Wind-vector estimation for RADARSAT-1 SAR images: Validation of wind-direction estimates based upon geometry diversity

In this letter, a new wind-vector algorithm is presented that uses radar backscatter sigma(0) measurements at two adjacent subscenes of RADARSAT-1 synthetic aperture radar (SAR) images, with each subscene having slightly different geometry. Resultant wind vectors are validated using in situ buoy measurements and compared with wind vectors determined from a hybrid wind-retrieval model using wind directions determined by spectral analysis of wind-induced image streaks and observed by colocated QuikSCAT measurements. The hybrid wind-retrieval model consists of CMOD-IFR2 [applicable to C-band vertical-vertical (W) polarization] and a C-band copolarization ratio according to Kirchhoff scattering. The new algorithm displays improved skill in wind-vector estimation for RADARSAT-1 SAR data when compared to conventional wind-retrieval methodology. In addition, unlike conventional methods, the present method is applicable to RADARSAT-1 images both with and without visible streaks. However, this method requires ancillary data such as buoy measurements to resolve the ambiguity in retrieved wind direction.

合成孔径雷达提取海面风、浪信息研究

传统海洋测量方法例如岸基观测站、船只和浮标等方式只能对海面进行单点观测，存在很多不足之处：观测点少，只能获得有限的点的资料，无法得到大范围的海面信息；费用较高，无法大范围的密集观测；传统方法受天气影响较大，无法长时间的连续观测；星载合成孔径雷达能全天候、全天时、高分辨率对海面成像，能实现多波段、多极化、多视角得观测海面，提供大范围、高精度的实时动态海面信息。同时，这种全天候、全天时和高分辨率观测海洋的优势是可见光和红外传感器所没有的。总之研究如何从SAR影像中有效地获取海面信息具有重要的科学和实用意义。 自从1978年Seasat卫星发射以来SAR图像就广泛的应用于海洋要素的反演如海面风场、波高、平均周期。其中对海面风场的反演研究的最多，一般的方法是首先对SAR图像进行快速傅立叶变换得到SAR图像谱，通过图像谱的峰值信息能够得到具有180°模糊的风向，利用SAR图像的条纹或者浮标等外部信息来消除180°模糊确定风向，其次把风向和图像的正交后向散射截面数值(NRCS)带入经验函数CMOD4通过迭代计算得到风速。这种方法得到的海面风场需要外部信息的辅助，不利于大范围的业务化的反演海面风场。因此本文试图直接通过SAR图像来反演海面风场而不利用其他外部信息的帮助。2002年3月1日发射的ENVISAT卫星所获得的ASAR图像具有不同于ERS系列的SAR图像的新特点：不同入射角的情况下获得ASAR图片、双极化数据以及可以获得更宽的带幅的图像。双极化数据可以同时地提供同极化和交错极化的影像，两种极化的数据能够增加分辨目标的能力和提供目标更多的信息，这方面的优势也使得利用双极化的数据来消除SAR图像反演风场中的180º模糊问题提供了可能。本文推导了利用双极化的ASAR图像反演海面风场的新函数，该函数是在同极化和交错极化函数的基础上推导出来的，有效的消除了只用一副ASAR图片反演风场所固有的180º模糊现象。风速和风向的反演结果与Quikscat数据之间的均方根误差分别为0.53 m/s和2.21º。该方法与传统方法的比较可以看出新方法与浮标数据以及Quikscat数据符合的更好。 利用SAR图像来反演海浪的波高也是一个研究的比较多的领域。比较传统的方法是由SAR图像得到的海浪谱计算出波高，目前国际上比较流行的SAR图像反演海浪谱的模式有两种：一种是Hasselmann（1991）提出的在Max－Planck Institute（MPI）发展起来的方法，Hasselmann（1996）进行了改进和完善。另一种是Mastenbroek和de Valk（2000）提出的半参数化反演方法。这两种方法的最主要的缺点就是需要引进外部信息(WAM模式结果或者散射计的风的信息)来消除图像所固有的180°方向模糊问题。而本文应用的经验函数方法(CWAVE)不需要引进任何外部信息而直接得到海浪的重要参数－有效波高(Hs)，输入该经验函数的参数主要有：波模式图像的雷达截面、图像方差以及由SAR谱得到的20个参数。这些参数的选择是通过逐步回归方法进行筛选的。CWAVE经验函数的系数是通过6000幅全球分布的ERS-2波模式的图像谱拟合同时同地点的WAM模式结果得到的。利用CWAVE经验函数反演了1998年9月到2000年11月两年多的全球接近一百万的ERS-2 SAR图像的有效波高，利用NOAA浮标数据对反演结果进行了验证，他们之间的相关系数为0.83, 均方根误差为0.61m，偏差为0.02 m；反演结果也与欧洲中长期预报中心(ECMWF)的ERA-40有效波高资料和高度计资料进行了比较，结果表明该方法是通过SAR图像反演Hs的一种有效方法。 SAR的反演结果也应用于有效波高的非线性统计分布的研究。结合动力学和随机统计学原理推导了海面高度、波高和有效波高的非线性统计分布，同时利用了NOAA浮标数据、华师大的波高数据以及SAR的有效波高数据分别对推导出的非线性统计分布函数进行了检验。 SAR的结果也用来反演新的风浪成长关系： ，与已有的风浪成长关系的比较表明该成长关系与已有的结果比较一致。

