Subset S3 WFR

Optical Toolbox
#1

I am trying to subset an image Sentinel 3 WFR using python and snappy with the following method:

def Crop(self, wkt):
if wkt != ‘nan’:
print (“Subset image with wkt”)
HashMap = jpy.get_type(‘java.util.HashMap’)
geom = WKTReader().read(wkt)
temp = self.product
parameters = self.HashMap()
parameters.put(“geoRegion”, geom)
self.product = GPF.createProduct(‘Subset’, parameters, temp)
else:
print(“Wkt was not defined. All the product will be used”)
return

However my result is cropped but not georeferenced and not with the correct resolution as it would be if i did subset from snap desktop.

Any ideas?

#2

How you define the subset looks correct.
How do you read the product and how do you write it? Which output format do you use?
The resolution should not change as long you don’t set the paramter subSampling[X/Y].