How to accurately measure the height of a 3D structure?

6 次查看（过去 30 天）
Goncalo Costa
Goncalo Costa 2023-3-10
编辑： Mathieu NOE 2023-3-17
I have a 3D data array consisting of a scan of a semi-circular sample on a flat surface. The purpose of this is to image and afterwards detect the height and width of each sample.
To do so, I wrote the following:
data_raw = load('sample.asc'); %3 columns with hundreds of rows
data_filt = imgaussfilt3(data_raw);
x = data_filt(:,1);
y = data_filt(:,2);
z = data_filt(:,3);
dotsize = 6;
scatter3(x, y, z, dotsize, z, 'filled'); %colobar linked to height (z-axis)
xlim([-6 14])
ylim([-20 -6]) %limits where the sample is
I wrote the limits above, as the original image goes from [(-60, 60), (-40, 40) ,(-10, 15)], and this included a lot of noise and unnecessary image space (especially when using the function plot3 instead of scatter3).
To try and measure the height of the sample and its width I thought of doing a simple max/min subtraction within a certain space array, but unfortunately there is enough noise that would make this not work. Is there another way of doing this?
The images below demonstrate the noise and the image before the array limits.
The image above shows all the noise around of the surface on which the sample laid on and why I applied the above limits.
The image below shows the noise still around the area of interest. How do I get the height and width (accurately) of the sample in this case?
Update:
I have attempted using imdistline and drawline the first one can only be used in a single dimension and I don't know how to obtain the lengthof the latter.
  5 个评论
显示 3更早的评论
Mathieu NOE
Mathieu NOE 2023-3-15
you could share it via internet (like Googledrive or alike) and share the link
Goncalo Costa
Goncalo Costa 2023-3-15
编辑：Goncalo Costa 2023-3-15
I have uploaded it up in the problem desciption. Thanks for the suggestion.

请先登录，再进行评论。

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

采纳的回答

Mathieu NOE
Mathieu NOE 2023-3-15
编辑：Mathieu NOE 2023-3-15
hello
this is my rough first trial
NB that I didn't use here any surface smoothing (needs still to be implemented) but I assume you can bring that bak in the code
(I cannot use imgaussfilt3 because I don't have the corresponding toolbox)
one last detail, I draw the lower boundary line once I could remove the "noise" islands from the main shape . I don't know if you prefer to take the mean or the min of the z points belonging to that line , please choose the best option for you
hope this helps !
data_raw = readmatrix('drying_1.asc','FileType','text'); %3 columns with hundreds of rows
x = data_raw(:,1);
y = data_raw(:,2);
z = data_raw(:,3);
% "zoom" into sample area
idx = (x>-6 & x<14) & (y>-20 & y<-6);
x = x(idx);
y = y(idx);
z = z(idx);
dotsize = 6;
figure
scatter3(x, y, z, dotsize, z, 'filled'); %colobar linked to height (z-axis)
% attempt to remove outliers once x,y coordinates are transformed in polar
% coordinates
x_center = mean(x);
y_center = mean(y);
[th,r,z] = cart2pol(x-x_center,y-y_center,z);
% remove "outliers" points => find the best radius for circle fit
M = 100;
[N,EDGES] = histcounts(r,M); % uses M bins.
idx = find(N<2); % find first bin of data with very small amount of samples => "gap" between valid data and noise islands
idx = idx(1);
r_threshold = EDGES(idx); % this should be the radius of the circle that contains the valid data (outside this radius is noise)
% remove points above r_threshold
idx = (r<= r_threshold);
% convert back to cartesian coordinates
[x,y,z] = pol2cart(th(idx),r(idx),z(idx));
figure
scatter3(x, y, z, dotsize, z, 'filled'); %colobar linked to height (z-axis)
% boundary ?
k = boundary(x,y,0.5);
hold on
plot3(x(k), y(k), z(k), 'r'); %colobar linked to height (z-axis)
% height distance
% z_low = mean(z(k));
z_low = min(z(k));
z_high = max(z);
distZ = z_high - z_low
  7 个评论
显示 5更早的评论
Goncalo Costa
Goncalo Costa 2023-3-16
The reference plane, is just a plane on which the sample being analyzed was place, so yeah it would be that line that you drew. Perfect!
Mathieu NOE
Mathieu NOE 2023-3-17
编辑：Mathieu NOE 2023-3-17
there are several "steps" so I am not quite sure which one is the true reference plane
you see my red line does not truly coincide with both sides too , I think the whole "basement" is not truly parallel to the XY plane (more noticeable seen from the side)

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Data Exploration 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by