c语言读取dem数据获取高程值计算填挖方不用第三方库dem列486行327格式tiff像元大小55类型浮点型深度32位

由于TIFF格式的DEM数据比较复杂，需要使用一些第三方库来读取，因此本回答将使用其他格式的DEM数据进行演示。

假设我们有一个名为"dem.txt"的文本文件，其中每一行代表一行DEM数据，每个数据以空格分隔，如下所示：

10.2 10.1 10.5 10.3 10.2
10.3 10.2 10.6 10.4 10.3
10.4 10.3 10.7 10.5 10.4
10.5 10.4 10.8 10.6 10.5

我们可以使用以下代码读取该文件并计算填挖方：

#include <stdio.h>

int main()
{
    // 打开文件
    FILE *fp = fopen("dem.txt", "r");
    if (fp == NULL)
    {
        printf("无法打开文件！");
        return -1;
    }

    // 读取DEM数据
    float dem[4][5];
    for (int i = 0; i < 4; i++)
    {
        for (int j = 0; j < 5; j++)
        {
            fscanf(fp, "%f", &dem[i][j]);
        }
    }

    // 计算填挖方
    float total_volume = 0;
    for (int i = 0; i < 3; i++)
    {
        for (int j = 0; j < 4; j++)
        {
            // 计算当前位置和右侧位置的高程差
            float height_diff = dem[i][j+1] - dem[i][j];

            // 如果高程差为正，则需要填方；否则需要挖方
            if (height_diff > 0)
            {
                total_volume += height_diff * 5 * 5;
            }
            else
            {
                total_volume -= height_diff * 5 * 5;
            }

            // 计算当前位置和下方位置的高程差
            height_diff = dem[i+1][j] - dem[i][j];

            // 如果高程差为正，则需要填方；否则需要挖方
            if (height_diff > 0)
            {
                total_volume += height_diff * 5 * 5;
            }
            else
            {
                total_volume -= height_diff * 5 * 5;
            }
        }
    }

    // 输出结果
    printf("填挖方量为：%.2f立方米", total_volume);

    // 关闭文件
    fclose(fp);

    return 0;
}

该程序首先打开文件"dem.txt"，然后使用双重循环读取DEM数据并计算填挖方。其中，每个像元的大小为5m×5m，因此需要将高程差乘以25来得到填挖方量。最后，程序输出填挖方量，并关闭文件。

请注意，由于DEM数据通常比较大，因此在实际应用中可能需要使用分块读取的方式来处理。此外，如果使用TIFF格式的DEM数据，则需要使用第三方库（如GDAL）来读取数据