From tpursimo at not.iac.es  Mon Dec  2 09:26:43 2019
From: tpursimo at not.iac.es (Tapio Pursimo)
Date: Mon, 2 Dec 2019 14:26:43 +0000 (WET)
Subject: [AstroPy] get_moon ... and the phase
In-Reply-To: <ED8FBB09-1A0A-4936-8EBD-CC2DC104B0F4@swarthmore.edu>
References: <alpine.DEB.2.20.1911291233180.25596@annie>
 <ED8FBB09-1A0A-4936-8EBD-CC2DC104B0F4@swarthmore.edu>
Message-ID: <alpine.DEB.2.20.1912021421280.8479@annie>

Hello Eric,

Thanks!
The moon illumination is not really coordinate stuff, however for 
an optical astronomer it would be very handy to have the moon coordinates 
and illumination in one go.

Cheers,

   T a p i o


> Hi Tapio, 
> 
> Here?s how I handled this in some of my code: 
> 
> # Moon illumination formula from Meeus, ?Astronomical
> # Algorithms".  Formulae 46.1 and 46.2 in the 1991 edition,
> # using the approximation cos(psi) \approx -cos(i).  Error
> # should be no more than 0.0014 (p. 316).
> 
> moon_illum =  0.5 * (1. - sin(dec_sun)*sin(dec_moon) -
> 				 cos(dec_sun)*cos(dec_moon)*
> 				 cos(ra_sun - ra_moon))
> 
> That gives the percent illumination of the moon, between 0 and 1. 
> 
>  If you can get the RA and Dec of the sun and moon, and the above 
> precision is good enough for your purposes, this is a simple way to do 
> it. 
> 
> Eric
> 
> 
> 
> > On Nov 29, 2019, at 7:38 AM, Tapio Pursimo <tpursimo at not.iac.es> wrote:
> > 
> > 
> > Hello,
> > 
> > Can anybody help me how to get the moon phase using astropy?
> > I can get the coordinates using get_moon, but the phase is apparently 
> > missing.
> > I would not like to use/install "ephem" nor "astroplan".
> > 
> > Best regards,
> > 
> > Tapio Pursimo
> > 
> > _______________________________________________
> > AstroPy mailing list
> > AstroPy at python.org
> > https://mail.python.org/mailman/listinfo/astropy
> 
> 

From p3y1i4n at gmail.com  Mon Dec  2 12:06:45 2019
From: p3y1i4n at gmail.com (Pey Lian Lim)
Date: Mon, 2 Dec 2019 12:06:45 -0500
Subject: [AstroPy] AstroPy Digest, Vol 159, Issue 1
In-Reply-To: <mailman.43.1575306004.24015.astropy@python.org>
References: <mailman.43.1575306004.24015.astropy@python.org>
Message-ID: <CAJt49iTqDQcVBemZ5u8=oo1GY7zpZdnz_9e348tyxhtowF8o4w@mail.gmail.com>

Hello. Looks like it is already implemented in astroplan --
https://github.com/astropy/astroplan/blob/master/astroplan/moon.py

On Mon, Dec 2, 2019 at 12:01 PM <astropy-request at python.org> wrote:
>
> Today's Topics:
>
>    1. Re: get_moon ... and the phase (Tapio Pursimo)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Mon, 2 Dec 2019 14:26:43 +0000 (WET)
> From: Tapio Pursimo <tpursimo at not.iac.es>
> To: Astronomical Python mailing list <astropy at python.org>
> Subject: Re: [AstroPy] get_moon ... and the phase
> Message-ID: <alpine.DEB.2.20.1912021421280.8479 at annie>
> Content-Type: text/plain; charset="utf-8"
>
> Hello Eric,
>
> Thanks!
> The moon illumination is not really coordinate stuff, however for
> an optical astronomer it would be very handy to have the moon coordinates
> and illumination in one go.
>
> Cheers,
>
>    T a p i o

From shbhuk at gmail.com  Mon Dec  2 18:00:31 2019
From: shbhuk at gmail.com (Shubham)
Date: Mon, 2 Dec 2019 18:00:31 -0500
Subject: [AstroPy] Getting Lunar phase
In-Reply-To: <mailman.30.1575046803.9048.astropy@python.org>
References: <mailman.30.1575046803.9048.astropy@python.org>
Message-ID: <CAKL_h6dH7W0qfOpteOaEGyRYSk=ZPfCvfC9dG4BYZ+4=fEGAoA@mail.gmail.com>

Try this -

from astroplan import Observer


observatory = Observer(location=location, name="", timezone=timezone)

moon_illumination = observatory.moon_illumination(obstime)


*Where location is an Astropy EarthLocation Object - Astropy EarthLocation
object, and time is an Astropy Time obkect. *

*Shubham Kanodia*


On Fri, 29 Nov 2019 at 12:01, <astropy-request at python.org> wrote:

> Send AstroPy mailing list submissions to
>         astropy at python.org
>
> To subscribe or unsubscribe via the World Wide Web, visit
>         https://mail.python.org/mailman/listinfo/astropy
> or, via email, send a message with subject or body 'help' to
>         astropy-request at python.org
>
> You can reach the person managing the list at
>         astropy-owner at python.org
>
> When replying, please edit your Subject line so it is more specific
> than "Re: Contents of AstroPy digest..."
>
>
> Today's Topics:
>
>    1. get_moon ... and the phase (Tapio Pursimo)
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Fri, 29 Nov 2019 12:38:45 +0000 (WET)
> From: Tapio Pursimo <tpursimo at not.iac.es>
> To: astropy at python.org
> Subject: [AstroPy] get_moon ... and the phase
> Message-ID: <alpine.DEB.2.20.1911291233180.25596 at annie>
> Content-Type: text/plain; charset=US-ASCII
>
>
> Hello,
>
> Can anybody help me how to get the moon phase using astropy?
> I can get the coordinates using get_moon, but the phase is apparently
> missing.
> I would not like to use/install "ephem" nor "astroplan".
>
> Best regards,
>
>  Tapio Pursimo
>
>
>
> ------------------------------
>
> Subject: Digest Footer
>
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
>
>
> ------------------------------
>
> End of AstroPy Digest, Vol 158, Issue 7
> ***************************************
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191202/40b05b04/attachment.html>

