C#+GDAL影像处理笔记06：获取两幅tiff影像的差值，并将大于某个数值的像素导出为shp文件

**C# + GDAL 影像处理笔记06**

**获取两幅 TIFF 影像的差值，并将大于某个数值的像素导出为 SHP 文件**

在前面的几篇笔记中，我们已经学习了如何使用 C# 和 GDAL 库进行影像处理。今天我们要讨论一个更复杂的问题：如何获取两幅 TIFF 影像的差值，并将大于某个数值的像素导出为 SHP 文件。

**问题描述**

假设我们有两幅 TIFF 影像，分别是 `image1.tif` 和 `image2.tif`。我们需要计算这两幅影像之间的差值，即 `image1.tif` 中的像素值减去 `image2.tif` 中的像素值。然后，我们需要将大于某个数值（例如100）的像素导出为 SHP 文件。

**解决方案**

为了解决这个问题，我们可以使用 GDAL 库提供的几种功能：

1. `GDAL.AllBandRead()`：读取多幅影像的所有波段。
2. `GDAL.Calculate()`：计算两幅影像之间的差值。
3. `GDAL.WriteRaster()`：将结果写入新的 TIFF 文件中。

我们还需要使用 C# 的 LINQ 库来过滤大于某个数值的像素。

**代码示例**

csharpusing System;
using System.IO;
using OSGeo.GDAL;

class Program{
 static void Main(string[] args)
 {
 //读取两幅 TIFF 影像 string image1Path = "image1.tif";
 string image2Path = "image2.tif";

 GDAL.AllBandRead(image1Path, out _);
 GDAL.AllBandRead(image2Path, out _);

 // 计算两幅影像之间的差值 double[] difference = new double[3];
 for (int i =0; i < 3; i++)
 {
 difference[i] = GDAL.Calculate(image1Path, image2Path, i);
 }

 // 过滤大于某个数值的像素 int threshold =100;
 var filteredPixels = new System.Collections.Generic.List<double>();
 for (int i =0; i < difference.Length; i++)
 {
 if (Math.Abs(difference[i]) > threshold)
 {
 filteredPixels.Add(difference[i]);
 }
 }

 // 将过滤后的像素写入 SHP 文件中 string shpPath = "filtered_pixels.shp";
 using (var writer = new ShapefileWriter(shpPath))
 {
 writer.Write(filteredPixels);
 }

 Console.WriteLine("完成!");
 }
}