滤波器卷积反射边界
filter convolution reflection boundary
我在进行卷积时如何设置边界遇到了大问题。我们得到了这部分代码,我已经实现了零填充,但现在我必须将像素值反映到边界。
bool ConvolutionFilter::Execute()
{
// First, create a valid output image.
// This fails, if no valid input image is available.
if( !CreateOutputImage() )
{
return false;
}
// Check if convolution kernel is set.
if( m_Kernel == NULL )
{
std::cout << "Error: No convolution kernel set!" << std::endl;
return false;
}
/* Execute the convolution filter. */
/*
* Hier sind schonmal ein paar Konstanten definiert, die Sie brauchen.
* Kommentieren Sie den Code aus, wenn Sie diese benutzen.
*/
const int kernelHalfSizeX = m_Kernel->GetHalfSizeX();
const int kernelHalfSizeY = m_Kernel->GetHalfSizeY();
//const int kernelSizeX = m_Kernel->GetSizeX();
const int imageSizeX = static_cast<int>(m_InputImage->GetSizeX());
const int imageSizeY = static_cast<int>(m_InputImage->GetSizeY());
int conv, value;
/*
* TODO: Aufgabe 3: Implementieren Sie die Faltung des Eingabe-Bildes m_InputImage
* mit der Filtermaske m_Kernel und schreiben sie das Ergebnis in das Ausgabe-Bild
* m_OutputImage.
* Die Faltung ist gegeben als:
* f_out(x,y)=sum from{j=-hY} to{hY} sum from{i=-hX} to{hX} f(x+i,y+j) * g(i,j)
*
* Aufgabe 3b: Verwenden Sie dabei 0-Randbedingungen, d.h. wenn (x+i,y+j) ausserhalb des Bildbereiches
* liegen, ist f(x+i,y+j) * g(i,j) = 0.
*
* Sie koennen mit m_Kernel->GetCoefficient(i,j) auf die Werte der Filtermaske zugreifen.
* Lesen Sie die Kommentare in der Klasse ConvolutionKernel für Details!
*
* Beachten Sie, dass i und j auch negative Werte annehmen, und deshalb nicht unsigned sein
* koennen. Evtl. muessen Sie (abhaengig vom Compiler) das Ergebnis von m_InputImage->GetSizeX()
* bzw. m_InputImage->GetSizeY() zunaechst von unsigned int nach signed int konvertieren, um
* Warnings bei einem Vergleich mit vorzeichenbehafteten Werten zu vermeiden, oder verwenden Sie
* die oben definierten Konstanten imageSizeX, imageSizeY.
*/
/*
* TODO: Aufgabe 3c: Veraendern Sie ihren Code aus Aufg 3b so, dass reflektierende Randbedingungen
* verwendet werden (siehe Aufgabenblatt).
* ACHTUNG: Kommentieren sie den Code aus Aufg 3b nur aus, und loeschen Sie diesen nicht!
*
*/
for(int x = 0; x<=imageSizeX; x++)
{
for(int y = 0; y<=imageSizeY; y++)
{
conv = 0;
for(int j = -kernelHalfSizeY; j<=kernelHalfSizeY; j++)
{
for(int i = -kernelHalfSizeX; i<=kernelHalfSizeX; i++)
{
if(x+i>0 || x+i<imageSizeX || y+j>0 || y+j<imageSizeY) {
value = m_InputImage->GetPixel(x+i,y+j);
}
else
{
value = m_InputImage->GetPixel();
//value = 0; //Zero-Padding
}
conv += value*m_Kernel->GetCoefficient(i,j);
}
}
printf("conv: %dn", conv);
m_OutputImage->SetPixel(x,y,conv);
}
}
return true;
}
那么我如何获取不同蒙版尺寸的像素值以将它们反射到图像的边界上呢?
int xi = x + i;
if ( xi < 0 ) xi = - xi;
else if ( xi > imageSizeX ) xi = imageSizeX - xi;
int yj = y + j;
... same ...
value = m_InputImage->GetPixel(xi,yj);
我认为这不需要解释...
相关文章:
- std::当在256字节边界上写入整数时,流的奇怪行为
- 光线跟踪器灯光反射错误
- 使用不变量来确定二分搜索中的边界条件
- 跨 DLL 边界访问虚拟方法是否安全/可能?
- 在 64 位边界上对齐C++结构数组?
- 如何有效地计算将单位立方体映射到自身的反射和旋转?
- std::数组边界检查如何工作?
- CGAL:如何创建填充边界曲线的曲面网格?
- 在 c++ 中是否允许创建具有运行时边界的数组?
- 是否可以跨 dll 边界返回常量引用/指向 std::vectors?
- 使用即将推出的C++反射工具打印类型的全名
- C++11:16 字节原子<>变量是否在 16 字节边界上自动对齐,从而允许CMPXCHG16B指令?
- 转换引用对象的边界框?
- 将迭代器作为 3 个元素的滑动窗口,可以超调边界(可能使用 Boost)
- 算法设计:用边界数字表示 2D 网格的最佳方式,以C++?
- 查找矩形的边界集
- Assimp 无法加载除 FBX 的漫反射纹理之外的其他纹理
- 在sfml中获取可转换的边界框
- SystemC sc_signal数组内存边界
- 滤波器卷积反射边界