ncl站点数据绘图的问题

zly4814624

本人是新手，现在需要讲一个站点数据进行绘图，数据的格式是TXT，第一行是站点序号，第二行是经度，第三是维度，第四行是运用到的温度数据，请问如何将我希望运用到的温度数据变成带有经纬度信息的二维数组呢。望各路大神帮帮忙解答一下。以下是我的原始数据：

xuebiz

http://www.ncl.ucar.edu/Applications/read_ascii.shtml
官网例子~~

可能这个例子比较像你的数据，供参考
stn_latlon.dat - a file with 980 rows and 10 columns of floating point data. (example)
http://www.ncl.ucar.edu/Applications/read_ascii.shtml#stn_latlon

stn_latlon.dat - a file with 980 rows and 10 columns of floating point data.

The first two methods show how to read this file if you know the exact number of rows and columns, and the third method shows how to read this file if you don't.

Method 1

; Read data into a 980 x 10 float array.
  nrows = 980
  ncols = 10
  data  = asciiread("stn_latlon.dat",(/nrows,ncols/),"float")
  printVarSummary(data)           ; Print information about file only.

; Two ways to print the data.
  print(data)                     ; Print data, one value per line
  write_matrix(data,ncols + "f7.2",0)   ; Formatted output

Method 2

This file is actually a file of latitude and longitude values, each dimensioned 70 x 70. The latitude values are written first on the file, followed by the longitude values. Given this information, here's another way to read in this file:

  nlat     = 70
  nlon     = 70
  latlon2d = asciiread("stn_latlon.dat",(/2,nlat,nlon/),"float")     ; 2 x 70 x 70
  lat2d    = latlon2d(0,:,:)    ; 70 x 70
  lon2d    = latlon2d(1,:,:)    ; 70 x 70

Method 3

Use the special contributed functions numAsciiCol and readAsciiTable function to first calculate the number of columns, and then to read the data into an array dimensioned nrows x ncols.

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

  filename = "stn_latlon.dat"

; Calculate the number of columns.
  ncols = numAsciiCol(filename)

; Given the # of columns, we can use readAsciiTable to read this file.
  data = readAsciiTable(filename,ncols,"float",0)

  nrows = dimsizes(data(:,0))    ; calculate # of rows

  print("'" + filename + "' has " + nrows + " rows and " + ncols + \
        " columns of data.")

xuebiz

下载要扣5个金币哦，不舍得，差不多就是这样子吧，(*^__^*) 嘻嘻……
下次可以贴几行数据看看就好  ------  守财奴留


参考下：
以下来自  ------  NCL常用气象格式的数据读取及可视化
（本论坛有电子版下载http://bbs.06climate.com/forum.php?mod=viewthread&tid=10831）

2.2文本数据

使用下面的 NCL 代码读取刚才生成的 test.txt 文件：
; NCL 脚本
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
begin
  fiASC = "./test.txt"
  nrow = numAsciiRow(fiASC)
  ncol = numAsciiCol(fiASC)
  data = asciiread(fiASC, (/nrow, ncol/), "integer")
  print(data)
end

zly4814624

xuebiz 发表于 2016-8-29 16:02
http://www.ncl.ucar.edu/Applications/read_ascii.shtml
官网例子~~

谢谢。学习了。

zongqing123

谢谢。学习了，很有用

ztf6863

谢谢楼主分享