topic Re: How to plot a set of unique coordinates on a map in Graphics Programming

How to plot a set of unique coordinates on a map

jserafy — Tue, 06 Feb 2024 22:28:24 GMT

I want to plot unique coordinates on the example map shown at: https://documentation.sas.com/doc/en/pgmsascdc/9.4_3.5/grmapref/n1reiydhvbvs7rn1lxrlps7gwt3l.htm

I would appreciate code to add to that below that places a point on the map generated by the code below.

data us48;
set mapsgfk.us_states;
if state not in (2 15 72);
run;

proc gproject data=us48
latlon
degrees eastlong
out=gulf
longmin=-98
longmax=-81
latmin=25
latmax=33;
where density<5;
id state;
run;

title "Northern Gulf Coast";
pattern value=mempty color=blue;
proc gmap map=gulf data=gulf all;
id state;
choro state / statistic=first nolegend levels=1;
run;
quit;

goptions reset=all;

Re: How to plot a set of unique coordinates on a map

Reeza — Tue, 06 Feb 2024 23:19:41 GMT

In GMAP you use annotate method, which is annoying IMO.

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

I would recommend SGMAP instead, set the bubble size to a fixed value in the data set to have just points.

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

You can find more examples at lexjansen.com and robslink

Re: How to plot a set of unique coordinates on a map

Ksharp — Wed, 07 Feb 2024 02:41:28 GMT

https://blogs.sas.com/content/graphicallyspeaking/2021/08/13/cary-nc-data-car-crashes/
https://blogs.sas.com/content/graphicallyspeaking/2021/08/17/cary-nc-data-fire-incidents/
https://blogs.sas.com/content/graphicallyspeaking/2021/08/31/cary-nc-data-parks-and-trails/
https://blogs.sas.com/content/graphicallyspeaking/2021/10/07/lets-dive-into-some-data-scuba-site-data/

Re: How to plot a set of unique coordinates on a map

jserafy — Wed, 07 Feb 2024 15:02:01 GMT

Thanks Reeza:
Unfortunately, I'm just not bright enough to modify the code shown at the
url's you have kindly provided for my purposes.
Could you possibly help me by further modifying the code below to place the
(water) locations on the map that emerges from the code below?
Sorry, but I just can't seem to make any progress.
j
*--------------------*
*Locations*

28.9787 -95.1672
29.1468 -94.8645
29.2313 -89.8375
29.4695 -94.0852
29.5513 -92.899
29.5282 -93.9293
29.016 -95.098
29.8865 -88.1793
29.5088 -93.9617
29.4177 -92.587
29.1007 -94.02
29.5347 -93.5135
25.2593 -81.9585
28.9765 -91.7018
28.1155 -93.6067
28.3387 -94.4748
28.9765 -91.7018
*SAS Code*

goptions reset=all border device=png;

data us48;
set mapsgfk.us_states;
if state not in (2 15 72);
run;

proc gproject data=us48
latlon
degrees eastlong
out=gulf
longmin=-98
longmax=-68
latmin=25
latmax=36.5;
where density<5;
id state;
run;

title "US Southern Atlantic and Gulf of Mexico";
pattern value=mempty color=blue;
proc gmap map=gulf data=gulf all;
id state;
choro state / statistic=first nolegend levels=1;
run;
quit;

goptions reset=all;
goptions reset=all border device=png;

data us48;
set mapsgfk.us_states;
if state not in (2 15 72);
run;

proc gproject data=us48
latlon
degrees eastlong
out=gulf
longmin=-98
longmax=-68
latmin=25
latmax=36.5;
where density<5;
id state;
run;

title "US Southern Atlantic and Gulf of Mexico";
pattern value=mempty color=blue;
proc gmap map=gulf data=gulf all;
id state;
choro state / statistic=first nolegend levels=1;
run;
quit;

goptions reset=all;

Re: How to plot a set of unique coordinates on a map

jserafy — Wed, 07 Feb 2024 15:09:01 GMT

Thanks Ksharp:
Unfortunately, I'm just not bright enough to modify the code shown at the
url's you have kindly provided for my purposes.
Could you possibly help me by further modifying the code below to place the
(water) locations on the map that emerges from the code below?
Sorry, but I just can't seem to make any progress.
j
*--------------------*
*Locations*

28.9787 -95.1672
29.1468 -94.8645
29.2313 -89.8375
29.4695 -94.0852
29.5513 -92.899
29.5282 -93.9293
29.016 -95.098
29.8865 -88.1793
29.5088 -93.9617
29.4177 -92.587
29.1007 -94.02
29.5347 -93.5135
25.2593 -81.9585
28.9765 -91.7018
28.1155 -93.6067
28.3387 -94.4748
28.9765 -91.7018

*SAS Code*

goptions reset=all border device=png;

data us48;
set mapsgfk.us_states;
if state not in (2 15 72);
run;

proc gproject data=us48
latlon
degrees eastlong
out=gulf
longmin=-98
longmax=-68
latmin=25
latmax=36.5;
where density<5;
id state;
run;

title "US Southern Atlantic and Gulf of Mexico";
pattern value=mempty color=blue;
proc gmap map=gulf data=gulf all;
id state;
choro state / statistic=first nolegend levels=1;
run;
quit;

goptions reset=all;
goptions reset=all border device=png;

data us48;
set mapsgfk.us_states;
if state not in (2 15 72);
run;

proc gproject data=us48
latlon
degrees eastlong
out=gulf
longmin=-98
longmax=-68
latmin=25
latmax=36.5;
where density<5;
id state;
run;

title "US Southern Atlantic and Gulf of Mexico";
pattern value=mempty color=blue;
proc gmap map=gulf data=gulf all;
id state;
choro state / statistic=first nolegend levels=1;
run;
quit;

goptions reset=all;

Re: How to plot a set of unique coordinates on a map

Ksharp — Fri, 09 Feb 2024 02:50:59 GMT

Sorry. I am not familiar with SAS/Graph.

Maybe @GraphGuy (a.k.a Robert.Allision) could give you a hand.

Re: How to plot a set of unique coordinates on a map

GraphGuy — Mon, 12 Feb 2024 14:28:27 GMT

Thanks for the heads-up @Ksharp!

Here's a chapter I wrote on this topic, that might help:

https://robslink.com/SAS/book2/Chapter_07_Annotating_on_Maps.pdf

Re: How to plot a set of unique coordinates on a map

GraphGuy — Mon, 12 Feb 2024 14:48:55 GMT

This code should get you started ...

data coordinates;
input lat long;
datalines;
28.9787 -95.1672
29.1468 -94.8645
29.2313 -89.8375
29.4695 -94.0852
29.5513 -92.899
29.5282 -93.9293
29.016 -95.098
29.8865 -88.1793
29.5088 -93.9617
29.4177 -92.587
29.1007 -94.02
29.5347 -93.5135
25.2593 -81.9585
28.9765 -91.7018
28.1155 -93.6067
28.3387 -94.4748
28.9765 -91.7018
;
run;



data us48;
set mapsgfk.us_states;
if state not in (2 15 72);
run;

proc gproject data=us48 out=gulf
latlon degrees eastlong
longmin=-98
longmax=-68
latmin=25
latmax=36.5
parmout=projparm;
where density<5;
id state;
run;

/* project the coordinates using the same projection parameters as were used to project the map */
proc gproject data=coordinates out=coordinates
latlon degrees eastlong
longmin=-98
longmax=-68
latmin=25
latmax=36.5
parmin=projparm;
id;
run;

/* convert the projected coordinates into annotate dataset */
data coordinates; set coordinates;
xsys='2'; ysys='2'; hsys='3'; when='a';
function='pie'; rotate=360; size=1.0; color='red';
run;

title "US Southern Atlantic and Gulf of Mexico";
pattern value=mempty color=blue;
proc gmap map=gulf data=gulf all anno=coordinates;
id state;
choro state / statistic=first nolegend levels=1;

Re: How to plot a set of unique coordinates on a map

jserafy — Mon, 12 Feb 2024 15:38:40 GMT

Thanks you so much, Graph Guy. Really helps novices like me figure this stuff out.

