Re: overlay shapefile on panel plots

Hi Mamadou,

You can use the gsn_add_shapefile_polyline function in your panel code.
 You need to call it right after the call to gsn_csm_contour_map_ce, but
before you panel the plots, and make sure that the variable on the left
side of the "=" sign is unique in every case.

You already have code that's calling "new" for each set of plots you need
to create, so you can add the unique variables there:

  plot1 = new(6,graphic)
  plot2 = new(4,graphic)
  plot3 = new(3,graphic)
  plot4 = new(2,graphic)
  shp1 = new(6,graphic)
  shp2 = new(4,graphic)
  shp3 = new(3,graphic)
  shp3 = new(2,graphic)

Then, after you are done calling "gsn_csm_contour_map_ce" for a group of
plots, you can use a do loop to add the shapefile outlines. Here's an
example for the "plot1" plots:

 do i=0,5
   shp1(i) = gsn_add_shapefile_polylines(wks,plot1(i),shapefile_name,lnre=
s)
 ​end do​

This calls for an enhancement to "gsn_add_shapefile_polylines" that would
allow you to pass in the whole array of plots, rather than having to call
it every time. I will add a ticket for this. You would then be able to
just have one call:

   shp1 = gsn_add_shapefile_polylines(wks,plot1,shapefile_name,lnres)

and you wouldn't need to create "shp1" beforehand.

--Mary



​



On Wed, Jul 9, 2014 at 11:53 AM, mamadoulamine.mbaye@ucad.edu.sn <
mamadoulamine.mbaye@ucad.edu.sn> wrote:

