interpol - interpolate (12-point 2-D cubic) the field into new size and dimension.

Calling Sequence

[field_out]=interpol(field_in,size_new[,new_number[,x_shift[,y_shift[,angle[,magnif]]]]])

Parameters

• field_in :input field. it must be diffractive wavefront.
• size_new : the new size after transformt.if size_new=0, then the size of the field will not changed.
• new_number : input the new number after transformt(integer). default is 0, if new_number=0, then the size of the field will not changed.
• x_shift : the shifts of the X direction in metre. ( default is 0.)
• y_shift : the shifts of the Y direction in metre. ( default is 0.)
• angle : the rotated angle in degree.( default is 0. )
• magnif : the coefficient of magnification ( default is 1.).
• field_out : the output field or diffractive wavefront.

Description

Interpol is the tool to manipulate the size and the dimension of the grid and to change the shift, rotation and the scale of the field distribution. besides the input field It accepts six command line arguments, the first is the new size of the grid, the second and others are optional. The second argument gives the new sampling, the third and fourth determine shifts in X and Y directions,the fifth gives the field rotation (after shifting). The last argument defines the magnification, its action is equivalent to passing the beam through a focal system with magnification M (without diffraction, but preserving the integral intensity).

For example if the field was propagated with FFT algorithm forvard, then the grid contains empty borders which are not necessary for further processing. Other way around, after forward we have to add some empty borders to continue with forvard. Interpol is useful for interpolating into grids with different size and samplings. Of course, it is not wise to interpolate a grid of 512x512 points into a grid of 4x4, and then back to 512x512, because all information about the initial field will be lost. The same is true for interpolating a distribution defined in the square of 1x1m to 1x1mm and back. When interpolating into a grid with larger size, for example from 1x1 to 2x2, the program puts zeros into added regions.

Examples of transformations (we suppose the input field is f):

interpol (f,1) - interpolate to a new grid with size 1

interpol (f,1,128) - interpolate to a new grid with size 1 and number of points 128

interpol (f,0) - don't change anything

interpol (f,1,128,0.2) - interpolate to a new grid with size 1 and number of points 128, then shift in X direction by 0.2

interpol (f,1,128,0.2,0.3) - interpolate to a new grid with size 1 and number of points 128, then shift in X direction by 0.2 and in Y direction by 0.3

interpol (f,1,128,0.2,0.3,23) - interpolate to a new grid with size 1 and number of points 128, then and shift in X direction by 0.2 and in Y direction by 0.3, then rotate by .

interpol (f,0,0,128,0.2,0.3,23,1.3) - interpolate into a new grid with size and sampling both inherited from the previous filter, shift in X direction by 0.2 and in Y direction by 0.3, then rotate and then magnify the field with a ratio of 1.3

interpol (f,0,0,0.,0.,0.,1.3) - inherit size and sampling, magnify (afocal) the field with a ratio of 1.3

Examples

	f1=begin(0.01,1e-6);
	f2=rect_ap(f1,0.005);
	f3=forvard(f2,0.2);
	file_ps(f3,"int1");
	f4=interpol(f3,0.0075,64);
	field_plot(f4);
 

Author

Chen jingyuan, Gan guangyong (IAPCM, Beijing,China.)

See Also

lp_interp ,   interpolate ,