Storage-efficient reconstruction framework for planar contours

A storage-efficient reconstruction framework for cartographic planar contours is developed. With a smaller number of control points, we aim to calculate the area and perimeter as well as to reconstruct a smooth curve. The input data forms an oriented contour, each control point of which consists of three values: the Cartesian coordinates (x, y) and tangent angle θ. Two types of interpolation methods are developed, one of which is based on an arc spline while the other one is on a cubic Hermite spline. The arc spline-based method reconstructs a G¹ continuous curve, with which the exact area and perimeter can be calculated. The benefit of using the Hermite spline-based method is that it can achieve G² continuity on most control points and can obtain the exact area, whereas the resulting perimeter is approximate. In a numerical experiment for analytically defined curves, more accurate computation of the area and perimeter was achieved with a smaller number of control points. In another experiment using a digital elevation model data, the reconstructed contours were smoother than those by a conventional method.     

三次参数样条曲线在等高线矢量化中的应用

主要阐述了如何应用参数样条曲线在AutoCAD中快速准确地矢量化等高线。通过对三次参数样条曲线不同边界条件所绘线形的分析,最终选定了用夹持端样条曲线来矢量化等高线,并在AutoCAD平台上和VC 6.0开发环境下,用Object ARX开发工具包开发完成该矢量化程序。实践证明,该法切实可行,提高了工作效率。     

精密卫星钟差的内插方法研究

分析了几种精密卫星钟差的内插方法,将插值结果和IGS提供的数据相比较,采用统计量和差值图表来显示不同方法的插值结果,得到0.11～0.38 ns的精度.结果表明,内插方案的最佳选择取决于不同卫星钟类型的钟差变化特性以及插值方法本身的特点.     

扰动引力快速确定的替代算法

扰动引力的点质量表达方法使得扰动引力的计算模式具有最为简单的结构 ,但它本身也具有计算量大的缺点 ,与实时计算的目标相矛盾。为解决这一矛盾 ,文中用三次等距B样条函数进行试验。结果表明 ,在已知预设弹道坐标的条件下 ,三次等距B样条函数方法是适于进行快速精确逼近弹道扰动引力的数值算法     

等轴自然硅在我国的发现

我们发现的等轴自然硅产于安徽贵池的天然河流重砂中,共生矿物有锡石、锆石、磁铁矿、金红石、(?)石、四方硅铁矿、角闪石和辉石等,自然硅与四方硅铁矿呈连生体。自然硅为灰黑色,金属光泽,不透明,颗粒大小0.074—0.015mm,显微硬度1018.9kg/mm~2,密度2.32g/cm~3六个不同颗粒自然硅电子探针分析结果:Si 99.11,Fe 0.138,Al 0.034,Sn 0.16。最强的衍射线为3.12(10,111),1.91(8,220),1.635(8,311),1.105(7,422)。晶胞参数α=5.430A,V=160.103A~3,Z=8.     

基于三次插值的大地电磁自适应有限元二维正演模拟

大地电磁(MT)数值模拟中通常使用有限单元法,通过伽辽金(Galerkin)法将微分方程转化为与其等价的泛函形式,对泛函求取极值并在单元上定义插值基函数,得到节点上电磁场值的线性方程组,最终形成大型复对称稀疏矩阵。要达到较高的有限元计算精度,一般采用密集的网格或高次插值的方法,这样做大大的减慢了正演的速度。结合两者的优点利用三次插值和h-型自适应相结合的有限元法来实现MT的正演算法。首先从一个粗网格出发并利用三次插值,通过后验误差估计方法局部加密网格,在计算量较小的情况获得较高的计算精度。这种方法可以针对目标区域和介质分界面发生突变处进行网格加密,不需要全局加密网格。最后通过对国际标准模型COMMEMI-2D1的模拟,分别比较二次插值与三次插值的自适应网格数量和数值模拟结果,证明了三次插值自适应有限元算法的可行性。     

一种“准拉格朗日法”和“欧拉法”统一算法时间积分方案

从欧拉算符出发,用泰勒级数展开,给出二阶时空微商余项预报方程,进而讨论一种数值分析新算法——双三次曲面拟合（插值）的准拉格朗日时间积分方案与数值模式。它是将大气运动描述成为非线性的“三次”运动,即是通过对原始大气运动方程中包括标量、矢量的压、温、湿、风、以及“旋转地球上广义牛顿力”加速度场和散度场等,做双三次曲面拟合,实现对各个大气运动变量场的二阶可导,即限定气块上游点在各个不同双三次曲面（片）、具有斜率、曲率和挠率的非线性“三次”变量场上活动,从而可对各个大气运动方程做时间离散积分,即为“双三次曲面拟合—时间步积分—双三次曲面拟合……”,实现成为一种新动力框架数值模式。由于双三次曲面具有数学定律“收敛性”和二阶可导“最优性”,故选用双三次曲面插值求算二阶余差“上游点”,具有“充分必要”的数学理由：它包含了大气运动变量场之斜率、曲率和挠率。因此,埃尔米特双三次曲面片具有对“网格”变量场二阶可导运算“等价性”、及其数学“收敛性”与最佳曲率“最优性”,并且将“准拉格朗日法”与“欧拉法”,以及柯朗 弗里德里希斯 列维判据统一起来。容易实现全球“网格”变量场的双三次曲面拟合,和可按双三次曲面变量场的斜率、曲率或挠率判断,作变量场局域或单点平滑,以此保持“三次”模式的时间积分稳定性。     

地层数字制图中纵、横分界线控制搜索新算法

针对目前地层数字连线中处理地层缺失和倒转困难的问题,考虑6种地层缺失和倒转情形,提出了纵、横向结合地层分界线控制点搜索算法。将剖面图绘制视为一个多层地层区域组合叠加的过程,在保持地层连续性的前提下兼顾各地层相互之间的影响,有效处理了地层倒转、缺失问题,在此基础上提出一种增加虚拟点、微调已有点的新算法。结合3次样条插值函数,对绘制出的剖面图进行修正,实现了地层光滑连线,达到了地层无缝、无重叠组合的效果。将提出的复杂地层连线新方法应用于广东梅河高速公路工程数字制图中,构造出的地层分布图像准确性好,比较符合实际,可以为工程的勘察、设计和施工提供参考。     

山东焦家金成矿带三维预测模型及靶区优选

为了实现焦家金成矿带资源的可持续开发利用,应用三维地质建模技术指导深部找矿工作迫在眉睫。文章从山东焦家金成矿带区域地质背景入手,全面收集研究区已有的地质、矿产、物探、化探等科研成果资料,结合前人的成矿理论研究,分析和总结有利成矿、控矿地质条件,并应用该区深部物化探信息资料,建立综合信息定量找矿模型,采用“立方体预测模型”及地质统计学方法推断深部成矿空间,实现真正意义的定位定量定概率的资源预测评价。预测结果表明,该成矿带深部找矿金资源量达842.88 t,具有较好的开发利用潜力。更重要的是,圈定了7个尚未开发的成矿有利靶区,并结合探矿工程模型,提出了具体的勘查建议,为深部找矿提供了科学依据。