> Hi,
>
> I would like to overlay the shapefile (toto2.png) on each plot of the
> panel (toto1.png). these are the separated scripts:
>
> 1-panel plots:
>
> wks = gsn_open_wks("png", "toto1")
>
> plot1 = new(6,graphic)
> plot2 = new(4,graphic)
> plot3=new(3,graphic)
> plot4=new(2,graphic)
>
> res = True
> res@gsnMaximize = True
> res@gsnDraw = False
> res@gsnFrame = False
> res@cnInfoLabelOn = False
> res@cnFillOn = True
> res@cnLinesOn = False
> res@cnLineLabelsOn = False
> res@gsnAddCyclic = False
> res@gsnSpreadColors = True
> res@cnFillMode = "RasterFill"
>
> res@lbLabelBarOn = False
> res@mpLimitMode="LatLon"
> res@mpMinLonF=-18
> res@mpMaxLonF=-6
> res@mpMinLatF=9
> res@mpMaxLatF=19
> res@tmXBTickSpacingF= 3
> res@tmYLTickSpacingF = 3
> ;res@tmXBLabelStride = 2
> res@gsnRightString=""
>
> res@tmXBLabelFontHeightF =0.025
> res@tmYLLabelFontHeightF =0.025
>
>
> res@cnLevelSelectionMode = "ManualLevels"
> res@cnMinLevelValF = 22.
> res@cnMaxLevelValF = 38.
> res@cnLevelSpacingF = 2.
>
>
> res@txFontHeightF = 0.03
> res@gsnLeftString = "UC_HISTO"
> plot1(0) = gsn_csm_contour_map_ce(wks,tas1,res)
>
> res@gsnLeftString = "BC_HISTO"
> plot1(1) = gsn_csm_contour_map_ce(wks,tas4,res)
>
> res@gsnLeftString = "UC_RCP45"
> plot1(2) = gsn_csm_contour_map_ce(wks,tas2,res)
>
> res@gsnLeftString = "BC_RCP45"
> plot1(3) = gsn_csm_contour_map_ce(wks,tas5,res)
>
> res@gsnLeftString = "UC_RCP85"
> plot1(4) = gsn_csm_contour_map_ce(wks,tas3,res)
>
> res@gsnLeftString = "BC_RCP85"
> plot1(5) = gsn_csm_contour_map_ce(wks,tas6,res)
>
>
> ;-----------------------------------------------Create
> Panel1-----------------------------------------------
> resP = True
> resP@gsnMaximize = True
> resP@lbAutoManage = False
> resP@gsnFrame = False
> resP@gsnPanelLabelBar = True ;
> add common colorbar
> resP@lbLabelFontHeightF = 0.008
> resP@gsnPanelTop = 0.97
> resP@gsnPanelBottom =0.43
> resP@gsnPanelYWhiteSpacePercent = 2
> resP@gsnPanelXWhiteSpacePercent = 2
> resP@lbTitleString ="(deg C)"
> resP@lbTitlePosition = "Right"
> resP@lbTitleFontHeightF = 0.008
> resP@lbTitleDirection = "Across"
> ;-----------------------------------------------Create
> Panel2-----------------------------------------------
> ;=========================
====== for
> baises=======================
======================
> res@cnLevelSelectionMode = "ManualLevels"
> res@cnMinLevelValF = 1.
> res@cnMaxLevelValF = 6.
> res@cnLevelSpacingF = 1.
> res@gsnLeftString = "Change_UC_RCP45"
> plot2(0) = gsn_csm_contour_map_ce(wks,tas7,res)
>
> res@gsnLeftString = "Change_BC_RCP45"
> plot2(1) = gsn_csm_contour_map_ce(wks,tas8,res)
>
> res@gsnLeftString = "Change_UC_RCP85"
> plot2(2) = gsn_csm_contour_map_ce(wks,tas9,res)
>
> res@gsnLeftString = "Change_BC_RCP85"
> plot2(3) = gsn_csm_contour_map_ce(wks,tas10,res)
>
> ;=========================
==Panel for biases
> =========================
=================
> resP2 =True
> resP2@gsnMaximize =True
> resP2@lbAutoManage =False
> resP2@gsnFrame =False
> resP2@gsnPanelTop =0.42
> resP2@gsnPanelBottom =0.042
> resP2@gsnPanelLabelBar =True
> resP2@txFontHeightF =0.012 ;
> add common colorbar
> resP2@lbLabelFontHeightF =0.0075
> resP2@gsnPanelYWhiteSpacePercent = 2
> resP2@gsnPanelXWhiteSpacePercent = 2
> resP2@lbTitleString ="(deg C)"
> resP2@lbTitlePosition = "Right"
> resP2@lbTitleFontHeightF = 0.008
> resP2@lbTitleDirection = "Across"
> ;-------------------------------------------------DELTA_1
> ------------------------------------------
> res@cnLevelSelectionMode = "ManualLevels"
> res@cnMinLevelValF = -2.
> res@cnMaxLevelValF = 1.
> res@cnLevelSpacingF = 1.
>
> res@gsnLeftString = "Delta_HISTO"
> plot3(0) = gsn_csm_contour_map_ce(wks,d10,res)
>
> res@gsnLeftString = "Delta_RCP45"
> plot3(1) = gsn_csm_contour_map_ce(wks,d20,res)
>
> res@gsnLeftString = "Delta_RCP85"
> plot3(2) = gsn_csm_contour_map_ce(wks,d30,res)
>
>
> ;=========================
==Panel delta_1
> =========================
=================
> resP3 =True
> resP3@gsnMaximize =True
> resP3@lbAutoManage =False
> resP3@gsnFrame =False
> resP3@gsnPanelLeft =0.69
> resP3@gsnPanelRight =0.88
> resP3@gsnPaneltop =0.98
> resP3@gsnPanelBottom =0.4
> resP3@gsnPanelLabelBar =True
> resP3@txFontHeightF =0.01 ;
> add common colorbar
> resP3@lbLabelFontHeightF =0.007
> resP3@gsnPanelYWhiteSpacePercent = 2
> resP3@gsnPanelXWhiteSpacePercent = 2
> resP3@lbTitleString ="(deg C)"
> resP3@lbTitlePosition = "Right"
> resP3@lbTitleFontHeightF = 0.008
> resP3@lbTitleDirection = "Across"
>
> ;-----------------------------------------DELTA_2-----------------------=
---------------------------
>
> res@cnLevelSelectionMode = "ManualLevels"
> res@cnMinLevelValF = -1.8
> res@cnMaxLevelValF = -0.4
> res@cnLevelSpacingF = 0.4
> res@cnFillColors = (/2,18,40,70,100,128/)
> res@gsnLeftString = "Delta_Change_RCP45"
>
> plot4(0) = gsn_csm_contour_map_ce(wks,d40,res)
>
> res@gsnLeftString = "Delta_Change_RCP85"
> plot4(1) = gsn_csm_contour_map_ce(wks,d50,res)
>
>
>
> ;=========================
==Panel delta_2
> =========================
=================
> resP4 =True
> resP4@gsnMaximize =True
> resP4@lbAutoManage =False
> resP4@gsnFrame =False
> resP4@gsnPanelLeft =0.044
> resP4@gsnPanelRight =0.46
> resP4@gsnPanelLabelBar =True
> resP4@txFontHeightF =0.012 ;
> add common colorbar
> resP4@lbLabelFontHeightF =0.0075
> resP4@gsnPanelLeft =0.69
> resP4@gsnPanelRight =0.88
> resP4@gsnPanelTop =0.46
> resP4@gsnPanelBottom =0.001
> resP4@gsnPanelYWhiteSpacePercent = 2
> resP4@gsnPanelXWhiteSpacePercent = 2
> resP4@lbTitleString ="(deg C)"
> resP4@lbTitlePosition = "Right"
> resP4@lbTitleFontHeightF = 0.008
> resP4@lbTitleDirection = "Across"
>
>
> gsn_panel(wks,plot1,(/3,2/),resP)
> gsn_panel(wks,plot2,(/2,2/),resP2)
> gsn_panel(wks,plot3,(/3,1/),resP3)
>
> gsn_panel(wks,plot4,(/2,1/),resP4)
>
> frame(wks)
> end
>
> 2-shapefile:
>
> wks = gsn_open_wks("png","toto2") ; Open X11 window for graphics
>
> res = True ; Set some plot options
>
> res@gsnMaximize = True ; maximize plot in frame
>
> res@mpMinLatF = 10.30
> res@mpMaxLatF =17.7
> res@mpMinLonF = -18.0
> res@mpMaxLonF = -6.5
>
> res@mpDataBaseVersion = "MediumRes" ; "LowRes", "HighRes"
>
> ;res@tiMainString = shapefile_name ; Main title
>
>
> plot = gsn_csm_map(wks,res)
>
> lnres = True ; Set some line options
>
> lnres@gsLineThicknessF = 2.0 ; default is 1.0
> lnres@gsLineColor = "NavyBlue" ; default is black
> ; lnres@gsLineDashPattern = 2 ; default is 0 (solid)
>
> id = gsn_add_shapefile_polylines(wks,plot,shapefile_name,lnres)
>
> draw(plot)
> frame(wks)
>
> end
>
>
> Thanks in advance
>
> Mamadou,
>
> Regards
>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

Received on Thu Jul 10 08:47:32 2014