Еще раз о картах Росреестра

[Adon73]

За "категории" спасибо.., а сослаться на неустаревший Перечень (№ 326 от 28.05.2007) увы публично не могу (он же ДСП), но интересующимся космоснимками можно обрадовать, что в дополнении к пункту в перечне о материалах ДДЗ (где указано, что уточнение порядка будет сделано дополнительно) в 2007 году появилось несколько писем: от МО в адрес Роскартографии от 18 апреля 2007 г. № 205/6687 п.2 "Космические снимки с разрешением на местности лучше 2 метров, полученные с иностранных аппаратов, являются несекретными", соответствующее письмо было направлено от Федерального космического агентства на ЗАО "Совзонд" (от 05.07.2007) о том, что космические снимки с разрешением на местности лучше 2 м, полученные с иностранных аппаратов явлются не секретными, аналогичное письмо было направлено от Росреестра в ОАО "Газпром космические системы" в 2010 г. Надо признать, что появление общедоступного сервиса Google со снимками высокого разрешения сдвинула в мертвой точки данный вопрос, который был актуальным уже много лет...

[vasketsov]

Но если решу делать, то собстно программулина в течение трех-четырех дней напишется и будет одинаково работать по всей карте и всем регионам

А что делать с подобными плодами жизнедеятельности тамошних криворуких баранов (это запрос Layer: Центр КК (ID: 7) с условием FullNumber like '59:37:010%') на предмет номера квартала?:
OBJECTID: 357838
KvartalId: 59371050106
ParentId: 5937
OkrugNumber: 59
RayonNumber: 37
KvartalNumber: 0010501
FullNumber: 59:37:010501
ShortName:
XMin: 6162777,37768555
XMax: 6369549,26608555
YMin: 8179061,44171143
YMax: 8342590,07331143
X_coord: 6257579,56848145
Y_coord: 8266758,64611816
ParcelsCount: 0
ErrorCode: 1

[PavelML]

А что делать с подобными плодами жизнедеятельности тамошних криворуких баранов
KvartalId: 59371050106

Ждать появления в этих местах планеты мыслящих форм жизни 8))))
Кадастровый квартал 59:37:1050106 официальный портал росреестра найти "не смог". Видимо геодезисты еще не вышли из тени.
Зато нашел кадастровый массив 59:37:0000000
Край Пермской области пересекается с краем Коми-пермяцкого автономного округа на глубину... 19 км (!)
Ну это же неплохо если у Камского водохранилища есть сразу два хозяина, верно? Хуже когда ни одного нету 8)))
Вообще-то... кажется наоборот... границы РАСХОДЯТСЯ на 19 км... ну тогда четыре острова получаются ничьи... надо будет съездить прихватизировать 8)))
Или уже проспали? Пермская область уже стала Пермским краем?

Тут такое дело:
KvartalId: 59371050106 - что соответствует кадастровому номеру 59:37:1050106
OkrugNumber: 59 - так, RayonNumber: 37 - так, KvartalNumber: 0010501 - не так
FullNumber: 59:37:010501 - совсем не так!
в кадастровом номере
первые две цифры (59) - номер кадастрового округа (то есть субъекта федерации)
вторые две цифры (37) - номер кадастрового района (обычно совпадает с административным районом области)
далее должны быть семь цифр номера кадастрового квартала, которые на карте росреестра могут быть указаны шестизначными, если левая цифра ноль, но в запросе обязательно нужно указывать все семь цифр,
и наконец, если запрашивается земельный участок - после номера кадастрового квартала следует номер участка, который в официальных документах пишется после последнего двоеточия БЕЗ незначащих нулей слева (в отличие от номера кадастрового квартала), но в запросе незначащие нули также нужны и это число должно быть кажется пятизначным.
Как пример:
В кадастровой выписке написанное 59:37:0150106:13 будет равно на публичной карте числам "59:37, 150106 и 13" и коду в запросе 5937015010600013

