An angular distortion matrix approach for joint wave-speed tomography and aberration correction in scattering media

Resumo Técnico: Abordagem de Matriz de Distorção Angular para Tomografia Conjunta de Velocidade de Onda e Correção de Aberração

Definição do Problema
Na imagem de reflexão (ultrassom, óptica, sismologia), a reconstrução precisa da refletividade do meio depende de um modelo de propagação de onda rigoroso, governado pela distribuição da velocidade da onda. No entanto, em meios de espalhamento complexos, a velocidade da onda raramente é conhecida a priori. Desajustes entre a velocidade assumida e a real introduzem aberrações de fase que degradam o foco, alargando a função de espalhamento de ponto (PSF) e distorcendo a localização aparente dos espalhadores. Embora essas aberrações codifiquem informações sobre as heterogeneidades subjacentes da velocidade da onda, separar essas informações da refletividade do meio é difícil, particularmente em meios dominados por speckle. Métodos existentes enfrentam um problema fundamental de acoplamento: corrigir aberrações requer um modelo de velocidade de onda preciso, mas recuperar esse modelo requer observáveis que sejam, elas próprias, não enviesadas pelas aberrações. Além disso, abordagens tradicionais de matriz de distorção assumem que as aberrações surgem ou do caminho incidente ou do caminho refletido, enquanto na realidade ambas as contribuições estão inerentemente misturadas, levando a deslocamentos de profundidade que não podem ser corrigidos sem atualizar o modelo de velocidade.

Metodologia
Os autores propõem um framework unificado que estima iterativamente a distribuição da velocidade da onda e corrige aberrações usando uma nova matriz de distorção angular ($D_{\theta\theta}$). A metodologia procede da seguinte forma:

1. Imagem de Matriz e Representação Frequência-Número de Onda:
   O processo começa com a aquisição de uma matriz de reflexão $R_{uu}$ (respostas transdutor-para-transdutor). Esta é projetada em uma base focada para formar uma matriz de reflexão focada $R_{xx}$. Para isolar as distorções, os autores introduzem uma representação frequência-número de onda das distorções de onda, $D_{kk}$. Este operador é construído removendo os componentes de fase geométrica esperados das respostas, como se o espalhamento originasse de um ponto arbitrário $x$. Isso permite a síntese de distorções para qualquer experimento virtual em pós-processamento.

2. Construção da Matriz de Distorção Angular:
   O operador $D_{kk}$ é projetado em uma base de feixe Gaussiano para formar a matriz de distorção angular $D_{\theta\theta}$. Esta base utiliza feixes Gaussianos direcionados, definidos pelos ângulos incidente ($\theta_{in}$) e refletido ($\theta_{out}$). Esta projeção aproveita o efeito de memória angular, onde as ondas exibem fortes correlações ao longo de direções angulares específicas, apesar do desordem. Ao analisar as antidiagonais de $D_{\theta\theta}$ (ângulo médio constante $\theta_m = (\theta_{in} + \theta_{out})/2$), o método isola aberrações de fase associadas a caminhos de propagação bem definidos.

3. Estimativa de Aberração Local:
   O campo de visão (FOV) é dividido em "patches isoplanáticos" sobrepostos. Dentro de cada patch, uma matriz de correlação local é computada a partir dos elementos da matriz de distorção angular. Sob a suposição de que as distorções de onda são invariantes sobre pequenas regiões, o primeiro autovetor desta matriz de correlação fornece uma estimativa da lei de aberração de fase local ($u_m$). Para evitar o aliasing de fase, o método utiliza diferenças de fase relativas entre sucessivos ângulos de meia-diferença ($\theta_d$).

4. Inversão Iterativa para Tomografia de Velocidade de Onda:
   As aberrações de fase estimadas são usadas para reconstruir a distribuição da velocidade do som. O problema é formulado como uma tarefa de otimização para minimizar o desajuste entre as diferenças de fase medidas e as previstas por um modelo de propagação. A inversão é realizada iterativamente:

   • Assume-se um modelo homogêneo inicial (ex: $c_0 = 1540$ m/s).
   • As aberrações são estimadas e corrigidas.
   • Um problema inverso linearizado (usando uma aproximação de raio reto para o Jacobiano e regularização de Tikhonov) resolve as perturbações de slowness ($\Delta s$).
   • O modelo de propagação de onda é atualizado, e o processo se repete até que as aberrações residuais se tornem negligenciáveis.

5. Tratamento de Espalhamento Múltiplo:
   Em meios altamente heterogêneos (ex: tecido mamário) onde o espalhamento múltiplo degrada o efeito de memória, um filtro confocal adaptativo é aplicado à matriz de reflexão antes de computar a matriz de distorção. Isso suprime o ruído de fundo incoerente enquanto preserva a energia diagonal associada ao espalhamento simples.

Principais Resultados
A abordagem foi validada usando experimentos de ultrassom em três cenários:

• Fantoma Simulador de Tecido: Um fantoma com estruturas de camadas conhecidas (gordura, músculo, fígado) e velocidades de som variando (1420–1550 m/s) foi utilizado. O método iterativo recuperou com sucesso o mapa de velocidade do som com erros médios dentro de 1% do valor real (ground truth). A razão de Strehl (uma medida da qualidade do foco) melhorou de 0,34 para 0,96, e as imagens de refletividade mostraram interfaces nítidas e contraste de speckle aprimorado.
• Imagem Mamária In Vivo: Aplicado a uma paciente com carcinoma invasivo. O método identificou uma região localizada de velocidade de som elevada (~1574 m/s), consistente com lesões malignas, distinguindo-a do tecido adiposo e fibroglandular circundante. A correção de aberração revelou que as imagens não corrigidas haviam enviesado a profundidade aparente de características anatômicas.
• Imagem de Fígado In Vivo: Aplicado a um voluntário saudável. O método resolveu a estrutura em camadas da parede abdominal (pele, gordura, músculo) e o parênquima hepático com alta resolução axial. As velocidades de som recuperadas coincidiram com valores da literatura para tecidos saudáveis, e a razão de Strehl aumentou de 0,61 para 0,93, nitidamente aguçando a imagem de refletividade e revelando estruturas anteriormente obscurecidas, como veias anecoicas.

Significância e Alegações
O artigo afirma que a abordagem da matriz de distorção angular fornece uma solução robusta para o problema acoplado de tomografia de velocidade de onda e correção de aberração. Sua principal significância reside em:

• Unificação: Integra conceitos de matriz de distorção (tipicamente usados para correção de aberração) com a tomografia de ultrassom computada em modo eco (CUTE) dentro de um único framework de imagem de matriz iterativo.
• Independência de Modelo: A formulação da matriz de distorção angular é independente do modelo específico de propagação de onda usado para o foco, permitando que o método refine progressivamente o próprio modelo.
• Alta Resolução: Ao projetar em uma base de feixe Gaussiano, o método pode estimar aberrações sobre patches isoplanáticos menores (aprox. 1 mm) em comparação com as abordagens clássicas de matriz de distorção, permitindo maior resolução espacial na tomografia.
• Utilidade Clínica: O método demonstra a capacidade de fornecer mapas quantitativos de velocidade do som (potenciais biomarcadores para patologia) e, simultaneamente, melhorar a qualidade da imagem estrutural para interpretação clínica em tecidos complexos e heterogêneos.
• Generalidade: Embora demonstrada com ultrassom, o framework é aplicável a qualquer modalidade de onda (óptica, sismologia) onde a imagem de matriz de reflexão possa ser implementada.