Parallel sets and morphological measurements of CT images of soil pore structure in a vineyard

Important physical and biological processes in soil-plant-microbial systems are dominated by the geometry of soil pore space, and a correct model of this geometry is critical for understanding them. We analyze the geometry of soil pore space with the X-ray computed tomography (CT) of intact soil columns. We present here some preliminary results of our investigation on Minkowski functionals of parallel sets to characterize soil structure. We also show how the evolution of Minkowski morphological measurements of parallel sets may help to characterize the influence of conventional tillage and permanent cover crop of resident vegetation on soil structure in a Spanish Mediterranean vineyard.

Volume, Surface, Connectivity and Size Distribution of Soil Pore Space in CT Images: Comparison of Samples at Different Depths from Nearby Natural and Tillage Areas

The study of soil structure, i.e., the pores, is of vital importance in different fields of science and technology. Total pore volume (porosity), pore surface, pore connectivity and pore size distribution are some (probably the most important) of the geometric measurements of pore space. The technology of X-ray computed tomography allows us to obtain 3D images of the inside of a soil sample enabling study of the pores without disturbing the samples. In this work we performed a set of geometrical measures, some of them from mathematical morphology, to assess and quantify any possible difference that tillage may have caused on the soil. We compared samples from tilled soil with samples from a soil with natural vegetation taken in a very close area. Our results show that the main differences between these two groups of samples are total surface area and pore connectivity per unit pore volume.

Fractal parameters of pore space from CT images of soils under contrasting management practices

Soil structure plays an important role in flow and transport phenomena, and a quantitative characterization of the spatial heterogeneity of the pore space geometry is beneficial for prediction of soil physical properties. Morphological features such as pore-size distribution, pore space volume or pore?solid surface can be altered by different soil management practices. Irregularity of these features and their changes can be described using fractal geometry. In this study, we focus primarily on the characterization of soil pore space as a 3D geometrical shape by fractal analysis and on the ability of fractal dimensions to differentiate between two a priori different soil structures. We analyze X-ray computed tomography (CT) images of soils samples from two nearby areas with contrasting management practices. Within these two different soil systems, samples were collected from three depths. Fractal dimensions of the pore-size distributions were different depending on soil use and averaged values also differed at each depth. Fractal dimensions of the volume and surface of the pore space were lower in the tilled soil than in the natural soil but their standard deviations were higher in the former as compared to the latter. Also, it was observed that soil use was a factor that had a statistically significant effect on fractal parameters. Fractal parameters provide useful complementary information about changes in soil structure due to changes in soil management. Read More: http://www.worldscientific.com/doi/abs/10.1142/S0218348X14400118?queryID=%24%7BresultBean.queryID%7D&

Efeitos da remoção do disco e cartilagem articular no crescimento e microarquitetura óssea da mandíbula de ratos: análise por microtomografia

Alterações na articulação emporomandibular (ATM) comumente geram desequilíbrios musculares que estão associados à alterações no tecido ósseo. Esta articulação pode sofrer a influência de traumas, fatores congênitos ou desordens de crescimento. Estudos sobre alterações de crescimento do complexo maxilomandibular decorrentes de problemas da ATM são escassos. O objetivo deste trabalho foi avaliar por meio da microtomografia os efeitos da remoção do disco articular e a remoção conjugada do disco e cartilagem articular no crescimento e na microarquitetura óssea da mandíbula de ratos. Trinta ratos da raça Wistar com um mês de idade foram divididos em três grupos: CTR (controle operado); RD (remoção de disco articular) e RDC (remoção conjugada do disco e cartilagem articular). Apenas o lado direito foi operado; o lado esquerdo permaneceu intacto. Após dois meses de acompanhamento, os ratos foram sacrificados e as hemimandíbulas escaneadas em microtomógrafo A remoção do disco articular e a remoção conjugada do disco e cartilagem articular alteram o volume e microestrutura do osso trabecular da mandíbula de ratos jovens. Estas duas intervenções provocaram uma queda na qualidade de parâmetros da microestrutura do trabeculado do processo angular e diminuição do crescimento da hemimandíbula do lado operado.

