JP7841786B2 - 視覚とimuとの融合に基づく海洋機器の水中損傷の3次元再構成方法 - Google Patents

視覚とimuとの融合に基づく海洋機器の水中損傷の3次元再構成方法

Info

Publication number
JP7841786B2
JP7841786B2 JP2025512740A JP2025512740A JP7841786B2 JP 7841786 B2 JP7841786 B2 JP 7841786B2 JP 2025512740 A JP2025512740 A JP 2025512740A JP 2025512740 A JP2025512740 A JP 2025512740A JP 7841786 B2 JP7841786 B2 JP 7841786B2
Authority
JP
Japan
Prior art keywords
imu
underwater
stereo camera
platform
laser sensor
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
JP2025512740A
Other languages
English (en)
Other versions
JP2025534947A (ja
Inventor
振民 王
鵬 遅
海鵬 廖
済語 田
▲チン▼ 張
Original Assignee
華南理工大学
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 華南理工大学 filed Critical 華南理工大学
Publication of JP2025534947A publication Critical patent/JP2025534947A/ja
Application granted granted Critical
Publication of JP7841786B2 publication Critical patent/JP7841786B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/80Analysis of captured images to determine intrinsic or extrinsic camera parameters, i.e. camera calibration
    • G06T7/85Stereo camera calibration
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • G01C21/1652Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments with ranging devices, e.g. LIDAR or RADAR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • G01C21/1656Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments with passive imaging devices, e.g. cameras
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/521Depth or shape recovery from laser ranging, e.g. using interferometry; from the projection of structured light
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/50Depth or shape recovery
    • G06T7/55Depth or shape recovery from multiple images
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10028Range image; Depth image; 3D point clouds

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Theoretical Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Optics & Photonics (AREA)
  • Length Measuring Devices By Optical Means (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Traffic Control Systems (AREA)

Description

本発明は、水中3次元再構成の技術分野に関し、具体的には、視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法に関する。
船舶、海上石油・ガスプラットフォーム、海上風力発電設備などの海洋機器は、長期間にわたり巨大な波、湿潤な環境、海水の腐食、衝突などの不利な要因の影響を受け、構造的な損傷が発生しやすい。従来の方法では、修理のために帰港したり、手作業で水中検査したりする必要があり、多大な時間や経済的コストがかかるとともに、多くの安全上のリスクも生じる。水中移動プラットフォームを使用して損傷の測位を行うことにより、上記問題をうまく解決することができる。水中移動プラットフォームの自律測位及び3次元再構成技術により、明確かつ正確な水中損傷モデルを確立して、自律修復システムと合わせて上記の修理作業を完了することができる。
現在、一般的に使用されている水中測位システムには、水中音響測位システム及び水中SLAM手法がある。水中音響測位システムには、超短基線測位、短基線測位、長基線測位が含まれており、機器が高価で、設置が困難である。一般的なSLAM手法には、ソナーによる手法とカメラによる手法とが含まれており、ソナー機器は、高価であり、音響手法なので解像度が低いため、深海測位に適しているが、カメラによる手法は、水中の光屈折を克服する必要があり、光線に対して敏感であるため、特徴点が不明瞭である場合に、測位失敗という問題が発生しやすい。
従来のカメラを用いた3次元再構成システムには、カメラ歪み補正と3次元再構成の2つの手順がある。カメラ歪み補正には、単一視点モデルに基づく方法と、キャリブレーションボード又は補助ハードウェアに基づく方法とが含まれる。単一視点モデルに基づく方法では、透視モデルのみが考慮され、水中の屈折モデルが考慮されていないため、精度が低くなる。キャリブレーションボード及び補助ハードウェアに基づく方法では、水中の屈折モデルが考慮されているため、精度が高くなる。3次元再構成は、主にカメラパラメータに基づいて3次元点群を直接又は間接的に取得することにより、測位データによって3次元点群の重ね合わせを行うが、上述したようにカメラのみを用いた水中測位の精度が低いため、3次元再構成の精度も低下する。
従来技術の欠点及び不足を克服するために、本発明の目的は、視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法を提供することにある。この方法は、損傷領域の高精度な3次元再構成結果を提供し、他の機器による自律修復を支援し、海洋機器の作業効率を向上させることができる。
上記目的を達成するために、本発明は、下記技術的手段により実現される。視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法は、水中損傷の3次元再構成システムによって実現され、前記水中損傷の3次元再構成システムは、水中移動プラットフォームと、ホストコンピュータとを備え、水中移動プラットフォームは、水中移動プラットフォーム本体と、水中移動プラットフォーム本体に搭載されたステレオカメラ、プラットフォーム用IMU、レーザセンサ、レーザ駆動システム、通信システム及び排水システムとを備え、レーザ駆動システムの各軸に駆動用IMUがそれぞれ設けられており、通信システムが水中移動プラットフォームとホストコンピュータとの間の通信に用いられ、
海洋機器の水中損傷の3次元再構成方法は、
S1、ステレオカメラ及びプラットフォーム用IMUをそれぞれ水中移動プラットフォームに固定し、水中環境において、ステレオカメラの左右カメラの内部パラメータキャリブレーション、及びステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションを行うステップと、
S2、レーザセンサを駆動システムの先端に固定し、駆動システムの各軸上の駆動用IMU座標系とレーザセンサ座標系の外部パラメータ行列をキャリブレーションするステップと、
S3、ステレオカメラ画像データ及びプラットフォーム用IMUデータを収集し、プラットフォーム用IMUデータ中の加速度データ及び角速度データを積分して、プラットフォーム用IMU座標系における位置及び姿勢観測を取得し、ステレオカメラ画像データに対応するフレームに対して、ステレオカメラ座標系における位置姿勢観測、損傷検出及び3次元点群生成をそれぞれ行い、位置姿勢観測結果を融合し、連続した3次元点群を重ね合わせて、水中3次元再構成点群を取得し、損傷検出結果の検証に使用し、位置姿勢観測及び損傷検出結果に基づいて、水中移動プラットフォームの最適な経路を計画し、3次元再構成点群に基づいて局所的な障害物回避を行い、水中移動プラットフォームを損傷領域の近傍に移動するように制御するステップと、
S4、排水システムの軌跡を計画し、排水システムにより損傷領域の排水を行い、S2で得られた外部パラメータ行列に基づき、駆動用IMUデータを用いてレーザセンサのレーザ位置を決定することにより、損傷領域の精細なレーザセンサデータの3次元再構成を実現するステップと、を含むことを特徴とする。
好ましくは、前記S1は、
S11、ステレオカメラを水中移動プラットフォーム本体の前端に固定し、視覚方向 を斜め下向きに10°~30°の間にし、プラットフォーム用IMUを水中移動プラットフォームの質量中心位置に相当する水中移動プラットフォーム本体の中部に固定するステップと、
S12、キャリブレーションボード及び水中移動プラットフォームを同時に水中に配置し、キャリブレーションボードがステレオカメラの左右カメラの視野に同時に現れるようにし、
キャリブレーションボードがステレオカメラの左右カメラの視野の各位置に分布するように、水中移動プラットフォームを移動させ、複数組のステレオカメラ画像データを記録し、通信システムにより複数組のステレオカメラ画像データをホストコンピュータに送信し、ホストコンピュータにより、関連するキャリブレーション計算、ステレオカメラの左右カメラの内部パラメータキャリブレーション、及びステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションを行うステップと、を含む。
好ましくは、前記S12では、
ステレオカメラの左右カメラの内部パラメータキャリブレーションは、
(式中、lは左カメラを表し、rは右カメラを表し、K,Kはそれぞれ左右カメラの内部パラメータ行列を表し、fxl,fyl,fxr,fyrはそれぞれ左右カメラのx軸方向及びy軸方向における焦点距離の長さをピクセル単位で表し、(u0l,vol),(u0r,v0r)はそれぞれ左右カメラの画像平面座標系における主点の実際のピクセル座標を表す)であり、
ステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションは、
プラットフォーム用IMU座標系を世界座標系とすると、ステレオカメラの左右カメラの画像点からプラットフォーム用IMU座標系への変換関係が、
(式中、
はそれぞれ左右カメラ座標系における2次元座標であり、
はプラットフォーム用IMU座標系における3次元座標であり、Rlr,Rriはそれぞれ右カメラから左カメラへ、左カメラからプラットフォーム用IMU座標系への3*3の回転行列であり、Tlr,Triはそれぞれ右カメラから左カメラへ、左カメラからプラットフォーム用IMU座標系への1*3の並進ベクトルである)である。
好ましくは、前記S2では、駆動システムの各軸上の駆動用IMU座標系とレーザセンサ座標系とを位置合わせすることは、
ステレオカメラ、排水システム及びレーザ駆動システムの位置関係に基づいて、水中移動プラットフォーム質量中心座標系から排水システム座標系への変換関係、排水システム座標系とレーザ駆動システム座標系との変換関係を取得し、
レーザセンサのレーザポイントをパラメータが既知のキャリブレーションボードを移動するように制御し、通信システムによりレーザセンサと駆動用IMUを接続してデータを取得し、ホストコンピュータに送信し、ホストコンピュータによりキャリブレーション計算を行い、レーザ駆動システム座標系とレーザセンサ座標系との変換関係を取得し、
レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心の4つの座標系を位置合わせする。
好ましくは、レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心の4つの座標系を位置合わせする方法は、
レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心のいずれか2つの座標系の回転行列及び並進ベクトルを含む外部パラメータ行列をキャリブレーションする。
(式中、A及びBはそれぞれ2つの座標系を表し、Xは4*4の外部パラメータ行列を表し、Rは3*3の回転行列を表し、Tは1*3の並進ベクトルを表す)
好ましくは、前記S3では、プラットフォーム用IMU座標系における姿勢観測は、
プラットフォーム用IMUデータを時刻kから時刻k+1まで積分して得られた速度V、並進ベクトルT、及び回転行列Rがそれぞれ
(式中、V,Vk+1はそれぞれ時刻k及び時刻k+1における速度であり、aは加速度であり、Δtは時間間隔であり、T,Tk+1はそれぞれ時刻k及び時刻k+1における並進ベクトルであり、R,Rk+1はそれぞれ時刻k及び時刻k+1における回転行列であり、ωは角速度であり、
はクロネッカー積である)で表される。
好ましくは、前記S3では、ステレオカメラ座標系における姿勢観測は、
ステレオカメラ画像データに対して特徴点抽出を行い、特徴点を中心として円領域を構成し、
θ=arctan(m01/m10
[式中、Cは円領域の質量中心を表し、θは特徴点の方向ベクトルを表し、mpqは円領域のモーメントを表し、
(式中、Rは円領域の半径を表し、x,yはx軸座標、y軸座標を表し、I(x,y)はグレースケール方程式を表す)で定義される]
連続した複数フレームのステレオカメラ画像データの特徴点を抽出してマッチングすることにより、マッチングされた画素点を用いてPnP問題を設定し、ステレオカメラの回転行列R及び並進ベクトルTを得る。
好ましくは、前記S3では、ステレオカメラ座標系における3次元点群生成は、
同一フレームのステレオカメラの左右カメラ画像に対して上記特徴点を抽出してマッチングし、グレースケール二乗和誤差アルゴリズムに基づいて視差計算を行い、
(式中、x,y,dはそれぞれx軸座標、y軸座標、視差であり、i,jはそれぞれx軸方向及びy軸方向の変化値であり、m,nはそれぞれx軸方向及びy軸方向の最大値であり、I(x,y),I(x,y)はグレースケール方程式を表す)
視差及び元の座標により3次元点群データを生成し、3次元座標が、
[式中、x,xはそれぞれ左右カメラに対応する横座標値であり、y,yはそれぞれ左右カメラの縦座標値であり、f,fはそれぞれ左右カメラの内部パラメータの対応する焦点距離であり、X,Y,Zはそれぞれ3次元座標であり、Dは深度値であり、下記式:
D=Bf/d
(式中、Bは基線長であり、fはカメラの焦点距離であり、dは左右画像の視差である)で算出される]で表される。
好ましくは、前記S4は、損傷領域の位置に基づいて、排水システムが損傷領域を被覆して排水して乾燥空間を形成するように、水中移動プラットフォームの移動軌跡を計画し、レーザ駆動システムを使用して、レーザセンサを乾燥空間で3次元走査を行うように制御し、レーザセンサデータと駆動用IMUデータとを通信システムを介してホストコンピュータに送信し、ホストコンピュータは、S2で得られた外部パラメータ行列に基づいて、駆動用IMUデータを用いてレーザ駆動システムの姿勢を取得し、変換してレーザセンサの位置を取得し、レーザセンサの位置及び点群データに基づいて、精細なレーザセンサデータの3次元再構成を得、3次元再構成結果に基づいて損傷位置を検出する。
好ましくは、前記S4では、レーザセンサデータの3次元再構成は、
レーザセンサは、固定周波数でレーザパルスを発射し、受信器により戻ってきた反射光を受光して距離を判断するとともに、反射強度に基づいて対象材質を大まかに区別し、距離測定式が、
L=tc/2
(式中、Lは目標距離であり、tは戻り時間であり、cは光速である)であり、
駆動用IMUを使用してレーザセンサの位置姿勢予測を行い、回転行列R及び並進ベクトルTに基づいてレーザセンサの3次元再構成結果を得る。
本発明は、従来技術に比べて、以下の利点及び有益な効果を有する。
1、本発明は、海洋機器の水中損傷の検出及び3次元再構成を自律的に行うことができるため、人件費や経済的コストを削減するとともに、安全性を向上させる。
2、本発明は、視覚とIMUとの融合により、測位精度及び水中3次元再構成精度を向上させ、画像及び点群の融合検証に基づく損傷検出方法により、水中海洋機器の損傷位置をより正確に測位することができる。
3、本発明は、損傷領域の近傍の水を正確に排水することができるため、高精度なレーザによる3次元再構成及び損傷識別を実現し、他の自律修復機器にも利便性を提供する。
本発明に係る視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法のフローチャートである。 本発明に係る水中損傷の3次元再構成システムの構造概略図である。 本発明に係る視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法の通信を示す図である。 本発明に係る視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法における座標系変換を示す図である。
以下、図面を参照して具体的な実施形態により本発明を詳細に説明する。
(実施例)
本実施例に係る視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法は、具体的なプロセスが図1に示されており、水中損傷の3次元再構成システムによって実現される。
水中損傷の3次元再構成システムは、水中移動プラットフォームと、ホストコンピュータとを備え、図2に示すように、水中移動プラットフォームは、水中移動プラットフォーム本体1と、水中移動プラットフォーム本体1に搭載されたステレオカメラ4、プラットフォーム用IMU6、レーザセンサ5、レーザ駆動システム3、通信システム及び排水システム2とを備える。
レーザ駆動システム3の各軸に駆動用IMUがそれぞれ設けられており、通信システムが水中移動プラットフォームとホストコンピュータとの間の通信に用いられる。具体的には、通信システムは、水中移動プラットフォーム本体の中部に固定され、ステレオカメラ画像データ及びIMUデータを収集してホストコンピュータに送信するとともに、水中移動プラットフォームを駆動するようにホストコンピュータの関連制御指令を受信するために使用される。
海洋機器の水中損傷の3次元再構成方法は、以下のステップS1~S4を含む。
S1:ステレオカメラ及びプラットフォーム用IMUをそれぞれ水中移動プラットフォームに固定し、水中環境において、ステレオカメラの左右カメラの内部パラメータキャリブレーション、及びステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションを行う。
S1は、
S11、ステレオカメラを水中移動プラットフォーム本体の前端に固定し、視覚方向を斜め下向きに10°~30°の間にし、プラットフォーム用IMUを水中移動プラットフォームの質量中心位置に相当する水中移動プラットフォーム本体の中部に固定するステップと、
S12、キャリブレーションボード及び水中移動プラットフォームを同時に水中に配置し、キャリブレーションボードがステレオカメラの左右カメラの視野に同時に現れるようにし、ステレオカメラの視野にキャリブレーションボードが同時に完全に含まれることを確認しながら、駆動用IMU3軸とのキャリブレーションを確保するように、各方向で可能な限り回転させ、このステップにおけるデータ記録時間は長くする必要がなく、ステレオカメラは15フレーム/秒以上で、駆動用IMUは100フレーム/秒以上であり、
キャリブレーションボードがステレオカメラの左右カメラの視野の各位置に分布するように、水中移動プラットフォームを移動させ、複数組のステレオカメラ画像データを記録し、通信システムにより複数組のステレオカメラ画像データをホストコンピュータに送信し、ホストコンピュータにより、関連するキャリブレーション計算、ステレオカメラの左右カメラの内部パラメータキャリブレーション、及びステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションを行うステップと、を含む。
S12では、ステレオカメラの左右カメラの内部パラメータキャリブレーションは、
(式中、lは左カメラを表し、rは右カメラを表し、K,Kはそれぞれ左右カメラの内部パラメータ行列を表し、fxl,fyl,fxr,fyrはそれぞれ左右カメラのx軸方向及びy軸方向における焦点距離の長さをピクセル単位で表し、(u0l,vol),(u0r,v0r)はそれぞれ左右カメラの画像平面座標系における主点の実際のピクセル座標を表す)であり、
ステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションは、
プラットフォーム用IMU座標系を世界座標系とすると、ステレオカメラの左右カメラの画像点からプラットフォーム用IMU座標系への変換関係が、
(式中、
はそれぞれ左右カメラ座標系における2次元座標であり、
はプラットフォーム用IMU座標系における3次元座標であり、Rlr,Rriはそれぞれ右カメラから左カメラへ、左カメラからプラットフォーム用IMU座標系への3*3の回転行列であり、Tlr,Triはそれぞれ右カメラから左カメラへ、左カメラからプラットフォーム用IMU座標系への1*3の並進ベクトルである)である。
S2:レーザセンサを駆動システムの先端に固定し、駆動システムの各軸上の駆動用IMU座標系とレーザセンサ座標系の外部パラメータ行列をキャリブレーションする。
具体的には、図4に示すように、ステレオカメラ、排水システム及びレーザ駆動システムの位置関係に基づいて、水中移動プラットフォーム質量中心座標系から排水システム座標系への変換関係、排水システム座標系とレーザ駆動システム座標系との変換関係を取得し、
レーザセンサのレーザポイントをパラメータが既知のキャリブレーションボードを移動するように制御し、通信システムによりレーザセンサと駆動用IMUを接続してデータを取得し、ホストコンピュータに送信し、ホストコンピュータによりキャリブレーション計算を行い、レーザ駆動システム座標系とレーザセンサ座標系との変換関係を取得し、
レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心の4つの座標系を位置合わせする。オフラインキャリブレーションが完了すると、図4におけるすべての座標系の変換関係が既知となっている。
レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心の4つの座標系を位置合わせする方法は、
レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心のいずれか2つの座標系の回転行列及び並進ベクトルを含む外部パラメータ行列をキャリブレーションする。
(式中、A及びBはそれぞれ2つの座標系を表し、Xは4*4の外部パラメータ行列を表し、Rは3*3の回転行列を表し、Tは1*3の並進ベクトルを表す)
S3:ステレオカメラ画像データ及びプラットフォーム用IMUデータをホストコンピュータにより収集し、前の内部パラメータキャリブレーション結果及び外部パラメータキャリブレーション結果を読み取り、次に、プラットフォーム用IMUデータ及びステレオカメラ画像データを融合して、左カメラ座標系における測位結果を得、ステレオカメラの検出原理に基づいて、左カメラ座標系において損傷情報を検出して現在フレーム画像の3次元点群を生成し、測位結果及び3次元点群情報を融合してフレーム毎の点群をフィルタリングして重ね合わせ、連続した3次元再構成結果を生成し、3次元再構成の点群に基づいてステレオカメラで検出された損傷位置を検証し、次に、左カメラ座標系における測位結果及び損傷領域の位置に基づいて、水中移動プラットフォームの全体移動経路を計画し、水中移動プラットフォーム質量中心座標系に変換し、通信バスを介して水中移動プラットフォームの通信システムに制御信号を送信し、水中移動プラットフォームが移動する際に、水中移動プラットフォームが損傷領域の近傍に移動するまで、リアルタイムの3次元再構成結果により保存された3次元情報に基づいて局所的な障害物回避を行う。
プラットフォーム用IMU座標系における姿勢観測は、
プラットフォーム用IMUデータを時刻kから時刻k+1まで積分して得られた速度V、並進ベクトルT、及び回転行列Rがそれぞれ
(式中、V,Vk+1はそれぞれ時刻k及び時刻k+1における速度であり、aは加速度であり、Δtは時間間隔であり、T,Tk+1はそれぞれ時刻k及び時刻k+1における並進ベクトルであり、R,Rk+1はそれぞれ時刻k及び時刻k+1における回転行列であり、ωは角速度であり、
はクロネッカー積である)で表される。
ステレオカメラ座標系における姿勢観測は、
ステレオカメラ画像データに対して特徴点抽出を行い、特徴点を中心として円領域を構成し、
θ=arctan(m01/m10
[式中、Cは円領域の質量中心を表し、θは特徴点の方向ベクトルを表し、mpqは円領域のモーメントを表し、
(式中、Rは円領域の半径を表し、x、yはx軸座標、y軸座標を表し、I(x,y)はグレースケール方程式を表す)で定義される]
連続した複数フレームのステレオカメラ画像データの特徴点を抽出してマッチングすることにより、マッチングされた画素点を用いてPnP問題を設定し、ステレオカメラの回転行列R及び並進ベクトルTを得る。
ステレオカメラ座標系における3次元点群生成は、
同一フレームのステレオカメラの左右カメラ画像に対して上記特徴点を抽出してマッチングし、グレースケール二乗和誤差アルゴリズムに基づいて視差計算を行い、
(式中、x,y,dはそれぞれx軸座標、y軸座標、視差であり、i,jはそれぞれx軸方向及びy軸方向の変化値であり、m,nはそれぞれx軸方向及びy軸方向の最大値であり、I(x,y),I(x,y)はグレースケール方程式を表す)
視差及び元の座標により3次元点群データを生成し、3次元座標が、
[式中、x,xはそれぞれ左右カメラに対応する横座標値であり、y,yはそれぞれ左右カメラの縦座標値であり、f,fはそれぞれ左右カメラの内部パラメータの対応する焦点距離であり、X,Y,Zはそれぞれ3次元座標であり、Dは深度値であり、下記式:
D=Bf/d
(式中、Bは基線長であり、fはカメラの焦点距離であり、dは左右画像の視差である)で算出される]で表される。
S4:排水システムの軌跡を計画し、排水システムを排水するように制御し、S2で得られた外部パラメータ行列に基づき、駆動用IMUデータを用いてレーザセンサのレーザ位置を決定することにより、損傷領域の精細なレーザセンサデータの3次元再構成を実現する。
具体的には、損傷領域の位置に基づいて、排水システムが損傷領域を被覆して排水して乾燥空間を形成するように、水中移動プラットフォームの移動軌跡を計画し、レーザ駆動システムを使用して、レーザセンサを乾燥空間で3次元走査を行うように制御し、レーザセンサデータと駆動用IMUデータとを通信システムを介してホストコンピュータに送信し、ホストコンピュータは、S2で得られた外部パラメータ行列に基づいて、駆動用IMUデータを用いてレーザ駆動システムの姿勢を取得し、変換してレーザセンサの位置を取得し、レーザセンサの位置及び点群データに基づいて、精細なレーザセンサデータの3次元再構成を得、3次元再構成結果に基づいて損傷位置を検出する。
レーザセンサデータの3次元再構成は、
レーザセンサは、固定周波数でレーザパルスを発射し、受信器により戻ってきた反射光を受光して距離を判断するとともに、反射強度に基づいて対象材質を大まかに区別し、距離測定式が、
L=tc/2
(式中、Lは目標距離であり、tは戻り時間であり、cは光速である)であり、
駆動用IMUを使用してレーザセンサの位置姿勢予測を行い、回転行列R及び並進ベクトルTに基づいてレーザセンサの3次元再構成結果を得る。このとき、精細な損傷位置の検出を実現し、誤差を0.2mm以内にし、他の自律修復機器に高精度な測位結果を提供することができる。
上記実施例は、本発明の好ましい実施形態であるが、本発明の実施形態は、上記実施例に限定されるものではなく、本発明の精神及び原理を逸脱することなく行われるその他の変更、修正、置換、組み合わせ及び簡略化は、いずれも同等の置換であるべきであり、全て本発明の保護範囲に含まれるものとする。
(付記)
(付記1)
視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法であって、
水中損傷の3次元再構成システムによって実現され、前記水中損傷の3次元再構成システムは、水中移動プラットフォームと、ホストコンピュータとを備え、水中移動プラットフォームは、水中移動プラットフォーム本体と、水中移動プラットフォーム本体に搭載されたステレオカメラ、プラットフォーム用IMU、レーザセンサ、レーザ駆動システム、通信システム及び排水システムとを備え、レーザ駆動システムの各軸に駆動用IMUがそれぞれ設けられており、通信システムが水中移動プラットフォームとホストコンピュータとの間の通信に用いられ、
海洋機器の水中損傷の3次元再構成方法は、
S1、ステレオカメラ及びプラットフォーム用IMUをそれぞれ水中移動プラットフォームに固定し、水中環境において、ステレオカメラの左右カメラの内部パラメータキャリブレーション、及びステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションを行うステップと、
S2、レーザセンサを駆動システムの先端に固定し、駆動システムの各軸上の駆動用IMU座標系とレーザセンサ座標系の外部パラメータ行列をキャリブレーションするステップと、
S3、ステレオカメラ画像データ及びプラットフォーム用IMUデータを収集し、プラットフォーム用IMUデータ中の加速度データ及び角速度データを積分して、プラットフォーム用IMU座標系における位置及び姿勢観測を取得し、ステレオカメラ画像データに対応するフレームに対して、ステレオカメラ座標系における位置姿勢観測、損傷検出及び3次元点群生成をそれぞれ行い、位置姿勢観測結果を融合し、連続した3次元点群を重ね合わせて、水中3次元再構成点群を取得し、損傷検出結果の検証に使用し、位置姿勢観測及び損傷検出結果に基づいて、水中移動プラットフォームの最適な経路を計画し、3次元再構成点群に基づいて局所的な障害物回避を行い、水中移動プラットフォームを損傷領域の近傍に移動するように制御するステップと、
S4、排水システムの軌跡を計画し、排水システムにより損傷領域の排水を行い、S2で得られた外部パラメータ行列に基づき、駆動用IMUデータを用いてレーザセンサのレーザ位置を決定することにより、損傷領域の精細なレーザセンサデータの3次元再構成を実現するステップと、を含む、ことを特徴とする視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
(付記2)
前記S1は、
S11、ステレオカメラを水中移動プラットフォーム本体の前端に固定し、視覚方向を斜め下向きに10°~30°の間にし、プラットフォーム用IMUを水中移動プラットフォームの質量中心位置に相当する水中移動プラットフォーム本体の中部に固定するステップと、
S12、キャリブレーションボード及び水中移動プラットフォームを同時に水中に配置し、キャリブレーションボードがステレオカメラの左右カメラの視野に同時に現れるようにし、
キャリブレーションボードがステレオカメラの左右カメラの視野の各位置に分布するように、水中移動プラットフォームを移動させ、複数組のステレオカメラ画像データを記録し、通信システムにより複数組のステレオカメラ画像データをホストコンピュータに送信し、ホストコンピュータにより、関連するキャリブレーション計算、ステレオカメラの左右カメラの内部パラメータキャリブレーション、及びステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションを行うステップと、を含む、ことを特徴とする付記1に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
(付記3)
前記S12では、
ステレオカメラの左右カメラの内部パラメータキャリブレーションは、
(式中、lは左カメラを表し、rは右カメラを表し、K,Kはそれぞれ左右カメラの内部パラメータ行列を表し、fxl,fyl,fxr,fyrはそれぞれ左右カメラのx軸方向及びy軸方向における焦点距離の長さをピクセル単位で表し、(u0l,vol),(u0r,v0r)はそれぞれ左右カメラの画像平面座標系における主点の実際のピクセル座標を表す)
であり、
ステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションは、
プラットフォーム用IMU座標系を世界座標系とすると、ステレオカメラの左右カメラの画像点からプラットフォーム用IMU座標系への変換関係が、
(式中、
はそれぞれ左右カメラ座標系における2次元座標であり、
はプラットフォーム用IMU座標系における3次元座標であり、Rlr,Rriはそれぞれ右カメラから左カメラへ、左カメラからプラットフォーム用IMU座標系への3*3の回転行列であり、Tlr,Triはそれぞれ右カメラから左カメラへ、左カメラからプラットフォーム用IMU座標系への1*3の並進ベクトルである)である、ことを特徴とする付記2に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
(付記4)
前記S2では、駆動システムの各軸上の駆動用IMU座標系とレーザセンサ座標系とを位置合わせすることは、
ステレオカメラ、排水システム及びレーザ駆動システムの位置関係に基づいて、水中移動プラットフォーム質量中心座標系から排水システム座標系への変換関係、排水システム座標系とレーザ駆動システム座標系との変換関係を取得し、
レーザセンサのレーザポイントをパラメータが既知のキャリブレーションボードを移動するように制御し、通信システムによりレーザセンサと駆動用IMUを接続してデータを取得し、ホストコンピュータに送信し、ホストコンピュータによりキャリブレーション計算を行い、レーザ駆動システム座標系とレーザセンサ座標系との変換関係を取得し、
レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心の4つの座標系を位置合わせする、ことを特徴とする付記1に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
(付記5)
レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心の4つの座標系を位置合わせする方法は、
レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心のいずれか2つの座標系の回転行列及び並進ベクトルを含む外部パラメータ行列をキャリブレーションする、ことを特徴とする付記4に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
(式中、A及びBはそれぞれ2つの座標系を表し、Xは4*4の外部パラメータ行列を表し、Rは3*3の回転行列を表し、Tは1*3の並進ベクトルを表す)
(付記6)
前記S3では、プラットフォーム用IMU座標系における姿勢観測は、
プラットフォーム用IMUデータを時刻kから時刻k+1まで積分して得られた速度V、並進ベクトルT、及び回転行列Rがそれぞれ
(式中、V,Vk+1はそれぞれ時刻k及び時刻k+1における速度であり、aは加速度であり、Δtは時間間隔であり、T,Tk+1はそれぞれ時刻k及び時刻k+1における並進ベクトルであり、R,Rk+1はそれぞれ時刻k及び時刻k+1における回転行列であり、ωは角速度であり、
はクロネッカー積である)で表される、ことを特徴とする付記1に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
(付記7)
前記S3では、ステレオカメラ座標系における姿勢観測は、
ステレオカメラ画像データに対して特徴点抽出を行い、特徴点を中心として円領域を構成し、
θ=arctan(m01/m10
[式中、Cは円領域の質量中心を表し、θは特徴点の方向ベクトルを表し、mpqは円領域のモーメントを表し、
(式中、Rは円領域の半径を表し、x,yはx軸座標、y軸座標を表し、I(x,y)はグレースケール方程式を表す)で定義される]
連続した複数フレームのステレオカメラ画像データの特徴点を抽出してマッチングすることにより、マッチングされた画素点を用いてPnP問題を設定し、ステレオカメラの回転行列R及び並進ベクトルTを得る、ことを特徴とする付記6に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
(付記8)
前記S3では、ステレオカメラ座標系における3次元点群生成は、
同一フレームのステレオカメラの左右カメラ画像に対して上記特徴点を抽出してマッチングし、グレースケール二乗和誤差アルゴリズムに基づいて視差計算を行い、
(式中、x,y,dはそれぞれx軸座標、y軸座標、視差であり、i,jはそれぞれx軸方向及びy軸方向の変化値であり、m,nはそれぞれx軸方向及びy軸方向の最大値であり、I(x,y),I(x,y)はグレースケール方程式を表す)
視差及び元の座標により3次元点群データを生成し、3次元座標が、
[式中、x,xはそれぞれ左右カメラに対応する横座標値であり、y,yはそれぞれ左右カメラの縦座標値であり、f,fはそれぞれ左右カメラの内部パラメータの対応する焦点距離であり、X,Y,Zはそれぞれ3次元座標であり、Dは深度値であり、下記式:
D=Bf/d
(式中、Bは基線長であり、fはカメラの焦点距離であり、dは左右画像の視差である)で算出される]で表される、ことを特徴とする付記7に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
(付記9)
前記S4は、損傷領域の位置に基づいて、排水システムが損傷領域を被覆して排水して乾燥空間を形成するように、水中移動プラットフォームの移動軌跡を計画し、レーザ駆動システムを使用して、レーザセンサを乾燥空間で3次元走査を行うように制御し、レーザセンサデータと駆動用IMUデータとを通信システムを介してホストコンピュータに送信し、ホストコンピュータは、S2で得られた外部パラメータ行列に基づいて、駆動用IMUデータを用いてレーザ駆動システムの姿勢を取得し、変換してレーザセンサの位置を取得し、レーザセンサの位置及び点群データに基づいて、精細なレーザセンサデータの3次元再構成を得、3次元再構成結果に基づいて損傷位置を検出する、ことを特徴とする付記1に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
(付記10)
前記S4では、レーザセンサデータの3次元再構成は、
レーザセンサは、固定周波数でレーザパルスを発射し、受信器により戻ってきた反射光を受光して距離を判断するとともに、反射強度に基づいて対象材質を大まかに区別し、距離測定式が、
L=tc/2
(式中、Lは目標距離であり、tは戻り時間であり、cは光速である)であり、
駆動用IMUを使用してレーザセンサの位置姿勢予測を行い、回転行列R及び並進ベクトルTに基づいてレーザセンサの3次元再構成結果を得る、ことを特徴とする付記9に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。

Claims (10)

  1. 視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法であって、
    水中損傷の3次元再構成システムによって実現され、前記水中損傷の3次元再構成システムは、水中移動プラットフォームと、ホストコンピュータとを備え、水中移動プラットフォームは、水中移動プラットフォーム本体と、水中移動プラットフォーム本体に搭載されたステレオカメラ、プラットフォーム用IMU、レーザセンサ、レーザ駆動システム、通信システム及び排水システムとを備え、レーザ駆動システムの各軸に駆動用IMUがそれぞれ設けられており、通信システムが水中移動プラットフォームとホストコンピュータとの間の通信に用いられ、
    海洋機器の水中損傷の3次元再構成方法は、
    S1、ステレオカメラ及びプラットフォーム用IMUをそれぞれ水中移動プラットフォームに固定し、水中環境において、ステレオカメラの左右カメラの内部パラメータキャリブレーション、及びステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションを行うステップと、
    S2、レーザセンサを駆動システムの先端に固定し、駆動システムの各軸上の駆動用IMU座標系とレーザセンサ座標系の外部パラメータ行列をキャリブレーションするステップと、
    S3、ステレオカメラ画像データ及びプラットフォーム用IMUデータを収集し、プラットフォーム用IMUデータ中の加速度データ及び角速度データを積分して、プラットフォーム用IMU座標系における位置及び姿勢観測を取得し、ステレオカメラ画像データに対応するフレームに対して、ステレオカメラ座標系における位置姿勢観測、損傷検出及び3次元点群生成をそれぞれ行い、位置姿勢観測結果を融合し、連続した3次元点群を重ね合わせて、水中3次元再構成点群を取得し、損傷検出結果の検証に使用し、位置姿勢観測及び損傷検出結果に基づいて、水中移動プラットフォームの最適な経路を計画し、3次元再構成点群に基づいて局所的な障害物回避を行い、水中移動プラットフォームを損傷領域の近傍に移動するように制御するステップと、
    S4、排水システムの軌跡を計画し、排水システムにより損傷領域の排水を行い、S2で得られた外部パラメータ行列に基づき、駆動用IMUデータを用いてレーザセンサのレーザ位置を決定することにより、損傷領域の精細なレーザセンサデータの3次元再構成を実現するステップと、を含む、ことを特徴とする視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
  2. 前記S1は、
    S11、ステレオカメラを水中移動プラットフォーム本体の前端に固定し、視覚方向を斜め下向きに10°~30°の間にし、プラットフォーム用IMUを水中移動プラットフォームの質量中心位置に相当する水中移動プラットフォーム本体の中部に固定するステップと、
    S12、キャリブレーションボード及び水中移動プラットフォームを同時に水中に配置し、キャリブレーションボードがステレオカメラの左右カメラの視野に同時に現れるようにし、
    キャリブレーションボードがステレオカメラの左右カメラの視野の各位置に分布するように、水中移動プラットフォームを移動させ、複数組のステレオカメラ画像データを記録し、通信システムにより複数組のステレオカメラ画像データをホストコンピュータに送信し、ホストコンピュータにより、関連するキャリブレーション計算、ステレオカメラの左右カメラの内部パラメータキャリブレーション、及びステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションを行うステップと、を含む、ことを特徴とする請求項1に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
  3. 前記S12では、
    ステレオカメラの左右カメラの内部パラメータキャリブレーションは、

    (式中、lは左カメラを表し、rは右カメラを表し、K,Kはそれぞれ左右カメラの内部パラメータ行列を表し、fxl,fyl,fxr,fyrはそれぞれ左右カメラのx軸方向及びy軸方向における焦点距離の長さをピクセル単位で表し、(u0l,vol),(u0r,v0r)はそれぞれ左右カメラの画像平面座標系における主点の実際のピクセル座標を表す)
    であり、
    ステレオカメラの左右カメラとプラットフォーム用IMUとの外部パラメータキャリブレーションは、
    プラットフォーム用IMU座標系を世界座標系とすると、ステレオカメラの左右カメラの画像点からプラットフォーム用IMU座標系への変換関係が、


    (式中、

    はそれぞれ左右カメラ座標系における2次元座標であり、

    はプラットフォーム用IMU座標系における3次元座標であり、Rlr,Rriはそれぞれ右カメラから左カメラへ、左カメラからプラットフォーム用IMU座標系への3*3の回転行列であり、Tlr,Triはそれぞれ右カメラから左カメラへ、左カメラからプラットフォーム用IMU座標系への1*3の並進ベクトルである)である、ことを特徴とする請求項2に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
  4. 前記S2では、駆動システムの各軸上の駆動用IMU座標系とレーザセンサ座標系とを位置合わせすることは、
    ステレオカメラ、排水システム及びレーザ駆動システムの位置関係に基づいて、水中移動プラットフォーム質量中心座標系から排水システム座標系への変換関係、排水システム座標系とレーザ駆動システム座標系との変換関係を取得し、
    レーザセンサのレーザポイントをパラメータが既知のキャリブレーションボードを移動するように制御し、通信システムによりレーザセンサと駆動用IMUを接続してデータを取得し、ホストコンピュータに送信し、ホストコンピュータによりキャリブレーション計算を行い、レーザ駆動システム座標系とレーザセンサ座標系との変換関係を取得し、
    レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心の4つの座標系を位置合わせする、ことを特徴とする請求項1に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
  5. レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心の4つの座標系を位置合わせする方法は、
    レーザセンサ、レーザ駆動システム、排水システム、及び水中移動プラットフォーム質量中心のいずれか2つの座標系の回転行列及び並進ベクトルを含む外部パラメータ行列をキャリブレーションする、ことを特徴とする請求項4に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。

    (式中、A及びBはそれぞれ2つの座標系を表し、Xは4*4の外部パラメータ行列を表し、Rは3*3の回転行列を表し、Tは1*3の並進ベクトルを表す)
  6. 前記S3では、プラットフォーム用IMU座標系における姿勢観測は、
    プラットフォーム用IMUデータを時刻kから時刻k+1まで積分して得られた速度V、並進ベクトルT、及び回転行列Rがそれぞれ



    (式中、V,Vk+1はそれぞれ時刻k及び時刻k+1における速度であり、aは加速度であり、Δtは時間間隔であり、T,Tk+1はそれぞれ時刻k及び時刻k+1における並進ベクトルであり、R,Rk+1はそれぞれ時刻k及び時刻k+1における回転行列であり、ωは角速度であり、

    はクロネッカー積である)で表される、ことを特徴とする請求項1に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
  7. 前記S3では、ステレオカメラ座標系における姿勢観測は、
    ステレオカメラ画像データに対して特徴点抽出を行い、特徴点を中心として円領域を構成し、

    θ=arctan(m01/m10
    [式中、Cは円領域の質量中心を表し、θは特徴点の方向ベクトルを表し、mpqは円領域のモーメントを表し、

    (式中、Rは円領域の半径を表し、x,yはx軸座標、y軸座標を表し、I(x,y)はグレースケール方程式を表す)で定義される]
    連続した複数フレームのステレオカメラ画像データの特徴点を抽出してマッチングすることにより、マッチングされた画素点を用いてPnP問題を設定し、ステレオカメラの回転行列R及び並進ベクトルTを得る、ことを特徴とする請求項6に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
  8. 前記S3では、ステレオカメラ座標系における3次元点群生成は、
    同一フレームのステレオカメラの左右カメラ画像に対して上記特徴点を抽出してマッチングし、グレースケール二乗和誤差アルゴリズムに基づいて視差計算を行い、

    (式中、x,y,dはそれぞれx軸座標、y軸座標、視差であり、i,jはそれぞれx軸方向及びy軸方向の変化値であり、m,nはそれぞれx軸方向及びy軸方向の最大値であり、I(x,y),I(x,y)はグレースケール方程式を表す)
    視差及び元の座標により3次元点群データを生成し、3次元座標が、

    [式中、x,xはそれぞれ左右カメラに対応する横座標値であり、y,yはそれぞれ左右カメラの縦座標値であり、f,fはそれぞれ左右カメラの内部パラメータの対応する焦点距離であり、X,Y,Zはそれぞれ3次元座標であり、Dは深度値であり、下記式:
    D=Bf/d
    (式中、Bは基線長であり、fはカメラの焦点距離であり、dは左右画像の視差である)で算出される]で表される、ことを特徴とする請求項7に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
  9. 前記S4は、損傷領域の位置に基づいて、排水システムが損傷領域を被覆して排水して乾燥空間を形成するように、水中移動プラットフォームの移動軌跡を計画し、レーザ駆動システムを使用して、レーザセンサを乾燥空間で3次元走査を行うように制御し、レーザセンサデータと駆動用IMUデータとを通信システムを介してホストコンピュータに送信し、ホストコンピュータは、S2で得られた外部パラメータ行列に基づいて、駆動用IMUデータを用いてレーザ駆動システムの姿勢を取得し、変換してレーザセンサの位置を取得し、レーザセンサの位置及び点群データに基づいて、精細なレーザセンサデータの3次元再構成を得、3次元再構成結果に基づいて損傷位置を検出する、ことを特徴とする請求項1に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
  10. 前記S4では、レーザセンサデータの3次元再構成は、
    レーザセンサは、固定周波数でレーザパルスを発射し、受信器により戻ってきた反射光を受光して距離を判断するとともに、反射強度に基づいて対象材質を大まかに区別し、距離測定式が、
    L=tc/2
    (式中、Lは目標距離であり、tは戻り時間であり、cは光速である)であり、
    駆動用IMUを使用してレーザセンサの位置姿勢予測を行い、回転行列R及び並進ベクトルTに基づいてレーザセンサの3次元再構成結果を得る、ことを特徴とする請求項9に記載の視覚とIMUとの融合に基づく海洋機器の水中損傷の3次元再構成方法。
JP2025512740A 2022-08-31 2023-08-30 視覚とimuとの融合に基づく海洋機器の水中損傷の3次元再構成方法 Active JP7841786B2 (ja)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN202211051347.X 2022-08-31
CN202211051347.XA CN115471570B (zh) 2022-08-31 2022-08-31 基于视觉和imu融合的海洋装备水下损伤三维重建方法
PCT/CN2023/115908 WO2024046390A1 (zh) 2022-08-31 2023-08-30 基于视觉和imu融合的海洋装备水下损伤三维重建方法

Publications (2)

Publication Number Publication Date
JP2025534947A JP2025534947A (ja) 2025-10-22
JP7841786B2 true JP7841786B2 (ja) 2026-04-07

Family

ID=84368940

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2025512740A Active JP7841786B2 (ja) 2022-08-31 2023-08-30 視覚とimuとの融合に基づく海洋機器の水中損傷の3次元再構成方法

Country Status (3)

Country Link
JP (1) JP7841786B2 (ja)
CN (1) CN115471570B (ja)
WO (1) WO2024046390A1 (ja)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115471570B (zh) * 2022-08-31 2025-06-20 华南理工大学 基于视觉和imu融合的海洋装备水下损伤三维重建方法
CN118298112B (zh) * 2024-06-04 2024-09-20 华南理工大学 一种基于悬浮式云台的水下目标三维重建方法
CN119024382B (zh) * 2024-09-18 2025-05-27 中国热带农业科学院农业机械研究所 基于边缘计算的异物收集装备定位轨迹模拟方法及系统
CN119354037A (zh) * 2024-10-16 2025-01-24 杭州电子科技大学 一种长基线深海双目测量相机
CN119693559B (zh) * 2024-12-24 2025-09-26 大连理工大学 一种基于双目视觉的海洋波浪场重构方法和装置
CN120194609B (zh) * 2025-03-20 2026-01-30 华中科技大学 一种基于六维传感器的大型构件测量系统及其标定方法
CN120846235B (zh) * 2025-09-22 2025-11-21 武汉科技大学 基于激光断面轮廓差动扫描的高动态形变测量系统与方法
CN121617074B (zh) * 2026-02-03 2026-04-03 上海海事大学 纹理增强轨迹与视觉slam融合的桥梁防撞预警方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170061632A1 (en) 2015-08-24 2017-03-02 Qualcomm Incorporated Systems and methods for depth map sampling
JP2018062810A (ja) 2016-10-14 2018-04-19 株式会社Jmuアムテック 水中調査装置および方法
JP2018151964A (ja) 2017-03-14 2018-09-27 朝日航洋株式会社 桟橋上部工下面の検査装置、桟橋上部工下面の検査システムおよび桟橋上部工下面の検査方法
JP2020507072A (ja) 2017-01-27 2020-03-05 カールタ インコーポレイテッド 実時間オンライン自己運動推定を備えたレーザスキャナ
CN113570715A (zh) 2021-07-23 2021-10-29 东北大学 基于传感器融合的旋转激光实时定位建模系统及方法
CN113777580A (zh) 2021-07-19 2021-12-10 长春理工大学 一种基于可旋转单线激光雷达的三维激光扫描仪
JP2022023592A (ja) 2020-07-27 2022-02-08 株式会社トプコン 測量システム、測量方法および測量用プログラム

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9430822B2 (en) * 2013-06-14 2016-08-30 Microsoft Technology Licensing, Llc Mobile imaging platform calibration
CN110763152B (zh) * 2019-10-09 2021-08-20 哈尔滨工程大学 一种水下主动旋转结构光三维视觉测量装置及测量方法
CN111951305B (zh) * 2020-08-20 2022-08-23 重庆邮电大学 一种基于视觉和激光雷达的目标检测和运动状态估计方法
CN112833816A (zh) * 2020-12-31 2021-05-25 武汉中观自动化科技有限公司 一种标志点定位与智能反向定位混合的定位方法和系统
CN113744337B (zh) * 2021-09-07 2023-11-24 江苏科技大学 一种融合视觉、imu与声纳的同步定位与建图方法
CN115471570B (zh) * 2022-08-31 2025-06-20 华南理工大学 基于视觉和imu融合的海洋装备水下损伤三维重建方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170061632A1 (en) 2015-08-24 2017-03-02 Qualcomm Incorporated Systems and methods for depth map sampling
JP2018062810A (ja) 2016-10-14 2018-04-19 株式会社Jmuアムテック 水中調査装置および方法
JP2020507072A (ja) 2017-01-27 2020-03-05 カールタ インコーポレイテッド 実時間オンライン自己運動推定を備えたレーザスキャナ
JP2018151964A (ja) 2017-03-14 2018-09-27 朝日航洋株式会社 桟橋上部工下面の検査装置、桟橋上部工下面の検査システムおよび桟橋上部工下面の検査方法
JP2022023592A (ja) 2020-07-27 2022-02-08 株式会社トプコン 測量システム、測量方法および測量用プログラム
CN113777580A (zh) 2021-07-19 2021-12-10 长春理工大学 一种基于可旋转单线激光雷达的三维激光扫描仪
CN113570715A (zh) 2021-07-23 2021-10-29 东北大学 基于传感器融合的旋转激光实时定位建模系统及方法

Also Published As

Publication number Publication date
JP2025534947A (ja) 2025-10-22
WO2024046390A1 (zh) 2024-03-07
CN115471570A (zh) 2022-12-13
CN115471570B (zh) 2025-06-20

Similar Documents

Publication Publication Date Title
JP7841786B2 (ja) 視覚とimuとの融合に基づく海洋機器の水中損傷の3次元再構成方法
CN110728715B (zh) 一种智能巡检机器人摄像机角度自适应调整方法
CN109283538B (zh) 一种基于视觉和激光传感器数据融合的海上目标大小检测方法
CN110033489B (zh) 一种车辆定位准确性的评估方法、装置及设备
US9007601B2 (en) Automatic measurement of dimensional data with a laser tracker
CN106289184B (zh) 一种无gnss信号和无控制点下无人机协同视觉形变监测方法
CN109859271B (zh) 一种水下摄像机与前视声纳的联合标定方法
CN112258590A (zh) 一种基于激光的深度相机外参标定方法、设备及其存储介质
CN105844692B (zh) 基于双目立体视觉的三维重建装置、方法、系统及无人机
CN109242918A (zh) 一种直升机机载双目立体视觉标定方法
CN109410234A (zh) 一种基于双目视觉避障的控制方法及控制系统
CN118628573B (zh) 一种针对水下视觉退化场景的机器人在线定位和建图方法
EP4257924A1 (en) Laser scanner for verifying positioning of components of assemblies
CN111640177A (zh) 一种基于水下声呐探测的三维建模方法及无人潜水器
Yin et al. Study on underwater simultaneous localization and mapping based on different sensors
CN116129669A (zh) 一种基于激光雷达的车位测评方法、系统、设备及介质
JP2000028332A (ja) 3次元計測装置及び3次元計測方法
CN109798877A (zh) 一种仿生水下机器鱼双目立体测距方法
CN115876099A (zh) 一种基于倾斜预警式宽范围双目相机船舶高度测量系统
CN120807603A (zh) 一种超高压输电线路杆塔倾斜度计算方法、装置及介质
CN114705682A (zh) 一种岩体结构智能视觉探测识别成像装置
CN120510590A (zh) 基于激光雷达与摄像头数据融合的码头感知方法及系统
CN110533702B (zh) 一种基于光场多层折射模型的水下双目视觉立体匹配方法
CN118671773A (zh) 一种基于双声纳数据融合的水下三维重建装置及方法
CN116086431A (zh) 一种移动机器人地图构建方法及装置

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20250327

TRDD Decision of grant or rejection written
A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20260227

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20260303

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20260318

R150 Certificate of patent or registration of utility model

Ref document number: 7841786

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150