EP1889092B1 - Verfahren und vorrichtung zum korrigieren von wetterdaten sowie computerprogrammprodukt - Google Patents

Description

• Die Erfindung betrifft ein Verfahren und eine Vorrichtung zum Korrigieren von Wetterdaten sowie ein Computerprogrammprodukt nach dem Oberbegriff des Anspruchs 1, 18 bzw. 20.
• Bei der Erfassung von Wetterdaten mittels eines bodengestützten Wetterradars beispielsweise zur Wetterbeobachtung in der Nähe von Flughäfen werden neben Wetterechos, d.h. Echos von Wolken und Niederschlag usw., auch Bodenechos, d.h. Echos von Hindernissen wie Gebäuden, Bergen, Bäumen usw. vom Wetterradar aufgenommen. Zum Herausfiltern der auch als Clutter bezeichneten Bodenechos aus der Intensitätsverteilung der Echos sind zwei Ansätze bekannt. Einerseits kann aus der Intensitätsverteilung einer mit Clutter verunreinigten Messung eine Cluttermap subtrahiert werden, in der eine Intensitätsverteilung nur mit Boden- und ohne Wetterechos gespeichert ist und die unter geeigneten Wetterbedingungen aufgenommen worden ist. Andererseits können durch Vergleich der Meßwerte mit Werten der Cluttermap diejenigen Bereiche, die hinreichend mit Clutter verunreinigt sind, bestimmt und durch Inter- oder Extrapolation aus Bereichen, die nicht hinreichend mit Clutter verunreinigt sind, ersetzt werden.
• Sowohl bei der Subtraktion als auch bei der Inter- oder Extrapolation sind diverse Schwellenwerte und komplexe Algorithmen zu verwenden, die zudem nur die Intensität des Echos korrigieren. Für die Wetterdaten differentielle Reflektivität sind keine geeigneten Korrekturmöglichkeiten unter Beibehaltung eines gültigen Werts bekannt.
• Aus J.E. Evans et al.:TERMINAL DOPPLER WEATHER RADAR CLUTTER CONTROL, IEEE International Radar Conference, 1990, S. 12-16, ist das Eliminieren bzw. Streichen von Wetterdaten bekannt. Die Wetterdaten werden in Form einer Windgeschwindigkeit unter Zuhilfenahme einer zuvor gespeicherten Cluttermap eliminiert.
• Aus L.P. Ligthart et al.: "An X-band solid-state FM-CW weather radar", IEEE Proceedings, Bd. 137, Heft 6, 1990, S. 418-426 ist ein FMCW-Radar bekannt, dessen Spektren korrigiert werden, woraus sich Geschwindigkeit und spektrale Breite berechnen lassen.
• Der Erfindung liegt daher die Aufgabe zugrunde, ein Verfahren und eine Vorrichtung zum Korrigieren von Wetterdaten sowie ein Computerprogrammprodukt nach dem Oberbegriff des Anspruchs 1, 18 bzw. 20 zu schaffen, mit denen die Wetterdaten differentielle Reflektivität zuverlässig und einfach unter Beibehaltung eines gültigen Werts korrigierbar sind.
• Diese Aufgabe wird entsprechend den Merkmalen der Ansprüche 1, 18 bzw. 20 gelöst.
• Hierdurch können auch Wetterdaten vom Typ Radialgeschwindigkeit, spektrale Breite und differentielle Reflektivität automatisiert von Clutter bereinigt werden. Die automatisierte Korrektur dieser Wetterdaten ermöglicht beispielsweise eine Verbesserung der Zuverlässigkeit der Unwettererkennung und bei hydrologischen Anwendungen, die bisher, d.h. wenn nur die Intensitätswerte korrigiert werden, häufig inkorrekte Ergebnisse liefern.
• Weitere Ausgestaltungen der Erfindung sind der nachfolgenden Beschreibung und den Unteransprüchen zu entnehmen.
• Die Erfindung wird nachstehend anhand von in den beigefügten Figuren dargestellten Ausführungsbeispielen näher erläutert.
• Fig. 1 illustriert ein bodengestütztes radarbasiertes Fernerkundungsgerät.
• Fig. 2A, 2B illustriert jeweils eine Interpolation.
• Fig. 3 illustriert Abschattungen.
• Fig. 4 illustriert die Berechnung eines Höhenunterschieds.
• Das in Fig. 1 dargestellte Fernerkundungsgerät in Form eines Wetterradars 1 auf der Erdoberfläche 2 umfaß einen dreh- und/oder schwenkbaren Sende-/Empfangsspiegel 3 für einen zweckmäßigerweise gepulsten elektromagnetischen Strahl einer Wellenlänge insbesondere im Mikrowellenbereich mit einer Hauptstrahlungskeule 4, die mit ihrer 3dB-Breite θ dargestellt ist, und mehreren Nebenstrahlungskeulen 5.
• Im Betrieb wird der dargestellte Wetterradar 1 um eine Hochachse 6 gedreht und/oder um eine sich zweckmäßigerweise mitdrehende Querachse 7 verschwenkt, um einen Raumwinkel oder gegebenenfalls einen Flächenwinkel abzufahren. Dabei werden über den Sende-/Empfangsspiegel 3 elektromagnetische Impulse ausgesendet, die von Objekten wie Wolken, Niederschlag, Gebäuden oder der Erdoberfläche 2 gegebenenfalls teilweise reflektiert und/oder absorbiert werden und am Sende-/Empfangsspiegel 3 mit Amplitude, Phase und/oder Polarisation detektierbare Echoimpule erzeugen. Die Echoimpulse werden als Funktion des Dreh- und/oder Schwenkwinkels sowie des aus der Laufzeit zwischen Aussenden eines Impulses und Empfangen des Echoimpulses berechneten Abstandes eines echogebenden Objekts erfaßt und zweckmäßigerweise für eine Auswertung gespeichert. Auf diese Weise werden Radarechodaten erhalten, die im Rahmen der erzielbaren Auflösung Informationen über jeden Punkt im abgefahrenen Raumbereich enthalten.
• Echogebend sind nicht nur für die Wettererkundung gewünschte Ziele wie Wolken und Niederschlag, sondern auch unerwünschte Ziele wie Gebäude oder die Erdoberfläche 2. Bodenechos erzeugende Bereiche sind für die Hautpstrahlungskeule 4 in Zeile A und für die Nebenstrahlungskeulen 5 in Zeile B der Fig. 1 schraffiert dargestellt. Zudem treten wegen der unerwünschten Ziele teilweise und vollständige Abschattungen auf, die in Zeilen C bzw. D der Fig. 1 schraffiert dargestellt sind. Für die Wettererkundung gewünschte Ziele werden bei der in Fig. 1 dargestellten Position des Wetterradars 1 nur im in Zeile E schraffierten Bereich weder durch Bodenechos gestört noch abgeschattet.
• Zum Korrigieren der vom Wetterradar 1 auf Basis von Radarechodaten gelieferten Wetterdaten Radialgeschwindigkeit V, spektrale Breite W und/oder differentielle Reflektivität ZDR wird zunächst eine Cluttermap Z_CI erstellt, d.h. eine Intensitätsverteilung von Bodenechos ohne Wetterechos. Diese kann unter wetterechofreien Bedingungen aufgenommen und gegebenenfalls per Hand oder automatisiert, insbesondere durch Mittelung verschiedener Aufnahmen, optimiert werden. Die Cluttermap Z_CI ist wie die Wetterdaten V, W und ZDR zweckmäßigerweise dreidimensional und eine Funktion des Abstands r, des Drehwinkels α und des Schwenkwinkels β, kann jedoch auch zweidimensional als Funktion des Abstands r und des Drehwinkels α oder nur eindimensional als Funktion des Abstands r sein. Gegebenenfalls werden kartesische Koordinaten verwendet. Die Koordinatenschrittgröße der Cluttermap Z_CI stimmt zweckmäßigerweise mit der Koordinatenschrittgröße der Wetterdaten V, W, ZDR überein.
• Die Wetterdaten werden wie folgt korrigiert.
• Eine korrigierte Radialgeschwindigkeit V_Wea wird entsprechend $${\mathrm{V}}_{\mathrm{Wea}}\mathrm{=}\mathrm{V}\left({\mathrm{Z}}_{\mathrm{Cl}}\mathrm{+}{\mathrm{Z}}_{\mathrm{Wea}}\right)\mathrm{/}{\mathrm{Z}}_{\mathrm{Wea}}$$

  

  
  berechnet, wobei V die noch nicht korrigierte Radialgeschwindigkeit, Z_CI die Intensität aus der Cluttermap und Z_Wea die Intensität des Wettersignals ist.
