Программа RUSSIA.GSB

Осталось теперь только понять, как этот файл сетки прикрутить в QGIS. Вроде по описаниям в сети должно быть всё просто, но только сами описания эти мутные-мутные.

 

В мапинфо наверное проще.

 

В ютубе можно поискать
ntv2 transformation qgis - YouTube

--- Сообщения объединены, Feb 23, 2023, Оригинальное время сообщения: Feb 23, 2023 ---

Spoiler

Немного не в тему, но для задумчивых про утонченность МСК - состояние реального ГСК-2011 на февраль 2023 года

 

Официальная документация:
Integrate an NTv2-transformation in QGIS

 

Но только я пока не увидел упоминания про 4 СК. Ведь сказано, что данная сетка каким то образом связывает сразу 4 СК: Программа RUSSIA.GSB #13 (хотя мне это очень непонятно). Но факт остаётся фактом. Либо я конкретно чего то не понимаю.

 

2 СК Из 42 в 84.... В кредо в обратном направлении не посчитать (у меня не вышло во всяком случае)

 

То есть стандартные параметры для СК42 + сетка. Ну так то понятнее. А то с 4-мя как то совсем не алё. А ГСК-2011 это тупо по ГОСТу через ПЗ90. Оки-доки.

PS: Только вот насчёт "стандартных параметров СК42" у меня внезапно "вкралось смутное сомнение". Хмм... По логике они вроде как ненужны, нужна только проекция. Но "это" моя логика, а не автора сетки.

wolodya, ты ещё не пробовал в параметры что то вроде "24.83,-130.97,-81.74,0,0,-0.13,-0.22" прикрутить?

 

Да нет.... Использую какие вбиты.

 

Ну так я навсидку, не вдаваясь в подробности. Значит параметры между Крассовским и WGS-84 всё-таки нужны. Только ежели Крассовский, то почему у тебя везде эллипсоид WGS84 указан? Не принципиально конечно, просто очередная странность.

 

В транскоре вообще то wgs не задается..... Но возможно используется...

 

Честно говоря без понятия, что там задаётся и используется. Но у тебя то вбит (Программа RUSSIA.GSB #28) эллипсоид WGS84. Так и должно быть? Странно это всё. Параметры сдвигов-вращений считай одного эллипсоида, а сам эллипсоид внезапно другой. Странно. Несостыковка какая то между проекциями и эллипсоидом возникает.

 

Как тяжело в этой стране жить, кругом одни кулибины ...

 

Спасибо конечно. Но лингвист из меня только в путь. Не осилить мне по-человечески "этот" момент на буржуйском, без нашего никак.

 

Внесу скромный лепт в вычисления.

Исходный файл a0.txt:

Code:

59d46'07.9482"N 30d19'38.4992"E 0 ГАО_РАН
67d15'57.74658"N 59d19'20.49591"E 0 Контроль_верх
67d10'18.87855"N 59d19'04.72041"E 0 Контроль_низ
57d11'11.0076"N 75d13'37.6140"E 0 Контроль

Считаем из СК-42 в WGS 84 через трансформацию RUSSIA.gsb:

Code:

cs2cs -r -W6 +proj=longlat +ellps=krass +nadgrids=@RUSSIA.gsb,null +to EPSG:4236 a0.txt > b0.txt

Выходной файл b0.txt:

Code:

59d46'07.938563"N 30d19'30.594478"E 0.000 ГАО_РАН
67d15'59.761229"N 59d19'13.167914"E 0.000 Контроль_верх
67d10'20.889600"N 59d18'57.421261"E 0.000 Контроль_низ
57d11'12.974931"N 75d13'34.437637"E 0.000 Контроль

 

А почему у меня данное колдунство не происходит?

Code:

cs2cs Rel. 7.2.1, January 1st, 2021

Файлы a0.txt и b0.txt за исключением табуляций идентичные. RUSSIA.gsb рядом с a0.txt.

Code:

sha1sum RUSSIA.gsb
cd18010c4dd32b95030f2d64464104506a033b87 RUSSIA.gsb

PS: И EPSG:4326.

Отбой. Сам разобрался:

Code:

cs2cs -r -W6 +proj=latlong +ellps=krass +nadgrids=./RUSSIA.gsb,null +to EPSG:4326 a0.txt
59d46'07.938563"N 30d19'30.594478"E 0.000 ГАО_РАН
67d15'59.761229"N 59d19'13.167914"E 0.000 Контроль_верх
67d10'20.889600"N 59d18'57.421261"E 0.000 Контроль_низ
57d11'12.974931"N 75d13'34.437637"E 0.000 Контроль

PS: Интересно, а где по умолчанию (без указания текущей директории "./") эти сетки ищутся?

 

Здесь:

Code:

/usr/share/proj

В M$ Windows здесь:

Code:

C:\OSGeo4W64\share\proj

Ссылки:
Where are PROJ resource files looked for ?
Integrate an NTv2-transformation in QGIS

 

Не получается у меня с буржуйским. До конца не разобрался с "этим" колдунством, но обратная строка такая:

Code:

cs2cs -W6 EPSG:4326 +to +proj=latlong +ellps=krass +nadgrids=./RUSSIA.gsb,null b0.txt
30d19'38.499200"E       59d46'07.948200"N 0.000 ГАО_РАН
59d19'20.495910"E       67d15'57.746580"N 0.000 Контроль_верх
59d19'04.720410"E       67d10'18.878550"N 0.000 Контроль_низ
75d13'37.614000"E       57d11'11.007600"N 0.000 Контроль

Как бы ещё Lat-Long сделать, вместо Long-Lat. Было бы вообще замечательно.

 

-s — меняет порядок координат на выходе (-r меняет порядок координат на входе)

 

Отлично!:

Code:

cs2cs -W6 -s EPSG:4326 +to +proj=longlat +ellps=krass +nadgrids=@./RUSSIA.gsb,null b0.txt
59d46'07.948200"N       30d19'38.499200"E 0.000 ГАО_РАН
67d15'57.746580"N       59d19'20.495910"E 0.000 Контроль_верх
67d10'18.878550"N       59d19'04.720410"E 0.000 Контроль_низ
57d11'11.007600"N       75d13'37.614000"E 0.000 Контроль

 

Что то у меня не получилось в Trimble RTX Джавадовский ринекс загрузить....
Может из за этого вышли большие ошибки по северу ХМАО.
Может кому интересно будет - файлы прилагаю.

https://disk.yandex.ru/d/CKEeuOMhkozmPg