Statistical Characterisation of Speckle in Clinical Echocardiographic Images with Pearson Family of Distributions

The statistical characterisation of gray level distribution of echocardiographic images is commonly done in terms of unimodal probability densities such as Rayleigh, Gamma, Weibull, Nakagami, and Lognormal. Amongst these distributions, the Gamma density is found to provide better empirical model that fits to real data sets. We propose to extend the class of probability distributions by exploring Pearson family to characterise blood and tissue in echocardiographic images. It is found that Pearson Type I characterises the tissue regions whereas Type I, Type IV and Type VI classify blood regions. The statistical measures viz. Jensen-Shannon (JS) divergence and Kolmogorov-Smirnov (KS) statistics reveal that Pearson family of curves outperforms the Gamma distribution.Defence Science Journal, 2011, 61(5), pp.473-478, DOI:http://dx.doi.org/10.14429/dsj.62.116

CAROTID PLAQUE AND INTIMA MEDIA THICKNESS IN THE ASSESSMENT OF CARDIOVASCULAR RISK

Dramatic advances have been demonstrated over the past decade in the prevention and treatment of cardiovascular disease (CVD). Despite these major strides, CVD continues to be our nation's most significant cause of morbidity and mortality. The risk status of persons without known CVD varies greatly and thus requires a range of intense screening and interventions. This dissertation focuses on subclinical CVD measures as well as a new methodology that will improve the evaluation of CVD in clinical trials and eventually improve primary prevention of CVD. There are three related projects in this dissertation, each of which uses noninvasive subclinical methodologies to assess cardiovascular risk. The first project focuses on a high-risk population, the elderly, and evaluates the association calcified carotid plaques with cardiovascular outcomes. Carotid plaque characterization is a new focus of research across the nation and what makes one plaque more dangerous than another is unclear. We do know that as plaques age, the plaques often become more complicated and often calcify. However, the significance of calcification in the carotid arteries is poorly understood. In this project, I assess if carotid calcification is predictive of cardiovascular outcomes.The second project focuses on another high cardiovascular risk population, women systemic lupus erythematosus (SLE). Women with SLE have a significantly high risk of myocardial infarction compared to women without SLE. The role that lupus-related risk factors have in cardiovascular disease progression above the traditional risk factors is unclear. Using carotid ultrasound, associations are evaluated between intima-media thickness and plaque with both cardiovascular and SLE-specific risk factors. The third and final project is the development of a protocol that will allow new computerized assessment of carotid artery plaques. Over the past decade both ultrasound technology and computerized assessment tools have improved. This creates opportunity for improved plaque assessment in vivo. This methodology characterized plaque components, possibly identifying plaques that may be dangerous. Plaque characterization software is now available for use with ultrasound and I have developed the protocols to execute this technique in the Ultrasound Research Laboratory. The final project outlines the software and testing process development, staff training, worksheet design, quality control processes, and a pilot study to evaluate the reproducibility of the measure. This research will contribute to public health through new cardiovascular risk assessment techniques and may lead to improved primary prevention and research methods

First-order statistical speckle models improve robustness and reproducibility of contrast-enhanced ultrasound perfusion estimates

Contrast-enhanced ultrasound (CEUS) permits the quantification and monitoring of adaptive tumor responses in the face of anti-angiogenic treatment, with the goal of informing targeted therapy. However, conventional CEUS image analysis relies on mean signal intensity as an estimate of tracer concentration in indicator-dilution modeling. This discounts additional information that may be available from the first-order speckle statistics in a CEUS image. Heterogeneous vascular networks, typical of tumor-induced angiogenesis, lead to heterogeneous contrast enhancement of the imaged tumor cross-section. To address this, a linear (B-mode) processing approach was developed to quantify the change in the first-order speckle statistics of B-mode cine loops due to the incursion of microbubbles. The technique, named the EDoF (effective degrees of freedom) method, was developed on tumor bearing mice (MDA-MB-231LN mammary fat pad inoculation) and evaluated using nonlinear (two-pulse amplitude modulated) contrast microbubble-specific images. To improve the potential clinical applicability of the technique, a second-generation compound probability density function for the statistics of two-pulse amplitude modulated contrast-enhanced ultrasound images was developed. The compound technique was tested in an antiangiogenic drug trial (bevacizumab) on tumor bearing mice (MDA-MB-231LN), and evaluated with gold-standard histology and contrast-enhanced X-ray computed tomography. The compound statistical model could more accurately discriminate anti-VEGF treated tumors from untreated tumors than conventional CEUS image. The technique was then applied to a rapid patient-derived xenograft (PDX) model of renal cell carcinoma (RCC) in the chorioallantoic membrane (CAM) of chicken embryos. The ultimate goal of the PDX model is to screen RCC patients for de novo sunitinib resistance. The analysis of the first-order speckle statistics of contrast-enhanced ultrasound cine loops provides more robust and reproducible estimates of tumor blood perfusion than conventional image analysis. Theoretically this form of analysis could quantify perfusion heterogeneity and provide estimates of vascular fractal dimension, but further work is required to determine what physiological features influence these measures. Treatment sensitivity matrices, which combine vascular measures from CEUS and power Doppler, may be suitable for screening of de novo sunitinib resistance in patients diagnosed with renal cell carcinoma. Further studies are required to assess whether this protocol can be predictive of patient outcome

Analysis of contrast-enhanced medical images.

Early detection of human organ diseases is of great importance for the accurate diagnosis and institution of appropriate therapies. This can potentially prevent progression to end-stage disease by detecting precursors that evaluate organ functionality. In addition, it also assists the clinicians for therapy evaluation, tracking diseases progression, and surgery operations. Advances in functional and contrast-enhanced (CE) medical images enabled accurate noninvasive evaluation of organ functionality due to their ability to provide superior anatomical and functional information about the tissue-of-interest. The main objective of this dissertation is to develop a computer-aided diagnostic (CAD) system for analyzing complex data from CE magnetic resonance imaging (MRI). The developed CAD system has been tested in three case studies: (i) early detection of acute renal transplant rejection, (ii) evaluation of myocardial perfusion in patients with ischemic heart disease after heart attack; and (iii), early detection of prostate cancer. However, developing a noninvasive CAD system for the analysis of CE medical images is subject to multiple challenges, including, but are not limited to, image noise and inhomogeneity, nonlinear signal intensity changes of the images over the time course of data acquisition, appearances and shape changes (deformations) of the organ-of-interest during data acquisition, determination of the best features (indexes) that describe the perfusion of a contrast agent (CA) into the tissue. To address these challenges, this dissertation focuses on building new mathematical models and learning techniques that facilitate accurate analysis of CAs perfusion in living organs and include: (i) accurate mathematical models for the segmentation of the object-of-interest, which integrate object shape and appearance features in terms of pixel/voxel-wise image intensities and their spatial interactions; (ii) motion correction techniques that combine both global and local models, which exploit geometric features, rather than image intensities to avoid problems associated with nonlinear intensity variations of the CE images; (iii) fusion of multiple features using the genetic algorithm. The proposed techniques have been integrated into CAD systems that have been tested in, but not limited to, three clinical studies. First, a noninvasive CAD system is proposed for the early and accurate diagnosis of acute renal transplant rejection using dynamic contrast-enhanced MRI (DCE-MRI). Acute rejection–the immunological response of the human immune system to a foreign kidney–is the most sever cause of renal dysfunction among other diagnostic possibilities, including acute tubular necrosis and immune drug toxicity. In the U.S., approximately 17,736 renal transplants are performed annually, and given the limited number of donors, transplanted kidney salvage is an important medical concern. Thus far, biopsy remains the gold standard for the assessment of renal transplant dysfunction, but only as the last resort because of its invasive nature, high cost, and potential morbidity rates. The diagnostic results of the proposed CAD system, based on the analysis of 50 independent in-vivo cases were 96% with a 95% confidence interval. These results clearly demonstrate the promise of the proposed image-based diagnostic CAD system as a supplement to the current technologies, such as nuclear imaging and ultrasonography, to determine the type of kidney dysfunction. Second, a comprehensive CAD system is developed for the characterization of myocardial perfusion and clinical status in heart failure and novel myoregeneration therapy using cardiac first-pass MRI (FP-MRI). Heart failure is considered the most important cause of morbidity and mortality in cardiovascular disease, which affects approximately 6 million U.S. patients annually. Ischemic heart disease is considered the most common underlying cause of heart failure. Therefore, the detection of the heart failure in its earliest forms is essential to prevent its relentless progression to premature death. While current medical studies focus on detecting pathological tissue and assessing contractile function of the diseased heart, this dissertation address the key issue of the effects of the myoregeneration therapy on the associated blood nutrient supply. Quantitative and qualitative assessment in a cohort of 24 perfusion data sets demonstrated the ability of the proposed framework to reveal regional perfusion improvements with therapy, and transmural perfusion differences across the myocardial wall; thus, it can aid in follow-up on treatment for patients undergoing the myoregeneration therapy. Finally, an image-based CAD system for early detection of prostate cancer using DCE-MRI is introduced. Prostate cancer is the most frequently diagnosed malignancy among men and remains the second leading cause of cancer-related death in the USA with more than 238,000 new cases and a mortality rate of about 30,000 in 2013. Therefore, early diagnosis of prostate cancer can improve the effectiveness of treatment and increase the patient’s chance of survival. Currently, needle biopsy is the gold standard for the diagnosis of prostate cancer. However, it is an invasive procedure with high costs and potential morbidity rates. Additionally, it has a higher possibility of producing false positive diagnosis due to relatively small needle biopsy samples. Application of the proposed CAD yield promising results in a cohort of 30 patients that would, in the near future, represent a supplement of the current technologies to determine prostate cancer type. The developed techniques have been compared to the state-of-the-art methods and demonstrated higher accuracy as shown in this dissertation. The proposed models (higher-order spatial interaction models, shape models, motion correction models, and perfusion analysis models) can be used in many of today’s CAD applications for early detection of a variety of diseases and medical conditions, and are expected to notably amplify the accuracy of CAD decisions based on the automated analysis of CE images

