/// <summary>
        /// 根据图层的列表计算土地利用变化历史趋势。
        /// </summary>
        /// <param name="listLayers"></param>
        private void CalculateHistoryTrends(List <string> listLayers)
        {
            //写入列。列名为图层名。
            DataTable dtResults = new DataTable();

            for (int i = 0; i < listLayers.Count; i++)
            {
                dtResults.Columns.Add(listLayers[i]);
            }
            //写入行。除了所有的土地利用类型的栅格数量,还有数量总计,土地开发强度,城市用地增长率。
            for (int i = 0; i < landUseTypesCount + 3; i++)
            {
                DataRow row = dtResults.NewRow();
                dtResults.Rows.Add(row);
            }
            //计算并写入数据
            int          urbanIndex = GetUrbanIndex();
            int          urbanCount = 0;
            IRasterLayer rasterLayer;
            IRasterUniqueValueRenderer rasterUniqueValueRenderer;

            for (int i = 0; i < listLayers.Count; i++)
            {
                rasterLayer = ArcGISOperator.GetRasterLayerByName(listAllLayersInMapName[i]);
                rasterUniqueValueRenderer = (IRasterUniqueValueRenderer)rasterLayer.Renderer;
                int[] counts = ArcGISOperator.GetLandUseTypesCount(rasterLayer);
                int   sum    = 0;
                for (int j = 0; j < landUseTypesCount; j++)
                {
                    dtResults.Rows[j][i] = counts[j];
                    sum += counts[j];
                    if (urbanIndex == j)
                    {
                        urbanCount = counts[j];
                    }
                }
                dtResults.Rows[landUseTypesCount][i]     = sum;
                dtResults.Rows[landUseTypesCount + 1][i] = (Convert.ToDecimal(urbanCount) * 100 / sum).ToString("0.00");
                if (i != 0)
                {
                    int oldUrbanCount = Convert.ToInt32(dtResults.Rows[urbanIndex][i - 1]);
                    int newUrbanCount = Convert.ToInt32(dtResults.Rows[urbanIndex][i]);
                    dtResults.Rows[landUseTypesCount + 2][i] = (Convert.ToDecimal(newUrbanCount - oldUrbanCount) * 100
                                                                / oldUrbanCount).ToString("0.00");
                }
            }
            dataGridViewResult.DataSource = dtResults;
            //设置行名
            rasterLayer = ArcGISOperator.GetRasterLayerByName(listAllLayersInMapName[0]);
            rasterUniqueValueRenderer = (IRasterUniqueValueRenderer)rasterLayer.Renderer;
            for (int i = 0; i < landUseTypesCount; i++)
            {
                dataGridViewResult.RowHeadersWidth          = 130;
                dataGridViewResult.Rows[i].HeaderCell.Value = rasterUniqueValueRenderer.get_Label(0, i);
            }
            dataGridViewResult.Rows[landUseTypesCount].HeaderCell.Value     = resourceManager.GetString("String150");
            dataGridViewResult.Rows[landUseTypesCount + 1].HeaderCell.Value = resourceManager.GetString("String151");
            dataGridViewResult.Rows[landUseTypesCount + 2].HeaderCell.Value = resourceManager.GetString("String152");
        }
Пример #2
0
 private void buttonCalConvertCells_Click(object sender, EventArgs e)
 {
     if (SimulationEndImageName != "")
     {
         IRasterLayer rasterLayer      = ArcGISOperator.GetRasterLayerByName(simulationStartImageName);
         int          startUrbanConuts = ArcGISOperator.GetUrbanCount(rasterLayer, landUseClassificationInfo);
         rasterLayer = ArcGISOperator.GetRasterLayerByName(simulationEndImageName);
         int endUrbanCounts = ArcGISOperator.GetUrbanCount(rasterLayer, landUseClassificationInfo);
         numericUpDownConvertCount.Value = Convert.ToDecimal(endUrbanCounts - startUrbanConuts);
     }
 }
Пример #3
0
        void t_Tick(object sender, EventArgs e)
        {
            try
            {
                if (VariableMaintainer.IsNeedRefresh)
                {
                    string       simulationLayerName  = VariableMaintainer.CurrentFoucsMap.get_Layer(0).Name;
                    IRasterLayer simulationImageLayer = new RasterLayerClass();
                    if (VariableMaintainer.CurrentModel == EnumCurrentModel.Simulation_CA_LogisticRegression)
                    {
                        simulationImageLayer.CreateFromFilePath(VariableMaintainer.CurrentFormLogisticCAWizard.OutputFolder + @"\" + simulationLayerName);
                        ArcGISOperator.WriteRaster(simulationImageLayer, VariableMaintainer.CurrentFormLogisticCAWizard.SimulationImage);
                    }
                    else if (VariableMaintainer.CurrentModel == EnumCurrentModel.Simulation_CA_ANN)
                    {
                        simulationImageLayer.CreateFromFilePath(VariableMaintainer.CurrentFormANNCAWizard.OutputFolder + @"\" + simulationLayerName);
                        ArcGISOperator.WriteRaster(simulationImageLayer, VariableMaintainer.CurrentFormANNCAWizard.SimulationImage);
                    }
                    else if (VariableMaintainer.CurrentModel == EnumCurrentModel.Simulation_CA_DT)
                    {
                        simulationImageLayer.CreateFromFilePath(VariableMaintainer.CurrentFormDTCAWizard.OutputFolder + @"\" + simulationLayerName);
                        ArcGISOperator.WriteRaster(simulationImageLayer, VariableMaintainer.CurrentFormDTCAWizard.SimulationImage);
                    }

                    IRasterLayer l       = ArcGISOperator.GetRasterLayerByName(simulationLayerName);
                    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);
                    VariableMaintainer.IsNeedRefresh = false;
                }

                if (VariableMaintainer.IsSimulationFinished)
                {
                    VariableMaintainer.CurrentTimer.Stop();
                    VariableMaintainer.IsSimulationFinished = false;
                }
            }
            catch (Exception ex)
            {
                //MessageBox.Show(ex.Message);
            }
        }
Пример #4
0
        /// <summary>
        /// 读取并填充土地利用类型信息。
        /// </summary>
        private void FillLanduseInfo(string layerName)
        {
            dataGridViewLandUse.Rows.Clear();
            //根据第一个图层的符号化设置读取土地利用类型信息。
            //注意:因为一般来说符号化设置的设置和属性表一致,因为不需要考虑顺序不一致的问题。
            IRasterLayer rasterLayer = ArcGISOperator.GetRasterLayerByName(layerName);
            IRasterUniqueValueRenderer  rasterUniqueValueRenderer  = (IRasterUniqueValueRenderer)rasterLayer.Renderer;
            IRasterRendererUniqueValues rasterRendererUniqueValues = (IRasterRendererUniqueValues)rasterUniqueValueRenderer;
            IUniqueValues uniqueValues     = rasterRendererUniqueValues.UniqueValues;
            List <object> listUniqueValues = new List <object>();

            for (int i = 0; i < uniqueValues.Count; i++)
            {
                listUniqueValues.Add(uniqueValues.get_UniqueValue(i));
            }
            int classCount = rasterUniqueValueRenderer.get_ClassCount(0);

            for (int i = 0; i < classCount; i++)
            {
                dataGridViewLandUse.Rows.Add();
                dataGridViewLandUse.Rows[i].Cells[0].Value = rasterUniqueValueRenderer.get_Label(0, i);
            }
        }
Пример #5
0
        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);
            }
        }
Пример #6
0
        /// <summary>
        /// 更新信息素,同时更新最优状态
        /// </summary>
        /// <param name="Ants"></param>
        private void UpdatePheromoneData(Ant[] Ants)
        {
            for (int i = 0; i < structRasterMetaData.RowCount; i++)
            {
                for (int j = 0; j < structRasterMetaData.ColumnCount; j++)
                {
                    if (suitabilityData[i, j] != structRasterMetaData.NoDataValue)
                    {
                        //计算剩下的信息素
                        pheromoneData[i, j] = (float)(pheromoneData[i, j] * (1 - rho));
                        if (goalUtilityInMicroSearch > goalUtilityInAllIterations)
                        {
                            goalUtilityAntsStatusInAllIterations[i, j] = 0;
                        }
                    }
                }
            }

            double onePheromon = q * goalUtilityInMicroSearch;

            int conNN = 5;

            if (currentIteration > 200 && currentIteration <= 700)
            {
                conNN = 5;                                                      /*Beita=3;*/
            }
            if (currentIteration > 700)
            {
                conNN = 3;                           /*Beita=3;*/
            }

            for (int i = 0; i < Ants.Length; i++)
            {
                int x = Ants[i].X;
                int y = Ants[i].Y;

                //-------------计算增加的信息素-----------------------------------------------
                pheromoneData[x, y] += (float)onePheromon;                                                                                 //更新信息素
                UpdateNeighbour(pheromoneData, conNN, onePheromon, structRasterMetaData.ColumnCount, structRasterMetaData.RowCount, y, x); //更新邻域信息素;
                //-------------------------------------------------------------------------

                //判断是否需要更新最优状态
                if (goalUtilityInMicroSearch > goalUtilityInAllIterations)
                {
                    goalUtilityAntsStatusInAllIterations[x, y] = 1;
                }
            }

            if (goalUtilityInMicroSearch > goalUtilityInAllIterations)
            {
                goalUtilityInAllIterations = goalUtilityInMicroSearch;
                isNeedUpdate             = true;
                needFinishIterationCount = 1;
                currentRefreshIteration += 1;
            }
            else
            {
                needFinishIterationCount += 1;
            }

            currentIteration += 1;
            dockableWindowOutput.AppendText("CurrentInteration=" + currentIteration + "     CurrentRefreshIteration=" + currentRefreshIteration + "     NeedUpdate=" + isNeedUpdate + "     NeedFinishIterationCount=" + needFinishIterationCount + "       GoalUtilityValueInAllIterations=" + goalUtilityInAllIterations + "\r\n");
            dockableWindowOutput.AppendText(("\n"));
            dockableWindowOutput.ScrollTextbox();

            if (currentIteration == 1)
            {
                dockableWindowGraphy.CreateGraph(currentIteration, goalUtilityInAllIterations, "");
                dockableWindowGraphy.RefreshGraph();
            }
            else
            {
                dockableWindowGraphy.UpdateData(currentIteration, goalUtilityInAllIterations, 0);
                dockableWindowGraphy.RefreshGraph();
            }

            float[,] currentOccpied = 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] == -9999f)
                    {
                        currentOccpied[i, j] = 255;
                    }
                    else
                    {
                        if (occupiedStatus[i, j] == 1)
                        {
                            currentOccpied[i, j] = Convert.ToSingle(occupiedStatus[i, j]);
                        }
                        else
                        {
                            currentOccpied[i, j] = 0;
                        }
                    }
                }
            }

            int mapRefreshInterval = VariableMaintainer.CurrentStructACOParameters.MapRefreshInterval;

            if (mapRefreshInterval == 0)
            {
                mapRefreshInterval = 5;
            }
            if (currentIteration % mapRefreshInterval == 0)
            {
                string       simulationLayerName  = VariableMaintainer.CurrentFoucsMap.get_Layer(0).Name;
                IRasterLayer simulationImageLayer = new RasterLayerClass();
                simulationImageLayer.CreateFromFilePath(VariableMaintainer.DefaultOutputFolder + @"\" + simulationLayerName);
                ArcGISOperator.WriteRaster(simulationImageLayer, currentOccpied);
                IRasterLayer l = ArcGISOperator.GetRasterLayerByName(simulationLayerName);
                ((IRasterEdit)l.Raster).Refresh();

                IActiveView activeView = VariableMaintainer.CurrentFoucsMap as IActiveView;
                activeView.PartialRefresh(esriViewDrawPhase.esriViewGeography, null, null);
            }

            //保存每次迭代的目标值,如果需要保存
            if (isNeedSave == true)
            {
                goalUtilityArray[currentIteration + 1] = currentIteration + "," + currentRefreshIteration + "," + isNeedUpdate + "," + needFinishIterationCount + "," + goalUtilityInAllIterations;
                if (currentIteration == totalItearationCount)
                {
                    string outResult = outputPath + outGoalName;
                    //CommonOperator.writeGoalData(outResult, goalSave);
                }
            }

            //判断是否需要把结果保存成文件
            if (isNeedSave == true && needFinishIterationCount == 1)
            {
                outFileName = outPreFileName + currentRefreshIteration + "_s" + currentIteration + ".txt";
                String outResult = outputPath + outFileName;

                dockableWindowOutput.AppendText("Writing file:" + outResult + "\r\n");
                dockableWindowOutput.ScrollTextbox();

                //CommonOperator.WriteData(outResult, allOptStatus, rows, cols, xllcorner, yllcorner, cellsize, NODATA_value);

                dockableWindowOutput.AppendText("Writing file finished\r\n");
                dockableWindowOutput.ScrollTextbox();
            }
        }
