特許 6018032 レーダ画像判読支援装置及びレーダ画像判読支援方法

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

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

(11)【特許番号】6018032

(24)【登録日】2016年10月7日

(45)【発行日】2016年11月2日

(54)【発明の名称】レーダ画像判読支援装置及びレーダ画像判読支援方法

(51)【国際特許分類】

   G06T 1/00 20060101AFI20161020BHJP

   G01S 13/90 20060101ALI20161020BHJP

【ＦＩ】

   G06T1/00 285

   G01S13/90

【請求項の数】9

【全頁数】20

(21)【出願番号】特願2013-214629(P2013-214629)

(22)【出願日】2013年10月15日

(65)【公開番号】特開2015-79296(P2015-79296A)

(43)【公開日】2015年4月23日

【審査請求日】2015年8月31日

(73)【特許権者】

【識別番号】000135771

【氏名又は名称】株式会社パスコ

(74)【代理人】

【識別番号】110000925

【氏名又は名称】特許業務法人信友国際特許事務所

(72)【発明者】

【氏名】下村 博之

(72)【発明者】

【氏名】鵜殿 俊昭

(72)【発明者】

【氏名】吉川 和男

(72)【発明者】

【氏名】三五 大輔

(72)【発明者】

【氏名】木村 詩織

【審査官】 村松 貴士

(56)【参考文献】

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

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

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

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

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

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

Ｇ０６Ｔ １／００

Ｇ０６Ｔ １１／６０

Ｇ０１Ｓ ７／００ − ７／４２

Ｇ０１Ｓ １３／００ − １３／９５

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

【請求項1】

幾何歪みなしの地図座標の画像を幾何歪みありの地図座標の画像に変換する歪みあり座標変換テーブルに基づいて、幾何歪みのない地図座標の参照画像を、幾何歪みありの地図座標の比較画像に変換する歪みあり座標変換部と、
飛翔体に搭載されたレーダ装置で観測された幾何歪みありのスラントレンジ座標の第１のレーダ画像を変換して得られる幾何歪みありの地図座標の第２のレーダ画像と、前記歪みあり座標変換部で得られた前記比較画像を出力する出力部と、
を備えるレーダ画像判読支援装置。

【請求項2】

幾何歪みありのスラントレンジ座標の前記第１のレーダ画像を、幾何歪みありの地図座標の前記第２のレーダ画像に変換する第１の地図投影部と、を更に備える
請求項１に記載のレーダ画像判読支援装置。

【請求項3】

スラントレンジ座標の幾何歪みありの前記第１のレーダ画像を、幾何歪みなしの地図座標の第３のレーダ画像に変換する第２の地図投影部と、を更に備え、
前記第１の地図投影部は、前記第１のレーダ画像を前記第２のレーダ画像に変換する過程で、前記第１のレーダ画像を前記第２のレーダ画像に変換する第１の変換テーブル及びその逆の変換を行う第１の逆変換テーブルを生成し、
前記第２の地図投影部は、前記第１のレーダ画像を前記第３のレーダ画像に変換する過程で、前記第１のレーダ画像を前記第３のレーダ画像に変換する第２の変換テーブル及びその逆の変換を行う第２の逆変換テーブルを生成し、
前記歪みあり座標変換テーブルは、前記第２の逆変換テーブルと前記第１の変換テーブルから構成される
請求項２に記載のレーダ画像判読支援装置。

【請求項4】

幾何歪みありの地図座標の画像を幾何歪みなしの地図座標の画像に変換する歪みなし座標変換テーブルに基づいて、前記第２のレーダ画像のユーザから指定された領域を、幾何歪みなしの地図座標に変換する歪みなし座標変換部と、を更に備える
請求項３に記載のレーダ画像判読支援装置。

【請求項5】

前記歪みなし座標変換部は、前記第２のレーダ画像又は前記比較画像の前記指定された領域を幾何歪みなしの地図座標に変換した後、地理情報として出力する
請求項４に記載のレーダ画像判読支援装置。

【請求項6】

前記歪みなし座標変換テーブルは、前記第１の逆変換テーブルと前記第２の変換テーブルから構成される
請求項４又は５に記載のレーダ画像判読支援装置。

【請求項7】

幾何歪みなしの地図座標の画像を、幾何歪みありのスラントレンジ座標の画像に変換する歪みあり座標変換テーブルに基づいて、幾何歪みのない地図座標の参照画像を、幾何歪みありのスラントレンジ座標の比較画像に変換する歪みあり座標変換部と、
飛翔体に搭載されたレーダ装置で観測された幾何歪みありのスラントレンジ座標のレーダ画像から、前記歪みあり座標変換部で得られた前記比較画像に対する変化域を抽出する抽出部と、
を備えるレーダ画像判読支援装置。

【請求項8】

レーダ画像判読支援装置の歪みあり座標変換部が、幾何歪みなしの地図座標の画像を幾何歪みありの地図座標の画像に変換する歪みあり座標変換テーブルに基づいて、幾何歪みのない地図座標の参照画像を、幾何歪みありの地図座標の比較画像に変換し、
出力部が、飛翔体に搭載されたレーダ装置で観測された幾何歪みありのスラントレンジ座標の第１のレーダ画像を変換して得られる幾何歪みありの地図座標のレーダ画像と、前記歪みあり座標変換部で得られた前記比較画像を出力する
レーダ画像判読支援方法。

【請求項9】

レーダ画像判読支援装置の歪みあり座標変換部が、幾何歪みなしの地図座標の画像を幾何歪みありのスラントレンジ座標の画像に変換する歪みあり座標変換テーブルに基づいて、幾何歪みのない地図座標の参照画像を、幾何歪みありのスラントレンジ座標の比較画像に変換し、
抽出部が、飛翔体に搭載されたレーダ装置で観測された幾何歪みありのスラントレンジ座標のレーダ画像から、前記歪みあり座標変換部で得られた前記比較画像に対する変化域を抽出する
レーダ画像判読支援方法。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、人工衛星等の飛翔体に搭載され、地表を広範囲にわたり観測するレーダ装置により得られるレーダ画像を判読するレーダ画像判読支援装置及びレーダ画像判読支援方法に関する。

【背景技術】

【0002】

従来、地表に照射した電波（マイクロ波パルス）の反射波を測定する能動型センサとしてレーダ装置がある。レーダ装置の一例として、例えば合成開口レーダ（SAR：Synthetic Aperture Radar）がある。合成開口レーダは、マイクロ波の特性を利用して、昼夜、天候に関係なく、広域を面的に観測することができる。

【0003】

合成開口レーダを用いて観測されたレーダ画像（ＳＡＲ画像）には、レイオーバやフォアショートニングなどの幾何歪みが含まれている。ＳＡＲ画像を用いて斜面崩壊地等を抽出する際に、地図投影された幾何歪みが補正されていないＳＡＲ画像を使用すると、判読性が良いことが知られている。地図投影された幾何歪みが補正されていないＳＡＲ画像として、例えばＴｅｒｒａＳＡＲ-Ｘ衛星の観測データから生成されるＳＡＲ画像を地図投影したＧＥＣ（Geocoded Ellipsoid Corrected）画像がある。

【0004】

ＧＥＣ画像は、ＳＡＲ画像特有の幾何歪みの補正（以下、「オルソ補正」と記す）が行われていないため、既存のＧＩＳ（Geographic Information System）情報や航空写真（オルソ補正済み画像）とのマッチング（位置の整合）に問題がある。

【0005】

一方、オルソ補正済み画像（以下、「オルソ画像」と略称する）は、幾何歪みが補正されているが、そのために判読性が著しく落ちるという問題がある。

【0006】

図１８ＡはＧＥＣ画像（ＳＡＲ画像）を示し、図１８Ｂはオルソ画像（ＳＡＲ画像）を示す。また図１９ＡはＧＥＣ画像（ＳＡＲ画像）を示し、図１９Ｂはオルソ画像（ＳＡＲ画像）を示す。例えば図１８Ａに示すＧＥＣ画像では、谷底と河川（太線）の位置が一致しないが、図１８Ｂに示すオルソ画像では、谷底と河川（太線）の位置が一致する。しかし、図１９ＡのＧＥＣ画像では、土砂崩れ発生地（破線部１０１）を判読可能であるが、図１９Ｂのオルソ画像では、土砂崩れ発生地に対応する地点（破線部１０２）は白みがかっており、判読性が悪化している。

【0007】

ＳＡＲ画像の判読性を向上させるために、ＳＡＲ画像の解析に際し、観測対象の地域が含まれる光学画像を参照画像として用いてＳＡＲ画像と光学画像を比較することが行われている。

【0008】

例えば、ＳＡＲ／光学観測衛星に搭載されたＳＡＲと光学センサとにより地上の同じ地域を観測し、ＳＡＲ／光学伝送装置によりＳＡＲの観測データと光学センサの観測データとを重畳させた重畳データを生成し、重畳データをダウンリンクすることが記載されている（例えば、特許文献１参照。）。地上センタでは、ＳＡＲ／光学再生装置により重畳データからＳＡＲ画像と光学画像とを再生する。

【先行技術文献】

【特許文献】

【0009】

【特許文献1】特開２０１０−２３６９７０号公報

【発明の概要】

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

【0010】

地図に投影された幾何歪みが補正されたＳＡＲ画像と光学画像を比較する場合、両画像は地図投影された画像であるためマッチングに問題はないが、斜面崩壊地等の災害発生地を抽出する際の判読性が落ちてしまう。

【0011】

上記の状況から、本発明は、幾何歪みありのレーダ画像上の地表の変化域を容易に判読できるようにすることを目的とする。

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

【0012】

