0

我对英语的理解不太好，而且我不是最好的 C++ 程序员，所以我决定在这里问一个问题来弥补这一点。

我正在尝试从我从该页面下载的 GeoTiff 图像中获取海拔高度：http: //terracolor.net/sample_imagery.html

我通过 GDAL lib 从图像中获取了一些信息，但我不确定这到底是什么。是颜色吗？这是坐标吗？那是什么团结？我也对我能读到的不同乐队感到困惑。也许过去使用过 GDAL 的人可以向我解释这是怎么回事。最后我喜欢的是给定像素的高度，例如米。

这是我的代码：

int ofApp::getAlt(int x,int y){ 
GDALDataset  *poDataset;

GDALAllRegister();

poDataset = (GDALDataset *) GDALOpen( "data/test.tif", GA_ReadOnly );
if( poDataset == NULL )
{
    cout << "no" << endl;
}else{
    GDALRasterBand  *poBand;
    int             nBlockXSize, nBlockYSize;
    int             bGotMin, bGotMax;
    double          adfMinMax[2];

    //printf( "Size is %dx%dx%d\n", poDataset->GetRasterXSize(), poDataset->GetRasterYSize(), poDataset->GetRasterCount() );

    poBand = poDataset->GetRasterBand( 3 );
    poBand->GetBlockSize( &nBlockXSize, &nBlockYSize );

    adfMinMax[0] = poBand->GetMinimum( &bGotMin );
    adfMinMax[1] = poBand->GetMaximum( &bGotMax );
    if( ! (bGotMin && bGotMax) )
        GDALComputeRasterMinMax((GDALRasterBandH)poBand, TRUE, adfMinMax);

    float *pafScanline;
    int   nXSize = poBand->GetXSize();

    pafScanline = (float *) CPLMalloc(sizeof(float)*nXSize);
    poBand->RasterIO( GF_Read, x, y, 1, 1, 
                      pafScanline, nXSize, 1, GDT_Float32, 
                      0, 0 );

    //cout << "vvv" << pafScanline[0] << endl;
    //printf( "value %f \n", pafScanline[0]);
    return pafScanline[0];
 }
 }
4

2 回答 2

0

你说的对。多么简单的错误。我使用的这些图像中没有像带这样的高程数据。我现在使用来自http://srtm.csi.cgiar.org/SELECTION/inputCoord.asp的图像，这些图像以我需要的方式提供所需的数据。然后我可以使用我的函数来读出它。

于 2014-05-28T09:14:14.733 回答
0

Terracolor 提供的这些 GeoTIFF 文件似乎包含颜色信息，而不是高程。这些文件在 RGB 空间中具有像素颜色，共有三个波段（1：红色，2：绿色，3：蓝色）。您可以使用与您的示例类似的代码处理 GeoTIFF 文件，但您需要阅读三个波段。数据类型是字节，而不是 GDT_Float32。在使用 GDAL 处理它们之前，我建议与数据提供者核实或使用 GIS 软件（如 QGIS）打开其中一些栅格图层。.tif 文件似乎带有 .txt 文件中的元数据。

于 2014-05-25T15:52:09.720 回答