1 算法

算法来自于GEMLink 5.2的帮助文档，这个文档基本解决了算法问题。

GemLink日变改正模块界面

1.1 概述

日变改正模块旨在执行磁力日变数据计算，而不用在仪器上进行日变（仪器是未经过校正的原始数据）。这个模块要求已经保存了Base和Rover文件，结果将另存为另外的文件，其中包含Rover文件的所有信息，以及附加在该行末尾的校正磁场值。

1.2 文件格式：

文件必须是ASCII文本文件。每行文本的末尾必须有一个CR(回车)或CRLF(回车/换行)。仅以LF终止的行不能使用，因为它们被视为一行。

字段或列之间的默认分隔符是SPACE。可以从“文件菜单”下的“字段分隔符”菜单中选择不同的分隔符。这个设置是临时的，每次GEMLinkW打开时，空格都是默认的分隔符。

注意，Base和Rover文件必须具有相同的分隔符。如果选择了不正确的分隔符，系统可能会产生意想不到的结果或错误。

1.3 时间戳：

自从GEMlinkW 2.4发布以来，日变改正模块在时间戳格式方面提供了更大的灵活性。Base和Rover文件接受以下24小时的时间戳格式:

HHMMSS;
HHMMSS.s;
HHMMSS.ss;
HHMMSS.sss;
HH: MM: SS;
HH: MM: SS.s;
HH: MM: SS;
HH:MM: ss.sss

不符合这些格式的条目将被拒绝。

此外，如果发现240000.0是表示午夜的时间戳，则将其视为000000.0。分和秒不能超过59。匹配这些格式，但hour、Minutes或Seconds值非法的条目也会被拒绝。

1.4 时间转换

为了便于比较时间值，所有时间都通过简单的小时、分钟、秒和秒的分数加权相加转换为十分之一秒(1/10秒)。

虽然小数点后三位的时间戳也被接受，但时间只与十分之一秒进行比较。小数点后第2位和第3位将被丢弃，不进行舍入。

时间戳0 =< |日变站时间-移动站时间|<= 0.099秒被认为是相等的。

1.5 同步问题

当然，为了进行日改，日变站和移动站必须在测量前进行时间同步。如果两个设备都配备了GPS，它们会自动同步到UTC，因此它们的读取周期由GPS定时控制。如果只有一台设备配备了GPS，您可能希望将该设备同步到UTC，然后使用同步电缆将UTC时间传输到另一台设备。如果没有一台设备配备GPS，在其中一台设备上手动设置时间，并使用同步电缆将时间传输到另一台设备上。

事实上，野外都是手动进行同步的，而且同步引起的误差很小。

1.6 日期

本程序不检查文件的日期。用户自己确保两个文件都是同一天就行。

文件时间不能超过午夜，因为程序并没有识别这种情况的功能。

1.7 未修正行的包含

如果在“日变改正”工具对话框中选中“包括未修正的行”(默认)，系统将存储无法修正的数据行而不作任何更改。结果是一个包含移动站所有读数的文件，然而，只有经过校正的行才会显示在行末附加的日改校正的磁场值。

如果选中“舍弃未更正的行”选项，则结果文件只包含已成功更正的行。

1.8 使用说明

打开Base或Rover的文件后，会弹出“设置”对话框。对话框里会显示文件的前50行，以验证是否为正确的文件。

然后，程序会检查所打开的文件，确定首行和尾行。首行定义为长度相同的3行中的第一行。尾行被定义为与首行具有相同长度和数量的字段分隔符的最接近文件末尾的行。

第一行数据将根据当前字段分隔符设置进行解析，并将值传递给Time和Magnetic field下拉菜单。下拉框的值用于确定行中每个字段的位置及其格式。

除了作为选项(来自首行数据)显示的值之外，不可能选择其他值。如果遇到其他值，您将看到一条错误消息。如果未使用选定的字段分隔符分隔行字段，则日变改正计算可能发生错误。

可以任意在Rover或Base设置窗口中输入基准值datum。

当准备好开始计算时(即当两个文件都打开时)，在“日变改正”对话框中按下“correct”按钮。

程序首先检查文件以确定它们是否有相同的读数。如果两个文件从开始到结束的时间间隔不重叠，则无法进行更正，进程将停止。

同样，请记住，文件只根据时间进行比较，日期不被考虑在内。确保两个文件都是同一天。

从Rover文件的第一行数据中提取的时间与Base文件的时间进行比较。

如果移动站的时间不是一个共同的读数(即它的时间不在重叠的时间间隔内)，则根据'Include' 'Discard'选项的设置，数据线要么不受影响地传递，要么被丢弃。

如果时间经检验后没问题，则通过以下两种方式之一进行修正：

（1）如果在Base文件中找到了完全相同的时间，那么磁场修正的计算方式为:

改正值Corrected Value = 移动站值Rover Field - 日变站值Base Field + 基准值Datum

修正后的值被附加在移动站的数据行的末尾，这条线将被标记为“i——”

（2）如果时间不完全匹配，则从Base文件的最近时间读数线性插值Base值。修正后的值被计算为:

改正值Corrected Value = 移动站值Rover Field - 插值后的日变站值Interpolated Base Field + 基准值Datum

为衡量计算结果准确度，用于插值的基线的时间间隔(以秒为单位)也添加在该行的末尾。例如，从间隔5秒的基本读数插入的校正字段将被标记为“i005”。

如果时间间隔超过999秒，该行将被标记为“i+++”。

对移动站的所有行都执行相同的计算，直到所有行都被校正完。

Rover文件中的无效数据行(例如位置标签、注释、不包含有效时间格式或数值分隔符的行)将不做任何改动的传递给新文件。如果该行包含有效的时间，但其中一些值不能提取，则会发出警告，显示无效的行，并且该行将传递给标记为“Error”的最终文件。

GEMLinkW的以前版本对日变文件中的无效数据行容忍度较低。但是，从GEMLnkW 2.4开始，您现在会收到一条显示无效行消息的警告，该过程将继续使用下一个有效行。为了避免看到这些消息，我们建议您不要手动编辑这些文件(意思就是说，从仪器导出来的文件不要编辑，否则可能会编辑错而导致出现错误的改正信息)。

在文件被改正完后，一个摘要消息显示了以下几项内容的改正量: