PIE创建带压缩的栅格数据集

这段时间我一直在研究如何用PIE创建带压缩的栅格数据集，由于我在比赛中使用的原始影像大小普遍都在300M以上，软件加载较慢，因此希望能对原始影像进行压缩，加快加载时间。

首先，该方法的关键是修改DataSource.DatasetFactory.CreateRasterDataset()方法中的一个参数。CreateRasterDataset()方法的参数说明请参考《PIESDK For Donet V5.0二次开发指南》。

由于《PIESDK For Donet V5.0二次开发指南》还没来得及更新，实际开发中CreateRasterDataset()方法最后还存在一个指南中未提及的参数mList，这个参数就是用来决定创建的栅格数据集是否是带压缩的。我们只需要将其设置为 "COMPRESS=PACKBITS"即可完成我们的需求。

经过对比，两幅影像不存在失真现象（左图为原始影像，右图为压缩后的影像）。

 

 

压缩后的大小仅为原影像的一半左右。

 

源代码如下：

 

/// <summary>
        /// 创建带压缩的栅格数据集事件
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void toolStripButton1_Click(object sender, EventArgs e)
        {
            string pathA = @"D:
onghe.tif";
            string[] options = new string[1];
            options[0] = "COMPRESS=PACKBITS";
            IRasterDataset rDataset = DatasetFactory.OpenRasterDataset(pathA, OpenMode.ReadOnly);
            int xSize = rDataset.GetRasterXSize();
            int ySize = rDataset.GetRasterYSize();
            int bandCount = rDataset.GetBandCount();
            IRasterBand rasterband = rDataset.GetRasterBand(bandCount - 1);
            PixelDataType type = rasterband.GetRasterDataType();

            int[] bandMap = new int[bandCount];
            for (int i = 0; i < bandCount; i++)
            {
                bandMap[i] = i + 1;
            }

            double[] geoTransA = rDataset.GetGeoTransform();

            UInt16[] buf = new UInt16[xSize * ySize * bandCount];

            bool ok = rDataset.Read(0,0, xSize, ySize, buf, xSize, ySize, type, bandCount, bandMap);

            ISpatialReference spGLL = rDataset.SpatialReference;
            string stGLL = spGLL.ExportToWkt();

            string rPathSaveGLL = @"D:CompressedImage.tiff";
            IRasterDataset rasterDatasetA = DatasetFactory.CreateRasterDataset(rPathSaveGLL, xSize, ySize, bandCount, type, "GTIFF", options);

            rasterDatasetA.Write(0, 0, xSize, ySize, buf, xSize, ySize, type, bandCount, bandMap);
            rasterDatasetA.SpatialReference = SpatialReferenceFactory.CreateSpatialReference(stGLL);
            rasterDatasetA.SetGeoTransform(geoTransA);
            (rasterDatasetA as IDisposable).Dispose();

            mapControlMain.AddLayerFromFile(pathA, 0);
            mapControlMain.AddLayerFromFile(rPathSaveGLL, 1);
            mapControlMain.PartialRefresh(ViewDrawPhaseType.ViewAll);
        }