0

HGT ファイルに保存されている標高データを読み込もうとしています。私の知る限り、バイナリ ファイルとして読み取ることができます。

このスレッドを見つけました:
How do I access .HGT SRTM files in C++?

その投稿に基づいて、私のサンプルコードは次のとおりです。

#include <iostream>
#include <fstream>

int main(int argc, const char * argv[])
{

std::ifstream::pos_type size;
char * memblock;

std::ifstream file ("N45W066.hgt", std::ios::in|std::ios::binary|std::ios::ate);

if (file.is_open())
{
    size = 2;
    memblock = new char [size];

    file.seekg(0, std::ios::beg);
    file.read(memblock, size);

    int srtm_ver = 1201;
    int height[1201][1021];

    for (int i = 0; i<srtm_ver; ++i){
        for (int j = 0; j < srtm_ver; ++j) {

            height[i][j] = (memblock[0] << 8 | memblock[1]);
            std::cout<<height[i][j]<<" ";
        }
        std::cout<<std::endl;
    }
}


return 0;
}

最初の実行後、たくさんのゼロが表示され、他には何も表示されません:| hgt ファイルは適切です。いくつかの種類のマップ ファイルを読み取ることができるアプリケーションでテストしました。必要な標高データが含まれています。

4

1 に答える 1

3

これにより、ファイルが読み取られ、配列が正しく読み込まれます。一度に 2 バイトを読み取るのは、一般的に最も効率的な方法ではありませんが、簡単です。別の方法は、ファイル全体を読み取ってから、後でバイトを交換することです。

Visual Studio の既定のスタック サイズでのスタック オーバーフローの問題を回避するために、高さの配列をメインの外に移動しました。スタックが十分に大きい場合は、元に戻すか、ヒープにメモリを動的に割り当てることができます。

#include <iostream>
#include <fstream>

const int SRTM_SIZE = 1201;
short height[SRTM_SIZE][SRTM_SIZE] = {0};

int main(int argc, const char * argv[])
{
    std::ifstream file("N45W066.hgt", std::ios::in|std::ios::binary);
    if(!file)
    {
        std::cout << "Error opening file!" << std::endl;
        return -1;
    }

    unsigned char buffer[2];
    for (int i = 0; i < SRTM_SIZE; ++i)
    {
        for (int j = 0; j < SRTM_SIZE; ++j) 
        {
            if(!file.read( reinterpret_cast<char*>(buffer), sizeof(buffer) ))
            {
                std::cout << "Error reading file!" << std::endl;
                return -1;
            }
            height[i][j] = (buffer[0] << 8) | buffer[1];
        }
    }

    //Read single value from file at row,col
    const int row = 500;
    const int col = 1000;
    size_t offset = sizeof(buffer) * ((row * SRTM_SIZE) + col);
    file.seekg(offset, std::ios::beg);
    file.read( reinterpret_cast<char*>(buffer), sizeof(buffer) );
    short single_value = (buffer[0] << 8) | buffer[1];
    std::cout << "values at " << row << "," << col << ":" << std::endl;
    std::cout << "  height array: " << height[row][col] << ", file: " << single_value << std::endl;

    return 0;
}
于 2013-04-30T11:39:45.133 に答える