From gabriel.rau at kit.edu  Tue Dec  3 03:22:50 2019
From: gabriel.rau at kit.edu (Rau, Gabriel (AGW))
Date: Tue, 3 Dec 2019 08:22:50 +0000
Subject: [AstroPy] Help needed: Amplitudes and phases from Lomb-Scargle
Message-ID: <85694d1eb05745d8924b53905a5e0b5e@kit-msx-29.kit.edu>

Hello,
Does anyone know how to recover harmonic component amplitudes and phases from the Astropy Lomb-Scargle implementation?
Many thanks,
Gabriel

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191203/9d097dec/attachment.html>

From bsipocz at gmail.com  Sat Dec  7 15:46:22 2019
From: bsipocz at gmail.com (Brigitta Sipocz)
Date: Sat, 7 Dec 2019 12:46:22 -0800
Subject: [AstroPy] astropy v4.0rc2 release candidate
Message-ID: <CAEjJ3fhzRxBX7y8v4LJ3M2-Hiq5jpZjtgFr5ebQTYefUU0coJg@mail.gmail.com>

Dear all,

The second release candidate for the new major astropy version has been
just released and available for testing.

https://pypi.python.org/pypi/astropy/4.0rc2

We would appreciate if you could try out this RC and report back any
success or failure via the following wiki page. If you maintain any
downstream packages please check whether they are compatible or there are
any blocking issues:

https://github.com/astropy/astropy/wiki/v4.0-RC-testing

Please also open issues for each different type of failure (please check
the issue tracker first to make sure no one has reported it earlier).

You should be able to install the RC by doing:

pip install astropy --pre

or

pip install astropy==4.0rc2

If no blocker issues come up, I plan to finalize the release next week
during the coordination meeting.

Cheers,
 Brigitta
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191207/411abb53/attachment.html>

From wtbsvs at gmail.com  Mon Dec  9 10:50:59 2019
From: wtbsvs at gmail.com (William Bridgman)
Date: Mon, 9 Dec 2019 10:50:59 -0500
Subject: [AstroPy] Using Astropy to plot skymaps...possible error in
 tutorial?
Message-ID: <CAKJeYxR6M+fOK3BWCTjVXEKhseh_f_KSxZSbiS_R3bEntuDJLA@mail.gmail.com>

I occasionally need to generate reference charts and plots for some of
my work producing content for general audiences and I was trying to
use the charting tutorial at the bottom of
http://learn.astropy.org/rst-tutorials/plot-catalog.html

However, if I plot the Yale Bright star catalog using this algorithm,
I get unrecognizable constellations.

ra,dec,magnitude = YaleStars.ExtractBasic()
coords_icrs = coord.SkyCoord(ra,dec,frame='icrs')
fig = plt.figure(figsize=(8,4.5))
axes = fig.add_subplot(111,projection='mollweide')
axes.scatter(ra.radian, dec.radian,s=(vMagThreshold - magnitude+1)**2)
axes.set_xticklabels(['14h','16h','18h','20h', '22h', '0h',
'2h','4h','6h','8h','10h'])
axes.grid(True)
plt.show()

Checking against my reference charts clearly shows the constellations
reversed (see ICRSflipped.pdf).

A change to the above fixed this issue (ICRS.pdf):

axes.scatter(-ra.radian, dec.radian,s=(vMagThreshold - magnitude+1)**2)
axes.set_xticklabels(['10h','8h','6h','4h', '2h', '0h',
'22h','20h','18h','16h','14h'])

Is there a legitimate convention for plotting skymaps reversed,
perhaps to match instrument view?

What's the procedure for revising tutorials?

Thanks,
Tom
-- 
Dr. William T."Tom" Bridgman
Scientific Visualization Studio
NASA/Goddard Space Flight Center
Greenbelt, MD  20771
http://svs.gsfc.nasa.gov
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ICRSflipped.pdf
Type: application/pdf
Size: 103673 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/astropy/attachments/20191209/4489fc37/attachment-0002.pdf>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: ICRS.pdf
Type: application/pdf
Size: 103657 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/astropy/attachments/20191209/4489fc37/attachment-0003.pdf>

From kellecruz at gmail.com  Mon Dec  9 12:25:31 2019
From: kellecruz at gmail.com (Kelle Cruz)
Date: Mon, 9 Dec 2019 12:25:31 -0500
Subject: [AstroPy] Using Astropy to plot skymaps...possible error in
 tutorial?
In-Reply-To: <CAKJeYxR6M+fOK3BWCTjVXEKhseh_f_KSxZSbiS_R3bEntuDJLA@mail.gmail.com>
References: <CAKJeYxR6M+fOK3BWCTjVXEKhseh_f_KSxZSbiS_R3bEntuDJLA@mail.gmail.com>
Message-ID: <CAApiRV7nL-o4+D+xam-VbaFzVmyk59BvvodYqnwWJ_Zc3E-=8g@mail.gmail.com>

Hi Tom,

Cool, glad to hear that tutorial is useful! And this tweak to the tutorial
would be very welcome! As far as I know, I don't think there is a
convention and I think pointing out that folks may want to think
deliberately about which way to do it is useful.

