圖1 假設的全球某成分分布
圖2 北半球和南半球球面立體投影座標系
以球面立體投影開始,這個在大氣痕量成分分布產品中應用得非常多,能清楚明了地看到世界各地的某高度的濃度或總量分布。
二、解壓,然後通過matlab設定路徑包含m_map
三、投影座標系繪製
1,首先定義投影方式
m_proj('stereographic','latitude',90,'radius',90,'rotagnle',0);
;投影名稱為stereographic
;中心緯度為北緯90度
;繪製緯度長度為90度,即從n90-0,如果radius為55的話,則繪製n90-n55的區域
;座標系旋轉角度為0,即垂直方向為經度0度
2,繪製海岸線
m_coast('patch',[.7 .7 .7], 'edgecolor', 'none');
;填充區域顏色為[.7 .7 .7],三個分量分別代表r,g,b,以小數表示,即255為1,0為0
;邊界無顏色
3,繪製邊界、座標軸等
m_grid;
四、繪製目標資料
1,資料準備
例如:臭氧ozone = [360*90],即有經度-180:180和緯度0:90的臭氧資料,即北半球資料
經度資料lon=0:180
緯度資料lat=0:90
將座標其轉化為網格
[lon,lat] = meshgrid(lat, lon); 注意這裡的lat和lon的順序
現在lon = [360*90],lat = [360*90];
2,座標轉換
將經緯度座標轉換成投影後座標
[x,y] = m_ll2xy(lon,lat);
3,繪製目標資料
hold on; 一定不要忘記此句,否則前面的將全部覆蓋
pcolor(lon,lat,ozone);
shading interp;
box on;
效果見圖1.
幾點注意的地方:
1,座標轉換時,給定的經緯度一定是投影範圍內的點,否則得不到有效值。
2,緯度的範圍是0-90,北緯為正,南緯為負。即-55表示南緯55度。
3,經度的範圍是0-180度,東經為正,西經為負。
4,繪製南半球時,設定中心為s90,即-90。要想0也在豎起時,要設定旋轉方向為180。
m_proj('stereographic','latitude',-90,'radius',90,'rotagnle',180);
5,繪製南半球時,要通過m_grid設定經度軸標在最外層,否則都聚焦在最裡。
m_grid('xaxislocation','top');
5、繪製目標資料之前,加上hold on;否則將繪製的座標系全部覆蓋了。
mmap 用mmap實現cp命令
mmap 可以把磁碟檔案的一部分直接對映到記憶體,這樣檔案中的位置直接就有對 應的記憶體位址,對檔案的讀寫可以直接用指標來做而不需要read write 函式。include void mmap void addr,size t len,int prot,int flag,int filedes,o...
mmap例子詳解
mmap函式是unix linux下的系統呼叫,來看 unix netword programming 卷二12.2節有詳細介紹。mmap系統呼叫並不是完全為了用於共享記憶體而設計的。它本身提供了不同於一般對普通檔案的訪問方式,程序可以像讀寫記憶體一樣對普通檔案的操作。而posix或系統v的共享記憶...
mmap詳細解釋
mmap系統呼叫並不是完全為了用於共享記憶體而設計的。它本身提供了不同於一般對普通檔案的訪問方式,程序可以像讀寫記憶體一樣對普通檔案的操作。而posix或系統v的共享記憶體ipc則純粹用於共享目的,當然mmap 實現共享記憶體也是其主要應用之一。mmap系統呼叫使得程序之間通過對映同乙個普通檔案實現...