Перепроецирование из EPSG: 3785 в WGS-84

[SS_Rebelious]

Ограбив один из ресурсов с помощью SASPlanet, получил растры в проекции +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs (EPSG: 3785).
К нему я привязал растр, но в его свойствах значится проекция +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=clrk66 +towgs84=-8,160,176,0,0,0,0 +units=m +no_defs; хотя во время привязки я указывал "EPSG: 3785".

Проблема в том, что при перепроецировании из обеих этих проекций в WGS-84 растры улетают на 20 километров на север. Интрересно, проблема преобразования связана именно с этими проекциями (встречал в интренете мнения, что использование +towgs84=0,0,0,0,0,0,0 не есть хорошо) или это SASPlanet так подкачал с привязкой (чего раньше не наблюдалось)?

[manisha]

Проблема с исходной проекцией, в том что меркатор на сфере. Например, ердас в версиях до 9-ки (до того, как в параметры проекции добавили выбор сфера/эллипсоид) тоже неправильно (с таким же улетом на север) пересчитывал в проекции на эллипсоиде, даже если задавать новый кастомный эллипосид с равной длины полуосями (для имитации сферы).

[SS_Rebelious]

Спасибо. Буду знать.

[Mavka]

+towgs84=0,0,0,0,0,0,0

Это 7 параметров, применяемых при переходе от одного эллипсоида к другому. В Google-СК используются одновременно 2 эллипсоида - сфера и WGS84, поэтому параметры нулевые. Точнее эллипсоид WGS84, а для проекции заданы другие значения полуосей.

Насколько помню, правильнее было использовать нулевую таблицу поправок "+nadgrids=@null".
Вот собственно так: http://lists.maptools.org/pipermail/pro ... 03122.html

[SS_Rebelious]

Насколько помню, правильнее было использовать нулевую таблицу поправок "+nadgrids=@null".
Вот собственно так: http://lists.maptools.org/pipermail/pro ... 03122.html

Большое спасибо, всё получилось!

P.S. Спасибо bim2010 за обсуждение в irc.

[wlad762]

Странно. В версии QGis 1.6 данная проблема проявлялась явно. Решал с помощью изменения системы координат в Свойствах растрового слоя с EPSG:3785 - Popular Visualisation CRS / Mercator на EPSG:900913 - Google Mercator . Считал это багом СасПланеты. Далее, по снимку делалась привязка плана, делалась векторизация с последующим перепроецированием векторных слоёв в нужное (для нас это WGS84).
Но в версии QGis 1.7 всё решается включением в свойствах проекта опции преобразования координат "на лету".
В скриншоте (555.JPG) растровый слой с СасПланеты, векторные - трек GPS. Всё легло друг на друга с точностью в 4 метра.
Был приятно удивлён перепроецированным растром. Жить стало легче, жить стало веселее.

[Mavka]

wlad762
Перепроецирование растра это, действительно, очень хорошая штука. А на счет СК я немного не понял - 3785 и 900913 это практически одно и тоже...

[wlad762]

Mavka
Перепроецирование растра это, действительно, очень хорошая штука. А на счет СК я немного не понял - 3785 и 900913 это практически одно и тоже...

Понимаю. Но всё же. Баг с расбросом в 20 км пропадал в QGis 1.6 только после вышеописанных действий. А доходить до этого пришлось методом научного тыка. Слава создателям за QGis 1.7 и прочих новых библиотек. Работать стало в разы комфортнее.

Кстати! При использовании растров СасПланеты в MapEdit всё ложилось как надо, правда, насколько помню, для него использовался другой вид привязки.