Vjylku

Getting extremely large arrows with tikzcd

Why were 5.25" floppy drives cheaper than 8"?

Is it possible to create a QR code using text?

Why was Sir Cadogan fired?

Knowledge-based authentication using Domain-driven Design in C#

What is the most common color to indicate the input-field is disabled?

How do conventional missiles fly?

Is there a hemisphere-neutral way of specifying a season?

My ex-girlfriend uses my Apple ID to login to her iPad, do I have to give her my Apple ID password to reset it?

What does the same-ish mean?

Does int main() need a declaration on C++?

Blending or harmonizing

Can compressed videos be decoded back to their uncompresed original format?

How does a dynamic QR code work?

Notepad++ delete until colon for every line with replace all

When handwriting 黄 (huáng; yellow) is it incorrect to have a disconnected 草 (cǎo; grass) radical on top?

Do Iron Man suits sport waste management systems?

OP Amp not amplifying audio signal

Does the Idaho Potato Commission associate potato skins with healthy eating?

Are British MPs missing the point, with these 'Indicative Votes'?

Mathematica command that allows it to read my intentions

Can I hook these wires up to find the connection to a dead outlet?

Is this draw by repetition?

How badly should I try to prevent a user from XSSing themselves?

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?

0

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

asked 8 mins ago

EcostriderEcostrider

62

add a comment | 

0

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

asked 8 mins ago

EcostriderEcostrider

62

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

asked 8 mins ago

EcostriderEcostrider

62

share

asked 8 mins ago

EcostriderEcostrider

62

share

asked 8 mins ago

EcostriderEcostrider

62

asked 8 mins ago

EcostriderEcostrider

62

asked 8 mins ago

EcostriderEcostrider

62