The present disclosure involves a computer-implemented unmixing algorithm that employs a least-square method involving image patches, and a computer-implemented unmixing or color deconvolution method that incorporates spatial smoothness and structure continuity constraints, for example, into a neighborhood graph regularizer, such as using a graph to enforce the pixel value similarities for those pixels in the same neighborhood.