NCL添加地区分界线-shp文件哪里可以下？

Soaring

shp_filename = ("$NCL_TUT/data/Shapefiles/DEU_adm/DEU_adm1.shp")
poly = gsn_add_shapefile_polylines(wks,plot,shp_filename,False)
以上命令是给图形加上地区轮廓分界线的，中国各个省的.shp文件在哪里可以下到呢？
请教各位大神！！

andrewsoong

本论坛有做好的，好像在MI板块

Soaring

andrewsoong 发表于 2017-3-20 12:37
本论坛有做好的，好像在MI板块

我找到啦，谢谢！！！

1649518749

http://bbs.06climate.com/forum.p ... &extra=page%3D1
感觉这个也适合你

Soaring

1649518749 发表于 2017-3-20 15:49
http://bbs.06climate.com/forum.php?mod=viewthread&tid=45252&extra=page%3D1
感觉这个也适合你

好的，我自己先研究研究吧！

Soaring

Soaring 发表于 2017-3-20 15:59
好的，我自己先研究研究吧！

好像不对，我想画北京市具体的图，你这个好像只是为了画中国地图而描绘的每个省的大致轮廓，那每个省具体的shp在哪里呢？

Soaring

andrewsoong 发表于 2017-3-20 12:37
本论坛有做好的，好像在MI板块

你好，还有一个小问题：
begin

  f=addfile("pres.mon.ltm.nc","r")
  pres= f->pres(0,:,:)

;---------------------------------------------------------------
  wks = gsn_open_wks("png","jiangsu")
  gsn_define_colormap( wks ,"gui_default")  

  res                            = True            
  res@gsnMaximize                = True
  res@gsnDraw                    = False
  res@gsnFrame                   = False

;---------------------------------------------------------------
; Mask Beijing plot

; set for the map
  res@mpMinLatF                  = 30
  res@mpMaxLatF                  = 36
  res@mpMinLonF                  = 115
  res@mpMaxLonF                  = 125

  res@mpFillOn                   = True
  res@mpDataSetName              = "Earth..4"
  res@mpDataBaseVersion          = "MediumRes" ; or "Ncarg4_1"

  res@mpAreaMaskingOn            = True
  res@mpOutlineSpecifiers        = "Jiangsu"
  res@mpMaskAreaSpecifiers       = "Jiangsu"

  res@mpLandFillColor            = "white"
  res@mpOceanFillColor           = "white"
  res@mpProvincialLineColor      = "red"
  res@mpProvincialLineThicknessF = 2.

  res@mpShapeMode                = "FreeAspect"
  res@mpFillDrawOrder            = "PostDraw"

;---------------------------------------------------------------
; set for the plot

  res@cnFillOn                   = True
  res@cnFillDrawOrder            = "PreDraw"
  res@cnLinesOn                  = False            
  res@cnLevelSpacingF            = 20.            
  res@gsnSpreadColors            = True   
  res@lbLabelAutoStride          = True
  res@pmTickMarkDisplayMode      = "Always"
  res@gsnRightString             = "hPa"
  res@gsnLeftString              = "Surface pressure in Beijing with NCL Mask"

  map4 = gsn_csm_contour_map(wks,pres,res)

;--- add shp file boundaries to test whether the border overlap---
  ShpDir = "C:/cygwin/usr/local/ncl/lib/ncarg/nclscripts/cnmap/" ;"$NCARG_ROOT/lib/ncarg/nclscripts/cnmap/"
  cnres                       = True
  cnres@gsLineThicknessF      = 1.5
  cnres@gsLineColor           = "blue"

  getvalues map4
    "mpMinLatF" : cnres@minlat
    "mpMaxLatF" : cnres@maxlat
    "mpMinLonF" : cnres@minlon
    "mpMaxLonF" : cnres@maxlon
  end getvalues

  plotshp = gsn_add_shapefile_polylines(wks,map4,ShpDir+"county_2004.shp",cnres)

;---------------------------------------------------------------
  draw(map4)
  frame(wks)

end  


怎么把江苏以外的中国轮廓线去掉呢？
谢谢！

Soaring

andrewsoong 发表于 2017-3-20 12:37
本论坛有做好的，好像在MI板块

已经解决啦！mpOutlineBoundarySets="NoBoundaries"
谢谢

1334335264

Soaring 发表于 2017-3-20 16:29
好像不对，我想画北京市具体的图，你这个好像只是为了画中国地图而描绘的每个省的大致轮廓，那每个省具体 ...

我每个省的时候shp文件。。

蜻蜓队长

Soaring 发表于 2017-3-20 21:41
你好，还有一个小问题：
begin

请问你这个.shp文件哪里找到的呢