91超碰碰碰碰久久久久久综合_超碰av人澡人澡人澡人澡人掠_国产黄大片在线观看画质优化_txt小说免费全本

溫馨提示×

溫馨提示×

您好,登錄后才能下訂單哦!

密碼登錄×
登錄注冊×
其他方式登錄
點擊 登錄注冊 即表示同意《億速云用戶服務條款》

NCL常用統計學函數怎么用

發布時間:2022-01-07 09:17:50 來源:億速云 閱讀:1565 作者:iii 欄目:大數據

這篇文章主要介紹“NCL常用統計學函數怎么用”,在日常操作中,相信很多人在NCL常用統計學函數怎么用問題上存在疑惑,小編查閱了各式資料,整理出簡單好用的操作方法,希望對大家解答”NCL常用統計學函數怎么用”的疑惑有所幫助!接下來,請跟著小編一起來學習吧!

NCL作為一門氣象專業語言,自帶了很多氣象屆常用的算法和命令, 比如計算相關系數、均方根誤差等統計量,相關系數顯著性檢驗,EOF,以及各種插值。這里將一些常用函數做了歸納:

1、相關系數
NCL中計算相關系數的函數有好幾個,常用的為escorc系列和regCoef系列。

escorc(x,y)是計算兩個數組最右側一列的相關系數。如果想指定計算的列數,比如想計算x(lat,lon)和y(lat,lon,year)在lat這一列上的相關系數,就要用到escorc_n: escorc_n(x,y,0)。

相關系數的檢驗方法包括t檢驗、f檢驗和顯著性檢驗,分別為ttest,ftest,rtest。

regCoef系列的命令比escorc強大很多,因為它除了返回相關系數rc,還捎帶手把兩個序列的標準差(rc@rstd)和t檢驗結果(rc@tval)都給算了。

2、均方根誤差
函數名為dim_rmsd。它是計算兩個變量在所有其他維度上最右側維度的均方根誤差。而對于想要計算指定維度的情況,則要用dim_rmsd_n。

3、EOF分解
計算EOF分解的函數為eofunc開頭的系列函數,一般用eofunc_Wrap計算空間模態,用eofunc_ts_Wrap計算時間系數。

以下為針對變量x計算EOF分解并繪圖的NCL子程序:

undef("draw_eof_plot")procedure draw_eof_plot(dir_plot,file_plot,type_plot,x,year,neof)begin;--EOF analysisoptEof = Trueeof    = eofunc_Wrap( x, neof, optEof)eof_ts = eofunc_ts_Wrap( x, eof, False)lat = x&latlon = x&lon;--Begin plotting section.wks  = gsn_open_wks(type_plot,dir_plot+file_plot)      ; Opens a ps file gsn_define_colormap(wks,"rainbow")   plot = new(neof,graphic)   res              = True                       ; plot mods desired;************************************************; original data;************************************************   res@gsnDraw                    = False        ; don't draw yet   res@gsnFrame                   = False        ; don't advance frame yet   res@gsnAddCyclic             = False;--map plot resources  res@mpFillOn                   = False        ; no grey continents  res@mpCenterLonF      = 180.  res@mpDataBaseVersion     = "MediumRes"    ; or "Ncarg4_1"  res@mpDataSetName="Earth..4"  res@mpOutlineSpecifiers=(/"China:states","Taiwan"/)  res@mpOutlineBoundarySets = "AllBoundaries"  res@mpMinLatF  = min(lat)            ; range to zoom in on  res@mpMaxLatF  = max(lat)  res@mpMinLonF  = min(lon)  res@mpMaxLonF  = max(lon)
;--contour resources  res@cnFillOn             = True              ; turn on contour fill  res@cnLineLabelsOn       = False              ; turn off contour  res@cnLinesOn            = False        ; add countor or not,True is default  res@gsnLeftString = " "  res@gsnRightString = " "
;--tickmark resources  res@tmXTOn            = False  res@tmYROn            = False  res@tmYLLabelFontHeightF =0.02  res@tmXBLabelFontHeightF =0.018  res@tmXTOn            = False  res@tmYROn            = False
;--labelbar resources  res@lbLabelBarOn         = False  res@lbLabelFontHeightF   = 0.02  res@lbOrientation       = "Vertical"          ; vertical label bar
 symMinMaxPlt(eof, 16, False, res); contributed.ncl; panel plot only resources  resP                     = True         ; modify the panel plot  resP@gsnMaximize         = True         ; large format  resP@gsnPanelLabelBar    = True         ; add common colorbar  resP@txString  = "EOF"do n=0,neof-1        res@gsnLeftString  = " EOF "+(n+1)        res@gsnRightString = sprintf("%5.1f", eof@pcvar(n)) +"%"        plot(n) = gsn_csm_contour_map(wks,eof(n,:,:),res)end do   gsn_panel(wks,plot,(/neof/1,2/),resP)     ; draw all 'neof' as one plot;*******************************************; time series (principal component) plot;*******************************************  eof_ts@long_name = "Amplitude"  rts           = True  rts@gsnDraw   = False       ; don't draw yet  rts@gsnFrame  = False       ; don't advance frame yet
; decide exactly where on the page to draw it.  rts@vpHeightF = 0.40        ; Changes the aspect ratio  rts@vpWidthF  = 0.85  rts@vpXF      = 0.10        ; change start locations  rts@vpYF      = 0.75        ; the plot  rts@gsnYRefLine           = 0.              ; reference line  rts@gsnAboveYRefLineColor = "red"           ; above ref line fill red  rts@gsnBelowYRefLineColor = "blue"          ; below ref line fill blue; panel plot only resources  rtsP                     = True             ; modify the panel plot  rtsP@gsnMaximize         = True             ; large format  rtsP@txString  = "EOF"  do n=0,neof-1        rts@gsnLeftString  = " EOF "+(n+1)        rts@gsnRightString = sprintf("%5.1f", eof@pcvar(n)) +"%"        plot(n) = gsn_csm_xy (wks,year,eof_ts(n,:),rts)  end do  gsn_panel(wks,plot,(/neof/1,2/),rtsP)        ; draw all 'neof' as one plotend

到此,關于“NCL常用統計學函數怎么用”的學習就結束了,希望能夠解決大家的疑惑。理論與實踐的搭配能更好的幫助大家學習,快去試試吧!若想繼續學習更多相關知識,請繼續關注億速云網站,小編會繼續努力為大家帶來更多實用的文章!

向AI問一下細節

免責聲明:本站發布的內容(圖片、視頻和文字)以原創、轉載和分享為主,文章觀點不代表本網站立場,如果涉及侵權請聯系站長郵箱:is@yisu.com進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。

ncl
AI

修水县| 兴国县| 修文县| 乐安县| 合江县| 永仁县| 阜新| 金昌市| 合肥市| 临江市| 乐安县| 全南县| 梨树县| 岫岩| 大田县| 墨江| 浦江县| 射洪县| 多伦县| 青州市| 仁怀市| 江川县| 泾阳县| 民勤县| 花莲市| 固阳县| 平度市| 寻乌县| 灌南县| 迁西县| 清水河县| 调兵山市| 乌海市| 金秀| 遂溪县| 光山县| 济源市| 罗源县| 博乐市| 竹山县| 沧源|