我正在尝试读取包含 1501×1501 的双精度矩阵的二进制文件,并将其插入特征矩阵。这是我的代码:
#include <fstream>
#include <iostream>
#include <Eigen/Dense>
#include <string>
using namespace std;
using namespace Eigen;
int main()
{
MatrixXd B(1501, 1501);
ifstream inputFile;
double toread;
inputFile.open("/path/to/bathymetry_S1000s2500s_E65d1000s65d2500s.bin",
ios::out | ios::in | ios::binary);
if (!inputFile) {
cout << "The file can't be opened.\n";
exit(10);
} else {
for (int i2 = 0; i2 < 1501; i2++) {
for (int i1 = 0; i1 < 1501; i1++) {
inputFile.read( reinterpret_cast<char*>( &toread ),
sizeof(toread) );
inputFile >> toread;
B(i1, i2) = toread;
}
}
inputFile.close();
}
cout << "Max value:" << B.maxCoeff() << endl; // Just to check the result
cout << "Mean Value:" << B.mean() << endl; // The same
}
我的问题是,当我运行代码时,我的矩阵 B 实际上只填充了 inputFile 的第一个值,即 -4502,然后由两个 cout 给出。(矩阵的所有元素都是-4502)。如何让编译器理解我希望它在上一个值之后继续读取 inputFile,而不是在每个循环步骤从头开始?