function [F]=Forcollect(H,W,POSP,th1,FS)
COLH=[]
COLW=[];
j=1;
for i=1:1:length(POSP)
[I,J]= find((W-POSP(i))>0)
if(I>th1)
COLH{j}=H(I(1)-2*th1:1:I(1)+2*th1);
j=j+1;
COLW=[COLW W(I(1)-1)];
end
end
R=[];
F=[];
for i=1:1:length(COLH)    
[M,N]=max(COLH{i});
if ((N>th1)&(COLH{i}(N)>0.1))    
%R=[R ((COLH{i}(N)/COLH{i}(N-th1))+(COLH{i}(N)-COLH{i}(N+th1)))/2];
temp=((COLH{i}(N)/COLH{i}(N-th1))+(COLH{i}(N)/COLH{i}(N+th1)))/2 ;
if(temp<1.414)
    F=[F (COLW(i)*FS)/(2*pi)];
else
    F=[F 0];
end
else
if(N>1)    
 temp=((COLH{i}(N)/COLH{i}(N-1))+(COLH{i}(N)/COLH{i}(N+1)))/2 ;
if(temp<(1.414/th1))
    F=[F (COLW(i)*FS)/(2*pi)];
else
    F=[F 0];
end     
else
    F=[F 0];
end
end
end
