Fully Operational?

[Grant Privett]

Hi there,

I've been using Python 3.7.4 under W10 to access Astrometry.net via astroquery.

I'm finding star positions (pixel space) using DAOphot and then creating arrays which I pass to Astrometry.net and employ: 

                # Supply coord list to astrometry.net to create a wcs  

                w=ast.solve_from_source_list( xc,yc,isx,isy,submission_id=submission_id,tweak_order=2,solve_timeout=timeout )

                

                # Pass coordinate to pixel coordinate to ra and dec conversion routine

                try:

                    ra,de=w.wcs_pix2world(xv,yv,0)

                    success=False

                except:

                    success=True

 

To get the RA/Dec of an asteroid in scene. This code worked fine for a week or two, but late Friday night (UTC) started returning WCS header files that would not process. Its now got to the stage where none are processing. I've checked that values of xc,yc,isx,isy plus timeout (1200s - I'm willing to wait) are sensible.

Is there any problem or am I doing something dumb? :)  It wouldn't be the first time?

Grant

[Dustin Lang]

Hi,

Can you please send an example submission ID?

We did have some issues over the weekend (temp disk filling up), but things should be fine now.

PS, isn't it success=True if the function completes and False if there's an exception?

cheers,

--dustin

--
You received this message because you are subscribed to the Google Groups "astrometry" group.
To unsubscribe from this group and stop receiving emails from it, send an email to astrometry+...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/astrometry/8db322c7-dbc2-445a-9ff5-d568120be018o%40googlegroups.com.

[Dustin Lang]

Again --

Please can you tell me the Submission ID for your successful (manual) solve.

thanks,

--dustin

On Mon, Oct 26, 2020 at 1:41 PM Grant Privett <g.pr...@virgin.net> wrote:

Hi Duston,

Thanks for replying.

The success== was a typo from when I tried to reduce something more complicated to its essentials...

I just tried running an image "NGC6729 191018 60s1_P.fit" via the normal manual interface and it worked fine - AstroArt7 showed RA/Dec for any pixel I passed the cursor over.

Then tried the same thing using the Python Astroquery interface - but supplying the star coordinates as derived by DAOphot. Its 120 or so stars.

The code was simplified from a few days ago and is shown below:

        # Kick off astroquery
        ast=AstrometryNet()
        ast.api_key=apikey
        ast.key=apikey
        wcs1=None

        # Use astrometry.net
        success=False
        nloops=0
        while success==False and nloops<3:
          
            # Short delay to avoid bombarding website if something goes wrong
            time.sleep(2)
           
            # Try to plate solve the image  
            try:
                # Supply stars x/y co-ord list plus image dimensions and set tweak to 2
                # with timeout=1200 for astrometry.net to create a wcs from 
                wcs1=ast.solve_from_source_list( \
                     xc,yc,isx,isy, \
                     tweak_order=tweak,solve_timeout=timeout )
                   
            except:
                # This attempt failed. Possibly connection closed, timeout or too few stars?
                print("\nFailed.")
       
            # Check if we have a real WCS by trying to use it
            try:
                # Attempt to get the ra/dec of central image pixel
                xv=int(isx/2)
                yv=int(isy/2)
                ra,de=wcs1.wcs_pix2world(xv,yv,0)
                success=True
               
            except:
                # Problems with WCS
                print("\nReceived WCS corrupt or unpopulated.")
                success=False
         
            # Increment attempts counter
            nloops+=1         
               
        # end of while construct
       
        # Inform user of result
        if nloops<3:
            print("\nCompleted successfully.")
        else:   
            print("\nProblems with WCS created.")
            return     
               

A few days ago I had 2 programs working fine and now I have not had a single solution since Sunday. Would like to resolve this.

Any thoughts very welcome.

Thanks

Grant - Perplexed of Salisbury :)

Virus-free. www.avast.com