• Eine korrigierte spektrale Breite W_Wea wird entsprechend $${\mathrm{W}}_{\mathrm{Wea}}\mathrm{=}\mathrm{W}\left({\mathrm{Z}}_{\mathrm{Cl}}\mathrm{+}{\mathrm{Z}}_{\mathrm{Wea}}\right)\mathrm{/}{\mathrm{Z}}_{\mathrm{Wea}}$$

  

  
  berechnet, wobei W die noch nicht korrigierte spektrale Breite, Z_CI die Intensität aus der Cluttermap und Z_Wea die Intensität des Wettersignals ist.
• Eine korrigierte differentielle Reflektivität ZDR_Wea wird entsprechend $${\mathrm{ZDR}}_{\mathrm{Wea}}\mathrm{=}\mathrm{ZDR}\left({\mathrm{Z}}_{\mathrm{Cl}}\mathrm{+}{\mathrm{Z}}_{\mathrm{Wea}}\right)\mathrm{/}{\mathrm{Z}}_{\mathrm{Wea}}\mathrm{-}{\mathrm{ZDR}}_{\mathrm{Cl}}{\mathrm{Z}}_{\mathrm{Cl}}\mathrm{/}{\mathrm{Z}}_{\mathrm{Wea}}$$

  

  
  berechnet, wobei ZDR die noch nicht korrigierte differentielle Reflektivität, Z_CI die Intensität aus der Cluttermap, Z_Wea die Intensität des Wettersignals und ZDR_CI die differentielle Reflektivität aus der Cluttermap ist.
• Wenn die Bedingung $${\mathrm{Z}}_{\mathrm{Wea}}\mathrm{<}{\mathrm{Z}}_{\mathrm{Cl}}$$

  

  
  und insbesondere $${\mathrm{Z}}_{\mathrm{Wea}}\ll {\mathrm{Z}}_{\mathrm{Cl}}$$

  

  
  erfüllt ist, können die Gleichungen (1) bis (3) instabile Ergebnisse liefern. In diesem Fall werden die jeweiligen Wetterdaten zweckmäßigerweise durch gewichtete Mittelwertbildung aus wenigstens zwei unabhängig erhaltenen korrigierten Werten wie folgt berechnet.
• Für die Radialgeschwindigkeit V können wenigstens zwei der folgenden Werte zweckmäßigerweise gewichtet kombiniert verwendet werden:
• Als ein Wert W₁ kann V_Wea gemäß Gleichung (1) verwendet werden, gegebenenfalls mit einer Gewichtungsfunktion g₁, die mit zunehmendem Faktor Z_Wea / Z_CI zunimmt.
• Weitere Werte W₂, W₃, W₄ können mittels Interpolation bestimmt werden. Zweckmäßigerweise werden die Werte W₂, W₃, W₄ aus benachbarten und von Clutter unbeeinträchtigten Daten für V entlang jeweils einer der drei Achsen des verwendeten Koordinatensystems, beispielsweise entlang des Drehwinkels α, des Schwenkwinkels β und/oder des Abstands r interpoliert.
• Die Interpolation ist beispielhaft in Fig. 2A für die Achse α illustriert. Der dort dargestellte zweidimensionale Koordinatenbereich 8 weist in der Cluttermap Einträge infolge Bodenechos auf. Der koordinatenmäßig entsprechende Bereich der gemessenen Radialgeschwindigkeit V wird ersetzt durch interpolierte Werte. Dargestellt ist beispielhaft der hier zwischen den gemessenen Werten V (r₁, α₀) und V (r₁, α₂) interpolierte Wert V_i (r₁, α₁): $${\mathrm{W}}_{\mathrm{2}}\mathrm{=}{\mathrm{V}}_{\mathrm{i}}\left({\mathrm{r}}_{\mathrm{1}},{\mathrm{\alpha }}_{\mathrm{1}}\right)\mathrm{=}\mathrm{f}\left({\mathrm{\alpha }}_{\mathrm{1}}\right)\mathrm{V}\left({\mathrm{r}}_{\mathrm{1}},{\mathrm{\alpha }}_{\mathrm{0}}\right)\mathrm{+}\left(\mathrm{1}\mathrm{-}\mathrm{f}\left({\mathrm{\alpha }}_{\mathrm{1}}\right)\right)\mathrm{V}\left({\mathrm{r}}_{\mathrm{1}},{\mathrm{\alpha }}_{\mathrm{2}}\right)\mathrm{.}$$

  
• Für den Interpolationsfaktor f(α₁) ist ein Wert zwischen 0 und 1 entweder fest oder entsprechend einer Funktion f(α) vorgebbar, wobei die Funktion f(α) zweckmäßigerweise linear zwischen f(α₀) = 1 und f(α₂) = 0 verläuft, so daß für α₁ = ½ (α₀ + α₂) der Interpolationsfaktor f(α₁) = ½ beträgt.
• Gegebenenfalls ist für den Wert W₂ eine Gewichtungsfunktion g₂ vorgesehen, die zweckmäßigerweise mit zunehmendem Abstand, z.B. Δ = α₂ - α₁, im dreidimensionalen Raum zwischen den Koordinatentripeln der zur Interpolation verwendeten Werte von V abnimmt.
• Für den Wert W₃ kann eine Interpolation entlang der Achse β analog zur Interpolation entlang der Achse α für den Wert W₂ durchgeführt werden. Dies ist zwar beispielsweise im in Fig. 1 in Zeile C dargestellten Bereich nicht möglich, da dort bodenseitig kein unbeeinträchtigter Bereich vorhanden ist. Können jedoch in Ausnahmefällen beispielsweise unterhalb sich quer erstreckender oder punktueller Hindernisse Meßwerte aufgenommen werden, kann der Wert W₃ sinnvoll berechnet und verwendet werden.
• Für den Wert W₄ kann wie in Fig. 2B dargestellt eine Interpolation entlang der Achse r durchgeführt werden. Der beispielhaft dargestellte Wert V_i (r₁, α₁) ist hier zwischen den gemessenen Werten V (r₀, α₁) und V (r₂, α₁) interpoliert: $${\mathrm{W}}_4\mathrm{=}{\mathrm{V}}_{\mathrm{i}}\left({\mathrm{r}}_{\mathrm{1}},{\mathrm{\alpha }}_{\mathrm{1}}\right)\mathrm{=}\mathrm{f}\left({\mathrm{r}}_{\mathrm{1}}\right)\mathrm{V}\left({\mathrm{r}}_{\mathrm{0}},{\mathrm{\alpha }}_{\mathrm{1}}\right)\mathrm{+}\left(\mathrm{1}\mathrm{-}\mathrm{f}\left({\mathrm{r}}_{\mathrm{1}}\right)\right)\mathrm{V}\left({\mathrm{r}}_{\mathrm{2}},{\mathrm{\alpha }}_{\mathrm{1}}\right)\mathrm{.}$$

  
• Für den Interpolationsfaktor f(r₁) ist ein Wert zwischen 0 und 1 entweder fest oder entsprechend einer Funktion f(r) vorgebbar, wobei die Funktion f(r) zweckmäßigerweise linear zwischen f(r₀) = 1 und f(r₂) = 0 verläuft, so daß für r₁ = ½ (r₀ + r₂) der Interpolationsfaktor f(r₁) = ½ beträgt.
• Gegebenenfalls ist für den Wert W₄ eine Gewichtungsfunktion g₄ vorgesehen, die zweckmäßigerweise mit zunehmendem Abstand, z.B. Δ = r₂ - r₀, im dreidimensionalen Raum zwischen den Koordinatentripeln der zur Interpolation verwendeten Werte von V abnimmt.
• Ein weiterer Wert W₅ kann basierend auf einer vertikalen Extrapolation berechnet werden, insbesondere bei einer Abschattung vertikal von oben nach unten wie für den Bereich in Zeile C der Fig. 1. Die vertikale Extrapolation wird zweckmäßigerweise entlang β oder gegebenenfalls entlang der kartesischen Koordinate z jeweils vom unbeeinträchtigten Bereich oberhalb der Abschattung nach unten in den abgeschatteten Bereich durch die einzelnen (α,r)-Ebenen, die durch die β-Koordinaten im abgetasteten Raum festgelegt sind, hindurch durchgeführt. Gegebenenfalls wird für den Wert W₅ eine Gewichtungsfunktion g₅ vorgesehen, die zweckmäßigerweise mit zunehmendem vertikalen Abstand im dreidimensionalen Raum zum zur Extrapolation verwendeten unbeeinträchtigten Bereich ab- und/oder mit einem Zuverlässigkeitsmaß eines zuvor bestimmten Vertikalprofils der Radialgeschwindigkeit zunimmt. Sofern ein Wert W₃ für eine vertikale Interpolation berechenbar ist, kann auf W₅ verzichtet bzw. die Gewichtungsfunktion g₅ auf einen kleinen Wert oder 0 gesetzt werden.
• Aus den unabhängig berechneten korrigierten Werten W₁ bis W₅ können eine beliebige Anzahl i linear oder mittels der jeweiligen Gewichtungsfunktion zu einem einzigen korrigierten Wert V_Wea für die Radialgeschwindigkeit an dem jeweiligen Punkt kombiniert werden: $${\mathrm{V}}_{\mathrm{Wea}}\mathrm{=}{\displaystyle \mathrm{\sum }_{\mathit{i}}}{\mathrm{g}}_{\mathrm{i}}{\mathrm{W}}_{\mathrm{i}}\mathrm{,}\mathrm{i}\mathrm{=}\mathrm{1}\mathrm{,}\mathrm{2}\mathrm{,}\mathrm{3}\mathrm{,}\mathrm{\dots }$$

  
• Für die weiteren Wetterdaten spektrale Breite und differentielle Reflektivität können als Wert W₁ W_Wea gemäß Gleichung (2) bzw. ZDR_Wea gemäß Gleichung (3) verwendet sowie analog zur Radialgeschwindigkeit berechnete Werte W₂ bis W₅ und Gewichtungsfunktionen g₁ bis g₅ verwendet werden.
• Die Gewichtungsfunktionen insbesondere für die Inter- und Extrapolation können zudem insbesondere diskontinuierlich auf 0 oder 1, gegebenenfalls 1/i, abfallen bzw. ansteigen, wenn empirisch vorgebbare Schwellenwerte über- bzw. unterschritten werden, beispielsweise der Abstand zwischen den zur Interpolation verwendeten Werten zu groß bzw. hinreichend klein ist.
• Im in Zeile C der Fig. 1 schraffierten Bereich ist die Hauptstrahlungskeule 4 teilweise abgeschattet. Ein Querschnitt durch die Hauptstrahlungskeule 4 für verschiedene Werte der Abschattung ist in Fig. 3 illustriert. Bei der Korrektur der Wetterdaten wird erfindungsgemäß der in Fig. 3, 4 illustrierte Höhenunterschied Δβ zwischen der mittleren Höhe β_m des verbleibenden Teils der teilweise abgeschatteten Hauptstrahlungskeule 4 und der Mitte β₀ der Hauptstrahlungskeule 4 derart berücksichtigt, daß ein Meßwert für ein Wetterdatum nicht der Koordinate β₀ sondern der Koordinate β_m = β₀ + Δβ zugeordnet wird.
• Die mittlere Höhe β_m des verbleibenden Teils der teilweise abgeschatteten Hauptstrahlungskeule 4 wird zweckmäßigerweise aus der Intensitätsverteilung I(α, β) des Radarstrahls berechnet. Hierbei wird die Intensitätsverteilung über den gesamten Querschnitt, in dem der Radarstrahl eine nicht vernachlässigbar geringe Intensität aufweist, berücksichtigt. Dieser Querschnitt ist größer als der dargestellte 3dB-Querschnitt und abhängig von der Geometrie des Sende-/Empfangsspiegels 3. Zur Illustration ist β_m grob angenähert mittig im verbleibenden Teil der teilweise abgeschatteten Hauptstrahlungskeule 4 dargestellt.
• Der Höhenunterschied Δβ kann für eine Extrapolation aus einem Vertikalprofil der jeweiligen Wetterdaten verwendet werden. Hierzu wird das gemessene Vertikalprofil um den Höhenunterschied Δβ nebst zugehörigem Wetterdatenwert nach unten verlängert, und basierend auf dem verlängerten Vertikalprofil wird eine Extrapolation vertikal nach unten durchgeführt. Gegebenenfalls wird der so erhaltene korrigierte Wert als weiterer Wert W₆ zur gewichteten Kombination mit den Werten W₁ bis W₅ verwendet. Hierbei kann eine Gewichtungsfunktion g₆ verwendet werden, die z.B. mit zunehmendem Grad der Abschattung oder mit zunehmendem Höhenunterschied Δβ abnimmt.
• Das Verfahren ist in einer dem Wetterradar nachgeschalteten Vorrichtung in Form einer Datenverarbeitungsanlage mit einer Dateneingabevorrichtung zum Empfangen der zu korrigierenden Wetterdaten und einer Datenausgabevorrichtung zum Ausgeben der korrigierten Wetterdaten durchführbar. Hierzu kann ein zum Durchführen des Verfahrens erforderliche Instruktionen und Daten enthaltendes Computerprogrammprodukt über einen Datenträger und/oder ein Datennetz in die Vorrichtung geladen werden.

