I am trying to subset a S2 image and i get the following error
Traceback (most recent call last):
File “C:\Python34\doriforikes\test.py”, line 30, in
ProductIO.writeProduct(sub_product,‘subset_output.dim’,‘BEAM-DIMAP’)
RuntimeError: java.lang.RuntimeException: Cannot construct DataBuffer.
The code i have used is
import snappy
import numpy
import math
from snappy import ProductIO
from snappy import GPF
from snappy import jpy
from snappy import HashMap
SubsetOp = snappy.jpy.get_type(‘org.esa.snap.core.gpf.common.SubsetOp’)
WKTReader = snappy.jpy.get_type(‘com.vividsolutions.jts.io.WKTReader’)
HashMap = snappy.jpy.get_type(‘java.util.HashMap’)
snappy.GPF.getDefaultInstance().getOperatorSpiRegistry().loadOperatorSpis()
product=ProductIO.readProduct(‘C:\Program Files\snap\S2A_MSIL2A_20170910T090021_N0205_R007_T35SMD_20170910T090725.SAFE\MTD_MSIL2A.xml’)
HashMap = jpy.get_type(‘java.util.HashMap’)
parameters = HashMap()
parameters.put(‘targetResolution’,10)
result = GPF.createProduct(‘Resample’, parameters, product)
WKTReader = snappy.jpy.get_type(‘com.vividsolutions.jts.io.WKTReader’)
wkt =‘POLYGON ((26.060354643380858 39.226140746346694,26.318622752332 39.22793788561884,26.32017059429704 39.06723050504229,26.062488995951703 39.065443565792464,26.060354643380858 39.226140746346694))’
geometry = WKTReader().read(wkt)
op = SubsetOp()
op.setSourceProduct(result)
op.setGeoRegion(geometry)
sub_product = op.getTargetProduct()
ProductIO.writeProduct(sub_product,‘subset_output.dim’,‘BEAM-DIMAP’)
Any ideas?