Change Detection in Sequences of Images by Multifractal Analysis

International audienceIn this work, we propose a multifractal approach to the problem of change detection in image sequences, such as registrated remotely sensed images of the same scene or sequences of medical images. We show that the multifractal analysis of images -- based on a modelisation of the two-dimensional signal as measure -- can be of great help if we want to detect changes without any a priori knowledge of the objects to be extracted. We first present a simple change detection method based on the classical multifractal analysis of images w.r.t. the Lebesgue measure. We then describe an improved method based on the analysis of images w.r.t. a reference measure, which in this case is the first image of the sequence. We finally show some results on real data

Fractal and multifractal analysis of PET-CT images of metastatic melanoma before and after treatment with ipilimumab

PET/CT with F-18-Fluorodeoxyglucose (FDG) images of patients suffering from metastatic melanoma have been analysed using fractal and multifractal analysis to assess the impact of monoclonal antibody ipilimumab treatment with respect to therapy outcome. Our analysis shows that the fractal dimensions which describe the tracer dispersion in the body decrease consistently with the deterioration of the patient therapeutic outcome condition. In 20 out-of 24 cases the fractal analysis results match those of the medical records, while 7 cases are considered as special cases because the patients have non-tumour related medical conditions or side effects which affect the results. The decrease in the fractal dimensions with the deterioration of the patient conditions (in terms of disease progression) are attributed to the hierarchical localisation of the tracer which accumulates in the affected lesions and does not spread homogeneously throughout the body. Fractality emerges as a result of the migration patterns which the malignant cells follow for propagating within the body (circulatory system, lymphatic system). Analysis of the multifractal spectrum complements and supports the results of the fractal analysis. In the kinetic Monte Carlo modelling of the metastatic process a small number of malignant cells diffuse throughout a fractal medium representing the blood circulatory network. Along their way the malignant cells engender random metastases (colonies) with a small probability and, as a result, fractal spatial distributions of the metastases are formed similar to the ones observed in the PET/CT images. In conclusion, we propose that fractal and multifractal analysis has potential application in the quantification of the evaluation of PET/CT images to monitor the disease evolution as well as the response to different medical treatments.Comment: 38 pages, 9 figure

Scaling analyses based on wavelet transforms for the Talbot effect

The fractal properties of the transverse Talbot images are analysed with two well-known scaling methods, the wavelet transform modulus maxima (WTMM) and the wavelet transform multifractal detrended fluctuation analysis (WT-MFDFA). We use the widths of the singularity spectra, Delta alpha=alpha_H-alpha_min, as a characteristic feature of these Talbot images. The tau scaling exponents of the q moments are linear in q within the two methods, which proves the monofractality of the transverse diffractive paraxial field in the case of these imagesComment: 9 pages, 6 figures, version accepted at Physica

Two-dimensional matrix algorithm using detrended fluctuation analysis to distinguish Burkitt and diffuse large B-cell lymphoma

Copyright © 2012 Rong-Guan Yeh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.A detrended fluctuation analysis (DFA) method is applied to image analysis. The 2-dimensional (2D) DFA algorithms is proposed for recharacterizing images of lymph sections. Due to Burkitt lymphoma (BL) and diffuse large B-cell lymphoma (DLBCL), there is a significant different 5-year survival rates after multiagent chemotherapy. Therefore, distinguishing the difference between BL and DLBCL is very important. In this study, eighteen BL images were classified as group A, which have one to five cytogenetic changes. Ten BL images were classified as group B, which have more than five cytogenetic changes. Both groups A and B BLs are aggressive lymphomas, which grow very fast and require more intensive chemotherapy. Finally, ten DLBCL images were classified as group C. The short-term correlation exponent α1 values of DFA of groups A, B, and C were 0.370 ± 0.033, 0.382 ± 0.022, and 0.435 ± 0.053, respectively. It was found that α1 value of BL image was significantly lower (P < 0.05) than DLBCL. However, there is no difference between the groups A and B BLs. Hence, it can be concluded that α1 value based on DFA statistics concept can clearly distinguish BL and DLBCL image.National Science Council (NSC) of Taiwan the Center for Dynamical Biomarkers and Translational Medicine, National Central University, Taiwan (also sponsored by National Science Council)

