ggk, rgis, eguar, smgis

Aleksandr Dezhin · 2011-02-10 20:18:34

Чтобы как-то систематизировать всяких генеральных планов и пр. фигни сделал штуку, которая умеет на лету делать gdal warp по контрольным точкам и тянуть нужные тайлы с проксируемого сервиса. На пробу сделал прокси для ggk.mos.ru, rgis.spb.ru и eguar.ru. Вот ссылки для josm imagery

Санкт-Петербург http://rgis.spb.ru
* привязка исходя из проекции описанной в вики

tms:http://osm.dezhin.net/tyler/tms/rgis-address/{zoom}/{x}/{y} - улицы + дома
tms:http://osm.dezhin.net/tyler/tms/rgis-boundary/{zoom}/{x}/{y} - границы

Екатеринбург http://eguar.ru
* максимальный зум больше дефолтного, можно поставить 21 в соотв колонку в настройках imagery

tms:http://osm.dezhin.net/tyler/tms/eguar-address/{zoom}/{x}/{y} - улицы + дома
tms:http://osm.dezhin.net/tyler/tms/eguar-boundary/{zoom}/{x}/{y} - границы районов, города и ГО

Москва (старая привязка) http://ggk.mos.ru
* присутствуют всякие косяки, типа эпизодически выпадающих тайлов и не очень контрастной цветовой палитры для светлых слоев, но вроде оно не сильно мешает, хотя со временем это поправится.

tms:http://osm.dezhin.net/tyler/tms/ggk/{zoom}/{x}/{y} - ggk.mos.ru как есть

Смоленская область http://admin.smolensk.ru/~gis_site/

tms:http://osm.dezhin.net/tyler/tms/smgis/{zoom}/{x}/{y}

Смоленск http://www.smoladmin.ru/map/mapr.html?mode=js&map=SmMap

tms:http://osm.dezhin.net/tyler/tms/smcity/{zoom}/{x}/{y}

Как скорректировать привязку?

1) переходим по ссылке http://osm.dezhin.net/tyler/layer/ggk и на вкладке находим интересующий нас фрагмент
2) стрелками на клавиатуре двигаем слой до совпадения с данными осм или с треками (слой с треками можно включить на последних трех zoom level-ах)
3) жмем add gcp, то что будет в поле под кнопкой нужно скопировать сюда или в личку и подождать когда привязка обновится

Last edited by Aleksandr Dezhin (2011-04-30 23:32:33)

Sergey Astakhov · 2011-02-10 20:57:08

Aleksandr Dezhin wrote:

Спб
* привязка исходя из проекции описанной в вики
tms:http://osm.dezhin.net/tyler/tms/rgi … m}/{x}/{y} - улицы + дома
tms:http://osm.dezhin.net/tyler/tms/rgi … m}/{x}/{y} - границы

Круть!
Правда сдвиг метров в 6 относительно данных OSM получается... Будем дальше проекцию подкручивать или так оставим?

Aleksandr Dezhin · 2011-02-10 21:13:29

Sergey Astakhov wrote:

Правда сдвиг метров в 6 относительно данных OSM получается... Будем дальше проекцию подкручивать или так оставим?

Я склонен к тому чтобы подкрутить. Сейчас сделано следующим образом: по проекции описанной в вики сгенерирована сетка из 64 точек. Сейчас можно добавлять точки прямо в нее в виде:
x y lon lat, где x, y координаты в rgis.spd.ru.
Через некоторое время заделаю чего-нибудь, чтобы можно было как-то по-удобней эти точки ставить.

dimuzz · 2011-02-11 08:01:08

Aleksandr Dezhin, супер!
По Екатеринбургу у меня получилось такое смещение относительно треков:
Проекция Меркатора Eguar Восточная-Декабристов 9.08290562762204E-8 -2.890015437184701E-6
Это смещение в окрестностях точки http://www.openstreetmap.org/?lat=56.82 … 8&layers=M

Fortress · 2011-02-11 08:17:07

