Розуміння значень мінімуму та максимуму за допомогою Mosaic до New Raster?

Мені доводиться поєднувати кілька растрових DEM в одну з цього джерела: http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp , тому я використовую Mosaic to New Raster tool.

Я налаштував інструмент:

Я ввожу растри (усі з одного джерела, однакового розміру і без проекції);
Налаштувати тип пікселя (32-бітна плаваюча точка як у оригінальних растрах);
Встановити розмір комірок як оригінальні растри;
Кількість смуг = 1 як вихідні растри; Мозаїчний оператор, я це зробив за допомогою BLEND та MEAN (я отримую той самий результат).

Проблема, яку я маю, полягає в тому, що отримана мозаїка показує інший діапазон максимальних і мінімальних значень, ніж максимальне і мінімальне значення для окремих растрових, наприклад, растр 1 (-5123,8, 23,25), растр 2 (-5974,6, 40,09), растр 3 (-57770,2, 38), растр 4 (-2534,3, 23,55) та кінцевий мозаїчний растр (-5975,8, 81,1).

Я думаю, що це рішення не є правильним, принаймні, я не сподівався цього отримати. У когось є ідея, якщо це нормально, а якщо це не так, як це вирішити та отримати правильний мозаїчний растр із правильними максимальними та мінімальними значеннями?

Я використовую ArcGIS 10.2.2 для робочого столу.

— Том
джерело

Не могли б ви сказати нам, як ви знайшли статистику для вхідних та вихідних растрах? Іноді, в інтересах ефективності, ArcGIS обманює і використовує лише підмножину даних для оцінки статистики. Зокрема, це означає, що це може бути відхиленням у звітності про мінімуми та максимуми.

— whuber

Дякуємо за ваш інтерес Я клацаю правою кнопкою миші на шарі (як початкові растри, так і мозаїка) і заглядаю у Властивості. Варто зазначити, що я перевірив переключення значень кожного растру в Symbology з Streched на Unique, і виявляється, що якщо я переключусь на значення Streched для всіх шарів, то воно відповідає мінімальному значенню Streched з одного шару і макс. підкреслене значення від іншого, до значень min та max, які я бачу при створенні мозаїки. Хтось міг би мені пояснити, чи все-таки правильна процедура, і ці відмінності між значеннями Streched та Uniques? Спасибі!

— Том

Я підозрюю, що вкладка "Властивості" обманює. Докладніше про це дивіться на gis.stackexchange.com/questions/8780 .

— whuber

Спасибі. Ніколи не думав про можливість того, що вкладка "Властивості" обманює. Можливо, це може бути уточнено, якщо ми знаємо, чому причина ArcGIS дає різні значення для Streched та Unique в Symbology для одного растрового шару. Хтось би знав, чому?

— Том

Як було сказано, часто статистичні дані, знайдені у властивостях растру, іноді є приблизними або застарілими. Вони є заздалегідь визначеними властивостями, які можуть вводити в оману фактичні значення растру.

Обчислили власні значення min / max від 100% фактичних даних за допомогою масивів NumPy. Див. Розділ Робота з NumPy в ArcGIS та RasterToNumPyArray (arcpy) . Наприклад:

import arcpy

inrast = r'C:\data\inRaster.tif'
my_array = arcpy.RasterToNumPyArray(inrast)
print((my_array.min(), my_array.max()))

Якщо у вас відсутні значення (NODATA), маскується масив, щоб отримати правильну статистику:

import numpy as np
my_array = arcpy.RasterToNumPyArray(inrast)
my_masked_array = np.ma.masked_equal(my_array, arcpy.Raster(inrast).noDataValue)
print((my_masked_array.min(), my_masked_array.max()))

Крім того, вам не потрібен ArcGIS, щоб читати растри як масиви NumPy; наприклад, GDAL або растеріо можуть зробити подібне .

— Майк Т
джерело

Як свідчать інші відповіді, статистика, ймовірно, застаріла. Якщо ви віддаєте перевагу використанню ArcGIS, спробуйте інструмент Обчислити статистику в панелі інструментів управління даними. Це має оновити статистику для вас.

— Мартін
джерело