Evaluation of deformable image registration for improved 4D CT-derived ventilation for image guided radiotherapy

Recent treatment planning studies have demonstrated the use of physiologic images in radiation therapy treatment planning to identify regions for functional avoidance. This image-guided radiotherapy (IGRT) strategy may reduce the injury and/or functional loss following thoracic radiotherapy. 4D computed tomography (CT), developed for radiotherapy treatment planning, is a relatively new imaging technique that allows the acquisition of a time-varying sequence of 3D CT images of the patient\u27s lungs through the respiratory cycle. Guerrero et al. developed a method to calculate ventilation imaging from 4D CT, which is potentially better suited and more broadly available for IGRT than the current standard imaging methods. The key to extracting function information from 4D CT is the construction of a volumetric deformation field that accurately tracks the motion of the patient\u27s lungs during the respiratory cycle. The spatial accuracy of the displacement field directly impacts the ventilation images; higher spatial registration accuracy will result in less ventilation image artifacts and physiologic inaccuracies. Presently, a consistent methodology for spatial accuracy evaluation of the DIR transformation is lacking. Evaluation of the 4D CT-derived ventilation images will be performed to assess correlation with global measurements of lung ventilation, as well as regional correlation of the distribution of ventilation with the current clinical standard SPECT. This requires a novel framework for both the detailed assessment of an image registration algorithm\u27s performance characteristics as well as quality assurance for spatial accuracy assessment in routine application. Finally, we hypothesize that hypo-ventilated regions, identified on 4D CT ventilation images, will correlate with hypo-perfused regions in lung cancer patients who have obstructive lesions. A prospective imaging trial of patients with locally advanced non-small-cell lung cancer will allow this hypothesis to be tested. These advances are intended to contribute to the validation and clinical implementation of CT-based ventilation imaging in prospective clinical trials, in which the impact of this imaging method on patient outcomes may be tested

In-house Implementation and Validation of the Mid-Position CT approach for the Treatment Planning of Respiration-induced Moving Tumours in Radiotherapy for Lung and Upper abdomen cancer

Tese mestrado integrado, Engenharia Biomédica e Biofísica (Engenharia Clínica e Instrumentação Médica) Universidade de Lisboa, Faculdade de Ciências, 2022A Radioterapia é uma das modalidades principais para tratamentos de foro oncológico que visa destruir a ação proliferativa das células cancerígenas e reduzir o volume tumoral. A sua ação terapêutica através do uso de radiação ionizante tem, subjacente, a máxima de irradiar o tumor com uma elevada dose, ao mesmo tempo que os órgãos de risco (OARs) adjacentes, são tanto quanto possível protegidos. Quando um tumor se localiza no pulmão ou abdómen superior, como no fígado ou pâncreas, o seu movimento devido à respiração pode alcançar até 4 cm, especialmente na direção crânio-caudal, aumentando as incertezas relativas à posição do tumor. No Centro Clínico Champalimaud (CCC), o planeamento convencional dos tratamentos de radioterapia faz uso de uma tomografia computadorizada (CT) que é adquirida aquando da respiração livre do doente e que, por isso, apresenta geralmente artefactos que podem ser uma fonte de erro durante o planeamento. Nos casos em que o movimento do tumor é considerável, é ainda adquirida uma tomografia computadorizada quadrimensional (4DCT) que consiste entre 8 e 10 CTs que representam fases do ciclo respiratório. Posteriormente, a 4DCT é utilizada para delinear o volume interno do alvo (ITV) que engloba toda a extensão do movimento do tumor. Apesar da estratégia do ITV garantir uma adequada cobertura do volume-alvo, os OARs ficam expostos a doses de radiação desnecessárias e a um maior risco de toxicidade. Este efeito é ainda mais preocupante em tratamentos hipofracionados, onde doses mais elevadas são administradas num número reduzido de frações. Nos últimos anos têm sido desenvolvidas estratégias que visam tornar os tratamentos de radioterapia mais eficazes. Uma delas é a reconstrução de uma CT que representa a posição média do doente ao longo do ciclo respiratório (Mid-P CT). Esta estratégia resulta em volumes de tratamento menores do que a estratégia do ITV, possibilitando o aumento da dose e maior controlo tumoral local. O primeiro passo para a reconstrução do Mid-P CT é o registo deformável de imagens (DIR) entre uma das fases da respiração (uma CT da 4DCT), definida como a fase de referência, e as restantes fases. Deste processo resultam campos vetoriais deformáveis (DVF) que contém informação do deslocamento dos tecidos. Os DVFs são subsequentemente utilizados para transformar cada uma das fases da respiração para a posição média. O método do Mid-P foi implementado com sucesso no Instituto do Cancro Holandês (NKI) em 2008. Apesar dos bons resultados clínicos, o número de centros de radioterapia que utiliza esta técnica é muito reduzido. Tal deve-se, por um lado, à inexistência de soluções comerciais com esta funcionalidade e, por outro, ao esforço necessário alocar para implementar e validar soluções desenvolvidas internamente. O presente projeto teve como principal objetivo implementar a estratégia do Mid-P no CCC (Portugal). Para tal, foi otimizado um módulo – RunMidP – desenvolvido para o software 3D Slicer, que calcula o Mid-P CT e estima a amplitude do movimento do tumor e OARs com base nos DVFs. Considerando que a precisão do módulo e a qualidade de imagem do Mid-P CT devem atender os requisitos para o planeamento em radioterapia, foram realizados testes para validar o módulo. Sempre que possível, a sua performance foi comparada com outras aplicações desenvolvidas para a implementação da técnica do Mid-P, nomeadamente com um protótipo desenvolvido pela empresa Mirada Medical Ltd. (Reino Unido) – Mirada – e com o software desenvolvido no NKI (Holanda) – Wimp. Os testes foram divididos em três estudos diferentes, cada um com um conjunto de dados diferente. No primeiro estudo (estudo A), foram utilizadas 4DCT de 2 fantomas digitais, cuja função respiratória e cardíaca foi modelada de forma simplificada, e de 18 doentes com tumores localizados no pulmão (N = 8), no fígado (N = 6) e no pâncreas (N = 4). Neste estudo, foram comparados dois algoritmos DIR disponíveis no software 3D Slicer, o Plastimatch e o Elastix, em termos da precisão do registo e da qualidade de imagem do Mid-P CT reconstruído. Foi ainda avaliado a capacidade dos softwares RunMidP e Mirada representarem corretamente a posição média do doente e as diferenças das amplitudes do movimento do tumor estimadas pelos dois softwares. No estudo B, foram realizados testes de verificação semelhantes aos supre mencionados, em imagens sintéticas provenientes de 16 doentes, desta vez com a vantagem de se conhecer o “verdadeiro” Mid-P CT e as “verdadeiras” amplitudes do movimento do tumor. Estes foram comparados com os resultados obtidos com os softwares RunMidP e Mirada. Ainda, as unidades de Hounsfield (HU) no Mid-P CT reconstruído por RunMidP e Mirada foram comparadas com as HU na fase de referência, de modo a verificar se os Mid P CTs produziriam diferenças dosimétricas relevantes. No último estudo (estudo C), a qualidade de imagem do Mid-P CT foi avaliada quantitativamente e qualitativamente. Durante a análise qualitativa, foi pedido a dois médicos especialistas que avaliassem a viabilidade dos Mid-P CTs, reconstruídos pelos três softwares (RunMidP, Mirada e Wimp), para o planeamento dos tratamentos. O tempo da reconstrução do Mid-P CT a partir da 4DCT foi de cerca de 1h. Ambos os algoritmos, Plastimach e Elastix, demonstraram ser adequados para DIR de imagens do pulmão e abdómen superior, com diferenças estatisticamente não significativas (p > 0.05) em termos da precisão do registo. Contudo, o Mid-P CT reconstruído com Elastix apresentou uma melhoria na qualidade de imagem, sendo assim o algoritmo DIR escolhido para ser implementado no RunMidP. Em termos de métricas aplicadas a contornos definidos manualmente, tais como a distância de Hausdorf (HD) e coeficiente de Dice (DSC), o erro do registo de imagem foi menor que 1 mm, dentro do contorno do tumor, e 2 mm no pulmão. Os Mid-P CTs reconstruídos com o RunMidP e Mirada apresentaram maiores diferenças, relativamente ao “verdadeiro” Mid-P CT, na região do diafragma e zonas de maior homogeneidade como, por exemplo, no ar presente no intestino. Contudo, para a maioria dos doentes do estudo B, o Mid-P CT reconstruído com o software Mirada apresentou maior índice de similaridade estrutural (SSIM) relativamente ao “verdadeiro” Mid-P CT. Estes resultados podem estar na origem do uso de diferentes algoritmos DIR, mas deveram-se principalmente a uma falha na aplicação das transformações deformáveis pelo módulo RunMiP que foi corrigida posteriormente. Ainda, as diferenças entre as amplitudes estimadas e previstas foram menores que 1 mm para 37 tumores (78,9%), que resultam em diferenças menores que 0.3mm quando convertidas em margens de planeamento. Para além disso, as diferenças nos valores de HU dos Mid-P CTs comparativamente à fase de referência foram, em média, de 1 HU no tumor e OARs. Foram também observadas melhorias na qualidade de imagem do Mid-P CT, nomeadamente um aumento da relação sinal-ruído (SNR) e diminuição dos artefactos. Estes resultados estão de acordo com a avaliação dos médicos que, em geral, consideraram que os Mid-P CTs reconstruídos pelos três softwares são adequados para o planeamento dos tratamentos. No entanto, os Mid-P CTs reconstruídos com dados 4DCT provenientes do CCC apresentaram classificações inferiores aos reconstruídos com dados 4DCT do NKI. Em suma, as modificações do algoritmo DIR Plastimach para Elastix e a correção do método para aplicar as transformações deformáveis, permitiram uma melhoria na qualidade de imagem do Mid P CT e melhor performance do algoritmo, respetivamente. O módulo RunMidP, neste projeto otimizado e validado, apresenta um forte potencial para a reconstrução e implementação da estratégia do Mid-P na clínica, com performance comparável a outras aplicações existentes (Mirada e Wimp). Atenção especial deve ser dada aos dados 4DCT de input que parecem afetar a qualidade de imagem final do Mid-P CT. No futuro, valerá a pena otimizar os parâmetros de aquisição e reconstrução da 4DCT de modo a melhorar a qualidade de imagem e, ainda, o módulo RunMidP pode potencialmente ser otimizado no que respeita ao tempo de reconstrução do Mid-P CT e à precisão do DIR.Radiotherapy for tumours in the thorax and upper abdomen is challenging since they move notably with breathing. To cover the whole extent of tumour motion, relatively large margins are added to treatment volumes, posing a higher risk of toxicity for surrounding organs-at-risk (OARs). The Mid Position (Mid-P) method accounts for breathing motion by using deformable image registration (DIR) to transform all phases of a 4DCT scan to a time-weighted average 3DCT scan (Mid-P CT). The Mid-P strategy results in smaller treatment volumes, potentially boosting the delivery of hypofractionated treatments. To bring the Mid-P approach to the Champalimaud Clinical Centre (CCC), an in-house Mid position software module – RunMidP – was optimized. The module reconstructs the Mid-P CT and estimates breathing motion amplitudes of tumours and relevant OARs. In addition, this project presents a set of experiments to evaluate the performance of the Mid-P method and its feasibility for clinical implementation. The experiments were conducted throughout three different studies using 4DCT data from 18 phantoms and 23 patients. In Study A, the accuracy and image quality of two DIR algorithms (Plastimatch and Elastix) were assessed using quantitative metrics applied on either warped images or manually delineated contours. The reproduction of the patient’s mean position by the Mid-P CT and the estimation of motion amplitudes were compared to a soon-to-be Mid-P commercial software developed by Mirada Medical Ltd. In Study B,similar experiments were performed, this time using a more rigorous reference – “true” Mid-P CT scans and “true” motion estimations. In Study C, the image quality of Mid P CT scans was assessed quantitatively and qualitatively. Both Plastimatch and Elastix registration showed comparable registration accuracy, although Elastix showed superior image quality of reconstructed Mid-P CTs. Based on contour metrics, the registration error was less than 2 mm. In-house Mid-P CTs showed a slightly lower match to ground truth Mid-P CTs than the ones reconstructed by the Mirada prototype due to differences in DIR methods and small shifts to the original image geometry. Higher image differences were found in the diaphragm lung interface, where the patient's anatomy moves faster due to breathing, and in homogeneous regions such as the air regions in the bowel. On the other hand, differences (estimated-predicted) in motion amplitudes smaller than 1 mm were observed in 37 moving tumours (78.7%), showing a good performance of the Mid-P algorithm. Regarding the image quality, improvements in the signal-to-noise ratio and removal of image artefacts in Mid-P CTs are great advantages for using them as the planning CT. Clinicians also gave a good assessment of the suitability of Mid-P CT scans for treatment planning. No significant differences were found in the performance of the RunMidP compared to other Mid-Position packages, although worse scores were given to the CCC dataset than the dataset from another hospital. The in-house Mid-position algorithm shows promising results regarding the use of the software module in radiotherapy for lung and upper abdomen cancer. Further exploration must be given to improve the registration accuracy, image quality of the input data, and speed up the reconstruction of the Mid-P CT scan

