Практическая работа по дисциплине «Обобщающий курс»

Практическая работа по дисциплине «Обобщающий курс»

Определение координат пункта IGS  в СК ITRF

Цель работы: освоение метода точного точечного ГНСС-позиционирования при определении координат пункта IGS в СК ITRF

Бюджет времени:6 часов аудиторной работы и 6 часов самостоятельной.

Исходные данные:

- дата и время наблюдения, постоянно действующая станция IGS (Приложение 1);

- свободное ПО RTKLib;

- презентация «Метод точного точечного позиционирования».

Остальные необходимые данные собираются студентами самостоятельно.

Содержаниеработы

1. Подготовка контрольных данных – получение координат XYZ с сайта ITRF в СКITRF2014 и на дату наблюдения

Ознакомиться с официальным сайтом ITRF http://itrf. ign. fr/.

Для получения координат пунктов ITRF на эпоху каталога 2014 г перейти по ссылке

http://itrf. ign. fr/ITRF_solutions/2014/doc/ITRF2014_GNSS. SSC. txt

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

Получение координат пунктов с сайта ITRFна дату наблюдения выполняется в порядке, приведенном в Приложении 2.

Результаты получения координат и скоростей пунктов IGS занести в таблицу.

Координаты и скорости пунктов IGS

Пункт	Координаты	X, м	mX, м	Y, м	mY, м	Z, м	mZ, м
…	ITRF 2014
Дата варианта
Скорости, м/год	VX	mX’	VY	mY’	VZ	mZ’
…	ITRF 2014
Дата варианта

Совершить переход от пространственных прямоугольных координат пунктов X, Y,Zк геодезическим эллипсоидальным LBH в геоцентрической системе координат с параметрами эллипсоида WGS-84, по формулам:

,

,

,        

rp=(X2+Y2)Ѕ, a= 6378237 м, e2 = 0,00669437999013.

Примечание: вместо вычислений по формулам можно грамотно воспользоваться геодезическим калькулятором или другим ПО (см. Интернет)

Результаты вычислений оформить в виде таблицы

Эллипсоидальные координаты пунктов  IGS в СК ITRF

Пункт	Дата	B, 0‘“	L, 0‘“	H, м
1	2014
Дата варианта
2	2014
Дата варианта

2. Обработка GPS-измерений в ПО Rtklib

программа – в архиве. После разархивации папки зайти в папку bin и запустить программу Rtklaunch.

При обработке пользоваться указаниями в презентации PPP. Презентация прилагается отдельным файлом.

2.1  Подготовка исходных данных с помощью утилиты RTKGET.

Заполнить таблицу для всех скачанных файлов

№ пп	Вид информации	Название файла	Размер файла, КБ
…

…

Примечание. Файл с приливными нагрузками можно не брать

2.2. Вычисление координат методом точного точечного позиционирования в ПО RTKLib. Постобработка RTKPOST.

2.2.1 Обработка в режиме «Статика».

2.2.2 Обработка в режиме «Кинематика».

3. Обработка GPS и ГЛОНАСС измерений.

3.1 Подготовка необходимых исходных данных с Российского сайта Системы высокоточного определения эфемеридно-временных поправок (СВО ЭВП).

http://www. glonass-svoevp. ru/index. php? option=com_content&view=article&id=55&Itemid=259&lang=ru

3.2 Обработка в режимах  статика и кинематика в RTKPOST, GPS+ГЛОНАСС, аналогично п.2.

4.  Обработка с помощью онлайн-службы GAPS.

Выполнить обработку данных с помощью онлайн-службы (см. презентацию).

5. Оценка точности и анализ результатов вычисления.

5.1Привести в отчет графики:

GPS статика (ITRF),

GPSкинематика (ITRF),

GPS+ГЛОНАСС статика (ПЗ-90.11),

GPS+ГЛОНАСС кинематика (ПЗ-90.11),

По каждому графику указать время сходимости, отклонение от среднего на участке после того, как графики “выровнялись”, - на восток Е, север N и вверх U.

5.2 Работа с файлами *.pos

По каждому методу обработки сохранить файл с расширением *.pos.

И файла привести в отчет по работе заголовок, несколько первых и несколько последних записей.

Из самой последней записи выписать координаты пункта X, Y,Z и их ошибки mx, my, mz.

Заполнить таблицу

