Как совместить карту и подоснову? Нужной проекции нет

[ali]

Здравствуйте,
посоветуйте, пожалуйста, как решить следующую проблему.
Использую карты в районе Москвы. Так как мне нужны метры, то перевожу исходные .shp файлы
с помощью Mapinfo в формат .mif проекцией Гаусса-Крюгера Пулково 1942 (зона 7)
(у моей программы есть ограничение - работает только с mif-файлами и системой
координат nonEarth).
Но у меня есть электронная векторная подоснова Москвы, которая не сходится с полученными
картами - выглядит чуть повёрнутой и сдвинутой. В организации-изготовителе этой
подосновы мне сказали, что используют "прямоугольную проекцию на эллипсоиде Бесселя".
В моём старом Mapinfo (2006 г) такой нет. Скачал помощь к Mapinfo 10, там есть
эллипсоид Бесселя, но нужной проекции в списке нет.

Есть ли какой-нибудь способ совместить карту и подоснову?
Заранее спасибо за совет.
Александр

[nadiopt]

это у них ЕГСК, она реально на эллипсоиде Бесселя
создайте вручную, эллипсоид Бесселя, а проекция, наверное, Гаусса-Крюгера, а потом экспортните в нормальную Пулково

[ali]

Мне не совсем понятен ответ, видимо, я знаю меньше, чем Вы думаете :)

Что значит, что у них "на эллипсоиде"? Я думал, что на поверхности должны быть градусы (широта/долгота),
а у них метры (хотя не уверен - я уже работаю с MIF-файлом, а что было в самом начале - не знаю).

Что значит "создать вручную эллипсоид Бесселя"? Чем и где создать? Какой файл получится?

Последнее предложение тоже не совсем ясно: Вы имеете в виду просто экспорт
в проекцию Гаусса-Крюгера (Пулково,1942), зона 7? Или слово "нормальная"
значит что-то особое?

С уважением,
Александр

[nadiopt]

грубо говоря, система координат состоит из эллипсоида и проекции
ЕГКО (пардон, не ЕГСК, это я системой семеноводческих станций перепутала) вот:
viewtopic.php?t=5859
http://www.dataplus.ru/ARCREV/Number_16/6_Moscow.htm
вроде бы там система координат на эллипсоиде Бесселя, проекция прямоугольная гаусса-крюгера, а вот где у них ноль - точно не знаю, по слухам, на крыше МИИГАИКа
а Пулково-42 базируется на эллипсоиде Красовского, а проекция там Гаусса-Крюгера
нету мапинфо под рукой, не могу показать, как проекцию создать. для арки вот как:

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

Projection: Gauss_Kruger
False_Easting: 0,000000
False_Northing: 0,000000
Central_Meridian: 0,000000
Scale_Factor: 1,000000
Latitude_Of_Origin: 0,000000
Linear Unit: Meter (1,000000)

Geographic Coordinate System: GCS_Bessel_1841
Angular Unit: Degree (0,017453292519943295)
Prime Meridian: Greenwich (0,000000000000000000)
Datum: D_Bessel_1841
  Spheroid: Bessel_1841
    Semimajor Axis: 6377397,155000000300000000
    Semiminor Axis: 6356078,962818188600000000
    Inverse Flattening: 299,152812799999990000

[ali]

Пока не понял
В приведённых Вами двух ссылках ничего не говорится о проекциях для ЕГКО.
Правильно ли я Вас понял, что слова о том, что проекция "Прямоугольная на эллипсоиде Бесселя"
означают именно проекцию Гаусса-Крюгера, потому что она прямоугольная? Других прямоугольных
не бывает?
Про ноль тоже непонятно: если Гаусса-Крюгера, то надо брать (как мне раньше говорили на этом
форуме) зону 7 для Москвы. Как может быть ноль в другом месте? Или я что-то не то говорю?
"Арки" в Вашем ответе означает ArcView? В его триальной версии можно будет использовать Ваш код?
Я попытался догадаться о смысле строчек кода. А где там что-то относящееся к Москве? И если ноль
("по слухам") - на крыше чего-то в Москве, то как это учитывается?

С уважением,
Александр

[nadiopt]

ух я вас совсем запутала
в ссылках нету про систему координат ЕГКО, я что-то вообще не нашла эти данные в сети. это просто для общего представления
прямоугольная может быть и не гаусса-крюгера, просто она в СССР широко использовалась, вот я и предположила, что в Москве тоже она.
про ноль - это может быть смещение начала координат False_Easting: 0,000000 (на восток)
False_Northing: 0,000000 (на север). а зона скорее всего действительно седьмая.
в арквью код использовать можно, я делала для девятки
завтра спрошу поточнее про ЕГКО, чтоб больше вас не путать

[ali]

> завтра спрошу поточнее про ЕГКО...

Это будет здорово, если получится. Большое Вам спасибо за помощь!
Я звонил в мосгоргеотрест и задавал им вопрос про проекцию, но они
берутся переводить в их проекцию только по договору и, естественно, за деньги, хотя почему бы
им просто не сказать - не понимаю: ведь это не работа, а просто знание.

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

