对做任意两点的时间斜剖面的一些想法

chnlhttom

首先是基于先找点，再写入新的文件，最后画图的思想。
首先一个时次，计算两点之间多个点的经纬度，然后再写入二进制文件，在画图。下面是我的gs和ctl文件以及一张实验图，不知道对不对。

1. 'reinit'
   
2. 'sdfopen E:\uwnd.nc'
   
3. 'set gxout fwrite '
   
4. 'set fwrite E:\uwnd.dat'
   
5. 'set grads off'
   
6. 'set zlog on'
   
7. 'set csmooth on'
   
8. 
   
9. 'set t 1 20'
   
10. lon1 = 70.0
    
11. lon2 = 120.0
    
12. lat1 = 50.0
    
13. lat2 = 25.0
    
14. lon = lon1
    
15. 
    
16. while (lon <= lon2)
    
17.   lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)
    
18. 'set lat 'lat''
    
19. 'set lon 'lon''
    
20. 'a=u'
    
21. 'set y 1'
    
22. 'set lon 'lon''
    
23. 'd a'
    
24.   lon = lon + 0.5
    
25. 
    
26. endwhile
    
27. 'disable  fwrite '
    
28. 'c'
    
29. 'reinit '
    
30. ;

复制代码

1. dset F:\Obliqueprofile\uwnd.dat
   
2. undef 9999
   
3. title grid data
   
4. xdef 101 linear 60 0.5
   
5. ydef 1 linear 20 0.5
   
6. zdef 1 linear 1000 10
   
7. tdef 20 linear jan2003 1yr
   
8. vars 1
   
9. xie 0 99 xie data
   
10. endvars

复制代码

river

'set y 1'
'set lon 'lon''
纬向为什么是固定到1，'set y 1'为什么就代表了北纬20度？

chnlhttom

river 发表于 2017-4-23 09:12
'set y 1'
'set lon 'lon''
纬向为什么是固定到1，'set y 1'为什么就代表了北纬20度？

读取的时候是读取的任意两点的一条斜线，写入时我就想写成一条水平线，这样ctl我才想的出怎么写，所以你明白我的那种意思吗，这样想能实现画图，但是我不知道对不对

river

chnlhttom 发表于 2017-4-23 15:27
读取的时候是读取的任意两点的一条斜线，写入时我就想写成一条水平线，这样ctl我才想的出怎么写，所以你 ...

我觉得你的思路是对的，但是set y 1不能代表set lat 20 吧，你直接用set lat 20看看结果一样不？

oydl1985

先占个座，看看大神怎么说！

18393810639

期待大神回复

jolincai

楼主在吗？是沿一条斜线？那横坐标应该同时标注经度和纬度啊？不然图看上去就是沿经或纬线的传播。不知道脚本是否正确 想用用

jolincai

while (lon <= lon2)

  lat = lat1 + (lat2-lat1)*(lon-lon1) / (lon2-lon1)

'set lat 'lat''

'set lon 'lon''（斜线上等分点的经纬度）

'a=u'

'set y 1'

'set lon 'lon''（可是又set  y 1为什么呢？这个y是否是上一步的lat值呢？如果不是的话那相当于提取沿着y 1的纬线提取数值了啊）

'd a'

  lon = lon + 0.5

chnlhttom

jolincai 发表于 2017-6-8 16:27
while (lon

不好意思，这几天忙着没上。这个坐标轴我还没有解决，你如果需要私信我你qq号我们可以交流一下

jolincai

chnlhttom 发表于 2017-6-14 20:53
不好意思，这几天忙着没上。这个坐标轴我还没有解决，你如果需要私信我你qq号我们可以交流一下

我后来用插值法  得到这条直线上的点 然后画出来了 用你的方法也画过 感觉是错的