МСК и точность привязки ПКК.

Случайно набрел на таблицу параметров местных систем координат: http://www.mapbasic.ru/msksolutions
хотя там есть не все регионы, может будет кому-то полезна.
По моим наблюдениям, ошибки в 10 метров, несовпадения границ участков измеренных
разными фирмами такого же порядка дело абсолютно обычное,
что легко демонстрируется наложением в openlayers прозрачного WMS слоя ПКК росреестра
на гугловский спутниковый basemap (хотя и с ним надо быть осторожным).

А как это дело превратить в смещение для джосма?
У меня есть свое смещение кадастра для Татарстана, интересно, на сколько я ошибся.

Насколько я в курсе мы сейчас не используем ПКК.
А так две подложки двигаются след. образом:

  1. Сначала по трекам двигаем к примеру Бинговскую
  2. Вторую двигаем, так что бы совпало с бинговской.

Согласно постановлению политбюро ЦК OSMF на координаты отдельных точек (nodes) никакие
лицензионные ограничения не распространяются (PD), т.е. если подвинуть административные
границы примерно по ПКК, то такое действие ODBL-совместимо.
Это краеугольный камень т.н. “перелицензирования”.
ПКК же является самым точным источником открытых официальных данных (по крайней мере в теории).

Если вы используете проекцию Меркатора, то одним смещением тут не обойтись.
Для того чтобы добавить все МСК в josm, надо как-то договориться о “приватных” EPSG номерах
для них, вроде 9000NN, где NN номер МСК, при этом все, кто будет предоставлять WxS сервисы,
должны будут следовать также этому соглашению. В любом случае для PicLayer это будет
очень удобная и полезная добавка.

Если хотите что-то от PicLayer - сформулируйте задание поконкретнее, возможно даже в теме плагина. Авось кто-то сделает, у кого время будет :wink:

Ок, но это подложка, а что вы можете сказать про данные получаемые приближением (до участков) нажатии на кнопку информация, и получается адрес ввиде **обл, *** город, ** улица, *** дом. И допустимо ли в комментарии к правке писать: Обрисовка Бинг, адреса по кадастру?

Мое личное мнение - такую аттрибутивную информацию копировать нельзя.

Я уже вроде писал в теме плагина, что (имхо) главная проблема с PicLayer это превращение квадратных пикселов из
world файла в какую-то жуть в .cal файле из-за использования варварских приближенных формул для вычисления расстояний.
Может быть это и не сказывается на конечном результате (если приближенные формулы точно обратимы),
но я бы точно не стал пользоваться PicLayer в современном виде для обработки точных измерений.
А использование МСК с PicLayer уже возможно и сейчас: выбрав ‘Custom projection’ и
добавив proj4 строку в определение параметров. Правда +ellps=krass не распознается валидатором,
и приходится пользоваться вместо него эквивалентным “полным” определением +a=6378245 +rf=298.2997381
например, для Санкт-Петербурга


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

Для остальных МСК нужно только подставить соответствующие числа lon_0, x_0 и y_0 из mapinfo параметров
в таблице, т.е. для Татарстана из


"МСК-16 зона 1", 8, 1001, 7, 49.033333333333, 0, 1, 1300000, -5709414.70

получится


+proj=tmerc +lat_0=0 +lon_0=49.033333333333 +k=1 +x_0=1300000 +y_0=-5709414.70 +a=6378245 +rf=298.2997381 +units=m +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1 +no_defs

При сравнении ПКК и использовании координат СК-63 зоны C2 (EPSG:3352 http://spatialreference.org/ref/epsg/3352/ , она же МСК-47 зона 1 со сдвигом 1000000;6200000)


+proj=tmerc +lat_0=0 +lon_0=27.95 +k=1 +x_0=250000 +y_0=-11057.626 +ellps=krass +units=m +towgs84=24,-123,-94,0.02,-0.25,-0.13,1.1

у меня получился сдвиг (в реальных метрах) -8.741 ; -1.568
Так что без независимых измерений координат правильно установленных межевых знаков (= в соответствии с законодательством http://www.rg.ru/2010/03/02/mezhevye-znaki-dok.html , а не той липы, которой обычно занимаются кадастровые конторки)
продвинуться в этом направлении невозможно.

Если использовать для towgs84 простое 3x параметрическое преобразование (ПЗ-90 → СК-95)


+proj=tmerc +lat_0=0 +lon_0=27.95 +k=1 +x_0=250000 +y_0=-11057.63 +a=6378245 +rf=298.3 +units=m +towgs84=25.90,-130.94,-81.76,0.000,-0.0,0.0,0.00 +no_defs

то горизонтальный сдвиг полностью изчезает, но остается вертикальный сдвиг на -2.38 метра.

А кто-нибудь может подсказать как из координат, полученных с главной, получить координаты в МСК (конкртено, Тутаевский район Ярославской области, соответственно, в МСК-76).

Берем для примера деревню Летешовка Тутаевского района Ярославской области:
http://www.openstreetmap.org/?lat=57.796863&lon=39.736436&zoom=18&layers=M


МСК-76 Ярославская область
"МСК-76 зона 1", 8, 1001, 7, 38.55, 0, 1, 1250000, -6014743.504
"МСК-76 зона 2", 8, 1001, 7, 41.55, 0, 1, 2250000, -6014743.504

3 градусная зона 1 с центром 38.55° (37.05° ↔ 40.05°), точка с lon=39.736436° находится в ней.

Устанавливаете в JOSM Edit->Preferences->Map Settings->Projection method = “Сustom projection” с параметрами


+proj=tmerc +lat_0=0 +lon_0=38.55 +k=1 +x_0=1250000 +y_0=-6014743.504 +a=6378245 +rf=298.3 +units=m +towgs84=25.90,-130.94,-81.76,0.000,-0.0,0.0,0.00 +no_defs

и Edit->Preferences->Map Settings->Display Coordinates as = Projected Coordinates
скачиваете векторную карту OSM с сервера, и видите в левом нижнем углу
(строго говоря примерные, т.к. точные цифры для +towgs84= для этого региона неизвестны)
координаты в “МСК-76 зона 1”.

Для точки latWGS84=57.796863 lonWGS84=39.736436 eleWGS84=0 у меня получается
EastingМСК-76з1=1320670.676474938 NorthingМСК-76з=394634.101737767 eleМСК-76з1=-6.757255170

Если опустить номер зоны (1) и высоту, то получится X=394634.1 Y=320670.7 (Y=Easting, X=Northing)


$ cs2cs -f "%.9f" +init=epsg:4326 +to +proj=tmerc +lat_0=0 +lon_0=38.55 +k=1 +x_0=1250000 +y_0=-6014743.504 +a=6378245 +rf=298.3 +units=m +towgs84=25.90,-130.94,-81.76,0.000,-0.0,0.0,0.00 +no_defs
39.736436 57.796863 0
1320670.676474938       394634.101737767 -6.757255170

Хм… а без JOSMа? есть какие-то формулы? Я бы даже общедоступный пхп-скриптик под это дело слепил бы, куда координаты с главной вводишь, а он выдаёт в в МСК-ХХ, скажем, только бы формулы где взять.

Вся математика основана на proj4 http://trac.osgeo.org/proj , и вроде как даже для php есть
http://proj4php.svn.sourceforge.net/viewvc/proj4php/trunk/src/proj4php
На форуме гислаба пишут, что все параметры МСК-43 официально опубликованы
http://www.to43.rosreestr.ru/kadastr/program/mapinfo