16acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn/*M/////////////////////////////////////////////////////////////////////////////////////// 26acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 36acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// IMPORTANT: READ BEFORE DOWNLOADING, COPYING, INSTALLING OR USING. 46acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 56acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// By downloading, copying, installing or using the software you agree to this license. 66acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// If you do not agree to this license, do not download, install, 76acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// copy or use the software. 86acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 96acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 106acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Intel License Agreement 116acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// For Open Source Computer Vision Library 126acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 136acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Copyright (C) 2000, Intel Corporation, all rights reserved. 146acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Third party copyrights are property of their respective owners. 156acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 166acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Redistribution and use in source and binary forms, with or without modification, 176acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// are permitted provided that the following conditions are met: 186acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 196acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// * Redistribution's of source code must retain the above copyright notice, 206acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// this list of conditions and the following disclaimer. 216acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 226acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// * Redistribution's in binary form must reproduce the above copyright notice, 236acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// this list of conditions and the following disclaimer in the documentation 246acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// and/or other materials provided with the distribution. 256acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 266acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// * The name of Intel Corporation may not be used to endorse or promote products 276acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// derived from this software without specific prior written permission. 286acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 296acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// This software is provided by the copyright holders and contributors "as is" and 306acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// any express or implied warranties, including, but not limited to, the implied 316acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// warranties of merchantability and fitness for a particular purpose are disclaimed. 326acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// In no event shall the Intel Corporation or contributors be liable for any direct, 336acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// indirect, incidental, special, exemplary, or consequential damages 346acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// (including, but not limited to, procurement of substitute goods or services; 356acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// loss of use, data, or profits; or business interruption) however caused 366acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// and on any theory of liability, whether in contract, strict liability, 376acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// or tort (including negligence or otherwise) arising in any way out of 386acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// the use of this software, even if advised of the possibility of such damage. 396acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 406acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn//M*/ 416acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn#include "_cv.h" 426acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 436acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn/* The function calculates center of gravity and central second order moments */ 446acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Rennstatic void 456acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennicvCompleteMomentState( CvMoments* moments ) 466acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn{ 476acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double cx = 0, cy = 0; 486acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double mu20, mu11, mu02; 496acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 506acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn assert( moments != 0 ); 516acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->inv_sqrt_m00 = 0; 526acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 536acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( fabs(moments->m00) > DBL_EPSILON ) 546acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 556acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double inv_m00 = 1. / moments->m00; 566acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn cx = moments->m10 * inv_m00; 576acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn cy = moments->m01 * inv_m00; 586acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->inv_sqrt_m00 = sqrt( fabs(inv_m00) ); 596acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 606acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 616acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* mu20 = m20 - m10*cx */ 626acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mu20 = moments->m20 - moments->m10 * cx; 636acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* mu11 = m11 - m10*cy */ 646acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mu11 = moments->m11 - moments->m10 * cy; 656acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* mu02 = m02 - m01*cy */ 666acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mu02 = moments->m02 - moments->m01 * cy; 676acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 686acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->mu20 = mu20; 696acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->mu11 = mu11; 706acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->mu02 = mu02; 716acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 726acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* mu30 = m30 - cx*(3*mu20 + cx*m10) */ 736acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->mu30 = moments->m30 - cx * (3 * mu20 + cx * moments->m10); 746acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mu11 += mu11; 756acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* mu21 = m21 - cx*(2*mu11 + cx*m01) - cy*mu20 */ 766acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->mu21 = moments->m21 - cx * (mu11 + cx * moments->m01) - cy * mu20; 776acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* mu12 = m12 - cy*(2*mu11 + cy*m10) - cx*mu02 */ 786acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->mu12 = moments->m12 - cy * (mu11 + cy * moments->m10) - cx * mu02; 796acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* mu03 = m03 - cy*(3*mu02 + cy*m01) */ 806acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->mu03 = moments->m03 - cy * (3 * mu02 + cy * moments->m01); 816acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn} 826acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 836acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 846acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Rennstatic void 856acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennicvContourMoments( CvSeq* contour, CvMoments* moments ) 866acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn{ 876acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int is_float = CV_SEQ_ELTYPE(contour) == CV_32FC2; 886acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 896acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( contour->total ) 906acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 916acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CvSeqReader reader; 926acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double a00, a10, a01, a20, a11, a02, a30, a21, a12, a03; 936acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double xi, yi, xi2, yi2, xi_1, yi_1, xi_12, yi_12, dxy, xii_1, yii_1; 946acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int lpt = contour->total; 956acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 966acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a00 = a10 = a01 = a20 = a11 = a02 = a30 = a21 = a12 = a03 = 0; 976acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 986acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn cvStartReadSeq( contour, &reader, 0 ); 996acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1006acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !is_float ) 1016acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 1026acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xi_1 = ((CvPoint*)(reader.ptr))->x; 1036acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yi_1 = ((CvPoint*)(reader.ptr))->y; 1046acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 1056acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn else 1066acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 1076acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xi_1 = ((CvPoint2D32f*)(reader.ptr))->x; 1086acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yi_1 = ((CvPoint2D32f*)(reader.ptr))->y; 1096acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 1106acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_NEXT_SEQ_ELEM( contour->elem_size, reader ); 1116acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1126acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xi_12 = xi_1 * xi_1; 1136acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yi_12 = yi_1 * yi_1; 1146acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1156acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn while( lpt-- > 0 ) 1166acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 1176acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !is_float ) 1186acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 1196acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xi = ((CvPoint*)(reader.ptr))->x; 1206acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yi = ((CvPoint*)(reader.ptr))->y; 1216acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 1226acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn else 1236acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 1246acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xi = ((CvPoint2D32f*)(reader.ptr))->x; 1256acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yi = ((CvPoint2D32f*)(reader.ptr))->y; 1266acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 1276acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_NEXT_SEQ_ELEM( contour->elem_size, reader ); 1286acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1296acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xi2 = xi * xi; 1306acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yi2 = yi * yi; 1316acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn dxy = xi_1 * yi - xi * yi_1; 1326acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xii_1 = xi_1 + xi; 1336acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yii_1 = yi_1 + yi; 1346acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1356acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a00 += dxy; 1366acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a10 += dxy * xii_1; 1376acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a01 += dxy * yii_1; 1386acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a20 += dxy * (xi_1 * xii_1 + xi2); 1396acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a11 += dxy * (xi_1 * (yii_1 + yi_1) + xi * (yii_1 + yi)); 1406acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a02 += dxy * (yi_1 * yii_1 + yi2); 1416acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a30 += dxy * xii_1 * (xi_12 + xi2); 1426acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a03 += dxy * yii_1 * (yi_12 + yi2); 1436acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a21 += 1446acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn dxy * (xi_12 * (3 * yi_1 + yi) + 2 * xi * xi_1 * yii_1 + 1456acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xi2 * (yi_1 + 3 * yi)); 1466acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a12 += 1476acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn dxy * (yi_12 * (3 * xi_1 + xi) + 2 * yi * yi_1 * xii_1 + 1486acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yi2 * (xi_1 + 3 * xi)); 1496acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1506acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xi_1 = xi; 1516acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yi_1 = yi; 1526acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn xi_12 = xi2; 1536acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn yi_12 = yi2; 1546acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 1556acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1566acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double db1_2, db1_6, db1_12, db1_24, db1_20, db1_60; 1576acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1586acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( fabs(a00) > FLT_EPSILON ) 1596acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 1606acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( a00 > 0 ) 1616acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 1626acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_2 = 0.5; 1636acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_6 = 0.16666666666666666666666666666667; 1646acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_12 = 0.083333333333333333333333333333333; 1656acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_24 = 0.041666666666666666666666666666667; 1666acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_20 = 0.05; 1676acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_60 = 0.016666666666666666666666666666667; 1686acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 1696acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn else 1706acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 1716acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_2 = -0.5; 1726acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_6 = -0.16666666666666666666666666666667; 1736acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_12 = -0.083333333333333333333333333333333; 1746acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_24 = -0.041666666666666666666666666666667; 1756acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_20 = -0.05; 1766acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn db1_60 = -0.016666666666666666666666666666667; 1776acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 1786acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1796acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* spatial moments */ 1806acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m00 = a00 * db1_2; 1816acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m10 = a10 * db1_6; 1826acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m01 = a01 * db1_6; 1836acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m20 = a20 * db1_12; 1846acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m11 = a11 * db1_24; 1856acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m02 = a02 * db1_12; 1866acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m30 = a30 * db1_20; 1876acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m21 = a21 * db1_60; 1886acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m12 = a12 * db1_60; 1896acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m03 = a03 * db1_20; 1906acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1916acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvCompleteMomentState( moments ); 1926acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 1936acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 1946acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn} 1956acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1966acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 1976acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn/* summarizes moment values for all tiles */ 1986acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Rennstatic void 1996acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennicvAccumulateMoments( double *tiles, CvSize size, CvSize tile_size, CvMoments * moments ) 2006acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn{ 2016acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int x, y; 2026acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2036acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn for( y = 0; y < size.height; y += tile_size.height ) 2046acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 2056acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn for( x = 0; x < size.width; x += tile_size.width, tiles += 10 ) 2066acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 2076acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double dx = x, dy = y; 2086acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double dxm = dx * tiles[0], dym = dy * tiles[0]; 2096acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2106acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m00 ( = m00' ) */ 2116acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m00 += tiles[0]; 2126acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2136acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m10 ( = m10' + dx*m00' ) */ 2146acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m10 += tiles[1] + dxm; 2156acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2166acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m01 ( = m01' + dy*m00' ) */ 2176acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m01 += tiles[2] + dym; 2186acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2196acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m20 ( = m20' + 2*dx*m10' + dx*dx*m00' ) */ 2206acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m20 += tiles[3] + dx * (tiles[1] * 2 + dxm); 2216acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2226acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m11 ( = m11' + dx*m01' + dy*m10' + dx*dy*m00' ) */ 2236acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m11 += tiles[4] + dx * (tiles[2] + dym) + dy * tiles[1]; 2246acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2256acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m02 ( = m02' + 2*dy*m01' + dy*dy*m00' ) */ 2266acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m02 += tiles[5] + dy * (tiles[2] * 2 + dym); 2276acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2286acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m30 ( = m30' + 3*dx*m20' + 3*dx*dx*m10' + dx*dx*dx*m00' ) */ 2296acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m30 += tiles[6] + dx * (3. * tiles[3] + dx * (3. * tiles[1] + dxm)); 2306acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2316acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m21 (= m21' + dx*(2*m11' + 2*dy*m10' + dx*m01' + dx*dy*m00') + dy*m20') */ 2326acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m21 += tiles[7] + dx * (2 * (tiles[4] + dy * tiles[1]) + 2336acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn dx * (tiles[2] + dym)) + dy * tiles[3]; 2346acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2356acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m12 (= m12' + dy*(2*m11' + 2*dx*m01' + dy*m10' + dx*dy*m00') + dx*m02') */ 2366acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m12 += tiles[8] + dy * (2 * (tiles[4] + dx * tiles[2]) + 2376acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn dy * (tiles[1] + dxm)) + dx * tiles[5]; 2386acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2396acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* + m03 ( = m03' + 3*dy*m02' + 3*dy*dy*m01' + dy*dy*dy*m00' ) */ 2406acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments->m03 += tiles[9] + dy * (3. * tiles[5] + dy * (3. * tiles[2] + dym)); 2416acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 2426acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 2436acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2446acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvCompleteMomentState( moments ); 2456acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn} 2466acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2476acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2486acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn/****************************************************************************************\ 2496acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn* Spatial Moments * 2506acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn\****************************************************************************************/ 2516acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 2526acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn#define ICV_DEF_CALC_MOMENTS_IN_TILE( __op__, name, flavor, srctype, temptype, momtype ) \ 2536acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Rennstatic CvStatus CV_STDCALL icv##name##_##flavor##_CnCR \ 2546acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn( const srctype* img, int step, CvSize size, int cn, int coi, double *moments ) \ 2556acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn{ \ 2566acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int x, y, sx_init = (size.width & -4) * (size.width & -4), sy = 0; \ 2576acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn momtype mom[10]; \ 2586acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 2596acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn assert( img && size.width && (size.width | size.height) >= 0 ); \ 2606acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn memset( mom, 0, 10 * sizeof( mom[0] )); \ 2616acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 2626acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( coi ) \ 2636acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn img += coi - 1; \ 2646acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn step /= sizeof(img[0]); \ 2656acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 2666acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn for( y = 0; y < size.height; sy += 2 * y + 1, y++, img += step ) \ 2676acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { \ 2686acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn temptype x0 = 0; \ 2696acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn temptype x1 = 0; \ 2706acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn temptype x2 = 0; \ 2716acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn momtype x3 = 0; \ 2726acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int sx = sx_init; \ 2736acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn const srctype* ptr = img; \ 2746acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 2756acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn for( x = 0; x < size.width - 3; x += 4, ptr += cn*4 ) \ 2766acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { \ 2776acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn temptype p0 = __op__(ptr[0]), p1 = __op__(ptr[cn]), \ 2786acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn p2 = __op__(ptr[2*cn]), p3 = __op__(ptr[3*cn]); \ 2796acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn temptype t = p1; \ 2806acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn temptype a, b, c; \ 2816acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 2826acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn p0 += p1 + p2 + p3; /* p0 + p1 + p2 + p3 */ \ 2836acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn p1 += 2 * p2 + 3 * p3; /* p1 + p2*2 + p3*3 */ \ 2846acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn p2 = p1 + 2 * p2 + 6 * p3; /* p1 + p2*4 + p3*9 */ \ 2856acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn p3 = 2 * p2 - t + 9 * p3; /* p1 + p2*8 + p3*27 */ \ 2866acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 2876acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a = x * p0 + p1; /* x*p0 + (x+1)*p1 + (x+2)*p2 + (x+3)*p3 */ \ 2886acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn b = x * p1 + p2; /* (x+1)*p1 + 2*(x+2)*p2 + 3*(x+3)*p3 */ \ 2896acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn c = x * p2 + p3; /* (x+1)*p1 + 4*(x+2)*p2 + 9*(x+3)*p3 */ \ 2906acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 2916acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn x0 += p0; \ 2926acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn x1 += a; \ 2936acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn a = a * x + b; /*(x^2)*p0+((x+1)^2)*p1+((x+2)^2)*p2+((x+3)^2)*p3 */ \ 2946acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn x2 += a; \ 2956acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn x3 += ((momtype)(a + b)) * x + c; /*x3 += (x^3)*p0+((x+1)^3)*p1 + */ \ 2966acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* ((x+2)^3)*p2+((x+3)^3)*p3 */ \ 2976acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } \ 2986acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 2996acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn /* process the rest */ \ 3006acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn for( ; x < size.width; sx += 2 * x + 1, x++, ptr += cn ) \ 3016acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { \ 3026acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn temptype p = __op__(ptr[0]); \ 3036acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn temptype xp = x * p; \ 3046acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 3056acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn x0 += p; \ 3066acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn x1 += xp; \ 3076acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn x2 += sx * p; \ 3086acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn x3 += ((momtype)sx) * xp; \ 3096acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } \ 3106acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 3116acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { \ 3126acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn temptype py = y * x0; \ 3136acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 3146acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[9] += ((momtype)py) * sy; /* m03 */ \ 3156acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[8] += ((momtype)x1) * sy; /* m12 */ \ 3166acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[7] += ((momtype)x2) * y; /* m21 */ \ 3176acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[6] += x3; /* m30 */ \ 3186acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[5] += x0 * sy; /* m02 */ \ 3196acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[4] += x1 * y; /* m11 */ \ 3206acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[3] += x2; /* m20 */ \ 3216acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[2] += py; /* m01 */ \ 3226acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[1] += x1; /* m10 */ \ 3236acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mom[0] += x0; /* m00 */ \ 3246acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } \ 3256acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } \ 3266acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 3276acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn for( x = 0; x < 10; x++ ) \ 3286acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moments[x] = (double)mom[x]; \ 3296acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn \ 3306acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn return CV_OK; \ 3316acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn} 3326acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3336acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3346acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennICV_DEF_CALC_MOMENTS_IN_TILE( CV_NOP, MomentsInTile, 8u, uchar, int, int ) 3356acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennICV_DEF_CALC_MOMENTS_IN_TILE( CV_NOP, MomentsInTile, 16u, ushort, int, int64 ) 3366acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennICV_DEF_CALC_MOMENTS_IN_TILE( CV_NOP, MomentsInTile, 16s, short, int, int64 ) 3376acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennICV_DEF_CALC_MOMENTS_IN_TILE( CV_NOP, MomentsInTile, 32f, float, double, double ) 3386acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennICV_DEF_CALC_MOMENTS_IN_TILE( CV_NOP, MomentsInTile, 64f, double, double, double ) 3396acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3406acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennICV_DEF_CALC_MOMENTS_IN_TILE( CV_NONZERO, MomentsInTileBin, 8u, uchar, int, int ) 3416acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennICV_DEF_CALC_MOMENTS_IN_TILE( CV_NONZERO, MomentsInTileBin, 16s, ushort, int, int ) 3426acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennICV_DEF_CALC_MOMENTS_IN_TILE( CV_NONZERO_FLT, MomentsInTileBin, 32f, int, int, int ) 3436acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennICV_DEF_CALC_MOMENTS_IN_TILE( CV_NONZERO_FLT, MomentsInTileBin, 64f, int64, double, double ) 3446acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3456acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn#define icvMomentsInTile_8s_CnCR 0 3466acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn#define icvMomentsInTile_32s_CnCR 0 3476acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn#define icvMomentsInTileBin_8s_CnCR icvMomentsInTileBin_8u_CnCR 3486acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn#define icvMomentsInTileBin_16u_CnCR icvMomentsInTileBin_16s_CnCR 3496acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn#define icvMomentsInTileBin_32s_CnCR 0 3506acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3516acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennCV_DEF_INIT_FUNC_TAB_2D( MomentsInTile, CnCR ) 3526acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennCV_DEF_INIT_FUNC_TAB_2D( MomentsInTileBin, CnCR ) 3536acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3546acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn////////////////////////////////// IPP moment functions ////////////////////////////////// 3556acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3566acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennicvMoments_8u_C1R_t icvMoments_8u_C1R_p = 0; 3576acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennicvMoments_32f_C1R_t icvMoments_32f_C1R_p = 0; 3586acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennicvMomentInitAlloc_64f_t icvMomentInitAlloc_64f_p = 0; 3596acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennicvMomentFree_64f_t icvMomentFree_64f_p = 0; 3606acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennicvGetSpatialMoment_64f_t icvGetSpatialMoment_64f_p = 0; 3616acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3626acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renntypedef CvStatus (CV_STDCALL * CvMomentIPPFunc) 3636acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn ( const void* img, int step, CvSize size, void* momentstate ); 3646acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3656acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennCV_IMPL void 3666acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RenncvMoments( const void* array, CvMoments* moments, int binary ) 3676acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn{ 3686acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn static CvFuncTable mom_tab; 3696acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn static CvFuncTable mombin_tab; 3706acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn static int inittab = 0; 3716acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double* tiles = 0; 3726acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn void* ippmomentstate = 0; 3736acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3746acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_FUNCNAME("cvMoments"); 3756acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3766acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __BEGIN__; 3776acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3786acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int type = 0, depth, cn, pix_size; 3796acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int coi = 0; 3806acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int x, y, k, tile_num = 1; 3816acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CvSize size, tile_size = { 32, 32 }; 3826acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CvMat stub, *mat = (CvMat*)array; 3836acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CvFunc2DnC_1A1P func = 0; 3846acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CvMomentIPPFunc ipp_func = 0; 3856acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CvContour contour_header; 3866acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CvSeq* contour = 0; 3876acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CvSeqBlock block; 3886acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3896acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( CV_IS_SEQ( array )) 3906acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 3916acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn contour = (CvSeq*)array; 3926acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !CV_IS_SEQ_POLYGON( contour )) 3936acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_ERROR( CV_StsBadArg, "The passed sequence is not a valid contour" ); 3946acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 3956acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 3966acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !inittab ) 3976acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 3986acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvInitMomentsInTileCnCRTable( &mom_tab ); 3996acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvInitMomentsInTileBinCnCRTable( &mombin_tab ); 4006acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn inittab = 1; 4016acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 4026acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4036acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !moments ) 4046acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_ERROR( CV_StsNullPtr, "" ); 4056acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4066acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn memset( moments, 0, sizeof(*moments)); 4076acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4086acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !contour ) 4096acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 4106acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_CALL( mat = cvGetMat( mat, &stub, &coi )); 4116acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn type = CV_MAT_TYPE( mat->type ); 4126acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4136acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( type == CV_32SC2 || type == CV_32FC2 ) 4146acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 4156acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_CALL( contour = cvPointSeqFromMat( 4166acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_SEQ_KIND_CURVE | CV_SEQ_FLAG_CLOSED, 4176acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mat, &contour_header, &block )); 4186acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 4196acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 4206acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4216acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( contour ) 4226acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 4236acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvContourMoments( contour, moments ); 4246acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn EXIT; 4256acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 4266acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4276acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn type = CV_MAT_TYPE( mat->type ); 4286acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn depth = CV_MAT_DEPTH( type ); 4296acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn cn = CV_MAT_CN( type ); 4306acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn pix_size = CV_ELEM_SIZE(type); 4316acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn size = cvGetMatSize( mat ); 4326acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4336acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( cn > 1 && coi == 0 ) 4346acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_ERROR( CV_StsBadArg, "Invalid image type" ); 4356acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4366acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( size.width <= 0 || size.height <= 0 ) 4376acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 4386acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn EXIT; 4396acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 4406acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4416acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( type == CV_8UC1 ) 4426acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn ipp_func = (CvMomentIPPFunc)icvMoments_8u_C1R_p; 4436acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn else if( type == CV_32FC1 ) 4446acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn ipp_func = (CvMomentIPPFunc)icvMoments_32f_C1R_p; 4456acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4466acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( ipp_func && !binary ) 4476acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 4486acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int matstep = mat->step ? mat->step : CV_STUB_STEP; 4496acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn IPPI_CALL( icvMomentInitAlloc_64f_p( &ippmomentstate, cvAlgHintAccurate )); 4506acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn IPPI_CALL( ipp_func( mat->data.ptr, matstep, size, ippmomentstate )); 4516acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 0, 0, 0, cvPoint(0,0), &moments->m00 ); 4526acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 1, 0, 0, cvPoint(0,0), &moments->m10 ); 4536acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 0, 1, 0, cvPoint(0,0), &moments->m01 ); 4546acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 2, 0, 0, cvPoint(0,0), &moments->m20 ); 4556acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 1, 1, 0, cvPoint(0,0), &moments->m11 ); 4566acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 0, 2, 0, cvPoint(0,0), &moments->m02 ); 4576acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 3, 0, 0, cvPoint(0,0), &moments->m30 ); 4586acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 2, 1, 0, cvPoint(0,0), &moments->m21 ); 4596acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 1, 2, 0, cvPoint(0,0), &moments->m12 ); 4606acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvGetSpatialMoment_64f_p( ippmomentstate, 0, 3, 0, cvPoint(0,0), &moments->m03 ); 4616acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvCompleteMomentState( moments ); 4626acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn EXIT; 4636acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 4646acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4656acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn func = (CvFunc2DnC_1A1P)(!binary ? mom_tab.fn_2d[depth] : mombin_tab.fn_2d[depth]); 4666acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4676acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !func ) 4686acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_ERROR( CV_StsBadArg, cvUnsupportedFormat ); 4696acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4706acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( depth >= CV_32S && !binary ) 4716acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn tile_size = size; 4726acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn else 4736acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn tile_num = ((size.width + tile_size.width - 1)/tile_size.width)* 4746acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn ((size.height + tile_size.height - 1)/tile_size.height); 4756acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4766acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_CALL( tiles = (double*)cvAlloc( tile_num*10*sizeof(double))); 4776acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4786acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn for( y = 0, k = 0; y < size.height; y += tile_size.height ) 4796acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 4806acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CvSize cur_tile_size = tile_size; 4816acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( y + cur_tile_size.height > size.height ) 4826acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn cur_tile_size.height = size.height - y; 4836acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4846acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn for( x = 0; x < size.width; x += tile_size.width, k++ ) 4856acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 4866acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( x + cur_tile_size.width > size.width ) 4876acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn cur_tile_size.width = size.width - x; 4886acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4896acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn assert( k < tile_num ); 4906acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4916acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn IPPI_CALL( func( mat->data.ptr + y*mat->step + x*pix_size, 4926acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mat->step, cur_tile_size, cn, coi, tiles + k*10 )); 4936acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 4946acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 4956acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4966acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvAccumulateMoments( tiles, size, tile_size, moments ); 4976acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 4986acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __END__; 4996acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5006acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( ippmomentstate ) 5016acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn icvMomentFree_64f_p( ippmomentstate ); 5026acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5036acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn cvFree( &tiles ); 5046acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn} 5056acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5066acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn/*F/////////////////////////////////////////////////////////////////////////////////////// 5076acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Name: cvGetHuMoments 5086acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Purpose: Returns Hu moments 5096acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Context: 5106acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Parameters: 5116acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// mState - moment structure filled by one of the icvMoments[Binary]*** function 5126acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// HuState - pointer to output structure containing seven Hu moments 5136acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Returns: 5146acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// CV_NO_ERR if success or error code 5156acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Notes: 5166acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn//F*/ 5176acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennCV_IMPL void 5186acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RenncvGetHuMoments( CvMoments * mState, CvHuMoments * HuState ) 5196acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn{ 5206acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_FUNCNAME( "cvGetHuMoments" ); 5216acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5226acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __BEGIN__; 5236acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5246acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !mState || !HuState ) 5256acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_ERROR_FROM_STATUS( CV_NULLPTR_ERR ); 5266acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5276acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 5286acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double m00s = mState->inv_sqrt_m00, m00 = m00s * m00s, s2 = m00 * m00, s3 = s2 * m00s; 5296acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5306acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double nu20 = mState->mu20 * s2, 5316acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn nu11 = mState->mu11 * s2, 5326acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn nu02 = mState->mu02 * s2, 5336acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn nu30 = mState->mu30 * s3, 5346acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn nu21 = mState->mu21 * s3, nu12 = mState->mu12 * s3, nu03 = mState->mu03 * s3; 5356acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5366acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double t0 = nu30 + nu12; 5376acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double t1 = nu21 + nu03; 5386acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5396acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double q0 = t0 * t0, q1 = t1 * t1; 5406acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5416acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double n4 = 4 * nu11; 5426acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double s = nu20 + nu02; 5436acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double d = nu20 - nu02; 5446acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5456acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn HuState->hu1 = s; 5466acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn HuState->hu2 = d * d + n4 * nu11; 5476acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn HuState->hu4 = q0 + q1; 5486acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn HuState->hu6 = d * (q0 - q1) + n4 * t0 * t1; 5496acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5506acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn t0 *= q0 - 3 * q1; 5516acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn t1 *= 3 * q0 - q1; 5526acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5536acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn q0 = nu30 - 3 * nu12; 5546acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn q1 = 3 * nu21 - nu03; 5556acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5566acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn HuState->hu3 = q0 * q0 + q1 * q1; 5576acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn HuState->hu5 = q0 * t0 + q1 * t1; 5586acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn HuState->hu7 = q1 * t0 - q0 * t1; 5596acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 5606acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5616acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __END__; 5626acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn} 5636acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5646acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5656acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn/*F/////////////////////////////////////////////////////////////////////////////////////// 5666acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Name: cvGetSpatialMoment 5676acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Purpose: Returns spatial moment(x_order, y_order) which is determined as: 5686acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// m(x_o,y_o) = sum (x ^ x_o)*(y ^ y_o)*I(x,y) 5696acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 0 <= x_o, y_o; x_o + y_o <= 3 5706acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Context: 5716acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Parameters: 5726acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// mom - moment structure filled by one of the icvMoments[Binary]*** function 5736acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// x_order - x order of the moment 5746acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// y_order - y order of the moment 5756acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Returns: 5766acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// moment value or large negative number (-DBL_MAX) if error 5776acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Notes: 5786acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn//F*/ 5796acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennCV_IMPL double 5806acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RenncvGetSpatialMoment( CvMoments * moments, int x_order, int y_order ) 5816acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn{ 5826acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int order = x_order + y_order; 5836acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double moment = -DBL_MAX; 5846acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5856acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_FUNCNAME( "cvGetSpatialMoment" ); 5866acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5876acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __BEGIN__; 5886acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5896acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !moments ) 5906acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_ERROR_FROM_STATUS( CV_NULLPTR_ERR ); 5916acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( (x_order | y_order) < 0 || order > 3 ) 5926acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_ERROR_FROM_STATUS( CV_BADRANGE_ERR ); 5936acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5946acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn moment = (&(moments->m00))[order + (order >> 1) + (order > 2) * 2 + y_order]; 5956acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5966acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __END__; 5976acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 5986acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn return moment; 5996acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn} 6006acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6016acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6026acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn/*F/////////////////////////////////////////////////////////////////////////////////////// 6036acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Name: cvGetCentralMoment 6046acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Purpose: Returns central moment(x_order, y_order) which is determined as: 6056acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// mu(x_o,y_o) = sum ((x - xc)^ x_o)*((y - yc) ^ y_o)*I(x,y) 6066acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 0 <= x_o, y_o; x_o + y_o <= 3, 6076acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// (xc, yc) = (m10/m00,m01/m00) - center of gravity 6086acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Context: 6096acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Parameters: 6106acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// mom - moment structure filled by one of the icvMoments[Binary]*** function 6116acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// x_order - x order of the moment 6126acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// y_order - y order of the moment 6136acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Returns: 6146acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// moment value or large negative number (-DBL_MAX) if error 6156acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Notes: 6166acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn//F*/ 6176acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennCV_IMPL double 6186acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RenncvGetCentralMoment( CvMoments * moments, int x_order, int y_order ) 6196acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn{ 6206acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int order = x_order + y_order; 6216acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double mu = 0; 6226acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6236acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_FUNCNAME( "cvGetCentralMoment" ); 6246acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6256acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __BEGIN__; 6266acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6276acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( !moments ) 6286acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_ERROR_FROM_STATUS( CV_NULLPTR_ERR ); 6296acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( (x_order | y_order) < 0 || order > 3 ) 6306acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_ERROR_FROM_STATUS( CV_BADRANGE_ERR ); 6316acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6326acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn if( order >= 2 ) 6336acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn { 6346acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mu = (&(moments->m00))[4 + order * 3 + y_order]; 6356acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn } 6366acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn else if( order == 0 ) 6376acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mu = moments->m00; 6386acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6396acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __END__; 6406acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6416acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn return mu; 6426acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn} 6436acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6446acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6456acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn/*F/////////////////////////////////////////////////////////////////////////////////////// 6466acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Name: cvGetNormalizedCentralMoment 6476acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Purpose: Returns normalized central moment(x_order,y_order) which is determined as: 6486acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// nu(x_o,y_o) = mu(x_o, y_o)/(m00 ^ (((x_o + y_o)/2) + 1)) 6496acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// 0 <= x_o, y_o; x_o + y_o <= 3, 6506acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// (xc, yc) = (m10/m00,m01/m00) - center of gravity 6516acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Context: 6526acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Parameters: 6536acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// mom - moment structure filled by one of the icvMoments[Binary]*** function 6546acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// x_order - x order of the moment 6556acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// y_order - y order of the moment 6566acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Returns: 6576acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// moment value or large negative number (-DBL_MAX) if error 6586acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn// Notes: 6596acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn//F*/ 6606acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RennCV_IMPL double 6616acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius RenncvGetNormalizedCentralMoment( CvMoments * moments, int x_order, int y_order ) 6626acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn{ 6636acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn int order = x_order + y_order; 6646acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double mu = 0; 6656acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn double m00s, m00; 6666acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6676acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_FUNCNAME( "cvGetCentralNormalizedMoment" ); 6686acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6696acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __BEGIN__; 6706acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6716acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mu = cvGetCentralMoment( moments, x_order, y_order ); 6726acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn CV_CHECK(); 6736acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6746acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn m00s = moments->inv_sqrt_m00; 6756acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn m00 = m00s * m00s; 6766acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6776acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn while( --order >= 0 ) 6786acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn m00 *= m00s; 6796acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn mu *= m00; 6806acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6816acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn __END__; 6826acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6836acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn return mu; 6846acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn} 6856acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6866acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn 6876acb9a7ea3d7564944e12cbc73a857b88c1301eeMarius Renn/* End of file. */ 688