基于matlab的可见光成像通信OOK解调

直接给出代码和视觉效果

clc
clear
close all

s=strcat('frame0010','.jpg'); 
A=imread(s);
Agray=rgb2gray(A);
Agray = Agray(780:870,895:988);
imshow(Agray);
[row,col] = size(Agray);
%%%a=row%%%%%行
%%%b=col%%%%%%列
Agray = double(Agray);
%col = Agray(783:865,923);
a=zeros(row,1);
for i=1:row
    a(i)=sum(Agray(i,:))/col;
end
col=a;
%% LED1
ratio  = 3.9;
q = 0:ratio:length(col)*ratio-ratio;
xi = 0:1:length(col)*ratio;
Yi=interp1(q,col,xi,'spline') ;
col = Yi; % interpolated signal


figure,imshow(uint8(Agray));
ss = 1:9:length(col);
figure(1),plot(col),hold on,plot(ss,col(ss),'*'),hold on,plot(1:length(col),ones(1,length(col))*80,'b--');