WO2020122025A1

WO2020122025A1 - 地表面推定方法、計測領域表示システムおよびクレーン

Info

Publication number: WO2020122025A1
Application number: PCT/JP2019/048133
Authority: WO
Inventors: 巖石川; 孝幸小阪; 諭窪田; 田中　成典; 中村　健二; 雄平山本; 匡哉中原
Original assignee: 株式会社タダノ; 学校法人関西大学
Priority date: 2018-12-10
Filing date: 2019-12-09
Publication date: 2020-06-18
Also published as: US12032657B2; JP2020094835A; JP7129894B2; US20220017336A1; EP3896391A4; EP3896391A1

Abstract

他の領域によって分断されている領域について地表面の推定を行う地表面推定方法、計測領域表示システムおよびクレーンを提供する。点群データ（Ｐ）を演算処理するデータ処理部（７０）によって、地表面（Ｆ）と推定された連続領域に隣接していない連続領域である非隣接連続領域（Ｊａ）毎に、非隣接連続領域（Ｊａ）におけるグリッド（ｇ）の代表点（ｐｒ）の平均標高値を算出し、地表面（Ｆ）と推定された連続領域の所定の範囲（Ｎａ）におけるグリッド（ｇ）の代表点（ｐｒ）の平均標高値を算出し、非隣接連続領域（Ｊａ）におけるグリッド（ｇ）の代表点（ｐｒ）の平均標高値と地表面（Ｆ）と推定された連続領域の所定の範囲（Ｎａ）におけるグリッド（ｇ）の代表点（ｐｒ）の平均標高値との差が領域間の閾値以下である非隣接連続領域（Ｊａ）を地表面（Ｆ）と推定する第二地表面推定工程を備える。

Description

地表面推定方法、計測領域表示システムおよびクレーン

　本発明は、地表面推定方法、計測領域表示システムおよびクレーンに関する。

　従来、モニタに表示されたガイド情報に基づいて、荷物や荷物の周辺の状態を確認し、クレーン作業の効率化を図る技術が知られている。斯かる技術は、例えば、以下の特許文献１に開示されている。

　特許文献１には、吊荷（荷物）周辺の物体高さを通知する高さ情報通知システムに係る技術が開示されている。特許文献１記載の吊荷周辺の物体の高さ情報通知システムでは、レーザ距離センサ、マイクロ波距離センサ、ステレオカメラ等の距離計測手段によって、ブームの先端から吊荷周辺までの距離を計測する。そして、距離の計測結果を用いて吊荷周辺の物体の位置を検出するとともに高さを算出し、カメラによって撮像（撮影）された撮像画像（カメラ画像）に吊荷周辺の物体の位置および高さを対応させた処理画像（ガイド情報）を作成（生成）し、その処理画像をモニタに表示する構成としている。

　特許文献１に記載された従来技術は、吊荷周辺の領域において地表面の推定を行う構成ではない。そのため、従来技術は、他の領域によって分断されている領域がある場合に、分断されている領域の高さの情報を提供するだけである。つまり、分断されている領域について提供される情報には、地表面か否かをオペレータが判断できる情報は含まれていない。このような情報がクレーン作業に用いられた場合、地表面が存在していても、表示された高さの地物が存在するとオペレータが誤って認識してしまう可能性がある。

特開２０１３－１２０１７６号公報

　本発明の目的は、他の領域によって分断されている領域について地表面の推定を行う地表面推定方法、計測領域表示システムおよびクレーンの提供を目的とする。

　本発明の解決しようとする課題は以上の如くであり、次にこの課題を解決するための手段を説明する。

　即ち、本発明においては、レーザスキャナによって、計測領域にレーザを照射して点群データを取得し、前記点群データから照射点の標高値を算出する点群データ取得工程と、前記点群データを演算処理するデータ処理手段によって、前記計測領域を複数のグリッドに分割し、前記グリッド毎に、当該グリッドの重心位置と当該グリッド内の前記点群データの平均標高値とを算出し、前記グリッド内の点群データの平均標高値における前記重心位置を前記グリッド毎の代表点の位置として設定するグリッド生成処理工程と、前記複数のグリッドのうち、一のグリッドの代表点と隣接する他のグリッドの代表点との標高値の差がグリッド間の閾値以下であれば、前記一のグリッドと前記他のグリッドとが連続する領域である連続領域として認識する連続領域認識工程と、前記連続領域のうち、最も前記グリッドの数が多い連続領域を地表面と推定する第一地表面推定工程と、前記地表面と推定された連続領域に隣接していない連続領域である非隣接連続領域毎に、前記非隣接連続領域におけるグリッドの代表点の平均標高値を算出し、前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値を算出し、前記非隣接連続領域におけるグリッドの代表点の平均標高値と前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値との差が領域間の閾値以下である前記非隣接連続領域を地表面と推定する第二地表面推定工程と、を備える、ことを特徴とする地表面推定方法である。

　本発明においては、前記地表面と推定された連続領域の所定の範囲は、前記非隣接連続領域と、前記地表面と推定された連続領域と、の間にある連続領域に隣接するグリッドから構成される範囲である、ことを特徴とする地表面推定方法である。

　本発明においては、前記地表面と推定された非隣接連続領域の所定の範囲における、前記連続領域の数に応じて前記計測領域を複数のグリッドに分割する際の前記グリッドの幅を変更する、ことを特徴とする地表面推定方法である。

　本発明においては、前記第二地表面推定工程は、前記非隣接連続領域と、前記地表面と推定された連続領域と、の間の距離が所定の距離以上であれば、前記非隣接連続領域を地表面と推定しない、ことを特徴とする地表面推定方法である。

　本発明においては、計測領域にレーザを照射して点群データを取得するレーザスキャナを備えたデータ取得部と、取得した点群データを演算処理するデータ処理部と、データ表示部と、を備え、前記データ処理部は、前記データ取得部から前記点群データを取得し、前記点群データから照射点の標高値を算出し、前記計測領域を複数のグリッドに分割し、前記グリッド毎に、当該グリッドの重心位置と当該グリッド内の前記点群データの平均標高値とを算出し、前記グリッド内の点群データの平均標高値における前記重心位置を前記グリッド毎の代表点の位置として設定し、前記複数のグリッドのうち、一のグリッドの代表点と隣接する他のグリッドの代表点との標高値の差がグリッド間の閾値以下であれば、前記一のグリッドと前記他のグリッドとが連続する領域である連続領域として認識し、前記連続領域のうち、最も前記グリッドの数が多い連続領域を地表面と推定し、前記地表面と推定された連続領域に隣接していない連続領域である非隣接連続領域毎に、前記非隣接連続領域におけるグリッドの代表点の平均標高値を算出し、前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値を算出し、前記非隣接連続領域におけるグリッドの代表点の平均標高値と前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値との差が領域間の閾値以下である前記非隣接連続領域を地表面と推定し、前記地表面と推定された連続領域を他の領域と区別して前記データ表示部に表示する、ことを特徴とする計測領域表示システムである。

　本発明においては、旋回台と、前記旋回台に設けられるブームと、前記ブームに取り付けられ、点群データを取得するレーザスキャナと、取得した点群データを演算処理する制御装置と、表示装置と、を備えるクレーンにおいて、前記旋回台の旋回操作と、前記ブームの伸縮操作と起伏操作と、に伴って前記レーザスキャナを移動させながらレーザを照射させることで、前記レーザスキャナのレーザ照射時の位置毎の点群データを取得し、前記制御装置は、前記レーザスキャナのレーザ照射時の位置毎の前記点群データを取得し、前記レーザスキャナのレーザ照射時の位置とその姿勢とに基づいてレーザ照射時の位置毎の前記点群データを重ね合わせ、照射点の標高値を算出し、前記レーザスキャナの計測領域を複数のグリッドに分割し、前記グリッド毎に、当該グリッドの重心位置と当該グリッド内の前記点群データの平均標高値とを算出し、前記グリッド内の点群データの平均標高値における前記重心位置を前記グリッド毎の代表点の位置として設定し、前記複数のグリッドのうち、一のグリッドの代表点と隣接する他のグリッドの代表点との標高値の差がグリッド間の閾値以下であれば、前記一のグリッドと前記他のグリッドとが連続する領域である連続領域として認識し、前記連続領域のうち、最も前記グリッドの数が多い連続領域を地表面と推定し、前記地表面と推定された連続領域に隣接していない連続領域である非隣接連続領域毎に、前記非隣接連続領域におけるグリッドの代表点の平均標高値を算出し、前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値を算出し、前記非隣接連続領域におけるグリッドの代表点の平均標高値と前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値との差が領域間の閾値以下である前記非隣接連続領域を地表面と推定し、前記地表面と推定された連続領域を他の領域と区別して前記表示装置に表示する、ことを特徴とするクレーンである。

　本発明は、以下に示すような効果を奏する。

　本発明によれば、非隣接連続領域の標高値が地表面の標高値と異なっていても閾値の範囲内の差であれば、非隣接連続領域が地表面と推定される。これにより、他の領域によって分断されている領域について地表面の推定を行える。

　本発明によれば、地表面と推定されている領域の中で非隣接連続領域に最も近接している部分の標高値で判定するので、標高値の差が大きい地表面の影響を受けない。これにより、他の領域によって分断されている領域について地表面の推定を行える。

　本発明によれば、複数の非隣接連続領域が密集する場合、計測領域を複数のグリッドに分割する際のグリッドの幅を小さくすることによって、密集した非隣接連続領域がより詳細に認識される。これにより、他の領域によって分断されている領域について地表面の推定を行える。

　本発明によれば、地表面と推定された連続領域からの距離が長くなるに応じて分断している領域にある地表面の標高値が閾値の範囲内となることの不確実性が高くなるため、非隣接連続領域と、地表面と推定された連続領域と、の間の距離が所定の距離以上の非隣接連続領域を地表面と推定しないようにすることで、誤った推定が抑制される。これにより、他の領域によって分断されている領域について地表面の推定を行える。

　本発明によれば、旋回台やブームの動作に伴ってクレーンの作業領域における点群データが取得され、非隣接連続領域の標高値が地表面の標高値と異なっていても範囲内の差であれば、非隣接連続領域が地表面と推定される。これにより、クレーンの作業領域において、他の領域によって分断されている領域について地表面の推定を行える。

