MATLAB用canny求得图像轮廓之后，根据列将轮廓分为左右两边，分布求左右两边曲率变化最大处

要将轮廓分为左右两边，并找到左右两边曲率变化最大处，可以使用以下代码：

% 读取图像
img = imread('example.png');
% 边缘检测
img_edge = edge(img, 'canny');
% 提取轮廓
[B, L] = bwboundaries(img_edge, 'noholes');
% 取第一个轮廓
boundary = B{1};
% 计算曲率
k = curvature(boundary(:,2), boundary(:,1));
% 找到曲率变化最大的位置
[~, max_idx] = max(abs(diff(k)));
max_point = boundary(max_idx, :);
% 将轮廓分成左右两边
left_boundary = boundary(boundary(:,2) < max_point(2), :);
right_boundary = boundary(boundary(:,2) >= max_point(2), :);
% 计算左右两边的曲率
k_left = curvature(left_boundary(:,2), left_boundary(:,1));
k_right = curvature(right_boundary(:,2), right_boundary(:,1));
% 找到左右两边曲率变化最大的位置
[~, max_idx_left] = max(abs(diff(k_left)));
[~, max_idx_right] = max(abs(diff(k_right)));
max_point_left = left_boundary(max_idx_left, :);
max_point_right = right_boundary(max_idx_right, :);

这段代码假设你已经读取了一张图像，然后进行了边缘检测和轮廓提取。然后它计算了轮廓的曲率，并找到了曲率变化最大的位置。接下来，它将轮廓分成左右两边，并计算了左右两边的曲率。最后，它找到了左右两边曲率变化最大的位置。最后的输出是两个二元素向量 max_point_left 和 max_point_right，它们分别包含了左右两边曲率变化最大的位置的 x 和 y 坐标。