用GrADS绘制可旋转椭圆【原创】

mofangbao

未经本人允许，该绘制椭圆脚本拒绝任何转载   
    昨天刚要睡觉看到一个人提问的帖子，说是如何用GrADS绘制椭圆，觉得还挺有意思，因为GrADS里面貌似没有自带这方面的东西，draw mark出来的也只能是圆，还看到他的帖子里面附带了一个脚本，于是我就看了看。
    那个帖子在这里：http://bbs.06climate.com/forum.php?mod=viewthread&tid=14278
    帖子里的脚本写的有点凌乱，看不下去了，于是直接试试效果，画出来是这么个图：


额。。。既然是画图的脚本，就不应该跟经纬度联系，而应该用屏幕坐标，要不然如果是各种剖面图就没办法画了。而且这个图是不是有点。。。

    于是，就决定自己写啦！在维基百科又复习了一遍高中数学关于椭圆的一些知识，感觉差不多了就动手写，下面就是我写的脚本画出来的东西：


    这才是画图嘛，这样不管是什么坐标系，什么维度设置，都能够在想要的地方画图椭圆（别告诉我你没发现这图里面有椭圆。。。）
    下载下面的drawEllipse脚本，放到GrADS安装目录的Contents/Scripts目录（有的版本是lib目录）
    话说这图可不是ps的，感兴趣的可以自己下载脚本试试吧！下面是脚本的主要用法：

1. drawEllipse -arg arg
   
2. 其中，-arg表示要设置的参数，arg表示参数的内容
   
3. 具体如下：
   
4. -a 设置椭圆的半长轴，单位英寸，默认1
   
5. -b 设置椭圆的半短轴，单位英寸，默认2
   
6. -cx 设置椭圆的中心x坐标，页面坐标，默认5
   
7. -cy 设置椭圆的中心y坐标，页面坐标，默认5
   
8. -type 设置椭圆的绘制方式，-1表示填色，0表示不填色，1/2/3/4/5分别表示用1-5的对应的draw mark的符号来绘制，默认0
   
9. -thick 表示绘制线条或者mark的粗细，默认1
   
10. -col 表示用来绘制（填色）的颜色，默认为3
    
11. -theta 表示将椭圆逆时针旋转多少度，单位是角度，如45，表示将椭圆逆时针旋转45度
    
12. -accur 表示绘制的精细程度，默认为1，该数值越大，绘制耗时越长，建议为1
    
13. 直接输入drawEllipse可以查看用法
    

drawEllipse -arg arg<br /> 其中，-arg表示要设置的参数，arg表示参数的内容<br /> 具体如下：<br /> -a 设置椭圆的半长轴，单位英寸，默认1<br /> -b 设置椭圆的半短轴，单位英寸，默认2<br /> -cx 设置椭圆的中心x坐标，页面坐标，默认5<br /> -cy 设置椭圆的中心y坐标，页面坐标，默认5<br /> -type 设置椭圆的绘制方式，-1表示填色，0表示不填色，1/2/3/4/5分别表示用1-5的对应的draw mark的符号来绘制，默认0<br /> -thick 表示绘制线条或者mark的粗细，默认1<br /> -col 表示用来绘制（填色）的颜色，默认为3<br /> -theta 表示将椭圆逆时针旋转多少度，单位是角度，如45，表示将椭圆逆时针旋转45度<br /> -accur 表示绘制的精细程度，默认为1，该数值越大，绘制耗时越长，建议为1<br /> 直接输入drawEllipse可以查看用法<br />

下面这个脚本就是用来绘制上面那幅图的（第二幅）：

1. 'reinit'
   
2. drawFlo(1,2,7,6,7,5)
   
3. drawFlo(1,6,7,2,7,5)
   
4. drawFlo(1,10,7,8,7,5)
   
5. drawRot(1,1.5,2,3,1,45,1)
   
6. drawRot(1,4.3,2,2,-1,45,0)
   
7. drawRot(1,7,2,7,-1,45,0)
   
8. drawRot(1,9.6,2,7,2,45,1)
   
9. drawFlo(0.5,3.8,5.5,6,7,7)
   
10. drawFlo(0.5,7.8,5.5,2,7,7)
    
11. ;
    
12. 
    
13. function drawRot(size,px,py,col,type,delta,round)
    
14. a=1*size
    
15. b=0.5*size
    
16. tmp='-a '%a%' -b '%b%' -cx '%px%' -cy '%py%' -type '%type%' -col '%col
    
17. i = 0
    
18. while(i<=180)
    
19.     args=tmp%' -theta '%i
    
20.     'drawEllipse.gs 'args''
    
21.     i=i+delta
    
22. endwhile
    
23. if(round=1)
    
24.     args='-a '%a%' -b '%a%' -cx '%px%' -cy '%py%' -type 0 -thick 3 -col '%col
    
25.     'drawEllipse.gs 'args''
    
26. endif
    
27. return
    
28. 
    
29. function drawFlo(size,px,py,col1,col2,col3)
    
30. a=1*size
    
31. b=0.65*size
    
32. tmp='-a '%a%' -b '%b
    
33. args=tmp%' -col '%col1%' -type -1 -theta 60'%' -cx '%px%' -cy '%py
    
