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

13top

感谢分享

lhjou

您好，您的 “[作图] （更新）最最最准确的中国地图（附带南海，精确到市界）” 这篇帖子对我受益匪浅，非常感谢有您这样的前辈给与我们后辈这些支持。
我有一问，希望能听取您的见解。我现在需要将地级市的界限显示出来，您2017.3.2 更新的市界地图与您帖子中的图不符，反而少了市里县级市的分界，比如我所需要的辽宁省-铁岭市-开原市的分界就没了。之后我也想了办法，在您的原始代码中添加了我自己寻找的辽宁省县级市的分界，但是加入之后并不能显示，代码是不报错的，这是否与.shp文件的具体内容有关呢？

尽头的尽头

lhjou 发表于 2019-11-19 14:29
您好，您的 “[作图] （更新）最最最准确的中国地图（附带南海，精确到市界）” 这篇帖子对我受益匪浅，非 ...

你完全可以用自己的底图替换附件里面的

尽头的尽头

lhjou 发表于 2019-11-19 14:29
您好，您的 “[作图] （更新）最最最准确的中国地图（附带南海，精确到市界）” 这篇帖子对我受益匪浅，非 ...

你完全可以用自己的底图替换附件里面的

flyapple1992

谢谢楼主，试一试

飞雪残剑

太棒啦，感谢楼主！

gleamlsc

这个对于写paper画图来说真的很有用！

方莹红

在画大范围地图的时候怎样才可以做到把中国的省分界线画出来但不管其他国家啊

尽头的尽头

方莹红 发表于 2020-3-11 17:09
在画大范围地图的时候怎样才可以做到把中国的省分界线画出来但不管其他国家啊

不画世界地图，就画中国区域的地图就好了

方莹红

尽头的尽头 发表于 2020-3-12 13:29
不画世界地图，就画中国区域的地图就好了

怎样修改呢？我试了一下这篇帖子分享的方法，发现报错了，也不知道是什么原因了，您能帮我看看吗？
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
  load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/diagnostics_cam.ncl"
  begin
  ;================read
f1= addfile("d:/changshui/1800.pressure.nc","r")
f2= addfile("d:/changshui/1800uv.nc","r")
f3= addfile("d:/changshui/1800.2m.t.nc","r")

sp=short2flt(f1->sp(28,:,:))
sp1=sp
sp1=sp/100
t2m=short2flt(f3->t2m(28,:,:))
t2m_1=t2m
t2m_1=t2m-273.15
u10=short2flt(f2->u10(28,:,:))
v10=short2flt(f2->v10(28,:,:))
;printVarSummary(t2m)
time=f1->time
tim= cd_calendar(time,-1)
;print(tim)
;printVarSummary(hgt)
;exit
;printMinMax (t2m, False)
;=================================================================plot
wks = gsn_open_wks("eps","dimian")
res0                   = True
res0@tiMainString     = " dimian"
;res0@gsnLeftString     = "sp"
res0@gsnRightString    = "Wind"
res0@gsnDraw            =False
res0@gsnFrame          =False  
res0@mpFillOn     = False
res0@mpOutlineOn  = True
res0@mpMinLonF         = 80
res0@mpMaxLonF         = 110
res0@mpMinLatF         =10
res0@mpMaxLatF         = 44.75
res0@gsnAddCyclic      =False
res0@cnLineThicknessF  =3
;res0@cnFillOn           = True
res0@cnLevelSelectionMode = "ManualLevels"
;res0@cnLineLabelInterval   = 1
res0@cnLevelSpacingF     =20


fill_area = (/"China: states"/)
res0@mpDateBaseVesion = "MediumRes"
res0@mpDateSetMame = "Earth..4"
res0@mpOutlineOn = True
res0@mpOutlineSpecifiers = fill_area

plot=gsn_csm_contour_map_ce(wks,sp1,res0)

;===============================================================res
res                      = True               ; plot mods desired
res@gsnDraw            =False
res@gsnFrame          =False
res@vcRefMagnitudeF           = 6.5             ; define vector ref mag
res@vcRefLengthF              =  0.035   ; define length of vec
res@vcRefAnnoOrthogonalPosF   = -1.2  
res@vcRefAnnoArrowLineColor   = "purple"         ; change ref vector color
res@vcRefAnnoArrowUseVecColor = False           ; don't use vec color for ref
res@vcMinDistanceF           = 0.03          ; thin vectors
res@vcLineArrowColor        = "purple"           ; change vector color
res@vcLineArrowThicknessF   =2.0              ; change vector thickness
res@vcLineArrowHeadMaxSizeF=0.01
res@vcRefAnnoString2On     =True
res@vcRefAnnoString2       ="m/s"

res@gsnLeftString     = ""
res@gsnRightString    = ""
plot1=gsn_csm_vector(wks,u10,v10,res)
overlay(plot,plot1)

res2     =True
res2@gsnDraw  =False
res2@gsnFrame  =False
res2@cnFillOn =True
res2@cnLineThicknessF  =1.5
res2@gsnRightString  ="~F35~J~F~C"
res2@cnLevelSelectionMode = "ManualLevels"
res2@cnLevelSpacingF     = 4
res2@cnFillColor ="red"
plot2=gsn_csm_contour(wks, t2m_1, res2)
overlay(plot, plot2)

draw(plot)    ; This will draw everything
frame(wks)

               
end