imread进来的影像资料
Clear all; close all;
DR = 60;
????? = imread(????); % imread進來的影像資料 data type 為 "uint8" ?????= rgb2gray(?????); % 轉換為gray scale
GrayIm = double(?????); % 將uint8的data type轉成"double
% gray to dB 由0-255的灰階轉成 dB
dBIm = GrayIm - min(min(GrayIm)); dBIm = dBIm/max(max(dBIm)); dBIm = dBIm*DR;
figure; % show B-mode image
image(?????);
colormap(gray(DR));
axis image;
colorbar;
title('?????');
xlabel('?????');
ylabel('?????');
% ------------------ speckle std ------------------
speckleStd = std2(dBIm(???:???,???:???)); % 將影像ROI均質的部份選出來計算 % -------------------- speckle histogram -----------------------
% 取影像上均質的部份出來計算histogram
% 計算histogram前,得先將dB資料,轉為原來的linear的資料格式 % dB = 20*log10(E), E: amplitude => E = 10^(dB/20)
% dB = 10*log10(I), I: intensity => I = 10^(dB/10)
LinearIm_I =?????? % intensity
LinearIm_E = ????? % amplitude
figure
hist(reshape(LinearIm_I,prod(size(LinearIm_I)),1),100);
% reshape將2-D data 轉成1-D,再使用"hist"來求出及畫出histgoram:機率分布: title('?????');
xlabel('?????');
ylabel('?????');
figure
hist(reshape(LinearIm_E,prod(size(LinearIm_E)),1),100);
title('?????');
xlabel('?????');
ylabel('?????');