干涉合成孔径雷达海浪遥感理论与应用研究

干涉合成孔径雷达 (Interferometric Synthetic Aperture Radar，InSAR) 是以合成孔径雷达复数据提取的相位为信息源获取地表三维信息和海表散射体运动信息的新型微波成像雷达。 InSAR通过两幅天线同时观测 (单轨模式)，或两次近平行的观测 (重复轨道模式)，获取地面或海面同一景观的复图像对。20世纪90年代以来，InSAR陆地和海洋研究成为微波遥感的热点，广泛应用于地表变形监测、南极冰流测量、地面或海面慢速运动目标检测等领域。 近年来，国际上逐渐应用机载顺轨或交轨干涉合成孔径雷达进行海表面流速测量以及海面波成像机制研究。相对于传统单天线合成孔径雷达 (Synthetic Aperture Radar, SAR)，双天线干涉合成孔径雷达 (InSAR) 测量海表面波有着独特的优势： (1) InSAR复图像的相位近似正比于海面散射体的径向速度，这种内在的成像机制提供了直接测量海表面动态运动的机会。 (2) 真实孔径雷达调制传递函数几乎对InSAR相位图像没有影响，而对传统SAR图像影响较大。 基于干涉合成孔径雷达测量海浪的优势，本文做了一些干涉合成孔径雷达海浪遥感理论与应用研究工作，主要内容大致可归纳如下： (1)基于新建立的顺轨干涉合成孔径雷达 (Along-Track Interferometric Synthetic Aperture Radar，ATI-SAR) 相位谱与海浪谱之间的非线性映射关系，通过数值模拟研究了不同雷达参数和海况参数对应的ATI-SAR相位谱。数值模拟结果表明：距离速度比率、雷达入射角、天线间距和有效波高和波长比率是影响ATI-SAR海浪成像的重要因素。进一步，利用机载X波段水平极化相位图像和机载C波段水平极化相位图像谱结合方向波骑士浮标测量的海浪方向谱验证了ATI-SAR相位谱与海浪谱之间的非线性映射关系。结果显示用前向映射关系计算的相位谱与实际观测的相位谱较为一致，二者相关系数总体大于0.6，而且对成像非线性不敏感. (2)建立了包含海表面高度和速度聚束的交轨干涉合成孔径雷达 (Across-Track Interferometric Synthetic Aperture Radar，XTI-SAR) 涌浪干涉相位模型，得到了涌浪成像的解析表达式，进一步研究了XTI-SAR沿方位向传播的涌浪成像机制。定义二次谐波振幅与基波振幅比率来表征成像非线性，通过比较XTI-SAR和ATI-SAR相位的二阶调和分量，分析不同海况和干涉SAR参数情况下的数值模拟，结果表明：当速度聚束弱时，XTI-SAR相位比ATI-SAR相位具有较强的非线性，ATI-SAR比XTI-SAR更适合测量海浪。当速度聚束强时，XTI-SAR相位比ATI-SAR相位具有较弱的非线性，XTI-SAR比ATI-SAR更适合测量海浪。 (3)基于包含海表面高度和速度聚束的交轨干涉合成孔径雷达 (XTI-SAR) 涌浪干涉相位模型，结合多维高斯变量的特征函数方法建立了新的XTI-SAR相位谱与海浪谱非线性积分变换。新积分变换不同于Bao (1999) 建立的积分变换，两者形式上区别在于新积分变换中包含了长波径向轨道速度一阶倒数项。数值模拟显示：通常情况下，长波径向轨道速度一阶倒数项不能忽略。进一步，我们针对不同雷达参数和海况结合新非线积分变换对XTI-SAR海浪成像进行了数值模拟，结果表明：同顺轨干涉合成孔径雷达 (ATI-SAR) 海浪成像一致，距离速度比率 和有效波高与波长比率 是影响XTI-SAR海浪成像的重要因子。 (4)基于新的ATI-SAR相位谱与海浪谱之间的非线性映射关系，发展了利用ATI-SAR相位图像反演海浪方向谱的参数化反演模式，并由此得到海浪波长、波向和有效波高。反演结果与现场浮标观测结果比较一致。相对于其它反演模式，参数化反演模式的优点在于：(1) 不需要任何附加信息如初猜海浪谱、散射计提供的风速风向等信息。 (2) 不需要对相位图像进行辐射定标，可以由反演的海浪谱直接计算有效波高。(3) 反演结束后还可以得到成像区域的局地风速信息。因此，参数化反演模式可以实现风、浪信息的联合反演。

波浪破碎散射模式的研究及应用

海面反射、散射或自发辐射的各个波段的电磁波携带着海表面温度、海平面高度、海表面粗糙度以及海水所含各种物质浓度的信息。波浪破碎是海洋的一种重要现象，波浪破碎支撑着海－气之间的热通量和气体通量的交换，这些通量对于天气和气候都有深远影响。同时，波浪破碎使得海表的粗糙度增加，这对于海表散射电磁波有极大的影响，因此在处理雷达散射问题时，考虑波浪破碎的影响十分必要。 本文分析了雷达散射的两种模式，并阐述了波浪破碎在海表电磁波散射过程中的重要性。论文的创新点在于利用Kudryavtsev03模式，反演得到了波浪破碎率与风速的两种新关系式，同时得到了较高精度的风速反演的隐含模式。针对波浪破碎微波散射模式问题，论文主要研究结果如下： 1.在Kudryavtsev03模式的基础上，以黄兴忠的白冠覆盖率公式替代破碎率，利用ERS-1散射计数据，对Bragg散射模式和Kudryavtsev03模式的结果进行了比较。结果表明：在相同条件下，考虑波浪破碎的Kudryavtsev03模式明显优于纯Bragg散射模式。 2 对Kudryavtsev03模式进行了修正。通过与实测ERS-1散射计数据比较，表明修正后，入射角在18-35°范围内，模式的结果具有较高的精度。 3 利用修正后的模式及散射计数据，在入射角为30-35°范围内，反演破碎率，并对反演结果进行了拟合，得到了两个新的破碎率与风速的关系式。 4 利用修正后的模式及新的破碎率公式，在入射角为30-35°范围内，反演风速，具有较高精度。 5 利用反演的风速，根据经验关系式，推算了各种海洋物理量。

交通监控系统中运动目标分类和跟踪研究

文章讲述了交通监控系统中应用视频图像流来跟踪运动目标并对目标进行分类的具体过程和原则.基于目标检测提出了双差分的目标检测算法,目标分类应用到了连续时间限制和最大可能性估计的原则,目标跟踪则结合检测到的运动目标图像和当前模板进行相关匹配.实验结果表明,该过程能够很好地探测和分类目标,去除背景信息的干扰,并能够在运动目标部分被遮挡、外观改变和运动停止等情况下连续地跟踪目标.

车辆检测毫米波雷达信号处理技术研究

道路交通信息采集在智能交通系统中，担负着提供准确可靠的信息源以使整个智能交通系统得以顺利准确运行的重任。毫米波雷达交通检测器作为获取路面信息的一种方式，因其特有的优势具有广阔的应用前景。本论文是围绕一个基于毫米波雷达的车辆检测系统的研究课题进行的。 本文的主要工作及创新之处在于： (1) 结合毫米波车辆检测雷达的应用背景，首先阐述毫米波车辆检测雷达在智能交通系统中的应用，介绍了智能交通系统中路面信息采集各种不同方式的优缺点，同时给出采用毫米波车辆检测雷达系统进行交通信息采集的优越性和必要性，说明了本文的选题背景及主要研究内容。 (2) 针对毫米波车辆检测雷达的具体应用需要，在提出车辆检测雷达的工作原理，阐述毫米波车辆检测雷达的系统构成和整体框架结构之后，分别从系统硬件设计与实现、数字脉冲压缩原理和基于一维距离像的车辆目标识别三个方面进行阐述。根据实际应用提出了车辆检测雷达系统的性能技术参数要求，确定了系统实际的工作频段和工作体制，分析了线性调频系统的组成原理，对系统主要参数进行了定量计算和误差估计。 (3) 分析了几种简单形状目标的散射特性，对典型车辆的雷达截面积进行了粗略的工程估算。对宽带雷达杂波特性进行了分析，采用幅度为韦布尔分布、功率谱为高斯分布的杂波模型I模拟地杂波，幅度为瑞利分布、功率谱为高斯分布的杂波模型II来模拟气象杂波。 (4) 介绍了一维距离像的基本概念和获得一维距离像的方法，深入分析了一维距离像的基本特性。通过对基于位置信息的全程相关检测方法和基于能量积累的几种滑窗检测方法的比较，分析了它们在宽带雷达信号检测中的性能特点和应用条件。在此基础上，提出了一种适合宽带雷达的检测方法，即相邻脉冲重复周期相关检测(IPCD)方法。当采用多脉冲积累后，其检测能力比能量相加脉间非相干积累EAD检测平均提高2dB。

Millstone Hill电离层F层等效风场的气候学研究

Neutral winds and electric fields in the ionospheric F layer play important roles in the variations of the ionosphere, and also affect the thermospheric circulation via the close coupling between the ionosphere and the thermosphere. By now, the neutral winds and electric drifts are generally observed with ground-based Fabry-Perot interferometers (FPI) and incoherent scatter radars (ISR), rockets, and satellite-borne instrument. Based on the servo theory, the ionospheric equivalent winds, which include the information of both the neutral winds and electric fields, can be derived from these characteristic parameters observed by ionosondes. This indirect derivation has potential values in climatological researches and space weather forecast. With the data set of the incoherent scatter radar observations at Millstone Hill, USA, from 1976 to 2006, we statistically analyzed the climatological variations of the vertical component of the equivalent winds (VEWs) over Millstone Hill, which are derived from the ionospheric key parameters (the peak electron number density and peak height of the F2 layer, NmF2 and hmF2) on the basis of the servo theory, Liu's method, and measurements from the ion line-of-sight velocity as well. The main results of this analysis are summarized as follows: (1) The values of VEWs over Millstone Hill during nighttime are stronger than in the daytime, and the upward drift dominates most of the day. In 1993, Hagan found that the component of the neutral winds in the magnetic meridion in daytime is weaker than during nighttime under both solar maximum and minimum conditions; he also found that the equatorward winds dominate most of the day. Both results suggest that the thermosphere in Millstone Hill is modulated by the aurorally driven high-latitude circulation cell; that is, during geomagnetic quiet periods, the average auroral activity is strong enough to drive thermospheric circulation equatorward for most of the day at Millstone Hill. Moreover, since ion drag is the strongest during daytime when F region densities are enhanced by photoionization, the wind speeds are smaller during the daytime than in the nighttime. (2) There is equinoctial symmetry in VEWs at Millstone Hill. The amplitudes and phases of VEWs in spring are quite similar to those in autumn. In contrast, the nighttime upward drift in winter is weaker than in summer and the difference becomes more significant with increasing solar activity. This solstice asymmetry indicates that, the aurorally driven circulation in the northern hemisphere at Millstone Hill latitude is weaker in winter due to arctic darkness, because the subsolar point is in the southern hemisphere. (3) The comparison of the VEWs derived from three methods, i.e., the servo theory, Liu's method, and the ISR ion line-of-sight velocity measurements, indicates that the amplitudes and main phase tendencies of these VEWs accord well with each other during nighttime hours. However, the case in the daytime is relatively worse. This daytime discrepancy can be explained in terms of the effects of photochemical processes and the choices of the servo constants. A larger servo constant gives a stronger plasma drift in daytime. Therefore, this study tells how important to choose a suitable constant for deriving VEWs at Millstone Hill.

Mapeamento da aptidão agrícola das terras do Estado do Maranhão.

2007

Aptidão agrícola das terras do nordeste do Estado de São Paulo.

Solos do nordeste do Estado de São Paulo; Sistema de Avaliação da Aptidão Agrícola das Terras; Mapeamento da Aptidão Agrícola das Terras do Nordeste do Estado de São Paulo.

Levantamento de reconhecimento de baixa intensidade dos solos e aptidao agricola das terras de parte da Regiao Geoeconomica de Brasilia.

O levantamento visou atender o CPAC que em convenio com a SUDECO se propos a estudar os recursos naturais da area Programa Especial da Regiao Geoeconomica de Brasilia, de maneira que sirvam de suporte para indicacao de sistemas alternativos de producao agropecuaria condizente com o potencial de uso das terras. O SNLCS se propos a executar o mapeamento de solos e a avaliacao da aptidao agricola das terras, a nivel de reconhecimento de baixa intensidade de uma area prioritaria indicada pelo CPAC de aproximadamente 68.000km2. Este estudo propiciou o conhecimento da potencialidade da area e a obtencao de indicadores seguros e economicos para a escolha de locais apropriados a merecerem maiores atencoes e recursos. A metodologia empregada foi basicamente a utilizadapelo SNLCS. Utilizou-se como material basico de mapeamento, imagens de radar na escla 1:250.000, plantas alti-planimetricas 1:100.000 e 1:250.000. A escala final dos mapas e de 1:500.000. Foram mapeados os seguintes solos : Latossolo Vermelho-Amarelo, Latossolo Vermelho-Escuro, Latossolo Roxo, Latossolo Petroplintico, Podzolico Vermelho-Amarelo, Podzolico Vermelho-Escuro, Cambissolo, Plintossolo, Podzol, Solos Aluviais, Solos Hidromorficos e Solos Litolicos. Quando necessario, foram subdivididos de acordo com a saturacao de bases, tipo de horizonte A, textura, ocorrencia de cascalhos e calhaus, drenagem e fasados de acordo com o relevo e vegetacao. A interpretacao da aptidao agricola das terras indicou que aproximadamente 36% das terras pertencem a classe ...

