HMI limb on STEREO map

[jcmo]

Hi Guys,

I would lie to plot the HMI limb onto a stereo image, but I don't know how to do that in sunpy. Could you help me with this? 

Cheers,

Juan

[DVD PS]

Hi Juan,

I'm not sure if this is the best way, and I don't know whether you are asking just to overlie the contour on top, or you are asking something more complicated as knowing where the limb of the sun seen from Earth is on as from seen from other point of view (eg, stereo).  If the last is what you are asking, then I think we still need to do some development with the WCS transformations.

But if it's the first what you meant, then you could do as it follows (AIA with EIT limb):

import sunpy
import matplotlib.pyplot as plt
aia  = sunpy.make_map(sunpy.AIA_171_IMAGE)
fig = plt.figure()
ax = plt.subplot(111)
aia.plot()
eit  = sunpy.make_map(sunpy.EIT_195_IMAGE)
plt.colorbar()
eit.draw_limb()
plt.show()


Cheers,

David

--
You received this message because you are subscribed to the Google Groups "SunPy" group.
To unsubscribe from this group and stop receiving emails from it, send an email to sunpy+un...@googlegroups.com.
To post to this group, send email to su...@googlegroups.com.
Visit this group at http://groups.google.com/group/sunpy?hl=en.
For more options, visit https://groups.google.com/groups/opt_out.
 
 

[Christe, Steven D. (GSFC-6710)]

Hi David and Juan, I actually think all of the base functionality is in sunpy to be able to do this. Give me a little bit of time to look into this.

[Jack Ireland]

Steven,

You need to know the location of STEREO in space.  We don't have that.  The IDL code pb0r.pro

http://hesperia.gsfc.nasa.gov/ssw/gen/idl/solar/pb0r.pro

has a branch

http://hesperia.gsfc.nasa.gov/ssw/gen/idl/solar/pb0r_stereo.pro

for calculating solar p, b0 and radius from STEREO which calls all sorts of STEREO-related routines via

http://darts.isas.jaxa.jp/pub/ssw/stereo/gen/idl/spice/get_stereo_lonlat.pro

and others.  SPICE is down there ultimately.

SunPy (or at least my branch of it) has a version pb0r.pro for the Sun as seen from Earth only.

We need to include pySPICE, and then we can work on supporting STEREO and Solar Orbiter, Probe and SOHO.

jack

[Christe, Steven D. (GSFC-6710)]

Jack, as far as I know that information should be included in the (thankfully) WCS-compliant STEREO fits file.

[Jack Ireland]

Looks that way.  Someone needs to try it through.

And that someone might be me.  Helioviewer could really use some convenient functions taking feature/event locations taken from along the Sun-Earth line and projecting them onto the disk as seen from STEREO, or from another WCS header compliant mission.  I'll pencil that in for June.

[jcmo]

Hi Guys,

Thanks for all your ideas. The way that I am doing this with SSWIDL is based on a tutorial that William Thomson wrote. These are the commands to "re-project" the limb:

wcs_limb, wcs_hmi, lon, lat

wcs_convert_to_coord, wcs_euvi, coord, 'hg', lon, lat

pix = wcs_get_pixel(wcs_euvi, coord)

In principle, I could use the output of this routines and continue doing my analysis with Sunpy. The output is in pixels coordinates, so I am not sure how to overplot the results on a map.

Thanks again,

Juan

[Christe, Steven D. (GSFC-6710)]

Juan, the sunpy wcs code is based on bill's code so all of the functionality you need should already be there and you should be able to translate this IDL code to sunpy fairly straightforwardly. Please check out the sunpy wcs functionality that is available.

From: jcmo <jua...@gmail.com>
Reply-To: "su...@googlegroups.com" <su...@googlegroups.com>
Date: Sunday, April 28, 2013 0:31 AM
To: "su...@googlegroups.com" <su...@googlegroups.com>

--

[Erica Lastufka]

Hi all,

Came across this old thread via Google. I want to do something similar, overplotting the AIA limb on STEREO. Was there any SunPy solution to this HMI/STEREO situation? Right now I'm playing with the longitude line objects but haven't managed to apply a coordinate transform to them yet.

Cheers,
~ Erica

[Stuart Mumford]

Hello Erica,

Thanks for the question, this is now really easy to do with the 0.7 coordinate framework and Astropy 1.3.

See my example here: https://gist.github.com/Cadair/de415eeeaad7f340c33cde1a8f618898

Thanks,

Stuart

Visit this group at https://groups.google.com/group/sunpy.
For more options, visit https://groups.google.com/d/optout.

[Stuart Mumford]

Hello,

There was an error in my original version below, I have now updated it to be compatible with SunPy 0.7.

Stuart

[Erica Lastufka]

Hi Stuart,

Thanks for that example, it's beautifully simple!

Cheers,

~ Erica

To unsubscribe from this group and stop receiving emails from it, send an email to sunpy+unsubscribe@googlegroups.com.

To post to this group, send email to su...@googlegroups.com.
Visit this group at https://groups.google.com/group/sunpy.
For more options, visit https://groups.google.com/d/optout.

--
You received this message because you are subscribed to the Google Groups "SunPy" group.

To unsubscribe from this group and stop receiving emails from it, send an email to sunpy+unsubscribe@googlegroups.com.

To post to this group, send email to su...@googlegroups.com.
Visit this group at https://groups.google.com/group/sunpy.
For more options, visit https://groups.google.com/d/optout.

--
You received this message because you are subscribed to a topic in the Google Groups "SunPy" group.
To unsubscribe from this topic, visit https://groups.google.com/d/topic/sunpy/Bti-K5OsLBc/unsubscribe.
To unsubscribe from this group and all its topics, send an email to sunpy+unsubscribe@googlegroups.com.