[ncl-talk] Subsetting a shapefile

[Lyndz]

Dear NCL-experts,

I am plotting a shapefile based on the basic commands provided in the NCL
webpage.

ncl 1> load = "../shapefile_utils.ncl"

ncl 2> sname

ncl 2> sname="PHL_climatetype.shp"

ncl 3> print_shapefile_info(sname)

(0)

(0) Filename: "PHL_climatetype.shp"

(0)   Geometry type: polygon

(0)   # of features: 4

(0)   Min/max lat:      4.59/  21.12

(0)   Min/max lon:    116.93/ 126.61

(0)   Variable names and their types:

(0)       geometry : integer

(0)       segments : integer

(0)       x : double

(0)       y : double

(0)       Id : integer

(0)       Type : double

(0)



I want to do the following:


[1] The "Type" in the shapefile are numbers 1 to 4.

I would like to select only Type 2 and plot the result.


My question is, is it possible to select only the Type 2 in the shapefile
and plot it using the plot_shapefile(sname) command?


[2] I want to save the selected shapefile in [1] as an independent
shapefile. How can I do this in NCL?



I'll appreciate any advise on this.



Sincerely,

Lyndz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20200406/02f51d52/attachment-0001.html>