JP2000121358A - Method and device for obtaining inspection point position - Google Patents

Method and device for obtaining inspection point position

Info

Publication number
JP2000121358A
JP2000121358A JP10297298A JP29729898A JP2000121358A JP 2000121358 A JP2000121358 A JP 2000121358A JP 10297298 A JP10297298 A JP 10297298A JP 29729898 A JP29729898 A JP 29729898A JP 2000121358 A JP2000121358 A JP 2000121358A
Authority
JP
Japan
Prior art keywords
image
plane
inspection
coordinates
coefficient
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
Application number
JP10297298A
Other languages
Japanese (ja)
Other versions
JP4170469B2 (en
Inventor
Kazuo Oda
和夫 織田
Takeshi Kondo
剛 近藤
Masayoshi Ohata
昌芳 尾幡
Takeshi Doihara
健 土居原
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Asia Air Survey Co Ltd
Original Assignee
Asia Air Survey Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Asia Air Survey Co Ltd filed Critical Asia Air Survey Co Ltd
Priority to JP29729898A priority Critical patent/JP4170469B2/en
Publication of JP2000121358A publication Critical patent/JP2000121358A/en
Application granted granted Critical
Publication of JP4170469B2 publication Critical patent/JP4170469B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Landscapes

  • Image Input (AREA)
  • Image Processing (AREA)
  • Image Analysis (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Abstract

PROBLEM TO BE SOLVED: To find the position of an arbitrary point in an inspection area included in the plane of a structure to be inspected on the plane. SOLUTION: An actual size coordinate system is set by installing four indexes on the plane structure 4. Two kinds of wide-area and narrow-area digital images of the plane structure are obtained by a camera 1 and a film reader 2. A data process part 3 calculates a first projection conversion coefficient from the pixel coordinates and actual size coordinates of the four indexes in the wide-area image. A second projection conversion coefficient is calculated on the basis of coordinate information on two corresponding places in the narrow-area and wide-area images. For the narrow-area image, an arbitrary point to be inspected is inputted and coordinate conversion is performed by using the first projection conversion coefficient and second projection conversion coefficient.

Description

【発明の詳細な説明】DETAILED DESCRIPTION OF THE INVENTION

【0001】[0001]

【発明の属する技術分野】本発明は、検査対象構造物の
平面に含まれる検査領域の任意点の位置を求める方法お
よびその装置に関する。
BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a method and an apparatus for determining a position of an arbitrary point in an inspection area included in a plane of a structure to be inspected.

【0002】[0002]

【従来の技術】長大平面構造物(典型的にはコンクリー
ト壁面等)の整備点検においては、一般的に外観が重要
な良否判断要素となる。例えば、ひび割れの長さ、幅や
水漏れなどの有無等は、構造物の外観から判読される。
近年では、平面構造物の外観画像を写真撮影した後、コ
ンピュータによりデジタル画像処理や解析測定を行い点
検する方法が開発されている。この方法を使うと、ディ
スプレイ画面を見ながら遠隔的に点検を行うことができ
る。
2. Description of the Related Art In the maintenance and inspection of a long and long planar structure (typically, a concrete wall or the like), appearance is generally an important quality judgment factor. For example, the length, width, presence or absence of water leakage, and the like of a crack can be read from the appearance of the structure.
In recent years, a method has been developed in which a photograph of an external appearance image of a planar structure is taken and then inspected by performing digital image processing and analysis measurement with a computer. By using this method, the inspection can be performed remotely while looking at the display screen.

【0003】このようなコンピュータ解析による点検方
法では、画像のA−D変換時の分解能が重要となる。例
えば、0.1mm 幅程度の細いひび割れを検出するために
は、デジタル画像上でのピクセル間距離で0.025mm 程度
の分解能が必要である。
In such an inspection method based on computer analysis, the resolution at the time of A / D conversion of an image is important. For example, in order to detect a small crack having a width of about 0.1 mm, a resolution of about 0.025 mm is required as a distance between pixels on a digital image.

【0004】この際に問題になるのはデジタル画像のサ
イズである。例えば、上記分解能で50cm×50cmをカバー
しようとすると、それだけで20000 ×20000 のピクセル
数が必要になる。1ピクセルあたり8ビットの階調でこ
の画像を保存する場合は、約400MByteの記憶容量が必要
となる。従って、数m〜数10mのオーダーの大きさの平
面構造物全体をこの解像度および階調で保存する場合、
ギガbyteオーダーのメモリが必要である。また画像処理
を実行するコンピュータにも高い処理能力が要求され
る。このような理由から従来は、構造物の外観から要点
検箇所の存在を確認し、その場所の画像だけを高い分解
能で取り込んでコンピュータ画像解析を行っていた。
[0004] At this time, a problem is the size of the digital image. For example, if it is intended to cover 50 cm × 50 cm with the above resolution, the number of pixels of 20000 × 20000 is required alone. If this image is stored at a gradation of 8 bits per pixel, a storage capacity of about 400 MByte is required. Therefore, when storing the entire planar structure having a size on the order of several meters to several tens of meters at this resolution and gradation,
A gigabyte order memory is required. In addition, a computer that executes image processing also requires high processing capability. For these reasons, conventionally, the presence of a point requiring inspection has been confirmed from the appearance of the structure, and only the image of that location has been captured with high resolution and computer image analysis has been performed.

【0005】[0005]

【発明が解決しようとする課題】上記のように従来のひ
び割れの画像解析は、ひび割れ部分の画像を高い解像度
で取得するものなので、ひび割れ自体の解析には適して
いるが、一方で、ひび割れが検査対象の構造物のどの位
置を占めるかが問題となることがある。
As described above, the conventional image analysis of a crack is suitable for analyzing the crack itself because the image of the crack is obtained at a high resolution. The position of the structure to be inspected may be a problem.

【0006】例えば、ダムにおいてひび割れらしきもの
が発見された場合、その位置次第では、実際に現地のダ
ム壁に登り、その部分の精細調査を行う必要性が生じ
る。しかし、このような精細調査には膨大な時間と労力
を要すため、予め、ひび割れのダム壁における位置を確
認しておきたいという要望があった。
[0006] For example, when a thing that looks like a crack is found in a dam, depending on its position, it is necessary to actually climb the dam wall on the site and conduct a detailed investigation of that part. However, since such a detailed survey requires an enormous amount of time and effort, there has been a demand to confirm the position of the crack on the dam wall in advance.

【0007】この要望に対し、概略の位置取得だけなら
ば、人間がダム全体を見てひび割れ部分の位置を見当づ
けることもできる。そのためにはダムを正面から撮影し
て正対画像を得るのが理想だが、撮影場所の確保の面か
らこのような撮影は容易ではなく、通常は斜めの角度か
ら撮影が行われるので、写真を見てひび割れの位置を見
当づけるのは困難であった。
In response to this demand, if only a rough position acquisition is performed, a person can look at the entire dam and find the position of the cracked portion. For this purpose, it is ideal to take a picture of the dam from the front and obtain a facing image.However, such a shooting is not easy from the viewpoint of securing a shooting place, and usually the shooting is performed from an oblique angle, so the photograph is taken. It was difficult to see the location of the crack.

【0008】そこで、この発明は上記従来の課題に鑑み
なされたものであり、検査対象構造物の平面に含まれた
検査領域の任意点の平面における位置を求める検査点位
置取得方法および装置を提供することを目的とする。
The present invention has been made in view of the above-mentioned conventional problems, and provides an inspection point position obtaining method and apparatus for obtaining a position on an arbitrary plane of an inspection area included in a plane of an inspection target structure. The purpose is to do.

【0009】[0009]

【課題を解決するための手段】この発明による検査点位
置取得方法は、検査対象構造物の平面に含まれた検査領
域の任意点の前記平面における位置を求める検査点位置
取得方法であって、前記検査領域が撮像された検査画像
を読み込むステップと、予め前記平面に定められた複数
の基準点と前記検査領域とが撮像された平面画像を読み
込むステップと、予め測定された前記複数の基準点の位
置データを設定するステップと、前記平面画像における
前記複数の基準点の座標を認識し、その認識した座標と
前記設定された位置データとに基づいて、前記平面画像
の任意の座標を前記平面における位置に変換するための
係数を算出するステップと、前記検査画像および前記平
面画像に共通な複数の共通撮像点に基づいて、前記検査
画像の任意の座標を前記平面画像の座標に変換するため
の係数を算出するステップと、前記検査画像の範囲で任
意点が指定されたときは、その指定された任意点の前記
検査画像の座標を、前記平面画像の座標に変換するため
の係数と前記平面における位置に変換するための係数と
を用いて、前記検査対象構造物の平面における位置に変
換するステップと、を有することを特徴とする。
An inspection point position obtaining method according to the present invention is an inspection point position obtaining method for obtaining an arbitrary point of an inspection area included in a plane of a structure to be inspected on the plane. Reading an inspection image in which the inspection area is imaged, reading a plane image in which the inspection area is imaged with a plurality of reference points determined in advance on the plane, and reading the plurality of reference points measured in advance. Setting the position data, and recognizing the coordinates of the plurality of reference points in the plane image, and, based on the recognized coordinates and the set position data, changing any coordinates of the plane image to the plane. Calculating a coefficient for converting to a position in the inspection image, and based on a plurality of common imaging points common to the inspection image and the planar image, arbitrary coordinates of the inspection image Calculating a coefficient for converting to the coordinates of the plane image, and when an arbitrary point is specified in the range of the inspection image, the coordinates of the inspection image at the specified arbitrary point are calculated as the coordinates of the plane image. Using a coefficient for converting to coordinates and a coefficient for converting to a position on the plane, converting to a position on the plane of the structure to be inspected.

【0010】この発明による検査点位置取得方法は、検
査対象構造物の平面に含まれた検査領域の任意点の前記
平面における位置を求める検査点位置取得方法であり、
従来のコンピュータ解析のための検査領域画像は、通常
平面の部分画像であるから、この画像に加えて、撮像領
域を大きくしながら段階的に平面全体の画像まで取得
し、これら画像間を座標変換係数で対応づけようとした
ものである。
[0010] An inspection point position acquisition method according to the present invention is an inspection point position acquisition method for finding the position on the plane of an arbitrary point in an inspection area included in the plane of a structure to be inspected,
Conventionally, the inspection area image for computer analysis is a partial image of a plane, so in addition to this image, an image of the entire plane is acquired stepwise while enlarging the imaging area, and coordinate conversion between these images is performed. It is intended to be associated with a coefficient.

【0011】先ずひび割れ等が撮像された検査画像を先
ず読み込み、また予め平面に定められた複数の基準点と
検査領域とが撮像された平面画像(全体画像)を読み込
む。また、予め測定された複数の基準点の位置データを
コンピュータシステムなどに設定する。
First, an inspection image in which cracks and the like are captured is first read, and a planar image (whole image) in which a plurality of reference points and an inspection region determined in advance on a plane are captured is read. In addition, position data of a plurality of reference points measured in advance is set in a computer system or the like.

【0012】次に、オペレータが平面画像において指定
した複数の基準点の座標をコンピュータシステムなどに
より認識し、その認識した座標と設定された位置データ
とに基づいて、平面画像の任意の座標を平面における位
置に変換するための係数、即ち座標変換係数を算出す
る。これは、いわゆる射影変換係数であり、平面画像と
実際の平面との関係において、次の要素を含む変換を可
能にする。即ち、拡大縮小、回転、平行移動および奥行
き方向の傾きをもった変換が可能になる。
Next, the coordinates of a plurality of reference points specified by the operator on the plane image are recognized by a computer system or the like, and arbitrary coordinates of the plane image are determined based on the recognized coordinates and the set position data. Is calculated, ie, a coordinate conversion coefficient. This is a so-called projective transformation coefficient, which enables transformation including the following elements in the relationship between the plane image and the actual plane. That is, conversion with scaling, rotation, parallel movement, and inclination in the depth direction becomes possible.

【0013】そして、オペレータが指定した検査画像お
よび平面画像に共通な複数の共通撮像点に基づいて、コ
ンピュータシステムなどで演算を行い、検査画像の任意
の座標を平面画像の座標に変換するための係数を算出す
る。この係数も前記した「位置に変換するための係数」
と同様、射影変換係数(座標変換係数)である。
Then, a computer system or the like performs an operation based on a plurality of common imaging points common to the inspection image and the plane image designated by the operator to convert arbitrary coordinates of the inspection image into coordinates of the plane image. Calculate the coefficient. This coefficient is also referred to as the "coefficient for converting to position".
Similarly to the above, they are projection transformation coefficients (coordinate transformation coefficients).

【0014】そして検査画像の範囲でオペレータが位置
を求めたい任意点を指定したときは、その指定された任
意点の検査画像の座標を、予め算出された検査画像の任
意の座標を平面画像の座標に変換するための係数と、平
面画像の任意の座標を平面における位置に変換するため
の係数とを用いて、座標変換を段階的に行って、検査対
象構造物の平面における位置に変換する。従って、検査
画像の範囲において検査点として任意点を指定すれば、
平面における位置を求めることができる。
When the operator designates an arbitrary point whose position is desired to be obtained in the range of the inspection image, the coordinates of the inspection image at the specified arbitrary point are replaced with the arbitrary coordinates of the inspection image calculated in advance. Using a coefficient for converting to coordinates and a coefficient for converting arbitrary coordinates of a planar image to a position on a plane, coordinate conversion is performed stepwise to convert to a position on the plane of the inspection target structure. . Therefore, if an arbitrary point is designated as the inspection point in the range of the inspection image,
The position on the plane can be determined.

【0015】また、この発明の検査点位置取得装置は、
検査対象構造物の平面に含まれた検査領域の任意点の前
記平面における位置を求める検査点位置取得装置であっ
て、前記検査領域が撮像された検査画像を読み込む検査
画像読込手段と、予め前記平面に定められた複数の基準
点と前記検査領域とが撮像された平面画像を読み込む平
面画像読込手段と、予め測定された前記複数の基準点の
位置データを設定する基準点設定手段と、前記平面画像
における前記複数の基準点の座標を認識し、その認識し
た座標と前記設定された位置データとに基づいて、前記
平面画像の任意の座標を前記平面における位置に変換す
るための係数を算出する位置変換係数算出手段と、前記
検査画像および前記平面画像に共通な複数の共通撮像点
に基づいて、前記検査画像の任意の座標を前記平面画像
の座標に変換するための係数を算出する座標変換係数算
出手段と、前記検査画像の範囲で任意点が指定されたと
きは、その指定された任意点の前記検査画像の座標を、
前記平面画像の座標に変換するための係数と前記平面に
おける位置に変換するための係数とを用いて、前記検査
対象構造物の平面における位置に変換する変換手段と、
を有することを特徴とするこの発明による検査点位置取
得装置は、検査対象構造物の平面に含まれた検査領域の
任意点の平面における位置を求めるために、上記説明し
た検査点位置取得方法を具現化した装置であり、撮像領
域が段階的に平面全体まで大きく撮像された複数画像を
読み込み、これら画像間を射影変換係数(座標変換係
数)で対応づけることで、検査画像の範囲における任意
点の平面における位置を求めることができるようにする
ものである。
Further, the inspection point position acquiring device of the present invention is
An inspection point position acquisition device for determining a position on the plane of an arbitrary point of an inspection region included in a plane of an inspection target structure, wherein inspection image reading means for reading an inspection image obtained by imaging the inspection region, Plane image reading means for reading a plane image in which the plurality of reference points determined on the plane and the inspection region are imaged, reference point setting means for setting position data of the plurality of reference points measured in advance, and Recognizing the coordinates of the plurality of reference points in the plane image, and calculating a coefficient for converting any coordinates of the plane image to a position in the plane based on the recognized coordinates and the set position data. And converting the arbitrary coordinates of the inspection image into the coordinates of the plane image based on a plurality of common imaging points common to the inspection image and the plane image. A coordinate conversion coefficient calculating means for calculating the coefficients of the fit, when any point in the range of the test image is specified, the coordinates of the inspection image of the specified arbitrary point,
Using a coefficient for converting to the coordinates of the plane image and a coefficient for converting to a position on the plane, a conversion unit for converting the position of the inspection target structure into a plane,
The inspection point position obtaining apparatus according to the present invention is characterized by having the inspection point position obtaining method described above in order to obtain the position on the plane of an arbitrary point of the inspection area included in the plane of the inspection target structure. This is a embodied device that reads a plurality of images in which the imaging area is imaged stepwise to the entire plane and associates these images with a projective transformation coefficient (coordinate transformation coefficient), so that an arbitrary point in the range of the inspection image is obtained. To determine the position on the plane.

【0016】検査画像読込手段および平面画像読込手段
は、フィルムリーダ等のA−D変換装置からなり、それ
ぞれ検査画像および複数の基準点と検査領域とが撮像さ
れた平面画像を読み込む。
Each of the inspection image reading means and the plane image reading means comprises an A / D converter such as a film reader, and reads an inspection image and a plane image obtained by capturing a plurality of reference points and an inspection area.

【0017】また、予め測定された複数の基準点の位置
データが、オペレータによって入力され、これを基準点
設定手段がメモリに保存(設定)する。
Further, position data of a plurality of reference points measured in advance is inputted by an operator, and these are stored (set) in a memory by reference point setting means.

【0018】位置変換係数算出手段は、オペレータが平
面画像において指定した複数の基準点の座標を認識し、
その認識した座標と設定された位置データとに基づい
て、平面画像の任意の座標を平面における位置に変換す
るための係数を算出する。これは、いわゆる射影変換係
数であり、平面画像と実際の平面との関係において、次
の要素を含む変換が可能になる。即ち、拡大縮小、回
転、平行移動および奥行き方向の傾きをもった変換が可
能になる。
The position conversion coefficient calculating means recognizes the coordinates of a plurality of reference points specified by the operator on the plane image,
Based on the recognized coordinates and the set position data, a coefficient for converting arbitrary coordinates of the plane image into a position on the plane is calculated. This is a so-called projective transformation coefficient, which enables transformation including the following elements in the relationship between the plane image and the actual plane. That is, conversion with scaling, rotation, parallel movement, and inclination in the depth direction becomes possible.

【0019】座標変換係数算出手段は、オペレータによ
って検査画像読込手段および平面画像読込手段によって
それぞれ読み込まれた検査画像および平面画像に共通な
複数の共通撮像点が入力されると、これら共通撮像点の
各画像における座標に基づいて、検査画像の任意の座標
を平面画像の座標に変換するための係数、即ち座標変換
のための射影変換係数を算出する。
When a plurality of common imaging points common to the inspection image and the plane image read by the inspection image reading unit and the plane image reading unit are input by the operator, the coordinate transformation coefficient calculating unit inputs the common imaging points. Based on the coordinates in each image, a coefficient for converting any coordinates of the inspection image into the coordinates of the planar image, that is, a projective conversion coefficient for coordinate conversion is calculated.

【0020】そして、変換手段は、オペレータによって
検査画像の範囲で位置を求めたい任意点が指定されたと
きは、その指定された任意点の前記検査画像の座標を、
予め算出された検査画像の任意の座標を平面画像の座標
に変換するための係数と、平面画像の任意の座標を平面
における位置に変換するための係数とを用いて、座標変
換を段階的に行って、検査対象構造物の平面における位
置に変換する。
When the operator specifies an arbitrary point whose position is desired to be obtained within the range of the inspection image, the conversion means converts the coordinates of the specified arbitrary point in the inspection image into
Using a coefficient for converting any coordinates of the previously calculated inspection image into coordinates of the plane image and a coefficient for converting any coordinates of the plane image into a position on the plane, the coordinate conversion is performed step by step. Then, the position is converted into a position on the plane of the inspection target structure.

【0021】また、この発明の検査点位置取得装置は、
前記平面画像読込手段は、予め前記検査領域を取り囲む
ように定められた4箇所の基準点が撮像された平面画像
を読み込み、前記基準点設定手段は、予め測定された前
記4箇所の基準点の位置データを設定し、前記位置変換
係数算出手段は、前記平面画像における前記4箇所の基
準点の座標を認識し、その認識した座標と前記設定され
た位置データとに基づいて、前記平面画像の任意の座標
を前記平面における位置に変換するための8係数からな
る射影変換係数を算出することを特徴とする。
Further, the inspection point position acquisition device of the present invention
The plane image reading unit reads a plane image in which four reference points determined in advance to surround the inspection area are captured, and the reference point setting unit reads the four reference points measured in advance. Position data is set, and the position conversion coefficient calculation unit recognizes the coordinates of the four reference points in the plane image, and, based on the recognized coordinates and the set position data, calculates the position of the plane image. It is characterized in that a projective transformation coefficient consisting of eight coefficients for transforming an arbitrary coordinate into a position on the plane is calculated.

【0022】このような構成とすることで、検査領域を
取り囲むように定められた4箇所の基準点から、8係数
からなる射影変換係数を算出できるので、平面画像の任
意の座標から平面における位置への変換において、回
転、平行移動、縮尺および奥行き方向の傾きの要素を含
んだ変換が可能となる。
With this configuration, a projection transformation coefficient composed of eight coefficients can be calculated from four reference points defined so as to surround the inspection area. In the conversion to, conversion including elements of rotation, translation, scale, and inclination in the depth direction becomes possible.

【0023】また、この発明の検査点位置取得装置は、
前記座標変換係数算出手段は、8係数からなる射影変換
係数を算出することを特徴とする。
Further, the inspection point position obtaining apparatus of the present invention
The coordinate transformation coefficient calculating means calculates a projective transformation coefficient composed of eight coefficients.

【0024】座標変換係数算出手段でも、8係数からな
る射影変換係数を算出することで、検査画像の任意の座
標から平面画像の座標への変換において、回転、平行移
動、縮尺および奥行き方向の傾きの要素を含んだ変換が
可能となる。
The coordinate conversion coefficient calculating means also calculates a projective conversion coefficient consisting of eight coefficients, so that in the conversion from the arbitrary coordinates of the inspection image to the coordinates of the plane image, rotation, translation, scale and inclination in the depth direction are performed. The conversion including the element of becomes possible.

【0025】更に、この発明の検査点位置取得装置は、
前記座標変換係数算出手段は、非線形最適化手法を利用
した画像処理により、前記射影変換係数を算出すること
を特徴とする。
Further, the inspection point position acquiring apparatus of the present invention
The coordinate transformation coefficient calculating means calculates the projective transformation coefficient by image processing using a non-linear optimization technique.

【0026】非線形最適化手法を利用した画像処理によ
り射影変換係数を最適化することができる。
The projective transformation coefficient can be optimized by image processing using a non-linear optimization technique.

【0027】[0027]

【発明の実施の形態】以下、本発明の検査点位置取得方
法および装置の実施の形態を、図1ないし図7を参照し
ながら詳細に説明する。
DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS An embodiment of a method and an apparatus for acquiring inspection point positions according to the present invention will be described below in detail with reference to FIGS.

【0028】図1は、本発明の検査点位置取得装置の第
1の実施の形態の機器構成ブロック図である。
FIG. 1 is a block diagram of a device configuration of a first embodiment of an inspection point position acquiring apparatus according to the present invention.

【0029】カメラ1は、平面構造物の外観撮影を行う
ためのアナログカメラである。
The camera 1 is an analog camera for photographing the appearance of a planar structure.

【0030】フィルムリーダ2は、撮影フィルムをスキ
ャンして読み込み、平面構造物の外観画像をディジタル
化する装置である。
The film reader 2 is a device that scans and reads a photographic film and digitizes an appearance image of a planar structure.

【0031】データ処理部3は、外観画像のディジタル
データを処理、解析して、本発明の目的である構造物に
おける検査対象点の実寸位置座標(位置)を求めるコン
ピュータシステムである。
The data processing unit 3 is a computer system that processes and analyzes the digital data of the appearance image to obtain the actual position coordinates (position) of the inspection target point on the structure, which is the object of the present invention.

【0032】この実施の形態では、作業者が平面構造物
をカメラ1で写真撮影し、その撮影フィルムは現像され
てフィルムリーダ2でディジタル画像に変換され、デー
タ処理部3において処理される。
In this embodiment, an operator takes a picture of a planar structure with a camera 1, and the photographed film is developed, converted into a digital image by a film reader 2, and processed in a data processing section 3.

【0033】図2は、データ処理部3のシステム構成図
である。このデータ処理部3において、CPU31は、
各種処理を行うための演算部31aと主記憶部31bと
から構成される。CPU31には種々の周辺機器が接続
される。まずデータバスを介して、前述したフィルムリ
ーダ2および外部記憶装置32が接続される。外部記憶
装置32は、取り外し可能な記憶媒体である磁気ディス
クMにプログラム記憶領域32aと保存するデータ保存
領域32bを設定し、解析処理を実行するためのプログ
ラム、処理の過程で生成される中間データや結果データ
等をデータをそれぞれ保存する。フィルムリーダ2から
は、ディジタル化された平面構造物の外観画像データが
転送されてくる。
FIG. 2 is a system configuration diagram of the data processing unit 3. In the data processing unit 3, the CPU 31
It comprises an operation unit 31a for performing various processes and a main storage unit 31b. Various peripheral devices are connected to the CPU 31. First, the above-described film reader 2 and external storage device 32 are connected via a data bus. The external storage device 32 sets a program storage area 32a and a data storage area 32b to be stored on a magnetic disk M, which is a removable storage medium, a program for executing an analysis process, and intermediate data generated during the process. And data such as result data. From the film reader 2, digitized appearance image data of the planar structure is transferred.

【0034】また、CPU31には、ディスプレイ装置
33が接続される。ディスプレイ装置33は、3箇所の
表示領域、即ち、低解像度画像領域33a、精細画像領
域33bおよびメッセージ領域33cを備えている。C
PU31には、また、マウス34およびキーボード35
が接続され、画像上でのデータ入力および数値データ入
力が可能となっている。
A display device 33 is connected to the CPU 31. The display device 33 includes three display areas, that is, a low-resolution image area 33a, a fine image area 33b, and a message area 33c. C
The PU 31 also has a mouse 34 and a keyboard 35.
Are connected to enable data input and numerical data input on an image.

【0035】このような構成のデータ処理部3では、演
算部31aがプログラム記憶領域32aに記憶された解
析プログラムを逐次主記憶部31bに読み込ませて、実
行する。処理の過程における中間データは、必要に応じ
て主記憶部31bやデータ保存領域32bに保存され
る。ポインタを連動させるマウス34やキーボード35
をオペレータが操作して、適宜入力が行われる。ディス
プレイ装置33は、低解像度画像領域33a、精細画像
領域33bにフィルムリーダ2から転送された画像を表
示させて、データ入力を補助するとともに、処理の結果
求められた検査対象点の実寸位置座標を、メッセージ領
域33cに表示する。またメッセージ領域33cにオペ
レータへの指示メッセージを表示して、対話式に作業を
進めることもできる。
In the data processing section 3 having such a configuration, the calculation section 31a sequentially reads the analysis program stored in the program storage area 32a into the main storage section 31b and executes it. Intermediate data in the process of processing is stored in the main storage unit 31b or the data storage area 32b as needed. Mouse 34 or keyboard 35 for linking the pointer
Is operated by the operator, and an input is appropriately performed. The display device 33 displays the image transferred from the film reader 2 in the low-resolution image area 33a and the fine image area 33b, assists in data input, and displays the actual position coordinates of the inspection target point obtained as a result of the processing. , In the message area 33c. Further, an instruction message to the operator can be displayed in the message area 33c, and the operation can be performed interactively.

【0036】このようにオペレータの入力操作に伴い、
逐次データ処理を実行することにより、実寸位置座標取
得のための解析処理を実施することができる。
As described above, with the input operation of the operator,
By executing the sequential data processing, it is possible to execute the analysis processing for obtaining the actual size position coordinates.

【0037】尚、この実施の形態では、平面構造物の画
像を取得するために、アナログカメラとフィルムリーダ
とを使用するが、これら装置に代えて、ディジタルカメ
ラを利用し、画像を直接ディジタル化してもよい。
In this embodiment, an analog camera and a film reader are used to acquire an image of a planar structure. Instead of these devices, a digital camera is used to directly digitize the image. Is also good.

【0038】図3は、この実施の形態の検査点位置取得
装置を用いた検査点位置取得方法の概念図である。以下
この図を参照しながらこの発明の概念を説明する。
FIG. 3 is a conceptual diagram of an inspection point position acquisition method using the inspection point position acquisition device of this embodiment. Hereinafter, the concept of the present invention will be described with reference to FIG.

【0039】<基準点設定のための指標の設置作業>平
面構造物の表面の四隅近傍に、この図に示すような市松
模様の指標を設置する。このような市松模様の指標とし
たのは、その図柄により中心点が視認しやすいからであ
る。これら中心点は、以下に説明する実寸座標系の基準
点となる。指標の代わりに平面構造物の目地の交差部分
や、ひび割れの分岐点を4点定めて基準点としてもよ
い。
<Indicator setting work for setting reference point> A checkerboard index as shown in this figure is set near four corners of the surface of the planar structure. The reason why such a checkered pattern is used is that the center point is easily recognized by the design. These center points serve as reference points in the actual size coordinate system described below. Instead of the index, four intersections of joints of the planar structure or four branch points of cracks may be determined and used as reference points.

【0040】<基準点の実寸座標の計測および実寸座標
系の設定>そして、これら4基準点の相対位置関係を測
量する。この位置関係をもとに平面構造物の表面を座標
平面と見立てて、実寸座標系(X,Y) を設定し、各指標の
中心点R1,R2,R3,R4(基準点)の実寸座標(X
i,Yi)(i=1,2,3,4)を決定しておく。基準点の相対位置関
係は実際に平面構造物において作業者が測量してもよい
し、3角測量により各基準点の3次元座標を計測し、そ
れらを平面構造物に投影する方法を用いてもよい。ま
た、平面構造物の目地の交差部分を基準点とする場合
は、設計図から取得してもよい。
<Measurement of Actual Size Coordinates of Reference Points and Setting of Actual Size Coordinate System> Then, the relative positional relationship between these four reference points is measured. Based on the positional relationship, the surface of the planar structure is regarded as a coordinate plane, an actual size coordinate system (X, Y) is set, and the actual size coordinates of the center points R1, R2, R3, R4 (reference points) of each index are set. (X
i, Yi) (i = 1, 2, 3, 4) is determined in advance. The relative positional relationship of the reference points may be actually measured by the worker on the planar structure, or a method of measuring the three-dimensional coordinates of each reference point by triangulation and projecting them on the planar structure. Is also good. When the intersection of the joints of the planar structure is used as the reference point, the intersection may be obtained from the design drawing.

【0041】<低解像度画像の取得>次に、基準点R1
〜R4を含む平面構造物のディジタル画像(通常は平面
すべて)を取得する。先ず、カメラ1を使って、アナロ
グ撮影を行い、現像後にフィルムリーダ2を使用してデ
ィジタル画像に変換する(以下、このディジタル画像を
低解像度画像と記す。)。撮影に際しては、歪みの少な
いレンズを使用するのが好ましい。デジタルカメラを用
いる場合は、なるべく画素数の多い高解像度カメラを使
用する。後述する精細画像の取得においても同様であ
る。
<Acquisition of Low Resolution Image> Next, the reference point R1
Acquire a digital image (usually all of the plane) of a planar structure containing RR4. First, analog photography is performed using the camera 1, and after development, it is converted into a digital image using the film reader 2 (hereinafter, this digital image is referred to as a low-resolution image). In photographing, it is preferable to use a lens with less distortion. When a digital camera is used, a high-resolution camera having as many pixels as possible is used. The same applies to the acquisition of a fine image described later.

【0042】<精細画像の取得>低解像度画像を得る一
方、ひび割れ部分Aの精細なディジタル画像を取得す
る。これは従来から行われている工程であり、ひび割れ
近傍の比較的狭い領域の画像を低解像度画像よりも高い
解像度で取得する(以下、このディジタル画像を精細画
像と記す。)。この精細画像も上記低解像度画像と同様
に、フィルムリーダ2を使用して現像後の撮影フィルム
を読み取って取得し、コンピュータに保存する。
<Acquisition of a Fine Image> While obtaining a low-resolution image, a fine digital image of the cracked portion A is acquired. This is a process that has been conventionally performed, and an image in a relatively narrow area near a crack is acquired at a higher resolution than a low-resolution image (hereinafter, this digital image is referred to as a fine image). As in the case of the low resolution image, the fine image is read and acquired by using the film reader 2 to read the developed photographic film, and stored in the computer.

【0043】<基準点の低解像度画像座標の計測>次
に、基準点R1〜4の低解像度画像の範囲における座標
(xi,yi)(i=1,2,3,4)を取得する。低解像度画像の範囲
における座標はピクセル座標で表現される(以下、低解
像度画像座標系と記す)。
<Measurement of low-resolution image coordinates of reference point> Next, the coordinates of the reference points R1 to R4 in the range of the low-resolution image will be described.
(xi, yi) (i = 1, 2, 3, 4) is obtained. The coordinates in the range of the low-resolution image are expressed by pixel coordinates (hereinafter, referred to as a low-resolution image coordinate system).

【0044】<第1射影変換係数h の算出>次に、基準
点の実寸座標および低解像度画像座標から第1射影変換
係数h を算出する。
<Calculation of First Projection Conversion Coefficient h> Next, the first projection conversion coefficient h is calculated from the actual size coordinates of the reference point and the low resolution image coordinates.

【0045】第1射影変換係数h は、平面構造物が凹凸
のない理想平面で、カメラ1をピンホールカメラと近似
的にみなしたときの座標変換係数であり、低解像度画像
座標系(x,y) と実寸座標系(X,Y) とを数式1で表したと
きの、係数h1,h2,...,h8からなる集合である。
The first projective transformation coefficient h is a coordinate transformation coefficient when the camera 1 is approximately regarded as a pinhole camera on an ideal plane having no irregularities, and has a low-resolution image coordinate system (x, y) and the actual coordinate system (X, Y) are represented by Equation 1, and are a set of coefficients h1, h2,..., h8.

【0046】[0046]

【数1】X=(h1x+h2y+h3)/(h7x+h8y+1) Y=(h4x+h5y+h6)/(h7x+h8y+1) 第1射影変換係数h は基準点R1〜4の実寸座標(xi,y
i)(i=1,2,3,4) および低解像度画像座標(xxi,yyi)(i=
1,2,3,4) を数式1に代入して8元1次連立方程式を立
て、これを解くことにより得ることができる。
X = (h1x + h2y + h3) / (h7x + h8y + 1) Y = (h4x + h5y + h6) / (h7x + h8y + 1) The first projective transformation coefficient h is the reference point R1-4. Actual coordinates (xi, y
i) (i = 1,2,3,4) and low resolution image coordinates (xxi, yyi) (i =
1,2,3,4) is substituted into Equation 1 to form an eight-dimensional linear simultaneous equation, which can be obtained by solving the equation.

【0047】<第2射影変換係数H の算出>上述の第1
射影変換係数h と同じように、精細画像における座標系
(xx,yy) (以下、精細画像座標系と記す。)から低解像
画像度画像座標系(x,y) への座標変換も、射影変換係数
を用いて表すことができる。
<Calculation of Second Projection Transform Coefficient H>
As with the projective transformation coefficient h, the coordinate system
The coordinate transformation from (xx, yy) (hereinafter, referred to as a fine image coordinate system) to the low-resolution image quality image coordinate system (x, y) can also be represented using a projective transformation coefficient.

【0048】すなわち、第2射影変換係数H(H1,H2,...,
H8)は数式2の変換を与える射影変換係数である。
That is, the second projective transformation coefficient H (H1, H2,...,
H8) is a projective transformation coefficient that gives the transformation of Equation 2.

【0049】[0049]

【数2】x=(H1xx+H2yy+H3)/(H7xx+H8yy+1) y=(H4xx+H5yy+H6)/(H7xx+H8yy+1) 第2射影変換係数H も理論的には第1射影変換係数h を
求めた方法のように、基準となる4点の座標情報を用い
て求めることができるが、精細画像には必ずしも4箇所
の指標が含まれるとは限らないので、この実施の形態で
は後述するLM法(Levenberg-Marquardt 法)という非
線形最適化手法を用いて第2射影変換係数H を求める。
X = (H1xx + H2yy + H3) / (H7xx + H8yy + 1) y = (H4xx + H5yy + H6) / (H7xx + H8yy + 1) The second projective transformation coefficient H is also theoretically As in the method of calculating the projective transformation coefficient h, it can be obtained by using the coordinate information of four reference points. However, since a fine image does not always include four indices, In the embodiment, the second projective transformation coefficient H is obtained by using a non-linear optimization technique called an LM method (Levenberg-Marquardt method) described later.

【0050】<射影変換係数の合成>既に求められた第
1射影変換係数h 、および第2射影変換係数H とを合成
して、精細画像座標系(xx,yy) から実寸座標系(X,Y) へ
の合成射影変換係数HHを算出する。これらの射影変換係
数は、いずれも8係数から構成されるので、これらを用
いて、回転、平行移動、縮尺および奥行き方向の傾きの
要素を含んだ変換を行うことができる。
<Synthesis of Projective Transformation Coefficients> The first projective transformation coefficient h and the second projective transformation coefficient H, which have already been obtained, are combined to convert the fine image coordinate system (xx, yy) to the actual size coordinate system (X, yy). Calculate the composite projective transformation coefficient HH to Y). Since all of these projective transformation coefficients are composed of eight coefficients, they can be used to perform transformation including the elements of rotation, translation, scale, and inclination in the depth direction.

【0051】<オブジェクトの指定>この実施の形態で
オブジェクトとは実寸座標が求められる検査対象点であ
り、オペレータが精細画像の範囲内で指定する。
<Designation of Object> In this embodiment, an object is an inspection target point at which actual size coordinates are determined, and is designated by an operator within the range of a fine image.

【0052】通常この段階の前に、ひび割れ写真解析な
どの他の技術を用いて解析点検が行われ、ある部分の構
造物における実寸座標が必要となった場合に、これ以降
の段階へと進む。ここで指定されたオブジェクトPobjの
精細画像座標をPobj(xxobj,yyobj) とする。
Normally, before this stage, an analysis and inspection is performed by using another technique such as crack photo analysis, and when the actual coordinates of a certain portion of the structure are required, the process proceeds to the subsequent stages. . The fine image coordinates of the object Pobj specified here are Pobj (xxobj, yyobj).

【0053】<オブジェクトの実寸座標の算出>最後
に、オブジェクトPobjの精細画像座標(xxobj,yyobj) を
射影変換係数HHで座標変換して、オブジェクトPobjの実
寸座標(Xobj,Yobj) 、つまり平面構造物における位置を
算出する。
<Calculation of Actual Size Coordinates of Object> Finally, the fine image coordinates (xxobj, yyobj) of the object Pobj are coordinate-transformed by the projection transformation coefficient HH, and the actual size coordinates (Xobj, Yobj) of the object Pobj, that is, the planar structure Calculate the position on the object.

【0054】以上この発明の概念を説明したが、この発
明は、射影変換係数を用いて指定されたオブジェクトの
実寸座標(位置)を段階的に求めるものである。ただ
し、このように2段階に限定するものではなく、3段階
以上の変換を行うことでも、同様に実寸座標を求めるこ
とができる。
The concept of the present invention has been described above. In the present invention, the actual size coordinates (position) of an object specified by using a projective transformation coefficient are obtained stepwise. However, the present invention is not limited to the two stages, and the actual size coordinates can be similarly obtained by performing three or more stages of conversion.

【0055】次に、この発明の検査点位置取得装置の第
1の実施の形態における、検査対象点の実寸座標取得動
作を説明する。
Next, the operation of obtaining the actual coordinates of the inspection target point in the first embodiment of the inspection point position obtaining apparatus of the present invention will be described.

【0056】図4は、第1射影変換係数h を求めるため
データ処理部3において実行される解析プログラムのフ
ローチャートである。
FIG. 4 is a flowchart of an analysis program executed in the data processing section 3 to obtain the first projective transformation coefficient h.

【0057】先ず、予め実測された基準点R1〜4の実
寸座標(Xi,Yi)(i=1,2,3,4)をキーボード35から数値入
力する(S41)。次にディスプレイ装置33の低解像
度画像表示部33aに表示された画像上で、マウス34
によりポインタを移動させ、クリックして基準点R1〜
4を指定する(S42)。演算部31aにて基準点R1
〜4の低解像度画像座標(xi,yi) (i=1,2,3,4) を算出す
る(S43)。8元1次連立方程式解法に基づく演算を
演算部31aにて実行し、第1射影変換係数hを算出
し、これを主記憶部31bおよび外部記憶装置32に装
着された磁気ディスクMのデータ保存領域32bに保存
する(S44)。すなわちデータ処理部3において、4
点の基準点の実寸座標と低解像度画像座標との関係を演
算処理することで第1射影変換係数h を求めることがで
きる。
First, the actual coordinates (Xi, Yi) (i = 1, 2, 3, 4) of the previously measured reference points R1 to R4 are numerically input from the keyboard 35 (S41). Next, on the image displayed on the low resolution image display section 33a of the display device 33, the mouse 34
The pointer is moved by clicking and the reference points R1 to R1 are clicked.
4 is designated (S42). The reference point R1 is calculated by the arithmetic unit 31a.
The low-resolution image coordinates (xi, yi) (i = 1, 2, 3, 4) are calculated (S43). The arithmetic unit 31a executes an operation based on the solution of the eight-dimensional linear simultaneous equation to calculate a first projective transformation coefficient h, and stores the data in the magnetic disk M mounted on the main storage unit 31b and the external storage device 32. It is stored in the area 32b (S44). That is, in the data processing unit 3, 4
The first projective transformation coefficient h can be obtained by calculating the relationship between the actual size coordinates of the reference point of the point and the low-resolution image coordinates.

【0058】次に、第2射影変換係数H を求めるための
解析プログラムについて説明するが、この解析プログラ
ムの実行に先立って、オペレータは、ディスプレイ装置
33の低解像度画像領域33aおよび精細画像領域33
bに表示された低解像度画像および精細画像を見て、両
画像に含まれた特徴点を2点任意に選択する。
Next, an analysis program for obtaining the second projective transformation coefficient H will be described. Prior to execution of this analysis program, the operator operates the low-resolution image area 33a and the fine image area 33 of the display device 33.
Looking at the low-resolution image and the fine image displayed in b, two feature points included in both images are arbitrarily selected.

【0059】図5は、特徴点P1,P2 を両画像において指
定するときのディスプレイ装置33に表示される画像の
模式図である。この図に示すように、ひび割れの分岐点
等の特徴的な形状部分を特徴点とすればよい。
FIG. 5 is a schematic diagram of an image displayed on the display device 33 when the feature points P1 and P2 are designated in both images. As shown in this figure, a characteristic shape portion such as a branch point of a crack may be set as a characteristic point.

【0060】図6は、第2射影変換係数H を求めるため
データ処理部3において実行される解析プログラムのフ
ローチャートである。
FIG. 6 is a flowchart of an analysis program executed in the data processing unit 3 to obtain the second projective transformation coefficient H.

【0061】オペレータは先ず、マウス34によってポ
インタを操作して、図5に示すように低解像度画像上
で、特徴点P1,P2 を指定する(S61)。これらの特徴
点は、最適化計算の初期値を単に与えるだけのものであ
るから、形状を正確にとらえなくてもよい。次にデータ
処理部3の演算部31aが、P1,P2 の低解像度画像座標
(xp1,yp1),(xp2,yp2) を算出する(S62)。そして、
精細画像においてもオペレータは特徴点P1,P2 を指定す
る(S63)。演算部31aが、P1,P2 の精細画像座標
(xxp1,yyp1),(xxp2,yyp2) を算出する(S64)。
First, the operator operates the pointer with the mouse 34 to designate feature points P1 and P2 on the low-resolution image as shown in FIG. 5 (S61). Since these feature points merely provide an initial value of the optimization calculation, it is not necessary to accurately grasp the shape. Next, the calculation unit 31a of the data processing unit 3 calculates the low-resolution image coordinates of P1 and P2.
(xp1, yp1), (xp2, yp2) are calculated (S62). And
The operator also specifies feature points P1 and P2 in the fine image (S63). The calculating unit 31a calculates the fine image coordinates of P1 and P2.
(xxp1, yyp1), (xxp2, yyp2) are calculated (S64).

【0062】次に、最適化計算のために必要な初期値H0
を算出する(S65)。
Next, the initial value H0 required for the optimization calculation
Is calculated (S65).

【0063】LM法では、解析計算を実行する前に第2
射影変換係数H の初期値H0を与える必要がある。この初
期値H0にはある程度の誤差が含まれていてもよい。
In the LM method, before executing the analysis calculation, the second
It is necessary to give the initial value H0 of the projective transformation coefficient H. This initial value H0 may include some error.

【0064】この実施の形態では、精細画像座標から低
解像度画像座標への変換を近似的にヘルマート変換(射
影変換の一種であり、回転・平行移動・スケール変換の
みによる変換)とみなして先ず次の制限を与えることに
より、初期値H0の未知係数を8から4に減らす。
In this embodiment, the conversion from the fine image coordinates to the low-resolution image coordinates is approximately regarded as a Helmert transform (a type of projective transformation, and is performed only by rotation, translation, and scale transformation). , The unknown coefficient of the initial value H0 is reduced from 8 to 4.

【0065】[0065]

【数3】H1=H5,H2=-H4,H7=0,H8=0 さらに、特徴点P1,P2 の低解像度画像座標および精細画
像座標を数式2に代入して得られる4元1次連立方程式
を解く演算を、演算部31aに実行させることで、射影
変換係数H の初期値H0を算出する。
## EQU3 ## H1 = H5, H2 = -H4, H7 = 0, H8 = 0 Further, a quaternary linear system obtained by substituting the low-resolution image coordinates and the fine image coordinates of the feature points P1 and P2 into Equation 2 The arithmetic unit 31a executes an operation to solve the equation, thereby calculating an initial value H0 of the projective transformation coefficient H.

【0066】また、低解像度画像と精細画像の概略の縮
尺比Sが分かっており、かつ両者の間にほとんど回転が
ない場合、さらに次の制限を加えることができる。
When the approximate scale ratio S between the low-resolution image and the fine image is known and there is almost no rotation between the two, the following restriction can be further applied.

【0067】[0067]

【数4】H1=S,H2=0 これにより、[数2]の未知係数は4から2に減り、2
元1次連立方程式を解くだけでよい。このときオペレー
タは、ただ一つの特徴点の低解像度画像座標および精細
画像座標のみが必要となる。
H1 = S, H2 = 0 As a result, the unknown coefficient of [Equation 2] is reduced from 4 to 2, and
It is only necessary to solve the original linear simultaneous equations. At this time, the operator needs only the low-resolution image coordinates and the fine image coordinates of only one feature point.

【0068】次に、初期値H0を基にして、LM法による
非線形最適化を演算部31aに実行させて、第2射影変
換係数H の最終値を求め、このデータを保存する(S6
6)。
Next, based on the initial value H0, the arithmetic unit 31a executes the non-linear optimization by the LM method, obtains the final value of the second projective transformation coefficient H, and saves this data (S6).
6).