Please consider posting an issue here in the tutorials repo:
https://github.com/astropy/astropy-tutorials/issues

If you wanted to make the tweak yourself, that would also be amazing. You
could fork the repo and modify the notebook or download it and edit it in
Google Colab and then share it with us. Both methods are described in the
Contributor instructions:
https://github.com/astropy/astropy-tutorials/blob/master/CONTRIBUTING.rst#procedure-for-contributing

And notebook is here:
https://github.com/astropy/astropy-tutorials/tree/master/tutorials/notebooks/plot-catalog

Thanks!
Kelle

--
Kelle Cruz, PhD
917.837.9748 ? Hunter: x16486 ? AMNH: x7930
Pronouns: she/her and they/them



On Mon, Dec 9, 2019 at 12:15 PM William Bridgman <wtbsvs at gmail.com> wrote:

> I occasionally need to generate reference charts and plots for some of
> my work producing content for general audiences and I was trying to
> use the charting tutorial at the bottom of
> http://learn.astropy.org/rst-tutorials/plot-catalog.html
>
> However, if I plot the Yale Bright star catalog using this algorithm,
> I get unrecognizable constellations.
>
> ra,dec,magnitude = YaleStars.ExtractBasic()
> coords_icrs = coord.SkyCoord(ra,dec,frame='icrs')
> fig = plt.figure(figsize=(8,4.5))
> axes = fig.add_subplot(111,projection='mollweide')
> axes.scatter(ra.radian, dec.radian,s=(vMagThreshold - magnitude+1)**2)
> axes.set_xticklabels(['14h','16h','18h','20h', '22h', '0h',
> '2h','4h','6h','8h','10h'])
> axes.grid(True)
> plt.show()
>
> Checking against my reference charts clearly shows the constellations
> reversed (see ICRSflipped.pdf).
>
> A change to the above fixed this issue (ICRS.pdf):
>
> axes.scatter(-ra.radian, dec.radian,s=(vMagThreshold - magnitude+1)**2)
> axes.set_xticklabels(['10h','8h','6h','4h', '2h', '0h',
> '22h','20h','18h','16h','14h'])
>
> Is there a legitimate convention for plotting skymaps reversed,
> perhaps to match instrument view?
>
> What's the procedure for revising tutorials?
>
> Thanks,
> Tom
> --
> Dr. William T."Tom" Bridgman
> Scientific Visualization Studio
> NASA/Goddard Space Flight Center
> Greenbelt, MD  20771
> http://svs.gsfc.nasa.gov
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191209/60f3b61c/attachment.html>

From Paolo.Tanga at oca.eu  Mon Dec  9 16:55:20 2019
From: Paolo.Tanga at oca.eu (Paolo Tanga)
Date: Mon, 9 Dec 2019 22:55:20 +0100
Subject: [AstroPy] Using Astropy to plot skymaps...possible error in
 tutorial?
In-Reply-To: <CAApiRV7nL-o4+D+xam-VbaFzVmyk59BvvodYqnwWJ_Zc3E-=8g@mail.gmail.com>
References: <CAKJeYxR6M+fOK3BWCTjVXEKhseh_f_KSxZSbiS_R3bEntuDJLA@mail.gmail.com>
 <CAApiRV7nL-o4+D+xam-VbaFzVmyk59BvvodYqnwWJ_Zc3E-=8g@mail.gmail.com>
Message-ID: <eb17141d-eb99-47bc-9a9f-a5e2d7ea4d04@oca.eu>

An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191209/321ce07a/attachment.html>

From wtbsvs at gmail.com  Tue Dec 10 08:31:26 2019
From: wtbsvs at gmail.com (William Bridgman)
Date: Tue, 10 Dec 2019 08:31:26 -0500
Subject: [AstroPy] Using Astropy to plot skymaps...possible error in
 tutorial?
In-Reply-To: <CAApiRV7nL-o4+D+xam-VbaFzVmyk59BvvodYqnwWJ_Zc3E-=8g@mail.gmail.com>
References: <CAKJeYxR6M+fOK3BWCTjVXEKhseh_f_KSxZSbiS_R3bEntuDJLA@mail.gmail.com>
 <CAApiRV7nL-o4+D+xam-VbaFzVmyk59BvvodYqnwWJ_Zc3E-=8g@mail.gmail.com>
Message-ID: <CAKJeYxQW-urQ9GeUvr6ufb2NwGVFbzULKy_bWq_Wx1aq_gaSPg@mail.gmail.com>

It appears most of this stuff is done as Jupyter notebooks - which is
really inconsistent with my workflows which require more command-line
style tools.
It will take some time for me to convert this stuff to a notebook,
especially if it needs to be totally self-contained/standalone.
Thanks for the pointers.
Tom

