多年平均气温数据空间化误差的尺度效应

属性数据空间化是利用矢量数据生成栅格数据产品的有效方法,它有助于不同来源、不同格式之间的数据的综合分析。空间化是一种必然有误差伴随的过程,为探讨空间化误差与数据源密度、空间化模型方法,以及空间化分辨率之间的关系,本文利用7种水平的气象站点密度、5种空间化模型方法和19种栅格分辨率分析多年平均气温数据空间化误差与这3类影响因子之间的关系。分析发现：（1）气象站点密度的降低导致多年平均气温数据的空间化误差增加;（2）在IDW、Kriging、Adjusted IDW、Regression和Anusplin 5种空间化模型方法中,Adjusted IDW、Regression、Anusplin比IDW、Kriging的精度高;（3）随着栅格分辨率的变粗,多年平均气温数据空间化误差增大;（4）在影响空间化精度的3类因子中,空间化模型方法对空间化精度的影响最大,栅格分辨率次之,气象站点密度的影响最小。通过多元回归分析,建立了多年平均气温数据空间化误差与这3类影响因子之间的定量模型,可为空间化技术方案的制定提供参考和依据。     

面向稀疏散布数据集的时空Kriging优化

时空Kriging法通过将变异函数向时空域进行扩展得到时空变异函数,有效地利用时空邻近的采样点综合进行插值,由于时空稀疏散布数据集具有单一时刻下样本点数量少以及时空分布不规律的特点,难以满足使用时空Kriging插值法的基本条件,导致插值精度不高,据此本文提出了优化方法：通过多时段叠置拟合空间变异函数的方法,综合利用时空邻域内的采样点以解决单一时刻下空间邻域内数量不足情况;控制时间变异对空间变异函数拟合的误差影响;采用积合式模型构建时空变异函数进行插值。最后使用Argo海温数据进行插值实验,在相同条件下与时空Kriging法以及时空权重法的交叉验证结果对比得出,该方法在保证拟合所需采样点数量要求的同时,有效削减了一般时空Kriging法中时间变异对空间变异函数拟合结果的干扰,插值结果的绝对误差均值从0.5降低至0.2以内,稳定性进一步增强,改善了时空Kriging法在稀疏散布数据条件下精度上的不足。     

局部SVT算法的遥感反演场数据恢复实验分析

 遥感反演场数据会由于云雾、地物的遮挡,传感器性能等原因造成部分区域数据的缺失而影响遥感反演场数据的应用。矩阵填充理论针对低秩矩阵,利用矩阵的低秩性,即数据的高相关性,可以高精度地对低秩矩阵中的缺值数值进行恢复,其中矩阵填充理论中的SVT(Singular Value Thresholding)算法可以对矩阵中缺失数值进行快速、高精度的估计,应用广泛。本文应用矩阵填充理论的SVT算法,以缺值点为中心,方差最小作为窗口尺度选择的标准,这样可以保证区域数据的高相关性,建立局部窗口,对窗口进行SVT算法填充。本文也针对相同缺值区域进行了距离反比加权插值、Kriging插值法插值和整体SVT算法插值,整体SVT算法插值即并未对缺值点进行相关性窗口判断,而是直接对整个区域进行SVT填充。并对这几种方法的精度进行比较,得到局部SVT算法的精度相比整体SVT算法和距离反比加权插值算法的精度要高,与Kriging算法相比,其精度变化趋势相似,在锋面区域局部SVT算法精度比Kriging方法要高。     

几种土壤属性制图方法的稳定性与影响因素分析

 本文以土壤CEC和土壤全氮为研究对象,在河南省黄河以北6个地市选取870个样点,并随机均分为两个数据集,进行了制图对比分析;同时研究了Kriging法、IDW法和以点代面法制图结果的稳定性及其与精度的关系,以及影响因素。结果表明:采用实测数据与预测数据的交叉验证并不能用来衡量制图结果的稳定性,验证精度高低并不表示采用不同数据集制图结果可能会重现的概率;而且采用不同的方式进行精度检验,结论也会有所不同。采用不同抽样集合制图,Kriging法和IDW的结果较稳定,两者较接近,其相对误差>0.3的区域均不超过20%,且相对差异度高值区的空间分布格局相对分散;而以点代面法制图结果不稳定,相对误差>0.3的区域达到54.21%,且差异度较高的图斑相对集中,呈大片状分布。制图结果稳定性受到实测数据分布特征和局部地区土壤的高度变异性的影响,其中Kriging方法制图结果的不稳定性受样点分布格局的影响较另两种方法要大,而IDW法和以点代面法受实测数据自身变异性影响更明显。     

