# An Innovative Direct NF-FF Transformation Technique with Helicoidal Scanning

Authors: Francesco D’Agostino; Flaminio Ferrara; Claudio Gennarelli; Rocco Guerriero; Massimo Migliozzi
Journal: International Journal of Antennas and Propagation
 (2012)

---

## Abstract

A direct near-field-far-field transformation with helicoidal scanning is developed. It is based on the nonredundant sampling representation of electromagnetic fields and uses a spherical antenna modelling to determine the number of helix turns. Moreover, the number of voltage samples on each of them is fixed by the maximum transverse dimension of the antenna, both to simplify the mechanical scanning and to reduce the computational effort. This technique allows the evaluation of the antenna far field directly from a minimum set of near-field data without interpolating them. Although the number of near-field data employed by the developed technique is slightly increased with respect to that required by rigorously applying the nonredundant sampling representation on the helix, it is still remarkably smaller than that needed by the standard near-field-far-field transformation with cylindrical scanning. The effectiveness of the technique is assessed by numerical and experimental results.


---

## Body

## 1. Introduction

The techniques for the reconstruction of antenna radiation patterns from near-field (NF) measurements have been widely investigated and used for applications ranging from cellular phone antennas to large-phased arrays and complex multibeam communication satellite antennas [1–5]. They have been proved to be efficient and attractive alternatives to conventional far-field (FF) and compact range measurements. In addition, the NF measurements have the advantage to be performed in a controlled environment, such as an anechoic chamber, thus overcoming those drawbacks due to weather conditions, electromagnetic (EM) interference, and so forth, which cannot be eliminated in FF outdoor measurements. In the recent years, the researchers interested in the antenna NF measurements and in the related NF-FF transformation techniques have spent many efforts to reduce the time required for the acquisition of the NF data. As a matter of fact, this time is currently very much greater than that needed to carry out the corresponding NF-FF transformation. In this context, a significant reduction of the number of required NF data (and, as a consequence, of the measurement time) has been obtained for all the conventional scannings [6–14] by applying the theoretical results on the nonredundant sampling representations of EM fields [15]. In particular, optimal sampling interpolation (OSI) formulas of central type [16] have been employed to efficiently recover the data required by the corresponding NF-FF transformation from the knowledge of the acquired nonredundant ones.A more convenient way of reducing the measurement time is the use of innovative spiral scanning techniques. They have been implemented, as suggested by Yaccarino et al. in [17], by means of continuous and synchronized movements of the positioning systems of the probe and antenna under test (AUT), and are based on the aforementioned nonredundant sampling representations and OSI expansions. In particular, NF-FF transformations using the helicoidal scanning [18–20], the planar [20, 21] and spherical [20, 22] spiral scannings have been developed, by considering the AUT as enclosed in the smallest sphere able to contain it. Then, more effective AUT modellings, that allow a further reduction of the required NF data in the case of elongated or quasiplanar antennas, have been adopted in [23–27] by properly employing the unified theory of spiral scans for nonspherical antennas [28]. These modellings allow one to consider measurement cylinders (planes) with a radius (distance) smaller than one half the AUT maximum size, thus reducing the error related to the truncation of the scanning zone. In all the cases, the NF data needed by the corresponding NF-FF transformation can be reconstructed by interpolating the nonredundant ones acquired on the spiral.A probe uncompensated NF-FF transformation technique with planar spiral scanning has been proposed in [29], to efficiently evaluate the antenna far field directly from the acquired NF data. It uses the convolution property of the radiation integral and the fast Fourier transform (FFT) algorithm to determine the antenna far-field pattern without requiring any interpolation step. However, since such an approach does not exploit the nonredundant representations of EM fields, it needs a useless large amount of measurements. This comment holds also for the NF-FF transformation technique with helicoidal scanning proposed by the same authors in [30]. On the contrary, a direct NF-FF transformation with cylindrical scanning, based on the aforementioned nonredundant sampling representations, has been developed in [10]. It allows the evaluation of the antenna far field in any cut plane directly from the nonredundant NF data without interpolating them.In this paper, the direct NF-FF transformation with cylindrical scanning [10] has been extended to the helicoidal one [18–20]. To this end, the approach in [10] has been reviewed in order to match the advantages of the direct cylindrical NF-FF transformation with those own of the fast helicoidal scanning. In particular, the number of helix turns is determined by the nonredundant sampling representation along a generatrix which makes use of the spherical AUT modelling, whereas the voltage samples on each of them is fixed by the AUT maximum transverse dimension.


## 2. Nonredundant Sampling Representations on a Cylinder from Helicoidal Samples

The nonredundant sampling representation of the probe voltage on a cylinder from NF data acquired along a helix, when adopting a spherical AUT modelling, is summarized in the following.Let us consider a nondirective probe which scans a helix with constant angular step lying on a cylinder of radiusd surrounding the AUT (see Figure 1) and adopt the spherical coordinate system (r,ϑ,φ) for denoting the observation point P in the NF region. Since the voltage V measured by this kind of probe has the same effective spatial bandwidth of the field [7], the theoretical results on the nonredundant representation of EM fields [15] can be applied to such a voltage. Accordingly, if the AUT is enclosed in a sphere of radius a (AUT ball) and the helix is described by a proper analytical parameterizationr̲=r̲(ξ), the probe “reduced voltage”(1)Ṽ(ξ)=V(ξ)ejγ(ξ)
can be closely approximated by a spatially band limited function [15], γ(ξ) being a phase function to be determined. The related band limitation error becomes negligible as the bandwidth exceeds a critical value Wξ [15], so that it can be effectively controlled by choosing a bandwidth equal to χ′Wξ, wherein the excess bandwidth factor χ′ is slightly greater than unity for an electrically large AUT.Figure 1
Helicoidal scanning.To obtain a sampling representation of the voltage on the cylinder from its nonredundant samples collected along a helix with constant angular step, it is necessary (a) to develop a nonredundant voltage representation on the helix; (b) to choose the helix step equal to the sample spacing needed to interpolate the voltage along a cylinder generatrix.The parametric equations of the helix, when imposing its passage through a fixed pointQ0 of the generatrix at φ=0, are x=dcos(ϕ-ϕi), y=dsin(ϕ-ϕi), z=dcotθ, where ϕ is the angular parameter describing the helix, ϕi is the value of ϕ at Q0, and θ=kϕ. The parameter k is such that the angular step, determined by the consecutive intersections Q(ϕ) and Q(ϕ+2π) with a generatrix, is Δθ=2πk. It is worth noting that the helix can be obtained by radially projecting on the measurement cylinder the spiral wrapping the AUT sphere with the same angular step.As shown in [19, 20], a nonredundant sampling representation of the voltage on the helix can be obtained by using the following expressions for the optimal phase function and parameterization:(2)γ=β∫0r1-a2r′2dr′=βr2-a2-βacos⁡⁡-1(ar),(3)ξ=βaWξ∫0ϕk2+sin⁡2kϕ′dϕ′,
where β is the wavenumber.Since the elevation step of the helix must be equal to the sample spacing required for the interpolation along a generatrix, then, according to [19, 20], Δθ=2π/(2N′′+1), with N′′=Int(χN′)+1 and N′=Int(χ′βa)+1, χ>1 being an oversampling factor controlling the truncation error [15, 16] and Int(x) denoting the integer part of x. It can be easily verified that k=1/(2N′′+1).According to (3), ξ is proportional to the curvilinear abscissa along the spiral wrapping the sphere modelling the source. Since such a spiral is a closed curve, it is convenient to choose the bandwidth Wξsuch that ξ covers a 2π range when the whole curve on the sphere is described. As a consequence,(4)Wξ=βaπ∫0(2N′′+1)πk2+sin⁡⁡2kϕ′dϕ′.According to these results, the OSI formula of central type to reconstruct the voltage at any pointQ of the helix is [20]:(5)Ṽ(ξ)=∑m=m0-p+1m0+pṼ(ξm)ΩM(ξ-ξm)DM′′(ξ-ξm),
where m0=Int[(ξ-ξ(ϕi))/Δξ] is the index of the sample nearest (on the left) to Q, 2p is the number of retained samples Ṽ(ξm), and(6)ξm=ξ(ϕi)+mΔξ=ξ(ϕi)+2πm(2M′′+1)
with M′′=Int(χM′)+1 andM′=Int(χ′Wξ)+1. Moreover,(7)DM′′(ξ)=sin⁡((2M′′+1)ξ/2)(2M′′+1)sin⁡⁡(ξ/2),ΩM(ξ)=TM[-1+2(cos⁡(ξ/2)/cos⁡⁡(ξ̅/2))2]TM[-1+2/cos⁡⁡2(ξ̅/2)]
are the Dirichlet and Tschebyscheff sampling functions, wherein TM(ξ) is the Tschebyscheff polynomial of degree M=M′′-M′ and ξ̅=pΔξ.The OSI formula (5) can be used to evaluate the “intermediate samples,” namely, the voltage values at the intersection points between the helix and the generatrix passing through the observation point P. Once these samples have been evaluated, the voltage at P can be reconstructed via the following OSI expansion:(8)Ṽ(ϑ,φ)=∑n=n0-q+1n0+qṼ(ϑn)ΩN(ϑ-ϑn)DN′′(ϑ-ϑn),
where N=N′′-N′, n0=Int[(ϑ-ϑ0)/Δϑ],2q is the number of the retained intermediate samples Ṽ(ϑn), and(9)ϑn=ϑn(φ)=ϑ(ϕi)+kφ+nΔϑ=ϑ0+nΔϑ.The described two-dimensional OSI algorithm can be properly applied to recover the NF data required by the NF-FF transformation technique [31] or [32].


## 3. Standard NF-FF Transformation with Cylindrical Scanning

For reader’s convenience, the key steps of the standard probe compensated NF-FF transformation technique with cylindrical scanning [31] are reported. According to such a technique, the modal coefficients aν and bν of the cylindrical wave expansion of the field radiated by the AUT are related to (i) the two-dimensional Fourier transforms Iν1 and Iν2 of the output voltage of the probe for two independent sets of measurements (the probe is rotated 90° about its longitudinal axis in the second set); (ii) the coefficients cm,dm, and cm′,dm′ of the cylindrical wave expansion of the field radiated by the probe and the rotated probe, respectively, when used as transmitting antennas. The key relations are(10)aν(η)=β2Λ2Δν(η)[Iν1(η)∑m=-∞∞dm′(-η)Hν+m(2)(Λd)-Iν2(η)∑m=-∞∞dm(-η)Hν+m(2)(Λd)],bν(η)=β2Λ2Δν(η)[Iν2(η)∑m=-∞∞cm(-η)Hν+m(2)(Λd)-Iν1(η)∑m=-∞∞cm′(-η)Hν+m(2)(Λd)],(11)Iν1,2(η)=∫-∞∞∫02πV1,2(φ,z)e-jνφejηzdφdz,(12)Δv(η)=∑m=-∞∞cm(-η)Hν+m(2)(Λd)×∑m=-∞∞dm′(-η)Hν+m(2)(Λd)-∑m=-∞∞cm′(-η)Hν+m(2)(Λd)×∑m=-∞∞dm(-η)Hν+m(2)(Λd),
where Λ=β2-η2, Hν(2)(·) is the Hankel function of second kind and order ν, and V1, V2 represent the output voltage of the probe and the rotated probe at the point of cylindrical coordinates (d,φ,z).Once the modal coefficients have been determined, the FF components of the electric field in the spherical coordinate system(R,Θ,Φ), adopted to denote an observation point in the far-field region, can be evaluated by(13)EΘ(R,Θ,Φ)=FΘ(Θ,Φ)e-jβRR=-j2βe-jβRRsin⁡⁡Θ∑ν=-∞∞jνbν(βcos⁡⁡Θ)ejνΦ,EΦ(R,Θ,Φ)=FΦ(Θ,Φ)e-jβRR=-2βe-jβRRsin⁡⁡Θ∑ν=-∞∞jνaν(βcos⁡⁡Θ)ejνΦ.


## 4. Direct NF-FF Transformation Technique