The lat-lon coordinates should be in the water versus on the land.

Is this a projection problem?

Thanks again.

Re: How to plot a set of unique coordinates on a map

jserafy — Mon, 12 Feb 2024 18:37:01 GMT

Hi GraphGuy:
I get the three warnings below (in blue).
j
---------------------
NOTE: There were 1097 observations read from the data set WORK.GULF.
WHERE resolution<=4;
NOTE: There were 17 observations read from the data set WORK.COORDINATES.
NOTE: PROCEDURE GMAP used (Total process time):
real time 57:55.91
cpu time 2.96 seconds

320 data coordinates;
321 input lat long;
322 datalines;

NOTE: The data set WORK.COORDINATES has 17 observations and 2 variables.
NOTE: DATA statement used (Total process time):
real time 0.00 seconds
cpu time 0.00 seconds

340 ;
341 run;
342
343
344
345 data us48;
346 set mapsgfk.us_states;
347 if state not in (2 15 72);
348 run;

NOTE: There were 33570 observations read from the data set
MAPSGFK.US_STATES.
NOTE: The data set WORK.US48 has 17692 observations and 13 variables.
NOTE: DATA statement used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds

349
350 proc gproject data=us48 out=gulf
351 latlon degrees eastlong
352 longmin=-98
353 longmax=-68
354 latmin=25
355 latmax=36.5
356 parmout=projparm;
357 where density<5;
358 id state;
359 run;

NOTE: PARALLEL1 = 27.875.
NOTE: PARALLEL2 = 33.625.
NOTE: MERIDIAN = -83.
NOTE: The data set WORK.PROJPARM has 1 observations and 9 variables.
NOTE: The data set WORK.GULF has 1976 observations and 13 variables.
NOTE: PROCEDURE GPROJECT used (Total process time):
real time 0.03 seconds
cpu time 0.03 seconds

360 /* project the coordinates using the same projection parameters as
were used to project the map
360! */
361 proc gproject data=coordinates out=coordinates
362 latlon degrees eastlong
363 longmin=-98
364 longmax=-68
365 latmin=25
366 latmax=36.5
367 parmin=projparm;
368 id;
369 run;

WARNING: No entry named COORDINATES found in parameter data set
WORK.PROJPARM.
WARNING: The map data are completely contained within the specified
clipping region.
WARNING: Data NOT clipped due to reasons noted above.
NOTE: PARALLEL1 = 26.4161.
NOTE: PARALLEL2 = 28.7297.
NOTE: MERIDIAN = -88.56285.
NOTE: The data set WORK.COORDINATES has 17 observations and 4 variables.
NOTE: PROCEDURE GPROJECT used (Total process time):
real time 0.00 seconds
cpu time 0.00 seconds

370 /* convert the projected coordinates into annotate dataset */
371 data coordinates; set coordinates;
372 xsys='2'; ysys='2'; hsys='3'; when='a';
373 function='pie'; rotate=360; size=1.0; color='blue';
374 run;

NOTE: There were 17 observations read from the data set WORK.COORDINATES.
NOTE: The data set WORK.COORDINATES has 17 observations and 12 variables.
NOTE: DATA statement used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds

375
376 title "US Southern Atlantic and Gulf of Mexico";
377 pattern value=mempty color=black;
378 proc gmap map=gulf data=gulf all anno=coordinates;
379 id state;
380 choro state / statistic=first nolegend levels=1; run;

NOTE: 50456 bytes written to C:\Users\JOE~1.SER\AppData\Local\Temp\1\SAS
Temporary
Files\_TD17660_SECML4380988JSE_\gmap7.png.

Re: How to plot a set of unique coordinates on a map

jserafy — Mon, 12 Feb 2024 19:35:28 GMT

Thanks so much! Sorry to be a pest, but the points on the map should be on water versus on land.

Perhaps this problem is related to the warnings (bolded) in the sas log below?

---------------

NOTE: There were 1097 observations read from the data set WORK.GULF.
WHERE resolution<=4;
NOTE: There were 17 observations read from the data set WORK.COORDINATES.
NOTE: PROCEDURE GMAP used (Total process time):
real time 57:55.91
cpu time 2.96 seconds

