关于#matlab#绘制圆锥的问题，如何解决？

Question

我想在DEM中找到一个圆，然后生成一个顶点位于中心点的理想圆锥，就是把一个穹隆或者凹坑理想化成一个正圆锥，我试了一下，只能单独地用matlab读取dem和绘制一个圆锥，不能把两者放在一起，让圆锥有dem的坐标系统。这要怎么写代码？

Answer 1

参考GPT和自己的思路：

你可以通过以下步骤来绘制一个顶点在中心的圆锥，并将其与DEM的坐标系统对齐：

1. 通过Matlab读取DEM数据，并确定圆锥的中心点。
2. 使用Matlab的plot3函数绘制一个圆锥，其中圆锥的底面圆是从中心点开始的固定半径的圆。
3. 可以使用surf函数将DEM数据映射为3D表面，并将其添加到圆锥上。需要将DEM数据从其原始坐标系转换为与圆锥相同的坐标系。
4. 将图形调整为最佳视角，以使DEM表面和圆锥更好地对齐。

下面是可能用到的一些代码示例：

读取DEM数据：

dem = imread('dem_data.tif');

确定圆锥的中心点：

center = [x,y,z]; % 坐标系统中圆锥的中心点

绘制圆锥：

r = 10; % 圆锥的半径
h = 20; % 圆锥的高度
theta = linspace(0, 2*pi, 50);
x = r*cos(theta);
y = r*sin(theta);
z = zeros(size(x));
zs = linspace(0, h, length(z));
zs(end) = h;
[X, Y, Z] = cylinder2P(x, y, z, zs);
figure;
plot3(X,Y,Z,'-','LineWidth',2,'Color','k');
axis equal;

将DEM数据映射到圆锥表面：

[X,Y]=meshgrid(1:size(dem,2),1:size(dem,1));
Z = double(dem);
% 坐标系偏移量，将DEM数据对齐到圆锥所在的坐标系
X = X - center(1) - size(dem,2)/2; 
Y = Y - center(2) - size(dem,1)/2;
Z = Z - center(3);
% 将DEM数据的坐标系旋转，以匹配圆锥坐标系
Rx = [1,0,0;0,cos(rad),-sin(rad);0,sin(rad),cos(rad)];
Ry = [cos(rad),0,-sin(rad);0,1,0;sin(rad),0,cos(rad)];
Rz = [cos(rad),-sin(rad),0;sin(rad),cos(rad),0;0,0,1];
R = Rx*Ry*Rz; % 旋转角度为rad
XYZ = R*[X(:)';Y(:)';Z(:)'];
X = reshape(XYZ(1,:),size(dem));
Y = reshape(XYZ(2,:),size(dem));
Z = reshape(XYZ(3,:),size(dem));
% 绘制DEM表面
hold on;
surf(X,Y,Z,'EdgeColor','none','FaceColor','texturemap','CData',dem);
axis equal;

最后，您可以使用视角调整函数来调整图形以使DEM表面和圆锥更好地对齐：

view(45,30);

希望这能帮到你！