[AstroPy] [SciPy-user] FITS images with header-supplied axes?
Bridgman, William T.
Wed Apr 2 11:52:42 CDT 2008
Here's my 'work python' list.
Mostly codes for organizing image and volumetric datasets.
I use PIL, PyFITS, scipy, numpy, VTK for 3-d data extraction and
image generation. I have some cases where I need to retrieve the
data in IDL to get a format that can be used on the python side.
I have a number of custom modules which co-register the images and 3-
d data into various geospatial coordinate systems. I also have
modules that extract spacecraft orbit information from sites like the
SSCweb (though the limitations in ZSI interfere with me migrated to
their newer web interface). The co-registration works with time-
dependent data and connects the data components at appropriate times
(however, it does not interpolate between actual datasets).
I use CGkit to map the images and objects created above onto 3-d
objects and generate files which are rendered into images by the
RenderMan computer graphics package.
On Apr 1, 2008, at 1:51 PM, Bridgman, William T. wrote:
> I agree there is some lack of coordination on some of these
> projects. Such is the nature of Open Source.
> I think there is also an issue where some of the projects are
> overkill for what individuals might need for production-pipeline/
> research/educational uses. Hence many of us start from scratch to
> keep the work compact. It's only later that we discover others
> working in similar directions.
> I obtained a copy of pywcs many months ago, trying to get
> heliographic coordinate systems installed. Other priorities
> intervened and by the time I got back to it, migration to numpy and
> changes in pyFITS seemed to break the previous version. I didn't get
> a chance to make any revisions operational.
> Would there be any interest in members of the list publishing a short
> description of what types of modules they are designing in their own
> work? It might be worthwhile for coordination & possible
> collaborations. My requirements for work projects are quite
> different from my recreational & educational python projects.
> On Apr 1, 2008, at 10:41 AM, email@example.com wrote:
>> I just want to discuss about a problem that our python astronomer
>> have and the
>> answer from Perry is a very good example of it. There are too many
>> which are doing exactly the same thing and essentially because
>> nothing is
>> centralized or worst, in the case of pywcs, not advertised.
>> For example Perry told us about pywcs developped by STSCI (it was
>> the first
>> time I saw any reference to this project) but Adam spoke about
>> astlib ( http://astlib.sourceforge.net/ ) which are different
>> package with
>> exactly the same goal and I did myself something similar (even if
>> it was fast
>> and dirtier).
>> I think it's time to try to identified the most important task that
>> astronomer are needing and try to centralized all the effort at the
>> place. The astropy mail list is probably a good start as the
>> astropy.scipy.org website.
>> Perhaps we can start to identified the need and desiderata. We will
>> have too
>> many but that some most important or urgent must be identified. As
>> we clearly needed the pywcs and thanks to STSCI we have it now. We
>> also need
>> a package to plot our data, images etc. Matplotlib is very good but
>> does have
>> a major problem (at least for my point of view), it's slow, very
>> slow for big
>> array or can't even produce the image if the image is too big. ( If I
>> remember a precedent discussion, the problem is mainly due to Agg.)
>> After we need to know, if possible, for which project if someone is
>> something, who and how to contact him/her. So interested people can
>> eventually help him/her to do it. It's seems that STSCI is doing
>> most of the
>> work but I'm pretty sure that other people are willingly to help
>> them to
>> extend python to be the ideal tool for astronomer.
>> Just my 2 cents,
>> Le Sunday 30 March 2008 10:54:04 Perry Greenfield, vous avez écrit :
>>> On Mar 28, 2008, at 8:28 PM, Keflavich wrote:
>>>> Is there any plotting routine in scipy / matplotlib that can plot a
>>>> fits image with correct WCS coordinates on the axes? I know pyfits
>>>> can load fits files, astLib has routines to interpret header
>>> coordinates, and I think you can make the axes different using
>>>> matplotlib transforms, but is there anything that puts all three
>>>> together currently available?
>>> Well, we (STScI) recently wrapped WCSLIB to obtain a mapping
>>> between pixel and sky coordinates for python (you can find it as
>>> pywcs in
>>> astrolib on scipy; that may have been what you were referring to).
>>> But I'm not sure you understand what you are asking for with
>>> regard to
>>> matplotlib. The new transforms stuff should make it much easier to
>>> the sky coordinates in the interactive display. The axis labeling
>>> is a
>>> different matter. Suppose your image (let's say it's 1Kx1K for the
>>> sake of
>>> discussion) is rotated 45 degrees with regard to north (either
>>> way, it
>>> doesn't really matter). What would you expect to see for axis
>>> labels? I
>>> don't think it is at all obvious how people would want labeling to
>>> be done
>>> along the edges of the image. I can imagine someone wanting axes
>>> or grids
>>> superimposed on the image itself, but that's not quite the same
>>> thing. Do
>>> you want the image rotated so that it is resampled on to RA and
>>> Dec and
>>> displayed that way?
>>> In any event, no we haven't yet done anything to try to integrate
>>> all three
>>> things. Among other things we wanted to make sure that the api for
>>> the wcs
>>> info was suitable before doing a lot with it (and in the meantime,
>>> Mike is
>>> working on rewriting drizzle which is taking a lot of his time).
>>> SciPy-user mailing list
>> AstroPy mailing list
> **** Please Update my e-mail address to my *NEW* OneNASA address
> below ****
> Dr. William T."Tom" Bridgman Scientific Visualization
> Global Science & Technology, Inc. NASA/Goddard Space Flight
> Email: William.T.Bridgman@nasa.gov Code 610.3
> Phone: 301-286-1346 Greenbelt, MD 20771
> FAX: 301-286-1634 http://svs.gsfc.nasa.gov/
> AstroPy mailing list
**** Please Update my e-mail address to my *NEW* OneNASA address
Dr. William T."Tom" Bridgman Scientific Visualization
Global Science & Technology, Inc. NASA/Goddard Space Flight
Email: William.T.Bridgman@nasa.gov Code 610.3
Phone: 301-286-1346 Greenbelt, MD 20771
FAX: 301-286-1634 http://svs.gsfc.nasa.gov/
More information about the AstroPy