CN116625240B - 一种组合式水下定位设备的标定方法 - Google Patents

一种组合式水下定位设备的标定方法 Download PDF

Info

Publication number
CN116625240B
CN116625240B CN202310894786.5A CN202310894786A CN116625240B CN 116625240 B CN116625240 B CN 116625240B CN 202310894786 A CN202310894786 A CN 202310894786A CN 116625240 B CN116625240 B CN 116625240B
Authority
CN
China
Prior art keywords
coordinate system
displacement meter
steps
coordinates
stay wire
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202310894786.5A
Other languages
English (en)
Other versions
CN116625240A (zh
Inventor
杨润来
锁旭宏
李允良
张乃受
成益品
管泽旭
王强
宁进进
于博
宋智强
安鸣赞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong Provincial Laboratory Of Artificial Intelligence And Digital Economy Shenzhen
CCCC First Harbor Engineering Co Ltd
No 2 Engineering Co Ltd of CCCC First Harbor Engineering Co Ltd
Original Assignee
Guangdong Provincial Laboratory Of Artificial Intelligence And Digital Economy Shenzhen
CCCC First Harbor Engineering Co Ltd
No 2 Engineering Co Ltd of CCCC First Harbor Engineering Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Guangdong Provincial Laboratory Of Artificial Intelligence And Digital Economy Shenzhen, CCCC First Harbor Engineering Co Ltd, No 2 Engineering Co Ltd of CCCC First Harbor Engineering Co Ltd filed Critical Guangdong Provincial Laboratory Of Artificial Intelligence And Digital Economy Shenzhen
Priority to CN202310894786.5A priority Critical patent/CN116625240B/zh
Publication of CN116625240A publication Critical patent/CN116625240A/zh
Application granted granted Critical
Publication of CN116625240B publication Critical patent/CN116625240B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/002Measuring arrangements characterised by the use of optical techniques for measuring two or more coordinates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01BMEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
    • G01B11/00Measuring arrangements characterised by the use of optical techniques
    • G01B11/26Measuring arrangements characterised by the use of optical techniques for measuring angles or tapers; for testing the alignment of axes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明属于定位设备标定技术领域,涉及一种组合式水下定位设备的标定方法。该方法包括:S1、建立基于仪器箱的坐标系、基于摄像位移计的坐标系、基于拉线位移计的坐标系;S2、设置参考点及多个标靶点;S3、用全站仪测量坐标系原点、于全站仪坐标系下坐标;用摄像位移计测量坐标系下坐标;用拉线位移计测量坐标系下坐标;S4、解算全站仪坐标系向坐标系的转换参数,以将于全站仪坐标系下坐标换算为于坐标系下坐标;S5、建立第一观测方程,解算坐标系向坐标系的转换参数;S6、建立第二观测方程,解算坐标系向坐标系的转换参数。本发明实现了定位设备的组合式设置及统一标定,提高了标定精度。

Description

一种组合式水下定位设备的标定方法
技术领域
本发明属于定位设备标定技术领域,具体涉及一种组合式水下定位设备的标定方法。
背景技术
沉管隧道施工中,常用的水下定位方法有测量塔法、声呐法等;但随着水下作业水深的逐渐增大,这些常规水下定位方法的定位精度逐渐降低,且成本较高。目前,拉线位移计、摄像位移计等新式定位装置已开始出现在沉管隧道施工的水下定位中,如CN216815272U公开的沉管隧道最终接头水下安装定位系统,但这些新式定位装置均独立安装使用,一方面多种定位装置分别安装难度较大、占用空间大、敷设线缆多,另一方面多种定位装置各自单独标定,标定误差较大。
发明内容
针对相关技术中存在的不足之处,本发明提供一种组合式水下定位设备的标定方法,旨在实现多种定位装置的组合式设置以及统一标定,提高组合式水下定位设备的标定精度。
本发明的组合式水下定位设备的标定方法,适用于组合式水下定位设备;组合式水下定位设备包括仪器箱、安装于仪器箱内的摄像位移计和拉线位移计;摄像位移计包括并排设置的两台相机;拉线位移计包括可转动的拉杆、贯穿拉杆且可伸缩的拉线、测定拉线长度的拉线编码器和测定拉线方向的激光测角装置;
组合式水下定位设备的标定方法包括以下步骤:
S1、分别建立基于仪器箱的坐标系、基于摄像位移计的/>坐标系、基于拉线位移计的/>坐标系;将/>坐标系的原点记作/>,/>位于仪器箱顶面的中心,将/>坐标系的主轴记作/>轴,/>轴朝向摄像位移计的拍摄方向;
S2、将仪器箱放置于地面上并使其保持水平状态;在轴上间隔设置两个参考点、/>,将其在/>坐标系下的坐标分别记作/>、/>;沿/>轴方向间隔设置多个标靶并在每一标靶上布设多个标靶点,对所有标靶点编号并记作/>,其中,/>或/>或…或/>
S3、架设全站仪以测量、/>、/>及/>于全站仪坐标系下的坐标,将其分别记作、/>、/>、/>;利用摄像位移计测量/>于/>坐标系下的坐标,将其记作/>;利用拉线位移计测量/>于/>坐标系下的坐标,将其记作/>
S4、通过、/>、/>于全站仪坐标系下的坐标和于/>坐标系下的坐标,解算全站仪坐标系向/>坐标系转换时的全站仪坐标系转换参数;根据全站仪坐标系转换参数,将/>于全站仪坐标系下的坐标/>换算为/>于/>坐标系下的坐标,将其记作
S5、利用和/>,建立标靶点/>的第一观测方程,以解算/>坐标系向/>坐标系转换时的摄像位移计坐标系转换参数/>
S6、利用和/>,建立标靶点/>的第二观测方程,以解算/>坐标系向/>坐标系转换时的拉线位移计坐标系转换参数/>
上述技术方案,解决了现有技术中多种定位装置独立安装使用、各自单独标定而导致标定误差较大的问题,实现了摄像位移计和拉线位移计的统一标定,提高了标定精度和标定效率,进而能够在工程应用中提高对待测物体的定位精度。
在其中一些实施例中,在步骤S4中,全站仪坐标系转换参数包括水平方向角和高度转换系数,全站仪坐标系转换参数的解算包括以下步骤:
S41、在全站仪坐标系下,将、/>、/>的坐标与/>坐标/>相减以归置零,得到/>、/>、/>归置零后的新坐标,将其分别记作/>、/>,此时/>
S42、计算水平方向角,/>
S43、根据式(1)计算转换到/>坐标系下的坐标/>,根据式(2)计算/>转换到/>坐标系下的坐标/>;根据式(3)计算第一高度转换系数/>和第二高度转换系数/>
(1);
(2);
(3)。
在其中一些实施例中,在步骤S4中,全站仪坐标系转换参数的解算包括以下步骤:
S44、全站仪坐标系转换精度验证:判断、/>是否位于/>范围内,/>为预设测量容差;判断/>是否等于/>、/>是否等于/>;判断/>和/>是否相等;若三项判断结果均为是,则验证通过;若三项判断结果中任一项为否,则验证未通过,重新架设全站仪测量/>、/>、/>及/>于全站仪坐标系下的坐标,然后重新执行步骤S4。
在其中一些实施例中,在步骤S4中,于/>坐标系下的坐标/>的解算包括以下步骤:
S45、在全站仪坐标系下,将的坐标/>与/>坐标/>相减以归置零,得到/>归置零后的新坐标,将其记作/>
S46、根据式(4)计算转换到/>坐标系下的坐标/>
(4)。
在其中一些实施例中,在步骤S1中,将坐标系的原点记作/>,/>位于摄像位移计的一相机的摄像中心,将/>坐标系的主轴记作/>轴,/>轴朝向该相机主光轴方向,将/>坐标系的横轴记作/>轴,/>轴朝向摄像位移计的另一相机的摄像中心;
在步骤S5中,摄像位移计坐标系转换参数的解算包括以下步骤:
S51、建立标靶点的第一观测方程,表示为式(5);
(5);
式(5)中,、/>、/>为/>坐标系与/>坐标系之间的平移参数;/>、/>、/>为/>坐标系与/>坐标系之间的旋转参数;/>、/>、/>为/>坐标系与/>坐标系之间三个坐标轴方向的尺度;/>,/>为第/>个标靶点的第一测量改正数;
,/>,/>,建立所有标靶点的第一观测误差方程,表示为式(6);
(6);
根据最小二乘法平差计算得到坐标系向/>坐标系转换时的摄像位移计坐标系转换参数/>,表示为式(7);
(7)。
在其中一些实施例中,在步骤S5中,摄像位移计坐标系转换参数的解算还包括以下步骤:
S52、将摄像位移计坐标系转换参数代入到式(6)中,得到/>;计算出/>的标准差/>
S53、根据三倍标准差法,剔除对应的标靶点数据;若剔除的标靶点数据量超过全部标靶点数据量的/>,则重新利用摄像位移计测量/>于/>坐标系下的坐标,然后重新执行步骤S5;若剔除的标靶点数据量未超过全部标靶点数据量的/>,则利用其它未被剔除的标靶点数据,重新执行步骤S51,以解算出修正后的摄像位移计坐标系转换参数
在其中一些实施例中,在步骤S1中,将坐标系的原点记作/>,/>位于拉线位移计的拉杆转动中心,将/>坐标系的主轴记作/>轴,/>轴朝向拉线位移计正前方;
在步骤S3中,在拉线位移计的测量数据中,/>为拉线编码器的读数,/>和/>分别为激光测角装置中光斑相机拍摄的激光光点于相片上像素点的横向位置和竖向位置;
在步骤S6中,拉线位移计坐标系转换参数的解算包括以下步骤:
S61、拉线位移计测量数据规正所需参数的计算,具体包括:
S611、在坐标系下,量取拉杆转动中心/>与/>坐标系原点/>之间的坐标差,该坐标差为常数,记作/>
S612、分别在拉线未拉出和最长拉出状态下读取拉线编码器的读数、/>,同时用卷尺量取拉线的拉出距离/>;根据式(8)计算拉线尺度/>
(8);
S613、在拉线未拉出时,量取拉杆转动中心至拉线的拉环内侧的距离/>;根据式(9)计算拉线读数常数/>
(9);
S614、将拉线拉出至接近最长距离,调准拉线方向使其与坐标系的/>轴保持平行,读取激光光点于相片上像素点的横向位置和纵向位置的读数,分别记作光线中心位置常数/>和/>
S615、在一标靶上选择同一水平位置上的两个水平标靶点,将该两个水平标靶点的间距记作,将拉杆转动中心/>到该标靶中心的距离记作/>,根据式(10)计算该两个水平标靶点处拉线的夹角/>;利用拉线位移计读取拉线位于该两个水平标靶点时的光斑相机读数,相减得到读数差/>,将光斑相机焦距记作/>、像素点大小记作/>,根据式(11)计算对应该读数差的夹角/>;然后根据式(12)计算横向张角系数/>
在该标靶上选择同一竖直位置上的两个竖直标靶点,将该两个竖直标靶点的间距记作,根据式(13)计算该两个竖直标靶点处拉线的夹角/>;利用拉线位移计读取拉线位于该两个竖直标靶点时的光斑相机读数,相减得到读数差/>,根据式(14)计算对应该读数差的夹角/>;然后根据式(15)计算纵向张角系数/>
(10);
(11);
(12);
(13);
(14);
(15);
S62、对拉线位移计的测量数据进行规正计算:
根据式(16)和式(17)计算出规正后的拉线位移计的测量数据,将其记作
(16);
(17)。
在其中一些实施例中,在步骤S61中,拉线位移计测量数据规正所需参数的计算还包括:
S616、拉线位移计规正参数验证:在另一位置处的标靶上选择同一水平位置上的两个水平标靶点和同一竖直位置上的两个竖直标靶点,参考步骤S615,重新计算和/>;判断/>和/>的两次计算结果的差值是否小于本次计算数值的20%,若是,则验证通过,若否,则验证未通过,重新执行步骤S61。
在其中一些实施例中,在步骤S6中,拉线位移计坐标系转换参数的解算还包括以下步骤:
S63、建立标靶点的第二观测方程,表示为式(18);
(18);
式(18)中,、/>、/>为/>坐标系与/>坐标系之间的平移参数;/>、/>为/>坐标系与/>坐标系之间的旋转参数;/>、/>、/>为/>坐标系与/>坐标系之间三个坐标轴方向的尺度;/>,/>为第/>个标靶点的第二测量改正数;
,/>,/>,建立所有标靶点的第二观测误差方程,表示为式(19);
(19);
根据最小二乘法平差计算得到坐标系向/>坐标系转换时的拉线位移计坐标系转换参数/>,表示为式(20);
(20)。
在其中一些实施例中,在步骤S6中,拉线位移计坐标系转换参数的解算还包括以下步骤:
S64、将拉线位移计坐标系转换参数代入到式(19)中,得到/>;计算出/>的标准差/>
S65、根据三倍标准差法,剔除对应的标靶点数据;若剔除的标靶点数据量超过全部标靶点数据量的/>,则重新利用拉线位移计测量/>于/>坐标系下的坐标,然后重新执行步骤S6;若剔除的标靶点数据量未超过全部标靶点数据量的/>,则利用其它未被剔除的标靶点数据,重新执行步骤S62-S63,以解算出修正后的拉线位移计坐标系转换参数/>
基于上述技术方案,本发明实施例中的组合式水下定位设备的标定方法,实现了多种定位装置的组合式设置以及统一标定,显著提高了组合式水下定位设备的标定精度,降低了标定误差,进而能够在工程应用中提高对待测物体的定位精度,提升施工效率和施工精度。
附图说明
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1为组合式水下定位设备的结构示意图(透视化显示);
图2为本发明的组合式水下定位设备的标定方法中各坐标系的示意图;
图3为本发明的组合式水下定位设备的标定方法的基本流程图;
图4为本发明的组合式水下定位设备的标定方法中的全站仪坐标系转换参数解算的步骤图;
图5为本发明的组合式水下定位设备的标定方法中的摄像位移计坐标系转换参数解算的步骤图;
图6为本发明的组合式水下定位设备的标定方法中的拉线位移计坐标系转换参数解算的步骤图。
具体实施方式
下面将结合本发明实施例中的附图,对实施例中的技术方案进行清楚、完整的描述。显然,所描述的实施例仅仅是本发明的一部分实施例,而非全部的实施例。基于本发明的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要理解的是,术语“中心”、“横向”、“纵向”、“上”、“下”、“顶”、“底”、“内”、“外”、“左”、“右”、“前”、“后”、“竖直”、“水平”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
术语“第一”、“第二”等仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”等的特征可以明示或者隐含地包括一个或者更多个该特征。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“安装”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是直接相连,也可以通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
本发明的组合式水下定位设备的标定方法,适用于组合式水下定位设备。参考图1-图2所示,组合式水下定位设备包括仪器箱、安装于仪器箱内的摄像位移计和拉线位移计。仪器箱内具有容纳腔。摄像位移计包括并排设置的两台相机,可以理解的是,摄像位移计的拍摄方向朝向待测物体。拉线位移计包括可转动的拉杆、贯穿拉杆且可伸缩的拉线、测定拉线长度的拉线编码器和测定拉线方向的激光测角装置;可以理解的是,拉杆转动中心位于仪器箱内,拉线与拉杆的转动状态一致。进一步说明,拉线的一端缠绕于一绕线盘上,另一端用于与待测物体连接,当拉线因待测物体的运动而改变方向时,拉杆随拉线同步改变方向;绕线盘用于收放拉线,拉线编码器设于绕线盘处以测量绕线盘的旋转角度,进而计算拉线长度。激光测角装置包括激光器和光斑相机;激光器安装于拉杆上以随拉杆同步转动,激光器用于向外发出激光;光斑相机面向激光器,以采集激光于相片上形成的像素点位置信息,进而计算拉杆也即拉线方向。进一步地,该组合式水下定位设备实现了多种定位装置的组合式集成设置,解决了现有技术中多种定位装置分别安装难度大、占用空间大、敷设线缆多的问题,由此能够提高施工效率、降低施工成本。
参考图1-图3所示,组合式水下定位设备的标定方法包括以下步骤:
S1、建立基于仪器箱的仪器坐标系,将其记作坐标系;建立基于摄像位移计的摄像位移计坐标系,将其记作/>坐标系;建立基于拉线位移计的拉线位移计坐标系,将其记作坐标系;
具体地,将坐标系的原点记作/>,/>位于仪器箱顶面的中心;将/>坐标系的主轴记作/>轴,/>轴朝向摄像位移计的拍摄方向;将/>坐标系的横轴记作/>轴,/>轴朝向右侧90°方向;构成左手坐标系的第三轴方向记作/>轴,即/>坐标系的竖轴,/>轴朝上;
具体地,将坐标系的原点记作/>,/>位于摄像位移计的左相机的摄像中心;将/>坐标系的主轴记作/>轴,/>轴朝向左相机主光轴方向;将/>坐标系的横轴记作/>轴,/>轴朝向摄像位移计的右相机的摄像中心,可以理解的是,左相机的摄像中心和右相机的摄像中心处于同一水平线上;构成左手坐标系的第三轴方向记作/>轴,即/>坐标系的竖轴;
具体地,将坐标系的原点记作/>,/>位于拉线位移计的拉杆转动中心;将/>坐标系的主轴记作/>轴,/>轴朝向拉线位移计正前方;将/>坐标系的横轴记作/>轴,/>轴朝向右侧90°方向;构成左手坐标系的第三轴方向记作/>轴,即/>坐标系的竖轴。
S2、将仪器箱放置于地面上并使其保持水平状态;
延长轴,可采用细线沿/>轴向前延伸,或在仪器箱顶面设置激光指示器,利用激光指示器发出的激光来指示/>轴的方向;在/>轴上间隔设置两个参考点/>、/>,将其在坐标系下的坐标分别记作/>、/>,可以理解的是,/>为/>至/>的直线距离,/>为/>至/>的直线距离;
沿轴方向间隔设置多个标靶并在每一标靶上布设多个标靶点,对所有标靶点编号并记作/>,其中,/>或/>或…或/>;具体地,每个标靶上的多个标靶点相对于标靶中心呈对称设置,且能遍布上下左右范围。
S3、架设全站仪以测量、/>、/>及/>于全站仪坐标系下的坐标,将其分别记作、/>、/>、/>;需要说明的是,全站仪设站时尽量使全站仪零方向与/>轴方向平行;
利用摄像位移计测量于/>坐标系下的坐标,将其记作/>
利用拉线位移计测量于/>坐标系下的坐标,将其记作/>;需要说明的是,/>为拉线编码器的读数,/>和/>分别为激光测角装置中光斑相机拍摄的激光光点于相片上像素点的横向位置和竖向位置。
S4、通过、/>、/>于全站仪坐标系下的坐标和于/>坐标系下的坐标,解算全站仪坐标系向/>坐标系转换时的全站仪坐标系转换参数;然后根据该全站仪坐标系转换参数,将于全站仪坐标系下的坐标/>换算为/>于/>坐标系下的坐标,将其记作
S5、利用步骤S4换算得到的于/>坐标系下的坐标/>和步骤S3中摄像位移计测量得到的/>于/>坐标系下的坐标/>,建立标靶点/>的第一观测方程,以解算/>坐标系向/>坐标系转换时的摄像位移计坐标系转换参数/>
S6、利用步骤S4换算得到的于/>坐标系下的坐标/>和步骤S3中拉线位移计测量得到的/>于/>坐标系下的坐标/>,建立标靶点/>的第二观测方程,以解算/>坐标系向/>坐标系转换时的拉线位移计坐标系转换参数/>
需要说明的是,通过该实施例能够得到摄像位移计坐标系向仪器坐标系的转换参数以及拉线位移计坐标系向仪器坐标系的转换参数;当该组合式水下定位设备应用于工程中时,能够将摄像位移计对待测物体的测量结果和拉线位移计对待测物体的测量结果均换算成同一仪器坐标系下的测量数据,因而,能够将两种定位装置换算后的测量数据进行相互校核,进而提高对待测物体的定位精度。
上述示意性实施例,解决了现有技术中多种定位装置独立安装使用、各自单独标定而导致标定误差较大的问题,实现了摄像位移计和拉线位移计的统一标定,提高了标定精度和标定效率,进而能够在工程应用中提高对待测物体的定位精度。
参考图1-图4所示,在一些实施例中,在步骤S4中,全站仪坐标系转换参数包括水平方向角和高度转换系数;全站仪坐标系转换参数的解算包括以下步骤:
S41、在全站仪坐标系下,将、/>、/>的坐标与/>坐标/>相减以归置零,得到/>、/>、/>归置零后的新坐标,将其分别记作/>、/>,可以理解的是此时/>
S42、计算水平方向角,/>
S43、根据式(1)计算归置零后的新坐标/>转换到/>坐标系下的坐标/>,根据式(2)计算/>归置零后的新坐标/>转换到/>坐标系下的坐标/>;根据式(3)计算第一高度转换系数/>和第二高度转换系数/>
(1);
(2);
(3)。
上述示意性实施例,实现了全站仪坐标系转换参数的解算。
参考图1-图4所示,在一些实施例中,在步骤S4中,
S44、全站仪坐标系转换精度验证:判断、/>是否位于/>范围内,/>为预设测量容差;判断/>是否等于/>、/>是否等于/>;判断/>和/>是否相等;若三项判断结果均为是,则验证通过;若三项判断结果中任一项为否,则验证未通过,说明计算有差错或全站仪测量数据有差错,重新架设全站仪测量/>、/>、/>及/>于全站仪坐标系下的坐标,然后重新执行步骤S4。
上述示意性实施例,实现了全站仪坐标系转换参数的转换精度的验证,确保全站仪坐标系转换参数的准确性。
参考图1-图4所示,在一些实施例中,在步骤S4中,标靶点于/>坐标系下的坐标的解算包括以下步骤:
S45、在全站仪坐标系下,将的坐标/>与/>坐标/>相减以归置零,得到/>归置零后的新坐标,将其记作/>。/>
S46、根据式(4)计算归置零后的新坐标/>转换到/>坐标系下的坐标/>
(4)。
上述示意性实施例,实现了将标靶点于全站仪坐标系下的坐标换算为/>于/>坐标系下的坐标。
参考图1-图5所示,在一些实施例中,在步骤S5中,摄像位移计坐标系转换参数的解算包括以下步骤:
S51、建立标靶点的第一观测方程,表示为式(5);
(5);
式(5)的计算方法为九参数转换法,其中,、/>、/>为/>坐标系与/>坐标系之间的平移参数;/>、/>、/>为/>坐标系与/>坐标系之间的旋转参数;/>、/>、/>为/>坐标系与/>坐标系之间三个坐标轴方向的尺度。需要说明的是,/>为第/>个标靶点的第一测量改正数;
,/>,/>,/>,建立所有标靶点的第一观测误差方程,表示为式(6);
(6);
根据最小二乘法平差计算得到坐标系向/>坐标系转换时的摄像位移计坐标系转换参数/>,表示为式(7);
(7)。
上述示意性实施例,实现了坐标系向/>坐标系转换时的摄像位移计坐标系转换参数/>的解算。
参考图1-图5所示,在一些实施例中,在步骤S5中,摄像位移计坐标系转换参数的解算还包括以下步骤:
S52、将摄像位移计坐标系转换参数代入到式(6)中,得到/>;计算出/>的标准差/>,标准差的计算方法为本领域技术人员所知晓,在此不做展开和赘述;
S53、根据三倍标准差法,剔除对应的标靶点数据;若剔除的标靶点数据量超过全部标靶点数据量的/>,则重新利用摄像位移计测量/>于/>坐标系下的坐标,然后重新执行步骤S5;若剔除的标靶点数据量未超过全部标靶点数据量的/>,则利用其它未被剔除的标靶点数据,重新执行步骤S51,以解算出修正后的摄像位移计坐标系转换参数/>
上述示意性实施例,实现了摄像位移计坐标系转换参数的修正计算,确保摄像位移计坐标系转换参数/>的准确性,进而提高标定精度。
参考图1-图4、图6所示,在一些实施例中,在步骤S6中,拉线位移计坐标系转换参数的解算包括以下步骤:
S61、拉线位移计测量数据规正所需参数的计算,具体包括:
S611、在坐标系下,量取拉杆转动中心/>与/>坐标系原点/>之间的坐标差,该坐标差为常数,记作/>
S612、分别在拉线未拉出和最长拉出状态下读取拉线编码器的读数、/>,同时用卷尺量取拉线的拉出距离/>;根据式(8)计算拉线尺度/>
(8);/>
S613、在拉线未拉出时,量取拉杆转动中心至拉线的拉环内侧的距离/>;根据式(9)计算拉线读数常数/>
(9);
S614、将拉线拉出至接近最长距离,调准拉线方向使其与坐标系的/>轴保持平行,读取激光光点于相片上像素点的横向位置和纵向位置的读数,分别记作光线中心位置常数/>和/>
S615、在一标靶上选择同一水平位置上的两个水平标靶点,将该两个水平标靶点的间距记作,将拉杆转动中心/>到该标靶中心的距离记作/>,根据式(10)计算该两个水平标靶点处拉线的夹角/>;利用拉线位移计读取拉线位于该两个水平标靶点时的光斑相机读数,相减得到读数差/>,将光斑相机焦距记作/>、像素点大小记作/>,根据式(11)计算对应该读数差的夹角/>;然后根据式(12)计算横向张角系数/>
在该标靶上选择同一竖直位置上的两个竖直标靶点,将该两个竖直标靶点的间距记作,根据式(13)计算该两个竖直标靶点处拉线的夹角/>;利用拉线位移计读取拉线位于该两个竖直标靶点时的光斑相机读数,相减得到读数差/>,根据式(14)计算对应该读数差的夹角/>;然后根据式(15)计算纵向张角系数/>
(10);
(11);
(12);
(13);
(14);
(15)。
S62、对拉线位移计的测量数据进行规正计算:
根据式(16)和式(17)计算出规正后的拉线位移计的测量数据,将其记作
(16);
(17)。
上述示意性实施例,实现了拉线位移计测量数据的规正计算,便于进行后续标定工作。
参考图1-图4、图6所示,在一些实施例中,在步骤S61中,拉线位移计测量数据规正所需参数的计算还包括:
S616、拉线位移计规正参数验证:在另一位置处的标靶上选择同一水平位置上的两个水平标靶点和同一竖直位置上的两个竖直标靶点,参考步骤S615,重新计算横向张角系数和纵向张角系数/>;判断/>和/>的两次计算结果的差值是否小于本次计算数值的20%;若是,则验证通过,以距离拉线位移计最近的标靶位置处计算的/>和/>的数值为准;若否,则验证未通过,重新执行步骤S61。
上述示意性实施例,实现了拉线位移计规正参数的验证,确保拉线位移计规正参数的准确性,进而确保拉线位移计测量数据的规正计算的准确性。
参考图1-图4、图6所示,在一些实施例中,在步骤S6中,拉线位移计坐标系转换参数的解算还包括以下步骤:
S63、建立标靶点的第二观测方程,表示为式(18);
(18);
式(18)的计算方法为九参数转换法,其中,、/>、/>为/>坐标系与/>坐标系之间的平移参数;/>、/>、/>为/>坐标系与/>坐标系之间的旋转参数;/>、/>、/>为/>坐标系与/>坐标系之间三个坐标轴方向的尺度。需要说明的是,,/>为第/>个标靶点的第二测量改正数;
,/>,/>,/>,建立所有标靶点的第二观测误差方程,表示为式(19);
(19);
根据最小二乘法平差计算得到坐标系向/>坐标系转换时的拉线位移计坐标系转换参数/>,表示为式(20);
(20)。
上述示意性实施例,实现了坐标系向/>坐标系转换时的拉线位移计坐标系转换参数/>的解算。
参考图1-图4、图6所示,在一些实施例中,在步骤S6中,拉线位移计坐标系转换参数的解算还包括以下步骤:
S64、将拉线位移计坐标系转换参数代入到式(19)中,得到/>;计算出/>的标准差/>,标准差的计算方法为本领域技术人员所知晓,在此不做展开和赘述;
S65、根据三倍标准差法,剔除对应的标靶点数据;若剔除的标靶点数据量超过全部标靶点数据量的/>,则重新利用拉线位移计测量/>于/>坐标系下的坐标,然后重新执行步骤S6;若剔除的标靶点数据量未超过全部标靶点数据量的/>,则利用其它未被剔除的标靶点数据,重新执行步骤S62-S63,以解算出修正后的拉线位移计坐标系转换参数/>
上述示意性实施例,实现了拉线位移计坐标系转换参数的修正计算,确保拉线位移计坐标系转换参数/>的准确性,进而提高标定精度。
综上所述,本发明的组合式水下定位设备的标定方法,实现了多种定位装置的组合式设置以及统一标定,显著提高了组合式水下定位设备的标定精度,降低了标定误差,进而能够在工程应用中提高对待测物体的定位精度,提升施工效率和施工精度。
最后应当说明的是:本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
以上实施例仅用以说明本发明的技术方案而非对其限制;尽管参照较佳实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者对部分技术特征进行等同替换;而不脱离本发明技术方案的精神,其均应涵盖在本发明请求保护的技术方案范围当中。

Claims (9)

1.一种组合式水下定位设备的标定方法,其特征在于,
所述组合式水下定位设备包括仪器箱、安装于所述仪器箱内的摄像位移计和拉线位移计;所述摄像位移计包括并排设置的两台相机;所述拉线位移计包括可转动的拉杆、贯穿所述拉杆且可伸缩的拉线、测定所述拉线长度的拉线编码器和测定所述拉线方向的激光测角装置;
所述组合式水下定位设备的标定方法包括以下步骤:
S1、分别建立基于所述仪器箱的坐标系、基于所述摄像位移计的/>坐标系、基于所述拉线位移计的/>坐标系;将/>坐标系的原点记作/>,/>位于仪器箱顶面的中心,将/>坐标系的主轴记作/>轴,/>轴朝向摄像位移计的拍摄方向;
S2、将所述仪器箱放置于地面上并使其保持水平状态;在轴上间隔设置两个参考点、/>,将其在/>坐标系下的坐标分别记作/>、/>;沿/>轴方向间隔设置多个标靶并在每一标靶上布设多个标靶点,对所有所述标靶点编号并记作/>,其中,/>或/>或…或/>
S3、架设全站仪以测量、/>、/>及/>于全站仪坐标系下的坐标,将其分别记作、/>、/>、/>;利用摄像位移计测量/>于/>坐标系下的坐标,将其记作/>;利用拉线位移计测量/>于/>坐标系下的坐标,将其记作/>
S4、通过、/>、/>于全站仪坐标系下的坐标和于/>坐标系下的坐标,解算全站仪坐标系向/>坐标系转换时的全站仪坐标系转换参数;根据所述全站仪坐标系转换参数,将/>于全站仪坐标系下的坐标/>换算为/>于/>坐标系下的坐标,将其记作;所述全站仪坐标系转换参数包括水平方向角和高度转换系数,全站仪坐标系转换参数的解算包括以下步骤:
S41、在全站仪坐标系下,将、/>、/>的坐标与/>坐标/>相减以归置零,得到/>、/>、/>归置零后的新坐标,将其分别记作/>、/>,此时/>
S42、计算水平方向角,/>
S43、根据式(1)计算转换到/>坐标系下的坐标/>,根据式(2)计算/>转换到/>坐标系下的坐标/>;根据式(3)计算第一高度转换系数/>和第二高度转换系数/>
(1);
(2);
(3);
S5、利用和/>,建立标靶点/>的第一观测方程,以解算/>坐标系向/>坐标系转换时的摄像位移计坐标系转换参数/>
S6、利用和/>,建立标靶点/>的第二观测方程,以解算/>坐标系向/>坐标系转换时的拉线位移计坐标系转换参数/>
2.根据权利要求1所述的组合式水下定位设备的标定方法,其特征在于,在步骤S4中,所述全站仪坐标系转换参数的解算包括以下步骤:
S44、全站仪坐标系转换精度验证:判断、/>是否位于/>范围内,/>为预设测量容差;判断/>是否等于/>、/>是否等于/>;判断/>和/>是否相等;若三项判断结果均为是,则验证通过;若三项判断结果中任一项为否,则验证未通过,重新架设全站仪测量/>、/>、/>及/>于全站仪坐标系下的坐标,然后重新执行步骤S4。
3.根据权利要求2所述的组合式水下定位设备的标定方法,其特征在于,在步骤S4中,于/>坐标系下的坐标/>的解算包括以下步骤:
S45、在全站仪坐标系下,将的坐标/>与/>坐标/>相减以归置零,得到/>归置零后的新坐标,将其记作/>
S46、根据式(4)计算转换到/>坐标系下的坐标/>
(4)。
4.根据权利要求3所述的组合式水下定位设备的标定方法,其特征在于,
在步骤S1中,将坐标系的原点记作/>,/>位于摄像位移计的一相机的摄像中心,将/>坐标系的主轴记作/>轴,/>轴朝向该相机主光轴方向,将/>坐标系的横轴记作/>轴,/>轴朝向摄像位移计的另一相机的摄像中心;
在步骤S5中,摄像位移计坐标系转换参数的解算包括以下步骤:
S51、建立标靶点的第一观测方程,表示为式(5);
(5);
式(5)中,、/>、/>为/>坐标系与/>坐标系之间的平移参数;/>、/>、/>为/>坐标系与/>坐标系之间的旋转参数;/>、/>、/>为/>坐标系与/>坐标系之间三个坐标轴方向的尺度;/>,/>为第/>个标靶点的第一测量改正数;
,/>,/>,建立所有标靶点的第一观测误差方程,表示为式(6);
(6);
根据最小二乘法平差计算得到坐标系向/>坐标系转换时的摄像位移计坐标系转换参数/>,表示为式(7);
(7)。
5.根据权利要求4所述的组合式水下定位设备的标定方法,其特征在于,在步骤S5中,摄像位移计坐标系转换参数的解算还包括以下步骤:
S52、将摄像位移计坐标系转换参数代入到式(6)中,得到/>;计算出/>的标准差/>
S53、根据三倍标准差法,剔除对应的标靶点数据;若剔除的标靶点数据量超过全部标靶点数据量的/>,则重新利用摄像位移计测量/>于/>坐标系下的坐标,然后重新执行步骤S5;若剔除的标靶点数据量未超过全部标靶点数据量的/>,则利用其它未被剔除的标靶点数据,重新执行步骤S51,以解算出修正后的摄像位移计坐标系转换参数/>
6.根据权利要求3所述的组合式水下定位设备的标定方法,其特征在于,
在步骤S1中,将坐标系的原点记作/>,/>位于拉线位移计的拉杆转动中心,将/>坐标系的主轴记作/>轴,/>轴朝向拉线位移计正前方;
在步骤S3中,在拉线位移计的测量数据中,/>为拉线编码器的读数,和/>分别为激光测角装置中光斑相机拍摄的激光光点于相片上像素点的横向位置和竖向位置;
在步骤S6中,拉线位移计坐标系转换参数的解算包括以下步骤:
S61、拉线位移计测量数据规正所需参数的计算,具体包括:
S611、在坐标系下,量取拉杆转动中心/>与/>坐标系原点/>之间的坐标差,该坐标差为常数,记作/>
S612、分别在拉线未拉出和最长拉出状态下读取拉线编码器的读数、/>,同时用卷尺量取拉线的拉出距离/>;根据式(8)计算拉线尺度/>
(8);
S613、在拉线未拉出时,量取拉杆转动中心至拉线的拉环内侧的距离/>;根据式(9)计算拉线读数常数/>
(9);
S614、将拉线拉出至接近最长距离,调准拉线方向使其与坐标系的/>轴保持平行,读取激光光点于相片上像素点的横向位置和纵向位置的读数,分别记作光线中心位置常数和/>
S615、在一标靶上选择同一水平位置上的两个水平标靶点,将该两个水平标靶点的间距记作,将拉杆转动中心/>到该标靶中心的距离记作/>,根据式(10)计算该两个水平标靶点处拉线的夹角/>;利用拉线位移计读取拉线位于该两个水平标靶点时的光斑相机读数,相减得到读数差/>,将光斑相机焦距记作/>、像素点大小记作/>,根据式(11)计算对应该读数差的夹角/>;然后根据式(12)计算横向张角系数/>
在该标靶上选择同一竖直位置上的两个竖直标靶点,将该两个竖直标靶点的间距记作,根据式(13)计算该两个竖直标靶点处拉线的夹角/>;利用拉线位移计读取拉线位于该两个竖直标靶点时的光斑相机读数,相减得到读数差/>,根据式(14)计算对应该读数差的夹角/>;然后根据式(15)计算纵向张角系数/>
(10);
(11);
(12);
(13);
(14);
(15);
S62、对拉线位移计的测量数据进行规正计算:
根据式(16)和式(17)计算出规正后的拉线位移计的测量数据,将其记作
(16);
(17)。
7.根据权利要求6所述的组合式水下定位设备的标定方法,其特征在于,在步骤S61中,拉线位移计测量数据规正所需参数的计算还包括:
S616、拉线位移计规正参数验证:在另一位置处的标靶上选择同一水平位置上的两个水平标靶点和同一竖直位置上的两个竖直标靶点,参考步骤S615,重新计算和/>;判断和/>的两次计算结果的差值是否小于本次计算数值的20%,若是,则验证通过,若否,则验证未通过,重新执行步骤S61。
8.根据权利要求7所述的组合式水下定位设备的标定方法,其特征在于,在步骤S6中,拉线位移计坐标系转换参数的解算还包括以下步骤:
S63、建立标靶点的第二观测方程,表示为式(18);
(18);
式(18)中,、/>、/>为/>坐标系与/>坐标系之间的平移参数;/>、/>、/>为/>坐标系与/>坐标系之间的旋转参数;/>、/>、/>为/>坐标系与/>坐标系之间三个坐标轴方向的尺度;/>,/>为第/>个标靶点的第二测量改正数;
,/>,/>,建立所有标靶点的第二观测误差方程,表示为式(19);
(19);
根据最小二乘法平差计算得到坐标系向/>坐标系转换时的拉线位移计坐标系转换参数/>,表示为式(20);
(20)。
9.根据权利要求8所述的组合式水下定位设备的标定方法,其特征在于,在步骤S6中,拉线位移计坐标系转换参数的解算还包括以下步骤:
S64、将拉线位移计坐标系转换参数代入到式(19)中,得到/>;计算出/>的标准差
S65、根据三倍标准差法,剔除对应的标靶点数据;若剔除的标靶点数据量超过全部标靶点数据量的/>,则重新利用拉线位移计测量/>于/>坐标系下的坐标,然后重新执行步骤S6;若剔除的标靶点数据量未超过全部标靶点数据量的/>,则利用其它未被剔除的标靶点数据,重新执行步骤S62-S63,以解算出修正后的拉线位移计坐标系转换参数/>
CN202310894786.5A 2023-07-20 2023-07-20 一种组合式水下定位设备的标定方法 Active CN116625240B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310894786.5A CN116625240B (zh) 2023-07-20 2023-07-20 一种组合式水下定位设备的标定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310894786.5A CN116625240B (zh) 2023-07-20 2023-07-20 一种组合式水下定位设备的标定方法

Publications (2)

Publication Number Publication Date
CN116625240A CN116625240A (zh) 2023-08-22
CN116625240B true CN116625240B (zh) 2023-09-19

Family

ID=87636845

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310894786.5A Active CN116625240B (zh) 2023-07-20 2023-07-20 一种组合式水下定位设备的标定方法

Country Status (1)

Country Link
CN (1) CN116625240B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115371639B (zh) * 2022-08-11 2023-04-18 深圳大学 一种水下摄影测量沉管管节对接测量方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204676554U (zh) * 2015-04-17 2015-09-30 武汉理工大学 一种沉管水下施工定位系统
CN110220505A (zh) * 2019-06-11 2019-09-10 江苏海事职业技术学院 适用于多种传感器的控制测量标靶及空间坐标配准方法
CN111189436A (zh) * 2020-01-08 2020-05-22 中国矿业大学(北京) 基于激光靶向跟踪的掘进机位姿测量系统
CN111829472A (zh) * 2019-04-17 2020-10-27 初速度(苏州)科技有限公司 利用全站仪测定传感器间相对位置的方法及装置
CN113295149A (zh) * 2021-05-17 2021-08-24 中铁第四勘察设计院集团有限公司 一种基于联合观测量的cpⅲ坐标计算方法、装置
WO2022016356A1 (zh) * 2020-07-21 2022-01-27 中国科学院长春光学精密机械与物理研究所 一种测绘相机高精度内、外方位元素的标定方法
CN114488093A (zh) * 2021-12-09 2022-05-13 苏州挚途科技有限公司 激光雷达与车辆之间的外参标定方法、装置和电子设备
CN216815272U (zh) * 2022-03-15 2022-06-24 中交第一航务工程局有限公司 沉管隧道最终接头水下安装定位系统
CN116051659A (zh) * 2023-03-30 2023-05-02 山东科技大学 一种线阵相机与2d激光扫描仪联合标定方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN204676554U (zh) * 2015-04-17 2015-09-30 武汉理工大学 一种沉管水下施工定位系统
CN111829472A (zh) * 2019-04-17 2020-10-27 初速度(苏州)科技有限公司 利用全站仪测定传感器间相对位置的方法及装置
CN110220505A (zh) * 2019-06-11 2019-09-10 江苏海事职业技术学院 适用于多种传感器的控制测量标靶及空间坐标配准方法
CN111189436A (zh) * 2020-01-08 2020-05-22 中国矿业大学(北京) 基于激光靶向跟踪的掘进机位姿测量系统
WO2022016356A1 (zh) * 2020-07-21 2022-01-27 中国科学院长春光学精密机械与物理研究所 一种测绘相机高精度内、外方位元素的标定方法
CN113295149A (zh) * 2021-05-17 2021-08-24 中铁第四勘察设计院集团有限公司 一种基于联合观测量的cpⅲ坐标计算方法、装置
CN114488093A (zh) * 2021-12-09 2022-05-13 苏州挚途科技有限公司 激光雷达与车辆之间的外参标定方法、装置和电子设备
CN216815272U (zh) * 2022-03-15 2022-06-24 中交第一航务工程局有限公司 沉管隧道最终接头水下安装定位系统
CN116051659A (zh) * 2023-03-30 2023-05-02 山东科技大学 一种线阵相机与2d激光扫描仪联合标定方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于激光标靶和捷联惯导系统组合位姿测量方法;郭庆尧;林嘉睿;任永杰;高扬;;激光与光电子学进展(01);全文 *
多站激光扫描点云配准算法研究;刘志超;顾广杰;;海洋测绘(05);全文 *
适合近岸管节水下施工的测量塔全站仪定位方法;吴瑞大;方长远;;中国港湾建设(04);全文 *

Also Published As

Publication number Publication date
CN116625240A (zh) 2023-08-22

Similar Documents

Publication Publication Date Title
CN116625240B (zh) 一种组合式水下定位设备的标定方法
CN116592768B (zh) 沉管隧道最终接头水下安装定位系统及定位方法
CN200986436Y (zh) 检测可见激光束位置偏差的装置
Fryer et al. On the calibration of underwater cameras
CN115810041B (zh) 一种基于多摄像头视觉的线缆尺寸检测方法及系统
CN111366902A (zh) 卫星热变形试验相对指向变化测量系统及方法
CN105865349A (zh) 一种大型建筑物位移监测方法
CN115205466A (zh) 一种基于结构光的输电通道三维重建方法和系统
CN101655378A (zh) 一种数字摄影测量相机的内方位参数测量方法
CN111156899A (zh) 一种视觉测量装置、方法及机器人
JP4472618B2 (ja) 線条の弛度測定方法および弛度測定装置
CN114813061A (zh) 一种近眼成像设备的光学参数检测方法及系统
CN113063399A (zh) 一种高程测量方法及系统
CN111336937B (zh) 一种提高圆柱体直径测量精度的2d设备及其测量方法
CN112802115A (zh) 一种多焦面拼接大视场离轴相机的几何标定方法及装置
CN216645327U (zh) 基于机器视觉的整套索股高差测量设备
CN113465480A (zh) 检测场地平整度的方法和电子设备
CN115290008A (zh) 一种影像测量平行光管的角度标定算法
CN109029404B (zh) 一种手机摄影标定地面点的方法
CN112507871B (zh) 一种巡检机器人及其检测方法
CN111197486A (zh) 一种基于机器视觉的顶管机导向方法
JP5243805B2 (ja) 水中距離測定方法および水中距離測定システム
WO2024017408A1 (zh) 沉管隧道最终接头水下安装定位系统及定位方法
CN114964734B (zh) 一种光纤光谱天文望远镜的光纤极坐标位置测量方法
CN114353765B (zh) 一种倾斜巷道腰线施放方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant