Как определить проекцию карты?

denizzz · 2011-03-26T08:19:55.000Z

Есть карта (точнее, план, нарисованный на кальке обводом карты). Известны точные координаты 6 точек (можно и больше) на этой карте.

При преобразовании карты в GeoTIFF и подключении её в merkaartor она неточно ложится на местность, хотя координаты первых трёх точек отлично совпадают.

Я так понимаю, нужно определить проекцию карты и задать её?

(Для привязки и преобразований использую qlandkartegt из Debian)

Zverik · 2011-03-26T10:02:49.000Z

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

denizzz · 2011-03-26T11:27:31.000Z

Что значит привязать по всем точкам?

usm78-gis · 2011-03-26T18:35:52.000Z

denizzz:

Есть карта (точнее, план, нарисованный на кальке обводом карты).
…
(Для привязки и преобразований использую qlandkartegt из Debian)

Каков примерный размер плана (10x10 km , 100x100km )
и можете ли вы примерно определить координаты центра ?
qlandkartegt (насколько я помню) использует систему координат


+proj=merc +ellps=WGS84 +datum=WGS84 +no_defs

также как яндегс, космоснимки и прочая. Для небольшого плана не самая точная, и
сложная для понимания из-за пересчета из “меркаторовских метров”.
Чтобы добиться максимальной точности, имхо, лучше всего задать свою
собственную местную систему координат на базе трансверсального меркатора.
Тут многое будет зависеть от того, внесли ли создатели плана сознательно
нелинейные искажения: в кадастре / картах основанных на ГШ и Местных CK их
нет (и быть не должно), а вот во всяких крупномасштабных планах городов
и т.п. советских времен встречаются сплошь и рядом.

andriano · 2011-03-27T08:27:34.000Z

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

mguzelevich · 2011-03-27T12:14:41.000Z

denizzz:

При преобразовании карты в GeoTIFF и подключении её в merkaartor она неточно ложится на местность, хотя координаты первых трёх точек отлично совпадают.

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

denizzz · 2011-03-27T13:36:13.000Z

Вот этот план:
http://rghost.ru/4953724
(3.5 Мб)

А вот точки для qlandkartegt:
http://rghost.ru/4953740
(там проекции заданы какие-то, но они не работают и нужны другие)

простой текстовый формат, точки можно просмотреть глазами:

-gcp 4019 465 94.6113474 56.1138328
-gcp 813 2943 94.5565941 56.1022981
-gcp 3099 1862 94.5922707 56.1037540
-gcp 2397 293 94.5873587 56.1191711
-gcp 3214 853 94.5975274 56.1121705
-gcp 1094 1761 94.5640222 56.1093609

denizzz · 2011-03-27T13:37:45.000Z

wildMan:

denizzz:

При преобразовании карты в GeoTIFF и подключении её в merkaartor она неточно ложится на местность, хотя координаты первых трёх точек отлично совпадают.

я правильно понял что вы делаете привязку по 6-ти точкам но в итоговом привязаном файле у вас совпали только 3 первых точки?
тогда вы просто неверно выполнили привязку.

Сначала их было 3, но результат выглядел криво. Потом я постепенно добавил ещё 3 точки, думал результат “загнётся” по нужному эллипсоеду но это не помогло.

denizzz · 2011-03-27T13:40:52.000Z

usm78-gis:

denizzz:

Есть карта (точнее, план, нарисованный на кальке обводом карты).
…
(Для привязки и преобразований использую qlandkartegt из Debian)

Каков примерный размер плана (10x10 km , 100x100km )

меньше чем 10 на 10 км

и можете ли вы примерно определить координаты центра ?

да

qlandkartegt (насколько я помню) использует систему координат
+proj=merc +ellps=WGS84 +datum=WGS84 +no_defs

там можно выбирать, вроде бы

также как яндегс, космоснимки и прочая. Для небольшого плана не самая точная, и
сложная для понимания из-за пересчета из “меркаторовских метров”.
Чтобы добиться максимальной точности, имхо, лучше всего задать свою
собственную местную систему координат на базе трансверсального меркатора.
Тут многое будет зависеть от того, внесли ли создатели плана сознательно
нелинейные искажения: в кадастре / картах основанных на ГШ и Местных CK их
нет (и быть не должно), а вот во всяких крупномасштабных планах городов
и т.п. советских времен встречаются сплошь и рядом.

искажений вроде бы нету, хотя глаз - не алмаз )

usm78-gis · 2011-03-27T17:11:06.000Z

denizzz:

искажений вроде бы нету, хотя глаз - не алмаз )

Карта поворачивается где-то на 30° по часовой стрелке,
план составлен в сильно самодельной CK.
Размер пиксела < 2 m, это конечно круто. А гугловские
картинки в этом месте правильно (= с метровой точностью)
привязаны ?

denizzz · 2011-03-27T17:34:25.000Z

usm78-gis:

denizzz:

искажений вроде бы нету, хотя глаз - не алмаз )

Карта поворачивается где-то на 30° по часовой стрелке,
план составлен в сильно самодельной CK.
Размер пиксела < 2 m, это конечно круто. А гугловские
картинки в этом месте правильно (= с метровой точностью)
привязаны ?

да, и бинг тоже

denizzz · 2011-03-29T12:11:21.000Z

Так что, на этом плане можно ставить крест из-за его кривизны?

usm78-gis · 2011-03-30T20:59:34.000Z

denizzz:

Так что, на этом плане можно ставить крест из-за его кривизны?

Почему? Натыкать точек из bing равномерно по всей площади картинки,
и использовать (например) gdalwarp -tps
Предварительно желательно уменьшить число цветов до 16 или даже 8
(например в gimp) и применть nearest neighbour интерполяцию.