SNR-提取波的初至时

能够精确地提取雷达波的初至时，进行层析成像-Radar can accurately extract the first to when tomography

function speed=SNR（speed）
%信噪比法
% clear;clc
% [datpathname]= uigetfile （‘*.out‘ ‘Load Ramac format radar file‘）;
% [HeaderFields]=gprmax（[pathnamedat]）;

% dt=1.0e9*Header.dt;
% nlen=Header.iterations;
% nt=Header.nrx;
% dx=1;
% 
% DAT13=Fields.ez‘;
% DAT=DAT13‘+mean（mean（DAT13‘））;
% 
% am1=0.25*max（max（abs（DAT13）））;
% for jj=1:nt
% DAT13（:jj）=（1+jj*0.3）*DAT13（:jj）;
% end
A=speed.A+mean（mean（speed.A））;


m=length（A（1:））;
n=length（A（:1））;
R（1:m）=0;
for i=1:m
    d=0;p=0;
    r=0;
    for j=11:n-11
        for t=j-5:j
            temp1（t）=（A（ti）.^2）*（exp（-0.5*（t-j））.^2）;
            d=d+temp1（t）;
        end
        for t=j:j+5
            temp2（t）=（A（ti）.^2）*（exp（-0.5*（j-t））.^2）;
            p=p+temp2（t）;
        end
        r（j）=（p-d）./d;
    end
   
    y=max（r）;
    k=find（r==y）;
    if（size（k）==1）
        R（i）=k;
    else
        R（i）=k（1）;
    end
end
figure

imagesc（A）;
% for ii=1:nt
%     sig（1:nlen）=A13（ii1:nlen）;
%     am（ii）=2*max（sig）;
%     h=plot（（ii）*dx+sig（1:nlen）/（am1*3）[1:nlen]*dt‘k-‘）;hold on;
%     set（h‘Linewidth‘1.5）
% end
% hold off
axis ij
axis tight
hold on;
% plot（R*dt）;
plot（R）;
% colorbar;
title（‘信噪比最大法‘）
xlabel（‘Trace number‘）;
ylabel（‘time [ns]‘）;
speed.Cmin=0;
speed.Cmax=0;
speed.amp=0;
p=length（R（:1））;
q=length（R（1:））;
if p>q
%     speed.time=C.*speed.dt+163.4;
    speed.time=R.*speed.dt;
else
%     speed.time=C‘.*speed.dt+163.4;
    speed.time=R‘.*speed.dt;
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1518  2010-04-08 09:47  SNR.m

----------- ---------  ---------- -----  ----

                 1518                    1