Гонсалес Р., Вудс Р., Эддинс С. Цифровая обработка изображений в среде Matlab (2006), страница 103

' ); епб ГипсСзоп [з/С, А, РРС)= со[отдгасЦГ, Т) %СО[ ОКСКАЯ СогприСея СЬе чесСот дгаЖепС о[ ап КСВ зтпаде. / [ЧС, ЧА, РРС) = СОЬОНСНАО(Р, Т) сошриСев СЬе чессог Егаб1епС, ЧС, Х алб соггевропб1пЕ ал81е аггау, ЧА, (1п таб1алв) о1 НОВ АшаЕе Х Р 1С а1во сошриСев РРС, СЬе рег-р1але сошровАСе Егаб1епС '/ оЬСа1пеб Ьу вишш1пЯ СЬе 2-0 дгаб1ептв о1 СЬе 1пб1ч1биа1 со1ог '/ р1апев.

1приС Т зв а СЬгевЬо1б Ап СЬе галЕе [О, 1). П 1С Ав '/ 1пс1ибеб 1п СЬе атЕишепС 11зС, СЬе ча1иев о1 ЧС алб РРС аге '/ СЬтевЬо1беб Ьу 1еСС1пЯ ЧС(х,у) = 0 1ог ча1иез (= Т алб ЧС(х,у) '/. ЧС(х,у) оСЬегя1ве. 81ш11аг сошшептв арр1у Со РРС. 11 Т 1в поС '/ Апс1ибеб гп СЬе агЯишепС 11зС СЬеп Т 1в веС Со О. ВоСЬ оисриС '/ дгаб[епсв ате зса1еб Со СЬе галде [О, 1) . 11 (пб1шз(1) = 3) ! (з1яе(1, 3) = 3) еггог('1приС 1шаЕе шияС Ье НОВ.'); епб Х Соврите СЬе х алб у бег1чаС1чев о1 СЬе СЬгее сошропепС АшаЯев '/ ив1пЯ БоЬе1 орегатогв.

вЬ = 1врес1а1('воЬе1'); вч = вЬ'; Нх = 1шт11Сег(боиЬ1е(т(:,:, 1)), вЬ, 'гер11саСе'); Ву = 1ш$11Сез (боиЫе(1(:, з, 1)), зч, 'тер11сасе'); Сх = Аш111Сег(боиЫе(1(:, :, 2)), вЬ, 'гер11сате'); Су = 1шт11Сет(боиЫе(т(:, :, 2)), зч, 'гер1зсате'); Вх = 1шт11Сег(боиЫе(т(:, :, 3)), зЬ, 'тер11сате'); Ву = 1шт11Сег(боиЬ1е(т(:,:, 3)), вч, 'гер11сасе'); Х Сошрите СЬе рагашесегя о1 СЬе чессот Егаб1епС. Яхх = Нх."2 + Сх."2 + Вх.

2; Еуу = Ну."2 + Су."2 + Ву."2; Яху = Нх.еНу + Сх.*бу + Вх.*Ву; А = 0.8*(атал(2*дху./(Ехх — Еуу + ерв))); С1 = О.бе((Ехх + дуу) + (дхх — Еуу).асов(2еА) + 2едху.ез1п(2еА)); '/ Мою гереаС 1ог аия1е + РА/2. ТЬеп зе1есС СЬе шах1шиш аС еасЬ ро1пС. П--..., В БДЗ Х 1пХСХа11хе 1 ая а со1пшп чессог. 1С в111 Ье гевЬаре6 1атег / шяо ап ХшаЕе. 1 хетов(МеМ, 1); Т чагагЯ1п13т; ш = нагагЕ1п14); ш = ш(:)'; '/ МаЬе вате СЬаС ш Хя а гов чесСог. 1й 1епЕСЬ(чагагЕш) == 4 шеСЬо6 = 'епс116еап'; е1ве11 1епЕСЬ(чагагЕ1п) == 5 шеСЬо6 = 'шаЬа1апоЫв'; е1ве еггог('ИгопЯ шппЬег ой Хприсв.'); еп6 яч1СсЬ шеСЬо6 саве 'епс116еап' '/ Сошрпсе СЬе Епс116еап 61всвлсе Ьесвееп а11 тоня о1 Х ап6 ш.

Яее '/ ЯесСХоп 12.2 о1 01РОМ 1ог зл ехр1влаС1оп о1 СЬе 1о11ов1пЯ '/ ехргевв1оп. 0(1) 1в СЬе Епс1гбеал 61яСалсе ЬеСвееп чессог Х(1,:) '/ ал6 чесСог ш. р = 1епЯСЬ(т); О = яс[тС(зпш(аЬя(1 — тершаС(ш, р, 1))."2, 2)); саяе 'шаЬа1алоЫв' С = чагагя1п(Я; О = шаЬа1апоЫв(1, С, ш); оСЬегч1ве еггог('ОпЬповп веЕшепсаС1оп шеСЬо6.') еп6 '/ О 1в а честог о1 в1хе МН-Ьу-1 сопса1п1пЯ СЬе 61всалсе сошрисаСХопв Х Хтош а11 СЬе со1от РХхе1я Со чесСог ш.

Р1п6 СЬе 61вСелсев <= Т. 3 = 11п6(О <= Т); '/. Яес СЬе ча1пез о1 1(О) Со 1. ТЬеве аге СЬе яеЯшепсе6 '/. со1ог РХхе1в. 1(О) = 1; Х ВезЬаре 1 шсо ап М-Ьу-И ХшаЯе. 1 = гевЬаре(1, М, Н); тппсС[оп с = соппесСро)у(х, у) %СОТно[ЕСТРОХЛ' СоппесСв негероев о( а ро1удоп. Х С = СОММЕСТРОЬУ(Х, У) соплесСя СЬе ро1пся вХСЬ соог61пасев ЯХнеп Х 1п Х ап6 У в1СЬ ясга1ЯЬС 11пез. ТЬеве ро1псв аге авяише6 Со Ъе а '/ веопепсе о1 ро1уЯоп нетСХсез отЕап1хе6 ш СЬе с1осЬв1ве ог Х сочлсетс1осЬв1ве 61гесС1оп. ТЬе опсриС, С, 1в СЬе веС о1 ро1псв '/ а1опЕ СЬе Ьопп6ату о1 СЬе ро1уЯоп ш СЬе Хотш о1 ап пг-Ъу-2 Х соог61пасе вес)пепсе Хп СЬе ваше 61гесСХоп ав СЬе 1прпС. ТЬе 1авС '/ рошС 1п СЬе ее~[палое Хв ес[па1 Со СЬе 11гзС.

ч = [х(:), у(:)); % С1ове ро1уЕоп. 11 1вес(па1(ч(еп6,:), н(1,:)) ч(еп6 + 1, :) = ч(1, :); еп6 % Соппесс негСАсез. зеЕшепсв = се11(1, 1епЕСЬ(ч) — 1); Хог 1 = 2:1епЕСЬ(ч) (х, у) = 1пС11пе(ч(1 - 1, 1), ч(1, 1), н(1 - 1, 2), ч(1, 2)); веЕшепсзП вЂ” 1) = (х, у); еп6 с = саС(1, веЕшепгвСВ)); 1ппсСюп в = ЖатпеСег(Т ) %ЕИАМЕТЕгс Меаяиге 6!апьеСег апд ге!аСей ргорег11ея оГ илаЕе геЕЬвпз.

% Я = РТАИЕТЕК(Ь) сошрпсея СЬе 61ашегег, СЬе ша)ог ах1в еп6ро1пгз, % СЬе ш1пог ах1з еп6ро1псз, ал6 СЬе Ьаз1с гесСапЕ1е оХ еасЬ 1аЬе1е6 % геЕАоп гп СЬе 1аЬе1 шаСг1х 1. Ров1С1не 1пгеЕег е1ешепгя о1 1 % соггезроп6 Со 61ггегепС геЕ1опв. рог ехашр1е, СЬе зеС оХ е1ешепСя % оХ 1 ес(па1 Со 1 соггезроп6з Со геЕ1оп 1; СЬе веС о1 е1ешепгв о1 1 % ецпа1 Со 2 соггевроп6в Со геЕАоп 2; ап6 во оп.

Я Ав а вСгпсспге % аггау о1 1епЕСЬ шах(Ь(:)). ТЬе 11е16в оХ СЬе вггпсгпге аггау % 1пс1п6е: % % Р1ашесег % Иа)огАх1з % ИАпогАх1в % Вав1сКесгапЕ1е % % ТЬе Р1ашегег 11е16, а вса1аг, Ав СЬе шах1шпш 61зсапсе Ьесчееп алу % Сяо р1хе1з 1п СЬе соггевроп61пЕ геЕ1оп. % % ТЬе Иа)огАх1в 11е16 1з а 2-Ьу-2 шагг1х. ТЬе гоня сопСа1п СЬе гоч % ап6 со1пшп соог61пасез 1ог СЬе еп6ро1пСв оХ СЬе ша)ог ах1з ой СЬе % соггеяроп61пЕ гед1оп. % % ТЬе ИАпогАх1з 11е16 1я а 2-Ьу-2 шагг1х.

