r.stream.extractgrass - Интернет в облаке

ПРОГРАММА:

ИМЯ

r.stream.extract - Выполняет извлечение потоковой сети.

КЛЮЧЕВЫЕ СЛОВА

растр, гидрология, водотоковая сеть

СИНТАКСИС

r.stream.extract
r.stream.extract --Помогите
r.stream.extract высота=имя [накопление=имя] [депрессии=имя] порог=плавать
[d8cut=плавать] [мэксп=плавать] [длина_потока=целое] [Память=целое]
[поток_растр=имя] [поток_вектор=имя] [направление=имя] [-затирать] [-помощь]
[-подробный] [-тихий] [-ui]

Флаги:
- перезаписать
Разрешить выходным файлам перезаписывать существующие файлы

--Помогите
Распечатать сводку использования

--подробный
Подробный вывод модуля

--тихий
Тихий выход модуля

--уи
Принудительный запуск диалогового окна GUI

Параметры:
высота=имя [требуется]
Имя входной растровой карты высот

накопление=имя
Имя входной растровой карты накопления
При извлечении потока будет использоваться полученное накопление вместо повторного расчета.

депрессии=имя
Название исходной растровой карты с реальными впадинами
Потоки не будут выведены из настоящих депрессий

порог=плавать [требуется]
Минимальное накопление потока для ручьев
Должно быть> 0

d8cut=плавать
Используйте SFD выше этого порога
Если накопление больше, чем d8cut, вместо MFD используется SFD. Применяется только если нет
дана карта накопления.
По умолчанию: бесконечность

мэксп=плавать
Показатель Монтгомери для наклона, отключен 0
Монтгомери: накопление умножается на pow (slope, mexp), а затем сравнивается с
порог
По умолчанию: 0

длина_потока=целое
Удалить сегменты потока короче ячеек stream_length
Применяется только к сегментам потока первого порядка (пружины / головки потока)
По умолчанию: 0

Память=целое
Максимальный объем используемой памяти (в МБ)
Размер кеша для растровых строк
По умолчанию: 300

поток_растр=имя
Имя выходной растровой карты с уникальными идентификаторами потоков

поток_вектор=имя
Имя выходной векторной карты с уникальными идентификаторами потоков

направление=имя
Имя выходной растровой карты с направлением потока

ОПИСАНИЕ

r.stream.extract извлекает потоки как в растровом, так и в векторном формате из требуемого ввода
высота карта и дополнительный ввод накопление карта.

ПРИМЕЧАНИЯ

Ячейки NULL (nodata) во входных данных высота map игнорируются, нулевые и отрицательные значения
действительные данные о высоте. Разрывы на карте высот, расположенные в районе
проценты должны быть заполнены заранее, например, р.fillnulls, чтобы избежать искажений.

Все ненулевые и ненулевые ячейки депрессии карту будем рассматривать как настоящие впадины.
Потоки не будут выведены из впадин. Если область отмечена как депрессия, но
Модель рельефа не имеет депрессии в этом месте, потоки на этом не остановятся. Если поток
Приведены карта накопления и карта с реальными депрессиями, карта накопления стока
должен соответствовать карте депрессии, чтобы поток не выходил за пределы указанного
депрессии. Рекомендуется использовать внутреннее вычисление накопления расхода, если
предоставляется карта депрессии.

Опция порог определяет минимальное (необязательно измененное) значение накопления потока, которое
инициирует новый поток. Если используется метод Монтгомери для инициации канала,
значение ячейки входной карты накопления умножается на (tan (локальный наклон)) mexp, а затем
в сравнении с порог. Если мэксп дан, чем метод Монтгомери и
Foufoula-Georgiou (1993), чтобы инициировать поток с этим значением. Значение ячейки
входная карта накопления умножается на (tan (местный наклон)) mexp и затем сравнивается с
порог. Если порог достигнут или превышен, инициируется новый поток. По умолчанию
значение 0 отключает Монтгомери. Монтгомери и Фуфула-Георгиу (1993) обычно рекомендуют
использовать 2.0 как показатель степени. мэксп значения ближе к 0 создадут потоки, более похожие на
потоки, извлеченные с отключенным Монтгомери. Больше мэксп значения уменьшают количество
ручьи на равнинных участках и увеличивают количество ручьев на крутых участках. Если вес is
учитывая, груз применяется первым.

