Mads Nielsen

Mammographic density is a well-established risk factor for breast cancer. We investigated the association between three different methods of measuring density or parenchymal pattern/texture on digitized film-based mammograms, and examined to what extent textural features independently and jointly with density can improve the ability to identify screening women at increased risk of breast cancer. The study included 121 cases and 259 age- and time matched controls based on a cohort of 14,736 women with negative screening mammograms from a population-based screening programme in Denmark in 2007 (followed until 31 December 2010). Mammograms were assessed using the Breast Imaging-Reporting and Data System (BI-RADS) density classification, Tabár’s classification on parenchymal patterns and a fully automated texture quantification technique. The individual and combined association with breast cancer was estimated using binary logistic regression to calculate Odds Ratios (ORs) and the area under the receiver operating characteristic (ROC) curves (AUCs). Cases showed significantly higher BI-RADS and texture scores on average than controls (p < 0.001). All three methods were individually able to segregate women into different risk groups showing significant ORs for BI-RADS D3 and D4 (OR: 2.37; 1.32–4.25 and 3.93; 1.88–8.20), Tabár’s PIII and PIV (OR: 3.23; 1.20–8.75 and 4.40; 2.31–8.38), and the highest quartile of the texture score (3.04; 1.63–5.67). AUCs for BI-RADS, Tabár and the texture scores (continuous) were 0.63 (0.57–0–69), 0.65 (0.59–0–71) and 0.63 (0.57–0–69), respectively. Combining two or more methods increased model fit in all combinations, demonstrating the highest AUC of 0.69 (0.63-0.74) when all three methods were combined (a significant increase from standard BI-RADS alone). Our findings suggest that the (relative) amount of fibroglandular tissue (density) and mammographic structural features (texture/parenchymal pattern) jointly can improve risk segregation of screening women, using information already available from normal screening routine, in respect to future personalized screening strategies.

Brain atrophy from structural magnetic resonance images (MRIs) is widely used as an imaging surrogate marker for Alzheimers disease. Their utility has been limited due to the large degree of variance and subsequently high sample size estimates. The only consistent and reasonably powerful atrophy estimation methods has been the boundary shift integral (BSI). In this paper, we first propose a tensor-based morphometry (TBM) method to measure voxel-wise atrophy that we combine with BSI. The combined model decreases the sample size estimates significantly when compared to BSI and TBM alone. © (2016) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.

Mammographic breast density and parenchymal patterns are well-established risk factors for breast cancer. We aimed to report inter-observer agreement on three different subjective ways of assessing mammographic density and parenchymal pattern, and secondarily to examine what potential impact reproducibility has on relative risk estimates of breast cancer. This retrospective case–control study included 122 cases and 262 age- and time matched controls (765 breasts) based on a 2007 screening cohort of 14,736 women with negative screening mammograms from Bispebjerg Hospital, Copenhagen. Digitised randomized film-based mammograms were classified independently by two readers according to two radiological visual classifications (BI-RADS and Tabár) and a computerized interactive threshold technique measuring area-based percent mammographic density (denoted PMD). Kappa statistics, Intraclass Correlation Coefficient (ICC) (equivalent to weighted kappa), Pearson’s linear correlation coefficient and limits-of-agreement analysis were used to evaluate inter-observer agreement. High/low-risk agreement was also determined by defining the following categories as high-risk: BI-RADS’s D3 and D4, Tabár’s PIV and PV and the upper two quartiles (within density range) of PMD. The relative risk of breast cancer was estimated using logistic regression to calculate odds ratios (ORs) adjusted for age, which were compared between the two readers. Substantial inter-observer agreement was seen for BI-RADS and Tabár (κ=0.68 and 0.64) and agreement was almost perfect when ICC was calculated for the ordinal BI-RADS scale (ICC=0.88) and the continuous PMD measure (ICC=0.93). The two readers judged 5% (PMD), 10% (Tabár) and 13% (BI-RADS) of the women to different high/low-risk categories, respectively. Inter-reader variability showed different impact on the relative risk of breast cancer estimated by the two readers on a multiple-category scale, however, not on a high/low-risk scale. Tabár’s pattern IV demonstrated the highest ORs of all density patterns investigated. Our study shows the Tabár classification has comparable inter-observer reproducibility with well tested density methods, and confirms the association between Tabár’s PIV and breast cancer. In spite of comparable high inter-observer agreement for all three methods, impact on ORs for breast cancer seems to differ according to the density scale used. Automated computerized techniques are needed to fully overcome the impact of subjectivity.

