Сложность при привязке растра

[indeec]

Друзья, добрый день!

Возникла сложность при привязке растра, которую никак не могу побороть. Собственно, проблема на скриншоте видна. Точки привязки находятся очень далеко друг от друга и пересекаются.
Я много раз выполнял привязку растра и с таким сталкиваюсь впервые. Конкретно это изображение пробовал привязать полностью с нуля дважды. В точности определения положения на местности уверен.

Привязка.png (80.42 КБ) 8122 просмотра

[Игорь Белов]

Cудя по всему, растр "Градзонирование целиком_modified.tif" уже привязан.

[indeec]

Игорь, спасибо!
Вы имеете в виду, что он уже на карте? Да, совершенно верно. Я открыл свою последнюю попытку. Точки были сохранены в файле.
Растр привязывается, но с большими погрешностями. Добавление точек ничего особо не меняет. Подозреваю, что дело в этой ошибке.

[tikhpetr]

indeec, А выложите Вы его сюда с точками привязки. Я попробую.

[indeec]

indeec, А выложите Вы его сюда с точками привязки. Я попробую.

https://cloud.mail.ru/public/AQsZ/3PYGyDBkG
Спасибо!

[tikhpetr]

Все забрал. Сделаю завтра, т.к., шибко медленно качает SAS. Но полигон по растру лег не прям вот в ЮАР. В смысле на OSM в SAS.

[indeec]

Все забрал. Сделаю завтра, т.к., шибко медленно качает SAS. Но полигон по растру лег не прям вот в ЮАР. В смысле на OSM в SAS.

В вашем сообщении меня пугает аббревиатура ЮАР) Это Московская область, а не Африка)
Тут нужно сказать, наверное, что я новичок в ГИС. Оказалось очень интересная, но и очень объемная по составу знаний область.
Спасибо большое!

[tikhpetr]

Не пугайтесь, это у меня по ночам юмор такой. Я cверил Ваши Tiff-ы со снимком ESRI и картой OSM. Все у Вас нормально привязано. Отлеты минимальные, и связаны они, что тип карты, очевидно отличается от OSM и, тем более от снимка. Я проверял в Global Mapper и ArcMap.
Вот ссылка на ЯДиск. Можете сами позырить.
https://yadi.sk/d/7Go4wOTL_7hWtQ

[indeec]

Не пугайтесь, это у меня по ночам юмор такой. Я cверил Ваши Tiff-ы со снимком ESRI и картой OSM. Все у Вас нормально привязано. Отлеты минимальные, и связаны они, что тип карты, очевидно отличается от OSM и, тем более от снимка. Я проверял в Global Mapper и ArcMap.
Вот ссылка на ЯДиск. Можете сами позырить.
https://yadi.sk/d/7Go4wOTL_7hWtQ

Спасибо большое! То есть, получается, это максимальная точность? Или можно ее еще как-то повысить?

[tikhpetr]

Только инструментальным/документальным - с тахеометром на местности отснять N точек на местности(привязавшись к известным репперам), или получив список координат(достаточно точных) для характерных точек. И то, и другое - ваще морока!
З.Ы. А что, 5-10 м Вас не устраивает?

[indeec]

Только инструментальным/документальным - с тахеометром на местности отснять N точек на местности(привязавшись к известным репперам), или получив список координат(достаточно точных) для характерных точек. И то, и другое - ваще морока!
З.Ы. А что, 5-10 м Вас не устраивает?

Понятно, спасибо! Что касается точности, то в целом устраивает. Но создает некоторые неудобства( Но ничего не поделаешь. Спасибо еще раз!

[tikhpetr]

Понимаете, когда карту привязывают, скажем так, графическим способом, должно быть что-то(объекты), к-рые должны правильно лечь на привязываемую карту. Напр., когда делают карты для навигаторов, критерием является то, чтобы ранее записанные треки и точки легли, как надо. У Вас есть такие критерии? К чему необходимо стремиться?

[gamm]

Или можно ее еще как-то повысить?

1) у вас точность координат (по крайней мере тех, что на первой картинке), порядка 10 метров
2) похоже, что вы пытаетесь привязать карту в градусах, что обычно приводит к искажениям внутри области, даже если точки сядут.

[Игорь Белов]

Ребята, так нельзя!

У вас карта в проекции, результат растеризации без разворота и прочих искажений. Необходимо выяснить, что это за проекция, и привязывать в ней. Выставляем линейное преобразование и перебираем подходящие проекции. В данном случае это оказалась МСК-50 зона 2. Вот правильный файл points:

Код: Выделить всё

mapX,mapY,pixelX,pixelY,enable
2195367.9188821641728282,428592.83977056108415127,265.49787030507690133,-4989.52365787629787519,1
2212402.52802088391035795,433481.52290798351168633,6188.99020412609388586,-3289.95794262701429034,1
2216621.48982810415327549,409760.53861769009381533,7655.53108740053721704,-11537.44991668019065401,1
2199220.88468059338629246,402933.96739552821964025,1605.27916140076240481,-13911.75623692229419248,1

В колонках mapX, mapY координаты в проекции, pixelX, pixelY - координаты в ПИКСЕЛАХ.

Вместо трансформации просим вывести world-файл. Получаем такое:

Код: Выделить всё

2.87591787110083086
0
0
-2.87596360580451194
2194604.22852857038378716
442942.83678602881263942

Теперь можно трансформировать растр в такую проекцию, с которой GeoTIFF хорошо знаком. В данном случае возьмём EPSG:32637 "WGS 84 / UTM zone 37N".

Геопривязанный и трансформированный растр

[tikhpetr]

Игорь Белов, Я, пардон, не понял. Я смотрел, как ложатся "modified tif" на соответствующие снимок и кусок OSM. Окрывал я сначала, ес-сно, tif. GM и Арка сами определили их как WGS. Я и решил, что ТС так их и вязал. Допустим, они изначально были в проекции. Да и хрен с ним. Зачем точки преобразовывать, а потом еще и растр трансформировать? ИМХО, нуно, чтобы ТС выяснил СК исходника, настриг бы потом точек со снимка и привязал по-новой в родной СК. Нет?