ArcGIS Engine 10 开发手册(8-9)示例:栅格数据分块

本文目录
  • 正文

我们知道栅格数据往往是很大的,为了提高处理的效率,我们有的时候需要将栅格数据分块,处理每 一个小块的数据,然后在需要的时候进行合并,当然 ArcGIS 的工具箱直接提供了这个工具,下面是我们 通过代码实现类似的功能。

/// <summary>

/// 分割栅格数据

/// </summary>

/// <param name="pRasterDataset"></param>

/// <param name="pOutputWorkspace"></param>

/// <param name="pWidth"></param>

/// <param name="pHeight"></param>

public void CreateTilesFromRasterDataset (IRasterDataset pRasterDataset, IWorkspace pOutputWorkspace, int pWidth, int pHeight)

{

  IRasterProps pRasterProps = (IRasterProps) pRasterDataset.CreateDefaultRaster ();

  double xTileSize = pRasterProps.MeanCellSize ().X * pWidth;
  double yTileSize = pRasterProps.MeanCellSize ().Y * pHeight;

  int xTileCount = (int) Math.Ceiling ((double) pRasterProps.Width / pWidth);
  int yTileCount = (int) Math.Ceiling ((double) pRasterProps.Height / pHeight);

  IEnvelope pExtent = pRasterProps.Extent;
  IEnvelope pTileExtent = new EnvelopeClass ();
  ISaveAs pSaveAs = null;

  for (int i = 0; i < xTileCount; i++)

  {

    for (int j = 0; j < yTileCount; j++)

    {

      pRasterProps = (IRasterProps) pRasterDataset.CreateDefaultRaster ();

      pTileExtent.XMin = pExtent.XMin + i * xTileSize;
      pTileExtent.XMax = pTileExtent.XMin + xTileSize;

      pTileExtent.YMin = pExtent.YMin + j * yTileSize;
      pTileExtent.YMax = pTileExtent.YMin + yTileSize;

      pRasterProps.Height = pHeight;
      pRasterProps.Width = pWidth;
      pRasterProps.Extent = pTileExtent;

    }

    pSaveAs = (ISaveAs) pRasterProps;

    pSaveAs.SaveAs ("tile_" + i + "_" + j + ".tif", pOutputWorkspace, "TIFF");

  }

}

ArcGIS Engine 10 开发手册全集

ArcGIS Engine 10 开发手册全集: ArcGIS Engine 10 开发手册

如果本文对你有所帮助,欢迎对我们团队进行打赏捐助,让我们在传播3S的路上可以走得更远。
赞赏支持
微信捐助麻辣GIS 支付宝捐助麻辣GIS

麻辣GIS-Sailor
作者:
GIS爱好者,学GIS,更爱玩GIS。

如您有疑问,可在文末留言,也可在 麻辣GIS•小分队 搜索回答或发帖 ,还可以在QQ群中提问。

本站QQ群(一):291616564 麻辣GIS

本站QQ群(二):166408035 麻辣GIS(二)

微信公众号:malagis,扫描右边二维码直接关注。

发表评论