# Near-Field/Far-Field Transformation with Helicoidal Scanning from Irregularly Spaced Data

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

---

## Abstract

A fast and accurate technique for the compensation of the probe positioning errors in the near-field/far-field transformation with helicoidal scanning is proposed in this paper. It relies on a nonredundant sampling representation using a spherical modelling of the antenna under test and employs an iterative scheme to evaluate the near-field data at the points fixed by the helicoidal nonredundant representation from the acquired irregularly distributed ones. Once these helicoidal data have been recovered, those required by a classical cylindrical near-field/far-field transformation are efficiently determined by using an optimal sampling interpolation algorithm. Some numerical tests assessing the effectiveness of the proposed approach and its stability with respect to random errors affecting the near-field data are shown.


---

## Body

## 1. Introduction

As well-known, far-field (FF) range size limitations, transportation, and mounting problems can make impossible or impractical the measurement of antenna radiation patterns on a conventional FF range. In these cases, it is convenient to exploit near-field (NF) measurements and recover the FF patterns by using NF-FF transformation techniques [1–3]. In addition, the NF measurements may be performed in a controlled environment, as an anechoic chamber, thus overcoming those drawbacks that cannot be eliminated in FF outdoor measurements. In this framework, the reduction of the time needed for acquiring the NF data is assuming an ever growing relevance for the antenna measurement community, since this time is currently very much greater than that required to perform the corresponding NF-FF transformation. Such a reduction can be achieved by decreasing the number of the NF data to be collected and/or by making faster the acquisition of each NF value. A significant reduction of the number of required NF data has been obtained for all the conventional scannings [4–8] by applying the nonredundant sampling representations of electromagnetic (EM) fields and the optimal sampling interpolation (OSI) expansions [9], whereas, the use of the modulated scattering technique employing arrays of scattering probes has been proposed in [10] to realize a very fast electronic scanning. However, antenna testing NF facilities based on such a technique are not very flexible as those employing mechanical scans. These last can be made faster by exploiting innovative spiral scannings which use, as suggested by  Yaccarino et al. in [11], continuous and synchronized movements of the positioning systems of the probe and antenna under test (AUT). In particular, accurate, stable, and efficient NF-FF transformations using the helicoidal scanning, the planar and spherical spiral scannings have been developed [12–19]. They are based on the aforementioned nonredundant representations and reconstruct the NF data needed by the NF-FF transformation with the corresponding classical scanning by interpolating, via appropriate OSI formulas, the nonredundant samples acquired on the spiral. The required two-dimensional algorithm has been obtained (a) by assuming the AUT enclosed in a proper convex domain bounded by a surface Σ with rotational symmetry; (b) by developing a nonredundant sampling representation of the voltage acquired by the probe on the spiral; (c) by choosing the spiral step equal to the spacing needed to interpolate the data along a meridian curve. In particular, the AUT has been assumed to be enclosed in the smallest sphere able to contain it in [12–15], whereas more effective AUT modellings, that allow a further reduction of required NF data when dealing with antennas having one or two predominant dimensions, have been adopted in [16–18] by properly applying the unified theory of spiral scannings for nonspherical antennas [19]. This last has been obtained by heuristically extending the rigorous approach for spherical antennas [15]. It is worth noting that these effective modellings allow one to consider measurement cylinders (planes) with a radius (distance) smaller than one half the antenna maximum size, thus reducing the error related to the truncation of the scanning zone.Unfortunately, it may be impossible to get regularly distributed NF data due to an inaccurate control of the positioning systems, but the measurements points position can be accurately read by optical devices. In addition, the finite resolution of the positioning systems and their imprecise synchronization do not allow one to exactly locate the probe at the points fixed by the sampling representation. In the light of the above considerations, the development of an accurate and stable reconstruction procedure from irregularly spaced data becomes relevant. In this context, an approach based on the conjugate gradient iteration method and using the unequally spaced fast Fourier transform [20, 21] has been proposed in the planar [22] and spherical [23] classical scannings. However, such an approach is not tailored for scannings exploiting the nonredundant sampling representations of EM fields, wherein the “a priori” information on the AUT and proper sampling interpolations are employed to recover the NF data required by the corresponding standard NF-FF transformation technique. The interpolation from nonuniform samples has been well investigated in the one-dimensional case of bandlimited functions defined over the real axis. Several sufficient conditions assuring the possibility of reconstructing a function from its nonuniform samples have been stated in [24]. The stability in a nonuniform sampling algorithm; that is, the requirement that small errors affecting the samples give rise to small errors in the reconstructed functions, has been exhaustively investigated in [25, 26], wherein it has been shown that a stable sampling cannot be accomplished at an “average” rate lower than the Nyquist one. Some closed-form expressions for the interpolation of bandlimited functions from nonuniform samples have been developed in [27]. However, they are valid only for particular sampling points arrangements, are cumbersome and not user friendly. Moreover, they become more and more unstable as the sampling points distribution deviates from the uniform one. The two-dimensional nonuniform sampling has not attracted an equal consideration. In any case, it has been shown [28] that, again, a stable sampling cannot be performed at a rate lower than the Nyquist one. As it has been clearly stressed in [29], wherein a more exhaustive discussion on this topic can be found, a nonuniform sampling algorithm useful for practical applications must be computationally manageable, accurate, and stable. Accordingly, it is more convenient to recover the uniform samples from the irregularly distributed ones than to resort to a direct interpolation formula. In fact, once the uniform samples have been determined, the value at any point of the scanning surface can be recovered by an accurate and stable OSI expansion.In this context, two different approaches have been proposed [29–31]. The former [29, 30] is based on an iterative technique which has been found convergent only if there exists a one-to-one correspondence associating at each uniform sampling point the nearest nonuniform one. The latter [31] exploits the singular value decomposition (SVD) method [32] and has been applied when the two-dimensional problem can be reduced to the research of the solution of two independent one-dimensional ones. This, f.i., occurs in a cylindrical near-field facility, wherein the nonuniformly distributed data can be realistically assumed to lie on not regularly spaced rings when the measurements are made by rings [31]. Such a hypothesis is no longer valid in the helicoidal scanning. Accordingly, the iterative technique will be here applied for reconstructing the uniformly spaced helicoidal samples from the acquired irregularly distributed data. Obviously, the SVD-based approach could be generalized to such an intrinsically two-dimensional problem, but the dimension of the involved matrix would become too large, thus requiring a huge computational effort.In the following, the helicoidal scanning based on the spherical AUT modelling will be considered since it is more simple and general than those using modellings tailored for elongated antennas.


## 2. Nonredundant Sampling Representation on a Cylinder

Let us consider an AUT and a nondirective probe scanning a helix with constant angular step lying on a cylinder of radiusd (Figure 1) and adopt the spherical coordinate system (r,ϑ,φ) to denote an observation point P in the NF region. Since, as shown in [33], the voltage measured by such a kind of probe has the same effective spatial bandwidth of the AUT field, the nonredundant sampling representations of EM fields [9] can be applied to it. Accordingly, by assuming the AUT as enclosed in the smallest sphere of radius a able to contain it and describing the helix by a proper analytical parameterization r̲=r̲(ξ), it is possible to consider the “reduced voltage”
(1)Ṽ(ξ)=V(ξ)ejγ(ξ),
where γ(ξ) is an optimal phase function to be determined. The error, occurring when Ṽ is approximated by a spatially bandlimited function, becomes negligible as the bandwidth exceeds a critical value Wξ [9], so that it can be effectively controlled by choosing a bandwidth equal to χ′Wξ, where χ′ is an enlargement factor, slightly greater than unity for electrically large antennas.Figure 1
Geometry of the problem.The parametric equation of the helix, when imposing its passage through a fixed pointP0 of the generatrix at φ=0, are(2)x=dcos(ϕ-ϕs),y=dsin(ϕ-ϕs),z=dcotθ,
where ϕ is the angular parameter describing the helix, ϕs is the value of ϕ at P0, and θ=kϕ. Such a helix can be constructed as intersection of the cylinder with the line from the origin to a point moving on a spiral which wraps the sphere enclosing the AUT. In order to allow the two-dimensional interpolation, the helix step Δθ (Figure 1) must be equal to the spacing required for the voltage interpolation along a generatrix [6]. Therefore, the parameter k is such that the angular step, determined by the consecutive intersections of the helix with a generatrix, is Δθ=2π/(2N′′+1), with N′′=Int(χN′)+1 and N′=Int(χ′βa)+1. Accordingly, being Δθ=2πk, it results that k=1/(2N′′+1). The function Int(x) denotes the integer part of x,β is the wavenumber, and χ>1 is an oversampling factor needed for controlling the truncation error [9].A nonredundant sampling representation of the voltage on the helix can be obtained by using the following expressions for the phase function and parameterization [15](3)γ=β∫0r1-a2r′2dr′=βr2-a2-βacos-1(ar),ξ=βaWξ∫0ϕk2+sin2kϕ′dϕ′.As can be seen, the optimal parameterξ is proportional to the curvilinear abscissa along the spiral wrapping the sphere modelling the AUT. 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 AUT sphere is described. Therefore,(4)Wξ=βaπ∫0(2N′′+1)πk2+sin2kϕ′dϕ′.In the light of these results, the reduced voltage at any point of the helix can be reconstructed via the OSI expansion:(5)Ṽ(ξ)=∑m=m0-p+1m0+pṼ(ξm)ΩM(ξ-ξm)DM′′(ξ-ξm),
where m0=Int[(ξ-ξs)/Δξ] is the index of the sample nearest (on the left) to the output point, 2p is the number of retained samples Ṽ(ξm), and (6)ξm=ξ(ϕs)+mΔξ=ξs+2πm(2M′′+1)
with M′′=Int(χM′)+1, M′=Int(χ′Wξ)+1, and M=M′′-M′. Moreover,(7)DM′′(ξ)=sin((2M′′+1)ξ/2)(2M′′+1)sin(ξ/2),(8)ΩM(ξ)=TM[2(cos(ξ/2)/cos(ξ̅/2))2-1]TM[2/cos2(ξ̅/2)-1]
are the Dirichlet and Tschebyscheff Sampling functions, respectively, TM(·) being the Tschebyscheff polynomial of degree M and ξ̅=pΔξ.The OSI expansion (5) can be used to evaluate the “intermediate samples”; that is, the reduced voltage values at the intersection points between the helix and the generatrix passing through P. Once these samples have been evaluated, the reduced voltage at P can be recovered via the following OSI formula:(9)Ṽ(ϑ,φ)=∑n=n0-q+1n0+qṼ(ϑn)DN′′(ϑ-ϑn)ΩN(ϑ-ϑn),
where n0=Int[(ϑ-ϑ0)/Δϑ], N=N′′-N′, Ṽ(ϑn) are the intermediate samples (10)ϑn=ϑn(φ)=ϑ(ϕs)+kφ+nΔϑ=ϑ0+nΔϑ
and the other symbols have the same meaning as in (5). It is so possible to recover the NF data required to perform the standard NF-FF transformation with cylindrical scanning [34], whose key steps are reported in the next section for reader’s convenience.


## 3. Probe Compensated NF-FF Transformation

As rigorously demonstrated in [34], the modal coefficients av and bv of the cylindrical wave expansion of the field radiated by the AUT are related to (a) the two-dimensional Fourier transforms Iν and Iν′ 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); (b) 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. In particular,(11)aν(η)=β2Λ2Δν(η)×[Iν(η)∑m=-∞∞dm′(-η)Hν+m(2)(Λd)-Iν′(η)∑m=-∞∞dm(-η)Hν+m(2)(Λd)],bν(η)=β2Λ2Δν(η)×[Iν′(η)∑m=-∞∞cm(-η)Hν+m(2)(Λd)-Iν(η)∑m=-∞∞cm′(-η)Hν+m(2)(Λd)],Iν(η)=∫-∞∞∫-ππV(φ,z)e-jνφejηzdφdz,Iν′(η)=∫-∞∞∫-ππV′(φ,z)e-jνφejηzdφdz,Δν(η)=∑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 V, V′ 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,Θ,Φ) can be evaluated by(12)EΘ=-j2βe-jβRRsinΘ∑ν=-∞∞jνbν(βcosΘ)ejνΦ,EΦ=-2βe-jβRRsinΘ∑ν=-∞∞jνaν(βcosΘ)ejνΦ.


