Пример #1
0
def findsources(dataid,datapath):
    d = DataSource.objects.filter(id=dataid)
    dfile = '%s%s' % (datapath,d[0].fits)
    #dfile = '%s%s/astrometry-solution.fits' % (path,url)
    dc = fits.open(dfile)
    #### Map WCS coords to pixel values for the first file in the dataset
    f = aplpy.FITSFigure(dc[0])
    head = dc[0].header
    w = wcs.WCS(head)  # This line extracts the WCS info from the header
    #head = dc[0].head[0].header
    print head['TELESCOP']
    ra= head['CRVAL1'] #291.75
    dec = head['CRVAL2'] #1.38
    if (head['TELESCOP'] == 'SQA-0m8' or head['TELESCOP'] == 'SBA-0m8'):
        xpix  = [0,1024,0,1024]
        ypix = [0,0,1024,1024]
        r = 0.2
    if (head['TELESCOP'] == 'Faulkes Telescope North' or head['TELESCOP'] == 'Faulkes Telescope South'):
        r = head['PIXSCALE']*head['CCDXIMSI']/3600. #0.0792178
    # Change to 0,0 1024,1024 if in Northern Hem
        xpix  = [0,head['CCDXIMSI'],0,head['CCDXIMSI']]
        ypix = [0,0,head['CCDXIMSI'],head['CCDXIMSI']]
    #ra1, dec1 = f.pixel2world(xpix,ypix)
    ra1, dec1 = w.wcs_pix2world(xpix, ypix, 1) # argument 1 starts from (1,1)
    ra_max = max(ra1)
    dec_max = max(dec1)
    ra_min = min(ra1)
    dec_min = min(dec1)
   # print head['CCDXIMSI'],head['CCDXIMSI'],ra_min, ra_max, dec_min,dec_max
    t = atpy.Table(catalog='USNO-A2', ra=ra, dec=dec, radius=2.*r, type='vo_conesearch')
    t4 = t.where((t.RA > ra_min) & (t.RA < ra_max) & (t.DEC > dec_min) & (t.DEC < dec_max) )
    print len(t4)
    f.add_grid()
    f.show_grayscale()
    if t4:
        f.show_markers(t4.RA,t4.DEC,c='r', s= 6, alpha=0.5,edgecolor='none')
    imagej = dfile.replace('.fits','.png')
    print imagej
    f.save(imagej)
    #x,y = f.world2pixel(t4.RA,t4.DEC)
    x,y = w.wcs_world2pix(t4.RA,t4.DEC)
    for i,val in enumerate(t4.id):
        #print i, val,int(x[i]),int(y[i]),ra,dec
        if (x[i] >0 and y[i] >0):
            cat = CatSource(name=val,
                            xpos=int(x[i]),
                            ypos=int(y[i]),
                            catalogue='USNO-A2.0',
                            data=d[0])
            try:
                cat.save()
                print "Saved %s" % val
            except:
                print "error on save %s" % val
    f.close()
Пример #2
0
    def handle(self, *args, **options):
        args = [u'wasp12b']
        try:
            planet = Event.objects.get(name=args[0])
            d = DataSource.objects.filter(id=planet.finder)
            print d
            dfile = '%s%s' % (DATA_LOCATION,d[0].fits)
            #dc = pyfits.open(dfile)
            dc = fits.open(dfile)
            #### Map WCS coords to pixel values for the first file in the dataset
            f = aplpy.FITSFigure(dc[0])
            head = dc[0].header
            print head['TELESCOP']
            ra= head['CRVAL1'] #291.75
            dec = head['CRVAL2'] #1.38
            if (head['TELESCOP'] == 'SQA-0m8' or head['TELESCOP'] == 'SBA-0m8'):
                xpix  = [0,1024,0,1024]
                ypix = [0,0,1024,1024]
                r = 0.2
            if (head['TELESCOP'] == 'Faulkes Telescope North' or head['TELESCOP'] == 'Faulkes Telescope South'):
                r = head['PIXSCALE']*head['CCDXIMSI']/3600. #0.0792178
                # Change to 0,0 1024,1024 if in Northern Hem
                xpix  = [0,head['CCDXIMSI'],0,head['CCDXIMSI']]
                ypix = [0,0,head['CCDXIMSI'],head['CCDXIMSI']]
            ra1, dec1 = f.pixel2world(xpix,ypix)
            #ra1, dec1 = w.wcs_pix2world(xpix, ypix, 1) # argument 1 starts from (1,1)
            ra_max = max(ra1)
            dec_max = max(dec1)
            ra_min = min(ra1)
            dec_min = min(dec1)
            print 'dec_min=', dec_min
            #print head['CCDXIMSI'],head['CCDXIMSI'],ra_min, ra_max, dec_min,dec_max
            #t = atpy.Table(catalog='USNO-A2', ra=ra, dec=dec, radius=2.*r, type='vo_conesearch')
            t = Vizier.query_region(coord.SkyCoord(ra=ra, dec=dec,
                                            unit = (u.deg, u.deg),
                                            frame='icrs'),
                                            radius=(2.*r)*u.deg,
                                            catalog='USNO-A2.0'
                                            )
            t_table = t['I/252/out']
            print 't_table is', t_table
            print 'colnames:', t_table.colnames
            #print "t['RA'] > ra_min is", (t['RA'] > ra_min)
            #t4 = np.where((t_table['RAJ2000'] > ra_min) & (t_table['RAJ2000'] < ra_max) & (t_table['DEJ2000'] > dec_min) & (t_table['DEJ2000'] < dec_max))
            #t4 = np.where((t['RA'] > ra_min) & (t['RA'] < ra_max) & (t['DEC'] > dec_min) & (t['DEC'] < dec_max))

            t4 = Table([[],[]],names=['RA','DEC'])
            for i in range(0,len(t_table)):
                print i
                if (t_table['RAJ2000'][i] > ra_min) & (t_table['RAJ2000'][i] < ra_max):
                    if (t_table['DECJ2000'][i] > dec_min) & (t_table['DECJ2000'][i] < dec_max):
                        t4.add_row([t_table['RAJ2000'][i],t_table['DECJ2000'][i]])

            print 't4 is', t4
            print len(t4)
            x,y = f.world2pixel(t4['RA'],t4['DEC'])
            #x,y = w.wcs_world2pix(t4.RA,t4.DEC)
            for i,val in enumerate(t4.id):
                if (x[i] >0 and y[i] >0):
                    cat = CatSource(name=val,
                                 xpos=int(x[i]),
                                 ypos=int(y[i]),
                                 catalogue='USNO-A2.0',
                                 data=d[0])
                    try:
                        cat.save()
                        print "Saved %s" % val
                    except:
                        print "error on save %s" % val
        except Exception,e:
            print "Could not find planet %s - %s" % (args[0],e)