Guest User

Untitled

a guest
Nov 19th, 2018
94
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.66 KB | None | 0 0
  1. output = path + folder + "\"
  2.  
  3. timestamp = folder.split("_")[4]
  4.  
  5. date = timestamp[:8]
  6.  
  7. sentinel_1 = ProductIO.readProduct(output)
  8.  
  9. print sentinel_1
  10.  
  11. pols = ['VH','VV']
  12.  
  13. for p in pols:
  14.  
  15. polarization = p
  16.  
  17. wkt = 'POLYGON((17.06138 98.4766, 17.07500 98.55027, 17.10500 98.46833, 17.11833 98.54166))'
  18.  
  19. geom = WKTReader().read(wkt)
  20.  
  21. parameters = HashMap()
  22.  
  23. parameters.put('geoRegion', geom)
  24.  
  25. parameters.put('outputImageScaleInDb', False)
  26.  
  27. subset = output+date+"_subset_" + polarization
  28.  
  29. target_1 = GPF.createProduct("Subset", parameters, sentinel_1)
  30.  
  31. ProductIO.writeProduct(target_1, subset, 'BEAM-DIMAP')
  32.  
  33. geom = WKTReader().read(wkt)
Add Comment
Please, Sign In to add comment