Reconstrução de imagens de tomografia por emissão de pósitrons com base em compressive sensing e informação a priori

Abstract

A Tomografia por Emissão de Pósitrons(PET), é um exame diagnóstico no âmbito do imageamento médico não invasivo. Ela utiliza, em pequenas quantidades, materiais radioativos denominados radiofármacos, para o diagnóstico de doenças por meio de imagens. Essa técnica combina informações anatômicas e metabólicas, que permitem a detecção de doenças em fases iniciais ao medir a atividade em nível celular no corpo humano. É amplamente usada no tratamento de doenças gastrointestinais, endócrinas, cardíacas e de diversos tipos de câncer. Contudo, o custo de exame PET e o tempo de aquisição dos dados limitam seu uso. Algoritmos de reconstrução alternativos aos tradicionais têm sido desenvolvidos com o propósito de diminuir a quantidade de medidas necessárias ou o tempo em que essas medidas são adquiridas pelo escâner de PET. Uma forma de buscar esse objetivo é utilizar técnicas de reconstrução de dados sub-amostrados baseadas em Compressive Sensing(CS).O CSpermite reconstruir as imagens a partir de uma quantidade de medidas inferior à definida pelo critério de Nyquist. Nessas condições, é fundamental que exista um domínio transformado conhecido em que o sinal seja esparso e que as medidas sejam em números suficientes, o que depende do grau de esparsidade. Além disso, é necessário que o domínio em que o sinal é esparso seja incoerente com relação ao domínio de medidas.Outra abordagem que também permite uma melhora da relação entre qualidade de imagem e o número de medidas é a Pré-filtragem, uma técnica que combinada com CSpermite reduzir o número de coeficientes necessários para a reconstrução. O princípio se baseia na geração de imagens pré-filtradas em uma primeira etapa, de forma a favorecer a esparsidade em cada versão reconstruída. Posteriormente, uma etapa de composição espectral gera a imagem objetivo a partir das versões filtradas. Além disso, a informação de suporte constitui outro método que usado junto aos anteriores também diminui o tempo de aquisição dos dados, com base nas informações de cortes ou quadros prévios. De fato, o uso de Pré-filtragem e o uso de Informação a priori têm melhorado a qualidade de reconstrução de imagens em algumas técnicas de imageamento, sobretudo na Ressonância Magnética. No entanto, não foram encontrados esses algoritmos como referência na literatura para imagens de PET, provavelmente pelas dificuldades adicionais de aplicar um processo de medição computacional e seu operador adjunto que é necessário para implementação de CS.Diante desses desafios, nesse trabalho foram implementados os algoritmos de CScom Pré-filtragem, para reconstrução de sinais com uso de informação a priori, para o caso de medidas de PET. Os algoritmos, uma vez implementados, foram utilizados também sem informação a priori, para efeito de comparação. A implementação exigiu programas de otimização (minimização de l1e de lp), realizados em Octave. Para avaliação desses algoritmos, foram utilizadas imagens de um banco de dados de Tomografia por Emissão de Pósitronsnomeado Laboratóriode Neuroimagem(LONI) Image and DATA Archive(IDA).Os valoresdas medidas de PET foram calculados a partir das imagens obtidas do banco de dados e então reconstruídas empregando os algoritmos clássicos para PET (retroprojeção filtrada), CScom Pré-filtragem sem informação a priori e, por fim, o CScom Pré-filtragem. Foi feita uma comparação da relação sinal erro das imagens reconstruídas a partir de cada método. Também foi realizada uma série de análises da qualidade de reconstrução de imagens à medida em que se aumenta o número de linhas radiais. Os resultados sugerem que as imagens reconstruídas pelo método propostos são associadas a uma melhor qualidade em termos de relação sinal erro (p=3.9874e-63) em comparação com a retroprojeção filtrada. Os testes estatísticos realizados sugerem que as imagens reconstruídas pelo método proposto (CScom Pré-filtragem e informações a priori) resultam em uma melhor qualidade média da imagem em comparação aos outros dois métodos. Foram realizados testes preliminares com informação a priori com sinais de domínio unidimensional que mostraram que esses resultados melhoram em termos de SNR (dB) a qualidade da imagem reconstruída. Além disso, foi investigado o impacto que essa informação a priori pode ter nas imagens PET. Foram consideradas sequências de imagens que representam diferentes quadros temporais (ouframes). A informação que pode ser extraída de um quadro é utilizada para auxiliar na reconstrução do próximo quadro. Ao final, foi avaliado o impacto que isso tem em termos da qualidade da reconstrução conforme medida pela SNR (dB), apresentando os melhores resultados em imagens de dimensões maiores.