計測自動制御学会東北支部第 289 回研究集会 (2014.6.24)
資料番号 289-6

夜間の熱赤外画像のための GCP としての人工ホットスポットの抽出
Extraction of the artificial hotspot as the GCP for nighttime thermal infrared image○春木裕真, 丹波澄雄
○Yuma Haruki, Sumio Tanba
弘前大学
Hirosaki University
キーワード : 人工ホットスポット(artficial hotspot), 熱赤外画像(thermal infrared image), 輝度(luminance), 閾値処理(threshold processing)
連絡先 : 〒036-8561 青森県弘前市文京町 3 弘前大学大学院 理工学研究科 理工学専攻 電子情報工学コース 丹波研究室
春木裕真, Tel : , Email : h14ms514@stu.hirosaki-u.ac.jp―――――――――――――――――――――――――――――――――――――――――――
1.序論
は一般に温度が低いので検出可能だが、霧などの人工衛星画像の精密幾何学補正にはシステム補低高度の雲は地表との温度差が小さいので検出が正だけで十分な精度が得られない場合がある。そ困難となる。また、陸域と水域の境界での温度差はの場合、地上基準点(GCP)と呼ばれる地図上の位季節、地域、天候によって多様に変化するため、温置が確定している地点と画像上の対応点の情報の度情報だけで決定することは困難となる。以上より、組みが必要となる。GCP は一般的に半島や岬の先夜間画像の GCP 取得精度は低くなってしまう。

端など視覚的に特徴のある地形が選ばれるが、衛研究室では中間赤外域(3.7μm)と熱赤外域(10-星に搭載されているセンサの地上分解能により選12μm)の観測波長帯を有しているセンサによる夜択される対象物が異なる。高分解能センサの場合

間画像のための GCP として、製鉄所や、セメント工は 人 工構 造物 も よく 利用 さ れ る。 NOAA 衛星 の場といった人工ホットスポット(人工熱源)の利用を提AVHRR セ ン サ の 地 上 分 解 能 は 衛 星 直 下 で 約案している。

1.1km であるので、人工構造物はほとんど利用する本研究では、日本と韓国全域の人工ホットスポットことができない。

候補地点の熱赤外画像での現れ方を調べ、夜間可視域や熱赤外域の波長帯を用いる場合、日中画像の GCP データとして整備し、利用するためのも夜間も晴天でなければ地表の観測が難しいため、 方法を検討した。
GCP の取得は困難になる。日中画像の場合は、雲域をほぼ確実に把握でき、近赤外域を用いることで2.使用した画像データ陸域と水域の境界も明瞭に認識できる。そのため、東海大学情報技術センターで受信・蓄積される日中であれば良質な GCP が取得可能となる。しかNOAA 衛星の AVHRR 画像データを、研究者有志し、日中に比べて夜間画像の場合、熱赤外域のセ

た ち が 開 発 、 公 開 し て い る PaNDA(Package forンサの放射温度の情報のみ使用可能となる。雲域

Noaa Data Analysis)を使用して、日本列島及び、

表 1 : NOAA 衛星号数と撮影月内訳
NOAA
15号
16号
17号
18号
19号
計

4月

5月
0
10
1
1
30
50
66
176月

7月
000
112

8月
0000
000
40
43112013年
2014年
9月 10月
11月
12月 1月 2月 3月 4月
000000000383834400000000
1488115671514911457212302030121518計
04306267172

朝鮮半島が比較的晴れていて、なおかつノイズが輝度で比較すると、ch.3B は 20000 倍に増加するないものを 2013 年 4 月～2014 年 4 月の夜間画像

が、ch.4、ch.5 は 50 倍程度の増加であり、感度が 2から 172 枚選択した。15 号と 17 号はノイズが酷い桁以上異なることがわかる。人工ホットスポットは周ものが多かったため、選択されることはなかった。

囲に比べ温度が高いため、ch.3B と ch.4 の輝度温

AVHRR(改良型高分解能放射計)は衛星直下で度差画像を作成すると差が大きく現れるので検出の空間分解能が約 1.1km、観測幅が約 2800km をすることができる。

有するセンサである。AVHRR センサの ch3 は日中

AVHRR セ ン サ の 地 上 分 解 能 は 衛 星 直 下 では ch.3A を、夜間では ch.3B を使用するようになっ1.1km 四方であり、通常の人工構造物はこの分解ている。日中では全ての波長帯が使用可能だが、能以下のサイズとなる。さらに高温部分(ホットスポッ

夜間の使用可能波長帯は ch.3B の中間赤外域とト)は通常その構造物の一部分であるため、AVHRRch.4、 ch.5 の熱赤外域のみである。表２は NOAA-の 1 画素は非常に占有面積の小さい高温部分とそ

15 号以降の AVHRR センサ諸元である。の他の低温部分の混合画素(ミクセル)となる。

表 2:AVHRR 諸元
Ch
観測波長帯(μm)
1 0.58〜0.68(可視)2 0.725〜1.0(近赤外)3A 1.58〜1.64(中間赤外)3B 3.55〜3.93(中間赤外)
4 10.3〜11.3(熱赤外)5 11.5〜12.5(熱赤外)

観測項目
雲、氷、雪
陸域植生
表面温度、雲表面温度、雲表面温度、雲表面温度、雲４．人工ホットスポットの出現特徴
過去の研究より、NOAA 衛星画像に現れる人工ホットスポット候補地点が 39 地点存在し、地図と照らしあわせた結果、人工ホットスポットとして現れる熱源を有しているのは製鉄所やセメント工場であることがわかった。
選択した 172 枚の夜間衛星画像では陸域画像と3.人工ホットスポット検出原理

海岸線画像のズレが生じているものがほとんどであ

物体は自身の温度に応じた輝度を波長毎に放射る。このままでは切り出しデータ(ダンプデータ)を作しており、その関係はプランク関数の黒体放射関数成する際、安定して人工ホットスポットの検出ができ(1)式で記述される。

ないため、目視で陸域と海岸線を可能な限り一致さB λ=2hc21･5(1)hcλexp()−1kλT

ここで、Bλ は黒体の分光放射輝度[W･m-2 ･sr-1 ･-1せた。その後、ch.3B と ch.4 の輝度温度差画像で39 地点の各座標を中心とした 21 画素×21 画素のダンプデータを作成し、人工ホットスポットの分布状況を観察した。その結果、人工ホットスポットの大きμm ]、 T は黒体の絶対温度[K]、λ は波長[μm]さは各工場で異なり、熱源の分布が時間と共に変、c は光速度(2.998×108[m・s-1])、h はプランク定数化する可能性がある。製鉄所は明るいが、人工ホッ-34(6.626×10 [J ・ s]) 、 k は ボ ル ツ マ ン 定 数

トスポットとなる工場が比較的集中しており、近傍に(1.380×10-23[J・K-1])である。

他の人工ホットスポットが存在することが多い。セメ

プランク関数において NOAA 衛星の AVHRR セント工場は暗いが 1、2 画素程度に収まり、さらに周ンサの ch.3B、ch.4、ch.5 を 300K と 1300K の時の囲に他の人工ホットスポットは殆ど存在しないことがわかった。これは日本、韓国共に似たような結果になった。さらに、同一ダンプ地点内に複数の人工ホットスポットがいくつか存在したため、それらから最も高い輝度温度を示した工場を使用することにした 。
その結果、人工ホットスポットは 39 地点から 35 地

1882 1883 1884 1885 1886 1887
1209
1210
121112121213

439 2225 3064 3064
1831 2260 886
2260 765886

図 2 : 閾値以上の数値例

点となった。

5.人工ホットスポット領域特定
ダンプデータを調べることで製鉄所やセメント工場の特徴を得られることができたが、それだけではGCP の分類には欠ける。各人工ホットスポットの明るさや広さをより詳細に調べるため、図 1 の手順で処理を実行した。 閾値は(標準偏差×3＋平均)のことである

図 3 : 領域標準偏差平均と最大輝度温度差による区分図

図 4:領域標準偏差正規分布

図 1:閾値処理のフローチャート
図 5:最大輝度温度正規分布

出力された閾値以上の数値データから人工ホットスポットの領域を求める。図 2 は閾値以上の数値デはそれぞれの正規分布により求めることができる。

ータの例であるが、このデータから領域を求める方

領域標準偏差正規分布グラフは頂点の左側が狭法はピクセル方向(横)とライン方向(縦)の標準偏差い、右側が広いとみなす。最大輝度温度正規分布を求め、2 つの標準偏差からから人工ホットスポット

グラフは頂点の左側が暗い、右側が明るいとみな領域の標準偏差を求めることができる。

す。

次に、各人工ホットスポットの明暗・広狭区分を行うために領域標準偏差の平均と最大輝度温度の平6.人工ホットスポット検出結果

均を用いてグラフを作成する。図 3 が作成されたグ人工ホットスポットの明暗・広狭区分はできたが、ラフであり、領域標準偏差と最大輝度温度の区分