Monitoramento hidrológico em voçoroca submetida a práticas mecânicas e vegetativas: aeroporto internacional do Rio de Janeiro Galeão/Antônio Carlos Jobim.

O desenvolvimento desta dissertação está vinculado ao projeto de RAD (recuperação de áreas degradadas) do Aeroporto Internacional do Rio de Janeiro, sendo seu objeto o processo de recuperação da voçoroca do Morro do Radar. Através da realização de uma investigação do comportamento hidrológico desta voçoroca após a implantação do projeto, buscou-se fornecer informações que possam constituir um subsídio à avaliação dos efeitos das práticas adotadas A metodologia utilizada compreendeu seis etapas: avaliação das variações da topografia local; monitoramento da precipitação; monitoramento da intercepção das chuvas promovida pelas copas arbóreas e pela serrapilheira; avaliação das propriedades fisico-hídricas dos solos em diferentes profundidades; avaliação dos potenciais matriciais e das cargas totais nos solos nas mesmas profundidades, através de gráficos e construção de mapas de equipotenciais; avaliação da relação entre os condicionantes locais e o desempenho das práticas de RAD. Os valores de intercepção pelas copas arbóreas e pela serrapilheira revelaram uma proteção eficaz do solo contra o impacto direto das águas pluviais. Os resultados de condutividade hidráulica mostraram uma boa relação com o comportamento das cargas de pressão e totais. Com relação à expansão e contração da zona de saturação, observou-se que a umidade antecedente está condicionando de forma decisiva o processo de expansão. Foi encontrada também ao longo de todas as discussões grande influência da topografia sobre os processos de recuperação, o que indica a importância do planejamento adequado da reconformação do talude. As áreas com menores declives apresentaram uma recuperação proporcionalmente mais rápida, o que se relacionou tanto ao estabelecimento da cobertura vegetal e à incorporação de serrapilheira, quanto à própria dinâmica hidro-erosiva inerente.

Scalable Multi-Relational Association Mining

Clare, A., Williams, H. E. and Lester, N. M. (2004) Scalable Multi-Relational Association Mining. In proceedings of the 4th International Conference on Data Mining ICDM '04.

Interaction between coronal mass ejections and the solar wind

Jones, R. A.; Breen, A. R.; Fallows, R. A.; Canals, A.; Bisi, M. M.; Lawrence, G. (2007). Interaction between coronal mass ejections and the solar wind, Journal of Geophysical Research, 112, Issue A8 RAE2008

Estimating 3D Hand Pose from a Cluttered Image

A method is proposed that can generate a ranked list of plausible three-dimensional hand configurations that best match an input image. Hand pose estimation is formulated as an image database indexing problem, where the closest matches for an input hand image are retrieved from a large database of synthetic hand images. In contrast to previous approaches, the system can function in the presence of clutter, thanks to two novel clutter-tolerant indexing methods. First, a computationally efficient approximation of the image-to-model chamfer distance is obtained by embedding binary edge images into a high-dimensional Euclide an space. Second, a general-purpose, probabilistic line matching method identifies those line segment correspondences between model and input images that are the least likely to have occurred by chance. The performance of this clutter-tolerant approach is demonstrated in quantitative experiments with hundreds of real hand images.

Segmenting Foreground Objects from a Dynamic Textured Background via a Robust Kalman Filter

The algorithm presented in this paper aims to segment the foreground objects in video (e.g., people) given time-varying, textured backgrounds. Examples of time-varying backgrounds include waves on water, clouds moving, trees waving in the wind, automobile traffic, moving crowds, escalators, etc. We have developed a novel foreground-background segmentation algorithm that explicitly accounts for the non-stationary nature and clutter-like appearance of many dynamic textures. The dynamic texture is modeled by an Autoregressive Moving Average Model (ARMA). A robust Kalman filter algorithm iteratively estimates the intrinsic appearance of the dynamic texture, as well as the regions of the foreground objects. Preliminary experiments with this method have demonstrated promising results.