Team:USTC/Software

In this section, we will introduce all the software program we used when processing our modeling and experimental data. Through our modeling, many complicated processdures are simply omitted by our convenient programs.

In general, our software programs are used in several sections,

  • Fringes Analysis, analysis of interference pattern.
  • Adhesion Dynamics, a series of programs relating to bacteria adhesion dynamics.
  • Rose Prediction, which calculate the advantages of our ROSE improvement.

All software of our project, you can download them at Github:2015USTCiGEM. All codes are based on

Fringes Pattern Simulation-Film I

This program is used to simulate fringes pattern deliverd by film I.

Code:

function drawndh(r,h,t,n) %r is radius of film,h is the distance of deformation,t is the drift angle of film,n is pixel number wl=6.5e-7; %wave length R=r.^2/(2*h); [x,y]=meshgrid(linspace(-0.004,0.004,n),linspace(-0.004,0.004,n)); r2=x.^2+y.^2; delta=pi*r2/(R*wl)-2*pi*y*sin(t)/wl; I=abs(sin(delta).^2); imshow(I); end

Demo
Figure 1: Film II deformation simulation result

After the formation of demo, we are able to recongize the real situation in our pre-experiment.

Fringes Pattern Simulation-Film II

This program is used to simulate fringes pattern deliverd by film II.

Code

function drawfr(a,b,h,n)%a is the angle between film and x axis, b is the angle between film and y axis, h is the distance of film deformation, n is the pixel number
[X,Y]=meshgrid(linspace(-0.003,0.003,n),linspace(-0.003,0.003,n));
wl=6.5e-7;
L=2*(X.*tan(0.01)+Y.*tan(0.01));
delta=2*pi*L/wl;
I=1+cos(delta);
imshow(I)
end

Demo
Figure 2 :Film II deformation simulation result

Film I Fringes Analysis

Using this code, we are able to capture some important parameters in Film I.

Code

function ms3=solvendh(X,nx,ny,n,dx,dy)%nx is the number of fringes on x axis,ny is the number of fringes on y axis,n is the number of fringes on diag,dx is the x length of CCD,dx is the x length of CCD.
x=X(1);
y=X(2);
z=X(3);
l=1.2e-4;
ms3(1)=(2*(z+nx)*l)^0.5-(2*z*l)^0.5-((x+dx)^2+y^2)^0.5+(x^2+y^2)^0.5;
ms3(2)=(2*(z+ny)*l)^0.5-(2*z*l)^0.5-(x^2+(y+dy)^2)^0.5+(x^2+y^2)^0.5;
ms3(3)=(2*(z+n)*l)^0.5-(2*z*l)^0.5-((x+dx)^2+(y+dy)^2)^0.5+(x^2+y^2)^0.5;

Demo
Figure 3 : Film I Fringes analysis running result

Film II Fringes Analysis-Fringes number analysis

This is the final program used to get our NDM calibration. Based on this software program, we are able to count the number of bright fringes. As a matter of fact, the variation of bright fringes is the key breakthrough to antibiotic analysis.

Code

function seeim(i)
name=num2str(i);
I=imread(name,'pnm');
G=rgb2gray(I);
imshow(G);

Demo
Figure 4: Checking figure program

function A=getfr(i)
name=num2str(i);
I=imread(name,'pnm');
G=rgb2gray(I);
G1=medfilt2(G);
[x,y]=size(G1);
r=zeros(x,1);
for m=1:20
    p=G1(:,x-m+1);
    p=double(p);
    r=r+p;
end
[a,b]=size(r);
for m=1:b-1
    r(m,1)=(r(m-1,1)+r(m,1)+r(m+1,1));
end
A=r;
function M=getfrs(i)
A=getfr(i);
t=1:480;
[pksa,locsa]=findpeaks(A,'minpeakdistance',3);plot(t,A,'b',locsa,pksa,'bo');
[a,k]=size(pksa);
M=a;

Running result
ans =
78

Figure 5: Fringes counting program running result

To get a series of results, this program is required.