3-D lung deformation and function from respiratory-gated 4-D x-ray CT images : application to radiation treatment planning.

Many lung diseases or injuries can cause biomechanical or material property changes that can alter lung function. While the mechanical changes associated with the change of the material properties originate at a regional level, they remain largely asymptomatic and are invisible to global measures of lung function until they have advanced significantly and have aggregated. In the realm of external beam radiation therapy of patients suffering from lung cancer, determination of patterns of pre- and post-treatment motion, and measures of regional and global lung elasticity and function are clinically relevant. In this dissertation, we demonstrate that 4-D CT derived ventilation images, including mechanical strain, provide an accurate and physiologically relevant assessment of regional pulmonary function which may be incorporated into the treatment planning process. Our contributions are as follows: (i) A new volumetric deformable image registration technique based on 3-D optical flow (MOFID) has been designed and implemented which permits the possibility of enforcing physical constraints on the numerical solutions for computing motion field from respiratory-gated 4-D CT thoracic images. The proposed optical flow framework is an accurate motion model for the thoracic CT registration problem. (ii) A large displacement landmark-base elastic registration method has been devised for thoracic CT volumetric image sets containing large deformations or changes, as encountered for example in registration of pre-treatment and post-treatment images or multi-modality registration. (iii) Based on deformation maps from MOFIO, a novel framework for regional quantification of mechanical strain as an index of lung functionality has been formulated for measurement of regional pulmonary function. (iv) In a cohort consisting of seven patients with non-small cell lung cancer, validation of physiologic accuracy of the 4-0 CT derived quantitative images including Jacobian metric of ventilation, Vjac, and principal strains, (V?1, V?2, V?3, has been performed through correlation of the derived measures with SPECT ventilation and perfusion scans. The statistical correlations with SPECT have shown that the maximum principal strain pulmonary function map derived from MOFIO, outperforms all previously established ventilation metrics from 40-CT. It is hypothesized that use of CT -derived ventilation images in the treatment planning process will help predict and prevent pulmonary toxicity due to radiation treatment. It is also hypothesized that measures of regional and global lung elasticity and function obtained during the course of treatment may be used to adapt radiation treatment. Having objective methods with which to assess pre-treatment global and regional lung function and biomechanical properties, the radiation treatment dose can potentially be escalated to improve tumor response and local control

Comparison of hyperpolarised gas MRI and CT-based surrogates of ventilation

Background: Non-contrast CT-based surrogates of regional ventilation derived from pulmonary images acquired at multiple inflation levels have been proposed as alternatives to established modalities. However, their physiological accuracy has yet to be validated prior to clinical translation. Purpose: To address the hypothesis that these surrogates can provide information comparable to a direct measure of ventilation from hyperpolarised gas MRI ventilation via: i. development of a methodology for registering CT and gas MRI. ii. comparison of these surrogates with gas MRI at the lobar level. iii. evaluation of the impact of inflation levels when comparing gas MRI and ventilation CT. iv. development of an image acquisition and analysis framework to facilitate spatial correlations of both techniques. v. assessment of the effect of using different gases on the correlation. Methods: i. A method to indirectly register gas MRI to CT via same-breath 1H-structural MR images was developed and its accuracy was assessed. ii. A ventilation model based on expansion of lobar CT segmentations was compared with gas MRI lobar ventilation measurements. iii. The spatial overlap of ventilation CT was compared to gas MRI acquired at two different inflation levels. iv. An image acquisition protocol was designed to minimise differences in acquisition settings between scans such as posture and breathing manoeuvre and analysis methods were developed to enable direct regional and voxel level correlations. v. The effect of using two different noble gases, namely, 3He and 129Xe, on correlation with ventilation CT was assessed. Results: i. The indirect method of registration was more accurate than direct registration. ii. Despite subtle differences, lobar ventilation measurements derived from CT and hyperpolarised gas MRI were comparable. iii. Comparison of ventilation CT and gas MRI varied with inflation state. iv. The spatial correlation between ventilation CT and gas MRI increased at coarser levels. v. A marked improvement in correlation was observed for 3He and 129Xe MRI in contrast to when ventilation CT was compared with either 3He and 129Xe MRI. Conclusion: Although CT-based surrogates of ventilation show promise for replacing established ventilation modalities such as hyperpolarised gas MRI, particularly at coarser levels, they cannot be assumed to be equivalent to the techniques they purport to replace

Computational methods for the analysis of functional 4D-CT chest images.

Medical imaging is an important emerging technology that has been intensively used in the last few decades for disease diagnosis and monitoring as well as for the assessment of treatment effectiveness. Medical images provide a very large amount of valuable information that is too huge to be exploited by radiologists and physicians. Therefore, the design of computer-aided diagnostic (CAD) system, which can be used as an assistive tool for the medical community, is of a great importance. This dissertation deals with the development of a complete CAD system for lung cancer patients, which remains the leading cause of cancer-related death in the USA. In 2014, there were approximately 224,210 new cases of lung cancer and 159,260 related deaths. The process begins with the detection of lung cancer which is detected through the diagnosis of lung nodules (a manifestation of lung cancer). These nodules are approximately spherical regions of primarily high density tissue that are visible in computed tomography (CT) images of the lung. The treatment of these lung cancer nodules is complex, nearly 70% of lung cancer patients require radiation therapy as part of their treatment. Radiation-induced lung injury is a limiting toxicity that may decrease cure rates and increase morbidity and mortality treatment. By finding ways to accurately detect, at early stage, and hence prevent lung injury, it will have significant positive consequences for lung cancer patients. The ultimate goal of this dissertation is to develop a clinically usable CAD system that can improve the sensitivity and specificity of early detection of radiation-induced lung injury based on the hypotheses that radiated lung tissues may get affected and suffer decrease of their functionality as a side effect of radiation therapy treatment. These hypotheses have been validated by demonstrating that automatic segmentation of the lung regions and registration of consecutive respiratory phases to estimate their elasticity, ventilation, and texture features to provide discriminatory descriptors that can be used for early detection of radiation-induced lung injury. The proposed methodologies will lead to novel indexes for distinguishing normal/healthy and injured lung tissues in clinical decision-making. To achieve this goal, a CAD system for accurate detection of radiation-induced lung injury that requires three basic components has been developed. These components are the lung fields segmentation, lung registration, and features extraction and tissue classification. This dissertation starts with an exploration of the available medical imaging modalities to present the importance of medical imaging in today’s clinical applications. Secondly, the methodologies, challenges, and limitations of recent CAD systems for lung cancer detection are covered. This is followed by introducing an accurate segmentation methodology of the lung parenchyma with the focus of pathological lungs to extract the volume of interest (VOI) to be analyzed for potential existence of lung injuries stemmed from the radiation therapy. After the segmentation of the VOI, a lung registration framework is introduced to perform a crucial and important step that ensures the co-alignment of the intra-patient scans. This step eliminates the effects of orientation differences, motion, breathing, heart beats, and differences in scanning parameters to be able to accurately extract the functionality features for the lung fields. The developed registration framework also helps in the evaluation and gated control of the radiotherapy through the motion estimation analysis before and after the therapy dose. Finally, the radiation-induced lung injury is introduced, which combines the previous two medical image processing and analysis steps with the features estimation and classification step. This framework estimates and combines both texture and functional features. The texture features are modeled using the novel 7th-order Markov Gibbs random field (MGRF) model that has the ability to accurately models the texture of healthy and injured lung tissues through simultaneously accounting for both vertical and horizontal relative dependencies between voxel-wise signals. While the functionality features calculations are based on the calculated deformation fields, obtained from the 4D-CT lung registration, that maps lung voxels between successive CT scans in the respiratory cycle. These functionality features describe the ventilation, the air flow rate, of the lung tissues using the Jacobian of the deformation field and the tissues’ elasticity using the strain components calculated from the gradient of the deformation field. Finally, these features are combined in the classification model to detect the injured parts of the lung at an early stage and enables an earlier intervention

A Heterogeneous Patient-Specific Biomechanical Model of the Lung for Tumor Motion Compensation and Effective Lung Radiation Therapy Planning

Radiation therapy is a main component of treatment for many lung cancer patients. However, the respiratory motion can cause inaccuracies in radiation delivery that can lead to treatment complications. In addition, the radiation-induced damage to healthy tissue limits the effectiveness of radiation treatment. Motion management methods have been developed to increase the accuracy of radiation delivery, and functional avoidance treatment planning has emerged to help reduce the chances of radiation-induced toxicity. In this work, we have developed biomechanical model-based techniques for tumor motion estimation, as well as lung functional imaging. The proposed biomechanical model accurately estimates lung and tumor motion/deformation by mimicking the physiology of respiration, while accounting for heterogeneous changes in the lung mechanics caused by COPD, a common lung cancer comorbidity. A biomechanics-based image registration algorithm is developed and is combined with an air segmentation algorithm to develop a 4DCT-based ventilation imaging technique, with potential applications in functional avoidance therapies

Retrospective analysis of the impact of respiratory motion in treatment margins for frameless lung SBRT based on respiratory-correlated CBCT data-sets.

To investigate the impact of respiratory motion in the treatment margins for lung SBRT frameless treatments and to validate our treatment margins using 4D CBCT data analysis. Two hundred and twenty nine fractions with early stage NSCLC were retrospectively analyzed. All patients were treated in frameless and free breathing conditions. The treatment margins were calculated according to van Herk equation in Mid-Ventilation. For each fraction, three 4D CBCT scans, pre- and postcorrection, and posttreatment, were acquired to assess target baseline shift, target localization accuracy and intra-fraction motion errors. A bootstrap analysis was performed to assess the minimum number of patients required to define treatment margins. The retrospectively calculated target-baseline shift, target localization accuracy and intra-fraction motion errors agreed with the literature. The best tailored margins to our cohort of patients were retrospectively computed and resulted in agreement with already published data. The bootstrap analysis showed that fifteen patients were enough to assess treatment margins. The treatment margins applied to our patient's cohort resulted in good agreement with the retrospectively calculated margins based on 4D CBCT data. Moreover, the bootstrap analysis revealed to be a promising method to verify the reliability of the applied treatment margins for safe lung SBRT delivery

Biomechanical Modeling for Lung Tumor Motion Prediction during Brachytherapy and Radiotherapy

A novel technique is proposed to develop a biomechanical model for estimating lung’s tumor position as a function of respiration cycle time. Continuous tumor motion is a major challenge in lung cancer treatment techniques where the tumor needs to be targeted; e.g. in external beam radiotherapy and brachytherapy. If not accounted for, this motion leads to areas of radiation over and/or under dosage for normal tissue and tumors. In this thesis, biomechanical models were developed for lung tumor motion predication in two distinct cases of lung brachytherapy and lung external beam radiotherapy. The lung and other relevant surrounding organs geometries, loading, boundary conditions and mechanical properties were considered and incorporated properly for each case. While using material model with constant incompressibility is sufficient to model the lung tissue in the brachytherapy case, in external beam radiation therapy the tissue incompressibility varies significantly due to normal breathing. One of the main issues tackled in this research is characterizing lung tissue incompressibility variations and measuring its corresponding parameters as a function of respiration cycle time. Results obtained from an ex-vivo porcine deflated lung indicated feasibility and reliability of using the developed biomechanical model to predict tumor motion during brachytherapy. For external beam radiotherapy, in-silico studies indicated very significant impact of considering the lung tissue incompressibility on the accuracy of predicting tumor motion. Furthermore, ex-vivo porcine lung experiments demonstrated the capability and reliability of the proposed approach for predicting tumor motion as a function of cyclic time. As such, the proposed models have a good potential to be incorporated effectively in computer assisted lung radiotherapy treatment systems