特許 5671184 断層図作成装置及び断層図作成方法並びに断層図作成プログラム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B1)

(11)【特許番号】5671184

(24)【登録日】2014年12月26日

(45)【発行日】2015年2月18日

(54)【発明の名称】断層図作成装置及び断層図作成方法並びに断層図作成プログラム

(51)【国際特許分類】

   G09B 29/00 20060101AFI20150129BHJP

   G01V 5/02 20060101ALI20150129BHJP

   G06T 17/05 20110101ALI20150129BHJP

【ＦＩ】

   G09B29/00 Z

   G01V5/02

   G06T17/05

【請求項の数】21

【全頁数】57

(21)【出願番号】特願2014-187849(P2014-187849)

(22)【出願日】2014年9月16日

【審査請求日】2014年9月17日

(31)【優先権主張番号】特願2013-230313(P2013-230313)

(32)【優先日】2013年11月6日

(33)【優先権主張国】JP

【新規性喪失の例外の表示】特許法第３０条第２項適用 ２０１３年１０月２４日〜２５日、日本応用地質学会 平成２５年度研究発表会 名古屋大学野依記念学術交流館にて、平川泰之及び吉村達朗が公開

【早期審査対象出願】

(73)【特許権者】

【識別番号】591074161

【氏名又は名称】アジア航測株式会社

(73)【特許権者】

【識別番号】393003505

【氏名又は名称】復建調査設計株式会社

(74)【代理人】

【識別番号】100083806

【弁理士】

【氏名又は名称】三好 秀和

(74)【代理人】

【識別番号】100100712

【弁理士】

【氏名又は名称】岩▲崎▼ 幸邦

(74)【代理人】

【識別番号】100101247

【弁理士】

【氏名又は名称】高橋 俊一

(74)【代理人】

【識別番号】100095500

【弁理士】

【氏名又は名称】伊藤 正和

(74)【代理人】

【識別番号】100098327

【弁理士】

【氏名又は名称】高松 俊雄

(72)【発明者】

【氏名】平川 泰之

(72)【発明者】

【氏名】吉村 辰朗

【審査官】 有家 秀郎

(56)【参考文献】

【文献】 特開２０１２−０５３１９１（ＪＰ，Ａ）

【文献】 特開２０１１−１３３９５２（ＪＰ，Ａ）

【文献】 特開２０１１−０４８４９５（ＪＰ，Ａ）

【文献】 特開平０８−２４８１４６（ＪＰ，Ａ）

【文献】 国際公開第２０１３／１８７５０２（ＷＯ，Ａ１）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０９Ｂ ２９／００

Ｇ０１Ｖ ５／００− ５／１４

Ｇ０６Ｔ １７／０５

(57)【特許請求の範囲】

【請求項1】

微細断層線ＤＬｊを検出可能な間隔で航空機から地域にレーザパルスを発射して得た数値標高モデルが記憶された第１の記憶手段と、
赤色立体化画像ＲＧｉが記憶される第２の記憶手段と、
前記赤色立体化画像ＲＧｉに基づいて決定した前記地域におけるγ線測定点Ｐｉにおける複数測定ポイントのγ線強度と測定間隔とに基づいて求められた前記微細断層線ＤＬｊの断層長Ｌ及び前記γ線測定点Ｐｉでの前記微細断層線ＤＬｊの方向が現地γ線測定情報Ｆγｉとして記憶された第３の記憶手段と
前記微細断層線ＤＬｊの断層帯幅ＦＷを検出可能な検索範囲を設定するための入力ボックスを表示部の画面に表示する断層用検出範囲設定部と、
前記第１の記憶手段の前記数値標高モデルと前記断層用検出範囲設定部によって設定された検索範囲に基づく前記赤色立体化画像ＲＧｉを生成して前記第２の記憶手段に記憶して前記画面に表示する赤色立体化画像作成部と、
前記画面の前記赤色立体化画像ＲＧｉ上に、前記微細断層線ＤＬｊを描かせるための面を表示し、この面におけるカーソルの移動軌跡を前記微細断層線ＤＬｊとし、これを所定の色で表示する手段と、
前記微細断層線ＤＬｊを描かせるための面の上に、前記γ線測定点Ｐｉを指定させるための面を表示し、この面に指定された点を前記γ線測定点Ｐｉとして前記微細断層線ＤＬｊとは異なる色で表示させる手段と、
前記微細断層線ＤＬｊを描かせるための面の上又は下にγ線現地測定結果直線ＬＧｊを表示するための面を表示し、この面に前記γ線測定点Ｐｉを指定させるための面に定義された前記γ線測定点Ｐｉを定義し、このγ線測定点Ｐｉから前記第３の記憶手段の前記現地γ線測定情報Ｆγｉの前記方向及び前記断層長Ｌの直線を描き、この直線を前記γ線現地測定結果直線ＬＧｊとして前記γ線測定点Ｐｉ、前記微細断層線ＤＬｊとは異なる色で表示する手段と、
前記γ線現地測定結果直線ＬＧｊに前記方向を割り付ける手段と、
前記γ線現地測定結果直線ＬＧｊがカーソルで指定された場合は、このγ線現地測定結果直線ＬＧｊに前記割り付けられている前記方向を読み込み、前記カーソルの移動に追従させて前記γ線現地測定結果直線ＬＧｊを前記読み込んだ前記方向で移動させる手段と、
を有することを特徴とする断層図作成装置。

【請求項2】

前記微細断層線ＤＬｊを描かせるための面の上に、前記微細断層線ＤＬｊを線で繋ぐための面を表示し、この面における前記カーソルの移動に追従して前記微細断層線ＤＬｊ、前記γ線測定点Ｐｉ、前記γ線現地測定結果直線ＬＧｊの色とは異なる色の線をハイブリッド断層線ＨＢＬｊとして前記微細断層線ＤＬｊを線で繋ぐための面に描画する手段と、
を有することを特徴とする請求項１記載の断層図作成装置。

【請求項3】

前記断層長Ｌ（Ｋｍ）は、
Ｌ＝０．３６ＦＷ＋０．６２
但し、ＦＷは断層幅（ｍ）
としていることを特徴とする請求項１記載の断層図作成装置。

【請求項4】

前記微細断層線ＤＬｊの幅は、１０ｍ以下であることを特徴とする請求項１記載の断層図作成装置。

【請求項5】

前記数値標高モデルは、
２．５ｍ以下のメッシュサイズであることを特徴とする請求項１記載の断層図作成装置。

【請求項6】

前記赤色立体化画像作成部は、
前記数値標高モデルの格子に標本点を設定し、この設定毎にこの標本点から設定された前記検索範囲の内に存在する対象格子までの地上開度、地下開度及び斜度を求めて、地上開度の値が大きさほどに明るい色を、地下開度の値が大きいさほどに暗い色を、斜度の値が大きいほどに赤が強調された色を割り付けて前記赤色立体化画像ＲＧｉを生成すること
を特徴とする請求項１乃至５のいずれかに記載の断層図作成装置。

【請求項7】

前記検索範囲は、
前記断層帯幅ＦＷに対して５倍〜１０倍程度の長さにされていることを特徴とする請求項１記載の断層図作成装置。

【請求項8】

前記検索範囲は、５０〜１００ｍ程度にされていることを特徴とする請求項７記載の断層図作成装置。

【請求項9】

微細断層線ＤＬｊを検出可能な間隔で航空機から地域にレーザパルスを発射して得た数値標高モデルが記憶された第１の記憶手段と、
赤色立体化画像ＲＧｉが記憶される第２の記憶手段と、
前記赤色立体化画像ＲＧｉに基づいて決定した前記地域におけるγ線測定点Ｐｉにおける複数測定ポイントのγ線強度と測定間隔とに基づいて求められた前記微細断層線ＤＬｊの断層長Ｌ及び前記γ線測定点Ｐｉでの前記微細断層線ＤＬｊの方向が現地γ線測定情報Ｆγｉとして記憶された第３の記憶手段とを用意し、
コンピュータが、
前記微細断層線ＤＬｊの断層帯幅ＦＷを検出可能な検索範囲を設定するための入力ボックスを表示部の画面に表示する断層用検出範囲設定部と、
前記第１の記憶手段の前記数値標高モデルと前記断層用検出範囲設定部によって設定された検索範囲に基づく前記赤色立体化画像ＲＧｉを生成して前記第２の記憶手段に記憶して前記画面に表示する赤色立体化画像作成ステップと、
前記画面の前記赤色立体化画像ＲＧｉ上に、前記微細断層線ＤＬｊを描かせるための面を表示し、この面におけるカーソルの移動軌跡を前記微細断層線ＤＬｊとし、これを所定の色で表示するステップと、
前記微細断層線ＤＬｊを描かせるための面の上に、前記γ線測定点Ｐｉを指定させるための面を表示し、この面に指定された点を前記γ線測定点Ｐｉとして前記微細断層線ＤＬｊとは異なる色で表示させるステップと、
前記微細断層線ＤＬｊを描かせるための面の上又は下にγ線現地測定結果直線ＬＧｊを表示するための面を表示し、この面に前記γ線測定点Ｐｉを指定させるための面に定義された前記γ線測定点Ｐｉを定義し、このγ線測定点Ｐｉから前記第３の記憶手段の前記現地γ線測定情報Ｆγｉの前記方向及び前記断層長Ｌの直線を描き、この直線を前記γ線現地測定結果直線ＬＧｊとして前記γ線測定点Ｐｉ、前記微細断層線ＤＬｊとは異なる色で表示するステップと、
前記γ線現地測定結果直線ＬＧｊに前記方向を割り付けるステップと、
前記γ線現地測定結果直線ＬＧｊがカーソルで指定された場合は、このγ線現地測定結果直線ＬＧｊに前記割り付けられている前記方向を読み込み、前記カーソルの移動に追従させて前記γ線現地測定結果直線ＬＧｊを前記読み込んだ前記方向で移動させるステップと、
を行うことを特徴とする断層図作成方法。

【請求項10】

前記コンピュータが、
前記微細断層線ＤＬｊを描かせるための面の上に、前記微細断層線ＤＬｊを線で繋ぐための面を表示し、この面における前記カーソルの移動に追従して前記微細断層線ＤＬｊ、前記γ線測定点Ｐｉ、前記γ線現地測定結果直線ＬＧｊの色とは異なる色の線をハイブリッド断層線ＨＢＬｊとして前記微細断層線ＤＬｊを線で繋ぐための面に描画するステップと
を行うことを特徴とする請求項９記載の断層図作成方法。

【請求項11】

