WKTReader = beampy.jpy.get_type('com.vividsolutions.jts.io.WKTReader')
if len(sys.argv) != 3:
print("usage: %s <file> <geometry-wkt>" % sys.argv[0])
print(" %s ./TEST.N1 \"POLYGON((15.786082 45.30223, 11.798364 46.118263, 10.878688 43.61961, 14.722727"
"42.85818, 15.786082 45.30223))\"" % sys.argv[0])
sys.exit(1)
file = sys.argv[1]
wkt = sys.argv[2]
# Uncomment if you receive errors of type com.sun.media.jai.util.ServiceConfigurationError, see
# http://www.brockmann-consult.de/beam-jira/browse/BEAM-1699
#beampy.SystemUtils.init3rdPartyLibs(None) # Initialise BEAM's third party Java libraries JAI and GeoTools.
geom = WKTReader().read(wkt)
print("Reading...")
product = ProductIO.readProduct(file)
op = SubsetOp()
op.setSourceProduct(product)
op.setGeoRegion(geom)
sub_product = op.getTargetProduct()
print("Writing...")
ProductIO.writeProduct(sub_product, "beampy_subset_output.dim", "BEAM-DIMAP")
print("Done.")
if len(sys.argv) != 3:
print("usage: %s <file> <geometry-wkt>" % sys.argv[0])
print(
" %s ./TEST.N1 \"POLYGON((15.786082 45.30223, 11.798364 46.118263, 10.878688 43.61961, 14.722727"
"42.85818, 15.786082 45.30223))\"" % sys.argv[0])
sys.exit(1)
file = sys.argv[1]
wkt = sys.argv[2]
# Uncomment if you receive errors of type com.sun.media.jai.util.ServiceConfigurationError, see
# http://www.brockmann-consult.de/beam-jira/browse/BEAM-1699
#beampy.SystemUtils.init3rdPartyLibs(None) # Initialise BEAM's third party Java libraries JAI and GeoTools.
geom = WKTReader().read(wkt)
print("Reading...")
product = ProductIO.readProduct(file)
op = SubsetOp()
op.setSourceProduct(product)
op.setGeoRegion(geom)
sub_product = op.getTargetProduct()
print("Writing...")
ProductIO.writeProduct(sub_product, "beampy_subset_output.dim", "BEAM-DIMAP")
print("Done.")
targetBand2.name = 'band_2'
targetBand2.type = 'float32'
targetBand2.expression = '(radiance_9 - radiance_6) / (radiance_9 + radiance_6)'
targetBands = jpy.array('org.esa.beam.gpf.operators.standard.BandMathsOp$BandDescriptor', 2)
targetBands[0] = targetBand1
targetBands[1] = targetBand2
parameters = HashMap()
parameters.put('targetBands', targetBands)
result = GPF.createProduct('BandMaths', parameters, product)
print("Writing...")
ProductIO.writeProduct(result, 'beampy_bmaths_output.dim', 'BEAM-DIMAP')
print("Done.")
"""
Please note: the next major version of beampy/jpy will be more pythonic in the sense that implicit data type
conversions are performed. The 'parameters' from above variable could then be given as a Python dict object:
parameters = {
'targetBands': [
{
'name': 'band_1',
'type': 'float32',
'expression': '(radiance_10 - radiance_7) / (radiance_10 + radiance_7)'
},
