A method of image processing, including: (a) calculating at least one pixel color feature (PCF) value for each pixel in a color medical image to generate a set of PCF data; and (b) filtering the PCF data with at least one spatial adaptive bandpass filter (ABPF) to sort the pixels into physiologically significant regions; wherein the at least one PCF value for at least one pixel depends on at least 2 color components of the medical image.