【0069】即ちこのステップでは、演算部31aが次
の評価関数F(H)を最小にするような最適化演算を行って
第2射影変換係数H の収束値(最終値)を求める。
That is, in this step, the calculation unit 31a performs an optimization calculation so as to minimize the next evaluation function F (H) to obtain a convergence value (final value) of the second projective transformation coefficient H.

【0070】[0070]

【数5】 ここでN は精細画像のピクセル総数、I(x,y)は低解像度
画像の座標(x,y) のピクセルの輝度、II(x,y) は精細画
像の座標(xx,yy) のピクセルの輝度である。
(Equation 5) Where N is the total number of pixels in the detail image, I (x, y) is the luminance of the pixel at coordinates (x, y) in the low-resolution image, and II (x, y) is the pixel at coordinates (xx, yy) in the detail image. Luminance.

【0071】このように仮定の射影変換係数を用いた射
影変換前後の画像において、輝度差のピクセル総和が最
小となるような非線形最適化演算を実行することより、
第2射影変換係数H を収束的に求めることができる。ま
た、座標変換をヘルマート変換とみなして第2射影変換
係数H の初期値H0を最適化演算処理に与えることで、初
期値を与えない場合に比べ最適化の誤差が少なくなり、
計算時間も短縮できる。
As described above, in the image before and after the projective transformation using the assumed projective transform coefficients, the nonlinear optimization operation is performed such that the pixel sum of the luminance difference is minimized.
The second projective transformation coefficient H can be obtained convergently. Further, by giving the initial value H0 of the second projective transformation coefficient H to the optimization calculation processing by regarding the coordinate transformation as a Helmert transformation, the error of optimization is reduced as compared with a case where no initial value is given,
Calculation time can also be reduced.

