How to plot gaussian curve for each data point of boundary ?

2 visualizaciones (últimos 30 días)
Charu
Charu el 2 de Jun. de 2023
Respondida: Image Analyst el 8 de Ag. de 2023
  3 comentarios
Charu
Charu el 6 de Jun. de 2023
Its an extracted boundary from an image . Actually i want to calculate RMSF (ROOT MEAN SQUARE FLUCTUATION ) at each point between two images (deviated one and non deviated one)..and want to plot as pdf (probability density function) as gaussian curves with the calculation of skewness and kurtosis.
Mann Baidi
Mann Baidi el 8 de Ag. de 2023
It would be helpful if you can describe more about the issue

Iniciar sesión para comentar.

Respuestas (1)

Image Analyst
Image Analyst el 8 de Ag. de 2023
If you have the Statistics and Machine Learning Toolbox, what I'd try is pdist2 and put in the boundaries of both curves. Then for each row you can find out which two points are closest to each other. Then get the average of those closest distances. Something like (untested)
distanceMatrix = pdist2(xyRef, xyTest);
[rows, columns] = size(distanceMatrix)
% Allocate arrays for closest distances.
closestColDistances = nan(1, columns);
closestRowDistances = nan(1, rows);
% Get rid of zeros (distance of point to itself)
distanceMatrix(distanceMatrix==0) = inf; % Set to infinity.
% Get distances of ref curve to test curve.
for row = 1 : rows
closestRowDistances = min(distanceMatrix(row, :));
end
% Get distances of test curve to ref curve.
for col = 1 : columns
closestColDistances = min(distanceMatrix(:, col));
end
% Get the average distance.
averageDistance = mean([closestRowDistances, closestColDistances])
You need to get the distances both ways because if you only look at, say, the distance from the ref to the test, there may be some points on the test curve that are not closest and then they would not contribute to the average distance.

Categorías

Más información sobre Statistics and Machine Learning Toolbox en Help Center y File Exchange.

Community Treasure Hunt

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

Start Hunting!

Translated by