SASGIS - SAS.Планета
View Issue Details
0003179SAS.ПланетаРефакторингpublic18-02-2017 15:2218-02-2017 20:38
RedRat 
zed 
normalminoralways
resolvedfixed 
Windows7Professional
160707 
181221181221 
0003179: Уточнение формулы перевода из WGS84 в прямоугольные координаты СК-42
В файле Proj4.pas используется старые параметры для преобразования координат, в соответствии с ГОСТ Р 51794-2008 параметры преобразований из WGS84 в СК-42 были уточнены. В настоящее время следующая строка параметров для proj4 даёт точность порядка нескольких миллиметров и преобразует координаты из WGS84 в прямоугольные СК-42 в один приём:

+proj=latlong +datum=WGS84 +to +proj=tmerc +ellps=krass +towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22 +lon_0=39 +x_0=7500000

(последние два параметра зависят от листа). Проверить можно через онлайн-калькулятор http://www.latlong.ru/sk.php или через программу PHOTOMOD GeoCalculator (http://www.racurs.ru), которые выполняют преобразования координат в соответствии с последним ГОСТ от 2008 года.
системы координат
related to 0000376resolved zed Отображение координат в СК42 
Issue History
18-02-2017 15:22RedRatNew Issue
18-02-2017 18:20zedNote Added: 0017818
18-02-2017 18:31RedRatNote Added: 0017819
18-02-2017 18:32zedNote Added: 0017820
18-02-2017 18:41RedRatNote Added: 0017821
18-02-2017 18:48zedNote Added: 0017822
18-02-2017 19:28RedRatNote Added: 0017823
18-02-2017 20:34zedStatusnew => resolved
18-02-2017 20:34zedFixed in Version => 181221
18-02-2017 20:34zedResolutionopen => fixed
18-02-2017 20:34zedAssigned To => zed
18-02-2017 20:35zedTarget Version => 181221
18-02-2017 20:35zedSummaryУточнение формулы перевода из WGS84 в прямоугольные координаты => Уточнение формулы перевода из WGS84 в прямоугольные координаты СК-42
18-02-2017 20:35zedDescription Updatedbug_revision_view_page.php?rev_id=7133#r7133
18-02-2017 20:36zedTag Attached: системы координат
18-02-2017 20:38zedRelationship addedrelated to 0000376

Notes
(0017818)
zed   
18-02-2017 18:20   
У нас сейчас используются вот такие параметры:

+towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12

вы предлагаете заменить их на

+towgs84=23.57,-140.95,-79.8,0,0.35,0.79,-0.22

Я правильно вас понял?
(0017819)
RedRat   
18-02-2017 18:31   
Да, это параметры из последнего действующего ГОСТ-а.

Там ещё в Proj4.pas при вызове proj.dll передаются лишние параметры, и само преобразование делается, как я понял, в два этапа. На всякий случай я привёл строку параметров для proj.dll, которые преобразуют WGS84 (из GPS) сразу в прямоугольные координаты СК-42 (как на генштабовских картах). С ней у меня расхождение получается в третьем-четвёртом знаке после запятой (миллиметры и доли миллиметров).
(0017820)
zed   
18-02-2017 18:32   
Эти параметры подходят для перехода в обе стороны или чтобы перейти из СК-42 в WGS84 нужны другие?
(0017821)
RedRat   
18-02-2017 18:41   
Да, в обе стороны работают, надо только -I в параметры proj добавить. У меня есть небольшой Unix-скрипт, который вытаскивает из аэрофотосъёмки координаты GPS и преобразовывает их в прямоугольные генштабовские, так что я неделю примерно потратил, чтобы получить идеальное совпадение точек. Два вышеуказанных калькулятора можно использовать для проверки преобразований, они оба поддерживают ГОСТ Р 51794-2008.
(0017822)
zed   
18-02-2017 18:48   
Исправил и залил ночнушку, тестируйте.
(0017823)
RedRat   
18-02-2017 19:28   
Спасибо! Прямоугольные координаты совпадают тютелька в тютельку (с учётом округления до метра в программе).