همپوشانی رستر

همپوشانی رستر چندین لایه رستری را در رستر واحد ترکیب می‌کند که در آن هر سلول از هر رستر به موقعیت جغرافیایی مشابهی اشاره می‌کند. بنابراین همپوشانی رستری مانند تقابل بردار شامل تقاطع هندسی نمی‌شود. این شامل عملیات روی مقادیر جداگانه سلول‌های مربوطه در لایه‌های ورودی برای ایجاد لایه رستری جدید با مقادیر سلول جدید است. بر اساس ترکیبات منحصر به فرد مقادیر چند رستر ورودی می‌توان به هر سلول در رستر خروجی مقداری اختصاص داد. معمولاً از همپوشانی رستری برای ترکیب ریاضی رسترهای اختصاص داده شده با مقادیر سلول عددی برای محاسبه مقدار جدید برای هر سلول در رستر خروجی استفاده می‌شود. شکل 4-17 نمونه ای از همپوشانی رستری را با افزودن نشان می‌دهد. دو رستر ورودی با هم جمع می‌شوند تا رستر خروجی با مقادیر هر سلول جمع شوند.

از همپوشانی رستری برای انجام عملیات حسابی، منطقی، رابطه ای، شرطی و آماری در رسترهای ورودی استفاده می‌شود. طبق جدول 4-5 هنگام اعمال عملگرهای منطقی در همپوشانی رستری اغلب ابتدا رسترهای ورودی را در رسترهای دوتایی با مقادیر سلول 1 یا 0 طبقه بندی می‌کنیم، سپس عملیات همپوشانی را انجام می‌دهیم همانطور که در تصویر شکل 4-18 نشان داده شده است. کادر 4-10 نحوه انجام همپوشانی رستری در ArcGIS را نشان می‌دهد.

شکل 4-17 همپوشانی رستر با افزودن

شکل 4-18 پوشش رستری با عملیات منطقی

 

کادر 4-10 همپوشانی رستر در ArcGIS

کاربردی

برای پیروی از این مثالArcMap  را راه اندازی کنید و کلاس ویژگی خاک‌ها را دانلود کنید و landcv_r رستر از مسیر زیر دانلود کنید. پسوند Spatial Analyst را فعال کنید.
C:\Databases\GIS4EnvSci\VirtualCatchment\Geodata.gdb.

پوشش رستری بر اساس ترکیب منحصر به فرد مقادیر

1) از ابزار Feature to Raster برای تبدیل خاک به لایه رستری استفاده کنید (به کادر 3-6 مراجعه کنید). شطرنج خروجی را grounds_r نامگذاری کنید، Soil_order  را به عنوان فیلد برای تخصیص مقادیر به رستر خروجی انتخاب کنید و اندازه سلول خروجی را 40 (متر) تنظیم کنید.

2) ArcToolBox را باز کنید. به مسیر Spatial Analyst Tools > Local بروید و روی Combine دوبار کلیک کنید.

3) در کادر محاوره ای Combine :
الف) خاکهای_r و landcv_r را به عنوان رسترهای ورودی انتخاب کنید.
ب) رستری خروجی را نام ببرید.
ج) روی OK کلیک کنید. رستر خروجی ایجاد شده و به نمای داده اضافه می‌شود. هر مقدار سلول در رستری خروجی نشان دهنده یک مقدار منحصر به فرد است.

ترکیب خاک و پوشش زمین

4) روی رستر خروجی در فهرست مطالب کلیک راست کرده و Open Attribute Table را کلیک کنید. در جدول ویژگی، فیلدهای SOILS_R و LANDCV_R مقادیری را که نشان دهنده انواع خاک و پوشش زمین مربوط به هر مقدار سلول در رستری خروجی است فهرست می‌کنند.

 پوشش رستر از طریق عملیات ریاضی

