WO2007022376A2 - Constant altitude plan position indicator display for multiple radars - Google Patents

Constant altitude plan position indicator display for multiple radars Download PDF

Info

Publication number
WO2007022376A2
WO2007022376A2 PCT/US2006/032175 US2006032175W WO2007022376A2 WO 2007022376 A2 WO2007022376 A2 WO 2007022376A2 US 2006032175 W US2006032175 W US 2006032175W WO 2007022376 A2 WO2007022376 A2 WO 2007022376A2
Authority
WO
WIPO (PCT)
Prior art keywords
radar
data
radars
overlapping
method
Prior art date
Application number
PCT/US2006/032175
Other languages
French (fr)
Other versions
WO2007022376A3 (en
Inventor
Pravas R. Mahapatra
Vishnu V. Makkapati
Original Assignee
Honeywell International Inc.
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
Priority to US11/206,696 priority Critical patent/US7821448B2/en
Priority to US11/206,696 priority
Application filed by Honeywell International Inc. filed Critical Honeywell International Inc.
Publication of WO2007022376A2 publication Critical patent/WO2007022376A2/en
Publication of WO2007022376A3 publication Critical patent/WO2007022376A3/en

Links

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/87Combinations of radar systems, e.g. primary radar and secondary radar
    • 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/02Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
    • G01S13/50Systems of measurement based on relative movement of target
    • G01S13/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • G01S13/589Velocity or trajectory determination systems; Sense-of-movement determination systems measuring the velocity vector
    • 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/95Radar or analogous systems specially adapted for specific applications for meteorological use
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/04Display arrangements
    • G01S7/06Cathode-ray tube displays or other two-dimensional or three-dimensional displays
    • G01S7/10Providing two-dimensional and co-ordinated display of distance and direction
    • G01S7/12Plan-position indicators, i.e. P.P.I.
    • 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
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers
    • G01S7/295Means for transforming co-ordinates or for evaluating data, e.g. using computers
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change.
    • Y02A90/12Specially adapted for meteorology, e.g. weather forecasting, climate modelling
    • Y02A90/17Weather surveillance systems using the reflection or reradiation of electromagnetic waves
    • Y02A90/18Radar-based

Abstract

A system and a method to receive radar data from a plurality of radars and to reduce the data from the plurality of radars to a common universal latitude-longitude coordinate frame. The plurality of radars may be positioned at separate geographical locations and may have one or more overlapping regions. The data may be combined in the one or more overlapping regions.

Description

CONSTANT ALTITlJDE PLAN POSITION INDICATOR DISPLAY FOR MULTIPLE RADARS

Cross Reference to Related Applications

This application claims priority to India Patent Application No. 151/KOL/2005, filed March 10, 2005, which is incorporated herein by reference.

Background

[O001] Constant Altitude Plan Position Indicator (CAPPl) is a form of data presentation in weather radars. For CAPPI scanning, a horizontal slice is taken through the radar volume scan data at a constant altitude above the earth surface. The radar volume scan data is extracted from full 3-D scan data, and is converted to a 2-D format for presentation in polar coordinates on a computer display, paper printout, or any other two-dimensional surface.

[0002] CAPPI is a variant of the more general Plan Position Indicator (PPI) that displays weather parameters along radials from the radar as function of an azimuth scan angle. A radar antenna transmits and receives pulses at different elevation angles φ and at different azimuth angles θ by both performing a rotating scan operation in the azimuth and by varying the elevation angle. PPI data is generated and recorded by scanning a beam circularly at a constant elevation angle. A volume scan consists of multiple constant-elevation azimuth scans. PPI volume scan data at multiple elevation angles is used to produce CAPPL

[0003] Near the radar site there is often ground clutter, which may interfere with obtaining a clean display of weather. In the beam position(s) with low elevation angle(s), clutter is often so strong that filtering the ground clutter also removes weather signals resulting in gaps in the weather display.

[0004] In PPI scanning, the radar beam may overshoot precipitation altitudes, for a part of the radial, and thereby not detect any precipitation at the corresponding ranges (i.e. distances from the radar). [0005] Weather radar systems often deliberately degrade the time resolution for observation in order to improve signal quality and also to reduce the data handling specifications over long observation periods. The PPI volume scanning mode also degrades the spatial resolution by skipping certain elevation angles to reduce the time for scanning the region around the radar.

[0006] The maximum elevation of scan may be limited to a value less than 90° (i.e. vertical pointing), leaving a conical 'blind zone' over the radar location. This causes a circular hole to appear in the CAPPI, the hole being larger at higher altitudes. [0007] Individual radars may be limited in their range coverage. To get a weather picture over a large geographical area, data may be combined from multiple radars that are spatially separated.

[0008] Combining CAPPI data from multiple radars may pose technical challenges. When a geographical area is covered by many radars, the coverage pattern is not uniform. Certain areas may not be covered at all (i.e. fall between coverage zones of individual radars), certain areas may have coverage only from a single radar, certain areas may have coverage from two radars, certain areas may be covered by three radars, and certain other areas may receive coverage from more than three radars. Such variability of coverage poses technical challenges in data combining.

Brief Description of the Drawings

[0009] Figure 1 illustrates a schematic radar system according to an example embodiment.

[0010] Figure 2 illustrates a schematic cross-section of Earth and a surface at a constant altitude above the earth's surface at which Constant Altitude Plan Position Indicator (CAPPI) data is generated according to an example embodiment. [0011] Figure 3 illustrates a schematic partial cross-section of CAPPI geometry in a vertical section above a plane tangent to the Earth at the radar location, as employed in an example embodiment.

[0012] Figure 4 illustrates a schematic plan view of an area around the radar showing a constant-range circle over which CAPPI data is generated according to an example embodiment.

[0013] Figure 5 illustrates a latitude-longitude coordinate system according to an embodiment.

[0014] Figure 6 illustrates a radar mosaic system at which Constant Altitude Plan Position Indicator (CAPPI) data is generated for overlapping radars according to an example embodiment.

[0015] Figure 7 illustrates a process of generating a mechanism, such as look up tables, according to an example embodiment.

[0016] Figures 8 A, 8B and 8C illustrate a process of generating CAPPI data using the mechanism generated in Figure 7 according to an embodiment. [0017] Figures 9 A and 9B illustrate a process of extracting scalar data according to an embodiment.

[0018] Figure 1 OA and 1 OB illustrate a process of extracting wind vector data according to an embodiment.

[0019] Figure 11 shows a diagrammatic representation of a machine in the form of a computer system within which a set of instructions, for causing the machine to perform any one or more of the methodologies discussed herein, may be executed, according to an example embodiment. Detailed Description CAPPI Data Generation

[0020] Figure 1 illustrates a schematic radar system 100 according to an example embodiment. The radar system 100 includes an antenna 101 coupled to a radar receiver 102, a processor 103 coupled to the radar receiver 102, and a Plan Position Indicator (PPI) display 104 coupled with the processor 103. The processor 103 receives volume scan spherical coordinate data from the radar receiver 102, converts the radar data to rectangular coordinate data, and provides the rectangular coordinate data to the display 104. Additionally or alternatively, the processor 103 may receive the volume scan data recorded from a previous scan in a computer memory 105.

[0021] Figure 2 illustrates a schematic cross-section 200 of Earth 201 and a surface 202 having a constant altitude H above the earth's surface at which Constant Altitude Plan Position Indicator (CAPPI) data is generated according to an example embodiment. An origin of the radar antenna 101 is represented at point D on the earth's surface. The center of the Earth is represented by O. The Earth's radius is represented by Re. Point C represents the point at which the CAPPI data may be generated at the altitude H. The angle at the earth's center between point D and point C is represented by β. The point E is directly above point D at the altitude H.

[0022] An elevation angle φ represents the angle at point D between the earth's surface (i.e. a tangent to the Earth's surface at point D) and the line CD. The specific elevation angle φ may have been skipped during the radar scan operation. In an example embodiment, gaps such as this may be filled using an interpolation scheme to potentially provide spatially continuous information of weather at the given altitude H. The elevation angle φ of the radar is computed for each increment in EC using equation set (2) and Figure 2. n EC β~ Re ± H CI= (Re + H)sin(β) OI= (Re + H)cos(β) ID= (R + H)cos(β)- r%

