Кадастр

Из закона Об информации, информационных технологиях и о защите информации:

Значит есть (могут быть) ограничения.

Могут быть, но информация кадастра, выложенная в общий доступ на сайт, является общедоступной. Не может использоваться там информация, которая не является общедоступной. Эта - является. Значит может использоваться.

кадастровые кварталы не совпадают с административным делением, по крайней мере я могу привести примеры, которые показывают
более сложные отношения (примеры общедоступные, можно найти поиском на google)

По этому поводу уже многие годы ломают копья.
Вот что например говорит GDAL/OGR (первая цифра - EPSG номер):

$ fgrep CS63 /usr/local/share/gdal/pcs.csv
2935,"Pulkovo 1942 / CS63 zone A1",9001,4284,18441,9807,1,0,8801,0.07,9110,8805,1,9201,8807,0,9001,8806,1300000,9001,8802,41.32,9110,,,,,,
2936,"Pulkovo 1942 / CS63 zone A2",9001,4284,18442,9807,1,0,8807,0,9001,8801,0.07,9110,8806,2300000,9001,8802,44.32,9110,8805,1,9201,,,,,,
2937,"Pulkovo 1942 / CS63 zone A3",9001,4284,18443,9807,1,0,8801,0.07,9110,8805,1,9201,8807,0,9001,8802,47.32,9110,8806,3300000,9001,,,,,,
2938,"Pulkovo 1942 / CS63 zone A4",9001,4284,18444,9807,1,0,8801,0.07,9110,8805,1,9201,8807,0,9001,8802,50.32,9110,8806,4300000,9001,,,,,,
2939,"Pulkovo 1942 / CS63 zone K2",9001,4284,18446,9807,1,0,8807,0,9001,8805,1,9201,8802,50.46,9110,8806,2300000,9001,8801,0.08,9110,,,,,,
2940,"Pulkovo 1942 / CS63 zone K3",9001,4284,18447,9807,1,0,8802,53.46,9110,8801,0.08,9110,8806,3300000,9001,8805,1,9201,8807,0,9001,,,,,,
2941,"Pulkovo 1942 / CS63 zone K4",9001,4284,18448,9807,1,0,8801,0.08,9110,8802,56.46,9110,8806,4300000,9001,8805,1,9201,8807,0,9001,,,,,,
3350,"Pulkovo 1942 / CS63 zone C0",9001,4284,18450,9807,1,0,8807,0,9001,8805,1,9201,8802,21.57,9110,8801,0.06,9110,8806,250000,9001,,,,,,
3351,"Pulkovo 1942 / CS63 zone C1",9001,4284,18451,9807,1,0,8807,0,9001,8805,1,9201,8806,1250000,9001,8802,24.57,9110,8801,0.06,9110,,,,,,
3352,"Pulkovo 1942 / CS63 zone C2",9001,4284,18452,9807,1,0,8801,0.06,9110,8806,2250000,9001,8807,0,9001,8802,27.57,9110,8805,1,9201,,,,,,

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

Пример городской округ Муром и Муромский муниципальный район.

Округ Муром включает в себя город Муром и девять сельских населенных пунктов. В Муромский муниципальный район городской округ Муром не входит. Как городской округ Муром так и Муромский муниципальный район имеют свои равноправные органы местного самоуправления.

Одно из отличий от того, что было раньше: закон обязал четко определить границы каждого муниципального образования. Лежит куча мусора - кто за нее отвечает? Нашли нефть - в бюджете какого муниципального образования будут дополнительные деньги.

Есть статья в википедии http://ru.wikipedia.org/wiki/%D0%9C%D1%83%D0%BD%D0%B8%D1%86%D0%B8%D0%BF%D0%B0%D0%BB%D1%8C%D0%BD%D0%B0%D1%8F_%D1%80%D0%B5%D1%84%D0%BE%D1%80%D0%BC%D0%B0_%D0%B2_%D0%A0%D0%BE%D1%81%D1%81%D0%B8%D0%B8_(2006)

Дополнительно к этой статье можно почитать какие вопросы находятся в ведении муниципальных образований: статьи 14, 15, 16 закона “Об общих принципах организации местного самоуправления в Российской Федерации” http://www.kodeks.ru/noframe/com-pus-FullLegRF?d&nd=901876063&nh=1

Административное деление – классификатор ОКАТО, обновляется регулярно, используется повсеместно.
Муниципальное деление – классификатор ОКТМО.

Может быть это и система координат такая, но координаты там похоже в метрах - по крайней мере расстояния между контрольными точками сходятся с точностью до 0.5%.
Пробежался по http://spatialreference.org/, но подходящей системы координат не нашел (

Мне кажется, что принадлежность какой-то территории к муниципальному образованию, мало что говорит, относительно того где это место находится, а говорит о том, как и кем эта территория управляется, что на мой взгляд, не представляет особого интереса для карты, хотя релейшенов в будущем наплодить можно.

Чтобы ответить на этот вопрос, надо знать границы административного района. После муниципальной реформы мне непонятна сущность такого понятия как административный район. Где он используется? Если только в ОКАТО, то где тогда используется ОКАТО? Границы муниципального образования описаны в его уставе, а вот где искать границы административного района, если непонятно как где это понятие сейчас используется?

Чтобы задать систему координат типа CS63 нужны 4 параметра LA0, LO0, X0, Y0 :

+proj=tmerc +lat_0=LA0 +lon_0=LO0 +k=1.000000 +x_0=X0 +y_0=Y0 +ellps=krass +towgs84=33.4,-146.6,-76.3,-0.359,-0.053,0.844,-0.84 +units=m +no_defs

В тех известных EPSG кодах значения X0 и Y0 стандартные, а LA0 и LO0 сдвинуты на “простую” дробь градуса.
На форуме gps-forum.ru было много постингов с координатными парами для Московской области, так что
не должно быть проблем с “теоретическим” вычислением LA0, LO0, X0, Y0 уже по двум известным парам.

Не похоже, а именно в метрах. Там так и написано “неземная система координат, метры”.
На сайте море подобных СК. Например EPSG:2466, EPSG:2554, ну и куча подобных могут быть довольно близкими (если сдвиг добавить).

Сдвиг определяется только изменением X0 и Y0. При таком подходе достаточно знать координаты 1 (!) точки чтобы сделать пересчет
любой другой точки. Не думаю, что создатели местных систем координат такие уж идиоты.

  1. Я сказал могут, а не в самом деле окажутся, например где-то могли взять одну из всеросийских СК и добавить туда эти самые X0, Y0 :wink:
  2. Кроме X0, Y0 может потребоваться Phi
  3. Оптимально найти базовую точку проекции (где плоскость касается сферы), и пересчитывать оттуда. Причем это не 0,0 так как по правилам создания МСК (которых не во всех районах придержитаются) требуется добавить дельту, чтобы все координаты в зоне действия СК была положительными.

    n. Скорее всего даже неправильные параметры перехода от этой СК дадут приемлимую для OSM точность (в пару метров).
    P.S. а еще в МСК могут как-то учитывать кривизну поверхности (не знаю как но упоминание о такой возможности встречал).

Там еще вот такое есть - несмотря на то что растр, векторизовать его можно довольно просто.
Границы областей там по-подробней тех что есть сейчас. Но насчет совпадания границ кадастрового и муниципального деления пока нет, но городские округа там отдельно, а муниципальные районы отдельно.
Ну и опять очередная система координат - хотя тут попроще должно быть.
Оно нам надо?

если б сырые данные получить, то было б неплохо.А векторизировать растр какой смысл…

может кто попросит у админа границы?

Поковырял MIF файлы для московской области, получилось что-то типа:

+proj=tmerc +lat_0=0 +lon_0=38.479571 +k=1 +x_0=3249848.090161 +y_0=-12151.181141 +ellps=bessel +units=m +no_defs (ср .ошибка 12 метров)
+proj=tmerc +lat_0=0 +lon_0=38.482197 +k=1 +x_0=3250017.613303 +y_0=-12900.660956 +ellps=krass +units=m +no_defs (ср. ошибка 16 метров)

Точек 7 штук - искать их очень мучительно оказалось.
Некоторые районы в этих файлах сдвинуты на тысячи км в разные стороны )

Спасибо, посмотрим.
То, что разные районы сдвинуты в разные стороны было очевидно (по координанам).
Кстати, в Тверской области написано:

