JP2020201713A - 画像選択方法、画像選択プログラム、および画像選択装置 - Google Patents

画像選択方法、画像選択プログラム、および画像選択装置 Download PDF

Info

Publication number
JP2020201713A
JP2020201713A JP2019108130A JP2019108130A JP2020201713A JP 2020201713 A JP2020201713 A JP 2020201713A JP 2019108130 A JP2019108130 A JP 2019108130A JP 2019108130 A JP2019108130 A JP 2019108130A JP 2020201713 A JP2020201713 A JP 2020201713A
Authority
JP
Japan
Prior art keywords
image
scale
crack
captured
blur
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
JP2019108130A
Other languages
English (en)
Other versions
JP7211271B2 (ja
Inventor
学 中尾
Manabu Nakao
学 中尾
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.)
Fujitsu Ltd
Original Assignee
Fujitsu 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 Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2019108130A priority Critical patent/JP7211271B2/ja
Publication of JP2020201713A publication Critical patent/JP2020201713A/ja
Application granted granted Critical
Publication of JP7211271B2 publication Critical patent/JP7211271B2/ja
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

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

Abstract

【課題】ひび割れの計測を容易にすることができる。【解決手段】コンピュータが、クラックスケールの目盛り線を含む第1撮影画像のボケ量を算出し、第1撮影画像にボケがないと仮定した場合の画像であるクラックスケール模擬画像にボケ量に応じたボケを印加して、第1撮影画像に対応するクラックスケール画像を生成し、複数の第2撮影画像に含まれる計測対象物のひび割れと、クラックスケール画像内の各目盛り線との類似度を算出し、類似度に基づいて、複数の第2撮影画像のうちの計測対象の撮影画像を選択する処理を実行する。【選択図】図12

Description