CD= ^CI2 + ID2 [0023] Because the Earth is curved and the scan elevation interval may be between a minimum value and a maximum value, a point such as point C in Figure 2 on the CAPPI surface may have an elevation angle φ outside the scanned elevation interval (or in between two scanned elevation angles). For altitudes on the CAPPI surface that are below the minimum scan elevation, the data may be taken from the lowest elevation angle available. Conversely, if a selected point C is above the maximum scan elevation, the data may be taken from the highest elevation angle available. [0024] Embodiments may account for bending of the radar beam. The radar beam may bend as it passes through layers of air with different refractive indices. Under standard atmospheric conditions, the bending of the radar beam has a radius of curvature about four times the radius of the Earth. Thus, under normal conditions, a radar beam emitted horizontally and at other elevations would take paths that curve slightly below straight line paths.

[0025] Figure 3 illustrates a schematic partial cross-section of CAPPI geometry 220 in a vertical section above a plane tangent to the Earth at the radar location D, as employed in an example embodiment. The horizontal line through D is the tangent to the Earth's surface at point D and the arc 225 represents the spherical surface, at a height H above the radar location D, at which Constant Altitude Plan Position Indicator (CAPPI) data is generated according to an additional example embodiment. [0026] In a volume scanning mode, an elevation angle φ of the antenna 101 is changed incrementally by a determined angle and a horizontally rotating scanning operation is performed along each incremented elevation angle. CAPPI data for point C may not be readily available in the volume scan data, and may be constructed from radar data gathered along other elevation angles and/or azimuth angles, in embodiments of the present invention. The CAPPI data at point C may be generated by interpolating the gathered radar data in elevation. Elevation angles φm and φm.i denote the elevation angles of the radar scan that are closest to angle φ, on either side, as shown in Figure 3. The weather parameter (e.g. reflectivity Z) at the point C, having elevation φ, may be calculated as the weighted average of radar resolution cell data from points A and B, having elevations §m.ι and φm, respectively, in an embodiment. [0027] In Figure 3, the distance from point D to point B (DB) may be represented as RB, the distance from point D to point A (DA) may be represented as RA, and the distance from point D to point C (CD) maybe represented as R0. Ih embodiments, the CAPPI construction includes two processes. The first process includes a mechanism to generate CAPPI scan line data of a weather parameter Z. For a specified CAPPI altitude, the two bounding radials φm-i and φm (i.e. the value of m), the interpolation weight between them, and the ranges CD, DB and DA depend on the length of the arc EC alone, and are independent of the azimuth angle θ.

[0028] In the example embodiments, the CAPPI scan line data is generated by keeping EC constant and varying the azimuth angle θ of the radar from 0 to 2π (or a certain θπ»n to θmax for a sectoral CAPPI) clockwise or counter-clockwise depending on the direction of the radar scan. As discussed in more detail with regard to the process of Figures 8A, 8B, SC, after completing the circular arc at a constant arc length EC, the arc length EC maybe incremented by 1 km (or any other determined step, increment or unit) and the azimuth angle θ may varied again from 0 to 2π (or θmin to θmax)- [0029] In the second step of the CAPPI construction, the CAPPI scan line data is presented to the PPI coordinate conversion, formatting and display mechanism. [0030] Figure 4 illustrates a schematic plan view 230 of an area around the radar showing a horizontal projection of a constant-range circle 235 over which CAPPI data is generated according to an example embodiment. The circle 235 includes a radius of (Rc)(COS φ). This radius is the horizontal component of the range shown as the line segment CD in Figure 3. The CAPPI scan data may display weather conditions occurring over any given constant-height surface within a half-spherical observation volume of the radar. The CAPPI data may be generated about the radar origin D over azimuth angle θ from 0° ≤ θ ≤ 360°. The azimuth angle θ may be incremented by 1 ° or by any other angular increment as specified in the radar data in the CAPPI generation. Latitude-Longitude Coordinate System

[0031] Figure 5 illustrates a latitude-longitude coordinate system 240 according to an embodiment. The latitude-longitude coordinate system offers a grid for mosaicing

CAPPI data from overlapping radars at various altitudes.

10032] The CAPPI radar data may be converted from the individual radar coordinates to a universal grid system as follows. Let Re denote the radius of the earth and if be the altitude above the Earth's surface 245 at which the mosaic is to be constructed.

Coordinate (Xu φi) may be the latitude-longitude location of the radar i on the earth's surface 245.

[0033] Coordinates

Figure imgf000008_0001
)) of the Minimum Bounding Rectangle

(MBR) of radar coverage (e.g. the bounds of the coverage of radar i) may be calculated as

Figure imgf000008_0002

An^ = A1 -AA; A1 +AA

Φmin, = Φl -ΔΦl.'Φmx,

Figure imgf000008_0003
+M

[0035] The coordinates ((A1n^ Jmin),(Amιamax)) ofthe MBRforΛTradars maybe calculated as arguments of the minimum and maximum values, respectively, where i — 1, 2, ... N.

Figure imgf000008_0004

[0036] ^ - -*~««^ (8)

^« = argmin,.=ti2 Nφmhtι

[0037] A distance d between two latitude-longitude points ( A1 , φx ) and (A22) may be calculated as

[0038] d((Aιι),(A22))= Re cos-ι(cos(Aι)cos(λ2)cos(φι - φ2)+ Sm(A1)Sm(A2)) (9) Mosaic for Multiple Radars

[0039] Figure 6 illustrates a radar mosaic 250 at which Constant Altitude Plan Position Indicator (CAPPI) data is generated for overlapping radars according to an example embodiment. The radar mosaic 250 includes overlapping radars I, π, II and IV in this illustrated embodiment. There may be any number N of radars associated with the mosaic, and may not be limited to 4 or may be less than 4.

