TA的每日心情 | 开心 2022-1-29 15:04 |
|---|
签到天数: 1 天 [LV.1]初来乍到
|
function main(), X9 N. |. q7 v. \7 W
clc; close all;
+ _) R# a7 s, F. Z' s3 Z# j% r. LI=imread('cameraman.tif'); %用的是原图8 ~" g( E) T0 Y$ P
I1=imnoise(I,'gaussian',0,0.01);
: s& c. n) [- f' t3 X5 c5 s. nfigure;
4 q+ i& z8 v$ z! s! Q- L' asubplot(1, 2, 1); imshow(I); title('原图像');$ i9 [7 G; a( f i0 L0 K+ U
subplot(1, 2, 2); imshow(I1); title('加噪声图像'); N L0 q, T/ |3 H
[PSNR, MSE] = psnr(I, I1)7 e) ^! S% A; E3 C9 _8 `
% 图像峰值信噪比PSNR、均方根误差MSE, ~6 Z4 T8 ^5 H7 z" X, d
function [PSNR, MSE] = psnr(X, Y)
5 a6 }4 H& d# k) u' e. _% 计算峰值信噪比PSNR、均方根误差MSE
I0 A0 e5 v2 J# h) i! v5 W% 如果输入Y为空,则视为X与其本身来计算PSNR、MSE
* i9 M: u( a- O% K5 k5 \5 x
' h* M/ F1 `3 ]1 W" n/ Bif nargin<2/ \1 [+ H/ F, T5 l
D = X;
- L8 Z' C/ { @, N. |5 N! m) e( Selse
3 P9 ~% o! G* ]4 ?# [. {; p if any(size(X)~=size(Y))% M( {0 k! o7 @) s& d; {
error('The input size is not equal to each other!');# P' f8 ?9 o' o7 ]2 ~" i
end
2 g$ g4 O; }4 D7 ~8 E3 [ D = X-Y;
& i; q2 p! v6 a2 dend
1 a2 L5 _8 ]2 X2 ~$ W' [MSE = sum(D(: ).*D(: ))/prod(size(X));
& j& V/ ~" G+ aPSNR = 10*log10(255^2/MSE);% C: d; G8 D8 }& _* [
function main()
' J1 y y: M/ ~. I% g! P( oclc; close all;7 O$ a$ Z. v. V7 A. `2 q* B9 E
I=imread('cameraman.tif'); %用的是原图
- A8 o& s; o& r: M4 j/ PI1=imnoise(I,'gaussian',0,0.01); s$ E9 k: ?/ ^: D4 q& D& y
figure;
% ?7 f+ Q% f# l, F6 tsubplot(1, 2, 1); imshow(I); title('原图像');) |: D' K5 U* Y6 D+ d4 v0 d' d- P
subplot(1, 2, 2); imshow(I1); title('加噪声图像'); ]8 E& ]4 r: S8 `& O* w
[PSNR, MSE] = psnr(I, I1)# S$ h6 D& a/ @5 g4 b1 }# [6 R
% 图像峰值信噪比PSNR、均方根误差MSE
% E* q$ g2 }# H) F, k+ c8 @ ufunction [PSNR, MSE] = psnr(X, Y)
- v4 Q$ k; b2 ]- @7 \% 计算峰值信噪比PSNR、均方根误差MSE
* }8 _+ K6 z" Y* E% 如果输入Y为空,则视为X与其本身来计算PSNR、MSE
' L0 t [3 N9 [; M2 p% @5 J7 h6 B2 x) b% F' s' `
if nargin<2
4 s, d; F6 o. H D = X;
9 e6 h: o, l3 f, k/ q# P W8 [else" I, T% D" S f0 z! l. b) n
if any(size(X)~=size(Y))
, F" Z6 p8 i" F1 [5 Y error('The input size is not equal to each other!'); L% O+ Q4 v; N1 V# n4 c+ ]$ h, t
end; d; W- t5 C6 c7 U7 z, B
D = X-Y;3 H3 |/ I/ ]* d( w, e' |) [
end
O3 A( A; [ t. z3 H% v; LMSE = sum(D(: ).*D(: ))/prod(size(X));
( V* a3 }9 D# l% p) G/ Z: p4 [PSNR = 10*log10(255^2/MSE);
, g* L7 b3 t# @& S用这个算法算下吧,I为原图 I1为加噪图像/ s1 x2 D& g8 H; U- S% r
|
|