ncl如何添加线shp文件_NCL画图个例讲解

要在NCL中添加线形状文件（.shp），您需要使用以下代码：

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"

begin
;-- Open NCAR Graphics
  wks = gsn_open_wks("png","plot")

;-- Read in the shapefile
  shp_file = addfile("file.shp", "r")

;-- Set the map projection
  mpProjection = "lcc"
  mpDataBaseVersion = "MediumRes"
  mpLimitMode = "LatLon"
  mpMinLatF = 24.
  mpMaxLatF = 50.
  mpMinLonF = -125.
  mpMaxLonF = -66.

  mpOutlineOn = False

  mpFillOn = True
  mpFillDrawOrder = "PostDraw"

;-- Set the color map
  cnFillPalette = "BlWhRe"

;-- Create the plot
  plot = gsn_csm_vector(wks, shp_file, "arc", "Line", False, False, mpProjection,\
                        mpDataBaseVersion, mpLimitMode, mpMinLatF, mpMaxLatF,\
                        mpMinLonF, mpMaxLonF, mpOutlineOn, mpFillOn, mpFillDrawOrder,\
                        cnFillPalette)

;-- Draw the plot
  draw(plot)

;-- Close NCAR Graphics
  frame(wks)
  end

其中，您需要更改以下内容：

- `wks`：指定绘图设备，这里我使用 `png` 作为输出格式，可以根据需要更改。
- `shp_file`：指定您要添加的 .shp 文件的路径和文件名。
- `mpProjection`：指定地图投影方式，这里我使用的是 lambert conformal，可以根据需要更改。
- `mpMinLatF`、`mpMaxLatF`、`mpMinLonF`、`mpMaxLonF`：指定地图的经纬度范围，根据您的数据范围进行更改。
- `cnFillPalette`：指定颜色表，这里我使用黑白红三色渐变，可以根据需要更改。

运行这段代码，您将得到一个包含 .shp 文件中所有线形状的地图。