ggk, rgis, eguar, smgis

Сдвиги направлены в стороны координатной сетки на которой они определены (и в единицах соответствующей координатной системы),
что несколько усложняет их практическое использование.
В примере dimuzz точка 60.637912E, 56.826083N, epsg:4326 сдвигается на


$ echo 9.08290562762204E-8 -2.890015437184701E-6 | LANG=C awk '{R=6378137; print $1*R,$2*R}'
0.57932 -18.4329

+0.57932 -18.4329 меркаторовских метра от


$ echo 60.637912 56.826083 0.0 | cs2cs -f "%.3f" +init=epsg:4326 +to +init=epsg:3857
6750181.487     7724654.312 0.000

ее координат в меркаторовских метрах: 6750181.487E, 7724654.312N, epsg:3857
Если учесть, что 60.637912E, 56.826083N, epsg:4326 находится в зоне UTM41 (epsg:32641)
с координатами


$ echo 60.637912 56.826083 0.0 | cs2cs -f "%.3f" +init=epsg:4326 +to +init=epsg:32641
355855.845      6300514.112 0.000

355855.845E, 6300514.112N, epsg:4326, то после добавления “меркаторовской поправки”
она станет иметь координаты


$ LANG=C awk 'END{printf("%.4f %.4f %.4f\n",6750181.487+0.57932,7724654.312-18.4329,0.000)}' < /dev/null | cs2cs -f "%.3f" +init=epsg:3857 +to +init=epsg:32641
355855.815      6300504.019 0.000

355855.815E, 6300504.019N, epsg:32641, т.е. сдвинется ~ на +3см -10м в UTM41.

Надо бы СК-63 добавить, чтобы и область не обижать. Где-то на gislab форуме я видел таблицу proj4 параметров для C3 и C4 (?):


# Pulkovo 1942 / CS63 zone C2
<3352> +proj=tmerc +lat_0=0.1 +lon_0=27.95 +k=1 +x_0=2250000 +y_0=0 +ellps=krass +units=m +no_defs  <>

(при этом не забыть правильный +towgs84=)

Edit: Вот гислабовская таблица из гуглокэша, хотя вроде как
+lat_0=0.1 это не совсем 100% +y_0=11057.626 (?)


C2=TM, 27.95, 0, 1, 250000, -11057.626, 0
C3=TM, 30.95, 0, 1, 250000, -11057.626, 0
C4=TM, 33.95, 0, 1, 250000, -11057.626, 0
P1=TM, 32.48333333333, 0, 1, 250000, -12900.566, 0
P2=TM, 35.48333333333, 0, 1, 250000, -12900.566, 0
P3=TM, 38.48333333333, 0, 1, 250000, -12900.566, 0
P4=TM, 41.48333333333, 0, 1, 250000, -12900.566, 0

С правкой смещений как-то грустно. В том смысле, что оно постоянно плавает… Проверял по бингу, его в свою очередь, выравнивая по трекам.
Да, я о Москве говорю.

Но в любом случае, штука очень полезная, т.к. пользоваться оригинальным ggk.mos.ru терпения не хватает.
Спасибо!

usm78-gis, я запутан. Какую проекцию мне использовать? Сейчас это:

+proj=tmerc +lat_0=0 +lon_0=30 +k=1 +x_0=95948.85 +y_0=-6552807.1 +ellps=krass +units=m +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1 +no_defs

Могу попробовать поставить любую другую, если вы укажете ее параметры в виде строки proj4.

В шапке топика появилось описание, как скорректировать привязку.

Пока ничего лучше нет. +towgs84= такой же как и в mapinfo:


OGRFeature(Datum8parameters):2
  Datum# (String) = 1001
  Datum_Name (String) = Pulkovo 1942
  Ellipsoid_Name (String) = ELLIPSOID_KRASSOVSKY
  ShiftX (Real) = 24
  ShiftY (Real) = -123
  ShiftZ (Real) = -94
  RotationX (Real) = -0.02
  RotationY (Real) = 0.25
  RotationZ (Real) = 0.13
  scale_ppm (Real) = 1.1
  Prime_Meridian (Real) = 0

а со сдвигами в открытых источниках можно только гадать.
Вот в этом документе http://www.vsevreg.ru/files/munzakaz/kot735.doc пишут о местных системах МСК-64 и МСК-47, гугл опять же дает ссылку нa гислаб (cp. c CK63 выше):


"СК 47 зона 2", 8, 1001, 7, 30.95, 0, 1, 2250000, -6211057.626

но это никак не совпадает с моим кадастром.
Так что ничего другого не остается как улучшать +x_0 и +y_0 с помощью линейной
регрессии по GCP с фиксированными +lon_0 и +towgs84.

Еще сильно не хватает былой прозрачности слоя.

А куда все прокси удалились? :frowning: очень РГис надо бы.

умер винт, кэш копируется, скоро всё починят.

ааа, ох. Ждём-с! Спасибо! Я бы зеркальнул РГис бы например на всякий случай.

