raster

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.

Context: this from @agila5: https://twitter.com/a_gilardi5/status/1487813588454678528

I have not had a chance to read the explanation and am not sure what we should say on the topic but guess it's worth a mention here. Any thoughts?

Is your feature request related to a problem?

Tests running on Ubuntu download the test data every time they run.

The URL is:

https://grass.osgeo.org/sampledata/north_carolina/nc_spm_full_v2alpha2.tar.gz

Describe the solution you'd like

GitHub Actions cache is used instead of a download inside the test script.

Describe alternatives you've considered

All these would

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

The docstring doesn't show it

Oops wrong tag.

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