SASGIS

Веб-картография и навигация

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

Материалы от РосКосмоса, ГГЦ, РосРеестра, КосмоСнимков

Модератор: Tolik

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

Сообщение PavelML » 30 мар 2011, 12:37

Adon73 писал(а):Почему используете ГОСТ 51794-2001, если он отменен и действует ГОСТ 51794-2008, там другие параметры...

потому что функция была написана ДО появления этого госта. А "другие" параметры - это другой датум. Зачем он мне - если в любом случае требуется подбор индивидуальных для каждой МСК?
По округлению до шести знаков (дискретность 11,2 см, а для долготы 6,3 см), имеет ли это значение если точность пересчета МСК-WGS84-МСК составляет 6 - 8 м, т.е. грубее на несколько порядков.

Мои функции дают точность пересчета для WGS-84 --> СК-13 --> WGS-84 (мордовия) примерно 3 мм по всей территории размером 150х300 км. Остальные пока не проверял. Точнее сказать - при подгоне датума я остановился на этой погрешности, не стал "догонять" лучше.
Статья о гостайне здесь не причем - т.к. мы говорим о кадастровой информации, которая по закону открыта в независимости от точности данных, т.к. режимные объекты, не наносятся на публичную кадастровую карту.

Ошибаетесь. К гостайне отнесень точное определение координат любой точки в одной из МИРОВЫХ систем координат (например WGS-84). Кадастровая информация дана СУГУБО В МЕСТНЫХ системах координат. Ваши просьбы были по поводу преобразования - то есть получения точных координат в МИРОВОЙ системе из местной.
PavelML
Заслуженный тролль ресурса
 
Сообщения: 104
Зарегистрирован: 20 фев 2010, 17:29
Благодарил (а): 0 раз.
Поблагодарили: 6 раз.

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

Сообщение PavelML » 30 мар 2011, 12:52

Shurik писал(а):
tae1980 писал(а):В кадастровом паспорте моей дачи заявлена погрешность 50 см. Если числовая вариативность данных на сайте меньше - нужно ли париться?


Объясняю. Заявленная погрешность - нормативная, для геодезической съемки. Скажем для городских съемок она минимальна - 10 см. А для колхозных сельхозугодий - 2,5 м. Для леса если не ошибаюсь - 5 м.
Но! Дело в том что после того как координаты поставлены на учет - все эти погрешности уже теряют силу. Данные зафиксированы и их нельзя менять. Дискретность цифрования в государственном кадастре недвижимости - 1 см. Поэтому - если Вы определите координату смежного участка с мизерной погрешностью в 1 см - ваш участок может оказаться пересекающимся с ранее учтенным смежным участком. И это будет основанием для отказа в постановке его на учет. Имеется в виду именно кадастровая информация, а не реальная на местности. Обычно дело когда малограмотный или малопорядочный "кадастровый инженер" ставит на учет скажем, участок под боксом гаража - который по координатам оказывается примерно в 70 метрах от самого бокса гаража! Так вот - в картах росреестра выдаются учтенные данные, а не фактические. Поэтому указанные допуски погрешности к ним неприменимы. Они должны быть получены точными.
PavelML
Заслуженный тролль ресурса
 
Сообщения: 104
Зарегистрирован: 20 фев 2010, 17:29
Благодарил (а): 0 раз.
Поблагодарили: 6 раз.

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

Сообщение PavelML » 30 мар 2011, 13:10

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


Налогообложение тут ни при чем, роснедвижимость никак не связана с налоговыми органами (по крайней мере юридически). Были случаи, когда кадастровая площадь оказывалась вдвое больше фактической (и налог соотвественно больше) - но чаще бывает наоборот, когда владелец использует больше земли чем у него прописано в документах.
Земельные участки без координат возникают по трем причинам:
1. Ранее учтенные в "докомпьютерную эпоху". Докомпьютерная эпоха в роснедвижимости - это не сильно большая древность, примерно 1995-2002 годы (в зависимости от региона России). Тогда границы участков рисовались "на кальке", положенной на планшеты, рисованные тушью по аэрофотоснимкам. Само собой координаты тогда уже существовали, но погрешность такой технологии не лезла ни в какие нормативы. А снимать участки геодезическими приборами (каждую точку) было дороже стоимости самих земельных участков. Собсстно, базы данных о земельных участках даже не предусматривали полей под ввод координат. Базы это сильно сказано для экселевских таблиц.
2. Массово поставленные на кадастровый учет администрациями населенных пунктов в переходной период уже имели координаты и были занесены в базы данных. Но - отсутствует правовая сторона в узаконивании координат границ - то есть согласование с "заинтересованными лицами". По закону - земельный участок остается учтенным, поскольку есть права на участок. Но остается без координат, поскольку нет прав на координаты. Скажем, в моем городе несколько лет назад по этой причине лишились координат примерно 70% всех земельных участков.
3. Отмежеванные земельные участки могли быть поставлены на учет в локальных системах координат (система координат в пределах одного наспункта). А еще они могли быть поставлены с ошибками - которые позже были оспорены другими межевиками. Такие участки тоже лишают координат в материалах росреестра.
PavelML
Заслуженный тролль ресурса
 
Сообщения: 104
Зарегистрирован: 20 фев 2010, 17:29
Благодарил (а): 0 раз.
Поблагодарили: 6 раз.

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

Сообщение Adon73 » 30 мар 2011, 14:51

Небольшое замечания - в перечне материалов, подлежащим засекречиванию по линии Раскартографии (сейчас Росреестра) говорится о государственных системах координат, а WGS-84 никогда государственной в РФ не считалась... и это требование распространяется на пункты государственной геодезической сети (даже не на пункты ОМС), но не на точки границ земельных участков (больше половины из которых не закреплены на местности).
По поводу сроков обновления данных на ПКК, посмотрите на ней есть информациия на какую дату были загружены материалы, а заявленное регулярное обновление похоже пока не работает...
... Аналогичная история было 2 зарегистрированных земельных участка с 2003 года, кадастровая палата по каждобу выдает выписки и паспорта, но на публичной кадастровой карте одного из них нет...
интересно другое - когда отсутствующий участок решили разделить на 2, подготовили межевой план - то получили отказ в кадастровой палате, т.к. оказывается участок имел наложение на смежный участок на 8 см...
Adon73
Новичок
 
Сообщения: 49
Зарегистрирован: 10 июл 2010, 16:22
Благодарил (а): 0 раз.
Поблагодарили: 7 раз.

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

Сообщение PavelML » 30 мар 2011, 15:25

Adon73 писал(а):Небольшое замечания - в перечне материалов, подлежащим засекречиванию по линии Раскартографии (сейчас Росреестра) говорится о государственных системах координат, а WGS-84 никогда государственной в РФ не считалась... и это требование распространяется на пункты государственной геодезической сети (даже не на пункты ОМС), но не на точки границ земельных участков (больше половины из которых не закреплены на местности)...


Большие замечания:
1. Мировая система координат засекречиванию подлежать не может по той причине - что это тогда не будет не то что "мировая" - а вовсе не система координат!
2. "Государственная" система координат может быть как мировой так и местной. Мировые системы координат могут быть как "государственными" так и неофициальными. И их кстати, насчитываются многие десятки.
3. Ваш "перечень материалов" не имеет отношения к тому о чем я говорил. Я говорил о том что гостайной является точное определение координат конкретной точки в ЛЮБОЙ из мировых систем координат. А не сама система координат.
Цитата из перечня сведений, подлежащих засекречиванию:
15. Сведения о местонахождении (координаты) геодезических пунктов и географических объектов, определенные с точностью 30 метров и точнее в государственной системе координат 1942 года и в геоцентрических системах координат, а также геодезические, картографические и другие материалы, позволяющие вычислить или уточнить эти координаты, кроме астрономических, астрофизических, геофизических обсерваторий и станций, других научных и измерительных пунктов, средств и объектов, включенных в международные сети и системы на основании международных договоров и соглашений Российской Федерации по спискам, согласованным с Генеральным штабом Вооруженных Сил Российской Федерации - секретно.
16. Сведения, содержащие координаты режимного объекта, независимо от их точности и системы координат - секретно.
17. Сведения о местонахождении (координаты) пунктов геодезических сетей на территории режимных объектов, независимо от системы координат и площади этой территории - секретно.
18. Каталоги и списки координат геодезических пунктов и географических объектов в местных системах координат и системе координат 1963 года, составленные на участки местности, превышающие территорию одного административного района субъекта Российской Федерации или площадью свыше 5000 кв. км - секретно.
19. Исходные данные (ключи) местных систем координат, системы координат 1963 года и материалы, раскрывающие переход от этих систем координат к государственной системе координат 1942 года - секретно.

