"""
/*********************************************************************************/
* The MIT License (MIT) *
* *
* Copyright (c) 2014 EOX IT Services GmbH *
* *
* Permission is hereby granted, free of charge, to any person obtaining a copy *
* of this software and associated documentation files (the "Software"), to deal *
* in the Software without restriction, including without limitation the rights *
* to use, copy, modify, merge, publish, distribute, sublicense, and/or sell *
* copies of the Software, and to permit persons to whom the Software is *
* furnished to do so, subject to the following conditions: *
* *
* The above copyright notice and this permission notice shall be included in *
* all copies or substantial portions of the Software. *
* *
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR *
* IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, *
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE *
* AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER *
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, *
* OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE *
* SOFTWARE. *
* *
*********************************************************************************/
"""
#!/usr/bin/env python
import sys
import os
import re
import logging
import xml.etree.ElementTree as ET
import geojson
import shapely.wkt
from zipfile import ZipFile
try:
from osgeo import gdal
except:
sys.exit('ERROR: cannot find GDAL/OGR modules, install gdal with python bindings')
log = logging.getLogger(__name__)
def get_manifest_zip_path(zipfile_path):
with ZipFile(zipfile_path) as z:
manifest_zip_path = None
for p in z.namelist():
if p.endswith('manifest.safe'):
manifest_zip_path = p
break
return "/vsizip/" + zipfile_path.rstrip('/') + '/' + manifest_zip_path
def extract_file_from_zip(endswith_str, zipfile_path, out_dir=None):
with ZipFile(zipfile_path) as z:
file_zip_path = None
for p in z.namelist():
if p.endswith(endswith_str):
file_zip_path = p
break
if file_zip_path is None:
return None
# extract it
if out_dir is None:
targetpath = z.extract(file_zip_path)
else:
targetpath = z.extract(file_zip_path, path=out_dir)
return targetpath
def extract_manifest_from_zip(zipfile_path, out_dir=None):
with ZipFile(zipfile_path) as z:
# look for manifest
manifest_zip_path = None
for p in z.namelist():
if p.endswith('manifest.safe'):
manifest_zip_path = p
break
if manifest_zip_path is None:
return None
# extract it
if out_dir is None:
targetpath = z.extract(manifest_zip_path)
else:
targetpath = z.extract(manifest_zip_path, path=out_dir)
return targetpath
class S1GDALReader:
def __init__(self, sentinel1_product_zip_path):
self.product_zip_path = sentinel1_product_zip_path
sentinel1_product_dir = os.path.dirname(sentinel1_product_zip_path)
sentinel1_product_zipname = os.path.basename(sentinel1_product_zip_path)
self.product_dir = sentinel1_product_dir
self.sentinel1_product_zipname = sentinel1_product_zipname
self.granule_identifier, _ = os.path.splitext(sentinel1_product_zipname)
manifest_path = extract_manifest_from_zip(sentinel1_product_zip_path)
self.manifest_tree = ET.parse(manifest_path)
try:
os.remove(manifest_path)
os.rmdir(os.path.dirname(manifest_path))
except:
log.warn("Cannot cleanup manifest directory")
#sentinel1_safe_pkg_path = "/vsizip/{}/{}.SAFE/manifest.safe".format(self.product_zip_path, self.granule_identifier)
#self.safe_package_path = sentinel1_safe_pkg_path
#self.datastore = gdal.Open(sentinel1_safe_pkg_path)
def get_metadata(self):
manifest_zip_path = get_manifest_zip_path(self.product_zip_path)
datastore = gdal.Open(manifest_zip_path)
metadata_dict = datastore.GetMetadata()
metadata_dict['NAME'] = self.granule_identifier
startTime = metadata_dict['ACQUISITION_START_TIME']
endTime = metadata_dict['ACQUISITION_STOP_TIME']
# round to milliseconds
m = re.search("(.*)(\d{3})(\d{3})$", startTime)
if m:
startTime = m.groups()[0] + m.groups()[1]
m = re.search("(.*)(\d{3})(\d{3})$", endTime)
if m:
endTime = m.groups()[0] + m.groups()[1]
metadata_dict['ACQUISITION_START_TIME'] = startTime + 'Z'
metadata_dict['ACQUISITION_STOP_TIME'] = endTime + 'Z'
return metadata_dict
def get_footprint(self):
GML_NS = "{http://www.opengis.net/gml}"
gml_coordinates = ""
for el in self.manifest_tree.iter(GML_NS + "coordinates"):
gml_coordinates = el.text
wkt_coordinates = gml_coordinates.replace(",", ";")
wkt_coordinates = wkt_coordinates.replace(" ", ",")
wkt_coordinates = wkt_coordinates.replace(";", " ")
wkt_coordinates = wkt_coordinates + "," + wkt_coordinates.split(",")[0]
s = "POLYGON ((" + wkt_coordinates + "))"
log.info("stringa s = " + s)
g1 = shapely.wkt.loads(s)
g2 = geojson.Feature(geometry=g1, properties={})
print type(g2.geometry)
return g2.geometry
def get_quicklook(self):
return extract_file_from_zip("quick-look.png", self.product_zip_path)
def get_preview_image(self):
return os.path.join(self.safe_package_path, "preview", "quick-look.png")
