Hfrhyu

Geometric mean and geometric standard deviation

Unexpected result with right shift after bitwise negation

How should I respond to a player wanting to catch a sword between their hands?

Fishing simulator

What to do with post with dry rot?

Is there folklore associating late breastfeeding with low intelligence and/or gullibility?

Blender game recording at the wrong time

Am I ethically obligated to go into work on an off day if the reason is sudden?

Can I throw a sword that doesn't have the Thrown property at someone?

What's the difference between (size_t)-1 and ~0?

I'm having difficulty getting my players to do stuff in a sandbox campaign

Why is "Captain Marvel" translated as male in Portugal?

Do working physicists consider Newtonian mechanics to be "falsified"?

When communicating altitude with a '9' in it, should it be pronounced "nine hundred" or "niner hundred"?

Classification of bundles, Postnikov towers, obstruction theory, local coefficients

Complexity of many constant time steps with occasional logarithmic steps

How can I make names more distinctive without making them longer?

Passing functions in C++

How do I keep my slimes from escaping their pens?

How can players take actions together that are impossible otherwise?

What kind of display is this?

What items from the Roman-age tech-level could be used to deter all creatures from entering a small area?

What would be Julian Assange's expected punishment, on the current English criminal law?

How to colour the US map with Yellow, Green, Red and Blue to minimize the number of states with the colour of Green

Why shapefiles created with sf package (in R) do not display same CRS in ArcMap?

.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;

0

I have created a polygon shapefile in ArcMap using 7 points(lat,lon) then read this shapefile into R to obtain the CRS. I then created a shape file in R using the sf package (and the 7 points) and used the CRS obtained from the ArcMap shapefile. When I bring the shapefile created in R into ArcMap I get a warning regarding different coordinate systems. Why would this happen?

My steps:

7 points (lon and lat) where point 1 = point 7 (closed loop). Create a polygon shapefile in ArcMap and display. Properties in ArcMap show : Geographic Coordinate System: GCS_North_American_1983 Datum: D_North_American_1983

In R I use rgdal to read in shapefile and obtain CRS

cashes <- rgdal::readOGR(dns="path",layer="shapefile") 
cashes@proj4string 
This results in : "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 "

Now I create the shapefile in R using the sf package (I could use %>%)

Rshapefile <- readxl::read_xlsx(path,sheet="cashes")
Rshapefile <- sf::st_as_sf(Rshapefile,coords=c("lon","lat"),crs="+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 ")
Rshapefile <- dplyr::summarize(Rshapefile,geometry = sf::st_combine(geometry))
Rshapefile <- sf::st_cast(Rshapefile,"POLYGON")

sf::st_write(Rshapefile,dns=path,layer="Rshapefile",driver="ESRI Shapefile")

I then bring this Rshapefile into ArcMap and the coordinate systems are different. When I view the properties in ArcMap they show: Geographic Coordinate System: GRS 1980(IUGG, 1980) Datum: D_unknown

Why is this? I must be doing something wrong but I can't figure out what!

coordinate-system arcmap r sf

share|improve this question

edited Apr 8 at 14:20

Kadir Şahbaz

4,69221531

asked Apr 8 at 13:54

andyandy

1

New contributor

andy is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  I can't help you with this question, but thank you for including the references to the packages for each function.
  
  – Paul H
  Apr 8 at 14:16
  
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  Have you checked the CRS at each stage of manipulating Rshapefile after first assigning the CRS?
  
  – Simbamangu
  Apr 8 at 16:07
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  dns=path should be dsn=path. Not sure that fixes it but at least it makes your code run properly.
  
  – Spacedman
  Apr 8 at 18:09
  

  
  
  
  

add a comment | 

0

I have created a polygon shapefile in ArcMap using 7 points(lat,lon) then read this shapefile into R to obtain the CRS. I then created a shape file in R using the sf package (and the 7 points) and used the CRS obtained from the ArcMap shapefile. When I bring the shapefile created in R into ArcMap I get a warning regarding different coordinate systems. Why would this happen?

My steps:

7 points (lon and lat) where point 1 = point 7 (closed loop). Create a polygon shapefile in ArcMap and display. Properties in ArcMap show : Geographic Coordinate System: GCS_North_American_1983 Datum: D_North_American_1983

In R I use rgdal to read in shapefile and obtain CRS

cashes <- rgdal::readOGR(dns="path",layer="shapefile") 
cashes@proj4string 
This results in : "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 "

Now I create the shapefile in R using the sf package (I could use %>%)

