韦伯尔分布拟合海杂波

最新推荐文章于 2023-07-08 07:12:35 发布

xiexienN

最新推荐文章于 2023-07-08 07:12:35 发布

阅读量773

点赞数

本文链接：https://blog.csdn.net/xiexienN/article/details/106972901

版权

clear
img=imread("SARShip-1.0-31.tif");
[m,n]=size(img);
img=im2uint8(img);
index=find(img==0);
img(index)=1;
[Freq,X]=hist(double(img(:)),180);
binWidth=X(2)-X(1);
bar(X,Freq/binWidth/sum(Freq));
% axis([0,65536,0,4*1e5]
img=double(img);
u1=mean(img,'all');
sigma21=mean(img.*img,'all')-u1*u1;
u=mean(log(img),'all');
sigma2=mean(log(img).^2,'all')-(mean(log(img),'all'))^2;
x=1:256;
f=1./(x.*sqrt(sigma2*2*pi)).*exp(-(log(x)-u).^2/(2*sigma2));
hold on;
plot(x,f,'LineWidth',2);
f1=1/(sqrt(2*pi*sigma21))*exp(-0.5*(x-u1).^2./(2*sigma21));
hold on;
plot(f1,'LineWidth',2);
mlog2=mean(log(img).*log(img),'all');
mlog=mean(log(img),'all');
p=(6/(pi*pi)*(m*n)/(m*n-1)*(mlog2-mlog^2))^(-0.5);
q=exp(mlog+0.5772*p^(-1));
x=0.1:256;
f=p/q*(x./q).^(p-1).*exp(-(x./q).^p);
hold on;
plot(x,f,'LineWidth',2);
sum(f1)
legend('直方图','对数正态','正态','韦布尔');