چکیده:
خسارات ناشی از سیلاب در سرتاسر جهان باعث شده است تا از این رخداد طبیعی به عنوان یکی از مهمترین بلایای طبیعی نیز یاد شود. امروزه همراه با پیشرفت فناوری، انسان قادر است مناطق مستعد سیلاب را تعیین و از جنبههای گوناگون مورد بررسی قرار دهد. در این تحقیق با بهرهگیری از تصاویر راداری ماهواره سنتینل-1 و روش آستانهگذاری Otsu، مناطق سیلزده و مستعد سیلاب با توجه به سیلابهای فروردین ماه 1398 در محدوده رودخانه کشکان تعیین گردید. حدود آستانه Otsu در این حقیق، در نرم افزار متلب محاسبه شدند. همچنین مراحل پیشپردازش و پردازش تصاویر راداری در نرم افزار SNAP صورت گرفت. براساس نتایج، میانگین بارش روزانه معادل 46/47 میلیمتر باعث شده تا 51/6 درصد از منطقه مورد مطالعه در کلاس مناطق سیلزده قرار بگیرد، این در حالی است که با حدود 31 میلیمتر بارندگی، 69/3 درصد از منطقه در کلاس سیلزده قرار گرفته است. همچنین نتایج نشان داد که محل پیچانرودها از مهمترین مناطق مستعد سیلاب در مسیر رودخانه کشکان به شمار میروند. حفظ حریم رودخانه و عدم ساختوساز و فعالیت انسانی در مناطق مئاندری میتواند حجم خسارات ناشی از سیلاب را کاهش دهد. یافتههای تحقیق نشان داد که تصاویر راداری ابزاری مفید و قدرتمند در تعیین مناطق سیلزده خصوصاً برای بررسیهای سریع و نزدیک به زمان وقوع پدیده (Near Real Time) هستند.
Determination of flood-prone areas using Sentinel-1 Radar images
(Case study: Flood on March 2019, Kashkan River, Lorestan Province)
Introduction
Although natural hazards occur in all parts of the world, their incidence is higher in Asia than in any other part of the world. Natural phenomena are considered as natural hazards when they cause damage or financial losses to human beings. Iran is also one of the high-risk countries in terms of floods. Until 2002, about 467 floods have been recorded by the country's hydrometric stations. In addition to natural factors such as rainfall, researchers consider human impacts such as destruction of vegetation cover, soil destruction, inefficient management, destruction of pastures and forests, and encroachment on the river are the most important factors for the occurrence and damage of floods in the country. One of the most efficient and emerging tools in flood surveys is the use of radar images. SAR images and flood maps produced by radar images provide researchers valuable and reliable information. Moreover, maps obtained from SAR images help officials to manage the crisis and take preventive measures against floods. The Sentinel-1 satellite is part of the Copernicus program, launched by the European Space Agency, and is widely used in mapping flood-prone areas. The contribution of Sentinel-1 to the application of flood mapping arises from the sensitivity of the backscatter signal to open water. This study aims to determine high-risk and flood-prone areas along the Kashkan River using Sentinel-1 radar images.
Data and Methods
The study area includes a part of the Kashkan river from Mamolan city to the connection point of this river to Seymareh river, after Pol-dokhtar city. The average annual discharge of the Kashkan river is 33.2 cubic meters per second based on the data of the Pole-Kashkan Station. The length of the river in the study area is about 100 km. To investigate flood-prone areas, we applied pre-processing and image-processing steps to each flood event including SAR images belonging to March 25th, 2019, March 31st 2019, and April 2nd, 2019. SAR images were acquired from ESA Copernicus Open Access Hub. climatic data was downloaded from power.larc.nasa.gov. To create meander cross-sections, the Digital Elevation Model of the studied area was utilized. Cross-sections were created using QGIS software. Pre-processing steps include: applying orbit data, removing SAR thermal noise, calibration of SAR images, de-speckling and topographic correction. In image processing, we applied the Otsu thresholding method to distinguish water pixels from land pixels. In thresholding methods, the histogram of each image is divided into two parts according to the amount of gray composition. The higher the amount of gray (i.e., the pixel tends to be darker), the more pixels represent water, and conversely, the lighter-toned pixels (i.e., pixels that tend to whiten) represent land. The Otsu thresholding method is a commonly used method for water detection in SAR images. It uses an image histogram to determine the correct threshold. The most important feature of the Otsu method is that it is capable of determining the threshold automatically. The Otsu algorithm was applied to all images using MATLAB.
Results
According to the flood maps, on March 25th, 6.51 percent of the study area was flooded, while on March 31th, only 3.96 percent was flooded. This is mainly due to less precipitation on the 31st. On March 25th the average daily precipitation was 47.46 mm while on 31st of March the average daily precipitation was 31.64 mm. On April 2nd, however, there was no rainfall, on the day before more than 63 mm of precipitation has occurred. This massive amount of precipitation on the previous day has led to more than 25km2 being flooded in the studied area.
Conclusion
Results showed that meanders and their surrounding areas are the most dangerous sections in terms of flooding. The meander's dynamic and the river's hydrologic processes are essential factors affecting flooding in those sections. Generally, various factors affect flooding and the damage caused by it. This study aimed to determine flooded and flood-prone areas (according to flooded areas in previous events) using new methods in a short time and with high accuracy to use this tool for more accurate zoning and efficient planning in the future. The results showed that radar images are practical, robust, and reliable tools for determining flooded areas, especially for rapid and near-real-time studies of flood events.
Keywords: Floods, Radar images, Sentinel-1Satelitte, Kashkan river
خلاصه ماشینی:
در اين تحقيق با بهره گيري از تصاوير راداري ماهواره سنتينل -١ و روش آستانه گذاري Otsu، مناطق سيل زده و مستعد سيلاب با توجه به سيلاب هاي فروردين ماه ١٣٩٨ در محدوده رودخانه کشکان تعيين گرديد.
يافته هاي تحقيق نشان داد که تصاوير راداري ابزاري مفيد و قدرتمند در تعيين مناطق سيل زده خصوصاً براي بررسيهاي سريع و نزديک به زمان وقوع پديده (Near Real Time) هستند.
در اين مطالعه ، براي تعيين مناطق سيل زده ، از تصاوير راداري ، GRDH٣ ماهواره سنتينل -١، با قدرت تفکيک مکاني ١٠ متر، استفاده شده است .
در اين تحقيق از نرم افزار SNAP به منظور اعمال مراحل پردازشي بر روي تصاوير راداري، نرم افزار MATLAB براي اعمال آستانه گذاري روش Otsu و نرم افزارهاي سيستم اطلاعات جغرافيايي از 1 European Space Agency 2 Copernicus 3 Ground Range Detected High-resolution جمله ARCMAP و QGIS براي تهيه نقشه مناطق سيل زده و مستعد سيلاب استفاده شد.
روش هاي آستانه گذاري در تعيين مناطق سيل زده توسط تصاوير SAR به طور گسترده توسط محققين مورد استفاده قرار گرفته است (٢٠١٢، Fan and Lei؛ ٢٠١٤، .
حدود آستانه تعيين شده ، مساحت و درصد مساحت مناطق سيل زده و ميزان بارش در هر زور (رجوع شود به تصویر صفحه) • پنجم فروردين ماه بارندگي در روز پنجم فروردين ١٣٩٨ منجر به وقوع سيلاب در محدوده مورد مطالعه شد.
در اين تحقيق نيز با استفاده از تصاوير GRDH ماهواره سنتينل -١ و با استفاده از روش آستانه گذاري Otsu، مناطق سيل زده در حوادث سيلاب فروردين ماه ١٣٩٨ در مسير رودخانه کشکان ، مورد بررسي قرار گرفت .