Search

Many complex physical systems exhibit a rich variety of discrete behavioural modes. Often, the system complexity limits the applicability of standard modelling tools. Hence, understanding the underlying physics of different behaviours and distinguishing between them is challenging. Although traditional machine learning techniques could predict and classify behaviour well, typically they do not provide any meaningful insight into the underlying physics of the system. In this paper we present a novel method for extracting physically meaningful clusters of discrete behaviour from limited experimental observations. This method obtains a set of physically plausible functions that both facilitate behavioural clustering and aid in system understanding. We demonstrate the approach on the V-shaped falling paper system, a new falling paper type system that exhibits four distinct behavioural modes depending on a few morphological parameters. Using just 49 experimental observations, the method discovered a set of candidate functions that distinguish behaviours with an error of 2.04%, while also aiding insight into the physical phenomena driving each behaviour. © 2019 Howison et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Author summary Conferences are great venues for disseminating algorithmic bioinformatics results, but they unfortunately do not offer an opportunity to make major revisions in the way that journals do. As a result, it is not possible for authors to fix mistakes that might be easily correctable but nevertheless can cause the paper to be rejected. As a reviewer, I wish that I had the opportunity to tell the authors, “Hey, you forgot to do this really important thing, without which it is hard to accept the paper, but if you could go back and fix it, you might have a great paper for the conference.” This lack of a back and forth can be especially problematic for first-time submitters or those from outside the field, e.g., biologists. In this article, I outline Ten Simple Rules to follow when writing an algorithmic bioinformatics conference paper to avoid having it rejected.

I review my scientific career in terms of eight areas and 16 papers. The first two areas are associated with childhood. The other six are associated with residual maximum likelihood (REML), canonical transformation, inbreeding in selected populations, average information residual maximum likelihood (AIREML), the computer program ASReml and sampling-based estimation.

The Semantic Link Network is a general semantic model for modeling the structure and the evolution of complex systems. Various semantic links play different roles in rendering the semantics of complex system. One of the basic semantic links represents cause-effect relation, which plays an important role in representation and understanding. This paper verifies the role of the Semantic Link Network in representing the core of text by investigating the contribution of cause-effect link to representing the core of scientific papers. Research carries out with the following steps: (1) Two propositions on the contribution of cause-effect link in rendering the core of paper are proposed and verified through a statistical survey, which shows that the sentences on cause-effect links cover about 65% of key words within each paper on average. (2) An algorithm based on syntactic patterns is designed for automatically extracting cause-effect link from scientific papers, which recalls about 70% of manually annotated cause-effect links on average, indicating that the result adapts to the scale of data sets. (3) The effects of cause-effect link on four schemes of incorporating cause-effect link into the existing instances of the Semantic Link Network for enhancing the summarization of scientific papers are investigated. The experiments show that the quality of the summaries is significantly improved, which verifies the role of semantic links. The significance of this research lies in two aspects: (1) it verifies that the Semantic Link Network connects the important concepts to render the core of text; and, (2) it provides an evidence for realizing content services such as summarization, recommendation and question answering based on the Semantic Link Network, and it can inspire relevant research on content computing.

Quantifying the impact of a scholarly paper is of great significance, yet the effect of geographical distance of cited papers has not been explored. In this paper, we examine 30,596 papers published in Physical Review C, and identify the relationship between citations and geographical distances between author affiliations. Subsequently, a relative citation weight is applied to assess the impact of a scholarly paper. A higher-order weighted quantum PageRank algorithm is also developed to address the behavior of multiple step citation flow. Capturing the citation dynamics with higher-order dependencies reveals the actual impact of papers, including necessary self-citations that are sometimes excluded in prior studies. Quantum PageRank is utilized in this paper to help differentiating nodes whose PageRank values are identical. Refereed/Peer-reviewed

