```r
# all_dorset_sssi<-st_read(conn,query="select s.* from sssi s, dorset d where st_intersects(s.geom,d.geometry)")
# 
# hh<- st_read(conn,query="select s.* from sssi s where sssi_name ILIKE 'Christchurch%arbour%'")
# bound<-as(hh, "Spatial")
# 
# dsm<- pgGetRast(conn, "dsm2m",boundary=bound)
# dtm<- pgGetRast(conn, "dtm2m",boundary=bound)
# mapview(dsm) +mapview(hh)
```