3

たとえば、長方形を表すベクトルを取得した場合

xy=[165  88; 
    401  88; 
    401 278; 
    165 278];

画像上。

次のベクトルを取得するにはどうすればよいですか

[165  88; % increase X -     hold y
 166  88; 
 167  88;
   ...  ;
 399  88;
 400  88;
 401  88; % hold     x - increase y
 401  89;
 401  90;
 401  91;
   ...  ;
 401 276;
 401 277;
 401 278; % decrease X -     hold y
 400 278;
 399 278;
 398 278;
   ...  ;
 167 278;
 166 278;
 165 278; % hold     x - decrease y
 165 277;
 165 276;
   ...  ;
 165  87];

MATLAB 組み込み関数を使用するか、FOR LOOPS を使用して記述する必要がありますか?

このアルゴリズムは、n 点と xy 座標を持つ一般的なベクトルに対して機能する必要があります。

4

3 に答える 3

4

画像処理ツールボックスがある場合は、多角形の画像を作成してから輪郭を見つけることでこれを行うことができます。

xy=[165 88; 401 88; 401 278; 165 278];
%# create the image - check the help for impolygon for how to make sure that
%# your line is inside the pixel
img = poly2mask(xy(:,1),xy(:,2),max(xy(:,1))+3,max(xy(:,2))+3);
figure,imshow(img) %# show the image

%# extract the perimeter. Note that you have to inverse x and y, and that I had to
%# add 1 to hit the rectangle - this shows one has to be careful with rectangular 
%# polygons
boundary = bwtraceboundary(logical(img),xy(1,[2,1])+1,'n',8,inf,'clockwise');

%# overlay extracted boundary
hold on, plot(boundary(:,2),boundary(:,1),'.r')

bwtraceboundary の使用方法を示し、長方形でのピクセル オフセットを警告するように編集されました。

于 2010-03-01T19:53:19.080 に答える
0

IND2SUB を使用した 1 つのソリューション:

xy=[165 88; 401 88; 401 278; 165 278];
xmin = min(xy(:,1))-1;
xmax = max(xy(:,1));
ymin = min(xy(:,2))-1;
ymax = max(xy(:,2));

ncol=xmax-xmin;
nrow=ymax-ymin;

[xn yn]=ind2sub([nrow ncol],1:nrow*ncol);
xypairs = [xn'+xmin yn'+ymin];
于 2010-03-01T19:16:47.080 に答える
0

画面外のマトリックスにストレートを描画する手っ取り早い方法は、数式を評価することa*X+b*Y=cです。

h と w をバッファーの幅と高さとします。

X = repmat([0:w-1], h, 1)
Y = repmat([0:h-1]', 1, w)

点 (x1,y1)->(x2,y2) のすべてのペアについて、a、b、および c は次のとおりです。

a = y2-y1
b = x1-x2
c = x1*y2-x2*y1

今ストレートを計算します：

st = a*X+b*Y-c
st(abs(st)>1) = 1
st = 1 - abs(st)

Matrixstは、点 (x1,y1) と (x2,y2) を通るアンチエイリアス処理された直線を含む w*h 行列です。不要な部分をマスクして、直線から直線に進みましょう。

[xs] = sort([x1 x2])
st = st .* [zeros(h, xs(1)) ones(h, xs(2)-xs(1)) zeros(h, w-xs(2))]
[ys] = sort([y1 y2])
st = st .* [zeros(ys(1), w) ; ones(ys(2)-ys(1), w) ; zeros(h-ys(2), w)]

明示的なループなしで、手動で 1 本の線を描画しました。ただし、コードの効率についての保証はありません:-)

最後に: 上記のすべての式に別のディメンションを追加します (読者の演習として残しています)。

于 2010-03-01T19:53:11.380 に答える