ANÁLISE ESTATÍSTICA DO POSICIONAMENTO RELATIVO A PARTIR DO AUSPOS

Statistical analysis of the relative positioning obtained from AUSPOS

Resumos

A disponibilidade de sistemas de processamento de dados GPS (Global Positioning System), online e gratuitos, tem aumentado nos últimos anos. Dentre esses destaca-se o AUSPOS (v.2.1) o qual emprega o programa BERNESE (v. 5.0) para realização do posicionamento relativo estático. O objetivo desse trabalho é verificar a qualidade do posicionamento planimétrico e altimétrico obtido com o uso do AUSPOS. Para isso foram utilizados arquivos contendo sessões de observações GPS de 4, 6 e 24 horas, no período de 22 dias consecutivos, no mês de setembro de 2013, das estações CRUZ, IMBT, IMPZ, MGIN e POAL, pertencentes à RBMC (Rede Brasileira de Monitoramento Contínuo dos Sistemas GNSS). Os resultados mostraram menor variabilidade posicional com o aumento da duração da sessão de rastreio e as acurácias planimétrica e altimétrica sub-centimétricas para o período de 24h e centimétricas para os períodos de 4h e 6h.

AUSPOS; Posicionamento Relativo GPS online; Acurácia

The availability of processing systems of GPS data, online and free, has increased in the last years. Among these, the AUSPOS(v.2.1) system, which employs the Bernese software (v.5.0) to realize of the static relative positioning, stands out. This paper aims at checking the quality of planimetric and altimetric positining obtained using the AUSPOS system. Files containing GPS observations of 4, 6 and 24 hour sessions were used in the period of 22 consecutive days in September 2013 of the CRUZ, IMBT, IMPZ, MGIN and POAL, stations belonging to RBMC (Brazilian Network for Continuous Monitoring of the GNSS systems). The results have shown lower positional variability with increase of the tracking session duration, and the planimetric and altimetric accuracies has reached subcentimeter level for 24 hours and centimetric level for 4 and 6 hours.

AUSPOS; GPS Relative Positioning; Accuracy

1. INTRODUÇÃO

Vários serviços online de processamento de dados GPS (Global Positioning System) têm sido desenvolvidos e disponibilizados gratuitamente por diferentes organizações (GHODDOUSI-FARD e DARE, 2006 ; OCALAN, ERGAN e TUNALIOGLU, 2013 ). Dentre esses serviços há os que permitem a realização do posicionamento relativo e Posicionamento por Ponto Preciso (PPP).As ferramentas gratuitas mais utilizadas que permitem realizar o PPP são: CSRS-PPP (Canadian Spatial Reference System - Precise Point Positioning), o IBGE-PPP (Instituto Brasileiro de Geografia e Estatística - Posicionamento por Ponto Preciso), o GAPS(GPS Analysis and Positioning Software) e o APPS (Automatic Precise Positioning Service); já as que permitem a realização do posicionamento relativosão: AUSPOS (Online GPS Processing Service), o OPUS (Online Positioning User Service) e o SCOUT (Scripps Coordinate Update Tool) ( DAWSON, GOVIND E MANNING, 2001 ; GHODDOUSI-FARD E DARE, 2006 ; ADIYANTO E ROBERTS, 2007 ; TSAKIRI, 2008 ; EL-MOWAFY, 2011 ; OCALAN, ERDGAN E TUNALIOGLU, 2013 ; GOMEZ, COGLIANO E TELLO,2013; IBGE, 2014 ).Os requisitos básicos para que usuários utilizem esses diferentes serviços são praticamente os mesmos: acesso à Internet e um endereço de e-mail válido ( GHODDOUSI-FARD E DARE, 2006 ). Para a realização dos experimentos deste trabalho foi utilizado o serviço de processamento de dados GPS - AUSPOS.

O método de posicionamento relativo estático é frequentemente utilizado na implantação de redes geodésicas, no transporte de coordenadas e na determinação de pontos temáticos. Como resultado do pós-processamento dos dados GPS são obtidos vetores linha de base os quais ligam as estações de coordenadas conhecidas àquelas a serem determinadas. Para a realização do pós-processamento o AUSPOS utiliza efemérides precisas do IGS (International GNSS Service), as quais estão vinculadas desde 07 de outubro de 2012, semana GPS 1709, até o presente, ao sistema de referência IGb08 ( REBISCHUNG, 2012 ; REBISCHUNG et al., 2012 ). Atualmente, todas as coordenadas determinadas pelo AUSPOS estão no GDA94 (Geocentric Datum of Australia 1994) e no ITRF2008, época da coleta dos dados.

Devido a adoção de diferentes realizações de sistemas de referência para determinação das velocidades e coordenadas das estações, e as mesmas estarem vinculadas a uma época específica, torna-se fundamental a compatibilização entre referenciais e épocas.

O objetivo deste trabalho é analisar a qualidade do posicionamento planimétrico e altimétrico obtido com uso do AUSPOS (v. 2.1) utilizando para isso, dados coletados por estações pertencentes à RBMC (Rede Brasileira de Monitoramento Contínuo dos Sistemas GNSS).

