Terrain awareness system

0001The present application claims priority of Provisional <patcit id="pcit0001" dnum="US02509196P" dnum-type="L"><text>U.S. Patent Application Serial No. 60/025,091, filed September 3, 1996</text></patcit>, which is a continuation-in-part of commonly assigned, co-pending patent application Serial No. <patcit id="pcit0002" dnum="US08509642B"><text>08/509,642</text></patcit>, entitled "Terrain Awareness System", filed July 31, 1995.

BACKGROUND OF THE INVENTION

1.

Field of the Invention

0002The present invention relates to a terrain awareness system and more particularly to a system for alerting a pilot of an aircraft of a dangerous flight condition based upon stored terrain data which provides relatively longer warning times than known ground proximity warning systems while minimizing nuisance warnings.

2.

Description of the Prior Art

0003Various systems are known in the art that provide warnings and advisory indications of hazardous flight conditions. Among such systems are systems generally known as ground proximity warning systems (GPWS) which monitor the flight conditions of an aircraft and provide a warning if flight conditions are such that inadvertent contact with terrain is imminent. Among the flight conditions normally monitored by such systems are radio altitude and rate, barometric altitude and rate, air speed, flap and gear positions. These parameters are monitored and an advisory signal and/or warning signal is generated when the relationship between the parameters is such that terrain impact is likely to occur. Typical examples of such systems are disclosed in <patcit id="pcit0003" dnum="US3715718A"><text>U.S. Patent Nos. 3,715,718</text></patcit>; <patcit id="pcit0004" dnum="US3936796A"><text>3,936,796</text></patcit>; <patcit id="pcit0005" dnum="US3958218A"><text>3,958,218</text></patcit>; <patcit id="pcit0006" dnum="US3944968A"><text>3,944,968</text></patcit>; <patcit id="pcit0007" dnum="US3947808A"><text>3,947,808</text></patcit>; <patcit id="pcit0008" dnum="US3947810A"><text>3,947,810</text></patcit>; <patcit id="pcit0009" dnum="US3934221A"><text>3,934,221</text></patcit>; <patcit id="pcit0010" dnum="US3958219A"><text>3,958,219</text></patcit>; <patcit id="pcit0011" dnum="US3925751A"><text>3,925,751</text></patcit>; <patcit id="pcit0012" dnum="US3934222A"><text>3,934,222</text></patcit>; <patcit id="pcit0013" dnum="US4060793A"><text>4,060,793</text></patcit>; <patcit id="pcit0014" dnum="US4030065A"><text>4,030,065</text></patcit>; <patcit id="pcit0015" dnum="US4215334A"><text>4,215,334</text></patcit>; and <patcit id="pcit0016" dnum="US4319218A"><text>4,319,218</text></patcit>, all assigned to the same assignee as the assignee of the present invention.

0004While the above-referenced systems do provide advisory and warning signals in the event of proximity to terrain, the warnings generated by such systems are based solely upon flight conditions of the aircraft and do not utilize any navigational information. Consequently, the sensitivity of such systems must be adjusted to provide adequate warnings when a hazardous flight condition exists without generating false or spurious warnings. However, such an adjustment can result in a compromise that may still result in nuisance warnings over terrain unique to particular geographic areas and shorter than desired warning times in yet other geographic areas.

0005Several attempts have been made to improve upon such ground proximity warning systems utilizing ground-based navigational information. For example, <patcit id="pcit0017" dnum="US4567483A"><text>U.S. Patent Nos. 4,567,483</text></patcit>; <patcit id="pcit0018" dnum="US4646244A"><text>4,646,244</text></patcit>; <patcit id="pcit0019" dnum="US4675823A"><text>4,675,823</text></patcit>; and <patcit id="pcit0020" dnum="US4914436A"><text>4,914,436</text></patcit> all disclose ground proximity warning systems which monitor the position of the aircraft relative to stored terrain data in order to provide modified ground proximity warnings. However, the utility of such systems is limited. For example, the systems disclosed in <patcit id="pcit0021" dnum="US4567483A"><text>U.S. Patent Nos. 4,567,483</text></patcit> and <patcit id="pcit0022" dnum="US4914436A"><text>4,914,436</text></patcit> disclose ground proximity warning systems which utilize navigational data to modify predetermined warning envelopes surrounding certain particular airports.

0006<patcit id="pcit0023" dnum="US4646244A"><text>U.S. Patent Nos. 4,646,244</text></patcit> and <patcit id="pcit0024" dnum="US4675823A"><text>4,675,823</text></patcit> disclose terrain advisory systems which utilize various ground-based navigational inputs and stored terrain data to provide various ground proximity warning systems based on the position of the aircraft. In order to conserve memory, the terrain data is modeled in various geometric shapes as a function of the elevation of the terrain within the area defined by the geometric shape. The warning signals are generated as a function of the position of the aircraft relative to the model. While such systems do provide adequate terrain warnings, such systems may cause nuisance warnings under certain conditions since the predicted trajectory of the aircraft is not taken into consideration.

0007Another problem with such systems is that they are based upon ground-based navigational systems, such as VOR/DME and LORAN. Such ground-based navigational systems are being replaced with relatively more accurate satellite systems, such as the global positioning system (GPS). Moreover, such ground-based systems are only able to provide two-dimensional position data. Thus, with such systems, the altitude of the aircraft must be supplied by an auxiliary device, such as an altimeter, which increases the complexity of the system as well as the cost.

0008Enhanced ground proximity warnings (EGPWS), windshear warnings including predictive (PrWS) and reactive (ReWS) and TCAS alerts are provided in <patcit id="pcit0025" dnum="US5220322A"><text>U.S. Patent Nos. 5,220,322</text></patcit>, <patcit id="pcit0026" dnum="US4433323A"><text>4,433,323</text></patcit>, <patcit id="pcit0027" dnum="US5187478A"><text>5,187,478</text></patcit> and <patcit id="pcit0028" dnum="US4684948A"><text>4,684,948</text></patcit>. Suitable PrWS and ReWS systems are described in <patcit id="pcit0029" dnum="US4905000A"><text>U.S. Patent Nos. 4,905,000</text></patcit>, <patcit id="pcit0030" dnum="US5059964A"><text>5,059,964</text></patcit>, <patcit id="pcit0031" dnum="US4725811A"><text>4,725,811</text></patcit>, <patcit id="pcit0032" dnum="US4947164A"><text>4,947,164</text></patcit>, <patcit id="pcit0033" dnum="US5153588A"><text>5,153,588</text></patcit> and <patcit id="pcit0034" dnum="US4891642A"><text>4,891,642</text></patcit>. Suitable commercially available TCAS systems and associated hardware are described in <patcit id="pcit0035" dnum="US5122808A"><text>U.S. Patent Nos. 5,122,808</text></patcit>, <patcit id="pcit0036" dnum="US5272725A"><text>5,272,725</text></patcit>, <patcit id="pcit0037" dnum="US4914733A"><text>4,914,733</text></patcit>, <patcit id="pcit0038" dnum="US5008844A"><text>5,008,844</text></patcit> and <patcit id="pcit0039" dnum="US4855748A"><text>4,855,748</text></patcit>, and in published documents: "Pilot's Guide TCAS II CAS 67A/81A Bendix/King Traffic Alert and Collision Avoidance Systems" and "Pilot's Guide CAS 66A Bendix/King TCAS I Collision Avoidance System". Performance specifications on TCAS may be found in TCAS I - Advisory Circular (AC) 20-TCAS draft dated 08115/95; TCAS II - Advisory Circular 20-131 dated -3/29/93; DO 181 A RTCA MOPS Mode S Transponder; DO 185 TCAS RTCA MOPS and DO 185A TCAS RTCA MOPS.

SUMMARY OF THE INVENTION

0009The present invention provides a warning system for aircraft comprising: a terrain awareness device for receiving and storing terrain information relative to a position including altitude of an aircraft; a visual display, coupled to the terrain awareness device, for displaying the terrain information to the pilot; a controller, coupled to the visual display, adapted to automatically determine when the terrain information should be displayed on the visual display and alert the pilot accordingly, the controller comprising a selection switch constructed to automatically switch from displaying other information to displaying terrain information on the visual display when a terrain threat is detected, characterised by a terrain selection lamp indicating that the terrain information has been selected.

0010This invention and a method of using the invention system can provide increased warning times to the pilot of an aircraft or hazardous flight condition while minimizing nuisance warnings. An important aspect of the invention is that the terrain background information as well as the terrain advisory and warning indications may be displayed on a navigational or weather display, normally existing within an aircraft, which reduces the cost of the system and obviates the need for extensive modification to existing displays, such as navigational and weather radar type display.

0011<patcit id="pcit0040" dnum="FR2689668A"><text>FR-A-2 689 668</text></patcit> discloses an aerial navigation aid which provides for an alarm and display of a terrain image as stored on an on-board memory.

0012In a specific aspect of the invention, a terrain warning system includes a terrain awareness device for receiving and storing terrain information relative to a position of an aircraft, a visual display for displaying the terrain information to the pilot and a controller for determining when the terrain information is displayed. The controller includes a terrain selection switch for switching between terrain information and other information, such as weather and navigational information. Preferably, the terrain information on the visual display changes color based on the severity of the terrain threat.

0013In one embodiment, the controller can be moved between a first state, in which the visual display depicts no information, a second state in which the visual display depicts terrain information and a third state in which the visual display depicts weather information. The terrain selection switch moves the controller between the second and third states. A master selection switch (i.e., a WXR switch) moves the controller between the first and second states. For example, when the controller is in the first or no information displayed state, toggling the WXR switch to ON moves the controller to the second state where weather is displayed. Toggling the terrain selection switch from the second state moves the controller to the third state where the terrain information is displayed. The controller moves from the second state back to the first state when the master selection switch is turned OFF.

0014The present invention may also include a lamp assembly for providing immediate information regarding the terrain threat and the state of the terrain awareness system to the pilot. In one embodiment, the terrain warning system includes a terrain selection lamp for indicating that the terrain information has been selected and a terrain display lamp for indicating that the terrain information is being displayed on the visual display. The warning system also includes a terrain alert lamp for indicating the terrain threat. In this embodiment, the pilot will notice the terrain alert lamp and toggle the terrain selection switch to view the terrain threat. Alternatively, the system may automatically display the terrain during a terrain alert.

BRIEF DESCRIPTION OF THE DRAWINGS

0015<ul id="ul0001" list-style="none" compact="compact"><li><figref idref="f0001">FIGS. 1A</figref> and <figref idref="f0002">1B</figref> represent a block diagram of a terrain awareness system (TAS) in accordance with the present invention;</li><li><figref idref="f0003">FIG. 2</figref> is a diagram of the world broken down into various latitude segments and longitude segments in accordance with the present invention;</li><li><figref idref="f0004">FIG. 3</figref> is a graphical representation of the various memory map and resolution sequencing in accordance with the present invention;</li><li><figref idref="f0003">FIG. 4A</figref> is an exemplary representation of the digital header and subsquare mask word for identifying geographical areas in accordance with the present invention;</li><li><figref idref="f0005">FIG. 4B</figref> is a terrain map alternative to <figref idref="f0004">FIG. 3</figref>;</li><li><figref idref="f0006">FIG. 5</figref> is a graphical illustration of the LOOK-AHEAD distance assuming the aircraft turns at a 30° angle in accordance with the present invention;</li><li><figref idref="f0006">FIG. 6</figref> is a graphical representation of a DH terrain floor boundary which forms the bases for a terrain advisory signal and a terrain warning signal in accordance with the present invention;</li><li><figref idref="f0007">FIG. 7</figref> is a graphical illustration of a terrain advisory signal for an aircraft relative to the terrain and an airport for a condition when the aircraft flight path angle is less than a first predetermined reference plane or datum in accordance with the present invention;</li><li><figref idref="f0007">FIG. 8</figref> is similar to <figref idref="f0007">FIG. 7</figref> but for a condition when the aircraft flight path angle is greater than the first predetermined reference plane;</li><li><figref idref="f0008">FIG. 9</figref> is a graphical illustration of a terrain warning signal for an aircraft relative to the terrain and an airport for a condition when the aircraft flight path angle is greater than a second reference plane or datum in accordance with the present invention;</li><li><figref idref="f0008">FIG. 10</figref> is similar to <figref idref="f0008">FIG. 9</figref> but for when the flight path angle of the aircraft is less than the second reference plane or datum;</li><li><figref idref="f0009">FIG. 11</figref> is a graphical illustration of a cut-off angle correction boundary for a level flight condition in accordance with the present invention;</li><li><figref idref="f0009">FIG. 12</figref> is similar to <figref idref="f0009">FIG. 11</figref> but for a condition when the flight path angle of the aircraft is greater than a predetermined reference plane or datum;</li><li><figref idref="f0009">FIG. 13</figref> is similar to <figref idref="f0009">FIG. 11</figref> but for a condition when the flight path angle of the aircraft is less than a predetermined reference plane or datum and also illustrates a BETA sink rate enhancement boundary in accordance with the present invention;</li><li><figref idref="f0010">FIG. 14</figref> is a graphical illustration of the terrain advisory and superimposed terrain warning signals for an aircraft relative to the terrain for a condition when the aircraft flight path angle is less than a first reference plane or datum;</li><li><figref idref="f0010">FIG. 15</figref> is similar to <figref idref="f0010">FIG. 14</figref> but for a condition when the flight path angle is between a first datum and a second datum;</li><li><figref idref="f0010">FIG. 16</figref> is similar to <figref idref="f0010">FIG. 14</figref> but for a condition when the flight path angle is less than the second datum;</li><li><figref idref="f0011">FIG. 17</figref> is a graphical illustration of LOOK-AHEAD/LOOK-DOWN terrain advisory and warning boundaries in accordance with the present invention for a condition when the aircraft is descending;</li><li><figref idref="f0011">FIG. 18</figref> is similar to <figref idref="f0011">FIG. 17</figref> but for a condition when an aircraft is climbing;</li><li><figref idref="f0012">FIG. 19</figref> is a graphical illustration of LOOK-UP terrain advisory and warning boundaries in accordance with the present invention;</li><li><figref idref="f0013">FIG. 20</figref> is a graphical illustration of an aircraft during a pull-up maneuver;</li><li><figref idref="f0014">FIG. 21</figref> is an exemplary block diagram for a system for generating a signal representative of the altitude loss due to pilot reaction time ALPT, as well as altitude loss due to a pull-up maneuver ALPU in accordance with the present invention;</li><li><figref idref="f0015">FIG. 22</figref> is a graphical illustration of an alternative methodology for generating a cut-off altitude boundary in accordance with the present invention;</li><li><figref idref="f0016">FIG. 23</figref> is a block diagram of the display system in accordance with the present invention;</li><li><figref idref="f0017">FIG. 24</figref> is a functional diagram of the display data format in accordance with the ARINC 708/453 standard for a weather radar display;</li><li><figref idref="f0018">FIG. 25</figref> is a plan view of a background terrain display illustrating the variable density dot pattern in accordance with the present invention;</li><li><figref idref="f0019">FIG. 26</figref> is similar to <figref idref="f0018">FIG. 25</figref> but additionally illustrates the terrain threat indications;</li><li><figref idref="f0020">FIG. 27</figref> is a block diagram illustrating the configuration of a flash ROM for a terrain data base and a RAM used with the present invention;</li><li><figref idref="f0021">FIG. 28</figref> illustrates a wedding-cake configuration for a RAM utilized in the present invention;</li><li><figref idref="f0022">FIG. 29</figref> is a diagram of the file format of the files in the terrain data base in accordance with the present invention;</li><li><figref idref="f0022">FIG. 30</figref> illustrates a terrain map of the terrain data base in accordance with the present invention illustrating the highest elevation in each of the cells within the map;</li><li><figref idref="f0023">FIG. 31</figref> is similar to <figref idref="f0022">FIG. 30</figref> and illustrates a method for correcting various cells in the vicinity of a runway;</li><li><figref idref="f0023">FIG. 32</figref> is a simplified elevational view of certain data illustrated in <figref idref="f0023">FIG. 31</figref>;</li><li><figref idref="f0023">FIG. 33</figref> is a diagram of a display illustrating the weather radar sweep spokes and the determination of the range increments;</li><li><figref idref="f0024">FIG. 34</figref> is a block diagram of the display system in accordance with the present invention illustrating the memory configuration and the data flow;</li><li><figref idref="f0025">FIG. 34A</figref> is a block diagram of a configuration alternative to that shown in <figref idref="f0024">FIG. 34</figref>;</li><li><figref idref="f0026">FIG. 35</figref> is a diagram of a terrain map illustrating the determination of the look-ahead vector array in accordance with the present invention;</li><li><figref idref="f0023">FIG. 36</figref> is a diagram of a display illustrating a radar sweep spoke and a method for determining the range increments of the radar sweep spoke;</li><li><figref idref="f0027">FIG. 37</figref> is a terrain map illustrating the computation of the hazardous terrain in accordance with the present invention;</li><li><figref idref="f0028">FIG. 38</figref> is a diagram of a terrain map superimposed on a display screen illustrating the sweep range and the determination of the incremental range increments;</li><li><figref idref="f0029">FIG. 39</figref> is a diagram of a pixel map of the display in accordance with the present invention;</li><li><figref idref="f0030">FIG. 40</figref> is a diagram of the different fractal patterns used to provide the variable density display of non-hazardous display indications in accordance with the present invention;</li><li><figref idref="f0031">FIG. 41</figref> is a diagram illustrating the determination of the X and Y increments of the display screen;</li><li><figref idref="f0031">FIG. 42</figref> is a diagram of the method for superimposing the fractal patterns on the displayed pixel map in accordance with the present invention;</li><li><figref idref="f0032">Fig 43</figref> is an illustration of a terrain data format in accordance with the present invention;</li><li><figref idref="f0033">Fig. 44</figref> is illustrative of steps to build the Terrain database from the original data in accordance with the present invention;</li><li><figref idref="f0034">Fig. 45</figref> is illustrative of steps to build the Terrain database from the original data in accordance with the present invention;</li><li><figref idref="f0035">Fig. 46</figref> is a flow diagram of terrain-weather selection states of one embodiment of the present invention;</li><li><figref idref="f0036">Fig. 47</figref> is a flow diagram illustrating a method for displaying terrain information according to one embodiment of the present invention; and</li><li><figref idref="f0037">Fig. 48</figref> is a block diagram of the terrain selection/deselection and annunciation according to the present invention.</li></ul>

