[ncl-talk] shapefiles
Vanúcia Schumacher
vanucia-schumacher at hotmail.com
Thu Jan 31 05:58:47 MST 2019
Unfortunately I do not have another shapefile, just this one referring to my work.
________________________________
De: Dave Allured - NOAA Affiliate <dave.allured at noaa.gov>
Enviado: terça-feira, 29 de janeiro de 2019 16:16
Para: Vanúcia Schumacher
Cc: ncl-talk at ucar.edu
Assunto: Re: [ncl-talk] shapefiles
Vanúcia,
If you can find another shape file of the same region from a different source, give that a try. This might be quicker than waiting for analysis of the original shape file. If you find a working shape file, please let us know.
--Dave
On Mon, Jan 28, 2019 at 10:50 AM Dave Allured - NOAA Affiliate <dave.allured at noaa.gov<mailto:dave.allured at noaa.gov>> wrote:
Thanks. I can't look at this today. I will get back to you in a day or two. Perhaps someone else would like to check out the shape file and the script.
On Mon, Jan 28, 2019 at 10:45 AM Vanúcia Schumacher <vanucia-schumacher at hotmail.com<mailto:vanucia-schumacher at hotmail.com>> wrote:
Is attached.
Thank you for the advance
________________________________
De: Dave Allured - NOAA Affiliate <dave.allured at noaa.gov<mailto:dave.allured at noaa.gov>>
Enviado: segunda-feira, 28 de janeiro de 2019 15:37
Para: Vanúcia Schumacher
Cc: ncl-talk at ucar.edu<mailto:ncl-talk at ucar.edu>
Assunto: Re: [ncl-talk] shapefiles
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<mailto: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<mailto:dave.allured at noaa.gov>>
Enviado: domingo, 27 de janeiro de 2019 18:23
Para: Vanúcia Schumacher
Cc: ncl-talk at ucar.edu<mailto: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/20190131/c4f14c1c/attachment.html>
More information about the ncl-talk
mailing list