X-Git-Url: https://hackdaworld.org/gitweb/?p=physik%2Fnlsop.git;a=blobdiff_plain;f=dft.c;h=1d26cee7a2a1e6ab25c5cb003fed73d9c189861c;hp=4bbbb46fbc694cda4aaec05f155724168b19cccf;hb=0dcd23ead2402b25e59a63668c16d375f031177d;hpb=e255e57617d57fef9960f1a572b7d707ea21139d

diff --git a/dft.c b/dft.c
index 4bbbb46..1d26cee 100644
--- a/dft.c
+++ b/dft.c
@@ -12,15 +12,16 @@ int main(int argc,char **argv) {
   t_fourier fourier;
   t_bmp src;
   t_bmp dst;
+  t_bmp cut;
   int x,y;
   int offy,offt;
   int scale;
-  int help;
   double max;
   double *mag;
 
   bmp_init(&src,1);
   bmp_init(&dst,1);
+  bmp_init(&cut,1);
   src.mode=READ;
   dst.mode=WRITE;
   strcpy(src.file,argv[1]);
@@ -32,12 +33,14 @@ int main(int argc,char **argv) {
 
   bmp_read_file(&src);
 
-  dst.width=src.info.width;
-  dst.height=src.info.height;
+  bmp_cut_grab_bottom(&cut,&src,src.info.width,GRAB);
+
+  dst.width=cut.width;
+  dst.height=cut.height;
   bmp_alloc_map(&dst);
 
-  fourier.data_len[0]=src.info.width;
-  fourier.data_len[1]=src.info.height;
+  fourier.data_len[0]=cut.width;
+  fourier.data_len[1]=cut.height;
   fourier_alloc_data(&fourier);
 
   mag=(double *)malloc(fourier.data_len[0]*fourier.data_len[1]*sizeof(double));
@@ -82,9 +85,11 @@ int main(int argc,char **argv) {
     offt=0;
     for(y=0;y<dst.height;y++) {
       for(x=0;x<dst.width;x++) {
+	// usual image processing order
         offy=((y<dst.height/2)?y+dst.height/2:y-dst.height/2);
         offy*=dst.width;
         offy+=((x<dst.width/2)?x+dst.width/2:x-dst.width/2);
+	//offy=offt;
         dst.map[offy].r=scale*mag[offt];
         dst.map[offy].g=dst.map[offy].r;
         dst.map[offy].b=dst.map[offy].r;