MeteoInfoLab脚本示例：160站数据绘图

MeteoInfo · 发表于 2016-5-10 11:01:31

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 MeteoInfo 于 2019-7-31 14:14 编辑

中国气候160站数据下载地址：http://ncc.cma.gov.cn/Website/index.php?ChannelID=43&WCHID=5，下面用到的示例数据：

t1601.txt (52.59 KB, 下载次数: 214)

2018-1-24 11:56 上传

点击文件名下载附件。站点数据按站、年的顺序排列，从1951年开始。

数据说明里有站点名称和经纬度信息，经过简单的编辑成此文件：

st_160.csv (4.75 KB, 下载次数: 413)

2016-5-10 10:55 上传

点击文件名下载附件。

读取2015年1月的温度数据（0.1度为单位）并绘制站点图，标注站点名：

#Get file names
datafn = 'D:/Temp/ascii/t1601.txt'
stfn = 'D:/Temp/ascii/st_160.csv'
#Get years
years = []
for i in range(1951, 2017):
years.append(i)
yearn = len(years)
stn = 160
#Read temperature data of 160 stations from data file
temp = asciiread(datafn, shape=(yearn,stn))
temp = temp * 0.1
#Get temperature data of 2015
year = 2015
yearidx = years.index(year)
temp = temp[yearidx,:]
#Read station name and lon/lat
table = readtable(stfn, delimiter=',', format='%i%2s%2f')
stnames = table['Name']
lat = table['LAT']
lon = table['LON']
#Plot
axesm()
geoshow('cn_province', edgecolor='lightgray')
bou1_layer = geoshow('cn_border', facecolor=(0,0,255))
layer = scatterm(lon, lat, temp, 20)
layer.addfield('Name', 'string', stnames)
layer.addlabels('Name', fontname=u'楷体', fontsize=14, yoffset=15)
title('Temperature (' + str(year) + '-01)')
colorbar(layer)
xlim(72, 136)
ylim(16, 55)
#Add south China Sea
sc_layer = bou1_layer.clone()
axesm(position=[0.14,0.18,0.15,0.2], axison=False, frameon=True)
geoshow(sc_layer, facecolor=(0,0,255))
xlim(106, 123)
ylim(2, 23)

#Get file names 
datafn = 'D:/Temp/ascii/t1601.txt' 
stfn = 'D:/Temp/ascii/st_160.csv' 
 
#Get years 
years = [] 
for i in range(1951, 2017): 
    years.append(i) 
yearn = len(years) 
stn = 160 
 
#Read temperature data of 160 stations from data file 
temp = asciiread(datafn, shape=(yearn,stn)) 
temp = temp * 0.1 
 
#Get temperature data of 2015 
year = 2015 
yearidx = years.index(year) 
temp = temp[yearidx,:] 
#Read station name and lon/lat 
table = readtable(stfn, delimiter=',', format='%i%2s%2f') 
stnames = table['Name'] 
lat = table['LAT'] 
lon = table['LON'] 
 
#Plot 
axesm() 
geoshow('cn_province', edgecolor='lightgray') 
bou1_layer = geoshow('cn_border', facecolor=(0,0,255)) 
layer = scatterm(lon, lat, temp, 20) 
layer.addfield('Name', 'string', stnames) 
layer.addlabels('Name', fontname=u'楷体', fontsize=14, yoffset=15) 
title('Temperature (' + str(year) + '-01)') 
colorbar(layer) 
xlim(72, 136) 
ylim(16, 55) 
 
#Add south China Sea 
sc_layer = bou1_layer.clone() 
axesm(position=[0.14,0.18,0.15,0.2], axison=False, frameon=True) 
geoshow(sc_layer, facecolor=(0,0,255)) 
xlim(106, 123) 
ylim(2, 23)

读取2015年1月的温度数据（0.1度为单位），插值为格点数据并绘制等值线填色图：

#Get file names
datafn = 'D:/Temp/ascii/t1601.txt'
stfn = 'D:/Temp/ascii/st_160.csv'
#Get years
years = []
for i in range(1951, 2017):
years.append(i)
yearn = len(years)
stn = 160
#Read temperature data of 160 stations from data file
temp = asciiread(datafn, shape=(yearn,stn))
temp = temp * 0.1
#Get temperature data of 2015
year = 2015
yearidx = years.index(year)
temp = temp[yearidx,:]
#Read station name and lon/lat
table = readtable(stfn, delimiter=',', format='%i%2s%2f')
stnames = table['Name']
lat = table['LAT']
lon = table['LON']
#To grid data
x = arange(70, 140, 0.5)
y = arange(15, 58, 0.5)
gtemp,gx,gy = griddata((lon, lat), temp, xi=(x, y), method='idw', pointnum=4)
#Plot
axesm()
geoshow('cn_province', edgecolor='lightgray')
geoshow('cn_border', facecolor=(0,0,255))
china_layer = geoshow('china', visible=False)
layer = contourfm(x, y, gtemp, 20)
slayer = scatterm(lon, lat, temp, colors='k', size=4)
masklayer(china_layer, [layer])
title('Temperature (' + str(year) + '-01)')
colorbar(layer)
xlim(72, 136)
ylim(16, 55)
#Add south China Sea
axesm(position=[0.14,0.18,0.15,0.2], axison=False, frameon=True)
geoshow('cn_border', facecolor=(0,0,255))
xlim(106, 123)
ylim(2, 23)

#Get file names 
datafn = 'D:/Temp/ascii/t1601.txt' 
stfn = 'D:/Temp/ascii/st_160.csv' 
 
#Get years 
years = [] 
for i in range(1951, 2017): 
    years.append(i) 
yearn = len(years) 
 
stn = 160 
#Read temperature data of 160 stations from data file 
temp = asciiread(datafn, shape=(yearn,stn)) 
temp = temp * 0.1 
 
#Get temperature data of 2015 
year = 2015 
yearidx = years.index(year) 
temp = temp[yearidx,:] 
 
#Read station name and lon/lat 
table = readtable(stfn, delimiter=',', format='%i%2s%2f') 
stnames = table['Name'] 
lat = table['LAT'] 
lon = table['LON'] 
 
#To grid data 
x = arange(70, 140, 0.5) 
y = arange(15, 58, 0.5) 
gtemp,gx,gy = griddata((lon, lat), temp, xi=(x, y), method='idw', pointnum=4) 
 
#Plot 
axesm() 
geoshow('cn_province', edgecolor='lightgray') 
geoshow('cn_border', facecolor=(0,0,255)) 
china_layer = geoshow('china', visible=False) 
layer = contourfm(x, y, gtemp, 20) 
slayer = scatterm(lon, lat, temp, colors='k', size=4) 
masklayer(china_layer, [layer]) 
title('Temperature (' + str(year) + '-01)') 
colorbar(layer) 
xlim(72, 136) 
ylim(16, 55) 
 
#Add south China Sea 
axesm(position=[0.14,0.18,0.15,0.2], axison=False, frameon=True) 
geoshow('cn_border', facecolor=(0,0,255)) 
xlim(106, 123) 
ylim(2, 23)

读取北京站点1951-2016年一月份温度变化数据并绘图：

#Get file names
datafn = 'D:/Temp/ascii/t1601.txt'
stfn = 'D:/Temp/ascii/st_160.csv'
#Get years
years = []
for y in range(1951, 2017):
years.append(y)
yearn = len(years)
stn = 160
#Read temperature data of 160 stations from data file
temp = asciiread(datafn, shape=(yearn,stn))
temp = temp * 0.1
#Read station name and lon/lat
table = readtable(stfn, delimiter=',', format='%i%2s%2f')
stnames = table['Name']
#Get 1951-2016 Jan. temperature data of 北京
stname = u'北京'
stidx = stnames.index(stname)
temp = temp[:,stidx]
#Plot
plot(years, temp, color='r', linewidth=2)
xlabel(u'年份', fontname=u'宋体')
ylabel(u'温度', fontname=u'宋体')
title(u'北京1951-2016年一月份温度年际变化', fontname=u'黑体')

river · 发表于 2016-5-10 11:56:50

支持王老师！棒棒哒！

Ashe · 发表于 2016-5-10 13:51:57

必须赞，王老师正是无私啊，谢谢

aas314 · 发表于 2016-5-10 22:32:05

必须学习大赞

wsjqxjy · 发表于 2016-5-11 08:53:38

好好学习，太给力了，谢谢分享

勿忘草 · 发表于 2016-5-11 08:56:10

太给力了，谢谢分享

我是彤-Tina · 发表于 2016-5-12 09:25:17

看起来就很棒的样子~

湖图淳淙淙 · 发表于 2016-5-12 16:13:18

谢谢，刚一找就找到了。。。。。

FRUITKID · 发表于 2016-5-13 10:59:42

前排撒花！

drkuang · 发表于 2016-5-17 09:46:02

请问有没有中文版的软件使用教程，入门使用和你发表的示例都列出来。我现安了软件，但如何用，还摸不着门。

		自动登录	找回密码
密码			立即注册