Dear @marpet,
I have recently upgraded from SNAP 6.0 to SNAP 7.0, along with all the latest modules (module_list.txt)
Unfortunately I can no longer subset Sentinel-5P NRTI files (nc format). The Subset operator returns an error stating that the product has no geocoding:
RuntimeError: java.lang.IllegalArgumentException: Product without geo-coding.
I use snappy to process the products, here is an snippet:
from snappy import GPF
from snappy import ProductIO
from snappy import HashMap
# open product
fname = 'S5P_NRTI_L2__SO2____20181118T112353_20181118T112853_05691_01_010102_20181118T122344.nc'
p = ProductIO.readProduct(fname)
# subset product
parameters = HashMap()
parameters.put('copyMetadata', 'true')
parameters.put('geoRegion', 'POLYGON ((33.748 2.977, 24.752 2.977, 24.733 -6.007, 33.767 -6.007, 33.748 2.977))')
#parameters.put('sourceBands', 'sulfurdioxide_total_vertical_column_1km')
p = GPF.createProduct('Subset', parameters, p)
It seems that the band are no longer recognized/accessed properly: with SNAP 6 the recognized bands had names like “sulfurdioxide_total_vertical_column_1km”, whereas with SNAP 7 the same band is recognized as “sulfurdioxide_total_vertical_column_1km_PRODUCT/SUPPORT_DATA/DETAILED_RESULTS”.
Is there an easy fix to this? I use snappy to manage an operational monitoring framework using Sentinel-1, Sentinel-2, and Sentinel-5P products, so I would deeply like to stick with SNAP rather than try new toolboxes (e.g., BEAT).
Thanks in advance for any advice or suggestion!