def parseCamera(policy):
"""
Make a CameraConfig from a policy
@param policy: Policy object to parse
@return CameraConfig parsed from the policy
"""
camPolicy = policy.get('Camera')
camConfig = CameraConfig()
camConfig.name = camPolicy.get('name')
# Using pixel scale 0.185 arcsec/pixel from:
# http://arxiv.org/pdf/0908.3808v1.pdf
camConfig.plateScale = 13.70 #arcsec/mm
# Radial distortion correction polynomial coeff.
conv_1 = 14805.4
conv_2 = 13619.3
conv_3 = 426637.0
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'inverted'
radialClass = afwGeom.xyTransformRegistry['radial']
tConfig.transform.active.transform.retarget(radialClass)
coeffs = [0., conv_1, conv_2, conv_3]
tConfig.transform.active.transform.coeffs = coeffs
tmc = afwGeom.TransformMapConfig()
tmc.nativeSys = FOCAL_PLANE.getSysName()
tmc.transforms = {PUPIL.getSysName():tConfig}
camConfig.transformDict = tmc
return camConfig
def parseCamera(policy, cameraname):
""" Parse a policy file for a Camera object
@param[in] policy pexPolicy.Policy object to parse
@param[in] cameraname name of camerea being used
@return afw.CameraGeom.CameraConfig object
"""
camPolicy = policy.get('Camera')
camConfig = CameraConfig()
camConfig.name = camPolicy.get('name')
camConfig.plateScale = PIXELSIZE
#Need to invert because this is stored with FOCAL_PLANE to PUPIL as the
#forward transform: only for HSC
if cameraname.lower() == 'hsc':
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'hsc'
else:
#I don't know what the PUPIL transform is for non-HSC cameras is
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'identity'
tpConfig = afwGeom.TransformConfig()
tpConfig.transform.name = 'affine'
tpConfig.transform.active.linear = [1/PIXELSIZE, 0, 0, 1/PIXELSIZE]
tmc = afwGeom.TransformMapConfig()
tmc.nativeSys = FOCAL_PLANE.getSysName()
tmc.transforms = {PUPIL.getSysName():tConfig, FOCAL_PLANE_PIXELS.getSysName():tpConfig}
camConfig.transformDict = tmc
return camConfig
def parseCamera(policy, cameraname):
""" Parse a policy file for a Camera object
@param[in] policy pexPolicy.Policy object to parse
@param[in] cameraname name of camerea being used
@return afw.CameraGeom.CameraConfig object
"""
camPolicy = policy.get('Camera')
camConfig = CameraConfig()
camConfig.name = camPolicy.get('name')
camConfig.plateScale = PIXELSIZE
#Need to invert because this is stored with FOCAL_PLANE to PUPIL as the
#forward transform: only for HSC
if cameraname.lower() == 'hsc':
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'hsc'
else:
#I don't know what the PUPIL transform is for non-HSC cameras is
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'identity'
tpConfig = afwGeom.TransformConfig()
tpConfig.transform.name = 'affine'
tpConfig.transform.active.linear = [1 / PIXELSIZE, 0, 0, 1 / PIXELSIZE]
tmc = afwGeom.TransformMapConfig()
tmc.nativeSys = FOCAL_PLANE.getSysName()
tmc.transforms = {
PUPIL.getSysName(): tConfig,
FOCAL_PLANE_PIXELS.getSysName(): tpConfig
}
camConfig.transformDict = tmc
return camConfig
def ReturnCamera(baseDir):
"""
This method reads in the files
baseDir/focalplanelayout.txt
baseDir/segmentation.txt
and returns an afw.cameraGeom object
Below is the original documentation of the function this code was copied from:
Create the configs for building a camera. This runs on the files distributed with PhoSim. Currently gain and
saturation need to be supplied as well. The file should have three columns: on disk amp id (R22_S11_C00), gain, saturation.
For example:
DetectorLayoutFile -- https://dev.lsstcorp.org/cgit/LSST/sims/phosim.git/plain/data/lsst/focalplanelayout.txt?h=dev
SegmentsFile -- https://dev.lsstcorp.org/cgit/LSST/sims/phosim.git/plain/data/lsst/segmentation.txt?h=dev
"""
defaultOutDir = 'scratch'
DetectorLayoutFile = os.path.join(baseDir, 'focalplanelayout.txt')
SegmentsFile = os.path.join(baseDir, 'segmentation.txt')
GainFile = None
phosimVersion='1.0'
ampTableDict = makeAmpTables(SegmentsFile, GainFile)
detectorConfigList = makeDetectorConfigs(DetectorLayoutFile, phosimVersion)
#Build the camera config.
camConfig = CameraConfig()
camConfig.detectorList = dict([(i,detectorConfigList[i]) for i in xrange(len(detectorConfigList))])
camConfig.name = 'LSST'
camConfig.plateScale = 2.0 #arcsec per mm
pScaleRad = afwGeom.arcsecToRad(camConfig.plateScale)
pincushion = 0.925
# Don't have this yet ticket/3155
#camConfig.boresiteOffset_x = 0.
#camConfig.boresiteOffset_y = 0.
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'inverted'
radialClass = afwGeom.xyTransformRegistry['radial']
tConfig.transform.active.transform.retarget(radialClass)
# According to Dave M. the simulated LSST transform is well approximated (1/3 pix)
# by a scale and a pincusion.
#this is ultimately used to convert from focal plane coordinates to pupil coordinates
#see the asgnment below to tmc.transforms
tConfig.transform.active.transform.coeffs = [0., 1./pScaleRad, 0., pincushion/pScaleRad]
#tConfig.transform.active.boresiteOffset_x = camConfig.boresiteOffset_x
#tConfig.transform.active.boresiteOffset_y = camConfig.boresiteOffset_y
tmc = afwGeom.TransformMapConfig()
tmc.nativeSys = FOCAL_PLANE.getSysName()
tmc.transforms = {PUPIL.getSysName():tConfig}
camConfig.transformDict = tmc
myCamera = makeCameraFromCatalogs(camConfig, ampTableDict)
return myCamera
def _getCameraFromPath(path):
"""Return the camera geometry given solely the path to the location
of that definition."""
config = CameraConfig()
config.load(os.path.join(path, "camera.py"))
return makeCameraFromPath(
cameraConfig=config,
ampInfoPath=path,
shortNameFunc=lambda name: name.replace(" ", "_"),
)
def getCamera(self):
"Grab camera information from camera/camera.py file."
path = os.path.join(getPackageDir("obs_vista"), "camera")
config = CameraConfig()
config.load(os.path.join(path, "camera.py"))
return makeCameraFromPath(
cameraConfig=config,
ampInfoPath=path,
shortNameFunc=lambda name: name.replace(" ", "_"),
)
def loadCamera(repoDir):
"""Load a camera given the path to the repository (the directory that contains the
description directory). I use this just in testing from the interpreter prompt.
In general, it's probably best to do butler.get('camera')
@param repoDir: path to the root of the camera description tree
"""
inputPath = os.path.join(repoDir, "description", "camera")
camConfigPath = os.path.join(inputPath, "camera.py")
camConfig = CameraConfig()
camConfig.load(camConfigPath)
lsstSimMapper = LsstSimMapper
return makeCameraFromPath(camConfig, inputPath, lsstSimMapper.getShortCcdName)
def loadCamera(repoDir):
"""Load a camera given the path to its description
I use this just in testing from the interpreter prompt.
In general, it's probably best to do butler.get('camera')
@param repoDir: path to the root of the camera description tree
"""
inputPath = os.path.join(repoDir, "description", "camera")
camConfigPath = os.path.join(inputPath, "camera.py")
camConfig = CameraConfig()
camConfig.load(camConfigPath)
lsstSimMapper = LsstSimMapper
return makeCameraFromPath(camConfig, inputPath,
lsstSimMapper.getShortCcdName)
def getCamera(self):
"""Retrieve the cameraGeom representation of HSC.
This is a temporary API that should go away once obs_ packages have
a standardized approach to writing versioned cameras to a Gen3 repo.
"""
path = os.path.join(getPackageDir("obs_subaru"), "hsc", "camera")
config = CameraConfig()
config.load(os.path.join(path, "camera.py"))
return makeCameraFromPath(
cameraConfig=config,
ampInfoPath=path,
shortNameFunc=lambda name: name.replace(" ", "_"),
pupilFactoryClass=HscPupilFactory)
def getCamera(self):
"""Retrieve the cameraGeom representation of HSC.
This is a temporary API that should go away once obs_ packages have
a standardized approach to writing versioned cameras to a Gen3 repo.
"""
path = os.path.join(getPackageDir("obs_subaru"), "hsc", "camera")
config = CameraConfig()
config.load(os.path.join(path, "camera.py"))
return makeCameraFromPath(
cameraConfig=config,
ampInfoPath=path,
shortNameFunc=lambda name: name.replace(" ", "_"),
pupilFactoryClass=HscPupilFactory
)
def ReturnCamera(baseDir):
"""
This method reads in the files
baseDir/focalplanelayout.txt
baseDir/segmentation.txt
and returns an afw.cameraGeom object
Below is the original documentation of the function this code was copied from:
Create the configs for building a camera. This runs on the files distributed with PhoSim. Currently gain and
saturation need to be supplied as well. The file should have three columns: on disk amp id (R22_S11_C00), gain, saturation.
For example:
DetectorLayoutFile -- https://dev.lsstcorp.org/cgit/LSST/sims/phosim.git/plain/data/lsst/focalplanelayout.txt?h=dev
SegmentsFile -- https://dev.lsstcorp.org/cgit/LSST/sims/phosim.git/plain/data/lsst/segmentation.txt?h=dev
"""
defaultOutDir = 'scratch'
DetectorLayoutFile = os.path.join(baseDir, 'focalplanelayout.txt')
SegmentsFile = os.path.join(baseDir, 'segmentation.txt')
GainFile = None
phosimVersion='1.0'
ampTableDict = makeAmpTables(SegmentsFile, GainFile)
detectorConfigList = makeDetectorConfigs(DetectorLayoutFile, phosimVersion)
#Build the camera config.
camConfig = CameraConfig()
camConfig.detectorList = dict([(i,detectorConfigList[i]) for i in range(len(detectorConfigList))])
camConfig.name = 'LSST'
camConfig.plateScale = 2.0 #arcsec per mm
pScaleRad = afwGeom.arcsecToRad(camConfig.plateScale)
pincushion = 0.925
# Don't have this yet ticket/3155
#camConfig.boresiteOffset_x = 0.
#camConfig.boresiteOffset_y = 0.
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'inverted'
radialClass = afwGeom.xyTransformRegistry['radial']
tConfig.transform.active.transform.retarget(radialClass)
# According to Dave M. the simulated LSST transform is well approximated (1/3 pix)
# by a scale and a pincusion.
#this is ultimately used to convert from focal plane coordinates to pupil coordinates
#see the asgnment below to tmc.transforms
tConfig.transform.active.transform.coeffs = [0., 1./pScaleRad, 0., pincushion/pScaleRad]
#tConfig.transform.active.boresiteOffset_x = camConfig.boresiteOffset_x
#tConfig.transform.active.boresiteOffset_y = camConfig.boresiteOffset_y
tmc = afwGeom.TransformMapConfig()
tmc.nativeSys = FOCAL_PLANE.getSysName()
tmc.transforms = {PUPIL.getSysName():tConfig}
camConfig.transformDict = tmc
myCamera = makeCameraFromCatalogs(camConfig, ampTableDict)
return myCamera
def makeCamera(name="SDSS", outputDir=None):
"""Make a camera
@param name: name of the camera
@param outputDir: If not None, write the objects used to make the camera to this location
@return a camera object
"""
camConfig = CameraConfig()
camConfig.name = name
camConfig.detectorList = {}
camConfig.plateScale = 16.5 # arcsec/mm
pScaleRad = afwGeom.arcsecToRad(camConfig.plateScale)
radialDistortCoeffs = [0.0, 1.0 / pScaleRad]
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'inverted'
radialClass = afwGeom.xyTransformRegistry['radial']
tConfig.transform.active.transform.retarget(radialClass)
tConfig.transform.active.transform.coeffs = radialDistortCoeffs
tmc = afwGeom.TransformMapConfig()
tmc.nativeSys = FOCAL_PLANE.getSysName()
tmc.transforms = {PUPIL.getSysName(): tConfig}
camConfig.transformDict = tmc
ccdId = 0
ampInfoCatDict = {}
for i in range(6):
dewarName = str(i + 1)
filters = "riuzg"
for j, c in enumerate(reversed(filters)):
ccdName = "%s%s" % (c, dewarName)
offsetPoint = afwGeom.Point2D(25.4 * 2.5 * (2.5 - i),
25.4 * 2.1 * (2.0 - j))
ccdInfo = makeCcd(ccdName, ccdId, offsetPoint)
ampInfoCatDict[ccdName] = ccdInfo['ampInfo']
camConfig.detectorList[ccdId] = ccdInfo['ccdConfig']
ccdId += 1
if outputDir is not None:
camConfig.save(os.path.join(outputDir, 'camera.py'))
for k in ampInfoCatDict:
ampInfoCatDict[k].writeFits(
os.path.join(outputDir, "%s.fits" % (k)))
return makeCameraFromCatalogs(camConfig, ampInfoCatDict)
def makeCamera(name="SDSS", outputDir=None):
"""Make a camera
@param name: name of the camera
@param outputDir: If not None, write the objects used to make the camera to this location
@return a camera object
"""
camConfig = CameraConfig()
camConfig.name = name
camConfig.detectorList = {}
camConfig.plateScale = 16.5 # arcsec/mm
pScaleRad = afwGeom.arcsecToRad(camConfig.plateScale)
radialDistortCoeffs = [0.0, 1.0/pScaleRad]
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'inverted'
radialClass = afwGeom.transformRegistry['radial']
tConfig.transform.active.transform.retarget(radialClass)
tConfig.transform.active.transform.coeffs = radialDistortCoeffs
tmc = TransformMapConfig()
tmc.nativeSys = FOCAL_PLANE.getSysName()
tmc.transforms = {FIELD_ANGLE.getSysName(): tConfig}
camConfig.transformDict = tmc
ccdId = 0
ampInfoCatDict = {}
for i in range(6):
dewarName = str(i+1)
filters = "riuzg"
for j, c in enumerate(reversed(filters)):
ccdName = "%s%s" % (c, dewarName)
offsetPoint = afwGeom.Point2D(25.4*2.5*(2.5-i), 25.4*2.1*(2.0 - j))
ccdInfo = makeCcd(ccdName, ccdId, offsetPoint)
ampInfoCatDict[ccdName] = ccdInfo['ampInfo']
camConfig.detectorList[ccdId] = ccdInfo['ccdConfig']
ccdId += 1
if outputDir is not None:
camConfig.save(os.path.join(outputDir, 'camera.py'))
for k in ampInfoCatDict:
ampInfoCatDict[k].writeFits(os.path.join(outputDir, "%s.fits"%(k)))
return makeCameraFromCatalogs(camConfig, ampInfoCatDict)
parser = argparse.ArgumentParser()
parser.add_argument("DetectorLayoutFile", help="Path to detector layout file")
parser.add_argument("SegmentsFile", help="Path to amp segments file")
parser.add_argument("OutputDir",
help = "Path to dump configs and AmpInfo Tables; defaults to %r" % (defaultOutDir,),
nargs = "?",
default = defaultOutDir,
)
parser.add_argument("--clobber", action="store_true", dest="clobber", default=False,
help=("remove and re-create the output directory if it already exists?"))
args = parser.parse_args()
ampTableDict = makeAmpTables(args.SegmentsFile)
detectorConfigList = makeDetectorConfigs(args.DetectorLayoutFile)
#Build the camera config.
camConfig = CameraConfig()
camConfig.detectorList = dict([(i,detectorConfigList[i]) for i in xrange(len(detectorConfigList))])
camConfig.name = 'DECAM'
#From DECam calibration doc
camConfig.plateScale = 17.575
pScaleRad = afwGeom.arcsecToRad(camConfig.plateScale)
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'radial'
nomWavelen = 0.625 #nominal wavelen in microns
coeff0 = 0
coeff1 = 1 - 2.178e-4 - 2.329e-4/nomWavelen + 4.255e-5/nomWavelen**2
coeff2 = 0
coeff3 = -6.174e-8 + 5.569e-9/nomWavelen
tConfig.transform.active.coeffs = [pScaleRad*coeff0, pScaleRad*coeff1, pScaleRad*coeff2,
pScaleRad*coeff3]
tmc = afwGeom.TransformMapConfig()
help = "Path to dump configs and AmpInfo Tables; defaults to %r" % (defaultOutDir,),
nargs = "?",
default = defaultOutDir,
)
parser.add_argument("--clobber", action="store_true", dest="clobber", default=False,
help=("remove and re-create the output directory if it already exists?"))
args = parser.parse_args()
ampTableDict = makeAmpTables(args.SegmentsFile, args.GainFile)
if args.phosimVersion is None:
phosimVersion = getPhosimVersion(defaultDataDir)
else:
phosimVersion = args.phosimVersion
detectorConfigList = makeDetectorConfigs(args.DetectorLayoutFile, phosimVersion)
# Build the camera config.
camConfig = CameraConfig()
camConfig.detectorList = dict([(i, detectorConfigList[i]) for i in xrange(len(detectorConfigList))])
camConfig.name = 'LSST'
camConfig.plateScale = 20.0
pScaleRad = afwGeom.arcsecToRad(camConfig.plateScale)
pincushion = 0.925
# Don't have this yet ticket/3155
# camConfig.boresiteOffset_x = 0.
# camConfig.boresiteOffset_y = 0.
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'inverted'
radialClass = afwGeom.xyTransformRegistry['radial']
tConfig.transform.active.transform.retarget(radialClass)
# According to Dave M. the simulated LSST transform is well approximated (1/3 pix)
# by a scale and a pincusion.
tConfig.transform.active.transform.coeffs = [0., 1./pScaleRad, 0., pincushion/pScaleRad]
parser = argparse.ArgumentParser()
parser.add_argument("DetectorLayoutFile", help="Path to detector layout file")
parser.add_argument("SegmentsFile", help="Path to amp segments file")
parser.add_argument("OutputDir",
help="Path to dump configs and AmpInfo Tables; defaults to %r" % (defaultOutDir,),
nargs="?",
default=defaultOutDir,
)
parser.add_argument("--clobber", action="store_true", dest="clobber", default=False,
help=("remove and re-create the output directory if it already exists?"))
args = parser.parse_args()
ampTableDict = makeAmpTables(args.SegmentsFile)
detectorConfigList = makeDetectorConfigs(args.DetectorLayoutFile)
# Build the camera config.
camConfig = CameraConfig()
camConfig.detectorList = dict([(i, detectorConfigList[i]) for i in range(len(detectorConfigList))])
camConfig.name = 'DECAM'
# From DECam calibration doc
camConfig.plateScale = 17.575
pScaleRad = afwGeom.arcsecToRad(camConfig.plateScale)
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'radial'
nomWavelen = 0.625 # nominal wavelen in microns
coeff0 = 0
coeff1 = 1 - 2.178e-4 - 2.329e-4/nomWavelen + 4.255e-5/nomWavelen**2
coeff2 = 0
coeff3 = -6.174e-8 + 5.569e-9/nomWavelen
tConfig.transform.active.coeffs = [pScaleRad*coeff0, pScaleRad*coeff1, pScaleRad*coeff2,
pScaleRad*coeff3]
tmc = afwGeom.TransformMapConfig()
help="Path to dump configs and AmpInfo Tables; defaults to %r" % (defaultOutDir,),
nargs="?",
default=defaultOutDir,
)
parser.add_argument("--clobber", action="store_true", dest="clobber", default=False,
help=("remove and re-create the output directory if it already exists?"))
args = parser.parse_args()
ampTableDict = makeAmpTables(args.SegmentsFile, args.GainFile)
if args.phosimVersion is None:
phosimVersion = getPhosimVersion(defaultDataDir)
else:
phosimVersion = args.phosimVersion
detectorConfigList = makeDetectorConfigs(args.DetectorLayoutFile, phosimVersion)
# Build the camera config.
camConfig = CameraConfig()
camConfig.detectorList = dict([(i, detectorConfigList[i]) for i in range(len(detectorConfigList))])
camConfig.name = 'LSST'
camConfig.plateScale = 20.0
pScaleRad = afwGeom.arcsecToRad(camConfig.plateScale)
pincushion = 0.925
# Don't have this yet ticket/3155
# camConfig.boresiteOffset_x = 0.
# camConfig.boresiteOffset_y = 0.
tConfig = afwGeom.TransformConfig()
tConfig.transform.name = 'inverted'
radialClass = afwGeom.transformRegistry['radial']
tConfig.transform.active.transform.retarget(radialClass)
# According to Dave M. the simulated LSST transform is well approximated (1/3 pix)
# by a scale and a pincusion.
tConfig.transform.active.transform.coeffs = [0., 1./pScaleRad, 0., pincushion/pScaleRad]