Aproximação facial forense: relação entre as estruturas ósseas e a porção inferior do nariz

Este trabalho teve o propósito de avaliar a relação entre a cartilagem alar e a abertura piriforme a partir de imagens de tomografia computadorizada cone-beam e a relação do nariz com o padrão esquelético vertical da face. A pesquisa foi realizada com 96 imagens de indivíduos de ambos os sexos (49 masculino e 47 feminino), com idades entre 18 e 65 anos classificados de acordo com sua tipologia facial. Para a realização das marcações e mensurações de interesse foi utilizado o software OsiriX. A tipologia facial foi acessada através de três metodologias: índice facial, ângulo goníaco e proporção entre as alturas faciais. No corte axial da imagem de tomografia computadorizada, foram aferidas quatro grandezas lineares: largura do nariz externo, distância entre as inserções alares, extensão da base da abertura piriforme e máxima largura da abertura piriforme. Todas as grandezas foram mensuradas por dois examinadores em concordância. A análise dos resultados obtidos mostrou que há correlação entre as larguras do nariz externo e da abertura piriforme e entre a base da cavidade nasal e a distância entre as inserções alares. A largura da abertura piriforme aumenta proporcionalmente mais que a do nariz. Além disso, a tipologia facial longa associa-se à largura nasal.

Using two dimensional sectional distributions to infer three dimensional volumetric distributions - Validation using tomography

The critical process parameter for mineral separation is the degree of mineral liberation achieved by comminution. The degree of liberation provides an upper limit of efficiency for any physical separation process. The standard approach to measuring mineral liberation uses mineralogical analysis based two-dimensional sections of particles which may be acquired using a scanning electron microscope and back-scatter electron analysis or from an analysis of an image acquired using an optical microscope. Over the last 100 years, mathematical techniques have been developed to use this two dimensional information to infer three-dimensional information about the particles. For mineral processing, a particle that contains more than one mineral (a composite particle) may appear to be liberated (contain only one mineral) when analysed using only its revealed particle section. The mathematical techniques used to interpret three-dimensional information belong, to a branch of mathematics called stereology. However methods to obtain the full mineral liberation distribution of particles from particle sections are relatively new. To verify these adjustment methods, we require an experimental method which can accurately measure both sectional and three dimensional properties. Micro Cone Beam Tomography provides such a method for suitable particles and hence, provides a way to validate methods used to convert two-dimensional measurements to three dimensional estimates. For this study ore particles from a well-characterised sample were subjected to conventional mineralogical analysis (using particle sections) to estimate three-dimensional properties of the particles. A subset of these particles was analysed using a micro-cone beam tomograph. This paper presents a comparison of the three-dimensional properties predicted from measured two-dimensional sections with the measured three-dimensional properties.

EFEITOS ESQUELÉTICOS DA TRAÇÃO REVERSA DA MAXILA UTILIZANDO TOMOGRAFIA COMPUTADORIZADA

Este estudo avaliou os efeitos esqueléticos da tração reversa da maxila utilizando imagens 2D (telerradiografia lateral) geradas a partir da tomografia de feixe cônico (imagens 3D). A amostra foi composta por 20 crianças (15 do gênero feminino, e 5 do masculino), com idade variando de 5,6 a 10,7 anos que apresentavam má-oclusão de Classe III de Angle. A tomografia foi realizada antes do tratamento (T1) e logo após o tratamento (T2). O tratamento foi realizado por meio da tração reversa da maxila utilizando-se o aparelho expansor Hyrax associado à máscara facial individualizada, com força de 600 a 800g de cada lado, durante 14 horas por dia. A correção da relação de caninos em Classe I ou com sua sobrecorreção em Classe II foi obtida após 4 a 8 meses de tratamento. Para verificar o erro sistemático e casual foi utilizado o teste t pareado e a fórmula de Dahlberg, respectivamente. O teste t pareado (p<0,05) mostrou diferença significante entre as medidas cefalométricas obtidas em T1 e T2. Na maxila houve aumento do SNA 2,2°, A-Nperp 1,47mm e em Co-A 2,58mm. Na mandíbula, SNB diminuiu -0,54° e P-Nperp, -1,45mm, enquanto Co-Gn aumentou 1,04mm. Houve melhora na relação maxilo-mandibular ANB 2,74° e Wits 4,23mm. As variáveis GoGn.SN, Gn.SN, FH.Md, Mx.Md, e AFAI aumentaram demonstrando que houve uma rotação da mandíbula no sentido horário. O plano palatino rotacionou no sentido anti-horário. Pode se concluir que o tratamento de tração reversa da maxila na idade precoce promoveu uma melhora na relação maxilo-mandibular devido a um avanço da maxila e um deslocamento da mandíbula para baixo e para trás.

Development and Optimization of Four-dimensional Magnetic Resonance Imaging (4D-MRI) for Radiation Therapy

A tenet of modern radiotherapy (RT) is to identify the treatment target accurately, following which the high-dose treatment volume may be expanded into the surrounding tissues in order to create the clinical and planning target volumes. Respiratory motion can induce errors in target volume delineation and dose delivery in radiation therapy for thoracic and abdominal cancers. Historically, radiotherapy treatment planning in the thoracic and abdominal regions has used 2D or 3D images acquired under uncoached free-breathing conditions, irrespective of whether the target tumor is moving or not. Once the gross target volume has been delineated, standard margins are commonly added in order to account for motion. However, the generic margins do not usually take the target motion trajectory into consideration. That may lead to under- or over-estimate motion with subsequent risk of missing the target during treatment or irradiating excessive normal tissue. That introduces systematic errors into treatment planning and delivery. In clinical practice, four-dimensional (4D) imaging has been popular in For RT motion management. It provides temporal information about tumor and organ at risk motion, and it permits patient-specific treatment planning. The most common contemporary imaging technique for identifying tumor motion is 4D computed tomography (4D-CT). However, CT has poor soft tissue contrast and it induce ionizing radiation hazard. In the last decade, 4D magnetic resonance imaging (4D-MRI) has become an emerging tool to image respiratory motion, especially in the abdomen, because of the superior soft-tissue contrast. Recently, several 4D-MRI techniques have been proposed, including prospective and retrospective approaches. Nevertheless, 4D-MRI techniques are faced with several challenges: 1) suboptimal and inconsistent tumor contrast with large inter-patient variation; 2) relatively low temporal-spatial resolution; 3) it lacks a reliable respiratory surrogate. In this research work, novel 4D-MRI techniques applying MRI weightings that was not used in existing 4D-MRI techniques, including T2/T1-weighted, T2-weighted and Diffusion-weighted MRI were investigated. A result-driven phase retrospective sorting method was proposed, and it was applied to image space as well as k-space of MR imaging. Novel image-based respiratory surrogates were developed, improved and evaluated.

Effective Dose Estimation in Fast-kV Switch Dual Energy Computed Tomography

Purpose

The objective of our study was to test a new approach to approximating organ dose by using the effective energy of the combined 80kV/140kV beam used in fast kV switch dual-energy (DE) computed tomography (CT). The two primary focuses of the study were to first validate experimentally the dose equivalency between MOSFET and ion chamber (as a gold standard) in a fast kV switch DE environment, and secondly to estimate effective dose (ED) of DECT scans using MOSFET detectors and an anthropomorphic phantom.

Materials and Methods

A GE Discovery 750 CT scanner was employed using a fast-kV switch abdomen/pelvis protocol alternating between 80 kV and 140 kV. The specific aims of our study were to (1) Characterize the effective energy of the dual energy environment; (2) Estimate the f-factor for soft tissue; (3) Calibrate the MOSFET detectors using a beam with effective energy equal to the combined DE environment; (4) Validate our calibration by using MOSFET detectors and ion chamber to measure dose at the center of a CTDI body phantom; (5) Measure ED for an abdomen/pelvis scan using an anthropomorphic phantom and applying ICRP 103 tissue weighting factors; and (6) Estimate ED using AAPM Dose Length Product (DLP) method. The effective energy of the combined beam was calculated by measuring dose with an ion chamber under varying thicknesses of aluminum to determine half-value layer (HVL).

