Главная страница ДонНТУ      Портал магистров ДонНТУ     Кафедра геоинформатики и геодезии    


в библиотеку

Автор статьи: А. М. Ицков, геодезист.  г. Иркутск

Источник: http://www.dataplus.ru/Support/ProjUtility/Datum.htm

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

Как известно, при вычислении координат съемочного обоснования,  в измеренные величины не вводятся поправки за редукцию с поверхности Земли на поверхность референц-эллипсоида и за редуцирование с эллипсоида на плоскость (далее – редукции и поправки за редукции). До появления компьютеров, ввод поправок за редукции занимал бы слишком много времени. Поэтому действующие инструкции (порядком устаревшие), не требуют ввода поправок за редукцию. В настоящее время расчеты, как правило, производятся при помощи специальных программ и учет редукций не усложняет камеральные работы. Однако считается, что в вводе поправок за редукции нет необходимости вследствие того, что они малы. Так ли это?
Рассмотрим этот вопрос.
В данной статье мы рассмотрим системы координат, параметры которых известны.

В России, в качестве государственной, принята система координат 95 года (СК-95г). В ней принят эллипсоид Красовского и проекция Гаусса-Крюгера. Данная проекция является частным случаем  Поперечно-цилиндрической проекции Меркартора. Для обеспечения единой геодезической основы больших и протяженных объектов, например нефтепровод «Восточная Сибирь – Тихий океан», иногда применяется система координат 95 года.

В настоящее время все еще широко используется система координат 63 года (СК-63г). В этой системе координат так же принят эллипсоид Красовского и проекция Гаусса-Крюгера. Ее отличие состоит в том, что осевые меридианы , смещения на север (No) и восток (Eo) сдвинуты относительно СК-95г. При чем, для каждой зоны смещения разные. Кроме того, в СК-95г зоны имеют ширину в 6 градусов, а в СК-63г ширина зоны 3 градуса. Вследствие этого параметры каждой отдельной зоны СК-63г неизвестны. Однако нам известны общие свойства этой системы координат.

Так же, для геодезической основы изысканий и строительства объектов, начали использовать систему координат WGS-84, UTM (Проект Сахалин-2, Береговые трубопроводы). В системе WGS-84, UTM принят эллипсоид WGS-84 и проекция Гаусса-Боага, так же являющаяся частным случаем поперечно-цилиндрической проекции Меркартора.

Параметры систем координат приведены в таблице 1.

Таблица1

Система координат СК-95г WGS-84, UTM
Проекция (Гаусса-Крюгера) (Гаусса-Боага)
Смещение на север No 0 0
Смещение на восток Eo 500000 500000
Масштаб по осевому мередиану ( mo ) 1 0,9996
Применяемый эллипсоид Красовского WGS-84
Большая полуось эллипсоида, м 6378245 6378137
Сжатие эллипсоида 1/298.3 1/298.257223563

Как мы видим, различие систем координат СК-95г и WGS-84, UTM состоит в несколько разных размерах большой полуоси исходных эллипсоидов и в довольно существенной разнице масштабов по осевому меридиану. Поскольку системы координат используют одну проекцию, то для съемочных сетей, формулы вычисления поправок за редукцию измеренных величин будут общими. Мы не будем приводить всех поправок, которые крайне малы и вводятся только при обработке высокоточных геодезических сетей. Рассмотрим только редукцию измеренных линий на эллипсоид и на плоскость.

Формулы вычисления редукций.

Редуцирование измеренных линий на поверхность референц- эллипсоида.

Можно записать:

Dэ=D*(1-H/R) (2)

где:   D    -   горизонтальное проложение измеренной линии на местности;

Dэ   -   значение линии, приведенной на референц-эллипсоид;

H     -   средняя высота линии над референц-эллипсоидом;

R     -   средний радиус кривизны эллипсоида в средней точке линии;

Для территории России принимаем среднее значение  R=6380000 м  [1]

Редуцирование линий на плоскость

Поскольку стороны в съемочном обосновании редко превышают 1000 м, то формулу 1 [1] можно записать в следующем виде:

Dп=Dэ*mo*(1+YmІ/2RІ)   (2)

