Abstract:
The present invention relates to a method for storing measurements of a given type made by a radar, each measurement having been made at a position in a given coordinate system. The method comprises a step of accessing a structured data storage space, each location of which is addressable by a pair of positive integers and each location of which is able to store at least one evolution profile of the measurements of the given type as a function of altitude. It also comprises a meshing step which associates a ground position with a pair of positive integers allowing a location to be addressed. The method also comprises a step of modifying the content of the location with the pair of integers as its address by storing there an evolution profile such that for each ground position sufficiently close to the position associated with the pair and for which a measurement has been made at a certain altitude, the evolution profile provides approximately the measurement made at this altitude.

Description:
RELATED APPLICATIONS 
       [0001]    The present application is based on, and claims priority from, French Application Number 0704283, filed Jun. 15, 2007, the disclosure of which is hereby incorporated by reference herein in its entirety. 
       FIELD OF THE INVENTION 
       [0002]    The present invention relates to a method for storing measurements of a given type made by a radar. It applies, for example, in the field of weather radars. 
       BACKGROUND OF THE INVENTION 
       [0003]    A weather radar allows location of precipitations such as rain, snow or hail, measurement of their intensity and possibly identification of dangerous phenomena. Most weather radar are installed on the ground and are often part of a larger weather surveillance network. But airborne applications are increasingly coming into being, air transport being particularly interested in meteorological phenomena. The concern is notably with bypassing cumulonimbi, enormous clouds much feared by pilots as they sometimes create violent storms. Even airliners change their route in order to avoid crossing the path of certain particularly threatening cumulonimbi. This is because the lightning, the hail and the strong wind shears inside the cloud add to the risk of ice accretion and may endanger the flight if the pilot attempts to pass through. 
         [0004]    A weather radar allows the detection of extended voluminous targets that clouds are, of which it must provide the position, size and dangerousness. To do that, a weather radar may, for example, emit a wave in the X-band. The distance to a cloud is deduced from the time necessary for the emitted pulse to carry out a return trip from the radar antenna to the cloud at the speed of light. This time corresponds simply to the time between emitting a pulse and receiving its echo. Estimation of the size of a cloud entails estimating its extent, i.e. the maximum horizontal distance it extends over, and estimating its elevation, i.e. the maximum vertical distance it extends over. Estimation of the extent profits notably from the azimuth scan of the radar beam. Estimation of the elevation profits notably from the elevation scan of the radar beam. By way of an indication, the elevation of a cumulonimbus often exceeds 10000 meters! It is the elevation which primarily characterizes the dangerousness of the cloud, for the higher a convective cloud is, the more dangerous it is. But the danger level of the cloud is also linked with its reflectivity factor, designated Z, which characterizes the concentration of hydrometeors suspended in a volume of air, in liquid or solid form. In a way, the reflectivity factor Z represents the intensity of the cloud. Once in a logarithmic scale, it is represented in dBZ. 
         [0005]    In concrete terms, a display console displays to the on-board crew a simplified representation of clouds based on elementary geometrical shapes such as parallelepipeds, the color of which characterizes the intensity, whether it involves rain, snow or hail. For example, the color black is often used for dry air, i.e. the absence of cloud. Green and yellow may be used for moderate humidity concentrations. Red is often used for areas with very high humidity concentration, i.e. the most dangerous areas which it is absolutely necessary to bypass. 
         [0006]    The display must be almost instantaneous, possibly on several screens used by various crew members throughout the flight. This simplified graphical representation of clouds is constructed on the screen using previously made intensity measurements, these measurements being stored for a greater or lesser period in a suitable memory space. Access to the data must therefore be fast and efficient, which is not without difficulties when the large quantity of data to be stored and therefore the memory space to be addressed is considered. In fact, it involves still making an intensity value available for every position (x, y, z) of space within range of the radar! The data must therefore be stored in a structured manner in order to optimize access to them while avoiding notably accessing them sequentially. This is one of the technical problems to which the present invention proposes to provide an innovative solution. 
         [0007]    One current solution consists in saving the intensity information in a 3D matrix with each dimension corresponding to a spatial dimension. Hence, a triplet of indices (i, j, k) corresponds to each position (x, y, z), the coefficient of the 3D matrix at the location (i, j, k) containing an intensity value associated with the position (x, y, z). But such a solution is costly in terms of memory space: whatever the quantity of relevant information saved, the purely spatial quantification in a matrix makes it necessary to reserve a large memory space. For the matrix may possibly be largely filled with zero or insignificant values over an entire area. Thus, the memory is needlessly burdened, notably by clear weather, which is particularly damaging in airborne systems, the resources of which are limited. Of course, the problem is accentuated when the display precision increases, i.e. when the surface covered by a coefficient (i, j, k) of the 3D matrix decreases while the weather area to be stored remains the same. 
       SUMMARY OF THE INVENTION 
       [0008]    The aim of the invention is notably to propose an advantageous compromise between the memory space used and the spatial precision permitted, and to avoid storing too many zero or insignificant values. For a given 2D geographical ground area, it involves notably not storing all the intensity values measured for all the altitude levels, but instead storing a parametric representation of these values. To this end, it is functions that are stored, each function representing the evolution of the intensity as a function of the altitude at a considered point on the ground. 
         [0009]    To this end, the subject of the invention is a method for storing measurements of a given type made by a radar, each measurement having been made at a position in a given coordinate system. The method comprises a step of accessing a structured data storage space, each location of which is addressable by a pair of positive integers and each location of which is able to store at least one evolution profile of the measurements of the given type as a function of the altitude. It also comprises a meshing step which associates a ground position with a pair of positive integers allowing a location to be addressed. It also comprises a step of modifying the content of the location with the pair of integers as its address by storing there an evolution profile such that for each ground position sufficiently close to the position associated with the pair and for which a measurement has been made at a certain altitude, the evolution profile provides approximately the measurement made at this altitude. 
         [0010]    In one embodiment, several evolution profiles may be stored at the location with the pair as its address. It is possible that a single one of said profiles may be nonzero for a given altitude. 
         [0011]    Advantageously, the profile may be stored in the form of a structure containing a p-tuple of values, p being a positive nonzero integer. The evolution profile may thus be a piecewise constant function assuming p distinct nonzero values, or a polynomial function of degree p−1 or again a spatial distribution defined by p−2 real coefficients, a mean and a standard deviation. It may also be linear by segment. 
         [0012]    For example, the measurements may be values characterizing the reflectivity of clouds or their danger level made by an airborne weather radar. 
         [0013]    A principal advantage of the invention is also that it enables space to no longer be wasted for storing useless values, such as the values corresponding to areas without cloud or without a weather danger. 
         [0014]    Still other objects and advantages of the present invention will become readily apparent to those skilled in the art from the following detailed description, wherein the preferred embodiments of the invention are shown and described, simply by way of illustration of the best mode contemplated of carrying out the invention. As will be realized, the invention is capable of other and different embodiments, and its several details are capable of modifications in various obvious aspects, all without departing from the invention. Accordingly, the drawings and description thereof are to be regarded as illustrative in nature, and not as restrictive. 
     
    
     
       DETAILED DESCRIPTION OF THE DRAWINGS 
         [0015]    The present invention is illustrated by way of example, and not by limitation, in the figures of the accompanying drawings, wherein elements having the same reference numeral designations represent like elements throughout and wherein: 
           [0016]      FIG. 1 , an illustration through a view of a screen on board an aircraft of an example where the parametric functions are simple piecewise constant functions; and 
           [0017]      FIG. 2 , an illustration through a diagram of an example of enriched structures to store the parametric functions. 
       
    
    
     DETAILED DESCRIPTION OF THE INVENTION 
       [0018]    The method according to the present invention proposes storing the measurements made by a radar, such as the reflectivity of clouds or their danger level in the case of a weather radar, in a way enabling the quantity of information saved to be significantly limited. Advantageously, the method may be simply implemented in the majority of computer systems for processing radar signals, whether or not they are airborne, and this without modification of their hardware or software configuration. To do this, a section of random access memory (RAM) is structured in a two-dimensional array, each cell of the array being accessible through a pair of positive integer indices (i,j). The array forms a 2D ground mesh, that is to say that a ground position (x i ,y j ,0) in a given coordinate system, x i  and y j  being real numbers, corresponds to each cell with indices (i,j). In addition, each cell with indices (i,j) contains a complex data item enabling deduction of the reflectivity of clouds or their danger level at the vertical of the position (x i ,y j ,0) as a function of the altitude z. In other words, the complex data item stored in the cell (i,j) enables deduction of the reflectivity of clouds, their danger level, or any other quantity in any position (x i ,y j ,z). For example, the complex data item may be a record E ij ={S 1 , S 2 , . . . ,S N } of N structures S k , N being a positive nonzero integer and k an integer between 1 and N, each structure S k  enabling description of a parametric function f which associates a cloud reflectivity, a danger level or any other quantity with a height h. As illustrated by the three following examples, the type of parametric function f may vary without departing from the principles of the present invention: f may be a piecewise constant function, a polynomial function or again a spatial distribution. The invention thus requires storing in each cell of the array with indices (i,j) at least one structure representing a parametric function f such that, for any position (x,y,0) sufficiently close to (x i ,y j ,0) below a certain threshold and for which a measurement M(x,y,z) has been made at an altitude z, the measurement M(x,y,z) is sufficiently close to f(z) below a certain threshold in the sense of a norm suited to the measurement type. Advantageously, several structures representing several parametric functions may be stored in a cell with indices (i,j). In this case, advantageously but not exclusively, a single one of the stored parametric functions is nonzero in z for a given altitude z. In another embodiment, f may be a linear segment function. 
         [0019]      FIG. 1  illustrates through a view of a screen on board an aircraft a first example where the parametric functions are simple piecewise constant functions.  FIG. 1  shows how the invention may, for example, allow the reflectivity of clouds or their danger level to be displayed simply by a color code, the measurements having been made by an airborne radar. The screen displays a typical cloud  1  represented in vertical section in a graph, the x-axis indicating the distance relative to the radar and the y-axis indicating the altitude in feet relative to the ground. On the right of the graph, the screen also displays a color palette allowing each reflectivity level to be linked with a color. For example a reflectivity level less than 20 dBZ is displayed in a color c 0 , for example black. A level between 20 dBZ and 30 dBZ is displayed in a color c 1 , for example green. A level between 30 dBZ and 40 dBZ is displayed in a color c 2 , for example yellow. Finally, a level higher than 40 dBZ is displayed in a color c 3 , for example red. Typically, in terms of mean reflectivity value and for gaps between colors of around 10 dB, a cloud is characterized by a succession of intensity measurements that grow going toward the center of the cloud, then decrease when moving away from the center of the cloud. In the illustration of the cloud  1  this typical phenomenon is brought out by the column structures  2 ,  3  and  4  corresponding respectively to distances of around 1.15, 1.25 and 1.3 relative to the radar. In column  2 , there is hence a transition from black to green at an altitude h 1  of around 1000 feet, from green to yellow at an altitude h 2  of around 2500 feet, from yellow to red at an altitude h 3  of around 5000 feet, from red to yellow at an altitude h 4  of around 8000 feet, from yellow to green at an altitude h 5  of around 9000 feet, and finally from green to black at an altitude h 6  of around 12000 feet. A structure S enabling representation of column  2  may be defined by a record of 6 heights S={h 1 ,h 2 , . . . h 6 } in which each height respectively corresponds to a height of a color change: h 1  for the transition from black to green, h 2  for the transition from green to yellow, h 3  for the transition from yellow to red, h 4  for the transition from red to yellow, h 5  for the transition from yellow to green, h 6  for the transition from green to black. It seems that by postulating that the measurement of the intensity of a cloud is graduated, it is possible to limit the quantity of saved information. Thus a very small number of parameters, such as the six parameters of the present example, can enable the display of a reflectivity level without error in quantifying the information. Notably, the use of information about the graduation of the measurement of intensity of a cloud limits a priori the quantity of saved information. 
         [0020]    More generally, each parametric function f may be defined by a piecewise constant function in the form of a p-tuple of heights {h 1 ,h 2 , . . . h p }, where p is a nonzero integer, such that: 
         [0000]    
       
         
           
             
               f 
                
               
                 ( 
                 h 
                 ) 
               
             
              
             
               : 
             
              
             
               { 
               
                 
                   
                     
                       
                         
                           h 
                           1 
                         
                         &gt; 
                         
                           h 
                            
                           
                             f 
                              
                             
                               ( 
                               h 
                               ) 
                             
                           
                         
                       
                       = 
                       
                         z 
                         0 
                       
                     
                   
                 
                 
                   
                     
                       
                         
                           h 
                           1 
                         
                         &lt; 
                         h 
                         &lt; 
                         
                           
                             h 
                             2 
                           
                            
                           
                             f 
                              
                             
                               ( 
                               h 
                               ) 
                             
                           
                         
                       
                       = 
                       
                         z 
                         1 
                       
                     
                   
                 
                 
                   
                     … 
                   
                 
                 
                   
                     
                       
                         
                           h 
                           
                             p 
                             - 
                             1 
                           
                         
                         &lt; 
                         h 
                         &lt; 
                         
                           
                             h 
                             p 
                           
                            
                           
                             f 
                              
                             
                               ( 
                               h 
                               ) 
                             
                           
                         
                       
                       = 
                       
                         z 
                         p 
                       
                     
                   
                 
                 
                   
                     
                       h 
                       &gt; 
                       
                         
                           
                             h 
                             p 
                           
                            
                           
                             f 
                              
                             
                               ( 
                               h 
                               ) 
                             
                           
                         
                          
                         
                           z 
                           0 
                         
                       
                     
                   
                 
               
             
           
         
       
     
         [0000]    where the values {z 0 ,z 1 , . . . z p } are mean reflectivity values corresponding respectively to the colors {c 0 ,c 1 , . . . c p } fixed for all the structures S k . In this case, a structure S k  is defined by a record {h 1 ,h 2 , . . . h p }. Thus from h 1  to h 2  a color c 1  is displayed, from h 2  to h 3  a color c 2  is displayed, and so on through to h p−1  to h p  where the color c p−1  is displayed. For the heights above h p  or below h 1 , the color c 0  is displayed. In another embodiment, a continuous linear segment parametric function may be defined in a similar manner by a p-tuple of heights {h 1 ,h 2 , . . . h p } defining the points at which the slope of the function changes. 
         [0021]      FIG. 2  illustrates through a diagram an example of structures enriched in relation to the structure S of the example in  FIG. 1 . The structures are represented over a two-dimensional matrix Π which stores them. The matrix Π comprises 4 rows indexed from i=1 to i=4 and 11 columns indexed from j=1 to j=11. In the interests of clarity, only the index of the row i=4 and the index of the column j=1 have been represented. For example, in a structure S 1  stored in the cell with the indices i=3 and j=3 and defined by a record {h 10 , h 11 , h 12 , h 13 , h 14 , h 15 , h 16 }, the cloud  1  touches the ground G at a height h 10 , the height h 10  providing the altitude of the ground and h 10  being equal to a height h 11  providing the base of the cloud  1 . In the structure S 1 , the intensity measurement reaches a maximum value Z max1  at an altitude h Zmax1 . As another example, in a structure S 2  stored in the cell with the indices i=2 and j=8, defined by the record {h 20 , h 21 , h 22 , h 23 , h 24 , h 25 , h 26 }, the color yellow does not appear between the green and the red, which is marked by the fact that the height h 22  is equal to the height h 23 . Similarly, green does not appear above the yellow, which is marked by the fact that the height h 26  is equal to the height h 25 . In the structure S 2 , the intensity measurement reaches a maximum value Z max2  at an altitude h Zmax2 . As another example, a structure S 3  stored in the cell with the indices i=3 and j=10, defined by the record {h 30 , h 31 , h 32 , h 33 , h 34 , h 35 , h 36 }, and a structure S 4  also stored in the cell with the indices i=3 and j=10, defined by the record {h 41 , h 42 , h 43 , h 44 , h 45 , h 46 }, both correspond to the same position on the ground. This is marked by the absence of a height h 40  in S 4  to represent the altitude of the ground. In the structure S 3 , the intensity measurement reaches a maximum value Z max3  at an altitude h Zmax3 . In the structure S 4 , the intensity measurement reaches a maximum value Z max4  at an altitude h Zmax4 . 
         [0022]    The function f may also be a polynomial function defined by a p-tuple (a 0 , a 1 , . . . ,a p−1 ) of coefficients: 
         [0000]        f ( h )= a   0   +a   1   h+a   2   h   2   + . . . +a   p−1   h   p−1    
         [0023]    In this case, a structure S k  is defined by the record (a 0 , a 1 , . . . ,a p−1 ). 
         [0024]    The parametric function f may also be a spatial distribution around a pixel (x 0 ,y 0 ): 
         [0000]    
       
         
           
             
               
                 f 
                 
                   
                     z 
                     0 
                   
                    
                   
                     y 
                     0 
                   
                 
               
                
               
                 ( 
                 h 
                 ) 
               
             
             = 
             
               
                 ∑ 
                 
                   i 
                   = 
                   0 
                 
                 
                   p 
                   - 
                   1 
                 
               
                
               
                 
                   k 
                   i 
                 
                  
                 
                   exp 
                    
                   
                     [ 
                     
                       - 
                       
                         
                           
                             ( 
                             
                               h 
                               - 
                               
                                 m 
                                 h 
                               
                             
                             ) 
                           
                           2 
                         
                         
                           σ 
                           h 
                           2 
                         
                       
                     
                     ] 
                   
                 
               
             
           
         
       
     
         [0025]    In this case, a structure S k  is defined by a record {k 0 , k 1 , . . . ,k p−1 , m h ,σ h } where k 0 , k 1 , . . . ,k p−1  are real coefficients, m h  is the mean of the heights at which measurements have been made on the vertical of (x 0 ,y 0 ) and σ h  is the standard deviation of the heights at which measurements have been made on the vertical of (x 0 ,y 0 ). 
         [0026]    It must be noted that the previously described invention is not applicable only to the measurement of the intensity or reflectivity of a cloud. Without departing from the principles of the invention it is also applicable to any measured or calculated quantity likely to vary as a function of altitude, such as, for example, the risk level, the standard deviation in speed or again the integral scale. Moreover, the previously described invention enables the rapid increase in calculation power in relation to the available memory space to be profited from fully, especially in airborne systems. This is because the parametric representation of the reflectivity level according to the invention allows the relevant information to be compressed and hence to save memory space. Notably, the use of information about the graduation of the measurement of intensity of a cloud limits a priori the quantity of saved information. Hence, the transfer rate necessary to transfer information characterizing a weather situation is considerably reduced. The disappearance of quantifying information in terms of altitude simplifies the display in a panoramic view at a constant angle of elevation, the type of display better known by the acronym PPI meaning “Plan Position Indicator”. The parametric nature of the representation according to the invention also allows calculations of the vertical tendency to be made more easily and more precisely. 
         [0027]    It will be readily seen by one of ordinary skill in the art that the present invention fulfills all of the objects set forth above. After reading the foregoing specification, one of ordinary skill in the art will be able to affect various changes, substitutions of equivalents and various aspects of the invention as broadly disclosed herein. It is therefore intended that the protection granted hereon be limited only by definition contained in the appended claims and equivalents thereof.