【0072】以上説明したようにして第1射影変換係数
h および第2射影変換係数H が求まり保存されれば、後
述する実寸座標取得の準備が整ったことになる。
As described above, the first projective transformation coefficient
When h and the second projective transformation coefficient H are obtained and stored, it means that preparations for acquiring actual size coordinates, which will be described later, are completed.

【0073】図7は、演算部31aにて実行される、実
寸座標の取得プログラムのフローチャートである。
FIG. 7 is a flowchart of a program for obtaining actual size coordinates, which is executed by the arithmetic unit 31a.

【0074】先ず、オペレータは精細画像でのオブジェ
クトPobjの指定をする(S71)。この指定は、精細画
像をディスプレイ装置33の精細画像領域33bに表示
させて、マウス34によりポインタを移動させてクリッ
クして行う。次に、オブジェクトPobjの精細画像座標(x
xobj,yyobj) を算出する(S72)。そして先に求めた
第2射影変換係数H を読み出して、精細画像座標(xxob
j,yyobj) を座標変換し、低解像度画像座標(xobj,yobj)
を算出する(S73)。さらに第1射影変換係数h を
用いて低解像度画像座標(xobj,yobj) を座標変換し、オ
ブジェクトPobjの実寸座標(Xobj,Yobj) を算出する(S
74)。ディスプレイ装置33のメッセージ領域33c
に、算出されたオブジェクトPobjの実寸座標(Xobj,Yob
j) を表示する(S75)。
First, the operator specifies the object Pobj in the fine image (S71). This designation is performed by displaying a fine image in the fine image area 33b of the display device 33, moving the pointer with the mouse 34, and clicking. Next, the fine image coordinates (x
xobj, yyobj) is calculated (S72). Then, the second projective transformation coefficient H determined earlier is read out, and the fine image coordinates (xxob
j, yyobj) is converted to low-resolution image coordinates (xobj, yobj)
Is calculated (S73). Further, the low-resolution image coordinates (xobj, yobj) are transformed using the first projective transformation coefficient h to calculate the actual size coordinates (Xobj, Yobj) of the object Pobj (S).
74). Message area 33c of display device 33
Is the calculated actual coordinates of the object Pobj (Xobj, Yob
j) is displayed (S75).

