Comparison of on-site AOT and AOT from GIS algorithm (4/2010–5/2011).
Abstract
This chapter focuses on the development of an image-based integrated method for determining and mapping aerosol optical thickness (AOT). Using the radiative transfer (RT) equation, a methodology is developed using nonvariant targets with the revised darkest pixel method to improve the accuracy to calculate AOT over urban areas. This method can be applied to create thematic maps using GIS that present AOT values in a given area. The methodology is applied to Landsat TM/ETM+ satellite images of Limassol, Cyprus over a period of time to visualize and assess the AOT levels over the urban area. An accuracy assessment of the method shows a strong correlation between the in-situ AOT values from the sunphotometers and the AOT values derived from the image-based integrated method.
Keywords
- AOT
- GIS
- atmospheric correction
- image-based mapping
- Landsat
- Kriging
1. Introduction
Aerosol optical thickness (AOT) is defined as the degree to which aerosols prevent the transmission of light through absorption or scattering of light. AOT can be retrieved by using radiometers on board satellites [1–3], ground-based sunphotometers, and satellite sensors, such as the moderate resolution imaging spectroradiometer(MODIS), which provide AOT measurements on a daily and monthly basis [4, 5]. However, the spatial resolution of the MODIS AOT data products is 10 × 10 km, which does not permit the identification of specific AOT distribution over urban areas or complex terrains [5–7]. Due to the variability of aerosols, atmospheric aerosol monitoring is difficult. Significant efforts to improve aerosol characterizations have included using
Research indicates that the determined AOT from satellite image data can be used as a tool to assess air pollution as well as identify the sources of local emissions [2, 7, 11–24]. AOT values can be used as a way of measuring air quality; determining AOT in large-scale pollution areas provides a synoptic, cost-effective means to further assess the air quality in such areas [16, 18–24]. Indeed, the AOT values derived from the atmospheric path radiance can be utilized to assess and monitor air quality and atmospheric pollution [8]. The image-based integrated method presented in this chapter can accurately calculate AOT values retrieved from satellite imagery by using radiative transfer (RT) equations and GIS. The method can be used to visually display AOT levels using thematic maps in order to identify concentrations of AOT over an urban area [11]. The image-based method can also be used with archived satellite images, thereby providing detailed information regarding spatial aerosol concentration overtime [11].
AOT is directly related to the atmospheric aerosol load, which is the main variable describing the effects of aerosols on radiative transfer in the Earth’s atmosphere. According to Guanter et al. [25], modeling atmospheric constituents and surface reflectance involves modeling the radiative transfer across the atmosphere. The key parameter for assessing atmospheric pollution is the aerosol optical thickness, which is also the most important unknown of every atmospheric correction algorithm for solving the radiative transfer equation and removing atmospheric effects from remotely sensed satellite images [8, 26–29]. Several researchers [5, 18, 26, 28, 29–37] found that using the radiative transfer and atmospheric modeling in conjunction with field measurements of aerosol optical thickness can yield more accurate atmospheric corrections instead of using simple image-based techniques.
The image-based integrated method discussed in this chapter combines the RT equation, satellite imagery, the modified darkest pixel (DP) method of atmospheric correction, and GIS to derive AOT measures. An example including 11 Landsat satellite images with
2. Methodology
The innovation of this methodology is the retrieval of AOT by using an image-based algorithm based on the radiative transfer equation to derive AOT using the reflectance values from the darkest pixel method of atmospheric correction. In the image-based integrated method, the RT equation retrieves AOT values from the satellite images with as much precision as possible. Also, the integrated use of field spectroscopy, GIS, and remote sensing for AOT retrieval is an important contribution to the field of remote sensing.
The first step of the methodology is the preprocessing of the satellite images, including geometric and radiometric correction. Also, the satellite images need to be atmospherically corrected. The darkest pixel atmospheric correction method can be applied, as it provides reasonable atmospheric correction, especially in cloud-free skies [28, 38, 39]. If there is a known nonvariant target or overpass
2.1. Detailed description of the image-based integrated method for determining AOT
According to Hadjimitsis and Clayton [40], RT equations can be used to retrieve AOT values from the path radiance. Therefore, the algorithm used in the image-based integrated method is formulated to solve for AOT, which is included in the equation for path radiance (
Path radiance is the sum of the atmospheric path radiance due to Rayleigh (
where
2.1.1. Calculating path radiance
The path radiance “‘
where
Eq. (3) can be used to calculate the at-satellite radiance (
The “‘
where
The “‘
Although several equations to calculate
where
2.1.2. Calculating Rayleigh scattering (L pr )
To calculate path radiance due to Rayleigh scattering(
where
The Rayleigh scattering phase function (
2.1.3. Calculating Mie scattering
To calculate Mie scattering, Eq. (9) is used [47]:
where
2.1.4. Steps for performing the image-based integrated method
In order to retrieve AOT using the above equation and solve for
2.2. Image-based integrated method for GIS analysis
In order to solve for AOT and conduct a GIS analysis, it may be necessary to simplify the equation due to the complexity of the RT equations and their logarithmic components. The equations for path radiance, Rayleigh scattering, and Mie scattering can be simplified and then combined into one equation, to be solved in an image-based software. After the AOT values are calculated for every pixel in the satellite image, the values can be exported into a GIS geospatial database, with the AOT value and coordinates. In order to display AOT on a GIS thematic map, the Kriging method can be used to interpolate the AOT values and then overlay the data into GIS vector maps or satellite images of the area of interest to display a spatial distribution of AOT. The maps can also be overlaid with vector data including roads, building lots, parks, green areas, stadiums, industrial areas, and municipal boundaries to identify AOT distribution over the area of interest.
3. AOT retrieval using the image-based integrated method: the example of Limassol, Cyprus
The example described in this chapter includes AOT and surface reflectance field measurements in Limassol Cyprus during Landsat-5 TM and Landsat-7 ETM+ satellite overpass. Only 11 Landsat satellite images with low cloud cover are used. These images were chosen due to their availability and clarity. The Landsat-5 TM images used are for 24/6/2010, 10/7/2010, 27/8/2010, and 28/9/2010, and the Landsat-7 ETM+ images used are for 13/4/2010, 29/4/2010, 31/5/2010, 16/6/2010, 7/11/2010, 9/12/2010, and 2/5/2011. During each satellite overpass, AOT measurements were taken using the Microtops II sunphotometer at five different locations in Limassol and the stationery Cimel sunphotometer, which is part of the AERONET network, located in the center of the study area. These AOT values are later compared to the AOT values derived from the algorithm using the satellite images. At the same time, the SVC HR-1024 spectroradiometer was used to measure the reflectance value of a large nonvariant target dark asphalt surface located near the center of the satellite image.
Geometric and radiometric correction is performed on the 11 Landsat-5 and Landsat-7 satellite images, as well as atmospheric correction using the modified DP method. The image-based integrated method is applied to the Landsat-5 and Landsat-7 images taken over specific days of 2010–2011. Eqs.(1–9) are used to retrieve AOT from bands 1–4 of the Landsat-5 and Landsat-7 satellite images, following the steps outlined. As the
In this example, the model is run in the ERDAS Imagine and MATLAB software to generate an image that is consists of AOT values. The model is then applied to every surface using the Landsat band 1 satellite image, since the AOT required are in the 500 nm wavelength. The resulting image features the AOT values for every pixel in the image, as featured in Figure 2. Areas with no data, due to cloud presence or negative AOT values, appear in white. All AOT values can be exported to ASCII files to create a georeferenced AOT dataset, which is then imported into the ArcGIS software.
Using the ArcGIS software, each pixel is converted to a point and each point is associated with the calculated AOT value. The white area of the image contains values of “no data” (Figure 3, left). To create a GIS thematic map showing the AOT distribution in the urban area of Limassol, the interpolation method is used in order to estimate the values from the “no data” regions. The ordinary Kriging interpolation tool is applied, where the unknown values of the “no data” area are interpolated using the weighted average of neighboring samples. Figure 3 (right) shows the results of the Kriging interpolation for the Limassol example, where the dots indicate the points with AOT values, while the colored sections are the resulting AOT values determined from Kriging interpolation.
Once the interpolation is completed, an AOT thematic map is created, which is classified according to AOT values. The AOT values are displayed in different colors according to the AOT concentration and a legend of the AOT values is created. This facilitates the ability to visualize the concentration and distribution of AOT values over an urban area. Figure 4 features a map of Limassol following Kriging interpolation, which provides synoptic views of the AOT distribution. As is evident, high AOT levels are present throughout the city, especially in the industrial estates and on busy streets, with lower AOT levels present in parks, stadiums, and outside the city.
In the example, the thematic map is overlaid with GIS vector data from the Lands and Surveys Department for Limassol to facilitate the identification of sources of AOT values within the area of interest (Figure 5). Each polygon in the thematic map is connected to a GIS database, where data such as plot number and area information can be identified. In this way, the GIS database, combined with the AOT data, can create thematic maps to illustrate trends in different areas.
In the example, the image-based integrated method is used to produce GIS maps indicating the AOT distribution over Limassol. In order to visualize high and low AOT levels, thematic maps are generated using colors ranging from blue to green to yellow to red for the specified AOT range. A separate thematic map is created for each of the 11 Landsat-5 and Landsat-7 satellite images. Figure 6 features the GIS thematic map using Landsat-7 satellite imagery from 24/6/2010.
The AOT values derived from the algorithm are compared with the AOT values measured with the Microtops and Cimel sun photometers. Figure 7 shows all the locations of the
Table 1 compares the AOT measurements taken on site using a sunphotometer against the AOT values derived from the GIS maps using the proposed algorithm by location and date, as well as the correlation coefficient (
4. Conclusions
In this chapter, the image-based integrated method is developed to retrieve AOT values from satellite images and display these values in an urban area, using the example of Limassol, Cyprus. The findings show that the proposed integrated image-based method is able to accurately obtain AOT measurements from satellite images. The accuracy assessment shows strong agreement between the
Acknowledgments
The author would like to thank the Cyprus University of Technology, Department of Civil Engineering and Geomatics, Remote Sensing Lab for use of the sunphotometers and AERONET data.
References
- 1.
Anderson, T., Charlson, R., Schwartz, S., Knutti, R., Boucher, O., Rodhe, H. & Heintzenberg, J. Climate forcing by aerosol—a hazy picture. Science . 2003;300 :1103–1104. - 2.
Kaufman, Y.J., Tanre, D. & Boucher, O. A satellite view of aerosols in the climate system. Nature . 2002;419 :215–223. DOI: 10.1038/nature 01091 - 3.
Holben, B.N., Smirnov, A., Eck, T.F., Slutsker, I., Abuhassan, N., Newcomb, W., Schafer, J.S., Tanre, D., Chatenet, B. & Lavenu, F. An emerging ground-based aerosol climatology: aerosol optical depth from AERONET. Journal of Geophysical Research . 2001;106 :12067–12098. - 4.
Yu, H., Dickinson, R.E., Chin, M., Kaufman, Y.J., Zhou, L., Tian, Y., Dubovik, O. & Holben, B.N. Direct radiative effect of aerosols as determined from a combination of MODIS retrievals and GOCART simulations. Journal of Geophysical Research . 2004;109 :1–13(D03206). DOI: 10.1029/2003JD003914 - 5.
Remer, L., Kaufman, Y.J. & Tanre, D. The MODIS aerosol algorithm, products, and validation. Journal of the Atmospheric Scien ces. 2005;62 (4):947–973. - 6.
Misra, A., Jayaraman, A. & Ganguly, D. Validation of MODIS derived aerosol optical depth over western India. Journal of Geophysical Research . 2008;113 :1–12 DOI: 10.1029/2007JD009075 - 7.
Guleria, R.P., Kuniyal, J.C., Rawat, P.S., Thakur, H.K., Sharma, M., Sharma, N.L., Singh, M., Chand, K., Sharma, P., Thakur, A.K., Dhyani, P.P. & Bhuyan, P.K. Aerosols optical properties in dynamic atmosphere in the northwestern part of the Indian Himalaya: a comparative study from ground and satellite based observations. AtmosphericResearch . 2011;101 (3):726–738. - 8.
Retalis, A., Cartalis, C. & Athanassiou, E. Assessment of the distribution of aerosols in the area of Athens with the use of LANDSAT Thematic Mapper data . International Journal of Remote Sensing . 1999;20 (5):939–945. - 9.
Hadjimitsis, D.G. Aerosol optical thickness (AOT) retrieval over land using satellite image-based algorithm. Air Quality Atmosphere and Health-An International Journal . 2009;2 (2):89–97. - 10.
Huang, D., Knyazikhin, Y., Wang, W., Deering, D.W. & Stenberg, P. Stochastic transport theory for investigating the three-dimensional canopy structure from space measurements. Remote Sensing of Environment . 2008;112 :35–50. - 11.
Themistocleous, K. & Hadjimitsis, D.G. Development of an image-based integrated method for determining andmapping aerosol optical thickness (AOT) over urban areas using the darkestpixel atmospheric correction method, RT equation and GIS: a case study ofthe Limassol area in Cyprus. ISPRS Journal of Photogrammetry and Remote Sensing . 2013;86 :1–10. DOI: 10.1016/j.isprsjprs.2013.08.011 - 12.
Zawadzka, O., Markowicz, K.M., Pietruczuk, A., Zielinski, T. & Jaroslawski, J. Impact of urban pollution emitted in Warsaw on aerosol properties. Atmospheric Environment . 2013;69 :15–28. - 13.
Yahi, H., Santer, R., Weill, A., Crepon, M. & Thiria, S. Exploratory study for estimating atmospheric low level particle pollution based on vertical integrated optical measurements. Atmospheric Environment . 2011;45 (23):3891–3902. - 14.
Wang, Z., Chen, L., Tao, J., Zhang, Y. & Su, L. Satellite-based estimation of regional particulate matter (PM) in Beijing using vertical-and-RH correcting method. Remote Sensing of Environment . 2010;114 (15):50–63. - 15.
Jiang, X., Liu, Y., Yu, B. & Jiang, M. Comparison of MISR aerosol optical thickness with AERONET measurements in Beijing metropolitan area. Remote Sensing of Environment . 2007;107 :45–53. - 16.
Hadjimitsis, D.G. & Clayton, C.R.I. Detecting air pollution from space using image-based method. In: International Conference: Protection and Restoration of the Environment , VIII, 313; Crete, Greece. 2006. - 17.
Koukouli, M.E., Balis, D.S., Amiridis, V., Kazadzis, S., Bais, A., Nickovic, S. & Torres, O. Aerosol variability over Thessaloniki using ground based remote sensing observations and the TOMS aerosol index. Atmospheric Environment . 2006;40 (28):5367–5378. - 18.
Hadjimitsis, D.G., Retalis, A. & Clayton, C.R.I. The assessment of atmospheric pollution using satellite remote sensing techniques in large cities in the vicinity of airports. Water, Air and Soil Pollution: Focus, An International Journal of Environmental Pollution . 2002;2 (5–6):631–640. - 19.
Tang, J., Xue, T., Yuc, T. & Guen, Y.N. Aerosol optical thickness determination by exploiting the synergy of TERRRA and AQUA MODIS. Remote Sensing and Environment . 2005;94 :327–334. - 20.
Kanaroglou, P.S., Soulakellis, N.A. & Sifakis, N.I. Improvement of satellite derived pollution maps with the use of a geostatistical interpolation method. Journal of Geographical Systems . 2002;4 :193–208. - 21.
Prospero, J.M., Ginoux, P., Torres, O., Nicholson, S. & Gill, T. Environmental characterization of global sources of atmospheric soil dust identified with the NIMBUS7 total ozone mapping spectrometer (TOMS) absorbing aerosol product. Review of Geophysics . 2002;40 (1):1002. DOI: 10.1029/2000RG000095 - 22.
Sifakis, N. & Deschamps, P.Y. Mapping of air pollution using SPOT satellite data. Photogrammetric Engineering and Remote Sensing . 1992;58 :1433–1437. - 23.
Leon, J., Chazette, P. & Dulac, F. Direct and indirect methods for correcting the aerosol effect on remote sensing. Remote Sensing and Environment . 1995;55 :65–79. - 24.
Sifakis, N., Soulakellis, N. & Paronis, D. Quantitative mapping of air pollution density using Earth observations: a new processing method and application to an urban area. International Journal of Remote Sensing . 1998;19 (17):3289–3300. - 25.
Guanter, L., Gómez-Chova, L. & Moreno, J. Coupled retrieval of aerosol optical thickness, columnar water vapor and surface reflectance maps from ENVISAT/MERIS data over land. Remote Sensing of Environment . 2008;112 :2898–2913. - 26.
Hadjimitsis, D.G., Themistocleous, K., Vryonides, P., Toulios, L. & Clayton, C.R.I. Applications of satellite remote sensing and GIS to urban air-quality monitoring: potential solutions and suggestions for the Cyprus area. In: 6th International Conference on Urban Air Quality , 144; Limassol, Cyprus. 2007. - 27.
Hadjimitsis, D.G. & Clayton, C.R.I. Determination of the aerosol optical thickness and assessment of atmospheric conditions using satellite image-based processing algorithm and radiative transfer theory. In: 7th Pan-Hellenic International Conference of Meteorology, Climatology and Atmospheric Physics ; University of Cyprus, Nicosia, Cyprus. 2004. - 28.
Hadjimitsis, D.G., Clayton, C.R.I. & Hope, V.S. An assessment of the effectiveness of atmospheric correction algorithms through the remote sensing of some reservoirs. International Journal of Remote Sensing . 2004;25 (18):3651–3674. - 29.
Forster, B.C. Derivation of atmospheric correction procedures for Landsat MSS with particular reference to urban data. International Journal of Remote Sensing . 1984;5 (5):799–817. - 30.
Istomina, L.G., von Hoyningen-Huene, W., Kokhanovsky, A.A., Schultz, E. & Burrows, J.P. Remote sensing of aerosols over snow using infrared AATSR observation. Atmospheric Measurement Techniques . 2011;4 :1133–1145. DOI: 10.51941.amt-4-1133-2011 - 31.
Lee, K.H. & Kim, Y.J. Satellite remote sensing of Asian aerosols: a case study of clean, polluted and dust storm days. Atmospheric Measurement Techniques Discussions . 2010;3 :2651–2680. DOI: 10.5194/amtd-3-2651-2010 - 32.
Xue, Y., Wan W., Li, Y.J., Guang, J., Bai, L.Y., Wang Y. & Ai, J.W. A data intensive scientific computing framework for quantitative retrieval of geophysical parameters using satellite data. IEEE Computer . 2008;41 (4):33–40. - 33.
Seidel, F., Schlapfer, D., Nieke, J. & Itten, K.I. Sensor performance requirements for the retrieval of atmospheric aerosols by airborne optical remote sensing. Sensors . 2008;8 :1901–1914. - 34.
Hadjimitsis, D.G. The application of atmospheric correction algorithms in the satellite remote sensing of reservoirs [dissertation]. Guildford, U.K.: University of Surrey, School of Engineering in the Environment, Department of Civil Engineering; 1999. - 35.
Kaufman, Y.J. & Sendra, C. Algorithm for automatic atmospheric corrections to visible and near-IR satellite imagery. International Journal of Remote Sensing . 1988;9 (8):1357–1381. - 36.
Kaufman, Y.J. & Tanre, D. Direct and indirect methods for correcting the aerosol effect on remote sensing . Remote Sensing of Environment . 1996;55 :65–79. - 37.
Kaufman, Y.J., Tanre, D., Remer, L.A., Vermote, E.F., Chu, A. & Holben, B. N. Operational remote sensing of tropospheric aerosol over land from EOS moderate resolution imaging spectroradiometer. Journal of Geophysical Research . 1997;102 :17051–17067. - 38.
Schott, J.R. Remote Sensing: The Image Chain Approach . Oxford, UK: Oxford University Press; 2007. - 39.
Cheng, K.S. & Lei, T.C. Reservoir trophic state evaluation using Landsat TM images. Journal of the American Water Resources Association . 2001;37 :1321–1334. - 40.
Hadjimitsis, D.G. & Clayton, C.R.I. Assessment of temporal variations of water quality in inland water bodies using atmospheric corrected remotely sensed image data. Environmental Monitoring and Assessment . 2009;159 (1–4):281–292. - 41.
deHaan, J.F., Hovenier, J.W., Kokke, J.M.M. & Van Stokkom, H.T.C. Removal of atmospheric influences on satellite-borne imagery: a radiative transfer approach. Remote Sensing of Environment . 1991;37 :1–21. - 42.
Hill, J. & Sturm, B. Radiometric correction of multi-temporal Thematic Mapper data for use in agricultural land cover classification and vegetation monitoring. International Journal of Remote Sensing . 1991;12 (7):1471–1491. - 43.
Deschamps, P.Y., Herman, M. & Tanré, D. Modeling of the atmospheric effects and its application to the remote sensing of ocean colour. Applied Optics . 1983a;23 :3751–3758. - 44.
Gilabert, M.A., Conese, C. & Maselli. An atmospheric correction method for the automatic retrieval of surface reflectances from TM image. International Journal of Remote Sensing . 1994;15 (10):2065–2086. - 45.
Turner, R.E. & Spencer, M.M. Atmospheric model for correction of spacecraft data. In: Eight International Symposium on Remote Sensing of the Environment ; 895–934; Ann Arbor, Michigan; 1972.. - 46.
Landsat-7 Science Data Users Handbook. National Aeronautics and Space Administration; 2011. Accessed at: http://landsathandbook.gsfc.nasa.gov/pdfs/Landsat7_ Handbook.pdf - 47.
Moller, F. Strahlung in der UnterenAtmosphare, Handbuch der Physik . New York: Springer-Verlag; 1957. - 48.
Waggoner, A.P., Weiss, R.E., Ahlquist, N.C., Covert, D.S., Wills, S. & Charlson, R.J. Optical characteristics of atmospheric aerosols. Atmospheric Environmen t. 1981;15 :1891–1909.