本発明の一実施形態に係るクレーンの全体構成を示す模式図。本発明の一実施形態に係る計測領域表示システムとクレーンとの制御構成を示す模式図。レーザスキャナによるレーザの照射状況の説明図、図３ＡはＹ軸方向視模式図、図３ＢはＺ軸方向視模式図。計測領域、個別計測領域、作業領域を説明する平面模式図。個別計測領域とレーザにより描かれる軌跡との関係を示す模式図。グリッド生成処理の説明図、図６Ａは点群データの分割、グリッドのスライス、平均標高値の算出の説明図、図６Ｂは代表点の設定の説明図。三次元地図更新処理の説明図、図７Ａは設定したグリッド毎の代表点と現在時刻よりも前に生成した三次元地図の代表点との説明図、図７Ｂは同一のグリッド内に複数の代表点が存在する場合における代表点の削除の説明図。グリッドの標高値を濃淡で表した三次元地図の表示状態を示す図。ラベリング処理の説明図、図９Ａは三次元地図のグリッドの模式図、図９Ｂはラベルの付与の説明図、図９Ｃはラベルの付与完了時の説明図、図９Ｄはラベルの上書き前の説明図、図９Ｅはラベルの上書き後の説明図。地物領域推定処理の説明図、図１０Ａは推定を行う連続領域と所定の範囲との説明図、図１０Ｂは標高値の平均値の説明図、図１０Ｃは地物と地表面との推定の説明図、図１０Ｄは所定の範囲の説明図。ガイド情報の表示状態を示す図、図１１Ａはカメラ画像を表示したデータ表示部を示す図、図１１Ｂはカメラ画像とガイド情報を重畳表示したデータ表示部を示す図。連続領域毎に領域を区別して表した三次元地図の表示状態を示す図。データ処理部によるデータ処理の流れを示すフロー図。グリッド生成処理を示すフロー図。三次元地図更新処理を示すフロー図。ラベリング処理を示すフロー図。地物領域推定処理を示すフロー図。グリッド生成処理においてグリッドの幅を変更する場合の説明図、図１８Ａは所定の範囲とその範囲内の連続領域との説明図、図１８Ｂはグリッドの幅を変更した後における地表面の領域の説明図。グリッド生成処理を示すフロー図。地物領域推定処理において非隣接連続領域を地表面と推定しない場合の説明図。地物領域推定処理を示すフロー図。

　以下に、本発明の第一実施形態に係る計測領域表示システム５０を備えたクレーン１について説明する。尚、本実施形態においては、クレーン１として移動式クレーン（ラフテレーンクレーン）について説明を行うが、トラッククレーン等でもよい。

　図１に示すように、クレーン１は、不特定の場所に移動可能な移動式クレーンである。クレーン１は、車両２、クレーン装置６を有する。

　車両２は、クレーン装置６を搬送する走行車両である。車両２は、複数の車輪３を有し、エンジン４を動力源として走行する。車両２には、アウトリガ５が設けられている。アウトリガ５は、車両２の幅方向両側に油圧によって延伸可能な張り出しビームと地面に垂直な方向に延伸可能な油圧式のジャッキシリンダとから構成されている。

　クレーン装置６は、荷物Ｗをワイヤロープによって吊り上げる作業装置である。クレーン装置６は、旋回台７、ブーム９、ジブ９ａ、メインフックブロック１０、サブフックブロック１１、起伏用油圧シリンダ１２、メインウインチ１３、メインワイヤロープ１４、サブウインチ１５、サブワイヤロープ１６、キャビン１７、制御装置１８（図２参照）等を具備する。

　旋回台７は、クレーン装置６を旋回可能に構成する駆動装置である。旋回台７は、円環状の軸受を介して車両２のフレーム上に設けられる。旋回台７は、円環状の軸受の中心を回転中心として回転自在に構成されている。旋回台７には、アクチュエータである油圧式の旋回用油圧モータ８が設けられている。

　旋回用油圧モータ８は、電磁比例切換弁である旋回用バルブ２３（図２参照）によって回転操作されるアクチュエータである。旋回用バルブ２３は、旋回用油圧モータ８に供給される作動油の流量を任意の流量に制御することができる。つまり、旋回台７は、旋回用バルブ２３によって回転操作される旋回用油圧モータ８を介して任意の旋回速度に制御可能に構成されている。旋回台７には、旋回台７の旋回位置（角度）と旋回速度とを検出する旋回角度検出手段である旋回用センサ２７（図２参照）が設けられている。

　ブーム９は、荷物Ｗを吊り上げ可能な状態にワイヤロープを支持する可動支柱である。ブーム９は、複数のブーム部材から構成されている。ブーム９は、各ブーム部材をアクチュエータである図示しない伸縮用油圧シリンダで移動させることで軸方向に伸縮自在に構成されている。ブーム９は、ベースブーム部材の基端が旋回台７の略中央に揺動可能に設けられている。

　図示しない伸縮用油圧シリンダは、電磁比例切換弁である伸縮用バルブ２４（図２参照）によって伸縮操作されるアクチュエータである。伸縮用バルブ２４は、伸縮用油圧シリンダに供給される作動油の流量を任意の流量に制御することができる。つまり、ブーム９は、伸縮用バルブ２４によって任意のブーム長さに制御可能に構成されている。ブーム９には、ブーム９の長さを検出する伸縮長さ検出手段である伸縮用センサ２８（図２参照）が設けられている。

　メインフックブロック１０とサブフックブロック１１とは、荷物Ｗを吊る吊り具である。メインフックブロック１０には、メインワイヤロープ１４が巻き掛けられる複数のフックシーブと、荷物Ｗを吊るメインフック１０ａとが設けられている。サブフックブロック１１には、荷物Ｗを吊るサブフック１１ａが設けられている。

　起伏用油圧シリンダ１２は、電磁比例切換弁である起伏用バルブ２５（図２参照）によって伸縮操作される。起伏用バルブ２５は、起伏用油圧シリンダ１２に供給される作動油の流量を任意の流量に制御することができる。つまり、ブーム９は、起伏用バルブ２５によって任意の起伏速度に制御可能に構成されている。ブーム９には、ブーム９の起伏角度を検出する旋回角度検出手段である起伏用センサ２９（図２参照）や荷物Ｗの重量を検出する重量センサ等が設けられている。

　メインウインチ１３とサブウインチ１５とは、メインワイヤロープ１４とサブワイヤロープ１６との繰り入れ（巻き上げ）および繰り出し（巻き下げ）を行う巻回装置である。メインウインチ１３は、メインワイヤロープ１４が巻きつけられるメインドラムがアクチュエータである図示しないメイン用油圧モータによって回転され、サブウインチ１５は、サブワイヤロープ１６が巻きつけられるサブドラムがアクチュエータである図示しないサブ用油圧モータによって回転されるように構成されている。

　メイン用油圧モータは、電磁比例切換弁であるメイン用バルブ２６ｍ（図２参照）によって回転操作される。メイン用バルブ２６ｍは、メイン用油圧モータに供給される作動油の流量を任意の流量に制御することができる。つまり、メインウインチ１３は、メイン用バルブ２６ｍによって任意の繰り入れおよび繰り出し速度に制御可能に構成されている。同様に、サブウインチ１５は、電磁比例切換弁であるサブ用バルブ２６ｓ（図２参照）によって任意の繰り入れおよび繰り出し速度に制御可能に構成されている。メインウインチ１３とサブウインチ１５とには、メインワイヤロープ１４とサブワイヤロープ１６との繰り出し量をそれぞれ検出する巻回用センサ３０（図２参照）が設けられている。

　キャビン１７は、運転座席を覆うものである。運転座席には、車両２を走行操作するための操作具やクレーン装置６を操作するための旋回操作具１９、伸縮操作具２０、起伏操作具２１、メインドラム操作具２２ｍ、サブドラム操作具２２ｓ等が設けられている（図２参照）。旋回操作具１９は、旋回用油圧モータ８を操作することができる。伸縮操作具２０は、伸縮用油圧シリンダを操作することができる。起伏操作具２１は、起伏用油圧シリンダ１２を操作することができる。メインドラム操作具２２ｍは、メイン用油圧モータを操作することができる。サブドラム操作具２２ｓは、サブ用油圧モータを操作することができる。

　このように構成されるクレーン１は、車両２を走行させることで任意の位置にクレーン装置６を移動させることができる。また、クレーン１は、起伏操作具２１の操作によって起伏用油圧シリンダ１２でブーム９を任意の起伏角度に起立させて、伸縮操作具２０の操作によってブーム９を任意のブーム長さに延伸させたりすることでクレーン装置６の揚程や作業半径を拡大することができる。また、クレーン１は、メインドラム操作具２２ｍ等によって荷物Ｗを吊り上げて、旋回操作具１９の操作によって旋回台７を旋回させることで荷物Ｗを搬送することができる。

　図２に示すように、制御装置１８は、各バルブを介してクレーン装置６のアクチュエータを制御する。制御装置１８は、キャビン１７内に設けられている。制御装置１８は、実体的には、ＣＰＵ、ＲＯＭ、ＲＡＭ、ＨＤＤ等がバスで接続される構成であってもよく、あるいはワンチップのＬＳＩ等からなる構成であってもよい。制御装置１８は、各アクチュエータ、センサ等の動作を制御するために種々のプログラムやデータが格納されている。

　制御装置１８は、旋回操作具１９、伸縮操作具２０、起伏操作具２１、メインドラム操作具２２ｍ、サブドラム操作具２２ｓに接続され、旋回操作具１９、伸縮操作具２０、起伏操作具２１、メインドラム操作具２２ｍ、サブドラム操作具２２ｓのそれぞれの操作量を取得することができる。

　制御装置１８は、旋回用バルブ２３、伸縮用バルブ２４、起伏用バルブ２５、メイン用バルブ２６ｍ、サブ用バルブ２６ｓに接続され、旋回用バルブ２３、伸縮用バルブ２４、起伏用バルブ２５、メイン用バルブ２６ｍ、サブ用バルブ２６ｓに制御信号を伝達することができる。

　制御装置１８は、旋回用センサ２７、伸縮用センサ２８、起伏用センサ２９、重量センサ、巻回用センサ３０に接続され、旋回台７の旋回位置、ブーム長さ、起伏角度、荷物Ｗの重量、メインワイヤロープ１４とサブワイヤロープ１６との繰り出し量を取得することができる。

　制御装置１８は、旋回台７の旋回位置、ブーム長さ、起伏角度、荷物Ｗの重量、メインワイヤロープ１４とサブワイヤロープ１６との繰り出し量等のクレーン１に関する情報を計測領域表示システム５０に伝達することができる。

　尚、本説明では、ブーム９の起伏支点の軸方向を基準として、図１に示すようなＸＹＺ座標系を規定している（以下の説明においても同様）。
　Ｘ軸方向（奥行方向とも呼ぶ）は、ブーム９の起伏支点の軸方向に対して垂直、かつ、水平な方向である。また、Ｙ軸方向（水平方向とも呼ぶ）は、ブーム９の起伏支点の軸方向に対して平行、かつ、水平な方向である。さらに、Ｚ軸方向は、鉛直下方向である。即ち、ＸＹＺ座標系は、ブーム９を基準としたローカル座標系として規定している（図４参照）。

　次に、本発明の一実施形態に係る計測領域表示システム５０について、説明する。
　クレーン１は、図２に示すような計測領域表示システム５０を備えている。
　計測領域表示システム５０は、本発明に係る計測領域表示システム５０の一例であり、図１に示すようなクレーン１による作業を効率よく、かつ、安全に行うことを可能にするために、計測領域表示システム５０が計測対象とする領域（以下、計測領域ＫＡという（図４参照））の情報（以下、ガイド情報という）を画像で表示し、オペレータに提示するためのシステムである。

　図２に示すように、計測領域表示システム５０は、データ取得部６０、データ処理部７０、データ表示部８０、データ入力部９０によって、構成されている。

　データ取得部６０は、計測領域ＫＡにおけるガイド情報を生成するために必要なデータを取得する部位であり、カメラ６１、レーザスキャナ６２、慣性計測装置（ＩＭＵ）６３、第一ＧＮＳＳ受信機６４ａ、第二ＧＮＳＳ受信機６４ｂを備えている。カメラ６１、レーザスキャナ６２、慣性計測装置（ＩＭＵ）６３、第一ＧＮＳＳ受信機６４ａは、フレーム体に対して固定され、一体に構成されたセンサユニット６５となっている。