【0075】尚、オブジェクトの精細画像座標の指定に
ついては、別に用意したひび割れの解析装置等により求
められた種々のデータ(ひび割れの位置や、これと相対
的に求められた巾や長さをラスタ又はベクトルデータと
したもの等)から演算して、オペレータを介さずに入力
しても良い。
For the designation of the fine image coordinates of the object, various data (crack positions, widths and lengths determined relative to the positions of the cracks, and the widths and lengths determined relative thereto) obtained by a separately prepared crack analysis device or the like are used. Alternatively, it may be calculated from vector data, etc., and input without the intervention of an operator.

【0076】このように、指定されたオブジェクトの精
細画像座標に対し、射影変換(座標変換)を2段階に実
行することで、オブジェクトの実寸座標を求めることが
できる。そして、画像データや射影変換係数データは保
存されているので、実寸座標の算出は、異なるオブジェ
クトに対して何回でも行うことができる。尚、図3に示
したように、合成した射影変換係数HHを予め求めてお
き、オブジェクトPobjの精細画像座標(xxobj,yyobj) か
ら、直接、実寸座標(Xobj,Yobj) を算出することもでき
る。
As described above, the projective transformation (coordinate transformation) is performed on the fine image coordinates of the designated object in two steps, whereby the actual size coordinates of the object can be obtained. Since the image data and the projection transformation coefficient data are stored, the actual coordinates can be calculated for different objects any number of times. As shown in FIG. 3, it is also possible to previously obtain the synthesized projective transformation coefficient HH and directly calculate the actual size coordinates (Xobj, Yobj) from the fine image coordinates (xxobj, yyobj) of the object Pobj. .

【0077】<精細画像と低解像度画像の解像度比の制
限>尚、精細画像と低解像度画像の解像度比が10以内
なら、LM法による第2射影変換係数H(H1,H2,...,H8)
の収束が精度よく行われるが、解像度比が10以上のと
きは、求められた第2射影変換係数H には比較的大きい
誤差が含まれることが経験により分かっている。誤差を
少なくするには、中間の解像度を有する画像(中解像度
画像)を取得して、精細画像と中解像度画像、および中
解像度画像と低解像度画像の、それぞれを対応させる射
影変換係数を精度よく求め、これら射影変換係数による
座標変換を行うことで、解像度比が10以内の場合のよ
うに精度よい実寸座標を求めることができる。また、コ
ンクリート構造物などの表面の特徴の乏しいものを検査
対象平面とする場合は、最適化の初期値計算前に、ヒス
トグラムの正規化等の画像処理を行うことで、第2射影
変換係数H に含まれる誤差を小さくすることができる。
<Restriction on the resolution ratio between the fine image and the low-resolution image> If the resolution ratio between the fine image and the low-resolution image is within 10, the second projection conversion coefficient H (H1, H2,. H8)
It is known from experience that when the resolution ratio is 10 or more, the obtained second projective transformation coefficient H 2 contains a relatively large error. In order to reduce the error, an image having an intermediate resolution (medium resolution image) is obtained, and the projective transformation coefficients corresponding to the fine image, the medium resolution image, and the medium resolution image and the low resolution image are accurately determined. Then, by performing coordinate conversion using these projective conversion coefficients, it is possible to obtain accurate actual size coordinates as in the case where the resolution ratio is within 10. In addition, when a plane having a poor surface feature such as a concrete structure is used as the inspection target plane, image processing such as normalization of a histogram is performed before calculation of the initial value of optimization, so that the second projective transformation coefficient H Can be reduced.

【0078】また、本実施の形態では、回転、平行移
動、縮尺および奥行き方向の傾きの要素を含んだ座標変
換を行うために、8係数からなる射影変換係数を用いた
が、本発明による検査点位置取得装置は、この実施の形
態に限るものではない。
Further, in the present embodiment, in order to perform coordinate transformation including elements of rotation, translation, scale and inclination in the depth direction, a projective transformation coefficient composed of eight coefficients is used. The point position acquisition device is not limited to this embodiment.

【0079】例えば、実寸座標系と低解像度画像座標系
との座標変換をヘルマート変換であらわすことで実用上
十分であれば、第1射影変換係数h は4係数から構成さ
れることとなり、この実施の形態で4箇所設けた指標を
2箇所に減らすこともできる。
For example, if it is practically sufficient to express the coordinate transformation between the actual size coordinate system and the low-resolution image coordinate system by the Helmert transformation, the first projective transformation coefficient h is composed of four coefficients. It is also possible to reduce the index provided at four places in the form of the above to two places.

【0080】また、精細画像座標系と低解像度画像座標
系との座標変換も、ヘルマート変換のみで十分であれ
ば、説明したLM法による最適化を省いてもよい。
If the coordinate transformation between the fine image coordinate system and the low-resolution image coordinate system is sufficient only by the Helmert transformation, the optimization by the LM method described above may be omitted.

【0081】以上説明のように、この発明の検査点位置
取得方法あるいは装置は、基準点および共通撮像点(特
徴点)から求めた係数を用いて座標変換を段階的に行う
ことにより、検査画像(狭域画像)の任意点(オブジェ
クト)から、その検査対象構造物の平面全体における位
置を求めることができる。
As described above, the inspection point position acquiring method or apparatus of the present invention performs the coordinate conversion stepwise using the coefficients obtained from the reference point and the common imaging point (feature point), thereby obtaining the inspection image. From an arbitrary point (object) of the (narrow area image), the position of the inspection target structure on the entire plane can be obtained.

【0082】また、8係数からなる射影変換係数を用い
たことにより、斜めの角度から撮像された画像(奥行き
方向の傾き要素を含む画像)からでも検査点の位置が正
確に求まる。また射影変換係数を用いたことで、フィル
ムリーダの読取り歪みも吸収でき、正確な座標変換がで
きるので実用上で優れた効果を奏する。
In addition, the use of the projective transformation coefficient composed of eight coefficients enables the position of the inspection point to be accurately obtained even from an image captured at an oblique angle (an image including a tilt element in the depth direction). In addition, the use of the projective transformation coefficient can also absorb the reading distortion of the film reader, and can perform accurate coordinate transformation.

【0083】[0083]

【発明の効果】この発明の検査点位置取得方法によれ
ば、基準点の位置データおよび共通撮像点から求めた係
数を用いて座標変換を段階的に行うことにより、指定さ
れた検査領域の任意点からその検査対象構造物の平面全
体における位置を求めることができるので、実用上で優
れた効果を奏する。特に大規模な平面構造物の点検業務
に際し、遠隔的に位置が取得できるので、現地調査の手
間が省け時間と労力を軽減できる。
According to the inspection point position acquisition method of the present invention, the coordinate transformation is performed stepwise using the position data of the reference point and the coefficient obtained from the common imaging point, so that the specified inspection area can be arbitrarily selected. Since the position of the inspection target structure in the entire plane can be obtained from the point, an excellent effect in practical use is obtained. Especially when inspecting a large-scale planar structure, the position can be obtained remotely, so that the time and labor required for on-site inspection can be reduced and the time and labor can be reduced.

【図面の簡単な説明】[Brief description of the drawings]

【図1】本発明の検査点位置取得装置の第1の実施の形
態の機器構成ブロック図である。
FIG. 1 is a block diagram of a device configuration of an inspection point position acquisition device according to a first embodiment of the present invention.

【図2】図1に示す装置のデータ処理部の構成図であ
る。
FIG. 2 is a configuration diagram of a data processing unit of the device shown in FIG.

【図3】本発明による検査点位置取得方法の概念図であ
る。
FIG. 3 is a conceptual diagram of an inspection point position acquisition method according to the present invention.

【図4】図1に示す装置における第1射影変換係数を求
めるためのプログラムのフローチャートである。
FIG. 4 is a flowchart of a program for obtaining a first projective transformation coefficient in the apparatus shown in FIG. 1;

【図5】図1に示す装置における特徴点を指定するとき
のディスプレイ表示画像の模式図である。
FIG. 5 is a schematic diagram of a display image when a feature point is designated in the apparatus shown in FIG. 1;

【図6】図1に示す装置における第2射影変換係数を求
めるためのプログラムのフローチャートである。
FIG. 6 is a flowchart of a program for obtaining a second projective transformation coefficient in the apparatus shown in FIG. 1;

【図7】図1に示す装置におけるオブジェクトの実寸座
標を求めるためのプログラムのフローチャートである。
FIG. 7 is a flowchart of a program for obtaining actual size coordinates of an object in the apparatus shown in FIG. 1;

【符号の説明】[Explanation of symbols]

1 カメラ 2 フィルムリーダ 3 データ処理部 31 CPU 31a 演算部 31b 主記憶部 32 外部記憶装置 32a プログラム記憶領域 32b データ保存領域 33 ディスプレイ装置 33a 低解像度画像領域 33b 精細画像領域 33c メッセージ領域 34 マウス 35 キーボード Reference Signs List 1 camera 2 film reader 3 data processing unit 31 CPU 31a calculation unit 31b main storage unit 32 external storage device 32a program storage area 32b data storage area 33 display device 33a low resolution image area 33b detail image area 33c message area 34 mouse 35 keyboard

───────────────────────────────────────────────────── フロントページの続き (72)発明者 尾幡 昌芳 東京都新宿区新宿4−2−18新宿光風ビル アジア航測株式会社内 (72)発明者 土居原 健 東京都新宿区新宿4−2−18新宿光風ビル アジア航測株式会社内 Fターム(参考) 2F065 AA03 BB02 BB28 CC14 DD03 EE08 FF01 FF04 FF26 FF66 FF67 JJ13 QQ03 QQ21 QQ24 QQ29 QQ32 RR02 SS13 UU05 5B057 CD20 DA03 DA07 DB02 DC05 5L096 AA11 BA03 CA02 EA26 FA09 FA69 9A001 DZ13 EE05 GG08 HH24 HH27 KK27 KK37 KZ42  ──────────────────────────────────────────────────続 き Continuing on the front page (72) Inventor Masayoshi Obata 4-2-18 Shinjuku, Shinjuku, Shinjuku-ku, Tokyo Shinjuku Kofu Building Inside Asia Air Survey Co., Ltd. (72) Inventor Ken Ken Doihara 4-2, Shinjuku, Shinjuku-ku, Tokyo 18 Shinjuku Kofu Building Asia Air Survey Co., Ltd. EE05 GG08 HH24 HH27 KK27 KK37 KZ42

Claims (5)

【特許請求の範囲】[Claims] 【請求項1】 検査対象構造物の平面に含まれた検査領
域の任意点の前記平面における位置を求める検査点位置
取得方法であって、 前記検査領域が撮像された検査画像を読み込むステップ
と、 予め前記平面に定められた複数の基準点と前記検査領域
とが撮像された平面画像を読み込むステップと、 予め測定された前記複数の基準点の位置データを設定す
るステップと、 前記平面画像における前記複数の基準点の座標を認識
し、その認識した座標と前記設定された位置データとに
基づいて、前記平面画像の任意の座標を前記平面におけ
る位置に変換するための係数を算出するステップと、 前記検査画像および前記平面画像に共通な複数の共通撮
像点に基づいて、前記検査画像の任意の座標を前記平面
画像の座標に変換するための係数を算出するステップ
と、 前記検査画像の範囲で任意点が指定されたときは、その
指定された任意点の前記検査画像の座標を、前記平面画
像の座標に変換するための係数と前記平面における位置
に変換するための係数とを用いて、前記検査対象構造物
の平面における位置に変換するステップと、を有するこ
とを特徴とする検査点位置取得方法。
1. An inspection point position obtaining method for obtaining a position on an arbitrary point of an inspection area included in a plane of an inspection target structure on the plane, the method comprising: reading an inspection image in which the inspection area is imaged; Reading a planar image in which the plurality of reference points and the inspection region determined in advance on the plane are captured; setting the position data of the plurality of reference points measured in advance; and Recognizing the coordinates of a plurality of reference points, based on the recognized coordinates and the set position data, calculating a coefficient for converting any coordinate of the planar image to a position on the plane, A coefficient for converting arbitrary coordinates of the inspection image into coordinates of the plane image is calculated based on a plurality of common imaging points common to the inspection image and the plane image. When an arbitrary point is specified in the range of the inspection image, the coordinates of the inspection image at the specified arbitrary point are converted into a coefficient for converting into coordinates of the plane image and a position on the plane. Converting the position of the inspection target structure into a position on a plane using a coefficient for performing the inspection.
【請求項2】 検査対象構造物の平面に含まれた検査領
域の任意点の前記平面における位置を求める検査点位置
取得装置であって、 前記検査領域が撮像された検査画像を読み込む検査画像
読込手段と、 予め前記平面に定められた複数の基準点と前記検査領域
とが撮像された平面画像を読み込む平面画像読込手段
と、 予め測定された前記複数の基準点の位置データを設定す
る基準点設定手段と、 前記平面画像における前記複数の基準点の座標を認識
し、その認識した座標と前記設定された位置データとに
基づいて、前記平面画像の任意の座標を前記平面におけ
る位置に変換するための係数を算出する位置変換係数算
出手段と、 前記検査画像および前記平面画像に共通な複数の共通撮
像点に基づいて、前記検査画像の任意の座標を前記平面
画像の座標に変換するための係数を算出する座標変換係
数算出手段と、 前記検査画像の範囲で任意点が指定されたときは、その
指定された任意点の前記検査画像の座標を、前記平面画
像の座標に変換するための係数と前記平面における位置
に変換するための係数とを用いて、前記検査対象構造物
の平面における位置に変換する変換手段と、を有するこ
とを特徴とする検査点位置取得装置。
2. An inspection point position acquisition device for determining a position on an arbitrary point of an inspection area included in a plane of an inspection target structure on the plane, wherein an inspection image reading apparatus reads an inspection image obtained by capturing the inspection area. Means, a plurality of reference points predetermined on the plane and a plane image reading means for reading a plane image of the inspection area, and a reference point for setting position data of the plurality of reference points measured in advance. Setting means for recognizing coordinates of the plurality of reference points in the plane image, and converting arbitrary coordinates of the plane image to positions in the plane based on the recognized coordinates and the set position data. Position conversion coefficient calculating means for calculating a coefficient for calculating the coordinates of the inspection image and the plane image based on a plurality of common imaging points common to the plane image. Coordinate conversion coefficient calculating means for calculating a coefficient for converting into coordinates, and when an arbitrary point is specified in the range of the inspection image, the coordinates of the inspection image at the specified arbitrary point are converted to the coordinates of the plane image. Using a coefficient for converting to coordinates and a coefficient for converting to a position on the plane, converting means for converting to a position on the plane of the structure to be inspected; apparatus.
【請求項3】 前記平面画像読込手段は、予め前記検査
領域を取り囲むように定められた4箇所の基準点が撮像
された平面画像を読み込み、 前記基準点設定手段は、予め測定された前記4箇所の基
準点の位置データを設定し、 前記位置変換係数算出手段は、前記平面画像における前
記4箇所の基準点の座標を認識し、その認識した座標と
前記設定された位置データとに基づいて、前記平面画像
の任意の座標を前記平面における位置に変換するための
8係数からなる射影変換係数を算出することを特徴とす
る請求項2記載の検査点位置取得装置。
3. The plane image reading means reads a plane image in which four reference points predetermined so as to surround the inspection area are imaged, and the reference point setting means reads the four points measured in advance. Setting the position data of the reference point of the position, the position conversion coefficient calculating means recognizes the coordinates of the four reference points in the plane image, and based on the recognized coordinates and the set position data. 3. The inspection point position acquiring apparatus according to claim 2, wherein a projective transformation coefficient composed of eight coefficients for transforming arbitrary coordinates of the plane image into a position on the plane is calculated.
【請求項4】 前記座標変換係数算出手段は、8係数か
らなる射影変換係数を算出することを特徴とする請求項
2または請求項3記載の検査点位置取得装置。
4. The inspection point position acquiring apparatus according to claim 2, wherein said coordinate transformation coefficient calculating means calculates a projective transformation coefficient composed of eight coefficients.
【請求項5】 前記座標変換係数算出手段は、非線形最
適化手法を利用した画像処理により、前記射影変換係数
を算出することを特徴とする請求項4項記載の検査点位
置取得装置。
5. The inspection point position acquiring apparatus according to claim 4, wherein said coordinate transformation coefficient calculating means calculates said projective transformation coefficient by image processing using a non-linear optimization technique.
JP29729898A 1998-10-19 1998-10-19 Inspection point position acquisition method and apparatus Expired - Fee Related JP4170469B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP29729898A JP4170469B2 (en) 1998-10-19 1998-10-19 Inspection point position acquisition method and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP29729898A JP4170469B2 (en) 1998-10-19 1998-10-19 Inspection point position acquisition method and apparatus