где   Dп   -  значение линии редуцированной на плоскость;

Ym   -  среднее значение ординаты (Расстояние от осевого меридиана) линии в проекции  (Y1+Y2)/2;

mo   -  значение масштабного коэффициента по осевому меридиану;

На основании формул 1 и 2 получаем формулу суммарной поправки в линии:

D=D*(1-H/R)*mo*(1+YmІ/2RІ) (3)

В СК-95г и WGS-84, UTM для исключения отрицательных значений, ордината имеет смещение на +500000м. Формула (3) будет иметь вид:

D=D*(1-H/R)*mo*(1+(Y-500000)І/2RІ) (4)

где: Y   -  значение ординаты из каталога;

Используя формулу ( 4 ) вычислим поправки за редукцию на эллипсоид и за редукцию на плоскость в системах координат СК-95г и WGS-84, UTM в линии длиной 1000 метров.

 

Таблица 2. Поправки в линию длиной 1000м в СК-95г и СК-63г.

Отметка, м Расстояние от осевого меридиана, км
0 50 100 125 150 200 250
0 0,000 0,031 0,123 0,192 0,276 0,491 0,768
250 -0,039 -0,008 0,084 0,153 0,237 0,452 0,729
500 -0,078 -0,048 0,044 0,114 0,198 0,413 0,689
750 -0,118 -0,087 0,005 0,074 0,159 0,374 0,650
1000 -0,157 -0,126 -0,034 0,035 0,120 0,335 0,611

 

Таблица 3. Поправки в линию длиной 1000м в СК WGS-84, UTM.

Отметка, м Расстояние от осевого меридиана, км
0 50 100 150 200 250
0 -0,400 -0,369 -0,277 -0,124 0,091 0,367
250 -0,439 -0,408 -0,316 -0,163 0,052 0,328
500 -0,478 -0,448 -0,356 -0,202 0,013 0,289
750 -0,517 -0,487 -0,395 -0,241 -0,026 0,250
1000 -0,557 -0,526 -0,434 -0,280 -0,066 0,211

 

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

Возникает вопрос: «Всегда ли надо вводить поправки?»

Среднеквадратическая погрешность линейных измерений при создании планового обоснования должна быть не менее 1/4000. Это среднеквадратическая ошибка 0.250м на один километр. Пренебрегаемо малыми можно считать ошибки, не превышающие 10 – 20% от среднеквадратической. Возьмем величину в 20% - получим 0.050м на километр.

В таблицах 2 и 3 желтым цветом выделены поправки менее 0.050 м на один километр. То есть, когда поправки можно не учитывать. Жирным курсивом выделены ошибки более среднеквадратической (более 0.25м на 1000м), когда поправки однозначно необходимо вводить в измеренные линии.

Система координат 95года.  Ширина зоны в СК-95г равна 6 градусам. На территории России максимальное расстояние от осевого меридиана до края зоны составляет порядка 250 км. На краю зоны поправки за редуцирование могут достигать 0.76м на 1000м. Однако стоит отметить, что в диапазоне: высот от 0 до 500м и на расстоянии от осевого меридиана от 0 до 50км, редукции весьма малы. На каждом отдельном объекте геодезистам необходимо оценивать поправки за редукции и самостоятельно принимать решение о целесообразности их введения в измеренные линии.

Система координат 63г.  Ширина зоны в СК-63г равна 3 градусам. Максимальное расстояние до края зоны в России составляет порядка 125 км. Обратимся к таблице 2. Учитывая, что в подавляющем большинстве случаев расстояние от осевого меридиана редко достигает 100км, можно сделать вывод о том, что в большинстве случаев поправки за редукцию измеренных линий будут небольшими. В том случае, если все необходимые параметры известны, то можно рассчитать и оценить величину поправки за редукцию.

Система координат  WGS-84, UTM.  Поправки значительные. Небольшие по своей величине поправки только при расстоянии 200км от осевого меридиана. При работе с системой координат WGS-84, UTM поправки за редуцирование линий необходимо вводить постоянно.

Уравнивание геодезических сетей по способу наименьших квадратов (МНК) производится при условии, что:

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

