# how to Identify the centre of atom

3 ビュー (過去 30 日間)
2024 年 4 月 12 日
コメント済み: DGM 2024 年 4 月 13 日
how to Identify the centre of gravity like the fllowing figure

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

### 採用された回答

DGM 2024 年 4 月 12 日
Here's a start.
% the original image
inpict = im2gray(inpict);
% crop out the subimages
A0 = imcrop(inpict,[31.51 30.51 277.98 277.98]);
B0 = imcrop(inpict,[326.51 30.51 277.98 277.98]);
C0 = imcrop(inpict,[621.51 30.51 277.98 277.98]);
% process A
A = imflatfield(A0,10);
A = imopen(A,strel('disk',3,0));
maskA = imextendedmax(A,50); % roughly get a mask of minimal regions
maskA = imclearborder(maskA); % get rid of partial blobs
% process B
maskB = imextendedmin(B0,80); % roughly get a mask of minimal regions
maskB = imclearborder(maskB); % get rid of partial blobs
% process C
C = imflatfield(C0,10);
maskC = imextendedmin(C,80); % roughly get a mask of minimal regions
maskC = imclearborder(maskC); % get rid of partial blobs
% for example, find the centroids of the blobs in one of the masks
centers = vertcat(Sc.Centroid);
% mark the centers on the mask
hold on;
plot(centers(:,1),centers(:,2),'x','linewidth',2)
##### 7 件のコメント5 件の古いコメントを表示5 件の古いコメントを非表示
Image Analyst 2024 年 4 月 13 日
@DGM in your first response, since he specifically wanted the center of gravity, you should have used "WeightedCentroid" instead of 'Centroid', though for this image it looks like they should be really close to the same location.
Maybe, after field flattening, you can use kmeans to classify a pixel as light, dark, or gray.
DGM 2024 年 4 月 13 日
That's a pretty good point. I'll throw in an example just for sake of completeness.
% the original image
inpict = im2gray(inpict);
% crop out the subimages (i'm just going to use C)
C0 = imcrop(inpict,[621.51 30.51 277.98 277.98]);
% process C (same as before)
C = imflatfield(C0,10);
maskC = imextendedmin(C,80); % roughly get a mask of minimal regions
maskC = imclearborder(maskC); % get rid of partial blobs
% for example, find the weighted centroids of the blobs in one of the masks
% in this call to regionprops, it's important to make sure that the second argument
% has a light-on-dark polarity, hence the use of imcomplement() in the case of B or C.
% A is already light-on-dark, so we wouldn't need to complement it.
centers = vertcat(Sc.WeightedCentroid);
plot(centers(:,1),centers(:,2),'x','linewidth',2)

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

### カテゴリ

Help Center および File ExchangeRegion and Image Properties についてさらに検索

### Community Treasure Hunt

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

Start Hunting!

Translated by