| Re: Limits of custom axis [message #48251] |
Tue, 04 April 2006 09:06 |
burkina
Messages: 32
Registered: February 2005
|
Member |
|
|
Ok, I reconstructed the limits of the image directly from my script,
using the fits headers.
Everything seems ok now.
Thanks,
Stefano
|
|
|
|
| Re: Limits of custom axis [message #48278 is a reply to message #48251] |
Mon, 03 April 2006 08:33  |
David Fanning
Messages: 11724
Registered: August 2001
|
Senior Member |
|
|
burkina writes:
> This seems close to what I want, but not quite.
> The values for !X.CRange are not in "data" units, which is what I
> need.... I would like them to be in arcsec, like the imcontour axes.
Oh, right. IMCONTOUR doing things in its own way. I
guess you will just have to figure out how it works
and return the data from the variables it creates internally.
Cheers,
David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
|
|
|
|
| Re: Limits of custom axis [message #48279 is a reply to message #48278] |
Mon, 03 April 2006 08:14 |
burkina
Messages: 32
Registered: February 2005
|
Member |
|
|
David Fanning writes:
> I think you are looking for ![XYZ].CRANGE.
>
> axis, xaxis=1, xrange=[!X.CRange[0]*dist, !X.CRange[1]*dist], $
> xstyle=1, xtitle="pc, etc.
This seems close to what I want, but not quite.
The values for !X.CRange are not in "data" units, which is what I
need.... I would like them to be in arcsec, like the imcontour axes.
Thanks,
Stefano
|
|
|
|
| Re: Limits of custom axis [message #48283 is a reply to message #48279] |
Mon, 03 April 2006 06:48 |
David Fanning
Messages: 11724
Registered: August 2001
|
Senior Member |
|
|
burkina writes:
> I'm plotting an astronomical image with imcontour, suppressing two
> axes:
>
> imcontour, alog10(sm), hdr, position=pos, Levels=vals,
> color=fsc_color('white'), TYPE=0,/keep_aspect, /noerase, xstyle=8,
> ystyle=8, XMID=xmid, YMID=ymid
>
> I want those two axes to show a distance in pc, instead of arcsec.
> That's why I use:
>
> axis, xaxis=1, xrange=[9.5*dist, -16.5*dist], xstyle=1, xtitle="pc",
> /save,color=fsc_color('black')
>
> where dist is the distance in pc corresponding to 1 arcsec for that
> source.
> The problem is that I get the numbers 9.5 and -16.5 directly by eye. I
> wish I had a command which would return exactly the lower and upper
> limit of the axis created by imcontour.
> How can I do that?
I think you are looking for ![XYZ].CRANGE.
axis, xaxis=1, xrange=[!X.CRange[0]*dist, !X.CRange[1]*dist], $
xstyle=1, xtitle="pc, etc.
> I would also like to plot the contours in white, the ticks also in
> white, but the ticknumbers in black.
> All I can do now is to plot the contours in white and then call another
> time imcontour with /nodata and plot the rest in black. But I would
> like to draw the ticks also in white....
Be careful what you wish for: you could probably do this in
object graphics. :-)
Cheers,
David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.dfanning.com/
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
