Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RM
Autor(a) principal: | |
---|---|
Data de Publicação: | 2019 |
Tipo de documento: | Tese |
Idioma: | por |
Título da fonte: | Biblioteca Digital de Teses e Dissertações da USP |
Texto Completo: | https://www.teses.usp.br/teses/disponiveis/76/76132/tde-13052020-102037/ |
Resumo: | A Espectroscopia por Ressonância Magnética (ERM) in vivo é uma técnica não invasiva e não ionizante que permite a quantificação da concentração dos metabólitos presentes no corpo humano. No domínio da frequência, as ressonâncias são representadas por linhas espectrais ou picos obtidos através da Transformada de Fourier do sinal adquirido no domínio do tempo. A concentração de cada metabólito é proporcional ao número de núcleos visíveis à técnica de RM e à intensidade de seus respectivos picos. O sinal de um mesmo metabólito pode ser composto por vários picos no domínio da frequência, que podem estar sobrepostos aos demais em função das suas larguras e da resolução espectral. Métodos que realizam ajuste do espectro através da combinação linear de bases espectrais in vitro ou simuladas; e métodos que representam o sinal através de funções matemáticas estão disponíveis na literatura. Pertencente a esta última família, o Krylov Basis Diagonalization Method (KBDM), que representa o sinal através da sobreposição de senoides exponencialmente amortecidas no domínio do tempo foi pouco estudado na literatura no contexto da sua aplicação em ERM in vivo. Portanto, seu posicionamento em relação aos demais métodos da mesma família não é claro. Entre estes, o HLSVD-PRO (HLSVD with partial reorthogonalization), baseado no formalismo de espaço de estados, representa o estado da arte na literatura em tarefas como a remoção do sinal residual da água. Além de realizar um estudo sistemático sobre o KBDM, este trabalho propõe o desenvolvimento de um novo método denominado LLC-KBDM (Line List Clustering Krylov Basis Diagonalization Method), que realiza a promediação das linhas espectrais obtidas através das soluções do KBDM em múltiplos truncamentos do mesmo sinal no espaço de parâmetros. O agrupamento automático dos picos é realizado através do uso de técnicas de aprendizado de máquina não supervisionado. Comparações utilizando sinais de 900 espectros simulados em 9 níveis de ruído foram realizadas entre os métodos LLC-KBDM, KBDM e HLSVD-PRO. O desvio padrão do resíduo medido a partir da diferença entre o sinal simulado sem ruído e o sinal estimado, bem como o número de picos corretamente estimados são utilizados como métricas de comparação. Os resultados mostram que o KBDM é ligeiramente superior ao HLSVD-PRO em todos os níveis de ruído considerados. Por outro lado, o LLC-KBDM se mostrou mais adequado para sinais os mais ruidosos, embora algumas instabilidades numéricas tenham sido observadas no caso contendo o maior nível de ruído. Para sinais in vivo, a tarefa da remoção do sinal residual da água é comparada entre os mesmos métodos. Um teste estatístico que utiliza a curtose e a obliquidade das distribuições para verificação da normalidade dos valores do resíduo na região de supressão foi aplicado. De acordo com o teste para o nível de significância α = 0.05, o sinal residual da água para 10 dos 14 espectros foram adequadamente suprimidos pelo LLC-KBDM; enquanto o KBDM e HLSVD-PRO conseguiram o mesmo feito para apenas 6 e 7 espectros, respectivamente. Uma implementação para o LLC-KBDM está disponível em: https://github.com/danilomendesdias/llckbdm. |
id |
USP_7f04b421da6f7a5dfaa7a7d15419a57d |
---|---|
oai_identifier_str |
oai:teses.usp.br:tde-13052020-102037 |
network_acronym_str |
USP |
network_name_str |
Biblioteca Digital de Teses e Dissertações da USP |
repository_id_str |
2721 |
spelling |
Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RMKrylov basis diagonalization method with spectral line clustering for noisy MRS data fittingAgrupamento de linhas espectraisEspectroscopia por RMHarmonic inversion problemKrylov basis diagonalization methodLLC-KBDMLLC-KBDMMétodo da diagonalização na base de KrylovMR spectroscopyProblema da inversão harmônicaSpectral line list clusteringA Espectroscopia por Ressonância Magnética (ERM) in vivo é uma técnica não invasiva e não ionizante que permite a quantificação da concentração dos metabólitos presentes no corpo humano. No domínio da frequência, as ressonâncias são representadas por linhas espectrais ou picos obtidos através da Transformada de Fourier do sinal adquirido no domínio do tempo. A concentração de cada metabólito é proporcional ao número de núcleos visíveis à técnica de RM e à intensidade de seus respectivos picos. O sinal de um mesmo metabólito pode ser composto por vários picos no domínio da frequência, que podem estar sobrepostos aos demais em função das suas larguras e da resolução espectral. Métodos que realizam ajuste do espectro através da combinação linear de bases espectrais in vitro ou simuladas; e métodos que representam o sinal através de funções matemáticas estão disponíveis na literatura. Pertencente a esta última família, o Krylov Basis Diagonalization Method (KBDM), que representa o sinal através da sobreposição de senoides exponencialmente amortecidas no domínio do tempo foi pouco estudado na literatura no contexto da sua aplicação em ERM in vivo. Portanto, seu posicionamento em relação aos demais métodos da mesma família não é claro. Entre estes, o HLSVD-PRO (HLSVD with partial reorthogonalization), baseado no formalismo de espaço de estados, representa o estado da arte na literatura em tarefas como a remoção do sinal residual da água. Além de realizar um estudo sistemático sobre o KBDM, este trabalho propõe o desenvolvimento de um novo método denominado LLC-KBDM (Line List Clustering Krylov Basis Diagonalization Method), que realiza a promediação das linhas espectrais obtidas através das soluções do KBDM em múltiplos truncamentos do mesmo sinal no espaço de parâmetros. O agrupamento automático dos picos é realizado através do uso de técnicas de aprendizado de máquina não supervisionado. Comparações utilizando sinais de 900 espectros simulados em 9 níveis de ruído foram realizadas entre os métodos LLC-KBDM, KBDM e HLSVD-PRO. O desvio padrão do resíduo medido a partir da diferença entre o sinal simulado sem ruído e o sinal estimado, bem como o número de picos corretamente estimados são utilizados como métricas de comparação. Os resultados mostram que o KBDM é ligeiramente superior ao HLSVD-PRO em todos os níveis de ruído considerados. Por outro lado, o LLC-KBDM se mostrou mais adequado para sinais os mais ruidosos, embora algumas instabilidades numéricas tenham sido observadas no caso contendo o maior nível de ruído. Para sinais in vivo, a tarefa da remoção do sinal residual da água é comparada entre os mesmos métodos. Um teste estatístico que utiliza a curtose e a obliquidade das distribuições para verificação da normalidade dos valores do resíduo na região de supressão foi aplicado. De acordo com o teste para o nível de significância α = 0.05, o sinal residual da água para 10 dos 14 espectros foram adequadamente suprimidos pelo LLC-KBDM; enquanto o KBDM e HLSVD-PRO conseguiram o mesmo feito para apenas 6 e 7 espectros, respectivamente. Uma implementação para o LLC-KBDM está disponível em: https://github.com/danilomendesdias/llckbdm.In vivo Magnetic Resonance Spectroscopy (MRS) is a noninvasive and non-ionizing technique that allows the quantification of human body metabolite concentrations. In the frequency domain, resonances are represented by spectral lines or peaks obtained through the Fourier Transform of the acquired signal in the time domain. The metabolite concentration is proportional to the number of nuclei that are MR-visible and the respective peak intensities. The signal of a single metabolite can be composed by multiple peaks in the frequency domain, that may be overlapped with other peaks due large line-widths and spectral resolution. Methods that are capable of fitting the spectrum through linear combination of simulated or in vitro acquired spectral basis; and methods that represent the signal through model functions are found in the literature. Belonging to this last family of methods, Krylov Basis Diagonalization Method (KBDM) makes the assumption that the signal can be modeled through the superposition of exponentially damped sinusoids in the time domain. In the context of in vivo MRS application, this method has been little studied in the literature. Among these methods, HLSVD-PRO (HLSVD with partial reorthogonalization), which is based on the space-state formalism, represents the state-of-art in the literature for tasks such as residual water removal. Therefore, its position in relation to the other methods of the same family is not clear. Beyond the systematic study of the KBDM, this work also proposes the development of a new technique, named LLC-KBDM (Line List Clustering Krylov Basis Diagonalization Method), that makes use of spectral line list averaging in the parameter space of KBDM solutions that are obtained through multiple signal truncations. The automatic peak clustering is performed through the use of non-supervised machine learning techniques. Comparisons between LLCKBDM, KBDM and HLSVD-PRO were performed using 900 simulated spectra with 9 levels of noise. The standard deviation of the residual computed through the difference between the simulated signal without noise and the estimated signal, as well the number of correctly estimated peaks were used as benchmark metrics. The results show that KBDM is slightly superior to HLSVD-PRO for all noise levels considered. On the other hand, LLC-KBDM has demonstrated to be suitable for noisy signals, although numerical instabilities have been observed in the higher level of noise that was considered. For in vivo signals, the task that constitutes in the removal of residual water signal was chosen to perform the benchmark between the same methods. A statistical test for normality that uses the kurtosis and skewness of the distribution in the suppressed region was used. According to the results at significance level of α = 0.05, the residual water signal was successfully removed from 10 out of 14 spectra by LLC-KBDM; while KBDM and HLSVD-PRO, respectively, only removed 6 and 7 water signal residuals from the same spectra set.. An implementation for LLC-KBDM is available at: https://github.com/danilomendesdias/llckbdm.Biblioteca Digitais de Teses e Dissertações da USPPaiva, Fernando FernandesSilva, Danilo Mendes Dias Delfino da2019-11-27info:eu-repo/semantics/publishedVersioninfo:eu-repo/semantics/doctoralThesisapplication/pdfhttps://www.teses.usp.br/teses/disponiveis/76/76132/tde-13052020-102037/reponame:Biblioteca Digital de Teses e Dissertações da USPinstname:Universidade de São Paulo (USP)instacron:USPLiberar o conteúdo para acesso público.info:eu-repo/semantics/openAccesspor2020-05-18T21:07:01Zoai:teses.usp.br:tde-13052020-102037Biblioteca Digital de Teses e Dissertaçõeshttp://www.teses.usp.br/PUBhttp://www.teses.usp.br/cgi-bin/mtd2br.plvirginia@if.usp.br|| atendimento@aguia.usp.br||virginia@if.usp.bropendoar:27212020-05-18T21:07:01Biblioteca Digital de Teses e Dissertações da USP - Universidade de São Paulo (USP)false |
dc.title.none.fl_str_mv |
Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RM Krylov basis diagonalization method with spectral line clustering for noisy MRS data fitting |
title |
Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RM |
spellingShingle |
Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RM Silva, Danilo Mendes Dias Delfino da Agrupamento de linhas espectrais Espectroscopia por RM Harmonic inversion problem Krylov basis diagonalization method LLC-KBDM LLC-KBDM Método da diagonalização na base de Krylov MR spectroscopy Problema da inversão harmônica Spectral line list clustering |
title_short |
Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RM |
title_full |
Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RM |
title_fullStr |
Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RM |
title_full_unstemmed |
Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RM |
title_sort |
Método da diagonalização na base de Krylov com agrupamento de linhas no ajuste de sinais ruidosos de espectroscopia por RM |
author |
Silva, Danilo Mendes Dias Delfino da |
author_facet |
Silva, Danilo Mendes Dias Delfino da |
author_role |
author |
dc.contributor.none.fl_str_mv |
Paiva, Fernando Fernandes |
dc.contributor.author.fl_str_mv |
Silva, Danilo Mendes Dias Delfino da |
dc.subject.por.fl_str_mv |
Agrupamento de linhas espectrais Espectroscopia por RM Harmonic inversion problem Krylov basis diagonalization method LLC-KBDM LLC-KBDM Método da diagonalização na base de Krylov MR spectroscopy Problema da inversão harmônica Spectral line list clustering |
topic |
Agrupamento de linhas espectrais Espectroscopia por RM Harmonic inversion problem Krylov basis diagonalization method LLC-KBDM LLC-KBDM Método da diagonalização na base de Krylov MR spectroscopy Problema da inversão harmônica Spectral line list clustering |
description |
A Espectroscopia por Ressonância Magnética (ERM) in vivo é uma técnica não invasiva e não ionizante que permite a quantificação da concentração dos metabólitos presentes no corpo humano. No domínio da frequência, as ressonâncias são representadas por linhas espectrais ou picos obtidos através da Transformada de Fourier do sinal adquirido no domínio do tempo. A concentração de cada metabólito é proporcional ao número de núcleos visíveis à técnica de RM e à intensidade de seus respectivos picos. O sinal de um mesmo metabólito pode ser composto por vários picos no domínio da frequência, que podem estar sobrepostos aos demais em função das suas larguras e da resolução espectral. Métodos que realizam ajuste do espectro através da combinação linear de bases espectrais in vitro ou simuladas; e métodos que representam o sinal através de funções matemáticas estão disponíveis na literatura. Pertencente a esta última família, o Krylov Basis Diagonalization Method (KBDM), que representa o sinal através da sobreposição de senoides exponencialmente amortecidas no domínio do tempo foi pouco estudado na literatura no contexto da sua aplicação em ERM in vivo. Portanto, seu posicionamento em relação aos demais métodos da mesma família não é claro. Entre estes, o HLSVD-PRO (HLSVD with partial reorthogonalization), baseado no formalismo de espaço de estados, representa o estado da arte na literatura em tarefas como a remoção do sinal residual da água. Além de realizar um estudo sistemático sobre o KBDM, este trabalho propõe o desenvolvimento de um novo método denominado LLC-KBDM (Line List Clustering Krylov Basis Diagonalization Method), que realiza a promediação das linhas espectrais obtidas através das soluções do KBDM em múltiplos truncamentos do mesmo sinal no espaço de parâmetros. O agrupamento automático dos picos é realizado através do uso de técnicas de aprendizado de máquina não supervisionado. Comparações utilizando sinais de 900 espectros simulados em 9 níveis de ruído foram realizadas entre os métodos LLC-KBDM, KBDM e HLSVD-PRO. O desvio padrão do resíduo medido a partir da diferença entre o sinal simulado sem ruído e o sinal estimado, bem como o número de picos corretamente estimados são utilizados como métricas de comparação. Os resultados mostram que o KBDM é ligeiramente superior ao HLSVD-PRO em todos os níveis de ruído considerados. Por outro lado, o LLC-KBDM se mostrou mais adequado para sinais os mais ruidosos, embora algumas instabilidades numéricas tenham sido observadas no caso contendo o maior nível de ruído. Para sinais in vivo, a tarefa da remoção do sinal residual da água é comparada entre os mesmos métodos. Um teste estatístico que utiliza a curtose e a obliquidade das distribuições para verificação da normalidade dos valores do resíduo na região de supressão foi aplicado. De acordo com o teste para o nível de significância α = 0.05, o sinal residual da água para 10 dos 14 espectros foram adequadamente suprimidos pelo LLC-KBDM; enquanto o KBDM e HLSVD-PRO conseguiram o mesmo feito para apenas 6 e 7 espectros, respectivamente. Uma implementação para o LLC-KBDM está disponível em: https://github.com/danilomendesdias/llckbdm. |
publishDate |
2019 |
dc.date.none.fl_str_mv |
2019-11-27 |
dc.type.status.fl_str_mv |
info:eu-repo/semantics/publishedVersion |
dc.type.driver.fl_str_mv |
info:eu-repo/semantics/doctoralThesis |
format |
doctoralThesis |
status_str |
publishedVersion |
dc.identifier.uri.fl_str_mv |
https://www.teses.usp.br/teses/disponiveis/76/76132/tde-13052020-102037/ |
url |
https://www.teses.usp.br/teses/disponiveis/76/76132/tde-13052020-102037/ |
dc.language.iso.fl_str_mv |
por |
language |
por |
dc.relation.none.fl_str_mv |
|
dc.rights.driver.fl_str_mv |
Liberar o conteúdo para acesso público. info:eu-repo/semantics/openAccess |
rights_invalid_str_mv |
Liberar o conteúdo para acesso público. |
eu_rights_str_mv |
openAccess |
dc.format.none.fl_str_mv |
application/pdf |
dc.coverage.none.fl_str_mv |
|
dc.publisher.none.fl_str_mv |
Biblioteca Digitais de Teses e Dissertações da USP |
publisher.none.fl_str_mv |
Biblioteca Digitais de Teses e Dissertações da USP |
dc.source.none.fl_str_mv |
reponame:Biblioteca Digital de Teses e Dissertações da USP instname:Universidade de São Paulo (USP) instacron:USP |
instname_str |
Universidade de São Paulo (USP) |
instacron_str |
USP |
institution |
USP |
reponame_str |
Biblioteca Digital de Teses e Dissertações da USP |
collection |
Biblioteca Digital de Teses e Dissertações da USP |
repository.name.fl_str_mv |
Biblioteca Digital de Teses e Dissertações da USP - Universidade de São Paulo (USP) |
repository.mail.fl_str_mv |
virginia@if.usp.br|| atendimento@aguia.usp.br||virginia@if.usp.br |
_version_ |
1815256934079528960 |