基于VT K 的D IC OM 医学图像三维重建研究

医学图像三维重建技术利用二维医学图像序列重建出三维模型,为医生提供直观、全面、准确的病灶和正常组织信息．是当今医学影像领域研究的热点之一。VTK是国际上广泛应用的可视化工具包,具有优秀的架构和运行机制。本文研究了DICOM3．0标准,提出了正确解读DICOM医学图像的方法;深入VTK内部机制,解决了VTK和DICOM医学图像读取模块间的数据接口问题;在三维重建过程中,为解决数据量巨大、成像时间漫长、阶梯效应、交互性不强等问题,重点剖析了VTK的数据处理机制。并给出相关优化方法。实验结果表明本文提出的解决方案和优化方法实用可靠,为进一步开发医学三雏图形系统打下了基础。     

早期胚胎全场OCT图像的边界去模糊研究

全场光学相干层析成像技术(全场OCT)是研究早期胚胎形态发育的最理想成像设备,然而所采集图像难免受噪声干扰.这些噪声可模糊早期胚胎内不同组织结构的边界,从而给基于图像边界的结构划分带来干扰.为解决这一问题,本文运用中值滤波、维纳滤波、各向异性扩散算法处理全场OCT获得的早期胚胎图像,并运用信噪比、均方误差、峰值信噪比和边缘保留等指标评价图像处理效果.结果表明:经各向异性扩散算法处理的早期胚胎图像,可完整地保留原始图像信息,且边界最清晰,视觉效果最好.     

On the Mathematical Reconstruction of Two Dimensional Plants

A set of 10, chosen medicinal plants (some of them with a reputation as remedies for tuberculosis) has been investigated through Partitioned Iterated Function Systems-Semi Fractals with Angle (PIFS-SFA) coding, Lempel, Ziv, Welch with quantization and noise (LZW-QN) compression, and surface density statistics (f(α)-SDS) discrimination techniques. The final outcomes of this quantitative analysis were, firstly: the linear ordering of the plants in question accompanied by the hope that it reflects their medical significance, secondly: the mathematical representation of each of the plants, and thirdly: the impressive compression achieved, leading to remarkable computer memory saving, and still permitting successful pattern recognition i.e., proper identification of the plant from the compressed image.     

Image Alignment for Tomography Reconstruction from Synchrotron X-Ray Microscopic Images

A synchrotron X-ray microscope is a powerful imaging apparatus for taking high-resolution and high-contrast X-ray images of nanoscale objects. A sufficient number of X-ray projection images from different angles is required for constructing 3D volume images of an object. Because a synchrotron light source is immobile, a rotational object holder is required for tomography. At a resolution of 10 nm per pixel, the vibration of the holder caused by rotating the object cannot be disregarded if tomographic images are to be reconstructed accurately. This paper presents a computer method to compensate for the vibration of the rotational holder by aligning neighboring X-ray images. This alignment process involves two steps. The first step is to match the “projected feature points” in the sequence of images. The matched projected feature points in the - plane should form a set of sine-shaped loci. The second step is to fit the loci to a set of sine waves to compute the parameters required for alignment. The experimental results show that the proposed method outperforms two previously proposed methods, Xradia and SPIDER. The developed software system can be downloaded from the URL, http://www.cs.nctu.edu.tw/~chengchc/SCTA or http://goo.gl/s4AMx.     

Correlation between Lamina Cribrosa Tilt Angles,Myopia and Glaucoma Using OCT with a Wide Bandwidth Femtosecond Mode-Locked Laser

Purpose

To measure horizontal and vertical lamina cribrosa (LC) tilt angles and investigate associated factors using prototype optical coherence tomography (OCT) with a broad wavelength laser light source.

Design

Cross sectional study.

Methods

Twenty-eight no glaucoma eyes (from 15 subjects) and 25 glaucoma eyes (from 14 patients) were enrolled. A total of 300 optic nerve head B-scans were obtained in 10 µm steps and the inner edge of Bruch''s membrane opening (BMO) was identified as the reference plane. The vertical and horizontal angles between BMO line and approximate the best-fitting line for the surface of the LC were measured and potential associated factors were estimated with univariate and multivariate logistic regression analyses.

Results

The median (interquartile range) horizontal and vertical tilt angles were 7.10 (2.43–11.45) degrees and 4.15 (2.60–6.85) degrees in eyes without glaucoma and 8.50 (4.40–14.10) degrees and 9.30 (6.90–14.15) degrees in glaucoma eyes, respectively. The refractive errors had a statistically significant association with horizontal LC tilt angles (coefficients, −1.53 per diopter) and glaucoma had a significant correlation with vertical tilt angles (coefficients, 6.56) using multiple logistic regression analysis (p<0.001).

Conclusions

OCT allowed evaluation of the internal tilting of the LC compared with the BMO. The horizontal internal LC tilt angle was correlated with refractive errors, corresponding to myopic physiological changes, and vertical internal LC tilt was correlated with glaucoma, corresponding to glaucomatous pathological changes. These parameters have important implications for investigation of the correlation between myopia, glaucoma and LC morphological features.     

High-throughput Computer Method for 3D Neuronal Structure Reconstruction from the Image Stack of the Drosophila Brain and Its Applications

Drosophila melanogaster is a well-studied model organism, especially in the field of neurophysiology and neural circuits. The brain of the Drosophila is small but complex, and the image of a single neuron in the brain can be acquired using confocal microscopy. Analyzing the Drosophila brain is an ideal start to understanding the neural structure. The most fundamental task in studying the neural network of Drosophila is to reconstruct neuronal structures from image stacks. Although the fruit fly brain is small, it contains approximately 100 000 neurons. It is impossible to trace all the neurons manually. This study presents a high-throughput algorithm for reconstructing the neuronal structures from 3D image stacks collected by a laser scanning confocal microscope. The proposed method reconstructs the neuronal structure by applying the shortest path graph algorithm. The vertices in the graph are certain points on the 2D skeletons of the neuron in the slices. These points are close to the 3D centerlines of the neuron branches. The accuracy of the algorithm was verified using the DIADEM data set. This method has been adopted as part of the protocol of the FlyCircuit Database, and was successfully applied to process more than 16 000 neurons. This study also shows that further analysis based on the reconstruction results can be performed to gather more information on the neural network.     

Analysis of OCT1, OCT2 and OCT3 gene polymorphisms among Type 2 diabetes mellitus subjects in Indian ethnicity,Malaysia

BackgroundType 2 Diabetes mellitus (T2DM) is a chronic metabolic disorder. It is a major non-communicable disease affecting 463 million people globally in 2019 and is expected to be double to about 700 million by 2045. The majority are Asians with Indian ethnicity in Malaysia reported as the highest prevalence of T2DM. Cardiovascular disease, renal failure, blindness and neuropathy, as well as premature death are the known morbidity and mortality resulted from T2DM. T2DM is characterized by the dysfunctional insulin physiology that causes reduction of glucose transport into the cells which lead to hyperglycaemia. Hence, one of the important treatments is an oral antidiabetic drug that lowers the serum glucose level in patients with T2DM. This drug will be transported across cell membranes by organic cation transporters (OCT). Therefore, it is important to identify the OCT candidate gene polymorphisms related to T2DM especially among the Indian ethnicity in Malaysia.MethodsBlood samples were collected from 132 T2DM patients and 133 controls. Genotyping of OCT1 (rs628031), OCT2 (rs145450955), OCT3 (rs3088442 and rs2292334) was performed using (PCR-RFLP).ResultsNo association was observed for genotypic and allelic distributions in all the gene polymorphisms of OCT genes (P > 0.05). However, a logistic regression analysis stratified by gender in a dominant model showed a significant difference for OCT3 among males with T2DM (P = 0.006). Significant association was also observed for OCT3 when stratified to subjects aged > 45 years old (P = 0.009).ConclusionBased on these findings, the association of OCT3 (rs2292334) could be considered as a possible genetic risk factor for the development of T2DM among Indian males alone.     

Three-Dimensional Reconstruction from Reduced Sets of Very Noisy Images Acquired Following a Single-Axis Tilt Schema: Application of a New Three-Dimensional Reconstruction Algorithm and Objective Comparison with Weighted Backprojection

In this work we propose a reconstruction algorithm (ART with blobs) that has not been previously used in electron Tomography and we compare it with the standard method in the field (weighted back projection, WBP). We assume that only a limited set of very noisy images, collected around a single axis tilt, is available; which is a typical situation in Electron Tomography. In general, the reconstruction problem is underdetermined (due to the limited number of projections) and the data are inconsistent (due to the high level of noise). The evaluation of the results is performed in a rigorous way by a task-oriented approach which makes use of numerical observers. ART with blobs outperforms WBP for a number of key tasks. Results are presented both for simplified line integral data and for realistic simulations of macromolecular structures embedded in amorphous ice.     