5) با تخصیص مقدار 1 به نوع پوشش زمین جنگل‌های سوزنی برگ سرد معتدل و 0 به سایر انواع پوشش زمین با استفاده از ابزار Reclassify، landcv_r را مجدداً طبقه بندی کنید (به کادر 4-3 مراجعه کنید). نوع را به عنوان قسمت reclass  انتخاب کنید. رستر خروجی را سرد کانیفر نامگذاری کنید.
6) ArcToolBox را باز کنید. به Spatial Analyst Tools > Math بروید و روی Times دوبار کلیک کنید.
7) در کادر محاوره ای  :  Times
الف) grounds_r  را به عنوان رستر ورودی 1 و coldConifer را به عنوان شطرنج ورودی 2 انتخاب کنید.
ب) رستری خروجی را نام ببرید.
ج) روی OK کلیک کنید. رستری خروجی ایجاد می‌شود که انواع خاک (کد 1-5) را در منطقه جنگلی مخروطی سرد معتدل نشان می‌دهد. مناطق دیگر مقدار صفر را به دست می‌آورند.
کادر 4-11 نمونه‌هایی از توابع محلی در ArcGIS
کاربردی
برای پیروی از این مثال، ArcMap را شروع کنید و رستر dem، gradeSlope و landcv_r را از مسیر زیر دانلود کنید.
C:\Databases\GIS4EnvSci\VirtualCatchment\Geodata.gdb
توابع محلی در ArcGIS را می‌توان از مجموعه ابزار محلی Spatial Analyst یا ماشین حساب رستر دسترسی داشت. در این مثال از Raster Calculator استفاده شده است. پسوند Spatial Analyst را فعال کنید.
تبدیل شیب درجه به درصد شیب

1) ArcToolBox را باز کنید. به ابزارهای تحلیل مکانی > جبر نقشه بروید. روی Raster Calculator  دوبار کلیک کنید. شبیه شکل 4-19 به نظر می‌رسد.

شکل 4-19 ماشین حساب رستر

2) در کادر محاوره‌ای  Raster Calculator:
الف) فهرست توابع/ابزار را به پایین اسکرول کنید تا تابع Tan را پیدا کنید، سپس روی آن دوبار کلیک کنید. Tan() با قرار دادن مکان نما در داخل پرانتز به کادر عبارت اضافه می‌شود. این یک تابع مثلثاتی است که مماس هر مقدار سلول را در رستر ورودی محاسبه می‌کند. این نیاز به اندازه گیری مقادیر سلول ورودی بر حسب رادیان دارد.
ب) در کادر فهرست لایه‌ها و متغیرها، روی degreeSlope دوبار کلیک کنید. عبارت Tan (″degreeSlope″) می‌شود. با این حالdegreeSlope  شیب را بر حسب درجه ثبت می‌کند. مقادیر در درجه شیب ابتدا باید بر ضریب تبدیل رادیان به درجه π  تقسیم بر 180 یا تقریباً 296/57 تقسیم شوند.
ج) روی دکمه عملگر تقسیم  کلیک کنید، سپس 296/57 را وارد کنید. عبارت ( 296/57 Tan ″degreeSlope″/) می‌شود.
د) مکان نما را به انتهای عبارت منتقل کنید. روی دکمه عملگر ضرب (*) کلیک کنید و 100 را وارد کنید. عبارتTan(″degreeSlope″/57.296)*100   می‌شود. این عبارت جبر نقشه است که برای تبدیل درجه شیب به درصد شیب استفاده می‌شود.
ه) رستر خروجی را Slope نامگذاری کنید. روی OK کلیک کنید. رستر خروجی ایجاد می‌شود.
شناسایی علفزارهای واقع در نوع شیب خاصی در محدوده معینی از ارتفاع
3) ماشین حساب Raster را باز کنید.
4) در کادر محاوره‌ای  Raster Calculator :
الف) عبارت  Con(″degreeSlope″ < 15, 1, 0)  را وارد کنید.
ب)  Conتابع شرطی است که به کاربران اجازه می‌دهد تا مقدار خروجی هر سلول را بر اساس درست یا نادرست بودن مقدار سلول در یک عبارت شرطی مشخص شده در داخل پرانتز، کنترل کنند. در عبارت بالا، اگر مقدار یک سلول در درجه شیب کمتر از 15 (درجه) باشد، آنگاه 1 به محل سلول (true) در رستری خروجی اختصاص داده می‌شود. در غیر این صورت، مقادیر سلولی بزرگتر یا مساوی 15 در رستری خروجی صفر (نادرست) تخصیص داده می‌شود. با اجرای این تابع، در واقع یک عملیات طبقه بندی مجدد برای ایجاد یک نقشه باینری انجام می‌دهید.
ج) خروجی رستری را شیب ملایم نام ببرید. روی OK کلیک کنید. رستری باینری خروجی ایجاد می‌شود.
5) ماشین حساب Raster را دوباره باز کنید، سپس بیان جبری نقشه را وارد کنید
Con((dem > 1000) & (dem < 3000), 1, 0).
این عبارت 1 را به مکان‌های سلولی که ارتفاع آنها بیشتر از 1000 و کمتر از 3000 است و 0 را به مکان‌های سلول دیگر اختصاص می‌دهد. & عملگر AND منطقی است. رستری خروجی را Mid Elevation نامگذاری کنید. برای ایجاد رستر خروجی روی OK کلیک کنید.
6) رستر جدید به نام meadow با استفاده از ماشین حساب رستر با عبارت زیر ایجاد کنید.
Con(″landcv_r″ != 4, 0, 1)
در این عبارت != عملگر رابطه ای NOT EQUAL و 4 کد کلاس کلاس پوشش زمین Meadow است. اگر مقدار یک سلول در landcv_r برابر با 4 نباشد، 0 به محل سلول در رستری خروجی اختصاص داده می‌شود. در غیر این صورت، 1 به آن مکان سلول در رستری خروجی اختصاص داده می‌شود.
7) ماشین حساب Raster را دوباره باز کنید، سپس عبارت جبر نقشه را وارد کنید
″meadow″ * ″midElevation″ * ″gentleSlope″
رستری خروجی را نام ببرید. عبارت بالا را اجرا کنید. رستری خروجی مراتع را نشان می‌دهد که در ارتفاع 1000 تا 3000 متری و در شیب‌های کمتر از 15 درجه قرار دارند. این عملیات در اصل یک پوشش رستری است.