2. Sistemas de Referência e Realizações

O ITRS (International Terrestrial Reference System - Sistema de Referência Terrestre Internacional) é um sistema de referência tridimensional que gira com a Terra em seu movimento diurno no espaço. Nesse sistema, os pontos fixos à crosta da Terra sofrem pequenas variações com o tempo devido aos efeitos geofísicos ( PETIT E LUZUM, 2010 ).

O ITRS é realizado e mantido pelo IERS (International Earth Rotation and Reference Systems Service) ( ALTAMIMI, COLLILIEUX E MÉTIVIER, 2012 ). A realização mais atual é denominada de ITRF2008 e foi construída com base em observações obtidas por técnicas geodésicas espaciais: VLBI (Very Long Baseline Interferometry), SLR (Satellite Laser Ranging), DORIS (Doppler Orbitography and Radiopositioning Integrated by Satellite) e GNSS (Global Navigation Satellite System) ( ITRF, 2014 ). O ITRF2008 é composto por 934 estações distribuídas em 580 localidades, na qual 463 estão localizadas no hemisfério norte e 117 no hemisfério sul ( ALTAMIMI, COLLILIEUX E MÉTIVIER, 2012 ).

A realização do ITRS pelo IGS é denominada de IGb08 e foi obtida com o emprego de observações GNSS. As efemérides precisas disponibilizadas pelo IGS estão, atualmente, vinculadas ao referencial IGb08 o qual é alinhado ao ITRF2008 ( IGS, 2014 ).

O Sistema de Referência Geocêntrico para as Américas (SIRGAS) é realizado por meio de uma rede de estações GNSS de funcionamento contínuo denominada de SIRGAS-CON. As coordenadas semanais das estações são dadas no referencial IGb08, época média do período da coleta dos dados. Há também a solução multianual (SIR11P01) da rede SIRGAS-CON, na qual são fornecidas as coordenadas no ITRF2008, época 2005,0 e velocidades cartesianas geocêntricas no ITRF2005 das estações( SIRGAS,2014 ).

2.1 Relacionamento entre Sistemas de Referência

a)Mudança de Referencial das Coordenadas das Estações e das Velocidades

A transformação de similaridade é utilizada para relacionar dois sistemas de referência por meio de três translações, três rotações e um fator de escala (PETIT E LUZUM, 2010 ). Na atualização dos sete parâmetros de transformação e na mudança do referencial das velocidades dos pontos de interesse é necessário considerar suas variações temporais ( ITRF, 2014 ).

De acordo com Petit e Luzum (2010 ),Altamimi, Collilieux e Métivier (2012) e ITRF (2014) , o modelo matemático utilizado para transformação das coordenadas do sistema referência de origem (o) para o sistema de referência destino (d) pode ser expresso por:

em que:

: vetor das coordenadas cartesianas geocêntricas no sistema de referência destino, dadas em metros;

: vetor das coordenadas cartesianas geocêntricas no sistema de referência origem, dadas em metros;

[TX, TY, TZ] ^T : vetor das translações da origem do sistema de origem no sistema de referência destino, dadas em metros;

R _X ,R _Y ,R _Z : rotações diferenciais em torno dos eixos X, Y, Z respectivamente para estabelecer o paralelismo da rede origem com a rede destino, dadas em radianos;

D: fator diferença de escala, adimensional.

É importante salientar que os sete parâmetros de transformação e as coordenadas dos pontos nos dois sistemas devem estar na mesma época.

De acordo com Altamimi, Collilieux e Métivier (2012 ) o modelo matemático para realizar a mudança de referencial das velocidades das estações pode ser expresso por:

em que:

: vetor das velocidades cartesianas no sistema destino;

: vetor das velocidades cartesianas no sistema origem;

: vetor das variações temporais das translações;

: variações temporais das rotações diferenciais em torno dos eixos X, Y, Z;

D: variação temporal do fator de escala.

Na mudança de referencial das velocidades do ITRF2008 para ITRF2005, apenas a velocidade na componente X é alterada, pois com exceção da variação temporal da translação em torno do eixo X, as demais são iguais a zero.

Os parâmetros de transformação entre realizações do ITRS podem ser obtidos no site http://itrf.ensg.ign.fr/doc_ITRF/Transfo-ITRF2008_ITRFs.txt, os quais permitem o relacionamento do ITRF2008 para realizações anteriores. Deste modo, para transformações como, por exemplo, do ITRF2005 para ITRF2008, basta trocar o sinal do parâmetro. Para atualizar os parâmetros da época 2000,0 para outra época de interesse, pode-se utilizar o modelo matemático expresso por ( ITRF, 2014 ):

em que:

t: época de interesse dos parâmetros, dados em anos e sua fração;

t ₀ : época de origem dos parâmetros, dados em anos e sua fração.

b)Atualização de Coordenadas

A atualização de coordenadas de uma época origem para uma época de interesse pode ser realizada com o emprego das equações SIRGAS (2014) :

em que:

t: época de interesse para a determinação das coordenadas cartesianas, dada em ano e sua fração;

