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

lvchabj

lvchabj 发表于 2013-7-25 16:32

                               
登录/注册后可看大图

嗯 ，是的，全是压缩包里面自带的。

您好，我按照上面我注释的，后面画了观测降水图 。但是这个里面有两个国界线

lvchabj

绘制的图片

lvchabj

云南那块儿 好像有问题？不知道是不是这样，而且我在后面加了省界线，也没有显示
wks_type = "png"
  wks = gsn_open_wks(wks_type,argu(1))
  gsn_define_colormap(wks,"precip4_11lev")
;  gsn_draw_colormap(wks)

  res = True
  res@gsnAddCyclic  = False
  res@mpDataSetName         = "Earth..4"
  res@cnFillDrawOrder         = "PreDraw"

  res@mpDataBaseVersion     = "MediumRes"
;  res@mpOutlineOn           = False
res@mpOutlineOn           = True
res@mpOutlineColor = 0
; res@mpOutlineSpecifiers   = (/"China:states","Taiwan"/)
  res@mpOutlineSpecifiers   = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)


  res@mpMinLatF             =  17          ; Asia limits
  res@mpMaxLatF             =  53
  res@mpMinLonF             =  75
  res@mpMaxLonF             = 135


  res@mpGeophysicalLineThicknessF= 1.
  res@mpNationalLineThicknessF= 1.

  res@mpProjection = "LambertConformal"   ;?录虏庐???露?掳
  res@mpLambertMeridianF = 105
  res@mpLimitMode = "LatLon"
  res@mpLambertParallel1F = 29
  res@mpLambertParallel2F = 29
res@mpAreaMaskingOn = True
; res@mpMaskAreaSpecifiers = (/"China:states","Taiwan"/)
  res@mpMaskAreaSpecifiers = (/"China","Taiwan","Disputed area between India and China","India:Arunachal Pradesh"/)
res@mpLandFillColor         = "white"

  res@mpOceanFillColor = 0
  res@mpInlandWaterFillColor = 0

  res@cnFillOn      = True
  res@cnLinesOn     = False
  res@cnLineLabelsOn = False
  res@cnFillDrawOrder = "PreDraw"

  res@cnLevelSelectionMode = "ExplicitLevels"

  res@cnLevels   =  (/0.1, 10.0, 25., 50.0, 100.0 ,250.0/)

   res@lbLabelBarOn = True
   res@lbLabelStrings = (/"0.1","10","25","50","100","250"/)
   res@lbOrientation =   "vertical"          ; vertical label bar

   res@tiMainFont            = "helvetica"
  res@tiMainOffsetYF        = 0.02  ;set place for main title along Y,offset
  res@tiMainFontHeightF     = 0.02   ;set main title font size

  res@tiMainString = argu(1)
;  res@tiMainString          = argu(0)

  plot = gsn_csm_contour_map(wks,data1,res)

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
cnres = True
cnres@china = True
cnres@river = True
cnres@province = True
cnres@nanhai = True
cnres@diqu = False
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
  
; plot = gsn_csm_contour_map(wks,data1,res)

      chinamap = add_china_map(wks,plot,cnres)
  ; ADD Yellow River  and Yangtze River
  ; res for following lines
      resline = True
      resline@gsLineThicknessF = 2.0           ; thickness of lines
      resline@gsLineColor = 6

    ; Yangtze River
      Yangtze = asciiread("Yangtze.ascii",(/73,2/),"float")
      dum1 = new(dimsizes(Yangtze(:,0)),graphic)
      n=dimsizes(Yangtze(:,0))
      do i = 0 ,n-2
         dum1(i) = gsn_add_polyline(wks,plot,Yangtze(i:i+1,0),Yangtze(i:i+1,1),resline)
      end do
    ; Yellow River
      Yellow = asciiread("Yellow.ascii",(/75,2/),"float")
      dum2 = new(dimsizes(Yellow(:,0)),graphic)
      n2=dimsizes(Yellow(:,0))
      do j = 0 ,n2-2         
         ; dum2(j) = gsn_add_polyline(wks,plot,Yellow(j:j+1,0),Yellow(j:j+1,1),resline)
          dum2(j) = gsn_add_polyline(wks,plot,Yellow(j:j+1,0),Yellow(j:j+1,1),resline)
      end do

lvchabj

  ;   chinamap = gsn_add_china_map(wks,plot,cnres)

      draw(plot)
      frame(wks)


end

尽头的尽头

lvchabj 发表于 2013-7-25 18:05

                               
登录/注册后可看大图

;   chinamap = gsn_add_china_map(wks,plot,cnres)

      draw(plot)

将res@mpOutlineOn  = True改为Fasle，NCL自带的国界线不被我们国家承认，在藏南争议地区

lvchabj

我改过这个设置，变成fasle后就没有边界了，不知道是什么原因，是不是还是跟我前面注释有关？
因为虽然运行后，有图形还是会报警告错误。

lvchabj

还有我觉得这个方法绘制的图形，在云南那块儿是不是对了一点出来啊？

东篱

收了，学习学习

xiao^^@n

w我太佩服楼主了，真是偶像

xiao^^@n

只是我水平不行，发不了贴，