Работа с программой Teqc

Ну, «раз пошла такая пьянка», давайте продолжим тему и попробуем общими усилиями решить следующую интересную задачку. На примере того же файла svtl0470.11o.

Выполним команду
teqc.exe -O.s G -R +ds 86400 -O.sum s svtl0470.11o > svtl0470_8.11o

где -O.s G - задаёт систему GPS, так как удаляем опцией –R все спутники ГЛОНАСС.
Хотя опция -O.s G и не влияет на вычисления, но для строгости оформления файла, так правильнее.

Получился файл только со спутниками GPS.
Анализируем, и видно, что в интервалах времени:

11 2 16 7 31 30.0000000 0 13G03G06G08G11G14G15G16G18G19G21G22G27
G28
11 2 16 7 32 30.0000000 0 13G03G06G08G11G14G15G16G18G19G21G22G27
G28


11 2 16 8 25 0.0000000 0 13G03G06G09G11G14G15G18G19G22G24G27G28
G32
11 2 16 8 36 0.0000000 0 13G03G06G09G11G14G15G18G19G22G24G27G28
G32


11 2 16 9 25 0.0000000 0 13G03G06G09G11G14G18G19G20G22G24G27G28
G32
11 2 16 9 27 0.0000000 0 13G03G06G09G11G14G18G19G20G22G24G27G28
G32


11 2 16 9 28 30.0000000 0 13G03G09G11G14G17G18G19G20G22G24G27G28
G32
11 2 16 9 41 30.0000000 0 13G03G09G11G14G17G18G19G20G22G24G27G28
G32


11 2 16 13 18 0.0000000 0 13G01G02G04G10G13G17G20G23G24G25G30G31
G32
11 2 16 13 20 30.0000000 0 13G01G02G04G10G13G17G20G23G24G25G30G31
G32


11 2 16 13 26 0.0000000 0 13G01G02G04G10G13G16G20G23G24G25G30G31
G32
11 2 16 13 32 0.0000000 0 13G01G02G04G10G13G16G20G23G24G25G30G31
G32



11 2 16 15 27 30.0000000 0 13G01G02G04G05G07G08G10G13G16G20G23G29
G30
11 2 16 15 31 0.0000000 0 13G01G02G04G05G07G08G10G13G16G20G23G29
G30


количество отнаблюдённых спутников GPS составляет 13 штук (расшифровка цифр: 11 – год 2011, 2 – месяц февраль, 16 – число месяца, 15 – часы, 31 – минуты, 0.0000000 – секунды, 0 – неизвестно, 13 – количество отнаблюдённых спутников, G01G02G04G05G07G08G10G13G16G20G23G29 G30 – номера спутников GPS).

(По спутникам ГЛОНАСС – всё в порядке, их общее количество на эпоху регистрации составляет максимум 10 штук, можно проверить, выполнив команду: teqc.exe -O.s R -G +ds 86400 -O.sum s svtl0470.11o > svtl0470_9.11o).


Казалось бы, это хорошо, чем больше наблюдаемых спутников, тем лучше.
Но, ТТС может одновременно за одну эпоху регистрации загружать только 24 спутника, из них 12 – GPS, и 12 – ГЛОНАСС. Остальные спутники ТТС просто отбрасывает и не подгружает их в обработку (интересно, а как это происходит в ПО других производителей?).

Но, к сожалению, отбрасывает спутники без анализа их возвышения над горизонтом, то есть, отбрасывает те, которые стоят в списке последними. А в списке последними стоят спутники с наибольшим по значению номером – по возрастанию нумерации, например, в нашем случае всё время отбрасываются спутники G28, G32 и G30. Хотя в эти эпохи регистрации они находятся гораздо выше над горизонтом, чем другие спутники – имеющие меньшую по значению нумерацию.
В итоге мы имеем картинку, где видно, что непрерывность наблюдений спутников G28, G32 и G30 необоснованно нарушена, что не есть хорошо. Возможно, это не очень повлияет на качество результатов вычислений, но всё же, правильнее было бы, если бы период их наблюдений был непрерывным.

Итак, вопрос: как с помощью TEQC или другой ПО, выполнить отсечку тех спутников, которые имеют низкое возвышение над горизонтом, чтобы общее число наблюдаемых GPS спутников на эпоху регистрации было не более 12 штук.

 