With the rapid development of information technology, scientific social networks (SSNs) have become the fastest and most convenient way for researchers to communicate with each other. Many published papers are shared via SSNs every day, resulting in the problem of information overload. How to appropriately recommend personalized and highly valuable papers for researchers is becoming more urgent. However, when recommending papers in SSNs, only a small amount of positive instances are available, leaving a vast amount of unlabelled data, in which negative instances and potential unseen positive instances are mixed together, which naturally belongs to One-Class Collaborative Filtering (OCCF) problem. Therefore, considering the extreme data imbalance and data sparsity of this OCCF problem, a hybrid approach of Social and Content aware One-class Recommendation of Papers in SSNs, termed SCORP, is proposed in this study. Unlike previous approaches recommended to address the OCCF problem, social information, which has been proved playing a significant role in performing recommendations in many domains, is applied in both the profiling of content-based filtering and the collaborative filtering to achieve superior recommendations. To verify the effectiveness of the proposed SCORP approach, a real-life dataset from CiteULike was employed. The experimental results demonstrate that the proposed approach is superior to all of the compared approaches, thus providing a more effective method for recommending papers in SSNs.

This article proposes an Enhanced Memetic Differential Evolution (EMDE) for designing digital filters which aim at detecting defects of the paper produced during an industrial process. Defect detection is handled by means of two Gabor filters and their design is performed by the EMDE. The EMDE is a novel adaptive evolutionary algorithm which combines the powerful explorative features of Differential Evolution with the exploitative features of three local search algorithms employing different pivot rules and neighborhood generating functions. These local search algorithms are the Hooke Jeeves Algorithm, a Stochastic Local Search, and Simulated Annealing. The local search algorithms are adaptively coordinated by means of a control parameter that measures fitness distribution among individuals of the population and a novel probabilistic scheme. Numerical results confirm that Differential Evolution is an efficient evolutionary framework for the image processing problem under investigation and show that the EMDE performs well. As a matter of fact, the application of the EMDE leads to a design of an efficiently tailored filter. A comparison with various popular metaheuristics proves the effectiveness of the EMDE in terms of convergence speed, stagnation prevention, and capability in detecting solutions having high performance.

Cyclic dominance of three species is a commonly occurring interaction dynamics, often denoted the rock-paper-scissors (RPS) game. Such a type of interactions is known to promote species coexistence. Here, we generalize recent results of Reichenbach et al. [Nature (London) 448, 1046 (2007)] of a four-state variant of the RPS game. We show that spiral formation takes place only without a conservation law for the total density. Nevertheless, in general, fast diffusion can destroy species coexistence. We also generalize the four-state model to slightly varying reaction rates. This is shown both analytically and numerically not to change pattern formation, or the effective wavelength of the spirals, and therefore not to alter the qualitative properties of the crossover to extinction.

Procedures for grouping students into homogeneous subsets have long interested educational researchers. The research reported in this paper is an investigation of a set of objective grouping procedures based on multivariate analysis considerations. Four multivariate functions that might serve as criteria for adequate grouping are given and discussed; a method for optimizing these functions is also described. The set of procedures is illustrated through application to data from two samples of students, each student with scores on either ten or eleven subtests of a criterion referenced mathematics inventory. The results indicate that the procedures discussed provide a promising means for grouping students to minimize classroom heterogeneity. (Author)

Research paper recommenders emerged over the last decade to ease finding publications relating to researchers' area of interest. The challenge was not just to provide researchers with very rich publications at any time, any place and in any form but to also offer the right publication to the right researcher in the right way. Several approaches exist in handling paper recommender systems. However, these approaches assumed the availability of the whole contents of the recommending papers to be freely accessible, which is not always true due to factors such as copyright restrictions. This paper presents a collaborative approach for research paper recommender system. By leveraging the advantages of collaborative filtering approach, we utilize the publicly available contextual metadata to infer the hidden associations that exist between research papers in order to personalize recommendations. The novelty of our proposed approach is that it provides personalized recommendations regardless of the research field and regardless of the user's expertise. Using a publicly available dataset, our proposed approach has recorded a significant improvement over other baseline methods in measuring both the overall performance and the ability to return relevant and useful publications at the top of the recommendation list.

Question: 1—Echocardiography images Images were acquired from end of diastole to end of systole at different phases (7 phases) within a cardiac cycle by a Vivid 3 GE Health echocardiography machine. To obtain worthy images at denoted phases, it would be dependent to operator (for the best clinical examination). Of course, having assumed images not only were acquired from A4C and A2C but also from short-axis views. Short axis views have not been used at your papers, so some radial, circumferential and sample rotatory data would be missed. Response: It is simply possible for us to record and analyze images of short-axis views. However, according to American Society of Echocardiography (ASE) [1, 2], in the modified Simpson’s rule which is a method for LV volume computation from 2-D echocardiography images, it suffices to acquire apical two and four chamber views (A4C, A2C) for left ventricular (LV) volume computation. Question: 2—Observed data yi′ s Echocardiography images paly role of observed data and they are symbolized by yi′ s. there is a sequence of observed data y1, y2,…,yN. N is the number of obtained images. In fact we have made a chain of displacement, rotation and pure strain or nonrigid transformation (deformation) by these observed data that is started from y1 and is ended up to yN. This chain has some conceptual interpretations of the elasticity of the global left ventricular motion/function and the regional LV fiber arrangement/movement. Response: Yes, it is correct. Question: 3—Embed observed data to a high dimensional space Response: Yes, it is correct. Question: 4—Image processing on yi′ s and geodesic distance between yi′ s Observed data y1, y2,…,yN are embed to a meshed surface sized number of pixels of an image. Translation, rotation and pure deformation of y1, y2,…,yN are occurred/and studied at this surface utilizing the mathematical elasticity theory. A distance metric is defined that computes distances between observed points. Medical interpretations of observed points over the time should be checked/and stated clearly at your manuscripts. Images are not usual images but also are images from the left ventricle. Motion (displacement and velocity), deformation (strain and strain rate) and torsion of each myocardial sample have to be extracted during a cardiac cycle in the mentioned surface. These are used on the creation of graph G referred to their papers. Response: In this stage, we have just computed the distance between the two images in a cardiac cycle using non-rigid image registration regardless of relationship between images over the time. Therefore by this method, the motion and longitudinal deformation (short-axis views are needed to get for radial and circumferential deformations) of LV myocardium between all of two images in a cardiac cycle are extracted (for example, between end systole image and end diastole image, etc.). Then, in other stage of our method (calculation of the matrix of pairwise geodesic distance and apply multidimensional scaling on the resulting geodesic distance matrix to construct the low dimensional data points xi), motion and longitudinal deformation of LV myocardium during a cardiac cycle are extracted. Question: 5—An Isomap f The main tool at this method is the function “f” which is an isometric map [2, 3]. I will be interested to know the structure/formula or relationship of this isometric function which has most probably a lot of practical information of the left ventricular function and structure. I mean, how we can get a good representation of such a this function “f” for heart science study? Is it “f” known as an algebraic function? What information has been coded to fibers (f−1) of this isometric function “f”? Response: Isomap algorithm is one of the most popular sets of nonlinear dimensionality reduction (NLDR) algorithms. This algorithm attempts to extract low-dimensional data points from input points in a high-dimensional space in such a way that pairwise geodesic distances (the distance between two points measured over the manifold) are preserved. So that nearby and far points in high-dimensional space map to nearby and far points in low-dimensional space. This method is frequently used for visualization of medial image set. Detail of this method is described in Tenenbaum et al. and Borg et al. [3, 4]. It should be noted that we do not change the structure or formula of isomap algorithm. We only define a new image distance function based on image registration in calculation of geodesic distance between yi′ s in isomap algorithm. With this modification, Isomap algorithm is specified for assessment of the left ventricular function and structure. However, precise mathematical descriptions of f and (f−1) applied to echocardiography images and using the other NLDR algorithms need further researches in the following of this article. Question: 6—A reconstructed curve on a 2D manifold space & hidden data xi′ s; f(yi) = xi Isometric function results in a reconstructive curve crossing from hidden data (xi′ s) in a surface. It’s natural that some information (maybe clinical information) would betransferred to these hidden points like the curve of LV volume changes and so on and so on. A main question is: what are these new points xi′ s exactly (medically point of view)? What the other data has been come out from these hidden points? I think all of these questions back to gain a good understanding of isometric function “f”. Response: These points xi′ s will probably enable physicians to diagnose and follow up many cardiac structures and functions that also open doors to much more research at imaging cardiology in the future. References 1. AlizadehSani Z, Shalbaf A, Behnam H, Shalbaf R: Automatic computation of left ventricularvolume changes over a cardiac cycle from echocardiography images by nonlinear dimensionality reduction. J Digit Imaging : July, 2014 2. Tenenbaum JB, de Silva V, Langford JC: global geometric framework for nonlinear dimensionalityreduction. Science 290:2319–2323, 2000. Reprint available online: http://web.mit.edu/cocosci/Papers/sci_reprint.pdf 3. Ledesma-Carbayo MJ, Kybic J, Desco M, et al.: Spatio-temporal nonrigid registration for ultrasoundcardiac motion estimation. IEEE Trans Med Imaging 24:1113–1126, 2005.

