Landsat8 DN to TOA Reflectance (nilai minus)

8

Saya belajar cara melakukan koreksi atmosfer secara manual, karenanya tidak menggunakan metode otomatis yang tersedia di ENVI

metadata (band 1) untuk adegan Landsat 8 saya adalah sebagai berikut;

RADIANCE_MULT_BAND_1 = 1.2344E-02

RADIANCE_ADD_BAND_1 = -61.72016

REFLECTANCE_MULT_BAND_1 = 2.0000E-05

REFLECTANCE_ADD_BAND_1 = -0.100000

SUN_ELEVATION = 48.12237303

Jadi ketika saya menggunakan band matematika untuk menghitung TOA Radiance saya ketik ini:

(0,012344 * B1) -61.72016

Dan itu menghasilkan statistik spektral yang terlihat benar Statistik Spectral

Ketika saya menggunakan salah satu rumus matematika band berikut;

TOA Reflectance (tidak termasuk penyesuaian ketinggian matahari)

(B1 * 0,00002) -0.1

TOA Reflectance (termasuk penyesuaian ketinggian matahari)

((B1 * 0,00002) -0.1) /0.74457226676389733207607359928648

Atau hitung langsung dari DN mentah ke TOA Reflectance (termasuk penyesuaian ketinggian matahari)

(((((0,012344 * B1) -61.72016) * 0.00002) -0.1) /0.74457226676389733207607359928648

Saya berakhir dengan sesuatu seperti ini:

Reflektansi TOA

yang terlihat salah.

Saya sudah googled dan istilah pencarian saya adalah sampah atau tidak ada yang menggunakan matematika band manual untuk melakukan ini.

[EDIT]

Jawaban Mikkel di bawah ini dengan perhitungan matematika band revisi (menghilangkan langkah-langkah menengah) menghasilkan statistik berikut:

masukkan deskripsi gambar di sini

yang terlihat agak lebih realistis.

Ingin tahu apakah saya harus menambahkan 13,4305 ke setiap band, atau mengabaikan minimum negatif?

[AKHIR EDIT]

[EDIT # 2]

Mikkel tepat lagi, dalam menyarankan bahwa nilai -13,4305 diisolasi ke area "tidak ada data" di sekitar gambar;

masukkan deskripsi gambar di sini

seperti yang Anda lihat, terisolasi hanya ke tepi utara tempat kejadian.

[AKHIR EDIT # 2]

Kyuss
sumber

Jawaban:

7

Satu kelemahan dalam pendekatan Anda. Anda tidak perlu melalui DN untuk bersinar.
Anda dapat langsung menuju DN ke reflektansi.
Tetap berpegang pada ((B1*0.00002)-0.1)/0.74457226676389733207607359928648.

Mikkel Lydholm Rasmussen
sumber
ini bekerja dengan sangat baik! Saya kira tidak pernah dijelaskan bahwa TOA radiance / TOA reflectance excl. ketinggian matahari tidak diperlukan langkah-langkah dalam urutan koreksi.
Kyy
per pos yang diedit di atas, apakah perlu menambahkan 13,4305 ke setiap band untuk membawanya hingga min = 0, atau itu tidak perlu ??
Kyuss
@kyuss, coba periksa di mana -0.134305 berada. Saya menduga bahwa itu adalah area tanpa data di sekitar tepi gambar.
Mikkel Lydholm Rasmussen
brilian, saya bahkan tidak memeriksa apa pun kecuali perbatasan selatan di mana tidak ada nilai negatif. terima kasih atas semua bantuan anda.
Kyuss
@kyuss sekarang langkah selanjutnya dalam proses Anda adalah "hanya" untuk melakukan koreksi atmosfer.
Mikkel Lydholm Rasmussen
2

Berikut ini adalah bacaan yang bagus tentang fisika koreksi atmosfer, dan mungkin membantu dengan masalah ini ... Saya seorang pria yang menekan tombol, saya merasakan sakit Anda, setidaknya Anda bisa mengklik tombol QUAC atau FLAASH yang indah begitu Anda berada di luar akademisi.

http://modis.gsfc.nasa.gov/data/atbd/atbd_mod08.pdf

geodranik
sumber
terima kasih Steven. Saya telah - seperti yang telah Anda simpulkan - belajar banyak selama studi saya tetapi dokumen yang modis mungkin membantu mengisi kesenjangan pemahaman saya kembali: transfer radiatif.
Kyuss