Вроде бы все должно работать.

Спасибо все супер!

На geodesist.ru есть программа преобразования координат TransLine_05++.exe, которая вроде как умеет следующее:


    * Преобразование выполняется между следующими системами:
         1. Система координат 1963г - зоны "С2", "С3", "С4" и "Р1", "Р2", "Р3";
         2. Система координат М64 - местная система Петербурга.
         3. Система координат 1942г - 5-я зона - UTM27 и 6-я зона - UTM33, названные так по градусу осевого меридиана
         4. Система геодезическких координат (B,L,H) на эллипсоиде Красовского "BLH(Kras)"
         5. Система геодезическких координат (B,L,H) на эллипсоиде WGS-84
         6. Региональная система координат 47-го региона - 1,2 и 3-я зоны - Ленинградская область
         7. Региональная система координат 53-го региона - 1,2 и 3-я зоны - Новгородская область

У меня windows нет, так что проверить не могу.

Я проверял M64->WGS84 - работает, причём правильно работает.

Эээ. Точка 60N 30E (WGS84) имеет там какие координаты в “M64” ?
Совпадают ли они с результатом формулы в wiki ?

Результат (с точностью до сантиметра ;)):

Ось X - вверх, ось Y - вправо.
А что получается по формулам в wiki?

Вот ещё пару точек (я их использовал при привязки планшетов 1:500, когда рисовал строящийся Суздальский пр.):
WGS-84 60,0680049 30,3040889; МСК-64 109000,113000; левый нижний угол планшета 2829-10-13
WGS-84 60,063547430,2905812; МСК-64 108500;112500; левый нижний угол планшета 2829-13-06
WGS-84 60,0591070 30,2681007; МСК-64 108000;111000; левый нижний угол планшета 2828-16-13
WGS-84 60,0591249 30,2591242; МСК-64 108000;110500; левый нижний угол планшета 2828-15-15

Во второй точке ошибка в 250 метров, а так они имеют систематический
сдвиг порядка 17 сантиметров по горизонтали :slight_smile:


МСК-64 113000 109000 -0.035 0.006
МСК-64 112500 108500 -250.03 -0.001
МСК-64 111000 108000 -0.023 -0.009
МСК-64 110500 108000 -0.024 -0.008

Если ваши точки определены с сантиметровой точностью, то надо использовать +x_0=95949 в wiki.

Да должно быть так:
WGS-84 60,0635474 30,2905812; МСК-64 108500;112250; левый нижний угол планшета 2829-13-06

А откуда видно 17 сантиметров?

Для 4-ой точки у меня получался коэффициент сжатия 0.998318863, что даёт тоже 17 см, но я его использовал не только по широте, но и по долготе. В добавок, планшет с этой же точкой надо повернуть на 0.2237226 градусов по часовой стрелке. Вот как выглядел файл привязки для PicLayer:


#\u0414\u0430\u043D\u043D\u044B\u0435 \u043A\u0430\u043B\u0438\u0431\u0440\u043E\u0432\u043A\u0438 \u043C\u043E\u0434\u0443\u043B\u044F JOSM PicLayer
#Fri Jun 24 01:52:39 MSD 2011
SHEARX=0.0
INITIAL_POS_y=8413031.911570232
SCALEY=17.64
SCALEX=17.64
POSITION_Y=8413031.709057873
INITIAL_POS_X=3368215.366136855
POSITION_X=3368215.703657454
INITIAL_SCALE=0.998318863
SHEARY=0.0
ANGLE=0.2237226

Если чисто формально преобразовать ваши lat/long WGS84 координаты
по формуле в wiki (при этом высота над WGS84 эллипсоидом берется =0),
то будет сдвиг порядка 17 сантиметров (и высота в Пулково-42 -15 метров, что уже само по себе некорректно).

Если действительно серьезно пытаться рисовать в ОСМ с дециметровой
точностью, то первое, от чего придется отказаться, это от сферического
меркатора, и добавить в josm все местные проекции (цифры EPSG/гислаба для CК-63-С и МСК-47 полностью совпадают с TransLine_05++) Тогда не придется делать ни сжатие, ни поворот. Но боюсь, что если все действительно делать правильно, то полезут на свет ошибки в системе кадастре (с его описаниями типа “пересечение канавы с канавой” и мифическими ОМЗ ) :slight_smile:

С СК-63 и МСК-47 вроде как все понятно, но остается одна интересная
загвоздка: реальные “кадастровые” цифры иногда не относятся ни
к одной из них (!). Вот конкретный пример:
http://www.gps-forum.ru/cgi-bin/forum/showpost.pl?Board=gpsgeneral&Number=108545&page=34&view=expanded&sb=2#Post108545
Таких координат в приозерском районе быть не должно.
Неудивительно, что данных ГКН для ЛО на сервере Росреестра до
сих пор нет :wink:

А можно попросить обновить кеш РГиса? А то ряд зданий там уже есть, а через wms не видать :frowning: