A method for content-based image retrieval for the classification of breast density from mammographic imagery is described. The breast density is characterized through the Fisher linear discriminants (FLD) extracted from the Principal Component Analysis (PCA). Unlike PCA, the FLD provides a very discriminative representation of the mammographic images in terms of the breast density. Various exemplary methods, systems and computer program products are also disclosed.