Вид решения	X, м	mx, м	Y, м	my, м	Z, м	mz, м
ITRF на дату
GPS статика (ITRF)
GPSкинематика (ITRF
GPS+ГЛОНАСС статика (ПЗ-90.11)
GPS+ГЛОНАСС кинематика
Онлайн-служба GAPS

5.3 Сравнение результатов вычисления координат

Сравнить полученные координаты с контрольными данными (вычислить разности вычисленных координат с контрольными значениями ITRF). Заполнить таблицу

Разности с контрольными значениями ITRF на дату	ΔX, мм	ΔY, мм	ΔZ, мм
GPS статика (ITRF)
GPSкинематика (ITRF
GPS+ГЛОНАСС статика (ПЗ-90.11)
GPS+ГЛОНАСС кинематика
Онлайн-служба GAPS

Сделать вывод, попадает ли разность вычисленных и контрольных координат в доверительный интервал, задаваемый СКП.

6. Перевычисление координат.

6.1 Перевычислить координаты с даты наблюдения t на эпоху ITRF2014:

X2014 = X(дата) + VX(t-2014)

Y2014 = Y(дата) + VY(t-2014)

Z2014 = Z(дата) + VZ(t-2014)

t = год + (количество суток от начала года)/365.

Заполнить таблицу

Вид решения	X, м	mx, м	Y, м	my, м	Z, м	mz, м
ITRF2014
GPS статика (ITRF)
GPSкинематика (ITRF
GPS+ГЛОНАСС статика (ПЗ-90.11)
GPS+ГЛОНАСС кинематика
Онлайн-служба GAPS

Заполнить таблицу

Разности с контрольными значениями ITRF2014	ΔX, мм	ΔY, мм	ΔZ, мм
GPS статика (ITRF)
GPSкинематика (ITRF
GPS+ГЛОНАСС статика (ПЗ-90.11)
GPS+ГЛОНАСС кинематика
Онлайн-служба GAPS

6.2 Привести формулы  перевычисления координат XYZ из ITRF в ПЗ-90.11

2017 г

Приложение 1 Список пунктов IGS Номер варианта – последние 2 цифры в зачетке

№ вар	IGS	дата	№ вар	IGS	дата	№ вар	IGS	дата
BADG	15.01	11	BADG	15.02	21	BADG	15.03
IRKJ	12	IRKJ	22	IRKJ
KZN2	13	KZN2	23	KZN2
MDVJ	14	MDVJ	24	MDVJ
MOBJ	15	MOBJ	25	MOBJ
MOBK	16	MOBK	26	MOBK
NVSK	17	NVSK	27	NVSK
SVTL	18	SVTL	28	SVTL
TIXI	19	TIXI	29	TIXI
ZECK	20	ZECK	30	ZECK

№ вар	IGS	дата	№ вар	IGS	дата	№ вар	IGS	дата
31	BADG	15.04	41	BADG	15.05	51	BADG	15.12
32	IRKJ	42	IRKJ	52	IRKJ
33	KZN2	13	KZN2	53	KZN2
34	MDVJ	44	MDVJ	54	MDVJ
35	MOBJ	45	MOBJ	55	MOBJ
36	MOBK	46	MOBK	56	MOBK
37	NVSK	47	NVSK	57	NVSK
38	SVTL	48	SVTL	58	SVTL
39	TIXI	49	TIXI	59	TIXI
40	ZECK	50	ZECK	60	ZECK

№ вар	IGS	дата	№ вар	IGS	дата	№ вар	IGS	дата
61	BADG	15.10	71	BADG	15.11	81	BADG	15.09
62	IRKJ	72	IRKJ	82	IRKJ
63	KZN2	73	KZN2	83	KZN2
64	MDVJ	74	MDVJ	84	MDVJ
65	MOBJ	75	MOBJ	85	MOBJ
66	MOBK	76	MOBK	86	MOBK
67	NVSK	77	NVSK	87	NVSK
68	SVTL	78	SVTL	88	SVTL
69	TIXI	79	TIXI	89	TIXI
70	ZECK	80	ZECK	90	ZECK

Приложение 2

Получение координат пунктов IGScсайт ITRF на дату наблюдения

Координаты станций IGS на дату наблюдения с сайта

http://itrf. ign. fr/

Определить номер домена или 4-х символьный идентификатор из списка станций:IERSNetwork – Networkdescription – получить текстовый файл с информацией о станцияхiers_sta_list. txt. Пример: станция arti – идентификатор – ARTU, номер домена – 12362M001

GetITRFcoord(получить координаты, слева на странице);selectedpoints(выбор точек); addpointstothecart(добавить точку в карту);  выбрать точку по номеру домена или по идентификатору;addtocart (добавить в карту);displayselectedpoints(вывод на экран информации о точке);Get ITRF coordinatesofthesepoints(получение координат этих точек); указать эпоху в требуемом формате (год, месяц, день);  согласиться ОК;  получить и распечатать/скопировать таблицу координат и скоростей точки на дату наблюдения.

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