Intersection between geoJson and product

[tilmahutli]

Hello marpet,

I have a question. I need to obtained intersection zone between an area saved in a .geojson file versus a sentinel 1 product using snappy python.

The file .geojson has the polygon, how I can obtained first the image polygon? and second, how I can obtain the intersection between both images or zones?

Do you have any idea?

Thanks marpet

[marpet]

Geojson is not supported by by SNAP, but you can convert it to a WKT (https://mygeodata.cloud/converter/geojson-to-wkt).

You can use this WKT with the subset operator. Here is an example.

Python subset

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_ty…

Instead of specifying the WKT in the code you can read it from a file.