基于VTK和MFC的医学图像三维重建研究与实现

VTK是医学可视化领域的主流工具,MFC是Windows平台下的应用程序框架。尝试将两者进行结合编程,以实现二维医学图像的三维重建。实现医学图像三维重建的主要方法是面绘制和体绘制。将利用多组医学图像数据进行三维重建研究,其中面绘制用移动立方体法,体绘制用光线投射法、最大密度投影法和合成体绘制法实现。最后比较两种绘制技术的结果并讨论了它们的特点。结果表明,VTK作为一种图像处理和三维可视化工具其功能是十分强大的。     

Reconstruction of 3D structures from protein contact maps

The prediction of the protein tertiary structure from solely its residue sequence (the so called Protein Folding Problem) is one of the most challenging problems in Structural Bioinformatics. We focus on the protein residue contact map. When this map is assigned it is possible to reconstruct the 3D structure of the protein backbone. The general problem of recovering a set of 3D coordinates consistent with some given contact map is known as a unit-disk-graph realization problem and it has been recently proven to be NP-Hard. In this paper we describe a heuristic method (COMAR) that is able to reconstruct with an unprecedented rate (3-15 seconds) a 3D model that exactly matches the target contact map of a protein. Working with a non-redundant set of 1760 proteins, we find that the scoring efficiency of finding a 3D model very close to the protein native structure depends on the threshold value adopted to compute the protein residue contact map. Contact maps whose threshold values range from 10 to 18 Ångstroms allow reconstructing 3D models that are very similar to the proteins native structure.     

Continuous Stroke Volume Estimation from Aortic Pressure Using Zero Dimensional Cardiovascular Model: Proof of Concept Study from Porcine Experiments

Introduction

Accurate, continuous, left ventricular stroke volume (SV) measurements can convey large amounts of information about patient hemodynamic status and response to therapy. However, direct measurements are highly invasive in clinical practice, and current procedures for estimating SV require specialized devices and significant approximation.

Method

This study investigates the accuracy of a three element Windkessel model combined with an aortic pressure waveform to estimate SV. Aortic pressure is separated into two components capturing; 1) resistance and compliance, 2) characteristic impedance. This separation provides model-element relationships enabling SV to be estimated while requiring only one of the three element values to be known or estimated. Beat-to-beat SV estimation was performed using population-representative optimal values for each model element. This method was validated using measured SV data from porcine experiments (N = 3 female Pietrain pigs, 29–37 kg) in which both ventricular volume and aortic pressure waveforms were measured simultaneously.

Results

The median difference between measured SV from left ventricle (LV) output and estimated SV was 0.6 ml with a 90% range (5^th–95^th percentile) −12.4 ml–14.3 ml. During periods when changes in SV were induced, cross correlations in between estimated and measured SV were above R = 0.65 for all cases.

Conclusion

The method presented demonstrates that the magnitude and trends of SV can be accurately estimated from pressure waveforms alone, without the need for identification of complex physiological metrics where strength of correlations may vary significantly from patient to patient.     

MR Image Reconstruction Using Block Matching and Adaptive Kernel Methods

An approach to Magnetic Resonance (MR) image reconstruction from undersampled data is proposed. Undersampling artifacts are removed using an iterative thresholding algorithm applied to nonlinearly transformed image block arrays. Each block array is transformed using kernel principal component analysis where the contribution of each image block to the transform depends in a nonlinear fashion on the distance to other image blocks. Elimination of undersampling artifacts is achieved by conventional principal component analysis in the nonlinear transform domain, projection onto the main components and back-mapping into the image domain. Iterative image reconstruction is performed by interleaving the proposed undersampling artifact removal step and gradient updates enforcing consistency with acquired k-space data. The algorithm is evaluated using retrospectively undersampled MR cardiac cine data and compared to k-t SPARSE-SENSE, block matching with spatial Fourier filtering and k-t ℓ₁-SPIRiT reconstruction. Evaluation of image quality and root-mean-squared-error (RMSE) reveal improved image reconstruction for up to 8-fold undersampled data with the proposed approach relative to k-t SPARSE-SENSE, block matching with spatial Fourier filtering and k-t ℓ₁-SPIRiT. In conclusion, block matching and kernel methods can be used for effective removal of undersampling artifacts in MR image reconstruction and outperform methods using standard compressed sensing and ℓ₁-regularized parallel imaging methods.     

Disease Progression in CNGA3 and CNGB3 Retinopathy; Characteristics of Slovenian Cohort and Proposed OCT Staging Based on Pooled Data from 126 Patients from 7 Studies

Achromatopsia has been proposed to be a morphologically predominately stable retinopathy with rare reports of progression of structural changes in the macula. A five-grade system of optical coherence tomography (OCT) features has been used for the classification of structural macular changes. However, their association with age remains questionable. We characterized the Slovenian cohort of 12 patients with pathogenic variants in CNGA3 or CNGB3 who had been followed up with OCT for up to 9 years. Based on observed structural changes in association with age, the following four-stage classification of retinal morphological changes was proposed: (I) preserved inner segment ellipsoid band (Ise), (II) disrupted ISe, (III) ISe loss and (IV) ISe and RPE loss. Data from six previously published studies reporting OCT morphology in CNGA3 and CNGB3 patients were additionally collected, forming the largest CNGA3/CNGB3 cohort to date, comprising 126 patients aged 1–71 years. Multiple regression analysis showed a significant correlation of OCT stage with age (p < 0.001) and no correlation with gene (p > 0.05). The median ages of patients with stages I–IV were 12 years, 23 years, 27 years and 48 years, respectively, and no patient older than 50 years had continuous ISe. Our findings suggest that achromatopsia presents with slowly but steadily progressive structural changes of the macular outer retinal layers. However, whether morphological changes in time follow the proposed four-stage linear pattern needs to be confirmed in a long-term study.     

Automated correlation of single particle tilt pairs for Random Conical Tilt and Orthogonal Tilt Reconstructions

One of the major methodological challenges in single particle electron microscopy is obtaining initial reconstructions which represent the structural heterogeneity of the dataset. Random Conical Tilt and Orthogonal Tilt Reconstruction techniques in combination with 3D alignment and classification can be used to obtain initial low-resolution reconstructions which represent the full range of structural heterogeneity of the dataset. In order to achieve statistical significance, however, a large number of 3D reconstructions, and, in turn, a large number of tilted image pairs are required. The extraction of single particle tilted image pairs from micrographs can be tedious and time-consuming, as it requires intensive user input even for semi-automated approaches. To overcome the bottleneck of manual selection of a large number of tilt pairs, we developed an algorithm for the correlation of single particle images from tilted image pairs in a fully automated and user-independent manner. The algorithm reliably correlates correct pairs even from noisy micrographs. We further demonstrate the applicability of the algorithm by using it to obtain initial references both from negative stain and unstained cryo datasets.     

Image Quality in Children with Low-Radiation Chest CT Using Adaptive Statistical Iterative Reconstruction and Model-Based Iterative Reconstruction

Objective

To evaluate noise reduction and image quality improvement in low-radiation dose chest CT images in children using adaptive statistical iterative reconstruction (ASIR) and a full model-based iterative reconstruction (MBIR) algorithm.

Methods

Forty-five children (age ranging from 28 days to 6 years, median of 1.8 years) who received low-dose chest CT scans were included. Age-dependent noise index (NI) was used for acquisition. Images were retrospectively reconstructed using three methods: MBIR, 60% of ASIR and 40% of conventional filtered back-projection (FBP), and FBP. The subjective quality of the images was independently evaluated by two radiologists. Objective noises in the left ventricle (LV), muscle, fat, descending aorta and lung field at the layer with the largest cross-section area of LV were measured, with the region of interest about one fourth to half of the area of descending aorta. Optimized signal-to-noise ratio (SNR) was calculated.

Result

In terms of subjective quality, MBIR images were significantly better than ASIR and FBP in image noise and visibility of tiny structures, but blurred edges were observed. In terms of objective noise, MBIR and ASIR reconstruction decreased the image noise by 55.2% and 31.8%, respectively, for LV compared with FBP. Similarly, MBIR and ASIR reconstruction increased the SNR by 124.0% and 46.2%, respectively, compared with FBP.

Conclusion