Results

The effective energy of the combined dual-energy beams was found to be 42.8 kV. After calibration, tissue dose in the center of the CTDI body phantom was measured at 1.71 ± 0.01 cGy using an ion chamber, and 1.73±0.04 and 1.69±0.09 using two separate MOSFET detectors. This result showed a -0.93% and 1.40 % difference, respectively, between ion chamber and MOSFET. ED from the dual-energy scan was calculated as 16.49 ± 0.04 mSv by the MOSFET method and 14.62 mSv by the DLP method.

Iterative reconstruction as a novel method of radiation dose reduction at computed tomography in patients with Crohn's disease

Prior work of our research group, that quantified the alarming levels of radiation dose to patients with Crohn’s disease from medical imaging and the notable shift towards CT imaging making these patients an at risk group, provided context for this work. CT delivers some of the highest doses of ionising radiation in diagnostic radiology. Once a medical imaging examination is deemed justified, there is an onus on the imaging team to endeavour to produce diagnostic quality CT images at the lowest possible radiation dose to that patient. The fundamental limitation with conventional CT raw data reconstruction was the inherent coupling of administered radiation dose with observed image noise – the lower the radiation dose, the noisier the image. The renaissance, rediscovery and refinement of iterative reconstruction removes this limitation allowing either an improvement in image quality without increasing radiation dose or maintenance of image quality at a lower radiation dose compared with traditional image reconstruction. This thesis is fundamentally an exercise in optimisation in clinical CT practice with the objectives of assessment of iterative reconstruction as a method for improvement of image quality in CT, exploration of the associated potential for radiation dose reduction, and development of a new split dose CT protocol with the aim of achieving and validating diagnostic quality submillisiever t CT imaging in patients with Crohn’s disease. In this study, we investigated the interplay of user-selected parameters on radiation dose and image quality in phantoms and cadavers, comparing traditional filtered back projection (FBP) with iterative reconstruction algorithms. This resulted in the development of an optimised, refined and appropriate split dose protocol for CT of the abdomen and pelvis in clinical patients with Crohn’s disease allowing contemporaneous acquisition of both modified and conventional dose CT studies. This novel algorithm was then applied to 50 patients with a suspected acute complication of known Crohn’s disease and the raw data reconstructed with FBP, adaptive statistical iterative reconstruction (ASiR) and model based iterative reconstruction (MBIR). Conventional dose CT images with FBP reconstruction were used as the reference standard with which the modified dose CT images were compared in terms of radiation dose, diagnostic findings and image quality indices. As there are multiple possible user-selected strengths of ASiR available, these were compared in terms of image quality to determine the optimal strength for this modified dose CT protocol. Modified dose CT images with MBIR were also compared with contemporaneous abdominal radiograph, where performed, in terms of diagnostic yield and radiation dose. Finally, attenuation measurements in organs, tissues, etc. with each reconstruction algorithm were compared to assess for preservation of tissue characterisation capabilities. In the phantom and cadaveric models, both forms of iterative reconstruction examined (ASiR and MBIR) were superior to FBP across a wide variety of imaging protocols, with MBIR superior to ASiR in all areas other than reconstruction speed. We established that ASiR appears to work to a target percentage noise reduction whilst MBIR works to a target residual level of absolute noise in the image. Modified dose CT images reconstructed with both ASiR and MBIR were non-inferior to conventional dose CT with FBP in terms of diagnostic findings, despite reduced subjective and objective indices of image quality. Mean dose reductions of 72.9-73.5% were achieved with the modified dose protocol with a mean effective dose of 1.26mSv. MBIR was again demonstrated superior to ASiR in terms of image quality. The overall optimal ASiR strength for the modified dose protocol used in this work is ASiR 80%, as this provides the most favourable balance of peak subjective image quality indices with less objective image noise than the corresponding conventional dose CT images reconstructed with FBP. Despite guidelines to the contrary, abdominal radiographs are still often used in the initial imaging of patients with a suspected complication of Crohn’s disease. We confirmed the superiority of modified dose CT with MBIR over abdominal radiographs at comparable doses in detection of Crohn’s disease and non-Crohn’s disease related findings. Finally, we demonstrated (in phantoms, cadavers and in vivo) that attenuation values do not change significantly across reconstruction algorithms meaning preserved tissue characterisation capabilities with iterative reconstruction. Both adaptive statistical and model based iterative reconstruction algorithms represent feasible methods of facilitating acquisition diagnostic quality CT images of the abdomen and pelvis in patients with Crohn’s disease at markedly reduced radiation doses. Our modified dose CT protocol allows dose savings of up to 73.5% compared with conventional dose CT, meaning submillisievert imaging is possible in many of these patients.

