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
