The surface positions of a calculated trajectory are now shown thicker and have a date.

lib/+artoa/+controller/+track/+trajectoryOutput/plot.m

@@ -52,6 +52,16 @@ for i = 1:length(trajectories)
         trajectories{i}, ...
         trajectoryColor ...
     );
+    % plot surface positions thicker than trajectory
+    [ ...
+        artoaGui.trajectoryOutput.trajectoryHandles{i}.scatterSurfacePositions, ...
+        artoaGui.trajectoryOutput.trajectoryHandles{i}.textSurfacePositions ...
+    ] = artoa.controller.track.trajectoryOutput.plotTrajectorySurfacePositions( ...
+        artoaGui.trajectoryOutput.axesTrajectoryOutput, ...
+        trajectories{i}, ...
+        artoaWorkspace.satData, ...
+        trajectoryColor ...
+    );
     % plot trajectory reference point
     [ ...
         artoaGui.trajectoryOutput.trajectoryHandles{i}.scatterReferencePoints, ...
@@ -83,7 +93,9 @@ for i = 1:length(trajectories)
             artoaGui.trajectoryOutput.trajectoryHandles{i}.linePositions, ...
             artoaGui.trajectoryOutput.trajectoryHandles{i}.textPositions, ...
             artoaGui.trajectoryOutput.trajectoryHandles{i}.scatterReferencePoints, ...
-            artoaGui.trajectoryOutput.trajectoryHandles{i}.textReferencePoints ...
+            artoaGui.trajectoryOutput.trajectoryHandles{i}.textReferencePoints, ...
+            artoaGui.trajectoryOutput.trajectoryHandles{i}.scatterSurfacePositions, ...
+            artoaGui.trajectoryOutput.trajectoryHandles{i}.textSurfacePositions ...
         }, ...
         ~trajectories{i}.hidden ...
     );

lib/+artoa/+controller/+track/+trajectoryOutput/plotTrajectorySurfacePositions.m

+function [scatterHandle, textHandles] = plotTrajectorySurfacePositions(pAxesHandle, pTrajectory, pSatData, pColor)
+%PLOTTRAJECTORYSURFACEPOSITIONS Summary of this function goes here
+%   Detailed explanation goes here
+
+%% Get required data
+satDates = artoa.convert.dmy2rd(pSatData.day_sat, pSatData.month_sat, pSatData.year_sat);
+trajectoryDates = pTrajectory.date;
+[~, ~, surfaceDates] = intersect(round(satDates), round(trajectoryDates));
+x = pTrajectory.longitude(surfaceDates);
+y = pTrajectory.latitude(surfaceDates);
+
+%% Hold axes
+hold(pAxesHandle, 'on');
+
+scatterHandle = scatter( ...
+    pAxesHandle, ...
+    x, ...
+    y, ...
+    50, ...
+    pColor, ...
+    'filled' ...
+);
+
+t = artoa.convert.rd2dmy(trajectoryDates(surfaceDates));
+trd = round(trajectoryDates(surfaceDates));
+textHandles = cell(1, length(t));
+for i = 1:length(t)
+    textHandles{i} = text( ...
+        x(i), ...
+        y(i), ...
+        {' '; [num2str(t(i, 3)) '-' pad(num2str(t(i, 2)), 2, 'left', '0') '-' pad(num2str(t(i, 1)), 2, 'left', '0') ', ' num2str(trd(i))]}, ...
+        'FontSize', 10, ...
+        'Color', pColor ...
+    );
+end
+
+%% Resume axes
+hold(pAxesHandle, 'off');
+
+end
+