<div dir="ltr">Olá pessoal!<div><br></div><div>Após a produção de shapefiles no R, eles apresentaram o seguinte datum:</div><div><br></div><div><div>Geographic Coordinate System: GCS_WGS_84_with_axis_order_normalized_for_visualization<br>Datum: D_WGS_1984</div><div><br></div><div>Sendo que o que preciso seria:</div><div><br></div><div>Geographic Coordinate System: GCS_WGS_1984<br>Datum: D_WGS_1984</div></div><div><br></div><div>Ou seja, preciso fazer a transformação sem gerar essa extensão de "with_axis_order_normalized_for_visualization", mas não consigo encontrar uma solução.</div><div><br></div><div>Resumindo o que fiz: <br>1. Download de dados municipais do IBGE (o datum deles é o SIRGAS2000)<br><br></div><div>2. Transformação dos dados do IBGE para o datum que preciso (WGS84) com a função 'st_transform':<br><br>dadosibge <- st_transform(dadosibge, 4326)<br><br></div><div>3. União destes dados do IBGE com minha planilha de dados, mais a transformação de objeto st em sf:<br><br></div><div>df <- inner_join(dadosibge, df, by="id") %>%<br>    st_as_sf(coords=c("lon","lat"), crs=4326)<br></div><div><br></div><div>4. Transformação do data frame em objeto espacial:</div><div><br></div><div>df <- as(df, "Spatial")<br></div><div><br></div><div>5. Exportação do shapefile:</div><div><br></div><div>writeOGR(df, dsn = "diretorio", layer = i, driver = 'ESRI Shapefile')<br></div><div> </div><div>Obs.: fiz com múltiplas planilhas, então o código está em um loop.</div><div><br></div><div>Segue o meu código:</div><div><br></div><div>###Pacotes<br>x<-c("geobr","ggplot2","dplyr","RColorBrewer","ggspatial","maps","mapdata","rworldmap","maptools","mapproj","ggmap","sf","raster","spData","spDataLarge", "rgdal")<br>lapply(x, require, character.only=T)<br><br>###Carregar dados IBGE<br>##All municipalities limits of all states<br>spatial <- read_municipality(code_muni = 'all', year = 2019)<br><br>##Check CRS - If this fuction return "NA", the CRS is unknown<br>st_is_longlat(spatial)<br><br>##Verify the DATUM. In this case, shows SIRGAS 2000<br>st_crs(spatial)<br><br>##Transform SGC SIRGAS 2000 to WGS84 <br>###WGS84 EPSG code: 4326<br>spatial <- st_transform(spatial, 4326)<br>st_crs(spatial)<br>class(spatial)<br><br>###Carregar diretorio dos meus dados<br>setwd('diretorio')<br>list.files()<br>alldata = list.files()<br><br>###Loop para exportar varios shapefiles<br>for (i in alldata) {<br>  df = read.csv(i, sep = ";", stringsAsFactors = FALSE) %>%<br>   <br>#Unindo meus dados com o df do IBGE<br>  df <- inner_join(spatial_f, df, by="id") %>%<br>    st_as_sf(coords=c("lon","lat"), crs=4326)<br>  <br>#Transformando df em dados espaciais<br>  df <- as(df, "Spatial")<br>  <br>  #Exportando shapefile<br>  writeOGR(df, dsn = "C:/Users/Bianca/Desktop/ipam_consult/shapes_test", layer = i, driver = 'ESRI Shapefile')<br>  <br>}<br><br>Agradeço desde já!</div><div><br></div></div>