Flou de défocalisation OpenCV avec noyau personnalisé?

Question

Je souhaite simuler le flou de défocalisation, l'intensité pour chaque pixel d'une image est la suivante:

1/(pi * r^2) pour un rayon r donné, si le pixel est dans sqrt (x^2 + y^2) et 0 sinon (voir le code pour une meilleure explication)

Ceci donne un noyau de flou/convolution de forme circulaire. J'ai essayé de le faire en OpenCV bu avait pas de chance: OpenCV seulement « pixélise » les bords de mon image: testimage http://www.bilderkiste.org/show/original/1131895735815/test_out.jpg

Je ne peux pas vraiment comprendre pourquoi ce qui se passe, voici mon code jusqu'à présent:

//includes. then: 
using namespace std; 
#define KERNELLENGTH 3 
#define PI 3.14159265 
int main() { 
    IplImage *src = 0; 
    IplImage *dst = 0; 
    src = cvLoadImage("test.bmp"); //create image matrixes.. 
    dst = cvLoadImage("test.bmp"); // 
    CvMat *filter; 
    double kernel[KERNELLENGTH * KERNELLENGTH]; //create an appropriate kernel 
    int r = KERNELLENGTH/2; //calculate the radius 
    double value = 1/(PI * KERNELLENGTH * KERNELLENGTH/(4 * r)); //calculate the defocus blur value 
    cout << "Kernel:" << "\n"; 
    for (int x = 0; x < KERNELLENGTH; x++) //calculate kernel (seems to work right!) 
    { 
     for (int y = 0; y < KERNELLENGTH; y++) { 
      if (sqrt((x - KERNELLENGTH/2) * (x - KERNELLENGTH/2) + (y 
        - KERNELLENGTH/2) * (y - KERNELLENGTH/2)) <= r) { 
       kernel[y * 4 + x] = value; //Wert zuweisen 
       cout << value << "\t"; 
      } else 
       cout << 0 << "\t"; 
     } 
     cout << "\n"; 
    } 

    filter = cvCreateMatHeader(KERNELLENGTH, KERNELLENGTH, CV_32FC1);//create the filter 
    cvSetData(filter, kernel, KERNELLENGTH * sizeof(kernel[0]));//link kernel and filter 
    cvFilter2D(src, //convolve filter and src, save to dst 
      dst, filter, cvPoint(-1, -1)); 

    cvSaveImage("test_out.bmp", dst); //save dst on disk 

    cvReleaseImage(&src); 
    cvReleaseImage(&dst); 
    return 0; 
} 

J'apprécierais vraiment de l'aide avec ceci, merci!

Answer 1

Il semble que le problème est dans #define KERNELLENGTH 3, que vous obtenez KERNELLENGTH/2 == 1 et le noyau est 3-par-3 quelque chose, que je n'appellerais pas un disque de défocalisation approprié.

Avez-vous testé avec e.g. #define KERNELLENGTH 10?