　センサユニット６５は、クレーン１のブーム９の先端部分に付設されており、荷物Ｗの真上に位置するブーム９の先端部分から真下の状況を捉えることができる状態で配置されている（図１参照）。尚、ここでいう荷物Ｗの「真上」は、荷物Ｗの鉛直上方の位置と、その位置を基準とした一定範囲（例えば、荷物Ｗの上面の範囲）の位置と、を含む概念である。

　センサユニット６５は、ブーム９の先端部分に対してジンバル６６（図１参照）を介して付設されており、旋回台７の旋回操作、ブーム９の起伏操作、伸縮操作が行われたときに、センサユニット６５の姿勢（Ｚ軸方向に向けた姿勢）を略一定に保持することができるように構成されている。これにより、カメラ６１とレーザスキャナ６２とを常に荷物Ｗに向けておくことができる。このため、センサユニット６５は、カメラ６１とレーザスキャナ６２とによって、荷物Ｗとその下方の地表面Ｆから、常にデータを取得することができる。また、荷物Ｗの下方の地表面Ｆに地物Ｅが存在する場合には、カメラ６１とレーザスキャナ６２とによって、地物Ｅのデータを取得することができる。

　カメラ６１は、センサユニット６５の下方の領域（以下、個別計測領域ｋａという（図４参照））の画像を撮影するためのデジタルビデオカメラであり、撮影したカメラ画像をリアルタイムで外部に出力する機能を有している。また、カメラ６１は、適切なガイド情報の生成に必要なデータ量を考慮した画素数、画角、フレームレート、画像伝送レートを有している。

　レーザスキャナ６２は、計測対象物にレーザを照射し、そのレーザの計測対象物における反射光を受光することによって、その反射点に係る情報を取得し、計測対象物の点群データを取得する装置である。レーザスキャナ６２は、センサユニット６５を介してブーム９の先端部分に取り付けられている。レーザスキャナ６２の計測対象物は、荷物Ｗ、地物Ｅ、地表面Ｆである。また、レーザスキャナ６２には、計測時刻を取得するための第三ＧＮＳＳ受信機６４ｃが接続されている。
　そして、レーザスキャナ６２は、個別計測領域ｋａに向けて照射されるレーザにより描かれる軌跡がＹ軸方向に対して平行となるように配置されている（図３Ｂ参照）。また、レーザスキャナ６２は、レーザの照射角度を変更する基準軸が、Ｘ軸方向に対して平行とされている。
　計測領域表示システム５０は、レーザスキャナ６２によって、リアルタイムに平面的な三次元点群データを取得する。

　慣性計測装置（Ｉｎｅｒｔｉａｌ　Ｍｅａｓｕｒｅｍｅｎｔ　Ｕｎｉｔ、以下ＩＭＵという）６３は、データ取得時におけるカメラ６１とレーザスキャナ６２との姿勢データを取得するための装置である。ＩＭＵ６３は、リアルタイムで姿勢角を計測することが可能であり、レーザスキャナ６２によって取得した点群データの補正に利用可能な計測精度を有している。また、ＩＭＵ６３には、計測時刻を取得するための第四ＧＮＳＳ受信機６４ｄが接続されている。

　第一ＧＮＳＳ受信機６４ａは、衛星から測距電波を受信し、座標である緯度、経度、標高値を算出するための装置である。データ処理部７０は、第一ＧＮＳＳ受信機６４ａとレーザスキャナ６２とが離間する距離が設定されており、設定された距離に基づいてリアルタイムでレーザスキャナ６２の座標を算出可能である。

　第二ＧＮＳＳ受信機６４ｂは、衛星から測距電波を受信し、座標である緯度、経度、標高値を算出するための装置である。第二ＧＮＳＳ受信機６４ｂは、旋回台７の旋回中心位置に配置されている。第二ＧＮＳＳ受信機６４ｂは、リアルタイムで旋回台７の旋回中心の座標を算出可能である。

　本実施形態では、第一ＧＮＳＳ受信機６４ａと第二ＧＮＳＳ受信機６４ｂとには、測定精度の高いＲＴＫ－ＧＰＳ（Ｒｅａｌ　Ｔｉｍｅ　Ｋｉｎｅｍａｔｉｃ　ＧＰＳ）測位方式を採用する。ＲＴＫ－ＧＰＳ測位方式を採用することにより、レーザスキャナ６２の位置と旋回台７の旋回中心の位置の測定精度を高めることができる。尚、ＲＴＫ－ＧＰＳ測位方式に限定されず、他の測位方法を採用してもよい。

　第一ＧＮＳＳ受信機６４ａと第二ＧＮＳＳ受信機６４ｂとを結ぶ直線と、レーザスキャナ６２の計測軸と、ＩＭＵ６３の計測軸と、が同一直線上にあるように、第一ＧＮＳＳ受信機６４ａが配置される。また、第一ＧＮＳＳ受信機６４ａが算出するレーザスキャナ６２の座標と、第二ＧＮＳＳ受信機６４ｂが算出する旋回台７の旋回中心の座標と、によって、ブーム９を基線としたＧＮＳＳコンパスを構成し、レーザスキャナ６２の向き（方位）を算出可能である。第一ＧＮＳＳ受信機６４ａと第二ＧＮＳＳ受信機６４ｂとは、レーザスキャナ６２によって取得した点群データの補正に利用可能な計測精度を有している。

　図３に示すように、レーザスキャナ６２は、合計１６個のレーザ送受信センサを備えており、同時に１６本のレーザを計測対象物に照射して、計測対象物の点群データを取得することができる装置である。図３Ａに示すように、レーザスキャナ６２の１６個の各レーザ送受信センサは、Ｘ軸方向において２°ずつ照射角度を異ならせて配置されており、計測対象物に対して、全体で３０°の拡がりを持ってレーザを照射可能に構成されている。また、レーザスキャナ６２の各レーザ送受信センサは、Ｘ軸回りに３６０°（全方位）回転可能に構成されている。図３Ｂに示すように、個別計測領域ｋａに向けて照射されるレーザにより描かれる軌跡は、Ｙ軸方向に対して平行であり、レーザスキャナ６２では、１６本の当該軌跡が同時に描かれる。

　尚、レーザスキャナ６２は、ブーム９の最高到達高さを考慮して、その最高到達高さ（例えば、約１００ｍ）から計測対象物の三次元形状を計測可能な機器が選択される。また、レーザスキャナ６２は、適切なガイド情報を生成するために必要なデータ量およびデータ精度を考慮して、計測スピード、計測ポイント数、計測精度等の各仕様について所定の性能を有する機器が選択される。

　尚、本実施形態では、合計１６個のレーザ送受信センサを備えたレーザスキャナ６２を用いる場合を例示しているが、本発明に係る計測領域表示システム５０は、レーザスキャナ６２を構成するレーザ送受信センサの個数によっては限定されない。即ち、本発明に係る計測領域表示システム５０では、クレーン１のブーム９（ジブ９ａ）の最高到達高さ等に応じて、最適な仕様のレーザスキャナ６２が適宜選択される。

　センサユニット６５によって個別計測領域ｋａにおいて取得するデータには、荷物Ｗ、荷物Ｗの下方の地表面Ｆ、荷物Ｗの下方の地表面Ｆに存在する地物Ｅをカメラ６１によって撮影した画像データが含まれる。また、センサユニット６５によって個別計測領域ｋａにおいて取得するデータには、荷物Ｗ、地表面Ｆ、地物Ｅをレーザスキャナ６２によってスキャンして取得した点群データが含まれる。尚、ここでいう地表面Ｆには、荷物Ｗの搬送元および搬送先となる面を広く含み、地上面のみならず、建物屋上の床面や屋根面等も含まれる。

　データ処理手段であるデータ処理部７０は、データ取得部６０で取得したデータを演算処理して、オペレータに提示するガイド情報を生成するための部位であり、本実施形態では、所定のデータ処理プログラムがインストールされた汎用のパーソナルコンピュータによって構成している。
　また、データ処理部７０は、クレーン１の制御装置１８と電気的に接続されており、制御装置１８から出力される旋回台７の旋回位置、ブーム長さ、起伏角度、荷物Ｗの重量、メインワイヤロープ１４とサブワイヤロープ１６との繰り出し量等のクレーン１に関する情報が入力される。

　データ表示部８０は、オペレータに提示するガイド情報を表示するための部位であり、データ処理部７０に接続された表示装置により構成される。
　データ表示部８０は、計測領域ＫＡの三次元地図と個別計測領域ｋａのカメラ画像とをリアルタイムに表示する。

　ここで、計測領域ＫＡとガイド情報とについて説明する。

　計測領域ＫＡとは、ガイド情報を生成する対象となる領域である。
　図４に示すように、計測領域表示システム５０は、センサユニット６５の下方の個別計測領域ｋａに存在する荷物Ｗ、地表面Ｆ、地物Ｅを計測する。センサユニット６５は、旋回台７の旋回操作、起伏操作、ブーム９の伸縮操作に応じて移動する。このため、個別計測領域ｋａも、センサユニット６５の移動に応じて移動する。計測領域ＫＡは、異なる位置から計測された個別計測領域ｋａを併せた領域である。計測領域ＫＡは、計測領域表示システム５０が計測可能な領域の全体を計測した場合、ブーム９の先端部分が稼動可能な領域である作業領域ＳＡを含む領域となる。

　ここでいうガイド情報は、オペレータがクレーン１によって荷物Ｗを搬送するときに、ブーム９の長さ・旋回位置・起伏角度、メインワイヤロープ１４とサブワイヤロープ１６との繰り出し量等の良否について、オペレータの判断を補助する情報である。ガイド情報には、計測領域ＫＡの三次元地図、カメラ画像情報、荷物Ｗと地物Ｅとの形状に係る情報、荷物Ｗの高さ情報、地物Ｅの高さ情報、荷物Ｗの動線に係る情報等が含まれる。

　図２に示すように、データ入力部９０は、データ処理部７０に対して、設定値等を入力するための部位であり、タッチパネル、マウス、キーボード装置等の入力装置により構成される。

　図１に示すように、データ表示部８０とデータ入力部９０は、キャビン１７内の運転座席の前方のオペレータが見やすい位置に配置する。データ処理部７０は、センサユニット６５の近傍に配置することが好ましい。尚、計測領域表示システム５０は、データ処理部７０、データ表示部８０、データ入力部９０をタブレットＰＣによって一体的に構成した場合には、データ処理部７０をキャビン１７内に配置する構成としてもよい。
　データ取得部６０とデータ処理部７０間のデータの伝送は、有線ＬＡＮによることが好ましい。尚、データ取得部６０とデータ処理部７０間のデータの伝送は、無線ＬＡＮを採用してもよく、あるいは、電力線通信を採用してもよい。

　ここで、データ取得部６０によるデータの取得状況を説明する。
　データ取得部６０は、カメラ６１によって、個別計測領域ｋａを連続的に撮影し、個別計測領域ｋａのカメラ画像を取得する。