توابع کانونی یک مقدار جدید را برای هر سلول در شطرنج خروجی به عنوان تابعی از مقادیر، فواصل و/یا جهت سلول‌های موجود در همسایگی آن در رستری ورودی محاسبه می‌کنند. سلولی که در حال حاضر پردازش می‌شود، سلول کانونی نامیده می‌شود. محله را می‌توان در اشکال و اندازه‌های مختلف تعریف کرد. رایج‌ترین شکل‌های مورد استفاده عبارتند از: مستطیل (معمولاً با عرض و ارتفاع آن در سلول‌های متمرکز در سلول کانونی، به عنوان مثال 3×3)، دایره (تعریف شده با شعاع آن در مرکز سلول کانونی)، گوه (تعریف شعاع امتدادی از سلول کانونی) سلول کانونی، زوایای شروع و پایان که در خلاف جهت عقربه‌های ساعت از شرق اندازه‌گیری می‌شوند) و حلقۀ حلقوی (تعریف شده توسط شعاع‌های داخلی و خارجی در مرکز سلول کانونی) (شکل 4-20 را ببینید). محله همچنین می‌تواند به عنوان یک شکل نامنظم به نام کرنل تعریف شود که می‌تواند برای تعیین سلولهای محله در عملیات یا اعمال وزنهای مختلف به سلولهای محله استفاده شود ( شکل 4-21 را ببینید). توابع کانونی عمدتاً برای آمارهای کانونی یا همسایگی استفاده می‌شوند. مطالعه موردی 7 فصل 10 نمونه ای از نحوه استفاده از توابع کانونی در مدل سازی آلودگی هوا را ارائه می‌دهد. آنها همچنین برای فیلتر کردن در پردازش تصویر دیجیتال بسیار مفید هستند (به بخش 6-3 مراجعه کنید).