On Mon, Dec 9, 2019 at 12:26 PM Kelle Cruz <kellecruz at gmail.com> wrote:
>
> Hi Tom,
>
> Cool, glad to hear that tutorial is useful! And this tweak to the tutorial would be very welcome! As far as I know, I don't think there is a convention and I think pointing out that folks may want to think deliberately about which way to do it is useful.
>
> Please consider posting an issue here in the tutorials repo:
> https://github.com/astropy/astropy-tutorials/issues
>
> If you wanted to make the tweak yourself, that would also be amazing. You could fork the repo and modify the notebook or download it and edit it in Google Colab and then share it with us. Both methods are described in the Contributor instructions: https://github.com/astropy/astropy-tutorials/blob/master/CONTRIBUTING.rst#procedure-for-contributing
>
> And notebook is here:
> https://github.com/astropy/astropy-tutorials/tree/master/tutorials/notebooks/plot-catalog
>
> Thanks!
> Kelle
>
> --
> Kelle Cruz, PhD
> 917.837.9748 ? Hunter: x16486 ? AMNH: x7930
> Pronouns: she/her and they/them
>
>
>
> On Mon, Dec 9, 2019 at 12:15 PM William Bridgman <wtbsvs at gmail.com> wrote:
>>
>> I occasionally need to generate reference charts and plots for some of
>> my work producing content for general audiences and I was trying to
>> use the charting tutorial at the bottom of
>> http://learn.astropy.org/rst-tutorials/plot-catalog.html
>>
>> However, if I plot the Yale Bright star catalog using this algorithm,
>> I get unrecognizable constellations.
>>
>> ra,dec,magnitude = YaleStars.ExtractBasic()
>> coords_icrs = coord.SkyCoord(ra,dec,frame='icrs')
>> fig = plt.figure(figsize=(8,4.5))
>> axes = fig.add_subplot(111,projection='mollweide')
>> axes.scatter(ra.radian, dec.radian,s=(vMagThreshold - magnitude+1)**2)
>> axes.set_xticklabels(['14h','16h','18h','20h', '22h', '0h',
>> '2h','4h','6h','8h','10h'])
>> axes.grid(True)
>> plt.show()
>>
>> Checking against my reference charts clearly shows the constellations
>> reversed (see ICRSflipped.pdf).
>>
>> A change to the above fixed this issue (ICRS.pdf):
>>
>> axes.scatter(-ra.radian, dec.radian,s=(vMagThreshold - magnitude+1)**2)
>> axes.set_xticklabels(['10h','8h','6h','4h', '2h', '0h',
>> '22h','20h','18h','16h','14h'])
>>
>> Is there a legitimate convention for plotting skymaps reversed,
>> perhaps to match instrument view?
>>
>> What's the procedure for revising tutorials?
>>
>> Thanks,
>> Tom
>> --
>> Dr. William T."Tom" Bridgman
>> Scientific Visualization Studio
>> NASA/Goddard Space Flight Center
>> Greenbelt, MD  20771
>> http://svs.gsfc.nasa.gov
>> _______________________________________________
>> AstroPy mailing list
>> AstroPy at python.org
>> https://mail.python.org/mailman/listinfo/astropy
>
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy



-- 
Dr. William T."Tom" Bridgman
Scientific Visualization Studio
NASA/Goddard Space Flight Center
Greenbelt, MD  20771
http://svs.gsfc.nasa.gov

From m.shemuni at gmail.com  Thu Dec 12 03:47:56 2019
From: m.shemuni at gmail.com (Mohammad Shameoni Niaei)
Date: Thu, 12 Dec 2019 11:47:56 +0300
Subject: [AstroPy] Alt-Az Splitting a Fisheye Image
Message-ID: <CAAXud8MbmCvReGHOfK6D_CZh0F1mLmPFWFpD_YR7oiP3U_M_XA@mail.gmail.com>

Hello the Astropy community

I have a project to detect cloud coverage from ground-based data (All-Sky
camera data).
I'm trying to train a classifier using SVM or CNN.

The problem:
I want to train the machine for parts of the sky. And for this, I want to
chop the image data into pieces. These pieces can be polygons(This part is
not the problemn) or Alt-Az(This is the problem) areas. And I want to use
the Numpy masks (The problem can be solved using a nested loop. But that
would be inefficient I think).
[image: part.png]
At-Az mask example(What I want to achieve.* Image created using gimp*)

I can create a full radius mask which I call pizza mask.
[image: pizza.png]
Full radius (Pizza) mask example

