Изменение значений интересующей области растра в ArcGIS: различия между версиями

Материал из GIS-Lab
Перейти к навигации Перейти к поиску
 
(не показаны 4 промежуточные версии этого же участника)
Строка 2: Строка 2:
{{Аннотация|Создание "заплатки" и склеивание мозаики}}
{{Аннотация|Создание "заплатки" и склеивание мозаики}}
===Введение===
===Введение===
На равнинной территории снимка высот (данные SRTM4) находится техногенный карьер, представляющий из себя отрицательную форму рельефа, вносящий искажения в статистику растра, а отсюда и в восприятие диапазона высот природной составляющей территории. Задача придания области карьера высот, характерных прилегающим территориям, сводиться к созданию т.н. "заплатки" и получению мозаики на основе исходного растра высот.<br />
На равнинной территории снимка высот (данные SRTM4) находится техногенный карьер, представляющий из себя отрицательную форму рельефа, вносящий искажения в статистику растра, а отсюда и в восприятие диапазона высот природной составляющей территории. Задача придания области карьера высот, характерных прилегающим территориям, сводиться к созданию т.н. "заплатки" и получению мозаики на основе исходного растра высот.[[Файл:Quarry.jpg|мини|справа|750|Рис. 1. Карьер. Landsat ETM+ (5,3,1,+8)]]<br />
В заметке уточнения обозначаются "'''*'''", а также присутствуют следующие сокращения:<br />
В заметке уточнения обозначаются "'''*'''", а также присутствуют следующие сокращения:<br />
'''srtm4.tif''' - исходные данные SRTM4<br />
'''srtm4.tif''' - исходные данные SRTM4<br />
Строка 11: Строка 11:


===Алгоритм действий===
===Алгоритм действий===
[[Файл:Quarry_fill.jpg|мини|справа|750|Рис. 2. Алгоритм создания заплатки]]
1. Создаём новый слой mask.shp в координатной системе такой же, как у srtm4.tif и обводим область, которую хотим убрать (наш карьер).<br />
1. Создаём новый слой mask.shp в координатной системе такой же, как у srtm4.tif и обводим область, которую хотим убрать (наш карьер).<br />
2. В атрибутивной таблице mask.shp создаём поле GridID и вписываем одно значение пикселя (основываясь на значениях прилегающих территорий).* Закрываем сессию редактирования с сохранением правок.
2. В атрибутивной таблице mask.shp создаём поле GridID и вписываем одно значение пикселя (основываясь на значениях прилегающих территорий).* Закрываем сессию редактирования с сохранением правок.
Строка 16: Строка 17:
3. Конвертируем функцией Feature to Raster, где в Field указываем GridID, а в Output cell size выбираем srtm4.tif.* В Output raster присваиваем название mask.tif.
3. Конвертируем функцией Feature to Raster, где в Field указываем GridID, а в Output cell size выбираем srtm4.tif.* В Output raster присваиваем название mask.tif.
* ''мы конвертируем в растровый формат, имеющий свойства srtm4.tif''
* ''мы конвертируем в растровый формат, имеющий свойства srtm4.tif''
4. ПКМ на mask.tif->Data->Export Data..., где указываем значение NoData as: такое же, как и у srtm4.tif*. Получаем mask1.tif.
4. ПКМ на mask.tif->Data->Export Data..., где указываем значение NoData as: такое же, как и у srtm4.tif.* Получаем mask1.tif.
* ''значение NoData для srtm4.tif можно посмотреть в свойствах слоя. Этот шаг окончательно унифицирует свойства mask.tif с srtm4.tif (их pixel depth и nodata value) и нужен для получения в следующем шагу elevation.tif без хаотичных дыр NoData по снимку.''
* ''значение NoData для srtm4.tif можно посмотреть в свойствах слоя. Этот шаг окончательно унифицирует свойства mask.tif с srtm4.tif (их pixel depth и nodata value) и нужен для получения в следующем шагу elevation.tif без хаотичных дыр NoData по снимку.''
5. Используем функцию Mosaic to new raster, первым добавляем mask1.tif, потом srtm4.tif. Устанавливаем pixel depth, band и выбираем mosaic operator - first (заплатка у нас будет наложена поверх). Получаем результирующий растр elevation.tif.<br />
5. Используем функцию Mosaic to new raster, первым добавляем mask1.tif, потом srtm4.tif. Устанавливаем pixel depth, band и выбираем mosaic operator - first (заплатка у нас будет наложена поверх). Получаем результирующий растр elevation.tif.<br />
6. Функцией Build raster attribute table строим атрибутивную таблицу для elevation.tif и получаем новую.<br />
6. Функцией Build raster attribute table строим атрибутивную таблицу для elevation.tif и получаем новую.<br />
В данном случае elevation.tif и будет растром представляющим диапазон высот природной составляющей территории.
В данном случае elevation.tif и будет растром представляющим диапазон высот природной составляющей территории.

Текущая версия от 15:00, 24 апреля 2013

Эта страница является черновиком статьи.


Создание "заплатки" и склеивание мозаики

Введение

На равнинной территории снимка высот (данные SRTM4) находится техногенный карьер, представляющий из себя отрицательную форму рельефа, вносящий искажения в статистику растра, а отсюда и в восприятие диапазона высот природной составляющей территории. Задача придания области карьера высот, характерных прилегающим территориям, сводиться к созданию т.н. "заплатки" и получению мозаики на основе исходного растра высот.

Рис. 1. Карьер. Landsat ETM+ (5,3,1,+8)


В заметке уточнения обозначаются "*", а также присутствуют следующие сокращения:
srtm4.tif - исходные данные SRTM4
mask.shp - шейпфайл заплатки
mask.tif - растр заплатки
mask1.tif - растр заплатки со свойствами исходных данных
elevation.tif - результирующий растр

Алгоритм действий

Рис. 2. Алгоритм создания заплатки

1. Создаём новый слой mask.shp в координатной системе такой же, как у srtm4.tif и обводим область, которую хотим убрать (наш карьер).
2. В атрибутивной таблице mask.shp создаём поле GridID и вписываем одно значение пикселя (основываясь на значениях прилегающих территорий).* Закрываем сессию редактирования с сохранением правок.

  • это подготовка к получению Constant Raster

3. Конвертируем функцией Feature to Raster, где в Field указываем GridID, а в Output cell size выбираем srtm4.tif.* В Output raster присваиваем название mask.tif.

  • мы конвертируем в растровый формат, имеющий свойства srtm4.tif

4. ПКМ на mask.tif->Data->Export Data..., где указываем значение NoData as: такое же, как и у srtm4.tif.* Получаем mask1.tif.

  • значение NoData для srtm4.tif можно посмотреть в свойствах слоя. Этот шаг окончательно унифицирует свойства mask.tif с srtm4.tif (их pixel depth и nodata value) и нужен для получения в следующем шагу elevation.tif без хаотичных дыр NoData по снимку.

5. Используем функцию Mosaic to new raster, первым добавляем mask1.tif, потом srtm4.tif. Устанавливаем pixel depth, band и выбираем mosaic operator - first (заплатка у нас будет наложена поверх). Получаем результирующий растр elevation.tif.
6. Функцией Build raster attribute table строим атрибутивную таблицу для elevation.tif и получаем новую.
В данном случае elevation.tif и будет растром представляющим диапазон высот природной составляющей территории.