raster

Please respect the following template to report a new issue, otherwise your issue will be automatically closed.

Before submitting, search through existing open or closed issues reports and through blender.stackexchange.com to check whether the issue already exists.

If the issue concerns an operator that grab geodata from an external web service (Basemaps, OSM, SRTM), then it can be a network

The gdal/swig/python/samples/ogr_layer_algebra.py should be promoted to an official Python script

Steps:

• move it to gdal/swig/python/scripts
• create a new test script for it in autotest/pyscripts/test_ogr_layer_algebra.py
• write a documentation page in doc/source/programs/ogr_layer_algebra (and link it to the index)

rasterio.warp.reproject takes a kwargs parameter, which is passed to the "transformation function". It is not by any means clear what this means. What transformation function? What are the allowed values for kwargs? Can this please be expanded on? Thanks.

Describe the bug
When recompressing data from one compression to another, everything works but you still get

[1]   Exit 1                  r.compress SVI_2000_log

To Reproduce

r.compress elevation
echo $?
1

Additional context
Looks like the problem is here:
https://github.com/OSGeo/grass/blob/master/raster/r.compress/main.c#L99
The logic should be rev

Describe the bug
Using ep.plot_rgb() with stretch=True fails when plotting a masked array when the original array contains NaNs. Plotting returns a black plot and gives the following warning "UserWarning: One or more intensity levels are NaN. Rescaling will broadcast NaN to the full image. Provide intensity levels yourself to avoid this. E.g. with np.nanmin(image), np.nanmax(image)."

I th

link: https://github.com/milesgranger/pyrus-cramjam

right now we have to install brotli dependency and use default middleware which is not customizable (ref developmentseed/titiler#268)

I think we could create a simple middleware which support zlib and brotli using pyrus-cramjam and also add a path filter as we've done in https://github.com/developmentseed/ti

https://carpentries-incubator.github.io/geospatial-python/ is a good course that uses rioxarray. Could add a link to it in https://corteva.github.io/rioxarray/stable/getting_started/getting_started.html

In some ways fasterize does diverge from rasterize (e.g., #13). These should be tersely described in the README/vignette and function documentation.