import os
import glob
from IPython.core.display import Image
from pymodis import downmodis

# destination foldert
dest = "/tmp"
# tiles to download
tiles = "h18v04,h19v04"
# starting day
day = "2014-08-14"
# number of day to download
delta = 1

modisDown = downmodis.downModis(destinationFolder=dest, tiles=tiles, today=day, delta=delta)

modisDown.connect()

modisDown.downloadsAllDay()

# create the list of files to use
files = glob.glob(os.path.join(dest, 'MOD11A1.A2014*.hdf'))
print files

from topng import hdf2png
pngfile = hdf2png(files[0])
Image(filename=pngfile)

from pymodis import parsemodis

modisParse = parsemodis.parseModis(files[0])

# bounding box of the tiles
modisParse.retBoundary()

# quality statistics
modisParse.retMeasure()

modisParse.retRangeTime()

parsemodis.parseModis?

modisMultiParse = parsemodis.parseModisMulti(files)

modisMultiParse.valBound()
modisMultiParse.boundary

modisMultiParse.writexml(os.path.join(dest, 'modismultiparse.xml'))

f = open(os.path.join(dest, 'modismultiparse.xml'))
lines = f.readlines()
p = [l.strip() for l in lines]
f.close()
print "\n".join(p)

from pymodis.convertmodis_gdal import createMosaicGDAL
# [daily temp, quality for daily, not used, not used, nightly temp, quality for nightly]
subset = [1,1,0,0,1,1]
output_pref = os.path.join(dest, 'MOD11A1.A2014226.mosaic')
output_tif = os.path.join(dest, 'MOD11A1.A2014226.mosaic.tif')

# the first parameter is a list with the original tiles,
# the second one is a list with the the subset to process,
# the last is the output format, in this case GeoTiff
mosaic = createMosaicGDAL(files, subset, 'GTiff')

mosaic.run(output_tif)

mosaic.write_vrt(output_pref)

from pymodis.convertmodis_gdal import convertModisGDAL
output_pref = os.path.join(dest, 'MOD11A1.A2014226.h18v04')

convertsingle = convertModisGDAL(hdfname=files[0], prefix=output_pref, subset=subset, res=1000, epsg=3035)
convertsingle.run()

vrtfiles = glob.glob(os.path.join(dest, 'MOD11A1.A2014*.vrt'))
print vrtfiles

for f in vrtfiles:
    base = os.path.basename(f).replace('.vrt', '_vrt')
    output = os.path.join(dest, base)
    convertsingle = convertModisGDAL(hdfname=f, prefix=output, subset= [1,1,1,1], res=1000, epsg=3035, vrt=True)
    convertsingle.run_vrt_separated()

vrttiffiles = glob.glob(os.path.join(dest, 'MOD11A1.A2014*_vrt.tif'))
print vrttiffiles

from topng import tif2png
pngfile = tif2png(vrttiffiles[1])
Image(filename=pngfile)