Demarcating Z-line and Gastric Folds Boundary Based on the Segmentation of the Lower Esophageal Sphincter Images.

Background and Objective: The endoscopic diagnosis of pathological changes in the gastroesophageal junction including esophagitis and Barrett's mucosa is based on the visual detection of two boundaries: mucosal color change between esophagus and stomach, and top endpoint of gastric folds. The presence and pattern of mucosal breaks in the gastroesophageal mucosal junction (Z line) classify esophagitis in patients and the distance between the two boundaries points to the possible columnar lined epithelium. Since visual detection may suffer from intra- and interobserver variability, our objective was to define the boundaries automatically based on image processing algorithms, which may enable us to measure the detentions of changes in future studies. Methods: To demarcate the Z-line, first the artifacts of endoscopy images are eliminated. In the second step, using SUSAN edge detector, Mahalanobis distance criteria, and Gabor filter bank, an initial contour is estimated for the Z-line. Using region-based active contours, this initial contour converges to the Z-line. Finally, by applying morphological operators and Gabor Filter Bank to the region inside of the Z-line, gastric folds are segmented. Results: To evaluate the results, a database consisting of 50 images and their ground truths were collected. The average dice coefficient and mean square error of Z-line segmentation were 0.93 and 3.3, respectively. Furthermore, the average boundary distance criteria are 12.3 pixels. In addition, two other criteria that compare the segmentation of folds with several ground truths, i.e., Sweet-Spot Coverage and Jaccard Index for Golden Standard, are 0.90 and 0.84, respectively. Conclusions: Considering the results, automatic segmentation of Z-line and gastric folds are matched to the ground truths with appropriate accuracy.

CNV-Net: Segmentation, Classification and Activity Score Measurement of Choroidal Neovascularization (CNV) Using Optical Coherence Tomography Angiography (OCTA).

This paper aims to present an artificial intelligence-based algorithm for the automated segmentation of Choroidal Neovascularization (CNV) areas and to identify the presence or absence of CNV activity criteria (branching, peripheral arcade, dark halo, shape, loop and anastomoses) in OCTA images. Methods: This retrospective and cross-sectional study includes 130 OCTA images from 101 patients with treatment-naïve CNV. At baseline, OCTA volumes of 6 × 6 mm2 were obtained to develop an AI-based algorithm to evaluate the CNV activity based on five activity criteria, including tiny branching vessels, anastomoses and loops, peripheral arcades, and perilesional hypointense halos. The proposed algorithm comprises two steps. The first block includes the pre-processing and segmentation of CNVs in OCTA images using a modified U-Net network. The second block consists of five binary classification networks, each implemented with various models from scratch, and using transfer learning from pre-trained networks. Results: The proposed segmentation network yielded an averaged Dice coefficient of 0.86. The individual classifiers corresponding to the five activity criteria (branch, peripheral arcade, dark halo, shape, loop, and anastomoses) showed accuracies of 0.84, 0.81, 0.86, 0.85, and 0.82, respectively. The AI-based algorithm potentially allows the reliable detection and segmentation of CNV from OCTA alone, without the need for imaging with contrast agents. The evaluation of the activity criteria in CNV lesions obtains acceptable results, and this algorithm could enable the objective, repeatable assessment of CNV features.

Correction to: Spike Sorting of Non-Stationary Data in Successive Intervals Based on Dirichlet Process Mixtures.

[This corrects the article DOI: 10.1007/s11571-022-09781-7.].

Spike Sorting of Non-Stationary Data in Successive Intervals Based on Dirichlet Process Mixtures.

This paper proposes a new automatic method for spike sorting and tracking non-stationary data based on the Dirichlet Process Mixture (DPM). Data is divided into non-overlapping intervals and mixtures are applied to individual frames rather than to the whole data. In this paper, we have used the information of the previous frame to estimate the cluster parameters of the current interval. Specifically, the means of the clusters in the previous frame are used for estimating the cluster means of the current one, and other parameters are estimated via noninformative priors. The proposed method is capable to track variations in size, shape, or location of clusters as well as detecting the appearance and disappearance of them. We present results in two-dimensional space of first and second principal components (PC1-PC2), but any other feature extraction method leading to the ability of modeling spikes with Normal or t-Student distributions can also be applied. Application of this approach to simulated data and the recordings from anesthetized rat hippocampus confirms its superior performance in comparison to a standard DPM that uses no information from previous frames.

Tropical cyclone simulations over Bangladesh at convection permitting 4.4 km & 1.5 km resolution.

High resolution simulations at 4.4 km and 1.5 km resolution have been performed for 12 historical tropical cyclones impacting Bangladesh. We use the European Centre for Medium-Range Weather Forecasting 5th generation Re-Analysis (ERA5) to provide a 9-member ensemble of initial and boundary conditions for the regional configuration of the Met Office Unified Model. The simulations are compared to the original ERA5 data and the International Best Track Archive for Climate Stewardship (IBTrACS) tropical cyclone database for wind speed, gust speed and mean sea-level pressure. The 4.4 km simulations show a typical increase in peak gust speed of 41 to 118 knots relative to ERA5, and a deepening of minimum mean sea-level pressure of up to -27 hPa, relative to ERA5 and IBTrACS data. The downscaled simulations compare more favourably with IBTrACS data than the ERA5 data suggesting tropical cyclone hazards in the ERA5 deterministic output may be underestimated. The dataset is freely available from https://doi.org/10.5281/zenodo.3600201 .

Graph Embedding Using Frequency Filtering.

The target of graph embedding is to embed graphs in vector space such that the embedded feature vectors follow the differences and similarities of the source graphs. In this paper, a novel method named Frequency Filtering Embedding (FFE) is proposed which uses graph Fourier transform and Frequency filtering as a graph Fourier domain operator for graph feature extraction. Frequency filtering amplifies or attenuates selected frequencies using appropriate filter functions. Here, heat, anti-heat, part-sine and identity filter sets are proposed as the filter functions. A generalized version of FFE named GeFFE is also proposed by defining pseudo-Fourier operators. This method can be considered as a general framework for formulating some previously defined invariants in other works by choosing a suitable filter bank and defining suitable pseudo-Fourier operators. This flexibility empowers GeFFE to adapt itself to the properties of each graph dataset unlike the previous spectral embedding methods and leads to superior classification accuracy relative to the others. Utilizing the proposed part-sine filter set, which its members filter different parts of the spectrum in turn, improves the classification accuracy of GeFFE method. Additionally, GeFFE resolves the cospectrality problem entirely in tested datasets.

Parameter Estimation in Synaptic Coupling Model Using a Point Process Modeling Framework.

Biophysical models are widely used to characterize temporal dynamics of the brain networks on different topological and spatial scales. In parallel, the state-space modeling framework with point process observations has been successfully applied in characterizing spiking activity of neuronal ensembles in response to different dynamical covariates. Parameter estimation in biophysical models is generally done heuristically, which hampers their applicability and interpretability. Heuristic parameter estimation becomes an intractable problem when the number of model parameters grows. Here, we propose an algorithm for estimating biophysical model parameters using point-process models and a state-space framework. The framework provides methods for parameter estimation as well as model validation. We demonstrate the application of this methodology to the problem of estimating the parameters of a dynamic synapse model. We generate simulation data for the dynamic synapse across a range of parameters values and assess the estimation accuracy of our method using a combination of goodness-of-fit measures. The proposed methodology can be applied broadly to parameter estimation problems across a broad range of biophysical models, including Hodgkin-Huxley models and network models.

Fully Automated Segmentation of Fluid/Cyst Regions in Optical Coherence Tomography Images With Diabetic Macular Edema Using Neutrosophic Sets and Graph Algorithms.

This paper presents a fully automated algorithm to segment fluid-associated (fluid-filled) and cyst regions in optical coherence tomography (OCT) retina images of subjects with diabetic macular edema. The OCT image is segmented using a novel neutrosophic transformation and a graph-based shortest path method. In neutrosophic domain, an image is transformed into three sets: (true), (indeterminate) that represents noise, and (false). This paper makes four key contributions. First, a new method is introduced to compute the indeterminacy set , and a new -correction operation is introduced to compute the set in neutrosophic domain. Second, a graph shortest-path method is applied in neutrosophic domain to segment the inner limiting membrane and the retinal pigment epithelium as regions of interest (ROI) and outer plexiform layer and inner segment myeloid as middle layers using a novel definition of the edge weights . Third, a new cost function for cluster-based fluid/cyst segmentation in ROI is presented which also includes a novel approach in estimating the number of clusters in an automated manner. Fourth, the final fluid regions are achieved by ignoring very small regions and the regions between middle layers. The proposed method is evaluated using two publicly available datasets: Duke, Optima, and a third local dataset from the UMN clinic which is available online. The proposed algorithm outperforms the previously proposed Duke algorithm by 8% with respect to the dice coefficient and by 5% with respect to precision on the Duke dataset, while achieving about the same sensitivity. Also, the proposed algorithm outperforms a prior method for Optima dataset by 6%, 22%, and 23% with respect to the dice coefficient, sensitivity, and precision, respectively. Finally, the proposed algorithm also achieves sensitivity of 67.3%, 88.8%, and 76.7%, for the Duke, Optima, and the university of minnesota (UMN) datasets, respectively.

Fully-automated segmentation of fluid regions in exudative age-related macular degeneration subjects: Kernel graph cut in neutrosophic domain.

A fully-automated method based on graph shortest path, graph cut and neutrosophic (NS) sets is presented for fluid segmentation in OCT volumes for exudative age related macular degeneration (EAMD) subjects. The proposed method includes three main steps: 1) The inner limiting membrane (ILM) and the retinal pigment epithelium (RPE) layers are segmented using proposed methods based on graph shortest path in NS domain. A flattened RPE boundary is calculated such that all three types of fluid regions, intra-retinal, sub-retinal and sub-RPE, are located above it. 2) Seed points for fluid (object) and tissue (background) are initialized for graph cut by the proposed automated method. 3) A new cost function is proposed in kernel space, and is minimized with max-flow/min-cut algorithms, leading to a binary segmentation. Important properties of the proposed steps are proven and quantitative performance of each step is analyzed separately. The proposed method is evaluated using a publicly available dataset referred as Optima and a local dataset from the UMN clinic. For fluid segmentation in 2D individual slices, the proposed method outperforms the previously proposed methods by 18%, 21% with respect to the dice coefficient and sensitivity, respectively, on the Optima dataset, and by 16%, 11% and 12% with respect to the dice coefficient, sensitivity and precision, respectively, on the local UMN dataset. Finally, for 3D fluid volume segmentation, the proposed method achieves true positive rate (TPR) and false positive rate (FPR) of 90% and 0.74%, respectively, with a correlation of 95% between automated and expert manual segmentations using linear regression analysis.

Analysis of Fundus Fluorescein Angiogram Based on the Hessian Matrix of Directional Curvelet Sub-bands and Distance Regularized Level Set Evolution.

This paper presents a new procedure for automatic extraction of the blood vessels and optic disk (OD) in fundus fluorescein angiogram (FFA). In order to extract blood vessel centerlines, the algorithm of vessel extraction starts with the analysis of directional images resulting from sub-bands of fast discrete curvelet transform (FDCT) in the similar directions and different scales. For this purpose, each directional image is processed by using information of the first order derivative and eigenvalues obtained from the Hessian matrix. The final vessel segmentation is obtained using a simple region growing algorithm iteratively, which merges centerline images with the contents of images resulting from modified top-hat transform followed by bit plane slicing. After extracting blood vessels from FFA image, candidates regions for OD are enhanced by removing blood vessels from the FFA image, using multi-structure elements morphology, and modification of FDCT coefficients. Then, canny edge detector and Hough transform are applied to the reconstructed image to extract the boundary of candidate regions. At the next step, the information of the main arc of the retinal vessels surrounding the OD region is used to extract the actual location of the OD. Finally, the OD boundary is detected by applying distance regularized level set evolution. The proposed method was tested on the FFA images from angiography unit of Isfahan Feiz Hospital, containing 70 FFA images from different diabetic retinopathy stages. The experimental results show the accuracy more than 93% for vessel segmentation and more than 87% for OD boundary extraction.

Automated acoustic analysis in detection of spontaneous swallows in Parkinson's disease.

Acoustic monitoring of swallow frequency has become important as the frequency of spontaneous swallowing can be an index for dysphagia and related complications. In addition, it can be employed as an objective quantification of ingestive behavior. Commonly, swallowing complications are manually detected using videofluoroscopy recordings, which require expensive equipment and exposure to radiation. In this study, a noninvasive automated technique is proposed that uses breath and swallowing recordings obtained via a microphone located over the laryngopharynx. Nonlinear diffusion filters were used in which a scale-space decomposition of recorded sound at different levels extract swallows from breath sounds and artifacts. This technique was compared to manual detection of swallows using acoustic signals on a sample of 34 subjects with Parkinson's disease. A speech language pathologist identified five subjects who showed aspiration during the videofluoroscopic swallowing study. The proposed automated method identified swallows with a sensitivity of 86.67 %, a specificity of 77.50 %, and an accuracy of 82.35 %. These results indicate the validity of automated acoustic recognition of swallowing as a fast and efficient approach to objectively estimate spontaneous swallow frequency.

Vesselness-guided Active Contour: A Coronary Vessel Extraction Method.

Vessel extraction is a critical task in clinical practice. In this paper, we propose a new approach for vessel extraction using an active contour model by defining a novel vesselness-based term, based on accurate analysis of the vessel structure in the image. To achieve the novel term, a simple and fast directional filter bank is proposed, which does not employ down sampling and resampling used in earlier versions of directional filter banks. The proposed model not only preserves the performance of the existing models on images with intensity inhomogeneity, but also overcomes their inability both to segment low contrast vessels and to omit non-vessel structures. Experimental results for synthetic images and coronary X-ray angiograms show desirable performance of our model.

On the quantity and composition of floating plastic debris entering and leaving the Tamar Estuary, Southwest England.

The majority of plastic debris found in the marine environment has land based sources and rivers are considered an important medium for transfer of this debris. Here we report on the quantity and composition of floating plastic debris collected from surface waters of the Tamar Estuary. This represents the first study of riverine transport of floating plastic debris into European waters during different tidal regimes. Plastics were found in a variety of forms and sizes and microplastics (<5 mm) comprised 82% of the debris. The most abundant types of plastic were Polyethylene (40%), Polystyrene (25%) and Polypropylene (19%). There was a significant difference in size frequency distribution between the spring and neap tides with more fragments of larger size observed during spring tides. While it is clear that debris has accumulated on shorelines within the estuary, during our study this river did not identify as a net source or sink.

Segmentation of dermoscopy images using wavelet networks.

This paper introduces a new approach for the segmentation of skin lesions in dermoscopic images based on wavelet network (WN). The WN presented here is a member of fixed-grid WNs that is formed with no need of training. In this WN, after formation of wavelet lattice, determining shift and scale parameters of wavelets with two screening stage and selecting effective wavelets, orthogonal least squares algorithm is used to calculate the network weights and to optimize the network structure. The existence of two stages of screening increases globality of the wavelet lattice and provides a better estimation of the function especially for larger scales. R, G, and B values of a dermoscopy image are considered as the network inputs and the network structure formation. Then, the image is segmented and the skin lesions exact boundary is determined accordingly. The segmentation algorithm were applied to 30 dermoscopic images and evaluated with 11 different metrics, using the segmentation result obtained by a skilled pathologist as the ground truth. Experimental results show that our method acts more effectively in comparison with some modern techniques that have been successfully used in many medical imaging problems.

A review of coronary vessel segmentation algorithms.

Coronary heart disease has been one of the main threats to human health. Coronary angiography is taken as the gold standard; for the assessment of coronary artery disease. However, sometimes, the images are difficult to visually interpret because of the crossing and overlapping of vessels in the angiogram. Vessel extraction from X-ray angiograms has been a challenging problem for several years. There are several problems in the extraction of vessels, including: weak contrast between the coronary arteries and the background, unknown and easily deformable shape of the vessel tree, and strong overlapping shadows of the bones. In this article we investigate the coronary vessel extraction and enhancement techniques, and present capabilities of the most important algorithms concerning coronary vessel segmentation.

An automated method for segmentation of epithelial cervical cells in images of ThinPrep.

We present an automated method for segmentation of epithelial cells in images taken from ThinPrep scenes by a digital camera in a cytology lab. The method covers both steps of localization of cell objects in low resolution and detection of cytoplasm and nucleus boundary in high resolution. The underlying method makes use of geometric active contours as a powerful tool of segmentation. We also provide the analysis of the connected cells. For this purpose an automatic circular decomposition method is incorporated and adapted to the application by changing its segmentation condition. The results are evaluated numerically and compared with those of previous work in literature.