Lines Matching refs:fourier_info
502 static MagickBooleanType ForwardFourier(const FourierInfo *fourier_info, in ForwardFourier() argument
545 magnitude_info=AcquireVirtualMemory((size_t) fourier_info->width, in ForwardFourier()
546 fourier_info->height*sizeof(*magnitude_pixels)); in ForwardFourier()
547 phase_info=AcquireVirtualMemory((size_t) fourier_info->width, in ForwardFourier()
548 fourier_info->height*sizeof(*phase_pixels)); in ForwardFourier()
561 (void) memset(magnitude_pixels,0,fourier_info->width* in ForwardFourier()
562 fourier_info->height*sizeof(*magnitude_pixels)); in ForwardFourier()
564 (void) memset(phase_pixels,0,fourier_info->width* in ForwardFourier()
565 fourier_info->height*sizeof(*phase_pixels)); in ForwardFourier()
566 status=ForwardQuadrantSwap(fourier_info->width,fourier_info->height, in ForwardFourier()
569 status=ForwardQuadrantSwap(fourier_info->width,fourier_info->height,phase, in ForwardFourier()
571 CorrectPhaseLHS(fourier_info->width,fourier_info->height,phase_pixels); in ForwardFourier()
572 if (fourier_info->modulus != MagickFalse) in ForwardFourier()
575 for (y=0L; y < (ssize_t) fourier_info->height; y++) in ForwardFourier()
576 for (x=0L; x < (ssize_t) fourier_info->width; x++) in ForwardFourier()
585 for (y=0L; y < (ssize_t) fourier_info->height; y++) in ForwardFourier()
587 q=GetCacheViewAuthenticPixels(magnitude_view,0L,y,fourier_info->width,1UL, in ForwardFourier()
591 for (x=0L; x < (ssize_t) fourier_info->width; x++) in ForwardFourier()
593 switch (fourier_info->channel) in ForwardFourier()
637 for (y=0L; y < (ssize_t) fourier_info->height; y++) in ForwardFourier()
639 q=GetCacheViewAuthenticPixels(phase_view,0L,y,fourier_info->width,1UL, in ForwardFourier()
643 for (x=0L; x < (ssize_t) fourier_info->width; x++) in ForwardFourier()
645 switch (fourier_info->channel) in ForwardFourier()
692 static MagickBooleanType ForwardFourierTransform(FourierInfo *fourier_info, in ForwardFourierTransform() argument
728 source_info=AcquireVirtualMemory((size_t) fourier_info->width, in ForwardFourierTransform()
729 fourier_info->height*sizeof(*source_pixels)); in ForwardFourierTransform()
737 memset(source_pixels,0,fourier_info->width*fourier_info->height* in ForwardFourierTransform()
741 for (y=0L; y < (ssize_t) fourier_info->height; y++) in ForwardFourierTransform()
743 p=GetCacheViewVirtualPixels(image_view,0L,y,fourier_info->width,1UL, in ForwardFourierTransform()
747 for (x=0L; x < (ssize_t) fourier_info->width; x++) in ForwardFourierTransform()
749 switch (fourier_info->channel) in ForwardFourierTransform()
783 forward_info=AcquireVirtualMemory((size_t) fourier_info->width, in ForwardFourierTransform()
784 (fourier_info->height/2+1)*sizeof(*forward_pixels)); in ForwardFourierTransform()
796 fftw_r2c_plan=fftw_plan_dft_r2c_2d(fourier_info->width,fourier_info->height, in ForwardFourierTransform()
811 gamma=PerceptibleReciprocal((double) fourier_info->width* in ForwardFourierTransform()
812 fourier_info->height); in ForwardFourierTransform()
813 for (y=0L; y < (ssize_t) fourier_info->height; y++) in ForwardFourierTransform()
814 for (x=0L; x < (ssize_t) fourier_info->center; x++) in ForwardFourierTransform()
829 if (fourier_info->modulus != MagickFalse) in ForwardFourierTransform()
830 for (y=0L; y < (ssize_t) fourier_info->height; y++) in ForwardFourierTransform()
831 for (x=0L; x < (ssize_t) fourier_info->center; x++) in ForwardFourierTransform()
838 for (y=0L; y < (ssize_t) fourier_info->height; y++) in ForwardFourierTransform()
839 for (x=0L; x < (ssize_t) fourier_info->center; x++) in ForwardFourierTransform()
858 fourier_info; in ForwardFourierTransformChannel() local
867 fourier_info.width=image->columns; in ForwardFourierTransformChannel()
868 fourier_info.height=image->rows; in ForwardFourierTransformChannel()
873 fourier_info.width=(extent & 0x01) == 1 ? extent+1UL : extent; in ForwardFourierTransformChannel()
875 fourier_info.height=fourier_info.width; in ForwardFourierTransformChannel()
876 fourier_info.center=(ssize_t) (fourier_info.width/2L)+1L; in ForwardFourierTransformChannel()
877 fourier_info.channel=channel; in ForwardFourierTransformChannel()
878 fourier_info.modulus=modulus; in ForwardFourierTransformChannel()
879 magnitude_info=AcquireVirtualMemory((size_t) fourier_info.width, in ForwardFourierTransformChannel()
880 (fourier_info.height/2+1)*sizeof(*magnitude_pixels)); in ForwardFourierTransformChannel()
881 phase_info=AcquireVirtualMemory((size_t) fourier_info.width, in ForwardFourierTransformChannel()
882 (fourier_info.height/2+1)*sizeof(*phase_pixels)); in ForwardFourierTransformChannel()
896 status=ForwardFourierTransform(&fourier_info,image,magnitude_pixels, in ForwardFourierTransformChannel()
899 status=ForwardFourier(&fourier_info,fourier_image,magnitude_pixels, in ForwardFourierTransformChannel()
1107 static MagickBooleanType InverseFourier(FourierInfo *fourier_info, in InverseFourier() argument
1141 magnitude_info=AcquireVirtualMemory((size_t) fourier_info->width, in InverseFourier()
1142 fourier_info->height*sizeof(*magnitude_pixels)); in InverseFourier()
1143 phase_info=AcquireVirtualMemory((size_t) fourier_info->width, in InverseFourier()
1144 fourier_info->height*sizeof(*phase_pixels)); in InverseFourier()
1145 inverse_info=AcquireVirtualMemory((size_t) fourier_info->width, in InverseFourier()
1146 (fourier_info->height/2+1)*sizeof(*inverse_pixels)); in InverseFourier()
1167 for (y=0L; y < (ssize_t) fourier_info->height; y++) in InverseFourier()
1169 p=GetCacheViewVirtualPixels(magnitude_view,0L,y,fourier_info->width,1UL, in InverseFourier()
1173 for (x=0L; x < (ssize_t) fourier_info->width; x++) in InverseFourier()
1175 switch (fourier_info->channel) in InverseFourier()
1209 status=InverseQuadrantSwap(fourier_info->width,fourier_info->height, in InverseFourier()
1211 (void) memcpy(magnitude_pixels,inverse_pixels,fourier_info->height* in InverseFourier()
1212 fourier_info->center*sizeof(*magnitude_pixels)); in InverseFourier()
1215 for (y=0L; y < (ssize_t) fourier_info->height; y++) in InverseFourier()
1217 p=GetCacheViewVirtualPixels(phase_view,0,y,fourier_info->width,1, in InverseFourier()
1221 for (x=0L; x < (ssize_t) fourier_info->width; x++) in InverseFourier()
1223 switch (fourier_info->channel) in InverseFourier()
1256 if (fourier_info->modulus != MagickFalse) in InverseFourier()
1259 for (y=0L; y < (ssize_t) fourier_info->height; y++) in InverseFourier()
1260 for (x=0L; x < (ssize_t) fourier_info->width; x++) in InverseFourier()
1268 CorrectPhaseLHS(fourier_info->width,fourier_info->height,phase_pixels); in InverseFourier()
1270 status=InverseQuadrantSwap(fourier_info->width,fourier_info->height, in InverseFourier()
1272 (void) memcpy(phase_pixels,inverse_pixels,fourier_info->height* in InverseFourier()
1273 fourier_info->center*sizeof(*phase_pixels)); in InverseFourier()
1279 if (fourier_info->modulus != MagickFalse) in InverseFourier()
1280 for (y=0L; y < (ssize_t) fourier_info->height; y++) in InverseFourier()
1281 for (x=0L; x < (ssize_t) fourier_info->center; x++) in InverseFourier()
1293 for (y=0L; y < (ssize_t) fourier_info->height; y++) in InverseFourier()
1294 for (x=0L; x < (ssize_t) fourier_info->center; x++) in InverseFourier()
1309 static MagickBooleanType InverseFourierTransform(FourierInfo *fourier_info, in InverseFourierTransform() argument
1337 source_info=AcquireVirtualMemory((size_t) fourier_info->width, in InverseFourierTransform()
1338 fourier_info->height*sizeof(*source_pixels)); in InverseFourierTransform()
1356 gamma=PerceptibleReciprocal((double) fourier_info->width* in InverseFourierTransform()
1357 fourier_info->height); in InverseFourierTransform()
1358 for (y=0L; y < (ssize_t) fourier_info->height; y++) in InverseFourierTransform()
1359 for (x=0L; x < (ssize_t) fourier_info->center; x++) in InverseFourierTransform()
1373 fftw_c2r_plan=fftw_plan_dft_c2r_2d(fourier_info->width,fourier_info->height, in InverseFourierTransform()
1379 for (y=0L; y < (ssize_t) fourier_info->height; y++) in InverseFourierTransform()
1383 q=GetCacheViewAuthenticPixels(image_view,0L,y,fourier_info->width > in InverseFourierTransform()
1384 image->columns ? image->columns : fourier_info->width,1UL,exception); in InverseFourierTransform()
1387 for (x=0L; x < (ssize_t) fourier_info->width; x++) in InverseFourierTransform()
1390 switch (fourier_info->channel) in InverseFourierTransform()
1443 fourier_info; in InverseFourierTransformChannel() local
1451 fourier_info.width=magnitude_image->columns; in InverseFourierTransformChannel()
1452 fourier_info.height=magnitude_image->rows; in InverseFourierTransformChannel()
1459 fourier_info.width=(extent & 0x01) == 1 ? extent+1UL : extent; in InverseFourierTransformChannel()
1461 fourier_info.height=fourier_info.width; in InverseFourierTransformChannel()
1462 fourier_info.center=(ssize_t) (fourier_info.width/2L)+1L; in InverseFourierTransformChannel()
1463 fourier_info.channel=channel; in InverseFourierTransformChannel()
1464 fourier_info.modulus=modulus; in InverseFourierTransformChannel()
1465 inverse_info=AcquireVirtualMemory((size_t) fourier_info.width, in InverseFourierTransformChannel()
1466 (fourier_info.height/2+1)*sizeof(*inverse_pixels)); in InverseFourierTransformChannel()
1475 status=InverseFourier(&fourier_info,magnitude_image,phase_image, in InverseFourierTransformChannel()
1478 status=InverseFourierTransform(&fourier_info,inverse_pixels,fourier_image, in InverseFourierTransformChannel()