JPWO2021044628A1 - Deformity detection system and deformation detection method - Google Patents
Deformity detection system and deformation detection method Download PDFInfo
- Publication number
- JPWO2021044628A1 JPWO2021044628A1 JP2019567758A JP2019567758A JPWO2021044628A1 JP WO2021044628 A1 JPWO2021044628 A1 JP WO2021044628A1 JP 2019567758 A JP2019567758 A JP 2019567758A JP 2019567758 A JP2019567758 A JP 2019567758A JP WO2021044628 A1 JPWO2021044628 A1 JP WO2021044628A1
- Authority
- JP
- Japan
- Prior art keywords
- unit
- data
- point cloud
- cloud data
- dimensional
- 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.)
- Granted
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 232
- 238000006243 chemical reaction Methods 0.000 claims abstract description 81
- 239000011159 matrix material Substances 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims description 169
- 239000013598 vector Substances 0.000 claims description 123
- 238000011156 evaluation Methods 0.000 claims description 17
- 238000000034 method Methods 0.000 description 57
- 230000008569 process Effects 0.000 description 34
- 238000010586 diagram Methods 0.000 description 27
- 238000004891 communication Methods 0.000 description 23
- 230000006870 function Effects 0.000 description 22
- 238000011161 development Methods 0.000 description 15
- 230000018109 developmental process Effects 0.000 description 15
- 230000015654 memory Effects 0.000 description 15
- 230000009466 transformation Effects 0.000 description 14
- 230000008859 change Effects 0.000 description 10
- 238000005259 measurement Methods 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 4
- 238000012935 Averaging Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000012447 hatching Effects 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000006866 deterioration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000008685 targeting Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/20—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/50—Depth or shape recovery
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Theoretical Computer Science (AREA)
- Length Measuring Devices By Optical Means (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
- Image Analysis (AREA)
Abstract
変状検出システム(100)は、取得部(21)と、データ変換部(22)と、段差検出部(23)と、を備える。取得部(21)は、構造物の表面を3次元点群で表すデータである3次元点群データと構造物の形状を示す形状データとを取得する。データ変換部(22)は、形状データに基づいて構造物の表面に対応する基準面を設定し、3次元点群データを、基準面に行列状に配列される複数の点から選択される1つの点の位置と基準面からのずれ量とで各々表される複数の要素点のデータを含む要素点群データへ変換する。段差検出部(23)は、要素点群データに含まれる複数の要素点に対してガウシアン差分フィルタをかけて得られる結果に基づいて、構造物の表面における段差の位置を検出する。The deformation detection system (100) includes an acquisition unit (21), a data conversion unit (22), and a step detection unit (23). The acquisition unit (21) acquires three-dimensional point cloud data, which is data representing the surface of the structure as a three-dimensional point cloud, and shape data indicating the shape of the structure. The data conversion unit (22) sets a reference plane corresponding to the surface of the structure based on the shape data, and selects the three-dimensional point cloud data from a plurality of points arranged in a matrix on the reference plane. It is converted into element point cloud data including data of a plurality of element points represented by the position of one point and the amount of deviation from the reference plane. The step detection unit (23) detects the position of the step on the surface of the structure based on the result obtained by applying a Gaussian difference filter to a plurality of element points included in the element point cloud data.
Description
本発明は、構造物の表面に生じる段差を検出することができる変状検出システムおよび変状検出方法に関する。 The present invention relates to a deformation detection system and a deformation detection method capable of detecting a step generated on the surface of a structure.
トンネル、橋梁、または道路などの構造物の点検を効率的に行うため、MMS(Mobile Mapping System)と呼ばれる走行型計測装置が用いられている。MMSは、走行しながら周囲の構造物をレーザスキャナ装置などで計測し、計測した結果を示すデータを生成する。かかるデータには、構造物の表面を複数の3次元点で表す3次元点群のデータである3次元点群データが含まれており、3次元点群データは構造物の点検に用いられる。構造物の点検では、例えば、構造物の表面に生じるひびまたは段差などの変状が検出される。 In order to efficiently inspect structures such as tunnels, bridges, and roads, a traveling type measuring device called MMS (Mobile Mapping System) is used. The MMS measures the surrounding structure with a laser scanner device or the like while traveling, and generates data showing the measurement result. Such data includes three-dimensional point cloud data which is data of a three-dimensional point cloud representing the surface of the structure by a plurality of three-dimensional points, and the three-dimensional point cloud data is used for inspection of the structure. In the inspection of the structure, for example, deformation such as cracks or steps generated on the surface of the structure is detected.
例えば、特許文献1には、堤防の表面を複数の3次元点で表す3次元点群データに基づいて、堤防の表面に形成された段差を検出する技術が提案されている。特許文献1に記載の技術は、堤防の表面を表す複数の3次元点のうち、近傍空間に他の3次元点が予め定められた数以上包含される3次元点を、段差に由来する特徴点として抽出し、抽出した特徴点の位置から段差の位置を検出する。
For example,
しかしながら、特許文献1に記載の技術は、近傍空間に他の3次元点が予め定められた数以上包含されることを条件として3次元点を、段差に由来する特徴点として抽出する。そのため、構造物の表面において段差がある領域が上記条件を満たさない場合、構造物の表面に生じる段差の位置を精度よく検出することができない。
However, the technique described in
本発明は、上記に鑑みてなされたものであって、構造物の表面を表す3次元点群から構造物の表面に生じる段差の位置を精度よく検出することができる変状検出システムを得ることを目的とする。 The present invention has been made in view of the above, and obtains a deformation detection system capable of accurately detecting the position of a step generated on the surface of a structure from a three-dimensional point cloud representing the surface of the structure. With the goal.
上述した課題を解決し、目的を達成するために、本発明の変状検出システムは、取得部と、データ変換部と、段差検出部と、を備える。取得部は、構造物の表面を3次元点群で表すデータである3次元点群データと構造物の形状を示す形状データとを取得する。データ変換部は、形状データに基づいて構造物の表面に対応する基準面を設定し、3次元点群データを、基準面に行列状に配列される複数の点から選択される1つの点の位置と基準面からのずれ量とで各々表される複数の要素点のデータを含む要素点群データへ変換する。段差検出部は、要素点群データに含まれる複数の要素点に対してガウシアン差分フィルタをかけて得られる結果に基づいて、構造物の表面における段差の位置を検出する。 In order to solve the above-mentioned problems and achieve the object, the deformation detection system of the present invention includes an acquisition unit, a data conversion unit, and a step detection unit. The acquisition unit acquires three-dimensional point cloud data, which is data representing the surface of the structure as a three-dimensional point cloud, and shape data indicating the shape of the structure. The data conversion unit sets a reference plane corresponding to the surface of the structure based on the shape data, and sets the three-dimensional point cloud data of one point selected from a plurality of points arranged in a matrix on the reference plane. It is converted into element point cloud data including data of a plurality of element points represented by the position and the amount of deviation from the reference plane. The step detection unit detects the position of the step on the surface of the structure based on the result obtained by applying a Gaussian difference filter to a plurality of element points included in the element point cloud data.
本発明によれば、構造物の表面を表す3次元点群から構造物の表面に生じる段差の位置を精度よく検出することができる、という効果を奏する。 According to the present invention, there is an effect that the position of a step generated on the surface of the structure can be accurately detected from the three-dimensional point cloud representing the surface of the structure.
以下に、本発明の実施の形態にかかる変状検出システムおよび変状検出方法を図面に基づいて詳細に説明する。なお、この実施の形態によりこの発明が限定されるものではない。 Hereinafter, the deformation detection system and the deformation detection method according to the embodiment of the present invention will be described in detail with reference to the drawings. The present invention is not limited to this embodiment.
実施の形態1.
図1は、本発明の実施の形態1にかかる変状検出システムの構成例を示す図である。図1に示すように、実施の形態1にかかる変状検出システム100は、構造物の表面を3次元点群で表す3次元点群のデータである3次元点群データDpgを生成する走行型計測装置1と、3次元点群データDpgに基づいて、構造物の表面に生じる段差を検出する段差検出装置2とを備える。以下においては、上記構造物の表面がトンネル5の内壁面5aであるものとして説明するが、構造物の表面はトンネル5の内壁面5aに限定されず、橋梁の表面、道路の表面、または建物の表面などであってもよい。
FIG. 1 is a diagram showing a configuration example of a deformation detection system according to a first embodiment of the present invention. As shown in FIG. 1, the
走行型計測装置1は、不図示のレーザスキャナ装置を有しており、走行しながらレーザスキャナ装置によって得られるデータからトンネル5の内壁面5aの3次元点群データDpgを生成する。走行型計測装置1は、MMSとも呼ばれる。3次元点群データDpgには、複数の3次元点の計測データが含まれており、各3次元点の計測データには、3次元直交座標系における3次元点の位置を示すデータが含まれる。
The traveling
段差検出装置2は、通信部10と、記憶部11と、処理部12とを備える。通信部10は、ネットワーク3を介して走行型計測装置1と通信可能に接続される。ネットワーク3は、例えば、インターネットなどのWAN(Wide Area Network)である。通信部10は、走行型計測装置1からネットワーク3を介して送信される3次元点群データDpgを受信する。
The
記憶部11は、トンネル5の形状データDfを記憶している。形状データDfは、例えば、トンネル5の内壁面5aの断面形状を示すデータである。内壁面5aの断面形状は、トンネル軸方向に直交する面でトンネル5の内壁を切断したと仮定した場合における内壁面5aを示す線の形状である。トンネル軸方向は、トンネル5が延伸する方向である。
The storage unit 11 stores the shape data Df of the
処理部12は、取得部21と、データ変換部22と、段差検出部23とを備える。取得部21は、通信部10を介して走行型計測装置1から3次元点群データDpgを取得し、記憶部11から形状データDfを取得する。データ変換部22は、取得部21によって取得された形状データDfに基づいてトンネル5の内壁面5aに対応する基準面を設定し、設定した基準面を用いて、取得部21によって取得された3次元点群データDpgを要素点群データDeへ変換する。要素点群データDeは、基準面に行列状に配列される複数の点から選択される1つの点の位置と基準面からのずれ量とで各々表される複数の要素点のデータを含む。
The
基準面に行列状に配列される複数の点は、基準面に格子を形成したと仮定した場合の格子点と言え、以下、基準面に行列状に配列される複数の点の各々を格子点と記載する場合がある。行方向に隣接する2つの格子点の間隔と、列方向に隣接する2つの格子点の間隔とは同じである。また、行方向および列方向の一方は、トンネル軸方向であり、他方はトンネル周方向である。なお、トンネル周方向は、トンネル軸方向に直交する面でトンネル5の内壁面5aを切断した場合の内壁面5aを示す線に沿った方向である。
A plurality of points arranged in a matrix on the reference plane can be said to be grid points when it is assumed that a grid is formed on the reference plane. Hereinafter, each of the plurality of points arranged in a matrix on the reference plane is a grid point. May be described as. The distance between two grid points adjacent to each other in the row direction is the same as the distance between two grid points adjacent to each other in the column direction. One of the row direction and the column direction is the tunnel axial direction, and the other is the tunnel circumferential direction. The tunnel circumferential direction is a direction along a line indicating the
段差検出部23は、要素点群データDeに含まれる複数の要素点に対してガウシアン差分フィルタをかけて得られる結果に基づいて、内壁面5aに生じている段差の位置を検出する。例えば、段差検出部23は、ガウシアン差分フィルタをかけた後の複数の要素点のうち隣接する2つの要素点同士を結ぶ線が基準面に交差する位置を特定し、特定した位置を内壁面5aに生じている段差の位置として検出する。
The
図2は、実施の形態1にかかる変状検出システムによる計測処理およびデータ変換処理の流れを説明するための図である。図2に示すように、走行型計測装置1は、トンネル5の内壁面5aを計測し、3次元点群データDpgを生成する(ステップS1)。段差検出装置2の取得部21は、走行型計測装置1から送信される3次元点群データDpgを取得する(ステップS2)。
FIG. 2 is a diagram for explaining a flow of measurement processing and data conversion processing by the deformation detection system according to the first embodiment. As shown in FIG. 2, the traveling
段差検出装置2のデータ変換部22は、記憶部11に記憶された形状データDfに基づいてトンネル5の内壁面5aに対応する基準面Swrを設定し、取得部21で取得された3次元点群データDpgを要素点群データDeへ変換するデータ変換処理を行う(ステップS3)。要素点群データDeは、基準面Swrに行列状に配列される複数の格子点Pgから選択される1つの格子点Pgの位置と基準面Swrからのずれ量とで各々表される複数の要素点Peのデータを含むデータである。基準面Swrは、形状データDfに基づいて得られるトンネル5の内壁面5aを示す面である。基準面Swrは、トンネル5の内壁面5aに凹凸がないと仮定した場合のトンネル5の内壁面5aであり、仮想的な面である。
The
図2において、要素点群データDeを表す図では、要素点Peの位置を分かりやすくするために、トンネル軸方向で隣接する要素点Pe間とトンネル周方向で隣接する要素点Pe間とが各々破線で接続されている。また、図2に示す基準面Swrは、トンネル軸方向における位置とトンネル周方向における位置と基準面Swrからのずれ量とを各々異なる軸とする直交座標系である展開図座標系の平面へ変換された後の基準面である。以下においては、説明の便宜上、展開図座標系へ変換前の基準面Swrを基準面Swr’と記載する場合がある。 In FIG. 2, in the diagram showing the element point cloud data De, in order to make the position of the element point Pe easy to understand, the element points Pe adjacent to each other in the tunnel axial direction and the element points Pe adjacent to each other in the tunnel circumferential direction are respectively. It is connected by a broken line. Further, the reference plane Swr shown in FIG. 2 is converted into a plane of the development view coordinate system which is an orthogonal coordinate system in which the position in the tunnel axial direction, the position in the tunnel circumferential direction, and the amount of deviation from the reference plane Swr are different axes. It is a reference plane after being made. In the following, for convenience of explanation, the reference plane Swr before conversion to the developed view coordinate system may be described as the reference plane Swr'.
ステップS3において、データ変換部22は、3次元点群データDpgに含まれる各3次元点の基準面Swr’に対する位置を算出する。算出される3次元点の位置は、トンネル軸方向の位置とトンネル周方向の位置と基準面Swr’からのずれ量とで表される。トンネル軸方向の位置およびトンネル周方向の位置は、基準面Swr’上の位置である。また、基準面Swr’からのずれ量は、基準面Swr’からの距離で大きさが表され、基準面Swr’と直交する方向であって基準面Swr’から見た方向が極性として表される。
In step S3, the
次に、データ変換部22は、基準面Swr’に対する各3次元点の位置と、展開図座標系へ変換後の基準面Swrに行列状に配列される複数の格子点Pgとに基づいて、3次元点群データDpgを展開図座標系の要素点群データDeへ変換する。
Next, the
例えば、データ変換部22は、3次元点群データDpgを展開図座標系のデータである展開図点群データDpg1へ変換する。データ変換部22は、3次元点群データDpgに含まれる各3次元点のトンネル周方向の位置を、展開図座標系においてトンネル周方向を示す一つの軸上の位置へ変換することによって、3次元点群データDpgを展開図点群データDpg1へ変換する。
For example, the
データ変換部22は、展開図点群データDpg1に含まれる複数の3次元点のうちトンネル軸方向における位置とトンネル周方向における位置とが1つの格子点Pgの位置から予め定められた範囲内にある1つ以上の3次元点を1つの格子点Pgに割り当てる。データ変換部22は、割り当てた1つの格子点Pgのトンネル軸方向における位置とトンネル周方向における位置とを、要素点Peのトンネル軸方向における位置とトンネル周方向における位置に決定する。
The
また、データ変換部22は、1つの格子点Pgに割り当てられた1つ以上の3次元点の基準面Swrからのずれ量から、要素点Peにおける基準面Swrからのずれ量を算出する。例えば、データ変換部22は、1つの格子点Pgに割り当てられた1つ以上の3次元点の基準面Swrからのずれ量を平均化して得られる値を要素点Peにおける基準面Swrからのずれ量にすることができる。
Further, the
データ変換部22は、このように1つ以上の3次元点を1つの格子点Pgに割り当てて要素点Peを算出する処理を格子点Pg毎に行うことで、展開図点群データDpg1を要素点群データDeへ変換する。
The
図3は、実施の形態1にかかる変状検出システムによる段差検出処理の流れを説明するための図である。図3において、要素点群データDeを表す図では、図2に示す要素点群データDeと同様に、要素点Peの位置を分かりやすくするために、要素点Pe間が破線で接続されている。 FIG. 3 is a diagram for explaining the flow of the step detection process by the deformation detection system according to the first embodiment. In the figure showing the element point cloud data De in FIG. 3, similarly to the element point cloud data De shown in FIG. 2, the element points Pe are connected by a broken line in order to make the position of the element point Pe easy to understand. ..
図3に示すように、段差検出部23は、データ変換部22によって変換された要素点群データDeに含まれる複数の要素点Peに対してガウシアン差分フィルタをかける処理であるガウシアン差分処理を行う。ガウシアン差分フィルタは、DoG(Difference of Gaussian)とも呼ばれる。ガウシアン差分処理には、第1フィルタ処理(ステップS4)と、第2フィルタ処理(ステップS5)と、差分演算処理(ステップS6)とが含まれる。
As shown in FIG. 3, the
段差検出部23は、ステップS4の第1フィルタ処理として、要素点群データDeに含まれる複数の要素点Peに対して第1ガウシアンフィルタをかける処理を行う。また、段差検出部23は、ステップS5の第2フィルタ処理として、要素点群データDeに含まれる複数の要素点Peに対して第2ガウシアンフィルタをかける処理を行う。第1ガウシアンフィルタおよび第2ガウシアンフィルタによって基準面Swrからの要素点Peのずれ量に対する平滑化処理が行われる。
The
第1ガウシアンフィルタと第2ガウシアンフィルタとは、互いに特性が異なるガウシアンフィルタであり、少なくとも一部のパラメータの値が互いに異なる。段差検出部23は、ステップS6の差分演算処理において、第1フィルタ処理を行って得られる第1ガウシアンフィルタの出力と第2フィルタ処理を行って得られる第2ガウシアンフィルタの出力との差分を算出する。
The first Gaussian filter and the second Gaussian filter are Gaussian filters having different characteristics from each other, and the values of at least some parameters are different from each other. The
複数の要素点Peに対してガウシアン差分フィルタをかけることで、複数の要素点Peの各々の基準面Swrからのずれ量の差分がガウシアン差分フィルタから出力される。基準面Swrからのずれ量の差分は、要素点Peにおける基準面Swrからのずれ量に対して第1ガウシアンフィルタをかけて得られる基準面Swrからのずれ量と、要素点Peにおける基準面Swrからのずれ量に対して第2ガウシアンフィルタをかけて得られる基準面Swrからのずれ量との差分である。 By applying the Gaussian difference filter to the plurality of element points Pe, the difference in the amount of deviation of the plurality of element points Pe from the reference plane Swr is output from the Gaussian difference filter. The difference in the amount of deviation from the reference surface Swr is the amount of deviation from the reference surface Swr obtained by applying the first Gaussian filter to the amount of deviation from the reference surface Swr at the element point Pe and the reference surface Swr at the element point Pe. It is a difference from the deviation amount from the reference plane Swr obtained by applying the second Gaussian filter to the deviation amount from.
基準面Swrからのずれ量の差分は、例えば、要素点Peに対して第2ガウシアンフィルタをかけて得られる基準面Swrからのずれ量から要素点Peに対して第1ガウシアンフィルタをかけて得られる基準面Swrからのずれ量を減算して得られる。ガウシアン差分フィルタをかけた後の要素点Peの位置は、トンネル軸方向の位置と、トンネル周方向の位置と、基準面Swrからのずれ量の差分とで表される。 The difference in the amount of deviation from the reference surface Swr is obtained, for example, by applying the first Gaussian filter to the element point Pe from the amount of deviation from the reference surface Swr obtained by applying the second Gaussian filter to the element point Pe. It is obtained by subtracting the amount of deviation from the reference plane Swr. The position of the element point Pe after applying the Gaussian difference filter is represented by the position in the tunnel axial direction, the position in the tunnel circumferential direction, and the difference in the amount of deviation from the reference plane Swr.
次に、段差検出部23は、要素点群データDeの複数の要素点Peに対してガウシアン差分フィルタをかけて得られる結果に基づいて、段差の位置を検出する(ステップS7)。ステップS7において、段差検出部23は、差分演算処理の結果として得られる要素点Pe毎の基準面Swrからのずれ量の差分に基づいて、零交差が生じる位置を特定する。零交差が生じる位置は、基準面Swrからのずれ量の差分が零になる位置であり、トンネル軸方向の位置とトンネル周方向の位置とで表される。
Next, the
具体的には、段差検出部23は、トンネル軸方向の位置とトンネル周方向の位置と基準面Swrからのずれ量の差分とを各々異なる軸とする直交座標系で要素点Peを扱う。段差検出部23は、トンネル軸方向またはトンネル周方向で隣接し且つガウシアン差分フィルタをかけて得られる基準面Swrからのずれ量の差分値の極性が異なる2つの要素点Pe同士を結ぶ線が零交差する位置である零交差点を特定する。零交差点は、トンネル軸方向の位置とトンネル周方向の位置とで表される。段差検出部23は、零交差点の位置をトンネル5の内壁面5aに生じている段差の位置として検出する。
Specifically, the
走行型計測装置1がトンネル5内を走行しながらトンネル5の内壁面5aをレーザスキャナ装置で計測する場合、トンネル5の内壁面5aに対して走行型計測装置1のレーザスキャナ装置からのレーザが正対して照射される。そのため、トンネル5の内壁面5aに生じる段差に対してレーザの照射方向が平行になり、走行型計測装置1の計測によって得られるトンネル5の3次元点群において段差による3次元点の疎密差が現われにくい。また、トンネル5の内壁面5aにおける段差は、浮きまたは剥離に伴う段差でありその高さは相対的に低いことから、このことからも、トンネル5の内壁面5aを表す3次元点群に段差による3次元点の疎密差が生じにくい。
When the traveling
段差検出装置2の処理部12は、複数の要素点Peに対してガウシアン差分フィルタをかけて得られる結果に基づいて、トンネル5の内壁面5aにおける段差の位置を検出する。そのため、処理部12は、トンネル5の内壁面5aを表す3次元点群において段差による3次元点の疎密差が少ない場合であっても、トンネル5の内壁面5aに生じる段差の位置を精度よく検出することができる。
The
また、トンネル軸方向およびトンネル周方向の各々で要素点Peは等間隔に配置され、各要素点Peに対して1以上の3次元点が割り当てられる。そのため、複数の要素点Peがトンネル軸方向およびトンネル周方向の各々で等間隔に配置される。これにより、処理部12は、トンネル5の内壁面5aの全体に亘って段差を同じ精度で検出することができる。
Further, the element points Pe are arranged at equal intervals in each of the tunnel axial direction and the tunnel circumferential direction, and one or more three-dimensional points are assigned to each element point Pe. Therefore, a plurality of element points Pe are arranged at equal intervals in each of the tunnel axial direction and the tunnel circumferential direction. As a result, the
また、処理部12は、複数の3次元点を一つの要素点Peに割り当てることができるため、3次元点群データDpgで表される3次元点群に含まれる3次元点の数に比べ、要素点群データDeで表される要素点群に含まれる要素点Peの数を減らすことができる。そのため、ガウシアン差分フィルタの処理における処理部12の負荷を軽減することができる。
Further, since the
以下、実施の形態1にかかる変状検出システム100についてさらに詳細に説明する。図1に示す走行型計測装置1は、レーザスキャナ装置に加え、2次元画像を撮像する撮像部を有しており、走行しながら撮像部によって得られる撮像データからトンネル5の内壁面5aの画像データDimgを生成する。画像データDimgは、撮像部によって撮像された画像である撮像画像を繋ぎ合せて得られるトンネル5の内壁面5aの全体画像を示すデータである。
Hereinafter, the
走行型計測装置1は、トンネル5の内壁面5aを定期的に計測し、計測を行う毎に3次元点群データDpgと画像データDimgを生成する。走行型計測装置1は、トンネル5の内壁面5aを計測する毎に、3次元点群データDpgと画像データDimgを段差検出装置2へネットワーク3を介して送信する。なお、3次元点群データDpgと画像データDimgとには各々走行型計測装置1によって計測された時刻である計測時刻の情報が含まれる。
The traveling
次に、段差検出装置2の構成についてさらに具体的に説明する。図4は、図1に示す段差検出装置をさらに詳細化した構成例を示す図である。図4に示すように、段差検出装置2は、図1に示す段差検出装置2の構成に加え、入力部13と、表示部14とを備える。
Next, the configuration of the
通信部10は、図1に示す走行型計測装置1によって生成されたデータをネットワーク3経由で受信する。具体的には、通信部10は、走行型計測装置1からネットワーク3経由で3次元点群データDpgと画像データDimgを受信する。
The
記憶部11は、形状データDfと、3次元点群データDpgと、画像データDimgと、変状データDdfとを記憶する。形状データDfは、トンネル5の断面形状を示すデータである。形状データDfは、設計図から得られる内壁面5aの断面形状のデータに代えて、3次元点群データDpgから得られる内壁面5aの断面形状のデータであってもよい。
The storage unit 11 stores the shape data Df, the three-dimensional point cloud data Dpg, the image data Dimg, and the deformation data Ddf. The shape data Df is data indicating the cross-sectional shape of the
処理部12は、例えば、3次元点群データDpgで示される内壁面5aの3次元点群に柱体を当てはめて内壁面5aの断面形状を推定することによって形状データDfを生成することができる。また、処理部12は、3次元点群データDpgに含まれる複数の3次元点のうちトンネル軸方向に直交する面である直交面から設定された距離内にある複数の3次元点を直交面に投影し、投影した複数の点を曲線近似することによって形状データDfを生成することもできる。
The
また、形状データDfは、トンネル5の内壁面5aの全体形状を示すデータであってもよい。処理部12は、例えば、3次元点群データDpgで示される内壁面5aの3次元点群に柱体を当てはめて内壁面5aの全体形状を推定することによって形状データDfを生成することができる。
Further, the shape data Df may be data indicating the overall shape of the
3次元点群データDpgおよび画像データDimgは、上述したように、走行型計測装置1によって生成され通信部10で受信されるデータであり、上述した時刻情報を含む。
As described above, the three-dimensional point cloud data Dpg and the image data Dimg are data generated by the traveling
変状データDdfは、内壁面5aに生じた変状の種別、範囲、および時刻などを示すデータを含む。変状の種別には、亀裂、段差、剥離、および漏水などがある。変状の範囲は、内壁面5aにおける変状の範囲であり、例えば、基準面Swr’における範囲または基準面Swrにおける範囲である。変状の時刻は、変状が検出された時刻であり、例えば、変状が検出された3次元点群データDpgに含まれる時刻情報で示される時刻である。
The deformation data Ddf includes data indicating the type, range, time, and the like of the deformation that occurred on the
処理部12は、取得部21と、データ変換部22と、段差検出部23と、段差評価部24と、属性設定部25と、座標変換部26と、表示処理部27とを備える。取得部21は、走行型計測装置1からネットワーク3を介して3次元点群データDpgおよび画像データDimgを取得し、取得した3次元点群データDpgおよび画像データDimgを記憶部11に記憶させる。なお、取得部21は、3次元点群データDpgおよび画像データDimgが記録された不図示の記録媒体から3次元点群データDpgおよび画像データDimgを取得することもできる。
The
データ変換部22は、形状データDfに基づいて内壁面5aに対応する基準面Swrを設定し、3次元点群データDpgを要素点群データDeへ変換する。要素点群データDeは、上述したように、基準面Swrに行列状に配列される複数の格子点Pgから選択される1つの格子点Pgの位置と基準面Swrからのずれ量とで各々表される複数の要素点Peのデータを含む。なお、データ変換部22は、3次元点群データDpgから形状データDfを生成する上述した処理を行うことができる。
The
図5は、実施の形態1にかかるデータ変換部による3次元点群データの要素点群データへの変換方法を説明するための図である。図6は、実施の形態1にかかるデータ変換部による要素点における基準面からのずれ量の第1の決定方法を説明するための図である。図7は、実施の形態1にかかるデータ変換部による要素点における基準面からのずれ量の第2の決定方法を説明するための図である。図8は、実施の形態1にかかるデータ変換部による要素点における基準面からのずれ量の第3の決定方法を説明するための図である。 FIG. 5 is a diagram for explaining a method of converting the three-dimensional point cloud data into the element point cloud data by the data conversion unit according to the first embodiment. FIG. 6 is a diagram for explaining a first method of determining the amount of deviation from the reference plane at the element point by the data conversion unit according to the first embodiment. FIG. 7 is a diagram for explaining a second method of determining the amount of deviation from the reference plane at the element point by the data conversion unit according to the first embodiment. FIG. 8 is a diagram for explaining a third method of determining the amount of deviation from the reference plane at the element point by the data conversion unit according to the first embodiment.
データ変換部22は、形状データDfに基づいて、図5に示すトンネル5の内壁面5aの全体形状を特定する。例えば、データ変換部22は、内壁面5aの断面形状をトンネル軸方向に連続させることで内壁面5aの全体形状を特定する。また、データ変換部22は、形状データDfが内壁面5aの全体形状を示す3次元データである場合、形状データDfに含まれる3次元データから内壁面5aの全体形状を特定することができる。
The
そして、データ変換部22は、特定した内壁面5aの全体形状を基準面Swr’に設定する。データ変換部22は、3次元点群データDpgに含まれる各3次元点を基準面Swr’に直交する方向から基準面Swr’へ投影する。データ変換部22は、基準面Swr’に投影された各3次元点の基準面Swr’上の位置と、投影前の各3次元点の基準面Swr’からのずれ量とを算出する。基準面Swr’からのずれ量の大きさは3次元点における基準面Swr’からの距離を表し、基準面Swr’からのずれ量の極性は基準面Swr’と直交する方向であって基準面Swr’から見た方向を表す。
Then, the
次に、データ変換部22は、各3次元点における投影後の基準面Swr’上の位置と投影前の基準面Swr’からのずれ量に基づいて、3次元点群データDpgを展開図座標系のデータである展開図点群データDpg1へ変換する。展開図座標系は、トンネル軸方向における位置とトンネル周方向における位置と基準面Swrからのずれ量とを各々異なる軸とする直交座標系である。展開図点群データDpg1は、トンネル軸方向における位置と、トンネル周方向における位置と、基準面Swrからのずれ量とで各々表される複数の3次元点のデータを含む。
Next, the
データ変換部22は、3次元点の投影後の基準面Swr’上の位置を3次元点の基準面Swr上の位置へ変換し、投影前の基準面Swr’からのずれ量を基準面Swrからのずれ量とすることで、3次元点群データDpgを展開図点群データDpg1へ変換する。基準面Swrからのずれ量の大きさは3次元点における基準面Swrからの距離を表し、基準面Swrからのずれ量の極性は基準面Swrと直交する方向であって基準面Swrから見た方向を表す。
The
次に、データ変換部22は、展開図点群データDpg1に含まれる複数の3次元点のうちトンネル軸方向における位置とトンネル周方向における位置とが1つの格子点Pgの位置から予め定められた範囲内にある1つ以上の3次元点を1つの格子点Pgに割り当てる。データ変換部22は、割り当てた1つの格子点Pgのトンネル軸方向における位置とトンネル周方向における位置とを、要素点Peのトンネル軸方向における位置とトンネル周方向における位置に設定する。
Next, in the
また、データ変換部22は、1つの格子点Pgに割り当てられた1つ以上の3次元点の基準面Swrからのずれ量から、要素点Peの基準面Swrからのずれ量を決定する。例えば、図6に示すように、データ変換部22は、格子点Pgに最も近い3次元点の基準面Swrからのずれ量を、要素点Peの基準面Swrからのずれ量に決定することができる。格子点Pgに最も近い3次元点は、格子点Pgとの直線距離が最も短い3次元点である。
Further, the
また、図7に示すように、データ変換部22は、格子点Pgから予め定められた範囲内にある3次元点の基準面Swrからのずれ量を平均して得られる値を、要素点Peの基準面Swrからのずれ量に決定することができる。図7に示す例では、格子点Pgから予め定められた範囲内にある格子点Pgにハッチングが付されている。図7に示すように、格子点Pgから予め定められた範囲内には7つの3次元点が存在するため、これらの7つの3次元点の基準面Swrからのずれ量の平均値が、要素点Peの基準面Swrからのずれ量に決定される。
Further, as shown in FIG. 7, the
また、図8に示すように、データ変換部22は、格子点Pgから予め定められた範囲内にある複数の3次元点を基準面Swrに直交する方向に並べた場合に、並び順の真ん中に位置する3次元点の値を、要素点Peの基準面Swrからのずれ量に決定することができる。図8に示す例では、格子点Pgから予め定められた範囲内にある格子点Pgにハッチングが付されている。図8に示すように、格子点Pgから予め定められた範囲内には7つの3次元点が存在し、7つの3次元点を基準面Swrに直交する方向に並べた場合に、並び順の4番目の位置になる3次元点の基準面Swrからのずれ量が、要素点Peの基準面Swrからのずれ量に決定される。
Further, as shown in FIG. 8, when a plurality of three-dimensional points within a predetermined range from the grid point Pg are arranged in a direction orthogonal to the reference plane Swr, the
データ変換部22は、このように1つ以上の3次元点を1つの格子点Pgに割り当てて要素点Peを算出する処理を、格子点Pg毎に行うことで、展開図点群データDpg1を要素点群データDeへ変換する。なお、データ変換部22は、例えば、3次元点群データDpgから展開図点群データDpg1への変換を行わずに、3次元点群データDpgから要素点群データDeへ直接変換する構成であってもよい。
The
データ変換部22は、基準面Swr’上に複数の格子点Pgを配置することもできる。この場合、データ変換部22は、基準面Swr’上に1つ以上の3次元点を1つの格子点Pgに割り当てる処理を格子点Pg毎に行った後、展開図座標系のデータへ変換することで、要素点群データDeを生成する。
The
図4に戻って、処理部12の説明を続ける。処理部12の段差検出部23は、データ変換部22によって変換された要素点群データDeに対してガウシアン差分フィルタをかける処理であるガウシアン差分処理を行う。図9は、実施の形態1にかかるガウシアン差分処理を説明するための図である。なお、図9においてステップ関数を示す図以外の図は、y軸を省略している。
Returning to FIG. 4, the description of the
ガウシアン差分フィルタは、下記式(1)に示すガウス関数を用いた2つの異なる第1ガウシアンフィルタおよび第2ガウシアンフィルタの各々の出力間の差分を出力するフィルタである。下記式(1)において、「x」は、xy座標系におけるx軸の座標であり、「y」は、xy座標系におけるy軸の座標であり、「σ」は、パラメータである。 The Gaussian difference filter is a filter that outputs the difference between the outputs of two different first Gaussian filters and the second Gaussian filter using the Gaussian function shown in the following equation (1). In the following equation (1), "x" is the coordinate of the x-axis in the xy coordinate system, "y" is the coordinate of the y-axis in the xy coordinate system, and "σ" is the parameter.
ここで、トンネル5の内壁面5aに生じる段差が下記式(2)に示すステップ関数s(x,y)で表されるとする。下記式(2)において、「a」は、ステップの大きさ、すなわち、段差の高さを示す。下記式(2)に示すステップ関数s(x,y)は、図9のように表される。また、下記式(2)において、x軸の座標が零である位置をエッジの位置としており、x軸の座標が零である位置が段差の位置である。x座標は、要素点Peにおけるトンネル軸方向の位置およびトンネル周方向の位置のうち一方の位置を示し、y座標は、要素点Peにおけるトンネル軸方向の位置およびトンネル周方向の位置のうち他方の位置を示す。以下においては、x座標がトンネル軸方向の位置であり、y座標がトンネル周方向の位置であるものとして説明する。
Here, it is assumed that the step generated on the
第1ガウシアンフィルタは、下記式(3)のように表される。下記式(3)において、「Sσ1(x,y)」は、xy座標系における任意の位置における第1ガウシアンフィルタの出力を示す。The first Gaussian filter is represented by the following equation (3). In the following equation (3), "S σ1 (x, y)" indicates the output of the first Gaussian filter at an arbitrary position in the xy coordinate system.
上記式(3)において、「*」はコンボリューションであり、「σ1」は第1ガウシアンフィルタのパラメータであり、「u」は、注目点になる要素点Pe以外の要素点Peにおけるトンネル軸方向の位置である。また、上記式(3)において、「v」は、注目点になる要素点Pe以外の要素点Peにおけるトンネル周方向の位置であり、「a」は、上記式(2)に示すステップ関数s(x,y)のステップの大きさである。In the above equation (3), "*" is a convolution, "σ 1 " is a parameter of the first Gaussian filter, and "u" is a tunnel axis at an element point Pe other than the element point Pe that is the point of interest. The position of the direction. Further, in the above equation (3), "v" is a position in the tunnel circumferential direction at an element point Pe other than the element point Pe which is the point of interest, and "a" is the step function s shown in the above equation (2). The size of the step (x, y).
第2ガウシアンフィルタは、下記式(4)のように表される。下記式(4)において「Sσ2(x,y)」は、xy座標系における任意の位置における第2ガウシアンフィルタの出力を示し、「σ2」は第2ガウシアンフィルタのパラメータである。第1ガウシアンフィルタと第2ガウシアンフィルタとはパラメータの値が互いに異なるため、図9に示すように、第1ガウシアンフィルタの特性と第2ガウシアンフィルタの特性とは互いに異なる。図9に示す例では、σ2>σ1である。また、図9において第1ガウシアンフィルタの特性を示す図は、横軸がx軸の座標を示し、縦軸が第1ガウシアンフィルタの出力を示す。また、図9において第2ガウシアンフィルタの特性を示す図は、横軸がx軸の座標を示し、縦軸が第2ガウシアンフィルタの出力を示す。The second Gaussian filter is expressed by the following equation (4). In the following equation (4), "S σ 2 (x, y)" indicates the output of the second Gaussian filter at an arbitrary position in the xy coordinate system, and "σ 2 " is a parameter of the second Gaussian filter. Since the parameter values of the first Gaussian filter and the second Gaussian filter are different from each other, the characteristics of the first Gaussian filter and the characteristics of the second Gaussian filter are different from each other as shown in FIG. In the example shown in FIG. 9, σ 2 > σ 1 . Further, in the figure showing the characteristics of the first Gaussian filter in FIG. 9, the horizontal axis shows the coordinates of the x-axis and the vertical axis shows the output of the first Gaussian filter. Further, in the figure showing the characteristics of the second Gaussian filter in FIG. 9, the horizontal axis shows the coordinates of the x-axis and the vertical axis shows the output of the second Gaussian filter.
そして、ガウシアン差分フィルタは、下記式(5)のように表される。下記式(5)において、「H(x,y)」は、ガウシアン差分フィルタの出力を示す。下記式(5)および図9に示すように、ガウシアン差分フィルタの特性は、ステップ関数s(x,y)の演算結果が「a」に変化する位置、すなわちx=0である場合に、零交差する。段差位置は、上述したようにx=0の位置であるため、段差検出部23は、零交差が生じる位置を段差の位置として検出する。図9においてガウシアン差分フィルタの特性を示す図は、横軸がx軸の座標を示し、縦軸がガウシアン差分フィルタの出力を示す。
Then, the Gaussian difference filter is expressed by the following equation (5). In the following equation (5), "H (x, y)" indicates the output of the Gaussian difference filter. As shown in the following equation (5) and FIG. 9, the characteristic of the Gaussian difference filter is zero when the operation result of the step function s (x, y) changes to "a", that is, x = 0. Cross. Since the step position is the position of x = 0 as described above, the
図4に示すように、段差検出部23は、段差位置検出部31と、段差高検出部32と、ベクトル生成部33とを備える。段差検出部23の段差位置検出部31は、要素点群データDeに対して上記式(5)で示されるガウシアン差分フィルタの演算を行う。具体的には、段差位置検出部31は、要素点群データDeに含まれる各要素点Peに対して、上記式(5)の演算を行う。そして、段差位置検出部31は、ガウシアン差分フィルタの出力に基づいて、零交差が生じる位置を段差の位置として決定する。
As shown in FIG. 4, the
また、図9に示されるガウシアン差分フィルタの特性を示す曲線において零交差が生じる位置での傾き∇H(x,y)は、下記式(6)の演算によって求めることができる。下記式(6)において、零交差が生じる位置をx=0とした場合、下記式(7)のように表される。 Further, the slope ∇H (x, y) at the position where the zero intersection occurs in the curve showing the characteristics of the Gaussian difference filter shown in FIG. 9 can be obtained by the calculation of the following equation (6). In the following equation (6), when the position where the zero intersection occurs is x = 0, it is expressed as the following equation (7).
そして、∇H(0,y)の絶対値は、下記式(8)のように表すことができる。下記式(8)へ∇H(0,y)の絶対値、「σ1」の値、および「σ2」の値を代入することによって、段差の高さである「a」を算出することができる。Then, the absolute value of ∇H (0, y) can be expressed by the following equation (8). By substituting the absolute value of ∇H (0, y), the value of "σ 1 ", and the value of "σ 2 " into the following formula (8), the height of the step "a" is calculated. Can be done.
図10は、実施の形態1にかかる段差位置検出部による段差の位置を検出する方法の一例を示す図である。図10では、ガウシアン差分フィルタをかけた後の4つの要素点Peの位置が黒丸で示されている。図10に示す要素点Peの位置は、トンネル軸方向の位置とトンネル周方向の位置と基準面Swrからのずれ量の差分とを各々異なる軸とする直交座標系の位置である。基準面Swrからのずれ量の差分は、上述したように、要素点Peに対して第1ガウシアンフィルタをかけて得られる基準面Swrからのずれ量と、要素点Peに対して第2ガウシアンフィルタをかけて得られる基準面Swrからのずれ量との差分である。 FIG. 10 is a diagram showing an example of a method of detecting the position of a step by the step position detecting unit according to the first embodiment. In FIG. 10, the positions of the four element points Pe after the Gaussian difference filter is applied are indicated by black circles. The position of the element point Pe shown in FIG. 10 is a position in a Cartesian coordinate system in which the position in the tunnel axial direction, the position in the tunnel circumferential direction, and the difference in the amount of deviation from the reference plane Swr are different axes. As described above, the difference in the amount of deviation from the reference surface Swr is the amount of deviation from the reference surface Swr obtained by applying the first Gaussian filter to the element point Pe and the second Gaussian filter with respect to the element point Pe. It is a difference from the amount of deviation from the reference plane Swr obtained by multiplying by.
図10に示す「H(xn,yn)」、「H(xn+1,yn)」、「H(xn,yn+1)」、および「H(xn+1,yn+1)」は、4つの要素点Peに対してガウシアン差分フィルタをかけて得られる基準面Swrからのずれ量の差分を示す。nは整数である。具体的には、「H(xn,yn)」は、トンネル軸方向の位置が「xn」であり且つトンネル周方向の位置が「yn」である要素点Pe(xn,yn)に対してガウシアン差分フィルタをかけて得られる基準面Swrからのずれ量の差分である。「H(xn+1,yn)」は、トンネル軸方向の位置が「xn+1」であり且つトンネル周方向の位置が「yn」である要素点Pe(xn+1,yn)に対してガウシアン差分フィルタをかけて得られる基準面Swrからのずれ量の差分である。 “H (x n , y n )”, “H (x n + 1 , y n )”, “H (x n , y n + 1 )”, and “H (x n + 1 , y n + 1 )” shown in FIG. 10 are The difference in the amount of deviation from the reference plane Swr obtained by applying a Gaussian difference filter to the four element points Pe is shown. n is an integer. Specifically, "H (x n , y n )" is an element point Pe (x n , y) whose position in the tunnel axial direction is "x n " and whose position in the tunnel circumferential direction is "y n". It is the difference of the amount of deviation from the reference plane Swr obtained by applying the Gaussian difference filter to n). “H (x n + 1 , y n )” is relative to the element point Pe (x n + 1 , y n ) whose position in the tunnel axial direction is “x n + 1 ” and whose position in the tunnel circumferential direction is “y n”. It is the difference of the amount of deviation from the reference plane Swr obtained by applying the Gaussian difference filter.
「H(xn,yn+1)」は、トンネル軸方向の位置が「xn」であり且つトンネル周方向の位置が「yn+1」である要素点Pe(xn,yn+1)に対してガウシアン差分フィルタをかけて得られる基準面Swrからのずれ量の差分である。「H(xn+1,yn+1)」は、トンネル軸方向の位置が「xn+1」であり且つトンネル周方向の位置が「yn+1」である要素点Pe(xn+1,yn+1)に対してガウシアン差分フィルタをかけて得られる基準面Swrからのずれ量の差分である。「xn」と「xn+1」とは、トンネル軸方向の位置が隣接する2つ要素点Peの位置を示し、「yn」と「yn+1」とは、トンネル周方向の位置が隣接する2つ要素点Peの位置を示す。“H (x n , y n + 1 )” is relative to the element point Pe (x n , y n + 1 ) whose position in the tunnel axial direction is “x n ” and whose position in the tunnel circumferential direction is “y n + 1”. It is the difference of the amount of deviation from the reference plane Swr obtained by applying the Gaussian difference filter. “H (x n + 1 , y n + 1 )” is for the element point Pe (x n + 1 , y n + 1 ) whose position in the tunnel axis direction is “x n + 1 ” and whose position in the tunnel circumferential direction is “y n + 1”. It is the difference of the amount of deviation from the reference plane Swr obtained by applying the Gaussian difference filter. “X n ” and “x n + 1 ” indicate the positions of two element points Pe whose positions in the tunnel axial direction are adjacent to each other, and “y n ” and “y n + 1 ” indicate positions in the tunnel circumferential direction are adjacent to each other. The positions of the two element points Pe are shown.
段差位置検出部31は、ガウシアン差分フィルタをかけて得られる複数の要素点Peのうちトンネル軸方向の位置が隣接する2つの要素点Pe同士を結ぶ線が零交差する位置である零交差点を判定する。例えば、段差位置検出部31は、トンネル軸方向の位置が隣接する要素点Pe(xn,yn)と要素点Pe(xn+1,yn)とで基準面Swrからのずれ量の差分の極性が異なるか否かを判定する。すなわち、段差位置検出部31は、「H(xn,yn)」と「H(xn+1,yn)」とで極性が異なるか否かを判定する。段差位置検出部31は、要素点Pe(xn,yn)における「H(xn,yn)」と要素点Pe(xn+1,yn)における「H(xn+1,yn)」とで極性が異なると判定した場合、要素点Pe(xn,yn)と要素点Pe(xn+1,yn)を内分して零交差点を求める。例えば、下記式(9)の演算によって、トンネル軸方向において零交差が生じる位置である零交差点の位置xcを算出することができる。The step position detection unit 31 determines the zero intersection, which is the position where the line connecting the two element points Pe having adjacent positions in the tunnel axis direction intersects zero among the plurality of element points Pe obtained by applying the Gaussian difference filter. do. For example, the step position detection unit 31 determines the difference in the amount of deviation from the reference plane Swr between the element points Pe (x n , y n ) and the element points Pe (x n + 1 , y n ) whose positions in the tunnel axial direction are adjacent to each other. Determine if the polarities are different. That is, the step position detection unit 31 determines whether or not the polarities of "H (x n , y n )" and "H (x n + 1 , y n)" are different. It stepped position detection unit 31, element point Pe (x n, y n) "H (x n, y n)" in the element point Pe (x n + 1, y n) in the "H (x n + 1, y n) " When it is determined that the polarities are different from each other, the element point Pe (x n , y n ) and the element point Pe (x n + 1 , y n ) are internally divided to obtain the zero intersection. For example, the position x c of the zero intersection, which is the position where the zero intersection occurs in the tunnel axis direction, can be calculated by the calculation of the following equation (9).
また、段差位置検出部31は、ガウシアン差分フィルタをかけて得られる複数の要素点Peのうちトンネル周方向の位置が隣接する2つの要素点Pe同士を結ぶ線が零交差する位置である零交差点を判定する。例えば、段差位置検出部31は、トンネル周方向の位置が隣接する要素点Pe(xn,yn)と要素点Pe(xn,yn+1)とで基準面Swrからのずれ量の差分の極性が異なるか否かを判定する。すなわち、段差位置検出部31は、「H(xn,yn)」と「H(xn,yn+1)」とで極性が異なるか否かを判定する。段差位置検出部31は、要素点Pe(xn,yn)における「H(xn,yn)」と要素点Pe(xn,yn+1)における「H(xn,yn+1)」とで極性が異なると判定した場合、要素点Pe(xn,yn)と要素点Pe(xn,yn+1)を内分して零交差点を求める。例えば、下記式(10)の演算によって、トンネル周方向において零交差が生じる位置である零交差点の位置ycを算出することができる。Further, the step position detection unit 31 is a zero intersection where the lines connecting the two element points Pe that are adjacent to each other in the tunnel circumferential direction among the plurality of element points Pe obtained by applying the Gaussian difference filter intersect at zero. To judge. For example, the step position detection unit 31 determines the difference in the amount of deviation from the reference plane Swr between the element points Pe (x n , y n ) and the element points Pe (x n , y n + 1 ) whose positions in the tunnel circumferential direction are adjacent to each other. Determine if the polarities are different. That is, the step position detection unit 31 determines whether or not the polarities of "H (x n , y n )" and "H (x n , y n + 1)" are different. It stepped position detection unit 31, element point Pe (x n, y n) "H (x n, y n)" in the element point Pe (x n, y n + 1) in the "H (x n, y n + 1) " When it is determined that the polarities are different, the element point Pe (x n , y n ) and the element point Pe (x n , y n + 1 ) are internally divided to obtain the zero intersection. For example, the position y c of the zero intersection, which is the position where the zero intersection occurs in the tunnel circumferential direction, can be calculated by the calculation of the following equation (10).
なお、段差位置検出部31は、基準面Swrからのずれ量の差分の極性が異なる2つの要素点Peのうち基準面Swrからのずれ量の差分の絶対値が小さい要素点Peの位置を段差の位置として検出することができる。このように、段差位置検出部31は、零交差点に近い位置を段差の位置として検出することもできる。 The step position detection unit 31 sets the position of the element point Pe having the smaller absolute value of the difference in the amount of deviation from the reference surface Swr among the two element points Pe having different polarities of the difference in the amount of deviation from the reference surface Swr. Can be detected as the position of. In this way, the step position detection unit 31 can also detect a position close to the zero intersection as the step position.
次に、図4に示す段差高検出部32について説明する。段差高検出部32は、内壁面5aに生じる段差の高さを検出する。図11は、実施の形態1にかかる段差高検出部による段差の高さを検出する方法の一例を示す図である。なお、図11では、図10と同様にガウシアン差分フィルタをかけた後の4つの要素点Peの位置が黒丸で示されている。
Next, the step height detection unit 32 shown in FIG. 4 will be described. The step height detecting unit 32 detects the height of the step generated on the
段差高検出部32は、ガウシアン差分フィルタの出力H(x,y)における零交差点での傾き∇H(x,y)を求める。段差高検出部32は、零交差点での傾き∇H(x,y)に基づいて、段差の高さを算出する。 The step height detection unit 32 obtains the slope ∇H (x, y) at the zero intersection at the output H (x, y) of the Gaussian difference filter. The step height detection unit 32 calculates the step height based on the slope ∇H (x, y) at the zero intersection.
図11に示す例では、段差高検出部32は、零交差点P(xc,yn)および零交差点P(xn+1,yc)での傾き∇H(x,y)を、零交差点P(xc,yn)と零交差点P(xn+1,yc)とを結んで得られる後述の段差ベクトル毎に算出する。例えば、段差高検出部32は、零交差点P(xc,yn)および零交差点P(xn+1,yc)での傾き∇H(x,y)として、要素点Pe(xn,yn)、要素点Pe(xn+1,yn)、および要素点Pe(xn+1,yn+1)とを含む面の零平面に対する傾きを求める。零平面は、ガウシアン差分フィルタの出力H(x,y)が零になる面である。In the example shown in FIG. 11, the step height detection unit 32 sets the inclination ∇H (x, y) at the zero intersection P (x c , y n ) and the zero intersection P (x n + 1 , y c) to the zero intersection P. It is calculated for each step vector described later obtained by connecting (x c , y n ) and the zero intersection P (x n + 1 , y c). For example, the step height detection unit 32 has an element point Pe (x n , y) as an inclination ∇H (x, y) at the zero intersection P (x c , y n ) and the zero intersection P (x n + 1 , y c). Find the slope of the plane containing n ), the element point Pe (x n + 1 , y n ), and the element point Pe (x n + 1 , y n + 1 ) with respect to the zero plane. The zero plane is the plane where the output H (x, y) of the Gaussian difference filter becomes zero.
また、段差高検出部32は、零交差点P(xc,yn)および零交差点P(xn+1,yc)での傾き∇H(x,y)として、要素点Pe(xn,yn)、要素点Pe(xn+1,yn)、要素点Pe(xn+1,yn+1)、および要素点Pe(xn+1,yn+1)とを含む面の零平面に対する傾きを求めることもできる。Further, the step height detection unit 32 has an element point Pe (x n , y) as an inclination ∇H (x, y) at the zero intersection P (x c , y n ) and the zero intersection P (x n + 1 , y c). It is also possible to obtain the slope of the surface including n ), element point Pe (x n + 1 , y n ), element point Pe (x n + 1 , y n + 1 ), and element point Pe (x n + 1 , y n + 1) with respect to the zero plane. ..
このように、段差高検出部32は、ガウシアン差分フィルタをかけた後の複数の要素点Peのうち隣接する3以上の要素点Peであって少なくとも1つの要素点Peと他の且つ残りの要素点Peとでガウシアン差分フィルタをかけて得られる値の極性が異なる3以上の要素点Peを含む面の傾きを傾き∇H(x,y)として、段差の高さを検出する。 As described above, the step height detection unit 32 is three or more adjacent element points Pe among the plurality of element points Pe after the Gaussian difference filter is applied, and is at least one element point Pe and other and remaining elements. The height of the step is detected by setting the slope of the surface including the three or more element points Pe having different polarities of the values obtained by applying the Gaussian difference filter with the point Pe as the slope ∇H (x, y).
また、段差高検出部32は、ガウシアン差分フィルタをかけて得られる値であるガウシアン差分フィルタの出力H(x,y)の極性が異なる2つの要素点Pe同士を結ぶ線の零平面に対する傾きを傾き∇H(x,y)として段差の高さを検出することができる。図11に示す例では、段差高検出部32は、零交差点P(xc,yn)での傾き∇H(x,y)として、要素点Pe(xn,yn)と要素点Pe(xn+1,yn)とを結ぶ線の零平面に対する傾きを傾き∇H(x,y)として求めることができる。Further, the step height detection unit 32 determines the slope of the line connecting the two element points Pe having different polarities of the output H (x, y) of the Gaussian difference filter, which is a value obtained by applying the Gaussian difference filter, with respect to the zero plane. The height of the step can be detected as the slope ∇H (x, y). In the example shown in FIG. 11, the step height detection unit 32 has an element point Pe (x n , y n ) and an element point Pe as an inclination ∇ H (x, y) at the zero intersection P (x c , y n). The slope of the line connecting (x n + 1 , y n ) with respect to the zero plane can be obtained as the slope ∇H (x, y).
また、図11に示す例では、段差高検出部32は、零交差点P(xn+1,yc)での傾き∇H(x,y)として、要素点Pe(xn+1,yn)と要素点Pe(xn+1,yn+1)とを結ぶ線の零平面に対する傾きを傾き∇H(x,y)として求めることができる。また、段差高検出部32は、零交差点P(xc,yn+1)および零交差点P(xn+1,yc)での傾き∇H(x,y)として、要素点Pe(xn,yn)と要素点Pe(xn+1,yn)とを結ぶ線の零平面に対する傾きと要素点Pe(xn+1,yn)と要素点Pe(xn+1,yn+1)とを結ぶ線の零平面に対する傾きとの平均値を傾き∇H(x,y)として求めることもできる。Further, in the example shown in FIG. 11, the step height detection unit 32 has an element point Pe (x n + 1 , y n ) and an element as an inclination ∇ H (x, y) at the zero intersection P (x n + 1 , y c). The slope of the line connecting the points Pe (x n + 1 , y n + 1 ) with respect to the zero plane can be obtained as the slope ∇H (x, y). Further, the step height detection unit 32 has an element point Pe (x n , y) as an inclination ∇H (x, y) at the zero intersection P (x c , y n + 1 ) and the zero intersection P (x n + 1 , y c). The slope of the line connecting n) and the element point Pe (x n + 1 , y n ) with respect to the zero plane, and the zero of the line connecting the element point Pe (x n + 1 , y n ) and the element point Pe (x n + 1 , y n + 1 ). The average value of the slope with respect to the plane can also be obtained as the slope ∇H (x, y).
そして、段差高検出部32は、下記式(11)の演算によって、段差の高さ「a」を検出する。下記式(11)は、上記式(8)に基づいて求めることができる。 Then, the step height detection unit 32 detects the step height “a” by the calculation of the following equation (11). The following formula (11) can be obtained based on the above formula (8).
次に、ベクトル生成部33について説明する。ベクトル生成部33は、内壁面5aに生じる段差の高さを検出する。図12は、実施の形態1にかかるベクトル生成部による段差ベクトルを生成する方法の一例を示す図である。なお、図12では、図10と同様にガウシアン差分フィルタをかけた後の4つの要素点Peの位置が黒丸で示されている。
Next, the vector generation unit 33 will be described. The vector generation unit 33 detects the height of the step generated on the
ベクトル生成部33は、隣接する4つの要素点Peを直線で結んで構成されるメッシュ内で、2以上の零交差点がある場合、零交差点を結ぶ線を段差ベクトルの単位ベクトルとして生成する。要素点Peがトンネル軸方向とトンネル周方向とに各々10ずつ配列されている場合、隣接する4つの要素点Peを直線で結んで構成されるメッシュの数は、81個である。ベクトル生成部33は、各メッシュに対して、2以上の零交差点がある場合に零交差点を結ぶ線を単位ベクトルとして生成する処理を行う。図12に示す例では、ベクトル生成部33は、零交差点P(xc,yn)と零交差点P(xn+1,yc)とを結ぶ線を単位ベクトルとして生成する。When there are two or more zero intersections in the mesh formed by connecting four adjacent element points Pe with a straight line, the vector generation unit 33 generates a line connecting the zero intersections as a unit vector of the step vector. When 10 element points Pe are arranged in each of the tunnel axial direction and the tunnel circumferential direction, the number of meshes formed by connecting four adjacent element points Pe with a straight line is 81. The vector generation unit 33 performs a process of generating a line connecting the zero intersections as a unit vector when there are two or more zero intersections for each mesh. In the example shown in FIG. 12, the vector generation unit 33 generates a line connecting the zero intersection P (x c , y n ) and the zero intersection P (x n + 1 , y c ) as a unit vector.
ベクトル生成部33は、単位ベクトルを生成したメッシュに隣接するメッシュである隣接メッシュのうち零交差点を共有する隣接メッシュの単位ベクトルを、連続する単位ベクトルとして接続する。これにより、ベクトル生成部33は、連続する段差に対応する段差ベクトルを生成することができる。なお、ベクトル生成部33は、いずれの隣接メッシュにも単位ベクトルがない場合、一つの単位ベクトルを段差ベクトルとして扱う。 The vector generation unit 33 connects the unit vectors of the adjacent meshes that share the zero intersection among the adjacent meshes that are adjacent to the mesh that generated the unit vector as continuous unit vectors. As a result, the vector generation unit 33 can generate a step vector corresponding to a continuous step. If there is no unit vector in any of the adjacent meshes, the vector generation unit 33 treats one unit vector as a step vector.
ベクトル生成部33は、段差ベクトルの節を減らすことで段差ベクトルを簡略化することもできる。例えば、ベクトル生成部33は、生成した段差ベクトルをDouglas−Peucker法などを用いて、段差ベクトルの節を減らすことができる。なお、段差ベクトルの節は、段差ベクトルを構成する連続する2つの単位ベクトルが為す角が0度以外の角度である場合に零交差点の位置に形成される。 The vector generation unit 33 can also simplify the step vector by reducing the nodes of the step vector. For example, the vector generation unit 33 can reduce the number of steps of the step vector by using the Douglas-Pucker method or the like for the generated step vector. The node of the step vector is formed at the position of the zero intersection when the angle formed by the two consecutive unit vectors constituting the step vector is an angle other than 0 degrees.
段差検出部23は、検出した段差の位置を示す情報、検出した段差の高さを示す情報、および生成した段差ベクトルを示す情報を含む段差データを、記憶部11に記憶される変状データDdfに加えることによって、変状データDdfを更新する。また、段差検出部23は、段差の検出に用いられた3次元点群データDpgに含まれる時刻情報を段差データに追加することができる。
The
図4に戻って、処理部12の説明を続ける。処理部12の段差評価部24は、段差検出部23によって検出された段差の高さおよび段差の広がりに基づいて、段差の状態をレベルで表す状態レベルを判定する。段差評価部24は、判定した段差レベルを示す情報を段差データに追加することができる。
Returning to FIG. 4, the description of the
段差評価部24は、基準面Swr上での広がりを段差の広がりとして算出する。段差ベクトルの広がりは、例えば、段差ベクトルの長さ、および段差ベクトルが存在する2次元範囲の大きさの少なくとも一方から算出される。段差ベクトルが存在する2次元範囲の大きさは、例えば、段差ベクトルを囲むことができる最も小さい矩形の面積である。 The step evaluation unit 24 calculates the spread on the reference surface Swr as the spread of the step. The spread of the step vector is calculated from, for example, at least one of the length of the step vector and the magnitude of the two-dimensional range in which the step vector exists. The magnitude of the two-dimensional range in which the step vector exists is, for example, the area of the smallest rectangle that can enclose the step vector.
段差評価部24は、段差の高さと段差の広がりとをパラメータとする演算式またはテーブルを用いて、段差の状態レベルを判定することができる。状態レベルは、例えば、段差の危険度または段差の劣化度を示す。段差評価部24は、例えば、段差検出部23によって生成された段差ベクトル毎に、第1レベル、第2レベル、および第3レベルからなる3段階のレベルから1つのレベルを判定することができる。第1レベルは、例えば、段差が危険な状態でないことを示すレベルである。第2レベルは、例えば、段差が危険な状態でないが注意を要する状態であることを示すレベルである。第3レベルは、例えば、段差が危険な状態であることを示すレベルである。
The step evaluation unit 24 can determine the state level of the step by using an arithmetic expression or a table in which the height of the step and the spread of the step are parameters. The state level indicates, for example, the degree of risk of the step or the degree of deterioration of the step. For example, the step evaluation unit 24 can determine one level from three levels including the first level, the second level, and the third level for each step vector generated by the
段差評価部24は、2段階のレベルから段差の状態レベルを判定することもでき、4段階以上のレベルから段差の状態レベルを判定することもできる。段差評価部24は、判定した段差の状態レベルを示す情報を段差データに追加することができる。この場合、段差データには、検出した段差の位置を示す情報、検出した段差の高さを示す情報、および生成した段差ベクトルを示す情報に加え、段差の状態レベルを示す情報が含まれる。 The step evaluation unit 24 can determine the state level of the step from two levels, and can also determine the state level of the step from four or more levels. The step evaluation unit 24 can add information indicating the state level of the determined step to the step data. In this case, the step data includes information indicating the position of the detected step, information indicating the height of the detected step, information indicating the generated step vector, and information indicating the state level of the step.
属性設定部25は、段差評価部24によって算出された段差ベクトルの長さおよび段差の広がりなどの情報を段差の属性情報として記憶部11に記憶されている段差データに追加する。また、属性設定部25は、入力部13によって入力される段差の属性情報を記憶部11に記憶されている段差データを追加する。入力部13によって入力される段差の属性情報は、例えば、識別番号およびコメントなどの段差属性を示す情報である。属性設定部25は、段差データに含まれる段差ベクトルの情報に段差の属性情報を関連付ける。入力部13は、例えば、キーボードなどである。このように、属性設定部25によって段差データに段差の属性情報を含めることができる。
The
座標変換部26は、記憶部11に記憶された変状データDdfの段差データに含まれる段差ベクトルの情報に基づいて、段差ベクトルを上述した展開図座標系の座標から上述した3次元直交座標系の座標へ変換する。座標変換部26は、変換した3次元直交座標系の段差ベクトルの座標を示す情報を3次元段差ベクトルデータとして記憶部11に記憶された段差データに追加することができる。これにより、段差データには、展開図座標系の段差ベクトルの情報と3次元直交座標系の段差ベクトルの情報とが含まれる。 The coordinate conversion unit 26 uses the step vector information included in the step data of the deformation data Ddf stored in the storage unit 11 to obtain the step vector from the coordinates of the above-mentioned developed view coordinate system to the above-mentioned three-dimensional orthogonal coordinate system. Convert to the coordinates of. The coordinate conversion unit 26 can add information indicating the coordinates of the step vector of the converted three-dimensional Cartesian coordinate system to the step data stored in the storage unit 11 as three-dimensional step vector data. As a result, the step data includes the step vector information of the developed view coordinate system and the step vector information of the three-dimensional Cartesian coordinate system.
また、座標変換部26は、記憶部11に記憶された変状データDdfと画像データDimgに基づいて、トンネル5の画像をトンネル5の内壁面5aの形状に沿うように座標変換して、トンネル5の3次元画像を生成することもできる。また、座標変換部26は、記憶部11に記憶された3次元点群データDpgに基づいて、トンネル5の内壁面5aの3次元点群を3次元座標系から展開図座標系へ変換し、変換した3次元点群を展開図平面へ投影し、投影した複数の点を2次元点群とする2次元点群データを記憶部11に記憶させることができる。
Further, the coordinate conversion unit 26 performs coordinate conversion of the image of the
表示処理部27は、記憶部11に記憶された形状データDfおよび変状データDdfに基づき、基準面Swrで表される展開図に対して、段差ベクトルを示す図形と段差属性を示す情報とを付加した図である変状展開図を示す2次元画像を表示部14に表示させることができる。表示処理部27は、記憶部11に記憶された形状データDf、変状データDdf、および画像データDimgに基づいて、トンネル5の内壁面5aの撮像画像に変状展開図を示す画像を重畳した2次元画像を表示部14に表示させることができる。また、表示処理部27は、記憶部11に記憶した形状データDf、変状データDdf、および2次元点群データに基づいて、変状展開図に2次元点群を重畳した2次元画像を表示部14に表示させることもできる。
Based on the shape data Df and the deformation data Ddf stored in the storage unit 11, the display processing unit 27 displays a figure showing the step vector and information showing the step attribute with respect to the developed view represented by the reference plane Swr. A two-dimensional image showing a modified development view, which is an added figure, can be displayed on the display unit 14. The display processing unit 27 superimposes an image showing a deformation development diagram on the captured image of the
図13は、実施の形態1にかかる変状展開図の一例を示す図である。図13に示す変状展開図80は、トンネル5の内壁面5aの外形を展開した外形枠81と、段差ベクトル82a,82b,82c,82dと、属性情報83a,83b,83c,83dを含む。段差ベクトル82a,82b,82c,82dは、段差ベクトルを示す2次元の図形である。
FIG. 13 is a diagram showing an example of a modified development diagram according to the first embodiment. The modified
属性情報83a,83b,83c,83dは、段差ベクトルの長さを示す情報である。属性情報83a,83b,83c,83dは、段差ベクトルの長さに代えてまたは加えて、入力部13から入力された識別番号およびコメントを含んでいてもよい。
The
また、表示処理部27は、記憶部11に記憶された3次元点群データDpgと3次元段差ベクトルデータとに基づいて、トンネル5の内壁面5aを表す3次元点群に、3次元の段差ベクトルを重畳した3次元画像を表示部14に表示させることができる。また、表示処理部27は、3次元点群データDpgと3次元段差ベクトルデータと画像データDimgに基づいて、トンネル5の内壁面5aを表す3次元点群に、トンネル5の3次元画像と3次元の段差ベクトルとを重畳した3次元画像を表示部14に表示させることもできる。
Further, the display processing unit 27 adds a three-dimensional step to the three-dimensional point cloud representing the
図14は、実施の形態1にかかるトンネルの内壁面を表す3次元点群に3次元の段差ベクトルを重畳した3次元画像の一例を示す図である。図14に示す3次元画像90は、トンネル5の内壁面5aを表す3次元点群91と、段差ベクトル92a,92b,92c,92dとを含む。段差ベクトル92a,92b,92c,92dは、段差ベクトルを示す3次元の図形である。なお、3次元画像90は、段差ベクトルの属性情報を含むこともできる。
FIG. 14 is a diagram showing an example of a three-dimensional image in which a three-dimensional step vector is superimposed on a three-dimensional point cloud representing the inner wall surface of the tunnel according to the first embodiment. The three-
また、表示処理部27は、段差検出部23によって検出された段差の情報を時刻情報に基づいて時系列に表示部14に表示させることができる。例えば、表示処理部27は、段差データに含まれる時刻情報に基づいて、各時刻で検出された段差ベクトルを検出時刻の順に切り替えながら、表示部14に表示させることができる。これにより、変状検出システム100の利用者は、段差の状態変化を把握することができる。例えば、表示処理部27は、各時刻で検出された段差ベクトル82a,82b,82c,82dを検出時刻の順に切り替えながら表示部14に表示させることができる。また、表示処理部27は、各時刻で検出された段差ベクトル92a,92b,92c,92dを検出時刻の順に切り替えながら表示部14に表示させることができる。
Further, the display processing unit 27 can display the step information detected by the
図15は、実施の形態1にかかる段差検出装置の処理部による処理手順の一例を示すフローチャートである。図15に示す処理は、処理部12によって繰り返し実行される。図15に示すように、段差検出装置2の処理部12は、段差検出タイミングになったか否かを判定する(ステップS10)。例えば、処理部12は、走行型計測装置1から3次元点群データDpgを取得した場合、または入力部13への特定操作があった場合に、段差検出タイミングになったと判定する。
FIG. 15 is a flowchart showing an example of a processing procedure by the processing unit of the step detection device according to the first embodiment. The process shown in FIG. 15 is repeatedly executed by the
処理部12は、段差検出タイミングになったと判定した場合(ステップS10:Yes)、段差検出処理を行う(ステップS11)。ステップS11の段差検出処理は、図16に示すステップS20からS23の処理であり、後で詳述する。次に、処理部12は、ステップS11で検出した段差の段差ベクトルを3次元直交座標系の段差ベクトルへ変換する座標変換処理を行う(ステップS12)。また、処理部12は、ステップS11で検出した段差の高さと段差の広がりとをパラメータとする演算式またはテーブルを用いて、段差の状態レベルを判定する段差評価処理を行う(ステップS13)。
When the
処理部12は、ステップS13の処理が終了した場合、または段差検出タイミングになっていないと判定した場合(ステップS10:No)、属性入力があるか否かを判定する(ステップS14)。ステップS14において、処理部12は、入力部13へ段差の属性情報が入力された場合に属性入力があると判定する。処理部12は、属性入力があると判定した場合(ステップS14:Yes)、入力された属性情報を記憶部11に記憶させる(ステップS15)。
When the processing of step S13 is completed, or when it is determined that the step detection timing has not been reached (step S10: No), the
処理部12は、ステップS15の処理が終了した場合、または属性入力がないと判定した場合(ステップS14:No)、表示要求があるか否かを判定する(ステップS16)。ステップS16において、処理部12は、例えば、入力部13への特定操作があった場合に、表示要求があると判定する。処理部12は、表示要求があると判定した場合(ステップS16:Yes)、表示要求の対象を表示部14に表示させる(ステップS17)。ステップS17において、処理部12は、表示要求が変状展開図の表示要求である場合、変状展開図を示す2次元画像を表示部14に表示させる。また、処理部12は、表示要求が時系列の表示要求である場合、段差検出部23によって検出された段差の情報を時刻情報に基づいて時系列に表示部14に表示させる。
When the processing of step S15 is completed or when it is determined that there is no attribute input (step S14: No), the
処理部12は、ステップS17の処理が終了した場合、または表示要求がないと判定した場合(ステップS16:No)、図15に示す処理を終了する。
When the process of step S17 is completed or when it is determined that there is no display request (step S16: No), the
図16は、実施の形態1にかかる段差検出装置の処理部による段差検出処理の一例を示すフローチャートである。図16に示すように、段差検出装置2の処理部12は、3次元点群データDpgを要素点群データDeへ変換する(ステップS20)。次に、処理部12は、要素点群データDeをガウシアン差分フィルタにかけるガウシアン差分処理を行う(ステップS21)。処理部12は、ガウシアン差分フィルタの出力に基づいて、段差を検出する(ステップS22)。そして、処理部12は、検出した段差の情報を示す段差データを記憶部11に記憶して(ステップS23)、図16に示す処理を終了する。
FIG. 16 is a flowchart showing an example of the step detection process by the processing unit of the step detection device according to the first embodiment. As shown in FIG. 16, the
図17は、実施の形態1にかかる段差検出装置のハードウェア構成の一例を示す図である。図17に示すように、段差検出装置2は、プロセッサ101と、メモリ102と、通信装置103と、入出力回路104と、表示装置105とを備えるコンピュータを含む。
FIG. 17 is a diagram showing an example of the hardware configuration of the step detection device according to the first embodiment. As shown in FIG. 17, the
プロセッサ101、メモリ102、通信装置103、入出力回路104、および表示装置105は、例えば、バス106によって互いにデータの送受信が可能である。通信部10は、通信装置103で実現される。記憶部11は、メモリ102によって実現される。入力部13は、入出力回路104によって実現される。表示部14は、表示装置105によって実現される。プロセッサ101は、メモリ102に記憶されたプログラムを読み出して実行することによって、取得部21、データ変換部22、段差検出部23、段差評価部24、属性設定部25、座標変換部26、および表示処理部27の機能を実行する。プロセッサ101は、例えば、処理回路の一例であり、CPU(Central Processing Unit)、DSP(Digital Signal Processer)、およびシステムLSI(Large Scale Integration)のうち一つ以上を含む。
The
メモリ102は、RAM(Random Access Memory)、ROM(Read Only Memory)、フラッシュメモリ、EPROM(Erasable Programmable Read Only Memory)、およびEEPROM(登録商標)(Electrically Erasable Programmable Read Only Memory)のうち一つ以上を含む。また、メモリ102は、コンピュータが読み取り可能なプログラムが記録された記録媒体を含む。かかる記録媒体は、不揮発性または揮発性の半導体メモリ、磁気ディスク、フレキシブルメモリ、光ディスク、コンパクトディスク、およびDVD(Digital Versatile Disc)のうち一つ以上を含む。なお、段差検出装置2は、ASIC(Application Specific Integrated Circuit)およびFPGA(Field Programmable Gate Array)などの集積回路を含んでいてもよい。
The
なお、上述した例では、段差検出装置2の段差検出部23によって段差が検出される構造物の表面は、トンネル5の内壁面5aであるが、段差検出装置2の段差検出部23は、トンネル5の内壁面5a以外に生じる段差を検出することができる。例えば、段差検出部23は、橋梁の表面、道路の表面、およびビルの表面を検出対象として、段差を検出することができる。また、段差検出部23は、3次元点群で表される構造物の表面が平面である場合、展開図座標系へ変換することなく、3次元点群データDpgを要素点群データDeへ変換することができる。
In the above example, the surface of the structure in which the step is detected by the
以上のように、実施の形態1にかかる変状検出システム100の段差検出装置2は、取得部21と、データ変換部22と、段差検出部23とを備える。取得部21はトンネル5の内壁面5aを3次元点群で表すデータである3次元点群データDpgとトンネル5の形状を示す形状データDfとを取得する。データ変換部22は、形状データDfに基づいてトンネル5の内壁面5aに対応する基準面Swrを設定し、3次元点群データDpgを、基準面Swrに行列状に配列される複数の点Pgから選択される1つの点Pgの位置と基準面Swrからのずれ量とで各々表される複数の要素点Peのデータを含む要素点群データDeへ変換する。段差検出部23は、要素点群データDeに含まれる複数の要素点Peに対してガウシアン差分フィルタをかけて得られる結果に基づいて、トンネル5の内壁面5aにおける段差の位置を検出する。トンネル5は、構造物の一例であり、トンネル5の内壁面5aは、構造物の表面の一例である。これにより、例えば、トンネル5の内壁面5aを表す3次元点群において段差による3次元点の疎密差が少ない場合であっても、トンネル5の内壁面5aに生じる段差の位置を検出することができる。
As described above, the
また、段差検出部23は、ガウシアン差分フィルタをかけた後の複数の要素点Peのうち互いに隣接し且つガウシアン差分フィルタをかけて得られる値の極性が異なる2つの要素点Pe同士を結ぶ線に零交差が生じる位置を特定し、特定した位置を段差の位置として検出する。これにより、ガウシアン差分フィルタから得られる結果から、トンネル5の内壁面5aに生じる段差の位置を精度よく検出することができる。
Further, the
また、段差検出部23は、ガウシアン差分フィルタをかけた後の複数の要素点Peのうち隣接する3以上の要素点Peであって少なくとも1つの要素点Peと残りの要素点Peとでガウシアン差分フィルタをかけて得られる値の極性が異なる3以上の要素点Peを含む面の傾きに基づいて、段差の高さを検出する段差高検出部32を備える。これにより、ガウシアン差分フィルタから得られる結果から、段差の高さを検出することができる。
Further, the
また、段差検出部23は、ガウシアン差分フィルタをかけて得られる結果から隣接する2つの要素点Pe同士を結ぶ線の零交差点の位置での傾きに基づいて段差の高さを検出する段差高検出部32を備える。これにより、ガウシアン差分フィルタから得られる結果から、段差の高さを精度よく検出することができる。
Further, the
また、段差検出部23は、段差の位置を結んで得られるベクトルである段差ベクトルを生成するベクトル生成部33を備える。これにより、連続する段差を表す段差ベクトルを生成することができる。
Further, the
また、段差検出装置2は、ベクトル生成部33で生成された段差ベクトルに基づいて段差の広がりを算出し、算出した段差の広がりと段差高検出部32によって検出された段差の高さとに基づいて、段差の状態をレベルで表す状態レベルを判定する段差評価部24を備える。これにより、段差の状態を容易に把握することができる。
Further, the
また、段差検出装置2は、ベクトル生成部33によって生成された段差ベクトルを表示部14に表示させる表示処理部27を備える。これにより、段差の状態を容易に把握することができる。
Further, the
また、表示処理部27は、トンネル5の展開図および撮像画像のうち少なくとも一つに段差ベクトル82a,82b,82c,82dを重畳した画像を表示部14に表示させる。これにより、例えば、段差の内壁面5a上の位置および段差の広がりなどを容易に把握することができる。
Further, the display processing unit 27 causes the display unit 14 to display an image in which the
また、段差検出装置2は、段差ベクトル82a,82b,82c,82dを3次元座標系のデータへ変換する座標変換部26を備える。表示処理部27は、3次元座標系のデータに変換された段差ベクトル92a,92b,92c,92dを表示部14に表示させる。これにより、段差の状態を3次元で把握することができる。
Further, the
また、表示処理部27は、トンネル5の内壁面5aを表す3次元点群の画像に3次元座標系のデータに変換された段差ベクトル92a,92b,92c,92dを重畳した画像を表示部14に表示させる。これにより、段差ベクトル92a,92b,92c,92dに対応する画像によりトンネル5の内壁面5aにおける段差の位置および段差の広がりなどを容易に把握することができる。
Further, the display processing unit 27 displays an image in which the
また、取得部21は、異なる時刻で各々トンネル5の内壁面5aを計測して得られた複数の3次元点群データDpgを取得する。複数の3次元点群データDpgは、トンネル5の内壁面5aを計測した時刻を示す時刻情報を含む。段差検出部23は、複数の3次元点群データDpgの各々から段差を検出し、検出した段差に時刻情報を関連付ける。表示処理部27は、段差検出部23によって検出された段差の情報を時刻情報に基づいて時系列に表示部14に表示させる。これにより、段差の状態を時系列に把握することができる。
Further, the
実施の形態2.
実施の形態2にかかる変状検出システムは、段差検出装置が段差の経年変化を検出すると共に段差の分布などを算出する点で、実施の形態1にかかる変状検出システム100と異なる。以下においては、実施の形態1と同様の機能を有する構成要素については同一符号を付して説明を省略し、実施の形態1の変状検出システム100と異なる点を中心に説明する。
The deformation detection system according to the second embodiment is different from the
図18は、本発明の実施の形態2にかかる変状検出システムにおける段差検出装置の構成例を示す図である。図18に示すように、実施の形態2にかかる変状検出システムの段差検出装置2Aは、処理部12に代えて、処理部12Aを備える点で、段差検出装置2と異なる。処理部12Aは、比較部28と、統計処理部29とをさらに備える点で、処理部12と異なる。
FIG. 18 is a diagram showing a configuration example of a step detection device in the deformation detection system according to the second embodiment of the present invention. As shown in FIG. 18, the step detection device 2A of the deformation detection system according to the second embodiment is different from the
比較部28は、走行型計測装置1において異なる時刻で各々計測されて得られる複数の3次元点群データDpgから各々段差検出部23によって検出される複数の段差データ同士を比較する。例えば、比較部28は、走行型計測装置1において異なる時刻で各々計測されて得られる複数の3次元点群データDpgから各々段差検出部23によって検出される段差の高さ同士を比較し、段差の高さの経年変化を判定する。また、比較部28は、異なる時刻で各々計測されて得られる複数の3次元点群データDpgから各々段差検出部23によって検出される段差の広がり同士を比較し、段差の広がりの経年変化を判定する。表示処理部27は、比較部28によって判定された段差の高さの経年変化と段差の広がりの経年変化とを示す情報を表示部14に表示させることができる。
The comparison unit 28 compares a plurality of step data detected by the
比較部28は、例えば、前回計測された3次元点群データDpgを用いて検出される段差の高さと、今回計測された3次元点群データDpgを用いて検出される段差の高さとを比較し、高さが変化した段差を判定する。また、比較部28は、例えば、前回計測された3次元点群データDpgを用いて検出される段差の広がりと、今回計測された3次元点群データDpgを用いて検出される段差の広がりとを比較し、広がりが変化した段差を判定する。表示処理部27は、比較部28によって高さが変化したと判定された段差および比較部28によって広がりが変化したと判定された段差とを強調表示して表示部14に表示させることができる。 The comparison unit 28 compares, for example, the height of the step detected using the previously measured 3D point cloud data Dpg with the height of the step detected using the 3D point cloud data Dpg measured this time. Then, the step whose height has changed is determined. Further, in the comparison unit 28, for example, the spread of the step detected using the previously measured three-dimensional point cloud data Dpg and the spread of the step detected using the three-dimensional point cloud data Dpg measured this time. To determine the level difference in which the spread has changed. The display processing unit 27 can highlight and display the step determined by the comparison unit 28 that the height has changed and the step determined by the comparison unit 28 that the spread has changed.
比較部28は、例えば、前回計測された3次元点群データDpgを用いて検出される段差の高さと、今回計測された3次元点群データDpgを用いて検出される段差の高さとの差分を算出することができる。また、比較部28は、例えば、前回計測された3次元点群データDpgを用いて検出される段差の広がりと、今回計測された3次元点群データDpgを用いて検出される段差の広がりとの差分を算出することができる。表示処理部27は、比較部28によって算出された段差の高さの差分を示す情報および段差の広がりの差分を示す情報を表示部14に表示させることができる。表示処理部27は、比較部28によって高さが変化したと判定された段差の高さの差分および比較部28によって広がりが変化したと判定された段差の広がりの差分とを強調表示して表示部14に表示させることもできる。 In the comparison unit 28, for example, the difference between the height of the step detected using the previously measured 3D point cloud data Dpg and the height of the step detected using the 3D point cloud data Dpg measured this time. Can be calculated. Further, in the comparison unit 28, for example, the spread of the step detected using the previously measured three-dimensional point cloud data Dpg and the spread of the step detected using the three-dimensional point cloud data Dpg measured this time. The difference between can be calculated. The display processing unit 27 can display the information indicating the difference in the height of the steps calculated by the comparison unit 28 and the information indicating the difference in the spread of the steps on the display unit 14. The display processing unit 27 highlights and displays the difference in the height of the step determined by the comparison unit 28 to change the height and the difference in the spread of the step determined by the comparison unit 28 to change the spread. It can also be displayed on the unit 14.
統計処理部29は、トンネル5の内壁面5aにおいて設定された領域毎に、段差ベクトルの数および長さの少なくとも1つを段差ベクトル情報として算出する統計処理を行う。例えば、統計処理部29は、トンネル5のスパンであるトンネルスパン毎に、トンネルスパンに含まれる段差ベクトルの数、トンネルスパンに含まれる段差ベクトルのトータルの長さ、トンネルスパンに含まれる段差ベクトルにおける段差の高さ、およびトンネルスパンに含まれる段差ベクトルの広がりなどを算出する。トンネルスパンに含まれる段差ベクトルにおける段差の高さは、例えば、複数の高さの段差から段差ベクトルが生成される場合、段差ベクトルの生成に用いられる複数の段差の高さの平均値、最大値、および最小値の少なくとも一つを含む。
The
また、統計処理部29は、トンネルスパン毎に、段差ベクトルの分布を算出することができる。段差ベクトルの分布とは、例えば、トンネルスパンがトンネル周方向に複数の領域に区分されたと仮定した場合において区分された各領域に存在する段差ベクトルの割合を意味する。
In addition, the
表示処理部27は、統計処理部29によって算出された結果を表示部14に表示させることができる。図19は、実施の形態2にかかる統計処理部によって算出された結果の一例を示す図である。図19に示すように、統計処理部29は、統計処理部29によって算出された各トンネルスパンの段差ベクトルの数を示すグラフを表示部14に表示させることができる。図19に示すグラフでは、2017年で発生した段差ベクトルの数、2018年で発生した段差ベクトルの数、および2019年で発生した段差ベクトルの数が、トンネルスパン毎に示されている。
The display processing unit 27 can display the result calculated by the
また、統計処理部29は、各トンネルスパンの段差ベクトルの数に代えて、各トンネルスパンの段差ベクトルのトータルの長さを示すグラフ、または各トンネルスパンの段差ベクトルの分布を示すグラフなどを表示部14に表示させることができる。
Further, the
図20は、実施の形態2にかかる段差検出装置の処理部による処理手順の一例を示すフローチャートである。図20のステップS30からS37の処理は、図15に示すステップS10からS17の処理と同じであるため、説明を省略する。図20に示す処理は、処理部12Aによって繰り返し実行される。 FIG. 20 is a flowchart showing an example of a processing procedure by the processing unit of the step detection device according to the second embodiment. Since the processes of steps S30 to S37 of FIG. 20 are the same as the processes of steps S10 to S17 shown in FIG. 15, the description thereof will be omitted. The process shown in FIG. 20 is repeatedly executed by the processing unit 12A.
図20に示すように、段差検出装置2Aの処理部12Aは、比較処理タイミングになったか否かを判定する(ステップS38)。例えば、処理部12Aは、入力部13への特定操作があった場合に、比較処理タイミングになったと判定する。処理部12Aは、比較処理タイミングになったと判定した場合(ステップS38:Yes)、異なる時刻で各々計測されて得られる複数の3次元点群データDpgから各々段差検出部23によって検出される複数の段差データ同士を比較する比較処理を行い、比較処理の結果を表示部14に表示させる表示処理を行う(ステップS39)。
As shown in FIG. 20, the processing unit 12A of the step detection device 2A determines whether or not the comparison processing timing has been reached (step S38). For example, the processing unit 12A determines that the comparison processing timing has come when there is a specific operation on the
処理部12Aは、ステップS39の処理が終了した場合、または比較処理タイミングになっていないと判定した場合(ステップS38:No)、統計処理タイミングになったか否かを判定する(ステップS40)。例えば、処理部12Aは、入力部13への特定操作があった場合に、統計処理タイミングになったと判定する。処理部12Aは、統計処理タイミングになったと判定した場合(ステップS40:Yes)、段差ベクトルの数および長さの少なくとも1つをトンネル5の内壁面5aの設定された領域毎に段差ベクトル情報として算出する統計処理を行い、統計処理の結果を表示部14に表示させる表示処理を行う(ステップS41)。処理部12Aは、ステップS41の処理が終了した場合、または統計処理タイミングになっていないと判定した場合(ステップS40:No)、図20に示す処理を終了する。
When the processing of step S39 is completed or when it is determined that the comparison processing timing has not been reached (step S38: No), the processing unit 12A determines whether or not the statistical processing timing has been reached (step S40). For example, the processing unit 12A determines that the statistical processing timing has come when there is a specific operation on the
実施の形態2にかかる段差検出装置2Aのハードウェア構成例は、図17に示す段差検出装置2のハードウェア構成と同じである。プロセッサ101は、メモリ102に記憶されたプログラムを読み出して実行することによって、比較部28および統計処理部29の機能を実行することができる。
The hardware configuration example of the step detection device 2A according to the second embodiment is the same as the hardware configuration of the
以上のように、実施の形態2にかかる変状検出システム100Aの段差検出装置2Aは、統計処理部29を備える。統計処理部29は、段差ベクトル82a,82b,82c,82dの数および長さの少なくとも1つをトンネル5の内壁面5aのトンネルスパン毎に段差ベクトル情報として算出する。トンネルスパンは、予め定められた領域の一例である。表示処理部27は、統計処理部29で算出された段差ベクトル情報をグラフ化して表示部14に表示させる。これにより、段差の状態をトンネルスパン毎に把握することができる。
As described above, the step detection device 2A of the deformation detection system 100A according to the second embodiment includes a
また、段差検出装置2Aは、複数の3次元点群データDpgから各々段差検出部23によって検出される複数の段差のデータ同士を比較する比較部28を備える。これにより、例えば、段差の経年変化を把握することができる。
Further, the step detection device 2A includes a comparison unit 28 that compares the data of the plurality of steps detected by the
実施の形態3.
実施の形態3にかかる変状検出システムは、変状サーバおよびデータサーバなどが用いられる点で、実施の形態1にかかる変状検出システム100と異なる。以下においては、実施の形態1にかかる変状検出システム100と同様の構成については同一符号を付して説明を省略し、実施の形態1にかかる変状検出システム100と異なる構成について主に説明するものとする。Embodiment 3.
The deformation detection system according to the third embodiment is different from the
図21は、本発明の実施の形態3にかかる変状検出システムの第1の構成例を示す図である。図21に示す変状検出システム100Bは、走行型計測装置1と、段差検出装置2Bと、変状サーバ4とを備える。変状サーバ4は、通信部40と、記憶部41と、処理部42とを備える。処理部42は、走行型計測装置1からネットワーク3および通信部40を介して3次元点群データDpgおよび画像データDimgを取得し、取得した3次元点群データDpgおよび画像データDimgを記憶部41に記憶させる。記憶部41は、取得部21Bによって取得される3次元点群データDpgおよび画像データDimgの他、形状データDfと変状データDdfとを記憶する。
FIG. 21 is a diagram showing a first configuration example of the deformation detection system according to the third embodiment of the present invention. The
段差検出装置2Bは、取得部21を有する処理部12に代えて、取得部21Bを有する処理部12Bを備える。取得部21Bは、変状サーバ4から、3次元点群データDpg、画像データDimg、形状データDf、および変状データDdfをネットワーク3および通信部10を介して取得する点で、取得部21と異なる。また、取得部21Bは、段差検出部23および段差評価部24によって生成された段差データを通信部10からネットワーク3を介して変状サーバ4へ送信する。変状サーバ4は、段差検出装置2Bから取得した段差データを記憶部41に記憶された変状データDdfへ追加して、変状データDdfを更新する。
The
このように、図21に示す変状検出システム100Bでは、データ量が多い3次元点群データDpgなどが変状サーバ4に記憶される。これにより、実施の形態1にかかる段差検出装置2の記憶部11の記憶容量に比べ、段差検出装置2Bの記憶部11の記憶容量を低減することができる。
As described above, in the
図22は、実施の形態3にかかる変状検出システムの第2の構成例を示す図である。図22に示す変状検出システム100Cは、走行型計測装置1と、段差検出装置2Cと、変状サーバ4Cとを備える。変状サーバ4Cは、通信部40と、記憶部41と、処理部42Cとを備える。処理部42Cは、データ変換部51を備えており、上述した処理部42の機能に加え、データ変換部51によって3次元点群データDpgを要素点群データDeへ変換する機能を有する。データ変換部51の機能は、図4に示すデータ変換部22の機能と同じである。
FIG. 22 is a diagram showing a second configuration example of the deformation detection system according to the third embodiment. The
段差検出装置2Cは、取得部21Bを有する処理部12Bに代えて、取得部21Cを有する処理部12Cを備える。取得部21Cは、変状サーバ4Cから、3次元点群データDpgおよび形状データDfに代えて、要素点群データDeをネットワーク3および通信部10を介して取得する点で、取得部21Bと異なる。段差検出装置2Cは、段差検出装置2Bと同様に、生成した段差データを変状サーバ4Cへ送信する。
The step detection device 2C includes a processing unit 12C having an acquisition unit 21C instead of the processing unit 12B having the acquisition unit 21B. The acquisition unit 21C differs from the acquisition unit 21B in that the element point cloud data De is acquired from the
このように、図22に示す変状検出システム100Cでは、データ量が多い3次元点群データDpgなどが変状サーバ4Cに記憶され、3次元点群データDpgが要素点群データDeへ変状サーバ4Cによって変換される。これにより、段差検出装置2Cの記憶部11の記憶容量を低減することができる。また、要素点群データDeは、3次元点群データDpgよりもデータ容量が小さいことから、変状検出システム100Cは、変状検出システム100Bに比べて、段差検出装置2Cと変状サーバ4Cとの間で通信量を削減することができる。
As described above, in the
図23は、実施の形態3にかかる変状検出システムの第3の構成例を示す図である。図23に示す変状検出システム100Dは、走行型計測装置1と、変状サーバ4Dと、端末装置6とを備える。変状サーバ4Dは、段差検出装置2の機能のうち表示機能および入力受付機能を除く機能が実行される。変状サーバ4Dは、データ変換部51と、段差検出部52と、段差評価部53と、属性設定部54と、座標変換部55とを備える処理部42Dを有する。
FIG. 23 is a diagram showing a third configuration example of the deformation detection system according to the third embodiment. The
データ変換部51は、データ変換部22と同様の機能を有しており、記憶部41に記憶された形状データDfおよび3次元点群データDpgに基づいて、3次元点群データDpgを要素点群データDeへ変換する。段差検出部52は、段差検出部23と同様に、要素点群データDeから段差を検出し、検出した段差に関する情報を含む段差データを生成する。段差検出部52は、生成した段差データを変状データDdfへ追加する。
The
段差評価部53は、段差評価部24と同様に、段差の状態レベルを判定し、判定した段差レベルを示す情報を段差データに追加する。属性設定部54は、端末装置6からの要求に基づいて、属性設定部25と同様に、段差データに段差属性を示す情報を追加する。座標変換部55は、座標変換部26と同様に、記憶部41に記憶された2次元座標のデータを3次元座標のデータへ変換し、記憶部41に記憶された3次元座標のデータを2次元座標のデータへ変換する。座標変換部55は、変換したデータを記憶部41へ記憶させる。
Similar to the step evaluation unit 24, the step evaluation unit 53 determines the state level of the step and adds information indicating the determined step level to the step data. Based on the request from the terminal device 6, the attribute setting unit 54 adds information indicating the step attribute to the step data, similarly to the
端末装置6は、通信部60と、処理部61と、入力部62と、表示部63とを備える。処理部61は、要求処理部70と、表示処理部71と、出力処理部72とを備える。要求処理部70は、入力部62への特定操作があった場合、特定操作に対応するデータの要求を変状サーバ4Dへ通信部60を介して行う。要求処理部70は、通信部60を介して変状サーバ4Dから要求に応じたデータを取得する。
The terminal device 6 includes a communication unit 60, a processing unit 61, an input unit 62, and a
表示処理部71は、要求処理部70によって取得されたデータを表示部63に表示させる。表示処理部71によって表示されるデータは、表示処理部27によって表示されるデータと同じである。これにより、端末装置6は、段差検出装置2と同様のデータを表示部63に表示させることができる。また、出力処理部72は、入力部62によって入力された段差の属性情報を変状サーバ4Dへ通信部60から送信する。これにより、変状サーバ4Dによって段差の属性情報が段差データに追加される。
The display processing unit 71 causes the
このように、変状検出システム100Dでは、変状サーバ4Dによって主な処理が行われることから、端末装置6の処理負荷を軽減することができる。
As described above, in the
図24は、実施の形態3にかかる変状検出システムの第4の構成例を示す図である。図24に示す変状検出システム100Eは、走行型計測装置1と、変状サーバ4Eと、端末装置6と、データサーバ7を備える。変状検出システム100Eでは、形状データDf、3次元点群データDpg、画像データDimg、および変状データDdfが変状サーバ4Eではなくデータサーバ7に記憶される点で変状検出システム100Dと異なる。
FIG. 24 is a diagram showing a fourth configuration example of the deformation detection system according to the third embodiment. The
変状サーバ4Eは、ネットワーク3を介してデータサーバ7との間で、形状データDf、3次元点群データDpg、画像データDimg、および変状データDdfの送受信を行う。このように、変状検出システム100Eは、変状サーバ4Eが処理を行うサーバとして用いられ、データサーバ7が変状サーバ4Eまたは端末装置6で用いられるデータを記憶するデータサーバとして用いられる。
The
以上のように、実施の形態3にかかる変状検出システム100B,100C,100D,100Eでは、実施の形態1にかかる段差検出装置2が有する複数の機能の一部または全部を複数の装置に分散させている。これにより、システム設計の自由度を高めることができる。また、同様に、実施の形態2にかかる段差検出装置2Aが有する比較部28および統計処理部29の各々の機能を変状サーバ4,4C,4D,4Eに設けることもできる。
As described above, in the
実施の形態3にかかる段差検出装置2B,2Cおよび端末装置6の各々のハードウェア構成例は、図17に示す段差検出装置2のハードウェア構成と同じである。プロセッサ101は、メモリ102に記憶されたプログラムを読み出して実行することによって、段差検出装置2B,2Cおよび端末装置6の各々の各機能を実行することができる。
Each of the hardware configuration examples of the
また、実施の形態3にかかる変状サーバ4,4C,4D,4Eおよびデータサーバ7の各々のハードウェア構成例は、図17に示す段差検出装置2のハードウェア構成のうち表示装置105を除いた構成である。プロセッサ101は、メモリ102に記憶されたプログラムを読み出して実行することによって、変状サーバ4,4C,4D,4Eおよびデータサーバ7の各々の各機能を実行することができる。
Further, in the hardware configuration examples of the
以上の実施の形態に示した構成は、本発明の内容の一例を示すものであり、別の公知の技術と組み合わせることも可能であるし、本発明の要旨を逸脱しない範囲で、構成の一部を省略、変更することも可能である。 The configuration shown in the above-described embodiment shows an example of the content of the present invention, can be combined with another known technique, and is one of the configurations without departing from the gist of the present invention. It is also possible to omit or change the part.
1 走行型計測装置、2,2A,2B,2C 段差検出装置、3 ネットワーク、4,4C,4D,4E 変状サーバ、5 トンネル、5a 内壁面、6 端末装置、7 データサーバ、10,40,60 通信部、11,41 記憶部、12,12A,12B,12C,42,42C,42D,61 処理部、13,62 入力部、14,63 表示部、21,21B,21C 取得部、22,51 データ変換部、23,52 段差検出部、24,53 段差評価部、25,54 属性設定部、26,55 座標変換部、27,71 表示処理部、28 比較部、29 統計処理部、31 段差位置検出部、32 段差高検出部、33 ベクトル生成部、70 要求処理部、72 出力処理部、80 変状展開図、81 外形枠、82a,82b,82c,82d,92a,92b,92c,92d 段差ベクトル、83a,83b,83c,83d 属性情報、90 3次元画像、91 3次元点群、100,100A,100B,100C,100D,100E 変状検出システム。 1 Travel type measuring device, 2,2A, 2B, 2C step detection device, 3 network, 4,4C, 4D, 4E deformation server, 5 tunnel, 5a inner wall surface, 6 terminal device, 7 data server, 10,40, 60 Communication unit, 11,41 Storage unit, 12,12A, 12B, 12C, 42, 42C, 42D, 61 Processing unit, 13,62 Input unit, 14,63 Display unit, 21,21B, 21C Acquisition unit, 22, 51 Data conversion unit, 23,52 Step detection unit, 24,53 Step evaluation unit, 25,54 Attribute setting unit, 26,55 Coordinate conversion unit, 27,71 Display processing unit, 28 Comparison unit, 29 Statistical processing unit, 31 Step position detection unit, 32 step height detection unit, 33 vector generation unit, 70 request processing unit, 72 output processing unit, 80 deformation development drawing, 81 outline frame, 82a, 82b, 82c, 82d, 92a, 92b, 92c, 92d step vector, 83a, 83b, 83c, 83d attribute information, 90 3D image, 91 3D point cloud, 100, 100A, 100B, 100C, 100D, 100E Deformity detection system.
Claims (14)
前記形状データに基づいて前記構造物の表面に対応する基準面を設定し、前記3次元点群データを、前記基準面に行列状に配列される複数の点から選択される1つの点の位置と前記基準面からのずれ量とで各々表される複数の要素点のデータを含む要素点群データへ変換するデータ変換部と、
前記要素点群データに含まれる前記複数の要素点に対してガウシアン差分フィルタをかけて得られる結果に基づいて、前記構造物の表面における段差の位置を検出する段差検出部と、を備える
ことを特徴とする変状検出システム。An acquisition unit that acquires three-dimensional point cloud data, which is data representing the surface of a structure as a three-dimensional point cloud, and shape data indicating the shape of the structure.
A reference plane corresponding to the surface of the structure is set based on the shape data, and the position of one point selected from a plurality of points arranged in a matrix on the reference plane for the three-dimensional point cloud data. And a data conversion unit that converts data into element point cloud data including data of a plurality of element points represented by the amount of deviation from the reference plane.
It is provided with a step detection unit that detects the position of a step on the surface of the structure based on the result obtained by applying a Gaussian difference filter to the plurality of element points included in the element point cloud data. Characteristic deformation detection system.
前記ガウシアン差分フィルタをかけた後の前記複数の要素点のうち互いに隣接し且つ前記ガウシアン差分フィルタをかけて得られる値の極性が異なる2つの要素点同士を結ぶ線に零交差が生じる位置を特定し、特定した位置を前記段差の位置として検出する
ことを特徴とする請求項1に記載の変状検出システム。The step detection unit
Specify the position where zero intersection occurs in the line connecting two element points that are adjacent to each other and have different polarities in the values obtained by applying the Gaussian difference filter among the plurality of element points after the Gaussian difference filter is applied. The deformation detection system according to claim 1, wherein the specified position is detected as the position of the step.
前記ガウシアン差分フィルタをかけた後の前記複数の要素点のうち隣接する3以上の要素点であって少なくとも1つの要素点と残りの要素点とで前記ガウシアン差分フィルタをかけて得られる値の極性が異なる3以上の要素点を含む面の傾きに基づいて、前記段差の高さを検出する段差高検出部を備える
ことを特徴とする請求項1または2に記載の変状検出システム。The step detection unit
Polarity of a value obtained by applying the Gaussian difference filter to at least one element point and the remaining element points of three or more adjacent element points among the plurality of element points after the Gaussian difference filter is applied. The deformation detection system according to claim 1 or 2, further comprising a step height detection unit that detects the height of the step based on the inclination of a surface including three or more element points having different values.
前記2つの要素点同士を結ぶ線の傾きに基づいて前記段差の高さを検出する段差高検出部を備える
ことを特徴とする請求項2に記載の変状検出システム。The step detection unit
The deformation detection system according to claim 2, further comprising a step height detecting unit that detects the height of the step based on the inclination of a line connecting the two element points.
前記段差の位置を結んで得られるベクトルである段差ベクトルを生成するベクトル生成部を備える
ことを特徴とする請求項1から4のいずれか一つに記載の変状検出システム。The step detection unit
The deformation detection system according to any one of claims 1 to 4, further comprising a vector generation unit that generates a step vector, which is a vector obtained by connecting the positions of the steps.
ことを特徴とする請求項5に記載の変状検出システム。The deformation detection system according to claim 5, further comprising a display processing unit that displays the step vector generated by the vector generation unit on the display unit.
前記構造物の展開図および撮像画像のうち少なくとも一つに前記段差ベクトルを重畳した画像を前記表示部に表示させる
ことを特徴とする請求項6に記載の変状検出システム。The display processing unit
The deformation detection system according to claim 6, wherein an image in which the step vector is superimposed on at least one of the developed view of the structure and the captured image is displayed on the display unit.
前記表示処理部は、
前記3次元座標系のデータに変換された前記段差ベクトルを前記表示部に表示させる
ことを特徴とする請求項6または7に記載の変状検出システム。A coordinate conversion unit that converts the step vector into data in a three-dimensional coordinate system is provided.
The display processing unit
The deformation detection system according to claim 6 or 7, wherein the step vector converted into data of the three-dimensional coordinate system is displayed on the display unit.
前記構造物の表面を表す3次元点群の画像に前記3次元座標系のデータに変換された前記段差ベクトルを重畳した画像を前記表示部に表示させる
ことを特徴とする請求項8に記載の変状検出システム。The display processing unit
The eighth aspect of claim 8, wherein the display unit displays an image in which the step vector converted into the data of the three-dimensional coordinate system is superimposed on the image of the three-dimensional point cloud representing the surface of the structure. Deformity detection system.
前記表示処理部は、
前記統計処理部で算出された前記段差ベクトル情報をグラフ化して前記表示部に表示させる
ことを特徴とする請求項9に記載の変状検出システム。A statistical processing unit that calculates at least one of the number and length of the step vector as step vector information for each predetermined region on the surface of the structure is provided.
The display processing unit
The deformation detection system according to claim 9, wherein the step vector information calculated by the statistical processing unit is graphed and displayed on the display unit.
異なる時刻で各々前記構造物の表面を計測して得られた複数の3次元点群データを取得し、
前記複数の3次元点群データは、
前記構造物の表面を計測した時刻を示す時刻情報を含み、
前記段差検出部は、
前記複数の3次元点群データの各々から前記段差を検出し、検出した前記段差に前記時刻情報を関連付け、
前記表示処理部は、
前記段差検出部によって検出された前記段差の情報を前記時刻情報に基づいて時系列に前記表示部に表示させる
ことを特徴とする請求項6から10のいずれか一つに記載の変状検出システム。The acquisition unit
A plurality of three-dimensional point cloud data obtained by measuring the surface of the structure at different times are acquired.
The plurality of three-dimensional point cloud data are
Includes time information indicating the time when the surface of the structure was measured.
The step detection unit
The step is detected from each of the plurality of three-dimensional point cloud data, and the time information is associated with the detected step.
The display processing unit
The deformation detection system according to any one of claims 6 to 10, wherein the information on the step detected by the step detection unit is displayed on the display unit in chronological order based on the time information. ..
ことを特徴とする請求項11に記載の変状検出システム。The deformation detection system according to claim 11, further comprising a comparison unit that compares the data of the plurality of steps detected by the step detection unit from the plurality of three-dimensional point cloud data.
ことを特徴とする請求項1から12のいずれか一つに記載の変状検出システム。One of claims 1 to 12, wherein a step evaluation unit for determining a state level representing the state of the step based on the state of the step detected by the step detection unit is provided. Deformation detection system described.
構造物の表面を3次元点群で表すデータである3次元点群データと前記構造物の形状を示す形状データとを取得する取得ステップと、
前記形状データに基づいて前記構造物の表面に対応する基準面を設定し、前記3次元点群データを、前記基準面に行列状に配列される複数の点から選択される1つの点の位置と前記基準面からのずれ量とで各々表される複数の要素点のデータを含む要素点群データへ変換するデータ変換ステップと、
前記要素点群データに含まれる前記複数の要素点に対してガウシアン差分フィルタをかけて得られる結果に基づいて、前記構造物の表面における段差の位置を検出する段差検出ステップと、を含む
ことを特徴とする変状検出方法。It ’s a computer-executed error detection method.
An acquisition step for acquiring 3D point cloud data, which is data representing the surface of a structure in a 3D point cloud, and shape data indicating the shape of the structure, and
A reference plane corresponding to the surface of the structure is set based on the shape data, and the position of one point selected from a plurality of points arranged in a matrix on the reference plane for the three-dimensional point cloud data. And a data conversion step of converting to element point cloud data including data of a plurality of element points represented by the amount of deviation from the reference plane.
It includes a step detection step of detecting the position of a step on the surface of the structure based on the result obtained by applying a Gaussian difference filter to the plurality of element points included in the element point cloud data. A characteristic deformation detection method.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2019/035261 WO2021044628A1 (en) | 2019-09-06 | 2019-09-06 | Deformation detection system and deformation detection method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP6664569B1 JP6664569B1 (en) | 2020-03-13 |
JPWO2021044628A1 true JPWO2021044628A1 (en) | 2021-09-27 |
Family
ID=70000328
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019567758A Active JP6664569B1 (en) | 2019-09-06 | 2019-09-06 | Deformation detection system and deformation detection method |
Country Status (2)
Country | Link |
---|---|
JP (1) | JP6664569B1 (en) |
WO (1) | WO2021044628A1 (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113673011A (en) * | 2021-06-04 | 2021-11-19 | 浙江华东测绘与工程安全技术有限公司 | Method for intelligently identifying tunnel invasion boundary in operation period based on point cloud data |
WO2023058098A1 (en) * | 2021-10-04 | 2023-04-13 | 日本電気株式会社 | Display device, display method, and non-transitory computer-readable medium |
WO2024195097A1 (en) * | 2023-03-23 | 2024-09-26 | 日本電気株式会社 | Information processing device, information processing method, and storage medium |
CN118640838B (en) * | 2024-08-15 | 2024-10-22 | 济冠钢铁(山东)有限公司 | Surface digital detection method and system for rolling mill sheet |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3208273B2 (en) * | 1995-02-14 | 2001-09-10 | 株式会社東芝 | Surface profile measuring device |
JP5161936B2 (en) * | 2010-08-11 | 2013-03-13 | 株式会社パスコ | Data analysis apparatus, data analysis method, and program |
JP5857397B2 (en) * | 2010-11-02 | 2016-02-10 | トヨタ自動車株式会社 | Step information acquisition system and step information acquisition method |
JP2018115937A (en) * | 2017-01-17 | 2018-07-26 | リコーエレメックス株式会社 | Inspection system |
-
2019
- 2019-09-06 WO PCT/JP2019/035261 patent/WO2021044628A1/en active Application Filing
- 2019-09-06 JP JP2019567758A patent/JP6664569B1/en active Active
Also Published As
Publication number | Publication date |
---|---|
WO2021044628A1 (en) | 2021-03-11 |
JP6664569B1 (en) | 2020-03-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
WO2021044628A1 (en) | Deformation detection system and deformation detection method | |
ES2973165T3 (en) | Method and device for filling terrain elevation model data regions | |
JP4917351B2 (en) | Calibration method in three-dimensional shape measuring apparatus | |
JP6359868B2 (en) | 3D data display device, 3D data display method, and 3D data display program | |
JP6565037B2 (en) | Displacement measuring device, displacement measuring method and program thereof | |
Gao et al. | An approach to combine progressively captured point clouds for BIM update | |
JP2008185375A (en) | 3d shape calculation device of sar image, and distortion correction device of sar image | |
CN105306922A (en) | Method and device for obtaining depth camera reference diagram | |
CN117195361B (en) | High slope stability monitoring method, device, equipment and readable storage medium | |
CN112308969B (en) | Three-dimensional map generation method and device | |
CN116152218A (en) | Intelligent detection method and device for construction quality | |
JP2015184061A (en) | Extracting device, method, and program | |
CN112348781A (en) | Method, device and equipment for detecting height of reference plane and storage medium | |
CN114897959A (en) | Phase unwrapping method based on light field multi-view constraint and related components | |
KR101814023B1 (en) | Apparatus and Method for Automatic Calibration of Finite Difference Grid Data | |
JP2019164109A (en) | Measuring device, measuring system, measuring method, and program | |
JP6991700B2 (en) | Information processing equipment, information processing method, program | |
JP5496033B2 (en) | Defect dimension measuring apparatus, defect dimension measuring method, and program | |
KR101808958B1 (en) | Method for obtaining shape information of structure and method for measuring deformation of structure | |
KR101621858B1 (en) | Apparatus and method for calculating horizontal distance between peak and structure point | |
JP2007264952A (en) | Ground-analyzing mesh generation method and ground-analyzing mesh generation program | |
JP7285678B2 (en) | Information processing device and control program | |
CN113421292A (en) | Three-dimensional modeling detail enhancement method and device | |
CN111145347B (en) | Method and device for correcting digital elevation data of terrain section and correction equipment | |
CN113219489B (en) | Point-to-point determination method, device, computer equipment and storage medium for multi-line laser |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20191206 |
|
A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20191206 |
|
A975 | Report on accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A971005 Effective date: 20200117 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20200121 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20200218 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6664569 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |