QVision: Qt's Image, Video and Computer Vision Library

00001 /*
00002  *      Copyright (C) 2011, 2012. PARP Research Group.
00003  *      <http://perception.inf.um.es>
00004  *      University of Murcia, Spain.
00005  *
00006  *      This file is part of the QVision library.
00007  *
00008  *      QVision is free software: you can redistribute it and/or modify
00009  *      it under the terms of the GNU Lesser General Public License as
00010  *      published by the Free Software Foundation, version 3 of the License.
00011  *
00012  *      QVision is distributed in the hope that it will be useful,
00013  *      but WITHOUT ANY WARRANTY; without even the implied warranty of
00014  *      MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00015  *      GNU Lesser General Public License for more details.
00016  *
00017  *      You should have received a copy of the GNU Lesser General Public
00018  *      License along with QVision. If not, see <http://www.gnu.org/licenses/>.
00019  */
00020 
00024 
00025 #include <QVVector>
00026 #include <QVKeypoint>
00027 #include <QFile>
00028 #include <QTextStream>
00029 #include <QVSiftFeature>
00030 
00031 bool saveSIFTDescriptorsInLoweSIFTFileFormat(const QString &fileName, const QList<QVSiftFeature> &features)
00032         {
00033         QFile file(fileName);
00034         if (not file.open(QIODevice::WriteOnly | QIODevice::Text))
00035                 return false;
00036 
00037         QTextStream stream(&file);
00038 
00039         stream.setRealNumberPrecision(16); // Default is 6, not enough precission to faithfully store certain maps.
00040 
00041         stream << features.count() << " 128\n";
00042         foreach(QVSiftFeature feature, features)
00043                 {
00044                 stream  << feature.keypoint.y() << " "
00045                         << feature.keypoint.x() << " "
00046                         << feature.keypoint.scale() << " "
00047                         << feature.keypoint.orientation();
00048 
00049                 const QVVector descriptor = 512 * feature.descriptor / feature.descriptor.norm2();
00050 
00051                 for(int i = 0; i < 128; i++)
00052                         {
00053                         if ( (i % 20) == 0 )
00054                                 stream << "\n";
00055                         stream << " " << int(descriptor[i]);
00056                         }
00057                 stream << "\n";
00058                 }
00059 
00060         file.close();
00061         return true;
00062         }
00063
src/qvip/qvsiftfeature.cpp
