Diagnosis of human bladder cancer in untreated tissue sections is achieved by using imaging data from desorption electrospray ionization mass spectrometry (DESI-MS) combined with multivariate statistical analysis. We use the distinctive DESI-MS glycerophospholipid (GP) mass spectral profiles to visually characterize and formally classify twenty pairs (40 tissue samples) of human cancerous and adjacent normal bladder tissue samples. The individual ion images derived from the acquired profiles correlate with standard histological hematoxylin and eosin (H&E)-stained serial sections. The profiles allow us to classify the disease status of the tissue samples with high accuracy as judged by reference histological data. To achieve this, the data from the twenty pairs were divided into a training set and a validation set. Spectra from the tumor and normal regions of each of the tissue sections in the training set were used for orthogonal projection to latent structures (O-PLS) treated partial least-square discriminate analysis (PLS-DA). This predictive model was then validated by using the validation set and showed a 5% error rate for classification and a misclassification rate of 12%. It was also used to create synthetic images of the tissue sections showing pixel-by-pixel disease classification of the tissue and these data agreed well with the independent classification that uses histological data by a certified pathologist. This represents the first application of multivariate statistical methods for classification by ambient ionization although these methods have been applied previously to other MS imaging methods. The results are encouraging in terms of the development of a method that could be utilized in a clinical setting through visualization and diagnosis of intact tissue.