def search(self):
""" Search for Landsat-8 scenes """
# Convert AOI to GeoJSON
aoi_file = tempdir() + '/aoi_geojson.geojson'
Module('v.out.ogr',
overwrite=True,
input=options['map'],
format='GeoJSON',
output=aoi_file)
# Reproject to espg 4326
aoi_file_re = tempdir() + '/aoi_reprojected.geojson'
os.system("ogr2ogr -t_srs epsg:4326 " + aoi_file_re + " " + aoi_file)
# Search for scenes
with open(aoi_file_re) as f:
aoi = json.dumps(json.load(f))
search = Search(date_from=options['date_from'], date_to=options['date_to'], satellite_name='Landsat-8',
intersects=aoi, cloud_from=0, cloud_to=options['clouds'])
self._scenes = search.scenes()
os.remove(aoi_file)
def test_search(self):
""" Perform simple query """
with open(os.path.join(self.path, 'aoi1.geojson')) as f:
aoi = json.dumps(json.load(f))
search = Search(datetime='2017-01-05', intersects=aoi)
self.assertEqual(search.found(), 1)
scenes = search.scenes()
self.assertTrue(isinstance(scenes[0], Scene))
def query_satapi(aoi, clip=False, **kwargs):
""" Search sat-api using sat-search """
if clip:
aoi = get_coastline(aoi)
search = Search(intersects=json.dumps(aoi), **kwargs)
return search.scenes()
