void SBTopHat::SBTopHatImpl::fillXImage(ImageView<double> im, double x0, double dx, int izero, double y0, double dy, int jzero) const { dbg<<"SBTopHat fillXImage\n"; dbg<<"x = "<<x0<<" + i * "<<dx<<", izero = "<<izero<<std::endl; dbg<<"y = "<<y0<<" + j * "<<dy<<", jzero = "<<jzero<<std::endl; const int m = im.getNCol(); const int n = im.getNRow(); double* ptr = im.getData(); int skip = im.getNSkip(); assert(im.getStep() == 1); // The columns to consider have -r0 <= y < r0 // given that y = y0 + j dy double absdx = std::abs(dx); double absdy = std::abs(dy); int j1 = std::max(0, int(std::ceil(-_r0/absdy - y0/dy))); int j2 = std::min(n, int(std::ceil(_r0/absdy - y0/dy))); y0 += j1 * dy; ptr += j1*im.getStride(); for (int j=j1; j<j2; ++j,y0+=dy,ptr+=skip) { double ysq = y0*y0; double xmax = std::sqrt(_r0sq - ysq); // Set to _norm all pixels with -xmax <= x < xmax // given that x = x0 + i dx. int i1 = std::max(0, int(std::ceil(-xmax/absdx - x0/dx))); int i2 = std::min(m, int(std::ceil(xmax/absdx - x0/dx))); int i=0; for (; i<i1; ++i) ++ptr; for (; i<i2; ++i) *ptr++ = _norm; for (; i<m; ++i) ++ptr; } }
void SBBox::SBBoxImpl::fillXImage(ImageView<double> im, double x0, double dx, int izero, double y0, double dy, int jzero) const { dbg<<"SBBox fillXImage\n"; dbg<<"x = "<<x0<<" + i * "<<dx<<", izero = "<<izero<<std::endl; dbg<<"y = "<<y0<<" + j * "<<dy<<", jzero = "<<jzero<<std::endl; const int m = im.getNCol(); const int n = im.getNRow(); double* ptr = im.getData(); int skip = im.getNSkip(); assert(im.getStep() == 1); x0 /= dx; double wo2 = _wo2 / std::abs(dx); y0 /= dy; double ho2 = _ho2 / std::abs(dy); // Fill the interior with _norm: // Fill pixels where: // x0 + i >= -width/2 // x0 + i < width/2 // y0 + j >= -width/2 // y0 + j < width/2 int i1 = std::max(0, int(std::ceil(-wo2 - x0))); int i2 = std::min(m, int(std::ceil(wo2 - x0))); int j1 = std::max(0, int(std::ceil(-ho2 - y0))); int j2 = std::min(n, int(std::ceil(ho2 - y0))); ptr += im.getStride() * j1 + i1; skip += m - (i2-i1); for (int j=j1; j<j2; ++j,ptr+=skip) { for (int i=i1;i<i2;++i) *ptr++ = _norm; } }