これは、Ubuntu Online、Fedora Online、Windows オンライン エミュレーター、MAC OS オンライン エミュレーターなど、複数の無料オンライン ワークステーションのいずれかを使用して、OnWorks 無料ホスティング プロバイダーで実行できるコマンド i.atcorrgrass です。
プログラム：
NAME
i.atcorr ・6Sアルゴリズムによる大気補正を行います。
6S - 太陽スペクトルにおける衛星信号の XNUMX 番目のシミュレーション。
KEYWORDS
画像、大気補正
SYNOPSIS
i.atcorr
i.atcorr - 助けて
i.atcorr [-イラク人] =名 [範囲=最小、最大] [標高=名] [視認性=名]
パラメータ=名 出力=名 [再スケール=最小、最大] [-上書きする] [-助けます] [-詳細]
[-静かな] [-ui]
フラグ：
-i
ラスター マップを整数として出力する
-r
反射率に変換された入力ラスター マップ (デフォルトは放射輝度)
-a
1 年 2000 月 XNUMX 日以降に撮影された ETM+ 画像からの入力
-b
1 年 2000 月 XNUMX 日より前に撮影された ETM+ 画像からの入力
-上書き
出力ファイルが既存のファイルを上書きできるようにする
- 助けて
使用状況の概要を印刷する
-詳細
冗長モジュール出力
- 静かな
静かなモジュール出力
--ui
GUIダイアログを強制的に起動する
パラメーター：
=名 [必要]
入力ラスターマップの名前
範囲=最小、最大
入力範囲
デフォルト： 0,255
標高=名
入力標高ラスター マップの名前 (m)
視認性=名
入力可視性ラスター マップの名前 (km)
パラメータ=名 [必要]
6S パラメータを含む入力テキスト ファイルの名前
出力=名 [必要]
出力ラスターマップの名前
再スケール=最小、最大
出力ラスター マップの縮尺変更
デフォルト： 0,255
DESCRIPTION
i.atcorr 6S アルゴリズムを使用して、入力ラスター マップで大気補正を実行します
(秒 of 衛星 シグナル in 太陽 スペクトラム）。 詳細なアルゴリズム
説明は、Land Surface Reflectance Science Computing Facility で入手できます。
ウェブサイトをご覧ください。
重要 注意事項： 電流プローブ 地域 設定 無視！ リージョンは、
大気補正が実行される前の入力ラスター マップ。 以前の設定は
その後復元。 このフラグが教えてくれます i.atcorr 計算の高速化を試みます。 でも、
このオプションにより、メモリ要件が増加します。
イフフラグ -r が使用される場合、入力ラスター データは次のように扱われます。 反射率. それ以外の場合、入力
ラスター データは次のように扱われます。 輝き での反射率に変換されます。
i.atcorr ランタイム。 出力データは常に反射率です。
衛星陸橋の時刻は、グリニッジ標準時 (GMT) で指定する必要があることに注意してください。
6S パラメータの例:
8 - 幾何学的条件=Landsat ETM+
2 19 13.00 -47.410 -20.234 - 月 日 hh.ddd 経度 緯度 (「hh.ddd」は XNUMX 進数の時間 GMT)
1 - 大気モード = トロピカル
1 - エアロゾル モデル = 大陸
15 - 可視性 [km] (エアロゾル モデル濃度)
-0.600 - 平均目標海抜 [km] (ここでは海抜 600m)
-1000 - センサーの高さ (ここでは、衛星に搭載されたセンサー)
64 - ETM+ ランドサット 4 の第 7 バンド
位置が緯度経度 (WGS84) で利用できない場合、 m.プロジェクト 変換
モジュールを使用して、別の投影から再投影できます。
6S CODE パラメータ 選択肢
A. 幾何学的 条件
CPコード 詳細説明 Details
1 隕石 月、日、XNUMX 進数の時間 (世界時 - hh.ddd) を入力します。 n. の
列、n。 行の。 (フルスケール5000×2500)
2 goes 東 月、日、XNUMX 進数の時間 (世界時 - hh.ddd) を入力します。 n. の
列、n。 行の。 (フルスケール 17000*12000)c
3 goes 西 月、日、XNUMX 進数の時間 (世界時 - hh.ddd) を入力します。 n. の
列、n。 行の。 (フルスケール17000×12000)
4 avhrr (PM noaa) 月、日、XNUMX 進数の時間を入力 (世界時 - hh.ddd) n. の
column(1-2048),xlonan,hna long.(xlonan) と高架
赤道のアセンダント ノードでの時間 (hna)
5 avhrr (AM noaa) 月、日、XNUMX 進数の時 (世界時 - hh.ddd) を入力します。 の
column(1-2048),xlonan,hna long.(xlonan) と高架
赤道のアセンダント ノードでの時間 (hna)
6 hrv (スポット) 月、日、hh.ddd、経度、緯度を入力します。 *
7 tm (landsat) 月、日、hh.ddd、long.、lat を入力します。 *
8 etm+ (landsat7) 月、日、hh.ddd、経度、緯度を入力します。 *
9 LISS (IRS 1C) 月、日、hh.ddd、経度、緯度を入力します。 *
10 アスター 月、日、hh.ddd、long.、lat を入力します。 *
11 アヴニール 月、日、hh.ddd、long.、lat を入力します。 *
12 イコノス 月、日、hh.ddd、long.、lat を入力します。 *
13 RapidEye 月、日、hh.ddd、long.、lat を入力します。 *
14 VGT1 (スポット4) 月、日、hh.ddd、long.、lat を入力します。 *
15 VGT2 (スポット5) 月、日、hh.ddd、long.、lat を入力します。 *
16 世界観 2 月、日、hh.ddd、long.、lat を入力します。 *
17 クイックバード 月、日、hh.ddd、long.、lat を入力します。 *
18 ランドサット 8 月、日、hh.ddd、long.、lat を入力します。 *
* 注意: HRV、TM、ETM+、LISS、および ASTER 実験の場合、経度と緯度は
シーン中心の座標。 緯度は、北半球では > 0、北半球では < 0 でなければなりません
南方の。 経度は、東半球では > 0、西半球では < 0 でなければなりません。
B. 大気の
CPコード 意味
0 ガス吸収なし
1 トロピカル
2 中緯度の夏
3 中緯度の冬
4 亜寒帯の夏
5 亜寒帯の冬
6us 標準 62
7 独自の大気モデルを次のセットとして定義します
測定ごとに 5 つのパラメータ: 高度 [km] 気圧
[mb] 温度 [k] h2o 密度 [g/m3] o3 密度 [g/m3]
例: それぞれに XNUMX つのラジオゾンデ測定値があります。
高度 0 ～ 25km を 1km ごとに XNUMX 回測定
25km刻みで高度50～5km、シングルXNUMX台
高度 70km および 100km の測定。 これで34になります
測定。 その場合、入力する値は 34*5 です。
8 独自の大気モデルを定義して、
水蒸気とオゾンの含有量: uw [g/cm2] uo3 [cm-atm]
プロフィールはus62から取っています。
C. エアロゾル
CPコード 意味 Details
0 エアロゾルなし
1大陸モデル
2 海事モデル
3 都市モデル
背景の砂漠エアロゾルの 4 シェル モデル
5 バイオマス燃焼
6成層圏モデル
7 独自のモデルを定義する 各コンポーネントの体積割合を入力します。 c（1）=
ダスト状の体積％ c(2) = 水溶性の体積%
c(3) = 海洋の体積 % c(4) = すすの体積 % すべて
0 から 1 の間の値。
8 独自のモデルを定義 サイズ分布関数: Multimodal Log Normal (最大 4
モード）。
9 独自のモデルを定義する サイズ分布関数: 修正ガンマ。
10 独自のモデルを定義 サイズ分布関数: Junge Power-Law。
11 独自のモデルを定義する 太陽光度計の測定値、最大 50 の値、次のように入力: r
および d V / d (logr) ここで、r は半径 [ミクロン]、V は
体積、d V / d (logr) [cm3/cm2/ミクロン]。 続いて: nr
nr と ni はそれぞれ、各波長の ni
屈折率の実部と虚部。
D. エアロゾル 濃度 (可視性)
気象パラメーターの可視性 v の推定値がある場合は、
v [km] の値 (エアロゾルの光学的深度 (AOD) は、標準的なエアロゾルから計算されます。
プロフィール）。
エアロゾルの光学的深さの推定値がある場合は、可視性に 0 を入力し、
次の行は、550 nm でのエアロゾルの光学的深さを入力します (iaer は入力の「i」を意味し、
エアロゾルの場合は「aer」)、例:
0 - 可視性
0.112 - エアロゾルの光学的深さ 550 nm
注: iaer が 0 の場合、可視性のために -1 を入力します。
E. ターゲット 高度 (xps)、 センサー プラットフォーム (xpp)
ターゲット高度 (xps、負の [km]): xps >= 0 は、ターゲットが海面にあることを意味します。
それ以外の場合、xps はターゲットの高度 (平均高度など) を [km] で表します。
負の値として
センサー プラットフォーム (xpp、マイナス [km] または -1000):
xpp = -1000 は、センサーが衛星に搭載されていることを意味します。
xpp = 0 は、センサーが地上にあることを意味します。
-100 < xpp < 0 は、[km] で表されるセンサーの高度を定義します。 この標高は
与えられた 相対 〜へ ターゲット 負の値としての高度。
航空機シミュレーションのみ (xpp は 0 でも -1000 でもありません): puw,po3
(航空機と表面の間の水蒸気含有量、オゾン含有量)
taerp (航空機と表面の間の 550nm でのエアロゾルの光学的厚さ)
これらのデータが利用できない場合は、すべて負の値を入力してください。 puw、po3はその後
地上での値に従って、us62 標準プロファイルから補間されます
レベル。 taerp は、エアロゾルの 2km 指数プロファイルに従って計算されます。
F. センサー バンド
2 つの可能性があります: 独自のスペクトル条件を定義する (コード -1、-0、XNUMX、
または 1) または、事前定義された衛星の XNUMX つの帯域を示すコードを選択します。
独自のスペクトル条件を定義します:
CPコード 意味
-2 wlinf、wlsup と入力します。 フィルタ関数は 1 になります
帯域全体（iwave = 0として）ですが、段階的な出力
印刷されます。
-1 wl を入力します (単色条件、ガス吸収を含む)。
0 wlinf、wlsup と入力します。 フィルタ関数は次のようになります
1 バンド全体で。
1 wlinf、wlsup、およびユーザーのフィルター関数 s(lambda) を次のように入力します。
0.0025マイクロメートルのステップ。
事前定義された衛星バンド:
CPコード 意味
2 隕石 vis バンド (0.350-1.110)
3 goes 東 バンド vis (0.490-0.900)
4 西側バンド vis (0.490-0.900)
5 avhrr (ノアア6) バンド 1 (0.550-0.750)
6 avhrr (noaa6) バンド 2 (0.690-1.120)
7 avhrr (ノアア7) バンド 1 (0.500-0.800)
8 avhrr (noaa7) バンド 2 (0.640-1.170)
9 avhrr (ノアア8) バンド 1 (0.540-1.010)
10 avhrr (noaa8) バンド 2 (0.680-1.120)
11 avhrr (ノアア9) バンド 1 (0.530-0.810)
12 avhrr (noaa9) バンド 1 (0.680-1.170)
13 avhrr (ノアア10) バンド 1 (0.530-0.780)
14 avhrr (noaa10) バンド 2 (0.600-1.190)
15 avhrr (ノアア11) バンド 1 (0.540-0.820)
16 avhrr (noaa11) バンド 2 (0.600-1.120)
17 hrv1 (スポット1) バンド 1 (0.470-0.650)
18 hrv1 (スポット 1) バンド 2 (0.600-0.720)
19 hrv1 (スポット 1) バンド 3 (0.730-0.930)
20 hrv1 (spot1) バンドパン (0.470-0.790)
21 hrv2 (スポット1) バンド 1 (0.470-0.650)
22 hrv2 (スポット 1) バンド 2 (0.590-0.730)
23 hrv2 (スポット 1) バンド 3 (0.740-0.940)
24 hrv2 (spot1) バンドパン (0.470-0.790)
25 tm (ランドサット5) バンド 1 (0.430-0.560)
26 tm (landsat5) バンド 2 (0.500-0.650)
27 tm (landsat5) バンド 3 (0.580-0.740)
28 tm (landsat5) バンド 4 (0.730-0.950)
29 tm (landsat5) バンド 5 (1.5025-1.890)
30 tm (landsat5) バンド 7 (1.950-2.410)
31 MSS (ランドサット5) バンド 1 (0.475-0.640)
32 ミリ秒 (landsat5) バンド 2 (0.580-0.750)
33 ミリ秒 (landsat5) バンド 3 (0.655-0.855)
34 ミリ秒 (landsat5) バンド 4 (0.785-1.100)
35 MAS (ER2) バンド 1 (0.5025-0.5875)
36 MAS (ER2) バンド 2 (0.6075-0.7000)
37 MAS (ER2) バンド 3 (0.8300-0.9125)
38 MAS (ER2) バンド 4 (0.9000-0.9975)
39 MAS (ER2) バンド 5 (1.8200-1.9575)
40 MAS (ER2) バンド 6 (2.0950-2.1925)
41 MAS (ER2) バンド 7 (3.5800-3.8700)
42 モディス バンド 1 (0.6100-0.6850)
43 MODIS バンド 2 (0.8200-0.9025)
44 MODIS バンド 3 (0.4500-0.4825)
45 MODIS バンド 4 (0.5400-0.5700)
46 MODIS バンド 5 (1.2150-1.2700)
47 MODIS バンド 6 (1.6000-1.6650)
48 MODIS バンド 7 (2.0575-2.1825)
49 avhrr (ノアア12) バンド 1 (0.500-1.000)
50 avhrr (noaa12) バンド 2 (0.650-1.120)
51 avhrr (ノアア14) バンド 1 (0.500-1.110)
52 avhrr (noaa14) バンド 2 (0.680-1.100)
53 干拓地 バンド 1 (0.4125-0.4775)
54 POLDER バンド 2 (無極性) (0.4100-0.5225)
55 POLDER バンド 3 (無極性) (0.5325-0.5950)
56 POLDER バンド 4 P1 (0.6300-0.7025)
57 POLDER バンド 5 (無極性) (0.7450-0.7800)
58 POLDER バンド 6 (無極性) (0.7000-0.8300)
59 POLDER バンド 7 P1 (0.8100-0.9200)
60 POLDER バンド 8 (無極性) (0.8650-0.9400)
61 etm+ (ランドサット7) バンド 1 (0.435-0.520)
62 etm+ (landsat7) バンド 2 (0.506-0.621)
63 etm+ (landsat7) バンド 3 (0.622-0.702)
64 etm+ (landsat7) バンド 4 (0.751-0.911)
65 etm+ (landsat7) バンド 5 (1.512-1.792)
66 etm+ (landsat7) バンド 7 (2.020-2.380)
67 etm+ (landsat7) バンド 8 (0.504-0.909)
68 LISS (IRC 1C） バンド 2 (0.502-0.620)
69 liss (IRC 1C) バンド 3 (0.612-0.700)
70 liss (IRC 1C) バンド 4 (0.752-0.880)
71 liss (IRC 1C) バンド 5 (1.452-1.760)
72 アスター バンド 1 (0.480-0.645)
73 アスターバンド 2 (0.588-0.733)
74 アスターバンド 3N (0.723-0.913)
75 アスターバンド 4 (1.530-1.750)
76 アスターバンド 5 (2.103-2.285)
77 アスターバンド 6 (2.105-2.298)
78 アスターバンド 7 (2.200-2.393)
79 アスターバンド 8 (2.248-2.475)
80 アスターバンド 9 (2.295-2.538)
81 アヴニール バンド 1 (0.390-0.550)
82 avnir バンド 2 (0.485-0.695)
83 avnir バンド 3 (0.545-0.745)
84 avnir バンド 4 (0.700-0.925)
85 イコノス グリーンバンド (0.350-1.035)
86 ikonos 赤帯 (0.350-1.035)
87 ikonos NIR バンド (0.350-1.035)
88 RapidEye ブルーバンド (0.438-0.513)
89 RapidEye グリーンバンド (0.463-0.594)
90 RapidEye レッドバンド (0.624-0.690)
91 RapidEye RedEdge バンド (0.500-0.737)
92 RapidEye NIR バンド (0.520-0.862)
93 VGT1 (スポット4) バンド 0 (0.400-0.500)
94 VGT1 (SPOT4) バンド 2 (0.580-0.782)
95 VGT1 (SPOT4) バンド 3 (0.700-1.030)
96 VGT1 (SPOT4) MIR バンド (1.450-1.800)
97 VGT2 (スポット5) バンド 0 (0.400-0.550)
98 VGT2 (SPOT5) バンド 2 (0.580-0.780)
99 VGT2 (SPOT5) バンド 3 (0.700-1.000)
100 VGT2 (SPOT5) MIR バンド (1.450-1.800)
101 WorldView 2 パンクロマティック バンド (0.447-0.808)
102 WorldView 2 Coastal Blue バンド (0.396-0.458)
103 WorldView 2 ブルーバンド (0.442-0.515)
104 WorldView 2 グリーンバンド (0.506-0.586)
105 WorldView 2 黄色のバンド (0.584-0.632)
106 WorldView 2 レッドバンド (0.624-0.694)
107 WorldView 2 レッドエッジバンド (0.699-0.749)
108 WorldView 2 NIR1 バンド (0.765-0.901)
109 WorldView 2 NIR2 バンド (0.856-0.1043)
110 クイックバード パンクロバンド (0.405-1.053)
111 QuickBird ブルーバンド (0.430-0.545)
112 QuickBird グリーン バンド (0.466-0.620)
113 QuickBird レッドバンド (0.590-0.710)
114 QuickBird NIR1 バンド (0.715-0.918)
115 ランドサット 8 沿岸エアロゾルバンド (0.427nm - 0.459nm)
116 Landsat 8 ブルーバンド (436nm - 527nm)
117 Landsat 8 グリーンバンド (512nm-610nm)
118 ランドサット 8 レッドバンド (625nm-691nm)
119 ランドサット 8 パンクロマティック バンド (488nm-692nm)
120 Landsat 8 NIR バンド (829nm-900nm)
121 ランドサット 8 シーラス バンド (1340nm-1409nm)
122 ランドサット 8 SWIR1 バンド (1515nm - 1697nm)
123 ランドサット 8 SWIR2 バンド (2037nm - 2355nm)
例
大気の 補正 of a LANDSAT-7 チャンネル
この例は、ノースカロライナ州のサンプル データセット (GMT -5 時間) に基づいています。 最初に設定します
衛星地図への計算領域 (例: チャネル 4):
g.region raster=lsat7_2002_40 -p
太陽の位置の利用可能なメタデータを確認することが重要です。
大気補正用に定義されています。 オプションは、衛星高架時間を確認することです
メタデータで報告されている太陽の位置。 ノースカロライナ州のサンプル データセットの場合、
チャネルごとに保存されており、次のように取得できます。
r.info lsat7_2002_40
この場合、SUN_AZIMUTH = 120.8810347、SUN_ELEVATION = 64.7730999 です。
太陽の位置のメタデータが利用できない場合は、陸橋から計算することもできます
時間は以下の通り（r.サンマスク SOLPOS を使用):
r.sunmask -s elev=標高 out=ダミー year=2002 month=5 day=24hour=10 min=42 sec=7 timezone=-5
# .. レポート: 太陽の方位角: 121.342461、水平線上の太陽の角度 (屈折補正): 65.396652
高架通過時間が不明な場合は、サテライト 高架通過予測を使用します。
を使用して、DN (デジタル数 = ピクセル値) を大気のトップ (TOA) での放射輝度に変換します。
式
Lλ = ((LMAXλ - LMINλ)/(QCALMAX-QCALMIN)) * (QCAL-QCALMIN) + LMINλ
どこ：
· Lλ = センサーのアパーチャでの分光放射輝度 (ワット/(メートル平方 * スター)
* µm)、サテライト センサーで見た見かけの放射輝度。
· QCAL = DN の量子化された校正済みピクセル値。
· LMINλ = QCALMIN にスケーリングされた分光放射輝度 (ワット/(メートル)
二乗 * スター * µm);
· LMAXλ = QCALMAX にスケーリングされた分光放射輝度 (ワット/(メートル)
二乗 * スター * µm);
· QCALMIN = 最小の量子化されたキャリブレーション済みピクセル値 (に対応)
LMINλ) DN;
· QCALMAX = 量子化されたキャリブレーション済みピクセルの最大値 (に対応)
LMAXλ) DN = 255 で。
LMINλ と LMAXλ は、DN の最小値と最大値に関連する放射輝度です。
各画像のメタデータ ファイルまたは表 1 で報告されます。 高ゲインまたは低ゲイン
ゲインは、各 Landsat 画像のメタデータ ファイルでも報告されます。 DN の最小値
(QCALMIN) は Landsat ETM+ 画像の場合は 1 です (Landsat ハンドブック、第 11 章を参照)。
最大 DN 値 (QCALMAX) は 255 です。QCAL は、各ピクセルの DN 値です。
ランドサット画像。
放射輝度マップを取得するために、係数を抽出して適用します。
ちゃん=4
r.info lsat7_2002_${CHAN}0 -h | tr '\n' ' ' | sed 's+ ++g' | tr ':' '\n' | grep "LMIN_BAND${CHAN}\|LMAX_BAND${CHAN}"
LMAX_BAND4=241.100,p016r035_7x20020524.met
LMIN_BAND4=-5.100,p016r035_7x20020524.met
QCALMAX_BAND4=255.0,p016r035_7x20020524.met
QCALMIN_BAND4=1.0,p016r035_7x20020524.met
放射輝度への変換 (この計算はバンド 4 に対して行われ、他のバンドでは、
斜体の数字は関連する値に置き換える必要があります):
r.mapcalc "lsat7_2002_40_rad = ((241.1 - (-5.1)) / (255.0 - 1.0)) * (lsat7_2002_40 - 1.0) + (-5.1)"
# 平均標高を見つける (制御ファイルで初期化値として使用される海抜のターゲット)
r.univar標高
メタデータに基づいて、チャネル 4 (NIR) の制御ファイル 'icnd.txt' を作成します。 高架用
時間、XNUMX 進数の時間を定義する必要があります。
10:42:07 NC 現地時間 = 10.70 42 進時 (100 進分: 60 * XNUMX / XNUMX)、つまり
15.70GMT:
8 - 幾何学的条件=Landsat ETM+
5 24 15.70 -78.691 35.749 - 月 日 hh.ddd 経度 緯度 (「hh.ddd」は GMT XNUMX 進時)
2 - 大気モード=中緯度の夏
1 - エアロゾル モデル = 大陸
50 - 可視性 [km] (エアロゾル モデル濃度)
-0.110 - 海抜の平均目標標高 [km]
-1000 - 衛星搭載センサー
64 - ETM+ ランドサット 4 の第 7 バンド
最後に、大気補正を実行します (反射率入力マップの場合は -r、日付 >XNUMX 月の場合は -a)。
2000）：
i.atcorr -r -a lsat7_2002_40_rad elev=標高パラメータ=icnd_lsat4.txt 出力=lsat7_2002_40_atcorr
「icnd_lsat4.txt」ファイルからの高度値が最初に読み取られることに注意してください
初期変換を計算します。 平均となる値を与える必要があります
標高モデルの値。 大気補正の場合はラスター標高
値はマップから使用されます。
このプロセスは計算集約的であることに注意してください。
また、次の点にも注意してください。 i.atcorr 太陽ではなく地平線上の太陽仰角を報告する
天頂角。
残り マニュアル 問題
1. ビジビリティ値またはマップの影響と重要性についても説明する必要があります。
550nm での可視性またはエアロゾルの光学的深さの推定値を取得する方法。
onworks.net サービスを使用してオンラインで i.atcorrgrass を使用する
