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