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 开发手册
手机阅读
赞赏支持
手机阅读

微信打赏

支付宝打赏
声明
1.本文所分享的所有需要用户下载使用的内容(包括但不限于软件、数据、图片)来自于网络或者麻辣GIS粉丝自行分享,版权归该下载资源的合法拥有者所有,如有侵权请第一时间联系本站删除。
2.下载内容仅限个人学习使用,请切勿用作商用等其他用途,否则后果自负。
3.如果本文对您有所帮助,请您对我们团队进行 打赏捐助,让我们在传播3S的路上可以走得更远,不胜感激。PS:如果你是学生党,请优先把经费用于购买学习资料 以及 与小哥哥/小姐姐约会上:-)