How do I retrieve the actual latitudes and longitudes on a National Park Service shapefile?
The shapefile contains x and y, which gives the leftbc, rightbc, bottombc, and topbc but not
the actual latitudes/longitudes.
Below I include a segment from my NCL script, a dump of the shapefile, and internet links
with sites to download the shapefile and get into on its dimensions/projection (Universal Transverse Mercator).
Thanks,
Michael
----------------------------------------
Where to download info on the shapefile:
http://pubs.usgs.gov/of/2005/1225/Shapefiles/park_boundary.shp.xml
Where to download the shapefile:
http://pubs.usgs.gov/of/2005/1225/Shapefiles/
----------------------------------------
NCL CODE SEGMENT:
f=addfile("park_boundary.shp","r")
segments = f->segments
geometry = f->geometry
segsDims = dimsizes(segments)
geomDims = dimsizes(geometry)
geom_segIndex = f@geom_segIndex
geom_numSegs = f@geom_numSegs
segs_xyzIndex = f@segs_xyzIndex
segs_numPnts = f@segs_numPnts
lines = new(segsDims(0),graphic)
numFeatures = geomDims(0)
lon = f->x
lat = f->y
------------------------------------------------
FILE DUMP:
Variable: f
Type: file
filename: park_boundary
path: park_boundary.shp
file global attributes:
layer_name : park_boundary
geometry_type : polygon
geom_segIndex : 0
geom_numSegs : 1
segs_xyzIndex : 0
segs_numPnts : 1
dimensions:
geometry = 2
segments = 2
num_features = 1 // unlimited
num_segments = 1
num_points = 2860
variables:
integer geometry ( num_features, geometry )
integer segments ( num_segments, segments )
double x ( num_points )
double y ( num_points )
integer OBJECTID ( num_features )
double AREA ( num_features )
double PERIMETER ( num_features )
integer PARKBNDARY ( num_features )
integer PARKBNDA_1 ( num_features )
string ROCKTYPE ( num_features )
double Shape_Leng ( num_features )
double Shape_Area ( num_features )
------------------------------------------------------
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Dec 21 14:32:57 2012
This archive was generated by hypermail 2.1.8 : Fri Jan 04 2013 - 15:32:29 MST