Выравнивание двух несовпадающих растровых сеток с равным разрешением в ArcGIS Desktop?

14

У меня есть два растра: цифровая модель рельефа (DEM) и группа Landsat. Разрешение сетки обоих составляет 30 м, но сетки не совпадают (см. Изображение: темный растр наверху - DEM; светлый растр внизу - полоса Landsat).

Я хочу выполнить некоторое преобразование полосы Landsat, чтобы ее сетка идеально совпала с сеткой матрицы высот. Я пытался пересэмплировать полосу Landsat до 30 м в ArcGIS, используя DEM в качестве Snap Raster в настройках Environment, но это не имеет никакого эффекта. Я также попытался обрезать полосу Landsat до размера векторного многоугольника (к которому также подрезана матрица высот), при этом установив матрицу высот как Snap Raster, но безрезультатно.

Я подозреваю, что это происходит, потому что оба растра уже имеют одинаковое разрешение (только на разных сетках), но мне действительно нужно выровнять их сетки, чтобы иметь возможность проводить дальнейший анализ.

Могу ли я сделать это в ArcGIS Desktop?

введите описание изображения здесь

Hendra
источник
2
Не все инструменты соответствуют настройкам среды Snap Raster (хотя Resample говорит, что это так, поэтому я не уверен, почему это не работает - возможно, потому что размер плитки не изменился (?)). Вы пробовали Copy Raster с включенной настройкой Snap Raster? Предположительно, он соблюдает настройки. Кстати, какую версию ArcGIS вы используете и какую версию GDAL?
om_henners
1
@om_henners: очевидно, проблема заключается в ошибке в ArcGIS 10. Смотрите мой ответ ниже.
Хендра

Ответы:

6

После дальнейших исследований я обнаружил, что неправильная настройка Snap Raster на самом деле была вызвана ошибкой в ​​ArcGIS 10. Похоже, что решение состоит в том, чтобы отключить фоновую геообработку (параметры геообработки \ параметры геообработки ... затем снимите флажок Включить в разделе Фоновая обработка ), и выполните операцию отсечения, с указанием Snap Raster, в ArcCatalog.

На форумах ArcGIS есть обсуждение здесь: http://forums.arcgis.com/threads/50808-Another-Snap-to-Raster-problem

Очевидно, ошибка исправлена ​​в ArcGIS SP 4. Я установил пакет обновления, но еще не тестировал.

Интересно (или, скорее, досадно), если я выполняю клип на моем полигоне, представляющем интерес, с указанным Snap Raster (Snap Raster - DEM) в соответствии с решением в сообщении на форуме выше, Snap Raster работает, но Landsat В band есть еще один столбец и еще одна строка, чем растр DEM, который также был обрезан в ArcGIS до полигона интересующей области. Это раздражает, потому что для дальнейшего анализа мне нужно, чтобы два растра имели одинаковое количество строк и столбцов. Простой обходной путь - выполнить вычисление с помощью Con()функции, которая создает новый растр с тем же количеством строк и столбцов, что и у матрицы высот, но который имеет значения полосы Landsat.

Hendra
источник
2

Я долго боролся с этой проблемой, и я нашел то, что считаю очень эффективным и довольно простым решением. Выберите один из ваших растров, который, по вашему мнению, находится в «правильном» месте, увеличьте настолько, насколько это позволяет ArcGIS, и создайте новый точечный шейп-файл прямо поверх пересечения четырех растровых ячеек. Единственный трюк в этом заключается в том, чтобы быть уверенным в том, что вы увеличены настолько, насколько это возможно, и что вы создаете новую точку прямо поверх «креста», образованного четырьмя растровыми ячейками. Затем запишите координаты точки, которую вы только что создали. Наконец, перепроектируйте ваш второй растр (но оставьте размер проекции и выходной ячейки без изменений) и используйте координаты вновь созданной точки в качестве «Контрольной точки». Это отлично выровняет ваши два растра. Нет ошибок,

Мэтью Лонг
источник
1

Вы можете попробовать использовать инструмент маскирования для наложения матрицы высот с помощью Landsat.

Возможно, вам придется создать окно области исследования, чтобы захватить пиксели, которые полностью падают друг на друга.

