user29a/fastrometry

Fast Automatic World Coordinate Solution Solver. See the github Wiki link below for more info.

fastrometry's Issues

Example for python usage?

How to call fastrometry in a Python script?

Thank you!

ValueError: Buffer has wrong number of dimensions (expected 2, got 3)

Hi all,
when running fastrometry on the file obtained from here: https://fits.gsfc.nasa.gov/nrao_data/samples/hst/w0bu0101t_c0h.fit.gz I get the error below.

% fastrometry w0bs0102t_c0h.fit -ra "RA" -dec "DEC" -scale 0.9 
.
.
.
Extracting sources from image...
Traceback (most recent call last):
  File "/Users/alex/Astronomy/fastrometry/.venv/bin/fastrometry", line 8, in <module>
    sys.exit(callFromCommandLine())
         
  File "/Users/alex/Astronomy/fastrometry/.venv/lib/python3.11/site-packages/fastrometry/fastrometry.py", line 885, in callFromCommandLine
    solution = findWCS(args.filename, ra=args.ra, dec=args.dec, scale=args.scale, scalebnds=args.scalebnds, rotation=args.rotation, rotationbnds=args.rotationbnds, buffer=args.buffer, fieldshape=args.fieldshape, catalogue=args.catalogue, filter=args.filter, pmepoch=args.pmepoch, npts=args.npts, nrefinepts=args.nrefinepts, pixsat=args.pixsat, kernelrad=args.kernelrad, sourcesep=args.sourcesep, vertextol=args.vertextol, nmatchpoints=args.nmatchpoints, nmatchpercent=args.nmatchpercent, wcsdiagnostics=args.wcsdiagnostics, overwrite=args.overwrite, save=args.save, load=args.load, verbosity=args.verbosity, debug=args.debug)

  File "/Users/alex/Astronomy/fastrometry/.venv/lib/python3.11/site-packages/fastrometry/fastrometry.py", line 784, in findWCS
    num_psesources = PSE(img, img_xmax, img_ymax, kernelrad, sourcesep, pixsat, npts, nrefinepts, pixelradius, fieldshape, srcindexmap_initial, srcindexmap_refine, pse_metadata, debug_report, filepath, debug, verbosity)
                     
  File "PSE.pyx", line 528, in fastrometry.cython_code.PSE.PSE
ValueError: Buffer has wrong number of dimensions (expected 2, got 3)

ValueError: Buffer dtype mismatch, expected 'int' but got 'long'

Hello, I keep getting the error below. Could you please advise? The example FITS file is from SDSS.

Thanks,
CS

$ fastrometry frame-g-000756-3-0095.fits -ra "RA" -dec "DEC" -scale 0.4
.
.
.
Extracting sources from image...
Traceback (most recent call last):
  File "/home/developer/.conda/envs/ngpsdev/bin/fastrometry", line 8, in <module>
    sys.exit(callFromCommandLine())
  File "/home/developer/.conda/envs/ngpsdev/lib/python3.10/site-packages/fastrometry/fastrometry.py", line 885, in callFromCommandLine
    solution = findWCS(args.filename, ra=args.ra, dec=args.dec, scale=args.scale, scalebnds=args.scalebnds, rotation=args.rotation, rotationbnds=args.rotationbnds, buffer=args.buffer, fieldshape=args.fieldshape, catalogue=args.catalogue, filter=args.filter, pmepoch=args.pmepoch, npts=args.npts, nrefinepts=args.nrefinepts, pixsat=args.pixsat, kernelrad=args.kernelrad, sourcesep=args.sourcesep, vertextol=args.vertextol, nmatchpoints=args.nmatchpoints, nmatchpercent=args.nmatchpercent, wcsdiagnostics=args.wcsdiagnostics, overwrite=args.overwrite, save=args.save, load=args.load, verbosity=args.verbosity, debug=args.debug)
  File "/home/developer/.conda/envs/ngpsdev/lib/python3.10/site-packages/fastrometry/fastrometry.py", line 784, in findWCS
    num_psesources = PSE(img, img_xmax, img_ymax, kernelrad, sourcesep, pixsat, npts, nrefinepts, pixelradius, fieldshape, srcindexmap_initial, srcindexmap_refine, pse_metadata, debug_report, filepath, debug, verbosity)
  File "src\fastrometry\cython_code\PSE.pyx", line 515, in fastrometry.cython_code.PSE.PSE
ValueError: Buffer dtype mismatch, expected 'int' but got 'long'

Not able to verify known good WCS solutions

Hello Joe, Hello Cameron,
I have spent some time getting to grips with Astrometry and the fastrometry package as my tool of choice. While working on some files of my own I ran into issues, namely fastrometry not being able to produce meaningful results after finding pse sources. The files used are of high enough quality to do meaningful astrometry, yet I did not manage to reproduce result from a different pipeline that was used to come up with a solution. In fact, every file from the source im working on failed.
Judging from the very useful debug outputs included in the package, things go wrong when setting the intermediate points.

You can find one of the files im struggling with here.

Any help is appreciated.
-Eyck

ImportError: cannot import name 'Mapping' from 'collections'

I get the following error when trying to run fastrometry from the terminal. I'm on Ubuntu 18 and installed Python 3.10 to run this, so I'm not sure if the problem is on my end. I also tried to run it inside a conda environment but get the same result.

Traceback (most recent call last): File "/home/orlando/.local/bin/fastrometry", line 5, in <module> from fastrometry.fastrometry import callFromCommandLine File "/home/orlando/.local/lib/python3.10/site-packages/fastrometry/__init__.py", line 1, in <module> from fastrometry.fastrometry import findWCS File "/home/orlando/.local/lib/python3.10/site-packages/fastrometry/fastrometry.py", line 7, in <module> from .catalogue import getIntermediateCoords File "/home/orlando/.local/lib/python3.10/site-packages/fastrometry/catalogue.py", line 4, in <module> from AstraCarta import astracarta File "/home/orlando/.local/lib/python3.10/site-packages/AstraCarta/__init__.py", line 1, in <module> from AstraCarta.AstraCarta import astracarta File "/home/orlando/.local/lib/python3.10/site-packages/AstraCarta/AstraCarta.py", line 1, in <module> from astroquery.simbad import Simbad File "/home/orlando/.local/lib/python3.10/site-packages/astroquery/simbad/__init__.py", line 34, in <module> from .core import Simbad, SimbadClass, SimbadBaseQuery File "/home/orlando/.local/lib/python3.10/site-packages/astroquery/simbad/core.py", line 8, in <module> import requests File "/usr/lib/python3/dist-packages/requests/__init__.py", line 43, in <module> import urllib3 File "/usr/lib/python3/dist-packages/urllib3/__init__.py", line 8, in <module> from .connectionpool import ( File "/usr/lib/python3/dist-packages/urllib3/connectionpool.py", line 29, in <module> from .connection import ( File "/usr/lib/python3/dist-packages/urllib3/connection.py", line 40, in <module> from .util.ssl_ import ( File "/usr/lib/python3/dist-packages/urllib3/util/__init__.py", line 3, in <module> from .connection import is_connection_dropped File "/usr/lib/python3/dist-packages/urllib3/util/connection.py", line 3, in <module> from .wait import wait_for_read File "/usr/lib/python3/dist-packages/urllib3/util/wait.py", line 1, in <module> from .selectors import ( File "/usr/lib/python3/dist-packages/urllib3/util/selectors.py", line 14, in <module> from collections import namedtuple, Mapping ImportError: cannot import name 'Mapping' from 'collections' (/usr/lib/python3.10/collections/__init__.py)