Пояснить реализацию алгоритма отсечения

@toypypp · Регистрация: 04.10.2016

Студворк — интернет-сервис помощи студентам

Добавьте,пожалуйста, немного поясняющих комментариев, что делает каждый из блоков. Нужно разобраться для доклада.

Кликните здесь для просмотра всего текста

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
n = 100;
Xmin = -4; Xmax = 4;
Ymin = -3; Ymax = 3;
 
figure; hold on;
 
P1x = rand(1,n)*20 - 10;  
P2x = rand(1,n)*20 - 10;
P1y = rand(1,n)*20 - 10;  
P2y = rand(1,n)*20 - 10;
P1code = zeros(n,4);
P2code = zeros(n,4);
 
for i = 1:n    
 
    if P1x(i) < Xmin
        P1code(i,1) = 1;
    end
    if P1x(i) > Xmax
        P1code(i,2) = 1;
    end
    if P1y(i) < Ymin
        P1code(i,3) = 1;
    end
    if P1y(i) > Ymax
        P1code(i,4) = 1;
    end    
 
    if P2x(i) < Xmin
        P2code(i,1) = 1;
    end
    if P2x(i) > Xmax
        P2code(i,2) = 1;
    end
    if P2y(i) < Ymin
        P2code(i,3) = 1;
    end
    if P2y(i) > Ymax
        P2code(i,4) = 1;
    end
 
    plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'b-');
end
hold off;
 
P_label = zeros(1,n);
figure; hold on;
for i = 1:n
    P_or = P1code(i,1:4) | P2code(i,1:4);
    if sum(P_or) == 0
        P_label(i) = 1;
        plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'r-');
    end
    P_and = P1code(i,1:4) & P2code(i,1:4);
    if sum(P_and) > 0
        P_label(i) = 2;
        plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'g-');
    end
end
hold off;
 
figure; hold on;
for i = 1:n
    if P_label(i) == 0
 
        P_or = P1code(i,1:4) | P2code(i,1:4);
        plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'g-');
 
        if P_or(1) == 1
            Py = P1y(i) + (Xmin-P1x(i))*(P2y(i)-P1y(i))/(P2x(i)-P1x(i));
            if P1x(i) < Xmin
                P1x(i) = Xmin; P1y(i) = Py;
            elseif P2x(i) < Xmin
                P2x(i) = Xmin; P2y(i) = Py;
            end            
        end
 
        if P_or(2) == 1
            Py = P1y(i) + (Xmax-P1x(i))*(P2y(i)-P1y(i))/(P2x(i)-P1x(i));
            if P1x(i) > Xmax
                P1x(i) = Xmax; P1y(i) = Py;
            elseif P2x(i) > Xmax
                P2x(i) = Xmax; P2y(i) = Py;
            end            
        end
 
        if P_or(3) == 1
            Px = P1x(i) + (Ymin-P1y(i))/(P2y(i)-P1y(i))*(P2x(i)-P1x(i));
            if P1y(i) < Ymin
                P1x(i) = Px; P1y(i) = Ymin;
            elseif P2y(i) < Ymin
                P2x(i) = Px; P2y(i) = Ymin;
            end            
        end
 
        if P_or(4) == 1
            Px = P1x(i) + (Ymax-P1y(i))/(P2y(i)-P1y(i))*(P2x(i)-P1x(i));
            if P1y(i) > Ymax
                P1x(i) = Px; P1y(i) = Ymax;
            elseif P2y(i) > Ymax
                P2x(i) = Px; P2y(i) = Ymax;
            end            
        end
        if P1x(i) >= Xmin & P1x(i) <= Xmax & ...
                P2x(i) >= Xmin & P2x(i) <= Xmax & ...
                P1y(i) >= Ymin & P1y(i) <= Ymax & ...
                P2y(i) >= Ymin & P2y(i) <= Ymax
            plot([P1x(i),P2x(i)],[P1y(i),P2y(i)],'m-','LineWidth',2);
        end
        plot([-4 4 4 -4 -4],[-3 -3 3 3 -3],'b-','LineWidth',2);
    end
end
hold off;

https://blog.csdn.net/soulmeetliang/article/details/79179350

@nuHrBuH · 14.05.2019, 19:42

Насколько хорошо должен работать код?

@AlMih · 14.05.2019, 21:21

toypypp, Правильно я понимаю, что ты где-то стырил код и хочешь его где-то представить как свой в виде доклада, но не понимаешь как он работает?

