中国shp格式底图，河流，省界，市界

黄小仙儿

从地理信息网上下的，需要的人可以来下
添加底图的例子可以参考NCL的shapefiles板块
这是链接http://www.ncl.ucar.edu/Applications/shapefiles.shtml这个是官网上wrf数据的例子，大家可以先参考一下
;----------------------------------------------------------------------
; This is a basic NCL template for creating contours over maps using
; WRF data, and attaching shapefile outlines. This template plots one
; timestep of HGT from a WRF file, and attaches shapefile outlines
; from the "USA_adm2.shp" and "CAN_adm2.shp" shapefiles, downloaded
; from http://www.gadm.org/country/
;----------------------------------------------------------------------
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/wrf/WRFUserARW.ncl"

begin
  filename = "wrfout_d01_2002-07.nc"
  a        = addfile(filename,"r")

  ter = a->HGT(0,:,:)                        ; Read the variable to memory

  wks = gsn_open_wks("x11","wrf_hgt_shapefile")
  gsn_define_colormap(wks,"OceanLakeLandSnow")  ; Change color map

  res                   = True               ; Use basic options for this field
  res@cnFillOn          = True               ; Create a color fill plot
  res@ContourParameters = (/1,1100,20/)

  contour = wrf_contour(a,wks,ter,res)

  pltres            = True   ; Set plot options
  pltres@PanelPlot  = True   ; Tells wrf_map_overlays not to remove contours

  mpres             = True   ; Set map options
  mpres@mpOutlineOn = False
  mpres@mpFillOn    = False

;---Create the contours over the WRF map (nothing will be drawn yet).
  plot = wrf_map_overlays(a,wks,contour,pltres,mpres)

;---Attach the shapefile polylines using files read off gadm.org/country.
  usa_shp_name    = "USA_adm/USA_adm2.shp"
  canada_shp_name = "CAN_adm/CAN_adm2.shp"

  lnres                  = True
  lnres@gsLineColor      = "gray25"
  lnres@gsLineThicknessF = 0.5   
  usa_id = gsn_add_shapefile_polylines(wks,plot,usa_shp_name,lnres)
  can_id = gsn_add_shapefile_polylines(wks,plot,canada_shp_name,lnres)

  draw(plot)       ; This will draw the map and the shapefile outlines.
  frame(wks)       ; Advance the frame

end
还有这个是ncl对grib,hdf,wrf数据处理及格式的讲解pdf,需要的也可以下载看看

syl136

为什么每次都提示.SHP文件打不开。。。

andrewsoong

顶~~~PS：要是把你自己的NCL脚本例子和画的图放上就更好了~~~

黄小仙儿

andrewsoong 发表于 2014-4-13 16:06
顶~~~PS：要是把你自己的NCL脚本例子和画的图放上就更好了~~~

恩，我现在添加了地区界的底图，但是省界那里是空的。不知道怎么在一张图上即设置省界，也设置市界。
现在还在请教其他人，等我把这个弄好了，就贴上来

黄小仙儿

这个是官网上wrf数据的例子，大家可以先参考一下
;----------------------------------------------------------------------
; This is a basic NCL template for creating contours over maps using
; WRF data, and attaching shapefile outlines. This template plots one
; timestep of HGT from a WRF file, and attaches shapefile outlines
; from the "USA_adm2.shp" and "CAN_adm2.shp" shapefiles, downloaded
; from http://www.gadm.org/country/
;----------------------------------------------------------------------
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/wrf/WRFUserARW.ncl"

begin
  filename = "wrfout_d01_2002-07.nc"
  a        = addfile(filename,"r")

  ter = a->HGT(0,:,:)                        ; Read the variable to memory

  wks = gsn_open_wks("x11","wrf_hgt_shapefile")
  gsn_define_colormap(wks,"OceanLakeLandSnow")  ; Change color map

  res                   = True               ; Use basic options for this field
  res@cnFillOn          = True               ; Create a color fill plot
  res@ContourParameters = (/1,1100,20/)

  contour = wrf_contour(a,wks,ter,res)

  pltres            = True   ; Set plot options
  pltres@PanelPlot  = True   ; Tells wrf_map_overlays not to remove contours

  mpres             = True   ; Set map options
  mpres@mpOutlineOn = False
  mpres@mpFillOn    = False

;---Create the contours over the WRF map (nothing will be drawn yet).
  plot = wrf_map_overlays(a,wks,contour,pltres,mpres)

;---Attach the shapefile polylines using files read off gadm.org/country.
  usa_shp_name    = "USA_adm/USA_adm2.shp"
  canada_shp_name = "CAN_adm/CAN_adm2.shp"

  lnres                  = True
  lnres@gsLineColor      = "gray25"
  lnres@gsLineThicknessF = 0.5   
  usa_id = gsn_add_shapefile_polylines(wks,plot,usa_shp_name,lnres)
  can_id = gsn_add_shapefile_polylines(wks,plot,canada_shp_name,lnres)

  draw(plot)       ; This will draw the map and the shapefile outlines.
  frame(wks)       ; Advance the frame

end

andrewsoong

黄小仙儿 发表于 2014-4-13 16:09
恩，我现在添加了地区界的底图，但是省界那里是空的。不知道怎么在一张图上即设置省界，也设置市界。
现 ...

kongfeng0824

谢谢分享   

rjys1219

非常好，学习中

陌上尘

这个怎么用啊，摸索一下！！！！

nuist2006

呵呵，学习学习

yrovl

感谢楼主！