LSSTApplications  16.0-10-g0ee56ad+5,16.0-11-ga33d1f2+5,16.0-12-g3ef5c14+3,16.0-12-g71e5ef5+18,16.0-12-gbdf3636+3,16.0-13-g118c103+3,16.0-13-g8f68b0a+3,16.0-15-gbf5c1cb+4,16.0-16-gfd17674+3,16.0-17-g7c01f5c+3,16.0-18-g0a50484+1,16.0-20-ga20f992+8,16.0-21-g0e05fd4+6,16.0-21-g15e2d33+4,16.0-22-g62d8060+4,16.0-22-g847a80f+4,16.0-25-gf00d9b8+1,16.0-28-g3990c221+4,16.0-3-gf928089+3,16.0-32-g88a4f23+5,16.0-34-gd7987ad+3,16.0-37-gc7333cb+2,16.0-4-g10fc685+2,16.0-4-g18f3627+26,16.0-4-g5f3a788+26,16.0-5-gaf5c3d7+4,16.0-5-gcc1f4bb+1,16.0-6-g3b92700+4,16.0-6-g4412fcd+3,16.0-6-g7235603+4,16.0-69-g2562ce1b+2,16.0-8-g14ebd58+4,16.0-8-g2df868b+1,16.0-8-g4cec79c+6,16.0-8-gadf6c7a+1,16.0-8-gfc7ad86,16.0-82-g59ec2a54a+1,16.0-9-g5400cdc+2,16.0-9-ge6233d7+5,master-g2880f2d8cf+3,v17.0.rc1
LSSTDataManagementBasePackage
Namespaces | Classes | Functions
lsst::ip::isr Namespace Reference

Namespaces

 assembleCcdTask
 
 crosstalk
 
 fringe
 
 isrFunctions
 
 isrQa
 
 isrTask
 
 linearize
 
 masking
 
 measureCrosstalk
 
 straylight
 
 version
 
 vignette
 

Classes

class  CountMaskedPixels
 

Functions

template<typename PixelT >
int applyLookupTable (afw::image::Image< PixelT > &image, ndarray::Array< PixelT, 1, 1 > const &table, PixelT indOffset)
 Add the values in a lookup table to an image, e.g. More...
 
template<typename PixelT >
size_t maskNans (afw::image::MaskedImage< PixelT > const &mi, afw::image::MaskPixel maskVal, afw::image::MaskPixel allow=0)
 Mask NANs in an image. More...
 
template<typename ImagePixelT , typename FunctionT >
void fitOverscanImage (std::shared_ptr< lsst::afw::math::Function1< FunctionT > > &overscanFunction, lsst::afw::image::MaskedImage< ImagePixelT > const &overscan, double ssize=1., int sigma=1)
 
 PYBIND11_MODULE (applyLookupTable, mod)
 
 PYBIND11_MODULE (isr, mod)
 
template int applyLookupTable< float > (afw::image::Image< float > &, ndarray::Array< float, 1, 1 > const &, float)
 
template int applyLookupTable< double > (afw::image::Image< double > &, ndarray::Array< double, 1, 1 > const &, double)
 
std::string between (std::string &s, char ldelim, char rdelim)
 
template void fitOverscanImage (std::shared_ptr< afw::math::Function1< double > > &overscanFunction, afw::image::MaskedImage< float > const &overscan, double ssize, int sigma)
 
template void fitOverscanImage (std::shared_ptr< afw::math::Function1< double > > &overscanFunction, afw::image::MaskedImage< double > const &overscan, double ssize, int sigma)
 
template size_t maskNans< float > (afw::image::MaskedImage< float > const &, afw::image::MaskPixel, afw::image::MaskPixel)
 
template size_t maskNans< double > (afw::image::MaskedImage< double > const &, afw::image::MaskPixel, afw::image::MaskPixel)
 

Function Documentation

◆ applyLookupTable()

template<typename PixelT >
int lsst::ip::isr::applyLookupTable ( afw::image::Image< PixelT > &  image,
ndarray::Array< PixelT, 1, 1 > const &  table,
PixelT  indOffset 
)

Add the values in a lookup table to an image, e.g.

for non-linearity correction

The algorithm is as follows: numOutOfRange = 0 For each i,j of the image: lookupInd = int(indOffset + image[i,j]) if lookupInd not in range [0, table.size() - 1]: set lookupInd to nearest edge and increment numOutOfRange image[i,j] += table[lookupInd] return numOutOfRange

Parameters
[in,out]imageimage to which to add the values; modified in place
[in]tablelookup table
[in]indOffsetscalar added to image value before truncating to lookup column
Returns
the number of pixels whose values were out of range

Definition at line 35 of file applyLookupTable.cc.

39  {
40  if (table.size() == 0u) {
41  throw LSST_EXCEPT(
42  pex::exceptions::LengthError,
43  "Lookup table has zero size."
44  );
45  }
46  int numOutOfRange = 0;
47  int const maxLookupCol = table.size() - 1;
48  for (int col = 0, imHeight = image.getHeight(); col < imHeight; ++col) {
49  for (auto imPtr = image.row_begin(col), end = image.row_end(col); imPtr != end; ++imPtr) {
50  int lookupCol = indOffset + *imPtr;
51  if (lookupCol < 0) {
52  lookupCol = 0;
53  ++numOutOfRange;
54  } else if (lookupCol > maxLookupCol) {
55  lookupCol = maxLookupCol;
56  ++numOutOfRange;
57  }
58  *imPtr += table[lookupCol];
59  }
60  }
61  return numOutOfRange;
62 }
col
int col
Definition: CR.cc:144
LSST_EXCEPT
#define LSST_EXCEPT(type,...)
Create an exception with a given type.
Definition: Exception.h:48
end
int end
Definition: BoundedField.cc:105

◆ applyLookupTable< double >()

template int lsst::ip::isr::applyLookupTable< double > ( afw::image::Image< double > &  ,
ndarray::Array< double, 1, 1 > const &  ,
double   
)

◆ applyLookupTable< float >()

template int lsst::ip::isr::applyLookupTable< float > ( afw::image::Image< float > &  ,
ndarray::Array< float, 1, 1 > const &  ,
float   
)

◆ between()

std::string lsst::ip::isr::between ( std::string &  s,
char  ldelim,
char  rdelim 
)

Definition at line 99 of file Isr.cc.

99  {
100  std::string::iterator b(s.begin());
101  std::string::iterator e(s.end());
102  std::string::iterator lp;
103  std::string::iterator rp;
104 
105  std::string result;
106 
107  if((lp = std::find(b, e, ldelim)) != e)
108  if((rp = std::find(++lp, e, rdelim)) != e)
109  result = std::string(lp, rp);
110 
111  return result;
112 }
b
table::Key< int > b
Definition: TransmissionCurve.cc:468
result
py::object result
Definition: schema.cc:284
std::string::end
T end(T... args)
std::string
STL class.
std::find
T find(T... args)
std::string::begin
T begin(T... args)

◆ fitOverscanImage() [1/3]

template<typename ImagePixelT , typename FunctionT >
void lsst::ip::isr::fitOverscanImage ( std::shared_ptr< lsst::afw::math::Function1< FunctionT > > &  overscanFunction,
lsst::afw::image::MaskedImage< ImagePixelT > const &  overscan,
double  ssize = 1.,
int  sigma = 1 
)

afw::geom::Box2I bbox = afw::geom::Box2I( afw::geom::Point2I(0, y), afw::geom::Point2I(0, width) ); The above was how this was defined before ticket #1556. As I understand it the following is the new way to do this

afw::geom::Box2I bbox = afw::geom::Box2I( afw::geom::Point2I(0, y), afw::geom::Point2I(0, width) ); The above was how this was defined before ticket #1556. As I understand it the following is the new way to do this

Definition at line 52 of file Isr.cc.

57  {
58  typedef afw::image::MaskedImage<ImagePixelT> MaskedImage;
59 
60 
61  const int height = overscan.getHeight();
62  const int width = overscan.getWidth();
63  std::vector<double> values(height);
64  std::vector<double> errors(height);
65  std::vector<double> positions(height);
66 
67  std::vector<double> parameters(overscanFunction->getNParameters(), 0.);
68  std::vector<double> stepsize(overscanFunction->getNParameters(), ssize);
69 
70  for (int y = 0; y < height; ++y) {
77  afw::geom::Box2I bbox = afw::geom::Box2I(afw::geom::Point2I(0,y), afw::geom::Point2I(width,y));
78  MaskedImage mi = MaskedImage(overscan, bbox);
79  afw::math::Statistics stats = afw::math::makeStatistics(*(mi.getImage()), afw::math::MEAN | afw::math::STDEV);
80 
81  values[y] = stats.getValue(afw::math::MEAN);
82  errors[y] = stats.getValue(afw::math::STDEV);
83  positions[y] = y;
84 
85  }
86  afw::math::FitResults fitResults = afw::math::minimize(
87  *overscanFunction,
88  parameters,
89  stepsize,
90  values,
91  errors,
92  positions,
93  sigma
94  );
95 
96  overscanFunction->setParameters(fitResults.parameterList);
97 }
lsst::afw::math::STDEV
estimate sample standard deviation
Definition: Statistics.h:68
y
int y
Definition: SpanSet.cc:49
lsst::afw::math::minimize
FitResults minimize(lsst::afw::math::Function1< ReturnT > const &function, std::vector< double > const &initialParameterList, std::vector< double > const &stepSizeList, std::vector< double > const &measurementList, std::vector< double > const &varianceList, std::vector< double > const &xPositionList, double errorDef)
Find the minimum of a function(x)
Definition: minimize.cc:177
lsst::afw::image::MaskedImage::getHeight
int getHeight() const
Return the number of rows in the image.
Definition: MaskedImage.h:1096
lsst::afw::math::makeStatistics
Statistics makeStatistics(lsst::afw::image::Image< Pixel > const &img, lsst::afw::image::Mask< image::MaskPixel > const &msk, int const flags, StatisticsControl const &sctrl=StatisticsControl())
Handle a watered-down front-end to the constructor (no variance)
Definition: Statistics.h:354
sigma
afw::table::Key< double > sigma
Definition: GaussianPsf.cc:50
bbox
table::Box2IKey bbox
Definition: Detector.cc:166
lsst::geom::Point2I
Point< int, 2 > Point2I
Definition: Point.h:321
lsst::afw::math::MEAN
estimate sample mean
Definition: Statistics.h:67
lsst::afw::image::MaskedImage::getWidth
int getWidth() const
Return the number of columns in the image.
Definition: MaskedImage.h:1094

◆ fitOverscanImage() [2/3]

template void lsst::ip::isr::fitOverscanImage ( std::shared_ptr< afw::math::Function1< double > > &  overscanFunction,
afw::image::MaskedImage< float > const &  overscan,
double  ssize,
int  sigma 
)

◆ fitOverscanImage() [3/3]

template void lsst::ip::isr::fitOverscanImage ( std::shared_ptr< afw::math::Function1< double > > &  overscanFunction,
afw::image::MaskedImage< double > const &  overscan,
double  ssize,
int  sigma 
)

◆ maskNans()

template<typename PixelT >
size_t lsst::ip::isr::maskNans ( afw::image::MaskedImage< PixelT > const &  mi,
afw::image::MaskPixel  maskVal,
afw::image::MaskPixel  allow = 0 
)

Mask NANs in an image.

NANs in the image or variance that are not already masked by the 'allow' value are masked with the 'maskVal'.

Returns
Number of pixels masked
Parameters
miInput image
maskValBit mask value to give a NaN
allowRetain NANs with this bit mask (0 to mask all NANs)

Definition at line 34 of file Isr.cc.

36 {
37  typedef typename afw::image::MaskedImage<PixelT>::x_iterator x_iterator;
38  size_t nPix = 0;
39  for (int y = 0; y != mi.getHeight(); ++y) {
40  for (x_iterator ptr = mi.row_begin(y), end = mi.row_end(y); ptr != end; ++ptr) {
41  if (!(ptr.mask() & allow) && (!std::isfinite(ptr.image()) ||
42  !std::isfinite(ptr.variance()))) {
43  nPix += 1;
44  ptr.mask() |= maskVal;
45  }
46  }
47  }
48  return nPix;
49 }
ptr
uint64_t * ptr
Definition: RangeSet.cc:88
y
int y
Definition: SpanSet.cc:49
std::isfinite
T isfinite(T... args)
end
int end
Definition: BoundedField.cc:105

◆ maskNans< double >()

template size_t lsst::ip::isr::maskNans< double > ( afw::image::MaskedImage< double > const &  ,
afw::image::MaskPixel  ,
afw::image::MaskPixel   
)

◆ maskNans< float >()

template size_t lsst::ip::isr::maskNans< float > ( afw::image::MaskedImage< float > const &  ,
afw::image::MaskPixel  ,
afw::image::MaskPixel   
)

◆ PYBIND11_MODULE() [1/2]

lsst::ip::isr::PYBIND11_MODULE ( applyLookupTable  ,
mod   
)

Definition at line 45 of file applyLookupTable.cc.

45  {
46  declareApplyLookupTable<float>(mod);
47  declareApplyLookupTable<double>(mod);
48 }

◆ PYBIND11_MODULE() [2/2]

lsst::ip::isr::PYBIND11_MODULE ( isr  ,
mod   
)

Definition at line 68 of file isr.cc.

68  {
69  declareAll<float>(mod, "F");
70  declareAll<double>(mod, "D");
71 }
Generated on Fri Feb 15 2019 09:18:51 for LSSTApplications by   doxygen 1.8.13