Publications (2)

Publication Number Publication Date
JP2000121358A true JP2000121358A (en) 2000-04-28
JP4170469B2 JP4170469B2 (en) 2008-10-22

Family

ID=17844708

Family Applications (1)

Application Number Title Priority Date Filing Date
JP29729898A Expired - Fee Related JP4170469B2 (en) 1998-10-19 1998-10-19 Inspection point position acquisition method and apparatus

Country Status (1)

Country Link
JP (1) JP4170469B2 (en)

Cited By (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002328096A (en) * 2001-04-27 2002-11-15 Tokyo Electric Power Co Inc:The Program, method, and system for detecting crack defect generated on structure
JP2007164258A (en) * 2005-12-09 2007-06-28 Kurabo Ind Ltd Image synthesizing device and method
JP2012504222A (en) * 2008-07-04 2012-02-16 シック アイヴィピー エービー calibration
JP2014035199A (en) * 2012-08-07 2014-02-24 Kurabo Ind Ltd Precise measurement method
KR101370523B1 (en) 2012-10-19 2014-03-27 주식회사 지오시스템 Method for measurement and accuracy control of ship and marine structures
CN103808271A (en) * 2012-11-15 2014-05-21 张文镒 Method and device for remotely measuring cracks
JP2016019053A (en) * 2014-07-04 2016-02-01 Kddi株式会社 Video transfer system, terminal, program, and method that are for displaying imaging area frame superposed on wide-area picture
JP2016205901A (en) * 2015-04-17 2016-12-08 株式会社フジタ Position recognition apparatus for inspected object
WO2018037689A1 (en) * 2016-08-22 2018-03-01 富士フイルム株式会社 Image processing device and image processing method
JP2020064566A (en) * 2018-10-19 2020-04-23 ファナック株式会社 Image processing apparatus using a plurality of imaging conditions
CN112767415A (en) * 2021-01-13 2021-05-07 深圳瀚维智能医疗科技有限公司 Chest scanning area automatic determination method, device, equipment and storage medium
CN115719492A (en) * 2022-11-29 2023-02-28 中国测绘科学研究院 Method, device and equipment for identifying width characteristics of planar elements and readable storage medium
JP7450502B2 (en) 2020-09-15 2024-03-15 株式会社安藤・間 Target spacer and image analysis method

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002328096A (en) * 2001-04-27 2002-11-15 Tokyo Electric Power Co Inc:The Program, method, and system for detecting crack defect generated on structure
JP2007164258A (en) * 2005-12-09 2007-06-28 Kurabo Ind Ltd Image synthesizing device and method
JP2012504222A (en) * 2008-07-04 2012-02-16 シック アイヴィピー エービー calibration
JP2014035199A (en) * 2012-08-07 2014-02-24 Kurabo Ind Ltd Precise measurement method
KR101370523B1 (en) 2012-10-19 2014-03-27 주식회사 지오시스템 Method for measurement and accuracy control of ship and marine structures
CN103808271A (en) * 2012-11-15 2014-05-21 张文镒 Method and device for remotely measuring cracks
JP2014098683A (en) * 2012-11-15 2014-05-29 National Applied Research Laboratories Method for remotely measuring crack
JP2016019053A (en) * 2014-07-04 2016-02-01 Kddi株式会社 Video transfer system, terminal, program, and method that are for displaying imaging area frame superposed on wide-area picture
JP2016205901A (en) * 2015-04-17 2016-12-08 株式会社フジタ Position recognition apparatus for inspected object
WO2018037689A1 (en) * 2016-08-22 2018-03-01 富士フイルム株式会社 Image processing device and image processing method
JPWO2018037689A1 (en) * 2016-08-22 2019-03-07 富士フイルム株式会社 Image processing device
JP2020064566A (en) * 2018-10-19 2020-04-23 ファナック株式会社 Image processing apparatus using a plurality of imaging conditions
JP7376227B2 (en) 2018-10-19 2023-11-08 ファナック株式会社 Image processing device using multiple imaging conditions
JP7450502B2 (en) 2020-09-15 2024-03-15 株式会社安藤・間 Target spacer and image analysis method
CN112767415A (en) * 2021-01-13 2021-05-07 深圳瀚维智能医疗科技有限公司 Chest scanning area automatic determination method, device, equipment and storage medium
CN115719492A (en) * 2022-11-29 2023-02-28 中国测绘科学研究院 Method, device and equipment for identifying width characteristics of planar elements and readable storage medium
CN115719492B (en) * 2022-11-29 2023-08-11 中国测绘科学研究院 Method, device and equipment for recognizing width characteristics of planar elements and readable storage medium

Also Published As

Publication number Publication date
JP4170469B2 (en) 2008-10-22

Similar Documents

Publication Publication Date Title
JP4232064B2 (en) Crack evaluation system using image processing
US7257249B2 (en) Discrete linear space sampling method and apparatus for generating digital 3D models
JP2005215917A (en) Working drawing creation support method and replacement model creation method
KR20130138247A (en) Rapid 3d modeling
JP4170469B2 (en) Inspection point position acquisition method and apparatus
JP6823719B2 (en) Image composition method, image composition device, and recording medium
JP2015228215A (en) Positional information processing method
JPWO2019021719A1 (en) Damage data editing device, damage data editing method, and program
WO2020136523A1 (en) System and method for the recognition of geometric shapes
JP3282516B2 (en) 3D image measurement method
JPH10122819A (en) Method and device for calibration
JP3251557B2 (en) Wall image / damage / repair information registration management system and method
JP7329951B2 (en) Image processing device and its control method
Habib et al. 3-d Modelling of Historical Sites using low-cost Digital Cameras
JP5218849B2 (en) Defect inspection equipment
JPH06195472A (en) Image-processing system
Ahmad et al. Photogrammetric capabilities of the Kodak DC40, DCS420 and DCS460 digital cameras
JP3924576B2 (en) Three-dimensional measurement method and apparatus by photogrammetry
JP2005174151A (en) Three-dimensional image display device and method
JPH01291380A (en) Analyzing method for crack image
JPH076777B2 (en) Pattern contour detection method and length measuring apparatus using this method
Adami et al. Digital Representations and analysis of deformations induced in map supporting materials
JP4796295B2 (en) Camera angle change detection method, apparatus and program, and image processing method, equipment monitoring method, surveying method, and stereo camera setting method using the same
JP3615509B2 (en) Wall image / damage / repair information registration management system
JPH10214349A (en) Method and instrument for measuring photographic parameter and recording medium

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20051005

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: 20080729

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

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20080807

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110815

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110815

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120815

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120815

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130815

Year of fee payment: 5

LAPS Cancellation because of no payment of annual fees