ТЬе гоня сопгагп СЬе гоч % ал6 со1пшп соог61пасез Хог СЬе епброгпсз оХ СЬе ш1пог ах1в о1 СЬе % соггевроп61пЕ геЕ1оп. % % ТЬе Ваз1сКессапЕ1е 11е16 Ав а 4-Ъу-2 шасг1х. ЕасЬ гоч сопса1пв % СЬе гоч ял6 со1ояш соог6гпасев о1 а согпег оХ СЬе % геЕ1оп-епс1ов1пЕ гесгялЕ1е 6е11пе6 Ьу СЬе ша)ог ап6 шгаог ахея.

П В 5В~~5~ % % Рог шоте 1пйогшаС1оп аЬоиС СЬеве шеазигешепСв, зее ЯесС1оп 11.2. 1 % о1 01Я1Са1 1шаде Ргосеввгпд, Ьу Оопха1ех ап6 'яоо6в, 2п6 е61С1оп, % РтепСАсе На11. в = геЕ1опргорв((., (' 1шаЕе', 'Вопп61пЕВох'г); 1ог Ь = 1: 1епЯСЬ(в) [в(Ь) .Р1ашесег, в(Ь) .Ма)огАх1в, рег1ш г, рет1ш с] сошриСе 61ашегег(в(Ь)); [в (К) .