ЗЫЖ кстати... кадастровый квартал 59:37:1050106 не нашел, а вот 59:37:1050105 - нашел у деревни Щекино... мож что в номере неправильно?

[vasketsov]

мож что в номере неправильно?

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

Похожая хрень (в смысле мегафееричности косяков) связана с районом
59:41
ShortName=ЗАТО Пермь-76
OBJECTID=941
XMin=6248219,95031738
XMax=6278025,46991738
YMin=7897096,9074707
YMax=7921200,8156707
X_coord=6263122,71008301
Y_coord=7905886,13427734
PointX=6263122,7101
PointY=7905886,1343
KvartalsCount=5
ParcelsCount=3816
ErrorCode=0

Для него запрос кварталов возвращает ровно одну запись
59:41:0000000
ShortName=
OBJECTID=356886
XMin=6248219,95031738
XMax=6278025,46991738
YMin=7897096,9074707
YMax=7921200,8156707
X_coord=6263122,71008301
Y_coord=7905886,13427734
PointX=6263122,7101
PointY=7905886,1343
KvartalNumber=0000000
KvartalId=59410000000
ParcelsCount=0
ErrorCode=0

и понимай как хочешь, то ли там 5 кварталов, то ли 1, и что мешало вернуть в списке районов в параметрах KvartalsCount=1 (то, что это ЗАТО, бардак не оправдывает).

С ErrorCode тоже не всё просто оказалось, местами хорошо бьётся, то если ErrorCode=0, то в KML будет граница, во всех прочих случаях не будет границы в KML, а местами в KML попадают и объекты с ErrorCode=пусто, причём не все, а как попало. Что 18-й, что 81-й округа, везде такой бардак, а иногда и с ErrorCode=0 и ненулевым номером квартала не попадает в KML (запрос границ). Кировская область вроде нормально получилась, но судя по всему нужно какую-то проверку данных выдумать, чтобы спать спокойно. Нижегородская и Марий Эл тоже косяков не замечал (но там не всё смотрел).

В общем пока автомат если и может работать, то только под пристальным вниманием автоматчика ))

[PavelML]

то, что это ЗАТО, бардак не оправдывает

Дык садица в камеру никто не желает...
как мы тут недавно нацитировали - режимные объекты нельзя координировать НИ В КАКОЙ системе координат! Секретно!
То есть даже просто нарисовать мелом на стене - и то нельзя!

По ходу у нас тоже ЗАТО имеется - небезызвестный ядреный центр Саров. Так вот - когда я по контракту выполнял координирование границ муниципальных образований - мне пришлось брать их тупо из гугля, по забору охранной зоны. 8)
Почему? Потому что размер этой зоны в реале оказался раза в три больше чем во всех официальных документах, в официальных документах граница ЗАТО нарисована аккурат по периметру городка Саров, в реальности... кто не ленивый может глянуть в гугле 8)
Кстати по публичной кадастровой карте росреестра ЗАТО (граница которого есстно не нарисована и сам Саров - входят в состав Русско-Караевского сельского поселения Республики Мордовия, но один Саров (в сильно погрызенном виде и без ЗАТЫ) при этом умудряется еще и входить в состав Нижегородской области! И все это безобразие - сельское поселение размером 40х30 км, город и пригород - имеет один единственный номер кадастрового квартала 0532001. А вот номера кадастрового округа и кадастрового района - двойные - 13:19 и 52:60... 8)

[PavelML]

Короче... способ получить из росреестра координаты высокой точности есть, но ручонками это очень канительно, надо программулинку делать.
Что я получил руками:
1. берем кадастровый номер нужного земельного участка (имеющегося в публичной кадастровой карте есстно)
2. Скармливаем браузеру строку запроса:

Код: Выделить всё

http://maps.rosreestr.ru/ArcGIS/rest/services/Cadastre/CadastreSelected/MapServer/1/query?text=&geometry=&geometryType=esriGeometryEnvelope&inSR=&spatialRel=esriSpatialRelIntersects&where=ParcelId+like+%271306031300300245%25%27&returnGeometry=true&outSR=4326&outFields=OBJECTID+%2CParcelId&f=kmz

