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:\ronghe.tif";
             string[] options = new string[];
             options[] = "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 - );
             PixelDataType type = rasterband.GetRasterDataType();

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

             double[] geoTransA = rDataset.GetGeoTransform();

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

             bool ok = rDataset.Read(,, 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(, , xSize, ySize, buf, xSize, ySize, type, bandCount, bandMap);
             rasterDatasetA.SpatialReference = SpatialReferenceFactory.CreateSpatialReference(stGLL);
             rasterDatasetA.SetGeoTransform(geoTransA);
             (rasterDatasetA as IDisposable).Dispose();

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