@toypypp · 16.05.2019, 13:40 **[ТС]**

Ну я даже ссылку указал откуда взял

@nuHrBuH · 16.05.2019, 14:54

Я сделал

Кликните здесь для просмотра всего текста

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
clc,clear
 
n = 20;                    % количество отрезков
    % границы области
% t = linspace(0,2*pi,5)+pi/4;
% xv = 8*cos(t);
% yv = 6*sin(t);
xv = [-4 -4 4 4 -4];
yv = [-3 3 3 -3 -3];
 
    % генерируем точки отрезков в границах -10 : 10
a = 20;
x1 = rand(1,n)*a - a/2;  
x2 = rand(1,n)*a - a/2;
y1 = rand(1,n)*a - a/2;  
y2 = rand(1,n)*a - a/2;
 
    % точки пересечения отрезков и границ области
CR = false(1,n);
[xc1,xc2,yc1,yc2] = deal( zeros(1,n) );
for i = 1:n
        % отрезок
    x = [x1(i) x2(i)];
    y = [y1(i) y2(i)];
    [xc,yc] = polyxpoly(x,y,xv,yv);
    if ~isempty(xc)            % если точки пересечения есть
        CR(i) = true;
        xc1(i) = xc(1);        % первая точка пересечения
        yc1(i) = yc(1);
        if length(xc) == 2     % если точек - две
            xc2(i) = xc(2);
            yc2(i) = yc(2);
        else                   % если точка одна, определяем 
                               % какая из точек отрезка 
                               % находится внутри области
            if inpolygon(x1(i),y1(i),xv,yv)
                xc2(i) = x1(i);
                yc2(i) = y1(i);
            else
                xc2(i) = x2(i);
                yc2(i) = y2(i);
            end
        end
    end
end
 
in1 = inpolygon(x1,y1,xv,yv);
in2 = inpolygon(x2,y2,xv,yv);
IN  = in1 & in2;                % отрезки находяющиемя внутри
OUT = ~(IN | CR);            % отрезки находяющиемя снаружи
 
plot(xv, yv, 'b-','LineWidth',2);                   % границы области
hold on
plot([x1; x2], [y1; y2],'g')                        % все отрезки
plot([xc1(CR); xc2(CR)], [yc1(CR); yc2(CR)],'r')    % пересекающие
plot([x1(IN);   x2(IN)], [y1(IN);   y2(IN)],'k')    % внутри
plot([x1(OUT); x2(OUT)], [y1(OUT); y2(OUT)],'b')    % снаружи
hold off

@toypypp · 18.05.2019, 10:32 **[ТС]**

nuHrBuH, спасибо большое! можете прислать результат? Не могу запустить без mapping toolbox, а в пробной версии matlab не устанавливается.

@nuHrBuH · 18.05.2019, 13:40

Странно, это из-за inpolygon() и polyxpoly()?

Измененный код:

Кликните здесь для просмотра всего текста

Matlab M
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
clc,clear
 
n = 50;                    % количество отрезков
    % границы области
% t = linspace(0,2*pi,5)+pi/4;
% xv = 8*cos(t);
% yv = 6*sin(t);
xv = [-4 -4 4 4 -4]';
yv = [-3 3 3 -3 -3]';
 
    % генерируем точки отрезков в границах -10 : 10
a = 20;
x1 = rand(n,1)*a - a/2;  
x2 = rand(n,1)*a - a/2;
y1 = rand(n,1)*a - a/2;  
y2 = rand(n,1)*a - a/2;
 
    % точки пересечения отрезков и границ области
CR = false(n,1);
[xc1,xc2,yc1,yc2] = deal( zeros(n,1) );
for i = 1:n
        % отрезок
    x = [x1(i) x2(i)];
    y = [y1(i) y2(i)];
    [xc,yc] = intersections(x,y,xv,yv);
    if ~isempty(xc)            % если точки пересечения есть
        CR(i) = true;
        xc1(i) = xc(1);        % первая точка пересечения
        yc1(i) = yc(1);
        if length(xc) == 2     % если точек - две
            xc2(i) = xc(2);
            yc2(i) = yc(2);
        else                   % если точка одна, определяем 
                               % какая из точек отрезка 
                               % находится внутри области
            if inpoly2([x1(i) y1(i)],[xv yv])
                xc2(i) = x1(i);
                yc2(i) = y1(i);
            else
                xc2(i) = x2(i);
                yc2(i) = y2(i);
            end
        end
    end
