1. A method of multi-spectral image processing, the method comprising: (a) calculating a pixel color feature (PCF) value for each pixel of a spectral component of a multi-spectral image to generate a set of PCF data; (b) applying one or more pairs of adaptive smoothing filters to each said spectral PCF value, wherein each pair of adaptive smoothing filters has fine and coarse smoothing parameters, thereby determining new pairs of PCF values for each said spectral PCF value at each pixel of a spectral component of the multi-spectral image; (c) differencing between said pairs of new PCF values thereby producing respective spatial adaptive bandpass filters (ABPFs) for each pixel of a spectral component, whereby said different smoothing parameters define degrees of local spectral contrast with respect to adjacent regions in said spectral component; (d) enhancing and detecting the multi-spectral image regions by combining a measured degree of the local spectral contrasts of said ABPFs from adjoining or overlapping regions from one or more spectral components of the multi-spectral image.