(57)【特許請求の範囲】
【請求項１】
地図画像を複数のメッシュ状の領域に分解し、当該領域各々の緯度経度及び高度からなる３次元位置情報の示す位置の線量率から、前記領域各々における線量率の放射線を放射する放射線源の放射能を算出し、算出した放射能から周辺の他の領域における線量率への３次元空間における寄与率を算出し、前記放射能、前記寄与率及び外部から供給される除染率により、当該除染率による除染後の前記領域毎の線量率を求める除染処理推定サーバと、
線量率測定器から供給される、当該線量率測定器の配置されている前記３次元位置情報の示す位置の前記線量率を、当該３次元位置情報に対応させて測定線量率データベースに書き込んで記憶させるデータベース管理サーバと
を備え、
前記除染処理推定サーバが、
前記測定線量率データベースの前記線量率の３次元位置情報における緯度経度が含まれる緯度経度範囲の前記領域の線量率として、当該領域の前記放射線源の放射能を算出し、前記地図画像の前記領域に対応する線量率が前記測定線量率データベースに記憶されていない場合、前記データベース管理サーバに対し、前記測定線量率データベースに記憶されていない前記領域の前記線量率を要求する線量率要求情報を出力し、
前記データベース管理サーバが、
インターネットにおける線量率を公開する専用サイトから、前記線量率要求情報の供給する緯度経度範囲における線量率を読み込み、公開線量率データベースに対し、緯度経度範囲に対応して書き込んで記憶させる
ことを特徴とする３次元線量評価マッピングシステム。
【請求項２】
前記除染処理推定サーバが、
前記測定線量率データベースと前記公開線量率データベースとから、前記地図画像における前記領域の線量率を読み出し、前記除染率が設定された領域の前記放射能の他の領域に対する寄与率である推定寄与率を求め、当該推定寄与率及び前記放射能により、前記地図画像における前記領域各々の除染後の線量率を算出する
ことを特徴とする請求項１に記載の３次元線量評価マッピングシステム。
【請求項３】
前記除染処理推定サーバが、
前記メッシュ状の前記領域に平面視で重なるメッシュからなる半透明のメッシュ画像を生成し、前記線量率の数値範囲毎に異なる色を前記メッシュに付与し、前記地図画像に対してこのメッシュ画像を重ねて表示部に表示させる
ことを特徴とする請求項１または請求項２に記載の３次元線量評価マッピングシステム。
【請求項４】
地図画像を複数のメッシュ状の領域に分解する過程と、
当該領域各々の緯度経度及び高度からなる３次元位置情報の示す位置の線量率から、前記領域各々における線量率の放射線を放射する放射線源の放射能を算出する過程と、
算出した放射能から周辺の他の領域における線量率への３次元空間における寄与率を算出する過程と、
前記放射能、前記寄与率及び外部から供給される除染率により、当該除染率による除染後の前記領域毎の線量率を求める過程と、
線量率測定器から供給される、当該線量率測定器の配置されている前記３次元位置情報の示す位置の前記線量率を、当該３次元位置情報に対応させて測定線量率データベースに書き込んで記憶させる過程と、
前記測定線量率データベースの前記線量率の３次元位置情報における緯度経度が含まれる緯度経度範囲の前記領域の線量率として、当該領域の前記放射線源の放射能を算出する過程と、
前記地図画像の前記領域に対応する線量率が前記測定線量率データベースに記憶されていない場合、前記測定線量率データベースに記憶されていない前記領域の前記線量率を要求する線量率要求情報を出力する過程と、
インターネットにおける線量率を公開する専用サイトから、前記線量率要求情報の供給する緯度経度範囲における線量率を読み込み、公開線量率データベースに対し、緯度経度範囲に対応して書き込んで記憶させる過程と
を含むことを特徴とする３次元線量評価マッピング方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、放射性物質が放射する放射線量を評価する３次元線量評価マッピングシステム及びその方法に関する。
【背景技術】
【０００２】
放射性物質である放射線源（以下、線源とする）から、この線源の放射能に対応した放射線量（以下、線量とする）の放射線が放射される。
そして、ある地域において、この線量を測定することにより、その地域での人体に対する放射線の線量当量（単位：シーベルト）を測定した線量から算出することができる。
この線量当量による放射線の人体の各組織及び臓器毎に受ける影響が知られており、国が決めた基準を超える線量当量となる場合、その線量の地域への立ち入りが制限されることになる。
【０００３】
従来においては、表計算ソフトなどを用いて、測定された線量データを、作業者が測定した地点を示す地点情報（例えば、メッシュ形状の緯度経度情報）に対応付けて、コンピュータに入力している。
この場合、表計算ソフトで用いることができる形式のデータファイルを用意して、データファイルから線量データを読み込み、線量分布を作成している。このとき、作業者がデータファイルにおいて地図の位置と線量データを特定した位置との対応付けを、地図データにおける縮尺を含めて行う必要がある。
また、上述した問題を解決するため、構造物の地点毎の線量をマッピングする放射線線量率マップ作成装置が開示されている（例えば、特許文献１）。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開平０６−１８６３４０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、特許文献１においては、放射性物質の除染を行った後の線量が推定できないため、放射性物質を除染した後に、この除染によりどの程度の線量低下があるかの効果を得ることができない。
また、線源分布を求めて、線量の推定は行われているが、この推定は距離による放射線の減衰補正が主体となっており、除染を行った地域周辺の３次元的な地形を考慮しているわけではない。
このため、実際の地形に即した線源分布による周辺の寄与が不明であり、正確な除染の効果を推定することができない。
【０００６】
本発明は、このような事情に鑑みてなされたもので、線源の分布を明確に推定して、除染後の線源の量に即した線量を、地形に対応させた３次元的な線源の配置から求める３次元線量評価マッピングシステム及びその方法を提供することを目的とする。
【課題を解決するための手段】
【０００７】
この発明は上述した課題を解決するためになされたもので、本発明の３次元線量評価マッピングシステムは、地図画像を複数のメッシュ状の領域に分解し、当該領域各々の緯度経度及び高度からなる３次元位置情報の示す位置の線量率から、前記領域各々における線量率の放射線を放射する放射線源の放射能を算出し、算出した放射能から周辺の他の領域における線量率への３次元空間における寄与率を算出し、前記放射能、前記寄与率及び外部から供給される除染率により、当該除染率による除染後の前記領域毎の線量率を求める除染処理推定サーバと、線量率測定器から供給される、当該線量率測定器の配置されている前記３次元位置情報の示す位置の前記線量率を、当該３次元位置情報に対応させて測定線量率データベースに書き込んで記憶させるデータベース管理サーバとを備え、 前記除染処理推定サーバが、前記測定線量率データベースの前記線量率の３次元位置情報における緯度経度が含まれる緯度経度範囲の前記領域の線量率として、当該領域の前記放射線源の放射能を算出し、前記地図画像の前記領域に対応する線量率が前記測定線量率データベースに記憶されていない場合、前記データベース管理サーバに対し、前記測定線量率データベースに記憶されていない前記領域の前記線量率を要求する線量率要求情報を出力し、前記データベース管理サーバが、インターネットにおける線量率を公開する専用サイトから、前記線量率要求情報の供給する緯度経度範囲における線量率を読み込み、公開線量率データベースに対し、緯度経度範囲に対応して書き込んで記憶させることを特徴とする。
【００１０】
本発明の３次元線量評価マッピングシステムは、前記除染処理推定サーバが、前記測定線量率データベースと前記公開線量率データベースとから、前記地図画像における前記領域の線量率を読み出し、前記除染率が設定された領域の前記放射能の他の領域に対する寄与率である推定寄与率を求め、当該推定寄与率及び前記放射能により、前記地図画像における前記領域各々の除染後の線量率を算出することを特徴とする。
【００１１】
本発明の３次元線量評価マッピングシステムは、前記除染処理推定サーバが、前記メッシュ状の前記領域に平面視で重なるメッシュからなる半透明のメッシュ画像を生成し、前記線量率の数値範囲毎に異なる色を前記メッシュに付与し、前記地図画像に対してこのメッシュ画像を重ねて表示部に表示させることを特徴とする。
【００１２】
本発明の３次元線量評価マッピング方法は、地図画像を複数のメッシュ状の領域に分解する過程と、当該領域各々の緯度経度及び高度からなる３次元位置情報の示す位置の線量率から、前記領域各々における線量率の放射線を放射する放射線源の放射能を算出する過程と、算出した放射能から周辺の他の領域における線量率への３次元空間における寄与率を算出する過程と、前記放射能、前記寄与率及び外部から供給される除染率により、当該除染率による除染後の前記領域毎の線量率を求める過程と、線量率測定器から供給される、当該線量率測定器の配置されている前記３次元位置情報の示す位置の前記線量率を、当該３次元位置情報に対応させて測定線量率データベースに書き込んで記憶させる過程と、前記測定線量率データベースの前記線量率の３次元位置情報における緯度経度が含まれる緯度経度範囲の前記領域の線量率として、当該領域の前記放射線源の放射能を算出する過程と、前記地図画像の前記領域に対応する線量率が前記測定線量率データベースに記憶されていない場合、前記測定線量率データベースに記憶されていない前記領域の前記線量率を要求する線量率要求情報を出力する過程と、インターネットにおける線量率を公開する専用サイトから、前記線量率要求情報の供給する緯度経度範囲における線量率を読み込み、公開線量率データベースに対し、緯度経度範囲に対応して書き込んで記憶させる過程とを含むことを特徴とする。
【発明の効果】
【００１３】
この発明によれば、測定した線量から、各地点毎の線源の３次元座標における分布を明確に推定した後、この線源の分布に対応して、除染後の線源の放射能に即した線量を推定するため、除染後の３次元空間における線量を正確に推定することができる。
【図面の簡単な説明】
【００１４】
【図１】この発明の一実施形態による３次元線量評価マッピングシステムの構成例を示す概略ブロック図である。
【図２】測定線量率データベース１３に記憶されている線量率テーブルの構成例を示す図である。
【図３】統合線量率データベース１５に記憶されている統合線量率テーブルの構成例を示す図である。
【図４】地図画像上に重ねられるメッシュ画像における各メッシュの線量率の大きさを色分けした例を示す図である。
【図５】放射線源データベース１６に記憶されている放射線源テーブルの構成例を示す図である。
【図６】本実施形態による３次元線量評価マッピングシステム１における、除染した地域とこの地域の周囲の他の地域との推定線量率を算出する動作例を示すタイミングチャートである。
【図７】測定線量率データベース１３の線量率テーブルに記憶されている３次元位置線量率データの測定位置を表示した図である。
【図８】公開線量率データベース１４に記憶されている３次元位置線量率データの測定位置を表示した図である。
【図９】図４に示すメッシュ画像と地図画像との対応関係によって、実際にシミュレーションで仮想的に求めた線量率をメッシュ画像に反映させ、地図画像に重ね合わせて表示した際の表示画像を示す図である。
【発明を実施するための形態】
【００１５】
以下、図面を参照して、本発明の実施の形態について説明する。図１は、この発明の一実施形態による３次元線量評価システムの構成例を示す概略ブロック図である。
本実施形態による３次元線量評価マッピングシステムは、３次元線量評価マッピングシステム１と、線量測定装置２＿１から２＿ｎと、データ収集装置３とを備えている。
線量測定装置２＿１から２＿ｎの各々は、それぞれの配置した位置における、放射線の線量率、緯度経度及び高度を含む線量データを測定する。
データ収集装置３は、予め設定された周期により、各地点に配置されている線量測定装置２＿１から２＿ｎの各々より３次元位置線量データ（後述）を収集し、収集した３次元位置線量データに対して、この３次元位置線量データを識別する測定識別情報、例えば読み込んだ順番を付加し、３次元線量評価マッピングシステム１へ送信する。
３次元線量評価マッピングシステム１は、データ収集装置３を介して供給される、外部装置である線量測定装置２＿１〜２＿ｎの３次元位置線量データを基に、地図画像を所定の分解率によって、複数のメッシュ上の領域に分割し、このメッシュ状の領域における線源の放射能を推定し、この線源の除染の程度に応じて、各メッシュにおける除染後の線量率を推定する。
【００１６】
また、線量測定装置２＿１から２＿ｎの各々は、線量率測定器２１とＧＰＳ（Ｇｌｏｂａｌ  Ｐｏｓｉｔｉｏｎｉｎｇ  Ｓｙｓｔｅｍ：全地球測位システム）・高度測定器２２とを備えている。
線量率測定器２１は、自身の配置されている地点の線量率を測定し、この線量率と、測定した測定日時と、この線量率の測定に要した時間の長さを示す測定時間と、自身の識別情報である測定器識別情報とともに、線量率データとして出力する。
ＧＰＳ・高度測定器２２は、自身の配置されている緯度経度と高度とを測定し、緯度経度及び高度を含む３次元位置情報とし、この３次元位置情報と上記線量率データとを合わせて３次元位置線量率データとして、３次元線量評価マッピングシステム１へ出力する。
【００１７】
また、３次元線量評価マッピングシステム１は、データベース管理サーバ１１、除染処理推定サーバ１２、測定線量率データベース１３、公開線量率データベース１４、統合線量率データベース１５、放射線源データベース１６及び表示部１７を備えている。
次に、図２は、データベース管理サーバ１１が記憶する線量率テーブルの構成を示す図である。この図２において、線量率テーブルは、測定識別情報と、測定日時と、３次元位置情報と、線量率と、測定時間と、測定器識別情報とから構成されている。
図１に戻り、データベース管理サーバ１１は、データ収集装置３から供給される３次元位置線量率データを、測定識別情報毎に、測定線量率データベース１３における線量率テーブルに順次書き込んで記憶させる。
【００１８】
除染処理推定サーバ１２は、表示部１７に表示される地図画像上に、この地図画像に重ねるための予め設定された分解率により地図画像の領域をメッシュ状の領域に分割する。
そして、除染処理推定サーバ１２は、この地図画像のメッシュ状の領域の各々に対して、それぞれが平面視で重なるメッシュからなるメッシュ画像を生成する。
また、除染処理推定サーバ１２は、この地図画像における緯度軽度範囲に重なるメッシュ画像（半透明状態としたフィルタ画像とする、色が付与されても半透明状態となる）に対し、このメッシュの位置に含まれる緯度経度に対応する線量率データを、３次元位置情報における緯度経度により測定線量率データベース１３から読み出す。
また、除染処理推定サーバ１２は、メッシュ画像における、対応する３次元位置情報のないメッシュの位置の線量率を、メッシュを識別するメッシュ識別情報と、このメッシュ識別情報の示すメッシュの緯度経度範囲とを付加した線量率要求情報を、データベース管理サーバ１１に対して送信する。
【００１９】
データベース管理サーバ１１は、除染処理推定サーバ１２からの線量率要求情報に対応して、要求されたメッシュに対する線量率を、このメッシュの緯度経度範囲によって、インターネット上の関連サイトで公開されている線量率データを読み込み、この線量率と、この線量率データの緯度経度を含むメッシュのメッシュ識別情報と、読み込んだ日時と、３次元位置情報（公開サイトから読み込んだ緯度経度及び高度のデータ）とを対応付けて、公開線量率データベース１４に対して書き込んで記憶させる。
【００２０】
次に、図３は、統合線量率データベース１５に記憶されている統合線量率テーブルの構成例を示す図である。この図において、統合線量率データベース１５の統合線量率テーブルには、メッシュ識別情報に対応付けられて、このメッシュ識別情報の示すメッシュにおける線量率と、測定された３次元位置情報と、測定された測定日時と、メッシュの緯度経度範囲とが書き込まれて記憶されている。
【００２１】
図１に戻り、除染処理推定サーバ１２は、地図画像におけるメッシュ状の領域の緯度軽度範囲に重なる、メッシュ画像におけるメッシュに対してメッシュ識別情報を付与し、統合線量率データベース１５における統合線量率テーブルに対し、メッシュ画像における全てのメッシュのメッシュ識別情報及びこのメッシュ識別情報の示すメッシュの緯度経度範囲とを書き込んで記憶させる。
また、データベース管理サーバ１１は、メッシュの位置に含まれる緯度経度に対応する線量率データを、３次元位置情報における緯度経度により、測定線量率データベース１３の線量率テーブルから読み出す。
【００２２】
また、データベース管理サーバ１１は、読み出したメッシュ毎の線量率データを、読み出したメッシュを示すメッシュ識別情報と対応付けて、統合線量率データベース１５の統合線量率テーブルに書き込んで記憶させる。
また、データベース管理サーバ１１は、メッシュの緯度経度範囲に含まれる緯度経度に対応する線量率データを、３次元位置情報における緯度経度により、公開線量率データベース１４から読み出す。
そして、データベース管理サーバ１１は、読み出したメッシュ毎の線量率データを、読み出したメッシュを示すメッシュ識別情報及びメッシュの緯度経度範囲と対応付けて、統合線量率データベース１５の統合線量率テーブルに書き込んで記憶させる。
【００２３】
次に、図４は、地図画像上に重ねられるメッシュ画像における各メッシュの線量率の大きさを色分けして分類した例を示す図である。この図４において、色Ａ、Ｂ、Ｃ、Ｄ、Ｅの順番に線量率が大きい領域であることを示している。また、図４においては、この色の違いをハッチングの種類を変えて表現している。また、地図画像を視認できるハッチングの種類を用いて、色の替わりに線量率の大きさをハッチングの違いによって表現しても良い。例えば、表示部１７に表示する地図画像において、最も小さい線量率と最も大きい線量率とを複数の範囲、本実施形態においては５個の範囲（段階）に分け、線量率の低い方の範囲から色Ａ、Ｂ、Ｃ、Ｄ、Ｅの各々を付与し、このメッシュ画像を地図画像に重ねて表示する構成になっている。これにより、各メッシュの重なる地図画像の領域における線量率の大小が明確に視認することができる。また、メッシュ画像の各メッシュは、この図４に示すように、地図画像のメッシュ状の領域に平面視で重なる大きさに形成されている。この構成により、地図画像のメッシュ状の領域における線量率が、重なった下部の地図画像が透けて見える半透明のメッシュ画像のメッシュの色により、地図画像におけるメッシュ状の領域毎にビジュアル的に確認することができる。
【００２４】
図１に戻り、除染処理推定サーバ１２は、地図画像に重ねられるメッシュ画像における全てのメッシュに対する線量率データが得られると、これらの線量率に対応して、最も小さい線量率と最も大きい線量率とを複数の範囲、本実施形態においては５個の範囲に分ける。
そして、除染処理推定サーバ１２は、線量率に応じて各メッシュに対して色を付与して、表示部１７に表示された地図画像の緯度経度と、それぞれの緯度経度範囲を対応させ、メッシュ画像を、この地図画像に重ねて表示する。
【００２５】
また、除染処理推定サーバ１２は、各メッシュの線量率から、それぞれのメッシュにおける線源の放射能を推定するための計算を行う。例えば、除染処理推定サーバ１２は、ＳＰＭ（Ｓａｍｐｌｅｄ  Ｐａｔｔｅｒｎ  Ｍａｔｃｈｉｎｇ）法による逆推定アルゴリズムを用い、場の測定データが与えられた際、この場を生じさせているソース（源）の位置分布を逆推定する。
ここで、除染処理推定サーバ１２は、各メッシュに対して、それぞれ放射性物質の堆積による単位放射線源を仮定する。そして、除染処理推定サーバ１２は、線量率を測定した場の分布に最も近い放射線源の放射能を逐次近似的に探索していく。この結果、除染処理推定サーバ１２は、統合線量率データベース１５の統合線量率テーブルにおけるメッシュ毎の線量率と予め設定された誤差範囲、例えば±１０％以内の誤差率において近似する線量率の分布が得られた放射能分布を、推定放射能分布として出力する。
【００２６】
また、除染処理推定サーバ１２は、この放射能分布を算出する際に求められた、各メッシュの放射能が周囲のメッシュの線量率に寄与する寄与率を、予め設定した範囲のメッシュに対して放射線源データベース１６に対して書き込んで記憶させる。ここで、寄与率とは、各メッシュにおける放射能が、自身と周囲のメッシュとの各々における線量率にどの程度の割合で影響しているかを示す率である。
【００２７】
次に、図５は、放射線源データベース１６に記憶されている放射線源テーブルの構成例を示す図である。この放射線源テーブルは、メッシュ識別情報に対応付けて、線量率を測定した３次元情報と、このメッシュ識別情報の示すメッシュに対応して算出した（推定した）放射線源の放射能と、周囲のメッシュに対する寄与率とから構成されている。
図１に戻り、除染処理推定サーバ１２は、逆推定アルゴリズムの計算で得られた放射能分布による各メッシュの放射能を、それぞれ対応するメッシュのメッシュ識別情報に対応付けて、３次元位置情報と、寄与率とともに、放射線源データベース１６の放射線源テーブルに書き込んで記憶させる。
【００２８】
また、除染処理推定サーバ１２は、操作者がいずれか一つあるいは複数のメッシュを選択し、その選択したメッシュに対して除染率を入力し、この除染率による除染後の線量率の計算を指示した場合、以下の処理を行う。ここで、除染率とは、除染によって除染前に比較して残存した放射性物質の割合を示すものであり、１以下に設定される。
すなわち、除染処理推定サーバ１２は、除染した放射線分布を放射線源とし、土地の高度、周囲の建物の分布、山、あるいは丘陵の位置などを考慮した３次元空間における放射線輸送により、各メッシュにおける線量率を推定する計算を行う。ここで、除染処理推定サーバ１２は、各メッシュの放射能と、放射線源に対する除染率とにより、モンテカルロ法（ＭＣＮＰ：Monte Carlo Transport for N Particle）を用いた３次元の放射線輸送計算コード （輸送計算プログラム）により、各メッシュの線量率を算出し、算出した結果を除染後の推定線量率として出力する。
【００２９】
また、除染処理推定サーバ１２は、除染後の各メッシュにおける推定線量率を、図４に示す線量率の範囲に分類し、分類された範囲に対応する色を各メッシュに付与する。
そして、除染処理推定サーバ１２は、それぞれ色分けされたメッシュからなるメッシュ画像を、表示部１７に表示されている表示画像に重ねて表示する。
これにより、本実施形態によれば、操作者が除染率を入力することにより、この除染率に対応して変化する、地図画像における各メッシュ単位の地域の推定線量率を、容易にビジュアル的に得ることができる。
【００３０】
次に、図１、図６、図７、図８及び図９を参照して、本実施形態による３次元線量評価マッピングシステム１における除染した地域と、この地域の周囲の他の地域との推定線量率を算出する処理を説明する。図６は、本実施形態による３次元線量評価マッピングシステム１における、除染した地域とこの地域の周囲の他の地域との推定線量率を算出する動作例を示すタイミングチャートである。また、図７は、測定線量率データベース１３の線量率テーブルに記憶されている３次元位置線量率データの測定位置を表示した図である。また、図８は、公開線量率データベース１４に記憶されている３次元位置の線量率のデータの位置を表示した図である。また、図９は、図４に示すメッシュ画像と地図画像との対応関係によって、実際にシミュレーションで仮想的に求めた線量率をメッシュ画像に反映させ、地図画像に重ね合わせて表示した際の表示画像を示す図である。
【００３１】
ステップＳ１：
データベース管理サーバ１１は、データ収集装置３を介して、線量測定装置２＿１〜２＿ｎの各々から、それぞれが配置されている位置の３次元位置線量率データが供給される。
次に、データベース管理サーバ１１は、供給された３次元位置線量率データに対し、それぞれ測定識別情報を付与し、この測定識別情報毎に対応付けて、３次元位置線量率データを、測定線量率データベース１３の線量率テーブルに書き込んで記憶させる。
そして、操作者が、線量率の分布を表示させる表示エリアを、表示部１７に表示されている地図画像において選択し、線量率の分布マップを表示する処理を図示しない入力手段により入力する。
【００３２】
ステップＳ２：
次に、除染処理推定サーバ１２は、選択された範囲を新たな地図画像として表示し、予め設定された緯度経度範囲にて、正方形のメッシュ状に分割する。
そして、除染処理推定サーバ１２は、この地図画像の正方形のメッシュ状に対し、平面視で重なるメッシュから構成されるメッシュ画像を生成する。このメッシュ画像におけるメッシュの各々は、それぞれ平面視で重なる地図画像におけるメッシュ状の範囲と同一の緯度経度の範囲情報を有している。
【００３３】
ステップＳ３：
次に、除染処理推定サーバ１２は、メッシュの緯度経度により、測定線量率データベース１３の線量率テーブルから、それぞれのメッシュと重なる地図画像上のメッシュ状の領域に含まれた３次元情報に対応する線量率を読み出し、その線量率を地図画像上のメッシュ状の領域における線量率データとする。
そして、除染処理推定サーバ１２は、統合線量率データベース１５の統合線量率テーブルにメッシュのメッシュ識別情報に対応させ、書き込んで記憶させる。
【００３４】
ステップＳ４：
次に、除染処理推定サーバ１２は、統合線量率データベース１５の統合線量率テーブルにおいて、線量率データが付与されていないメッシュ識別情報を検索し、このメッシュ識別情報と、このメッシュ識別情報に対応するメッシュの緯度経度範囲とを含む線量率要求情報を、データベース管理サーバ１１に対して送信する。
これにより、データベース管理サーバ１１は、対応する緯度経度範囲の線量率データを、インターネットの専用サイトから読み込み、公開線量率データベース１４に対し、対応する緯度経度範囲のメッシュ識別情報と、この緯度経度範囲の線量率データとを組にして書き込んで記憶させる。
【００３５】
そして、除染処理推定サーバ１２は、公開線量率データベース１４から、メッシュ識別情報に対応して記憶されている線量率データを読み込み、統合線量率データベース１５の統合線量率テーブルにおける対応するメッシュ識別情報の線量率データに追加して書き込み記憶させる。
すなわち、上述した公開されている線量率により、除染処理推定サーバ１２は、線量率データが測定されないメッシュに対して、線量率データの内挿（補間）の処理を行う。
【００３６】
ステップＳ５：
次に、除染処理推定サーバ１２は、各メッシュに対して、それぞれの線量率が含まれる数値の数値範囲に対応した色を付与し、この色を付与したメッシュからなるメッシュ画像を、表示部１７に表示されている図９に示すように、地図画像に対して重ねて表示する。
これにより、操作者は、表示部１７に表示される図９の画像により、選択した表示エリア（現在、表示部１７に表示されている地図画像の領域）における各メッシュ状の領域の現在の線量率をビジュアル的に観察することができる。
【００３７】
ステップＳ６：
次に、除染処理推定サーバ１２は、ＳＰＭ法による推定アルゴリズムにより、メッシュの線量率に対応した、このメッシュと重なる地図画像のメッシュ状の領域における放射線源の放射能を順次算出する。
【００３８】
ステップＳ７：
次に、除染処理推定サーバ１２は、放射線源データベース１６の放射線源テーブルにおける、表示エリアに対応したメッシュ画像の全てのメッシュに放射能が書き込まれているか否かの判定を行う。
そして、除染処理推定サーバ１２は、メッシュ画像の全てのメッシュに放射能が放射線源テーブルに書き込まれている場合、メッシュ画像の全てのメッシュにおける放射能の算出が終了していると判定し、処理をステップＳ８へ進める。
一方、除染処理推定サーバ１２は、メッシュ画像のメッシュにおいて放射能が放射線源テーブルに書き込まれていないメッシュが存在する場合、メッシュ画像の全てのメッシュにおける放射能の算出が終了していないと判定し、処理をステップＳ６へ進め、放射能の算出がされていないメッシュを選択し、このメッシュの放射能を算出する処理を行わせる。
【００３９】
ステップＳ８及びＳ９：
次に、除染処理推定サーバ１２は、放射線源データベース１６の放射線源テーブルから、各メッシュの放射能を読み出す。
そして、除染処理推定サーバ１２は、メッシュ画像の各メッシュの放射能が自身及び周辺の他のメッシュにおける線量率の放射線の線量に寄与する寄与率を、緯度経度に高度を加えた３次元における放射線輸送計算を行って算出する。
また、除染処理推定サーバ１２は、算出した各メッシュの寄与率を、メッシュ識別情報に対応させ、放射線源データベース１６における放射線源テーブルに書き込み記憶させる。
【００４０】
ステップＳ１０：
次に、操作者が表示部１７に表示されている表示エリアのなかから、除染を行おうとする地図画像上のメッシュ状の領域を、マウスなどによりクリックすることにより、除染エリアとして選択する。
これにより、除染処理推定サーバ１２は、表示部１７において選択された地図画像上のメッシュ状の領域の緯度経度を含む緯度経度範囲のメッシュ画像のメッシュを、放射線源データベース１６から検索する。
そして、除染処理推定サーバ１２は、表示部１７において、選択した地図画像のメッシュ状の領域に対応するメッシュ画像のメッシュの枠を、メッシュ画像の他の枠に対して太く表示する。
また、除染処理推定サーバ１２は、選択した地図画像上のメッシュ状の領域の除染後の除染率の入力欄を、表示部１７に対して表示する。
【００４１】
ステップＳ１１：
次に、操作者が表示部１７に表示された入力欄に対して、除染率をキーボード等により入力する。
これにより、除染処理推定サーバ１２は、入力欄に入力された除染率を読み込み、自身の内部記憶部に除染率が設定されたメッシュのメッシュ識別情報とともに書き込んで、除染率設定テーブルとして記憶させる。
【００４２】
ステップＳ１２：
次に、除染処理推定サーバ１２は、メッシュ画像の全てのメッシュの放射能を、メッシュのメッシュ識別情報とともに、放射線源データベース１６の放射線源テーブルから読み込み、自身の内部記憶部に対して書き込んで記憶させる。
そして、除染処理推定サーバ１２は、読み出したメッシュから除染率の設定されたメッシュを、除染率設定テーブルからメッシュ識別情報により抽出し、この抽出された除染率の設定されたメッシュの放射能の寄与率に対して、メッシュに設定された除染率を乗算し、除染後の推定寄与率を算出する。ここで、除染処理推定サーバ１２は、除染率設定テーブルにおいて、推定寄与率の算出が終了したメッシュのメッシュ識別情報に対しフラグを設定する。
【００４３】
ステップＳ１３：
次に、除染処理推定サーバ１２は、自身の内部記憶部の除染率設定テーブルにおいて、フラグが設定されていないメッシュ識別情報の有無の判定を行う。
そして、除染処理推定サーバ１２は、内部記憶部の除染率設定テーブルの全てのメッシュ識別情報に対してフラグが設定されている場合、除染率が設定されたメッシュの推定寄与率の算出が終了していると判定し、処理をステップＳ１４へ進める。
一方、除染処理推定サーバ１２は、内部記憶部の除染率設定テーブルのメッシュ識別情報にフラグが設定されていないメッシュ識別情報が存在する場合、除染率が設定されたメッシュの推定寄与率の算出が終了していない判定し、処理をステップＳ１４へ進め、フラグが設定されていなメッシュを選択し、このメッシュの推定寄与率の算出処理を行わせる。
【００４４】
ステップＳ１４：
次に、除染処理推定サーバ１２は、各メッシュの放射能と、各メッシュの寄与率（除染されない領域のメッシュに対して寄与率、一方、除染された領域のメッシュに対して推定寄与率）とを用い、地図画像上の各メッシュ状の領域における線量率を算出する。
【００４５】
ステップＳ１５：
次に、除染処理推定サーバ１２は、除染後の線量率を算出することで推定した後、この推定した地図画像上のメッシュ状の領域の線量率が含まれる数値の数値範囲に対応した色を、このメッシュ状の領域に重なるメッシュに付与し、この色を付与したメッシュからなるメッシュ画像を、表示部１７に表示されている地図画像に対して重ねて表示する。このときの、色と数値範囲との関係は、ステップＳ５における表示処理の際の図９のメッシュ画像と同様である。
これにより、操作者は、選択した表示エリア（現在、表示部１７に表示されている地図画像の領域）における各メッシュ状の領域の除染後の線量率をビジュアル的に観察することができる。
【００４６】
また、操作者が除染前の状態との比較処理を、３次元線量評価マッピングシステム１に対して指示する制御を行う。
この比較処理が指示されると、除染処理推定サーバ１２は、除染前の線量率に対応した色分布（メッシュの色分布）のメッシュ画像が重ねられた地図画像と、除染後の線量率に対応した色分布のメッシュ画像が重ねられた地図画像との各々を、表示部１７の画面における左右領域に列べて表示する。
これにより、操作者は、選択した表示エリア（現在、表示部１７に表示されている地図画像の領域）における除染前の線量率の分布と、除染後の線量率の分布とをビジュアル的に比較することができ、除染処理の効果を明確に把握することができる。
【００４７】
また、操作者が測定線量率データベース１３の線量率テーブルに記憶されている３次元位置線量率データの測定位置の表示処理を、３次元線量評価マッピングシステム１に対して指示する制御を行った場合、データベース管理サーバ１１は、操作者の指定した地図画像における緯度経度範囲内の３次元位置線量率データの測定位置を、表示部１７に対して図７に示すように表示する。
また、操作者が公開線量率データベース１４に記憶されている３次元位置の線量率のデータの測定位置の表示処理を、３次元線量評価マッピングシステム１に対して指示する制御を行った場合、データベース管理サーバ１１は、操作者の指定した地図画像における緯度経度範囲内の３次元位置の線量率のデータの測定位置を、表示部１７に対して図８に示すように表示する。
また、本実施形態においては、メッシュの形状を正方形として説明したが、正方形に限らず直方形で構成するようにしても良い。
【００４８】
また、図１におけるデータベース管理サーバ１１と、除染処理推定サーバ１２との各々の機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより除染後の推定線量率を算出して推定する処理を行ってもよい。なお、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
【００４９】
また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。
【００５０】
以上、この発明の実施形態を図面を参照して詳述してきたが、具体的な構成はこの実施形態に限られるものではなく、この発明の要旨を逸脱しない範囲の設計等も含まれる。
【符号の説明】
【００５１】
１…３次元線量評価マッピングシステム
２＿１，２＿ｎ…線量率測定装置
３…データ収集装置
１１…データベース管理サーバ
１２…除染処理推定サーバ
１３…測定線量率データベース
１４…公開線量率データベース
１５…統合線量率データベース
１６…放射線源データベース
１７…表示部
２１…線量率測定器
２２…ＧＰＳ・高度測定器
【図１】
0005839280.tif 000002
【図２】
0005839280.tif 000003
【図３】
0005839280.tif 000004
【図４】
0005839280.tif 000005
【図５】
0005839280.tif 000006
【図６】
0005839280.tif 000007
【図７】
0005839280.tif 000008
【図８】
0005839280.tif 000009
【図９】
0005839280.tif 000010