LSST Applications  21.0.0+04719a4bac,21.0.0-1-ga51b5d4+f5e6047307,21.0.0-11-g2b59f77+a9c1acf22d,21.0.0-11-ga42c5b2+86977b0b17,21.0.0-12-gf4ce030+76814010d2,21.0.0-13-g1721dae+760e7a6536,21.0.0-13-g3a573fe+768d78a30a,21.0.0-15-g5a7caf0+f21cbc5713,21.0.0-16-g0fb55c1+b60e2d390c,21.0.0-19-g4cded4ca+71a93a33c0,21.0.0-2-g103fe59+bb20972958,21.0.0-2-g45278ab+04719a4bac,21.0.0-2-g5242d73+3ad5d60fb1,21.0.0-2-g7f82c8f+8babb168e8,21.0.0-2-g8f08a60+06509c8b61,21.0.0-2-g8faa9b5+616205b9df,21.0.0-2-ga326454+8babb168e8,21.0.0-2-gde069b7+5e4aea9c2f,21.0.0-2-gecfae73+1d3a86e577,21.0.0-2-gfc62afb+3ad5d60fb1,21.0.0-25-g1d57be3cd+e73869a214,21.0.0-3-g357aad2+ed88757d29,21.0.0-3-g4a4ce7f+3ad5d60fb1,21.0.0-3-g4be5c26+3ad5d60fb1,21.0.0-3-g65f322c+e0b24896a3,21.0.0-3-g7d9da8d+616205b9df,21.0.0-3-ge02ed75+a9c1acf22d,21.0.0-4-g591bb35+a9c1acf22d,21.0.0-4-g65b4814+b60e2d390c,21.0.0-4-gccdca77+0de219a2bc,21.0.0-4-ge8a399c+6c55c39e83,21.0.0-5-gd00fb1e+05fce91b99,21.0.0-6-gc675373+3ad5d60fb1,21.0.0-64-g1122c245+4fb2b8f86e,21.0.0-7-g04766d7+cd19d05db2,21.0.0-7-gdf92d54+04719a4bac,21.0.0-8-g5674e7b+d1bd76f71f,master-gac4afde19b+a9c1acf22d,w.2021.13
LSST Data Management Base Package
_maskedImageContinued.py
Go to the documentation of this file.
1 # This file is part of afw.
2 #
3 # Developed for the LSST Data Management System.
4 # This product includes software developed by the LSST Project
5 # (https://www.lsst.org).
6 # See the COPYRIGHT file at the top-level directory of this distribution
7 # for details of code ownership.
8 #
9 # This program is free software: you can redistribute it and/or modify
10 # it under the terms of the GNU General Public License as published by
11 # the Free Software Foundation, either version 3 of the License, or
12 # (at your option) any later version.
13 #
14 # This program is distributed in the hope that it will be useful,
15 # but WITHOUT ANY WARRANTY; without even the implied warranty of
16 # MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17 # GNU General Public License for more details.
18 #
19 # You should have received a copy of the GNU General Public License
20 # along with this program. If not, see <https://www.gnu.org/licenses/>.
21 import numpy as np
22 
23 from lsst.utils import TemplateMeta
24 from ._maskedImage import MaskedImageI, MaskedImageF, MaskedImageD, MaskedImageU, MaskedImageL
25 
26 from ..image._fitsIoWithOptions import imageReadFitsWithOptions, exposureWriteFitsWithOptions
27 from ..image._slicing import supportSlicing
28 from ..image._disableArithmetic import disableImageArithmetic
29 
30 __all__ = ["MaskedImage", "VariancePixel"]
31 
32 
33 VariancePixel = np.float32
34 
35 
36 class MaskedImage(metaclass=TemplateMeta):
37 
38  def set(self, value, mask=None, variance=None):
39  """Assign a tuple of scalars to the entirety of all three planes.
40  """
41  if mask is None and variance is None:
42  try:
43  value, mask, variance = value
44  except TypeError:
45  pass
46  if mask is None:
47  mask = 0
48  if variance is None:
49  variance = 0.0
50  self.image.set(value)
51  self.mask.set(mask)
52  self.variance.set(variance)
53 
54  def _set(self, index, value, origin):
55  """Set the pixel at the given index to a triple (value, mask, variance).
56 
57  Parameters
58  ----------
59  index : `geom.Point2I`
60  Position of the pixel to assign to.
61  value : `tuple`
62  A tuple of (value, mask, variance) scalars.
63  origin : `ImageOrigin`
64  Coordinate system of ``index`` (`PARENT` or `LOCAL`).
65  """
66  self.image[index, origin] = value[0]
67  self.mask[index, origin] = value[1]
68  self.variance[index, origin] = value[2]
69 
70  def _get(self, index, origin):
71  """Return a triple (value, mask, variance) at the given index.
72 
73  Parameters
74  ----------
75  index : `geom.Point2I`
76  Position of the pixel to assign to.
77  origin : `ImageOrigin`
78  Coordinate system of ``index`` (`PARENT` or `LOCAL`).
79  """
80  return (self.image[index, origin],
81  self.mask[index, origin],
82  self.variance[index, origin])
83 
84  def getArrays(self):
85  """Return a tuple (value, mask, variance) numpy arrays."""
86  return (self.image.array if self.image else None,
87  self.mask.array if self.mask else None,
88  self.variance.array if self.variance else None)
89 
90  def convertF(self):
91  return MaskedImageF(self, True)
92 
93  def convertD(self):
94  return MaskedImageD(self, True)
95 
96  def __reduce__(self):
97  from lsst.afw.fits import reduceToFits
98  return reduceToFits(self)
99 
100  def __str__(self):
101  string = "image={},\nmask={}, maskPlaneDict={}\nvariance={}, bbox={}"
102  return string.format(self.image.array,
103  self.mask.array,
104  self.mask.getMaskPlaneDict(),
105  self.variance.array,
106  self.getBBox())
107 
108  def __repr__(self):
109  return "{}.{}=({})".format(self.__module__, self.__class__.__name__, str(self))
110 
111  readFitsWithOptions = classmethod(imageReadFitsWithOptions)
112 
113  writeFitsWithOptions = exposureWriteFitsWithOptions
114 
115 
116 MaskedImage.register(np.int32, MaskedImageI)
117 MaskedImage.register(np.float32, MaskedImageF)
118 MaskedImage.register(np.float64, MaskedImageD)
119 MaskedImage.register(np.uint16, MaskedImageU)
120 MaskedImage.register(np.uint64, MaskedImageL)
121 MaskedImage.alias("I", MaskedImageI)
122 MaskedImage.alias("F", MaskedImageF)
123 MaskedImage.alias("D", MaskedImageD)
124 MaskedImage.alias("U", MaskedImageU)
125 MaskedImage.alias("L", MaskedImageL)
126 
127 for cls in set(MaskedImage.values()):
128  supportSlicing(cls)
def format(config, name=None, writeSourceLine=True, prefix="", verbose=False)
Definition: history.py:174