محاسبه مساحت مناطق سیل زده با تصاویر سنتینل 1 و 3 – در سامانه Google Earth Engine
۳۶۹,۰۰۰ تومان
- مدرس: امیرحسین احراری
- تخصص: کارشناس سنجش از دور و پردازش تصاویر ماهواره ای
- نوع محصول: دانلودی
- داده های تمرینی: دارد.
- پیش نیاز: آشنایی مقدماتی با گوگل ارث انجین
- نرم افزار: سامانه Google Earth Engine
- موضوع: کاربرد سنجش از دور در مخاطرات طبیعی
توضیحات
محصول آموزشی برآورد مساحت مناطق سیل زده با استفاده از تصاویر سنتینل ۱ و ۳ در سامانه Google Earth Engine منتشر شد. در این محصول چگونگی استفاده از تصاویر ماهواره ای سنتینل ۱ و ۳ در سامانه گوگل ارث انجین در راستای شناسایی مناطق سیل زده آموزش داده شده است. در این سامانه بدون نیاز به دانلود داده ها و با سرعتی بسیار بالا امکان پردازش تمامی مراحل این پروژه به سادگی فراهم شده است.
توضیحات محصول
در این محصول آموزشی، امیرحسین احراری روش محاسبه مساحت مناطق سیل زده را ابتدا با استفاده از تصاویر راداری و سپس با استفاده از تصاویر اپتیکی در گوگل ارث انجین، آموزش داده است. در همین راستا از تصاویر راداری ماهواره سنتینل ۱ و اپتیکی سنتینل ۳ در سامانه پردازشی گوگل ارث انجین استفاده شده است. استفاده از تصاویر ماهواره ای متنوع برای انجام محاسبات متفاوت، افزایش دقت و کاهش عدم قطعیت در نتایج سنجش از دور را به دنبال دارد.
عناوین آموزشی
نکات آموزش داده شده در سامانه گوگل ارث انجین و ماهواره های سنتینل ۱ و ۳ عبارت اند از:
- روش فراخوانی تصاویر
- روش انتخاب منطقه مورد مطالعه
- روش برش زدن تصاویر ورودی
- روش فراخوانی تصاویر سری زمانی
- روش اعمال فرمول در تصاویر سری زمانی
- روش محاسبه تصویر کمینه
- روش محاسبه تصویر میانه
- روش شناسایی پهنه های آبی
- روش شناسایی محدوده های سیل زده
- روش خروجی گرفتن از نتایج
- روش محاسبه مساحت مناطق سیل زده
ماهواره Sentinel-1
سنتینل ۱، ماهواره ای کاربردی است که قابلیت تصویربرداری در محدوده مایکرویو فعال (رادار) را دارا است. این ماهواره با توان تفکیک مکانی ۱۰ متر و توان تفکیک زمانی ۵ تا ۶ روز از نقش موثری در مطالعه پدیده های مختلف سطح زمین برخوردار است. ماهواره سنتینل ۱ در محدوده باند c تصویر برداری کرده و برای مطالعات کاربری اراضی، تولید مدل های رقومی زمین، کشاورزی، رطوبت خاک و محاطرات زیست محیطی بسیار کاربردی است.
ماهواره Sentinel-3
یکی از سنجنده های چندگانه ماهواره سنتینل ۳، OLCI نام دارد. این سنجنده از ۲۱ باند در محدوده مرئی و مادون قرمز نزدیک تشکیل شده که برای مطالعات کشاورزی، پوشش گیاهی و آب بسیار مطلوب و کاربردی است. توان تفکیک زمانی ۲ روزه به همراه توان تفکیک مکانی ۳۰۰ متری آن امکان مطالعه و پایش تغییرات پهنه های آبی و گیاهی سطح زمین را با دقت بالایی فراهم ساخته است.
کاربرد سنجش از دور
سنجش از دور ابزاری مهم در مطالعات مخاطرات محیطی است. با استفاده از سنجش از دور وقوع برخی از مخاطرات مانند طوفان ها قابل پیش بینی است. از سوی دیگر سنجش از دور امکان پایش تغییرات ناشی از مخاطرات را فراهم می سازد. در همین راستا تصاویر ماهواره ای می توانند برای برآورد خسارت های ناشی ار مخاطرات طبیعی و انسانی نیز مورد استفاده قرار بگیرند.
نقش سنجش از دور در مخاطره سیل نیز حائز اهمیت است. به عنوان مثال خرابی ها و خسارت های ناشی از وقوع سیل به خوبی در تصاویر ماهواره ای قابل ردیابی است. درهمین راستا، شما می توانید مساحت محدوده سیل زده را با استفاده از تصاویر ماهواره ای محاسبه کنید. به علاوه، دقت بالای زمانی و مکانی تصاویر ماهواره ای نقش مهمی در افزایش دقت برآورد مساحت محدوده سیل زده خواهد داشت.
منطقه مورد مطالعه
استان خوزستان به عنوان منطقه مورد مطالعه، در این محصول آموزشی در نظر گرفته شده است. در گام نخست با استفاده از تصاویر چند زمانه راداری ماهواره سنیتنل ۱ محدوده سیل زده در این استان شناسایی و سپس مساحت آن محاسبه شد. در شکل مقابل محدوده سیل زده در تصویر راداری به خوبی قابل تشخیص است.
در گام دوم، تصاویر اپتیکی ماهواره سنتینل ۳ سنجنده OLCI برای شناسایی مناطق سیل زده با استفاده از شاخص های طیفی کاربردی در این زمینه، مورد استفاده قرار گرفت. در شکل روبروز مناطق سیل زده استان خوزستان به رنگ روشن و بصورت واضح قابل تشخیص است.
لینک های مفید
- کارگاه رایگان گوگل ارث انجین
- کاربرد تصاویر سنتینل ۳ در جزیره حرارتی شهری
- آشکارسازی مناطق سیل زده در نرم افزار SNAP
- آشکارسازی مناطق سیل زده در نرم افزار ENVI




مصطفی –
سلام جناب احراری، وقت بخیر. سوال من در مورد این آموزش از شما این است که آیا روشی برای اعتبار سنجی این دو روشی که در این آموزش به آنها پرداخته اید وجود دارد؟؟
امیرحسین احراری –
با سلام و احترام
در این آموزش مساحت سطح مناطق سیل زده قابل استخراج است. با مقایسه عدد آن و مقایسه گزارش سطح مناطقی که تحت تاثیر سیل قرار گرفته اند می توانید اعتبار سنجی کنید.
موفق باشید
امیر –
با سلام
در قسمت آخر که تصویر رو خروجی گرفتید اعداد scale و maxpixels رو بر چه اساسی تعیین کردید؟ کلا میخوام بدونم برای کارهای مختلف این دو عدد رو باید چه جوری تعیین کنم؟ و اینکه من میخوام فقط مناطق سیل زده رو خروجی بگیرم و ببرم تو arcgis ولی نمیدونم چجوری این کارو انجام بدم. ممنون میشم راهنمایی کنید.
امیرحسین احراری –
با سلام و احترام
scale بر اساس توان تفکیک مکانی تعیین میشه. اگر تصویر ۱۰ متری هست باید عدد ۱۰ وارد شود. البته من در آموزش ها برای اینکه خروجی ها رو سریع تر بگیرم عدد بالاتری برای این کار انتخاب میکنم.
maxPixels ماکزیموم پیکسل هایی است که ارث انجین برای هر تصویر خروجی در نظر میگیره. این عدد یک عدد ثابت است 1e13
موفق باشید
شهرام شریفی هشجین –
سلام
آیا تو این محصول بین خروجی راداری و اپتیک در محیط ارث انجین مقایسه یا ارزیابی دقت صورت گرفته؟ مثلا به شرط مهیا بودن شرایط، تصویر سیل مستخرج از سنتینل ۱ با شاخص نرمالیزه آب مستخرج از اپتیک مورد ارزیابی قرار گرفته؟
با تشکر
امیرحسین احراری –
با سلام و احترام
در این محصول آموزشی ارزیابی دقت بدلیل عدم در اختیار داشتن داده های زمینی مناسب انجام نشده است.
موفق باشید
mohamad –
سلام
من یک دیدگاه ثبت کرده بودم ولی متاسفانه تایید نشده
لطف میکنید تایید کنید که آقای مهندس لطف کنند پاسخ بدند
امیرحسین احراری –
با سلام و احترام
همه دیدگاه ها امروز تایید میشه و مدرسین پاسخ میدهند. اگر از بنده پرسیدید امروز همه سوالات رو جواب میدم.
موفق باشید
mohamad –
سلام جناب آقای احراری
یک سئال دارم از خدمت شما، شما این محاسبات را بر این اساس انجام دادید که دی اده هاراداری مناطق شهری و مناطق کوهستانی را روشن نمایش می دهد و آب را تیره.
سوال بنده این جاست اگر بخواهیم flood detection را برای یک منطقه شهری انجام بدیم به هدف این که ببینیم در کدام نواحی شهر سیل بیشتر بوده حال در این صورت ایا میشود از این الگوریتمی که ارائه دادید استفاده کرد زیرا کلا شهر را روشن نشان میدهد.
آیا به طور کلی برای flood detection در مناطق شهری الگوریتم معتبری وجود دارد
ممنئن
امیرحسین احراری –
با سلام و احترام
همین روش را می توانید برای مناطق شهری استفاده کنید. البته من خودم در شهر تا کنون تست نکردم. اما این روش آشکارسازی تغییراتی که آموزش دادم تابعی از منطقه مورد مطالعه نیست و یک تکنیک مرسوم و کلی است و برای همه جای دنیا قابل استفاده است.
به نظرم برای شهر باید قابل استفاده باشه. بصورت پیش فرض مشکلی نباید داشته باشه چون بازپراکنش شهر آب کاملا متمایز و جدا از همدیگر هست.
موفق باشید
mohamad –
بله بنده این الگوریتم را برای یکی از شهرستان ها پیاده کردم ولی متاسفانه اصلا خروجی مطلوبی حاصل نشد
امیرحسین احراری –
البته این روهم عرض کنم که منطقه مورد مطالعه و مقیاس رخداد نیز بسیار مهم است.
این سیلی که بنده انتخاب کردم خیلی وسیع بوده و کارایی روش به خوبی دیده شده. در صورتی که وسعت منطقه و رخداد کوچک باشد قطعا در نتایج مشکل ساز خواهد شد.
موفق باشید
امیر –
با سلام
تخفیف هایی که گذاشته میشه خیلی کمه نهایتا تا سقف ۳۵ درصد است. کلا فکر میکنم فقط یکبار تخفیف ۵۰ درصدی گذاشتید اونم برای افتتاح سایت جدیدتون بود. ممنون میشم تخفیف هارو بیشتر کنید که همه بتونن از محصولات خوبتون استفاده کنند.
با تشکر
حسین زاده –
سلام و وقت شما به خیر،
با تشکر از آموزشهای مفیدتون، بنده این کدهارو برای یک تاریخ واحد (یک روز مشخص) می خواستم اما وقتی به جای مجموعه تصاویر، آی دی مربوط به یک تصویر واحدرو داخل کد می نویسم اجرا نمیشه و خطا میده، فیلتر مربوط به زمان و حتی تک تک فیلترهارو یک یک حذف کردم باز هم اجرا نمیشه، ممنون میشم اگر راجع به فراخوانی تصویر سیل و قبل از سیل (هر کدام تنها یک تصویر با آی دی مشخص) راهنمایی بفرمایید.
همچنین وقتی برای یک روز مشخص مثلا ۱۳ ژانویه، ۳ تصویر راداری موجود هست و هدف ما فراخوانی تنها یک تصویر است از چه فیلترهایی باید استفاده کنیم تا این ۳ تصویر را به یک تصویر تبدیل کند؟
ممنونم از زحمات شما.
FerSa –
سلام آقای احراری
خیلی ممنون از تهیه این بسته آموزشی.
برای سالهای ۲۰۰۰ تا ۲۰۰۵ دنبال داده راداری جهت شناسایی سیل هستم. شما پیشنهاد خاصی دارید؟
امیرحسین احراری –
با سلام و احترام
متاسفانه برای این بازه زمانی اطلاعات دقیقی ندارم.
موفق باشید
محمد –
سلام جناب احراری عزیز با تشکر از شما
چرا امکان دانلود تصاویر سنتینل ۲ از سایت Copernicus Open Access Hub نیست و برروی تصاویر offline میزنه
با تشکر
امیرحسین احراری –
با سلام و احترام
تصاویر نسبتا قدیمی حتی برای یک سال پیش بصورت آفلاین شده است. باید روی گزینه Add to card کلیک کرده و سپس چند ساعت صبر کنید. کمتر از ۲۴ ساعت بعد داخل اکانتتون بشید و سپس وارد بخش card شوید تا امکان دانلود آن فراهم شده باشد.
موفق باشید
hamidreza rahimi –
با سلام و خسته نباشید. ببخشید برای محاسبه مساحت مناطق سیل زده از با snap بهتره محاسبه کنیم یا گوگل ارث انجین؟
امیر –
سلام و با تشکر. جناب احراری بعد از اینکه شما از تصاویر خروجی گرفتید، وقتی که داخل نرم افزار اسنپ باز میکنیم فقط مقادیر ۰ یا ۱ را نشان میدهد و مقادیر اعشاری نشون نمیده. لطفا اون کدی که برای خروجی نوشته میشه و مقادیر اعشاری هم شاملش میشه بفرمایید. ممنون از شما
امیرحسین احراری –
با سلام و احترام
راستش بنده دقیقا از همین کد استفاده کردم. غیر از آنجه که در ویدیو توضیح دادم و کدش هم ضمیمه شده کد دیگری برای این کار ندارم. احتمالا در بخشی از محاسبات اشتباهی داره صورت میگیره. تنها راهکار این هست که با کد بنده مقایسه کنید.
موفق باشید
امیر –
باسلام اشتباهی صورت گرفته فقط چون اعداد بصورت float نیست اعداد اعشاری نشون نمیده. شماهم میتونید اون تصویری که اخذ کردید ببینید که فقط مقادیر ۰ و ۱ نشون میده. برای اخذ تصاویر که اعداد اعشاری هم نشون بده بایستی در قسمت export یک کد دیگه نوشته بشه. با تشکر
امیرحسین احراری –
تصویر خروجی ما یک تصویر ماسک شده ۰ و ۱ هست و قرار نیست بصورت اعشاری باشه. چون حدآستانه رو قبلا تعیین کردیم و خروجی فقط یک تصویر ۰ و ۱ هست.
موفق باشید
محمد –
سلام جناب احراری با تشکر از محصول خوبتون. سوالی که دارم این هست که مثلا اگه سال ۲۰۲۰ سیل کمتر از سال ۲۰۱۹ باشه و اینو به ترتیب تایم ۱ ئ ۲ قرار بدیم و اون ndci اجرا کنیم تصویره درستی به ما نمیده و حتی جابجا کردم. به نظر شما حتما بایستی pre flood سالش عقب تر از postflood باشه؟ ممنون از شما
امیرحسین احراری –
با سلام و احترام
ترتیب باید به این صورت باشد مگر اینکه در یک نمونه عملی خلاف آن تست شود و نتیجه درست بگیرید. می توانید حتی یک سال بدون سیل را با یک سال سیل زده مقایسه کنید. لزوما نباید دقیقا سال های پشت سرهم باشد.
موفق باشید
اسماعیلی –
با سلام
جناب آقای مهندس احراری
برای پروژه ای تصاویر لندست ۵، ۷ و ۸ را از طریق GEE دریافت کرده و بصورت stack آن ها را export کردم. برای کار با این تصاویر در انوی مشکلی دارم به این صورت که مانند تصاویر معمولی با باندهای کامل، از این تصویر استک شده نمی توانم استفاده کنم. لطفا راهنمایی فرمایید.
با تشکر
امیرحسین احراری –
با سلام و احترام
دقیقا چه مشکلی با این داده ها دارید. ممنون میشم دقیق تر توضیح بدید که بتونم بررسی کنم.
موفق باشید
مشفقی –
با عرض سلام خدمت استاد احراری بزرگوار
من برای دانلود تصاویر ماهواره ای جهت سیل در google earth engine با مشکل مواجه شدم.
کد را مطابق آموزش شما وارد کردم اما باز ایراد می گیرد. اگر محبت بفرمایید کنترل و اشکال را اطلاع بدهید ممنون می شوم.
Map.addLayer(table);
Map.centerObject(table);
// Load the Sentinel-1 ImageCollection.
var sentinel1 = ee.ImageCollection(‘COPERNICUS/S1_GRD’)
.filterBounds(table)
.filterDate(‘2019-01-01′,’2019-02-01’)
.filter(ee.Filter.listContains(‘transmitterReceiverPolarisation’,’VV’));
.filter(ee.Filter.eq(‘instrumentMode’,’IW’);
.select(‘VV’)
.min()
.divide(100)
.clip(table)
.rename(‘time1’)
// Display as a preflood.
Map.addLayer(pre);
پیشاپیش عذرخواهی میکنم که گویا نشد تصویر خطا را خدمتتان ارسال کنم.
امیرحسین احراری –
با سلام و احترام
لطفاً با استفاده از گزینه get link آدرس کدتون رو در اینجا برای بنده کامنت بزارید و مشکل را مطرح کنید تا بنده آن را اصلاح کنم.
موفق باشید
سعید کرمانشاهی –
سلام و وقت بخیر. توی آموزشهاتون کار با دیتاستهای اماده گوگل ارث انجین جهت شناسایی زمینهای کشاورزی هم وجود داره؟
امیرحسین احراری –
با سلام و احترام
یک آموزش آقای کاکویی در این رابطه تولید کردند البته علاوه بر پروداکت های آماده از داده های لندست و سنتینل نیز استفاده کردند.
https://girs.ir/cropland-in-gee/
موفق باشید
سیدامیرحسین عقیلی –
سلام روز بخیر
لینک کد :https://code.earthengine.google.com/b2dd3e6670d38f2854332c656c6b0819
من چون شیپ فایل رو نتوستم ایمپورت کنم مجبور شدم rectangle بکشم که آموزش رو جلو برم ولی بازم به مشکل خورد
شیپ فایل ها هم که ایمپورت نمیشدند
امیرحسین احراری –
با سلام و احترام
من کد شما را بررسی کردم. لینک کد اصلاح شده ضمیمه شده است. مشکل کد شما از انتخاب منطقه مورد مطالعه نیست پارامترهای transmitter و instrument را با دیکته اشتباه تایپ کرده بودید.
https://code.earthengine.google.com/1938fe4792a79a4fbdb1018b9aa4b41b
برای لود شیپ فایل ها دقیقا چه خطایی به شما نمایش میده؟ عبارت دقیق خطا را در این جا کامنت بزارید تا بررسی کنم.
سیدامیرحسین عقیلی –
سلام و روز بخیر
خیلی ممنون از وقتی که گذاشتید
پیغام خطا رو در کامنت های بالا گذاشتم ولی لازم نیست چندبار امتحان کردم دوباره و بلاخره فعلا درست شد
امیرحسین احراری –
با سلام و احترام
خواهش میکنم. موفق باشید
سیدامیرحسین عقیلی –
سلام وقت بخیر
من وقتی فایل زیپ شده خوزیتانی که قرار دادید رو ایمپورت می کنم چنین پیغامی میده وبا مشکل مواجه میشه
URI “gs://earthengine-uploads/GJMtdnpjfpXLSJ1N/khozestan.zip” is invalid. All URIs must have a “.shp” extension.
امیرحسین احراری –
با سلام و احترام
مشکل از شیپ فایل ورودی هست. لطفا شیپ فایلتون رو برای بنده ارسال کنید تا بررسی کنم مشکلش رو.
موفق باشید
سیدامیرحسین عقیلی –
سلام ایمیل کردم
رسید بهتون؟
امیرحسین احراری –
سلام بله رسیده. تعداد سوالات زیاد هست هنوز نرسیدم داده های شما رو بررسی کنم. بهتون اطلاع میدم.
موفق باشید
سیدامیرحسین عقیلی –
سلام ایمیل فرستادم
ایمیل بهتون رسید
اگر نرسید لطفا ایمیلتون رو بفرمایید
امیرحسین احراری –
با سلام و احترام
بله تعداد سوالات بسیار زیاد هست. هنوز نرسیدم داده های شما رو بررسی کنم. بهتون اطلاع خواهم داد.
موفق باشید
هادی امامی –
سلام مهندس احراری خسته نباشین. آیا با همین محصول میشه برف رو هم برآورد کرد؟؟
یعنی مشابه هم هستن روش هاش؟
Amirhossein.Ahrari@Gmail.com –
با سلام و احترام
اگر شاخص برف رو جایگزینش کنید بله قابل استفاده هست.
موفق باشید
hamidreza rahimi –
با سلام و خسته نباشید. ببخشید شاخص برف رو چطور جایگزین کنیم؟
امیرحسین احراری –
با سلام و احترام
داخل دستور select اسم باند شاخص برف یعنی NDSI را تایپ کنید.
موفق باشید
آیسان –
نظرتون راجع به موزاییک کردنشون چیه؟ آیا کدی هست که من بتونم اونا رو موزاییک کنم …؟
Amirhossein.Ahrari@Gmail.com –
گوگل ارث انجین خروجی موازییک شده میده اما چون حجمش زیاد بوده به چند فایل اونها رو تفسیم بندی کرده. پیشنهاد میکنم از نرم افزار انوی برای موازییک کردن آنها استفاده کنید.
موفق باشید
آیسان –
با سلام و احترام
بله من به کمک انوی تونستم موزاییک کنم ….متشکرم از راهنمایی تون.
میتونم ازتون یه سوال کنم …شیپ فایلی که شما استفاده کردین برای خوزستان یه مقدار با اون چیزی که من معرفی کردم به GEE فرق میکنه در حاشیه جنوبی…این موضوع خیلی مهم هست یا نه؟
و یه نکته دیگه اینکه من کدی خدمتتون ایمیل کرده بودم، فک کنم اینقد مشغله دارین نرسیدین زحمت بررسی اش رو بکشین.
با تشکر
Amirhossein.Ahrari@Gmail.com –
با سلام و احترام
بسیار خوب. هیچ مشکلی از بابت منطقه مورد مطالعه نیست. حتی میتونید منطقه مورد مطالعه را بصورت دستی تعیین کنید. منطقه مورد مطالعه تابعی از کار شما است نه آموزش بنده.
کدهای ارسالی برای بررسی زیاد هست و شرمنده خیلی از کاربران هم شدیم که نرسیدیم کدشون رو چک کنیم. ممنون میشم دوباره ایمیل ارسال کنید. سعی میکنم در کوتاه مدت بررسی کنم.
موفق باشید
آیسان –
با سلام و احترام
ممنونم از شما و این همه پشتیبانی
من مجددا ایمیل کردم .
Amirhossein.Ahrari@Gmail.com –
با سلام و احترام
ممنون. اگر پاسخ داده نشد مجددا ایمیل کنید چون اینقدر لینک و ایمیل های زیادی دریافت کردم واقعا نمیدونم الان کدوم ایمیل رو شما ارسال کردید.
موفق باشید
آیسان –
سلام یه سوال داشتم از حضورتون
اگه بخوایم از یه مرحله ای به بعد تصاویر رو خروجی گیریم و در محیط انوی یا هر محیط دیگه ای یه سری پردازش هایی رو تصویر انجام بدیم چه کار کنیم که خروجی در قالب یک تیف باشه …آیا به خاطر مساحت بالای منطقه من با این مشکل مواجه میشم….استان خوزستان هست منطقه مطالعاتی ام
ممنونم
Amirhossein.Ahrari@Gmail.com –
با سلام و احترام
باید از همان دستور export image استفاده کنید که در آموزش ها وجود دارد. درصورتی که وسعت منطقه زیاد هست میتوانید maxpixels رو مقدارش رو افزایش بدید مثلا ۱e20
برای اینکه خروجی گرفتن هم راحت تر بشه پیشنهاد میکنم توان تفکیک مکانی رو به ۲۰ تا ۳۰ متر در هنگام ذخیره سازی کاهش دهید تا فرایند آن راحت تر انجام بگیرد.
موفق باشید
آیسان –
ممنونم از پاسختون
ولی ۲۰ رو از من قبول نکرد و ماکزیمم ۱۳ رو تو خطا گفت میشه زد منم ۱e13 رو وارد کردم با تفکیک مکانی ۳۰متر با این وجود باز هم در قالب ۴ تیف به من خروجی داد ….راه حل دیگه ای شما به نظرتون نمیرسه؟
تشکر
Amirhossein.Ahrari@Gmail.com –
بسیار خوب. بدلیل حجم زیاد فایل هست که چهار فایل به شما داده که میتونید در نرم افزار انوی آنها را موازییک کنید.
موفق باشید