画青藏高原轮廓错误

当时年少春衫薄

想画轮廓，不知道为什么会是阴影

denghaijun

这是因为原始的shp文件导致的，那不是阴影，是边线太多造成的，
重新下载这个边界，应该就没问题了，http://www.geodoi.ac.cn/WebCn/doi.aspx?Id=135

当时年少春衫薄

请问各位这是shp文件的问题还是代码有错误

当时年少春衫薄

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"
begin
  wks = gsn_open_wks("png","coast")   ;   Open X11 window
dir = "/home/fanweiwei/NCEP/heat/TP/shp/Qinghai-Tibet_Plateau/"
filename =dir + "Qinghai-Tibet_Plateau.shp"
  res               = True
  res@gsnMaximize   = True
  res@gsnDraw       = False
  res@gsnFrame      = False

  res@mpOutlineOn   = False
  res@mpFillOn      = False  

;---Turn on fancier tickmark labels.


;---Zoom in on Australia
  res@mpLimitMode           = "LatLon"
  res@mpMinLatF             = 20
  res@mpMaxLatF             =  60
  res@mpMinLonF             = 60
  res@mpMaxLonF             = 110

;---Create map
  map = gsn_csm_map(wks,res)

;---Attach three different sets of shapefile outlines to map
  lnres                  = True
  lnres@gsLineColor      = "brown"
  lnres@gsLineThicknessF = 1.0

  p = gsn_add_shapefile_polylines(wks,map,filename,lnres)

  lnres@gsLineColor      = "purple"
  lnres@gsLineThicknessF = 2.0


  draw(map)
  frame(wks)
end

御风者

我觉得可以是你的shapefile本身存在问题，你的shapefile来源是什么？

当时年少春衫薄

御风者 发表于 2016-6-5 15:37
我觉得可以是你的shapefile本身存在问题，你的shapefile来源是什么？

你好我这个是在论坛里下载的，请问还有别的下载途径吗

当时年少春衫薄

denghaijun 发表于 2016-6-5 23:16
这是因为原始的shp文件导致的，那不是阴影，是边线太多造成的，
重新下载这个边界，应该就没问题了，http: ...

哦哦好的非常感谢

AnneChen

我也想要青藏高原边界图呢。感谢分享讨论

denghaijun

当时年少春衫薄 发表于 2016-6-6 06:57
哦哦好的非常感谢

客气，多多交流。

自行车大大爷

denghaijun 发表于 2016-6-5 23:16
这是因为原始的shp文件导致的，那不是阴影，是边线太多造成的，
重新下载这个边界，应该就没问题了，http: ...

我也遇到了这个问题，非常感谢您的回答！