求问一个散点图加上对角线的问题

stepdance

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"

procedure add_line(wks,plot,x,y)
local res_line, str
begin
  res_line = True
  res_line@xyDashPatterns         = 16                  ; dashed line for 2nd
  res_line@xyLineThicknesses      = 3                  ; thicker line

  str = unique_string("polyline")   ; "unique_string" will return a unique
                                    ; string every time it is called from
                                    ;  within a single NCL session.
;
; You can then use this unique string as an attribute variable name
; that gets attached to the plot variable. This ensures that this
; value will live for the duration of the script.
;
  plot@$str$ = gsn_add_polyline(wks, plot, x, y, res_line)
end
;---------------仿照例子定义一个add_line的procedure，用来画对角线，输入的参数是wks, plot, x和y
begin

Rn_path="/home/zhaojc/work/863/Rn/"
f_Rn_cam=addfile(Rn_path+"Rn_cam_dot.nc","r")
var_Rn_cam=f_Rn_cam->Rn_cam_dot
f_Rn_rs=addfile(Rn_path+"Rn_rs.nc","r")
var_Rn_rs=f_Rn_rs->Rn_rs

x_Rn=var_Rn_cam
y_Rn=var_Rn_rs
;--------------------------------------将要化的数据读取出来---------------------------------------------
xx_Rn = ispan (60,72,1)
yy_Rn = xx_Rn
；--------------------------------------生成准备画对角线的数据------------------------------------------------------
wks   = gsn_open_wks ("pdf","scatter")            ; open workstation
  res_Rn = True
  res_Rn@xyMarkLineMode    = "Markers"                ; choose to use markers
  res_Rn@xyMarkers         =  16                      ; choose type of marker
  res_Rn@xyMarkerColor     = "Black"               ; Marker color
  res_Rn@xyMarkerSizeF     = 0.01                     ; Marker size (default 0.01)
  res_Rn@trXMinF = 60
  res_Rn@trYMinF = 60
  res_Rn@trXMaxF = 72
  res_Rn@trYMaxF = 72
  res_Rn@gsnFrame = False                  ; don't advance frame yet
  res_Rn@tiXAxisString = "ERA-40"
  res_Rn@tiXAxisFont = 25
  res_Rn@tiYAxisString = "calculation"
  res_Rn@tiYAxisFont = 25

   plot  = gsn_csm_xy (wks,x_Rn,y_Rn,res_Rn)
；---------------------------------------------------画散点图------------------------------------------------------------------
  add_line(wks,plot,xx_Rn,yy_Rn)
；---------------------------------------------------调用add_line画对角线------------------------------------------------
  frame(wks)
end

但这样画出来，只出来散点图，没有对角线，求问错在什么地方？
多谢！

talkd

hi，完全可以参照http://www.ncl.ucar.edu/Applications/scatter.shtml 中例4将趋势线换为你要的对角线就行了

stepdance

talkd 发表于 2013-11-21 12:04

                               
登录/注册后可看大图

hi，完全可以参照http://www.ncl.ucar.edu/Applications/scatter.shtml 中例4将趋势线换为你要的对角线就行 ...

非常感谢，我试试去

abc258886

subtropical

talkd 发表于 2013-11-21 12:04
hi，完全可以参照http://www.ncl.ucar.edu/Applications/scatter.shtml 中例4将趋势线换为你要的对角线就行 ...

谢谢您，刚好要用到

subtropical

talkd 发表于 2013-11-21 12:04
hi，完全可以参照http://www.ncl.ucar.edu/Applications/scatter.shtml 中例4将趋势线换为你要的对角线就行 ...

您好，有个问题想请教一下，我画了一个因子先简称A和台风数的散点图,结果表达得是A因子越强台风数越少，
然后用ncl里面的第四个例子http://www.ncl.ucar.edu/Applications/Scripts/scatter_4.ncl添加了一条直线regline，请问我写图注的时候怎么写，就说A和台风数的散点图及其回归线吗？我在文献中没找到这种图，所以图注不太会写。。。求大神指教

talkd

subtropical 发表于 2016-10-19 17:34
您好，有个问题想请教一下，我画了一个因子先简称A和台风数的散点图,结果表达得是A因子越强台风数越少，
...

应该有很多这种图吧，Caption可以大概这样写：Scatter plot of A against the Number of Typhoon. The dash line indicates the regression line. 仅供参考...

subtropical

talkd 发表于 2016-10-21 10:17
应该有很多这种图吧，Caption可以大概这样写：Scatter plot of A against the Number of Typhoon. The da ...

好的灰常感谢{:eb511:}

橙子鲜

talkd 发表于 2013-11-21 12:04
hi，完全可以参照http://www.ncl.ucar.edu/Applications/scatter.shtml 中例4将趋势线换为你要的对角线就行 ...

你好，我已经使用了例子中介绍的方法，画出散点跟对角线图，现在想再画出x轴y轴一横一竖的图，该怎么添加呢？

talkd

橙子鲜 发表于 2016-11-18 15:45
你好，我已经使用了例子中介绍的方法，画出散点跟对角线图，现在想再画出x轴y轴一横一竖的图，该怎么添加 ...

设置XRefline和YRefline就行了