画像内の余分なところを除いて新しい行列にしたい
显示 更早的评论
采纳的回答
Takuji Fukumoto
2017-1-24
画像の周囲と中央の境界になっている部分を検出してトリミングするサンプルです。
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 个评论
thuru
2017-1-24
Takuji Fukumoto
2017-1-25
编辑：Takuji Fukumoto
2017-1-25
各処理後の領域の最大最小をあらかじめ求めてから、全部同じ座標で切り出すのではいかがでしょうか。
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
画像の読み出しと出力のところは編集してご利用ください。
thuru
2017-1-26
更多回答（2 个）
michio
2017-1-23
上記の「添字による代入の次元が一致しません。」というエラーですが、
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) = [];
% 一列に並んだデータを２次元に並べ替え (ここで、事前に中心部分の画像サイズが必要）
imag2 = reshape(copy,3,2)
下記を実行して、変数 row や clm などの値を確認してみてください。
また、
copy(:,row) = []; % 要素がすべて１の行を削除
など、= [] で対応する要素を削除することができることを利用します。
% サンプルデータ（周囲が値100で囲まれた 3x2 の画像）
imag1 = 100*ones(5,5);
imag1(2:4,2:3) = rand(3,2)
copy = imag1;
% 100の値を持つ要素位置に１をもった論理配列を作成
idx = copy == 100
row = all(idx,1); % 要素が全て１の行を検出
clm = all(idx,2); % 要素が全て１の列を検出
copy(:,row) = []; % 要素がすべて１の行を削除
copy(clm, :) = [];% 要素がすべて１の列を削除
copy
类别
在 帮助中心 和 File Exchange 中查找有关 Graphics Performance 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
