sentinel-1/subset.py at master · MapTailor/sentinel-1 · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
import snappy
from snappy import ProductIO as product
from snappy import HashMap as hash
from snappy import GPF as gpf

import os, sys

gpf.getDefaultInstance().getOperatorSpiRegistry().loadOperatorSpis()
snappy.HashMap = snappy.jpy.get_type('java.util.HashMap')

path = "C:\\"

for folder in os.listdir(path):

    filename = folder + ".SAFE"

    dataset = os.path.join(path, folder, filename)
    dataset = dataset.format(sys.platform, os.sep)

    print dataset

    timestamp = folder.split("_")[5]
    date = timestamp[:8]

    parameters = hash()
    sources = hash()

    sentinel = product.readProduct(dataset + "//manifest.safe")

    WKTReader = snappy.jpy.get_type('com.vividsolutions.jts.io.WKTReader')

    wkt = "POLYGON((" \
          "14.1223022771676 36.116795674787561, " \
          "14.623741501140653 36.119420063224389, " \
          "14.625412269082306 35.766175579973442, " \
          "14.126199181086793 35.763584814742948, " \
          "14.1223022771676 36.116795674787561" \
          "))"
    geom = WKTReader().read(wkt)

    pols = ["VV", "VH"]

    for pol in pols:

        parameters.put('sourceBands', "Amplitude_" + pol)
        parameters.put('geoRegion', geom)
        parameters.put('copyMetadata', True)

        output = date + "_subset_" + pol
        checks = output + ".tif"

        writefile = os.path.join(path, folder, filename, output)
        checkfile = os.path.join(path, folder, filename, checks)

        if os.path.exists(checkfile) == True:
            pass

        else:

            item = gpf.createProduct("Subset", parameters, sentinel)
            product.writeProduct(item, writefile, 'GeoTIFF')