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 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 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.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()
tmc.nativeSys = FOCAL_PLANE.getSysName()
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()
tmc.nativeSys = FOCAL_PLANE.getSysName()
