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

х = х(:)„ у=у(:); ( (е.еса30 + е.еса12) 2 — Зв(е.еса21 + е.еса03) 2 ) (Зве.еса12 — е.еса30) в (е.еСа21 + е.еСа03) ( Зв(е.еса30 + е.еса12) "2 — (е.еса21 + е.еса03) "2 ); ГипсС1оп [х, у~ = пипрегро[у(В, се11з[зе) %МТХРЕНРОЬУ СогприСея СЬе ппшпппп регппеСег ро[удоп. % [Х, У) = МТНРЕВРОЬУ(Р, СЕЬЬЯ12Е) сошриСез СЬе чегС1сев 1п [Х, У) % оу СЬе шйп1шиш регйшесег ро1уЕоп оу а з1пЯ1е Ь1пагу геЕ1оп ог % Ьоипйагу 1п 1шаЕе В. ТЬе ргосейиге 1в Ьавей оп 31апвЬу'в % зЬг1пЫпЕ гиЬЬег Ьапй арргоасЬ.

Рагашесег СЕЬЬЯ12Е йесегш1пев СЬе % я1зе о1 СЬе в~[пахе сеПв СЬаС епс1озе СЬе Ьоипйагу о1 СЬе геЯ1оп % 1п В. СЕЬЬБ12Е шивС Ье а попхего 1псеЕег Егеасег СЬап 1. % % ТЬе а1Еог1СЬш 1я арр11саЫе оп1у Со Ьоипйагйев СЬаС аге поС % ве1г-1псегвесС1пЕ апй СЬаС йо поС Ьаче опе-рйхе1-СЫсЬ % ргосгив1опв. 11 сеПз1хе <= 1 еггог('СЕЕЬБ12Е шизС Ье ап 1псеяег > 1.'); епй % РШ В гп саяе СЬе 1приС чая ргонййей ая а Ьоипйагу. Ьасег % сье Ьоипйагу чгП Ье ехсгассей ч1СЬ 4-соппесС1чйСу, чЫсЬ % [в ге~)и1гей Ьу СЬе а1Еог1СЬш. ТЬе иве оу Ьчрег1ш авзигев % СЬаС 4-соппесСйчйсу 1в ргевегчей аС СЫв ро1пС.

В = йш11П(В, 'Ьо1ев'); В = Ьчрегйш(В); [М, М) = вйхе(В); % Тпсгеаве 1шаЯе в1хе во СЬаС СЬе ХшаЕе йв о1 в1хе К-Ьу-К % ч1СЬ (а) К >= шах(М,И) апй (Ь) К/сеПя1хе = а рочег о1 2. К = пехсроч2(шах(М, И)/сеПвйхе); К = (2 К)всеПв1хе; % Хпсгеаве 1шаЕе в1зе Со пеагевС хпСеЯег рочег о1 2, Ьу % аррепй1пЯ хегов Со СЬе епй о1 СЬе йшаЕе.

ТЫз чйП аПоч % «[иайсгее йесошроя1С[опз ав вшаП ав сеПв о1 в1хе 2-Ьу-2, % чЫсЬ 1в СЬе вшаПевС аПочей ча1ие о1 сеПзйхе, М=К-М; И=К-И; В = райаггау(В, [М М), 'розС'); % 1 1в поч о1 в1хе К-Ьу-К % ЦиайСгее йесошров1С1оп. 0 = ~[Сйесошр(В, О, сеПвйхе); % СеС аП СЬе виЫшаЯев оу в1яе сеПвйхе-Ьу-сеПв1хе. [ча1в, г, с) = ЧСЕеСЫЬ(В, Ц, сеПв1хе); % Сес аП СЬе яиЬйшаЯев СЬаС сопСа1п аС Теавс опе Ь1асЬ П В ° 9~9 % апд во оп. ТЬе 1авС е1ешепС оХ сЬе соде 1в СЬе Сгапв1С1оп Ьесчееп '/ СЬе 1авС апд 1вС е1ешепгз ой В. ТЬе е1ешепСв оХ В йогш а сю % веииепсе (вее аЪоче), зо че изе 'ваше' 1от СЬе ддгесС1оп гп % йшсС1оп ХСЬсоде. соде = 1сЬсоде(В, 4, 'ваше'); соде = соде.1сс; % Го11оч сЬе соде весиепсе со ехсгасс сЬе В1асК 0осв, В0, (сопчех % согпегз) апд 'еЬ1Се РоСз, 'н0, (сопсаче согпегв).

ТЬе Сгапв1С1опв аге % ав 1о11онв: О-Со-1=н0; О-Со-3=ВО; 1-Со-О=ВО; 1-Со-2=И0; 2-Со-1=ВО; '/ 2-Со-З=в0; 3-Со-О='н0; 3-Со-2=доС. ТЬе 1огши1а С = 2>11гвС вЂ” зесопд % х1чев СЬе Хо11оч1пз ипгоие на1иев 1ог СЬеве Сгапв1С1опз; -1, -3, 2, % О, 3, 1, б, 4. ТЬеве аге арр11саЬ1е со Сгаче1 тп СЬе сю д1гесгдоп. % ТЬе И0'в аге дгвр1асед опе-Ьа11 а д1азопа1 1гош СЬе В0'в Со йогш % СЬе Ьа11-се11 ехрапв1оп ге~(и1тед Тп СЬе а1зог1СЬш.

% ЧегС1сез ч111 Ье сошрисед ав аггау >чеггдсев> ой д1шепвдоп пн-Ьу-З, % нЬеге пч 1в СЬе пишЬег о1 негС1сев. ТЬе 11гвС Сно е1ешепСв ой елу % гоч о1 агтау негС1сев аге СЬе (х,у) соогд1паСев ой СЬе чегсех % соггевропд1пз Со СЬаС гоч, апд СЬе СЬггд е1ешепС гв 1 11 СЬе % чегСех 1з сопчех (В0) ог 2 1Х 1С 1в сопсане (И0). ТЬе 11гвС четСех '/ 1в Кпонп Со Ье сопчех, во 1С 1в ЫасК. негС1сез = [х1, у1, 1],' и=1; К=1; Хог К = 2т1епзСЬ(соде) 11 соде(К вЂ” 1) = соде(К) и = и+ 1; С = 2>соде(К-1) — соде(К); % С = ча1ие оХ Хогши1а. 11 С == -3 ~ С == г ~ С == 3 ) С == 4 % Сопчех: В1асК 0осв. чегС1сез(п, 1:3) = [х(К), у(К), 1]; е1ве11 С == -1 ( С == О ) С == 1 ( С == В '/ Оопсаче: МЬ1Се 0огз. 1Х С == -1 чегС1сез(п, 1:3) = [х(К) — сеПз1хе, у(К) — се11вгхе,2]; е1зе11 С==О чегС1сев(п, 1:3) = [х(К) + се11в1хе, у(К) - се11з1хе,2]; е1веН С =1 четС1сев(п, 1:3) = [х(К) + се11в1ге, у(К) + сеПв1хе,2]; е1зе негС1сев(п, 1:3) = [х(К) — се11в1хе, у(К) + се11в1хе,2]; епд е1ве '/ ИоСЬ1пВ Со до Ьете.

епд епд епд % ТЬе гевС оХ ш1прегро1у.ш ргосеввев СЬе чегС1сев Со % агтдче аС СЬе ИРР. (600 и В П1аК = 1' чЬт1е 11ав / Оесетш1пе юЬХсЬ четСХсев 11е оп от 1пв1йе СЬе % ро1у8оп чЬове четС1сев аге СЬе 81асК Оосв. Ое1есе а11 Х оСЬег ро1псв. 1 = 11пб(четС1сев(:, 3) == 1); хч = четСХсев(1, 1); '/ Соотд1пасев оХ СЬе В1асК Оосв. уч = четС1сев(1, 2); Х = четС1сев(:, 1); '/ Соот61пасев оХ а11 четС1сев. У четС1сев(:, 2); 1И = 1про1у8оп(Х, У, хч, уч); 1 = 11пй(1И = 0); четСХсев = четСХсев(1, :); '/ Моч сЬесК азот елу В1асК Посв СЬаС шау Ьаче Ьееп Ситпей 1псо '/ сопсаче четС1сев аХСег СЬе ртечХоив йе1есХоп всер. Ое1ете '/ апу висЬ В1асК Оосв ал6 тесошрисе СЬе ро1у8оп ав та СЬе '/ ртеч1оив весС1оп ой сойе.

ЫЬеп по шоте сЬелзев оссит, веС '/ Х1а8 Со О, чЬХсЬ салаев СЬе 1оор Со Сетш1пате. х = четС1сев(:, 1); у = нетСХсев(:, 2); апд1ев = ро1уапд1ев(х, у); '/ РХпй а11 СЬе таСет1от ап81ев. 1 = 11пй(ал81ев > 180 й нетСХсев(:, 3) == 1); 1Х ХвешрСу(1) 11а8 = 0; е1ве 3 = 1:1епВСЬ(нетС1сев); 1от К = 1:1еиВСЬ(1) К = Нп8(,) = 1(К)); Л =,)(К); епй четСХсев = четСХсев(3, :); еж$ еп6 '/ Ртпа1 рава Со йе1есе СЬе нетСХсев ч1СЬ ап81ев оу 180 йе8теев.

х = нетС1сев(:, 1); у = нетСХсев(:, 2); аи81ев = ро1уап81ев(х, у); 1 = Х1пб(ап81ев = 180); Х ЧетСХсев оХ СЬе ИРР: х = четСХсев(1, 1); у = четСХсев(1, 2); щ 8 6~~~~1 ГппсС1оп  — р1хе[йпр(А, тп, п) %Р1ХЕ[ 01)Р ОирйсаСев р1хе[в от" ап ипаце 1п ЬоСЬ сИгесС)опя. % В = Р1ХЕЫ()Р(А, И, М) оир11саСев еасЬ ргхе1 о1 А И С1шея 1п СЬе % нетс1са1 61гесс1оп апд И САшея Ап СЬе Ьог1гопса1 Ыгесс1оп. % Рагашесегв И апб М шияС Ье АпСеЕетя. П И 1я пос 1пс1иоей, Ас % бе1аи1ся со И. % СЬесК 1присв. 11 патЕ1п ( 2 егтог('Ас 1еаяс сюо Априся ате тес(и1геб.'); епо гЕ пагЕАп == 2 п=ш; епй % Сепегаге а несСот в1СЬ е1ешепСз 1:в1ге(А, 1) .

и = 1:вгге(А, 1); % 0ир11сасе еасЬ е1ешепс о1 сЬе нессот ш с1шев. ш = гоипб(ш); % Ргосесс аЕагпвс попгпсегЕегя. и = и(опев(1, ш), :); и = и(:); % Моя гереас 1от СЬе оСЬег й1гесс1оп. и = 1:зсге(А, 2); и = гоип6(п); и = н(опея(1, п), :); и = и(:); В = А(и, н); Гипс11оп апЕ1ея =. ро1уапЕ[ев[х, у) %РО?ХАХС? ЕЯ Сотприсев 1пСетпа1 ро1удоп апд[ев. % АИСЬЕЯ = РОЬУАИСЬЕЯ(Х, У) сошрисев СЬе гпсетгот апЕ1ев (1п % йеЕгеев) о1 ап агЬАСгагу ро1уЕоп нЬове негсгсея аге ЕАнеп Ап % [Х, У], огйегей 1п а с1осхнгве шаппег. ТЬе ртоЯгаш е11шгпаСез % бир11сасе аб]асепс гоня 1п [Х У], ехсерС СЬаС СЬе 11твС гон шау % еииа1 СЬе 1авС, во СЬас СЬе ро1уЕоп 1в с1овеб. % Рге11шспаг1ез.

[х у] = йирЕопе(х, у); % Е11ш1паСе бир11сасе негс1сез. ху = [х(:) у(:)]; 11 1вешрСу(ху) % Мо негс1сев! апЕ1ея = гегоя(0, 1); гесигп; еаза 11 ягге(ху, 1) == 1 ~ 1яе~[иа1(ху(1, :), ху(епб, :)) % С1ове СЬе ро1уЕоп ху(епд в 1, ;) = ху(1, :); еш( С602 и В % Ргесошрисе воша с)иапС1С1ев. 6 = 6111(ху, 1); ч1 = -6(1:епй, :); ч2 = [6(2:епй, :); 6(1, :)); ч1 йоС ч2 = виш(ч1 .е ч2, 2); шад н1 = загс(виш(ч1."2, 2)); шах ч2 = зцтс(виш(н2. 2, 2)); % РгоСесС ава1пяС пеаг1у йир11сасе негСАсев; оиСриС апв1е ч111 Ье 90 % йехтеев Хог висЬ савев. ТЬе сгеа1> ТигсЬег ртосессв а8а1пвс % ровв1Ь1е зша11 1шац1пагу апя1е сошропепсв 1п СЬове савве.

шах ч1( шах н1) = ерв; шах ч2( шах ч2) = еря; апв1ев = геа1(асов(ч1 йоС н2 ./ шах н1 ./ ша8 ч2) е 180 / р1); % ТЬе 11гвС апх1е сошрисей чая 1ог СЬе весопй чегсех, апй СЬе % 1азС чав 1ог СЬе 11гзС чегСех. Бсто11 опе ровАСАоп йочп Со % ша1се СЬе 1авС четсех Ье СЬе 11гвС. апв1ев = сйгсвЬ11С(апд1ев, [1, 0]); % Боч йесегшйпе 11 алу чегсйсев ате сопсаче аий ай)ивС СЬе апх1ев % ассогй1пв1у. ввп = сопчех апв1е Севс(ху); % Апу е1ешепс ой в8п сЬас'в -1 1п61сасев сЬас сЬе апв1е 1в % сопсане.

ТЬе соггеяропйгп8 апв1ев Ьане Со Ье зиЬСтасСей % 1гош 360. 1 = 11пй(зяп == -1)' апх1ев(1) = 380 — апв1ев(1); /- -/ ХипсС1оп в8п = сопчех апв1е Севс(ху) % ТЬе гочв о1 агтау ху аге отйетей четсгсев о1 а ро1увоп. 11 СЬе % ЬСЬ апв1е 1в сопнех ()О вл6 <= 180 йе8тевв) СЬеп ввп(Ь) % 1. ОСЬегч1ве ввп(К) = -1. ТЬйв 1ипссйоп аввишев СЬаС СЬе ХггвС % нетсех гп СЬе 11вС 1з сопнех, апй СЬаС по оСЬет чегсех Ьав а % вша11ег ча1ие оХ х-соотй1пасе. ТЬеяе Счо сопйАС1опв аге Стае гп % сЬе 11твС чегсех вепетасей Ьу сЬе ИРР а1вог1СЬш.

А1яо СЬе % чегС1сев аге авзишей Со Ье огйегей 1п а с1осЬч1ве вес(иепсе, впй % СЬеге сап Ье по йир11сасе чегсгсев. у % ТЬе Севс Ав Ьавей оп СЬе Тасс СЬаС ечегу сопчех негсех Ая оп СЬе % ровАСАне вгйе о1 СЬе 1гпе разв1пв СЬгоивЬ СЬе Счо чегсйсев % Ашшеййасе1у То11оч1пв еасЬ чегСех Ье1пх сопв16егей. 11 а чегСех % Ав сопсане СЬев 1С 11ев оп СЬе педасйче в16е о1 СЬе 11пе )о1п1пя % СЬе пехС Счо чегсйсея. ТЬАя ргорегсу Ая Стае а1зо 11 роз1С1не ап6 % пе8аСАне аге 1псегсЬвлвей 1п СЬе ргесей1пв Счо вепсепсев. % ?с гв аввише6 СЬаС СЬе ро1увоп 1в с1овей.

?1 поС, с1ове 1С. И в1хе(ху, 1) == 1 ! 1зе9иа1(ху(1, :), ху(еп6, :)) ху(еп6 + 1, :) = ху(1, :); епй Л В 603)) % 818п сопчептйоп: вдп = 1 1ог сопчех негС1сев (1.е, 1птетбог % ап81е > 0 апб <= 180 беЕгеев), вЕп = -1 1ог сопсаче четсйсев. '/ Ехтгеше ротпсв Со Ье ивеб гп СЬе 1о11оч1пЕ 1оор. А 1 1в аррепбеб % со ретХогш СЬе йппет (бес) ргобисС вйСЬ в, вййсЬ 1в 1-Ьу-3 (вее '/.