Вав1сНесгяпЯ1е, в (К) . М1погАх1в] сошрисе Ьав1с гесгапЯ1е(в(Ь), рег1ш г, рег1ш с); еп6 '/ -% 1ппсС1оп [6, ша]огах1в, г, с] = сошрпге 61ашегег(в) % [О, МАООВАХТБ, В, С] = СОМРОТЕ ОХАМЕТЕН(Б) сошрпсев СЬе 61ашеСет '/ ап6 ша]ог ахгз 1от СЬе геЕАоп гергезепсе6 Ьу СЬе всгпсСпге Я. Я % шпвС сопСагп СЬе 11е16в 1шаде ап6 Вопп61пЯВох. СОМРОТЕ ОХАМЕТЕК '/ а1зо гегпгпв СЬе гоя ап6 со1пшп соог6тпаСез (К ап6 С) о1 СЬе % рег1шеСет РАхе1з о1 в.1шаде. '/ Сошрпсе гоя ап6 со1ошп соот6гпасез о1 рег1шесег РАхе1в. [г, с] = 11п6(Ьярет1ш(в.

1шаяе) ); г=г(:); с = с(:); [гр, ср] = ргипе р1хе1 11вС(т, с); пош р1хе1в = 1епЯСЬ(гр); вюАСсЬ ппш р1хе1з саве 0 6 = -1п1; ша]огах1в = опев(2, 2); саяе 1 6=0; ша)огах1в = [гр ср; гр ср]; саве 2 6 = (гр(2) — гр(1))"2 + (ср(2) - ср(1))"2; ша]огах1в = [гр ср]; оСЬегюгяе % Оепегаге а11 сошЬАпаС1опв ой 1:пош р1хе1в Са)сеп Сяо аС аС С1ше. % МесЬо6 впЯЯевсе6 Ьу Ресег АСЬ)аш. [16х(:, 2) 16х(:, 1)] = 11п6(Сг11(спев(ппш РАхе1з), -1)); гг = гр(16х); сс = ср(16х); 6гвС вс)паге6 = (гг(:, 1) — тт(:, 2)). 2 + (сс(:, 1) — сс(:, 2))."2; [шах п1яС вЧпаге6, 16х] = шах(61вС зчпаге6); ша]отах1в = [гг(16х,:)' сс(16х,:)']; 6 = вчгС(шах 6гвС вчиаге6); баррет 1шаяе тоя = в.Вопп61пяВох(2) е 0.5; ~~66 О В 1еХС 1шаЯе со1 = в.Яопп61пЯВох(1) + 0.5; ша)огах1в(:, 1) = ша)огах1в(:, 1) + пррег гшаЕе гов - 1; ша5огах1в(:, 2) = ша)огах1в(:, 2) + 1ейС 1шаЯе со1 - 1; епВ '/-— -/.