CoordSys Earth Projection 8, 1001, "m", 1.5, 0, 1, 1250000, 0 Bounds (-6999281.53901, -10002137.4978) (9499281.53901, 10002137.4978)

Что вообще должно быть Земной СК с проекцией “Transverse Mercator”!!!

В порядке эксперимента провел по полученным данным провел границы Сергиево-Посадского и Талдомского районов и совпадающую с ней границу МО, что всяко лучше чем GNS который там был и были НП, которые по карте были во Владимирской области, а на самом деле они в МО. Пока рисовал границы - появились мысли как обозначать, например такой момент, что для way-я границы, которая собирается в relation с type=boundary почти никогда не надо указывать admin_level, так как непонятно к чему он относится.

Мини-конвертер mif2osm собрал, чтобы в JOSM-е смотреть, так что если кому интересно, могу поделится предварительными версиями osm файлов, районов, которых не разметало по вселенной. Готовые есть для Дмитровкого, Талдомского и Сергиево-Посадского, другие можно быстренько сгенерить. Мне кажется оттуда много чего можно будет натаскать в OSM - у меня прилегающая к даче территория очень не дурно легла.

Ну и появился скриптик, который симлексами подбирает параметры проекции для proj4 - может пригодится в твери.
Хотя если честно, я совсем недавно узнал про всякие проекции, и не совсем понял, что нам даст, тот факт, что там написано, что это земная СК )

Естественно все на горячо любимом нами perl-е.

Вообще-то можно указываеть максимальный (т.е. минимальное число) уровень границы. Например для куска границы НП, совпадающей с границей области указывать “admin_level=4”, и т.д., поскольку мало кто по гос.границе еще ведет границу соответсвующего района. Тогда и простым рендерерам (которые не смотрят на отношения) будет полегче они просто проведут нужную в данном месте границу.

А на сам конвертер (и скриптик) посмотреть нельзя?

По идее это должно дать отсутствие надобности в подобных скриптах. Т.е. здесь имеет быть проекция с известными параметрами перехода в тот же WGS-84, в отличии от МСК, где ключи перехода лежат с 1-ых отделах и, скорее всего, тем, кто хочет иметь дело с OSM, их лучше не видеть.
К сожалению мои попытки запустить “ogr2ogr -t_srs “EPSG:4326” out Query1.MIF” дали что-то “размазанное по вселенной”. Или шапка файла не соответсвует действительности, или данные 1001-ого датума беруться неизвестно откуда, или, что наверняка и происходит, я не умею готовить огров.

Кстати, погрешность порядка 10-м вполне соответсвует точности бытовых GPS (+точности установки точек…) и точности данных в OSM (что, к тому же, неплохо соотносится с параноидальными устремлениями нашего государства). Отсюда же следует непригодность импорта данных из OSM в кадастры и прочие серьезные БД, где требуется субметромая точность .

Да, полностью согласен. Сам потом включил голову и к такому-же выводу пришел. (right|left):region предлагаю заполнять в виде: “Московская область, Шатурский район” записывая туда четные уровни из admin_level. Кстати простановку и *:region и admin_level для вэя границы можно поручить боту. Я так понял admin_level нужен всем рендерам ).

http://nopaste.info/b54dd43960.html

./mif2osm --file ../cadastr/50/group/X05-0.MIF --proj "+proj=tmerc +lat_0=0 +lon_0=38.479596 +k=1 +x_0=3249849.890208 +y_0=-12150.598631 +ellps=bessel +units=m +no_defs" > ~/tmp/05.osm

Ну тогда подберем, если надо будет.

У меня закрадывается подозрение, что у yahoо местами ошибка привязки метров 20-25 - в чистом поле треки лежат кучно, но на дорогу по снимку не попадают. Попробую сделать контрольные точки не по снимку, а по osm карте. То что есть у меня сейчас лучше ложится на яховские снимки, нежели на osm карту.

У снимков что, Yahoo, что google, имеется примерно подобрая погрешность. Причем она плавает от точки к точке. Подозреваю, что это погрешность орторектификации. Когда снимок, полученный под углом к нормали, разворачивают в прямоугольный, то приходится учитывать рельеф (точки с разной высотой отклоняются в стороны по разному), и при не очень точной 3-d модели поверхности получаются разные погрешности (посмотри насколько крыши домов отъезжают от фундамента).