　図５に示すように、データ取得部６０は、レーザスキャナ６２によって、個別計測領域ｋａを連続的にスキャンし、個別計測領域ｋａにおける計測対象物の点群データを取得する。以下では、レーザスキャナ６２によって取得する点群データを、点群データＰと呼ぶ。点群データＰは、点データｐの集合であり、点データｐは、個別計測領域ｋａに存在する地表面Ｆ、荷物Ｗ、地物Ｅに位置する点を表している。

　データ取得部６０は、レーザスキャナ６２によって点群データＰを取得すると同時に、第三ＧＮＳＳ受信機６４ｃによって複数の測位衛星から時間情報を受信する。そして、データ処理部７０は、点データｐに対して、該点データｐの取得時間ｔｐに係る情報を付与する。

　また、データ取得部６０は、レーザスキャナ６２によって、点群データＰを取得すると同時に、ＩＭＵ６３によって所定の周期でレーザスキャナ６２の姿勢データＱを取得し、第一ＧＮＳＳ受信機６４ａと第二ＧＮＳＳ受信機６４ｂとによってレーザスキャナ６２の位置データＲと方位データＫとを取得する。尚、データ処理部７０は、第二ＧＮＳＳ受信機６４ｂによって取得される旋回台７の旋回中心の位置と、制御装置１８から出力される旋回台７の旋回位置、ブーム長さ、起伏角度に基づいて、レーザスキャナ６２の位置と向きとを算出してもよい。

　姿勢データＱには、レーザスキャナ６２のＸ・Ｙ・Ｚ軸の各軸方向に対する角度と加速度に係る情報が含まれる。尚、ＩＭＵ６３による姿勢データＱの取得周期は、レーザスキャナ６２による点群データＰの取得周期よりも短くする。姿勢データＱは、計測周期ごとに計測される個別姿勢データｑの集合である。

　データ取得部６０は、ＩＭＵ６３によって姿勢データＱを取得すると同時に、第四ＧＮＳＳ受信機６４ｄによって、複数の測位衛星から時間情報を受信する。データ処理部７０は、個別姿勢データｑに対して、該個別姿勢データｑの取得時間に係る情報として取得時間ｔｑを付与する。

　位置データＲには、レーザスキャナ６２のＸＹＺ座標系における位置に係る情報が含まれる。尚、第一ＧＮＳＳ受信機６４ａと第二ＧＮＳＳ受信機６４ｂとによる位置データＲの取得周期は、レーザスキャナ６２による点群データＰの取得周期よりも短くする。位置データＲは、計測周期ごとに計測される個別位置データｒの集合である。データ処理部７０は、個別位置データｒに基づいて個別方位データｋを算出する。方位データＫは、計測周期ごとに計測される個別方位データｋの集合である。

　第一ＧＮＳＳ受信機６４ａと第二ＧＮＳＳ受信機６４ｂとは、位置データＲと方位データＫとを取得すると同時に、複数の測位衛星から時間情報を受信する。データ処理部７０は、位置データＲと方位データＫとに対して、該位置データＲと該方位データＫとの取得時間に係る情報として取得時間ｔｒを付与する。

　次に、データ処理部７０によるデータの処理状況を説明する。

　本発明の一実施形態に係る地表面推定方法は、点群データ取得工程、グリッド生成処理工程、連続領域認識工程、第一地表面推定工程、第二地表面推定工程を備える。

（点群データ取得工程）
　データ処理部７０は、計測データを取得する。具体的には、データ処理部７０は、点群データＰのストリームデータから、１フレーム分の点群データＰを切り出して出力する。１フレーム分の点群データＰは、レーザスキャナ６２によるレーザの照射方向がＸ軸回りに１周する間に取得する点データｐの集合である。また、データ処理部７０は、カメラ６１のカメラ画像、ＩＭＵ６３の姿勢データＱ、第一ＧＮＳＳ受信機６４ａと第二ＧＮＳＳ受信機６４ｂとによって取得する位置データＲおよび方位データＫを計測データとして取得する。

　データ処理部７０は、１フレーム分の点群データＰに含まれる点データｐを、ＩＭＵ６３によって取得した姿勢データＱと、第一ＧＮＳＳ受信機６４ａと第二ＧＮＳＳ受信機６４ｂとによって取得した位置データＲおよび方位データＫと、時間によって対応付ける（図１３の時刻対応処理（ＳＴＥＰ－１０１）参照）。
　具体的には、データ処理部７０は、姿勢データＱについて、個々の点データｐにおいて、その点データｐの取得時間ｔｐに最も近い個別姿勢データｑの取得時間ｔｑを探索し、該取得時間ｔｑにおける個別姿勢データｑを、その点データｐに対応付ける。また、データ処理部７０は、位置データＲおよび方位データＫについて、個々の点データｐにおいて、その点データｐの取得時間ｔｐに最も近い位置データＲおよび方位データＫの取得時間ｔｒを探索し、該取得時間ｔｒにおける個別位置データｒおよび個別方位データｋを、その点データｐに対応付ける。

　このようにしてデータ処理部７０は、個別姿勢データｑ、個別位置データｒ、個別方位データｋとが時間によって対応付けられた点データｐを出力する。

　データ処理部７０は、時間によって対応付けられた点データｐ、個別姿勢データｑ、個別位置データｒ、個別方位データｋの組み合わせに対して剛体変換処理を行う（図１３の剛体変換処理（ＳＴＥＰ－１０２）参照）。具体的には、データ処理部７０は、個別姿勢データｑ、個別位置データｒ、個別方位データｋとに基づいて点群データＰのアフィン変換を行って、点群データＰを平面直角座標系に変換する。これにより、傾き、位置、向きが補正された点群データＰが出力される。点データｐの標高値は、傾き、位置、向きが補正された点データｐのＺ座標の値となる。つまり、レーザスキャナ６２の照射点の標高値は、点データｐのＺ座標の値として算出されている。但し、Ｚ軸は、鉛直下方向の軸であるため、Ｚ座標の値が大きくなるに応じて、標高が低い位置となる。

（グリッド生成処理工程）
　図６に示すように、データ処理部７０は、補正された点群データＰの全体を一定の大きさのグリッドｇで分割し、グリッドｇの代表点ｐｒを設定する（図１３のグリッド生成処理Ａ（ＳＴＥＰ－１０３））。具体的には、まず、データ処理部７０は、補正された点群データＰを重ね合わせる。図６Ａに示すように、次に、データ処理部７０は、重ね合わせた点群データＰの全体を一定の大きさのグリッドｇで分割する。

　重ね合わせた点群データＰには、地物Ｅや地表面Ｆの上面だけではなく、地物Ｅや地表面Ｆの側面にレーザが照射され取得された点データｐが含まれる可能性がある。地物Ｅや地表面Ｆの上面を対象としてガイド情報を生成するため、まず、データ処理部７０は、グリッドｇ毎に、グリッドｇを任意の間隔で水平の方向にスライスする。そして、データ処理部７０は、グリッドｇ毎に、スライスされて形成される層Ｌ１・Ｌ２・Ｌ３・・・のうち、点データｐが存在し、かつ、最も標高が高い位置にある層Ｌ１の点データｐ（点群データＰ）の標高値の平均値である平均標高値Ｈを算出する。また、データ処理部７０は、グリッドｇ毎に、グリッドｇの重心位置を算出する。図６Ｂに示すように、最後に、データ処理部７０は、グリッドｇ毎に、算出した平均標高値Ｈにおけるグリッドｇの重心位置をグリッドｇの代表点ｐｒの位置として設定する。これにより、データ処理部７０は、地表面Ｆに凹凸があり、地表面Ｆの凹凸の側面にレーザが照射された場合でも、地表面Ｆの凹凸の上面にレーザが照射され取得されたと推定される点群データＰの代表点ｐｒを設定できる。

　尚、平均標高値Ｈを算出する層は、層Ｌ１に限定されず、各層における、点データｐの数、点データｐの位置、点データｐの密度に基づいて選択される層であってもよい。例えば、各層のうち、点データｐの数が最も多い層が選択されてもよい。或いは、任意の数以上の点データｐが存在する層のうち、最も高い層が選択されてもよい。或いは、Ｚ方向視で均等に点データｐが存在する層は、地物Ｅの上面に照射された点データｐである可能性が高いため、Ｚ方向視で均等に点データｐが存在する層が選択されてもよい。

　図７に示すように、データ処理部７０は、設定したグリッドｇ毎の代表点ｐｒを用いて現在時刻よりも前に生成した三次元地図全体を更新する（図１３の三次元地図更新処理Ｂ（ＳＴＥＰ－１０４））。図７Ａ左図は、設定したグリッドｇ毎の代表点ｐｒを代表点ｐｒ１として、グリッドｇに代表点ｐｒ１のみを示した図である。図７Ａ右図は、現在時刻よりも前に生成した三次元地図の代表点ｐｒを代表点ｐｒ２として、グリッドｇに代表点ｐｒ２のみを示した図である。具体的には、まず、データ処理部７０は、代表点ｐｒ１と代表点ｐｒ２とをグリッドｇ毎に取得する。図７Ｂに示すように、データ処理部７０は、同一のグリッドｇ内に複数の代表点ｐｒが存在する場合は、代表点ｐｒ２を削除する。データ処理部７０は、削除後の全ての代表点ｐｒ毎にグリッドｇの大きさの面を生成し、計測領域ＫＡの三次元地図として出力する。図８に示すように、三次元地図Ｍをデータ表示部８０に表示した場合、グリッドｇの代表点ｐｒの標高値が高くなるに応じて面Ｓの色が濃く表示される。

（連続領域認識工程）
　図９に示すように、データ処理部７０は、三次元地図を用いて連続領域毎にラベルを付与する（図１３のラベリング処理Ｃ（ＳＴＥＰ－１０５）参照）。連続領域とは、隣接するグリッドｇの代表点ｐｒの標高値の差がグリッドｇ間の閾値であるグリッド閾値以下となる領域である。グリッド閾値は、レーザスキャナ６２のＺ軸方向の計測誤差に基づいて設定されており、例えば、０．１ｍとしている。ラベルは、１番から順番に番号で付与される。以下では、説明のため、グリッドｇの代表点ｐｒの標高値は、標高が最も低い位置に存在する代表点ｐｒの標高値を基準として、基準とした代表点ｐｒからの高さとしている。

　図９Ａ及び図９Ｂに示すように、具体的には、まず、データ処理部７０は、三次元地図のグリッドｇから注目グリッドｇａを選択する。注目グリッドｇａの代表点ｐｒの標高値と、その左上、上、右上、左の４近傍（以下、単に４近傍という）にある各グリッドｇの代表点ｐｒの標高値と、の比較が行われる。注目グリッドｇａが選択される順番は、三次元地図の左上を起点に左端から右に、右端に到達すると一つ下の段の左端から右に選択される。

　データ処理部７０は、注目グリッドｇａの代表点ｐｒの標高値とその４近傍にある各グリッドｇの代表点ｐｒの標高値との差のいずれかがグリッド閾値以下の場合は、標高値の差がグリッド閾値以下であるグリッドｇと同一のラベルを注目グリッドｇａに付与し、注目グリッドｇａの代表点ｐｒの標高値とその４近傍にある各グリッドｇの代表点ｐｒの標高値との差が全てグリッド閾値よりも大きい場合は、新たなラベルを付与する。新たなラベルとは、この時点までに付与されたラベルの番号のうち、最も大きい番号を１だけ増加させた番号のラベルである。但し、データ処理部７０は、注目グリッドｇａの代表点ｐｒの標高値とその４近傍にある各グリッドｇの代表点ｐｒの標高値との差がグリッド閾値以下であるグリッドｇのラベルの種類が複数の場合は、標高値の差がグリッド閾値以下であるグリッドｇのラベルのうち、最も番号が小さいラベルを付与する。これらの処理を全てのグリッドｇに対して行うことで、注目グリッドｇａとその左上、上、右上、左、右、左下、下、右下の８近傍にあるグリッドｇとの比較ができる。

