Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- import = Drop[Import["Fermi2", "Table"], 1];
- ra = Table[(Take[import[[i]], {2, 4}][[1]] +
- Take[import[[i]], {2, 4}][[2]]/60. +
- Take[import[[i]], {2, 4}][[3]]/3600.)/24 360, {i, 1,
- Length[import]}];
- dec = Table[Sign[Take[import[[i]], {5, 7}][[
- 1]]] (Abs[Take[import[[i]], {5, 7}][[1]]] +
- Take[import[[i]], {5, 7}][[2]]/60. +
- Take[import[[i]], {5, 7}][[3]]/3600.), {i, 1, Length[import]}];
- data = Transpose[{dec, ra}];
- GeoGraphics[{Black, Point@GeoPosition@data},
- GeoRange -> {All, {0, 360}}, PlotRangePadding -> Scaled@.01,
- GeoGridLinesStyle -> Directive[Black, Dashed],
- GeoProjection -> "Sinusoidal", GeoGridLines -> Automatic,
- GeoBackground -> White, Frame -> True,
- ImagePadding -> {{60, 40}, {40, 15}}, ImageSize -> 800,
- FrameTicks -> {Table[{N[i Degree], Row[{i/15 + 12, " h"}]}, {i, -180, 180, 30}],
- Table[{N[i Degree], Row[{i, " [Degree]"}]}, {i, -90, 90, 30}]},
- Background -> White, FrameStyle -> Black, TicksStyle -> 15,
- FrameLabel -> {{"DEC", None}, {"RA", None}}]
- frame = GeoGraphics[{White, Point@GeoPosition@data},
- GeoRange -> {All, {0, 360}}, PlotRangePadding -> Scaled@.01,
- GeoGridLinesStyle -> Directive[White, Dashed],
- GeoProjection -> "Sinusoidal", GeoGridLines -> Automatic,
- GeoBackground -> White, Frame -> True,
- ImagePadding -> {{60, 40}, {40, 15}}, ImageSize -> 800,
- FrameTicks -> {Table[{N[i Degree], Row[{i/15 + 12, " h"}]}, {i, -180, 180, 30}], Table[{N[i Degree], Row[{i, " [Degree]"}]}, {i, -90, 90, 30}]},
- Background -> White, FrameStyle -> Black, TicksStyle -> 15,
- FrameLabel -> {{"DEC", None}, {"RA", None}}];
- sky = GeoGraphics[{Black, Point@GeoPosition@data},
- GeoRange -> {All, {0, 360}}, PlotRangePadding -> Scaled@.01,
- GeoGridLinesStyle -> Directive[Black, Dashed],
- GeoProjection -> "Mollweide", GeoGridLines -> Automatic,
- GeoBackground -> White, Frame -> True,
- ImagePadding -> {{60, 40}, {40, 15}}, ImageSize -> 800,
- FrameTicks -> {Table[{N[i Degree], Row[{i/15 + 12, " h"}]}, {i, -180, 180, 30}], Table[{N[i Degree], Row[{i, " [Degree]"}]}, {i, -90, 90, 30}]},
- Background -> White, FrameStyle -> White, TicksStyle -> 15,
- FrameLabel -> {{"DEC", None}, {"RA", None}}];
- sky3 = Show[Graphics @@ frame, Graphics @@ sky, PlotRange -> All]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement