Access bands with snapista

development python
#1

Hi,

i’m trying to process Sentinel-1 GRD data with snapista.
Following are my added nodes:

g.add_node(
    operator=Operator("Read", file=filename, formatName="SENTINEL-1", bandNames='Band_IW_VH'),
    node_id="read",
)

g.add_node(
    operator=Operator("Apply-Orbit-File", continueOnFail="true"),
    node_id="apply-orbit-file",
    source="read",
)

g.add_node(
    operator=Operator("Remove-GRD-Border-Noise", borderLimit="2000", trimThreshold="0.2"),
    node_id="border-noise-removal",
    source="apply-orbit-file",
)

g.add_node(
    operator=Operator("ThermalNoiseRemoval"),
    node_id="thermal-noise-removal",
    source="border-noise-removal")

g.add_node(
    operator=Operator("Calibration"),
    node_id="calibration",
    source="border-noise-removal"
)

g.add_node(
    operator=Operator('Speckle-Filter'),
    node_id="speckle-filtering",
    source="calibration"
)

g.add_node(
    operator=Operator("Terrain-Correction", mapProjection=crs),
    node_id="terrain-correction",
    source="calibration"
)

g.add_node(
    operator=Operator("SetNoDataValue"),
    node_id="set-no-data",
    source="terrain-correction")

g.add_node(
    operator=Operator("Write", file="test.tif", formatName='GDAL-GTiff-WRITER'),
    source="set-no-data",
    node_id="write")

Specifically, i want to use the cross-polarization data but i’m not sure if it works like that. The band name is from the abstracted metadata in SNAP. The resulting Geotiff contains 2 bands while i’m expecting just 1 band with the backscatter intensity. Further processing is carried out with rasterio, numpy etc.

I would be very grateful for any advice!