Rendering Transformation can cause bad rendering on reprojection to 3995

I'm using the Jiffle rendering transformation to remove some pixels, and the result is then reprojected from 4326 to 3995.

The result of the reprojection with a style that uses the map algebra is different from the result of the reprojection without the map algebra. It looks like the resolution got lower.

Here is how it looks without the map algebra:

Here is how it looks with the map algebra:

In a multiband NetCDF file (attached) one band is used to decide whether or not render pixels from another band.
The Jiffle script is this:

A colormap is then applied; the full style is attached.
The layer preview image looks correct but when using different projections, artifacts come up.

To reproduce:

• Get the latest platform independent binary

• Add NetCDF plugin

• Add WPS plugin

• Create a NetCDF store using the attached .nc file

• Create a CoverageView with the following bands in this order:

  • sea_surface_temperature

  • quality_level

• Apply the attached style to it

• Issue a request for a tile in EPSG:3995 : http://localhost:8080/geoserver/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&FORMAT=image%2Fpng&TRANSPARENT=true&LAYERS=cite%3Asingle_file_coverage_view&STYLES=jiffle&WIDTH=512&HEIGHT=512&CRS=EPSG%3A3995&BBOX=-1128800%2C-1128800%2C10442400%2C10442400