　例えば、データ処理部７０は、図９Ｂに示す注目グリッドｇａが選択されている場合、注目グリッドｇａの代表点ｐｒの標高値とその４近傍にある各グリッドｇの代表点ｐｒの標高値との差を、全て０．５ｍと算出する。データ処理部７０は、注目グリッドｇａの代表点ｐｒの標高値とその４近傍にある各グリッドｇの代表点ｐｒの標高値との差が全てグリッド閾値０．１ｍよりも大きいため、１番のラベルの番号を１だけ増加させた２番のラベルを注目グリッドｇａに付与している。

　図９Ｃ、図９Ｄ、及び図９Ｅに示すように、データ処理部７０は、全てのグリッドｇにラベルが付与された後に、グリッドｇから注目グリッドｇａを再び選択する。注目グリッドｇａが選択される順番は、ラベルを付与する際と同様に、三次元地図の左上を起点に左端から右に、右端に到達すると一つ下の段の左端から右に選択される。

　データ処理部７０は、注目グリッドｇａの代表点ｐｒの標高値とその４近傍にある各グリッドｇの代表点ｐｒの標高値とを比較する。データ処理部７０は、注目グリッドｇａの４近傍のグリッドｇに、標高値の差がグリッド閾値以下で、かつ、注目グリッドｇａと異なるラベルが付与されている場合は、標高値の差がグリッド閾値以下で、かつ、注目グリッドｇａと異なるラベルが付与されているグリッドｇのラベルを注目グリッドｇａのラベルで上書きする。

　例えば、データ処理部７０は、図９Ｄに示す注目グリッドｇａが選択されている場合、その左と右上にある各グリッドｇの代表点ｐｒとの標高値の差を、それぞれ０ｍと算出し、その左上と上にある各グリッドｇの代表点ｐｒとの標高値の差を、それぞれ２ｍと算出する。さらに、注目グリッドｇａの左にあるグリッドｇは、５番のラベルが付与されており、注目グリッドｇａに付与されている３番のラベルと異なる。従って、注目グリッドｇａの左にあるグリッドｇは、標高値の差がグリッド閾値０．１ｍ以下で、かつ、注目グリッドｇａと異なるラベルが付与されている。図９Ｅに示すように、このため、データ処理部７０は、注目グリッドｇａの左にあるグリッドｇのラベルを３番のラベルで上書きしている。

　最後に、データ処理部７０は、グリッドｇにラベルが付与された三次元地図を出力する。このように、データ処理部７０は、隣接するグリッドｇの間で代表点ｐｒの標高値を比較して、傾斜の有無によらず、同一のラベルが付与されたグリッドｇを、連続領域として認識している。

（第一地表面推定工程）
　データ処理部７０は、ラベルが付与された三次元地図の中から、地表面Ｆ、荷物Ｗとそれら以外の領域を示すラベルを推定する（図１３の同一領域推定処理（ＳＴＥＰ－１０６））。具体的には、データ処理部７０は、地表面Ｆの領域が最も広いと考えられるため、グリッドｇの数が最も多い同一のラベルを地表面Ｆのグリッドｇに付与されたラベルとする。つまり、データ処理部７０は、連続領域のうち、最もグリッドｇの数が多い連続領域を地表面Ｆと推定する。図９Ｅに示すように、例えば、データ処理部７０は、１番のラベルのグリッドｇが最も多いため、１番のラベルの連続領域を地表面Ｆと推定する（白抜き部分参照）。データ処理部７０は、三次元地図Ｍをデータ表示部８０に表示した場合、地表面Ｆの領域を他の領域と色や濃淡を区別して表示する（図１２参照）。

　また、データ処理部７０は、手動で入力された荷物Ｗの範囲（図１１Ｂ参照）にあるグリッドｇの代表点ｐｒの中で最も標高値の低い代表点ｐｒのグリッドｇのラベルを取得し、荷物Ｗのグリッドｇに付与されたラベルとする。これは、グリッドｇの中に荷物Ｗのグリッドｇとフックやワイヤのグリッドｇが含まれる場合、フックやワイヤロープのグリッドｇにおける代表点ｐｒの標高値よりも荷物Ｗのグリッドｇにおける代表点ｐｒの標高値の方が低いためである。尚、荷物Ｗの高さ（Ｚ軸方向の長さ）は、ブーム９の先端部分からでは荷物Ｗの側面や下部を計測できないため、実測されてデータ処理部７０に入力されている。

（第二地表面推定工程）
　図１０に示すように、データ処理部７０は、グリッドｇに付与されたラベルを用いて、同一領域推定処理で推定された地表面Ｆと荷物Ｗ以外の各連続領域について、地表面Ｆまたは地物Ｅの推定を実施する（図１３の地物領域推定処理Ｄ（ＳＴＥＰ－１０７））。具体的には、地表面Ｆと推定された連続領域に隣接しない連続領域である非隣接連続領域Ｊａと、地表面Ｆと推定された連続領域に一部でも隣接する連続領域である隣接連続領域Ｊｂと、について実施する。ここで、「隣接する」とは、上下左右または斜めに隣り合うことをいう。

　まず、データ処理部７０は、推定を行う非隣接連続領域Ｊａにおけるグリッドｇの代表点ｐｒの平均標高値（以下、単に「非隣接連続領域Ｊａの平均標高値」という）を算出する。次に、データ処理部７０は、推定を行う非隣接連続領域Ｊａに対する所定の範囲Ｎａにおけるグリッドｇの代表点ｐｒの平均標高値（以下、単に「所定の範囲Ｎａの平均標高値」という）を算出する。所定の範囲Ｎａは、地表面Ｆと推定された連続領域内の範囲であるとともに、推定を行う非隣接連続領域Ｊａと、地表面Ｆと推定された連続領域と、の間にある隣接連続領域Ｊｂに隣接するグリッドｇで構成される範囲である（図１０Ａ参照）。但し、所定の範囲Ｎａは、この範囲に限定されず、地表面Ｆの状況に応じて任意の範囲に設定可能である。

　データ処理部７０は、非隣接連続領域Ｊａの平均標高値と所定の範囲Ｎａの平均標高値との差が領域間の閾値である領域閾値以下の場合、推定を行う非隣接連続領域Ｊａを地表面Ｆと推定する。データ処理部７０は、非隣接連続領域Ｊａの平均標高値と所定の範囲Ｎａの平均標高値の差が領域閾値よりも大きい場合、推定を行う非隣接連続領域Ｊａを地物Ｅと推定する。領域閾値は、グリッド閾値と同様にレーザスキャナ６２のＺ軸方向の計測誤差に基づいて設定されており、例えば、０．１ｍとしている。グリッド閾値と領域閾値とが個別に設定されるものとして説明するが、同一の閾値として設定されてもよい。

　また、データ処理部７０は、推定を行う隣接連続領域Ｊｂにおけるグリッドｇの代表点ｐｒの平均標高値（以下、単に「隣接連続領域Ｊｂの平均標高値」という）を算出する。次に、データ処理部７０は、推定を行う隣接連続領域Ｊｂに対する所定の範囲Ｎｂにおけるグリッドｇの代表点ｐｒの平均標高値（以下、単に「所定の範囲Ｎｂの平均標高値」という）を算出する。所定の範囲Ｎｂは、地表面Ｆと推定された連続領域内の範囲である。同時に、所定の範囲Ｎｂは、推定を行う隣接連続領域Ｊｂに隣接する他の隣接連続領域Ｊｂがない場合、隣接連続領域Ｊｂに隣接するグリッドｇで構成される範囲である（図１０Ａ参照）。所定の範囲Ｎｂは、推定を行う隣接連続領域Ｊｂに隣接する他の隣接連続領域Ｊｂがある場合、推定を行う隣接連続領域Ｊｂを含み、互いに隣接し合う複数の隣接連続領域Ｊｂを隣接領域Ｊｃとして、隣接領域Ｊｃに隣接するグリッドｇで構成される範囲である（図１０Ｄ参照）。但し、所定の範囲Ｎｂは、この範囲に限定されず、地表面Ｆの状況に応じて任意の範囲に設定可能である。

　データ処理部７０は、隣接連続領域Ｊｂの平均標高値と所定の範囲Ｎｂの平均標高値との差が領域閾値以下の場合、推定を行う隣接連続領域Ｊｂを地表面Ｆと推定し、隣接連続領域Ｊｂの平均標高値と所定の範囲Ｎｂの平均標高値との差が領域閾値よりも大きい場合、推定を行う隣接連続領域Ｊｂを地物Ｅと推定する。

　図１０Ａ及び図１０Ｂに示すように、例えば、同一領域推定処理において、１番のラベルの連続領域が地表面Ｆと推定されており、４番のラベルと６番のラベルとの連続領域が地表面Ｆまたは荷物Ｗと推定されていないものとする。６番のラベルの連続領域は、地表面Ｆと推定された連続領域に隣接していないため、非隣接連続領域Ｊａとなる。４番のラベルの連続領域は、地表面Ｆと推定された連続領域に隣接しているため、隣接連続領域Ｊｂとなる。６番のラベルの連続領域に対する所定の範囲Ｎａと４番のラベルの連続領域に対する所定の範囲Ｎｂとは、地表面Ｆと推定された連続領域において、４番のラベルの連続領域に隣接するグリッドｇで構成される範囲（図１０Ａ斜線部）となる。

　データ処理部７０は、６番のラベルの連続領域の平均標高値を、０．３ｍと算出している。データ処理部７０は、４番のラベルの連続領域の平均標高値を、１．７５ｍと算出している。データ処理部７０は、所定の範囲Ｎａと所定の範囲Ｎｂとの平均標高値を、０．３５ｍと算出している。図１０Ｃに示すように、データ処理部７０は、６番のラベルの連続領域の平均標高値と所定の範囲Ｎａの平均標高値との差が領域閾値０．１ｍより小さいため、６番のラベルの連続領域を地表面Ｆと推定する。データ処理部７０は、４番のラベルの連続領域の平均標高値と所定の範囲Ｎｂの平均標高値との差が領域閾値０．１ｍより大きいため、４番のラベルの連続領域を地物Ｅと推定する。

　最後に、データ処理部７０は、グリッドｇに付与されたラベルに基づいて、荷物Ｗ、地表面Ｆ、地物Ｅ毎に、グリッドｇの代表点ｐｒの集合を出力する。このように、データ処理部７０は、周囲の地表面Ｆから地物Ｅによって分断された領域を地表面Ｆとして推定できる。データ処理部７０は、このような地表面Ｆの推定をリアルタイムで自動的に行うことができる。

