GDALDEM(1) GDAL GDALDEM(1)

gdaldem - ابزارهایی برای تحلیل و مصورسازی مدل‌های رقومی ارتفاع (DEM)

gdaldem [--help] [--help-general] <mode> <input> <output> <options>

از هر رستر ارتفاعی پشتیبانی‌شده در GDAL:

تولید یک نقشه سایه‌روشن برجستگی (shaded relief):

gdaldem hillshade <input_dem> <output_hillshade>
            [-z <zfactor>] [[-s <scale>] | [-xscale <xscale> -yscale <yscale>]]
            [-az <azimuth>] [-alt <altitude>]
            [-alg ZevenbergenThorne] [-combined | -multidirectional | -igor]
            [-compute_edges] [-b <Band>] [-of <format>] [-co <NAME>=<VALUE>]... [-q]

تولید یک نقشه شیب (slope map):

gdaldem slope <input_dem> <output_slope_map>
            [-p] [[-s <scale>] | [-xscale <xscale> -yscale <yscale>]]
            [-alg ZevenbergenThorne]
            [-compute_edges] [-b <band>] [-of <format>] [-co <NAME>=<VALUE>]... [-q]

تولید یک نقشه جهت شیب (aspect map)، یک رستر اعشاری ۳۲ بیتی با مقادیر پیکسل بین 0 تا 360 که نشان‌دهنده آزیموت است خروجی می‌دهد:

gdaldem aspect <input_dem> <output_aspect_map>
            [-trigonometric] [-zero_for_flat]
            [-alg ZevenbergenThorne]
            [-compute_edges] [-b <band>] [-of format] [-co <NAME>=<VALUE>]... [-q]

تولید یک نقشه رنگی برجستگی (color relief):

gdaldem color-relief <input_dem> <color_text_file> <output_color_relief_map>
             [-alpha] [-exact_color_entry | -nearest_color_entry]
             [-b <band>] [-of format] [-co <NAME>=<VALUE>]... [-q]
که در آن color_text_file شامل خطوطی با قالب "elevation_value red green blue [alpha]" است. در صورت وجود ستون alpha می‌توان آن را با '-alpha' فعال کرد.

تولید یک نقشه شاخص ناهمواری زمین (TRI):

gdaldem TRI input_dem output_TRI_map
            [-alg Wilson|Riley]
            [-compute_edges] [-b Band (default=1)] [-of format] [-q]

تولید یک نقشه شاخص موقعیت توپوگرافی (TPI):

gdaldem TPI <input_dem> <output_TPI_map>
            [-compute_edges] [-b <band>] [-of <format>] [-co <NAME>=<VALUE>]... [-q]

تولید یک نقشه زبری سطح (roughness map):

gdaldem roughness <input_dem> <output_roughness_map>
            [-compute_edges] [-b <band>] [-of <format>] [-co <NAME>=<VALUE>]... [-q]

دستور gdaldem با استفاده از الگوریتم‌های مختلف برای تحلیل عوارض زمین، رسترهای خروجی تولید می‌کند.

به طور کلی، این ابزار فرض می‌کند که واحدهای x و y و z یکسان هستند. با این حال، برای محاسبه سایه‌روشن (hillshade) و شیب (slope)، از نگارش GDAL 3.11 به بعد، در صورتی که هیچ‌یک از گزینه‌های -scale، -xscale و -yscale مشخص نشده باشند و CRS یک CRS جغرافیایی یا تصویرشده باشد، gdaldem به طور خودکار نسبت مناسب را بر اساس واحدهای CRS و همچنین مقدار احتمالی واحدهای باند رستر (که توسط GDALRasterBand::GetUnitType() برگردانده می‌شود، اگر متر، فوت بین‌المللی یا فوت نقشه‌برداری آمریکا باشد) تعیین می‌کند. توجه داشته باشید که برای CRS جغرافیایی، نتیجه برای مجموعه داده‌های ورودی در عرض‌های جغرافیایی بالا ممکن است نادرست باشد و ممکن است بازتصویرسازی پیشین به یک تصویر قطبی لازم باشد.

