Hello,

I’m trying to calculate the curvature of each surface point from a morphology picture. However, when using surfature(), I do not get the results I want. A lot of points have a gaussian curvature of 0, which is not plausible for the example I use.
I tried a lot of different methods for calculating and ploting, but never got a good result…

I wrote the following:
% reading in picture
data = imread("Morphology.PNG");

% if not yet, convert to gray image
if size(data, 3) == 3
data = rgb2gray(data);
end

% smoothing image
data = imgaussfilt(data, 4);

% converting to double
data = double(data);
z = data;

% generating 2D arrays for X and Y with size of data-dimensions
[Rows, Cols] = size(data);
[x, y] = meshgrid(1:Cols, 1:Rows);

% calculating gaussian and normal curvature
[K, H] = surfature(x, y, data);

% display surface
figure;
surf(x, y, z, "EdgeColor","none");
title(‘Surface Morphology’);
colorbar;

% ploting the results
figure;
surf(x, y, z, K, "EdgeColor","none");
title(‘3D Surface Gaussian Curvature Map’);
colorbar;

disp(K);Hello,

I’m trying to calculate the curvature of each surface point from a morphology picture. However, when using surfature(), I do not get the results I want. A lot of points have a gaussian curvature of 0, which is not plausible for the example I use.
I tried a lot of different methods for calculating and ploting, but never got a good result…

I wrote the following:
% reading in picture
data = imread("Morphology.PNG");

% if not yet, convert to gray image
if size(data, 3) == 3
data = rgb2gray(data);
end

% smoothing image
data = imgaussfilt(data, 4);

% converting to double
data = double(data);
z = data;

% generating 2D arrays for X and Y with size of data-dimensions
[Rows, Cols] = size(data);
[x, y] = meshgrid(1:Cols, 1:Rows);

% calculating gaussian and normal curvature
[K, H] = surfature(x, y, data);

% display surface
figure;
surf(x, y, z, "EdgeColor","none");
title(‘Surface Morphology’);
colorbar;

% ploting the results
figure;
surf(x, y, z, K, "EdgeColor","none");
title(‘3D Surface Gaussian Curvature Map’);
colorbar;

disp(K); Hello,

I’m trying to calculate the curvature of each surface point from a morphology picture. However, when using surfature(), I do not get the results I want. A lot of points have a gaussian curvature of 0, which is not plausible for the example I use.
I tried a lot of different methods for calculating and ploting, but never got a good result…

I wrote the following:
% reading in picture
data = imread("Morphology.PNG");

% if not yet, convert to gray image
if size(data, 3) == 3
data = rgb2gray(data);
end

% smoothing image
data = imgaussfilt(data, 4);

% converting to double
data = double(data);
z = data;

% generating 2D arrays for X and Y with size of data-dimensions
[Rows, Cols] = size(data);
[x, y] = meshgrid(1:Cols, 1:Rows);

% calculating gaussian and normal curvature
[K, H] = surfature(x, y, data);

% display surface
figure;
surf(x, y, z, "EdgeColor","none");
title(‘Surface Morphology’);
colorbar;

% ploting the results
figure;
surf(x, y, z, K, "EdgeColor","none");
title(‘3D Surface Gaussian Curvature Map’);
colorbar;

disp(K); surfature, surface, curvature, image processing MATLAB Answers — New Questions

​

I would analyze the data of MAXM86146I would analyze the data of MAXM86146 I would analyze the data of MAXM86146 sofware, maxm86146 MATLAB Answers — New Questions

​

Hi all,

simply, how to make the simple hollow 3D triangle depicted below?

Is there some contraint method I could use instead? I use the following command to get the surface:
tri = triangulation(delaunay(x,y),x,y,z);Hi all,

simply, how to make the simple hollow 3D triangle depicted below?

Is there some contraint method I could use instead? I use the following command to get the surface:
tri = triangulation(delaunay(x,y),x,y,z); Hi all,

simply, how to make the simple hollow 3D triangle depicted below?

Is there some contraint method I could use instead? I use the following command to get the surface:
tri = triangulation(delaunay(x,y),x,y,z); delaunay, triangulation MATLAB Answers — New Questions

​

Hi,
My computations are stored in a Clickhouse databse server on a Linux (Ubuntu) virtual machine (VM). I would like to access the data from the Clickhouse server on the VM to MATLAB in my local machine on windows OS, so that I can process the data locally.
I tried to explore the database tool box, but since I am new to this, I haven’t been able to connect to the database on the VM.
I would appreciate pointers from the experts.
Thanks in advance.Hi,
My computations are stored in a Clickhouse databse server on a Linux (Ubuntu) virtual machine (VM). I would like to access the data from the Clickhouse server on the VM to MATLAB in my local machine on windows OS, so that I can process the data locally.
I tried to explore the database tool box, but since I am new to this, I haven’t been able to connect to the database on the VM.
I would appreciate pointers from the experts.
Thanks in advance. Hi,
My computations are stored in a Clickhouse databse server on a Linux (Ubuntu) virtual machine (VM). I would like to access the data from the Clickhouse server on the VM to MATLAB in my local machine on windows OS, so that I can process the data locally.
I tried to explore the database tool box, but since I am new to this, I haven’t been able to connect to the database on the VM.
I would appreciate pointers from the experts.
Thanks in advance. database, clickhouse MATLAB Answers — New Questions

​

I am working with a finite element model and want to decompose my system using the eigenvectors of the system. The eigenvalue problem is of the form: K*V = M*V*D, where V is the eigenvector matrix and D is the eigenvalue matrix. Both M and K are block diagonal symmetric matrices. For the generalized eigenvalue problem it should be the case that V’*M*V = a matrix whose only nonzero terms are along the diagonal, but I am getting nonzero terms on either side of these. Using the code included below for something like N =5, I get (after rounding to nearest 3rd decimal to remove near zero terms)

abs( Orthgonal_Check_eig) = [ 78.8050 0 0 0 0
0 43.4540 0 0 0
0 0 0 34.9310 0
0 0 34.9310 0 0
0 0 0 0 29.6430]
Where I am taking the absolute value because this case has complex eigenvectors (though I believe I’ve seen it with real eigenvector matrices as well). If I rerun the code it does not always happen, but I am trying to understand why this can occur for both eig() and eigs() so that I can determine if it is the conditioning of the matrices or if it is inherent to the approximations inside of eig() and eigs(). Is there a better one of the two to use to avoid this issue?

N = 5;
M = zeros(N,N);
K = zeros(N,N);

for ii = 1:N
% Diagonal Terms
M(ii,ii) = rand*50;
K(ii,ii) = rand*500;

if ii ~=N
% Off Diagonal Terms
off_diagonal = rand*50;
M(ii,ii+1) = off_diagonal;
K(ii,ii+1) = off_diagonal*10;
M(ii+1,ii) = off_diagonal;
K(ii+1,ii) = off_diagonal*10;
end
end
M= round(M);
K= round(K);
[Veig,~] = eig(K,M);
[Veigs,~] = eigs(K,M,3);

Orthgonal_Check_eig = Veig’*M*Veig;
Orthgonal_Check_eigs =Veigs’*M*Veigs;I am working with a finite element model and want to decompose my system using the eigenvectors of the system. The eigenvalue problem is of the form: K*V = M*V*D, where V is the eigenvector matrix and D is the eigenvalue matrix. Both M and K are block diagonal symmetric matrices. For the generalized eigenvalue problem it should be the case that V’*M*V = a matrix whose only nonzero terms are along the diagonal, but I am getting nonzero terms on either side of these. Using the code included below for something like N =5, I get (after rounding to nearest 3rd decimal to remove near zero terms)

abs( Orthgonal_Check_eig) = [ 78.8050 0 0 0 0
0 43.4540 0 0 0
0 0 0 34.9310 0
0 0 34.9310 0 0
0 0 0 0 29.6430]
Where I am taking the absolute value because this case has complex eigenvectors (though I believe I’ve seen it with real eigenvector matrices as well). If I rerun the code it does not always happen, but I am trying to understand why this can occur for both eig() and eigs() so that I can determine if it is the conditioning of the matrices or if it is inherent to the approximations inside of eig() and eigs(). Is there a better one of the two to use to avoid this issue?

N = 5;
M = zeros(N,N);
K = zeros(N,N);

for ii = 1:N
% Diagonal Terms
M(ii,ii) = rand*50;
K(ii,ii) = rand*500;

if ii ~=N
% Off Diagonal Terms
off_diagonal = rand*50;
M(ii,ii+1) = off_diagonal;
K(ii,ii+1) = off_diagonal*10;
M(ii+1,ii) = off_diagonal;
K(ii+1,ii) = off_diagonal*10;
end
end
M= round(M);
K= round(K);
[Veig,~] = eig(K,M);
[Veigs,~] = eigs(K,M,3);

Orthgonal_Check_eig = Veig’*M*Veig;
Orthgonal_Check_eigs =Veigs’*M*Veigs; I am working with a finite element model and want to decompose my system using the eigenvectors of the system. The eigenvalue problem is of the form: K*V = M*V*D, where V is the eigenvector matrix and D is the eigenvalue matrix. Both M and K are block diagonal symmetric matrices. For the generalized eigenvalue problem it should be the case that V’*M*V = a matrix whose only nonzero terms are along the diagonal, but I am getting nonzero terms on either side of these. Using the code included below for something like N =5, I get (after rounding to nearest 3rd decimal to remove near zero terms)

abs( Orthgonal_Check_eig) = [ 78.8050 0 0 0 0
0 43.4540 0 0 0
0 0 0 34.9310 0
0 0 34.9310 0 0
0 0 0 0 29.6430]
Where I am taking the absolute value because this case has complex eigenvectors (though I believe I’ve seen it with real eigenvector matrices as well). If I rerun the code it does not always happen, but I am trying to understand why this can occur for both eig() and eigs() so that I can determine if it is the conditioning of the matrices or if it is inherent to the approximations inside of eig() and eigs(). Is there a better one of the two to use to avoid this issue?

N = 5;
M = zeros(N,N);
K = zeros(N,N);

for ii = 1:N
% Diagonal Terms
M(ii,ii) = rand*50;
K(ii,ii) = rand*500;

if ii ~=N
% Off Diagonal Terms
off_diagonal = rand*50;
M(ii,ii+1) = off_diagonal;
K(ii,ii+1) = off_diagonal*10;
M(ii+1,ii) = off_diagonal;
K(ii+1,ii) = off_diagonal*10;
end
end
M= round(M);
K= round(K);
[Veig,~] = eig(K,M);
[Veigs,~] = eigs(K,M,3);

Orthgonal_Check_eig = Veig’*M*Veig;
Orthgonal_Check_eigs =Veigs’*M*Veigs; eigenvalue problem, eig, eigs, orthogonality MATLAB Answers — New Questions

​

Hello, I can´t answer the last task, task 9, of the section 8.1 of the simscape onramp course, I do exactly what I´m requested to do but I always get it wrong, I´ve been for several hours trying it but I simply don´t understand what I´m getting wrong. I´ve looked for the answer on Youtube but literally nobady has upload it, at least I haven´t found anybody, does anybody knows how to obtain the answer or can anybody tell me where can I find the way to solve it?, thank you!Hello, I can´t answer the last task, task 9, of the section 8.1 of the simscape onramp course, I do exactly what I´m requested to do but I always get it wrong, I´ve been for several hours trying it but I simply don´t understand what I´m getting wrong. I´ve looked for the answer on Youtube but literally nobady has upload it, at least I haven´t found anybody, does anybody knows how to obtain the answer or can anybody tell me where can I find the way to solve it?, thank you! Hello, I can´t answer the last task, task 9, of the section 8.1 of the simscape onramp course, I do exactly what I´m requested to do but I always get it wrong, I´ve been for several hours trying it but I simply don´t understand what I´m getting wrong. I´ve looked for the answer on Youtube but literally nobady has upload it, at least I haven´t found anybody, does anybody knows how to obtain the answer or can anybody tell me where can I find the way to solve it?, thank you! simscape MATLAB Answers — New Questions

​

Hello, I am developing some control logic in Virtual Vehicle Composer (VVC). There are some example scenarios like the one below. However, when I use the "Double lane change" scenario, the steering response is a bit weird (too sudden compared to other simulations such as CarSim). So, I want to know how to change the scenario or manually replace the driver model.
I tried once, but it was difficult to consider the related signals from other vehicles. If anyone has attempted this or if I could get some help regarding this, I would appreciate it.Hello, I am developing some control logic in Virtual Vehicle Composer (VVC). There are some example scenarios like the one below. However, when I use the "Double lane change" scenario, the steering response is a bit weird (too sudden compared to other simulations such as CarSim). So, I want to know how to change the scenario or manually replace the driver model.
I tried once, but it was difficult to consider the related signals from other vehicles. If anyone has attempted this or if I could get some help regarding this, I would appreciate it. Hello, I am developing some control logic in Virtual Vehicle Composer (VVC). There are some example scenarios like the one below. However, when I use the "Double lane change" scenario, the steering response is a bit weird (too sudden compared to other simulations such as CarSim). So, I want to know how to change the scenario or manually replace the driver model.
I tried once, but it was difficult to consider the related signals from other vehicles. If anyone has attempted this or if I could get some help regarding this, I would appreciate it. simulink, virtualvehiclecomposer MATLAB Answers — New Questions

​

I have been assigned a project which works by linking Aspen Plus with Matlab. What my work is to extract multiple data points to matlab using actx server. For example, I was able to get temperature & pressure data as they are single variable whereas, for Particle Size Distribution(PSD) data it cosists of an array of elements. We have mess sizes & mole fraction of crystals. How do I tell matlab to get that data from Aspen Plus? Please help me with this doubt. Currently using Aspen V12 version.I have been assigned a project which works by linking Aspen Plus with Matlab. What my work is to extract multiple data points to matlab using actx server. For example, I was able to get temperature & pressure data as they are single variable whereas, for Particle Size Distribution(PSD) data it cosists of an array of elements. We have mess sizes & mole fraction of crystals. How do I tell matlab to get that data from Aspen Plus? Please help me with this doubt. Currently using Aspen V12 version. I have been assigned a project which works by linking Aspen Plus with Matlab. What my work is to extract multiple data points to matlab using actx server. For example, I was able to get temperature & pressure data as they are single variable whereas, for Particle Size Distribution(PSD) data it cosists of an array of elements. We have mess sizes & mole fraction of crystals. How do I tell matlab to get that data from Aspen Plus? Please help me with this doubt. Currently using Aspen V12 version. matlab, aspenplus, actxserver MATLAB Answers — New Questions

​

I have a variable that varies from 0 to 1 and its format is 32 bit fixdt(1,32,30). I want to convert this number 0 to 1 as 0-4095. Help is required for converting fixdt(1,32,30) to 12 bits integer value ranging from 0 to 4095.I have a variable that varies from 0 to 1 and its format is 32 bit fixdt(1,32,30). I want to convert this number 0 to 1 as 0-4095. Help is required for converting fixdt(1,32,30) to 12 bits integer value ranging from 0 to 4095. I have a variable that varies from 0 to 1 and its format is 32 bit fixdt(1,32,30). I want to convert this number 0 to 1 as 0-4095. Help is required for converting fixdt(1,32,30) to 12 bits integer value ranging from 0 to 4095. data type conversion MATLAB Answers — New Questions

​

strange behavior:
if I work in my matlab source with fft and ifft and then call in the same source my mex-function which works with fftw_threads and of course fftw,ifft, it works very good, but after closing matlab, I’ll get "Segmentation violation detected" with something like fftw_plan_awake….
If I do this process the other way around, first call of my mex function and afterwards fft,ifft in my matlab source, closing matlab, I have no "Segmentation violation detected".
As well, if I run my mex-function without fftw_threads everything is fine!! (irrelevant if I call matlab fft then mex fct or mex fct then matlab fft)
In both source, matlab and c-mex-source, I only make some easy standard fft, nothing special. New is the try to work with fftw_threads.

My Question: what happens in the background so that matlab crashes after closing in some constellations?
Thank you!strange behavior:
if I work in my matlab source with fft and ifft and then call in the same source my mex-function which works with fftw_threads and of course fftw,ifft, it works very good, but after closing matlab, I’ll get "Segmentation violation detected" with something like fftw_plan_awake….
If I do this process the other way around, first call of my mex function and afterwards fft,ifft in my matlab source, closing matlab, I have no "Segmentation violation detected".
As well, if I run my mex-function without fftw_threads everything is fine!! (irrelevant if I call matlab fft then mex fct or mex fct then matlab fft)
In both source, matlab and c-mex-source, I only make some easy standard fft, nothing special. New is the try to work with fftw_threads.

My Question: what happens in the background so that matlab crashes after closing in some constellations?
Thank you! strange behavior:
if I work in my matlab source with fft and ifft and then call in the same source my mex-function which works with fftw_threads and of course fftw,ifft, it works very good, but after closing matlab, I’ll get "Segmentation violation detected" with something like fftw_plan_awake….
If I do this process the other way around, first call of my mex function and afterwards fft,ifft in my matlab source, closing matlab, I have no "Segmentation violation detected".
As well, if I run my mex-function without fftw_threads everything is fine!! (irrelevant if I call matlab fft then mex fct or mex fct then matlab fft)
In both source, matlab and c-mex-source, I only make some easy standard fft, nothing special. New is the try to work with fftw_threads.

My Question: what happens in the background so that matlab crashes after closing in some constellations?
Thank you! mex fftw_threads, crash after closing matlab MATLAB Answers — New Questions

​