　図１０Ｄに示すように、隣接連続領域Ｊｂの少なくとも一部が地表面Ｆと推定された連続領域に隣接していない場合の所定の範囲Ｎｂの一例について説明する。例えば、７番のラベルと８番のラベルとの連続領域が地表面Ｆまたは荷物Ｗと推定されていないものとする。７番のラベルと８番のラベルとの連続領域は、地表面Ｆと推定された連続領域に隣接しているため、隣接連続領域Ｊｂとなる。但し、８番のラベルの連続領域における２つのグリッドｇが地表面Ｆと推定された連続領域に隣接していない。７番のラベルの連続領域と８番のラベルの連続領域とに対する所定の範囲Ｎｂは、地表面Ｆと推定された連続領域において、７番のラベルと８番のラベルとの連続領域に隣接するグリッドｇで構成される範囲（図１０Ｄ斜線部）となる。

　このような地表面推定方法を実施する計測領域表示システム５０は、非隣接連続領域Ｊａの平均標高値と所定の範囲Ｎａの平均標高値との差が閾値の範囲内の差であれば、非隣接連続領域Ｊａが地表面Ｆと推定される。また、計測領域表示システム５０は、隣接連続領域Ｊｂの平均標高値と所定の範囲Ｎｂの平均標高値との差が閾値の範囲内の差であれば、隣接連続領域Ｊｂが地表面Ｆと推定される。所定の範囲Ｎａと所定の範囲Ｎｂとを上述の実施形態の範囲にすると、地表面Ｆと推定されている領域の中で非隣接連続領域Ｊａまたは隣接連続領域Ｊｂに最も近接している部分の標高値で判定するので、標高値の差が大きい地表面Ｆの影響を受けない。これにより、他の領域によって分断されている領域についても地表面Ｆの推定を行える。

　本実施形態において、計測領域表示システム５０がクレーン１から独立した構成として説明したが、データ取得部６０が備える各種センサ６１～６４ｄ、データ表示部８０が備える表示装置、データ入力部９０が備える入力装置をクレーン１が具備し、データ処理部７０が行うデータ処理を制御装置１８が行う構成であってもよい。

　このようなクレーン１においては、旋回台７やブーム９の動作に伴ってクレーン１の作業領域ＳＡにおける点群データＰが取得され、非隣接連続領域Ｊａの平均標高値と所定の範囲Ｎａの平均標高値との差が閾値の範囲内の差であれば、非隣接連続領域Ｊａが地表面Ｆと推定される。また、クレーン１においては、隣接連続領域Ｊｂの平均標高値と所定の範囲Ｎｂの平均標高値との差が閾値の範囲内の差であれば、隣接連続領域Ｊｂが地表面Ｆと推定される。所定の範囲Ｎａと所定の範囲Ｎｂとを上述の実施形態の範囲にすると、地表面Ｆと推定されている領域の中で非隣接連続領域Ｊａまたは隣接連続領域Ｊｂに最も近接している部分の標高値で判定するので、標高値の差が大きい地表面Ｆの影響を受けない。これにより、クレーン１の作業領域ＳＡにおいて、他の領域によって分断されている領域についても地表面Ｆの推定を行える。尚、計測領域表示システム５０は、クレーン１以外のその他の作業車両に用いることもできる。例えば、高所作業車等にも用いることができる。

　データ処理部７０はグリッドｇの代表点ｐｒの集合を用いて、荷物Ｗ、地物Ｅ毎に大きさと高さをカメラ画像上に可視化する（図１３の領域可視化処理（ＳＴＥＰ－１０８））。具体的には、まず、データ処理部７０は、第一ＧＮＳＳ受信機６４ａと第二ＧＮＳＳ受信機６４ｂとの計測データから現在時刻のレーザスキャナ６２の位置と向き（方位）とを算出する。次に、データ処理部７０は、算出したレーザスキャナ６２の位置と向きとに基づいて座標軸を決定し、荷物Ｗや地物Ｅ等の代表点ｐｒのカメラ画像上における位置を合わせて表示する。

　図１１に示すように、データ処理部７０は、荷物Ｗや地物Ｅ等の代表点ｐｒのカメラ画像Ｉ上における位置を用いて、荷物Ｗや地物Ｅをカメラ画像Ｉ上で内包する外形線を算出し、外形線が成す外形図形をガイド情報ＧＤ１としてカメラ画像Ｉ上に表示する。最後に、データ処理部７０は、ガイド情報ＧＤ１の上に荷物Ｗや地物Ｅの標高値と地表面Ｆの標高値の差分値をガイド情報ＧＤ２（高さの情報）として表示する。尚、ガイド情報ＧＤ１の描画色は、荷物Ｗや地物Ｅの高さに応じて濃淡を変更した所定の色で表示する。また、荷物Ｗや地物Ｅの高さの有効桁数は、小数第一位とする。

　図１２に示すように、データ処理部７０は、荷物Ｗ、地表面Ｆ、地物Ｅ毎にラベルが付与された計測領域ＫＡの三次元地図を用いて、それぞれの三次元空間における位置関係と大きさを可視化する（図１３の三次元地図可視化処理（ＳＴＥＰ－１０９））。具体的には、まず、データ処理部７０は、荷物Ｗ、地表面Ｆ、地物Ｅ毎の代表点ｐｒの位置と標高値とに基づいて、代表点ｐｒを重心とする面Ｓを生成する。このとき、面Ｓの幅は代表点ｐｒの生成時に用いたグリッドｇの幅とする。そして、データ処理部７０は、荷物Ｗ、地表面Ｆ、地物Ｅ毎に面Ｓに色付けを行い、三次元空間における三次元地図Ｍとしてデータ表示部８０に可視化する。

　尚、クレーン１や計測領域ＫＡの状況に合わせてグリッドｇの幅が変更されて、三次元地図Ｍの可視性が調整される。例えば、データ処理部７０は、ブーム９の先端部分の高さが高くなるに応じてグリッドｇの幅を小さくすると、高い位置から計測された広範囲の計測領域ＫＡを、より細かな面Ｓの三次元地図Ｍで可視化できる。データ処理部７０は、計測領域ＫＡにおいて高さの差がグリッド閾値以内である複数の地物Ｅが互いに近接している場合、グリッドｇの幅を大きくすると、複数の地物Ｅを一つの地物Ｅとして三次元地図Ｍに可視化できる。

　以下に、計測領域表示システム５０の制御態様について具体的に説明する。

　図１３に示すように、データ処理部７０は、計測データが入力され、ＳＴＥＰ－１０１に移行する。

　データ処理部７０は、ＳＴＥＰ－１０１において、時刻対応処理を行い、ＳＴＥＰ－１０２に移行する。

　データ処理部７０は、ＳＴＥＰ－１０２において、剛体変換処理を行い、ＳＴＥＰ－１０３に移行する。

　データ処理部７０は、ＳＴＥＰ－１０３において、グリッド生成処理Ａを開始し、ＳＴＥＰ－２０１に移行する（図１４参照）。そして、グリッド生成処理Ａが終了するとＳＴＥＰ－１０４に移行する。

　データ処理部７０は、ＳＴＥＰ－１０４において、三次元地図更新処理Ｂを開始し、ＳＴＥＰ－３０１に移行する（図１５参照）。そして、三次元地図更新処理Ｂが終了するとＳＴＥＰ－１０５に移行する。尚、データ処理部７０は、現在時刻よりも前に三次元地図が生成されていない場合は、新に三次元地図を生成する。

　データ処理部７０は、ＳＴＥＰ－１０５において、ラベリング処理Ｃを開始し、ＳＴＥＰ－４０１に移行する（図１６参照）。そして、ラベリング処理Ｃが終了するとＳＴＥＰ－１０６に移行する。

　データ処理部７０は、ＳＴＥＰ－１０６において、同一領域推定処理を行い、ＳＴＥＰ－１０７に移行する。

　データ処理部７０は、ＳＴＥＰ－１０７において、地物領域推定処理Ｄを開始し、ＳＴＥＰ－５０１に移行する（図１７参照）。そして、地物領域推定処理Ｄが終了するとＳＴＥＰ－１０８に移行する。

　データ処理部７０は、ＳＴＥＰ－１０８において、領域可視化処理を行い、ＳＴＥＰ－１０９に移行する。

　データ処理部７０は、ＳＴＥＰ－１０９において、三次元地図可視化処理を行い、可視化結果を出力する。可視化結果とは、三次元地図Ｍやカメラ画像情報等である。

　図１４に示すように、データ処理部７０は、ＳＴＥＰ－２０１において、グリッド生成処理Ａを開始し、補正された点群データＰを重ね合わせて、ＳＴＥＰ－２０２に移行する。

　データ処理部７０は、ＳＴＥＰ－２０２において、重ね合わせた点群データＰを一定の大きさのグリッドｇで分割し、ＳＴＥＰ－２０３に移行する。

　データ処理部７０は、ＳＴＥＰ－２０３において、点群データＰが存在するグリッドｇを任意の間隔で水平の方向にスライスし、複数の層Ｌ１・Ｌ２・Ｌ３・・・を形成し、ＳＴＥＰ－２０４に移行する。

　データ処理部７０は、ＳＴＥＰ－２０４において、各層のうち、点データｐが存在し、かつ、最も標高が高い位置にある層Ｌ１の点データｐの平均標高値Ｈを算出し、ＳＴＥＰ－２０５に移行する。

　データ処理部７０は、ＳＴＥＰ－２０５において、グリッドｇの重心位置を算出し、ＳＴＥＰ－２０６に移行する。

　データ処理部７０は、ＳＴＥＰ－２０６において、算出した平均標高値Ｈにおけるグリッドｇの重心位置をグリッドｇの代表点ｐｒとして設定し、ＳＴＥＰ－２０７に移行する。

　データ処理部７０は、ＳＴＥＰ－２０７において、点群データＰが存在する全てのグリッドｇで代表点ｐｒの設定が完了したか否かを判定する。
　その結果、データ処理部７０は、点群データＰが存在する全てのグリッドｇで代表点ｐｒの設定が完了した場合、グリッド生成処理Ａを終了し、ＳＴＥＰ－１０４に移行する。
　一方、データ処理部７０は、点群データＰが存在する全てのグリッドｇで代表点ｐｒの設定が完了していない場合、ＳＴＥＰ－２０３に移行して、代表点ｐｒの設定が行われていないグリッドｇに対してＳＴＥＰ－２０３を行う。

　図１５に示すように、データ処理部７０は、ＳＴＥＰ－３０１において、三次元地図更新処理Ｂを開始し、グリッドｇにおける、設定した代表点ｐｒと、現在時刻より前に生成されている三次元地図の代表点ｐｒと、を取得する。

　データ処理部７０は、ＳＴＥＰ－３０２において、同一のグリッドｇ内に複数の代表点ｐｒが存在するか否かを判定する。
　その結果、データ処理部７０は、同一のグリッドｇ内に複数の代表点ｐｒが存在する場合、ＳＴＥＰ－３０３に移行する。
　一方、データ処理部７０は、同一のグリッドｇ内に複数の代表点ｐｒが存在しない場合、ＳＴＥＰ－３０４に移行する。

　データ処理部７０は、ＳＴＥＰ－３０３において、現在時刻よりも前に生成されている三次元地図の代表点ｐｒを削除し、ＳＴＥＰ－３０４に移行する。