برای محاسبه سایه‌روشن و شیب، اگر واحدهای x (شرق-غرب) و y (شمال-جنوب) یکسان باشند اما واحدهای z (ارتفاع) متفاوت باشند، می‌توان از گزینه -scale (یا -xscale و -yscale) برای تنظیم نسبت واحدهای عمودی به افقی استفاده کرد. برای سیستم‌های مختصات طول و عرض جغرافیایی (LatLong) در نزدیکی خط استوا، که در آن‌ها واحدهای عرض جغرافیایی و طول جغرافیایی مشابه هستند، واحدهای ارتفاع (z) را می‌توان با استفاده از scale=370400 (اگر ارتفاع بر حسب فوت باشد) یا scale=111120 (اگر ارتفاع بر حسب متر باشد) برای سازگاری تبدیل کرد. برای مکان‌هایی که نزدیک خط استوا نیستند، بهتر است پیش از استفاده از gdaldem، شبکه رستر خود را با استفاده از gdalwarp بازتصویرسازی کنید.

<mode>
که در آن <mode> یکی از هفت حالت موجود است:
•
hillshade
تولید یک نقشه سایه‌روشن برجستگی از هر رستر ارتفاعی پشتیبانی‌شده توسط GDAL.
•
slope
تولید یک نقشه شیب از هر رستر ارتفاعی پشتیبانی‌شده توسط GDAL.
•
aspect
تولید یک نقشه جهت شیب از هر رستر ارتفاعی پشتیبانی‌شده توسط GDAL.
•
color-relief
تولید یک نقشه رنگی برجستگی از هر رستر ارتفاعی پشتیبانی‌شده توسط GDAL.
•
TRI
تولید یک نقشه شاخص ناهمواری زمین (Terrain Ruggedness Index) از هر رستر ارتفاعی پشتیبانی‌شده توسط GDAL.
•
TPI
تولید یک نقشه شاخص موقعیت توپوگرافی (Topographic Position Index) از هر رستر ارتفاعی پشتیبانی‌شده توسط GDAL.
•
roughness
تولید یک نقشه زبری سطح از هر رستر ارتفاعی پشتیبانی‌شده توسط GDAL.

نکته:

معادل‌ها در واسط خط فرمان جدید "gdal":
  • gdal raster aspect
  • gdal raster color-map
  • gdal raster hillshade
  • gdal raster slope
  • gdal raster tpi
  • gdal raster tri

گزینه‌های عمومی زیر در دسترس هستند:

نمایش این پیام راهنما و خروج.
--help-general
نمایش یک پیام کاربرد مختصر برای گزینه‌های عمومی خط فرمان GDAL و خروج.
<input_dem>
رستر DEM ورودی که باید پردازش شود.
<output_xxx_map>
رستر خروجی تولیدشده.
انتخاب قالب خروجی.
در صورت عدم تعیین، قالب از پسوند حدس زده می‌شود (قبلاً GTiff -- GeoTIFF File Format بود). از نام کوتاه قالب استفاده کنید.
انجام محاسبات در لبه‌های رستر و نزدیک مقادیر nodata.
انتخاب یک باند ورودی برای پردازش. شماره‌گذاری باندها از 1 شروع می‌شود.
بسیاری از قالب‌ها یک یا چند گزینه اختیاری ایجاد دارند که می‌توان از آن‌ها برای کنترل ویژگی‌های خاص فایل ایجادشده استفاده کرد. برای نمونه، درایور GeoTIFF از گزینه‌های ایجاد برای کنترل فشرده‌سازی و تعیین اینکه آیا فایل باید تایل‌بندی شود یا خیر پشتیبانی می‌کند.

گزینه‌های ایجاد موجود بسته به درایور قالب متفاوت است و برخی قالب‌های ساده هیچ گزینه ایجادی ندارند. فهرستی از گزینه‌های پشتیبانی‌شده برای یک قالب را می‌توان با گزینه خط فرمان --format مشاهده کرد، اما مستندات مربوط به همان قالب مرجع قطعی اطلاعات درباره گزینه‌های ایجاد درایور است. برای گزینه‌های مجاز ایجاد هر قالب، به مستندات اختصاصی قالب در Raster drivers مراجعه کنید.

حالت بی‌صدا؛ جلوگیری از نمایش پیام‌های پیشرفت و سایر خروجی‌های غیرخطا.

