转换文本文件格点数据为NetCDF格式

MeteoInfo · 发表于 2012-2-1 13:52:30

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

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

x

这个脚本演示了如何将文本文件格点数据转换为NetCDF格式。示例脚本程序中的文本文件是Surfer ASCII 格点数据，是一个全球气压场，是从GrADS的示例文件model.ctl (model.dat)输出的。脚本程序中有注释，这里就不细说了。

文本文件：

pressure_19870102.grd (54.99 KB, 下载次数: 57)

2012-2-1 13:45 上传

点击文件名下载附件

转换后的NetCDF文件：

pressure.nc (27.77 KB, 下载次数: 36)

2012-2-1 13:51 上传

点击文件名下载附件

转换后的NetCDF文件在MeteoInfo中打开：

脚本程序：

#--------------------------------------------------------
# Author: Yaqiang Wang
# Date: 2012-2-1
# Purpose: Convert Surfer ASCII grid data to NetCDF data format
# Note: Sample
#-----------------------------------------------------------
import clr
clr.AddReferenceByPartialName("System.Windows.Forms")
clr.AddReferenceByPartialName("System.Drawing")
from System import *
from System.Windows.Forms import *
from System.Drawing import *
clr.AddReference("MeteoInfoC.dll")
from MeteoInfoC import *
from MeteoInfoC.Data import *
from MeteoInfoC.Data.MeteoData import *
import os.path
from datetime import *
import time
#---- Set data folder
DataDir = "F:\\Temp\"
#---- Open Surfer ASCII Data file
infile = open(DataDir + "pressure_19870102.grd")
i = 0
lonNum = 0
latNum = 0
xMin = 0.0
xMax = 0.0
yMin = 0.0
yMin = 0.0
for aline in infile:
if i == 1:
lonNum = int(aline.split()[0])
latNum = int(aline.split()[1])
elif i == 2:
xMin = float(aline.split()[0])
xMax = float(aline.split()[1])
elif i == 3:
yMin = float(aline.split()[0])
yMax = float(aline.split()[1])
elif i == 4:
break
i = i + 1
lonArray = Array.CreateInstance(Double, lonNum)
latArray = Array.CreateInstance(Double, latNum)
dataArray = Array.CreateInstance(Double, lonNum * latNum)
dx = (xMax - xMin) / (lonNum - 1)
dy = (yMax - yMin) / (latNum - 1)
print xMin, xMax
print yMin, yMax
print dx, dy
for i in range(lonNum):
lonArray[i] = xMin + dx * i
for i in range(latNum):
latArray[i] = yMin + dy * i
i = 0
d = 0
for aline in infile:
if i >= 0:
dataList = aline.split()
# print dataList
for dstr in dataList:
dataArray[d] = float(dstr)
d = d + 1
i = i + 1
infile.close()
unDef = -2.56E33
#---- Create netCDF data info object
outfilepath = DataDir + "pressure.nc"
outDataInfo = NetCDFDataInfo()
outDataInfo.fileName = outfilepath
outDataInfo.UNDEF = unDef
outDataInfo.unlimdimid = 2
#---- Add dimensions: lon, lat, time
lonDim = outDataInfo.AddDimension("lon", lonNum)
latDim = outDataInfo.AddDimension("lat", latNum)
timeDim = outDataInfo.AddDimension("time", -1)
#---- Add variables
outDataInfo.AddVariable("lon", NetCDF4.NcType.NC_DOUBLE, Array[DimStruct]([lonDim]))
outDataInfo.AddVariable("lat", NetCDF4.NcType.NC_DOUBLE, Array[DimStruct]([latDim]))
outDataInfo.AddVariable("time", NetCDF4.NcType.NC_DOUBLE, Array[DimStruct]([timeDim]))
outDataInfo.AddVariable("pressure", NetCDF4.NcType.NC_DOUBLE, Array[DimStruct]([timeDim, latDim, lonDim]))
#---- Add variable attributes
outDataInfo.AddVariableAttribute("lon", "units", "degrees_east")
outDataInfo.AddVariableAttribute("lon", "long_name", "longitude")
outDataInfo.AddVariableAttribute("lat", "units", "degrees_north")
outDataInfo.AddVariableAttribute("lat", "long_name", "latitude")
outDataInfo.AddVariableAttribute("time", "units", "days since 1-1-1 00:00:00")
outDataInfo.AddVariableAttribute("time", "long_name", "time")
outDataInfo.AddVariableAttribute("pressure", "units", "hPa")
outDataInfo.AddVariableAttribute("pressure", "long_name", "Pressure")
outDataInfo.AddVariableAttribute("pressure", "missing_value", unDef)
#---- Add global attributes
outDataInfo.AddGlobalAttribute("title", "Pressure data")
outDataInfo.AddGlobalAttribute("description", "script sample data")
#---- Create netCDF file
outDataInfo.CreateNCFile(outfilepath)
#---- Write lon, lat, time data
outDataInfo.WriteVar("lon", lonArray)
outDataInfo.WriteVar("lat", latArray)
timeArray = Array.CreateInstance(Double, 1)
t = datetime(1987, 1, 2)
st = datetime(1, 1, 1)
days = (t - st).days
timeArray[0] = days
#print timeArray[0]
startArray = Array.CreateInstance(int, 1)
countArray = Array.CreateInstance(int, 1)
startArray[0] = 0
countArray[0] = 1
outDataInfo.WriteVara("time", startArray, countArray, timeArray)
#---- Write data
startArray = Array.CreateInstance(int, 3)
countArray = Array.CreateInstance(int, 3)
startArray[0] = 0
startArray[1] = 0
startArray[2] = 0
countArray[0] = 1
countArray[1] = latDim.dimLen
countArray[2] = lonDim.dimLen
outDataInfo.WriteVara("pressure", startArray, countArray, dataArray)
#---- Close netCDF file
outDataInfo.CloseNCFile()

复制代码

mofangbao · 发表于 2012-2-1 17:05:36

赞一个，不知道楼主有没有把这个功能集成到meteoinfo中的打算呢

MeteoInfo · 发表于 2012-2-1 17:17:52

mofangbao 发表于 2012-2-1 17:05

登录/注册后可看大图

赞一个，不知道楼主有没有把这个功能集成到meteoinfo中的打算呢

用脚本程序是因为其灵活性，稍加修改可以转换别的文本格式数据。如果放在桌面软件中只能是一个比较死板的功能。

mofangbao · 发表于 2012-2-1 17:50:40

MeteoInfo 发表于 2012-2-1 17:17

登录/注册后可看大图

用脚本程序是因为其灵活性，稍加修改可以转换别的文本格式数据。如果放在桌面软件中只能是一个比较死板的 ...

恩说明楼主的dll很强大

孤蓝et · 发表于 2012-2-1 21:55:54

脚本强大啊。。。。。就是不擅长python

孤蓝et · 发表于 2012-2-1 21:56:23

脚本强大啊。。。。。就是不擅长python

♂雨已~○ · 发表于 2012-2-3 23:00:59

楼主，厉害~~~真正学习二次开发中！

浪迹天涯 · 发表于 2012-2-7 16:48:02

普通的ascii数据具体的格式要怎么做？

MeteoInfo · 发表于 2012-2-7 17:33:05

浪迹天涯发表于 2012-2-7 16:48

登录/注册后可看大图

普通的ascii数据具体的格式要怎么做？

随便，自己知道怎么读文件就可以了。

浪迹天涯 · 发表于 2012-2-8 14:04:06

MeteoInfo 发表于 2012-2-7 17:33

登录/注册后可看大图

随便，自己知道怎么读文件就可以了。

嗯，了解~

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