METHOD AND APPARATUS FOR ESTIMATING SLANT VISIBILITY IN FOG

Abstract: Apparatus for estimating the slant visibility in obscured atmospheres such as fog, whereby an indication is derived as to when, during an aircraft landing approach in such conditions, the aircrew are likely to be able to see the runway lights. The apparatus comprises a transmissometer array arranged to take visibility readings over a slanted distance from different heights not far from a runway threshold, and a computer arranged to compare the readings with a plurality of fog profile models, to determine the definitive model and hence to derive and indicate a visual acquisition estimate.

Description: METHOD AND APPARATUS FOR ESTIMATING SLANT VISIBILITY IN FOG The present invention relates to apparatus for estimating the slant visibility in obscured atmospheres such as fog, falling snow and rain, and in particular for providing an indication for use by aircrew and air traffic controllers as to when, during an aircraft landing approach in fog conditions, the pilot is likely to be able to see the runway lights and the extent to which the lighting pattern or marking will be in view at any instant. The mode of development of this extent during the course, say, of an aircraft landing run, is called a visual sequence. Hitherto operational measurements of low visibility for assisting air traffic control have almost always been made by human observers or instruments, such as transmissometers, close to the ground. These measurements are reported as 'runway visual range*. This technique cannot usually give a reliable indication of any vertical element to visibility, and this vertical element has a profound effect upon what is called the slant visual range, that is for example, the distance an aircraft pilot can actually see when looking toward the ground, be that anywhere between directly below and the horizon. The ground measurement derived runway visual range then is used to estimate the probability of a landing success. The shortcomings of this method tend to be in favour of safety, that is, to prevent visual landing attempts more often than might have been possible rather than less. This is because the majority of low visibility situations are those of deep mature fog or low cloud, and since fog density increases with altitude (up to a reletively instantaneous ceiling), once a pilot has made visual contact with a runway lighting pattern .he will be able to maintain it down to the landing. This method does not, OMPI ., of course, endeavour to imply any actual estimate either of the visual contact height or of the visual sequence. Moreover it can make no allowance for shallow, developing, dispersing or lifting fogs, which although being perhaps less prevalent than deep mature fogs are nevertheless frequently encountered. Drawbacks also exist in the measuring technique per se. This h customarily comprised deploying transmissometers with horizontal measuring paths at a standard measuring height of 5 u above the groun (5 ft in the UK). Shallow ground radiation fog the ceiling of which is below transmissometer level will not be detected. A fog ceiling slightly above transmissometer level will result in an unduly severe fog indication and an unnecessary refusal being issued to the pilots of tall aircraft, (eg those where the cockpit is over 10 m from the ground). Furthermore a fog ceiling intermittently penetrating the measurement path of the instrument will result in large and quite illusory fluctuations of the reported values. It will also be appreciated that there are, to say the least, constraints pertaining to the deployment in the vicinity of airfields of instrumentation operating at significant heights, especially heights where at instrumentation could be employed to simulate meaningfully slant visual range. What has now been discovered is that by making a small number of slant or vertical measurements (typically three in number) from vario heights to the ground sufficient information is available which by the judicious application of predetermined models of mean extinction coefficient profile allow estimation of the actual profile over a ran of heights not less than that of the measurements with a useful accuracy. (Extinction coefficient (r) is the reciprocal of the distance over which light falls, due exclusively to obscuration, to 1/e times its original intensity. Accordingly, in those circumstances where the fog or cloud extend considerably above the ground (several hundreds of feet) then the mean extinction coefficient profile in the bottom one hundred or so feet may be adequately described by the use of a polynomial of order two or less such that OMPI /., IPO where -T is the mean integral of extinction coefficient between the n ι ground and a height = h, o" is the extinction coefficient at a height h a, b and β are constants accordingly 3 values of σ" at the heights h , h_, h., will determine a, b and c t In those cases where the upper extent of the fog or cloud is below the upper limit of measurement one may proceed by determining the height of the nominal fog top in relation to the heights of the slant or vertical measurements and then to model separately the profile of the mean integral of extinction coefficient both below and above the fog top. Although the three measurements are apportioned between the two now independent models of mean integral of extinction coefficient, it has been found that models of a lesser complexity are' of value in those cases, viz: δ L = J'h + k below the fog top and above the fog top, where j, k, p, q and r are constants. Where the fog top lies between the upper two fog measurements only qne point is available for the calculation of <? and an assumption u can be made about r which generally improves the performance of the model compared to removing r from the model. In the case where the fog top lies between the lower pair of measurement heights (upper extent of individual slant or vertical) only one value is available for the calculation of σ * and in this case h L the term jh is ignored. The cusp resulting from the intersection of these two models at the nominal fog top may additionally be abated. According to one aspect of the present invention therefore • apparatus for estimating a fog structure comprises a transmissometer array disposed for slant operation at different heights above the ground, and means for employing signals derived from the transmissomete array in estimating slant visual range. By slant operation is meant that the points between which the extinction coefficient are measured are one above, but not necessarily vertically above, the other. A transmissometer generally comprises a light source and a detector of light intensity over a given distance from the source. Rather than have, however, a source in one unit and a detector in anoth it is convenient if the said other unit is a reflector. In this manner the source/detector unit can readily be arranged so that the detector unit compares the intensity of both the emitted and the reflected light. This arrangement is particularly suitable in the present invention, as it permits a transmissometer source/detector unit array to be mounted close to the ground, where it can readily be serviced, and reflectors to be mounted at different heights. The source/detector unit array may comprise a single unit, a composite of several units, or several (in particular three) distinct units. It is an important feature of the present invention, from the point of view particularly of cost , that the apparatus may be arranged to provide readings at three different heights, and preferably comprises three reflectors, one at each of the heights, and a source/detector unit at ground level. It is preferred not to use a singxe reflector traversing the heights, for simultaneous readings are preferred and the traverse movement will usually give rise to eddies which will affect the readings. Means for keeping the source and/or detector clean may also be provided. These preferably comprise air curtains and a shield tube, though washing and wiping means may additionally or alternatively be employed. It may be preferred, for the purpose of both realism and equipment hygiene, for the light paths not to be vertical. According to a feature of the invention the means for employing signals derived from the transmissometer array in estimating the slant OMP /,, WIP visual range may comprise means for converting the transmissometer • signals into signals representing the mean integral of extinction coefficient (<? ) over the associated distance, means for comparing the &, signals and selecting an appropriate model from a predetermined set of models, means for determining the values of the constants in the selected model, and thence the definitive models, and means for deriving an estimate of the mean integral of extinction coefficient profile from the definitive models. The means for converting the transmissometer array signals into cr, signals may employ the function: where R is the distance over which the measurement is taken, I, is the signal from the transmissometer, and I is a clear day datum transmissometer signal. The predetermined set of models may comprise: - 2 - - - ' O" . = ah + bh + c for use when cr, > σ, > σ, , n h h 2 h and j(h_ or h ) + k for use when ?, < σ > <f qh, + r •> ' 2 V h, h_ h„ , and j i> \ and k when σ, < σ, < , . q(h, or h + r h 3 h 2 h 1 3 2 Typically, heights from the ground, over which the I, signals can be derived, may be k, 12, 30 metres, so that the reflector units may be mounted on a tower 30 metres high, which, in accordance with the above- mentioned current regulations, may be sited 370 m laterally from the centreline of any airfield runway. The <^ profile may be extrapolated to a small amount, eg 10-20 above the maximum height (h_). A lower measuring base may be sited closer to the runway for deriving an I, signal at, say, metres height to provide some confirmation of fog extent and consistency and to improve the estimate of the o " profile, h OMPI /-, WIPO ,* The apparatus may further have means for determining slant visua range and visual sequence from the cr profile, that is to say the distance from and position with respect to, the runway lighting array at which the pilot will see some part of that array, how much he shou be able to see, and how his acquisition will develop in the course of the landing (visual sequence). The determination of a visual contact height (the height at whic a descending aircraft pilot will first see the landing lights) and a visual sequence, requires consideration of, in addition to the fog structure, the lighting intensity, pattern and direction, the pilot's visual threshold, the structure and configurat on of the aircraft, th amount and constitution of day or night light, and the aircraft cabin lighting. It may be possible to derive for certain airfields a function: where H is the visual contact height, c o " is the mean integral slant extinction coefficient at that he o Y and ι are empirical constants, related to flight path and lighting pattern characteristics respectively, and EL, is the eye illuminance threshold (measured in lux). This function was derived from Allard's law which states: where I is the intensity of a point source of illumination, and R is the visual range to extinction, by determining a number of visual sequence patterns for a particular flight path in different fog and eye illuminance threshold situations. The apparatus may therefore have means for determining the inter section of the H function and the definitive model, this intersectio c occurring at H . For circumstances where the eye illuminance c threshold is not known the apparatus may have means for deriving Log 0 (I ) from the function: - ■ siEA OMPI /., WIPO Log 1Q (^) = u log 1Q β - v where β is background brightness in nits and u and v empirical constants typical values for which are 0.675 and 5 * 7 respectively. β may be derived from a detector mounted at the top of the transmissometer tower. In these circumstances the apparatus may be arranged to feed only Log (E_) * - w into the H function, w being an accustomation factor or night use, typically 6.1. The apparatus may additionally have means for storing aircraft configuration information, in particular information on the extent to which a particular aircraft in a landing run configuration obstructs the downward view of a pilot, applying this information to the derived H c value, and modifying H accordingly. C It may be more appropriate in certain cases for the apparatus to have means for determining the slant visual range and the visual sequence rather than or in addition to H . A suitable method of determining slant visual range is iterative, involving applying the derived definitive cr model, the particular lighting pattern, direction and intensity and visual threshold values at a sequence of aircraft positions during a landing run. This method is described in Royal Aircraft Establishment Tech Memo Avionics 18 (BLEU) entitled "A Method of Computing the Visual Range of Lighting Patterns During the Approach to Landing" by A Puffett, August 1 7--+. The apparatus may therefore have means for carrying out this method. A second aspect of the invention provides processes for estimating a fog structure, and perhaps visual ' contact height, slant visual range and visual sequence, by operation of the means specified above. By virtue of the present invention, then, a much more accurate assessment can be obtained than hitherto available of the slant visual range and visual sequence of the pilot of an aircraft approaching a runway in order to land. Moreover it has been found that significant discrepancies between the estimated values and the values in practice occur only briefly during transient conditions such as fog decay. Thus it should be possible, by means of the invention, to permit aircraft to land much more frequently in fog conditions than hitherto. OMPI /., WIPO . Moreover the apparatus can readily be adapted to give visibility • indications or such other obscuration conditions as snow-and rain. Apparatus for estimating slant visibility in fog will now be described by way of example with reference to the accompanying drawings, of which: Figure 1 is a schematic arrangement of a transmissometer array, Figure 2 is a plan view of a runway threshold environment including the transmissometer array, Figure 3 is a flow chart depicting operation of the apparatus, Figure k is a graph comparing experimentally derived and model extinction coefficient curves in shallow fog, Figure ka. is a graph demonstrating the deviation of the shallow fog model from the empirical values, Figure 5 is a graph comparing experimentally derived and model extinction coefficient curves in deep fog, Figure 5a is a graph demonstrating the deviation of the.deep fog model from the empirical values, Figure 6 is a generalised computer program flow chart for determini visual sequence, and Figure 7 is a tabulated output example of the program of figure 6. Figure 1 shows the principal data source, the Slant Visual Range (SVR) transmissometer array tower 10 carrying three transmissometer reflectors 11 at k, 11 and 29 metres respectively, and a background brightness detector 12 at the top of the tower 10. Based on the ground a short distance from the tower 10 is a transmissometer 12 of the source detector unit type. The tower also carries obligatory warning lights 1 As illustrated in figure 2 the tower 10 is located 370 metres laterally from the centre of a runway 20 and 350 metres ahead of the threshold thereof. A runway visual range (RVR) transmissometer array 21 is located 190 metres from the runway centre line and behind an approach lighting array 22. It is disposed for level operation at 5 metres height. A OMPI WIPO computer for processing signals from the SVR and RVR detectors is located in a control tower 23. The process by which the computer derives visibility information using the RVR and SVR inputs is illustrated in figures 3 to 5« As shown in figure 3 SVR transmissometer signals are converted to extinction coefficient signals f (k m), cr (11 m) cr (29 m) in or element 3 of the computer. The computer has been preprogrammed with extinction coefficient (σ;n) function models as follows: -cr - h2 + bh + c for use when -cr > o mZ- > σ -m- , h h, 2 h ■ , ■ "■ and j(h_ or h . ) + k for use when cr < cr > cr and qh^ + r 2 1 h . h 2 1^ ' —q(rrh or-P h, ≤\ + r and k when cr.h,■< σ:h ≥ < cr.. For deep fog conditions, such as illustrated in figure k, when the second order polynominal applies, a signal cr -^ derived (at 32) from the RVR detector 21 is employed to determine the constant -c'. Figure 5 clearly illustrates the second of the above conditions, shallow fog, when two different expressions for cr are applied. In the comparator 31 the computer compares crn with cΛ ' and selects the worst case cr,h ,' selects in accordance with the relation wc one to another of cr, , σ, and σ, , which of the above models is h h_ h.., wc 2 3 appropriate, and determines the remainder of the constants (including assuming a value for the constant r in cases where the fog ceiling is between 11 and 29 ω) and hence a definitive c profile. From the background brightness detector 12 is derived a signal which is converted (at 33) into a function of eye illuminance threshold EL according to the formula: Log 1Q ^ = u log 10 β - v OMPI The inventor used the value O.675 and 5.7 for the empirical constants u and v respectively, but others skilled in the art are known to use higher values. A discriminator 3**+ then determines whether Log..J- should be i passed on as such or limited to a minimum value w quantifying night time accustomisation factor. The inventor uses w = -6.1 but others skilled in the art may prefer a different value. The chosen signal for ^o ,.-. Ϊ , is fed into a visual contact height H ffuunction determining element 35» The visual contact height 0 function is: YClog 0 E^ * H = " as mentioned earlier in the specification, C oh c where and « ■ < , are empirical constants representing a particular flight path and lighting pattern (37)• I one context their values were 288.7 and 1.3 respectively. Together with the definitive σ profile from the comparator element 5 31 the visual contact height function H is fed into a visual contact c height estimating element 38-1 and the intersect, which occurs, at the definitive H estimate, is fed to an information output 39 . thus c providing a rapid, relatively inexpensive, and in some cases sufficient indication of when, during an approach to land, an aircraft pilot can 0 expect to see the runway lighting. As an adjunct or alternative to estimating a visual contact height using the function for Hc and the elements 35 and 38 as described above, the apparatus has a visual sequence estimator element hθ for providing an estimate of the development of a pilot's visual acquisition of a runway and its lighting pattern -during the course of a landing run. In one example the estimator element kO comprises an Elliot ^100 computer to which can be configured a visual sequence program written in Algol 6θ as now described below. In order to simulate the lighting pattern and the conditions under which it is viewed, a number of models are needed to define the various parameters. These fall into two categories; those representing variables which are normally continuous functions such as the distributi of light in the beams of the different airfield lights, the vertical - ■ RE OMP ∑ l profile of fog density, and the path in space traced by the pilot's eyes. Each of these variables is modelled by a series of straight lines as detailed in Appendix D below.. However it is necessary to use discrete values to describe the configuration of the ground lighting system, ie the disposition of the various types of lamps on the ground, and their associated setting angles. The adjustment of the intensity of portions of the lighting system has to be similarly described. In the present version of the program the lights are assumed to have two optical axes of symmetry and the distribution of light to be such that it may be represented in an isocandle diagram by a series of concentric ellipses of constant eccentricity. The consequence of these simplifying assumptions is that a lamp is fully described by one half of its vertical-axis light distribution and the ratio of the major and minor axes of the ellipse. Alternatively other beam characteristics such as may be represented by rectangles or ellipses of varying eccentricity can be easily incorporated. In the present scheme the light output characteristics .of each type of lamp is defined, as described in Appendix A, by five segments and in the examples considered this has adequately described the light distribution down to approximately two orders of magnitude below the peak value. This level of light has been found to be equivalent to off-axis, viewing angles far in excess of those which are desirable in practice. No restriction has been placed on the number of lamp types so that a virtually unlimited number of types may be read into the program and called upon either in mixed lighting patterns or in differing successive cases. Since this program is concerned with the conspicuity of the lamp fittings rather than with the information contained within the pattern itself, only a single line of lights is considered at any one time. That is to say a single line of lights lying in the ground plane and parallel to the centre line or X axis. However, where a comparison is required between the visual sequences produced by lights at differing lateral positions (ie runway centreline and runway edge lights) this may be achieved by the execution of a succession of cases, as is made clear below. The various parameters associated with the lamps are their type, ■ elevation and azimuth setting angles, and the intensity factor. This latter parameter allows the simulation of sub-standard sections of the lighting pattern as will be found necessary in practical examples. However, all the parameters are independent in terms of their applica¬ bility to a given longitudinal section of the pattern. Consequently each parameter must be separately defined over the whole length of the pattern. This is achieved by defining the downwind end of the pattern and assigning a value to the parameter. This value is then assumed to hold for increasing upwind positions until restated, or until the upwind end of the pattern is reached. In this way a series of blocks of values for a parameter may be built up, the only restriction being that the minimum distance over which a parameter can be assigned a particular value is 10 metres. This distance is the size of the step currently in use for scanning along the lights, thus any block of smaller length might be passed over without its presence being recognise Whilst smaller step intervals may be used the resulting increase of computational time may prove unacceptable. If it should be required to introduce gaps into the lighting pattern this may be achieved in several ways, one being by the definitio of a dummy light model having zero intensity; another method is that of modifying one of the existing lighting models by a zero intensity factor. The former method may be preferred since this is more readily identified in the subsequent output data. The major visibility conditions w ich dictate the range at which a given light can be seen, are the illuminance threshold of the pilot's eyes and the effective mean extinction coefficient of the atmosphere between the pilot and the ground. In the present version of the program the illuminance threshold is single-valued, whereas in practice this is a variable dependent not only on the quantity of fog between the pilot and the ground in the direction of view, but also on the direction of view in terms of the veiling illumination due to the sun by day and to the ground lighting in the pilot's near field by night. OMPI The parameter which represents the attenuation of the atmosphere • is that of the mean extinction coefficient which is derived from the definitive σ, profile. The path of the pilot's eyes is defined in terms of height above the ground and lateral displacement from the centreline. Both the models are of the segmented-line type as mentioned above and described & Appendix A. In addition, limits for the excursion in longitudinal position are given. These enable the pilot's motion to be restricted to that portion of the approach or landing which gives information relevant to the case in hand. As a matter of convenience these models are organised with the elements commencing at the upwind and running end back towards the approach, so that in effect the pilot flies backwards up the approach path. When it is required to run cases involving offset viewing or curved flight paths sign conventions assume some importance. In order to facilitate the construction and interpretation of cases all spatial coordinates are referred to ground axes. The choice of the origin is somewhat arbitrary, however, the use of the intersection of the runway centreline with either the runway threshold or glidepath origin will generally be the most convenient. The convention adopted is that the downwind end of the pattern, a lateral displacement to the right (when viewed from the downwind end), and the pilot eye height are all positive. When considering lamp setting angles a change to local axis is effected (for convenience of implementation). The convention is such that azimuth angles are referred to a line parallel to the centreline and positive to the right when looking downwind. Whilst the elevation angle is considered positive, it should be noted that the lamp is set in azimuth before the elevation angle is applied. The described sign conventions follow the approach normally adopted by engineers in considering airfield lighting. As has been stated above, a number of models must be defined and variables specified before the program can be run. In order to minimise the work necessary to assemble the data the models are organised in the following order: lamp, models, lighting pattern, illuminance threshold, fog model, pilot eye path models and finally program controls. Thus, for instance, changes can be made to the fog model without having to re-enter the lamp and pattern models. This is achieved by the provision of a number of entry points in the progra so that for a succeeding case only the data for the models requiring change and those following in the list need to be entered. The 'constructions and format of input data for a single run is shown in Appendix B. This may be compared with actual input data, an example of which is shown in Appendix C. It should be noted that apart from the caption, only the numbers in Appendix C have any significance for the computer and that the headings merely serve to break up the number into more recognisable groups. Examination of Appendix C will reveal the numbers to be of both REAL and INTEGER typr. The variables have been defined by type in Appendix B, since with the exception of L . ™ in the lighting pattern, integers are unsubscripted. A further point of interest is the use of the variable R in the lighting pattern, which serves to close all the parameter blocks in the lighting pattern on their upwind side. It will be seen that both Appendix B and Appendix C include a section labelled CASE, a description of which is below: If it is assumed that input data has been read then the general method of determining the position of the limit of vision is as follows: The upwind limit of X which is the pilot's initial longitudinal position is used to determine the lateral position and height from the eye path models. The height is then used to determine the mean extinction coefficient to the ground by means of the fog model. For this initial position of the pilot a search is made along the ground from the upwind end of the pattern and parallel to the centreline at 30 m intervals. At each such point on the ground the appropriate lig is positioned and given the required setting angles. The cartesian coordinates of the pilot are then calculated with respect to the angled and tilted lamp axes, so that the pilot's bearings may be determined as azimuth and elevation angles relative to the principal axes of the lamp. By this means the effective intensit -t O li- W directed towards the pilot can be calculated using the relevant lamp model. This value is then factored as required and substituted together with the slant range and the mean extinction coefficient into the equation for Allard's Law. The resulting eye illuminance is then tested against the stated illuminance threshold. Once a ground position has been found where the eye illumination exceeds the threshold value, ie a lamp is visible, then the step size is reduced and the direction of the search reversed. This procedure is repeated until two successive points are within 0.2 m. The ground position is now moved 0.3 m in the direction of increased eye illumination and the various parameters calculated as before, these values being the ones used as output data. The pilot's eye is now moved to a new position 10 m further downwind and the procedure is repeated. Slight modification of the procedure occurs when visual contact was achieved for the pilot's previous position (which is the normal situation) in that the initial search step is reduced from 30 m to 10 m and the search commences 100 m upwind of the previous ground position. This allows a generally better resolution of the visual sequence at initial contact without any increase in the computation time. The present version of the program gives a solution of Allard's law to within 0.3 in after approximately 20 iterations. This occupies a computation time of • approximately two thirds of a second. It should be noted that the ground positions determined do not, in general, relate to the actual position of lamps, rather where such a lamp should be placed for threshold eye detection. However, the computed results do allow the actual positions of lamps to be considered retrospectively. Measures are taken to combat a number of anomalous situations that may arise during the search procedure. Such conditions include the visual range being dictated by the extent of the lighting pattern rather than by the visibility, and a complete absence of visual contact because the pilot is either too high, or too far downwind for the prevailing conditions, or upwind of all lights. In addition, an alternative mode Of operation permits the deter¬ mination of the nearest visible light. Once visual contact has been established as has been described the search is continued with the OMPI ,, WIPO » discriminants reversed so that threshold eye illumination is sought for a nearer lamp. In this way a more, realistic picture can be bu up of the extent to which the lighting pattern is visible since bot the far and near points of vision ar known. To further elucidate the operation of the program a generalise flow chart is shown in Fig 6. This figure outlines the complete process from the reading of the input data to the termination of th program. In order to prevent confusion some labelled entry points have been substituted in place of flow paths. It should be pointe out that the name of the logical variable 'GO' has no intrinsic meaning, but the logical variable 'USE* is employed to indicate the suitability of data for output purposes. Explanations of the code letters used in some of the boxes may be desirable since particular meaning has been given to the words in the key. Read is used to denote the reading of the appropriate input data, while find means that a value is read directly from the stored model. This latter contrasts with the use of determine which implies that some calcula tion is required in order to obtain the required value from the mod Finally calculate alludes to the evaluation of algebraic expression not directly linked with the input models. Facilities exist for both the tabulation and plotting of certa of the variables. However, to avoid possible confusion, each case produces a tabulation of all the defining input models and paramete (Appendix D), irrespective of the current form of output employed. The tabulated output data includes all basic parameters and, i required, the current states of all models at the completion of eac search (Fig 7)» Briefly, the interpretation of the columns in Fig is as follows: 1. Height of the pilot's eyes above the ground plane (metres) 2. X component of separation between the lamp and the pilot (metres). 3- X position of the lamp ground axes (metres). k. X position of the pilot _OM 5. Lateral separation (Y) between the pilot and the line of lights (metres). in lamp axes (degrees). 8. Angle resulting from (6) and (7) uses to determine intensity (degrees). 9. Value of intensity from lamp model by application of (8) (candles). 10. Intensity expressed as a percentage of the peak value of that lamp. 11. Mean extinction coefficient along the vertical path between the pilot and the ground (km ). The remaining columns indicate which defining blocks and elements in Appendix C or Appendix D apply to the particular solution, whilst the column marked TEST uses a zero to indicate the satisfactory solution of Allard's law. In addition output is suppressed for those pilot positions when no visual contact was achieved. Finally, the right hand column shows the number of iterations required to achieve the resolution of 0.5 o. The examples shewn in Appendix D and Fig 7 and derived in part from Appendix C apply to a simple case where both the line of lights and the pilot are placed on the centreline. As a consequence no results have been obtained for the azimuth offset angle (col 7 in Fig 7 )• However, persons skilled in the art will, in view of the above, readily perceive how they can be obtained. APPENDIX A A number of models are required within the program and a standard technique has been applied. Parameters which may require . approximations to be made to continuous functions are the polar distribution of light about the principal axis of a lamp, the vertical gradient of fog density, and the pilot's eye height and lateral position. An approach to modelling which is both easy to comprehend and implement, is that of the segmented straight line approximation. Since data sources for the models will inevitably be diverse in nature ranging perhaps from a few discrete points to highly accurate continuous functions it is suggested that the segments are fitted by eye to graphical plots of the source data. It should be noted that when calculating slopes of segments the should be defined by the coordinates of intersection even if these a only estimated, rather than estimates of the slopes themselves, ie: ~ 3 - - - ~ -ή This prevents the occurrence of disruption of the models at the intersection of the segments. For the purposes of input data to the program the segments have been defined in terms of the lower value of X at which the segment becomes effective, the slop directed towards increasing X, and the value of Y appropriate to the lower value of X. In the case of both lamp and fog models the initial value of X will be zero, whilst for lamps the value of both the slope and Y for the final segment will normally be zero. Failure to observe this practice in cases of normal light distributions which fall off with angle from the principal axis may give rise to the anamolous condition of negative intensity. Similarly in defining vertical fo density which decreases with height care must be taken to prevent th occurrence of negative for (clearer than cleari). During a recent study of visual sequences it was found that the lamps were adequately described by five segments, however the facili for a variable number of segments has been retained for both fog and pilot's eye path models. APPENDIX B The various models and program controls are read in the followin order: ENTRY POINTS 1. LAMPS (light distribution 2. LIGHTING PATTERN (configuration) 3. ET (visual threshold) ' 4. FOG (vertical model) 5. PILOT (eye path) 6. CASE (output and program control). The form of the data required to establish the above models is set out below: (ENTRY POINT 1) LAMPS (refer Appendix A) NL number of different lamp models to be read (the first model read is subsequently known as lamp 1) R,. ratio of AZ EL axes in isocandle diagram. (First lamp) is the angle at the break point in degrees M is the slopeof the segment in candela/degree I is the intensity at the break point in candela NB: five elements describing one half of the vertical light distribution. (NL th lamp) NL θ 1NL = °-° M 1NL " 1NL 3NL M. 3NL 3NL 5NL M '5.NL X 5NL =°'° (ENTRY POINT 2) LIGHTING PATTERN LP Lateral position of line of lights NT Total number of blocks of all lamp types. >1 NE Total number of blocks of all lamp elevation settings. NA Total number of blocks of all lamp pzimuth settings. NI Total number of blocks of all lamp _intensity factors. >1 The following positions are those of the downwind end of blocks. POT LT. 1 NT position pf block 1 for lamp _type lamp type in block Λ_ tt it P0T NT LT NT position pf block NT for lamp type lamp type in block NT P0E„ E. 1 1 position pf block 1 for plevation lamp plevation in block _1. NE NE position pf block NE for elevation lamp plevation in block N P0A 1 1 position pf block 1 for azimuth lamp azimuth in block 1 It II 11 II P0A„„ LA NA ^NA position pf block NA for azimuth lamp azimuth in block NA P0I 1 * U 1 position pf block 1 for intensity _1 lamp _intensity factor in factor block 1_ II II II II OMP /,, WWIIPP POI LI NI NI position pf block NI for intensity lam-p ^Intensity factor in factor block NI R p upwind end of last block in the line of lights. (ENTRY POINT k) FOG MODEL NUMBER specifying type of model in use may assume values of 1 , 2 or 3« 1 Continuous statistical model of mean extinction coefficient. 2 Segmented function of mean extinction coefficient. 3 Segmented function of extinction coefficient. HN number of elements in segmented fog models 2. Note: No illustration of the use of segmented fog models is given, if such models are required the technique described in Appendix A should be followed. For the continuous model, HN should assume a dummy value of 1. The vertical profile is specified by two parameters S„ and C, where S is the extinction coefficient (km ) at a height of 2 m and C is a constant related to the probability of occurrence. % This is used to obtain the mean extinction coefficient S by means of an equation within the program as set out below. where the height H is given in metres. (ENTRY POINT 5) PILOT NH number of elements in height position model * 2 NL number of elements in ^ lateral position model > 2 OMPI /,, WIPO « 22 X or XX is range at the break point in metres S or SS is slope of segment in metres/metre NH ~NH Z and Y are the height and lateral displacement X. SSm at the break point in metres. J W upwind and d_ownwind limits of the pilot's eye pa "V Note: The upwind limit should be downwind of both the extreme values (ie X^ and X ^-) of the eye path models. (ENTRY POINT 6) CASE HTU a three digit number controlling the form of the output. Hundreds Program mode control 0 far point computed 1 near point computed Tens Tabulation output control >1 basic parameter tabulation on line printer >2 additional program progress and state indicators Units Digital plotter control <5 generate new axes for plotter >5 superimpose next plot on existing axes 0 plotter output suppressed 1 or 6 height/visual ground range 2 or 7 height/elevation offset angle 3 or 8 height/azimuth offset angle • + or 9 (not allocated) RUN a serial or date number up to 6 digits (appears on all tabulations and plots). INSTRING caption for all tabulations (and plots when new axes are generated), up to 120 characters within apostrophies. NC number dictating entry point for subsequent pase. 1 LAMPS π 2 LIGHTING 3 ET (entry at a given point requires that and 4 FOG all following models to be re-read). 5 PILOT 6 CASE N number exceeding 6 terminates program. In the above, code letters have been used to represent variables. These have been chosen to act as shorthand descriptions of the variables and users assembling input data may find their use as annotations eases the task of maintaining the correct data format. APPENDIX C LAMPS 2 2.0 0.0 -1000.0 7500.0 0.5 -2333.3 7000.0 2.6 -615.4 2100.0 5.2 -100.0 500.0 10.2 0.0 0.0 2.4 0.0 -1666.7 30000.0 0.6 -5239.4 29000.0 4.15 -3125.0 10400.0 6.55 -1054.5 2900.0 9.3 0.0 0.0 OMPI /,, WIPO <ι LIGHTING PATTERN 0.0 2 5 1 2 900.0 2 0.0 1 900.0 6.0 650.0 5.5 450.0 5.0 300.0 4.5 0.0 3.0 900.0 0.0 900.0 1 .0 0.0 0.33 -I5OO.O ET - .G FOG 1 1 8.76 1 .50 PILOT 6 . 2 -1100.0 0.0 8.9 -6O5.O 0.01808 8.9 - 45.O 0.02263 13.6 -I 5.O 0.03053 17.9 -6O.O 0.04188 20.8 +100.0 0.05241 .27-5 -1500.0 0.0 0.0 0.0 0.0 0.0 OMP Λ. WIP LIMITS ■620.0 +300.0 CASE 021 100574 'B 47 3 DEG AUTOLAND ICAO LIGHTING DAY P OF FOG GRADIENT = 0.5 5/6/7*10 2 M RVR ' 6 106 2 •NEAR PT ' 4 FOG 1 1 6.69 1 .50 PILOT 6 2 -1100.0 0.0 8.9 -605.0 0.01808 8.9 -345.0 0.02263 13.9 -155.0 0.03053 17.9 -60.0 0.04188 20.8 +100.0 0.05241 27.5 ■1500.0 0.0 0.0 0.0 0.0 0.0 LIMITS -62O.O ^400.0 CASE 006 3 26 'FAR PT» 6 106 4 'NEAR FT ' (The remainder has been omitted for brevity) APPENDIX D RUN IDENT 100574 B 473 DEG AUTOLAND ICAO LIGHTING DAY P OF FOG GRADIENT = 0.5 -5/β/7*ΛO 2 M RVR LAMP 1 AZIMUTH FACTOR 2.00 Angle Slope Intensity 0.00 -1000.0 7500.0 0.50 -2333-3 7000.0 2.60 -615.4 2100.0 5.20 -100.0 500.0 10.20 0.0 0.0 LAMP 2 AZIMUTH FACTOR 2.4θ Angle Slope Intensity 0.00 -1666.7 30000.0 0.60 -5239.4 29OOO.O 4.15 -3125.O 10400.0 6.55 -IO54. 2900.0 9.30 0.0 0.0 LATERAL POSITION 0.0 LIGHTING PATTERN Position Angle Type 900.0 2 0.0 1 ELEVATION 900.0 6.00 650.0 5.50 450.0 5.00 300.0 4.50 0.0 3.00 AZIMUTH 900.0 0.00 INTENSITY FACTOR 900.0 1.00 0.0 0.33 RUNWAY STOP END -1500.0 LUMINANCE THRESHOLD -4.00 FOG EXTINCTION COEFFICIENT Height Slope Sigma 8.76 1.50 PILOT EYE PATH A/C Position Slope Height -1100.0 0.0000 8.9 -605.O 0.0181 8.9 -345.O 0.0226 13.6 -15 .O 0.0305 17.9 -6O.O 0.0419 20.8 100.0 0.0524 27.5 A/C Position Slope Lateral -I5OO.O 0.0000 0.0 0.0 0.0000 0.0 Computation Limits Upwind -620 Downwind 300 Output Option 21 In shallow fog conditions a tall aircraft override facility 41 gives information as to the minimum height a pilot must be from the ground when the aircraft is landed for uninterrupted visual acquisition of runway lighting. Figures 4 and 5 illustrate the variation of slant extinction coefficient with height in actual shallow and deep fog situations together with model profiles fitted to accord with the three transmissometer measurements. As shown in figures 4a and 5a, the model and the empirical profiles deviate one from another by no mor than +5%, a correlation which has been found consistent throughout all deep and shallow fog conditions examined by the inventor and which conditions are by far the more common in any fog timescale. Thus is validated his discovery that a reliable estimate can be obtained from just three slant readings at different heights. _OMPl

Claims: CLAIMS Claim 1. Apparatus for estimating a fog structure comprising a transmissometer array disposed for slant operation at different heights above the ground, and means for employing signals derived from the transmissometer array in estimating slant visual range. Claim 2. Apparatus as claimed in claim 1 and wherein the means for employing signals derived from the transmissometer array in estimating the slant visual range comprises means for converting the transmissometer signals into signals representing the mean integral of extinction coefficient (°\ ) over the associated distance, means for h comparing the &, signals and selecting an appropriate model from a predetermined set of models, means for determining the values of the constants in the selected model, and thence the definitive models, and means for deriving a mean integral of extinction coefficient profile from the definitive models. Claim 3« Apparatus as claimed in claim 2 and wherein the means for converting thhee ttrr<ansmissometer array signals into o" signals employs the function: where R is the distance over which the measurement is taken, I, is the signal from the transmissometer, and I is a clear day datum transmissometer signal. Claim 4. Apparatus as claimed in claim 2 or claim 3 and wherein the predetermined set of models comprises: = ah + bh + c for use when cr > σ, > or h 3 h 2 h 1 and (h 2 or h ) + k for use when <τ, < cr, > <τ, , and qh, + r h 3 h 2 h ι and k when cr, < σ, < - . q(h or h ) + r h 3 h 2 Claim 5« Apparatus as claimed in any one of claims 2 to 4 and having means for determining slant visual range and visual sequence from the cr profile. Claim 6. Apparatus as claimed in any one of claims 2 to 5 and having means for determining visual contact height comprising means for determining the intersection of the function: Visual contact height (H and the definitive model. Claim 7« Apparatus as claimed in claim 6 and having means for deriving Log_ |n (E_) from the function: Log 1Q (E τ = u log 1Q β - v Claim 8. Apparatus as claimed in claim 7 and arranged to feed Log 1 _(E r _ [1 ) > - w into the H function. Claim 9« Apparatus as claimed in any one of claims 6 -to 8 and having means for storing aircraft configuration information, applyin this information to the derived H value, and modifying H according c c Claim 10. Apparatus as claimed in any one of claims 2 to 5 and wherein the means for determining slant visual range is iterative, comprising means for processing the derived definitive cr model, the particular lighting patterns, direction and intensity and visual threshold values at a sequence of aircraft positions, during a landing, run. Claim 11. Apparatus as claimed in any one of the preceding claims and wherein the transmissometer array is disposed for operati at three different heights above the ground. Claim 12. Apparatus as claimed in any one of the preceding claims and wherein the transmissometer array comprises a source/ detector transmissometer unit mounted at or close to ground level and a plurality of reflectors mounted at the different heights. OMP /., WIP &RM

Inventors: PUFFETT A
Applicants: PUFFETT A, US GOVERNMENT
Application year: 1980
Main classification: G01W1
IPC codes: B64C13, B64F1, G01N21, G01W1, G08G5