This article is concerned with the development of discrimination procedure for seismic data based on feature extraction. These features are the maximum and minimum (maximin) values of Nondecimated Wavelets Transform (NDWT), which were extracted from the P component of the original earthquake and explosion traces. The method provides small error rates in both discrimination and clustering seismic wave datasets. It has been shown that our procedure performs as well as other procedures and in some cases better than other classifications methods. Also, the performance of our procedure is evaluated via a simulation study.