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

Автор работы: Пользователь скрыл имя, 12 Марта 2015 в 13:35, курсовая работа

Краткое описание

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

Содержание

ВВЕДЕНИЕ…………………………………………………………………………3
ГЛАВА 1. ВИДЫ ГЕОДЕЗИЧЕСКИХ РАБОТ, ВЫПОЛНЯЕМЫХ ДЛЯ ВЕДЕНИЯ ЗЕМЕЛЬНОГО КАДАСТРА………………………………………….5
1.1 МЕЖЕВАНИЕ ЗЕМЕЛЬ………………………………………………….……..6
1.2 СОЗДАНИЕ ОПОРНОЙ МЕЖЕВОЙ СЕТИ…………………………….……8
1.3 ТОПОГРАФИЧЕСКИЕ СЪЕМКИ………………………………………..…..11
1.4 СОЗДАНИЕ ЦИФРОВЫХ КАРТ………………………………………..……15

ГЛАВА 2.НОРМАТИВНО - МЕТОДИЧЕСКАЯ БАЗА ВЫПОЛНЕНИЯ ГЕОДЕЗИЧЕСКИХ РАБОТ ДЛЯ ВЕДЕНИЯ ЗЕМЕЛЬНОГО КАДАСТРА…...16
2.1 СОСТАВЛЕНИЕ КАРТЫ (ПЛАНА) ОБЪЕКТА ЗЕМЛЕУСТРОЙСТВА ИЛИ КАРТЫ (ПЛАНА) ГРАНИЦ ОБЪЕКТА ЗЕМЛЕУСТРОЙСТВА……..….23
2.2 КОНТРОЛЬ ЗА ПРОВЕДЕНИЕМ МЕЖЕВАНИЯ………………………..…28
2.3 НОРМАТИВНО – МЕТОДИЧЕСКАЯ БАЗА СОЗДАНИЯ ОПОРНО-МЕЖЕВОЙ СЕТИ……………………………………………………………….…29
2.4 КЛАССИФИКАЦИЯ ОПОРНОЙ МЕЖЕВОЙ СЕТИ И ЕЕ ТОЧНОСТЬ.…30

ГЛАВА 3. СОЗДАНИЕ ОПОРНЫХ МЕЖЕВЫХ СЕТЕЙ МЕТОДАМИ GPS…………………………………………………………………………………..36
3.1 ИСТОЧНИКИ ОШИБОК………………………………………………………38
3.2 РЕЖИМЫ РАБОТЫ GPS ПРИЕМНИКА…………………………………….40
3.3 ОБРАБОТКА РЕЗУЛЬТАТОВ GPS ИЗМЕРЕНИЙ………………………….43
3.4 ИМПОРТ ДАННЫХ ИЗ “GeoGenius” В GIS ”MapInfo”…………………….47

ГЛАВА 4.ПРОЛОЖЕНИЕ ОПОРНЫХ ТЕОДОЛИТНЫХ ХОДОВ, НАКОПЛЕНИЕ ДАННЫХ В ПОЛЯХ БАЗ ДАННЫХ, ОБРАБОТКА ИЗМЕРЕНИЙ В GIS MAPINFO…………………………………………...………50
4.1 РАБОТА С ТАБЛИЦАМИ MapInfo ………………………………………….50
4.2 РАБОТА С ТАХЕОМЕТРОМ-ПОЛУАВТОМАТОМ “Geodimeter 510”………………………………………………..………………………………..55
4.3 ПЕРЕДАЧА ДАННЫХ ИЗ ПАМЯТИ Geodimeter 510 В СРЕДУ GIS MapInfo………………………………………………..……………………………61
4.4 ОБРАБОТКА РЕЗУЛЬТАТОВ СЪЕМКИ В СРЕДЕ MapInfo. РАСЧЕТ ТЕОДОЛИТНЫХ ХОДОВ, ПРЯМОЙ ЗАСЕЧКИ, ОБРАТНОЙ ЗАСЕЧКИ, ЛИНЕЙНОЙ ЗАСЕЧКИ В СРЕДЕ MapInfo……………..………………………64
ГЛАВА 5.СОЗДАНИЕ ЦИФРОВЫХ КАРТ В СРЕДЕ MAPINFO…..…………81
5.1 СЛОИ КАРТЫ………………………………………………………..………..81
5.2 РАБОЧИЕ НАБОРЫ…………………………………………………….…….85
5.3 ВЫБОР ДЕТАЛЕЙ. ВЫБОР НА ЭКРАНЕ. ВЫБОР SQL-ЗАПРОСАМИ……………………………………………………………………..86
ТЕМАТИЧЕСКОЕ КАРТОГРАФИРОВАНИЕ…………………………….87
РЕДАКТИРОВАНИЕ И ИЗМЕНЕНИЕ СЛОЕВ КАРТЫ…………………92
ПОДПИСЫВАНИЕ КАРТЫ…………………………………………….…..93
РАСТРОВЫЕ ИЗОБРАЖЕНИЯ…………………………………………..…95
ВВОД ДАННЫХ - РАБОТА С ДИГИТАЙЗЕРОМ………………………...96
ГЕОГРАФИЧЕСКИЙ АНАЛИЗ……………………………………….…….98
СОЗДАНИЕ ОТЧЕТА……………………………………………………….101

ГЛАВА 6. ОХРАНА ТРУДА……………………………………………………..103
6.1. АНАЛИЗ ОПАСНЫХ И ВРЕДНЫХ ПРОИЗВОДСТВЕННЫХ ФАКТОРОВ ПРИ ВЫПОЛНЕНИИ ЗЕМЛЕУСТРОИТЕЛЬНЫХ РАБОТ НА ПЭВМ
6.2. РАЗРАБОТКА ТЕХНИЧЕСКИХ МЕРОПРИЯТИЙ ПО ОБЕСПЕЧЕНИЮ ЗДОРОВЫХ И БЕЗОПАСНЫХ УСЛОВИЙ ТРУДА ПРИ РАБОТЕ НА ПЭВМ ……………………………………………..106

ВЫВОДЫ И ПРЕДЛОЖЕНИЯ…………………………………………………..110
СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ……

Скачать полностью (334.25 Кб) Сколько стоит заказать работу?
Прикрепленные файлы: 1 файл

diplom v. 0.7.3.doc

— 923.00 Кб (Скачать документ)

s22=(Хр-Хв)2+(Ур-Ув)2                      (4.4.2.5)

Расхождения в координатах при первом и втором решениях должны удовлетворять условию:

r=3Mr, где                                   (4.4.2.6)


r=     (x’-x’’)2+(y’-y’’)2                       (4.4.2.7)

Mr – среднее квадратическое расхождение в положении пункта Р из двух решений

 Mr=     M12+M22,            , где M1 M2  – средние квадратические ошибки положения пункта Р из первого и второго решения.


Если расхождение r окажется допусимым, то за окончательное значение координат пункта Р берут среднее арифметическое, которое будет иметь ошибку М= Mr/2

При значительных расхождениях между М1 и М2 для достижения более точного окончательного результата можно воспользоваться формулами среднего весового

Хср = (X’P1 + X’’P2)/(P1+P2)                                     (4.4.2.8)

Yср = (Y’P1 +Y’’P2)/(P1+P2),  где P1 =1/M12      P2=1/M22              (4.4.2.9)

 

Средняя квадратическая ошибка окончательного положения пункта в этом случае найдется по формуле М =1/     Р1+Р2                      (4.4.2.10)


 

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

 

Declare Sub Main                                                              ‘Объявление

Declare Sub pr_zasechka                                                   ‘процедур

Include "Mapbasic.def"                                                     ‘

Global x(10),y(10),b(10),s(10),ro, M(10) As Float          ‘объявление

Global m_m,xp(10),yp(10), x1(5),y1(5),y2(3) As Float   ‘глобальных переменных

 

 

Sub Main                                                                                ‘тело главной процедуры

Set CoordSys NonEarth Units "m" Bounds(0,0)(1000000,1000000)

Dialog                                                                                      ‘открытие диалогового окна

Control StaticTExt                                                                  ‘

Title "x1 ="                                                                             ‘описание элементов

  Control EditText                                                                   ‘управления диалогом

  Into x(1)                                                                                ‘и занесение исходных

     Control StaticText                                                              ‘значений в соответствующие

     Title "y1 ="                                                                        ‘переменные

        Control EditText

        Into y(1)

     Control StaticText

     Title "b1 ="

        Control EditText

        Into b(1)

Control StaticText

Title "x2 ="

   Control EditText

   Into x(2)

      Control StaticText

      Title "y2 ="

         Control EditText

         Into y(2)

     Control StaticText

     Title "b2 ="

        Control EditText

        Into b(2)

Control StaticText

Title "x3 ="

   Control EditText

   Into x(3)

      Control StaticText

      Title "y3 ="

         Control EditText

         Into y(3)

     Control StaticText

     Title "b3 ="

        Control EditText

        Into b(3)

  Control StaticText

     Title "Ошибка измерения углов ="

        Control EditText

        Into m_m

Control OKButton                                                ‘описание элемента управления кнопка ОК

Calling pr_Zasechka                                             ‘и вызов процедуры pr_Zasechka  по нажатии                                              

Title "OK"                                                            ‘этой кнопки

Position 8,25

 

End Sub                                                                ‘конец тела главной процедуры

 

Sub pr_zasechka                                                     ‘тело процедуры pr_Zasechka

 

Dim  c_tgLeft, c_tgRight, c_tgY2,p(5),xfin,yfin, mfin As Float   ‘объявление локальных

Dim i As integer                                                                             ‘переменных

ro =206000                                                                                    

for i=1 to 2 Step 1                                                  ‘поскольку решать засечку надо дважды

y2(i)=180-b(i)-b(i+1)                                             ‘целесообразно это делать в цикле.

c_tgLeft = 1/tan(b(i)*DEG_2_RAD)                     ‘вычисление ctg левого,

c_tgRight =1/tan(b(i+1)*DEG_2_RAD)               ‘правого углов

c_tgY2 = 1/tan(y2(i)*DEG_2_RAD)                    ‘и угла    g

xp(i)=x(i)*c_tgRight+x(i+1)*c_tgLeft -y(i)+y(i+1)/(c_tgRight+c_tgLeft)     ‘расчет Xp(i), Yp(i)

yp(i)=(y(i)*c_tgRight+y(i+1)*c_tgLeft +x(i)-x(i+1))/(c_tgRight+c_tgLeft)  ‘

x1(i)=(x(i+1)*c_tgY2+xp(i)*c_tgRight-y(i+1)+yp(i))/(c_tgY2+c_tgRight)  ‘контрольн. расчет

y1(i)=(y(i+1)*C_tgY2+yp(i)*c_tgRight+x(i+1)-xp(i))/(c_tgY2+c_tgRignt)  ‘координат

create point(x1(i),y1(i))                        ‘исходного пункта  и его построение

s(i)=sqr((xp(i)-x(i))^2+(yp(i)-y(i))^2)                        ‘расчет расстояний

s(i+1)=sqr((xp(i)-x(i+1))^2+(yp(i)-y(i+1))^2)           ‘AP, BP, CP

m(i)=m_m*sqr(s(i)^2+s(i+1)^2)/ro*sin(Y2(i)) ‘расчет ср. квадр. ошибки для обоих решений

p(i)=1/(m(i)^2)                                                   ‘расчет P для каждого решения

next

xfin = (xp(1)*p(1)+xp(2)*p(2))/(p(1)+p(2))      ‘расчет окончательных координат

yfin = (yp(1)*p(1)+yp(2)*p(2))/(p(1)+p(2))      ‘определяемого пункта Р

Create Point (xfin, yfin)                                     ‘построение этого пункта

mfin=1/(sqr(p(1)+p(2)))                                     ‘расчет окончательного значения М и

Note(mfin)                                                          ‘вывод его значения на экран

End Sub

 

Блок-схему для решения этой задачи см. в Приложении 1

 

4.4.3 Расчет обратной засечки  в среде GIS MapInfo

    Обратная засечка(задача  Потенота) находит широкое применение  при создании съемочных сетей, привязке аэроснимков, выносе проектов  в натуру и в других случаях.

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

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

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

    Координаты пунктов A,B,C,D и измеренные углы мы знаем. Координаты пункта Р необходимо вычислить дважды, один раз с исползованием пунктов А,B,C и при следующем решении с использованием пунктов B,C,D

 

           Рис.4.4.2.  Схема расчета обратной засечки

Решение задачи производится по формулам Кнейсля:

Х’в=Хв-Ха    У’в=Ув-Уа                      (4.4.3.1)


Х’с=Хс-Ха    У’с=Ус-Уа

Х’d=Хd-Ха    У’d=Уd-Уа                                             (4.4.3.2)

k1 =ctgg У’в - Х’в              k2 =ctgg Х’в + У’в                        (4.4.3.3)


k3 =ctgg У’с - Х’с              k4 =ctgg Х’с + У’с

k5 =ctgg У’d - Х’d              k6 =ctgg Х’d + У’d

C1 =k2-k4/k1-k2 =ctg(AP)’                        (4.4.3.4)


C2 =k2-k6/k1-k5 =ctg(AP)’’

 

 

Y’1 =k2-C1k1/C12+1          Y’2 = k2-C2k1/C22+1                             (4.4.3.5)

X’1=C1Y’1                          X’2 = C2Y’2                                        (4.4.3.6)

X1 = Xa + X’1                       X2= Xa + X’2                                    (4.4.3.7)

Y1 = Ya + Y’1                       Y2 = Ya + Y’2                                   (4.4.3.8)

Формула окончательных координат пункта Р из двух решений:

Xp = X1+X2/2                       Yp = Y1+Y2/2                                    (4.4.3.9)

 

Declare Sub Main                                                  ‘объявление процедур

Declare Sub Obr_zasechka                                    ‘

Include "Mapbasic.def"                                          ‘ 

 

Global x(6),y(6),b(6),s(6),ro, M(6),m_m As Float     ‘объявление глобальных переменных

 

Sub Main                                                                ‘тело главной процедуры

Set CoordSys NonEarth Units "m" Bounds(0,0)(1000000,1000000)

 

Dialog                                                                      ‘открытие диалогового окна  

 

Control StaticTExt                                                   ‘описание элементов 

Title "x1 ="                                                              ‘управления диалогом

  Control EditText                                                    ‘и занесение исходных

  Into x(1)                                                                 ‘данных

     Control StaticText                                              ‘в соответствующие переменные

     Title "y1 ="

        Control EditText

        Into y(1)

     Control StaticText

     Title "b1 ="

        Control EditText

        Into b(1)

Control StaticText

Title "x2 ="

   Control EditText

   Into x(2)

      Control StaticText

      Title "y2 ="

         Control EditText

         Into y(2)

     Control StaticText

     Title "b2 ="

        Control EditText

        Into b(2)

Control StaticText

Title "x3 ="

   Control EditText

   Into x(3)

      Control StaticText

      Title "y3 ="

         Control EditText

         Into y(3)

     Control StaticText

     Title "b3 ="

        Control EditText

        Into b(3)

  Control StaticText

     Title "Ошибка измерения углов ="

        Control EditText

        Into m_m

Control OKButton                                                   ‘описание элемента управления

Calling obr_Zasechka                                              ‘кнопки ОК и вызов по нажатии

Title "OK"                                                               ‘кнопки процедуры obr_Zasechka

Position 8,25

 

End Sub                                                                   ‘конец тела главной процедуры

 

Sub obr_zasechka                                                    ‘тело процедуры obr_Zasechka

Dim  x1(10),y1(10),k(20),C(4),U(4),x2(3),y2(3), x3(3), xp,yp As Float ’объявление

Dim i As integer                       ‘локальных переменных

ro =206000

for i=1 to 2 Step 1                            ‘поскольку засечка решается дважды

                                                         ‘целесообразно для ее решения  использовать цикл

x1(i)=x(i+1)-x(i)                                                                 ‘Вычисление

x1(i+1)=x(i+2)-x(i)                                                             ‘X’в, Y’в

x1(i+2)=x(i+3)-x(i)                                                             ‘X’с, Y’с

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