[ncl-talk] shapefiles
Vanúcia Schumacher
vanucia-schumacher at hotmail.com
Mon Jan 28 10:45:16 MST 2019
Is attached.
Thank you for the advance
________________________________
De: Dave Allured - NOAA Affiliate <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
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/20190128/3db0ef77/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cuencasx.cpg
Type: application/octet-stream
Size: 5 bytes
Desc: cuencasx.cpg
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/3db0ef77/attachment.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cuencasx.dbf
Type: application/x-dbf
Size: 802 bytes
Desc: cuencasx.dbf
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/3db0ef77/attachment.bin>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cuencasx.prj
Type: application/octet-stream
Size: 409 bytes
Desc: cuencasx.prj
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/3db0ef77/attachment-0001.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cuencasx.sbn
Type: application/octet-stream
Size: 156 bytes
Desc: cuencasx.sbn
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/3db0ef77/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cuencasx.sbx
Type: application/octet-stream
Size: 116 bytes
Desc: cuencasx.sbx
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/3db0ef77/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cuencasx.shp
Type: application/octet-stream
Size: 20308 bytes
Desc: cuencasx.shp
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/3db0ef77/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cuencasx.shp.xml
Type: text/xml
Size: 1510 bytes
Desc: cuencasx.shp.xml
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/3db0ef77/attachment.xml>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: cuencasx.shx
Type: application/octet-stream
Size: 132 bytes
Desc: cuencasx.shx
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/3db0ef77/attachment-0005.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: shape21.ncl
Type: application/octet-stream
Size: 11064 bytes
Desc: shape21.ncl
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190128/3db0ef77/attachment-0006.obj>
More information about the ncl-talk
mailing list