libs/interpolate/fresize.h

#ifndef __FRESIZE_H
#define __FRESIZE_H
#include <complex.h>
#include <fftw3.h>

struct fresize_struct {
  int method;
  int nsub;
  int hwidth;
  float *ker,*kerx,*kery;
  fftwf_complex *helpc,*fkernel,*fkernely;
  float *helpin,*helpout;
  fftwf_plan plan1,plan2,plan1y,plan2y;
  int nxin,nyin,nxinp,nyinp;
};

int init_fresize_sample(
  struct fresize_struct *pars,
  int nsub // Distance between sampled points
);

int init_fresize_bin(
  struct fresize_struct *pars,
  int nsub // Binsize
);

int init_fresize_boxcar(
  struct fresize_struct *pars,
  int hwidth, // Half width of boxcar. Full is 2*hwidth+1.
  int nsub // Distance between sampled points
);

int init_fresize_boxcar_fft(
  struct fresize_struct *pars,
  int hwidth, // Half width of boxcar. Full is 2*hwidth+1.
  int nsub, // Distance between sampled points
  int nxin, // Array size
  int nyin // Array size
);

int init_fresize_gaussian( // Simple square truncated Gaussian
  struct fresize_struct *pars,
  float sigma, // Shape is exp(-(d/sigma)^2/2)
  int hwidth, // Half (truncation) width of kernel. Full is 2*hwidth+1.
  int nsub // Distance between sampled points
);

int init_fresize_gaussian_fft( // Simple square truncated Gaussian. FFT version.
  struct fresize_struct *pars,
  float sigma, // Shape is exp(-(d/sigma)^2/2)
  int hwidth, // Half (truncation) width of kernel. Full is 2*hwidth+1.
  int nsub, // Distance between sampled points
  int nxin, // Array size
  int nyin // Array size
);

int init_fresize_sinc( // Sinc filter
  struct fresize_struct *pars,
  float wsinc, /* Shape is sinc(d/wsinc)*ap(d)
                  wsinc is the amount by which the Nyquist is reduced.
                  May want wsinc=nsub. */
  int hwidth, // Half width of kernel. Full is 2*hwidth+1.
  int iap, /* Apodization method. Always ap=0 for d>nap*wsinc.
              iap=0 means no apodization ap=1
              iap=1 uses parabola ap=1-(d/(nap*wsinc))^2
              iap=2 uses sinc ap=sinc(d/(nap*wsinc))
              all other cases give ap=1 (not guaranteed) */
  int nap, /* Sinc apodization width in units of wsinc.
              Normally hwidth=nap*wsinc,
              but hwidth=nap*wsinc-1 works for integer */
  int nsub // Distance between sampled points
);

int init_fresize_sinc_fft( // Sinc filter. FFT version.
  struct fresize_struct *pars,
  float wsinc, /* Shape is sinc(d/wsinc)*ap(d)
                  wsinc is the amount by which the Nyquist is reduced.
                  May want wsinc=nsub. */
  int hwidth, // Half width of kernel. Full is 2*hwidth+1.
  int iap, /* Apodization method. Always ap=0 for d>nap*wsinc.
              iap=0 means no apodization ap=1
              iap=1 uses parabola ap=1-(d/(nap*wsinc))^2
              iap=2 uses sinc ap=sinc(d/(nap*wsinc))
              all other cases give ap=1 (not guaranteed) */
  int nap, /* Sinc apodization width in units of wsinc.
              Normally hwidth=nap*wsinc,
              but hwidth=nap*wsinc-1 works for integer */
  int nsub, // Distance between sampled points
  int nxin, // Array size
  int nyin // Array size
);

int init_fresize_gaussian2( // Circularly truncated Gaussian
  struct fresize_struct *pars,
  float sigma, // Shape is exp(-(d/sigma)^2/2)
  float rmax, // Truncation radius. Probably rmax<=hwidth.
  int hwidth, // Half (truncation) width of kernel. Full is 2*hwidth+1.
  int nsub // Distance between sampled points
);