function N=getpkss(i,j)
B=zeros(j-i+1,1);
for m=i:j
    A=getfr(m);
    [pksa,locsa]=findpeaks(A,'minpeakdistance',3);
    [a,k]=size(pksa);
    B(m-i+1,1)=a;
end
N=B;

Demo

ans =

78
80
77
76
75
79
81
80
75
80
82
82
81
78
79
82
81
81
83
81
79

During analysis of ahdesion assay results, we are quite annoyed about counting total bacteria number and moving bacteria number, which is really important for mechanism analysis and results determination. Consequently, this part explains how we efficiently accomplished these tough work based on our software.

Bacteria counting program

This program is used to count the total number of bacteria.

Code

function a=shujun(i)
name=num2str(i);
I=imread(name,'jpg');
L=graythresh(I);%get the self-adapting threshold
G=im2bw(I,0.3);%image binaryzation with special threshold
W1=~G;
W2=bwmorph(W1,'majority',20);%mathematical morphology operations
W3=medfilt2(W2,[10,10]);%filtering
W4=bwareaopen(W3,40);%delete small area to reduce error noises
[labeled,num]=bwlabel(W4,4);%counting number of objects
a=num;
end

Demo

input image:
Figure 6: Adhesion Assay image
output result
ans =
865

function A=mulshu(m)
A=zeros(m,1);
for i=1:m
    A(i,1)=shujun(i);
end

Moving Bacteria Counting Program

This program is used to count the total number of moving bacteria. Along with bacterial counting program, we finally got the mechanism of polylysine interaction.

Code

function b=act(m,n)
s1=num2str(m);
s2=num2str(n);
I1=imread(s1,'jpg');
I2=imread(s2,'jpg');
G=I1-I2;
b=shujun(G);
end
function A=mulact(m)
A=zeros(m,1);
for i=1:m
    name1=num2str(i);
    name2=num2str(i+1);
    I1=imread(name1,'jpg');
    I2=imread(name2,'jpg');
    c=act(I1,I2);
    A(i,1)=c;
end

Demo

ans =
405

Bacteria total number-time curve

Coding

function drawBT1(Ka,Kd,C,Vz,sigma0)%Ka is the adhesive rate of bacteria, Kd is the drop rate of bacteria, C is the density of bacteria solution,Vz is the velocity of bacteria,sigma0 is the maximum density of baacteria on surface
t=linspace(0,100,101);
K=Ka*C*Vz/(Kd*sigma0+Ka*C*Vz);
sigma=K*sigma0*(1-exp(-Ka*C*Vz*t/(K*sigma0));
plot(t,sigma)
end

Running result
Figure 7: Movement number simulation result

Bacteria movement number-Time curve

Coding

function drawBMT(c,b,m0,k)
t=linspace(0,100,101);
M=c*(1-exp(-b*t))*m0.*exp(-k*t);
plot(t,M)
xlabel('Time(s)');
ylabel('Movement bacteria number');
title('PAO1-PLL-0 Movement number-time');
end

Running result
Figure8:Bacteria movement number-Time simulation result

This program is used to predict the advantages of ROSE construction. We successfully demonstrate our ROSE will be able to amplify fluorescence siganl.

Code

clear all
R=R,k=k,k1=k1,k2=k2,k3=k3,k4=k4,k5=k5,k6=k6,k7=k7,k0=k0; %define the constant
t=linspace(0,10); %detecting time range
F=k1*S;
A=k*dsolve('Dx=k2*F','x(0)=0','t');
RA=k3*R*A;
X=k4*RA;
Cl=dsolve('Dx=k5*X','x(0)=0','t');
lambda=1-k6*cl;
t0=fzero(lambda,t); %find the zero point of lambda
if lambda>=0
    G=dsolve('Dx=k7*lambda','x(0)=0','t');
else
    Gm=G(t0);
end
%finish the equation solving step
plot(G,t); %print the G-t picture
hold on
plot(k0*t,t); %print the traditional strategy result in the same image

Demo

Figure 9: ROSE Dynamics simulation result

Contact Us

University of Science and Technology of China, No.96, JinZhai Road Baohe District,Hefei,Anhui, 230026,P.R.China.

Links