این توابع آمار خلاصه ای از مقادیر سلول را در همسایگی از پیش تعریف شده هر سلول محاسبه می‌کنند. در طول عملکرد تابع کانونی، سلول کانونی از یک سلول به سلول دیگر در جهت چپ به راست و از بالا به پایین منتقل می‌شود تا زمانی که همه سلول‌ها بازدید شوند. بنابراین، همسایگی سلول کانونی درست مانند یک پنجره متحرک است (شکل 4-22). توابع آماری فهرست شده در جدول 4-5 همگی در توابع کانونی استفاده می‌شوند. شکل 4-23 چند مثال ارائه می‌کند که در آنها از یک همسایگی مستطیلی 3×3 استفاده شده است، اما فقط مقادیر خروجی در سلول‌های سایه دار زرد نشان داده شده است. تابع کانونی «میانگین» میانگین مقادیر سلول‌ها را در همسایگی محاسبه می‌کند، به طور موثری تغییرات مکان را کاهش می‌دهد و داده‌ها را هموار می‌کند. اغلب به عنوان یک فیلتر پایین گذر برای صاف کردن یک تصویر دیجیتال برای کاهش نویزهای تصویر استفاده می‌شود. شکل 4-24 بخشی از مدل رقومی ارتفاع را قبل و بعد از هموارسازی با استفاده از میانگین نشان می‌دهد. تابع کانونی “محدوده” تفاوت بین حداکثر و حداقل مقدار سلول در همسایگی را محاسبه می‌کند. یک مقدار محدوده بزرگ مرز احتمالی بین کلاس‌ها یا دسته‌های مختلف را نشان می‌دهد. بنابراین تابع محدوده می‌تواند به عنوان یک فیلتر تشخیص لبه برای بهبود لبه در پردازش تصویر دیجیتال استفاده شود، به ویژه برای متمایز کردن ویژگی‌های خط و مرزهای ویژگی‌های منطقه مفید است. تابع کانونی “اکثریت” مقدار سلولی را که اغلب در همسایگی رخ می‌دهد تعیین می‌کند. بخش 6.3 بحث کاملی در مورد فیلترهای مکانی در پردازش تصویر دیجیتال ارائه می‌دهد. علاوه بر آمار کانونی، توابع شیب، جنبه و انحنای سطح (به بخش 7-4 مراجعه کنید) معمولاً در توابع کانونی گروه بندی می‌شوند، زیرا محاسبات آنها بر اساس مقادیر سلول‌های همسایگی است. کادر 4-12 نمونه‌هایی از استفاده از توابع کانونی در ArcGIS را ارائه می‌دهد.

شکل 4-20 انواع مجاورت : (الف) مستطیل، (ب) گوه، (ج) دایره و (د) حلقوی.

شکل 4-21 هسته ها

شکل 4-22 پنجره متحرک

شکل 4-23 نمونه‌هایی از آمار محله

شکل 4-24 بخشی از DEM : (الف) قبل و (ب) پس از صاف شدن

 

کادر 4-12 نمونه‌هایی از توابع کانونی در ArcGIS

کاربردی

برای پیروی از این مثال ArcMap  را شروع کنید و رستر dem و landcv_r را از مسیر زیر دانلود کنید.
C:\Databases\GIS4EnvSci\VirtualCatchment\Geodata.gdb
پسوند Spatial Analyst را فعال کنید. تنوع پوشش زمین در محله مشخص شده را محاسبه کنید
1) ArcToolBox را باز کنید. به ابزارهای تحلیل مکانی > محله بروید. روی Focal Statistics دوبار کلیک کنید.
2) در کادر محاوره ای Focal Statistics  :
الف) landcv_r  را به عنوان رستر ورودی انتخاب کنید.
ب) رستری خروجی را نام ببرید.
ج) Rectangle را به عنوان شکل همسایگی انتخاب کنید و ارتفاع و عرض آن را روی 5 خانه قرار دهید – یعنی اندازه محله 5×5 خانه است.
د) Variety را به عنوان نوع آمار انتخاب کنید، که تعداد مقادیر منحصربه‌فرد سلول‌های مجاورت را محاسبه می‌کند.
ه) روی OK کلیک کنید. رستری خروجی ایجاد می‌شود. رستری مناطق با بیش از یک نوع پوشش زمین را نشان می‌دهد. جدول ویژگی‌ها را باز کنید. این فهرست چند سلول وجود دارد که حاوی چندین نوع پوشش زمین هستند.

