def test_pcolormesh(self):
''' Should use Nansatmap.pcolormesh '''
n = Nansat(self.test_file_stere, logLevel=40)
b1 = n[1]
nmap = Nansatmap(n)
nmap.pcolormesh(b1)
tmpfilename = os.path.join(ntd.tmp_data_path, 'nansatmap_pcolormesh.png')
nmap.save(tmpfilename)
self.assertTrue(os.path.exists(tmpfilename))
def test_pcolormesh(self):
''' Should use Nansatmap.pcolormesh '''
n = Nansat(self.test_file_stere, logLevel=40)
b1 = n[1]
nmap = Nansatmap(n)
nmap.pcolormesh(b1)
tmpfilename = os.path.join(ntd.tmp_data_path,
'nansatmap_pcolormesh.png')
nmap.save(tmpfilename)
self.assertTrue(os.path.exists(tmpfilename))
nMap = Nansatmap(n) # draw filled contour plot nMap.contourf(w, v=range(4, 22, 2)) # draw black smooth contour plot with labels nMap.contour(w, smooth=True, fontsize=8, colors='k') # add colorbar nMap.add_colorbar(fontsize=10) # add geocoordinates nMap.drawgrid() # save to a file nMap.save(oFileName + '01_contourf_contour.png', landmask=False) # Create Nansatmap object from nansat (domain) object nMap = Nansatmap(n, resolution='l') # pseudo-color plot over the map nMap.pcolormesh(w) # quiver plot nMap.quiver(v, u, quivectors=20) # save to a file nMap.save(oFileName + '02_pcolormesh_quiver.png') # use Nansatmap for converting lon/lat into x/y # 1. Create domain over the area of interest in stereographic projection extentString = '-lle -10 50 20 70 -tr 1000 1000' srsString = '+proj=stere +lon_0=10 +lat_0=60 +k=1 +ellps=WGS84 +datum=WGS84 +no_defs' d = Domain(srsString, extentString) # 2. Create nansatmap object from tha domain nmap = Nansatmap(d) # 3. Use the created object to convert from lon/lat into x/y: x, y = nmap([0, 2, 4], [63, 64, 65]) # 4 or from x/y into lon/lat
# Create a Nansat object (n)
n = Nansat(iFileName)
# Get data from 1st, 2nd and 3rd bands as numpy array (u,v and w)
u = n[1]; v = n[2]; w = n[3]
nMap = Nansatmap(n)
# draw filled contour plot
nMap.contourf(w)
# draw black smooth contour plot with labels
nMap.contour(w, smooth=True, fontsize=8, colors='k')
# add colorbar
nMap.add_colorbar(fontsize=10)
# add geocoordinates
nMap.drawgrid()
# save to a file
nMap.save(oFileName+'contourf_contour.png', landmask=False)
nMap = Nansatmap(n, resolution='h')
# pseudo-color plot over the map
nMap.pcolormesh(w)
# quiver plot
nMap.quiver(u, v, quivectors=20)
nMap.save(oFileName+'pcolormesh_quiver.png')
#n.write_map(grid='w', contour='w', smooth=true, quiver=('u','v'), step=5)
