Here is my code.
(* Create empty lists *)
xlist={};
ylist={};
zlist={};
(* Take inputs for ECEF *)
(Label[startinput];x = Input["What is the x coordinate?"];
y = Input["What is the y coordinate?"];
z = Input["What is the z coordinate?"];
AppendTo [xlist, x]
AppendTo [ylist, y]
AppendTo [zlist, z]
(Label[errorcontrol];
again=InputString["Enter more coords? Y for yes, N for no"];
If [again=="Y", Print["Okay!"];Goto[startinput]];
If [again!="N", Print["Sorry, this is an invalid input here!"];Goto[errorcontrol]];
If [again=="N", Print["Okay."]]));
(* Debugging *)
Print[xlist]
Print[ylist]
Print[zlist]
Print[x]
Print[y]
Print[z]
(* Combine Lists *)
list=Transpose[{xlist,ylist,zlist}]
(* More Debugging *)
Print[list]
(* Put in Coordinate Form *)
GeoPositionXYZ[{x, y, z}, "ITRF00"];
(* Convert to LLA *)
GeoPosition[%]
(* Display Map *)
GeoGraphics[GeoMarker[GeoPosition[%]],GeoRange -> "World", GeoProjection -> "Robinson"]
(* Display 3D Plot *)
ListPointPlot3D[list, DataRange -> {{-1000000, 1000000}, {-1000000, 1000000}}]
The issue is: The map only displays the very last point entered. It should show all the points but I have absolutely no idea how to convert an entire list of cartesian coordinates to latitude longitude and then display them all on the map. Thanks for the help.