## 4. Uniform Samples Reconstruction

Let us now suppose that the samples are irregularly distributed (Figure1) and denote with (ϑ̅i,φ̅i) the nonuniform sampling point corresponding to the nearest uniform one ξi lying on the helix. By expressing the reduced voltage at each nonuniform sampling point as a function of the unknown values at the nearest uniform ones via the following two-dimensional OSI formula (obtained by properly merging expansions (5) and (9)), it results in(13)Ṽ(ϑ̅i,φ̅i)=∑n=n0-q+1n0+q[∑m=m0-p+1m0+pDN′′(ϑ̅i-ϑn)ΩN(ϑ̅i-ϑn)·∑m=m0-p+1m0+pṼ(ξm)ΩM(ξ(ϑn)-ξm)DM′′(ξ(ϑn)-ξm)∑m=m0-p+1m0+p(ϑ̅i-ϑn)]i=1,2,…,Q,
where Q is the number of samples. The system (13) can be recast in the matrix form(14)A̲̲x̲=b̲,
where A̲̲ is a Q×Q sparse banded matrix whose elements are given by(15)Aim=DN′′(ϑ̅i-ϑn)ΩN(ϑ̅i-ϑn)×DM′′(ξ(ϑn)-ξm)ΩM(ξ(ϑn)-ξm),
where b̲ is the vector of the collected nonuniform data, and x̲ is the vector of the unknown uniform samples.By splittingA̲̲ in its diagonal and nondiagonal parts, A̲̲D and Δ̲̲, respectively, multiplying both members of (14) by A̲̲D-1 and rearranging the terms, it results in(16)x̲=A̲̲D-1b̲-A̲̲D-1Δ̲̲x̲.The following iterative scheme is so obtained:(17)x̲(ν)=A̲̲D-1b̲-A̲̲D-1Δ̲̲x̲(ν-1)=x̲(0)-A̲̲D-1Δ̲̲x̲(ν-1),
where x̲(ν) is the vector of the uniform samples estimated at the νth step. Necessary conditions for the convergence of such a scheme [29, 30] are that Aii≠0,foralli, and |Aii|≥|Aim|,forallm≠i. These conditions are surely verified in the assumed hypothesis of biunique correspondence between each uniform sampling point and the “nearest” nonuniform one.By straightforward evaluations from (17), it results
(18)Ṽ(ν)(ξi)=1Aii{∑n=n0-q+1n0+q[ΩN(ϑ̅i-ϑn)DN′′(ϑ̅i-ϑn)Ṽ(ϑ̅i,ϕ̅i)-∑n=n0-q+1n0+q[∑m=m0-p+1m0+pΩN(ϑ̅i-ϑn)DN′′(ϑ̅i-ϑn)·∑m=m0-p+1m≠im0+pṼ(ξm)(ν-1)ΩM(ξ(ϑn)-ξm)×DM′′(ξ(ϑn)-ξm)∑m=m0-p+1m0+p⇳()\\\(ϑ̅i-ϑn)]∑1D45B=n0-q+1n0+q}i=1,2,…,Q,
where(19)Aii=DN′′(ϑ̅i-ϑni)ΩN(ϑ̅i-ϑni)×DM′′(ξ(ϑn)-ξm)ΩM(ξ(ϑn)-ξm)ni=Nint[(ϑ̅i-ϑ0)/Δϑ] being the index of the intermediate sampling point nearest to the uniform one ξi.


## 5. Numerical Results

The effectiveness and robustness of the proposed algorithm for compensating the probe positioning errors in the NF-FF transformation with helicoidal scanning have been assessed by many numerical tests. The reported simulations refer to a uniform circular array lying on the planey = 0 (see Figure 1), symmetric with respect to the plane z=0 and having radius a=16λ, λ being the wavelength. Its elements are elementary Huygens sources polarized along the z axis and are radially and azimuthally spaced of λ/2. An open-ended WR-90 rectangular waveguide, operating at the frequency of 10 GHz, is chosen as probe. The considered helix wraps a cylinder with radius d=20λ and height 2h=140λ. The irregularly distributed samples have been generated by imposing that the distances in ξ and ϑ between the position of each nonuniform sample and the associate uniform one are random variables uniformly distributed in (-0.3Δξ,0.3Δξ) and (-0.3Δϑ,0.3Δϑ). Note that this represents a pessimistic occurrence in a real scanning procedure.Figures2 and 3 show a representative reconstruction example of the output voltage V (the most significant one) on the generatrix at φ=90°, obtained by 0 and 8 iterations, respectively. As it can be seen, only 8 iterations are enough to get a very good reconstruction. The evaluation of the maximum and mean-square errors (normalized to the maximum value of the output voltage V on the cylinder) in the reconstruction of the uniform samples assesses more quantitatively the effectiveness of the proposed algorithm. They have been obtained by comparing the reconstructed uniform samples and the exact ones. As can be seen (Figures 4 and 5), on increasing the number of iterations, the errors decrease quickly until a constant saturation value is reached. Such a value decreases on increasing the retained samples number.Even better results are to be expected when the distances between the nonuniform samples and the uniform ones are smaller. Moreover, Figure 6 shows the normalized maximum and mean-square errors in the reconstruction of the NF data needed to carry out the NF-FF transformation [34] both when using the directly collected nonuniform samples and the recovered uniform helicoidal ones. As can be seen, in this last case an increased accuracy of about 60 dB can be obtained.Figure 2
Amplitude of the probe voltageV on the generatrix at φ=90°. Solid line: exact. Crosses: recovered from irregularly spaced NF data at the iteration 0.Figure 3
Amplitude of the probe voltageV on the generatrix at φ=90°. Solid line: exact. Crosses: recovered from irregularly spaced NF data at the iteration 8.Figure 4
Normalized maximum error in the reconstruction of the uniform helicoidal samples.Figure 5
Normalized mean-square error in the reconstruction of the uniform helicoidal samples.Figure 6
Normalized maximum (■) and mean-square  (●) errors in the reconstruction of the NF data needed for the NF-FF transformation. Dotted lines: when using the nonuniform samples. Solid lines: when using the recovered helicoidal uniform samples.The robustness of the algorithm has been assessed (see Figure7) by corrupting the exact samples with random errors. Both a background noise (bounded to Δa in amplitude and with arbitrary phase) and uncertainties on the data of ±Δar in amplitude and ±Δα in phase have been simulated.Figure 7
Amplitude of the probe voltageV on the generatrix at φ=90°. Solid line: exact. Crosses: recovered at the iteration 8 from error affected irregularly spaced NF data.At last, Figures8 and 9 show the antenna FF pattern in the principal planes E and H reconstructed from the irregularly distributed helicoidal samples. As can be seen, the exact and recovered patterns are practically indistinguishable, thus providing an overall assessment of the proposed iterative technique.Figure 8
E-plane pattern. Solid line: exact. Crosses: recovered from the irregularly spaced helicoidal NF data.Figure 9
H-plane pattern. Solid line: exact. Crosses: recovered from the irregularly spaced helicoidal NF data.It is worth noting that the number of employed samples (guard samples included) for reconstructing the NF data on the considered cylinder is 26 817, significantly less than that (71 680) required by the standard cylindrical scanning and by the helicoidal scanning technique [35].



---

*Source: 859396-2010-09-26.xml*