Magnetic resonance (MR) imaging systems produce spatial distribution estimates of proton density, relaxation time, and flow, in a two dimensional matrix form that is analogous to that of the image data obtained from multispectral imaging satellites. Advanced NASA satellite image processing offers sophisticated multispectral analysis of MR images. Spin echo and inversion recovery pulse sequence images were entered in a digital format compatible with satellite images and accurately registered pixel by pixel. Signatures of each tissue class were automatically determined using both supervised and unsupervised classification. Overall tissue classification was obtained in the form of a theme map. In MR images of the brain, for example, the classes included CSF, gray matter, white matter, subcutaneous fat, muscle, and bone. These methods provide an efficient means of identifying subtle relationships in a multi-image MR study.