end
 
in1 = inpoly2([x1 y1],[xv yv]);
in2 = inpoly2([x2 y2],[xv yv]);
IN  = in1 & in2;                % отрезки находяющиемя внутри
OUT = ~(IN | CR);            % отрезки находяющиемя снаружи
 
plot(xv, yv, 'b-','LineWidth',2);                   % границы области
hold on
plot([x1 x2]', [y1 y2]','g')                        % все отрезки
plot([xc1(CR) xc2(CR)]', [yc1(CR) yc2(CR)]','r')    % пересекающие
plot([x1(IN)   x2(IN)]', [y1(IN)   y2(IN)]','k')    % внутри
plot([x1(OUT) x2(OUT)]', [y1(OUT) y2(OUT)]','b')    % снаружи
hold off

Скачать аналог inpolygon(): inpoly2()
Скачать аналог polyxpoly(): intersections()

Добавлено через 1 минуту
Медленно считает, што поделать

Новые блоги и статьи Все статьи Все блоги /
SDL3 для Desktop (MinGW): Создаём пустое окно с нуля для 2D-графики на SDL3, Си и C++ 8Observer8 10.03.2026 Содержание блога Финальные проекты на Си и на C++: hello-sdl3-c. zip hello-sdl3-cpp. zip Результат:	Установка CMake и MinGW 13.1 для сборки С и C++ приложений из консоли и из Qt Creator в EXE 8Observer8 10.03.2026 Содержание блога MinGW - это коллекция инструментов для сборки приложений в EXE. CMake - это система сборки приложений. Здесь описаны базовые шаги для старта программирования с помощью CMake и. . .	Как дизайн сайта влияет на конверсию: 7 решений, которые реально повышают заявки Neotwalker 08.03.2026 Многие до сих пор воспринимают дизайн сайта как “красивую оболочку”. На практике всё иначе: дизайн напрямую влияет на то, оставит человек заявку или уйдёт через несколько секунд. Даже если у вас. . .	Модульная разработка через nuget packages DevAlt 07.03.2026 Сложившийся в . Net-среде способ разработки чаще всего предполагает монорепозиторий в котором находятся все исходники. При создании нового решения, мы просто добавляем нужные проекты и имеем. . .
Модульный подход на примере F# DevAlt 06.03.2026 В блоге дяди Боба наткнулся на такое определение: В этой книге («Подход, основанный на вариантах использования») Ивар утверждает, что архитектура программного обеспечения — это структуры,. . .	Управление камерой с помощью скрипта OrbitControls.js на Three.js: Вращение, зум и панорамирование 8Observer8 05.03.2026 Содержание блога Финальная демка в браузере работает на Desktop и мобильных браузерах. Итоговый код: orbit-controls-threejs-js. zip. Сканируйте QR-код на мобильном. Вращайте камеру одним пальцем,. . .	SDL3 для Web (WebAssembly): Синхронизация спрайтов SDL3 и тел Box2D 8Observer8 04.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-sync-physics-sprites-sdl3-c. zip На первой гифке отладочные линии отключены, а на второй включены:. . .	SDL3 для Web (WebAssembly): Идентификация объектов на Box2D v3 - использование userData и событий коллизий 8Observer8 02.03.2026 Содержание блога Финальная демка в браузере. Итоговый код: finish-collision-events-sdl3-c. zip Сканируйте QR-код на мобильном и вы увидите, что появится джойстик для управления главным героем. . . .

@nuHrBuH 483 / 427 / 205 Регистрация: 04.03.2011 Сообщений: 1,259
	14.05.2019, 19:42
	Насколько хорошо должен работать код? 0

@AlMih 416 / 200 / 69 Регистрация: 20.01.2019 Сообщений: 714
	14.05.2019, 21:21
	toypypp, Правильно я понимаю, что ты где-то стырил код и хочешь его где-то представить как свой в виде доклада, но не понимаешь как он работает? 0

@toypypp 0 / 0 / 6 Регистрация: 04.10.2016 Сообщений: 110
	16.05.2019, 13:40 [ТС]
	Ну я даже ссылку указал откуда взял 0

@toypypp 0 / 0 / 6 Регистрация: 04.10.2016 Сообщений: 110
	18.05.2019, 10:32 [ТС]
	nuHrBuH, спасибо большое! можете прислать результат? Не могу запустить без mapping toolbox, а в пробной версии matlab не устанавливается. 0

Пояснить реализацию алгоритма отсечения

Решение