Kalibrasi Radiometrik: Mengubah Digital Number (DN) ke Radiance dan/atau Reflectance
Posted on: 10 Dec, 2013, by : lmjaelani

Koreksi Radiometrik (Radiometric Correction) dikelompokkan menjadi dua:

  1. Kalibrasi Radiometrik (radiometric calibration) (bagian ini akan dijelaskan di sini)
  2. Koreksi Atmosferik (atmospheric correction) (bagian ini akan dijelaskan kemudian, sementara bisa membaca paper saya di http://lmjaelani.com/2013/11/24/evaluation-of-four-meris-atmospheric-correction-algorithms-in-lake-kasumigaura-japan/ )

Kalibrasi Radiometrik  merupakan langkah pertama yang harus dilakukan saat kita mengolah data citra satelit. Tujuan utama dari Kalibrasi radiometrik ini adalah untuk mengubah data pada citra yang (pada umumnya) disimpan dalam bentuk Digital Number (DN) menjadi radiance dan/atau reflectance, bisa juga ke brightness temperature (untuk kanal Termal Infra Red)

Terkait dengan Kalibrasi ini, ada istilah yang perlu diperhatikan, yakni Resolusi Radiometrik (radiometric resolution) yang menunjukkan `Berapa banyak bit yang digunakan dalam satu pixel?` lihat penjelasannya di  4 Resolusi penting dalam INDERAJA

Sebelum melakukan kalibrasi radiometrik, penting untuk mengetahui resolusi radiometrik dari citra yang kita gunakan. Misal untuk Landsat 7,  resolusi radiometriknya adalah 8 bits atau setara dengan ( 2 pangkat 8 ) 256 pixel value, (atau Digital Number, atau Digital Count atau dikenal juga dengan istilah grayscale),  berarti data yang kita gunakan memiliki gradasi grayscale dari 0 sampai 255. Contoh lainnya adalah Landsat 8 yang mengggunakan 16 bits, artinya digital number (DN) terletak dalam rentang 0 sampai 2 pangkat 16 (==65536).

Sekarang, kita lanjutkan dengan kalibrasi radiometriknya. Ada dua cara yang umum digunakan (tergantung data yang tersedia):

  1. Menggunakan Gain dan Offset, data yang diperlukan adalah radiance atau reflectance  multiple rescalling factor (GAIN) dan additive rescalling factor (OFFSET)
  2. Menggunakan nilai radiance atau reflectance maksimum dan minimum.

Sebagai contoh, kalibrasi radiometrik ini menggunakan data Landsat 8 (ingat, resolusi spasialnya 16 bits)

CARA PERTAMA : GAIN-OFFSET

a). Konversi DN  ke  TOA Radiance (catatan: TOA = Top of atmosphere)

Rumus yang digunakan untuk mengubah DN ke radiance (L) adalah sebagai berikut:

Lλ = MLQcal + AL

dimana:

  • Lλ          = TOA spectral radiance (Watts/( m2 * srad * μm))
  • ML         = Band-specific multiplicative rescaling factor from the metadata (RADIANCE_MULT_BAND_x, where x is the band number)
  • AL          = Band-specific additive rescaling factor from the metadata (RADIANCE_ADD_BAND_x, where x is the band number)
  • Qcal        = Quantized and calibrated standard product pixel values (DN)

Radiance rescaling factors  di atas bisa ditemukan di file metadata (yang selalu menyertai setiap data citra satelit)

b). Konversi DN ke  TOA Reflectance (ρ)
Sama dengan point (a) di atas, konversi DN bisa juga dilakukan langsung menjadi reflectance, gunakan rumus berikut:

ρλ = MρQcal + Aρ

dimana;

  • ρλ          = TOA planetary reflectance, without correction for solar angle.  Note that ρλ’ does not contain a correction for the sun angle. — untuk mengoreksi reflectance terhadap sudut matahari  lihat CARA KEDUA bagian b
  • Mρ         = Band-specific multiplicative rescaling factor from the metadata (REFLECTANCE_MULT_BAND_x, where x is the band number)
  • Aρ          = Band-specific additive rescaling factor from the metadata (REFLECTANCE_ADD_BAND_x, where x is the band number)
  • Qcal        = Quantized and calibrated standard product pixel values (DN)

http://landsat.usgs.gov/Landsat8_Using_Product.php

CARA KEDUA : MAX-MIN

a). Konversi DN  ke  TOA Radiance

Lλ={(Lmax-Lmin)/(Qcalmax-Qcalmin)} * (Qcal-Qcalmax)

Dimana:
Lλ = Spectral Radiance in watts/(meter squared * ster * μm)
LMAXR = Max Detected Radiance Level
LMINR = Min Detected Radiance Level
QCALMAX = Max Pixel Value
QCALMIN = Min Pixel Value
QCAL = Digital Number

b). Konversi DN ke  TOA Reflectance (ρ) cara 1

ρλ`={(ρmax-ρmin)/(Qcalmax-Qcalmin)} * (Qcal-Qcalmax)

Dimana:
ρλ = Spectral reflectance (tanpa koreksi solar angle)
ρMAXR = Max Detected Reflectance Level
ρMINR = Min Detected Reflectance Level
QCALMAX = Max Pixel Value
QCALMIN = Min Pixel Value
QCAL = Digital Number

Berikut langkah tambahan untuk mengoreksi reflectance dengan solar angle

ρλ* =ρλ/cos(</em><em>θ<sub>SZ</sub></em>) atau <em>ρλ/sin(θSE)

dimana:

ρλ*          = TOA planetary reflectance
θSE         = Local sun elevation angle. The scene center sun elevation angle in degrees is provided in the metadata (SUN_ELEVATION). Untuk penghitungan yang akurat, solar angle harus pada pixel yang dihitung, namun sayang sekali di metadata hanya disediakan sudut matahari untuk pusat scene saja.
θSZ         = Local solar zenith angle;  θSZ = 90° – θSE

c). Konversi DN ke  TOA Reflectance (ρ) cara 2 (tidak usah dipakai, merepotkan)

ρ=(π*Lλ *d^2)/(Esunλ *cosθ)

dimana:
ρ = Unitless Planetary Reflectance
Lλ = Spectral Radiance at The Sensor’s Aperture
d2 = Earth-Sun Distance in Astronomical Units
ESUNλ = Mean Solar Exoatmospheric Irradiances
θ = Solar Zenith Angle in Degrees

ESUN Landsat 8
ESUN Landsat 8

CONTOH ISI METADATA (yang terkait dengan kalibrasi radiometrik)

GROUP = MIN_MAX_RADIANCE –>LMAXR dan LMINR
RADIANCE_MAXIMUM_BAND_1 = 775.47809
RADIANCE_MINIMUM_BAND_1 = -64.03923
RADIANCE_MAXIMUM_BAND_2 = 790.78357
RADIANCE_MINIMUM_BAND_2 = -65.30316
RADIANCE_MAXIMUM_BAND_3 = 724.09100
RADIANCE_MINIMUM_BAND_3 = -59.79567
RADIANCE_MAXIMUM_BAND_4 = 613.26074
RADIANCE_MINIMUM_BAND_4 = -50.64327
RADIANCE_MAXIMUM_BAND_5 = 372.13345
RADIANCE_MINIMUM_BAND_5 = -30.73090
RADIANCE_MAXIMUM_BAND_6 = 93.75905
RADIANCE_MINIMUM_BAND_6 = -7.74265
RADIANCE_MAXIMUM_BAND_7 = 30.49965
RADIANCE_MINIMUM_BAND_7 = -2.51867
RADIANCE_MAXIMUM_BAND_8 = 690.79944
RADIANCE_MINIMUM_BAND_8 = -57.04644
RADIANCE_MAXIMUM_BAND_9 = 152.92506
RADIANCE_MINIMUM_BAND_9 = -12.62860
RADIANCE_MAXIMUM_BAND_10 = 22.00180
RADIANCE_MINIMUM_BAND_10 = 0.10033
RADIANCE_MAXIMUM_BAND_11 = 22.00180
RADIANCE_MINIMUM_BAND_11 = 0.10033
END_GROUP = MIN_MAX_RADIANCE

GROUP = MIN_MAX_REFLECTANCE –>ρMAXR, ρMINR
REFLECTANCE_MAXIMUM_BAND_1 = 1.210700
REFLECTANCE_MINIMUM_BAND_1 = -0.099980
REFLECTANCE_MAXIMUM_BAND_2 = 1.210700
REFLECTANCE_MINIMUM_BAND_2 = -0.099980
REFLECTANCE_MAXIMUM_BAND_3 = 1.210700
REFLECTANCE_MINIMUM_BAND_3 = -0.099980
REFLECTANCE_MAXIMUM_BAND_4 = 1.210700
REFLECTANCE_MINIMUM_BAND_4 = -0.099980
REFLECTANCE_MAXIMUM_BAND_5 = 1.210700
REFLECTANCE_MINIMUM_BAND_5 = -0.099980
REFLECTANCE_MAXIMUM_BAND_6 = 1.210700
REFLECTANCE_MINIMUM_BAND_6 = -0.099980
REFLECTANCE_MAXIMUM_BAND_7 = 1.210700
REFLECTANCE_MINIMUM_BAND_7 = -0.099980
REFLECTANCE_MAXIMUM_BAND_8 = 1.210700
REFLECTANCE_MINIMUM_BAND_8 = -0.099980
REFLECTANCE_MAXIMUM_BAND_9 = 1.210700
REFLECTANCE_MINIMUM_BAND_9 = -0.099980
END_GROUP = MIN_MAX_REFLECTANCE

GROUP = MIN_MAX_PIXEL_VALUE –>QCALMAX, QCALMIN
QUANTIZE_CAL_MAX_BAND_1 = 65535
QUANTIZE_CAL_MIN_BAND_1 = 1
QUANTIZE_CAL_MAX_BAND_2 = 65535
QUANTIZE_CAL_MIN_BAND_2 = 1
QUANTIZE_CAL_MAX_BAND_3 = 65535
QUANTIZE_CAL_MIN_BAND_3 = 1
QUANTIZE_CAL_MAX_BAND_4 = 65535
QUANTIZE_CAL_MIN_BAND_4 = 1
QUANTIZE_CAL_MAX_BAND_5 = 65535
QUANTIZE_CAL_MIN_BAND_5 = 1
QUANTIZE_CAL_MAX_BAND_6 = 65535
QUANTIZE_CAL_MIN_BAND_6 = 1
QUANTIZE_CAL_MAX_BAND_7 = 65535
QUANTIZE_CAL_MIN_BAND_7 = 1
QUANTIZE_CAL_MAX_BAND_8 = 65535
QUANTIZE_CAL_MIN_BAND_8 = 1
QUANTIZE_CAL_MAX_BAND_9 = 65535
QUANTIZE_CAL_MIN_BAND_9 = 1
QUANTIZE_CAL_MAX_BAND_10 = 65535
QUANTIZE_CAL_MIN_BAND_10 = 1
QUANTIZE_CAL_MAX_BAND_11 = 65535
QUANTIZE_CAL_MIN_BAND_11 = 1
END_GROUP = MIN_MAX_PIXEL_VALUE

GROUP = RADIOMETRIC_RESCALING –>ML , AL, Mρ,  Aρ
RADIANCE_MULT_BAND_1 = 1.2810E-02
RADIANCE_MULT_BAND_2 = 1.3063E-02
RADIANCE_MULT_BAND_3 = 1.1962E-02
RADIANCE_MULT_BAND_4 = 1.0131E-02
RADIANCE_MULT_BAND_5 = 6.1474E-03
RADIANCE_MULT_BAND_6 = 1.5488E-03
RADIANCE_MULT_BAND_7 = 5.0383E-04
RADIANCE_MULT_BAND_8 = 1.1412E-02
RADIANCE_MULT_BAND_9 = 2.5262E-03
RADIANCE_MULT_BAND_10 = 3.3420E-04
RADIANCE_MULT_BAND_11 = 3.3420E-04
RADIANCE_ADD_BAND_1 = -64.05204
RADIANCE_ADD_BAND_2 = -65.31623
RADIANCE_ADD_BAND_3 = -59.80763
RADIANCE_ADD_BAND_4 = -50.65340
RADIANCE_ADD_BAND_5 = -30.73705
RADIANCE_ADD_BAND_6 = -7.74420
RADIANCE_ADD_BAND_7 = -2.51917
RADIANCE_ADD_BAND_8 = -57.05785
RADIANCE_ADD_BAND_9 = -12.63113
RADIANCE_ADD_BAND_10 = 0.10000
RADIANCE_ADD_BAND_11 = 0.10000

REFLECTANCE_MULT_BAND_1 = 2.0000E-05
REFLECTANCE_MULT_BAND_2 = 2.0000E-05
REFLECTANCE_MULT_BAND_3 = 2.0000E-05
REFLECTANCE_MULT_BAND_4 = 2.0000E-05
REFLECTANCE_MULT_BAND_5 = 2.0000E-05
REFLECTANCE_MULT_BAND_6 = 2.0000E-05
REFLECTANCE_MULT_BAND_7 = 2.0000E-05
REFLECTANCE_MULT_BAND_8 = 2.0000E-05
REFLECTANCE_MULT_BAND_9 = 2.0000E-05
REFLECTANCE_ADD_BAND_1 = -0.100000
REFLECTANCE_ADD_BAND_2 = -0.100000
REFLECTANCE_ADD_BAND_3 = -0.100000
REFLECTANCE_ADD_BAND_4 = -0.100000
REFLECTANCE_ADD_BAND_5 = -0.100000
REFLECTANCE_ADD_BAND_6 = -0.100000
REFLECTANCE_ADD_BAND_7 = -0.100000
REFLECTANCE_ADD_BAND_8 = -0.100000
REFLECTANCE_ADD_BAND_9 = -0.100000
END_GROUP = RADIOMETRIC_RESCALING

TAMBAHAN:

Kita bisa mengubah DN ke brightness temperature menggunakan rumus:

T=K2/{ln(k1/Lλ +1)}

dimana:

T           = At-satellite brightness temperature (K)
Lλ          = TOA spectral radiance (Watts/( m2 * srad * μm))
K1          = Band-specific thermal conversion constant from the metadata (K1_CONSTANT_BAND_x, where x is the band number, 10 or 11)
K2          = Band-specific thermal conversion constant from the metadata (K2_CONSTANT_BAND_x, where x is the band number, 10 or 11)

Data yang diperlukan juga ada di metadata, contohnya sebagai berikut:

GROUP = TIRS_THERMAL_CONSTANTS
K1_CONSTANT_BAND_10 = 774.89
K1_CONSTANT_BAND_11 = 480.89
K2_CONSTANT_BAND_10 = 1321.08
K2_CONSTANT_BAND_11 = 1201.14
END_GROUP = TIRS_THERMAL_CONSTANTS

Brightness Temperature (T) bukanlah suhu di permukaan bumi. Untuk mendapatkan suhu tersebut, T harus diubah ke Surface Temperature dengan memanfaatkan data Emmisivity (yang diturunkan dari data NDVI). –>

https://www.youtube.com/watch?v=7W4IwlvPLbQ

 

Semoga ada manfaatnya

BEBERAPA TOOL:

– earth-solar distance d

-date to DOY conversion : http://quake.geo.berkeley.edu/convert/dd.html

46 thoughts on “Kalibrasi Radiometrik: Mengubah Digital Number (DN) ke Radiance dan/atau Reflectance

  1. mas mau tanya jadi klo mau koreksi radiometric untuk landsat 8 itu harus lengkap ya, dn to rad, rad to ref, ref to sudut elevasi, trus band yang di pake juga band hasil proses bertahap itu kan ?

  2. Betul pak, tapi bisa juga langsung dari DN ke Ref (sesuai rumus di atas, mohon periksa satu satu). Semua band dipakai bertahap (band 1-9)

    Maaf, telat reply, masuk spam

  3. Ref nya ada dua, ada yang rumusnya sudah terkoreksi oleh sub elevation ada yang belum. Jika belum tinggal di koreksi dengan :

    refterkoreksi = refbelumterkoreksi/ cos (sudut elevasi matahari)

  4. assalamualaikum pak
    nama saya jeddah
    saya ingin bertanya yang mencari harga TI pada software ENVI 4.7, untuk fungsi rumus “ln” itu penulisannya di band math bagaimana ya pak ?
    mohon konfirmasinya pak
    terima kasih
    Wasalam

  5. Terima kasih banyak atas jawabannya pak
    saya ingin bertanya lagi untuk menghitung rumus dengan akar pangkat 3, cara penulisannya di ENVI 4.7 bagaimana ya pak ?
    dan kalau misalnya di suatu rumus ada konstanta 255 untuk Landsat 7, untuk Landsat 8 apakah ttp 255 atau tidak pak ?
    mohon konfirmasinya pak
    terima kasih
    Wasalam

  6. Terima kasih untuk tutorialnya mas, sangat bermanfaat. Saya mau bertanya mas. Apa bedanya reflectance dan radiance? Maksud saya, misalnya saya ingin melakukan klasifikasi citra landsat, kemudian mana yang harus saya pilih,menggunakan nilai reflectance ataukah radiance? Apa bedanya mas?

    Terima kasih, salam

    1. Kalau sekedar untuk klasifikasi citra, tidak perlu dikalibrasi pak. cukup dalam format DN.

      Radiance: total energi yang sampai di sensor per luas area tertentu (dari objek), satuannya W/(m2 * sr * µm).
      Reflectance: rasio dari energi yang dipantulkan oleh objek (dan diterima oleh sensor) dengan energi yang sampai di objek itu, satuannya tidak ada (dimensionless)

  7. maaf mas, saya bingung banget ni. masih mau coba soalnya, boleh tunjukin perhitungan angka yang sudah dimasukkan ke rumus? punya saya landsat 8, terimakasih

  8. Terimakasih infonya mas.
    kalau boleh tau caranya mengubah nilai DN ke Radian itu kan pakai rumus Lλ = MLQcal + AL.
    untuk melihat nilai ML,Qcal dan AL itu di lihat dari mananya mas?
    hehe maaf baru belajar ginian mas soalny.

  9. Assalamualaikum Pak jaelani, saya mau tanya koreksi radiometrik yg tepat untuk algoritma NDVI itu apakah yg merubah nilai DN ke radiance atau ke yg reflectance? Trims

    1. Wa’alaikumsalam wr. wb

      untuk NDVI, idealnya menggunakan remote sensing reflectance (Rrs), artinya reflektan yang sudah diperoleh dari tahap kalibrasi di atas, harus dikoreksi atmosfer.

        1. Karena data yang akan diekstrak berada di permukaan bumi (di bawah atmosfer), sementara data yang direkam oleh sensor adalah data yang sudah berinteraksi dengan atmosfer

  10. Assalamualaikum pak Jaelani mohon izin bertanya, untuk mengetahui sebaran terumbu karang lebih baik menggunakan koreksi radiometrik yang reflektan atau yang radian pak? terima kasih

    1. Untuk mendeteksi terumbu karang, diperlukan algoritma koreksi kolom air (Lyzenga). Algoritma ini, memerlukan input berupa radian permukaan (BOA-Radiance).
      Untuk mendapatkan radian permukaan tersebut, kita memerlukan reflektan permukaan (BOA-reflectance), sehingga perlu kalibrasi dari BOA reflectance ke BOA Radiance. Untuk mendapatkan BOA Reflectance dari TOA Reflectance, diperlukan koreksi atmosfer. TOA Reflectance ini yang dihasilkan dari kalibrasi radiometerik (DN ke ‘TOA’ reflectance) di judul artikel ini.

  11. selamat siang mas.
    maaf mas, apakah ada tutorial untuk mengkonversikannya? citra yang yang saya pakai landsat 8. mohon bantuan nyaa mas.
    terimakasih.

  12. pak mau tanya apa nanti nilai yg dihasilkan dr DN ke REF benar minus? krn nilai REFLECTANCE_ADD_BAND_1 = -0.100000 (minus)

  13. Sy hitung scr manual di excel tp hasilnya negatif pak. Apa nilai negatifnya dijadikan positif? Bgmn yah? Mohon bantuannya.

Leave a Reply

%d bloggers like this: