|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"|
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?