INSTITUT FÜR INFORMATIK UND PRAKTISCHE MATHEMATIK
|
|
- Jane Wilkerson
- 6 years ago
- Views:
Transcription
1 INSTITUT FÜR INFORMATIK UND PRAKTISCHE MATHEMATIK Structure Multivector for Local Analysis of Images Michael Felsberg and Gerald Sommer Bericht Nr Februar 2000 CHRISTIAN-ALBRECHTS-UNIVERSITÄT KIEL
2 Institut für Informatik und Praktische Mathematik der Christian-Albrechts-Universität zu Kiel Olshausenstr. 40 D 24098Kiel Structure Multivector for Local Analysis of Images Michael Felsberg and Gerald Sommer Bericht Nr Februar mfe@ks.informatik.uni-kiel.de Dieser Bericht ist als persönliche Mitteilung aufzufassen. This work has been supported by German National Merit Foundation and by DFG Graduiertenkolleg No. 357(M. Felsberg) and by DFG Grant So (G. Sommer).
3 Abstract The structure multivector is a new approach for analyzing the local properties of a two-dimensional signal(e.g. images). It combines the classical concepts of the structuretensorandtheanalyticsignalinanewway.thishasbeenmadepossible by using a representation in the algebra of quaternions. The resulting method is linear and of low complexity. The filter-response includes local phase, local amplitude and local orientation of intrinsically one-dimensional neighborhoods in the signal. As for the structure tensor, the structure multivector field can be used to apply special filters to it for detecting two-dimensional features like corners. Experiments and comparisons with other approaches have been made. 1
4 Contents 1 Introduction 3 2 A new approach for the 2D analytic signal Fundamentals The2Dphaseconcept Thesphericalanalyticsignal Properties of the Spherical Analytic Signal Thespatialrepresentation Thestructuremultivector Detectionofkeypoints Experiments and Discussion Implementationandtests Discussionofresults Futurework
5 1 Introduction In image and image sequence processing, different paradigms of interpreting the signals exist. Regardless of they are following a constructive or an appearance based strategy, they all need a capable low-level preprocessing scheme. One of the surely most discussed topics is the detection of key points, or points of higher intrinsic dimension (in the context of appearance based approaches one should better talk of representation than of detection). Key points are needed for solving the correspondence problem occuring in 3D reconstruction, for applying steerable filters(see e.g.[16, 20]), for all kinds of model based representations, etc.. Sincethepreprocessingisonlythefirstlinkinalongchainofoperations,itisuseful to have a linear approach, because otherwise it would be nearly impossible to design the higher-level processing steps. On the other hand, we need a rich representation if we want to treat as much as possible in the preprocessing. Furthermore, the representation of the signal during the different operations should be complete, in order to prevent a loss of information. These constraints enforce us to use the framework of geometric algebra which is also advantageous if we combine image processing with neural computing and robotics(see[18]). The representation of points with special intrinsic dimensionality needs a mathematical treatment of the relevant data, which means that structural information has to be separated from textures and noise. In the one-dimensional case, quadrature filters are a frequently used approach for processing data. But the standard extension to two dimensions is not very efficient because it closely depends on the preference direction of the filter. The structure tensor(see e.g.[13]) is a capable approach for detecting the existence and orientation of local, intrinsic one-dimensional neighborhoods. From the tensor field the orientation vector field can be constructed and by a normalized or differential convolution special symmetries can be detected(see e.g.[14]). The structure tensor itself can be computed with quadrature filters(see[10]) but the tensor itself does not possess the typical properties of a quadrature filter. Especially the linearity and the split of the identity is lost, because the phase is neglected. Inthisreport,weintroduceanewapproachforthe2Danalyticsignalwhichenables ustosubstitutethestructuretensorbyanentitywhichislinear,preservesthesplitofthe identity and has a geometrically meaningful representation: the structure multivector. 3
6 2 Anewapproachforthe2Danalyticsignal Inthefirstsection,werepeatsomefundamentals,inordertomakethenewconcepts easier to understand. In the following sections, we introduce the new approach for the analytic signal in two steps. Firstly, we define the coordinate system, in which we work in the sequel and which also follows directly from geometric constraints. Secondly, the spherical analytic signal is defined using the introduced coordinate system. 2.1 Fundamentals Sinceweworkonimages,whichcanbetreatedassampledintervalsof Ê ¾,weusethe geometricalgebra Ê ¼¾ whichisisomorphictothealgebraofquaternions À.Thewhole complex signal theory naturally embeds in the algebra of quaternions, e.g. complex numbers are considered as a subspace of quaternions here. The basis of the quaternions reads ½ whilethebasisofthecomplexnumbersreads ½. Normally,the basis vector ½ is omitted. Throughout this document, we use the following notations: incontrasttomultivectors,vectors(1dand2d)areboldface Ü Ü ½ Ü ¾, Ù Ù ½ Ù ¾, Ø Ø,etc.and2Dvectorshavethedirections ³ Ü, ³ Ù,etc. thefouriertransformofthe ÒDsignal ܵisdenoted ܵ Æ Ùµ Ê Ò Üµ ÜÔ ¾Ù ܵ Ü therealpart,the -part,the -part,andthe -partofaquaternion Õisobtainedby ÊÕ, ÁÕ,  Õ,and ÃÕ,respectively theautomorphismsofthequaternionsaredenotedby «Õµ Õ, Õµ Õ, and Õµ Õ,wheretheconjugationofthecomplexnumberscanbeexpressed eitherby Þµorby Þµ(Þ ÊÞ ÁÞ) The1Danalyticsignalcanbedefinedasfollows(seee.g.[13]).Thesignalwhichis obtainedfrom ܵbyaphaseshiftof ¾iscalledtheHilberttransform À ܵof ܵ. Since À ܵshallbereal-valued,thespectrummusthaveanoddsymmetry.Therefore, thetransferfunctionhastheform 1 À Ùµ Ò Ùµ (1) 1 Sinceweusevectornotationfor1Dfunctions,wehavetoredefinesomereal-valuedfunctionsaccordingto Ò Ùµ Ò Ùµ,where Ù Ù. 4
7 If we combine a signal and its Hilbert transform corresponding to ܵ ܵ À ܵ (2) wegetacomplex-valuedsignal,whichiscalledtheanalyticsignalof ܵ. Accordingto(1),theFouriertransformoftheanalyticsignal ܵislocatedinthe right half-space of the frequency domain: ܵ Æ Ùµ ¾ Ùµ Ò Ùµ Ùµ ¼µ ¼ if Ù ¼ if Ù ¼ if Ù ¼ Hence, we have recalled all facts of the one-dimensional analytic signal which are relevant for this report. More information can be found in any good textbook about signaltheory. Sinceweusethealgebraofquaternionsasageometricalgebra,wehavetointroduce some terms of geometric algebra(for a complete introduction, see e.g. [12]). If the algebraofquaternionsisidentifiedwiththegeometricalgebra Ê ¼¾,therealpartisalso calledthescalar,theentityconsistingof -and -partsiscalledthevectorandthe -part iscalledthebivector( ). Wehavealreadyusedtwoproducts:theinnerproduct,e.g. Ü Ù,andthegeometric product,e.g. Ù. Thegeometricproductisevaluatedbyapplyingtheassociativelaw and the multiplication table of the algebra. The inner product is well known for vectors, but it can be extended to elements of all grades(geometric algebras are graded algebras) by (3) Ö (4) where isthegradeoperatorand Ö and.theinnerproductwitha scalar is always zero. Theouterproductisdefinedby Ö (5) whichyieldstheidentity ÜÙ ÜÙ Ü Ùforvectors.Ifatleastonefactorisascalar,the outer product is identical to the geometric product respective the field multiplication. Arotationin2Dvectorspacecanbeexpressedbyaspinor,inthecaseofquaternions forexample,wehaveascalar,avector Ü Ü ½ Ü ¾ andabivector.together,they formthequaternion Õ Ü.Ifwewanttorotatethevector(andonlythevector) 5
8 by ³,thiscanbedonebythespinor ÜÔ ³¾µ: ÜÔ ³¾µÕ ÜÔ ³¾µ ÜÔ ³¾µÜ ÜÔ ³¾µ Ó ³¾µÜ ½ Ò ³¾µÜ ¾ µ Ó ³¾µÜ ¾ Ò ³¾µÜ ½ µµ ÜÔ ³¾µ Ó ¾ ³¾µÜ ½ ¾ Ó ³¾µ Ò ³¾µÜ ¾ Ò ¾ ³¾µÜ ½ µ ¾ Ó ³¾µ Ò ³¾µÜ ½ Ò ¾ ³¾µÜ ¾ Ó ¾ ³¾µÜ ¾ µ Ü ½ Ó ³µ Ü ¾ Ò ³µµ Ü ½ Ò ³µ Ü ¾ Ó ³µµ 2.2 The 2D phase concept Wewanttodevelopanew2Danalyticsignalforintrinsically1Dsignals(incontrastto the 2D analytic signal in[4] which is designed for intrinsically 2D signals), which shall include three properties: local amplitude, local phase and local orientation. Compared tothe1danalyticsignalweneedoneadditionalphase. Wecannotchoosethisphase withoutconstraints:ifthesignalisrotatedby,weobtainthesameanalyticsignal,but conjugated. Therefore, we have the following relationship: negation of the local phase isidenticaltoarotationof. One straightforward consequence is that points where the analytic signal has zero phase(i.e.itisreal-valued)haveorientation ³ Ù and ³ Ù atthesametime-therefore, these points have no orientation at all. Note the difference between direction and orientationinthiscontext;thedirectionisavaluein ¼ ¾µandtheorientationisavaluein ¼ µ(see[9]). Asecondconsequenceisthefollowing. Anyvalueofthe2Danalyticsignalcan beunderstoodasa3dvector. Theamplitudefixesthesphereonwhichthevalueis located.thelocalphasecorrespondstorotationsonagreatcircleonthissphere.tobe consistent,arotationofthesignalmustthencorrespondtoarotationonasmallcircle. An obvious parameterization is the following: assuming that we have a complexvalued function over a two-dimensional domain, we map the imaginary unit onto the normalized position vector(the direction vector). The resulting function is quaternionvalued and has the following properties: 6
9 Ifwechooseonefixedorientation ³ Ù intherighthalf-space,theimaginaryunitis changedinto Ó ³ Ù µ Ò ³ Ù µµ Ò Ù ½ µ. For every orientation, the function lies on a different complex plane(note that everyalgebraover ½ Õ ¼, Õ ¼ apureunitquaternion,isisomorphictothealgebra ofcomplexnumbers).thephaseinthisplaneisdenoted ³. The values of the function include information about the orientation of the function. Thecoordinatesystemdefinedinthiswayisdisplayedinfigure1. Itisthesame -part ³Ù -part ³ real part Figure 1: Coordinate system of the 2D phase approach asin[10],butgranlundandknutssonusethe2dphaseinthecontextoforientation adaptive filtering(see also[11]). Consider a given quaternion Õ, how can we extract the two phases, described just above? In fact, the answer of this question defines the introduced coordinate system. 7
10 Definition1Let Ö Þµbethefunctionwhichevaluatesthecomplexphaseof Þintheinterval ¼ ¾µandlet Õ beaquaternionwithzero part.then,wegettheorientationphase ³ Ù by ³ Ù Ö ÁÕ Â Õ µ ¾ µ¾ (6) andthephaseinthecomplexsubspacespannedby ½and Ó ³ Ù µ Ò ³ Ù µisobtainedby ³ Ö ÊÕ ÁÕ Â Õ µ ÜÔ ³ Ù µµ Ö ÊÕ Â Õ ÁÕ µ ÜÔ ³ Ù µµ (7) Note that these definitions of quaternionic phases are totally different from those of the quaternionic Fourier transform(see e.g. [4]). The reason why we use a different phasewillbeexplainedattheendofsection The spherical analytic signal Now,havingaphaseconceptwhichisrichenoughtocodealllocalpropertiesofintrinsically one-dimensional signals, we construct a generalized Hilbert transform and an analyticsignalforthe2dcase,whichmakeuseofthenewembedding. The following definition of the spherical Hilbert transform is motivated by theorem 1, which establishes a correspondence between the Hilbert transform and the spherical Hilbert transform. Definition2Let ÙµbetheFouriertransformof ܵ.ThenthesphericalHilberttransformof ܵisdefinedintheFourierdomainby À Ùµ À Ùµ Ùµ Ù Ùµ (8) Ù and À Ùµ Æ À ܵ. Notethattheoperator À Ùµhastobeappliedfromtheleft.Insection3.1,itwillbe shown that the components of the spherical Hilbert transform are identical to the Riesz transform. 8
11 Example 1: ܵ Ó ¾Ù ¼ ܵ À Ùµ À Ùµ Æ ¼ Ù Ù ¼ µ Æ ¼ Ù Ù ¼ µ ÙÆ ¼ Ù Ù ¼ µ ÙÆ ¼ Ù Ù ¼ µ ¾Ù Ù ¼Æ ¼ Ù Ù ¼ µ Ù ¼ Æ ¼ Ù Ù ¼ µ ¾Ù ¼ Ù Æ ¼ Ù ¼ ¼ Ù Ù ¼ µ Æ ¼ Ù Ù ¼ µ ¾ ¾ sothat À Ùµ Æ ÜÔ ³ ¼ µ Ò ¾Ù ¼ ܵ,where ³ ¼ Ö Ù ¼ µ. Example 2: ܵ Ò ¾Ù ¼ ܵ À Ùµ Ù Æ ¼ Ù ¼ ¼ Ù Ù ¼ µ Æ ¼ Ù Ù ¼ µ ¾ sothat À Ùµ Æ ÜÔ ³ ¼ µ Ó ¾Ù ¼ ܵ. Obviously, the spherical Hilbert transform yields functions which are identical to the 1D Hilbert transforms of the cosine and sine function, respectively, except for an additional exponential function. Up to now, we have only considered two special examples, but what about general signals? Whatkindofsignalscanbetreatedwiththisapproach? Theanswercanbe found easily: the orientation phase must be independent of the frequency coordinate. This sounds impossible, but in fact, the orientation phase is constant, if the spectrum is located on a line through the origin. Signals which have a spectrum of this form are intrinsically one-dimensional(i.e. theyareconstantinonedirection).thisisexactlytheclassoffunctionsforwhichthe structure tensor has been designed. For these functions we have the following theorem: Theorem1Let صbeaone-dimensionalfunctionwiththeHilberttransform À ص. Then, thesphericalhilberttransformofthetwo-dimensionalfunction ¼ ܵ Ü Òµµreads ¼ À ܵ Ò À Ü Òµµ,where Ò Ó µ Ò µisanarbitraryunitvector. 9
12 Proof: TheFouriertransformof ¼ ܵreads Ù ÒµµÆ ¼ ÙÒµif ص Æ µ.corresponding to(8)weobtain Ù Ù Ù ÒµµÆ ¼ Ù Òµ Ò Ò Ù Òµ Ù ÒµµÆ ¼ Ù Òµ andwecantaketheterm ÒintheinverseFouriertransformoutsidetheintegral: and therefore, Ò Ò Ù Òµ Ù ÒµµÆ ¼ Ù Òµ ÜÔ ¾Ù ܵ Ù Ê ¾ Ò Ê ¾ Ò Ò µ µ ÜÔ ¾Ò ܵ Ê Ò À Ò Üµµ Ò Ù Òµ Ù ÒµµÆ ¼ Ù Òµ ÜÔ ¾Ù ܵ Ù ¼ À ܵ Ò À Ü Òµµ ÜÔ µ À Ü Òµµ Nowwesimplyadapt(2)forthe2Dcaseandobtain: Definition3Thesphericalanalyticsignalofa2Dsignal ܵisobtainedby where À ܵisthesphericalHilberttransformof ܵ. ܵ ܵ À ܵ (9) Using this definition, we obtain for our two examples: Ó ¾Ù ¼ ܵ ÜÔ ³ ¼ µ Ò ¾Ù ¼ ܵ Ó ¾Ù ¼ ܵ Ù¼ Ù ¼ Ò ¾Ù ¼ ܵ ÜÔ Ù¼ Ù ¼ ¾Ù ¼ ܵ Ò ¾Ù ¼ ܵ ÜÔ ³ ¼ µ Ó ¾Ù ¼ ܵ Ò ¾Ù ¼ ܵ Ù ¼ and Ù ¼ Ó ¾Ù ¼ ܵ ÜÔ Ù¼ Ù ¼ ¾Ù ¼ Ü ¾µµ So the spherical analytic signal uses the phase concept, which has been defined in section 2.2. According to theorem 1, the spherical analytic signal of an intrinsically onedimensionalsignal ¼ ܵ Ü Òµµreads ¼ ܵ Ü Òµµ Ò À Ü Òµµ (10) 10
13 Of course, the spherical analytic signal can be computed for all functions which are Fourier transformable. However, for signals which do not have an intrinsic dimension ofone 2,thecorrespondencetothe1Danalyticsignalislost. Independently of the intrinsic dimensionality of the signal, the analytic signal can alsobecalculatedinadifferentway. The1Danalyticsignalcanbeobtainedinthe Fourierdomainbythetransferfunction ½ Ò Ùµ. Forthesphericalanalyticsignal wehavethesameresultifwemodifythefouriertransform. Definition4TheinversesphericalFouriertransformof Ùµisdefinedby ܵ ÜÔ ³ Ù ¾µ Ùµ ÜÔ ¾Ù ܵ ÜÔ ³ Ù ¾µ Ù (11) Ê ¾ Before we investigate the inverse of this transform and some of its properties, we establish a relation between the spherical analytic signal and the inverse spherical Fourier transform: Theorem2Thesphericalanalyticsignalof ܵistheinversesphericalFouriertransformof Ùµ. Proof: Ù TheinverseFouriertransformof ½ Ù µ Ùµcanbecalculatedbyanintegrationover Ù or Ù, the integrals are identical. Therefore, we have ܵ ¾ ½ Ù µ Ùµ ÜÔ ¾Ù ܵ Ù Ê Ù ½ ¾ ¾ ½ Ù µ Ùµ ÜÔ ¾Ù ܵ ½ Ù µ Ùµ ÜÔ ¾Ù ܵ Ù Ê Ù Ù Ù Ê Ùµ ÜÔ ¾Ù ܵ Á Ùµ ÜÔ ¾Ù ܵ Ù Ê ¾ Ù Ê ¾ ÜÔ ³ Ù ¾µ Ùµ ÜÔ ¾Ù ܵ ÜÔ ³ Ù ¾µ ٠ܵ 2 Thecaseofintrinsicdimensionzero(i.e.aconstantsignal)isirrelevant,becausetheHilberttransform iszeroinbothcases. 11
14 Above, we mentioned that the spherical analytic signal can be obtained by a similar transferfunctionasinthe1dcaseifthemodifiedfouriertransformisused.uptonow, we used the whole spectrum, but since ÜÔ ³ Ù ¾µ Ùµ ÜÔ ¾Ù ܵ ÜÔ ³ Ù ¾µ ÜÔ ³ Ù µ¾µ Ùµ ÜÔ ¾Ù ܵµ µ ÜÔ ³ Ù µ¾µ ÜÔ ³ Ù µ¾µ Ùµ ÜÔ ¾Ù ܵµ µ ÜÔ ³ Ù µ¾µ ÜÔ ³ Ù µ¾µ Ùµ ÜÔ ¾Ù ܵ ÜÔ ³ Ù µ¾µ i.e., the integrand is symmetric, we can also integrate over the half domain and multiply theintegralbytwo.therefore,wecanuseanytransferfunctionoftheform ½ Ò Ù Òµ withoutchangingtheintegral.themostworkablecasesarethosewhere Ò or Ò. By simply omitting half of the data, the redundancy in the representation is removed. Theinverseofthetransformdefinedby(11)isthetransformwhichtakesaspherical analytic signal and returns the classical complex spectrum. Since the original signal is simplytherealpartoftheanalyticsignal,wehave Ê Üµ Æ Ùµ (12) In order to calculate the energy of the spherical analytic signal, we firstly need the transferfunction,whichchanges Ùµto Ùµ:itisobtainedfrom(8)and(9)andreads Ù ½ Ù ½ Ó ³ Ùµ Ò ³ Ù µ.theenergyofthesphericalanalyticsignalis ½ Ó ³ Ù µ Ò ³ Ù µµ Ùµ ¾ Ù Ê ¾ Ê ¾ ½ Ó ³ Ù µµ ¾ Ùµ ¾ Ò ¾ ³ Ù µ Ùµ ¾ Ù ¾ ½ Ó ³ Ù µµ Ùµ ¾ Ù Ê ¾ andsince Ùµ Ùµ ¾ ¾ ½ Ó ³ Ù µµ Ùµ ¾ ½ Ê Ê ¾ Ùµ ¾ Ù i.e.itistwotimestheenergyoftheoriginalsignal 3. Ó ³ Ù µµ Ùµ ¾ Ù From the group of similarity transformations(i.e. shifts, rotations and dilations) only the rotation really affects the spherical analytic signal; the orientation phase is changed according to the rotation. If we interpret the spherical analytic signal as a vector field 3 ThisisonlyvalidforDCfreesignals.TheenergyoftheDCcomponentisnotdoubledasinthecase of the 1D analytic signal. 12
15 in3d(seealsosection3.2),thegroupof2dsimilaritytransforms 4 evencommuteswith the operator which yields the spherical analytic signal. The reader might ask, why do we use a quaternion-valued Fourier transform which differsfromtheqft(see[4,6,7]andforamoremathematicaldiscussion[2])? The reason is not obvious. Therefore, we will go into some detail. The QFT covers more symmetry concepts than the complex Fourier transform. The classical transform maps a reflection at the origin onto the conjugation operator. The QFTmapsareflectionatoneoftheaxesontooneofthealgebraautomorphisms.This canbeshownbytheisomorphismbetween Å and ¾ : InthecalculationoftheQFT,onlyproductsoftheform,,, and areused, theproducts,, and neverappear.therefore,wecanusethealgebra Å insteadof ÉtocalculatetheQFT.Moreover,wehave Å ¾ (seee.g.[5,8])and consequently, we can calculate the QFT of a real signal by two complex transformations using the formula Õ Ùµ Ùµ ½ ¾ Ù ½ Ù ¾ µ ½ ¾ (13) andthesymmetrywrt.theaxesisobvious. In this report, we want to present an isotropic approach which means that symmetry wrt. theaxesisnotsufficient. Therefore,wehadtodesignthenewtransform. The designofisotropicdiscretefiltersisaquiteoldtopic,seee.g.[3],butstillcurrent(e.g. [15]). 4 Notethatinthecontextofa3Dembedding,thegroupof2Dsimilaritytransformsisthesubgroup ofthe3dtransformsrestrictedtoshiftsin or direction,rotationsaroundtherealaxisandadilationof the and axes. 13
16 3 Properties of the Spherical Analytic Signal In this section we describe some global and local aspects of the spherical analytic signal. 3.1 The spatial representation The definition of the spherical Hilbert transform in the frequency domain can be transformed into a spatial integral transform. The transfer function(8) can be split into two functions: Ù ½ Ù and Ù ¾ Ù.TheonlythingleftistocalculatetheinverseFouriertransformof thesetwofunctions.in[19,17]thetransformpairscanbefound: Ü ½ ¾Ü Æ Ü ¾ ¾Ü Æ Ù ½ Ù and (14) Ù ¾ Ù (15) These Fourier correspondences can also be found by direct calculation. First we try tofindtheinversetransformof ½ Ù.ItcanbeobtainedusingtheHankeltransform,see [1]: ½ ÜÔ ¾Ù Ü Ù Ê ¾ Ù substitute Ü Ö Ó µ Ò µµand Ù Õ Ó ³µ Ò ³µµ ½ ¾ Õ ½ ÜÔ ¾ÖÕ Ó ³ µõ Õ ³ ¼ ¼ whichisthehankeltransformof Õ ½ Ö ½ Ü ½ Now, using the derivative theorem of the Fourier transform we obtain: and since we obtain(14). Ü ½ Ü ½ Æ ¾ Ù ½ Ù Ü ½ Ü ½ Ü The functions ½ ¾Ü and Ü ¾ ¾Ü arethekernelsofthe2driesztransform([19]).from a mathematician s point of view, the Riesz transform is the multidimensional generalization of the Hilbert transform. From(14) and(15) we obtain Ü ½ Ü Ü ½ ¾Ü Æ Ù½ Ù and (16) Ü ¾ ¾Ü Æ Ù¾ Ù (17) 14
17 by the linearity of the Fourier transform Ü ¾Ü Æ Ù Ù (18) Consequently,thesphericalHilberttransform(ortheRiesztransform)ofasignal ܵ isobtainedbyconvolvingthesignalwith Ü ¾Ü : À ܵ Ê ¾ Ø Ü Øµ Ø (19) ¾Ø The spherical analytic signal is directly obtained by the convolution ܵ Ê ¾ Æ ¼ ص Ø ¾Ø Ü Øµ Ø (20) Thegraphinfigure2sumsupallwaystocalculatethesphericalanalyticsignalfrom theprecedingsections.theinversesphericalfouriertransformisdenoted ½. Æ ¼ ܵ Ü ¾Ü µ ½ Ù Ù µ ½ ½ Figure 2: Three ways to calculate the spherical analytic signal 15
18 3.2 The structure multivector Normally, images are intrinsically two-dimensional, so the concepts described in section 2.3 cannot be applied globally. On the other hand, large areas of images are intrinsically one-dimensional, at least on a certain scale. Therefore, a local processing would take advance of the new approach. Compared to local properties of 1D signals, intrinsic one-dimensional signals only have one additional property: the local orientation. Therefore, it is quite reasonable to extend the local description by this parameter. The classical approach of the analytic signal has its local counterpart in the quadrature filters. A pair of quadrature filters(or a complex quadrature filter) is characterized bythefactthattheimpulseresponseisananalyticsignal. Ontheotherhand,both impulse responses can be bounded(i.e. they act on a local neighborhood), so that the problem of the unlimited impulse response of the Hilbert filter is circumvented. Anexampleoftheoutputofa1Dquadraturefiltercanbefoundinfigure3.While Figure 3: Upper left: impulse, upper right: magnitude of filter output, lower left: real part of filter output, lower right: imaginary part of filter output therepresentationinfigure3isverycommon,wewillintroducenowadifferentrep- 16
19 resentation. One-dimensional signals can be interpreted as surfaces in 2D space. If we assigntherealaxistothesignalvaluesandtheimaginaryaxistotheabscissa,weobtain a representation in the complex plane. The analytic signal can be embedded in the sameplane itcorrespondstoavectorfieldwhichisonlynon-zeroontheimaginary axis(see figure 4) Figure 4: Representation of the 1D analytic signal as a vector field The vector field representation of the analytic signal looks like particles which roll fromahill(impulse).therotationisdirectedawayfromthepeakandthespeed(amplitude of vectors) decreases with increasing distance. Based on the spherical analytic signal we now introduce the spherical quadrature filters. They are defined according to the 1D case. Definition 5 A hypercomplex filter whose impulse response is a spherical analytic signal is called a spherical quadrature filter. In practical applications, the convolution(20) cannot be applied directly either, because of the unlimited impulse response of the expression(20). Of course, one could create optimized filters for signals which are supposed to be band-limited, but it is better workable to use radial bandpass filters to band-limit the signal actively. The spatial representation of a radial bandpass can be used to design the spherical quadrature filter in the spatial domain. The spatial representation looks similar to a Bessel function(see figure 5). In fact, it is a linear combination of Bessel functions with different scales. It is remarkable that the spherical quadrature filters have isotropic energy and exactlychoosethefrequencybandstheyaredesignedfor.infigure6itcanbeseenthat theenergyisisotropicandthatitismaximalfortheradius,whichcorrespondsto afrequencyof ½.Theusedbandpasshasacenterfrequencyof ½. ¾ ½ 17
20 Figure 5: Bessel function Figure 6: Siemens-star convolved with a spherical quadrature filter(magnitude) 18
21 Infigure7theoutputofasphericalquadraturefilterappliedtoanimpulse-lineis displayed. The lower right image shows the combined imaginary parts which means Figure 7: Upper left: impulse-line, upper right: amplitude of filter output, lower left: real part of filter output, lower right: combined imaginary parts of filter output thatitissimplytheimaginarypartofthe correct imaginaryplane(thatonewhichis obtained from the orientation phase). Sameasfor1Dsignals,2Dsignalscanbeembeddedasasurfacein3Dspace.The signalvaluesareassignedtotherealaxisandthespatialcoordinatestothe -andthe -axis. Thesphericalanalyticsignalcanberepresentedinthesameembedding. It correspondstoavectorfield,whichisonlynon-zerointheplanespannedby and (see figure 8). Asinthe1Dcase,thevectorfieldlookslikeparticleswhichrollfromaridge(impulse-line). The rotation is directed orthogonal away from the ridge and the speed decreases with increasing distance. The result of filtering a signal with a spherical quadrature filter is a quaternionvaluedfield. Thoughthebivector(or )componentofthefieldisalwayszero,we denote this field as a multivector field or the structure multivector of the signal. As already the name induces, the structure multivector is closely related to the structure tensor. The structure tensor as defined in[10] includes the following information (in fact, only the orientation vector is considered): amplitude measurement for the existence of local structure 19
22 Figure 8: Representation of the spherical analytic signal as a vector field orientation orientation of the local structure Jähne[13] extracts an additional information: the coherence. The coherence is the relationship between the oriented gradients and all gradients, so it is a measurement for thedegreeoforientationinastructureanditiscloselyrelatedtothevarianceofthe orientation. Thevarianceisasecondorderproperty.Itincludesaproductoftheargumentsand therefore, it is not linear. Consequently, the coherence cannot be measured by a linear approach like the structure multivector. Two structures with different orientations simply yield the vector sum of both multivector fields. The structure multivector consists of three independent components(local phase, local orientation and local amplitude) and it codes three properties. Consequently, there is no additional information possible. The structure tensor possesses three degrees of freedom(it is a symmetric tensor) and codes only two properties. Therefore, one can additionally extract a third information, the coherence. 20
23 3.3 Detection of key points The current status of our approach of detecting key points is not totally satisfying, becauseitisstillnonlinear.thekeyideaistocalculatetheouterenergyofthemultivector field. Outer energy means being contrary to the inner energy, which is calculated by the inner product of the multivectors. Obviously, this cannot be done pointwise, because theouterproductofamultivectorwithitselfisalwayszero 5. Therefore, we consider a local neighborhood and in this local neighborhood we combine those areas which are assumed to have similar orientations. This is achieved as follows: the first vector for the outer product is obtained by convolving the signal with acosinemaskandthesecondoneisobtainedbyconvolvingwithasinemask: ܵ À ܵ Ó ³ Ü µö ܵµ (21) Ë Üµ À ܵ Ò ³ Ü µö ܵµ (22) where Ö Üµisafunctionwhichiszerofor Ü ¼andforlarge Ü(e.g. ashifted Gaussian function). Ifwehaveanintrinsically1Dstructure,bothsignalshavethesamephases. Ifwe have a rectangular corner, they are orthogonal. Consequently, the outer product(cross product 6 ) iszerointhefirstcaseandismaximizedinthesecondcase.  ܵ ܵ Ë Üµ (23) 5 Noteinthiscontext,theinnerproductstandsforthesquaredmagnitudeandtheouterproductisthe cross product between the 3D component vectors obtained from the multivectors. 6 Intheformula,weusethenotationofthecrossproductinordertodistinguishbetweenthealgebras weworkon(seealsofootnote5). 21
24 4 Experiments and Discussion This section shows some experiments and an implementation in MatLab. The results arecomparedtothosewhichcanbefoundin[10]forexample. 4.1 Implementation and tests Throughout this section we used the implementations which can be found in the appendix. The function which computes the structure multivector uses a multi-scale approach, becausetheshiftofthegaussianbandpassiscoupledwiththevarianceasforthegabor wavelets. Both functions are first implementations, which should only underline the practicability of the approach. For testing the implementations, we chose some synthetic examples with letters as gray level or textured images. Figure 9: Structure multivector applied to an image without texture. Upper left: original,upperright:amplitude,lowerleft: ³ -phase,lowerright: ³ Ù -phase Infigure9itcanbeseenthatthestructuremultivectorrespondsonlyattheedges. Therefore,theamplitudeisameasureforthepresenceofstructure. The ³ -phaseis 22
25 Figure 10: Structure multivector applied to an image with one texture. Upper left: original,upperright:amplitude,lowerleft: ³ -phase,lowerright: ³ Ù -phase Figure 11: Structure multivector applied to an image with two superposed textures. Upperleft:original,upperright:amplitude,lowerleft: ³ -phase,lowerright: ³ Ù -phase 23
26 Figure 12: Structure multivector applied to an image with two superposed textures and texturedbackground.upperleft:original,upperright:amplitude,lowerleft: ³ -phase, lowerright: ³ Ù -phase linear,whichcanonlybeguessedinthisrepresentation. Butitcanbeseenthatthe ³ -phaseismonotonicmoduloamaximalinterval(whichisinfact ¾).The ³ Ù -phase represents the orientation of the edge. Note that the highest and the lowest gray level (standingfor andzero,respectively)representthesameorientation. Infigure10itcanbeseenthatthestructuremultivectorrespondsalsoinsidethe object. The amplitude is nearly constant, which corresponds to a texture with constant energy.ofcourse,thispropertyislostifthewrongscaleisconsidered.the ³ -phase islinear,seenotesabove. The ³ Ù -phaserepresentstheorientationofthetexture. A constant gray level corresponds to a constant estimated orientation. The small spikes in the figure are produced by the extraction of the local orientation angle. The underlying quaternion-valued field does not show these artifacts. Infigure11itcanbeseenthatthestructuremultivectorrespondsonlywithrespect to the dominant texture. The magnitude of the response is modulated with that componentoftheweakertexturethatisnormaltothedominanttexture.thiseffectiseven moreobviousinfigure12. The ³ -phaseisalwaysdirectedparalleltothedominant texture.the ³ Ù -phaserepresentstheorientationofthedominanttextureineachcase. The described 2D detection algorithm works fine for images without texture. The amplitude of the detection filter response can be seen in figure 13. Obviously, the response increases with increasing sharpness of the corner. 24
27 Figure 13: Result of the 2D key-point detection applied to the multivector field from figure Discussion of results In contrast to the structure tensor which becomes zero for ideal 2D points(e.g. rectangular corner) the structure multivector responds at a high energy level. Therefore, it is easier to design a detector for those points, because the energy is higher. Furthermore, a concept like the outer energy does not produce artificial double lines beneath the original edges. These double lines makes it necessary to use a normalized convolution. Our approach does not need such tricks to obtain clear cut results. 4.3 Futurework Inourfuturework,wewilltrytolinearizethe2Dkeypointdetection.Thiswillbedone in an 8D algebra. Furthermore, we will try to separate texture information from edge information in a second filtering step, using the same algebra. Up to now, textures which are not parallel to the edge produce artificial 2D key points. The long-term objective will be the design of filters which covers the entire hierarchy of intrinsic dimensionality, consisting of(1d and 2D) textures, edges, corners and planes. Another extension of the presented theory will be the consideration of higher dimensional cases, especially the image sequence analysis. The approach of the spherical analytic signal can be extended to higher dimensions easily. The theoretical background will be further investigated with regard to the results from Clifford analysis, especially the monogenic function and other multidimensional generalizations of the analytic function. 25
28 Appendix MatLab function for computing the structure multivector ÙÒØÓÒÅ ØÖÑÛ Ê Ñ µ ± ÓÑÔÙØ Ø ØÖÙØÙÖ ÑÙÐØÚØÓÖ Ó ÖÐ ÒÐ Ê ÛØ ± ÔÖÐ ÕÙÖØÙÖ ÐØÖ ÛØ Ù Ò ÒÔ ± Ѿ ÚÖÒ ÖØÓ Ó ÛÚ ÒÖÒ ¼º Ò Þ Ê ½µ»¾ ¾¾ Ñ Ù Ú Ñ Ö ¹ ¹½ ¹ ¹½ µ ÕÖØ Ùº¾ Úº¾µ ÜÔ ¹¾ Ô» ѹµµº¾µ Ñ ÕÖØ ¾ Ôµ ½ ½µ¼ ÆØ Ø Ø¾ Ø Ø µµµ ÆÆ ½¹ ¾ ¾ ½¹ ¾ ¾µ ½ ½µ½ º» ÙÙº ÚÚº ÍØ Ø Ø¾ Ø Ø Ùµµµ ÎØ Ø Ø¾ Ø Ø Úµµµ ÍÍ ½¹ ¾ ¾ ½¹ ¾ ¾µ ÎÎ ½¹ ¾ ¾ ½¹ ¾ ¾µ Å ¾µÓÒÚ¾ Ê ¹Ñ ͵ ³ ѳµ Å µóòú¾ Ê ¹Ñ ε ³ ѳµ Å ½µÓÒÚ¾ Ê ÖРƵ ³ ѳµ 26
29 MatLab function for detecting intrinsically 2D structures ÙÒØÓÒ Øؾ Å Ñ Ûµ ± ØØÓÒ Ó ¾ ݹÔÓÒØ ÛØ ÖÐ ÐÓ¹ÒÓÖÑÐ ÛØÒ Ù Ú Ñ Ö ¹ ¹ µ Ñ ÕÖØ Ùº¾ Úº¾µ Ñ µ½ ÔÜÔ ¹Û ÐÓ Ñ»Ñµº¾»ÐÓ ¾µµº»Ñ Ô µ¼ ÙÙº Ô ÚÚº Ô ÓÒÚ¾ Å ½µ Ù ³ ѳµ ÓÒÚ¾ Å ½µ Ú ³ ѳµ ÓÒÚ¾ Å ¾µ Ù ³ ѳµ ÓÒÚ¾ Å ¾µ Ú ³ ѳµ ÓÒÚ¾ Å µ Ù ³ ѳµ ÓÒÚ¾ Å µ Ú ³ ѳµ º ¹ º º ¹º º ¹ º º¾ º¾ º¾ 27
30 References [1] BRACEWELL, R. N. Two-Dimensional Imaging. Prentice Hall signal processing series. Prentice Hall, Englewood Cliffs, [2] BRACKX, F., DELANGHE, R., AND SOMMEN, F. Clifford Analysis. Pitman, Boston, [3] BRADY, J. M., AND HORN, B. M. P. Rotationally symmetric operators for surface interpolation. Computer Vision, Graphics, and Image Processing 22, 1(April 1983), [4] BÜLOW, T. Global and Local Hypercomplex Spectral Signal Representations for Image Processing and Analysis. PhD thesis, Christian-Albrechts-University of Kiel, [5] DAVENPORT, C. M. A Commutative Hypercomplex Algebra with Associated Function Theory. In Clifford Algebras with Numeric and Symbolic Computations, R. Ablamowicz, Ed. Birkhäuser Bosten, 1996, pp [6] ELL, T. A. Hypercomplex Spectral Transformations. PhD thesis, University of Minnesota, [7] ERNST, R. R., AUE, W. P., BACHMANN, P., KARHAN, J., KUMAR, A., AND MÜLLER, L. Two-dimensional NMR spectroscopy. In Proc. 4th Ampère Int. Summer School, Pula, Yugoslavia(1976). [8] FELSBERG, M., BÜLOW, T., AND SOMMER, G. Commutative Hypercomplex Fourier Transforms of Multidimensional Signals. In Geometric Computing with Clifford Algebra, G. Sommer, Ed. Springer, Heidelberg, to appear. [9] GRANLUND, G. H. Hierarchical computer vision. In Proc. of EUSIPCO-90, Fifth European Signal Processing Conference, Barcelona(1990), L. Torres, E. Masgrau, and M. A. Lagunas, Eds., pp [10] GRANLUND, G. H., AND KNUTSSON, H. Signal Processing for Computer Vision. Kluwer Academic Publishers, Dordrecht, [11] HAGLUND, L. Adaptive Multidimensional Filtering. PhD thesis, Linköping University, [12] HESTENES, D., AND SOBCZYK, G. Clifford algebra to geometric calculus, A Unified Language for Mathematics and Physics. Reidel, Dordrecht, [13] JÄHNE, B. Digitale Bildverarbeitung. Springer, Berlin, [14] KARLHOLM, J. Local Signal Models for Image Sequence Analysis. PhD thesis, Linköping University,
31 [15] MERRON, J., AND BRADY, M. Isotropic gradient estimation. In IEEE Computer Vision and Pattern Recognition(1996), pp [16] MICHAELIS, M., AND SOMMER, G. A Lie group approach to steerable filters. Pattern Recognition Letters 16(1995), [17] NABIGHIAN, M. N. Toward a three-dimensional automatic interpretation of potential field data via generalized Hilbert transforms: Fundamental relations. Geophysics 49, 6(June 1984), [18] SOMMER, G. The global algebraic frame of the perception-action cycle. In Handbook of Computer Vision and Applications(1999), B. Jähne, H. Haußecker, and P. Geissler, Eds., vol. 3, Academic Press, San Diego, pp [19] STEIN, E., AND WEISS, G. Introduction to Fourier Analysis on Euclidean Spaces. Princeton University Press, New Jersey, [20] YU, W., DANIILIDIS, K., BEAUCHEMIN, S., AND SOMMER, G. Detection and characterization of multiple motion points. In 18th IEEE conference on Computer Vision and Pattern Recognition, Fort Collins, Colorado(1999). 29
Linköping University Post Print. The monogenic signal
Linköping University Post Print The monogenic signal Michael Felsberg Gerald Sommer N.B.: When citing this work, cite the original article. 2009 IEEE. Personal use of this material is permitted. However,
More informationThe GET Operator. LiTH-ISY-R-2633
The GET Operator Michael Felsberg Linköping University, Computer Vision Laboratory, SE-58183 Linköping, Sweden, mfe@isy.liu.se, http://www.isy.liu.se/cvl/ October 4, 004 LiTH-ISY-R-633 Abstract In this
More informationA Probabilistic Definition of Intrinsic Dimensionality for Images
A Probabilistic Definition of Intrinsic Dimensionality for Images Michael Felsberg 1 and Norbert Krüger 2 1 Computer Vision Laboratory, Dept. of Electrical Engineering, Linköping University SE-58183 Linköping,
More informationA Novel Approach to the 2D Analytic Signal? Thomas Bulow and Gerald Sommer. Christian{Albrechts{Universitat zu Kiel
A Novel Approach to the 2D Analytic Signal? Thomas Bulow and Gerald Sommer Christian{Albrechts{Universitat zu Kiel Institute of Computer Science, Cognitive Systems Preuerstrae 1{9, 24105 Kiel Tel:+49 431
More informationHypercomplex Signals A Novel Extension of the Analytic Signal to the Multidimensional Case
2844 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 49, NO. 11, NOVEMBER 2001 Hypercomplex Signals A Novel Extension of the Analytic Signal to the Multidimensional Case Thomas Bülow and Gerald Sommer Abstract
More informationDemystification of the Geometric Fourier Transforms
Demystification of the Geometric Fourier Transforms Roxana Buack, Gerik Scheuermann and Eckhard Hitzer Universität Leipzig, Institut für Informatik, Abteilung für Bild- und Signalverarbeitung, Augustuplatz
More informationAn Uncertainty Principle for Quaternion Fourier Transform
An Uncertainty Principle for Quaternion Fourier Transform Mawardi Bahri a, Eckhard S. M. Hitzer a Akihisa Hayashi a Ryuichi Ashino b, a Department of Applied Physics, University of Fukui, Fukui 9-857,
More informationA UNIFIED THEORY FOR STEERABLE AND QUADRATURE FILTERS
A UNIFIED THEORY FOR STEERABLE AND QUADRATURE FILTERS Kai Krajsek J. W. Goethe University, Visual Sensorics and Information Processing Lab Robert Mayer Str.2-4, D-60054 Frankfurt am Main, Germany krajsek@vsi.cs.uni-frankfurt.de
More informationOptical Flow Estimation from Monogenic Phase
Optical Flow Estimation from Monogenic Phase Michael Felsberg Linköping University, Computer Vision Laboratory, SE-58183 Linköping, Sweden, mfe@isy.liu.se, http://www.isy.liu.se/cvl/ Abstract. The optical
More informationGeometry and Kinematics with Uncertain Data
Geometry and Kinematics with Uncertain Data Christian Perwass, Christian Gebken, and Gerald Sommer Institut für Informatik, CAU Kiel Christian-Albrechts-Platz 4, 24118 Kiel, Germany {chp,chg,gs}@ks.informatik.uni-kiel.de
More informationApplications of Discrete Mathematics to the Analysis of Algorithms
Applications of Discrete Mathematics to the Analysis of Algorithms Conrado Martínez Univ. Politècnica de Catalunya, Spain May 2007 Goal Given some algorithm taking inputs from some set Á, we would like
More informationTwo-Dimensional Clifford Windowed Fourier Transform
Two-Dimensional Clifford Windowed Fourier Transform Mawardi Bahri, Eckhard M. S. Hitzer and Sriwulan Adji Abstract Recently several generalizations to higher dimension of the classical Fourier transform
More informationCONVEX OPTIMIZATION OVER POSITIVE POLYNOMIALS AND FILTER DESIGN. Y. Genin, Y. Hachez, Yu. Nesterov, P. Van Dooren
CONVEX OPTIMIZATION OVER POSITIVE POLYNOMIALS AND FILTER DESIGN Y. Genin, Y. Hachez, Yu. Nesterov, P. Van Dooren CESAME, Université catholique de Louvain Bâtiment Euler, Avenue G. Lemaître 4-6 B-1348 Louvain-la-Neuve,
More informationAGEOMETRICALGEBRA APPROACH TO SOME PROBLEMS OF ROBOT VISION
AGEOMETRICALGEBRA APPROACH TO SOME PROBLEMS OF ROBOT VISION Gerald Sommer Institut für Informatik und Praktische Mathematik Christian-Albrechts-Universität zu Kiel, Kiel, Germany gs@ks.informatik.uni-kiel.de
More information2 Hallén s integral equation for the thin wire dipole antenna
Ú Ð Ð ÓÒÐ Ò Ø ØØÔ»» Ѻ Ö Ùº º Ö ÁÒغ º ÁÒ Ù ØÖ Ð Å Ø Ñ Ø ÎÓк ÆÓº ¾ ¾¼½½µ ½ ¹½ ¾ ÆÙÑ Ö Ð Ñ Ø Ó ÓÖ Ò ÐÝ Ó Ö Ø ÓÒ ÖÓÑ Ø Ò Û Ö ÔÓÐ ÒØ ÒÒ Ëº À Ø ÑÞ ¹Î ÖÑ ÞÝ Ö Åº Æ Ö¹ÅÓ Êº Ë Þ ¹Ë Ò µ Ô ÖØÑ ÒØ Ó Ð ØÖ Ð Ò Ò
More informationOn Urquhart s C Logic
On Urquhart s C Logic Agata Ciabattoni Dipartimento di Informatica Via Comelico, 39 20135 Milano, Italy ciabatto@dsiunimiit Abstract In this paper we investigate the basic many-valued logics introduced
More informationThe Monogenic Wavelet Transform Sofia C. Olhede and Georgios Metikas
3426 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 57, NO. 9, SEPTEMBER 2009 The Monogenic Wavelet Transform Sofia C. Olhede and Georgios Metikas Abstract This paper extends the 1-D analytic wavelet transform
More informationGET: The Connection Between Monogenic Scale-Space and Gaussian Derivatives
GET: The Connection Between Monogenic Scale-Space and Gaussian Derivatives Michael Felsberg 1 and Ullrich Köthe 2 1 Linköping University, Computer Vision Laboratory, SE-58183 Linköping, Sweden, mfe@isy.liu.se,
More informationSelf-Testing Polynomial Functions Efficiently and over Rational Domains
Chapter 1 Self-Testing Polynomial Functions Efficiently and over Rational Domains Ronitt Rubinfeld Madhu Sudan Ý Abstract In this paper we give the first self-testers and checkers for polynomials over
More informationLund Institute of Technology Centre for Mathematical Sciences Mathematical Statistics
Lund Institute of Technology Centre for Mathematical Sciences Mathematical Statistics STATISTICAL METHODS FOR SAFETY ANALYSIS FMS065 ÓÑÔÙØ Ö Ü Ö Ì ÓÓØ ØÖ Ô Ð ÓÖ Ø Ñ Ò Ý Ò Ò ÐÝ In this exercise we will
More informationEnergy Tensors: Quadratic, Phase Invariant Image Operators
Energy Tensors: Quadratic, Phase Invariant Image Operators Michael Felsberg and Erik Jonsson Linköping University, Computer Vision Laboratory, SE-58183 Linköping, Sweden Abstract. In this paper we briefly
More informationDense Optical Flow Estimation from the Monogenic Curvature Tensor
Dense Optical Flow Estimation from the Monogenic Curvature Tensor Di Zang 1, Lennart Wietzke 1, Christian Schmaltz 2, and Gerald Sommer 1 1 Department of Computer Science, Christian-Albrechts-University
More informationRandomized Simultaneous Messages: Solution of a Problem of Yao in Communication Complexity
Randomized Simultaneous Messages: Solution of a Problem of Yao in Communication Complexity László Babai Peter G. Kimmel Department of Computer Science The University of Chicago 1100 East 58th Street Chicago,
More informationLecture 16: Modern Classification (I) - Separating Hyperplanes
Lecture 16: Modern Classification (I) - Separating Hyperplanes Outline 1 2 Separating Hyperplane Binary SVM for Separable Case Bayes Rule for Binary Problems Consider the simplest case: two classes are
More informationThe structure tensor in projective spaces
The structure tensor in projective spaces Klas Nordberg Computer Vision Laboratory Department of Electrical Engineering Linköping University Sweden Abstract The structure tensor has been used mainly for
More informationObservations on the Stability Properties of Cooperative Systems
1 Observations on the Stability Properties of Cooperative Systems Oliver Mason and Mark Verwoerd Abstract We extend two fundamental properties of positive linear time-invariant (LTI) systems to homogeneous
More informationA Language for Task Orchestration and its Semantic Properties
DEPARTMENT OF COMPUTER SCIENCES A Language for Task Orchestration and its Semantic Properties David Kitchin, William Cook and Jayadev Misra Department of Computer Science University of Texas at Austin
More informationRiesz-Transforms vs. Derivatives: On the Relationship Between the Boundary Tensor and the Energy Tensor
Riesz-Transforms vs. Derivatives: On the Relationship Between the Boundary Tensor and the Energy Tensor Ullrich Köthe 1 and Michael Felsberg 1 Cognitive Systems Group, University of Hamburg koethe@informatik.uni-hamburg.de
More informationChapter 1 Numerical evaluation of Versors with Clifford Algebra
This is page 1 Printer: Opaque this Chapter 1 Numerical evaluation of Versors with Clifford Algebra Christian B. U. Perwass, Gerald Sommer 1 ABSTRACT This paper has two main parts. In the first part we
More informationGET: The Connection Between Monogenic Scale-Space and Gaussian Derivatives
GET: The Connection Between Monogenic Scale-Space and Gaussian Derivatives Michael Felsberg 1 and Ullrich Köthe 2 1 Linköping University, Computer Vision Laboratory, SE-58183 Linköping, Sweden, mfe@isy.liu.se,
More informationAnalyzing Real Vector Fields with Clifford Convolution and Clifford Fourier Transform
Analyzing Real Vector Fields with Clifford Convolution and Clifford Fourier Transform Wieland Reich and Gerik Scheuermann Abstract Post-processing in computational fluid dynamics as well as processing
More informationLocal textures parameters of images in space domain obtained with different analytical approaches
European Journal of Mathematics and Computer Science Vol. 4 No. 1, 017 Local textures parameters of images in space domain obtained with different analytical approaches Robert Goutte University of Lyon,
More informationImage enhancement. Why image enhancement? Why image enhancement? Why image enhancement? Example of artifacts caused by image encoding
13 Why image enhancement? Image enhancement Example of artifacts caused by image encoding Computer Vision, Lecture 14 Michael Felsberg Computer Vision Laboratory Department of Electrical Engineering 12
More informationMacLane s coherence theorem expressed as a word problem
MacLane s coherence theorem expressed as a word problem Paul-André Melliès Preuves, Programmes, Systèmes CNRS UMR-7126, Université Paris 7 ÑÐÐ ÔÔ ºÙ ÙºÖ DRAFT In this draft, we reduce the coherence theorem
More informationFast Fourier Transform Solvers and Preconditioners for Quadratic Spline Collocation
Fast Fourier Transform Solvers and Preconditioners for Quadratic Spline Collocation Christina C. Christara and Kit Sun Ng Department of Computer Science University of Toronto Toronto, Ontario M5S 3G4,
More informationand multiplication as q + p = (q 0 + iq 1 + jq 2 + kq 3 ) + (p 0 + ip 1 + j p 2 + kp 3 )
Ben Witten* and Jeff Shragge, Stanford University SUMMAY Hypercomlex numbers are primarily used for pattern recognition, offer many useful applications to geophysics. Image disparity estimation is a hypercomplex,
More informationA Generalization of Principal Component Analysis to the Exponential Family
A Generalization of Principal Component Analysis to the Exponential Family Michael Collins Sanjoy Dasgupta Robert E. Schapire AT&T Labs Research 8 Park Avenue, Florham Park, NJ 7932 mcollins, dasgupta,
More informationMotivation towards Clifford Analysis: The classical Cauchy s Integral Theorem from a Clifford perspective
Motivation towards Clifford Analysis: The classical Cauchy s Integral Theorem from a Clifford perspective Lashi Bandara November 26, 29 Abstract Clifford Algebras generalise complex variables algebraically
More informationMargin Maximizing Loss Functions
Margin Maximizing Loss Functions Saharon Rosset, Ji Zhu and Trevor Hastie Department of Statistics Stanford University Stanford, CA, 94305 saharon, jzhu, hastie@stat.stanford.edu Abstract Margin maximizing
More informationArchiv der Mathematik Holomorphic approximation of L2-functions on the unit sphere in R^3
Manuscript Number: Archiv der Mathematik Holomorphic approximation of L-functions on the unit sphere in R^ --Manuscript Draft-- Full Title: Article Type: Corresponding Author: Holomorphic approximation
More informationSUBOPTIMALITY OF THE KARHUNEN-LOÈVE TRANSFORM FOR FIXED-RATE TRANSFORM CODING. Kenneth Zeger
SUBOPTIMALITY OF THE KARHUNEN-LOÈVE TRANSFORM FOR FIXED-RATE TRANSFORM CODING Kenneth Zeger University of California, San Diego, Department of ECE La Jolla, CA 92093-0407 USA ABSTRACT An open problem in
More informationThe orthogonal planes split of quaternions and its relation to quaternion geometry of rotations 1
The orthogonal planes split of quaternions and its relation to quaternion geometry of rotations 1 Eckhard Hitzer Osawa 3-10-, Mitaka 181-8585, International Christian University, Japan E-mail: hitzer@icu.ac.jp
More informationFinal exam: Automatic Control II (Reglerteknik II, 1TT495)
Uppsala University Department of Information Technology Systems and Control Professor Torsten Söderström Final exam: Automatic Control II (Reglerteknik II, TT495) Date: October 22, 2 Responsible examiner:
More informationf k j M k j Among them we recognize the gradient operator on polynomial-valued functions which is part of the complete gradient (see [9]) u j xj.
Advances in Applied Clifford Algebras 4, No. 1 (1994) 65 FUNCTIONS OF TWO VECTOR VARIABLES F. Sommen* and N. Van Acker Department of Mathematical Analysis, University of Gent, Galglaan 2 B-9000 Gent, Belgium
More informationOn the Equivariance of the Orientation and the Tensor Field Representation Klas Nordberg Hans Knutsson Gosta Granlund Computer Vision Laboratory, Depa
On the Invariance of the Orientation and the Tensor Field Representation Klas Nordberg Hans Knutsson Gosta Granlund LiTH-ISY-R-530 993-09-08 On the Equivariance of the Orientation and the Tensor Field
More informationContinuous-time Fourier Methods
ELEC 321-001 SIGNALS and SYSTEMS Continuous-time Fourier Methods Chapter 6 1 Representing a Signal The convolution method for finding the response of a system to an excitation takes advantage of the linearity
More informationCambridge University Press The Mathematics of Signal Processing Steven B. Damelin and Willard Miller Excerpt More information
Introduction Consider a linear system y = Φx where Φ can be taken as an m n matrix acting on Euclidean space or more generally, a linear operator on a Hilbert space. We call the vector x a signal or input,
More informationNarrowing confidence interval width of PAC learning risk function by algorithmic inference
Narrowing confidence interval width of PAC learning risk function by algorithmic inference Bruno Apolloni, Dario Malchiodi Dip. di Scienze dell Informazione, Università degli Studi di Milano Via Comelico
More informationFrequency domain representation and singular value decomposition
EOLSS Contribution 643134 Frequency domain representation and singular value decomposition AC Antoulas Department of Electrical and Computer Engineering Rice University Houston, Texas 77251-1892, USA e-mail:
More informationThe Fourier series are applicable to periodic signals. They were discovered by the
3.1 Fourier Series The Fourier series are applicable to periodic signals. They were discovered by the famous French mathematician Joseph Fourier in 1807. By using the Fourier series, a periodic signal
More informationA quick introduction to Clifford algebras
A quick introduction to Clifford algebras Paul Leopardi paul.leopardi@unsw.edu.au Presented at School of Mathematics, University of New South Wales, 2003-06-05. Real representation matrix R 2 8 16 8 16
More informationarxiv:math/ v1 [math.dg] 6 Jul 1998
arxiv:math/9807024v [math.dg] 6 Jul 998 The Fundamental Theorem of Geometric Calculus via a Generalized Riemann Integral Alan Macdonald Department of Mathematics Luther College, Decorah, IA 520, U.S.A.
More informationarxiv: v1 [math.dg] 19 Mar 2012
BEREZIN-TOEPLITZ QUANTIZATION AND ITS KERNEL EXPANSION XIAONAN MA AND GEORGE MARINESCU ABSTRACT. We survey recent results [33, 34, 35, 36] about the asymptotic expansion of Toeplitz operators and their
More informationAutomatic Control III (Reglerteknik III) fall Nonlinear systems, Part 3
Automatic Control III (Reglerteknik III) fall 20 4. Nonlinear systems, Part 3 (Chapter 4) Hans Norlander Systems and Control Department of Information Technology Uppsala University OSCILLATIONS AND DESCRIBING
More informationLecture Notes 5: Multiresolution Analysis
Optimization-based data analysis Fall 2017 Lecture Notes 5: Multiresolution Analysis 1 Frames A frame is a generalization of an orthonormal basis. The inner products between the vectors in a frame and
More informationA Tutorial on Data Reduction. Principal Component Analysis Theoretical Discussion. By Shireen Elhabian and Aly Farag
A Tutorial on Data Reduction Principal Component Analysis Theoretical Discussion By Shireen Elhabian and Aly Farag University of Louisville, CVIP Lab November 2008 PCA PCA is A backbone of modern data
More informationAn Introduction to Optimal Control Applied to Disease Models
An Introduction to Optimal Control Applied to Disease Models Suzanne Lenhart University of Tennessee, Knoxville Departments of Mathematics Lecture1 p.1/37 Example Number of cancer cells at time (exponential
More informationEstimators for Orientation and Anisotropy in Digitized Images
Estimators for Orientation and Anisotropy in Digitized Images Lucas J. van Vliet and Piet W. Verbeek Pattern Recognition Group of the Faculty of Applied Physics Delft University of Technolo Lorentzweg,
More informationScale Space Smoothing, Image Feature Extraction and Bessel Filters
Scale Space Smoothing, Image Feature Extraction and Bessel Filters Sasan Mahmoodi and Steve Gunn School of Electronics and Computer Science, Building 1, Southampton University, Southampton, SO17 1BJ, UK
More informationNew Foundations in Mathematics: The Geometric Concept of Number
New Foundations in Mathematics: The Geometric Concept of Number by Garret Sobczyk Universidad de Las Americas-P Cholula, Mexico November 2012 What is Geometric Algebra? Geometric algebra is the completion
More informationVARIATIONAL CALCULUS IN SPACE OF MEASURES AND OPTIMAL DESIGN
Chapter 1 VARIATIONAL CALCULUS IN SPACE OF MEASURES AND OPTIMAL DESIGN Ilya Molchanov Department of Statistics University of Glasgow ilya@stats.gla.ac.uk www.stats.gla.ac.uk/ ilya Sergei Zuyev Department
More informationHyperbolic Conformal Geometry with Clifford Algebra 1)
MM Research Preprints, 72 83 No. 18, Dec. 1999. Beijing Hyperbolic Conformal Geometry with Clifford Algebra 1) Hongbo Li Abstract. In this paper we study hyperbolic conformal geometry following a Clifford
More informationAn EM algorithm for Batch Markovian Arrival Processes and its comparison to a simpler estimation procedure
An EM algorithm for Batch Markovian Arrival Processes and its comparison to a simpler estimation procedure Lothar Breuer (breuer@info04.uni-trier.de) University of Trier, Germany Abstract. Although the
More informationarxiv: v1 [math.ca] 6 Feb 2015
The Fourier-Like and Hartley-Like Wavelet Analysis Based on Hilbert Transforms L. R. Soares H. M. de Oliveira R. J. Cintra Abstract arxiv:150.0049v1 [math.ca] 6 Feb 015 In continuous-time wavelet analysis,
More informationIndex. p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96
p, lip, 78 8 function, 107 v, 7-8 w, 7-8 i,7-8 sine, 43 Bo,94-96 B 1,94-96 M,94-96 B oro!' 94-96 BIro!' 94-96 I/r, 79 2D linear system, 56 2D FFT, 119 2D Fourier transform, 1, 12, 18,91 2D sinc, 107, 112
More informationSignal Modeling Techniques in Speech Recognition. Hassan A. Kingravi
Signal Modeling Techniques in Speech Recognition Hassan A. Kingravi Outline Introduction Spectral Shaping Spectral Analysis Parameter Transforms Statistical Modeling Discussion Conclusions 1: Introduction
More information4.3 Laplace Transform in Linear System Analysis
4.3 Laplace Transform in Linear System Analysis The main goal in analysis of any dynamic system is to find its response to a given input. The system response in general has two components: zero-state response
More informationThe New Graphic Description of the Haar Wavelet Transform
he New Graphic Description of the Haar Wavelet ransform Piotr Porwik and Agnieszka Lisowska Institute of Informatics, Silesian University, ul.b dzi ska 39, 4-00 Sosnowiec, Poland porwik@us.edu.pl Institute
More informationRandom Matrices and Transfer Matrix Techniques for the Integer Quantum Hall Effect
Random Matrices and Transfer Matrix Techniques for the Integer Quantum Hall Effect Dissertation zur Erlangung des akademischen Grades eines Doktors der Naturwissenschaften an der Mathematisch-Naturwissenschaftlichen
More informationSpherical Diffusion. ScholarlyCommons. University of Pennsylvania. Thomas Bülow University of Pennsylvania
University of Pennsylvania ScholarlyCommons Technical Reports (CIS) Department of Computer & Information Science 12-30-2001 Spherical Diffusion Thomas Bülow University of Pennsylvania Follow this and additional
More informationAppendix A Wirtinger Calculus
Precoding and Signal Shaping for Digital Transmission. Robert F. H. Fischer Copyright 0 2002 John Wiley & Sons, Inc. ISBN: 0-471-22410-3 Appendix A Wirtinger Calculus T he optimization of system parameters
More informationDiscrete Fourier Transform
Last lecture I introduced the idea that any function defined on x 0,..., N 1 could be written a sum of sines and cosines. There are two different reasons why this is useful. The first is a general one,
More informationCovariance Tracking Algorithm on Bilateral Filtering under Lie Group Structure Yinghong Xie 1,2,a Chengdong Wu 1,b
Applied Mechanics and Materials Online: 014-0-06 ISSN: 166-748, Vols. 519-50, pp 684-688 doi:10.408/www.scientific.net/amm.519-50.684 014 Trans Tech Publications, Switzerland Covariance Tracking Algorithm
More informationTOPICS IN HARMONIC ANALYSIS WITH APPLICATIONS TO RADAR AND SONAR. Willard Miller
TOPICS IN HARMONIC ANALYSIS WITH APPLICATIONS TO RADAR AND SONAR Willard Miller October 23 2002 These notes are an introduction to basic concepts and tools in group representation theory both commutative
More informationLOCALIZED analyses in one dimension(al) (1-D) have
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 55, NO. 3, MARCH 2007 921 Multiple Multidimensional Morse Wavelets Georgios Metikas Sofia C. Olhede Abstract This paper defines a set of operators that localize
More informationConnecting spatial and frequency domains for the quaternion Fourier transform
Connecting spatial and frequency domains for the quaternion Fourier transform Ghent University (joint work with N. De Schepper, T. Ell, K. Rubrecht and S. Sangwine) MOIMA, Hannover, June, 2016 Direct formulas
More informationUsing Skew Gabor Filter in Source Signal Separation and Local Spectral Multi-Orientation Analysis
Using Skew Gabor Filter in Source Signal Separation and Local Spectral Multi-Orientation Analysis Weichuan Yu Gerald Sommer Kostas Daniilidis 3 Dept. of Diagnostic Radiology Institut für Informatik 3 GRASP
More informationLearnability and the Doubling Dimension
Learnability and the Doubling Dimension Yi Li Genome Institute of Singapore liy3@gis.a-star.edu.sg Philip M. Long Google plong@google.com Abstract Given a set of classifiers and a probability distribution
More informationAdvanced Signals and Systems
TECHNISCHE FAKULTÄT DER CHRISTIAN-ALBRECHTS-UNIVERSITÄT ZU KIEL DIGITALE SIGNALVERARBEITUNG UND SYSTEMTHEORIE DSS Advanced Signals and Systems Exam SS 08 Examiner: Prof. Dr.-Ing. Gerhard Schmidt Date:
More informationGeometric algebra: a computational framework for geometrical applications (part I: algebra)
Geometric algebra: a computational framework for geometrical applications (part I: algebra) Leo Dorst and Stephen Mann Abstract Geometric algebra is a consistent computational framework in which to define
More information14 Singular Value Decomposition
14 Singular Value Decomposition For any high-dimensional data analysis, one s first thought should often be: can I use an SVD? The singular value decomposition is an invaluable analysis tool for dealing
More informationLinear Operators and Fourier Transform
Linear Operators and Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013
More informationFeature Extraction and Image Processing
Feature Extraction and Image Processing Second edition Mark S. Nixon Alberto S. Aguado :*авш JBK IIP AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO
More informationCorrelation at Low Temperature: I. Exponential Decay
Correlation at Low Temperature: I. Exponential Decay Volker Bach FB Mathematik; Johannes Gutenberg-Universität; D-55099 Mainz; Germany; email: vbach@mathematik.uni-mainz.de Jacob Schach Møller Ý Département
More informationClifford Fourier Transform for Color Image Processing
Clifford Fourier Transform for Color Image Processing Thomas Batard, Michel Berthier and Christophe Saint-Jean Abstract The aim of this paper is to define a Clifford Fourier transform that is suitable
More informationPeriodic monopoles and difference modules
Periodic monopoles and difference modules Takuro Mochizuki RIMS, Kyoto University 2018 February Introduction In complex geometry it is interesting to obtain a correspondence between objects in differential
More informationMultiscale Image Transforms
Multiscale Image Transforms Goal: Develop filter-based representations to decompose images into component parts, to extract features/structures of interest, and to attenuate noise. Motivation: extract
More informationImproving the Berlekamp algorithm for binomials x n a
Improving the Berlekamp algorithm for binomials x n a Ryuichi Harasawa Yutaka Sueyoshi Aichi Kudo Nagasaki University July 19, 2012 1 2 3 Idea Example Procedure after applying the proposed mehod 4 Theoretical
More informationArbeitstagung: Gruppen und Topologische Gruppen Vienna July 6 July 7, Abstracts
Arbeitstagung: Gruppen und Topologische Gruppen Vienna July 6 July 7, 202 Abstracts ÁÒÚ Ö Ð Ñ Ø Ó Ø¹Ú ÐÙ ÙÒØ ÓÒ ÁÞØÓ Ò ÞØÓ º Ò ÙÒ ¹Ñ º ÙÐØÝ Ó Æ ØÙÖ Ð Ë Ò Ò Å Ø Ñ Ø ÍÒ Ú Ö ØÝ Ó Å Ö ÓÖ ÃÓÖÓ ½ ¼ Å Ö ÓÖ ¾¼¼¼
More informationOrientation Estimation in Ambiguous Neighbourhoods Mats T. Andersson & Hans Knutsson Computer Vision Laboratory Linköping University 581 83 Linköping, Sweden Abstract This paper describes a new algorithm
More informationLie Algebra of Unit Tangent Bundle in Minkowski 3-Space
INTERNATIONAL ELECTRONIC JOURNAL OF GEOMETRY VOLUME 12 NO. 1 PAGE 1 (2019) Lie Algebra of Unit Tangent Bundle in Minkowski 3-Space Murat Bekar (Communicated by Levent Kula ) ABSTRACT In this paper, a one-to-one
More informationLars Schmidt-Thieme, Information Systems and Machine Learning Lab (ISMLL), Institute BW/WI & Institute for Computer Science, University of Hildesheim
Course on Information Systems 2, summer term 2010 0/29 Information Systems 2 Information Systems 2 5. Business Process Modelling I: Models Lars Schmidt-Thieme Information Systems and Machine Learning Lab
More informationF O R SOCI AL WORK RESE ARCH
7 TH EUROPE AN CONFERENCE F O R SOCI AL WORK RESE ARCH C h a l l e n g e s i n s o c i a l w o r k r e s e a r c h c o n f l i c t s, b a r r i e r s a n d p o s s i b i l i t i e s i n r e l a t i o n
More informationA New Scheme for Anisotropic Diffusion: Performance Evaluation of Rotation Invariance, Dissipativity and Efficiency no author given no address given A
A New Scheme for Anisotropic Diffusion: Performance Evaluation of Rotation Invariance, Dissipativity and Efficiency Hanno Scharr 1 and Joachim Weickert 1 Interdisciplinary Center for Scientific Computing,
More informationThe Spectral Basis of a Linear Operator
The Spectral Basis of a Linear Operator Garret Sobczyk Universidad de Las Americas-P Cholula, Mexico http://www.garretstar.com Thursday Jan. 10, 2013, 2PM AMS/MAA Joint Math Meeting San Diego Convention
More informationOperator Exponentials for the Clifford Fourier Transform on Multivector Fields
Operator Exponentials for the Clifford Fourier Transform on Multivector Fields David Eelbode and Echard Hitzer Mathematics Subject Classification 000). Primary 4A38; Secondary 15A66,34L40. Keywords. Fourier
More informationSCALE and directionality are key considerations for visual
636 IEEE TRANSACTIONS ON IMAGE PROCESSING, VOL. 19, NO. 3, MARCH 2010 Wavelet Steerability and the Higher-Order Riesz Transform Michael Unser, Fellow, IEEE, and Dimitri Van De Ville, Member, IEEE Abstract
More informationA General Geometric Fourier Transform
A General Geometric Fourier Transform Roxana Bujack, Gerik Scheuermann and Eckhard Hitzer Abstract. The increasing demand for Fourier transforms on geometric algebras has resulted in a large variety. Here
More informationRECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK
RECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK TRNKA PAVEL AND HAVLENA VLADIMÍR Dept of Control Engineering, Czech Technical University, Technická 2, 166 27 Praha, Czech Republic mail:
More informationMerging Mathematical Technologies by Applying the Reverse bra-ket Method By J.A.J. van Leunen
Merging Mathematical Technologies by Applying the Reverse bra-ket Method By J.A.J. van Leunen http://www.e-physics.eu Abstract Last modified: 28 September 2016 Hilbert spaces can store discrete members
More information