前記断層長Ｌ（Ｋｍ）は、
Ｌ＝０．３６ＦＷ＋０．６２
但し、ＦＷは断層幅（ｍ）
としていることを特徴とする請求項９記載の断層図作成方法。

【請求項12】

前記微細断層線ＤＬｊの幅は、１０ｍ以下であることを特徴とする請求項９記載の断層図作成方法。

【請求項13】

前記数値標高モデルは、
２．５ｍ以下のメッシュサイズであることを特徴とする請求項９記載の断層図作成方法。

【請求項14】

前記赤色立体化画像作成ステップは、
前記数値標高モデルの格子に標本点を設定し、この設定毎にこの標本点から設定された前記検索範囲の内に存在する対象格子までの地上開度、地下開度及び斜度を求めて、地上開度の値が大きさほどに明るい色を、地下開度の値が大きいさほどに暗い色を、斜度の値が大きいほどに赤が強調された色を割り付けて前記赤色立体化画像ＲＧｉを生成することを特徴とする請求項９記載の断層図作成方法。

【請求項15】

前記検索範囲は、
前記断層帯幅ＦＷに対して５倍〜１０倍程度の長さにされていることを特徴とする請求項９記載の断層図作成方法。

【請求項16】

前記検索範囲は、５０〜１００ｍ程度にされていることを特徴とする請求項１５記載の断層図作成方法。

【請求項17】

微細断層線ＤＬｊを検出可能な間隔で航空機から地域にレーザパルスを発射して得た数値標高モデルが記憶された第１の記憶手段と、
赤色立体化画像ＲＧｉが記憶される第２の記憶手段と、
前記赤色立体化画像ＲＧｉに基づいて決定した前記地域におけるγ線測定点Ｐｉにおける複数測定ポイントのγ線強度と測定間隔とに基づいて求められた前記微細断層線ＤＬｊの断層長Ｌ及び前記γ線測定点Ｐｉでの前記微細断層線ＤＬｊの方向が現地γ線測定情報Ｆγｉとして記憶された第３の記憶手段とを用いて、
コンピュータを、
前記微細断層線ＤＬｊの断層帯幅ＦＷを検出可能な検索範囲を設定するための入力ボックスを表示部の画面に表示する断層用検出範囲設定手段、
前記第１の記憶手段の前記数値標高モデルと前記断層用検出範囲設定手段によって設定された検索範囲に基づく前記赤色立体化画像ＲＧｉを生成して前記第２の記憶手段に記憶して前記画面に表示する赤色立体化画像手段、
前記画面の前記赤色立体化画像ＲＧｉ上に、前記微細断層線ＤＬｊを描かせるための面を表示し、この面におけるカーソルの移動軌跡を前記微細断層線ＤＬｊとし、これを所定の色で表示する手段、
前記微細断層線ＤＬｊを描かせるための面の上に、前記γ線測定点Ｐｉを指定させるための面を表示し、この面に指定された点を前記γ線測定点Ｐｉとして前記微細断層線ＤＬｊとは異なる色で表示させる手段、
前記微細断層線ＤＬｊを描かせるための面の上又は下にγ線現地測定結果直線ＬＧｊを表示するための面を表示し、この面に前記γ線測定点Ｐｉを指定させるための面に定義された前記γ線測定点Ｐｉを定義し、このγ線測定点Ｐｉから前記第３の記憶手段の前記現地γ線測定情報Ｆγｉの前記方向及び前記断層長Ｌの直線を描き、この直線を前記γ線現地測定結果直線ＬＧｊとして前記γ線測定点Ｐｉ、前記微細断層線ＤＬｊとは異なる色で表示する手段、
前記γ線現地測定結果直線ＬＧｊに前記方向を割り付ける手段、
前記γ線現地測定結果直線ＬＧｊがカーソルで指定された場合は、このγ線現地測定結果直線ＬＧｊに前記割り付けられている前記方向を読み込み、前記カーソルの移動に追従させて前記γ線現地測定結果直線ＬＧｊを前記読み込んだ前記方向で移動させる手段、
としての機能を実行させるための断層図作成プログラム。

【請求項18】

前記コンピュータを、
前記微細断層線ＤＬｊを描かせるための面の上に、前記微細断層線ＤＬｊを線で繋ぐための面を表示し、この面における前記カーソルの移動に追従して前記微細断層線ＤＬｊ、前記γ線測定点Ｐｉ、前記γ線現地測定結果直線ＬＧｊの色とは異なる色の線をハイブリッド断層線ＨＢＬｊとして前記微細断層線ＤＬｊを線で繋ぐための面に描画する手段、
としての機能を実行させるための請求項１７記載の断層図作成プログラム。

【請求項19】

前記断層長Ｌ（Ｋｍ）は、
Ｌ＝０．３６ＦＷ＋０．６２
但し、ＦＷは断層幅（ｍ）
としていることを特徴とする請求項１７に記載の断層図作成プログラム。

【請求項20】

前記微細断層線ＤＬｊの幅は、１０ｍ以下であることを特徴とする請求項１７記載の断層図作成プログラム。

【請求項21】

前記数値標高モデルは、
２．５ｍ以下のメッシュサイズであることを特徴とする請求項１７記載の断層図作成プログラム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、航空機で取得した地域の画像上に、現地の計測結果に基づいた断層線を精度よく容易に繋げさせることができる断層図作成装置に関する。

【背景技術】

【0002】

内陸地震がしばしば発生し、断層の位置や活動性の把握が重要になっている。内陸地震の原因として活断層があるが、従来知られていない活断層が近年新たに発見されることがしばしばある。

【0003】

このような断層の位置把握は航空機でステレオ写真を撮影し、このステレオ写真を立体視化して断層を検出する方式を主としている。

【0004】

しかし、植生に覆われる斜面が広く分布する地域では１０ｍ以下の微細な断層変位地形（以下微細断層という）を認定することは困難である。

【0005】

これに対して、航空機からレーザを発射して反射光を取得して地表面の形状を表現するＤＥＭ（数値標高モデル）を生成し、このＤＥＭを用いて地表面の凹凸を強調する地図作成装置が特許文献１に開示されている。

【0006】

特許文献１の地形図作成装置は、数値標高モデルの標高データに基づいて地表面の微地形を表現すべく微細な凹凸を強調する。つまり、高低差を増幅することで凹凸を強調している。そして、この凹凸の彩度を標高に比例させて決定している。また、標高によって色相を決定して、これらを合成したカラー画像を得る。

【先行技術文献】

【特許文献】

【0007】

【特許文献1】特開２００８−２４２２９８号公報

【発明の概要】

【発明が解決しようとする課題】

【0008】

しかしながら、活断層と思われる微細断層というのは、幅数メートル程度が一般的である。ところが、特許文献１の地形図作成装置は、標高の局所的な平均値を維持しつつ微細な凹凸の標高差を増幅して強調するものであるから、幅数メートル程度の微細断層の検出に不向きである。

【0009】

また、断層とする箇所を写真から判読して地図上に描いたとしても、断層は数キロ以上にも渡るのが一般的である。また、断層は複雑に入り組んでいる。

【0010】

このため、ある断層がどの断層と繋がって、どの方向にどの程度の長さに渡っているかは容易には判読できない。

【0011】

さらに、断層長というのは、断層幅と地盤の活動性に依存すると言われている。このため、現地で測定した断層幅、活動性、微細断層の方向を考慮して断層線同士を繋げるのが望ましい。

【0012】

本発明は以上の課題を解決するためになされたもので、断層線を描く対象地域の画像上に、現地の計測結果を反映した断層長で、かつ断層の方向で断層線同士を容易に精度よく繋げさせることができる断層線図を得ることができる断層図作成装置を得ることを目的とする。

【課題を解決するための手段】

【0013】

本発明に係る断層図作成装置は、微細断層線ＤＬｊを検出可能な間隔で航空機から地域にレーザパルスを発射して得た数値標高モデルが記憶された第１の記憶手段と、
赤色立体化画像ＲＧｉが記憶される第２の記憶手段と、
前記赤色立体化画像ＲＧｉに基づいて決定した前記地域におけるγ線測定点Ｐｉにおける複数測定ポイントのγ線強度と測定間隔とに基づいて求められた前記微細断層線ＤＬｊの断層長Ｌ及び前記γ線測定点Ｐｉでの前記微細断層線ＤＬｊの方向が現地γ線測定情報Ｆγｉとして記憶された第３の記憶手段と、

【0014】

前記微細断層線ＤＬｊの断層帯幅ＦＷを検出可能な検索範囲を設定するための入力ボックスを表示部の画面に表示する断層用検出範囲設定部と、
前記第１の記憶手段の前記数値標高モデルと前記断層用検出範囲設定部によって設定された検索範囲に基づく前記赤色立体化画像ＲＧｉを生成して前記第２の記憶手段に記憶して前記画面に表示する赤色立体化画像作成部と、
前記画面の前記赤色立体化画像ＲＧｉ上に、前記微細断層線ＤＬｊを描かせるための面を表示し、この面におけるカーソルの移動軌跡を前記微細断層線ＤＬｊとし、これを所定の色で表示する手段と、
前記微細断層線ＤＬｊを描かせるための面の上に、前記γ線測定点Ｐｉを指定させるための面を表示し、この面に指定された点を前記γ線測定点Ｐｉとして前記微細断層線ＤＬｊとは異なる色で表示させる手段と、
前記微細断層線ＤＬｊを描かせるための面の上又は下にγ線現地測定結果直線ＬＧｊを表示するための面を表示し、この面に前記γ線測定点Ｐｉを指定させるための面に定義された前記γ線測定点Ｐｉを定義し、このγ線測定点Ｐｉから前記第３の記憶手段の前記現地γ線測定情報Ｆγｉの前記方向及び前記断層長Ｌの直線を描き、この直線を前記γ線現地測定結果直線ＬＧｊとして前記γ線測定点Ｐｉ、前記微細断層線ＤＬｊとは異なる色で表示する手段と、
前記γ線現地測定結果直線ＬＧｊに前記方向を割り付ける手段と、
前記γ線現地測定結果直線ＬＧｊがカーソルで指定された場合は、このγ線現地測定結果直線ＬＧｊに前記割り付けられている前記方向を読み込み、前記カーソルの移動に追従させて前記γ線現地測定結果直線ＬＧｊを前記読み込んだ前記方向で移動させる手段とを備えたことを要旨とする。

【発明の効果】

【0015】

以上のように、本発明によれば、画面に表示された赤色立体化画像ＲＧｉから微細断層線ＬＧｊを判読し、この微細断層線ＬＧｊが存在する現地のポイントにおいて複数ポイントでγ線測定して得た現地γ線測定情報Ｆγｉを記憶しておく。
そして、画面の赤色立体画像ＲＧｉｉ上でγ線測定点Ｐｉが指定された場合は、このγ線測定点Ｐｉを基準として現地γ線測定情報Ｆγｉに含まれている断層長Ｌを現地γ線測定情報Ｆγｉに含まれている現地で得た方向で表示し、このγ線現地測定結果直線ＬＧｊを前述の方向に対して逆方向に移動させて表示させる。

【0016】

このため、赤色立体化画像ＲＧｉ上に表示されている小まぎれな各々の微細断層線ＤＬｊに対して、γ線現地測定結果直線ＬＧｊがほぼ平行になるように表示されるので、一連の断層線が一つの線で繋がる断層であるかを容易に判断できる。

【図面の簡単な説明】

【0017】

【図1】実施の形態１の断層図作成装置の概略構成図である。

【図2】ＤＥＭデータの説明図である。

【図3】ハイブリッド断層線を得るまでの作業手順を説明するフローチャートである。

【図4】赤色立体化画像ＲＧｉの説明図である。

【図5】断層の説明図である。

【図6】赤色立体化画像ＲＧｉで見られる微細断層線ＤＬｊの説明図である。

【図7】断層の説明図である。

【図8】等高線図と赤色立体化画像ＲＧｉとの比較説明図である。

【図9】等高線図と赤色立体化画像ＲＧｉとの比較説明図である。

【図10】図３の第２工程（２）の説明図である。

【図11】図３の第３工程（３）の説明図である。

【図12】現地の測定状況の説明図である。

【図13】現地の測定状況のγ線測定の測定図である。

【図14】現地の測定状況のγ線測定の測定方法の説明図である。

【図15】現地の測定状況のγ線測定の測定方法の説明図である。

【図16】図３の第４工程（４）の説明図である。

【図17】図３の第５工程の処理を説明するフローチャートである。

【図18】図３の第５工程の処理を説明するフローチャートである。

【図19】γ線測定点Ｐｉの表示の説明図である。

【図20】γ線測定点・断層線描画部２４のダイアログボックスＢＰｉの説明図である。

【図21】γ線測定点Ｐｉの関連付情報の説明図である。

【図22】γ線現地測定結果直線ＬＧｊの説明図である。

【図23】γ線現地測定結果直線ＬＧｊの移動の説明図である。

【図24】ハイブリッド断層線ＨＢＬｊの作成の説明図である。

【図25】γ線現地測定結果直線ＬＧｊの削除の説明図である。

【図26】実施の形態２の断層図作成装置の概略構成図である。

【図27】実施の形態２の断層図作成装置の表示画面説明図である。

【図28】地上開度及び地下開度の原理説明図である。

【図29】地上開度及び地下開度の主要パターン説明図である。

【図30】地上開度及び地下開度の立体的説明図である。

【図31】地上開度及び地下開度の標本点及び距離の説明図である。

【図32】実施の形態３の断層図作成装置の概略構成図である。

【図33】レイヤ順位設定部の階層の説明図である。

【図34】メモリ８０のデータ構成の説明図である。

【図35】γ線現地測定結果直線ＬＧｊを得る過程の説明図である。

【図36】γ線現地測定結果直線ＬＧｊの移動の説明図である。

【図37】ハイブリッド断層線ＨＢＬｊの属性情報の説明図である。

【図38】ハイブリッド断層線ＨＢＬｊの描画の説明図である。

【図39】第１層をγ線現地測定結果直線ワーキング用メモリ７６にした場合のレイヤ順位設定部の説明図である。

【図40】第１層をγ線現地測定結果直線ワーキング用メモリ７６にした場合のγ線現地測定結果直線ＬＧｊの説明図である。

【発明を実施するための形態】

【0018】

以下に示す実施の形態は、発明の技術的思想を具体化するための装置や方法を例示するものであって、発明の技術的思想は、構造、配置は下記のものに特定するものではない。
発明の技術的思想は、特許請求の範囲に記載された技術的範囲内において、種々の変更を加えることができる。図面は模式的なものであり、装置やシステムの構成等は現実のものとは異なることに留意すべきである。

【0019】

一般に、活断層と思われる微細断層というのは、数キロにも渡る断層長であり、この断層長は断層幅に依存しているといわれている。また、断層長というのは、地盤の活動性が影響するといわれている。この活動性の検出にγ線測定装置を用いている。

【0020】

つまり、実施の形態１においては、航空レーザ計測データから得た赤色立体化画像ＲＧｉを画面に表示させて微細断層を判読させ、判読した微細断層を描画させて、現地で測定したγ線測定結果を用いた断層長、及び方向に基づいて微細断層同士を容易に精度よく繋げられるようにしている。

【0021】

＜実施の形態１＞
図１は実施の形態１の断層図作成装置の概略構成図である。実施の形態１の断層図作成装置は、図１に示すように、パソコン本体部１００に以下のプログラム構成を備えて表示部２００の画面にハイブリッド断層図ＨＨＧｉを表示する。

【0022】

パソコン本体部１００は、赤色立体化画像作成部１０と、断層用検索範囲設定部１２と、画像化部１４と、範囲定義部１６と、既存・赤色判読描画部１８と、γ線測定点・断層線描画部２４と、線調整部２６等のプログラム構成を備えている。

【0023】

さらに、パソコン本体部１００は、記憶手段として、図１に示すように、数値標高モデル（ＤＥＭ：Ｄｉｇｉｔａｌ Ｅｌｅｖａｔｉｏｎ Ｍｏｄｅｌ）を記憶したＤＥＭ用データベース３０と、赤色立体化画像ＲＧｉが記憶されるメモリ３２と、画像化部１４のメモリ３９と、赤色立体化画像ＲＧｉに基づいて判読した断層帯幅が１０ｍ以下の微細断層線（以下微細断層線ＤＬｊという）を描画するための微細断層線定義用メモリ３４と、既存の断層線（以下既存断層線ＫＬｊという）が描画される既存断層線用メモリ３６と、ハイブリッド断層線ＨＢＬｊが描画されるメモリ３８と、現地γ線測定情報Ｆγｉが記憶された現地測定結果情報用メモリ４０等を備えている。本実施の形態では、微細断層線ＤＬｊは、断層線の幅が約１０ｍ以下をいう。

【0024】

前述のＤＥＭ用データベース３０に記憶された数値標高モデル（Ｄｉｇｉｔａｌ Ｅｌｅｖａｔｉｏｎ Ｍｏｄｅｌ）とは、計測地域全体に所望の格子間隔ｄ（例えば、０．２ｍ、０．５ｍ、１ｍ又は１，５ｍ、２．５ｍ等）の格子構造を設定し、レーザ反射パルスのうち主に最後に返ってきたパルス（ラストパルス）によって計測された標高データから、地表面以外の建物や樹木などを取り除くフィルタリングを行い、標高値内挿補間法によって得た地盤の格子状の標高データである。

【0025】

図２に具体例を示す。図２に示すように、ＤＥＭデータＤｉは、格子番号ｉ（ｉ＝１，２，・・・，ｎ）を付与した各格子の中心点のＸ座標（経度Ｘｉ）、Ｙ座標（緯度Ｙｉ）、Ｚ座標（地盤標高値Ｚｇｉ）によって構成される。

【0026】

標高値内挿補間法の例としては、同じ標高値を結んだ等高線図を作成し、この等高線図に対して不整三角形網（ＴＩＮ）を作成して地盤を復元し、ＴＩＮと各格子の交点の高さを求める方法がある。

【0027】

赤色立体化画像作成部１０は、ＤＥＭ用データベース３０の数値標高モデルに標本格子点を順次設定し、この標本格子点毎に断層用検索範囲設定部１２によって設定された断層を検出するため検索範囲を設定し、この断層検索範囲内において、ＤＥＭ（Digital Elavation Model）データＤｉに基づいて斜度と、地上開度と、地下開度との３つのパラメータを求め、その平面分布をグレイスケール画像として保存する。
なお、断層用検索範囲設定部１２は、前述の微細断層線ＤＬｊの断層帯幅ＦＷを検出可能な検索範囲を設定するための入力ボックス（図示せず）を表示する。

【0028】

そして、地上開度と地下開度の差分画像にはグレイを、傾斜を赤のチャンネルにいれて、擬似カラー画像を作成することにより、尾根や山頂部分を白っぽく、また谷や窪地を黒っぽく表現し、傾斜が急な部分ほど赤く表現した赤色立体化地図ＲＧｉを生成して、メモリ３２に記憶する。この赤色立体化画像作成部１０について図を用いて後述する。

【0029】

画像化部１４は、出力された各種画像を合成（重ねて）して表示部２００に表示する。このとき、画像に付加されている色情報で画面に表示する。

【0030】

範囲定義部１６は、赤色立体化画像作成部１０を介してメモリ３２に記憶された赤色立体化画像ＲＧｉの平面範囲（Ｘ、Ｙ）を読み、これらの範囲を微細断層線定義用メモリ３４、既存断層線用メモリ３６、メモリ３８に定義する。この微細断層線定義用メモリ３４に定義された平面を赤色判読結果用平面と称し、既存断層線用メモリ３６に定義された平面を既存用平面と称し、メモリ３８に定義された平面図をハイブリッド用平面と称する。

【0031】

すなわち、微細断層線定義用メモリ３４の赤色判読結果用平面は、赤色立体化画像ＲＧｉの原点座標を有し、かつ赤色立体化画像ＲＧｉの解像度のメッシュサイズにされて定義されている。

【0032】

また、既存断層線用メモリ３６の既存用平面は赤色立体化画像ＲＧｉの原点座標を有し、かつ赤色立体化画像ＲＧｉの解像度のメッシュサイズ又はこの解像度より細かいメッシュサイズにされて定義されている。

【0033】

また、メモリ３８のハイブリッド用平面は、赤色立体化画像ＲＧｉの原点座標を有し、かつ赤色立体化画像ＲＧｉの解像度のメッシュサイズにされて定義されている。

【0034】

既存・赤色判読描画部１８は、画面の赤色立体化画像ＲＧｉ上に微細断層線ＤＬｊを描く指示がある場合は、画面上のカーソルの移動に伴って、このカーソルの移動に伴う座標を微細断層線定義用メモリ３４の赤色判読結果用平面の座標に変換して定義し、この平面に定義された線を微細断層線ＤＬｊとして画像化部１４によって表示させる。

【0035】

