On 13/02/12 23:50, Elaina Ann Hyde wrote:
Hi, I'm working on a routine for reading in the Schlegel dustmaps. I
have an ascii table with values, Ra-Dec and I'm trying to convert 2
columns to l,b and get the dust values out, this is a 2 part problem as
I need to first convert to l,b, keep those values for other uses
(calculating VLSR), and get extinctions for each l,b value,
That sounds very interesting but almost totally meaningless to most of
us on this list I suspect! :-)
import astropysics
from astropysics.coords import ICRSCoordinates,GalacticCoordinates
import asciitable
From those lines alone I suspect you might have more luck on either a
specialised forum or on the general Python mailing list/newsgroup.
It looks like a much more specialised query than is usual on this list
which is for beginning Python programmers.
However....
dustmap='SFD_dust_4096_ngp.fits'
EB_V=astropysics.obstools.get_SFD_dust(l,b,dustmap,interpolate=True)
-------
this gives the error 'Only length-1 arrays can be converted to Python
scalars' ... however, I cannot do:
Always post the full error text not just the last line or a summary.
It usually contains a lot of useful information, especially what code
Python actually saw when it tried to execute your code.
gcoords=ICRSCoordinates(dat['ra-drad'],dat['dec-drad'],radians=True).convert(GalacticCoordinates)
without the same error. Any ideas would be appreciated, thanks!
On the surface I have no idea because that looks like an issue
related to the use of the astrophysics module. But the error
message may reveal a more fundamental issue underneath.
--
Alan G
Author of the Learn to Program web site
http://www.alan-g.me.uk/
_______________________________________________
Tutor maillist - Tutor@python.org
To unsubscribe or change subscription options:
http://mail.python.org/mailman/listinfo/tutor
