```r
# library(rpostgis)
# arne<-st_read(conn,query="select * from sssi where sssi_name = 'Arne'")
# bound<-as(arne, "Spatial")
# dsm<- pgGetRast(conn, "dsm2m",boundary=bound)
# dtm<- pgGetRast(conn, "dtm2m",boundary=bound)
# mapview(dsm) +mapview(arne)
# 
```