```r
hh_habitat<-"select main_habit mhabitat, 
s1habtype,s2habtype,s3habtype habitat,
st_intersection(h.geom,s.geom) geom from 
(select * from local_reserves where lnr_name = 'Hengistbury Head') s,
ph_v2_1 h 
where st_intersects(h.geom,s.geom)"
habitats<-st_read(conn, query=hh_habitat)
qtm(habitats, "mhabitat") +tm_style("classic")
```

<img src="OS_GIS_files/figure-html/unnamed-chunk-87-1.png" width="672" />