Umeå University's logo

Manual segmentations are considered the gold standard for ground truth in machine learning applications. Such tasks are tedious and time-consuming, albeit necessary to train reliable models. In this work, we present a dataset with expert segmentations of the prostatic zones and urethra for 200 randomly selected patients from the PROSTATEx dataset. Notably, independent duplicate segmentations were performed for 40 patients, providing inter-reader variability data. This results in a total of 240 segmentations. This dataset can be used to train machine learning models or serve as an external test set for evaluating models trained on private data, thereby addressing a current gap in the field. The delineated structures and terminology adhere to the latest Prostate Imaging Reporting and Data Systems v2.1 guidelines, ensuring consistency.

The purpose of this work is to investigate spatial statistical modelling approaches to improve contrast agent quantification in dynamic contrast enhanced MRI, by utilising the spatial dependence among image voxels. Bayesian hierarchical models (BHMs), such as Besag model and Leroux model, were studied using simulated MRI data. The models were built on smaller images where spatial dependence can be incorporated, and then extended to larger images using the maximum a posteriori (MAP) method. Notable improvements on contrast agent concentration estimation were obtained for both smaller and larger images. For smaller images: the BHMs provided substantial improved estimates in terms of the root mean squared error (rMSE), compared to the estimates from the existing method for a noise level equivalent of a 12-channel head coil at 3T. Moreover, Leroux model outperformed Besag models with two different dependence structures. Specifically, the Besag models increased the estimation precision by 27% around the peak of the dynamic curve, while the Leroux model improved the estimation by 40% at the peak, compared with the existing estimation method. For larger images: the proposed MAP estimators showed clear improvements on rMSE for vessels, tumor rim and white matter.

Striatal dopamine is involved in facilitation of motor action as well as various cognitive and emotional functions. Positron emission tomography (PET) is the primary imaging method used to investigate dopamine function in humans. Previous PET studies have shown striatal dopamine release during simple finger tapping in both the putamen and the caudate. It is likely that dopamine release in the putamen is related to motor processes while dopamine release in the caudate could signal sustained cognitive component processes of the task, but the poor temporal resolution of PET has hindered firm conclusions. In this study we simultaneously collected [11C]Raclopride PET and functional Magnetic Resonance Imaging (fMRI) data while participants performed finger tapping, with fMRI being able to isolate activations related to individual tapping events. The results revealed fMRI-PET overlap in the bilateral putamen, which is consistent with a motor component process. Selective PET responses in the caudate, ventral striatum, and right posterior putamen, were also observed but did not overlap with fMRI responses to tapping events, suggesting that these reflect non-motor component processes of finger tapping. Our findings suggest an interplay between motor and non-motor-related dopamine release during simple finger tapping and illustrate the potential of hybrid PET-fMRI in revealing distinct component processes of cognitive functions.

Purpose: To enable the evaluation of locoregional disease in the on-going RECTOPET (REctal Cancer Trial on PET/MRI/CT) study; a methodology to match mesorectal imaging findings to histopathology is presented, along with initial observations.

Methods: FDG-PET/MRI examinations were performed in twenty-four consecutively included patients with rectal adenocarcinoma. In nine patients, of whom five received neoadjuvant treatment, a postoperative MRI of the surgical specimen was performed. The pathological cut-out was performed according to clinical routine with the addition of photo documentation of each slice of the surgical specimen, meticulously marking the location, size, and type of pathology of each mesorectal finding. This allowed matching individual nodal structures from preoperative MRI, via the specimen MRI, to histopathology.

Results: Preoperative MRI identified 197 mesorectal nodal structures, of which 92 (47%) could be anatomically matched to histopathology. Of the matched nodal structures identified in both MRI and histopathology, 25% were found to be malignant. These malignant structures consisted of lymph nodes (43%), tumour deposits (48%), and extramural venous invasion (9%). One hundred eleven nodal structures (55%) could not be matched anatomically. Of these, 97 (87%) were benign lymph nodes, and 14 (13%) were malignant nodal structures. Five were malignant lymph nodes, and nine were tumour deposits, all of which had a short axis diameter < 5 mm.

Conclusions: We designed a method able to anatomically match and study the characteristics of individual mesorectal nodal structures, enabling further research on the impact of each imaging modality. Initial observations suggest that small malignant nodal structures assessed as lymph nodes in MRI often comprise other forms of mesorectal tumour spread.

Trial registration: Clinical Trials Identifier: NCT03846882.

Haralick texture features are common texture descriptors in image analysis. To compute the Haralick features, the image gray-levels are reduced, a process called quantization. The resulting features depend heavily on the quantization step, so Haralick features are not reproducible unless the same quantization is performed. The aim of this work was to develop Haralick features that are invariant to the number of quantization gray-levels. By redefining the gray-level co-occurrence matrix (GLCM) as a discretized probability density function, it becomes asymptotically invariant to the quantization. The invariant and original features were compared using logistic regression classification to separate two classes based on the texture features. Classifiers trained on the invariant features showed higher accuracies, and had similar performance when training and test images had very different quantizations. In conclusion, using the invariant Haralick features, an image pattern will give the same texture feature values independent of image quantization.

PurposeThe role of hybrid imaging using F-18-fluoro-2-deoxy-D-glucose positron-emission tomography (FDG-PET), computed tomography (CT) and magnetic resonance imaging (MRI) to improve preoperative evaluation of rectal cancer is largely unknown. To investigate this, the RECTOPET (REctal Cancer Trial on PET/MRI/CT) study has been launched with the aim to assess staging and restaging of primary rectal cancer. This report presents the study workflow and the initial experiences of the impact of PET/CT on staging and management of the first patients included in the RECTOPET study.MethodsThis prospective cohort study, initiated in September 2016, is actively recruiting patients from Region Vasterbotten in Sweden. This pilot study includes patients recruited and followed up until December 2017. All patients had a biopsy-verified rectal adenocarcinoma and underwent a minimum of one preoperative FDG-PET/CT and FDG-PET/MRI examination. These patients were referred to the colorectal cancer multidisciplinary team meeting at Umea University Hospital. All available data were evaluated when making management recommendations. The clinical course was noted and changes consequent to PET imaging were described; surgical specimens underwent dedicated MRI for anatomical matching between imaging and histopathology.ResultsTwenty-four patients have so far been included in the study. Four patients were deemed unresectable, while 19 patients underwent or were scheduled for surgery; one patient was enrolled in a watch-and-wait programme after restaging. Consequent to taking part in the study, two patients were upstaged to M1 disease: one patient was diagnosed with a solitary hepatic metastasis detected using PET/CT and underwent metastasectomy prior to rectal cancer surgery, while one patient with a small, but metabolically active, lung nodulus experienced no change of management. PET/MRI did not contribute to any recorded change in patient management.ConclusionsThe RECTOPET study investigating the role of PET/CT and PET/MRI for preoperative staging of primary rectal cancer patients will provide novel data that clarify the value of adding hybrid to conventional imaging, and the role of PET/CT versus PET/MRI.Trial registrationNCT03846882.

Aim: To evaluate if quantifying proton density fat fraction (PDFF) would be useful in separating lipoma, atypical lipomatous tumour (ALT) and liposarcoma in the extremities and trunk. In addition, differentiating ALT versus non-classical lipomas using magnetic resonance imaging (MRI)-based fatty acidcomposition (FAC) and three-dimensional (3D) texture analysis was tested.

Material and methods: This prospective study (undertaken between 2014–2017; comprising 20 women, 21 men) was approved by the Regional Ethical Review Board and informed consent was obtained from all participants. For PDFF and FAC 3D spoiled gradient multi-echo images were acquired. PDFF was analysed in 16 lipomas (25–76 years), 14 ALTs (42–78 years) and 11 myxoid liposarcomas (31–68 years). The difference of mean PDFF was tested with one-way analysis of variance. A support vector machine algorithm was used to find the separating mean PDFF values.

Results: Mean PDFF for lipomas was 90% (range 76–98%), for ALT 83% (range 62–91%), and for liposarcoma 4% (range 0–21%). The difference of mean PDFF for liposarcomas versus ALT and lipoma was significant (p=0.0001, for both), and for ALT versus lipoma (p=0.021). The optimal threshold for separating liposarcoma from ALT and lipoma was 41.5%, and for ALT and lipoma 85%. Texture analysis could not separate ALT and non-classical lipomas, while the difference for FAC unsaturation degree was significant (p=0.013).

Conclusion: Measuring PDFF is a promising complement to standard MRI, to separate liposarcomas from ALT and lipomas. Lipomas that are not solely composed of fat cannot confidently be separated from ALT using PDFF, FAC, or texture analysis.

The Haralick texture features are common in the image analysis literature, partly because of their simplicity and because their values can be interpreted. It was recently observed that the Haralick texture features are very sensitive to the size of the GLCM that was used to compute them, which led to a new formulation that is invariant to the GLCM size. However, these new features still depend on the sample size used to compute the GLCM, i.e. the size of the input image region-of-interest (ROI).

The purpose of this work was to investigate the performance of density estimation methods for approximating the GLCM and subsequently the corresponding invariant features.

Three density estimation methods were evaluated, namely a piece-wise constant distribution, the Parzen-windows method, and the Gaussian mixture model. The methods were evaluated on 29 different image textures and 20 invariant Haralick texture features as well as a wide range of different ROI sizes.

The results indicate that there are two types of features: those that have a clear minimum error for a particular GLCM size for each ROI size, and those whose error decreases monotonically with increased GLCM size. For the first type of features, the Gaussian mixture model gave the smallest errors, and in particular for small ROI sizes (less than about 20×20).

In conclusion, the Gaussian mixture model is the preferred method for the first type of features (in particular for small ROIs). For the second type of features, simply using a large GLCM size is preferred.

Purpose: Simultaneous collection of PET and MR data for radiotherapy purposes are useful for, for example, target definition and dose escalations. However, a prerequisite for using PET/MR in the radiotherapy workflow is the ability to image the patient in treatment position. The aim of this work was to adapt a GE SIGNA PET/MR scanner to image patients for radiotherapy treatment planning and evaluate the impact on signal-to-noise (SNR) of the MR images, and the accuracy of the PET attenuation correction. Method: A flat tabletop and a coil holder were developed to image patients in the treatment position, avoid patient contour deformation, and facilitate attenuation correction of flex coils. Attenuation corrections for the developed hardware and an anterior array flex coil were also measured and implemented to the PET/MR system to minimize PET quantitation errors. The reduction of SNR in the MR images due to the added distance between the coils and the patient was evaluated using a large homogenous saline-doped water phantom, and the activity quantitation errors in PET imaging were evaluated with and without the developed attenuation corrections. Result: We showed that the activity quantitation errors in PET imaging were within ±5% when correcting for attenuation of the flat tabletop, coil holder, and flex coil. The SNR of the MRI images were reduced to 74% using the tabletop, and 66% using the tabletop and coil holders. Conclusion: We present a tabletop and coil holder for an anterior array coil to be used with a GE SIGNA PET/MR scanner, for scanning patients in the radiotherapy work flow. Implementing attenuation correction of the added hardware from the radiotherapy setup leads to acceptable PET image quantitation. The drop in SNR in MR images may require adjustment of the imaging protocols.