Rshapefile <- readxl::read_xlsx(path,sheet="cashes")
Rshapefile <- sf::st_as_sf(Rshapefile,coords=c("lon","lat"),crs="+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 ")
Rshapefile <- dplyr::summarize(Rshapefile,geometry = sf::st_combine(geometry))
Rshapefile <- sf::st_cast(Rshapefile,"POLYGON")

sf::st_write(Rshapefile,dns=path,layer="Rshapefile",driver="ESRI Shapefile")

I then bring this Rshapefile into ArcMap and the coordinate systems are different. When I view the properties in ArcMap they show: Geographic Coordinate System: GRS 1980(IUGG, 1980) Datum: D_unknown

Why is this? I must be doing something wrong but I can't figure out what!

coordinate-system arcmap r sf

share|improve this question

edited Apr 8 at 14:20

Kadir Şahbaz

4,69221531

asked Apr 8 at 13:54

andyandy

1

New contributor

andy is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  I can't help you with this question, but thank you for including the references to the packages for each function.
  
  – Paul H
  Apr 8 at 14:16
  
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  Have you checked the CRS at each stage of manipulating Rshapefile after first assigning the CRS?
  
  – Simbamangu
  Apr 8 at 16:07
  

  
  
  
  


• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  dns=path should be dsn=path. Not sure that fixes it but at least it makes your code run properly.
  
  – Spacedman
  Apr 8 at 18:09
  

  
  
  
  

add a comment | 

I have created a polygon shapefile in ArcMap using 7 points(lat,lon) then read this shapefile into R to obtain the CRS. I then created a shape file in R using the sf package (and the 7 points) and used the CRS obtained from the ArcMap shapefile. When I bring the shapefile created in R into ArcMap I get a warning regarding different coordinate systems. Why would this happen?

My steps:

7 points (lon and lat) where point 1 = point 7 (closed loop). Create a polygon shapefile in ArcMap and display. Properties in ArcMap show : Geographic Coordinate System: GCS_North_American_1983 Datum: D_North_American_1983

In R I use rgdal to read in shapefile and obtain CRS

cashes <- rgdal::readOGR(dns="path",layer="shapefile") 
cashes@proj4string 
This results in : "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 "

Now I create the shapefile in R using the sf package (I could use %>%)

Rshapefile <- readxl::read_xlsx(path,sheet="cashes")
Rshapefile <- sf::st_as_sf(Rshapefile,coords=c("lon","lat"),crs="+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 ")
Rshapefile <- dplyr::summarize(Rshapefile,geometry = sf::st_combine(geometry))
Rshapefile <- sf::st_cast(Rshapefile,"POLYGON")

sf::st_write(Rshapefile,dns=path,layer="Rshapefile",driver="ESRI Shapefile")

I then bring this Rshapefile into ArcMap and the coordinate systems are different. When I view the properties in ArcMap they show: Geographic Coordinate System: GRS 1980(IUGG, 1980) Datum: D_unknown

Why is this? I must be doing something wrong but I can't figure out what!

coordinate-system arcmap r sf

share|improve this question

edited Apr 8 at 14:20

Kadir Şahbaz

4,69221531

asked Apr 8 at 13:54

andyandy

1

New contributor

andy is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

cashes <- rgdal::readOGR(dns="path",layer="shapefile") 
cashes@proj4string 
This results in : "+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 "

Rshapefile <- readxl::read_xlsx(path,sheet="cashes")
Rshapefile <- sf::st_as_sf(Rshapefile,coords=c("lon","lat"),crs="+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0 ")
Rshapefile <- dplyr::summarize(Rshapefile,geometry = sf::st_combine(geometry))
Rshapefile <- sf::st_cast(Rshapefile,"POLYGON")

sf::st_write(Rshapefile,dns=path,layer="Rshapefile",driver="ESRI Shapefile")

share|improve this question

edited Apr 8 at 14:20

Kadir Şahbaz

4,69221531

asked Apr 8 at 13:54

andyandy

1

New contributor

andy is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

share|improve this question

edited Apr 8 at 14:20

Kadir Şahbaz

4,69221531

asked Apr 8 at 13:54

andyandy

1

New contributor

andy is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

edited Apr 8 at 14:20

Kadir Şahbaz

4,69221531

edited Apr 8 at 14:20

Kadir Şahbaz

4,69221531

asked Apr 8 at 13:54

andyandy

1

New contributor

andy is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.

asked Apr 8 at 13:54

andyandy

1