Форумы GIS-Lab.info

Добрый день,

В документации Mapinfo читаем:
To create your own coordinate system using a particular projection, add an entry to the MAPINFOW.PRJ file.
The parameters of a coordinate system are (in this order):
l Coordinate System Name
l Projection Type
l Datum
l Units l Original Longitude
l Original Latitude
l Standard Parallel 1 l Standard Parallel 2
l Azimuth
l Scale Factor
l False Easting l False Northing
l Range

Есть СК Выборга:
"Выборг", 1008, 1001, 7, 27.95, 0, 1, 250000, -11057.626, 7, 0.999477237, 0.031204196, -485909.751, -0.031204196, 0.999477237, -6699985.630

Начальные параметры нормально укладываются: "Выборг", 1008, 1001, 7, 27.95, 0, 1, 250000, -11057.626

6 аффинных параметров в конце строки соответствуют: 0.999477237, 0.031204196, -485909.751, -0.031204196, 0.999477237, -6699985.630

Какому параметру в описании WKT будет соответствовать выделенная жирным и подчеркнутая 7 ? Это единицы измерения для аффинных коэффициентов?

Zorgis писал(а): ↑

14 фев 2019, 15:45

Какому параметру в описании WKT будет соответствовать выделенная жирным и подчеркнутая 7 ? Это единицы измерения для аффинных коэффициентов?

Совершенно верно, имя этого параметра "Affine units".

спасибо.

Игорь Белов писал(а): ↑

14 фев 2019, 16:04

Совершенно верно, имя этого параметра "Affine units".

не подскажите, есть в интернетах полная спецификация формата Mapinfow.proj со всеми возможными параметрами?
и еще вопрос: возможен вариант пользовательской ск 9999 c одновременным добавлением аффинных параметров (на месте 9999 должно быть 3008 или 1008)?
Делаю программный парсер для этого формата и нужны все возможные варианты написания.

Zorgis писал(а): ↑

14 фев 2019, 16:15

есть в интернетах полная спецификация формата Mapinfow.proj со всеми возможными параметрами?

Всё-превсё найдёте в руководствах пользователя на официальном сайте. Например, версия 12.0, версия 15.0.

Zorgis писал(а): ↑

14 фев 2019, 16:15

возможен вариант пользовательской ск 9999 c одновременным добавлением аффинных параметров (на месте 9999 должно быть 3008 или 1008)?

Разумеется, возможен. Не путайте номер проекции 8 с идентификатором датума, и вопрос отпадёт сам собой.

Делаю программный парсер для этого формата и нужны все возможные варианты написания.

уже всё есть

trir писал(а): ↑

14 фев 2019, 18:46

уже всё есть

пишу для своего проекта, в любом случае нужно разобраться с параметрами. за код спасибо, посмотрю.

Игорь Белов писал(а): ↑

14 фев 2019, 18:19

Всё-превсё найдёте в руководствах пользователя на официальном сайте. Например, версия 12.0, версия 15.0.

руководства смотрел. не нашел состав параметров для пользовательского датума.
правильно понимаю, что указанный в пользовательском датуме Основной Меридиан соответствует PRIMEM в WKT?

это есть в моей статье

Описание датума имеет следующую структуру:
9999, НомерЭллипсоида, dX, dY, dZ, EX, EY, EZ, m, ОсновнойМеридиан
Или упрощённая форма:
999, НомерЭллипсоида, dX, dY, dZ

что указанный в пользовательском датуме Основной Меридиан соответствует PRIMEM в WKT?

чего?
а точно, совсем забыл. там указывается долгота нулевого мередиана

Zorgis писал(а): ↑

14 фев 2019, 19:51

руководства смотрел. не нашел состав параметров для пользовательского датума.

Действительно, эта информация выкинута из документации, и уже давно. Ищите руководства к старым версиям. Здесь найдёте версию 5.0. Смотрите приложение I на странице 501.

Игорь Белов писал(а): ↑

15 фев 2019, 08:25

Ищите руководства к старым версиям. Здесь найдёте версию 5.0. Смотрите приложение I на странице 501.

наконец, конкретика.
Видимо, пользовательские датумы актуальны только для СНГ. Аж 1995 года руководство.

Не совсем понятно, в параметрах toWGS84 для proj.4 параметры с 4 по 6 меняют знак на противоположный. Для формата WKT это тоже актуально или у него знаки как в Mapinfo?

Zorgis писал(а): ↑

15 фев 2019, 15:05

Не совсем понятно, в параметрах toWGS84 для proj.4 параметры с 4 по 6 меняют знак на противоположный.

Если коротко, есть два метода трансформации 7 параметров, Position Vector и Coordinate Frame Rotation. В PROJ.4 принят первый, в MapInfo второй.
Переход от одной системы координат к другой - методы трансформации
Переход от одной системы координат к другой - реализация в программах

Zorgis писал(а): ↑

15 фев 2019, 15:05

Для формата WKT это тоже актуально или у него знаки как в Mapinfo?

WKT 1 или WKT 2? В спецификации WKT 1 о методе за TOWGS84 я не нашёл ни слова. Термин "Бурша-Вольфа" обозначает формулу и относится как к Position Vector, так и к Coordinate Frame Rotation. Ну а GDAL использует параметры PROJ.4 как есть.

В WKT 2 использовать TOWGS84 не рекомендуют. Там всё несколько иначе.

в примере стабильной версии geotools
http://docs.geotools.org/stable/javadoc ... ml#TOWGS84
параметр присутствует. парсер библиотеки проходит.
если он есть в форматах proj.4 и MapInfo, мне кажется, ничто не должно препятствовать использовать TOWGS84 в рамках библиотеки geotools как параметр WKT. Или могут быть подводные камни?
в предыдущих сообщениях на форуме склоняются к мысли, что знаки для WKT в TOWGS84 должны быть как в MapInfo.

Zorgis писал(а): ↑

15 фев 2019, 18:26

в предыдущих сообщениях на форуме склоняются к мысли, что знаки для WKT в TOWGS84 должны быть как в MapInfo

Между тем GDAL/OGR конвертирует TOWGS84 из PROJ.4 в WKT и обратно без изменений.
Раз стандарт не определяет тип трансформации жёстко, разработчики софта могут интерпретировать знаки так, как им удобно.

Игорь Белов писал(а): ↑

15 фев 2019, 22:30

Между тем GDAL/OGR конвертирует TOWGS84 из PROJ.4 в WKT и обратно без изменений.
Раз стандарт не определяет тип трансформации жёстко, разработчики софта могут интерпретировать знаки так, как им удобно.

действительно, утилита gdalsrsinfo при конвертации WKT в proj4 не меняет знаки в toWGS84. Получается нужно менять на противоположные знакам MapInfo.