# NOT RUN {
tr <- opq ("trentham australia") %>% osmdata_sf ()
coliban <- tr$osm_lines [which (tr$osm_lines$name == "Coliban River"),]
pts <- osm_points (tr, rownames (coliban)) # all points of river
waterfall <- pts [which (pts$waterway == "waterfall"),] # the waterfall point
# }
Run the code above in your browser using DataLab