Объясните ламеру, куда вставлять ссылку? http://tms:http://osm.dezhin.net/tyler/tms/eguar-address/{zoom}/{x}/{y}
Если в графу "Служебная ссылка" а потом "Получить слои" то джосм ругается "Неверная служебная ссылка" Если в графу "Ссылка на WMS" то вместо тайлов "Error". И ещё: 2 раза повторяющееся http: это правильно? Пробовал и так: "tms:http://osm.dezhin.net/tyler/tms/eguar-address/{zoom}/{x}/{y}" тоже не работает. josm-latest

Sergey Astakhov · 2011-02-11 08:41:52

Fortress wrote:

Объясните ламеру, куда вставлять ссылку? http://tms:http://osm.dezhin.net/tyler/tms/eguar-address/{zoom}/{x}/{y}
Если в графу "Служебная ссылка" а потом "Получить слои" то джосм ругается "Неверная служебная ссылка" Если в графу "Ссылка на WMS" то вместо тайлов "Error".

Вставлять в "Ссылка на WMS"

Fortress wrote:

И ещё: 2 раза повторяющееся http: это правильно?

Нет

Fortress wrote:

Пробовал и так: "tms:http://osm.dezhin.net/tyler/tms/eguar-address/{zoom}/{x}/{y}" тоже не работает. josm-latest

Только что попробовал - внутри Екатеринбурга работает. Версия JOSM 3880

Fortress · 2011-02-11 10:16:40

Sergey Astakhov wrote:

Вставлять в "Ссылка на WMS"
Fortress wrote:
И ещё: 2 раза повторяющееся http: это правильно?
Нет
Fortress wrote:
Пробовал и так: "tms:http://osm.dezhin.net/tyler/tms/eguar-address/{zoom}/{x}/{y}" тоже не работает. josm-latest
Только что попробовал - внутри Екатеринбурга работает. Версия JOSM 3880

Спасибо! Ргис работает. П/С: ссылки нужно поправить...

Aleksandr Dezhin · 2011-02-12 06:00:06

dimuzz wrote:

По Екатеринбургу у меня получилось такое смещение относительно треков:
Проекция Меркатора Eguar Восточная-Декабристов 9.08290562762204E-8 -2.890015437184701E-6

В общем я не одолел перевод этих смещений в какой-либо человеческий вид - если это меркатор, то это должны быть метры, но чего-то маловато для метров…
Привязал Екатеринбург по большему кол-ву точек, опираясь на ключевые дороги, должно получше быть. Возможно потребуется сбросить руками кеш тайлов imagery, flush tile cache похоже ничего не делает.

Last edited by Aleksandr Dezhin (2011-02-12 06:01:19)

Zverik · 2011-02-12 08:54:59

Это градусы в wgs84, скорее всего.

Sergey Astakhov · 2011-02-12 12:17:55

Aleksandr Dezhin wrote:

dimuzz wrote:
По Екатеринбургу у меня получилось такое смещение относительно треков:
Проекция Меркатора Eguar Восточная-Декабристов 9.08290562762204E-8 -2.890015437184701E-6
В общем я не одолел перевод этих смещений в какой-либо человеческий вид - если это меркатор, то это должны быть метры, но чего-то маловато для метров…

Upliner говорил что "Для меркатора это радианы по долготе и логарифм тангенса по широте". Вроде (как упоминал usm78-gis), достаточно просто домножить на 6378137 для получения обычных значений меркатора.

usm78-gis · 2011-02-12 13:44:41

Sergey Astakhov wrote:

Aleksandr Dezhin wrote:
dimuzz wrote:
По Екатеринбургу у меня получилось такое смещение относительно треков:
Проекция Меркатора Eguar Восточная-Декабристов 9.08290562762204E-8 -2.890015437184701E-6
В общем я не одолел перевод этих смещений в какой-либо человеческий вид - если это меркатор, то это должны быть метры, но чего-то маловато для метров…
Upliner говорил что "Для меркатора это радианы по долготе и логарифм тангенса по широте". Вроде (как упоминал usm78-gis), достаточно просто домножить на 6378137 для получения обычных значений меркатора.

Сдвиги направлены в стороны координатной сетки на которой они определены (и в единицах соответствующей координатной системы),
что несколько усложняет их практическое использование.
В примере 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.

Last edited by usm78-gis (2011-02-12 16:21:58)

usm78-gis · 2011-02-12 13:59:01

Sergey Astakhov wrote:

Будем дальше проекцию подкручивать или так оставим?

Надо бы СК-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

Last edited by usm78-gis (2011-02-12 16:18:13)

esaulenka · 2011-02-12 17:42:35

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

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

Aleksandr Dezhin · 2011-02-14 03:08:51

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.

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

Last edited by Aleksandr Dezhin (2011-02-14 10:04:12)

usm78-gis · 2011-02-14 20:15:04

Aleksandr Dezhin wrote:

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

Пока ничего лучше нет. +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.

wowik · 2011-03-29 18:03:33

Aleksandr Dezhin wrote:

Присутствуют всякие косяки

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

GaM · 2011-08-24 12:37:58

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

Zverik · 2011-08-24 13:10:35

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

GaM · 2011-08-24 13:50:32

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

Last edited by GaM (2011-08-24 13:50:48)

Aleksandr Dezhin · 2011-08-24 23:16:09

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

GaM · 2011-08-25 09:05:29

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

usm78-gis · 2011-08-28 14:48:37

usm78-gis wrote:

Вот гислабовская таблица из гуглокэша

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

На 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 нет, так что проверить не могу.

Дмитрий Кудряшов · 2011-08-28 15:28:46

usm78-gis wrote:

На geodesist.ru есть программа преобразования координат TransLine_05++.exe, которая вроде как умеет следующее:
    * Преобразование выполняется между следующими системами:
...         2. Система координат М64 - местная система Петербурга.
...
   
У меня windows нет, так что проверить не могу.

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

usm78-gis · 2011-08-28 18:38:47

Дмитрий Кудряшов wrote:

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

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

Дмитрий Кудряшов · 2011-08-28 19:26:36

usm78-gis wrote:

Дмитрий Кудряшов wrote:
Я проверял 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

#1 2011-02-10 20:18:34

ggk, rgis, eguar, smgis

#2 2011-02-10 20:57:08

Re: ggk, rgis, eguar, smgis

#3 2011-02-10 21:13:29

Re: ggk, rgis, eguar, smgis

#4 2011-02-11 08:01:08

Re: ggk, rgis, eguar, smgis

#5 2011-02-11 08:17:07

Re: ggk, rgis, eguar, smgis

#6 2011-02-11 08:41:52

Re: ggk, rgis, eguar, smgis

#7 2011-02-11 10:16:40

Re: ggk, rgis, eguar, smgis

#8 2011-02-12 06:00:06

Re: ggk, rgis, eguar, smgis

#9 2011-02-12 08:54:59

Re: ggk, rgis, eguar, smgis

#10 2011-02-12 12:17:55

Re: ggk, rgis, eguar, smgis

#11 2011-02-12 13:44:41

Re: ggk, rgis, eguar, smgis

#12 2011-02-12 13:59:01

Re: ggk, rgis, eguar, smgis

#13 2011-02-12 17:42:35

Re: ggk, rgis, eguar, smgis

#14 2011-02-14 03:08:51

Re: ggk, rgis, eguar, smgis

#15 2011-02-14 20:15:04

Re: ggk, rgis, eguar, smgis

#16 2011-03-29 18:03:33

Re: ggk, rgis, eguar, smgis

#17 2011-08-24 12:37:58

Re: ggk, rgis, eguar, smgis

#18 2011-08-24 13:10:35

Re: ggk, rgis, eguar, smgis

#19 2011-08-24 13:50:32

Re: ggk, rgis, eguar, smgis

#20 2011-08-24 23:16:09

Re: ggk, rgis, eguar, smgis

#21 2011-08-25 09:05:29

Re: ggk, rgis, eguar, smgis

#22 2011-08-28 14:48:37

Re: ggk, rgis, eguar, smgis

#23 2011-08-28 15:28:46

Re: ggk, rgis, eguar, smgis

#24 2011-08-28 18:38:47

Re: ggk, rgis, eguar, smgis

#25 2011-08-28 19:26:36

Re: ggk, rgis, eguar, smgis

Board footer