Пример #7
0
        //--------------------------------------------------------------------------------------------------------------

        //------------------------------------------------方法----------------------------------------------------------

        /// <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";
        }
Пример #8
0
        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);
        }
Пример #9
0
        /// <summary>
        /// 当向导页的新一页出现前可以做的该页初始化工作。
        /// </summary>
        /// <param name="sender"></param>
        /// <param name="e"></param>
        private void wizardLogistcCA_AfterSwitchPages(object sender, CristiPotlog.Controls.Wizard.AfterSwitchPagesEventArgs e)
        {
            //切换至逻辑回归计算页时,应首先读取地图图层信息,填写部分控件
            if (e.OldIndex == 1 && e.NewIndex == 2)
            {
                //如果不是第一次打开窗体,则再次读取地图数据
                application             = ArcMap.Application;
                document                = application.Document as IMxDocument;
                map                     = document.FocusMap;
                ArcGISOperator.FoucsMap = document.FocusMap;
                //如果图层选择框为空
                if (comboBoxTrainingStartImage.Items.Count == 0)
                {
                    for (int i = 0; i < map.LayerCount; i++)
                    {
                        string layerName = map.get_Layer(i).Name;
                        comboBoxTrainingStartImage.Items.Add(layerName);
                        comboBoxTrainingEndImage.Items.Add(layerName);
                        comboBoxSimStartImage.Items.Add(layerName);
                        comboBoxSimEndImage.Items.Add(layerName);
                    }
                }
            }
            //切换到土地利用类型设置页时首先获取图层信息,然后填充土地利用类型
            else if (e.OldIndex == 2 && e.NewIndex == 3)
            {
                listVariableLayersName.Clear();
                trainingStartImageName = comboBoxTrainingStartImage.SelectedItem.ToString();
                trainingEndImageName   = comboBoxTrainingEndImage.SelectedItem.ToString();
                for (int i = 0; i < dataGridViewVariableDatas.Rows.Count; i++)
                {
                    listVariableLayersName.Add(dataGridViewVariableDatas.Rows[i].Cells[0].Value.ToString());
                }

                //如果已经完成了回归,则不再进行
                if (isRegressed)
                {
                    return;
                }

                dataGridViewLandUse.Rows.Clear();
                IRasterLayer rasterLayer = ArcGISOperator.GetRasterLayerByName(trainingStartImageName);
                IRasterUniqueValueRenderer  rasterUniqueValueRenderer  = (IRasterUniqueValueRenderer)rasterLayer.Renderer;
                IRasterRendererUniqueValues rasterRendererUniqueValues = (IRasterRendererUniqueValues)rasterUniqueValueRenderer;
                IUniqueValues uniqueValues     = rasterRendererUniqueValues.UniqueValues;
                List <object> listUniqueValues = new List <object>();
                for (int i = 0; i < uniqueValues.Count; i++)
                {
                    listUniqueValues.Add(uniqueValues.get_UniqueValue(i));
                }
                int classCount = rasterUniqueValueRenderer.get_ClassCount(0);
                for (int i = 0; i < classCount; i++)
                {
                    dataGridViewLandUse.Rows.Add();
                    dataGridViewLandUse.Rows[i].Cells[0].Value = listUniqueValues[i];
                    dataGridViewLandUse.Rows[i].Cells[1].Value = rasterUniqueValueRenderer.get_Label(0, i);

                    ISymbol     symbol     = rasterUniqueValueRenderer.get_Symbol(0, i);
                    IFillSymbol fillSymbol = symbol as IFillSymbol;
                    //ISimpleFillSymbol simpleFillSymbol = (ISimpleFillSymbol)symbol;
                    //IColor esriColor = simpleFillSymbol.Color;
                    IColor    esriColor = fillSymbol.Color;
                    IRgbColor rgbColor  = new RgbColorClass();
                    rgbColor.CMYK = esriColor.CMYK;
                    Color color = Color.FromArgb(rgbColor.Red, rgbColor.Green, rgbColor.Blue);
                    dataGridViewLandUse.Rows[i].Cells[3].Value =
                        GeneralOpertor.GetBitmap(15, 15, color);
                }
                ColumnLUCharacter.Items.Clear();
                ColumnLUCharacter.Items.Add(resourceManager.GetString("String6"));
                ColumnLUCharacter.Items.Add(resourceManager.GetString("String7"));
                ColumnLUCharacter.Items.Add(resourceManager.GetString("String8"));
            }
            //切换到逻辑回归计算页面时进行计算
            if (e.OldIndex == 3 && e.NewIndex == 4)
            {
                if (isRegressed)
                {
                    return;
                }

                dataGridViewLogisticRegressionCoeff.Rows.Clear();
                labelCalculte.Text           = resourceManager.GetString("String9");
                labelCalculte.Visible        = true;
                progressBarCalculate.Visible = true;
                Application.DoEvents();

                //1.读取起始图像、终止图像、各变量图像为数组
                IRasterLayer   rasterLayerStartImage = ArcGISOperator.GetRasterLayerByName(comboBoxTrainingStartImage.SelectedItem.ToString());
                IRasterLayer   rasterLayerEndImage   = ArcGISOperator.GetRasterLayerByName(comboBoxTrainingEndImage.SelectedItem.ToString());
                RasterSampling rasterSampling        = new RasterSampling();
                labelCalculte.Text = resourceManager.GetString("String12");
                Application.DoEvents();
                List <int> notNullRows;
                List <int> notNullColumns;
                trainingStartImage = ArcGISOperator.ReadRasterAndGetNotNullRowColumn(rasterLayerStartImage, out structRasterMetaData,
                                                                                     out notNullRows, out notNullColumns);
                //获取最小空间范围
                List <IRasterLayer> listVariablesLayers = new List <IRasterLayer>();
                for (int i = 0; i < listVariableLayersName.Count; i++)
                {
                    listVariablesLayers.Add(ArcGISOperator.GetRasterLayerByName(listVariableLayersName[i]));
                }
                ArcGISOperator.GetSmallestBound(rasterLayerStartImage, rasterLayerEndImage, listVariablesLayers, ref structRasterMetaData);

                labelCalculte.Text = resourceManager.GetString("String13");
                Application.DoEvents();
                trainingEndImage     = ArcGISOperator.ReadRaster(rasterLayerEndImage, structRasterMetaData.NoDataValue);
                float[,] changeImage = GeneralOpertor.GetBinaryImageByTwoImages(trainingStartImage, trainingEndImage, structRasterMetaData.RowCount,
                                                                                structRasterMetaData.ColumnCount, structRasterMetaData.NoDataValue);
                //GeneralOpertor.WriteDataFloat(GetOutputFolder() + @"\lr_change.txt", changeImage, structRasterMetaData.RowCount, structRasterMetaData.ColumnCount);

                foreach (string layerName in listVariableLayersName)
                {
                    labelCalculte.Text = resourceManager.GetString("String14") + layerName + ".....";
                    Application.DoEvents();
                    varibaleImages.Add(ArcGISOperator.ReadRaster(ArcGISOperator.GetRasterLayerByName(layerName), structRasterMetaData.NoDataValue));
                }
                //2.按设置的比例进行随机抽样,得到抽样结果,然后计算
                //1.进行抽样
                labelCalculte.Text = resourceManager.GetString("String11");
                Application.DoEvents();
                int samplingCellsCount = (int)(structRasterMetaData.RowCount * structRasterMetaData.ColumnCount * numericUpDownSamplingPrecent.Value / 100);
                float[,] datas = new float[samplingCellsCount, listVariableLayersName.Count + 1];
                //抽取真值的比例默认设置为一半
                datas = rasterSampling.SamplingData(varibaleImages, changeImage, (double)numericUpDownSamplingPrecent.Value,
                                                    50, structRasterMetaData, notNullRows, notNullColumns);
                //2.进行计算
                labelCalculte.Text = resourceManager.GetString("String9");
                Application.DoEvents();
                //输出抽样结果数据
                //GeneralOpertor.WriteDataFloat(GetOutputFolder() + @"\lrSamplingData.txt", datas,
                //    rasterSampling.SamplingCellsCount, listVariableLayersName.Count + 1);

                //3.使用Accord.NET进行逻辑回归计算
                labelCalculte.Text = resourceManager.GetString("String9");
                Application.DoEvents();
                LogisticRegression lr = new LogisticRegression();
                double[]           odds;
                double[]           stde;
                double[]           min;
                double[]           max;
                lr.LogisticRegresssionAanlysis(datas, rasterSampling.SamplingCellsCount, listVariableLayersName.Count + 1,
                                               out coef, out odds, out stde, out min, out max);

                for (int i = 0; i < listVariableLayersName.Count + 1; i++)
                {
                    dataGridViewLogisticRegressionCoeff.Rows.Add();
                    if (i == 0)
                    {
                        dataGridViewLogisticRegressionCoeff.Rows[i].Cells[0].Value = resourceManager.GetString("String15");
                    }
                    else
                    {
                        dataGridViewLogisticRegressionCoeff.Rows[i].Cells[0].Value = listVariableLayersName[i - 1];
                    }
                    dataGridViewLogisticRegressionCoeff.Rows[i].Cells[1].Value = coef[i];
                    dataGridViewLogisticRegressionCoeff.Rows[i].Cells[2].Value = odds[i];
                    dataGridViewLogisticRegressionCoeff.Rows[i].Cells[3].Value = stde[i];
                    dataGridViewLogisticRegressionCoeff.Rows[i].Cells[4].Value = max[i];
                    dataGridViewLogisticRegressionCoeff.Rows[i].Cells[5].Value = min[i];
                }

                labelCalculte.Visible        = false;
                progressBarCalculate.Visible = false;
                Application.DoEvents();
                isRegressed = true;
            }
            //切换到模拟参数设置页时给定控件初始值,同时保存模拟数据信息
            else if (e.OldIndex == 4 && e.NewIndex == 5)
            {
                if (comboBoxSimEndImage.SelectedIndex == -1)
                {
                    buttonCalConvertCells.Enabled = false;
                }
                else
                {
                    buttonCalConvertCells.Enabled = true;
                }
            }
            //切换到模拟过程参数设置页时给定控件初始值,同时保存模拟参数信息
            else if (e.OldIndex == 5 && e.NewIndex == 6)
            {
                simulationStartImageName = comboBoxSimStartImage.SelectedItem.ToString();
                if (comboBoxSimEndImage.SelectedIndex != -1)
                {
                    simulationEndImageName = comboBoxSimEndImage.SelectedItem.ToString();
                }

                convertCount         = Convert.ToInt32(numericUpDownConvertCount.Value);
                simulationIterations = Convert.ToInt32(numericUpDownIterations.Value);
                delta = Convert.ToInt32(numericUpDownDelta.Value);

                numericUpDownRefresh.Value     = 10;
                numericUpDownOutputImage.Value = 10;
                outputFolder             = GetOutputFolder();
                textBoxOutputFolder.Text = outputFolder;
            }
            //切换到完成页时填充摘要,同时保存模拟过程输出参数的信息
            else if (e.OldIndex == 6 && e.NewIndex == 7)
            {
                refreshInterval     = Convert.ToInt32(numericUpDownRefresh.Value);
                outputImageInterval = Convert.ToInt32(numericUpDownOutputImage.Value);
                if (radioButtonOutput.Checked)
                {
                    isOutput = true;
                }
                else
                {
                    isOutput = false;
                }
                outputFolder = textBoxOutputFolder.Text;

                textBoxSummay.Text = WriteSummay();
            }
        }