---
title: "Rhinolophus, Palm civet and site map"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = FALSE)
library(giscourse)
con<-sconnect11("clima2")
query<-"select genus,binomial,geometry from iucnmammals where genus = 'Rhinolophus'"
library(sf)
d<-st_read(con,query=query)
library(mapview)
unique(d$binomial)
ch<-c("Rhinolophus ferrumequinum","Rhinolophus sinicus","Rhinolophus affinis")
d<-d[d$binomial %in% ch,]

query<-"select genus,binomial,geometry from iucnmammals where binomial = 'Paradoxurus hermaphroditus'"
d2<-st_read(con,query=query)
d<-rbind(d,d2)
query<-"select genus,binomial,geometry from iucnmammals where binomial = 'Manis javanica'"
d2<-st_read(con,query=query)
d<-rbind(d,d2)

mp<-mapview(d,zcol="binomial",burst=TRUE, legend=FALSE)
dd<-data.frame(name=c("Wuhan", "Tongguan"),X=c(114.3,110.24),Y=c(30.6,34.54))
Sites=st_as_sf(dd,coords=c("X","Y"),crs=4326)
mp<-mp+mapview(Sites, legend=FALSE)
save(mp,file="rmap.rda")
```

```{r}
mp
```