The plan is to have 3 pizza masks for different radii and apply logical
comparison to achieve At-Az mask. (I'm open to better ideas. With this
method I'll need to calculate 3 pizza masks and it takes time.)

Here comes the catch: I need to know what is the linear radius of the given
altitude angle.(Sounds like *someone* didn't attend to
trigonometry/spherical astronomy courses)
The fisheye projection looks like I need to use Sky2Pix_ZenithalEqualArea.
But how?

-- 
Mohammad SHAMEONI NIAEI
Astronomer / Software Specialist
Astrofizik Ara?t?rma ve Uygulama Merkezi (ATASAM),
Yakutiye, ERZURUM/T?RK?YE
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191212/ca53aacc/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: part.png
Type: image/png
Size: 141879 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/astropy/attachments/20191212/ca53aacc/attachment-0002.png>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: pizza.png
Type: image/png
Size: 129012 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/astropy/attachments/20191212/ca53aacc/attachment-0003.png>

From stefanv at berkeley.edu  Thu Dec 12 15:44:36 2019
From: stefanv at berkeley.edu (Stefan van der Walt)
Date: Thu, 12 Dec 2019 12:44:36 -0800
Subject: [AstroPy] Alt-Az Splitting a Fisheye Image
In-Reply-To: <CAAXud8MbmCvReGHOfK6D_CZh0F1mLmPFWFpD_YR7oiP3U_M_XA@mail.gmail.com>
References: <CAAXud8MbmCvReGHOfK6D_CZh0F1mLmPFWFpD_YR7oiP3U_M_XA@mail.gmail.com>
Message-ID: <20191212204436.nxbh2mvshgbhopeq@aurelius.localdomain>

Hi Mohammad,

On Thu, 12 Dec 2019 11:47:56 +0300, Mohammad Shameoni Niaei wrote:
> The problem:
> I want to train the machine for parts of the sky. And for this, I want to
> chop the image data into pieces. These pieces can be polygons(This part is
> not the problemn) or Alt-Az(This is the problem) areas. And I want to use
> the Numpy masks (The problem can be solved using a nested loop. But that
> would be inefficient I think).

You may be interested in the skimage draw module, that does some of the
mask generation:

https://scikit-image.org/docs/dev/api/skimage.draw.html#skimage.draw.polygon2mask

Although, from what you write above it sounds like you have challenges
in generating the polygon coordinates from coordinates; I hope someone
else here can help out with that part!

Best regards,
St?fan

From bruce.simonson at gmail.com  Sat Dec 21 14:12:42 2019
From: bruce.simonson at gmail.com (Bruce Simonson)
Date: Sat, 21 Dec 2019 11:12:42 -0800
Subject: [AstroPy] Is the solar position calculation
Message-ID: <CAM7j2eDSR=hx1d8YTRh+XKR4YivM3yeuwH3s81Y+yNW2afVQCg@mail.gmail.com>

Greetings,

I'm new to astropy, and the mailing list.  I'm not sure (yet) how to search
previous posts, or if there's an appropriate forum to use to post questions.

I'm wondering if a particular solar position calculator has been
implemented in astropy.  Here is the reference:

Reda, I.; Andreas, A. (2003). Solar Position Algorithm for Solar Radiation
Applications. 55 pp.; NREL Report No. TP-560-34302, Revised January 2008.
<http://www.nrel.gov/docs/fy08osti/34302.pdf>

Here is an implementation, from NREL, online:

https://midcdmz.nrel.gov/solpos/spa.html

I would like to compute solar positions, based on (most of) the parameters
available in this solar position calculator -- ideally from a python
module, especially astropy, if the calculator has been implemented.

Sorry for the post, if this has been discussed elsewhere.

Thanks,
-Bruce
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191221/b035dd53/attachment.html>

From kellecruz at gmail.com  Sat Dec 21 14:58:47 2019
From: kellecruz at gmail.com (Kelle Cruz)
Date: Sat, 21 Dec 2019 14:58:47 -0500
Subject: [AstroPy] Is the solar position calculation
In-Reply-To: <CAM7j2eDSR=hx1d8YTRh+XKR4YivM3yeuwH3s81Y+yNW2afVQCg@mail.gmail.com>
References: <CAM7j2eDSR=hx1d8YTRh+XKR4YivM3yeuwH3s81Y+yNW2afVQCg@mail.gmail.com>
Message-ID: <CAApiRV71+tFE=QZH0cDZgXUQfnwkT8gfS+hFK+gyfR0wbKtjfw@mail.gmail.com>

Maybe it?s in sunpy?
https://docs.sunpy.org/en/stable/


On Sat, Dec 21, 2019 at 2:14 PM Bruce Simonson <bruce.simonson at gmail.com>
wrote:

> Greetings,
>
> I'm new to astropy, and the mailing list.  I'm not sure (yet) how to
> search previous posts, or if there's an appropriate forum to use to post
> questions.
>
> I'm wondering if a particular solar position calculator has been
> implemented in astropy.  Here is the reference:
>
> Reda, I.; Andreas, A. (2003). Solar Position Algorithm for Solar Radiation
> Applications. 55 pp.; NREL Report No. TP-560-34302, Revised January 2008.
> <http://www.nrel.gov/docs/fy08osti/34302.pdf>
>
> Here is an implementation, from NREL, online:
>
> https://midcdmz.nrel.gov/solpos/spa.html
>
> I would like to compute solar positions, based on (most of) the parameters
> available in this solar position calculator -- ideally from a python
> module, especially astropy, if the calculator has been implemented.
>
> Sorry for the post, if this has been discussed elsewhere.
>
> Thanks,
> -Bruce
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
>
-- 
kelle (via iOS)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191221/a1745c28/attachment.html>

From juanlu001 at gmail.com  Sat Dec 21 16:13:20 2019
From: juanlu001 at gmail.com (Juan Luis Cano)
Date: Sat, 21 Dec 2019 22:13:20 +0100
Subject: [AstroPy] Is the solar position calculation
In-Reply-To: <CAApiRV71+tFE=QZH0cDZgXUQfnwkT8gfS+hFK+gyfR0wbKtjfw@mail.gmail.com>
References: <CAM7j2eDSR=hx1d8YTRh+XKR4YivM3yeuwH3s81Y+yNW2afVQCg@mail.gmail.com>
 <CAApiRV71+tFE=QZH0cDZgXUQfnwkT8gfS+hFK+gyfR0wbKtjfw@mail.gmail.com>
Message-ID: <CAFmivUXqbT-+7zqZRZJr1Pi2ZL+eZ5F9+7r1VHXu_Y2TGF3NUA@mail.gmail.com>

It seems that it's implemented in pvlib:

https://pvlib-python.readthedocs.io/en/stable/generated/pvlib.solarposition.get_solarposition.html

Best,

On Sat, Dec 21, 2019 at 8:59 PM Kelle Cruz <kellecruz at gmail.com> wrote:

> Maybe it?s in sunpy?
> https://docs.sunpy.org/en/stable/
>
>
> On Sat, Dec 21, 2019 at 2:14 PM Bruce Simonson <bruce.simonson at gmail.com>
> wrote:
>
>> Greetings,
>>
>> I'm new to astropy, and the mailing list.  I'm not sure (yet) how to
>> search previous posts, or if there's an appropriate forum to use to post
>> questions.
>>
>> I'm wondering if a particular solar position calculator has been
>> implemented in astropy.  Here is the reference:
>>
>> Reda, I.; Andreas, A. (2003). Solar Position Algorithm for Solar
>> Radiation Applications. 55 pp.; NREL Report No. TP-560-34302, Revised
>> January 2008. <http://www.nrel.gov/docs/fy08osti/34302.pdf>
>>
>> Here is an implementation, from NREL, online:
>>
>> https://midcdmz.nrel.gov/solpos/spa.html
>>
>> I would like to compute solar positions, based on (most of) the
>> parameters available in this solar position calculator -- ideally from a
>> python module, especially astropy, if the calculator has been implemented.
>>
>> Sorry for the post, if this has been discussed elsewhere.
>>
>> Thanks,
>> -Bruce
>> _______________________________________________
>> AstroPy mailing list
>> AstroPy at python.org
>> https://mail.python.org/mailman/listinfo/astropy
>>
> --
> kelle (via iOS)
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
>


-- 
Juan Luis Cano
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191221/4227bd81/attachment.html>

From adrianmpw at gmail.com  Sun Dec 22 09:24:43 2019
From: adrianmpw at gmail.com (Adrian Price-Whelan)
Date: Sun, 22 Dec 2019 09:24:43 -0500
Subject: [AstroPy] Is the solar position calculation
In-Reply-To: <CAM7j2eDSR=hx1d8YTRh+XKR4YivM3yeuwH3s81Y+yNW2afVQCg@mail.gmail.com>
References: <CAM7j2eDSR=hx1d8YTRh+XKR4YivM3yeuwH3s81Y+yNW2afVQCg@mail.gmail.com>
Message-ID: <CAEUL7mGu6CvJhGuBp24snVHWRBAAgK0=QMeG_NRRqAJGDR_PfA@mail.gmail.com>

Hi Bruce --

This functionality and most of these parameters are available within
Astropy, but they are spread between a few different classes. The relevant
pieces are:
http://docs.astropy.org/en/latest/api/astropy.coordinates.get_sun.html#astropy.coordinates.get_sun
http://docs.astropy.org/en/latest/api/astropy.coordinates.EarthLocation.html#astropy.coordinates.EarthLocation
http://docs.astropy.org/en/latest/api/astropy.coordinates.builtin_frames.AltAz.html#astropy.coordinates.builtin_frames.AltAz

Here is a small demo (as a Jupyter notebook) comparing the azimuth and
altitude values computed on a grid of times in Astropy to the NREL form
output:
https://gist.github.com/9e58e1ba422c1df9b8d1fe28d790071b

There are some small differences in altitude (possibly because of the
approximate refraction estimate we are using from SOFA/ERFA), but otherwise
they agree quite well.

Let me know if you need any help tracking down various parameters and
figuring out where to pass them in to Astropy!

best,
Adrian

On Sat, Dec 21, 2019 at 2:13 PM Bruce Simonson <bruce.simonson at gmail.com>
wrote:

> Greetings,
>
> I'm new to astropy, and the mailing list.  I'm not sure (yet) how to
> search previous posts, or if there's an appropriate forum to use to post
> questions.
>
> I'm wondering if a particular solar position calculator has been
> implemented in astropy.  Here is the reference:
>
> Reda, I.; Andreas, A. (2003). Solar Position Algorithm for Solar Radiation
> Applications. 55 pp.; NREL Report No. TP-560-34302, Revised January 2008.
> <http://www.nrel.gov/docs/fy08osti/34302.pdf>
>
> Here is an implementation, from NREL, online:
>
> https://midcdmz.nrel.gov/solpos/spa.html
>
> I would like to compute solar positions, based on (most of) the parameters
> available in this solar position calculator -- ideally from a python
> module, especially astropy, if the calculator has been implemented.
>
> Sorry for the post, if this has been discussed elsewhere.
>
> Thanks,
> -Bruce
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
>


-- 
Adrian M. Price-Whelan
Flatiron Institute, NYC
http://adrn.github.io
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191222/c87eac13/attachment-0001.html>

From gordon.b.dennis at gmail.com  Sat Dec 28 07:54:11 2019
From: gordon.b.dennis at gmail.com (Gordon Dennis)
Date: Sat, 28 Dec 2019 12:54:11 -0000
Subject: [AstroPy] Continuum subtraction
Message-ID: <000001d5bd7d$e787b550$b6971ff0$@gmail.com>

Good afternoon all

I am observing low-Z spiral galaxies through H-alpha filters as tracers of
SFR.  Are there any Python libraries available to perform continuum
subtraction on such images?

Gordon Dennis

 

+44(0)1934 838096 / +44(0)7774628110 / Skype: gordon.dennis

Sent from my Analytical Engine

Description: analyticalEngine

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191228/7f1a4aa1/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.jpg
Type: image/jpeg
Size: 3301 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/astropy/attachments/20191228/7f1a4aa1/attachment.jpg>

From teuben at astro.umd.edu  Sat Dec 28 14:25:35 2019
From: teuben at astro.umd.edu (Peter Teuben)
Date: Sat, 28 Dec 2019 14:25:35 -0500
Subject: [AstroPy] Continuum subtraction
In-Reply-To: <000001d5bd7d$e787b550$b6971ff0$@gmail.com>
References: <000001d5bd7d$e787b550$b6971ff0$@gmail.com>
Message-ID: <bf5f6f35-6f91-ef45-9676-01ab930c94a5@astro.umd.edu>

both pyspeckit and specutils have methods to do this. but look on 
ascl.net, where i saw some other candidates


 ?-peter

,On 12/28/19 7:54 AM, Gordon Dennis wrote:
>
> Good afternoon all
>
> I am observing low-Z spiral galaxies through H-alpha filters as 
> tracers of SFR.? Are there any Python libraries available to perform 
> continuum subtraction on such images?
>
> Gordon Dennis
>
> +44(0)1934 838096 / +44(0)7774628110 / Skype: gordon.dennis
>
> Sent from my Analytical Engine
>
> Description: analyticalEngine
>
>
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191228/58b4220c/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.jpg
Type: image/jpeg
Size: 3301 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/astropy/attachments/20191228/58b4220c/attachment.jpg>

From pdzwig at summaventures.com  Sun Dec 29 12:32:17 2019
From: pdzwig at summaventures.com (Peter Dzwig)
Date: Sun, 29 Dec 2019 17:32:17 +0000
Subject: [AstroPy] Continuum subtraction
In-Reply-To: <000001d5bd7d$e787b550$b6971ff0$@gmail.com>
References: <000001d5bd7d$e787b550$b6971ff0$@gmail.com>
Message-ID: <08f670c3-41f0-87c8-e34b-582cc3089351@summaventures.com>

Thanks for asking this.  I have been looking at something a bit similar
in the optical and have also run into issues.

IMHO the specutils tools are OK/fine until the data gets noisy. Looking
at medium redshift (z~1.0) QSO spectra in particular (even woth a
reasonable S/N).

I too would be interested in any comments that people have on the
relative merits of the packages and recommendations.

Peter Dzwig

PS I also think that the specutils demos could make more use of real
data in this area and suggest ways of handling it.

On 28/12/2019 12:54, Gordon Dennis wrote:
> Good afternoon all
> 
> I am observing low-Z spiral galaxies through H-alpha filters as tracers
> of SFR.? Are there any Python libraries available to perform continuum
> subtraction on such images?
> 
> Gordon Dennis
> 
> ?
> 
> +44(0)1934 838096 / +44(0)7774628110 / Skype: gordon.dennis
> 
> Sent from my Analytical Engine
> 
> Description: analyticalEngine
> 
> ?
> 
> 
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
> 

-- 

Dr. Peter Dzwig

From pdzwig at summaventures.com  Sun Dec 29 12:31:29 2019
From: pdzwig at summaventures.com (Peter Dzwig)
Date: Sun, 29 Dec 2019 17:31:29 +0000
Subject: [AstroPy] Continuum subtraction
In-Reply-To: <000001d5bd7d$e787b550$b6971ff0$@gmail.com>
References: <000001d5bd7d$e787b550$b6971ff0$@gmail.com>
Message-ID: <6b92cd25-5938-414a-379d-8375241214c6@summaventures.com>

Thanks for asking this.  I have been looking at something a bit similar
in the optical and have also run into issues.

IMHO the specutils tools are OK/fine until the data gets noisy. Looking
at medium redshift (z~1.0) QSO spectra in particular (even woth a
reasonable S/N).

I too would be interested in any comments that people have on the
relative merits of the packages and recommendations.

Peter Dzwig

PS I also think that the specutils demos could make more use of real
data in this area and suggest ways of handling it.

On 28/12/2019 12:54, Gordon Dennis wrote:
> Good afternoon all
> 
> I am observing low-Z spiral galaxies through H-alpha filters as tracers
> of SFR.? Are there any Python libraries available to perform continuum
> subtraction on such images?
> 
> Gordon Dennis
> 
> ?
> 
> +44(0)1934 838096 / +44(0)7774628110 / Skype: gordon.dennis
> 
> Sent from my Analytical Engine
> 
> Description: analyticalEngine
> 
> ?
> 
> 
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
> 

-- 

Dr. Peter Dzwig

From rfinn at siena.edu  Mon Dec 30 11:47:46 2019
From: rfinn at siena.edu (Finn, Rose)
Date: Mon, 30 Dec 2019 11:47:46 -0500
Subject: [AstroPy] Continuum subtraction
In-Reply-To: <000001d5bd7d$e787b550$b6971ff0$@gmail.com>
References: <000001d5bd7d$e787b550$b6971ff0$@gmail.com>
Message-ID: <CAF=o2jF6V2nSr4s_fThTka+4AZoryemy8wcEB7DN77xiHUehsQ@mail.gmail.com>

Hi Gordon,

I have some code to do continuum subtraction of H-alpha images
https://github.com/rfinn/halphagui/wiki  It's tailored to the projects I am
working on, but I could help you adapt it for your needs if you are
interested.

Take care,
Rose

On Sat, Dec 28, 2019 at 7:51 AM Gordon Dennis <gordon.b.dennis at gmail.com>
wrote:

> Good afternoon all
>
> I am observing low-Z spiral galaxies through H-alpha filters as tracers of
> SFR.  Are there any Python libraries available to perform continuum
> subtraction on such images?
>
> Gordon Dennis
>
>
>
> +44(0)1934 838096 / +44(0)7774628110 / Skype: gordon.dennis
>
> Sent from my Analytical Engine
>
> [image: Description: analyticalEngine]
>
>
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
>


-- 
Rose A. Finn, PhD
Department of Physics & Astronomy
Siena College
Loudonville, NY  12211
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/astropy/attachments/20191230/a0f19b72/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.jpg
Type: image/jpeg
Size: 3301 bytes
Desc: not available
URL: <http://mail.python.org/pipermail/astropy/attachments/20191230/a0f19b72/attachment.jpg>

From gordon.b.dennis at gmail.com  Tue Dec 31 16:00:20 2019
From: gordon.b.dennis at gmail.com (Gordon Dennis)
Date: Tue, 31 Dec 2019 21:00:20 -0000
Subject: [AstroPy] AstroPy Digest, Vol 159, Issue 11
In-Reply-To: <mailman.24.1577725203.31322.astropy@python.org>
References: <mailman.24.1577725203.31322.astropy@python.org>
Message-ID: <000901d5c01d$5081dc10$f1859430$@gmail.com>

Many thanks Peter and Rose.  I'll take a good look at specutils and also
Rose's GitHub entry.  
Happy New Year
Gordon

-----Original Message-----
From: AstroPy [mailto:astropy-bounces+gordon.b.dennis=gmail.com at python.org]
On Behalf Of astropy-request at python.org
Sent: Monday, December 30, 2019 17:00
To: astropy at python.org
Subject: AstroPy Digest, Vol 159, Issue 11

Send AstroPy mailing list submissions to
	astropy at python.org

To subscribe or unsubscribe via the World Wide Web, visit
	https://mail.python.org/mailman/listinfo/astropy
or, via email, send a message with subject or body 'help' to
	astropy-request at python.org

You can reach the person managing the list at
	astropy-owner at python.org

When replying, please edit your Subject line so it is more specific than
"Re: Contents of AstroPy digest..."


Today's Topics:

   1. Re: Continuum subtraction (Peter Dzwig)
   2. Re: Continuum subtraction (Peter Dzwig)
   3. Re: Continuum subtraction (Finn, Rose)


----------------------------------------------------------------------

Message: 1
Date: Sun, 29 Dec 2019 17:32:17 +0000
From: Peter Dzwig <pdzwig at summaventures.com>
To: Gordon Dennis <gordon.b.dennis at gmail.com>, astropy at python.org
Subject: Re: [AstroPy] Continuum subtraction
Message-ID: <08f670c3-41f0-87c8-e34b-582cc3089351 at summaventures.com>
Content-Type: text/plain; charset=windows-1252

Thanks for asking this.  I have been looking at something a bit similar in
the optical and have also run into issues.

IMHO the specutils tools are OK/fine until the data gets noisy. Looking at
medium redshift (z~1.0) QSO spectra in particular (even woth a reasonable
S/N).

I too would be interested in any comments that people have on the relative
merits of the packages and recommendations.

Peter Dzwig

PS I also think that the specutils demos could make more use of real data in
this area and suggest ways of handling it.

On 28/12/2019 12:54, Gordon Dennis wrote:
> Good afternoon all
> 
> I am observing low-Z spiral galaxies through H-alpha filters as 
> tracers of SFR.? Are there any Python libraries available to perform 
> continuum subtraction on such images?
> 
> Gordon Dennis
> 
> ?
> 
> +44(0)1934 838096 / +44(0)7774628110 / Skype: gordon.dennis
> 
> Sent from my Analytical Engine
> 
> Description: analyticalEngine
> 
> ?
> 
> 
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
> 

-- 

Dr. Peter Dzwig


------------------------------

Message: 2
Date: Sun, 29 Dec 2019 17:31:29 +0000
From: Peter Dzwig <pdzwig at summaventures.com>
To: Gordon Dennis <gordon.b.dennis at gmail.com>, astropy at python.org
Subject: Re: [AstroPy] Continuum subtraction
Message-ID: <6b92cd25-5938-414a-379d-8375241214c6 at summaventures.com>
Content-Type: text/plain; charset=windows-1252

Thanks for asking this.  I have been looking at something a bit similar in
the optical and have also run into issues.

IMHO the specutils tools are OK/fine until the data gets noisy. Looking at
medium redshift (z~1.0) QSO spectra in particular (even woth a reasonable
S/N).

I too would be interested in any comments that people have on the relative
merits of the packages and recommendations.

Peter Dzwig

PS I also think that the specutils demos could make more use of real data in
this area and suggest ways of handling it.

On 28/12/2019 12:54, Gordon Dennis wrote:
> Good afternoon all
> 
> I am observing low-Z spiral galaxies through H-alpha filters as 
> tracers of SFR.? Are there any Python libraries available to perform 
> continuum subtraction on such images?
> 
> Gordon Dennis
> 
> ?
> 
> +44(0)1934 838096 / +44(0)7774628110 / Skype: gordon.dennis
> 
> Sent from my Analytical Engine
> 
> Description: analyticalEngine
> 
> ?
> 
> 
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
> 

-- 

Dr. Peter Dzwig


------------------------------

Message: 3
Date: Mon, 30 Dec 2019 11:47:46 -0500
From: "Finn, Rose" <rfinn at siena.edu>
To: Astronomical Python mailing list <astropy at python.org>
Subject: Re: [AstroPy] Continuum subtraction
Message-ID:
	<CAF=o2jF6V2nSr4s_fThTka+4AZoryemy8wcEB7DN77xiHUehsQ at mail.gmail.com>
Content-Type: text/plain; charset="utf-8"

Hi Gordon,

I have some code to do continuum subtraction of H-alpha images
https://github.com/rfinn/halphagui/wiki  It's tailored to the projects I am
working on, but I could help you adapt it for your needs if you are
interested.

Take care,
Rose

On Sat, Dec 28, 2019 at 7:51 AM Gordon Dennis <gordon.b.dennis at gmail.com>
wrote:

> Good afternoon all
>
> I am observing low-Z spiral galaxies through H-alpha filters as 
> tracers of SFR.  Are there any Python libraries available to perform 
> continuum subtraction on such images?
>
> Gordon Dennis
>
>
>
> +44(0)1934 838096 / +44(0)7774628110 / Skype: gordon.dennis
>
> Sent from my Analytical Engine
>
> [image: Description: analyticalEngine]
>
>
> _______________________________________________
> AstroPy mailing list
> AstroPy at python.org
> https://mail.python.org/mailman/listinfo/astropy
>


--
Rose A. Finn, PhD
Department of Physics & Astronomy
Siena College
Loudonville, NY  12211
-------------- next part --------------
An HTML attachment was scrubbed...
URL:
<http://mail.python.org/pipermail/astropy/attachments/20191230/a0f19b72/atta
chment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.jpg
Type: image/jpeg
Size: 3301 bytes
Desc: not available
URL:
<http://mail.python.org/pipermail/astropy/attachments/20191230/a0f19b72/atta
chment-0001.jpg>

------------------------------

Subject: Digest Footer

_______________________________________________
AstroPy mailing list
AstroPy at python.org
https://mail.python.org/mailman/listinfo/astropy


------------------------------

End of AstroPy Digest, Vol 159, Issue 11
****************************************