برای تمام الگوریتم‌ها، به جز color-relief، اگر حداقل یک پیکسل با مقدار nodata در پنجره 3x3 متمرکز بر هر پیکسل مبدأ یافت شود، یک مقدار nodata در مجموعه داده هدف درج خواهد شد. پیامد این امر این است که یک حاشیه 1 پیکسلی در اطراف هر تصویر با مقدار nodata تنظیم می‌شود.

اگر -compute_edges مشخص شده باشد، gdaldem مقادیر را در لبه‌های تصویر یا در صورت یافتن مقدار nodata در پنجره 3x3 با میان‌یابی مقادیر گمشده محاسبه می‌کند.

این دستور یک رستر 8 بیتی با جلوه سایه‌روشن برجستگی زیبا تولید می‌کند. این ابزار برای تجسم عوارض زمین بسیار مفید است. شما می‌توانید به صورت اختیاری آزیموت و زاویه ارتفاع منبع نور، ضریب اغراق عمودی و ضریب مقیاس‌بندی را برای در نظر گرفتن تفاوت بین واحدهای عمودی و افقی مشخص کنید.

مقدار 0 به عنوان مقدار nodata در خروجی استفاده می‌شود.

گزینه‌های اختصاصی زیر در دسترس هستند:

منابع علمی نشان می‌دهند که فرمول Zevenbergen & Thorne برای مناظر هموارتر مناسب‌تر است، در حالی که فرمول Horn روی عوارض ناهموارتر عملکرد بهتری دارد.
ضریب اغراق عمودی که برای پیش‌ضرب در ارتفاعات استفاده می‌شود.
نسبت واحدهای عمودی به واحدهای افقی. اگر واحد افقی DEM مبدأ درجه باشد (مانند تصویر Lat/Long WGS84)، در صورتی که واحدهای عمودی متر باشند می‌توانید از scale=111120 استفاده کنید (یا اگر بر حسب فوت باشند از scale=370400).

از نگارش GDAL 3.11 به بعد، اگر هیچ‌یک از گزینه‌های -scale، -xscale و -yscale مشخص نشده باشند و CRS از نوع جغرافیایی یا تصویرشده باشد، gdaldem به طور خودکار نسبت مناسب را بر اساس واحدهای CRS و همچنین مقدار احتمالی واحدهای باند رستر (که توسط GDALRasterBand::GetUnitType() برگردانده می‌شود، اگر متر، فوت بین‌المللی یا فوت نقشه‌برداری آمریکا باشد) تعیین می‌کند. توجه داشته باشید که برای CRS جغرافیایی، نتیجه برای مجموعه داده‌های مبدأ در عرض‌های جغرافیایی بالا ممکن است نادرست باشد و ممکن است بازتصویرسازی پیشین به یک تصویر قطبی لازم باشد.

اثر این گزینه مشابه مشخص کردن -xscale و -yscale با همان مقدار -scale است. گزینه -scale با گزینه‌های -xscale و -yscale ناسازگار است.

افزوده‌شده در نگارش 3.11.

نسبت واحدهای عمودی به واحدهای افقی محور X. اگر واحد افقی DEM مبدأ درجه باشد (مانند تصویر Lat/Long WGS84)، در صورتی که واحدهای عمودی متر باشند می‌توانید از scale=111120 استفاده کنید (یا اگر بر حسب فوت باشند از scale=370400).

از نگارش GDAL 3.11 به بعد، اگر هیچ‌یک از گزینه‌های -scale، -xscale و -yscale مشخص نشده باشند و CRS از نوع جغرافیایی یا تصویرشده باشد، gdaldem به طور خودکار نسبت مناسب را بر اساس واحدهای CRS و همچنین مقدار احتمالی واحدهای باند رستر (که توسط GDALRasterBand::GetUnitType() برگردانده می‌شود، اگر متر، فوت بین‌المللی یا فوت نقشه‌برداری آمریکا باشد) تعیین می‌کند. توجه داشته باشید که برای CRS جغرافیایی، نتیجه برای مجموعه داده‌های مبدأ در عرض‌های جغرافیایی بالا ممکن است نادرست باشد و ممکن است بازتصویرسازی پیشین به یک تصویر قطبی لازم باشد.

اگر -xscale مشخص شود، -yscale نیز باید مشخص گردد.

افزوده‌شده در نگارش 3.11.