320 data coordinates;
321 input lat long;
322 datalines;

NOTE: The data set WORK.COORDINATES has 17 observations and 2 variables.
NOTE: DATA statement used (Total process time):
real time 0.00 seconds
cpu time 0.00 seconds

340 ;
341 run;
342
343
344
345 data us48;
346 set mapsgfk.us_states;
347 if state not in (2 15 72);
348 run;

NOTE: There were 33570 observations read from the data set MAPSGFK.US_STATES.
NOTE: The data set WORK.US48 has 17692 observations and 13 variables.
NOTE: DATA statement used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds

349
350 proc gproject data=us48 out=gulf
351 latlon degrees eastlong
352 longmin=-98
353 longmax=-68
354 latmin=25
355 latmax=36.5
356 parmout=projparm;
357 where density<5;
358 id state;
359 run;

NOTE: PARALLEL1 = 27.875.
NOTE: PARALLEL2 = 33.625.
NOTE: MERIDIAN = -83.
NOTE: The data set WORK.PROJPARM has 1 observations and 9 variables.
NOTE: The data set WORK.GULF has 1976 observations and 13 variables.
NOTE: PROCEDURE GPROJECT used (Total process time):
real time 0.03 seconds
cpu time 0.03 seconds

360 /* project the coordinates using the same projection parameters as were used to project the map
360! */
361 proc gproject data=coordinates out=coordinates
362 latlon degrees eastlong
363 longmin=-98
364 longmax=-68
365 latmin=25
366 latmax=36.5
367 parmin=projparm;
368 id;
369 run;

WARNING: No entry named COORDINATES found in parameter data set WORK.PROJPARM.
WARNING: The map data are completely contained within the specified clipping region.
WARNING: Data NOT clipped due to reasons noted above.
NOTE: PARALLEL1 = 26.4161.
NOTE: PARALLEL2 = 28.7297.
NOTE: MERIDIAN = -88.56285.
NOTE: The data set WORK.COORDINATES has 17 observations and 4 variables.
NOTE: PROCEDURE GPROJECT used (Total process time):
real time 0.00 seconds
cpu time 0.00 seconds

370 /* convert the projected coordinates into annotate dataset */
371 data coordinates; set coordinates;
372 xsys='2'; ysys='2'; hsys='3'; when='a';
373 function='pie'; rotate=360; size=1.0; color='blue';
374 run;

NOTE: There were 17 observations read from the data set WORK.COORDINATES.
NOTE: The data set WORK.COORDINATES has 17 observations and 12 variables.
NOTE: DATA statement used (Total process time):
real time 0.01 seconds
cpu time 0.01 seconds

375
376 title "US Southern Atlantic and Gulf of Mexico";
377 pattern value=mempty color=black;
378 proc gmap map=gulf data=gulf all anno=coordinates;
379 id state;
380 choro state / statistic=first nolegend levels=1; run;

NOTE: 50456 bytes written to C:\Users\JOE~1.SER\AppData\Local\Temp\1\SAS Temporary
Files\_TD17660_SECML4380988JSE_\gmap7.png.

Re: How to plot a set of unique coordinates on a map

GraphGuy — Tue, 13 Feb 2024 03:37:05 GMT

Ahh - I forgot one thing in the gproject. I had specified the parmin= dataset, but I forgot to tell gproject which entry in that dataset to use. Here's the new/fixed code:

data coordinates;
input lat long;
datalines;
28.9787 -95.1672
29.1468 -94.8645
29.2313 -89.8375
29.4695 -94.0852
29.5513 -92.899
29.5282 -93.9293
29.016 -95.098
29.8865 -88.1793
29.5088 -93.9617
29.4177 -92.587
29.1007 -94.02
29.5347 -93.5135
25.2593 -81.9585
28.9765 -91.7018
28.1155 -93.6067
28.3387 -94.4748
28.9765 -91.7018
;
run;



data us48;
set mapsgfk.us_states;
if state not in (2 15 72);
run;