Compared with FBP and ASIR, overall image quality and noise reduction were significantly improved by MBIR. MBIR image could reconstruct eligible chest CT images in children with lower radiation dose.     

Semi-automated 3D Leaf Reconstruction and Analysis of Trichome Patterning from Light Microscopic Images

Trichomes are leaf hairs that are formed by single cells on the leaf surface. They are known to be involved in pathogen resistance. Their patterning is considered to emerge from a field of initially equivalent cells through the action of a gene regulatory network involving trichome fate promoting and inhibiting factors. For a quantitative analysis of single and double mutants or the phenotypic variation of patterns in different ecotypes, it is imperative to statistically evaluate the pattern reliably on a large number of leaves. Here we present a method that enables the analysis of trichome patterns at early developmental leaf stages and the automatic analysis of various spatial parameters. We focus on the most challenging young leaf stages that require the analysis in three dimensions, as the leaves are typically not flat. Our software TrichEratops reconstructs 3D surface models from 2D stacks of conventional light-microscope pictures. It allows the GUI-based annotation of different stages of trichome development, which can be analyzed with respect to their spatial distribution to capture trichome patterning events. We show that 3D modeling removes biases of simpler 2D models and that novel trichome patterning features increase the sensitivity for inter-accession comparisons.     

Adaptive Iterative Dose Reduction Using Three Dimensional Processing (AIDR3D) Improves Chest CT Image Quality and Reduces Radiation Exposure

Objective

To assess the advantages of Adaptive Iterative Dose Reduction using Three Dimensional Processing (AIDR3D) for image quality improvement and dose reduction for chest computed tomography (CT).

Methods

Institutional Review Boards approved this study and informed consent was obtained. Eighty-eight subjects underwent chest CT at five institutions using identical scanners and protocols. During a single visit, each subject was scanned using different tube currents: 240, 120, and 60 mA. Scan data were converted to images using AIDR3D and a conventional reconstruction mode (without AIDR3D). Using a 5-point scale from 1 (non-diagnostic) to 5 (excellent), three blinded observers independently evaluated image quality for three lung zones, four patterns of lung disease (nodule/mass, emphysema, bronchiolitis, and diffuse lung disease), and three mediastinal measurements (small structure visibility, streak artifacts, and shoulder artifacts). Differences in these scores were assessed by Scheffe''s test.

Results

At each tube current, scans using AIDR3D had higher scores than those without AIDR3D, which were significant for lung zones (p<0.0001) and all mediastinal measurements (p<0.01). For lung diseases, significant improvements with AIDR3D were frequently observed at 120 and 60 mA. Scans with AIDR3D at 120 mA had significantly higher scores than those without AIDR3D at 240 mA for lung zones and mediastinal streak artifacts (p<0.0001), and slightly higher or equal scores for all other measurements. Scans with AIDR3D at 60 mA were also judged superior or equivalent to those without AIDR3D at 120 mA.

Conclusion

For chest CT, AIDR3D provides better image quality and can reduce radiation exposure by 50%.     

共聚焦图像数据的三维重建和显示

本文介绍利用激光扫描共聚焦显微镜获得的共聚焦图像的三维重建和显示方法,并以ＡＣＡＳ Ｕｌｔｉｎａ３１２激光扫描共聚焦显微镜系统为例,分析了ＳＰＦ算法、投影算法和深度阴影算法等共聚焦图像数据的三维重建和图像显示方法的特点。     

浙江发现无尾蹄蝠及其3D头骨图像重建

2019年11月23日,在浙江省丽水市莲都区白云山森林公园1废弃矿洞内（28°29′51″N,119°54′34″E,海拔306 m）采集到蝙蝠标本1号。此蝙蝠主要特征为,前臂长35.4 mm;耳壳半透明,呈圆形漏斗;无尾,股间膜内凹呈“∧”形;背毛基部黑褐色,毛尖赤褐色;翼膜黑褐色。基于形态及分子数据（Cyt b）鉴定为无尾蹄蝠（Coelops frithii）,为浙江省兽类分布新记录种。本研究应用CT扫描和计算机辅助三维图像重建技术,在保证标本完整性的前提下,对其头骨进行CT扫描与三维（3D）图像复原,高清显示了其头部骨骼形态,该技术的应用,是在小型兽类骨骼形态学研究与描述上的有益探索和尝试。