00001 /* 00002 * Software License Agreement (BSD License) 00003 * 00004 * Copyright (c) 2010, Willow Garage, Inc. 00005 * All rights reserved. 00006 * 00007 * Redistribution and use in source and binary forms, with or without 00008 * modification, are permitted provided that the following conditions 00009 * are met: 00010 * 00011 * * Redistributions of source code must retain the above copyright 00012 * notice, this list of conditions and the following disclaimer. 00013 * * Redistributions in binary form must reproduce the above 00014 * copyright notice, this list of conditions and the following 00015 * disclaimer in the documentation and/or other materials provided 00016 * with the distribution. 00017 * * Neither the name of Willow Garage, Inc. nor the names of its 00018 * contributors may be used to endorse or promote products derived 00019 * from this software without specific prior written permission. 00020 * 00021 * THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS 00022 * "AS IS" AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT 00023 * LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS 00024 * FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL THE 00025 * COPYRIGHT OWNER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, 00026 * INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES (INCLUDING, 00027 * BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; 00028 * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER 00029 * CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT 00030 * LIABILITY, OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN 00031 * ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE 00032 * POSSIBILITY OF SUCH DAMAGE. 00033 */ 00034 00039 #include <iostream> 00040 #include <map> 00041 00042 namespace pcl { 00043 00044 inline float Narf::getDescriptorDistance(const Narf& other) const 00045 { 00046 float ret = L1_Norm(descriptor_, other.descriptor_, descriptor_size_); 00047 //float ret = Sublinear_Norm(descriptor_, other.descriptor_, descriptor_size_); 00048 ret /= descriptor_size_; 00049 //std::cerr << ret<<", ";; 00050 return ret; 00051 } 00052 00053 inline void Narf::copyToNarf36(Narf36& narf36) const 00054 { 00055 if (descriptor_size_ != 36) 00056 { 00057 std::cerr << __PRETTY_FUNCTION__ << ": descriptor size is not 36!\n"; 00058 return; 00059 } 00060 getTranslationAndEulerAngles(getInverse(transformation_), narf36.x, narf36.y, narf36.z, narf36.roll, narf36.pitch, narf36.yaw); 00061 memcpy(narf36.descriptor, descriptor_, 36*sizeof(*descriptor_)); 00062 } 00063 00064 //inline float Narf::getDescriptorDistance(const Narf& other) const 00065 //{ 00066 //float middle_value = 0.1f; 00067 //float normalization_factor1 = 1.0f/middle_value, 00068 //normalization_factor2 = 1.0f/(1.0f-middle_value); 00069 //const float* descriptor1_ptr = descriptor_; 00070 //const float* descriptor2_ptr = other.getDescriptor(); 00071 //float ret = 0; 00072 //for (int i=0; i<descriptor_size_; ++i) { 00073 //float diff = fabsf(*(descriptor2_ptr++) - *(descriptor1_ptr++)); 00075 //if (diff < middle_value) 00076 //{ 00077 //diff = diff*normalization_factor1; 00078 //diff = 0.5f*diff*diff; 00080 //} 00081 //else 00082 //{ 00083 //diff = (diff - middle_value)*normalization_factor2; 00084 //diff = 0.5f + 0.5f*diff; 00087 //} 00089 //ret += diff; 00090 //} 00091 //ret /= descriptor_size_; 00092 //return ret; 00093 //} 00094 00095 //inline float Narf::getDescriptorDistance(const Narf& other) const 00096 //{ 00097 //float max_diff_between_cells = 0.25; 00098 00099 //const float* descriptor1_ptr = descriptor_; 00100 //const float* descriptor2_ptr = other.getDescriptor(); 00101 //float ret = 0; 00102 //for (int i=0; i<descriptor_size_; ++i) { 00103 //ret += (std::min)(max_diff_between_cells, fabsf(*(descriptor2_ptr++) - *(descriptor1_ptr++))); 00104 //} 00105 //ret /= descriptor_size_*max_diff_between_cells; 00106 //return ret; 00107 //} 00108 00109 } // namespace end