نسبت واحدهای عمودی به واحدهای افقی محور Y. اگر واحد افقی DEM مبدأ درجه باشد (مانند تصویر Lat/Long WGS84)، در صورتی که واحدهای عمودی متر باشند می‌توانید از scale=111120 استفاده کنید (یا اگر بر حسب فوت باشند از scale=370400).

از نگارش GDAL 3.11 به بعد، اگر هیچ‌یک از گزینه‌های -scale، -xscale و -yscale مشخص نشده باشند و CRS از نوع جغرافیایی یا تصویرشده باشد، gdaldem به طور خودکار نسبت مناسب را بر اساس واحدهای CRS و همچنین مقدار احتمالی واحدهای باند رستر (که توسط GDALRasterBand::GetUnitType() برگردانده می‌شود، اگر متر، فوت بین‌المللی یا فوت نقشه‌برداری آمریکا باشد) تعیین می‌کند. توجه داشته باشید که برای CRS جغرافیایی، نتیجه برای مجموعه داده‌های مبدأ در عرض‌های جغرافیایی بالا ممکن است نادرست باشد و ممکن است بازتصویرسازی پیشین به یک تصویر قطبی لازم باشد.

اگر -yscale مشخص شود، -xscale نیز باید مشخص گردد.

زاویه آزیموت (سمت) تابش نور بر حسب درجه. مقدار 0 بدین معنی است که نور از بالای رستر می‌تابد، 90 از شرق و... مقدار پیش‌فرض، 315، به ندرت باید تغییر داده شود زیرا مقداری است که معمولاً برای تولید نقشه‌های سایه‌روشن استفاده می‌شود.
زاویه ارتفاع تابش نور بر حسب درجه. مقدار 90 بدین معنی است که نور از بالای سر DEM می‌تابد و 0 یعنی نور مماس (افقی) می‌تابد.
سایه‌زنی ترکیبی، ترکیبی از شیب و سایه‌زنی مورب (oblique shading).
سایه‌زنی چندجهته، ترکیبی از سایه‌روشن‌های تابیده‌شده از آزیموت‌های 225 درجه، 270 درجه، 315 درجه و 360 درجه.
سایه‌زنی که تلاش می‌کند اثرات منفی روی سایر عوارض نقشه در زیر لایه را به حداقل برساند. با گزینه -alt قابل استفاده نیست.

افزوده‌شده در نگارش 3.0.

سایه‌روشن چندجهته از فرمول ارائه‌شده در http://pubs.usgs.gov/of/1992/of92-422/of92-422.pdf استفاده می‌کند.

سایه‌روشن Igor از فرمول Maperitive http://maperitive.net/docs/Commands/GenerateReliefImageIgor.html استفاده می‌کند.

این دستور یک رستر DEM را دریافت کرده و یک رستر اعشاری ۳۲ بیتی با مقادیر شیب خروجی می‌دهد. شما این امکان را دارید که نوع مقدار شیب مورد نظر خود را مشخص کنید: درجه یا درصد شیب. در مواردی که واحدهای افقی با واحدهای عمودی متفاوت هستند، می‌توانید ضریب مقیاس‌بندی را نیز مشخص کنید.

مقدار -9999 به عنوان مقدار nodata خروجی استفاده می‌شود.

گزینه‌های اختصاصی زیر در دسترس هستند:

منابع علمی نشان می‌دهند که فرمول Zevenbergen & Thorne برای مناظر هموارتر مناسب‌تر است، در حالی که فرمول Horn روی عوارض ناهموارتر عملکرد بهتری دارد.
در صورت مشخص شدن، شیب بر حسب درصد بیان می‌شود. در غیر این صورت، بر حسب درجه محاسبه می‌گردد.
نسبت واحدهای عمودی به واحدهای افقی. اگر واحد افقی DEM مبدأ درجه باشد (مانند تصویر Lat/Long WGS84)، در صورتی که واحدهای عمودی متر باشند می‌توانید از scale=111120 استفاده کنید (یا اگر بر حسب فوت باشند از scale=370400).