Тем не менее, считается, что поправки при уравнивании компенсируют поправки за редукции. Поэтому при создании планового обоснования нет необходимости в учете этих редукций.  Рассмотрим этот вопрос на примере двух теодолитных ходов.

Уравнивание теодолитных ходов.

Воспользуемся программой CREDO DAT. Создадим проект теодолитного хода с нулевой невязкой, с учетом параметров СК-42г. Координаты возьмем в проекции Гаусса-Крюгера, ордината – 150 км, отметка 50м над уровнем моря. Суммарная поправка 0.269м на 1000 метров. Получим таким образом «истинные» координаты. После этого установим систему координат «Местная» и отключим учет поправок. Смоделируем уравнивание хода без ошибок в измерениях, но и без редукций. Сравним полученные результаты.

Пример первый.

Ход в виде буквы П.  Схема хода на рис.1.  Журнал измерений на рис.2.

 

Таблица 4. Сравнение координат точек хода из двух вариантов уравнивания.

Точки Координаты точек ходов, в м Ошибки в определении координат, м
С учетом поправок за редукцию Без учета поправок
Хи Х Y H Xи-X Yи-Y S
4 6000600,169 650000,016 50,000 6000600,008 650000,043 50,000 0,161 -0,027 0,163
5 6000900,250 650000,023 50,000 6000900,008 650000,063 50,000 0,242 -0,040 0,245
6 6000900,242 650300,104 50,000 6000900,000 650300,066 50,000 0,242 0,038 0,245
7 6000600,161 650300,096 50,000 6000600,000 650300,071 50,000 0,161 0,025 0,163

 

Таблица 5. Характеристики теодолитного хода.

Длина N Fb факт. Fb доп. Fx Fy Fs [S]/Fs
2100 8  0°00'00" 0°02'50" 0,000 -0,080 0,080 25926

 

 Рис. 1  Рис. 2
Схема хода Журнал измерений

 

Линейная невязка, в характеристиках хода (табл.5), составляет всего 8 см. Однако погрешности по линии 2-6 и по линии 7-1 компенсируют друг друга и на конечную плановую невязку хода не влияют. Полные погрешности вычисления координат точек в слабом месте хода ( точки 5 и 6 ) составляют 0,245м.  Много это или мало? Длина хода 2100м. Согласно СП 11-104-97 ход такой протяженности может служить плановым обоснованием для топографической съемки масштаба 1:1000 ( табл. 5.1, прим.1 [3] ).

Допустимая средняя ошибка определения координат точек хода составляет 0.1 мм в масштабе плана. В нашем случае – 0.10  м. Соответственно средняя квадратическая ошибка – 0.14 м и предельная ошибка – 0.28 м.

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

 

Пример второй.

Ход с координатной привязкой. Схема хода на рис. 3.  Журнал измерений на рис. 4.

 

Таблица 6. Сравнение координат точек из двух вариантов уравнивания.


Точки Координаты точек ходов, в м. Ошибки в определении координат, м
С учетом поправок за редукцию Без учета поправок
X Y H Xи-X Yи-Y S
5 6000000,008 650000,004 50,000 5999999,949 649999,837 50,000 0,059 0,167 0,177
6 6000300,089 650000,004 50,000 6000300,010 649999,816 50,000 0,079 0,188 0,204
7 6000600,169 650000,003 50,000 6000600,070 649999,825 50,000 0,099 0,178 0,204
8 6000900,250 650000,003 50,000 6000900,131 649999,871 50,000 0,119 0,132 0,178
9 6001200,330 650000,003 50,000 6001200,192 649999,962 50,000 0,138 0,041 0,144

 

Таблица 7. Сравнение дирекционных углов сторон хода

Сторона Дирекционный угол стороны хода Аи-Ак
С учетом редукций Без учета редукций
9-10 0°00'00" 0°01'38" 01'38"
10-11 300°00'00" 300°02'15 02'15"
11-12 240°00'00" 240°02'34" 02'34"

 

Таблица 8. Характеристики теодолитного хода.

Длина N Fb факт. Fb доп. Fx Fy Fs [S]/Fs
2250 11 -0°03'19 0°03'00" -0,52 0,26 0,58 3876

 