[0040] In this embodiment, there are 13 separate regions in the mosaic 250. The regions along with the overlapping radars are shown in Table I. Table I is indexed by region and shows the corresponding radars in that particular region. The region may or may not have overlapping radars. Ih the instance where there are at least two radars covering the region, there are overlapping radars. In the regions covered by three radars (e.g. regions 4, 5, 8, 9) or more than three radars (e.g. region 7), a three dimensional wind vector may be constructed. A two dimensional wind vector may be constructed for a region with two overlapping radars (e.g. regions 2, 6, 10, 12).

TABLE I

Figure imgf000009_0001
Flowcharts

[0041] Figure 7 illustrates a process 300 of generating a mechanism, such as look up tables, according to an example embodiment. The process 300 includes the first of the two processes discussed with regard to Figure 3.

[0042] At block 305 of Figure 7, the process begins. At block 310, a value of M is selected such that M-I equals the number of elevation scans performed by the radar 102.

In embodiments of the process, set φo equal to zero, set ΦM equal to 90°, and the angles φi through ΦM-I equal respectively to the radar scan elevation angles from the lowest to the highest along the determined increments.

[0043] At block 320, an arc length EC (φm) may be calculated for each elevation angle φra of the radar. The interval for the elevation number m may be: 0^n<M. The elevation number m may be incremented, by 1, for example, in the interval. The arc length EC may be calculated for each elevation angle φm with the formula:

EC (φm)

Figure imgf000010_0001
- φπj.

[0044] At block 330, a first mechanism of EC(φm) values in terms of the elevation number m of the radar may be generated using the values generated at block 320. The first mechanism may be a look up table indexed by elevation number m and/or elevation angle φm, a graph, an algorithm, a chart, and/or any other possible mechanism.

[0045] At block 335, the process of generating the first mechanism ends.

[0046] Figures 8A, 8B, and 8C illustrate a process 400 of generating CAPPI data using the mechanisms generated in Figure 5 according to an example embodiment.

[0047] At block 401, the process 400 of Figure 8A begins.

[0048] At block 405, the process 400 sets the value of the arc length EC, illustrated in

Figure 3, to zero. The value of M is set to one more than the number of total scan elevation angles φ in the radar volume scan data set.

[0049] At block 410, the length EC is incremented by one step, which may be one (1) km or any other chosen value.

[0050] At block 415, the range CD and the corresponding angle φ may be calculated using equation (2). The range AC may be calculated using the equation: AC = CD x sin(φ - ΦM-I)/(COS(ΦM-I + β). The range DA may be calculated by solving for DA in the equation: (DA)2 cos2M<) + (DA) SUn(^-1) = H -(AC) . Optionally, for a given

radar and scan cycle, instead of calculating the range DA, the values of range DA may be read from a pre-stored mechanism, such as apre-stored look-up table LUTR, of range

DA indexed by arc length EC. The pre-stored mechanism may be from a previous computation of the range DA for the same CAPPI altitude.

[0051] The azimuth angle θ may be incremented from θmin to θmax in steps of Δθ, where Δθ may be any degree, such as 1 degree. The radar data ZA may be retrieved at range DA, from the original scan data at elevation ΦM-I, for all θ values. The radar data

ZA may be stored in a CAPPI scan line data buffer indexed by arc length EC and azimuth angle θ.

[0052] At block 420, the process 400 is queried as to whether EC is less than or equal to EC(ΦM-I). If the answer to the query is 'yes' then the process may proceed to block

410. If the answer to the query is 'no', then the process may proceed to block 425.

[0053] At block 425, m may be initialized to M-I and the process may proceed to block 430 in Figure 8B.

[0054] At block 430, the weights BC and AC may be calculated using the equations

BC = CD x sin(φm - φ)/(cos(φm + β) and AC = CD x sin(φ - φm-iyCcosCφm-i + β), where

CD and φ are obtained using equation (2). Optionally, for a given radar and scan cycle, the values of weights BC and AC may be read from a pre-stored mechanism, such as a pre-stored look-up table LUTW, from a previous computation of BC and AC for the same CAPPI altitude.

[0055J At block 435, determine ranges DB and DA solving for DB and DA in the equations: ^-(DB)2 cos2m) + (DB)sm(φm) = H + (BC) and

— — (DA)2 cos2m.i) + (DA)sva(φm_ι) = H - (AC) , respectively. Optionally, for a

given radar and scan cycle, the values of ranges DB and DA may be read from a pre- stored mechanism, such as the look-up table LUTR indexed by arc length EC, &om a previous computation of ranges DB and DA for the same CAPPI altitude. [0056] At block 440, CAPPI data is calculated at point C. Weather data ZA may be retrieved from the radar volume scan data for elevation φm.j and range DA. Weather data ZB may be retrieved for elevation φm and range DB from the radar volume scan data. The weather data Zc at CAPPI data point C is the interpolation of ZA and ZB with weights BC

((BCX Z M (ACX Z11)) and AC, and may be calculated according to the equation Zc jz— — — .

[0057] At block 445, the length EC is incremented by one step.

[0058] At block 450, the process is queried as to whether EC is less than or equal to

ECCφm-O. If the answer to the query is 'yes' then the process may proceed to block 430.

If the answer to the query is 'no', then the process may proceed to block 455.

[0059] At block 455, m is decremented by 1.

[0060] At block 460, the process is queried as to whether m is equal to 1. If the answer to the query is 'no' then the process may proceed to block 430. If the answer to the query is 'yes', then the process proceeds to block 465. The process may proceed to block 465 in Figure 8C.

[0061] At block 465, for a given arc length EC, calculate weight BC using the equation BC = CD x sin(φi - φ)/(cos(φi + β), then solve for range DB using the equation: (DBf cos2 tø, ) + (DB) sintø, ) = H + (BC) . Angle θ may be incremented from θmin

to θmax in steps of Δθ. The radar data ZB at range DB may be retrieved from the scan data at elevation φi, for all θ values. The radar data ZB may be stored in CAPPI scan line data buffer indexed by length EC and angle θ. Optionally, for a given radar and scan cycle, the values of DB may be read from the pre-stored mechanism, such as the pre-stored look-up table LUTR indexed by EC, from a previous computation of range DB for the same CAPPI altitude.

[0062] At block 470, the length EC is incremented by one step.

[0063] At block 475, the process is queried as to whether EC is less than or equal to

EC(φo) at φo = 0°. If the answer to the query is 'yes' then the process proceeds to block

465. If the answer to the query is 'no', then the process proceeds to block 480.

[0064] At block 480, CAPPI scan line data generation ends. [0065] At block 485 , PPI display scan conversion occurs, as described herein. PPI scan converted CAPPI scan line data maybe obtained and stored. [0066] At block 490, the generated CAPPI radar data Zc along the surface 202 at the altitude H above the earth's surface may be displayed on the display 104 and/or printed.

Superposition of Scalar Data

[0067] Figures 9A and 9B illustrate a process 600 of extracting scalar data according to an embodiment. The overlapping radars at a point, as shown in Figure 6, may be determined dynamically. For scalar data, the average or maximum of the values for the overlapping radars at the point may be calculated and displayed at that point.

[0068] The process 600 for constructing the mosaic of scalar data with data from N radars may begin at block 605.

[0069] At block 610 the Minimum Bounding Rectangle (MBR) coordinates

Figure imgf000013_0001
max, )) may be determined for radar i (l ≤ i ≤N) using equation (7).

