本发明专利技术涉及一种基于SEGY文件直接重构SEGY数据坐标的方法及处理终端，包括：步骤1：获得SEGY文件，SEGY文件中存储有各个测点的大地坐标；步骤2：将各个测点的大地坐标转换为平面坐标；步骤3：根据当前测点的坐标和上一个测点的坐标确定当前测点对应的调查船方位角A；步骤4：根据当前测点对应的调查船方位角A，以及反射点相对于当前测点的偏移量，得到反射点的平面坐标，将反射点的坐标作为测点的校正后的平面坐标；步骤5：将测点的校正后的平面坐标转换为和步骤1相同格式的大地坐标，完成重构SEGY数据坐标。本发明专利技术无需导航文件即可对SEGY数据坐标校正，占用资源更少，易用性和实用性高。高。高。
【技术实现步骤摘要】
一种基于SEGY文件直接重构SEGY数据坐标的方法及处理终端
[0001]本专利技术涉及SEGY数据坐标校正
，具体涉及一种基于SEGY文件直接重构SEGY数据坐标的方法及处理终端。
技术介绍
[0002]目前，实际测量中浅地层剖面(简称浅剖)、单道地震记录的坐标为导航定位设备传输的位置，此位置是GPS天线头或导航参考点的位置，并不是信号反射点(简称反射点)的位置。由于记录的坐标位置不准，存在误差，导致后续对地震剖面解释带来了困难，造成根据剖面推断的地下地质情况与实际情况不符，使得给后续应用造成麻烦甚至无法使用。对于固定在船上的浅剖和拖缆单道地震设备来说，反射点位置由设备相对位置和船首向决定，通常反射点与GPS天线头之间存在几十米的偏移量，对于海砂调查、井场调查、海底管线探测等高精度海底探测来说，几十米的偏移量校正也是必不可少的。通常浅剖和单道地震的数据存储格式为SEGY格式，现有的SEGY数据处理软件，如OMEGA、SU、RadExpro等，不能进行偏移量校正。
[0003]而在目前对SEGY数据坐标进行校正方法中，大都是需要利用导航文件，所需数据种类较多，操作过程较繁杂。在具体校正操作过程中，需要将导航读入内存，内存占用大，对大体积量的数据处理不利。更重要的是，实际数据测量过程中，浅剖或单道比导航上线早、下线晚，造成导航文件中未记录测线开始和结束的部分位置，对于这部分缺失导航信息的位置，还是需要特殊处理。更近一步，对于很多使用浅剖或单道数据的人来说，未必能获得导航文件，对于缺失导航文件的情况现有方法无法进行坐标偏移量校正。因此，需要一种不依赖于导航文件甚至不依赖任何其他文件能够对SEGY数据坐标进行校正的技术。
[0004]相关的参考文献如下：
[0005]文献1：公布号为CN 111078953 A、专利技术名称为“基于导航文件直接重构SEGY数据坐标的方法”的中国申请专利技术专利。
技术实现思路
[0006]针对现有技术的不足，本专利技术的目的之一提供一种基于SEGY文件直接重构SEGY数据坐标的方法，其能够解决不依赖其他文件数据对SEGY数据坐标校正的问题；
[0007]本专利技术的目的之二提供一种处理终端，其能够解决不依赖其他文件数据对SEGY数据坐标校正的问题。
[0008]实现本专利技术的目的之一的技术方案为：一种基于SEGY文件直接重构SEGY数据坐标的方法，包括如下步骤：
[0009]步骤1：获得SEGY文件，SEGY文件中存储有各个测点的大地坐标；
[0010]步骤2：将各个测点的大地坐标转换为平面坐标，其中，第i个测点的平面坐标记为(x
i
,y
i
)，i≥1，第i+1个测点的平面坐标记为(x
i+1
,y
i+1
)；
[0011]步骤3：根据当前测点的坐标和上一个测点的坐标确定当前测点对应的调查船方
位角A；
[0012]步骤4：根据当前测点对应的调查船方位角A，以及反射点相对于当前测点的偏移量，得到反射点的平面坐标，将反射点的平面坐标作为当前测点的平面坐标，从而得到各个测点的校正后的平面坐标；
[0013]步骤5：将各个测点的校正后的平面坐标转换为和步骤1中测点相同格式的大地坐标，从而得到校正后SEGY坐标数据，完成重构SEGY数据坐标。
[0014]进一步地，所述步骤1中，若SEGY文件中存在异常点或重复坐标的测点，则对异常点或重复坐标进行处理，以使得SEGY数据中不存在异常点和重复坐标后的测点后，再执行步骤2。
[0015]进一步地，所述步骤2中，采用UTM投影将各个测点的大地坐标转换为平面坐标。
[0016]进一步地，所述根据当前测点的坐标和上一个测点的坐标确定当前测点对应的调查船方位角A，其具体实现过程包括：
[0017]当y
i+1
>y
i
时，
[0018]若x
i+1
≥x
i
，则A
i+1
＝α，否则，A
i+1
＝
‑
α，α表示调查船中心轴线与Y轴之间的夹角，A
i+1
表示第i+1个测点对应的调查船方位角，
[0019]当y
i+1
＝y
i
时，
[0020]若x
i+1
≥x
i
，则A
i+1
＝90
°
，否则，A
i+1
＝270
°
。
[0021]当y
i+1
<y
i
时，
[0022]若x
i+1
≥x
i
，则A
i+1
＝180
°‑
α，否则，A
i+1
＝180
°
+α
[0023]其中，A2＝A1。
[0024]进一步地，所述根据当前测点对应的调查船方位角A，以及反射点相对于当前测点的偏移量，得到反射点的坐标，将反射点的坐标作为当前测点的平面坐标，从而得到各个测点的校正后的平面坐标，其具体实现过程包括：
[0025]按公式
①
计算得到反射点的平面坐标：
[0026][0027]式中，(x
′
i+1
,y
′
i+1
)表示第i+1个测点的校正后的平面坐标，(dx,dy)表示反射点相对测点的偏移量。
[0028]实现本专利技术的目的之二的技术方案为：一种处理终端，其包括：
[0029]存储器，用于存储程序指令；
[0030]处理器，用于运行所述程序指令，以执行所述于SEGY文件直接重构SEGY数据坐标的方法的步骤。
[0031]本专利技术的有益效果为：本专利技术解决了在缺失导航文件的情况下，可对浅剖或单道数据进行坐标偏移量校正的问题，在实际操作中，由于不需要将其他导航文件等其他数据导入，不需要占用过多的内存，效率更高，占用资源更少。另外，通过对比，本专利技术计算结果和其他方法计算结果基本相同，但本方法使用更少的数据类型和数据量，操作步骤更简单，计算效率更高。以及支持不同道长、不同采样率的多文件批量处理，提高了本专利技术的易用性
和实用性。
附图说明
[0032]图1为较佳实施例的流程示意图；
[0033]图2是反射点与GPS天线头相对位置示意图；
[0034]图3为采用本专利技术和其他现有方法得到的重构SEGY数据坐标和原始测得位置对比示意图；
[0035]图4为使用SEGY处理软件查询到的本方法处理结果的坐标校正后SEGY文件的道头信息示意图；
[0036]图5为处理终端的示意图。
具体实施方式
[0037]下面，结合附图以及具体实施方案，对本专利技术做进一步描述。<本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种基于SEGY文件直接重构SEGY数据坐标的方法，其特征在于，包括如下步骤：步骤1：获得SEGY文件，SEGY文件中存储有各个测点的大地坐标；步骤2：将各个测点的大地坐标转换为平面坐标，其中，第i个测点的平面坐标记为(x
i
,y
i
)，i≥1，第i+1个测点的平面坐标记为(x
i+1
,y
i+1
)；步骤3：根据当前测点的坐标和上一个测点的坐标确定当前测点对应的调查船方位角A；步骤4：根据当前测点对应的调查船方位角A，以及反射点相对于当前测点的偏移量，得到反射点的平面坐标，将反射点的平面坐标作为当前测点的平面坐标，从而得到各个测点的校正后的平面坐标；步骤5：将各个测点的校正后的平面坐标转换为和步骤1中测点相同格式的大地坐标，从而得到校正后SEGY坐标数据，完成重构SEGY数据坐标。2.根据权利要求1所述的基于SEGY文件直接重构SEGY数据坐标的方法，其特征在于，所述步骤1中，若SEGY文件中存在异常点或重复坐标的测点，则对异常点或重复坐标进行处理，以使得SEGY数据中不存在异常点和重复坐标的测点后，再执行步骤2。3.根据权利要求1所述的基于SEGY文件直接重构SEGY数据坐标的方法，其特征在于，所述步骤2中，采用UTM投影将各个测点的大地坐标转换为平面坐标。4.根据权利要求1所述的基于SEGY文件直接重构SEGY数据坐标的方法，其特征在于，所述根据当前测点的坐标和上一个测点的坐标确定当前测点对应的调查船方位角A，其具体实现过程包括：当y
i+1
>y
i
时，若x
i+1
≥x
i
，则A
i+1
＝α，否则，A
...
【专利技术属性】
技术研发人员：徐岩，侯建梅，曾凡祥，邓贤泽，
申请(专利权)人：南方海洋科学与工程广东省实验室广州，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。