<div dir="ltr"><div>Bia,</div><div><br></div><div>Acredito que não seja um problema do R ou pacote 'sf'.<br></div><div><br></div><div>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().</div><div><br></div><div>Além disso, você pode exportar teus dados espaciais usando o próprio pacote sf com o comando st_write.<br></div><div><br></div><div>library(sf)<br></div><div><br></div><div>dadosibge <- st_transform(municipios, 4326)<br>st_crs(municipios)<br>st_crs(dadosibge)<br><br>df <- inner_join(dadosibge, df, by= c("geocodigo"="cod_municipio")) %>% st_as_sf(coords=c("lon","lat"), crs = 4326)<br></div><div><br></div><div>st_crs(df)<br><br>st_write(df, dsn = getwd(),<br>         layer = "teste", driver = 'ESRI Shapefile')<br><br>teste <- st_read(dsn = getwd(),<br>                 layer = "teste")</div><div><br>st_crs(teste)<br></div><div><br></div><div>
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.</div><div><br>

</div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr">José Lucas Safanelli</div><div dir="ltr"><a href="https://zecojls.netlify.app/" target="_blank">https://zecojls.netlify.app/</a><br></div></div></div></div></div></div></div>