Доступно с лицензией Spatial Analyst.
Инструмент Топо в растр (Topo to Raster) – это метод интерполяции, специально разработанный для создания гидрологически корректных цифровых моделей рельефа (ЦМР). Он основан на программе ANUDEM, разработанной Майклом Хатчинсоном (Michael Hutchinson – 1988, 1989, 1996, 2000, 2011). Литература о применении программы ANUDEM для построения ЦМР на территорию континентов: Hutchinson and Dowling (1991); ANU Fenner School of Environment and Society and Geoscience Australia (2008) Применения ЦМР для моделирования окружающей среды – Hutchinson and Gallant (2000); Hutchinson (2008). Дальнейшее развитие ANUDEM – Hutchinson et al. (2009, 2011). Текущая версия ANUDEM, которая используется в ArcGIS, – 5.3.
Инструмент Топо в растр (Topo to Raster) предназначен для получения растра путем интерполяции значений высот при введении ограничений, которые обеспечивают:
- Связанную дренажную структуру
- Корректное представление хребтов и водотоков из входных данных изолиний
Как таковой, он является просто интерполятором ArcGIS, специально разработанным для разумной работы с входными данными изолиний.
Инструмент Топо в растр по параметрам (Topo to Raster by File) удобен для неоднократного запуска инструмента Топо в растр (Topo to Raster), т. к. обычно легче поменять одну запись в файле параметров и перезапустить инструмент, чем каждый раз заполнять диалоговое окно инструмента.
Процесс интерполяции
Процедура интерполяции была разработана таким образом, чтобы в максимальной степени воспользоваться преимуществами обычно используемых входных данных, а также известными характеристиками поверхностей высот. Этот метод использует методику интерполяции, известную как итеративная конечная разница. Она оптимизирована таким образом, чтобы сочетать в себе вычислительную эффективность локальных методов интерполяции, например, интерполяции по методу обратно взвешенных расстояний (ОВР), без потери непрерывности поверхности, обеспечиваемой глобальными методами интерполяции, такими как Кригинг и Сплайн. По сути, это дискретизированная методика плоского сплайна (thin plate spline) (Wahba, 1990), для которого был изменен параметр поправки за шероховатость (roughness penalty), что позволяет подобранной ЦМР следовать резким изменениям в земной поверхности (таким как ущелья, горные хребты, крутые обрывы).
Вода является первичной эрозионной силой, определяющей общую форму многих ландшафтов. По этой причине большинство ландшафтов имеют много вершин холмов (локальных максимумов) и несколько локальных понижений (локальных минимумов), что приводит к связанной дренажной модели. Инструмент Топо в растр (Topo to Raster) использует эти знания о поверхностях и накладывает ограничения на процесс интерполяции, что приводит к связанной структуре речного бассейна и корректному представлению водотоков и горных хребтов. Это введенное дренажное условие производит поверхности более высокой точности с меньшим количеством входных данных. Количество входных данных может быть на порядок меньше, чем обычно требуется для адекватного описания поверхности с оцифрованными изолиниями, что в дальнейшем минимизирует затраты на получение достоверных ЦМР. Общее дренажное условие также практически исключает необходимость редактирования или повторной обработки для удаления ложных локальных понижений в созданной поверхности.
При удалении ложных понижений программа действует консервативно, она не будет накладывать условие дренажа на те местоположения, где существуют противоречия с входными данными по высотам. Такие местоположения обычно возникают в файле диагностики как локальные понижения. Используйте эту информацию для исправления ошибок данных, особенно при обработке больших наборов данных.
Процесс принудительного заполнения
Целью процесса принудительного заполнения является удаление всех точек локальных понижений в выходной ЦМР, которые не были идентифицированы как локальные понижения в выходном наборе классов объектов локальных понижений. Программой предполагается, что все не идентифицированные локальные понижения – ошибки, т.к. локальные понижения обычно редко встречаются в естественных ландшафтах (Goodchild and Mark, 1987).
Алгоритм принудительного заполнения предпринимает попытку очистить ложные локальные понижения путем изменения ЦМР, проводя дренажные линии или линии стока через самую низкую точку седловины в области стока, окружающей каждое ложное локальное понижение. Он не стремится удалить реальные локальные понижения, как это делала функция Локальное понижение (Sink). Поскольку очистка локальных понижений зависит от допуска в высотах, при попытках удалить ложные понижения программа ведет себя консервативно. Другими словами, она не удаляет ложные локальные понижения, которые противоречат входным данным высот больше, чем на значение Допуска 1.
В качестве дополнительных данных при формировании принудительного стока могут быть также использованы данные линий водотоков. Это полезно, если требуется более точное размещение водотоков. Распределение потоков моделируется тем, что каждая ячейка может иметь до двух направлений стока.
Опция формирования принудительного стока может быть отключена; в этом случае процесс удаления локальных понижений будет проигнорирован. Эта возможность может быть полезной в тех случаях, когда у вас есть данные по изолиниям для иного, чем высоты, показателя (например, температуры), и вы хотите построить поверхность по этим данным.
Использование данных изолиний
Изолинии были первоначально наиболее распространённым методом хранения и представления информации о высотах. К сожалению, этот метод также является наиболее трудным для правильного использования общими методами интерполяции. Этот недостаток определяется отсутствием информации между изолиниями, особенно в областях с понижениями рельефа.
В начале процесса интерполяции инструмент Топо в растр использует информацию, заданную изолиниями, для построения первичной генерализованной дренажной модели. Для этого определяются точки локального максимума кривизны каждой изолинии. Далее при помощи изначального растра высот рассчитывается сеть криволинейных водотоков и горных хребтов, проходящих через эти точки (Hutchinson, 1988). Положение этих линий итеративно корректируется по мере обновления высот ЦМР. Эта информация используется для обеспечения надлежащих гидрогеоморфологических свойств выходной ЦМР, а также может использоваться для проверки точности выходной ЦМР.
Точки изолиний также используются при интерполяции значений высот для каждой ячейки. Вся информация по изолиниям считывается и обощается. В пределах каждой ячейки с изолиний может быть считано до 100 точек, при этом уникальное значение высоты ячейки, через которую проходит изолиния, равно среднему значению высот этих точек. При каждом разрешении ЦМР для каждой ячейки используется только одна критическая точка. По этой причине, данные с такой плотностью горизонталей, когда выходную ячейку пересекает несколько изолиний, являются избыточными.
После определения общей морфологии поверхности данные изолиний также используются в интерполяции значений высот в каждой ячейке.
Если данные изолиний используются для интерполяции информации высот, все данные изолиний считываются и обощаются. Из этих изолиний в пределах каждой ячейки считываются максимум 50 точек данных. В конечном разрешении только одна критическая точка используется для каждой ячейки. По этой причине, данные с такой плотностью горизонталей, когда выходную ячейку пересекает несколько изолиний, являются избыточными.
Использование данных озер
В ранних версиях реализации Топо в растр полигоны-озера были простыми масками, принимающими за высоту поверхности каждого озера наименьшее значение высоты из всех значений ЦМР, непосредственно окружающих озеро. Алгоритм границ озер был обновлен, и сейчас он включает автоматическое определение высоты озера, в полном согласии со связанными линиями водотоков и соседними значениями высот.
Новый метод границы водоема все так же принимает каждую береговую линию за изолинию с неизвестной высотой и итеративно определяет высоту этой изолинии исходя из значений ячеек вдоль нее. В то же время при расчете высоты береговой линии каждого озера учитываются высоты других озер, расположенных вверх и вниз по течению. Высота границы каждого водоема также приводится в соответствие с соседними значениями ЦМР. В результате работы алгоритма значения ячеек, расположенных сразу за пределами озера, превышают высоту границы озера, тогда как значения ячеек в пределах озера являются меньше значений границы.
Границы озер могут включать в себя острова, а также озера внутри таких островов. Все значения ЦМР в пределах озера, задаваемыми полигонами границ водоема, принимают значение высоты, рассчитанное для границы озера.
Использование данных обрывов
Линии обрывов и скал служат для обозначения полного разрыва непрерывности значений ячеек, прилегающих с обеих сторон к линии обрыва, при кодировке этих значений в выходной растр. Линии обрывов должны иметь направление, при этом низкая сторона обрыва находится слева от линии, а высокая сторона – справа.. Это дает возможность удаления точек данных высот, расположенных на несоответствующей стороне от обрыва, при кодировке растра, а также более аккуратного размещения обрывов относительно линий водотоков.
Было выявлено, что незначительные сдвиги в положении водотоков и обрывов при их объединении в растре могут привести к появлению фиктивных взаимодействий между ними. Как следствие, был разработан автоматический метод для внесения небольших корректировок в положение как линий водотоков, так и линий обрывов, для того чтобы свести эти фиктивные отношения к минимуму.
Использование данных береговой линии
Ячейки конечной выходной ЦМР, лежащие вне полигонов, указанных в данном классе полигональных объектов, принимают внутренне заданное определенное значение, которое меньше минимального предела высоты, указанного пользователем. В результате, если в состав входных данных включить полные очертания береговой линии, они будут автоматически обрезаны под обрабатываемый экстент.
Интерполяция для нескольких разрешений
Программа использует метод интерполяции с несколькими разрешениями; выполнение алгоритма начинается с растра с более грубым разрешением, а затем, постепенно, алгоритм переходит к более высокому, заданному пользователем, разрешению. При каждом разрешении происходит принудительное формирование дренажных условий, выполняется интерполяция, а количество оставшихся локальных понижений записывается в выходной файл диагностики.
Обработка данных водотоков
Для инструмента Топо в растр требуется, чтобы дуги в сети водотоков были направлены вниз по склону, и чтобы в сети не было полигонов (озер).
Данные водотоков должны иметь вид древовидной системы, состоящей из одиночных дуг, при этом все параллельные берега рек, полигоны озер и т.д. должны быть удалены путем интерактивного редактирования. При редактировании полигонов водоемов, не входящих в дренажную сеть, формируемая дуга, которая огибает соответствующую область, должна начинаться и заканчиваться в одной точке. Дуга должна следовать по пути исторического русла, если оно известно и существует. Если известен рельеф озера, полигон озера и его рельеф могут использоваться в качестве входных данных типа CONTOUR.
Чтобы отобразить направление линейных секторов, измените символы на Символы со стрелкой на конце. Это позволит показать отдельные сегменты линий, при этом стрелка на конце отрезка будет указывать на направление линии.
Создание смежных растров и мозаики из них
Иногда необходимо создать ЦМР из смежных листов входных данных. Обычно это происходит, когда входные объекты получены из ряда листов карты или когда, из-за ограничений памяти, входные данные должны быть обработаны в несколько частей.
Процесс интерполяции использует входные данные из прилегающих областей для определения морфологии поверхности и дренажного стока и интерполирует выходные значения. Однако значения ячеек на границах любой выходной ЦМР не так надёжны, как в центральной области, т.к. они интерполируются с половиной информации.
Чтобы сделать наиболее точные прогнозы на краях исследуемой области, экстент входных наборов данных должен быть больше исследуемой области. Параметр Область расширения (в ячейках) (Margin in cells) отвечает за метод обрезки краев выходных ЦМР на основе заданного пользователем расстояния. Края перекрывающихся областей должны быть шире, по крайней мере на 20 ячеек.
В смежных областях должно быть небольшое перекрытие входных данных, если несколько выходных ЦМР будут комбинированы в один растр. Без этого перекрытия края объединённых ЦМР не могут быть сглажены. Экстенты входных наборов данных из каждой интерполяции должны иметь даже большую площадь, чем если бы предстояло выполнить единственную интерполяцию для одного набора данных, чтобы таким образом удостовериться, что края можно расчитать как можно точнее.
После построения нескольких ЦМР их можно объединить, используя инструмент геообработки Мозаика (Mosaic) с опциями BLEND или MEAN. Эта функция предоставляет опции для обработки перекрывающихся областей, чтобы сгладить переход между наборами данных.
Оценка выходных данных
Каждая созданная поверхность должна быть оценена для обеспечения того, что данные и параметры, поставляемые в программе, приводят к реалистичному представлению поверхности. Существует много способов оценки качества выходной поверхности в зависимости от типа входных данных, доступных для создания поверхности.
Стандартный метод оценки – создание изолиний по новой поверхности с помощью инструмента Изолинии (Contour) и сравнение их с входными изолиниями. Для сравнения результатов лучше строить эти новые изолинии с половинным интервалом исходного интервала изолиний. Одновременное отображение исходных изолиний и вновь созданных изолиний может помочь в выявлении ошибок интерполяции.
Другой метод визуального сравнения – сравнение дополнительного выходного дренажного покрытия с известными водотоками и хребтами. Класс объектов понижений содержит водотоки и хребты, созданные программой в процессе принудительного заполнения понижений. Эти водотоки и хребты должны совпадать с известными потоками и хребтами на местности. Если в качестве входных данных был использован класс водотоков, выходные водотоки должны почти полностью совпадать с входными, хотя они могут быть немного более генерализованными.
Общий метод оценки качества созданной поверхности – исключение некоторой части входных данных из процесса интерполяции. После построения поверхности высоты в этих известных точках могут быть вычтены из построенной поверхности, и может быть оценено, насколько хорошо новая поверхность представляет истинную поверхность. Эти различия могут использоваться для вычисления измерения ошибки для поверхности, например, среднеквадратичной ошибки (RMS).
Инструмент Топо в растр (Topo to Raster) имеет полный набор процедур для оценки качества корректированной ЦМР, для оптимизации разрешения ЦМР и для обнаружения ошибок во входных данных.
Дополнительный Выходной файл диагностики (Output diagnostic file) может использоваться для оценки того, насколько эффективны заданные параметры допуска при удалении локальных понижений во входных данных. Уменьшение значений допусков может заставить программу вести себя более консервативно при удалении локальных понижений.
Класс объектов Остаточные точки локальных понижений (Output remaining sink point feature) содержит местоположения всех остаточных ложных локальных понижений. Его следует использовать в сочетании с выходными объектами полилиний водотоков для обнаружения ошибок во всех входных топографических данных.
Класс выходных точечных объектов невязок содержит местоположения всех больших невязок данных высот в масштабе локальной ошибки дискретизации. Большие масштабные невязки указывают на конфликты между входными данными высот и данными линий водотоков. Они также могут быть связаны с ошибками, выявленными при автоматическом принудительном заполнении речного бассейна. Эти конфликты могут быть устранены при помощи дополнительной линии водотока и/или дополнительных значений высот точек, но сначала следует проверить и исправить ошибки в существующих входных данных. Большие немасштабные невязки обычно указывают на ошибки в высотах входных данных.
Класс выходных точечных объектов ошибок изолиний содержит местоположения точек, лежащих на входных изолиниях, со значительными невязками относительно корректированной ЦМР. Значение ошибки, равное 1, чаще всего указывает на местоположения точек, в которых соединяются изолинии с различной высотой – это верный индикатор ошибки в подписях изолиний.
Класс выходных точечных объектов ошибок в водотоках и обрывах служит важным показателем качества линий водотоков и обрывов, в частности, он указывает на ошибки в их направлениях, и этот класс всегда следует подвергать анализу.
Данный класс объектов имеет следующие коды:
1. Настоящий контур в сети данных водотока.
2. Контур в сети водотока, закодированный в выходном растре.
3. Контур в сети водотока через связующие озера.
4. Точка рукавов.
5. Водоток выше скалы (водопад).
6. Точки, указывающие на множество точек выхода водотока из озер.
7. Код не используется.
8. Точки около скал, чья высота не соответствует направлению скал.
9. Код не используется.
10. Кольцевой рукав удален.
11. Рукав без впадающего водотока.
12. Растеризованный рукав в выходной ячейке, отличной от той, в которой разделилась линия водотока.
13. Дополнительные условия ошибок обработки – индикатор очень сложных данных водотока.
Класс выходных полилинейных объектов водотоков содержит все ограничения речного бассейна, налагаемые инструментом Топо в растр, которые определяются входными линиями водотоков, линиями водотоков и обрывов, построенными на основе изолиний, а также линиями водотоков, определенными при автоматическом принудительном заполнении речного бассейна. Эта информация может использоваться для проверки ошибок в местоположении входных линий водотоков, а также для проверки данных на соответствие ограничениям, связанным со входными линиями водотоков и с автоматическим принудительным заполнением бассейна. Каждому типу производной линии водотока присвоен свой код. Линии водотока, перекрещивающиеся с линиями обрывов, обозначаются короткими по длине линиями водотока, каждая ячейка имеет отдельный код. Данный класс объектов также включает в себя линии, помечающие большие области входных данных рельефа, превышающие второй допуск по высоте, посредством связанных линий водотоков и озер. Эта информация является удобным показателем ошибок в исходных данных.
Полилинейные объекты имеют следующие коды:
1. Входная линия водотока не выше скалы.
2. Входная линия водотока выше скалы (водопад).
3. Принудительное заполнение, очищающее случайный сток.
4. Линия водотока, определенная по углу изолинии.
5. Линия гребня, определенная по углу изолинии.
6. Код не используется.
7. Дополнительные условия данных линии водотока.
8. Код не используется.
9. Линия, указывающая на большой промежуток в данных высоты.
Смещение изолиний
В алгоритме интерполяции существует небольшое смещение, которое приводит к тому, что входные изолинии оказывают более сильное влияние на выходную поверхность непосредственно в месте изолинии. Это смещение может привести к небольшому сглаживанию выходной поверхности, т.к. она пересекает изолинию. Это может привести к ошибочным результатам при вычислении кривизны профиля выходного растра, а в остальном это не заметно.
Вероятные причины проблем с инструментом Topo в растр
Ниже приведены разъяснения и методы решения наиболее распространенных проблем, возникающих при работе инструмента Топо в растр.
- Не хватает системных ресурсов. Алгоритм инструмента Топо в растр стремится в процессе обработки хранить как можно больше информации в оперативной памяти. Это обеспечивает одновременный доступ к данным точек, изолиний, локальных понижений, водотоков и озёр. Чтобы облегчить обработку больших наборов данных, рекомендуется закрыть ненужные приложения до запуска инструмента, чтобы освободить физическую RAM. Также важно иметь достаточно свободного места на диске.
- Плотность изолиний или точек может быть избыточной для заданного размера ячейки. Если одна выходная ячейка покрывает несколько входных изолиний или точек, алгоритм, возможно, не сможет установить значение для этой ячейки. Чтобы решить эту проблему, попробуйте выполнить один из следующих шагов:
- Уменьшите размер ячейки, затем, после работы инструмента Топо в растр, произведите пересчет растра обратно к большему размеру ячейки.
- Переведите входные данные в растр по частям, используя параметры Выходной экстент и Область расширения (в ячейках). Полученные составляющие части растра соедините при помощи инструмента Мозаика (Mosaic).
- Разрежьте входные данные на секции с перекрытиями и запустите инструмент Топо в растр отдельно для каждой секции. Полученные составляющие части растра соедините при помощи инструмента Мозаика (Mosaic).
- Применение метода интерполяции поверхности может не соотноситься с входным набором данных. Например, если есть входные данные локальных понижений с большим числом точек, чем было бы ячеек в выходном растре, инструмент не будет работать. Источники данных, имеющие высокую плотность точек (например, данные лазерного сканирования), могут привести к сходным проблемам. В этом случае может помочь использование опции NO_ENFORCE , но для предотвращения неправильного применения важно правильно понимать работу методов интерполяции.
Литература
ANU Fenner School of Environment and Society and Geoscience Australia, 2008. GEODATA 9 Second DEM and D8 Digital Elevation Model and Flow Direction Grid, User Guide. Geoscience Australia, 43 pp. Веб-сайт: http://www.ga.gov.au/image_cache/GA11644.pdf.
Goodchild, M. F., and D. M. Mark. 1987. The fractal nature of geographic phenomena. Annals of Association of American Geographers. 77 (2): 265–278.
Hutchinson, M. F. 1988. Calculation of hydrologically sound digital elevation models. Статья представлена на Третьем международном симпозиуме по работе с пространственными данными (Third international Symposium on Spatial Data Handling) в г. Сидней, Австралия.
Hutchinson, M. F. 1989. A new procedure for gridding elevation and stream line data with automatic removal of spurious pits. Journal of Hydrology, 106: 211–232.
Hutchinson, M. F., and T. I. Dowling. 1991. A continental hydrological assessment of a new grid-based digital elevation model of Australia. Hydrological Processes 5: 45–58.
Hutchinson, M. F. 1993. Development of a continent-wide DEM with applications to terrain and climate analysis. In Environmental Modeling with GIS, ed. M. F. Goodchild et al., 392–399. New York: Oxford University Press.
Hutchinson, M. F. 1996. A locally adaptive approach to the interpolation of digital elevation models. In Proceedings, Third international Conference/Workshop on integrating GIS and Environmental Modeling. Santa Barbara, CA: National Center for Geographic Information and Analysis. Веб-сайт: http://www.ncgia.ucsb.edu/conf/SANTA_FE_CD-ROM/sf_papers/hutchinson_michael_dem/local.html.
Hutchinson, M.F. 2000. Optimising the degree of data smoothing for locally adaptive finite element bivariate smoothing splines. ANZIAM Journal 42(E): C774–C796.
Hutchinson, M.F. and Gallant, J.C. 2000. Digital elevation models and representation of terrain shape. In: J.P. Wilson and J.C. Gallant (eds) Terrain Analysis. Wiley, New York, pp. 29–50.
Hutchinson, M.F. 2008. Adding the Z-dimension. In: J.P. Wilson and A.S. Fotheringham (eds), Handbook of Geographic Information Science, Blackwell, pp 144–168.
Hutchinson, M.F., Stein, J.A., Stein, J.L. and Xu, T. 2009. Locally adaptive gridding of noisy high resolution topographic data. In Anderssen, R.S., R.D. Braddock and L.T.H. Newham (eds) 18th World IMACS Congress. Modelling and Simulation Society of Australia and New Zealand and international Association for Mathematics and Computers in Simulation, July 2009, pp. 2493–2499. Веб-сайт: http://www.mssanz.org.au/modsim09/F13/hutchinson.pdf.
Hutchinson, M.F., Xu, T. and Stein, J.A. 2011. Recent Progress in the ANUDEM Elevation Gridding Procedure. In: Geomorphometry 2011, edited by T. Hengel, I.S. Evans, J.P. Wilson and M. Gould, pp. 19–22. Redlands, California, USA. Веб-сайт: http://geomorphometry.org/HutchinsonXu2011.
Wahba, G. 1990. Spline models for Observational data. Статья представлена на очередной региональной конференции CBMS-NSF по прикладной математике. Philadelphia: Soc. Ind. Appl. Maths.