Это в некоторой степени теоретический вопрос, вытекающий из некоторых дискуссий с коллегами по теме последствий разграничения бассейнов с прогнозируемыми (например, равными по площади Альберса) и непроектированными (NAD 83) данными, полученными из 10-метровой матрицы высот, находящейся в NAD 83.
Некоторые утверждают, что это не проблема, поскольку значения, рассчитанные на основе непроецированных данных, просто корректируются, если вы решите спроецировать.
Я не уверен, что это так, поскольку существуют различия между данными в географической системе координат и проецируемыми данными. Я попробовал один пример прохождения процедуры, начиная с непроецированных данных матрицы высот, а затем протестировал тот же сайт с прогнозными данными матрицы высот. Шаги, предпринятые для обоих, были выполнены (вся работа выполнена в ArcGIS 9.3.1) с использованием данных высотой 10 м.
Один запуск был выполнен с использованием матрицы высот в NAD 83, а второй - путем проецирования той же матрицы высот в USA_Contiguous_Albers_Equal_Area_Conic_USGS_version.
- определить направление потока с помощью инструмента геообработки FlowDirection_sa
- получить накопление потока с помощью инструмента FlowDirection_sa
- установить точку застывания на расстоянии 50 метров
- определить водораздел с помощью инструмента Watershed_sa
При сравнении этих двух я мог заметить визуальную разницу между отображением сеток направления потока.
ПРИМЕЧАНИЕ. После дальнейших исследований я полагаю, что эффект чередования обусловлен не использованием повторной выборки CUBIC, а ошибочным выбором значения по умолчанию NEAREST в инструменте ArcGIS Project Raster. Я не верю, что это дает какое-то разрешение этой дискуссии, хотя ...
Направления потока с использованием незапроектированной матрицы высот
Направления потока с использованием спроектированной матрицы высот
Я понимаю, что визуальное сравнение не на 100% научно, но может быть хорошей отправной точкой.
Соответственно, была разница между температурой застывания и тем, как она щелкала при каждом прогоне. И было определенное различие в производных водоразделах, учитывая, как инструмент мгновенного определения точки текучести решил выполнить привязку на основе соответствующих спроецированных / непроецированных наборов данных. Водораздел, показанный зеленым, является водоразделом, полученным с использованием прогнозируемой матрицы высот и последующих прогнозируемых производных данных высоты. Водораздел, показанный фиолетовым контуром, является водоразделом, полученным с использованием не спроецированных данных матрицы высот.
Водораздел
Я сталкивался с этими двумя другими ветками форума ГИС (ссылки ниже), которые обсуждают эту проблему на старых форумах ESRI, но мне все еще неясно, как работает инструмент «Направление потока» по отношению к прогнозируемым и непроектированным данным (я понимаю, концепция гидрологического потока и направления потока, хотя). Если каждая ячейка по-прежнему имеет одно и то же значение высоты в проецируемой ЦМР по сравнению с непроецированной ЦМР (верно ли это?), Почему существует разница в растре направления потока, полученном из проецируемых данных, и растре, полученном из данных ЦМР в NAD83?
http://forums.esri.com/Thread.asp?c=93&f=995&t=292503
http://forums.esri.com/Thread.asp?c=93&f=995&t=290652
Кроме того, не будет ли теоретически различий меньше, если проводить разграничения в более высоких широтах, таких как национальный парк Шенандоа в Вирджинии, по сравнению с разграничениями в штате Техас?
Я говорил с одним экспертом по картографии, который подумал, что искажение восток-запад, которое вы получаете при удалении от экватора, может быть проблемой (например, как на некоторых картах Канада сильно раздута и искажена), в том случае, если вы более На 10 градусов широты от экватора они думали, что проецируемые данные - это путь, если вы заинтересованы в точности.
Одним из главных неизвестных является уровень неопределенности с бассейнами, очерченными с использованием непроецированных данных, с которыми мы пытаемся справиться. Есть разница, но какова величина?
Спасибо всем, кто может дать прямой ответ на эту дискуссию или просто полезную информацию об этом.
редактировать
Основной вопрос, который нас интересует / который волнует, заключается в том, возникнут ли проблемы точности с обозначенными водоразделами в результате запуска процесса с использованием незапроектированной ЦМР.
Итак, если я понимаю ответ, то обозначенные бассейны должны быть хорошими с точки зрения представления площади дренажа для точки застывания? Кажется, однако, что, если направления потока неправильны, это приведет к некоторой ошибке в конечном очерченном водоразделе.
Это очень интересная и действительно важная тема - мне еще предстоит увидеть отчет или документацию, подтверждающую, что можно использовать прогнозируемые ООН данные для определения водоразделов. Я провел технические конференции ESRI User Conference во главе с ведущим инженером-разработчиком расширения Spatial Analyst (в котором находятся инструменты гидрологии), где они сказали, что вы должны также использовать проекцию равных площадей (например, равную площадь Альберса).
Кроме того, кажется, что не существует какого-либо авторитетного «библейского» стандарта для того, как это сделать - просто кажется, что это практически общепризнанный подход де-факто для проецирования данных до вычисления ваших производных высот.
Нигде я не смог найти краткого и прямого ответа о том, как это влияет на расчет направления потока и, следовательно, на определение водораздела.
И, если вы в конечном итоге работаете с водоразделами, очерченными с использованием не спроецированных данных ЦМР, и затем проектируете эти водоразделы, разве нет еще неточности (например, с точки зрения определения площади водосбора или любых других характеристик, таких как пропорции земного покрова и т. Д.)?
Кроме того, я предполагаю, что проецирование растра направления потока, который был получен из не спроектированной матрицы высот, также не исправляет ошибки, так как исходные данные не были спроецированы ....
спасибо - оцените любую дополнительную информацию, которую вы можете предоставить
РЕДАКТИРОВАТЬ - 20110331
@whuber:
спасибо за это обширное обсуждение. Мы исследовали эту проблему больше и на самом деле столкнулись с некоторыми рекомендациями, которые предполагают, что на самом деле лучше не проецировать ЦМР, прежде чем получать поток управления, накопление потока и разграничение.
Один ответ по электронной почте от анонимного источника (но который является довольно авторитетным человеком), когда был поставлен вопрос 1.) проекта DEM 2.) производства деривативов ИЛИ 1.) производства деривативов 2.) проекта DEM сказал:
В двух словах, это зависит от производной. Для непрерывных производных, которые будут визуализироваться, вы должны получить и затем спроецировать - это снижает риск того, что артефакты границ плитки будут улучшены или введены (с помощью алгоритма проецирования), а затем переданы производной, если вы сначала спроектировали DEM. Исключением является случай, когда вы также используете расстояние или площадь в качестве основы для своего производного расчета. Это, конечно, зависит от того, насколько велики расстояния / площади и как далеко вы можете приемлемо уйти от экватора. Итак, представьте, что для производных, таких как склон или склон, которые зависят от размера клетки, есть последствия. Эти производные будут наиболее точными на экваторе, и точность значительно снизится после 60 градусов севернее или югу. В обоих случаях я предполагаю, что ЦМР покрывает очень большую площадь (шире, чем 1,5 зоны UTM) и традиционный подход на основе плиток, где листы либо произвольны, либо соответствуют существующим стандартам, таким как границы четырехугольных листов USGS. Таким образом, сказанное подразумевает, что большая часть этого мышления предшествует наборам данных мозаики, которые я не могу комментировать. Главной заботой для меня было бы желание узнать, насколько хорошо соответствуют плитки DEM. Если они хорошо согласованы (например, NED), то я ожидаю, что все будет работать хорошо, производные будут получены из плиток (как функции, применяемые к набору данных мозаики), а затем они будут отображаться на лету. Если они не соответствуют друг другу, то мусор на входе, мусор на выходе. Возвращаясь к вашему первоначальному вопросу, я думаю, что это просто границы водораздела,
Они продолжали говорить:
Причина, по которой я придерживаюсь непрогнозируемой методологии, заключается в том, что мы используем растры, которые сами по себе являются производными от матрицы высот (которую мы обычно не имеем, но думаем, что облако точек LiDAR). Для растров, которые покрывают очень большие области, например, континенты с относительно хорошим уровнем разрешения, проецирование на что-то вроде Альберса приведет к потере или введению информации, когда растр использует ячейки обычного размера (как это делают растры Эсри). Это означает, что такие инструменты, как Flow Accumulation, будут давать результаты на основе частичной или интерполированной информации. По сути, все алгоритмы проецирования, применяемые к растрам, будут вызывать проблемы, как только происходит расширение или сжатие, превышающее расстояние ширины пикселя (проекции, подобные Альберсу, могут вносить ошибку, вводя новые пиксели между двумя старыми). Из этого следует, что вероятность кумулятивной ошибки высока.
Кажется, это говорит об обратном: проецирование создает больше шума, если вы не достигли 60 градусов широты.
Мы также сталкивались с некоторыми опубликованными источниками, в которых предлагалось подобное, что незапроектированный является приемлемым подходом для небольших водосборов (последние 2 абзаца раздела 1.6) из Распределенного гидрологического моделирования для ГИС (Vieux, 2004): http: //www.springerlink. ком / содержание / x877238532533g20 / fulltext.pdf
Итак, в конце концов, сводится ли это к 1), где вы выполняете работу на поверхности земли, 2) масштабам, на которых вы работаете, и 3) ли шум, создаваемый проекцией это лучше сохранит атрибуты, которые влияют на алгоритм направления потока, меньше, чем искажение, вносимое непроецированными данными (выгода возрастает по мере вашего движения к полюсам), чтобы определить, следует ли вам проецировать что-то вроде конформности или это не имеет значения?
Когда вы начинаете копаться в этой теме, кажется, что больший консенсус заключается в том, чтобы спроектировать, но есть некоторые, которые, кажется, говорят, что это не жесткое и быстрое правило.