Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035 #ifndef OPENMS_FORMAT_SPECARRAYFILE_H
00036 #define OPENMS_FORMAT_SPECARRAYFILE_H
00037
00038 #include <OpenMS/DATASTRUCTURES/String.h>
00039 #include <OpenMS/CONCEPT/Exception.h>
00040 #include <OpenMS/CONCEPT/LogStream.h>
00041 #include <OpenMS/KERNEL/Feature.h>
00042 #include <OpenMS/FORMAT/TextFile.h>
00043
00044 #include <fstream>
00045 #include <vector>
00046
00047 namespace OpenMS
00048 {
00061 class OPENMS_DLLAPI SpecArrayFile
00062 {
00063 public:
00065 SpecArrayFile();
00067 virtual ~SpecArrayFile();
00068
00077 template <typename FeatureMapType>
00078 void load(const String & filename, FeatureMapType & feature_map)
00079 {
00080
00081 TextFile input(filename, false);
00082
00083
00084 FeatureMapType fmap;
00085 feature_map = fmap;
00086
00087 for (Size i = 1; i < input.size(); ++i)
00088 {
00089 String line = input[i];
00090
00091 std::vector<String> parts;
00092 line.split('\t', parts);
00093
00094 if (parts.size() < 5)
00095 {
00096 throw Exception::ParseError(__FILE__, __LINE__, __PRETTY_FUNCTION__, "", String("Failed to convert line") + String(i + 1) + "not enough columns (expected 5 or more, got " + String(parts.size()) + ")");
00097 }
00098
00099 Feature f;
00100 try
00101 {
00102 f.setMZ(parts[0].toDouble());
00103 f.setRT(parts[1].toDouble() * 60.0);
00104 f.setMetaValue("s/n", parts[2].toDouble());
00105 f.setCharge(parts[3].toInt());
00106 f.setIntensity(parts[4].toDouble());
00107 }
00108 catch (Exception::BaseException )
00109 {
00110 throw Exception::ParseError(__FILE__, __LINE__, __PRETTY_FUNCTION__, "", String("Failed to convert value into a number (line '") + (i + 1) + ")");
00111 }
00112 feature_map.push_back(f);
00113 }
00114 }
00115
00123 template <typename SpectrumType>
00124 void store(const String & filename, const SpectrumType & spectrum) const
00125 {
00126 std::cerr << "Store() for SpecArrayFile not implemented. Filename was: " << filename << ", spec of size " << spectrum.size() << "\n";
00127 throw Exception::NotImplemented(__FILE__, __LINE__, __PRETTY_FUNCTION__);
00128 }
00129
00130 };
00131 }
00132
00133 #endif // OPENMS_FORMAT_SPECARRAYFILE_H