[go: up one dir, main page]

DE102006009121B4 - Verfahren zur Verarbeitung und Darstellung von mittels Synthetik-Apertur-Radarsystemen (SAR) gewonnen Bodenbildern - Google Patents

Verfahren zur Verarbeitung und Darstellung von mittels Synthetik-Apertur-Radarsystemen (SAR) gewonnen Bodenbildern Download PDF

Info

Publication number
DE102006009121B4
DE102006009121B4 DE102006009121A DE102006009121A DE102006009121B4 DE 102006009121 B4 DE102006009121 B4 DE 102006009121B4 DE 102006009121 A DE102006009121 A DE 102006009121A DE 102006009121 A DE102006009121 A DE 102006009121A DE 102006009121 B4 DE102006009121 B4 DE 102006009121B4
Authority
DE
Germany
Prior art keywords
sar
image
azimuth
speed
focusing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
DE102006009121A
Other languages
English (en)
Other versions
DE102006009121A1 (de
Inventor
Dr. Eineder Michael
Prof. Dr. Bamler Richard
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Deutsches Zentrum fuer Luft und Raumfahrt eV
Original Assignee
Deutsches Zentrum fuer Luft und Raumfahrt eV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Deutsches Zentrum fuer Luft und Raumfahrt eV filed Critical Deutsches Zentrum fuer Luft und Raumfahrt eV
Priority to DE102006009121A priority Critical patent/DE102006009121B4/de
Priority to IT000006A priority patent/ITSV20070006A1/it
Priority to IT000007A priority patent/ITSV20070007A1/it
Priority to CA2579898A priority patent/CA2579898C/en
Publication of DE102006009121A1 publication Critical patent/DE102006009121A1/de
Application granted granted Critical
Publication of DE102006009121B4 publication Critical patent/DE102006009121B4/de
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/91Radar or analogous systems specially adapted for specific applications for traffic control
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9023SAR image post-processing techniques combined with interferometric techniques
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/89Radar or analogous systems specially adapted for specific applications for mapping or imaging
    • G01S13/90Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
    • G01S13/9021SAR image post-processing techniques
    • G01S13/9029SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S13/00Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
    • G01S13/88Radar or analogous systems specially adapted for specific applications
    • G01S13/91Radar or analogous systems specially adapted for specific applications for traffic control
    • G01S13/92Radar or analogous systems specially adapted for specific applications for traffic control for velocity measurement

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

Verfahren zur Verarbeitung und zur Darstellung von mittels flugzeug- oder weltraumgestützten Synthetik-Apertur-Radarsensor-Systemen (SAR) gewonnenen Bodenbildern unter Nutzung eines bekannten Verkehrswegeverlaufs als Vorwissen, wobei eine hypothetische Fahrzeug-Radialgeschwindigkeit vr jedes als Fahrzeug in Betracht kommenden Bildpunktes aus seinem Abstand zu einem Verkehrsweg in Azimutrichtung und die zugehörigen Fahrzeugbildposition durch Azimutprojektion auf den betreffenden Verkehrsweg ermittelt wird und eine hypothetische Fahrzeug-Azimutgeschwindigkeit va aus der hypothetischen Fahrzeug-Radialgeschwindigkeit vr, dem Einfallswinkel Θ des Radarstrahls auf die Erdoberfläche und dem zwischen dem betreffenden Verkehrsweg und der Flugbahn des Radarsensors vorliegenden Winkel α mittels der Funktionabgeleitet wird, und dass der durch geographische Koordinaten und Höhen gegebene Verkehrswegeverlauf einer Karte in Azimut/Range-Koordinaten des SAR-Radarbildes konvertiert wird, dass aus der ermittelten hypothetischen Fahrzeug-Radialgeschwindigkeit vr in Abhängigkeit der Radarparameter ein erster, die Verschiebung des Azimut-Spektrums betreffender Fokussierungsparameter ΔfDC abgeleitet wird, dass aus der hypothetischen Fahrzeug-Azimutgeschwindigkeit va unter Zuhilfenahme der Flugbahn des Radarsensors und der Radarparameter ein zweiter, die Änderung der...

Description

  • Die Erfindung betrifft ein Verfahren zur Verarbeitung und zur Darstellung von mittels flugzeug- oder weltraumgestützten Synthetik-Apertur-Radar-Systemen gewonnenen Bodenbildern.
  • Synthetik-Apertur-Radar-(SAR-)Systeme werden zur Fernerkundung aus dem Weltraum oder aus dem Flugzeug eingesetzt. Dabei wird im zivilen Bereich meistens die ruhende Erdoberfläche abgebildet. Zukünftige weltraumgestützte Radarsysteme haben eine genügend hohe Auflösung, um auch einzelne sich bewegende Objekte abbilden zu können. Falls zwei oder mehrere in Flugrichtung ortlich getrennte Antennen verwendet werden, wird des Objekt zu leicht unterschiedlichen Zeiten beobachtet und es erschließt sich die Möglichkeit, ein sich bewegendes Objekt zu erkennen und seine Geschwindigkeit zu erfassen. Diese kann aus der interferometrischen Phasendifferenz
    Figure 00010001
    der aus den beiden Antennen gewonnenen SAR-Bildern abgeleitet werden. Dabei sind:
  • λ
    die Wellenlange des Radars,
    BATI
    der Abstand der beiden Antennen in Flugrichtung,
    VSAR
    die Geschwindigkeit des SAR-Systems, und
    vr
    die Geschwindigkeit des Objektes in radialer Richtung zur Flugbahn des Radars.
  • Wegen der Beschränkung der beobachteten Phase auf des Intervall [–180°, 180°] kann die Geschwindigkeit allerdings nur mehrdeutig gemessen werden. Beide Effekte, die Phase und die später noch beschriebene geschwindigkeitsabhängige Verschiebung, ergeben zusammen gute Hinweise auf die Geschwindigkeit des Objektes.
  • Allerdings können die von sich bewegenden Objekten reflektierten Signale vergleichsweise schwach sein und deshalb kann deren Erkennung und Messung vor dem Hintergrund (Clutter) entsprechend schwierig sein.
  • Bei einen Synthetik-Apertur-Radar wird die Bewegung der Antenne entlang einer bekannten Flugbahn genutzt, um eine großere Apertur und damit eine hohere Auflösung zu erreichen. Bewegt sich nun das abzubildende Objekt während der Beleuchtungszeit durch die Antenne, so hat dies, wie in 1 schematisch gezeigt ist, mehrere Störeinflusse zur Folge:
    • a) Eine Bewegungskomponente Vr des Objektes in Richtung der Verbindungslinie zwischen Sensor und Objekt erzeugt im SAR-Bild eine Verschiebung Δa des Objektes in Flugrichtung. Ein sich auf einer Straße bewegendes Fahrzeug wird deshalb abseits der Straße abgebildet. Da das Ruckstreusignal der dort vorherrschenden Gebäude, Walder und Felder in der Größenordnung des Signals des Fahrzeugs liegen kann, ist dieses nur schwer zu finden. Die Verschiebung kann auch durch ein verschobenes Dopplerspektrum des bewegten Objektes in Bezug auf das Spektrum der ruhenden Objekte erklart werden. Wegen der fehlenden Bandbegrenzung und der Abtastung des Dopplerspektrums mit der Radar-Pulswiederholfrequenz werden verschiedene Geschwindigkeiten mehrdeutig auf dem gleichen Ort abgebildet.
    • b) Wegen des verschobenen Dopplerspektrums kann es bei einer fur ruhende Objekte eingestellten SAR-Fokussierung zudem leicht zu fehlerhaften ”Geisterbildern” (Ambiguities) kommen. Auch dieser Effekt bewirkt eine Signalabschwachung, da sich die empfangene Energie des Fahrzeugs auf zwei weit auseinanderliegende Bildpunkte verteilt. Entfernt sich das Fahrzeug während der Antennenbeleuchtungszeit um mehr als eine Auflosungszelle, so wird wegen der zusätzlichen ”Range Migration” die Azimut- und Range-Fokussierung verschlechtert und das Objekt erscheint noch dunkler.
    • c) Eine Bewegungskomponente va des Objektes in Flugrichtung der Antenne verändert die für die SAR-Fokussierung wichtige Azimut-Frequenzmodulationsrate FM und bewirkt somit eine Verschmierung der Punktantwort in Flugrichtung.
  • 1 zeigt diesbezüglich schematisch die Fahrzeug-Fehlabbildungen in einem gewöhnlichen herkömmlichen SAR-Bild. Danach wird ein bewegtes Fahrzeug bei der Abbildung durch ein Synthetik-Apertur-Radar gegenüber dem ruhenden Hintergrund in Flugrichtung versetzt, im Azimut verschmiert und moglicherweise mehrfach abgebildet (Ambiguity).
  • Bei Radarsystemen, die fur die Abbildung ruhender Objekte konstruiert sind, ist aus obigen Gründen die Verarbeitung der Daten und die Erkennung von Fahrzeugen mit erheblichen Schwierigkeiten verbunden.
  • Deshalb und wegen der relativ geringen Auflosung von heutigen zivilen Radarsystemen spielt die Detektion von bewegten Fahrzeugen, die auch GMTI (Ground Moving Target Indication) genannt wird, bisher nur im militärischen Bereich eine Rolle. Dort werden fur GMTI optimierte, flugzeuggetragene Radarsysteme eingesetzt. Um die Entdeckung von bewegten Objekten zu erleichtern, haben diese Systeme beispielsweise mehrere Antennen, eine stark erhöhte Pulswiederholfrequenz und eine starkere Sendeleistung.
  • Um das im Vergleich zum Hintergrund schwache Signal eines Fahrzeugs in einem Radarbild zu finden, muss das Fahrzeug optimal fokussiert sein, d. h. seine gesamte Energie sollte auf einen Bildpunkt konzentriert sein. Wegen der zuvor genannten Effekte mussten dazu allerdings Ort, Geschwindigkeit und Richtung des Fahrzeugs bereits vor der Verarbeitung bekannt sein.
  • Eine verhaltnismäßig naheliegende, aber sehr aufwandige Losung fur dieses Problem bestunde darin, in einem ersten Schritt jeweils ein Bild für alle möglichen Kombinationen aus Radial- und Azimutgeschwindigkeiten (vr und va) zu erzeugen und in dieser Bildermenge in einem zweiten Schritt die Fahrzeuge zu suchen. Allerdings kann die Zahl der dazu notwendigen Bilder durchaus bei mehreren hundert liegen. Deren Erzeugung ist sehr rechenaufwändig und auch die numerische Weiterverarbeitung, z. B. die Fahrzeugdetektion in diesem Bildstapel, ist kompliziert. Auch eine visuelle Auswertung einer so großen Bildermenge ist unrealistisch.
  • Die Grundlagen der Verschiebung sowie Verschmierung in der Darstellung von Bewegtzielen bei der SAR-Verarbeitung sind aus dem Artikel von Hinz, S. u. a. (Spaceborn Traffic Monitoring with Dual Channel Synthetic Aperture Radar-Theory and Experiments, In: Proceedings of the 2005 IEEE Computer Society Conference an Computer Vison and Pattern Recognition CVPR'05, 2005) bekannt. Bei dem bekannten Verfahren wird ein bekannter Verkehrswegeablauf als a priori Wissen genutzt und die Radialgeschwindigkeit vr eines Bildpunktes mit einem Abstand zu einem Verkehrsweg bestimmt. Weiter kann bei dem bekannten Verfahren durch einfache Trigonometrie die Azimutgeschwindigkeit va bestimmt werden, wenn der Bildpunkt dem entsprechenden Verkehrsweg zugeordnet ist.
  • Ferner ist aus dem zuvor genannten Artikel bekannt, dass eine Azimutfokussierung von SAR Bildern mittels eines signalangepassten Filters erfolgt. Um das signalangepasste Filter zu berechnen, muss die tatsächliche Entfernungshistorie eines jeden Ziels simuliert werden. Hierzu ist ein a priori Wissen über Sensor, Streupunkt und Bewegung erforderlich.
  • Aus den Artikeln von Thompson, A. A. and Livingstone, C. E. (Moving Target Performance for RADARSAT 2, in: Geoscience and Remote Sensing Symposium, 2000, Proceedings, IGARSS 2000, IEEE 2000 International. Volume 6, 24–28 July 2000, S. 2599–2601) und Suchandt, S. u. a. (Results from a Airborne SAR GMTI Experiment supporting TerraSAR-X Traffic Processor Development, in: Geoscience and Remote Sensing Symposium, 2005, IGARSS '05, Proceedings, 2005 IEEE International, Volume 4, 25–29 July, S. 2949–2952) ist es dem Fachmann ferner bekannt, dass eine Geschwindigkeitskomponente des Streupunktes in radialer Richtung zu einer Verschiebung der Doppelfrequenz und eine Bewegung des Streupunktes in Azimutrichtung zu einer Änderung der Erequenzmodulationsrate führt.
  • Aus DE 44 23 899 C1 ist ein Verfahren zur Detektion, Lokalisierung und Geschwindigkeitsbestimmung von Bewegtzielen aus Radar-Rohdaten eines von einem Träger mitgeführten, kohärenten, ein- oder mehrkanaligen SAR Abbildungssystemen bekannt. Hierbei wird mit einem Abbildungsradar die Landschaft einer Straße aufgenommen und anhand einer Dopplerate-Karte bzw. einer Frequenzversatz-Karte ein Bewegtzielen detektiert und lokalisiert.
  • In dem Artikel von Sharma, J. (The influence of Target Acceleration on Dual-Channel SAR-GMTI (Synthetic Aperture Radar Ground Moving Target Indication) Data, UCGE Reports number 20206, October 2004, University of Calgary, S. 186–189) wird im Zusammenhang mit der Auswertung von SAR-Bildern von sich bewegenden Objekten auf die Problematik von Beschleunigungseffekten hingewiesen und deren Bestimmung aus den gewonnenen Daten mit Hilfe einer Filterbank, d. h. durch Ausprobieren, in Aussicht gestellt.
  • Zusammenfassend ist demnach festzustellen, dass die Erkennung solcher sich bewegender Objekte in SAR-Bildern sehr schwierig ist, weil sie gegenüber ruhenden Objekten verschoben erscheinen und weil sie aufgrund ihrer Eigenbewegung unscharf fokussiert werden.
  • Aufgabe der vorliegenden Erfindung ist es daher, bewegte Fahrzeuge, insbesondere Straßen und Schienenfahrzeuge mit Hilfe von flugzeug- oder weltraumgestützten Synthetik-Apertur-Radar-Systemen (SAR) besser und dabei mit weniger technischem und personellem Aufwand zu erfassen.
  • Gemäß der Erfindung, wird diese Aufgabe gelöst durch ein Verfahren zur Verarbeitung und zur Darstellung von mittels flugzeug- oder weltraum-gestützten Synthetik-Apertur-Radarsensor-Systemen (SAR) gewonnenen Bodenbildern mit den Merkmalen des Anspruchs 1.
  • Die Erfindung nutzt das Vorwissen des Verkehrsweges für eine prazise und effiziente Fokussierung sowie eine Darstellungsform, welche die Interpretation der SAR-Radaraufnahmen wesentlich erleichtert.
  • Die Erfindung zeichnet sich dadurch aus, dass digitale Verkehrswegekarten, also vor allem digitale Straßenkarten, Schienenkarten oder auch Flüsseverlaufkarten, in die SAR-Fokussierung einbezogen werden. Der Verkehrswegeverlauf wird dabei in die Radar-Aufnahmegeometrie transformiert. Es erfolgt dann eine Bestimmung der Fahrzeug-Radialgeschwindigkeit jedes Bildpunktes aus seinem Abstand zum Verkehrsweg in Azimutrichtung und der zugehörigen Fahrzeugposition durch Azimut-Projektion auf die Straße. Im Anschluss daran wird die Azimutgeschwindigkeit aus obiger Radialgeschwindigkeit und dem Straßenwinkel ermittelt. Schließlich erfolgt aus ermittelter Radial- und Azimutgeschwindigkeit die präzise und effiziente Fokussierung des einzelnen Bildpunktes, der einem bewegten Objekt, also einem sich bewegenden Fahrzeug, entspricht.
  • Die Vorteile des gemäß der vorliegenden Erfindung arbeitenden Verfahrens bestehen vor allem darin, dass nur ein SAR-Bild erzeugt und nach Fahrzeugen durchsucht werden muss, dass also eine erhebliche Zeitersparnis und Vereinfachung stattfinden. Jedes Fahrzeug wird optimal fokussiert, was die Fahrzeug-Detektionswahrscheinlichkeit steigert. Der Hintergrund (Clutter) für schnell bewegte Objekte wird kontrastärmer und für manche Geschwindigkeitsbereiche dunkler, was ebenfalls die Detektionswahrscheinlichkeit steigert. Schließlich kann bei Anwendung des Verfahrens nach der vorliegenden Erfindung das Bild leicht und schnell von einem menschlichen Auswerter interpretiert werden.
  • Vorteilhafte und zweckmäßige Weiterbildungen der Erfindung sind in den sich auf den Anspruch 1 unmittelbar oder mittelbar rückbeziehenden Ansprüchen angegeben.
  • Die Fokussierung des einzelnen Bildpunktes oder seiner Umgebung unter Berücksichtigung der Azimut- und Radialgeschwindigkeit kann vorteilhaft im Zeitbereich mit einem adaptiven Zeitbereichskorrelator erfolgen.
  • Alternativ kann die Fokussierung mehrerer Bilder mit verschiedenen Geschwindigkeitsparametern (vr, va) durch Standard-SAR-Prozessoren auf Vorrat vorgenommen werden, wobei nur Bildbereiche in der Nahe von Straßen fokussiert werden. Dann erfolgt eine Kombination zu einem Bild, wobei jeder Bildpunkt aus dem Bild entnommen wird, das der hypothetischen Radial- und Azimutgeschwindigkeit entspricht.
  • In einer anderen Alternative kann eine Vorfokussierung des Bildes mit einem modifiziertem Standard-SAR-Prozessor und eine Nachfokussierung mit ortsadaptivem Korrekturfilter im Zeitbereich durchgeführt werden. Die Modifikation eines Standard-SAR-Prozessors besteht in der internen Replikation des Azimut-Spektrums.
  • Eine Weiterbildung der Erfindung besteht in der Darstellung nur des Straßenverlaufes, also nicht des ganzen Bildes, mit Geschwindigkeitskonturen zur visuellen menschlichen Auswertung.
  • Eine farbige Darstellung des interferometrischen SAR-Bildes mit überlagerten farbigen Geschwindigkeitskonturen für Mehkanalsysteme (ATI; Along-Track-Interferometrie) ist in vorteilhafter Weise möglich.
  • Auch eine farbige Darstellung des SAR-Bildes mit überlagerten farbigen Geschwindigkeitskonturen für DPCA(Displaced Phase Center Antenna)-Systeme, bei denen es sich um Systeme mit einer Antenne mit zwei oder mehr versetzten Phasenzentren handelt, lässt sich in vorteilhafter Weise beim Verfahren nach der vorliegenden Erfindung realisieren.
  • Eine Verwendung des Verfahrens nach der vorliegenden Erfindung besteht in besonders vorteilhafter Weise bei der Erfassung und dem Einsatz sowie der Vermarktung von Verkehrsdaten für Wissenschaft, Wirtschaft und Sicherheit und bei der Erstellung von übersichtlichen Bildern von Verkehrsflüssen.
  • Die Erfindung wird nachfolgend anhand von Zeichnungen und Diagrammen erläutert. Es zeigen:
  • 1 die bereits einleitend beschriebene schematische Darstellung zur Erläuterung der bei üblichen herkömmlichen SAR-Systemen bei der Abbildung von bewegten Fahrzeugen auftretenden Schwierigkeiten, die darin bestehen, dass die die Fahrzeuge darstellenden Bildpunkte gegenuber dem ruhenden Hintergrund in Flugrichtung versetzt, im Azimut verschmiert und möglicherweise mehrfach (Ambiguity) abgebildet werden,
  • 2 eine schematische Darstellung zur Erläuterung der kartengesteuerten geschwindigkeitsabhängigen SAR-Fokussierung gemäß dem Verfahren nach der vorliegenden Erfindung,
  • 3 die schematische Darstellung eines entsprechend dem Verfahren nach der vorliegenden Erfindung geschwindigkeitsabhängig fokussierten SAR-Bildes mit Iso-Geschwindigkeiten,
  • 4 das Ablaufdiagramm einer entsprechend dem Verfahren nach der vorliegenden Erfindung ablaufenden geschwindigkeitsadaptiven SAR-Prozessierung mit nachfolgender SAR-Verarbeitung in einer ersten Variante (Variante a) unter Verwendung einer Visualisierung mit einem adaptiven Zeitbereichskorrelator,
  • 5 das Ablaufdiagramm einer entsprechend dem Verfahren nach der vorliegenden Erfindung ablaufenden geschwindigkeitsadaptiven SAR-Prozessierung mit nachfolgender SAR-Verarbeitung in einer zweiten Variante (Variante b) unter Verwendung einer Visualisierung mit Hilfe von Standard-SAR-Prozessoren,
  • 6 das Ablaufdiagramm einer entsprechend dem Verfahren nach der vorliegenden Erfindung ablaufenden geschwindigkeitsadaptiven SAR-Prozessierung mit nachfolgender SAR-Verarbeitung in einer dritten Variante (Variante c) unter Verwendung einer Visualisierung mit Hilfe eines Standard-SAR-Prozessors und adaptiver Nachfokussierung, und
  • 7 vier Bildbeispiele eines Prototyp-Prozessors, wobei 7a ein herkömmliches SAR-Bild mit Straße und eingeblendeten Iso-Geschwindigkeitslinien, 7b ein herkömmliches SAR-Bild des relevanten Bereichs mit Straße und eingeblendeten Iso-Geschwindigkeitslinien, 7c ein geschwindigkeitsadaptiv fokussiertes SAR-Bild mit Straße und eingeblendeten Iso-Geschwindigkeitslinien und 7d ein geschwindigkeitsadaptiv fokussiertes SAR-Bild mit eingeblendeter interferometrischer Phase, Straße und Iso-Geschwindigkeitslinien mit farbig eingeblendeter interferometrischer Phase darstellen.
  • Anhand von 2 wird nun das Prinzip des Verfahrens zur Verarbeitung und zur Darstellung der SAR-Bilder nach der vorliegenden Erfindung veranschaulicht und erläutert. Die Verarbeitung wird durch die neuartige Kombination von zwei Lösungsansätzen stark vereinfacht.
  • Das Straßennetz ist bekannt. Damit kann gemäß dem ersten Lösungsansatz der Lösungsraum stark eingeschränkt werden. Dies geschieht gemäß der vorliegenden Erfindung, indem der durch geografische Koordinaten und Höhen {(x, y, z), ...} gegebene Straßenverlauf in die Azimut/Range-Koordinaten des Radarbildes {(a, r), ...} konvertiert wird.
  • Gemäß dem zweiten Lösungsansatz wird die SAR-Fokussierung durch Rückwärtsprojektion gesteuert. Dabei wird für jeden in Frage kommenden Bildpunkt die Hypothese zu Grunde gelegt, dass er das Abbild eines Fahrzeugs ist. Zunächst wird der Abstand da des Punktes zum angenommenen Verkehrsweg (z. B. Straße, Schiene) in Flugrichtung des flugzeug- oder weltraumgestützten SAR-Sensors bestimmt. Daraus lässt sich sowohl die Position des Fahrzeugs auf der Straße als auch seine Radialgeschwindigkeit Vr ableiten. Dabei ist R der Abstand zwischen Antenne und Punkt und vSAR die Geschwindigkeit der SAR-Antenne: vr = Δα / RvSAR
  • Aus der hypothetischen Radialgeschwindigkeit vr lässt sich für die Fokussierung in Abhängigkeit der Radarparameter leicht der proportionale Korrekturparameter ΔfDC, die Verschiebung des Azimut-Spektrums im Vergleich zu dem der festen Erde, ableiten. Aus dem Einfallswinkel θ des Radarstrahls auf die Erdoberfläche und dem Winkel α zwischen Straße und Flugbahn des Radarsensors kann aus der Radialgeschwindigkeit vr auch die hypothetische Azimutgeschwindigkeit va und daraus unter Zuhilfenahme der Flugbahn und der Radarparameter der zweite Fokussierungsparameter, die Änderung der Frequenzmodulationsrate ΔFM bestimmt werden. Dies gelingt allerdings nicht für Bewegungen nahezu parallel zur Flugbahn mit α ≅ 0.
  • Die beiden aus den hypothetischen Radial- und Azimutgeschwindigkeiten vr und va abgeleiteten Parameter ΔfDC und ΔFM werden in einer Matrix abgelegt und dazu verwendet, die SAR-Fokussierung exakt auf das bewegte Objekt abzustimmen. Damit ist sichergestellt, dass das hypothetische Fahrzeug optimal fokussiert wird. Falls die Hypothese nicht zutrifft, der Bildpunkt also kein Fahrzeug enthält, wird der ruhende Hintergrund defokussiert, was aber für, die Detektion des Fahrzeugs unerheblich ist. In der Karte ist ebenfalls hinterlegt, welche Bereiche überhaupt nicht verarbeitet werden müssen, weil der Versatz zu der nächstliegenden Straße einer unrealistischen Geschwindigkeit entsprechen würde.
  • Die Erzeugung der Fokussiergparameterkarte erfolgt entsprechend den in 4 bis 6 skizzierten Ablaufdiagrammen und wird nachfolgend erklärt:
    • • Der Straßenverlauf wird in Form von Segmenten oder einer 3-dimensionalen Punktefolge {x, y, z} aus einer digitalen Straßendatenbank extrahiert.
    • • Die Segmente bzw. die Punktefolge {(x, y, z), ...} wird in das Azimut/Range-Radarkoordinatensystem {(as, rs,) ...} transformiert.
    • • Aus den Segmenten wird ein nahezu kontinuierlicher Verlauf für jedes Pixel des Radarbildes interpoliert. Dieser Straßenverlauf wird zur visuellen Auswertung zum Schluss zusätzlich in das Radarbild eingeblendet.
    • • Für jedes Pixel in Range wird nun mit Hilfe einer angenommenen maximalen Fahrzeuggeschwindigkeit vmax der maximale Verschiebungsbereich in Azimutrichtung bestimmt. Dieser hängt vom Einfallswinkel θ und vom lokalen Straßenwinkel α ab:
      Figure 00140001
      Für jede Fahrbahnstelle {as, rs} kann das Fahrzeug nur in den Bereich {[as – Δamax, as + Δamax], rs} abgebildet werden.
    • • Für jeden möglichen Bildpunkt {a, r} im Bereich {[as – Δamax, as + Δamax], rs} werden nun die Radialgeschwindigkeit vr und die Azimutgeschwindigkeit va berechnet:
      Figure 00140002
  • Aus diesen Parametern ergeben sich die Fokussierungsfilterparameter ΔfDC und ΔEM, die den SAR-Prozessor für die erwartete Geschwindigkeit justieren. Dieses Fokussierungsparameterpaar wird für jeden in Frage kommenden Bildpunkt in der Nähe der Straße berechnet.
  • Falls das Fahrzeug sich nicht geradlining mit konstanter Geschwindigkeit bewegt, also in Radialrichtung zur Radarantenne beschleunigt ist, so ergeben sich weitere Einflüsse auf die Frequenzmodulationsrate ΔFM mit der Folge, dass das Bild des Fahrzeugs in Azimut defokussiert wird. Wenn solch eine Beschleunigungskomponente durch Kurvenfahrt mit konstanter Geschwindigkeit verursacht wird, so kann die resultierende Defokussierung mit einer leichten Modifikation des Verfahrens kompensiert werden. Hierzu wird durch numerische Ableitung der Radialkomponente der Geschwindigkeit die Radialbeschleunigung errechnet
  • Figure 00150001
  • Diese Radialbeschleuningung bewirkt eine weitere Korrektur der Frequenzmodulationsrate ΔFMacc. Die bei der Fokussierung zu berücksichtigende Frequenzmodulationsratenkorrektur setzt sich dann aus der Korrektur für die Azimutgeschwindigkeit und der Korrektur für die Radaialbeschleunigung zusammen: ΔFMkorr = ΔFM + ΔFMacc
  • Gleichzeitig ergeben sich für den gesamten Straßenverlauf die Extremwerte ΔfDC–min, ΔfDC–max, ΔFM–min, und ΔFM–max.
  • Für die Realisierung der nachfolgenden SAR-Verarbeitung werden gemäß der Erfindung drei Varianten a), b) und c) vorgeschlagen, die nachfolgend anhand von 4 bis 6 beschrieben werden.
  • Variante a): Es wird ein spezieller Zeitbereichskorrelator verwendet, der zur Fokussierung jedes einzelnen Bildpunktes die SAR-Rohdaten mit einem zweidimensionalen Korrelationskern korreliert. Gemäß der Erfindung wird dieser Korrelationskern mit Hilfe der beiden Fokussierungsparameter ΔfDC und ΔFM. für jeden Bildpunkt adaptiert und ermöglicht so ein optimales bewegungsadaptives Bild der Straße und der auf ihr fahrenden Fahrzeuge. Aus den in Frage kommenden Geschwindigkeiten und Positionen wird eine Steuermaske berechnet. Das Bild wird nur für die in der Steuermaske markierten Bereiche fokussiert, um Zeit zu sparen. Dieses Bild ist mathematisch exakt und daher optimal fokussiert. Diese Variante ist in 4 im dortigen Ablaufdiagramm skizziert. Die Berechnung kann unter Umständen länger als die nachfolgend unter Variante b) beschriebene Erzeugung mehrerer Bilder durch Standard-SAR-Prozessoren dauern.
  • Variante b): Hierfür werden Standard-SAR-Prozessoren verwendet. Der Bereich [ΔfDCmin, ΔfDCmax] wird in N Intervalle und der Bereich [ΔFMmin, ΔFMmax] in M Intervalle diskretisiert, wobei die Intervallgrenzen so gewählt sind, dass die Verarbeitungsfehler durch die Diskretisierung tolerabel sind. Die Mittelwerte der Intervalle werden orthogonal in einer 2-dimensionalen Diskretisierungstabelle eingetragen. Für jedes Element der Z = N × M großen Diskretisierungstabelle wird die Matrix der Fokussierungsparameter nach Wertepaaren untersucht, die in dieses Diskretisierungsintervall fallen. Das Ergebnis sind Z Steuermasken, die den SAR-Prozessor dahingehend anweisen, welche Bildteile überhaupt zu verarbeiten sind. Wie in 5 skizziert ist, werden in Z parallel geschalteten oder im Rechner nacheinander ausgeführten Standard-SAR-Prozessoren Z Bilder verarbeitet, die alle benötigten Z Kombinationen von Fokussierungsparametern abdecken. Es ist vorteilhaft, einen Standard-SAR-Prozessor leicht zu modifizieren. Ein solcher modifizierter Standard-SAR-Prozessor führt diejenigen Verarbeitungsschritte, die unabhängig von den Parametern ΔfDC und ΔFM sind, nur einmal durch und alle folgenden Verarbeitungsschritte Z-mal mit adaptierten Parametern. Er verarbeitet für jedes Bild gemäß der Erfindung nur die in der Steuermaske markierten Bildpunkte und spart auf diese Weise sehr viel Rechenzeit. Es liegen nun Z Teilbilder vor, die anhand der Z Steuermasken wieder zu einem Bild zusammengesetzt werden. Das Ergebnis ist ein Bild der Straße und ihrer Umgebung wie bei der Variante a) mit dem Zeitbereichskorrelator.
  • Variante c): Ein Standard-SAR-Prozessor verarbeitet die Szene mit den mittleren Fokussierungsparametern ΔfDC und ΔFM, d. h. normalerweise denjenigen der ruhenden Szene, aus der Diskretisierungstabelle. In diesem ersten Schritt werden die bewegungsbedingten Fokussierungsartefakte in Kauf genommen. Zusätzlich wird eine Modifikation vorgenommen, die darin besteht, dass das Signalspektrum vor der Fokussierung in Azimut soweit repliziert wird, so dass auch bewegte Objekte bezüglich ihres Spektrums und ihrer Range-Migration korrekt vom Prozessor verarbeitet werden. Die Azimut-Bandbreite des Prozessors wird für die folgende Verarbeitung dementsprechend breit dimensioniert. Das fokussierte Bild enthält daher das Bild von Fahrzeugen und Clutter sowie deren spektral und örtlich verschobene Versionen. In einem ortsadaptiven Nachverarbeitungsschritt werden dann die verbleibenden Bewegungsunschärfen ΔFM in Azimutrichtung korrigiert und das Azimut-Spektrum des erwarteten Bildes wird herausgefiltert. Die Filterung kann im Zeitbereich oder durch schnelle Faltung erfolgen. Auch in dieser Variante wird eine Steuermaske verwendet, um nur die notwendigen Bildteile zu verarbeiten. In 6 ist ein schematisches Ablaufdiagramm dieser entsprechend dem Verfahren nach der Erfindung ablaufenden, geschwindigkeitsadaptiven SAR-Prozessierung mit nachfolgender SAR-Verarbeitung gemäß der Variante c) unter Verwendung einer Visualisierung mit Hilfe eines Standard-SAR-Prozessors und adaptiver Nachfokussierung veranschaulicht.
  • Für interferometrische Mehr-Antennen-Systeme wird die Berechnung der Steuermasken und der geschwindigkeitsabhängigen Parameter nur einmal ausgeführt, die eigentliche SAR-Prozessierung für jeden Antennen-Kanal des Systems getrennt.
  • Eine vorteilhafte Weiterbildung des Verfahrens nach der vorliegenden Erfindung besteht in einer neuartigen Darstellung des Bildes zum Zweck der visuellen Auswertung. Dazu werden sowohl der Straßenverlauf als auch die Linien gleicher Geschwindigkeit in das nach der vorstehenden Beschreibung erzeugte SAR-Bild eingeblendet. Die Azimut-Verschiebung für eine konstante Geschwindigkeit viso wird für jeden Punkt der Straße in Range aus der nachfolgenden Gleichung berechnet:
    Figure 00180001
  • Diese Linien gleicher Geschwindigkeit (Iso-Geschwindigkeitslinien) werden wie Höhenlinien in einer topographischen Karte in das SAR-Bild eingeblendet. Eine schematische Veranschaulichung eines geschwindigkeitsabhängig fokussierten SAR-Bildes mit Iso-Geschwindigkeitslinien ist in 3 gezeigt.
  • Im Falle von interferometrischem SAR entsprechend der Along-Track-Interferometrie (ATI) werden die Linien vorteilhaft farbig entsprechend der erwarteten interferometrischen Phase eingezeichnet. Das SAR-Bild wird aus den beiden vorhandenen Antennen-Kanälen generiert und entsprechend der interferometrischen Phasendifferenz zwischen den beiden Kanälen eingefärbt. Die visuelle Auswertung geschieht nun folgendermaßen: Helle Punkte werden als Fahrzeuge angenommen und deren Geschwindigkeit wird aus den nächstliegenden Iso-Geschwindigkeitslinien geschätzt. Im Falle von ATI-Verarbeitung wird zusätzlich die Farbe (= Phase) des Punktes mit der Farbe der nächstliegenden Iso-Geschwindigkeitslinie verglichen. Wenn somit die aus dem Versatz und die aus der Farbe abgeleitete Geschwindigkeit übereinstimmen, dann ist die Wahrscheinlichkeit sehr hoch, dass tatsächlich eine Bewegung mit dieser Geschwindigkeit vorliegt. Zusätzlich kann das Problem der Mehrdeutigkeiten gelöst werden, wenn die Geschwindigkeitsmehrdeutigkeitsintervalle der Phase und der Pulswiederholfrequenz verschieden sind.
  • Ein alternatives Verfahren zur Along-Track-Interferometrie (ATI) ist die Subtraktion der Ergebnisse der beiden Antennenkanale, die auch als DPCA (Displaced Phase Center Antenna) bezeichnet wird. Hier wird nicht die Phasendifferenz, sondern die Differenz der komplexen Bilder beider Antennen gebildet. Idealerweise wird dabei der ruhende Hintergrund (Clutter) eliminiert und die bewegten Fahrzeuge bleiben als helle Punkte erhalten. Auch dieses Bild wird zweckmäßig mit eingeblendetem Straßenverlauf und Linien gleicher Geschwindigkeit dargestellt.
  • Das so erzeugte Bild eignet sich wesentlich besser zur Interpretation als ein mit konstanten Parametern fokussiertes Standard-SAR-Bild oder gar Stapel von Bildern, da nur ein Band entlang der Straße in einem einzigen Bild ausgewertet werden muss und die visuelle Auswertung durch die Geschwindigkeitskonturen unterstützt wird.
  • 7 zeigt vier Bildbeispiele aus Daten der Shuttle-Radar-Topographie-Mission (SRTM). Im herkömmlichen Standard-SAR-Bild nach 7a sind die Fahrzeuge und die Straße nur schwer zu erkennen. Die eingeblendeten Iso-Geschwindigkeitslinien für Vielfache von 50 km/h helfen bereits bei der Interpretation des SAR-Bildes. Noch übersichtlicher ist die Variante nach 7b, bei der nur der relevante Bereich fokussiert wurde und außerdem Rechenzeit gespart wurde. 7c zeigt das gesamte SAR-Bild rein geschwindigkeitsadaptiv fokussiert und mit eingeblendeten Iso-Geschwindigkeitslinien. Der ruhende Clutter im Bereich der Straße erscheint daher mehrdeutig. Für bestimmte Geschwindigkeitsbereiche außerhalb des Clutter-Spektrums sind die Fahrzeuge jetzt viel heller zu sehen. In 7d sind das SAR-Bild und die Iso-Geschwindigkeitslinien mit ihrer interferometrischen Phase eingefärbt. So können auch Mehrdeutigkeiten aufgelöst werden.

