# Copyright (c) 2012-2023 by the GalSim developers team on GitHub
# https://github.com/GalSim-developers
#
# This file is part of GalSim: The modular galaxy image simulation toolkit.
# https://github.com/GalSim-developers/GalSim
#
# GalSim is free software: redistribution and use in source and binary forms,
# with or without modification, are permitted provided that the following
# conditions are met:
#
# 1. Redistributions of source code must retain the above copyright notice, this
# list of conditions, and the disclaimer given in the accompanying LICENSE
# file.
# 2. Redistributions in binary form must reproduce the above copyright notice,
# this list of conditions, and the disclaimer given in the documentation
# and/or other materials provided with the distribution.
#
from .extra import ExtraOutputBuilder, RegisterExtraOutput
from ..image import ImageS
# The badpix extra output type is currently just a placeholder for when we eventually add
# defects, saturation, etc. Now it always just builds an Image with all 0's.
[docs]class BadPixBuilder(ExtraOutputBuilder):
"""This builds a bad pixel mask image to go along with each regular data image.
There's not much here currently, since GalSim doesn't yet have any image artifacts that
would be appropriate to do something with here. So this is mostly just a placeholder for
when we eventually add defects, saturation, etc.
"""
# The function to call at the end of building each stamp
def processStamp(self, obj_num, config, base, logger):
# Note: This is just a placeholder for now. Once we implement defects, saturation, etc.,
# these features should be marked in the badpix mask. For now though, all pixels = 0.
if base['do_noise_in_stamps']:
badpix_im = ImageS(base['current_stamp'].bounds, wcs=base['wcs'], init_value=0)
self.scratch[obj_num] = badpix_im
# The function to call at the end of building each image
def processImage(self, index, obj_nums, config, base, logger):
image = ImageS(base['image_bounds'], wcs=base['wcs'], init_value=0)
if len(self.scratch) > 0.:
# If we have been accumulating the variance on the stamps, build the total from them.
# Make sure to only use the stamps for objects in this image.
for obj_num in obj_nums:
stamp = self.scratch[obj_num]
b = stamp.bounds & image.bounds
if b.isDefined(): # pragma: no branch
# This next line is equivalent to:
# image[b] |= stamp[b]
# except that this doesn't work through the proxy. We can only call methods
# that don't start with _. Hence using the more verbose form here.
image.setSubImage(b, image.subImage(b) | stamp[b])
else:
# Otherwise, build the bad pixel mask here.
# Again, nothing here yet.
pass
self.data[index] = image
# Register this as a valid extra output
RegisterExtraOutput('badpix', BadPixBuilder())