REDD+

Reducing Emission from Deforestation
and Forest Degradation-plus

平成2５年度基礎講習

第７章
リモートセンシングを用いたリモートセンシングを用いた森林面積の把握手法

一般社団法人 日本森林技術協会
鈴木 圭
1

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

リモートセンシングデータを用いた土地利被覆区分の時系列解析の流れ
前処理

単時点の分類

時系列解析

・幾何補正
・大気補正
・地形補正
・オブジェクト分類
・ピクセル分類・目視判読・分類結果の差によるもの
・２時点の画像間の変化を直接検出する方法
2

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

変化抽出技術の分類[１] 各時点でそれぞれ森林タイプ分類図を作成して差分をとる方法
オブジェクトベース分類

教師付き分類*ピクセルベース分類

教師なし分類*

自動分類

手動分類

ルールベース分類

目視判読*：ピクセルベース分類でも[２] ２時点の画像間の変化を直接検出する方法
画像間の差バンド間の差画像間の積指数画像の差

ＮＤＶＩ

自動検出
主成分分析

手動検出

・・・・・・・・
目視判読
3REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

変化抽出技術の特徴（１）長所

短所[１] 各時点でそれぞれ森 作成した分類図から炭素林タイプ分類図を作成し 量への換算が比較的容易て差分をとる方法

分類精度が高くないと、変化がない箇所を誤抽出する可能性がある[２] ２時点の画像間の変 変化のある箇所の検出精化を直接検出する方法 度は比較的高い

炭素量への換算が困難

4

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

変化抽出技術の特徴（2）『[１]各時点でそれぞれ森林タイプ分類図を作成し差分をとる方法』の各分類手法
長所

自動分類
高分解能の衛星画像の- オブジェクトベース分類 分類に適しているため、詳細な分類図を得られる

短所

分類のためのパラメータ設定が複雑

自動分類- ピクセルベース分類

作業者の技術力によらず、 高分解能の衛星画像の比較的均質な成果を得ら 分類には適さない
れる

手動分類- 目視判読分類精度が高い経験にもとづく技術力が必要であり、また作業量が多くなる

5

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

[１]各時点でそれぞれ森林タイプ分類図を作成し差分をとる方法

オブジェクトベース分類

セグメンテーションと呼ばれる処理により、スペクトル情報や形状情報に基づいたオブジェクト（ピクセルの集合）が生成される

6

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

[１]各時点でそれぞれ森林タイプ分類図を作成し差分をとる方法

ピクセルベース分類

類似したスペクトル情報を持つピクセルをまとめることにより分類する。空間分解能の高くない衛星画像の分類に適す

7

REDD+Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

[１]各時点でそれぞれ森林タイプ分類図を作成し差分をとる方法

目視判読

対象物の持つ『色調』『形状』『大きさ』『きめ』『模様』などを手がかりにした作
業者の判断により分類する

8

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

[１]各時点でそれぞれ森林タイプ分類図を作成し差分をとる方法

３つの分類方法の比較
オブジェクトベース分類

ピクセルベース分類

目視判読9

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus[２] ２時点の画像間の変化を直接検出する方法

バンド間の差分による検出【 バンド７ 散布図 】

・各バンドは地表面の様々な特徴をとらえている。
期末・森林の変化検出にはLANDSATバンド７（短波長赤外線）がよく用いられる。
・一般的に、植林による変化より伐採による変化のほうが検出しやすい。
RESTECホームページより引用

期首10

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

[２] ２時点の画像間の変化を直接検出する方法

バンド間の差分による検出
期首画像
LANDSA 1992年
バンド７差分画像（緑：植林、赤：伐採）

期首画像
LANDSA 2002年

11

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

[２] ２時点の画像間の変化を直接検出する方法

NDVI画像の差分による検出
NDVI：正規化植生指標

NDVI

NIR
NIR

赤バンド（R)

R
R

NIR : 近赤外バンド
R : 赤バンド

近赤外バンド（NIR)

NDVI画像12

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus[２] ２時点の画像間の変化を直接検出する方法NDVI画像の見え方（例）
赤バンド(R）
正規化植生指標NDVI = (NIR – R) / (NIR + R)

近赤外バンド（NIR）

13

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

画像解析ソフトの紹介
オブジェクト分類に対応した主なソフトウェア
ソフトウェア

開発会社

eCognition
Feature Analyst
for ArcGIS
ENVI EX
PicassoTrimbleOverwatch Systems
EXELIS（株）つくばアグリサイエンス

オブジェクト分類に対応したソフトウェアの特徴としては、Segmentation処理を行うことにより画像分類をおこなう機能を持っている。
Segmentation？
14

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

画像解析ソフトの紹介
eCognitionの画像分割の仕組み• Segmentation：隣接画素との均質性を基に画素をグループ化– 色調や形から算出した異質性（Heterogeneity）の指標を用いて、分割是非を判断する。
異質性(ｆ )は色調（f）と形状（）の合計（2.1式）。

w * hcolor (1 w) * hshape(2.1)

hcolorとhshapeはwによって重み付けされ合計が１となる。

hshapewcmpct * hcmpct (1 - wcmpct) * hsmooth(2.2)

形状（）はとの2因子の合計（2.1式）。
(2.1式)に(2.2式)を代入すると次式を得る。

fw * hcolor (1 w) * {w cmpct * hcmpct (1 - wcmpct) * hsmooth } (2.3)
15

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

画像解析ソフトの紹介
eCognitionの画像分割の仕組み• 異質性（Heterogeneity）の算出– 色調の異質性（hcolor）：輝度値分散の総和とみなす。

hwc *(2.4)

cc

σ：レイヤーCのポリゴン内輝度値の分散
w：そのレイヤーCの重み付けした値

SegmentationではObject１とObject2を統合する場合、それぞれの標準偏差を合計し、統合後のObjectの標準偏差との差をとる。これにレイヤーの重み付けを乗算した値がhcolorである（2.5式）。

hcolorwc(nMerge *Mergec- (nObj1 *Obj1cnObj2 *

Obj 2c)) (2.5)

c16REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

画像解析ソフトの紹介
eCognitionの画像分割の仕組み• 異質性（Heterogeneity）の算出– 形状の異質性（hsmooth）：周長と短辺の比によって形状を評価する。

hlb(2.6)

l：ポリゴンの周長、ｂ：ポリゴンの短辺

Object１とObject2を統合する場合、 (2.6式)から個々に算出したhの値を合計し、統合後のObjectのhとの差をとった値がhsmoothである（2.7式）。

hsmoothlMergenMerge *bMerge

lObj1( nObj1 *
bObj1

lObj2
nObj2 *)bObj2

( 2 .7 )17

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

画像解析ソフトの紹介
eCognitionの画像分割の仕組み• 異質性（Heterogeneity）の算出– 形の異質性（ hcmpct ）：周長と画素数の平方根の比によって形状を評価する。

hln(2.8)

l：ポリゴンの周長、n：ポリゴン内の画素数

Object１とObject2を統合する場合、 (2.8式)から個々に算出したhの値を合計し、統合後のObjectのhとの差をとった値がhcmpctである（2.9式）。

hcmpctlMergenMerge *nMergelObj1(nObj1 *
nObj1

lObj2
nObj2 *)
nObj2

(2.9)
18

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

画像解析ソフトの紹介
eCognitionの画像分割の仕組み• 異質性（Heterogeneity）の算出– （2.5式）、(2.7式)、（2.9式）によって算出されたhcolor、hsmooth、hcmpctの3因子を（2.3式）に代入することで異質性(ｆ )が求められる。

• Segmentation処理では異質性(ｆ )はとして設定する。
– 値を小さくすると異質性の許容範囲が狭まり画像をより細かく分割するにとどまる結果となり、大きくすると粗い分割結果となる。そのほかに先に述べた3つの因子の重み付けや、Segmentation処理時に使用する各レイヤー（例えばRed、Blue、Nirなど）の重み付けによりSegmentation結果は左右される。
19

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

画像解析ソフトの紹介
eCognitionの画像分割の仕組み（異なるScale parameterによる区画線）
SP=50

区画されていない

SP=200

SP=100

まだ細かいが必要な区画はされている

空中写真の目視判読（正）20REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

画像解析ソフトの紹介
eCognitionの画像分割の仕組み（異なる林相における区画線）
SP=80
スギとヒノキヒノキ

異質性が低いのでSP80では区分されない。

スギと広葉樹異質性が高いのでSP80でも区分される。

SP=30スギ
スギ

広葉樹スギ

21

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

画像解析ソフトの紹介
eCognitionの画像分割の仕組み（画像分類結果の比較）

空中写真 – 目視判読

IKONOS - eCognition

平成14年度森林資源モニタリング調査データ地理解析事業（リモートセンシング資源解析事業）報告書より
22

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

画像解析ソフトの紹介
ENVI EXの画像分割の仕組み• 作業の流れ23

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

画像解析ソフトの紹介
ENVI EXの画像分割の仕組み• Segmentation
明るさ、テクスチャ、色などの数値が似ている近隣ピクセルをScale Levelを設定することでグループ化し、画像を分割する。
Scale Level 20

Scale Level 60

境界を作りたいところに線がある

土地利用の境界に線がない箇所が多い
24

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

森林－非森林変化のモニタリング（GOFC-GOLD / COP18）

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/25

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

モニタリングアプローチの選択と実施
2.1.2.4 Selection and implementation of a monitoring approach - deforestation
モニタリングアプローチの選択と実施－森林減少
•••
•
•
•

Step1: Selection of the forest definition
森林の定義の選択
Step2: Designation of forest area for acquiring satellite data
衛星データを必要とする森林域
Step3: Selection of satellite imagery and coverage
衛星画像と適用範囲
Step4: Decisions for sampling versus wall to wall coverage
サンプリングvs全域の決定
Step5: Proccess and analyze the satellite data
衛星データの処理と解析
Step6:Accuracy assessment
精度検証

2.1.2.5 Monitoring of increases in forest area – forestation
森林回復のモニタリング
GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/26

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

Step1: Selection of the forest definition
森林の定義の選択

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/27

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

Step2: Designation of forest area for acquiring satellite data

衛星データを必要とする森林域・ 国土（対象エリア）のすべての森林を含める必要がある。
・ 現存している全森林は，評価期間中すべての森林変化をモニタリングする必要がある。

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/28

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

Step3: Selection of satellite imagery and coverage
衛星画像と適用範囲
・光学センサーの解像度と利用状況

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/29

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

Step3: Selection of satellite imagery and coverage
衛星画像と適用範囲
・中解像度光学
センサの利用現況

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/30

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

【資料】 代表的な光学衛星センサのスペックおよびデータ価格（REDD-plus Cook Book）

REDD-plus Cook Book（森林総合研究所）

31

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

Step4: Decisions for samplingversus wall to wall coverage
サンプリングvs全域の決定

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

32

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

Step5: Proccess and analyze
the satellite data
衛星データの処理と解析

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

33

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

○解像度に応じた適切な区分（解析）手法

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

34

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

○解像度に応じた適切な区分（解析）手法

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

35

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

Step6:Accuracy assessment
精度検証
不確実性評価のキーとなる要因について（森林タイプ面積、炭素蓄積量の推定）
・計測、測定によるエラー（人為的な要因）
・データのキャリブレーション（調整）に起因するエラー
・モデリングによるエラー（現象を説明しきれないモデルを作成してしまう）
・サンプリングによるエラー（不適切なサンプリング間隔の設計）
・不適切な分類項目や定義に起因するエラー（カウントもれ、ダブルカウント）
GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/36

REDD+平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

Step6:Accuracy assessment
精度検証
森林タイプ面積の不確実性について（単時点の地図評価）
・ アクティビティーデータとしてリモセン由来の森林タイプ面積を使用する場合→ 分類（classification）エラーや、判読者によるバイアスが含まれる。
→ パラメータのチューニングによるバイアスの修正（恣意性が含まれる。）→ どの程度の恣意性が含まれるか評価することが重要
・ 高精度なリファレンスデータが入手できる場合→ 統計的にしっかりしたデータ（層化サンプリング、システマティック…）→ オリジナルデータをこれらのデータでキャリブレート（調整）する。
→ IPCC GPG ５章にいくつか方法論のリコメンデーション有り（複数年度の地図の経年変化評価）
・アクティビティーデータとしてリモセン由来の森林タイプ面積を使用する場合→ 土地利用・被覆面積の量的な不確実性の評価ができるような工夫が必要→ 主題となる分類項目（森林タイプ）に着目した精度計測、信頼区間算定→ 既存の統計情報やバイアスを数値化したものでキャリブレーションする等
GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

37

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

Step6:Accuracy assessment
精度検証
衛星画像データに関するエラーの要因、考慮事項について（考慮する事項）
・ 衛星画像データの質、特徴（空間情報、スペクトラル情報、撮影周期、雨季、乾季のようなフェノロジー的な要素 ・・・）
・ 衛星画像入手時には、センサーの処理レベルに応じて地理的な位置精度や画質のチェック・ 異なるセンサー間でのデータの取扱い
・ 幾何補正（地理的な歪みの取除き）、大気補正等
・ 主題図作成の基準（土地分類カテゴリーや最小作図単位など）
・ 画像判読の手順（分類アルゴリズム、目視判読基準など）
・ 地図作成後の処理（後処理・・・ベクタ／ラスタ変換、０値の取扱いなど）
・ 精度評価のためのレファレンスデータの入手（グランドトルース、キャリブ
レーション用データなど）
・ 地図作成時には、主題図作成の統一的かつ透明性のある基準に従って専門家による適切な判読、モニタリングが為されるべき。→ 「Consistent mapping」→ 判読キーなどの文書化が必須（統一性、透明性）→ トレーニングデータの取得方法、検査方法の手順化（統一性、透明性）→ 雲や雲陰などのNo dataの取扱い、処理方法も手順化する必要有り。
38

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

2.1.2.5 Monitoring of increases
in forest area – forestation
森林回復のモニタリング• 森林の回復スピードは遅いので，森林減少よりも見つけにくい。
• スタンダードな手法はない。
• 樹冠疎密度で推定などが考えられる。

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

39

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

森林内での変化モニタリング
2.2.2.1 Direct approach to monitor selective logging
択伐モニタリングの直接的アプローチ• Step1:Define the spatial resolution
利用する空間分解能• Step2:Enhance the image
画像強調の方法• Step3:Select the mapping feature and methods
マッピング（図化）機能と方法の選択
2.2.2.2 Indirect approach to monitor forest degradation
森林劣化モニタリングへの非直接的アプローチ

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/40

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

2.2.2.1 Direct approach to monitor selective logging
択伐モニタリングの直接的アプローチ• 様々な土地被覆の状態で見えるので，択伐地域の抽出は難しい。
• 高分解能衛星の利用が必要

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

41

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

2.2.2.1 Direct approach to monitor selective logging
択伐モニタリングの直接的アプローチ
Step 1. 利用する分解能を決める

Step 2. 画像強調の方法を決める

高分解能（5m以下）
中分解能（10～60m）

大気補正
ヒストグラムマッチング
テクスチャフィルタ
ミクセル分析

Step 3-1. 図化対象を決める

樹冠ギャップの抽出
領域での抽出Step 3-2. 抽出方法を決める

目視判読
自動分類

結果を確認するGOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

42

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

Step1:Define the spatial resolution
利用する空間分解能

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/43

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

Step2:Enhance the image
画像強調の方法

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

44

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

Step3:Select the mapping feature and methods
マッピング（図化）機能と方法の選択
・ 森林劣化の程度は?
・ 解析対象領域の広さは?
・ 解析手法は?

GOFC-GOLDホームページ
http://www.gofcgold.wur.nl/redd/

45

REDD+

平成25年度 基礎講習

Reducing Emission from Deforestation and Forest Degradation-plus

2.2.2.2 Indirect approach to monitor forest degradation
森林劣化モニタリングへの非直接的アプローチ

自然林（Intact Forest）の抽出• Positive approach‐ 画像解析で自然林（IntactIntact ForestForest）を抽出• Negative approach
1. GISデータなどで人為改変地域および周辺地域を除く。
2 .目視判読でFine Shapingする。

GOFC-GOLDホームページ http://www.gofcgold.wur.nl/redd/

46

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

リモートセンシングを用いた森林面積の推定（出典：REDD-plus COOKBOOK）

47

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

COOKBOOKの該当目次

48

REDD+Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

COOKBOOKの該当レシピ

49

REDD+Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

雲の除去（Recipe-T06)
処理前

