محاسبه سری زمانی پیشرفته مساحت سطحی دریاچه با تصاویر MODIS – در گوگل ارث انجین

272 بازدید

اخبار ، مقالات و تحقیقات گروهی را دنبال کنید.

آموزش رایگان نرم افزار ENVI

در ENVI حرفه ای شو 🦾

0 تا 100 مفاهیم سنجش از دور

کاملا رایگان شروع کن 🧩

آموزش رایگان ArcGIS pro

واقعاااااا رایگاااااانه 🤗

آموزش گوگل ارث انجین

شدیدا توصیه شده 🌏

دانلود لندیوز 10 متری ایران و …

به کمک گوگل ارث انجین 🚀

کتاب موتور پردازش مجازی

اولین کتاب آموزش گوگل ارث انجین 📗

طبقه بندی تصاویر ماهواره ای

آشنایی با انواع روش ها 😲

پیش پردازش تصاویر ماهواره ای

هندسی – رادیومتریکی – اتمسفری 🀄

دانلود داده های آلتیمتری

ترازسنجی برای دریاچه ها 🧿

پوشش اراضی 10 متری

تولید شده توسط ESRI 🧐

نمایش سه بعدی در ArcGIS

مدل رقومی ارتفاع 🗻

مطالب پربازدید

محصول آموزشی محاسبه سری زمانی پیشرفته مساحت سطحی دریاچه با تصاویر MODIS در گوگل ارث انجین منتشر شد! در این محصول آموزشی روش آشکارسازی و شناسایی پهنه های آبی و جداسازی آن از سایر پوشش ها به همراه محاسبه مساحت و پایش تغییرات آنها آموزش داده شده است. با استفاده از این محصول آموزشی شما با روش های پیشرفته سری…

محصول آموزشی محاسبه سری زمانی پیشرفته مساحت سطحی دریاچه با تصاویر MODIS در گوگل ارث انجین منتشر شد!

در این محصول آموزشی روش آشکارسازی و شناسایی پهنه های آبی و جداسازی آن از سایر پوشش ها به همراه محاسبه مساحت و پایش تغییرات آنها آموزش داده شده است. با استفاده از این محصول آموزشی شما با روش های پیشرفته سری زمانی مساحت سطحی آب در گوگل ارث انجین آشنا می شوید. در این محصول، پایش ۲۰ ساله تغییرات مساحت سطحی دریاچه ارومیه آموزش داده شده و برای سایر دریاچه ها نیز کاربردی است. 


مشخصات محصول آموزشی 

  • قیمت: ۸۹/۹۰۰ تومان
  • مدرس: امیرحسین احراری
  • تخصص: کارشناس سنجش از دور و پردازش تصاویر ماهواره ای
  • موضوع: سنجش از دور آب
  • نرم افزار: سامانه گوگل ارث انجین
  • داده های ماهواره ای: سنجنده MODIS
  • پیش نیاز: ندارد.

توضیحات محصول آموزشی 

در این آموزش موضوع محاسبه مساحت سطحی پهنه آبی دریاچه ارومیه به همراه پایش تغییرات آن در ۲۰ سال گذشته با استفاده از داده های سنجنده مادیس تدریس شده است.

با استفاده از این محصول آموزشی شما با پایش تغییرات ۸ روزه و ماهیانه دریاچه ارومیه با استفاده از پروداکت بازتاب سطحی سنجنده مادیس MOD09 با توان تفکیک مکانی ۵۰۰ متری آشنا خواهید شد. از سوی دیگر شما با روش محاسبه روند خطی تغییرات به همراه آشکارسازی تغییرات فصلی آن نیز آشنا می شوید. محاسبه روند خطی تغییرات در محاسبات سری زمانی به شما کمک خواهد کرد تا بتوانید روند ثابت تغییرات را شناسایی نمایید. 


عناوین آموزشی 

عناوین آموزش داده شده در گوگل ارث انجین بر اساس داده های سنجنده MODIS عبارت اند از:

  1. روش فراخوانی پروداکت بازتاب سطحی 
  2. روش ایجاد برش طیفی 
  3. روش ایجاد برش مکانی 
  4. روش محاسبه شاخص آب 
  5. روش جداسازی آب از سایر پوشش ها
  6. روش محاسبه مساحت سطحی 
  7. روش پایش تغییرات مساحت سطحی 
  8. روش محاسبه روند تغییرات خطی 
  9. روش محاسبه روند تغییرات خطی با مدل هارمونیک
  10. روش تبدیل پروداکت ۸ روزه به ۱ ماهه
  11. روش پایش تغییرات ۸ روزه
  12. روش پایش تغییرات ماهیانه

منظور از مساحت سطحی پیشرفته چیست؟

در این محصول آموزشی موضوع مساحت سطحی پهنه های آبی بصورت پیشرفته مورد توجه قرار گرفته است.

منظور از این عنوان چیست؟ 🙄 در حقیقت با استفاده از گوگل ارث انجین شما می توانید پردازش های سری زمانی را به راحتی انجام دهید. در این محصول آموزشی ما کمی موضوع را فراتر بررسی کردیم. با استفاده از این آموزش شما می توانید علاوه بر مساحت سطحی و سری زمانی آن تغییرات خطی مساحت سطحی دریاچه را نیز بررسی کنید.

این موضوعی است که تا کنون در هیچ یک از محصولات آموزشی آن را بررسی نکردیم و ایده های جالبی را به شما می دهد. استفاده از این روش به شما کمک خواهد کرد تا بسیاری از خطاهای موجود در سری زمانی استخراج شده را از بین ببرید. 


مساحت سطحی ۸ روزه سنجنده MODIS