Many results have been obtained when studying scientific papers citations databases in a network perspective. Articles can be ranked according to their current in-degree and their future popularity or citation counts can even be predicted. The dynamical properties of such networks and the observation of the time evolution of their nodes started more recently. This work adopts an evolutionary perspective and proposes an original algorithm for the construction of genealogical trees of scientific papers on the basis of their citation count evolution in time. The fitness of a paper now amounts to its in-degree growing trend and a "dying" paper will suddenly see this trend declining in time. It will give birth and be taken over by some of its most prevalent citing "offspring". Practically, this might be used to trace the successive published milestones of a research field., SCOPUS: ar.j, info:eu-repo/semantics/published

We study the oscillatory dynamics in the generic three-species rock-paper-scissors games with mutations. In the mean-field limit, different behaviors are found: (a) for high mutation rate, there is a stable interior fixed point with coexistence of all species; (b) for low mutation rates, there is a region of the parameter space characterized by a limit cycle resulting from a Hopf bifurcation; (c) in the absence of mutations, there is a region where heteroclinic cycles yield oscillations of large amplitude (not robust against noise). After a discussion on the main properties of the mean-field dynamics, we investigate the stochastic version of the model within an individual-based formulation. Demographic fluctuations are therefore naturally accounted and their effects are studied using a diffusion theory complemented by numerical simulations. It is thus shown that persistent erratic oscillations (quasi-cycles) of large amplitude emerge from a noise-induced resonance phenomenon. We also analytically and numerically compute the average escape time necessary to reach a (quasi-)cycle on which the system oscillates at a given amplitude., Comment: 25 pages, 9 figures. To appear in the Journal of Theoretical Biology

We study the cyclic dominance of three species in two-dimensional constrained Newman-Watts networks with a four-state variant of the rock-paper-scissors game. By limiting the maximal connection distance $R_{max}$ in Newman-Watts networks with the long-rang connection probability $p$, we depict more realistically the stochastic interactions among species within ecosystems. When we fix mobility and vary the value of $p$ or $R_{max}$, the Monte Carlo simulations show that the spiral waves grow in size, and the system becomes unstable and biodiversity is lost with increasing $p$ or $R_{max}$. These results are similar to recent results of Reichenbach \textit{et al.} [Nature (London) \textbf{448}, 1046 (2007)], in which they increase the mobility only without including long-range interactions. We compared extinctions with or without long-range connections and computed spatial correlation functions and correlation length. We conclude that long-range connections could improve the mobility of species, drastically changing their crossover to extinction and making the system more unstable., Comment: 6 pages, 7 figures

