Трансформация MapInfo.prj в WKT и proj

[trir]

нет, я работаю в AutoCAD Civil 3D

[Stasss]

Я показывал на конференции, как это делается

Не подскажете, как найти запись конференции?

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

нужно переписать proj

По ссылке kbevers предлагает бесполезную проекцию:

+proj=affine +A=... +B=... +D=... +E=... +xoff=... +yoff=...

Если же подумать головой, следует вместо этого предложить добавить к списку стандартных параметров проекций параметр дополнительного преобразования:

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

+affine=A,D,B,E,C,F

Не подскажете, как найти запись конференции?

Не знаю. На сайте конференции ссылка на материалы ведёт в никуда. Надеюсь в скором времени переработать старую статью, чтобы соответствовала. Это точно будет лучше, чем видео.

[Ariki]

По ссылке kbevers предлагает бесполезную проекцию:

Начиная с версии 5.0 proj поддерживает цепочки преобразований, так что в предлагаемом виде аффинную трансформацию можно будет комбинировать с проекциями.

https://proj4.org/operations/pipeline.html#pipeline

Хотя в виде дополнительного параметра в большинстве случаев было бы удобнее.

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

Начиная с версии 5.0 proj поддерживает цепочки преобразований, так что в предлагаемом виде аффинную трансформацию можно будет комбинировать с проекциями.

Просто досадно, что концепция +proj кардинально изменилась, превратишись из описания системы координат в произвольную элементарную трансформацию. Не хватило ума придумать какое-нибудь новое служебное слово вроде +op или +transformation.

[Ariki]

Меня в этих нововведениях больше всего смущает, что раньше, когда проекция описывалась декларативно, все имена параметров в её описании были уникальны, и их порядок не имел значения. А теперь это не так, и проекции, описанные через pipeline, будут несовместимы со многими программами, которые полагались на старое поведение (например, хранили параметры в хэш-таблице). В том числе это касается популярных питоновских интерфейсов к GDAL: rasterio и Fiona.

[MaksymLapshyk]

добрый день. Есть параметры описания системы координат в proj:
+proj=tmerc +lat_0=0 +lon_0=36 +k=1 +x_0=130000 +y_0=-5343999.94 +ellps=krass +towgs84=25,-141,-78.5,0,0.35,0.736,0 +units=m +no_defs +type=crs

записал параметры в MAPINFOW.PRJ следующей строкой:
"название", 8, 9999, 3, 25, -141, -78.5, 0, 0.35, 0.736, 0.000, 0, 7, 36, 0, 1, 130000, -5343999.94

в итоге участки улетают метров на десять.
Помогите, что неправильно записал в MAPINFOW.PRJ???

[lam]

С чем сравниваете?

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

Не забывайте менять знаки параметров вращения:
25, -141, -78.5, 0, -0.35, -0.736, 0

[MaksymLapshyk]

С чем сравниваете?

с тем как участки ложатся на снимки ЕSRI
если делаю в СК63 зона Х, 5 зона , то участок ложиться нормально.

проверил улетает примерно метров на 35.

[MaksymLapshyk]

Не забывайте менять знаки параметров вращения:
25, -141, -78.5, 0, -0.35, -0.736, 0

поменял улетело на 200 м

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

если делаю в СК63 зона Х, 5 зона , то участок ложиться нормально

Вот с этого места подробнее. Как в Вашем MAPINFOW.PRJ выглядит "СК63 зона Х, 5 зона"? И как связана с ней безымянная система координат?

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

поменял улетело на 200 м

Т. е. в QGIS на снимки ESRI ложилось нормально, а в MapInfo улетает?

если делаю в СК63 зона Х, 5 зона , то участок ложиться нормально

Как в Вашем MAPINFOW.PRJ выглядит "СК63 зона Х, 5 зона"? И как связана с ней безымянная система координат?

[lam]

с тем как участки ложатся на снимки ЕSRI

участки улетают метров на десять.

Да, может быть такая погрешность

[MaksymLapshyk]

Т. е. в QGIS на снимки ESRI ложилось нормально, а в MapInfo улетает?

да

Как в Вашем MAPINFOW.PRJ выглядит "СК63 зона Х, 5 зона"?

"1963_X3 зона 5", 8, 1001, 7, 35.5, 0, 1, 5300000, -9214.69

И как связана с ней безымянная система координат?

это описание системы локальной системы координат г. Днепр, которую использует местная архитектура для ведения градкадастра