（更新）最最最准确的中国地图（附带南海，精确到市界）

LPR

试一下                           

LPR

试一下                                      

LPR

试一下                                      

BestRiven

1. ; These files are loaded by default in NCL V6.2.0 and newer
   
2. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
   
3. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
   
4. load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
   
5. load "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.ncl"
   
6. ;**************************************************************
   
7. begin
   
8. ;***************************************************************
   
9. ; User Input
   
10. ;***************************************************************
    
11. 
    
12.   ;  diri    = ""                            ; input directory
    
13.    ; fili    = ""
    
14. 
    
15.    shpname = "$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/cnmap.shp"
    
16. 
    
17.    cnres           = True
    
18.    cnres[url=home.php?mod=space&uid=46658]@china[/url]     = True       ;draw china map or not
    
19.    cnres[url=home.php?mod=space&uid=15693]@river[/url]     = True       ;draw changjiang&huanghe or not
    
20.    cnres@province  = True       ;draw province boundary or not
    
21.    cnres@nanhai    = True       ;draw nanhai or not
    
22.    cnres@diqu      = False       ; draw diqujie or no
    
23. ; End User Input
    
24. ;***************************************************************
    
25. ; Read
    
26. ;***************************************************************
    
27.    f      = addfile ("/run/media/hua/lsh/HM/results/CCTM_hm12_Linux2_x86_64.JJJ27_12.CGRID.2012183.ncf", "r")
    
28. ;   x      = f->NH3(0, 0, :, :)
    
29.     a1 = f->AECI(:,0,:,:)
    
30.     a3 = f->AECJ(:,0,:,:)
    
31.     a4 = f->ASO4I(:,0,:,:)
    
32.     a5 = f->ASO4J(:,0,:,:)
    
33.     a6 = f->ANO3I(:,0,:,:)
    
34.     a7 = f->ANO3J(:,0,:,:)
    
35.     a8 = f->AORGPAI(:,0,:,:)
    
36.     a9 = f->AORGPAJ(:,0,:,:)
    
37.     a10 = f->A25J(:,0,:,:)
    
38.     a11 = f->AH2OJ(:,0,:,:)
    
39.     a12 = f->ANAJ(:,0,:,:)
    
40.     a13 = f->ACLJ(:,0,:,:)
    
41.     a14 = f->ANH4I(:,0,:,:)
    
42.     a15 = f->ANH4J(:,0,:,:)
    
43.    x1 = a1 + a3 + a4 +a5 +a6 +a7 + a8 + a9 + a10 + a11 + a12 + a13 + a14 + a15
    
44.     x  = dim_avg_n(x1,0)
    
45. ;   x@_FillValue    = -1
    
46. ;   x@missing_value = -1
    
47. 
    
48.    dimx   = dimsizes(x)
    
49. ;*****************************************************
    
50. ; Create netCDF coordinate variables. The lat are N->S
    
51. ;*****************************************************
    
52. 
    
53.     XCENT = f@XCENT
    
54.     XORIG = f@XORIG
    
55.     C1 = XCENT + XORIG/108000                        ; limit bnd
    
56.     C2 = XCENT - XORIG/108000
    
57.     YCENT = f@YCENT
    
58.     YORIG = f@YORIG
    
59.     R1 = YCENT + YORIG/108000
    
60.     R2 = YCENT - YORIG/108000
    
61. 
    
62.    nlat       = dimx(0)
    
63.    mlon       = dimx(1)
    
64.    lat        = fspan(R1, R2, nlat)
    
65.    lon        = fspan(C1, C2, mlon)
    
66. 
    
67.    lat@long_name = "latitude"
    
68.    lat@units  = "degrees_north"
    
69.    lat!0      = "lat"
    
70.    lat&lat    =  lat
    
71. 
    
72.    lon@long_name = "longitude"
    
73.    lon@units  = "degrees_east"
    
74.    lon!0      = "lon"
    
75.    lon&lon    =  lon
    
76. 
    
77.    printMinMax(lat, True)
    
78.    printMinMax(lon, True)
    
79. 
    
80. ;***************************************************************
    
81. ; Associate the spatial coordinates with variables
    
82. ;***************************************************************
    
83.    x!0      = "lat"                   ; 1st ... name the dimensions
    
84.    x!1      = "lon"
    
85.    x&lat    =  lat                    ; create coordinate variable
    
86.    x&lon    =  lon
    
87. 
    
88. 
    
89. 
    
90. ;-- define the workstation (plot type and name)
    
91.     wks = gsn_open_wks("png","1PM2.5" )
    
92. ;-- set resources
    
93. 
    
94.     res              = True
    
95.     res@tiMainString = "1PM2.5"
    
96. ;   res@vpHeightF    =
    
97.    ; res@mpDataBaseVersion = "MediumRes"
    
98.     res@mpProjection        = "CylindricalEquidistant"
    
99.     res@gsnAddCyclic = False
    
100. 
     
101.     res@mpLimitMode                = "LatLon"
     
102.     res@mpMinLatF                  = 0
     
103.     res@mpMaxLatF                  = 60
     
104.     res@mpMinLonF                  = 70
     
105.     res@mpMaxLonF                  = 140
     
106.     ;res@mpCenterLonF               = f@P_GAM
     
107. 
     
108.     res@mpOutlineOn   = False
     
109.     res@gsnDraw       =  False
     
110.     res@gsnFrame      =  False
     
111.     res@gsnMaximize   =  True
     
112. 
     
113. 
     
114.     res@cnFillOn      =  True
     
115.     res@cnFillPalette = "BlAqGrYeOrRe"
     
116.     res@cnLinesOn     = False
     
117.     res@cnLineLabelsOn= False
     
118.    ; res@cnFillMode    = "CellFill"           ; Cell Mode
     
119.     res@cnFillMode    = "RasterFill"         ; Raster Mode
     
120.     ; res@cnFillMode    = "AreaFill"
     
121.     res@cnSmoothingOn       =  True
     
122.     res@cnLevelSelectionMode = "ManualLevels"     ; set manual contour levels
     
123.     res@cnMinLevelValF       = 0               ; set min contour level
     
124.     res@cnMaxLevelValF       = 200                ; set max contour level
     
125.     res@cnLevelSpacingF      = 20     ; set contour spacing
     
126. ;-- draw the contour map
     
127.     plot = gsn_csm_contour_map(wks, x, res)
     
128. 
     
129. ;-- set resources for the polylines
     
130.     plres               =  True
     
131.     plres@gsLineColor   = "black"
     
132.     plres@gsLineThicknessF = 2.5
     
133.     poly                = gsn_add_shapefile_polylines(wks,plot,shpname,plres)
     
134. ;-- draw the plot
     
135.     draw(plot)
     
136.     frame(wks)
     
137. 
     
138. end
     
139. 
     

复制代码

博主，我没理解您的教程，我想直接用shapename调用您叠加好的中国底图该怎么操作啊，您看我脚本那一行shapename该怎么填啊

明明~~

谢谢楼主分享~~

yaotianci

因为分享而强大。

wade

精细化地图数据很适用

smile2me

厉害了厉害了，这个超赞的。一直在找地图终于找到了

taotao

膜拜大神！谢谢分享！

陈牛奶

好棒好棒，太感谢啦~~~