Форумы GIS-Lab.info

Посоветуйте по или модуль к qgis который поможет привязать по точкам вектор SHP.
Или можно ли привязать как с помощью файла привязки растров????.
Все варианты рассмотрю)

Я вижу только плагин qgsAffine. Но для него нужно самому рассчитывать коэффициенты.

dvd_xaker писал(а):Посоветуйте по или модуль к qgis который поможет привязать по точкам вектор SHP.
Или можно ли привязать как с помощью файла привязки растров????.

В GRASS есть модуль GEORECTIFY TOOL, который позволяет привязывать, как растры, так и вектор. Единственно что - в отличии от привязки растров, при привязке вектора можно использовать трансформацию только первого порядка.

0) забыть о "волшебной кнопке", которая все сделает сама.
1) создать в QGIS новый слой, и наполнить его линиями сдвига (Xold,Yold)->(Xnew,Ynew). Сохранить.
2) загрузить этот слой в R и создать модель преобразования координат, например аффинную. Посчитать ошибки. Выбрать хорошую модель.
3) загрузить слой, который нужно корежить, в R, и применить к нему модель.
4) сохранить трансформированный файл.

KolesovDmitry писал(а):

dvd_xaker писал(а):Посоветуйте по или модуль к qgis который поможет привязать по точкам вектор SHP.
Или можно ли привязать как с помощью файла привязки растров????.

В GRASS есть модуль GEORECTIFY TOOL, который позволяет привязывать, как растры, так и вектор. Единственно что - в отличии от привязки растров, при привязке вектора можно использовать трансформацию только первого порядка.

Наверно самый быстрый путь) так как векторов более 400)
Все спасибо.

Было бы хорошо тему эту осветить на странице вопросов))))))) Я думаю не я один столкунлся с этим)

Поясните как добавить и привязать... что то вообще до меню с самой привязкой не могу дойти(((.
К примеру у меня слой называется karta2_4.shp(+к нему shx, dbf,prj, qpg). Что мне надо ввести в меню? данном или как именно, если не сложно обьясните попорядку..... с Grass вообще не работал(((

GRASS по-нормальному работает только со своим форматом данных, поэтому:
1) Скачиваете какой-нибудь готовый набор (лучше всего spearfish), разархивируете его, помещаете в папку на диске (без пробелов и русского языка в названии)
2) Заходите в область spearfish, набор PERMANENT (как это сделать, например, из интерфейса wxGUI, смотрите сюда). Запускаете модуль v.in.ogr с параметром "-е" и указываете ваш SHP-файл --> и т.д.
3) Предыдущим шагом вы создали новую область GRASS с системой координат и проекцией вашего SHP-файла, по пути его туда импортировали
4) Так же импортируете в отдельную область/набор точки, по которым надо привязывать вектор
5) Заходите в область, где точки. Отображаете точки на экране. Запускаете File --> Georectify, указываете область/набор, где ваш вектор, создаёте группу с вектором --> появляются окна для привязки с вектором --> Привязка...

Как-то так. Если будут вопросы, пож-та, сделайте отд. тему на форуме GRASS.

dvd_xaker писал(а): Наверно самый быстрый путь) так как векторов более 400)

тогда тем более R - один раз запрограммировать, и пусть машина работает

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

Вопрос к gamm

У меня есть 8 наборов векторных данных (непривязанных) для разных территорий. В каждом наборе по 3 слоя, которые нужно привязать в единую систему координат. Ваш подход к привязке векторных данных меня оч. заинтересовал ввиду того, что необходимо эти данные привязать в кратчайшие сроки.

Расскажите, пожалуйста, поподробней про модель преобразования в R. Сам я когда-то сталкивался с R, но вплотную не работал.

rhot писал(а):Вопрос к gamm

У меня есть 8 наборов векторных данных (непривязанных) для разных территорий. В каждом наборе по 3 слоя, которые нужно привязать в единую систему координат. Ваш подход к привязке векторных данных меня оч. заинтересовал ввиду того, что необходимо эти данные привязать в кратчайшие сроки.

Расскажите, пожалуйста, поподробней про модель преобразования в R. Сам я когда-то сталкивался с R, но вплотную не работал.

модель преобразования - это просто функция, например X_new=a11+a12*X_old+a13*Y_old, и т.д. - в общем, что вы предполагаете имеет место. Выбираете базовую систему (или "землю"), в итоге у вас получаются либо пары координат "набор-набор", которые связывают два набора (обычно два снимка, так и буду говорить), и "набор-земля", когда для точек на снимке известны фактические координаты. После этого составляете систему уравнений, включающую преобразования, которая все это увязывает, и пишите функционал метода наименьших квадратов (либо квадраты разностей преобразованных координат "набор-набор", либо квадрат разности преобразованных координат "набор-земля"). Начальное приближение делаете, идя последовательно от "земли", особой точности не нужно.

Лучше, если все координаты будут масштабированы вычитанием среднего и делением на СКО, чтобы не иметь проблем с вычислительной дисперсией.

Запускаете для минимизации квадратов сначала метод Нелдера-Мида- optim(,method="Nelder-Mead"), потом "доводите" квазиньютоновским ("BFGS").

"Остальное подогнать по месту и обработать напильником" (с) Анекдот про автомат Калашникова

А вот систему уравнений на основании чего составлять? того нового слоя с линиями сдвига?

Есть ли примеры применения такого интересного способа привязки векторов? Хочется наглядно посмотреть на эту математику в действии. =)

rhot писал(а):А вот систему уравнений на основании чего составлять? того нового слоя с линиями сдвига?

Есть ли примеры применения такого интересного способа привязки векторов? Хочется наглядно посмотреть на эту математику в действии. =)

этот "интересный способ" называется блочная фототриангуляция, и описан книжках по фотограмметрии много десятилетий назад ... система уравнения составляется из равенства преобразованных координат между собой (для связки наборов между собой - "набор-набор" ), либо равенства преобразованных координат известным (с "землей"). Вместо равенства разности координат нулю берется, как обычно, сумма квадратов разностей и минимизируется. Можно с весами, если каким-то координатам веры больше, каким-то меньше.

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

У меня скорее всего будет "набор-земля". Т.е., например, есть ландсатовский снимок и есть 3 оцифрованных слоя на определённую территорию. Хочется для одного слоя указать линии сдвига, а остальные сдвинуть в R. Это ведь возможно? Я имею ввиду, получится ли это быстрее ручной привязки всех трёх слоёв?..

rhot писал(а):Я имею ввиду, получится ли это быстрее ручной привязки всех трёх слоёв?..

Врядли, ведь вам все-так же придется находить реперы и координаты к ним + скриптинг на R. Количество слоев вообще значения не имеет, если они друг с другом увязаны, хоть 100 их будет, хоть 1.