Claims (21)

  1. Verfahren zur Verarbeitung und zur Darstellung von mittels flugzeug- oder weltraumgestützten Synthetik-Apertur-Radarsensor-Systemen (SAR) gewonnenen Bodenbildern unter Nutzung eines bekannten Verkehrswegeverlaufs als Vorwissen, wobei eine hypothetische Fahrzeug-Radialgeschwindigkeit vr jedes als Fahrzeug in Betracht kommenden Bildpunktes aus seinem Abstand zu einem Verkehrsweg in Azimutrichtung und die zugehörigen Fahrzeugbildposition durch Azimutprojektion auf den betreffenden Verkehrsweg ermittelt wird und eine hypothetische Fahrzeug-Azimutgeschwindigkeit va aus der hypothetischen Fahrzeug-Radialgeschwindigkeit vr, dem Einfallswinkel Θ des Radarstrahls auf die Erdoberfläche und dem zwischen dem betreffenden Verkehrsweg und der Flugbahn des Radarsensors vorliegenden Winkel α mittels der Funktion
    Figure 00210001
    abgeleitet wird, und dass der durch geographische Koordinaten und Höhen gegebene Verkehrswegeverlauf einer Karte in Azimut/Range-Koordinaten des SAR-Radarbildes konvertiert wird, dass aus der ermittelten hypothetischen Fahrzeug-Radialgeschwindigkeit vr in Abhängigkeit der Radarparameter ein erster, die Verschiebung des Azimut-Spektrums betreffender Fokussierungsparameter ΔfDC abgeleitet wird, dass aus der hypothetischen Fahrzeug-Azimutgeschwindigkeit va unter Zuhilfenahme der Flugbahn des Radarsensors und der Radarparameter ein zweiter, die Änderung der Frequenzmodulationsrate betreffender Fokussierungsparameter ΔFM abgeleitet wird und dass die beiden aus den hypothetischen Radial- und Azimutgeschwindigkeiten vr und va abgeleiteten Fokussierungsparameter ΔfDC und ΔFM in einer Matrix abgelegt und dazu verwendet werden, die SAR-Fokussierung exakt auf das bewegte Objekt abzustimmen.
  2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass in der Karte hinterlegt ist, welche Bereiche nicht verarbeitet werden müssen, weil der Versatz zum nächstliegenden Verkehrsweg einer unrealistischen Geschwindigkeit von Fahrzeugen entspricht.
  3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass der Verkehrswegeverlauf in Form von Segmenten oder einer 3-dimensionalen Punktefolge {x, y, z} aus einer digitalen Verkehrswege-Datenbank extrahiert wird, dass die Segmente bzw. die Punktefolge {(x, y, z), ...} in das Azimut/Range-SAR-Radarkoordinatensystem {(as, rs,) ...} transformiert werden, dass aus den Segmenten bzw. der Punktefolge ein nahezu kontinuierlicher Verlauf für jedes Pixel des Radarbildes interpoliert und dieser Verkehrswegeverlauf zur visuellen Auswertung schließlich zusätzlich in das Radarbild eingeblendet wird, dass für jedes Pixel in der SAR-Range-Koordinate (r) dann mit Hilfe einer angenommenen maximalen Fahrzeuggeschwindigkeit vmax der maximale Verschiebungsbereich in Azimutrichtung bestimmt wird, der entsprechend der Gleichung
    Figure 00220001
    vom Einfallswinkel B und vom zwischen dem betreffenden Verkehrsweg und der Flugbahn des Radarsensors vorliegenden Winkel α abhängig ist, wobei VSAR die Fluggeschwindigkeit und R der Abstand des hypothetischen Fahrzeugbildpunktes vom Radarsensor ist, dass für jede Verkehrswegestelle {as, rs} das betreffende Fahrzeug nur in den Bereich {[as – Δamax, as – Δamax], rs} abgebildet wird und dann für jeden möglichen Bildpunkt {a, r} im Bereich {[as – Δamax, as + Δamax], rs} nun die hypothetische Radialgeschwindigkeit
    Figure 00220002
    und die hypothetische Azimutgeschwindigkeit
    Figure 00220003
    berechnet werden, dass sich aus diesen Geschwindigkeitsparametern die beiden Fokussierungsparameter ΔfDC und ΔFM ergeben, die den SAR-Prozessor für die erwartete Geschwindigkeit justieren, und dass dieses Fokussierungsparameterpaar für jeden in Betracht kommenden Bildpunkt in der Nähe des Verkehrsweges berechnet wird.
  4. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass zur Realisierung der nachfolgenden SAR-Verarbeitung ein spezieller Zeitbereichskorrelator verwendet wird, der zur Fokussierung jedes einzelnen Bildpunktes die SAR-Rohdaten mit einem zweidimensionalen Korrelationskern korreliert, und dass dieser Korrelationskern mit Hilfe der beiden Fokussierungsparameter ΔfDC und ΔFM für jeden Bildpunkt adaptiert wird und damit ein bewegungsadaptives Bild des Verkehrsweges und der auf ihm fahrenden Fahrzeuge erzeugt wird, das mathematisch exakt und daher optimal fokussiert ist.
  5. Verfahren nach Anspruch 4, dadurch gekennzeichnet, dass aus den in Frage kommenden Geschwindigkeiten und Positionen eine Steuermaske berechnet wird und zur Zeitersparnis das Bild nur für die in der Steuermaske markierten Bereiche fokussiert wird.
  6. Verfahren nach einem der Ansprüche 1 bis 3, dadurch gekennzeichnet, dass zur Realisierung der nachfolgenden SAR-Verarbeitung unter Verwendung von Standard-SAR-Prozessoren der Bereich zwischen den Extremwerten des ersten Fokussierungsparameters [ΔfDCmin, ΔfDCmax] in N Intervalle und der Bereich zwischen den Extremwerten des zweiten Fokussierungsparameters [ΔFMmin, ΔFMmax] in M Intervalle diskretisiert werden, wobei die Intervallgrenzen so gewählt sind, dass die Verarbeitungsfehler durch die Diskretisierung tolerabel sind, dass die Mittelwerte der Intervalle orthogonal in einer 2-dimensionalen Diskretisierungstabelle eingetragen werden, dass für jedes Element der Z = N × M großen Diskretisierungstabelle die Matrix der Fokussierungsparameter nach Wertepaaren untersucht wird, die in dieses Diskretisierungsintervall fallen, dass das Ergebnis Z Steuermasken sind, die den SAR-Prozessor anweisen, welche Bildteile zu verarbeiten sind, dass in 2 parallel geschalteten oder im SAR-Rechner nacheinander ausgeführten Standard-SAR-Prozessoren Z Bilder verarbeitet werden, die alle benötigten Z Kombinationen von Fokussierungsparametern abdecken, und dass, nachdem nun Z Teilbilder vorliegen, diese anhand der Z Steuermasken wieder zu einem Bild zusammengesetzt werden, so dass als Ergebnis ein Bild des Verkehrswegs und seiner Umgebung erzeugt wird.
  7. Verfahren nach Anspruch 6, gekennzeichnet durch leichte Modifizierung eines Standard-SAR-Prozessurs, derart, dass ein solcher modifizierter SAR-Prozessor diejenigen Verarbeitungsschritte, die unabhängig von den beiden Fokussierungsparametern ΔfDC und ΔFM sind, nur einmal und alle folgenden Verarbeitungsschritte 2-mal mit adaptierten Fokussierungsparametern durchführt, so dass er in Rechenzeit sparender Weise für jedes Bild nur die in der Steuermaske markierten Bildpunkte verarbeitet.
  8. Verfahren nach einem der Ansprüche 1 bis 3, dadurch gekennzeichnet, dass zur Realisierung der nachfolgenden SAR-Verarbeitung unter Verwendung eines Standard-SAR-Prozessors die Szene mit den Fokussierungsparametern ΔfDC und ΔFM der ruhenden Szene aus einer Diskretisierungstabelle verarbeitet wird, dass in diesem ersten Schritt die bewegungsbedingten Fokussierungsartefakte in Kauf genommen werden und zusätzlich eine Modifikation vorgenommen wird, die darin besteht, dass das Signalspektrum vor der Fokussierung in der SAR-Azimut-Dimension mehrmals repliziert wird, so dass auch bewegte Objekte bezüglich ihres Spektrums und ihrer Range-Migration korrekt vom Prozessor verarbeitet werden, dass die Doppler-Bandbreite des SAR-Prozessors für die folgende Verarbeitung dementsprechend breit dimensioniert wird, so dass das fokussierte Bild daher das Bild von Fahrzeugen und Clutter sowie deren spektral und örtlich verschobene Versionen enthält, und dass in einem ortsadaptiven Nachverarbeitungsschritt dann die verbleibenden, wegen der Annahme eines mittleren Wertes des zweiten Fokussierungsparameters ΔFM entstandenen Bewegungsunschärfen in Azimutrichtung korrigiert werden und das Azimut-Spektrum des erwarteten Bildes herausgefiltert wird.
  9. Verfahren nach Anspruch 8, dadurch gekennzeichnet, dass eine Steuermaske verwendet wird, um nur die notwendigen Bildteile zu verarbeiten.
  10. Verfahren nach einem der Ansprüche 8 bis 9, dadurch gekennzeichnet, dass die Filterung im Zeitbereich erfolgt.
  11. Verfahren nach einem der Ansprüche 8 bis 9, dadurch gekennzeichnet, dass die Filterung durch schnelle Faltung vorgenommen wird.
  12. Verfahren nach einem der Ansprüche 8 bis 11, gekennzeichnet durch leichte Modifizierung eines Standard-SAR-Prozessors, derart, dass ein solcher modifizierter Standard-SAR-Prozessor diejenigen Verarbeitungsschritte, die unabhängig von den beiden Fokussierungsparametern ΔfDC und ΔFM sind, nur einmal und alle folgenden Verarbeitungsschritte Z-mal mit adaptierten Fokussierungsparametern durchführt, so dass er in Rechenzeit sparender Weise für jedes Bild nur die in der Steuermaske markierten Bildpunkte verarbeitet.
  13. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass für interferometrische Mehr-Antennen-Systeme die Berechnung der Masken und der geschwindigkeitsabhängigen Parameter nur einmal ausgeführt wird und die eigentliche SAR-Prozessierung für jeden Antennen-Kanal des Systems getrennt vorgenommen wird.
  14. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass eine Darstellung des Bildes zum Zweck der visuellen Auswertung darin besteht, dass sowohl der Verkehrswegeverlauf als auch die Linien gleicher Geschwindigkeit in das erzeugte SAR-Bild eingeblendet werden.
  15. Verfahren nach Anspruch 14, dadurch gekennzeichnet, dass die Azimut-Verschiebung für eine konstante Geschwindigkeit viso für jeden Punkt des Verkehrswegs in Range gemäß der Gleichung
    Figure 00250001
    berechnet wird, wobei diese Linien gleicher Geschwindigkeit als Iso-Geschwindigkeitslinien wie Höhenlinien in einer topographischen Karte in das SAR-Bild eingeblendet werden.
  16. Verfahren nach einem der Ansprüche 14 oder 15, dadurch gekennzeichnet, dass im Falle von interferometrischem SAR mit Along-Track-Interferometrie (ATI) die Linien farbig entsprechend der erwarteten interferometrischen Phase eingezeichnet werden.
  17. Verfahren nach Anspruch 16, dadurch gekennzeichnet, dass das SAR-Bild aus beiden Antennen-Kanälen generiert und entsprechend der interferometrischen Phasendifferenz beider Kanäle eingefärbt wird, dass die visuelle Auswertung in der Weise vorgenommen wird, dass helle Punkte als Fahrzeuge angenommen werden und deren Geschwindigkeit aus den nächstliegenden Iso-Geschwindigkeitslinien geschätzt wird, und dass im Falle von ATI-Verarbeitung zusätzlich die Farbe, d. h. Phase, des Punktes mit der Farbe der nächstliegenden Iso-Geschwindigkeitslinie verglichen wird, und dass dann, wenn die aus dem Versatz und die aus der Farbe abgeleiteten Geschwindigkeiten übereinstimmen, angenommen wird, dass eine Bewegung mit dieser Geschwindigkeit vorliegt, da die Wahrscheinlichkeit sehr hoch ist, dass tatsächlich Bewegung mit dieser Geschwindigkeit vorliegt.
  18. Verfahren nach einem der Ansprüche 14 oder 15, dadurch gekennzeichnet, dass im Falle von mehrkanaligem SAR mit Subtraktion der Ergebnisse zweier Antennenkanäle die Differenz der komplexen Bilder zweier Antennen gebildet wird.
  19. Verfahren nach Anspruch 18, dadurch gekennzeichnet, dass der ruhende Hintergrund eliminiert wird und die bewegten Fahrzeuge als helle Punkte erhalten bleiben und dass das Bild mit eingeblendetem Verkehrswegeverlauf und Linien gleicher Geschwindigkeit dargestellt wird.
  20. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass eine Defokussierung, die aus einer durch Kurvenfahrt des Fahrzeugs mit konstanter Geschwindigkeit verursachten Beschleunigungskomponente in Radialrichtung zur SAR-Sensorantenne resultiert, durch eine Modifikation kompensiert wird, die darin besteht, dass durch numerische mathematische Ableitung der Radialkomponente vr der Geschwindigkeit nach der Zeit t die Radialbeschleunigung v .r gemäß der Gleichung
    Figure 00270001
    berechnet wird, dass aus dieser berechneten Radialbeschleunigung v .r ein Korrekturwert der Frequenzmodulationsrate ΔFMacc fur den zweiten Fokussierungsparameter abgeleitet wird und dass sich die bei der Fokussierung zu berücksichtigende Frequenzmodulationsratenkorrektur, d. h. der gesamte korrigierte zweite Fokussierungsparameter ΔFMkorr, aus der Korrektur für die Azimutgeschwindigkeit, d. h. dem zweiten Fokussierungsparameter ΔFM, und dem Korrekturwert der Frequenzmodulationsrate, d. h. dem berechneten Korrekturwert ΔFMacc, gemäß der Gleichung ΔFMkorr = ΔFM + ΔFMacc zusammensetzt.
  21. Verwendung des Verfahrens nach einem der vorhergehenden Ansprüche bei der Erfassung und dem Einsatz von Verkehrsdaten für Wissenschaft, Wirtschaft und Sicherheit und der Erstellung von übersichtlichen Bildern von Verkehrsflüssen.
