m map投影 科學網 基於m

2021-10-12 14:39:54 字數 1394 閱讀 2652

圖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系統呼叫使得程序之間通過對映同乙個普通檔案實現...