صاف کردن DEM با تابع کانونی MEAN

3) در ArcToolBox، به Spatial Analyst Tools > Neighborhood بروید. روی Focal Statistics دوبار کلیک کنید.
4) در کادر محاوره ای Focal Statistics  :
الف) dem  را به عنوان رستر ورودی انتخاب کنید.
ب) رستری خروجی را نام ببرید.
ج) Rectangle را به عنوان شکل محله انتخاب کنید و ارتفاع و عرض آن را روی 10 خانه قرار دهید. به طور کلی، هر چه اندازه محله بزرگتر باشد، نتیجه هموارتر است.
د) میانگین را به عنوان نوع آمار انتخاب کنید.
ه) روی OK کلیک کنید. رستری خروجی ایجاد می‌شود. بخشی از dem و خروجی صاف شده آن در شکل 4-24 نشان داده شده است.

توابع منطقه ای مقدار جدید را برای هر سلول در رستری خروجی به صورت a محاسبه می‌کنند. تابع مقادیر موجود سلول‌های متعلق به همان منطقه است. ناحیه گروهی از سلولها با ارزش یکسان است. این می‌تواند منطقه واحد یا چندین منطقه قطع شده که ارزش یکسانی دارند را تشکیل دهد. توابع منطقه ای آمار خلاصه تمام مقادیر سلول را در هر ناحیه محاسبه می‌کنند و آمار را در هر سلول منطقه در رستری خروجی می‌نویسند. آمار خلاصه شامل تمام موارد ذکر شده در جدول 4-5 است که مناطق معمولاً از لایه دیگری جدا از رستر مشخص می‌شوند که مقادیر سلول آنها باید به صورت آماری خلاصه شود. شکل 4-25 دو نمونه از آمار منطقه ای را نشان می‌دهد که میانگین و حداکثر ارتفاع را در هر نوع جنگل محاسبه می‌کند. برخی از بسته‌های GIS دارای توابع هندسی منطقه ای هستند که می‌توان از آنها برای محاسبه مساحت، محیط، ضخامت و سایر خصوصیات هندسی برای هر منطقه در یک رستر استفاده کرد. کادر 4-13 نحوه دسترسی به توابع منطقه ای در ArcGIS را نشان می‌دهد.

 

شکل 4-25 نمونه‌هایی از آمار منطقه ای

 

 

کادر 4-13 نمونه‌هایی از توابع ناحیه ای در ArcGIS

کاربردی

برای پیروی از این مثال، ArcMap  را شروع کنید و رستر dem و landcv_r را از مسیر زیر دانلود کنید.
C:\Databases\GIS4EnvSci\VirtualCatchment\Geodata.gdb
پسوند Spatial Analyst را فعال کنید.

محاسبه میانگین ارتفاع در هر نوع پوشش زمین

1)ArcToolBox  را باز کنید. به Spatial Analyst Tools > Zonal بروید. روی Zonal Statistics دوبار کلیک کنید.
2) در کادر محاوره ای  Zonal Statistics :
الف) landcv_r را به عنوان رستری ورودی و Value را به عنوان قسمت ناحیه انتخاب کنید. در اینجا از هر نوع پوشش زمین برای تعریف یک منطقه استفاده می‌شود.
ب) dem را به عنوان رستری مقدار ورودی انتخاب کنید.
ج) رستری خروجی را نام ببرید.
د) میانگین را به عنوان نوع آمار انتخاب کنید.
ه) روی OK کلیک کنید. رستری خروجی ایجاد می‌شود. میانگین ارتفاع برای هر منطقه پوشش زمین را نشان می‌دهد.
3) در ArcToolBox، به Spatial Analyst Tools > Zonal بروید.
روی Zonal Statistics as Table دوبار کلیک کنید.
4)  در کادر محاوره ای  Zonal Statistics as Table:
الف) landcv_r  را به عنوان رستری ورودی و Type را به عنوان قسمت zone انتخاب کنید.
ب) dem را به عنوان رستری مقدار ورودی انتخاب کنید.
ج) جدول خروجی را نام ببرید.
د) ALL را به عنوان نوع آمار انتخاب کنید و روی OK کلیک کنید. جدولی ایجاد می‌شود که شامل هفت آمار خلاصه ارتفاع در هر نوع پوشش زمین شامل حداقل، حداکثر، محدوده، میانگین، انحراف معیار و مجموع است.

