public IRaster CreateRasterSlope(IRaster raster)
{
try
{
RasterSurfaceOpClass class2 = new RasterSurfaceOpClass();
IGeoDataset geoDataset = raster as IGeoDataset;
double zFactor = Convert.ToDouble(textBoxX3.Text);
if (zFactor <= 0)
{
MessageBox.Show("请指定大于0的正数！", "信息", MessageBoxButtons.OK, MessageBoxIcon.Information);
return(null);
}
//float zFactor = 1;
object o = zFactor;
esriGeoAnalysisSlopeEnum geoType = esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopeDegrees;
IGeoDataset pGeoDataset = class2.Slope(geoDataset, geoType, ref o);
IRaster pRaster = (IRaster)pGeoDataset;
ISaveAs2 pSaveAs = pRaster as ISaveAs2;
IDataset pDataset = pSaveAs.SaveAs(m_ResultPath + ".tif", null, "TIFF");
//System.Runtime.InteropServices.Marshal.ReleaseComObject(pDataset);
MessageBox.Show("坡度计算完毕！", "信息", MessageBoxButtons.OK, MessageBoxIcon.Information);
return(pRaster);
}
catch (SystemException e)
{
MessageBox.Show(e.Message, "提示", MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
return(null);
}
}
/// <summary>
/// 计算DEM坡度
/// </summary>
/// <param name="rasterLayer">输入DEM图层</param>
/// <param name="outputFullPath">输出坡度图路径，类型为IMAGE，文件名为img</param>
public void CreateRasterSlope(ILayer rasterLayer, string outputFullPath)
{
if (File.Exists(outputFullPath))
{
MessageBox.Show(outputFullPath + "已经存在，请重命名");
return;
}
IRasterLayer rasterLyr = rasterLayer as IRasterLayer;
IRaster raster = rasterLyr.Raster;
RasterSurfaceOpClass rasterOpCls = new RasterSurfaceOpClass();
object o = 1;
IGeoDataset geoDataset = raster as IGeoDataset;
esriGeoAnalysisSlopeEnum geoType = esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopeDegrees;
IGeoDataset outGeodataset = rasterOpCls.Slope(geoDataset, geoType, ref o);
IRaster pRaster = outGeodataset as IRaster;
FileInfo fi = new FileInfo(outputFullPath);
string fileDir = fi.Directory.FullName;
string name = fi.Name;
if (!Directory.Exists(fileDir))
{
Directory.CreateDirectory(fileDir);
}
IWorkspaceFactory pWSF = new RasterWorkspaceFactoryClass();
IWorkspace pWS = pWSF.OpenFromFile(fileDir, 0);
ISaveAs pSaveAs = pRaster as ISaveAs;
pSaveAs.SaveAs(name, pWS, "IMAGINE Image");
IRasterLayer rasterLyr2 = new RasterLayerClass();
rasterLyr2.CreateFromRaster(pRaster);
rasterLyr2.Name = name;
m_mapControl.Map.AddLayer(rasterLyr2 as ILayer);
}
/// <summary>
/// 计算DEM坡度
/// </summary>
/// <param name="rasterLayer">输入DEM图层</param>
/// <param name="outputFullPath">输出坡度图路径，类型为IMAGE，文件名为img</param>
public void CreateRasterSlope(ILayer rasterLayer, string outputFullPath)
{
if (File.Exists(outputFullPath))
{
MessageBox.Show(outputFullPath + "已经存在，请重命名");
return;
}
IRasterLayer rasterLyr = rasterLayer as IRasterLayer;
IRaster raster = rasterLyr.Raster;
RasterSurfaceOpClass rasterOpCls = new RasterSurfaceOpClass();
object o = 1;
IGeoDataset geoDataset = raster as IGeoDataset;
esriGeoAnalysisSlopeEnum geoType = esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopeDegrees;
IGeoDataset outGeodataset = rasterOpCls.Slope(geoDataset, geoType, ref o);
IRaster pRaster = outGeodataset as IRaster;
FileInfo fi = new FileInfo(outputFullPath);
string fileDir = fi.Directory.FullName;
string name = fi.Name;
if (!Directory.Exists(fileDir))
Directory.CreateDirectory(fileDir);
IWorkspaceFactory pWSF = new RasterWorkspaceFactoryClass();
IWorkspace pWS = pWSF.OpenFromFile(fileDir, 0);
ISaveAs pSaveAs = pRaster as ISaveAs;
pSaveAs.SaveAs(name, pWS, "IMAGINE Image");
IRasterLayer rasterLyr2 = new RasterLayerClass();
rasterLyr2.CreateFromRaster(pRaster);
rasterLyr2.Name = name;
m_mapControl.Map.AddLayer(rasterLyr2 as ILayer);
}
private void Btn_Slope_Click(object sender, EventArgs e)
{
IRasterLayer pFromRasterLayer = m_selectedLayer as IRasterLayer;
if (pFromRasterLayer == null)
{
MessageBox.Show("非栅格图层");
return;
}
IRaster pFromRaster = pFromRasterLayer.Raster;
IRasterBandCollection pFromRasterBandCollection = pFromRaster as IRasterBandCollection;
IRasterBand pFromRasterBand = pFromRasterBandCollection.Item(0);
IRasterDataset pFromRasterDataset = pFromRasterBand as IRasterDataset;
IWorkspaceFactory pWorkspaceFactory = new RasterWorkspaceFactoryClass();
IRasterWorkspace pRasterWorkspace = pWorkspaceFactory.OpenFromFile(@"DB", 0) as IRasterWorkspace;
ISurfaceOp pSurfaceOp = new RasterSurfaceOpClass();
IRasterAnalysisEnvironment pRasterAnalysisEnvironment;
pRasterAnalysisEnvironment = pSurfaceOp as IRasterAnalysisEnvironment;
pRasterAnalysisEnvironment.OutWorkspace = pWorkspaceFactory as IWorkspace;
object zFactor = new object();
IGeoDataset pGeoDataset, pRasterGetDataset;
pRasterGetDataset = pFromRasterDataset as IGeoDataset;
pGeoDataset = pSurfaceOp.Slope(pRasterGetDataset, esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopePercentrise, ref zFactor);
IRasterBandCollection pOutRasterBandCollection = pGeoDataset as IRasterBandCollection;
string strOutRasterName = pFromRasterLayer.Name.Split('.')[0] + "_slope.tif";
pOutRasterBandCollection.SaveAs(strOutRasterName, pRasterWorkspace as IWorkspace, "TIFF");
IRasterDataset pNewRasterDataset = pRasterWorkspace.OpenRasterDataset(strOutRasterName);
IRasterLayer pNewRasterLayer = new RasterLayerClass();
pNewRasterLayer.CreateFromDataset(pNewRasterDataset);
ILayer pNewLayer = pNewRasterLayer as ILayer;
Ctrl_Map.AddLayer(pNewLayer);
}
private void btnGO_Click(object sender, EventArgs e)
{
string strInFileName;
string strOutFileName;
int startX, endX;
string strOutDir;
try
{
if (bDataPath == true)
{
strInFileName = txtOutPath.Text;
strOutDir = strInFileName.Substring(0, strInFileName.LastIndexOf("\\"));
startX = strInFileName.LastIndexOf("\\");
endX = strInFileName.Length;
strOutFileName = strInFileName.Substring(startX + 1, endX - startX - 1);
}
else
{
strOutDir = txtOutPath.Text;
strOutFileName = "slope";
}
if (File.Exists(strOutDir + "\\" + strOutFileName + ".aux") == true)
{
MessageBox.Show("�ļ�" + strOutFileName + "�Ѿ����ڣ�������������");
return;
}
if (m_pRasterLyr != null)
{
double dCellSize = Convert.ToDouble(txtCellSize.Text);
ISurfaceOp pRasterSurfaceOp = new RasterSurfaceOpClass();
double dZFactor = Convert.ToDouble(txtZFactor.Text);
object objZFactor = dZFactor;
IRaster pInRaster = m_pRasterLyr.Raster;
IGeoDataset pGeoDs = null;
IRasterLayer pOutRasterLayer = new RasterLayerClass();
pOutRasterLayer.Name = strOutFileName;
if (radioDegree.Checked == true)
{
pGeoDs = pRasterSurfaceOp.Slope(pInRaster as IGeoDataset, esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopeDegrees, ref objZFactor);
}
else
{
pGeoDs = pRasterSurfaceOp.Slope(pInRaster as IGeoDataset, esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopePercentrise, ref objZFactor);
}
IRaster pOutRaster = (IRaster)pGeoDs;
pOutRasterLayer.CreateFromRaster(pOutRaster);
pOutRasterLayer.Name = strOutFileName;
IWorkspaceFactory pWSF = new RasterWorkspaceFactoryClass();
IWorkspace pRWS = pWSF.OpenFromFile(strOutDir, 0);
IRasterBandCollection pRasBandCol = (IRasterBandCollection)pGeoDs;
pRasBandCol.SaveAs(strOutFileName, pRWS, "GRID");//"IMAGINE Image"
pOutRasterLayer = UtilityFunction.SetStretchRenderer(pOutRasterLayer.Raster);
pOutRasterLayer.Name = strOutFileName;
pMap.AddLayer(pOutRasterLayer);
}
}
catch (Exception ex)
{
MessageBox.Show(ex.Message);
}
}
private IRasterLayer calWSCons(string _outPath)
{
//��դ�� pMapAlgebraOp
try
{
if (pMapAlgebraOp==null)
{
pMapAlgebraOp = new RasterMapAlgebraOpClass();
}
string strExp = "";
string sRainPath = "";
//��ֱ���ý�����ʴ�����ݻ��ǽ������ݼ���
IRaster pRasterR = null;
if (rbtnR.Checked)
{
pRasterR = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbR.Text)) as IRasterLayer).Raster;
pMapAlgebraOp.BindRaster(pRasterR as IGeoDataset, "R");
}
else
{
sRainPath = this.txtPcpPath.Text;
//
IWorkspaceFactory pWsF = new RasterWorkspaceFactoryClass();
IRasterWorkspace pRWs;
IRasterDataset pRDs = new RasterDatasetClass();
IRasterLayer pRLyr = new RasterLayerClass();
pRWs = pWsF.OpenFromFile(sRainPath, 0) as IRasterWorkspace;
string sFileName="";
for (int i = 1; i < 13; i++)
{
sFileName = this.txtPcpPrefix.Text + i.ToString() + this.txtPcpSuffix.Text;
pRDs = pRWs.OpenRasterDataset(sFileName);
pMapAlgebraOp.BindRaster(pRDs as IGeoDataset, sFileName);
strExp = strExp + "[" + sFileName + "] + ";
}
//cal total pcp
strExp = strExp.TrimEnd("+ ".ToCharArray());
IGeoDataset pGeoDsPcp = pMapAlgebraOp.Execute(strExp);
pMapAlgebraOp.BindRaster(pGeoDsPcp, "Pcp");
//cal Ri
strExp = "";
string sFileNameE = "",sDsName="",strExpR="";
for (int i = 1; i < 13; i++)
{
sFileName = this.txtPcpPrefix.Text + i.ToString() + this.txtPcpSuffix.Text;
sFileNameE="["+sFileName+"]";
strExp = "1.735 * pow(10,1.5 * log10((" + sFileNameE + " * " + sFileNameE + ") / [Pcp]) - 0.08188)";
IGeoDataset pGeoDsRi = pMapAlgebraOp.Execute(strExp);
sDsName = "Pcp" + i.ToString();
pMapAlgebraOp.BindRaster(pGeoDsRi, sDsName);
pMapAlgebraOp.UnbindRaster(sFileName);
strExpR = strExpR + "[" + sDsName + "] + ";
}
//cal R
strExpR = strExpR.TrimEnd("+ ".ToCharArray());
IGeoDataset pGeoDsR = pMapAlgebraOp.Execute(strExpR);
pMapAlgebraOp.BindRaster(pGeoDsR, "R");
for (int i = 1; i < 13; i++)
{
sDsName = "Pcp" + i.ToString();
pMapAlgebraOp.UnbindRaster(sDsName);
}
}
//����K����
IRaster pRasterSclay = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbSoilClay.Text)) as IRasterLayer).Raster;
pMapAlgebraOp.BindRaster(pRasterSclay as IGeoDataset, "clay");
IRaster pRasterSsand = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbSoilSand.Text)) as IRasterLayer).Raster;
pMapAlgebraOp.BindRaster(pRasterSsand as IGeoDataset, "sand");
IRaster pRasterSslit = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbSoilSlit.Text)) as IRasterLayer).Raster;
pMapAlgebraOp.BindRaster(pRasterSslit as IGeoDataset, "slit");
IRaster pRasterSOrg = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbSoilOrganic.Text)) as IRasterLayer).Raster;
pMapAlgebraOp.BindRaster(pRasterSOrg as IGeoDataset, "org");
//cal K
//K =(0.2 + 0.3 * exp(-0.0256 * soil_sand * (1.0 - soil_silt / 100.0))) * pow((soil_silt * 1.0 / (soil_clay * 1.0 + soil_silt * 1.0)),0.3) * (1.0 - 0.25 * soil_oc / (soil_oc * 1.0 + exp(3.72 - 2.95 * soil_oc))) * (1.0 - (0.7 * ksd) / (ksd + exp(-5.51 + 22.9 * ksd)))
strExp = "(0.2 + 0.3 * Exp(-0.0256 * [sand] * (1.0 - [slit] / 100.0))) * Pow(([slit] * 1.0 / ([clay] * 1.0 + [slit] * 1.0)), 0.3) * (1.0 - 0.25 * [org] * 0.58 / ([org] * 0.58 + Exp(3.72 - 2.95 * [org] * 0.58))) * (1.0 - (0.7 * (1.0 - [sand] / 100.0)) / ((1.0 - [sand] / 100.0) + Exp(-5.51 + 22.9 * (1.0 - [sand] / 100.0))))";
IGeoDataset pGeoDsK = pMapAlgebraOp.Execute(strExp);
pMapAlgebraOp.BindRaster(pGeoDsK, "K");
pMapAlgebraOp.UnbindRaster("clay");
pMapAlgebraOp.UnbindRaster("sand");
pMapAlgebraOp.UnbindRaster("slit");
pMapAlgebraOp.UnbindRaster("org");
//cal L*S
IHydrologyOp pHydrologyOp = new RasterHydrologyOpClass();
object objZLimit = System.Type.Missing;
IRaster pRasterDem = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbDem.Text)) as IRasterLayer).Raster;
//Fill Dem
IGeoDataset pGeoDsDemFill =pHydrologyOp.Fill(pRasterDem as IGeoDataset,ref objZLimit);
// cal FlowDirection
IGeoDataset pGeoDsFlowDir = pHydrologyOp.Fill(pGeoDsDemFill, ref objZLimit);
//cal FlowAccumulation
IGeoDataset pGeoDsFlowAcc = pHydrologyOp.Fill(pGeoDsFlowDir, ref objZLimit);
pMapAlgebraOp.BindRaster(pGeoDsFlowAcc, "FlowAcc");
//cal Slope with Deg
ISurfaceOp pSurfaceOp = new RasterSurfaceOpClass();
object objZFactor = System.Type.Missing;
IGeoDataset pGeoDsSlope = pSurfaceOp.Slope(pGeoDsDemFill, esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopeDegrees, ref objZFactor);
// bind raster data "Slope"
pMapAlgebraOp.BindRaster(pGeoDsSlope, "Slope");
//cal S
strExp = "Con([Slope] < 5,10.8 * Sin([Slope] * 3.14 / 180) + 0.03,[Slope] >= 10,21.9 * Sin([Slope] * 3.14 / 180) - 0.96,16.8 * Sin([Slope] * 3.14 / 180) - 0.5)";
IGeoDataset pGeoDsS = pMapAlgebraOp.Execute(strExp);
pMapAlgebraOp.BindRaster(pGeoDsS, "S");
//cal m
strExp = "Con([Slope] <= 1,0.2,([Slope] > 1 & [Slope] <= 3),0.3,[Slope] >= 5,0.5,0.4)";
IGeoDataset pGeoDsM = pMapAlgebraOp.Execute(strExp);
pMapAlgebraOp.BindRaster(pGeoDsM, "m");
//cal ls
strExp = "[S] * Pow(([FlowAcc] * "+this.txtCellSize.Text+" / 22.1),[m])";
IGeoDataset pGeoDsLS = pMapAlgebraOp.Execute(strExp);
pMapAlgebraOp.BindRaster(pGeoDsLS, "LS");
pMapAlgebraOp.UnbindRaster("m");
pMapAlgebraOp.UnbindRaster("S");
pMapAlgebraOp.UnbindRaster("Slope");
pMapAlgebraOp.UnbindRaster("FlowAcc");
IRaster pRasterC = null;
if (rbtnVegCover.Checked)
{
pRasterC = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbC.Text)) as IRasterLayer).Raster;
pMapAlgebraOp.BindRaster(pRasterC as IGeoDataset, "C");
}
else
{
//cal vegetation cover
IRaster pRasterNDVI = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbC.Text)) as IRasterLayer).Raster;
IRasterBandCollection pRasterBandCollection = pRasterNDVI as IRasterBandCollection;
IRasterBand pRasterBand = pRasterBandCollection.Item(0);
pRasterBand.ComputeStatsAndHist();
IRasterStatistics pRasterStatistics = pRasterBand.Statistics;
double dMax = pRasterStatistics.Maximum;
double dMin = pRasterStatistics.Minimum;
pMapAlgebraOp.BindRaster(pRasterNDVI as IGeoDataset, "NDVI");
if (dMin < 0)
{
dMin = 0;
}
//veg%yr% = (ndvi%yr% - ndvi%yr%min) / (ndvi%yr%max - ndvi%yr%min)
strExp = "([NDVI] - "+dMin+") / ("+dMax +" - "+dMin+")";
IGeoDataset pGeoDsC = pMapAlgebraOp.Execute(strExp);
pMapAlgebraOp.BindRaster(pGeoDsC, "C");
pMapAlgebraOp.UnbindRaster("NDVI");
}
//����P����
IRaster pRasterP = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbP.Text)) as IRasterLayer).Raster;
pMapAlgebraOp.BindRaster(pRasterP as IGeoDataset, "P");
if (_outPath != "")
{
strExp = _outPath + " = [R] * [K] * [LS] * (1 - [C]) * [P]";
}
else
{
strExp = "[R] * [K] * [LS] * (1 - [C]) * [P]";
}
IGeoDataset pGeoDsSr = pMapAlgebraOp.Execute(strExp);
IRaster pOutRaster = pGeoDsSr as IRaster;
IRasterLayer pOutRasterLayer = new RasterLayerClass();
pOutRasterLayer.CreateFromRaster(pOutRaster);
//��������
//string strOutDir = _outPath.Substring(0, _outPath.LastIndexOf("\\"));
int startX = _outPath.LastIndexOf("\\");
int endX = _outPath.Length;
string sRLyrName = _outPath.Substring(startX + 1, endX - startX - 1);
pOutRasterLayer.Name = sRLyrName;
pMapAlgebraOp = null;
return pOutRasterLayer;
}
catch (Exception ex)
{
MessageBox.Show(ex.Message.ToString());
return null;
}
}
private void btnGO_Click(object sender, EventArgs e)
{
string strInFileName;
string strOutFileName;
int startX, endX;
string strOutDir;
try
{
if (bDataPath == true)
{
strInFileName = txtOutPath.Text;
strOutDir = strInFileName.Substring(0, strInFileName.LastIndexOf("\\"));
startX = strInFileName.LastIndexOf("\\");
endX = strInFileName.Length;
strOutFileName = strInFileName.Substring(startX + 1, endX - startX - 1);
}
else
{
strOutDir = txtOutPath.Text;
strOutFileName = "slope";
}
if (File.Exists(strOutDir + "\\" + strOutFileName + ".aux") == true)
{
MessageBox.Show("文件" + strOutFileName + "已经存在，请重新命名！");
return;
}
if (m_pRasterLyr != null)
{
double dCellSize = Convert.ToDouble(txtCellSize.Text);
ISurfaceOp pRasterSurfaceOp = new RasterSurfaceOpClass();
double dZFactor = Convert.ToDouble(txtZFactor.Text);
object objZFactor = dZFactor;
IRaster pInRaster = m_pRasterLyr.Raster;
IGeoDataset pGeoDs = null;
IRasterLayer pOutRasterLayer = new RasterLayerClass();
pOutRasterLayer.Name = strOutFileName;
if (radioDegree.Checked == true)
{
pGeoDs = pRasterSurfaceOp.Slope(pInRaster as IGeoDataset, esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopeDegrees, ref objZFactor);
}
else
{
pGeoDs = pRasterSurfaceOp.Slope(pInRaster as IGeoDataset, esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopePercentrise, ref objZFactor);
}
IRaster pOutRaster = (IRaster)pGeoDs;
pOutRasterLayer.CreateFromRaster(pOutRaster);
pOutRasterLayer.Name = strOutFileName;
IWorkspaceFactory pWSF = new RasterWorkspaceFactoryClass();
IWorkspace pRWS = pWSF.OpenFromFile(strOutDir, 0);
IRasterBandCollection pRasBandCol = (IRasterBandCollection)pGeoDs;
pRasBandCol.SaveAs(strOutFileName, pRWS, "GRID");//"IMAGINE Image"
pOutRasterLayer = UtilityFunction.SetStretchRenderer(pOutRasterLayer.Raster);
pOutRasterLayer.Name = strOutFileName;
pMap.AddLayer(pOutRasterLayer);
}
}
catch (Exception ex)
{
MessageBox.Show(ex.Message);
}
}
private IRasterLayer calTerrain(string _outPath)
{
//绑定栅格到 pMapAlgebraOp
try
{
if (pMapAlgebraOp == null)
{
pMapAlgebraOp = new RasterMapAlgebraOpClass();
}
string strExp = "";
//dem
IRaster pRasterDem = (pMap.get_Layer(LayerOprate.getLayerIndexByName(pMap, this.cmbDem.Text)) as IRasterLayer).Raster;
IGeoDataset pGeoDsDem = pRasterDem as IGeoDataset;
// bind raster data "Dem"
pMapAlgebraOp.BindRaster(pGeoDsDem, "Dem");
//cal Slope with Deg
ISurfaceOp pSurfaceOp = new RasterSurfaceOpClass();
object objZFactor = System.Type.Missing;
IGeoDataset pGeoDsSlope = pSurfaceOp.Slope(pRasterDem as IGeoDataset, esriGeoAnalysisSlopeEnum.esriGeoAnalysisSlopeDegrees, ref objZFactor);
// bind raster data "Slope"
pMapAlgebraOp.BindRaster(pGeoDsSlope, "Slope");
//cal Slo dem
string demStart = this.txtDemStart.Text;
string demEnd = this.txtDemEnd.Text;
string sloStart = this.txtSloStart.Text;
string sloEnd = this.txtSloEnd.Text;
//strExp = "Con([Dem] < " + demStart + ",0,[Dem] > " + demEnd + ",0,1)";
//IGeoDataset pGeoDsDemC = pMapAlgebraOp.Execute(strExp);
//pMapAlgebraOp.BindRaster(pGeoDsDemC, "DemC");
////限制区为1，非限制区为0
//strExp = "Con([Slope] < " + sloStart + ",0,[Slope] > "+sloEnd+",0,1)";
//IGeoDataset pGeoDsSloC = pMapAlgebraOp.Execute(strExp);
//pMapAlgebraOp.BindRaster(pGeoDsSloC, "SloC");
//cal in together
strExp = "Con((([Slope] > " + sloStart + " & [Slope] < " + sloEnd + ") & ([Dem] > " + demStart + " & [Dem] < " + demEnd + ")),1,0)";
IGeoDataset pGeoDsC = pMapAlgebraOp.Execute(strExp);
pMapAlgebraOp.BindRaster(pGeoDsC, "TerrainRedL");
pMapAlgebraOp.UnbindRaster("Dem");
pMapAlgebraOp.UnbindRaster("Slope");
IRaster pOutRaster = pGeoDsC as IRaster;
IRasterLayer pOutRasterLayer = new RasterLayerClass();
pOutRasterLayer.CreateFromRaster(pOutRaster);
//数据名称
//string strOutDir = _outPath.Substring(0, _outPath.LastIndexOf("\\"));
int startX = _outPath.LastIndexOf("\\");
int endX = _outPath.Length;
string sRLyrName = _outPath.Substring(startX + 1, endX - startX - 1);
pOutRasterLayer.Name = sRLyrName;
pMapAlgebraOp = null;
return pOutRasterLayer;
}
catch (Exception ex)
{
MessageBox.Show(ex.Message.ToString());
return null;
}
}
