1 /* Create masks and filter with them.
2 *
3 * Copyright: N. Dessipris 1991,
4 * Written on: Nov 1991
5 * Updated on: Dec 1991
6 * 20/9/95 JC
7 * - modernised
8 * 22/3/10
9 * - gtkdoc
10 */
11
12 /*
13
14 This file is part of VIPS.
15
16 VIPS is free software; you can redistribute it and/or modify
17 it under the terms of the GNU Lesser General Public License as published by
18 the Free Software Foundation; either version 2 of the License, or
19 (at your option) any later version.
20
21 This program is distributed in the hope that it will be useful,
22 but WITHOUT ANY WARRANTY; without even the implied warranty of
23 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
24 GNU Lesser General Public License for more details.
25
26 You should have received a copy of the GNU Lesser General Public License
27 along with this program; if not, write to the Free Software
28 Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA
29 02110-1301 USA
30
31 */
32
33 /*
34
35 These files are distributed with VIPS - http://www.vips.ecs.soton.ac.uk
36
37 */
38
39 #ifdef HAVE_CONFIG_H
40 #include <config.h>
41 #endif /*HAVE_CONFIG_H*/
42 #include <vips/intl.h>
43
44 #include <stdio.h>
45 #include <math.h>
46 #include <stdarg.h>
47
48 #include <vips/vips.h>
49 #include <vips/vips7compat.h>
50 #include <vips/internal.h>
51
52 /* Make a mask image.
53 */
54 static int
build_freq_mask(IMAGE * out,int xs,int ys,ImMaskType flag,va_list ap)55 build_freq_mask( IMAGE *out, int xs, int ys, ImMaskType flag, va_list ap )
56 {
57 /* May be fewer than 4 args ... but extract them all anyway. Should be
58 * safe.
59 */
60 double p0 = va_arg( ap, double );
61 double p1 = va_arg( ap, double );
62 double p2 = va_arg( ap, double );
63 double p3 = va_arg( ap, double );
64 double p4 = va_arg( ap, double );
65
66 VipsImage *t;
67
68 switch( flag ) {
69 case IM_MASK_IDEAL_HIGHPASS:
70 if( vips_mask_ideal( &t, xs, ys, p0,
71 NULL ) )
72 return( -1 );
73 break;
74
75 case IM_MASK_IDEAL_LOWPASS:
76 if( vips_mask_ideal( &t, xs, ys, p0,
77 "reject", TRUE,
78 NULL ) )
79 return( -1 );
80 break;
81
82 case IM_MASK_BUTTERWORTH_HIGHPASS:
83 if( vips_mask_butterworth( &t, xs, ys, p0, p1, p2,
84 NULL ) )
85 return( -1 );
86 break;
87
88 case IM_MASK_BUTTERWORTH_LOWPASS:
89 if( vips_mask_butterworth( &t, xs, ys, p0, p1, p2,
90 "reject", TRUE,
91 NULL ) )
92 return( -1 );
93 break;
94
95 case IM_MASK_GAUSS_HIGHPASS:
96 if( vips_mask_gaussian( &t, xs, ys, p0, p1,
97 NULL ) )
98 return( -1 );
99 break;
100
101 case IM_MASK_GAUSS_LOWPASS:
102 if( vips_mask_gaussian( &t, xs, ys, p0, p1,
103 "reject", TRUE,
104 NULL ) )
105 return( -1 );
106 break;
107
108 case IM_MASK_IDEAL_RINGPASS:
109 if( vips_mask_ideal_ring( &t, xs, ys, p0, p1,
110 NULL ) )
111 return( -1 );
112 break;
113
114 case IM_MASK_IDEAL_RINGREJECT:
115 if( vips_mask_ideal_ring( &t, xs, ys, p0, p1,
116 "reject", TRUE,
117 NULL ) )
118 return( -1 );
119 break;
120
121 case IM_MASK_BUTTERWORTH_RINGPASS:
122 if( vips_mask_butterworth_ring( &t, xs, ys, p0, p1, p2, p3,
123 NULL ) )
124 return( -1 );
125 break;
126
127 case IM_MASK_BUTTERWORTH_RINGREJECT:
128 if( vips_mask_butterworth_ring( &t, xs, ys, p0, p1, p2, p3,
129 "reject", TRUE,
130 NULL ) )
131 return( -1 );
132 break;
133
134 case IM_MASK_GAUSS_RINGPASS:
135 if( vips_mask_gaussian_ring( &t, xs, ys, p0, p1, p2,
136 NULL ) )
137 return( -1 );
138 break;
139
140 case IM_MASK_GAUSS_RINGREJECT:
141 if( vips_mask_gaussian_ring( &t, xs, ys, p0, p1, p2,
142 "reject", TRUE,
143 NULL ) )
144 return( -1 );
145 break;
146
147 case IM_MASK_FRACTAL_FLT:
148 if( vips_mask_fractal( &t, xs, ys, p0,
149 NULL ) )
150 return( -1 );
151 break;
152
153 case IM_MASK_IDEAL_BANDPASS:
154 if( vips_mask_ideal_band( &t, xs, ys, p0, p1, p2,
155 NULL ) )
156 return( -1 );
157 break;
158
159 case IM_MASK_IDEAL_BANDREJECT:
160 if( vips_mask_ideal_band( &t, xs, ys, p0, p1, p2,
161 "reject", TRUE,
162 NULL ) )
163 return( -1 );
164 break;
165
166 case IM_MASK_BUTTERWORTH_BANDPASS:
167 if( vips_mask_butterworth_band( &t, xs, ys, p0, p1, p2, p3, p4,
168 NULL ) )
169 return( -1 );
170 break;
171
172 case IM_MASK_BUTTERWORTH_BANDREJECT:
173 if( vips_mask_butterworth_band( &t, xs, ys, p0, p1, p2, p3, p4,
174 "reject", TRUE,
175 NULL ) )
176 return( -1 );
177 break;
178
179 case IM_MASK_GAUSS_BANDPASS:
180 if( vips_mask_gaussian_band( &t, xs, ys, p0, p1, p2, p3,
181 NULL ) )
182 return( -1 );
183 break;
184
185 case IM_MASK_GAUSS_BANDREJECT:
186 if( vips_mask_gaussian_band( &t, xs, ys, p0, p1, p2, p3,
187 "reject", TRUE,
188 NULL ) )
189 return( -1 );
190 break;
191
192 default:
193 im_error( "im_freq_mask", "%s", _( "unimplemented mask type" ) );
194 return( -1 );
195 }
196
197 if( im_copy( t, out ) ) {
198 g_object_unref( t );
199 return( -1 );
200 }
201 g_object_unref( t );
202
203 return( 0 );
204 }
205
206 int
im_flt_image_freq(IMAGE * in,IMAGE * out,ImMaskType flag,...)207 im_flt_image_freq( IMAGE *in, IMAGE *out, ImMaskType flag, ... )
208 {
209 IMAGE *mask = im_open_local( out, "tempmask", "p" );
210 va_list ap;
211
212 if( !mask )
213 return( -1 );
214
215 /* Generate mask.
216 */
217 va_start( ap, flag );
218 if( build_freq_mask( mask, in->Xsize, in->Ysize, flag, ap ) )
219 return( -1 );
220 va_end( ap );
221
222 if( im_freqflt( in, mask, out ) )
223 return( -1 );
224
225 return( 0 );
226 }
227
228 int
im_create_fmask(IMAGE * out,int xsize,int ysize,ImMaskType flag,...)229 im_create_fmask( IMAGE *out, int xsize, int ysize, ImMaskType flag, ... )
230 {
231 va_list ap;
232
233 va_start( ap, flag );
234 if( build_freq_mask( out, xsize, ysize, flag, ap ) )
235 return( -1 );
236 va_end( ap );
237
238 return( 0 );
239 }
240