proc gproject data=us48 out=gulf
latlon degrees eastlong
longmin=-98
longmax=-68
latmin=25
latmax=36.5
parmout=projparm;
where density<5;
id state;
run;
/* project the coordinates using the same projection parameters as were used to project the map */
proc gproject data=coordinates out=coordinates
latlon degrees eastlong
longmin=-98
longmax=-68
latmin=25
latmax=36.5
parmin=projparm parmentry=gulf;
id;
run;
/* convert the projected coordinates into annotate dataset */
data coordinates; set coordinates;
xsys='2'; ysys='2'; hsys='3'; when='a';
function='pie'; rotate=360; size=1.0; color='red';
run;

title "US Southern Atlantic and Gulf of Mexico";
pattern value=mempty color=blue;
proc gmap map=gulf data=gulf all anno=coordinates;
id state;
choro state / statistic=first nolegend levels=1;

Re: How to plot a set of unique coordinates on a map

jserafy — Tue, 13 Feb 2024 15:31:16 GMT

Fantastic! Your help is deeply appreciated. Thanks so much. Now that I have working code , I can try to understand it!

Re: How to plot a set of unique coordinates on a map

GraphGuy — Sat, 17 Feb 2024 15:31:18 GMT

I have a couple of similar examples on this page, that might interest you. (they're on the last row of thumbnails - bottom/left):

https://robslink.com/SAS/democd_mpg/aaaindex.htm

Re: How to plot a set of unique coordinates on a map

jserafy — Tue, 05 Mar 2024 21:45:03 GMT

Hi again, GraphGuy:

Sorry to bother you again, but I can't figure out how to have the dots on the map reflect presence/absence information (0=absence, 1=presence).

All my attempts have failed thus far.

The data and code are below.

--------------------

data coordinates;
input lat long present;
datalines;
28.9787 -95.1672 0
29.1468 -94.8645 0
29.2313 -89.8375 1
29.4695 -94.0852 0
29.5513 -92.899 1
29.5282 -93.9293 1
29.016 -95.098 1
29.8865 -88.1793 1
29.5088 -93.9617 0
29.4177 -92.587 0
29.1007 -94.02 1
29.5347 -93.5135 0
25.2593 -81.9585 1
28.9765 -91.7018 1
28.1155 -93.6067 1
28.3387 -94.4748 0
28.9765 -91.7018 1
;
run;

data us48;
set mapsgfk.us_states;
if state not in (2 15 72);
run;

proc gproject data=us48 out=gulf
latlon degrees eastlong
longmin=-98
longmax=-68
latmin=25
latmax=36.5
parmout=projparm;
where density<5;
id state;
run;
/* project the coordinates using the same projection parameters as were used to project the map */
proc gproject data=coordinates out=coordinates
latlon degrees eastlong
longmin=-98
longmax=-68
latmin=25
latmax=36.5
parmin=projparm parmentry=gulf;
id;
run;
/* convert the projected coordinates into annotate dataset */
data coordinates; set coordinates;
xsys='2'; ysys='2'; hsys='3'; when='a';
function='pie'; rotate=360; size=1.0; color='red';
run;

title "US Southern Atlantic and Gulf of Mexico";
pattern value=mempty color=blue;
proc gmap map=gulf data=gulf all anno=coordinates;
id state;
choro state / statistic=first nolegend levels=1;
run;

Re: How to plot a set of unique coordinates on a map

ballardw — Tue, 05 Mar 2024 22:20:02 GMT

Basic approach is to use If/Then statements in a data step to set a desired color /symbol based on the value of such a variable.

Then use Legend statements to provide description text to match the color /symbol you used.

Re: How to plot a set of unique coordinates on a map

GraphGuy — Wed, 06 Mar 2024 15:46:43 GMT

You could use the value of the 'present' variable to control the color of the dots, such as ...

data coordinates; set coordinates;

length color $12;
xsys='2'; ysys='2'; hsys='3'; when='a';
function='pie'; rotate=360; size=1.0;

if present=0 then color='red';

else if present=1 then color='green';
run;

Re: How to plot a set of unique coordinates on a map

jserafy — Thu, 07 Mar 2024 21:14:35 GMT

Thanks again. Really apreciate the help!