Опция d8cut определяет минимальный объем сухопутного стока (накопления), когда SFD (D8) будет
используется вместо MFD (FD8) для расчета накопления потока. Применяется только при отсутствии накопления
карта предоставляется. Установка на 0 полностью отключает МФД.

Опция длина_потока определяет минимальную длину потока в количестве ячеек для первого порядка
(головной / родниковый) сегменты ручья. Все сегменты потока первого порядка короче, чем длина_потока
будет удален.

Результат направление растровая карта содержит направление потока для всех ячеек, отличных от NULL, во входных данных
высота. Направление потока имеет тип D8 с диапазоном от 1 до 8. Умножение значений на
45 - градус против часовой стрелки с востока. Направление потока регулировали во время прореживания, принимая
ярлыки и пропуск ячеек, которые были устранены процедурой прореживания.

Поток добыча
Если нет накопление предоставлена ​​входная карта, накопление потока определяется с помощью
гидрологический анализ, аналогичный правый водораздел. Алгоритм - MFD (FD8) по Холмгрену.
1994, что касается правый водораздел, порог опция определяет количество потоков и детализацию
потоковых сетей. Когда накопление потока достигает порог, новый поток запущен
и прослеживается вниз по потоку до точки выхода. Что касается правый водораздел, накопление потока
рассчитывается как количество ячеек, истекающих через ячейку.

If накопление дается, чем накопленные значения предоставленных накопление карта
используется и не рассчитывается на основе ввода высота карта. В этом случае высота карта должна
быть точно такой же картой, которая использовалась для расчета накопление. Если накопление был рассчитан
р.террафлоу, заполненный выход высот р.террафлоу должны быть использованы. Далее
текущий регион должен быть выровнен по накопление карта. Направление потока первое
рассчитывается из высота а затем настроен на накопление. Нет необходимости
обеспечивать накопление как количество ячеек, оно также может быть отрегулировано по желанию или
общая взвешенная площадь в квадратных метрах или в любых других единицах. Когда оригинальный поток
карта накопления скорректирована или взвешена, юстировка или взвешивание не должны преобразовываться
допустимые значения накопления до значений NULL (nodata).

Взвешенный поток накопление
Сначала можно рассчитать накопление потока, например, с помощью правый водораздел, а затем изменен до
используя его как вход для r.stream.extract. В общем виде карта взвешенного накопления имеет вид
генерируется путем создания сначала карты взвешивания, а затем умножения карты накопления на
карта взвешивания с использованием r.mapcalc. Настоятельно рекомендуется оценить взвешенный поток.
карту накопления, прежде чем использовать ее в качестве входных данных для r.stream.extract.

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

Другая возможность - ограничить начало канала долинами, определенными на местности.
морфология. Долины можно определить с помощью r.param.scale param = crosc (поперечное сечение или
тангенциальная кривизна). Значения кривизны <0 указывают на вогнутые элементы, то есть впадины. В
размер окна обработки определяет, будут ли узкие или широкие впадины
идентифицированы (см. пример ниже).

Определяющий a область of интерес
Процедура извлечения потока может быть ограничена определенной интересующей областью, например
суббассейн, задав расчетную область с г. регион и / или создание МАСКИ. Такой
представляющая интерес область должна быть полностью обслуживаемой, в том смысле, что
Включена вся площадь перед выходным отверстием и буферизуется, по крайней мере, с одной ячейкой.

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

Результат поток_растр растровая карта хранит извлеченные потоки. Значения ячеек кодируют уникальный идентификатор
для каждого сегмента потока.

