导读：各位师兄师弟，本文档仅供学习交流，请不要过多传阅，因为有的地方不一定适合实际情况，如有不对的地方请指正。其中，数据处理中遇到的软件MAG3.0可以用其他软件（师弟们可以用RGIS）代替，比如在日变改正的时候，平滑滤波和计算正常场的时候等等，但结果是相似的。工作中，要注意细节，比如数据的顺序千万要注意，其次，实际测量可能存在跳点漏点，这些点要单独处理。再就是一些“异常”点的处理，比如受房屋电线干扰后异常突变的点最后一定要有相应处理。最后祝各位师弟将李老师的5507发扬光大，不辜负李老师的教导。汪洋1：导数据本次磁测数据是适用捷克PMG-1质子磁力仪采集，坐标数据采集采用的是集思宝GPS，各自数据的导出见各自说明书。2：日变数据处理圆滑，将日变原始数据“2012-7-21.dat”拷贝一份出来，修改名称为“y2012-7-21.dat”，打开这个数据文件，删除除除测得场值外的所有行列，如图：，打开mag3.0，选择如图项进行圆滑，如图：其中，分解尺度表示圆滑滤波程度，一般填3-5，滤波结果如图，其中，逼近文件“y2012-7-21A.dat”就是我们需要的滤波结果，数据如图所示：打开“y2012-7-21A.dat”和“2012-7-21.dat”，用“y2012-7-21A.dat”的数据将“2012-7-21.dat”中的数据列替换掉，重命名为“yh2012-7-21.dat”，其数据如图:这就是所要的圆滑后的日变文件。3：日变改正打开mag3.0，如图改正如下其中，文件格式按所用的磁力仪型号选择，因本次所用磁力仪是捷克PMG-1，故而选PMG，总场T0值为分基点T0值（分基点所测值的算术平均值）。改正以后结果文件为“0261_C”，数据格式如下：其中J列grad即为改正以后的结果。将上面的结果另存为xls，将上面的J列减去总基点T0值，即为日变改正后的结果（L列）。4：数据整理将上面的表格进行整理。首先将pos列修改成磁测设计坐标值（一般写后四位，注意坐标值一定要与后面的实测值对应，不能反了，其次要注意实际测量中如果有跳点，一定要将所跳的点的位置空出来）如图：其中，5、6两行即为跳点的缺失项。5：填表建立如图1所示的表格（图示为现成的表格）：图1将上面的数据按设计坐标顺序填入表格（如果上面步骤4里的数据顺序不对，经调整以后填入），当前主要能填的是原始数据、日变改正值。6：GPS数据处理（以集思宝GPS为例）将导出的1830线的GPS数据用集思宝的软件GIS-office打开，选择文件-导出，弹出对话框如图：设置也如图，点OK，得到结果文件，打开，结果如下在表格中，将B、C两列的字母N和E替换掉，只要数字，如图同样，将B、C、D三列填入图一所示表格实测x、实测y和高程列。注意：1，同上，请注意填入坐标的顺序，应该与设计坐标顺序一致；2，从GPS导出的x坐标小数点前仅6位，而实际坐标应该为8位，如499100应该为34499100，故而写入图一实测x列时应该加上34000000，即499100+34000000=34499100。最终结果如图一。7：计算正常场写如图所示的一个数据文件“总基点经纬度.dat”其中A表示野外总基点测量数据时间，B、C表示总基点经纬度，D表示总基点高程。打开mag3.0，选择，设置如图计算所得正常场文件“总基点经纬度A.dat”文件如图：其中X、Y、Z表示各向分量，F表示正常场，I为倾角，D为偏角。8：高度改正和水平梯度改正高度改正：高度改正系数∂T0∂R=-3T0R，其中T0即为上面算出来的总基点正常场，R为地球半径=6371000m。水平梯度改正：利用上面所得X、Y、Z分量，求得H2=X2+Y2，T02=H2+Z2，而梯度改正系数即为∂T0∂x=3ZH2RT0。高度改正如图：其中，高程值G列为实测高程值，改正公式为：改正值=（实测高程值-总基点高程值）*高度改正系数，此处为（1930-1796）*0.002265=0.30。水平梯度改正如图：其中，改正公式为：改正值（实测y值-总基点y值）*梯度改正系数，此处为（2911831-2910635）*0.005591=-6.6840405。9：△T值的求取将上面图一所得的3列相加，即可得△T值。数据预处理完毕。