In this paper a reliable method is found for approximating the value of the Bivariate Normal Correlation Coefficient, rho, given values of the joint probability and the normal deviates, h and k, or the related areas. This technique finds useful application in the computation of the tetrachoric correlation coefficient, r, when the underlying distributions may be assumed to be normal. (Author)

We implemented three recently proposed approaches to the identification of overlapping and hierarchical substructures in graphs and applied the corresponding algorithms to a network of 492 information-science papers coupled via their cited sources. The thematic substructures obtained and overlaps produced by the three hierarchical cluster algorithms were compared to a content-based categorisation, which we based on the interpretation of titles and keywords. We defined sets of papers dealing with three topics located on different levels of aggregation: h-index, webometrics, and bibliometrics. We identified these topics with branches in the dendrograms produced by the three cluster algorithms and compared the overlapping topics they detected with one another and with the three pre-defined paper sets. We discuss the advantages and drawbacks of applying the three approaches to paper networks in research fields., Comment: 18 pages, 9 figures

Purpose Preference-based measures are essential for producing quality-adjusted life years (QALYs) that are widely used for economic evaluations. In the absence of such measures, mapping algorithms can be applied to estimate utilities from disease-specific measures. This paper aims to develop mapping algorithms between the MacNew Heart Disease Quality of Life Questionnaire (MacNew) instrument and the English and the US-based EQ-5D-5L value sets. Methods Individuals with heart disease were recruited from six countries: Australia, Canada, Germany, Norway, UK and the US in 2011/12. Both parametric and non-parametric statistical techniques were applied to estimate mapping algorithms that predict utilities for MacNew scores from EQ-5D-5L value sets. The optimal algorithm for each country-specific value set was primarily selected based on root mean square error (RMSE), mean absolute error (MAE), concordance correlation coefficient (CCC), and r-squared. Leave-one-out cross-validation was conducted to test the generalizability of each model. Results For both the English and the US value sets, the one-inflated beta regression model consistently performed best in terms of all criteria. Similar results were observed for the cross-validation results. The preferred model explained 59 and 60% for the English and the US value set, respectively. Linear equating provided predicted values that were equivalent to observed values. Conclusions The preferred mapping function enables to predict utilities for MacNew data from the EQ-5D-5L value sets recently developed in England and the US with better accuracy. This allows studies, which have included the MacNew to be used in cost-utility analyses and thus, the comparison of services with interventions across the health system.

Purpose Structure-guided total variation is a recently introduced prior that allows reconstruction of images using knowledge of the location and orientation of edges in a reference image. In this work, we demonstrate the advantages of a variant of structure-guided total variation known as directional total variation (DTV), over traditional total variation (TV), in the context of compressed-sensing reconstruction and super-resolution. Methods We compared TV and DTV in retrospectively undersampled ex vivo diffusion tensor imaging and diffusion spectrum imaging data from healthy, sham, and hypertrophic rat hearts. Results In compressed sensing at an undersampling factor of 8, the RMS error of mean diffusivity and fractional anisotropy relative to the fully sampled ground truth were 44% and 20% lower in DTV compared with TV. In super-resolution, these values were 29% and 14%, respectively. Similarly, we observed improvements in helix angle, transverse angle, sheetlet elevation, and sheetlet azimuth. The RMS error of the diffusion kurtosis in the undersampled data relative to the ground truth was uniformly lower (22% on average) with DTV compared to TV. Conclusion Acquiring one fully sampled non-diffusion-weighted image and 10 diffusion-weighted images at 8× undersampling would result in an 80% net reduction in data needed. We demonstrate efficacy of the DTV algorithm over TV in reducing data sampling requirements, which can be translated into higher apparent resolution and potentially shorter scan times. This method would be equally applicable in diffusion MRI applications outside the heart.

