ITU-R REPORT P 2090-2006 Measuring the input parameters for the radiative energy transfer model of vegetation attenuation《植被衰减辐射能量转移模型的输入值测量》.pdf

1、 Rep. ITU-R P.2090 1 REPORT ITU-R P.2090 Measuring the input parameters for the radiative energy transfer model of vegetation attenuation (2006) 1 Introduction The radiative energy transfer (RET) theory allows an analytical calculation of the attenuation of a scattering medium, such as vegetation. R

2、ecommendation ITU-R P.833-5 details the RET-based model and gives the necessary input parameters for several tree species. However, if the vegetation is different to that found in Recommendation ITU-R P.833-5, one may derive the input parameters for RET from measurements. This Report details how to

3、derive the RET parameters from a measurement. 2 RET input parameters The RET requires six input parameters, four of which are specific for the scattering medium and two are geometry and system descriptors. The four medium dependent parameters need to be established experimentally and are dependent o

4、n the type of vegetation, frequency and state of foliation. i) , the ratio of the forward scattered power to the total scattered power; ii) , the beamwidth of the phase function; iii) , the combined absorption and scatter cross-section (a+ S); iv) W, the albedo; v) R, the beamwidth of the receiving

5、antenna; vi) d, the distance into the vegetation (m). Parameters and are best determined by measuring the scatter function usually referred to as the “phase” function of the vegetation medium. The phase function characterizes the scattering behaviour of the medium at short distances after the air/ve

6、getation interface. When the four RET parameters are extracted from the vegetation attenuation data only, the resulting parameters may be lacking in real physical properties of vegetation scattering in most of cases, even though the parameters could fit the measurement data well. So, it is strongly

7、recommended that the phase function measurement must be performed whenever possible. Paragraph 3 describes the phase function measurements and the method of determining the parameters and in more detail. Parameters and W can be determined from attenuation curve measurements into the vegetation mediu

8、m as described in 4. 3 Phase function measurements The phase function forms an integral part in the determination of the parameters needed for the RET theory. This Report presents a guide outlining how to conduct phase function measurements. The main purpose of such measurements will be to extract p

9、arameters and required for the derivation of the RET model. The method of determining parameters and from measurement data is also described in this Report. 2 Rep. ITU-R P.2090 3.1 Choosing the measurement site Phase function measurements are necessary to derive parameters and for the vegetation med

10、ium, which in turn enables RET modelling of excess attenuation in the medium. In many cases the phase function measurement will therefore be conducted in conjunction with measurements of the excess attenuation experienced in the vegetation medium as a function of vegetation depth. An ideal site will

11、 comprise of a group of trees, behind which receiver access at various depth locations is possible. The site requires a reasonable foreground clearance to enable the tree group to be illuminated by the transmitter, thus fulfilling the condition of a half space of vegetation, illuminated with an inci

12、dent plane wave signal. Furthermore, the site needs to be free of objects that can cause spurious reflected signals to interfere with the scattered signal from the trees to the receiver location. Phase function measurements combined with excess attenuation measurements, performed for a line of trees

13、 made up of the same tree species can jointly provide the RET parameters. The results can be used to develop a prediction model for the behaviour of excess attenuation as a function of vegetation depth. 3.2 Measurement geometry The RET describes the energy transfer in a medium exhibiting both scatte

14、ring and absorbing properties. This is treated by the RET as homogenous infinite half space. Both transmit and receive antennas are arranged to have a common volume irradiated in the scattering medium. Under these conditions the RET postulates a phase function with Gaussian shaped forward lobe and a

15、n isotropic background. The two parameters and , define the exact shape of the phase function, which for the vegetation medium is dependent on factors such as the species and density of vegetation and foliation state. The location at which the phase function is measured is also important. The RET pr

16、edicts that the forward lobe broadens as the observation point moves inside the medium. This occurs due to the attenuation of the coherent component and the relative increase in the level of the isotropic scatter. This suggests that there is an optimum region within which the measurement should be c

17、onducted. The following conditions are proposed for the measurement geometry: i) The region in which the phase function measurement is carried out needs to lie within that defined by the initial predominately linear (dB scale) slope of the excess attenuation curve. A location two to three trees deep

18、 into the medium is recommended. This allows sufficient amount of scattering medium after the interface to be present to result in an isotropic back-scatter coherent forward lobe, without the widening of the forward lobe yet taking place. ii) The assumption of homogeneity of the medium in the RET is

19、 more easily satisfied in the canopy region of the tree rather than that containing the trunk. A further benefit of largely confining the common volume to the canopy region accrues from the avoidance of strong diffracted and reflected signal components emanating from trunks and the ground. In single

20、 tree measurements of the phase function it was found that illuminating around 70% to 90% of the tree canopy yields good results. When groups of trees are considered, the transmit antenna vegetation distance should be fixed carefully. This distance, in combination with the transmit antenna elevation

21、 beamwidth, is utilized to ensure that ground reflections and diffraction components over the top of the canopy are minimized. Rep. ITU-R P.2090 3 iii) The phase function forward lobe is aligned with the direction of propagation of the incident signal. It is therefore preferable to place the transmi

22、t and receive antennas at the same height. This also ensures normal incidence conditions. If the heights differ the phase function will be tilted with respect to the horizontal and the direction of incidence. The orientation of the receiver antenna as well as that of the axis of rotation needs to be

23、 adjusted accordingly. The beamwidth of the phase function is relatively narrow, so that tilts of a few degrees with respect to the receiver alignment will lead to significant reduction of the magnitude of the forward lobe and consequently the determination of . 3.3 Antenna characteristics 3.3.1 Tra