DESCRIPTION OF THE PREFERRED EMBODIMENT

0016Referring to <figref idref="f0001 f0002">FIG. 1</figref>, a terrain awareness system (TAS), generally identified with the reference number 20, is illustrated. As will be discussed in more detail below, the TAS 20 utilizes inputs from a satellite-based navigational system, such as global positioning system (GPS) 22 (longitude, latitude, altitude, groundtrack, ground speed), and/or an FMS/IRS navigational system which may be updated by the GPS and/or DME/DME, a terrain data base 24, an airport data base 26 and corrected barometric altitude, for example from a barometric altimeter 28 to provide a LOOK AHEAD warning system which provides relatively longer warning times while minimizing nuisance warnings. Since the TAS 20 does not require a radio altimeter input, the system can be used with certain aircraft, such as commuter aircraft, which are not normally equipped with a radio altimeter.

0017The current longitude and latitude of the aircraft from the GPS 22 are applied to an Airport and Terrain Search Algorithm, indicated by a block 29, which includes location search logic for determining the terrain data, as well as the airport data surrounding the aircraft. Such search logic is described in detail in <patcit id="pcit0041" dnum="US4675823A"><text>U.S. Patent Nos. 4,675,823</text></patcit> and <patcit id="pcit0042" dnum="US4914436A"><text>4,914,436</text></patcit> assigned to the same assignee as the present invention. The GPS inputs, along with terrain and airport data surrounding the aircraft from the search algorithm, indicated by the block 29, are applied to the LOOK-AHEAD warning generator 30, which provides both terrain advisory and terrain warning signals based upon the position and projected flight path of the aircraft. The LOOK-AHEAD warning generator 30 may provide both an aural warning by way of a voice warning generator 32 and speaker 34, and/or a visual warning by way of a map or a display 36, as discussed below.

GLOBAL POSITIONING SYSTEM

0018The primary positional information for the aircraft is provided by the GPS 22, such as disclosed in <patcit id="pcit0043" dnum="US4894655A"><text>U.S. Patent Nos. 4,894,655</text></patcit>; <patcit id="pcit0044" dnum="US4903212A"><text>4,903,212</text></patcit>; <patcit id="pcit0045" dnum="US4912645A"><text>4,912,645</text></patcit>; <patcit id="pcit0046" dnum="US4954959A"><text>4,954,959</text></patcit>; <patcit id="pcit0047" dnum="US5155688A"><text>5,155,688</text></patcit>; <patcit id="pcit0048" dnum="US5257195A"><text>5,257,195</text></patcit>; <patcit id="pcit0049" dnum="US5265025A"><text>5,265,025</text></patcit>; <patcit id="pcit0050" dnum="US5293163A"><text>5,293,163</text></patcit>; <patcit id="pcit0051" dnum="US5293318A"><text>5,293.318</text></patcit>; and <patcit id="pcit0052" dnum="US5337242A"><text>5,337,242</text></patcit>. The GPS 22 includes a GPS receiver 37 and a combination GPS receiver and monitor 38. The positional inputs from the GPS 22 are processed by a dead reckoning algorithm, indicated by the block 40, for example as discussed in detail in <patcit id="pcit0053" dnum="US5257195A"><text>U.S. Patent No. 5,257,195</text></patcit>, to compensate for temporary outages of the GPS 22, such as loss or masking of satellites. The dead reckoning algorithm can be replaced by an additional navigational system, such as an FMS or an IRS, which can be used during GPS outages.

0019A differential GPS receiver 42 can also be used to provide an indication of a GPS altitude. As long as the differential GPS information is received and four more satellites are visible, the differential GPS information from the differential GPS receiver 42 is sufficiently accurate for the TAS 20. If the differential GPS receiver 42 is not available, a compound altitude signal may be generated by the block 44 and connected to the LOOK-AHEAD warning generator 30 by way of a single-pole, double-throw switch 46, under the control of a signal "DIFF GPS AVAILABLE", which indicates when the differential GPS receiver 42 is available. In a first mode when differential GPS information is available, the switch 46 connects the GPS altitude signal from the GPS 22 to the LOOK-AHEAD warning generator 30. In this position, the GPS altitude signal is generated either from the GPS 22 or from the differential GPS receiver 42, as discussed above. When a differential GPS receiver 42 is not available, a compound altitude signal is provided by the block 44 to the LOOK-AHEAD warning generator 30. During this condition, the switch 46, under the control of a "DIFF GPS AVAILABLE" signal, connects the block 44 to the LOOK-AHEAD warning generator 30 to enable the compound altitude signal to be provided to the LOOK-AHEAD warning generator 30. Within the Differential GPC Receiver 42, the barometric altitude is compared with the GPS altitude to generate a compound altitude signal which could be used to limit the maximum difference between the GPS output altitude signal and the barometric altitude to the maximum expected GPS altitude error. Such a configuration would reduce the effect of erroneous pressure correction settings on the barometric altimeter 28.

TERRAIN AND AIRPORT DATA BASES

0020As mentioned above, the airport data and terrain data are separated into two different data bases 26 and 24, respectively. Such a configuration allows either of the data bases 24 or 26 to be updated without the need to update the other.

0021The airport data base 26 contains various types of information relating to airports, such as runway midpoint coordinates, runway length, runway heading, runway elevation and airport/runway designation data. Additional data relating to obstacles along the approach path (i.e., high-rise hotels adjacent Heathrow Airport in England) and a nominal final approach slope could also be included in the airport data base 26. A separate obstacle data base may also be provided, structured similar to the airport data base.

0022In order for the TAS 20 to be useful, all airports, for example, which commuter planes without radio altimeters can land, would have to be included in the airport data base 26. In areas where the airport data base 26 is not complete, the search algorithm 29 would indicate NO DATA. Since the airport data base 26 is organized in a similar fashion as the terrain data base 24 as described below, additional airport data could be added in a piecewise fashion when available. The estimated size of the airport data base 26 is contemplated to be relatively small, compared to the terrain data base 24, for example less than 200 kilobytes depending on the amount of data stored per airport.

0023The terrain data base 24 may require up to 40 megabytes of storage space, which makes it compatible with available flash erasable programmable read only memory (ROM) devices, such as a flash EEPROM, or mini hard disk drives, which could take advantage of known data compression techniques to reduce the amount of storage space required and also promote easy data retrieval by the search algorithm 29.

0024An important aspect of the invention is that the terrain data base 24 is structured to provide varying resolutions of terrain data as a function of the topography of the terrain, as well as distance to airports. For example, a relatively high resolution can be provided close to the airport on the order of 1/4 to 1/8 nautical miles and a medium resolution, for example 1/2 to 1 nautical miles within a 30-mile radius of the airport. Outside of the 30-mile radius from the airport, a coarser resolution is sufficient as will be described below.

0025<figref idref="f0003 f0004 f0005">FIGS. 2 through 4</figref> illustrate the organization of the terrain data base. Referring first to <figref idref="f0003">FIG. 2</figref>, the world is divided into a plurality of latitude bands 50, for example each about 4' wide. Each latitude band 50 is then divided into a plurality of longitudinal segments 52, which are about 4' wide around the equator, such that each longitudinal segment 52 is about 256 x 256 nautical miles. In order to maintain a relatively constant segment size, the number of longitudinal segments 52 per latitude band 50 are reduced closer to the poles.

0026Since the number of latitude bands 50 and the number of longitudinal segments 52 per latitude band 50 is fixed, determination of the particular segment corresponding to the current aircraft position is readily determined. For example, the aircraft current latitude X is used to determine the latitude band number. Next, the number of longitudinal segments 52 associated with that band 50 is determined either by way of a LOOK-UP table or by calculation. Once the number of longitudinal segments 52 in the particular latitude band 50 of interest is determined, the current longitudinal position Y of the aircraft can be used to readily determine the longitudinal segment number.

0027The minimum data associated with each longitudinal segment 52 corresponds to the highest altitude within that segment. As mentioned above, each of the segments 52 is approximately 256 x 256 nautical miles. As shown in <figref idref="f0004">FIG. 3</figref>, these longitudinal segments 52 can be broken down into various subsquares to provide varying levels of resolution. For example, each segment 52 may be broken down into a plurality of subsquares 54, each subsquare 54 being 64 x 64 nautical miles to provide a very coarse resolution. The subsquares 54, can, in turn, be further subdivided into a number of subsquares 56, for example, each 16 x 16 nautical miles, to provide a coarse resolution. These subsquares 56, in turn, can be subdivided into a plurality of subsquares 58, for example, each 4 x 4 nautical miles, to provide a medium resolution. The subsquares 58 can also be broken down into a plurality of subsquares 60, for example 1 x 1 nautical miles, to provide a fine resolution. In the vicinity around airports, it may even be desirable to break down the subsquares 60 down into smaller subsquares 62 to provide even finer resolution, for example 1/4 x 1/4 nautical miles.

0028As shown in <figref idref="f0003">FIG. 4A</figref>, the minimum data associated with each longitudinal segment 52 consists of a header 65 which includes a multiple data byte 66 which includes the reference altitude which corresponds to the highest altitude for all the subsquares within the segment, which assumes that the highest elevation in a square is representative of the entire square elevation. The multiple data byte 66 may also include a flag to indicate when no further subdividing is required for certain geographical areas. For example, for segments representing the ocean, all subsquares would have the same maximum altitude and thus no further subdividing would be required. In order to enable the data base to be created and updated on a piecewise basis, the multiple data byte 66 may also contain a flag bit or code indicating that data corresponding to further subdivisions does not exist for a particular segment 52 containing a code indicating that no map data exists in the segment.

0029For geographical areas, such as mountainous areas and areas in the vicinity of an airport, the longitudinal segments 52 are subdivided as discussed above. In such a situation, the stored data would include a 2-byte mask word 68, which points to the subsquares having different altitudes. For example, as shown in <figref idref="f0003 f0005">FIG. 4</figref>, the 2-byte mask word 68 is illustrated with a pointer in the box 15 within the subsquare 54 within the longitudinal segment 52, which represents a different altitude which points to the next finer layer of resolution as discussed above. The subsquares each include a header 70 including a mask as described above to enable a pointer for each subsquare having a different altitude to point to the next finer layer until the finest resolution level desire is reached as shown.

0030The data base structure provides for flexibility as to resolution required and can be constructed on a piecewise basis at different times. As such, the structure contains the entire framework necessary to add new data or modify existing data without changing any of the software.

0031The size of the data base is proportional to the resolution and the size of the area to be covered. In order to cover the entire earth's surface, about 149 million square nautical miles, only about 2,500 to 3,500 longitudinal segments 52 are required. Thus, the overhead for the minimum header sizes is relatively small.

0032Alternatively, the world may be divided into a plurality of, for example, 1° x 1° map files or cells 51 as shown in <figref idref="f0005">FIG. 4B</figref>. The various map files 51 may be formed with variable resolution. More particularly, each of the map files may be subdivided into a plurality of subcells 53 with the highest terrain elevation in each subcell indicated, with the resolution of all of the subcells within a particular map file being of the same resolution. As indicated above, the maps are generally 1° x 1°, which may be combined for coarser resolution (i.e., 2° x 2° or 10° x 10°).

0033The size of the subcells 53 is varied to provide varying degrees of resolution. For example, for relatively high resolution, the size of the subcells 53 may be selected as either 15 x 15 arc seconds or 30 x 30 arc seconds. The size of the subcells 53 may also be varied as a function of the terrain in particular geographic areas. For example, the size of the subcells 53 may be selected as 30 x 60 arc seconds at relatively far north longitudes and as 30 x 120 arc seconds at longitudes even further north. The size of the subcells 53 may also be varied as a function of the distance to the nearest airport. For example, for distances greater than 64 miles from an airport, the size of the subcells may be selected as 60 x 60 arc seconds while at distances greater than, for example, 128 miles from the airport, size of the subcell 53 may be selected as 120 x 120 arc seconds. For a course resolution, for example for use during a cruise mode as discussed below, the size of subcells may be selected as 5 x 5 arc minutes.

0034As discussed below and as illustrated in <figref idref="f0020">FIG. 27</figref>, the data base may be compressed by a compression algorithm and stored in a flash read only memory (ROM), configured as a wedding cake as illustrated in <figref idref="f0021">FIG. 28</figref>. Various compression algorithms are known and suitable for this application. The compression algorithm does not form a part of the present invention.

0035In order to simplify computation, each of the map files includes a header 55, which includes, <i>inter alia,</i> the resolution of the particular map file 51 as well as the location of one or more corners, for example, the northwest and southeast corners. The headers 55 as well as the map files 51 are discussed in more detail below and illustrated in more detail in <figref idref="f0022">FIG. 29</figref>.

0036For certain terrain, for example, terrain adjacent a runway, the resolution of the map files 51 (i.e., relatively large size of the subcells 53) may result in errors in runway elevation. More particularly, since each subcell 53 contains the highest elevation within the subcell 53, certain errors relating to runway elevations can occur, depending on the resolution of the map file 51. Correction of the elevations adjacent the runways is discussed in more detail below in conjunction with <figref idref="f0022 f0023">FIGS. 30-33</figref>. Systematic errors may also be introduced by the maps, which are digitized, for example, by the Defense Mapping Agency, using known averaging algorithms, which are known to fill in valleys and cut-off the tops of sharp peaks of the terrain.

LOOK-AHEAD WARNING GENERATOR

0037As mentioned above, the LOOK-AHEAD warning generator 30 generates both a terrain advisory signal and a terrain warning signal based upon the position and trajectory of the aircraft relative to stored terrain data. There are two aspects of the terrain advisory and terrain warning signals: LOOK-AHEAD distance/direction; and terrain threat boundaries.

LOOK-AHEAD DISTANCE/DIRECTION

0038The LOOK-AHEAD direction for detecting threatening terrain is along the groundtrack of the aircraft. In order to prevent nuisance warnings, the LOOK-AHEAD distance is limited, as will be discussed below. Otherwise, potentially threatening terrain along the current flight path of the aircraft relatively far from its current position could produce nuisance warnings.

0039Two different LOOK-AHEAD distances (LAD) are utilized. The first LAD is used for a terrain advisory signal (also referred to as a yellow alert LAD). A second LAD is used for terrain warning signals which require immediate evasive action (also referred to as a red-alert LAD).

0040The LAD for a terrain advisory condition is considered first in determining the LAD because it is assumed that the pilot could make a 30° bank turn at any time at a turning radius R. The LAD is equal to the product of the speed of aircraft and the total LOOK-AHEAD time. As shown in <figref idref="f0006">FIG. 5</figref>, the total LOOK-AHEAD time is equal to the sum of the LOOK-AHEAD time T₁ for a single turning radius R; the LOOK-AHEAD time T₂ for terrain clearance at the top of the turn (i.e., point 73) plus a predetermined reaction T₃.

0041The terrain clearance at the top 73 of the turn is provided to prevent inadvertent terrain contact as a result of the turn. This terrain clearance may be selected as a fixed distance or may be made equal to the turning radius R.

0042As shown in equation (1), the turning radius R is proportional to the square of the speed of the aircraft and inversely proportional to the bank angle TG (ROLL). <maths id="math0001" num="(1)"><math display="block"><mi mathvariant="normal">R</mi><mo mathvariant="normal">=</mo><msup><mi mathvariant="normal">V</mi><mn mathvariant="normal">2</mn></msup><mo mathvariant="normal">/</mo><mi mathvariant="normal">G</mi><mo mathvariant="normal">×</mo><mi>TG</mi><mfenced><mi>Roll</mi></mfenced></math><img file="EP0925518B1_D0001.tif" /></maths>

0043For a 30° bank angle, the turning radius R in nautical miles (nM) as a function of the speed in knots is given by equation (2). <maths id="math0002" num="(2)"><math display="block"><mi mathvariant="normal">R</mi><mo mathvariant="normal">=</mo><mn mathvariant="normal">0.000025284</mn><mo mathvariant="normal">*</mo><msup><mi mathvariant="normal">V</mi><mn mathvariant="normal">2</mn></msup></math><img file="EP0925518B1_D0002.tif" /></maths>

0044The LOOK-AHEAD time T ₁ for a single turning radius is given by equation (3). <maths id="math0003" num="(3)"><math display="block"><msub><mi mathvariant="normal">T</mi><mn mathvariant="normal">1</mn></msub><mo mathvariant="normal">=</mo><mi mathvariant="normal">R</mi><mo mathvariant="normal">/</mo><mi mathvariant="normal">V</mi></math><img file="EP0925518B1_D0003.tif" /></maths>

0045Substituting R from equation (1) into equation (3) yields the LOOK-AHEAD time T₁ for a single turn radius as shown in equation (4) as a function of the speed of the aircraft and the bank angle. <maths id="math0004" num="(4)"><math display="block"><mi mathvariant="normal">T</mi><mo mathvariant="normal">=</mo><mi mathvariant="normal">V</mi><mo mathvariant="normal">/</mo><mi mathvariant="normal">G</mi><mo mathvariant="normal">×</mo><mi>TG</mi><mfenced><mi>Roll</mi></mfenced></math><img file="EP0925518B1_D0004.tif" /></maths>

0046TABLE 1 provides various LOOK-AHEAD times T₁ at various turning radii and ground speeds. <tables id="tabl0001" num="0001"><table frame="all"><title><u style="single">TABLE 1</u></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="16mm" /><colspec colnum="2" colname="col2" colwidth="16mm" /><colspec colnum="3" colname="col3" colwidth="13mm" /><thead><row rowsep="0"><entry align="center" valign="top"><b>Radius</b></entry><entry align="center" valign="top"><b>Speed</b></entry><entry align="center" valign="top"><b>T₁</b></entry></row><row><entry align="center" valign="top"><b>(nM)</b></entry><entry align="center" valign="top"><b>(knots)</b></entry><entry align="center" valign="top"><b>(sec)</b></entry></row></thead><tbody><row><entry align="center">1/2</entry><entry align="center">140</entry><entry align="center">13</entry></row><row><entry align="center">1</entry><entry align="center">200</entry><entry align="center">18</entry></row><row><entry align="center">2</entry><entry align="center">280</entry><entry align="center">26</entry></row><row><entry align="center">3</entry><entry align="center">345</entry><entry align="center">31</entry></row><row><entry align="center">4</entry><entry align="center">400</entry><entry align="center">38</entry></row><row><entry align="center">5</entry><entry align="center">445</entry><entry align="center">40</entry></row><row><entry align="center">6</entry><entry align="center">490</entry><entry align="center">44</entry></row></tbody></tgroup></table></tables>

0047By assuming that the fixed terrain clearance at the top 73 of the turn is equal to one turning radius R, the total LOOK-AHEAD time for two turn radii (i.e., T₁ + T₂) is simply twice the time for a single turning radius. Thus, the total LOOK-AHEAD time is 2 * T₁ plus the predetermined reaction time T₃, for example 10 seconds, as given by equation (5). <maths id="math0005" num="(5)"><math display="block"><msub><mi mathvariant="normal">T</mi><mi>TOTAL</mi></msub><mo mathvariant="normal">=</mo><mn mathvariant="normal">2</mn><mo mathvariant="normal">*</mo><msub><mi mathvariant="normal">T</mi><mn mathvariant="normal">1</mn></msub><mo mathvariant="normal">+</mo><msub><mi mathvariant="normal">T</mi><mn mathvariant="normal">3</mn></msub></math><img file="EP0925518B1_D0005.tif" /></maths>

0048By substituting R from equation (2) into equation (4), the LOOK-AHEAD time for a single turn radius is given by equation (6). <maths id="math0006" num="(6)"><math display="block"><msub><mi mathvariant="normal">T</mi><mn>1</mn></msub><mo mathvariant="normal">=</mo><mn mathvariant="normal">0.000025284</mn><mo mathvariant="normal">*</mo><mi mathvariant="normal">V</mi></math><img file="EP0925518B1_D0006.tif" /></maths>

0049Thus, if it is assumed that the fixed clearance X at the top of the turn 73 is also equal to the turning radius R, the total LOOK-AHEAD time for two turning radii is given by equation (7) below, simply twice the value determined in equation (6). <maths id="math0007" num="(7)"><math display="block"><msub><mi mathvariant="normal">T</mi><mn>1</mn></msub><mo>+</mo><msub><mi mathvariant="normal">T</mi><mn>2</mn></msub><mo mathvariant="normal">=</mo><mn mathvariant="normal">2</mn><mo mathvariant="normal">*</mo><mn mathvariant="normal">0.000025284</mn><mo mathvariant="normal">*</mo><mi mathvariant="normal">V</mi></math><img file="EP0925518B1_D0007.tif" /></maths>

0050For a predetermined reaction time T₃, for example 10 seconds, the LOOK-AHEAD distance (LAD) in nautical miles for a terrain advisory signal can be determined simply by multiplying the speed of the aircraft (V) by the total time T_TOTAL as shown in equation (8). <maths id="math0008" num="(8)"><math display="block"><mi>LAD</mi><mo mathvariant="normal">=</mo><mi mathvariant="normal">V</mi><mo mathvariant="normal">*</mo><mfenced><mn mathvariant="normal">0.0000278</mn><mo mathvariant="normal">+</mo><mn mathvariant="normal">2</mn><mo mathvariant="normal">*</mo><mn mathvariant="normal">0.000252854</mn><mo mathvariant="normal">*</mo><mi mathvariant="normal">V</mi></mfenced><mo mathvariant="normal">+</mo><mi>fixed distance</mi><mo mathvariant="normal">*</mo><mi mathvariant="normal">K</mi><mo mathvariant="normal">,</mo></math><img file="EP0925518B1_D0008.tif" /></maths> where K is a constant, for example 0.

0051However, in order to provide optimal conditions during a terrain advisory condition, the LAD is limited with an upper limit and a lower limit. The lower limit may be a configurable amount, for example, either 0.75, 1 or 1-1/2 nautical miles at relatively low speeds (i.e., speeds less than 150 knots) and limited to, for example, 4 nautical miles at speeds, for example, greater than 250 knots. The LAD may also be limited to a fixed amount regardless of the speed when the distance to the runway is less than a predetermined amount, for example 2 nautical miles, except when the aircraft altitude is greater than 3,000 feet, relative to the runway.

0052The terrain warning LAD is given by equation (9) below. <maths id="math0009" num="(9)"><math display="block"><mtable columnalign="left"><mtr><mtd><mi>LAD</mi><mo mathvariant="normal">=</mo></mtd></mtr><mtr><mtd><mspace width="10em" /><mtable columnalign="left"><mtr><mtd><msub><mi mathvariant="normal">k</mi><mn mathvariant="normal">1</mn></msub><mo mathvariant="normal">*</mo><mi>LAD</mi><mrow><mo mathvariant="normal">(</mo><mi>terrain for LOOK</mi><mo mathvariant="normal">-</mo><mi>DOWN</mi><mo mathvariant="normal">/</mo><mi>LOOK</mi><mo mathvariant="normal">-</mo><mi>AHEAD advisory indication</mi><mo mathvariant="normal">)</mo><mo mathvariant="normal">,</mo></mrow></mtd></mtr><mtr><mtd><msub><mi mathvariant="normal">k</mi><mn mathvariant="normal">2</mn></msub><mo mathvariant="normal">*</mo><mi>LAD</mi><mrow><mo mathvariant="normal">(</mo><mi>terrain for LOOK</mi><mo mathvariant="normal">-</mo><mi>DOWN</mi><mo mathvariant="normal">/</mo><mi>LOOK</mi><mo mathvariant="normal">-</mo><mi>AHEAD warning</mi><mo mathvariant="normal">)</mo><mo mathvariant="normal">,</mo></mrow></mtd></mtr><mtr><mtd><msub><mi mathvariant="normal">k</mi><mn mathvariant="normal">3</mn></msub><mo mathvariant="normal">*</mo><mi>LAD</mi><mo>⁢</mo><mfenced><mi>terrain for LOOK</mi><mo mathvariant="normal">-</mo><mi>UP advisory</mi></mfenced><mo mathvariant="normal">,</mo></mtd></mtr></mtable></mtd></mtr></mtable></math><img file="EP0925518B1_D0009.tif" /></maths> where k₁ = 1.5, except when the LAD is limited at its lower limit, in which case k₁ = 1, k₃ = 1 and where k₂ = 0.5, k₃ = 2.

TERRAIN THREAT BOUNDARIES

0053The terrain threat boundaries include a terrain floor boundary, terrain advisory boundaries (yellow alert), as well as terrain warning boundaries (red alert) and are provided along the groundtrack of the aircraft.

TERRAIN FLOOR BOUNDARY

0054The terrain floor boundary is the basis for the terrain threat boundaries and is similar to the terrain floor developed for the GPWS. The terrain floor relates to a distance DH below the aircraft and is proportional to the distance to the closest runway to prevent nuisance warnings when the aircraft is taking off and landing, while providing adequate protection in other modes of operation. As illustrated in equation (10), the terrain floor boundary below the aircraft is essentially based upon providing 100 feet clearance per nautical mile from the runway, limited to 800 feet. <maths id="math0010" num="(10)"><math display="block"><mi>DH</mi><mo>=</mo><mn>100</mn><mrow><mo>(</mo><mi>ft</mi><mo>/</mo><mi>nM</mi><mo>)</mo><mo>*</mo><mrow><mo>(</mo><mi>distance to runway threshold</mi><mo>-</mo><mi>an offset</mi><mo>)</mo><mo>+</mo><mn>100</mn><mrow><mo>(</mo><mi>ft</mi><mo>/</mo><mi>nM</mi><mo>)</mo><mo>*</mo><mfenced><mi>distance to runway center</mi><mo>-</mo><mn>12</mn><mspace width="1em" /><mi>nautical miles</mi></mfenced></mrow></mrow></mrow></math><img file="EP0925518B1_D0010.tif" /></maths>

0055Equation (10) is shown graphically in <figref idref="f0006">FIG. 6</figref>. Referring to <figref idref="f0006">FIG. 6</figref>, the horizontal axis represents the distance from the runway, while the vertical axis represents the DH terrain floor boundary beneath the aircraft. The first segment 71 represents the runway length, while the point 72 represents the runway center. After a small offset 74, the next segment 78, which starts at 0 feet, slopes at 100 feet per nautical mile distance from the runway threshold, identified by the point 80, up to a maximum of 500 feet. The 500 foot maximum continues along a segment 82 until the distance to the runway center 72 is a predetermined distance D, for example 12 nautical miles. The next segment 84 slopes up at a 100 feet per nautical mile distance until the distance rises 300 feet from the segment 82 for a total of 800 feet.

0056The terrain floor DH boundary beneath the aircraft is limited such that the segment 78 begins at 0 and the segment 82 never goes above a predetermined maximum, for example 500 feet. In addition, the terrain floor DH boundary is also limited so that the segment 84 never goes below the vertical height of the segment 82 and rises no more than 300 feet relative to the segment 82. During conditions when the aircraft is 3,000 feet above the terrain, the terrain floor DH boundary is maintained at the upper limit of 800 feet.

TERRAIN ADVISORY BOUNDARIES

0057Two terrain advisory boundaries are shown graphically in <figref idref="f0007">FIGS. 7 and 8. FIG. 7</figref> represents a condition when the aircraft is descending while <figref idref="f0007">FIG. 8</figref> represents a condition when the aircraft is ascending. As will be discussed in more detail below, the terrain advisory boundaries are based upon the relationship between the flight path angle g and a first configurable datum, THETA1.

0058Referring first to <figref idref="f0007">FIG. 7</figref>, the terrain advisory (yellow alert) boundaries are shown. The first segment of the terrain advisory boundary, identified with the reference numeral 92, corresponds to the DH terrain floor boundary. As mentioned above, the DH terrain floor boundary is a function of the distance from the runway 94. In order to determine the bottom segment 100 of the terrain advisory boundary, it is necessary to determine whether the flight path angle g is more or less than a configurable datum, identified as THETA1. As shown and described herein, THETA1 is set for 0. However, it should be clear that other angles for THETA1 are within the broad scope of the present invention. For the condition illustrated in <figref idref="f0007">FIG. 7</figref>, the flight path angle g is less than THETA1. Thus, the segment 100 extends from the DH terrain floor boundary at the angle THETA1 to the LOOK-AHEAD distance for a terrain advisory. The final segment 102 is generally parallel to the segment 92 and extends in a vertical direction along the LAD for the terrain advisory boundary.

0059If the flight path angle g is greater than THETA 1 as shown in <figref idref="f0007">FIG 8</figref>., then different terrain advisory boundaries are provided. The segment 92 will be the same as shown in <figref idref="f0007">FIG. 7</figref> and described above. However, when the flight path angle g is greater than THETA1 as shown in <figref idref="f0007">FIG. 8</figref>, a bottom segment 106 is formed by extending a line segment from the bottom of the segment 92, which represents the DH terrain floor, along a direction parallel with the flight path angle g up to the LAD. A vertical segment 108 extends from the bottom segment 106 along the LAD to form the terrain advisory boundaries.

0060Thus, as shown, if the airplane is climbing, only terrain which penetrates the upward-sloping terrain advisory boundaries will cause an advisory signal, such that if the aircraft clears the terrain by a safe margin, no advisory signal is given, even if the terrain ahead is at or above the present aircraft altitude. However, if terrain penetrates the terrain advisory boundaries, an aural warning, such as "CAUTION, TERRAIN" aural advisory may be given. In addition, visual warnings may be provided.

TERRAIN WARNING BOUNDARIES

0061The terrain warning boundaries are shown in <figref idref="f0008">FIGS. 9 and 10</figref> and, in general, indicate to the pilot of an aircraft conditions when evasive action is required to avoid terrain contact. The terrain warning boundaries are based on the relationship of the flight path angle of the aircraft relative to a second configurable datum, THETA2. The datum THETA2 may be selected with an upslope of, for example, 6°, which is equal to the average climb capability of airliners. The datum THETA2 could be modified, taking into consideration aircraft type, configuration, altitude and time for takeoff. However, because of the longer LOOK-AHEAD distance as discussed above for the terrain warning boundaries, a terrain warning could occur before a terrain advisory indication for extreme terrain conditions.

0062Referring to <figref idref="f0008">FIGS. 9 and 10</figref>, two different terrain warning boundaries are illustrated. In particular, <figref idref="f0008">FIG. 9</figref> illustrates a condition when the flight path angle g is greater than the second datum THETA2. <figref idref="f0008">FIG. 10</figref> illustrates a condition when the flight path angle g is less than the second datum THETA2. In both conditions, the first segment 114 of the terrain warning boundary is extended below the aircraft for a distance 1/2 of the DH terrain floor, which, as mentioned above, is a function of the distance of the aircraft from a runway. As mentioned above, the terrain warning boundaries are dependent on the relationship between the flight path angle g of the aircraft and the datum THETA2.

0063Referring first to <figref idref="f0008">FIG. 9</figref>, since the flight path g is less than the second datum THETA2, the bottom segment 118 of the terrain warning boundary is formed by extending a line from the segment 114 at an angle equal to the flight path angle g up to the LOOK-AHEAD distance for a terrain warning as discussed above. A vertical segment 120 extends from the lower segment 118 along the LAD for a terrain warning, forming the terrain warning boundaries illustrated in <figref idref="f0008">FIG. 9</figref>.

0064If the flight path angle of the aircraft is less than the slope of the second datum THETA2 as shown in <figref idref="f0008">FIG. 10</figref>, a bottom segment 122 extends from the segment 114 up to the LOOK-AHEAD distance at an angle equal to the flight path angle g. A segment 124 extends upwardly along the LAD from the segment 122. If any terrain along the groundtrack of the aircraft out to the LAD for the terrain warning penetrates the envelope, aural and/or visual terrain ahead warnings may be given.

CUT-OFF BOUNDARIES

0065In order to avoid spurious warnings when the aircraft overflies a ridge at relatively low altitudes, the warning boundaries may include cut-off boundaries, for example, as illustrated in <figref idref="f0009">FIGS. 11, 12 and 13</figref>. Without the cut-off boundaries, warnings would be given, although the terrain is virtually below the aircraft and no terrain is visible ahead. Referring first to <figref idref="f0009">FIG. 11</figref>, the cut-off boundary 126 begins at a predetermined cut-off offset 128 below the aircraft and extends in a direction in front of the aircraft at a predetermined envelope cut-off angle 130. The envelope cut-off angle 130 is equal to the flight path angle g plus a configurable predetermined cut-off angle, described and illustrated as -6. For level flight as shown in <figref idref="f0009">FIG. 11</figref>, the cut-off boundary 126 extends from the cut-off offset 128 in the direction of the envelope cut-off angle 130 toward the front of the aircraft to a point 132 where it intersects a terrain advisory boundary or terrain warning boundary, identified with the reference numeral 134. For level flight, as shown in <figref idref="f0009">FIG. 11</figref>, the flight path angle g is zero. Thus, the cut-off boundary 126 illustrated in <figref idref="f0009">FIG. 11</figref> will extend from the cut-off offset 128 along an angle equal to the cut-off angle, which, as mentioned above, is selected as -6° for illustration. As mentioned above, the cut-off boundary 126 extends from the cut-off offset 128 to the point 132 where it intersects the terrain advisory boundary 134 discussed above. The warning boundary is then selected to be the highest of the terrain advisory boundary 134 and the envelope cut-off boundary 126. Thus, for the example illustrated in <figref idref="f0009">FIG. 11</figref>, the terrain advisory boundary would consist of the cut-off boundary 126 up to the point 132, where the envelope cut-off boundary 124 intersects the warning envelope 126. From the point 132 forward, the normal terrain advisory boundary 134, corresponding, for example, to a THETA1 slope, is utilized. Thus, if either a terrain advisory boundary or terrain warning boundary is below the cut-off boundary 126, the cut-off boundary 126 becomes the new boundary for the advisory or warning signal.

0066The cut-off angle 130 is limited such that if the flight path angle g is a relatively high climb angle, the cut-off angle 130 does not exceed a predetermined limit, for example a configurable limit between 0 to 6° as will be illustrated in <figref idref="f0009">FIG. 12</figref>. In particular, referring to <figref idref="f0009">FIG. 12</figref>, the flight path angle g is shown at a relatively high climb rate (i.e., g = 9°). An unlimited cut-off envelope angle 130 is illustrated by the partial segment 136. As shown, this segment 132 is 6° below the flight path angle g However, in order to decrease the sensitivity of the system 20 during conditions when the flight path angle is relatively large, for example 9°, indicative of a fairly steep climb, the cut-off envelope angle 130 is limited. In the example illustrated in <figref idref="f0009">FIG. 12</figref>, the cut-off angle limit is selected, for example, to be 1.4°. Thus, at flight path angles g greater than 7.4°, the cut-off angle will result in an envelope cut-off boundary 140, which is more than 6° below the flight path angle g. As shown, since the cut-off angle limit was selected as 1.4°, the cut-off angle will be limited to 1.4° anytime the flight path angle is greater than 7.4°, even though the cut-off angle will be equal to or more than 6° less than the flight path angle g. Thus, for the conditions illustrated in <figref idref="f0009">FIG. 12</figref>, the segment 140 begins at the cut-off offset 128 and continues along a line that is 1.4° from the horizon until it intersects the terrain advisory or terrain warning boundary 142 at the point 144. The actual advisory or warning boundary will be formed by the limited cut-off boundary 140 up to the point 144 and will continue with the normal warning boundary 142 beyond the point 144.

0067<figref idref="f0009">FIG. 13</figref> illustrates a situation when the aircraft is descending and thus the flight path angle g is less than 0. During conditions when the flight path angle g is less than the normal descent angle, the warning boundary is lowered past the cut-off boundary for high sink-rate situations in order to prevent nuisance warnings during step-down approaches. In particular, a BETA sink rate enhancement is provided and is effective anytime the flight path g is below a predetermined configurable descent bias GBBIAS, for example, -4°. During such conditions, an angle BETA is added to the terrain advisory datum (THETA1) as indicated below. In particular, when the flight path angle g is greater than the descent bias angle GBBIAS or the flight path angle g is less than zero, BETA is selected as 0. When the flight path angle g is greater than zero, BETA is given by equation (11). <maths id="math0011" num="(11)"><math display="block"><mi>BETA</mi><mo mathvariant="normal">=</mo><mi mathvariant="normal">g</mi><mo mathvariant="normal">-</mo><mi>THETA</mi><mo>⁢</mo><mn mathvariant="normal">1</mn><mo mathvariant="normal">;</mo><mi>for g</mi><mo mathvariant="normal">></mo><mn mathvariant="normal">0</mn></math><img file="EP0925518B1_D0011.tif" /></maths>

0068When the flight path angle g is less than GBBIAS, the value for BETA is selected in accordance with equation (12) below. <maths id="math0012" num="(12)"><math display="block"><mi>BETA</mi><mo mathvariant="normal">=</mo><mi mathvariant="normal">K</mi><mo mathvariant="normal">*</mo><mfenced><mi mathvariant="normal">g</mi><mo mathvariant="normal">-</mo><mi>GBBIAS</mi></mfenced><mo>,</mo></math><img file="EP0925518B1_D0012.tif" /></maths> where k = 0.5, g = flight path angle and GBBIAS is selected as -4°.

0069<figref idref="f0009">FIG. 13</figref> illustrates a condition when the flight path angle g is -8°. During this condition, the cut-off angle 130 is selected to be 6° below the flight path angle g and thus defines an envelope cut-off boundary 150 that extends from the cut-off offset 128 along a line which is -14° from the horizontal axis. Since the flight path angle g is less than GBBIAS, the angle BETA is added to the normal warning boundary according to equation (12) above. Thus, for a flight path angle of -8°, the BETA angle will be 2°. Thus, the terrain advisory or warning boundary, indicated by the boundary 151, instead of extending along the THETA1 slope, will extend along an angle which is 2° below the THETA1 slope to define the enhanced boundary 152. The warning cut-off boundary 152 thus will extend to the point 156 where it intersects the normal terrain advisory or terrain warning boundary 152. Since the cut-off angle boundary 150 is higher than boundary 152 up to the point 156, the cut-off boundary 150 will be utilized as the warning boundary up to the point 156. Beyond point 156, the BETA-enhanced warning boundary 152 will define the warning envelope.

0070<figref idref="f0010">FIGS. 14, 15 and 16</figref> illustrate the composite terrain advisory and terrain warning boundaries for different conditions. The cross-hatched area illustrates terrain in the vicinity of the aircraft. Three conditions are illustrated. In particular, <figref idref="f0010">FIG. 14</figref> represents a condition when the flight path angle g of the aircraft is less than the first datum THETA1. <figref idref="f0010">FIG. 15</figref> represents a condition when the flight path angle g is between the first datum THETA1 and the second datum THETA2. <figref idref="f0010">FIG. 16</figref> represents a condition when the flight path angle g is greater than the second datum THETA2. The terrain advisory boundaries are indicated in solid line while the terrain warning boundaries are illustrated by dashed lines.

0071Referring first to <figref idref="f0010">FIG. 14</figref>, the terrain advisory boundaries include a boundary 155, which extends from the aircraft to the DH terrain floor. Normally, the lower terrain advisory boundary 157 would extend from the segment 155 along an angle equal to the THETA1 slope. However, in this situation, the flight path angle g is illustrated to be greater than the GBBIAS angle; therefore, a BETA-enhanced boundary 159, as discussed above, is utilized which extends along an angle that is equal to the THETA1 angle plus the BETA angle, up to the LOOK-AHEAD distance for a terrain advisory. A vertical boundary 158 extends from the first datum THETA1 vertically upwardly along the LOOK-AHEAD distance for a terrain advisory. Since the cut-off boundary 157 is higher than the terrain advisory boundary 157 up to the point 169 where the two boundaries intersect, the cut-off boundary 157 becomes the terrain advisory boundary up to the point 169. Beyond the point 169, the boundary 159 is utilized as the warning boundary.

0072The terrain warning boundary includes a boundary 162 extending to a point equal to 1/2 of the terrain floor as discussed above. A boundary 165 extends from the boundary 162 at the THETA2 slope. A vertical boundary 167 extends from the boundary 165 along the LOOK-AHEAD distance for a terrain warning as discussed above.

0073<figref idref="f0010">FIG. 15</figref> illustrates a situation when the flight path angle g is greater than the THETA1 angle and less than or equal to the THETA2 angle. During this condition, the boundary for the terrain warning extends along a segment 180, which is at 1/2 of the DH terrain floor. A sloping boundary 184 extends from the segment 180 to the LOOK-AHEAD distance for a terrain warning (point 186). From there, a vertical boundary 188 extends vertically upwardly to define the terrain warning boundaries.

0074The terrain advisory boundary includes a boundary 185 which extends downwardly to the terrain floor, i.e., point 187. From point 187, the terrain advisory boundary extends along a segment 189 at the flight path angle g up to a point 191, the LOOK-AHEAD distance for a terrain advisory. A vertical segment 194 extends along the LAD from the boundary 189 at the point 191.

0075In this illustration, an envelope cut-off boundary 193 cuts off a portion of the terrain advisory indication boundary since the flight path angle g of the aircraft is greater than a predetermined limit. Thus, the portions of the warning boundary 189 below the cut-off angle boundary 193 are ignored; and the segment 193 becomes the effective terrain advisory indication up to the point 195, where the cut-off boundary 193 intersects the terrain advisory indication boundary 189. After the point 195, the boundary 189 will be the effective warning boundary.

0076<figref idref="f0010">FIG. 16</figref> illustrates a situation when the flight path angle g is greater than THETA2. In this situation, a terrain advisory boundary 200 extends from the aircraft to the boundary terrain floor. A lower boundary 202 extends at the flight path angle g from the boundary 200 up to the LOOK-AHEAD distance for a terrain advisory. A vertical terrain advisory boundary 204 extends upwardly along the LOOK-AHEAD distance for the terrain advisory. Since the flight path angle g is illustrated to the greater than GBBIAS, a segment 211, the cut-off boundary, defines the effective warning boundary up to a point 213 where the terrain advisory boundary intersects the cut-off warming boundary. After that point, segment 202, which is higher, becomes the warning boundary.

0077The terrain warning boundary extends along a segment 208 down to a point of 1/2 the terrain floor. A lower boundary 210 extends along an angle equal to flight path angle g since it is greater than the THETA2 slope. A vertical terrain warning boundary 212 extends upwardly from the boundary 210 along the LOOK-AHEAD distance for a terrain warning.

0078In order to provide additional sensitivity of the warning system, terrain advisory signals within a predetermined portion of the LAD for a terrain advisory, for example 1/2 LAD, may be treated as terrain warnings. Thus, as shown in <figref idref="f0010">FIGS. 14, 15 and 16</figref>, vertical boundaries 161, 181 and 201, respectively, are shown, for example, at 1/2 LAD. Based upon the flight path angles shown, these segments define effective terrain advisory envelopes 166, 182 and 203.

ALTERNATIVE TERRAIN THREAT BOUNDARIES

0079Alternative terrain threat boundaries are illustrated in <figref idref="f0011 f0012 f0013 f0014 f0015">FIGS. 17-22</figref>. Similar to the terrain threat boundary discussed above and illustrated in <figref idref="f0006 f0007 f0008 f0009 f0010">FIGS. 5-16</figref>, the alternative terrain threat boundaries include a terrain floor boundary, terrain advisory boundaries (yellow alert) and terrain warning boundaries (red alert). In the alternative embodiment, the advisory and warning boundaries are further divided into two parts, a LOOK-AHEAD/LOOK-DOWN boundary for detecting terrain ahead or below the aircraft (<figref idref="f0011">FIGS. 17 and 18</figref>) and a LOOK-UP boundary for detecting precipitous high terrain ahead of the aircraft which may be difficult to clear (<figref idref="f0012">FIG. 19</figref>).

LOOK-AHEAD/LOOK-DOWN TERRAIN ADVISORY AND WARNING BOUNDARIES

0080As mentioned above, the LOOK-AHEAD/LOOK-DOWN boundaries are illustrated in <figref idref="f0011">FIGS. 17 and 18</figref>. Referring first to <figref idref="f0011">FIG. 17</figref>, LOOK-AHEAD/LOOK-DOWN advisory and warning boundaries are illustrated for a condition when the aircraft is descending (i.e., g &lt; 0). During such a configuration, the first segment of the LOOK-AHEAD/LOOK-DOWN terrain advisory boundary, identified with the reference numeral 300, corresponds to the DH terrain floor boundary. As discussed above, the DH terrain floor boundary is a function of the aircraft from the runway. In order to determine the bottom segment 302 of the LOOK-AHEAD/LOOK-DOWN terrain advisory boundary, the flight path angle g is compared with a configurable datum, THETA1, for example 0°. During descent conditions, the flight path angle g will thus be less than THETA1. Thus, the LOOK-AHEAD/LOOK-DOWN terrain advisory boundary segment will extend from the DH terrain floor boundary segment 300 along the angle THETA1 to the look-ahead distance for a terrain advisory (LAD). The final segment 304 extends vertically upwardly from the segment 302 along the LAD.

0081The LOOK-AHEAD/LOOK-DOWN terrain advisory boundary may also be modified by the BETA sink rate enhancement as discussed above. In this embodiment, the BETA sink rate enhancement ensures that an advisory indication always precedes a warning indication when the aircraft descends into or on top of terrain. The BETA sink rate enhancement is determined as a function of the flight path angle g and two (2) configurable constants KBETA and GBIAS. The BETA sink rate enhancement BETA1 for a LOOK-AHEAD/LOOK-DOWN terrain advisory is provided in equation (13) below for a condition when the aircraft is descending. <maths id="math0013" num="(13)"><math display="block"><mi>BETA</mi><mo>⁢</mo><mn>1</mn><mo mathvariant="normal">=</mo><mi>KBETA</mi><mo mathvariant="normal">*</mo><mfenced><mi mathvariant="normal">g</mi><mo mathvariant="normal">-</mo><mi>GBIAS</mi></mfenced><mo>,</mo></math><img file="EP0925518B1_D0013.tif" /></maths> where GBIAS is a configurable constant, selected for example, to be zero (0) and KBETA is also a configurable constant selected, for example, to be 0.5.

0082Referring to <figref idref="f0011">FIG. 17</figref>, the BETA sink rate enhancement BETA1 for the LOOK-AHEAD/LOOK-DOWN terrain advisory boundary provides an advisory warning at a distance less than 1/2 LAD, unlike the terrain advisory boundary described above). More particularly, for the values defined above for the configurable constants KBETA and GBIAS, the BETA sink rate enhancement BETA1 results in a segment 306 which extends from the DH terrain floor segments 300 at an angle equal to g/2 up to 1/2 of the LAD. Beyond 1/2 LAD, a segment 308 extends at the angle THETA1 to a distance equal to the LAD. A vertical segment 310 extends along the LAD to connect the segments 308 to the segment 304.

0083As discussed above, a cut-off boundary may be provided to reduce nuisance warnings. Referring to <figref idref="f0011">FIG. 17</figref>, the cut-off boundary is identified with the reference numeral 312. The cut-off boundary 312 extends from a vertical distance 314 below the aircraft along a cut-off angle up to the point of intersection 316 with the terrain advisory boundary. For distances less than the intersection 316, the cut-off boundary 312 forms the terrain advisory boundary. For distances beyond the point of intersection 316, the boundaries 306 and 308 form the terrain advisory boundaries.

0084The terrain warning boundary includes the segment 300 extending from the aircraft along the DH terrain floor. A bottom segment 318 connects to the segment 300 and extends along a BETA sink rate enhancement angle BETA2. The BETA sink rate enhancement angle is determined as a function of the flight path angle g and a configurable constant KBETA2 as well as the constant GBIAS discussed above and provided in equation (14) below. <maths id="math0014" num="(14)"><math display="block"><mi>BETA</mi><mo>⁢</mo><mn>2</mn><mo mathvariant="normal">=</mo><mi>KBETA</mi><mo>⁢</mo><mn>2</mn><mo mathvariant="normal">*</mo><mfenced><mi>GAMMA</mi><mo mathvariant="normal">-</mo><mi>GBIAS</mi></mfenced><mo>,</mo></math><img file="EP0925518B1_D0014.tif" /></maths> where GBIAS is a configurable constant selected, for example, to be 0 and KBETA2 is also a configurable constant selected, for example, to be 0.25.

0085For such values of the constants KBETA2 and GBIAS, the BETA enhancement angle KBETA2 will be 1/4 * g. Thus, the segment 318 extends from the segment 300 at an angle equal to 1/4 * g up to 1/2 the LAD. A vertical segment 320 extends along a distance equal to 1/2 * LAD from the segment 318 to define the terrain warning boundary.

0086Similar to the above, the terrain warning boundaries are also limited by the cut-off boundary 312. Thus, the cut-off boundary 312 forms the terrain warning boundary up to a point 322, where the cut-off boundary 312 intersects the lower terrain warning boundary 318. At distances beyond the point of intersection 322, the segment 318 forms the lower terrain warning boundary up to a distance equal to 1/2 of the LAD.

0087The terrain advisory and terrain warning boundaries for a condition when the aircraft is climbing (i.e., g > 0) is illustrated in <figref idref="f0011">FIG. 18</figref>. During such a condition, the BETA sink rate enhancement angles BETA1 and BETA2 are set to a configurable constant, for example, zero (0).

0088The terrain advisory boundary during a climbing condition is formed by extending a vertical segment 324 from the aircraft for a distance below the aircraft equal to the DH terrain floor. During a climbing condition, a segment 326 is extended from the segment 324 to the LAD at an angle equal to the flight path angle g. At a point 328 where the segment 326 intersects a position equal to 1/2 of the LAD, a vertical segment 330 is extended up from the segment 326, forming a first vertical boundary for the terrain advisory condition. The line segment 326 from the point 328 to the LAD forms the lower terrain advisory boundary while a line segment 332 extending vertically upward from the line segment 326 along the LAD forms a second vertical boundary.

0089For the exemplary condition illustrated, a cut-off boundary 334 does not intersect the terrain advisory boundaries. Thus, the terrain advisory boundaries for the exemplary condition illustrated is formed by the segments 330 and 332 and that portion of the line segment 326 between the line segments 330 and 332.

0090The terrain warning boundaries for a condition when the aircraft is climbing includes the vertical segment 324 (<figref idref="f0011">FIG. 18</figref>) which extends from the aircraft to vertical distance equal to the DH terrain floor below the aircraft forming a first vertical boundary. For a condition when the aircraft is climbing, the line segment 326 extends from the segment 324 at the flight path angle g to form the lower terrain warning boundary. The vertical segment 330 at a distance equal to 1/2 of the LAD forms the second vertical terrain warning boundary.

0091The cut-off boundary 334 limits a portion of the terrain warning boundary. In particular, the cut-off boundary 334 forms the lower terrain warning boundary up to a point 340, where the cut-off boundary 334 intersects the line segment 326. Beyond the point 340, a portion 342 of the line segment 340 forms the balance of the lower terrain warning boundary up to a distance equal to 1/2 of the LAD.

LOOK-UP TERRAIN ADVISORY AND WARNING BOUNDARIES

0092The LOOK-UP terrain advisory and terrain warning boundaries are illustrated in <figref idref="f0012">FIG. 19</figref>. As will be discussed in more detail below, the LOOK-UP terrain advisory and warning boundaries start at altitudes DHYEL2 and DHRED2, respectively, below the aircraft. These altitudes DHYEL2 and DHRED2 are modulated by the instantaneous sink rate (i.e., vertical speed, HDOT of the aircraft). The amount of modulation is equal to the estimated altitude loss for a pull-up maneuver, for example, at 1/4G (i.e., 8 ft/sec²) at the present sink rate. The altitudes DHRED2 and DHYEL2 are dependent upon the altitude loss during a pull-up maneuver (A1PU) and the altitude loss due to reaction time (ALRT) discussed below.

0093The estimated altitude loss due to a pull-up maneuver ALPU is best understood with reference to <figref idref="f0013">FIG. 20</figref>, wherein the vertical axis relates to altitude in feet and the horizontal axis relates to time in seconds. The trajectory of the aircraft from a time shortly before pull-up is initiated to a time when the aircraft has recovered is shown and identified with the reference numeral 344.

0094Assuming an advisory or warning indication is generated at time T₁ at a point 346 along the trajectory 344, the altitude loss due to the reaction time of the pilot (ALRT) is given by equation (15) below. <maths id="math0015" num="(15)"><math display="block"><mi>ALRT</mi><mo>=</mo><mi>HDOT</mi><mo>*</mo><msub><mi mathvariant="normal">T</mi><mi mathvariant="normal">R</mi></msub><mo>,</mo></math><img file="EP0925518B1_D0015.tif" /></maths> where HDOT equals the vertical acceleration of the aircraft in feet/sec and T_R equals the total reaction time of the pilot in seconds.

0095Assuming a pull-up maneuver is initiated at time T₂, the altitude loss due to the pull-up maneuver ALPU may be determined by integrating the vertical velocity HDOT with respect to time as set forth in equation (16) below. <maths id="math0016" num="(16)"><math display="block"><mi>HDOT</mi><mfenced><mi mathvariant="normal">t</mi></mfenced><mo mathvariant="normal">=</mo><mi mathvariant="normal">a</mi><mo mathvariant="normal">*</mo><mi mathvariant="normal">t</mi><mo mathvariant="normal">+</mo><msub><mi>HDOT</mi><mn mathvariant="normal">0</mn></msub><mo mathvariant="normal">,</mo></math><img file="EP0925518B1_D0016.tif" /></maths> where "a" equals the pull-up acceleration and HDOT₀ is a constant.

0096Integrating both sides of equation (16) yields the altitude loss as a function of time H(t) as provided in equation (17) below. <maths id="math0017" num="(17)"><math display="block"><mi mathvariant="normal">H</mi><mfenced><mi mathvariant="normal">t</mi></mfenced><mo mathvariant="normal">=</mo><mn mathvariant="normal">1</mn><mo mathvariant="normal">/</mo><mn mathvariant="normal">2</mn><mo mathvariant="normal">*</mo><mi mathvariant="normal">a</mi><mo mathvariant="normal">*</mo><msup><mi mathvariant="normal">t</mi><mn mathvariant="normal">2</mn></msup><mo mathvariant="normal">+</mo><msub><mi>HDOT</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><mi mathvariant="normal">t</mi></math><img file="EP0925518B1_D0017.tif" /></maths>

0097Assuming a constant acceleration during the pull-up maneuver, the time t until vertical speed reaches zero (point 348) is given by equation (18). <maths id="math0018" num="(18)"><math display="block"><mi mathvariant="normal">t</mi><mo mathvariant="normal">=</mo><mo mathvariant="normal">-</mo><msub><mi>HDOT</mi><mn mathvariant="normal">0</mn></msub><mo mathvariant="normal">/</mo><mi mathvariant="normal">a</mi></math><img file="EP0925518B1_D0018.tif" /></maths>

0098Substituting equation (18) into equation (17) yields equation (19). <maths id="math0019" num="(19)"><math display="block"><mi>ALPU</mi><mo mathvariant="normal">=</mo><mo mathvariant="normal">-</mo><mrow><mo mathvariant="normal">(</mo><msub><mi>HDOT</mi><mn mathvariant="normal">0</mn></msub><mo>⁢</mo><msup><mrow><mo mathvariant="normal">)</mo></mrow><mn mathvariant="normal">2</mn></msup></mrow><mo mathvariant="normal">/</mo><mfenced><mn mathvariant="normal">2</mn><mo mathvariant="normal">*</mo><mi mathvariant="normal">a</mi></mfenced></math><img file="EP0925518B1_D0019.tif" /></maths> Equation (19) thus represents the altitude loss during the pull-up maneuver.

0099An exemplary block diagram for generating the signals ALTR and ALPU is illustrated in <figref idref="f0014">FIG. 21</figref>. In particular, a signal representative of the vertical velocity of the aircraft HDOT, available, for example, from a barometric altimeter rate circuit (not shown), is applied to a filter 350 in order to reduce nuisance warnings due to turbulence. The filter 350 may be selected with a transfer function of 1/(TAHDOT * S + 1); where TAHDOT is equal to one second. The output of the filter 350 is a signal HDOTf, which represents the filtered instantaneous vertical speed; positive during climbing and negative during descent.

0100In order to generate the altitude loss due to reaction time signal ALTR, the signal HDOTf is applied to a multiplier 352. Assuming a pilot reaction time Tr, for example, of 5 seconds, a constant 354 equal to 5 seconds is applied to another input of the multiplier 352. The output of the multiplier 352 represents the signal ALTR, which is positive when HDOTf is negative and set to zero if the signal HDOTf is positive, which represents a climbing condition. More particularly, the signal HDOTf is applied to a comparator 356 and compared with a reference value, for example, zero. If the comparator 356 indicates that the signal HDOTf is negative, the signal HDOTf is applied to the multiplier 352. During climbing conditions, the signal HDOTf will be positive. During such conditions, the comparator 356 will apply a zero to the multiplier 352.

0101The altitude loss due to the pull-up maneuver signal ALPU is developed by a square device 358, a divider 360 and a multiplier 362. The filtered instantaneous vertical speed signal HDOTf is applied to the square device 358. Assuming a constant acceleration during the pull-up maneuver, for example, equal to 8 feet/ sec² (0.25 G), a constant is applied to the multiplier 362 to generate the signal 2a. This signal, 2a, is applied to the divider 360 along with the output of the square device 350. The output of the divider 360 is a signal (HDOT f)²/2a, which represents the altitude loss during a pull-up maneuver signal ALPU.

0102These signals ALRT and ALPU are used to modulate the distance below the aircraft where terrain advisory and terrain warning boundaries begin during a LOOK-UP mode of operation. More particularly, during such a mode of operation, DH terrain floor segment of the terrain advisory boundary DHYEL2 is modulated by the signals ALRT and ALPU while the DH terrain floor segment of the terrain warning boundary DHRED2 is modulated by the signal ALPU as indicated in equations (20) and (21), respectively. <maths id="math0020" num="(20)"><math display="block"><mi>DHYEL</mi><mo>⁢</mo><mn>2</mn><mo>=</mo><mn>3</mn><mo>/</mo><mn>4</mn><mo>*</mo><mi>DH</mi><mo>+</mo><mi>ALRT</mi><mo>+</mo><mi>ALPU</mi></math><img file="EP0925518B1_D0020.tif" /></maths><maths id="math0021" num="(21)"><math display="block"><mi>DHRED</mi><mo>⁢</mo><mn mathvariant="normal">2</mn><mo mathvariant="normal">+</mo><mn mathvariant="normal">1</mn><mo mathvariant="normal">/</mo><mn mathvariant="normal">2</mn><mo mathvariant="normal">*</mo><mi mathvariant="normal">D</mi><mo mathvariant="normal">+</mo><mi>ALPU</mi><mo mathvariant="normal">,</mo></math><img file="EP0925518B1_D0021.tif" /></maths> where DH represents the terrain floor as discussed above.

0103Thus, referring to <figref idref="f0012">FIG. 19</figref>, the LOOK-UP terrain advisory warning begins at a point 364 below the aircraft; equal to DHYEL2. If the flight path angle g is less than a configurable datum THETA2, a terrain advisory boundary 366 extends from the point 364 to the advisory LAD at an angle equal to THETA2. Should the flight path angle g be greater than THETA2, the lower advisory boundary, identified with the reference numeral 368, will extend from the point 364 at an angle equal to the flight path angle g.

0104Similarly, the LOOK-UP terrain warning boundary begins at point 370 below the aircraft; equal to DHRED2. If the flight path angle g is less than the angle THETA2, a warning boundary 372 extends from the point 370 at angle THETA2 to the warning LAD. Should the flight path angle g be greater than THETA2, a warning boundary 374 will extend at an angle equal to-the flight path angle g between the point 370 and the warning LAD.

CUT-OFF ALTITUDE

0105The cut-off altitude is an altitude relative to the nearest runway elevation; set at, for example, 500 feet. Altitudes below the cut-of altitude are not displayed and are ignored by the terrain advisory and terrain warning computations.

0106The use of a cut-off altitude declutters the display discussed below around airports, especially during a final approach when the aircraft approaches the ground. A second advantage is that nuisance warnings on a final approach, due to altitude errors, terrain data base resolution and accuracy errors are minimized.

0107However, the use of a cut-off altitude during certain conditions, such as an approach to an airport on a bluff (i.e., Paine Field) at a relatively low altitude or even at an altitude below the airport altitude, may compromise system performance. More particularly, during such conditions, the use of a cut-off altitude may prevent a terrain warning from being generated because the threatening terrain may be below the cut-off altitude. In order to account of such situations, the system selects the lower of two cut-off altitudes; a nearest runway cut-off altitude (NRCA) and a cut-off altitude relative to aircraft (CARA). The NRCA is a fixed cut-off altitude, relative to the nearest runway. The CARA is an altitude below the instantaneous aircraft altitude (ACA) by an amount essentially equivalent to DH terrain floor boundary for an advisory condition, which, as discussed above, is a function of the aircraft distance to the nearest runway.

0108Equations (22) and (23) below set forth the NRCA and CARA. As mentioned above, the absolute cut-off altitude (ACOA) is the lower of the NRCA and CARA as set forth in equation (24). <maths id="math0022" num="(22)"><math display="block"><mi>NRCA</mi><mo>=</mo><mi>COH</mi><mo>+</mo><mi>RE</mi><mo>,</mo></math><img file="EP0925518B1_D0022.tif" /></maths> where COH relates to the cut-off height and is a fixed configurable value, initially set between 400 feet and 500 feet; and RE equals the runway elevation. <maths id="math0023" num="(23)"><math display="block"><mi>CARA</mi><mo>=</mo><mi>ACA</mi><mo>-</mo><mi>DH</mi><mo>-</mo><mi>DHO</mi><mo>,</mo></math><img file="EP0925518B1_D0023.tif" /></maths> where ACA is the instantaneous aircraft altitude; DH is the terrain floor which is proportional to the distance to the runway as discussed above; and DHO is a configurable bias, set to, for example, 50 feet. <maths id="math0024" num="(24)"><math display="block"><mi>ACOA</mi><mo>=</mo><mi>lower of CARA</mi><mo>,</mo><mi>NRCA</mi><mo>,</mo></math><img file="EP0925518B1_D0024.tif" /></maths> except where DH DH1, in which case ACOA is always equal to NRCA.

0109DH1 is a point at which the ACOA is forced to be equal to NRCA independent of the aircraft altitude. The point DH1 is related to COH, DH and DHO such that on a nominal three (3) degree approach slope, CARA is equal to NRCA when the aircraft is at a distance equal to a distance DH1 from the airport, as illustrated in TABLE 2 below. <tables id="tabl0002" num="0002"><table frame="all"><title><u style="single">TABLE 2</u></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="13mm" /><colspec colnum="2" colname="col2" colwidth="13mm" /><colspec colnum="3" colname="col3" colwidth="31mm" /><thead><row rowsep="0"><entry align="center" valign="top"><b>COH</b></entry><entry align="center" valign="top"><b>DH1</b></entry><entry align="center" valign="top"><b>DISTANCE TO</b></entry></row><row><entry align="center" valign="top"><b>(feet)</b></entry><entry align="center" valign="top"><b>(feet)</b></entry><entry align="center" valign="top"><b>RUNWAY (n mile)</b></entry></row></thead><tbody><row><entry align="center">300</entry><entry align="center">50</entry><entry align="center">1</entry></row><row><entry align="center">400</entry><entry align="center">100</entry><entry align="center">1.5</entry></row><row><entry align="center">500</entry><entry align="center">150</entry><entry align="center">2</entry></row></tbody></tgroup></table></tables>

0110The point DH1 forces the cut-off altitude (COH) above the runway whenever the aircraft is close to the runway to ensure robustness against nuisance warnings caused by altitude errors and terrain data base resolution margins to disable the terrain advisory and warning indications when the aircraft is within the aircraft perimeter. There are trade-offs between nuisance warnings and legitimate warnings. In particular, the lower COH, the closer the advisory and warning indications are given, making the system more nuisance prone. As indicated above, for a COH of 400, terrain advisory and warning indications are effectively disabled when the aircraft is closer than 1.5 n miles from the airport runway.

0111<figref idref="f0015">FIG. 22</figref> illustrates the operation of the alternative cut-off altitude boundaries. In particular, <figref idref="f0015">FIG. 22</figref> illustrates a condition when the COH is set to 300 feet with DH1 equal to 50 feet. The cut-off altitude substantially for area from the runway, for example, greater than 4 nautical miles, is 300 feet, as indicated by the segment 378 when the glide slope angle is less than a predetermined angle, for example, 3°. As the aircraft gets closer to the runway, the COH is lowered, as illustrated by the segment 388, until the aircraft is within one (1) nautical mile of the runway, at which point the COH is forced to be 300 feet, which effectively disables any terrain advisory and warning indications within the aircraft is closer than one (1) nautical mile from the runway, as represented by the segment 382.

0112During a condition when the aircraft is on a, for example, 3° glide slope angle, the ACOA is forced to be the NRCA. As shown, the NRCA is illustrated by the segment 382.

DISPLAY SYSTEM

0113A display system, generally identified with the reference numeral 400, is illustrated in <figref idref="f0016 f0017 f0018 f0019 f0020 f0021 f0022 f0023 f0024 f0025 f0026 f0027 f0028 f0029 f0030 f0031">FIGS. 23-42</figref>. The display system 400 is used to provide a visual indication of the terrain advisory and terrain warning indications discussed above as a function of the current position of the aircraft. Background terrain information is also provided which provides an indication of significant terrain relative to the current position of the aircraft.

0114In order to declutter the display, background information may be displayed in terms of predetermined dot patterns whose density varies as a function of the elevation of the terrain relative to the altitude of the aircraft. Terrain advisory and terrain warning indications may be displayed in solid colors, such as yellow and red, relative to the groundtrack of the aircraft.

0115An important aspect of the invention is that the terrain background information as well as the terrain advisory and warning indications may be displayed on a navigational or weather display, normally existing within an aircraft, which reduces the cost of the system and obviates the need for extensive modifications to existing displays, such as a navigational and weather radar type display. In particular, the terrain data is converted to a weather radar format in accordance with the standard ARINC 708/453 for digital buses on civil aircraft. The terrain data, masked as weather data, can then be easily displayed on an existing navigational or a dedicated weather radar display that conforms with the ARINC 708/453 serial interface standard.

0116Referring to <figref idref="f0016">FIG. 23</figref>, the display system 400 includes a weather radar display or navigational display 36 connected to a control bus 404 conforming to an ARINC 429 standard and a terrain data bus switch 406. The terrain data bus switch 406 enables selection between terrain data and weather data for display. More particularly, the data bus switch 406 includes a common pole 408 connected to the display 36 by way of a serial display bus 410. One contact 412 on the data bus switch 406 is connected to a serial weather data bus 414 and, in turn, to a weather radar R/T unit 416, which transmits weather radar data over the weather data bus 414 in conformance with the ARINC 708/453 standard.

0117An antenna 418 provides weather radar data to the weather radar R/T unit 416. The weather radar R/T unit 416 is also connected to the control bus 404, which is used for range and mode data.

0118In accordance with an important aspect of the invention, the terrain advisory and terrain warning indications, as well as the background information discussed above, are converted to "RHO/THETA" format. The converted terrain data is then applied to a terrain data serial bus 420, which conforms to the ARINC 453 standard, and connected to a contact 422 on the data bus switch 406 to enable selective display of either weather radar or terrain data on the display 36.

0119The system for converting the terrain advisory and warning indications as well as background terrain data to an ARINC 708/453 standard serial format is indicated as a functional block 424 in <figref idref="f0016">FIG. 23</figref> and is described in more detail below. By converting such data into an ARINC 708/453 standard serial format, the data can be displayed on an existing display 36, which conforms to the ARINC 453 standard without the need for relatively complex and expensive modifications to the display's symbol generator.

0120Referring to <figref idref="f0017">FIG. 24</figref>, each word 425, in accordance with the ARINC 453/708 standard, is 1600 bits long and represents one spoke or radial 426 of the display 36. Normally, 512 spokes or radials 426 are transmitted for each complete sweep of the display 36. At a transmission rate of, for example, 1 megahertz (MHz), each word 425 takes about 1.6 milliseconds (msec) to transmit. For 512 spokes, one complete sweep will thus take about 4 seconds.

0121The 1600-bit words 425 include a header 428 and 512 range bins 430. The header 428 includes the control data on the control bus 404 and also includes a 12-bit code representing the antenna scan angle; the angle of the spoke relative to the aircraft heading equivalent to an indicator UP direction. The range bins 430 cover the distance from the antenna 418 to the maximum selected. Each of the 512 range bins includes three intensity bits, encoded into colors as indicated in TABLE 3 below. <tables id="tabl0003" num="0003"><table frame="all"><title><u style="single">TABLE 3</u></title><tgroup cols="2"><colspec colnum="1" colname="col1" colwidth="35mm" /><colspec colnum="2" colname="col2" colwidth="21mm" /><thead><row><entry align="center" valign="top"><b>3-BIT RANGE CODE</b></entry><entry align="center" valign="top"><b>COLOR</b></entry></row></thead><tbody><row><entry align="center">000</entry><entry align="center">BLACK</entry></row><row><entry align="center">100</entry><entry align="center">GREEN</entry></row><row><entry align="center">010</entry><entry align="center">YELLOW</entry></row><row><entry align="center">110</entry><entry align="center">RED</entry></row><row><entry align="center">001</entry><entry align="center">-----</entry></row><row><entry align="center">101</entry><entry align="center">CYAN</entry></row><row><entry align="center">011</entry><entry align="center">MAGENTA</entry></row><row><entry align="center">111</entry><entry align="center">-----</entry></row></tbody></tgroup></table></tables>

0122As will be discussed in more detail below, the display system 400 is capable of displaying background terrain information as well as terrain threat indications as a function of the current position of the aircraft. The threat detection algorithm is run along the groundtrack of the aircraft. Thus, terrain threat indications are typically displayed along the groundtrack while background terrain is displayed relative to the heading of the aircraft.