Open-source virtual bronchoscopy for image guided navigation

This thesis describes the development of an open-source system for virtual bronchoscopy used in combination with electromagnetic instrument tracking. The end application is virtual navigation of the lung for biopsy of early stage cancer nodules. The open-source platform 3D Slicer was used for creating freely available algorithms for virtual bronchscopy. Firstly, the development of an open-source semi-automatic algorithm for prediction of solitary pulmonary nodule malignancy is presented. This approach may help the physician decide whether to proceed with biopsy of the nodule. The user-selected nodule is segmented in order to extract radiological characteristics (i.e., size, location, edge smoothness, calcification presence, cavity wall thickness) which are combined with patient information to calculate likelihood of malignancy. The overall accuracy of the algorithm is shown to be high compared to independent experts' assessment of malignancy. The algorithm is also compared with two different predictors, and our approach is shown to provide the best overall prediction accuracy. The development of an airway segmentation algorithm which extracts the airway tree from surrounding structures on chest Computed Tomography (CT) images is then described. This represents the first fundamental step toward the creation of a virtual bronchoscopy system. Clinical and ex-vivo images are used to evaluate performance of the algorithm. Different CT scan parameters are investigated and parameters for successful airway segmentation are optimized. Slice thickness is the most affecting parameter, while variation of reconstruction kernel and radiation dose is shown to be less critical. Airway segmentation is used to create a 3D rendered model of the airway tree for virtual navigation. Finally, the first open-source virtual bronchoscopy system was combined with electromagnetic tracking of the bronchoscope for the development of a GPS-like system for navigating within the lungs. Tools for pre-procedural planning and for helping with navigation are provided. Registration between the lungs of the patient and the virtually reconstructed airway tree is achieved using a landmark-based approach. In an attempt to reduce difficulties with registration errors, we also implemented a landmark-free registration method based on a balanced airway survey. In-vitro and in-vivo testing showed good accuracy for this registration approach. The centreline of the 3D airway model is extracted and used to compensate for possible registration errors. Tools are provided to select a target for biopsy on the patient CT image, and pathways from the trachea towards the selected targets are automatically created. The pathways guide the physician during navigation, while distance to target information is updated in real-time and presented to the user. During navigation, video from the bronchoscope is streamed and presented to the physician next to the 3D rendered image. The electromagnetic tracking is implemented with 5 DOF sensing that does not provide roll rotation information. An intensity-based image registration approach is implemented to rotate the virtual image according to the bronchoscope's rotations. The virtual bronchoscopy system is shown to be easy to use and accurate in replicating the clinical setting, as demonstrated in the pre-clinical environment of a breathing lung method. Animal studies were performed to evaluate the overall system performance.

Utilisation de la tomodensitométrie à deux énergies pour le calcul de dose en curiethérapie à bas débit

