You are not logged in.
- Topics: Active | Unanswered
Announcement
#1 2011-02-10 20:18:34
- Aleksandr Dezhin
- Member
- From: Msk
- Registered: 2008-05-25
- Posts: 2,109
- Website
ggk, rgis, eguar, smgis
Чтобы как-то систематизировать всяких генеральных планов и пр. фигни сделал штуку, которая умеет на лету делать 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)
[ Проверка границ и НП | Дампы | Валидатор рек | Кое-что ещё ]
- Нет! Не надо пытаться. Делай. Или не делай. Не надо пытаться. (c) Йода
Offline
#2 2011-02-10 20:57:08
- Sergey Astakhov
- Member
- From: St.Petersburg, Russia
- Registered: 2009-11-13
- Posts: 5,779
Re: ggk, rgis, eguar, smgis
Спб
* привязка исходя из проекции описанной в вики
tms:http://osm.dezhin.net/tyler/tms/rgi … m}/{x}/{y} - улицы + дома
tms:http://osm.dezhin.net/tyler/tms/rgi … m}/{x}/{y} - границы
Круть!
Правда сдвиг метров в 6 относительно данных OSM получается... Будем дальше проекцию подкручивать или так оставим?
Offline
#3 2011-02-10 21:13:29
- Aleksandr Dezhin
- Member
- From: Msk
- Registered: 2008-05-25
- Posts: 2,109
- Website
Re: ggk, rgis, eguar, smgis
Правда сдвиг метров в 6 относительно данных OSM получается... Будем дальше проекцию подкручивать или так оставим?
Я склонен к тому чтобы подкрутить. Сейчас сделано следующим образом: по проекции описанной в вики сгенерирована сетка из 64 точек. Сейчас можно добавлять точки прямо в нее в виде:
x y lon lat, где x, y координаты в rgis.spd.ru.
Через некоторое время заделаю чего-нибудь, чтобы можно было как-то по-удобней эти точки ставить.
[ Проверка границ и НП | Дампы | Валидатор рек | Кое-что ещё ]
- Нет! Не надо пытаться. Делай. Или не делай. Не надо пытаться. (c) Йода
Offline
#4 2011-02-11 08:01:08
- dimuzz
- Member
- From: Екатеринбург
- Registered: 2009-09-10
- Posts: 1,843
Re: ggk, rgis, eguar, smgis
Aleksandr Dezhin, супер!
По Екатеринбургу у меня получилось такое смещение относительно треков:
Проекция Меркатора Eguar Восточная-Декабристов 9.08290562762204E-8 -2.890015437184701E-6
Это смещение в окрестностях точки http://www.openstreetmap.org/?lat=56.82 … 8&layers=M
Offline
#5 2011-02-11 08:17:07
- Fortress
- Member
- From: Russia, Saint-Petersburg
- Registered: 2009-07-09
- Posts: 1,109
Re: ggk, rgis, eguar, smgis
Объясните ламеру, куда вставлять ссылку? 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
Offline
#6 2011-02-11 08:41:52
- Sergey Astakhov
- Member
- From: St.Petersburg, Russia
- Registered: 2009-11-13
- Posts: 5,779
Re: ggk, rgis, eguar, smgis
Объясните ламеру, куда вставлять ссылку? http://tms:http://osm.dezhin.net/tyler/tms/eguar-address/{zoom}/{x}/{y}
Если в графу "Служебная ссылка" а потом "Получить слои" то джосм ругается "Неверная служебная ссылка" Если в графу "Ссылка на WMS" то вместо тайлов "Error".
Вставлять в "Ссылка на WMS"
И ещё: 2 раза повторяющееся http: это правильно?
Нет
Пробовал и так: "tms:http://osm.dezhin.net/tyler/tms/eguar-address/{zoom}/{x}/{y}" тоже не работает. josm-latest
Только что попробовал - внутри Екатеринбурга работает. Версия JOSM 3880
Offline
#7 2011-02-11 10:16:40
- Fortress
- Member
- From: Russia, Saint-Petersburg
- Registered: 2009-07-09
- Posts: 1,109
Re: ggk, rgis, eguar, smgis
Вставлять в "Ссылка на WMS"
Fortress wrote:И ещё: 2 раза повторяющееся http: это правильно?
Нет
Fortress wrote:Пробовал и так: "tms:http://osm.dezhin.net/tyler/tms/eguar-address/{zoom}/{x}/{y}" тоже не работает. josm-latest
Только что попробовал - внутри Екатеринбурга работает. Версия JOSM 3880
Спасибо! Ргис работает. П/С: ссылки нужно поправить...
Offline
#8 2011-02-12 06:00:06
- Aleksandr Dezhin
- Member
- From: Msk
- Registered: 2008-05-25
- Posts: 2,109
- Website
Re: ggk, rgis, eguar, smgis
По Екатеринбургу у меня получилось такое смещение относительно треков:
Проекция Меркатора Eguar Восточная-Декабристов 9.08290562762204E-8 -2.890015437184701E-6
В общем я не одолел перевод этих смещений в какой-либо человеческий вид - если это меркатор, то это должны быть метры, но чего-то маловато для метров…
Привязал Екатеринбург по большему кол-ву точек, опираясь на ключевые дороги, должно получше быть. Возможно потребуется сбросить руками кеш тайлов imagery, flush tile cache похоже ничего не делает.
Last edited by Aleksandr Dezhin (2011-02-12 06:01:19)
[ Проверка границ и НП | Дампы | Валидатор рек | Кое-что ещё ]
- Нет! Не надо пытаться. Делай. Или не делай. Не надо пытаться. (c) Йода
Offline
#10 2011-02-12 12:17:55
- Sergey Astakhov
- Member
- From: St.Petersburg, Russia
- Registered: 2009-11-13
- Posts: 5,779
Re: ggk, rgis, eguar, smgis
dimuzz wrote:По Екатеринбургу у меня получилось такое смещение относительно треков:
Проекция Меркатора Eguar Восточная-Декабристов 9.08290562762204E-8 -2.890015437184701E-6В общем я не одолел перевод этих смещений в какой-либо человеческий вид - если это меркатор, то это должны быть метры, но чего-то маловато для метров…
Upliner говорил что "Для меркатора это радианы по долготе и логарифм тангенса по широте". Вроде (как упоминал usm78-gis), достаточно просто домножить на 6378137 для получения обычных значений меркатора.
Offline
#11 2011-02-12 13:44:41
- usm78-gis
- Member
- Registered: 2008-04-21
- Posts: 2,662
Re: ggk, rgis, eguar, smgis
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)
Offline
#12 2011-02-12 13:59:01
- usm78-gis
- Member
- Registered: 2008-04-21
- Posts: 2,662
Re: ggk, rgis, eguar, smgis
Будем дальше проекцию подкручивать или так оставим?
Надо бы СК-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)
Offline
#13 2011-02-12 17:42:35
- esaulenka
- Member
- From: Москва
- Registered: 2010-08-21
- Posts: 1,194
- Website
Re: ggk, rgis, eguar, smgis
С правкой смещений как-то грустно. В том смысле, что оно постоянно плавает... Проверял по бингу, его в свою очередь, выравнивая по трекам.
Да, я о Москве говорю.
Но в любом случае, штука очень полезная, т.к. пользоваться оригинальным ggk.mos.ru терпения не хватает.
Спасибо!
Offline
#14 2011-02-14 03:08:51
- Aleksandr Dezhin
- Member
- From: Msk
- Registered: 2008-05-25
- Posts: 2,109
- Website
Re: ggk, rgis, eguar, smgis
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)
[ Проверка границ и НП | Дампы | Валидатор рек | Кое-что ещё ]
- Нет! Не надо пытаться. Делай. Или не делай. Не надо пытаться. (c) Йода
Offline
#15 2011-02-14 20:15:04
- usm78-gis
- Member
- Registered: 2008-04-21
- Posts: 2,662
Re: ggk, rgis, eguar, smgis
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.
Offline
#16 2011-03-29 18:03:33
- wowik
- Member
- From: Zelenograd
- Registered: 2009-09-29
- Posts: 9,310
Re: ggk, rgis, eguar, smgis
Присутствуют всякие косяки
Еще сильно не хватает былой прозрачности слоя.
Offline
#17 2011-08-24 12:37:58
- GaM
- Member
- From: Saint-Petersburg/Russia
- Registered: 2011-01-10
- Posts: 3,187
- Website
Re: ggk, rgis, eguar, smgis
А куда все прокси удалились? очень РГис надо бы.
Чем меньше ты флудишь, тем больше нарисуешь!
Offline
#19 2011-08-24 13:50:32
- GaM
- Member
- From: Saint-Petersburg/Russia
- Registered: 2011-01-10
- Posts: 3,187
- Website
Re: ggk, rgis, eguar, smgis
ааа, ох. Ждём-с! Спасибо! Я бы зеркальнул РГис бы например на всякий случай.
Last edited by GaM (2011-08-24 13:50:48)
Чем меньше ты флудишь, тем больше нарисуешь!
Offline
#20 2011-08-24 23:16:09
- Aleksandr Dezhin
- Member
- From: Msk
- Registered: 2008-05-25
- Posts: 2,109
- Website
Re: ggk, rgis, eguar, smgis
Вроде бы все должно работать.
[ Проверка границ и НП | Дампы | Валидатор рек | Кое-что ещё ]
- Нет! Не надо пытаться. Делай. Или не делай. Не надо пытаться. (c) Йода
Offline
#21 2011-08-25 09:05:29
- GaM
- Member
- From: Saint-Petersburg/Russia
- Registered: 2011-01-10
- Posts: 3,187
- Website
Re: ggk, rgis, eguar, smgis
Спасибо все супер!
Чем меньше ты флудишь, тем больше нарисуешь!
Offline
#22 2011-08-28 14:48:37
- usm78-gis
- Member
- Registered: 2008-04-21
- Posts: 2,662
Re: ggk, rgis, eguar, smgis
Вот гислабовская таблица из гуглокэша
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 нет, так что проверить не могу.
Offline
#23 2011-08-28 15:28:46
- Дмитрий Кудряшов
- Member
- From: Санкт-Петербург, Пушкин
- Registered: 2010-06-15
- Posts: 444
Re: ggk, rgis, eguar, smgis
На geodesist.ru есть программа преобразования координат TransLine_05++.exe, которая вроде как умеет следующее:
* Преобразование выполняется между следующими системами: ... 2. Система координат М64 - местная система Петербурга. ...
У меня windows нет, так что проверить не могу.
Я проверял M64->WGS84 - работает, причём правильно работает.
Offline
#24 2011-08-28 18:38:47
- usm78-gis
- Member
- Registered: 2008-04-21
- Posts: 2,662
Re: ggk, rgis, eguar, smgis
Я проверял M64->WGS84 - работает, причём правильно работает.
Эээ. Точка 60N 30E (WGS84) имеет там какие координаты в "M64" ?
Совпадают ли они с результатом формулы в wiki ?
Offline
#25 2011-08-28 19:26:36
- Дмитрий Кудряшов
- Member
- From: Санкт-Петербург, Пушкин
- Registered: 2010-06-15
- Posts: 444
Re: ggk, rgis, eguar, smgis
Дмитрий Кудряшов 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
Offline