data earthquakes;
  length Type $ 10;
  infile "https://support.sas.com/rnd/app/data/earthquakes.txt" url;
  input Latitude Longitude Depth Magnitude dNearestStation 
        RootMeanSquareTime Type $;
  run;


data hotsprings;
  length Type $ 10;
  infile "https://support.sas.com/rnd/app/data/hotsprings.txt" url;
  input Latitude 1-6 Longitude 8-15 TemperatureFarenheit $ 17-19 
        TemperatureCelsius $ 21-23;
  Type = "hotspring"; 
  run;


data QuakesAndSprings;
   set earthquakes hotsprings;
   where (type in ('hotspring','earthquake'));
run;


data GProjectIn;
   set QuakesAndSprings;
   ID = _N_;
   rename latitude=y longitude=x;
run;

proc gproject data=GProjectIn out=GProjectOut degrees;
   id ID;
run;

data GProjectOut;
   set GProjectOut;
   format _character_;
   informat _character_;
run;


 ods graphics on;
 proc spp data=GProjectOut plots(unpack)=(all observ(attr=mark));
    process p = (X,Y / mark=type)
                / G cross=types('hotspring','earthquake') maxdist=max;
 run;