処理後

雲、ヘイズ（もや）及びそれらの影の抽出

前後の時期の画像から雲のない部分を抽出し、モザイク処理
50

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

季節性の調整（Recipe-T06)
・熱帯落葉樹林では、着葉期と落葉期があることから、同季節の画像を収集することが基本となる。
・同季節の画像の入手が困難な場合は、季節性の影響を画像処理で調整してから分類するなどの処理が必要となる。
処理前

処理後

季節変化の影響の除去の例
51

REDD+

Reducing Emission from Deforestation and Forest Degradation-plus

平成25年度 基礎講習

グランドトゥルース（Recipe-T0８)
グランドトゥルース：リモートセンシングデータと観測対象物（森林）との対応関係を明らかにするために観測・測定・収集した地上の実態に関する情報。
現地調査が難しい場合は、より高分解能の衛星画像や空中写真の判読で代用する場合もある。
グランドトゥルースの位置情報：現地調査で取得する際には、ＧＰＳや既存の地図情報、空中写真の画像判読などを利用する。
グランドトゥルースの抽出方法：抽出方法は単純無作為抽出、層化抽出、系統的抽出など様々、統計的に適切に抽出する必要がある。
サンプル数：教師つき分類の教師データ（分類に用いるデータ）として必要なクラ
スごとのサンプル数は、用いる説明変数の10倍（例えばバンド数）といわれている。また、精度検証に用いるサンプルは、分類に用いるデータと異なるサンプルである必要があり、クラス毎に60か所は必要といわれる。

52