<div dir="ltr"><div><div><div>The current STC standard contains the following language in Section 4.5.1.3 on this subject:<br><br>
<h4 style="mso-list:l0 level4 lfo1"><a name="_Toc181531823"></a><span style="mso-bookmark:_Toc181531823"><span style="mso-fareast-font-family:Arial;mso-bidi-font-family:Arial"><span style="mso-list:Ignore"><span style="font:7pt "Times New Roman"">
</span></span></span>Ellipse</span></h4>
<p class="MsoNormal">The Ellipse (2-dimensional) is similar to the Circle but
has, in addition, a minor radius and a position angle. Position angles are
measured following the definition in Section 4.4.1.2.5
and refer to the first axis. The definition of an ellipse in a Cartesian
coordinate system is unambiguous, but this is not the case for spherical
coordinates. In a spherical coordinate system the ellipse shall be defined as
the intersection of an elliptical cone with the unit sphere, where the axes and
position angle describe the geometry of the cone.</p>
<br><br></div>Which is the second option you mention.<br>Considering that we have agreed on that in the past, I would suggest that we stick with that.<br><br></div>Cheers,<br><br></div> - Arnold<br></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr">-------------------------------------------------------------------------------------------------------------<br>Arnold H. Rots Chandra X-ray Science Center<br>Smithsonian Astrophysical Observatory tel: +1 617 496 7701<br>60 Garden Street, MS 67 fax: +1 617 495 7356<br>Cambridge, MA 02138 <a href="mailto:arots@cfa.harvard.edu" target="_blank">arots@cfa.harvard.edu</a><br>USA <a href="http://hea-www.harvard.edu/~arots/" target="_blank">http://hea-www.harvard.edu/~arots/</a><br>--------------------------------------------------------------------------------------------------------------<br><br></div></div></div>
<br><div class="gmail_quote">On Mon, May 9, 2016 at 6:29 AM, Walter Landry <span dir="ltr"><<a href="mailto:wlandry@caltech.edu" target="_blank">wlandry@caltech.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Everyone,<br>
<br>
I noticed the question in DALI 1.1 about whether we should add an<br>
Ellipse type. We already sort of support ellipses in our TAP service<br>
because we have users that need it. So it would be nice to add it to<br>
the spec.<br>
<br>
However, I think we might have a problem with how to define it.<br>
pgSphere defines an ellipse as<br>
<br>
If the center of any spherical ellipse is the North Pole, the<br>
perpendicular projection into the x-y-plane gives an ellipse as in<br>
two-dimensional space.<br>
<br>
<a href="http://pgsphere.projects.pgfoundry.org/types.html#dt.sellipse" rel="noreferrer" target="_blank">http://pgsphere.projects.pgfoundry.org/types.html#dt.sellipse</a><br>
<br>
On the other hand, q3c and our tinyhtm-based service define an ellipse<br>
as an intersection of an elliptical cone with the unit sphere.<br>
Essentially, it is the difference between projecting an ellipse<br>
perpendicularly or radially.<br>
<br>
To distinguish between these two choices, we can look at how the two<br>
implementations answer various questions about ellipses.<br>
<br>
1) Does ellipse E CONTAIN point P?<br>
<br>
This is expressible in closed form for both formulations. If you<br>
are feeling adventurous, you can write the expressions in ADQL.<br>
<br>
2) Does an ellipse E1 INTERSECT or CONTAIN ellipse E2?<br>
<br>
pgSphere uses an iterative algorithm. It is not in closed form.<br>
In contrast, I think (but have not implemented) that q3c/tinyhtm can<br>
compute it by finding the null spaces of 3x3 matrices. This is, in<br>
principle, expressible in closed form. But for numerical<br>
stability, you probably want to use something like QR or SVD.<br>
<br>
3) Does a polygon intersect/overlap with ellipse E?<br>
<br>
No one implements this.<br>
<br>
Just to make things more complicated, we also use PostGIS, which does<br>
not have ellipses at all.<br>
<br>
Cheers,<br>
Walter Landry<br>
</blockquote></div><br></div>