مدلهای رقومی ارتفاع (DEM) پیکربندی سطح زمین را به تصویر میکشند و در بسیاری از مناطق در علوم زمین و محیطزیست کاربرد دارند. در این مطالعه، دقت مدل سطح دیجیتال جهان سه بعدی ماهوارهای رصد زمین پیشرفته نسخه 2.1 (ALOS W3D30)، مدل مدل ارتفاعی دیجیتال توپوگرافی رادار شاتل نسخه 3.0 (SRTM30) و گسیل گرمایی و انعکاس رادیومتر پیشرفته نسخه DEM 2.0 (ASTER GDEM2) با استفاده از داده های نظرسنجی GPS با دقت بالا مورد ارزیابی آماری قرار گرفت. خطاهای ریشه-میانگین مربع ~5.40 متر، ~7.47 متر و ~20.03 متر به ترتیب برای ALOS W3D30، SRTM30 و ASTER GDEM2 به دست آمد. در تجزیه و تحلیل های بیشتر، ما متوجه شدیم که ALOS W3D30 و SRTM30 در مناطقی که فواصل ارتفاع در 201 متر – 400 متر و بیش از 801 متر است، بسیار دقیق تر هستند. ثابت شد که ALOS W3D30 دقیقترین DEM است که توپوگرافی سطح زمین را به بهترین شکل نشان میدهد و میتواند برای برخی کاربردهای زمینی و محیطی در نیجریه استفاده شود. ما توصیه می کنیم که این مطالعه باید به عنوان راهنمای استفاده از هر یک از این DEM ها برای کاربردهای زمین و محیط زیست در نیجریه باشد.
کلید واژه ها
ALOS W3D30 ، ASTER GDEM2 ، SRTM30 ، نیجریه ، DEMs ، ارزیابی دقت ، خطای ریشه میانگین مربع
1. مقدمه
مدلهای دیجیتال ارتفاع (نمایش سه بعدی سطح زمین) منابع اصلی اطلاعات ارتفاع هستند که در بسیاری از رشتهها کاربرد زیادی دارند. بسیاری از مناطقی که DEM ها در آنها اعمال می شود عبارتند از: مدل سازی طغیان سیل [ 1 ]; نقشه برداری پوشش گیاهی [ 2 ] [ 3 ] ; نقشه برداری از محیط های صخره مرجانی [ 4 ] ; توسعه مدل های جهانی ژئوپتانسیل [ 5 ] ; ارزیابی تغییر حجم یخچال های طبیعی [ 6 ] ; سیستم های ناوبری برای هوانوردی تجاری [ 7 ] ; مدلسازی اقلیمی [ 8 ] ; باستان شناسی [ 9 ] ; تغییر سطح یخچال طبیعی [ 10 ] ; تجزیه و تحلیل هیدرولوژیکی و شبیه سازی [ 11] ؛ علم خاک و زمین شناسی [ 12 ] ; ژئومورفولوژی و هیدرولوژی حوضه [ 13 ] ; و نظارت بر فرسایش و رسوبات ساحلی [ 14 ]. در مطالعه دیگری، نویسنده بیشتر مناطق مختلفی را که در آن DEM های جهانی یا نزدیک به جهانی می توان اعمال کرد، دسته بندی کرد [ 15 ].
درست است که DEM ها به منابع بسیار مفید داده برای طیف وسیعی از کاربردها در علوم زمین و محیط زیست تبدیل شده اند [ 16 ]، اما علیرغم مفید بودن آنها، منابع بسیاری از خطاهای ذاتی در آنها وجود دارد [ 17 ]] . با توجه به این واقعیت و در نتیجه بهبودهای اخیر و انتشار نسخه های جدیدتر مدل سطح دیجیتال سه بعدی ماهواره جهان پیشرفته رصد زمین نسخه 2.1 (ALOS W3D30). شاتل رادار توپوگرافی مدل ارتفاعی دیجیتال نسخه 3.0 (SRTM30) و رادیومتر گسیل حرارتی پیشرفته و انعکاس فضایی پیشرفته DEM نسخه 2.0 (ASTER GDEM2)، ارزیابی و مقایسه کیفیت این داده ها بسیار مهم است تا ببینید چقدر خوب است. DEM ها با داده های محلی موجود یا به دست آمده مطابقت دارند. این به تعیین اندازه خطاهای آنها در یک منطقه کاربردی کمک می کند.
چندین محقق [ 15 ] [ 17 ] – [ 24 ] از مناطق مختلف جهان، دقت یا تناسب این DEM ها را با توجه به داده های محلی موجود یا به دست آمده آنها ارزیابی کردند و دقیق ترین DEM را در منطقه خود گزارش کردند. جدول 1 نتایج به دست آمده توسط این نویسندگان از کشورهای مختلف جهان را ارائه می دهد. اگرچه ALOS W3D30 بهعنوان دقیقترین DEM در تقریباً همه مطالعات گزارش شده است، نیاز به تعیین سطح واقعی دقت آن در هر منطقه مورد مطالعه وجود دارد زیرا مطالعات سطوح مختلفی از دقت را نشان میدهند. این به این دلیل است که دقت DEM ها به منطقه مورد مطالعه، ماهیت محیط، روش های توسعه الگوریتم، داده های ورودی، پردازش داده ها و وضوح حسگر بستگی دارد. به عنوان مثال ازجدول 1 ، RMSE ALOS W3D30 در فیلیپین، کامرون و روسیه به ترتیب 5.68 متر، 13.06 متر و 7.87 متر است.
در حال حاضر، هیچ نقشه توپوگرافی به راحتی در دسترس نیست که بتواند به راحتی اطلاعات توپوگرافی را برای کاربردهای علمی مختلف در نیجریه ارائه دهد و این یک واقعیت شناخته شده است که به دست آوردن زمینی داده های مکانی، پر زحمت، زمان بر و بسیار پرهزینه تر از انجام همین کار از راه دور است. اگرچه مطالعات متعددی در مورد ارزیابی دقت DEM ها در نقاط مختلف جهان انجام شده است، اما هیچ مطالعه جامعی در مورد دقت عمودی این DEM های آزاد در دسترس در نیجریه وجود ندارد. این در حالی است که این DEM ها به عنوان منابع اصلی اطلاعات توپوگرافی برای کاربردهای متعدد در علوم زمین و محیط زیست مورد استفاده قرار می گیرند. بنابراین این مطالعه هدف آن ارزیابی دقت این DEM ها در نیجریه با استفاده از داده های نظرسنجی GPS (سیستم موقعیت یاب جهانی) به دست آمده از زمین است زیرا روشی مستقل برای ارزیابی کیفیت این سه DEM در نیجریه ارائه می دهد. این اعتبار همچنین بهعنوان بازخوردی برای گروههای تحقیقاتی و/یا سازمانهای دولتی که این DEMها را توسعه دادهاند، عمل میکند و در نظر گرفته شده است که تمام مطالعات دیگری را که در کشورهای دیگر برای ارزیابی کیفیت آنها انجام شده است، تکمیل کند.
2. مواد و روش
2.1. مجموعه داده ها
مجموعه داده های مورد استفاده در این مطالعه عبارتند از: شصت و پنج (65) نقطه GPS، ALOS W3D30، SRTM30 و ASTER GDEM2. شصت و پنج (65) نقطه GPS مختصات ژئودتیکی هستند که بخشی از شبکه ژئودتیک نیجریه را تشکیل می دهند. ارتفاعات بیضی شکل از 22.84 متر تا 1793.41 متر متغیر است. خطاهای ریشه میانگین مربع (RMSE) ارتفاعات بیضی شکل در دوره مرجع، (01. JAN.2012)، از 0.00101 متر تا 0.0244 متر برای شصت و پنج (65) نقطه GPS متغیر است. این نقاط GPS از دفتر نقشه بردار کل فدراسیون (OSGoF) در نیجریه به دست آمده است. شکل 1 توزیع نقاط GPS را در نیجریه نشان می دهد.
ما بخشهای مرتبط ALOS W3D30 DSM [ 25 ]، SRTM30 DEM و ASTER GDEM2 [ 26 ] را در نیجریه دانلود کردیم. مقادیر رستری (ارتفاع) این DEM ها به مختصات نقاط GPS استخراج شد. این ارتفاعات به عنوان ارتفاعات به دست آمده از هر یک از DEM ها گفته می شود. جدول 2 خلاصه ای از ویژگی های DEM های مورد استفاده را نشان می دهد.
شکل 1 . توزیع نقاط GPS در نیجریه
2.2. روشهای ارزیابی دقت
همانطور که در جدول 2 ارائه شده است ، ارتفاعات به دست آمده از این سه DEM به صورت عمودی به مدل گرانشی زمین 1996 (EGM96) ارجاع داده شده است و این منجر به تبدیل ارتفاعات بیضی شکل با استفاده از موج های زمینی محاسبه شده از EGM96 [ 31 ] شد. از نظر ریاضی، رابطه (2) رابطه بین ارتفاع بیضی شکل و ارتفاع مشتق شده EGM96 را نشان می دهد:
ساعتجی پی اس–نEGM96=اچارتوhGPS−NEGM96=HOrtho(1)
جایی که: ساعتجی پی اس= ارتفاع بیضی شکلhGPS=Ellipsoidal height، نEGM96= موج دار شدن ژئوئید برگرفته از EGM 96 ، اچارتو= ارتفاعNEGM96=Geoid undulation derived from EGM96, HOrtho=heightبرگرفته از EGM96.
میانگین خطا (معادله (2))، خطای انحراف معیار (معادله (3))، خطای ریشه میانگین مربع (معادله (4)) و ضریب همبستگی (معادله (5)) ابزارهای آماری هستند که در ارزیابی دقت عمودی ارتفاعات بهدستآمده از این سه DEM که توسط سایر محققین اتخاذ شده است [ 32 ] [ 33 ] [ 34 ] [ 35 ]. تفاوت در ارتفاع به دست آمده از هر یک از نقاط DEM و GPS به عنوان “خطا” نامیده می شود، زیرا نقاط GPS به دست آمده از زمین از دقت بالاتری برخوردار هستند.
ME ( MeanError ) =∑نi = 1( E )نME(MeanError)=∑i=1N(E)N(2)
جایی که؛ E = خطا = اچجی پی اس–اچDEMHGPS−HDEM، اچجی پی اسHGPS= ارتفاعات حاصل از EGM96 از داده های نظرسنجی GPS، اچDEMHDEM= ارتفاع به دست آمده از هر یک از DEM ها، N = تعداد نقاط آزمون.
STDE (خطای انحراف استاندارد) = ∑نi = 1( E – ME )2N – 1———-√∑i=1N(E−ME)2N−1(3)
RMSE (خطای ریشه-میانگین مربع) = ∑نi = 1(E2)ن——-√∑i=1N(E2)N(4)
هرچه مقدار RMSE به صفر نزدیکتر باشد، ارتفاعهای بهدستآمده از DEMها دقیقتر است، در حالی که هر چه مقدار RMSE از صفر دورتر باشد، ارتفاعهای بهدستآمده از DEMها دقت کمتری دارند.
کورل ( X ، ی) =∑را( x- _ایکس¯) ( y–y¯)∑را( x- _ایکس¯)2∑را( y–y¯)2√Correl (X,Y)=∑(x−x¯)(y−y¯)∑(x−x¯)2∑(y−y¯)2(5)
جایی که؛ X = ارتفاعات حاصل از EGM96 از داده های نظرسنجی GPS، Correl = ضریب همبستگی، Y = ارتفاعات به دست آمده از هر یک از DEM ها، ایکس¯ و y¯ نمونه ها .x¯ and y¯ aresamplemeans.
هرچه مقدار ضریب همبستگی به 1± نزدیکتر باشد، سطح توافق ارتفاعات بهدستآمده از هر یک از DEMها به ارتفاعات EGM96 از دادههای نظرسنجی GPS بیشتر است و بالعکس.
علاوه بر این، خطاهای خطی (LE) هر یک از سه DEM با این فرض که خطاهای عمودی معمولاً توزیع می شوند و خطاهای خطی مستقیماً با خطاهای انحراف معیار [ 35 ] [ 36 ] متناسب هستند.
LE @ 90 ٪= 1.6449 × STDELE@90%=1.6449×STDE(6)
LE @ 95 ٪= 1.9000 × STDELE@95%=1.9000×STDE(7)
LE @ 99.73 ٪= 3.0000 × STDELE@99.73%=3.0000×STDE(8)
شصت و پنج (65) نقطه GPS و هر یک از ارتفاعات مربوطه به دست آمده از سه DEM در فواصل ارتفاع 200 متر برای ارزیابی عملکرد فشرده تر DEMs طبقه بندی شدند. فواصل ارتفاع 0 – 200 متر، 201 – 400 متر، 401 – 600 متر، 601 – 800 متر و > 800 متر است. نتایج آماری بهدستآمده از هر یک از این کلاسها برای ارزیابی تأثیر زمینهای موجدار بر دقت عمودی هر یک از سه DEM استفاده شد.
3. نتایج و بحث
موج های زمینی مشتق شده از EGM96 در شصت و پنج (65) ایستگاه GPS در شکل 2 نشان داده شده است. این مقادیر از ارتفاعات بیضی شکل در هر یک از شصت و پنج (65) نقطه GPS کم شد تا ارتفاعاتی که به عنوان ارتفاعات مشتق شده از EGM96 نامیده می شوند به دست آید.
ارتفاعات حاصل از EGM96 از داده های بررسی GPS و ارتفاعات به دست آمده از هر یک از DEM ها در شکل 3 نشان داده شده است.
شکل 2 . موجهای ژئوییدی
(الف)
(ب)
(ج)
شکل 3 . (ac): ارتفاعات (الف) GPS در مقابل ALOS W3D30 (ب) GPS در مقابل SRTM30 (ج) GPS در مقابل ASTER GDEM2.
نگاهی دقیق تر به شکل 3 نشان می دهد که تفاوت شدیدی در ایستگاه های 42، 1 و 40 به ترتیب برای ALOS W3D30، SRTM30 و ASTER GDEM2 وجود دارد. اینها ایستگاه هایی هستند که حداکثر اختلاف ارتفاع را بین هر DEM و نقاط GPS دارند. تفاوت (خطاهای) به دست آمده در ارتفاعات بین هر یک از DEM ها و شصت و پنج (65) نقطه GPS در شکل 4 نشان داده شده است. نتایج آماری این خطاها در جدول 3 نشان داده شده است. ضرایب همبستگی بین ارتفاعات به دست آمده از EGM96 از داده های نظرسنجی GPS و ALOS W3D30، SRTM30 و ASTER GDEM2 به ترتیب 0.9999، 0.9998 و 0.9993 است، به این معنی که هر یک از DEM ها دارای همبستگی بالایی با ارتفاع های GPS9-EGM هستند. امتیاز اما ALOS W3D30 بالاترین سطح توافق را دارد.
از نتایج نشان داده شده در جدول 3 ، می توان استنباط کرد که ALOS W3D30 دقیق تر از STRM30 و ASTER GDEM2 در نیجریه است. به طور کلی، RMSE به دست آمده از ALOS W3D30 40 سانتی متر با 5 متر مشخص شده متفاوت است [ 28 ]. قابل توجه است که SRTM30 به مراتب بهتر از دقت عمودی مشخص شده 16 متر عمل کرد.
علاوه بر این، ارتفاعات به فواصل ارتفاع 200 متری طبقه بندی شدند تا فاصله ارتفاعی را که بهتر با داده های نظرسنجی GPS مشاهده شده محلی مطابقت دارد تشخیص دهند. نتایج آماری مربوط به خطاها در هر بازه ارتفاعی در شکل 5 ارائه شده است در حالی که خطاهای خطی در جدول 4 فهرست شده اند .
در تمام فواصل ارتفاع، ALOS W3D30 در دقت عملکرد بهتری داشت، پس از آن SRTM30 بسیار عقب بود و ASTER GDEM2 بسیار عقب است. همانطور که توسط مطالعات دیگر تایید شد [ 15 ] [ 18 ] [ 22 ] [ 23 ] ALOS W3D30 از نظر دقت بهتر از دو DEM دیگر ارزیابی شده در نیجریه است. این بدان معناست که ALOS W3D30 توپوگرافی سطح زمین را در منطقه مورد مطالعه بهویژه در مناطقی که ارتفاع آنها بیش از 801 متر است، به بهترین شکل نشان میدهد.
افزایش مداوم دقت ALOS W3D30 و SRTM30 در فواصل ارتفاع 0 تا 200 متر وجود دارد. 201 – 400 متر و افزایش شدید در 601 – 800 متر و > 801 متر. این نشان میدهد که مناطق یا ایالتهایی که ارتفاع آنها بین 201 تا 400 متر است، مدلسازی بهتری از زمین توسط ALOS W3D30 نسبت به مناطق بین 0 تا 200 متر دارند. 401 – 600 متر و 601 – 800 متر. بر اساس تمام مقادیر آماری بهدستآمده در فاصله ارتفاعی 601 – 800 متر، ALOS W3D30 و ASTER GDEM2 عملکرد ضعیفی داشتند و ALOS W3D30 بیش از دو برابر دقت مورد انتظار داشت در حالی که ASTER GDEM تقریبا دو برابر دقت مورد انتظار خود را داشت.
اگرچه پنج (5) نقطه GPS در بازه ارتفاع >801 متر قرار داشتند، ALOS W3D30 دو و نیم برابر بهتر از دقت مورد انتظار خود عمل کرد، SRTM30 بیش از سه برابر بهتر از دقت مورد انتظار خود عمل کرد در حالی که ASTRE GDEM2 در دقت مورد انتظار خود باقی می ماند. بدیهی است و همانطور که توسط نویسندگان دیگر [ 15 ] [ 18 ] [ 22 ] [ 23 ] تأیید شده است که دقت عمودی DEM های جهانی به شدت تحت تأثیر شیب زمین قرار می گیرد. خطاهای خطی ( جدول 4) محاسبه شده از هر بازه ارتفاع، برتری ALOS W3D30 را نسبت به دو DEM دیگر ارزیابی شده در این مطالعه تایید کرد و این به وضوح نشان می دهد که ALOS W3D30 می تواند به تنهایی یا در ترکیب با داده های زمینی برای برخی از برنامه های زمینی و محیطی استفاده شود.
شکل 4 . خطاهای ارتفاع به دست آمده از هر یک از DEM ها.
(الف)
(ب)
(ج)
شکل 5 . (الف) نتایج آماری فواصل ارتفاع (الف) خطاهای ریشه – میانگین – مربع (ب) خطاهای انحراف استاندارد (ج) خطاهای میانگین.
4. نتیجه گیری
این مطالعه با هدف ارزیابی کیفیت DEMS جهانی یا نزدیک به جهانی، از چندین ابزار آماری برای تعیین دقت ارتفاع بهدستآمده از ALOS W3D30، ASTER GDEM2 و SRTM30 با استفاده از دادههای پیمایش GPS با دقت بالا در نیجریه استفاده کرد. در تمام تحلیلها، ALOS W3D30 دقیقترین DEM است که میتواند توپوگرافی سطح زمین را در نیجریه به طور نسبی به تصویر بکشد. ما کشف کردیم که مناطق یا ایالتها، در منطقه مورد مطالعه، که در آن فواصل ارتفاع بیش از 801 متر است، نتایج آماری بهتری نسبت به سایرین دارند در حالی که مناطق در فاصله ارتفاعی 601 متر تا 800 متر نتایج آماری بدتری هنگام استفاده از ALOS W3D30 در نیجریه دارند.
ارزیابی دقت DEM ها در علوم زمین و محیط زیست از اهمیت بالایی برخوردار است، زیرا نشان می دهد که چگونه DEM ها به بهترین شکل به سطح زمین دینامیکی تقریب می کنند. به طور کلی، دقت DEM ها به منطقه مورد مطالعه، ماهیت محیط، روش های توسعه الگوریتم، داده های ورودی، پردازش داده ها و وضوح حسگر بستگی دارد. این سطوح مختلف دقت ثبت شده توسط هر یک از DEM ها را توضیح می دهد. بنابراین، ما توصیه می کنیم که این مطالعه باید به عنوان راهنمای استفاده از هر یک از این DEM ها برای کاربردهای زمین و محیط زیست در نیجریه باشد.
بدون دیدگاه