DE102006009121A 2006-02-24 2006-02-24 Verfahren zur Verarbeitung und Darstellung von mittels Synthetik-Apertur-Radarsystemen (SAR) gewonnen Bodenbildern Expired - Fee Related DE102006009121B4 (de)

Priority Applications (4)

Application Number Priority Date Filing Date Title
DE102006009121A DE102006009121B4 (de) 2006-02-24 2006-02-24 Verfahren zur Verarbeitung und Darstellung von mittels Synthetik-Apertur-Radarsystemen (SAR) gewonnen Bodenbildern
IT000006A ITSV20070006A1 (it) 2006-02-24 2007-02-13 Procedimento per l'elaborazione e la rappresentazione di immagini di terra ottenute mediante sistemi radar ad apertura sintetica (sar)
IT000007A ITSV20070007A1 (it) 2006-02-24 2007-02-14 Procedimento per l'elaborazione e la rappresentazione di immagini di terra ottenute mediante sistemi radar ad apertura sintetica (sar)
CA2579898A CA2579898C (en) 2006-02-24 2007-02-23 Method for the processing and representing of ground images obtained by synthetic aperture radar systems (sar)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102006009121A DE102006009121B4 (de) 2006-02-24 2006-02-24 Verfahren zur Verarbeitung und Darstellung von mittels Synthetik-Apertur-Radarsystemen (SAR) gewonnen Bodenbildern