از نگارش GDAL 3.11 به بعد، اگر هیچ‌یک از گزینه‌های -scale، -xscale و -yscale مشخص نشده باشند و CRS از نوع جغرافیایی یا تصویرشده باشد، gdaldem به طور خودکار نسبت مناسب را بر اساس واحدهای CRS و همچنین مقدار احتمالی واحدهای باند رستر (که توسط GDALRasterBand::GetUnitType() برگردانده می‌شود، اگر متر، فوت بین‌المللی یا فوت نقشه‌برداری آمریکا باشد) تعیین می‌کند. توجه داشته باشید که برای CRS جغرافیایی، نتیجه برای مجموعه داده‌های مبدأ در عرض‌های جغرافیایی بالا ممکن است نادرست باشد و ممکن است بازتصویرسازی پیشین به یک تصویر قطبی لازم باشد.

اثر این گزینه مشابه مشخص کردن -xscale و -yscale با همان مقدار -scale است. گزینه -scale با گزینه‌های -xscale و -yscale ناسازگار است.

افزوده‌شده در نگارش 3.11.

نسبت واحدهای عمودی به واحدهای افقی محور X. اگر واحد افقی DEM مبدأ درجه باشد (مانند تصویر Lat/Long WGS84)، در صورتی که واحدهای عمودی متر باشند می‌توانید از scale=111120 استفاده کنید (یا اگر بر حسب فوت باشند از scale=370400).

از نگارش GDAL 3.11 به بعد، اگر هیچ‌یک از گزینه‌های -scale، -xscale و -yscale مشخص نشده باشند و CRS از نوع جغرافیایی یا تصویرشده باشد، gdaldem به طور خودکار نسبت مناسب را بر اساس واحدهای CRS و همچنین مقدار احتمالی واحدهای باند رستر (که توسط GDALRasterBand::GetUnitType() برگردانده می‌شود، اگر متر، فوت بین‌المللی یا فوت نقشه‌برداری آمریکا باشد) تعیین می‌کند. توجه داشته باشید که برای CRS جغرافیایی، نتیجه برای مجموعه داده‌های مبدأ در عرض‌های جغرافیایی بالا ممکن است نادرست باشد و ممکن است بازتصویرسازی پیشین به یک تصویر قطبی لازم باشد.

اگر -xscale مشخص شود، -yscale نیز باید مشخص گردد.

افزوده‌شده در نگارش 3.11.

نسبت واحدهای عمودی به واحدهای افقی محور Y. اگر واحد افقی DEM مبدأ درجه باشد (مانند تصویر Lat/Long WGS84)، در صورتی که واحدهای عمودی متر باشند می‌توانید از scale=111120 استفاده کنید (یا اگر بر حسب فوت باشند از scale=370400).

از نگارش GDAL 3.11 به بعد، اگر هیچ‌یک از گزینه‌های -scale، -xscale و -yscale مشخص نشده باشند و CRS از نوع جغرافیایی یا تصویرشده باشد، gdaldem به طور خودکار نسبت مناسب را بر اساس واحدهای CRS و همچنین مقدار احتمالی واحدهای باند رستر (که توسط GDALRasterBand::GetUnitType() برگردانده می‌شود، اگر متر، فوت بین‌المللی یا فوت نقشه‌برداری آمریکا باشد) تعیین می‌کند. توجه داشته باشید که برای CRS جغرافیایی، نتیجه برای مجموعه داده‌های مبدأ در عرض‌های جغرافیایی بالا ممکن است نادرست باشد و ممکن است بازتصویرسازی پیشین به یک تصویر قطبی لازم باشد.

اگر -yscale مشخص شود، -xscale نیز باید مشخص گردد.

این دستور یک رستر اعشاری ۳۲ بیتی با مقادیری بین 0° تا 360° خروجی می‌دهد که نشان‌دهنده جهت (آزیموت) شیب‌ها است. تعریف آزیموت به این صورت است: 0° بدین معنی است که شیب رو به شمال است، 90° رو به شرق، 180° رو به جنوب و 270° رو به غرب (به شرطی که بالای رستر ورودی شما رو به شمال باشد). مقدار جهت شیب -9999 به عنوان مقدار nodata برای نشان دادن جهت شیب تعریف‌نشده در مناطق مسطح با شیب صفر استفاده می‌شود.

گزینه‌های اختصاصی زیر در دسترس هستند:

