Skip to content

Commit

Permalink
Merge pull request sunpy#7352 from nabobalis/example
Browse files Browse the repository at this point in the history
Add text from an email from Albert - Lazy edition
  • Loading branch information
nabobalis authored Jan 5, 2024
2 parents b4a7b2e + e1302c9 commit d524b1e
Show file tree
Hide file tree
Showing 2 changed files with 10 additions and 3 deletions.
1 change: 1 addition & 0 deletions changelog/7352.doc.rst
Original file line number Diff line number Diff line change
@@ -0,0 +1 @@
Added explanation text to :ref:`sphx_glr_generated_gallery_plotting_plotting_blank_map.py` about the offset between "(0, 0)" in helioprojective coordinates and the heliographic equator.
12 changes: 9 additions & 3 deletions examples/plotting/plotting_blank_map.py
Original file line number Diff line number Diff line change
Expand Up @@ -54,11 +54,17 @@
yc = [0, 100, 400] * u.arcsec

################################################################################
# Place and mark coordinates on the plot.
# Plot the blank map with the specified coordinates. Note that the marker for
# (0, 0) in helioprojective coordinates is at the center of the solar disk, yet
# the heliographic equator (zero degrees latitude) does not go through the disk
# center and instead curves below it. The reason for that is the observer,
# specified as Earth in this example, is almost always at non-zero heliographic
# latitude, and disk center as seen by such an observer will have that same
# heliographic latitude. The :func:`~sunpy.coordinates.sun.B0` function returns
# Earth's heliographic latitude at a specified time.

coords = SkyCoord(xc, yc, frame=blank_map.coordinate_frame)
p = ax.plot_coord(coords, 'o')
# Set title.
ax.plot_coord(coords, 'o')
ax.set_title('Plotting fixed points on a blank map')

plt.show()

0 comments on commit d524b1e

Please sign in to comment.