ArcGIS Engine 10 开发手册(8-10)示例：模拟栅格计算器

     ArcGIS 栅格数据分析的强大不仅仅在于 ArcGIS 直接提供了众多的分析工具，而且还提供了一个栅格 计算器，利用栅格计算器我们可以结合自己的模型，通过一定的脚本去执行，这让我们的分析更加灵活， ArcGIS Engine 要实现这样一个功能，只需要 IMapAlgebraOp 接口，该接口有以下几个属性和方法：

     其中 BindRaster 和 UnbindRaster 是一组相反的操作，分别用于将栅格数据和一个符号绑定和分离，下面，而 Execute 方法是用来执行我们的那个计算表达式，下面我们实现一个将栅格数据的数值除以 1000 的计算.

    代码如下

public static void UsingRasterMapAlgebra (string pFileName, string pValue)

{

  string pFolderName = Path.GetDirectoryName (pFileName);

  IRasterDataset pRaster = OpenFileRasterDataset (pFolderName, Path.GetFileNameWithoutExtension (pFileName));

  IMapAlgebraOp pMapAlgebraOp;

  pMapAlgebraOp = new RasterMapAlgebraOpClass ();

  IRasterAnalysisEnvironment pEnv = default (IRasterAnalysisEnvironment);

  pEnv = (IRasterAnalysisEnvironment) pMapAlgebraOp;

  IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass ();

  IWorkspace pWorkspace = workspaceFactory.OpenFromFile (pFolderName, 0);

  pEnv.OutWorkspace = pWorkspace;

  pMapAlgebraOp.BindRaster ((IGeoDataset) pRaster, "Ras01");

  string sOut = "[Ras01] " + "/ " + pValue;

  IRaster rasOut = (IRaster) pMapAlgebraOp.Execute (sOut);

  //保存 "[Ras01] / 1000" ISaveAs2 pSaveAs;

  pSaveAs = (ISaveAs2) rasOut;

  pSaveAs.SaveAs (pFileName + DateTime.Now.Date.Minute.ToString (), pWorkspace, "GRID");

}