рис. 3. Схема хода Рис. 4.  Журнал измерений хода
Схема хода Журнал измерений

В данном примере полные погрешности вычисления координат точек в слабом месте хода (точки 6 и 7) составляют 0.204 м. Длина хода 2250 м. Согласно СП 11-104-97, ход такой протяженности может служить плановым обоснованием для топографической съемки масштаба 1:1000 (табл.5.1, прим.1 [3]). Допустимая ошибка определения координат точек хода составляет 0.1 мм в масштабе плана. В нашем случае – 0.1 м. Соответственно средняя квадратическая ошибка – 0.14 м и предельная ошибка – 0.28 м.

Полученная погрешность составляет 70% от предельной погрешности.

Вследствии того, что ход с координатной привязкой, искажаются не только линии, но и азимуты сторон. В таблице 7 приведены дирекционные углы сторон хода. Искажения достигают 02'34". Так же показателен тот факт, что вычисленная фактическая угловая ошибка хода ( табл. 8 ) превышает допустимую погрешность.

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

 

Рекомендации по обработке геодезических измерений.

1. При помощи программы CREDODAT.

Создание различных систем координат в проекте CREDO DAT 3.0 детально описано в главе 3 описаний к программе и в статье «Создание пользовательской системы координат в CREDO_DAT»  [ 1 ]. Параметры проекций Гаусса-Крюгера и UTM приведены в таблице 1  данной статьи. Однако необходимо добавить, что для правильного редуцирования измеренных линий на эллипсоид, необходимо вводить отметки точек. При совместной обработке нивелирования и плановых сетей проблем не возникает. Если же уравнивается только плановое обоснование, то можно вводить отметки для каждой точки. Либо можно ввести среднюю отметку на объекте.  При правильной настройке проекта все поправки вносятся в измеренные величины автоматически.

2. При помощи программы «TRIMBLE GEOMATIC OFFICE» (далее TGO).

В TGO необходимо установить систему координат. Можно воспользоваться готовой или создать свою. Так же, после установления системы координат, поправки определяются автоматически. В программе TGO при решении обратной задачи рассчитывается длина линии в проекции (Crid Distance), на эллипсоиде (Ellipsoid Distance), на Земной поверхности (Ground Distance). См. рис 5. Это весьма удобно. Можно сразу оценить разницу между расстоянием измеренным на местности (Ground Distance) и расстоянием в проекции (Crid Distance). В примере на рис.5, расстояние на местности равно 1070.241м , а расстояние в проекции равно 1069.972м. Разница составляет 0.269м на 1070 метров.

 

Рис. 5

 

3. Ввод поправок в измерения за редукции «вручную» при помощи таблиц.

Можно составить таблицы поправок для конкретного объекта, рассчитав поправки по формуле ( 3 ). Выбрать интервал значений Y (от 50км до 65 например) и среднюю отметку для объекта. По одной оси откладывать Y, по другой длины линий.  Интервал линий рекомендуется брать через 10метров (210м, 220м, 230м и т.д.). Интервал значений Y через 10 км.

 

Выводы и рекомендации.

  1. При обработке планового съемочного обоснования, необходимо вводить поправки за редукцию измеренных величин с поверхности Земли на поверхность референц-эллипсоида и за редуцирование измеренных величин с эллипсоида на плоскость.
  2. Для обработки геодезических сетей рекомендуется выбирать программное обеспечение, в котором поправки вычисляются и вводятся в измеренные величины автоматически, не усложняя процесс обработки геодезических сетей.
  3. При работе с электронными тахеометрами в режиме вычисления координат, а так же при выносе в натуру, рекомендуется учитывать параметры проекции и отметку. Либо вводить усредненный масштабных коэффициент для района работ.

 

Список использованной литературы:

  1. Создание пользовательской системы координат в CREDO_DAT. Журнал «Автоматизированные технологии изысканий и проектирования» №13/2004, стр 58.
  2. Яковлев Н.В. Высшая геодезия. – Москва.: Недра, 1989.
  3. Инженерно-геодезические изыскания для строительства» СП 11-104-97, дата введения 01.01.1998г
Магистерская работа Cсылки Отчет о поиске Автобиография Индивидуальное задание