def setUp(self): # Load sample input from disk testDir=os.path.dirname(__file__) # Set up local astrometry_net_data helper.setupAstrometryNetDataDir('photocal', verbose=True) self.srcSet = afwTable.SourceCatalog.readFits(os.path.join(testDir, "v695833-e0-c000.xy.fits")) self.bbox = afwGeom.Box2I(afwGeom.Point2I(0, 0), afwGeom.Extent2I(2048, 4612)) # approximate self.exposure = afwImg.ExposureF(os.path.join(testDir, "v695833-e0-c000-a00.sci.fits")) config = measAstrom.ANetBasicAstrometryConfig() logLevel = Log.INFO #logLevel = Log.DEBUG self.astrom = measAstrom.ANetBasicAstrometryTask(config) self.astrom.log.setThreshold(logLevel)
def setUp(self): # Set up local astrometry_net_data self.datapath = setupAstrometryNetDataDir('photocal') self.config = LoadAstrometryNetObjectsTask.ConfigClass() self.bbox = afwGeom.Box2I(afwGeom.Point2I(0, 0), afwGeom.Extent2I(3001, 3001)) self.ctrPix = afwGeom.Point2I(1500, 1500) metadata = PropertySet() metadata.set("RADECSYS", "FK5") metadata.set("EQUINOX", 2000.0) metadata.set("CTYPE1", "RA---TAN") metadata.set("CTYPE2", "DEC--TAN") metadata.set("CUNIT1", "deg") metadata.set("CUNIT2", "deg") metadata.set("CRVAL1", 215.5) metadata.set("CRVAL2", 53.0) metadata.set("CRPIX1", self.ctrPix[0] + 1) metadata.set("CRPIX2", self.ctrPix[1] + 1) metadata.set("CD1_1", 5.1e-05) metadata.set("CD1_2", 0.0) metadata.set("CD2_2", -5.1e-05) metadata.set("CD2_1", 0.0) self.wcs = makeWcs(metadata) self.desNumStarsInPixelBox = 270 self.desNumStarsInSkyCircle = 410
def setUp(self): # Set up local astrometry_net_data self.datapath = helper.setupAstrometryNetDataDir('photocal') self.config = measAstrom.LoadAstrometryNetObjectsTask.ConfigClass() self.bbox = afwGeom.Box2I(afwGeom.Point2I(0, 0), afwGeom.Extent2I(3001, 3001)) self.ctrPix = afwGeom.Point2I(1500, 1500) metadata = dafBase.PropertySet() metadata.set("RADECSYS", "FK5") metadata.set("EQUINOX", 2000.0) metadata.set("CTYPE1", "RA---TAN") metadata.set("CTYPE2", "DEC--TAN") metadata.set("CUNIT1", "deg") metadata.set("CUNIT2", "deg") metadata.set("CRVAL1", 215.5) metadata.set("CRVAL2", 53.0) metadata.set("CRPIX1", self.ctrPix[0] + 1) metadata.set("CRPIX2", self.ctrPix[1] + 1) metadata.set("CD1_1", 5.1e-05) metadata.set("CD1_2", 0.0) metadata.set("CD2_2", -5.1e-05) metadata.set("CD2_1", 0.0) self.wcs = afwImage.makeWcs(metadata) self.desNumStarsInPixelBox = 270 self.desNumStarsInSkyCircle = 410
def setUp(self): # Set up local astrometry_net_data self.datapath = helper.setupAstrometryNetDataDir('photocal') self.bbox = afwGeom.Box2I(afwGeom.Point2I(0, 0), afwGeom.Extent2I(3001, 3001)) self.ctrPix = afwGeom.Point2I(1500, 1500) metadata = dafBase.PropertySet() metadata.set("RADECSYS", "FK5") metadata.set("EQUINOX", 2000.0) metadata.set("CTYPE1", "RA---TAN") metadata.set("CTYPE2", "DEC--TAN") metadata.set("CUNIT1", "deg") metadata.set("CUNIT2", "deg") metadata.set("CRVAL1", 215.5) metadata.set("CRVAL2", 53.0) metadata.set("CRPIX1", self.ctrPix[0] + 1) metadata.set("CRPIX2", self.ctrPix[1] + 1) metadata.set("CD1_1", 5.1e-05) metadata.set("CD1_2", 0.0) metadata.set("CD2_2", -5.1e-05) metadata.set("CD2_1", 0.0) self.tanWcs = afwImage.cast_TanWcs(afwImage.makeWcs(metadata)) self.exposure = afwImage.ExposureF(self.bbox) self.exposure.setWcs(self.tanWcs) self.exposure.setFilter(afwImage.Filter("r", True))
def setUp(self): self.conf = measAstrom.AstrometryConfig() # Load sample input from disk testDir=os.path.dirname(__file__) self.srcCat = afwTable.SourceCatalog.readFits(os.path.join(testDir, "data", "v695833-e0-c000.xy.fits")) # The .xy.fits file has sources in the range ~ [0,2000],[0,4500] # which is bigger than the exposure self.bbox = afwGeom.Box2I(afwGeom.Point2I(0,0), afwGeom.Extent2I(2048, 4612)) smallExposure = afwImage.ExposureF(os.path.join(testDir, "data", "v695833-e0-c000-a00.sci.fits")) self.exposure = afwImage.ExposureF(self.bbox) self.exposure.setWcs(smallExposure.getWcs()) self.exposure.setFilter(smallExposure.getFilter()) self.exposure.setCalib(smallExposure.getCalib()) # Set up local astrometry_net_data helper.setupAstrometryNetDataDir('photocal', rootDir=testDir)
def setUp(self): # Set up local astrometry_net_data datapath = helper.setupAstrometryNetDataDir('photocal') self.conf = measAstrom.ANetBasicAstrometryConfig() # Load andConfig2.py rather than the default. confpath = os.path.join(datapath, 'andConfig2.py') self.andConfig = measAstrom.AstrometryNetDataConfig() self.andConfig.load(confpath)
def getAstrom(self): andpath = helper.setupAstrometryNetDataDir('photocal', rootDir=self.mypath) andcfn = os.path.join(andpath, 'andConfigOpenFiles.py') andconfig = measAstrom.AstrometryNetDataConfig() andconfig.load(andcfn) conf = measAstrom.ANetBasicAstrometryConfig() return measAstrom.ANetBasicAstrometryTask(config=conf, andConfig=andconfig,)
def setUp(self): self.conf = ANetBasicAstrometryConfig() # Load sample input from disk testDir = os.path.dirname(__file__) self.srcCat = afwTable.SourceCatalog.readFits(os.path.join(testDir, "v695833-e0-c000.xy.fits")) # The .xy.fits file has sources in the range ~ [0,2000],[0,4500] self.bbox = afwGeom.Box2I(afwGeom.Point2I(0, 0), afwGeom.Extent2I(2048, 4612)) # approximate self.exposure = afwImg.ExposureF(os.path.join(testDir, "v695833-e0-c000-a00.sci.fits")) # Set up local astrometry_net_data self.an_data_dir = setupAstrometryNetDataDir('photocal')
def setUp(self): self.conf = measAstrom.AstrometryConfig() # Load sample input from disk testDir = os.path.dirname(__file__) self.srcCat = afwTable.SourceCatalog.readFits( os.path.join(testDir, "data", "v695833-e0-c000.xy.fits")) self.srcCat["slot_ApFlux_fluxSigma"] = 1 self.srcCat["slot_PsfFlux_fluxSigma"] = 1 # The .xy.fits file has sources in the range ~ [0,2000],[0,4500] # which is bigger than the exposure self.bbox = afwGeom.Box2I(afwGeom.Point2I(0, 0), afwGeom.Extent2I(2048, 4612)) smallExposure = afwImage.ExposureF( os.path.join(testDir, "data", "v695833-e0-c000-a00.sci.fits")) self.exposure = afwImage.ExposureF(self.bbox) self.exposure.setWcs(smallExposure.getWcs()) self.exposure.setFilter(smallExposure.getFilter()) self.exposure.setCalib(smallExposure.getCalib()) # Set up local astrometry_net_data helper.setupAstrometryNetDataDir('photocal', rootDir=testDir)
def setUp(self): helper.setupAstrometryNetDataDir('t2155')
def setUp(self): helper.setupAstrometryNetDataDir('t2155')
import unittest import lsst.afw.table as afwTable import lsst.utils.tests as utilsTests import lsst.afw.geom as afwGeom import lsst.meas.astrom as measAstrom import lsst.meas.astrom.sip as sip import lsst.meas.astrom.sip.genDistortedImage as distort import lsst.meas.astrom.sip.cleanBadPoints as cleanBadPoints import testFindAstrometryNetDataDir as helper ############################ # Set up local astrometry_net_data helper.setupAstrometryNetDataDir('cfhttemplate') class CreateWcsWithSipCase(unittest.TestCase): def setUp(self): self.config = measAstrom.ANetBasicAstrometryConfig() self.config.defaultFilter = "r" self.astrom = measAstrom.ANetBasicAstrometryTask(config=self.config) testDir=os.path.dirname(__file__) self.filename=os.path.join(testDir, "cat.xy.fits") self.tolArcsec = .4 self.tolPixel = .1 def tearDown(self): del self.config