reconstructor - construct a wavefront reconstructor.
construct a least squares reconstructors. it aims to represent reconstructors that yield zonal residuals as well as some number of zernike mode ( here tip and tilt ) residuals.
this function construct a least squares reconstructor from an aperture, an ideal deformable mirror, and a square lenslet array. The influence function of the ideal deformable mirror is assumed to be pyramidal.
NOTE - currently only Fried geometry is supported
The subaperture illumination threshold argument to this function allows one to zero geometry matrix elements corresponding to actuator influence functions and subapertures that overlap by less than the threshold value. This value should lie between zero and one.
The eigenvalue threshold argument to this function allows one to remove modes from the reconstructor with eigenvalues lower than the threshold. This argument is interpreted as the ratio between the threshold eigenvalue and the largest eigenvalue, and its value should lie between zero and one.
nsubaps=16; areal_threshold=0.5; eigenvalue_threshold=1e-6; lenslet_x_axes=16;lenslet_y_axes=16; focal_length = .012246;lenslet_pitch = .000252; final_wf_pix_per_lenslet = 32;final_wf_pix_per_xform = 32; LA=lnslt_array(lenslet_x_axes,lenslet_y_axes,focal_length,lenslet_pitch,final_wf_pix_per_lenslet,final_wf_pix_per_xform); D = 16*lenslet_pitch; AP=circle(D); pitch = D/nsubaps; actuator_velocity = 100; dm_AP=circle(2*D); act_dim=nsubaps+1; act=actuator_array(act_dim,act_dim,pitch,actuator_velocity); dm=ideal_dm(dm_AP,act); projection_modes=znk_mod(1); projection_modes=set_znk_cos_coef(projection_modes,1,1,1); projection_modes=set_znk_sin_coef(projection_modes,1,1,1); Reconstructor= reconstructor(AP,dm,LA,projection_modes,areal_threshold,eigenvalue_threshold); reconstructor_fits(Reconstructor,"reconstructor");
Chen jingyuan, Gan guangyong (IAPCM, Beijing,China.)
reconstructor_fits , zernike_residuals , zonal_residuals , reconstructor ,