Результат поток_вектор Векторная карта хранит извлеченные сегменты и точки потока. Очки
записывается в начале каждого сегмента потока и на выходе из потока
сеть. В слое 1 категории представляют собой уникальные идентификаторы, идентичные значению ячейки растра.
выход. Таблица атрибутов для уровня 1 содержит информацию о типе потока.
сегмент: начальный сегмент или промежуточный сегмент с притоками. Столбцы - это cat int,
stream_type varchar (), type_code int. Кодировка для type_code: 0 = начало, 1 =
средний. На уровне 2 категории идентичны type_code на уровне 1 с дополнительными
категория 2 = выход для точек выхода. Баллы с категорией 1 = средний уровень 2
находятся в месте слияния.

ПРИМЕР

Этот пример основан на карте высот «elev_ned_30m» из образца Северной Каролины.
набор данных и использует долины, определенные с помощью r.param.scale взвесить карту накопления
производится с правый водораздел.
# установить регион
g.region -p raster = elev_ned_30m @ ПОСТОЯННЫЙ
# вычислить накопление потока
р. водораздел ele = elev_ned_30m @ PERMANENT acc = elevation.10m.acc
# кривизна для получения узких впадин
r.param.scale input = elev_ned_30m @ PERMANENT output = tangential_curv_5 size = 5 param = crosc
# кривизна, чтобы сделать долины немного шире
r.param.scale input = elev_ned_30m @ PERMANENT output = tangential_curv_7 size = 7 param = crosc
# кривизна, чтобы получить широкие впадины
r.param.scale input = elev_ned_30m @ PERMANENT output = tangential_curv_11 size = 11 param = crosc
# создать карту веса
r.mapcalc "weight = if (tangential_curv_5 <0, -100 * tangential_curv_5, \
if (tangential_curv_7 <0, -100 * tangential_curv_7, \
if (tangential_curv_11 <0, -100 * tangential_curv_11, 0.000001))) "
# карта накопления веса
r.mapcalc expr = "elev_ned_30m.acc.weighed = elev_ned_30m.acc * вес"
# копируем таблицу цветов из исходной карты накопления
r.colors map = elev_ned_30m.acc.weighed raster = elev_ned_30m.acc
Отображение исходной и взвешенной карты накопления. Сравните их и продолжайте, если
карта взвешенного накопления имеет смысл.
# извлечь потоки
r.stream.extract elevation = elev_ned_30m @ ПОСТОЯННЫЙ \
накопление = высота_нед_30м. взвеш. \
порог = 1000 \
stream_rast = elev_ned_30m.streams
# извлекаем потоки, используя исходную карту накопления
r.stream.extract elevation = elev_ned_30m @ ПОСТОЯННЫЙ \
накопление = отметка_нед_30м.acc \
порог = 1000 \
stream_rast = elev_ned_30m.streams.noweight
Теперь отобразите обе карты потоков и решите, какая из них более реалистична.

Ссылки

· Эльшлегер, К. (1989). . AT Поиск Алгоритм в Разработка Гидрологический
Модели от Digital Высота Данные, Производство of международный Географический
Информация системы (ИГИС) КОНФЕРЕНЦИЯ ПО СИНЕСТЕЗИИ. МОСКВА, XNUMX-XNUMX ОКТЯБРЯ, XNUMX '89, стр. 275-281 (Балтимор, Мэриленд, 18-19 марта
1989). URL-адрес: http://faculty.wiu.edu/CR-Ehlschlaeger2/older/IGIS/paper.html

· Холмгрен, П. (1994). Многочисленные поток направление алгоритмы для сток моделирование in
сетка основанный высота модели: An эмпирический оценка. гидрологический Процессы Vol
8(4), стр. 327-334. DOI: 10.1002 / hyp.3360080405

· Монтгомери, Д. Р., Фуфула-Георгиу, Э. (1993). Канал сеть источник
представление через Интернет высота модели. Вода Ресурсы Исследование Vol
29(12), с. 3925-3934.

Используйте r.stream.extractgrass онлайн с помощью сервисов onworks.net