00001 /****************************************************************************** 00002 * __ _ _ _ _ _ _____ * 00003 * | \ / | | | | | | | | | | ___| * 00004 * | \/ | | | | | | | | | | |_ * 00005 * | | | | | | | | | | | _| * 00006 * | |\/| | | | | |__ | |__ | | | |__ * 00007 * |_| |_| |_| |____| |____| |_| |____| * 00008 * __________________________________________________________________________ * 00009 * Multifunctional Library For Image Processing * 00010 * * 00011 * * 00012 * * 00013 * (c) Copyright 2007 by Humbert Florent * 00014 * * 00015 * This program is free software; you can redistribute it and/or modify * 00016 * it under the terms of the GNU General Public License as published by * 00017 * the Free Software Foundation; only version 2 of the License. * 00018 * * 00019 * This program is distributed in the hope that it will be useful, * 00020 * but WITHOUT ANY WARRANTY; without even the implied warranty of * 00021 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * 00022 * GNU General Public License for more details. * 00023 * * 00024 * You should have received a copy of the GNU General Public License * 00025 * along with this program; if not, write to the Free Software * 00026 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA * 00027 * 02111-1307, USA. * 00028 ******************************************************************************/ 00029 00032 #include "HypersurfaceFlow.hpp" 00033 #include "operator/DifferentialOperators.hpp" 00034 #include "MillieStd.hpp" 00035 00036 namespace Millie 00037 { 00038 00039 HypersurfaceFlow::HypersurfaceFlow() 00040 {} 00041 00042 float HypersurfaceFlow::computeFlow(const Image & image, int i, int j, int canal) const 00043 { 00044 00045 float opttempo = 1.0f + diffSquareGradN(image, i,j, canal); 00046 float dnn = diffDnnC(image, i, j, canal); 00047 float dnt = diffDntntC(image, i, j, canal); 00048 float cnt = 1.0f / (2.0f * std::sqrt(opttempo)); 00049 float cn = 1.0f / (2.0f * opttempo * std::sqrt(opttempo)); 00050 00051 return((cn * dnn) + (cnt * dnt)); 00052 } 00053 00054 HypersurfaceFlow * HypersurfaceFlow::clone() const 00055 { 00056 return new HypersurfaceFlow(*this); 00057 } 00058 00059 00060 00061 00062 } 00063 00064
1.5.1