йшсС1оп [Ьав[стесС, ш1погах1в] = сошрисе Ьав[с тессапЯ1е(в, рег1ш г, рег[ш с) / [ВАЯТСНЕСТ,МТМОНАХТЯ] = СОМРОТЕ ВАЯТС НЕСТАИОЬЕ(Я, РЕН?М Н, Х РЕНТМ С) сошрисев СЬе Ьав1с гесСапЯ1е апй СЬе ш1пот ах[в '/ епд-ро1псв 1ог СЬе геЯ1оп тергевепсей Ьу СЬе вСгиссиге Я. Я шпвС '/ сопса[п СЬе 11е16в ТшаЯе, ВоппЖЬЕВох, Ма]огАх1в, апй '/ 01ашесег. РЕНТМ Н ап6 РЕНТМ С аге СЬе гов апй со1пшп соог61пасев '/ о1 рег[шесег ой в.1шаЯе.

ВАЯТСНЕСТ 1в а 4-Ьу-2 шаст[х, еасЬ тою '/ о1 вЬАсЬ сопСа1пв СЬе гов апд со1ишп соогВ1пасев ой опе сотпег о1 '/ СЬе Ъав1с тессапЯ1е. '/. Сошрпве СЬе ог[епсаС[оп о1 СЬе ша)от ахтв. СЬеса = аСап2(в.Ма)огАх[в(2, 1) — в.Ма]огАх[в(1, 1), в. Ма) огАх1в (2, 2) — в . Ма) огАх[в (1, 2) ); '/ Гогш гоСаС1оп шаСг[х. Т = [сов(СЬеса) в1п(СЬеса); -в1п(СЬеса) сов(СЬеСа)]; Х НоСасе рег1шеСег р1хе1в. р = [рег1ш с рет1ш г]; р = р * Т ; Х Са1сп1асе ш1п1шпш епй шах1шпш х- апй у-соог61патев 1ог СЬе тосасед '/ рег1шеСег р[хе1в. х=р(:, 1); у=р(;,2); ш1п х = ш1п(х); шах х = шах(х); ш1п у = ш1п(у),' шах у = шах(у); согпегв х = [ш1п х шах х шах х ш1п х]'; согпегв у = [шга у ш1п у шах у шах у]'; '/ НоСасе согпегв о1 СЬе Ьав1с гессапя1е. сотпегв = [согпегв„х сотпегв у] е Т; '/ Тгапв1асе ассог61пЯ Со СЬе геЕАоп'в Ьопп61пЯ Ьох. баррет 1шаЯе тою = в.Вопп61пЯВох(2) + 0.5; 1еХС 1шаЕе со1 = в.ВошЫ1пЯВох(1) + О.б; Ьав[стесС = [согпегв(:, 2) + баррет 1шаЯе гов — 1, согпегв(:, 1) + 1еХС АшаЯе со1 — 1]; Х Сошрпсе ш1пог ах1в епй-ро1псв, госасей.

х = (ш[п х + шах х) / 2; у1 = шгп у; у2 = шах у; епйро1псв = [х у1; х у2]; '/ НотаСе ш1пог ах[в епй-ро1псв ЬасН. П В Б~Д7 епйро1псв = еибро1псв е Т; % Ттапя1асе ассог81пЕ Со СЬе теЕйоп'в Ьоип61пЕ Ьох. ш1погахйв = [епйро1псв(:, 2) е иррег 1шаЕе тою — 1, епйро1псв(:, 1) е 1е1С 1шаЕе со1 — 1); %-— ХипсС1оп [г, с) = ргипе рйхе1 11вС(г, с) % [Е, С) = РЮМЕ РТХЕЕ ЕТЯТ(Е, С) гешочев р1хе1в Хгош СЬе чессогв % К ап6 С СЬаС саппоС Ье епбро1пся оК СЬе ша)ог ахйя.

ТЬ[в % е11ш1паС1оп 1в Ьавед оп ЕеошеСтйса1 сопясгатпСв Йезсг1Ьеб 1п % Нияз, ТшаЕе Ргосезв1пЕ НапйпооИ, СЬарсег 8. Сор = ш1п(г); ЬоССош = шах(т); Хе1С = ш1п(с); г18ЬС = шах(с); % 'вЬ1сЬ роЫСв аге ЫвЫе СЬе иррег сйгс1е? х = (1е1С е г1ЕЬС)/2; у = Сор; га81ив = ЬоССош — Сор; тпяЫе иррег = ( (с — х). 2 е (г - у). 2 ) < гад1ив 2; % еЬ1сЬ ро1пся аге 1пвЫе СЬе 1овег с1гс1е? у = Ьоссош; йпзЫе 1овег = ( (с — х)."2 е (г — у)."2 ) < гаЖив"2; '/ еЬ1сЬ ротпсв аге 1пвЫе СЬе 1е1С сйгс1е? х = Хе1С; у = (Сор е Ьоссош)/2; га81ив = гйЕЬС вЂ” 1еХС; 1пвЫе 1е1с = ( (с — х) ."2 е (г — у) ."2 ) < гаЫия"2; % 'еЬ1сЬ ро1пся аге 1пвЫе СЬе гйЕЬС сйгс1е? х = г1ЕЬС; 1пяЫе г1ЕЬС = ( (с — х). 2 е (г — у).

2 ) < гаЖия"2; % Е11штпасе ротпся СЬаС аге тпв16е а11 Тоит с1гс1ев. ое1есе Ых = 11пй(1пяЫе 1е1С й 1пвЫе т1ЕЬС й 1пвЫе иррег й 1пвЫе 1овет); гИе1есе Ых) = [); сИе1есе Ых) = [); [ипсС(оп с = [сЬсос[е[Ь, оопп, й[г) %ЕСНСО[)Е СотриСея СЬе г1еетпап сЬа[п соде оГ а Ьоипдагу.

% С = РСНСОВЕ(В) сошрисев СЬе 8-соппессей Ртеешап сЬа1п сойе о1 а % веС о1 2-0 соог61пасе ра1гя сопсайпей 1п В, ап пр-Ьу-2 аггау. С % 1в а вСгиссиге ч[СЬ СЬе Хо11ою1пЕ 11еЫв: / с.Хсс = Ртеешеп сЬайп сойе (1-Ьу-пр) Х с.д1тт = ГАтвс д111етепсе о1 соде с.1сс (1-Ьу-пр) Х с.шш = 1псеНег о1 штпйшиш шаЕп1сиде 1тош с.1сс (1-Ьу-пр) Х с.д1ттшш = РйтзС д1ттетепсе о1 соде с.шш (1-Ьу-пр) Х с.хОуО = СоотдйпаСев вЬете СЬе соде зСаття (1-Ъу-2) Х Х С = РСНСОРЕ(В, СОИМ) ртодисев Сье ваше оитритз ав аЬоче, ЬиС Х в1СЬ СЬе соде соппесс1ч1су зрес111ед 1п СОИМ. СОНИ сап Ье 8 1от Х ап 8-соппесСед сЬа1п соде, от СОИМ сал Ье 4 1от а 4-соплессед Х сЬа1п соде. Ярес11у1пН СОИМ=4 Ав на1тд оп1у 11 СЬе 1приС / яеоиепсе, В, сопса1пя Сталв1С1опв юйсЬ ча1иев О, 2, 4, апд б, Х ехс1ивйче1у.