МАГИСТЕРСКАЯ РАБОТА

Информационные технологии обеспечения пространственной
привязки объектов кадастровой съемки

Павельчак Татьяна Валентиновна

Руководитель: д.т.н., проф. Могильный Сергей Георгиевич


Как показывает опыт наиболее развитых стран, одним из ключевых моментов эффективной экономической политики является научно обоснованное землепользование, средством достижения которого является земельная реформа. Для эффективного управления земельными ресурсами и развития земельного рынка необходимо располагать достоверной и надежной информацией о земле. Ядром этой системы должен быть земельный кадастр, как методически упорядоченный государственный учет сведений об участках земли с уникальными интересами владения. В последнее время в нашей стране создана система областных и городских кадастровых центров, основной задачей которых является создание полноценной и достоверной кадастровой системы.

Препятствием для эффективной работы в этой области является проблема однозначной пространственной привязки объектов землепользования и недвижимости, то есть определение или восстановление в пространстве границ прав и интересов владельцев.

Для каждого такого участка землепользования проводится кадастровая съемка. Традиционно она происходит по такой схеме:

построение опорной геодезической сети ®

® построение съемочной геодезической сети ®

® координирование точек поворота границ объекта кадастровой съемки.

Это есть так называемый “принцип от общего к частному”, при котором точность измерений снижается от этапа к этапу съемки (от измерения съемочного обоснования к определению координат точек поворота границ объекта). Тем самым снижаются трудозатраты, а, следовательно, и денежные расходы.

Попытаемся выделить недостатки вышеупомянутого принципа при проведении кадастровых работ. Развитие, сохранение и восстановление пунктов геодезического обоснования достаточной плотности – очень дорогое и сложное дело. А в большинстве случаях на сегодняшний день пункты геодезического обоснования нарушены в силу различных причин и не отвечают требованиям инструкции. В условиях недостаточного финансирования геодезической отрасли нашей страны эта задача вообще является проблематичной. На практике очень часто, в зависимости от потребностей заказчиков, проводится не сплошная кадастровая съемка, а выборочная съемка отдельных участков землепользования. В таких условиях неэффективно развитие дорогого геодезического обоснования. Поэтому координаты объекта кадастровой съемки часто определяют в условной системе, где основное внимание уделяется взаимному положению точек поворота границ этого объекта. Предполагается, что положение данного участка землепользования в общепринятой системе координат будет определено в последующем, по мере поступления средств и при дальнейшей съемке соседних участков. Такой подход весьма неэффективен.

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

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

Можно назвать следующие современные средства измерений:

Эти средства позволяют нам практически одинаково точно (1-2 см) определить как точки хода обоснования, так и отдельные пикетные точки, что не соответствует “принципу от общего к частному”.

В данной работе предлагается нетрадиционное решение вышеописанных проблем пространственной привязки участков кадастровой съемки.

Чтобы обосновать предлагаемый метод, нужно выделить два предположения, при которых он имеет смысл:

  1. Измерения производятся современными приборами с точностью порядка 1-2 см;
  2. Точки поворота границ землепользования не являются геометрическими точками, а есть некоторые тела или части тел, и их идентификация на местности при измерениях производится с погрешностью, значительно превышающей точность измерений и достигающей 5-10 см.

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

Сущность предлагаемой схемы пространственной привязки показана на рисунке 1 .

Рисунок 1

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

Для опис ;ания теории ура внивания таких сетей необходи мо представить ситуацию: имеют ся четыре участ ка i, j, m, z, коо рдинаты точек п оворота границ участков определены в разных системах, имеются опорные точки.

Пусть для точки А координаты : (ХjА,YjА), (ХiА,YiА), (ХmА,YmА), (ХzА,YzА).

Рисунок 2

Из-за погрешности измерений эти точки могут не совпадать. Тогда возникают перекрытия площадей участков (между участками i и j) или пробелы (между участками j и m, i и z). Таким образом, возникает задача перевычислить координаты в одну систему, и уравновесить их значения, а следовательно, сбалансировать площади участков. Необходимо также найти, какие координаты точек определены с грубой ошибкой.

Рассмотрим решение данной задачи.

Для перехода из одной системы координат в другую известны соотношения для точки :

Xk =Xo +t(Xn cosθ+Yn sinθ); (1)
Yk =Yo +t(Yn cosθ-Xn sinθ),

где Xo , Yo - начало координат системы n относительно системы k ;

t-масштаб;

θ- угол поворота системы n относительно системы k .

Используем данное соотношение для точки А:

