%修正p值gm(1,1) 检测
clc;clear;
m2=imread('lena_gray.jpg');
m=double(m2)./255;
mm2=m;
p=0.5;
[n1,n2]=size(m);
for i=1:n1-1 %水平方向
for j=1:n2-2
AA2=[m(i,j);m(i,j)+m(i,j+1);m(i,j)+m(i,j+1)+m(i,j+2);m(i,j)+m(i,j+1)+m(i,j+2)+m(i+1,j)];
%A2=[AA2(1),AA2(1)+AA2(2),AA2(1)+AA2(2)+AA2(3),AA2(1)+AA2(2)+AA2(3)+AA2(4)];
B2=[-0.5*(AA2(1)+AA2(2)),1;-0.5*(AA2(2)+AA2(3)),1;-0.5*(AA2(3)+AA2(4)),1];
F2=[m(i,j+1);m(i,j+2);m(i+1,j)]; %AA2为累加过的,这里要不累加的
E2=inv(B2'*B2)*B2'*F2; %E2(1)会变得非常小
% G2=(m(i,j)-E2(2)/E2(1))*exp(-4*E2(1)/p)+E2(2)/E2(1);
% T2=(m(i,j)-E2(2)/E2(1))*exp(-3*E2(1)/p)+E2(2)/E2(1);
% S2=(m(i,j)-E2(2)/E2(1))*exp(-2*E2(1)/p)+E2(2)/E2(1);
% mm2(i+1,j+1)=G2-2*T2+S2;
mm2(i+1,j+1) = ( (1-exp(E2(1)*p)) * (m(i,j) - E2(2)/E2(1)) * exp(-E2(1)*4/p) ); %此为化简后的式子,减少不必要的计算
% R2=uint8(round(mm2*255)); %可以提取到循环外部
end
end
R2=uint8(round(mm2*255));
imwrite(R2,'xiesong2.jpg');figure;imshow(R2);
fangli.rar_GM(1_1)
版权申诉
16 浏览量
2022-09-22
18:36:23
上传
评论
收藏 642B RAR 举报
四散
- 粉丝: 54
- 资源: 1万+
最新资源
- nursery.data
- EasyCodeConfig-Final
- whisper-small-ct2 OPENAI ct2 语言识别模型实例,python环境请自行下载
- FPGA项目verilog实例资料带项目源码设计文档基于FPGA的序列检测器
- Allegro Free viewer
- 独库公路 kmz文件(含沿途主要地点)
- ROS机器人操作系统Robot Operating System
- 170700352842914FusionApp_2.0.0-beta4.apk
- Python高级的、通用的、解释型的编程语言
- FPGA项目verilog实例资料带项目源码设计文档基于FPGA的频率计与串口通信
资源上传下载、课程学习等过程中有任何疑问或建议,欢迎提出宝贵意见哦~我们会及时处理!
点击此处反馈