NCL如何批量读取txt数据画轨迹!

对方的凤飞飞

学了ncl也有快两年时间了,这是我从来没涉足的问题.单纯的如果读一个文件进行循环照葫芦画瓢我还行.现在是读多个文件,文件行数不一样,列是一样的话,画不同文件的轨迹,我就犯了难了,一下是我个人的困难和思考,希望大家予以指导!
画一个文件的轨迹图

1. 
   
2. begin
   
3. 
   
4.    wks = gsn_open_wks("png","traj1")               ; send graphics to PNG file
   
5.    data=asciiread("1.txt", (/38,13/), "float")
   
6.    lon  = data(:,2)
   
7.    lat  = data(:,3)
   
8. 
   
9.   
   
10.    res             = True                              ; plot mods desired
    
11.    res@gsnFrame    = False                             ; don't advance yet
    
12.    res@gsnMaximize = True
    
13.    res@gsnDraw          = False   ; so we can add poly stuff
    
14.    res@mpDataBaseVersion     =  "MediumRes" ;
    
15.    res@mpPerimOn             = False          ;-- don't draw the box around the plot
    
16.    res@gsnPolar="SH"
    
17.    res@mpCenterLonF          = 0              ;-- center at lon=10
    
18.    res@mpMaxLatF          = -30   
    
19.    res@mpGridLineColor       = "grey30"         ;-- set grid line color
    
20.    res@tiMainString    = ""  ;
    
21.    res@gsnCenterString = ""   ; center string
    
22.    res@mpFillOn               = False    ; 打开之后线才会显现
    
23.    res@mpOutlineOn            = True    ; turn on continental outline打开后图变清晰了
    
24. 
    
25. 
    
26.    map = gsn_csm_map(wks,res)                        ; gsn_csm_map_ce
    
27.    pres                  = True               ; polyline resources
    
28.    pres@gsLineThicknessF = 2.0                ; line thickness
    
29.   
    
30.    pres@gsLineColor = "red"                        
    
31.    gsn_polyline(wks,map,lon,lat,pres)      ; draw the traj
    
32.    draw(map)                                          
    
33.    frame(wks)                                         
    
34.   
    
35. end
    
36. 
    

复制代码

如果读多个文件呢,行数不一样,如何利用循环进行数据的可视化.目前只找到一个代码

1. fils = systemfunc ("ls *.txt")
   
2. ; print(fils)
   
3. ; print(dimsizes(fils))
   
4. 
   
5. cols=13
   
6. do i=0,dimsizes(fils)-1
   
7.     print(i+"Now reading..."+fils(i))   ;过程监控
   
8. <div class="blockcode"><blockquote>data = asciiread(fils(i),-1,"float")      ;读成一整列
   
9.     rows=dimsizes(data)/cols
   
10.     nfile=onedtond(data,(/rows,cols/))

复制代码

    delete(data)
    delete(nfile)
    delete(rows)

end do
请问中间核心的三行代码是什么意思呢?不是非常的懂这个循环该如何写,思考时间睡觉都在思考!谢谢!

随缘

fils = systemfunc ("ls *.txt")
; print(fils)
; print(dimsizes(fils))
cols=13
do i=0,dimsizes(fils)-1
    print(i+"Now reading..."+fils(i))   ;过程监控
   data = asciiread(fils(i),-1,"float")      ;读成一整列
    rows=dimsizes(data)/cols
    nfile=onedtond(data,(/rows,cols/))

    delete(data)
    delete(nfile)
    delete(rows)
end do
这个写的很清楚了，先把列定了cols=13，，这是所有文件共同的，然后先顺序读成一维数据，然后再通过onedtond函数转成2维的最终数据，方法很好啊，很适合行数不一样,列一样这种情况

oydl1985

列是一样的话,不同文件合并这块建议楼主用Python处理，处理完成后再用NCL绘画呀！或者直接PyNCL就画啦！仅供参考

对方的凤飞飞

oydl1985 发表于 2020-9-25 12:55
列是一样的话,不同文件合并这块建议楼主用Python处理，处理完成后再用NCL绘画呀！或者直接PyNCL就画啦！仅 ...

请问能推荐什么函数呢,我是正有这个意思.

对方的凤飞飞

随缘 发表于 2020-9-25 15:07
fils = systemfunc ("ls *.txt")
; print(fils)
; print(dimsizes(fils))

我在家园找到的，我研究研究，大神，请问按照这样的画图的时候我的数据就存放在了nfile(i)对吧,利用这个数据来写循环对吗?

随缘

对方的凤飞飞 发表于 2020-9-25 15:31
我在家园找到的，我研究研究，大神，请问按照这样的画图的时候我的数据就存放在了nfile(i)对吧,利用这个 ...

是的，每个文件的数据就在nfile里，是一个rows,cols的二维数组

小拉拉

请问楼主，数据列不一样，怎么读成想要的二位数组？

对方的凤飞飞

小拉拉 发表于 2020-10-17 09:56
请问楼主，数据列不一样，怎么读成想要的二位数组？

我就是这个不懂呀

ZJN1998

想问一下楼主大大，追踪的是什么呀，我最近在学习做水汽追踪，可以分享一下lz画轨迹的脚本嘛，感谢

贫道敬孔

print(i+"Now reading..."+fils(i))   ;过程监控这个思路很有意思