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.

Fix one or more of the issues reported by Flake8.

This issue is not something which will be fixed by one PR. Instead, it is servers as a pointer to the list of ignored warnings in the source code. Pick one or more issues or files, fix them, and submit a PR.

Getting the existing issues to fix

1. Go to the source code.
2. Find one of the .flake8 files in the source tree.
3. See the `

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

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.