The here proposed approach for evaluating the antenna FF pattern directly from the acquired helicoidal samples is described in the following. In particular, such a NF-FF transformation, unlike that developed in [18, 20], does not explicitly employ the OSI expansions (5) and (8), since there is no need to interpolate the nonredundant acquired data for reconstructing those needed by the classical ones [31, 32].The starting point is the evaluation of the two-dimensional Fourier transforms of the probe output voltage in (11). By taking into account the OSI formula (8), we can rearrange the integrals (11) in the form:(14)Iν1,2(η)=∑n∈Nr∫02π{∫-∞∞DN"(ϑ(z)-ϑn)Q(ϑ(z)-ϑn)e-jγ(z)ejηzdzṼ1,2(ϑn,φ)e-jνφ×∫-∞∞DN′′(ϑ(z)-ϑn)Q(ϑ(z)-ϑn)e-jγ(z)ejηzdz}dφ,
where Nr is the set of indexes of all considered NF turns, Ṽ1,2(ϑn,φ) are the reduced voltages samples at the intersection points between the generatrix at φ and the scanning helix, and Q=ΩN, if|ϑ(z)-ϑn|≤qΔϑ, or Q = 0, otherwise. Let us first consider the integration over z. For any fixed φ, it results
(15)Gnη(φ)=∫zizfDN′′(ϑ(z)-ϑn)ΩN(ϑ(z)-ϑn)e-jγ(z)ejηzdz,
where zi=z(ϑn+qΔϑ) and zf=z(ϑn-qΔϑ). By taking into account (15), the relation (14) can be rewritten as follows:(16)Iν1,2(η)=∑n∈Nr∫02πṼ1,2(ϑn,φ)Gnη(φ)e-jνφdφ.Such a relation involves an integration overφ that must be solved numerically for any n∈Nr, by using, from time to time, the voltage samples on the considered helix turn, whose number varies, according to the nonredundant sampling representation on the helix (see Section 2), when moving from the central to the peripheral turns and does not ensure the convergence of the numerical integration. Moreover, the Gnη(φ) values depend not only on the measurement cylinder and on the radius of the sphere enclosing the AUT, but also on the points at which the generatrix passing through φ intersects the scanning helix. It can be easily recognized that integration over φ can be efficiently solved by using the FFT algorithm, provided that the number of the voltage samples on each helix turn be always the same and equal to the smallest integer MH, product of powers of 2, 3, and 5 equal or greater than 2[Int(χ′βρmax)+1], ρmax being the AUT maximum transverse dimension. In such a way, the samples lying on the helix at φm=mΔφ=2πm/MH with m=0,…,MH-1 are all aligned. Accordingly, we get(17)∫02πṼ1,2(ϑn,φ)Gnη(φ)e-jνφdφ=2πMH∑m=0MH-1Ṽ1,2(ϑm,n,φm)Gnη(φm)e-j2πmν/MH,
where(18)ϑm,n=ϑn(φm)=ϑ(ϕi)+kφm+nΔϑ.The summation in (17) can be, obviously, efficiently performed via a direct FFT algorithm.It is worthy to note that theGnη(φm) values can be calculated (once and for all) for given sets of antennas.According to the above results, the determination of the far field is obtained through the following steps.(a)
The samples of the probe and rotated probe voltages acquired at the constant stepΔφ along the scanning helix specified in Section 2 are multiplied for the phase factor ejγ.(b)
For each required value of the polar angleΘ, fixing the corresponding value of η by the relation η=βcosΘ, the Gnη(φm) values are computed or read, if already calculated.(c)
For each helix turn, the FFT of the sequenceṼ1,2(ϑm,n,φm)Gnη(φm) is performed.(d)
The probe voltage Fourier transformsIν1,2(η) are evaluated by performing the summation over n∈Nr. The corresponding values of the cylindrical modal expansion coefficients aν and bν can be then evaluated.(e)
The evaluation of the far field at the considered angleΘ is finally effectively achieved by performing the summations in (13) again with the FFT.From the efficiency viewpoint, it is convenient to use this method to evaluate only the FF samples necessary to recover the antenna pattern via the following far-field OSI expansion tailored for even numbers of samples along the meridians and parallels [10]:(19)FΘ,Φ(Θ,Φ)=2NF′′-12NF′′∑n=n0-q+1n0+q{∑m=m0-p+1m0+pΩNF(Θ-Θn)DNF′′-1(Θ-Θn)2Mn′′-12Mn′′×∑m=m0-p+1m0+pFΘ,Φ(Θn,Φm,n)×ΩMn(Φ-Φm,n)×DMn′′-1(Φ-Φm,n)∑m=m0-p+1m0+p},
wherein n0=Int[Θ/ΔΘ] and m0=Int[Φ/ΔΦn] are the indexes of the sample nearest (on the left) to the output point, and(20)Θn=nΔΘ=nπNF′′,Φm,n=mΔΦn=mπMn′′,(21)NF′′=2[Int(χN′2)+1],NF=NF′′-N′,(22)Mn′′=2i≥Int(χMn′)+1,Mn′=Int[χ*βρmax⁡⁡sin⁡⁡Θn]+1,Mn=Mn′′-Mn′,χ*=1+(χ′-1)[sin⁡⁡Θn]-2/3.The need of an OSI expansion tailored for an even number of samples along the parallels is due to the employment of an efficient power of two FFT algorithms for computing (13), whereas NF′′ has been chosen according to (21) in order to have FF samples on the equator.It is worthy to note that there is no need to extract the phase factor from the far-field expression, since it is constant on the far-field sphere.


## 5. Numerical and Experimental Results

Many numerical tests have been performed to assess the validity of the presented technique. In particular, three sets of figures are reported.The first one refers to a uniform planar circular array (see Figure1) symmetric with respect to the plane z=0 and having radius a=16λ,λ being the wavelength. Its elements are elementary Huygens sources linearly polarized along the z-axis, radially and azimuthally spaced of 0.6λ. An open-ended WR-90 rectangular waveguide, operating at the frequency of 10 GHz, is considered as probe. The radius d of the cylinder wrapped by the scanning helix is 25λ, and its height h is 150λ.Figures2 and 3 show the FF pattern reconstruction in the principal planes obtained by using the direct helicoidal NF-FF transformation described in Section 4. As can be seen, the exact and recovered fields are practically indistinguishable, thus assessing the effectiveness of the technique. Moreover, a further example showing the reconstructed FF pattern in the plane at Φ=60° is reported in Figure 4, in order to show the capability of the technique to accurately reconstruct the far field in any cut-plane.Figure 2
Far-field pattern in theE-plane. Solid line: exact field. Crosses: reconstructed via the direct helicoidal NF-FF transformation.Figure 3
Far-field pattern in theH-plane. Solid line: exact field. Crosses: reconstructed via the direct helicoidal NF-FF transformation.Figure 4
Far-field pattern in the plane atΦ=60°. Solid line: exact field. Crosses: reconstructed via the direct helicoidal NF-FF transformation.The second set of figures is relevant to a smaller antenna. It is again a uniform planar circular array of elementary Huygens sources linearly polarized along thez-axis. Its radius a is now 4.2λ, and its elements are radially and azimuthally spaced of 0.4λ. The scanning helix wraps a cylinder with height h equal to 80λ and radius d=15λ. In order to highlight another peculiar feature of the here proposed direct helicoidal NF-FF transformation, the FF patterns in the E-plane obtained by applying such a NF-FF transformation and the standard cylindrical probe-compensated one [31] are compared with the exact FF pattern in Figures 5 and 6, respectively. As can be seen, the reconstruction obtained via the proposed transformation is clearly more accurate than that achieved by the standard one, since the characteristic ripple caused by the discontinuity of the near field at the edge of the scanning surface is not present. This is due to the different method employed to evaluate the two-dimensional Fourier transforms Iν1 and Iν2 of the probe voltage. As a matter of fact, in the standard transformation technique [31], they are computed via FFT by taking into account only the NF data falling in the measurement area so that the integration over z results truncated to the cylinder height. When the evaluation of Iν1 and Iν2 is carried out according to the proposed technique, the effect of each NF sample is considered (see (15)) in the range [zi,zf], so that the peripheral samples can affect the evaluation even far from the scanning area. Thus, the proposed technique intrinsically eliminates the field discontinuity at the edges of the scanning surface, without using any windowing function, that, as well known, reduces the ripple as well as the zone of good far-field reconstruction.Figure 5
Far-field pattern in theE-plane. Black line: exact field. Red line: reconstructed via the direct helicoidal NF-FF transformation.Figure 6
Far-field pattern in theE-plane. Black line: exact field. Red line: reconstructed via the standard cylindrical NF-FF transformation [31].The third set of figures refers to the experimental validation of the direct helicoidal NF-FF transformation. Such a validation has been carried out in the anechoic chamber of the UNISA Antenna Characterization Lab. The chamber is equipped with a rotating table and a vertical scanner supplied by MI Technologies. The amplitude and phase measurements are performed by means of a vectorial network analyzer. An open-ended WR90 rectangular waveguide is used as probe. The considered antenna is a MI-12-8.2 standard gain horn with aperture19.4cm×14.4cm, located on the plane y=0 of the adopted reference system and operating at the frequency of 10 GHz. Such an AUT has been considered as enclosed in a spherical surface having radius equal to 12.6 cm. The probe output voltages have been collected along a helix lying on cylinder with d= 45 cm and h= 240 cm.It must be stressed that the proposed NF-FF transformation incorporates the probe characterization. Therefore, first of all, the employed probe has been characterized according to [33] as done in the software package MI-3000 implementing the standard probe compensated NF-FF transformation with cylindrical scanning [32]. Then, it has been verified that practically identical results are obtained when the same NF data are transformed by using the MI package and the so developed version of the probe compensated NF-FF transformation [31]. Afterwards, the FF pattern in the principal planes E and H reconstructed by using the direct helicoidal NF-FF transformation has been compared with that obtained via the MI software (reference 1). These reconstructions are shown in Figures 7 and 8. As can be seen, there is a good agreement save for the zones characterized by very low field levels, wherein the discrepancies are due to the fact that two different NF data sets have been employed. As expected, a better agreement results (see Figures 9 and 10) when the NF data used by the direct helicoidal NF-FF transformation are employed to obtain via a proper interpolation expansion those needed by the software package MI-3000, thus getting a new FF pattern (reference 2). The last reconstruction example (Figure 11) is relevant to the FF pattern in the plane at Φ = 60°.Figure 7
Far-field pattern in theE-plane. Black line: reference 1. Red line: reconstructed via the direct helicoidal NF-FF transformation.Figure 8
Far-field pattern in theH-plane. Black line: reference 1. Red line: reconstructed via the direct helicoidal NF-FF transformation.Figure 9
Far-field pattern in theE-plane. Black line: reference 2. Red line: reconstructed via the direct helicoidal NF-FF transformation.Figure 10
Far-field pattern in theH-plane. Black line: reference 2. Red line: reconstructed via the direct helicoidal NF-FF transformation.Figure 11
Far-field pattern in the plane atΦ=60°. Black line: reference 2. Red line: reconstructed via the direct helicoidal NF-FF transformation.As can be seen, the experimental results relevant to the FF reconstructions in theE-plane and in the plane at Φ=60° confirm the particular property of the developed technique to eliminate the ripple caused by the discontinuity of the near field at the edges of the scanning surface.It can be interesting to compare the number of NF data used by the direct NF-FF transformation (31 044 for the first example and 2 417 for the other two) with those (67 725 for the first example and 11 592 for the other two) needed by the MI software and by the helicoidal NF-FF transformation technique [30].For what concerns the time needed for the NF data acquisition, the technique in [30] is certainly quicker than the traditional one, since the acquisition is performed by continuous and synchronized movements of the positioning systems. Accordingly, it is enough to compare the here developed technique only with the one in [30]. In such a case, it can be easily realized that the measurement times are directly proportional to the numbers of needed samples.


## 6. Conclusions

A direct NF-FF transformation with helicoidal scanning, which allows one to evaluate the antenna far field in any cut plane directly from a minimum set of NF data without interpolating them, has been developed in this paper. To this end, the approach in [10] has been reviewed in order to match the advantages of the direct cylindrical NF-FF transformation with those own of the fast helicoidal scanning. In particular, the NF sampling arrangement is obtained by determining the number of helix turns according to the nonredundant sampling representation along a generatrix which makes use of the spherical AUT modelling and the voltage samples on each of them according to the AUT maximum transverse dimension. Many numerical tests have assessed the validity of the presented technique and its capability to eliminate the ripple caused by the discontinuity of the near field at the edges of the scanning surface. Moreover, the technique effectiveness has been further confirmed by the experimental results. At last, it has been shown that, although the number of NF data employed by the developed technique is slightly increased with respect to that required by rigorously applying the nonredundant sampling representation on the helix, it remains still remarkably smaller than that needed by the standard NF-FF transformation with cylindrical scan. It is worthy to note that the proposed technique, as well as all those using the spiral scannings, allow also the frequency extension of existing near-field ranges, since both measurement axes are used simultaneously hence reducing the effect of backlash.



---

*Source: 912948-2012-05-09.xml*