　データ処理部７０は、ＳＴＥＰ－３０４において、全てのグリッドｇでＳＴＥＰ－３０１、３０２が完了したか否かを判定する。
　その結果、データ処理部７０は、全てのグリッドｇでＳＴＥＰ－３０１、３０２が完了した場合、ＳＴＥＰ－３０５に移行する。
　一方、データ処理部７０は、全てのグリッドｇでＳＴＥＰ－３０１、３０２が完了していない場合、ＳＴＥＰ－３０１に移行して、ＳＴＥＰ－３０１、３０２が行われていないグリッドｇに対してＳＴＥＰ－３０１を行う。

　データ処理部７０は、ＳＴＥＰ－３０５において、グリッドｇと同じ大きさの面を生成し、計測領域ＫＡの三次元地図を出力して、三次元地図更新処理Ｂを終了し、ＳＴＥＰ－１０５に移行する。

　図１６に示すように、データ処理部７０は、ＳＴＥＰ－４０１において、ラベリング処理Ｃを開始し、注目グリッドｇａの代表点ｐｒの標高値とその４近傍にある各グリッドｇの代表点ｐｒの標高値との差を算出し、ＳＴＥＰ－４０２に移行する。

　データ処理部７０は、ＳＴＥＰ－４０２において、算出した標高値の差のいずれかがグリッド閾値以下か否かを判定する。尚、データ処理部７０は、４近傍のグリッドｇがない（注目グリッドｇａが三次元地図の左上にある）場合は、ＳＴＥＰ－４０４に移行して注目グリッドｇａに１番のラベルを付与する。
　その結果、データ処理部７０は、算出した標高値の差のいずれかがグリッド閾値以下である場合、ＳＴＥＰ－４０３に移行する。
　一方、データ処理部７０は、算出した標高値の差が全てグリッド閾値よりも大きい場合、ＳＴＥＰ－４０４に移行する。

　データ処理部７０は、ＳＴＥＰ－４０３において、標高値の差がグリッド閾値以下のグリッドｇの中で最も番号が小さいラベルを注目グリッドｇａに付与し、ＳＴＥＰ－４０５に移行する。

　データ処理部７０は、ＳＴＥＰ－４０４において、注目グリッドｇａに新たなラベルを付与し、ＳＴＥＰ－４０５に移行する。

　データ処理部７０は、ＳＴＥＰ－４０５において、全てのグリッドｇに対してラベルの付与が完了したか否かを判定する。
　その結果、データ処理部７０は、全てのグリッドｇに対してラベルの付与が完了した場合、ＳＴＥＰ－４０６に移行する。
　一方、データ処理部７０は、全てのグリッドｇに対してラベルの付与が完了していない場合、ＳＴＥＰ－４０１に移行し、注目グリッドｇａの右隣または次の行のグリッドｇを次の注目グリッドｇａとしてＳＴＥＰ－４０１を行う。

　データ処理部７０は、ＳＴＥＰ－４０６において、注目グリッドｇａの代表点ｐｒの標高値とその４近傍にある各グリッドｇの代表点ｐｒの標高値との差を算出し、ＳＴＥＰ－４０７に移行する。

　データ処理部７０は、ＳＴＥＰ－４０７において、注目グリッドｇａの４近傍のグリッドｇに、標高値の差がグリッド閾値以下で、かつ、注目グリッドｇａと異なるラベルが付与されているか否かを判定する。
　その結果、データ処理部７０は、注目グリッドｇａの４近傍のグリッドｇに、標高値の差がグリッド閾値以下で、かつ、注目グリッドｇａと異なるラベルが付与されている場合、ＳＴＥＰ－４０８に移行する。
　一方、データ処理部７０は、注目グリッドｇａの４近傍のグリッドｇに、標高値の差がグリッド閾値以下で、かつ、注目グリッドｇａと異なるラベルが付与されていない場合、ＳＴＥＰ－４０９に移行する。

　データ処理部７０は、ＳＴＥＰ－４０８において、標高値の差がグリッド閾値以下で、かつ、注目グリッドｇａと異なるラベルが付与されているグリッドｇのラベルを注目グリッドｇａのラベルで上書きし、ＳＴＥＰ－４０９に移行する。

　データ処理部７０は、ＳＴＥＰ－４０９において、全てのグリッドｇに対してＳＴＥＰ－４０６、４０７の処理が完了したか否かを判定する。
　その結果、データ処理部７０は、全てのグリッドｇに対してＳＴＥＰ－４０６、４０７の処理が完了した場合、ＳＴＥＰ－４１０に移行する。
　一方、データ処理部７０は、全てのグリッドｇに対してＳＴＥＰ－４０６、４０７の処理が完了していない場合、ＳＴＥＰ－４０６に移行し、注目グリッドｇａの右隣または次の行のグリッドｇを次の注目グリッドｇａとしてＳＴＥＰ－４０６を行う。

　データ処理部７０は、ＳＴＥＰ－４１０において、グリッドｇにラベルが付与された三次元地図を出力して、ラベリング処理Ｃ処理を終了し、ＳＴＥＰ－１０６に移行する。

　図１７に示すように、データ処理部７０は、ＳＴＥＰ－５０１において、地物領域推定処理Ｄを開始し、推定を行う連続領域が非隣接連続領域Ｊａか否かを判定する。
　その結果、データ処理部７０は、推定を行う連続領域が非隣接連続領域Ｊａである場合、ＳＴＥＰ－５０２に移行する。
　一方、データ処理部７０は、推定を行う連続領域が非隣接連続領域Ｊａでない場合、ＳＴＥＰ－５０６に移行する。

　データ処理部７０は、ＳＴＥＰ－５０２において、非隣接連続領域Ｊａの平均標高値と所定の範囲Ｎａの平均標高値とを算出し、ＳＴＥＰ－５０３に移行する。

　データ処理部７０は、ＳＴＥＰ－５０３において、非隣接連続領域Ｊａの平均標高値と所定の範囲Ｎａの平均標高値との差が領域閾値以下か否かを判定する。
　その結果、データ処理部７０は、非隣接連続領域Ｊａの平均標高値と所定の範囲Ｎａの平均標高値との差が領域閾値以下である場合、ＳＴＥＰ－５０４に移行する。
　一方、データ処理部７０は、非隣接連続領域Ｊａの平均標高値と所定の範囲Ｎａの平均標高値との差が領域閾値よりも大きい場合、ＳＴＥＰ－５０５に移行する。

　データ処理部７０は、ＳＴＥＰ－５０４において、推定を行う非隣接連続領域Ｊａを地表面Ｆと推定し、ＳＴＥＰ－５１０に移行する。

　データ処理部７０は、ＳＴＥＰ－５０５において、推定を行う非隣接連続領域Ｊａを地物Ｅと推定し、ＳＴＥＰ－５１０に移行する。

　データ処理部７０は、ＳＴＥＰ－５０６において、隣接連続領域Ｊｂの平均標高値と所定の範囲Ｎｂの平均標高値とを算出し、ＳＴＥＰ－５０７に移行する。

　データ処理部７０は、ＳＴＥＰ－５０７において、隣接連続領域Ｊｂの平均標高値と所定の範囲Ｎｂの平均標高値との差が領域閾値以下か否かを判定する。
　その結果、データ処理部７０は、隣接連続領域Ｊｂの平均標高値と所定の範囲Ｎｂの平均標高値との差が領域閾値以下である場合、ＳＴＥＰ－５０８に移行する。
　一方、データ処理部７０は、隣接連続領域Ｊｂの平均標高値と所定の範囲Ｎｂの平均標高値との差が領域閾値よりも大きい場合、ＳＴＥＰ－５０９に移行する。

　データ処理部７０は、ＳＴＥＰ－５０８において、推定を行う隣接連続領域Ｊｂを地表面Ｆと推定し、ＳＴＥＰ－５１０に移行する。

　データ処理部７０は、ＳＴＥＰ－５０９において、推定を行う隣接連続領域Ｊｂを地物Ｅと推定し、ＳＴＥＰ－５１０に移行する。

　データ処理部７０は、ＳＴＥＰ－５１０において、地表面Ｆと荷物Ｗ以外の全て連続領域で地表面Ｆまたは地物Ｅの推定が完了したか否かを判定する。
　その結果、データ処理部７０は、地表面Ｆと荷物Ｗ以外の全て連続領域で地表面Ｆまたは地物Ｅの推定が完了した場合、ＳＴＥＰ－５１１に移行する。
　一方、データ処理部７０は、地表面Ｆと荷物Ｗ以外の全て連続領域で地表面Ｆまたは地物Ｅの推定が完了していない場合、ＳＴＥＰ－５０１に移行し、地表面Ｆまたは地物Ｅの推定が行われていない連続領域に対してＳＴＥＰ－５０１を行う。

　データ処理部７０は、ＳＴＥＰ－５１１において、荷物Ｗ、地表面Ｆ、地物Ｅ毎に、グリッドｇの代表点ｐｒの集合を出力し、地物領域推定処理Ｄを終了し、ＳＴＥＰ－１０８に移行する。

　以下に、非隣接連続領域Ｊａの所定の範囲内における、連続領域の数に応じて計測領域ＫＡを複数のグリッドｇに分割する際のグリッドｇの幅を変更する地表面推定方法について説明する。

　図１８に示すように、データ処理部７０は、非隣接連続領域Ｊａの所定の範囲Ｎｃにおける、連続領域の数に応じて計測領域ＫＡを複数のグリッドｇに分割する際のグリッドｇの幅を変更する。具体的には、まず、データ処理部７０は、非隣接連続領域Ｊａのグリッドｇ上で所定の範囲Ｎｃを移動させていく。所定の範囲Ｎｃは、地物Ｅの状況に応じて設定されており、例えば、縦横が２グリッドの長さの正方形の範囲としている。次に、データ処理部７０は、所定の範囲Ｎｃの連続領域の数に応じてグリッドｇの幅を変更して、再びラベリング処理Ｃを行う。グリッドｇの幅の変化率は、任意に設定されている。

　このように、データ処理部７０は、非隣接連続領域Ｊａの所定の範囲Ｎｃにおける、連続領域の数に応じてグリッドｇの幅を変更する。これは、非隣接連続領域Ｊａの所定の範囲Ｎｃにおける、連続領域の数が多くなるに応じて地物Ｅ等の領域が密集し、より詳細に非隣接連続領域Ｊａを認識する必要がある可能性があるためである。

　例えば、データ処理部７０が非隣接連続領域Ｊａのグリッドｇ上で所定の範囲Ｎｃを移動させていくと、所定の範囲Ｎｃに３番から５番のラベルの連続領域があるとする。データ処理部７０は、所定の範囲Ｎｃに連続領域が３つある場合、グリッドｇの幅を半分に減少させる設定がされているものとする。データ処理部７０は、所定の範囲Ｎｃに３つの連続領域があるため、グリッドｇの幅を半分に減少させて、再びラベリング処理Ｃを行っている。図１８Ｂに示すように、データ処理部７０は、グリッドｇの幅を減少させることによって、３番のラベルと４番のラベルとの連続領域をより詳細に認識している。尚、データ処理部７０は、グリッドｇの幅を小さくしても連続領域の形状が変わらなかった場合、小さくする前のグリッドｇの幅に戻してもよい。

　図１９に示すように、データ処理部７０は、ＳＴＥＰ－２０８において、非隣接連続領域Ｊａの所定の範囲Ｎｃにおける、連続領域の数に応じてグリッドｇの幅を変更し、ＳＴＥＰ－２０２に移行する。

