[ncl-talk] shapefiles
Dave Allured - NOAA Affiliate
dave.allured at noaa.gov
Mon Jan 28 10:37:39 MST 2019
Vanúcia,
Shape files come in different varieties. It looks like example
shapefiles_21 may need an upgrade to be able to read your shape file. If
you can send a copy of the shape file, and also your current script
shape21.ncl, I will test it on my system. Please do not include the
scientific data file; I do not need it for testing the shape file. The NCL
support team might need to fix this.
--Dave
On Mon, Jan 28, 2019 at 9:47 AM Vanúcia Schumacher <
vanucia-schumacher at hotmail.com> wrote:
> Thanks for the tip.
> I'm trying to use the script from example shapefiles_21, but I'm having
> this error in function. I can not fix it.
>
> fatal:Subscript out of range, error in subscript #0
> fatal:An error occurred reading features
> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 56 in
> file shape21.ncl
>
> fatal:["Execute.c":8640]:Execute: Error occurred at or near line 250 in
> file shape21.ncl
>
>
> The lines correspond to:
>
>
> function get_areas_of_interest(shp_file_name,shp_var_name,opt[1]:logical)
> begin
> ;---Open the shapefile
> f = addfile(shp_file_name,"r")
> features = f->$shp_var_name$
>
> if(opt.and.isatt(opt,"areas_to_exclude")) then
> features at _FillValue = default_fillvalue(typeof(features))
> do na=0,dimsizes(opt at areas_to_exclude)-1
> ii := ind(features.eq.opt at areas_to_exclude(na))
> if(.not.any(ismissing(ii))) then
> features(ii) = features at _FillValue
> end if
> end do
> return(features(ind(.not.ismissing(features)))) ------> line 56
> else
> return(features)
> end if
> end
>
> opt = True
> opt at areas_to_exclude = areas_of_interest
> areas_to_fill =
> get_areas_of_interest(shp_filename1,shape_var_name,opt) ----> line 250
>
> Infos:
>
> 0) Terrain Height (m) : min=0 max=5070.39
> (0) ======================================================================
> (0) Filename: "cuencasx.shp"
> (0) Geometry type: polygon
> (0) # of features: 4
> (0) Min/max lat: -34.77/ -34.23
> (0) Min/max lon: -70.53/ -70.13
> (0) Variable names and their types:
> (0) geometry : integer
> (0) segments : integer
> (0) x : double
> (0) y : double
> (0) COD_CUEN : string
> (0) COD_SUBC : string
> (0) COD_SSUBC : string
> (0) NOMBRE : string
> (0) AREAKM2 : double
> (0) ======================================================================
>
> Variable: features
> Type: string
> Total Size: 32 bytes
> 4 values
> Number of Dimensions: 1
> Dimensions and sizes: [num_features | 4]
> Coordinates:
> Number Of Attributes: 0
> (0) Cuenca Cipreses
> (1) Cuenca Coton
> (2) Cuenca Cortaderal
> (3) Cuenca Universidad
> (0) ==================================================
> (0) Shapefile: cuencasx.shp
> (0) Areas of interest: Cuenca Cipreses,Cuenca Coton,Cuenca
> Cortaderal,Cuenca Universidad
> (0) min_lat_chk: -37.4943
> (0) max_lat_chk: -31.9474
> (0) min_lon_chk: -73.6526
> (0) max_lon_chk: -66.7254
> (0) min_lat_data: -37.4943
> (0) max_lat_data: -31.9474
> (0) min_lon_data: -73.6526
> (0) max_lon_data: -66.7254
> (0) 4761 data values originally
> (0) Will keep data values inside given shapefile areas
> (0) ==================================================
> (0) 13 data values kept
> (0) shapefile_mask_data: elapsed time: 0.00758901 CPU seconds.
> (0) ==================================================
>
>
> *De:* Dave Allured - NOAA Affiliate <dave.allured at noaa.gov>
> *Enviado:* domingo, 27 de janeiro de 2019 18:23
> *Para:* Vanúcia Schumacher
> *Cc:* ncl-talk at ucar.edu
> *Assunto:* Re: [ncl-talk] shapefiles
>
> Vanúcia,
>
> Your two data files are on different grids. If you only need to make
> plots that show the same region on a map, then use graphics masking, not
> data masking. See the third plot in example shapefiles_21 on this page:
>
> https://www.ncl.ucar.edu/Applications/shapefiles.shtml
>
> If you need to perform calculations between the two data sets, then you
> will need to regrid the data from one file to match the other grid. Please
> see this documentation for regridding:
>
> https://www.ncl.ucar.edu/Applications/regrid.shtml
> https://www.ncl.ucar.edu/Applications/ESMF.shtml
>
> --Dave
>
>
> On Sun, Jan 27, 2019 at 12:40 PM Vanúcia Schumacher <
> vanucia-schumacher at hotmail.com> wrote:
>
> Hi all,
>
> I need support to find the problem with my script, with the purpose of
> cutting two different data in a shapefile, but keeping the information
> such as time lat and lon.
> My script (attachment) is "cutting" the region from shapefile different
> for each input (data) (see Figure attachment).
> I'd like to select the shapefile region of the same size independent of
> the input file, and keep the time, lat, and lon information.
>
> I appreciate any help
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/60534b37/attachment.html>
More information about the ncl-talk
mailing list