In [4]: a=fromfile("pr.dec",float32).byteswap().reshape(220,250) In [5]: a Out[5]: array([[-999., -999., -999., ..., -999., -999., -999.], [-999., -999., -999., ..., -999., -999., -999.], [-999., -999., -999., ..., -999., -999., -999.], ..., [-999., -999., -999., ..., -999., -999., -999.], [-999., -999., -999., ..., -999., -999., -999.], [-999., -999., -999., ..., -999., -999., -999.]], dtype=float32) In [6]: imshow(a) Out[6]: In [7]: from matplotlib.toolkits.basemap import Basemap In [8]: M=Basemap(urcrnrlat=50,urcrnrlon=150,llcrnrlat=25,llcrnrlon=125, resolution='l') In [9]: M1=Basemap(urcrnrlat=50,urcrnrlon=150,llcrnrlat=25,llcrnrlon=125, resolution='c') In [10]: M.drawcoastlines() In [11]: show() In [12]: M.imshow(a);M.drawco M.drawcoastlines M.drawcountries In [12]: M.imshow(a);M.drawcoastlines() Out[12]: In [13]: M.imshow(flipud(a));M.drawcoastlines() Out[13]: In [14]: M.imshow(flipud(a));M.drawcoastlines();show() Out[14]: In [15]: M1.imshow(flipud(a));M.drawcoastlines();show() \Out[15]: In [16]: M1.imshow(flipud(a));M.drawcoastlines();show() Out[16]: In [17]: M1.imshow(flipud(a));M1.drawcoastlines();show() Out[17]: