Hi,<div><br></div><div>I uploaded a script here</div><div><br></div><div><a href="https://gist.github.com/4078531">https://gist.github.com/4078531</a><br></div><div><br></div><div>I made it a gist, since it sounds like it may need to be refined.</div>
<div><br></div><div>cheers,</div><div>chris</div><div class="gmail_extra"><br><br><div class="gmail_quote">On Thu, Nov 15, 2012 at 3:29 AM, David Berry <span dir="ltr">&lt;<a href="mailto:d.berry@jach.hawaii.edu" target="_blank">d.berry@jach.hawaii.edu</a>&gt;</span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="im">On 14 November 2012 22:01, Demitri Muna &lt;<a href="mailto:demitri.muna@gmail.com">demitri.muna@gmail.com</a>&gt; wrote:<br>

&gt; Hi,<br>
&gt;<br>
&gt; We are in the process of testing the new astropy coordinates package. One test that I think would be useful is this (which Tom also mentioned):<br>
&gt;<br>
&gt; * generate 1000 ra/dec J2000 random points on the sky<br>
&gt; * convert them to, in steps,<br>
&gt;     - galactic<br>
&gt;     - B1950<br>
&gt;     - ecliptic<br>
&gt;     - equatorial J2000<br>
&gt;<br>
&gt; Can someone do this in IDL? The result should be<br>
&gt;<br>
&gt; * the file containing 1000 points<br>
&gt; * the resulting points at the end of the conversions<br>
&gt; * the IDL script<br>
&gt;<br>
&gt; The aim is to write a script using the coordinates package and see if there are any discrepancies. We will add this to the testing suite.<br>
&gt;<br>
&gt; I thought I&#39;d cc the astropy list as I&#39;m sure there are lots of people who are proficient in IDL and can probably do this pretty quickly.<br>
<br>
</div>1) Different conversion algorithms are always (or nearly always) going<br>
to produce different answers. So do you have a target accuracy for the<br>
astropy coordinates package? Do you want arcseconds accuracy or<br>
milli-arcsecond accuracy? The more accuracy you want, the more<br>
sophisticated the algorithms need to be.<br>
<br>
2) Why does it need to be IDL? It&#39;s proprietary software, so many<br>
people will not use it. There are many tools in many freely available<br>
languages for creating tables of corresponding positions in different<br>
coordinate systems.   For instance, I could create such a table<br>
quickly using PyAST, which since we&#39;re talking Python seems an obvious<br>
choice.<br>
<br>
3) You include ecliptic in the list. Does the coordinates package<br>
support ecliptic? I could not see it in the code.<br>
<br>
4) Your list of coordinate systems starts at &quot; ra/dec J2000&quot; and ends<br>
at &quot;equatorial J2000&quot; - are you making a distinction between these?<br>
And I presume you mean FK5 J2000 (there is a common usage of &quot;J2000<br>
RA/Dec&quot; which uses the mean dynamical equator and equinox of the J2000<br>
epoch as reference, which is different to FK5)?<br>
<br>
5) Was the omission of ICRS intentional?<br>
<br>
6) No mention of the epoch of observation. As I mentioned in another<br>
message you need to specify the epoch of observation, to be able to<br>
convert accurately to and from an FK4 RA/Dec system (which I presume<br>
is what you mean by &quot;B1950&quot;). Strictly, you also need the epoch of<br>
observation when converting between FK5 and ICRS, albeit the<br>
dependency on the epoch is much weaker.<br>
<br>
7) You need to specify an equinox for the ecliptic coordinates.<br>
<span class="HOEnZb"><font color="#888888"><br>
David<br>
</font></span><div class="HOEnZb"><div class="h5">_______________________________________________<br>
AstroPy mailing list<br>
<a href="mailto:AstroPy@scipy.org">AstroPy@scipy.org</a><br>
<a href="http://mail.scipy.org/mailman/listinfo/astropy" target="_blank">http://mail.scipy.org/mailman/listinfo/astropy</a><br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>************************************<br>Chris Beaumont<br>Graduate Student<br>Institute for Astronomy<br>University of Hawaii at Manoa<br>2680 Woodlawn Drive<br>
Honolulu, HI 96822<br><a href="http://www.ifa.hawaii.edu/~beaumont">www.ifa.hawaii.edu/~beaumont</a><br>************************************<br>
</div>