87

日本惑星科学会誌Vol.13.No.2,2004

 
る[9]. しかし, 非専門家によるクレータ抽出には上記の判断基準の問題があり, 注意が必要である. このため, 自動クレータ抽出アルゴリズムを利用するとい近年, 世界各国で月や火星, 小惑星に向けて探査機う選択も考える必要がある.

が打ち上げられ, また SELENE(SELenological and E

自動抽出したクレータ情報の利用分野として, クレーNgineering Explorer)[1]等の大型プロジェクトも計画

タ年代学以外にも Pieters[10]のようにクレータの中央,されている. 現在既に処理しきれないほどの画像デー

壁, リム等のスペクトル情報から, その場所の地殻のタが蓄積され, 今後も増え続けることが必至である.

鉛直構造を推定することがあげられる. 我々は, 自動この莫大なデータからどのように欲しい情報を抽出す抽出したクレータに対し. スペクトル情報を用いて地ればよいのだろうか？ここではクレータに焦点を絞る.

殻の鉛直構造を自動解析するシステム開発を視野に入クレータは, 隕石の衝突や火山によって形成されたれている[11]. その他にも小さなクレータが自動抽出

円形状の窪み地形を言う. 月においては, そのほとん可能であれば, 2次クレータ研究の効率化が図れる.

どが隕石起源であるといわれている. クレータの情報

また, 天体表面への隕石フラックスの直径分布をさらを利用した分野として, クレータ年代学[2-5]がある.に小さな領域まで広げて推測できる可能性がある.

クレータ年代学とは, クレータの面積密度からその惑

そこで, 本論文では目視によるクレータ抽出結果が星の表面が形成された年代を推定するもので, 目視判人によってどれだけ違うのか, 複数の人の同じ画像に読によるクレータ抽出結果を用いて発展してきた分野対する結果を比較し, 議論する. また同じ画像に, 著である. 現在クレータの自動,または半自動抽出手法

者が提案している自動クレータ抽出アルゴリズム[12]の研究[6-8]が行われているが, 未だ実用化されていを適用した結果との相違点を明らかにする. それらをない. そのため, 現在でもクレータ抽出は手動により踏まえ, 今後増大する月惑星画像に対するクレータ抽行われている.

出の一つの方向を提案したい.

手動クレータ抽出の際に, 画像上の輝度パターンをクレータとして認識する判断基準を明確に定義し, か 

つそれを作業中一定に保つことは困難である. そのため, 従来では米国地質調査所等の限られた専門家によ複数の作業者が, Apollo, Lunar Orbiter, Clementineる結果のみが, クレータ年代学の議論等では用いられ画像について目視によるクレータ抽出を行った. ここてきた. しかし, それでは増大するデータに対応しきではクレータを画像上の円形構造と定義する. 目視判れない. それを受け, 非専門家により, 人海戦術で多読の基準として,くのデータからクレータを抽出しようという動きがあ  

・抽出するクレータの半径は２画素以上複数の作業者によるクレータ手動抽出結果とクレータ自動抽出への展望／沢辺, 松永, 六川88  ! 

"



#


#



#
$




! 


"


#


%&&#

'
"


()
&
*&&&&+'"
*&

%&&##



・疑わしいものは抽出しないという2点を定め, 後はそれぞれの作業者の判断にまかせた. この半径2画素以上とは, 画像上で認識可能

月のほぼ全球を光学センサにより観測した衛星は,な最小円に対応し, 自動アルゴリズム[12]においても現在のところ Clementine 衛星, Lunar Orbiter 衛星で最小円を半径2画素としている. 複数の作業者は, 2名ある. また, Apollo ミッションで月周回中に撮影さの専門家 (月画像を扱っている研究者をここでは専門れた写真がある. そこで, 本研究では Clementine デー家と定義. ただし月惑星地質図作成業務に携わっていタ, Apollo データ, Lunar Orbiter データを使用した.

ない), 4名の非専門家 (月画像は扱っていないが, 画使用データの詳細を表1に示す.

像を用いた研究を行っているリモートセンシング関係の研究者) からなる. ただし, 全員が全ての画像に対 

して作業をしたわけではない. また, 実際の作業は,1968年から1972年にかけて, Apollo 計画として8機の自作の手動クレータカウント補助ソフトを使用し, パ探査機が打ち上げられた. 本研究ではそのうちの Apo

ソコン画面の画像上にマウスを使って円を描くというllo16号機が月周回時に撮影したものを利用した. こ方法を取った. このソフトは画像の拡大・縮小, 画像

れは, Metric Camera によって高度約110km で取得さのスクロール, コントラストの調整 (今回は不使用),れた写真をデジタル化した画像である. 本研究では空仮入力した円の移動機能が付いており, 抽出した円の間分解能の違いによる議論をしないため, Clementine中心位置, 大きさを表データとして保存する.

画像に合わせて100m 空間分解能にした海を撮影した画像を用いた.

  !
"#$

著者が開発した自動クレータ抽出アルゴリズムを目

1966年から1967年にかけて Lunar Orbiter 探査機5機視判読に使用した画像に適用し, 手動抽出結果と比較が打ち上げられた. この5回のミッションによって月する. 本アルゴリズムは, 複数のアプローチを利用す面のおよそ99％が撮影された. 本研究で使用したデーることにより, クレータの大きさ, 形状, 場所 (海やタは, Lunar Orbiter4号機に搭載された Medium Resol高地), カメラ, 太陽高度によらず, 安定して高い精ution カメラによる画像である. この画像を Bowker a度でクレータを抽出することが可能である.

nd Huges[13]により, 空間分解能を約300m にデジタル化したデータを使用した.89

日本惑星科学会誌Vol.13.No.2,2004

  !"#Clementine は, NASA が1994年に打ち上げた月・小惑星探査機である. 搭載された観測機器は, UVVIS(紫外可視・近赤外カメラ) 等である.

UVVIS は,5つの狭帯域バンドと一つの広帯域バンドを持つ. 本
研究では, USGS (アメリカ地質調査所) により測光補正, モザイク処理がすでに施されている UVVIS raw image collection の750nm バンドのデータを利用し!$#た. Clementine データは, 高い太陽高度の下で撮影されており一般にクレータ抽出には不向きである. しかし, 本論文では太陽高度が高い場合でも抽出が可能かどうか検討するため, また今後抽出したクレータに対しスペクトル情報を用いて属性を付加していくため,Clementine データを利用した.
!%#  
Apollo 画像に対して, 専門家(1) (白) と非専門家(1)(黒) が手動抽出を行った結果の例を図1.(a)に示す.
一致率は, 全員が抽出したのべ個数に対し71％である.
また, 専門家同士の違いを見るため, 専門家(1)(白)と(2)(黒)の抽出結果を図1.(b)に示す. 専門家(2)が抽出したクレータ数は, 専門家(1)の319個と比べ, 19 !"# !$%& ' 
()*+,-./0! 1&23*+,./0

7個と少なく, 一致率は65％となった. ただし, 直径2〜 (4) は, 明暗パターンがはっきりとしているクレー

km 以上のクレータは24個中20個一致しており, 不一タを抽出し, あいまいなものや小さいものはクレータ致はより小さいクレータのみであった. 小さなクレーとして抽出しない傾向にあった.

タでは, 専門家同士の抽出結果は, 専門家, 非専門家

また, 6名が全員作業した Apollo 画像について, その抽出結果よりも一致しない. これは, クレータと認の結果の違いを図1. (c) に示す. ここで, 白は6人と識する対象も異なっていたためである. 専門家 (1)も抽出したクレータ, 黒は5人以下が抽出したクレーは, クレータの明暗パターンが非常に薄い (よく見れタを表している. 6人全員が抽出したクレータは, 全

ばクレータに見えるかもしれない) ものをクレータと員が抽出したクレータののべ個数に対し 38％と非常して認識する傾向にあり, 専門家 (2) はほとんどそに低い.

ういったものを抽出していない. また, 専門家 (2)専門家, 非専門家による抽出結果の頻度分布を図3のみが抽出したクレータは非常に小さく, 明暗が比較(a)に, 抽出個数を表4に示す. これより, 直径2km 以的はっきりしているという傾向がある. 非専門家 (1)上までは非常に良く一致しているが, それより小さい複数の作業者によるクレータ手動抽出結果とクレータ自動抽出への展望／沢辺, 松永, 六川   90                        !" #$%&'(クレータでは, かなり個人差が見られることが分かる.

     ')*

+,-た. これは, 用いた地域の地形が複雑で, クレータか地表面の起伏かの判別が非常に難しいためと考えられ  る. この傾向は Clementine の高地画像 (図2 (b)) で同様に, Lunar Orbiter, Clementine 画像について,も見られる. Clementine の海画像 (図2 (a)) は, 太

目視判読結果を図2に, それぞれの抽出結果を表2に記

陽高度が高いこともあり, クレータのリムがはっきりす. 一画像あたりのクレータ数が減ったせいもあり,と見えずクレータの認識が難しい. 特に白く点状に見一致率はさらに低い. 特に Lunar Orbiter の高地画像えるものの取り扱いによって結果に大きく差がでた.

(図2 (d)) では, 全員が一致したのは9個 (2％) となっ

なお, 同一人物による抽出結果から, 同じ特徴をクレー91

日本惑星科学会誌Vol.13.No.2,2004

 !$%&' !" #()*+,-./0! 1&23*+,-./0

タとする判断基準が一定ではない様子も一部見られた.

出に適していたか述べることは難しい. 今後は, 同じ
以上より,クレータであると人の目が判断するにはカメラ, 同じ地域, 異なる太陽高度の画像により検討

個人差があることが, 非常に良く示された.

する必要がある.

  

この節では自動クレータ抽出結果の信頼性及び人にClementine 高地と Lunar Orbiter 高地の画像はほぼよる抽出結果の相関について検討する.

同じ地域を異なる太陽高度で撮影した画像である. 全員が抽出したクレータは, 両画像で同じである. しか

表3に, 手動抽出結果と自動抽出結果の一致率を示し, 画像中央の最大クレータ内部に存在するクレータす. 作業者全員が抽出したクレータを自動アルゴリズは, Clementine 画像では抽出できているが, Lunar Orムも抽出した割合は平均で94％と非常に高い. また,biter 画像では明度が高すぎて抽出できていない. 通半数以上の人が抽出したクレータを正解とすると, 自常太陽高度が低い画像 (ここでは Lunar Orbiter 画像)

動アルゴリズムの抽出率は平均91％となった. ここでの方がクレータ抽出に適しているといわれるが, この半数以上としたのは, 非専門家のみが抽出した不確かように明暗がはっきりしすぎてクレータが埋もれてしな抽出結果に大きく左右されるのを防ぐためである.

まう場合がある. 小さなクレータに関しては, Lunar半数以上とすると, 今回用いた画像については, ほぼOrbiter の空間分解能が低いこともあり, どちらが抽少なくとも一人の専門家が抽出していた.









$
$
"
"
$
$

$
$
"
"
$
$
     
 
!
%%
$



"%


%
"
"


%

45

"

"
"


%
67

-./0

#

$
"


#


$
%
!




#
#
#
#
#
#


複数の作業者によるクレータ手動抽出結果とクレータ自動抽出への展望／沢辺, 松永, 六川



92




 !"#$%&'()*+,-.*+,/ 012!3456'7)*+, / 01-812!3456'9)*+, / 01-81:53456

Apollo 画像に対する専門家による手動抽出と自動今回の手動抽出では画像のコントラスト変換機能を使抽出の累積ヒストグラムを図3(b)に, 差分ヒストグラ用せずに作業を行ったが, 自動抽出のみが抽出したクムを(c)に, 抽出個数を表4に示す. これより, 半径4レータは, コントラストを変えたところ目視による認画素以上の結果は非常に良く一致している. 一方クレー識ができるものがほとんどであった. このことから,タの数と大きさは指数関数的な関係があると一般に言自動抽出のみが抽出をしたクレータは必ずしも誤認識われているが, 手動抽出結果は半径2, 3画素になるとではないことが分かる.

