ダムMatlab問題

Question

私はこれが何か簡単だと確信していますが、私はむしろそれに困惑しています。私は、それぞれ異なる露出時間で撮影された一連の画像を取り込み、出力の平均値（DN）と信号のRMSを測定しています。全体のポイントは、光子伝達曲線を構築することです。私が実行している問題は、ウィンドウの信号またはRMSの平均を取ろうとするとき、または1つのピクセルさえも平均して保存しようとすると、格納された値がわかりにくいということです。これはむしろまっすぐ進むはずなので、それはとてもイライラしています。どんな洞察にも感謝します。

pix 
    16 
    18 
    21 
    19 
    22 
    17 
    19 
    18 
    16 
    20 
    18 
    21 
    23 
    17 
    20 

pixmean 
    19 

pixrms 
frame(t) 
    19 

framerms(t) 
    18.8880 

mean(frame) 
    0.4404 

frameMean(i) 
    0.4404 

********** 
frameRMS(i) 
    0.4404 

********** 

は、あなたがこのいずれかを明確にしたい場合は私に知らせてください：

stacklength = 25; 
frame = []; 
frameMean = []; 
framerms = []; 
frameRMS = []; 
t = 1; 
stacknum = 99; 
stackdepth = 15; 
imagewidth=1944; 
imageheight=2592; 
uStart = 1001; 
uEnd = 1001; 
vStart = 1001; 
vEnd = 1001; 
tic; 
for i=1:stacklength 
    filename=sprintf('PTLstack_%d.mat' ,i); %name files one at a time 
    load(filename); 
    for u = uStart:uEnd 
     for v = vStart:vEnd 
      pix = double(squeeze(pixel(u,v,:))); 
      pixmean = mean(pix); 
      pixrms = rms(pix); 
      frame(t) = pixmean; 
      framerms(t) = abs(pixrms-pixmean); 
      t = t+1; 
     end 
    end 
    frameMean(i) = mean(frame); 
    frameRMS(i) = mean(framerms); 
    t = 1; 
    disp(i); 
end 
toc; 

そして、ここでは、私はそれが私のために表示されていたいくつかの値です：

は、ここに私のコードです。私はコードの塊を見て、その中で起こっていること全てを理解することが期待されるのは面倒かもしれません。

ありがとう、10^6。

Answer 1

結果には、frame変数にいくつかの迷子ゼロがあり、平均値が低くなることが示唆されています。これらがどこに来るかは私には分かりませんが、これはこの種のコードを扱うMatlabスタイルの方法ではないからです。あなたはまた、いくつかの偽のデータを有していてもよく、予備のゼロを持っている場合は

frameMean(i) = mean(frame(frame~=0)); 

：

ショートカット、洗練ものの、これを修正する唯一の非ゼロの平均を取るために、それを伝えることであろうそこに他の並べ替え。

もっと良い解決策は、これらのループと増分変数などを完全に取り除くことです。あなたのイメージのサイズの画像スタック（X、Y、15）としてロードされていると仮定すると、あなたが意味かかることがありますが、単一行くに望む地域全体のrms：

pix = pix(ustart:uend,vstart:vend,:); % take window area of choice 
frame = mean(pix,3); % make mean 
framerms = rms(pix,3); 
framerms = abs(framerms-frame); 

ここでの出力は、2D行列でありますフレーム全体にわたって。フレーム全体に渡って平均を取るのが最も簡単です。 mean(frame(:))。