Optimization of Computer Aided Detection systems: an evolutionary approach

Computer Aided Diagnosis (CAD) systems are designed to aid the radiologist in interpreting medical images. They are usually based on lesion detection and segmentation algorithms whose performance depends on a large number of parameters. While time consuming and sub-optimal, manual adjustment is still widely used to adjust parameter values. Genetic or evolutionary algorithms (GA) are effective optimization methods that mimic biological evolution. Genetic algorithms have been shown to efficiently manage complex search spaces, and can be applied to all kinds of objective functions, including discontinuous, nondifferentiable, or highly nonlinear ones. In this study, we have adopted an evolutionary approach to the problem of parameter optimization. We show that the genetic algorithm is able to effectively converge to a better solution than manual optimization on a case study for digital breast tomosynthesis CAD. Parameter optimization was framed as a constrained optimization problem, where the function to be maximized was defined as weighted sum of sensitivity, false positive rate and segmentation accuracy. A modified Dice coefficient was defined to assess the segmentation quality of individual lesions. Finally, all viable solutions evaluated by the GA were studied by means of exploratory data analysis techniques, such as association rules, to gain useful insight on the strength of the influence of each parameter on overall algorithm performance. We showed that this combination was able to identify multiple ranges of viable solutions with good segmentation accuracy

ADNet : diagnóstico assistido por computador para doença de Alzheimer usando rede neural convolucional 3D com cérebro inteiro

Orientadores: Anderson de Rezende Rocha, Marina WeilerDissertação (mestrado) - Universidade Estadual de Campinas, Instituto de ComputaçãoResumo: Demência por doença de Alzheimer (DA) é uma síndrome clínica caracterizada por múltiplos problemas cognitivos, incluindo dificuldades na memória, funções executivas, linguagem e habilidades visuoespaciais. Sendo a forma mais comum de demência, essa doença mata mais do que câncer de mama e de próstata combinados, além de ser a sexta principal causa de morte nos Estados Unidos. A neuroimagem é uma das áreas de pesquisa mais promissoras para a detecção de biomarcadores estruturais da DA, onde uma técnica não invasiva é usada para capturar uma imagem digital do cérebro, a partir da qual especialistas extraem padrões e características da doença. Nesse contexto, os sistemas de diagnóstico assistido por computador (DAC) são abordagens que visam ajudar médicos e especialistas na interpretação de dados médicos, para fornecer diagnósticos aos pacientes. Em particular, redes neurais convolucionais (RNCs) são um tipo especial de rede neural artificial (RNA), que foram inspiradas em como o sistema visual funciona e, nesse sentido, têm sido cada vez mais utilizadas em tarefas de visão computacional, alcançando resultados impressionantes. Em nossa pesquisa, um dos principais objetivos foi utilizar o que há de mais avançado sobre aprendizagem profunda (por exemplo, RNC) para resolver o difícil problema de identificar biomarcadores estruturais da DA em imagem por ressonância magnética (IRM), considerando três grupos diferentes, ou seja, cognitivamente normal (CN), comprometimento cognitivo leve (CCL) e DA. Adaptamos redes convolucionais com dados fornecidos principalmente pela ADNI e avaliamos no desafio CADDementia, resultando em um cenário mais próximo das condições no mundo real, em que um sistema DAC é usado em um conjunto de dados diferente daquele usado no treinamento. Os principais desafios e contribuições da nossa pesquisa incluem a criação de um sistema de aprendizagem profunda que seja totalmente automático e comparativamente rápido, ao mesmo tempo em que apresenta resultados competitivos, sem usar qualquer conhecimento específico de domínio. Nomeamos nossa melhor arquitetura ADNet (Alzheimer's Disease Network) e nosso melhor método ADNet-DA (ADNet com adaptação de domínio), o qual superou a maioria das submissões no CADDementia, todas utilizando conhecimento prévio da doença, como regiões de interesse específicas do cérebro. A principal razão para não usar qualquer informação da doença em nosso sistema é fazer com que ele aprenda e extraia padrões relevantes de regiões importantes do cérebro automaticamente, que podem ser usados para apoiar os padrões atuais de diagnóstico e podem inclusive auxiliar em novas descobertas para diferentes ou novas doenças. Após explorar uma série de técnicas de visualização para interpretação de modelos, associada à inteligência artificial explicável (XAI), acreditamos que nosso método possa realmente ser empregado na prática médica. Ao diagnosticar pacientes, é possível que especialistas usem a ADNet para gerar uma diversidade de visualizações explicativas para uma determinada imagem, conforme ilustrado em nossa pesquisa, enquanto a ADNet-DA pode ajudar com o diagnóstico. Desta forma, os especialistas podem chegar a uma decisão mais informada e em menos tempoAbstract: Dementia by Alzheimer's disease (AD) is a clinical syndrome characterized by multiple cognitive problems, including difficulties in memory, executive functions, language and visuospatial skills. Being the most common form of dementia, this disease kills more than breast cancer and prostate cancer combined, and it is the sixth leading cause of death in the United States. Neuroimaging is one of the most promising areas of research for early detection of AD structural biomarkers, where a non-invasive technique is used to capture a digital image of the brain, from which specialists extract patterns and features of the disease. In this context, computer-aided diagnosis (CAD) systems are approaches that aim at assisting doctors and specialists in interpretation of medical data to provide diagnoses for patients. In particular, convolutional neural networks (CNNs) are a special kind of artificial neural network (ANN), which were inspired by how the visual system works, and, in this sense, have been increasingly used in computer vision tasks, achieving impressive results. In our research, one of the main goals was bringing to bear what is most advanced in deep learning research (e.g., CNN) to solve the difficult problem of identifying AD structural biomarkers in magnetic resonance imaging (MRI), considering three different groups, namely, cognitively normal (CN), mild cognitive impairment (MCI), and AD. We tailored convolutional networks with data primarily provided by ADNI, and evaluated them on the CADDementia challenge, thus resulting in a scenario very close to the real-world conditions, in which a CAD system is used on a dataset differently from the one used for training. The main challenges and contributions of our research include devising a deep learning system that is both completely automatic and comparatively fast, while also presenting competitive results, without using any domain specific knowledge. We named our best architecture ADNet (Alzheimer's Disease Network), and our best method ADNet-DA (ADNet with domain adaption), which outperformed most of the CADDementia submissions, all of them using prior knowledge from the disease, such as specific regions of interest of the brain. The main reason for not using any information from the disease in our system is to make it automatically learn and extract relevant patterns from important regions of the brain, which can be used to support current diagnosis standards, and may even assist in new discoveries for different or new diseases. After exploring a number of visualization techniques for model interpretability, associated with explainable artificial intelligence (XAI), we believe that our method can be actually employed in medical practice. While diagnosing patients, it is possible for specialists to use ADNet to generate a diversity of explanatory visualizations for a given image, as illustrated in our research, while ADNet-DA can assist with the diagnosis. This way, specialists can come up with a more informed decision and in less timeMestradoCiência da ComputaçãoMestre em Ciência da Computaçã

Computational Methods For The Diagnosis of Rheumatoid Arthritis With Diffuse Optical Tomography

Diffuse optical tomography (DOT) is an imaging technique where near infrared (NIR) photons are used to probe biological tissue. DOT allows for the recovery of three-dimensional maps of tissue optical properties, such as tissue absorption and scattering coefficients. The application of DOT as a tool to aid in the diagnosis of rheumatoid arthritis (RA) is explored in this work. Algorithms for improving the image reconstruction process and for enhancing the clinical value of DOT images are presented in detail. The clinical data considered in this work consists of 99 fingers from subjects with RA and 120 fingers from healthy subjects. DOT scans of the proximal interphalangeal (PIP) joint of each finger is performed with modulation frequencies of 0, 300, and 600 MHz. A computer-aided diagnosis (CAD) framework for extracting heuristic features from DOT images and a method for using these same features to classify each joint as affected or not affected by RA is presented. The framework is applied to the clinical data and results are discussed in detail. Then, an algorithm for recovering the optical properties of biological media using the simplified spherical harmonics (SPN) light propagation model is presented. The computational performance of the algorithm is analyzed and reported. Finally, the SPN reconstruction algorithm is applied to clinical data of subjects with RA and the resulting images are analyzed with the CAD framework. As the first part of the CAD framework, heuristic image features are extracted from the absorption and the scattering coefficient images using multiple compression and dimensionality reduction techniques. Overall, 594 features are extracted from the images of each joint. Then, machine-learning techniques are used to evaluate the ability to discriminate between images of joints with RA and images of healthy joints. An evolution-strategy optimization algorithm is developed to evaluate the classification strength of each feature and to find the multidimensional feature combination that results in optimal classification accuracy. Classification is performed with k-nearest neighbors (KNN), linear (LDA) and quadratic discriminate analysis (QDA), self-organizing maps (SOM), or support vector machines (SVM). Classification accuracy is evaluated based on diagnostic sensitivity and specificity values. Strong evidence is presented that suggest there are clear differences between the tissue optical parameters of joints with RA and joints without RA. It is first shown that data obtained at 600 MHz leads to better classification results than data obtained at 300 and 0 MHz. Analysis of each extracted feature shows that DOT images of subjects with RA are statistically different (p < 0.05) from images of subjects without RA for over 90% of the features. Evidence shows that subjects with RA that do not have detectable signs of erosion, effusion, or synovitis (i.e. asymptomatic subjects) in MRI and US images have optical profiles similar to subjects who do have signs of erosion, effusion, or synovitis; furthermore, both of these cohorts differ from healthy controls subjects. This shows that it may be possible to accurately identify asymptomatic subjects with DOT scans. In contrast, these subjects remain difficult to identify from MRI and US images. The implications of these results are profound, as they suggest it may be possible to identify RA with DOT at an earlier stage compared to standard imaging techniques. Results from the feature-selection algorithm show that the SVM algorithm (with a third order polynomial kernel) achieves 100.0% sensitivity and 97.8% specificity. Lower bounds for these results (at 95.0% confidence level) are 96.4% and 93.8%, respectively. Image features most predictive of RA are from the spatial variation of optical properties and the absolute range in feature values. The optimal classifiers are low dimensional combinations (< 7 features). Robust cross- validation is performed to ensure the generalization of these classification results. The SPN -based reconstruction algorithm uses a reduced-Hessian sequential quadratic programming (rSQP) PDE-constrained optimization approach to maximize computational efficiency. The complex-valued forward model, or frequency domain SPN equations (N = 1, 3), is discretized using the finite-volume method and solved on unstructured computational grids using the restarted GMRES algorithm. The image reconstruction algorithm is presented in detail and its performance benchmarked against the ERT algorithm. The algorithm is subsequently used to recover the absorption and scattering coefficient images of joints scanned in the RA clinical study. While the SPN model is inherently less accurate than the ERT model, it is nevertheless shown that the images obtained with the SP3-based reconstruction algorithm are sufficiently accurate and allow for the diagnosis of RA at clinically relevant sensitivity [87.9% (78.1%, 100.0%)] and specificity [92.9% (84.6%, 100.0%)] values (the 95.0% confidence interval is specified in brackets). In contrast to results obtained with the SP3 model, the images generated with the SP1 algorithm yield significantly lower sensitivity [66.7% (46.6%, 100.0%)] and specificity [81.0% (64.8%, 100.0%)] values. While some numerical accuracy is sacrificed by selecting the SP3 model over the ERT model, the superior computational performance of the SP3 algorithm allows for computation of the absorption and the scattering coefficient images in under 15 minutes and requires less than 200 MB of RAM per finger (compared to the over 180 minutes and over 6 GB of RAM needed by the ERT-based algorithm). Overall, results indicate that the SP3-based reconstruction algorithm provides computational advantages over the ERT-based algorithm without sacrificing significant classification accuracy. In contrast, the SP1 model provides computational advantages compared to the ERT at the expense of classification accuracy. This indicates that the frequency-domain SP3 model is an ideal light propagation model for use in DOT scanning of finger joints with RA. Altogether, the results presented in this dissertation underscore the high potential for DOT to become a clinically useful diagnostic tool. The algorithms and framework developed as part of this dissertation can be directly used on future data to help further validate the hypotheses presented in this work and to further establish DOT imaging as a valuable diagnostic tool