CINXE.COM

<!DOCTYPE html> <html lang="en"> <head> <title>iawn: II: Orbit Fitting-Least Squares Technique for Orbit Estimation</title>   <script async src="https://www.googletagmanager.com/gtag/js?id=G-DQPHPE6HZ3"></script> <script> window.dataLayer = window.dataLayer || []; function gtag(){dataLayer.push(arguments);} gtag('js', new Date()); gtag('config', 'G-DQPHPE6HZ3'); </script>  <meta charset="UTF-8"> <meta name="viewport" content="initial-scale = 1.0, maximum-scale = 1.0, user-scalable = no, width = device-width" />  <link rel="stylesheet" href="/css-n-inc/style.css" media="screen" />  <link rel="stylesheet" href="/css-n-inc/style.responsive.css" media="all" /> <link rel="stylesheet" type="text/css" href="http://fonts.googleapis.com/css?family=Merriweather&subset=latin" /> <script src="/css-n-inc/jquery.js"> </script> <script src="/css-n-inc/script.js"> </script> <script src="/css-n-inc/script.responsive.js"> </script> <script src="/css-n-inc/sorttable.js"> </script> </head> <body> <div id="art-main"> <div class="art-sheet clearfix"> <header class="art-header"> <div class="art-shapes"> <div class="art-object1140172892"></div> </div> </header>    <div class="art-layout-wrapper"> <div class="art-content-layout"> <div class="art-content-layout-row"> <div class="art-layout-cell art-sidebar1"> <div class="art-vmenublock clearfix"> <div class="art-vmenublockcontent"> <ul class="art-vmenu"> <li><a href="/index.shtml">Home</a></li> <li><a href="/about.shtml">ABOUT</a> <ul> <li><a href="/about/charter.shtml">Charter</a></li> <li><a href="/about/members.shtml">Members</a></li> <li><a href="/about/assets.shtml">Assets</a></li> </ul></li> <li><a href="/features.shtml" class="active">FEATURES</a> <ul class="active"> <li><a href="/features/SquirrelValley.shtml">Squirrel Valley Obs</a></li> <li><a href="/features/MAP.shtml">MAP NEO Disc. Project</a></li> <li><a href="/features/orbfitting.shtml" class="active">Least Squares Technique</a></li> <li><a href="/features/orbits.shtml">Intro. to Ellip. Orbits</a></li> <li><a href="/features/million.shtml">1 Million Orbits</a></li> <li><a href="/features/units.shtml">New Units for Impact Energies?</a></li> <li><a href="/features/comets.shtml">Comets</a></li> <li><a href="/features/phas.shtml">What is a PHA?</a></li> <li><a href="/features/impacts.shtml">Meteorite Impacts & Recovery</a></li> <li><a href="/features/welcome.shtml">Welcome New Members</a></li> <li><a href="/features/whatrneos.shtml">What are NEOs?</a></li> <li><a href="/features/astrometric.shtml">Astrometric Follow-up</a></li> <li><a href="/features/case-study.shtml">Case Study</a></li> <li><a href="/features/record.shtml">Record</a></li> <li><a href="/features/tenagra.shtml">Tenagra</a></li> </ul></li> <li><a href="/neo-news.shtml">NEO News</a></li> <li><a href="/meetings.shtml">MEETINGS</a> <ul> <li><a href="/meetings/19-steering-cmte.shtml">19th Steering Cmte</a></li> <li><a href="/meetings/18-steering-cmte.shtml">18th Steering Cmte</a></li> <li><a href="/meetings/17-steering-cmte.shtml">17th Steering Cmte</a></li> <li><a href="/meetings/16-steering-cmte.shtml">16th Steering Cmte</a></li> <li><a href="/meetings/15-steering-cmte.shtml">15th Steering Cmte</a></li> <li><a href="/meetings/14-steering-cmte.shtml">14th Steering Cmte</a></li> <li><a href="/meetings/13-steering-cmte.shtml">13th Steering Cmte</a></li> <li><a href="/meetings/12-steering-cmte.shtml">12th Steering Cmte</a></li> <li><a href="/meetings/11-steering-cmte.shtml">11th Steering Cmte</a></li> <li><a href="/meetings/10-steering-cmte.shtml">10th Steering Cmte</a></li> <li><a href="/meetings/09-steering-cmte.shtml">9th Steering Cmte</a></li> <li><a href="/meetings/08-steering-cmte.shtml">8th Steering Cmte</a></li> <li><a href="/meetings/07-steering-cmte.shtml">7th Steering Cmte</a></li> <li><a href="/meetings/06-steering-cmte.shtml">6th Steering Cmte</a></li> <li><a href="/meetings/05-steering-cmte.shtml">5th Steering Cmte</a></li> <li><a href="/meetings/04-steering-cmte.shtml">4th Steering Cmte</a></li> <li><a href="/meetings/smpag.shtml">SMPAG</a></li> <li><a href="/meetings/03-steering-cmte.shtml">3rd Steering Cmte</a></li> <li><a href="/meetings/02-steering-cmte.shtml">2nd Steering Cmte</a></li> <li><a href="/meetings/communications.shtml">Communications</a></li> <li><a href="/meetings/01-steering-cmte.shtml">1st Steering Cmte</a></li> </ul></li> <li><a href="/close-approaches.shtml">CLOSE APPROACHES</a> <ul> <li><a href="/close-approaches/ca-table.shtml">CA-table</a></li> <li><a href="/close-approaches/ca-table_multi-opp.shtml">Multi-Opp. CAs</a></li>  </ul></li> <li><a href="/campaigns.shtml">CAMPAIGNS</a> <ul> <li><a href="/obscamp/2023DZ2/" onclick="this.target='_blank'">2023 DZ2 Rapid response</a></li> <li><a href="/obscamp/2005LW3/" onclick="this.target='_blank'">2005 LW3 Timing Assessment</a></li> <li><a href="/obscamp/2019XS/" onclick="this.target='_blank'">2019 XS Timing Assessment</a></li> <li><a href="/obscamp/Apophis/" onclick="this.target='_blank'">Apophis</a></li> <li><a href="/obscamp/1999KW4/" onclick="this.target='_blank'">1999 KW4</a></li> <li><a href="http://2012tc4.astro.umd.edu/" onclick="this.target='_blank'">2012 TC4</a></li> </ul></li> <li><a href="/sup-docs.shtml">SUPPORT DOCUMENTS</a></li> <li><a href="/gallery.shtml">GALLERY</a></li> <li><a href="/misc/faqs.shtml">FAQs</a></li> <li><a href="/misc/links.shtml">Links</a></li> <li><a href="/misc/contacts.shtml">Contacts</a></li> </ul> </div> </div> </div> <div class="art-layout-cell art-content"> <article class="art-post art-article"> <h1 class="art-postheader">Understanding Orbital Uncertainties in the Context of Impact Predictions</h1> <div class="art-postcontent art-postcontent-0 clearfix">  <h3>II: Orbit Fitting-Least Squares Technique for Orbit Estimation</h3> <p>This is the second in a series of articles detailing the concept of orbital uncertainties as it relates to impact predictions. Celestial mechanics can be a complicated subject, but here we'll try to distill some of these difficulties and render the subject a bit more transparent.</p> <div class="image-caption-wrapper" style="width: 40%; float: left"><img src="/gallery/article_pics/290px-Orbit1.svg.png" style="width: 100%; max-width: 450 px;" alt="??" class="art-lightbox" /> <p class="caption">Orbital elements <br />Credit: (<a href="https://commons.wikimedia.org/wiki/File:Orbit1.svg">Lasunncty at the English Wikipedia</a>, <a href="http://creativecommons.org/licenses/by-sa/3.0/">CC BY-SA 3.0</a>, via Wikimedia Commons</p> </div> <p>Our previous article described some of the quantities required for determining an orbit, and thus for predicting where objects will be at any point in the future. We discussed the need for 6 orbital elements and presented 3: eccentricity, inclination, and semi-major axis. The remaining 3 involve the orientation of the ellipse and the time when the object is at its closest point to the sun. </p> <p>In practice we never really measure these orbital elements-rather we measure positions in space and time, and when we measure enough of these we can compute the full orbit. The actual mathematics is a little complicated and really beyond what we need to discuss here, but the point is that solutions exist for these problems and the calculations are routine on modern computers. All of these calculations rely on the concept of least squares, a mathematical technique that computes the best fit to a set of points. The least squares method was first published by Legendre in 1805, but Gauss had used a form of it for the orbital solution for the first minor planet Ceres , found in 1801.</p> <div class="image-caption-wrapper" style="width: 40%; float: right"><img src="/gallery/article_pics/988px-Linear_least_squares.svg.png" style="width: 100%; max-width: 450 px;" alt="??" class="art-lightbox" /> <p class="caption">Illustration of linear least squares. <br />Credit: (<a href="https://commons.wikimedia.org/wiki/File:Linear_least_squares.svg">Oleg Alexandrov at the English Wikipedia</a> via Wikimedia Commons</p> </div> <p>The least squares method works best when we have sets of constraints on our mathematical problem. For example, if you know you have an equation that represents a line in an X-Y graph, then you can assume that particular form for the solution.</p> <p>For orbital elements, our assumption is that the orbits are elliptical and this reduces the complexity of the least squares solution dramatically. Note that even if the orbit is hyperbolic, this is still a conic section and the math is quite similar to the elliptical calculations and one can use the same methods. This means it was just as easy to compute the orbit for 2I/Borisov on a strongly hyperbolic orbit as it was to compute the orbit for 175P/Hergenrother, a short-period comet that orbits the sun every ~7 years.</p> <p>In the next article in this series, we will discuss 'astrometry', or the measurement of positions and times using astronomical coordinate systems. Accurate and precise astrometry is the currency by which all orbits are effectively computed.</p>  </div> </article> </div> </div> </div> </div> <footer class="art-footer">  <div class="art-content-layout"> <div class="art-content-layout-row"> <div class="art-layout-cell layout-item-0" style="width: 33%"> <p class="art-page-footer">Copyright © 2023. All Rights Reserved. <br /><a href="https://umd.edu/web-accessibility" onclick="this.target='_blank'">Web Accessibility</a> <br /><a href="http://www.artisteer.com/" onclick="this.target='_blank'">Web Template</a> created with Artisteer.</p> </div> <div class="art-layout-cell layout-item-0" style="width: 34%"> <a href="/misc/contacts.shtml">Contact</a> | <a href="/misc/faqs.shtml">FAQs</a> | <a href="/misc/links.shtml">Links</a> | <a href="/misc/credits.shtml">Credits</a> | <a href="/misc/site-map.shtml">Sitemap</a> | <a href="/misc/search.shtml">Search</a> <br /><a href="https://pdssbn.astro.umd.edu/" title="Small Bodies Node">PDS-SBN</a> | <a href="https://www.minorplanetcenter.net/" title="Minor Planet Center">MPC</a> | <a href="https://www.nasa.gov/planetarydefense" title="Planetary Defense Coordination Office">PDCO</a> | <a href="https://www.cosmos.esa.int/web/smpag" title="Space Mission Planning Advisory Group">SMPAG</a> </div> <div class="art-layout-cell layout-item-0" style="width: 33%"> <p class="art-page-footer">IAWN Mgr: <a href='mailto:tsp%61hr4%34%40g%6D%61%69%6C%2Ec%6Fm'>Dr. Tim Spahr</a> <br />SBN PI: <a href='mailto:gerbsb@astr%6F.umd.%65d%75'>Dr. Gerbs Bauer</a> <br />Web: <a href='mailto:warnerem@%61s%74ro%2Eumd.ed%75'>Elizabeth Warner</a></p> </div> </div> </div>  <p class="art-page-footer">Updated: 30-Oct-2024</p> </footer> </div> </div> </body> </html>