- ����
- 22699
- ����
-
- ����
- ����ʱ��
- Сʱ
- ע��ʱ��
- 2011-7-23
- ����¼
- 1970-1-1
|
��¼��鿴���ྫ������~
����Ҫ ��¼ �ſ������ػ�鿴��û���ʺţ�����ע��
x
��������� ƽ������ܲ� �� 2022-11-16 21:24 �༭
Python���ƺ���һ������ô�������⣬��������������������Ѿ��õ��������Ľ��������Basemap��pyshp������ģ�飬ʵ�ֶ�����shp�ļ���������а����ܾ��հ�ݣ�˦��Խ����ͼ����maskout������������Ǻ�shp�߽硣
�����˼·�ǣ���ȡ�������ݣ��ճ���ͼ���õ�matplotlib plot��һ��instance����Ӧδ������ͼ��Ȼ����йؼ���һ������ȡshapefile��shp�ļ�������meteoinfo�ĵ�ͼ���ļ�country1.shp�������������ͼ������ȡshapefile��ĸ���Ȥ���������й������γ�һ��clip�������õ�����ͼ��
Ч��ͼ���£�һ�ź�ƽ����ѹ����ͼ�����й���ĵ���maskout�ˡ�
�Ŵ�֮��������û�о�ݣ������Ǻϡ�
��Ȼ��Ҳ�����ؼ���ͼ���ӵ�����ȥ��ֻҪ����clip����Ҫ����ͬʱ���Ӹ��Ӿ�ϸ���ؼ���ͼ�Ϳ����ˡ�
��һ��ͼ�Ľű����£�����ͼ��ͬС�죩��- #coding=utf-8
- import maskout
- from mpl_toolkits.basemap import Basemap
- import matplotlib.pyplot as plt
- from matplotlib.font_manager import FontProperties
- import numpy
- import netCDF4 as nc
- font = FontProperties(fname=r"C:\\Windows\\Fonts\\impact.ttf")
- font1 = FontProperties(fname=u"C:\\Windows\\Fonts\\simkai.ttf")
- fig = plt.figure(figsize=(16,9.6))
- ax = fig.add_subplot(111)
- # ncdata=nc.Dataset(r'2004_500_grtuvw.nc')
- ncdata=nc.Dataset('pres.mon.ltm.nc')
- # data=ncdata.variables['t'][0,:,:]
- data=ncdata.variables['pres'][0,:,:]
- # exit()
- # lat=ncdata.variables['latitude'][:]
- # lon=ncdata.variables['longitude'][:]
- lat=ncdata.variables['lat'][:]
- lon=ncdata.variables['lon'][:]
- lon1,lon2=lon[0],lon[-1]
- lat1,lat2=lat[-1],lat[0]
- nx=data.shape[1];ny=data.shape[0]
- m = Basemap(llcrnrlon=lon1,llcrnrlat=lat1,urcrnrlon=lon2,urcrnrlat=lat2,projection = 'cyl')
- xx,yy=m.makegrid(nx,ny)
- yy=yy[::-1]
- m.readshapefile('country1','whatevername',color='gray')
- minval,maxval=int(numpy.amin(data)),int(numpy.amax(data))+1
- cs = m.contourf(xx,yy,data,range(minval,maxval),cmap= plt.cm.get_cmap('jet') )
- bar=m.colorbar(cs)
- bar.set_ticks(range(minval-1,maxval,40))
- bar.set_ticklabels(range(minval-1,maxval,40))
- clip=maskout.shp2clip(cs,ax,'country1','China')
- plt.title(u'Python Super Mask',fontproperties=font,fontsize=40)
- lon1,lon2=ax.set_xlim(70,140);lat1,lat2=ax.set_ylim(15,55)
- #signature=u"by ƽ������ܲ�"
- #plt.text(lon2-(lon2-lon1)*3.0/10,lat2+(lat2-lat1)*0.1/10,signature,fontproperties=font1,fontsize=25)
- # plt.show()
- plt.savefig('111.png')
���ƴ��� ���������clip=maskout.shp2clip(cs,ax,'country1','China') ������maskout���Լ������һ��ģ�飬ʵ�ֵĹ��ܾ�����ؼ��Ķ���ͼ���м��У���Ϊģ�������ֲ�����������У������Ժ�ʹ�á�
�ѵ�һ��ͼ�IJ������ݡ���ͼ�ͽű�һ�����ˣ������Զ���maskoutģ�������ӭ����Ȥ���������أ�
Python������_�������ݼ��ű�.rar
(1.41 MB, ���ش���: 4638)
|
����
-
�鿴ȫ������
|