閾値処理の際、人工ホットスポット以外の地点が領域に含まれることが幾度なくあった。Google Map や上に雲がかかっていても識別できるように、雲を除国土地理院のデータと照らし合わせることでその地去するアルゴリズムの作成を予定している。

点に何が在るか、またどのようなものが領域に含まれるかを調査した。
・造船所、発電所、製油所、空港、自衛隊基地
・製粉工場、ゴム工場、精米工場、ガス工場
・石油化学工場、ステンレス工場
・電気化学工場、機電工場、石灰採石場
・水田、海上、山
その結果、以上のものが領域に含まれることがわかった。この中で検出頻度の高いのが造船所、発電所、電気系、化学系の工場であった。特に発電所の輝度温度は高く、セメント工場の付近にある場
合、セメント工場の輝度温度を上回る可能性もある。
これらの中で異分子なのが水田、海上、山頂の 3つだが、水田に関しては検出したのが夏や秋だったので藁焼きをしていたと考えられる。海上に関しても、検出されたのが人工ホットスポット付近だったため資材を運ぶための大型の船などが検出されたと考えられる。
さらに人工ホットスポットの位置が内陸か海沿いかによる領域標準偏差の違いを検証したが、先に述べた人工ホットスポット以外の地点が存在するか否かで領域標準偏差の値が大きく変わるためこの検証は意味をなさなかった。しかし、人工ホットスポッ
ト周辺に水域がなく、工場などの閾値処理で検出されるような建造物もない場合、検出されるのは人工

表 3:人工ホットスポット分類結果
番号1234567891011121314151617181920212223242526272829303132333435

工場名
太平洋セメント 上磯工場
太平洋セメント 熊谷工場
太平洋セメント 埼玉工場
太平洋セメント 藤原工場
太平洋セメント 大分工場
住友大阪セメント 栃木工場
住友大阪セメント 赤穂工場
宇部興産 宇部セメント工場
宇部興産 伊佐セメント工場
宇部興産 苅田セメント工場
トクヤマ
双龍セメント 東海工場
双龍セメント工業 蔚山工場
東洋セメント 釜山工場
東洋セメント 三陟工場
ラファージハルラセメント港湾工場
太平洋金属 八戸工場
新日鉄八幡製鐵所
新日鉄室蘭製鐵所
JFEスチール 東日本製鉄所京浜地区
住友金属工業 鹿島工場
JFEスチール 東日本製鉄所千葉地区
新日本製鐵 君津工場
新日本製鐵 名古屋工場
住友金属工業 和歌山工場
JFEスチール 西日本製鉄所倉敷地区
JFEスチール 西日本製鉄所福山地区
新日本製鐵 大分工場神戸製鋼 神戸製鉄所
新日本製鐵 広畑製鉄所
日本精鉱 呉製鉄所
住友金属鉱山株式会社 播磨事業所
POSCO 光陽製鉄所
現代製鉄 唐津製鉄所
現代製鉄 本社

位置
海岸
内陸
内陸
内陸
海岸
内陸
海岸
海岸
海岸
内陸
海岸
内陸
内陸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸
海岸

LA HA
狭 暗
広 暗
狭 暗
狭 暗
狭 暗狭 暗狭 暗広 暗狭 暗狭 暗狭 暗広 明狭 暗狭 暗広 暗狭 暗狭 明
広 暗広 明
広 明
狭 明
広 明
広 明
狭 明
狭 明
広 明
狭 明
狭 明
狭 暗
広 暗
狭 明
狭 明
広 明
広 明
広 暗

ホットスポットのみとなるため、そのような人工ホットスポットは非常に良質な GCP と呼べる。表 3 は人工ホットスポットの分類結果となる。ここで、LA は領参考文献

域標準偏差の平均、HA は最大輝度温度の平均を[1] 佐々木靖彦「夜間の衛星画像データ処理の占めす。

ための人工ホットスポットの利用可能性に関する研究」2013 年度学士論文

7.まとめと今後の方針[2] 佐藤佑樹「人工ホットスポットを用いた衛NOAA 衛星の AVHRR 画像内で利用可能な人工星画像の精密幾何補正」2011 年度修士論文

ホットスポットを対象とした GCP のデータセットが作[3] 日本測量協会 リモートセンシング読本成できた。これにより、明るく広い GCP を用いて 2,3画素程度の誤差でズレを修正し、暗く狭い GCP を用いて 1 画素程度までズレを修正する二段階補正が可能となった。
今後は、ノイズが多い NOAA-15 号、17 号での人工ホットスポットの高精度検出や人工ホットスポット