NCL批量读取台站数据中某一值

蔷薇花烬 · 发表于 2017-4-20 15:46:39

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

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

x

本帖最后由蔷薇花烬于 2017-4-20 16:01 编辑

ncl新手小白求各位大神帮助！！！

{:5_275:}
最近在用ncl做毕业论文，需要处理台站数据资料，求得6、7、8月的降水量的平均值画图，台站数据是.txt格式的，网上查找了相关的读取处理方法，但还是不知道怎么样提取出具体月份降水量这一数据，跪求各位大大能不能提供类似的例子参考一下，数据的格式，前两行都是表头，然后两行是一组。月份排列有的是乱的，没按照顺序来，真心不知道怎么办了。。。。

蔷薇花烬 · 发表于 2017-4-20 15:48:46

重新上传一次

yesyesnotion · 发表于 2017-4-20 16:20:18

打开数据的这个软件换一个试试呢，排列的错位了都

蔷薇花烬 · 发表于 2017-4-20 16:26:02

我是直接用linux系统下自带的文本编辑器打开的，我还以为这数据就长这样呢，也是一脸的懵逼

，那兄台有没有相关的程序可以分享一下呢？真的是搞的晕头转向的。。。。

梁苏洁 · 发表于 2017-4-28 14:43:56

我只会这样：把表头去掉，asc_file = "/cygdrive/*/50136.txt"
data = asciiread(asc_file,-1,"string")，这样把所有数据读入

梁苏洁 · 发表于 2017-4-28 14:45:00

ncols = numAsciiCol(asc_file)
nrows = numAsciiRow(asc_file) ; Don't include first row which doesn't have numerical data
设置行列数，你这个应该是19列吧，行数不知道，没下载

梁苏洁 · 发表于 2017-4-28 14:45:29

values = asciiread(asc_file,(/nrows,ncols/),"float")
这样不知道可不可以啊，试试？

梁苏洁 · 发表于 2017-4-28 14:47:25

data = onedtond(asc_file,(/nrows,ncols/)) ;
忘记，onedtond 命令啥意思了

梁苏洁 · 发表于 2017-4-28 14:48:56

iyear = data(:,0)
imon = data(:,1)

梁苏洁 · 发表于 2017-4-28 14:50:23

  fname  = "/cygdrive/*/50136.txt"
; Read all data into a one-dimensional variable.
  dummy  = asciiread(fname,-1,"float")
; Count the number of fields, just to show it can be done.
  ncol = 19                               ; # of columns
  npts = dimsizes(dummy)/ncol          ; # of points
data01 = onedtond(dummy,(/npts,ncol/)) ; npts x ncol
iyear    = data(:,0)
...

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

NCL批量读取台站数据中某一值

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

回帖奖励 +1 金钱

浏览过的版块