我试图读取存储在HGT文件中的高程数据。 据我所知,他们可以理解为二进制文件。
我发现这个线程:
如何访问在C ++ .HGT SRTM文件?
基于这个职位,我的示例代码:
#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文件是很好的,我已经有可以读取多种类型的地图文件的应用进行了测试,它包含了我所需要的高程数据。
这将读取该文件,并正确填充阵列。 每次读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;
}