RE: [ncl-talk] overlay shapefile on panel plots

Hi Mary,

It works now.

Thank you very much

Mamadou

________________________________
De : Mary Haley <haley@ucar.edu>
Envoyé : jeudi 10 juillet 2014 13:47
À : mamadoulamine.mbaye@ucad.edu.sn
Cc : ncl-talk@ucar.edu USERS
Objet : Re: [ncl-talk] 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,lnres)
 ​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<mailto:mamadoulamine.mbaye@ucad.edu.sn> <mamadoulamine.mbaye@ucad.edu.sn<mailto: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 Mon Jul 14 01:27:50 2014