├── OCTSimulationBooks ├── 827B3d01.pdf └── MasterThesis_PatriciaCarvalho.pdf ├── README.md ├── interference.asv ├── interference.m ├── rough.asv └── rough.m /OCTSimulationBooks/827B3d01.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/ajaygunalan/octSimulation/a28944e2a807aca0ea0918aca8a4b8a152456a6d/OCTSimulationBooks/827B3d01.pdf -------------------------------------------------------------------------------- /OCTSimulationBooks/MasterThesis_PatriciaCarvalho.pdf: -------------------------------------------------------------------------------- https://raw.githubusercontent.com/ajaygunalan/octSimulation/a28944e2a807aca0ea0918aca8a4b8a152456a6d/OCTSimulationBooks/MasterThesis_PatriciaCarvalho.pdf -------------------------------------------------------------------------------- /README.md: -------------------------------------------------------------------------------- 1 | # oct 2 | Simulation of Optical Coherence Tomography (OCT) 3 | 4 | 23 Dec 2020 5 | 6 | ## References 7 | 8 | * [Optical-Coherence-Tomography-using-Matlab](https://github.com/shbkd/Optical-Coherence-Tomography-using-Matlab) 9 | * [Optical Coherence Tomography – System and Simulation (Zemax)](https://www.youtube.com/watch?v=MvCzGLSMqNM) 10 | * [OPTICAL COHERENCE TOMOGRAPHY Layout Simulation using MATLAB](https://estudogeral.sib.uc.pt/bitstream/10316/31107/1/MasterThesis_PatriciaCarvalho.pdf) 11 | * [DTU Ph.D. thesis on Optical Coherence Tomography: Modeling and Applications](https://www.google.com/url?sa=t&rct=j&q=&esrc=s&source=web&cd=&ved=2ahUKEwjkqazyteTtAhXD-6QKHdQ4CPAQFjACegQIAxAC&url=https%3A%2F%2Forbit.dtu.dk%2Ffiles%2F7810773%2F827B3d01.pdf&usg=AOvVaw0M_E68vavnHRyagOEP9NwA) 12 | * [OCT System design (Matlab Simulation)](https://www.youtube.com/watch?v=6JgmXwkjkmE&list=PLAJchV67pirTQqxZUYJv_PHDE4muOOTbn) 13 | * [Michelson Interferometer Fringe pattern](https://in.mathworks.com/matlabcentral/fileexchange/73482-michelson-interferometer-fringe-pattern) 14 | * [JOptics](http://www.ub.edu/javaoptics/index-en.html) 15 | * [Best OCT explanation](https://www.youtube.com/playlist?list=PL-dIBMwXD0RXVArXRGrRyPcUZSF4BNSCp) 16 | * [Dark Filed & Phase Contrast](https://www.youtube.com/watch?v=14oH3jMdCFM&feature=emb_title) 17 | * [Basics of Microscopy (Dark Field, Phase Contrast, DIC & DHM)- Part 2 by Dr Kedar Khare, IIT-Delhi](https://www.youtube.com/watch?v=mETbbR_ciR8&feature=emb_title) 18 | -------------------------------------------------------------------------------- /interference.asv: -------------------------------------------------------------------------------- 1 | clc 2 | clear 3 | 4 | no_elemen = 100; 5 | x = -0.5:1/(N-1):0.5; -------------------------------------------------------------------------------- /interference.m: -------------------------------------------------------------------------------- 1 | clc 2 | clear 3 | 4 | no_elements = 100; 5 | x = -0.5:1/(no_elements-1):0.5; 6 | y = x; 7 | 8 | lambda = 0.633; %wavelength of red laser in micrometer 9 | k = 2*pi/lambda; 10 | theta = pi/3; 11 | 12 | [X, Y] = meshgrid(x, y); 13 | E1 = exp(i*0); 14 | E2 = exp(i*k*x*cos(theta)); 15 | 16 | E = (E1 + E2).*conj(E1+E2); 17 | 18 | figure(); 19 | imagesc(E); 20 | %colormap gray; 21 | colorbar; -------------------------------------------------------------------------------- /rough.asv: -------------------------------------------------------------------------------- 1 | 2 | x = -1:0.5:1; 3 | y = 1:0.5:1; 4 | [X, Y] = meshgrid(x, y); 5 | plot -------------------------------------------------------------------------------- /rough.m: -------------------------------------------------------------------------------- 1 | % constants 2 | c = 3.00e+08; 3 | % gaussian source 4 | lc = 800e-09; 5 | l = linspace(lc-(lc/3), lc+(lc/3), 1024); 6 | fwhm = 50e-09; 7 | sigma = fwhm./sqrt(8*log(2)); 8 | wc = (2*pi*c)/lc; 9 | w = linspace(wc-(wc/3), wc+(wc/3), 1024); 10 | fwhm_w = (2*pi*c*fwhm)./(lc.^2); 11 | sigma2 = fwhm_w./sqrt(8*log(2)); 12 | y = gaussmf(w, [sigma2, wc]); 13 | figure (), plot (w, y), axis tight, title ('Light Source Spectrum') 14 | xlabel ('Optical Frequency (rad.s-1)'), ylabel ('Amplitude (a.u.)') 15 | % sample 16 | n = [1.00 1.30 1.50 1.00]; 17 | z = [5.00e-06 15.00e-06 30.00e-06 0.00e-06]; 18 | s1 = 0; 19 | h = 0; 20 | for i = 1:3 21 | rj = (n(i+1)-n(i))/(n(i+1)+n(i)); 22 | s1 = s1 + n(i)*z(i); 23 | h = h + rj*exp(1i*2.*(w./c)*s1); 24 | end 25 | % time domain 26 | x = linspace (0, 100e-6, 1024); 27 | T1 = zeros(size(x)); 28 | for j = 1:length(x) 29 | for jj = 1:length(w) 30 | ph = cos(2*x(j)*w(jj)/c); 31 | T1(j) = T1(j) + real(0.5*(y(jj)*h(jj)*ph)); 32 | end 33 | end 34 | 132 35 | figure (), plot(x./1e-6, T1), axis tight, title ('TD-OCT Interferogram') 36 | xlabel ('Mirror displacement (um)'), ylabel ('Amplitude (a.u.)') 37 | --------------------------------------------------------------------------------