本発明は、画像選択方法、画像選択プログラム、および画像選択装置に関する。
橋などのコンクリート構造物のひび割れを、ひび割れ周辺部分を撮影して得られた撮影画像を用いて点検する技術が用いられている。例えば、ひび割れにスケールを重ねた画像を表示して、ひび割れの幅を計測する処理が行われる。
関連する技術として、複数の部分変換画像を合成し、検査領域画像を生成し、その検査領域画像を出力する技術が提案されている(例えば、特許文献1を参照)。
また、関連する技術として、線状の被写体を撮像して得られた画像を入力し、線状の被写体の幅方向の濃度分布に対応する確率分布関数に基づいて、被写体の幅を示す指標を生成する技術が提案されている(例えば、特許文献2を参照)。
また、関連する技術として、構造物の画像におけるひび割れの位置及び方向に応じた位置及び方向に、生成したスケール画像を重ね合わせて表示する技術が提案されている(例えば、特許文献3を参照)。
また、関連する技術として、撮像対象物を撮像することによって得られた撮像画像を取得し、クラックスケール画像が撮像画像上に重ねて表示されるように、合成画像データを生成する技術が提案されている(例えば、特許文献4を参照)。
特開2016−218762号公報 国際公開第2017/187966号 国際公開第2017/122641号 特開2009−85785号公報
ひび割れ周辺を撮影して得られた撮影原画像は、調査範囲と比べて狭く、ひび割れの全体を把握することができない可能性がある。その場合、ひび割れの全体を把握して、計測対象のひび割れ損傷を選択するためには、ユーザが撮影原画像を複数枚見る必要があり、時間がかかる。
そこで、複数の撮影原画像を繋いで合成したオルソ画像を、ひび割れの計測に用いることが考えられる。合成したオルソ画像を用いると、ユーザは、1枚の画像でひび割れを俯瞰的に見ることができるので、計測対象のひび割れの選択および計測を効率的に行うことができる。しかし、オルソ画像内における撮影原画像の繋ぎ合わせ部分では、画像がボケて画質が低下する場合があり、撮影原画像をそのまま用いる場合と比べて、ひび割れの計測精度が低い。
よって、ひび割れ幅の計測に関しては、撮影画像を合成したオルソ画像ではなく、ひび割れが鮮明に映っている撮影原画像を使うことが望ましい。しかし、撮影原画像にも、撮影時のカメラのブレ等のため、ひび割れが大きくボケた画像が存在する可能性がある。ひび割れが映っている複数の撮影原画像の中から計測に用いる画像をユーザが選択した場合、選択に時間がかかる。
1つの側面として、本発明は、ひび割れの計測を容易にすることを目的とする。
1つの態様では、画像選択方法において、コンピュータが、クラックスケールの目盛り線を含む第1撮影画像のボケ量を算出し、前記第1撮影画像にボケがないと仮定した場合の画像であるクラックスケール模擬画像に前記ボケ量に応じたボケを印加して、前記第1撮影画像に対応するクラックスケール画像を生成し、複数の第2撮影画像に含まれる計測対象物のひび割れと、前記クラックスケール画像内の各目盛り線との類似度を算出し、前記類似度に基づいて、前記複数の第2撮影画像のうちの計測対象の撮影画像を選択する処理を実行する。
1つの側面によれば、ひび割れの計測を容易にすることができる。
実施形態のシステムの全体構成の一例を示す図である。 クラックスケールの撮影画像の輝度分布の一例を示す図である。 ボケ量の算出方法の一例を示す図である。 クラックスケール画像の生成方法の一例を示す図である。 オルソ画像の生成方法の一例を示す図である。 第2撮影画像の抽出の一例を示す図である。 第2撮影画像内のひび割れとクラックスケール画像の目盛り線との類似度算出の一例を示す。 第2撮影画像にクラックスケールを重ねた表示例を示す図である。 実施形態の処理の一例を示すフローチャートである。 図9のステップS102の処理の一例を示すフローチャートである。 図9のステップS108の処理の一例を示すフローチャートである。 情報処理装置の他の例を示す図である。 情報処理装置のハードウェア構成の一例を示す図である。
以下、図面を参照して、実施形態について説明する。図1は、実施形態のシステムの全体構成の一例を示す図である。実施形態のシステムは、情報処理装置1と表示装置2とを含む。なお、情報処理装置1内に表示装置2が含まれていてもよい。情報処理装置1は、コンピュータの一例である。情報処理装置1は、取得部11と、第1算出部12と、第1生成部13と、第2生成部14と、第3生成部15と、抽出部16と、第2算出部17と、選択部18と、表示制御部19と、記憶部20とを含む。
取得部11は、クラックスケールの目盛り線を含む撮影画像(第1撮影画像)を取得する。クラックスケールは、複数の異なる幅の目盛り線を含み、各目盛りの幅を示す数値が目盛りの近傍に表示されたスケールである。取得部11は、計測対象物(例えば、コンクリート構造物)を網羅的に撮影して得られた複数の撮影原画像(第2撮影画像)を取得する。このクラックスケールの目盛り線の幅と長さは既知であるとする。
第1算出部12は、第1撮影画像と、クラックスケールの目盛り線の模擬画像とに基づいて、第1撮影画像のボケ量を示す値σを算出する。この模擬画像は、第1撮影画像にボケがないと仮定した場合の画像である。第1算出部12は、クラックスケールの模擬画像に、複数の異なるボケ量σを印加して得られた複数のボケ印加画像Iσj(x)と、第1撮影画像との類似度をそれぞれ算出する。そして、第1算出部12は、第1撮影画像との類似度が最も高いボケ印加画像Iσj(x)に用いられたボケ量σを、クラックスケール画像の生成に用いるボケ量として決定する。
第1生成部13は、模擬画像に、第1算出部12が算出したボケ量に応じたボケを印加して、第1撮影画像に対応するクラックスケール画像を生成する。第1生成部13は、例えば、算出されたボケ量を示す値σを用いたガウス分布関数G(x,σ)と目盛り線の一次元模擬画像との畳み込み演算を行い、演算結果に基づく1次元画像を生成する。そして、第1生成部13は、生成された1次元画像を連結して、第1撮影画像に対応するクラックスケール画像を生成する。
第2生成部14は、第2撮影画像を合成し、3Dモデルを生成する。なお、第2撮影画像と第1撮影画像の撮影環境(例えば、撮影装置、撮影距離等)は同じであるとする。
第3生成部15は、生成された3Dモデルを変換して、計測対象物を含むオルソ画像を生成する。生成したオルソ画像に対して、ユーザが、ひび割れ線に沿った線および計測点を描画する。
抽出部16は、計測点が描画されたオルソ画像を参照し、オルソ画像上の計測点の座標に基づいて、複数の第2撮影画像のうち計測点を含む第2撮影画像を抽出する。
第2算出部17は、抽出された複数の第2撮影画像内に含まれる計測対象物の計測点近傍のひび割れと、クラックスケール画像の各目盛り線との類似度を、複数の第2撮影画像毎、目盛り線毎、目盛り線の位置および傾き毎に算出する。第2算出部17は、類似度として、例えば、複数の第2撮影画像内の計測点近傍のひび割れと、クラックスケール画像の各目盛り線の輝度値との相関係数を用いる。
選択部18は、算出された類似度に基づいて、複数の第2撮影画像のうちの計測対象の撮影画像と、ひび割れの計測に用いるクラックスケール画像内の目盛り線とを選択し、目盛り線の位置および傾きを選択する。選択部18は、例えば、第2算出部17が算出した類似度が最も高くなる、第2撮影画像、目盛り線、目盛り線の位置および傾きを選択する。
表示制御部19は、選択された第2撮影画像とクラックスケール画像とを重ね合わせた画像を表示装置2に表示させる。表示制御部19は、例えば、選択された第2撮影画像内のひび割れの近傍に、選択された目盛り線が、選択された位置および傾きで配置されるように、第2撮影画像とクラックスケール画像の目盛り線とを重ね合わせた画像を表示装置2に表示させる。
表示制御部19は、目盛り線を選択可能なメニューを表示装置に表示させてもよい。そして、表示制御部19は、そのメニューから、ひび割れの最も近くに配置されている目盛り線と異なる目盛り線が選択された場合、そのメニューから選択された目盛り線がひび割れの最も近くに配置されるように、クラックスケール画像の表示位置を変更する。
記憶部20は、情報処理装置1の処理に用いられる各種データを記憶する。記憶部20は、例えば、取得部11が取得した第1撮影画像および第2撮影画像、第1生成部13が生成したクラックスケール画像、第2生成部14が生成した3Dモデル、第3生成部15が生成したオルソ画像等を記憶する。
図2は、クラックスケールの撮影画像の輝度分布の一例を示す図である。図2に示す撮影画像は、取得部11が取得する第1撮影画像である。第1撮影画像には、ピントのずれ、レンズ収差、カメラの動き等により画像にボケが加わる。そのため、例えば、分解能(1画素に映る撮影対象の大きさ)が0.3mmであり、目盛り線の幅が0.1mmであったとしても、目盛り線による輝度変化が1画素内だけでなく複数画素に発生する。計測対象のひび割れを含む撮影原画像にもボケがあるため、撮影原画像と同程度のボケ量に応じたボケを印加したクラックスケール画像を用いて計測することが好ましい。しかし、クラックスケールの撮影画像は、ノイズ等が含まれる可能性がある。
そのため、情報処理装置1は、撮影画像のボケ量を算出し、算出されたボケ量に応じたボケを印加したクラックスケール画像を生成し、そのクラックスケール画像を、ひび割れ幅を計測するためのスケールとして用いる。
図3は、ボケ量の算出方法の一例を示す図である。取得部11は、クラックスケール画像の目盛り線の撮影画像(第1撮影画像)を取得する。なお、第1算出部12は、以下に説明する処理において、第1撮影画像内の複数の目盛り線のうちのいずれかの目盛り線を、ボケ量の算出に用いる。
第1算出部12は、第1撮影画像内のいずれかの目盛り線の長さと、目盛り線の画素数とに基づいて、以下の式(1)のように、分解能を算出する。なお、目盛り線の長さは既知であるとする。
(分解能[mm/画素])=(目盛り線の長さ[mm])/(目盛り線の画素数)(1)
また、第1算出部12は、第1撮影画像内の目盛り線を含む、目盛り線と垂直方向の1列の画素で形成された1次元画像(I(x))を生成する。I(x)は、この1次元画像の画素xにおける輝度を示す。以下、この1次元画像をI(x)と称することがある。
第1算出部12は、予め計測された、クラックスケール画像の背景輝度と、目盛り線部分の輝度と、式(1)を用いて算出された分解能とに基づいて、クラックスケールの目盛り線の模擬画像(I(x))を生成する。I(x)は、模擬画像の画素xにおける輝度を示す。以下、目盛り線の1次元の模擬画像をI(x)と称する。I(x)は、ボケがないと仮定した場合の、クラックスケールの目盛り線の1次元画像(I(x))に相当する。
I(x)のうち、目盛り線部分を含む画素kの輝度I(k)は、以下の式(2−1)のように表される。背景輝度は、第1撮影画像における目盛り線以外の領域の輝度である。暗部輝度は、第1撮影画像における目盛り線部分の輝度である。背景輝度、暗部輝度には、予め計測された値が用いられる。また、分解能は、式(1)で算出された分解能である。また、目盛り線幅は、I(x)に含まれる目盛り線の幅である。
Figure 2020201713
なお、目盛り線部分を含まない画素xの輝度I(x)には、背景輝度がそのまま適用される。
例えば、画素x=0において、背景輝度=200,暗部輝度=50,分解能=0.3,目盛り線幅=0.1である場合、輝度I(0)は、以下の式(2−2)のように算出される。
Figure 2020201713
第1算出部12は、生成した模擬画像(I(x))とガウス分布関数G(x,σ)との畳み込み演算を、複数のσ(j=1〜J)を用いて行う。この畳み込み演算の演算結果をIσj(x)(ボケ印加画像)とする。ガウス分布関数G(x,σ)、演算結果Iσj(x)は、以下の式(3)、(4)のように表される。式(4)における「*」は、畳み込み演算を示す記号である。
Figure 2020201713
Figure 2020201713
第1算出部12は、I(x)と、算出された複数のIσj(x)との類似度を算出する。第1算出部12は、例えば、類似度として、以下の式(5)に示す相関係数rを用いる。式(5)において、M=I(x)であり、l=Iσj(x)であるとする。また、Nは、I(x)、Iσj(x)内の全画素数であるとする。
Figure 2020201713
第1算出部12は、各Iσj(x)(j=1〜J)について上記式(5)の計算を行い、類似度(相関係数r)が最も大きい場合のσを、第1撮影画像のボケ量を示す値として決定する。
図4は、クラックスケール画像の生成方法の一例を示す図である。第1生成部13は、第1算出部12が算出した類似度が最大となるときのσを用いたガウス分布関数G(x,σ)と、クラックスケールの一次元模擬画像I(x)との畳み込み演算を行い、演算結果の1次元画像Iσj(x)を生成する。そして、第1生成部13は、生成された1次元画像Iσ(x)を複製して縦方向に連結することにより、目盛り線画像を生成する。
第1生成部13は、クラックスケールにおける全目盛り線について同様の処理を行うことにより、第1撮影画像に対応するクラックスケール画像を生成することができる。
以上のように、情報処理装置1は、第1画像のボケ量を示す値σを算出し、模擬画像I(x)と、そのσを用いたガウス分布関数G(x,σを)との畳み込み演算を行うことにより、クラックスケール画像を生成する。これにより、生成されたクラックスケール画像の輝度は、図2に示す実際のクラックスケール画像(第1撮影画像)の輝度に近い輝度となる。また、生成されたクラックスケール画像は、模擬画像を基に生成した画像であるため、クラックスケールの実際の撮影画像(第1撮影画像)よりノイズが少ない画像となる。
図5は、オルソ画像の生成方法の一例を示す図である。第2生成部14は、計測対象物(例えば、コンクリート構造物)を網羅的に撮影して得られた複数の撮影原画像(第2撮影画像)を合成し、3Dモデルを生成する。第2生成部14は、例えば、3Dモデルの生成処理に、Structure from Motion(SfM)等の技術を用いる。そして、第3生成部15は、生成された3Dモデルを変換してオルソ画像を生成する。
表示制御部19は、生成されたオルソ画像を表示装置2に表示し、ユーザのスケッチ描画を受け付ける。ユーザは、ひび割れに沿った線と計測点を、表示装置2の画面上に描画する。
図6は、第2撮影画像の抽出の一例を示す図である。抽出部16は、オルソ画像上に描画された計測点の座標を算出する。そして、抽出部16は、オルソ画像上の計測点の座標に基づいて、計測点の3次元座標を算出する。そして、抽出部16は、計測点の3次元座標に基づいて、複数の第2撮影画像の座標系における計測点の座標を算出する。例えば、複数の第2撮影画像A〜Eの座標系における計測点の座標が、以下のように表されるとする。なお、各撮影原画像の座標の単位は、画素である。
第2撮影画像Aの座標系における計測点の座標=(−200,500)
第2撮影画像Bの座標系における計測点の座標=(300,50)
第2撮影画像Cの座標系における計測点の座標=(200,400)
第2撮影画像Dの座標系における計測点の座標=(250,200)
第2撮影画像Eの座標系における計測点の座標=(1000,100)
例えば、各第2撮影画像のサイズが500×500[画素]であるとする。そして、抽出部16は、各第2撮影画像内の計測点の座標が(0,0)〜(499,499)の範囲内であれば、計測点が撮影原画像内に存在すると判定する。上記の例では、第2撮影画像A〜Eのうち、第2撮影画像B、C、Dの座標が(0,0)〜(499,499)の範囲内であるため、抽出部16は、第2撮影画像B、C、Dを抽出する。
図7は、第2撮影画像内のひび割れとクラックスケール画像の目盛り線との類似度算出の一例を示す。第2算出部17は、複数の第2撮影画像内の計測点近傍のひび割れと、クラックスケール画像の各目盛り線の類似度を、抽出された複数の第2撮影画像および各目盛り線の全組み合わせについて算出する。第2算出部17は、類似度を算出する際、目盛り線の位置および傾きの変更と類似度の算出を繰り返す。すなわち、第2算出部17は、複数の第2撮影画像内の計測点近傍におけるひび割れと、クラックスケール画像の各目盛り線の類似度を、複数の第2撮影画像毎、目盛り線毎、目盛り線の位置および傾き毎に算出する。
第2算出部17は、類似度として、例えば、相関係数を用いる。第2算出部17は、第2撮影画像内の計測点近傍のひび割れを含む所定範囲の画素xの輝度値をMとし、比較対象の目盛り線を含む所定範囲の画素xの輝度値をlとして、式(5)の相関係数rを算出し、相関係数rを類似度として用いる。
そして、選択部18は、例えば、第2算出部17が算出した類似度が最も高くなる、第2撮影画像、目盛り線、目盛り線の位置および傾きを選択する。
このように、情報処理装置1は、表示対象の第2撮影画像、目盛り線、目盛り線の位置および傾きを自動選択するため、ユーザが第2撮影画像を選択する手間と、目盛り線の位置合わせをする手間を省き、ひび割れの計測を容易にすることができる。
図8は、第2撮影画像にクラックスケールを重ねた表示例を示す図である。図8に示す表示画面31は、表示装置2の表示画面である。表示制御部19は、選択された第2撮影画像とクラックスケール画像34とを重ね合わせた画像32を、表示装置2の表示画面31に表示させる。この際、表示制御部19は、第2撮影画像内の計測点近傍のひび割れ35に、選択された目盛り線が選択された位置および傾きで配置されるように、第2撮影画像とクラックスケール画像34とを重ね合わせた画像32を表示させる。
ボケ量の算出に用いた第1撮影画像と、ひび割れを含む第2撮影画像の撮影環境は同じである。そして、クラックスケール画像は、第1撮影画像と同程度のボケ量に応じたボケが印加されているため、第2撮影画像と同程度のボケ量に応じたボケが印加されていることとなる。すなわち、情報処理装置1は、ひび割れを含む第2画像と同程度のボケが印加されたクラックスケール画像を計測用スケールとして表示させることにより、計測をし易くし、計測精度を向上することができる。
さらに、表示制御部19は、ひび割れの近傍に表示させる目盛り線を選択可能なプルダウンメニュー33を、表示装置2の表示画面31に表示させる。ユーザは、選択された目盛り線がひび割れの幅と合っていないと判断した場合、プルダウンメニュー33を用いて、ひび割れの幅に近い目盛り線を選択することができる。プルダウンメニュー33により、ひび割れの最も近くに配置されている目盛り線と異なる目盛り線が選択された場合、表示制御部19は、プルダウンメニュー33から選択された目盛り線がひび割れの最も近くに配置されるように、クラックスケール画像の表示位置を変更する。
このように、情報処理装置1は、目盛り線を選択可能なプルダウンメニュー33を表示させることにより、ひび割れの幅と目盛り線の幅が合っていない場合に、ユーザの選択によりひび割れの近傍に表示させる目盛り線を変更することができる。
図9は、実施形態の処理の一例を示すフローチャートである。取得部11は、クラックスケールの目盛り線を含む撮影画像(第1撮影画像)と、計測対象物を網羅的に撮影して得られた複数の撮影原画像(第2撮影画像)を取得する(ステップS101)。
第1算出部12は、第1撮影画像と、クラックスケールの目盛り線の模擬画像とに基づいて、第1撮影画像のボケ量を示す値σを算出する(ステップS102)。ステップS102について、詳細は後述する。
第1生成部13は、算出されたボケ量を示す値σを用いたガウス分布関数G(x,σ)と目盛り線の一次元模擬画像I(x)との畳み込み演算を行い、演算結果に基づく1次元画像Iσj(x)を生成する(ステップS103)。第1生成部13は、各目盛り線について、ステップS103の処理を行う。
第1生成部13は、生成された1次元画像Iσj(x)を連結して、第1撮影画像に対応するクラックスケール画像を生成する(ステップS104)。
第2生成部14は、計測対象物を網羅的に撮影して得られた複数の撮影原画像(第2撮影画像)を合成し、3Dモデルを生成する(ステップS105)。第3生成部15は、3Dモデルを変換して、計測対象物を含むオルソ画像を生成する(ステップS106)。生成したオルソ画像に対して、ユーザがひび割れ線に沿った線とおよび計測点を描画する。抽出部16は、計測点が描画されたオルソ画像を参照し、オルソ画像上の計測点の座標に基づいて、複数の第2撮影画像のうち計測点を含む第2撮影画像を抽出する(ステップS107)。
なお、情報処理装置1は、ステップS105〜S107の処理を、ステップS101〜S104の処理の前、またはステップS102〜S104の処理と並列に実行してもよい。
第2算出部17は、抽出された複数の第2撮影画像内に含まれる計測対象物の計測点近傍のひび割れと、クラックスケール画像の各目盛り線の類似度を、複数の第2撮影画像毎、目盛り線毎、目盛り線の位置および傾き毎に算出する(ステップS108)。ステップS108の処理について、詳細は後述する。
選択部18は、ステップS108の処理で最も類似度が高くなる、第2撮影画像、クラックスケールの目盛り線、目盛り線の位置および傾きを選択する(ステップS109)。
表示制御部19は、選択された第2撮影画像とクラックスケール画像の目盛り線とを重ね合わせた画像を表示装置2に表示させる(ステップS110)。この際、表示制御部19は、選択された第2撮影画像内のひび割れの近傍に、選択された目盛り線が、選択された位置および傾きで配置されるように、第2撮影画像とクラックスケール画像の目盛り線とを重ね合わせた画像を表示装置2に表示させる。
以上のように、情報処理装置1は、複数の第2撮影画像に含まれる計測対象物のひび割れと、クラックスケール画像内の各目盛り線との類似度を算出して、計測に用いる第2撮影画像と目盛り線を自動選択する。そのため、情報処理装置1は、ユーザが第2撮影画像と目盛り線を選択する手間を省き、ひび割れの計測を容易にすることができる。また、情報処理装置1は、クラックスケールの目盛り線を含む第1撮影画像のボケ量を算出し、そのボケ量に応じたボケを印加したクラックスケール画像を生成する。これにより、情報処理装置1は、計測対象のひび割れと目盛り線のボケを一致させて計測をし易くし、計測精度を向上することができる。
図10は、図9のステップS102の処理の一例を示すフローチャートである。第1算出部12は、上記式(1)を用いて、第1撮影画像の分解能を算出する(ステップS201)。また、第1算出部12は、第1撮影画像内の目盛り線を含む、目盛り線と垂直方向の1列の画素で形成された1次元画像(I(x))を生成する(ステップS202)。
第1算出部12は、予め計測された、クラックスケール画像の背景輝度と、目盛り線部分の輝度と、算出された分解能とに基づいて、I(x)に対応するクラックスケールの目盛り線の模擬画像(I(x))を生成する(ステップS203)。第1算出部12は、例えば、目盛り線部分を含む画素kの輝度I(k)を上記式(2−1)を用いて算出し、目盛り線部分を含まない画素xの輝度I(x)には、背景輝度をそのまま適用する。
第1算出部12は、生成した模擬画像(I(x))とガウス分布関数G(x,σ)との畳み込み演算を、複数のσを用いて行い、演算結果をIσj(x)(ボケ印加画像)とする(ステップS203)。ガウス分布関数G(x,σ)、演算結果Iσj(x)は、上記の式(3)、(4)のように表される。
第1算出部12は、I(x)と、算出された複数のIσj(x)との類似度を算出する(ステップS205)。第1算出部12は、例えば、類似度として、上記の式(5)に示す相関係数rを用いる。そして、第1算出部12は、I(x)と類似度が最も高いIσj(x)に用いられたσを、ボケ量を示す値として決定する(ステップS206)。
図11は、図9のステップS108の処理の一例を示すフローチャートである。第2算出部17は、第2撮影画像内の計測点近傍におけるひび割れと、クラックスケール画像の目盛り線の類似度を算出する(ステップS301)。算出部17は、ステップS301の処理を、抽出された複数の第2撮影画像毎、クラックスケールの目盛り線毎、目盛り線の位置および傾き毎に繰り返す。
図12は、情報処理装置1の他の例を示す図である。図12に示す情報処理装置1は、第2生成部14と、第3生成部15と、表示制御部19と、記憶部20とを含まない点で、図1に示す情報処理装置1と異なる。第2生成部14と、第3生成部15と、表示制御部19と、記憶部20は、情報処理装置1と通信可能な他の情報処理装置等に含まれているとする。他の情報処理装置は、一つであってもよいし、複数であってもよい。
この場合、例えば、取得部11は、第2生成部14および第3生成部15を有する他の情報処理装置から、生成されたオルソ画像を取得する。また、情報処理装置1は、選択部18が選択した、第2撮影画像、目盛り線、目盛り線の位置および傾きを、表示制御部19を有する他の情報処理装置に送信し、その情報処理装置が表示装置2に対して表示制御を行う。
また、情報処理装置1は、実施形態の処理に用いる各種データを、記憶部20を有する他の情報処理装置に送信し、記憶させる。
次に、情報処理装置1のハードウェア構成の一例を説明する。図13は、情報処理装置1のハードウェア構成の一例を示す図である。図13の例に示すように、情報処理装置1において、バス100に、プロセッサ111とメモリ112と補助記憶装置113と通信インタフェース114と媒体接続部115と入力装置116と出力装置117とが接続される。
プロセッサ111は、メモリ112に展開されたプログラムを実行する。実行されるプログラムには、実施形態における処理を行う画像選択プログラムが適用されてもよい。
メモリ112は、例えば、Random Access Memory(RAM)である。補助記憶装置113は、種々の情報を記憶する記憶装置であり、例えばハードディスクドライブや半導体メモリ等が適用されてもよい。補助記憶装置113に実施形態の処理を行う画像選択プログラムが記憶されていてもよい。
通信インタフェース114は、Local Area Network(LAN)、Wide Area Network(WAN)等の通信ネットワークに接続され、通信に伴うデータ変換等を行う。
媒体接続部115は、可搬型記録媒体118が接続可能なインタフェースである。可搬型記録媒体118には、光学式ディスク(例えば、Compact Disc(CD)またはDigital Versatile Disc(DVD)等)、半導体メモリ等が適用されてもよい。可搬型記録媒体118に実施形態の処理を行う画像選択プログラムが記録されていてもよい。
入力装置116は、例えば、キーボード、ポインティングデバイス等であり、ユーザからの指示及び情報等の入力を受け付ける。
出力装置117は、例えば、表示装置、プリンタ、スピーカ等であり、ユーザへの問い合わせ又は指示、及び処理結果等を出力する。図1に示す表示装置2には、図12に示す出力装置117が適用されてもよい。
図1に示す記憶部20は、メモリ112、補助記憶装置113または可搬型記録媒体118等により実現される。図1に示す取得部11、第1算出部12、第1生成部13、第2生成部14、第3生成部15、抽出部16、第2算出部17、選択部18、表示制御部19は、メモリ112に展開された画像選択プログラムをプロセッサ111が実行することにより実現される。
メモリ112、補助記憶装置113および可搬型記録媒体118は、コンピュータが読み取り可能であって非一時的な有形の記憶媒体であり、信号搬送波のような一時的な媒体ではない。
なお、情報処理装置1が図13に示す全ての構成要素を含んでいなくてもよく、一部の構成要素が省略されていてもよい。また、一部の構成要素が情報処理装置1の外部装置に存在し、情報処理装置1が外部装置に接続して、外部装置内の構成要素を利用してもよい。
本実施形態は、以上に述べた実施の形態に限定されるものではなく、本実施形態の要旨を逸脱しない範囲内で様々な変更、追加、省略が適用可能である。
1 情報処理装置
2 表示装置
11 取得部
12 第1算出部
13 第1生成部
14 第2生成部
15 第3生成部
16 抽出部
17 第2算出部
18 選択部
19 表示制御部
20 記憶部
100 バス
111 プロセッサ
112 メモリ
113 補助記憶装置
114 通信インタフェース
115 媒体接続部
116 入力装置
117 出力装置
118 可搬型記録媒体

