netcdf - How to plot smooth line in python? -


ncep data website. want plot picture this: wind field , heights or one(this 1 add trough lines): wind field , heights

my data different them, content different. but, method should same.

i dont't know how smooth line. result: my_picture

here code:

import numpy np import scipy import matplotlib.pyplot plt netcdf4 import dataset scipy.interpolate import rbf scipy.ndimage import zoom  mpl_toolkits.basemap import basemap  m=basemap(projection='cyl',llcrnrlat=20,urcrnrlat=50,llcrnrlon=90,urcrnrlon=130)  chnshp = 'd:\python\shapefile\data\chn_adm_shp\chn_adm1' m.readshapefile(chnshp,'chn',drawbounds = false) twnshp = 'd:\python\shapefile\data\twn_adm_shp\twn_adm0' m.readshapefile(twnshp,'twn',drawbounds = false)  info, shape in zip(m.chn_info, m.chn):         x, y = zip(*shape)          m.plot(x, y, marker=none,color='k',linewidth = 0.5) info, shape in zip(m.twn_info, m.twn):         x, y = zip(*shape)          m.plot(x, y, marker=none,color='k',linewidth = 0.5)  parallels = np.arange(-90.,91.,10.) m.drawparallels(parallels,labels=[1,0,0,1],linewidth=0.5,xoffset=1.2) meridians = np.arange(-180.,181.,10.) m.drawmeridians(meridians,labels=[1,0,0,1],linewidth=0.5,yoffset=1.2)  u=dataset(r'd:\python\try\ncep\uwnd.2016.nc','r') v=dataset(r'd:\python\try\ncep\vwnd.2016.nc','r') hgt_data=dataset(r'd:\python\try\ncep\hgt.2016.nc','r')  uwnd=u.variables['uwnd'][728][2][:] vwnd=v.variables['vwnd'][728][2][:] hgt=hgt_data.variables['hgt'][728][2][:] lat=u.variables['lat'][:] lon=u.variables['lon'][:]  index1=np.logical_and(lon>=90,lon<=130);index2=np.logical_and(lat>=20,lat<=50) lons=lon[index1];lats=lat[index2] u1=uwnd[index2,:];u2=u1[:,index1] v1=vwnd[index2,:];v2=v1[:,index1] hgt1=hgt[index2,:];hgt2=hgt1[:,index1]  nx,ny=np.meshgrid(lons,lats) x,y=m(nx,ny) q = m.quiver(x,y,u2,v2,scale=250,width=0.003) qk = plt.quiverkey(q, 0.85, -0.12, 20, '20 m/s', labelpos='n')  rbf = scipy.interpolate.rbf(x, y, hgt2) zi = rbf(x, y) plt.contour(x,y,zi,color='k')  plt.show()

update:

lons = zoom(lons,3,order=3) lats = zoom(lats,3,order=3) x,y  = np.meshgrid(lons,lats,copy=false) hgt2 = zoom(hgt2,3,order=3) cs = plt.contour(x,y,hgt2,levels=levels,colors='k',linewidths=0.7)

have @ examples contour() function on matplotlib site https://matplotlib.org/examples/pylab_examples/contour_demo.html

heres how generate x , y coordinates plot:

delta = 0.025 x = np.arange(-3.0, 3.0, delta) y = np.arange(-2.0, 2.0, delta) x, y = np.meshgrid(x, y)

what need increase resolution lons , lats fields use in meshgrid() function in own program.

higher resolution -> smoother lines


Comments

Post a Comment

Popular posts from this blog

inversion of control - Autofac named registration constructor injection -

does autofac support specifying registration name in components' constructors? example: ninject's namedattribute . you need use autofac.extras.attributed package on top achieve this so let's have 1 interface , 2 classes: public interface ihello { string sayhello(); } public class englishhello : ihello { public string sayhello() { return "hello"; } } public class frenchhello : ihello { public string sayhello() { return "bonjour"; } } then have consumer class, in want select instance injected: public class helloconsumer { private readonly ihello helloservice; public helloconsumer([withkey("en")] ihello helloservice) { if (helloservice == null) { throw new argumentnullexception("helloservice"); } this.helloservice = helloservice; } public string sayhello() { return this.helloservice.sayhello()...
Read more

verilog - Systemverilog dynamic casting issues -

i've code snippet following in testbench function void write_to_port( my_data_type_base data ); my_data_type_extended data_ext; if(!$cast(data_ext, data)); `uvm_error(get_type_name(), "failed cast"); `uvm_info(get_name(), $psprintf("data_ext :\n%s", data_ext.sprint()), uvm_medium) // write data_ext out port.... endfunction when run it, i'm getting uvm_error "failed cast." i'm not quire sure why $cast not returning 1. can see, i'm printing out extended class data item after casting uvm_info. can see it's being cast properly. if don't use $cast if condition, don't runtime error. isn't coding practice use if dynamic cast check if $cast returning 1 ? what might reason behind cast not returning 1 in above case? i think semicolon on line 'if' not belong? i think consumes if statement, , uvm_error executes regardless of how if evaluates: if(!$cast(data_ext, data)); <- no semico...
Read more

ios - Change Storyboard View using Seague -

i have following code checks connection server, , depending on want change views if no connection found. nsstring *connect = [nsstring stringwithcontentsofurl:[nsurl urlwithstring:@"http://myserver.com"] encoding:nsutf8stringencoding error:nil]; if (connect == null) { [self performseguewithidentifier:@"network_failure" sender:self]; } else { // other code here } my problem not change views. have added storyboard segue between 2 views identifier of network_failure . not have navigation controller or that, problem? thanks help! no not problem. first of all, in objectivec exists nil object oriented null . must use that. so: if(connect == nil) or better: if(!connect) so due fact nil not same of null , probably code does't enter in if scope . you have been able discover breakpoint. anyway: be sure set trigger segue (line) 2 view controller , , name network_failure (that not best name. should instead view1toview2 ) as, m...
Read more