The aim of this study was to develop an automated segmentation approach for small gross tumor volumes (GTVs) in 3D planning computed tomography (CT) images using dense V-networks (DVNs) that offer more advantages in segmenting smaller structures than conventional V-networks. Regions of interest (ROI) with dimensions of 50 × 50 × 6–72 pixels in the planning CT images were cropped based on the GTV centroids when applying stereotactic body radiotherapy (SBRT) to patients. Segmentation accuracy of GTV contours for 192 lung cancer patients [with the following tumor types: 118 solid, 53 part-solid types and 21 pure ground-glass opacity (pure GGO)], who underwent SBRT, were evaluated based on a 10-fold cross-validation test using Dice’s similarity coefficient (DSC) and Hausdorff distance (HD). For each case, 11 segmented GTVs consisting of three single outputs, four logical AND outputs, and four logical OR outputs from combinations of two or three outputs from DVNs were obtained by three runs with different initial weights. The AND output (combination of three outputs) achieved the highest values of average 3D-DSC (0.832 ± 0.074) and HD (4.57 ± 2.44 mm). The average 3D DSCs from the AND output for solid, part-solid and pure GGO types were 0.838 ± 0.074, 0.822 ± 0.078 and 0.819 ± 0.059, respectively. This study suggests that the proposed approach could be useful in segmenting GTVs for planning lung cancer SBRT.

Purpose Multi-phase PCASL has been proposed as a means to achieve accurate perfusion quantification that is robust to imperfect shim in the labelling plane. However, there exists a bias in the estimation process that is a function of noise in the data. In this work, this bias is characterised and then addressed in animal and human data. Methods The proposed algorithm to overcome bias uses the initial biased voxel-wise estimate of phase tracking error to cluster regions with different off-resonance phase shifts, from which a high-SNR estimate of regional phase offset is derived. Simulations were used to predict the bias expected at typical SNR. Multi-phase PCASL in three rat strains (n = 21) at 9.4 T was considered, along with 20 human subjects previously imaged using ASL at 3 T. The algorithm was extended to include estimation of arterial blood flow velocity. Results Based on simulations, a perfusion estimation bias of 6–8% was expected using 8-phase data at typical SNR. This bias was eliminated when a high-precision estimate of phase error was available. In the preclinical data, the bias-corrected measure of perfusion (107±14 mL/100g/min) was lower than the standard analysis (116±14 mL/100g/min), corresponding to a mean observed bias across strains of 8.0%. In the human data, bias correction resulted in a 15% decrease in the estimate of perfusion. Conclusion Using a retrospective algorithmic approach, it was possible to exploit common information found in multiple voxels within a whole region of the brain, offering superior SNR and thus overcoming the bias in perfusion quantification from multi-phase PCASL.

Motivated by the current COVID-19 health crisis, we consider data analysis for quantitative polymerase chain-reaction (qPCR) measurements. We derive a theoretical result specifying the conditions under which all qPCR amplification curves (including their plateau phases) are identical up to an affine transformation, i.e. a multiplicative factor and horizontal shift. We use this result to develop a data analysis procedure for determining when an amplification curve exhibits characteristics of a true signal. The main idea behind this approach is to invoke a criterion based on constrained optimization that assesses when a measurement signal can be mapped to a master reference curve. We demonstrate that this approach: (i) can decrease the fluorescence detection threshold by up to a decade; and (ii) simultaneously improve confidence in interpretations of late-cycle amplification curves. Moreover, we demonstrate that the master curve is transferable reference data that can harmonize analyses between different labs and across several years. Application to reverse-transcriptase qPCR measurements of a SARS-CoV-2 RNA construct points to the usefulness of this approach for improving confidence and reducing limits of detection in diagnostic testing of emerging diseases. Graphical Abstract Left: a collection of qPCR amplification curves. Right: Example of data collapse after affine transformation.

Background Margin reflex distance 1 (MRD1), margin reflex distance 2 (MRD2), and levator muscle function (LF) are crucial metrics for ptosis evaluation and management. However, manual measurements of MRD1, MRD2, and LF are time-consuming, subjective, and prone to human error. Smartphone-based artificial intelligence (AI) image processing is a potential solution to overcome these limitations. Objective We propose the first smartphone-based AI-assisted image processing algorithm for MRD1, MRD2, and LF measurements. Methods This observational study included 822 eyes of 411 volunteers aged over 18 years from August 1, 2020, to April 30, 2021. Six orbital photographs (bilateral primary gaze, up-gaze, and down-gaze) were taken using a smartphone (iPhone 11 Pro Max). The gold-standard measurements and normalized eye photographs were obtained from these orbital photographs and compiled using AI-assisted software to create MRD1, MRD2, and LF models. Results The Pearson correlation coefficients between the gold-standard measurements and the predicted values obtained with the MRD1 and MRD2 models were excellent (r=0.91 and 0.88, respectively) and that obtained with the LF model was good (r=0.73). The intraclass correlation coefficient demonstrated excellent agreement between the gold-standard measurements and the values predicted by the MRD1 and MRD2 models (0.90 and 0.84, respectively), and substantial agreement with the LF model (0.69). The mean absolute errors were 0.35 mm, 0.37 mm, and 1.06 mm for the MRD1, MRD2, and LF models, respectively. The 95% limits of agreement were –0.94 to 0.94 mm for the MRD1 model, –0.92 to 1.03 mm for the MRD2 model, and –0.63 to 2.53 mm for the LF model. Conclusions We developed the first smartphone-based AI-assisted image processing algorithm for eyelid measurements. MRD1, MRD2, and LF measures can be taken in a quick, objective, and convenient manner. Furthermore, by using a smartphone, the examiner can check these measurements anywhere and at any time, which facilitates data collection.

Purpose To obtain a fast and robust fat-water separation with simultaneous estimation of water T1 , fat T1 , and fat fraction maps. Methods We modified an MR fingerprinting (MRF) framework to use a single dictionary combination of a water and fat dictionary. A variable TE acquisition pattern with maximum TE = 4.8 ms was used to increase the fat-water separability. Radiofrequency (RF) spoiling was used to reduce the size of the dictionary by reducing T2 sensitivity. The technique was compared both in vitro and in vivo to an MRF method that incorporated 3-point Dixon (DIXON MRF), as well as Cartesian IDEAL with different acquisition parameters. Results The proposed dictionary-based fat-water separation technique (DBFW MRF) successfully provided fat fraction, water, and fat T1 , B0 , and B1+ maps both in vitro and in vivo. The fat fraction and water T1 values obtained with DBFW MRF show excellent agreement with DIXON MRF as well as with the reference values obtained using a Cartesian IDEAL with a long TR (concordance correlation coefficient: 0.97/0.99 for fat fraction-water T1 ). Whereas fat fraction values with Cartesian IDEAL were degraded in the presence of T1 saturation, MRF methods successfully estimated and accounted for T1 in the fat fraction estimates. Conclusion The DBFW MRF technique can successfully provide T1 and fat fraction quantification in under 20 s per slice, intrinsically correcting T1 biases typical of fast Dixon techniques. These features could improve the diagnostic quality and use of images in presence of fat.

Significance: The Mueller matrix decomposition method is widely used for the analysis of biological samples. However, its presumed sequential appearance of the basic optical effects (e.g., dichroism, retardance, and depolarization) limits its accuracy and application. Aim: An approach is proposed for detecting and classifying human melanoma and non-melanoma skin cancer lesions based on the characteristics of the Mueller matrix elements and a random forest (RF) algorithm. Approach: In the proposal technique, 669 data points corresponding to the 16 elements of the Mueller matrices obtained from 32 tissue samples with squamous cell carcinoma (SCC), basal cell carcinoma (BCC), melanoma, and normal features are input into an RF classifier as predictors. Results: The results show that the proposed model yields an average precision of 93%. Furthermore, the classification results show that for biological tissues, the circular polarization properties (i.e., elements m44, m34, m24, and m14 of the Mueller matrix) dominate the linear polarization properties (i.e., elements m13, m31, m22, and m41 of the Mueller matrix) in determining the classification outcome of the trained classifier. Conclusions: Overall, our study provides a simple, accurate, and cost-effective solution for developing a technique for classification and diagnosis of human skin cancer.

