Determine the condition to find elements inside semi circular region in a rectangular mesh

6 ビュー (過去 30 日間)
Abhishek Saini
Abhishek Saini 2018 年 11 月 6 日
コメント済み: Abhishek Saini 2018 年 12 月 4 日
I have a rectangular mesh of 500000x500000. I extracted the coordinates of the mesh elements. Now I would like to remove elements in the semicircular form from this rectangular mesh. I have the function to remove elements. But I am struggling to set the condition for the semicircular region to determine all the elements inside the semicircular region. For example if ex and ey are the coordinates of elements in mesh, then I am using
find(ex<=R*cos(theta) & ey <=R*sin(theta)).
How to set this condition according to the grid points in the rectangular mesh. I have also attached the image for more clarity.
Thanks.

採用された回答

KSSV
KSSV 2018 年 11 月 6 日
編集済み: KSSV 2018 年 11 月 6 日
Read about inpolygon. This should help you.
M = 100 ;
N = 100 ;
x = linspace(0,1,M) ;
y = linspace(0,1,N) ;
[X,Y] = meshgrid(x,y) ;
% semi circle
O = [0.5 0. ]; % origin of circle
R = 0.3 ; % RAdius of circle
th = linspace(0,2*pi) ;
xc = O(1)+R*cos(th) ;
yc = O(2)+R*sin(th) ;
% Get points lying inside circle
idx = inpolygon(X(:),Y(:),xc,yc) ;
% plot
figure
hold on
plot(X,Y,'.r') ;
plot(xc,yc,'b')
plot(X(idx),Y(idx),'.k')
  11 件のコメント
Abhishek Saini
Abhishek Saini 2018 年 12 月 4 日
Whether this question will be considered or not?

サインインしてコメントする。

その他の回答 (0 件)

タグ

Community Treasure Hunt

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

Start Hunting!

Translated by