The mosaic procedure demonstrated here is simmilar to the SDSS example. For more detailed information see the SDSS example. DSS is different from SDSS in that the photographic plates downloaded are large files so most of the mosaic consists of only cropping the downloaded image data.
cd ../pipeline/
/Users/dorislee/Desktop/GSoC2014/rc3-sdss/pipeline
from dss import DSS
from rc3 import RC3
WARNING: AstropyDeprecationWarning: ConfigurationItem has been deprecated in astropy 0.4. Use ConfigItem objects as members of ConfigNamespace subclasses instead. See ConfigNamespace for an example. [astropy.config.configuration] WARNING:astropy:AstropyDeprecationWarning: ConfigurationItem has been deprecated in astropy 0.4. Use ConfigItem objects as members of ConfigNamespace subclasses instead. See ConfigNamespace for an example.
pgc2 = RC3(0.0075,47.2744444444,0.0399805487303,2)
dss = DSS()
rfits = pgc2.mosaic_band('r',0.0075,47.2744444444,3*0.0399805487303,0.0399805487303,2,dss)
------------------mosaic_band---------------------- Now mosaic_band on 2 (0.0075, 47.2744444444, 0.1199416461909) POSSIR POSSIR ['raw_r.fits'] dss_out: raw_r_2 Only one field in region of interest
update = pgc2.source_info(rfits,dss)
1th iteration ------------------source_info---------------------- Source info for DSS_r_2.fits 1 [['PGC 2', 0.0079000000000000008, 47.2744]] ra,dec of catalog sources rc3_data: [[2, 0.0079, 47.2744]] Source is Obvious Biggest Galaxy with radius 194.695274724 pixels! new_ra and new_dec: 0.0064810 , +47.2747672 Radii: 194.695274724 pixel Radii: 0.0214080111832 degrees rc3: 0.0075 , updated: 0.0064810 rc3: 47.2744444444 , updated: +47.2747672 rc3: 0.0399805487303 , updated: 0.0214080111832 Mosaic_all on 2 ------------------mosaic_all_bands---------------------- Working on 2 ------------------mosaic_band---------------------- Now mosaic_band on 2 ('0.0064810', '+47.2747672', 0.1199416461909) POSSIR POSSIR POSSIB ['raw_1b.fits'] dss_out: raw_1b_2 Only one field in region of interest ------------------mosaic_band---------------------- Now mosaic_band on 2 ('0.0064810', '+47.2747672', 0.1199416461909) POSSIR POSSIR POSSIR ['raw_1r.fits'] dss_out: raw_1r_2 Only one field in region of interest ------------------mosaic_band---------------------- Now mosaic_band on 2 ('0.0064810', '+47.2747672', 0.1199416461909) POSSIR POSSIR POSSIIB ['raw_2b.fits'] dss_out: raw_2b_2 Only one field in region of interest ------------------mosaic_band---------------------- Now mosaic_band on 2 ('0.0064810', '+47.2747672', 0.1199416461909) POSSIR POSSIR POSSIIR ['raw_2r.fits'] dss_out: raw_2r_2 Only one field in region of interest ------------------mosaic_band---------------------- Now mosaic_band on 2 ('0.0064810', '+47.2747672', 0.1199416461909) POSSIR POSSIR POSSIIIR ['raw_2ir.fits'] dss_out: raw_2ir_2 Only one field in region of interest Completed Mosaic
/Users/dorislee/.virtualenvs/py27dev/lib/python2.7/site-packages/astroquery/vizier/core.py:560: UserWarning: VOTABLE parsing raised exception: warnings.warn("VOTABLE parsing raised exception: {0}".format(ex))
update
[0.006481, 47.2747672, 0.1199416461909, 0.021408011183161203, 2]
cd ../2
/Users/dorislee/Desktop/GSoC2014/rc3-sdss/pipeline/2
ls
DSS_1b_2.fits DSS_2_BEST.tiff DSS_2b_2.fits DSS_2r_2.fits DSS_1r_2.fits DSS_2_LOW.tiff DSS_2ir_2.fits
from IPython.core.display import Image
import os
os.system("gm convert DSS_2_BEST.tiff DSS_2_BEST.jpg") #for ipynb display purpose only
Image(filename="DSS_2_BEST.jpg")
os.system("gm convert DSS_2_LOW.tiff DSS_2_LOW.jpg") #for ipynb display purpose only
Image(filename="DSS_2_LOW.jpg")