Publications (2)

Publication Number Publication Date
DE102006009121A1 DE102006009121A1 (de) 2007-08-30
DE102006009121B4 true DE102006009121B4 (de) 2012-05-31

Family

ID=38319907

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102006009121A Expired - Fee Related DE102006009121B4 (de) 2006-02-24 2006-02-24 Verfahren zur Verarbeitung und Darstellung von mittels Synthetik-Apertur-Radarsystemen (SAR) gewonnen Bodenbildern

Country Status (3)

Country Link
CA (1) CA2579898C (de)
DE (1) DE102006009121B4 (de)
IT (2) ITSV20070006A1 (de)

Families Citing this family (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7889133B2 (en) 1999-03-05 2011-02-15 Itt Manufacturing Enterprises, Inc. Multilateration enhancements for noise and operations management
US7777675B2 (en) 1999-03-05 2010-08-17 Era Systems Corporation Deployable passive broadband aircraft tracking
US7782256B2 (en) 1999-03-05 2010-08-24 Era Systems Corporation Enhanced passive coherent location techniques to track and identify UAVs, UCAVs, MAVs, and other objects
US8203486B1 (en) 1999-03-05 2012-06-19 Omnipol A.S. Transmitter independent techniques to extend the performance of passive coherent location
US8446321B2 (en) 1999-03-05 2013-05-21 Omnipol A.S. Deployable intelligence and tracking system for homeland security and search and rescue
US7570214B2 (en) 1999-03-05 2009-08-04 Era Systems, Inc. Method and apparatus for ADS-B validation, active and passive multilateration, and elliptical surviellance
US7612716B2 (en) 1999-03-05 2009-11-03 Era Systems Corporation Correlation of flight track data with other data sources
US7739167B2 (en) 1999-03-05 2010-06-15 Era Systems Corporation Automated management of airport revenues
US7908077B2 (en) 2003-06-10 2011-03-15 Itt Manufacturing Enterprises, Inc. Land use compatibility planning software
US7667647B2 (en) 1999-03-05 2010-02-23 Era Systems Corporation Extension of aircraft tracking and positive identification from movement areas into non-movement areas
US7965227B2 (en) 2006-05-08 2011-06-21 Era Systems, Inc. Aircraft tracking using low cost tagging as a discriminator
DE102009007704A1 (de) 2009-02-05 2010-08-12 Daimler Ag Verfahren und Vorrichtung zur Ermittlung einer voraussichtlichen Fahrspur eines Fahrzeugs
DE102009013768A1 (de) 2009-03-18 2010-10-14 Eads Deutschland Gmbh Verfahren und Einrichtung zum Ermitteln von Aspektwinkeln
CN103377558A (zh) * 2012-04-26 2013-10-30 鸿富锦精密工业(深圳)有限公司 交通流量管控系统及方法
CN108305274B (zh) * 2018-03-08 2021-11-23 中国民航大学 PolSAR图像多特征融合的飞机目标检测方法
CN108549063B (zh) * 2018-06-29 2021-11-26 北京航空航天大学 一种星载合成孔径雷达图像中运动目标航迹向速度快速估计方法
CN110389366B (zh) * 2019-08-05 2021-03-30 中国人民解放军军事科学院国防科技创新研究院 一种基于多源sar卫星的海上目标运动估计方法
KR102653404B1 (ko) * 2020-05-08 2024-04-02 주식회사 에이치엘클레무브 차량의 제어 장치 및 제어 방법
DE102020207879A1 (de) * 2020-06-25 2021-12-30 Robert Bosch Gesellschaft mit beschränkter Haftung Verfahren zum Betrieb eines Radarsensors in einem Kraftfahrzeug
CN113960599B (zh) * 2021-09-17 2024-05-31 北京工商大学 一种扫描模式sar成像重聚焦方法
CN114296080A (zh) * 2021-12-27 2022-04-08 天津天地伟业信息系统集成有限公司 一种交通数据采集方法及计算机可读存储介质

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4423899C1 (de) * 1994-07-08 1995-11-02 Deutsche Forsch Luft Raumfahrt Verfahren zur Detektion, Lokalisierung und Geschwindigkeitsbestimmung von Bewegtzielen aus Radar-Rohdaten eines von einem Träger mitgeführten, kohärenten, ein- oder mehrkanaligen Abbildungssystems

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4423899C1 (de) * 1994-07-08 1995-11-02 Deutsche Forsch Luft Raumfahrt Verfahren zur Detektion, Lokalisierung und Geschwindigkeitsbestimmung von Bewegtzielen aus Radar-Rohdaten eines von einem Träger mitgeführten, kohärenten, ein- oder mehrkanaligen Abbildungssystems

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
HINZ, S. [u.a.]: Spaceborn Traffic Monitoring with Dual Channel Synthetic Aperture Radar – Theory and Experiments. In: Proceedings of the 2005 IEEE Computer Society Conference on Computer Vison and Pattern Recognition CVPR’05, 2005.
HINZ, S. [u.a.]: Spaceborn Traffic Monitoring with Dual Channel Synthetic Aperture Radar - Theory and Experiments. In: Proceedings of the 2005 IEEE Computer Society Conference on Computer Vison and Pattern Recognition CVPR'05, 2005. *
SHARMA, J.: The influence of Target Acceleration on Dual-Channel SAR-GMTI (Synthetic Aperture Radar Ground Moving Target Indication) Data. UCGE Reports Number 20206, Oktober 2004, University of Calgary, S. 186, 189. *
SUCHANDT, S., [u.a.]: Results from an Airborne SAR GMTI Experiment supporting TerraSAR-X Traffic Processor Development. In: Geoscience and Remote Sensing Symposium, 2005, IGARSS '05, Proceedings, 2005 IEEE International, Volume 4, 25-29 July 2005, S. 2949 - 2952. *
SUCHANDT, S., [u.a.]: Results from an Airborne SAR GMTI Experiment supporting TerraSAR-X Traffic Processor Development. In: Geoscience and Remote Sensing Symposium, 2005, IGARSS '05, Proceedings, 2005 IEEE International, Volume 4, 25-29 July 2005, S. 2949 – 2952.
THOMPSON, A. A., LIVINGSTONE, C. E.: Moving Target Performance for RADARSAT-2. In: Geoscience and Remote Sensing Symposium, 2000, Proceedings, IGARSS 2000, IEEE 2000 International, Volume 6, 24-28 July 2000, S. 2599 - 2601. *
THOMPSON, A. A., LIVINGSTONE, C. E.: Moving Target Performance for RADARSAT-2. In: Geoscience and Remote Sensing Symposium, 2000, Proceedings, IGARSS 2000, IEEE 2000 International, Volume 6, 24-28 July 2000, S. 2599 – 2601.

Also Published As

Publication number Publication date
CA2579898C (en) 2012-05-22
CA2579898A1 (en) 2007-08-24
ITSV20070007A1 (it) 2007-08-25
ITSV20070006A1 (it) 2007-08-25
DE102006009121A1 (de) 2007-08-30

Similar Documents

Publication Publication Date Title
DE102006009121B4 (de) Verfahren zur Verarbeitung und Darstellung von mittels Synthetik-Apertur-Radarsystemen (SAR) gewonnen Bodenbildern
EP0406877B1 (de) Verfahren zur Extraktion von Bewegungsfehlern eines ein kohärentes Abbildungsradarsystem mit führenden Trägers aus Radar-Rohdaten und Einrichtung zur Durchführung des Verfahrens
DE19757309C1 (de) Verfahren zur Verarbeitung von Spotlight SAR-Rohdaten
DE102013104443B4 (de) Verkehrsüberwachungssystem zur Geschwindigkeitsmessung und Zuordnung von bewegten Fahrzeugen bei einem Mehrziel-Aufnahmemodul
DE69615587T2 (de) Verfahren zur momentabschätzung und ermittlung von merkmalen von spektren als entfernungs- oder zeitfunktion bestimmenden vorrichtungen
EP4211490A1 (de) Verfahren, radarsystem und fahrzeug zur signalverarbeitung von radarsignalen
DE3922427C1 (de)
DE102019114723A1 (de) Abschwächen von schwingungen in einem radarsystem auf einer beweglichen plattform
EP0784213A2 (de) Verfahren zur Bestimmung des Abstandes und/oder der Differenzgeschwindigkeit zwischen einem Radarsensor und einem oder mehreren Objekten sowie Anordnung zum Durchführen des Verfahrens
DE102018202864A1 (de) Verfahren und System für Synthetische-Apertur-Radarsignalverarbeitung
DE102018202294A1 (de) Schätzung von kartesischen Geschwindigkeiten von ausgedehnten Radarobjekten mit einem Radarsensor
DE102018000517A1 (de) Verfahren zur radarbasierten Vermessung und/oder Klassifizierung von Objekten in einer Fahrzeugumgebung
EP2409174B1 (de) Verfahren und einrichtung zum ermitteln von aspektwinkelverläufen
DE4311754C1 (de) Verfahren zur Extraktion von Bewegungsfehlern eines ein kohärentes Abbildungsradarsystem mitführenden Trägers aus Radar-Rohdaten und Einrichtung zur Durchführung des Verfahrens
DE202019105270U1 (de) System zur Herleitung von Seestatusparametern
DE10341893A1 (de) Verfahren zur Verringerung des Dopplerzentroids bei einem kohärenten impuls-Radarsystem
DE102022208465A1 (de) Verfahren zur Korrektur eines Radarsignals zum Bestimmen einer synthetischen Apertur, Computerprogramm, Vorrichtung und Fahrzeug
DE102019106293B3 (de) Verfahren zum Detektieren eines Flugobjekts und Passives Radardetektionssystem zum Detektieren eines Flugobjekts
EP3936892A1 (de) Vorrichtung und verfahren zum kompensieren von doppleraufweitungen
DE102011016337B4 (de) Verfahren und Vorrichtung zur Abbildung von mit Clutter überlagerten bewegten Objekten mit Radar oder Sonar
DE102006051091B4 (de) Verfahren zur Objekterkennung von Fahrzeugen mittels Nahbereichsdetektion
DE102022133858B4 (de) Verfahren zur rechnergestützten Verarbeitung von SAR-Rohdaten
DE102013213304B3 (de) Verfahren zur rechnergestützten Verarbeitung von SAR-Daten
DE102022117088A1 (de) Verfahren zur Erkennung und Auflösung von Phasenmehrdeutigkeiten in interferometrischen SAR-Daten
DE102023121112A1 (de) Erkennung von Radarreflektionen mit gefalteten neuronalen Netzwerk-Kerneln, die an einen Reflexions-Geist-Offset angepasst sind

Legal Events

Date Code Title Description
OP8 Request for examination as to paragraph 44 patent law
R016 Response to examination communication
R016 Response to examination communication
R018 Grant decision by examination section/examining division
R020 Patent grant now final

Effective date: 20120901

R119 Application deemed withdrawn, or ip right lapsed, due to non-payment of renewal fee