Acessibilidade / Reportar erro

PPP COM INTEGRAÇÃO DE DADOS GPS/GLONASS: FUNDAMENTOS ENVOLVIDOS E ANÁLISE DE ACURÁCIA

PPP with GPS/GLONASS data Integration: Fundamentals Involved and Accuracy Assessment

Resumos

Atualmente o GNSS possui o GPS e GLONASS operacionais, o que possibilita a integração destes sistemas no posicionamento geodésico. A integração de dados GPS/GLONASS fornece melhor geometria e maior confiabilidade na estimativa dos parâmetros envolvidos. Ao aplicar a integração no posicionamento é necessário realizar a compatibilização do sistema de tempo e também do Sistema Geodésico de Referência no caso da utilização de órbitas transmitidas. Em termos de posicionamento geodésico especial atenção tem sido dada ultimamente para o método PPP, o qual requer um único receptor no nível de usuário, órbitas e correções precisas de relógios dos satélites, além da modelagem matemática que permita corrigir os diversos efeitos envolvidos com a propagação dos sinais GNSS e efeitos geodinâmicos que afetam a estação. Neste artigo apresentam-se os modelos matemáticos do ajustamento no modo PPP com integração de dados GPS/GLONASS, além da análise de acurácia posicional envolvendo dados coletados por receptores em estações brasileiras. Os processamentos de dados foram realizados a partir do software GPSPPP do NRCAN com soluções diárias e a geração de série temporal anual. As análises foram realizadas em termos de Erro Médio Quadrático diário e foram observadas melhorias de até 25% na estimativa da posição ao aplicar os dados GPS/GLONASS no PPP ao invés de somente dados GPS.

Palavras-chave:
GNSS; Integração GPS e GLONASS; PPP


Currently GNSS has GPS and GLONASS operating, which enables the integration of these systems in geodetic positioning. The integration of GPS/GLONASS data provides better geometry and greater reliability in the estimation of the involved parameters. By applying the integration in the positioning it is necessary to perform the compatibility of the time system and also of the Geodetic Reference System, in case of applying broadcasted orbits. In terms of geodetic positioning special attention has been given nowadays for the PPP method that requires a single receiver at the user level, precise orbit and clock corrections for satellites, as well as mathematical modeling capable to correct for several effects involved in the GNSS signals propagation and the geodynamic effects in the station. In this paper it will be presented the adjustment mathematical model for PPP with integration of GPS/GLONASS data, as well as analysis of positional accuracy involving data collected by receivers in Brazilian stations. The data processing was carried out by using the GPSPPP NRCan software considering daily solutions with generation of annual time series. Analyses were performed in terms of daily Root Mean Square and improvements were observed reaching up to 25% in the position estimation when applying GPS/GLONASS in the PPP instead of only GPS data

Keywords:
GNSS; GPS and GLONASS integration; PPP


1. Introdução

