Hi All:
I wrote a script to export and reproject NDVI and LST from S3 SLSTR L2, but the NDVI band is exported with wrong values: the only values for all pixels is -32768.000000. The LST band is exported with correct values.
The workflow I use is: subset the S3 images, reproject the result in EPSG 3857 and save the reprojected dataset in Geotiff format.
Does anyone have similar export errors?
The export from SNAP software returns correct values for NDVI.
Script:
subset_op = snappy.jpy.get_type(‘org.esa.snap.core.gpf.common.SubsetOp’)()
subset_op.setSourceProduct(product_s3)
subset_op.setGeoRegion(WKTReader().read(overlap_WKT))
subset_op.setBandNames([“NDVI”])
target_op = subset_op.getTargetProduct()
reproject_op = snappy.jpy.get_type("org.esa.snap.core.gpf.common.reproject.ReprojectionOp")()
reproject_op.setSourceProduct(target_op)
reproject_op.setParameter("crs","EPSG:3857")
reproject_op.setParameter("resampling", "Nearest")
reproject_p = reproject_op.getTargetProduct()
ProductIO.writeProduct(reproject_p, r"path_to_file.tif", "GeoTIFF-BigTIFF")
Thank you