imread进来的影像资料★[doc]★

imread进来的影像资料★[doc]★ 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('?????');