A microwave imaging (MWI) methodology for early diagnosis of breast cancer is presented. Instead of generating a tomographic image of the breast, the proposed technique aims to reconstruct a map of malignant tumours inside the breast by adopting an extended form of factorization method. The implementation of factorization method requires (i) two multi-static scattered field measurements around the breast, which correspond to two different states of the breast, and (ii) the inhomogeneous Green's function associated to the breast. For this purpose, the paper proposes the use of contrast agents, which selectively increase the dielectric properties of the malign tissues. Two multi-static field measurements are collected before and after the administration of contrast agents. Later, the inhomogeneous Green's function of the breast is estimated by back-propagating the scattered field measurements, which are taken before the contrast agent usage. The feasibility and efficiency of the proposed technique are demonstrated with numerical examples that are performed on a slice of the realistic breast phantoms, which are derived from real three dimensional magnetic resonance imaging (3D-MRI) measurements.