BLUR.C
Upload User: gzsenex
Upload Date: 2013-02-24
Package Size: 25984k
Code Size: 1k
Category:

Graph program

Development Platform:

C/C++

  1. #define MAX
  2. #include "fftlib.h"
  3. int
  4. main (int argc, char *argv[])
  5. {
  6. int n=0, nr, nc;
  7. IMAGE im1, p;
  8. COMPLEX_IMAGE fft1, fft2;
  9. float x, xi;
  10. int i,j;
  11. printf ("Image Restoration System module 3 - Blurn");
  12. /* Try to open standard image files.  */
  13. if (argc < 4)
  14. {
  15. printf ("BLUR input psf outputnn");
  16. exit (0);
  17. }
  18. im1 = Input_PBM (argv[1]);
  19. if (im1 == 0)
  20. {
  21.   printf ("No such file as '%s'n", argv[1]);
  22.   exit(1);
  23. }
  24. p = Input_PBM (argv[2]);
  25. if (p == 0)
  26. {
  27.   printf ("No such file as '%s'n", argv[2]);
  28.   exit(1);
  29. }
  30. if (im1->info->nr != im1->info->nc)
  31. {
  32. printf ("Input image is not square. Use ALIGN.n");
  33. exit(0);
  34. }
  35.         if (p->info->nr != p->info->nc)
  36.         {
  37.                 printf ("PSF image is not square. Use ALIGN.n");
  38.                 exit(0);
  39.         }
  40. if (im1->info->nr != p->info->nr)
  41. {
  42. printf ("Input and PSF images differ in size.n");
  43. exit (0);
  44. }
  45. nr = im1->info->nr; nc = im1->info->nc;
  46. image_fftoc (im1, &fft1);
  47. normalize_set();
  48. image_fftoc (p,   &fft2);
  49. normalize_clear();
  50. for (i=0; i<nr; i++)
  51.   for (j=0; j<nc; j++)
  52.   {
  53.     x = fft1[i][j]; xi = fft1[i][j+nc];
  54.     cprod (fft2[i][j], fft2[i][j+nc], &x, &xi);
  55.     fft1[i][j] = x; fft1[i][j+nc] = xi;
  56.   }
  57. freecomplex (fft2); fft2 = (COMPLEX_IMAGE)0;
  58. image_fftinvoc (fft1, &fft2);
  59. freecomplex (fft1); fft1 = (COMPLEX_IMAGE)0;
  60. filt_toint (fft2, p, 0);
  61. Output_PBM (p, argv[3]);
  62. }