Dans la pratique actuelle de la curiethérapie à bas débit, l'évaluation de la dose dans la prostate est régie par le protocole défini dans le groupe de travail 43 (TG-43) de l'American Association of Physicists in Medicine. Ce groupe de travail suppose un patient homogène à base d'eau de même densité et néglige les changements dans l'atténuation des photons par les sources de curiethérapie. En considérant ces simplifications, les calculs de dose se font facilement à l'aide d'une équation, indiquée dans le protocole. Bien que ce groupe de travail ait contribué à l'uniformisation des traitements en curiethérapie entre les hôpitaux, il ne décrit pas adéquatement la distribution réelle de la dose dans le patient. La publication actuelle du TG-186 donne des recommandations pour étudier des distributions de dose plus réalistes. Le but de ce mémoire est d'appliquer ces recommandations à partir du TG-186 pour obtenir une description plus réaliste de la dose dans la prostate. Pour ce faire, deux ensembles d'images du patient sont acquis simultanément avec un tomodensitomètre à double énergie (DECT). Les artéfacts métalliques présents dans ces images, causés par les sources d’iode, sont corrigés à l'aide d’un algorithme de réduction d'artefacts métalliques pour DECT qui a été développé dans ce travail. Ensuite, une étude Monte Carlo peut être effectuée correctement lorsque l'image est segmentée selon les différents tissus humains. Cette segmentation est effectuée en évaluant le numéro atomique effectif et la densité électronique de chaque voxel, par étalonnage stoechiométrique propre au DECT, et en y associant le tissu ayant des paramètres physiques similaires. Les résultats montrent des différences dans la distribution de la dose lorsqu'on compare la dose du protocole TG-43 avec celle retrouvée avec les recommandations du TG-186.

Improving the Accuracy of 3-D Reconstruction in Robotic Vision Applications

Thesis (Ph.D.)--University of Washington, 2016-08

Evaluating deviations in prostatectomy patients treated with IMRT

Aim - To evaluate the deviations in prostatectomy patients treated with IMRT in order to calculate appropriate margins to create the PTV. Background - Defining inappropriate margins can lead to underdosing in target volumes and also overdosing in healthy tissues, increasing morbidity. Material and methods - 223 CBCT images used for alignment with the CT planning scan based on bony anatomy were analyzed in 12 patients treated with IMRT following prostatectomy. Shifts of CBCT images were recorded in three directions to calculate the required margin to create PTV. Results and discussion - The mean and standard deviation (SD) values in millimetres were −0.05 ± 1.35 in the LR direction, −0.03 ± 0.65 in the SI direction and −0.02 ± 2.05 the AP direction. The systematic error measured in the LR, SI and AP direction were 1.35 mm, 0.65 mm, and 2.05 mm with a random error of 2.07 mm; 1.45 mm and 3.16 mm, resulting in a PTV margin of 4.82 mm; 2.64 mm, and 7.33 mm, respectively. Conclusion - With IGRT we suggest a margin of 5 mm, 3 mm and 8 mm in the LR, SI and AP direction, respectively, to PTV1 and PTV2. Therefore, this study supports an anisotropic margin expansion to the PTV being the largest expansion in the AP direction and lower in SI.

Oncological patient anxiety in imaging studies: the PET/CT example

The purpose of this study was to investigate the subjective perception of anxiety pre- and post-procedure, and explore the relationship between demographic, clinical variables and cancer patients' anxiety during a positron emission tomography/computed tomography (PET/CT) scan. Two hundred and thirty-two oncological out patients, with clinical indication for performing an (18)F-2-fluoro-2-deoxy-D-glucose ((18)F-FDG) PET/CT scan and attending a nuclear medicine (NM) department, participated in the study. Patients' anxiety and subjective experience of PET/CT were examined using two self-report questionnaires. The pre-procedure questionnaire focused on demographic information, level of knowledge regarding the scan and subjective perception of anxiety before the procedure. The post-procedure questionnaire included the subjective perception anxiety after the procedure, information adequacy and satisfaction with the NM department. The self-reported data indicate that patients were anxious during PET/CT. Furthermore, our data revealed a significant difference between the anxiety pre-procedure and post-procedure (z = -3909, p < 0.05), in which the anxiety pre-procedure has significantly higher values. No significant correlation was found between anxiety and age of the patients, education levels, adequacy of information or satisfaction with the NM Department. Perception of anxiety post-procedure differs between gender (U = 5641, p = 0.033). In conclusion, PET/CT generated anxiety levels in oncological patients, especially before the procedure. Although patients seemed to be satisfied with information delivered by staff and with the NM Department, attention has to be focused on effective interventions strategies that help patients to reduce anxiety.