Заливка (Fill) (Spatial Analyst)
Резюме
Заполняет локальные понижения в растре поверхности для удаления всех небольших ошибок и неточностей, присущих данным.
Использование
-
Локальное понижение – это ячейка с неопределенным направлением стока; вокруг нее не существует ячеек с меньшей, чем у этой ячейки, высотой. Точка устья – это пограничная ячейка с наименьшей высотой для водосборной области локального понижения. Если бы локальные понижения были заполнены водой, в этих точках вода «утекала бы» с поверхности.
-
Ограничение по z задает максимальную допустимую разницу между глубиной точки локального понижения и устья и определяет, какое понижение будет заполнено, а какое останется незатронутым. Ограничение по z не является значением максимальной глубины, до которой будет заполнено локальное понижение.
Например, рассмотрим область локального понижения, где точка устья находится в 210 футах по высоте, а самая глубокая точка локального понижения – 204 фута (то есть разница составляет 6 футов). Если ограничение по z равно 8, это конкретное локальное понижение будет заполнено. Однако, если ограничение по z равно 4, это локальное понижение не будет заполнено, поскольку глубина понижения превышает разницу и оно будет считаться допустимым понижением.
-
Все локальные понижения, которые по глубине меньше, чем ограничение по z, и ниже, чем самая низкая из прилегающих соседних ячеек, будут заполнены до высоты своих точек устьев.
-
Запуск инструмента Заполнение (Fill) может быть требователен к процессору, оперативной памяти и занимаемому месту на диске. Он может требовать до четырех раз большего объема дискового пространства, чем входной растр.
-
Число локальных понижений, определяемых с использованием ограничения по z, будет задавать время, требуемое на обработку. Чем больше локальных понижений имеется, тем больше потребуется времени на обработку.
-
Инструмент Локальное понижение (Sink) может быть использован перед тем, как запустить инструмент Заполнение (Fill), чтобы найти число локальных понижений и помочь определить их глубину. Знание глубины локальных понижений может помочь в определении подходящего значения ограничения по z.
-
Инструмент Заполнение локальных понижений может быть также использован для удаления пиков. Пик – это ячейка, вокруг которой нет ячеек выше нее. Для удаления пиков входной растр поверхности должен быть инвертирован. Эту операцию можно выполнить с помощью инструмента Минус (Minus). Задайте наибольшее значение растра поверхности в качестве первых входных данных инструмента Вычитание (Minus) и растр поверхности в качестве вторых входных данных. Выполните заполнение. Инвертируйте результаты с тем, чтобы получить поверхность, имеющую исходные значения растра поверхности с удаленными пиками. Ограничение по z может быть применено и к этому процессу. Если для ограничения по z не задано никакого значения, будут удалены все пики. Если значение определено, в тех случаях, когда разница в значении z между пиком и самой высокой соседней ячейкой больше, чем ограничение по z, этот пик удален не будет.
Синтаксис
| Параметр | Объяснение | Тип данных |
in_surface_raster | Входной растр непрерывной высотной поверхности. | Raster Layer |
z_limit (дополнительно) | Максимальная разница в высотах между заполняемым локальным понижением и его точкой устья. Если разница в z-значениях между локальным понижением и его точкой устья больше, чем значение z_limit, это локальное понижение не будет заполнено. Значение ограничения по z должно быть больше нуля. Если для этого параметра не задано значение, все локальные понижения будут заполнены, независимо от их глубины. | Double |
Возвращено значение
| Имя | Объяснение | Тип данных |
| out_surface_raster |
Выходной набор растровых данных поверхности после заполнения локальных понижений. | Raster |
Пример кода
В этом примере заполняются локальные понижения входного растра поверхности рельефа Grid.
import arcpy
from arcpy import env
from arcpy.sa import *
env.workspace = "C:/sapyexamples/data"
outFill = Fill("elevation")
outFill.save("C:/sapyexamples/output/outfill01")
В этом примере заполняются локальные понижения входного растра поверхности рельефа Grid с применение ограничения по z.
# Name: Fill_Ex_02.py
# Description: Fills sinks in a surface raster.
# Requirements: Spatial Analyst Extension
# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *
# Set environment settings
env.workspace = "C:/sapyexamples/data"
# Set local variables
inSurfaceRaster = "elevation"
zLimit = 3.28
# Check out the ArcGIS Spatial Analyst extension license
arcpy.CheckOutExtension("Spatial")
# Execute FlowDirection
outFill = Fill(inSurfaceRaster, zLimit)
# Save the output
outFill.save("C:/sapyexamples/output/outfill02")