```r
# query<- "drop table if exists calc_grass;create table calc_grass as select * from ph_v2_1 WHERE main_habit ILIKE '%calc%grass%'"
# dbSendQuery(conn,query)
# query<- "CREATE INDEX cg_ind ON calc_grass USING GIST (geom) "
# dbSendQuery(conn,query)
# query<-"select * from calc_grass"
# d<-st_read(conn,query=query)
# # library(mapview)
# mapview(d)
```