また、既存・赤色判読描画部１８は、画面の赤色立体化画像ＲＧｉ上に既存断層線ＫＬｊを描く指示がある場合は、画面上のカーソルの移動に伴って、この座標を既存断層線用メモリ３６の既存用平面面の座標に変換して定義し、この平面に定義された線を既存断層線ＫＬｊとして画像化部１４によって表示させる。

【0036】

また、現地測定結果情報用メモリ４０には、画面の赤色立体化画像ＲＧｉに基づいて決定した現地におけるγ線測定点Ｐｉにおける複数測定ポイントのγ線強度と測定間隔に基づいて求められたγ線現地測定結果直線ＬＧｊの断層長Ｌ及び地域のγ線測定点Ｐｉでのγ線現地測定結果直線ＬＧｊの方向Ｌθｉが記憶されている。これらを総称して現地γ線測定情報Ｆγｉと称する。この現地γ線測定情報Ｆγｉには、γ線現地測定結果直線ＬＧｊの線番号ＬＧＢｉを生成してγ線測定点Ｐｉに関連付けている。

【0037】

γ線測定点・断層線描画部２４は、表示部２００の画面上に表示された赤色立体化画像ＲＧｉにカーソルで指示された点をγ線測定点Ｐｉとして読み込み、このγ線測定点Ｐｉをメモリ３８のハイブリッド用平面の座標に変換して定義する。そして、現地測定結果情報用メモリ４０に記憶されている現地γ線測定情報Ｆγｉに含まれている方向Ｌθｉで断層長Ｌの長さを有する線をγ線現地測定結果直線ＬＧｊとして表示させる。そして、このγ線現地測定結果直線ＬＧｊがカーソルで指定され、このカーソルの移動に伴ってこのカーソルの座標を微細断層線定義用メモリ３４の赤色判読結果用平面の座標に変換して表示させる。この移動は方向Ｌθｉに対して逆方向に移動させている。

【0038】

線調整部２６は、画面上においてγ線現地測定結果直線ＬＧｊと既存断層線ＫＬｊとが比較されて、ハイブリッド断層線ＨＢＬｊの作成指示の入力がある場合は、カーソルの座標を赤色判読結果用平面の座標に定義する。
そして、このカーソルの移動に伴って、この移動座標を微細断層線定義用メモリ３４の赤色判読結果用平面の座標に変換してハイブリッド断層線ＨＢＬｊ（ピンク）として表示させる。

【0039】

つまり、画面に表示している微細断層線ＤＬｊ同士がなぞられながら線で繋げられた場合は、これらの線及び該繋げられた微細断層線ＤＬｊをハイブリッド断層線ＨＢＬｊとして異なる色（ピンク）で表示させる。

【0040】

（詳細説明）
図３はハイブリッド断層線を得るまでの作業手順を説明するフローチャートである。

【0041】

図３に示すように、航空機によって対象地域にレーザを発射して２．５ｍ以下（例えば1m）の メッシュのＤＥＭをＤＥＭ用データベース３０に得る（Ｓ１００）。

【0042】

そして、赤色立体化画像ＲＧｉを赤色立体化画像作成部１０によって生成する（Ｓ２００）。図４に赤色立体化画像作成部１０によって生成した赤色立体化画像ＲＧｉを示す。この赤色立体化画像作成部１０によって生成を第１工程（１）とも称する。

【0043】

この赤色立体化画像ＲＧｉの生成において検索範囲を小さくとる（５０〜１００ｍ程度）。これは、赤色立体化画像ＲＧｉは単画像で地形が立体的に見えるように加工した画像であり、画像解像度は航空レーザデータのメッシュサイズと同じ1m（例えば）となるため、例えば活断層判読のキーとなる微地形も十分に判読できる。

【0044】

ここで、赤色立体化画像ＲＧｉによって微細幅の断層線（微細断層線ＤＬｊ（太い黒点線））を検出できる理由を説明する。図５の（ａ）〜（ｄ）に示すように、断層は線状に延びた段差地形や、尾根・谷等の地形がある線を境にして一斉にずれている場合には、断層であることが疑われる。

【0045】

ただし、段差地形が線状に延びていても、その方向が斜面の最大傾斜方向や河川流下方向に平行である場合には、断層ではなく侵食地形である可能性が高い。また、地盤を形成する岩石が、ある線を境にして異なっており、両者の耐侵食性が異なる場合にも、段差地形が形成されることがある。断層と疑われる地形があった場合、周辺を含めた地形形成過程を何パターンか想像して比較し、その地形要因が断層であると考えた方がより容易に説明できる場合に、断層と推定する。

【0046】

断層と疑われる地形を赤色立体化画像ＲＧｉで見た場合には、赤色立体化画像ＲＧｉは、例えば１ｍメッシュのＤＥＭを用いて検索範囲５０〜１００ｍ程度で演算するので、図６に示すように周辺と特徴的に異なる色調が、線形に並ぶ事になる。これは傾斜などの地形量が、断層地形の部分だけ周辺と異なるためである。つまり、図７に示すように、図６のＡ−Ｂ間の断面は地盤がずれている部分を検出できる。

【0047】

ここで、等高線図と比較して説明する。

【0048】

等高線図では、等高線どうしの間に何の情報もない（図８（ｂ）、図９（ｂ）参照）。
このため、段差などの地形の連続性を追うことが難しいので結果として地形スケール（段差の高さや幅）が微小であればあるほど難しくなる。

【0049】

また、断層による鉛直変位（小さい場合には１ｍ以下）が等高線間隔（一般に縮尺１/２５０００ で１０ｍ、１/２５００ で２ｍ程度）よりも小さい場合（図８（ｂ）、図９（ｂ）参照）には、等高線では表現不可能である。

【0050】

これに対して、１ｍメッシュ程度（０.２ｍ、０.５ｍでもかまわない）のＤＥＭを用いて５０〜１００ｍ範囲の検索範囲で赤色立体化画像ＲＧｉを作成しているので図８（ａ）、図９（ａ）に示すように、断層線（微細断層線ＤＬｊ）が明らかに判読できる。

【0051】

また、既存資料（地図、地質図等）を読み込み（Ｓ３００）、かつ赤色立体化画像ＲＧｉにおいて断層を判読する（Ｓ４００）。

【0052】

そして、既存資料から断層（既存断層線ＫＬｊ）と分かっているところは、赤色立体化画像ＲＧｉ上に移写（図１０（ａ）参照）する一方、赤色立体化画像ＲＧｉにおいて断層（微細断層線ＤＬｊ）と判読できる箇所は線（図１０（ｂ）参照）で描く（Ｓ５００）。これを第２工程（２）とも称する。

【0053】

つまり、赤色立体化画像ＲＧｉにより、従来手法（空中写真判読）では判読できなかった微細な断層を判読して、その判読結果を描いている。

【0054】

しかし、断層「と考えられる地形」を判読することはできるが、断層でないリニアメント（例えば地質境界に沿って発達した地形ギャップ）を含む可能性がある。

【0055】

例えば、 堆積土砂で覆われている場合や、逆に渓流の侵食や人工改変によって断層が後から分断されている場合、断層の連続性が判定できない。

【0056】

このため、現地でγ線測定によってこれらを調査する（Ｓ６００）。これを第３工程（３）とも称する。

【0057】

現地でγ線測定によってこれらを調査する前に図１１に示すように、工程（２）で判読された断層に対して、例えば、厚い堆積土砂に覆われていない尾根部や鞍部、岩が地表近くに存在する地点で、かつ車道から近い等アプローチの良い箇所を選定する。但し、図１１において等高線を重ねて示している。

【0058】

これによって、例えば 工程（２）の赤色立体地図の判読によって断層の存在・位置・方向などを推定できるだけでなく、地形・地質等の条件面からもγ線測定（現地調査）の適地を選定することができる。また、現地調査の効率化に繋がる。さらに、より多くのポイントで現地調査を行うことができる。

【0059】

また、 地表に土砂が厚く堆積していない場所など適切な測定場所を選定することにより、１０ｃｍ オーダーの精度で断層破砕帯の境界位置を決めることができる。

【0060】

このようにして決めた測定場所における測定状況を以下に説明する。図１２（ａ）は現地の測定場所で用いるγ線測定装置の写真図である。図１２（ｂ）は測定状況の写真図である。なお、γ線探査箇所数は計６ カ所である。γ線測定装置は、シンチレーションサーベイメータ（全計数法）を用いている。

【0061】

測点間隔は例えば3m である。そして、γ線強度異常値が出現した地点では１０ｃｍ 間隔で測定し、γ線強度異常値範囲(破砕帯)を詳細に求める。

【0062】

図１２（ｂ）に示すようにγ線測定装置は破砕帯の方向に平行で、測線方向に直角に地面に設置される。

【0063】

つまり、推定される断層破砕帯の方向と直交方向に測線を設ける。次に、表層土砂をはぎ取る。次に、測線上を移動しながら地上のγ線強度を計測し、図１３に示すように記録していく。

【0064】

図１３において、測線上のある区間で、その両側に比較して明らかにγ線強度が急に変化していれば、その区間が断層破砕帯であると推定できる。

【0065】

そして、並行する測線を設定し、同様の計測を行って、断層破砕帯の位置を特定する。２つの測線のγ線計測から破砕帯の幅・方向が特定できる。

【0066】

図１４及び図１５を用いて計測の方法を具体的に説明する。図１４に示すように、ＳＴＥＰ「０」において、第３工程（３）で決定した地盤の調査地点に測線Ａ（ロープ）を設定する。

【0067】

ＳＴＥＰ「１−１」において、粗い間隔（例えば５ｍ）でγ線測定装置を配置してそれが、異常値を示すか正常値を示すかをプロットする。

【0068】

ＳＴＥＰ「１−２」において、間隔を狭める（例えば０.５〜２ｍ）。ＳＴＥＰ「１−３」において、異常値の左右境界を例えば０.１ｍまで狭める。つまり、図１４における（ａ１）と（ａ２）とが異常区間の境界である。

【0069】

そして、図１５の「ＳＴＥＰ２−１」において、測線Ａ（ロープ）と平行に測線Ｂ（ロープ）を地面に設置する。

【0070】

そして、「ＳＴＥＰ２−２」において、測線Ｂでのγ線強度を測定する（図１４と同様な方法）。そして、「ＳＴＥＰ２−３」において、測線Ｂにおけるγ線異常境界位置ｂ１、ｂ２を決定する。次に、「ＳＴＥＰ２−４」において、測線Ａの異常境界位置ａ１、ａ２及び測線Ｂの異常境界位置ｂ１、ｂ２を結ぶ。そして、「ＳＴＥＰ２−４」において測線Ａと測線Ｂとを直角に結ぶ線の長さを断層帯幅ＦＷと決定する。

【0071】

そして、このようにして現地で決定した断層帯幅ＦＷを用いてγ線測定結果の整理・解析を行う（Ｓ７００）。これを本実施の形態では第４工程（４）と称する。

【0072】

第４工程（４）について図１６を用いて説明する。第４工程（４）は第３工程（３）において現地で取得した断層帯幅ＦＷをパソコンに入力して断層長Ｌを求めさせる。

【0073】

例えば、現地での測線Ａ、Ｂでのγ線強度測定結果を図１６（ａ）に示すようにグラフに描いた場合に、断層帯幅ＦＷが８．８ｍの場合は、図１６（ｂ）に示すように断層長Ｌは、Ｌ＝０．３６ＦＷ＋０．６２として求められる。

【0074】

本実施の形態では断層帯幅ＦＷが８．８ｍであるので断層長Ｌは３．８Ｋｍと求められる。なお、断層長Ｌを求めるための係数「０．３６」、「０．６２」は一例である。

【0075】

そして、これらの計測データを図１に示すＬｒＧ法を用いた断層図作成装置に入力して図３に示すように、断層線（γ線現地測定結果直線ＬＧｊ）を自動作成する（Ｓ８００）。これを本実施の形態においては第６工程（６）と称する。

【0076】

この断層線（γ線現地測定結果直線ＬＧｊ）の自動作成処理については図を用いて後述する。

【0077】

また、断層線（γ線現地測定結果直線ＬＧｊ）の自動作成に当たっては図３に示すよう赤色立体化画像ＲＧｉにおいて断層と考えられる箇所には断層線（微細断層線ＤＬｊ）を描いておく。

【0078】

そして、γ線現地測定結果直線ＬＧｊに基づいて第６工程で得られた微細断層線ＤＬｊを修正して（Ｓ９００）、最終的にハイブリッド断層線ＨＢＬｊを得る（Ｓ１０００）。

【0079】

以下に第５工程で行うＬｒＧ法を用いた断層図作成装置における断層線の自動作成処理（γ線測定点・断層線描画部２４）について図１７のフローチャートを用いて説明する。

【0080】

図１７に示すように初めに初期設定を行う（Ｓ８０１）。

【0081】

この初期設定は、断層用検索範囲設定部１２によって検索範囲（例えば５０ｍ）が赤色立体化画像作成部１０に設定されて、赤色立体化画像作成部１０が図４に示す赤色立体化画像ＲＧｉを表示部２００に表示する。

【0082】

前述の検索範囲を５０ｍ程度にする理由は、通常の赤色立体地図では、微地形と大地形を同時に把握可能なように表現するため、探索範囲を数百メートルとする場合が多い。

【0083】

しかし目的を断層地形判読に特化する場合には、探索範囲を５０ｍ程度とすることによって、より微地形のみを認識しやすい画像とすることができ、かつ計算時間を短縮することができるためである。

【0084】

また、範囲定義部１６によって赤色判読結果用平面が定義され、既存断層線用メモリ３６に既存用平面が定義され、メモリ３８にハイブリッド用平面が定義されているとする。

【0085】

さらに、現地調査によって断層長Ｌが得られ、かつ断層長Ｌの方向Ｌθｉが得られているとする。

【0086】

また、第２工程によって図１０（ａ）に示すように既存資料による既存断層線ＫＬｊ（黒）が赤色立体化画像ＲＧｉに描かれて、かつ図１０（ｂ）に示すように赤色立体地図ＲＧｉに微細断層線ＤＬｊ（黒点線）が描かれているとする。

【0087】

このような状態において、図１７に示すように、γ線測定点・断層線描画部２４が表示部２００に表示している赤色立体地図ＲＧｉ上にγ線測定点Ｐｉが設定されたかどうかを判定する（Ｓ８０３）。

【0088】

ステップＳ８０３において、赤色立体地図ＲＧｉ上にγ線測定点Ｐｉが設定された判定した場合は、γ線測定点・断層線描画部２４は、このγ線測定点Ｐｉをメモリ３８のハイブリッド用平面に定義すると共に、図１９に示すようにγ線測定点Ｐｉを色表示（例えばイエロー）する（Ｓ８０５）。

【0089】

そして、図２０に示すようにγ線測定点・断層線描画部２４は距離と角度を要求するダイアログボックスＢＰｉを表示する（Ｓ８０７）。ダイアログボックスＢＰｉには、例えば計算式から求められた断層長Ｌ＝３８００ｍを入力する。また、例えば、断層の方向（走行Ｎ６４°Ｗ）を、東を０°として反時計回りに図った角度に換算した「１５４°」を入力する。

【0090】

次に、γ線測定点・断層線描画部２４はこのダイアログボックスＢＰｉに入力された距離（３８００ｍ）と角度（１５４°）を読み込んでラインツールに渡す（Ｓ８０９）。つまり、メモリ３８のハイブリッド用平面にγ線測定点Ｐｉが定義されると、図２１に示すように、γ線測定点・断層線描画部２４がγ線測定点ＰｉにラインツールとダイアログボックスＢＰｉに入力された距離（３８００ｍ）と角度（１５４°）とを関連付ける。

【0091】

前述のラインツールは、距離（３８００ｍ）と角度（１５４°）とを読み込み、γ線測定点Ｐｉを原点（基準）にして角度（１５４°）で距離（３８００ｍ）のγ線現地測定結果直線ＬＧｊをメモリ３８のハイブリッド用平面に描くと共に（Ｓ８１１）、画像化部１４によって図２２に示すように、γ線現地測定結果直線ＬＧｊ（黒点線）を赤色立体化画像ＲＧｉ上に重ね表示する（Ｓ８１３）。

【0092】

図２２に示すように、γ線測定点Ｐｉを原点としているので、西方向のみにγ線現地測定結果直線ＬＧｊ（イエロー）が引かれている。

【0093】

そして、図１８に示すように、γ線測定点・断層線描画部２４はγ線現地測定結果直線ＬＧｊ（黒点線）の調整かどうかを判断する（Ｓ８１５）。

【0094】

具体的にはγ線現地測定結果直線ＬＧｊがカーソルで指定(クリック)された時点でγ線現地測定結果直線ＬＧｊの調整と判断している。

【0095】

そして、ステップＳ８１５において、γ線現地測定結果直線ＬＧｊの調整と判断した場合は、γ線現地測定結果直線ＬＧｊの調整を受け付けて、図２３に示すように、新たに指定された点までγ線現地測定結果直線ＬＧｊを移動させる（Ｓ８１７）。

【0096】

具体的にはγ線測定点Ｐｉを通って、角度（１５４°）を維持したままクリックがはなされるまでγ線現地測定結果直線ＬＧｊを移動させる。つまり、メモリ３８のハイブリッド用平面のγ線測定点Ｐｉを通って、角度（１５４°）を維持したままクリックがはなされるまでγ線現地測定結果直線ＬＧｊを移動させ、これを画像化部１４によって表示させる。
つまり、断層面が傾斜している場合、断層面と地表面の交わる線＝断層線は完全な直線形状とはならない。このため、線分の近傍を通り連続性が良い「既存資料による断層」や「判読による断層と考えられる地形」を、一連の断層であると判断し、それらをトレースした線を描いている。

【0097】

これによって、図２３に示すように図２２におけるγ線現地測定結果直線ＬＧｊがγ線測定点Ｐｉを通って角度（１５４°）で平行に移動させられているから、図２２においては赤色立体化画像ＲＧｉにおいて３８００ｍのγ線現地測定結果直線ＬＧｊの全体が見えていなかったのが、図２３に示すように、γ線測定点Ｐｉを通るであろうと推定される全ての微細断層線ＤＬｊに対応付けることになるので、各々の微細断層線ＤＬｊがどのように繋がっているのかを容易に把握できることになる。

【0098】

そして、図２４に示すようにγ線現地測定結果直線ＬＧｊと微細断層線ＤＬｊ（点線で示している：青）との関係を目視で確認しながら各々の微細断層線ＤＬｊが繋がる断層線であるかどうかを判断し、繋がると判断した場合は、オペレータはハイブリッド断層線ＨＢＬｊ（ピンク）の作成指示を行って図２４に示すように微細断層線ＤＬｊを接続（太線：色はピンク）する（Ｓ８１９）。つまり、メモリ３８のハイブリッド用平面には太線が定義される。

【0099】

そして、図２５に示すように、γ線現地測定結果直線ＬＧｊを消去して（Ｓ８２１）、終了かどうかを判断する（Ｓ８２３）。終了と判断した場合は、画像化部１４によって画面の画像をハイブリッド断層図ＨＨＧｉとしてメモリ３９に記憶させる（Ｓ８２５）。

【0100】

従って、断層の連続性は活動度評価に大きな影響を及ぼすので、赤色立体地図からの判読だけでは不明であった断層の連続性を、γ線測定からのγ線現地測定結果直線ＬＧｊによって補うことができている。

【0101】

さらに、これまで見えなかった断層破砕帯の分布と、それらの特性（幅、長さ等）を広域に把握することができる。こうして得られた情報をもとに、活断層の活動度評価（断層断層帯幅・長さや地質年代から評価）することも可能である。

【0102】

従って、鉄道・道路等インフラストラクチャに対する地震や土砂災害のリスク評価、深層崩壊等の土砂災害に関する発生・影響予測などに活用することができる。

【0103】

＜実施の形態２＞
図２６は実施の形態２の断層図作成装置の概略構成図である。図２６においては図１と同一符号のものについては説明を省略する。

【0104】

図２６に示すように実施の形態２の断層図作成装置は、赤色立体化画像ＲＧｉと同じ範囲のオルソフォト画像が記憶されたデータベース４１と、赤色立体化画像ＲＧｉと同じ範囲の等高線図が記憶されたデータベースと、切換部４４とを備えている。

【0105】

切換部４４は、オペレータによってオルソフォト又は等高線図が指定されると、データベース４１のオルソフォト画像又はデータベース４２の等高線図を読み出して画像化部１４を介して表示させる。

【0106】

例えば、図２７（ａ）に示すように実施の形態１によってハイブリッド断層図を表示している場合に、オペレータがオルソフォトを入力した場合は、切換部４４はデータベース４１のオルソフォト画像を画像化部１４に出力する。画像化部１４は切換部４４からオルソフォト画像が出力された場合は、メモリ３２からの赤色立体化画像ＲＧｉの読み込みを停止して、赤色立体化画像ＲＧｉに代えて切換部４４からオルソフォト画像を図２７（ｂ）に示すように表示する。

【0107】

