站点等值填色图

ANDYKYLE

fn='E:/24di.txt'
ncol=numasciicol(fn)#获取文件的列数
nrow=numasciirow(fn)#获取文件的行数
a=asciiread(fn,shape=(nrow,ncol));#读取文件的行、列数
lon=a[:,0]
lat=a[:,1]
diw=a[:,2]

#设置坐标
x=arange(92,109,0.05)
y=arange(32,43,0.05)
gtemp,gx,gy=griddata((lon,lat),diw,xi=(x,y),method='idw', radius=3)

#画图
axesm()
bou1_layer=shaperead('E:/gansu.shp')
mlayer=shaperead('E:/gansu.shp')
geoshow(bou1_layer,edgecolor='gray')
geoshow(mlayer,visible=False)

levs=[0.1,0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9]
cols = [(255,255,255),(0,255,0),(127,255,0),(255,255,0),(255,215,0),(255,128,0),(255,97,0), \
    (255,0,0),(176,23,31),(135,38,87),(255,0,255)]
#layer = contourfm(gx,gy, gtemp,20)
layer=contourfm(x,y,gtemp,levs,colors=cols)
masklayer(mlayer, [layer])
xlim(92,109)
ylim(32,43)
title(u'24小时最低温度预报准确率',fontname=u'黑体',fontsize=20,bold=False,color='blue')
#text(113.3, 23.9, u'2017年5月6日13时', fontname=u'黑体', fontsize=16)
colorbar(layer)
#show()