Weekend GIS Fun: Part 2

See part 1 for the relevant shape file and description of what’s going on! In the last post, we identified the map codes of interest to subset the shape file. Here’s some code for working with it in R. library(rgdal) bedrock <- readOGR(dsn=”~/Arrowhead_South_Bedrock_Geology”,”Arrowhead_South_Bedrock_Geology”) basalt = subset(bedrock,map_label==”Mnd” | map_label==”Mna”) writeOGR(basalt,dsn=”.”,layer=”basalt”,driver=”ESRI Shapefile”) So now we have a […]

Continue reading..Weekend GIS Fun: Part 2