Comparing Atmospheric Correction Methods For Landsat Oli Data
Comparing Atmospheric Correction Methods For Landsat Oli Data
Comparing Atmospheric Correction Methods For Landsat Oli Data
Abstract. Landsat data used for monitoring activities to land cover because it has spatial resolution
and high temporal. To monitor land cover changes in an area, atmospheric correction is needed to be
performed in order to obtain data with precise digital value picturing current condition. This study
compared atmospheric correction methods namely Quick Atmospheric Correction (QUAC), Dark Object
Subtraction (DOS) and Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH). The
correction results then were compared to Surface Reflectance (SR) imagery data obtained from the
United States Geological Survey (USGS) satelite. The three atmospheric correction methods were
applied to Landsat OLI data path/row126/62 for 3 particular dates. Then, sample on vegetation, soil
and bodies of water (waterbody) were retrieved from the image. Atmospheric correction results were
visually observed and compared with SR sample on the absolute value, object spectral patterns, as
well as location and time consistency. Visual observation indicates that there was a contrast change
on images that had been corrected by using FLAASH method compared to SR, which mean that the
atmospheric correction method was quite effective. Analysis on the object spectral pattern, soil,
vegetation and waterbody of images corrected by using FLAASH method showed that it was not good
enough eventhough the reflectant value differed greatly to SR image. This might be caused by certain
variables of aerosol and atmospheric models used in Indonesia. QUAC and DOS made more
appropriate spectral pattern of vegetation and water body than spectral library. In terms of average
value and deviation difference, spectral patterns of soil corrected by using DOS was more compatible
than QUAC.
dioxide, oxygen, ozone, as well as light historical data, and suggest DOS-Iteration
effect and transmittance on the since the method gathered broader
atmosphere (Zhang et al., 2010). spectral range. Lee et al. (2000) who
To eliminate noise which might cause examined Landsat 5 Data with FLAASH,
image distortion, sun correction and Foster and 6S models found that FLAASH
atmospheric correction needed to be done. was the best model. Nguyen et al (2015)
In relation to position of the sun. Sun who reviewed QUAC, FLAASH and 6S
correction was performed by converting models with Landsat 5 Data stated that
digital number (DN) to reflectance value. 6S was the best model which resembled
Then Atmospheric correction was done to on ground condition. Afterall, this
reduce or eliminate atmospheric effect and research aimed to analyze a suitable
to obtain reflectant value on the suface. atmospheric correction model for Landsat
Some atmospheric correction methods OLI data in order to retrieve Indonesian
were, for example, digital number (DN) land cover information which
used directly for transformation reflection corresponded to actual coverage.
by Smith et al (1999), Dark Object To achieve the above result, an
Subtraction (DOS) by (Zhang et al., 2010; analysis was performed on some
Trisakti et al., 2014), Fast Line-of-sight atmospheric correction methods, i.e.
Atmospheric Analysis of Spectral FLAASH, QUAC and DOS for Landsat OLI
Hypercubes (FLAASH), as well as Quick path/row 126/62 in certain area of study
Atmospheric Correction (QUAC), 6S, and compared to USGS’s Surface Reflectance
ATCOR 2-3. data which had corrected. This research
This research utilized atmospheric was aimed to find the most accurate and
correction method which compatible with consistent method in accordance with
data quality, data availability, and actual condition.
research purposes. Some previous studies
on the correction method had been 2. MATERIALS AND METHODOLOGY
performed. Such as Somdatta et al. (2011) 2.1 Data
who implemented both FLAASH and Remote sensing data used in this
QUACC methods on hyperion data. They study were Landsat 8 level LT1 which was
conclude FLAASH was better than QUAC recorded by LAPAN Remote Sensing
in terms of atmospheric correction result. Center, Parepare and Rumpin with 3
Guo et al. (2012) who studied the spot different dates and Surface Reflectance
data with both methods found that data from USGS as the comparative data
FLAASH was more effective to reduce (Table 2-1).
noise than QUAC. Furthermore, Nazeer et
al. (2001) who implemented 6S, FLAASH, Table 2-1: Satelite images used in this research
ATCOR, DOS, and ELM methods on
Landsat 7 Data and compared the results Num Data Explanation Path/Row
with Surface Reflectance (SR) data from ber Acquisition
USGS, found that DOS and 6S generated 1 June 21, 2014 Landsat 8 126/62
the most appropriate and consistent 2 July 26, 2015 Landsat 8 126/62
result. Cui et al. (2001) also conducted a 3 May 9, 2016 Landsat 8 126/62
study on Landsat TM 1, Landsat TM 2, 4 July 26, 2015 USGS 126/62
Landsat TM 3, Landsat TM 4, Landsat TM Surface
5, and Landsat TM 7 with DOS, Reflectance
DOSCOST, DOS-Constant, and DOS-
Iteration. At last, Cui et al choose DOS as For spectral analysis purpose,
the best model to acquire quantitative and sampling process was carried out on
106 International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016
Comparing Atmospheric Correction Methods for Landsat OLI Data
vegetation (forests), water and soil. The surface temperature and physical surface
data used were acquired in May, June, reflection.
July or during dry season with 1 to 7 FLAASH began with a standard
spectral bands. The spectral bands were equation of spectral glow from single pixel
choosen since they were often used in received by standard Lambertian planar
classifying land cover. The area analyzed (the closest Lambertian planar), which
was Lake Kerinci surroundings, which was accepted by a sensor based on solar
included the lake itself, sea, agricultural spectrum (excluding thermal radiation),
lands, and soil. with the following formula:
Sample observed in this study were
vegetation (forest), body of water and soil, (2-1)
each was selected from three sample
locations (Figure 2-1). Sample taken were
then compared to Surface Reflectance L was the light for single pixel received by
data. Reference value of Landsat OLI spectral sensor; ρ was pixel surface reflectance; ρe
for vegetation, water and settlement were was average surface reflectance for the
shown on Figure 2-2 below. pixel and its surrounding; S was the
accumulation of sunlight reflection and
diffusion by atmospheric particles; La was
the radiance when atmospheric radiation
penetrated sensor. A and B were
coefficient determined by atmospheric and
geometric conditions of the underlying
surface but had nothing to do with
surface reflectance. (Aρ / (1-ρe S))
Figure 2-1: Sampling on water (a), soil (b),
represented radiation energy which
vegetation (c) on Landsat imagery penetrated directly into the sensor from
the target’s surface, which indicated two
cases: reflection occurs when the sun
illuminates target surface; or the
surrounding surface reflected through the
atmosphere and shined in the target
surface target before another reflection.
(Bρe / (1-ρe S)) showed total radiation that
went into sensor from the surface through
the atmosphere. ρ and ρe explained
"proximity effect" of mixed pixel near the
radiation caused by atmospheric
scattering, with the assumption ρ = ρe.
Figure 2-2: Landsat 8 object profile of water, soil However, significant error might occur due
and vegetation Source: Journal
Optics Express, 2014 to fog or strong contrast on the surfaces.
According to Equation 1, surface
2.2 Methods reflection could be calculated pixel by
2.2.1 FLAASH Model pixel. FLAASH used average radian
FLAASH was an atmospheric spatial, ignoring "proximity effect", to get
correction tool to fix visible band, NIR, estimate equation (2) and to predict
SWIR. This model eliminated air and light spatial reflectance. While, Le was the
influences, as well as removed interference average radiant image generated from
due to reflectivity parameter, emissivity, imagery radian and spatial function.
International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016 107
Esthi Kurnia Dewi and Bambang Trisakti
108 International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016
Comparing Atmospheric Correction Methods for Landsat OLI Data
2.2.3 DOS Model dark gray object pixel from grey value of
DOS was an image-based each pixel on the image.
atmospheric correction. Chavez (1996)
stated a basic assumption that few pixels
radiances on covered by cloud image "# = "min + "%& % (2-6)
could be accepted by the satellite due to
atmospheric scattering (path radiance).
Considering the fact that there was very Which:
few target on the Earth surface was Lmin = the light that corresponded to
absolute black, it was assumed that one digital value which was the total
percent of minimum reflektansi was better of all pixel which digital sum
than zero percent. was lower or min DN value.
Landsat 8 data which was LD01% = radiance of dark object which
radiometrically corrected including ToA was assumed to have 0.01
Reflectance dan sunlight correction. ToA reflectance value
Reflectance was corrected by converting
the DN value to reflectance value. Based To calculate Landsat Imagery
on USGS (2014), ToA Reflectance equation
was:
"() = *" ∗ %,() + " (2-7)
'=M + (2-4)
To calculate Object Radiance
Which:
' = ToA Reflectance, without correction
(2-8)
of sunlight angle.
Mρ = REFLECTANCE_MULT_BAND_x, in
which x is Band number So to calculate path radiance of the dark
= REFLECTANCE_ADD_BAND_x, in object
which x is Band number
Qcal=Digital Number (DN) value
(2-9)
Then, image correction was conducted
due to sunlight angle change, to eliminate Assumed that
DN value difference, with the following Tv =1
equation: Tz =1
Edown = 0
(2-5)
= '/(cos(θSZ))= '/(sin(θSE))
Untuk mencari ESUN Landsat 8 OLI
Which: adalah sebagai berikut
! = ToA Reflectance
-./,
SE = sun elevation = (0 ∗ 1 )
SZ = zenith angle of the sun, where SZ = 2 1) 3_* 5)(6( (2-10)
90 ° - SE ∗
237 3 8 3_* 5)(6(
International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016 109
Esthi Kurnia Dewi and Bambang Trisakti
110 International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016
Comparing Atmospheric Correction Methods
M for Landsat OLI Data
Band (-)
Figure 3-5: Visual Analysis of Surface Reflectance of Water in Image July 16, 2015
112 International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016
Comparing Atmospheric Correction Methods for Landsat OLI Data
Figure 3-7: Visual Analysis of Surface Reflectance of Soil in Image July 16, 2015
Deviation differences to SR (-)
Average differences to SR (-)
International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016 113
Esthi Kurnia Dewi and Bambang Trisakti
Reflectance (-)
Band (-)
Figure 3-9: Visual Analysis of Surface Reflectance of Vegetation in Image on July 26, 2015
114 International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016
Comparing Atmospheric Correction Methods for Landsat OLI Data
pattern, but higher absolute value (Figure Inter (spatial) location observation
3-15). Meanwhile, when spectral profile of here was different to research conducted
vegetation (forest) in 3 different sample by Bongetal (2015) which found that
locations were reviewed, all had almost FLAASH might generate better results
identical value to SR; which mean that than QUAC. This might happen since the
sampling locations were correct and type of aerosol variable and atmospheric
consistent (Figure 3-16). conditions used in Indonesia differed in
previous research areas.
Reflectance (-)
Band (-)
Figure 3-11: Average Surface Reflectance of Water in Three Sample Locations
1 1
2 2
Reflectance (-)
Reflectance (-)
3 3
Band (-)
Figure 3-13: Average Surface Reflectance of Soil in Three Sampling Locations
International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016 115
Esthi Kurnia Dewi and Bambang Trisakti
1 1
Reflectance (-)
Reflectance (-)
2
3 3
Band (-)
Figure 3-15: Average Surface Reflectance of Vegetation in Three Sampling Locations
1 1
Reflectance (-)
2
Reflectance (-)
2
3 3
116 International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016
Comparing Atmospheric Correction Methods for Landsat OLI Data
3.2.3 Temporal spectral curve almost the same and the spectral pattern
The temporal analysis of images was in accordance with the spectral
corrected using DOS method showed that library, which mean that temporal
absolute value of water, settlement (soil), sampling was consistent and the value
and vegetation (forest) sampling were was correct (Figure 3-17, 3-18, 3-19).
2014
Reflectance (-)
2015
2016
Band (-)
Figure 3-17: Water Spectral Profiles using DOS in Different Years
2014
Reflectance (-)
2015
2016
Band (-)
Figure 3-18: Soil Spectral Profiles using DOS in Different Years
International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016 117
Esthi Kurnia Dewi and Bambang Trisakti
Reflectance (-)
2014
2015
2016
Band (-)
Figure 3-19: Vegetation Spectral Profiles using DOS in Different Years
118 International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016
Comparing Atmospheric Correction Methods for Landsat OLI Data
Land cover data from Landsat single-date Nazeer M., Nichol JE, Yung Y., (2014), Evaluation
imagery: an approach integrating pixel- of atmospheric correction models and
based and object-based classifiers. Landsat surface reflectance product in an
European Journal of Remote Sensing 46 : urban coastal environment. International
699 – 717. Journal of Remote Sensing 35(16): 6271 –
Chavez PS, (1988), An improved dark-object 6291. doi: 10.1080/01431161.2014. 951742.
subtraction technique for atmospheric Nguyen HC, Jung J., Lee J., Choi S., Hong S.,
scattering correction of multispectral data. Heo J., (2015), Optimal Atmospheric
Remote sensing of Environment 24 (3): Correction for Above-Ground Forest
459-479. doi: http://dx.doi.org/10.1016/ Biomass Estimation with the ETM+ Remote
0034-4257(88)90019-. Sensor. Sensors 15(8): 18865–18886. doi:
Chavez PS, (1996), Image-based atmospheric 10.3390/s150818865.
corrections-revisited and improved. Rahayu, Candra DS, (2014), Koreksi Radiometrik
Photogrammetric engineering and remote Citra Landsat-8 Kanal Multispektral
sensing 62 (9): 1025-1036. Available at: Menggunakan Top of Atmosphere (Toa)
http://www.unc.edu/courses/2008spring/ untuk Mendukung Klasifikasi Penutup
geog/577/001/www/Chavez96-PERS.pdf/ Lahan. Seminar Nasional Penginderaan
(lastaccessed: 24/10/2016). Jauh.
Cui L., Li G., Ren H., He L., Liao H., Ouyang N., Smith GM, Milton EJ, (1999), The use of the
Zhang Y., (2014), Assessment of empirical line method to calibrate remotely
atmospheric correction methods for sensed data to reflectance. International
historical Landsat TM images in the coastal Journal of Remote Sensing 20 (13): 2653-
zone: A case study in Jiangsu, China. 2662,doi: http://dx.doi.org/10.1080/0143
European Journal of Remote Sensing 47: 11699211994.
701-716. doi: 10.5721/EuJRS20144740 Somdatta C., Chakrabarti S., (2011), Pre-
Ginting AY, Latifah S., Rahmawaty, (2012), processing of Hyperspectral Data: A case
Analisis Perubahan Tutupan Lahan study of Henry and Lothian Islands in
Kabupaten Karo (Analysis of Karo Regency Sunderban Region, West Bengal, India.
Land Cover Changes). Peronema Forestry International Journal of Geomatics And
Science Journal 1(1). Geosciences 2(2).
Guo Y., Zeng F., (2012), Atmospheric Correction Tampubolon T., Jeddah Y., (2015), Aplikasi
Comparison of Spot-5 Image Based on Pemanfaatan Citra Satelit Landsat untuk
Model Flaash and Model Quac. Mengidentifikasi Perubahan Lahan Kritis di
International Archives of the Photogrammetry, Kota Medan danSekitarnya. Spektra:
Remote Sensing and Spatial Information Jurnal Fisika dan Aplikasinya 16(2): 15-19.
Sciences, XXXIX-B7: 7-11 Trisakti B., Suwarnana N., Cahyono JS, (2014),
Lee SB, La HP, Eo YD, Pyeon MW, (2015), Pemanfaatan Data Penginderan Jauh
Generation of Simulated Image from untuk Memantau Parameter Status
Atmospheric Corrected Landsat TM Images. Ekosistem Perairan Danau (Studi Kasus:
Journal of the Korean Society of Surveying Danau Rawa Pening). Seminar Nasional
Geodesy Photogrammetry and Cartography Penginderaan Jauh.
33(1): 1-9. Tuni MS, Barus B., Iskandar, (2013), Prediksi
Moran E., (2002), Assessment of atmospheric Perubahan Tutupan Lahan dan
correction methods for Landsat TM data Perencanaan Penggunaan Lahan
applicable to Amazon basin LBA research. Pascatambang Nikel di Kabupaten
ACT Publication No 02-06. Halmahera Timur. Globe 15(2): 146 – 152.
International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016 119
Esthi Kurnia Dewi and Bambang Trisakti
Vincent RK, (1972), An ERTS multispectral Sistem Informasi Geografis. Jurnal Ilmu
scanner experiment for mapping iron dan Teknologi KelautanTropis, 6 (2): 311-
compounds. Proceedings of the Eighth 318.
International Symposium on Remote Zhang X., Yang H., Shuai T., (2010). Comparision
Sensing of Environment, Ann Arbor, of FLAASH versus Empirical Line Approach
Michigan, II: 1239-1247. Available at: for Atmospheric Correction of OMIS-II
http://ntrs.nasa.gov/archive/nasa/casi.nt Imagery. Journal Chinese Academy of
rs.nasa.gov/19730001633.pdf/(last accessed: Sciences, Beijing.
24/09/2014). Zhang Z., He G., Wang X., (2010), A practical
Yulius, Tanto TA, Ramdhan M., Putra A., Salim DOS model-based atmospheric correction
HL, (2014), Perubahan Tutupan Lahan di algorithm. International Journal of Remote
Pesisir Bungus Teluk Kabung, Sumatra Sensing 31(11): 2837-2852, doi: http://
Barat Tahun 2003-2013 Menggunakan dx.doi.org/10.1080/01431160903124682.
120 International Journal of Remote Sensing and Earth Science Vol. 13 No. 2 December 2016