Fires in the UK and Croatia

library(giscourse)
library(dplyr)
conn<-sconnect11()

library(raster)
library(mapview)

query1<-"select * from nat_earth_countries_hi where name ilike '%roatia%'"
query2<-"select * from fires"


query<-sprintf("select f.* from (%s) c, (%s) f where st_intersects(f.geom,c.geometry)",query1,query2)

st_read(conn,query=query) -> croatia_fires
st_read(conn,query=query1) -> croatia

## type is 51

croatia_fires$area<-56594

lc<-get_lc(x="Zadar",dist=350)
land_cover<-data.frame(id=raster::extract(lc$lc,as(croatia_fires,"Spatial")))
land_cover<-base::merge(land_cover,lc$legend)
data.frame(table(land_cover$land_cover)) %>% DT::datatable()

xy<-st_coordinates(croatia_fires)
library(MASS)
k <- kde2d(xy[,1],xy[,2],h=.1,n=1000)
r <- raster(k)
r[r<0.1]<-NA
library(RColorBrewer)

pal<-RColorBrewer::brewer.pal(11,"Spectral")[11:1]
plot(croatia$geometry)
raster::plot(r,col=pal, add=TRUE) 

mapview(r, na.color="transparent") %>% giscourse::global_landcover_wms()