code: def dump(dataset, query, out_file, bounds):
      
      table = bcdata.validate_name(dataset)
      data = bcdata.get_data(table, query=query, bounds=bounds)
      if out_file:
          with open(out_file, "w") as f:
              json.dump(data.json(), f)
      else:
          sink = click.get_text_stream("stdout")
          sink.write(json.dumps(data))
docs: Write DataBC features to stdout as GeoJSON feature collection.
  
      \b
        $ bcdata dump bc-airports
        $ bcdata dump bc-airports --query "AIRPORT_NAME='Victoria Harbour (Shoal Point) Heliport'"
        $ bcdata dump bc-airports --bounds xmin ymin xmax ymax
  
      The values of --bounds must be in BC Albers.
  
       It can also be combined to read bounds of a feature dataset using Fiona:
      \b
        $ bcdata dump bc-airports --bounds $(fio info aoi.shp --bounds)