高程值的二维数组生成tiff栅格文件格式

把一个含有高程值的二维数组保存为tiff栅格文件格式

在C#中，将含有高程值的二维数组（假设这些数组是从某种VCT格式或其他格式中获取的Grid数据）保存为带有坐标参考系统的TIFF栅格文件格式，你可以使用GDAL库。不过，需要注意的是，VCT本身并不是GDAL直接支持的一种格式，但我们可以假设你已经有了从VCT或其他格式中提取的高程数据和相应的坐标信息。

以下是一个简化的步骤和示例代码，说明如何使用GDAL在C#中将高程值的二维数组保存为带有坐标参考系统的TIFF文件：

安装GDAL.NET：首先，确保你的C#项目中已经安装了GDAL.NET NuGet包。
准备高程数据和坐标信息：你需要有一个二维数组来存储高程值，以及关于这些高程值如何映射到地理空间的信息（如坐标参考系统和地理变换）。
创建TIFF文件：使用GDAL的Driver来创建一个新的TIFF文件，并设置其坐标参考系统和地理变换。
写入高程数据：将高程数据写入TIFF文件的第一个波段。
保存并关闭文件。
以下是一个示例代码：

using OSGeo.GDAL;  
using OSGeo.OSR;  
using System;  
  
class Program  
{  
    static void Main(string[] args)  
    {  
        // 假设你有一个二维数组，代表高程值  
        float[,] heightData = new float[100, 100];  
        // ... 填充heightData数组 ...  
  
        // 坐标参考系统（CRS），这里以WGS 84为例  
        SpatialReference srs = new SpatialReference(null);  
        srs.ImportFromEPSG(4326); // WGS 84  
  
        // 创建TIFF文件  
        Driver driver = Gdal.GetDriverByName("GTiff");  
        string filename = "elevation.tif";  
        Dataset dataset = driver.CreateCopy(filename, null, 0, null, 100, 100, 1, DataType.GDT_Float32, null);  
  
        // 设置坐标参考系统  
        dataset.SetProjection(srs.ExportToWkt());  
  
        // 设置地理变换信息  
        // 假设每个像素代表1度（这只是一个示例，实际值应根据你的数据确定）  
        // 并且图像的左上角是(0, 90)（北极附近的某个点，注意Y坐标是向下的）  
        double[] adfGeoTransform = new double[6];  
        adfGeoTransform[0] = 0;    // 左上角x  
        adfGeoTransform[1] = 1;    // x像素大小  
        adfGeoTransform[2] = 0;    // x旋转（通常为0）  
        adfGeoTransform[3] = 90;   // 左上角y（注意：TIFF中的Y坐标是向下的）  
        adfGeoTransform[4] = 0;    // y旋转（通常为0）  
        adfGeoTransform[5] = -1;   // y像素大小（负值因为y坐标向下增加）  
        dataset.SetGeoTransform(adfGeoTransform);  
  
        // 写入高程数据到第一个波段  
        Band band = dataset.GetRasterBand(1);  
        band.WriteRaster(0, 0, 100, 100, heightData, 100, 100, 0, 0);  
  
        // 刷新并关闭Dataset  
        dataset.FlushCache();  
        dataset.Dispose();  
  
        Console.WriteLine("TIFF文件已保存，包含坐标参考系统。");  
    }  
}

注意：

上面的代码示例中，坐标变换信息（adfGeoTransform）和坐标参考系统（CRS）都是简化的。你需要根据你的实际数据来设置这些值。特别是，adfGeoTransform中的像素大小和偏移量应该反映你的高程数据的实际地理覆盖范围。
坐标参考系统（CRS）在这里被设置为WGS 84（EPSG:4326），但你应该使用与你的数据相匹配的CRS。
TIFF文件中的Y坐标是向下的，因此当你设置左上角Y坐标时，你需要注意这一点。
确保在不再需要Dataset对象时调用Dispose()方法来释放资源。

本blog地址：https://blog.csdn.net/hsg77