[ncl-talk] shapefiles

[Vanúcia Schumacher]

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<mailto: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/414aca00/attachment.html>