当前位置:首页 > 代码相关 > C++ > 正文内容

C++读取ASCII格式DEM数据

admin12个月前 (10-03)C++318

要使用OpenGL显示DEM,就把tif的DEM转成了ASCII格式的,用的是ArcMap工具箱里的栅格转ASCII,转出来的格式是这样的:用了6行来描述数据,如下:

第一行:数据的列数ncols

第二行:数据的行数 nrows

第三行:数据左下角的X值xllcorner

第四行:数据左下角的Y值yllcorner

第五行:像元大小cellsize

第六行:无数据的值NODATA_value

第七行开始就是数据了,按照行和列的位置排列,值代表像元值的大小。

关于一些DEM格式的转换可以看这里

下面是ASCII格式DEM的截图:

image.png

这里的思路是,先读取第2个和第4个内容也就是列和行到int变量里,然后跳过后面的8个内容,从第13个内容开始就是数据了,开始读到float数组里,当然这是看的人家的代码,应该还有更好的方式来读取。这里的方式也相应的存在一些问题:如果数据没有空间参考,那么出来的ASCII文件没有数据左上角的坐标,读取时数据的起始位置会发生偏差,会出现问题。数据量较大的时候这种方式也会出现问题:需要开辟的数组太大了。当然这里先不管那些,问题留到以后解决。

代码:

#include <iostream>
#include <fstream>
#include <string>

using namespace std;

int main()
{
    float *H;
    float Z;
    int col, row;

    ifstream infile;
    infile.open("demclip1.txt");
    if (!infile)
        cout << "error" << endl;

    string str;

    infile >> str >> col >> str >> row;

    int datanum = col * row;
    H = new float[datanum];
    for (int i = 0; i < 8; i++)
    {
        infile >> str;
        cout << str << endl;
    }

    for (int i = 0; i < col * row; i++)
    {
        infile >> Z;
        if (Z <= -9999 || Z >= 9999)
            H[i] = 0;
        else
            H[i] = Z;
    }

    cin.get();
    return 0;
}

 

 

版权声明:本文由cyhu's essay发布,如需转载请注明出处。

本文链接:https://lovedm.club/?id=91

相关文章

OpenCV中简单的添加噪声

随机数生成器用到的是C++11中的新方法。以下代码来自OpenCV计算机视觉编程攻略第三版2.2.2节void salt(Mat image, int n)...

C++ vector

vector是一个容器,类似于数组,使用的顺序存储结构,可以看做是可以动态拓展的数组,使用vector的优点是其以常数时间进行索引,缺点是在中间进行插入和删除操作代价比较高,但是如果插入和删除发生在尾...

C++ 函数指针

1、什么是函数指针函数是实现特定功能的程序代码的集合,函数代码在内存中也要占据一段存储空间,这段空间的起始地址称为函数的入口地址。C++规定函数的入口地址为函数的指针,即函数名既代表函数,又是函数的指...

C++ 拷贝构造函数

拷贝构造函数是一种特殊的构造函数,它在创建对象时,是使用同一类中之前创建的对象来初始化新创建的对象。通常用于:通过使用另一个同类型的对象来初始化新创建的对象。复制对象把它作为参数传递给函数。复制对象,...

C++计算程序运行的时间

最近在程序中有个读取文件的操作,想知道耗费的时间,查找之后大家都在用C++11中的库chrono,找了个能运行的现成写好的类:链接代码如下:#include <iostream>...

C++常量

关于常量的一些形式以及含义在此记录一下:#include <iostream> using namespace std; int main(...