Итак... п.15 говорит как о "геодезических пунктах", так и о и "географических объектах" - которыми являются не только земельные участки, но и любые, однозначно привязанные к территории объекты - ручьи, колодцы, здания, мертвецки пьяные геодезисты...

Выделенное "и в геоцентрических системах координат" - как раз и означает в "мировых" в нашем контексте. Я могу отказаться от слова "мировых", чтобы не спорить о терминологии, но упомянутая система WGS-84 в любом случае является именно геоцентрической.

И наконец о "незакрепленных на местности" точках границ земельных участков. В перечне не указывается что секретными могут являться только "закрепленные". Будет очень прикольно, если мой участок станет секретным сразу как только я его обнесу забором... а когда меня арестуют я скажу жене, чтобы сосед сломал забор и тогда меня выпустят из кутузки... 8 ))
PavelML
Заслуженный тролль ресурса
 
Сообщения: 104
Зарегистрирован: 20 фев 2010, 17:29
Благодарил (а): 0 раз.
Поблагодарили: 6 раз.

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

Сообщение vasketsov » 31 мар 2011, 00:29

PavelML писал(а):господа, будете смеяться, но я после бесплодных мучений...
... получил... векторный прямоугольник запрошенного окна в WGS-84 и с координатами точностью в 13 знаков после запятой!
А теперь буду пытаться скрестить ежика с ужиком...

Смеяться конечно вряд ли, но на всякий случай напишу вот что, может будет полезно.
Итак, идём сюда и крутим вниз до конца страницы, видим экспорт карты и тыкаемся туда.
Вылазит примерно стандартное для этого портала табло для ввода информации для поиска.
Там пишем что-нибудь типа 4158174.33871359,7528741.53797672,4163066.30852384,7523849.56816647 в адрес.
А теперь самое интересное.
Ниже указываем image format = svg и format = image.
Жмакаем на экспорт имидж и получаем файло с именем export, внутри которого сидит, на минутку, честный svg квадрата (конкретно этот квадрат - северозапад мск-область на 14-м зуме в сасе).
Открывается, например, в опере и ещё много где, в педивикии по svg дофига ссылок кому надо - можно дальше рыть в этом направлении, например, в направлении превращения его в kml для саса, чтобы типа слоя вики отображалось.
Можно ли вернуть svgz - не знаю, может быть и можно, теоретически портал умеет выдавать в этом формате данные, но в списке его нету, комбинация svg+kmz ни к чему хорошему не привела, приходит файл с границами на 484 байта.
vasketsov
Специалист
 
Сообщения: 901
Зарегистрирован: 25 июл 2009, 21:15
Благодарил (а): 0 раз.
Поблагодарили: 198 раз.

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

Сообщение tae1980 » 31 мар 2011, 08:50

PavelML писал(а):Лично я конвертирую координаты своими функциями, базирующиммися на ГОСТ 51794-2001. Который нужно перед использованием весьма творчески переработать, поскольку датум туда вкрячен насмерть и он "типовой" - то есть ни на что непригодный.

А можно на это взглянуть на функции? Хотя бы в качестве учебного материала.
PavelML писал(а):К сожалению для кадастровых целей неприемлемо, там нужно иметь еще пару порядков

Дополнительные догон "ручками" еще ни кто не отменял, особенно если есть точные координаты пары кварталов или земельных участков.

Вы видимо невнимательно прочли. Какой "догон" при ОГРУБЛЕНИИ НЕИЗВЕСТНЫХ координат на пару порядков?
Если бы у меня были точные координаты земельных участков - а тогда зачем мне было бы их скачивать? Гы...[/quote]
зная координаты одного земельного участка, можно закоординировать весь квартал. Точных измерений проводить конечно нельзя, но факт что рядом стоит уже учтенный земельный участок, и было бы не плохо заказать выписку на него для избежания приостановки из-за наложения - можно.
tae1980
Новичок
 
Сообщения: 21
Зарегистрирован: 29 окт 2010, 13:17
Благодарил (а): 0 раз.
Поблагодарили: 0 раз.

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

Сообщение Adon73 » 31 мар 2011, 08:51

К сожалению, вы ссылаетесь на устаревший перечень материалов... , см.
"Перечень сведений, подлежащий засекречиванию. Утвержден приказом Министерства транспорта РФ от 10 ноября 2007 г. (с изменениями от 26.05.2008)"
(который сам является ДСП).
Adon73
Новичок
 
Сообщения: 49
Зарегистрирован: 10 июл 2010, 16:22
Благодарил (а): 0 раз.
Поблагодарили: 7 раз.

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

Сообщение PavelML » 31 мар 2011, 14:25

vasketsov писал(а):...Итак, идём сюда и крутим вниз до конца страницы, видим экспорт карты и тыкаемся туда...


Да, я тоже туда забегал... не хватает времени детально все изучить, надо копать дальше...
PavelML
Заслуженный тролль ресурса
 
Сообщения: 104
Зарегистрирован: 20 фев 2010, 17:29
Благодарил (а): 0 раз.
Поблагодарили: 6 раз.

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

Сообщение PavelML » 31 мар 2011, 14:40

tae1980 писал(а):
PavelML писал(а):Лично я конвертирую координаты своими функциями, базирующиммися на ГОСТ 51794-2001. Который нужно перед использованием весьма творчески переработать, поскольку датум туда вкрячен насмерть и он "типовой" - то есть ни на что непригодный.

А можно на это взглянуть на функции? Хотя бы в качестве учебного материала.

Можно. Но функции на васике. Предупреждаю, что выкладываю функции с неточным датумом:
Код: Выделить всё
'Набор функций для преобразований координат на основе ГОСТ Р 51794-2001
'версия от 01.05.2009, нуждается в подгонке смещений по X и Y
Const Pi As Double = 3.14159265358979 ' Число Пи
Const ro As Double = 206264.8062 ' Число угловых секунд в радиане

' Эллипсоид Красовского
Const aP As Double = 6378245 ' Большая полуось
Const alP As Double = 1 / 298.3 ' Сжатие
Const e2P As Double = 2 * alP - alP ^ 2 ' Квадрат эксцентриситета

' Эллипсоид WGS84 (GRS80, эти два эллипсоида сходны по большинству параметров)
Const aW As Double = 6378137 ' Большая полуось
Const alW As Double = 1 / 298.257223563 ' Сжатие
Const e2W As Double = 2 * alW - alW ^ 2 ' Квадрат эксцентриситета

' Вспомогательные значения для преобразования эллипсоидов
Const a As Double = (aP + aW) / 2
Const e2 As Double = (e2P + e2W) / 2
Const da As Double = aW - aP
Const de2 As Double = e2W - e2P

' Линейные элементы трансформирования, в метрах
'Const dx As Double = 20 '20
'Const dy As Double = -141 '-141
'Const dz As Double = -80 '-80
' Угловые элементы трансформирования, в секундах
'Const wx As Double = 0 '0
'Const wy As Double = -0.35 '-0.35
'Const wz As Double = -0.82 '-0.82
' Дифференциальное различие масштабов
'Const ms As Double = 0 '0

' Линейные элементы трансформирования, в метрах
Public dx As Double '20
Public dy As Double '-141
Public dz As Double '-80
' Угловые элементы трансформирования, в секундах
Public wx As Double '0
Public wy As Double '-0.35
Public wz As Double '-0.82
' Дифференциальное различие масштабов
Public ms As Double  '0

Public Function XY13_NW(X As Double, Y As Double, Optional G As Byte = 0, Optional H As Double = 200) As Double
'вычисление широты WGS84 по плоским прямоугольным координатам в системе СК-13 методом последовательного приближения, третий параметр = 1 для гугловской подложки
If IsNull(X) Or IsNull(Y) Then XY13_NW = 0: Exit Function 'данные недостаточны функция отменена
If X < 300000 Or X > 510000 Or Y < 1000000 Or Y > 1400000 Then XY13_NW = 0: Exit Function 'данные некорректны функция отменена
Dim N As Double: Dim E As Double: Dim X0 As Double: Dim Y0 As Double: Dim Cnt As Long
E = (1643130.458 + Y) * 0.000015398: N = (5639557.797 + X) * 0.0000089838: X0 = 0: Y0 = 0: Cnt = 0
Do Until (Abs(X - X0) < 0.0000001) And (Abs(Y - Y0) < 0.0000001)
X0 = W13_X(N, E, G, H): Y0 = W13_Y(N, E, G, H): Cnt = Cnt + 1
N = N + (X - X0) * 0.0000089838: E = E + (Y - Y0) * 0.000015398
Loop
XY13_NW = N
End Function

Public Function XY13_EW(X As Double, Y As Double, Optional G As Byte = 0, Optional H As Double = 200) As Double
'вычисление долготы WGS84 по плоским прямоугольным координатам в системе СК-13 методом последовательного приближения, третий параметр = 1 для гугловской подложки
If IsNull(X) Or IsNull(Y) Then XY13_EW = 0: Exit Function 'данные недостаточны функция отменена
If X < 300000 Or X > 510000 Or Y < 1000000 Or Y > 1400000 Then XY13_EW = 0: Exit Function 'данные некорректны функция отменена
Dim N As Double: Dim E As Double: Dim X0 As Double: Dim Y0 As Double: Dim Cnt As Long
E = (1643130.458 + Y) * 0.000015398: N = (5639557.797 + X) * 0.0000089838: X0 = 0: Y0 = 0: Cnt = 0
Do Until (Abs(X - X0) < 0.0000001) And (Abs(Y - Y0) < 0.0000001)
X0 = W13_X(N, E, G, H): Y0 = W13_Y(N, E, G, H): Cnt = Cnt + 1
N = N + (X - X0) * 0.0000089838: E = E + (Y - Y0) * 0.000015398
Loop
XY13_EW = E
End Function

Public Function W13_X(NW As Double, EW As Double, Optional G As Byte = 0, Optional HW As Double = 200) As Double
'параграф ГОСТ 4.3. Преобразование геодезических координат в плоские прямоугольные
'вычисление координаты X в системе СК-13 по геодезическим координатам WGS84 (Пулково-42 при выключении вызова функции) по шести-(трех-)градусной зоне
' Линейные элементы трансформирования, в метрах
dx = 19.9555 '20
dy = -139.673 '-141
dz = -77 '-80
' Угловые элементы трансформирования, в секундах
wx = 0 '0
wy = -0.35 '-0.35
wz = -0.82 '-0.82
' Дифференциальное различие масштабов
ms = -0.00000012 '0
If IsNull(NW) Or IsNull(EW) Then W13_X = 0: Exit Function 'данные недостаточны функция отменена
If NW < 53.5 Or NW > 55.5 Or EW < 42 Or EW > 47 Then W13_X = 0: Exit Function 'данные некорректны функция отменена
Dim NR As Double: Dim ER As Double: Dim A1 As Double: Dim A2 As Double: Dim A3 As Double: Dim A4 As Double 'геодезические координаты широта и расстояние от центрального меридиана, выраженные в радианах (градусы, деленные на 57.2957795130824)
Dim Cos1NR As Double: Dim Sin1NR As Double: Dim Sin2NR As Double: Dim Sin4NR As Double: Dim Sin6NR As Double 'вычисленные косинус и синусы широты в разной степени
Dim Sin1NR2 As Double: Dim ER2 As Double: Dim SmZonY As Double: Dim SmZonX As Double: Dim CMer As Double 'вычисленные синус удвоенной широты и квадрат расстояния широты от центрального меридиана, радиан в квадрате
N = WGS84_SK42_N(NW, EW, HW)
E = WGS84_SK42_E(NW, EW, HW)
CMer = 44.55: SmZonX = -5614743.5: SmZonY = 1250000 ': If E < 43.05 Then CMer = 41.55
NR = N / 57.29577951: ER = (E - CMer) / 57.29577951: ER2 = ER * ER
Cos1NR = Cos(NR): Sin1NR = Sin(NR): Sin2NR = Sin1NR * Sin1NR: Sin4NR = Sin2NR * Sin2NR: Sin6NR = Sin4NR * Sin2NR: Sin1NR2 = Sin(2 * NR)
A1 = ER2 * (109500 - (574700 * Sin2NR) + (863700 * Sin4NR) - (398600 * Sin6NR))
A2 = ER2 * (278194 - (830174 * Sin2NR) + (572434 * Sin4NR) - (16010 * Sin6NR) + A1)
A3 = ER2 * (672483.4 - (811219.9 * Sin2NR) + (5420 * Sin4NR) - (10.6 * Sin6NR) + A2)
A4 = ER2 * (1594561.25 + (5336.535 * Sin2NR) + (26.79 * Sin4NR) + (0.149 * Sin6NR) + A3)
W13_X = (6367558.4968 * NR) - (Sin1NR2 * (16002.89 + (66.9607 * Sin2NR) + (0.3515 * Sin4NR) - A4)) + SmZonX
End Function

Public Function W13_Y(NW As Double, EW As Double, Optional G As Byte = 0, Optional HW As Double = 200) As Double
'параграф ГОСТ 4.3. Преобразование геодезических координат в плоские прямоугольные
'вычисление координаты Y в системе СК-13 по геодезическим координатам WGS84 (Пулково-42 при выключении вызова функции) по шести-(трех-)градусной зоне
' Линейные элементы трансформирования, в метрах
dx = 19.9555 '20
dy = -139.673 '-141
dz = -77 '-80
' Угловые элементы трансформирования, в секундах
wx = 0 '0
wy = -0.35 '-0.35
wz = -0.82 '-0.82
' Дифференциальное различие масштабов
ms = -0.00000012 '0
If IsNull(NW) Or IsNull(EW) Then W13_Y = 0: Exit Function 'данные недостаточны функция отменена
If NW < 53.5 Or NW > 55.5 Or EW < 42 Or EW > 47 Then W13_Y = 0: Exit Function 'данные некорректны функция отменена
Dim NR As Double: Dim ER As Double: Dim A1 As Double: Dim A2 As Double: Dim A3 As Double: Dim A4 As Double 'геодезические координаты широта и расстояние от центрального меридиана, выраженные в радианах (градусы, деленные на 57.2957795130824)
Dim Cos1NR As Double: Dim Sin1NR As Double: Dim Sin2NR As Double: Dim Sin4NR As Double: Dim Sin6NR As Double 'вычисленные косинус и синусы широты в разной степени
Dim Sin1NR2 As Double: Dim ER2 As Double: Dim SmZonY As Double: Dim SmZonX As Double: Dim CMer As Double 'вычисленные синус удвоенной широты и квадрат расстояния широты от центрального меридиана, радиан в квадрате
N = WGS84_SK42_N(NW, EW, HW)
E = WGS84_SK42_E(NW, EW, HW)
CMer = 44.55: SmZonX = -5614743.5: SmZonY = 1250000 ': If E < 43.05 Then CMer = 41.55
NR = N / 57.29577951: ER = (E - CMer) / 57.29577951: ER2 = ER * ER
Cos1NR = Cos(NR): Sin1NR = Sin(NR): Sin2NR = Sin1NR * Sin1NR: Sin4NR = Sin2NR * Sin2NR: Sin6NR = Sin4NR * Sin2NR: Sin1NR2 = Sin(2 * NR)
A1 = ER2 * (79690 - (866190 * Sin2NR) + (1730360 * Sin4NR) - (945460 * Sin6NR))
A2 = ER2 * (270806 - (1523417 * Sin2NR) + (1327645 * Sin4NR) - (21701 * Sin6NR) + A1)
A3 = ER2 * (1070204.16 - (2136826.66 * Sin2NR) + (17.98 * Sin4NR) - (11.99 * Sin6NR) + A2)
W13_Y = (ER * Cos1NR * (6378245 + (21346.1415 * Sin2NR) + (107.159 * Sin4NR) + (0.5977 * Sin6NR) + A3)) + SmZonY
End Function
Public Function WGS84_SK42_N(Bd, Ld, H As Double) As Double
WGS84_SK42_N = Bd - dN(Bd, Ld, H) / 3600
End Function

Public Function SK42_WGS84_N(Bd, Ld, H As Double) As Double
SK42_WGS84_N = Bd + dN(Bd, Ld, H) / 3600
End Function

Public Function WGS84_SK42_E(Bd, Ld, H As Double) As Double
WGS84_SK42_E = Ld - dE(Bd, Ld, H) / 3600
End Function

Public Function SK42_WGS84_E(Bd, Ld, H As Double) As Double
SK42_WGS84_E = Ld + dE(Bd, Ld, H) / 3600
End Function

Function dN(Bd, Ld, H As Double) As Double
Dim B, L, M, N As Double
B = Bd * Pi / 180
L = Ld * Pi / 180
M = a * (1 - e2) / (1 - e2 * Sin(B) ^ 2) ^ 1.5
N = a * (1 - e2 * Sin(B) ^ 2) ^ -0.5
dN = ro / (M + H) * (N / a * e2 * Sin(B) * Cos(B) * da + (N ^ 2 / a ^ 2 + 1) * N * Sin(B) * Cos(B) * de2 / 2 - (dx * Cos(L) + dy * Sin(L)) * Sin(B) + dz * Cos(B)) - wx * Sin(L) * (1 + e2 * Cos(2 * B)) + wy * Cos(L) * (1 + e2 * Cos(2 * B)) - ro * ms * e2 * Sin(B) * Cos(B)
End Function

Function dE(Bd, Ld, H As Double) As Double
Dim B, L, N As Double
B = Bd * Pi / 180
L = Ld * Pi / 180
N = a * (1 - e2 * Sin(B) ^ 2) ^ -0.5
dE = ro / ((N + H) * Cos(B)) * (-dx * Sin(L) + dy * Cos(L)) + Tan(B) * (1 - e2) * (wx * Cos(L) + wy * Sin(L)) - wz
End Function

Public Function WGS84_H(Bd, Ld, H As Double) As Double
Dim B, L, N, dH As Double
B = Bd * Pi / 180
L = Ld * Pi / 180
N = a * (1 - e2 * Sin(B) ^ 2) ^ -0.5
dH = -a / N * da + N * Sin(B) ^ 2 * de2 / 2 + (dx * Cos(L) + dy * Sin(L)) * Cos(B) + dz * Sin(B) - N * e2 * Sin(B) * Cos(B) * (wx / ro * Sin(L) - wy / ro * Cos(L)) + (a ^ 2 / N + H) * ms
WGS84_H = H + dH
End Function


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

Мы с Вами разговариваем о разных вещах, Вы не заметили? Меня не интересует проблема геодезической съемки, ее для меня не существует. Меня интересует получение КАДАСТРОВЫХ КООРДИНАТ ИЗ ПУБЛИЧНОЙ КАДАСТРОВОЙ КАРТЫ. Ну да ладно...
PavelML
Заслуженный тролль ресурса
 
Сообщения: 104
Зарегистрирован: 20 фев 2010, 17:29
Благодарил (а): 0 раз.
Поблагодарили: 6 раз.

Пред.След.

Вернуться в РосКосмос, ГГЦ, РосРеестр, КосмоСнимки

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и гости: 3