t ₀ : época de origem das coordenadas cartesianas, dada em ano e sua fração;

X(t),Y(t),Z(t):coordenadas cartesianas geocêntricas de uma estação na época de destino;

X(t ₀ ),Y(t ₀ ),Z(t ₀ ):coordenadas cartesianas geocêntricas de uma estação na época de origem.

Este procedimento permite que seja considerado o deslocamento da estação no intervalo de tempo, dado entre a época de determinação da coordenada e a época de interesse, devido ao movimento da placa litosférica.

3. Serviço online de processamento de dados GPS - AUSPOS

AUSPOS é um serviço de pós-processamento de dados GPS, gratuito, fornecido pela agência australiana de geociências que pode ser acessado em < http://www.ga.gov.au/earth-monitoring/geodesy/auspos-online-gps-processing-service.html>. A versão 2.1 do AUSPOS foi lançada em primeiro de outubro de 2013.

A determinação posicional no ITRF em qualquer parte da Terra pelo AUSPOS é feita utilizando produtos IGS tais como: efemérides precisas, parâmetros de orientação da Terra, coordenadas e observáveis GPS das estações ativas da rede IGS ( GEOSCIENCE AUSTRALIA, 2014 ).

O pós-processamento dos dados GPS é feito utilizando o método relativo estático com o emprego do software Bernese GNSS, desenvolvido pela AIUB (Astronomical Institute of the University of Bern).Informações sobre esse software podem ser obtidas em (http://www.bernese.unibe.ch/ ).

Atualmente, para determinação das coordenadas no referencial ITRF2008, época média dos dados coletados, deve-se enviar, via upload ou ftp (File Transfer Protocol, as observáveis,no formato RINEX (Receiver INdependent EXchange), coletadas no modo estático, por receptores geodésicos de dupla frequência. Quando o arquivo de observações RINEX é submetido, 15 estações ativas da rede IGS e APREF (Asia-Pacific Reference Frame) mais próximas ao ponto a ser determinado poderão ser utilizadas como estações de referência para o posicionamento relativo e ajustamento ( GEOSCIENCE AUSTRALIA, 2014 ).

Nos processamentos, os efeitos das fontes de erros, tais como: erro do relógio dos receptores, troposfera e ionosfera são levados em consideração por meio da modelagem ou estimação dos parâmetros relacionados. Todos os cálculos são efetuados de acordo com as convenções do IERS< http://www.iers.org/ IERS/EN/DataProducts/Conventions/conventions.html?__nnn=true >. Salienta-se que a correção do efeito devido a carga oceânica não é levada em consideração nos processamentos.

Para a realização do processamento o usuário deverá submeter ao serviço apenas os dados de observação GPS, o tipo e altura da antena utilizada e um e-mail válido para o recebimento do relatório de processamento. Embora as observações da RBMC sejam armazenadas com intervalo de 15 segundos, são empregados intervalos de 30 segundos para o processamento.

OAUSPOS utiliza a nomenclatura de antena adotada pelo IGS e parâmetros de calibração NGS (National Geodetic Survey). Podem ser submetidos arquivos de observação de até 20 estações. No relatório estão presentes as coordenadas no GDA94 e no ITRF2008.

4. Metodologia

4.1 Obtenção dos dados

Os dados GPS das estações da RBMC utilizados compreendem o período de 08/09/2013 a 29/09/2013 e correspondem aos dias do ano de 251 a 272, perfazendo um total de 22 dias consecutivos. Foram escolhidas cinco estações da RBMC que contribuíram para solução multianual da rede SIRGAS-CON. As estações utilizadas são denominadas de: CRUZ (Cruzeiro do Sul / AC); IMBT (Imbituba / SC); IMPZ (Imperatriz / MA); MGIN (Inconfidentes / MG) e POAL (Porto Alegre / RS). As observações GPS dessas estações, no formato RINEX, foram adquiridas no site do IBGE mantenedor do projeto RBMC. Os períodos de observações utilizadas são da 0h à 4h, das 0h às 06h e das 0h às 24h.

4.2 Determinação das Coordenadas pelo AUSPOS

No processamento dos dados GPS, o AUSPOS utilizou 15 estações de controle pertencentes à rede IGS e efemérides precisas finais. No ajustamento, as coordenadas das estações IGS foram injuncionadas com incertezas de 0,001m para as componentes horizontais e de 0,002m para vertical. Todas as coordenadas foram determinadas em ITRF2008, época 2013,7.

No Quadro 1 encontram-se os principais modelos e estratégias utilizadas para o processamento das linhas de base pelo software Bernese no serviço AUSPOS.

Quadro 1:
Estratégias e modelos utilizadas no processamento pelo AUSPOS versão 2.1.

A Figura 1 ilustra a distribuição espacial aproximada das estações IGS e da RBMC utilizadas no processamento. Os triângulos representam as estações de referência da rede do IGS e os círculos representam as estações da RBMC utilizadas para essa pesquisa.

Figura 1:
Distribuição espacial aproximada das estações utilizadas nos processamentos. Fonte: Adaptado de ( GEOSCIENCE AUSTRALIA, 2014 ).

Estações de referência do IGS utilizadas: BOGT, BRFT, BRMU, CHPI, CRO1, GLPS, GODE, GUAT, ISPA, MANA, OHI2, PALM, PARC, UNSA e USNO.

4.3 Determinação das Coordenadas de Referência

Para verificar a acurácia do posicionamento proporcionado pelo AUSPOS, as coordenadas obtidas para cada estação foram comparadas com suas posições, provenientes da solução multianual (SIR11P01) da rede SIRGAS-CON, e tomadas como referência. As coordenadas das estações dessa solução são vinculadas ao ITRF2008, época 2005,0. A partir das velocidades das estações no ITRF2008foi possível atualizar suas coordenadas para a época média do período dos dados coletados (2013,7).Foi adotada a época média dos dados coletados para fins de comparação das coordenadas, pois em 11 dias, os deslocamentos das estações, devido ao movimento da placa litosférica,foram menoresque0,001 m.

No Quadro 2 encontram-se as velocidades cartesianas geocêntricas das estações obtidas da solução SIRG11P01, as quais podem ser encontradas no site ( http://www.sirgas.org/fileadmin/docs/SIR11P01.VEL ).

Quadro 2 :
velocidades cartesianas geocêntricas no ITRF2005.

Como as velocidades das estações são dadas no ITRF2005 e as coordenadas de origem estão no ITRF2008, foi necessário realizar a mudança de referencial das velocidades do ITRF2005 pra ITRF2008. Para isso, mudou-se o referencial das coordenadas das estações no ITRF2008 (2005,0) para ITRF2005 (2005,0) a partir das equações (1) e (3) com o uso dos parâmetros de transformação disponíveis no site do ITRF. Para mudança do referencial das velocidades do ITRF2005 para ITRF2008 foi utilizado o modelo matemático (2). A partir da compatibilidade dos referenciais das velocidades e das coordenadas foi possível por meio das equações (4), (5) e (6) atualizar as coordenadas ITRF2008, da época 2005,0 para a época média 2013,7.

As análises foram efetuadas utilizando as coordenadas planas no sistema Universal Transversa de Mercator (UTM) e as altitudes geométricas (elipsoide GRS80). Para isso, todas as coordenadas cartesianas geocêntricas em ITRF2008 (2013,7) foram convertidas para coordenadas planas UTM com o uso do aplicativo ProGriD, o qual pode ser obtido no site: ftp://geoftp.ibge.gov.br/aplicativos/ transformacao_coordenadas/progrid.zip

Na Tabela 1 encontram-se as coordenadas de referência,no sistema UTM,das estações utilizadas neste trabalho.

Tabela 1 :
UTM, ITRF2008 (2013,7), das estações.

4.4 Análises Estatísticas

Para a construção do Intervalo de Confiança (IC) para a média populacional μ, desconhecendo a variância populacional ^2 e para n<30, a variável z da distribuição normal é substituída pela variável t da distribuição de t Student. Para sua aplicação, foi necessário verificar o pressuposto de normalidade da população com o uso do teste não paramétrico de Shapiro-Wilk, pois o tamanho da amostra é menor que 50.As hipóteses foram enunciadas da seguinte forma:

H0: a amostra provém de uma população normal.

Ha: a amostra provém de uma população não normal.

Foi utilizado nível de significância = 5% e a decisão do teste foi tomada com base na probabilidade de significância (p-valor). Esse pode ser interpretado como o grau de concordância entre os dados e a hipótese nula. Deste modo, quanto menor for o p-valor, menor será a consistência entre os dados e a hipótese nula( MONTGOMERY E RUNGER, 2009 ). Conhecendo-se o p-valor pode-se testar a hipótese, sem impor, formalmente, um nível pré-selecionado de significância.

A decisão do teste de hipótese é tomada comparando-se o p-valor com o nível de significância fixado.

Se p-valor ≤ → rejeita-se H0

Se p-valor > → não se rejeita H0

Como foi selecionado = 5% se o p-valor for menor que 0,05 rejeita-se a hipótese nula.

A partir das coordenadas de referência e as obtidas via AUSPOS foram verificadas as acurácias, dado pelo Erro Quadrático Médio das componentes E(EQM _E)e N(EQM _N ); o Erro Planimétrico (EP), e as acurácias planimétrica e altimétrica (EQM _p e EQM _h ).

A acurácia das componentes planimétricas (E, N) é expressa matematicamente por (FGDC, 1998 ):

em que:

EAUSPOS,i,NAUSPOS,i: são as coordenadas da estação obtidas pelo AUSPOS;

ERef.,i, NRef.,i: são as coordenadas de referência da estação;

n: número de coordenadas;

i: é um inteiro que varia de 1 a n.

O Erro Planimétrico (EP) para uma estação i é definido como (FGDC, 1998):

A acurácia planimétrica (EQM _p ) é dada por ( FGDC, 1998 ):

ou

O Erro Altimétrico (EA) no ponto i dado por:

A acurácia altimétrica é dada por ( FGDC, 1998 ):

em que:

h _AUPOS,i , h _AUSPOS,i : é a altitude geométrica da estação obtida pelo AUSPOS;

h _Ref.,i , h _Ref.,i : é a altitude geométrica de referência da estação.

A média e variância amostral para cada componente posicional pode ser calculada por( MONTGOMERY E RUNGER, 2009 ):

em que:

n: número de dias de coleta de dados;

X _i : componentes (E, N, h) para cada dia, obtidas pelo AUSPOS, na época média adotada;

: média para cada componente (E, N, h);

: variância das componentes (E, N, h).

Outra forma de avaliar a acurácia proposta por Gauss, denominada EQM (Erro Quadrático Médio), a qual permite quantificar os valores de tendência e precisão, pode ser representada matematicamente pela expressão ( MONICO et al., 2009 ); MIKHAIL E ACKERMANN, 1976 ):

em que:

b: viés (tendência) (média das discrepâncias entre as coordenadas obtidas diariamente pelo AUSPOS em relação às de referência);

: variância amostral das discrepâncias;

De acordo com Monico et al. (2009 ), a expressão (16) para grandes amostras é praticamente igual à média quadrática dos erros .

As discrepâncias obtidas a partir da diferença entre os valores estimados com o uso do AUSPOS e os de referências são dadas por:

em que:

E _i ,N _i ,h _i : componente posicional de cada estação obtida pelo AUSPOS para cada dia;

E _R ,N _R ,h _R : coordenadas de referência para cada estação.

: discrepâncias das componentes (E, N, h).

O valor da tendência para cada componente de uma estação foi calculado a partir da média das discrepâncias ( ) obtidas em 22 dias. Sua precisão foi obtida a partir do desvio padrão da média das discrepâncias.

Para verificar se a tendência é significativa ou não foi construído um Intervalo de Confiança (IC) de 95% para a média populacional , no qual foi verificado se o intervalo construído continha a coordenada de referência.

A expressão matemática utilizada para o IC para a média populacional "" _"X" de uma população normalmente distribuída pode ser dada por ( MONTGOMERY E RUNGER, 2009 ):

em que:

: média das componentes (E, N, h) para cada estação analisada;

: média populacional para cada componente;

n: número de dias de coleta de dados;

: valor teórico da estatística t para com (n-1) graus de liberdade.

Se o valor de referência estiver inserido nesse intervalo, considera-se que não houve tendência ao nível de probabilidade "1- ", ou seja, o IC contém a verdadeira média com 95% de confiança ( MONICO et al., 2009 ; MONTGOMERY E RUNGER, 2009 ).

5. Resultados e Discussão

Para a determinação das estatísticas e realização dos testes de hipóteses foi empregado o software livre R v. 3.1.0.

A Tabela 2 contém as porcentagens médias das ambiguidades resolvidas (% MAR) e os comprimentos médios das linhas de base (CMLB) envolvidos na determinação posicional das estações CRUZ, IMBT, IMPZ, MGIN e POAL para cada dia utilizado neste trabalho.

Tabela 2 :
Ambiguidades (% MAR) e comprimentos das linhas de base(CMLB).

Verifica-se na Tabela 2 que a % MAR para duração da sessão de observação de 4 h tem valor máximo de 59,1, mínimo 21,7 e médio de 40,9. Para o CMLB foram encontrados valor máximo de 1.821,4 Km, mínimo de 1.673,2 Km e médio de 1.760,4 Km. A % MAR para duração da sessão de observação de 6 h tem valor máximo de 66,8, mínimo 35,0 e médio de 53,8. Para o CMLB foram encontrados valor máximo de 1.848,2 Km, mínimo de 1.696,6 Km e médio de 1.790,4 Km. A % MAR para duração da sessão de observação de 24 h valor máximo de 83,1, mínimo 72,6 e médio de 78,3. Para o CMLB foram encontrados valor máximo de 1.839,0 Km, mínimo de 1.640,5 Km e médio de 1.779,7 Km.

De acordo com os resultas apresentados na Tabela 2 , verifica-se que o aumento da duração da sessão de observação implica na ampliação da porcentagem média de resolução das ambiguidades das linhas de base. Assim, espera-se resultados mais confiáveis na determinação posicional das estações.

A Tabela 3 contém as coordenadas planas (E, N) no sistema UTM, fuso 18, Hemisfério Sul, e altitude geométrica (h) para a estação CRUZ no ITRF2008, época 2013,7.

Tabela 3 :
Coordenadas planas e altitude geométrica.

A Tabela 4 apresenta as estatísticas descritivas para as amostras das componentes (E, N, h) no ITRF2008, época 2013,7.

Tabela 4
: Estatísticas para as coordenadas (E, N, h) obtidas nas diferentes estações para sessões de observação de 4h, 6h e 24h.

Em que: Média (M), Amplitude Total (AT) e Desvio Padrão (DP) em metros; Coeficiente de Variação (CV) em porcentagem. Os valores das estatísticas foram calculados considerando os valores obtidos para cada componente posicional no período de 22 dias.

A partir das AT disponíveis na Tabela 4 , considerando as cinco estações podem-se obter as variabilidades médias das altitudes geométricas (h), as quais são: 0,055 m para 4h, 0,042 m para 6h e 0,019 m para 24h. Desse modo, verifica-se que aumentando a duração de coleta de dados de 4h para 6h e 24h houve diminuição de 0,013 m para 6h e 0,036 m para 24h da variabilidade média dada pela AT.

Observa-se a partir do desvio padrão médio das altitudes geométricas, diminuição da variabilidade média de 0,003 m para 6h e 0,009 m para 24h em relação a 4h de coleta de dados. Em relação aos CV, nota-se que na média, não houve uma melhora na mudança de 4h para 6h de coleta de dados, mas detecta-se diminuição de 0,014% para 0,006% quando da utilização de 24 horas de dados.

As AT médias encontradas para a componente planimétrica E foram 0,046 m para 4h, 0,020 m para 6h e 0,008 para 24h. Desse modo, observa-se uma diminuição da AT média de 0,026 m para 6h e de 0,038 m para 24h em relação a 4h de dados coletados. Em relação aos desvios padrão médio, observam-se valores de 0,010 m para 4h, 0,005 m para 6h e 0,002 m para 24h.

As AT médias encontradas para a componente planimétrica N foram 0,018 m para 4h, 0,017 m para 6h e 0,006 m para 24h. Observa-se que a AT média diminui de 0,001 m com a mudança de duração da sessão de observação de 4h para 6h e 0,012 m para 24h. O desvio padrão médio foi de 0,004 m para 4h e 6h e de 0,002 m para 24h. Nota-se a partir dos resultados, que na média, a variabilidade da componente E é maior do que a componente N.

A Tabela 5 apresenta os resultados do teste de normalidade com o emprego de Shapiro-Wilk.

Tabela 5:
P-valor para cada componente posicional.

*Significativo para = 0,05: população não segue uma distribuição normal.

Pode-se verificar que nem todas as componentes (E, N, h) das estações são provenientes de uma distribuição normal a 5% de significância com a aplicação do teste de Shapiro-Wilk. Porém, a 1% de significância aumenta o número de componentes advindas de população normal.

Na Tabela 6 são apresentados os valores da acurácia planimétrica dada pelo (EQM _P ), o Erro Planimétrico Médio (EPM), o desvio padrão do erro planimétrico (DPEP), o erro altimétrico médio (EAM) e o desvio padrão do erro altimétrico (DPEA).

Tabela 6 :
Valores, em metros, de acurácia e medidas de precisão.

Comparando os valores médios das acurácias planimétricas (EQM _P ) para as cinco estações, verifica-se que aumentando o período de coleta dos dados de 4h para 6h e 24h há uma melhora na qualidade do posicionamento de 28,36% (0,004 m) e 47,76% (0,006 m), respectivamente. Nota-se melhora 27,08% (0,003 m) quando da utilização de 24h de dados em relação à 6h.

Os erros planimétricos máximo e mínimo encontrados foram, respectivamente, 0,015 m e 0,008 m para 4 h, 0,011 m e 0,006 m para 6 h e de 0,010 m e 0,004 m para 24 h. Para 4h de dados, o erro planimétrico médio foi de 0,011 m, 0,008 m para 6h e 0,006 m para 24h.

A Figura 2 ilustra os valores de acurácia planimétrica e altimétrica obtidas para os tempos de coletas de 4, 6 e 24 horas, para as estações CRUZ, IMBT, IMPZ, MGIN e POAL.

A partir da Figura 2 , nota-se uma melhora na acurácia planimétrica em todas as estações com o aumento da duração da coleta dos dados. O mesmo não ocorreu com a acurácia altimétrica.

As acurácias planimétricas máxima e mínima foram, respectivamente, 0,016 m e 0,010 m para 4h, 0,012 m e 0,007 m para 6h e 0,011 m e 0,005 m para 24h. A acurácia planimétrica média para 4 h de coleta de dados foi de 0,013 m, 0,010 m para 6h e de0,007 m para 24h, considerando as acurácias planimétricas média de cada estação.

As acurácias altimétricas máximas e mínimas para 4h foram, respectivamente, 0,019 m e 0,011m; 0,014 m e 0,009 m para 6 h e de0,014 m e 0,005 m para 24h. A acurácia altimétrica média para 4 h de coleta de dados foi de 0,014 m, para 6h foi de 0,012 m e para 24h foi de 0,008 m.

Figura 2:
Acurácias planimétrica e altimétrica

A Figura 3 ilustra os erros planimétricos médios e os desvios padrão dos erros planimétricos em cada estação e duração da coleta dos dados.

Figura 3 :
Erro planimétrico médio e desvio padrão do erro planimétrico.

Percebe-se diminuição da variabilidade do erro planimétrico com o aumento do tempo de coleta dos dados.

Os Intervalos de Confiança de 95%, em metros, para as médias populacionais das componentes (E, N, h) das cinco estaçõespara coleta de dados de 4h, 6h e 24h e os resultados da análise de tendência estão apresentados na Tabela 7 .

Tabela 7 :
IC de 95% para as médias populacionais

A tendência pôde ser avaliada verificando se os valores das componentes posicionais tomadas como referência, disponíveis na Tabela 1 , estavam contidos em cada intervalo de confiança de 95% para a média populacional.

Pode-se verificar na Tabela 7 que aproximadamente 13,3 % das componentes (E, N, h) apresentaram tendência não significativa para = 0,05 para 24h de coleta de dados.

Nos casos onde as tendências não foram significativas para = 0,05 de significância, os valores de precisão se igualam ao valores calculados para as acurácias. Para as componentes que apresentaram tendência significativa o valor da REQM representa, de fato, a acurácia posicional, na qual envolve não apenas efeitos aleatórios, mas também efeitos sistemáticos dados pelas tendências.

Na Tabela 8 sãoapresentadas duas outras formas de representar os valores da acurácia, dados por (m) e EQM obtidos para as componentes (E, N, h), com o uso da equação (16) (com extração da raiz quadrada) e com as equações (7), (8) e (13).

Tabela 8:
Medidas de acurácia, em metros, para cada componente posicional das estações sob diferentes tempos de coleta de dados.

A partir dos valores presentes na Tabela 8 , percebe-se que na medida que aumenta a duração da coleta dos dados, os valores diminuem, isto é, aumenta a acurácia no posicionamento.

6. Considerações Finais e Conclusões

Para que o AUSPOS seja utilizado para posicionamento no Brasil deve-se realizar a transformação de referencial das coordenadas obtidas no ITRF2008 (IGb08) para o SIRGAS2000, o qual é o sistema de referência adotado pelo SGB. Recomenda-se determinar as velocidades das estações com o uso do modelo de velocidade VEMOS2009, adotado pelo IBGE, para atualizar das coordenadas da época do rastreio para a época 2000,4.

Para aplicação do teste t de Student é necessário que as observações que compõem a amostra sejam independentes e provenientes de uma população com distribuição gaussiana (normal). Foi verificado que nem todas as populações são gaussianas com a utilização do teste não paramétrico de Shapiro-Wilk(para n≤50).Isso passa a ser uma limitação na aplicação dessa ferramenta.

Pôde-se verificar que o aumento do tempo de coleta proporciona melhora na medida de acurácia das componentes do vetor posição, o que reflete diretamente na qualidade da determinação posicional.

Os resultados indicam que o posicionamento com o uso do AUSPOS possibilita determinação de acurácias planimétrica e altimétrica centimétrica para 4h e 6h de coleta de dados e sub-centimétrica para 24h. Assim, verifica-se que esse serviço poderá ser utilizado para densificação de redes geodésicas e determinação de pontos de controle para execução de levantamentos georreferenciados.

Recomenda-se a realização de novos experimentos considerando uma amostra com tamanho maior que a utilizada (22). E verificação da normalidade e independência das discrepâncias posicionais das componentes (E, N, h) de modo que seja possível a determinação das acurácias planimétricas e altimétricas com nível de confiança de 95%, conforme especificações encontradas em FGDC (1998) .

• ADIYANTO, F.H., ROBERTS, C.A., 2007. Temporary CORS networks for land reconstruction in Aceh, Sumatera.IGNSS2007 Symp. on GPS/GNSS, Sydney, Australia, 4-6 December, paper 54, CD-ROM procs. Disponível em <http://www.gmat.unsw.edu.au/snap/about/publications_author.htm> Acesso em: 22 jan. 2014.
  » http://www.gmat.unsw.edu.au/snap/about/publications_author.htm
• ALTAMIMI, Z.; COLLILIEUX, X.; MÉTIVIER, L. Analysis and results of ITRF2008 - IERS Technical Note No. 37. IERS. Frankfurt am Main: Verlag des Bundesamts für Kartographie und Geodäsie, p. 54. 2012.
• AUSPOS - Online GPS Processing Service. Australia Geoscience: Applying geoscience to Australia's most important challenges, 2014. Disponivel em: <http://www.ga.gov.au/earth-monitoring/geodesy/auspos-online-gps-processing-service.html>. Acesso em: 22 jan. 2014.
  » http://www.ga.gov.au/earth-monitoring/geodesy/auspos-online-gps-processing-service.html
• DAWSON, J., GOVIND, R., MANNING, J. Application of the AUSLIG Online GPS Processing System (AUSPOS) to Antarctica, Proceeding of Satnav2001, Canberra, July 2001. Disponível em <http://www.ga.gov.au/image_cache/ GA5057.pdf>. Acesso em: 22 jan. 2014.
  » http://www.ga.gov.au/image_cache/ GA5057.pdf
• EL-MOWAFY, A. Analysis of web-based gnss post-processing services for static and kinematic positioning using short data spans. Survey Review, 43, October 2011. 535-549. Disponivel em: <http://www.maneyonline.com/loi/sre>. Acesso em: 22 jan. 2014.
  » http://www.maneyonline.com/loi/sre
• FEDERAL GEOGRAPHIC DATA COMMITTEE - FGDC. Geospatial Positioning Accuracy Standards Part 3: National Standard for Spatial Data Accuracy. United States. 1998.
• GHODDOUSI-FARD, R.; DARE, P. Online GPS processing services: an initial study. GPS Solutions, 0, 2006. 10 - 20. Disponivel em: <http://folk.uio.no/treiken/GEO4530/Online_GPS_posisitoning_services.pdf>. Acesso em: 22 jan. 2014.
  » http://folk.uio.no/treiken/GEO4530/Online_GPS_posisitoning_services.pdf
• GEOSCIENCE AUSTRALIA. Australian Government -<http://www.ga.gov.au/earth-monitoring/geodesy/auspos-online-gps-processing-service.html>. Acesso em: 22 jan. 2014.
  » http://www.ga.gov.au/earth-monitoring/geodesy/auspos-online-gps-processing-service.html
• GOMEZ, M. E.; COGLIANO, D.; TELLO, M. N. R. Análisis de la performance de tres servicios de posicionamiento GPS online en Argentina. Geoacta, Argentina, 38, n. 2, 2013. 95-106. Disponivel em: <http://ppct.caicyt.gov.ar/ index.php/geoacta/index> Acesso em: 22 jan. 2014.
  » http://ppct.caicyt.gov.ar/ index.php/geoacta/index
• IBGE. Instituto Brasileiro de Geografia e Estatística, 2014. Disponivel em: <http://www.ibge.gov.br/home/default.php>. Acesso em: 22 jan. 2014.
  » http://www.ibge.gov.br/home/default.php
• IERS. International Earth Rotation and Reference Systems Service, 2014. Disponivel em: <http://www.iers.org/IERS/EN/IERSHome/home.html ?__nnn=true>. Acesso em: 22 jan. 2014.
  » http://www.iers.org/IERS/EN/IERSHome/home.html ?__nnn=true
• IGS. International GNSS Service, 2014. Disponível em <http://www.igs.org/>. Acesso em: 15 fev. 2014.
  » http://www.igs.org/
• ITRF. International Terrestrial Referece Frame, 2014. Disponivel em: <http://itrf.ensg.ign.fr/>. Acesso em: 22 jan. 2014.
  » http://itrf.ensg.ign.fr/
• MIKHAIL, E. M.; ACKERMANN, F. Observation and Least Squares. New York IEP, Dun-Donelly Publisher, 1976.
• MONICO, J.F.M.; PÓZ, A.P.D.; GALO, M.; SANTOS, M.C.; OLIVEIRA, L.C. Acúrácia e Precisão: Revendo os conceitos de forma acurada. Bol. Ciênc. Geod., séc.Comunicações, Curitiba, vol. 15, nº 3, p.469-483, jul-set, 2009.
• MONTGOMERY, D. C.; RUNGER, G. C. Estatística aplicada e probabilidade para engenheiros. Tradução de Verônica Caldo. 4. ed. Rio de Janeiro: LTC, 2009. 490 p.
• OCALAN, T.; ERDGAN,; TUNALIOGLU, N. Analysis of web-based online services for gps relative and precise point positioning techniques. Boletim de Ciências Geodésicas, Curitiba, 19, n. 2, abr-jun 2013. 191-207. Disponivel em: <http://ojs.c3sl.ufpr.br/ojs2/index.php/bcg>. Acesso em: 22 jan. 2014.
  » http://ojs.c3sl.ufpr.br/ojs2/index.php/bcg
• PETIT, G.; LUZUM, B. IERS Conventions (2010) - IERS Technical Note 36. IERS. Frankfurt am Main: Verlag des Bundesamts für Kartographie und Geodäsie, p. 179. 2010.
• R Development Core Team (2009). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0, URL <http://www.R-project.org>
  » http://www.R-project.org
• REBISCHUNG, P. et al. IGS08: the IGS realization of ITRF2008., GPS Solutions 16, 2012. 483-494.
• REBISCHUNG, P. IGb08: an update on IGS08. Disponivel em: <http://igscb.jpl.nasa.gov/pipermail/igsmail/2012/007853.html>. Acesso em: 22 fev. 2014.
  » http://igscb.jpl.nasa.gov/pipermail/igsmail/2012/007853.html
• SIRGAS. Sistema de Referência Geocêntrico para as Américas, 2014. Disponivel em: < http://www.sirgas.org/index.php?id=15> Acesso em: 22 jan. 2014.
  » http://www.sirgas.org/index.php?id=15
• TSAKIRI, M. GPS Processing using online services. Journal of Surveying Engineering, 134, n. 4, November 2008. 115-125. Disponivel em: <http://ascelibrary.org/doi/abs/10.1061/%28ASCE%290733-9453%282008%29134%3A4%28115%29>. Acesso em: 22 jan. 2014
  » http://ascelibrary.org/doi/abs/10.1061/%28ASCE%290733-9453%282008%29134%3A4%28115%29

Datas de Publicação

• Publicação nesta coleção
  Mar 2015

Histórico

• Recebido
  Jun 2014
• Aceito
  Out 2014