基于MATLAB的基本矩阵随机采样鲁棒估计

基本矩阵的估算可以分为线性算法、非线性优化迭代算法、鲁棒估计算法三大类,而以鲁棒估计效率、效果最为突出,但是实现比较复杂。本文通过数码相机获取的两张未标定照片,采样互相关系数法和RANSAC算法,在MATLAB中实现了基本矩阵的自动拟合,实验表明拟合精度较高、速度快,为下一步三维场景的自动重建和量测奠定了基础。     

Decumbent development: Urban sprawl in the Guwahati Metropolitan Area,India

Urban sprawl has become a global phenomenon as an outcome of growing population and rapid urbanization. Previous studies have addressed the rising incidence of uncontrollable urban development, particularly in peri-urban areas of cities, leading to chronic urban sprawl. The city of Guwahati, a million city in north east India, has expanded significantly in recent years. In this article, the links between population and growth of built-up areas were examined using geo-spatial techniques and remotely sensed datasets. The results indicate that the sprawl has accentuated in recent years. The intensity of land use remained uneven due to marked variations in the distribution of built-up areas, plausibly an outcome of unplanned urban growth. If current trends are anything to go by, future urban sprawl could pose serious threats to the vulnerable eco-sensitive and peri-urban areas of Guwahati. Secondary cities have unfortunately received scant attention in urban policy research, and Guwahati, epitomizes urban woes in a developing country.     

Wall-to-wall spatial prediction of growing stock volume based on Italian National Forest Inventory plots and remotely sensed data

Spatial predictions of forest variables are required for supporting modern national and sub-national forest planning strategies, especially in the framework of a climate change scenario. Nowadays methods for constructing wall-to-wall maps and calculating small-area estimates of forest parameters are becoming essential components of most advanced National Forest Inventory (NFI) programs. Such methods are based on the assumption of a relationship between the forest variables and predictor variables that are available for the entire forest area. Many commonly used predictors are based on data obtained from active or passive remote sensing technologies. Italy has almost 40% of its land area covered by forests. Because of the great diversity of Italian forests with respect to composition, structure and management and underlying climatic, morphological and soil conditions, a relevant question is whether methods successfully used in less complex temperate and boreal forests may be applied successfully at country level in Italy.For a study area of more than 48,657 km² in central Italy of which 43% is covered by forest, the study presents the results of a test regarding wall-to-wall, spatially explicit estimation of forest growing stock volume (GSV) based on field measurement of 1350 plots during the last Italian NFI. For the same area, we used potential predictor variables that are available across the whole of Italy: cloud-free mosaics of multispectral optical satellite imagery (Landsat 5 TM), microwave sensor data (JAXA PALSAR), a canopy height model (CHM) from satellite LiDAR, and auxiliary variables from climate, temperature and precipitation maps, soil maps, and a digital terrain model.Two non-parametric (random forests and k-NN) and two parametric (multiple linear regression and geographically weighted regression) prediction methods were tested to produce wall-to-wall map of growing stock volume at 23-m resolution. Pixel level predictions were used to produce small-area, province-level model-assisted estimates. The performances of all the methods were compared in terms of percent root mean-square error using a leave-one-out procedure and an independent dataset was used for validation. Results were comparable to those available for other ecological regions using similar predictors, but random forests produced the most accurate results with a pixel level R² = 0.69 and RMSE_% = 37.2% against the independent validation dataset. Model-assisted estimates were more precise than the original design-based estimates provided by the NFI.     

雷达-雨量计联合估测区域降水量方法研究Ⅰ——模式集成法的提出

天气雷达估测区域降水量是其应用研究的重要领域,不断研究改进雷达估测区域降水量精度和稳定度是雷达气象工作的重要目标。文章回顾了集合预报理论基础及其发展概况,总结分析集合预报方法、产品表现形式等。对前期雷达估测区域降水量的8种方法,从基本原理和特点等方面进行了系统分析。对该理论在雷达估测区域降水领域的适用性进行了分析,阐述了集合预报理论在该领域应用的概念模型,并提出了两种集成估测区域降水量的新方法,为下一步的应用研究提供了理论支撑。     

基于网络时空核密度的交通事故多发点鉴别方法

交通事故多发点是道路交通安全管理的重要治理对象,如何利用空间统计方法对其进行高效鉴别是研究热点。以华东某地为研究区域,以2013~2015年该研究区域的道路交通事故数据为研究对象,以时空道路网络为视角,通过路网匹配和路网裁剪形成事故时空子路段,提出一种基于交通事故场景的网络时空核密度估计值作为鉴别指标,用累计频率法和零膨胀的负二项回归模型确定鉴别阈值的事故多发点鉴别方法。     

建筑物侧面轮廓线约束的高分辨率遥感影像建筑物高度估算方法

针对现有利用阴影长度法提取建筑物高度时存在的阴影间相互遮挡问题,提出了一种基于建筑物侧面轮廓线进行建筑物高度估算的新方法。首先,利用RPC模型计算建筑物像点位移的方向与卫星成像角度,再将遥感影像进行旋转,使建筑物像点位移沿水平方向;然后,利用Canny算法进行轮廓检测,并构建一定长度的矩形形态学结构元素,对轮廓图像进行形态学开运算,以提取侧面轮廓线,再利用Hough变换与建筑物角点约束,对所提取的轮廓线进一步筛选;最后,根据卫星侧视成像时建筑物高度与像点位移的几何关系进行建筑物的高度估算。利用实际的高分辨率卫星影像对本文方法进行了验证,并与阴影法估算建筑物高度进行了对比。试验结果证明,利用建筑物侧面轮廓线进行建筑物高度估算平均误差可以达到0.7 m,且实际精度优于使用阴影法进行建筑物高度估算。     

基于椭球不确定性的平差模型与算法

测量平差模型中的参数通常存在一些不确定的附加信息或先验信息,充分利用它们可以对部分参数进行约束,从而保证参数解的唯一性和稳定性。本文利用椭球集合描述不确定性,建立了一个新的带有椭球不确定性的平差模型。以两个椭球交集的外接椭球的特征矩阵的迹最小平差准则,分析了不确定度的传播规律,给出了带有椭球不确定性的平差方法。最后,通过算例验证了算法的有效性,说明了平差解与带权混合估计的关系。     

Estimation of Grassland Production in Central and Eastern Mongolia from 2006 to 2015 via Remote Sensing

Mongolia is an important part of the Belt and Road Initiative “China-Mongolia-Russia Economic Corridor” and a region that has been severely affected by global climate change. Changes in grassland production have had a profound impact on the sustainable development of the region. Our study explored an optimal model for estimating grassland production in Mongolia and discovered its temporal and spatial distributions. Three estimation models were established using a statistical analysis method based on EVI, MSAVI, NDVI, and PsnNet from Moderate Resolution Imaging Spectroradiometer (MODIS) remote sensing data and measured data. A model evaluation and accuracy comparison showed that an exponential model based on MSAVI was the best simulation (model accuracy 78%). This was selected to estimate the grassland production in central and eastern Mongolia from 2006 to 2015. The results show that the grassland production in the study area had a significantly fluctuating trend for the decade study; a slight overall increasing trend was observed. For the first five years, the grassland production decreased slowly, whereas in the latter five years, significant fluctuations were observed. The grassland production (per unit yield) gradually increased from the southwest to northeast. In most provinces of the study area, the production was above 1000 kg ha ^-1, with the largest production in Hentiy, at 3944.35 kg ha ^-1. The grassland production (total yield) varied greatly among the provinces, with Kent showing the highest production, 2341.76×10 ⁴ t. Results also indicate that the trend in grassland production along the China-Mongolia railway was generally consistent with that of the six provinces studied.     

Ekman模型时间变化风应力系数的伴随参数估计

本文基于时间分布参数设置,利用伴随同化方法,反演了Ekman模型中随时间变化的风应力拖曳系数,并在孪生实验和实际实验中对该方法进行了验证。在孪生实验中,研究了参数反演结果对不同影响因素的响应,包括:风速分布、风应力系数分布、风应力系数初始猜测值、风应力系数独立变量个数、观测数据误差和观测的深度。孪生实验结果验证了伴随同化方法反演Ekman模型中时变风应力系数的有效性,具体包括如下五个方面结论:1)不同风速分布下均能成功反演出不同风应力拖曳系数分布; 2)反演结果对初始猜测值较为敏感,风应力系数初始猜测值越接近给定值,反演结果越好;3)风应力系数独立点个数的选取会显著影响反演结果,合理的选择有利于提高反演效率及减小观测数据误差;4)观测误差能够影响反演结果,观测数据误差在20%以下时能取得合理的反演结果; 5)反演结果对观测数据的表层和次表层流速更为敏感,这是由Ekman流的物理性质决定的。实际实验,利用百慕大锚系试验平台的风速和流速数据,去除周期性潮流和地转流成分后得到Ekman流成分,并作为观测输入到该同化模型,反演出了适用于该区域和该时段的随时间变化的风应力系数。通过比较模拟流速和观测流速,证明利用伴随同化方法能从实测数据中反演出合理的时变风应力系数,对于海洋模型风应力系数的确定是一项有益的尝试。     

Sea state estimation using multiple ships simultaneously as sailing wave buoys

The article presents initial ideas towards a network-based approach for sea state estimation used for marine operations and other maritime applications. In principle, all available means, ranging from in situ buoys, fleet of ships to remote sensing by satellite and aircraft, could be considered, emphasising that each means and any combinations among may act simultaneously. This study focuses on just one of the means; the use of ships as sailing wave buoys. The article introduces the wave buoy analogy, i.e. ship-as-a-wave-buoy, and it makes a proposal on how to impose (different) weights to the single ship-specific wave spectrum estimates obtained from multiple ships. Moreover, the work includes a discussion about the importance to associate a measure to reflect the (un)certainty of the wave spectrum estimate. The article presents a numerical case study, where multiple ships act simultaneously as wave spectrum-estimators. The case study relies on numerical motion simulations, as appropriate full-scale data is not yet available. In the analysis, it is shown that the use of simultaneous data from multiple ships leads to more accurate wave spectrum estimations.