lsst.meas.extensions.astrometryNet.anetAstrometry Namespace Reference

Classes
class	ANetAstrometryConfig
class	ANetAstrometryTask
	Use astrometry.net to match input sources with a reference catalog and solve for the Wcs. More...

Functions
def	showAstrometry (exposure, wcs, allMatches, useMatches, frame=0, title=None, pause=False)
	Show results of astrometry fitting. More...

Function Documentation

showAstrometry()

def lsst.meas.extensions.astrometryNet.anetAstrometry.showAstrometry	(		exposure,
			wcs,
			allMatches,
			useMatches,
			frame = 0,
			title = None,
			pause = False
	)

Show results of astrometry fitting.

Parameters

exposure	Image to display
wcs	Astrometric solution
allMatches	List of all astrometric matches (including rejects)
useMatches	List of used astrometric matches
frame	Frame number for display
title	Title for display
pause	Pause to allow viewing of the display and optional debugging?

• Matches are shown in yellow if used in the Wcs solution, otherwise red
  • +: Detected objects
  • x: Catalogue objects

Definition at line 456 of file anetAstrometry.py.

def showAstrometry(exposure, wcs, allMatches, useMatches, frame=0, title=None, pause=False):
    r"""!Show results of astrometry fitting

    \param exposure Image to display
    \param wcs Astrometric solution
    \param allMatches List of all astrometric matches (including rejects)
    \param useMatches List of used astrometric matches
    \param frame Frame number for display
    \param title Title for display
    \param pause Pause to allow viewing of the display and optional debugging?

    - Matches are shown in yellow if used in the Wcs solution, otherwise red
     - +: Detected objects
     - x: Catalogue objects
    """
    import lsst.afw.display as afwDisplay
    disp = afwDisplay.Display(frame=frame)
    disp.mtv(exposure, title=title)

    useIndices = set(m.second.getId() for m in useMatches)

    radii = []
    with disp.Buffering():
        for i, m in enumerate(allMatches):
            x, y = m.second.getX(), m.second.getY()
            pix = wcs.skyToPixel(m.first.getCoord())

            isUsed = m.second.getId() in useIndices
            if isUsed:
                radii.append(np.hypot(pix[0] - x, pix[1] - y))

            color = afwDisplay.YELLOW if isUsed else afwDisplay.RED

            disp.dot("+", x, y, size=10, ctype=color)
            disp.dot("x", pix[0], pix[1], size=10, ctype=color)

    radii = np.array(radii)
    print("<dr> = %.4g +- %.4g pixels [%d/%d matches]" % (radii.mean(), radii.std(),
                                                          len(useMatches), len(allMatches)))

    if pause:
        import sys
        while True:
            try:
                reply = input("Debugging? [p]db [q]uit; any other key to continue... ").strip()
            except EOFError:
                reply = ""

            if len(reply) > 1:
                reply = reply[0]
            if reply == "p":
                import pdb
                pdb.set_trace()
            elif reply == "q":
                sys.exit(1)
            else:
                break

---

Generated on Thu Aug 8 2019 20:21:22 for LSSTApplications by   1.8.13