0123Referring to <figref idref="f0018">FIG. 25</figref>, the terrain background information is shown on the display 36. As will be discussed in more detail below, the elevation of the highest terrain relative to the altitude of the aircraft is shown as a series of dot patterns whose density varies as a function of the distance between the aircraft and the terrain. For example, a relatively dense dot pattern 432 may be used to indicate terrain that is, for example, 500 feet or less below the aircraft. A medium dense dot pattern 434 may be used to represent terrain that is 1000 feet or less below the aircraft, while a lightly dotted pattern 436 may be used to indicate terrain 2000 feet or less below the aircraft. In order to declutter the display 402, terrain more than, for example, 2000 feet below the aircraft, is not shown. The dots may be displayed in one of the colors indicated in TABLE 3 above, for example, yellow (amber). In addition, in the vicinity of an airport, a runway graph 438, for example, in green, may also be provided. The apparent display resolution may be increased by using more than one color along with the variable density dot patterns, for example, as shown in TABLE 4 below. <tables id="tabl0004" num="0004"><table frame="all"><title><u style="single">TABLE 4</u></title><tgroup cols="4"><colspec colnum="1" colname="col1" colwidth="41mm" /><colspec colnum="2" colname="col2" colwidth="18mm" /><colspec colnum="3" colname="col3" colwidth="17mm" /><colspec colnum="4" colname="col4" colwidth="19mm" /><tbody><row><entry morerows="2" align="center">DOT DENSITY PATTERN</entry><entry align="center">HIGH</entry><entry align="center">3000</entry><entry align="center">0</entry></row><row><entry align="center">MEDIUM</entry><entry align="center">4000</entry><entry align="center">1000</entry></row><row><entry align="center">LOW</entry><entry align="center">5000</entry><entry align="center">2000</entry></row><row><entry align="center">COLOR</entry><entry align="center" /><entry align="center">GREEN</entry><entry align="center">YELLOW</entry></row></tbody></tgroup></table></tables>

0124The display of a terrain threat indication is illustrated in <figref idref="f0019">FIG. 26</figref>. As shown, the terrain threat indication may be displayed contemporaneously with the background terrain information. Terrain advisory and terrain warning indications are displayed in solid shapes 440 and 442, respectively, for example, "squares"; the displayed terrain map cells which represent a threat painted solid yellow or red. More particularly, at an aspect ratio for the display of, for example, 3 x 4 (vertical to horizontal), the terrain cells will appear "square", particularly at relatively low-range settings. Colors are used to distinguish between terrain advisory and terrain warning indications. For example, red may be used to represent a terrain warning indication 442 while yellow or amber is used to represent a terrain advisory indication 440. By using colored shapes for terrain threat indications and dot patterns of variable density for the terrain background information, clutter of the display 402 is minimized.

0125The terrain data from the terrain data base 24, for example, as illustrated in <figref idref="f0005">FIG. 4B</figref>, is used to compute the terrain threat indications as well as background information for the display 36. As discussed above, the terrain data base 24 may be stored in a flash ROM 444 (<figref idref="f0020">FIG. 27</figref>) for convenience of updates. In order to facilitate update of the display 36 as a function of the position of the aircraft, terrain data from the flash ROM 444 is transferred to a random access memory (RAM) A, identified with the reference numeral 446. As will be discussed in more detail below, the RAM A is mapped with terrain data from the flash ROM 444, centered about the current position of the aircraft, as illustrated in <figref idref="f0020">FIG. 27</figref>.

0126In order to provide a relatively large display range, for example, 160 nautical miles, while minimizing processing time as well as the size of the RAM A, the RAM A is configured as a "wedding cake" as generally illustrated in <figref idref="f0021">FIG. 28</figref> and configured such that the highest terrain resolution is near the position of the aircraft and the terrain farthest from the aircraft position has the lowest resolution. For example, the RAM A may be provided with four (4) layers. As shown, the instantaneous aircraft position is shown generally centered in layer 1, having the highest resolution. The resolution of layer 2 lies between the resolution of layer 3 and layer 1. Layer 4 (not shown) is provided with a relatively coarse resolution for use in a cruise mode.

0127Such a variable resolution display, configured as a layered wedding cake, minimizes the memory size of RAM A as well as the processing time. In order to improve the performance of the system, the flash ROM 444 may also be configured similar to a "wedding cake", as generally illustrated in <figref idref="f0020">FIG. 27</figref>.

0128As discussed above, the terrain data base 24 may be provided with variable resolution with a relatively fine resolution of, for example, 30 x 20 arc sec (i.e. 1/2 x 1/2 arc min). The relatively coarse resolution of layers 2 and 3 can be formed by combining finer resolution maps. With such a resolution, the configuration of RAM A may be selected in accordance with TABLE 5. <tables id="tabl0005" num="0005"><table frame="all"><title><u style="single">TABLE 5</u></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="16mm" /><colspec colnum="2" colname="col2" colwidth="26mm" /><colspec colnum="3" colname="col3" colwidth="32mm" /><thead><row rowsep="0"><entry align="center" valign="top"><b>LAYER</b></entry><entry align="center" valign="top"><b>RESOLUTION</b></entry><entry align="center" valign="top"><b>LAYER SIZE</b></entry></row><row><entry align="center" valign="top" /><entry align="center" valign="top"><b>(MINUTES)</b></entry><entry align="center" valign="top"><b>(MINUTES)</b></entry></row></thead><tbody><row><entry align="center">1</entry><entry align="center">1/2 x 1/2</entry><entry align="center">32 x 32</entry></row><row><entry align="center">2</entry><entry align="center">1 x 1</entry><entry align="center">64 x 64</entry></row><row><entry align="center">3</entry><entry align="center">2 x 2</entry><entry align="center">128 x 128</entry></row><row><entry align="center">4</entry><entry align="center">5 x 5</entry><entry align="center">REMAINING AREA</entry></row></tbody></tgroup></table></tables>

0129Each layer is generally centered relative to the aircraft position as shown in <figref idref="f0021">FIG. 28</figref> at the time of the update of RAM A. Each layer of the "wedding cake" has its own update boundary. For example, layer 1 of the RAM A includes an update boundary 448. As the aircraft passes through the update boundary 448, the RAM A is updated with new terrain data from the flash ROM 444 (<figref idref="f0020">FIG. 27</figref>) as will be discussed in more detail below.

0130Because of the refresh rates required for the display 36, the terrain data may be stored in files alternative to the format illustrated in <figref idref="f0003">FIG. 4A</figref>. For example, the file format of the terrain data may include an index file or header 55, stored in one or more 128 K byte storage blocks and a plurality of data files 51, as generally shown in <figref idref="f0022">FIG. 29</figref>.

0131As shown, each index file 55 may include pointers to specific data files 51, the length of the map file, as well as the position of one or more of the corner boundaries of the file 51, for example, northwest and southeast corners of each map file 51, which facilitates updates of the display 36 as a function of the current aircraft position. The index file 55 may contain information regarding the length of the map file 51, a start of file or altitude offset, the resolution of the map file 51, and whether and how the file 51 is compressed to enable use of various compression algorithms.

0132The data files 51 may include a file header 454 as well as data blocks 456 which, as shown, can be grouped according to resolution size. The file header 454 indicates various information about the file. The format of the header 454 may be structured as follows: <pre listing-type="program-listing"> <MAJOR VERSION BYTE> < MINOR VERSION BYTE> < FILE STATUS BYTE> < ATTRIB BYTE > <FILE NAME (8 CHARACTERS) > < EXTENSION (4 CHARACTERS) > <FILE LENGTH LONGWORD> <TIME STAMP LONGWORD > <CRC LONGWORD> <SPME LONGWORD></pre>

0133The major and minor version bytes relate to the revision status of the data in the data files 51. The minor version is set as the inverse of the byte value of the major version such that a default erase will indicate a zero version. The file status byte 457 provides various information regarding the status of the file as discussed below, while the attribute byte is reserved. The file name, extension file length, and spare data are self-explanatory. The time stamp relates to the operational time of the system (i.e., length of time in use) while the CRC data relates to the data files 51.

0134The status byte may be used to provide the following information regarding the data blocks 456, for example, as set forth in TABLE 6 below. <tables id="tabl0006" num="0006"><table frame="all"><title><u style="single">TABLE 6</u></title><tgroup cols="2"><colspec colnum="1" colname="col1" colwidth="34mm" /><colspec colnum="2" colname="col2" colwidth="33mm" /><thead><row><entry align="center" valign="top"><b>FLAG BYTE VALUE</b></entry><entry align="center" valign="top"><b>INFORMATION</b></entry></row></thead><tbody><row><entry align="center">FF</entry><entry align="center">EMPTY</entry></row><row><entry align="center">FE</entry><entry align="center">DOWNLOADING</entry></row><row><entry align="center">FO</entry><entry align="center">VALID</entry></row><row><entry align="center">EO</entry><entry align="center">READY TO ERASE</entry></row><row><entry align="center">·</entry><entry align="center">ESERVED</entry></row><row><entry align="center">·</entry><entry align="center">RESERVED</entry></row><row><entry align="center">00</entry><entry align="center">RESERVED</entry></row></tbody></tgroup></table></tables>

0135The EMPTY status (FF) is set by erasing the data blocks 456, while the DOWNLOADING status (FE) is set as part of the image. The VALID status (FO) is set when the file is loaded. Lastly, the READY TO ERASE (EO) is set by the file system for obsolete blocks. The file structure may also be as described in Appendix 1.

0136As indicated below, various corrections of the data files 51 may be required in areas near runways.

0137<figref idref="f0022">FIG. 30</figref> illustrates a portion of the terrain data base 24 around a particular airport, Boeing Field, with a resolution of, for example, 1/2 x 1/2 minute subcells 455, shown with a runway 456 drawn in. The elevations, as discussed above, are selected to be the highest elevations within the cells 455, rounded up to the next highest resolution increment, for example 100 feet. While such a system provides useful and conservative terrain data for use with the terrain data base 24, such a system in the vicinity of an airport can cause the elevations of the runway 456 to be too high or too low. For example, the elevation of the runway 456 at Boeing Field is known to be 17/15 feet; rounded off to zero (0) feet. <figref idref="f0022">FIG. 30</figref> shows the elevation of runway 456 to be 200/300 feet, for example, by the averaging algorithm discussed above. Such errors in the runway elevation can degrade system performance, and, in particular, cause improper performance of the terrain advisory and terrain warning cut-off boundaries below 500 feet, for example, as illustrated in <figref idref="f0015">FIG. 22</figref>.

0138In order to solve this problem, the actual elevation of the runway 456, rounded to the nearest resolution, is used for cells 455 around a perimeter of the runway 456, as generally shown in <figref idref="f0023">FIG. 31</figref>. In particular, a first perimeter, illustrated by the dashed box 458, is selected to be about ±1 cell (i.e., 1/2 minute) from the centerline of the runway 456. The elevations in all of the cells 455 intersected by the dashed box 458 are set to the actual runway elevation and rounded to the nearest resolution. For Boeing Field, the elevation for such cells is 0 feet.

0139Consideration is also given for a narrow approach and a wide approach. The dashed box 460 represents a wide approach while the dashed box 462 represents a narrow approach. The cells 464 and 466, intersected by the wide approach box 460, are set on the actual runway elevation, rounded to 0 feet. Inside the approach perimeter, as indicated by the arc 467, the elevation for the cells is selected to be either the runway elevation or the DH terrain floor (as discussed above). For example, if the corner of the cell 468, closest to the runway 456, is 1.386 nM from the runway 456, the DH terrain floor will be 86 feet: rounded to nearest 100 feet = 100 feet. Thus, the elevation of cell 468 is selected as 100 feet.

0140<figref idref="f0023">FIG. 32</figref> illustrates the differences between the elevations in selected cells 468, 470 and 472 before and after correction. Before correction, the elevation of cells 468, 470 and 472 was 300 feet. After correction, the elevation of cells 470 and 472 are corrected to the runway elevation and rounded to the nearest resolution (i.e., 0 feet). The cell 468 is corrected to the DH terrain floor elevation of 100 feet.

0141<figref idref="f0023">FIG. 33</figref> illustrates another situation where the elevation of cells, due to the rounding-off method discussed above, can result in a situation where the cell elevations adjacent a runway, represented as a staircase 474, are below the runway elevation. The correction method, discussed above, can also be used to correct the various cells around the runway 456 up to the runway elevation in a manner as generally discussed above.

0142<figref idref="f0024">FIG. 34</figref> represents a simplified block diagram for implementation for the display system 400 in accordance with the present invention. The display system 400 may include a microprocessor, for example, an Intel type 80486, 25 MHz type microprocessor ("486") and an Analog Devices type digital signal processor (DSP). The DSP is primarily used for calculating the RHO/THETA conversions to off load the 486.

0143The display system 400 may have a normal mode and a cruise mode. In a normal mode, the display system 400 displays terrain out to, for example, 100 nautical miles, while in cruise mode out to about 320 nautical miles. In the cruise mode, it is contemplated that the background terrain display levels can be lowered by the pilot such that terrain, for example 10,000 feet or 20,000 feet below the aircraft, can be displayed. Alternatively, the display level could automatically be adjusted to the highest terrain within the selected display range; the value of this level set forth in amber, for example.

0144The 486 runs a "file picker" routine which searches the terrain data base 24 for terrain data covering the area relative to the current position of the aircraft. For example, in a normal mode, the file picker routine searches the terrain data base 24 for terrain data covering an area approximately 256 x 256 minutes relative to the current position of the aircraft. As indicated above, the terrain data base 24 includes index files 55 (<figref idref="f0005">FIGS. 4B</figref> and <figref idref="f0022">29</figref>) which may include the northwest and southeast corners of each map file 51 in order to facilitate the search. A list of files covering the area is assembled and maintained.

0145As indicated above, the index file 55 also includes information as to whether and how the data files 51 are compressed. The files may be compressed/decompressed by a modified Hoffman code, such as PKZIP as manufactured by PKWARE, INC.

0146If the data files 51 are compressed, the files are decompressed in the background because of the relatively long decompression time and stored in a decompressed file RAM (DFR) 476, for example 256 K bytes. The decompressed files are then transferred to RAM A, identified with the reference number 446. On power-up, one or two files around the aircraft are decompressed and displayed to allow immediate display of a limited range, relative to the position of the aircraft. As more terrain data files are decompressed, the additional range is displayed.

0147The RAM A may be configured in a plurality of cells 479 (<figref idref="f0024">FIG. 34</figref>), representative of a terrain area 1/2 x 1/2 minutes, with true north appearing at the top of RAM A as shown. Each cell 479 is two bytes wide and includes data regarding the highest altitude in the cell as well as the identity of each cell 479. As such, the RAM A is addressable by X and Y coordinates and thus represents a map with an altitude in each cell 479. The identity of the cell may be used to distinguish between terrain data, for example, from the terrain data base 24 and obstacle data. As discussed above, in some applications, an obstacle data base may be utilized.

0148As indicated above, the RAM A may be configured as a wedding cake with various layers, as generally shown in <figref idref="f0021">FIG. 28</figref> with the top layer (i.e., layer 1) having the highest resolution (i.e., 1/2 x 1/2 minute cells). Terrain data is loaded into the RAM A such that the aircraft is roughly centered in the top layer, where the terrain advisory and warning detection algorithms are run. Once the aircraft moves outside the update boundary 448 (<figref idref="f0021">FIG. 28</figref>), the top layer of the RAM A is reloaded with terrain data with the aircraft centered. The lower layers may also be provided with update boundaries, which enable the lower layers to be updated independently and less frequently than the top layer.

0149As discussed above and illustrated in <figref idref="f0020">FIG. 27</figref>, the flash ROM 444 for the terrain data base 24 may also be configured as a wedding cake as illustrated in <figref idref="f0021">FIG. 28</figref>. In such an implementation, the highest resolution is provided near the airports while the coarser resolution is provided away from the airports as generally shown in <figref idref="f0020">FIG. 27</figref>. For example, as illustrated in <figref idref="f0025">FIG. 34A</figref>, the terrain data base 24 may be configured with 1° x 1°, 2° x 2° or 10' x 10° map files. Some of the map files may overlap, similar to a "wedding cake". As shown, the terrain maps closest to airports are provided with the highest resolution, for example 15 x 15 arc seconds. Maps further away from the airport may be provided with resolutions of, for example, 1 x 1 minute or 2 x 2 minutes while a coarse resolution (i.e., 2° x 2' or 10° x 10°) may be provided even further from the airport.

0150An upload routine may be used to upsample or downsample terrain or obstacle data into the RAM A from the DFR RAM 476 as shown in <figref idref="f0024">FIG. 34</figref>. Such a routine is well within the ordinary skill in the art. The top three layers (i.e., layers 1, 2 and 3) of the RAM A (<figref idref="f0021">FIG. 28</figref>) are updated from the map file source from the flash ROM 444 for the terrain data base 24 having the highest resolution, except for the 5 x 5 minute layer, which is updated only from a 5 x 5 minute file source to reduce processing time.

0151Alternatively, as shown in <figref idref="f0025">FIG. 34A</figref>, an additional RAM A*, identified with the reference numeral 477. may be provided and structured similar to the RAM A 446, forming a "ping-pong" memory arrangement. In this embodiment, one terrain data file 51 is decompressed one file at a time and loaded into the DFR RAM 476 and then upsampled to either the RAM A or RAM A*. For example, while the RAM A is being loaded with additional decompressed files, the terrain threat and display algorithms are executed from the RAM A*. Once the RAM A is updated, the RAM A is then used for the computations while the other RAM A* is updated with new data from the DFR RAM 476 to account for changes in the aircraft position.

0152The 5 x 5 minute resolution layer is primarily intended for use in a cruise mode to provide a look-ahead range of about 320 nM. As discussed above and illustrated in <figref idref="f0022">FIG. 29</figref>, each file includes an index file 55 which contains data relative to the resolution of the terrain data in a particular map file 51. Such information can be used to vary the resolution of the map files as a function of longitude and latitude position. For example, the latitude resolutions of the map files can be decreased at higher latitudes in a northward direction. Moreover, longitude resolution in the RAM A can be varied and provided with, for example, a one degree latitude hysteresis to prevent flipping back and forth as an aircraft cruises northwardly. Once the terrain data from the terrain data base 24 is loaded into the RAM A, the terrain advisory and warning algorithms ("detection algorithms") are run once per second as discussed below.

0153The terrain threat algorithms discussed above are based on a single vector along the groundtrack. The single vector provides acceptable results because the relatively crude resolution of the terrain data base 24 is based upon the highest elevation per cell and the natural noisiness of the groundtrack angle and position information. However, the single vector approach does not take into account various errors, due to, for example, the terrain data base 24, GPS longitude and latitude errors as well as the published track angle accuracy. As such, the single vector along the groundtrack may be substituted with an array of vectors, as generally shown in <figref idref="f0026">FIG. 35</figref>, which takes such errors into account. In particular, the groundtrack vector is identified with the reference numeral 482 originating at the instantaneous aircraft position X₀, Y₀. A pair of vectors 484 and 486 are disposed at a configurable distance D_OFF, for example 0.2 nM along a segment 488, perpendicular to the vector 482 at points X₁, Y₁ and X_r, Y_r. A pair of outside vectors 490 and 492, originating at points X₁, Y₁ and X_r, Y_r, respectively, are directed at a configurable angle D_ALPHA, relative to the vectors 484 and 486. The configurable angle D_ALPHA may be selected as 3 degrees, equivalent to the groundtrack accuracy of the GPS.

0154The coordinates X₁, Y₁ and X_r, Y_r for the vectors 490 and 492, respectively, are determined. More particularly, the incremental X and Y vectors, DX_OFF and DY_OFF, are determined as set forth in equations (25) and (26) as set forth below. <maths id="math0025" num="(25)"><math display="block"><msub><mi>DX</mi><mi>OFF</mi></msub><mo mathvariant="normal">=</mo><mfenced open="{" close="}"><msub><mi mathvariant="normal">D</mi><mi>OFF</mi></msub><mo mathvariant="normal">/</mo><mi>COS</mi><mfenced><mi>LAT</mi></mfenced></mfenced><mo mathvariant="normal">*</mo><mi>COS</mi><mfenced><mi>GROUNDTRACK</mi></mfenced><mo mathvariant="normal">,</mo></math><img file="EP0925518B1_D0025.tif" /></maths> where D_OFF is a configurable distance between X₀, Y₀ and X₁, Y₁, for example 0.2 nm and the factor 1/COS(LAT) is used to convert DX_OFF to nautical miles when DX_OFF, DY_OFF is in minutes and D_OFF is in nautical miles. <maths id="math0026" num="(26)"><math display="block"><msub><mi>DY</mi><mi>OFF</mi></msub><mo mathvariant="normal">=</mo><msub><mi mathvariant="normal">D</mi><mi>OFF</mi></msub><mo mathvariant="normal">*</mo><mi>SIN</mi><mfenced><mi>GROUNDTRACK</mi></mfenced></math><img file="EP0925518B1_D0026.tif" /></maths> The coordinates X_r, Y_r can be determined in a similar manner.

0155The threat detection algorithms are initiated by calculating the range increments DELTA X_NS and DELTA Y_NS along a north-south coordinate system along the groundtrack. Referring to <figref idref="f0023">FIG. 36</figref>, each spoke 426 has 512 size range bins as discussed above, or 256 double range bins. In order to convert from an X-Y coordinate system to the north-south coordinate system, the range increments DELTA X_NS, DELTA Y_NS are determined in accordance with equations (27) and (28). (The term 1/COS (LATITUDE) is used to convert from arc minutes to nautical miles.) <maths id="math0027" num="(27)"><math display="block"><mi>DELTA</mi><mspace width="1em" /><msub><mi mathvariant="normal">X</mi><mi>NS</mi></msub><mo mathvariant="normal">=</mo><mi>DETECTION STEP SIZE</mi><mo mathvariant="normal">*</mo><mi>SIN</mi><mfenced><mi>GROUNDTRACK</mi></mfenced><mo mathvariant="normal">*</mo><mn mathvariant="normal">1</mn><mo mathvariant="normal">/</mo><mi>COS</mi><mfenced><mi>LATITUDE</mi></mfenced></math><img file="EP0925518B1_D0027.tif" /></maths><maths id="math0028" num="(28)"><math display="block"><mi>DELTA</mi><mspace width="1em" /><msub><mi mathvariant="normal">Y</mi><mi>NS</mi></msub><mo mathvariant="normal">=</mo><mi>DETECTION STEP SIZE</mi><mo mathvariant="normal">*</mo><mi>COS</mi><mfenced><mi>GROUNDTRACK</mi></mfenced></math><img file="EP0925518B1_D0028.tif" /></maths> After the range increments DELTA X_NS and DELTA Y_NS are determined, the threat detection algorithm computation is initiated at the instantaneous position X₀, Y₀ of the aircraft (<figref idref="f0027">FIG. 37</figref>). The position is then incremented by the range increments DELTA X_NS and DELTA Y_NS and the range is incremented by a configurable step size, for example 1/8 unit, as discussed above. After the new coordinates are determined, the current terrain altitude at the new coordinate is determined. The terrain detection algorithms, as discussed above, are then computed for the current advisory LAD, as indicated by the arc 494, defining a yellow threat detection vector 496. In order to provide additional warning time for a terrain warning, the terrain warning indication can be computed out to a distance equal to 1.5 * LAD, as indicated by the arc 498 to define a red threat detection vector 500. The above steps are repeated for the look-ahead vector arrays 490 and 492 (<figref idref="f0026">FIG. 35</figref>). The entire computation is done at a rate of once per second.

0156Whenever the threat detection algorithm detects terrain along the look-ahead vector arrays 490 and 492 (<figref idref="f0026">FIG. 35</figref>), the offending terrain is painted in an expanded cone around the groundtrack as shown in <figref idref="f0027">FIG. 37</figref>. In particular, the expanded cone includes a plurality of vectors, spaced, for example, 4° apart (configurable value), for ±90° (configurable value) on each side of the groundtrack. For an expanded terrain warning, the cone may be expanded out distances, for example as set forth in TABLE 7 below. <tables id="tabl0007" num="0007"><table frame="all"><title><u style="single">TABLE 7</u></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="32mm" /><colspec colnum="2" colname="col2" colwidth="47mm" /><colspec colnum="3" colname="col3" colwidth="21mm" /><colspec colnum="4" colname="col4" colwidth="47mm" /><colspec colnum="5" colname="col5" colwidth="21mm" /><thead><row><entry namest="col1" nameend="col5" align="center" valign="top">(For an expansion factor of 2)</entry></row><row><entry align="center" valign="top" /><entry namest="col2" nameend="col3" align="center" valign="top"><b>NORMAL</b></entry><entry namest="col4" nameend="col5" align="center" valign="top"><b>EXPANDED</b></entry></row><row><entry align="center" valign="top">TYPE OF THREAT</entry><entry align="center" valign="top">LOOK-AHEAD/LOOK-DOWN</entry><entry align="center" valign="top">LOOK-UP</entry><entry align="center" valign="top">LOOK-AHEAD/LOOK-DOWN</entry><entry align="center" valign="top">LOOK-UP</entry></row></thead><tbody><row rowsep="0"><entry align="center">ADVISORY</entry><entry align="center">1.0*</entry><entry align="center">2.0 *</entry><entry align="center">2.0 *</entry><entry align="center">4.0 *</entry></row><row><entry align="center" /><entry align="center">LAD</entry><entry align="center">LAD</entry><entry align="center">LAD</entry><entry align="center">LAD</entry></row><row rowsep="0"><entry align="center">WARNING</entry><entry align="center">0.5 *</entry><entry align="center">1.5 *</entry><entry align="center">1.0 *</entry><entry align="center">3.0 *</entry></row><row><entry align="center" /><entry align="center">LAD</entry><entry align="center">LAD</entry><entry align="center">LAD</entry><entry align="center">LAD</entry></row></tbody></tgroup></table></tables> The terrain threat is painted on the display 36 by calculating the starting vector angle DISPA as set forth in equation (29) when a terrain threat along the groundtrack has been detected. <maths id="math0029" num="(29)"><math display="block"><mi>DISPA</mi><mo>=</mo><mi>GROUNDTRACK</mi><mo>+</mo><mn>30</mn><mspace width="1em" /><mi>degrees</mi></math><img file="EP0925518B1_D0029.tif" /></maths>

0157The range increments are then determined as set forth in equations (30) and (31). <maths id="math0030" num="(30)"><math display="block"><mi>DELTA</mi><mspace width="1em" /><msub><mi mathvariant="normal">X</mi><mi>NS</mi></msub><mo mathvariant="normal">=</mo><mi>DISPLAY STEP SIZE</mi><mo mathvariant="normal">*</mo><mi>SIN</mi><mfenced><mi>DISPA</mi></mfenced><mo mathvariant="normal">*</mo><mn mathvariant="normal">1</mn><mo mathvariant="normal">/</mo><mi>COS</mi><mfenced><mi>LATITUDE</mi></mfenced></math><img file="EP0925518B1_D0030.tif" /></maths><maths id="math0031" num="(31)"><math display="block"><mi>DELTA</mi><mspace width="1em" /><msub><mi mathvariant="normal">Y</mi><mi>NS</mi></msub><mo mathvariant="normal">=</mo><mi>DISPLAY STEP SIZE</mi><mo mathvariant="normal">*</mo><mi>COS</mi><mfenced><mi>DISPA</mi></mfenced><mo>,</mo></math><img file="EP0925518B1_D0031.tif" /></maths> where the display step size is a configurable constant, for example 1/2 unit.

0158After the range increments DELTA X_NS and DELTA Y_NS are determined, the instantaneous position of the aircraft X₀, Y₀ is incremented by the same and the range is incremented by the step size. The altitude is then determined for the current position on the terrain map. Subsequently, the threat detection algorithms are run for the distances as discussed above. Should a threat be detected, a color is loaded into a RAM B, identified with the reference numeral 504 (<figref idref="f0024">FIG. 34</figref>). At the end of the yellow and red threat vectors 496 and 500, the vector angle DISPA is incremented by 4 degrees with the above steps being repeated until all vectors are done.

0159In addition to the terrain detection and display algorithms, a background terrain algorithm is also executed relative to the terrain data in RAM A (<figref idref="f0024">FIG. 34</figref>) once every screen update, about 4 seconds. The background terrain algorithm simply compares the current altitude of the aircraft with the various terrain altitude data in the RAM A. Any of the dot patterns 432, 434 and 436 (<figref idref="f0018">FIG. 25</figref>) resulting from the computation are stored in the RAM B (<figref idref="f0024">FIG. 34</figref>) organized in a wedding-cake configuration (<figref idref="f0021">FIG. 28</figref>), similar to the RAM A (<figref idref="f0024">FIG. 34</figref>). In addition, the RAM B is organized in single byte cells containing the color/dot pattern as discussed above. The closest runway may also be stored in the RAM B as either a circular or square pattern in, for example, a green color.

0160The color/dot pattern in the RAM B (<figref idref="f0024">FIG. 34</figref>) is transferred to the DSP memory RAM B*, identified with the reference numeral 506, once per screen update (i.e., 4 sec) and synchronized to the DSP screen painting period. If a threat is detected as discussed above, a new RAM B image is transferred to the DSP RAM B* 206 immediately. The image transfer time typically fits between two spokes. If not, the DSP will send an old spoke stored in its spoke memory 508, 510 twice.

0161The DSP then performs the RHO/THETA conversion on the data in RAM B* 506 using the updated aircraft position, selected range and heading data from the 486, as discussed below. Changes in the selected range are implemented on the fly (i.e., the next spoke) without restarting the sweep.

0162As discussed above, the ARINC 708/453 weather radar format consists of 512 spokes covering a configurable scan range of typically ±90 degrees and a range selected, for example, on a weather data control panel. For ARINC systems, each spoke is calculated twice, once for the left selected range and once for the right selected range. The spokes are then transmitted by way of a serial port 512, interleaved left spoke and right spoke to the bus switch 406 (<figref idref="f0016">FIG. 23</figref>) by way of a pair of drivers 514 and 516 (<figref idref="f0024">FIG. 34</figref>).

0163The resolution of the scan conversion is 256 range bins; each range bin actually being a double range bin per the ARINC 708/453 standard format. Such a resolution is sufficient for the dot patterns 518 illustrated in <figref idref="f0018">FIGS. 25</figref> and <figref idref="f0019">26</figref>.

0164The RHO/THETA conversion is best understood with reference to <figref idref="f0028">FIG. 38</figref>. In particular, once per sweep (i.e., every 4 seconds) the range increment R1 is determined. The range increment is the selected range divided by 256. After the range increment is determined, the starting spoke angle is determined (i.e., the aircraft true heading plus 90 degrees). Lastly, the incremental spoke angle DELTA ALPHA is determined. The incremental spoke angle DELTA ALPHA is the sweep range divided by 512; the total number of spokes per sweep.

0165Once per spoke (i.e., 512 times per 4 second sweep), the incremental range coordinates DELTA X_ns, DELTA Y_ns are computed and the spoke angle is incremented by DELTA ALPHA, as discussed above. The incremental range increment coordinates X_ns and Y_ns are computed in accordance with equations (32) and (33), respectively. <maths id="math0032" num="(32)"><math display="block"><mi>DELTA</mi><mspace width="1em" /><msub><mi mathvariant="normal">X</mi><mi>ns</mi></msub><mo mathvariant="normal">=</mo><mi>RI</mi><mo mathvariant="normal">*</mo><mi>SIN</mi><mfenced><mi>SPOKE ANGLE</mi></mfenced><mo mathvariant="normal">*</mo><mn mathvariant="normal">1</mn><mo mathvariant="normal">/</mo><mi>COS</mi><mfenced><mi>LATITUDE</mi></mfenced></math><img file="EP0925518B1_D0032.tif" /></maths><maths id="math0033" num="(33)"><math display="block"><mi>DELTA</mi><mspace width="1em" /><msub><mi mathvariant="normal">Y</mi><mi>ns</mi></msub><mo mathvariant="normal">=</mo><mi>RI</mi><mo mathvariant="normal">*</mo><mi>COS</mi><mfenced><mi>SPOKE ANGLE</mi></mfenced></math><img file="EP0925518B1_D0033.tif" /></maths> Once per range increment (i.e., 512 x 256 times per 4 second sweep), the instantaneous aircraft position X₀, Y₀ is incremented by DELTA X_ns and DELTA Y_ns. The color/pattern for the new position is looked up and output to the display 36. These steps are repeated for the rest of the spoke.

0166With the use of the DSP, two computations can be done in parallel. In particular, the color/pattern of each range bin including the conversion of the longitude to nautical miles (i.e., multiplication by 1/COS (LAT)). With reference to <figref idref="f0028">FIG. 38</figref>, this computation is performed relative to the true north map coordinate system with the aircraft true heading pointing up and the spoke angles referenced to it. The second computation addresses each display indicator pixel in order to paint the required dot pattern, relative to the display screen coordinate system.

0167The generation of dot pattern is illustrated in <figref idref="f0029 f0030 f0031">FIGS. 39-42</figref>. In general, the scan conversion is simulated by overlaying each RHO/THETA range bin over a calculated indicator pixel map which includes corrections for display aspect ratio using fractal 8 x 8 dot patterns.

0168The fractal patterns corresponding to the variable density dot patterns 432, 434 and 436 (<figref idref="f0018">FIG. 25</figref>) are illustrated in <figref idref="f0030">FIG. 40</figref>. The fractal pattern size is selected to be 8 x 8 pixels. The densest pattern 432 includes 1/2, or 32, of the 64 pixels filled. In the medium density pattern 434, 16, or 1/4, of the pixels are filled. In the least dense pattern 4 of 64, or 1/16, of the pixels are filled.

0169Referring to <figref idref="f0029">FIG. 39</figref>, the display screen is divided into a pixel map, consisting of, for example, 256 x 256 pixels 520. As mentioned above, the 8 x 8 fractal patterns are overlaid onto the pixel map to provide the variable density dot patterns 432, 434 and 436 illustrated in <figref idref="f0018">FIG. 25</figref>.

0170The computation of the dot patterns is best understood with reference to <figref idref="f0029">FIGS. 39</figref>, <figref idref="f0031">41 and 42</figref>. Using an X-Y coordinate system, each spoke is determined at the bottom center of the display 36. For a 256 x 256 pixel display, the starting point corresponds to 120, 0. Next, the incremental X and Y coordinates relative to the scan angle ALPHA are determined according to equations (34) and (35). <maths id="math0034" num="(34)"><math display="block"><mi>DELTA X</mi><mo>=</mo><mi>RI</mi><mo>*</mo><mi>SIN</mi><mfenced><mi>ALPHA</mi></mfenced></math><img file="EP0925518B1_D0034.tif" /></maths><maths id="math0035" num="(35)"><math display="block"><mi>DELTA Y</mi><mo>=</mo><mi>RI</mi><mo>*</mo><mi>COS</mi><mfenced><mi>ALPHA</mi></mfenced><mo>,</mo></math><img file="EP0925518B1_D0035.tif" /></maths> where RI = a double range increment and ALPHA = the scan angle. These increments are added to the starting point 128, 0. After each increment, the selected fractal pattern is looked up in the RAM B * to determine whether to send the last range increment as either block or colored.

0171Referring to <figref idref="f0031">FIG. 42</figref>, the increments DELTA X and DELTA Y are applied to a pair of adders 522 and 524, which, in turn, are applied to an X-counter 526 and a Y-counter 528. As discussed above, the X counter 526 starts at 128 and counts up if DELTA X > 0. If DELTA X < 0, the X-counter counts down, indicating a pixel on the left (<figref idref="f0029">FIG. 39</figref>) of the starting position. Since the starting pixel is at the bottom of the display, the Y-counter 528 starts at zero and always counts up. A range counter 528 counts up to 256; the maximum number of double range bins. The range counter 528 is driven by an increment clock such that the pattern/color is painted once per range increment. The output indicates whether the pixel is to be painted left black.

0172As discussed above, the background terrain algorithm determines the particular fractal pattern to be used. This data is applied to the DSP to enable the pixels on the pixel map (<figref idref="f0022">FIG. 29</figref>) to be painted or left blank. In particular, the three least significant bits from the X counter 526 and Y counter 528 are used to select the fractal pattern as shown.

0173The scan conversion simulation is performed in parallel with the RHO/THETA conversion which determines the color/dot pattern of the range bin to be transmitted. The color/dot pattern provided by the RHO/THETA conversion is then overlaid with the fractal dot pattern to determine whether the range bin color should be black or painted to provide the proper dot pattern on the display.

0174Referring now to <figref idref="f0035 f0036 f0037 f0038">Figs. 46-48</figref>, another embodiment of the present invention incorporating a pop-up terrain display based on the aircraft's will now be described.. The system may be selected by the pilot, or it may automatically assert the visual terrain display whenever a terrain threat is present if the display mode is compatible. In an exemplary embodiment, the display changes color according to the relative severity of the terrain threat. The display is most commonly shared with an existing display such as weather radar, and the pilot or alerting system must decide which information is displayed. Of course, it will be recognized that the terrain display may be separate from existing displays such that the terrain information can be displayed at all times during flight.

0175Referring to <figref idref="f0035">Fig. 46</figref>, the terrain information system includes a WXR switch (not shown) and a terrain momentary selection switch (not shown) for switching between three different terrain/weather selection states. In a first state 806, terrain and weather information are not displayed, and the switching relays (not shown) are deasserted. This state is true if the WXR switch 802 is in the OFF position or deasserted. In this state 806, toggling the terrain momentary selection switch 804 has no effect on the state of the system. Turning the WXR switch 802 ON moves the system into a second state 808 in which the weather display is enabled and displayed when the EFIS display modes are selected on the Mode Control Panel. The switching relays are not energized in this state 808. In the second state 808, toggling or pressing the WXR switch 802 back to the OFF position moves the system back into the first state 806. Toggling or pressing the terrain selection switch 804 moves the system into a third state 810 in which the terrain display is enabled and terrain is displayed if the EFIS display modes are selected on the Mode Control Panel. In the third state 810, the switching relays are energized and the terrain ON lamp is illuminated. This third state can be asserted only from the second state 808 by asserting terrain selection switch 804. Pressing or toggling the WXR switch 802 to an OFF position returned the system to the first state 806.

0176The navigation (NAV) display modes for display of weather or terrain are: center map, map, EXP. VOR/ILS, Exp. NAV. Neither weather or terrain is displayable in the following modes: Plan, Full NAV, Full VOR/ILS.

0177<figref idref="f0037">Fig. 48</figref> illustrates a block diagram of the terrain display selection annunciation and fault logic according to one embodiment of the present invention. As shown, the system includes a Captain's Control Panel 900 with a Captain's Terrain Selection Logic 902 and a First Officer's Control Panel 904 with a First Officer's Selection Logic 906. The system further includes first and second relay assemblies 908, 909 (e.g., 453 bus relays) coupled to the selection logics 902, 904 for activating a number of lamps in the Captain's and First Officer's lamp assemblies 910, 912, respectively. Of course, it will be recognized by those in the art that a single selection logic and relay assembly may be used for both control panels and lamp assemblies.

0178Lamp assemblies 910, 912 each include a terrain display select/deselect push button/Lamp combination switch 914, 916 that provides a momentary input to terrain selection logics 902, 906. Switches 914, 916 allow the pilot to manually select the terrain display, as discussed above. Control mode data is provided to selection logics 902, 906 from a 429 bus 920 in each Control Panel 904, 906.

0179Lamp assemblies 910, 912 further include an terrain ON lamp 922, a terrain displayed lamp 926, a terrain alert lamp 930 and a terrain not available lamp 940. Terrain ON lamp 922 indicates that the terrain display has been selected, and is driven directly from a switching relay 924 within relay assemblies 908, 909. When relay 924 is energized, the ON lamp 922 is illuminated. Terrain ON lamp 922 function is analogous to the WXR On lamp when the weather display is activated. The terrain ON lamp 922 does not indicate that actual terrain information is being displayed unless the proper NAV display modes have been selected. Terrain displayed lamp 926 indicates to the pilot that the terrain has been selected and is actually being displayed. Terrain displayed lamp 926 is the analog of the WXR mode annunciation on the NAV display when weather is displayed. Terrain displayed lamp 926 is illuminated when the terrain display 932 is valid, the correct NAV display modes (received via 429 bus label 272/273) have been selected and a relay 928 from relay assemblies 908, 909 has been energized. The terrain displayed lamp 926 is preferably wired through a spare relay contact to ensure relay failures are indicated correctly.

0180Terrain alert lamp 930 indicates that a terrain alert is present so that the pilot is aware that the terrain information should be displayed. As shown in <figref idref="f0037">Fig. 48</figref>, terrain alert lamp 930 receives terrain input signals based on the terrain display input parameters 936, which may include latitude and longitude of the airplane, groundspeed, track angle, roll angle, vertical speed, barometric altitude, etc (discussed in detail above). Terrain not available lamp 940 is illuminated if any output signals required for the terrain display and terrain alerting system are invalid, or if the NAV solution is inaccurate, or if a EGPWS hardware failure affecting the terrain display is present. In addition, terrain not available lamp 940 is illuminated when the EGPWS power is OFF.

0181The terrain system further includes a number of lamps associated with the aircraft's ground proximity warning system (GPWS), which is described in greater detail above. A GPWS INOP lamp 950 is illuminated if the GPWS computer has failed or if any input signals are invalid. A GPWS ALERT lamp 952 is illuminated during a GPWS alert, and during a terrain alert. A GPWS PULL UP lamp 954 is illuminated during a GPWS warning, and during a terrain warning.

0182The present invention further includes a terrain display fault indication 960, as shown in <figref idref="f0037">Fig. 48</figref>. In this embodiment, any detected fault affecting terrain display (e.g., terrain display invalid), causes the terrain not available lamp 940 to be illuminated (a flow diagram of this logic is shown in <figref idref="f0036">Fig. 47</figref>). Preferably, the relay state is not changed without pilot action so that a blank display 962 is shown when the terrain display 932 is selected (i.e., neither a red "X" or a moving marker is displayed). The terrain display 932 is invalidated if any of the required terrain display input signals 936 are invalid or if the internal EGPWS fault monitor 964 detects a hardware fault. Of course, it will be recognized that some of the input signals may have default values if not selected, which would prevent invalidation of the terrain display. In an exemplary embodiment, for example, the selected range inputs are defaulted to 10 nautical miles if one range input is invalid or not present (i.e., no fault indication is given).

0183The terrain selection logics 904, 906 each include a latch 980 for moving the relay assemblies 908, 909 into energized and deenergized states. As shown in <figref idref="f0037">Fig. 48</figref>, logics 904, 906 include an AND gate 982 which has a positive output if the WXR control is ON, and if the terrain selection switch 914, 916 are pressed. If the latch 980 is in the deenergized state, AND gate 984 will send a positive signal to set the latch 980 to the energized state. In this state, the relay assembly 908 is energized, and the terrain selected and terrain displayed lamps 922, 926 are illuminated. Similarly, if the power is ON, or the WXR is OFF or if the relay has been energized and the terrain selection switch 914, 916 has been pressed, the latch 980 will be reset to the deenergized state. In this state, the terrain selected and terrain displayed lamps 922, 926 are turned OFF.

0184Obviously, many modifications and variations of the present invention are possible in light of the above teachings. Thus, it is to be understood that, within the scope of the appended claims, the invention may be practiced otherwise than as specifically described above.

<u style="single">APPENDIX 1</u>

TERRAIN DATA STRUCTURE

0185The terrain data format is made up of four distinct parts, the file indirection table, the group index, the data index, and the compressed elevation data itself. See <figref idref="f0032">Figure 43</figref>. Notes: All multi-byte numbers are stored Little Endian (LSB...MSB).

File Indirection Table

0186A file indirection table is provided at the top of the index file block. The file index field is each data index is used to index into this table to get the ASCII filename for file system usage. This filename is then mapped be the file system to the correct 128kbyte physical memory block. A time stamp is provided for each file, this can be used to match to the time stamp in the file header for each file. This helps check for consistency between the index files and the data files. <tables id="tabl0008" num="0008"><table frame="none"><tgroup cols="2" colsep="0" rowsep="0"><colspec colnum="1" colname="col1" colwidth="25mm" /><colspec colnum="2" colname="col2" colwidth="64mm" /><tbody><row><entry>SIZE:</entry><entry>16-bits - unsigned</entry></row><row><entry /><entry>Number of filenames defmed:</entry></row><row><entry namest="col1" nameend="col2" align="left">Repeated for each filename</entry></row><row><entry>FILENAME:</entry><entry>16-bytes - ASCII</entry></row><row><entry /><entry>ASCII filename with the extension included</entry></row><row><entry /><entry><Drive>:<File 8>.<Ext 3></entry></row><row><entry>TIME STAMP:</entry><entry>32 bits - unsigned</entry></row></tbody></tgroup></table></tables>

0187The number of seconds since January 1, 1980 (DOS's time base).

Group Index

0188The group index is used to reduce the number of data index's which must be searched to find the area of interest. The group indexes are built by dividing the earth into TBD bands or geographic areas, current thinking is 10 degrees bands. These bands need not be of equal size. The group index points to the corresponding data index block to use for this area. The area covered by the group must be a rectangle. Note: the group index and data index occupy the same file.

Definitions:

0189<dl id="dl0001"><dt>SIZE:</dt><dd>16-bits - unsigned Number of groups</dd></dl> The following is repeated for each group defined: INDEX:  32-bits - unsigned <tables id="tabl0009" num="0009"><table frame="all"><tgroup cols="2" colsep="0" rowsep="0"><colspec colnum="1" colname="col1" colwidth="18mm" /><colspec colnum="2" colname="col2" colwidth="17mm" /><tbody><row><entry colsep="1">12 - bits</entry><entry colsep="1">20 - bits</entry></row><row><entry colsep="1">File Index</entry><entry colsep="1">offset</entry></row></tbody></tgroup></table></tables> File Index: Always 0 for group index offset: offset in bytes to start of index file for this group LAT MAX:  24-bits - signed NW corner latitude Scaling: LSB = 180 / 2^23 Sign: (+) = North, (-) = South LAT MIN:  24-bits - signed SE corner latitude Scaling: LSB = 180 / 2^23 Sign: (+) = North, (-) = South LON MIN:  24-bits - signed NW corner longitude Scaling: LSB = 180 / 2"23 Sign: (+) = East, (-) = West LON MAX:  24-bits - signed SE corner longitude Scaling: LSB = 180 / 2^23 Sign: (+) = East, (-) = West Memory Usage:  16 bytes per group Assume 60 groups worst case = 9600 bytes

Data Index

0190The data index contains information defining the coverage, resolution and storage location of each individual terrain data file. The data indexes are grouped together per the group index defined above (i.e., group 2's data indexes immediately follow group 1's data indexes).

Definitions:

0191<dl id="dl0002"><dt>SIZE:</dt><dd>16-bits - unsigned Number of index's</dd></dl> The following is repeated for each index defined: INDEX:  32-bits - unsigned <tables id="tabl0010" num="0010"><table frame="all"><tgroup cols="2" colsep="0" rowsep="0"><colspec colnum="1" colname="col1" colwidth="18mm" /><colspec colnum="2" colname="col2" colwidth="17mm" /><tbody><row><entry colsep="1">12 - bits</entry><entry colsep="1">20 - bits</entry></row><row><entry colsep="1">File Index</entry><entry colsep="1">offset</entry></row></tbody></tgroup></table></tables> File Index: Filename index for data offset: offset in bytes to start of index file for this group LAT MAX:  24-bits - signed NW corner latitude Scaling: LSB = 180 / 2^23 Sign: (+) = North, (-) = South LAT RES:  8-bits - unsigned Number of rows per 1 degree LAT MIN:  24-bits - signed SE corner latitude Scaling: LSB = 180 / 2"23 LON MIN:  24-bits - signed NW corner longitude Scaling: LSB = 180 / 2^23 Sign: (+) = East, (-) = West LON RES:  8-bits - unsigned Number of columns per 1 degree LON MAX:  24-bits - signed SE corner longitude Scaling: LSB = 180 / 2^23 Memory Usage:  18 bytes index (terrain file) Assume 8000 files = 144k-bytes

Data

0192The following defines the format for the actual elevation data. The index files above are used to locate the desired elevation file, and also contain information on the location, resolution etc.

0193Note: for all data types below, the return elevation values at 16 bit singed values.

0194First value is always compression type as defined below: <ul id="ul0002" list-style="none" compact="compact"><li>COMPRESSION TYPE:8 bits - unsigned Compression: <ul id="ul0003" list-style="none" compact="compact"><li>0 = no compression</li><li>1 = HUFFCODE - type 1</li><li>2 = HUFFCODE - type 2 (currently not used)</li><li>3 = Run length (currently not used)</li><li>4 = water or constant elevation</li></ul></li></ul> For type 0, no compression, the uncompressed elevation values are stored directly.

0195The following definition is for the HUFFCODE compressed data format: <ul id="ul0004" list-style="none" compact="compact"><li>The top of the data contains the decode tables needed to uncompress the terrain data.</li></ul> Note: Decompression needs the Number of X and Y values of decoded file (from index file). Decode tables: <dl id="dl0003"><dt>FIRSTCODE SIZE:</dt><dd>1 byte - unsigned Size in bytes of firstcode and symbol pointer tables (max 32 bytes)</dd><dt>FIRSTCODE TAB:</dt><dd>n bytes - unsigned Firstcode table</dd><dt>SYMBOL POINTER TABLE:</dt><dd>n words - unsigned Used with Firstcode table to index into symbol table.</dd><dt>SYMBOL SIZE:</dt><dd>1 word - unsigned Size in bytes of symbol table (max 256 words)</dd><dt>SYMBOL TAB:</dt><dd>n words - signed Symbol table. Contains decoded symbol</dd><dt>FIRSTVALUE:</dt><dd>1 word - signed First elevation value. Used by linear predictor</dd></dl> Note: typical words decode table size is less than 150 bytes. <dl id="dl0004"><dt>ELEVATIONS:</dt><dd>Compressed elevation values in binary</dd></dl>

0196Note: padded to align on byte boundary currently. For Compression type 4, (water/same elevation value), only 3 bytes are needed. The compression type byte defined above, and the 16 bit elevation value. This is used for files that contain only one elevation value. This will typically be used for water and "no-elevation" areas.

Build System

0197<figref idref="f0033">Figures 44</figref> and <figref idref="f0034">45</figref> show a preliminary look at the steps to build the Terrain database from the original source data. Some key points: <ol id="ol0001"><li>1. Currently basic input file is 1 degree at 30 arc-seconds. This matches NOAA data. System should be built to handle any filesize.</li><li>2. Build system uses a control file automatically generated from runway database and a NOAA provided worldwide 1 degree terrain file. Control file contains a list of every 1 degree quad which is needed (i.e., on land), and what the desired resolution of that 1 degree file should be. Resolution is determined by proximity to an airport and the maximum elevation value of that one degree file. 1 of 3 resolutions are currently output, 30, 60, and 120 arc-seconds. In addition, the control file will indicate which 1 degree quads can be combined into one file. This is done for lower resolution quads in an attempt to keep the filesizes even and reduce the number of files.</li><li>3. Enroute data is at 5 minute resolution. Enroute data is separate from hi-resolution data above (overlaps), and is grouped in 10 degree by 10 degree blocks.</li><li>4. All files have the runway extrusion algorithm run on them. Runway extrusion should be run at the resolution that the data will be used in the system, it does not need to be run on the enroute data.</li><li>5. Some files will need a special extrusion run on them. This will be defined in a separate control file.</li></ol>

0198Grouping and packaging control files are needed to ensure that the build system groups the files into the correct structure.

Notes on "Terrain Data Embedded Format"

0199Maps are read from the file system, decompressed, and then used to update "RAM A" as the aircraft moves across it. Due to the slow speed of the decompression, two RAM A's are used. The detection and display algorithm uses one RAM A while the file picker updates the other RAM A with the latest maps. After the file picker has updated its RAM A, the RAM's are swapped and the process starts over.

0200The file picker process starts by the TAD asking the file system to open a file with a known, fixed name ("TERRAIN.IDX"?). This file contains the file indirection table, the group index table, and the data index tables. The TAD then periodically scans the group index table to find map groups that cover the general vicinity of the aircraft. This scan is done once for each source set of maps (currently two: hi-res and en-route). This step quickly reduces the number of maps that need to be considered for decompression from several thousand down to a few hundred. It identifies candidate groups by comparing the group bounding box stored in the group index table entry against the aircraft's current position. The number of group index table entries is stored at the top of the table, and all the entries immediately follow.

0201When a candidate group is found, the TAD uses the < fileindex > in the group index table entry to locate the file containing the data index tables. If the number is zero (the usual case), it indicates that the data is in the main index file (the one we are already reading). A future option would be to make the number non-zero. It would then identify a secondary index file named in the file indirection table. Once the file containing the data index table is opened, the TAD seeks to the < offset > in the group index table entry and starts reading the data index table for that group.

0202Each data index table begins with its number of entries. The TAD scans each of the data index table entries looking for any of the maps needed for the current aircraft position.

0203If a data index table entry is found to describe one of the maps we need, the TAD uses the < fileindex > to get the name of the map package file from the file indirection table. A map package file is a concatenation of a couple hundred compressed maps. From there it reads the compression type and calls the correct decompression algorithm to decompress the file. The decompressed data is then up-sampled or down-sampled into the appropriate tiers of RAM A. When up-sampling (filing higher resolution tiers with lower resolution data), the highest resolution data is used. To support this TAD must keep track of the current resolution used to fill each cell of RAM A. For down-sampling the highest elevation is always used. The enroute tier of RAM A is always filled with only enroute data. When all map files within the aircraft's area have been processed RAM A is released.