[R-br] Problema na conversão de datum: SIRGAS2000 para WGS84

[José Lucas Safanelli]

Bia,

Acredito que não seja um problema do R ou pacote 'sf'.

Digo isso porque fazendo as mesmas operações que você, percebi que não
houve nenhuma alteração do CRS. Você pode verificar isso utilizando a
função st_crs().

Além disso, você pode exportar teus dados espaciais usando o próprio pacote
sf com o comando st_write.

library(sf)

dadosibge <- st_transform(municipios, 4326)
st_crs(municipios)
st_crs(dadosibge)

df <- inner_join(dadosibge, df, by= c("geocodigo"="cod_municipio")) %>%
st_as_sf(coords=c("lon","lat"), crs = 4326)

st_crs(df)

st_write(df, dsn = getwd(),
         layer = "teste", driver = 'ESRI Shapefile')

teste <- st_read(dsn = getwd(),
                 layer = "teste")

st_crs(teste)

Gostaria de saber como você verificou essa inconsistência do CRS, foi
utilizando um SIG externo? Penso que talvez tenha acontecido isso porque a
nomenclatura descreve o termo visualization, então pode ser que teu SIG
tenha trocado automático para talvez otimizar a visualização, não!? Acho
que não somente o nome, mas seria bom conferir se os parâmetros elipsoidais
se alteram com essa pequena mudança. Enfim, precisa investigar melhor.

José Lucas Safanelli
https://zecojls.netlify.app/
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20201204/8347ae3b/attachment.html>