Advertisement
Guest User

Untitled

a guest
Aug 1st, 2015
217
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 2.17 KB | None | 0 0
  1. import = Drop[Import["Fermi2", "Table"], 1];
  2.  
  3. ra = Table[(Take[import[[i]], {2, 4}][[1]] +
  4. Take[import[[i]], {2, 4}][[2]]/60. +
  5. Take[import[[i]], {2, 4}][[3]]/3600.)/24 360, {i, 1,
  6. Length[import]}];
  7.  
  8. dec = Table[Sign[Take[import[[i]], {5, 7}][[
  9. 1]]] (Abs[Take[import[[i]], {5, 7}][[1]]] +
  10. Take[import[[i]], {5, 7}][[2]]/60. +
  11. Take[import[[i]], {5, 7}][[3]]/3600.), {i, 1, Length[import]}];
  12.  
  13. data = Transpose[{dec, ra}];
  14.  
  15. GeoGraphics[{Black, Point@GeoPosition@data},
  16. GeoRange -> {All, {0, 360}}, PlotRangePadding -> Scaled@.01,
  17. GeoGridLinesStyle -> Directive[Black, Dashed],
  18. GeoProjection -> "Sinusoidal", GeoGridLines -> Automatic,
  19. GeoBackground -> White, Frame -> True,
  20. ImagePadding -> {{60, 40}, {40, 15}}, ImageSize -> 800,
  21. FrameTicks -> {Table[{N[i Degree], Row[{i/15 + 12, " h"}]}, {i, -180, 180, 30}],
  22. Table[{N[i Degree], Row[{i, " [Degree]"}]}, {i, -90, 90, 30}]},
  23. Background -> White, FrameStyle -> Black, TicksStyle -> 15,
  24. FrameLabel -> {{"DEC", None}, {"RA", None}}]
  25.  
  26. frame = GeoGraphics[{White, Point@GeoPosition@data},
  27. GeoRange -> {All, {0, 360}}, PlotRangePadding -> Scaled@.01,
  28. GeoGridLinesStyle -> Directive[White, Dashed],
  29. GeoProjection -> "Sinusoidal", GeoGridLines -> Automatic,
  30. GeoBackground -> White, Frame -> True,
  31. ImagePadding -> {{60, 40}, {40, 15}}, ImageSize -> 800,
  32. 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}]},
  33. Background -> White, FrameStyle -> Black, TicksStyle -> 15,
  34. FrameLabel -> {{"DEC", None}, {"RA", None}}];
  35.  
  36. sky = GeoGraphics[{Black, Point@GeoPosition@data},
  37. GeoRange -> {All, {0, 360}}, PlotRangePadding -> Scaled@.01,
  38. GeoGridLinesStyle -> Directive[Black, Dashed],
  39. GeoProjection -> "Mollweide", GeoGridLines -> Automatic,
  40. GeoBackground -> White, Frame -> True,
  41. ImagePadding -> {{60, 40}, {40, 15}}, ImageSize -> 800,
  42. 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}]},
  43. Background -> White, FrameStyle -> White, TicksStyle -> 15,
  44. FrameLabel -> {{"DEC", None}, {"RA", None}}];
  45.  
  46. sky3 = Show[Graphics @@ frame, Graphics @@ sky, PlotRange -> All]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement