Assessment of regional non-linear tissue deformation and air volume change of human lungs via image registration

We evaluate the non-linear characteristics of the human lung via image registration-derived local variables based on volumetric multi-detector-row computed tomographic (MDCT) lung image data of six normal human subjects acquired at three inflation levels: 20% of vital capacity (VC), 60% VC and 80% VC. Local variables include Jacobian (ratio of volume change) and maximum shear strain for assessment of lung deformation, and air volume change for assessment of air distribution. First, the variables linearly interpolated between 20% and 80% VC images to reflect deformation from 20% to 60% VC are compared with those of direct registration of 20% and 60% VC images. The result shows that the linearly-interpolated variables agree only qualitatively with those of registration (P<0.05). Then, a quadratic (or linear) interpolation is introduced to link local variables to global air volumes of three images (or 20% and 80% VC images). A sinusoidal breathing waveform is assumed for assessing the time rate of change of these variables. The results show significant differences between two-image and three-image results (P<0.05). The three-image results for the whole lung indicate that the peak of the maximum shear rate occurs at about 37% of the maximum volume difference between 20% and 80% VC, while the peaks for the Jacobian and flow rate occur at 50%. This is in agreement with accepted physiology whereby lung tissues deform more at lower lung volumes due to lower elasticity and greater compliance. Furthermore, the three-image results show that the upper and middle lobes, even in the recumbent, supine posture, reach full expansion earlier than the lower lobes.     

Respiratory motion of adrenal gland metastases: Analyses using four-dimensional computed tomography images

PurposeTo evaluate the respiratory motion of adrenal gland metastases in three-dimensional directions using four-dimensional computed tomography (4DCT) images.MethodsFrom January 2013 to May 2016, 12 patients with adrenal gland metastases were included in this study. They all underwent 4DCT scans to assess respiratory motion of adrenal gland metastases in free breathing state. The 4DCT images were sorted into 10 image series according to the respiratory phase from the end inspiration to the end expiration, and then transferred to FocalSim workstation. All gross tumor volumes (GTVs) of adrenal gland metastases were drawn by a single physician and confirmed by a second. Relative coordinates of adrenal gland metastases were automatically generated to calculate adrenal gland metastases motion in different axial directions.ResultsThe average respiratory motion of adrenal gland metastases in left-right (LR), cranial-caudal (CC), anterior-posterior (AP), 3-dimensional (3D) vector directions was 3.4 ± 2.2 mm, 9.5 ± 5.5 mm, 3.8 ± 2.0 mm and 11.3 ± 5.3 mm, respectively. The ratios were 58.6% ± 11.4% and 63.2% ± 12.5% when the volumes of GTV_In0% and GTV _In100% were compared with volume of IGTV_10phase. The volume ratio of IGTV_10phase to GTV_3D was 1.73 ± 0.48.ConclusionsAdrenal gland metastasis is a respiration-induced moving target, and an internal target volume boundary should be provided when designing the treatment plan. The CC motion of adrenal gland metastasis is predominant and >5 mm, thus motion management strategies are recommended for patients undergoing external radiotherapy for adrenal gland metastasis.     

Registration methods in radiotherapy

PurposeThe aim of this study is to present a short and comprehensive review of the methods of medical image registration, their conditions and applications in radiotherapy. A particular focus was placed on the methods of deformable image registration.MethodsTo structure and deepen the knowledge on medical image registration in radiotherapy, a medical literature analysis was made using the Google Scholar browser and the medical database of the PubMed library.ResultsChronological review of image registration methods in radiotherapy based on 34 selected articles. A particular attention was given to show: (i) potential regions of the application of different methods of registration, (ii) mathematical basis of the deformable methods and (iii) the methods of quality control for the registration process.ConclusionsThe primary aim of the medical image registration process is to connect the contents of images. What we want to achieve is a complementary or extended knowledge that can be used for more precise localisation of pathogenic lesions and continuous improvement of patient treatment. Therefore, the choice of imaging mode is dependent on the type of clinical study. It is impossible to visualise all anatomical details or functional changes using a single modality machine. Therefore, fusion of various modality images is of great clinical relevance. A natural problem in analysing the fusion of medical images is geographical errors related to displacement. The registered images are performed not at the same time and, very often, at different respiratory phases.     

Lung clearance index in adults with non-cystic fibrosis bronchiectasis

Background

Lung clearance index (LCI) is a measure of abnormal ventilation distribution derived from the multiple breath inert gas washout (MBW) technique. We aimed to determine the clinical utility of LCI in non-CF bronchiectasis, and to assess two novel MBW parameters that distinguish between increases in LCI due to specific ventilation inequality (LCI_vent) and increased respiratory dead space (LCI_ds).

