The number of biomedical image analysis challenges organized per year is steadily increasing. These international competitions have the purpose of benchmarking algorithms on common data sets, typically to identify the best method for a given problem. Recent research, however, revealed that common practice related to challenge reporting does not allow for adequate interpretation and reproducibility of results. To address the discrepancy between the impact of challenges and the quality (control), the Biomedical Image Analysis ChallengeS (BIAS) initiative developed a set of recommendations for the reporting of challenges. The BIAS statement aims to improve the transparency of the reporting of a biomedical image analysis challenge regardless of field of application, image modality or task category assessed. This article describes how the BIAS statement was developed and presents a checklist which authors of biomedical image analysis challenges are encouraged to include in their submission when giving a paper on a challenge into review. The purpose of the checklist is to standardize and facilitate the review process and raise interpretability and reproducibility of challenge results by making relevant information explicit.

PURPOSE: Pathologists often look at whole slide images (WSIs) at low magnification to find potentially important regions and then zoom in to higher magnification to perform more sophisticated analysis of the tissue structures. Many automated methods of WSI analysis attempt to preprocess the down-sampled image in order to select salient regions which are then further analyzed by a more computationally intensive step at full magnification. Although it can greatly reduce processing times, this process may lead to small potentially important regions being overlooked at low magnification. We propose a texture analysis technique to ease the processing of H&E stained WSIs by triaging clinically important regions. METHOD: Image patches randomly selected from the whole tissue area were divided into smaller tiles and Gaussian-like texture filters were applied to them. Texture filter responses from each tile were combined together and statistical measures were derived from their histograms of responses. Bag of visual words pipeline was then employed to combine extracted features from tiles to form one histogram of words per every image patch. A support vector machine classifier was trained using the calculated histograms of words to be able to distinguish between clinically relevant and irrelevant patches. RESULT: Experimental analysis on 5151 image patches from 10 patient cases (65 tissue slides) indicated that our proposed texture technique out-performed two previously proposed colour and intensity based methods with an area under the ROC curve of 0.87. CONCLUSION: Texture features can be employed to triage clinically important areas within large WSIs.

This paper presents a robust atlas-based segmentation (ABS) algorithm for segmentation of the breast boundary in 3-D MR images. The proposed algorithm combines the well-known methodologies of ABS namely probabilistic atlas and atlas selection approaches into a single framework where two configurations are realized. The algorithm uses phase congruency maps to create an atlas which is robust to intensity variations. This allows an atlas derived from images acquired with one MR imaging sequence to be used to segment images acquired with a different MR imaging sequence and eliminates the need for intensity-based registration. Images acquired using a Dixon sequence were used to create an atlas which was used to segment both Dixon images (intra-sequence) and T1-weighted images (inter-sequence). In both cases, highly accurate results were achieved with the median Dice similarity coefficient values of 94% ±4% and 87 ±6.5%, respectively.

Assessing tumor response to therapy is a crucial step in personalized treatments. Pharmacokinetic (PK) modeling provides quantitative information about tumor perfusion and vascular permeability that are associated with prognostic factors. A fundamental step in most PK analyses is calculating the signal that is generated in the tumor vasculature. This signal is usually inseparable from the extravascular extracellular signal. It was shown previously using in vivo and phantom experiments that independent component analysis (ICA) is capable of calculating the intravascular time-intensity curve in dynamic contrast enhanced (DCE)-MRI. A novel adaptive complex independent component analysis (AC-ICA) technique is developed in this study to calculate the intravascular time-intensity curve and separate this signal from the DCE-MR images of tumors. The use of the complex-valued DCE-MRI images rather than the commonly used magnitude images satisfied the fundamental assumption of ICA, i.e., linear mixing of the sources. Using an adaptive cost function in ICA through estimating the probability distribution of the tumor vasculature at each iteration resulted in a more robust and accurate separation algorithm. The AC-ICA algorithm provided a better estimate for the intravascular time-intensity curve than the previous ICA-based method. A simulation study was also developed in this study to realistically simulate DCE-MRI data of a leaky tissue mimicking phantom. The passage of the MR contrast agent through the leaky phantom was modeled with finite element analysis using a diffusion model. Once the distribution of the contrast agent in the imaging field of view was calculated, DCE-MRI data was generated by solving the Bloch equation for each voxel at each time point. The intravascular time-intensity curve calculation results were compared to the previously proposed ICA-based intravascular time-intensity curve calculation method that applied ICA to the magnitude of the DCE-MRI data (Mag-ICA) using both simulated and experimental tissue mimicking phantoms. The AC-ICA demonstrated superior performance compared to the Mag-ICA method. AC-ICA provided more accurate estimate of intravascular time-intensity curve, having smaller error between the calculated and actual intravascular time-intensity curves compared to the Mag-ICA. Furthermore, it showed higher robustness in dealing with datasets with different resolution by providing smaller variation between the results of each datasets and having smaller difference between the intravascular time-intensity curves of various resolutions. Thus, AC-ICA has the potential to be used as the intravascular time-intensity curve calculation method in PK analysis and could lead to more accurate PK analysis for tumors.

Early detection of breast cancer is one of the most important factors in determining prognosis for women with malignant tumors. Dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) has been shown to be the most sensitive modality for screening high-risk women. Computer-aided diagnosis (CAD) systems have the potential to assist radiologists in the early detection of cancer. A key component of the development of such a CAD system will be the selection of an appropriate classification function responsible for separating malignant and benign lesions. The purpose of this study is to evaluate the effects of variations in temporal feature vectors and kernel functions on the separation of malignant and benign DCE-MRI breast lesions by support vector machines (SVMs). We also propose and demonstrate a classifier visualization and evaluation technique. We show that SVMs provide an effective and flexible framework from which to base CAD techniques for breast MRI, and that the proposed classifier visualization technique has potential as a mechanism for the evaluation of classification solutions.

Contrast enhanced ultrasound imaging provides a real-time tool for evaluating vasculature in the liver. Primary liver cancer is known to be perfused exclusively by blood from the hepatic artery, whereas normal liver is also supplied by the portal vein. Visual separation of two different phases of enhancement from the independent feeding vessels is important for diagnosis but remains a challenge. This paper presents a method of using factor analysis for extracting distinct time-intensity curves. A key component to this extraction is the clustering of measured bolus curves and their projection onto a positivity domain to obtain nonnegative curves. This technique provides complementary images representing spatial loadings on each curve. As little as 1% of the data is required to contain unmixed signals to extract time-intensity curves that correlate well with true curves. A method of combining this information to display a regional hepatic perfusion image is proposed, and results are tested on a set of 10 patients. Region of interest analysis suggests it is possible to detect changes in the hepatic perfusion index of liver lesions relative to normal liver parenchyma using contrast ultrasound.