航磁数据三维插值的实现和插值方法优选
论文--三维电阻率数据插值加密及数字图像插值的研究报告
.摘要本文主要分析比较了不同的插值方法对于已有的电阻率数据运用插值加密后的效果并进行了一系列的评价。
针对问题一，对于三维空间的插值加密可以运用的有线性插值法、三次样条插值法、三次多项式插值法、最邻近插值法等，对各个插值方法进行的深入分析，理解各个一维插值公式的推导过程以及一维至三维的衍生原理，得出保证极值大小及空间位置不变的两种方法，继而通过直接证明以及间接证明，通过计算机对于给定数据的运算，验证两种方法的可行性。
最后采用得出的两种方法，最邻近插值法以及线性插值法，进行三维分析并建立数学模型，得出最后给定点的电阻率数据。
针对问题二，通过第一题的给出的线性插值法以及最邻近插值法，分别计算出网格大小1m*1m*1m上任意点的电阻率数据，对所得数据进行三维建模，即可得出在网格加密后的电阻率数据三维成像结果。
对于插值方法的由一维至三维的延伸，比较在一维延伸至三维时各个方法的计算量，评估出网格加密过程中计算复杂度及计算量大小。
同时运用1m*1m*1m以及10m*10m*10m网格大小上所有点的电阻率得出各种情况下的平均值及标准差。
针对问题三，用颜色图展示网格加密后的直观效果。
将最小值置为纯蓝色（RGB 为(0,0,255)），中间值置为纯绿色（RGB为(0,255,0)），最大值置为纯红色（RGB 为(255,0,0)），对于中间数值，运用RGB成色三维矩阵，找出电阻率数值大小与颜色坐标的对应函数关系，用计算机作出与要求完全相同的Z=40对应的原图像。
继而绘制出对切片Z=0,50分别给出原数据，两种方法加密网格后数据的颜色对比图。
以及对切片X=82，Y=47, Z=88 ，两种加密方法得到数据的颜色图。
针对问题四，对于不同插值方法的加密效果，可以通过加密后切片图像的清晰度来表示其加密效果的好坏。
将RGB三维矩阵转换为灰度数值，而对灰度进行若干种定量的表示即可反映数字图像清晰度。
得某一指定切片图像的灰度值，通过4种不同方法对灰度处理得出相应定量指标以评价图像清晰度，进而反映两种插值加密法的效果。
航空摄影测量三维数据获取方法及研究
航空摄影测量三维数据获取方法及研究以航空摄影测量三维数据获取方法及研究为题，本文将详细介绍航空摄影测量在三维数据获取方面的应用和研究进展。
一、航空摄影测量的基本原理航空摄影测量是一种以航空摄影为手段，对地面目标进行三维空间位置测量的技术。
其基本原理为：通过航空相机对地面目标进行拍摄，再通过对航空影像的解译，测量出目标在三维空间中的位置和形态。
二、航空摄影测量三维数据获取方法1.立体像对法立体像对法是航空摄影测量中最常用的方法之一。
该方法是通过两张或多张相互重叠的航空影像进行立体测量，从而得出地面目标的三维坐标。
该方法的优点在于能够快速、准确地获取三维数据，并且适用范围广泛。
2.数字正射影像法数字正射影像法是一种利用数字影像进行三维坐标测量的方法。
该方法通过对数字正射影像进行解译，得出地面目标的高程信息，从而实现三维坐标的获取。
该方法的优点在于能够快速、准确地获取高程信息，并且不受地形起伏的影响。
3.激光雷达法激光雷达法是一种利用激光技术进行三维坐标测量的方法。
该方法通过激光雷达对地面目标进行扫描，从而得出地面目标的三维坐标。
该方法的优点在于能够快速、准确地获取三维数据，并且不受地形起伏的影响。
三、航空摄影测量三维数据获取的应用1.地形测绘航空摄影测量在地形测绘方面有着广泛的应用。
通过航空摄影测量技术，可以对地面的地形进行高精度的测绘，从而为城市规划、地质勘探等领域提供基础数据支撑。
2.资源调查航空摄影测量在资源调查方面也有着广泛的应用。
通过航空摄影测量技术，可以对林地、水资源等进行高精度的测量，从而为资源的管理和保护提供基础数据支持。
3.城市管理航空摄影测量在城市管理方面也有着广泛的应用。
通过航空摄影测量技术，可以对城市的道路、建筑等进行高精度的测绘，从而为城市规划和管理提供基础数据支持。
四、航空摄影测量三维数据获取的研究进展随着科学技术的不断发展，航空摄影测量三维数据获取的研究也在不断深入。
目前，研究者们主要关注以下几个方面：1.精度提升精度是航空摄影测量三维数据获取的重要指标之一。
三维插值算法原理
三维插值算法原理三维插值算法是一种用于在三维空间中估计未知数据点的方法。
它基于已知数据点的数值和位置，通过插值计算出未知数据点的数值。
三维插值算法广泛应用于地理信息系统、计算机图形学、物理模拟等领域。
三维插值算法的原理是利用已知数据点之间的关系，通过数学方法推断出未知数据点的数值。
常用的三维插值算法包括线性插值、三次样条插值和Kriging插值等。
线性插值是最简单的三维插值算法之一。
它假设未知数据点的数值与已知数据点之间的关系是线性的，并通过线性方程进行插值计算。
线性插值适用于数据点之间变化较为均匀的情况，但对于数据点之间的非线性关系无法准确估计。
三次样条插值是一种更为精确的三维插值算法。
它利用已知数据点之间的关系构建出一组三次样条函数，并通过这些函数进行插值计算。
三次样条插值可以更好地逼近数据点之间的曲线变化，能够较好地处理非线性关系。
Kriging插值是一种基于统计学原理的三维插值算法。
它假设未知数据点的数值是具有空间相关性的随机变量，并通过已知数据点之间的空间相关性进行插值计算。
Kriging插值能够考虑数据点之间的空间分布和变异性，适用于非均匀分布的数据点。
在三维插值算法中，选择合适的插值方法是十分重要的。
根据实际数据的特点和需求，可以选择不同的插值方法。
线性插值适用于数据点变化较为均匀的情况，速度较快但精度较低；三次样条插值能够更好地逼近数据点之间的曲线变化，精度较高但计算复杂度较大；Kriging插值考虑了数据点之间的空间相关性，适用于非均匀分布的数据点，但计算量较大。
除了选择合适的插值方法，还需要注意插值参数的选择。
参数的选择会直接影响插值结果的准确性和稳定性。
例如，在Kriging插值中，需要选择合适的半方差函数和变异性参数来描述数据点之间的空间相关性。
参数的选择通常需要根据实际数据的特点进行调整，以获得最佳的插值效果。
三维插值算法是一种通过已知数据点之间的关系推断未知数据点数值的方法。
无人机航磁测量工作内容
无人机航磁测量工作内容一、引言无人机航磁测量是利用无人机搭载磁力仪器对地面磁场进行测量的一种技术方法。
相比传统的地面磁测方法，无人机航磁测量具有高效、快速、灵活等优势，因此在地质勘探、资源调查、环境监测等领域得到广泛应用。
本文将介绍无人机航磁测量的工作内容及相应的技术要点。
二、飞行计划在进行无人机航磁测量前，需要制定详细的飞行计划。
首先，根据测区的地理特点和测量目的确定飞行区域。
然后，根据飞行区域的大小和形状，确定航线的起点、终点和间距，保证测量的全覆盖。
同时，考虑到无人机的续航能力，合理安排航线的长度和飞行高度。
最后，根据飞行计划确定起飞点和降落点，并考虑到无人机的安全飞行。
三、无人机设备准备在进行航磁测量前，需要对无人机设备进行准备。
首先，确保无人机具备搭载磁力仪器的能力，包括传感器的安装和校准。
其次，检查无人机的飞行控制系统和通信设备是否正常工作，以确保飞行的安全和数据的传输。
最后，对无人机进行全面的功能测试和预飞检查，确保其各项性能符合要求。
四、飞行操作在进行无人机航磁测量时，需要进行飞行操作。
首先，进行起飞前的准备工作，包括确认飞行控制系统和传感器的工作状态，检查飞行区域的环境条件，确保飞行的安全。
然后，根据飞行计划进行起飞、航线飞行和降落等操作。
在飞行过程中，要注意保持与地面站的通信，及时接收和传输数据。
同时，要根据实际情况进行飞行高度和速度的调整，以获取高质量的磁场数据。
五、数据处理与分析无人机航磁测量的核心是数据处理与分析。
首先，对采集到的原始数据进行预处理，包括数据的滤波、去噪和校正等操作，以提高数据的质量和准确性。
然后，根据测区的特点和测量目的，选择合适的数据处理方法，如数据插值、反演等，以获得地下磁场的空间分布和特征。
最后，对处理后的数据进行分析和解释，得出相应的结论和成果。
六、质量控制与评估在无人机航磁测量过程中，需要进行质量控制与评估，以确保数据的准确性和可靠性。
首先，对数据采集设备进行定期的检查和校准，以保证数据的一致性和稳定性。
详细讲解三次样条插值法及其实现方法
样条函数的定义 定义4.1 设区间[a,b]上给定一个节点划分
a=x0<x1<……<xn-1<xn=b 如果存在正整数k使得[a,b]上的分段函数s(x)满足 如下两条： (1)在[a,b]上有直到k-1阶连续导数。 (2)在每个小区间[xi,xi+1]上是次数不大于k的多项式。 则称分段函数s(x)是以(2.6)为节点集的k次样条函数。
x xi i i 1 hi
) mi1hi
( ) xi1x
1 hi
x [xi , xi1], hi xi1 xi , i 0,1,, n 1
(x) (2x 1)( x 1)2,1(x) x(x 1)2 13
对Si (x)求二阶导数 ,并整理后得
Si( x)
6( xi
xi 1 hi3
2x)
因为分段三次Hermite插值多项式已经至少是一阶连续 可导了，为了让它成为三次样条函数只需确定节点处 的一阶导数使这些节点处的二阶导数连续即可！
S(xi 0) S(xi 0), i 1,, n 1
S(x)
y ( xxi i 0 hi
)
y ( ) m h ( xi1x i1 0 hi
( yn
yn 1 )
2 hn1
(mn1
2mn )
立即可得下式：
21
其中：
nm1 nmn1 2mn gn
n
h0
h0 hn1
, n
hn1 h0 hn1
1 n
gn
3 n
y1 y0 h0
n
yn
yn1 hn1
联合基本方程得一个广义三对角或周期三对角方程组：
2 1
1
1
2
航空磁测数据常用转换处理方法
246管理及其他M anagement and other航空磁测数据常用转换处理方法田 亮，李亚南（中国冶金地质总局地球物理勘查院，河北 保定 071051）摘 要：航磁△T 原始数据是由各种不同空间位置、不同形态的磁性地质体磁场信息综合叠加反映的结果。
为了更好的提取有用的磁场信息，更好的进行磁异常的地质解释，航磁数据通常会采用不同的处理方法。
本文以北半球中高纬度某研究区为例，阐述了航空磁测数据常用的位场转换处理方法，如化极、上延、剩余异常提取、垂向和水平方向一阶导数、解析信号位等，以及应用此方法的目的和意义。
关键词：航磁测量；数据处理；化极；剩余异常；地质解释中图分类号：P631.2 文献标识码：A 文章编号：11-5004（2020）22-0246-2收稿日期：2020-11作者简介：田亮，男，生于1986年，汉族，河北唐县人，本科，工程师，研究方向：航空物探。
航空物探测量所得到的数据是地下所有地质体的综合反映，△T 原始数据是由各种不同空间位置、不同形态的磁性地质体磁场信息综合叠加反映的结果。
为了更好提取有用的磁场信息，提高磁异常解释地质效果，根据研究区域航磁异常和地质特点及地质解释的需要，使用地球物理专业软件，做平面数据处理，已达到解决不同的地质问题的效果。
常用的航磁△T 测量原始数据处理方法有原观测面化极、上延、剩余异常提取、垂向和水平方向一阶导数、解析信号等几种位场转换处理[1]。
下面以北半球中纬度某区域为例进行阐述。
1 航磁△T化极处理本区处于北半球中高纬度地区，由于倾斜磁化的影响，可能造成磁异常正值范围不是正好对应磁性地质体的正上方，而是相对于磁性地质体位置向南产生一定的偏移，这给磁性地质体的地面位置、形态及范围的确定均带来一定影响。
为了消除倾斜磁化对异常造成的影响，将实测的斜磁化△T 化到垂直磁化的垂直分量磁异常，简称化极，这样磁异常与场源的空间位置关系更为直观，有利于准确的确定异常场源的位置[2]，提高异常源的定位精度。
三维插值算法原理
三维插值算法原理引言：在地理信息系统（GIS）、计算机图形学和计算机辅助设计（CAD）等领域中，三维插值算法被广泛应用于空间数据的处理和分析。
三维插值算法可以通过已知的离散点数据，推断出未知位置的数值，从而生成连续、平滑的三维表面模型。
本文将介绍三维插值算法的原理及应用。
一、三维插值算法的基本原理三维插值算法的基本原理是通过已知的离散点数据，推断出未知位置的数值。
其关键在于确定未知位置的数值是通过哪些已知点来计算得出的，以及如何计算这些数值。
1.1 插值方法在三维插值算法中，常用的插值方法包括最邻近插值、反距离加权插值、克里金插值等。
最邻近插值方法简单直观，将未知位置的数值设置为与其最近的已知点的数值。
反距离加权插值方法则通过计算未知位置与已知点之间的距离，并根据距离的倒数进行加权平均得到未知位置的数值。
克里金插值方法则基于地理学中的克里金变异函数理论，利用空间相似性进行插值计算。
1.2 插值计算在三维插值算法中，插值计算的关键是确定插值点的权重。
权重的计算可以使用不同的函数或公式，例如高斯函数、指数函数等。
插值计算的结果将决定插值点的数值。
二、三维插值算法的应用三维插值算法在地理信息系统（GIS）、计算机图形学和计算机辅助设计（CAD）等领域中有广泛的应用。
以下是三维插值算法的一些典型应用案例。
2.1 地形表面模拟在地理信息系统（GIS）中，通过已知的地形高程数据点，可以使用三维插值算法生成地形表面模型。
这对于地形分析、地形可视化以及地形相关的工程应用非常重要。
2.2 气象数据插值在气象学中，三维插值算法可以用于推断空间上未观测到的气象数据，例如温度、湿度等。
通过已知的气象观测点数据，可以使用三维插值算法生成连续的气象场。
2.3 遥感图像处理在遥感图像处理中，三维插值算法可以用于填充或恢复图像中的缺失区域。
通过已知的图像数据点，可以使用三维插值算法推断出缺失区域的像素值。
2.4 计算机辅助设计在计算机辅助设计（CAD）中，三维插值算法可以用于生成平滑的曲面。
测控设备引导跟踪数据插值方法
测控设备引导跟踪数据插值方法庞岳峰;吴小东;牛攀峰【摘要】航天测控设备在引导跟踪时需要将转换后的方位、俯仰角度进行插值,在工程中尽可能采用易软件实现且不影响插值精度的插值方法.文中在介绍目前常用的Lagrange插值、Newton插值、Neville插值和Aitken插值4种方法原理的基础上,分析了4种插值下待插值点位置对插值结果的影响,通过实际算例得到结论.并讨论了4种方法在测控设备引导跟踪数据插值方面的优缺点,得出在航天测控设备引导跟踪插值时使用Neville算法最好的结论.【期刊名称】《电子科技》【年(卷),期】2016(029)011【总页数】4页(P118-121)【关键词】测控设备;引导跟踪;Neville插值;Lagrange插值【作者】庞岳峰;吴小东;牛攀峰【作者单位】酒泉卫星发射中心指挥控制站,甘肃酒泉732750;酒泉卫星发射中心指挥控制站,甘肃酒泉732750;酒泉卫星发射中心指挥控制站,甘肃酒泉732750【正文语种】中文【中图分类】TP391.9航天测控设备在航天器发射任务中主要承担目标跟踪、测量、控制任务，程序引导跟踪是测控设备的跟踪方式之一[1]。
程序引导跟踪原理是事先将目标理论弹道装订于测控设备跟踪单元，收到起飞信号后跟踪单元依据理论弹道转换为设备跟踪的A-E-R(方位-俯仰-距离)，最后控制天线按此A-E角运转。
目前测控设备装订的弹道都是1组/s，设备在进行引导时采用Lagrange算法插值为20组/s或者40组/s。
插值点在已知节点中间的称为内插，在节点外面的称为外推。
插值的方法有多种，目前常用的有Lagrange插值、Newton插值、Neville插值和Aitken插值等[2-7]，这些方法在GPS卫星精密星历差值方面的研究较为广泛[8-10]，在航天测控设备引导跟踪数据的插值研究较少。
1.1 Lagrange插值在数值计算方面，关于多项式插值的方法有多种，其中拉格朗日插值法是既简单而又实用的方法。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
航磁数据三维插值的实现和插值方法优选摘要：以铜陵铜官山矿区为试验区，研究了航磁数据在Surpac软件中不同三维插值方法的实现。
并针对不同方法插值的结果，从统计学角度正向分析和矿产预测角度逆向分析铜官山地区航磁数据最佳三维插值方法，以插值结果的均方差和磁有利成矿区间内矿体块数作为评价因子，从正逆双向进行铜官山矿区三维航磁插值方法的优选，最终得出距离幂次反比法是铜官山地区三维航磁数据插值的最佳方法的结论。
Abstract：关键词：三维插值，SURPAC，有利成矿，铜官山Keywords:1.研究区和研究数据介绍航空物探是一种最为快速和最有效的地质调查和地质找矿的方法，处理地质勘查中采集的航磁数据。
地球物理工作以及地质勘探中，为确保可实现性和可操作性，往往对研究区进行离散观测，然后对局部离散的数据经过一系列插值方法的处理以形成全区范围内时间或空间上离散数据体。
由于采样数据特点不同，各种插值方法的适用性也存在差异。
有效的插值方法是保证插值后数据准确性的关键，插值方法的不同直接影响到成矿预测结果的正确性与合理性。
二维空间的插值方法就是对研究区在平面上进行网格划分，再利用已知点值进行二维平面插值的方法。
而三维空间下的插值主要是对研究区进行立方块划分，再进行三维空间插值计算。
在在地球物理研究领域，采集到的数据通常是反映的是立体空间上的地质特征，数据的形式为（X,Y,Z,属性1，属性2，……），而目前对采集到的局部三维地质体数据进行插值处理的方法通常是将三维转化成二维，即设定在某一高程一定的情况下进行逐高程面分层插值，其数据形式为[(Z1,(X,Y,属性1，属性2，……)）,(Z2，(X,Y,属性1，属性2，……))，……（Zx，（X,Y,……））]。
文章选取以发育典型的矽卡岩型铜矿床而著称的安徽铜陵铜官山矿区作为试验区来研究处理该地区空间航磁数据最有效的插值方法。
铜陵铜官山矿区位于铜陵矿集区的西部，是中国东部重要的铜多金属矿产基地,已经发现大中型矿床20余处。
主要有天马山金（硫）矿床、铜官山铜（硫）矿床和金口岭铜（金）矿床等。
其范围为X：582807.553-587207.553，Y：3419535.407-3422635.407，Z：-3000-600。
2.三维航磁数据插值方法介绍根据试验区的范围建立的块体模型其最小单元尺寸为100×100×50，赋矿体为其属性。
根据铜陵矿集区航磁勘探采样布置情况，航磁数据的采样点坐标信息空间三向间隔均为200m，因此根据采样点的布置，我们在试验区的范围内建立最小单元尺寸为200×200×200的航磁块体模型，赋航磁值为其属性。
按照最小单元块尺寸对分别对试验区地质体模型和磁模型进行网格划分：将矿区地质体模型划分为98208块，其中试验区内矿体块数为537块，将矿区航磁模型划分为5888块，每一单元块体都对应每一个采样点的航磁值。
研究对选取的铜官山试验区的插值计算，就是将5888块的航磁块体模型单元块体的属性值采用最有效的插值方法赋值到98208块试验区地质体模型中去，使试验区98208块单元块体都对应一个航磁属性值。
常用的地统计插值方法有很多：反距离权重法、径向基函数插值法、距离幂次反比法、最邻近法和克里格插值法等，但Surpac 中用于块体模型插值的方法只有距离幂次反比法、最邻近法和克里金法。
另外，由于文章中选取的试验区块体模型最小单元尺寸与航磁块体模型最小单元的尺寸大小呈现比例关系，于是通过单元格的转换也可以实现由局部的航磁离散数据通过插值计算从而得到空间上离散数据体。
文章将主要分析距离幂次反比法（2次）、最邻近法和单元格转换法在应用到铜官山地区航磁数据插值处理上的有效性，具体插值方法原理介绍如下：2.1 最近距离法插值最近邻点插值法（NN ），又称泰森多边形方法，最近邻插值法是局部插值法中最为基础的方法，它以空间距离为基础，选取距离插值点最近的采样点，将其属性值赋予待插值点。
即插值点的属性值使用距离它最近的节点属性值表示：{}0111k 01231123123min 1123+11230,1,2k 123min ,,,,0,1,2k(,,p ),(,,p )(,,p )(,,p ),(,,p )(,,p )min ,i i i j j i j i j i j i k j A A A A A A i A p p p A p p p A p p p T A p p p A p p p A p p p T S S S S S S -+-===⎧⎫⎪⎪=⎨⎬⎪⎪⎩⎭=……………………………………,……………………，………… 式中,123(,,p )j p p p ……表示0,1,2k i A =……的属性值，0,1,2k i A i S =……表示其余点到i A 的距离值，min T 为其余点到i A 的距离值的最小值。
2.2距离幂次反比法距离幂次反比(IDW )法和地质统计学的插值方法一样,其方法原理是建立在区域化变量理论的基本原理基础之上的,即假设待插值的空间点属性值在一定的研究范围内具有相关性。
IDW 法的基本原理就是:假定区域化变量之间存在相关性并且这种相关性可以定量地表示为:样点与待估点之间的距离的幂次成反比,其公式如下：1()()ni i i Z B Z x λ==∑式中, ()Z B 为待估点的属性值；()i Z x 为已知的采样点的属性值; i λ为已知点的权重。
依据IDW 法的基本思想,确定权重i λ的方法为：111/n i k k i i i d d λ=⎡⎤⎡⎤=⎢⎥⎢⎥⎣⎦⎣⎦∑ 式中, i d 为待估点与已知点之间的距离;ｋ作为i d 的幂指数,其取值由具体的研究情况确定,通常它可以取1,2,3等整数。
IDW 法在利用样本信息对未采样点进行估值的过程中,比较充分、恰当地利用了样点与待估点的距离即采样点与待估点之间的空间位置信息,符合数据点属性分布的一般规律。
2.3单元格转换法在对研究区进行立方块划分中，通常采样点的准确坐标对应为最小单元块的质心坐标。
因此在勘查布置中，研究区块体模型的尺寸和根据勘查采样点分布构建的块体模型的尺寸不匹配是勘查数据有效运用普遍面临的难题，单元格分解法是最简单并且应用范围十分狭隘的一种插值方法。
当研究区的块体模型尺寸大小与勘查布置中采样点数据尺寸大小呈整数倍数比例时，几个较小尺寸的单元格组合后的尺寸刚好与较大尺寸的单元格吻合，这样就可以将较大尺寸的单元格属性赋给组合之后能与之尺寸匹配的若干小单元格，实现从局部区域到全区数据的插值转换。
比如：文章中试验区的地质体块体模型和航磁块体模型的范围是一致的，航磁块体模型的最小单元块数的尺寸是200×200×200，而地质体块体模型的最小单元尺寸是100×100×50，即某一位置上相邻的16块（2块×2块×4块）地质体模型的小块体就组合形成相同位置上200×200×200航磁块体模型最小单元块体，就可以将该块航磁模型块体的属性值赋给相同位置上16块地质体单元块体，使16块小尺寸的地质体块体模型的属性与较大尺寸的航磁块体单元属性保持一致。
即通过这种单元格转换法将局部的航磁属性值插值计算到试验区全区中去。
单元格转换法原理简单，其表达式如下：11112111(,,)(,,)(,,)(,,)i i i i A x y z a x y z a x y z a x y z =++若：…… 12j 112j 212j i 12j p (,,)(,,)(,,)=(,,)A a a a ==p p p p p p p p p p p 则：………………………… 式中，(,,)A x y z 表示的是以(,,)x y z 坐标为质心的较大尺寸立方块单元，(,,)i i i i a x y z 表示是以(,,)i i i x y z 为坐标为质心的较小尺寸立方块单元，12j (,,)p p p ……为立方块单元的属性值。
3.Surpac 中航磁数据插值的实现Surpac 软件支持下的以上块体模型插值方法主要运用于研究区内矿体品位的插值，对局部的矿体品位数据选取合理的插值方法进行插值计算，从而获得接近真实情况的全区矿体品位数据，再得以计算出准确的资源储量。
文章参照surpac 软件中矿体模型品位插值的方法对试验区的地质体块体模型进行航磁插值计算。
在进行矿体模型品位插值计算时，矿体的品位都是正值。
而试验区的航磁值范围是-100-300，航磁块体模型中存在大量的单元块体其属性值是负值。
为了消除负值的航磁数据对插值结果的影响，在对航磁采样点数据进行预处理时，我们将全区所有的航磁数据加上+100，将区内所有的航磁负值全部转换为正值进行插值计算，最后将插值计算的结果再同时减去100，即可得到符合实际情况的全区航磁插值数据。
为了对比铜官山地区航磁数据不同插值方法的优越性，必须确保各种插值方法中相关参数（搜索类型、样品数、搜索半径、各向异性比率）的一致性。
在Surpac 软件中，对建立好的试验区地质体块体模型新建三个“属性”（单元格转换法磁值、最邻近距离法磁值、距离幂次反比法磁值）。
对于最邻近距离法、距离幂次反比法插值，需要在Surpac 软件中将建立的最小单元尺寸为200×200×200航磁块体模型“导出”为通用的.csv 格式的文本数据格式，其数据格式为：（X,Y,Z ，磁值），每个坐标位置的磁值为以此坐标为质心的块体属性值，即地质勘探布置时航磁采样点的磁值。
以导出的.csv 文件为基础再“导入”Surpac 中存为.str 格式的线串文件作为矿区地质体模型的插值线要素。
然后在试验区块体模型中选入“.str格式线要素”进行分别进行插值计算，插值的结果将分别存在试验区模型的“最邻近距离法磁值”属性和“距离幂次反比法磁值”属性之中。
对于单元格转换法，则按照单元格转换法的原理，参照航磁块体模型和试验区地质体块体模型的最小单元格尺寸大小，将航磁块体模型某一位置上的属性值赋给相同位置上16块试验区地质体单元块体。
使这16块小尺寸的地质体块体模型的属性与较大尺寸的航磁块体单元属性保持一致。
4.Surpac中航磁数据插值方法优选不同的插值方法对应不同的计算原理，其方法的优越性也因应用的对象不同而存在差异。
在地质勘查、矿产预测等地质工作中，应用不同方法插值后数据的准确性将直接影响到对地质认识的正确性。
对以上三种方法插值的结果进行分析，对单元格转换法、最邻近距离法、距离幂次反比法三种插值方法在铜陵铜官山矿区三维航磁数据处理中的有效性和优越性作出评价，以分析出适用于该地区三维航磁数据处理最优越的插值方法。