code: def pointquery(features, raster, band, indent, nodata,
                 interpolate, property_name, sequence, use_rs):
      
  
      results = gen_point_query(
          features,
          raster,
          band=band,
          nodata=nodata,
          interpolate=interpolate,
          property_name=property_name,
          geojson_out=True)
  
      if sequence:
          for feature in results:
              if use_rs:
                  click.echo(b'\x1e', nl=False)
              click.echo(json.dumps(feature))
      else:
          click.echo(json.dumps(
              {'type': 'FeatureCollection',
               'features': list(results)}))
docs: Queries the raster values at the points of the input GeoJSON Features.
      The raster values are added to the features properties and output as GeoJSON
      Feature Collection.
  
      If the Features are Points, the point geometery is used.
      For other Feauture types, all of the verticies of the geometry will be queried.
      For example, you can provide a linestring and get the profile along the line
      if the verticies are spaced properly.
  
      You can use either bilinear (default) or nearest neighbor interpolation.