两幅影像波段计算时，莫宁奇妙失败（同样的数据在arcgis中可以）。
经过艰难的长期排查（都是心酸泪--------）
发现两幅影像坐标系不一样！
ae也不给个提示
//pRasterBanda与pRasterBandb坐标系必须一致
public static IGeoDataset BandMath(IGeoDataset pRasterBanda, IGeoDataset pRasterBandb, string gongshi)
{
IMapAlgebraOp rsalgebra = new RasterMapAlgebraOpClass();
rsalgebra.BindRaster(pRasterBanda, "band1");
rsalgebra.BindRaster(pRasterBandb, "band2");
IGeoDataset pGeoDataset = rsalgebra.Execute(gongshi);//gongshi=”[band1] / [band2]”
rsalgebra.UnbindRaster("band1");
rsalgebra.UnbindRaster("band2");
return pGeoDataset;
}
----------------------------------------------
栅格处理 异常
最近做的DEM图像处理的小程序，遇到异常，整理如下：
（1）栅格求差，可以运行，但是结果很诡异，完全不是求差。多次检查栅格计算的脚本，以为是栅格计算的语句出错。多次检查运行语句，多次试验，结果还是还是很诡异。
后来发现，是栅格赋值出错。导致的结果很诡异。诡异的错误由简单的错误引起
（2）运行时出现异常，提示为Error 010327 : Unable to set analysis window 还有一个错误 是无法执行栅格计算的语句。刚开始以为是语句出错，多次修改栅格计算的语句，还是出现这个错误。后来搜索Error 010327 : Unable to set analysis window，根据官网的错误处理方法，设置了栅格计算时的cellsize,Extent等属性，运行时还是出现奇怪的错误。后来修改了其中一幅栅格图的投影坐标，两幅待处理的栅格图投影坐标不一致，但是在ArcMap里能够正常显示，是同一个处理区域，且能正常栅格计算。然后将两幅栅格图的投影坐标统一，再处理，结果正常。
http://www.cnblogs.com/DayDreamEveryWhere/archive/2012/09/02/2667754.html