使用vc++6.0编写程序：（1）能够读取tif格式的数字高程数据；（2）对读取的数字高程数据进行填方挖方算法

好的，下面是一个简单的示例代码，用于读取tif格式的数字高程数据并进行填方挖方算法：

#include <iostream>
#include <gdal_priv.h>
#pragma comment(lib, "gdal_i.lib")

int main()
{
    // 打开tif文件
    GDALAllRegister();
    GDALDataset* poDataset = (GDALDataset*)GDALOpen("your_tif_file.tif", GA_ReadOnly);
    if (poDataset == NULL) {
        std::cout << "Failed to open tif file." << std::endl;
        return 1;
    }

    // 获取tif文件的宽度和高度
    int width = poDataset->GetRasterXSize();
    int height = poDataset->GetRasterYSize();

    // 读取tif文件的高程数据
    float* pBuf = new float[width * height];
    GDALRasterBand* poBand = poDataset->GetRasterBand(1);
    poBand->RasterIO(GF_Read, 0, 0, width, height, pBuf, width, height, GDT_Float32, 0, 0);

    // 进行填方挖方算法
    float sum = 0.0;
    for (int i = 0; i < width * height; i++) {
        sum += pBuf[i];
    }
    float avg = sum / (width * height);
    float volumn = 0.0;
    for (int i = 0; i < width * height; i++) {
        if (pBuf[i] > avg) {
            volumn += pBuf[i] - avg;
        } else {
            volumn -= pBuf[i] - avg;
        }
    }
    std::cout << "The volume is: " << volumn << std::endl;

    // 释放资源
    GDALClose(poDataset);
    delete[] pBuf;

    return 0;
}

需要注意的是，上述代码中使用了GDAL库来读取tif文件中的高程数据。在使用时，需要在项目的工程属性中添加相应的库文件和头文件，并将GDAL的dll文件拷贝到可执行文件所在目录下。

填方挖方算法的实现方式可以根据实际需求进行调整，上述代码中仅提供了一个简单的示例。