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
