SASGIS - SAS.Планета |
View Issue Details |
|
ID | Project | Category | View Status | Date Submitted | Last Update |
0003179 | SAS.Планета | Рефакторинг | public | 18-02-2017 15:22 | 18-02-2017 20:38 |
|
Reporter | RedRat | |
Assigned To | zed | |
Priority | normal | Severity | minor | Reproducibility | always |
Status | resolved | Resolution | fixed | |
Platform | Windows | OS | 7 | OS Version | Professional |
Product Version | 160707 | |
Target Version | 181221 | Fixed in Version | 181221 | |
|
Summary | 0003179: Уточнение формулы перевода из WGS84 в прямоугольные координаты СК-42 |
Description | В файле 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 года. |
Steps To Reproduce | |
Additional Information | |
Tags | системы координат |
Relationships | related to | 0000376 | resolved | zed | Отображение координат в СК42 |
|
Attached Files | |
|
Issue History |
Date Modified | Username | Field | Change |
18-02-2017 15:22 | RedRat | New Issue | |
18-02-2017 18:20 | zed | Note Added: 0017818 | |
18-02-2017 18:31 | RedRat | Note Added: 0017819 | |
18-02-2017 18:32 | zed | Note Added: 0017820 | |
18-02-2017 18:41 | RedRat | Note Added: 0017821 | |
18-02-2017 18:48 | zed | Note Added: 0017822 | |
18-02-2017 19:28 | RedRat | Note Added: 0017823 | |
18-02-2017 20:34 | zed | Status | new => resolved |
18-02-2017 20:34 | zed | Fixed in Version | => 181221 |
18-02-2017 20:34 | zed | Resolution | open => fixed |
18-02-2017 20:34 | zed | Assigned To | => zed |
18-02-2017 20:35 | zed | Target Version | => 181221 |
18-02-2017 20:35 | zed | Summary | Уточнение формулы перевода из WGS84 в прямоугольные координаты => Уточнение формулы перевода из WGS84 в прямоугольные координаты СК-42 |
18-02-2017 20:35 | zed | Description Updated | bug_revision_view_page.php?rev_id=7133#r7133 |
18-02-2017 20:36 | zed | Tag Attached: системы координат | |
18-02-2017 20:38 | zed | Relationship added | related 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
|
|
Спасибо! Прямоугольные координаты совпадают тютелька в тютельку (с учётом округления до метра в программе). |
|