基于距离图像的分块曲面LiDAR点云滤波方法

机载激光雷达点云滤波是获取高精度数字高程模型的关键,也是国内外研究的热点.这里采用区块索引机制对激光雷达点云进行分块曲面滤波;利用距离图像的像素位置信息指导分块数的确定.试验结果表明,该方法可有效滤除非地面点,从而保留地面部分点云,具有一定实用价值.     

一次范数最小估计的两种算法

测量数据服从正态分布情况下,最小二乘估计具有最优统计性质,如果测量数据受到粗差干扰而偏离正态分布,就会使最小二乘估计的最优性受到严重冲击,而且估计失实。然而,若采用一次范数最小估计,则所得结果比较好。介绍一次范数最小估计(L1估计)的两种算法——选权迭代法和线性规划法,并通过Matlab模拟实验比较这两种方法的优缺点,证明L1估计的良好抗差性。     

高精度曲面模型解算改进的Gauss-Seidel法

为了降低HASM的时间复杂度,采用一种改进Gauss-Seidel(GS)算法(MGS)解算HASM方程组。首先,从理论上分析了MGS算法收敛速度快于GS算法,然后以高斯合成曲面作为研究对象,用四组模拟试验表明,相同的网格数、达到相同的计算精度,MGS算法计算时间小于GS算法,且两种算法时间差与模拟区域网格数呈二次线性相关;固定网格数,使用相同的内迭代或者外迭代次数,MGS算法精度高于GS算法,但增加内迭代或者外迭代次数,GS算法同样收敛;MGS算法计算时间与网格数呈线性相关。MGS算法能够有效解决HASM模拟大区域的计算时间瓶颈,提高HASM运算速度。以甘肃省董志塬某测区SRTM3作为研究对象,基于MGS的HASM用于模拟DEM表明,HASM精度要高于传统的插值方法。     

大规模影像数据库的设计与实现

本文对现有遥感影像数据管理方法进行了讨论,针对其存在的不足,沿用传统数据库技术成果,如SQL技术、索引技术,并利用SQL Server 2005简单易用,且支持二进制大对象(BLOB),提出分幅分块存储管理方式,解决数据的物理存储问题,结合金字塔技术、索引技术、JPEG压缩技术,解决了数据的快速加载问题。工程应用实际表明,该方法能够满足海量遥感影像数据的管理需求。     

基于欧氏距离测度的激光点云配准

针对激光扫描测量系统得到的多视角点云数据,提出用离散对应特征和迭代最近点(ICP)算法相结合的方法,对点云模型多视配准技术进行了研究。首先给出单位四元数和旋转矩阵的关系以及线性最小二乘法原理;然后利用基于离散对应特征的方法求出刚体变换的一个初值,并用迭代最近点(ICP)算法精确估计刚体变换参数;最后用工件36副真实点云模型的配准结果证明此方法的有效性。     

计算Legendre函数导数的非奇异方法

引力场关于经度和纬度方向的梯度在两极附近会产生奇异性现象,这将会给诸如重力场和静态洋流探索(GOCE,Gravity field and stesdy-state Oceam Circulation Explorer)数据处理等引力场的研究工作带来诸多不便和困难。这里首先分析了该奇异性产生的原因,即目前采用的球坐标系自身在两极处是奇异的;然后利用Legendre函数的性质推导了一组不含任何奇异性的计算引力场梯度的计算公式;最后与常用的迭代方法进行了实例计算比较,结果表明所导出的公式不仅计算精度大大提高,而且计算用时也不会增加。     

感兴趣区域CT图像重建方法及模拟实验

虽然CT技术已经发展得相当成熟,但在大物体成像、高分辨率成像和减少辐射剂量等方面现有的CT成像方法仍然存在较大的困难。事实上,很多工程应用中并不要求对完整的物体进行全局CT成像,只需要获得某些感兴趣区域的物体图像即可,特别是医疗临床诊断中只要能够实现对可疑病灶部位的成像即可。因此,本文研究了针对感兴趣区域的CT图像重建方法,以及X射线束的视野只覆盖感兴趣区域的扫描方法设计,感兴趣区域CT成像研究的关键是如何在投影数据截断情况下实现断层图像的精确、稳定重建;本文介绍了三种投影数据在不同截断方式下的图像重建方法,并设计了相应的数值模拟实验,给出了数值模拟实验结果。最后,对感兴趣区域CT成像技术在工程中应用潜力做了展望。     

基于敏感性迭代拟合算法反演汶川Mw 8.0地震InSAR同震形变场

基于单一断层模型,运用敏感性迭代拟合算法反演了汶川地震InSAR同震形变场,获得了断层滑动分布及部分震源参数.结果表明,倾角线性变换的单断层模型模拟的同震形变场与InSAR形变场吻合较好,且残差较小,平均残差为0.11 m;反演的滑动分布主要集中于地下0～20km,其中,汶川地区在震源深度附近有较大的滑动量,10～15 km深度最高可达7 m,地表滑动量却相对较小,平均值仅为2.5 m,平均滑动角均值约为121°;北川地区最大滑动量可达到10 m,平均滑动角均值约为109°;而青川地区10～15 km深度滑动量最高也达8 m,平均滑动角均值约为135°;滑动矢量图显示,沿SW-NE走向断层面的滑动方向以强烈的逆冲为主,兼有一定右旋走滑分量.反演矩张量为1.0×1021 N·m,矩震级达Mw8.0.     

初至波地震层析技术及其在四川复杂地区的应用

近地表速度场的重建是获得准确静校正量的关键。根据实际地震数据的初至波时间建立初始速度模型;采用改进的最短路径射线追踪算法进行正演;利用改进的约束带阻尼联合迭代重建技术(CDSIRT)求解大型、稀疏的方程组来反演近地表的速度场。理论模型和实际资料的反演结果表明,该方法能稳定、快速、准确地重建复杂地区的近地表速度场,建立准确的近地表速度场,从而可计算准确的静校正量,保证复杂地区的构造成像。     

声波测井中的广义反射、透射系数方法

本文将广义反射、透射系数矩阵方法推广用于研究声法测井中的声波传播问题.结果表明,测井中的声场解可以表示成无数条广义声线的叠加.任意一条广义声线的积分表达式,则可根据每一条声线段的性质,由各反射界面上的广义反射、透射系数直接给出.这种表示方法特别适用于研究和识别声波记录图中各种“震相”的性质.当计算全波声场时,由于本文将声场计算中的被积函数形式上分解成一系列（2×2）子矩阵的乘积,并给出了相应的迭代公式,从而可以用较快的速度进行计算.