اخبار ، مقالات و تحقیقات گروهی را دنبال کنید.
محاسبه سری زمانی پیشرفته مساحت سطحی دریاچه با تصاویر MODIS – در گوگل ارث انجین
محصول آموزشی محاسبه سری زمانی پیشرفته مساحت سطحی دریاچه با تصاویر MODIS در گوگل ارث انجین منتشر شد! در این محصول آموزشی روش آشکارسازی و شناسایی پهنه های آبی و جداسازی آن از سایر پوشش ها به همراه محاسبه مساحت و پایش تغییرات آنها آموزش داده شده است. با استفاده از این محصول آموزشی شما با روش های پیشرفته سری…
محصول آموزشی محاسبه سری زمانی پیشرفته مساحت سطحی دریاچه با تصاویر MODIS در گوگل ارث انجین منتشر شد!
در این محصول آموزشی روش آشکارسازی و شناسایی پهنه های آبی و جداسازی آن از سایر پوشش ها به همراه محاسبه مساحت و پایش تغییرات آنها آموزش داده شده است. با استفاده از این محصول آموزشی شما با روش های پیشرفته سری زمانی مساحت سطحی آب در گوگل ارث انجین آشنا می شوید. در این محصول، پایش ۲۰ ساله تغییرات مساحت سطحی دریاچه ارومیه آموزش داده شده و برای سایر دریاچه ها نیز کاربردی است.
مشخصات محصول آموزشی
- مدرس: امیرحسین احراری
- تخصص: کارشناس سنجش از دور و پردازش تصاویر ماهواره ای
- موضوع: سنجش از دور آب
- نرم افزار: سامانه گوگل ارث انجین
- داده های ماهواره ای: سنجنده MODIS
- پیش نیاز: ندارد.
تهیه محصول آموزشی
برای خرید این محصول آموزش بر روی دکمه زیر کلیک کنید. بلافاصله پس از خرید می توانید نسبت به دانلود آن اقدام نمایید.
توضیحات محصول آموزشی
در این آموزش موضوع محاسبه مساحت سطحی پهنه آبی دریاچه ارومیه به همراه پایش تغییرات آن در ۲۰ سال گذشته با استفاده از داده های سنجنده مادیس تدریس شده است.
با استفاده از این محصول آموزشی شما با پایش تغییرات ۸ روزه و ماهیانه دریاچه ارومیه با استفاده از پروداکت بازتاب سطحی سنجنده مادیس MOD09 با توان تفکیک مکانی ۵۰۰ متری آشنا خواهید شد. از سوی دیگر شما با روش محاسبه روند خطی تغییرات به همراه آشکارسازی تغییرات فصلی آن نیز آشنا می شوید. محاسبه روند خطی تغییرات در محاسبات سری زمانی به شما کمک خواهد کرد تا بتوانید روند ثابت تغییرات را شناسایی نمایید.
عناوین آموزشی
عناوین آموزش داده شده در گوگل ارث انجین بر اساس داده های سنجنده MODIS عبارت اند از:
- روش فراخوانی پروداکت بازتاب سطحی
- روش ایجاد برش طیفی
- روش ایجاد برش مکانی
- روش محاسبه شاخص آب
- روش جداسازی آب از سایر پوشش ها
- روش محاسبه مساحت سطحی
- روش پایش تغییرات مساحت سطحی
- روش محاسبه روند تغییرات خطی
- روش محاسبه روند تغییرات خطی با مدل هارمونیک
- روش تبدیل پروداکت ۸ روزه به ۱ ماهه
- روش پایش تغییرات ۸ روزه
- روش پایش تغییرات ماهیانه
منظور از مساحت سطحی پیشرفته چیست؟
در این محصول آموزشی موضوع مساحت سطحی پهنه های آبی بصورت پیشرفته مورد توجه قرار گرفته است.
منظور از این عنوان چیست؟ 🙄 در حقیقت با استفاده از گوگل ارث انجین شما می توانید پردازش های سری زمانی را به راحتی انجام دهید. در این محصول آموزشی ما کمی موضوع را فراتر بررسی کردیم. با استفاده از این آموزش شما می توانید علاوه بر مساحت سطحی و سری زمانی آن تغییرات خطی مساحت سطحی دریاچه را نیز بررسی کنید.
این موضوعی است که تا کنون در هیچ یک از محصولات آموزشی آن را بررسی نکردیم و ایده های جالبی را به شما می دهد. استفاده از این روش به شما کمک خواهد کرد تا بسیاری از خطاهای موجود در سری زمانی استخراج شده را از بین ببرید.
مساحت سطحی ۸ روزه سنجنده MODIS
در شکل زیر نمودار تغییرات مساحت سطحی دریاچه ارومیه برای ۲۰ سال اخیر با استفاده از داده های پروداکت بازتاب سطحی نمایش داده شده است. واحد تغییرات مساحتی نمایش داده شده برحسب کیلومتر مربع است.
مساحت سطحی ماهیانه سنجنده MODIS
در شکل زیر تغییرات مساحت سطحی دریاچه ارومیه برای ۲۰ سال اخیر بصورت ماهیانه نمایش داده شده است. در این محصول روش تبدیل نمودار های سری زمانی ۸ روزه به ماهیانه نیز آموزش داده شده است. این تکنیک بسیار کاربردی بوده و نقش موثری در کاهش میزان خطای نتایج دارد.
روند خطی تغییرات مساحت سطحی
برای محاسبه روند خطی تعییرات مساحت سطحی تکنیک های مهمی در این محصول آموزش داده شده که جزو مطالب پیشرفته گوگل ارث انجین است و در هیچ یک از آموزش های سری زمانی قبلی تا کنون مورد اشاره قرار نگرفته است. با استفاده از روند های خطی، به خوبی روند صعودی و نزولی نمودار ها در طول زمان قابل ارزیابی است.
تهیه محصول آموزشی
برای خرید این محصول آموزش بر روی دکمه زیر کلیک کنید. بلافاصله پس از خرید می توانید نسبت به دانلود آن اقدام نمایید.
لینک های مفید
نوشته های مرتبط :
آموزش های رایگان پیشنهادی :
آموزش ثبت نام در USGS با جزئیات کامل
1,947 بازدید
خروجی گرفتن نقشه در جی آی اس
2,296 بازدید
وبینار جامع کاربردهای زبان R در سنجش از دور
1,202 بازدید
71 دیدگاه. ارسال دیدگاه جدید
سلام استاد محترم، در کدوم نرم فزار میشه نقشه های تولید شده گوگل ارث انجین را از حالت استک خارج کرد و تک تک سیو کرد؟
با سلام و احترام
در نرم افزار ENVI. با استفاده از دستور resize data و بخش spectral subset.
موفق باشید
سلام
ضمن عرض خسته نباشید
در روش فوق میتوان بجای تصاویر ماهواره مادیس از ماهوارههای دیگر مانند لندست یا سنتینل۲ استفاده کرد!؟
با تشکر
با سلام و احترام
بله قابل استفاده است. کدها بصورت پویا طراحی شده و برای تمامی ماهواره ها می توانید استفاده کنید.
موفق باشید
سلام و خسته نباشید
کدهای مربوطه که فرمودید ذخیره کردید در فایل زیپ موجود نمیباشد
همچنین در فایل زیپ ساخت فیلم متحرک نیز موجود نبود.
باتشکر
با سلام و احترام
زمانی که لینک دانلود محصول به ایمیل شما ارسال شد در آن لینک را قرار داده بودیم. به هر صورت شماره نام آموزش ها رو بفرمایید که example های چند هستن تا بنده براتون ایمیل کنم.
موفق باشید
سلام
در آموزش داده شده از فیلتر ابر( بعنوان مثال کمتر از ۱۰ درصد) میشود استفاده نمود!؟
تشکر از پاسخهای خوبتان
با سلام و احترام
اون فیلتر برای تصاویر لندست کاربردی است. برای مادیس باید کدهای نسبتا پیچیده تری برای این کار نوشت. در حال حاضر کد مناسبی متاسفانه ندارم. با این حال پروداکت های چند روزه میزان ابر بسیار کمتری دارند. یا تابع محاسبه پروداکت ماهیانه که آموزش دادم درصد ابر رو به حداقل می رساند.
موفق باشید
سلام
اگر امکانش هست تابع محاسبه پروداکت فصلی و سالانه را نیز توضیح فرمائید
تشکر فراوان
با سلام و احترام
بنده در همین آموزش روش تولید پروداکت ماهیانه رو با استفاده از تابع monthly product آموزش دادم. آیا آن را مشاهده کردید؟
موفق باشید
سلام
بله ماهیانه را کاملا فرا گرفتم،《 البته آموزش مانند همیشه عالی بود》، اما درخواستم آموزش جهت تولید پروداکت فصلی و همچنین تولید پروداکت سالیانه میباشد
البته اگر امکانش میباشد.
با تشکر
همان تابعی که برای monthly product نوشتیم فقط interval رو بزارید ۹۰ روز میشه قصلی. ۳۶۵ روز بزارید میشه سالیانه.
موفق باشید
سلام
خیلی ممنون از راهنمایی
کمی فرمول را بررسی و دستکاری کردم متوجه شدم در بخش monthly product در قسمت interval عدد ۱ گذاشته، اگر در بخش units واژه ‘month’ بنویسیم خروجی ماهیانه شده و دقیقا میانگین محاسبه شده را در اول هر ماه نمایش میدهد همچنین جهت نمایش خروجی سالانه نیز در بخش units واژه ‘year’ بگذاریم میانگین محاسبه شده دقیقا در اول هر سال نمایش داده میشود.
در خصوص نمایش خروجی فصلی واحدی برای قسمت units در حال حاضر پیدا نکردم لذا از روشی که جنابعالی توضیح فرمودید باید استفاده نمود
با تشکر
سلام مجدد
در خصوص نمایش خروجی فصلی نیز میتوان در قسمت units از واژه ‘month’ استفاده کرده، در بخش interval عدد ۳ را وارد نمود بر همین اساس برای خروجی شش ماهه نیز بجای عدد ۳ در بخش interval عدد ۶ را باید وارد نمود.
البته عذر خواهی میکنم زیاد پیام میگذارم. مطالب ذکر شده توسط بنده را درس پس دادن در محضر استاد تلقی فرمائید.
یا علی
بسیار عالی. خیلی خوب هست که به راهکار راحت تری رسیدید. ممنون از اینکه به اشتراک گذاشتید.
موفق باشید
سلام جناب احراری، یک سوال در مورد تصاویر رادار داشتم. برای محاسبه میزان فرونشست از طریق سری زمانی آیا میشود از سامانه گوگل ارث انجین استفاده نمود؟
با سلام و احترام
امکان استفاده هست اما بصورت تابع آماده دستورات آن نیست و باید تمامی مراحل رو فرمول نویسی کنید که مقرون به صرفه نیست. به همین دلیل پیشنهاد نمیکنم برای این کار از آن استفاده کنید.
موفق باشید
با سلام و احترام
جناب احراری بنظر بنده در کدنویسی شما یک اشکال وجود دارد که باعث می شود مساحت محاسبه شده کمتر از مقدار واقعی بدست آید.
در فانکشن نوشته شده در ابتدای برنامه باید
var mask= index. updateMask(thr
به
var mask=thr.updateMask(thr
تغییر پیدا کند زیرا شما در ادامه مساحت هر پیکسل رو در مقدار داخل هر پیکسل ضرب می کند در کد شما متغییر index با thr اپدیت شده که نتیج آن پیکسل های کمتر از ۰٫۳ حذف می شوند ولی مقدار سلول های بزرگ تر از ۰٫۳ یک نمی شود لذا در محا سبه مساحت هر پیکسل در عدد هر پیکسل که بین ۰٫۳ و ۱ هست ضرب می شود. اما اگر متغییر thr با خود thr اپدیت ماسک شود پیکسل های صفر حذف شده و فقط پیکس های یک باقی می ماند که در ادامه در مساحت پیکسل ها ضرب می شوند و در نهایت جمع می شوند.
ممنون از سایت خوبتون
با سلام و احترام خدمت شما.
ممنون از پیام شما. حتما بررسی میکنم.
موفق باشید
سلام کدش خیلی طولانیه میشه لطفا برای من هم ایمیل کنید example56-1 من هر اموزشی رو خریدم کد داخلش نیست و فقط فیلمه ممنون میشم کدها رو هم بزارید . بازم ممنون
سلام ممنون از زحماتتون کدش خیلی طولانیه میشه لطفا برای من هم ایمیل کنید example56-1 من هر اموزشی رو خریدم کد داخلش نیست و فقط فیلمه ممنون میشم کدها رو هم بزارید . بازم ممنون
ببخشید لینک داخل ایمیل بود من همیشه از همون صفحه خرید دانلود میکردم و ایمیل را باز نمیکردم الان که فرمودید ایمیل را چک کردم داخلش بود.
بله درسته. لینک کد ضمیمه شده است. الان پس به کد اصلی نیز دسترسی دارید.
موفق باشید
بله ممنونم . من همیشه تو اموزشها خودم کد رو از روی فیلمی که میگفتید تایپ میکردم ممنون که راهنمایی کردید.
سلام مجدد ببخشید میشه نام کد محصولی از سنتینل که میشه با اون مساحت پهنه آبی را به دست آورد بهمون بگید منطقه من خیلی کوچیکه و مادیس با پیکسل سایز ۵۰۰ متر خیلی براش خطا داره. ممنونم
با سلام و احترام
در این زمینه آموزشی فعلا تولید نشده. نزدیک ترین محصول به هدف شما آموزشی است که برای شناسایی سیل در تصاویر ماهواره ای تولید کردیم. در بخش جستجوی سایت عبارت سیل را وارد کنید تا آن را نمایش دهد.
موفق باشید
باسلام و خسته نباشید
در این اموزش در اخر بعد از محاسبه نمودارها و تغییرات مساحت نیاز به صحت سنجی داده ها می باشد برای این کار ایا روش یا فیلمی را پیشنهاد می کنید لطفا
با تشکر
با سلام و احترام
برای صحت سنجی داده های داخل اکسل رو باید با مقادیر زمینی مساحت و یا تراز مقایسه کرده و سپس RMSE و Correlation را بگیرید. این دو روش مرسوم ترین روش های اعتبارسنجی این دسته از تکنیک ها هستند.
موفق باشید
سلام
شما برای استفاده از NDWI از باند آبی استفاده کردید،در کنار باندnir
میخواستم بدونم که تو مقالات اکثرا از باند سبز استفاده کردند یا ویا SWIR و این که منبعی برای استفاده از این باندی که تو این محصول استفاده کردید معرقی کنید ممنون میشم.
باتشکر
در واقع باند آبی منظورم هستش که برای این باند اگه منبعی معرفی کنید ممنون میشم.
با سلام و احترام
بنده بصورت تجربی از این رابطه نتیجه گرفتم. اگر به منابع بخواهیم توجه کنیم فرمایش شما درست است و باند سبز و مادون قرمز نزدیک بیشتر استفاده شده است. اگر می خواهید رفرنس دهید من پیشنهاد میکنم از رابطه باند سبز و مادون قرمز نزدیک استفاده کنید که در مقالات با عنوان شاخص mc feeter شناخته شده است.
موفق باشید
سلام آقای احراری
من یک تصویر ndwi ایجاد کردم که ترکیبی از ndwi های یک دوره زمانیه. یعنی تصویری دارم چند بانده، هر باند ndwi یک تایم رو دربرمی گیره. به هر باندش هم تاریخ مربوطه را اختصاص دادم با ()rename. حالا توی arcGIS اکسپورتش کردم ولی اسم باندها رو نشون نمیده و این طوری هم کار کردن با تصویر خیلی سخته. راهی رو پیشنهاد میدین که بتونم توی arcGIS نام باندهای یک تصویر رو ببینم؟
ممنون
با سلام و احترام
اگر داخل qgis باز کنید اسم باندها نشان داده میشه. در غیر این صورت هریک از باندها باید بصورت جداگانه ذخیره سازی شود و به هر تصویر تاریخ متناظر آن را تعریف کرد.
موفق باشید
سلام
کد استفاده شده در این درس (example 56.1) رو میشه ارسال کنید؟ در بسته خریداری شده موجود نبود
با سلام و احترام
لینک کد مورد نظر برای شما ارسال شد.
موفق باشید
ا ستاد بزگوار احراری!
https://code.earthengine.google.com/3984214f1adca455915c8496c0d630b3
من این محصول را ا ز شما خریداری کردم و دقسمت سنجش سطح آب های سحظی در این بخش از آموزش از تصاویر مادیس استفاده نمودید.
ایا امکان دارد که در این بخش تصاویر لندسات را به جای مادیس استفاده کنم. و اگر امکان داشته باشد کدام قسمت های کود را باید تغیر بدهم.
لطفا همکاری تانرا خواهانم.
البته این سوال را من در ایمیل از شما هم پرسدم.
تشکر
با سلام و احترام
لینک ضمیمه شده را ملاحظه کنید:
https://code.earthengine.google.com/054d2aecd2390b30cc5b89b622c3b819
در کد بالا نمونه لندست محاسبه شده است. متناسب با این نمونه سایر بخش های کد را می توانید بر اساس لندست بازنویسی کنید.
لندست برای برخی از تاریخ ها عدد ۰ نشان میدهد. باید داخل اکسل این ۰ ها را حذف کنید.
محاسبه درست انجام شده اما ارث انجین برای برخی تاریخ ها از تصاویر لندست دوتا تصویر میدهد که یکی معادل با صفر و دیگری معادل با مساحت درست است.
داخل اکسل باید اعداد صفر را حذف کنید.
سلام و درود!
جناب استاد تشکر . چند نکته را که میخواستم بدانم.
۱) Harmonic Model: original and fitted values این را هم بالا دیتا های لندسات میتوانیم اجرا کنیم یا خیر؟
ممنون
بله می توانید اجرا کنید.
موفق باشید
استاد احراری عزیز!
کتاب جدید گوگل ارث را چه زمانی بدست آورده میتوانم و ایا شکل پی دی اف ان موجود هست یا خیر؟
تشکر
سلام
در روز وبینار همه جزئیات شرح داده می شود
با سلام و احترام
از روز ۱۸ دی ماه ساعت ۹ شب می توانید در سایت این کتاب را خریداری کنید. لطفا حتما در وبینار رونمایی از چاپ جدید کتاب شرکت کنید.
https://girs.ir/gee-book-unvelling/
موفق باشید
استاد احراری!
اگر در ویبنار شرکت کرده نتوانستم ولی حتمن کتاب را باید تهیه کنم چون خیلی منتظر این کتاب تان بودم.
تشکر
با سلام و احترام
خواهش میکنم هر طور صلاح میدونید. انشاءلله براون مفید باشه.
موفق باشید
سلام علیکم
ببخشید که باز هم در مورد کود هارمونیک مادل میپرسم.
// harmonic model
var modis = ee.ImageCollection(‘LANDSAT/LE07/C01/T1_SR’);
// dependent variable name
var dependent = ‘NDWI’;
// harmonic cycles
var harmonics = 1;
var harmonicFrequencies = ee.List.sequence(1, harmonics);
// band names
var constructBandNames = function(base, list){
return ee.List(list).map(function(i){
return ee.String(base).cat(ee.Number(i).int());
})};
// harmonic model name
var cosName = constructBandNames(‘cos_’,harmonicFrequencies);
var sinName = constructBandNames(‘sin_’,harmonicFrequencies);
// dependent variable list
var NDWI = function(img){
var scale = img.multiply(0.0001);
var clip = scale.clip(geometry);
var index = clip.normalizedDifference([‘B2′,’B5’]);
var thr = index.gte(0.2);
//var mask = index.updateMask(thr);
var area = thr.multiply(ee.Image.pixelArea().divide(1000000.0));
return area.rename(‘NDWI’).float()
.copyProperties(img,[‘system:time_start’,’system:time_end’]);
};
// add time bands
var addDependents = function(img){
var years = img.date().difference(‘2000-01-01′,’year’);
var timeRadians = ee.Image(years.multiply(2 * Math.PI)).rename(‘t’);
var constant = ee.Image(1);
return img.addBands(constant).addBands(timeRadians.float());
};
// sin and cos for time
var addHarmonics = function(freq){
return function(img){
var frequencies = ee.Image.constant(freq);
var time = ee.Image(img).select(‘t’);
var cosines = time.multiply(frequencies).cos().rename(cosName);
var sines = time.multiply(frequencies).sin().rename(sinName);
return img.addBands(cosines).addBands(sines);
};
};
// dependent variables finalized
var harmonicModis = modis
.filterBounds(geometry)
.filterDate(‘2000-01-01′,’2020-01-01’)
.map(NDWI)
.map(addDependents)
.map(addHarmonics(harmonicFrequencies));
// independents
var independents = ee.List([‘constant’,’t’])
.cat(cosName).cat(sinName);
var harmonicTrend = harmonicModis
.select(independents.add(dependent))
.reduce(ee.Reducer.linearRegression(independents.length(),1));
// finalized coefficient
var harmonicTrendCoefficient = harmonicTrend.select(‘coefficients’)
.arrayProject([0])
.arrayFlatten([independents]);
// compute fitted value
var fittedHarmonic = harmonicModis.map(function(image){
return image.addBands(
image.select(independents)
.multiply(harmonicTrendCoefficient)
.reduce(‘sum’)
.rename(‘fitted’)
);
});
آیا کود بالا آلان برای لندست هفت درست هست یا خیر. چون در آموزش برای مادیس نوشته شده اگر لطف کنید یک بار بیبیند که تغیرات ضروری که اوردم برای لندسات کافی هست یا خیر!
تشکر
با سلام و احترام
لطفا همین کد را در ارث انجین ران کنید اگرخطار داد لینک کد را برای بنده در این جا کامنت بزارید تا بررسی کنم.
موفق باشید
سلام جناب استاد.
لطفا لینک که برای شما میفرستم . من قصد اجرای یک مادل هارمونیک را با اندیکس NDWI با استفاده از دیتا های لندسات هفت را اجرا کنم اما در مورد نتیجه اش مطمین نیستم یعین اینکه این کود را من درست نوشته ام یا خیر میخواستم شما یک بار کود را بیبیند وتایید کنید.
https://code.earthengine.google.com/01cfb0d12825ce869da8ee4d24a27ec1
تشکر
سلام استاد احرارای!
در این لینک توجه شما را به نکته جلب میکنم که ممکن اشتباه باشد. مساحت ماهانه در این lake بند امیر افغانستان را محاسبه می کردم . اینجا یکی از جاذبه های گردشگری هست. اما مساحت اب این lake در ماه های جنوری ۵ کیلومتر مربع و لی در ماه های که باید اب این جا اضافه باشد زیر صفر نشان میدهد. مثلی ماه های جون جولای و اگست و سپتمبر تماما تحت صفر محاسبه میشود. لطفا یک بار ای کود ر ا ران کنید و اکلسل ماهانه مساحتش را بیبینید.
https://code.earthengine.google.com/a8fd7f04f3250a3f6f8350832beb12f9
سلام استاد احرارای!
در این لینک توجه شما را به نکته جلب میکنم که ممکن اشتباه باشد. مساحت ماهانه در این lake بند امیر افغانستان را محاسبه می کردم . اینجا یکی از جاذبه های گردشگری هست. اما مساحت اب این lake در ماه های جنوری ۵ کیلومتر مربع و لی در ماه های که باید اب این جا اضافه باشد زیر صفر نشان میدهد. مثلی ماه های جون جولای و اگست و سپتمبر تماما تحت صفر محاسبه میشود. لطفا یک بار ای کود ر ا ران کنید و اکلسل ماهانه مساحتش را بیبینید.
https://code.earthengine.google.com/a8fd7f04f3250a3f6f8350832beb12f9
با سلام و احترام
در بخش mask تغییری رو ایجاد کردم. در حقیقت thr باید updatemask بر اساس thr شود.
https://code.earthengine.google.com/efadba61b5f1a7a6709fd69df2796f08
موفق باشید
سلام جناب استاد.
لطفا لینک که برای شما میفرستم . من قصد اجرای یک مادل هارمونیک را با اندیکس NDWI با استفاده از دیتا های لندسات هفت را اجرا کنم اما در مورد نتیجه اش مطمین نیستم یعین اینکه این کود را من درست نوشته ام یا خیر میخواستم شما یک بار کود را بیبیند وتایید کنید.
https://code.earthengine.google.com/01cfb0d12825ce869da8ee4d24a27ec1
با سلام و احترام
در بخش mask تغییری رو ایجاد کردم. در حقیقت thr باید updatemask بر اساس thr شود.
https://code.earthengine.google.com/efadba61b5f1a7a6709fd69df2796f08
موفق باشید
جناب استاد!
من تشکر می کنم از اینکه سوالات مارا جواب میدهید.
شما یک لطف کنید همین کود را یک ران کنید در ماه های که ساحه مورد مطالعه پوشیده از برف می باشد مساحت ۵ کیلو متر ولی ما های که برف موجود نیست مساحت محاصبه شد ۰٫۰۲ کیلو متر میباشد . البته این تغیرات در اکسل که برای تغیرات مساحت ماهانه تولید می کند میتوانید بیبینید. اگر گو گل ارث انجنین دارای این نوع خطا باشد باید بدانیم و اگر خطا از محاسبات که در کود موجود هست هم باید بدانیم. لینک زیر اکسل فایل ان را برای تغیرات ماهانه دانلود کنید و تغیرات را مشاهده کنید.
https://code.earthengine.google.com/efadba61b5f1a7a6709fd69df2796f08
شاید مشکل از فرمول شاخص و حدآستانه شما باشد.
بنده این کد محاسبه مساحت سطحی آب رو برای سدها و دریاچه ها (با فرمولی متفاوت از فرمول شما) استفاده کردم و دقت RMSE آن کمتر از ۰٫۵ بوده. از لحاظ ساختار کد مطمئن هستم چون در حال حاضر هم در کارهام استفاده میکنم. مثلا در فرمول شما از باند سبز و مادون قرمز طول موج کوتاه استفاده شده که بنده هیچ وقت چنین فرمولی را در آموزش ها اشاره نکردم.
موفق باشید
استاد محترم!
For Landsat 7 data, NDWI = (Band 2 – Band 5) / (Band 2 + Band 5)
من این فورمول را مطابق تدریس شما استفاده نمودم شما در سنجش مساحت آب با استفاده از شاخص ان دی دبلیو آی . اگر مهربانی کنید یک بار دیگر در مورد تغیرات ماهانه سطح آب دریاجه بند امیر که چی مشکال زا خواهد بود.
تشکر استاد
عزیز در ماهواره لندست ۷ باند ۵ معادل با SWIR1 هست. بنده برای آب هیچ گاه از این فرمول استفاده نکردم. باند ۵ رو برای ماهواره لندست ۸ استفاده میکنم که معادل با NIR هست.
برای لندست ۷ به جای باند۵ باید از باند ۴ استفاده کنید. باند ۴ در لندست ۷ معادل با NIR است.
موفق باشید
استاد احراری عزیز!
MNDWI = (Green – SWIR) / (Green + SWIR)
For Landsat 7 data, NDWI = (Band 2 – Band 5) / (Band 2 + Band 5)
این هم مشخصات باند های لندسات هفت دیتا های سر فیس رفلیکتانیس
B2 ۰٫۰۰۰۱ ۰٫۵۲-۰٫۶۰ μm
Band 2 (green) surface reflectance
B5 ۰٫۰۰۰۱ ۱٫۵۵-۱٫۷۵ μm
Band 5 (shortwave infrared 1) surface reflectance
https://developers.google.com/earth-engine/datasets/catalog/LANDSAT_LE07_C01_T1_SR#bands
یعنی الان من از کدام باندها استفاده کنم تا تغیرات ماهانه را صفر نشان ندهد.
تشکر از همکاری تان
خب این فرمول را بنده پیشنهاد نداده ام.
هر فرمولی در مکان های مختلف نتایج گوناگونی دارند. من پیشنهاد میکنم به جای SWIR از باند NIR استفاده کنید. در لندست ۸ معادل باند ۵ و در لندست ۷ معادل با باند ۴ است.
موفق باشید
استاد احراری عزیز!
من اگر با باند ۴ هم کار کردم یعنی ترکیب از باند ۲ و ۴ نتیجه تغیر پیدا نمود ولی هیچ نتیچه منطقی ارایه نداد. ماه هایکه باید آب زیاد باشد کم هست و ما های که آب باید کم باشد زیاد نشان میدهد. به هر صورت مشکل پا بر جا هست تشکر از همکاری تان.
احمد
عرض سلام و ادب
من عین نوشته های جنابعالی رو می زنم ولی متأسفانه Eroor میده. ممنون میشم راهنمایی بفرمایید.
https://code.earthengine.google.com/377dcda8648378e72c362ab4e46ea027
با سلام و احترام
لینک کد اصلاح شده ضمیمه شده است.
https://code.earthengine.google.com/e19524a4784851369626730489bf1e2f
موفق باشید
سلام و عرض ادب
ببخشيد من اين آموزش رو براي داده هاي لندست 8 كالكشن 2(LANDSAT/LC08/C02/T1_L2) تغيير دادم ولي براي اطمينان از تغيير درست و نتيجه ي به دست آمده ممنون ميشم كد زير رو بررسي كنيد و ببينيد با توجه به خروجي نمودار هارمونيك پاياني آيا صحيح است يا خير؟
تشكر
https://code.earthengine.google.com/5b2dfd852959d0eaf6cd6bf6054bd6ad
با سلام و احترام
بدلیل تعدد زیاد سوالات فقط فرصت بررسی خطاهای موجود در کدها رو داریم. نتیجه حاصل از محاسبه را با نتیجه بدست أمده در آموزش مقایسه کنید اگر انطباق داشت مشکلی ندارد.
موفق باشید
سلام و عرض ادب
ببخشيد من كد زير رو براي داده هاي لندست 8 كالكشن 2 تغيير دادم ممنون ميشم بررسي كنيد آيا تغييرات درست اعمال شده يا خير؟
https://code.earthengine.google.com/52d8292fd24698c095db771db4f2e5ba
با سلام و احترام
بدلیل تعدد زیاد سوالات فقط فرصت بررسی خطاهای موجود در کدها رو داریم. نتیجه حاصل از محاسبه را با نتیجه بدست أمده در آموزش مقایسه کنید اگر انطباق داشت مشکلی ندارد.
موفق باشید
سلام و عرض ادب
من اين محصول آموزشي را براي داده هاي لندست 8 كالكشن 2(LANDSAT/LC08/C02/T1_L2) تغيير دادم، براي صحت تغيير كد ممنون ميشم كد زير رو بررسي كنيد و ببينيد با توجه به خروجي آيا تغييرات درست اعمال شده يا خير؟
تشكر
https://code.earthengine.google.com/52d8292fd24698c095db771db4f2e5ba
با سلام و احترام
بدلیل تعدد زیاد سوالات فقط فرصت بررسی خطاهای موجود در کدها رو داریم. نتیجه حاصل از محاسبه را با نتیجه بدست أمده در آموزش مقایسه کنید اگر انطباق داشت مشکلی ندارد.
موفق باشید
سلام مهندس
وقتی می خواهم تغییرات سطح آب دریاچه را محاسبه کنم با توجه به آموزش شما در قسمت ;()var modisWaterBody=modis.map(waterBody).sum وقتی خروجی می گیرم run نمی شود و ارور می دهد ولی برای درصد موردی ندارد می شود راهنمایی کنید. باسپاس
https://code.earthengine.google.com/bc8c9e44b73190f8cf0ce4b77066ddf2
با سلام و احترام
لینک کد اصلاح شده ضمیمه شده است.
https://code.earthengine.google.com/d81abf919236950bdd7f926c40080cf8
موفق باشید
سلام من برای یک بازه زمانی 20 ساله میخوام کد بنویسم تا اختلاف مساحت سد ها و دریاچه ها رو به دست بیارم آیا با این با آموزش و کد های نوشته شده امکانش هست؟ برای رودخانه های و تلاب های فصلی نیز کاربرد دارد؟ شما چه بازه زمانی رو برای تغییر مساحت پهنه های آبی توصیه می کنید
با سلام و احترام
آموزش های زیر مناسب هست
https://girs.ir/gee-dam-monitoring/
https://girs.ir/water-area-time-series-in-google-earth-engine/
موفق باشید