本発明の第１の態様のレーダ画像判読支援装置は、歪みあり座標変換部と出力部を備える。
歪みあり座標変換部は、幾何歪みなしの地図座標の画像を幾何歪みありの地図座標の画像に変換する歪みあり座標変換テーブルに基づいて、幾何歪みのない地図座標の参照画像を、幾何歪みありの地図座標の比較画像に変換する。
出力部は、飛翔体に搭載されたレーダ装置で観測された幾何歪みありのスラントレンジ座標の第１のレーダ画像を変換して得られる幾何歪みありの地図座標の第２のレーダ画像と、歪みあり座標変換部で得られた比較画像を出力する。

【0013】

本発明の第２の態様のレーダ画像判読支援装置は、歪みあり座標変換部と抽出部を備える。
歪みあり座標変換部は、幾何歪みなしの地図座標の画像を、幾何歪みありのスラントレンジ座標の画像に変換する歪みあり座標変換テーブルに基づいて、幾何歪みのない地図座標の参照画像を、幾何歪みありのスラントレンジ座標の比較画像に変換する。
抽出部は、飛翔体に搭載されたレーダ装置で観測された幾何歪みありのスラントレンジ座標のレーダ画像から、歪みあり座標変換部で得られた比較画像に対する変化域を抽出する。

【0014】

本発明の第１の態様及び第２の態様は、幾何歪みのない地図座標の参照画像が、幾何歪みあり且つ判読対象の幾何歪みありのレーダ画像と同じ座標の比較画像に変換される。この比較画像を用いて、判読対象の幾何歪みありのレーダ画像の判読が行われる。

【発明の効果】

【0015】

本発明の少なくとも一つの態様によれば、幾何歪みありのレーダ画像上の地表の変化域を容易に判読できる。

【図面の簡単な説明】

【0016】

【図1】本発明の第１の実施の形態に係るデータ解析システムの概要を示す構成図である。

【図2】本発明の第１の実施の形態に係るデータ解析の概要を示す説明図である。

【図3】図３Ａは衛星ＳＡＲの飛行方向と観測方向の説明図、図３Ｂはスラントレンジ座標とグランドレンジ座標を示す説明図である。

【図4】本発明の第１の実施の形態に係るデータ解析装置のハードウェアの構成例を示すブロック図である。

【図5】本発明の第１の実施の形態に係るデータ解析装置の機能ブロック図である。

【図6】再生画像（ＳＡＲ画像）から地図投影された幾何歪みありのレーダ画像又はオルソ画像への変換を示す説明図である。

【図7】図７Ａはスラントレンジ座標のＳＡＲ画像を示し、図７Ｂはグランドレンジ座標のＳＡＲ画像である。

【図8】図８Ａは南行軌道を通る衛星ＳＡＲから取得したＳＡＲ画像を地図投影したときのＳＡＲ画像を示し、図８Ｂは北行軌道を通る衛星ＳＡＲから取得したＳＡＲ画像を地図投影したときのＳＡＲ画像を示している。

【図9】ＳＡＲ画像からオルソ画像への変換を示す説明図である。

【図10】本発明の第１の実施の形態に係るデータ解析装置のＣＰＵによるデータ解析の手順を示すフローチャートである。

【図11】本発明の第１の実施の形態に係るデータ解析の手順を示す説明図である。

【図12】本発明の第２の実施の形態に係るデータ解析の概要を示す説明図である

【図13】本発明の第２の実施の形態に係るデータ解析装置の機能ブロック図である。

【図14】本発明の第２の実施の形態に係るデータ解析装置のＣＰＵによるデータ解析の手順を示すフローチャートである。

【図15】本発明の第２の実施の形態に係るデータ解析の手順を示す説明図である。

【図16】図１６ＡはＳＡＲ画像（オルソ画像）を示し、図１６Ｂは航空写真（オルソ画像）を示す。

【図17】ＧＥＣ画像と比較画像が表示された表示画面の一例である。

【図18】図１８ＡはＧＥＣ画像を示し、図１８Ｂはオルソ画像を示す。

【図19】図１９ＡはＧＥＣ画像を示し、図１９Ｂはオルソ画像を示す。

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

【0017】

以下、本発明を実施するための形態の例について、添付図面を参照しながら説明する。なお、各図において共通の構成要素には、同一の符号を付して重複する説明を省略する。

【0018】

＜１．第１の実施の形態＞
第１の実施の形態は、地図投影された幾何歪みありのレーダ画像を、幾何歪みのない地図座標の参照画像（例えば、光学画像又は地理情報）を地図投影された幾何歪みありの画像に変換した比較画像と比較することで、地図投影された幾何歪みありのレーダ画像を判読する。以下、レーダ画像の一例としてＳＡＲ画像を用いて説明する

【0019】

［データ解析システムの全体構成］
図１は、本発明の第１の実施の形態に係るデータ解析システムの概要を示す構成図である。
図１に示すデータ解析システムは、人工衛星に搭載されたレーダ装置の一例である合成開口レーダ（ＳＡＲ）を用いて地表３のＳＡＲ画像データを取得し、取得したＳＡＲ画像データをデータ解析センタ２で解析する。以下、人工衛星に搭載した合成開口レーダを「衛星ＳＡＲ」と記す。ここでは、レーダ装置を搭載する飛翔体として人工衛星を例示しているが、航空機でもよい。

【0020】

衛星ＳＡＲ１は、所定の軌道を周回しながら定期的に地表３へマイクロ波パルスを照射し、その反射波を受信して得たデータをＳＡＲ画像データとしてデータ解析センタ２へ送信する。

【0021】

データ解析センタ２は、観測計画を行い、その観測計画に基づいた観測指示を含む無線信号を、地上アンテナ４を介して衛星ＳＡＲ１に送信する。また、衛星ＳＡＲ１が取得したＳＡＲ画像データを、地上アンテナ４を介して受信する。これをデータ解析装置５（レーダ画像判読支援装置の一例）で合成開口処理することにより再生画像を生成する。

【0022】

データ解析装置５による一連の処理は、ハードウェアにより実行することもできるし、全部又は一部をソフトウェアにより実行させることもできる。一連の処理の全部又は一部をソフトウェアにより実行させる場合には、そのソフトウェアを構成するプログラムが専用のハードウェアに組み込まれているコンピュータ、または、各種の機能を実行するためのプログラムをインストールしたコンピュータにより、実行可能である。例えば汎用のパーソナルコンピュータなどに所望のソフトウェアを構成するプログラムをインストールして実行させてもよい。

【0023】

［データ解析の概要］
図２は、データ解析装置５により実行されるデータ解析の概要を示す説明図である。
まずデータ解析装置５は、衛星ＳＡＲ１又はその他から、オルソ補正されていないスラントレンジ座標のレーダ画像１１（ＳＡＲ画像）を取得する。このレーダ画像１１を地図座標に変換（地図投影）し、幾何歪みありの地図座標のレーダ画像１２を得る。

【0024】

一方、データ解析装置５は、オルソ補正された幾何歪みなしの地図座標の既存の情報１３を、スラントレンジ座標に変換（スラントレンジ化）した後、地図投影し、幾何歪みありの地図座標の既存の情報１４を得る。既存の情報１３として、例えばＧＩＳ情報等の地理情報、航空写真等の光学画像が用いられる。

【0025】

データ解析装置５は、レーダ画像１２と既存の情報１４を不図示の表示装置に出力する。ユーザは、表示画面に表示されたレーダ画像１２を既存の情報１４と比較して、レーダ画像１２上の地表３の変化域を判読する。そして、ユーザが入力部を操作してデータ解析装置５に判読結果１５を入力すると、データ解析装置５は、幾何歪みありの地図座標の判読結果１５をスラントレンジ化及びオルソ補正し、幾何歪みなしの地図座標の判読結果１６を得る。

【0026】

すなわち本実施の形態では、スラントレンジ座標のレーダ画像１１を地図投影したレーダ画像１２をユーザが判読する際に、既存の情報１３をスラントレンジ化及び地図投影して得られた既存の情報１４を参照して判読性を向上させている。

【0027】

図３において、図３Ａは衛星ＳＡＲ１の飛行方向と観測方向の説明図、図３Ｂはスラントレンジ座標とグランドレンジ座標を示す説明図である。図３Ａ，図３Ｂの原点Ｏは、人工衛星の鉛直直下にある。
人工衛星の飛行方向（進行方向）に対して直角方向をレンジ方向という。このレンジ方向にはスラントレンジ（視線方向）とグランドレンジがある。
スラントレンジは、人工衛星の飛行方向から直角に対象物に向かう方向である。スラントレンジ座標は、人工衛星（衛星ＳＡＲ１）と対象物との距離すなわち視線距離を表している。
グランドレンジは、スラントレンジ座標の地表３（楕円体）への投影距離を表している。図３Ｂでは、スラントレンジ座標の点Ａがグランドレンジ座標の点Ａ´に投影されている。

【0028】

［データ解析装置の構成］
図４は、データ解析装置５のハードウェアの構成例を示すブロック図である。この例では、データ解析装置５として汎用のコンピュータを用いている。

【0029】

データ解析装置５のＣＰＵ（Central Processing Unit）２１は、ＲＯＭ（Read Only Memory）２２、または記録部２８に記録されているプログラムに従って、上記一連の処理の他、各種の処理を実行する。ＲＡＭ（Random Access Memory）２３には、ＣＰＵ２１が実行するプログラムやデータなどが適宜記憶される。これらのＣＰＵ２１、ＲＯＭ２２、およびＲＡＭ２３は、バス２４により相互に接続されている。

【0030】

ＣＰＵ２１にはまた、バス２４を介して入出力インタフェース２５が接続されている。入出力インタフェース２５には、キーボード、マウス、マイクロホンなどよりなる入力部２６、ディスプレイ、スピーカなどよりなる出力部２７が接続されている。ＣＰＵ２１は、入力部２６から入力される指令に対応して各種の処理を実行する。そして、ＣＰＵ２１は、処理の結果を出力部２７に出力する。

【0031】

入出力インタフェース２５に接続されている記録部２８は、例えばハードディスクからなり、ＣＰＵ２１が実行するプログラムや各種のデータを記録する。

【0032】

通信部２９は、インターネットやローカルエリアネットワークなどのネットワークを介して外部の装置と通信する。この通信部２９を介してプログラムを取得し、記録部２８に記録してもよい。

【0033】

入出力インタフェース２５に接続されているドライブ３０は、磁気ディスク（フレキシブルディスクを含む）、光ディスク（ＣＤ−ＲＯＭ（Compact Disc - Read Only Memory）、ＤＶＤ（Digital Versatile Disc），光磁気ディスク等を含む）、あるいは半導体メモリなどのリムーバブルメディア３１が装着されたとき、それらを駆動し、そこに記録されているプログラムやデータなどを取得する。取得されたプログラムやデータは、必要に応じて記録部２８に転送され、記録される。

【0034】

コンピュータにインストールされ、コンピュータによって実行可能な状態とされるプログラムを格納するプログラム記録媒体は、磁気ディスク、光ディスク、もしくは半導体メモリなどよりなるパッケージメディアであるリムーバブルメディア３１、または、プログラムが一時的もしくは永続的に格納されるＲＯＭ２２や、記録部２８を構成するハードディスクなどにより構成される。プログラム記録媒体へのプログラムの格納は、必要に応じてルータ、モデムなどのインタフェースである通信部２９を介して、ローカルエリアネットワーク、インターネット、デジタル衛星放送といった、有線または無線の通信媒体を利用して行ってもよい。

【0035】

図５は、第１の実施の形態に係るデータ解析装置５の機能ブロック図である。
データ解析装置５のＣＰＵ２１は、第１の地図投影部２１１と、第２の地図投影部２１２と、歪みあり座標変換部２１３と、歪みなし座標変換部２１４を有している。

【0036】

第１の地図投影部２１１は、通信部２９を介して衛星ＳＡＲ１からレーダ画像データ（ＳＡＲ画像データ）を取得する。そして、このレーダ画像データから得られるスラントレンジ座標のレーダ画像を地図投影する。すなわち、第１の地図投影部２１１は、幾何歪みありのスラントレンジ座標のレーダ画像（第１のレーダ画像）を、幾何歪みありの地図座標のレーダ画像（第２のレーダ画像）に変換する。第１の地図投影部２１１は、幾何歪みありの地図座標のレーダ画像を、記録部２８（図４）に保存し、又は出力部２７に送る。

【0037】

また第１の地図投影部２１１は、スラントレンジ座標のレーダ画像を地図投影する過程で、幾何歪みありのスラントレンジ座標のレーダ画像を、幾何歪みありの地図座標のレーダ画像に変換する第１の変換テーブル２１１ａと、その逆の変換を行う第１の逆変換テーブル２１１ｂを生成する。第１の地図投影部２１１は、生成した第１の変換テーブル２１１ａを歪みあり座標変換部２１３へ供給し、第１の逆変換テーブル２１１ｂを歪みなし座標変換部２１４へ送る。なお、１つのテーブルで、変換および逆変換を実施しても構わない。

【0038】

図６は、再生画像（ＳＡＲ画像）から地図投影された幾何歪みありのレーダ画像（ＧＥＣ画像、右方向）又はオルソ画像（左方向）への変換を示す説明図である。
図６に示すように、第１の地図投影部２１１は、衛星ＳＡＲ１から受信したレーダ画像データ（ＳＡＲ画像データ）に対応する再生画像（ＳＡＲ画像４１）のスラントレンジ座標を、第１の変換テーブル２１１ａに基づいて、グランドレンジ座標に変換した後、さらに地図座標に変換する。これにより、幾何歪みありのスラントレンジ座標のＳＡＲ画像４１が、幾何歪みありの地図座標のレーダ画像（ＧＥＣ画像４２）に変換される。ＧＥＣ画像４２は、地図投影されているが、幾何歪みについて未補正のため、光学画像や地理情報とは位置関係が一致しない。

【0039】

図７Ａはスラントレンジ座標のＳＡＲ画像を示し、図７Ｂはグランドレンジ座標のＳＡＲ画像を示す。
スラントレンジ座標のＳＡＲ画像５５（図７Ａ）の座標Ａ、座標Ｂ、座標Ｃは、グランドレンジ化により、グランドレンジ座標のＳＡＲ画像５６（図７Ｂ）の座標Ａ´、座標Ｂ´、座標Ｃ´に変換されている。

【0040】

図８Ａは南行軌道を通る衛星ＳＡＲ１から取得したＳＡＲ画像を地図投影したときのＳＡＲ画像を示し、図８Ｂは北行軌道を通る衛星ＳＡＲ１から取得したＳＡＲ画像を地図投影したときのＳＡＲ画像を示している。

【0041】

人工衛星は基準面（北極から南極を通る面）に対し軌道傾斜角を持っているため、人工衛星の飛行方向に応じて衛星ＳＡＲ１の観測領域は経線に対して傾斜する。人工衛星が北から南（南行軌道）へ飛行した場合には、グランドレンジ座標のＳＡＲ画像５６を地図座標へ変換したＳＡＲ画像５７ａ（オルソ画像又はＧＥＣ画像）の観測領域は、飛行時間とともに西寄りになる（図８Ａ）。一方、人工衛星が南から北（北行軌道）へ飛行した場合には、グランドレンジ座標のＳＡＲ画像５６を地図座標へ変換したＳＡＲ画像５７ｂの観測領域は、飛行時間とともに西寄りになる（図８Ｂ）。図６に例示したＧＥＣ画像４２及びオルソ画像４３は、北行軌道を通った衛星ＳＡＲ１から取得したＳＡＲ画像４１を地図投影した画像である。

【0042】

人工衛星（衛星ＳＡＲ１）の軌道情報（飛行高度、飛行方向、飛行速度、マイクロ波の照射角度及び飛行時刻等）がわかると、衛星ＳＡＲ１と地表３の観測領域との視線距離（スラントレンジ座標）が求められる。このスラントレンジ座標に基づいて地表３の観測領域のグランドレンジ座標を求める。そして、地表３の観測領域のグランドレンジ座標から地表３の観測領域の緯度及び経度が求まる。最後に、グランドレンジ座標上の地表３の観測領域の頂点の緯度及び経度を、地図座標上の緯度及び経度に変換し、観測領域の全域をグランドレンジ座標から地図座標に変換する。第１の地図投影部２１１は、このスラントレンジ座標からオルソ補正を伴わない地図座標への変換処理を反映した第１の変換テーブル２１１ａと、その逆変換を行う第１の逆変換テーブル２１１ｂを生成する。生成された第１の変換テーブル２１１ａと、第１の逆変換テーブル２１１ｂは、記録部２８（図４）に保存される。

【0043】

図５の説明に戻る。第２の地図投影部２１２は、通信部２９を介して衛星ＳＡＲ１から受信したレーダ画像データ（ＳＡＲ画像データ）に対応するレーダ画像を、オルソ補正及び地図投影する。すなわち、第２の地図投影部２１２は、幾何歪みありのスラントレンジ座標のレーダ画像（第１のレーダ画像）を、幾何歪みなしの地図座標のレーダ画像（第３のレーダ画像）に変換する。第２の地図投影部２１２は、幾何歪みなしの地図座標のレーダ画像を、記録部２８（図４）に保存し、又は出力部２７に送る。

【0044】

また第２の地図投影部２１２は、スラントレンジ座標のレーダ画像をオルソ補正及び地図投影する過程で、幾何歪みありのスラントレンジ座標のレーダ画像を、幾何歪みなしの地図座標のレーダ画像に変換する第２の変換テーブル２１２ａと、その逆の変換を行う第２の逆変換テーブル２１２ｂを生成する。第２の地図投影部２１２は、生成した第２の変換テーブル２１２ａを歪みなし座標変換部２１４へ供給し、第２の逆変換テーブル２１２ｂを歪みあり座標変換部２１３へ供給する。

【0045】

図６に示すように、第２の地図投影部２１２は、衛星ＳＡＲ１から受信したレーダ画像データ（ＳＡＲ画像データ）に対応する再生画像（ＳＡＲ画像４１）のスラントレンジ座標を、第２の変換テーブル２１２ａに基づいて、オルソ補正を行った後、地図座標に変換する。これにより、幾何歪みありのスラントレンジ座標のＳＡＲ画像４１が、幾何歪みなしの地図座標のレーダ画像（オルソ画像４３）に変換される。また第２の地図投影部２１２は、第２の逆変換テーブル２１２ｂに基づいて、オルソ画像４３を再生画像であるＳＡＲ画像４１に変換する。

【0046】

図９は、ＳＡＲ画像４１からオルソ画像への変換を示す説明図である。この例では、オルソ画像の一例として数値標高モデル（ＤＥＭ：Digital Elevation Model）を用いて説明している。図９では、数値標高モデルをＤＥＭと記している。

【0047】

人工衛星（衛星ＳＡＲ１）の軌道情報と数値標高モデル６１から得られる仮想的な地表を用いて、衛星ＳＡＲ１で観測されたＳＡＲ画像４１の観測領域に対応するスラントレンジ座標のシミュレーション画像６２を生成する。数値標高モデル６１（オルソ画像）の地図座標とシミュレーション画像６２のスラントレンジ座標とを対応づける（処理Ｐ１）。

【0048】

次に、生成したシミュレーション画像６２の各部のスラントレンジ座標とＳＡＲ画像４１の各部のスラントレンジ座標とを対応づける（処理Ｐ２）。

【0049】

処理Ｐ１により数値標高モデル６１（オルソ画像）とシミュレーション画像６２が対応づけられ、処理Ｐ２によりシミュレーション画像６２とＳＡＲ画像４１がそれぞれ対応付けられる。それにより、ＳＡＲ画像４１のスラントレンジ座標と数値標高モデル６１（オルソ補正済み）の地図座標とが対応づけられる。すなわち、ＳＡＲ画像４１のスラントレンジ座標とオルソ画像４３（図６）の地図座標とが対応づけられる（処理Ｐ３）。第２の地図投影部２１２は、このスラントレンジ座標からオルソ補正済みの地図座標への変換処理を反映した第２の変換テーブル２１２ａと、その逆変換を行う第２の逆変換テーブル２１２ｂを生成する。生成された第２の変換テーブル２１２ａと第２の逆変換テーブル２１２ｂは、記録部２８（図４）に保存される。

【0050】

図５の説明に戻る。歪みあり座標変換部２１３は、幾何歪みなしの地図座標の画像を幾何歪みありの地図座標の画像に変換する歪みあり座標変換テーブルに基づいて、幾何歪みのない地図座標の参照画像２１５（図２の既存の情報１３に相当）を、幾何歪みありの地図座標の比較画像に変換する。参照画像２１５として、例えば光学画像又は地理情報が用いられる。すなわち、光学画像又は地理情報のラスタ情報又はベクタ情報が、幾何歪みありの地図座標に変換される。歪みあり座標変換テーブルは、例えば記録部２８（図４）に保存されており、第１の変換テーブル２１１ａと第２の逆変換テーブル２１２ｂの要素を有する。この例では、歪みあり座標変換部２１３は、第１の地図投影部２１１で生成された第１の変換テーブル２１１ａと、第２の地図投影部２１２で生成された第２の逆変換テーブル２１２ｂから歪みあり座標変換テーブルを生成している。歪みあり座標変換部２１３は、当該比較画像を、記録部２８に保存し、又は出力部２７に送る。

【0051】

ユーザは、例えば出力部２７としての表示装置の表示画面に表示された幾何歪みありの地図座標のレーダ画像を比較画像と比較し、幾何歪みありの地図座標のレーダ画像（ＧＥＣ画像４２）を判読する。そして、ユーザは、入力部２６を操作して、幾何歪みありの地図座標のレーダ画像上の災害発生地等（判読結果）を歪みなし座標変換部２１４へ入力する。

【0052】

歪みなし座標変換部２１４は、幾何歪みありの地図座標の画像を幾何歪みなしの地図座標の画像に変換する歪みなし座標変換テーブルに基づいて、入力部２６より入力された判読結果に対応する幾何歪みありの地図座標の画像を、幾何歪みなしの地図座標の画像に変換する。歪みなし座標変換テーブルは、例えば記録部２８（図４）に保存されており、第１の逆変換テーブル２１１ｂと第２の変換テーブル２１２ａの要素を有する。この例では、歪みなし座標変換部２１４は、第１の地図投影部２１１で生成された第１の逆変換テーブル２１１ｂと、第２の地図投影部２１２で生成された第２の変換テーブル２１２ａから歪みあり座標変換テーブルを生成している。歪みなし座標変換部２１４は、上記判読結果を、記録部２８に保存し、又は出力部２７に送る。

【0053】

なお、参照画像２１５を、記録部２８（図４）又はリムーバブルディスク３１から取得してもよいし、通信部２９を介してデータ解析装置５の外部から取得してもよい。

【0054】

［レーダ画像の判読を支援する手順］
次に、レーダ画像の判読を支援する手順を、図１０及び図１１を参照しながら具体的に説明する。
図１０は、データ解析装置５のＣＰＵ２１によるデータ解析の手順を示すフローチャートである。図１１は、データ解析の手順を示す説明図である。

【0055】

まずＣＰＵ２１は、データ解析装置５が取得したレーダ画像（ＳＡＲ画像４１）を、幾何歪みあり又は幾何歪みなしの地図座標のレーダ画像に変換するための変換テーブルを生成する（ステップＳ１）。

【0056】

すなわち、ＣＰＵ２１の第１の地図投影部２１１は、通信部２９を介して衛星ＳＡＲ１からレーダ画像データ（ＳＡＲ画像データ）に基づくスラントレンジ座標のＳＡＲ画像４１を取得する。そして、第１の地図投影部２１１は、スラントレンジ座標のＳＡＲ画像４１を地図座標に変換してＧＥＣ画像４２を生成するとともに、第１の変換テーブル２１１ａ及び第１の逆変換テーブル２１１ｂを生成する。以下の説明では、第１の変換テーブル２１１ａ及び第１の逆変換テーブル２１１ｂを、ＧＥＣ画像変換テーブル５０ａ及びＧＥＣ画像逆変換テーブル５０ｂと記している。

