In [ ]:

    

import numpy as np


    

In [ ]:

    

import matplotlib as mpl
import matplotlib.pyplot as plt
import matplotlib._cntr as cntr
%matplotlib inline
mpl.rcParams["figure.facecolor"] = "white"
mpl.rcParams["axes.facecolor"] = "white"
mpl.rcParams["savefig.facecolor"] = "white"


    

In [ ]:

    

thefile = '/Users/savoie/tmp/geotiff/nt_201409_f17_v01_n.bins'


    

In [ ]:

    

import rockbag as rb


    

In [ ]:

    

data = rb.suckgrid(thefile, gpd='n3b', header=300)


    

In [ ]:

    

plt.imshow(data)


    

In [ ]:

    

reproj_fn = '/Users/savoie/tmp/geotiff/nt_201409_f17_v01_n.wgs84.intermediate.tif'


    

In [ ]:

    

from PIL import Image
import numpy as np


    

In [ ]:

    

y = Image.open(reproj_fn)


    

In [ ]:

    

img = np.array(y)


    

In [ ]:

    

plt.imshow(img)


    

In [ ]:

    

data


    

In [ ]:

    

img


    

In [ ]:

    

diff = img - data


    

In [ ]:

    

np.any(diff != 0)


    

In [ ]: