| Date: | Wed, 4 Apr 2007 13:33:14 -0700 |
| Reply-To: | "Ferguson, Deron (OFM)" <Deron.Ferguson@OFM.WA.GOV> |
| Sender: | "SAS(r) Discussion" <SAS-L@LISTSERV.UGA.EDU> |
| From: | "Ferguson, Deron (OFM)" <Deron.Ferguson@OFM.WA.GOV> |
| Subject: | Proc mapimport question |
| Content-Type: | text/plain; charset="us-ascii" |
Hi all,
I'm importing unprojected (GCS83) shapefiles using proc mapimport and
I've come against a bug -- perhaps someone knows of a remedy:
If I use proc gmap with the imported shapefile, I get a decent looking
map, except that there are a few boundaries that stretch across the map
as straight lines, i.e., as if one of the lat/longs defining a polygon
was not imported correctly. I tried repairing the shapefile geometry in
ArcGIS and re-importing the shapefile to SAS, but that had no effect
(there were a couple of errors found and fixed in the repair).
Of note: When I look at the imported dataset, there are a number of x /
y pairs that are null. I can't tell off hand if these are the source of
this problem, and if so, how to remove them (I'm not too concerned about
the loss of accuracy of a single vertex here an there in a given
polygon). I have tried deleting the rows with null x/y sets, but that
makes the map come out even worse with proc gmap.
Anyone else have this problem?
Thanks,
Deron
|