また、図２７（ａ）に示すように実施の形態１によってハイブリッド断層図を表示している場合に、オペレータが等高線図を入力した場合は、切換部４４はデータベース４２の等高線図を画像化部１４に出力する。画像化部１４は切換部４４から等高線図が出力された場合は、メモリ３２からの赤色立体化画像ＲＧｉの読み込みを停止して、赤色立体化画像ＲＧｉに代えて切換部４４から等高線図を図２７（ｃ）に示すように表示する。

【0108】

従って、オルソフォト画像又は等高線図からハイブリッド断層線を確認できることになる。

【0109】

＜赤色立体化画像＞
以下に赤色立体化画像ＲＧｉの生成方法を説明する。

【0110】

赤色立体化画像ＲＧｉの生成においては、開度という概念を用いている。開度は当該地点が周囲に比べて地上に突き出ている程度及び地下に食い込んでいる程度を数量化したものである。つまり、地上開度は、図２８（ａ）に示すように、着目する標本地点から距離ＬＲの範囲内で見える空の広さを表しており、また地下開度は図２８（ｂ）に示すように、逆立ちをして地中を見渡す時、距離ＬＲの範囲における地下の広さを表している。

【0111】

開度は距離ＬＲと周辺地形に依存している。図２９は９種の基本地形についての地上開度及び地下開度を、方位毎の地上角及び地下角の８角形グラフで示したものである。

【0112】

一般に地上開度は周囲から高く突き出ている地点ほど大きくなり、山頂や尾根では大きな値をとり窪地や谷底では小さい。逆に地下開度は地下に低く食い込んでいる地点ほど大きくなり、窪地や谷底では大きな値をとり山頂や尾根では小さい。実際には、距離ＬＲの範囲内でも種々の基本地形が混在しているために、地上角及び地下角の８角形グラフは変形され開度も種々の値をとることが多い。

【0113】

また、開度図は計算距離の指定によって、地形規模に適合した情報抽出が可能であり、方向性及び局所ノイズに依存しない表示が可能である。

【0114】

つまり、尾根線及び谷線の抽出に優れており、豊富な地形・地質情報が判読できるものであり、図３０に示すように、一定範囲のＤＥＭデータ上（地表面：立体：図３０の（ａ））において、設定した当該地点Ａから８方向のいずれか一方を見たときに最大頂点となる点Ｂを結ぶ直線Ｌ１と、水平線とがなす角度ベクトルθｉを求める。

【0115】

この角度ベクトルの求め方を８方向に渡って実施し、これらを平均化したものを地上開度と称し、一定範囲のＤＥＭデータ上（地表面：立体）に空気層を押し当てた立体（図３０の（ｂ））を裏返した反転ＤＥＭデータ（図３０の（ｃ））の当該地点Ａから８方向のいずれか一方を見たときに最大頂点となる点Ｃ（一番深い所に相当する）を結ぶ直線Ｌ２と、水平線とがなす角度を求める。この角度を８方向に渡って求めて平均化したのを地下開度ψｉと称している。

【0116】

すなわち、地上開度は、着目点から一定距離までの範囲に含まれるＤＥＭデータ上において、８方向毎に地形断面を生成し、それぞれの地点と着目点を結ぶ線（図３０の（ａ）のＬ１）の傾斜の最大値（鉛直方向から見たとき）を求める。このような処理を８方向に対して行う。

【0117】

傾斜の角度は天頂からの角度（平坦なら９０度、尾根や山頂では９０度以上、谷底や窪地では９０度以下）である。

【0118】

また、反転ＤＥＭデータの着目点から一定距離までの範囲において、８方向毎に地形断面を生成し、それぞれの地点と着目点を結ぶ線の傾斜の最大値（図３０の（ａ）の地表面の立体図において鉛直方向からＬ２を見たときには最小値）を求める。このような処理を８方向に対して行う。

【0119】

図３０の（ａ）の地表面の立体図において鉛直方向からＬ２を見たときの角度（地下開度ψｉ）は平坦なら９０度、尾根や山頂では９０度以下、谷底や窪地では９０度以上である。

【0120】

つまり、地上開度と地下開度は、図３１に示すように、２つの標本地点Ａ（ｉA，ｊA，ＨA）とＢ（ｉB，ｊB，ＨB）を考える。標本間隔が１ｍであることからＡとＢの距離は、
P = {(iA − iB)2 + (jA − jB)2}1/2 …(1)
となる。

【0121】

図３１の（ａ）は標高０ｍを基準として、標本地点のＡとＢの関係を示したものである
。標本地点Ａの標本地点Ｂに対する仰角θは
θ＝ｔａｎ-1｛（ＨB −ＨA ）／Ｐ
で与えられる。θの符号は（１）ＨA＜ＨB の場合には正となり、（２）ＨA＞ＨB の場合には負となる。

【0122】

着目する標本地点から方位Ｄの距離ＬＲの範囲内にある標本地点の集合を DＳL と記述して、これを「着目する標本地点のＤ−Ｌ集合」を呼ぶことにする。ここで、
DβL ：着目する標本地点の DＳL の各要素に対する仰角のうちの最大値
DδL ：着目する標本地点の DＳL の各要素に対する仰角のうちの最小値
として（図３１の（ｂ）参照）、次の定義をおこなう。

【0123】

定義１：着目する標本地点のＤ−Ｌ集合の地上角及び地下角とは、各々
DφL ＝９０− DβL
及び
DψL ＝９０＋ DδL
を意味するものとする。

【0124】

DφL は着目する標本地点から距離ＬＲ以内で方位Ｄの空を見ることができる天頂角の最大値を意味している。一般に言われる地平線角とはＬを無限大にした場合の地上角に相当している。また、DψL は着目する標本地点から距離ＬＲ以内で方位Ｄの地中を見ることができる天底角の最大値を意味している。Ｌを増大させると、 DＳL に属する標本地点の数は増加することから、 DβL に対して非減少特性を持ち、逆に DδL は非増加特性を持つ。したがって DφL 及び Dψ1 は共にＬに対して非増加特性を持つことになる。

【0125】

測量学における高角度とは、着目する標本地点を通過する水平面を基準にして定義される概念であり、θとは厳密には一致しない。また地上角及び地下角を厳密に議論しようとすれば、地球の曲率も考慮しなければならず、定義１は必ずしも正確な記述ではない。定義１はあくまでもＤＥＭを用いて地形解析をおこなうことを前提として定義された概念である。

【0126】

地上角及び地下角は指定された方位Ｄについての概念であったが、これを拡張したものとして、次の定義を導入する。

【0127】

定義II：着目する標本地点の距離ＬＲの地上開度及び地下開度とは、各々
ΦL＝（0φL ＋45φL ＋90φL ＋135φL ＋180φL ＋225φL ＋270φL ＋315φL ）／８
及び
ΨL＝（０ψL ＋45ψL ＋90ψL ＋135ψL ＋180ψL ＋225ψL ＋270ψL ＋315ψL ）／８
を意味するものとする。

【0128】

地上開度は着目する標本地点から距離ＬＲの範囲内で見える空の広さを表しており、また地下開度は逆立ちをして地中を見渡す時、距離ＬＲの範囲における地下の広さを表している（図３０参照）。

【0129】

そして、ＤＥＭデータを正方形にメッシュ化し、このメッシュ上の着目点と隣接する正方形の面の平均傾斜を求める。隣接する正方形は４通り存在しており、いずれか一つを着目正方形とする。そして、この着目正方形の４隅の高度と平均傾斜とを求める。平均傾斜は最小二乗法を用いて４点から近似した面の傾きである。

【0130】

そして、凸部強調により、尾根、谷底を明るさで表現するための第１のグレイスケールを備え、地上開度（着目点からＬの範囲を８方向見たときの、平均角度：高いところにいるかを判定するための指標）を求める毎に、この地上開度θｉの値に対応する明るさ（明度）を算出する。

【0131】

例えば、地上開度の値が４０度から１２０度程度の範囲に収まる場合は、５０度から１１０度を第１のグレイスケールに対応させ、２５５諧調に割り当てる。

【0132】

つまり、尾根の部分（凸部）の部分ほど地上開度の値が大きいので、色が白くなる。

【0133】

そして、着目点（座標）を有するメッシュ領域（ＤＥＭデータの同じＺ値を繋いだ等高線を正方形でメッシュ化し（例えば１ｍ）、このメッシュの４隅のいずれかの点を着目点としている場合）に、第１のグレイスケールに基づく色データを割り付け、これを保存する。次に、諧調補正によってこの保存された地上開度画像データの色諧調を反転させた地上開度レイヤーを保存する。つまり、尾根が白くなるように調整した地上開度レイヤーを得ている。

【0134】

次に、谷底、尾根を明るさで表現するための第２のグレイスケールを備え、地下開度ψｉ（着目点から８方向の平均）を求める毎に、この地上開度の値に対応する明るさを算出する。

【0135】

例えば、地下開度の値が４０度から１２０度程度の範囲に収まる場合は、５０度から１１０度を第２のグレイスケールに対応させ、２５５諧調に割り当てる。つまり、谷底の部分（凹部）の部分ほど地下開度の値が大きいので、色が黒くなることになる。

【0136】

そして、地下開度画像データを読み、着目点（座標）を有するメッシュ領域（ＤＥＭデータの同じＺ値を繋いだ等高線を正方形でメッシュ化し（例えば１ｍ）、このメッシュの４隅のいずれかの点を着目点としている場合）に、第２のグレイスケールに基づく色データを割り付け、これを保存する。次に、諧調補正処理によって地下開度画像データの色諧調を補正する。

【0137】

色が黒くなり過ぎた場合は、トーンカーブを補正した度合いの色にする。これを地下開度レイヤーと称して保存する。

【0138】

次に、傾斜の度合いを明るさで表現するに応じたで表現するための第３のグレイスケールを備え、傾斜度（着目点から４方向の平均）を求める毎に、この傾斜度の値に対応する第３のグレイスケールの明るさ（明度）を算出する。

【0139】

例えば、斜度αｉの値が０度から７０度程度の範囲に収まる場合は、０度から５０度を第３のグレイスケールに対応させ、２５５諧調に割り当てる。つまり、０度が白、５０度以上が黒。傾斜α（斜度αｉともいう）の大きい地点ほど色が黒くなる。

【0140】

そして、地下開度画像データと地上開度画像データとの差画像を斜度画像として保存する。

【0141】

このとき、着目点（座標）を有するメッシュ領域（ＤＥＭデータの同じＺ値を繋いだ等高線を正方形でメッシュ化し（例えば１ｍ）、このメッシュの４隅のいずれかの点を着目点としている場合）に、第３のグレイスケールに基づく色データを割り付ける。次に、赤色処理がＲＧＢカラーモード機能でＲを強調する。つまり、傾斜が大きいほど赤が強調された傾斜強調画像を得る。

【0142】

そして、地上開度レイヤーと地下開度レイヤーとを乗算して合成した合成画像を得る。このとき、谷の部分が潰れないように両方のバランスを調整する。

【0143】

前述の「乗算」というのは、フォトショップ（photoshop）上のレイヤーモードの用語で、数値処理上はＯＲ演算となる。

【表1】

【0144】

このバランス調整は、地上開度と地下開度の値の配分は、ある地点を中心として一定の半径（Ｌ／２）の地表面を切り取る。

【0145】

空全体が一様な明るさの場合に地表面から見上げる空の広さが地面の明るさを与える。

【0146】

つまり、地上開度が明るさとなる。しかし、光が回り込むことまで考えると、地下開度の値も考慮するべきである。

【0147】

この両者の比をどのようにするべきかで、地形の尾根の部分を強調したり、任意に変化させることができる。谷の中の地形を強調したいときはｂの値を大きくする。

【0148】

明るさの指標＝ａ×地上開度−ｂ×地下開度
但し、ａ＋ｂ＝１
すなわち、地上開度レイヤー（尾根を白強調）と地下開度レイヤー（底を黒く強調）と乗算合成した灰色の諧調表現の合成画像を得る。

【0149】

一方、傾斜強調画像と合成画像と合成した尾根が赤色で強調された赤色立体地図ＲＧｉを得て、表示部２００に表示する。

【0150】

すなわち、尾根や山頂部分が白っぽく、谷や窪地が黒っぽく表現し、傾斜が急な部分ほど赤く表現した立体感のある画像が生成できる。このため、一目で凹凸の高低の度合い及び傾斜の度合いを把握させることができる。

【0151】

つまり、ＤＥＭデータを連結する地表面の着目点のＤＥＭデ-タから複数方向毎に、一定範囲内までの最大頂点となるＤＥＭデータと水平線とがなす角度ベクトルをそれぞれ求めて平均化した地上開度を求め、この地上開度の値の大きさほどに明るい色を割りあてた地上開度画像を得る。

【0152】

そして、一定範囲のＤＥＭデータ上に空気層を押し当てた立体を裏返した反転ＤＥＭデータの着目点のＤＥＭデータから複数方向毎に、一定範囲内までの最大頂点となるＤＥＭデータと水平線とがなす角度ベクトルをそれぞれ求めて平均化して地下開度を求め、この地下開度の値の大きさほどに暗い色を割りあてた地下開度画像を得る。

【0153】

そして、地上開度画像と地下開度画像とを重み付け合成し、この値に応じて階調表現した合成画像を赤色立体化画像ＲＧｉとして得ている。

【0154】

＜実施の形態３＞
図３２は実施の形態３の断層図作成装置の概略構成図である。図３２においては図１と同一符号のものは説明を省略する。
図３２に示すように、実施の形態３の断層図作成装置は、モード判定部５０と、メモリ定義部５２と、既存・赤色判読断層線描画部５４と、レイヤ順位設定部５６と、γ線測定点定義部５８と、推定断層線描画部６０と、ハイブリッド断層線描画部６２と、画像合成出力部６４等からなるプログラム構成を備えている。

【0155】

なお、画像化部１４は、画像メモリ（図示せず）を備え、一定時間毎にこの画像メモリに保存された各種画像を合成（重ねて）して表示部２００に表示する。このとき、画像に付加されている色情報、線情報等で画面に表示する。

【0156】

モード判定部５０は、画面に表示する画面を入力させるためのダイアログボックス（図示せず）を表示させて、このダイアログボックスに入力された情報に基づいて該当する各部を起動させる。
このダイアログボックスは、例えば、メモリ３２の赤色立体地図ＲＧｉを表示させるボタン、既存断層線用メモリ３６の既存断層線ＫＬｊを表示させるためのボタン、微細断層線定義用メモリ３４の微細断層線ＤＬｊを表示させるためのボタン、ハイブリッド断層線ＨＢＬｊを表示させるためのボタンを備えている。

【0157】

また、γ線測定点Ｐｉを定義させるモードを起動させるボタン、γ線現地測定結果直線ＬＧｊを作成するモードを起動させるためのボタン、γ線現地測定結果直線ＬＧｊを移動させるモードを起動させるためのボタン、ハイブリッド断層線ＨＢＬｊを作成するためのモードを起動させるためのボタン等を備えている。
さらに、後述する既存断層線ワーキング用メモリ７１の既存断層線ＫＬｊを削除するめのボタン、後述する微細断層線ワーキング用メモリ７３の微細断層線ＤＬｊを削除するためのボタン等を備えている。

【0158】

メモリ定義部５２は、起動に伴って赤色立体化画像ＲＧｉの原点を読み込み、この原点を有する微細断層線ワーキング用メモリ７３、既存断層線ワーキング用メモリ７１、γ線現地測定地点定義用ワーキングメモリ７５、γ線現地測定結果直線ワーキング用メモリ７６、ハイブリッド断層線ワーキング用メモリ７８をワーキング用メモリ部７０に生成して、各々が前述の原点を有するＸ―Ｙ平面を定義する。
つまり、メモリ定義部５２は、範囲定義部１６における赤色立体化画像ＲＧｉの原点座標を各平面に定義している。

【0159】

微細断層線ワーキング用メモリ７３に定義された平面を赤色判読結果微細断層線ワーキング用平面Ｈａと称し、既存断層線ワーキング用メモリ７１に定義された平面を既存断層線ワーキング用平面Ｈｂと称し、γ線現地測定地点定義用ワーキングメモリ７５に定義された平面をγ線現地測定地点定義ワーキング用平面Ｈｃと称し、γ線現地測定結果直線ワーキング用メモリ７６に定義された平面をγ線現地測定結果直線描画ワーキング用平面Ｈｄと称し、ハイブリッド断層線ワーキング用メモリ７８に定義された平面をハイブリッド断層線描画ワーキング用平面Ｈｅと称する。

【0160】

既存・赤色判読断層線描画部５４は、γ線微細断層線ワーキング用メモリ７３の生成に伴って、微細断層線定義用メモリ３４に記憶されている赤色立体化画像ＲＧｉを画面に表示して断層と判断された微細断層線ＤＬｊ（ＤＬ１、ＤＬ２・・・）を微細断層線ワーキング用メモリ７３に複写（コピー）する。
また、既存・赤色判読断層線描画部５４は、既存断層線ワーキング用メモリ７１の生成に伴って、既存断層線用メモリ３６に記憶されている既存用平面に描かれている既存断層線ＫＬｊ（ＫＬ１、ＫＬ２）を既存断層線ワーキング用メモリ７１に複写する。

【0161】

レイヤ順位設定部５６は、入力された順位に基づいて画面に表示させる画像の階層（レイヤー）の順番を図示しないメモリに設定して画像化部１４によってその階層の順番で画像を重ね表示させる。

【0162】

この階層は、図３３に示すように、例えば第１層（最上位）がγ線現地測定地点定義用ワーキングメモリ７５（γ線現地測定地点定義ワーキング用平面Ｈｃ）にされ（具体的にはメモリアドレスを関連付けるのが好ましい）、第２層がγ線現地測定結果直線ワーキング用メモリ７６（γ線現地測定結果直性描画ワーキング用平面Ｈｄ）にされ（具体的にはメモリアドレスを関連付けるが好ましい）、第３層がハイブリッド断層線ワーキング用メモリ７８（ハイブリッド断層線描画ワーキング用平面Ｈｅ）にされ（（具体的にはメモリアドレスを関連付けるのが好ましい）、第４層がγ線微細断層線ワーキング用メモリ７３（赤色判読結果微細断層線ワーキング用平面Ｈａ）にされ（具体的にはメモリアドレスを関連付けるのが好ましい）、第５層が既存断層線ワーキング用メモリ７１（既存断層線ワーキング用平面Ｈｂ）にされ（具体的にはメモリアドレスが好ましい）ている。

【0163】

そして、第６層は赤色立体化画像ＲＧｉ（３２）にされている（具体的にはメモリアドレスが好ましい）。
γ線測定点定義部５８は、起動に伴って、γ線測定点Ｐｉの色、記号種類（例えば〇）を入力させる画面（図示せず）を表示し、γ線測定点Ｐｉの識別番号（以下γ線測定点番号ＰＢｉという）を生成し、これらを関連付けてメモリ８０に記憶する（図３４参照）。

【0164】

そして、γ線測定点定義部５８はレイヤ順位設定部５６にγ線現地測定地点定義用ワーキングメモリ７５のメモリアドレスを読み込み、このメモリアドレスのメモリをワーキング用メモリ部７０から引き当てる。このγ線現地測定地点定義用ワーキングメモリ７５に定義されているγ線現地測定地点定義ワーキング用平面Ｈｃを画像化部１４に出力して画面（γ線測定点Ｐｉが定義されていない場合は透明）に表示させる。

【0165】

そして、画面上において、カーソルが位置され、これが例えばクリックされると、γ線測定点Ｐｉの決定と判断して、この画面上のカーソルの位置をγ線現地測定地点定義ワーキング用平面Ｈｃの座標に変換して現地測定点座標ＰＺｉとして定義（計算）する。
つまり、γ線測定点定義部５８は、γ線測定点・断層線描画部２４における
赤色立体化画像ＲＧｉにカーソルで指示された点をγ線測定点Ｐｉとして読み込み、このγ線測定点Ｐｉをγ線現地測定地点定義ワーキング用平面Ｈｃに定義する機能を行っている。
このとき、現地測定点座標ＰＺｉに入力されたγ線測定点Ｐｉの色、記号種類（例えば〇）、γ線測定点番号ＰＢｉ等を割り付ける（図３５参照）。

【0166】

そして、γ線測定点定義部５８は画像化部１４に、γ線現地測定地点定義ワーキング用平面Ｈｃに定義されたγ線測定点Ｐｉの現地測定点座標ＰＺｉと色（例えばイエロー）、大きさ、記号種類（例えば〇）等を組みにして出力する。
従って、図１９に示すように画面にγ線測定点Ｐｉが赤色立体化画像ＲＧｉ上に表示される。図１９は既に既存断層線ＫＬｊ、微細断層線ＤＬｊが表示されている例である。また、γ線現地測定地点定義用ワーキングメモリ７５が最上位にされている例である。

【0167】