不自然に減少していることが図3(c)から良く分かる.

Apollo 画像に対して, 自動抽出アルゴリズムは直それに対し, 自動抽出ではプロットが左肩上がりとな

径1.6km より大きいクレータについては専門家によるり, 非常に小さいものも良く抽出できている. このこ抽出と同程度抽出でき, 0.8〜1.6km のクレータは非とは, 自動抽出によって, 従来手動クレータ抽出の対専門家が抽出した程度には抽出可能であり, さらに小象とならない小さなクレータも抽出できる可能性を示

さな0.8km 以下のクレータでは専門家よりも多く抽出唆している.

可能であった. なお, 手動抽出の際, コントラスト調

また表4から, 直径16から32画素のクレータでは,専門家による手動抽出個数に比べ, 自動抽出個数が少ないことが分かる. この大きさで専門家のみが抽出したクレータは, クレータの明暗がはっきりとせず, エッ整の有無が抽出結果に及ぼす影響は, 今回の画像では微少であった.
以上のことから, 自動抽出アルゴリズムの精度は実用に耐えうるものと期待できる.

ジ抽出が困難であるものがほとんどであった. また,自動抽出の個数は非専門家が抽出した個数とほぼ同じである. このことから, 非専門家が抽出したクレータについては自動抽出アルゴリズムにより抽出可能であっ複数の作業者による手動クレータ抽出を同じ画像にたが, 専門家のみが抽出した陰日向のコントラストが対して行ったところ, 作業者全員が同じクレータを抽小さいクレータは抽出できなかったと言える.

出した割合は, 平均で３割であった. また, 専門家同自動抽出のみが抽出したクレータ (表2, 一致人数0士であっても同じクレータを抽出した一致率は６割での列) は, ほとんどが小さなクレータである. また,ある. 一致率の低下は, 直径16ピクセル以下のクレー93

日本惑星科学会誌Vol.13.No.2,2004

タで目立った.

自動抽出アルゴリズムの結果は, 半数以上の作業者が手動抽出した結果に対し, 9割以上の一致率となった.

2129.
[11] 沢辺頼子, et. al., 2003, 第35回リモートセンシ

今後の大量データに対するクレータ抽出には, 客観的な自動クレータ抽出手法の重要性が示された. また自動抽出アルゴリズムは, 専門家による限定された範囲の高品質なクレータ抽出を補完する役割を果たせる可能性を示した.


本研究を進めるにあたり, JAXA,[10] Pieters, C. M., 1991, Geophys. Res. Lett., 18, 11,大阪大学の佐伯和

人様より Apollo 画像をご提供いただきました. また,JAXA の春山純一様, 東大の小林史歩様, 東工大の小林秀樹様, 筑波大の石原光則様にはクレータ抽出の作業にご協力いただきました. 査読, 編集担当の皆様には適切なご指摘, 議論をいただきました. ここに感謝の意を記します. また, 著者の一人 (沢辺) は21世紀COE 機械システムイノベーションによる補助を受けました.

[1] 松永恒雄, et. al., 1999, 日本リモートセンシン
グ学会誌, 19, 5, 152.
[2] Hartmann, W. K., et. al., 1972, Astrophy. SpaceScience, 17, 48.
[3] Neukum, G. and Koenig, B., 1975, The Moon, 12,201.
[4] Lissauer, J. J., et. al., 1988, JGR, 93, B11, 13776.
[5] Hiesinger, H., et. al., 2002, Geophys. Res. Lett., 29, 8, 10.1029.
[6] Honda, R., et. al., 2002, LNAI 2281, 395.
[7] Burl, M. C., et. al., 2001, Proc. on Art. Intell.
Robot. Auto. Space 6th, AM118.
[8] Magee, M., et. al., 2003. LPSC 34th.
[9] http://www.cirklo.com/space/seleneathome.html

ング学会学術講演会論文集.
[12] 沢辺頼子, et. al., 2004, リモートセンシング学
会誌 (投稿中)[13] Gilles, J., et. al., 1999, LPSC 30th, 1770.