در شکل زیر نمودار تغییرات مساحت سطحی دریاچه ارومیه برای ۲۰ سال اخیر با استفاده از داده های پروداکت بازتاب سطحی نمایش داده شده است. واحد تغییرات مساحتی نمایش داده شده برحسب کیلومتر مربع است.


مساحت سطحی ماهیانه سنجنده MODIS

در شکل زیر تغییرات مساحت سطحی دریاچه ارومیه برای ۲۰ سال اخیر بصورت ماهیانه نمایش داده شده است. در این محصول روش تبدیل نمودار های سری زمانی ۸ روزه به ماهیانه نیز آموزش داده شده است. این تکنیک بسیار کاربردی بوده و نقش موثری در کاهش میزان خطای نتایج دارد.


روند خطی تغییرات مساحت سطحی 

برای محاسبه روند خطی تعییرات مساحت سطحی تکنیک های مهمی در این محصول آموزش داده شده که جزو مطالب پیشرفته گوگل ارث انجین است و در هیچ یک از آموزش های سری زمانی قبلی تا کنون مورد اشاره قرار نگرفته است. با استفاده از روند های خطی، به خوبی روند صعودی و نزولی نمودار ها در طول زمان قابل ارزیابی است.

 


تهیه محصول آموزشی 

برای خرید این محصول آموزش بر روی لینک زیر کلیک کنید. بلافاصله پس از خرید می توانید نسبت به دانلود آن اقدام نمایید. 


لینک های مفید

نوشته های مرتبط :

آموزش های رایگان پیشنهادی :

۵۹ دیدگاه. ارسال دیدگاه جدید

  • مهرانه خدامرادپور
    ۱۳۹۸-۰۵-۲۸ ۱۱:۱۸ ب٫ظ

    سلام استاد محترم، در کدوم نرم فزار میشه نقشه های تولید شده گوگل ارث انجین را از حالت استک خارج کرد و تک تک سیو کرد؟

    پاسخ
  • مهدی احمدی
    ۱۳۹۸-۰۷-۰۵ ۴:۳۰ ب٫ظ

    سلام
    ضمن عرض خسته نباشید

    در روش فوق میتوان بجای تصاویر ماهواره مادیس از ماهواره‌های دیگر مانند لندست یا سنتینل۲ استفاده کرد!؟

    با تشکر

    پاسخ
  • مهدی احمدی
    ۱۳۹۸-۰۷-۰۶ ۲:۰۳ ب٫ظ

    سلام و خسته نباشید

    کدهای مربوطه که فرمودید ذخیره کردید در فایل زیپ موجود نمیباشد

    همچنین در فایل زیپ ساخت فیلم متحرک نیز موجود نبود.

    باتشکر

    پاسخ
    • با سلام و احترام

      زمانی که لینک دانلود محصول به ایمیل شما ارسال شد در آن لینک را قرار داده بودیم. به هر صورت شماره نام آموزش ها رو بفرمایید که 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 این را هم بالا دیتا های لندسات میتوانیم اجرا کنیم یا خیر؟
    ممنون

    پاسخ
  • استاد احراری عزیز!
    کتاب جدید گوگل ارث را چه زمانی بدست آورده میتوانم و ایا شکل پی دی اف ان موجود هست یا خیر؟
    تشکر

    پاسخ
  • سلام علیکم
    ببخشید که باز هم در مورد کود هارمونیک مادل میپرسم.
    // 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

    پاسخ
  • سلام جناب استاد.
    لطفا لینک که برای شما میفرستم . من قصد اجرای یک مادل هارمونیک را با اندیکس NDWI با استفاده از دیتا های لندسات هفت را اجرا کنم اما در مورد نتیجه اش مطمین نیستم یعین اینکه این کود را من درست نوشته ام یا خیر میخواستم شما یک بار کود را بیبیند وتایید کنید.
    https://code.earthengine.google.com/01cfb0d12825ce869da8ee4d24a27ec1

    پاسخ
  • جناب استاد!
    من تشکر می کنم از اینکه سوالات مارا جواب میدهید.
    شما یک لطف کنید همین کود را یک ران کنید در ماه های که ساحه مورد مطالعه پوشیده از برف می باشد مساحت ۵ کیلو متر ولی ما های که برف موجود نیست مساحت محاصبه شد ۰٫۰۲ کیلو متر میباشد . البته این تغیرات در اکسل که برای تغیرات مساحت ماهانه تولید می کند میتوانید بیبینید. اگر گو گل ارث انجنین دارای این نوع خطا باشد باید بدانیم و اگر خطا از محاسبات که در کود موجود هست هم باید بدانیم. لینک زیر اکسل فایل ان را برای تغیرات ماهانه دانلود کنید و تغیرات را مشاهده کنید.
    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 استفاده کنید. در لندست ۸ معادل باند ۵ و در لندست ۷ معادل با باند ۴ است.

      موفق باشید

      پاسخ
      • استاد احراری عزیز!
        من اگر با باند ۴ هم کار کردم یعنی ترکیب از باند ۲ و ۴ نتیجه تغیر پیدا نمود ولی هیچ نتیچه منطقی ارایه نداد. ماه هایکه باید آب زیاد باشد کم هست و ما های که آب باید کم باشد زیاد نشان میدهد. به هر صورت مشکل پا بر جا هست تشکر از همکاری تان.

        احمد

        پاسخ

دیدگاهتان را بنویسید

نشانی ایمیل شما منتشر نخواهد شد.

این فیلد را پر کنید
این فیلد را پر کنید
لطفاً یک نشانی ایمیل معتبر بنویسید.

فهرست