//--------------------------------------------------------------------------------------------------------------
//------------------------------------------------方法----------------------------------------------------------
/// <summary>
/// 方法
/// </summary>
#region 方法
/// <summary>
/// 初始化方法
/// </summary>
#region ----------------初始化方法------------------------
public void Initialize(StructACOParameters structACOParameters, int microSearchIterationCountValue,
Boolean isNeedSaveValue, DockableWindowOutput currentDockableWindowOutput,
DockableWindowGraphy currentDockableWindowGraphy, Boolean isUseCompactnessValue)
{
q = structACOParameters.Q; //信息素强度
rho = structACOParameters.Rho; //挥发因子
alpha = structACOParameters.Alpha; //信息素权重
beta = structACOParameters.Beta; //启发函数权重
weightSutiable = structACOParameters.WeightSuitable; //阻抗性权重
weightCompact = structACOParameters.WeightCompact; //紧凑性权重
totalItearationCount = structACOParameters.InterationCount; //总迭代次数
currentIteration = 0; //当前迭代次数
currentRefreshIteration = 0; //更新次数
needFinishIterationCount = 0;
microSearchIterationCount = microSearchIterationCountValue; //小搜索次数
goalUtilityAntsStatusInMicroSearch = new Ant[structACOParameters.AntCount]; //小搜索内最优状态
currentIterationInMicroSearch = 0;
goalUtilityInMicroSearch = 0;
isNeedSave = isNeedSaveValue; //判断是否需要保存结果
isNeedUpdate = false; //初始化判断是否需要更新最优状态
isInitialStatusSave = false; //判断是否已经保存初始状态数据
//rows = rrows; //网格行数
//cols = ccols; //网格列数
//xllcorner = xxllcorner; //输出文件左上角横坐标
//yllcorner = yyllcorner; //输出文件左上角纵坐标
//cellsize = ccellsize; //栅格大小
//noDataValue = nnodata_value; //无数据值
outputPath = structACOParameters.OutputFolder; //结果输出路径
outInitName = "_Init.txt"; //结果输出初始状态文件名
outFileName = ""; //结果输出文件名
outPreFileName = "g"; //结果输出文件名前面部分
outGoalName = "GoalFunction.txt"; //结果输出目标函数名
antsCount = structACOParameters.AntCount; //智能体个数
//cellSize = bbianChang; //网格边长
//cellArea = cellSize * cellSize * antsCount; //智能体面积
goalUtilityInAllIterations = 0;
randomBound = 1;
goalUtilityArray = new string[totalItearationCount + 2]; //保存目标函数值数据
suitabilityData = GetSuitabilityData(VariableMaintainer.OptimizationExpression);
InitializeOccupiedAndGoalUtilityAntsStatusInAllIterations(); //初始化挖空
InitializeAnts(); //初始化智能体
InitializePheromone(); //初始化信息素
InitializeHeuristicFunction(); //初始化启发函数
string simulationLayerName = VariableMaintainer.CurrentFoucsMap.get_Layer(0).Name;
IRasterLayer simulationImageLayer = ArcGISOperator.GetRasterLayerByName(simulationLayerName);
string dateTime = GeneralOpertor.GetDataTimeFullString(DateTime.Now);
string rasterName = "ACOArea" + dateTime + ".img";
IRasterDataset rst = ArcGISOperator.CreateRasterDataset(VariableMaintainer.DefaultOutputFolder,
rasterName, simulationImageLayer, structRasterMetaData, occupiedStatus, 0);
IRasterLayer n = new RasterLayerClass();
n.CreateFromDataset(rst);
IColorRamp colorRamp = new RandomColorRampClass();
ArcGISOperator.UniqueValueRenderer(colorRamp, n, "Value");
ArcGISOperator.FoucsMap.AddLayer((ILayer)n);
double pixelWidth = ((IRasterDefaultProps)n.Raster).DefaultPixelWidth;
double pixelHeight = ((IRasterDefaultProps)n.Raster).DefaultPixelHeight;
cellArea = pixelWidth * pixelHeight * antsCount;
cellSize = pixelWidth;
random = new Random();
dockableWindowOutput = currentDockableWindowOutput;
dockableWindowGraphy = currentDockableWindowGraphy;
stringBuilderMessage = new StringBuilder();
isUseCompactness = isUseCompactnessValue;
dockableWindowGraphy.GraphTitle = "Goal Utility";
dockableWindowGraphy.XAxisTitle = "Iteration";
dockableWindowGraphy.YAxisTitle = "Goal Utility Value";
}
protected override void OnClick()
{
try
{
ResourceManager resourceManager = VariableMaintainer.CurrentResourceManager;
if (VariableMaintainer.CurrentModel == EnumCurrentModel.Null)
{
MessageBox.Show(resourceManager.GetString("String168"), resourceManager.GetString("String2"),
MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
return;
}
Thread thread = VariableMaintainer.SimulationThread;
if (thread != null)
{
if (thread.ThreadState != ThreadState.Stopped)
{
MessageBox.Show(resourceManager.GetString("String39"), resourceManager.GetString("String2"),
MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
return;
}
}
StructRasterMetaData structRasterMetaData = new StructRasterMetaData();
string model = "";
if (VariableMaintainer.CurrentModel == EnumCurrentModel.Simulation_CA_LogisticRegression)
{
model = resourceManager.GetString("String104");
}
else if (VariableMaintainer.CurrentModel == EnumCurrentModel.Simulation_CA_ANN)
{
model = resourceManager.GetString("String102");
}
else if (VariableMaintainer.CurrentModel == EnumCurrentModel.Simulation_CA_DT)
{
model = resourceManager.GetString("String103");
}
else if (VariableMaintainer.CurrentModel == EnumCurrentModel.Optimization_Area)
{
model = resourceManager.GetString("String167");
}
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String82")
+ model + resourceManager.GetString("String101") + DateTime.Now.ToShortTimeString());
VariableMaintainer.CurrentDockableWindowOutput.AppendText("\n");
VariableMaintainer.CurrentDockableWindowOutput.AppendText("\n");
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String58"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String58"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText("\n");
if (VariableMaintainer.CurrentModel == EnumCurrentModel.Simulation_CA_LogisticRegression)
{
//ArcGIS操作在另一线程中执行效率降低很多，因此把数据读取及刷新操作在线程外执行
//首先读取模拟起始时刻和/或终止时刻影像
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String12"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String12") +
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationStartImageName + ".....\n");
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationStartImage =
ArcGISOperator.ReadRasterAndGetMetaData(ArcGISOperator.GetRasterLayerByName(
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationStartImageName), out structRasterMetaData);
if (VariableMaintainer.CurrentFormLogisticCAWizard.SimulationEndImageName != "")
{
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String13"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String13") +
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationEndImageName + ".....\n");
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationEndImage = ArcGISOperator.ReadRaster(ArcGISOperator.GetRasterLayerByName(
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationEndImageName), -9999f);
}
//然后读取各变量影像
if (VariableMaintainer.CurrentFormLogisticCAWizard.IsUsingDefault)
{
foreach (string layerName in VariableMaintainer.CurrentFormLogisticCAWizard.ListVariableLayersName)
{
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String14") + layerName + ".....");
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String14") + layerName + ".....\n");
VariableMaintainer.CurrentFormLogisticCAWizard.VaribaleImages.Add(
ArcGISOperator.ReadRaster(ArcGISOperator.GetRasterLayerByName(layerName), -9999f));
}
//获取最小空间范围
IRasterLayer rasterLayerStartImage = ArcGISOperator.GetRasterLayerByName(VariableMaintainer.CurrentFormLogisticCAWizard.SimulationStartImageName);
IRasterLayer rasterLayerEndImage = ArcGISOperator.GetRasterLayerByName(VariableMaintainer.CurrentFormLogisticCAWizard.SimulationEndImageName);
List <IRasterLayer> listVariablesLayers = new List <IRasterLayer>();
for (int i = 0; i < VariableMaintainer.CurrentFormLogisticCAWizard.ListVariableLayersName.Count; i++)
{
listVariablesLayers.Add(ArcGISOperator.GetRasterLayerByName(VariableMaintainer.CurrentFormLogisticCAWizard.ListVariableLayersName[i]));
}
ArcGISOperator.GetSmallestBound(rasterLayerStartImage, rasterLayerEndImage, listVariablesLayers, ref structRasterMetaData);
VariableMaintainer.CurrentFormLogisticCAWizard.CurrentStructRasterMetaData = structRasterMetaData;
}
//最后用模拟起始时刻影像创建模拟影像
string dateTime = GeneralOpertor.GetDataTimeFullString(DateTime.Now);
string rasterName = "sim" + dateTime + ".img";
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String37"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String37") + rasterName + ".....\n");
IRasterLayer simulationStartImageLayer = ArcGISOperator.GetRasterLayerByName(
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationStartImageName);
IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass();
IWorkspace workspace = workspaceFactory.OpenFromFile(VariableMaintainer.CurrentFormLogisticCAWizard.OutputFolder, 0);
ISaveAs2 saveAs2 = (ISaveAs2)simulationStartImageLayer.Raster;
saveAs2.SaveAs(rasterName, workspace, "IMAGINE Image");
IRasterLayer simulationImageLayer = new RasterLayerClass();
simulationImageLayer.CreateFromFilePath(VariableMaintainer.CurrentFormLogisticCAWizard.OutputFolder + @"\" + rasterName);
simulationImageLayer.Renderer = simulationStartImageLayer.Renderer;
ArcGISOperator.FoucsMap.AddLayer((ILayer)simulationImageLayer);
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationImage = ArcGISOperator.ReadRaster(simulationImageLayer, -9999f);
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationLayerName = rasterName;
LogisticRegreesionCA lrCA = new LogisticRegreesionCA();
lrCA.DockableWindowGraphy = VariableMaintainer.CurrentDockableWindowGraphy;
lrCA.DockableWindowOutput = VariableMaintainer.CurrentDockableWindowOutput;
thread = new Thread(new ThreadStart(lrCA.DoSimulation));
}
else if (VariableMaintainer.CurrentModel == EnumCurrentModel.Simulation_CA_ANN)
{
//今后应使用父类和接口进行重构
//首先读取模拟起始时刻和/或终止时刻影像
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String12"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String12") +
VariableMaintainer.CurrentFormANNCAWizard.SimulationStartImageName + ".....\n");
VariableMaintainer.CurrentFormANNCAWizard.SimulationStartImage =
ArcGISOperator.ReadRasterAndGetMetaData(ArcGISOperator.GetRasterLayerByName(
VariableMaintainer.CurrentFormANNCAWizard.SimulationStartImageName), out structRasterMetaData);
//VariableMaintainer.CurrentANNCAWizard.CurrentStructRasterMetaData = structRasterMetaData;
if (VariableMaintainer.CurrentFormANNCAWizard.SimulationEndImageName != "")
{
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String13"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String13") +
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationEndImageName + ".....\n");
VariableMaintainer.CurrentFormANNCAWizard.SimulationEndImage = ArcGISOperator.ReadRaster(ArcGISOperator.GetRasterLayerByName(
VariableMaintainer.CurrentFormANNCAWizard.SimulationEndImageName), -9999f);
}
//最后用模拟起始时刻影像创建模拟影像
string dateTime = GeneralOpertor.GetDataTimeFullString(DateTime.Now);
string rasterName = "sim" + dateTime + ".img";
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String37"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String37") + rasterName + ".....\n");
IRasterLayer simulationStartImageLayer = ArcGISOperator.GetRasterLayerByName(
VariableMaintainer.CurrentFormANNCAWizard.SimulationStartImageName);
IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass();
IWorkspace workspace = workspaceFactory.OpenFromFile(VariableMaintainer.CurrentFormANNCAWizard.OutputFolder, 0);
ISaveAs2 saveAs2 = (ISaveAs2)simulationStartImageLayer.Raster;
saveAs2.SaveAs(rasterName, workspace, "IMAGINE Image");
IRasterLayer simulationImageLayer = new RasterLayerClass();
simulationImageLayer.CreateFromFilePath(VariableMaintainer.CurrentFormANNCAWizard.OutputFolder + @"\" + rasterName);
simulationImageLayer.Renderer = simulationStartImageLayer.Renderer;
ArcGISOperator.FoucsMap.AddLayer((ILayer)simulationImageLayer);
VariableMaintainer.CurrentFormANNCAWizard.SimulationImage = ArcGISOperator.ReadRaster(simulationImageLayer, -9999f);
VariableMaintainer.CurrentFormANNCAWizard.SimulationLayerName = rasterName;
ANNCA annCA = new ANNCA();
annCA.DockableWindowGraphy = VariableMaintainer.CurrentDockableWindowGraphy;
annCA.DockableWindowOutput = VariableMaintainer.CurrentDockableWindowOutput;
thread = new Thread(new ThreadStart(annCA.DoSimulation));
}
else if (VariableMaintainer.CurrentModel == EnumCurrentModel.Simulation_CA_DT)
{
//今后应使用父类和接口进行重构
//首先读取模拟起始时刻和/或终止时刻影像
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String12"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String12") +
VariableMaintainer.CurrentFormDTCAWizard.SimulationStartImageName + ".....\n");
VariableMaintainer.CurrentFormDTCAWizard.SimulationStartImage =
ArcGISOperator.ReadRasterAndGetMetaData(ArcGISOperator.GetRasterLayerByName(
VariableMaintainer.CurrentFormDTCAWizard.SimulationStartImageName), out structRasterMetaData);
//VariableMaintainer.CurrentDTCAWizard.CurrentStructRasterMetaData = structRasterMetaData;
if (VariableMaintainer.CurrentFormDTCAWizard.SimulationEndImageName != "")
{
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String13"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String13") +
VariableMaintainer.CurrentFormLogisticCAWizard.SimulationEndImageName + ".....\n");
VariableMaintainer.CurrentFormDTCAWizard.SimulationEndImage = ArcGISOperator.ReadRaster(ArcGISOperator.GetRasterLayerByName(
VariableMaintainer.CurrentFormDTCAWizard.SimulationEndImageName), -9999f);
}
//20170619添加限制层数据
//VariableMaintainer.CurrentDockableWindowOutput.AppendText("Add Restriction Layer " +
// "ACOres" + ".....\n");
//VariableMaintainer.RestrictImage = ArcGISOperator.ReadRaster(ArcGISOperator.GetRasterLayerByName(
// "ACOres"), -9999f);
//最后用模拟起始时刻影像创建模拟影像
string dateTime = GeneralOpertor.GetDataTimeFullString(DateTime.Now);
string rasterName = "sim" + dateTime + ".img";
ArcMap.Application.StatusBar.set_Message(0, resourceManager.GetString("String37"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(resourceManager.GetString("String37") + rasterName + ".....\n");
IRasterLayer simulationStartImageLayer = ArcGISOperator.GetRasterLayerByName(
VariableMaintainer.CurrentFormDTCAWizard.SimulationStartImageName);
IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass();
IWorkspace workspace = workspaceFactory.OpenFromFile(VariableMaintainer.CurrentFormDTCAWizard.OutputFolder, 0);
ISaveAs2 saveAs2 = (ISaveAs2)simulationStartImageLayer.Raster;
saveAs2.SaveAs(rasterName, workspace, "IMAGINE Image");
IRasterLayer simulationImageLayer = new RasterLayerClass();
simulationImageLayer.CreateFromFilePath(VariableMaintainer.CurrentFormDTCAWizard.OutputFolder + @"\" + rasterName);
simulationImageLayer.Renderer = simulationStartImageLayer.Renderer;
ArcGISOperator.FoucsMap.AddLayer((ILayer)simulationImageLayer);
VariableMaintainer.CurrentFormDTCAWizard.SimulationImage = ArcGISOperator.ReadRaster(simulationImageLayer, -9999f);
VariableMaintainer.CurrentFormDTCAWizard.SimulationLayerName = rasterName;
DecisionTreeCA dtCA = new DecisionTreeCA();
dtCA.DockableWindowGraphy = VariableMaintainer.CurrentDockableWindowGraphy;
dtCA.DockableWindowOutput = VariableMaintainer.CurrentDockableWindowOutput;
thread = new Thread(new ThreadStart(dtCA.DoSimulation));
}
else if (VariableMaintainer.CurrentModel == EnumCurrentModel.Optimization_Area)
{
if (!VariableMaintainer.IsACOUtilitySet || !VariableMaintainer.IsACOParametersSet)
{
MessageBox.Show(resourceManager.GetString("String166"), resourceManager.GetString("String2"),
MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
return;
}
AreaOptimizationACO areaOptimizationACO = new AreaOptimizationACO();
areaOptimizationACO.Initialize(VariableMaintainer.CurrentStructACOParameters, 1, false,
VariableMaintainer.CurrentDockableWindowOutput, VariableMaintainer.CurrentDockableWindowGraphy, true);
thread = new Thread(new ThreadStart(areaOptimizationACO.Run));
}
else
{
MessageBox.Show(resourceManager.GetString("String168"), resourceManager.GetString("String2"),
MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
return;
}
VariableMaintainer.SimulationThread = thread;
thread.SetApartmentState(ApartmentState.STA);
thread.Start();
VariableMaintainer.CurrentTimer.Interval = 500;
VariableMaintainer.CurrentTimer.Tick += new EventHandler(t_Tick);
VariableMaintainer.CurrentTimer.Enabled = true;
VariableMaintainer.CurrentTimer.Start();
ArcMap.Application.CurrentTool = null;
}
catch (Exception ex)
{
MessageBox.Show(ex.Message);
}
}
public float[,] GetSuitabilityData(string expression)
{
string[] splitedExpressionArray;
bool isSplited = SplitExpression(expression, out splitedExpressionArray);
if (isSplited == false)
{
MessageBox.Show("Suitibility Function Express Error", VariableMaintainer.CurrentResourceManager.GetString("String2"),
MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
return(null);
}
//根据获取的文件名字符串去读取相应的栅格文件
List <float[, ]> listLayers = new List <float[, ]>();
structRasterMetaData = new StructRasterMetaData();
for (int i = 0; i < splitedExpressionArray.Length; i++)
{
if (i > 0 && i % 2 > 0)
{
string layerName = splitedExpressionArray[i];
if (i == 1)
{
listLayers.Add(ArcGISOperator.ReadRasterAndGetMetaData(
ArcGISOperator.GetRasterLayerByName(layerName), out structRasterMetaData));
}
else
{
listLayers.Add(ArcGISOperator.ReadRaster(ArcGISOperator.GetRasterLayerByName(layerName),
structRasterMetaData.NoDataValue));
}
}
}
if (listLayers.Count == 0)
{
MessageBox.Show("Please select layers for calculation!", VariableMaintainer.CurrentResourceManager.GetString("String2"),
MessageBoxButtons.OK, MessageBoxIcon.Exclamation);
return(null);
}
//读取全部文件后进行计算，生成适宜性图层数据，并进行归一化
float[,] suitabilityData = new float[structRasterMetaData.RowCount, structRasterMetaData.ColumnCount];
float min = 1, max = 0;
Expression ex = new Expression("1+1");
//MathParser mathParser = new MathParser();
for (int i = 0; i < structRasterMetaData.RowCount; i++)
{
for (int j = 0; j < structRasterMetaData.ColumnCount; j++)
{
bool isNullData = false;
for (int k = 0; k < listLayers.Count; k++)
{
if (listLayers[k][i, j] == structRasterMetaData.NoDataValue)
{
isNullData = true;
break;
}
}
if (isNullData)
{
suitabilityData[i, j] = structRasterMetaData.NoDataValue;
}
else
{
string calculateExpression = "";
for (int w = 0; w < splitedExpressionArray.Length; w++)
{
if (w > 0 && w % 2 == 1)
{
calculateExpression += listLayers[w / 2][i, j];
}
else
{
calculateExpression += splitedExpressionArray[w];
}
}
ex = new Expression(calculateExpression);
suitabilityData[i, j] = Convert.ToSingle(ex.Evaluate());
//suitabilityData[i, j] = Convert.ToSingle(mathParser.Parse(calculateExpression, false));
if (suitabilityData[i, j] < min)
{
min = suitabilityData[i, j];
}
if (suitabilityData[i, j] > max)
{
max = suitabilityData[i, j];
}
}
}
}
float[,] suitabilityData2 = new float[structRasterMetaData.RowCount, structRasterMetaData.ColumnCount];
for (int i = 0; i < structRasterMetaData.RowCount; i++)
{
for (int j = 0; j < structRasterMetaData.ColumnCount; j++)
{
if (suitabilityData[i, j] == structRasterMetaData.NoDataValue)
{
suitabilityData2[i, j] = structRasterMetaData.NoDataValue;
}
else
{
suitabilityData2[i, j] = (suitabilityData[i, j] - min) / (max - min);
}
}
}
//将其他图层关闭，以免影响蚁群优化的刷新效果
for (int i = 0; i < VariableMaintainer.CurrentFoucsMap.LayerCount; i++)
{
if (VariableMaintainer.CurrentFoucsMap.get_Layer(i).Visible)
{
VariableMaintainer.CurrentFoucsMap.get_Layer(i).Visible = false;
}
}
//显示适宜性图层
string dateTime = GeneralOpertor.GetDataTimeFullString(DateTime.Now);
string rasterName = "suitability" + dateTime + ".img";
ArcMap.Application.StatusBar.set_Message(0, VariableMaintainer.CurrentResourceManager.GetString("String169"));
VariableMaintainer.CurrentDockableWindowOutput.AppendText(VariableMaintainer.CurrentResourceManager.GetString("String169") + rasterName + ".....\n");
VariableMaintainer.CurrentDockableWindowOutput.AppendText("\n");
IRasterLayer templateImageLayer = ArcGISOperator.GetRasterLayerByName(
splitedExpressionArray[1].ToString());
IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass();
IWorkspace workspace = workspaceFactory.OpenFromFile(VariableMaintainer.DefaultOutputFolder, 0);
ISaveAs2 saveAs2 = (ISaveAs2)templateImageLayer.Raster;
saveAs2.SaveAs(rasterName, workspace, "IMAGINE Image");
IRasterLayer sutiabilityImageLayer = new RasterLayerClass();
sutiabilityImageLayer.CreateFromFilePath(VariableMaintainer.DefaultOutputFolder + @"\" + rasterName);
//sutiabilityImageLayer.Renderer = templateImageLayer.Renderer;
ArcGISOperator.FoucsMap.AddLayer((ILayer)sutiabilityImageLayer);
IRasterLayer simulationImageLayer = new RasterLayerClass();
simulationImageLayer.CreateFromFilePath(VariableMaintainer.DefaultOutputFolder + @"\" + rasterName);
ArcGISOperator.WriteRaster(simulationImageLayer, suitabilityData2);
IRasterLayer l = ArcGISOperator.GetRasterLayerByName(rasterName);
IRaster2 raster = simulationImageLayer.Raster as IRaster2;
IPnt fromPnt = new PntClass();
fromPnt.SetCoords(0, 0);
IPnt blockSize = new PntClass();
blockSize.SetCoords(simulationImageLayer.ColumnCount, simulationImageLayer.RowCount);
IPixelBlock pixelBlock = ((IRaster)raster).CreatePixelBlock(blockSize);
l.Raster.Read(fromPnt, pixelBlock);
((IRasterEdit)l.Raster).Refresh();
IActiveView activeView = VariableMaintainer.CurrentFoucsMap as IActiveView;
activeView.PartialRefresh(esriViewDrawPhase.esriViewGeography, null, null);
return(suitabilityData2);
}