int init_fresize_gaussian2_fft( // Circularly truncated Gaussian. FFT version.
  struct fresize_struct *pars,
  float sigma, // Shape is exp(-(d/sigma)^2/2)
  float rmax, // Truncation radius. Probably rmax<=hwidth.
  int hwidth, // Half (truncation) width of kernel. Full is 2*hwidth+1.
  int nsub, // Distance between sampled points
  int nxin, // Input size
  int nyin // Input size
);

int init_fresize_airy( // 2D Airy filter
  struct fresize_struct *pars,
  float cdown, /* cdown is the amount by which the Nyquist is reduced. */
  int hwidth, /* Half width of kernel. Full is 2*hwidth+1.
                 Set to <0 to make routine set appropriate value */
  int iap, /* Apodization method. Always ap=0 for d>Z_nap, where Z_nap os
              the position of the nap'th zero.
              iap=0 means no apodization ap=1
              iap=1 uses parabola ap=1-(d/Z_nap)^2
              iap=2 uses sinc ap=sinc(d/(Z_nap))
              iap=3 uses Airy with first zero at Z_nap
              all other cases give ap=1 (not guaranteed) */
  int nap, /* Apodizes to nap'th zero */
  int nsub // Distance between sampled points
);

int init_fresize_airy_fft( // 2D Airy filter. FFT version.
  struct fresize_struct *pars,
  float cdown, /* cdown is the amount by which the Nyquist is reduced. */
  int hwidth, /* Half width of kernel. Full is 2*hwidth+1.
                 Set to <0 to make routine set appropriate value */
  int iap, /* Apodization method. Always ap=0 for d>nap*cdown.
              iap=0 means no apodization ap=1
              iap=1 uses parabola ap=1-(d/(nap*cdown))^2
              iap=2 uses sinc ap=sinc(d/(nap*cdown))
              iap=3 uses Airy with first zero at nap'th zero of main
              all other cases give ap=1 (not guaranteed) */
  int nap, /* Apodizes to nap'th zero */
  int nsub, // Distance between sampled points
  int nxin, // Input size
  int nyin // Input size
);

int init_fresize_user(
  struct fresize_struct *pars,
  int hwidth, // Half width of kernel. Full is 2*hwidth+1.
  int nsub, // Distance between sampled points
  float *user_ker // User specified kernel to convolve with.
                  // Must be of size (2*hwidth+1) x (2*hwidth+1).
                  // Kernel need not be and will not be normalized.
);

int free_fresize(
  struct fresize_struct *pars
);

int fresize(
  struct fresize_struct *pars, // Must have been initialized by init_fresize_XXX
  float *image_in,
  float *image_out,
  int nxin, // Size of input image
  int nyin, // Size of input image
  int nleadin, // Leading dimension of input image. nleadin>=nxin
  int nxout, // Size of xin, yin and image_out
  int nyout, // Size of xin, yin and image_out
  int nleadout, // Leading dimension. nlead>=nx
  int xoff, // Offset in x direction
  int yoff, // Offset in y direction
  float fillval // Value to use if outside area
);

int fsample(
  float *image_in,
  float *image_out,
  int nxin,
  int nyin,
  int nleadin,
  int nxout,
  int nyout,
  int nleadout,
  int nsub,
  int xoff,
  int yoff,
  float fillval
);

int fbin(
  float *image_in,
  float *image_out,
  int nxin,
  int nyin,
  int nleadin,
  int nxout,
  int nyout,
  int nleadout,
  int nsub,
  int xoff,
  int yoff,
  float fillval
);

