画像内の余分なところを除いて新しい行列にしたい
6 Ansichten (letzte 30 Tage)
Ältere Kommentare anzeigen
thuru
am 19 Jan. 2017
Kommentiert: thuru
am 26 Jan. 2017
添付してある画像のように、解析に用いたい部分が中央にあるのですが、そのまま使うと余計な数値が入ってしまい困っています。 添付画像は256×256 uint16です。周りのグレーの行列の値は32768です。32768の部分が邪魔になっています。トリミングみたいに、余計な部分を省いて新しい行列を作成したいです。
3 Kommentare
Akzeptierte Antwort
Takuji Fukumoto
am 24 Jan. 2017
画像の周囲と中央の境界になっている部分を検出してトリミングするサンプルです。
I = imread('grayscale2.png');
Ig = rgb2gray(I);
Ie = edge(Ig,'Sobel',0.00001);
Ie(end-3:end,:)= 0;
[row,col] = find(Ie);
rmax = max(row);rmin = min(row);cmax = max(col);cmin = min(col);
I2 = I(rmin:rmax,cmin:cmax,:);
figure,imshow(I2)
画像のトリミングでは直交する行列方向での切り出しになります。 領域周囲の直線が斜めになっているところは残ってしまいます。
領域の解析だけもう少し厳密にするのであれば、
h = impoly(gca, [188,30; 189,142; 93,141; 13,41]);
BW = createMask(h);
I(~BW) =0;
によってマスク以外のところはゼロにすることができるので、 Ieの頂点座標を使うか、マウス操作でimpolyの頂点を動かして領域を決定して用いることが出来ます。
3 Kommentare
Takuji Fukumoto
am 25 Jan. 2017
Bearbeitet: Takuji Fukumoto
am 25 Jan. 2017
各処理後の領域の最大最小をあらかじめ求めてから、全部同じ座標で切り出すのではいかがでしょうか。
for n=1:67
In = imread('Filename');
Ig = rgb2gray(In);
Ie = edge(Ig,'Sobel',0.00001);
Ie(end-3:end,:)= 0;
[row,col] = find(Ie);
rmax(n) = max(row);rmin(n) = min(row);cmax(n) = max(col);cmin(n) = min(col);
end
rmaxall = max(rmax);
rminall = min(rmin);
cmaxall = max(cmax);
cminall = min(cmin);
for k=1:67
In = ('Filename');
Inout = In(rminall:rmaxall,cminall:cmaxall,:);
end
画像の読み出しと出力のところは編集してご利用ください。
Weitere Antworten (2)
michio
am 23 Jan. 2017
上記の「添字による代入の次元が一致しません。」というエラーですが、
grayscale(:,:,k) = imread(tifFiles(k).name);
の実行において grayscale は 67×67×numfiles の配列ですが、imread で読み取られる画像サイズが 256×256 となっており、配列のサイズが合致していないことが原因です。
for k = 2:numfiles
img = imread(tifFiles(k).name);
grayscale(:,:,k) = img(58:124, 110:176);
end
と配列のサイズを合わせて grayscale に代入してやってください。
また、もし余分なところを取り除いた後に残る中心部分の画像のサイズが事前に分かっていれば、下記の方法が使えるかなと。
5x5の2次元配列から、100である部分を除いて3x2の配列を取り出しているサンプルになります。
% サンプルデータ(周囲が値100で囲まれた 3x2 の画像)
imag1 = 100*ones(5,5);
imag1(2:4,2:3) = rand(3,2)
copy = imag1;
% 100の値を持つ要素を削除
% 中心部分のデータが抽出できますが、すべての要素が
% 一列に並んだベクトルとして(ここでは 1x6 のベクトル)でてきます。
copy(copy==100) = [];
% 一列に並んだデータを2次元に並べ替え (ここで、事前に中心部分の画像サイズが必要)
imag2 = reshape(copy,3,2)
0 Kommentare
michio
am 23 Jan. 2017
Bearbeitet: michio
am 23 Jan. 2017
下記を実行して、変数 row や clm などの値を確認してみてください。
また、
copy(:,row) = []; % 要素がすべて1の行を削除
など、= [] で対応する要素を削除することができることを利用します。
% サンプルデータ(周囲が値100で囲まれた 3x2 の画像)
imag1 = 100*ones(5,5);
imag1(2:4,2:3) = rand(3,2)
copy = imag1;
% 100の値を持つ要素位置に1をもった論理配列を作成
idx = copy == 100
row = all(idx,1); % 要素が全て1の行を検出
clm = all(idx,2); % 要素が全て1の列を検出
copy(:,row) = []; % 要素がすべて1の行を削除
copy(clm, :) = [];% 要素がすべて1の列を削除
copy
Siehe auch
Kategorien
Mehr zu Matrix Indexing finden Sie in Help Center und File Exchange
Produkte
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!