Remote Sensing (RS)
Hosein Nesari; Reza Shah-Hosseini; Amirreza Goodarzi; Soheil Sobhan Ardakani; Saeed Farzaneh
Abstract
Extended Abstract
Introduction
Atmospheric aerosols are a colloid of solid particles or liquid droplets suspended in the atmosphere. Their diameter is between 10-2 to 10-3 micrometers. They directly and indirectly affect the global climate by absorbing and scattering solar radiation, and they also ...
Read More
Extended Abstract
Introduction
Atmospheric aerosols are a colloid of solid particles or liquid droplets suspended in the atmosphere. Their diameter is between 10-2 to 10-3 micrometers. They directly and indirectly affect the global climate by absorbing and scattering solar radiation, and they also have a serious impact on human health by emitting harmful substances. In addition, high concentrations of aerosols on a local scale due to natural or human activities have adverse effects on human health, including cancers, pulmonary inflammation, and cardiopulmonary mortality. Monitoring the temporal and spatial variability of high concentrations of aerosols requires regular measurement of their optical properties such as aerosol optical depth (AOD).
Materials & Methods
Algeria is a large country with little knowledge of the spatial and temporal diversity of AOD, and the low spatial resolution of existing products makes it very difficult to predict aerosols (airborne particles) at the local scale, especially in arid southern regions. As a result, AOD recovery with data with higher spatial resolution is crucial for determining air pollution and air quality information. Several AERONET stations have been installed in Algeria. The Tamanrasset_INM station has been selected based on its location and the availability of historical AOD data for the period (2015-2016).
In this study, Landsat-8 / OLI image from tile 192/44 was used for satellite images. To this end, 23 TOA-corrected L1G-level Landsat-8 / OLI cloudless scenes were downloaded from January 2015 to December 2016 in the study area. DN values are converted to TOA reflections using the scaling factor coefficients in the OLI Landsat-8 metadata file. In this study, the minimum monthly reflectance technique was used to recover AOD in this area. As a result, LSR images were used in the recovery process in different months of 2015 and 2016. The process of selecting reference LSRs was initially based on the selection of clear, foggy / cloudless sky images. The selected images were then used to construct artificial images in which each pixel corresponds to the second lowest surface reflection of all selected monthly images to be the LSR pixel for the respective month. The AOD retrieval method developed in this study is based on a LUT, using the 6S radiative transfer model. The advantage of using the 6S model is its ability to estimate direct components and scattering using a limited number of inputs for each spectral band in the entire solar domain. The effect of the viewing angle is limited because Landsat data are usually obtained with a fixed viewing angle. Surface reflectance can be estimated from a pre-calculated LSR database. The accuracy of AOD recovery depends on the use of the appropriate aerosol model. A continental model was selected from the available aerosol models. Other atmospheric parameters such as ozone, carbon dioxide, carbon monoxide and water vapor are considered by default. The AOD values used to make LUT are set as follows: 0.0, 0.05, 0.1, 1.5, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.2 and 1.5. The zenith angles of the sun and the sensor range from 0 to 70 degrees with a step of 5 degrees and the range of azimuth angles from 0 to 180 degrees with a step of 12 degrees. Using these parameters, the radiative transfer equation was run in forward to obtain the TOA reflection. Different combinations of input and TOA output parameters are stored in LUT. AOD retrieval is based on a comparison between the TOAs estimated with the model and the observed items using the best fit approach. Using such an approach, the estimated AODs are simulated in accordance with those used in the production of TOAs, using a competency function that minimizes the distance.
Results & Discussion
In this study, the AODs recovered at 550 nm in a 5-by-5-pixel window around the AERONET site were averaged. The considered AERONET values are the average of all measurements taken within ± 30 minutes of image acquisition time. Observation regression results (AOD from Landsat 8 images and AERONET stations) showed that the correlation coefficient is about 84%. This study shows a good fit of the model on the research data and shows the high capability of the model. This study showed a strong recovery of AOD against AERONET data of more than 70% at . The differences can be attributed to a limited number of points or hypotheses related to the aerosol model used in this study. The assumption of using a pre-calculated LSR does not limit the accuracy of this method because we have shown that in arid regions where the change in land cover in different months of the year is small, a pre-calculated LSR image can be representation used the share of surface reflection in the radiative transfer model throughout the month.
Conclusion
In this study, an AOD derived from a high-resolution satellite at an urban scale was produced in the city of Tamanrasset, Algeria. The developed method assumes that the change in land cover is minimal and the temporal change in LSR is not significant. A pre-calculated LSR image is created to show the surface reflection in the retrieval process. Based on the 6S radiative transfer model, an LUT was constructed to simulate the TOA reflection of the built-in LSRs and a set of geometric and atmospheric parameters. The retrieved AODs were compared with the AERONET ground data. The results show that this approach can achieve reasonable accuracy in AOD recovery, which reaches about 70.9% at . In addition, this approach is suitable for estimating AOD in urban areas compared to existing AOD products with low spatial resolution. The results of this study show a 4% improvement compared to the results of Omari et al. (2019). The results of this study showed that ignoring the monthly changes in LSR values leads to good results in AOD recovery.
Remote Sensing (RS)
Somayeh Aslani Katouli; Reza Shah-Hosseini; Hamid Bagheri
Abstract
Extended Abstract
Introduction
A flood is a widespread and dramatic natural disaster that affects the life, infrastructure, economy, and local ecosystems of the world. In this paper, a method for flood detection in urban (and suburban) environments using the intensity and coherence of SAR based on ...
Read More
Extended Abstract
Introduction
A flood is a widespread and dramatic natural disaster that affects the life, infrastructure, economy, and local ecosystems of the world. In this paper, a method for flood detection in urban (and suburban) environments using the intensity and coherence of SAR based on a convolutional neural network is introduced, and from the time series of SAR intensity and coherence to draw flood without obstruction (e.g. Flooded bare soils and short vegetation) are used. Non-cohesive areas blocked by floods (e.g., flooded vegetation) and cohesive areas with flood-blocked areas (e.g., frequently constructed flooded areas) are distinguished.
This method is flexible according to the time period of the data sequences (at least one pair of pre-event and event intensities and one pair of pre-event and in-event coherence are required). The increasing number of SAR missions in orbit that have a fixed viewing scenario with a short retry time increases the chances of seeing a flood event, while also having a good pre-event scene achieved by the same sensor. This makes this method desirable for operational emergency responses.
Materials & Methods
CNN algorithm is a multilayer perceptron that is designed to identify two-dimensional information of images and includes: input layer, convolution layer, sample layer, and output layer. The CNN algorithm has two main processes: collection and sampling.
The convolution process involves the use of a trainable Fx filter, deconvolution of the input image (the first step of image input, input after image convolution, is the feature of each layer called Feature Map), then by adding bx can be hand convolution of the CX layer Found. Sampling process: n pixels are collected from each neighborhood to form a pixel, then weighted with a scalar weight of Wx + 1 and a bx + 1 bias is added, then a map of The Narrow n times feature map properties are generated.
Three images of Sentinel-1A VV polarization, wide width interference (IW), and mode (SLC) data were used in this study. Intensity images were pre-processed with radiometric calibration, noise reduced with a spell-filter (window size 5.5 pixels), and converted from linear units to decibels. Coherent images were obtained with a pair of consecutive images with a window of 7.28 (range _ azimuth). Validation data set due to the lack of other data in two separate sections of ground data in the urban area of GonbadKavous that have been collected to identify homes damaged by floods and terrestrial reality data from gamma image thresholds for output validation were extracted.
Results & Discussion
In this section, the results of the study are qualitatively and quantitatively analyzed. Because the simultaneous display of SAR data over time in the form of RGB compounds is widely used in the qualitative interpretation of land cover and surface dynamics, RGB compounds are used to provide evidence of flood magnitude in terms of intensity and coherence. For both cases, the results of combining intensity and coherence and intensity alone and coherence alone are quantitatively analyzed. Overall accuracy (OA), kappa correlation coefficient, false-positive rate (FPR), precision (e.g., correctly predicted positive patterns out of the total predicted patterns in a positive class), recall (e.g., a fraction of properly classified positive patterns), and an F1 score (ie the harmonic mean between precision and recall). Flood reference and ground data are mentioned and reported based on the reference.
Conclusion
In this paper, a method for mapping floods in urban environments based on SAR intensity and interferometry coherence was introduced. A combination of intensity and coherence extracts flood information in different types of land cover and outlet. This method was tested on the KavousGonbad flood incident obtained by various SAR sensors and the flood maps were confirmed by the flood reference resulting from thresholding and ground harvesting and satisfactory results were shown in this case study. The findings of this experiment show that the shared use of SAR intensity and coherence provides more reliable information than the use of SAR intensity and coherence alone in urban areas with different landscapes. In particular, flood detection in less cohesive / non-cohesive areas (e.g., bare soils, vegetation, vegetated areas) relies heavily on multi-temporality, while multi-temporal coherence provides more comprehensive flood information in areas Create coherence (e.g., mostly built-up areas). However, some flood-specific situations, such as flooded parking lots and flooded dense building blocks, are still challenging in terms of intensity and coherence. Also, since the proposed method is sensor and scene independent, with very frequent and regular observations of SAR missions such as Sentinel-1 and RADARSAT (RCM), there are opportunities to map global floods on a global scale, especially in small countries. Provides income.
Remote Sensing (RS)
Moslem Darvishi; Reza Shah-Hosseini
Abstract
Extended Abstract IntroductionWith the expansion of the urban limits, some of the lands that were used for gardening years ago have been located within the urban limits. The difference between the value of garden land use and urban land use, such as residential and commercial, encourages gardeners to ...
Read More
Extended Abstract IntroductionWith the expansion of the urban limits, some of the lands that were used for gardening years ago have been located within the urban limits. The difference between the value of garden land use and urban land use, such as residential and commercial, encourages gardeners to change their land use. Urban managers try to prevent this change of use by enforcing strict rules.Assessing the success of such plans requires examining land-use change in the urban over a long periodof time. The main purpose of this study is to detect abandoned urban gardens using Landsat satellite imagery. The second goal is to determine the extent of changes in urban gardens in the study area over the past 30 years. In this study, based on Landsat satellite images in 2018 and 1988 for the northern slope of Alvand Mountain in Hamadan province and the city of Hamadan, the normalized differential index of vegetation (NDVI) along with land surface temperature (LST) in 9 time periods per year was extracted. The results indicated a 4/75 ° C increase in LSTfor the region over 30 years. Also, the inverse relationship of LST with NDVI is confirmed. Based on the separation of urban gardens, a comparison was made between 2018 and 1988, which showed a decrease of 175 hectares of urban gardens in the study area, which is equivalent to a 49% reduction in urban gardens. In the main part of the research, based on the behavioral evaluation of urban gardens, in these two characteristics, a differentiation index for active and abandoned gardens is presented. Examination of the results based on ground truth data including 25 active gardens and 25 abandoned gardens suggested that the proposed method had an overall accuracy of 82% and a Kappa coefficient of 0/64.Materials & MethodsThe study area includes a part of the northern slope of Alvand Mountain, which is limited to the southern part of Hamedan and has a latitude of 34 degrees and 45 minutes to 34 degrees and 48 minutes north and a longitude of 48 degrees and 27 minutes to 48 degrees and 31 minutes east. Ground truth data including 25 active gardens and 25 abandoned gardens were collected as field visits using a Garmin GPSMAP 62s handheld navigator so that coordinates were collected by attending the location of abandoned and active gardens. The satellite data used in this study concern the time series data of Landsat 8 satellite OLI and TIRS sensors for 2018 and Landsat 5 satellite TM sensor for 1988.To achieve the first objective and separate active and abandoned gardens in 2018, the land surface temperature (LST) and the normalized difference vegetation index (NDVI) are calculated and the behavior pattern of these two components is examined during the year for active and abandoned gardens in nine periods according to the proposed method, a final index for separating active and abandoned gardens is presented based on the NDVI behavior pattern throughout the year. The time series of NDVI for each year is evaluated in 9 periods and garden maps are extracted in 1988 and 2018 to achieve the second objective and prepare the maps of 30-year changes in active gardens in the study area. The rate of change of area and the percentage of changes in the class of gardens are obtained by comparing the maps.Results & DiscussionSince this study is conducted mainly to identify abandoned gardens in urban space, two criteria for assessing user accuracy and errors of commission in the abandoned garden class are very important. In other words, in this problem, the number of gardens that are properly divided into the abandoned garden class is important, and the proposed method provides an accuracy of 86%. The most important issue is the number of abandoned gardens that the proposed method has mistakenly labeled as active gardens, which is 14% in this method. Both accuracies provided are evaluated as acceptable. The overall accuracy of the proposed method is estimated at 82%, which is acceptable, indicating the efficiency of the proposed method.ConclusionOne of the problems facing human societies today is the reduction of forests and gardens. Given the important role that trees play in improving the quality of human life, protecting them is one of the inherent duties of rulers. Various factors cause the destruction of trees, one of which is the development of urban areas in the vicinity of forests and gardens. Traditional methods of conserving natural resources and monitoring their changes have failed in practice. For example, in the study area, 49% of the tree-covered areas have declined over the past 30 years. However, the ban on construction in the area has always been emphasized by city managers in the years under study, and the inefficiency of the methods used has been proven by the statistics provided. New methods of monitoring changes based on satellite image processing can be alternatives to traditional methods due to their high accuracy and speed and significant cost reduction. The proposed index is recommended to be evaluated to separate active and abandoned gardens in other areas facing this problem using images with higher spatial resolution. In different cases of threshold limit, the overall accuracy of the proposed method is examined based on the ground truth data of the evaluator. At best, the separation of active and abandoned gardens is associated with an overall accuracy of 82%.
Remote Sensing (RS)
Mahdiyeh Fathi; Reza Shah-Hosseini
Abstract
Extended AbstractIntroductionRice is an important crop and the main food of more than half of the world’s population, which needs water and heat to grow. Thus, mapping and monitoring rice fields with efficient means such as remote sensing technology is necessary for food security and the lack of ...
Read More
Extended AbstractIntroductionRice is an important crop and the main food of more than half of the world’s population, which needs water and heat to grow. Thus, mapping and monitoring rice fields with efficient means such as remote sensing technology is necessary for food security and the lack of water sources. The phenology extracted from the time series of vegetation indices is used for monitoring and mapping the area under rice cultivation. In addition to the phenological curve, the LST time series map, which is calculated from Landsat 8 images and is related to the phenomenon of evaporation and transpiration of irrigated crops, can cause the separation of rice cultivation from rainfed crops, summer crops, water, etc. Therefore, in this study, the effect of the LST time series map is investigated map for improving the accuracy of rice field identification.Materials & MethodsSince the planting to harvest period of rice is from May to October, in this study, the time-series maps of LST and NDVI for the 3rd of April, 21st of May, 6th of John, 22nd of John, 8th of July, 24th of July, 9th of August, 12th of October, and 28th of October have been calculated after download the Landsat-8 time-series in 2020 The ground truth map of the study area has been obtained from the US Department of Agriculture. To identify rice fields and calculate the LST and NDVI using the Landsat-8 images, initial pre-processing including radiometric and geometric corrections has been applied to these images first. After initial corrections and the calculation of NDVI and LST maps, to identify rice fields in the study area, machine learning algorithms such as Support Vector Machine, K-Nearest Neighborhood, Multilayer Perspective, Logistic Regression, and Decision Tree, have been proposed. Results & DiscussionThe results of the proposed method at the state of California showed that using the time series map of Land surface temperature (LST) with the time-series map of Normalized Difference vegetation Index, improved the results of identifying rice fields (the average Overall Accuracy= + 3/572% and the average kappa coefficient= +7/112%). Visual results showed that some cultivation such as tomato, corn, cucumber, fallow, and water were removed from the rice final map when using the LST time-series map with the NDVI time-series map. According to the numerical results, the Support Vector Machine algorithm (Overall Accuracy 94/28 and Kappa Coefficient 88/29), the Multilayer Perceptron algorithm (Overall Accuracy 94/26 and Kappa Coefficient 88/21), and the K-Nearest Neighborhood algorithm (Overall Accuracy 93/71 and Kappa Coefficient 87/08) showed the highest Overall Accuracy and Kappa Coefficient compared to the Logistic Regression algorithm (overall accuracy 91/96 and kappa coefficient 83/54) and the Decision Tree algorithm (Overall Accuracy 91/34 and Kappa Coefficient 81/97), respectively.ConclusionAlthough, many methods have been proposed to identify rice fields from satellite images. But, the similarity of rice class with other classes is one of the main challenges related to rice identification. In this research, the effect of LST time series maps to improve the identification accuracy of rice fields in Landsat-8 time-series images was investigated. In this study, the effect of the time series map of land surface temperature index extracted from Landsat-8 images on improving the accuracy of identifying rice fields from other rice fields due to the evapotranspiration process using machine learning algorithms was investigated. The results showed the effectiveness of the proposed index in improving the identification accuracy of rice fields. One of the reasons for improving the accuracy of identifying rice fields is to extract the characteristics of the thermal growing season from the Earth's surface temperature time series (LST) maps along with the rice phenology curve. The results showed that due to the flooding of rice fields when using the NDVI time series map, water class and fields summer crops were identified as rice class. But, water and summer crops classes were removed from the rice final map using a land surface temperature time-series map with the extraction of thermal growth season characteristics. Therefore, the results showed that there was a direct relationship between LST time-series maps and rice cultivation.
Mohsen Abedi; Mohammad SaadatSeresht; Reza Shahhoseini
Abstract
Extended Abstract
Introduction
Nowadays, updating information collected from urban areas is of great importance, since it provides the basis for many fields of study such as land cover changes and environmental studies. Remote sensing provides an opportunity to obtain information from urban areas ...
Read More
Extended Abstract
Introduction
Nowadays, updating information collected from urban areas is of great importance, since it provides the basis for many fields of study such as land cover changes and environmental studies. Remote sensing provides an opportunity to obtain information from urban areas at different levels of accuracy while widely used in various change detection applications. Detecting changes in buildings as one of the most important features in urban areas is of particular importance. Powerful and expensive processing systems are the only way to process large volume of remote sensing and photogrammetry data generated by the ever increasing number of sources to which laymen do not have access. The present study has applied deep learning methods and high computational volume of data processing in free clouds to make this possible for the public.
Materials & Methods
Two case studies have been selected in the present study. The first includes DSM and Orthophoto images captured by drones from Mashhad in 2011 and 2016. DSM and Orthophoto images in the second case study has been collected by drones from Aqda in Yazd province in 2015 and 2018. In accordance with the type of data used and high computational volume used for processing, the present study has applied fuzzy clustering method to detect buildings with a high computational speed and deep learning method to detect their changes. Object-based method and fuzzy logic theory have been used in the first step to classify features and detect buildings. In the second step, deep learning method and DSM differentiation method were also used to detect changes in buildings and evaluate results obtained from deep learning method. In the first step, buildings were detected using descriptors extracted from terrestrial and non-terrestrial features, and related decisions were made using fuzzy logic. In the second step, DSM differentiation method has applied the masks extracted from buildings in both epochs on the related DSMs to find their difference and detects changes using an elevation threshold. In deep learning method, a convolutional neural network model was trained to detect changes in buildings during both epochs. Using the DSM of buildings in both epochs and a part of their interface, the network input layers were generated for training. Changes detected in the buildings by the differentiation method were also introduced as the output layer. Following the training and introducing the entire interface in both epochs as the input layer, the trained neural network has detect changes in the buildings. The same process was performed once more using the difference between two DSMs. In other words, a single input layer was used in the network and the rest of the process was the same as before. Finally, changes detected by the neural network was compared with changes detected in the DSM differentiation method
Results & Discussion
In the first step, buildings were detected and images were classified in accordance with the fuzzy logic. The overall accuracy of the first epoch classification in Mashhad equaled 94.6% indicating higher acuracy of object-based methods as compared to pixel-based methods. The overall accuracy of first epoch in Aqda equaled 95.5%. Neural network method detected changes in buildings with an overall accuracy of 90%. In accordance with the ground truth used in network training (both using DSMs as the input layer and the difference between the epochs as the input layer), results indicated that deep learning method is highly accurate in one-dimensional convolution mode. Moreover, the second step has applied the difference between DSMs in the two epochs and thus, many areas lacking a change in height were removed in both epochs and the network was trained more appropriately and accurately.
Conclusion
Necessity of extracting features, especially urban features such as buildings and identifying their changes over time have been investigated in the present study. Due to the high computational volume of modern remote sensing and photogrammetry data and highly expensive systems required for their processing, a new method was presented in the present study to solve this problem. Considering the type of data used and the complexity of features, object-based methods were selected instead of pixel-based methods to identify features and buildings. Deep learning method was used to detect changes in buildings. The method was also compared with DSM differentiation method. A one-dimensional convolutional neural network was used in the deep learning method. Two different modes were used in the network to train and predict changes. In the first, DSMs extracted from the buildings in each epoch were used as the input layer, while in the second one, the difference between DSMs were introduced as a single input layer to the network and the network was trained in accordance with the ground truth collected from areas with and without change obtained from the DSM differentiation method. Following the training process, changes were predicted using the trained network. Much better results were obtained from the second mode in which the difference between DSMs were used.
Sara Khanbani; Reza Shahhoseini
Abstract
Extended AbstractIntroductionChange detection (CD) from remote sensing image considered an important topic among scientist because of its application in monitoring urban and non-urban area, environmental issue, damage assessment, etc. Presenting an efficient CD method from a high-resolution image can ...
Read More
Extended AbstractIntroductionChange detection (CD) from remote sensing image considered an important topic among scientist because of its application in monitoring urban and non-urban area, environmental issue, damage assessment, etc. Presenting an efficient CD method from a high-resolution image can face a different challenge; most of the CD method from a high-resolution image requires training procedure to overcome this challenge. In this paper, an unsupervised (without needing training process) CD algorithm proposed from the high-resolution image. In this method spatial and spectral features extracted from bi-temporal images of the studied area. Difference images generated from high information content features. Then generated different images mapped into spherical space. The Primary change map created using implemented multi-thresholding method on created spherical space and the second change map created using hierarchical clustering regularized by Markov random field method. The final change map created by integrating the result of primary and secondary change maps. The final change map shows an overall accuracy of 92.56% in the studied area. Data and methodsThe data used in this paper is a subset of the main data with dimensions of 2000 * 2000 from an urban area in the city of Mashhad. These images corresponded to the two periods of 1390 and 1395 and were taken with UAV. The orthoimage is related to the first time with a spatial resolution of 6 cm and the second image is taken with a pixel size of 10 cm. In this paper, in order to detect of change of high-resolution images, first, the input images are registered in terms of spectral and spatial, and then feature images are extracted from each input image separately. In the next step, the differences images corresponding to high information content feature images are calculated. . The optimal difference images are mapped to the spherical space using selected statistical methods and in order to better analysis of the results. Otsu multi-thresholding method implemented on r component of sphere space. In the next step, the optimal difference image mapped to a spherical space is divided into non-overlapping blocks with the same dimensions; a cumulative hierarchical clustering method is applied for each block separately. In this case, the computational volume and space proposed in the hierarchical clustering method are reduced. The results of the cumulative clustering of the blocks are merged together and then the Markov random field method is used in order to regularize the results of the cluster in order to reduce noise.In final clustering, the class values below the lowest Otsu threshold are known as unchanged pixels with high reliability and the values above the maximum threshold are determined as changed pixels. The class of middle interval is unknown. For determining, the class of middle interval the corresponded output of hierarchy clustering regularized with a random Markov field is used. In the last step, a vegetation and shadow mask is used for final post-processing. Results and discussionIn order to an accurate assessment of the proposed method on the mentioned study area, a ground truth image with 11073 pixels has been used as a ground test image. The proposed method has shown an overall accuracy of 92.56 in the study area. The accuracy of detecting changed pixels shows 81.61% and the accuracy of detection unchanged pixels shows 92.77%. The false alarm percentage is 0.21 percent and the missed alarm accuracy is 0.0723 percent. For comparative evaluation, the proposed method is compared with the change vector analysis algorithm. In this section, the selected features in the feature extraction section are entered in the change analysis algorithm, and then the multi thresholding algorithm and shadow analysis used to create the final change map. This method has shown increasing the alarm in comparison with the proposed method. The accuracy of changed and un-changed pixels in the change vector analysis method is equal to 52.98 and 89.24%, respectively. Comparing these results with the results of the proposed method shows the efficiency of the proposed method. ConclusionIn this paper, the new unsupervised change detection method presented based on the combination of multi thresholding and the hierarchical clustering algorithm. Compared to supervised methods that require training data, this method does not require training data. In this method, textural and spatial-spectral features are extracted from images with high spatial resolution, which covers the discussion of the importance of neighborhoods in images with high spatial resolution. In the next step, the extracted features that have a high information content are selected, which helps to reduce the redundancy of the information. The contrast images of the features with high information content are created to differentiate the location of the changes. Spherical computing space is considered as the basic computing space. In order to create a binary change map, two analyzes have been performed on the spherical computational space. First, the Otsu multi-thresholding method has been applied. The values of the smaller and larger thresholds have definite classes. But the value of the middle interval needs to be further analyzed using the hierarchical clustering method. In this section, the middle pixel class is examined, and then a final adjustment is performed using Markov field and shadow and vegetation analysis in order to post-process and prevent false changes. In this paper, the parameters of changed accuracy – unchanged accuracy - overall accuracy - false and missed alarms have been used to evaluate the accuracy of the proposed method with a ground accuracy map. In order to make a comparative study, the proposed method is compared with the change vector analysis method of the created feature space. The results show the efficiency of the proposed method.
Arastou Zarei; Reza Shahhoseini; Ronak Ghanbari
Abstract
Extended Abstract
Introduction
As a key parameter describing physics of land surface processes on local and global scales, land Surface Temperature (LST) is the result of all interactions and energy flows between land surface and the atmosphere. Temperature changes rapidly on temporal ...
Read More
Extended Abstract
Introduction
As a key parameter describing physics of land surface processes on local and global scales, land Surface Temperature (LST) is the result of all interactions and energy flows between land surface and the atmosphere. Temperature changes rapidly on temporal and spatial scales, and thus a complete description of LST require measurements involving spatial and temporal frequencies. Hence, climatological, meteorological, and hydrogeological studies require having access to wide scale information about spatial changes of air temperature. Since the LST product of SLSTR uses linear split-window algorithm, the present study has used nonlinear split-window algorithm to estimate LST in Sentinel-3 images. Linearity of the radiation transfer equation in linear algorithm and some approximations used in split-window algorithms (such as transfer approximation as a linear function of vapor value) result in considerable errors because of which nonlinear algorithm is used in the present study. Using linear split-window algorithm to estimate LST in tropical climates also leads to a high level of error. The present study seeks to estimate LST using a nonlinear split-window algorithm and data retrieved from Sentinel-3 in different seasons of 2018 and 2019. The results are also evaluated using temperature product of MODIS and SLSTR.
Materials & Method
A time series of sentinel-3 images retrieved from 2018 to 2019 was used as research data. Data were collected by Sentinel-3 SLSTR sensors operated by the European Space Agency (ESA). Obviously, images shall be radio-metrically corrected before calculating physical land surface parameters such as temperature, emissivity, reflectance and radiance, albedo, and etc. To reach this goal, it is necessary to omit or minimize the effect of atmosphere, epipolar geometry of sensor, sunlight, topography, and surface characteristics while estimating surface parameters in these images. The current study seeks to estimate LST applying a nonlinear split-window algorithm on Sentinel-3 data collected during different seasons of 2018 and 2019 and to evaluate the results using temperature product of MODIS, SLSTR, and in-situ data. Pearson Correlation Coefficient and Root Mean Square Error (RMSE) were also used as relative and quantitative criteria to evaluate the accuracy of the proposed method and determine the deference between temperature calculated by the proposed method and temperature product of MODIS and SLSTR sensor. Hence, four frames of LST product collected by MODIS, and SLSTR in April, June, and October, 2018 and January, 2019 were used to evaluate the proposed method.
Results & Discussion
The proposed method was also indirectly evaluated using temperature products of MODIS and SLSTR sensor. Applying parameters of mean and root mean square error, the evaluation has shown that the results obtained from the proposed method in the one-year reference period were more similar to the results obtained from MODIS sensor. Comparing nonlinear Split-Window algorithm and MODIS products, RMSE ranged from 1.21 to 2.46 and the highest and lowest accuracy belonged to winter and summer, respectively. Comparing this algorithm with the SLSTR product, RMSE ranged from 0.76 to 2.24 and the highest and lowest accuracy belonged to winter and summer, respectively. Proper performance of the algorithm in winter is due to the relative balance of atmospheric water vapour in this season. Comparing nonlinear modelling of atmospheric water vapour in the non-linear algorithm of a Split-window and the linear algorithm in SLSTR and MODIS products, the small difference between temperature calculated by the algorithm and the products can be justified. However, due to temperature fluctuations in summer, results obtained by the proposed method were not reliable enough compared to both temperature products. Generally, results obtained from the proposed method showed a higher correlation with the temperature product of SLSTR sensor, which is due to the similar spectral bands used in calculating the surface temperature. Relative comparison of the Split-Window and the MODIS product’s nonlinear algorithm showed a coefficient of determination ranging from 0.76 to 0.96, while comparing this algorithm with the SLSTR product showed a determination coefficient of 0.80 to 0.98. Comparing temperature obtained from the nonlinear Split-Window algorithm with SLSTR and MODIS temperature products, the proposed algorithm was relatively stable no matter which season was taken into account.
Conclusion
The present study seeks to estimate Land Surface Temperature using a nonlinear Split-Window algorithm and Sentinel-3 data collected in different seasons. Values obtained from the algorithm were validated using in-situ dataset retrieved from the meteorological station. They were also evaluated using temperature product of MODIS and SLSTR. To increase the accuracy level, temperature product of MODIS and SLSTR were also evaluated and compared with the in-situ dataset and provided good results. Generally, there is a significant difference between temperature values estimated by the NSW algorithm for different seasons especially summer. However, a similar trend was observed in temperature changes reported by SLSTR and MODIS, and the proposed algorithm in different seasons of the study area. Although, the nonlinear Split-Window algorithm showed a higher accuracy in spring and winter, overall results indicated that the proposed method was relatively stable no matter which season was taken into account. It can be concluded that LST estimation with nonlinear Split-window method and Sentinel-3 satellite data has an acceptable level of accuracy and thus, can be used in large scale environmental crises such as climate changes.
Saeed Farzaneh; Reza Shahhoseini; Iman Kordpour
Abstract
Introduction Drought is considered to be one of the most widespread natural disasters, ranking second in terms of damages. Due to the complex relationship between hydrological cycle parameters and atmospheric observations, predicting or modeling drought lacks the necessary precision. One of the most ...
Read More
Introduction Drought is considered to be one of the most widespread natural disasters, ranking second in terms of damages. Due to the complex relationship between hydrological cycle parameters and atmospheric observations, predicting or modeling drought lacks the necessary precision. One of the most significant problems in drought monitoring is lack of proper spatial coverage for the collected data (due to unavailibility of field data in some regions) and also lack of a suitable time scale (observations and thus drought estimation is not always possible). Since satellite observations do not face challenges like lack of spatial scale which is quite common in field observations, remote sensing satellites can provide a better estimate of droughts. However, satellite observations alone are not capable of accurately estimating the occurrence of droughts. Therefore, a combination of field and satellite observations has been used recentely to reach a better estimate of hydrological problems. Materials & Methods Temporal and spatial complexity of droughts have made a new global index combining ground-based and satellite-based observations quite necessary. Given the kind of data used in MDI index, we cannot expect it to be global. However, its performance is still acceptable in similar environments and climates, and thus it has been used in the United States (Texas). Datasets selected for the present study have different temporal and spatial scales and thus, a common scale must be found before calculating the index. Data received from GRACE satellite and MODIS sensor were downloaded monthly, but precipitation data were collected on a daily basis. Thus, aritmatic mean of precipitation data was calculated to reach a monthly avarage. Regarding the spatial scale, one-degree precipitation data were received from GRACE and MODIS while precipitation data extracted from synoptic stations had a point-based nature. Therefore, Inverse Distance Weighting (IDW) method was used to produce a one-degree network. Three types of observations were used in the present study including data received from synoptic stations of Iran meteorological organization, GRACE mission satellite-based gravity data and MODIS remote sensing satellite-based data. These were selected to identify droughts over a 14-year time series. Results & Discussion The present study has calculated MDI drought index on a one-degree spatial scale and monthly temporal scale for 168 months using Precipitation, NDVI, and TWS data. Severe droughts in northwestern and central areas of Iran from 2004 to 2014 have led to a shortage of water in reservoirs. In addition to drought, too much water harvesting in northwestern Iran has resulted in a decrease in groundwater level and thus, increased water harvesting from rivers and canals leading to the Urmia Lake and reduced water level in this lake. The results of MDI drought index calculated for Iran over the period of 2000 to 2014 show a high correlation with the results of standardized precipitation-evapotranspiration drought index. According to the type of data used to calculate MDI index, it is expected to have a strong correlation with PDSI index due to its sensitivity to precipitation, area temperature and soil moisture content. Since GRACE and MODIS satellite-based data, and data received from synoptic stations were used, a strong correlation with MDI is also expected. It should be noted that PDSI index is higher than MDI index in Iran, although both show the drought trends accurately. For example according to PDSI index, the worst drought of the last two decades in Iran has occurred in 2008, and MDI index shows the same year. Conclusion The present study has introduced a new drought index using a combination of precipitation data, GRACE_TWS and NDVI. These data were selected because of their high sensitivity to drought. GRACE_TWS observations monitor hydrological drought and include surface and subsurface water sources. NDVI observations are mostly used to identify photosynthetic activities of vegetation cover and are therefore very useful for detecting agricultural drought. Precipitation value shows the amount of surface water in the study area. Precipitation can have relatively rapid effects and is therefore useful for monitoring meteorological drought. MDI index has identified several droughts in each region of the country in the period of 2003 to 2016. These identified droughts have generally covered the country over time. However, each drought has had a different impact on ecosystem. In Iran, the most severe droughts have occurred during 2008 to 2009 and 2011 to 2012. Since MDI correlates well with PDSI, both show a drought in these years. In order to develop the proposed algorithm, the effect of different zoning of the study area on MDI index can be studied.
Ramin Mokhtari Dehkordi; Reza Shahhoseini
Abstract
Extended Abstract Introduction Nowadays, studying urban expansion is very importantin developing countries. Rapid growth of cities has devastating environmental impacts,and irreparable economic and social consequences. Moreover, studying urban expansion is of great importance for managers and planners ...
Read More
Extended Abstract Introduction Nowadays, studying urban expansion is very importantin developing countries. Rapid growth of cities has devastating environmental impacts,and irreparable economic and social consequences. Moreover, studying urban expansion is of great importance for managers and planners of a society. Land surface temperature (LST) is one of the important parameters in urban-regional planning.Urban heat, which is usually referred to as urban heat island, can affect human health, theecosystem, surrounding air, air pollution, urban planning, and energy management. The phenomenon of urban heat island (UHI) is closely related toland-use changes in urban areas, especially when natural surfaces turn intoimpermeable urban surfaces, and increases heat flux and reduces latent heat. Materials & Methods In this study, a collection of Landsat-5 multi-temporal satellite images received in 1986, 1989, 1993, 1998, 2001, 2008, and Landsat 8 multi-temporal satellite images received in 2013, 2015 and 2017, was used along with night images of the MODIS sensor recieved in 2001, 2008, 2013, 2015, 2017 (on the same day Landsat-5 and Landsat-8 satellite images were received). In order to classify land cover and calculate land surface temperature usingLandsat 5, Landsat 8 and MODIS sensorsatellite images, initial pre-processing (radiometric and geometric corrections)was performed.In order to classifyland cover in the study area, training areas were selected using Google Earth andthen, land cover classification was carried outusing Neural Network Algorithm. Since, classifying urban areas wasthe priority ofthe present study, Normalized Difference Built-up Index (NDBI) was also used.Ultimately, pixelidentified by classification algorithm and NDBI index was allocated tourban areas. A simple relationship suggested by the United States Geological Survey (USGS) was used to estimate land surface temperature from Landsat-5 imageries.Split-window algorithm was also used to estimate land surface temperature from Landsat-8 and MODIS imageries. Since, Landsat-8 and MODIS imageries were collectedwith only afew hours (or less than that)time difference, and their thermal bands’spectral rangeswere close to each other, Landsat-8 thermal bands’emissivity coefficient with a higher spatial resolution (30 m) was used to calculate land surface temperature from MODIS images. Results & Discussion Classifying land cover in Shahr-e Kordusing Landsat-5 and Landsat-8 imageries received in 1986, 1989, 1993, 1998, 2001, 2008, 2013, 2015, and 2017 indicated that in this31-year time period,residential areas were approximately duplicatedand reached from 1004 hectares to 2112 hectares. Analysis of land surface temperature maps using Landsat 5, and Landsat 8 imageries indicated that urban areas and areas with dense vegetation had lower surface temperatures compared to areas with thin vegetation cover. Therefore, land surface temperature of urban areas is lower than the surrounding areas. However, land surface temperature obtained from MODIS imageries indicated that land surface temperature of urban areas is higher at nights. Therefore, urban heat islands in this city occur at nights. Results indicated that with increasingexpansion of urban areas, urban heat islands also intensifyat nights. Conclusion Although, Shahr-ekordis a less developed urban area as compared to other Iranian metropolises,expansion of its constructed areas can stillhave negative effects on the environment and climate of the region. The present study investigates urban growth, and itsinfluence on land surface temperature and occurrence of urban heat island. Thermal maps produced in the present study indicated that daytime air temperature of this city was relatively lower than other regions. But this is not the case at nights: compared to other areas,residential areas have a higher temperature at nights. This indicates the existence of a heat island in the city, and possibly have adverse and devastating effects on humidity, reduces precipitation, changes local winds and the climate. Results also indicate that urban expansion have directlyaffected urban heat islands. Thus, urban heat islandshave intensified and expanded during this time period. Therefore, it is concluded that there is a direct relationship between land surface temperature and land use type.
Arash Karimi Zarchi; Reza Shahhoseini
Abstract
Extended Abstract Introduction Heat island phenomenon occurs when the land surface temperature and the air temperature in urban areas are higher than that of the surrounding areas. This temperature difference is shown as the urban heat islands on thermal maps. Information obtained from the urban heat ...
Read More
Extended Abstract Introduction Heat island phenomenon occurs when the land surface temperature and the air temperature in urban areas are higher than that of the surrounding areas. This temperature difference is shown as the urban heat islands on thermal maps. Information obtained from the urban heat islands can be a useful source in urban planning applications. The availability of reliable information about the urban heat islands plays an important role in predicting and preventing the occurrence of many heating risks in urban areas. One of the common methods of calculating heat islands intensity in urban areas is the use of two temperature sensors installed in the city and around it. Given the limited temperaturemeasuring stations, there is no accurate estimate of the urban heat islands. With the introduction of Remote Sensing technology into the space arena, and with the help of satellite images processing, a precise map can be produced for the land surface temperature, i.e. a precise estimation of the urban heat islands is obtained by calculating the pixels temperature difference at the urban areas and around them. Therefore, one of the important issues in such studies is to detect the urban and non-urban pixels and to separate them from each other. Materials&Methods The most important reason for the occurrence of the heat island phenomenon is the change in land use from rural to urban, which is well exhibited in the urban cover index maps.In this paper, in order to measurethe intensity of surface urban heat islands, a method based on generating the urban percentage map was proposed by combining the Land Surface Temperature (LST) map, the Normalized Difference Built-up Index (NDBI) map and the Normalized Difference Vegetation Index (NDVI) map.Considering the relationship between the land surface temperature and the land cover type, it can be said that the relationship between the land surface temperature and the urban percentage map follows a linear function which can be fitted to the land surface temperature graph in terms of land cover type. Finally, the Urban Heat Island Intensity (UHII) map was calculatedfrom the slope of the fitted line.In order to evaluate the strengths and weaknesses of the proposed method, a classification-based method was used to separate the urban and non-urban pixels and to calculate the urban heat island intensity. The proposed method was implemented on the Landsat-7 ETM + satellite data in the city of Rasht and on the Landsat-8 OLI / TIR satellite data in the city of Langroud. Results&Discussion The results of the classification-based method indicated a large difference between the maximum and the minimum temperature of the urban areas, which led to a high-temperature changein all land cover typesin the study area. Therefore, the use of the average temperature of each class to calculate the heat island intensity is not a suitable method and the accuracy of the heat islands maps is not high and they cannot be used in applications that require high precision.Although, this problem can be solved by increasing the number of classes, increasing the number of classes requires more training data and a sensor with higher spatial resolution. By contrast, the results indicated that the proposed method (based on the urban percentage map) had a high accuracy for calculating the urban heat island intensity which was similar for both study areas. Also, fitting a linear function to the values of land surface temperature and the urban percentage map led to decreasing the effect of suspicious pixels (noisy pixels) on the overall accuracy of the estimation of the urban heat island intensity. Meanwhile, the results obtained on two datasets indicated that this method did not require any training data or any other background information about the study area and it can be applied for many satellite images having thermal band with any spatial resolution. However, because of the ineffectiveness of urban cover indicators in desert areas, the heat islands intensity in these regions was underestimated. Conclusion In applications that do not require high accuracy in calculating the urban heat island intensity, and there are high spatial resolution satellite imagery and sufficient training data in a region, the use of a classification-based approach seems to be suitable. Since the collection of such data and information is costly, a new method based on the urban percentage map was proposed in this paper by fitting a line to the LST parameter diagram in terms of the NDBI index for measuring the heat island intensity. The results indicated the higher efficiency and accuracy of the proposed method compared to the conventional classification-based methods for calculating the urban heat island intensity.