C#调用GDAL进行栅格镶嵌

示例代码:

using System;
using OSGeo.GDAL;

namespace GDAL_test
{
    class Program
    {
        static void Main(string[] args)
        {
            GdalConfiguration.ConfigureGdal();
            // 注册所有的驱动
            Gdal.AllRegister();

            //打开输入数据集
            Dataset inputDataset1 = Gdal.Open("H:/DEM/test1.tif", Access.GA_ReadOnly);
            Dataset inputDataset2 = Gdal.Open("H:/DEM/test2.tif", Access.GA_ReadOnly);

            // 获取栅格大小
            int xSize = inputDataset1.RasterXSize + inputDataset2.RasterXSize;
            int ySize = Math.Max(inputDataset1.RasterYSize, inputDataset2.RasterYSize);
            int bandCount = inputDataset1.RasterCount;
            //int dataType = (int)inputDataset1.GetRasterBand(1).DataType;

            //创建输出数据集
            Driver driver = Gdal.GetDriverByName("GTiff");
            Dataset outputDataset = driver.Create("H:/DEM/output.tif", xSize, ySize, bandCount, DataType.GDT_Int16, null);

            //设置输出数据集的投影和地理变换
            string projection = inputDataset1.GetProjectionRef();
            outputDataset.SetProjection(projection);
            double[] geoTransform = new double[6];
            inputDataset1.GetGeoTransform(geoTransform);
            outputDataset.SetGeoTransform(geoTransform);

           //将数据从输入数据集复制到输出数据集
            for (int i = 1; i <= bandCount; i++)
            {
                Band inputBand1 = inputDataset1.GetRasterBand(i);
                Band outputBand = outputDataset.GetRasterBand(i);
                int[] buffer = new int[inputDataset1.RasterXSize];
                for (int j = 0; j < inputDataset1.RasterYSize; j++)
                {
                    //从第一个输入数据集中复制数据
                    inputBand1.ReadRaster(0, j, inputDataset1.RasterXSize, 1, buffer, inputDataset1.RasterXSize, 1, 0, 0);
                    outputBand.WriteRaster(0, j, inputDataset1.RasterXSize, 1, buffer, inputDataset1.RasterXSize, 1, 0, 0);
                }
            }

            for (int i = 1; i <= bandCount; i++)
            {
                   //从第二个输入数据集中复制数据
                Band inputBand2 = inputDataset2.GetRasterBand(i);
                Band outputBand = outputDataset.GetRasterBand(i);
                int[] buffer = new int[inputDataset2.RasterXSize];
                for (int j = 0; j < inputDataset2.RasterYSize; j++)
                {
                    inputBand2.ReadRaster(0, j, inputDataset2.RasterXSize, 1, buffer, inputDataset2.RasterXSize, 1, 0, 0);
                    outputBand.WriteRaster(inputDataset1.RasterXSize, j, inputDataset2.RasterXSize, 1, buffer, inputDataset2.RasterXSize, 1, 0, 0);
                }
            }
        }
    }
}

GIS结果展示:

 说明:首先要承认自己写这部分的代码是不太好的。有些细节没有考虑,我是用两个横排的数据进行镶嵌,如果换成两个竖排数据,那代码还要重新修改,没有做到一力破万法。

可参考的代码片段也就是设置输出数据集的投影和地理变换以及将数据从输入数据集复制到输出数据集。

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

三千思丶

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值