Injury induced neuroplasticity and cell specific targeting of the lumbar enlargement for gene therapy.

This dissertation is an examination of spinal cord injury induced neuroplasticity and tests whether noninvasive gene therapy can successfully target neurons in the lumbar spinal cord. It begins with an overview of neural control of locomotion and a brief summary of therapeutics that are used and/or in development for treating spinal anatomically characterize s subset of neurons in the spinal cord, long ascending propriospinal neurons, that are involved in interlimb coordination. Characterization of these neurons allows for subsequent evaluation of their potential involvement in injury induced neuroplasticity. This dissertation is divided into five chapters, covering spinal cord injury and therapeutics. Chapter One gives background on locomotor control, propriospinal neurons, spinal cord injury, and therapeutics. Chapter Two develops and characterizes viral tracing methods for spinal cord anatomy. Chapter Three then uses these methods to characterize long ascending propriospinal neurons and evaluate their involvement in injury induced plasticity. Chapter Four then focuses on the development of noninvasive delivery of gene transfer to the lumbar enlargement. This involves optimizing focused ultrasound and intravenous microbubble delivery to focally and transiently permeabilize the blood spinal cord barrier of the lumbar spinal cord. This optimization then allows for successful gene transfer in neurons in the lumbar spinal cord following intravenous delivery of viral vector. Lasty, Chapter Five discusses the implications for all of these findings and how these findings have contributed to our understanding spinal cord anatomy and injury, and how the proof-of-concept in Chapter 4 provides a promising new avenue for spinal cord injury therapeutics

Navigation system based in motion tracking sensor for percutaneous renal access

Tese de Doutoramento em Engenharia BiomédicaMinimally-invasive kidney interventions are daily performed to diagnose and treat several renal diseases. Percutaneous renal access (PRA) is an essential but challenging stage for most of these procedures, since its outcome is directly linked to the physician’s ability to precisely visualize and reach the anatomical target. Nowadays, PRA is always guided with medical imaging assistance, most frequently using X-ray based imaging (e.g. fluoroscopy). Thus, radiation on the surgical theater represents a major risk to the medical team, where its exclusion from PRA has a direct impact diminishing the dose exposure on both patients and physicians. To solve the referred problems this thesis aims to develop a new hardware/software framework to intuitively and safely guide the surgeon during PRA planning and puncturing. In terms of surgical planning, a set of methodologies were developed to increase the certainty of reaching a specific target inside the kidney. The most relevant abdominal structures for PRA were automatically clustered into different 3D volumes. For that, primitive volumes were merged as a local optimization problem using the minimum description length principle and image statistical properties. A multi-volume Ray Cast method was then used to highlight each segmented volume. Results show that it is possible to detect all abdominal structures surrounding the kidney, with the ability to correctly estimate a virtual trajectory. Concerning the percutaneous puncturing stage, either an electromagnetic or optical solution were developed and tested in multiple in vitro, in vivo and ex vivo trials. The optical tracking solution aids in establishing the desired puncture site and choosing the best virtual puncture trajectory. However, this system required a line of sight to different optical markers placed at the needle base, limiting the accuracy when tracking inside the human body. Results show that the needle tip can deflect from its initial straight line trajectory with an error higher than 3 mm. Moreover, a complex registration procedure and initial setup is needed. On the other hand, a real-time electromagnetic tracking was developed. Hereto, a catheter was inserted trans-urethrally towards the renal target. This catheter has a position and orientation electromagnetic sensor on its tip that function as a real-time target locator. Then, a needle integrating a similar sensor is used. From the data provided by both sensors, one computes a virtual puncture trajectory, which is displayed in a 3D visualization software. In vivo tests showed a median renal and ureteral puncture times of 19 and 51 seconds, respectively (range 14 to 45 and 45 to 67 seconds). Such results represent a puncture time improvement between 75% and 85% when comparing to state of the art methods. 3D sound and vibrotactile feedback were also developed to provide additional information about the needle orientation. By using these kind of feedback, it was verified that the surgeon tends to follow a virtual puncture trajectory with a reduced amount of deviations from the ideal trajectory, being able to anticipate any movement even without looking to a monitor. Best results show that 3D sound sources were correctly identified 79.2 ± 8.1% of times with an average angulation error of 10.4º degrees. Vibration sources were accurately identified 91.1 ± 3.6% of times with an average angulation error of 8.0º degrees. Additionally to the EMT framework, three circular ultrasound transducers were built with a needle working channel. One explored different manufacture fabrication setups in terms of the piezoelectric materials, transducer construction, single vs. multi array configurations, backing and matching material design. The A-scan signals retrieved from each transducer were filtered and processed to automatically detect reflected echoes and to alert the surgeon when undesirable anatomical structures are in between the puncture path. The transducers were mapped in a water tank and tested in a study involving 45 phantoms. Results showed that the beam cross-sectional area oscillates around the ceramics radius and it was possible to automatically detect echo signals in phantoms with length higher than 80 mm. Hereupon, it is expected that the introduction of the proposed system on the PRA procedure, will allow to guide the surgeon through the optimal path towards the precise kidney target, increasing surgeon’s confidence and reducing complications (e.g. organ perforation) during PRA. Moreover, the developed framework has the potential to make the PRA free of radiation for both patient and surgeon and to broad the use of PRA to less specialized surgeons.Intervenções renais minimamente invasivas são realizadas diariamente para o tratamento e diagnóstico de várias doenças renais. O acesso renal percutâneo (ARP) é uma etapa essencial e desafiante na maior parte destes procedimentos. O seu resultado encontra-se diretamente relacionado com a capacidade do cirurgião visualizar e atingir com precisão o alvo anatómico. Hoje em dia, o ARP é sempre guiado com recurso a sistemas imagiológicos, na maior parte das vezes baseados em raios-X (p.e. a fluoroscopia). A radiação destes sistemas nas salas cirúrgicas representa um grande risco para a equipa médica, aonde a sua remoção levará a um impacto direto na diminuição da dose exposta aos pacientes e cirurgiões. De modo a resolver os problemas existentes, esta tese tem como objetivo o desenvolvimento de uma framework de hardware/software que permita, de forma intuitiva e segura, guiar o cirurgião durante o planeamento e punção do ARP. Em termos de planeamento, foi desenvolvido um conjunto de metodologias de modo a aumentar a eficácia com que o alvo anatómico é alcançado. As estruturas abdominais mais relevantes para o procedimento de ARP, foram automaticamente agrupadas em volumes 3D, através de um problema de optimização global com base no princípio de “minimum description length” e propriedades estatísticas da imagem. Por fim, um procedimento de Ray Cast, com múltiplas funções de transferência, foi utilizado para enfatizar as estruturas segmentadas. Os resultados mostram que é possível detetar todas as estruturas abdominais envolventes ao rim, com a capacidade para estimar corretamente uma trajetória virtual. No que diz respeito à fase de punção percutânea, foram testadas duas soluções de deteção de movimento (ótica e eletromagnética) em múltiplos ensaios in vitro, in vivo e ex vivo. A solução baseada em sensores óticos ajudou no cálculo do melhor ponto de punção e na definição da melhor trajetória a seguir. Contudo, este sistema necessita de uma linha de visão com diferentes marcadores óticos acoplados à base da agulha, limitando a precisão com que a agulha é detetada no interior do corpo humano. Os resultados indicam que a agulha pode sofrer deflexões à medida que vai sendo inserida, com erros superiores a 3 mm. Por outro lado, foi desenvolvida e testada uma solução com base em sensores eletromagnéticos. Para tal, um cateter que integra um sensor de posição e orientação na sua ponta, foi colocado por via trans-uretral junto do alvo renal. De seguida, uma agulha, integrando um sensor semelhante, é utilizada para a punção percutânea. A partir da diferença espacial de ambos os sensores, é possível gerar uma trajetória de punção virtual. A mediana do tempo necessário para puncionar o rim e ureter, segundo esta trajetória, foi de 19 e 51 segundos, respetivamente (variações de 14 a 45 e 45 a 67 segundos). Estes resultados representam uma melhoria do tempo de punção entre 75% e 85%, quando comparados com o estado da arte dos métodos atuais. Além do feedback visual, som 3D e feedback vibratório foram explorados de modo a fornecer informações complementares da posição da agulha. Verificou-se que com este tipo de feedback, o cirurgião tende a seguir uma trajetória de punção com desvios mínimos, sendo igualmente capaz de antecipar qualquer movimento, mesmo sem olhar para o monitor. Fontes de som e vibração podem ser corretamente detetadas em 79,2 ± 8,1% e 91,1 ± 3,6%, com erros médios de angulação de 10.4º e 8.0 graus, respetivamente. Adicionalmente ao sistema de navegação, foram também produzidos três transdutores de ultrassom circulares com um canal de trabalho para a agulha. Para tal, foram exploradas diferentes configurações de fabricação em termos de materiais piezoelétricos, transdutores multi-array ou singulares e espessura/material de layers de suporte. Os sinais originados em cada transdutor foram filtrados e processados de modo a detetar de forma automática os ecos refletidos, e assim, alertar o cirurgião quando existem variações anatómicas ao longo do caminho de punção. Os transdutores foram mapeados num tanque de água e testados em 45 phantoms. Os resultados mostraram que o feixe de área em corte transversal oscila em torno do raio de cerâmica, e que os ecos refletidos são detetados em phantoms com comprimentos superiores a 80 mm. Desta forma, é expectável que a introdução deste novo sistema a nível do ARP permitirá conduzir o cirurgião ao longo do caminho de punção ideal, aumentado a confiança do cirurgião e reduzindo possíveis complicações (p.e. a perfuração dos órgãos). Além disso, de realçar que este sistema apresenta o potencial de tornar o ARP livre de radiação e alarga-lo a cirurgiões menos especializados.The present work was only possible thanks to the support by the Portuguese Science and Technology Foundation through the PhD grant with reference SFRH/BD/74276/2010 funded by FCT/MEC (PIDDAC) and by Fundo Europeu de Desenvolvimento Regional (FEDER), Programa COMPETE - Programa Operacional Factores de Competitividade (POFC) do QREN