notebook.community

Edit and run



In [1]:

    
%pylab notebook
from charistools.readers import read_tile
from charistools.readers import ModisTileCube
import glob
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
import re
pd.options.display.max_rows = 200









    



/Users/mj/miniconda3/envs/charis/lib/python2.7/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.
  warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')






    



Populating the interactive namespace from numpy and matplotlib



In [6]:

    
%cd /Users/mj/projects/CHARIS/basins/basin_MODIS_tiles/
%ls
mask_filename = '/Users/mj/projects/CHARIS/basins/basin_MODIS_tiles/GA_Langtang_at_Kyanjin.basin_mask.h25v06.tif'
#mask_filename = '/Users/mj/projects/CHARIS/basins/basin_MODIS_tiles/IN_Hunza_at_DainyorBridge.basin_mask.h23v05.tif'
mask_filename









    



/Users/mj/projects/CHARIS/basins/basin_MODIS_tiles
00notes.txt
AM_Vakhsh_at_Komsomolabad.basin_mask.h23v05.tif
AM_Vakhsh_at_Komsomolabad.basin_mask.h23v05.tif.aux.xml
AM_Vakhsh_at_Komsomolabad.basin_mask.h23v05.tif_SJSK_20160225
AM_Zeravshan_at_Dupuli.basin_mask.h23v05.tif
GA_Arun_at_Turkeghat.basin_mask.h25v06.tfw
GA_Arun_at_Turkeghat.basin_mask.h25v06.tif
GA_Ganga_at_Rishikesh.basin_mask.h24v05.tif
GA_Ganga_at_Rishikesh.basin_mask.h24v06.tif
GA_Karnali_at_Benighat.basin_mask.h24v05.tif
GA_Karnali_at_Benighat.basin_mask.h25v05.tif
GA_Karnali_at_Benighat.basin_mask.h25v06.tif
GA_Langtang_at_Kyanjin.basin_mask.h25v06.tif
GA_Narayani_at_Devghat.basin_mask.h25v06.tfw
GA_Narayani_at_Devghat.basin_mask.h25v06.tif
GA_SaptaKosi_at_Chatara.basin_mask.h25v06.tif
GA_Tulsi_at_Betrawati.basin_mask.h25v06.tfw
GA_Tulsi_at_Betrawati.basin_mask.h25v06.tif
IN_Chenab_at_Akhnoor.basin_mask.h24v05.tfw
IN_Chenab_at_Akhnoor.basin_mask.h24v05.tif
IN_Chenab_at_Marala.basin_mask.h24v05.tif
IN_Gilgit_at_AlamBridge.basin_mask.h23v05.tfw
IN_Gilgit_at_AlamBridge.basin_mask.h23v05.tif
IN_Gilgit_at_AlamBridge.basin_mask.h24v05.tfw
IN_Gilgit_at_AlamBridge.basin_mask.h24v05.tif
IN_Gilgit_at_Gilgit.basin_mask.h23v05.tfw
IN_Gilgit_at_Gilgit.basin_mask.h23v05.tif
IN_Gilgit_at_Gilgit.basin_mask.h24v05.tfw
IN_Gilgit_at_Gilgit.basin_mask.h24v05.tif
IN_Hunza_at_DainyorBridge.basin_mask.h23v05.tfw
IN_Hunza_at_DainyorBridge.basin_mask.h23v05.tif
IN_Hunza_at_DainyorBridge.basin_mask.h24v05.tfw
IN_Hunza_at_DainyorBridge.basin_mask.h24v05.tif
IN_Jhelum_at_Chinari.basin_mask.h24v05.tif
IN_UpperIndus_at_Besham.basin_mask.h23v05.tfw
IN_UpperIndus_at_Besham.basin_mask.h23v05.tif
IN_UpperIndus_at_Besham.basin_mask.h24v05.tfw
IN_UpperIndus_at_Besham.basin_mask.h24v05.tif
IN_UpperIndus_at_Besham.basin_mask.h24v05.tif.aux.xml
rastert_langtang_basin_mask_rastert_watrshsb_sinu_h25v06.tif






    Out[6]:





'/Users/mj/projects/CHARIS/basins/basin_MODIS_tiles/GA_Langtang_at_Kyanjin.basin_mask.h25v06.tif'



In [ ]:

    
mask = read_tile(mask_filename)



In [ ]:

    
%cd /Users/mj/projects/CHARIS/snow_cover/modice.v0.4/min05yr_nc
#%ls
modice_filename = '/Users/mj/projects/CHARIS/snow_cover/modice.v0.4/min05yr_nc/MODICE.v0.4.h25v06.3strike.min05yr.mask.nc'
modice_filename



In [ ]:

    
modice = read_tile(modice_filename, 'modice_min_year_mask')



In [ ]:

    
print(modice.shape)
print(np.amin(modice),np.amax(modice))



In [ ]:

    
def display_basin_mask(file):
    data = read_tile(filename=file)
    fig, ax = plt.subplots()
    print(np.amin(data), np.amax(data))
    ax.imshow(data, cmap='Greys', vmin=np.amin(data), vmax=np.amax(data), interpolation='None')
    ax.set_title(file)
    plt.axis('off')



In [ ]:

    
display_basin_mask(filename)



In [ ]:

    
%cd /Users/mj/projects/CHARIS/forcing_data/Downscaled/ERA-Interim/yearly/v3/h25v06
%ls



In [ ]:

    
temp_filename = "ERA_Interim_downscale_uncorrected_tsurf.v0.3.h25v06_2001.h5"



In [ ]:

    
from charistools.readers import ModisTileCube



In [ ]:

    
cube = ModisTileCube(temp_filename, varname='tsurf')



In [ ]:

    
d = cube.read(doy=1)



In [ ]:

    
d.shape



In [ ]:

    
np.amin(d), np.amax(d)



In [ ]:

    
fig, ax = plt.subplots()
ax.imshow(d * mask2, cmap='Greys', vmin=np.amin(d), vmax=np.amax(d), interpolation='None')
ax.set_title(temp_filename)
plt.axis('off')



In [3]:

    
mask = read_tile(mask_filename)



In [4]:

    
print(mask.shape)
print(np.amin(mask), np.amax(mask))









    



(2400, 2400)
(-2147483648, 1)



In [5]:

    
fig, ax = plt.subplots()
ax.imshow(mask, cmap='Greys', vmin=np.amin(mask), vmax=np.amax(mask), interpolation='None')
ax.set_title(filename)
plt.axis('off')









    














    











    



---------------------------------------------------------------------------
NameError                                 Traceback (most recent call last)
<ipython-input-5-8d6acdc0f218> in <module>()
      1 fig, ax = plt.subplots()
      2 ax.imshow(mask, cmap='Greys', vmin=np.amin(mask), vmax=np.amax(mask), interpolation='None')
----> 3 ax.set_title(filename)
      4 plt.axis('off')

NameError: name 'filename' is not defined



In [ ]:

    
mask1 = np.zeros(shape=(2400,2400))
mask1[mask==1] = 1



In [ ]:

    
fig, ax = plt.subplots()
ax.imshow(mask1, cmap='Greys', vmin=np.amin(mask1), vmax=np.amax(mask1), interpolation='None')
ax.set_title("mask is 1")
plt.axis('off')



In [ ]:

    
mask2 = np.zeros(shape=(2400,2400))
mask2[mask==2] = 1



In [ ]:

    
fig, ax = plt.subplots()
ax.imshow(mask2, cmap='Greys', vmin=np.amin(mask2), vmax=np.amax(mask2), interpolation='None')
ax.set_title("mask is 2")
plt.axis('off')



In [ ]:

    
dem_filename = '/Users/mj/projects/CHARIS/elevation_data/SRTMGL3_version2_SIN/CHARIS_DEM.v2.0.h25v06.tif'
dem = read_tile(dem_filename)



In [ ]:

    
print(dem.shape)
print(np.amin(dem), np.amax(dem))



In [ ]:

    
fig, ax = plt.subplots()
ax.imshow(dem * mask2, cmap='Greys', vmin=np.amin(dem), vmax=np.amax(dem), interpolation='None')
ax.set_title(dem_filename)
plt.axis('off')



In [ ]: