DEM 文件填挖方体积计算 C 语言代码

#include <stdio.h> #include <stdlib.h>

#define PIXEL_SIZE 30.0 // 像素大小 #define NO_DATA_VALUE -9999 // 无效像素值

int main() { FILE *fp_dem; char *filename = "dem.tif"; int width, height, i, j; float *elevation; float fill_volume = 0.0, cut_volume = 0.0;

// 打开 DEM.tif 文件
fp_dem = fopen(filename, "rb");
if (fp_dem == NULL) {
    printf("Error: cannot open file %s!

", filename); return 1; }

// 读取 DEM.tif 文件头信息
fseek(fp_dem, 0x12, SEEK_SET);
fread(&width, 4, 1, fp_dem);
fread(&height, 4, 1, fp_dem);

// 计算 DEM.tif 文件大小
fseek(fp_dem, 0, SEEK_END);
long filesize = ftell(fp_dem);
long datasize = filesize - 0x44;

// 分配内存空间
elevation = (float *)malloc(datasize);

// 读取 DEM.tif 数据
fseek(fp_dem, 0x44, SEEK_SET);
fread(elevation, sizeof(float), datasize/sizeof(float), fp_dem);

// 遍历像素值，计算填方和挖方的体积
for (i = 0; i < height; i++) {
    for (j = 0; j < width; j++) {
        float elev1 = elevation[i*width+j];
        if (elev1 == NO_DATA_VALUE) {
            continue;
        }
        float elev2 = NO_DATA_VALUE;
        if (i > 0) {
            elev2 = elevation[(i-1)*width+j];
            if (elev2 != NO_DATA_VALUE) {
                float diff = elev1 - elev2;
                if (diff > 0) {
                    cut_volume += diff * PIXEL_SIZE * PIXEL_SIZE;
                } else {
                    fill_volume += (-diff) * PIXEL_SIZE * PIXEL_SIZE;
                }
            }
        }
        if (i < height-1) {
            elev2 = elevation[(i+1)*width+j];
            if (elev2 != NO_DATA_VALUE) {
                float diff = elev1 - elev2;
                if (diff > 0) {
                    cut_volume += diff * PIXEL_SIZE * PIXEL_SIZE;
                } else {
                    fill_volume += (-diff) * PIXEL_SIZE * PIXEL_SIZE;
                }
            }
        }
        if (j > 0) {
            elev2 = elevation[i*width+j-1];
            if (elev2 != NO_DATA_VALUE) {
                float diff = elev1 - elev2;
                if (diff > 0) {
                    cut_volume += diff * PIXEL_SIZE * PIXEL_SIZE;
                } else {
                    fill_volume += (-diff) * PIXEL_SIZE * PIXEL_SIZE;
                }
            }
        }
        if (j < width-1) {
            elev2 = elevation[i*width+j+1];
            if (elev2 != NO_DATA_VALUE) {
                float diff = elev1 - elev2;
                if (diff > 0) {
                    cut_volume += diff * PIXEL_SIZE * PIXEL_SIZE;
                } else {
                    fill_volume += (-diff) * PIXEL_SIZE * PIXEL_SIZE;
                }
            }
        }
    }
}

// 输出结果
printf("Fill volume = %.2f m3\n", fill_volume);
printf("Cut volume = %.2f m3\n", cut_volume);

// 释放内存空间
free(elevation);

// 关闭文件
fclose(fp_dem);

return 0;

}

invalid file size 的解决方法:

如果程序运行时出现 "invalid file size" 的错误，可能是因为读取的文件大小不正确。可以尝试以下几种方法：

检查文件路径和文件名是否正确。
检查读取文件的方式是否正确。例如，如果文件是以二进制方式存储的，应该使用 "rb" 模式打开文件。
检查读取文件时使用的数据类型和数据大小是否正确。例如，如果文件中存储的是4字节的浮点数，应该使用 "sizeof(float)" 作为数据大小。
如果文件大小确实不正确，可以尝试重新生成文件或者修改文件大小。例如，可以使用图像处理软件重新生成 DEM.tif 文件，或者使用文件编辑器修改文件大小。