Importação das Redes Geodésicas do IBGE
(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")
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)
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ː
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)
- url=* (concatenar a url de relatórios do IBGE (http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=<*>) com ref = * da estação) - notaː todos links testados por amostragem aleatória (cerca de 100 testados) estão funcionando
- 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
- source=IBGE/INDE/ + < http://www.visualizador.inde.gov.br/VisualizaCamada/ * (* cf. camada fonte do SHP)
- note=Não deslocar o ponto (coordenadas originais do IBGE)
- note:en=Please: do not displace the point (official coordinates)
- - - - - - - - - -
- Quando contiver "ele_mapgeo"ː adicionar sourceːele = IBGE/MAPGEO2015 ; nestes houve conversão de valores de altitude Ortométrica para Geométricaː "h = H + N" (a partir de H = h - N) conforme Modelo de Ondulação Geoidal/IBGE (http://www.ibge.gov.br/home/geociencias/geodesia/calculo_do_modelo.shtm).
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.ː
- ele:orto=23.9010
- ele:orto=23.9092
- ele=15.181
- ele=15.189
- ref=2721T ("Chapa cravada no canto norte da ...antena parabólica")
- ref=2721X ("Chapa cravada no canto leste da ...antena parabólica")
- url=http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=2721T
- url=http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=2721X
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"ː
- description=IBGE-Rede Altimétrica-RN
- description=IBGE-SAT GPS
- description=IBGE-VTriangulacao
- ele:orto=1168.381 (1) (mesmo dado abaixo, apenas precisão em 3 casas)
- ele:orto=1168.3814 (2)
- ele=1155.874 (3) (mesmo dado)
- ref=1609
- ref=9501V
- ref=96123
- url=http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=1609
- url=http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=9501V
- url=http://www.bdg.ibge.gov.br/bdg/pdf/relatorio.asp?L1=96123
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 =