(57)【要約】【課題】気象衛星の画像データを用いて地震の発生を的確に検出可能な画像解析装置およびその方法、並びにプログラムを提供する。【解決手段】気象衛星による気象画像データを解析する気象画像解析装置１０は、地形データに取得した気象画像データを載せ、気象画像データを複数の明度値を設定することによりこの設定値より明度の低い気象画像データを削除してスキャン処理を行い、地震の有無を、明度差を利用して表出させることが可能な処理部１７を有する。【選択図】図１
2008209417.tif 000002
【特許請求の範囲】
【請求項１】
気象衛星による気象画像データを解析する画像解析装置であって、
地形画像データに取得した気象画像データを載せ、当該気象画像データを複数の明度値を設定することにより当該設定値より明度の低い気象画像データを削除してスキャン処理を行い、地震の有無を明度差を利用して表出させることが可能な処理部
を有する画像解析装置。
【請求項２】
前記処理部は、
取得した前記気象画像データのデータ間を線形補間してデータ数を増やす処理を行う
請求項１記載の画像解析装置。
【請求項３】
前記処理部は、
高速フーリエ変換によりフーリエ級数近似曲線を作成し、前記気象画像データの量を増やす
請求項２記載の画像解析装置。
【請求項４】
前記処理部は、
前記気象画像データ中の明度の高い少なくとも経線、緯線、地形境界線を削除して、当該線の周辺の雲または隣り合う明度に置き換えて画像に修正して３次元的に表出させる  請求項１から３のいずれか一に記載の画像解析装置。
【請求項５】
前記処理部は、
前記地形画像データから海域と湖とを示すデータを削除して３次元的に表出させる
請求項１から４のいずれか一に記載の画像解析装置。
【請求項６】
前記処理部は、
気象画像データを地形画像データ上に載せ、当該画像上をウォークスルーしたときの所定時刻の画像を表示可能である
請求項１記載の画像解析装置。
【請求項７】
気象衛星による気象画像データを解析する気象画像解析方法であって、
前記地形画像データに取得した気象画像データを載せ、
当該気象画像データを複数の明度値を設定することにより当該設定値より明度の低い気象画像データを削除してスキャン処理を行い、
地震の有無を、明度差を利用した表出する
画像解析方法。
【請求項８】
地形画像データに比べてピクセル数の少ない気象画像データのデータ間を線形補間し高速フーリエ変換によりフーリエ級数近似曲線を作成し、前記気象画像データの量を増やす  請求項７記載の画像解析方法。
【請求項９】
前記気象画像データ中の明度の高い少なくとも経線、緯線、地形境界線を削除して、当該線の周辺の雲または隣り合う明度に置き換えて画像に修正して３次元的に表出する
請求項７または８記載の画像解析方法。
【請求項１０】
前記地形画像データから海域と湖とを示すデータを削除して３次元的に表出する
請求項７から９のいずれか一に記載の画像解析方法。
【請求項１１】
気象衛星による気象画像データを解析する処理をコンピュータに実行させるプログラムであって、
地形画像データに取得した気象画像データを載せ、当該気象画像データを複数の明度値を設定することにより当該設定値より明度の低い気象画像データを削除してスキャン処理を行い、地震の有無を、明度差を利用して表出させる処理をコンピュータに実行させるプログラム。
【請求項１２】
前記地形画像データに比べてピクセル数の少ない気象画像データのデータ間を線形補間し高速フーリエ変換によりフーリエ級数近似曲線を作成し、前記気象画像データの量を増やす処理をコンピュータに実行させる
請求項１１記載のプログラム。
【請求項１３】
前記気象画像データ中の明度の高い少なくとも経線、緯線、地形境界線を削除して、当該線の周辺の雲または隣り合う明度に置き換えて画像に修正して３次元的に表出する処理をコンピュータに実行させる
請求項１１または１２記載のプログラム。
【請求項１４】
前記地形画像データから海域と湖とを示すデータを削除して３次元的に表出する処理をコンピュータに実行させる
請求項１１から１３のいずれか一に記載のプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】  本発明は、地震の発生を検出可能な画像解析装置およびその方法、並びにプログラムに関するものである。
【背景技術】
【０００２】  最近、台風による突風や異常気象による竜巻が数多く発生している。竜巻発生時にはその規模により、特殊な気象画像が確認されている。
気象衛星の画像は、地形データに重ね合わせて用いられるが、地形データに対応したいずれの位置に雲が存在するのか、また、台風級の大きな勢力の低気圧に伴う雲の流れや台風の目の存在、気圧配置等を確認することができる。【０００３】  このような、気象衛星の画像データは、所定時間ごと（現在は３０分ごと）に更新され、天気予報や台風の進路予想等に使用されている。【０００４】  気象衛星の画像としては、可視画像と赤外線画像がある。赤外線画像は温度観測により、雲頂の高度を見積もり、可視画像とは異なる。ただし、夜間の観測も可能である。
ある地域の一定時間継続した特殊な気象画像では、その地域で際立って明度の差が生じている。ここでは、そのような特殊な画像を異常気象画像と称する。
２００７年３月２５日９時４１分に能登半島沖でマグニチュード６．９の地震が発生した。この地域のその時刻の前後２時間で同様な異常気象画像が確認された。【０００５】
【非特許文献１】国土地理院：数値地図２５０ｍメッシュ（標高）、平成９年７月
【非特許文献２】杉本智彦：カシミール３Ｄ，実業之日本社，Ｄｅｃ，２００３．
【非特許文献３】世界を見渡す３Ｄソフトウェア，ｈｔｔｐ：／／ｅａｒｔｈ．ｇｏｏｇｌｅ．ｃｏ．ｊｐ
【非特許文献４】活断層研究会編、新編日本の活断層−分布図と資料、東京大学出版会、４．１９９５
【発明の開示】
【発明が解決しようとする課題】
【０００６】  ところで、前述したように、気象衛星の画像は、地形データに対応した雲の位置や、大きな勢力の低気圧に伴う雲の流れや台風の目の存在を確認することができるが、地震等の発生は、その画像の特質から確認、あるいは予測することが極めて困難である。【０００７】  地震の検出のために、気象衛星の画像を３次元表示することにより、解明の糸口を得ることができる可能性があるものと考えられる。
しかしながら、たとえば大きな勢力の台風を含む気象衛星の可視画像は、九州や四国の全体を覆うような厚い雲が存在する場合、明度が高く明度差が少ないため、３次元表示をすることは困難である。【０００８】  ３次元表示がより明確な赤外線画像を利用することも可能であるが、地形データに比べて、気象データはそのピクセル数が少なく粗いために、現状では気象異常状態の発生検出に用いることは難しい。【０００９】  本発明は、気象衛星の画像データを用いて地震等の異常状態の発生を的確に検出可能な画像解析装置およびその方法、並びにプログラムを提供することにある。
【課題を解決するための手段】
【００１０】  本発明の第１の観点は、気象衛星による気象画像データを解析する画像解析装置であって、地形画像データに取得した気象画像データを載せ、当該気象画像データを複数の明度値を設定することにより当該設定値より明度の低い気象画像データを削除してスキャン処理を行い、地震の有無を明度差を利用して表出させることが可能な処理部を有する。【００１１】  好適には、前記処理部は、取得した前記気象画像データのデータ間を線形補間してデータ数を増やす処理を行う。【００１２】  好適には、前記処理部は、高速フーリエ変換によりフーリエ級数近似曲線を作成し、前記気象画像データの量を増やす。【００１３】  好適には、前記処理部は、前記気象画像データ中の明度の高い少なくとも経線、緯線、地形境界線を削除して、当該線の周辺の雲または隣り合う明度に置き換えて画像に修正して３次元的に表出させる。【００１４】  好適には、前記処理部は、前記地形画像データから海域と湖とを示すデータを削除して３次元的に表出させる。【００１５】  好適には、前記処理部は、気象画像データを地形画像データ上に載せ、当該画像上をウォークスルーしたときの所定時刻の画像を表示可能である。【００１６】  本発明の第２の観点は、気象衛星による気象画像データを解析する画像解析方法であって、前記地形画像データに取得した気象画像データを載せ、当該気象画像データを複数の明度値を設定することにより当該設定値より明度の低い気象画像データを削除してスキャン処理を行い、地震の有無を、明度差を利用して表出する。【００１７】  本発明の第３の観点は、気象衛星による気象画像データを解析する処理をコンピュータに実行させるプログラムであって、地形画像データに取得した気象画像データを載せ、当該気象画像データを複数の明度値を設定することにより当該設定値より明度の低い気象画像データを削除してスキャン処理を行い、気象異常の有無を、明度差を利用して表出させる処理をコンピュータに実行させるプログラムである。
【発明の効果】
【００１８】  本発明によれば、気象衛星の画像データの明度差を利用することにより地震の発生を的確に検出することができる。
【発明を実施するための最良の形態】
【００１９】  以下、本発明の実施形態を図面に関連付けて説明する。【００２０】  図１は、本発明の実施形態に係る気象画像解析装置（気象異常検出装置）の構成例を示すブロック図である。【００２１】  画像解析装置（異常検出装置）１０は、図１に示すように、地形データ記憶部１１、気象画像データ提供部１２、気象画像データ取得部１３、記憶部１４、操作部１５、表示部１６、および処理部１７を主構成要素として有している。【００２２】  地形データ記憶部１１は、５ｍから１ｋｍメッシュまで種々の数値地図として国土地理院から公開（たとえば非特許文献１参照）されている地形データのうち、たとえば１ｋｍメッシュの地形データを記憶する。
地形データ記憶部１１は、たとえばＣＤＲＯＭ等の記憶媒体として形成され、処理部１７によりアクセス可能な図示しない光記録再生装置等にセットされ、処理部１７により記憶データが読み出される。
この読み出される地形データの数値データは、地形標高の他に標高０ｍの海域は−９９９９、湖は０００００で表示されように形成されている。【００２３】  気象画像データ提供部１２は、気象衛星を通して撮像された気象画像データを記憶するデータベースとして形成され、気象画像データ取得部１３を介して取得可能に構成されている。
気象画像データ提供部１２は、気象衛星を通して赤外線画像あるいは可視画像として撮像され、３０分ごとに更新され、所定期間のデータが蓄積される。【００２４】  気象画像データ取得部１３は、処理部１７による気象画像データ取得の指示に応答して、気象画像データ提供部１２にデータ要求等の送受によるネゴシエーション（データの交換手続）等を行って気象画像データを取得する。
気象画像データ取得部１３により取得された気象画像データは、記憶部１４に格納される。【００２５】  記憶部１４は、処理部１７のメインプログラムや気象画像の解析プログラム、計算プログラム等を記憶する。
また、記憶部１４は、画像解析のために取得した気象画像データ、地形データ、あるいは解析処理中のデータ等が一時的にあるいは保存データとして記憶される。
記憶部１４は、処理部１７によりアクセスされ、データが読み出されたり、書き込まれたりする。【００２６】  操作部１５は、キーボードやマウス等により構成され、処理部１７において、所望のプログラムを実行するためのコマンドや設定値を入力し、処理部１７に指示するために使用される。【００２７】  表示部１６は、たとえば、ＣＲＴ（Ｃａｔｈｏｄｅ  Ｒａｙ  Ｔｕｂｅ）や液晶表示パネル（ＬＤＣ）等の各種表示装置により構成され、表示部１６のモニタ部には、取得した画像解析のために取得した気象画像データ、地形データ、解析処理のデータ、解析結果が表示される。
また、表示部１６のモニタ部には、処理部１７への処理の指示のためのＧＵＩ（Ｇｒａｐｈｉｃａｌ  Ｕｓｅｒ  Ｉｎｔｅｒｆａｃｅ）画面が表示される。【００２８】  処理部１７は、解析装置１０全体の各部の制御を行う。
処理部１７は、操作部１５の操作入力に応答して気象画像解析のため、地形データ記憶部１１による地形データの取得、気象画像データ取得部１３を通しての気象画像データの取得、取得データの記憶部１４への格納処理、読み出し処理を行う。
さらに、処理部１７は、以下に詳述するような、取得したデータのモディファイ処理、合成処理、気象画像データの複数の明度値をしきい値としたスキャン処理等を行って、気象画像データの解析処理を行い、取得した画像解析のために取得した気象画像データ、地形データ、解析処理のデータ、解析結果等の表示部１６への３次元表示制御等を行う。【００２９】  以下、本実施形態に係る気象画像データの解析処理について、処理部１７の処理を中心に、図面に関連付けて説明する。【００３０】  まず、処理部１７の処理の大略を説明する。【００３１】＜処理の概略＞
最近、地形データの３次元表示が数多く試みられている（非特許文献２、あるいは非特許文献３参照）。
その地形データは５ｍから１ｋｍメッシュまで種々の数値地図として国土地理院から公開（前述の非特許文献１）されている。
本実施形態においては、たとえば、１ｋｍメッシュを用いて地形を３次元表示する。この数値データは地形標高の他に標高０ｍの海域は−９９９９、湖は０００００で表示されている。まず、これらの海域と湖を地形データから削除して３次元表示する。海域と湖を省くことで地形の表示がより明確になり、画像表示時間も短縮できる。【００３２】  次に、本実施形態においては、気象衛星ひまわりの気象画像データを地形データ上に載せて、その明度差を利用して３次元表示する。その際、明度の低いデータを削除して表示する。これにより、明度の高い厚い雲等が表示可能になる。なお、画像の明度は、０〜２５５の２５６の値をとることが可能であり、明度値２５５に近いほど画像は白く表示され、明度値０に近いほど画像は黒く表示される。
また、地形画像データと気象画像データの同じ領域において、地形画像データに比べて、気象画像データはそのピクセル数が少なく粗いために、後で詳述するように、データ間を線形補間して高速フーリエ変換によりフーリエ級数近似曲線を作成し、そのデータ数を増やす処理を行う。【００３３】  以下の処理部１７における気象画像データの解析処理のより具体的な処理について、順を追って説明する。
なお、ここでは、２００６年９月の台風１３号について１７日１１時から２３時の間の九州地形上の気象画像データを解析し、地形上のその変動について考察する。【００３４】＜地形データの海域と湖の削除＞
上述したように、１ｋｍメッシュの数値地図は地形標高の他に標高０ｍの海域は−９９９９、湖は０００００で表示されている。その行列数は８０×８０で、そのユニットを連結して求める地形が表示できる。
ここでは、そのユニットを連結して九州の地形データの行列を作成する。地形標高データを３次元表示するために、そのデータをＤＸＦ（Ｄｒａｗｉｎｇ  Ｉｎｔｅｒｃｈａｎｇｅ  Ｆｉｌｅ）変換して３次元ＣＡＤで読み込み、地形メッシュを作成する。
その際、湖と海域の標高０ｍ以下のメッシュを削除する。【００３５】  図２（Ａ），（Ｂ）は地形データのメッシュにその番号と節点番号、さらに座標値を付加した例を示す図であって、図２（Ａ）はメッシュの削除前を示し、図２（Ｂ）はメッシュの削除後を示している。【００３６】  地形データは、図２（Ａ）に示すように、メッシュにその番号と節点番号、さらに座標値を付加したものである。その後、標高０の節点を含むメッシュ、たとえば、図２（Ａ）で番号６，７，１０，１１のメッシュを削除する。
図２（Ｂ）に示すように、削除されたメッシュ番号を除いて、メッシュを通し番号に並べ替え、その座標値を入れ替える。さらに節点番号を通し番号に並べ替える。
これにより、節点番号が１から昇順に付加されてＤＸＦ変換が可能になる。【００３７】  図３（Ａ），（Ｂ）は、海域と湖の削除前後の九州の地形画像を示す図であって、図３（Ａ）が海域と湖の削除前の画像を示す図であり、図３（Ｂ）が海域と湖の削除後の画像を示す図である。
図４は、海域と湖の削除後の九州の３次元表示画像を示す図である。【００３８】  図３および図４は、行列数８０×８０のユニットを５行４列連結したもので、海域と湖のメッシュ削除前と削除後の表示である。
カラー表示を行う場合、たとえば標高の低所から高所の順に、青から赤色の順で着色して表示することが可能である。【００３９】＜気象画像データの３次元表示＞
図５（Ａ），（Ｂ）は、経線、緯線の削除前後の赤外線画像を示す図であって、図５（Ａ）が経線、緯線の削除前の画像を示す図であり、図５（Ｂ）が経線、緯線の削除後の画像を示す図である。
気象衛星の画像には雲画像の他に、明度の高い経線、緯線、地形境界線などが含まれているので、それらの線を削除して雲の明度のみを表示する必要がある。
そのため、図５（Ａ），（Ｂ）に示すように、それらの線の周辺の雲または隣り合う明度に置き換えて、余計なデータを削除した画像に修正する。
また、気象衛星の画像は可視画像と赤外線画像がある。赤外線画像は温度観測により、雲頂の高度を見積もり、可視画像とは異なる。ただし、夜間の観測も可能である。【００４０】  ２００６年９月の台風１３号はその勢力が強く、最低気圧は９１９ｈＰａ（１６日３〜９時）で、その勢力はあまり衰えずに９５０ｈＰａ（１７日１８時）で九州北部へ上陸した。その可視画像は九州全体を覆い明度差が少ないため、明度差を利用した３次元表示ができない。
そこで、ここでは３次元表示がより明確な赤外線画像を利用する。そして、図５（Ｂ）に示すように、経線、緯線の削除後の赤外線画像を利用する。【００４１】  また、地形と気象画像の同じ領域において、地形データに比べて、気象画像データはそのピクセル数が少なく表示が粗くなる。
そこで、本実施形態においては、データ間を線形補間して高速フーリエ変換によりフーリエ級数近似曲線を作成し、データ数を増やして表示する。【００４２】  次に、地形データの海域と湖の削除と同様に、明度の低い気象画像データを削除して表示する。
これにより、明度の高い厚い雲等が表示可能になる。
図６は線形補間前後の赤外線画像の３次元表示例である。【００４３】  図６（Ａ），（Ｂ）は、線形補間前後の赤外線画像表示を示す図であって、図６（Ａ）が画像のピクセル数の少ない線形補間前（７６行６２列）の赤外線画像を示す図であり、図６（Ｂ）が線形補間後（３倍，６０％）の赤外線画像を示した図である。なお、図６（Ａ），（Ｂ）は、明度値２３０以下を削除して表示した画像を示している。
図６（Ａ），（Ｂ）に示すように、ピクセル数の少ない線形補間前の画像を線形補間すると、線形補間後の画像はその３倍（級数の項数６０％）のピクセル数で、細密で滑らかな曲面になっている。【００４４】  本実施形態においては、気象画像データを線形補間し、たとえば等間隔にサンプリングされている離散画像データ群、あるいは、不等間隔にサンプリングされている離散画像データ群を元にして、データ数が画像データ群のデータ数の所定倍の回帰用データを生成し、この回帰用データを、項数が画像データ群のデータ数よりも小さいフーリエ級数に回帰させて元の画像データ群を近似する。その結果、近似により得られたデータ群の各データは滑らかに変化するようになり、かつ、近似の精度が向上する。【００４５】  すなわち、本実施形態においては、画像データ群をフーリエ級数に回帰させることによる近似を行う。
離散データ群をフーリエ級数へ回帰させる手法について、以下に簡単に説明する。【００４６】  フーリエ級数を用いると、回帰曲線により再現しようとする離散データ群の曲線を、三角関数を用いて１つの式によって表現することができる。
等間隔にＮ個のデータがサンプリングされる関数ｘ（ｔ）をフーリエ級数で表わすとすると、Ａ０，Ａ１，Ａ２，…，ＡｋおよびＢ０，Ｂ１，Ｂ２，…，Ｂｋを適宜決められる定数としたとき、ｘ（ｔ）＝Ａ０＋Ａ１ｃｏｓｔ＋Ａ２ｃｏｓ２ｔ＋…＋Ａｋｃｏｓｋｔ＋…＋Ｂ０＋Ｂ１ｓｉｎｔ＋Ｂ２ｓｉｎ２ｔ＋…＋Ｂｋｓｉｎｋｔ＋…として表わすことができる。
上式をまとめると下記式（１）のようになる。【００４７】
【数１】
2008209417.tif 000003
【００４８】  式（１）において、Δｔがサンプリング間隔を示し、ＮΔｔが関数ｘ（ｔ）の継続区間を示している。
式（１）はｋについて０から無限大まで総和している無限級数であるが、ｋ＝Ｎ／２までで打ち切ると、下記式（２）で表わされる有限三角級数となる。【００４９】
【数２】
2008209417.tif 000004
【００５０】  上式（２）によって表わされる関数を、サンプリング間隔ΔｔでサンプリングされたＮ個のサンプル値ｘｍ（ｍ＝０，１，２，…，Ｎ−１）を通るような関数に直すと、各サンプル値ｘｍは以下の式（３）のように表わすことができる。【００５１】
【数３】
2008209417.tif 000005
【００５２】  上式（３）における係数Ａｋ、Ｂｋはそれぞれ下記式（４）、（５）で表わすことができる。【００５３】
【数４】
2008209417.tif 000006
【００５４】  上式（３）で表わされる関数は、各サンプル値ｘｍを全て正確に通過する。したがって、上式（３）を用いれば、画像データ群をフーリエ級数に回帰させて、近似曲線を得ることができる。【００５５】  ただし、式（３）の関数をそのまま用いて画像データ群を全て通過するように回帰させると、フーリエ級数に回帰させて得られた近似曲線は図７（Ａ）に示すグラフＧＤＣ１のように細かい増減を繰返す曲線となってしまう。
フーリエ級数近似後のデータの値が細かく増減することを防止して滑らかな近似曲線が得られるようにするために、本実施形態においては、式（３）の関数の項数、すなわち係数Ａｋ、Ｂｋの数を、ステップＳＴ１において得られた画像データ群のデータ数よりも小さい所定の数にして画像データ群をフーリエ級数に回帰させる。【００５６】  フーリエ級数の項数の数については、画像データ群のデータ数の約数十％程度にすることが、等高線等の特定のためには好ましい。
画像データ群のデータ数が少な過ぎると、十分な近似精度が得られなくなるほど項数ｐの数が少なくなる。したがって、必要な近似精度にもよるが、項数ｐを増加してある程度の近似精度を確保するために、元となる画像データ群のデータ数もある程度必要である。画像データ群のデータ数が多いほど近似精度が良くなる。【００５７】  項数ｐが少なくなるほど、式（３）の関数の計算時間は短くなり、得られる近似曲線は滑らかになる傾向にある。
ただし、項数ｐが少なくなると画像データ群の後ろの部分において近似値との間の誤差が大きくなってくる。そこで、近似の精度を高めるために、本実施形態においては画像データ群に付加的な付加データ群を加えて、フーリエ級数に回帰させるデータ数を画像データ群の所定倍に増やす。【００５８】  処理部１７は、フーリエ級数に回帰させるデータ数を増やすために、得られた画像データ群に加える付加データ群を生成する。
フーリエ級数は周期関数であるため、この性質を利用して近似の精度を高めるために、付加データ群は、その個々のデータが、その大きさについて画像データ群の最後のデータを中心として画像データ群の各々のデータの大きさに対称となるように配列されたデータ群として生成する。
本実施形態においては、付加データ群の個々のデータを破線により繋いだグラフＧＤＣ１’が、図７（Ａ）に示すように画像データ群のグラフＧＤＣ１の最後の点を通る縦線に対してグラフＧＤＣ１に線対称になるように付加データ群を生成する。
グラフＧＤＣ１’がグラフＧＤＣ１の最後の点を中心としてグラフＧＤＣ１に点対称になるように付加データ群を生成してもよいが、近似精度向上の観点からは、図７（Ａ）に示すようにグラフが線対称となるようなデータを付加データ群とすることが好ましい。【００５９】  以上のように、画像データ群の最後のデータに対称となるように配列した付加データ群を画像データ群に付加することより、式（３）におけるサンプル値の数をＮ個ではなくαＮ個とすることができる。αＮ個となったデータ群を、回帰用データ群と呼ぶ。
ここで、係数αは、データ数を増やすという観点から自然数とする。特に、フーリエ級数の周期的性質を利用して近似の精度を向上させるという観点からは、係数αを偶数とすることが好ましい。係数αを偶数（たとえば２）とすれば、元の画像データ群の全てのデータが、画像データ群の最後のデータを中心として対称な大きさを持つデータとして全て用いられるためである。
ただし、データ数を増やすという観点からは、係数αを正の小数とすることもできる。
係数αを２以上にして最初のＮ個のデータの２倍以上のデータ数にするためには、最初の画像データ群に基づいて生成した回帰用データ群を第１のデータ群と考え、このデータ群の最後のデータを中心として対称となるように新たなデータ群を生成すればよい。【００６０】  以上のように生成して得た回帰用データ群を、処理部１７がフーリエ級数に回帰させる。【００６１】  回帰用データ群をフーリエ級数に回帰させて得られる近似曲線のグラフＦＧＤＣ１を図７（Ｂ）に示す。図７（Ｂ）の横軸は図７（Ａ）と同じ測定軸ＭＡＬ上の基準位置Ａからの距離［ｍｍ］を表わしており、縦軸は輝度の大きさを表わしている。
回帰用データ群をフーリエ級数に回帰させると回帰用データ群全体についての近似曲線を得ることができるが、等色線の位置特定のためには測定軸ＭＡＬ上の画像データ群に関する解析結果のみが必要であるため、図７（Ｂ）には元の画像データ群に関する部分の近似曲線のグラフＦＧＤＣ１のみを示している。【００６２】  グラフＦＧＤＣ１から明らかなように、項数ｐをある程度小さくしてフーリエ級数に回帰させることにより、各データの値は滑らかに変化するようになる。
また、付加データ群を含めて数が増えたデータを用いて回帰させているため、フーリエ級数への近似後の各データは滑らかに変化しつつも元の画像データ群との誤差は小さくなり、近似の精度が向上している。
元の画像データ群を用いることにより、生成した回帰用データ群が周期的なデータになるほど、近似の精度は向上すると考えられる。ただし、得られる近似曲線の滑らかさと近似の精度は、画像データ群の変化の特性や項数ｐ、回帰用データ群のデータ数等のパラメータの組み合わせにより複雑に変化し、一概に決まるものではない。【００６３】  そして、このようにして得られたフーリエ近似曲線による輝度値をＤＸＦ変換して３次元のメッシュ画像として表示させた。なお、３次元のメッシュ画像の表示にはＡｒｃｈｉＣＡＤ７（グラフィクソフト社）等を使用することが可能であるが、他の３次元表示ソフトを使用しても構わない。【００６４】  次に、地形画像と気象画像の３次元表示について説明する。【００６５】＜地形画像と気象画像の３次元表示＞
気象画像の基点として日本地形の区域では、東経１３０度、北緯３０度、東経１４０度、北緯３０度、東経１３０度、北緯４０度および東経１４０度、北緯４０度の４交点が利用可能である。
図８は、気象画像の基点として日本地形の区域での４交点を示す図である。
処理部１７は、この交点を基点として地形画像と気象画像を一致させて両画像を表示部１６に表示させる。
九州の地形ではその交点の東経１３０度、北緯３０度を基点として利用する。また、地形データに対して気象データを約１．７倍拡大して大きさを地形と一致させている。【００６６】  以下、上記実施の形態を実際の地震の検出に適用した具体例を述べる。【００６７】能登半島沖の地震
２００７年０３月２５日９時４１分頃、石川県で地震が発生した。震源地は能登半島沖（北緯３７．３°、東経１３６．５°）で震源の深さは約１１ｋｍ、地震の規模（マグニチュード）は６．９と推定されている。七尾市田鶴浜町で最大震度は６強である。気象庁発表の各地の震度に関する情報、地震情報を図９、１０に示す。震源地の茨城県沖を基点に、震度１、２が広範囲に分布している。【００６８】  石川県周辺の活断層、地震分布を図１１に示す（非特許文献４）。石川県、新潟県付近に活断層が分布している。【００６９】  ９時前後の気象画像について、能登半島付近の画像を詳細に検討する。各時刻の気象画像を図１２に示す。図９の震度分布付近で、画像に明度の変化が生じている。【００７０】  図１３は石川県付近を拡大した画像で、右は経線等を削除後の画像である。気象画像と地形との対応位置として、２個の基点を←で示す。【００７１】  図１４〜２１は石川県付近の地形上に各時刻の気象画像を載せた画像である。その平面（上）と断面で、↑は視線方向を示す。
図１４、１５は７時３０分の表示で平面は明度値１９０以下を削除した画像である。削除されたＡ地域の断面の明度値の最低は１５９、最高は１８１で、その明度差は２２である。ここで最高は最低値周辺の最高値を東西南北求めて平均した。その値は東１７１、西１９１、南１６６、北１９７である。断面は南北の方向について表示している。【００７２】  同様に図１６、１７は８時００分の表示で平面は明度値１９５以下を削除した画像である。明度の低い地域がＡ、Ｂ、Ｃ、Ｄの４ヶ所ある。東西断面▲１▼と南北断面▲２▼、▲３▼で、その地域が表示されている。Ｂ地域の断面の明度値の最低は１４９、最高は１９７で、その明度差は４８である。【００７３】  図１８は８時３０分の表示で平面は明度値１９０以下を削除した画像である。Ｆ地域の断面の明度値の最低は１４０、最高は１９９で、その明度差は５９である。【００７４】  図１９は９時００分の表示で平面は明度値１９０以下を削除した画像である。Ａ地域の断面の明度値の最低は１３２、最高は２１６で、その明度差は８４である。この時刻でＡ地域の明度差が最大になる。東西断面で約１３９ｋｍの範囲で、明度値の急激な変化がこの地域に生じている。この時刻の４１分後頃に、Ａ地域の海岸沖で地震が発生している。【００７５】  図２０は１０時００分の表示で平面は明度値１８０以下を削除した画像である。Ｂ地域の断面の明度値の最低は１３８、最高は２０３で、その明度差は６５である。その値は９時００分の最大値から減少している。【００７６】  図２１は９時００分の３Ｄ表示でＡ地域を表示している。図上は南西から北東方向の視線である。下は南東から北西方向の視線である。Ａ地域は地形の海岸沿上にある。【００７７】  図２２〜２７は９時００分の平面の表示で明度値１３３〜１８０以下を削除した画像である。明度の低い地域が能登半島の内陸部から震源地の海岸へと広がっている。【００７８】  以上の分析から明らかなように、能登半島沖の地震では、震源地付近の明度差が大で、その範囲が広い地域の震度は大きいことが判った。
このように気象画像データの明度差から地震の位置を特定（検出）することができると判った。【００７９】茨城県南部の地震
２００７年６月４日１３時３４分頃、茨城県で地震が発生した。震源地は南部（北緯３６．１°、東経１３９．９°）で震源の深さは約５０ｋｍ、地震の規模（マグニチュード）は４．４と推定されている。足利市大正町で最大震度は３である。気象庁発表の各地の震度に関する情報、地震情報を図２８に示す。震源地の茨城県沖を基点に、震度１、２が広範囲に分布している。【００８０】  茨城県周辺の活断層、地震分布を図２９に示す。茨城県南部付近に活断層が分布している。【００８１】  １３時前後の気象画像について、茨城県南部付近の画像を詳細に検討する。各時刻の気象画像を図３０に示す。図２８の震度分布付近で、画像に明度の変化が生じている。【００８２】  図３１は茨城県付近を拡大した画像で、右は経線等を削除後の画像である。気象画像と地形との対応位置として、２個の基点を←で示す。【００８３】  図３２〜３８は茨城県付近の地形上に各時刻の気象画像を載せた画像である。その平面（上）と断面で、↑は視線方向を示す。
図３２は１１時３０分の表示で平面は明度値５１以下を削除した画像である。削除されたＡ地域の断面の明度値の最低は２６、最高は６１で、その明度差は３５である。ここで最高は最低値周辺の最高値を東西南北求めて平均した。その値は東５６、西５４、南４９、北８５である。断面は東西、南北の２方向について表示している。【００８４】  同様に図３３は１２時００分の表示で平面は明度値５２以下を削除した画像である。明度の低い地域がＡ、Ｂ、Ｃ、Ｄ、Ｅの５ヶ所ある。東西断面▲１▼と南北断面▲２▼、▲３▼で、その地域が表示されている。Ｂ地域の断面の明度値の最低は２４、最高は６２で、その明度差は３８である。【００８５】  図３４は１２時３０分の表示で平面は明度値５３以下を削除した画像である。Ａ地域の断面の明度値の最低は２４、最高は６２で、その明度差は３８である。【００８６】  図３５は１３時００分の表示で平面は明度値５０以下を削除した画像である。Ａ地域の断面の明度値の最低は２３、最高は６３で、その明度差は４０である。【００８７】  図３６は１３時３０分の表示で平面は明度値５６以下を削除した画像である。Ａ地域の断面の明度値の最低は２２、最高は６６で、その明度差は４４である。この時刻でＡ地域の明度差が最大になる。東西断面で約１１４ｋｍ、南北断面で約１１３ｋｍの範囲で、明度値の急激な変化がこの地域に生じている。この時刻の４分後頃に、Ａ地域の東部で地震が発生している。【００８８】  図３７は１４時００分の表示で平面は明度値５０以下を削除した画像である。Ａ地域の断面の明度値の最低は２５、最高は６７で、その明度差は４２である。その値は１３時３０分の最大値から減少している。【００８９】  図３８は１３時３０分の３Ｄ表示でＡ地域を表示している。図上は北東から南西方向の視線である。図下は南東から北西方向の視線である。Ａ地域は地形の平野上にある。【００９０】  図３９〜４４は１３時３０分の平面の表示で明度値２２〜５１以下を削除した画像である。明度の低い地域が茨城県南部の内陸部から震源地の東部へと広がっている。【００９１】  以上から明らかなように、茨城県南部の地震ではマグニチュード４、震度３程度の地域の明度解析が可能である。【００９２】北海道根室支庁北部の地震
２００７年０７月０１日１３時１２分頃、北海道東部で地震が発生した。震源地は根室支庁北部（北緯４３．５°、東経１４４．９°）で震源の深さは約１３０ｋｍ、地震の規模（マグニチュード）は５．６と推定されている。釧路市幸町で最大震度は４である。
気象庁発表の各地の震度に関する情報、地震情報を図４５、４６に示す。震源地の根室支庁北部を基点に、震度１、２が広範囲に分布している。【００９３】  北海道周辺の活断層、地震分布を図４７に示す。石狩支庁、十勝支庁および根室支庁付近に活断層が分布している。【００９４】  １３時前後２時間の気象画像について、北海道の画像を詳細に検討する。各時刻の気象画像を図４８に示す。図４６の震度分布付近で、画像に明度の変化が生じている。【００９５】  図４９は北海道付近を拡大した画像で、右は経線等を削除後の画像である。気象画像と地形との対応位置として、２個の基点を←で示す。また１１時３０分から１４時３０分までの３０分毎の平均画像を示す。【００９６】  図５０〜５９は北海道の地形上に各時刻の気象画像を載せた画像である。その平面（上）と断面で、↑は視線方向を示す。
図５０は１１時３０分の表示で平面は明度値４８以下を削除した画像である。削除されたＡ地域の断面の明度値の最低は５０、最高は７４で、その明度差は２４である。ここで最高は最低値周辺の最高値を北、南求めて平均した。その値は北７０、南７７である。同様にＢ地域、Ｃ地域について各値を求めた。Ａ地域付近は石狩支庁、Ｂ地域は十勝支庁、Ｃ地域は震源地の根室支庁である。【００９７】  同様に図５１は１２時００分の表示で平面は明度値４８以下を削除した画像である。Ｂ地域の明度差は３７である。【００９８】  図５２は１２時３０分の表示で平面は明度値４７以下を削除した画像である。この時刻でＢ地域の明度差が最大値３９になる。約５７ｋｍの範囲で、明度値の急激な変化がこの地域に生じている。【００９９】  図５３は１３時００分の表示で平面は明度値４６以下を削除した画像である。Ｂ地域の明度差は３８である。１３時１２分頃にＣ地域で地震が発生している。【０１００】  図５４は１３時００分の平面の２方向について、その断面を求めた。断面１の明度差が断面２と比較して大きい。【０１０１】  図５５は１３時３０分の表示で平面は明度値４５以下を削除した画像である。平面の東西、南北の２方向について断面を求めた。南北断面では根室支庁北部に相当する左部分の明度差が大きい。Ｂ地域の明度差は３３である。【０１０２】  図５６は１４時００分の表示で平面は明度値４６以下を削除した画像である。Ｂ地域の明度差は３５である。【０１０３】  図５７は１４時３０分の表示で平面は明度値４８以下を削除した画像である。Ｂ地域の明度差は３６である。【０１０４】  図５８は１１時３０分〜１４時３０分の各明度値を平均した。平面は明度値５０以下を削除した画像である。Ｂ地域の明度差は３６である。震源地のＣ地域は２９でＢ地域の明度差よりも小さい。【０１０５】  図５９は１３時００分の平面の表示で明度値６１以下を削除した画像である。石狩支庁、十勝支庁、震源地の根室支庁付近が削除されている。【０１０６】  図６０は１３時００分の３Ｄ表示でＡ、Ｂ、Ｃ地域を表示している。図上は南西から北東方向の視線である。下は南東から北西方向の視線である。各地域は地形の平野上にある。【０１０７】  図６１〜６７は１３時００分の平面の表示で明度値４６〜５９以下を削除した画像である。明度の低い地域が石狩支庁付近から始まって十勝支庁、根室支庁へと広がっている。【０１０８】  以上の通り、北海道根室支庁北部の地震では、活断層の分布付近で明度が低く、その明度差が大である。【０１０９】新潟中越沖地震
２００７年０７月１６日１０時１３分頃、新潟県で地震が発生した。震源地は上中越沖（北緯３７．５°、東経１３８．６°、新潟の南西６０ｋｍ付近）で震源の深さは約１０ｋｍ、地震の規模（マグニチュード）は６．８と推定されている。飯綱町芋川で最大震度は６強である。気象庁発表の各地の震度に関する情報、地震情報を図６８、６９に示す。震源地の上中越沖を基点に、震度１、２および３が広範囲に分布している。【０１１０】  新潟県周辺の活断層、地震分布を図７０に示す。新潟県、石川県付近に活断層が分布している。【０１１１】  １０時前後の気象画像について、新潟県付近の画像を詳細に検討する。各時刻の気象画像を図７１に示す。図６８、６９の震度分布付近で、画像に明度の変化が生じている。【０１１２】  図７２は新潟県付近を拡大した画像で、右は経線等を削除後の画像である。気象画像と地形との対応位置として、２個の基点を←で示す。【０１１３】  図７３〜７７は新潟県付近の地形上に各時刻の気象画像を載せた画像である。その平面（上）と断面で、↑は視線方向を示す。
図７３は８時３０分の表示で平面は明度値８０以下を削除した画像である。削除されたＡ地域の断面の明度値の最低は７１、最高は８５で、その明度差は１４である。ここで最高は最低値周辺の最高値を北南西東求めて平均した。その値は北８９、南７３、西８９、東８８である。断面は東西、南北の２向について表示している。【０１１４】  同様に図７４は９時００分の表示で平面は明度値８１以下を削除した画像である。明度の低い地域がＡ、Ｂ、Ｃの３ヶ所ある。南北断面１と東西断面２、３で、その地域が表示されている。Ａ地域の断面の明度値の最低は６５、最高は８３で、その明度差は１８である。【０１１５】  図７５は１０時００分の表示で平面は明度値明度値８０以下を削除した画像である。Ａ地域の断面の明度値の最低は６１、最高は９３で、その明度差は３２である。この時刻でＡ地域の明度差が最大になる。南北断面で約７３ｋｍの範囲で、明度値の急激な変化がこの地域に生じている。この時刻の１３分後頃に、Ｂ地域の海岸で地震が発生している。【０１１６】  図７６は１０時３０分の表示で平面は明度値８０以下を削除した画像である。Ｃ地域の断面の明度値の最低は６０、最高は８１で、その明度差は２１である。【０１１７】  図７７は１１時００分の表示で平面は明度値８０以下を削除した画像である。Ｂ地域の断面の明度値の最低は６３、最高は９４で、その明度差は３１である。【０１１８】  図７８は１０時００分の３Ｄ表示でＡ地域を表示している。図上は東から西方向の視線である。図下は北西から南東方向の視線である。Ａ地域は地形の平野上にある。【０１１９】  図７９〜８５は１０時００分の平面の表示で明度値６２〜８０以下を削除した画像である。明度の低い地域が新潟県の内陸部から震源地の北西海岸へと広がっている。【０１２０】新潟県における余震
２００７年０７月１６日１５時３７分頃、新潟県で余震が発生した。震源地は上中越沖（北緯３７．５°、東経１３８．７°、新潟の南西６０ｋｍ付近）で震源の深さは約１０ｋｍ、地震の規模（マグニチュード）は５．６と推定されている。長岡市小島谷で最大震度は６弱である。気象庁発表の各地の震度に関する情報、地震情報を図８６、８７に示す。震源地の上中越沖を基点に、震度１、２および３が広範囲に分布している。【０１２１】  １６時前後の気象画像について、新潟県付近の画像を詳細に検討する。各時刻の気象画像を図８８に示す。図８６、８７の震度分布付近で、画像に明度の変化が生じている。【０１２２】  図８９は新潟県付近を拡大した画像で、右は経線等を削除後の画像である。気象画像と地形との対応位置として、２個の基点を←で示す。【０１２３】  図９０〜９２は新潟県付近の地形上に各時刻の気象画像を載せた画像である。その平面（上）と断面で、↑は視線方向を示す。
図９０は１５時００分の表示で平面は明度値１２６以下を削除した画像である。削除されたＢ地域の断面の明度値の最低は７３、最高は１４４で、その明度差は７１である。この時刻でＡ地域の明度差が最大になる。断面は東西、南北の２向について表示している。南北断面で約８２ｋｍの範囲で、明度値の急激な変化がこの地域に生じている。この時刻の３７分後頃に、Ｃ地域の海岸で地震が発生している。【０１２４】  同様に図９１は１６時００分の表示で平面は明度値１２８以下を削除した画像である。Ｂ地域の断面の明度値の最低は１０１、最高は１５３で、その明度差は５２である。【０１２５】  図９２は１６時３０分の表示で平面は明度値１３５以下を削除した画像である。Ｂ地域の断面の明度値の最低は９２、最高は１４９で、その明度差は５７である。【０１２６】  図９３は１５時００分の３Ｄ表示でＡ地域を表示している。図上は南西から北東方向の視線である。図下は北西から南東方向の視線である。Ｂ地域は地形の海岸上にある。【０１２７】  図９４〜９９は１５時００分の平面の表示で明度値７６〜１２６以下を削除した画像である。明度の低い地域が新潟県の南西海岸から震源地の北東へと広がっている。【０１２８】  以上から判る通り、新潟県中越沖の地震では気象条件によっては、その地域の震度が大きくても明度差が大、小異なる。【０１２９】  以上のことから明らかなように、地震時の気象画像データを３次元表示し、その明度差から震源地を特定することができた。
すなわち、具体例から、地震時気象画像解析結果から、地震発生約２時間前から異常画像の明度解析が可能である。また各地震時画像解析から以下のことが言える。
（１）能登半島沖の地震では震源地付近の明度差が大で、その範囲が広い地域の震度は大きい。
（２）茨城南部の地震ではマグニチュード４、震度３程度の地域の明度解析が可能である。
（３）北海道の地震では活断層分布付近で明度が低く、その明度差が大である。
（４）新潟中越沖の地震では気象条件によってはその地域の震度が大きくても明度差が大、小異なる。【０１３０】  なお、以上詳細に説明した方法は、上記手順に応じたプログラムとして形成し、ＣＰＵ等のコンピュータで実行するように構成することも可能である。
また、このようなプログラムは、半導体メモリ、磁気ディスク、光ディスク、フロッピー（登録商標）ディスク等の記録媒体、この記録媒体をセットしたコンピュータによりアクセスし上記プログラムを実行するように構成可能である。
【図面の簡単な説明】
【０１３１】
【図１】本発明の実施形態に係る画像解析装置（気象異常検出装置）の構成例を示すブロック図である。
【図２】地形データのメッシュにその番号と節点番号、さらに座標値を付加した例を示す図であって、（Ａ）はメッシュの削除前を示す図であり、（Ｂ）はメッシュの削除後を示す図である。
【図３】海域と湖の削除前後の九州の地形画像を示す図であって、（Ａ）が海域と湖の削除前の画像を示す図であり、（Ｂ）が海域と湖の削除後の画像を示す図である。
【図４】海域と湖の削除後の九州の３次元表示画像を示す図である。
【図５】経線、緯線の削除前後の赤外線画像を示す図であって、（Ａ）が経線、緯線の削除前の画像を示す図であり、（Ｂ）が経線、緯線の削除後の画像を示す図である。
【図６】線形補間前後の赤外線画像表示を示す図であって、（Ａ）が画像のピクセル数の少ない線形補間前（７６行６２列）の赤外線画像を示す図であり、（Ｂ）が線形補間後（３倍，６０％）の赤外線画像を示した図である。
【図７】画像データ群を示すグラフであり、（Ａ）はフーリエ級数に回帰させるデータ群のグラフであり、（Ｂ）は（Ａ）に示すデータ群をフーリエ級数に回帰させて得られる近似曲線のグラフである。
【図８】気象画像の基点として日本地形の区域での４交点を示す図である。
【図９】気象庁発表の各地の震度に関する情報、地震情報を示す図である。
【図１０】気象庁発表の各地の震度に関する情報、地震情報を示す図である。
【図１１】石川県周辺の活断層、地震分布を示す図である。
【図１２】各時刻の気象画像を示す図である。
【図１３】石川県付近を拡大した画像で、右は経線等を削除後の画像を示す図である。
【図１４】石川県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図１５】石川県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図１６】石川県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図１７】石川県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図１８】石川県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図１９】石川県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図２０】石川県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図２１】石川県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図２２】９時００分の平面の表示で明度値１３３〜１８０以下を削除した画像を示す図である。
【図２３】９時００分の平面の表示で明度値１３３〜１８０以下を削除した画像を示す図である。
【図２４】９時００分の平面の表示で明度値１３３〜１８０以下を削除した画像を示す図である。
【図２５】９時００分の平面の表示で明度値１３３〜１８０以下を削除した画像を示す図である。
【図２６】９時００分の平面の表示で明度値１３３〜１８０以下を削除した画像を示す図である。
【図２７】９時００分の平面の表示で明度値１３３〜１８０以下を削除した画像を示す図である。
【図２８】気象庁発表の各地の震度に関する情報、地震情報を示す図である。
【図２９】茨城県周辺の活断層、地震分布を示す図である。
【図３０】各時刻の気象画像を示す図である。
【図３１】茨城県付近を拡大した画像で、右は経線等を削除後の画像を示す図である。
【図３２】茨城県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図３３】茨城県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図３４】茨城県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図３５】茨城県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図３６】茨城県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図３７】茨城県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図３８】１３時３０分の３Ｄ表示でＡ地域を表示した図である。
【図３９】１３時３０分の平面の表示で明度値２２〜５１以下を削除した画像を示す図である。
【図４０】１３時３０分の平面の表示で明度値２２〜５１以下を削除した画像を示す図である。
【図４１】１３時３０分の平面の表示で明度値２２〜５１以下を削除した画像を示す図である。
【図４２】１３時３０分の平面の表示で明度値２２〜５１以下を削除した画像を示す図である。
【図４３】１３時３０分の平面の表示で明度値２２〜５１以下を削除した画像を示す図である。
【図４４】１３時３０分の平面の表示で明度値２２〜５１以下を削除した画像を示す図である。
【図４５】気象庁発表の各地の震度に関する情報、地震情報を示す図である。
【図４６】気象庁発表の各地の震度に関する情報、地震情報を示す図である。
【図４７】北海道周辺の活断層、地震分布を示す図である。
【図４８】各時刻の気象画像を示す図である。
【図４９】北海道付近を拡大した画像で、右は経線等を削除後の画像を示す図である。
【図５０】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図５１】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図５２】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図５３】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図５４】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図５５】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図５６】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図５７】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図５８】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図５９】北海道の地形上に各時刻の気象画像を載せた画像を示す図である。
【図６０】１３時００分の３Ｄ表示でＡ、Ｂ、Ｃ地域を表示している図である。
【図６１】１３時００分の平面の表示で明度値４６〜５９以下を削除した画像を示す図である。
【図６２】１３時００分の平面の表示で明度値４６〜５９以下を削除した画像を示す図である。
【図６３】１３時００分の平面の表示で明度値４６〜５９以下を削除した画像を示す図である。
【図６４】１３時００分の平面の表示で明度値４６〜５９以下を削除した画像を示す図である。
【図６５】１３時００分の平面の表示で明度値４６〜５９以下を削除した画像を示す図である。
【図６６】１３時００分の平面の表示で明度値４６〜５９以下を削除した画像を示す図である。
【図６７】１３時００分の平面の表示で明度値４６〜５９以下を削除した画像を示す図である。
【図６８】気象庁発表の各地の震度に関する情報、地震情報を示す図である。
【図６９】気象庁発表の各地の震度に関する情報、地震情報を示す図である。
【図７０】新潟県周辺の活断層、地震分布を示す図である。
【図７１】各時刻の気象画像を示す図である。
【図７２】新潟県付近を拡大した画像で、右は経線等を削除後の画像を示す図である。
【図７３】新潟県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図７４】新潟県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図７５】新潟県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図７６】新潟県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図７７】新潟県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図７８】１０時００分の３Ｄ表示でＡ地域を表示している図である。
【図７９】１０時００分の平面の表示で明度値６２〜８０以下を削除した画像を示す図である。
【図８０】１０時００分の平面の表示で明度値６２〜８０以下を削除した画像を示す図である。
【図８１】１０時００分の平面の表示で明度値６２〜８０以下を削除した画像を示す図である。
【図８２】１０時００分の平面の表示で明度値６２〜８０以下を削除した画像を示す図である。
【図８３】１０時００分の平面の表示で明度値６２〜８０以下を削除した画像を示す図である。
【図８４】１０時００分の平面の表示で明度値６２〜８０以下を削除した画像を示す図である。
【図８５】１０時００分の平面の表示で明度値６２〜８０以下を削除した画像を示す図である。
【図８６】気象庁発表の各地の震度に関する情報、地震情報を示す図である。
【図８７】気象庁発表の各地の震度に関する情報、地震情報を示す図である。
【図８８】各時刻の気象画像を示す図である。
【図８９】新潟県付近を拡大した画像で、右は経線等を削除後の画像を示す図である。
【図９０】新潟県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図９１】新潟県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図９２】新潟県付近の地形上に各時刻の気象画像を載せた画像を示す図である。
【図９３】１５時００分の３Ｄ表示でＡ地域を表示している図である。
【図９４】１５時００分の平面の表示で明度値７６〜１２６以下を削除した画像を示す図である。
【図９５】１５時００分の平面の表示で明度値７６〜１２６以下を削除した画像を示す図である。
【図９６】１５時００分の平面の表示で明度値７６〜１２６以下を削除した画像を示す図である。
【図９７】１５時００分の平面の表示で明度値７６〜１２６以下を削除した画像を示す図である。
【図９８】１５時００分の平面の表示で明度値７６〜１２６以下を削除した画像を示す図である。
【図９９】１５時００分の平面の表示で明度値７６〜１２６以下を削除した画像を示す図である。
【符号の説明】
【０１３２】  １０・・・画像解析装置（気象異常検出装置）、１１・・・地形データ記憶部、１２・・・気象画像データ提供部、１３・・・気象画像データ取得部、１４・・・記憶部、１５・・・操作部、１６・・・表示部、１７・・・処理部。
【図１】
2008209417.tif 000007
【図２】
2008209417.tif 000008
【図３】
2008209417.tif 000009
【図４】
2008209417.tif 000010
【図５】
2008209417.tif 000011
【図６】
2008209417.tif 000012
【図７】
2008209417.tif 000013
【図８】
2008209417.tif 000014
【図９】
2008209417.tif 000015
【図１０】
2008209417.tif 000016
【図１１】
2008209417.tif 000017
【図１２】
2008209417.tif 000018
【図１３】
2008209417.tif 000019
【図１４】
2008209417.tif 000020
【図１５】
2008209417.tif 000021
【図１６】
2008209417.tif 000022
【図１７】
2008209417.tif 000023
【図１８】
2008209417.tif 000024
【図１９】
2008209417.tif 000025
【図２０】
2008209417.tif 000026
【図２１】
2008209417.tif 000027
【図２２】
2008209417.tif 000028
【図２３】
2008209417.tif 000029
【図２４】
2008209417.tif 000030
【図２５】
2008209417.tif 000031
【図２６】
2008209417.tif 000032
【図２７】
2008209417.tif 000033
【図２８】
2008209417.tif 000034
【図２９】
2008209417.tif 000035
【図３０】
2008209417.tif 000036
【図３１】
2008209417.tif 000037
【図３２】
2008209417.tif 000038
【図３３】
2008209417.tif 000039
【図３４】
2008209417.tif 000040
【図３５】
2008209417.tif 000041
【図３６】
2008209417.tif 000042
【図３７】
2008209417.tif 000043
【図３８】
2008209417.tif 000044
【図３９】
2008209417.tif 000045
【図４０】
2008209417.tif 000046
【図４１】
2008209417.tif 000047
【図４２】
2008209417.tif 000048
【図４３】
2008209417.tif 000049
【図４４】
2008209417.tif 000050
【図４５】
2008209417.tif 000051
【図４６】
2008209417.tif 000052
【図４７】
2008209417.tif 000053
【図４８】
2008209417.tif 000054
【図４９】
2008209417.tif 000055
【図５０】
2008209417.tif 000056
【図５１】
2008209417.tif 000057
【図５２】
2008209417.tif 000058
【図５３】
2008209417.tif 000059
【図５４】
2008209417.tif 000060
【図５５】
2008209417.tif 000061
【図５６】
2008209417.tif 000062
【図５７】
2008209417.tif 000063
【図５８】
2008209417.tif 000064
【図５９】
2008209417.tif 000065
【図６０】
2008209417.tif 000066
【図６１】
2008209417.tif 000067
【図６２】
2008209417.tif 000068
【図６３】
2008209417.tif 000069
【図６４】
2008209417.tif 000070
【図６５】
2008209417.tif 000071
【図６６】
2008209417.tif 000072
【図６７】
2008209417.tif 000073
【図６８】
2008209417.tif 000074
【図６９】
2008209417.tif 000075
【図７０】
2008209417.tif 000076
【図７１】
2008209417.tif 000077
【図７２】
2008209417.tif 000078
【図７３】
2008209417.tif 000079
【図７４】
2008209417.tif 000080
【図７５】
2008209417.tif 000081
【図７６】
2008209417.tif 000082
【図７７】
2008209417.tif 000083
【図７８】
2008209417.tif 000084
【図７９】
2008209417.tif 000085
【図８０】
2008209417.tif 000086
【図８１】
2008209417.tif 000087
【図８２】
2008209417.tif 000088
【図８３】
2008209417.tif 000089
【図８４】
2008209417.tif 000090
【図８５】
2008209417.tif 000091
【図８６】
2008209417.tif 000092
【図８７】
2008209417.tif 000093
【図８８】
2008209417.tif 000094
【図８９】
2008209417.tif 000095
【図９０】
2008209417.tif 000096
【図９１】
2008209417.tif 000097
【図９２】
2008209417.tif 000098
【図９３】
2008209417.tif 000099
【図９４】
2008209417.tif 000100
【図９５】
2008209417.tif 000101
【図９６】
2008209417.tif 000102
【図９７】
2008209417.tif 000103
【図９８】
2008209417.tif 000104
【図９９】
2008209417.tif 000105