�߾����������ݸ����������Ƚ�

?о?3????????????????????Kriging????????????????????????????????????????????3????????????????μ??????????????????????3????????????????????????????????仯??????????????????Kriging????????????CoKriging?????Shepard????     

Kriging插值和序贯高斯条件模拟算法的对比分析

本文对Kriging插值与序贯高斯条件模拟值的算法联系进行了推导,并将两种计算结果和原始数据的统计参数作了对比。结果表明,Monto-carlo方法求得的序贯高斯条件模拟值经数学变换后等同于已知数据和此前模拟数据共同参与的Kriging插值结果与一个随机偏差的和,该随机偏差的均值为0,方差为Kriging误差方差。最优性条件导致Kriging插值结果的方差较原始数据降低了1个Kriging误差方差,造成Kriging平滑效应,其空间变异函数值降低,但自协方差函数值不变。序贯高斯条件模拟避免了平滑效应,其方差、变异函数和自协方差函数均不变,而其模拟值的误差方差较Kriging误差方差增加了1倍,表明1次随机模拟值的误差比Kriging插值大。然而,多次随机模拟值的平均值与Kriging插值的地理制图效果近似,可以弥补局部估值误差大的不足。因此,在应用中,Kriging插值是提供局部最优估计的方法,但却低估了全局的空间变异。而序贯高斯条件模拟的优点,在于提供若干等可能概率的模拟结果以进行估值的不确定性评价,并再现全局的空间可变性。     

普通Kriging法的参数设置及变异函数模型选择方法——以福建省一月均温空间内插为例

通过对普通Kriging法理论的深入理解,在ARCGIS和GS+软件平台上,以福建省一月均温为例,进行了深入的探索性空间数据分析,从数据的特征和规律出发探讨了普通Kriging插值法主要参数设置的依据和变异函数模型选择的方法。并用研究所获得的参数和变异函数模型对福建省一月均温的空间分布进行模拟,为合理的空间插值提供了一种较为清晰的思路。     

基于Kriging算法的三棱柱地质体三维可视化实现

采用了一种基于Kriging算法的三棱柱三维可视化方法对地质体进行了可视化,其基本思想是首先用Kriging算法将空间任意的离散的数据点生成规则的点数据,将点数据划分成多个三棱柱,然后以三棱柱为基本单元进行地质体可视化。本系统是在Visual C++/OpenGL平台上开发而成,具有精度高、交互好的特点,经过实践运行,切实可行。     

结合TRMM数据的区域降水高精度曲面建模研究

高精度曲面建模方法（High Accuracy Surface Modeling, HASM）,从理论上解决了传统方法在插值过程中峰值削平和边界震荡等问题。其模拟精度相对于经典插值方法有很大提高,已成功应用于人口密度、土壤属性,以及气候要素等领域的空间制图。然而,由于地面气象站点数量和分布的限制,使得HASM仅依靠站点数据难以得到高精度的空间降水估计数据,因此,本文以地貌与气候类型复杂多样的我国中西部地区2010年年降水量空间分布模拟为例,采用混合插值法进行HASM区域降水模拟。结果表明,TRMM作为背景场的HASM模拟的年降水量精度,在全局和局部明显优于IDW、Spline和Kriging等经典插值方法的结果,作为背景场的HASM模拟精度,MAE和RMSE分别为125.15 mm和155.80 mm,其他方法最好的模拟结果比其误差值分别高出53.6%和54.5%;其模拟误差在不同子区域都较小;各种方法在平原的精度都高于山区的精度。     

机载LiDAR点云密度和插值方法对DEM及地表粗糙度精度影响分析

机载LiDAR点云是获取高质量数字高程模型(Digital Elevation Model, DEM)的主要数据源,而地表粗糙度作为DEM的主要派生产品,在地学研究中发挥了重要作用,但点云密度和插值方法对DEM及地表粗糙度精度影响程度并没有明确结论。为此,本文利用不同地形条件下的林区机载LiDAR点云为实验对象,将原始点云随机缩减为不同的采样密度,利用5种常用插值方法(克里金(Ordinary Kriging, OK),径向基函数(Radial Basis Function, RBF),不规则三角网(Triangulated Irregular Network, TIN),自然邻域(Natural Neighbor, NN)和反距离加权(Inverse Distance Weighting, IDW))构建各个测区不同采样密度条件下的DEM,并通过空间特征和统计特征两方面对DEM及其地表粗糙度精度分析。结果表明:(1) DEM插值算法的精度随点云密度缩减而降低,且数据量缩减至原始数据量的30%后,不同算法精度区别较为明显,其中,RBF和OK精度最优,IDW精度最低;(2) DEM误差与...     

