Das deutsche Python-Forum

Code: Alles auswählen

1 cities = ['London ', 'New York ', 'Madrid ', 'Cairo ', 'Moscow ', 'Delhi ', 'Dakar ']
lat = [51.507778 , 40.716667 , 40.4 , 30.058 , 55.751667 , 28.61 , 14.692778]
 lon = [ -0.128056 , -74, -3.683333 , 31.229 , 37.617778 , 77.23 , -17.446667]

 m = Basemap ( projection ='ortho ', lat_0 =45 , lon_0 =10)
 m. drawmapboundary ()
 m. drawcoastlines ()
 m. fillcontinents ()
 x, y = m(lon , lat )
 plt . plot (x, y, 'ro ')
 for city , xc , yc in zip(cities , x, y):
 plt. text (xc +250000 , yc -150000 , city , bbox = dict ( facecolor ='yellow ', alpha =0.5))

Code: Alles auswählen

#Map create
# projection, lat/lon extents and resolution of polygons to draw
# resolutions: c - crude, l - low, i - intermediate, h - high, f - full
print('Creat a Basemap instance')

llcrnrlon=-29.	#Lon. definition left border
llcrnrlat=1.	#Lat. bottom border
urcrnrlon=29.	#lon. right
urcrnrlat=39.	#Lat. top
map = Basemap(projection='merc',llcrnrlon=llcrnrlon,llcrnrlat=llcrnrlat,urcrnrlon=urcrnrlon,urcrnrlat=urcrnrlat,resolution='i')

print('Location: \nEast:',llcrnrlon,'\nWest:',urcrnrlon,'\nNorth: ',urcrnrlat,'\nSouth: ',llcrnrlat)

#add the states, counties, coasts, national borders and a land-sea colored mask
map.drawcoastlines()
map.drawstates()
map.drawcountries()
map.drawlsmask(land_color='Linen', ocean_color='#CCFFFF') # can use HTML names or codes for colors
map.drawcounties() # you can even add counties (and other shapefiles!)

#plt.show()

#add Lat/Lon
parallels = np.arange(0,40,5.) # make latitude lines ever 5 degrees from 30N-50N
meridians = np.arange(-30,30,5.) # make longitude lines every 5 degrees from 95W to 70W
map.drawparallels(parallels,labels=[1,0,0,0],fontsize=10)
map.drawmeridians(meridians,labels=[0,0,0,1],fontsize=10)

#Cities definition for Visualisation
cities = ['Niamey', 'Bamako', 'Dakar', 'Bobo-Dioulasso','Marrakesch', 'Tunis', 'Abuja']
lat = [13.511596 ,  12.617098 ,  14.716677, 11.164922 ,  31.629472 , 36.806495,  9.076479]
lon = [ 2.125385 ,  -7.981084,  -17.467686 ,  -4.305154 ,  -7.981084 , 10.181532,  7.398574]

#Plot the data on top of the map
lon,lat = np.meshgrid(lons,lats)
x,y = map(lon,lat)

# display blue marble image ( from NASA )
# as map background
map.bluemarble()
for city , xc , yc in zip(cities , x, y):
	plt.text(xc, yc , city , bbox=dict(
	facecolor='yellow', alpha =0.5))

Code: Alles auswählen

In [27]: for city , xc , yc in zip(cities , x, y):
    ...: ^Iplt.text(xc, yc , city , bbox=dict(
    ...: ^Ifacecolor='yellow', alpha =0.5))
    ...:  

In [28]: plt.show()
---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/backends/backend_qt5agg.py in __draw_idle_agg(self, *args)
    174     def __draw_idle_agg(self, *args):
    175         try:
--> 176             FigureCanvasAgg.draw(self)
    177             self.update()
    178         finally:

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/backends/backend_agg.py in draw(self)
    472 
    473         try:
--> 474             self.figure.draw(self.renderer)
    475         finally:
    476             RendererAgg.lock.release()

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/artist.py in draw_wrapper(artist, renderer, *args, **kwargs)
     59     def draw_wrapper(artist, renderer, *args, **kwargs):
     60         before(artist, renderer)
---> 61         draw(artist, renderer, *args, **kwargs)
     62         after(artist, renderer)
     63 

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/figure.py in draw(self, renderer)
   1157         dsu.sort(key=itemgetter(0))
   1158         for zorder, a, func, args in dsu:
-> 1159             func(*args)
   1160 
   1161         renderer.close_group('figure')

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/artist.py in draw_wrapper(artist, renderer, *args, **kwargs)
     59     def draw_wrapper(artist, renderer, *args, **kwargs):
     60         before(artist, renderer)
---> 61         draw(artist, renderer, *args, **kwargs)
     62         after(artist, renderer)
     63 

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/axes/_base.py in draw(self, renderer, inframe)
   2322 
   2323         for zorder, a in dsu:
-> 2324             a.draw(renderer)
   2325 
   2326         renderer.close_group('axes')

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/artist.py in draw_wrapper(artist, renderer, *args, **kwargs)
     59     def draw_wrapper(artist, renderer, *args, **kwargs):
     60         before(artist, renderer)
---> 61         draw(artist, renderer, *args, **kwargs)
     62         after(artist, renderer)
     63 

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/text.py in draw(self, renderer)
    747 
    748         with _wrap_text(self) as textobj:
--> 749             bbox, info, descent = textobj._get_layout(renderer)
    750             trans = textobj.get_transform()
    751 

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/text.py in _get_layout(self, renderer)
    332         of a rotated text when necessary.
    333         """
--> 334         key = self.get_prop_tup()
    335         if key in self._cached:
    336             return self._cached[key]

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/text.py in get_prop_tup(self)
    901         need to know if the text has changed.
    902         """
--> 903         x, y = self.get_unitless_position()
    904         return (x, y, self.get_text(), self._color,
    905                 self._verticalalignment, self._horizontalalignment,

/home/samothkociok/.conda/envs/thomypy/lib/python3.4/site-packages/matplotlib/text.py in get_unitless_position(self)
    883         # This will get the position with all unit information stripped away.
    884         # This is here for convienience since it is done in several locations.
--> 885         x = float(self.convert_xunits(self._x))
    886         y = float(self.convert_yunits(self._y))
    887         return x, y

TypeError: only length-1 arrays can be converted to Python scalars

Code: Alles auswählen

#Load in the netCDF file
nc = NetCDFFile(ROOT+'LE_v2-1-000_conc-sfc_20170501.nc') 

# Extract data from NetCDF file (Global Parameter)
lats = nc.variables['latitude'][:]		#lat; Unit degrees_east
lons = nc.variables['longitude'][:]		#lon; Unit degrees_north
lev = nc.variables['level'][:]			#Layer; Hight
time = nc.variables['time'][:]			#Time; Sevonds since 2012-01-01
	
#Concentrations
accum_dust = nc.variables['tdust'][1,:,:,:]			#total dust Aerosol
accum_dust_units = nc.variables['tdust'].units		#kg m-3

#Map create
# projection, lat/lon extents and resolution of polygons to draw
# resolutions: c - crude, l - low, i - intermediate, h - high, f - full
print('Creat a Basemap instance')

llcrnrlon=-29.	#Lon. definition left border
llcrnrlat=1.	#Lat. bottom border
urcrnrlon=29.	#lon. right
urcrnrlat=39.	#Lat. top
map = Basemap(projection='merc',llcrnrlon=llcrnrlon,llcrnrlat=llcrnrlat,urcrnrlon=urcrnrlon,urcrnrlat=urcrnrlat,resolution='i')

print('Location: \nEast:',llcrnrlon,'\nWest:',urcrnrlon,'\nNorth: ',urcrnrlat,'\nSouth: ',llcrnrlat)

#add the states, counties, coasts, national borders and a land-sea colored mask
map.drawcoastlines()
map.drawstates()
map.drawcountries()
map.drawlsmask(land_color='Linen', ocean_color='#CCFFFF') # can use HTML names or codes for colors
map.drawcounties() # you can even add counties (and other shapefiles!)

#plt.show()

#add Lat/Lon
parallels = np.arange(0,40,5.) # make latitude lines ever 5 degrees from 30N-50N
meridians = np.arange(-30,30,5.) # make longitude lines every 5 degrees from 95W to 70W
map.drawparallels(parallels,labels=[1,0,0,0],fontsize=10)
map.drawmeridians(meridians,labels=[0,0,0,1],fontsize=10)

#Cities definition for Visualisation
cities = ['Niamey', 'Bamako', 'Dakar', 'Bobo-Dioulasso','Marrakesch', 'Tunis', 'Abuja']
lati = [13.511596 ,  12.617098 ,  14.716677, 11.164922 ,  31.629472 , 36.806495,  9.076479]
loni = [ 2.125385 ,  -7.981084,  -17.467686 ,  -4.305154 ,  -7.981084 , 10.181532,  7.398574]

#Plot the data on top of the map
lon,lat = np.meshgrid(lons,lats)
x,y = map(lon,lat)
xc,yc=map(lati,loni)

# display blue marble image ( from NASA )
# as map background
#map.bluemarble()
for city , xc , yc in zip(cities , x, y):
	plt.text(xc, yc , city , bbox=dict(facecolor='yellow', alpha =0.5))

accum_dust_corr= abs(accum_dust)

#Plot Temperature cmap=plt.cm.RdBu_r für Temp
tdustplot= map.pcolor(x,y,np.squeeze(accum_dust_corr), 
norm=colors.LogNorm(0,10),cmap=plt.cm.RdBu_r)

plt.title(nc.variables['tdust'].long_name +' for 01.05.2017 0100 UTC' )
cbar = map.colorbar(tdustplot, location='right', pad="10%")
cbar.set_label(accum_dust_units)

plt.show()