利用温度链对南海内波振幅的观测

内波的实测资料一般为走航或者浮标布站观测得到的一个关于时间和空间的大面资料,将此资料进行滤波后,进行简单的数值处理,即可得到内波的振幅。但是这种资料较少,布控站位难度较大,成本也较高,是一种不易获取的资料。本文利用的温度链资料是一种获取简单,成本也较低的实测资料。利用温度链资料基于内波带动水质点上下波动的物理过程,进行内波振幅的计算,将会有较大的适用性。     

基于特征选取及广义傅立叶分形的SAR图像海洋溢油检测算法

针对海上溢油SAR图像中油膜与类油膜的识别问题,提出了一种结合傅立叶分形与特征提取的检测算法。由于分形特征可以具有无穷多的细节,并在不同的研究尺度存在自仿射特性。这与油膜及类油膜表面的几何形貌特征非常吻合。该算法通过计算样本的傅立叶分形特征,组成油膜与类油膜的特征空间。然后,应用基于差分进化的特征选取方法将利于分类的重要特征值筛选出来。再利用重要特征值对原有样本进行分类。实验表明,经特征选取的分形特征向量能够以100%的准确率将两类样本区分开。该算法在选取重要特征的同时实现了对高维特征空间降维的目的,该思想可以应用于其他的基于高维特征的识别系统中,具有普遍的适用性。     

高分三号SAR影像双阈值变化检测

双阈值合成孔径雷达SAR(Synthetic Aperture Radar)变化检测算法具有在发现变化区域的同时还能确定地表发生后向散射变化类型的优点。针对广义高斯双阈值最小误差法D-GKIT(Dual Generalized Kittler and Illingworth Thresholding)在进行阈值选取时直方图中不同类别像素灰度级重叠严重时,分割结果容易在尖峰单侧选取出双阈值而导致无法正确分割差异图的问题,本文提出一种结合归一化最大类间方差和广义高斯最小误差法GKIT(Generalized Kittler and Illingworth Thresholding)的双阈值SAR变化检测方法。首先,提出以归一化最大类间方差值作为灰度级重叠程度的判别参数,确定阈值的选取顺序及两个候选区间;然后,利用GKIT在候选区间内进行分割,获取单侧阈值及非变化类拟合函数;最后,提出利用非变化类拟合函数更新后的直方图作为另一侧阈值选取基础进行分割,得到对应分割阈值。以宁波地区高分三号(GF-3)SAR卫星影像作为试验研究数据,结果表明:本文方法能较好地解决灰度级重叠时D-GKIT无法进行正确分割的问题,具有良好的变化检测效果和更强的鲁棒性且达到了利用研究区数据验证利用GF-3号SAR卫星影像进行变化检测研究可行性的目的。     

利用SIFT与粗差探测进行SAR影像配准

SAR影像特殊的成像机理使得SIFT在SAR影像配准中的错误率较高,加入相干系数的辅助只能在一定程度上削弱SIFT的错误配准。由于仅考虑灰度的配准策略对精度的提高有限,本文将摄影测量中粗差探测和剔除的方法与SIFT算法相结合。在二次多项式平差的过程中,将错误同名点视为粗差,利用粗差剔除的方法提高配准精度。实验证明了此方法的可行性。     

北极海域SAR海浪方向谱反演及其与中法海洋卫星CFOSAT/SWIM数据的比较

The interaction between ocean waves and sea ice in the Arctic ocean has received significant attention. However, the study on this issue is significantly limited due to the lack of observation data. Synthetic Aperture Radar (SAR) plays an important role in the research of ocean wave and sea-ice interaction because of its unique capability of imaging sea surface in two dimensions. Sentinel-1 (S1), which consists of Sentinel-1A (S1A) and Sentine-1B (S1B), can cover the entire Arctic Area within 2 days. The Interferometric Wide Swath (IW) mode, one of the main imaging modes of S1 in the Arctic, has been providing SAR images with a high resolution of 10 m. The ocean wave spectra are derived using the S1 IW data. The spectra are likely to provide vital observation for studying the interaction between sea-ice and waves as they present the distribution of wave energy and their variations in different frequencies and directions. Meanwhile, the retrieved ocean wave spectra are an excellent validation data source for SWIM, which is onboard the China-France Oceanography Satellite and provides measurements of ocean wave spectra in the global ocean.In this study, we use the sub-images from the S1 IW image with a size of 1024 pixels × 1024 pixels (10.24 km × 10.24 km) to retrieve ocean wave spectra by using a nonlinear retrieval method (i.e., MPI scheme). The retrieved ocean wave spectra and wave parameters are compared with SWIM measurements acquired at an incidence angle of 10°, and the Significant Wave Height (SWH) is measured at nadir. The SWIM spectrum covers a large area by approximately 70 km × 90 km. To make a comparison, all the SAR sub-image spectra within a SWIM beam coverage are averaged to calculate a new observed SAR spectrum, which is inputted into the MPI inversion scheme. Then, the retrieved SAR ocean wave spectrum is compared with the SWIM spectrum. The footprint size of the SWIM nadir beam is 18 km, which is comparable to the S1 sub-image size. Accordingly, a SAR sub-image is extracted at SWIM nadir, and the corresponding ocean wave spectrum is retrieved. The SWH is calculated by integrating the retrieved wave spectra to compare with the SWIM nadir measurements of SWH. The experiment was carried out using the data acquired in September 2020 in the Greenland Sea and Norwegian Sea, where the ocean waves generated in the North Atlantic and propagating vast distances to the ice-covered area in the Arctic ocean can be frequently observed. Fifty-four ocean wave spectra were retrieved from 25 IW data and are compared with the SWIM slope spectra. The comparison shows that the SAR-retrieved spectra are consistent with SWIM spectra in terms of structure and energy distribution. Good agreements are also found between the integral parameters of the SAR ocean wave spectra and SWIM slope spectra. The comparison yields a bias and an RMSE of 0.11 and 0.71 m for SWH and a bias and an RMSE of -0.52 s and 0.62 s for mean wave period. The comparison of the dominant wave parameters yield a bias of -7.74° and an RMSE of 15.75° for the dominant wave direction and a bias of -0.56 m and an RMSE of 52.73 m for the dominant wavelength. Furthermore, 5075 data pairs of S1-retrieved SWH and SWIM nadir SWH were collocated and compared. The comparison result yields a bias of 0.03 m, an RMSE of 0.48 m, and a correlation of 0.95. The comparison between the S1 retrieved results and the SWIM measurements suggests that ocean wave information can be effectively retrieved from S1 IW data by using the MPI method in the Arctic ocean. Although the MPI method relies on prior information, it is still an effective method for obtaining ocean wave spectra in high resolution. The spectra retrieved from S1 are likely to show the energy attenuation of ocean waves in different frequencies and directions when propagating toward an ice-covered area. This finding will be of great support for the further study on the interaction between sea ice and ocean waves. © 2023 National Remote Sensing Bulletin. All rights reserved.     

雷达卫星SAR与防卫气象卫星SSM/I对渤海海冰的观测研究

用雷达卫星(RADARSAT)合成孔径雷达(SAR)和防卫气象卫星(DMSP)特别微波成像辐射计(SSM/I)在1999年1月23日同一天对中国渤海区域海冰的观测数据，进行了主动SAR与被动SSM/I的组合研究.用一层海冰电磁散射辐射的建模与数值模拟，分析了中纬度渤海沿岸海冰的特征性变化，提出用SSM/I辐射亮度温度的散射指数、极化指数、极化比率来识别渤海海冰；用RADARSAT高分辨率SAR数据(水平极化后向散射系数)对渤海海冰物理特征进行识别与分类.这些特征指数在时间与空间尺度上的相关变化可有效地应用于渤海海冰的监察.     

Implications of broad dose distributions obtained with the single-aliquot regenerative-dose method on quartz fine-grains from loess

Six samples were collected from a section of Peoria Loess in Eustis, North America, for optically stimulated luminescence (OSL) dating of quartz, and all except one (LV90) produced narrow dose distributions. A comparative study was conducted on this sample and on a ‘well-behaved’ sample (LV91), involving other dating methods and examination of the quartz OSL. These investigations revealed differences in the quartz OSL growth with dose, OSL response to thermal treatments and the range of components within the OSL signals. An ultra-fast component was found in LV90 that displayed a higher rate of sensitivity change than the fast component and this had a malign influence on the determination of the equivalent dose. The distinctive luminescence characteristics of LV90 imply either a change in wind dynamics and/or the source area for the silt.     

InSAR复数影像配准方法探讨

复数影像配准是卫星合成孔径雷达干涉数据处理中的重要步骤之一,本文对目前影像配准最常用的方法进行了研究分析,提出了基于幅度影像的配准方案,并证明该方案配准精度高,适应能力强,运算效率高.     

UHV/VHF波段SAR对次地表层土壤湿度探测与反演

首先采用一层非球形粒子植被模型,计算各波段矢量辐射传输方程Mueller矩阵一阶散射解,对比微扰法所得各波段地表粗糙面直接后向散射解,结果证明L波段植被层的散射对观测结果仍有影响,与下垫土壤粗糙表面的散射不易分离.因此,宜采用更低频率的UHF和VHF波段,对地表和次地表层能有较大的渗透深度,并可忽略植被层影响.接着,运用矢量辐射传输的3层土壤全极化Mueller矩阵解,计算UHF/VHF波段分层土壤的散射与传输,分析该两波段探测深度的差异,证实UHF波段可探测大致10-60cm深处的土壤湿度,而VHF波段探测深度能更大一些.根据第3层中土壤体湿度变化0.1时能否引起土壤表面观测的后向散射系数变化0.1dB这一判据,分析VHF波段反演第3层土壤体湿度的必要条件,证实当第2层的体湿度较小时(<0.25)才能反演层3的体温度.基于UHF/VHF两波段探测深度的差异,耗散土壤层的贡献有不同的权重,先后采用UHF和VHF,迭代法实现3层土壤湿度廓线反演.误差分析表明,该方法是有意义的.     

地震损毁建筑物的高分辨率SAR图像模拟与分析

为了分析汶川地震震后高分辨率合成孔径(SAR)图像的城区建筑物特征,基于实际获取的机载X波段SAR图像,采用电磁模拟方法进行分析和研究.通过对城区的完整建筑和毁损建筑进行三维建模,采用射线跟踪的电磁计算方法和图像域积分的成像模拟方法得到不同受灾程度的建筑物SAR模拟图像.与真实SAR图像对比分析,提出的算法能够分析建筑物结构变化对SAR图像的影响,模拟主要的强散射点,能有效辅助SAR图像进行城区特征分析.