【0057】

一方、第２の地図投影部２１２は、通信部２９を介して衛星ＳＡＲ１からレーダ画像データ（ＳＡＲ画像データ）に基づくスラントレンジ座標のＳＡＲ画像４１を取得する。そして、第２の地図投影部２１２は、スラントレンジ座標のＳＡＲ画像４１をオルソ補正した後に地図座標に変換してオルソ画像４３を生成するとともに、第２の変換テーブル２１２ａ及び第２の逆変換テーブル２１２ｂを生成する。以下の説明では、第２の変換テーブル２１２ａ及び第２の逆変換テーブル２１２ｂを、地図座標変換テーブル５１ａ及び地図座標逆変換テーブル５１ｂと記している。

【0058】

次に、ＣＰＵ２１の歪みあり座標変換部２１３は、歪みあり座標変換テーブル４５に基づいて、光学画像又は地理情報等の歪みなしの地図座標の参照画像４４を、歪ありの地図座標の比較画像４６に変換する（ステップＳ２）。参照画像４４は、既存の情報１３（図２）に相当し、例えば記録部２８に保存されている、又は通信部２９を介して外部から供給される。歪みあり座標変換テーブル４５は、地図座標逆変換テーブル５１ｂとＧＥＣ画像変換テーブル５０ａから構成されている。

【0059】

次に、ＣＰＵ２１は、ＧＥＣ画像４２と比較画像４６を出力部２７に出力する。出力部２７としてのディスプレイは、その表示画面にＧＥＣ画像４２と比較画像４６を並べて又は重ね合わせて表示する。ユーザは、表示画面に表示されたＧＥＣ画像４２を比較画像４６と比較することにより、ＧＥＣ画像４２に災害発生地等の変化域があるかどうかを判読する（ステップＳ３）。

【0060】

ユーザは、入力部２６を操作してＧＥＣ画像４２上の災害発生地を指定し、判読結果４７としてＧＥＣ画像４２上の災害発生地を歪みなし座標変換部２１４へ入力する。歪みなし座標変換部２１４は、判読結果４７が入力されると、歪みなし座標変換テーブル４８に基づいて、ＧＥＣ画像４２上の指定された地点（領域）を、幾何歪みなしの地図座標（判読結果４９）に変換する（ステップＳ４）。歪みなし座標変換テーブル４８は、ＧＥＣ画像逆変換テーブル５０ｂと地図座標変換テーブル５１ａから構成されている。

【0061】

次に、ＣＰＵ２１は、ＧＥＣ画像４２上の指定された地点（判読結果４７）を幾何歪みなしの地図座標に変換した情報（判読結果４９）を、地理情報（例えばＧＩＳ情報）として出力する（ステップＳ５）。ＣＰＵ２１は、この判読結果４９を地理情報として出力後、この処理を終了する。

【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】

また、幾何歪みありの地図座標のレーダ画像としてＧＥＣ画像を例示したが、これに限られない。幾何歪みありの地図座標のレーダ画像としては、Cosmo-SkyMed衛星で用いられるLevel 1C、RADARSAT-2で用いられるSSG（Map Image），SPG（Precision Map Image）、ALOS，ALOS-2で用いられるLevel 1.5などが挙げられる。

【符号の説明】

【0094】

１…衛星ＳＡＲ、 ２…データ解析センタ、 ３…地表、 ４…地上アンテナ、 ５…データ解析装置、 １１，１２…レーダ画像、 １３，１４…既存の情報、 １５，１６…判読結果、 ２１，２１Ａ…ＣＰＵ、 ４１…ＳＡＲ画像、 ４２…ＧＥＣ画像、 ４３…オルソ画像、 ４４…参照画像、 ４５…歪みあり座標変換テーブル、 ４６，４６Ａ…比較画像、 ４７，４９…判読結果、 ４７Ａ，４９Ａ…抽出結果、 ４８…歪みなし座標変換テーブル、 ５０ａ…地図座標変換テーブル、 ５０ｂ…地図座標逆変換テーブル、 ７１…レーダ画像、 ７２…既存の情報、 ７３，７４…抽出結果、 ８０…地物情報判読画面、 ８７…比較画像、 ８８…ＳＡＲ画像、 ２１１…第１の地図投影部、 ２１１ａ…第１の変換テーブル、 ２１１ｂ…第１の逆変換テーブル、 ２１２…第２の地図投影部、 ２１２ａ…第２の変換テーブル、 ２１２ｂ…第２の逆変換テーブル、 ２１３…歪みあり座標変換部、 ２１４…歪みなし座標変換部、 ２１５…参照画像

【図1】

【図2】

【図3】

【図4】

【図5】

【図7】

【図8】

【図9】

【図10】

【図11】

【図12】

【図13】

【図14】

【図15】

【図6】

【図16】

【図17】

【図18】

【図19】