Spatial 2

setwd("J:/TX_CRIS and LADOTD/TX/SpecialCases")
library(data.table)
library(ggplot2)
dat0 <- fread("CRIS_2015_2019_LL_Year_Sev.csv")
#head(dat0)


#dat1 <- subset(dat0, Year==2018)
#dat2= head(dat1, 20000)

library(hexbin)
ggplot(dat0, aes(x=Latitude, y=Longitude)) + geom_hex(bins=100) +
  facet_wrap(~ Year, nrow = 2)+theme_bw(base_size=14)+
  scale_fill_gradientn(colours=c("#D3DDDC","blue"),name = "Frequency",na.value=NA)

library(hexbin)
ggplot(dat0, aes(x=Latitude, y=Longitude)) + geom_hex(bins=40) +
  facet_wrap(~ Year, nrow = 2)+theme_bw(base_size=14)+
  scale_fill_gradientn(colours=c("#D3DDDC","blue"),name = "Frequency",na.value=NA)

ggplot(dat0, aes(x=Latitude, y=Longitude)) + geom_hex(bins=25) +
  facet_wrap(~ Year, nrow = 2)+theme_bw(base_size=14)+
  scale_fill_gradientn(colours=c("#D3DDDC","blue"),name = "Frequency",na.value=NA)

ggplot(dat0, aes(x=Latitude, y=Longitude)) + geom_hex(bins=15) +
  facet_wrap(~ Year, nrow = 2)+theme_bw(base_size=14)+
  scale_fill_gradientn(colours=c("#D3DDDC","blue"),name = "Frequency",na.value=NA)

ggplot(dat0, aes(x=Latitude, y=Longitude)) + geom_hex(bins=10) +
  facet_wrap(~ Year, nrow = 2)+theme_bw(base_size=14)+
  scale_fill_gradientn(colours=c("#D3DDDC","blue"),name = "Frequency",na.value=NA)