[0070] At block 620, coordinates associated with the MBR for N radars

((λmk'ΦmtΛ&mox'Φmx)) may be determined using equation (8).

[0071 ] At block 630, grid resolution parameters (dimensions of latitude-longitude cell) δλ and δφ, may be set to any specified value. The latitude-longitude cell of dimension δλ and δφ at location (λmjnmm) is chosen as the first of the latitude-longitude cells x within the MBR for N radars, obtained in block 620.

[0072] At block 640, a sum s of individual CAPPI data at x, from all weather radars observing cell x, is set to zero. A number of overlapping radars n is set to zero. The radar index i is set to 1. The index of cell x is set to 1.

[0073] At block 650, the process 600 queries as to whether designated conditions are met for radar i. If the answer to the query is no, the process proceeds to block 680 with s and n held at their current values. If the answer to the query is yes, the process proceeds to block 670.

[0074] For the latitude-longitude cell x at coordinates (λ, φ), the designated conditions include first checking if the inequality An^1 ≤ λ ≤ X^1 is true, if yes then checking if the inequality φ^ ≤ φ ≤ φnm/ is satisfied, and if yes, then further checking if the inequality d((λ, φ),{λt,φ^< EC1^f) is satisfied, where ECj(O0) maybe retrieved from the first mechanism of CAPPI for the ift radar generated at block 330 in Figure 7.

[0075] At block 670, the number of overlapping radars n is set to n+1. The sum of weather data s is set to (s + PPI scan-converted CAPPI scan line data from the 1th radar at the cell x). The PPI scan-converted CAPPI scan line data from the ift radar at the cell x may be obtained by first computing the east-west distance cήawi and north-south distance dnsi using the formulae dm\ = <*(0W)>(Λ>&)) and ^NSi = d((λ,^.),(^,^)) respectively, of x from the ith radar location. The PPI scan converted data obtained and stored at block 485 of Figure 8C, corresponding to the coordinates closest to dewi and

O4 NSi, may be read. Ih case of reflectivity data, the data is converted from a logarithmic dBZ scale to a linear scale Z given by Z = 10dBZ/10.

[0076] At block 680, the process 600 is queried as to whether radar index i is less than or equal to the number of radars N. If the answer to the query is yes, the process proceeds to block 690. If the answer to the query is no, the process proceeds to block

700.

[0077] At block 690, the radar index / may be set to i+ 1. The process then proceeds to block 650 such that determinations may be made for each radar i.

[0078] At block 700, the process is queried as to whether the number of overlapping radars n is equal to zero. If the answer to the query is no, the process proceeds to block

720. If the answer to the query is yes, the process proceeds to block 730.

[0079] At block 720, weather data at cell x having coordinates (λ, <j>), may be set to the sum s of individual radar CAPPI data, divided by the number of overlapping radars n.

The process may proceed to block 735.

[0080] At block 730, there are no overlapping radars at cell x (λ, <j>), and the weather/mosaic data at cell x (λ, φ), may be set to zero. ( The process may proceed to block 735.

[0081] At block 735, weather data for cell x, at (λ, φ), may displayed against the coordinates λ and φ. In case of reflectivity data, the data is converted from the linear Z scale to the logarithmic dBZ scale by using the formula dBZ = 10 logto (Z). The process may proceed to block 740. [0082] At block 740, the process is queried as to whether each cell x of the MBR for the N radars has been determined. If the answer to the query is no, the process proceeds to block 745. If the answer to the query is yes, the process proceeds to block 750. [0083] At block 745, the cell x may be set to cell x+1 , and the process proceeds to block 640. [0084] At block 750, the process ends.

Superposition of Vector Data

[0085] Figures 1OA and 1OB illustrate a process 800 of extracting wind vector data according to an embodiment. With vector data, a wind velocity vector may be constructed from the values observed by the overlapping radars. The wind velocity may be extracted in two or three dimensions, depending on the combination of radars for a given element in space. For example, if there are two overlapping radars, the velocity vector may be extracted in two dimensions, and in three dimensions for three or more overlapping radars at that given element in space. The vector is usually represented by an arrow, with its head pointing towards the direction of the wind vector and the length proportional to the magnitude of the wind vector. The vector may not be displayed at every point if coverage by two or more radars is not available at every point. A single vector may be displayed for a region, since displaying an arrow for each pixel on a display would make the display very crowded. For convenience and neatness of display, each vector may be displayed over a square or rectangular block region (called a display block) as discussed in more detail below.

[0086] In the process 800 described in Figure 1OA, let kxp be the size of the rectangular display block chosen for constructing vector data. If k is chosen to be equal to p then the display block may be square. The process to extract the wind vectors with data fromiVradars begins at block 810.

[0087] At block 820, the Minimum Bounding Rectangle (MBR) coordinates (UΛ, ' Φmh, ) (Am*, - Φmax, )) ∞ay be determined for radar i (1 < i ≤ N) using equation (7). [0088] At block 830, coordinates ((λm!nmin),(λmaxmax)) of the MBR for N radars may be determined using equation (8). [0089] At block 840, the latitude and longitude interval between the coordinates {(Km > Φmm )> {Kox • Φmax )) may be divided into "D" rectangular display blocks of size k x p, i.e. there are "D" rectangular display blocks in the grid. Further, the grid resolution parameters (dimensions of latitude-longitude cell) δλ and δφ, may be set to any specified value.

[0090] At block 850, an index B of the current rectangular display block may be set to 0. At block 860, set B = 1 + B.

[0091] At block 870, the largest area within the display block B that is simultaneously observed by a given number of radars M may be determined. The wind vector for the display block B in Figure 6 may be determined. However, a display block may contain multiple regions, as illustrated by the display block B in Figure 6, which includes regions 7, 8, 9 and 12. These regions are determined by considering latitude- longitude cells of dimension δλ and δφ within the display block and determining the radar(s) which can observe each cell by applying the designated conditions. For a cell located at coordinates (λ, φ), the designated conditions include first checking if the inequality Xn^1 ≤ λ ≤ A10^1 is true, if yes then checking if the inequality φtώXs ≤φ ≤ φm&ι is satisfied, and if yes, then further checking if the inequality d((λ, φ),(Λ., φf)) < ECt(θ°) is satisfied, where ECj(O0) may be retrieved from the first mechanism of CAPPI for the iΛ radar generated at block 330 in Figure 7. The region with maximum area within that display block may be determined by counting the number of latitude-longitude cells that are observed by a given combination of radars, e.g. region 12 of Figure 6 is observed by two (2) overlapping radars. The wind vector for the entire display block B may be constructed using the data from radars corresponding to that region, e.g. a 2-dimensional wind vector.

[0092] At block 880, the PPI scan-converted CAPPI scan line data from the 1th radar for the cell located at (λ,φ) is obtained by first computing the east-west distance dEWi and north-south distance dNsi from the formulae dEwi - ^(0WXOW/)) 3^ dNsi = d((λ>Φi),(λι,Φi)) respectively, of x from the ith radar location. The PPI scan converted Doppler wind velocity data obtained and stored at block 485 of Figure 8C, corresponding to the coordinates closest to dswi and d^si, may be read. The average A(i) or At of the Doppler wind velocity values for each radar i corresponding to the largest region within B may be calculated independently. For example, the average of the Doppler wind velocity values for each radar (i.e. radars III and IV) covering the part of region 12 enclosed in display block B of Figure 6 may be calculated.

[0093] At block 900, the process 800 is queried as to whether the overlapping radars M may be equal to either 2 or 3. If the answer to the query is yes, the process proceeds to block 910. If the answer to the question is no, the process proceeds to block 920. [0094] At block 910, the average wind velocity V may be obtained by solving the dot-product equation V ■ P(i) = A(i), where P(i) is the unit position vector from the origin of the radar i to the center of the block B, with i taking on the values corresponding to the radars covering the largest region within B. The process may proceed to block 950. [0095] At block 920, the process is queried as to whether the overlapping radars M is less than 2. If the answer to the query is yes, the process proceeds to block 930. If the answer to the question is no, the process proceeds to block 940. [0096] At block 930, there is not enough data to determine a 2-dimensional or 3- dimensional wind vector at block B. The vector data for this block B is not displayed. [0097] At block 940, the average wind velocity V may be obtained by solving the dot-product equation V • P(i) = A(i) using the least squares method, where P(i) is the unit position vector from the origin of the radar i to the center of the block B, with i taking on the values corresponding to the radars covering the largest region within display block B. The process then proceeds to block 950.

[0098] At block 950, the wind vector is displayed as an arrow with its center located at the latitude and longitude coordinates corresponding to the center of the display block B. The process may proceed to block 960.

[0099] At block 960, the process is queried as to whether the wind vector for each block B has been determined. If the answer to the query is yes, the process proceeds to block 970. If the answer to the question is no, the process proceeds to block 860. [00100] At block 970, the process ends.

[00101] The wind velocity may be extracted using radial velocities from two or more separately located radars scanning overlapping areas. The 3-dimensional velocity components of the wind may be obtained from the combination of three Doppler radars. [00102] The position vector from the origin of the radar to the center of the block B may be calculated for the radars corresponding to the largest coverage region within display block B. The average of the values may be calculated independently for each radar covering that region. The dot product of the unit position vector and wind vector may .give the radial component of velocity. However, because the vector is displayed over a display block, the radial component may be taken as the average of the values. These equations may be solved to obtain the actual vector in two or three dimensions depending on the number of overlapping radars. If the set of equations is written in matrix notation (equation (10)), then the least squares solution may be found by solving equation (11) where Aτ denotes the transpose of A. [001031 AX = C (10)

[00104] A7AX = ATC (11)

[00105] In equations (10) and (11) X may represent a column vector for the wind velocity at display block B, A may be a matrix whose rows represent the unit position vectors of the radars observing display block B (i.e. the largest area within display block B), and C may be a column vector of the average values of the Doppler velocities for each radar over the display block B (i.e. the largest area within display block B).

Combined Scalar Data and Vector. Data from Multiple Radars [00106] Grids of multiple weather radars may generate data on scalar (reflectivity, spectrum width) as well as vector (wind velocity) parameters of the atmosphere. Because of the curvature of the earth's surface, the radar grid may not be co-planar. Scalar and vector data from multiple radars may be combined in a grid and displayed.

Computer System

[00107] Figure 11 illustrates a diagrammatic representation of a machine in the example form of a computer system 1001 within which a set of instructions, for causing the machine to perform any one or more of the methodologies discussed herein, may be executed. In alternative embodiments, the machine operates as a standalone device or may be connected (e.g., networked) to other machines. In a networked deployment, the machine may operate in the capacity of a server or a client machine in server-client network environment, or as a peer machine in a peer-to-peer (or distributed) network environment. The machine may be a server computer, a client computer, a personal computer (PC), a tablet PC, a set-top box (STB), a Personal Digital Assistant (PDA), a cellular telephone, a mobile device, a palmtop computer, a laptop computer, a desktop computer, a control system, or any machine capable of executing a set of instructions (sequential or otherwise) that specify actions to be taken by that machine. [00108] Further, while a single machine is illustrated, the term "machine" shall also be taken to include any collection of machines that individually or jointly execute a set (or multiple sets) of instructions to perform any one or more of the methodologies discussed herein. .

[00109] The example computer system 1001 includes the processor 103/1002 (e.g., a central processing unit (CPU), a graphics processing unit (GPU), or both), a main memory 1004 and a static memory 1006, which communicate with each other via a bus 1008. The computer system 1001 may further include a video display unit 1010 (e.g., a liquid crystal display (LCD) or a cathode ray tube (CRT)). The computer system 1001 also includes an input device 1012 (e.g., a keyboard), a cursor control device 1014 (e.g., a mouse), a disk drive unit 1016, a signal generation device 1018 (e.g., a speaker) and a network interface device 1020.

[00110] The disk drive unit 1016 includes a machine-readable medium 1022 on which is stored one or more sets of instructions (e.g., software 1024) embodying any one or more of the methodologies or functions described herein. Thainstructions 1024 may also reside, completely or at least partially, within the main memory 1004, the static memory 1006, and/or within the processor 1002 during execution thereof by the computer system 1001. The main memory 1004 and the processor 103/1002 also may constitute machine- readable media. In an additional embodiment, the mechanism (such as the look up table) is not stored, but rather a processor or additional processor is used to generate the weights substantially in real-time. This additional embodiment may be useful, e.g. where processing speed is more readily available as compared with memory. [00111] The term "machine-readable medium" shall also be taken to include any medium that is capable of storing, encoding or carrying a set of instructions for execution by the machine and that cause the machine to perform any one or more of the methodologies of the present invention. The term "machine-readable medium" shall accordingly be taken to include, but not be limited to, solid-state memories, optical and magnetic media, and carrier wave signals.

[00112] The instructions 1024 may further be transmitted or received over a network 1026 via the network interface device 1020.

[00113] Embodiments may be utilized by running the program on a digital computer that receives the appropriate type of data (scalar and vector) from a number of overlapping radar fields to generate a CAPPI display on the monitor in quasi-real-time. Embodiments are also useful for generating CAPPI display from pre-recorded radar data available from public or private data archives. The display unit 1010 may be of any resolution and embodiments maybe implemented without addition of any special hardware to a computer.

[00114] The mechanism(s) may also be implemented on a Digital Signal Processing (DSP) chip or any other computer board. The mechanism(s) is(are) implemented in a high level programming language for ease of coding, though they may also be implemented in other types of programming languages, e.g. in assembly or machine languages to achieve higher processing speed and reduced memory overheads. Due to the versatile nature of the mechanism, the technique may be embedded in hardware for CAPPI display of data from radar receiver during real-time operation. [00115] Applications that may include the apparatus and systems of various embodiments broadly include a variety of electronic and computer systems. Some embodiments implement functions in two or more specific interconnected hardware modules or devices with related control and data signals communicated between and through the modules, or as portions of an application-specific integrated circuit. Thus, the example system is applicable to software, firmware, and hardware implementations. [00116] Using the embodiments of methods and systems described herein, a horizontal slice is taken through the radar volume scan data at constant altitude above the earth surface using Constant Altitude Plan Position Indicator (CAPPI). Weather radar systems often deteriorate the spatial sampling of the scan volume by skipping certain elevation angles to reduce time used for scanning the region around the radar. This results in gaps in the CAPPI data, which have been filled in the embodiment using an interpolation scheme. Other CAPPI schemes may utilize a nearest neighbor approach, which may introduce artifacts into the display, including sharp jumps or boundaries between annular segments of the display. Embodiments determine the display parameter at the CAPPI surface by making a weighted interpolation of data from adjacent scan elevations to obtain a smoother and more accurate CAPPI display. The CAPPI scan line data constructed along each azimuth angle may be scan converted to computer display coordinates.

[001171 In systems described herein, the elevation from the ground of the CAPPI surface maybe first specified, and CAPPI data for increments in EC between the vertical and the highest elevation angle of radar is taken from the radar volume scan data for the highest elevation angle. The CAPPI data for increments in EC between the lowest elevation angle and 0° (i.e. horizontal) is taken from the radar volume scan data for the lowest elevation angle. The length of the arc EC for each elevation angle is computed. These pre-computed arc lengths may be used for determining the adjacent elevation angles corresponding to each increment in EC; the scheme may avoid computation of adjacent elevation angles for each increment in EC. Further, interpolation may add little overhead due to calculation of weights and weighted average of weather data from adjacent elevation angles.

[00118] Coverage by multiple radars may help in improving the quality of scalar radar data such as reflectivity, which corresponds to rainfall intensity. Multiple Doppler radar coverage may also help in retrieving the wind vectors which may not be observed directly by single Doppler radars. The multi-radar data combining method or system may handle scalar data as well as vector data.

[00119] Any number of radar fields overlapping at a given point may be handled using the method(s) described herein. The number of overlapping radars at a given element may be determined dynamically using a distance measure bounded by a value depending on the altitude of the element from the surface of the Earth. For scalar data, the data from the overlapping radars may be averaged and displayed on a CAPPI (Constant Altitude Plan Position Indicator) projection. For vector data, the projections corresponding to the individual radar radial directions may be used to reconstruct a true vector at the radar resolution cell. The vector data may be averaged over a display box and may be represented by an arrow for each box. For two overlapping radars, the vector is calculated as a planar (2-D) vector. With 3 overlapping radars, an exact solution of the 3- D vector is obtained. If the overlapping radars are more than 3 then a 3-D vector is reconstructed on a least-squares basis. The software may automatically determine the number of overlapping radars and employ the appropriate reconstruction algorithm. [00120] Individual radar data may be generated and archived in a radar-centric coordinate system. When combining data from multiple radars, the data may be converted to a universal coordinate system that may increase functionality and efficiency. In embodiments, a universal latitude-longitude grid for displaying the data generated at constant altitude maybe used.

[00121] In embodiments described herein, a mechanism, such as a graph, a look-up table (LUT), a set of equations, and/or an algorithm may be used for calculating or storing the arc lengths EC for each radar beam elevation φ, and weights and ranges for each increment in EC for a given altitude H from ground. For ease of explanation, the mechanism is a look up table (LUT) in these embodiments. The LUT for ranges and weights (i.e. the second and third mechanisms) is implemented as one dimensional arrays' indexed by the number of the increments in EC from the vertical, while the LUT for lengths of EC for each elevation angle is indexed by the elevation number. The LUT can be computed with minimal processing overhead and can be accessed quickly. The technique may be efficient because the memory used for storing the LUTs may be minimal and an efficient calculation of the lengths may add less processing overhead. This technique may not recompute the LUT unless the altitude or the elevation angles of the radar data change. Since the elevation angles of the radar data are less likely to change, multiple LUTs may be computed and stored for different altitudes. The use of LUTs enables high performance with minimal memory overhead and memory may be less expensive as compared to the increase of the processing speed of hardware. It also enables the CAPPI data computation to be carried on processors of lower capability such as those in airborne computers.

[00122] CAPPI may avoid issues associated with PPI by picking constant altitude data from different elevation scans. However, because CAPPI includes data from all elevation scans in the radar volume, processing in some embodiments may take considerably longer than PPI displays.

[00123] As described herein, the CAPPI data generation process described may be compatible with any PPI display algorithm or device that may handle any user-specified rotation, zooming, magnification, distance interval, and sector selection. This scheme may provide additional features for displaying CAPPI. The CAPPI data generated using these embodiments are indexed by the horizontal distance (along constant-height arc) from the vertical through the radar, and the azimuth angle, and may be converted to any other coordinate systems, e.g. latitude-longitude, Cartesian. Embodiments are suitable for generating CAPPI data for use with techniques superposing CAPPI data from overlapping radars, where the radar data should be converted from the individual radar coordinates to a universal grid system. CAPPI, generated using this interpolation scheme, suits viewing specifications and may be more accurate and continuous than data generated using another scheme, for example, a nearest neighbor scheme. [00124] Embodiments described herein may be used with the presentation and display of weather radar data by users such as meteorologists, air traffic controllers, pilots, TV weather broadcasters, and disaster monitors. Embodiments can be licensed to companies working on general radars, weather radars, imaging radars, meteorological data products and PPI displays.

[00125] Overlapping radar sites may increase reliability and may fill the void of non- operational sites. Such an overlapping network may offer a superior quality of coverage and the observed information may be more precise. For example, a single radar may have a blind azimuth that an adjacent radar may cover. Multiple radars may provide continuity and redundancy in case of a failure. Where a single radar surveillance system may experience catastrophic failure, multiple radar systems may provide for graceful degradation of overall performance.

[00126] A large number of radar sites are distributed across the world, e.g. US Weather Surveillance Radar - 1988 Doppler (WSR-88D), also known as Next Generation Radar (NEXRAD) network. These sites operate independently, but the ranges of these radar sites may overlap. The NEXRAD network may provide information associated with monitoring severe weather and issuing storm warnings, flash flood warnings. [00127] Gridded data may allow various WSR-88D users to benefit from a wide- variety of products and displays (flexible horizontal or vertical cross sections) that may be easily extracted from multiple radar analysis grids. Further, data fusion, i.e. the process of combining radar data with information from other sources e.g. satellite, may be performed for gridded radar data. The gridded data may provide a more complete depiction of storm and precipitation events than using a single radar. [O0128] A mosaic of data combined with radar images may be constructed from overlapping radars into a single image to potentially render a more accurate display. [00129] The following issues may be addressed using embodiments herein: Variability of the number of radar coverage volumes overlapping at a given location. These may range from 0 (no radar coverage) or 1 (coverage by a single radar) to coverage by several radars; variability of the size and orientation of the radar resolution volumes of multiple radars at a given location; the conical nature of the individual radar beam scanning surfaces for different elevation angles, and the non-overlap of these conical surfaces from adjacent radars; non-contiguous stacking of the scanning radar beams, i.e. presence of significant gaps between the scanning elevations of a radar, especially at the higher elevation angles; the curvature of the earth, which makes even the flat, zero- elevation scanning planes of adjacent radars to be non-overlapping, and uses additional formulation for transforming the coordinates from one radar to another, or to a common set of coordinates.

[00130] The following description includes terms, such as "up", "down", "upper", "lower", "first", "second", etc. that are used for descriptive purposes only and are not to be construed as limiting. The elements, materials, geometries, dimensions, and sequence of operations may all be varied to suit particular applications. Parts of some embodiments may be included in, or substituted for, those of other embodiments. While the foregoing examples of dimensions and ranges are considered typical, the various embodiments are not limited to such dimensions or ranges. [00131] The Abstract is provided to comply with 37 C.F.R. § 1.74(b) to allow the reader to quickly ascertain the nature and gist of the technical disclosure. The Abstract is submitted with the understanding that it will not be used to interpret or limit the scope or meaning of the claims. [00132] The illustrations of embodiments described herein are intended to provide a general understanding of the structure of various embodiments, and they are not intended to serve as a complete description of all the elements and features of apparatus and systems that might make use of the structures described herein. Many other embodiments will be apparent to those of skill in the art upon reviewing the above description. Other embodiments may be utilized and derived therefrom, such that structural and logical substitutions and changes may be made without departing from the scope of this disclosure. The figures are merely representational and may not be drawn to scale. Certain proportions thereof may be exaggerated, while others may be minimized. Accordingly, the specification and drawings are to be regarded in an illustrative rather than a restrictive sense.

[00133J In the foregoing Detailed Description, various features are grouped together in a single embodiment for the purpose of streamlining the disclosure. This method of disclosure is not to be interpreted as reflecting an intention that the claimed embodiments have more features than are expressly recited in each claim. Thus the following claims are hereby incorporated into the Detailed Description, with each claim standing on its own as a separate embodiment. One skilled in the art will recognize that other configurations are available and other methods of manufacture may function as well without exceeding the scope of the disclosed subject matter.

[00134] While particular embodiments have been illustrated and described, they are merely examples and a person skilled in the art may make variations and modifications to the embodiments described herein without departing from the spirit and scope of the presently disclosed subject matter.

Claims

Claims:
1. A system comprising: radar data received from a plurality of radars, the plurality of radars at separate geographical locations and having one or more overlapping radar regions; and a processor to reduce the data from the plurality of radars to a common universal latitude-longitude coordinate frame and to combine the data in the one or more overlapping radar regions.
2. The system of claim 1 further comprising a single image display of a group selected from combined radar data in the one or more overlapping regions and individual radar data in the one or more overlapping regions.
3. The system of claim 1 wherein the processor is to divide a latitude-longitude interval, spanning fields of view of the plurality of radars, into a plurality of blocks.
4. The system of claim 3 wherein a wind velocity is obtained, at one of the blocks, in three dimensions when there are at least three overlapping radars observing one of the plurality of blocks.
5. The system of claim 4 wherein when there are more than three overlapping radars observing one of the plurality of blocks, the processor uses least squares method to calculate the wind velocity, at the one of the blocks, in three dimensions.
6. The system of claim 1 wherein the plurality of radars includes a satellite.
7. The system of claim 1 wherein there are one or more coordinates associated with the common universal latitude-longitude coordinate frame, wherein at one of the one or more coordinates there are no overlapping radar regions.
8. The system of claim 1 wherein the radar data includes scalar data, and the scalar data obtained at one radar resolution cell is averaged when there are overlapping radar regions at the one radar resolution cell.
9. A system comprising: means for receiving radar data from a plurality of radars, the plurality of radars having one or more overlapping radar coverage regions; means for determining latitude-longitude coordinates associated with the radar data; and means for combining the data, in latitude-longitude coordinates, in the one or more overlapping regions.
10. The system of claim 9 wherein the means for combining the data in latitude- longitude coordinates includes a processor.
11. A method comprising: receiving radar data from a plurality of radars, the plurality of radars having one or more overlapping radar regions; determining latitude-longitude coordinates associated with the radar data; and combining the data in latitude-longitude coordinates in the one or more overlapping regions.
•12. The method of claim 11 further comprising displaying a single mosaiced image of a group selected from combined radar data in the one or more overlapping regions and individual radar data in the one or more overlapping regions.
13. The method of claim 11 further comprising dividing an interval of the latitude- longitude coordinates associated with the radar data into a plurality of blocks.
14. The method of claim 13 further comprising obtaining, at one of the blocks, a wind velocity in three dimensions when there are at least three overlapping radars observing one of the blocks.
15. The method of claim 14 further comprising using least squares method to calculate the wind velocity, at the one of the blocks, in three dimensions, when there are more than three overlapping radars observing one of the blocks.
16. The method of claim 13 further comprising obtaining, at one of the blocks, a wind velocity in two dimensions when there. are at least two overlapping radars observing one of the blocks.
17. The method of claim 11 wherein the radar data includes scalar data, the method further comprising averaging the scalar data obtained at one of the coordinates, when there are overlapping radar observations at the one of the coordinates.
18. A computer-readable medium having program instructions to perform the method of claim 11.
19. A method comprising: automatically determining a number of radars observing a given location in a multi-radar observation situation.
20. The method of claim 19 further comprising: automatically determining an algorithm to be used for combining multi-radar vector data: to yield radial velocity alone at a cell, when an observation at the cell is by a single radar; to yield a two-dimensional wind vector, when the observation at the cell is by two separate radars; to yield a three-dimensional wind vector, when the observation is by three separate radars; and to yield a least-squares three-dimensional wind vector, when the observation is by more than three radars.
21. The method of claim 19 further comprising: determining a position vector, in a latitude-longitude coordinate system, from a given radar of the number of radars to a center of a display block.
22. The method of claim 21 further comprising: calculating coordinates of a minimum bounding rectangle for radar i,
Figure imgf000029_0001
'Φmax, ))> ™d calculating coordinates of the minimum bounding rectangle coverage region for all the N radars, ((λaitt , φmin ), {λmax , φmax )) .
23. The method of claim 19 further comprising: displaying vector wind fields from data associated with the number of radars by using display blocks and arrows centered on the display blocks.
24. The method of claim 23 further comprising: combining vector data from the number of radars by prior-averaging individual radial components of the data before synthesizing wind vector over the display block.
25. The method of claim 22 further comprising: calculating a distance metric
Figure imgf000029_0002
and using designated conditions
Figure imgf000029_0003
may be retrieved from a generated first mechanism of CAPPI for the ith radar to test for coverage of a given coordinate by radar i.
26. The method of claim 22 further comprising: computing an east-west distance *fewi and an north-south distance ^NSΪ from the formulae d&m -
Figure imgf000030_0001
, respectively, of a latitude-longitude location (λ,φ) from the radar i; and retrieving PPI scan converted CAPPI scan line data from the radar i corresponding to coordinates closest to dΕwi and c?Nsi-
27. The method of claim 19 further comprising: linearly averaging reflectivity data available in dBZ units.
28. The method of claim 21 further comprising: j automatically detecting a variability of coverage pattern within a display block among the number of radars within individual display blocks.
29. The method of claim 28 further comprising: automatically handling a variability of coverage pattern within a display block by considering a largest region with a specific coverage pattern.
PCT/US2006/032175 2005-03-10 2006-08-17 Constant altitude plan position indicator display for multiple radars WO2007022376A2 (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
US11/206,696 US7821448B2 (en) 2005-03-10 2005-08-18 Constant altitude plan position indicator display for multiple radars
US11/206,696 2005-08-18

Publications (2)

Publication Number Publication Date
WO2007022376A2 true WO2007022376A2 (en) 2007-02-22
WO2007022376A3 WO2007022376A3 (en) 2007-09-07

Family

ID=37622030

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2006/032175 WO2007022376A2 (en) 2005-03-10 2006-08-17 Constant altitude plan position indicator display for multiple radars

Country Status (1)

Country Link
WO (1) WO2007022376A2 (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010130286A1 (en) * 2009-05-12 2010-11-18 Raytheon Anschütz Gmbh Combining data from multiple radar signals on a single plan position indicator (ppi) display
RU2574598C1 (en) * 2014-12-01 2016-02-10 ОТКРЫТОЕ АКЦИОНЕРНОЕ ОБЩЕСТВО "НИИ измерительных приборов - Новосибирский завод имени Коминтерна" /ОАО "НПО НИИИП-НЗиК"/ Radar space scanning method (versions)
US20160266249A1 (en) * 2015-03-13 2016-09-15 Honeywell International Inc. Weather radar integrating system combining ground-based and aircraft-based weather radar data

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3660844A (en) * 1970-05-04 1972-05-02 Sierra Research Corp Radar detector and identifier
DE3132009A1 (en) * 1981-08-13 1983-03-03 Licentia Gmbh Method for correcting systematic errors in multi-radar target tracking
US5554990A (en) * 1995-08-22 1996-09-10 Northrop Grumman Corporation Airspace management system and method
DE19614798A1 (en) * 1996-04-15 1997-10-16 Hannig Lutz Dipl Meteor Geographical, multistatic three=dimensional measuring platform
WO2001055742A2 (en) * 2000-01-12 2001-08-02 Honeywell International Inc. Projection of multi-sensor ray based data histories onto planar grids
EP1229348A1 (en) * 2001-01-30 2002-08-07 Meteo-France Method for detecting and characterising obstacles in a radar system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3660844A (en) * 1970-05-04 1972-05-02 Sierra Research Corp Radar detector and identifier
DE3132009A1 (en) * 1981-08-13 1983-03-03 Licentia Gmbh Method for correcting systematic errors in multi-radar target tracking
US5554990A (en) * 1995-08-22 1996-09-10 Northrop Grumman Corporation Airspace management system and method
DE19614798A1 (en) * 1996-04-15 1997-10-16 Hannig Lutz Dipl Meteor Geographical, multistatic three=dimensional measuring platform
WO2001055742A2 (en) * 2000-01-12 2001-08-02 Honeywell International Inc. Projection of multi-sensor ray based data histories onto planar grids
EP1229348A1 (en) * 2001-01-30 2002-08-07 Meteo-France Method for detecting and characterising obstacles in a radar system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
MAHAPATRA P R ET AL: "Vector reconstruction and mosaicing from multiple doppler weather radar velocity data" GEOSCIENCE AND REMOTE SENSING SYMPOSIUM, 2004. IGARSS '04. PROCEEDINGS. 2004 IEEE INTERNATIONAL ANCHORAGE, AK, USA 20-24 SEPT. 2004, PISCATAWAY, NJ, USA,IEEE, vol. 6, 20 September 2004 (2004-09-20), pages 3868-3871, XP010751598 ISBN: 0-7803-8742-2 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010130286A1 (en) * 2009-05-12 2010-11-18 Raytheon Anschütz Gmbh Combining data from multiple radar signals on a single plan position indicator (ppi) display
US9075145B2 (en) 2009-05-12 2015-07-07 Raytheon Anschutz Gmbh Combining data from multiple radar signals on a single plan position indicator (PPI) display
RU2574598C1 (en) * 2014-12-01 2016-02-10 ОТКРЫТОЕ АКЦИОНЕРНОЕ ОБЩЕСТВО "НИИ измерительных приборов - Новосибирский завод имени Коминтерна" /ОАО "НПО НИИИП-НЗиК"/ Radar space scanning method (versions)
US20160266249A1 (en) * 2015-03-13 2016-09-15 Honeywell International Inc. Weather radar integrating system combining ground-based and aircraft-based weather radar data
US10241203B2 (en) * 2015-03-13 2019-03-26 Honeywell International Inc. Weather radar integrating system combining ground-based and aircraft-based weather radar data

Also Published As

Publication number Publication date
WO2007022376A3 (en) 2007-09-07

Similar Documents

Publication Publication Date Title
US9874649B2 (en) Device and method for displaying full azimuth angle domain image data
CA2848216C (en) Device and method for 3d sampling with avian radar
Cheney et al. Fundamentals of radar imaging
Liu et al. A cloud and precipitation feature database from nine years of TRMM observations
Emery et al. Evaluating radial current measurements from CODAR high-frequency radars with moored current meters
Ford et al. Venus topography and kilometer‐scale slopes
Ye et al. Landslide monitoring in the Three Gorges area using D-InSAR and corner reflectors
Corbosiero et al. The effects of vertical wind shear on the distribution of convection in tropical cyclones
US20140063000A1 (en) Method and apparatus of taking aerial surveys
JP5726516B2 (en) Three-dimensional radar monitoring system and method
Korolev et al. Evaluation of the accuracy of PMS optical array probes
Trevett Imaging radar for resources surveys
Wright et al. Hurricane directional wave spectrum spatial variation in the open ocean
Richards A beginner's guide to interferometric SAR concepts and signal processing [aess tutorial iv]
Graham Synthetic interferometer radar for topographic mapping
US6930715B1 (en) Method, system and program product for augmenting an image of a scene with information about the scene
US6829536B2 (en) System and method for processing weather information
US20120200606A1 (en) Detail-in-context lenses for digital image cropping, measurement and online maps
Lakshmanan et al. The warning decision support system–integrated information
Bernstein Digital image processing of earth observation sensor data
Brock et al. Basis and methods of NASA airborne topographic mapper lidar surveys for coastal studies
Anagnostou et al. The use of TRMM precipitation radar observations in determining ground radar calibration biases
US6795019B2 (en) Method and system for calibrating wireless location systems
US5805098A (en) Method and system for forming image by backprojection
Barrick et al. Ocean surface currents mapped by radar

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application
NENP Non-entry into the national phase in:

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 06813510

Country of ref document: EP

Kind code of ref document: A2