Load the azimuth and elevation data from a logfile generated by an Adafruit® GPS satellite sensor. The data provided in this example contains the azimuth and elevation of each satellite and the pseudorandom noise (PRN) codes. Store these values as vectors.
load('gpsHWInfo','hwInfo')
satAz = hwInfo.SatelliteAzimuths;
satEl = hwInfo.SatelliteElevations;
prn = hwInfo.SatellitePRNs;
Separate the satellites based on the PRN codes. To correlate each position with a group, create a categorical array. PRNs greater than 32 are part of the wide area augmentation system (WAAS).
isWAAS = (prn > 32);
constellationGroup = categorical(isWAAS,[false,true],{'GPS','WAAS'});
Visualize the satellites and specify the categorical groups in the GroupData name-value argument. Specify the PRN as the label for each point. Show the legend.
skyplot(satAz,satEl,prn,'GroupData',constellationGroup)
legend('GPS','WAAS')