MeteoInfoLab脚本示例：旋转格点场

MeteoInfo

这个例子演示了2维水平格点场的旋转。脚本主要有以下几个步骤：

1、用arange()函数创建格点长的经纬度坐标，再用meshgrid()函数将lon, lat数组变为2维，此时创建的格点场是东西、南北方向的“正矩形”。

2、旋转格点场：利用了矩阵旋转方法（具体可以网上搜索了解）。
2.1 构建一个包含旋转角度正、余玄的2*2矩阵（也就是2维数组）。
2.2 创建一个包含所有格点经纬度坐标的矩阵（xy = zeros([n, 2]），并将经纬度值赋入。由于是希望格点场以左下角为中心旋转，构建的矩阵的左下角要平移到（0, 0）。
2.3 利用矩阵乘积函数dot()将2.1和2.2中构建的两个矩阵相乘即可得到旋转后的矩阵。旋转后的矩阵需要进行和2.2中相反方向的平移。

脚本程序：

1. #Plot
   
2. axesm(tickfontsize=12)
   
3. geoshow('country')
   
4. 
   
5. #Original matrix
   
6. lat = arange(10, 41, 5)
   
7. lon = arange(70, 131, 5)
   
8. lon, lat = meshgrid(lon, lat)
   
9. scatterm(lon, lat, facecolor='b', size=8, edge=False)
   
10. 
    
11. #Rotate matrix
    
12. angle = 20 * pi / 180
    
13. rotate = array([[cos(angle),sin(angle)], [-sin(angle),cos(angle)]])
    
14. n = lon.size
    
15. xy = zeros([n, 2])
    
16. xy[:,0] = lon - 70
    
17. xy[:,1] = lat - 10
    
18. xy = dot(xy, rotate)    #Matrix multiplication
    
19. xy[:,0] = xy[:,0] + 70
    
20. xy[:,1] = xy[:,1] + 10
    
21. scatterm(xy[:,0], xy[:,1], facecolor='r', size=8, edge=False)
    
22. 
    
23. xticks(arange(0, 361, 10))
    
24. yticks(arange(-90, 91, 10))
    
25. title('Matrix rotate sample')
    
26. xlim(50, 150)
    
27. ylim(0, 70)

#Plot<br /> axesm(tickfontsize=12)<br /> geoshow('country')<br /> <br /> #Original matrix<br /> lat = arange(10, 41, 5)<br /> lon = arange(70, 131, 5)<br /> lon, lat = meshgrid(lon, lat)<br /> scatterm(lon, lat, facecolor='b', size=8, edge=False)<br /> <br /> #Rotate matrix<br /> angle = 20 * pi / 180<br /> rotate = array([[cos(angle),sin(angle)], [-sin(angle),cos(angle)]])<br /> n = lon.size<br /> xy = zeros([n, 2])<br /> xy[:,0] = lon - 70<br /> xy[:,1] = lat - 10<br /> xy = dot(xy, rotate)    #Matrix multiplication<br /> xy[:,0] = xy[:,0] + 70<br /> xy[:,1] = xy[:,1] + 10<br /> scatterm(xy[:,0], xy[:,1], facecolor='r', size=8, edge=False)<br /> <br /> xticks(arange(0, 361, 10))<br /> yticks(arange(-90, 91, 10))<br /> title('Matrix rotate sample')<br /> xlim(50, 150)<br /> ylim(0, 70)

aas314

{:5_213:}