Для Teqc:
1. teqc +qc ***.11o (teqc +qc -set_mask # ***.11o, где # - угол отсечки);
2. teqc -G# ***.11o > new***.11o ,
где # - номер спутника с наименьшим количеством обсерваций над углом отсечки из файла ***.11S.

 

Гы, указанная задача, если я правильно понял, несколько сложнее.

В общем виде:
1. Выкинуть наблюдения со спутников ГЛОНАСС из файла.
2. Выбрать эпохи измерений в которых число наблюдаемых GPS спутников больше 12.
3. В выбранных эпохах обрезать спутники с наименьшим углом возвышения, так, чтобы их общее число осталось 12.
4. По возможности соблюсти непрерывность наблюдений.
5. Вставить наблюдения с ГЛОНАСС спутников обратно в файл .

Все, кроме пункта 4.,вроде можно сделать силами teqc. Надо попробовать.

Блин, teqc пирожки печь случайно не умеет?

 

Значит и решение ее будет сложнее
Например,
1. teqc +qc ***.11o (teqc +qc -set_mask # ***.11o, где # - угол отсечки);
В полученном файле ***.11S выбираем эпохи с числом GPS спутников больше 12.
2. Выделяем необходимую(ые) эпохи в отдельный файл(ы):
teqc -e str1 ***.11o > gap1***.11o
teqc -st str1 -e str2 ***.11o > gap2***.11o
teqc -st str2 ***.11o > gap***.11o ,
где str1 и str2 - время начала и конца эпохи с числом GPS спутников больше 12.
3. teqc +qc gap2***.11o (teqc +qc -set_mask # gap2***.11o, где # - угол отсечки);
4. teqc -G# gap2***.11o > newgap2***.11o ,
где # - номер спутника с наименьшим количеством обсерваций над углом отсечки из файла gap2***.11S
5. Объединяем полученные файлы:
teqc gap1***.11o newgap2***.11o gap3***.11o > new***.11o

А здесь какая собака зарыта ???

 

Я имел в ввиду, что в соседних эпохах лучше резать спутник(ки) с одинаковыми PRN.


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

I:\tmp>findstr /b /r ..............................[12][3456789] svtl0470_8.11o

Пример ниже:
(Добавление)

 

Пример для исправленной строки

 

А у меня чего не так при удалении спутника ?

Изящно...
Ну и приведите сразу Ваш вариант решения до конца...

Оффтоп

У меня в предыдущем посте фигурирует слово "эпоха", его нужно читать как "период времени (с эпохами содержащими более 12 спутников)".

 

У вас все верно. Я сразу не допер. ::tired:: Показалось, что -set_mask режет измерения по маске, а это же qc режим. Сперва просматриваем, потом выключаем не понравившееся.

И все же лучше так ::concordant::
findstr /b /r ..............................1[3456789] svtl0470_8.11o

 

Да, действительно, изящно, Viet Nam.

Кое-что уже получается, напрашивается создание bat-файла.

Например, Poisk_bolsche12.bat, с выводом в отдельный файл bolsche12-GPS.txt результатов поиска.


Итак, оставляем в папке только исходный файл svtl0470.11o и teqc.exe.
Запускаем bat-файл, получаем три файла

svtl0470_G.11o
svtl0470_R.11o
bolsche12-GPS.txt

В файле bolsche12-GPS.txt в начале строк записаны номера строк, где встречаются искомая информация.

20658: 11 2 16 7 31 30.0000000 0 13G03G06G08G11G14G15G16G18G19G21G22G27
20686: 11 2 16 7 32 0.0000000 0 13G03G06G08G11G14G15G16G18G19G21G22G27
20714: 11 2 16 7 32 30.0000000 0 13G03G06G08G11G14G15G16G18G19G21G22G27
23316: 11 2 16 8 25 0.0000000 0 13G03G06G09G11G14G15G18G19G22G24G27G28
23344: 11 2 16 8 25 30.0000000 0 13G03G06G09G11G14G15G18G19G22G24G27G28
и т.д.

Следующей строчки, где записан 13-й спутник - нет, но это и не нужно.

Как теперь использовать информацию из файла bolsche12-GPS.txt для дальнейшего продвижения к цели?

 

Прежде чем вплотную начать играться с +qc -set_mask #, тупой вопрос:

А как teqc узнает данные об угле места спутника, имея только файл OBS? Или каким то образом рядом должен находится NAV файл?

Посмотрел на вывод +qc ::shock:: Пока ушел читать ::rtfm.gif::

 

Рядом должен быть навигационный файл !

 

Используйте RTKPLOT из комплекта программ RTKLIB для визуализации QC .
TEQC интегрирован в ПО RTKLIB.

 

Лично меня не впечатлили возможности RTKPLOT для контроля качества первичных измерений

 

Viet Nam, сколько можно читать?
За прошедший период времени кое-что удалось «накопать».

Сделан ещё один шаг к достижению цели.
Получился bat-файл версии 1.1_

Чтобы проверить его работу, выполняем следующее:
1. Создаём новую папку.
2. Скачиваем упакованный файл.
3. Распаковываем в папку.
4. Запускаем bat-файл.
5. Откидываемся на спинку стула, расслабляемся, любуемся картинками

Game Over, переходим к анализу информации.

Полученные файлы:
_QC\svtl0470_G.ele
_SKYPLOT\svtl0470_G.elr
содержат искомую информацию о возвышениях спутников над горизонтом. Первый файл создан программой TEQC, другой – SKYPLOT, почему-то файлы немного отличаются.

Их отличия можно посмотреть здесь:
_ele-QC.txt
_ele-SkyPlot.txt

и сравнить с
_bolsche12-GPS.txt

Что дальше, предлагайте. http://geodesist.ru/forum/tools.php?action=smiles#

P.S. Навигационного файла для svtl0470 найти не удалось, поэтому взял brct0470.11o.
Если необходим анализ наблюдений для Ваших других файлов - можете их временно переименовать (файл****.11o - svtl0470.11o, файл****.11n - brct0470.11n) и запустить bat-файл на выполнение в новой папке. Туда-же нужно скопировать и файлы *.exe.

 

Оффтоп

Если позволите буду без сценария батника, так сказать только идеей...

Rodik,
В предыдущий Ваш bat-файл добавьте :
1. Поиск в файле ***.11o (исходный) строки с "INTERVAL", из которой получаем интервал записи эпох (а можно и жестко задать как 30 сек);
2. В цикле из каждой строки файла bolsche12-GPS.txt берем секунды и сравниваем с секундами предыдущей строкой. Их значения не должны быть больше интервала записи эпох. Т.о. получаем интервалы времени с эпохами, в которых количество спутников было больше 12.
3. Разбиваем исходный файл на полученные интервалы...

Пока так. Это реально сделать ?

 

Извиняюсь, но мой больничный закончился, а текучка уже заела

Если все же немножко уйти от teqc, но продолжать в духе минимализма, то можно использовать Skyplot из набора GPS-toolbox
Зная нужный интервал формируем файл skyplot.inp

Код:

11  2 16  7 31 30.0
11  2 16 15 31  0.0 
brdc0470.11n
svtl 2730153.4307  1562363.9908  5529988.4491
10

запускаем skyplot.exe & skyplot.bat
и оцениваем картинку skyplot.ps

 

Зная нужный интервал прежде всего необходимо будет выделить его (их) из исходного файла (т.е. разбить исходный файл).
После чего в файле ***.11S (для полученного файла выделенного интервала) достаточно найти спутник с наименьшим количеством записей и выключить его...

 

Так наверное все же не прокатит. Рассмотрим пример.
Провели проверку и выделили непрерывный интервал (в кавычках конечно, т.е интервал где разница во времени между соседними эпохами не превышает, в нашем случае, 30 сек)


---------------------------Восход G#1+++++++++++++++++++
++++++++++++++++ЗаходG#2-----------------------------------

t0_____________________t1__________________________t2


Выключаем G#2 (допустим у него на интервале наименьшее количество записей)
Остается интервал t1-t2, где число наблюдаемых спутников больше 12.

(извините без цитат, на работе скрипты отключены)
(Добавление)
Блин, проглядел,
со Skyplot Rodik опередил, ::смущен::
Финишный батник еще не тестировал, обязательно сделаю это на выходных.