In [5]:
import time
from shapely.geometry import box
from gbdxtools import Interface
gbdx = Interface()
bbox = [-109.84, 43.19, -109.59, 43.34]
aoi = box(*bbox).wkt
query = 'item_type:Landsat8 AND attributes.cloudCover_int:0'
res = gbdx.vectors.query(aoi, query)
print len(res), 'Images'
_id = res[0]['properties']['attributes']['catalogID']
print _id
In [6]:
from gbdxtools import LandsatImage
img = LandsatImage(_id, bbox=bbox)
img.plot(w=10,h=10, bands=[3,2,1])
In [ ]: