Importação das Redes Geodésicas do IBGE

From OpenStreetMap Wiki
Jump to navigation Jump to search

(EM DESENVOLVIMENTO)


Descrição e Procedimento Gerais

Descrição

"Proposta de Importação dos pontos (nós) das Redes Geodésicas do IBGE".
PropósitoːTer no OSM pontos com coordenadas e altitude precisas, para aferição de imagens, limites de fronteiras, localidades não mapeadas, etc.


Situação atual - histórico

Discussão inicialː no canal do Telegram em 06/03/2017ː https://t.me/OSMBrasil_Comunidade/18192

Publicado no talk-br em: https://lists.openstreetmap.org/pipermail/talk-br/2017-March/012063.html

Inicialmente importado emː http://www.openstreetmap.org/changeset/46795156 - 27/02/2017

Problema técnico apontado pelo DWG no changesetː "relação com excesso de membros, 37.948 objetos (nós), causou problema em osm2pgsql".

Removido emː
http://www.openstreetmap.org/changeset/46798238
http://www.openstreetmap.org/changeset/46798106
http://www.openstreetmap.org/changeset/46797916
http://www.openstreetmap.org/changeset/46797703
http://www.openstreetmap.org/changeset/46797541
http://www.openstreetmap.org/changeset/46797326
http://www.openstreetmap.org/changeset/46796572


Atualmenteː
Recuperado pacote dos objetos originais (re-reversão), salvo, corrigido ("removida a relação"ː sem mais relação agrupando os nós), feita adição de dados de altitude por conversão cf. IBGE para os nós que faltavam estes dados, e armazenado para eventual reposição.

O que ainda faltaː
Fazer "merge" entre os nós "antigos" recuperados (com histórico original) e os nós "novos" com dados complementados para as mesmas posições.

(Atualmente em desenvolvimento porː http://www.openstreetmap.org/user/smaprs. Com o valioso apoio de discussões prévias com a comunidade no Telegram, e mais recentes na talk-br.)


Questões iniciais sobre manutenção dos dados e eventuais problemas

Quando necessário os nós podem (teoricamente) ser alterados em bloco, ou se pode fazer uma reposição geral: isto é, quando houver significativa alteração de coordenadas.

As tags "description" são propostas para identificar os grupos e poder recuperar os objetos quando necessário, conforme os 6 tipos de estação.

Possibilidade de fazer "merge" paraː

  • substituir nós com alteração de coordenadas;
  • adicionar ou trocar dados de cada estação nas mesmas coordenadas ou alteradas;

PROBLEMAː o validator do JOSM consegue detectar todos os nós a serem subsituídos/alterados, como sobrepostos ( "nodes at same position" );
porém, somente executa "merge" (para fundir dados em um só nó), 01 nó de cada vez, o que atualmente torna inviável uma evetual manutenção/substituição de dados no pacote todo.
Solução possívelː fazer merge através de script próprio.

Notaː Considerando as medições para a "deriva continental" ("Continental Drift") no Brasil, a velocidade de movimentação das placas litosféricas é de cerca de 2cm/ano (http://www.ibge.gov.br/home/geociencias/geodesia/default_sirgas_int.shtm?c=11). Assim, daqui uns 25 anos talvez a alteração nas coordenadas atinja o deslocamento de 1 pixel, considerando as imagens atuais de no máximo 50cm/pixel, quando pode ser necessária uma revisão geral nos pontos conforme novos dados do IBGE.


Fonte dos dados: IBGE

Fonte: IBGE (acesso via Visualizador do INDEː http://www.visualizador.inde.gov.br/ - Menu "Tema"ː Redes Geodésicas)

Redes Geodésicas - originaisː Arquivos SHP baixados já em WGS84


Rede Geodésica Planimétrica:

  • Rede Brasileira de Monitoramento Contínuo - RBMC/GNSS Permanente: RBMCPoint.SHP / 4 KB / OBJECTS: 129
  • Estações SAT DOPPLER: BDG_DOPPoint.SHP / 28 KB / OBJECTS: 1006
  • Estações SAT GPS: BDG_GPSPoint.SHP / 82 KB / OBJECTS: 2977
  • Estações de Poligonal: BDG_EPPoint.SHP / 32 KB / OBJECTS: 1133
  • Vértices de Triangulação: BDG_VTPoint.SHP / 100 KB / OBJECTS: 3643

Sub-Total Planimétrica: 7882 OBJECTS


Rede Geodésica Altimétrica:

  • Referências de Nível: BDG_RNPoint.SHP / 613 MB / OBJECTS: 69887

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Total Original (Rede Planimétrica + Altimétrica): 77769 OBJECTS

TOTAL SELECIONADO A SER IMPORTADO = 45622 OBJECTS (somente as estações com campo "SITUACAO=BOM")



Redes Geodésicas-INDE.jpg


Sobre as referências (CRS) dos dados do IBGE

1) Exemplo conforme o RBMC-IBGE "Relatório de Informação de Estação"ː Estação POAL - Porto Alegreː ftp://geoftp.ibge.gov.br/informacoes_sobre_posicionamento_geodesico/rbmc/relatorio/Descritivo_POAL.pdf

"3. Coordenadas oficiais / 3.1. SIRGAS2000 (Época 2000.4)ː Coordenadas Geodésicas...Alt. Elip.: 76,745 m Sigma: 0,002 m"

2) Conforme SHP via INDE do mesmo ponto (RBMC-POAL): ALTGEOM=76,745 / DATUM_ALT=imbituba

3) Conforme IBGEː TRANSFORMAÇÃO DE COORDENADASː (http://www.ibge.gov.br/home/geociencias/geodesia/pmrg/faq.shtm#11)

"Atualmente não existem parâmetros de transformação entre SIRGAS2000 e WGS 84 porque eles são praticamente iguais, ou seja, DX = 0, DY = 0 e DZ = 0."

Portanto:

  • 1) Altitude Elipsoidal (e demais coordenadas) está em SIRGAS2000;
  • 2) Altitude Elipsoidal = ALTGEOM (DATUM_ALT=imbituba) - campo do SHP
  • 3) ALTGEOM: conversão direta para ele = * em WGS84 (DZ=0)

Redes Geodésicas-QGIS.jpg


Outros documentos importantes do IBGEː

Banco de Dados Geodésicos Onlineː

http://www.ibge.gov.br/home/geociencias/geodesia/sgb.shtm

Tutorial do Banco de Dados Geodésicosː

ftp://geoftp.ibge.gov.br/metodos_e_outros_documentos_de_referencia/outros_documentos_tecnicos/bdg/bdg_tutorial_fev_2015.pdf


Procedimento gerais de conversão para todos SHP originais

1) Conversão CRS: já baixados convertidos em WGS84 nos originais;


2) Seleção dos objetos: somente aqueles com classificação de qualidade no campo "SITUACAO"=BOM ;


3) Seleção e conversão de CAMPOS para TAGS úteis para o OSM:


3.1) QGISː Conversão das tags dos DADOS ORIGINAISː Identificação, Localidade, Altitude, URL, etc:

  • ref=* (identificar o campo compatível, como: nome da estação, etc)
  • ref:loc=* (identificar o campo compatível com referência da localidade)
  • ele=* (identificar os campo compatíveis com altitude em WGS84 a converter - diretamente ou via Altura Geoidal "N" ; alterar separação decimal e de milharː suprimir "." ponto de milhar; trocar separação decimal de "," para "." e vice-versa)
  • description=IBGE-<TEMAː RBMC/GNSS; SAT DOPPLER; SAT GPS; EPoligonal; VTriangulação; Rede Altimétrica-RN>


3.2) JOSMː Adição das TAGS próprias do OSM (para todos os pontos):

  • man_made=survey_point
  • note=Não deslocar o ponto (coordenadas originais do IBGE)
  • note:en=Please: do not displace the point (official coordinates)

- - - - - - - - - -


Particularidades da conversão e validação de cada grupo no JOSM

1. Rede Brasileira de Monitoramento Contínuo - RBMC/GNSS Permanente

  • NOTAː Foram removidos do pacote total de nós os 129 nós da rede "IBGE-RBMC/GNSS", pois correspondem "aos mesmos 129 nós" dentre os da camada "IBGE-SAT GPS" (abaixo), sendo nestes adicionada a tag "ref:RBMC=*".
  • NÃO SERÁ MAIS UTILIZADA ESTA CAMADA "RBMC/GNSS Permanente".


Arquivo Original: BDG RBMCPoint.SHP / 4 KB

Website originalː http://www.visualizador.inde.gov.br/VisualizaCamada/674

ORIGINAL - Nº OBJETOSː 129

SELECIONADOS (SITUACAO=BOM)ː 129

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

TABELA DE ATRIBUTOS E CONVERSÃOː

(exemplo de linha de valores; campos com (*) indicam conversão para .OSM; demais a remover)

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

(*)ESTACAO=99614 >> converter p/ "ref=99614"

(*)TEMA=GPS >> converter p/ "description=Rede Geodésica do IBGE - Estação RBMC-GNSS/GPS"

UF=MS

(*)MUNICIPIO=Aquidauana >> converter p/ "refːloc=Aquidauana"

GEOCODIGO=5001102

VISITA=28/05/2016

(*)SITUACAO=BOM >> selecionar somente objetos com "SITUACAO=BOM"; ao final, remover o campo

CONEX_GPS=

CONEX_VT=

CONEX_EP=

CONEX_DOP=

CONEX_RN=

CONEX_EG=

LATITUDE=-20,4529433333000000000000000000000

LONGITUDE=-55,6711646389000000000000000000000

LATGMS=- 20 27 10,5960

LONGMS=- 55 40 16,1927

SIGMALAT=0,002

SIGMALON=0,002

FONTECOORD=GPS GEODESICO

(*)ALTGEOM=188,113 >> em SIRGAS2000 (conversão para WGS84=zero) - converter p/ "ele=188,113" (129 de 129)

ADIÇÃO de tagː "sourceːele =IBGE/ALTGEOM"

SIGMALTGEO=0,008

FONTALTGEO=GPS GEODESICO

MEDICAOPLA=28/05/2016

DATCALCPLA=15/06/2016

DATUM_PLAN=SIRGAS2000

ALTORTO=

SIGMALT_RN=

DATUM_ALT=

FONTALTORT=

MEDICAOALT=

CALCULOALT=

GRAVIDADE=

MEDICAOGRA=

CALCULOGRA=

DATUM_GRAV=

(*)URL_LINK=ftp://geoftp.ibge.gov.br/RBMC/relatorio/Descritivo_MSAQ.pdf ː(Links originais estão inativos)

>> substituir por novo url = * ː

concatenar adicionando o código da estação ao final, ex.ː http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=99614

NOTA=

(*)SG_RBMC=MSAQ >> converter p/ "refːRMBC=MSAQ"

NR_DOMES=

NM_RBMC=AQUIDAUANA

- - - - - - - - - - - - - - - - - - - - -

VALIDAÇÃO NO JOSMː

  • ISOLADA (somente os 129)ː 0 Errors / 0 Warnings
  • GERAL (com todos os 45622 objetos)ː 129 (todos RBMC) com "Warning = "nodes at same position" (em relação a nós de outras redes)

2. Estações SAT DOPPLER

Arquivo Original: BDG_DOPPoint.SHP / 28 KB

Website originalː http://www.visualizador.inde.gov.br/VisualizaCamada/661

ORIGINAL - Nº OBJETOSː 1006

SELECIONADOS (SITUACAO=BOM)ː 905

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

TABELA DE ATRIBUTOS E CONVERSÃOː

(exemplo de linha de valores; campos com (*) indicam conversão para .OSM; demais a remover)

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

(*)ESTACAO=90001 >> converter p/ "ref=90001"

(*)TEMA=DOPPLER >> converter p/ "description=Rede Geodésica do IBGE - Estação SAT DOPPLER"

UF=RO

(*)MUNICIPIO=Nova Mamoré >> converter p/ "refːloc=Nova Mamoré"

GEOCODIGO=1100338

VISITA=25/06/2013

(*)SITUACAO=BOM >> selecionar somente objetos com "SITUACAO=BOM"; ao final, remover o campo

CONEX_GPS=

CONEX_VT=

CONEX_EP=

CONEX_DOP=

CONEX_RN=9520R

CONEX_EG=

LATITUDE=-10,0778088889000000000000000000000

LONGITUDE=-65,3164000000000000000000000000000

LATGMS=- 10 04 40,112

LONGMS=- 65 18 59,040

SIGMALAT=

SIGMALON=

FONTECOORD=DOPPLER

(*)ALTGEOM=121,620 >> em SIRGAS2000 (conversão para WGS84=zero) - converter p/ "ele=121,620" (757 de 905)

ADIÇÃO de tagː "sourceːele =IBGE/ALTGEOM"