We evaluated the performance of a newly developed three-dimensional (3D) model-based global-to-local registration of multiple organs, by comparing it with a 3D model-based global registration in the prostate region. This study included 220 prostate cancer patients who underwent intensity-modulated radiotherapy or volumetric-modulated arc therapy. Our registration proceeded sequentially, i.e. global registration including affine and piece-wise affine transformation followed by local registration. As a local registration, Laplacian-based and finite element method-based registration was implemented in Algorithm A and B, respectively. Algorithm C was for global registration alone. The template models for the prostate, seminal vesicles, rectum and bladder were constructed from the first 20 patients, and then three different registrations were performed on these organs for the remaining 200 patients, to assess registration accuracy. The 75th percentile Hausdorff distance was 1 mm in Algorithm B, except for the prostate; and 3.9 mm for the prostate and >7.8 mm for other organs in Algorithm C. The median computation time to complete registration was

Purpose Multi-echo spin-echo (MSE) transverse relaxometry mapping using multi-component models is used to study disease activity in neuromuscular disease by assessing the T-2 of the myocytic component (T-2water). Current extended phase graph algorithms are not optimized for fat fractions above 50% and the effects of inaccuracies in the T-2fat calibration remain unexplored. Hence, we aimed to improve the performance of extended phase graph fitting methods over a large range of fat fractions, by including the slice-selection flip angle profile, a through-plane chemical-shift displacement correction, and optimized calibration of T-2fat.Methods Simulation experiments were used to study the influence of the slice flip-angle profile with chemical-shift and T-2fat estimations. Next, in vivo data from four neuromuscular disease cohorts were studied for different T-2fat calibration methods and T-2water estimations.Results Excluding slice flip-angle profiles or chemical-shift displacement resulted in a bias in T-2water up to 10 ms. Furthermore, a wrongly calibrated T-2fat caused a bias of up to 4 ms in T-2water. For the in vivo data, one-component calibration led to a lower T-2fat compared with a two-component method, and T-2water decreased with increasing fat fractions.Conclusion In vivo data showed a decline in T-2water for increasing fat fractions, which has important implications for clinical studies, especially in multicenter settings. We recommend using an extended phase graph-based model for fitting T-2water from MSE sequences with two-component T-2fat calibration. Moreover, we recommend including the slice flip-angle profile in the model with correction for through-plane chemical-shift displacements.

Quantitative metrics in lung computed tomography (CT) images have been widely used, often without a clear connection with physiology. This work proposes a patient-independent model for the estimation of well-aerated volume of lungs in CT images (WAVE). A Gaussian fit, with mean (Mu.f) and width (Sigma.f) values, was applied to the lower CT histogram data points of the lung to provide the estimation of the well-aerated lung volume (WAVE.f). Independence from CT reconstruction parameters and respiratory cycle was analysed using healthy lung CT images and 4DCT acquisitions. The Gaussian metrics and first order radiomic features calculated for a third cohort of COVID-19 patients were compared with those relative to healthy lungs. Each lung was further segmented in 24 subregions and a new biomarker derived from Gaussian fit parameter Mu.f was proposed to represent the local density changes. WAVE.f resulted independent from the respiratory motion in 80% of the cases. Differences of 1%, 2% and up to 14% resulted comparing a moderate iterative strength and FBP algorithm, 1 and 3 mm of slice thickness and different reconstruction kernel. Healthy subjects were significantly different from COVID-19 patients for all the metrics calculated. Graphical representation of the local biomarker provides spatial and quantitative information in a single 2D picture. Unlike other metrics based on fixed histogram thresholds, this model is able to consider the inter-and intra-subject variability. In addition, it defines a local biomarker to quantify the severity of the disease, independently of the observer., 31 pages, 7 figures, accepted in EJMP