問題
利用GMT繪制小區(qū)域圖件時，經(jīng)緯度不利于表達圖內(nèi)部要素的相對方位沪么，但是可以很直觀的表達圖所處的絕對位置上忍。UTM平面坐標則可以直接表達圖內(nèi)的位置關(guān)系。那么如何才能讓經(jīng)緯度和UTM同時顯示在一個圖件中诱咏？達到以下這樣的效果，如圖示：
- 上缴挖、右軸顯示經(jīng)緯度標注
-
左袋狞、下軸顯示UTM投影平面坐標
目標圖樣
思路
這張圖的本質(zhì)是要用投影平面坐標系作為繪圖邊框，而非常規(guī)GMT的經(jīng)緯度坐標映屋。投影類型應(yīng)該選用UTM投影的-Ju苟鸯，但是該投影并不支持標注平面坐標。并且經(jīng)實驗后發(fā)現(xiàn)棚点，使用該投影后邊框有兩種不同形式：
- 邊框為帶弧度的地理坐標系
-R后接常規(guī)的經(jīng)緯度區(qū)域范圍早处，如-R103/104/33/34時，那么邊框就是地理坐標系下的帶一定弧度的邊框瘫析。 - 邊框為橫平豎直的平面坐標系
-R后如果用網(wǎng)格數(shù)據(jù)指范圍砌梆，且網(wǎng)格數(shù)據(jù)格式為UTM時，如-Rutm.grd贬循，那么邊框就是橫平豎直的平面坐標系了咸包。
了解這些后，我們應(yīng)這么處理：
- 左甘有、下軸
采用線性投影-Jx诉儒，讓左、下軸顯示UTM平面坐標亏掀。 - 右、上軸
右泛释、上軸采用-Ju投影滤愕，并和-Jx保持相同的比例尺。問題在于邊框怜校，這里圖示一下上述的地理坐標系邊框和平面坐標系邊框的區(qū)別：-
-R103/104/33/34
會導(dǎo)致右间影、上軸的為經(jīng)緯度邊框，這樣和左茄茁、下軸的平面坐標邊框是不匹配的魂贬。見代碼和示意圖：
-
gmt begin test pdf,png
gmt basemap -Jx1:500000 -Rutm.grd -Bag -BSWrt --MAP_FRAME_PEN=1p,blue --MAP_GRID_CROSS_SIZE_PRIMARY=5p
#gmt basemap -Ju48S/1:500000 -Rutm.grd -Bag -BNElb --MAP_FRAME_PEN=0.4p,red
gmt basemap -Ju48S/1:500000 -R103/104/33/34 -Bag -BNEls --MAP_FRAME_PEN=0.4p,red
gmt end
紅色邊框為使用-R103/104/33/34的結(jié)果巩割，紅線邊框為經(jīng)緯度弧形邊框，和藍線的平面坐標邊框有較大差距
-
-Rutm.grd
這樣右付燥、上軸是平面坐標系的邊框宣谈，和左、下軸匹配键科，見代碼和示意圖闻丑。
gmt begin test pdf,png
gmt basemap -Jx1:500000 -Rutm.grd -Bag -BSWrt --MAP_FRAME_PEN=1p,blue --MAP_GRID_CROSS_SIZE_PRIMARY=5p
gmt basemap -Ju48S/1:500000 -Rutm.grd -Bag -BNElb --MAP_FRAME_PEN=0.4p,red
#gmt basemap -Ju48S/1:500000 -R103/104/33/34 -Ba -BNEls --MAP_FRAME_PEN=0.4p,red
gmt end
-R使用平面坐標，結(jié)果紅線網(wǎng)格和藍線網(wǎng)格完全吻合
- 如何獲取UTM格式的網(wǎng)格數(shù)據(jù)
一般下載得到的數(shù)據(jù)都是基于地理坐標的勋颖，但grdproject可以對網(wǎng)格數(shù)據(jù)做地圖變換和逆變換嗦嗡。代碼如下：
gmt grdproject lon-lat.grd -Gutm.grd -Ju48S/1:1 -C -F
后期繪圖
如果是經(jīng)緯度數(shù)據(jù)，就將其放在-Ju那行代碼下面進行相關(guān)繪制饭玲。
如果是UTM坐標數(shù)據(jù)侥祭，就將其放在-Jx那行代碼后面進行繪制。