#endif
Revision:	1.3
Committed:	Fri Jan 22 03:47:55 2016 UTC (7 years, 8 months ago) by arta
Content type:	text/plain
Branch:	MAIN
CVS Tags:	Ver_9-1, Ver_LATEST, Ver_9-3, Ver_9-41, Ver_9-2, Ver_9-5, Ver_9-4, Ver_8-11, Ver_8-12, Ver_9-0, HEAD
Changes since 1.2:	+3 -0 lines
Log Message:	Add Keh-Chengs lev1_dcon program to CVS. Had to do a bunch of clean-up to get rid of multiple definitions in proj/lev0 and proj/lev1.5_hmi.
#	Content
1	#ifndef __FRESIZE_H
2	#define __FRESIZE_H
3	#include <complex.h>
4	#include <fftw3.h>
5
6	struct fresize_struct {
7	int method;
8	int nsub;
9	int hwidth;
10	float ker,kerx,*kery;
11	fftwf_complex helpc,fkernel,*fkernely;
12	float helpin,helpout;
13	fftwf_plan plan1,plan2,plan1y,plan2y;
14	int nxin,nyin,nxinp,nyinp;
15	};
16
17	int init_fresize_sample(
18	struct fresize_struct *pars,
19	int nsub // Distance between sampled points
20	);
21
22	int init_fresize_bin(
23	struct fresize_struct *pars,
24	int nsub // Binsize
25	);
26
27	int init_fresize_boxcar(
28	struct fresize_struct *pars,
29	int hwidth, // Half width of boxcar. Full is 2*hwidth+1.
30	int nsub // Distance between sampled points
31	);
32
33	int init_fresize_boxcar_fft(
34	struct fresize_struct *pars,
35	int hwidth, // Half width of boxcar. Full is 2*hwidth+1.
36	int nsub, // Distance between sampled points
37	int nxin, // Array size
38	int nyin // Array size
39	);
40
41	int init_fresize_gaussian( // Simple square truncated Gaussian
42	struct fresize_struct *pars,
43	float sigma, // Shape is exp(-(d/sigma)^2/2)
44	int hwidth, // Half (truncation) width of kernel. Full is 2*hwidth+1.
45	int nsub // Distance between sampled points
46	);
47
48	int init_fresize_gaussian_fft( // Simple square truncated Gaussian. FFT version.
49	struct fresize_struct *pars,
50	float sigma, // Shape is exp(-(d/sigma)^2/2)
51	int hwidth, // Half (truncation) width of kernel. Full is 2*hwidth+1.
52	int nsub, // Distance between sampled points
53	int nxin, // Array size
54	int nyin // Array size
55	);
56
57	int init_fresize_sinc( // Sinc filter
58	struct fresize_struct *pars,
59	float wsinc, /* Shape is sinc(d/wsinc)*ap(d)
60	wsinc is the amount by which the Nyquist is reduced.
61	May want wsinc=nsub. */
62	int hwidth, // Half width of kernel. Full is 2*hwidth+1.
63	int iap, /* Apodization method. Always ap=0 for d>nap*wsinc.
64	iap=0 means no apodization ap=1
65	iap=1 uses parabola ap=1-(d/(nap*wsinc))^2
66	iap=2 uses sinc ap=sinc(d/(nap*wsinc))
67	all other cases give ap=1 (not guaranteed) */
68	int nap, /* Sinc apodization width in units of wsinc.
69	Normally hwidth=nap*wsinc,
70	but hwidth=napwsinc-1 works for integer /
71	int nsub // Distance between sampled points
72	);
73
74	int init_fresize_sinc_fft( // Sinc filter. FFT version.
75	struct fresize_struct *pars,
76	float wsinc, /* Shape is sinc(d/wsinc)*ap(d)
77	wsinc is the amount by which the Nyquist is reduced.
78	May want wsinc=nsub. */
79	int hwidth, // Half width of kernel. Full is 2*hwidth+1.
80	int iap, /* Apodization method. Always ap=0 for d>nap*wsinc.
81	iap=0 means no apodization ap=1
82	iap=1 uses parabola ap=1-(d/(nap*wsinc))^2
83	iap=2 uses sinc ap=sinc(d/(nap*wsinc))
84	all other cases give ap=1 (not guaranteed) */
85	int nap, /* Sinc apodization width in units of wsinc.
86	Normally hwidth=nap*wsinc,
87	but hwidth=napwsinc-1 works for integer /
88	int nsub, // Distance between sampled points
89	int nxin, // Array size
90	int nyin // Array size
91	);
92
93	int init_fresize_gaussian2( // Circularly truncated Gaussian
94	struct fresize_struct *pars,
95	float sigma, // Shape is exp(-(d/sigma)^2/2)
96	float rmax, // Truncation radius. Probably rmax<=hwidth.
97	int hwidth, // Half (truncation) width of kernel. Full is 2*hwidth+1.
98	int nsub // Distance between sampled points
99	);
100
101	int init_fresize_gaussian2_fft( // Circularly truncated Gaussian. FFT version.
102	struct fresize_struct *pars,
103	float sigma, // Shape is exp(-(d/sigma)^2/2)
104	float rmax, // Truncation radius. Probably rmax<=hwidth.
105	int hwidth, // Half (truncation) width of kernel. Full is 2*hwidth+1.
106	int nsub, // Distance between sampled points
107	int nxin, // Input size
108	int nyin // Input size
109	);
110
111	int init_fresize_airy( // 2D Airy filter
112	struct fresize_struct *pars,
113	float cdown, /* cdown is the amount by which the Nyquist is reduced. */
114	int hwidth, /* Half width of kernel. Full is 2*hwidth+1.
115	Set to <0 to make routine set appropriate value */
116	int iap, /* Apodization method. Always ap=0 for d>Z_nap, where Z_nap os
117	the position of the nap'th zero.
118	iap=0 means no apodization ap=1
119	iap=1 uses parabola ap=1-(d/Z_nap)^2
120	iap=2 uses sinc ap=sinc(d/(Z_nap))
121	iap=3 uses Airy with first zero at Z_nap
122	all other cases give ap=1 (not guaranteed) */
123	int nap, /* Apodizes to nap'th zero */
124	int nsub // Distance between sampled points
125	);
126
127	int init_fresize_airy_fft( // 2D Airy filter. FFT version.
128	struct fresize_struct *pars,
129	float cdown, /* cdown is the amount by which the Nyquist is reduced. */
130	int hwidth, /* Half width of kernel. Full is 2*hwidth+1.
131	Set to <0 to make routine set appropriate value */
132	int iap, /* Apodization method. Always ap=0 for d>nap*cdown.
133	iap=0 means no apodization ap=1
134	iap=1 uses parabola ap=1-(d/(nap*cdown))^2
135	iap=2 uses sinc ap=sinc(d/(nap*cdown))
136	iap=3 uses Airy with first zero at nap'th zero of main
137	all other cases give ap=1 (not guaranteed) */
138	int nap, /* Apodizes to nap'th zero */
139	int nsub, // Distance between sampled points
140	int nxin, // Input size
141	int nyin // Input size
142	);
143
144	int init_fresize_user(
145	struct fresize_struct *pars,
146	int hwidth, // Half width of kernel. Full is 2*hwidth+1.
147	int nsub, // Distance between sampled points
148	float *user_ker // User specified kernel to convolve with.
149	// Must be of size (2hwidth+1) x (2hwidth+1).
150	// Kernel need not be and will not be normalized.
151	);
152
153	int free_fresize(
154	struct fresize_struct *pars
155	);
156
157	int fresize(
158	struct fresize_struct *pars, // Must have been initialized by init_fresize_XXX
159	float *image_in,
160	float *image_out,
161	int nxin, // Size of input image
162	int nyin, // Size of input image
163	int nleadin, // Leading dimension of input image. nleadin>=nxin
164	int nxout, // Size of xin, yin and image_out
165	int nyout, // Size of xin, yin and image_out
166	int nleadout, // Leading dimension. nlead>=nx
167	int xoff, // Offset in x direction
168	int yoff, // Offset in y direction
169	float fillval // Value to use if outside area
170	);
171
172	int fsample(
173	float *image_in,
174	float *image_out,
175	int nxin,
176	int nyin,
177	int nleadin,
178	int nxout,
179	int nyout,
180	int nleadout,
181	int nsub,
182	int xoff,
183	int yoff,
184	float fillval
185	);
186
187	int fbin(
188	float *image_in,
189	float *image_out,
190	int nxin,
191	int nyin,
192	int nleadin,
193	int nxout,
194	int nyout,
195	int nleadout,
196	int nsub,
197	int xoff,
198	int yoff,
199	float fillval
200	);
201
202	#endif