(2)

X(j)A = Xjo +aj XjA+ bj YjA;

Y(j)A = Yjo -bj XjA+ aj YjA,

X(i)A = Xio +ai XiA+ bi YiA;

Y(i)A = Yio –bi XiA+ ai YiA,

X(m)A = Xmo +am XmA+ bm YmA;

Y(m)A = Ymo –bm XmA+ am YmA.

X(z)A = Xzo +az XzA+ bz YzA;

Y(z)A = Yzo –bz XzA+ az YzA.

Неизвестные параметры для точки А : Xjo , Yjo, aj, bj,

Xio , Yio, ai, bi,

Xzo , Yzo, az, bz,

Xmo , Ymo, am, bm.

Если бы координаты точки А были определены без о шибок, то:

X(j)A = X(i)A = X(i)A = X(z)A ;

Y(j)A = Y(i)A = Y(i)A = Y(z)A.

Но фактически возникают невязки:

(3)

X(j)A -X(i)A = ν XA(i j);

Y(j)A -Y(i)A = νYA(i j),

X(j)A -X(z)A XA(j z);

Y(j)A -Y(z)A = νYA(j z),

X(j)A -X(m)A XA(j m);

Y(j)A -Y(m)A = νYA(j m).

Если точка А является общей для n участков, то таких систем можно составить (n-1). Если точка А является опорной и значения ее координат определены в собственных системах координат n смежных участков, то можно составить n систем вида:

X(q)A –XopA XA(q); (4)
Y(q)A –YopA = νYA(q),

где, q=1..n..

Запишем окончательный вид уравнений. Возьмем точку В, которая является общей для участка i и участка j. Для точки В можно записать следующие параметрические уравнения:

(5)

Запишем эти уравнения через приближенные значения параметров и поправки к ним:

(6)

В матричной форме:

- матрица коэффициентов уравнений поправок;

- вектор поправок в приближенные значения параметров;

- вектор приближенных значений параметров;

- невязки в уравнениях;

где, t - номер участка, К- имя точки.

Для точки В:

(7)

Можно составить такие уравнения для каждой рассматриваемой точки. Получится общая система уравнений в общем виде:

ν=A*ΔT+L. (8)

Данную систему можно решить двумя способами:

  • -методом наименьших квадратов [2] при условии
  • [ν x νx+ ν yν y]=min (метод Гаусса); (9)

  • - методом наименьших модулей [3] при условии
  • [|νx |+ |νy|]=min (метод Лапласа). (10)

    Решая систему (8) методом наименьших квадратов при условии (9), получаем систему нормальных уравнений, которую в матричной форме можно записать [2]:

    NΔT+ATPL=0, (11)

    где, N=ATPA (12) - матрица коэффициентов нормальных уравнений,

    Р - матрица весов.

    Матрица весов принимается единичной.

    Полученную систему нормальных уравнений решаем итерационным способом. Первоначально приближенные значения параметров берем равными нулю.

    При решении системы (8) методом наименьших модулей при условии (10) используется вариационный метод решения [3]. Составляется такая же система нормальных уравнений, как при использовании метода наименьших квадратов. Матрица весов при первой итерации принимается единичной. При первой итерации система решается методом наименьших квадратов. По полученным параметрам вычисляются невязки ν(0) в каждом уравнении. Составляется матрица весов 2x2 для каждой системы вида (7), где диагональные элементы – веса уравнения, которые принимаются равными:

    ; (13)
    ,

    где, i - номер или имя пары точек, для которых составлены параметрические уравнения вида (7);

    j - номер итерации,

    µ - средняя невязка по системе уравнений.

    Полученная матрица весов и поправки к параметрам d Xto, d Yto, d at, d bt используются в следующей итерации. Решается система нормальных уравнений для итерации (j+1):

    (14)

    Находят матрицу весов и поправки, но уже к следующей итерации. Итерации продолжаются до тех пор, пока поправки к параметрам a и b не будут меньше допустимой величины.

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

    В результате получаем неизвестные параметры Xo , Yo, a, b для каждой точки, по которым можно вычислить уравненные координаты в одной общей системе координат. Осталось привести значения координат общих точек в разных участках к единым значениям.

    Это могут быть средние арифметические значения из уравненных координат для конкретной точки. При использовании метода наименьших квадратов:

    ; (15)
    ,

    где, n-количество участков, на которых расположена эта точка.

    Если же использовался метод наименьших модулей, то определяется такие координаты и из всех уравненных значений, при которых:

    (16)

    где, n-количество участков, на которых расположена эта точка.

    и , удовлетворяющие условию (15), и будут необходимыми значениями координат. Если количество участков n – четное, то может возникнуть ситуация, что условию (15) соответствуют несколько значений , тогда берется среднее из них.

    Согласно выше описанной теории, для пространственной привязки большого количества соседних участков требуется оперировать с матрицами большого порядка. Современная компьютерная техника позволяет достаточно легко решить такую задачу. Разработана программа AreaEqual, которая позволяет выполнять уравнивание сетей и приводить координаты съемок отдельных участков в единую систему. Если нет данных об опорных точках, то программа перевычисляет все полигоны в общую для всех локальную систему координат.

    Исходными данными в программе для каждого участка землепользования служат:

    • локальные координаты точек, расположенных при обходе полигона по часовой стрелке;
    • описание соседства: имя соседнего полигона, номера начальной и конечной точек общей границы.
    • координаты опорных точек.

    Программа AreaEqual может быть представлена в виде выполнения следующих этапов:

      Согласно выше описанной теории, для пространственной привязки большого количества соседних участков требуется оперировать с матрицами большого порядка. Современная компьютерная техника позволяет достаточно легко решить такую задачу. Разработана программа AreaEqual, которая позволяет выполнять уравнивание сетей и приводить координаты съемок отдельных участков в единую систему. Если нет данных об опорных точках, то программа перевычисляет все полигоны в общую для всех локальную систему координат.

      Исходными данными в программе для каждого участка землепользования служат:

      • локальные координаты точек, расположенных при обходе полигона по часовой стрелке;
      • описание соседства: имя соседнего полигона, номера начальной и конечной точек общей границы.
      • координаты опорных точек.

      Программа AreaEqual может быть представлена в виде выполнения следующих этапов:

      1. Подготовка исходных данных к обработке.
      2. Нужно сказать, что данный модуль может дополняться новыми процедурами в зависимости от того, в каком виде представлены исходные данные, и какими методами они были получены. Например, если данные были получены в результате тахеосъемки точек поворота границ участка, то необходимо составить процедуру, которая бы вычисляла координаты углов границ участка в локальной системе координат участка.

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

      3. Расположение участков в таком порядке, при котором ресурсов оперативной памяти во время обработки для хранения данных отводится наименьшее количество.
      4. Эта задача возникает из-за того, что программа работает с большим количеством неизвестных и порядок нормальной матрицы также велик. Конкретный участок связан соседством не со всеми участками, только с некоторыми. Таким образом, полученная система нормальных уравнений является слабо заполненной (большинство коэффициентов равны “нулю”).

        Необходимо найти такую рациональную последовательность участков, чтобы матрица коэффициентов нормальных уравнений имела вид, показанный на рисунке 3:

        Рисунок 3

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

      5. Заполнение системы нормальных уравнений данными.

      Система нормальных уравнений разбивается на субматрицы, такая субматрица выступает как единица хранения. Это можно показать на следующем примере.

      Возьмем точку В, которая является общей для участка i и участка j.

      Для точки В:

      Заполняются матрицы AiB, и AjB, по исходным данным. Согласно N=ATPA вычисляются субматрицы для этих участков и заносятся в матрицу коэффициентов системы нормальных уравнений по адресам, которые соответствуют участкам i и j. Таким же образом вычисляется вектор свободных членов системы нормальных уравнений, согласно ATPL.

      Заполнение субматриц происходит в зависимости от выбранного метода решения (МНК или МНМ) согласно вышеописанной теории.

      1. Решение систем нормальных уравнений.

      Система нормальных уравнений уравнивается методом Гаусса.

      1. Перевычисление координат согласно новым параметрам.
      2. Определение единых координат общих точек поворота границ соседних участков.
      3. Определение сбалансированных площадей участков землепользования по окончательным координатам.

      Данную программу можно применять в следующих целях:

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

      2) Допустим, при инвентаризации участков землепользования измерения производятся в собственных системах координат для каждого из них. Причем в качестве начала координат (0,0) в участке была взята любая удобная точка, а ось Х направлена в сторону одной из

      вершин участка. Назовем эту систему координат “локальной”. Измерения проводились, например, электронным тахеометром (или другим способом). Можно легко вычислить координаты углов поворота границ участка в локальной системе координат по полученным измерениям, решая прямую геодезическую задачу. Достаточно иметь на территории инвентаризируемого района 3-4 опорных точки с координатами в государственной системе и не прокладывать опорную геодезическую сеть для определения значений координат углов поворотов границ участков в единой системе. Таким образом, в зависимости от того, каким техническим средством и методом были собраны исходные данные, нужно дополнить программу модулем. Этот модуль должен из данных измерений вычислять координаты поворотных точек участков землепользования (их можно рассматривать как приближенные координаты) и приводить полученные координаты в необходимый формат и структуру данных для работы программы. Имея опорные точки и координаты углов поворота в локальной системе координат, можно получить сбалансированные координаты в государственной или местной системе координат. А по ним можно определить площади участков землепользования.

      3) Пусть имеются значения координат углов поворота границ для системы участков. Некоторые из координат определены с грубой ошибкой. Если воспользоваться методом наименьших модулей в программе, то можно локализовать эту ошибку. Если же воспользоваться методом наименьших квадратов, то ошибка распределится на все соседние участки и исказит координаты соседних точек. Таким образом, программа позволяет успешно решать задачу отбраковки грубых измерений.

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

      1. Сеть состоит из 25 полигонов, образующих квадраты. В четырех углах общего контура расположены опорные точки (рисунок 4).

      Рисунок 4

      Методом статистического моделирования исследовано накопление ошибок в сети. В локальные координаты каждой точки модели вводились случайные ошибки, распределенные по нормальному закону с математическим ожиданием, равным нулю, и стандартом s =0.05м. Модель такой сети уравнивалась по программе AreaEqual, вычислялись оставшиеся истинные ошибки координат точек. По 100 подобным экспериментам найдены средние линейные ошибки ml положения каждой точки. Величины ml изображены в виде поверхности в перспективной проекции на рисунке 5.

      Рисунок 5 Видно, что наибольшие ошибки имеют граничные точки. Внутри контура ошибки находятся примерно на уровне 2s и достигают ± 10 см.

      Статистические свойства погрешностей локальных координат таковы, что метод наименьших квадратов (МНК) не всегда обеспечивает оптимальное уравнивание, поэтому в программе AreaEqual предусмотрена возможность решения системы нормальных уравнений методом наименьших модулей Лапласа (МНМ). МНМ позволяет эффективно выявить и локализовать грубые ошибки в исходных данных

      Наглядная демонстрация вышесказанного представлена на рисунке 4. В координаты точки A' искусственно вводится грубые ошибки. В целях эксперимента полученная модель уравнивается методом наименьших квадратов и методом наименьших модулей. На рисунке 4 четко видно, как исказились координаты точек всех соседних полигонов при использовании МНК. При использовании МНМ модель не исказилась, проявились лишь грубые ошибки lоординат точки A' .

       

      2.Теперь продемонстрируем работу программы на реальном объекте. Взят контур реально существующего села, состоящий из пяти кварталов и 114 участков землепользования ( рисунок 6), занимающих площадь » 45га. Известны локальные координаты каждого участка землепользования. Все кварталы можно легко увязать между собой. Для этого в исходные данные одного из участков n квартала А вносятся координаты угла близлежащего участка g, расположенного в соседнем квартале В, в системе координат участка n. Эти данные одновременно получают при съемке участка n.

      В локальные координаты каждой точки участков вводятся случайные ошибки, распределенные по нормальному закону с математическим ожиданием, равным нулю, и стандартом s =0.05м.

      Рисунок 6

      При наличии четырех опорных точек на все пять кварталов, расположенных по краям общего контура, линейные ошибки положения каждой точки накапливаются в пределе от ± 0.05м до ± 0.45м. Такой точности не достаточно для практических целей. Следовательно, необходимо повысить количество опорных точек.

      Накопление линейных ошибок координат менее существенно при 8-9 опорных точках . Накопление ошибок зависит от расположения и количества опорных точек. Желательно иметь опорную точку внутри контура и в его слабых местах на границе. При девяти опорных точках накопление линейных ошибок в большей части сети не превышает ± 0.20м при погрешностях измерения координат ± 0.05м. Таким образом, снижение накопления ошибок возможно за счет повышения количества опорных точек и точности определения координат. Целью дальнейшего исследования является более подробная разработка рекомендаций по оптимизации геометрических конструкций сетей пространственной привязки кадастровых съемок.


      Список литературы:
      1. Могильный С.Г., Фотограмметрия. – К.: 1985.
      2. Папазов М.Г., Могильный С.Г., Теория ошибок и способ наименьших квадратов. – М.: 1968.
      3. В.И. Мудров, В.Л. Кушко, Методы обработки измерений. Квазиправдоподобные оценки. – М.: 1983.




      Главная страница