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()