Учитывая, что у вас есть LIDAR DEM, вы должны использовать потоки, полученные из него. Это гарантирует безупречную регистрацию.
Суть идеи состоит в том, чтобы оценить средние уклоны с точки зрения возвышений на концах отрезков.
Одной из самых простых процедур является «взорвать» потоковую сеть на неразветвленные дуги ее компонентов. Преобразуйте коллекцию в слой «маршрут» на основе расстояния, сделав его «измеримым». Теперь просто сгенерировать набор «событий» маршрута на основе таблицы этапов (например, с интервалами 100 м) для каждой дуги и извлечь значения высот ЦМР из этих точек событий. Последовательные перепады высот вдоль каждой дуги, разделенные на 100 м, оценивают средние уклоны сегмента.
На следующем рисунке показаны дуги потоков, полученные в результате анализа накопления потока на 7-минутной ЦМР США (часть округа Хайленд, Вирджиния). Это около 10 км в поперечнике (6 миль).
Поскольку вы ищете остаточную плотину, на что может указывать изменение градиента всего за несколько десятков метров (для очень маленькой плотины), рассмотрите возможность использования еще меньших сегментов. Если набор данных слишком грубый, чтобы обеспечить четкие сигналы, вы можете легко отфильтровать его позже (с помощью скользящих средних или иным образом, например, сплайновая графики высот и дифференциация сплайна). По сути, этот подход помещает вас в область анализа временных рядов, где интересующей переменной является высота, а не градиент, и вы ищете шаблоны, состоящие из коротких участков уровня, за которыми следуют внезапные изменения.
Это график высот матрицы высот, наблюдаемых с интервалами 100 м вдоль большинства (не всех) изображенных сегментов потока. (Размер ячейки составляет 30 м.) При необходимости, дуги были переориентированы для уменьшения высоты, как правило, слева направо. (Если вы посмотрите внимательно, вы увидите, где я пропустил одну: она поднимается слева направо.)
Эта деталь дуги 16 (длинный сегмент в верхней части карты) показывает, что вы можете получить, когда потоки не полностью зарегистрированы в ЦМР: в местах, где поток кажется восходящим. Тем не менее, сегменты, предлагающие характеристики пула и падения, легко идентифицируются, особенно после вех 1800 (метры вдоль сегмента), 4000, 4600 и 6500. Эта идентификация может быть автоматизирована различными способами, особенно после очистки ряда возвышений (сглаживанием Это).
Вы можете видеть, что используемый здесь 100-метровый интервал отбора проб недостаточно хорош для определения характеристик, длина которых меньше 400-500 метров. Таким образом, чтобы найти небольшую остаточную плотину, вы, вероятно, захотите взять пробы с интервалом 10-25 м на своей лидере.
Кстати, то, что делает сегмент потока «слишком маленьким» для такого рода работы, - это не короткая длина и не большой размер ячейки, хотя оба играют роль в принятии решения. «Слишком маленький» зависит от того, как вы будете использовать предполагаемые уклоны и насколько неопределенными могут быть эти оценки. Для некоторой работы может иметь смысл оценивать градиенты с интервалами 10 метров по сетке 10 метров!