Application of wavelet packet in gravity anomaly processing

Based on the advantages of the wavelet to separate regional field and local anomalies in MATLAB environment,a high-precision regional-residual separation was finally realized. Analytical continuation and trend surface analysis are conventional methods for gravity anomaly separation. But the wavelet packet analysis in analyzing gravity data can make the gravity anomaly to be computed at a higher precision. In this paper,wavelet packet method is used to process gravity anomaly data obtained in Laos,and the separation result is good. Daubechies wavelet series has a higher precision in the wavelet packet.     

۰®���Ѵ����������������о�

????????????????????????????????о??????????4???4??????????????????????????????????????????????????????????????????????????????????????????????????????????????????????:??????????????????????Ч????á?     

实测重力异常与布格重力异常在高精度水准测量改正中的精度比较

按照一等水准测量要求在水准测点上进行重力测量。在水准测量数据处理中进行相同的水准标尺改正、正常水准面不平行改正、水准标尺尺长改正等，但在重力异常改正中采用实测重力异常和布格异常两种方法进行重力异常改正。比较两种方法对水准点高程改正的影响，并进行精度分析。     

���������ݶ����ݵ����񻯷����о�

??????????????????????????????????????????????????????????????淨????????????????Ч??????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????Ч????????????????????????????????????????????????????????     

�󻪱�����̬������������������о�

????μ????????о??????????????????????????????????????????????????????D1= 1.612????????????r=44 km??????????????????????????????????????????????????????????????????Ms5.0????????????????仯?????????????????????????????Ms5.0?????????н???????????     

小波多尺度分解特征分析

介绍了小波多尺度分析的数学原理，以三峡13期流动重力资料为例进行了小波分解。讨论了网格距对三峡流动重力资料分解结果的影响，提出了一种计算平均点距的方法；讨论了小波阶数与重力变化波长的关系．根据网格距及小波理论计算了各阶小波及逼近对应的波长范围。讨论了重力变化区移动与数据延拓、网格加密、网形完整性以及小波计算起始点之间的关系。结果表明．如重力网形不完整，分解结果中重力变化区随小波阶增加的移动没有规律．如果网形完整．分解结果普遍存在随小波阶增加而有规律移动的负重力变化区。重力变化区移动与小波分解起始点密切相关。     

Comparative analysis of fracture identification methods based on gravity anomaly in Heilongjiang

The identification of fractures is of great importance in gravity and magnetic data processing and interpretation.In this study,four fracture identification methods are applied,and widely used in processing and analysis of the gravity anomaly,including vertical second derivative method,tilt derivative method,theta map method and normalized differential method,for gravity data acquired in a given area in Heilongjiang.By comparing the distribution of the zero contour or maximum contour,we summarize the application effects,and both advantages and disadvantages of each method.It is found that tilt derivative method and normalized differential method provide better effects than other two methods:the narrower anomaly gradient belt and higher identification precision of fracture or geological boundary.The inferred fractures and geological boundaries have a great match with the results obtained from geologic map and remote sensing data interpretation.Those study results have definitely provided the theoretical foundation for identifying faults and the geological boundaries.     

�����������ƽ�������쳣��������г����

??????г????????????????????????????????????????????????????????????????????????????????????????????????????????????????????5????????????????г?????????????Stokes?????????￡????????????????????????г????????????????????????????????????????????     

宜川-泰安剖面重力异常与构造意义

利用宜川-泰安剖面高精度重力观测和GPS数据测量结果,对观测数据进行各项改正和处理,得到剖面的布格重力异常,并采用高斯低通滤波获得剩余重力异常。以诸城-宜川深地震测深（DSS）速度剖面为约束,对沿该剖面得到的高精度布格重力异常数据进行拟合,获得剖面地壳密度结构。结果表明,趋势异常主要反映该剖面地壳厚度的变化,太行山重力梯级带是地壳厚度的陡变带。结合前人相关研究,对剩余重力异常结合剖面经过的隆起、凹陷等地质问题进行讨论,进一步揭示出,剩余重力异常主要反映浅部地质构造和低密度地质体,剩余重力异常与地质构造不仅具有良好的对应关系, 而且其幅值能体现地质体的基本特征,比布格重力异常反映得更清晰、更细致。