Claims (20)

1. Verfahren zum Korrigieren von Wetterdaten, die aus von einem bodengestützten, radarbasierten Fernerkundungsgerät (1) zur Messung atmosphärischer Bedingungen aufgenommenen, auswertbare Wetter- und störende Bodenechos enthaltenden Radarechodaten gewonnen worden sind, dadurch gekennzeichnet, dass als korrigierte Wetterdaten die aus den Radarechodaten gewonnene differentielle Reflektivität durch Herausrechnen der störenden Bodenechos aus den gewonnen Wetterdaten unter Verwendung einer zuvor gespeicherten Cluttermap mit einer Intensitätsverteilung von Radarechodaten, die die Bodenechos im wesentlichen ohne Wetterechos enthalten, neu berechnet wird.
2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass die korrigierte differentielle Reflektivität ZDR_Wea entsprechend ZDR_Wea = ZDR (Z_CI + Z_Wea) / Z_Wea - ZDR_CI Z_CI / Z_Wea neu berechnet wird, wobei ZDR die aus den Radarechodaten gewonnene differentielle Reflektivität, Z_CI die Intensität aus der Cluttermap, Z_Wea die Intensität des Wettersignals und ZDR_CI die differentielle Reflektifität aus der Cluttermap ist.
3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass die Radialgeschwindigkeit V_Wea entsprechend V_Wea = V (Z_CI + Z_Wea) / Z_Wea neu berechnet wird, wobei V die aus den Radarechodaten gewonnene Radialgeschwindigkeit, Z_CI die Intensität aus der Cluttermap und Z_Wea die Intensität des Wettersignals ist.
4. Verfahren nach einem der Ansprüche 1 bis 3, dadurch gekennzeichnet, dass die spektrale Breite W_Wea entsprechend W_Wea = W (Z_CI + Z_Wea) / Z_Wea neu berechnet wird, wobei W die aus den Radarechodaten gewonnene spektrale Breite, Z_CI die Intensität aus der Cluttermap und Z_Wea die Intensität des Wettersignals ist.
5. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daß die korrigierte Radialgeschwindigkeit, spektrale Breite bzw. differentielle Reflektivität durch gewichtete Mittelwertbildung aus wenigstens zwei unabhängig erhaltenen Werten berechnet wird, wenn die Bedingung Z_Wea < Z_CI erfüllt ist, wobei Z_CI die Intensität aus der Cluttermap und Z_Wea die Intensität des Wettersignals ist.
6. Verfahren nach Anspruch 5, dadurch gekennzeichnet, daß als einer der Werte V_Wea gemäß Anspruch 3, W_Wea gemäß Anspruch 4 bzw. ZDR_Wea gemäß Anspruch 2 verwendet wird.
7. Verfahren nach Anspruch 6, dadurch gekennzeichnet, daß zum Gewichten eine Gewichtungsfunktion verwendet wird, die mit zunehmendem Faktor Z_Wea/Z_CI zunimmt.
8. Verfahren nach einem der Ansprüche 5 bis 7, dadurch gekennzeichnet, daß einer der Werte mittels radialer Interpolation aus von Clutter unbeeinträchtigen Bereichen (E) berechnet wird.
9. Verfahren nach Anspruch 8, dadurch gekennzeichnet, daß zum Gewichten eine Gewichtungsfunktion verwendet wird, die mit zunehmendem Abstand zwischen den nächstgelegenen von Clutter unbeeinträchtigten Bereichen (E) abnimmt.
10. Verfahren nach Anspruch 9, dadurch gekennzeichnet, daß ein azimutaler Abstand und/oder ein radialer Abstand verwendet werden.
11. Verfahren nach einem der Ansprüche 5 bis 10, dadurch gekennzeichnet, daß einer der Werte durch vertikale Extrapolation aus einem Vertikalprofil der Radialgeschwindigkeit, der spektralen Breite bzw. der differentiellen Reflektivität verwendet wird.
12. Verfahren nach Anspruch 11, dadurch gekennzeichnet, daß zum Gewichten eine Gewichtungsfunktion verwendet wird, die mit zunehmendem vertikalen Abstand zum nächstgelegenen von Clutter unbeeinträchtigten Bereich abnimmt.
13. Verfahren nach Anspruch 11 oder 12, dadurch gekennzeichnet, daß zum Gewichten eine Gewichtungsfunktion verwendet wird, die mit zunehmender Zuverlässigkeit des Vertikalprofils zunimmt.
14. Verfahren nach einem der Ansprüche 5 bis 10, dadurch gekennzeichnet, daß die mittlere Höhe (β_m) des verbleibenden Teils eines durch ein Hindernis teilweise abgeschatteten Radarstrahls (4) bestimmt und einer der Werte aus einem Vertikalprofil der jeweiligen aus den Radarechodaten gewonnenen Wetterdaten unter Berücksichtigung der mittleren Höhe (β_m) berechnet wird.
15. Verfahren nach Anspruch 14, dadurch gekennzeichnet, daß die Abschattung berechnet und zum Gewichten eine Gewichtungsfunktion verwendet wird, die mit zunehmender Abschattung abnimmt.
16. Verfahren nach einem der Ansprüche 8 bis 15, dadurch gekennzeichnet, daß als von Clutter unbeeinträchtige Bereiche (E) Bereiche verwendet werden, in denen das Verhältnis aus den Intensitäten der Wetter- und Bodenechos oberhalb eines vorgebbaren Schwellenwerts liegt.
17. Verfahren nach einem der Ansprüche 8 bis 16, dadurch gekennzeichnet, daß als von Clutter unbeeinträchtige Bereiche (E) Bereiche verwendet werden, in denen die Abschattung unterhalb eines vorgebbaren Schwellenwerts liegt.
18. Vorrichtung zum Korrigieren von Wetterdaten, die aus von einem bodengestützten, radarbasierten Fernerkundungsgerät (1) zur Messung atmosphärischer Bedingungen aufgenommenen, auswertbare Wetter- und störende Bodenechos enthaltenden Radarechodaten gewonnen worden sind, dadurch gekennzeichnet, daß die Vorrichtung zum Berechnen korrigierter Wetterdaten der differentiellen Reflektivität durch Herausrechnen der störenden Bodenechos aus den gewonnen Wetterdaten unter Verwendung einer zuvor gespeicherten Cluttermap mit einer Intensitätsverteilung von Radarechodaten, die die Bodenechos im wesentlichen ohne Wetterechos enthalten, ausgestaltet ist.
19. Vorrichtung nach Anspruch 18, dadurch gekennzeichnet, daß sie zum Durchführen des Verfahrens nach einem der Ansprüche 2 bis 17 ausgestaltet ist.
20. Computerprogrammprodukt mit darauf gespeicherten oder zumindest teilweise über ein Datennetz abrufbaren Instruktionen, dadurch gekennzeichnet, daß bei Ausführung auf einer Vorrichtung insbesondere nach Anspruch 18 oder 19 das Verfahren nach einem der Ansprüche 1 bis 17 durchführbar ist.

Images