<div dir="ltr">Walter wanted polar coordinates, by which I presume he wants zero radial distance at the centre. Something perhaps like the linked screen shots show below - basic polar coords on the left and a (slightly) logarithmic radial coord on the right. These were produced using AST, which can do this sort of thing.<div>
<br></div><div>David<br><div class="gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:#f5f5f5;padding:5px;color:#222;font-family:arial;font-style:normal;font-weight:bold;font-size:13px;border:1px solid #ddd">
<a href="https://docs.google.com/file/d/0B4jEb03bhCeOZi15d1JHQmxhemc/edit?usp=drive_web" target="_blank" style="display:inline-block;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0px;border:none;width:100%"><img style="vertical-align: bottom; border: none;" src="https://ssl.gstatic.com/docs/doclist/images/icon_11_image_list.png"> <span dir="ltr" style="color:#15c;text-decoration:none;vertical-align:bottom">logpolar.png</span></a></div>
 </div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 12 December 2013 15:26, Tom McGlynn <span dir="ltr"><<a href="mailto:Thomas.A.McGlynn@nasa.gov" target="_blank">Thomas.A.McGlynn@nasa.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">I guess I'm a little unclear on the actual requirement here.  Are we<br>
trying to represent the plane of some theoretical disk, or are we<br>
viewing a disk from some other location and looking at its emission?<br>
Walter's original message suggested the second since the radius had<br>
units of degrees from the center of the disk.  If we wanted the first<br>
presumably the units of radius should be some physical dimension.<br>
<br>
The former requires some tweaks to represent in standard WCS<br>
projections since they implicitly assume that we are projecting from<br>
the sphere to a plane, whereas this would be a plane->plane<br>
projection.  However we can just say that we are far away from the<br>
disk such that it subtends a small angle in the sky and is<br>
perpendicular to our line of sight.  Since the maximum angular radius<br>
is small (as small as we like) we can use any of the azimuthal<br>
projections using the center of the disk as our reference pixel.<br>
There will be a linear transformation of the angular radius to the<br>
physical radius.<br>
<br>
If we are modelling the emission of the disk as seen from some finite<br>
radius, then again we can set the CRVALn to the center of the disk and<br>
everything falls out pretty straightfowardly.  Now the ARC projection<br>
is the one that linearly renders the angular radius (by definition)<br>
but other projections could be used as appropriate.  [If we are not on<br>
the axis perpendicular to the center of the disk, then it may be more<br>
convenient to choose the closest point of the disk as the coordinate<br>
origin rather than the center of the disk.]<br>
<br>
In the first case we'd have NAXISn x | CDELTn | << 1 whereas in the<br>
second case we'd have NAXISn x | CDELT n | ~ 1 since we presume the<br>
disk is subtending a significant fraction of the sky.  The CRPIX's<br>
would be NAXISn/2 (assuming we are trying to cover the disk) and<br>
CRVALn=(0,90).  I don't know if LON/LAT--ARC is legal and if not I<br>
guess one has to use RA/DEC[-]--ARC.<br>
<br>
        Regards,<br>
        Tom<br>
<div class="HOEnZb"><div class="h5"><br>
David Berry wrote:<br>
> On 12 December 2013 13:48, Phil Hodge <<a href="mailto:hodge@stsci.edu">hodge@stsci.edu</a>> wrote:<br>
>> Walter,<br>
>><br>
>> Aside from the log, isn't this the ARC projection with CRVAL1 and CRVAL2<br>
>> set to the right ascension and declination at the center of the<br>
>> (theoretical) disk?<br>
><br>
> I'm not sure how that would work. Can you give an example? Typing the following:<br>
><br>
> NAXIS1  = 1000<br>
> NAXIS2  = 1000<br>
> CTYPE1  = 'RA---ARC'<br>
> CTYPE2  = 'DEC--ARC'<br>
> CRPIX1  = 500<br>
> CRPIX2  = 500<br>
> CDELT1  = -0.001<br>
> CDELT2  = 0.001<br>
> CRVAL1  = 45.0<br>
> CRVAL2  = 45.0<br>
><br>
> into the "Headers" box at<br>
> <a href="http://starlink.jach.hawaii.edu/cgi-bin/ast/fits-plotter" target="_blank">http://starlink.jach.hawaii.edu/cgi-bin/ast/fits-plotter</a>, putting<br>
> "grid=1" into the "Plot settings" box, and pressing the "Submit"<br>
> button just produces a roughly recti-linear grid as you might expect.<br>
><br>
> I'm not sure how using an ARC projection gets round the issue that<br>
> FITS-WCS can only describe latitude, not co-latitude.<br>
><br>
> David<br>
><br>
</div></div><div class="HOEnZb"><div class="h5">> _______________________________________________<br>
> fitsbits mailing list<br>
> <a href="mailto:fitsbits@listmgr.cv.nrao.edu">fitsbits@listmgr.cv.nrao.edu</a><br>
> <a href="http://listmgr.cv.nrao.edu/mailman/listinfo/fitsbits" target="_blank">http://listmgr.cv.nrao.edu/mailman/listinfo/fitsbits</a><br>
><br>
<br>
_______________________________________________<br>
fitsbits mailing list<br>
<a href="mailto:fitsbits@listmgr.cv.nrao.edu">fitsbits@listmgr.cv.nrao.edu</a><br>
<a href="http://listmgr.cv.nrao.edu/mailman/listinfo/fitsbits" target="_blank">http://listmgr.cv.nrao.edu/mailman/listinfo/fitsbits</a><br>
</div></div></blockquote></div><br></div>