NPS shapefile

From: Michael Notaro <mnotaro_at_nyahnyahspammersnyahnyah>
Date: Fri Dec 21 2012 - 14:32:34 MST

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