require(rgdal) require(raster) require(ggplot2) require(dplyr) require(rasterVis) setwd("C:/DATA/Blog/raster/NLCD2011_LC_N39W072") li <- raster("NLCD2011_LC_N39W072.tif") plot(li)
shape <- readOGR(dsn = "C:/DATA/Blog/raster/lishape/government_units_NRCSCNTY_ny_3451580_01/government_units", layer = "county_nrcs_a_ny") suffolk <- subset(shape, COUNTYNAME == "Suffolk") plot(suffolk)
suffolk_c <- crop(li, extent(suffolk)) suffolk_c2 <- mask(suffolk_c, suffolk)dist <- data.frame(freq(suffolk_c2))
g <- ggplot(dist, aes(x=X.1, y=count2, fill=X.1)) + geom_bar(stat = "identity") + coord_flip()
Advertisements