　このような地表面推定方法は、複数の非隣接連続領域Ｊａが密集する場合、計測領域ＫＡを複数のグリッドｇに分割する際のグリッドｇの幅を小さくすることによって、密集した非隣接連続領域Ｊａがより詳細に認識される。これにより、他の領域によって分断されている領域について地表面Ｆの推定を行える。

　以下に、非隣接連続領域Ｊａと、地表面Ｆと推定された連続領域と、の間の距離（以下、領域間距離Ｌａという）が所定の距離以上であれば、非隣接連続領域Ｊａを地表面Ｆと推定しない地表面推定方法について説明する。

　図２０に示すように、データ処理部７０は、領域間距離Ｌａが所定の距離Ｌｂ以上であれば、非隣接連続領域Ｊａを地表面Ｆと推定しない。具体的には、まず、データ処理部７０は、領域間距離Ｌａを算出する。次に、データ処理部７０は、領域間距離Ｌａが所定の距離Ｌｂ以上の場合、非隣接連続領域Ｊａを地表面Ｆまたは地物Ｅと推定せず、領域間距離Ｌａが所定の距離Ｌｂよりも小さい場合、非隣接連続領域Ｊａを地表面Ｆまたは地物Ｅと推定する。所定の距離Ｌｂは、地表面Ｆの状況に応じて設定されており、例えば、２グリッドの距離としている。領域間距離Ｌａは、非隣接連続領域Ｊａと、地表面Ｆと推定された連続領域と、の間の最短距離とするが、任意の地点間の距離でもよい。そして、データ処理部７０は、地表面Ｆまたは地物Ｅと推定しない非隣接連続領域Ｊａをガイド情報として表示する場合、地表面Ｆまたは地物Ｅと推定されていないことをオペレータが認識できるようにデータ表示部８０に表示する。

　このように、データ処理部７０は、領域間距離Ｌａが所定の距離Ｌｂ以上であれば、非隣接連続領域Ｊａを地表面Ｆと推定しない。これは、領域間距離Ｌａが長くなるに応じて非隣接連続領域Ｊａを地表面Ｆまたは地物Ｅと推定する判定の不確実性が高くなるためである。具体的に説明すると、データ処理部７０は、地表面Ｆの領域を分断している領域（例えば地物Ｅの領域）に存在する地表面Ｆの標高値の変動幅が領域閾値の範囲内になることを前提として、非隣接連続領域Ｊａを地表面Ｆまたは地物Ｅと推定する判定を行っている。しかし、地表面Ｆの領域を分断している領域には、地表面Ｆの傾斜や凹凸がある可能性がある。このため、領域間距離Ｌａが長くなるに応じて地表面Ｆの領域を分断している領域にある地表面Ｆの標高値の変動幅が領域閾値の範囲内となることの不確実性が高くなる。つまり、データ処理部７０は、領域間距離Ｌａが長くなるに応じて非隣接連続領域Ｊａを地表面Ｆまたは地物Ｅと推定する判定の不確実性も高くなるため、領域間距離Ｌａが所定の距離以上の場合、非隣接連続領域Ｊａを地表面Ｆまたは地物Ｅと推定しない。

　例えば、領域間距離Ｌａが３グリッドの距離とする。データ処理部７０は、領域間距離Ｌａが所定の距離Ｌｂ以上であるため、３番のラベルの連続領域を地表面Ｆまたは地物Ｅと推定しない。尚、領域間距離Ｌａが長くなるに応じて領域閾値を大きく設定するようにしてもよい。これは、地表面Ｆが傾斜して地表面Ｆと推定された連続領域から非隣接連続領域Ｊａに続いている場合、領域間距離Ｌａが長くなるに応じて、地表面Ｆと推定された連続領域と非隣接連続領域Ｊａとにおける、地表面Ｆの標高値の差が大きくなるためである。

　図２１に示すように、データ処理部７０は、ＳＴＥＰ－５１２において、領域間距離Ｌａを算出し、ＳＴＥＰ－５１３に移行する。

　データ処理部７０は、ＳＴＥＰ－５１３において、領域間距離Ｌａが所定の距離Ｌｂ以上か否かを判定する。
　その結果、データ処理部７０は、領域間距離Ｌａが所定の距離Ｌｂ以上である場合、ＳＴＥＰ－５１０に移行する。
　一方、データ処理部７０は、領域間距離Ｌａが所定の距離Ｌｂより小さい場合、ＳＴＥＰ－５０２に移行する。

　このような地表面推定方法は、地表面Ｆと推定された連続領域からの距離が長くなるに応じて分断している領域にある地表面Ｆの標高値が閾値の範囲内となることの不確実性が高くなるため、非隣接連続領域Ｊａと、地表面Ｆと推定された連続領域と、の間の距離が所定の距離Ｌｂ以上の非隣接連続領域Ｊａを地表面Ｆと推定しないようにすることで、誤った推定が抑制される。これにより、他の領域によって分断されている領域について地表面Ｆの推定を行える。

　上述の実施形態は、代表的な形態を示したに過ぎず、一実施形態の骨子を逸脱しない範囲で種々変形して実施することができる。さらに種々なる形態で実施し得ることは勿論のことであり、本発明の範囲は、特許請求の範囲の記載によって示され、さらに特許請求の範囲に記載の均等の意味、および範囲内のすべての変更を含む。

　本発明は、地表面推定方法、計測領域表示システムおよびクレーンに利用可能である。

　１　　　クレーン
　６１　　カメラ
　６２　　レーザスキャナ
　７０　　データ処理部
　Ｆ　　　地表面
　ｇ　　　グリッド
　Ｈ　　　平均標高値
　Ｊａ　　非隣接連続領域
　ＫＡ　　計測領域
　ｐ　　　点データ
　Ｐ　　　点群データ
　ｐｒ　　代表点

Claims

　レーザスキャナによって、計測領域にレーザを照射して点群データを取得し、前記点群データから照射点の標高値を算出する点群データ取得工程と、
　前記点群データを演算処理するデータ処理手段によって、
　前記計測領域を複数のグリッドに分割し、前記グリッド毎に、当該グリッドの重心位置と当該グリッド内の前記点群データの平均標高値とを算出し、前記グリッド内の点群データの平均標高値における前記重心位置を前記グリッド毎の代表点の位置として設定するグリッド生成処理工程と、
　前記複数のグリッドのうち、一のグリッドの代表点と隣接する他のグリッドの代表点との標高値の差がグリッド間の閾値以下であれば、前記一のグリッドと前記他のグリッドとが連続する領域である連続領域として認識する連続領域認識工程と、
　前記連続領域のうち、最も前記グリッドの数が多い連続領域を地表面と推定する第一地表面推定工程と、
　前記地表面と推定された連続領域に隣接していない連続領域である非隣接連続領域毎に、前記非隣接連続領域におけるグリッドの代表点の平均標高値を算出し、前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値を算出し、前記非隣接連続領域におけるグリッドの代表点の平均標高値と前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値との差が領域間の閾値以下である前記非隣接連続領域を地表面と推定する第二地表面推定工程と、
　を備える、ことを特徴とする地表面推定方法。
　前記地表面と推定された連続領域の所定の範囲は、前記非隣接連続領域と、前記地表面と推定された連続領域と、の間にある連続領域に隣接するグリッドから構成される範囲である、ことを特徴とする請求項１に記載の地表面推定方法。
　前記地表面と推定された非隣接連続領域の所定の範囲における、前記連続領域の数に応じて前記計測領域を複数のグリッドに分割する際の前記グリッドの幅を変更する、ことを特徴とする請求項１または請求項２に記載の地表面推定方法。
　前記第二地表面推定工程は、
　前記非隣接連続領域と、前記地表面と推定された連続領域と、の間の距離が所定の距離以上であれば、前記非隣接連続領域を地表面と推定しない、ことを特徴とする請求項１から請求項３のいずれか一項に記載の地表面推定方法。
　計測領域にレーザを照射して点群データを取得するレーザスキャナを備えたデータ取得部と、
　取得した点群データを演算処理するデータ処理部と、
　データ表示部と、を備え、
　前記データ処理部は、
　前記データ取得部から前記点群データを取得し、前記点群データから照射点の標高値を算出し、
　前記計測領域を複数のグリッドに分割し、前記グリッド毎に、当該グリッドの重心位置と当該グリッド内の前記点群データの平均標高値とを算出し、前記グリッド内の点群データの平均標高値における前記重心位置を前記グリッド毎の代表点の位置として設定し、
　前記複数のグリッドのうち、一のグリッドの代表点と隣接する他のグリッドの代表点との標高値の差がグリッド間の閾値以下であれば、前記一のグリッドと前記他のグリッドとが連続する領域である連続領域として認識し、
　前記連続領域のうち、最も前記グリッドの数が多い連続領域を地表面と推定し、
　前記地表面と推定された連続領域に隣接していない連続領域である非隣接連続領域毎に、前記非隣接連続領域におけるグリッドの代表点の平均標高値を算出し、前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値を算出し、前記非隣接連続領域におけるグリッドの代表点の平均標高値と前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値との差が領域間の閾値以下である前記非隣接連続領域を地表面と推定し、
　前記地表面と推定された連続領域を他の領域と区別して前記データ表示部に表示する、ことを特徴とする計測領域表示システム。
　旋回台と、
　前記旋回台に設けられるブームと、
　前記ブームに取り付けられ、点群データを取得するレーザスキャナと、
　取得した点群データを演算処理する制御装置と、
　表示装置と、を備えるクレーンにおいて、
　前記旋回台の旋回操作と、前記ブームの伸縮操作と起伏操作と、に伴って前記レーザスキャナを移動させながらレーザを照射させることで、前記レーザスキャナのレーザ照射時の位置毎の点群データを取得し、
　前記制御装置は、
　前記レーザスキャナのレーザ照射時の位置毎の前記点群データを取得し、前記レーザスキャナのレーザ照射時の位置とその姿勢とに基づいてレーザ照射時の位置毎の前記点群データを重ね合わせ、照射点の標高値を算出し、
　前記レーザスキャナの計測領域を複数のグリッドに分割し、前記グリッド毎に、当該グリッドの重心位置と当該グリッド内の前記点群データの平均標高値とを算出し、前記グリッド内の点群データの平均標高値における前記重心位置を前記グリッド毎の代表点の位置として設定し、
　前記複数のグリッドのうち、一のグリッドの代表点と隣接する他のグリッドの代表点との標高値の差がグリッド間の閾値以下であれば、前記一のグリッドと前記他のグリッドとが連続する領域である連続領域として認識し、
　前記連続領域のうち、最も前記グリッドの数が多い連続領域を地表面と推定し、
　前記地表面と推定された連続領域に隣接していない連続領域である非隣接連続領域毎に、前記非隣接連続領域におけるグリッドの代表点の平均標高値を算出し、前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値を算出し、前記非隣接連続領域におけるグリッドの代表点の平均標高値と前記地表面と推定された連続領域の所定の範囲におけるグリッドの代表点の平均標高値との差が領域間の閾値以下である前記非隣接連続領域を地表面と推定し、
　前記地表面と推定された連続領域を他の領域と区別して前記表示装置に表示する、ことを特徴とするクレーン。