منابع علمی نشان می‌دهند که فرمول Zevenbergen & Thorne برای مناظر هموارتر مناسب‌تر است، در حالی که فرمول Horn روی عوارض ناهموارتر عملکرد بهتری دارد.
بازگرداندن زاویه مثلثاتی به جای آزیموت. بنابراین 0° یعنی شرق، 90° شمال، 180° غرب و 270° جنوب.
بازگرداندن مقدار 0 برای مناطق مسطح با شیب صفر، به جای -9999.

با استفاده از این 2 گزینه، جهت شیب بازگردانده‌شده توسط gdaldem aspect با مقدار r.slope.aspect در سامانه GRASS یکسان خواهد بود. در غیر این صورت، با ابزار aspect.cpp اثر Matthew Perry یکسان است.

این دستور یک رستر ۳ بانده (RGB) یا ۴ بانده (RGBA) با مقادیری که از ارتفاع و یک فایل پیکربندی متنی رنگ محاسبه می‌شوند خروجی می‌دهد؛ این فایل حاوی نگاشت بین مقادیر مختلف ارتفاع و رنگ مورد نظر متناظر است. به طور پیش‌فرض، رنگ‌های بین مقادیر ارتفاع داده‌شده به آرامی ترکیب می‌شوند و نتیجه یک DEM رنگ‌آمیزی‌شده زیبا است. گزینه‌های -exact_color_entry یا -nearest_color_entry را می‌توان برای جلوگیری از میان‌یابی خطی برای مقادیری که با هیچ نمایه‌ای از فایل پیکربندی رنگ مطابقت ندارند استفاده کرد.

گزینه‌های اختصاصی زیر در دسترس هستند:

فایل متنی پیکربندی رنگ.
افزودن یک کانال آلفا به رستر خروجی.
استفاده از تطابق دقیق هنگام جستجو در فایل پیکربندی رنگ. اگر هیچ ورودی رنگ منطبقی یافت نشود، از مقدار چهارتایی RGBA به صورت "0,0,0,0" استفاده خواهد شد.
استفاده از مقدار چهارتایی RGBA متناظر با نزدیک‌ترین ورودی در فایل پیکربندی رنگ.

حالت color-relief تنها حالتی است که از VRT به عنوان قالب خروجی پشتیبانی می‌کند. در این حالت، فایل پیکربندی رنگ به عناصر مناسب جدول تطبیق رنگ (LUT) تبدیل می‌شود. توجه داشته باشید که ارتفاع‌های مشخص‌شده به صورت درصد به مقادیر مطلق ترجمه می‌شوند، که در صورت تفاوت آماری رستر مبدأ با آنچه در ساخت VRT استفاده شده است، باید مد نظر قرار گیرد.

فایل متنی پیکربندی رنگ عموماً شامل ۴ ستون در هر خط است: مقدار ارتفاع و مؤلفه‌های متناظر قرمز، سبز و آبی (بین 0 تا 255). مقدار ارتفاع می‌تواند هر عدد اعشاری، یا کلمه کلیدی nv برای مقدار nodata باشد. ارتفاع همچنین می‌تواند به صورت درصد بیان شود: 0% کمترین مقدار موجود در رستر و 100% بیشترین مقدار است.

جداکننده‌های مختلف فیلد پذیرفته می‌شوند: کاما، تب (tab)، فاصله و ':'.

رنگ‌های متداول مورد استفاده در GRASS را می‌توان به جای سه‌تایی RGB، با نام آن‌ها نیز مشخص کرد. فهرست پشتیبانی‌شده عبارت است از: white, black, red, green, blue, yellow, magenta, cyan, aqua, grey/gray, orange, brown, purple/violet و indigo.

فایل‌های پالت GMT .cpt نیز پشتیبانی می‌شوند (فقط COLOR_MODEL = RGB).

نکته: ساختار دستوری فایل پیکربندی رنگ برگرفته از ساختار پشتیبانی‌شده توسط ابزار r.colors در GRASS است. فایل‌های جدول رنگ HDR شرکت ESRI با پسوند (.clr) نیز با این ساختار مطابقت دارند. مؤلفه آلفا و پشتیبانی از tab و کاما به عنوان جداکننده، گسترش‌های اختصاصی GDAL هستند.

برای مثال:

3500   white
2500   235:220:175
50%   190 185 135
700    240 250 150
0      50  180  50
nv     0   0   0   0

