霍夫变换提取机场跑道区域的matlab代码

以下是使用霍夫变换提取机场跑道区域的 Matlab 代码：

%% 读取图像并进行预处理
im = imread('airport.jpg'); % 读取图像
im_gray = rgb2gray(im); % 转换为灰度图像
im_edge = edge(im_gray, 'canny'); % 进行边缘检测

%% 进行霍夫变换
[H,theta,rho] = hough(im_edge); % 进行霍夫变换，得到霍夫空间 H、极角 theta 和极径 rho
P = houghpeaks(H, 5); % 在霍夫空间中寻找峰值点
lines = houghlines(im_edge,theta,rho,P,'FillGap',20,'MinLength',50); % 根据峰值点提取直线段

%% 绘制直线段
figure, imshow(im), hold on
for k = 1:length(lines)
   xy = [lines(k).point1; lines(k).point2];
   plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');
end

%% 提取跑道区域
% 首先找出直线段中与 x 轴平行的直线
hor_lines = [];
for k = 1:length(lines)
    theta = lines(k).theta;
    if (theta > 85 && theta < 95) || (theta > -95 && theta < -85)
        hor_lines = [hor_lines; lines(k)];
    end
end

% 然后计算这些直线段的交点，得到跑道区域
corners = [];
for i = 1:length(hor_lines)-1
    for j = i+1:length(hor_lines)
        p1 = hor_lines(i).point1;
        p2 = hor_lines(i).point2;
        p3 = hor_lines(j).point1;
        p4 = hor_lines(j).point2;
        
        % 计算交点
        x = (p1(2) - p3(2) + (p3(1)-p1(1))*(p4(2)-p3(2))/(p4(1)-p3(1))) / ((p2(2)-p1(2))/(p2(1)-p1(1)) - (p4(2)-p3(2))/(p4(1)-p3(1)));
        y = (x - p1(1))*(p2(2)-p1(2))/(p2(1)-p1(1)) + p1(2);
        corners = [corners; x y];
    end
end

% 绘制跑道区域
line([corners(1,1) corners(2,1)], [corners(1,2) corners(2,2)], 'LineWidth', 2, 'Color', 'red');
line([corners(3,1) corners(4,1)], [corners(3,2) corners(4,2)], 'LineWidth', 2, 'Color', 'red');

注意：此代码仅提供参考，具体实现可能需要根据实际情况进行调整。