где "1306031300300245" - кадастровый номер земельного участка, в миру он такой: 13:06:0313003:245
3. Браузер вываливает файло с расширением kmz и каким то абракадабрским названием, меняем его на человеческое и сохраняем.
4. Расзиповвыввем kmz в kml и выковыриваем оттуда (можно блокнотом) координаты точек земельного участка в системе WGS-84 с округлением до шестого знака. Практическая погрешность этих координат по всей территории субъекта у меня получилась не более ±8 см.
5. Кому надо точнее, тот дальше возится с парой координат каждой точки, формируя из них запрос для браузера:

Код: Выделить всё

http://maps.rosreestr.ru/ArcGIS/rest/services/Cadastre/Cadastre/MapServer/export?dpi=96&transparent=true&format=png32&bbox=46.536612%2C54.375690%2C46.536614%2C54.375692&bboxSR=4326&imageSR=102100&size=2048%2C2048&f=image

где "46.536612%2C54.375690%2C46.536614%2C54.375692" координаты рамки растра, полученные из округленных координат точки 46.536613 (долгота) и 54.375691 (широта), путем уменьшения и увеличения младшей цифры на единичку, разделенные "запятыми" в виде кода %2C.
6. После скармливания запроса браузеру он вываливает в окно (иногда требует пинка) растр 2048х2048 пикселя, отображающий диапазон строго ±0,000001 градуса широты и столько же долготы вокруг округленного значения искомой координаты.
7. Остается как-то попроще отвекторизировать две пересекающиеся прямые линии на растре с разрешением что-то около 100 микрон на пиксель и получить высокоточные координаты каждой точки в системе WGS-84 с погрешностью в доли миллиметра. Далее конвертировать в местную систему при наличии конвертора и все - можно выписки не покупать.

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

[PavelML]

По описанной выше методике получения точных координат в WGS-84 из карты росреестра решил уточнить свои датумы по своему региону 150х300 км. На входе - кадастровые выписки земельных участков в прямоугольной плоской метрической проекции в местной системе координат (как эталон), и полученные вышеуказанным способом точные координаты этих же участков в геоцентрической системе координат WGS-84.
Для калибровки взяты семь зон - 4 по углам субъекта и три - юг, центр и север по среднему меридиану. В каждой зоне взято по 8 точек для усреднения погрешностей.
Результат такой - датум подобрать удалось с высокой точностью, но есть нюансы.
Оказалось что координаты X "север" в МСК получены с точностью, которая максимально возможна в моих условиях, то есть усредненные отклонения по зонам не превышают ±0.2 мм (!), а максимальное отклонение координаты X любой отдельной точки составляет не более ±0,6 мм... примерно с такой точностью я и брал мышкой координаты из растров 2000х2000 пикселей, которые в геодезических координатах имеют размер 0,000002х0,000002 угловых градуса широты-долготы.
А вот с координатой Y "восток" получить удалось только среднюю погрешность в зонах не более ±2.0 мм, при этом разброс по отдельным точкам получился в пределах ±12 мм.

При этом никакой закономерности в погрешностях горизонтальных координат выявить не удалось, они носят абсолютно случайный характер. Были подозрения на линейность работы своей фунцкии (выложено в теме ранее) - не подтвердились. Проверил вероятность расхождений в исходных координатах сервера росреестра - скажем, если по закону координаты в выписках всегда округлены до целых сантиметров, они могли вываливаться в кадастровую карту неокругленными. Но это не подтвердилось - на примере участка, который я сам ставил на учет в конце 90-х с миллиметровой дискретностью и в базу росреестра они попали в миллиметрах - оказалось что правильно бьются только округленные до целых сантиметров координаты.
Остается сделать вывод о нелинейности серверной функции, формирующей координаты долготы. Это может быть связано либо с использованием где-то внутри нее переменной одинарной точности вместо двойной (возможно для повышения скорости обработки), либо с неизбежно итерационным характером функции вычисления геоцентрических координат из плоских прямоугольных - количество итераций уменьшено (точность грублена) - тоже для скорости вычислений.

Датум, обеспечивший мне указанную точность на территории 13-го региона:
dX 23.860000 м
dY -141.270000 м
dZ -80.900000 м
wX 0.000000"
wY -0.292000"
wZ -0.849000"
m 0.000000000000
смещ X -5 614 705.365502 м
смещ Y 1 250 001.591651 м
меридиан 44.550000000000º
исходный эллипсоид (WGS-84, мировая система)
6 378 137.000000 1/298.2572235630
результирующий эллипсоид (Красовского, референсная система)
6 378 245.000000 1/298.3000000000

Да, чуть не забыл... сотрудники росреестра с медвежъей грациозностью увернулись от уголовной статьи за секретность тем - что выдаваемые их сервером координаты WGS-84 в реальности смещены относительно истинных координат примерно на 40 метров по широте и примерно на 3 метра по долготе. Это имеет значение, если использовать этот датум для работы с гуглом, САСом и GPS-приемниками. При этом нужно сказать - что данные смещения не должны меняться по всей территории России, по крайней мере резко. Если бы они менялись - я бы на территории протяженностью 300 км при точности измерений на уровне долей миллиметра обязытельно бы это обнаружил.

[Shurik]

выдаваемые их сервером координаты WGS-84 в реальности смещены относительно истинных координат примерно на 40 метров по широте и примерно на 3 метра по долготе

А я-то давно пытаюсь понять, почему у меня углы участка, померенные бытовым навигатором, отклоняются от публичного кадастра существенно больше, чем на погрешность прибора...

UPD: а вот ни фига, у меня, если смотреть на глаз, скорее 20 метров на восток и 3 на юг. Или надо широту с долготой было понимать наоборот и 40 умножать на cos(lat)?

[PavelML]

выдаваемые их сервером координаты WGS-84 в реальности смещены относительно истинных координат примерно на 40 метров по широте и примерно на 3 метра по долготе

А я-то давно пытаюсь понять, почему у меня углы участка, померенные бытовым навигатором, отклоняются от публичного кадастра существенно больше, чем на погрешность прибора...

UPD: а вот ни фига, у меня, если смотреть на глаз, скорее 20 метров на восток и 3 на юг. Или надо широту с долготой было понимать наоборот и 40 умножать на cos(lat)?

Если сравнивать с навигатором - то он может абсолютные значения сам довольно сильно искажать. Причем это сложно определить сразу - показания смещаются по мере ухода за горизонт одних спутников и выхода других, то есть пара-тройка часов на одном месте - будет приводить к уводам координат в пределах 3-15 метров. Это не связано с точностью навигатора, такова физика атмосферных прохождений радиоволн и прочих факторов. Скажем такие же уводы, но несколько меньше - характерны и для геодезических приемников, одночастотник легко дает показания на три метра в сторону от вчерашних на этой же точке, и это несмотря на то что в режиме дифференциальной съемки его точность составляет 5мм в плане и 10 мм по вертикали. Лучшие модели геодезических двух-пятичастотников после сложной математической обработки эфемерид с компенсацией атмосферных задержек позволяют определить абсолютное положение примерно с точностью 0,3 метра.

Если сравнивать со спутниковой фотосъемкой... ну вот самые последние снимки - лета прошлого года, по нашему региону - посажены довольно точно. Тем не менее даже два соседних снимка, отснятых с разницей в один месяц и выложенных почти одновременно - разъезжаются между собой примерно на 6 метров, в основном по долготе. Если сравнивать с дигитал глобом и гугловскими прошлых лет - то там есть и 20 метров и 30 метров легко.

С учетом сказанного - во первых указанные мной 40 метров по широте и 3 метра по долготе не являются истиной, я это определял по очень ненадежным источникам - нафигатору и спутниковым снимкам гугла - причем только в одном месте, во вторых - я об этом сказал - что они "приблизительные".

Наша проблема в том - что мы не можем даже с помощью перепрограммирования САС-планеты - скомпенсирвоать эти сдвиги спутниковых снимков чтобы перепривязать их правильно - потому что насколько я знаю гугль выдает именованные тайлы, которые уже привязаны к сетке координат.
Но если говорить о карте росреестра - то там ситуация как раз получше, там запрос идет по координатам и мы можем сдвинуть всю карту росреестра у себя на компьютере - вставив в zmp-параметрах смещения в строке запроса для координат широты и долготы.
Проблема только в том - что это можно сделать точно с привязкой только к одному, самому важному снимку, остальные будут разъезжаться неминуемо. Вот я у себя почти сразу внес эти смещения и привязал карту росреестра к гугловскому снимку г. Саранска от 10.04.2010г. Запрос (в файле GetUrlScript.txt) выглядит так:

Код: Выделить всё

ResultURL:=GetURLBase+'bbox={"xmin":'+RoundEx(topLeftM.x,8)+',"ymin":'+RoundEx(topLeftM.y-69,8)+',"xmax":'+RoundEx(bottomrightM.x,8)+',"ymax":'+RoundEx(bottomrightM.y-69,8)+',"spatialReference":{"wkid":102113}}&bboxSR=102113&imageSR=102113&size=256,256&f=image';

добавлены смещения по широте, по долготе не потребовалось. В результате земельные участки очень четко совпали с этим конкретно снимком. Если интересно можете проверить по г. Саранску, как они подойдут другим снимкам, не знаю.
Но оговариваюсь - после ввода этих смещений всю карту придется перезаливать заново. Или сделать другой zmp с другим кэшем.
Как апогей фанатизма, можно сотворить столько клонов карт росреестра со своими собственными папками кэша и описателями zmp - сколько имеется актуальных спутниковых снимков высокого разрешения в ареале обитания юзверя... 8 )

[Shurik]

Если сравнивать с навигатором - то он может абсолютные значения сам довольно сильно искажать.

Но не настолько сильно, погрешность нормальных бюджетных навигаторов гораздо ниже. Я снимок подстраиваю под навигатор, а не наоборот, поэтому "разные сдвиги снимков" - это не ко мне.

Кстати, что за SR=102113? Вроде портал выдаёт по умолчанию всегда 102100, и оно в целом соответствует 900913. Может, у каждой 102* свой сдвиг просто, в этом и вся разница?

скрытый текст: показать

Код: Выделить всё

stdClass Object
(
    [attributes] => stdClass Object
        (
            [OBJECTID] => 12921992
            [ParcelId] => 402401708100000601
            [ParentId] => 40240170810
            [OkrugNumber] => 40
            [RayonNumber] => 24
            [KvartalNumber] => 0170810
            [ParcelNumber] => 00006
            [TypeCode] => 01
            [FullAddress] => Калужская область р-н Юхновский д Мочалово ул Дачная 28
            [Area008_Value] => 800
            [Area008_Unit] => кв. м.
            [Area009_Value] => 800
            [Area009_Unit] => кв. м.
            [Area010_Value] => 799.95
            [Area010_Unit] => кв. м.
            [CadastralPrice_Value] => 93936
            [CadastralPrice_Unit] => Рубль
            [Utilization_ReferenceValues] => Земли под личными подсобными хозяйствами|
            [Utilization_Fact] =>
            [Utilization_ByDocument] => для ведения личного подсобного хозяйства
            [Category] => Земли поселений (земли населенных пунктов)
            [X_coord] => 3913218.1823
            [Y_coord] => 7302345.3459
            [XMin] => 3913189.3414
            [XMax] => 3913247.0231
            [YMin] => 7302305.2695
            [YMax] => 7302385.7996
            [ErrorCode] => 0
            [StateCode] => 01
            [FullNumber] => 40:24:0170810:00006
            [RaionID] => 4024
            [LoadDate] => 1.2833856E+12
        )

    [geometry] => stdClass Object
        (
            [x] => 3913218.1823
            [y] => 7302345.3459
        )

)