现在正在学NCL，用它绘制了一样中国地图

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin

;---Area to zoom in on.

minlat = 17

maxlat = 55

minlon = 72

maxlon = 136

wks = gsn_open_wks("png","China_Map")

mpres = True

mpres@gsnMaximize = True

mpres@gsnDraw = False

mpres@gsnFrame = False

mpres@mpDataBaseVersion = "LowRes"

mpres@mpOutlineOn = False

mpres@mpLandFillColor = "white"

mpres@mpOceanFillColor = "white"

mpres@mpInlandWaterFillColor = "white"

;---Zoom in on area of interest

mpres@mpLimitMode = "LatLon"

mpres@mpMinLatF = minlat

mpres@mpMaxLatF = maxlat

mpres@mpMinLonF = minlon

mpres@mpMaxLonF = maxlon

map_plot = gsn_csm_map(wks,mpres) ; Create map, but don't draw it yet.

;****************************************************************************

; sections for shapefiles

;****************************************************************************

;---适合中国的国界地图

shpfn1="bou1_4l.shp"

lnres = True

lnres@minlat = minlat

lnres@maxlat = maxlat

lnres@minlon = minlon

lnres@maxlon = maxlon

shp_plot1 = gsn_add_shapefile_polylines(wks,map_plot,shpfn1,lnres)

;---加上省级区域地图

shpfn2="bou2_4p.shp"

gres = True

gres@minlat = minlat

gres@maxlat = maxlat

gres@minlon = minlon

gres@maxlon = maxlon

gres@gsEdgesOn =True

gres@gsEdgeColor ="blue"

shp_plot2 = gsn_add_shapefile_polygons(wks,map_plot,shpfn2,gres)

;---加上一级河流

shpfn3="hyd1_4l.shp"

lnres3 = True

lnres3@minlat = minlat

lnres3@maxlat = maxlat

lnres3@minlon = minlon

lnres3@maxlon = maxlon

shp_plot3 = gsn_add_shapefile_polylines(wks,map_plot,shpfn3,lnres3)

maximize_output(wks,False)

end效果如下：