محاسبه هندسه پوشش زمین  

5) در ArcToolBox، به Spatial Analyst Tools > Zonal بروید. روی Zonal Geometry دوبار کلیک کنید.
6)در کادر محاوره ای Zonal Statistics:
الف) landcv_r را به عنوان رستری ورودی و Value را به عنوان قسمت ناحیه انتخاب کنید.
ب) رستری خروجی را نام ببرید.
ج) Area را به عنوان نوع هندسه و 40 را به عنوان اندازه سلول خروجی انتخاب کنید.
د) روی OK کلیک کنید. رستری خروجی ایجاد می‌شود. به هر سلول در رستری خروجی، یک مقدار مساحت برای نوع پوشش زمینی که به آن تعلق دارد، اختصاص داده می‌شود.
7) در ArcToolBox، به Spatial Analyst Tools > Zonal بروید. روی Zonal Geometry as Table دوبار کلیک کنید.
8) در کادر محاوره ای  Zonal Geometry as Table:
الف) landcv_r را به عنوان رستری ورودی و Value را به عنوان قسمت ناحیه انتخاب کنید.
ب) جدول خروجی را نام ببرید.
ج) اندازه سلول پردازش را 40 تنظیم کنید.
د) روی OK کلیک کنید. جدول خروجی تولید می‌شود که تمام ویژگی‌های هندسی اصلی هشت نوع پوشش زمین شامل مساحت، محیط، ضخامت و ویژگی‌های بیضی را فهرست می‌کند.

توابع جهانی مقدار جدیدی را برای هر سلول در رستر خروجی به عنوان تابعی از مقادیر موجود همه سلول‌ها در رستر ورودی محاسبه می‌کنند. توابع سطح مسافت و کم هزینه ، که قبلاً مورد بحث قرار گرفتند و توابع تجزیه و تحلیل دید (بخش 7.5 را ببینید) توابع جهانی هستند.

در ArcGIS همه ابزارها، اپراتورها و توابع Spatial Analyst را می‌توان با جبر نقشه استفاده کرد. علاوه بر توابع محلی، کانونی، منطقه ای و جهانی که در بالا مورد بحث قرار گرفت، جبر نقشه در ArcGIS می‌تواند از توابع خاص کاربردی در تحلیل مکانی ، از جمله هیدرولوژی، آب‌های زیرزمینی، تحلیل سطح و توابع تابش خورشیدی استفاده و پیاده سازی کند. دو روش برای استفاده از جبر نقشه در ArcGIS وجود دارد. یکی استفاده از ابزار Raster Calculator، همانطور که در کادر 4-11 نشان داده شده است. مورد دوم استفاده از بسته سایت ArcPy با پردازش جغرافیایی است که محیط یکپارچه پایتون را برای کاربران فراهم می‌کند تا بتوانند عبارات جبری نقشه را با اسکریپت‌های پایتون بنویسند و از ماژول‌ها و کتابخانه‌های پایتون و شخص ثالث پایتون استفاده کنند. استفاده از محیط یکپارچه پایتون ArcGIS مستلزم آشنایی با پایتون است که از حوصله این کتاب خارج است.

برگرفته از کتاب کاربرد GISدر محیط زیست

ترجمه:سعید جوی زاده،شهناز تیموری،فاطمه حسین پور فرزانه

درخواست مشاوره رایگان و شرکت در دوره

برخی از تالیفات دکتر سعید جوی زاده

دوره های پربازید و کاربردی

بدون دیدگاه

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

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