# Import the espatools module for raster I/O
import espatools
# Use matplotlib for simple plotting
import matplotlib.pyplot as plt

xmlfile='/Users/bane/Documents/OpenGeoVis/Data/data-testing/Readers/Landsat/Golden/LC080340322018062701T1-SC20180922164434/LC08_L1TP_034032_20180627_20180704_01_T1.xml'
reader = espatools.RasterSetReader(filename=xmlfile)

raster = reader.read()#cast=True, allowed=['sr_aerosol', 'sr_band7'])
bands = raster.bands
bands

for n,b in bands.items():
    print(b.data.dtype)

b = bands.get('sr_band3')

plt.figure(figsize=(10,10))
plt.imshow(b.data, vmin=0., vmax=2000)
plt.title(b.name)
plt.savefig('aer.png')
plt.show()

color = raster.get_rgb('false_a')

plt.figure(figsize=(10,10))
plt.imshow(color)
plt.title('RGB Color')
plt.savefig('RGB.png')
plt.show()

raster.global_metadata.satellite

raster.nsamps

raster.nlines

color.shape

raster.global_metadata.projection_information.corner_point[0].x

b.valid_range.min

b.fill_value

color.shape