我想在DEM中找到一个圆，然后生成一个顶点位于中心点的理想圆锥，就是把一个穹隆或者凹坑理想化成一个正圆锥，我试了一下，只能单独地用matlab读取dem和绘制一个圆锥，不能把两者放在一起，让圆锥有dem的坐标系统。这要怎么写代码？
参考GPT和自己的思路：
你可以通过以下步骤来绘制一个顶点在中心的圆锥，并将其与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);
希望这能帮到你！
不知道你这个问题是否已经解决, 如果还没有解决的话: