- 积分
- 810
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-5-12
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2015-1-12 22:29:33
|
显示全部楼层
- ;----------------------------------------------------------------------
- ; station_2.ncl
- ;----------------------------------------------------------------------
- ;
- ; Concepts illustrated:
- ; - Drawing markers on a map indicating the locations of station data
- ; - Generating dummy data using "random_uniform"
- ; - Drawing markers of different sizes and colors on a map
- ; - Drawing a custom legend outside of a map plot
- ; - Attaching a custom labelbar to a plot
- ;
- ; This example creates some dummy station data, and then plots each
- ; value by coloring and sizing it depending on which range of values
- ; it falls in.
- ;
- ; It creates two plots: one with a legend with markers and text,
- ; and one with a labelbar.
- ;----------------------------------------------------------------------
- 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"
- ;----------------------------------------------------------------------
- ; This procedure draws a legend with markers and text at the bottom
- ; of the screen
- ;----------------------------------------------------------------------
- procedure draw_legend(wks,lat[*][*]:numeric,lon[*][*]:numeric,\
- arr[*]:numeric,colors)
- local gsres, txres, xleg, xtxt, yleg, ytxt, i, labels, nitems
- begin
- narr = dimsizes(arr)
- nmarkers = narr+1
- labels = new(nmarkers,string)
- ;---Generate the labels for each marker.
- do i = 0, nmarkers-1
- if (i.eq.0) then
- labels(i) = "x < " + arr(0)
- end if
- if (i.eq.nmarkers-1) then
- labels(i) = "x >= " + max(arr)
- end if
- if (i.gt.0.and.i.lt.nmarkers-1) then
- labels(i) = arr(i-1) + " <= x < " + arr(i)
- end if
- end do
- ;
- ; Create logical variables to hold the marker and text resources.
- ; These markers are different than the XY markers, because they are not
- ; associated with an XY plot. You can put these markers on any plot.
- ;
- gsres = True
- gsres@gsMarkerIndex = 16 ; Use filled dots for markers.
- txres = True
- txres@txFontHeightF = 0.015
- ;
- ; Loop through each grouping of markers, and draw them one set at
- ; a time, assigning the proper color and size with gsn_marker.
- ;
- ; At the same time, draw a legend showing the meaning of the markers.
- ;
- xleg = (/0.07,0.07,0.32,0.32,0.57,0.57,0.82,0.82/) ; Location of
- xtxt = (/0.16,0.16,0.41,0.41,0.66,0.66,0.91,0.91/) ; legend markers
- yleg = (/0.22,0.17,0.22,0.17,0.22,0.17,0.22,0.17/) ; and text
- ytxt = (/0.22,0.17,0.22,0.17,0.22,0.17,0.22,0.17/) ; strings.
- do i = 0, dimsizes(lat(:,0))-1
- if (.not.ismissing(lat(i,0)))
- gsres@gsMarkerColor = colors(i)
- gsres@gsMarkerThicknessF = 0.7*(i+1)
- ;---Add marker and text for the legend.
- gsn_polymarker_ndc(wks, xleg(i),yleg(i),gsres)
- gsn_text_ndc (wks,labels(i),xtxt(i),ytxt(i),txres)
- end if
- end do
- end
- ;----------------------------------------------------------------------
- ; This function attaches a labelbar to the given plot.
- ;----------------------------------------------------------------------
- function attach_labelbar(wks,map,arr[*]:numeric,colors[*])
- local lbres, vph, vpw, nboxes
- begin
- getvalues map
- "vpHeightF" : vph
- "vpWidthF" : vpw
- end getvalues
- nboxes = dimsizes(colors)
-
-
- lbres = True ; labelbar only resources
- lbres@lbAutoManage = False ; Necessary to control sizes
- lbres@lbFillColors = colors
- lbres@vpWidthF = 0.7 * vpw ; labelbar width
- lbres@vpHeightF = 0.1 * vph ; labelbar height
- lbres@lbMonoFillPattern = True ; Solid fill pattern
- lbres@lbLabelFontHeightF = 0.01 ; font height. default is small
- lbres@lbOrientation = "horizontal"
- ;lbres@lbOrientation = "vertical"
- ;lbres@lbLabelPosition ="Right"
- lbres@lbPerimOn = False
- lbres@lbLabelAlignment = "InteriorEdges"
- lbid = gsn_create_labelbar(wks,nboxes,""+toint(arr),lbres)
- amres = True
- amres@amJust = "TopCenter"
- amres@amParallelPosF = 0.0 ; Center
- amres@amOrthogonalPosF = 0.6 ; Bottom
- annoid = gsn_add_annotation(map,lbid,amres)
- return(annoid)
- end
- ;----------------------------------------------------------------------
- ; Main code
- ;----------------------------------------------------------------------
- filename = new(3,string)
- filename(0) = "flood_T"
- filename(1) = "flood_S"
- filename(2) = "flood_TS"
- begin
- ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;Tide;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
- ;---Set some needed arrays
- arr = ispan(0,450,15)
- narr = dimsizes(arr)
- labels = new(narr+1,string) ; Labels for legend.
- nrow_T = numAsciiRow("flood_T.dat")
- ncol_T = numAsciiCol("flood_T.dat")
- print("flood_T.dat")
- data_T = asciiread("flood_T.dat",(/nrow_T,ncol_T/),"float")
- lat_T = data_T(:,1)
- lon_T = data_T(:,0)
- lat_T@units= "degrees-north"
- lon_T@units= "degrees-east"
- R_T = data_T(:,2)*100
- npts_T = nrow_T ; Number of points.
- ;
- ; Create X and Y arrays to hold the points for each range and initialize
- ; them to missing values. We want to use num_distinct_markers
- ; different colors, so we need num_distinct_markers sets of X and
- ; Y points.
- ;
- num_distinct_markers = dimsizes(arr)+1 ; number of distinct markers
- lat_new_T = new((/num_distinct_markers,dimsizes(R_T)/),float,-999)
- lon_new_T = new((/num_distinct_markers,dimsizes(R_T)/),float,-999)
- ;---Group the points according to which range they fall in.
- do i = 0, num_distinct_markers-1
- if (i.eq.0) then
- indexes = ind(R_T.lt.arr(0))
- end if
- if (i.eq.num_distinct_markers-1) then
- indexes = ind(R_T.ge.max(arr))
- end if
- if (i.gt.0.and.i.lt.num_distinct_markers-1) then
- indexes = ind(R_T.ge.arr(i-1).and.R_T.lt.arr(i))
- end if
- ;
- ; Now that we have the set of indexes whose values fall within
- ; the given range, take the corresponding lat/lon values and store
- ; them, so later we can color this set of markers with the appropriate
- ; color.
- ;
- if (.not.any(ismissing(indexes))) then
- npts_range = dimsizes(indexes) ; # of points in this range.
- lat_new_T(i,0:npts_range-1) = lat_T(indexes)
- lon_new_T(i,0:npts_range-1) = lon_T(indexes)
- end if
- delete(indexes) ; Necessary b/c "indexes" may be a different
- ; size next time.
- end do
- ;----------------------------------------------------------------------
- ; Begin plotting section.
- ;----------------------------------------------------------------------
- wks = gsn_open_wks("ps","flood_T") ; Open a workstation and
- ;---Set up some map resources.
- mpres = True
-
-
-
- mpres@gsnMaximize = True ; Maximize plot in frame.
- mpres@gsnFrame = False ; Don't advance the frame
- mpres@gsnDraw = False ; Don't advance the frame
-
- mpres@tmXTBorderOn = True ; Don't draw top axis.
- mpres@tmXTOn = True
- mpres@tmXBMode = "Explicit" ; Set tick mark mode.
- mpres@tmXBValues = ispan(10940,11090,30)/100.0
- mpres@tmXBLabels = ispan(10940,11090,30)/100.0
- mpres@tmYLMode = "Explicit" ; Set tick mark mode.
- mpres@tmYLValues = ispan(202,215,2)/10.0
- mpres@tmYLLabels = ispan(202,215,2)/10.0
- ;---Zoom in on United States.
- mpres@mpMinLatF = 20.2
- mpres@mpMaxLatF = 21.5
- mpres@mpMinLonF = 109.40
- mpres@mpMaxLonF = 110.9
- mpres@mpCenterLonF = 180
- mpres@mpDataBaseVersion = "HighRes" ; High resolution database
- mpres@mpDataResolution = "Finest"
- mpres@mpDataSetName = "Earth..3" ;
- mpres@mpFillColors = (/"transparent","transparent","gray","transparent"/)
- mpres@tiMainString = "FloodPlain Scatters "
- map = gsn_csm_map(wks,mpres)
-
- gsres = True
- gsres@gsMarkerIndex = 1 ; Use filled dots for markers.
- ;---Get nice spacing through color map for marker colors
- gsn_define_colormap(wks,"Rainbow")
- getvalues wks
- "wkColorMapLen" : clen ; number of colors in color map
- end getvalues
- nstep = (clen-2)/narr
- colors = ispan(2,clen-1,nstep)
- ;---Loop through each "bin" and attach the markers to the map.
- do i = 0, num_distinct_markers-1
- if (.not.ismissing(lat_new_T(i,0)))
- gsres@gsMarkerColor = colors(i)
- dumstr = unique_string("marker")
- map@$dumstr$ = gsn_add_polymarker(wks,map,lon_new_T(i,:),lat_new_T(i,:),gsres)
- end if
- end do
- ;----------------------------------------------------------------------
- ; First version of this plot has a legend at the bottom
- ;----------------------------------------------------------------------
- ;draw(map) ; Drawing the map will draw the markers
- ;draw_legend(wks,lat_new,lon_new,arr,colors) ; Draw a legend at the bottom
- ;frame(wks) ; Advance the frame
- ;----------------------------------------------------------------------
- ; Second version of this plot has a labelbar added.
- ;----------------------------------------------------------------------
- lbid = attach_labelbar(wks,map,arr,colors) ; Attach a labelbar
- draw(map) ; Drawing the map will draw everything
-
- frame(wks) ; Advance the frame.
- ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;Surge;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
- ;---Set some needed arrays
- arr = ispan(0,450,15)
- narr = dimsizes(arr)
- labels = new(narr+1,string) ; Labels for legend.
- nrow_S = numAsciiRow("flood_S.dat")
- ncol_S = numAsciiCol("flood_S.dat")
- print("flood_S.dat")
- data_S = asciiread("flood_S.dat",(/nrow_S,ncol_S/),"float")
- lat_S = data_S(:,1)
- lon_S = data_S(:,0)
- lat_S@units= "degrees-north"
- lon_S@units= "degrees-east"
- R_S = data_S(:,2)*100
- npts_S = nrow_S ; Number of points.
- ;
- ; Create X and Y arrays to hold the points for each range and initialize
- ; them to missing values. We want to use num_distinct_markers
- ; different colors, so we need num_distinct_markers sets of X and
- ; Y points.
- ;
- num_distinct_markers = dimsizes(arr)+1 ; number of distinct markers
- lat_new_S = new((/num_distinct_markers,dimsizes(R_S)/),float,-999)
- lon_new_S = new((/num_distinct_markers,dimsizes(R_S)/),float,-999)
- ;---Group the points according to which range they fall in.
- do i = 0, num_distinct_markers-1
- if (i.eq.0) then
- indexes = ind(R_S.lt.arr(0))
- end if
- if (i.eq.num_distinct_markers-1) then
- indexes = ind(R_S.ge.max(arr))
- end if
- if (i.gt.0.and.i.lt.num_distinct_markers-1) then
- indexes = ind(R_S.ge.arr(i-1).and.R_S.lt.arr(i))
- end if
- ;
- ; Now that we have the set of indexes whose values fall within
- ; the given range, take the corresponding lat/lon values and store
- ; them, so later we can color this set of markers with the appropriate
- ; color.
- ;
- if (.not.any(ismissing(indexes))) then
- npts_range = dimsizes(indexes) ; # of points in this range.
- lat_new_S(i,0:npts_range-1) = lat_S(indexes)
- lon_new_S(i,0:npts_range-1) = lon_S(indexes)
- end if
- delete(indexes) ; Necessary b/c "indexes" may be a different
- ; size next time.
- end do
- ;----------------------------------------------------------------------
- ; Begin plotting section.
- ;----------------------------------------------------------------------
- wks = gsn_open_wks("ps","flood_S") ; Open a workstation and
- ;---Set up some map resources.
- mpres = True
-
-
-
- mpres@gsnMaximize = True ; Maximize plot in frame.
- mpres@gsnFrame = False ; Don't advance the frame
- mpres@gsnDraw = False ; Don't advance the frame
-
- mpres@tmXTBorderOn = True ; Don't draw top axis.
- mpres@tmXTOn = True
- mpres@tmXBMode = "Explicit" ; Set tick mark mode.
- mpres@tmXBValues = ispan(10940,11090,30)/100.0
- mpres@tmXBLabels = ispan(10940,11090,30)/100.0
- mpres@tmYLMode = "Explicit" ; Set tick mark mode.
- mpres@tmYLValues = ispan(202,215,2)/10.0
- mpres@tmYLLabels = ispan(202,215,2)/10.0
- ;---Zoom in on United States.
- mpres@mpMinLatF = 20.2
- mpres@mpMaxLatF = 21.5
- mpres@mpMinLonF = 109.40
- mpres@mpMaxLonF = 110.9
- mpres@mpCenterLonF = 180
- mpres@mpDataBaseVersion = "HighRes" ; High resolution database
- mpres@mpDataResolution = "Finest"
- mpres@mpDataSetName = "Earth..3" ;
- mpres@mpFillColors = (/"transparent","transparent","gray","transparent"/)
- mpres@tiMainString = "FloodPlain Scatters "
- map = gsn_csm_map(wks,mpres)
-
- gsres = True
- gsres@gsMarkerIndex = 1 ; Use filled dots for markers.
- ;---Get nice spacing through color map for marker colors
- gsn_define_colormap(wks,"Rainbow")
- getvalues wks
- "wkColorMapLen" : clen ; number of colors in color map
- end getvalues
- nstep = (clen-2)/narr
- colors = ispan(2,clen-1,nstep)
- ;---Loop through each "bin" and attach the markers to the map.
- do i = 0, num_distinct_markers-1
- if (.not.ismissing(lat_new_S(i,0)))
- gsres@gsMarkerColor = colors(i)
- dumstr = unique_string("marker")
- map@$dumstr$ = gsn_add_polymarker(wks,map,lon_new_S(i,:),lat_new_S(i,:),gsres)
- end if
- end do
- ;----------------------------------------------------------------------
- ; First version of this plot has a legend at the bottom
- ;----------------------------------------------------------------------
- ;draw(map) ; Drawing the map will draw the markers
- ;draw_legend(wks,lat_new,lon_new,arr,colors) ; Draw a legend at the bottom
- ;frame(wks) ; Advance the frame
- ;----------------------------------------------------------------------
- ; Second version of this plot has a labelbar added.
- ;----------------------------------------------------------------------
- lbid = attach_labelbar(wks,map,arr,colors) ; Attach a labelbar
- draw(map) ; Drawing the map will draw everything
-
- frame(wks) ; Advance the frame.
- ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;Tide_Surge;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
- ;---Set some needed arrays
- arr = ispan(0,450,15)
- narr = dimsizes(arr)
- labels = new(narr+1,string) ; Labels for legend.
- nrow_TS = numAsciiRow("flood_TS.dat")
- ncol_TS = numAsciiCol("flood_TS.dat")
- print("flood_TS.dat")
- data_TS = asciiread("flood_TS.dat",(/nrow_TS,ncol_TS/),"float")
- lat_TS = data_TS(:,1)
- lon_TS = data_TS(:,0)
- lat_TS@units= "degrees-north"
- lon_TS@units= "degrees-east"
- R_TS = data_TS(:,2)*100
- npts_TS = nrow_TS ; Number of points.
- ;
- ; Create X and Y arrays to hold the points for each range and initialize
- ; them to missing values. We want to use num_distinct_markers
- ; different colors, so we need num_distinct_markers sets of X and
- ; Y points.
- ;
- num_distinct_markers = dimsizes(arr)+1 ; number of distinct markers
- lat_new_TS = new((/num_distinct_markers,dimsizes(R_TS)/),float,-999)
- lon_new_TS = new((/num_distinct_markers,dimsizes(R_TS)/),float,-999)
- ;---Group the points according to which range they fall in.
- do i = 0, num_distinct_markers-1
- if (i.eq.0) then
- indexes = ind(R_TS.lt.arr(0))
- end if
- if (i.eq.num_distinct_markers-1) then
- indexes = ind(R_TS.ge.max(arr))
- end if
- if (i.gt.0.and.i.lt.num_distinct_markers-1) then
- indexes = ind(R_TS.ge.arr(i-1).and.R_TS.lt.arr(i))
- end if
- ;
- ; Now that we have the set of indexes whose values fall within
- ; the given range, take the corresponding lat/lon values and store
- ; them, so later we can color this set of markers with the appropriate
- ; color.
- ;
- if (.not.any(ismissing(indexes))) then
- npts_range = dimsizes(indexes) ; # of points in this range.
- lat_new_TS(i,0:npts_range-1) = lat_TS(indexes)
- lon_new_TS(i,0:npts_range-1) = lon_TS(indexes)
- end if
- delete(indexes) ; Necessary b/c "indexes" may be a different
- ; size next time.
- end do
- ;----------------------------------------------------------------------
- ; Begin plotting section.
- ;----------------------------------------------------------------------
- wks = gsn_open_wks("ps","flood_TS") ; Open a workstation and
- ;---Set up some map resources.
- mpres = True
-
-
-
- mpres@gsnMaximize = True ; Maximize plot in frame.
- mpres@gsnFrame = False ; Don't advance the frame
- mpres@gsnDraw = False ; Don't advance the frame
-
- mpres@tmXTBorderOn = True ; Don't draw top axis.
- mpres@tmXTOn = True
- mpres@tmXBMode = "Explicit" ; Set tick mark mode.
- mpres@tmXBValues = ispan(10940,11090,30)/100.0
- mpres@tmXBLabels = ispan(10940,11090,30)/100.0
- mpres@tmYLMode = "Explicit" ; Set tick mark mode.
- mpres@tmYLValues = ispan(202,215,2)/10.0
- mpres@tmYLLabels = ispan(202,215,2)/10.0
- ;---Zoom in on United States.
- mpres@mpMinLatF = 20.2
- mpres@mpMaxLatF = 21.5
- mpres@mpMinLonF = 109.40
- mpres@mpMaxLonF = 110.9
- mpres@mpCenterLonF = 180
- mpres@mpDataBaseVersion = "HighRes" ; High resolution database
- mpres@mpDataResolution = "Finest"
- mpres@mpDataSetName = "Earth..3" ;
- mpres@mpFillColors = (/"transparent","transparent","gray","transparent"/)
- mpres@tiMainString = "FloodPlain Scatters "
- map = gsn_csm_map(wks,mpres)
-
- gsres = True
- gsres@gsMarkerIndex = 1 ; Use filled dots for markers.
- ;---Get nice spacing through color map for marker colors
- gsn_define_colormap(wks,"Rainbow")
- getvalues wks
- "wkColorMapLen" : clen ; number of colors in color map
- end getvalues
- nstep = (clen-2)/narr
- colors = ispan(2,clen-1,nstep)
- ;---Loop through each "bin" and attach the markers to the map.
- do i = 0, num_distinct_markers-1
- if (.not.ismissing(lat_new_TS(i,0)))
- gsres@gsMarkerColor = colors(i)
- dumstr = unique_string("marker")
- map@$dumstr$ = gsn_add_polymarker(wks,map,lon_new_TS(i,:),lat_new_TS(i,:),gsres)
- end if
- end do
- ;----------------------------------------------------------------------
- ; First version of this plot has a legend at the bottom
- ;----------------------------------------------------------------------
- ;draw(map) ; Drawing the map will draw the markers
- ;draw_legend(wks,lat_new,lon_new,arr,colors) ; Draw a legend at the bottom
- ;frame(wks) ; Advance the frame
- ;----------------------------------------------------------------------
- ; Second version of this plot has a labelbar added.
- ;----------------------------------------------------------------------
- lbid = attach_labelbar(wks,map,arr,colors) ; Attach a labelbar
- draw(map) ; Drawing the map will draw everything
-
- frame(wks) ; Advance the frame.
- end
|
|
窥视卡
雷达卡
发表于 2014-10-8 21:50:53
提升卡
置顶卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
发表于 2014-10-9 11:11:42
~
{:eb502:}{:eb502:}好高大上的感觉