[ncl-talk] How to avoid Interpolation outside shapefile
Rick Brownrigg
brownrig at ucar.edu
Wed Jan 8 16:19:59 MST 2020
Can you either i) provide the shapefile (and I'll point out that means the
.shp, .shx, .dbf files), or perhaps show us the results of
ncl_filedump District_Boundary.shp
Rick
On Wed, Jan 8, 2020 at 7:24 AM Ajay Bankar via ncl-talk <ncl-talk at ucar.edu>
wrote:
> Dear NCL Users,
> I'm trying to make spatial plot of rainfall
> from station data available in csv file. I tried different options in
> resources but data is not masking with shapefile.
> I have attached plot and data file for reference.
>
> Thanks for any help!
>
> Here is my script
>
> load "/shp.ncl"
>
> begin
> ;-- shapefile name
> shp_filename =
> "/scratch/Administrative_Boundary_Headquarters/District_Boundary.shp"
>
> infile = "station_data.csv"
> lines = asciiread(infile, -1, "string")
>
> rain = tofloat( str_get_field(lines(1:), 7, ","))
> rain at lat1d = tofloat( str_get_field(lines(1:), 5, ","))
> rain at lon1d = tofloat( str_get_field(lines(1:), 6, ","))
> ;print(rain)
>
> ;-- define sub-region, here for Karnataka
> minlat = 11
> maxlat = 19
> minlon = 73.5
> maxlon = 79.1
>
> ;-- open workstation
> wks_type = "pdf"
> wks = gsn_open_wks(wks_type, "station_plot")
> ;-- resource settings
> res = True
> res at gsnDraw = False ;-- don't
> draw plot yet
> res at gsnFrame = False ;-- don't
> advance frame yet
> res at gsnMaximize = True ;--
> maximize plot in frame
>
> ; res at tfDoNDCOverlay = True
>
> cmap := read_colormap_file("BlAqGrYeOrReVi200")
> cmap(0,:) = (/0,0,0,0/) ; make first color fully transparent
> res at cnFillOn = True
> res at cnLinesOn = False
> res at cnLineLabelsOn = False
> res at cnFillPalette = cmap
> ;"BlueYellowRed"
> res at cnLevelSelectionMode = "ExplicitLevels" ;--
> manually set contour levels
> res at cnLevels = (/5,20,35,50,65,80,95,110,125,140,155/)
> res at gsnAddCyclic = False
>
> res at mpGridLineColor = 0 ;"grey40"
> res at mpFillOn = False
> res at mpOutlineOn = False
> res at mpGeophysicalLineColor = "black"
> res at mpLimitMode = "LatLon"
> res at mpMinLatF = minlat
> res at mpMaxLatF = maxlat
> res at mpMinLonF = minlon
> res at mpMaxLonF = maxlon
> ;res at mpAreaMaskingOn = True
>
> res at lbBoxMinorExtentF = 0.1 ;-- decrease
> height of labelbar boxes
> res at pmLabelBarOrthogonalPosF = 0.08 ;-- move
> labelbar to the left side of plot
> res at lbBoxMinorExtentF = 0.2 ;-- decrease
> height of labelbar boxes
>
> res at pmTickMarkDisplayMode = "Always"
> res at tmXTOn = False
> res at tiMainString = "STATION DATA PLOT"
>
> ;-- create plot
> plot = gsn_csm_contour_map(wks,rain,res)
> var_mask = shapefile_mask_data(rain,shp_filename,True)
>
> ;-- create contours of masked data
> plot_mask = gsn_csm_contour_map(wks,var_mask,res)
>
> ;-- add plot_mask to plot
> lnres = True
> lnres at gsLineColor = "black"
> lnres at gsLineThicknessF = 1.0
>
> plot at lines = gsn_add_shapefile_polylines(wks, plot_mask,
> shp_filename, lnres)
>
> ;-- draw the plot and advance the frame
> draw(plot_mask)
> frame(wks)
> end
> --
>
> *Thanks & Regards,*
>
> *Ajay*
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200108/057abcf6/attachment.html>
More information about the ncl-talk
mailing list