SAR image formation toolbox for MATLAB
|
|
- Kevin Moris Quinn
- 6 years ago
- Views:
Transcription
1 SAR image formation toolbox for MATLAB LeRoy A. Gorham and Linda J. Moore Air Fore Researh Laboratory, Sensors Diretorate 2241 Avionis Cirle, Bldg 62, WPAFB, OH ABSTRACT While many syntheti aperture radar (SAR) image formation tehniques exist, two of the most intuitive methods for implementation by SAR novies are the mathed filter and bakprojetion algorithms. The mathed filter and (non-optimized) bakprojetion algorithms are undeniably omputationally omplex. However, the bakprojetion algorithm may be suessfully employed for many SAR researh endeavors not involving onsiderably large data sets and not requiring time-ritial image formation. Exeution of both image reonstrution algorithms in MATLAB is expliitly addressed. In partiular, a manipulation of the bakprojetion imaging equations is supplied to show how ommon MATLAB funtions, ifft and interp1, may be used for straight-forward SAR image formation. In addition, limits for sene size and pixel spaing are derived to aid in the seletion of an appropriate imaging grid to avoid aliasing. Example SAR images generated though use of the bakprojetion algorithm are provided given four publily available SAR datasets. Finally, MATLAB ode for SAR image reonstrution using the mathed filter and bakprojetion algorithms is provided. Keywords: SAR, image formation, mathed filter, bakprojetion algorithm 1. INTRODUCTION Over the past four deades, numerous image formation algorithms have been developed with varying levels of omplexity and auray. However, these algorithms have a steep learning urve for SAR novies. In this paper, we introdue two algorithms - mathed filter and bakprojetion - implemented in MATLAB, whih are very simple to understand. A benefit of mathed filter image reonstrution is that the method offers optimization of the signal-to-noise ratio (SNR). 1 Use of the mathed filter requires a hypothesis regarding a target s harateristis; these assumptions yield a speifially designed filter that will ideally math to the target. While the mathed filter algorithm is intuitive and straight-forward, the generation of a 2D SAR image using this tehnique requires O(N 4 ) operations. For most pratial situations, the omputational omplexity of the mathed filter algorithm renders use of the method unrealisti. However, understanding of this fundamental imaging algorithm is extremely benefiial, espeially for the SAR novie. While several tehniques exist for SAR image generation, two ommonly employed approahes worth noting here are the polar format and bakprojetion algorithms. These algorithms differ in regards to the tradeoff between omputational omplexity and image quality. 2,3 The polar format algorithm 4,5 for SAR image reonstrution is similar to omputer-aided tomography imaging tehniques. 6,7,8 Polar format reonstrution employs a diret 2D Fourier transform, whih provides the opportunity to take advantage of the extremely effiient fast Fourier transform. Polar format image reonstrution requires O(N 2 log N) operations given an NxN SAR image. While more omputationally effetive than the bakprojetion algorithm, the polar format imaging method involves an approximation of the reeived signal phase. This approximation introdues errors into the final reonstruted SAR image. Several orretion tehniques have been suggested as a means of improving the auray of SAR images generated via polar format reonstrution. 9,1 Another frequently used image formation tehnique is the bakprojetion algorithm. 7,11 A distint advantage of the bakprojetion algorithm is the ability to form SAR images as phase history is olleted, pulse by pulse, and to integrate newly obtained information into the SAR image as it beomes available. The bakprojetion algorithm is omputationally expensive, however, as it requires O(N 3 ) operations for an NxN SAR image. Fortunately, the bakprojetion algorithm lends itself naturally to parallel proessing. Bakprojetion has reently Further author information: ( leroy.gorham@wpafb.af.mil and linda.moore2@wpafb.af.mil)
2 been implemented on graphis proessing units (GPUs), 12,13 whih provide inexpensive platforms for parallel omputing. Several methods have been proposed as a means of reduing the omputational omplexity of the bakprojetion algorithm while retaining the advantages of this image formation approah. 14,15,16,17 Many of these fast bakprojetion algorithms are flexible, providing manners in whih image quality may be aeptably sarified in order to further improve data proessing time. These image formation algorithms, several of whih are favorable given wide bandwidths and large syntheti apertures, offer dereased omputational omplexities of O(N 5/2 ) 14 and O(N 2 log N). 15,16,17 Non-optimized bakprojetion imaging may be impratial for exeptionally large data sets with time-ritial image formation requirements; however, knowledge of this simple SAR imaging tehnique is undoubtedly valuable for many researh endeavors. The remainder of this paper is outlined as follows. In Setion 2, we define our terminology and the SAR signal model. We also derive limits for sene size and pixel spaing to aid in seleting an imaging grid that avoids aliasing. In Setion 3, we derive and desribe the mathed filter algorithm for image formation. In Setion 4, we derive a bakprojetion algorithm, whih is based on the mathed filter algorithm. We desribe how to utilize basi MATLAB funtions like ifft and interp1 to implement the algorithm. In Setion 5, we apply the bakprojetion algorithm to four publily available SAR datasets. Finally, MATLAB ode for the two imaging algorithms, mathed filter and bakprojetion, are supplied in Appendix A.1 and A.2. The authors will be glad to provide opies of the MATLAB ode used to generate the figures provided in this paper in response to requests via SIGNAL MODEL The following SAR signal model 18 is used for the development of MATLAB implementations of the mathed filter and bakprojetion imaging algorithms. The SAR sensor travels along a flight path suh that the antenna phase enter has a three-dimensional spatial loation denoted r a (τ) suh that r a (τ) = [x a (τ), y a (τ), z a (τ)] T (1) where τ denotes the syntheti aperture, or slow time, domain. All oordinates are defined with respet to a sene origin, whih orresponds to the SAR motion ompensation point. The distane from the antenna phase enter to the origin, r a (τ), is denoted by A target is speified at a loation d a (τ) = x 2 a(τ) + y 2 a(τ) + z 2 a(τ). (2) r(τ) = [x(τ), y(τ), z(τ)] T. (3) In general, this target an have any arbitrary motion, but in this paper, we will assume that the target is stationary. Thus, the dependene on τ is dropped for the remainder of the derivation. The target has a radar ross setion whih is dependent on frequeny and aspet angle. The distane from the antenna phase enter to the target is denoted by d a (τ) = (x a (τ) x) 2 + (y a (τ) y) 2 + (z a (τ) z) 2. (4) At periodi intervals, the radar transmits a pulse that reflets off satterers in the sene. Some of the energy is then reeived by the radar. In a given syntheti aperture, there are N p pulses used to form the image. The transmission time of eah pulse is denoted by the sequene {τ n n = 1, 2,..., N p }. The output of the reeiver at a given time, τ n, is a sequene of band-limited frequeny samples delayed with respet to the time of pulse transmission by the round-trip time to the target. There are K frequeny samples per pulse, and the assoiated frequeny values are represented by the sequene {f k k = 1, 2,..., K}. The reeiver output from the target loated at r is j4πfk R(τ n ) S(f k, τ n ) = A(f k, τ n ) exp. (5)
3 The amplitude, A(f k, τ n ), is related to the radar ross setion of the target while the phase is dependent on the frequeny of eah sample and on the differential range, R(τ n ), given by R(τ n ) = d a (τ n ) d a (τ n ). (6) Equation (5) assumes that eah pulse is motion ompensated suh that a satterer at the sene origin will have zero phase for all f k and τ n. The atual reeiver output is thus the sum of the ontributions of all satterers in the sene. The frequeny samples, {f k }, have a minimum value denoted by f 1, a median value denoted by f, a maximum value denoted by f K, and a step size denoted by f. The frequeny step size is inversely related to the maximum alias-free range extent of the image, W r, by W r = 2 f. (7) Thus, the frequeny step size is hosen to math the size of the sene to be imaged. The total bandwidth, B, of the reeived pulse is B = (K 1) f. The range resolution is thus δ r = 2B = 2(K 1) f. (8) In a similar manner as above, the azimuth angle traversed during the syntheti aperture determines the ross-range resolution, and the azimuth angle from pulse to pulse determines the maximum alias-free ross-range extent of the image, W x. Given an azimuth step size of θ, W x = λ min 2 θ where λ min is the minimum wavelength suh that λ min = /f K. The total azimuth angle, θ a, traversed during the syntheti aperture is θ a = (N p 1) θ. Thus, the ross-range resolution, δ x, is given by (9) δ x = λ λ = 2θ a 2(N p 1) θ (1) where λ is the enter wavelength suh that λ = /f. To utilize the imaging algorithms desribed in this paper, one an selet any arbitrary pixel loations. However, one should be areful when seleting these loations to avoid aliasing of the image or the frequeny support. The pixel spaing should be finer than the resolution defined in Equations (8) and (1) and the overall sene size should be less than the maximum sene size defined in Equations (7) and (9). 3. MATCHED FILTER ALGORITHM The most straightforward method for forming a SAR image is to perform a mathed filter. One an build the mathed filter to any kind of satterer, but here we will assume an isotropi point satterer. The reeived signal from a point satterer at loation r is given in Equation (5). An isotropi satterer will have a onstant amplitude and thus A(f k, τ n ) = A. Therefore, the mathed filter response, denoted by I(r), at loation r is given by I(r) = 1 N p K assuming a single satterer in the sene. N p n=1 +j4πfk R(τ n ) S(f k, τ n ) exp = A, (11) To form an image using this method, Equation (11) is applied for eah pixel in the image. This requires alulation of the differential range, R(τ n ), for every pixel for every pulse. The algorithm has a omputational omplexity of O(N 4 ) for 2D images, whih makes it impratial for most appliations. However, Equation (11) forms the basis for the derivation of the bakprojetion algorithm in Setion 4. MATLAB ode for the mathed filter algorithm is provided in Appendix A.1.
4 MATLAB Type Symbol Desription Units data.deltaf salar f Frequeny step size Hz data.minf N p vetor f 1 Minimum frequeny for every pulse Hz data.nfft salar N fft Length of FFT (bakprojetion only) data.x mat Arbitrary r X oordinate of every pixel in the image m data.y mat Arbitrary r Y oordinate of every pixel in the image m data.z mat Arbitrary r Z oordinate of every pixel in the image m data.antx N p vetor r a Antenna loation for every pulse m data.anty N p vetor r a Antenna loation for every pulse m data.antz N p vetor r a Antenna loation for every pulse m data.r N p vetor d a (τ n ) Distane to the mo-omp point m data.phdata (KxN p ) array S(f k, τ n ) Complex phase history data, Eq. (5) Table 1. Neessary fields in the MATLAB data struture for image formation y (m) 25 y (m) x (m) (a) mathed filter image x (m) (b) bakprojetion image Figure 1. Comparison of SAR images (in db) of 3 point targets, loated at (,,), (-3,2,), and (1,4,), generated using the (a) mathed filter and (b) bakprojetion algorithms. 5 In order to use the MATLAB funtion, mfbasi, the inputs are provided in a struture alled data. The fields of this struture are defined in Table 1. The image pixel loations, r, are speified in the pixel loation matries data.x mat, data.y mat and data.z mat whih ontain the (x, y, z) loations of eah pixel. The output image is stored in data.im final, whih has the same dimensions as the pixel loation matries. While these pixels an have arbitrary positions, it is often desirable to form an image onto a regular 2D or 3D grid. In this ase, the MATLAB funtion meshgrid is useful to build the pixel loation matries. In all the image examples provided in this paper, meshgrid was used to build data.x mat and data.y mat, while data.z mat was filled with zeros. If a Digital Elevation Map (DEM) is available, data.z mat an be filled with interpolated height values for every pixel to improve image quality. Figure 1(a) shows an example of an image formed using the mathed filter algorithm. Phase history data was simulated for 3 point targets using Equation (5) where A(f k, τ n ) = 1. Here, N p = 128 pulses were simulated with K = 512 frequeny samples per pulse, a enter frequeny of 1 GHz and a 6 MHz bandwidth. A irular flight path was used with a 3 degree depression angle and a slant range of 1 km. A 3 degree integration angle was used with a enter azimuth angle of 5 degrees. The sene extent was 1 m x 1 m with 2 m pixel spaing in eah dimension.
5 4. BACKPROJECTION ALGORITHM The bakprojetion algorithm offers an intuitive imaging tehnique for SAR novies. While the bakprojetion algorithm is omputationally expensive, its implementation is not unreasonable for many SAR researh ventures. In this setion, bakprojetion imaging equations are manipulated in order to illustrate how the algorithm an be exeuted through use of ommon MATLAB funtions. 4.1 Derivation of Effiient Calulation of Range Profiles The mathed filter response shown in Equation (11) an be used to ompute the target response at a disrete range bin, m. Given SAR phase history, S(f k, τ n ), olleted by N p pulses over a range of K frequenies, the range profile at range bin m given a reeived pulse at slow time τ n is s(m, τ n ) = +j4πfk R(m, τ n ) S(f k, τ n ) exp. (12) By substituting the frequeny values f k = (k 1) f + f 1 into Equation (12), the range profile may be rewritten as +j4π((k 1) f + f1 ) R(m, τ n ) s(m, τ n ) = S(f k, τ n ) exp ( +j4π f R(m, τn )(k 1) = S(f k, τ n ) exp + +j4πf ) 1 R(m, τ n ) +j4πf1 R(m, τ n ) = S(f k, τ n ) exp[φ( R(m, τ n )) (k 1)] exp where phase funtion Φ( R(m, τ n )) = (+j4π f R(m, τ n ))/. In order to implement the bakprojetion algorithm in MATLAB, Equation (13) must be rewritten in terms of MATLAB s inverse disrete Fourier transform (ifft) funtion. In MATLAB, the definitions of the disrete Fourier transforms between X(k) and x(m) are given as 19 (13) X(k) = M m=1 x(m) = 1 K x(m) ω (m 1)(k 1) K X(k) ω (m 1)(k 1) K (14) where ω K = exp( j2π/k). Therefore, the inverse disrete Fourier transform of X(k) is given as x(m) = ifft(x(k)) = 1 K +j2π(m 1) X(k) exp (k 1) = 1 K K X(k) exp (Θ(m) (k 1)) (15) where Θ(m) = (j2π(m 1))/K. In order to write Equation (13) in terms of MATLAB s ifft funtion, the previously defined phase funtion, Φ( R(m, τ n )), from Equation(13), must equal Θ(m) from Equation (15). To satisfy this requirement, the following equality must hold: R(m, τ n ) = (m 1) K 2 f = (m 1) K W r. (16) Note that R(m, τ n ) represents a sampling aross range and the maximum unambiguous range ours when (m 1) = K (reall Equation (7)).
6 Inserting the result from Equation (16) into Equation (13) yields +j4π f(k 1) +j4πf1 R(m, τ n ) s(m, τ n ) = S(f k, τ n ) exp R(m, τ n ) exp ( ) +j4π f(k 1) m 1 = S(f k, τ n ) exp K +j4πf1 R(m, τ n ) exp 2 f +j2π(m 1) +j4πf1 R(m, τ n ) = S(f k, τ n ) exp (k 1) exp. K (17) Finally, realling the expression for the inverse disrete Fourier transform from Equation (15), the range profile may be expressed in terms of MATLAB s ifft. +j4πf1 R(m, τ n ) s(m, τ n ) = K ifft(s(f k, τ n )) exp (18) The onstant K appears in Equation (18) beause by definition, the inverse disrete Fourier transform ontains a saling fator of 1/K. In addition, it s important to note that the MATLAB implementation of the inverse disrete Fourier transform assumes that K is even. However, for values of K that are powers of 2, the omputational omplexity of the inverse disrete Fourier transform is optimized. By inserting the expression for R(m, τ n ) from Equation (16), an alternative expression for the implementation of a range profile in MATLAB is given by j2πf1 (m 1) s(m, τ n ) = K ifft(s(f k, τ n )) exp. (19) K f One final onsideration must be addressed. Eah pulse is motion ompensated suh that a satterer at the sene origin has zero phase and will appear in the zero frequeny bin in the range profile. This bin orresponds to the enter of the range profile omputed in Equation (19). By default, the funtion ifft omputes the values from 1 m K where m = 1 orresponds to the zero frequeny bin. To put m = 1 at the enter of the range profile, the ommand fftshift is applied to the output of the ifft ommand. This operation ensures that the output is orretly ordered following the inverse disrete Fourier transform operation and that the zero frequeny omponent orresponds to the enter of the output vetor. As a result, range index, m, aquires values between K/2 + 1 and K/2. Therefore, R(m, τ n ) of Equation (16) has values between W r /2 and W r /2 W r /K and Equation (19) beomes j2πf1 (m 1) s(m, τ n ) = K fftshift{ifft(s(f k, τ n ))} exp. (2) K f 4.2 Image Formation Proess For implementation of the bakprojetion imaging algorithm in MATLAB, Equation (6) is used to ompute the differential range, R(τ n ), for eah pixel, for eah pulse, where the pixel (x, y, z) oordinates are inserted into Equation (4). In order to use Equation (2) to form a SAR image, an interpolation step is required due to the fat that the disrete values of R(m, τ n ) do not orrespond exatly to the R(τ n ) values alulated for every pixel. There are several methods for performing this interpolation, 2,11 but here we will simply use the MATLAB interp1 funtion to implement linear interpolation. Sine s(m, τ n ) is a band-limited signal, the ideal interpolator is a sin interpolator. This an be approximated by zero-padding the ifft omputation and then performing linear interpolation on the inverse disrete Fourier transform output. A good rule of thumb is that the length of the ifft, denoted N fft, should be 1 times the length of the data, K. Also, the ifft funtion is most effiient when N fft is a power of 2. In the bakprojetion algorithm, N fft is provided as an input to the imaging funtion, as shown in Table 1.
7 The first step in the image formation proess is to implement Equation (2) for every pulse, zero-padding the data suh that the length of the ifft is N fft. Thus, j2πf1 (m 1) s(m, τ n ) = N fft fftshift{ifft(s(f k, τ n ))} exp. (21) N fft f where S(f k, τ n ) = for all k > K. To find the image response for a pixel at loation r, given pulse n, R(τ n ) is alulated and used to find an interpolated value of s(m, τ n ). This is denoted as s int (r, τ n ). The final image response, I(r), is simply the summation of these values for every pulse. N p I(r) = s int (r, τ n ). (22) n=1 Figure 1(b) shows an example of an image formed using the bakprojetion algorithm. The image was generated using the same data and equivalent pixel loations as those employed in Figure 1(a). 5. IMAGE EXAMPLES In the last few years, the Air Fore Researh Laboratory (AFRL) has released several datasets of both syntheti and measured SAR data. In this setion, the bakprojetion algorithm introdued in Setion 4 is applied to data from the Bakhoe Data Dome, the 2D/3D Volumetri Challenge Problem, and the SAR-based Ground Moving Target Indiator (GMTI) Challenge Problem. In addition, the imaging ode was implemented on the Civilian Vehile Radar Data Domes syntheti data set from The Ohio State University. 5.1 Bakhoe Data Dome In 24, AFRL released a synthetially generated data dome of a bakhoe target. 21 The sattered field data for the bakhoe was omputed on a 2π-steradian data dome over the target for a frequeny bandwidth of 5.9 GHz at a enter frequeny of 1 GHz y (m) x (m) Figure 2. Bakprojetion image (in db) using the Bakhoe Data Dome syntheti dataset. The bakprojetion algorithm was performed on this dataset and a resulting image is given in Figure 2. The data imaged onsisted of all 36 aspet angles from a single elevation of 1 using VV polarization. The full 5.9 GHz bandwidth was used to form the image. The sene extent is 1 m x 1 m with 2 m pixel spaing, resulting in a 51 x 51 pixel image. No windowing was applied to the data before imaging. The alulated maximum sene size is W r = m; W x = 9.27 m. 7
8 5.2 2D/3D Volumetri Challenge Problem In 27, AFRL released a hallenge problem for 2D/3D imaging of targets from a volumetri data set in an urban environment. 22 The data was olleted of a sene onsisting of numerous ivilian vehiles and alibration targets. The bakprojetion algorithm was performed on this dataset and a resulting image is given in Figure 3. The data imaged onsisted of Pass 1 with HH polarization. An integration angle, θ, of 4 entered at 4 azimuth was used. The sene extent is 1 m x 1 m with 2 m pixel spaing, resulting in a 51 x 51 pixel image. No windowing was applied to the data before imaging. The alulated maximum sene size is W r = m; W x = m and the alulated resolution is δr =.24 m; δx =.23 m y (m) x (m) Figure 3. Bakprojetion image (in db) using the 2D/3D Volumetri Challenge Problem dataset SAR-based GMTI Challenge Problem In 29, AFRL released a hallenge problem for SAR-based GMTI in urban environments. 23 This data onsists of a 71-seond portion of phase history data from a radar operating in irular SAR mode. The data has been range-gated around the known loation of a moving vehile, severely limiting the range swath of the data. However, during the first 25 seonds of the senario, the vehile is waiting in line to ross a busy intersetion. The bakprojetion algorithm was performed on this dataset and a resulting image is provided in Figure 4. The sene extent is 2 m x 2 m with 2 m pixel spaing, resulting in a 11 x 11 pixel image. The first 8 pulses (almost 4 seonds) were used, and no windowing was applied to the data before imaging. The alulated maximum sene size is W r = m; W x = 2, m, and the alulated resolution is δr =.23 m; δx =.25 m. The relatively small range extent aounts for the blak areas in the image. 5.4 Civilian Vehile Radar Data Domes In 21, The Ohio State University released a set of synthetially generated data domes of ivilian targets. 24 The sattered field data for the targets was omputed for full azimuth overage of elevation angles from 3 to 6. The data is full polarimetri with a frequeny bandwidth of 5.35 GHz at a enter frequeny of 9.6 GHz. The bakprojetion algorithm was implemented on the 1993 Jeep target, and a resulting image is shown in Figure 5. The data imaged onsisted of all 36 aspet angles from a single elevation of 3 using VV polarization. The full bandwidth was used to form the image. The sene extent is 1 m x 1 m with 2 m pixel spaing, resulting in a 51 x 51 pixel image. No windowing was applied to the data before imaging. The alulated maximum sene size is W r = 14.3 m; W x = m.
9 5 1 2 y (m) x (m) Figure 4. Bakprojetion image (in db) using the SAR-based GMTI dataset y (m) x (m) Figure 5. Bakprojetion image (in db) using the Civilian Vehile Radar Data Domes syntheti dataset CONCLUSIONS In this paper, we developed two simple SAR imaging algorithms: a mathed filter algorithm and a bakprojetion algorithm. The bakprojetion algorithm was derived from the mathed filter algorithm by utilizing a fast Fourier transform to ompute range profiles. Both algorithms were implemented in MATLAB, and image examples were provided using four publily available SAR datasets.
10 A.1 Mathed Filter Algorithm APPENDIX A. MATLAB CODE 1 funtion data = mfbasi(data) 2 3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 4 % This funtion performs a mathed filter operation. The following % 5 % fields need to be populated: % 6 % % 7 % data.deltaf: Step size of frequeny data (Hz) % 8 % data.minf: Vetor ontaining the start frequeny of eah pulse (Hz) % 9 % data.x mat: The x position of eah pixel (m) % 1 % data.y mat: The y position of eah pixel (m) % 11 % data.z mat: The z position of eah pixel (m) % 12 % data.antx: The x position of the sensor at eah pulse (m) % 13 % data.anty: The y position of the sensor at eah pulse (m) % 14 % data.antz: The z position of the sensor at eah pulse (m) % 15 % data.r: The range to sene enter (m) % 16 % data.phdata: Phase history data (frequeny domain) % 17 % Fast time in rows, slow time in olumns % 18 % % 19 % The output is: % 2 % data.im final: The omplex image value at eah pixel % 21 % % 22 % Written by LeRoy Gorham, Air Fore Researh Laboratory, WPAFB, OH % 23 % leroy.gorham@wpafb.af.mil % 24 % Date Released: 8 Apr 21 % 25 % % 26 % Gorham, L.A. and Moore, L.J., "SAR image formation toolbox for % 27 % MATLAB," Algorithms for Syntheti Aperture Radar Imagery XVII % 28 % 7669, SPIE (21). % 29 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 3 31 % Define speed of light (m/s) 32 = ; % Determine the size of the phase history data 35 data.k = size(data.phdata,1); % The number of frequeny bins per pulse 36 data.np = size(data.phdata,2); % The number of pulses % Determine the azimuth angles of the image pulses (radians) 39 data.antaz = unwrap(atan2(data.anty,data.antx)); 4 41 % Determine the average azimuth angle step size (radians) 42 data.deltaaz = abs(mean(diff(data.antaz))); % Determine the total azimuth angle of the aperture (radians) 45 data.totalaz = max(data.antaz) min(data.antaz); % Determine the maximum sene size of the image (m) 48 data.maxwr = /(2*data.deltaF); 49 data.maxwx = /(2*data.deltaAz*mean(data.minF)); 5 51 % Determine the resolution of the image (m) 52 data.dr = /(2*data.deltaF*data.K); 53 data.dx = /(2*data.totalAz*mean(data.minF)); % Display maximum sene size and resolution 56 fprintf('maximum Sene Size: %.2f m range, %.2f m ross range\n',data.maxwr,data.maxwx); 57 fprintf('resolution: %.2fm range, %.2f m ross range\n',data.dr,data.dx); % Initialize the image with all zero values 6 data.im final = zeros(size(data.x mat)); % Set up a vetor to keep exeution times for eah pulse (se) 63 t = zeros(1,data.np); % Loop through every pulse 66 for ii = 1:data.Np % Display status of the imaging proess 69 if ii > 1 7 t sofar = sum(t(1:(ii 1))); 71 t est = (t sofar*data.np/(ii 1) t sofar)/6; 72 fprintf('pulse %d of %d, %.2f minutes remaining\n',ii,data.np,t est); 73 else
11 74 fprintf('pulse %d of %d\n',ii,data.np); 75 end 76 ti % Calulate differential range for eah pixel in the image (m) 79 dr = sqrt((data.antx(ii) data.x mat).ˆ (data.anty(ii) data.y mat).ˆ (data.antz(ii) data.z mat).ˆ2) data.r(ii); % Calulate the frequeny of eah sample in the pulse (Hz) 84 freq = data.minf(ii) + (:(data.k 1)) * data.deltaf; % Perform the Mathed Filter operation 87 for jj = 1:data.K 88 data.im final = data.im final + data.phdata(jj,ii) * exp(1i*4*pi*freq(jj)/*dr); 89 end 9 91 % Determine the exeution time for this pulse 92 t(ii) = to; 93 end return A.2 Bakprojetion Algorithm 1 funtion data = bpbasi(data) 2 3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 4 % This funtion performs a basi Bakprojetion operation. The % 5 % following fields need to be populated: % 6 % % 7 % data.nfft: Size of the FFT to form the range profile % 8 % data.deltaf: Step size of frequeny data (Hz) % 9 % data.minf: Vetor ontaining the start frequeny of eah pulse (Hz) % 1 % data.x mat: The x position of eah pixel (m) % 11 % data.y mat: The y position of eah pixel (m) % 12 % data.z mat: The z position of eah pixel (m) % 13 % data.antx: The x position of the sensor at eah pulse (m) % 14 % data.anty: The y position of the sensor at eah pulse (m) % 15 % data.antz: The z position of the sensor at eah pulse (m) % 16 % data.r: The range to sene enter (m) % 17 % data.phdata: Phase history data (frequeny domain) % 18 % Fast time in rows, slow time in olumns % 19 % % 2 % The output is: % 21 % data.im final: The omplex image value at eah pixel % 22 % % 23 % Written by LeRoy Gorham, Air Fore Researh Laboratory, WPAFB, OH % 24 % leroy.gorham@wpafb.af.mil % 25 % Date Released: 8 Apr 21 % 26 % % 27 % Gorham, L.A. and Moore, L.J., "SAR image formation toolbox for % 28 % MATLAB," Algorithms for Syntheti Aperture Radar Imagery XVII % 29 % 7669, SPIE (21). % 3 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% % Define speed of light (m/s) 33 = ; % Determine the size of the phase history data 36 data.k = size(data.phdata,1); % The number of frequeny bins per pulse 37 data.np = size(data.phdata,2); % The number of pulses % Determine the azimuth angles of the image pulses (radians) 4 data.antaz = unwrap(atan2(data.anty,data.antx)); % Determine the average azimuth angle step size (radians) 43 data.deltaaz = abs(mean(diff(data.antaz))); % Determine the total azimuth angle of the aperture (radians) 46 data.totalaz = max(data.antaz) min(data.antaz); % Determine the maximum sene size of the image (m) 49 data.maxwr = /(2*data.deltaF); 5 data.maxwx = /(2*data.deltaAz*mean(data.minF)); 51
12 52 % Determine the resolution of the image (m) 53 data.dr = /(2*data.deltaF*data.K); 54 data.dx = /(2*data.totalAz*mean(data.minF)); % Display maximum sene size and resolution 57 fprintf('maximum Sene Size: %.2f m range, %.2f m ross range\n',data.maxwr,data.maxwx); 58 fprintf('resolution: %.2fm range, %.2f m ross range\n',data.dr,data.dx); 59 6 % Calulate the range to every bin in the range profile (m) 61 data.r ve = linspae( data.nfft/2,data.nfft/2 1,data.Nfft)*data.maxWr/data.Nfft; % Initialize the image with all zero values 64 data.im final = zeros(size(data.x mat)); % Set up a vetor to keep exeution times for eah pulse (se) 67 t = zeros(1,data.np); % Loop through every pulse 7 for ii = 1:data.Np % Display status of the imaging proess 73 if ii > 1 74 t sofar = sum(t(1:(ii 1))); 75 t est = (t sofar*data.np/(ii 1) t sofar)/6; 76 fprintf('pulse %d of %d, %.2f minutes remaining\n',ii,data.np,t est); 77 else 78 fprintf('pulse %d of %d\n',ii,data.np); 79 end 8 ti % Form the range profile with zero padding added 83 r = fftshift(ifft(data.phdata(:,ii),data.nfft)); % Calulate differential range for eah pixel in the image (m) 86 dr = sqrt((data.antx(ii) data.x mat).ˆ (data.anty(ii) data.y mat).ˆ (data.antz(ii) data.z mat).ˆ2) data.r(ii); 89 9 % Calulate phase orretion for image 91 phcorr = exp(1i*4*pi*data.minf(ii)/*dr); % Determine whih pixels fall within the range swath 94 I = find(and(dr > min(data.r ve), dr < max(data.r ve))); % Update the image using linear interpolation 97 data.im final(i) = data.im final(i) + interp1(data.r ve,r,dr(i),'linear').* phcorr(i); % Determine the exeution time for this pulse 1 t(ii) = to; 11 end return REFERENCES [1] Soumekh, M., [Syntheti aperture radar signal proessing with MATLAB algorithms], Wiley-Intersiene (1999). [2] Jakowatz Jr, C., Wahl, D., Yoky, D., Bray, B., Bow Jr, W., and Rihards, J., Comparison of algorithms for use in real-time spotlight-mode SAR image formation, in [Proeedings of SPIE], 5427, 18 (24). [3] Jakowatz Jr, C. and Doren, N., Comparison of polar formatting and bak-projetion algorithms for spotlight-mode SAR image formation, in [Proeedings of SPIE], 6237, 6237H (26). [4] Carrara, W., Goodman, R., and Majewski, R., [Spotlight Syntheti Aperture Radar - Signal Proessing Algorithms], Norwood, MA: Arteh House (1995). [5] Jakowatz, C., Wahl, D., Eihel, P., Ghiglia, D., and Thompson, P., [Spotlight-mode syntheti aperture radar: a signal proessing approah], Kluwer Aademi Pub (1996). [6] Stark, H., Woods, J., Paul, I., and Hingorani, R., An investigation of omputerized tomography by diret fourier inversion and optimum interpolation, IEEE Transations on Biomedial Engineering, (1981). [7] Munson Jr, D., O Brien, J., and Jenkins, W., A tomographi formulation of spotlight-mode syntheti aperture radar, Proeedings of the IEEE 71(8), (1983).
13 [8] Gorham, L. A., Rigling, B. D., and Zelnio, E. G., A omparison between imaging radar and medial imaging polar format algorithm implementations, in [Proeedings of SPIE], 6568, 6568K (27). [9] Doren, N., Jakowatz Jr, C., Wahl, D., and Thompson, P., General formulation for wavefront urvature orretion in polar-formatted spotlight-mode SAR images using spae-variant post-filtering, in [Proeedings of the 1997 International Conferene on Image Proessing (ICIP 97) 3-Volume Set-Volume 1-Volume 1], IEEE Computer Soiety Washington, DC, USA (1997). [1] Preiss, M., Gray, D., and Stay, N., Spae Variant Filtering of Polar Format Spotlight SAR Images for Wavefront Curvature Corretion and Interferometri Proessing, in [IGARSS 22: IEEE International Geosiene and Remote Sensing Symposium (24th: 22: Toronto, Ontario)], IEEE: Institute of Eletrial and Eletronis Engineers (22). [11] Desai, M. and Jenkins, W., Convolution bakprojetion image reonstrution for spotlight mode syntheti aperture radar, IEEE Transations on Image Proessing 1(4), (1992). [12] Hartley, T., Fasih, A., Berdanier, C., Özgüner, F., and Çatalyürek, U., Investigating the use of GPUaelerated nodes for SAR image formation, in [Cluster Computing and Workshops, 29. CLUSTER 9. IEEE International Conferene on], 1 8 (Aug 31 Sep 4 29). [13] Rogan, A. and Carande, R., Improving the fast bakprojetion algorithm through massive parallelizations, Radar Sensor Tehnology XIV 7669(1), SPIE (21). [14] Yegulalp, A., Fast bakprojetion algorithm for syntheti aperture radar, in [Radar Conferene, The Reord of the 1999 IEEE], 6 65 (1999). [15] MCorkle, J. and Rofheart, M., An order N 2 log(n) bakprojetor algorithm for fousing wide-angle widebandwidth arbitrary-motion syntheti aperture radar, in [Pro. of SPIE Conferene on Radar Sensor Tehnology], 2747, (1996). [16] Ulander, L., Hellsten, H., and Stenstrom, G., Syntheti-aperture radar proessing using fast fatorized bak-projetion, IEEE Transations on Aerospae and eletroni systems 39(3), (23). [17] Wahl, D. E., Yoky, D. A., and Jakowatz Jr., C. V., An implementation of a fast bakprojetion image formation algorithm for spotlight-mode SAR, in [Proeedings of SPIE], 697, 697H (28). [18] Rigling, B. and Moses, R., Taylor expansion of the differential range for monostati SAR, IEEE Transations on Aerospae and Eletroni Systems 41(1), 6 64 (25). [19] The Mathworks, In, MATLAB Doumentation: fft - Disrete Fourier Transform. Version (R29b). [2] Crohiere, R. and Rabiner, L., Interpolation and deimation of digital signals - a tutorial review, Proeedings of the IEEE 69(3), (1981). [21] Naidu, K. and Lin, L., Data dome: full k-spae sampling data for high-frequeny radar researh, Algorithms for Syntheti Aperture Radar Imagery XI 5427(1), 2 27, SPIE (24). [22] Casteel Jr, C. H., Gorham, L. A., Minardi, M. J., Sarborough, S. M., Naidu, K. D., and Majumder, U. K., A hallenge problem for 2D/3D imaging of targets from a volumetri data set in an urban environment, Algorithms for Syntheti Aperture Radar Imagery XIV 6568(1), 6568D, SPIE (27). [23] Sarborough, S. M., Casteel Jr, C. H., Gorham, L., Minardi, M. J., Majumder, U. K., Judge, M. G., Zelnio, E., Bryant, M., Nihols, H., and Page, D., A hallenge problem for SAR-based GMTI in urban environments, Algorithms for Syntheti Aperture Radar Imagery XVI 7337(1), 7337G, SPIE (29). [24] Dungan, K. E., Austin, C., Nehrbass, J., and Potter, L. C., Civilian vehile radar data domes, Algorithms for Syntheti Aperture Radar Imagery XVII 7699(1), SPIE (21).
PARAMETRIC SAR IMAGE FORMATION - A PROMISING APPROACH TO RESOLUTION-UNLIMITED IMAGING. Yesheng Gao, Kaizhi Wang, Xingzhao Liu
20th European Signal Proessing Conferene EUSIPCO 2012) Buharest, Romania, August 27-31, 2012 PARAMETRIC SAR IMAGE FORMATION - A PROMISING APPROACH TO RESOLUTION-UNLIMITED IMAGING Yesheng Gao, Kaizhi Wang,
More informationRANGE DOPPLER ALGORITHM FOR BISTATIC SAR PROCESSING BASED ON THE IMPROVED LOFFELD S BISTATIC FORMULA
Progress In Eletromagnetis Researh Letters, Vol. 27, 161 169, 2011 RANGE DOPPLER ALGORITHM FOR ISTATIC SAR PROCESSING ASED ON THE IMPROVED LOFFELD S ISTATIC FORMULA X. Wang 1, * and D. Y. Zhu 2 1 Nanjing
More informationDetecting Moving Targets in Clutter in Airborne SAR via Keystoning and Multiple Phase Center Interferometry
Deteting Moving Targets in Clutter in Airborne SAR via Keystoning and Multiple Phase Center Interferometry D. M. Zasada, P. K. Sanyal The MITRE Corp., 6 Eletroni Parkway, Rome, NY 134 (dmzasada, psanyal)@mitre.org
More informationFocusing Translational Variant Bistatic Forward-Looking SAR Data Based on Two-Dimensional Non-Uniform FFT
Progress In Eletromagnetis Researh M, Vol. 37, 1 10, 2014 Fousing Translational Variant Bistati Forward-Looking SAR Data Based on Two-Dimensional Non-Uniform FFT Chan Liu 1, Shunsheng Zhang 1, *, Chunyang
More informationAn Extended Wavenumber-Domain Algorithm Combined with Two-Step Motion Compensation for Bistatic Forward-Looking SAR
Progress In Eletromagnetis Researh Letters, Vol. 63, 85 92, 2016 An Extended Wavenumber-Domain Algorithm Combined with Two-Step Motion Compensation for Bistati Forward-Looking SAR Yuebo Zha 1, * and Wei
More informationNONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION. Ken Sauer and Charles A. Bouman
NONLINEAR BACK PROJECTION FOR TOMOGRAPHIC IMAGE RECONSTRUCTION Ken Sauer and Charles A. Bouman Department of Eletrial Engineering, University of Notre Dame Notre Dame, IN 46556, (219) 631-6999 Shool of
More informationIMPROVEMENTS to processing methods for synthetic
1 Generalized SAR Proessing and Motion Compensation Evan C. Zaugg Brigham Young University Mirowave Earth Remote Sensing Laboratory 459 Clyde Building, Provo, UT 846 81-4-4884 zaugg@mers.byu.edu Abstrat
More informationA Challenge Problem for 2D/3D Imaging of Targets from a Volumetric Data Set in an Urban Environment
A Challenge Problem for 2D/3D Imaging of Targets from a Volumetric Data Set in an Urban Environment Curtis H. Casteel, Jr,*, LeRoy A. Gorham, Michael J. Minardi, Steven M. Scarborough, Kiranmai D. Naidu,
More informationCONSIDERING the high computational load of synthetic
IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, VOL. 5, NO. 4, OCTOBER 2008 573 Doppler Keystone Transform: An Approah Suitable for Parallel Implementation of SAR Moving Target Imaging Gang Li, Member, IEEE,
More informationTHE time-domain backprojection (BP) algorithm has already
1460 IEEE GEOSCIENCE AND REMOTE SENSING LETTERS, VOL. 11, NO. 9, SEPTEMBER 014 A Fast BP Algorithm With Wavenumber Spetrum Fusion for High-Resolution Spotlight SAR Imaging Lei Zhang, Member, IEEE, Hao-lin
More informationAnd, the (low-pass) Butterworth filter of order m is given in the frequency domain by
Problem Set no.3.a) The ideal low-pass filter is given in the frequeny domain by B ideal ( f ), f f; =, f > f. () And, the (low-pass) Butterworth filter of order m is given in the frequeny domain by B
More informationDrawing lines. Naïve line drawing algorithm. drawpixel(x, round(y)); double dy = y1 - y0; double dx = x1 - x0; double m = dy / dx; double y = y0;
Naïve line drawing algorithm // Connet to grid points(x0,y0) and // (x1,y1) by a line. void drawline(int x0, int y0, int x1, int y1) { int x; double dy = y1 - y0; double dx = x1 - x0; double m = dy / dx;
More informationA Novel Validity Index for Determination of the Optimal Number of Clusters
IEICE TRANS. INF. & SYST., VOL.E84 D, NO.2 FEBRUARY 2001 281 LETTER A Novel Validity Index for Determination of the Optimal Number of Clusters Do-Jong KIM, Yong-Woon PARK, and Dong-Jo PARK, Nonmembers
More informationAZIMUTH NONLINEAR CHIRP SCALING INTEGRATED WITH RANGE CHIRP SCALING ALGORITHM FOR HIGHLY SQUINTED SAR IMAGING
Progress In Eletromagnetis Researh, Vol. 143, 165 185, 2013 AZIMUTH NONLINEAR CHIRP SCALING INTEGRATED WITH RANGE CHIRP SCALING ALGORITHM FOR HIGHLY SQUINTED SAR IMAGING Qinglin Zhai *, Wei Wang, Jiemin
More informationA radiometric analysis of projected sinusoidal illumination for opaque surfaces
University of Virginia tehnial report CS-21-7 aompanying A Coaxial Optial Sanner for Synhronous Aquisition of 3D Geometry and Surfae Refletane A radiometri analysis of projeted sinusoidal illumination
More informationRotation Invariant Spherical Harmonic Representation of 3D Shape Descriptors
Eurographis Symposium on Geometry Proessing (003) L. Kobbelt, P. Shröder, H. Hoppe (Editors) Rotation Invariant Spherial Harmoni Representation of 3D Shape Desriptors Mihael Kazhdan, Thomas Funkhouser,
More informationDATA ACQUISITION AND PROCESSING OF PARALLEL FREQUENCY SAR BASED ON COMPRESSIVE SENSING
Progress In Eletromagnetis Researh, Vol. 133, 199 215, 2013 DATA ACQUISITION AND PROCESSING OF PARALLEL FREQUENCY SAR BASED ON COMPRESSIVE SENSING Y. N. You, H. P. Xu *, C. S. Li, and L. Q. Zhang Shool
More informationPlot-to-track correlation in A-SMGCS using the target images from a Surface Movement Radar
Plot-to-trak orrelation in A-SMGCS using the target images from a Surfae Movement Radar G. Golino Radar & ehnology Division AMS, Italy ggolino@amsjv.it Abstrat he main topi of this paper is the formulation
More informationSmooth Trajectory Planning Along Bezier Curve for Mobile Robots with Velocity Constraints
Smooth Trajetory Planning Along Bezier Curve for Mobile Robots with Veloity Constraints Gil Jin Yang and Byoung Wook Choi Department of Eletrial and Information Engineering Seoul National University of
More informationGraph-Based vs Depth-Based Data Representation for Multiview Images
Graph-Based vs Depth-Based Data Representation for Multiview Images Thomas Maugey, Antonio Ortega, Pasal Frossard Signal Proessing Laboratory (LTS), Eole Polytehnique Fédérale de Lausanne (EPFL) Email:
More informationPerformance of Histogram-Based Skin Colour Segmentation for Arms Detection in Human Motion Analysis Application
World Aademy of Siene, Engineering and Tehnology 8 009 Performane of Histogram-Based Skin Colour Segmentation for Arms Detetion in Human Motion Analysis Appliation Rosalyn R. Porle, Ali Chekima, Farrah
More informationParticle Swarm Optimization for the Design of High Diffraction Efficient Holographic Grating
Original Artile Partile Swarm Optimization for the Design of High Diffration Effiient Holographi Grating A.K. Tripathy 1, S.K. Das, M. Sundaray 3 and S.K. Tripathy* 4 1, Department of Computer Siene, Berhampur
More informationA 3D Polar Processing Algorithm for Scale Model UHF ISAR Imaging
A 3D Polar Proessing Algorithm for Sale Model UH ISAR Imaging Christopher J. Beaudoin* a, Andrew J. Gatesman a, Robert H. Giles a, Jerry Waldman a, and William E. Nixon b a Submillimeter-Wave Tehnology
More informationChapter 2: Introduction to Maple V
Chapter 2: Introdution to Maple V 2-1 Working with Maple Worksheets Try It! (p. 15) Start a Maple session with an empty worksheet. The name of the worksheet should be Untitled (1). Use one of the standard
More informationCompressed Sensing mm-wave SAR for Non-Destructive Testing Applications using Side Information
Compressed Sensing mm-wave SAR for Non-Destrutive Testing Appliations using Side Information Mathias Bequaert 1,2, Edison Cristofani 1,2, Gokarna Pandey 2, Marijke Vandewal 1, Johan Stiens 2,3 and Nikos
More informationCONVOLUTION BACK-PROJECTION IMAGING ALGO- RITHM FOR DOWNWARD-LOOKING SPARSE LINEAR ARRAY THREE DIMENSIONAL SYNTHETIC APER- TURE RADAR
Progress In Eletromagnetis Researh, Vol. 129, 287 313, 2012 CONVOLUTION BACK-PROJECTION IMAGING ALGO- RITHM FOR DOWNWARD-LOOKING SPARSE LINEAR ARRAY THREE DIMENSIONAL SYNTHETIC APER- TURE RADAR X. M. Peng
More informationDetection and Recognition of Non-Occluded Objects using Signature Map
6th WSEAS International Conferene on CIRCUITS, SYSTEMS, ELECTRONICS,CONTROL & SIGNAL PROCESSING, Cairo, Egypt, De 9-31, 007 65 Detetion and Reognition of Non-Oluded Objets using Signature Map Sangbum Park,
More informationA {k, n}-secret Sharing Scheme for Color Images
A {k, n}-seret Sharing Sheme for Color Images Rastislav Luka, Konstantinos N. Plataniotis, and Anastasios N. Venetsanopoulos The Edward S. Rogers Sr. Dept. of Eletrial and Computer Engineering, University
More informationTime delay estimation of reverberant meeting speech: on the use of multichannel linear prediction
University of Wollongong Researh Online Faulty of Informatis - apers (Arhive) Faulty of Engineering and Information Sienes 7 Time delay estimation of reverberant meeting speeh: on the use of multihannel
More informationPipelined Multipliers for Reconfigurable Hardware
Pipelined Multipliers for Reonfigurable Hardware Mithell J. Myjak and José G. Delgado-Frias Shool of Eletrial Engineering and Computer Siene, Washington State University Pullman, WA 99164-2752 USA {mmyjak,
More informationThe Mathematics of Simple Ultrasonic 2-Dimensional Sensing
The Mathematis of Simple Ultrasoni -Dimensional Sensing President, Bitstream Tehnology The Mathematis of Simple Ultrasoni -Dimensional Sensing Introdution Our ompany, Bitstream Tehnology, has been developing
More informationVideo Data and Sonar Data: Real World Data Fusion Example
14th International Conferene on Information Fusion Chiago, Illinois, USA, July 5-8, 2011 Video Data and Sonar Data: Real World Data Fusion Example David W. Krout Applied Physis Lab dkrout@apl.washington.edu
More information3-D IMAGE MODELS AND COMPRESSION - SYNTHETIC HYBRID OR NATURAL FIT?
3-D IMAGE MODELS AND COMPRESSION - SYNTHETIC HYBRID OR NATURAL FIT? Bernd Girod, Peter Eisert, Marus Magnor, Ekehard Steinbah, Thomas Wiegand Te {girod eommuniations Laboratory, University of Erlangen-Nuremberg
More informationIntroduction to Seismology Spring 2008
MIT OpenCourseWare http://ow.mit.edu 1.510 Introdution to Seismology Spring 008 For information about iting these materials or our Terms of Use, visit: http://ow.mit.edu/terms. 1.510 Leture Notes 3.3.007
More informationGray Codes for Reflectable Languages
Gray Codes for Refletable Languages Yue Li Joe Sawada Marh 8, 2008 Abstrat We lassify a type of language alled a refletable language. We then develop a generi algorithm that an be used to list all strings
More informationINTERPOLATED AND WARPED 2-D DIGITAL WAVEGUIDE MESH ALGORITHMS
Proeedings of the COST G-6 Conferene on Digital Audio Effets (DAFX-), Verona, Italy, Deember 7-9, INTERPOLATED AND WARPED -D DIGITAL WAVEGUIDE MESH ALGORITHMS Vesa Välimäki Lab. of Aoustis and Audio Signal
More informationUsing Augmented Measurements to Improve the Convergence of ICP
Using Augmented Measurements to Improve the onvergene of IP Jaopo Serafin, Giorgio Grisetti Dept. of omputer, ontrol and Management Engineering, Sapienza University of Rome, Via Ariosto 25, I-0085, Rome,
More informationMeasurement of the stereoscopic rangefinder beam angular velocity using the digital image processing method
Measurement of the stereosopi rangefinder beam angular veloity using the digital image proessing method ROMAN VÍTEK Department of weapons and ammunition University of defense Kouniova 65, 62 Brno CZECH
More information1. Introduction. 2. The Probable Stope Algorithm
1. Introdution Optimization in underground mine design has reeived less attention than that in open pit mines. This is mostly due to the diversity o underground mining methods and omplexity o underground
More informationWe P9 16 Eigenray Tracing in 3D Heterogeneous Media
We P9 Eigenray Traing in 3D Heterogeneous Media Z. Koren* (Emerson), I. Ravve (Emerson) Summary Conventional two-point ray traing in a general 3D heterogeneous medium is normally performed by a shooting
More informationParallelization and Performance of 3D Ultrasound Imaging Beamforming Algorithms on Modern Clusters
Parallelization and Performane of 3D Ultrasound Imaging Beamforming Algorithms on Modern Clusters F. Zhang, A. Bilas, A. Dhanantwari, K.N. Plataniotis, R. Abiprojo, and S. Stergiopoulos Dept. of Eletrial
More informationApproximate logic synthesis for error tolerant applications
Approximate logi synthesis for error tolerant appliations Doohul Shin and Sandeep K. Gupta Eletrial Engineering Department, University of Southern California, Los Angeles, CA 989 {doohuls, sandeep}@us.edu
More informationExtracting Partition Statistics from Semistructured Data
Extrating Partition Statistis from Semistrutured Data John N. Wilson Rihard Gourlay Robert Japp Mathias Neumüller Department of Computer and Information Sienes University of Strathlyde, Glasgow, UK {jnw,rsg,rpj,mathias}@is.strath.a.uk
More informationAccommodations of QoS DiffServ Over IP and MPLS Networks
Aommodations of QoS DiffServ Over IP and MPLS Networks Abdullah AlWehaibi, Anjali Agarwal, Mihael Kadoh and Ahmed ElHakeem Department of Eletrial and Computer Department de Genie Eletrique Engineering
More informationFUZZY WATERSHED FOR IMAGE SEGMENTATION
FUZZY WATERSHED FOR IMAGE SEGMENTATION Ramón Moreno, Manuel Graña Computational Intelligene Group, Universidad del País Vaso, Spain http://www.ehu.es/winto; {ramon.moreno,manuel.grana}@ehu.es Abstrat The
More informationFitting conics to paracatadioptric projections of lines
Computer Vision and Image Understanding 11 (6) 11 16 www.elsevier.om/loate/viu Fitting onis to paraatadioptri projetions of lines João P. Barreto *, Helder Araujo Institute for Systems and Robotis, Department
More informationOn - Line Path Delay Fault Testing of Omega MINs M. Bellos 1, E. Kalligeros 1, D. Nikolos 1,2 & H. T. Vergos 1,2
On - Line Path Delay Fault Testing of Omega MINs M. Bellos, E. Kalligeros, D. Nikolos,2 & H. T. Vergos,2 Dept. of Computer Engineering and Informatis 2 Computer Tehnology Institute University of Patras,
More informationA DYNAMIC ACCESS CONTROL WITH BINARY KEY-PAIR
Malaysian Journal of Computer Siene, Vol 10 No 1, June 1997, pp 36-41 A DYNAMIC ACCESS CONTROL WITH BINARY KEY-PAIR Md Rafiqul Islam, Harihodin Selamat and Mohd Noor Md Sap Faulty of Computer Siene and
More information特集 Road Border Recognition Using FIR Images and LIDAR Signal Processing
デンソーテクニカルレビュー Vol. 15 2010 特集 Road Border Reognition Using FIR Images and LIDAR Signal Proessing 高木聖和 バーゼル ファルディ Kiyokazu TAKAGI Basel Fardi ヘンドリック ヴァイゲル Hendrik Weigel ゲルド ヴァニーリック Gerd Wanielik This paper
More informationHEXA: Compact Data Structures for Faster Packet Processing
Washington University in St. Louis Washington University Open Sholarship All Computer Siene and Engineering Researh Computer Siene and Engineering Report Number: 27-26 27 HEXA: Compat Data Strutures for
More informationAbstract. Key Words: Image Filters, Fuzzy Filters, Order Statistics Filters, Rank Ordered Mean Filters, Channel Noise. 1.
Fuzzy Weighted Rank Ordered Mean (FWROM) Filters for Mixed Noise Suppression from Images S. Meher, G. Panda, B. Majhi 3, M.R. Meher 4,,4 Department of Eletronis and I.E., National Institute of Tehnology,
More informationSupplementary Material: Geometric Calibration of Micro-Lens-Based Light-Field Cameras using Line Features
Supplementary Material: Geometri Calibration of Miro-Lens-Based Light-Field Cameras using Line Features Yunsu Bok, Hae-Gon Jeon and In So Kweon KAIST, Korea As the supplementary material, we provide detailed
More informationAcoustic Links. Maximizing Channel Utilization for Underwater
Maximizing Channel Utilization for Underwater Aousti Links Albert F Hairris III Davide G. B. Meneghetti Adihele Zorzi Department of Information Engineering University of Padova, Italy Email: {harris,davide.meneghetti,zorzi}@dei.unipd.it
More informationA scheme for racquet sports video analysis with the combination of audio-visual information
A sheme for raquet sports video analysis with the ombination of audio-visual information Liyuan Xing a*, Qixiang Ye b, Weigang Zhang, Qingming Huang a and Hua Yu a a Graduate Shool of the Chinese Aadamy
More informationLearning Convention Propagation in BeerAdvocate Reviews from a etwork Perspective. Abstract
CS 9 Projet Final Report: Learning Convention Propagation in BeerAdvoate Reviews from a etwork Perspetive Abstrat We look at the way onventions propagate between reviews on the BeerAdvoate dataset, and
More informationChromaticity-matched Superimposition of Foreground Objects in Different Environments
FCV216, the 22nd Korea-Japan Joint Workshop on Frontiers of Computer Vision Chromatiity-mathed Superimposition of Foreground Objets in Different Environments Yohei Ogura Graduate Shool of Siene and Tehnology
More informationWide Angle, Staring Synthetic Aperture Radar
88 ABW-12-0578 Wide Angle, Staring Synthetic Aperture Radar Feb 2012 Ed Zelnio Sensors Directorate Air Force Research Laboratory Outline Review SAR Focus on Wide Angle, Staring SAR (90%) Technology Challenges
More informationAn Enhanced Forward-Looking SAR Imaging Algorithm Based on Compressive Sensing
Progress In Eletromagnetis Researh M, Vol. 78, 69 81, 2019 An Enhaned Forward-Looking SAR Imaging Algorithm Based on Compressive Sensing Bo Pang, Hao Wu *, Shiqi Xing, Dahai Dai, Yongzhen Li, and Xuesong
More informationDetecting Outliers in High-Dimensional Datasets with Mixed Attributes
Deteting Outliers in High-Dimensional Datasets with Mixed Attributes A. Koufakou, M. Georgiopoulos, and G.C. Anagnostopoulos 2 Shool of EECS, University of Central Florida, Orlando, FL, USA 2 Dept. of
More informationCluster-Based Cumulative Ensembles
Cluster-Based Cumulative Ensembles Hanan G. Ayad and Mohamed S. Kamel Pattern Analysis and Mahine Intelligene Lab, Eletrial and Computer Engineering, University of Waterloo, Waterloo, Ontario N2L 3G1,
More informationAlgorithms for External Memory Lecture 6 Graph Algorithms - Weighted List Ranking
Algorithms for External Memory Leture 6 Graph Algorithms - Weighted List Ranking Leturer: Nodari Sithinava Sribe: Andi Hellmund, Simon Ohsenreither 1 Introdution & Motivation After talking about I/O-effiient
More informationExploiting Enriched Contextual Information for Mobile App Classification
Exploiting Enrihed Contextual Information for Mobile App Classifiation Hengshu Zhu 1 Huanhuan Cao 2 Enhong Chen 1 Hui Xiong 3 Jilei Tian 2 1 University of Siene and Tehnology of China 2 Nokia Researh Center
More informationFast image-formation algorithm for ultrahigh-resolution airborne squint spotlight synthetic aperture radar based on adaptive sliding receivewindow
Fast image-formation algorithm for ultrahigh-resolution airborne squint spotlight syntheti aperture radar based on adaptive sliding reeivewindow tehnique Wei Yang Hong-heng Zeng Jie Chen Peng-bo Wang Fast
More informationAnalysis of input and output configurations for use in four-valued CCD programmable logic arrays
nalysis of input and output onfigurations for use in four-valued D programmable logi arrays J.T. utler H.G. Kerkhoff ndexing terms: Logi, iruit theory and design, harge-oupled devies bstrat: s in binary,
More informationarxiv: v1 [cs.db] 13 Sep 2017
An effiient lustering algorithm from the measure of loal Gaussian distribution Yuan-Yen Tai (Dated: May 27, 2018) In this paper, I will introdue a fast and novel lustering algorithm based on Gaussian distribution
More informationCompressive Sensing Ultrasound Beamformed Imaging In Time and Frequency Domain
Compressive Sensing Ultrasound Beamformed Imaging In Time and Frequeny Domain Pradeep Kumar Mishra, R Bharath, P Rajalakshmi, Uday B Desai Department of Eletrial Engineering, Indian Institute of Tehnology,
More informationA PROTOTYPE OF INTELLIGENT VIDEO SURVEILLANCE CAMERAS
INTERNATIONAL JOURNAL OF INFORMATION AND SYSTEMS SCIENCES Volume 1, 3, Number 1, 3, Pages 1-22 365-382 2007 Institute for Sientifi Computing and Information A PROTOTYPE OF INTELLIGENT VIDEO SURVEILLANCE
More informationA Load-Balanced Clustering Protocol for Hierarchical Wireless Sensor Networks
International Journal of Advanes in Computer Networks and Its Seurity IJCNS A Load-Balaned Clustering Protool for Hierarhial Wireless Sensor Networks Mehdi Tarhani, Yousef S. Kavian, Saman Siavoshi, Ali
More informationKERNEL SPARSE REPRESENTATION WITH LOCAL PATTERNS FOR FACE RECOGNITION
KERNEL SPARSE REPRESENTATION WITH LOCAL PATTERNS FOR FACE RECOGNITION Cuiui Kang 1, Shengai Liao, Shiming Xiang 1, Chunhong Pan 1 1 National Laboratory of Pattern Reognition, Institute of Automation, Chinese
More informationStable Road Lane Model Based on Clothoids
Stable Road Lane Model Based on Clothoids C Gakstatter*, S Thomas**, Dr P Heinemann*, Prof Gudrun Klinker*** *Audi Eletronis Venture GmbH, **Leibniz Universität Hannover, ***Tehnishe Universität Münhen
More informationHIGHER ORDER full-wave three-dimensional (3-D) large-domain techniques in
FACTA UNIVERSITATIS (NIŠ) SER.: ELEC. ENERG. vol. 21, no. 2, August 2008, 209-220 Comparison of Higher Order FEM and MoM/SIE Approahes in Analyses of Closed- and Open-Region Eletromagneti Problems Milan
More informationTrajectory Tracking Control for A Wheeled Mobile Robot Using Fuzzy Logic Controller
Trajetory Traking Control for A Wheeled Mobile Robot Using Fuzzy Logi Controller K N FARESS 1 M T EL HAGRY 1 A A EL KOSY 2 1 Eletronis researh institute, Cairo, Egypt 2 Faulty of Engineering, Cairo University,
More informationIterative Image Formation using Fast (Re/Back)-projection for Spotlight-mode SAR
Iterative Image Formation using Fast (Re/Back)-projection for Spotlight-mode SAR Shaun I. Kelly, Gabriel Rilling, Mike Davies and Bernard Mulgrew Institute for Digital Communications The University of
More informationDefect Detection and Classification in Ceramic Plates Using Machine Vision and Naïve Bayes Classifier for Computer Aided Manufacturing
Defet Detetion and Classifiation in Cerami Plates Using Mahine Vision and Naïve Bayes Classifier for Computer Aided Manufaturing 1 Harpreet Singh, 2 Kulwinderpal Singh, 1 Researh Student, 2 Assistant Professor,
More informationA Fast Decimation-in-image Back-projection Algorithm for SAR
A Fast Decimation-in-image Back-projection Algorithm for SAR Shaun I. Kelly and Mike E. Davies Institute for Digital Communications The University of Edinburgh email: {Shaun.Kelly, Mike.Davies}@ed.ac.uk
More informationFuzzy Meta Node Fuzzy Metagraph and its Cluster Analysis
Journal of Computer Siene 4 (): 9-97, 008 ISSN 549-3636 008 Siene Publiations Fuzzy Meta Node Fuzzy Metagraph and its Cluster Analysis Deepti Gaur, Aditya Shastri and Ranjit Biswas Department of Computer
More informationFlow Demands Oriented Node Placement in Multi-Hop Wireless Networks
Flow Demands Oriented Node Plaement in Multi-Hop Wireless Networks Zimu Yuan Institute of Computing Tehnology, CAS, China {zimu.yuan}@gmail.om arxiv:153.8396v1 [s.ni] 29 Mar 215 Abstrat In multi-hop wireless
More informationMultiple-Criteria Decision Analysis: A Novel Rank Aggregation Method
3537 Multiple-Criteria Deision Analysis: A Novel Rank Aggregation Method Derya Yiltas-Kaplan Department of Computer Engineering, Istanbul University, 34320, Avilar, Istanbul, Turkey Email: dyiltas@ istanbul.edu.tr
More informationThe Implementation of RRTs for a Remote-Controlled Mobile Robot
ICCAS5 June -5, KINEX, Gyeonggi-Do, Korea he Implementation of RRs for a Remote-Controlled Mobile Robot Chi-Won Roh*, Woo-Sub Lee **, Sung-Chul Kang *** and Kwang-Won Lee **** * Intelligent Robotis Researh
More informationDubins Path Planning of Multiple UAVs for Tracking Contaminant Cloud
Proeedings of the 17th World Congress The International Federation of Automati Control Dubins Path Planning of Multiple UAVs for Traking Contaminant Cloud S. Subhan, B.A. White, A. Tsourdos M. Shanmugavel,
More informationCluster Centric Fuzzy Modeling
10.1109/TFUZZ.014.300134, IEEE Transations on Fuzzy Systems TFS-013-0379.R1 1 Cluster Centri Fuzzy Modeling Witold Pedryz, Fellow, IEEE, and Hesam Izakian, Student Member, IEEE Abstrat In this study, we
More informationLARGE SCENE SAR IMAGE FORMATION
LARGE SCENE SAR IMAGE FORMATION A dissertation submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy by LEROY A. GORHAM A.S., Hudson Valley Community College, 1995
More informationA Partial Sorting Algorithm in Multi-Hop Wireless Sensor Networks
A Partial Sorting Algorithm in Multi-Hop Wireless Sensor Networks Abouberine Ould Cheikhna Department of Computer Siene University of Piardie Jules Verne 80039 Amiens Frane Ould.heikhna.abouberine @u-piardie.fr
More informationBSPLND, A B-Spline N-Dimensional Package for Scattered Data Interpolation
BSPLND, A B-Spline N-Dimensional Pakage for Sattered Data Interpolation Mihael P. Weis Traker Business Systems 85 Terminal Drive, Suite Rihland, WA 995-59-946-544 mike@vidian.net Robert R. Lewis Washington
More informationAn Optimized Approach on Applying Genetic Algorithm to Adaptive Cluster Validity Index
IJCSES International Journal of Computer Sienes and Engineering Systems, ol., No.4, Otober 2007 CSES International 2007 ISSN 0973-4406 253 An Optimized Approah on Applying Geneti Algorithm to Adaptive
More informationThe Minimum Redundancy Maximum Relevance Approach to Building Sparse Support Vector Machines
The Minimum Redundany Maximum Relevane Approah to Building Sparse Support Vetor Mahines Xiaoxing Yang, Ke Tang, and Xin Yao, Nature Inspired Computation and Appliations Laboratory (NICAL), Shool of Computer
More informationCalculation of typical running time of a branch-and-bound algorithm for the vertex-cover problem
Calulation of typial running time of a branh-and-bound algorithm for the vertex-over problem Joni Pajarinen, Joni.Pajarinen@iki.fi Otober 21, 2007 1 Introdution The vertex-over problem is one of a olletion
More informationAn Interactive-Voting Based Map Matching Algorithm
Eleventh International Conferene on Mobile Data Management An Interative-Voting Based Map Mathing Algorithm Jing Yuan* University of Siene and Tehnology of China Hefei, China yuanjing@mail.ust.edu.n Yu
More informationFinding the Equation of a Straight Line
Finding the Equation of a Straight Line You should have, before now, ome aross the equation of a straight line, perhaps at shool. Engineers use this equation to help determine how one quantity is related
More informationGradient based progressive probabilistic Hough transform
Gradient based progressive probabilisti Hough transform C.Galambos, J.Kittler and J.Matas Abstrat: The authors look at the benefits of exploiting gradient information to enhane the progressive probabilisti
More informationDiscrete Contours in Multiple Views: Approximation and Recognition
Disrete Contours in Multiple Views: Approximation and Reognition M. Pawan Kumar, Saurabh Goyal, Sujit Kuthirummal, C. V. Jawahar, and P. J. Narayanan Centre for Visual Information ehnology International
More informationWe don t need no generation - a practical approach to sliding window RLNC
We don t need no generation - a pratial approah to sliding window RLNC Simon Wunderlih, Frank Gabriel, Sreekrishna Pandi, Frank H.P. Fitzek Deutshe Telekom Chair of Communiation Networks, TU Dresden, Dresden,
More informationVolume 3, Issue 9, September 2013 International Journal of Advanced Research in Computer Science and Software Engineering
Volume 3, Issue 9, September 2013 ISSN: 2277 128X International Journal of Advaned Researh in Computer Siene and Software Engineering Researh Paper Available online at: www.ijarsse.om A New-Fangled Algorithm
More informationPreliminary investigation of multi-wavelet denoising in partial discharge detection
Preliminary investigation of multi-wavelet denoising in partial disharge detetion QIAN YONG Department of Eletrial Engineering Shanghai Jiaotong University 8#, Donghuan Road, Minhang distrit, Shanghai
More informationComprehensive Quantification of Signal-to-Noise Ratio and g-factor for Image-Based and k-space-based Parallel Imaging Reconstructions
Magneti Resonane in Mediine 60:895 907 (2008) Comprehensive Quantifiation of Signal-to-Noise Ratio and g-fator for Image-Based and k-spae-based Parallel Imaging Reonstrutions Philip M. Robson, 1 * Aaron
More informationA RAY TRACING SIMULATION OF SOUND DIFFRACTION BASED ON ANALYTIC SECONDARY SOURCE MODEL
19th European Signal Proessing Conferene (EUSIPCO 211) Barelona, Spain, August 29 - September 2, 211 A RAY TRACING SIMULATION OF SOUND DIFFRACTION BASED ON ANALYTIC SECONDARY SOURCE MODEL Masashi Okada,
More informationthe data. Structured Principal Component Analysis (SPCA)
Strutured Prinipal Component Analysis Kristin M. Branson and Sameer Agarwal Department of Computer Siene and Engineering University of California, San Diego La Jolla, CA 9193-114 Abstrat Many tasks involving
More informationNew Fuzzy Object Segmentation Algorithm for Video Sequences *
JOURNAL OF INFORMATION SCIENCE AND ENGINEERING 24, 521-537 (2008) New Fuzzy Obet Segmentation Algorithm for Video Sequenes * KUO-LIANG CHUNG, SHIH-WEI YU, HSUEH-JU YEH, YONG-HUAI HUANG AND TA-JEN YAO Department
More informationMigration based SAR imaging for ground penetrating radar systems
Migration based SAR imaging for ground penetrating radar systems K. Gu, G. Wang and J. Li Abstrat: The authors onsider migration based syntheti aperture radar (SAR imaging of surfae or shallowly buried
More informationA Unified Subdivision Scheme for Polygonal Modeling
EUROGRAPHICS 2 / A. Chalmers and T.-M. Rhyne (Guest Editors) Volume 2 (2), Number 3 A Unified Subdivision Sheme for Polygonal Modeling Jérôme Maillot Jos Stam Alias Wavefront Alias Wavefront 2 King St.
More information