0

我正在尝试读取包含 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,而不是在每个循环步骤从头开始?

4

0 回答 0