基于地理探测器的大熊猫生境适宜度评价模型及验证

动物生境适宜度评价对于野生动物生境保护十分重要。基于物种活动点来建模的生态位模型是目前应用最广泛的动物生境评价方法,但该方法不能直接表达生境适宜度与环境因子间具有生态学意义的数量关系。本文以雅安地区为例,提出一种新的大熊猫（Ailuropoda melanoleuca）生境适宜度评价方法,选取海拔、坡度、坡向、地形指数、距水源距离、植被类型、主食竹及距公路距离8个环境因子,引入地理探测器,在分别基于MAXENT模型和层次分析法（The Analytic Hierarchy Process,AHP）所构建生境适宜度模型的基础上,通过4个地理探测器（风险探测器、因子探测器、生态探测器和交互作用探测器）探寻大熊猫生境与各环境因子间的关系以及环境因子对大熊猫生境的影响机理,并将其预测结果与单一MAXENT模型和AHP法进行对比。结果表明：（1）AHP、AHP-Geogdetector、MAXENT和MAXENT-Geogdetector模型总体评价精度分别为85.6%、86.5%、91.3%和94.2%,kappa系数分别为0.699、0.718、0.821和0.882,AUC值分别为0.902、0.928、0.949和0.966,模型所预测的适宜和较适宜区与实际分布区重叠比分别为63.66%、61.30%、76.70%和90.10%,说明AHP-Geogdetector和MAXENT-Geogdetector模型精度均比相应的单一模型有所提高,且MAXENT-Geogdetector模型精度最高;（2）基于地理探测器的大熊猫生境适宜度评价模型能以“生境适宜度和环境因子间具有生态学意义的数量关系”的形式直接体现环境因子对动物生境利用的生态学作用,具有较好的生态学可解释性。因此,用地理探测器进行大熊猫生境适宜度评价具有较好的可行性。     

西藏羊湖二厂电站厂房区泥石流危险性评价

泥石流的活动特征及其危险性评价对电站的施工及运行有重要意义.根据羊湖电站厂房区勘查所获取的资料,对厂房区后山5条泥石流沟的成因进行了分析,在此基础上将层次分析法(AHP)应用到泥石流危险性评价研究中.对该厂房区5条泥石流沟的危险度评价表明,5号沟泥石流危险度最高,4号沟次之,1号沟、2号和3号沟危险度相对较低.泥石流危险度与流域面积的大小呈正比关系.     

基于内容的离散余弦变换域自适应遥感图像数字水印算法

提出一种适于遥感图像知识产权保护和机密信息传递的离散余弦变换域自适应数字水印嵌入算法,该算法能够通过数字水印置乱、嵌入位置自适应选取、嵌入深度智能调节等措施,将二值水印图像信息安全地嵌入到遥感图像(纹理区)的离散余弦变换域中高频系数内.实验结果表明:该遥感图像数字水印算法不仅具有较好的透明性,而且对诸如叠加噪声、JPEG压缩、平滑滤波、几何剪切、图像增强、马赛克等攻击均具有较好的抗差性,整体性能优于Ho等提出的数字水印嵌入方案.     

Discrete displacement analysis for geographic linear features and the application to glacier termini

The extent of glacier terminus displacement is instrumental in investigations of natural or artificial geographic changes. Its importance to earth science and engineering is reflected in the considerable efforts that have been devoted to the development of several boundary displacement analysis methods. Among the methods, the buffering-based approach compares favorably with other approaches in objectivity and robustness. However, it does not consider the relative positions of boundaries, because its buffering operation cannot determine features' relative directions. This limitation incurs inaccurate calculation results – underestimation of mean shifts and overestimation of shape variations, especially when the two compared boundaries intersect. Discrete displacement analysis (DDA), an alternative method that considers given geographic objects as a set of a finite number of points, is proposed here. In a series of tests carried out, including Jakobshavn glacier's calving front, DDA was found to correctly calculate mean shift and shape variation even in cases where the conventional buffering-based method failed. Moreover, this approach is independent of the dimension of space in which it is implemented, and thus is expected to be utilized for analysis of 3D geographic object displacement.     

The biaxial compression mechanical properties of crushed rock

Crushed rock subgrade, as one of the roadbed-cooling methods, has been widely used in the Qinghai-Tibet Railway. Much attention has been paid on the cooling effect of crushed rock; however, the mechanical properties of crushed rock are somehow neglected. Based on the discrete element method, biaxial compression test condition for crushed rock is compiled in FISH language in PFC^2D, and the natural shape of crushed rock is simulated with super particle "cluster". The effect of particle size, crushed rock strength and confining pressure level on overall mechanical properties of the crushed rock aggregate are respectively analyzed. Results show that crushed rock of large particle size plays an essential framework role, which is mainly responsible for the deformation of crushed rock aggregate. The strength of gravel has a great influence on overall mechanical properties which means that strength attenuation caused by the freeze thaw cycles cannot be ignored. The stress-strain curves can be divided into two stages including shear contraction and shear expansion at different confining pressures.     

On-line estimation of unit hydrographs using the wavelet-based LMS algorithm / Estimation en ligne des hydrogrammes unitaires grâce à l'algorithme des moindres carrés moyens à base d'ondelettes

Abstract

This study applies the discrete wavelet transform (DWT) to decompose the unit hydrograph, thereby generating parsimonious reparameterizations of the unit hydrograph. A model compression method is then employed to significantly compress the unit hydrograph requiring that fewer coefficients be estimated. Moreover, a wavelet-based linearly constrained least mean squares (WLCLMS) algorithm is also used to estimate on-line the wavelet coefficients of the unit hydrograph. The updated wavelet coefficients of the unit hydrograph, convoluted with effective rainfall input in the wavelet domain, allow for accurate prediction of one-step-ahead runoff in the time domain. The proposed approach allows the unit hydrographs to vary in time and accurately predicts runoff from a basin in Taiwan, thus making it highly promising for flood forecasting.     

基于地震信号波形形态差异的面波噪声稀疏优化分离方法

实际地震信号通常可表示为具有波形特征差异的多种基本波形信号的线性组合,如叠前道集中的工频干扰噪声与有效波信号、面波噪声与体波信号等.选择单一数学变换方法,往往不易实现地震信号的稀疏表示.近年来发展的形态成分分析理论,通过联合多种数学变换,可实现对复杂信号的稀疏表示.本文根据单道地震记录中面波与体波信号波形结构特征的差异性,提出一种基于形态成分分析的面波噪声衰减方法.针对面波的低频、窄带以及频散特性选择一维平稳小波变换作为其稀疏表示字典,而针对体波波形的局部相关特性选择局部离散余弦变换作为其稀疏表示字典,建立基于双波形字典的形态成分分析模型,通过求解该稀疏优化问题获得最终的信噪分离结果.理论模型和实际地震资料处理证实该方法不仅能够衰减单炮地震记录中的强面波干扰噪声,同时能够更好地保护有效信号的波形特征与频谱带宽,为地震资料的后续处理和分析提供良好的数据基础.     

全球六边形离散格网数据的三维可视化方法

针对全球六边形离散格网的三维显示化方法开展研究,设计了一种六边形格网的空间层次结构(hexagonalquaternary balanced structure,HQBS),采用四位码元对格网单元进行编码,定义并实现了格网向量的基本运算,利用这些运算可以方便地实现格网单元的空间索引。在此基础上还研究了全球离散格网的动态生成与显示算法、可视化区域裁剪等相关内容。试验表明:全球格网动态生成的效率110～370单元/ms之间,加载空间数据后,格网数据和空间数据逐层加载的时间在300 ms左右,能够保证加载空间数据后的显示刷新率在20帧/s左右。     

利用数值投影变换构建全球六边形离散格网

本文针对全球离散格网构建过程中从平面格网到球面格网的关键步骤进行讨论,提出了一种新型的评价球面离散格网几何属性最优化的目标函数,使用遗传算法优化,得到了球面上有限层次内最优化条件下的直接剖分格网。利用有限层次内最优化格网提供的控制点数据,结合数值投影变换理论,成功地构建了几何属性更加均匀的全球六边形离散格网系统。实验表明,相对于现有Snyder等积投影建立的全球格网,在格网单元的均匀度上更优;在运算效率方面,速度大约是Snyder投影的2.5-3倍。     

Comparative study on two methods for calculating the gravitational potential of a prism

The determination of the gravitational potential of a prism plays an important role in physical geodesy and geophysics. However, there are few literatures that provide accurate approaches for determining the gravitational potential of a prism. Discrete element method can be used to determine the gravitational potential of a prism, and can approximate the true gravitational potential values with sufficient accuracy (the smaller each element is, the more accurate the result is). Although Nagy’s approach provided a closed expression, one does not know whether it is valid, due to the fact that this approach has not been confirmed in literatures. In this paper, a study on the comparison of Nagy’s approach with discrete element method is presented. The results show that Nagy’s formulas for determining the gravitational potential of a prism are valid in the domain both inside and outside the prism.