v.surf.bsplinegrass - Интернет в облаке

ПРОГРАММА:

ИМЯ

v.surf.bspline - Выполняет бикубическую или билинейную сплайн-интерполяцию с Тихоновым
регуляризация.

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

вектор, поверхность, интерполяция, лидар

СИНТАКСИС

v.surf.bspline
v.surf.bspline --Помогите
v.surf.bspline [-ce] вход=имя [слой=string] [обзор=имя] [разреженный_вход=имя]
[выходной=имя] [растр_выход=имя] [маска=имя] [ew_step=плавать] [нс_степ=плавать]
[метод=string] [лямбда_i=плавать] [решающее устройство=имя] [максит=целое] [ошибка=плавать]
[Память=целое] [-затирать] [-помощь] [-подробный] [-тихий] [-ui]

Флаги:
-c
Найдите лучший регуляризирующий параметр Тихонова с помощью перекрестной проверки «исключение одного»
метод

-e
Оцените плотность точек и расстояние
Оцените плотность точек и расстояние до точек входного вектора в текущем
регион расширяется и уходит

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

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

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

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

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

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

слой=string
Номер или название слоя
Векторные объекты могут иметь значения категорий в разных слоях. Это число определяет
какой слой использовать. При использовании с прямым доступом к OGR это имя слоя.
По умолчанию: 1

обзор=имя
Имя столбца атрибутов со значениями, которые будут использоваться для аппроксимации
Если не указано и вводится трехмерная векторная карта, используются z-координаты.

разреженный_вход=имя
Имя входной векторной карты с разреженными точками
Или источник данных для прямого доступа к OGR

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

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

маска=имя
Растровая карта, используемая для маскирования (применяется только к выходным растровым изображениям)
Только ячейки, которые не равны NULL и не равны нулю, интерполируются

ew_step=плавать
Длина каждого шага сплайна в направлении восток-запад
По умолчанию: 4

нс_степ=плавать
Длина каждого шага сплайна в направлении север-юг
По умолчанию: 4

метод=string
Алгоритм сплайновой интерполяции
Опции: билинейный, бикубический
По умолчанию: билинейной
билинейной: Билинейная интерполяция
бикубический: Бикубическая интерполяция

лямбда_i=плавать
Параметр регуляризации Тихонова (влияет на сглаживание)
По умолчанию: 0.01

решающее устройство=имя
Тип решателя, который должен решать систему симметричных линейных уравнений
Опции: холецкий, cg
По умолчанию: Cholesky

максит=целое
Максимальное количество итераций, используемых для решения системы линейных уравнений
По умолчанию: 10000

ошибка=плавать
Критерии прерывания ошибки для итеративного решателя
По умолчанию: 0.000001

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

ОПИСАНИЕ

v.surf.bspline выполняет билинейную / бикубическую сплайн-интерполяцию с Тихоновым
регуляризация. В вход 2D или 3D вектор пунктов карта. Значения для интерполяции могут быть
значения z трехмерных точек или значения в столбце заданного пользователем атрибута в 3D или 2D
векторная карта. Вывод может быть растровым (растр_выход) или вектор (выходной) карта. По желанию
Можно ввести векторную карту "разреженных точек", которая указывает местоположение выходной вектор
пунктов.

ПРИМЕЧАНИЯ

С теоретической точки зрения процедура интерполяции состоит из двух частей:
во-первых, оценка линейных коэффициентов сплайн-функции выводится из
точки наблюдения с использованием регрессии наименьших квадратов; второй - вычисление
интерполированная поверхность (или интерполированные векторные точки). Здесь используются двухмерные шлицы.
кусочно-ненулевые полиномиальные функции, вычисляемые в ограниченной двумерной области. Длина
шага каждого сплайна определяется ew_step для направления восток-запад и нс_степ для
направление север-юг. Для оптимальной работы длина шага шлицев не должна быть
меньше, чем расстояние между точками наблюдения. Наблюдение за каждой векторной точкой
моделируется как линейная функция ненулевых сплайнов в области наблюдения.
Регрессия наименьших квадратов предсказывает коэффициенты этих линейных функций.
Регуляризация позволяет избежать необходимости иметь одно наблюдение и один коэффициент для каждого
сплайн (во избежание нестабильности).

При равномерно распределенных точках данных шаг сплайна, соответствующий максимальному
расстояние между двумя точками как в восточном, так и в северном направлениях достаточно. Но часто
точки данных не распределяются регулярно и требуют статистической регуляризации или
предварительный расчет. В таких случаях v.surf.bspline попытается минимизировать градиент
билинейные сплайны или кривизна бикубических сплайнов в областях, где отсутствуют точечные наблюдения.
Как правило, длина шага шлица должна быть больше среднего расстояния между
точки наблюдения (удвоенное расстояние между точками - хорошая отправная точка). Отдельный
аргументы длины шага сплайна восток-запад и север-юг позволяют пользователю учитывать некоторые
степень анизотропии в распределении точек наблюдения. Короткая длина шага шлицевого соединения
- особенно длины шага шлицев, которые меньше расстояния между наблюдениями
баллов - может значительно увеличить время обработки.

Более того, максимальное количество шлицев для каждого направления в каждый момент времени фиксировано,
независимо от длины шага шлица. По мере увеличения общего количества используемых шлицев (т. Е.
с малой длиной шага сплайна), область автоматически разбивается на подобласти для
интерполяция. Каждая подобласть может содержать не более 150x150 шлицев. Чтобы избежать субрегиона
граничные задачи, подобласти создаются так, чтобы частично перекрывать друг друга. Средневзвешенное значение
наблюдений, основанных на точках, рассчитывается в пределах каждой подобласти.

Параметр регуляризации Тихонова (лямбда_i) действует для сглаживания интерполяции. С
небольшой лямбда_i, интерполированная поверхность точно следует за точками наблюдения; больший
value приведет к более плавной интерполяции.

Входными данными может быть двухмерная или трехмерная векторная карта точек. Если вход 2D и обзор не дано чем
z-координаты используются для интерполяции. Параметр обзор требуется при вводе 2D
векторная карта.

v.surf.bspline может произвести растр_выход ИЛИ выходной (но НЕ одновременно). Обратите внимание, что
топология не строится для выходной векторной точечной карты. При необходимости можно построить топологию.
by v.сборка.

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

Для определения оптимального
лямбда_i значение, которое дает интерполяцию, которая наилучшим образом соответствует исходному наблюдению
данные. Чем больше точек используется для перекрестной проверки, тем больше времени требуется для
вычисление. Эмпирическое тестирование показывает, что максимальное количество баллов составляет 100 баллов.
рекомендуемые. Обратите внимание, что перекрестная проверка может выполняться очень медленно, если более 100 наблюдений
используются. Выходные отчеты перекрестной проверки значить и RMS остатков от истинного
значение точки и оценка по интерполяции для фиксированного ряда лямбда_i
ценности. Если выбрана перекрестная проверка, выходной векторный или растровый файл создаваться не будет.

ПРИМЕРЫ

Базовый интерполяция
v.surf.bspline input = point_vector output = interpolate_surface method = bicubic
Будет сделана интерполяция бикубическим сплайном и карта векторных точек с оценкой (т. Е.
интерполированные) значения будут созданы.

Базовый интерполяция и растр выходной a дольше сплайн шаг
v.surf.bspline input = point_vector raster = interpolate_surface ew_step = 25 ns_step = 25
Билинейная сплайн-интерполяция будет сделана с длиной шага сплайна 25 единиц карты. An
Будет создана интерполированная растровая карта с текущим разрешением региона.

оценка of лямбда_i параметр a пересекать Проверка процесс
v.surf.bspline -c input = point_vector

оценка on редкий пунктов
v.surf.bspline input = point_vector sparse = sparse_points output = interpolate_surface
Будет создана выходная карта векторных точек, соответствующая разреженной векторной карте,
с интерполированными значениями.

. атрибут ценности вместо z-координаты
v.surf.bspline input = point_vector raster = interpolate_surface layer = 1 \
column = attrib_column
Интерполяция будет выполняться с использованием значений в attrib_column, в таблице, связанной
со слоем 1.

север каролина расположение пример через z-координаты для интерполяция
g.region region = сельский_1м res = 2 -p
v.surf.bspline input = elev_lid792_bepts raster = elev_lid792_rast \
ew_step = 5 ns_step = 5 method = бикубическая лямбда_i = 0.1

KNOWN ВОПРОСЫ

Известные проблемы:

Во избежание проблем с оперативной памятью необходима вспомогательная таблица для записи некоторых
промежуточные расчеты. Это требует GROUP BY Используется функция SQL, которой нет
поддерживается драйвером DBF. По этой причине выходные данные векторной карты (выходной) не допускается
с драйвером DBF. Нет проблем с выводом растровой карты из драйвера DBF.

Ссылки

· Бровелли М.А., Канната М. и Лонгони Ю.М., 2004, LIDAR Data Filtering и DTM
Интерполяция в GRASS, Транзакции в ГИС, апрель 2004 г., т. 8, вып. 2, стр.
155-174(20), Blackwell Publishing Ltd.

· Бровелли М.А. и Канната М., 2004, Реконструкция цифровой модели местности в городе.
области по данным воздушного лазерного сканирования: метод и пример для Павии
(Северная Италия). Компьютеры и науки о Земле 30, стр. 325-331.

· Brovelli M. A e Longoni UM, 2003, Программное обеспечение для фильтрации данных LIDAR,
Rivista dell'Agenzia del Territorio, n. 3-2003, с. 11-22 (ISSN 1593-2192)

· Антолин Р. и Бровелли М.А., 2007, Фильтрация данных LiDAR с помощью GRASS GIS для
Определение цифровых моделей местности. Труды Jornadas de SIG Libre,
Жирона, Испания. CD ISBN: 978-84-690-3886-9

Используйте v.surf.bsplinegrass в Интернете с помощью сервисов onworks.net