I am currently reprojecting SLSTR and OLCI data and exporting them to geotiffs but running into two issues while re-projecting.
I would like to use the bicubic resampling instead of nearest neighbor, however when doing so NoData values (-327680 for the S1_radiacne for example) are being interpolated to nearby reprojected pixels causing the data to not get masked out and become extremely negative.
Is there any alternative to just using nearest neighbor that will yield nicely resampled image while somehow preventing pixels neighboring NoData values from being interpolated?
Also if I use NaN as the NoData value I can interpolate and everything seems to work as intended, but I do not want to do this for bands which are integer based such as S1-6 radiance since the integer values get converted to float.
I would like the preserve the datatype of the bands during the process of reprojecting -> writing as a geotiff. However the geotiff output seems to be float32 by default. Is there some implicit casting happening during reproject, and can I prevent this from happening?