NOTAː Sendo que em 148 das estações não consta este dado, para "estas que faltam" (somente) será convertido de ALTORTO (H) com o programa "MAPGEO2015_v1_0.exe", do IBGE, gerando valores de (N) para obter (h), no campo " ele_mapgeo = * " (nome a ser convertido para "ele" posteriormente no JOSM .

SIGMALTGEO=

FONTALTGEO=DOPPLER

MEDICAOPLA=15/09/1973

DATCALCPLA=14/01/2005

DATUM_PLAN=SIRGAS2000

(*)ALTORTO=100,4959 >> manter os dados originais na tag "eleːorto" ; converter com o programa "MAPGEO2015_v1_0.exe" para (h) em " ele = * " (com 3 casas), somente para as 148 que não possuem valor em ALTGEOM. ADIÇÃO de tagː "sourceːele=IBGE/MAPGEO2015(ALTORTO + N)"


SIGMALT_RN=

DATUM_ALT=Imbituba

FONTALTORT=NIVELAMENTO GEOMÉTRICO

MEDICAOALT=15/11/1974

CALCULOALT=15/06/2011

GRAVIDADE=

MEDICAOGRA=

CALCULOGRA=

DATUM_GRAV=

(*)LINK=* >> converter p/ "url=*"

(concatenar adicionando o código da estação ao final, ex.ː http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=90001)

NOTA=A Referência de Nível (RN) pertence a um ramal, isto é, uma linha de nivelamento geométrico que não forma circuito, cujos desníveis, portanto, não podem ser ajustados. Consequentemente, as respectivas altitudes são calculadas mediante simples transporte

- - - - - - - - - - - - - - - - - - - - -

VALIDAÇÃO NO JOSMː

  • ISOLADA (somente os 905)ː 0 Errors / 0 Warnings
  • GERAL (com todos os 45622 objetos)ː 231 com "Warning = "nodes at same position" (em relação a nós de outras redes)

3. Estações SAT GPS

Arquivo Original: BDG BDG_GPSPoint.SHP / 82 KB

Website originalː http://www.visualizador.inde.gov.br/VisualizaCamada/664

ORIGINAL - Nº OBJETOSː 2977

SELECIONADOS (SITUACAO=BOM)ː 2823

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

TABELA DE ATRIBUTOS E CONVERSÃOː

(exemplo de linha de valores; campos com (*) indicam conversão para .OSM; demais a remover)

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

(*)ESTACAO=99541 >> converter p/ "ref=99541"

(*)TEMA=GPS >> converter p/ "description=Rede Geodésica do IBGE - Estação SAT GPS"

UF=MG

(*)MUNICIPIO=Uruana de Minas >> converter p/ "refːloc=Uruana de Minas"

GEOCODIGO=3170479

VISITA=10/09/2011

(*)SITUACAO=BOM >> selecionar somente objetos com "SITUACAO=BOM"; ao final, remover o campo

CONEX_GPS=

CONEX_VT=

CONEX_EP=

CONEX_DOP=

CONEX_RN=

CONEX_EG=

LATITUDE=-16,0587896944000000000000000000000

LONGITUDE=-46,2555510833000000000000000000000

LATGMS=- 16 03 31,6429

LONGMS=- 46 15 19,9839

SIGMALAT=0,001

SIGMALON=0,001

FONTECOORD=GPS GEODESICO

(*)ALTGEOM=533,478 >> em SIRGAS2000 (conversão para WGS84=zero) - converter p/ "ele=533,478" (2817 de 2823)

ADIÇÃO de tagː "sourceːele =IBGE/ALTGEOM"

SIGMALTGEO=0,008

FONTALTGEO=GPS GEODESICO

MEDICAOPLA=10/09/2011

DATCALCPLA=01/11/2011

DATUM_PLAN=SIRGAS2000

(*)ALTORTO= * >> manter os dados originais na tag "eleːorto" ; Não será possível converter para (h)ː as mesmas 06 que não possuem ALTGEOM também não possuem ALTORTO

SIGMALT_RN=

DATUM_ALT=

FONTALTORT=

MEDICAOALT=

CALCULOALT=

GRAVIDADE=

MEDICAOGRA=

CALCULOGRA=

DATUM_GRAV=

(*)LINK=* >> converter p/ "url=*"

(concatenar o código da estação ao final, ex.ː http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=99541)

NOTA=

- - - - - - - - - - - - - - - - - - - - -

VALIDAÇÃO NO JOSMː

  • ISOLADA (somente os 2823)ː 0 Errors / 0 Warnings
  • GERAL (com todos os 45622 objetos)ː 1098 com "Warning = "nodes at same position" (em relação a nós de outras redes)

4. Estações de Poligonal

Arquivo Original: BDG BDG_EPPoint.SHP / 32 KB

Website originalː http://www.visualizador.inde.gov.br/VisualizaCamada/663

ORIGINAL - Nº OBJETOSː 1133

SELECIONADOS (SITUACAO=BOM)ː 1006

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

TABELA DE ATRIBUTOS E CONVERSÃOː

(exemplo de linha de valores; campos com (*) indicam conversão para .OSM; demais a remover)

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

(*)ESTACAO=10025 >> converter p/ "ref=10025"

(*)TEMA=EP >> converter p/ "description=Rede Geodésica do IBGE - Estação de Poligonal"

UF=AP

(*)MUNICIPIO=Ferreira Gomes >> converter p/ "refːloc=Ferreira Gomes"

GEOCODIGO=1600238

VISITA=24/09/2014

(*)SITUACAO=BOM >> selecionar somente objetos com "SITUACAO=BOM"; ao final, remover o campo

CONEX_GPS=

CONEX_VT=

CONEX_EP=

CONEX_DOP=

CONEX_RN=9533S

CONEX_EG=

LATITUDE=0,9399059444000000000000000000000

LONGITUDE=-51,1461391389000000000000000000000

LATGMS=+ 00 56 23,6614

LONGMS=- 51 08 46,1009

SIGMALAT=0,262

SIGMALON=0,409

FONTECOORD=POLIGONAÇÃO

(*)ALTGEOM=* >> em SIRGAS2000 (conversão para WGS84=zero) - converter p/ "ele=*" (32 de 1006)

ADIÇÃO de tagː "sourceːele =IBGE/ALTGEOM"

NOTAː Sendo que em 974 das estações não consta este dado, para "estas que faltam" (somente) será convertido de ALTORTO (H) com o programa "MAPGEO2015_v1_0.exe", do IBGE, gerando valores de (N) para obter (h), no campo " ele_mapgeo = * " (nome a ser convertido para "ele" posteriormente no JOSM .

SIGMALTGEO=

FONTALTGEO=

MEDICAOPLA=15/06/1961

DATCALCPLA=06/03/2006

DATUM_PLAN=SIRGAS2000

(*)ALTORTO=73,1136 >> manter os dados originais na tag "eleːorto" ; converter com o programa "MAPGEO2015_v1_0.exe" para (h) em " ele = * " (com 3 casas), somente para as demais 974 que não possuem valor em ALTGEOM. ADIÇÃO de tagː "sourceːele =IBGE/MAPGEO2015(ALTORTO + N)"

SIGMALT_RN=

DATUM_ALT=Santana

FONTALTORT=NIVELAMENTO GEOMETRICO

MEDICAOALT=18/02/1981

CALCULOALT=15/06/2011

GRAVIDADE=

MEDICAOGRA=

CALCULOGRA=

DATUM_GRAV=

(*)LINK=* >> converter p/ "url=*"

(concatenar adicionando o código da estação ao final, ex.ː http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=10025)

NOTA=A Referência de Nível (RN) pertence a um ramal, isto é,...

- - - - - - - - - - - - - - - - - - - - -

VALIDAÇÃO NO JOSMː

  • ISOLADA (somente os 1006)ː 0 Errors / 0 Warnings
  • GERAL (com todos os 45622 objetos)ː 169 com "Warning = "nodes at same position" (em relação a nós de outras redes)

5. Vértices de Triangulação

Arquivo Original: BDG_VTPoint.SHP / 100 KB

Website originalː http://www.visualizador.inde.gov.br/VisualizaCamada/666

ORIGINAL - Nº OBJETOSː 3643

SELECIONADOS (SITUACAO=BOM)ː 2814

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

TABELA DE ATRIBUTOS E CONVERSÃOː

(exemplo de linha de valores; campos com (*) indicam conversão para .OSM; demais a remover)

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

(*)ESTACAO=1 >> converter p/ "ref=1"

(*)TEMA=VT >> converter p/ "description=Rede Geodésica do IBGE - Vértice de Triangulação"

UF=RS

(*)MUNICIPIO=Mampituba >> converter p/ "refːloc=Mampituba"

GEOCODIGO=4311734

VISITA=15/12/1950

(*)SITUACAO=BOM >> selecionar somente objetos com "SITUACAO=BOM"; ao final, remover o campo

CONEX_GPS=

CONEX_VT=

CONEX_EP=

CONEX_DOP=

CONEX_RN=

CONEX_EG=

LATITUDE=-29,2425334722000000000000000000000

LONGITUDE=-49,9130385278000000000000000000000

LATGMS=- 29 14 33,1205

LONGMS=- 49 54 46,9387

SIGMALAT=0,389

SIGMALON=0,379

FONTECOORD=TRIANGULAÇÃO

(*)ALTGEOM=* >> em SIRGAS2000 (conversão para WGS84=zero) - converter p/ "ele=*" (83 de 2814).

ADIÇÃO de tagː "sourceːele =IBGE/ALTGEOM"

NOTAː Sendo que em 2731 das estações não consta este dado, para "estas que faltam" (somente) será convertido de ALTORTO (H) com o programa "MAPGEO2015_v1_0.exe", do IBGE, gerando valores de (N) para obter (h), no campo " ele_mapgeo = * " (nome a ser convertido para "ele" posteriormente no JOSM .

SIGMALTGEO=

FONTALTGEO=

MEDICAOPLA=15/12/1950

DATCALCPLA=06/03/2006

DATUM_PLAN=SIRGAS2000

(*)ALTORTO=235,26 >> manter os dados originais na tag "eleːorto" ; converter com o programa "MAPGEO2015_v1_0.exe" para (h) em " ele = * " (com 3 casas), somente para as demais 2731 que não possuem valor em ALTGEOM. ADIÇÃO de tagː "sourceːele =IBGE/MAPGEO2015(ALTORTO + N)"

SIGMALT_RN=

DATUM_ALT=Imbituba

FONTALTORT=NIVELAMENTO TRIGONOMÉTRICO

MEDICAOALT=15/12/1950

CALCULOALT=15/12/1950

GRAVIDADE=

MEDICAOGRA=

CALCULOGRA=

DATUM_GRAV=

(*)LINK=* >> converter p/ "url=*"

(concatenar adicionando o código da estação ao final, ex.ː http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=1)

NOTA=

- - - - - - - - - - - - - - - - - - - - -

VALIDAÇÃO NO JOSMː

  • ISOLADA (somente os 2814)ː 0 Errors / 0 Warnings
  • GERAL (com todos os 45622 objetos)ː 284 com "Warning = "nodes at same position" (em relação a nós de outras redes)

6. Referências de Nível (Rede Geodésica Altimétrica)

Arquivo Original: BDG_RNPoint.SHP / 613 MB

Website originalː http://www.visualizador.inde.gov.br/VisualizaCamada/665

ORIGINAL - Nº OBJETOSː 69.887

SELECIONADOS (SITUACAO=BOM)ː 37.945

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Problema encontradoː

  • Maioria dos pontos NÃO possui valor de "ALTGEOM" (h = Altitude Geométrica / GPS / WGS84).
  • Mas TODOS possuem valor de "ALTORTO" (H = Altitude Ortométrica / Geoidal).

Soluçãoː

Usar o cálculo "H = h - N", conforme IBGE / Modelo de Ondulação Geoidal do Brasil (em "http://www.ibge.gov.br/home/geociencias/geodesia/calculo_do_modelo.shtm"), para gerar "ALTGEOM (h)" a partir de "ALTORTO (H)" (h = H + N);

  • Utilizado o programa "MAPGEO2015_v1_0.exe", do IBGE (fornecido no mesmo site acima) para gerar valores de "N" para cada um dos pontos;
  • Concatenados valores de "N" no DBF original;
  • Adicionado novo campo "ALTGEOM-CONVERT" no SHP, com o cálculo dos valores de (h)ː "h = H + N".
  • Convertido campo "ALTGEOM-CONVERT" para tag .OSM " ele = * ".

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

TABELA DE ATRIBUTOS E CONVERSÃOː

(2 exemplos de linhas de valores; campos com (*) indicam conversão para .OSM; demais a remover)

- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

(*)ESTACAO = 144I  ; 146Q >> converter p/ "ref=*"

(*)TEMA = RN  ; RN >> converter p/ "description=Rede Geodésica do IBGE - Referência de Nível"

UF = MG  ; MG

(*)MUNICIPIO = São João Evangelista  ; Patos de Minas >> converter p/ "refːloc=*"

GEOCODIGO = 3162807  ; 3148004

VISITA = 04/05/2001  ; 13/07/2009

(*)SITUACAO = BOM  ; BOM >> selecionar somente objetos com "SITUACAO=BOM"; ao final, remover o campo

CONEX_GPS = <null>  ; 96056

CONEX_VT = <null>  ; <null>

CONEX_EP = <null>  ; <null>

CONEX_DOP = <null>  ; <null>

CONEX_RN = <null>  ; <null>

CONEX_EG = <null>  ; <null>

LATITUDE = -18,5838888889000000000000000000000  ; -18,5894766111000000000000000000000

LONGITUDE = -42,7836111111000000000000000000000  ; -46,5169277778000000000000000000000

LATGMS = - 18 35 02  ; - 18 35 22,1158

LONGMS = - 42 47 01  ; - 46 31 00,9400

SIGMALAT = <null>  ; 0,002

SIGMALON = <null>  ; 0,003

FONTECOORD = GPS NAVEGAÇÃO  ; GPS GEODESICO

(*)ALTGEOM = <null>  ; 824,346 >>

maioria não possui este dadoː precisa ser convertido a partir de ALTORTO (H)ː "h = H + N"

ADIÇÃO de tagː "sourceːele =IBGE/ALTGEOM"

NOTAː Sendo que em 36.929 das 37.945 estações não consta este dado, para "estas 36.929 que faltam" (somente) será convertido de ALTORTO (H) com o programa "MAPGEO2015_v1_0.exe", do IBGE, gerando valores de (N) para obter (h), no campo " ele_mapgeo = * " (nome a ser convertido para "ele" posteriormente no JOSM . ADIÇÃO de tagː "sourceːele =IBGE/ALTGEOM"

Notaː há uma pequena diferença em função da atualização do Modelo 2015/IBGEː o valor calculado em função do (N) do Modelo 2015 (N=-9.56) para a estação acima é 824,460; difere em relação ao mesmo valor calculado em 15/06/2011 (824,346) por modelo anterior, cf. data abaixo. RESULTADOː 1358 pontos com altitude negativa (ele = variando de -0.0001 a -21.801), entre os 37945. Convertido com precisão de 3 casas decimais (exː-21.801).

SIGMALTGEO = <null>  ; 0,017

FONTALTGEO = <null>  ; GPS GEODESICO

MEDICAOPLA = 04/05/2001  ; 04/10/2008

DATCALCPLA = <null>  ; 04/12/2008

DATUM_PLAN = SIRGAS2000  ; SIRGAS2000

(*)ALTORTO = 727,9531  ; 834,0203 >> todos possuem este dado ; manter os dados originais na tag "eleːorto" ; converter com o programa "MAPGEO2015_v1_0.exe" para (h) em " ele = * " (com 3 casas), somente para as 1016 demais que não possuem valor em ALTGEOM. ADIÇÃO de tagː "sourceːele =IBGE/MAPGEO2015(ALTORTO + N)"

manter os dados originais de (N) na tag "refːN_mapgeo = * "

SIGMALT_RN = 0,08  ; 0,07

DATUM_ALT = Imbituba  ; Imbituba

FONTALTORT = NIVELAMENTO GEOMETRICO  ; NIVELAMENTO GEOMETRICO

MEDICAOALT = 18/05/1951  ; 15/07/1951

CALCULOALT = 15/06/2011  ; 15/06/2011

GRAVIDADE = <null>  ; <null>

MEDICAOGRA = <null>  ; <null>

CALCULOGRA = <null>  ; <null>

DATUM_GRAV = <null>  ; <null>

(*)LINK=* >> converter p/ "url=*"

(concatenar adicionando o código da estação ao final, ex.ː http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=144I)

NOTA = <null>  ; <null>

- - - - - - - - - - - - - - - - - - - - -

VALIDAÇÃO NO JOSMː

  • ISOLADA (somente os 37.945)ː Errors (0) / Warnings (3) - Nodes at same position (3)ː
  • 2 nodes: 2680Z, 2781P
  • 4 nodes: 2721V, 2721X, 2721U,2721T
  • 4 nodes: 3018J, 3018S, 3018T, 3018L

NOTAː Relatórios originais indicam pontos com mesmas coordenadas; caso se pretenda resolver fazendo "merge nodes" (JOSM>M), perde-se os seguintes dados individuais e exclusivos, p.ex.ː

CONCLUSÃOː Até agora mantidos como estão.

  • GERAL (com todos os 45622 objetos)ː 1203 com "Warning = "nodes at same position" (em relação a nós de outras redes)

Validação Geral no JOSM e objeto .OSM

Validação internaː

Somente entre o grupo todo dos 6 tipos de Redes Geodésicas (45.622 nós) - sem baixar dados atuais do OSMː

VALIDAÇÃO NO JOSMː Errors (0) / Warnings (1579) - Nodes at same position (1579)ː

Tratam-se de estações das 6 diferentes redes, referentes a distintos SHP, sendo 1579 casos de mesmas coordenadas, num total de 3240 nós, sobre o mesmo elemento físico. Por exemploː

NOTAː Relatórios originais indicam pontos com mesmas coordenadas; caso se pretenda resolver fazendo "merge nodes" (JOSM>M), perde-se os seguintes dados individuais e exclusivos, p.ex. estes 3 nós em "Warning - Nodes at same position"ː

Questãoː manter pontos originais distintos? ou forçar "merge" dos dados em 1 só ponto? Se fundir, não permite que alguém encontre com uma simples busca pela "ref = *" a rede. Como identificar se fizer "merge"?

Propostaː manter como está nos originais, sem fusão (merge) de pontos.

Conclusãoː Até agora mantidos como estão.


Validação externaː

Querry Overpass no Brasil todo porː "man_made=survey_point and typeːnode"

Foram apontados 17 casos de objetos previamente existentes no OSM que apresentam "Warningː Nodes at same position" com os novos pontos da Rede Geodésica. A verificar possibilidade de fusão com "replace geometry" no JOSM, conservando histórico dos objetos, verificando as tags "ref=*" ː

Objetos existentes em conflito (CSV - via overpass)ː

  • @id;ref;name;source
  • node/3595547047;91950;;IBGE
  • node/3577081800;3279M;3279M;
  • node/3577081803;3279X;3279X;
  • node/3577081793;3279R;3279R;
  • node/3595097351;SAT 94048;RSAL;IBGE
  • node/3577081802;3279V;3279V;
  • node/3595534060;91938;;IBGE
  • node/3577081801;PR23;3279U;
  • node/3577081796;3092U;3092U;
  • node/3577081804;3279Z;3279Z;
  • node/3634429654;;Estação Geodésica SAT 93500;
  • node/3577081798;3278A;3278A;
  • node/3577081795;3279T;3279T;
  • node/3577081092;PR24;PR24;
  • node/3577081797;3092X;3092X;
  • node/3935936558;;SAT 93501;
  • node/3577081794;3279S;3279S;

ID Changeset

Proposta - a ser importado pela conta dedicadaː http://www.openstreetmap.org/user/smaprs_import (deː http://www.openstreetmap.org/user/smaprs)

Número total de objetosː 45.622

Sourceː IBGE/INDE

Changesetsː

1) Redes Geodésicas do IBGE - Rede Brasileira de Monitoramento Contínuo - GNSS (129 Objetos)

Comment = "Redes Geodésicas do IBGE - Rede Brasileira de Monitoramento Contínuo - GNSS - Documented atː https://wiki.openstreetmap.org/wiki/Importa%C3%A7%C3%A3o_das_Redes_Geod%C3%A9sicas_do_IBGE"

Changeset =

2) Redes Geodésicas do IBGE - Estações SAT DOPPLER (905 Objetos)

Comment = "Importː Redes Geodésicas do IBGE - Estações SAT DOPPLER - Documented atː https://wiki.openstreetmap.org/wiki/Importa%C3%A7%C3%A3o_das_Redes_Geod%C3%A9sicas_do_IBGE"

Changeset =

3) Redes Geodésicas do IBGE - Estações SAT GPS (2823 Objetos)

Comment = "Importː Redes Geodésicas do IBGE - Estações SAT GPS - Documented atː https://wiki.openstreetmap.org/wiki/Importa%C3%A7%C3%A3o_das_Redes_Geod%C3%A9sicas_do_IBGE"

Changeset =

4) Redes Geodésicas do IBGE - Estações de Poligonal (1006 Objetos)

Comment = "Importː Redes Geodésicas do IBGE - Estações de Poligonal - Documented atː https://wiki.openstreetmap.org/wiki/Importa%C3%A7%C3%A3o_das_Redes_Geod%C3%A9sicas_do_IBGE"

Changeset =

5) Redes Geodésicas do IBGE - Vértices de Triangulação (2814 Objetos)

Comment = "Importː Redes Geodésicas do IBGE - Vértices de Triangulação - Documented atː https://wiki.openstreetmap.org/wiki/Importa%C3%A7%C3%A3o_das_Redes_Geod%C3%A9sicas_do_IBGE"

Changeset =

6) Redes Geodésicas do IBGE - Referências de Nível (37.945 Objetos)

Comment = "Importː Redes Geodésicas do IBGE - Referências de Nível - Documented atː https://wiki.openstreetmap.org/wiki/Importa%C3%A7%C3%A3o_das_Redes_Geod%C3%A9sicas_do_IBGE"

Changeset =