２次元グラフの面積取得方法

27 次查看（过去 30 天）
H.O
H.O 2023-3-30
评论： H.O 2023-4-9
２次元の座標データをもっていて
点に囲まれた座標の面積を求めたいです．
polyareaを使ってみましたが，形状が変わってしまいます．　（sample.png)
左図：もとめたい形状
右図：plot した画像　→　polyareaはこの形状の面積を求めている？
左図の形状を保った状態で，面積を取得する方法はありますか？
座標の分解能は保持したいです．
  6 个评论
显示 4更早的评论
Atsushi Ueno
Atsushi Ueno 2023-3-31
移动：Atsushi Ueno 2023-4-2
Ran in:
領域境界のタイプ - MATLAB & Simulink - MathWorks 日本
convhull関数はboundary関数の係数sが0の時であると知りました。同係数sを0⇒1に変化させていくと、2次元点群をラバーバンドで括った形状が徐々に締め付けられていくのが判ります。
n = 100;
for k = 0:9 % ランダム値でサンプル点群データを生成
x(k*n+1:(k+1)*n) = rand(100,1)*n + randi(10,1,1)*n;
y(k*n+1:(k+1)*n) = rand(100,1)*n + randi(10,1,1)*n;
end
plot(x,y,'*')
hold on
h = plot(x(k),y(k));
for s = 0:0.01:1
[k,av] = boundary(x',y',s); % 境界を計算%[k,av] = convhull(x,y); % 凸包を計算
%av % 面積
h.XData = x(k);
h.YData = y(k);
drawnow
end
H.O
H.O 2023-4-2
移动：Atsushi Ueno 2023-4-2
ありがとうございます．
実はアルファシェイプは最初に試したのですが，
各領域の面積や領域間との距離が小さいためか，アルファシェイプだと形状の変形が大きかったので，
ほかの方法を探していました．
申し訳ありません．

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Hiroshi Iwamura
Hiroshi Iwamura 2023-4-6
Ran in:
外枠いらなければ、最後のループも不要ですね。
I = imread("https://jp.mathworks.com/matlabcentral/answers/uploaded_files/1345389/sample_c.png");
BW = imbinarize(I(:,:,2));
se = strel('disk',2);
BW2 = imclose(BW,se);
BW3 = imopen(BW2,se);
montage({I,BW,BW2,BW3})
stats = regionprops(BW3,'basic');
L = bwlabel(BW3);
T = struct2table(stats);
imshow(L,[],Colormap=jet)
idx = (T.Area > 100);
text(T.BoundingBox(idx,1),T.Centroid(idx,2),num2str(T.Area(idx)),'Color','white','FontSize',10);
fprintf('Total Area = %d\n',sum(T.Area(T.Area>10)))
Total Area = 4289
  2 个评论
Atsushi Ueno
Atsushi Ueno 2023-4-8
的確な回答ですが、問題は要の regionprops 関数等に ImageProcessingToolbox が必要な点です。
と思ったのですが、以前の@H.Oさんの質問でImageProcessingToolboxを使う回答にOKのメッセージがあるので、良いみたいですね。
Hiroshi Iwamura
Hiroshi Iwamura 2023-4-8
なるほど
Computer Vision Toolbox もお持ちのようですね。

请先登录，再进行评论。

更多回答（1 个）

Hiroshi Iwamura
Hiroshi Iwamura 2023-4-4
编辑：Hiroshi Iwamura 2023-4-4
Ran in:
なかなか正確に求めるのは難しく調整が必要になりますが、Image Processing Toolbox をお持ちであればモフォロジーを使う手はあります。Simulink (Computer Vision Toolbox) でやった方が色々と調整が簡単かもしれません。
I = imread("sample_c.png");
BW = imbinarize(I(:,:,2));
se = strel('disk',2);
BW2 = imclose(BW,se);
BW3 = imopen(BW2,se);
montage({I,BW,BW2,BW3})
stats = regionprops(BW3,'basic');
L = bwlabel(BW3);
T = struct2table(stats);
T = sortrows(T,'Area','descend');
% imshow(I)
imshow(L,[],Colormap=jet)
hold on
n = 1;
while T.Area(n) > 100
rectangle('Position',T.BoundingBox(n,:),EdgeColor=[1 0.2 0])
text(T.BoundingBox(n,1),T.Centroid(n,2),num2str(T.Area(n)),'Color','white','FontSize',10)
n = n + 1;
end
hold off
fprintf('Total Area = %d\n',sum(T.Area(T.Area>10)))
Total Area = 4289
  1 个评论
H.O
H.O 2023-4-9
ありがとうございます。 必要なtoolboxも運良く所有していました。 挙げて頂いたスクリプトで実行、面積取得(pixel)できました。

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 境界領域 的更多信息

标签

产品


版本

R2022b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!