С уважением,
Александр

[Aleksandr Dezhin]

+proj=tmerc +lat_0=55.6666666667 +lon_0=37.5 +x_0=0 +y_0=0 +k_0=1. +a=6377397 +rf=299.15 +towgs84=396,165,557.7,-0.05,0.04,0.01,0 +no_defs

Как сформулировать это в терминах MapInfo я не особо в курсе, но можно конвертнуть слой в шейп, прописать ему эту проекцию Qgis-ом или ogr2ogr, например, и открыть его в MapInfo. По-идее он должен в этом случае сам подхватить правильную СК.

Такая вот СК записанная в формате proj4 более менее подходит для данных ЕГКО, по крайней мере приблизительно, если сравнивать с данными ОСМ то у меня получалась ошибка не больше 10 метров, более точных данных для сравнения у меня не было. lat_0 и lon_0 и параметры эллипсоида нашлись в каких-то документах, towgs84 подбирались из имеющихся...

[Донецков]

Поищите на этом форуме, уже были темы про систему МГГТ, а еще советую посмотреть на сайте ГЕОДЕЗИСТ.RU, там это встречалось неоднократно, даже есть конвертор по пересчету в МГГТ из СК-42 7 зона...
В MapInfo проблема следующая, там нет возможности создать свой эллипсоид... Поэтому приходится подставлять имеющийся эллипсоид (Бесселя 1941), который в последних версиях программы показан как - Эстония 1937 г., но он немного различается с параметрами, которые указаны в ГКИНП, в котором приведены параметры МГГТ...

[ali]

Не могли бы Вы переформулировать ответ, а то я ничего не понял:
- что означает цитата? Это какой-то код? Куда его надо вставлять?
- что значит "конвертнуть слой"? Какой слой? в какой программе?
- что значит "прописать ему эту проекцию Qgis-ом или ogr2ogr"? Это какие-то отдельные программы?
- что значит "он должен в этом случае сам подхватить правильную СК"? что такое СК? ОСМ?

С уважением,
Александр

[ali]

Поищите на этом форуме, уже были темы про систему МГГТ, а еще советую посмотреть на сайте ГЕОДЕЗИСТ.RU, там это встречалось неоднократно, даже есть конвертор по пересчету в МГГТ из СК-42 7 зона...

К сожалению, по моему вопросу ничего не нашёл. Вроде бы кратко мою проблему теперь могу сформулировать таким образом: как перевести систему координат из Пулково 42, зона 7 в координаты Мосгоргеотреста (МГГТ)?
На сайте ГЕОДЕЗИСТ.RU нашёл программу перевода в обратную сторону (из МГГТ в Пулково-42), но даже она не работает под Windows 7.

В MapInfo проблема следующая, там нет возможности создать свой эллипсоид... Поэтому приходится подставлять имеющийся эллипсоид (Бесселя 1941), который в последних версиях программы показан как - Эстония 1937 г., но он немного различается с параметрами, которые указаны в ГКИНП, в котором приведены параметры МГГТ...

Не понял, как это помогает решить мою проблему. Не могли бы Вы пояснить подробнее.

С уважением,
Александр

[nadiopt]

ну к сказанному специалистами в части параметров ЕГКО мне добавить нечего, изучайте матчасть. а зачем вам карта в ЕГКО, сделайте в Пулково

[ali]

ну к сказанному специалистами в части параметров ЕГКО мне добавить нечего, изучайте матчасть.

К сожалению, я не понял, что сказали специалисты
Матчасть изучить готов, с чего надо начать? что почитать о действиях по переводу в нужную проекцию (или её созданию)?

а зачем вам карта в ЕГКО, сделайте в Пулково

У нас подложка от МГГТ. Поэтому подгоняем всё под неё. Можно и перевести МГГТ в Пулково, но на данный момент не понял, как это проделать.

С уважением,
Александр

[nadiopt]

На сайте ГЕОДЕЗИСТ.RU нашёл программу перевода в обратную сторону (из МГГТ в Пулково-42), но даже она не работает под Windows 7.

вот и делайте. попробуйте режим совместимости с windows XP или комп с установленной XP

[Донецков]

Странно - у меня программа работала под 7 и под Vist-ой (возможно проблема - запуск от имени администратора)

В MapInfo в файле PRJ добавляете строку:
"Москва (Бесселя 1941)", 8, 122, 7, 37.5, 55.66666666666666, 1, 0, 0

для пересчета нужно учесть:
план схема - > проекция GK-7 < - > МГГТ - > план схема

как переводить из план схему в проекцию и обратно неоднократно обсуждалось на этом форуме...

Если используете QGIS, то параметры уже были приведены:
+proj=tmerc +lat_0=55.6666666667 +lon_0=37.5 +x_0=0 +y_0=0 +k_0=1. +a=6377397 +rf=299.15 +towgs84=396,165,557.7,-0.05,0.04,0.01,0 +no_defs

Под слоем можно понимать таблицу (в терминах MapInfo)....

OSM - набор данных http://gis-lab.info/projects/osm_shp.html, представлен в десятичных градусах WGS84