LJF
источник
0

Измерьте зазор между одним растром и другим, а затем сместите один растр соответственно.

Вы можете использовать инструмент Shift из панели инструментов управления данными.

nickves
источник
0

У меня такая же проблема. Я попробовал решение Louis Longchamps , которое прекрасно работало на небольшом тестовом образце, но мой набор данных был слишком большим, чтобы использовать этот метод практически. После долгого времени и множества способов сделать это, я наконец-то по прихоти использовал «Извлечь по прямоугольнику» после повторной выборки, которая работала так же, как должна была работать функция привязки растра в инструменте клипа. Я настроил и попытался использовать построитель моделей, но это дало мне ошибку, так что в итоге пришлось делать это по одному шагу за раз, но в итоге все мои растры были согласованы.

Роберт Махол
источник
0

Мне не нравилась идея смещения растра, потому что, по-моему, вы меняете информацию. Чтобы избежать этого, я использовал последовательность преобразований.

Вот шаги:

  1. Для растра, который вы хотите использовать в качестве шаблона, зайдите в Properties/Sourceи запишите Cellsizeна листе бумаги все доступные десятичные дроби.
  2. Преобразуйте этот растр в точки, используя Conversion Tools/Raster to Point.
  3. Извлеките значения из всех других растров, которые вы хотите использовать Spatial Analyst Tools/Extraction/Extract Multi Values to Points, указав только что созданный точечный объект (на шаге 2) в качестве входного точечного объекта и всех растров, которые вы хотите выбрать во входных растрах.
  4. Создайте растры из этих точек по очереди Conversion Tools/To Raster/Point to Raster, указав, какой слой вы создаете в поле Значение, и используя соответствующее имя файла И Cellsizeзначение, которое вы записали на листе бумаги.

Они должны идеально совмещаться, и вы можете перенести их, например, в R и иметь матрицы и векторы одинаковых размеров.

Луи Лоншам
источник
1
Правильно выполненное растровое «смещение» просто меняет два числа: координаты x и y своего начала. Это делает его чрезвычайно эффективной надежной операцией, которая не изменяет никакие данные в растре.
whuber
0

Мы были так обеспокоены нашим DEM-изображением, что оно также не было выровнено с контролируемым классифицированным растровым изображением. Мы также использовали метод «Мэтью Лонг» .

мы добавили снимки не выровненных и выровненных изображений.

После метода Мэтью Лонга наше DEM-изображение было полностью выровнено с нашим контролируемым классифицированным растровым изображением.

Это изображение инструмента, который мы использовали для выравнивания нашего изображения DEM, в ARCGIS 10.2.1 в инструментах управления данными:

введите описание изображения здесь

Это изображение перед выравниванием, черное изображение является ЦМР, а цветное изображение контролируется классифицированным растровым изображением.

введите описание изображения здесь

Это изображение после выравнивания, черное изображение является ЦМР, а цветное изображение контролируется классифицированным растровым изображением.

введите описание изображения здесь

Анила и Набила
источник
0

Я решаю свою проблему с помощью метода «Мэтью Лонг». Однако сделать это еще проще и 100% точность. Без какого-либо увеличения, просто используйте верхнюю левую координату экстента хорошо выровненного растра.

Для этого: щелкните правой кнопкой мыши, свойство, панель «Источник» под свойством экстента

Кроме того, используя верхнюю левую координату экстента, можно создать решение Python без ручного вмешательства.

Алекс С.
источник
-1

У меня была похожая проблема с двумя файлами из разных источников. Во-первых, растровый файл для контрольного уровня обезлесения на страновом уровне (Перу). Во-вторых, консолидированный файл на платформе GEE загружен и перепроецирован в систему первого элемента. В заключение, наложение обоих показало мне, что было смещение между ячейками (30 * 30 ячейка) не более 5 метров. Чтобы выровнять и настроить второй файл по первому, мне пришло в голову установить * расширенную основу первого файла, которая может охватывать оба файла * и поддерживать выравнивание по пикселям первого файла с его системой проекции и другими функциями , Таким образом, во время перепроектирования моего второго файла или любого побочного продукта с аналогичным смещением,

Л. Кампос
источник