O posicionamento geodésico atualmente é realizado em sua maioria a partir do GNSS (Global Navigation Satellite System), podendo-se destacar o GPS (Global Positioning System) e o GLONASS (GLObal'naya NAvigatsionnaya Sputnikkovaya Sistema), os quais se encontram atualmente em sua forma operacional. Outros sistemas, tais como o GALILEO (European Global Navigation Satellite System) e o COMPASS/BeiDou (CNSS - Compass Navigation Satelllite System) se encontram em fase de desenvolvimento. Dentre os sistemas de navegação, o GPS ainda é o mais utilizado na atualidade considerando que este tem sido aplicado como sistema de navegação em sua forma operacional completa por aproximadamente três décadas. Contudo, o GLONASS atualmente também se encontra em sua forma totalmente operacional com uma constelação de satélites completa.

Atualmente, o GPS e o GLONASS passam por uma fase de modernização, a qual envolve o segmento espacial e o de controle. A modernização do GLONASS contempla a transmissão de sinais com base na técnica CDMA (Code Division Multiple Access) em três portadoras e a modernização do GPS envolve a inclusão do novo código civil na portadora L2 (L2C), o código civil na L1 (L1C), além da nova portadora L5 e sinais militares. (Hofmann-Wellenhof et al., 2008Hofmann-Wellenhof, Bernhard with Lichtenegger, H. and Wasle, E. GNSS: Global Navigation Satellite Systems Springer Wien New York: 2008.; ICD-GLONASS, 2008ICD-GLONASS Interface Control Document: Navigation radiosignal in bands L1, L2 Edition 5.1. 2008. Moscow.). O sinal L2C e a portadora L5 já estão disponíveis a partir dos satélites GPS do bloco IIR e IIF. Estes fatores favorecem a interoperabilidade dos sistemas de navegação abrindo oportunidades seja para aplicações de ordem prática ou de pesquisa científica.

Em relação ao posicionamento destaca-se método PPP (Posicionamento por Ponto Preciso), o qual tem sido investigado pela comunidade cientifica nos últimos anos. Este método quando aplicado em sua forma convencional requer o uso de órbitas e correções precisas de relógios dos satélites, aplicação da observável livre da ionosfera (ion-free) para a pseudodistância e fase da onda portadora, além da modelagem matemática dos efeitos envolvidos visando obtenção de coordenadas com acurácia centimétrica. A estimativa das coordenadas pode ser no modo PPP pós-processado ou em quase tempo real (latência de poucos segundos) (Marques et al., 2014Marques, Haroldo. A. with Monico, João F. G., Shimabukuro, Milton H., Oyama, R. T. and Wentz, Jens P. "PPP Em Tempo Real: Fundamentos, Implementação Computacional e Análises de Resultados no Modo Estático e Cinemático." Revista Brasileira de Cartografia 66, no 6 (2014): 1331-45.). Os efeitos envolvidos estão relacionados principalmente com a propagação dos sinais GNSS, como por exemplo, aqueles provocados pela ionosfera e troposfera, multicaminho, erro da órbita, erros de relógios do receptor e satélites, efeitos de relatividade, atraso de hardware, centro de fase da antena do satélite e do receptor, entre outros. (Seeber, 2003Seeber, Gunter. Satellite geodesy: foundations, methods, and applications Berlin, New York: Walter de Gruyter, 2003.; Hofmann-Wellenhof et al., 2008Hofmann-Wellenhof, Bernhard with Lichtenegger, H. and Wasle, E. GNSS: Global Navigation Satellite Systems Springer Wien New York: 2008.; Monico, 2008Monico, João F. G. Posicionamento pelo GNSS: descrição, fundamentos e aplicações. São Paulo: Unesp, 2008.).

Neste trabalho, apresentam-se os modelos matemáticos envolvidos com integração GPS/GLONASS no PPP e análise da acurácia posicional envolvendo dados GNSS coletados em estações da RBMC (Rede Brasileira de Monitoramento Contínuo dos sistemas GNSS). Para análise de acurácia foram realizados processamentos de dados no modo PPP com o software GPSPPP do NRCAN (Natural Resources Canada), onde, primeiramente foram utilizadas somente medidas GPS e em seguida medidas GPS/GLONASS com soluções diárias e a geração de série temporal anual. As coordenadas estimadas, seja aplicando somente GPS ou GPS/GONASS, foram comparadas com as coordenadas conhecidas das estações e análises foram realizadas em termos de Erro Médio Quadrático (EMQ). Análise da geometria posicional também foi realizada e, em geral, verificam-se melhorias de até 25% na acurácia posicional ao aplicar os dados GPS/GLONAS no PPP ao invés de somente GPS.

2. Modelo Matemático do PPP com integração GPS/GLONASS

Os satélites GPS transmitem sinais gerados a partir de uma frequência fundamental (f0) de 10,23 MHz, cuja multiplicação por números inteiros, 154, 120 e 115 permite formar a frequência de cada uma das três portadoras denominadas de L1 (1575,42 MHz), L2 (1227,60 MHz) e L5 (1176,45 MHz). Os sinais são transmitidos com base na técnica CDMA, no qual todos os satélites podem transmitir numa mesma frequência e serem identificados pelo seu código PRN (PseudoRandom Noise). Diferentemente do GPS, cada satélite GLONASS transmite sinais em frequências específicas. Desta forma, a identificação do satélite se dá pela frequência do sinal, técnica denominada FDMA (Frequency Division Multiple Access) (Seeber, 2003Seeber, Gunter. Satellite geodesy: foundations, methods, and applications Berlin, New York: Walter de Gruyter, 2003.; Leick, 2004Leick, Alfred. "GPS satellite surveying". New York: John Wiley & Sons, 2004.; Monico, 2008Monico, João F. G. Posicionamento pelo GNSS: descrição, fundamentos e aplicações. São Paulo: Unesp, 2008.). Contudo, o GLONASS já possui atualmente dois satélites (GLONASS-K) transmitindo dados com base na técnica CDMA. As frequências em cada satélite GLONASS são definidas por (Roosbach, 2001Rossbach, U. Positioning and Navigation Using the Russian Satellite System GLONASS. PhD Thesis. 2001. Universität der Bundeswehr München, Neubiberg.; Hofmann-Wellenhof et al., 2008Hofmann-Wellenhof, Bernhard with Lichtenegger, H. and Wasle, E. GNSS: Global Navigation Satellite Systems Springer Wien New York: 2008.):

As frequências f1, f2 e f3 estão, respectivamente, associadas às portadoras L1, L2 e L3. O valor de k diferencia os canais de frequência. Os termos Δf1, Δf2 e Δf3 são fatores que denotam incrementos na frequência de sinais dos satélites em dois canais adjacentes (Hofmann-Wellenhof et al., 2008Hofmann-Wellenhof, Bernhard with Lichtenegger, H. and Wasle, E. GNSS: Global Navigation Satellite Systems Springer Wien New York: 2008.). As equações de pseudodistância (PD) e de fase da onda portadora para o GPS e GLONASS em unidades de metro podem ser escritas como (Seebeer, 2003Seeber, Gunter. Satellite geodesy: foundations, methods, and applications Berlin, New York: Walter de Gruyter, 2003.; Monico, 2008Monico, João F. G. Posicionamento pelo GNSS: descrição, fundamentos e aplicações. São Paulo: Unesp, 2008.):

Nas Equações 4 a 7, os símbolos G e R são usados para representar respectivamente o GPS e o GLONASS. A distância geométrica entre o satélite e o receptor é representada por ρ e os erros dos relógios do receptor e do satélite são representados, respectivamente, por dtrec e dtsat . A ambiguidade da portadora é denotada por e os erros da ionosfera, troposfera, órbita e multicaminho são representados, respectivamente, por ILi e T, Orb e dm. Os termos , ,representam erros aleatórios ou não modelados. O subscrito Li se refere à frequência da portadora.

No caso do PPP, a posição e o erro do relógio do satélite são, geralmente, obtidos a partir das efemérides precisas e arquivos de correções de relógios, os quais atualmente são disponibilizados pelos centros do IGS tanto para o GPS quanto para o GLONASS. Neste caso, o SGR (Sistema Geodésico de Referência) ao qual a órbita está vinculada é um dos ITRFs (International Reference Frame), atualmente o ITRF2008. No caso da utilização de efemérides transmitidas, é importante destacar que é necessária a contabilização do SGR, uma vez que a órbita transmitida do GPS é dada no WGS84 e a do GLONASS no PZ90. Enquanto a órbita do GPS é disponibilizada em forma de elementos Keplerianos, a do GLONASS é dada em forma de coordenadas cartesianas, velocidades e acelerações, requerendo um processo de integração numérica do modelo dinâmico de órbita (Stewart; Tsakiri, 1998Stewart, Mike and Tsakiri, Maria. "GLONASS broadcast orbit computation." GPS Solutions2, (1998): 16-27.).

O erro do relógio do receptor é calculado como a diferença entre o instante de recepção (tr) no receptor e sistema de tempo GPS. O mesmo ocorre para o erro do relógio no caso do GLONASS, porém, este é obtido com base na diferença entre o tempo de recepção no receptor em relação ao sistema de tempo GLONASS (ver Equações 6 e 7). Desta forma, o erro do relógio do receptor pode ser representado por (Cai; Gao, 2007Cai, Changseng and Gao, Yang. "Precise Point Positioning Using Combined GPS and GLONASS Observations." Journal of Global Positioning System 6, (2007): 13-22.):

O termo tsys denota o sistema de tempo GPS (TGPS) para observações GPS ou o sistema de tempo GLONASS (TGLONASS) para observações GLONASS. Assim, em termos de ajustamento no modo PPP com integração GPS/GLONASS, comparecem dois parâmetros relacionados ao erro de relógio do receptor. Por esta razão, utiliza-se a estratégia de estimar o atraso entre os sistemas de tempo GPS e GLONASS, o qual é denotado por dtsys. Fazendo uso da Equação 8 para o caso do erro do relógio do receptor em relação ao GLONASS, tem-se (Cai; Gao, 2007Cai, Changseng and Gao, Yang. "Precise Point Positioning Using Combined GPS and GLONASS Observations." Journal of Global Positioning System 6, (2007): 13-22.):

Verifica-se na Equação 9 que o erro do relógio do receptor em relação ao tempo GLONASS () é calculado em função do erro do relógio em relação ao tempo GPS () e em relação ao atraso entre os sistemas de tempo GPS e GLONASS (dtsys). Aplicando a Equação 9 nas Equações 6 e 7, tem-se:

Para usuários com receptores de dupla frequência, os efeitos ionosféricos de primeira ordem podem ser praticamente eliminados através da combinação ion-free (livre da ionosfera). Os efeitos ionosféricos de ordem superior podem ser corrigidos conforme descrito em Marques et al., (2011Marques, Haroldo. A. with Monico, João F. G.and Aquino, M. "RINEX_HO: Second and third-order ionospheric corrections for RINEX observation files." GPS Solutions 15, (2011): 305-314.). As equações da combinação ion-free para a pseudodistância e fase do GPS e GLONASS são dadas por:

Nas Equações 12 a 15 têm-se como parâmetros as coordenadas da estação, erro do relógio do receptor em relação ao tempo GPS, diferença do sistema de tempo (dtsys), as ambiguidades da fase e a componente úmida da troposfera (Tzw) na direção zenital sendo que mf representa sua função de mapeamento (Monico, 2008Monico, João F. G. Posicionamento pelo GNSS: descrição, fundamentos e aplicações. São Paulo: Unesp, 2008.; Marques et al., 2014Marques, Haroldo. A. with Monico, João F. G., Shimabukuro, Milton H., Oyama, R. T. and Wentz, Jens P. "PPP Em Tempo Real: Fundamentos, Implementação Computacional e Análises de Resultados no Modo Estático e Cinemático." Revista Brasileira de Cartografia 66, no 6 (2014): 1331-45., Sapucci, 2001Sapucci, Luíz. F. "Estimativa do vapor d'água atmosférico e avaliação da modelagem do atraso zenital troposférico utilizando GPS". diss., Universidade Estadual Paulista, 2001.). As órbitas e erros de relógios dos satélites advêm das efemérides precisas e outros erros envolvidos devem ser modelados. Desta forma, as equações de pseudodistância e fase em sua forma linearizada (Gemael, 1994Gemael, Camilo. Introdução ao ajustamento de observações: Aplicações Geodésicas Curitiba: Editora UFPR, 1994.) podem ser escritas como:

Onde o termo Δ é usado para representar as diferenças entre valores observados e calculados, ou seja, o vetor LXS. Os termos representam as coordenadas dos satélites, seja GPS ou GLONASS (i=G ou R); ΔXrec, ΔYrec e ΔZrec representam as correções das coordenadas aproximadas das estações(). Designando as derivadas parciais relacionadas às coordenadas por podem-se reescrever as Equações 16 a 19 em forma matricial:

A Matriz de Variância Covariância (MVC) das observações (∑Lb) no caso do PPP com integração GPS e GLONASS pode ser escrita da seguinte forma:

No caso da Equação 21, consideram-se as observações não correlacionadas. Estratégias de ponderação podem ser adotadas conforme descritas em Wang, (1999Wang, Jinling. "Modelling and quality control for precise GPS and GLONASS satellite positioning" PhD diss., Curtin University of Technology, 1999.) e Silva, (2009)Silva. Heloísa A. "Avaliação de modelos estocásticos no posicionamento GNSS." diss., Universidade Estadual Paulista, 2009., como por exemplo, em função do ângulo de elevação do satélite. Os experimentos e análises realizadas serão apresentados na próxima seção.

3. Resultados e análises

Para a realização dos experimentos e análises foram utilizados dados GNSS coletados em estações pertencentes à RBMC e distribuídas em diferentes latitudes ao longo do Brasil. A estações utilizadas, foram BOAV (Boa Vista - RR), GOJA (Jataí - GO), PPTE (Pres. Prudente - SP), RECF (Recife-PE) e UFPR em Curitiba - PR, cuja localização aproximada pode ser vista na Figura 1.

Figura 1
Localização aproximada das estações utilizadas nos experimentos

Os dados foram processados com o software GPSPPP do NRCAN licenciado para o DECART (Departamento de Engenharia Cartografia da UFPE) para fins de pesquisa. Além disto, utilizou-se o software TEQC (Translation Edition Quality Control) (UNAVCO, 2015Unavco. "Teqc - Translation Edition Quality Control.", accessed April, 2015. http://www.unavco.org/software/data-processing/teqc/teqc.html
http://www.unavco.org/software/data-proc...
) para a retirada de observações GLONASS do arquivo RINEX (Receiver Independent Exchange Format) com a geração de um novo arquivo. Assim primeiramente, foram processados os dados somente GPS e em seguida GPS/GLONASS.

Para a realização dos experimentos, as precisões adotadas para as medidas de pseudodistância e de fase foram, respectivamente de 1 m e 0,01 m. Além disto, foi adotado fator de variância a priori unitário. Foram utilizadas órbitas precisas finais (formato .sp3) e correções finais dos relógios dos satélites (formato .clk com atualização em intervalos de 30 s). Além disto, foram aplicadas correções absolutas de variação de centro de fase das antenas, carga de marés oceânicas, entre diversas outras correções disponíveis no software GPSPPP. A solução do PPP foi realizada no modo forward do ajustamento recursivo. Ao final do processamento PPP, seja usando somente dados GPS ou usando GPS/GLONASS, as coordenadas estimadas foram comparadas com as coordenadas de referência obtidas a partir da rede SIRGAS-CON (solução em ITRF 2008) com a devida atualização para a época da coleta dos dados através da solução multianual SIR1101P01 e modelos de velocidades VEMOS2009 (Velocity Model for SIRGAS 2009). Os erros em coordenadas cartesianas foram transformados para o Sistema Geodésico Local (SGL) com base nas equações apresentadas em Torge, (2001Torge, Wolfgang. Geodesy Berlim, New York: Walter de Gruyter, 2001.) e Seeber. (2003Seeber, Gunter. Satellite geodesy: foundations, methods, and applications Berlin, New York: Walter de Gruyter, 2003.). Neste caso, o sistema de eixos do SGL é centrado na latitude e longitude média de todo o período de processamento e as análises foram realizadas em termos de discrepância nas coordenadas referenciadas ao SGL, ou seja, DE, DN e DU. As discrepâncias são doravante denominadas de 'erros' no SGL. A análise de acurácia foi realizada em termos de EMQ, o qual pode ser tomado como um medidor de acurácia como pode ser visto em Monico et al., (2009Monico, João F. G. with Póz, Aluir P. D., Galo, Maurício, Santos, Marcelo C. and Oliveira, Leonardo C. "Acurácia e Precisão: Revendo os Conceitos de Forma Acurada." Boletim de Ciências Geodésicas15, no 3 (2009): 469-83.). Os dados das cinco estações (Figura 1) foram processados em soluções diárias com a geração de séries temporais das componentes do SGL para o ano de 2013. No caso de análise anual, a origem do sistema de eixos do SGL (DE, DN e DU) é centrada na coordenada média de todo o período de processamento para cada estação.

3.1 Análise de precisão e DOPs com integração GPS/GLONASS

Considerando que se tenha disponível as coordenadas da estação e as órbitas dos satélites, pode-se calcular a MVC de parâmetros dada por:, em que é o fator de variância a posteriori. A matriz Jacobiana (A) e a matriz dos pesos (P) podem ser obtidas com base nas Equações 20 e 21. Assumindo que o ajustamento seja aceito no teste Qui-quadrado, o fator de variância a posteriori pode ser considerado estatisticamente igual ao fator de variância a priori, o qual pode ser arbitrado no ajustamento. Assim, com base na MVC de parâmetros pode-se analisar a priori a geometria (DOP - Diluition of Precision) do GPS e GLONASS e a precisão estimadas dos parâmetros.

Para o caso deste experimento utilizou-se as coordenadas da estação RECF e a posição dos satélites foi obtida a partir das efemérides precisas considerando uma época escolhida de forma arbitrária (0h 0m 15 seg GPS do dia 09/06/2014). Foram utilizados dados de órbita de oito satélites GPS (PRN: 01, 08, 11, 15, 19, 23, 27 e 32) e seis satélites GLONASS (Número: 02, 06, 13, 15, 18 e 24). Após a obtenção da MVC foram calculados os valores de PDOP (Positioning DOP), GDOP (Geometry DOP) e ADOP (Ambiguity DOP). No caso do posicionamento por ponto o cálculo de DOP pode ser realizado com base na MVC propagada para o SGL, cuja formulação pode ser encontrada, por exemplo, em Seeber, (2003Seeber, Gunter. Satellite geodesy: foundations, methods, and applications Berlin, New York: Walter de Gruyter, 2003.) ou Hofmann-Wellenhof et al., (2008Hofmann-Wellenhof, Bernhard with Lichtenegger, H. and Wasle, E. GNSS: Global Navigation Satellite Systems Springer Wien New York: 2008.). No caso do ADOP quanto menor seu valor maior a probabilidade de se obter a solução inteira das ambiguidades (Teunissen et al.; 1999Teunissen, Peter. J. G with Odijk, Dennis and De Jong, C. D. "Ambiguity Dilution of Precision: An additional tool for GPS quality control." Department of Mathematical Geodesy and Positioning. Thijsseweg, 1999.; Monico et al., 2014Monico, João F. G. with Rofatto, Vinícius F. and Mendonça, Marco A. M. "Introdução e aplicação diluição da precisão das ambiguidades GNSS - ADOP"., Boletim de Ciências Geodésicas 20 no 2 (2014): 279-99.). Os valores das precisões obtidas para os parâmetros envolvidos são apresentados na Tabela 1 juntamente com o nível de melhoria em termos percentuais ao aplicar GPS/GLONASS ao invés de somente GPS.

Tabela 1:
Estimativa das precisões (em metros) e melhorias percentuais

Verifica-se na Tabela 1 que ao utilizar dados GPS/GLONASS a melhoria percentual na precisão das coordenadas foi de aproximadamente 25%, 37% e 22%, respectivamente para X, Y e Z. As precisões das ambiguidades GPS também apresentaram melhorias com a aplicação dos dados GPS/GLONASS sendo o máximo de aproximadamente 11% para o caso do PRN08. Os valores de PDOP, GDOP e ADOP obtidos a partir da MVC calculada para o caso apenas GPS e GPS/GLONASS são apresentados na Tabela 2, juntamente com o nível de melhorias ao usar a integração dos sistemas.

Tabela 2:
DOPs calculados e melhorias esperadas ao usar a integração

Na Tabela 2 verifica que os valores de PDOP, GDOP e ADOP ao inserir dados da constelação GLONASS juntamente com GPS apresentam, respectivamente, melhorias de aproximadamente 30%, 14% e 3%. Com base neste experimento, pode-se dizer que o nível máximo de melhoria esperado nas precisões da estimativa das coordenadas ao inserir dados da constelação GLONASS é de aproximadamente 30%. No caso das ambiguidades, o nível de melhoria obtido foi no máximo de aproximadamente 11%, contudo é importante destacar que o experimento foi realizado para uma única época.

Com o objetivo de avaliar a variação diária do GDOP no PPP, os dados GPS e GLONASS coletados na estação RECF foram processados numa solução diária aplicando a estratégia de processamento usando somente dados GPS e em seguida GPS/GLONASS. A Figura 2 apresenta os valores de GDOP e número de satélites para cada um dos processamentos realizados.

Figura 2:
GDOP e nº de satélites disponíveis nas soluções GPS e GPS/GLONASS

A partir da Figura 2 verifica que os valores de GDOP para o caso GPS/GLONASS apresentam, em geral, valores melhores, principalmente, para os picos próximos das 18 e 21 horas UTC onde o número de satélites GPS disponíveis foi de aproximadamente seis e atingindo 14 satélites com a inserção das medidas GLONASS.

A convergência do PPP (tempo para as precisões dos parâmetros se estabilizarem) pode ser analisada com base na série temporal das precisões estimadas das coordenadas. No caso, foram processados dados da estação RECF para cinco dias da semana GPS 1830 (Fevereiro de 2015) considerando o caso do PPP cinemático. Foram analisadas somente as primeiras quatro horas de dados (961 épocas), seja com solução GPS ou GPS/GLONASS. A série temporal das precisões resultantes (3D) neste caso pode ser vista na Figura 3.

Figura 3:
Precisões 3D no PPP cinemático com solução GPS e GPS/GLONASS

Verifica-se na Figura 3 que, em geral, ao aplicar a integração de dados GPS/GLONASS no PPP, as precisões das coordenadas apresentaram melhor convergência nas primeiras épocas de processamentos. Considerando os dados coletados e processados relativos às quatro primeiras horas de dados, o tempo de convergência obtido foi em média de 2 horas (480 épocas) quando utilizados apenas dados GPS, enquanto que com a utilização de dados da integração GPS/GLONASS, o tempo de convergência médio foi de 1h e 10 min (280 épocas).

3.2 Análise da acurácia no PPP com integração GPS/GLONASS

Os dados GNSS das cinco estações foram processados no modo PPP estático com solução diária para o ano de 2013, cuja série temporal dos erros no SGL pode ser vista na Figura 4.

Figura 4:
Discrepâncias no SGL (ano 2013) para as cinco estações analisadas

Analisando a Figura 4 verificam-se em geral valores menores nas componentes do SGL quando aplicada a integração GPS/GLONASS em relação à aplicação de somente GPS. Geralmente a componente altimétrica é a mais afetada e nota-se um comportamento cíclico ao longo do ano nas séries temporais com maiores valores próximos ao período dos solstícios, 21 de Junho (dia 172) e 21 de Dezembro (dia 355). O efeito é mais intenso para a estação BOAV, situada em Roraima, região amazônica. As altitudes das estações da referida região sofrem uma consideravelmente variação anual atingindo a ordem de centímetros devido à força exercida pelo acúmulo de água na bacia Amazônica. Este efeito sazonal também pode ser observado em estações das regiões Centro-oeste e Sudeste do País (Costa et al., 2012Costa, Sonia M. A. with Silva, Alberto L., Lima, Marco A. A. and Júnior, Newton J. M. "Centro de análise SIRGAS - IBGE: Novas estratégias de processamento e combinação, e influencia da mudança de referencial global nos resultados." Boletim de Ciências Geodésicas, 18, no 1 (2012): 65-83.).

Em geral, os erros nas componentes planimétrica e altimétrica não ultrapassaram valores de 0,02 m e 0,04 m, respectivamente. Na Figura 5 são apresentadas as séries temporais dos erros em planimetria (2D) e resultante (3D) calculados para cada uma das estações utilizadas nos experimentos.

Figura 5:
Erros 2D e 3D (ano 2013) para as cinco estações analisadas

A partir da Figura 5 verifica-se que os erros 2D e 3D, em geral apresentaram menores valores ao aplicar a integração GPS/GLONASS no PPP. No caso 2D, os valores máximos nas duas soluções atingem a ordem de 0,02 m e no caso 3D aproximadamente 0,04 m. As diferenças nas séries temporais ao aplicar a integração em relação ao uso de somente GPS são mais aparentes para a estação BOAV. Com base nos valores diários obtidos, foi calculado o EMQ considerando todo o período do ano de 2013 (EMQ anual), seja para a solução PPP usando somente GPS ou GPS/GLONASS. Os valores de EMQ anual são apresentados na Tabela 3, juntamente com o nível de melhoria ao aplicar a integração GPS/GLONAS em relação ao uso de somente GPS.

Tabela 3:
EMQ anual e melhorias para cada uma das estações

Ao analisar a Tabela 3, nota-se que o EMQ anual ao aplicar a integração de dados GPS/GLONASS no PPP apresenta melhorias em relação à aplicação de somente dados GPS. A máxima melhoria em EMQ 3D foi para a estação BOAV com valor de aproximadamente 24,5% e a mínima ocorreu para a estação GOJA com valor de 0,74%. No caso da estação GOJA não houve melhorias na componente planimétrica. O nível de melhoria com a componente altimétrica (EMQ 3D), para cada estação, foi aproximadamente de, 24,5%, 0,7%, 16,8%, 10,9% e 17,4%, respectivamente para as estações BOAV, GOJA, PPTE, RECF e UFPR.

4. Conclusões

Este trabalho teve como principal objetivo apresentar os fundamentos matemáticos envolvidos com o ajustamento de dados no modo PPP com a integração de dados do GPS/GLONASS, além de avaliar a acurácia posicional ao aplicar a integração. Para o ajustamento dos dados no PPP, verifica-se que comparece um atraso entre os sistemas de tempo, o qual deve fazer parte dos parâmetros a serem estimados.

As análises de geometria e da precisão esperada para o caso GPS e GPS/GLONASS foram realizadas a princípio com base na MVC dos parâmetros, a qual requer o conhecimento da órbita dos satélites e da posição de uma estação para que se construa a matriz A. Verificou-se que a inserção de medidas GLONASS no posicionamento pode proporcionar melhorias na estimativa das precisões da ordem máxima de até 36% (ver Tabela 1), haja vista que o número de satélites, geralmente, aumenta em média aproximadamente 60% por época. Ao aplicar a integração, os valores de PDOP, GDOP e ADOP ficaram melhores com valores, respectivamente de aproximadamente, 30%, 14% e 3%, (ver Tabela 2), sendo importante destacar que se trata de análise utilizada em planejamento de coletas de dados GNSS.

Em termos de análise de acurácia do posicionamento, foram processados dados de cinco estações da RBMC distribuídas em diferentes latitudes ao longo do Brasil com solução diária no modo PPP estático e geração de séries temporais para o ano de 2013. Verificou-se um comportamento cíclico anual nas componentes do SGL em cada uma das estações, principalmente para a estação BOAV, a qual é bastante afetada pelas cargas provocadas pelo peso da massa d'água devido à cheia dos rios nos períodos de chuva que ocorrem na região. É importante destacar neste caso que ao aplicar a integração GPS/GLONASS, o efeito cíclico diminui consideravelmente, principalmente na componente altimétrica. Em relação aos erros no SGL verifica em geral valores menores nas componentes do SGL ao aplicar a integração GPS/GLONASS. No caso do EMQ 3D anual, a máxima melhoria foi para a estação BOAV com valor de 24,5% e a mínima ocorreu para a estação GOJA com valor de 0,7% e na média geral a melhoria em EMQ 3D anual foi de 14,06%. Com base nos experimentos e análises realizadas, pode-se dizer que em geral, a inserção de dados GLONASS juntamente com dados GPS no PPP proporciona melhorias em termos de acurácia. A integração dos dados no ajustamento requer técnicas especificas e é necessário considerar que com a modernização de ambos os sistemas, novos sinais estarão disponíveis abrindo oportunidades de desenvolvimento e pesquisa.

AGRADECIMENTOS

Os autores agradecem ao CNPq pelo apoio financeiro em forma de projeto Universal (Processo: 475775/2013-9); a CAPES pelo apoio em forma de bolsa de mestrado junto ao Programa de Pós-Graduação em Ciências Geodésicas e Tecnologias de Geoinformação da UFPE; ao projeto P&D UFPE/UNESP/EMBRAER (Processo: VDT0044-12) e ao projeto P&D ANEEL/CHESF/UFPE (92.2011.5510.00)

REFERÊNCIAS BIBLIOGRÁFICAS

  • Cai, Changseng and Gao, Yang. "Precise Point Positioning Using Combined GPS and GLONASS Observations." Journal of Global Positioning System 6, (2007): 13-22.
  • Costa, Sonia M. A. with Silva, Alberto L., Lima, Marco A. A. and Júnior, Newton J. M. "Centro de análise SIRGAS - IBGE: Novas estratégias de processamento e combinação, e influencia da mudança de referencial global nos resultados." Boletim de Ciências Geodésicas, 18, no 1 (2012): 65-83.
  • Gemael, Camilo. Introdução ao ajustamento de observações: Aplicações Geodésicas Curitiba: Editora UFPR, 1994.
  • Hofmann-Wellenhof, Bernhard with Lichtenegger, H. and Wasle, E. GNSS: Global Navigation Satellite Systems Springer Wien New York: 2008.
  • ICD-GLONASS Interface Control Document: Navigation radiosignal in bands L1, L2 Edition 5.1. 2008. Moscow.
  • Leick, Alfred. "GPS satellite surveying". New York: John Wiley & Sons, 2004.
  • Marques, Haroldo. A. with Monico, João F. G., Shimabukuro, Milton H., Oyama, R. T. and Wentz, Jens P. "PPP Em Tempo Real: Fundamentos, Implementação Computacional e Análises de Resultados no Modo Estático e Cinemático." Revista Brasileira de Cartografia 66, no 6 (2014): 1331-45.
  • Marques, Haroldo. A. with Monico, João F. G.and Aquino, M. "RINEX_HO: Second and third-order ionospheric corrections for RINEX observation files." GPS Solutions 15, (2011): 305-314.
  • Monico, João F. G. Posicionamento pelo GNSS: descrição, fundamentos e aplicações. São Paulo: Unesp, 2008.
  • Monico, João F. G. with Póz, Aluir P. D., Galo, Maurício, Santos, Marcelo C. and Oliveira, Leonardo C. "Acurácia e Precisão: Revendo os Conceitos de Forma Acurada." Boletim de Ciências Geodésicas15, no 3 (2009): 469-83.
  • Monico, João F. G. with Rofatto, Vinícius F. and Mendonça, Marco A. M. "Introdução e aplicação diluição da precisão das ambiguidades GNSS - ADOP"., Boletim de Ciências Geodésicas 20 no 2 (2014): 279-99.
  • Rossbach, U. Positioning and Navigation Using the Russian Satellite System GLONASS. PhD Thesis. 2001. Universität der Bundeswehr München, Neubiberg.
  • Sapucci, Luíz. F. "Estimativa do vapor d'água atmosférico e avaliação da modelagem do atraso zenital troposférico utilizando GPS". diss, Universidade Estadual Paulista, 2001.
  • Seeber, Gunter. Satellite geodesy: foundations, methods, and applications Berlin, New York: Walter de Gruyter, 2003.
  • Silva. Heloísa A. "Avaliação de modelos estocásticos no posicionamento GNSS." diss, Universidade Estadual Paulista, 2009.
  • Stewart, Mike and Tsakiri, Maria. "GLONASS broadcast orbit computation." GPS Solutions2, (1998): 16-27.
  • Teunissen, Peter. J. G with Odijk, Dennis and De Jong, C. D. "Ambiguity Dilution of Precision: An additional tool for GPS quality control." Department of Mathematical Geodesy and Positioning. Thijsseweg, 1999.
  • Torge, Wolfgang. Geodesy Berlim, New York: Walter de Gruyter, 2001.
  • Unavco. "Teqc - Translation Edition Quality Control.", accessed April, 2015. http://www.unavco.org/software/data-processing/teqc/teqc.html
    » http://www.unavco.org/software/data-processing/teqc/teqc.html
  • Wang, Jinling. "Modelling and quality control for precise GPS and GLONASS satellite positioning" PhD diss, Curtin University of Technology, 1999.
  • Errata

    Na página 224, onde se lia:
    Tabela 1:
    Estimativa das precisões (em metros) e melhorias percentuais
    Leia-se:
    Tabela 1:
    Estimativa das precisões (em metros) e melhorias percentuais
    Boletim de Ciências Geodésicas (2016) 22(3):574

Datas de Publicação

  • Publicação nesta coleção
    Jun 2016

Histórico

  • Recebido
    Maio 2015
  • Aceito
    Out 2015
Universidade Federal do Paraná Centro Politécnico, Jardim das Américas, 81531-990 Curitiba - Paraná - Brasil, Tel./Fax: (55 41) 3361-3637 - Curitiba - PR - Brazil
E-mail: bcg_editor@ufpr.br