Methods

Forty-three patients with non-CF bronchiectasis and 18 healthy control subjects underwent MBW using the sulphur hexafluoride wash-in technique, and data from 40 adults with CF were re-analysed. LCI_vent and LCI_ds were calculated using a theoretical two-compartment lung model, and represent the proportional increase in LCI above its ideal value due to specific ventilation inequality and increased respiratory dead space, respectively.

Results

LCI was significantly raised in patients with non-CF bronchiectasis compared to healthy controls (9.99 versus 7.28, p < 0.01), and discriminated well between these two groups (area under receiver operating curve = 0.90, versus 0.83 for forced expiratory volume in one second [% predicted]). LCI, LCI_vent and LCI_ds were repeatable (intraclass correlation coefficient > 0.75), and correlated significantly with measures of spirometric airflow obstruction.

Conclusion

LCI is repeatable, discriminatory, and is associated with spirometric airflow obstruction in patients with non-CF bronchiectasis. LCI_vent and LCI_ds are a practical and repeatable alternative to phase III slope analysis and may allow a further level of mechanistic information to be extracted from the MBW test in patients with severe ventilation heterogeneity.     

Effect of diagnostic cone-beam computed tomography protocols on image quality,patient dose,and lesion detection

ObjectiveTo evaluate the effect of cone-beam computed tomography (CBCT) image acquisition protocols on image quality, lesion detection, delineation, and patient dose.Methods100-patients and a CTDI phantom combined with an electron density phantom were examined using four different CBCT-image acquisition protocols during image-guided transarterial chemoembolization (TACE). Protocol-1 (time: 6 s, tube rotation: 360°), protocol-2 (5 s, 300°), protocol-3 (4 s, 240°) and protocol-4 (3 s, 180°) were used. The protocols were first investigated using a phantom. The protocols that were found to be clinically appropriate in terms of image quality and radiation dose were then assessed on patients. A higher radiation dose and/or a poor image quality were inappropriate for the patient imaging. Patient dose (patient-entrance dose and dose-area product), image quality (Hounsfield Unit, noise, signal-to-noise ratio and contrast-to-noise ratio), and lesion delineation (tumor-liver contrast) were assessed and compared using appropriate statistical tests. Lesion detectability, sensitivity, and predictive values were estimated for CBCT-image data using pre-treatment patient magnetic resonance imaging.ResultsThe estimated patient dose showed no statistical significance (p > 0.05) between protocols-2 and -3; the assessed image quality between these protocols manifested insignificant difference (p > 0.05). Two other phantom protocols were not considered for patient imaging due to significantly higher dose (protocols-1) and poor image quality (protocol-4). Lesion delineation and detection were insignificant (p > 0.05) between protocols-2 and -3. Lesion sensitivities generated were 81–89% (protocol-2) and 81–85% (protocol-3) for different lesion types.ConclusionData acquisition using protocols-2 and -3 provided good image quality, lesion detection and delineation with acceptable patient dose during CBCT-imaging mainly due to similar frame numbers acquired.     

智能期相联合SSF技术改善冠状动脉CTA图像运动伪影

目的:评价智能期相选择(Smart phase)和自动追踪冻结(Snap-Shot-Freeze, SSF)技术联合应用在前门控冠状动脉CT血管成像(Coronary CT angiography, CCTA)改善图像运动伪影中的价值。方法:94例疑似心血管疾病的患者行CCTA检查。所有患者图像均采用smart phase、SSF及smart phase+SSF两种技术联合进行后处理。将原始图像(A组)、smart phase图像(B组)、SSF图像(C组)和smart phase+SSF图像(D组)导入GE AW4.6工作站。对四组患者冠状动脉三大分支(右冠状动脉、左前降支及左旋支)图像质量进行主观评价(5分法)和测量每组图像主动脉根部及前胸壁肌肉CT值和标准差,并计算图像的噪声(Noise,N)、信噪比(Signal-to-noise ratio,SNR)和对比噪声比(Contrast-to-noise ratio,CNR)和比较。结果:经smart phase处理后,A组冠状动脉三大分支图像评分均得到了明显提高(P<0.05)。经SSF矫正的A组和B组冠状动脉三大分支图像质量评分均得到明显提高(P<0.05)。Smart Phase联合SSF技术得到的图像评分在所有组中评分最高,与其他组均存在显著差异(P<0.05)。结论:Smart phase联合SSF技术可明显改善前门控CCTA图像运动伪影,值得临床推广应用。     

基于恒定空间分辨率离散序列小波变换的眼底血管造影图像拼接方法

