使用DEM将图片提取单位坡度的matlab代码

在MATLAB中，你可以使用Digital Elevation Model (DEM)数据来计算单位坡度。以下是一个简单的示例代码片段，它假设你已经有了一个DEM矩阵(`dem`)和对应的行向量(`row`)和列向量(`col`)：

% 加载或创建DEM数据
% dem = imread('your_dem_image.tif'); % 替换为实际的图像路径
% 或者如果你已经有一个数字矩阵，直接使用 dem

% 计算每个像素点的水平和垂直梯度
[row, col] = meshgrid(1:size(dem,1), 1:size(dem,2));
dx = diff(dem, 0, 2); % 水平梯度（右 - 左）
dy = diff(dem, 0, 1); % 垂直梯度（下 - 上）

% 计算单位坡度 (如果需要对齐为正北)
unit_slope = sign(dy);
slope_angle_rad = atan2(sign(dx), unit_slope); % 角度转换到弧度，使得0对应于北方
slope_angle_deg = degrees(slope_angle_rad); % 转换成度数

% 如果你想得到单位坡度，可以取绝对值并除以最大值
unit_slope_normalized = abs(slope_angle_deg) ./ max(abs(slope_angle_deg));

% 查看第一个像素的坡度信息
disp(['First pixel slope angle: ', num2str(unit_slope_normalized(1,1)) ' degrees']);

%