topic Re: Map data points in SAS Programming

Map data points

sascode — Fri, 15 Jul 2022 16:47:02 GMT

Hello ,

attached are two datasets, first one(Baltimore city) contains longitude and latitude for a city in Baltimore, and in second one are few data points with location information which i need to place on map and color by tower Cell.
I am able to do it in R but not in SAS.

Thank you.

Re: Map data points

PaigeMiller — Fri, 15 Jul 2022 16:51:55 GMT

Many of us will not download Excel files as they are security threats. The best way to provide data is via SAS data step code (and not in any other format). Instructions

But also, show us what you have tried. Have you looked at the PROC SGMAP documentation?

Re: Map data points

Ksharp — Sat, 16 Jul 2022 10:42:03 GMT

Calling @Rick_SAS and @GraphGuy

Re: Map data points

Rick_SAS — Sat, 16 Jul 2022 11:18:17 GMT

As others have suggested, start with the doc. Here is an example that overlays points and lines on a map:

https://go.documentation.sas.com/doc/en/pgmsascdc/9.4_3.5/grmapref/n0x9p71lic5q5an15063k0qjo8er.htm

You should also look at the dozens of examples in the SAS blogs. For example, Google quickly located this one:

https://blogs.sas.com/content/sastraining/2017/03/29/plotting-markers-on-a-map-at-zip-code-locations-using-gmap-or-sgplot/

Re: Map data points

sascode — Tue, 19 Jul 2022 02:29:17 GMT

Thank you for your suggestions.
Much appreciated it.

Re: Map data points

GraphGuy — Wed, 20 Jul 2022 02:23:53 GMT

As PaigeMiller mentioned, it would be interesting to see what you've tried.

That being said, here is one solution ...

proc import datafile="/home/u59782200/datapoints.xlsx" out=datapoints dbms=xlsx;
run;

proc import datafile="/home/u59782200/balt_city.xlsx" out=balt_city dbms=xlsx;
run;

proc sgmap mapdata=balt_city plotdata=datapoints;
openstreetmap;
choromap / mapid=subregion lineattrs=(color=red thickness=2px);
scatter x=long y=lat / group=tower 
 markerattrs=(symbol=circlefilled size=12px) 
 /* More current versions let you specify colors via colormodel=(cxd7191c yellow) */
 ;
run;