24、nsmit antenna The beamwidth of the transmit antenna needs to be wide enough to illuminate sufficient width and height of vegetation in order to provide, together with the receive antenna, sufficient excited common volume for the measurement. However the beamwidth must be, relatively speaking, narrow

25、 so as not to give rise to diffracted and reflected signal components, which will influence the measured phase function pattern. 3.3.2 Receive antenna The receive antenna ideally has a very small beamwidth, in order to measure the relatively narrow beamwidth of the phase function itself. In practice

26、 a receive antenna of a few degrees beamwidth gives good results. 3.4 Receiver antenna displacement The most practical and most accurate method to measure the phase function is to place the receiver antenna inside the vegetation medium at a distance from the air vegetation interface specified in 3.2

27、. At this location the receiver antenna needs to be rotated about its vertical axis in discrete increments; after each increment the signal level needs to be recorded. Since the beamwidth of the phase function is considered to be relatively narrow the antenna only needs to be rotated over a section

28、of a circle, e.g. 40 with respect to the line receiver transmitter. The principal set-up is sketched in Fig. 1. FIGURE 1 Principal measurement set up for phase function measurement 4 Rep. ITU-R P.2090 3.5 Link budget considerations Apart from the free space loss due to the distance between transmitt

29、er and receiver the excess attenuation of the vegetation as well as the difference between forward lobe and isotropic back-scatter region of the phase function have to be considered. At millimetre and near millimetre frequencies the attenuation resulting from one large dense tree can be considerable

30、, i.e. in the region of a few tens of dB. Even at 11 GHz a large horse chestnut tree in full leaf was found to give as much as 35 dB attenuation. The difference between forward lobe and back-scatter regions of the phase function, is expected to be between 20 and 50 dB depending on the values of . Th

31、e measurement system must have sufficient dynamic range to overcome the strong initial excess attenuation due to vegetation, and the difference in level between the forward lobe signal and the back-scatter. As a guide the dynamic range of the receiver needs to be about 80 dB to yield the phase funct

32、ion shape and parameters. 3.6 Determination of parameters and Parameters and can be determined from the measurement data via best-fit curve fitting. In order to extract parameters and , the following model of phase function needs to be calculated: )1()()( += fp (1a) where: 2)/(22)(= ef (1b) Figure 2

33、 illustrates the fitting of predicted phase function curve to measurement data. 4 Excess attenuation curve measurements Excess attenuation curve measurements are obtained by measuring the signal level at various depths into the vegetation medium. The RET predicts a graph with an initially linear par

34、t (dB) and a transition region into a final attenuation rate. The slope of the initial slope is used to obtain parameter . Whereas the transition region and the final attenuation rate are used to find parameter W (albedo) via best-fit curve fitting. In the case when phase function data are not avail

35、able, parameters W, and must be determined by curve fitting. 4.1 Determination of parameter using initial slope The first term Iriof the RET equation describes the initial slope of the attenuation function with distance z (m). zePP=0(2) where: P: measured signal power inside the vegetation 0: power

36、received at the incident air-vegetation interface. Rep. ITU-R P.2090 5 The RET describes excess attenuation due to vegetation, so that the signal decay, due to free space loss, is not included in P. FIGURE 2 Predicted and measured phase function vs. rotation angle If the attenuation curve is display

37、ed in dB, equation (2) can be written as: )loglog(10log100ezPP= (3) The initial slope of the graph (dB/m) is given by: =43.4zA(4) where A is the attenuation (dB). 4.2 Determination of parameter W (albedo) by curve fitting Figure 3 illustrates the fitting of predicted attenuation curves with differen

38、t albedos to measured excess attenuation. 6 Rep. ITU-R P.2090 FIGURE 3 Predicted and measured excess attenuation vs. optical density for different albedos In order to model the RET attenuation curves the RET equation needs to be calculated. The RET equation is as follows: )(1112)()(!14)(2212112maxIs

39、eAPeIqqWmeqeeIePPkNNNkskNRMmMmmMRriRk+=+=(5) where: PR: received power by the receiving antenna with its gain pattern max: received signal strength received in the absence of vegetation R : beamwidth of the receiving antenna m: order of the term I1. Rep. ITU-R P.2090 7 The term I1is more accurately

40、evaluated for higher values of m, however it will not change significantly for m 10. The following relations also apply: zsa+= )( defines the optical density and z is the distance (m). )2(sin2NPN= (6) = )1( W The line in the RHS of equation (5) represents Iri, which is usually referred to as the “fi

41、rst” term, whereas the sum of the second and third line give I1, referred to as the “second” term. These appear in explicit form. The fourth line represents the “third” term I2of the RET equation, which contains the attenuation coefficients skand the amplitude factors Ak. These need to be determined

42、 numerically using the following expressions: The attenuation coefficients skare determined by the characteristic equation: =NnnnsPW0112(7) where: )1.,1(,sinsin = = NnNnNPnWWWW1)1(= is referred to as the reduced albedo. The amplitude factors Akare determined by a system of linear equations given by:

43、 NnNNkknkPsA =+=211for NNn 21+= where: =Nnncos 0=nfor Nn and: 1=nfor Nn = (8) The antenna gain beamwidth and the scatter function beamwidth in equation (5) have been assumed to be Gaussian, in which the main lobe of the power radiation pattern is reduced by 1/e. This angle is related to the 3 dB bea

44、mwidth by db36.0 =R. 8 Rep. ITU-R P.2090 4.3 Determination of parameter W, and by curve fitting, when phase function data are not available When phase function data are not available, three parameters W (albedo), and must be determined from the measured excess attenuation data via best-fit curve fitting of equation (5). In this case, the extracted parameters may lack physical meaning as mentioned before.