[求助]各位处理过IPIX雷达数据的请进！

数据可以读到matlab里！现在的问题是按说明《how to process the IPIX radar files in MATLAB》前3步都没有问题！现在的问题是在第4步Loading the data into the MATLAB workspace出现问题，我已经将网站ipixload.m文件下载放到F:\netcdf\install files（ipixinfo.m文件也放在这个路径里且好用），当我按说明书执行命令[I,Q,meanIQ,stdIQ,inbal] = ipixload(nc,pol,rangebin,mode)，matlab提示问题??? Undefined function or variable 'pol'.matlab当前工作路径也为：F:\netcdf\install files！这一步不成功也不能得到海杂波处理结果？这问题可能是出在哪里呢？谢谢各位指点！！附件：是按说明书《how to process the IPIX radar files in MATLAB》的第四步4. Loading the data into the MATLAB workspaceThe actual IPIX radar data is stored in the variable adc_data. This array has three or four dimensions, depending on the number of transmit polarizations used, as indicated by the global attribute TX_polarization. This attribute can have the values 'H' (all pulses have horizontal polarization), 'V' (all pulses have vertical polarization) or 'A' (alternating). The dimensions of the variable adc_data are: 1. nsweep - the number of sweeps stored in the array. Each sweep contains two pulses if the transmit polarization is alternating, and one otherwise. The number of sweeps measured per second can thus be derived from the variable PRF, the Pulse Repetition Frequency. 2. ntxpol- the number of different transmit polarizations used, as indicated by the global attribute TX_polarization. This dimension only exists if the transmit polarization alternates (TX_polarization = 'A'). 3. nrange - the number of range bins measured. The variable range contains the starting range of each range bin. 4. nadc - the number of analog-to-digital-converted numbers stored in the array. This number is always 4, because the radar receiver measures, for each radar return, the horizontally and vertically polarized components. These two components are further divided in their in-phase and quadrature components. Download the m-file ipixload.m. The m-file ipixload.m pre-processes the data, and loads all sweeps of a particular polarization and rangebin into MATLAB. The syntax is: [I,Q,meanIQ,stdIQ,inbal] = ipixload(nc,pol,rangebin,mode)The inputs are: • nc - pointer to netCDF file • pol - polarization specification, with the following four possible values: 'hh' (horizontal transmit, horizontal receive) 'hv' (horizontal transmit, vertical receive) 'vv' (vertical transmit, vertical receive) 'vh' (vertical transmit, horizontal receive) • rangebin - range bin to load (between 1 and nrange) • mode - specifies the type of pre-processing to be used: 'raw' (no pre-processing) 'auto' (automatic pre-processing, see below) 'dartmouth' (pre-processing for dartmouth files containing land, see below) The outputs are: • I - nsweep-by-1 matrix, containing the in-phase component after pre-processing • Q - nsweep-by-1 matrix, containing the quadrature component after pre-processing • meanIQ - 1-by-2 matrix, containing the mean of I and Q before pre-processing • stdIQ - 1-by-2 matrix, containing the standard deviation of I and Q before pre-processing • inbal - 1-by-1 matrix, containing the phase-inbalance (in degrees) between I and Q before pre-processing - Automatic Pre-processing -The automatic preprocessing performed by ipixload.m consists of two steps: 1. Remove mean and standard deviation from the I and Q channels separately. 2. Remove the phase inbalance. A phase inbalance occurs when the I and Q channels are not orthogonal to each other, due to hardware imperfections. Let A and φ be the amplitude and phase of the received radar echo, and let β be the phase inbalance in radians, equal to inbal*π/180. Ideally, the I and Q channels measure I=A cos(φ) and Q=A sin(φ). With a phase inbalance of β radians, the I and Q channels are I=A cos(X) and Q=A sin(X β). Given a long data record and assuming that the phase is uniformly distributed between 0 and 2π, the inbalance β is estimated by sin(β)=2(I,Q)/A2, where (I,Q) is the inner product of I and Q. The phase inbalance is removed by rotating the I component: Irot=[I-Q*sin(β)]/sqrt(1-sin(β)2). Derivation: E[ A cos(φ) A sin(φ β) ] = A2 E[ 1/2 sin(β) 1/2 sin(2φ β) ] = A2/2 sin(β), since the expectation of sin(2 φ β) is zero for uniformly distributed φ. - Files containing land clutter -The phase inbalance removal procedure assumes that the phase is uniformly distributed between 0 and 2π. This assumption is true if the distance between the radar antenna and the major scatterers is randomly varying with time, with a variance much larger than the radar wavelength of 3 cm. For ocean wave scatterers this certainly holds, but for land scatterers, the phase may be nearly constant. Therefore, to apply the phase-inbalance removal to scanning radar files which contain land clutter, the phase inbalance is estimated using only that part of the data that is known not .. 5p}ri,Y<  
L/Ytkag  

未注册仅能浏览部分内容，查看全部内容及附件请先 登录 或 注册