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

y1(i)=y(i+1)-y(i)                                                                 ‘X’d Y’d

y1(i+1)=y(i+2)-y(i)                                                             ‘по формулам

y1(i+2)=y(i+3)-y(i)                                                             ‘ 

k(i)=(1/tan(b(i)*DEG_2_RAD))*y1(i)-x1(i)                          ‘вычисление к1,к2,к3,к4,к5,к6

k(i+1)=(1/tan(b(i)*DEG_2_RAD))*x1(i)+y1(i)                     ‘по формулам

k(i+2)=(1/tan(b(i+1)*DEG_2_RAD))*y1(i+1)-x1(i+1)         ‘ 

k(i+3)=(1/tan(b(i+1)*DEG_2_RAD))*x1(i+1)+y1(i+1)        ‘  

k(i+4)=(1/tan(b(i+2)*DEG_2_RAD))*y1(i+2)-x1(i+2)         ‘ 

k(i+5)=(1/tan(b(i+2)*DEG_2_RAD))*x1(i+2)+y1(i+2)        ‘

c(i)=(k(i+1)-k(i+3))/(k(i)-k(i+2))                                       ‘вычисление с1,с2

c(i+1)=(k(i+1)-k(i+4))/k(i)-k(i+4)                                     ‘по формулам

U(i)=(k(i+1)-c(i)*k(i))/(c(i)^2+1)                                               ‘вычислеине У’1, У‘2

U(i+1)=(k(i+1)-c(i+1)*k(i))/(c(i+1)^2+1)                                  ‘по формулам

x2(i)=c(i)*u(i)                                                     ‘вычисление Х’1,  Х‘2

x2(i+1)=c(i+1)*u(i+1)                                        ‘

x3(i)=x(i)+x2(i)                                           ‘вычисление Х1, Х2

x3(i+1)=x(i)+x2(i+1)                                  ‘

y2(i)=y(i)+U(i)                                   ‘вычисление У1, У2

y2(i+1)=y(i)+U(i+1)                          ‘

xp=(x3(i)+x3(i+1))/2         ‘вычисление окончательных

yp=(y2(i)+y2(i+1))/2         ‘координат пункта Р

Create Point (xp,yp)

next

 

End Sub

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

 

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

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

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

Решение задачи, осуществим дважды, с использованием исходных пунктов А, В, С и  В, С, D.

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

 

(АВ) = arc tg (Yв-Ya)/(Xв-Xa)                        (4.4.4.1)

Sав =      (Хв-Ха)2 +(Ув-Уа)2                            (4.4.4.2)

(ВС) = arc tg (Yc-Yв)/(Xc-Xв)                         (4.4.4.3)

Sвc =      (Хc-Хв)2 +(Уc-Ув)2                            (4.4.4.4)

 

На основании теоремы косинусов определяем углы в1 и в2

 

в1 =arc cos(S2aв+S12-S22)/2SавS1                       (4.4.4.5)

в2 =arc cos(S2вс+S22-S32)/2SвсS2                       (4.4.4.6)

 

После этого, вычисляем дирекционные углы линий АР и ВР

(АР) = (АВ)+в1                      (4.4.4.7)

(ВР) = (ВС)+в2                       (4.4.4.8)

 

После нахождения дирекционных углов (АР) и (ВР) мы можем определить координаты пункта Р из двух решений

X’ = Xa + S1cos(AP)                       (4.4.4.9)

X’’ = Xв + S2cos(ВР)                      (4.4.4.10)

Y’ = Ya + S1sin(AP)                        (4.4.4.11)

Y’’ = Yв + S2sin(BP)                       (4.4.4.12)

Вычисляем углы засечки g1 g2

g1 =arc sin Sавsinв1/S2                       (4.4.4.13)

g2 =arc sin Sвcsinв2/S3                       (4.4.4.14)

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

М1 = msÖ2/sing1                       (4.4.4.15)

М12= msÖ2/sing2                       (4.4.4.16)

Вычисляем Mr = ÖM12 + M22,  r=Ö(X’-X’’)2 + (Y’-Y’’)2 и проверяем выполнение условия r <= 3Mr                       (4.4.4.17)

Вычисляем окончательные значения координат пункта Р и среднюю квадратическую ошибку его положения

Xp = (X’ +X’’)/2        Yp = (Y’ + Y’’)/2                   M=Mr/2       (4.4.4.18-19-20)

 

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

Declare Sub lin_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 "s1 ="

        Control EditText

        Into s(1)

 

 

Control StaticText

Title "x2 ="

   Control EditText

   Into x(2)

      Control StaticText

      Title "y2 ="

         Control EditText

         Into y(2)

     Control StaticText

     Title "s2 ="

        Control EditText

        Into s(2)

 

Control StaticText

Title "x3 ="

   Control EditText

   Into x(3)

      Control StaticText

      Title "y3 ="

         Control EditText

         Into y(3)

     Control StaticText

     Title "s3 ="

        Control EditText

        Into s(3)

 

  Control StaticText

     Title "Ошибка измерения сторон ="

        Control EditText

        Into m_m

 

Control OKButton                             ‘объявление элемента управления диалогом

Calling lin_Zasechka                         ‘кнопка ОК и вызов по нажатии этой кнопки

Title "OK"                                         ‘процедуры lin_Zasechka

Position 8,25                                      ‘    

 

 

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

 

 

Sub lin_zasechka                               ‘ тело процедуры lin_Zasechka

 

Dim l(10),l1(10),s1(10),chisl,chisl2,x1(10),y1(10),U(3),mr,xfin,yfin,mfin As Float    ‘

Dim i As integer                            ‘объявление локальных переменных

ro =206000

for i=1 to 2 Step 1

 

l(i)=atn((y(i+1)-y(i))/(x(i+1)-x(i)))*RAD_2_DEG       ‘вычисление дирекционных углов

s1(i)=sqr((x(i+1)-x(i))^2+(y(i+1)-y(i))^2)                    ‘и длин исходных линий

chisl=s1(i)^2+s(i+1)^2-s(i+1)^2                                   ‘ 

b(i)=acos(chisl/(2*s1(i)*s(i)))                                       ‘вычисление углов в1 и в2

l1(i)=(l(i)+b(1))*DEG_2_RAD                     ‘ вычисление дирекц.углов линий АР и ВР

x1(i)=x(i)+s(i)*cos(l1(i))                               ‘определение Х‘ и У‘

y1(i)=y(i)+s(i)*sin(l1(i))                                ‘для обоих решений засечки

 

chisl2=s1(i)*sin(b(i))                                                    ‘вычисление углов засечки g и

u(i)=asin(chisl2/s(i+1))                                                 ‘ср.квадратических ошибок

m(i)=(m_m*(sqr(2)))/sin(u(i))                                      ‘ для обоих решений

 

next

mr=sqr(m(1)^2+m(2)^2)                                              ‘вычислеине Mr

xfin=(x1(1)+x1(2))/2                                                    ‘вычисление окончательных

yfin=(y1(1)+y1(2))/2                                                    ‘координат пункта Р

mfin=mr/2                                                                     ‘вычислеине ср.квадратич.ошибки

Create Point(xfin,yfin)                                                  ‘для окончательного положения

Note(mfin)                                                                    ‘пункта Р и вывод его значения на

End Sub                                                                        ‘экран

 

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

 

 

 

 

 

 

 

 

 

 

ГЛАВА 5. СОЗДАНИЕ ЦИФРОВЫХ КАРТ В СРЕДЕ GIS MapInfo

Человечеству доступны огромные объемы информации, которая хранится на бумаге и в памяти компьютеров. По разным оценкам, до 80% всех баз данных содержат какую-либо географическую информацию. MapInfo – это мощное средство анализа данных. Вы можете придать графический вид статистическим и прочим данным. К данным применяются графические операторы, а доступ к ним можно оформлять через запросы. MapInfo совмещает преимущества обработки данных, которыми обладают системы управления базами данных (СУБД), включая мощный язык запросов SQL и наглядность карт, схем и графиков.

Возможности MapInfo.

• просмотр данных в удобном виде;
• доступ к файлам, созданным другими программами;
• доступ к отдаленным базам данных;
• многослойность карт;
• анализ данных с помощью тематических карт;
• составление запросов различной сложности;
• сохранение Рабочих наборов;
• OLE-вложение окон MapInfo в документы других программ;
• мощный набор функций изменения вида других программ (рисование, редактирование);
• наборы готовых карт и функции создания своих карт;
• создание печатных отчетов;
• изменение проекций карт;
• справочная система.

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

Все, что вы будете делать в MapInfo, начинается с открытия одной или нескольких таблиц.

Таблица состоит, по крайней мере, из двух различных файлов; первый содержит данные, а второй - описание структуры данных:

• <имя файла>.TAB. Этот файл содержит описание структуры данных таблицы. Он представляет собой текстовый файл, описывающий формат того файла, который содержит данные.
• <имя файла>.DAT или <имя файла>.WKS,.XLS. Этот файл содержит табличные данные. Если Вы работаете с файлами dBase/FoxBASE, ASCII с разделителями, Lotus 1-2-3 или Microsoft Excel, таблица MapInfo будет состоять из файла с расширением .TAB и либо файла данных, либо файла электронной таблицы. Таблицы, содержащие растровые изображения, хранят данные в файлах-компонентах форматов BMP, TIF или GIF.

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

• <имя файла>.MAP. Этот файл описывает графические объекты.
• <имя файла>.ID. Этот файл содержит список указателей (индексов) на графические объекты, позволяющий MapInfo быстро находить объекты на карте.

Таблица может содержать еще и индексный файл. Индексный файл позволяет проводить поиск объектов на карте с помощью команды Find. Если Вы хотите искать дома, улицы, города, области и т.д. используя команду Find, соответствующие поля таблицы должны быть проиндексированы. Индекс хранится в файле: <имя файла>.IND.

Для того чтобы запустить MapInfo, необходимо нажать на кнопку «Пуск» и в меню «Программы» выбрать MapInfo/MapInfo Professional или запустить ее с рабочего стола (при наличии файла запуска программы на нем). Через мгновение появится диалог начала сеанса работы с MapInfo.

В нем Вы можете указать, с чего следует начать. Если Вы уже работали ранее с MapInfo, Вы можете восстановить то состояние, которое было на экране в конце последнего сеанса работы, выбрав «Restore Previous Session», или открыть последний использовавшийся Рабочий набор. Для того чтобы открыть таблицу, выберите «Table» и нажмите кнопку ОК. Появится диалог «Open Table».

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