def test_imshow(self):
''' Should use Nansatmap.imshow '''
n = Nansat(self.test_file_stere, logLevel=40)
b1 = n[1]
nmap = Nansatmap(n)
nmap.imshow(b1, cmap='ak01')
tmpfilename = os.path.join(ntd.tmp_data_path, 'nansatmap_imshow.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))
def test_imshow(self):
''' Should use Nansatmap.imshow '''
n = Nansat(self.test_file_stere, logLevel=40)
b1 = n[1]
nmap = Nansatmap(n)
nmap.imshow(b1, cmap='ak01')
tmpfilename = os.path.join(ntd.tmp_data_path, 'nansatmap_imshow.png')
nmap.save(tmpfilename)
self.assertTrue(os.path.exists(tmpfilename))
def test_imshow_random(self):
''' Should use Nansatmap.imshow '''
n = Nansat(self.test_file_stere, logLevel=40)
b1 = n[1]
nmap = Nansatmap(n)
nmap.imshow(b1/5, cmap='random')
nmap.add_colorbar()
tmpfilename = os.path.join(ntd.tmp_data_path, 'nansatmap_imshow_random.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))
def test_imshow_random(self):
''' Should use Nansatmap.imshow '''
n = Nansat(self.test_file_stere, logLevel=40)
b1 = n[1]
nmap = Nansatmap(n)
nmap.imshow(b1 / 5, cmap='random')
nmap.add_colorbar()
tmpfilename = os.path.join(ntd.tmp_data_path,
'nansatmap_imshow_random.png')
nmap.save(tmpfilename)
self.assertTrue(os.path.exists(tmpfilename))
def test_add_labels(self):
size, npo = 100, 10
xy = np.random.randint(0, size, npo*2).reshape(npo, 2)
z = np.random.randint(0, size, npo)
xg, yg = np.meshgrid(range(size), range(size))
zg = griddata(xy, z, np.dstack([xg, yg]), method='nearest')
dstDomain = Domain(NSR().wkt, '-te -10 -10 10 10 -ts 100 100')
nmap = Nansatmap(dstDomain)
nmap.imshow(zg, cmap='random')
nmap.add_zone_labels(zg, fontsize=10)
tmpfilename = os.path.join(ntd.tmp_data_path, 'nansatmap_zonelables.png')
nmap.save(tmpfilename)
self.assertTrue(os.path.exists(tmpfilename))
def test_add_labels(self):
size, npo = 100, 10
xy = np.random.randint(0, size, npo * 2).reshape(npo, 2)
z = np.random.randint(0, size, npo)
xg, yg = np.meshgrid(range(size), range(size))
zg = griddata(xy, z, np.dstack([xg, yg]), method='nearest')
dstDomain = Domain(NSR().wkt, '-te -10 -10 10 10 -ts 100 100')
nmap = Nansatmap(dstDomain)
nmap.imshow(zg, cmap='random')
nmap.add_zone_labels(zg, fontsize=10)
tmpfilename = os.path.join(ntd.tmp_data_path,
'nansatmap_zonelables.png')
nmap.save(tmpfilename)
self.assertTrue(os.path.exists(tmpfilename))
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] # Create Nansatmap object from nansat (domain) object 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)
# 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)
