Hfrhyu

"to be prejudice towards/against someone" vs "to be prejudiced against/towards someone"

How is the claim "I am in New York only if I am in America" the same as "If I am in New York, then I am in America?

How is it possible to have an ability score that is less than 3?

Why was the small council so happy for Tyrion to become the Master of Coin?

How do we improve the relationship with a client software team that performs poorly and is becoming less collaborative?

Is it possible to do 50 km distance without any previous training?

Today is the Center

LaTeX closing $ signs makes cursor jump

Characters won't fit in table

Why can't I see bouncing of a switch on an oscilloscope?

What defenses are there against being summoned by the Gate spell?

Why, historically, did Gödel think CH was false?

Test if tikzmark exists on same page

How much RAM could one put in a typical 80386 setup?

Modeling an IP Address

Why not use SQL instead of GraphQL?

The use of multiple foreign keys on same column in SQL Server

What typically incentivizes a professor to change jobs to a lower ranking university?

tikz: show 0 at the axis origin

Can an x86 CPU running in real mode be considered to be basically an 8086 CPU?

Why Is Death Allowed In the Matrix?

Is it unprofessional to ask if a job posting on GlassDoor is real?

How do I create uniquely male characters?

What are these boxed doors outside store fronts in New York?

gIntersection error: linearRing not created

Reprojecting Google Maps (latlong) to Irish Grid TM65 using spTransform in R?Transforming and Projecting Shapefiles in R & ArcGIS for Desktop?Crop a raster by a grid shapefile with different extensions in RProblem points into a shape polygonCRS Projection problems prevent spatial merge in RGetting “Error in par(pin = p)” when using spTransform in RMissing 1.6% of expected area after buffering SpatialPoints from vector of widthsUnexpected output from rasterToPolygon in R raster packageCompute a distance matrix in R with geosphereAddressing Z-dimension discarded when reading in shapefiles from ArcMap to R?

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

2

I am trying to do something relatively simple but with large shapefiles. I have a shapefile of protected areas in the USA:

and watershed boundaries:

I want to calculate the percentage and proportion of protected areas in each watershed. I want to retain the designation of each protected area as well (i.e., federal, state, etc.)

To do this I started to try to intersect the two shapefiles:

pas <- readOGR(dsn = ".", layer = "PADUS1_4Combined") # protected areas

huc <- readOGR(dsn = ".", layer = "HUC2_USA_PRJ") # hydrologic units (2 digit, USGS)

Checked the CRS

crs(huc)

crs(pas)

crs(huc)
CRS arguments:
+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0

crs(pas)
CRS arguments:
+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0

Intersection:

test <- gIntersection(pas, huc)

I get this error:

Error in RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, drop_lower_td, unaryUnion_if_byid_false, :
rgeos_crdMat2LinearRing: linearRing not created

I searched for this error and cannot find anything. I am dealing with a very large dataset here so I am wondering if this is the issue.

Is there a workaround using simple features (sf)?

r sp rgdal sf rgeos

share|improve this question

asked Apr 2 at 22:16

EcostriderEcostrider

162

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  It could be the size of the dataset (you've not said how big it is) or it could be a problem with the geometry of a feature - perhaps a polygon not valid or something. Try it with subsets. Also try sf and see what happens.
  
  – Spacedman
  Apr 3 at 7:12
  
  

  
  
  
  

add a comment | 

2

I am trying to do something relatively simple but with large shapefiles. I have a shapefile of protected areas in the USA:

and watershed boundaries:

I want to calculate the percentage and proportion of protected areas in each watershed. I want to retain the designation of each protected area as well (i.e., federal, state, etc.)

To do this I started to try to intersect the two shapefiles:

pas <- readOGR(dsn = ".", layer = "PADUS1_4Combined") # protected areas

huc <- readOGR(dsn = ".", layer = "HUC2_USA_PRJ") # hydrologic units (2 digit, USGS)

Checked the CRS

crs(huc)

crs(pas)

crs(huc)
CRS arguments:
+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0

crs(pas)
CRS arguments:
+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0

Intersection:

test <- gIntersection(pas, huc)

I get this error:

Error in RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, drop_lower_td, unaryUnion_if_byid_false, :
rgeos_crdMat2LinearRing: linearRing not created

I searched for this error and cannot find anything. I am dealing with a very large dataset here so I am wondering if this is the issue.

Is there a workaround using simple features (sf)?

r sp rgdal sf rgeos

share|improve this question

asked Apr 2 at 22:16

EcostriderEcostrider

162

• 
  
  
  

  
  

  
  
  
  
  
  

  
  
  It could be the size of the dataset (you've not said how big it is) or it could be a problem with the geometry of a feature - perhaps a polygon not valid or something. Try it with subsets. Also try sf and see what happens.
  
  – Spacedman
  Apr 3 at 7:12
  
  

  
  
  
  

add a comment | 

I am trying to do something relatively simple but with large shapefiles. I have a shapefile of protected areas in the USA:

and watershed boundaries:

I want to calculate the percentage and proportion of protected areas in each watershed. I want to retain the designation of each protected area as well (i.e., federal, state, etc.)

To do this I started to try to intersect the two shapefiles:

pas <- readOGR(dsn = ".", layer = "PADUS1_4Combined") # protected areas

huc <- readOGR(dsn = ".", layer = "HUC2_USA_PRJ") # hydrologic units (2 digit, USGS)

Checked the CRS

crs(huc)

crs(pas)

crs(huc)
CRS arguments:
+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0

crs(pas)
CRS arguments:
+proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0

Intersection:

test <- gIntersection(pas, huc)

I get this error:

Error in RGEOSBinTopoFunc(spgeom1, spgeom2, byid, id, drop_lower_td, unaryUnion_if_byid_false, :
rgeos_crdMat2LinearRing: linearRing not created

I searched for this error and cannot find anything. I am dealing with a very large dataset here so I am wondering if this is the issue.

Is there a workaround using simple features (sf)?

r sp rgdal sf rgeos

share|improve this question

asked Apr 2 at 22:16

EcostriderEcostrider

162

share|improve this question

asked Apr 2 at 22:16

EcostriderEcostrider

162

share|improve this question

asked Apr 2 at 22:16

EcostriderEcostrider

162

asked Apr 2 at 22:16

EcostriderEcostrider

162

asked Apr 2 at 22:16

EcostriderEcostrider

162