Methodology, with a programmable-processor imaging system, for control and determination of breast lesion viscoelastic properties with the use of a creep-like test. Two dimensional reconstruction maps are used for different parameters of a linear viscoelastic model. Description of different aspects of the test used on live subjects and suitability of a 1-D inversion model in capturing different viscoelasticity parameters. An automated methodology for the selection of a region of interest derived only from the appearance of the breast lesion on pre-compressed B-mode images. Based on the ROI and estimated viscoelasticity parameters, contrast values are determined that facilitate the enhanced differentiation of breast mass. Employing the methodology in a large group of patients provides better understanding of variations of different viscoelasticity parameters in different types of breast lesion and helps to identify new biomarkers for enhanced differentiation of benign from malignant cases.