25 Years of Self-Organized Criticality: Numerical Detection Methods

The detection and characterization of self-organized criticality (SOC), in both real and simulated data, has undergone many significant revisions over the past 25 years. The explosive advances in the many numerical methods available for detecting, discriminating, and ultimately testing, SOC have played a critical role in developing our understanding of how systems experience and exhibit SOC. In this article, methods of detecting SOC are reviewed; from correlations to complexity to critical quantities. A description of the basic autocorrelation method leads into a detailed analysis of application-oriented methods developed in the last 25 years. In the second half of this manuscript space-based, time-based and spatial-temporal methods are reviewed and the prevalence of power laws in nature is described, with an emphasis on event detection and characterization. The search for numerical methods to clearly and unambiguously detect SOC in data often leads us outside the comfort zone of our own disciplines - the answers to these questions are often obtained by studying the advances made in other fields of study. In addition, numerical detection methods often provide the optimum link between simulations and experiments in scientific research. We seek to explore this boundary where the rubber meets the road, to review this expanding field of research of numerical detection of SOC systems over the past 25 years, and to iterate forwards so as to provide some foresight and guidance into developing breakthroughs in this subject over the next quarter of a century.Comment: Space Science Review series on SO

Scale-Free and Multifractal Time Dynamics of fMRI Signals during Rest and Task

Scaling temporal dynamics in functional MRI (fMRI) signals have been evidenced for a decade as intrinsic characteristics of ongoing brain activity (Zarahn et al., 1997). Recently, scaling properties were shown to fluctuate across brain networks and to be modulated between rest and task (He, 2011): notably, Hurst exponent, quantifying long memory, decreases under task in activating and deactivating brain regions. In most cases, such results were obtained: First, from univariate (voxelwise or regionwise) analysis, hence focusing on specific cognitive systems such as Resting-State Networks (RSNs) and raising the issue of the specificity of this scale-free dynamics modulation in RSNs. Second, using analysis tools designed to measure a single scaling exponent related to the second order statistics of the data, thus relying on models that either implicitly or explicitly assume Gaussianity and (asymptotic) self-similarity, while fMRI signals may significantly depart from those either of those two assumptions (Ciuciu et al., 2008; Wink et al., 2008). To address these issues, the present contribution elaborates on the analysis of the scaling properties of fMRI temporal dynamics by proposing two significant variations. First, scaling properties are technically investigated using the recently introduced Wavelet Leader-based Multifractal formalism (WLMF; Wendt et al., 2007). This measures a collection of scaling exponents, thus enables a richer and more versatile description of scale invariance (beyond correlation and Gaussianity), referred to as multifractality. Also, it benefits from improved estimation performance compared to tools previously used in the literature. Second, scaling properties are investigated in both RSN and non-RSN structures (e.g., artifacts), at a broader spatial scale than the voxel one, using a multivariate approach, namely the Multi-Subject Dictionary Learning (MSDL) algorithm (Varoquaux et al., 2011) that produces a set of spatial components that appear more sparse than their Independent Component Analysis (ICA) counterpart. These tools are combined and applied to a fMRI dataset comprising 12 subjects with resting-state and activation runs (Sadaghiani et al., 2009). Results stemming from those analysis confirm the already reported task-related decrease of long memory in functional networks, but also show that it occurs in artifacts, thus making this feature not specific to functional networks. Further, results indicate that most fMRI signals appear multifractal at rest except in non-cortical regions. Task-related modulation of multifractality appears only significant in functional networks and thus can be considered as the key property disentangling functional networks from artifacts. These finding are discussed in the light of the recent literature reporting scaling dynamics of EEG microstate sequences at rest and addressing non-stationarity issues in temporally independent fMRI modes