برای پیاده‌سازی یک حالت گرد کردن به پایین، مقدار ارتفاع را می‌توان با مقداری جدید که کمی بالاتر از آستانه است تکرار کرد. برای نمونه برای داشتن قرمز در [0,10]، سبز در ]10,20] و آبی در ]20,30]:

0       red
10      red
10.001  green
20      green
20.001  blue
30      blue

این دستور یک رستر تک‌بانده با مقادیر محاسبه‌شده از ارتفاع خروجی می‌دهد. TRI مخفف شاخص ناهمواری زمین (Terrain Ruggedness Index) است که تفاوت بین یک پیکسل مرکزی و سلول‌های اطراف آن را اندازه‌گیری می‌کند.

مقدار -9999 به عنوان مقدار nodata خروجی استفاده می‌شود.

گزینه زیر در دسترس است:

از نگارش GDAL 3.3 به بعد، الگوریتم Riley (رجوع کنید به Riley, S.J., De Gloria, S.D., Elliot, R. (1999): A Terrain Ruggedness that Quantifies Topographic Heterogeneity. Intermountain Journal of Science, Vol.5, No.1-4, pp.23-27) در دسترس بوده و مقدار پیش‌فرض جدید است. این الگوریتم از جذر مجموع مربعات اختلاف بین یک پیکسل مرکزی و سلول‌های اطراف آن استفاده می‌کند. این گزینه برای سناریوهای خشکی پیشنهاد می‌شود.

الگوریتم Wilson (رجوع کنید به Wilson et al 2007, Marine Geodesy 30:3-35) از میانگین اختلاف بین یک پیکسل مرکزی و سلول‌های پیرامون آن استفاده می‌کند. این گزینه برای سناریوهای ژرفاسنجی پیشنهاد می‌شود.

این دستور یک رستر تک‌بانده با مقادیر محاسبه‌شده از ارتفاع خروجی می‌دهد. TPI مخفف شاخص موقعیت توپوگرافی (Topographic Position Index) است که به عنوان تفاوت بین یک پیکسل مرکزی و میانگین سلول‌های پیرامون آن تعریف می‌شود (رجوع کنید به Wilson et al 2007, Marine Geodesy 30:3-35).

مقدار -9999 به عنوان مقدار nodata خروجی استفاده می‌شود.

هیچ گزینه اختصاصی ندارد.

این دستور یک رستر تک‌بانده با مقادیر محاسبه‌شده از ارتفاع خروجی می‌دهد. زبری سطح (Roughness) بزرگ‌ترین اختلاف بین‌سلولی یک پیکسل مرکزی و سلول‌های پیرامون آن است، طبق تعریف Wilson et al (2007, Marine Geodesy 30:3-35).

مقدار -9999 به عنوان مقدار nodata خروجی استفاده می‌شود.

هیچ گزینه اختصاصی ندارد.

برنامه در صورت موفقیت کد وضعیت 0 و در صورت بروز خطا مقدار غیرصفر برمی‌گرداند (خطاهای غیرمسدودکننده که به صورت هشدار صادر می‌شوند، به عنوان اجرای موفقیت‌آمیز در نظر گرفته می‌شوند).

این ابزار از طریق C نیز با تابع GDALDEMProcessing() قابل فراخوانی است.

Matthew Perry perrygeo@gmail.com, Even Rouault even.rouault@spatialys.com, Howard Butler hobu.inc@gmail.com, Chris Yesson chris.yesson@ioz.ac.uk

برگرفته از کدهای نوشته‌شده توسط Michael Shapiro, Olga Waupotitsch, Marjorie Larson, Jim Westervelt: U.S. Army CERL, 1993. GRASS 4.1 Reference Manual. U.S. Army Corps of Engineers, Construction Engineering Research Laboratories, Champaign, Illinois, 1-425.

مستندات ابزارهای مرتبط در GRASS:

https://grass.osgeo.org/grass84/manuals/r.slope.aspect.html

https://grass.osgeo.org/grass84/manuals/r.relief.html

https://grass.osgeo.org/grass84/manuals/r.colors.html

Matthew Perry <perrygeo@gmail.com>, Even Rouault <even.rouault@spatialys.com>, Howard Butler <hobu.inc@gmail.com>, Chris Yesson <chris.yesson@ioz.ac.uk>

1998-2026

August 18, 2026