Fast Eikonal Phase Retrieval for High-Throughput Beamlines

Resumo Técnico: Recuperação de Fase Eikonal Rápida para Beamlines de Alto Rendimento

Definição do Problema
A microtomografia baseada em contraste de fase por propagação (PPC-µCT) em fontes de síncrotron de quarta geração permite a imagem de alta sensibilidade de características de fraca refração ao utilizar distâncias de propagação amostra-detector consideráveis. No entanto, este regime introduz desafios significativos para a reconstrução de imagens. Métodos padrão de recuperação de fase de distância única (ex: filtros de Paganin) baseiam-se em suposições de transporte de intensidade linearizadas. Em espécimes heterogêneos contendo tanto tecidos moles fracos quanto interfaces fortes (ex: osso–ar ou osso–fluido), gradientes de fase íngremes geram modulações de intensidade fortes que violam estas suposições lineares. Isso leva a artefatos não lineares característicos, como estrias de longo alcance, que degradam a qualidade da imagem e comprometem a análise subsequente.

Embora formulações completas de óptica de ondas de Fresnel possam modelar estes efeitos, elas exigem a representação explícita pixel a pixel do campo de onda complexo. Isso impõe restrições rigorosas de amostragem; quando as variações de fase são rápidas (como causado por interfaces nítidas), a propagação de Fresnel padrão torna-se computacionalmente proibitiva devido ao aliasing e aos requisitos massivos de memória. Por outro lado, o algoritmo de Recuperação de Fase Eikonal (EPR) introduzido anteriormente modelou com sucesso a refração forte usando uma aproximação eikonal, mas sofreu com altos custos computacionais (dezenas de horas para grandes conjuntos de dados), tornando-se incompatível com as taxas de aquisição de alto rendimento das modernas beamlines.

Metodologia
Os autores introduzem uma formulação de "Recuperação de Fase Eikonal Rápida" (Fast EPR) que acelera o algoritmo EPR original em mais de duas ordens de magnitude, mantendo uma precisão controlada. A metodologia é construída sobre uma expansão assintótica de segunda ordem do propagador de Fresnel na distância de propagação $L$, complementada pela correção de onda de ótica WKB (Wentzel–Kramers–Brillouin) de primeira ordem.

A abordagem emprega um framework unificado com dois solvers complementares para lidar com diferentes regimes físicos:

1. Solver Local (Deslocamentos sub-pixel): Para regimes onde o deslocamento transversal induzido pelo eikonal $\Delta(r)$ é menor que um pixel ($\eta < 1$), o método utiliza um fechamento local $O(L^2)$. Isto envolve uma expansão diferencial de segunda ordem do transporte de intensidade, incluindo o termo eikonal quadrático (WKB0) e o termo de pressão de difração líder (WKB1). Este modelo local é pré-condicionado por operadores inversos dependentes da energia e FFT-diagonais (filtros do tipo Paganin) para acelerar a convergência.
2. Solver Não-local (Deslocamentos multi-pixel): Quando os deslocamentos excedem um pixel ($\eta > 1$), a expansão de Taylor local falha. O método muda para uma formulação não-local onde a intensidade é propagada através de um mapeamento de raio eikonal explícito. Isto é implementado como uma redistribuição bilinear de conservação de massa na grade do detector. Os resíduos do detector são transferidos de volta para a grade do objeto usando o operador adjunto correspondente (transposto), implementado como interpolação bilinear, antes de aplicar o pré-condicionador FFT-diagonal.

Suporte Policromático
O framework suporta dados policromáticos através de uma discretização espectral compacta. O espectro é modelado como uma soma incoerente de $N_s$ componentes espectrais. Todas as quantidades por linha (parâmetros de fase/absorção e deslocamentos) são tratadas de forma dependente da energia, permitendo o transporte e a inversão dependentes da energia enquanto mantém a eficiência de FFT.

Principais Contribuições

• Aceleração Algorítmica: A nova formulação reduz o custo computacional da recuperação de fase em aproximadamente 580 vezes (2,8 ordens de magnitude) em comparação com a implementação EPR original, quando normalizada para os mesmos recursos de GPU.
• Estratégia de Solver Híbrido: A introdução de uma estratégia de dual-solver (solver local $O(L^2)$ para deslocamentos sub-pixel e mapeamento de raio não-local para deslocamentos multi-pixel) garante robustez através de uma ampla gama de condições experimentais, desde tecidos moles de fraca refração até interfaces fortes.
  em modelos de primeira ordem ($L^1$), fornecendo uma renderização mais limpa do contraste de tecido mole próximo a gradientes fortes.
• Física de Segunda Ordem: Ao manter todo o conteúdo $O(L^2)$ da expansão WKB (incluindo o termo de pressão de difração), o método suprime artefatos de estrias não lineares que persistem em modelos lineares de primeira ordem.
• Eficiência Policromática: O método permite a reconstrução precisa de dados policromáticos sem sacrificar a eficiência computacional necessária para fluxos de trabalho de alto rendimento.

Resultos
Os autores validaram o método utilizando um conjunto de dados HiP-CT de cabeça de ovelha e um teste de estresse de amostra de bambu:

• Redução de Artefatos: Comparações no conjunto de dados da cabeça de ovelha demonstraram que o modelo de segunda ordem ($L^2$) reduz significativamente os artefatos de propagação não lineares residuais em comparação com modelos de primeira ordem ($L^1$), proporcionando uma renderização mais limpa do contraste de tecido mole próximo a gradientes fortes.
• Convergência: O algoritmo converge rapidamente; a maior parte da melhoria é alcançada na primeira iteração uma vez que a física de segunda ordem é habilitada.
• Robustez: Na amostra de bambu, onde os deslocamentos de raio excederam o tamanho do pixel do detector, o solver local falhou (divergindo após a primeira iteração), enquanto o solver não-local permaneceu estável e preciso.
• Desempenho: Para um conjunto de dados de 8.000 radiografias ($256 \times 3104$ pixels), o método Fast EPR completou a recuperação de fase em aproximadamente 1,1 a 1,9 minutos em uma única GPU (para dados monocromáticos) e 1,47 a 1,9 minutos para dados policromáticos (5 linhas espectrais). Isto contrasta com o benchmark EPR original, que exigiu aproximadamente 870 minutos para um volume comparável.

Significância
O artigo afirma que esta abordagem unificada permite a recuperação de fase precisa e computacionalmente eficiente através de condições de propagação relevantes para experimentos de PPC-µCT de alto rendimento. Ao reduzir os tempos de reconstrução para serem compatíveis com as taxas de aquisição de beamlines (um exemplo a cada poucas horas), o método facilita o uso rotineiro de recuperação de fase não linear em fluxos de trabalho de alto rendimento. Além disso, a disponibilidade de um modelo de campo próximo de segunda ordem rápido com correção explícita de óptica de ondas permite que pesquisadores operem mais próximos do limite do campo próximo e explorem regimes com gradientes de fase mais fortes (ex: menor energia ou maior resolução) onde abordagens lineares padrão falham, sem incorrer nos custos proibitivos da propagação de Fresnel total. Os autores observam que a complexidade algorítmica do novo método permite que ela seja ocultada atrás da retroprojeção tomográfica em configurações de hardware concorrentes, integrando-a ainda mais em pipelines de reconstrução práticos.