Inversion of limited-aperture Fresnel experimental data using orthogonality sampling method with single and multiple sources

Probleemstelling

Het artikel adresseert het probleem van de inverse verstrooiing (inverse scattering) met een beperkte apertuur (limited-aperture). Het doel is het snel en stabiel identificeren van de aanwezigheid, locatie en omtrek van kleine objecten (diëlektrische cilinders) op basis van gemeten verstrooide velddata.

• Context: De meeste bestaande algoritmen zijn iteratief (zoals Newton-type methoden) en vereisen een goede initiële schatting, wat computatie-intensief is. Niet-iteratieve methoden, zoals de Orthogonal Sampling Method (OSM), zijn sneller en stabieler.
• Uitdaging: Hoewel OSM effectief is voor synthetische data, zijn er beperkingen bij het toepassen op experimentele data (van het Fresnel-instituut). Er ontbreekt een adequate wiskundige theorie om te verklaren waarom de beeldkwaliteit sterk afhankelijk is van de positie van de bron (emitter) en de werkfrequentie. Bovendien is het moeilijk om meerdere objecten te onderscheiden bij lage of zeer hoge frequenties met een enkele bron.

Methodologie

De auteur analyseert en verbetert de OSM door een wiskundige relatie te leggen tussen de indicatorfunctie en Besselfuncties (Bessel functions).

1. OSM met enkele bron (Single Source):

   • De auteur leidt een asymptotische expansie af voor de indicatorfunctie $F_{OSM}(r', a_m)$ bij gebruik van één emitter $a_m$.
   • Wiskundige structuur: De indicatorfunctie wordt uitgedrukt als een integraal die de Besselfunctie van de eerste soort van orde nul ($J_0$) en een oneindige reeks van Besselfuncties van niet-nul gehele orde bevat.
   • Analyse: De term $J_0(k_b|r' - r|)$ zorgt voor de pieken op de locatie van de objecten, maar de term $E(r', r, a_m)$ (een oneindige reeks) veroorzaakt artefacten. De amplitude hangt af van de Green-functie $G(r, a_m)$, wat betekent dat de beeldkwaliteit sterk varieert met de positie van de emitter en de frequentie.

2. OSM met meerdere bronnen (Multiple Sources - MOSM):

   • Om de beperkingen van de enkele bron te overwinnen, wordt een nieuwe indicatorfunctie ontworpen die gebruikmaakt van data van meerdere emitters ($M$ bronnen).
   • Nieuwe indicatorfunctie: In plaats van een simpele som van indicatorfuncties, wordt een nieuwe vectoriële benadering gebruikt die de orthogonaliteit tussen de data en de testvectoren maximaliseert.
   • Wiskundige structuur: De nieuwe indicatorfunctie $F_{MOSM}(r')$ wordt uitgedrukt als een integraal die de kwadraten van de Besselfuncties bevat ($J_0^2$ en een reeks van $J_p^2$).
   • Theoretisch voordeel: Het kwadrateren van de Besselfuncties vermindert de amplitude van de artefacten veroorzaakt door de oneindige reekstermen, wat leidt tot een scherpere en betrouwbaardere detectie.

Belangrijkste Bijdragen

• Wiskundige Karakterisering: Voor het eerst wordt een expliciete wiskundige relatie afgeleid tussen de OSM-indicatorfunctie en oneindige reeksen van Besselfuncties voor zowel enkele als meerdere bronnen.
• Verklaring van Fenomenen: De theorie verklaart waarom:
  • De detectie bij een enkele bron sterk afhankelijk is van de emitterpositie.
  • Zeer hoge frequenties leiden tot verwaarloosbare signalen (door de afname van de Green-functie).
  • Lage frequenties leiden tot onvoldoende resolutie om objecten te onderscheiden (afhankelijk van de golflengte $\lambda$).
• Ontwerp van een Verbeterde Algoritme: De auteur introduceert een nieuwe indicatorfunctie voor meerdere bronnen die wiskundig bewezen is om unieke bepaling (unique determination) van kleine objecten mogelijk te maken, onafhankelijk van de emitterpositie.
• Validatie met Experimentele Data: De theorie wordt getoetst aan echte meetdata van het Fresnel-instituut, in plaats van alleen aan synthetische data.

Resultaten

De resultaten zijn gebaseerd op simulaties met experimentele data van twee kleine diëlektrische cilinders:

• Enkele Bron (OSM):
  • De prestaties zijn zeer variabel. Bij sommige emitterposities en frequenties (bijv. 2 GHz) zijn objecten zichtbaar, maar bij andere posities of frequenties (bijv. 1 GHz of $\ge$ 6 GHz) falen de detectie of ontstaan er sterke artefacten.
  • De vorm van de objecten kan vaak niet nauwkeurig worden bepaald.
• Meerdere Bronnen (Nieuwe MOSM):
  • De nieuwe indicatorfunctie $F_{MOSM}$ levert consistente resultaten op, ongeacht de positie van de emitters.
  • Bij frequenties van 2, 3 en 4 GHz worden zowel de locatie als de omtrek van de objecten succesvol en duidelijk gereconstrueerd.
  • De Jaccard-index (een maat voor overlap tussen de gereconstrueerde en werkelijke vorm) toont aan dat de nieuwe methode iets beter presteert dan de traditionele som-methode voor meerdere bronnen.
  • Bij zeer hoge frequenties (6-8 GHz) is detectie mogelijk, maar de beeldkwaliteit neemt af (vergelijkbaar met de enkele bron), wat bevestigt dat er een optimale frequentieband is.

Betekenis en Conclusie

Dit onderzoek biedt een fundamentele theoretische basis voor het gebruik van de Orthogonal Sampling Method in praktische, experimentele situaties met beperkte apertuur.

• Het bewijst dat de beperkingen van de OSM met een enkele bron (afhankelijkheid van bronpositie en frequentie) wiskundig kunnen worden verklaard en opgelost door het gebruik van meerdere bronnen.
• De nieuwe indicatorfunctie met meerdere bronnen biedt een robuustere, snellere en nauwkeurigere oplossing voor het identificeren van kleine objecten, wat relevant is voor toepassingen zoals medische beeldvorming, defectdetectie in beton en landmijnopsporing.
• De auteur stelt dat uitbreiding naar 3D Fresnel-gegevensbanken een logische volgende stap is voor toekomstig onderzoek.