1. A method of multi-spectral image processing, the method comprising: (a) receiving spectral values for each pixel of a spectral component of a multi-spectral image as digital data; (b) applying one or more pairs of adaptive smoothing filters to each spectral value, wherein each pair of adaptive smoothing filters has fine and coarse smoothing parameters, 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; (c) calculating a pixel color feature output (PCFO) value for each filtered pixel of the multi-spectral image to generate a set of PCFO data; (d) enhancing and detecting the multi-spectral image regions by combining a measured degree of the local spectral contrasts of said PCFOs from adjoining or overlapping regions from one or more spectral components of the multi-spectral image.