Claims (7)

  1. コンピュータが、
    クラックスケールの目盛り線を含む第1撮影画像のボケ量を算出し、
    前記第1撮影画像にボケがないと仮定した場合の画像であるクラックスケール模擬画像に前記ボケ量に応じたボケを印加して、前記第1撮影画像に対応するクラックスケール画像を生成し、
    複数の第2撮影画像に含まれる計測対象物のひび割れと、前記クラックスケール画像内の各目盛り線との類似度を算出し、
    前記類似度に基づいて、前記複数の第2撮影画像のうちの計測対象の撮影画像を選択する
    処理を実行することを特徴とする画像選択方法。
  2. 前記コンピュータが、
    前記クラックスケール模擬画像に複数の異なるボケ量に応じたボケを印加して得られた複数のボケ印加画像と、前記第1撮影画像との類似度をそれぞれ算出し、
    前記第1撮影画像との類似度が最も高い前記ボケ印加画像に用いられた前記ボケ量を、前記クラックスケール画像の生成に用いる前記ボケ量として決定する
    処理を実行することを特徴とする請求項1に記載の画像選択方法。
  3. 前記コンピュータが、
    選択された前記第2撮影画像と前記クラックスケール画像とを重ね合わせた画像を表示装置に表示させる
    処理を実行することを特徴とする請求項1または2に記載の画像選択方法。
  4. 前記コンピュータが、
    前記計測対象物のひび割れと、前記クラックスケール画像内の各目盛り線との類似度を、前記複数の第2撮影画像毎、前記目盛り線毎、前記目盛り線の位置および傾き毎に算出し、
    前記類似度が最も高くなる、前記第2撮影画像、前記目盛り線、前記目盛り線の位置および傾きを選択し、
    選択された前記第2撮影画像内に、選択された前記目盛り線が選択された位置および傾きで配置されるように、前記第2撮影画像と前記クラックスケール画像とを重ね合わせた画像を前記表示装置に表示させる
    処理を実行することを特徴とする請求項3に記載の画像選択方法。
  5. 前記目盛り線を選択可能なメニューを前記表示装置に表示させ、
    前記メニューから、前記ひび割れの最も近くに配置されている前記目盛り線と異なる目盛り線が選択された場合、前記メニューから選択された目盛り線が前記ひび割れの最も近くに配置されるように、前記クラックスケール画像の表示位置を変更する
    処理を実行することを特徴とする請求項4に記載の画像選択方法。
  6. クラックスケールの目盛り線を含む第1撮影画像のボケ量を算出し、
    前記第1撮影画像にボケがないと仮定した場合の画像であるクラックスケール模擬画像に前記ボケ量に応じたボケを印加して、前記第1撮影画像に対応するクラックスケール画像を生成し、
    複数の第2撮影画像に含まれる計測対象物のひび割れと、前記クラックスケール画像内の各目盛り線との類似度を算出し、
    前記類似度に基づいて、前記複数の第2撮影画像のうちの計測対象の撮影画像を選択する
    処理をコンピュータに実行させるための画像選択プログラム。
  7. クラックスケールの目盛り線を含む第1撮影画像のボケ量を算出する第1算出部と、
    前記第1撮影画像にボケがないと仮定した場合の画像であるクラックスケール模擬画像に前記ボケ量に応じたボケを印加して、前記第1撮影画像に対応するクラックスケール画像を生成する生成部と、
    複数の第2撮影画像に含まれる計測対象物のひび割れと、前記クラックスケール画像内の各目盛り線との類似度を算出する第2算出部と、
    前記類似度に基づいて、前記複数の第2撮影画像のうちの計測対象の撮影画像を選択する選択部と
    を備えることを特徴とする画像選択装置。
JP2019108130A 2019-06-10 2019-06-10 画像選択方法、画像選択プログラム、および画像選択装置 Active JP7211271B2 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2019108130A JP7211271B2 (ja) 2019-06-10 2019-06-10 画像選択方法、画像選択プログラム、および画像選択装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2019108130A JP7211271B2 (ja) 2019-06-10 2019-06-10 画像選択方法、画像選択プログラム、および画像選択装置

Publications (2)

Publication Number Publication Date
JP2020201713A true JP2020201713A (ja) 2020-12-17
JP7211271B2 JP7211271B2 (ja) 2023-01-24

Family

ID=73743381

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2019108130A Active JP7211271B2 (ja) 2019-06-10 2019-06-10 画像選択方法、画像選択プログラム、および画像選択装置

Country Status (1)

Country Link
JP (1) JP7211271B2 (ja)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009085785A (ja) * 2007-09-28 2009-04-23 Sanyo Electric Co Ltd クラック幅計測システム、操作装置、クラック幅計測方法、及びクラック幅計測プログラム
JP2015034445A (ja) * 2013-08-09 2015-02-19 ジビル調査設計株式会社 構造物の点検装置
JP2017227466A (ja) * 2016-06-20 2017-12-28 富士フイルム株式会社 計測支援装置及び計測支援方法
WO2019021729A1 (ja) * 2017-07-25 2019-01-31 富士フイルム株式会社 損傷図作成方法、損傷図作成装置、損傷図作成システム、及び記録媒体

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009085785A (ja) * 2007-09-28 2009-04-23 Sanyo Electric Co Ltd クラック幅計測システム、操作装置、クラック幅計測方法、及びクラック幅計測プログラム
JP2015034445A (ja) * 2013-08-09 2015-02-19 ジビル調査設計株式会社 構造物の点検装置
JP2017227466A (ja) * 2016-06-20 2017-12-28 富士フイルム株式会社 計測支援装置及び計測支援方法
WO2019021729A1 (ja) * 2017-07-25 2019-01-31 富士フイルム株式会社 損傷図作成方法、損傷図作成装置、損傷図作成システム、及び記録媒体

Also Published As

Publication number Publication date
JP7211271B2 (ja) 2023-01-24

Similar Documents

Publication Publication Date Title
JP6716996B2 (ja) 画像処理プログラム、画像処理装置、及び画像処理方法
JP7058277B2 (ja) 再構成方法および再構成装置
JP2007000205A (ja) 画像処理装置及び画像処理方法並びに画像処理プログラム
JP6319972B2 (ja) 画像処理装置、撮像装置、画像処理方法、およびプログラム
WO2019093457A1 (ja) 情報処理装置、情報処理方法及びプログラム
JP6594170B2 (ja) 画像処理装置、画像処理方法、画像投影システムおよびプログラム
JP6552256B2 (ja) 画像処理装置及び画像処理装置の制御方法
JP2021163190A (ja) 情報処理装置、情報処理方法、およびプログラム
JP2015524946A (ja) 画像解像度が改善された超解像画像を形成するための方法及び測定装置
US11941749B2 (en) Image display method, display control device, and recording medium for displaying shape image of subject and coordinates estimated from two-dimensional coordinates in reference image projected thereon
JPWO2019003796A1 (ja) 画像合成方法、画像合成装置、及び記録媒体
JP2004191856A5 (ja)
US20220405878A1 (en) Image processing apparatus, image processing method, and image processing program
JP2021128592A (ja) 画像処理装置、画像処理方法、学習済みモデルの生成方法、およびプログラム
JP7118729B2 (ja) 情報処理装置、情報処理方法及びプログラム
JP6494402B2 (ja) 画像処理装置、撮像装置、画像処理方法、プログラム
JP6371515B2 (ja) X線画像処理装置、x線画像処理方法、及びプログラム
JP6409769B2 (ja) 3次元形状計測装置、3次元形状計測方法及び3次元形状計測プログラム
WO2017187935A1 (en) Information processing apparatus, information processing method, and program
KR102076635B1 (ko) 산재된 고정 카메라를 이용한 파노라마 영상 생성 장치 및 방법
JP7211271B2 (ja) 画像選択方法、画像選択プログラム、および画像選択装置
JP2008217330A (ja) 速度推定方法および速度推定プログラム
JP2005141655A (ja) 3次元モデリング装置及び3次元モデリング方法
US20210304417A1 (en) Observation device and observation method
JP6355206B2 (ja) カメラキャリブレーション方法、カメラキャリブレーション装置及びカメラキャリブレーションプログラム

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20211208

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20221202

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20221226

R150 Certificate of patent or registration of utility model

Ref document number: 7211271

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150