Introduction: Breast density has been established as a major risk factor for breast cancer. We have previously demonstrated that mammographic texture resemblance (MTR), recognizing the local texture patterns of the mammogram, is also a risk factor for breast cancer, independent of percent breast density. We examine if these findings generalize to another population. Methods: Texture patterns were recorded in digitalized pre-diagnosis (3.7 years) film mammograms of a nested case–control study within the Dutch screening program (S1) comprising of 245 breast cancers and 250 matched controls. The patterns were recognized in the same study using cross-validation to form resemblance scores associated with breast cancer. Texture patterns from S1 were examined in an independent nested case–control study within the Mayo Mammography Health Study cohort (S2) of 226 cases and 442 matched controls: mammograms on average 8.5 years prior to diagnosis, risk factor information and percent mammographic density (PD) estimated using Cumulus were available. MTR scores estimated from S1, S2 and S1 + S2 (the latter two as cross-validations) were evaluated in S2. MTR scores were analyzed as both quartiles and continuously for association with breast cancer using odds ratios (OR) and adjusting for known risk factors including age, body mass index (BMI), and hormone usage. Results: The mean ages of S1 and S2 were 58.0 ± 5.7 years and 55.2 ± 10.5 years, respectively. The MTR scores on S1 showed significant capability to discriminate cancers from controls (area under the operator characteristics curve (AUC) = 0.63 ± 0.02, P <0.001), which persisted after adjustment for PD. S2 showed an AUC of 0.63, 0.61, and 0.60 based on PD, MTR scores trained on S2, and MTR scores trained on S1, respectively. When adjusted for PD, MTR scores of S2 trained on S1 showed an association with breast cancer for the highest quartile alone: OR in quartiles of controls as reference; 1.04 (0.59 to 1.81); 0.95 (0.52 to 1.74); 1.84 (1.10 to 3.07) respectively. The combined continuous model with both PD and MTR scores based on S1 had an AUC of 0.66 ± 0.03. Conclusions: The local texture patterns associated with breast cancer risk in S1 were also an independent risk factor in S2. Additional textures identified in S2 did not significantly improve risk segregation. Hence, the textural patterns that indicated elevated risk persisted under differences in X-ray technology, population demographics, follow-up time and geography.

In fields ranging from computer vision to signal processing and statistics, increasing computational power allows a move from classical linear models to models that incorporate non-linear phenomena. This shift has created interest in computational aspects of differential geometry, and solving optimization problems that incorporate non-linear geometry constitutes an important computational task. In this paper, we develop methods for numerically solving optimization problems over spaces of geodesics using numerical integration of Jacobi fields and second order derivatives of geodesic families. As an important application of this optimization strategy, we compute exact Principal Geodesic Analysis (PGA), a non-linear version of the PCA dimensionality reduction procedure. By applying the exact PGA algorithm to synthetic data, we exemplify the differences between the linearized and exact algorithms caused by the non-linear geometry. In addition, we use the numerically integrated Jacobi fields to determine sectional curvatures and provide upper bounds for injectivity radii.

Many classification/segmentation tasks in medical imaging are particularly challenging for machine learning algorithms because of the huge amount of training data required to cover biological variability. Learning methods scaling badly in the number of training data points may not be applicable. This may exclude powerful classifiers with good generalization performance such as standard non-linear support vector machines (SVMs). Further, many medical imaging problems have highly imbalanced class populations, because the object to be segmented has only few pixels/voxels compared to the background. This article presents a two-stage classifier for large-scale medical imaging problems. In the first stage, a classifier that is easily trainable on large data sets is employed. The class imbalance is exploited and the classifier is adjusted to correctly detect background with a very high accuracy. Only the comparatively few data points not identified as background are passed to the second stage. Here a powerful classifier with high training time complexity can be employed for making the final decision whether a data point belongs to the object or not. We applied our method to the problem of automatically segmenting tibial articular cartilage from knee MRI scans. We show that by using nearest neighbor (kNN) in the first stage we can reduce the amount of data for training a non-linear SVM in the second stage. The cascaded system achieves better results than the state-of-the-art method relying on a single kNN classifier.

Aortic calcification is a major risk factor for death from cardiovascular disease. We investigated the relationship between mortality and the composite markers of number, size, morphology and distribution of calcified plaques in the lumbar aorta.