将恒定空间分辨率离散序列小波变换（discrete sequence wavelet transform,DSWT)应用于眼底吲哚青绿血管造影（indocyanine green angiography,ICGA)图像的拼接,解决了传统基为2的DSWT会导致分解结果的空间分辨率下降的问题。提出对图像小波分解细节逼近和平滑逼近分别使用加权平均拼接和直接平均拼接进行处理的策略,以得到兼顾视觉效果和保真性的拼接结果。并且针对眼底图像背景光照不一致,提出在小波域进行处理的策略。实验结果表明拼接算法效果良好。     

Impact of inflammation,emphysema, and smoking cessation on V/Q in mouse models of lung obstruction

Background

Chronic obstructive pulmonary disease (COPD) is known to greatly affect ventilation (V) and perfusion (Q) of the lung through pathologies such as inflammation and emphysema. However, there is little direct evidence regarding how these pathologies contribute to the V/Q mismatch observed in COPD and models thereof. Also, little is known regarding how smoking cessation affects V/Q relationships after inflammation and airspace enlargement have become established. To this end, we have quantified V/Q on a per-voxel basis using single photon emission computed tomography (SPECT) in mouse models of COPD and lung obstruction.

Methods

Three distinct murine models were used to investigate the impact of different pathologies on V/Q, as measured by SPECT. Lipopolysaccharide (LPS) was used to produce neutrophilic inflammation, porcine pancreatic elastase (PPE) was used to produce emphysema, and long-term cigarette smoke (CS) exposure and cessation were used to investigate the combination of these pathologies.

Results

CS exposure resulted in an increase in mononuclear cells and neutrophils, an increase in airspace enlargement, and an increase in V/Q mismatching. The inflammation produced by LPS was more robust and predominantly neutrophilic, compared to that of cigarette smoke; nevertheless, inflammation alone caused V/Q mismatching similar to that seen with long-term CS exposure. The emphysematous lesions caused by PPE administration were also capable of causing V/Q mismatch in the absence of inflammation. Following CS cessation, inflammatory cell levels returned to those of controls and, similarly, V/Q measures returned to normal despite evidence of persistent mild airspace enlargement.

Conclusions

Both robust inflammation and extensive airspace enlargement, on their own, were capable of producing V/Q mismatch. As CS cessation resulted in a return of V/Q mismatching and inflammatory cell counts to control levels, lung inflammation is likely a major contributor to V/Q mismatch observed in the cigarette smoke exposure model as well as in COPD patients. This return of V/Q mismatching to control values also took place in the presence of mild airspace enlargement, indicating that emphysematous lesions must be of a larger volume before affecting the lung significantly. Early smoking cessation is therefore critical before emphysema has an irreversible impact on gas exchange.     

Influence of the type of imaging on the delineation process during the treatment planning

AimThe aim of this study was to compare the intra- and interobserver contouring variability for structures with density of organ at risk in two types of tomography: kilovoltage computed tomography (KVCT) versus megavoltage computed tomography (MVCT). The intra- and interobserver differences were examined on both types of tomography for structures which simulate human tissue or organs.Materials and methodsSix structures with density of the liver, bone, trachea, lung, soft tissue and muscle were created and used. For the measurements, the special water phantom with all structures was designed. To evaluate interobserver variability, five observers delineated the structures in both types of computed tomography (CT).ResultsIntraobserver variability was in the range of 1–14% and was the largest for the liver. The observers segmented larger volumes on MVCT compared with KVCT for the trachea (79.56 ccm vs.74.91 ccm), lung (87.61 vs. 82.50), soft tissue (154.24 vs. 145.47) and muscle (164.01 vs. 157.89). For the liver (98.13 vs. 99.38) and bone (51.86 vs. 67.97), the volume on MVCT was smaller than KVCT. The statistically significant differences between observers were observed for structures with density of the liver, bone and soft tissue on KVCT and for the liver, lung and soft tissue on MVCT. For the structures with density of the trachea and muscles, there were no significant differences for both types of tomography.ConclusionsDuring the contouring process the interobserver and intraobserver contouring uncertainty was larger on MVCT, especially for structures with HU near 80, compared with KVCT.     

基于配准的插值

断层间图像插值是三维重建的一个关键步骤,因为图像上像素之间的间隔常常小于断层图像之间的距离,而在三维重建需要它们有一致的分辨率.由于是同模态断层图像层间插值,对于解决同模态弹性配准问题,Thirion的demons算法比较适合.所以配准采用Demons方法.Demons算法先判断出待配准图像上各个象素点的运动方法,通过对各个象素点的移动来实现非刚性配准.在这个算法中,每张图像都被视为同灰度值轮廓的集合.该算法可以应用于精度要求比较高的体数据插值重建过程.