asm_reconstruction.m File Reference

Go to the source code of this file.

Functions
	fclose ('all')
	imwrite (I_rec, sprintf('OPH_point_3_%.2fmm.jpg', zz/mm+20))
	title (sprintf('ASM 3 point Reconst %.2f mm', zz/mm))
	saveas (fig, sprintf('ASM 3 point Reconst %.2f mm.jpg', zz/mm))
	pause (0.2)

Variables
	clc
unit	mm = 1e-3
	um = mm*mm
	nm = um*mm
command	asm_kernel = @(f, wvl, x, y, res, pp) exp(-2j*pi*f .* sqrt( wvl^-2 - ((x/pp - 0.5)./pp/res).^2 - ((y/pp - 0.5)./pp/res).^2))
angular spectrum method kernel	norm2int = @(img) uint8(255.*img./max(max(img)))
parameters	h = 2160
	v = 3840
	wvl = 525 * nm
green	pp = 3.6*um
pixel pitch	res = v/2
	ReconstLength = 1
input file name	fim_py = 'File name_IM.bmp'
	fre_py = 'File name_RE.bmp'
Convert to complex image phase part	im_py = imread(fim_py)
crop	phase = (im_py - 0.5) .* (2*pi)
amplitude part	re_py = imread(fre_py)
adapt ASM	r = (-pp*h/2 + pp/2):pp:(pp*h/2 - pp/2)
	c = r
Complex image	ch = real + 1j*imag
	A = fftshift(fft2(fftshift(ch)))
	p = asm_kernel(zz, wvl, C, R, h, pp)
Reconstruction part	figure
for	zz
	Az1 = A .* p
	EI = fftshift(ifft2(fftshift(Az1)))
	I_rec = EI .* conj(EI)
	fig = imagesc(I_rec)

Function Documentation

fclose()

fclose

(

'all'

)

imwrite()

imwrite	(	I_rec	,
		sprintf( 'OPH_point_3_%.2fmm.jpg', zz/mm+20)
	)

pause()

pause

(

0.

2

)

saveas()

saveas	(	fig	,
		sprintf( 'ASM 3 point Reconst %.2f mm.jpg', zz/mm)
	)

title()

title

(

sprintf( 'ASM 3 point Reconst %.2f mm', zz/mm)

)

Variable Documentation

A

A = fftshift(fft2(fftshift(ch)))

Definition at line 45 of file asm_reconstruction.m.

asm_kernel

command asm_kernel = @(f, wvl, x, y, res, pp) exp(-2j*pi*f .* sqrt( wvl^-2 - ((x/pp - 0.5)./pp/res).^2 - ((y/pp - 0.5)./pp/res).^2))

inline

Definition at line 10 of file asm_reconstruction.m.

Az1

Az1 = A .* p

Definition at line 54 of file asm_reconstruction.m.

c

c = r

Definition at line 39 of file asm_reconstruction.m.

ch

ch = real + 1j*imag

Definition at line 43 of file asm_reconstruction.m.

clc

clc

Definition at line 1 of file asm_reconstruction.m.

EI

EI = fftshift(ifft2(fftshift(Az1)))

Definition at line 55 of file asm_reconstruction.m.

fig

fig = imagesc(I_rec)

Definition at line 61 of file asm_reconstruction.m.

figure

Reconstruction part figure

Definition at line 50 of file asm_reconstruction.m.

fim_py

input file name fim_py = 'File name_IM.bmp'

Definition at line 22 of file asm_reconstruction.m.

fre_py

fre_py = 'File name_RE.bmp'

Definition at line 23 of file asm_reconstruction.m.

h

parameters h = 2160

Definition at line 14 of file asm_reconstruction.m.

I_rec

normalize term I_rec = EI .* conj(EI)

Definition at line 56 of file asm_reconstruction.m.

im_py

im_py = imread(fim_py)

Definition at line 27 of file asm_reconstruction.m.

mm

unit mm = 1e-3

Definition at line 5 of file asm_reconstruction.m.

nm

nm = um*mm

Definition at line 7 of file asm_reconstruction.m.

norm2int

angular spectrum method kernel norm2int = @(img) uint8(255.*img./max(max(img)))

Definition at line 11 of file asm_reconstruction.m.

p

p = asm_kernel(zz, wvl, C, R, h, pp)

Definition at line 47 of file asm_reconstruction.m.

phase

crop phase = (im_py - 0.5) .* (2*pi)

Definition at line 30 of file asm_reconstruction.m.

pp

green pp = 3.6*um

Definition at line 16 of file asm_reconstruction.m.

r

adapt ASM r = (-pp*h/2 + pp/2):pp:(pp*h/2 - pp/2)

Definition at line 38 of file asm_reconstruction.m.

re_py

re_py = imread(fre_py)

Definition at line 33 of file asm_reconstruction.m.

ReconstLength

ReconstLength = 1

Definition at line 18 of file asm_reconstruction.m.

res

pixel pitch res = v/2

Definition at line 17 of file asm_reconstruction.m.

um

um = mm*mm

Definition at line 6 of file asm_reconstruction.m.

v

v = 3840

Definition at line 14 of file asm_reconstruction.m.

wvl

wvl = 525 * nm

Definition at line 15 of file asm_reconstruction.m.

zz

for zz

Initial value:

= (400:10:980)*mm
%for zz = ReconstLength * mm
    p = asm_kernel(zz, wvl, C, R, h, pp)

Definition at line 51 of file asm_reconstruction.m.