34. 'drawEllipse.gs 'args''
    
35. args=tmp%' -col '%col2%' -type 0 -theta 60 -thick 2'%' -cx '%px%' -cy '%py
    
36. 'drawEllipse.gs 'args''
    
37. args=tmp%' -col '%col1%' -type -1 -theta 120'%' -cx '%(px-a)%' -cy '%py
    
38. 'drawEllipse.gs 'args''
    
39. args=tmp%' -col '%col2%' -type 0 -theta 120 -thick 2'%' -cx '%(px-a)%' -cy '%py
    
40. 'drawEllipse.gs 'args''
    
41. tmp='-a '%(a*0.66)%' -b '%(b*0.7)
    
42. args=tmp%' -col '%col1%' -type -1 -theta 90 -thick 2'%' -cx '%(px-a*0.5)%' -cy '%(py-b*0.4)
    
43. 'drawEllipse.gs 'args''
    
44. xl=px-b-0.05*size
    
45. yt=py-a+0.1*size
    
46. xr=px-b+0.35*size
    
47. yb=yt-2*size
    
48. 'set line 'col3
    
49. 'draw recf 'xl' 'yb' 'xr' 'yt
    
50. return
    
51. 
    

'reinit'<br /> drawFlo(1,2,7,6,7,5)<br /> drawFlo(1,6,7,2,7,5)<br /> drawFlo(1,10,7,8,7,5)<br /> drawRot(1,1.5,2,3,1,45,1)<br /> drawRot(1,4.3,2,2,-1,45,0)<br /> drawRot(1,7,2,7,-1,45,0)<br /> drawRot(1,9.6,2,7,2,45,1)<br /> drawFlo(0.5,3.8,5.5,6,7,7)<br /> drawFlo(0.5,7.8,5.5,2,7,7)<br /> ;<br /> <br /> function drawRot(size,px,py,col,type,delta,round)<br /> a=1*size<br /> b=0.5*size<br /> tmp='-a '%a%' -b '%b%' -cx '%px%' -cy '%py%' -type '%type%' -col '%col<br /> i = 0<br /> while(i<=180)<br />     args=tmp%' -theta '%i<br />     'drawEllipse.gs 'args''<br />     i=i+delta<br /> endwhile<br /> if(round=1)<br />     args='-a '%a%' -b '%a%' -cx '%px%' -cy '%py%' -type 0 -thick 3 -col '%col<br />     'drawEllipse.gs 'args''<br /> endif<br /> return<br /> <br /> function drawFlo(size,px,py,col1,col2,col3)<br /> a=1*size<br /> b=0.65*size<br /> tmp='-a '%a%' -b '%b<br /> args=tmp%' -col '%col1%' -type -1 -theta 60'%' -cx '%px%' -cy '%py<br /> 'drawEllipse.gs 'args''<br /> args=tmp%' -col '%col2%' -type 0 -theta 60 -thick 2'%' -cx '%px%' -cy '%py<br /> 'drawEllipse.gs 'args''<br /> args=tmp%' -col '%col1%' -type -1 -theta 120'%' -cx '%(px-a)%' -cy '%py<br /> 'drawEllipse.gs 'args''<br /> args=tmp%' -col '%col2%' -type 0 -theta 120 -thick 2'%' -cx '%(px-a)%' -cy '%py<br /> 'drawEllipse.gs 'args''<br /> tmp='-a '%(a*0.66)%' -b '%(b*0.7)<br /> args=tmp%' -col '%col1%' -type -1 -theta 90 -thick 2'%' -cx '%(px-a*0.5)%' -cy '%(py-b*0.4)<br /> 'drawEllipse.gs 'args''<br /> xl=px-b-0.05*size<br /> yt=py-a+0.1*size<br /> xr=px-b+0.35*size<br /> yb=yt-2*size<br /> 'set line 'col3<br /> 'draw recf 'xl' 'yb' 'xr' 'yt<br /> return<br /> <br />

绘制椭圆脚本下载（下载后请放到GrADS的脚本目录下）：谢谢@银狐的提醒，文件名里多了一个空格，请自己去掉~ drawEllipse .gs (3.06 KB, 下载次数: 125)

2013-5-10 19:16 上传

点击文件名下载附件
测试脚本下载（上面脚本放好后，直接运行该脚本即可，无需打开任何数据。如果找不到GrADS脚本目录，可以使用绝对路径来执行drawEllipse.gs脚本）
testDraw.gs (1.37 KB, 下载次数: 50)

2013-5-10 18:03 上传

点击文件名下载附件

海豚∩_∩浩

很牛拜啊，学习了

随缘

又出极品！膜拜！

topmad

明天周末啦....画个小图娱乐一下

天目神眉

哈哈，膜拜一下大大，

2010fuqiang

~~~

mofangbao

topmad 发表于 2013-5-10 18:48

                               
登录/注册后可看大图

明天周末啦....画个小图娱乐一下

我还真是这么想的

topmad

mofangbao 发表于 2013-5-10 19:23

                               
登录/注册后可看大图

我还真是这么想的

刚刚画出来了一个想了好几天的图...高兴一下
明天去周游一遭

river

原创好帖啊，除了能画好看的图，还很有参考价值啊

kongfeng0824

挺有意思的！！！