Development of an image processing software for quantification of histological calcification staining images

Quantification of the histological staining images gives important insights in biomedical research. In wet lab, it is common to have some stains off the target to become unwanted noisy stains during the generation of histological staining images. The current tools designed for quantification of histological staining images do not consider such situations; instead, the stained region is identified based on assumptions that the background is pure and clean. The goal of this study is to develop a light software named Staining Quantification (SQ) tool which could handle the image quantification job with features for removing a large amount of unwanted stains blended or overlaid with Region of Interest (ROI) in complex scenarios. The core algorithm was based on the method of higher order statistics transformation, and local density filtering. Compared with two state-of-art thresholding methods (i.e. Otsu's method and Triclass thresholding method), the SQ tool outperformed in situations such as (1) images with weak positive signals and experimental caused dirty stains; (2) images with experimental counterstaining by multiple colors; (3) complicated histological structure of target tissues. The algorithm was developed in R4.0.2 with over a thousand in-house histological images containing Alizarin Red (AR) and Von Kossa (VK) staining, and was validated using external images. For the measurements of area and intensity in total and stained region, the average mean of difference in percentage between SQ and ImageJ were all less than 0.05. Using this as a criterion of successful image recognition, the success rate for all measurements in AR, VK and external validation batch were above 0.8. The test of Pearson's coefficient, difference between SQ and ImageJ, and difference of proportions between SQ and ImageJ were all significant at level of 0.05. Our results indicated that the SQ tool is well established for automatic histological staining image quantification.