このとき、メモリ８０に現地測定点座標ＰＺｉをγ線測定点番号ＰＢｉに関連付けて記憶する（図３４参照）。
推定断層線描画部６０は、起動に伴って、γ線現地測定結果直線ＬＧｊを描かせるための色、線種（太さ）等を入力させる画面（図示せず）を表示してγ線現地測定結果直線ＬＧｊの線番号ＬＧＢｉ（ユニークコード：年月日、時刻、座標を含む）を生成してγ線測定点番号ＰＢｉに関連付けてメモリ８０に記憶する（図３４参照）。

【0168】

そして、推定断層線描画部６０はγ線測定点番号ＰＢｉに関連付けられている現地測定点座標ＰＺｉをγ線現地測定結果直線ワーキング用メモリ７６（γ線現地測定結果直線描画ワーキング用平面Ｈｄ）に定義する（図３５参照）。
そして、メモリ８０に記憶されたγ線測定点Ｐｉのγ線測定点番号ＰＢｉ又は線番号ＬＧＢｉを有する現地γ線測定情報Ｆγｉを現地測定結果情報用メモリ４０から読み込み、この現地γ線測定情報Ｆγｉに含まれている方向Ｌθｉ、断層長Ｌ、幅等を読み込む。

【0169】

そして、推定断層線描画部６０は、レイヤ順位設定部５６に設定されているγ線現地測定結果直線ワーキング用メモリ７６（γ線現地測定結果直線描画ワーキング用平面Ｈｄ）のメモリアドレスを読み込み、このメモリアドレスのメモリをワーキング用メモリ部７０から引き当てる。
そして、このγ線現地測定結果直線ワーキング用メモリ７６のγ線現地測定結果直線描画ワーキング用平面Ｈｄ上のγ線測定点Ｐｉの現地測定点座標ＰＺｉを起点として方向Ｌθｉ、断層長Ｌ（例えば３８００ｍ）、幅Ｗのγ線現地測定結果直線ＬＧｊを描画（計算）する（図３５参照）。このとき、γ線現地測定結果直線ＬＧｊには線番号ＬＧＢｉと方向Ｌθｉと断層長Ｌ（例えば３８００ｍ）と幅Ｗ、線種、色、現地測定点座標ＰＺｉ等を関連付ける。

【0170】

そして、推定断層線描画部６０は、このγ線現地測定結果直線描画ワーキング用平面Ｈｄに描かれたγ線現地測定結果直線ＬＧｊに関連付けられている方向Ｌθｉと断層長Ｌ（例えば３８００ｍ）と幅Ｗ、線種、色、現地測定点座標ＰＺｉ等を線属性情報ＺＬｉとして画像化部１４に出力する。従って、図２２に示すように表示される。

【0171】

さらに、推定断層線描画部６０は、γ線現地測定結果直線ＬＧｊを移動させるモードの場合は、レイヤ順位設定部５６からγ線現地測定結果直線ＬＧｊが定義されているメモリアドレスを読み込む。つまり、γ線現地測定結果直線ワーキング用メモリ７６（γ線現地測定結果直線描画ワーキング用平面Ｈｄ）のメモリアドレスを読み込む。 そして、このγ線現地測定結果直線描画ワーキング用平面Ｈｄのγ線現地測定結果直線ＬＧｊに関連図蹴られている線属性情報ＺＬｉを読み込む。

【0172】

そして、推定断層線描画部６０は、カーソルの移動に伴って読み込んだ線属性情報ＺＬｉの方向Ｌθｉと逆方向にγ線現地測定結果直線ＬＧｊを移動させる（図３６参照）。このとき、一定時間毎にγ線現地測定結果直線描画ワーキング用平面Ｈｄのγ線現地測定結果直線ＬＧｊを画像化部１４に出力している。
つまり、推定断層線描画部６０は、γ線測定点・断層線描画部２４におけるγ線現地測定結果直線ＬＧｊを表示させ、カーソルの移動に伴ってγ線現地測定結果直線ＬＧｊの移動を方向Ｌθｉに対して逆方向に移動させる機能を行っている。

【0173】

従って、図２３に示すように、γ線現地測定結果直線ＬＧｊが表示されることになる。

【0174】

ハイブリッド断層線描画部６２は、起動に伴って、ハイブリッド断層線ＨＢＬｊを描くための色（例えばピンク）、線の太さ等を入力させる画面（図示せず）を表示させる。そして、ハイブリッド断層線ＨＢＬｊの識別番号（以下ハイブリッド断層線番号ＨＢｂｉという）を生成し、これらを関連付けてメモリ９０に記憶する（図３７参照）。そして、レイヤ順位設定部５６のハイブリッド断層線ＨＢＬｊが定義されているメモリアドレスを読み込む。つまり、ハイブリッド断層線ワーキング用メモリ７８（ハイブリッド断層線描画ワーキング用平面Ｈｅ）を引き当てる。

【0175】

そして、このハイブリッド断層線描画ワーキング用平面Ｈｅを画像化部１４によって表示させる（透明）。
そして、画面のカーソルの移動（微細断層線ＤＬｊ上）に伴って、このカーソルの座標をハイブリッド断層線描画ワーキング用平面Ｈｅの座標に変換してハイブリッド断層線描画ワーキング用平面Ｈｅに定義（計算）していき（図３８参照）、この座標と色（ピンク）と太さ等を画像化部１４に出力する。
つまり、ハイブリッド断層線描画部６２は、線調整部２６におけるハイブリッド断層線ＨＢＬｊ（ピンク）として表示させる機能を行っている。

【0176】

従って、図２４に示すように、微細断層線ＤＬｊ同士がなぞられながら線で繋げられていくことになる。

【0177】

画像合成出力部６４は、メモリ３２の赤色立体地図ＲＧｉに、γ線微細断層線ワーキング用メモリ７３の赤色判読結果微細断層線ワーキング用平面Ｈａ、既存断層線ワーキング用メモリ７１の既存断層線ワーキング用平面Ｈｂ、γ線現地測定地点定義用ワーキングメモリ７５のγ線現地測定地点定義ワーキング用平面Ｈｃ、γ線現地測定結果直線ワーキング用メモリ７６のγ線現地測定結果直線描画ワーキング用平面Ｈｄ、ハイブリッド断層線ワーキング用メモリ７８のハイブリッド断層線描画ワーキング用平面Ｈｅのいずれか又は全ての画像を重ね合わせてこれをプリンタ（図示せず）に出力して印刷させる。

【0178】

（ハイブリッド断層図ＨＨＧｉ）
ハイブリッ断層図ＨＨＧｉは、既存断層線ＫＬｊを削除するめのボタン、微細断層線ＤＬｊを削除するためのボタン、γ線測定点Ｐｉを削除するボタン、γ線現地測定結果直線ＬＧｊを削除するボタンの選択によって得る。

【0179】

すなわち、γ線測定点定義部５８は、γ線測定点Ｐｉの削除モードの場合は、γ線現地測定地点定義用ワーキングメモリ７５を引き当て、画面のγ線測定点Ｐｉが指定されるとγ線現地測定地点定義ワーキング用平面Ｈｃに定義されているγ線測定点Ｐｉを削除する。

【0180】

既存・赤色判読断層線描画部５４は、既存断層線ＫＬｊの削除モードの場合は、既存断層線ワーキング用メモリ７１を引き当てる。そして、カーソルで画面の既存断層線ＫＬｊがなぞられるに伴って、この座標に該当する既存断層線ワーキング用メモリ７１の既存断層線ワーキング用平面Ｈｂ上の既存断層線ＫＬｊを削除する。

【0181】

また、既存・赤色判読断層線描画部５４は、微細断層線ＤＬｊの削除モードの場合は、γ線微細断層線ワーキング用メモリ７３を引き当てる。そして、カーソルで画面の微細断層線ＤＬｊがなぞられるに伴って、この座標に該当するγ線微細断層線ワーキング用メモリ７３の赤色判読結果微細断層線ワーキング用平面Ｈａ上の微細断層線ＤＬｊを削除する。

【0182】

また、推定断層線描画部６０は、γ線現地測定結果直線ＬＧｊの削除モードの場合は、γ線現地測定結果直線ワーキング用メモリ７６を引き当てる。そして、カーソルで画面のγ線現地測定結果直線ＬＧｊがなぞられるに伴って、この座標に該当するγ線現地測定結果直線ワーキング用メモリ７６のγ線現地測定結果直線ＬＧｊを削除する。
このような処理によって、図２５に示すハイブリッド断層図ＨＨＧｉを得ている。

【0183】

また、レイヤ順位設定部５６において図３９に示すように、第１層をγ線現地測定結果直線ワーキング用メモリ７６（γ線現地測定結果直線描画ワーキング用平面Ｈｄ）とした場合は、図４０に示すようにγ線測定点Ｐｉ上をγ線現地測定結果直線ＬＧｊが通ることになる。

【0184】

なお、上記実施の形態ではオルソフォト画像、等高線図にハイブリッド断層線を重ね表示するとしたが、高度断彩図、陰影図であってもかまわない。

【0185】

また、上記の各実施の形態においては微細断層として説明したが、幅が１５ｍ、２０ｍ、・・・の断層であってもかまわない。
また、各線の色は任意であるが、各々が異なる色又は太さ若しくは線種であることが好ましい。

【符号の説明】

【0186】

１０ 赤色地図作成部
１２ 断層検索範囲設定部
１４ 画像化部
１６ 範囲定義部
１８ 既存・赤色判読描画部
２４ γ線測定点・断層線描画部
２６ 線調整部

【要約】

【課題】幅数メートル程度の断層を画面から容易に判読して、この判読結果に基づく断層線図を容易に得ることができる断層図作成装置を得る。
【解決手段】断層図作成装置は、赤色立体化画像作成部１０と、断層用検索範囲設定部１２と、画像化部１４と、範囲定義部１６と、既存・赤色判読描画部１８と、γ線測定点・断層線描画部２４と、線調整部２６等を備えて、現地での測定結果に基づくγ線現地測定結果直線ＬＧｊ（黄色）が赤色立体化画像ＲＧｉ上から判読してこの赤色立体化画像ＲＧｉに描いた複数の微細断層線ＤＬｊが並ぶ線にほぼ平行で、かつこれらの微細断層線ＤＬｊに渡る区間に容易に対応させることができるγ線現地測定結果直線ＬＧｊを表示する。
【選択図】図１

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図16】

【図17】

【図18】

【図19】

【図20】

【図21】

【図22】

【図23】

【図24】

【図25】

【図26】

【図27】

【図28】

【図29】

【図30】

【図31】

【図32】

【図33】

【図34】

【図35】

【図36】

【図37】

【図38】

【図39】

【図40】