public IRasterDataset3 getRasterDataset(string catalogpath)
{
IRasterDataset3 rasterDataset = null;
if (exists(catalogpath))
{
rasterDataset = (IRasterDataset3)gpUtl.OpenRasterDatasetFromString(catalogpath);
}
return(rasterDataset);
}
private void simpleButton4_Click(object sender, EventArgs e)
{
if (this.comboBoxInputRaster.Text.ToString() == "" || this.comboBoxContrastRaster.Text.ToString() == "" || this.textBoxOutputRaster.Text.ToString() == "")
{
MessageBox.Show("请输入参数！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}//有任意一个空为空则弹出对话框提示：请输入参数。
string strResultsDBPath = m_pGDBHelper.GetResultsDBPath();//声明变量获取数据库结果路径
string ToolName = "ToolName";
if (comboBoxMethod.Text == "主成分差异法")
{
ToolName = "NewPCADiff";
}
if (comboBoxMethod.Text == "多波段主成分变换法")
{
ToolName = "NewBandPCA";
}
if (comboBoxMethod.Text == "差异主成分法")
{
ToolName = "NewDiffPCA";
}
if (Utilities.GDBUtilites.CheckNameExist(strResultsDBPath, this.textBoxOutputRaster.Text.ToString()) == false)//如果没有同名文件
{
//执行PCA分析
ILayer pInputLayer1 = null;
ILayer pInputLayer2 = null;
for (int i = 0; i < m_pMapCtrl.LayerCount; i++)
{
if (m_pMapCtrl.get_Layer(i).Name == this.comboBoxInputRaster.Text.ToString())
{
pInputLayer1 = m_pMapCtrl.get_Layer(i);
}
else if (m_pMapCtrl.get_Layer(i).Name == this.comboBoxContrastRaster.Text.ToString())
{
pInputLayer2 = m_pMapCtrl.get_Layer(i);
}
}
if (pInputLayer1 != null && pInputLayer2 != null)
{
IRasterLayer pRasterLayer1 = pInputLayer1 as IRasterLayer;
IRasterLayer pRasterLayer2 = pInputLayer2 as IRasterLayer;
string rasterpath1 = pRasterLayer1.FilePath;
string rasterpath2 = pRasterLayer2.FilePath;
IGeoProcessor2 pGP = new GeoProcessorClass();
string GPPath = m_pGDBHelper.GetToolboxPath();
pGP.AddToolbox(GPPath);
IVariantArray gpParameters = new VarArrayClass();
//gpParameters.Add(pRasterLayer1.Raster);
//gpParameters.Add(pRasterLayer2.Raster);
gpParameters.Add(rasterpath1);
gpParameters.Add(rasterpath2);
gpParameters.Add(strResultsDBPath + "\\" + this.textBoxOutputRaster.Text);
IGeoProcessorResult pGeoProcessorResult = null;
try
{
pGeoProcessorResult = pGP.Execute(ToolName, gpParameters, null);
}
catch (Exception ex)
{
if (ex is System.Runtime.InteropServices.COMException)
{
int errorCode = (ex as System.Runtime.InteropServices.COMException).ErrorCode;
MessageBox.Show(errorCode.ToString());
}
}
if (pGeoProcessorResult.Status == esriJobStatus.esriJobSucceeded)
{
if (this.checkBoxAdd.Checked)
{
IWorkspaceFactory2 pWKF = new FileGDBWorkspaceFactoryClass();
IRasterWorkspaceEx pRasterWKEx = (IRasterWorkspaceEx)pWKF.OpenFromFile(m_pGDBHelper.GetResultsDBPath(), 0);
IRasterDataset3 pRasterDataset = pRasterWKEx.OpenRasterDataset(this.textBoxOutputRaster.Text) as IRasterDataset3;
Utilities.MapUtilites.AddRasterLayer(m_pMapCtrl.ActiveView, pRasterDataset, null);
}
MessageBox.Show("PCA分析完成！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
}
}
this.DialogResult = DialogResult.OK;
this.Close();
}
else
{
MessageBox.Show("已存在同名数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
}
public void DoSimulation()
{
try
{
//0准备开始模拟
ResourceManager resourceManager = VariableMaintainer.CurrentResourceManager;
FormANNCAWizard formANNCAWizard = VariableMaintainer.CurrentFormANNCAWizard;
Random random = new Random();
Stopwatch stopWatch = new Stopwatch();
stopWatch.Start();
//0.1获得数据
int rowCount = formANNCAWizard.CurrentStructRasterMetaData.RowCount;
int columnCount = formANNCAWizard.CurrentStructRasterMetaData.ColumnCount;
float[,] simulationStartImage = formANNCAWizard.SimulationStartImage;
float[,] simulationEndImage = formANNCAWizard.SimulationEndImage;
float[,] simulationImage = formANNCAWizard.SimulationImage;
//0.2计算初始每种土地利用类型的单元数量
List <StructLanduseInfoAndCount> listLanduseInfoAndCount; //记录每种土地利用类型的单元数
listLanduseInfoAndCount = new List <StructLanduseInfoAndCount>();
StructLanduseInfoAndCount landuseInfoAndCount;
foreach (StructLanduseInfo structLanduseInfo in formANNCAWizard.LandUseClassificationInfo.AllTypes)
{
landuseInfoAndCount = new StructLanduseInfoAndCount();
landuseInfoAndCount.structLanduseInfo = structLanduseInfo;
landuseInfoAndCount.LanduseTypeCount = 0;
listLanduseInfoAndCount.Add(landuseInfoAndCount);
}
for (int i = 0; i < rowCount; i++)
{
for (int j = 0; j < columnCount; j++)
{
CommonLibrary.GeneralOpertor.ChangeLandUseCount(simulationStartImage[i, j], -10000, listLanduseInfoAndCount);
}
}
//0.3显示输出结果窗体和图表窗体
dockableWindowGraphy.GraphTitle = resourceManager.GetString("String40");
dockableWindowGraphy.XAxisTitle = resourceManager.GetString("String41");
dockableWindowGraphy.YAxisTitle = resourceManager.GetString("String42");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String43"));
Application.DoEvents();
//0.4绘制初始的图表
List <string> listPointPairListName = new List <string>();
List <string> notToDrawList = new List <string>();
notToDrawList.Add(resourceManager.GetString("String44"));
dockableWindowGraphy.CreateGraph(listLanduseInfoAndCount, notToDrawList, out listPointPairListName);
dockableWindowGraphy.RefreshGraph();
int convertedCellCount = 0; //模拟中总共转换的元胞数量
int randomRow, randomColumn; //Monte Carlo方法选取的随机行和随机列
int convertCountInOneIteration = formANNCAWizard.ConvertCount / formANNCAWizard.Iterations; //每次迭代应转换的数量
int convertCountOnce = 0; //每次迭代已经转换的数量
float oldValue, newValue; //每次转换前土地利用类型的新值和旧值
int iteration = 0; //迭代的次数
List <float> listLanduseValues = new List <float>(); //记录土地利用类型值的List
for (int j = 0; j < formANNCAWizard.LandUseClassificationInfo.AllTypesCount; j++)
{
listLanduseValues.Add(formANNCAWizard.LandUseClassificationInfo.AllTypes[j].LanduseTypeValue);
}
int neiWindowSize = formANNCAWizard.NeighbourWindowSize;
int inputNeuronsCount = formANNCAWizard.InputNeuronsCount;
int outputNeuronsCount = formANNCAWizard.OutputNeuronsCount;
//2.开始进行转换
while (convertedCellCount < formANNCAWizard.ConvertCount)
{
convertCountOnce = 0;
//2.3选择元胞与随机比较进行转换
//完成一次迭代
while (convertCountOnce < convertCountInOneIteration)
{
//随机选择一个栅格进行计算
randomRow = random.Next(0, rowCount);
randomColumn = random.Next(0, columnCount);
//计算逻辑为：
//这里需要首先获取当前元胞的输入值
//然后通过神经网络计算每种土地利用类型的概率，乘以随机量
//然后选出最高的概率
//再查看是否是同样的土地利用类型，如果不是则转变
//如果不是并且概率大于一个值，则进行转变（目前先不用）
//如果是空值，则不进行计算
if (simulationImage[randomRow, randomColumn] == -9999f)
{
continue;
}
double[] tempInputsArray = new double[inputNeuronsCount];
double[] tempOutputsArray = new double[outputNeuronsCount];
//获取该栅格的空间变量值
for (int i = 0; i < formANNCAWizard.ListVaribaleImages.Count; i++)
{
tempInputsArray[i] = formANNCAWizard.ListVaribaleImages[i][randomRow, randomColumn];
}
//获取该栅格的邻域值
NeighbourOperator neighbourOperator = new NeighbourOperator();
float[] counts = neighbourOperator.GetNeighbourCount(simulationImage, randomRow, randomColumn, neiWindowSize,
rowCount, columnCount, formANNCAWizard.LandUseClassificationInfo.AllTypesCount, listLanduseValues);
for (int z = 0; z < counts.Length; z++)
{
tempInputsArray[formANNCAWizard.ListVaribaleImages.Count + z] = counts[z] / (neiWindowSize * neiWindowSize - 1);
}
//获取该栅格的土地利用类型
tempInputsArray[inputNeuronsCount - 1] = simulationImage[randomRow, randomColumn];
//计算神经网络输出的概率值
double[] output = formANNCAWizard.ANN.ANNActivationNetwork.Compute(tempInputsArray);
double maxOutput = 0;
int maxIndex = -1;
Random r = new Random(1);
double gamma;
double disturbance;
double probability;
double[] probabilityArray = new double[output.Length];
for (int k = 0; k < output.Length; k++)
{
gamma = random.NextDouble();
disturbance = 1 + Math.Pow((-Math.Log(gamma)), formANNCAWizard.Delta);
probabilityArray[k] = disturbance * output[k];
if (maxOutput < probabilityArray[k])
{
maxOutput = probabilityArray[k];
maxIndex = k;
}
}
probability = maxOutput; //计算最高值的最终概率
//未使用邻域因素模拟较为散，所以加入邻域因素。但可能导致难以达到需要的转换数量，模拟时间过长。
double neiValue = 0f;
for (int t = 0; t < listLanduseValues.Count; t++)
{
if (listLanduseValues[t] == formANNCAWizard.LandUseClassificationInfo.UrbanValues[0].LanduseTypeValue)
{
neiValue = Convert.ToDouble(counts[t]) / (neiWindowSize * neiWindowSize - 1);
break;
}
}
probability = probability * neiValue;
//得到应转变为的土地利用类型值，以及当前土地利用类型值
newValue = formANNCAWizard.LandUseClassificationInfo.AllTypes[maxIndex].LanduseTypeValue;
oldValue = simulationImage[randomRow, randomColumn];
//oldValue有时取到0值等不正确的值
bool isOldValueCorrect = false;
for (int w = 0; w < formANNCAWizard.LandUseClassificationInfo.AllTypes.Count; w++)
{
if (oldValue == formANNCAWizard.LandUseClassificationInfo.AllTypes[w].LanduseTypeValue)
{
isOldValueCorrect = true;
}
}
if (!isOldValueCorrect)
{
continue;
}
//根据转换矩阵进行判断是否可以转换
DataTable dtMatrix = formANNCAWizard.DTMatrix;
int oldValueIndex = -1;
for (int k = 0; k < formANNCAWizard.LandUseClassificationInfo.AllTypesCount; k++)
{
if (formANNCAWizard.LandUseClassificationInfo.AllTypes[k].LanduseTypeValue == oldValue)
{
oldValueIndex = k;
break;
}
}
string canConvert = dtMatrix.Rows[oldValueIndex][maxIndex].ToString();
if (canConvert == "0")
{
continue;
}
//与设定的阈值进行比较
double convertThreshold = formANNCAWizard.ConvertThreshold;
if ((oldValue != newValue) && (probability >= convertThreshold))
{
simulationImage[randomRow, randomColumn] = newValue;
CommonLibrary.GeneralOpertor.ChangeLandUseCount(newValue, oldValue, listLanduseInfoAndCount);
//以城市用地的转变量为总转变量的标准
for (int w = 0; w < formANNCAWizard.LandUseClassificationInfo.UrbanValues.Count; w++)
{
if (newValue == formANNCAWizard.LandUseClassificationInfo.UrbanValues[w].LanduseTypeValue)
{
convertCountOnce++;
convertedCellCount++;
}
}
//System.Diagnostics.Debug.WriteLine(convertedCellCount + " - Old: " + oldValue + " New: " + newValue);
}
}
iteration++;
System.Diagnostics.Debug.WriteLine("iteration: " + iteration);
//2.4.刷新外部界面并输出中间结果数据
if (convertedCellCount == 1 || (iteration % formANNCAWizard.RefreshInterval == 0 && convertedCellCount != 0))
{
//刷新图像
formANNCAWizard.SimulationImage = simulationImage;
VariableMaintainer.IsNeedRefresh = true;
//刷新图表窗体
string landuseTypeName = "";
for (int k = 0; k < listLanduseInfoAndCount.Count; k++)
{
for (int l = 0; l < listPointPairListName.Count; l++)
{
if (System.Globalization.CultureInfo.CurrentCulture.Parent.Name == "zh-CHS")
{
landuseTypeName = listLanduseInfoAndCount[k].structLanduseInfo.LanduseTypeChsName;
}
else if (System.Globalization.CultureInfo.CurrentCulture.Parent.Name == "zh-CHT")
{
landuseTypeName = listLanduseInfoAndCount[k].structLanduseInfo.LanduseTypeChtName;
}
else
{
landuseTypeName = listLanduseInfoAndCount[k].structLanduseInfo.LanduseTypeEnName;
}
if (landuseTypeName == listPointPairListName[l])
{
dockableWindowGraphy.UpdateData(iteration, listLanduseInfoAndCount[k].LanduseTypeCount, l);
}
}
}
dockableWindowGraphy.RefreshGraph();
//刷新输出结果窗体
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String45") + iteration.ToString()
+ resourceManager.GetString("String46"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String47") + convertedCellCount.ToString());
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.ScrollTextbox();
Application.DoEvents();
}
//输出中间结果
if (formANNCAWizard.IsOutput && (iteration % formANNCAWizard.OutputImageInterval == 0))
{
GeneralOpertor.WriteDataFloat(formANNCAWizard.OutputFolder + @"\" + GeneralOpertor.GetNowString()
+ "_ann_iterate_" + iteration.ToString() + @".txt", simulationImage, rowCount, columnCount);
}
}
//3.完成模拟，输出结果。
stopWatch.Stop();
VariableMaintainer.IsSimulationFinished = true;
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String48"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String49") +
GeneralOpertor.GetElapsedTimeString(stopWatch.Elapsed));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText("\n");
//修改结果栅格的属性表
IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass();
IRasterWorkspace rasterWorkspace = workspaceFactory.OpenFromFile(
VariableMaintainer.CurrentFormANNCAWizard.OutputFolder, 0) as IRasterWorkspace;
IRasterDataset rasterDataset = rasterWorkspace.OpenRasterDataset(VariableMaintainer.CurrentFoucsMap.get_Layer(0).Name);
IRasterDatasetEdit3 rasterDatasetEdit3 = rasterDataset as IRasterDatasetEdit3;
rasterDatasetEdit3.BuildAttributeTable();
IRasterDataset3 rasterDataset3 = rasterDataset as IRasterDataset3;
rasterDataset3.Refresh();
if (formANNCAWizard.SimulationEndImageName != "")
{
//GeneralOpertor.WriteDataFloat(formLogisticCAWizard.OutputFolder + @"\CA_ANN_Reslut" + GeneralOpertor.GetNowString() + ".txt",
// simulationImage, rowCount,columnCount);
DataTable dtMatrixNumber = GeneralOpertor.GetMultiTypesMatrix(
simulationImage, simulationEndImage, rowCount, columnCount, formANNCAWizard.LandUseClassificationInfo);
double overallAccuracy = 0d;
double kappa = 0d;
GeneralOpertor.GetMultiTypesAccuracy(dtMatrixNumber, ref overallAccuracy, ref kappa, formANNCAWizard.LandUseClassificationInfo);
FormConfusionMatrix formConfusionMatrix = new FormConfusionMatrix();
formConfusionMatrix.DataTableValues = dtMatrixNumber;
DataTable dtMatrixPercent = dtMatrixNumber.Clone();
GeneralOpertor.CopyDataTableValues(dtMatrixNumber, dtMatrixPercent);
formConfusionMatrix.DataTablePercents = dtMatrixPercent;
formConfusionMatrix.DataGridViewConfusionMatrix.DataSource = dtMatrixNumber;
formConfusionMatrix.LabelOverallAccuracy.Text = (overallAccuracy * 100).ToString("0.00") + " %";
formConfusionMatrix.LabelKappa.Text = kappa.ToString("0.000");
float[] fomValues = GeneralOpertor.GetMultiFoMAccuracy(simulationStartImage, simulationEndImage, simulationImage, rowCount, columnCount);
formConfusionMatrix.LabelFoMValues.Text = "A: " + fomValues[0] + "\nB: " + fomValues[1] +
"\nC: " + fomValues[2] + "\nD: " + fomValues[3];
formConfusionMatrix.LabelFoM.Text = fomValues[4].ToString("0.000");
formConfusionMatrix.LabelPA.Text = fomValues[5].ToString("0.000");
formConfusionMatrix.LabelUA.Text = fomValues[6].ToString("0.000");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String84"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(GeneralOpertor.WriteCoufusionMatrix(dtMatrixPercent));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String83"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String85") + (overallAccuracy * 100).ToString("0.00") + " %");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String86") + kappa.ToString("0.000"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String87") + fomValues[4].ToString("0.000"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String88") + fomValues[5].ToString("0.000"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String89") + fomValues[6].ToString("0.000"));
dockableWindowOutput.AppendText("\n");
formConfusionMatrix.Text = resourceManager.GetString("String105") + " - " + formANNCAWizard.SimulationLayerName;
formConfusionMatrix.ShowDialog();
}
dockableWindowOutput.AppendText("-------------------------------------------");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.ScrollTextbox();
Application.DoEvents();
}
catch (Exception ex)
{
MessageBox.Show(ex.Message);
}
}
public void DoSimulation()
{
//0准备开始模拟
ResourceManager resourceManager = VariableMaintainer.CurrentResourceManager;
FormLogisticCAWizard formLogisticCAWizard = VariableMaintainer.CurrentFormLogisticCAWizard;
Random random = new Random();
Stopwatch stopWatch = new Stopwatch();
stopWatch.Start();
//0.1获得数据
int rowCount = formLogisticCAWizard.CurrentStructRasterMetaData.RowCount;
int columnCount = formLogisticCAWizard.CurrentStructRasterMetaData.ColumnCount;
float[,] simulationStartImage = formLogisticCAWizard.SimulationStartImage;
float[,] simulationEndImage = formLogisticCAWizard.SimulationEndImage;
float[,] simulationImage = formLogisticCAWizard.SimulationImage;
//0.2计算初始每种土地利用类型的单元数量
List <StructLanduseInfoAndCount> listLanduseInfoAndCount; //记录每种土地利用类型的单元数
listLanduseInfoAndCount = new List <StructLanduseInfoAndCount>();
StructLanduseInfoAndCount landuseInfoAndCount;
foreach (StructLanduseInfo structLanduseInfo in formLogisticCAWizard.LandUseClassificationInfo.AllTypes)
{
landuseInfoAndCount = new StructLanduseInfoAndCount();
landuseInfoAndCount.structLanduseInfo = structLanduseInfo;
landuseInfoAndCount.LanduseTypeCount = 0;
listLanduseInfoAndCount.Add(landuseInfoAndCount);
}
for (int i = 0; i < rowCount; i++)
{
for (int j = 0; j < columnCount; j++)
{
CommonLibrary.GeneralOpertor.ChangeLandUseCount(simulationStartImage[i, j], -10000, listLanduseInfoAndCount);
}
}
//0.3显示输出结果窗体和图表窗体
dockableWindowGraphy.GraphTitle = resourceManager.GetString("String40");
dockableWindowGraphy.XAxisTitle = resourceManager.GetString("String41");
dockableWindowGraphy.YAxisTitle = resourceManager.GetString("String42");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String43"));
Application.DoEvents();
//0.4绘制初始的图表
List <string> listPointPairListName = new List <string>();
List <string> notToDrawList = new List <string>();
notToDrawList.Add(resourceManager.GetString("String44"));
dockableWindowGraphy.CreateGraph(listLanduseInfoAndCount, notToDrawList, out listPointPairListName);
dockableWindowGraphy.RefreshGraph();
//1.计算Pg
float[,] pg = new float[rowCount, columnCount];
for (int i = 0; i < rowCount; i++)
{
for (int j = 0; j < columnCount; j++)
{
if (simulationStartImage[i, j] == formLogisticCAWizard.LandUseClassificationInfo.NullValue.LanduseTypeValue)
{
pg[i, j] = formLogisticCAWizard.LandUseClassificationInfo.NullValue.LanduseTypeValue;
}
else
{
pg[i, j] += Convert.ToSingle(formLogisticCAWizard.Coefficents[0]);
for (int k = 0; k < formLogisticCAWizard.ListVariableLayersName.Count; k++)
{
pg[i, j] += Convert.ToSingle(formLogisticCAWizard.Coefficents[k + 1]) *
formLogisticCAWizard.VaribaleImages[k][i, j];
}
}
pg[i, j] = 1 / (1 + Convert.ToSingle(Math.Exp(-1 * pg[i, j])));
}
}
//GeneralOpertor.WriteDataFloat(formLogisticCAWizard.OutputFolder + @"\pg.txt", pg, rowCount, columnCount);
int convertedCellCount = 0; //模拟中总共转换的元胞数量
List <float> neighbours; //Moore领域的元胞
float omega; //扩散系数
float conSuitable; //适应度函数
float maxPct; //每次迭代时最大的Pct值
float[,] pct;
float[,] pst;
float sumPtt = 0; //Ptt的总和
int randomRow, randomColumn; //Monte Carlo方法选取的随机行和随机列
int convertCountInOneIteration = formLogisticCAWizard.ConvertCount / formLogisticCAWizard.Iterations; //每次迭代应转换的数量
int convertCountOnce = 0; //每次迭代已经转换的数量
float oldValue, newValue; //每次转换前土地利用类型的新值和旧值
int[,] changeCells = new int[rowCount, columnCount]; //记录每次转换为城市的单元。
int iteration = 0; //迭代的次数
//2.开始进行转换
while (convertedCellCount < formLogisticCAWizard.ConvertCount)
{
if (convertedCellCount % formLogisticCAWizard.OutputImageInterval == 0 && convertedCellCount != 0)
{
changeCells = new int[rowCount, columnCount];
}
convertCountOnce = 0;
//2.1每次迭代先得到Pct和MaxPct
maxPct = float.MinValue;
pct = new float[rowCount, columnCount];
for (int i = 0; i < rowCount; i++)
{
for (int j = 0; j < columnCount; j++)
{
neighbours = CommonLibrary.GeneralOpertor.GetMooreNeighbors(i, j, simulationImage, rowCount, columnCount);
omega = Convert.ToSingle(CommonLibrary.GeneralOpertor.GetUrbanCount(neighbours, formLogisticCAWizard.LandUseClassificationInfo.UrbanValues[0].LanduseTypeValue)) / (neighbours.Count - 1);
if (formLogisticCAWizard.LandUseClassificationInfo.IsExistInConvertValue(simulationImage[i, j]))
{
conSuitable = 1f;
}
else
{
conSuitable = 0f;
}
pct[i, j] = pg[i, j] * conSuitable * omega;
if (pct[i, j] > maxPct)
{
maxPct = pct[i, j];
}
}
}
//CommonOperator.WriteData(Application.StartupPath + @"\\OutputData\pct.txt", pct, structAsciiImageFileMetaData);
//2.2再得到Pst
pst = new float[rowCount, columnCount];
sumPtt = 0;
//2.2.1先得到ptt和Ptt的总和
for (int i = 0; i < rowCount; i++)
{
for (int j = 0; j < columnCount; j++)
{
pst[i, j] = pct[i, j] * Convert.ToSingle(Math.Exp(-1f * formLogisticCAWizard.Delta * (1f - pct[i, j] / maxPct)));
sumPtt += pst[i, j];
}
}
for (int i = 0; i < rowCount; i++)
{
for (int j = 0; j < columnCount; j++)
{
pst[i, j] = convertCountInOneIteration * pst[i, j] / sumPtt;
}
}
//CommonOperator.WriteData(Application.StartupPath + @"\\OutputData\pst.txt", pst, structAsciiImageFileMetaData);
//2.3选择元胞与随机比较进行转换
//完成一次迭代
while (convertCountOnce < convertCountInOneIteration)
{
randomRow = random.Next(0, rowCount);
randomColumn = random.Next(0, columnCount);
if (pst[randomRow, randomColumn] > Convert.ToSingle(random.NextDouble()))
{
oldValue = simulationImage[randomRow, randomColumn];
newValue = formLogisticCAWizard.LandUseClassificationInfo.UrbanValues[0].LanduseTypeValue;
simulationImage[randomRow, randomColumn] = newValue;
CommonLibrary.GeneralOpertor.ChangeLandUseCount(newValue, oldValue, listLanduseInfoAndCount);
convertCountOnce++;
convertedCellCount++;
changeCells[randomRow, randomColumn] = 1;
}
}
iteration++;
//2.4.刷新外部界面并输出中间结果数据
if (convertedCellCount == 1 || (iteration % formLogisticCAWizard.RefreshInterval == 0 && convertedCellCount != 0))
{
//刷新图像
formLogisticCAWizard.SimulationImage = simulationImage;
VariableMaintainer.IsNeedRefresh = true;
//刷新图表窗体
string landuseTypeName = "";
for (int k = 0; k < listLanduseInfoAndCount.Count; k++)
{
for (int l = 0; l < listPointPairListName.Count; l++)
{
if (System.Globalization.CultureInfo.CurrentCulture.Parent.Name == "zh-CHS")
{
landuseTypeName = listLanduseInfoAndCount[k].structLanduseInfo.LanduseTypeChsName;
}
else if (System.Globalization.CultureInfo.CurrentCulture.Parent.Name == "zh-CHT")
{
landuseTypeName = listLanduseInfoAndCount[k].structLanduseInfo.LanduseTypeChtName;
}
else
{
landuseTypeName = listLanduseInfoAndCount[k].structLanduseInfo.LanduseTypeEnName;
}
if (landuseTypeName == listPointPairListName[l])
{
dockableWindowGraphy.UpdateData(iteration, listLanduseInfoAndCount[k].LanduseTypeCount, l);
}
}
}
dockableWindowGraphy.RefreshGraph();
//刷新输出结果窗体
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String45") + iteration.ToString()
+ resourceManager.GetString("String46"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String47") + convertedCellCount.ToString());
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.ScrollTextbox();
Application.DoEvents();
}
//输出中间结果
if (formLogisticCAWizard.IsOutput && (iteration % formLogisticCAWizard.OutputImageInterval == 0))
{
GeneralOpertor.WriteDataFloat(formLogisticCAWizard.OutputFolder + @"\" + GeneralOpertor.GetNowString()
+ "_lr_iterate_" + iteration.ToString() + @".txt", simulationImage, rowCount, columnCount);
}
}
//3.完成模拟，输出结果。
stopWatch.Stop();
VariableMaintainer.IsSimulationFinished = true;
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String48"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String49") +
GeneralOpertor.GetElapsedTimeString(stopWatch.Elapsed));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText("\n");
//修改结果栅格的属性表
IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass();
IRasterWorkspace rasterWorkspace = workspaceFactory.OpenFromFile(
VariableMaintainer.CurrentFormLogisticCAWizard.OutputFolder, 0) as IRasterWorkspace;
IRasterDataset rasterDataset = rasterWorkspace.OpenRasterDataset(VariableMaintainer.CurrentFoucsMap.get_Layer(0).Name);
IRasterDatasetEdit3 rasterDatasetEdit3 = rasterDataset as IRasterDatasetEdit3;
rasterDatasetEdit3.BuildAttributeTable();
IRasterDataset3 rasterDataset3 = rasterDataset as IRasterDataset3;
rasterDataset3.Refresh();
if (formLogisticCAWizard.SimulationEndImageName != "")
{
//GeneralOpertor.WriteDataFloat(formLogisticCAWizard.OutputFolder + @"\CA_LR_Reslut" + GeneralOpertor.GetNowString() + ".txt",
// simulationImage, rowCount, columnCount);
StructBinaryConfusionMatrix structConfusionMatrix = GeneralOpertor.GetBinaryAccuracy(
simulationImage, simulationEndImage, rowCount, columnCount, formLogisticCAWizard.LandUseClassificationInfo);
string accuracyString = GeneralOpertor.GetBinaryAccuracyReportString(structConfusionMatrix, convertedCellCount);
dockableWindowOutput.AppendText(accuracyString);
DataTable dtMatrixNumber = GeneralOpertor.GetMultiTypesMatrix(
simulationImage, simulationEndImage, rowCount, columnCount, formLogisticCAWizard.LandUseClassificationInfo);
double overallAccuracy = 0d;
double kappa = 0d;
GeneralOpertor.GetMultiTypesAccuracy(dtMatrixNumber, ref overallAccuracy, ref kappa, formLogisticCAWizard.LandUseClassificationInfo);
FormConfusionMatrix formConfusionMatrix = new FormConfusionMatrix();
formConfusionMatrix.DataTableValues = dtMatrixNumber;
DataTable dtMatrixPercent = dtMatrixNumber.Clone();
GeneralOpertor.CopyDataTableValues(dtMatrixNumber, dtMatrixPercent);
formConfusionMatrix.DataTablePercents = dtMatrixPercent;
formConfusionMatrix.DataGridViewConfusionMatrix.DataSource = dtMatrixNumber;
formConfusionMatrix.LabelOverallAccuracy.Text = (overallAccuracy * 100).ToString("0.00") + " %";
formConfusionMatrix.LabelKappa.Text = kappa.ToString("0.000");
float[] fomValues = GeneralOpertor.GetBinaryFoMAccuracy(simulationStartImage, simulationEndImage, simulationImage,
rowCount, columnCount, formLogisticCAWizard.LandUseClassificationInfo.UrbanValues[0].LanduseTypeValue);
formConfusionMatrix.LabelFoMValues.Text = "A: " + fomValues[0] + "\nB: " + fomValues[1] +
"\nC: " + fomValues[2] + "\nD: " + fomValues[3];
formConfusionMatrix.LabelFoM.Text = fomValues[4].ToString("0.000");
formConfusionMatrix.LabelPA.Text = fomValues[5].ToString("0.000");
formConfusionMatrix.LabelUA.Text = fomValues[6].ToString("0.000");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String84"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(GeneralOpertor.WriteCoufusionMatrix(dtMatrixPercent));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String83"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String85") + (overallAccuracy * 100).ToString("0.00") + " %");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String86") + kappa.ToString("0.000"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String87") + fomValues[4].ToString("0.000"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String88") + fomValues[5].ToString("0.000"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String89") + fomValues[6].ToString("0.000"));
dockableWindowOutput.AppendText("\n");
formConfusionMatrix.Text = resourceManager.GetString("String105") + " - " + formLogisticCAWizard.SimulationLayerName;
formConfusionMatrix.ShowDialog();
dockableWindowOutput.ScrollTextbox();
Application.DoEvents();
}
dockableWindowOutput.AppendText("-------------------------------------------");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.ScrollTextbox();
Application.DoEvents();
}
public void DoSimulation()
{
//0准备开始模拟
ResourceManager resourceManager = VariableMaintainer.CurrentResourceManager;
FormDTCAWizard formDTCAWizard = VariableMaintainer.CurrentFormDTCAWizard;
Random random = new Random();
Stopwatch stopWatch = new Stopwatch();
stopWatch.Start();
//0.1获得数据
int rowCount = formDTCAWizard.CurrentStructRasterMetaData.RowCount;
int columnCount = formDTCAWizard.CurrentStructRasterMetaData.ColumnCount;
float[,] simulationStartImage = formDTCAWizard.SimulationStartImage;
float[,] simulationEndImage = formDTCAWizard.SimulationEndImage;
float[,] simulationImage = formDTCAWizard.SimulationImage;
//0.2计算初始每种土地利用类型的单元数量
List <StructLanduseInfoAndCount> listLanduseInfoAndCount; //记录每种土地利用类型的单元数
listLanduseInfoAndCount = new List <StructLanduseInfoAndCount>();
StructLanduseInfoAndCount landuseInfoAndCount;
foreach (StructLanduseInfo structLanduseInfo in formDTCAWizard.LandUseClassificationInfo.AllTypes)
{
landuseInfoAndCount = new StructLanduseInfoAndCount();
landuseInfoAndCount.structLanduseInfo = structLanduseInfo;
landuseInfoAndCount.LanduseTypeCount = 0;
listLanduseInfoAndCount.Add(landuseInfoAndCount);
}
for (int i = 0; i < rowCount; i++)
{
for (int j = 0; j < columnCount; j++)
{
CommonLibrary.GeneralOpertor.ChangeLandUseCount(simulationStartImage[i, j], -10000, listLanduseInfoAndCount);
}
}
//0.3显示输出结果窗体和图表窗体
dockableWindowGraphy.GraphTitle = resourceManager.GetString("String40");
dockableWindowGraphy.XAxisTitle = resourceManager.GetString("String41");
dockableWindowGraphy.YAxisTitle = resourceManager.GetString("String42");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String43"));
Application.DoEvents();
//0.4绘制初始的图表
List <string> listPointPairListName = new List <string>();
List <string> notToDrawList = new List <string>();
notToDrawList.Add(resourceManager.GetString("String44"));
dockableWindowGraphy.CreateGraph(listLanduseInfoAndCount, notToDrawList, out listPointPairListName);
dockableWindowGraphy.RefreshGraph();
int convertedCellCount = 0; //模拟中总共转换的元胞数量
int randomRow, randomColumn; //Monte Carlo方法选取的随机行和随机列
int convertCountInOneIteration = formDTCAWizard.ConvertCount / formDTCAWizard.Iterations; //每次迭代应转换的数量
int convertCountOnce = 0; //每次迭代已经转换的数量
float oldValue, newValue; //每次转换前土地利用类型的新值和旧值
int iteration = 0; //迭代的次数
//float gamma = 1f / formDTCAWizard.Iterations;
//float gamma = 0.4f;
DecisionTree decisionTree = formDTCAWizard.CurrentDecisionTree;
List <float> listUrbanValues = new List <float>();
for (int j = 0; j < formDTCAWizard.LandUseClassificationInfo.UrbanValues.Count; j++)
{
listUrbanValues.Add(formDTCAWizard.LandUseClassificationInfo.UrbanValues[j].LanduseTypeValue);
}
int neiWindowSize = formDTCAWizard.NeighbourWindowSize;
try
{
//2.开始进行转换
while (convertedCellCount < formDTCAWizard.ConvertCount)
{
convertCountOnce = 0;
while (convertCountOnce < convertCountInOneIteration)
{
//随机选择一个栅格进行计算
randomRow = random.Next(0, rowCount);
randomColumn = random.Next(0, columnCount);
//计算逻辑为：
//这里需要首先获取当前元胞的输入值
//然后通过决策树计算是否发生转变
//然后看随机数是否小于等于预定值，则进行转变
//如果是空值，则不进行计算
if (simulationImage[randomRow, randomColumn] == -9999f)
{
continue;
}
//如果模拟影像该栅格是空值，也不进行计算
if (formDTCAWizard.ListVaribaleImages[formDTCAWizard.ListVaribaleImages.Count - 1][randomRow, randomColumn] == -9999f)
{
continue;
}
//如果已经是城市用地，则不计算
if (simulationImage[randomRow, randomColumn] == formDTCAWizard.LandUseClassificationInfo.UrbanValues[0].LanduseTypeValue)
{
continue;
}
double[] tempInputsArray = new double[formDTCAWizard.ListVaribaleImages.Count];
//获取该栅格的各空间变量值和其他输入变量
for (int i = 0; i < formDTCAWizard.ListVaribaleImages.Count; i++)
{
tempInputsArray[i] = formDTCAWizard.ListVaribaleImages[i][randomRow, randomColumn];
}
//获取该栅格的邻域中城市用地数量
tempInputsArray[formDTCAWizard.ListVaribaleImages.Count - 2] =
CommonLibrary.GeneralOpertor.GetNeighbors(randomRow, randomColumn, simulationImage,
rowCount, columnCount, neiWindowSize, listUrbanValues);
//获取该栅格当前的土地利用类型值
tempInputsArray[formDTCAWizard.ListVaribaleImages.Count - 1] = simulationImage[randomRow, randomColumn];
//计算决策树输出的转换值
int output = formDTCAWizard.CurrentDecisionTree.Compute(tempInputsArray);
//如果计算为不转换，则不进行转换
if (output == 0)
{
continue;
}
//如果可以转换为城市用地
if (!formDTCAWizard.LandUseClassificationInfo.IsExistInConvertValue(simulationImage[randomRow, randomColumn]))
{
continue;
}
Random r = new Random(1);
double convertThreshold = random.NextDouble();
//得到应转变为的土地利用类型值，以及当前土地利用类型值
newValue = formDTCAWizard.LandUseClassificationInfo.UrbanValues[0].LanduseTypeValue;
oldValue = simulationImage[randomRow, randomColumn];
//未使用邻域因素模拟较为散，所以加入邻域因素
double neighbourValue = tempInputsArray[formDTCAWizard.ListVaribaleImages.Count - 2] / (neiWindowSize * neiWindowSize - 1);
//double probability = convertThreshold * neighbourValue;
//if ((oldValue != newValue) && (convertThreshold <= gamma))
//if ((oldValue != newValue) && (probability >= gamma))
//20170619添加限制层数据
// if ((VariableMaintainer.RestrictImage[randomRow, randomColumn] == 2) || (VariableMaintainer.RestrictImage[randomRow, randomColumn] == 1))
//if (VariableMaintainer.RestrictImage[randomRow, randomColumn] == 1)
// continue;
//引入随机因素，如果邻域因素大于随机数
if ((oldValue != newValue) && (neighbourValue > convertThreshold))
{
simulationImage[randomRow, randomColumn] = newValue;
CommonLibrary.GeneralOpertor.ChangeLandUseCount(newValue, oldValue, listLanduseInfoAndCount);
convertCountOnce++;
convertedCellCount++;
//System.Diagnostics.Debug.WriteLine(convertedCellCount + " - Old: " + oldValue + " New: " + newValue);
}
}
iteration++;
//2.4.刷新外部界面并输出中间结果数据
if (convertedCellCount == 1 || (iteration % formDTCAWizard.RefreshInterval == 0 && convertedCellCount != 0))
{
//刷新图像
formDTCAWizard.SimulationImage = simulationImage;
VariableMaintainer.IsNeedRefresh = true;
//刷新图表窗体
string landuseTypeName = "";
for (int k = 0; k < listLanduseInfoAndCount.Count; k++)
{
for (int l = 0; l < listPointPairListName.Count; l++)
{
if (System.Globalization.CultureInfo.CurrentCulture.Parent.Name == "zh-CHS")
{
landuseTypeName = listLanduseInfoAndCount[k].structLanduseInfo.LanduseTypeChsName;
}
else if (System.Globalization.CultureInfo.CurrentCulture.Parent.Name == "zh-CHT")
{
landuseTypeName = listLanduseInfoAndCount[k].structLanduseInfo.LanduseTypeChtName;
}
else
{
landuseTypeName = listLanduseInfoAndCount[k].structLanduseInfo.LanduseTypeEnName;
}
if (landuseTypeName == listPointPairListName[l])
{
dockableWindowGraphy.UpdateData(iteration, listLanduseInfoAndCount[k].LanduseTypeCount, l);
}
}
}
dockableWindowGraphy.RefreshGraph();
//刷新输出结果窗体
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String45") + iteration.ToString()
+ resourceManager.GetString("String46"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String47") + convertedCellCount.ToString());
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.ScrollTextbox();
Application.DoEvents();
}
//输出中间结果
if (formDTCAWizard.IsOutput && (iteration % formDTCAWizard.OutputImageInterval == 0))
{
GeneralOpertor.WriteDataFloat(formDTCAWizard.OutputFolder + @"\" + GeneralOpertor.GetNowString()
+ "_ann_iterate_" + iteration.ToString() + @".txt", simulationImage, rowCount, columnCount);
}
}
}
catch (Exception ex)
{
MessageBox.Show("DTCA: " + ex.Message);
}
//3.完成模拟，输出结果。
stopWatch.Stop();
//isFinished = true;
VariableMaintainer.IsSimulationFinished = true;
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String48"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String49") +
GeneralOpertor.GetElapsedTimeString(stopWatch.Elapsed));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText("\n");
//修改结果栅格的属性表
IWorkspaceFactory workspaceFactory = new RasterWorkspaceFactoryClass();
IRasterWorkspace rasterWorkspace = workspaceFactory.OpenFromFile(
VariableMaintainer.CurrentFormDTCAWizard.OutputFolder, 0) as IRasterWorkspace;
IRasterDataset rasterDataset = rasterWorkspace.OpenRasterDataset(VariableMaintainer.CurrentFoucsMap.get_Layer(0).Name);
IRasterDatasetEdit3 rasterDatasetEdit3 = rasterDataset as IRasterDatasetEdit3;
rasterDatasetEdit3.BuildAttributeTable();
IRasterDataset3 rasterDataset3 = rasterDataset as IRasterDataset3;
rasterDataset3.Refresh();
if (formDTCAWizard.SimulationEndImageName != "")
{
//GeneralOpertor.WriteDataFloat(formLogisticCAWizard.OutputFolder + @"\CA_ANN_Reslut" + GeneralOpertor.GetNowString() + ".txt",
// simulationImage, rowCount,columnCount);
StructBinaryConfusionMatrix structConfusionMatrix = GeneralOpertor.GetBinaryAccuracy(
simulationImage, simulationEndImage, rowCount, columnCount, formDTCAWizard.LandUseClassificationInfo);
string accuracyString = GeneralOpertor.GetBinaryAccuracyReportString(structConfusionMatrix, convertedCellCount);
dockableWindowOutput.AppendText(accuracyString);
DataTable dtMatrixNumber = GeneralOpertor.GetMultiTypesMatrix(
simulationImage, simulationEndImage, rowCount, columnCount, formDTCAWizard.LandUseClassificationInfo);
double overallAccuracy = 0d;
double kappa = 0d;
GeneralOpertor.GetMultiTypesAccuracy(dtMatrixNumber, ref overallAccuracy, ref kappa, formDTCAWizard.LandUseClassificationInfo);
FormConfusionMatrix formConfusionMatrix = new FormConfusionMatrix();
formConfusionMatrix.DataTableValues = dtMatrixNumber;
DataTable dtMatrixPercent = dtMatrixNumber.Clone();
GeneralOpertor.CopyDataTableValues(dtMatrixNumber, dtMatrixPercent);
formConfusionMatrix.DataTablePercents = dtMatrixPercent;
formConfusionMatrix.DataGridViewConfusionMatrix.DataSource = dtMatrixNumber;
formConfusionMatrix.LabelOverallAccuracy.Text = (overallAccuracy * 100).ToString("0.00") + " %";
formConfusionMatrix.LabelKappa.Text = kappa.ToString("0.000");
float[] fomValues = GeneralOpertor.GetBinaryFoMAccuracy(simulationStartImage, simulationEndImage, simulationImage,
rowCount, columnCount, formDTCAWizard.LandUseClassificationInfo.UrbanValues[0].LanduseTypeValue);
formConfusionMatrix.LabelFoMValues.Text = "A: " + fomValues[0] + "\nB: " + fomValues[1] +
"\nC: " + fomValues[2] + "\nD: " + fomValues[3];
formConfusionMatrix.LabelFoM.Text = fomValues[4].ToString("0.000");
formConfusionMatrix.LabelPA.Text = fomValues[5].ToString("0.000");
formConfusionMatrix.LabelUA.Text = fomValues[6].ToString("0.000");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String84"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(GeneralOpertor.WriteCoufusionMatrix(dtMatrixPercent));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String83"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String85") + (overallAccuracy * 100).ToString("0.00") + " %");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String86") + kappa.ToString("0.000"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String87") + fomValues[4].ToString("0.000"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String88") + fomValues[5].ToString("0.000"));
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText(resourceManager.GetString("String89") + fomValues[6].ToString("0.000"));
dockableWindowOutput.AppendText("\n");
formConfusionMatrix.Text = resourceManager.GetString("String105") + " - " + formDTCAWizard.SimulationLayerName;
formConfusionMatrix.ShowDialog();
}
dockableWindowOutput.AppendText("-------------------------------------------");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.AppendText("\n");
dockableWindowOutput.ScrollTextbox();
Application.DoEvents();
}
private void simpleButton4_Click(object sender, EventArgs e)
{
if (this.comboBoxEditInputRaster1.Text.ToString() == "" || this.comboBoxEditInputRaster2.Text.ToString() == "" || this.comboBoxEditOutputRaster.Text.ToString() == "")
{
MessageBox.Show("请输入参数！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
string strResultsDBPath = m_pGDBHelper.GetResultsDBPath();
if (Utilities.GDBUtilites.CheckNameExist(strResultsDBPath, this.comboBoxEditOutputRaster.Text.ToString()) == false)
{
//执行PCA分析
ILayer pInputLayer1 = null;
ILayer pInputLayer2 = null;
for (int i = 0; i < m_pMapCtrl.LayerCount; i++)
{
if (m_pMapCtrl.get_Layer(i).Name == this.comboBoxEditInputRaster1.Text.ToString())
{
pInputLayer1 = m_pMapCtrl.get_Layer(i);
}
else if (m_pMapCtrl.get_Layer(i).Name == this.comboBoxEditInputRaster2.Text.ToString())
{
pInputLayer2 = m_pMapCtrl.get_Layer(i);
}
}
if (pInputLayer1 != null && pInputLayer2 != null)
{
IRasterLayer pRasterLayer1 = pInputLayer1 as IRasterLayer;
IRasterLayer pRasterLayer2 = pInputLayer2 as IRasterLayer;
string RasterPath1 = pRasterLayer1.FilePath;
string RasterPath2 = pRasterLayer2.FilePath;
IGeoProcessor2 pGP = new GeoProcessorClass();
pGP.AddToolbox(m_pGDBHelper.GetToolboxPath());
IVariantArray gpParameters = new VarArrayClass();
//gpParameters.Add(pRasterLayer1.Raster);
//gpParameters.Add(pRasterLayer2.Raster);
gpParameters.Add(RasterPath1);
gpParameters.Add(RasterPath2);
gpParameters.Add(String.Format("{0}\\{1}", m_pGDBHelper.GetResultsDBPath(), this.comboBoxEditOutputRaster.Text));
string ToolName = "ReseePACBandMath";
ToolName = "NewDiffPCA";
IGeoProcessorResult pGeoProcessorResult = pGP.Execute(ToolName, gpParameters, null);
if (pGeoProcessorResult.Status == esriJobStatus.esriJobSucceeded)
{
if (this.checkBox1.Checked)
{
IWorkspaceFactory2 pWKF = new FileGDBWorkspaceFactoryClass();
IRasterWorkspaceEx pRasterWKEx = (IRasterWorkspaceEx)pWKF.OpenFromFile(m_pGDBHelper.GetResultsDBPath(), 0);
IRasterDataset3 pRasterDataset = pRasterWKEx.OpenRasterDataset(this.comboBoxEditOutputRaster.Text.ToString()) as IRasterDataset3;
Utilities.MapUtilites.AddRasterLayer(m_pMapCtrl.ActiveView, pRasterDataset, null);
}
MessageBox.Show("PCA分析完成！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
}
}
this.DialogResult = DialogResult.OK;
this.Close();
}
else
{
MessageBox.Show("已存在同名数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
}
/// <summary>
/// Sets band information on items in a mosaic dataset
/// </summary>
/// <param name="md">The mosaic dataset with the items</param>
private static void SetMosaicDatasetItemInformation(IMosaicDataset md)
{
// Get the Attribute table from the Mosaic Dataset.
IFeatureClass featureClass = md.Catalog;
ISchemaLock schemaLock = (ISchemaLock)featureClass;
IRasterDataset3 rasDs = null;
try
{
// A try block is necessary, as an exclusive lock might not be available.
schemaLock.ChangeSchemaLock(esriSchemaLock.esriExclusiveSchemaLock);
// Get an update cursor going through all the rows in the Moasic Dataset.
IFeatureCursor fcCursor = featureClass.Update(null, false);
// Alternatively, a read cursor can be used if the item does not need to be changed.
// featureClass.Search(null, false);
// For each row,
IRasterCatalogItem rasCatItem = (IRasterCatalogItem)fcCursor.NextFeature();
while (rasCatItem != null)
{
// get the functionrasterdataset from the Raster field.
IFunctionRasterDataset funcDs = (IFunctionRasterDataset)rasCatItem.RasterDataset;
if (funcDs != null)
{
// Check if the 'BandName' property exists in the dataset.
bool propertyExists = false;
for (int bandID = 0; bandID < funcDs.RasterInfo.BandCount; ++bandID)
{
object bandNameProperty = null;
bandNameProperty = GetBandProperty((IDataset)funcDs, "BandName", bandID);
if (bandNameProperty != null)
{
propertyExists = true;
}
}
if (propertyExists == false && funcDs.RasterInfo.BandCount > 2)
{
// If the property does not exist and the dataset has atleast 3 bands,
// set Band Definition Properties for first 3 bands of the dataset.
SetBandProperties((IDataset)funcDs);
funcDs.AlterDefinition();
rasDs = (IRasterDataset3)funcDs;
// Refresh the dataset.
rasDs.Refresh();
}
}
fcCursor.UpdateFeature((IFeature)rasCatItem);
rasDs = null;
rasCatItem = (IRasterCatalogItem)fcCursor.NextFeature();
}
rasCatItem = null;
fcCursor = null;
featureClass = null;
}
catch (Exception exc) { Console.WriteLine("Exception Caught in SetMosaicDatasetItemInformation: " + exc.Message); }
finally
{
// Set the lock to shared, whether or not an error occurred.
schemaLock.ChangeSchemaLock(esriSchemaLock.esriSharedSchemaLock);
}
}
private void simpleButton4_Click(object sender, EventArgs e)
{
if (this.comboBoxInputRaster.Text.ToString() == "" || this.comboBoxContrastRaster.Text.ToString() == "" || this.comboBoxOutputRaster.Text.ToString() == "")
{
MessageBox.Show("请输入参数！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}//有任意一个空为空则弹出对话框提示：请输入参数。
string strResultsDBPath = m_pGDBHelper.GetResultsDBPath();//声明变量获取数据库结果路径
//string ToolName = "变化向量角5";
if (Utilities.GDBUtilites.CheckNameExist(strResultsDBPath, this.comboBoxOutputRaster.Text.ToString()) == false)//如果没有同名文件
{
//执行光谱角分析
ILayer pInputLayer1 = null;
ILayer pInputLayer2 = null;
for (int i = 0; i < m_pMapCtrl.LayerCount; i++)
{
if (m_pMapCtrl.get_Layer(i).Name == this.comboBoxInputRaster.Text.ToString())
{
pInputLayer1 = m_pMapCtrl.get_Layer(i);
}
else if (m_pMapCtrl.get_Layer(i).Name == this.comboBoxContrastRaster.Text.ToString())
{
pInputLayer2 = m_pMapCtrl.get_Layer(i);
}
}
if (pInputLayer1 != null && pInputLayer2 != null)
{
IRasterLayer pRasterLayer1 = pInputLayer1 as IRasterLayer;
IRasterLayer pRasterLayer2 = pInputLayer2 as IRasterLayer;
string rasterpath1 = pRasterLayer1.FilePath;
string rasterpath2 = pRasterLayer2.FilePath;
IGeoProcessor2 pGP = new GeoProcessorClass();
string GPPath = m_pGDBHelper.GetToolboxPath();
pGP.AddToolbox(GPPath);
IVariantArray gpParameters = new VarArrayClass();
//gpParameters.Add(pRasterLayer1.Raster);
//gpParameters.Add(pRasterLayer2.Raster);
gpParameters.Add(rasterpath1);
gpParameters.Add(rasterpath2);
gpParameters.Add(strResultsDBPath + "\\" + this.comboBoxOutputRaster.Text);
IGeoProcessorResult pGeoProcessorResult = null;
pGeoProcessorResult = pGP.Execute("changevector", gpParameters, null);
if (pGeoProcessorResult.Status == esriJobStatus.esriJobSucceeded)
{
if (this.checkBox1.Checked)
{
IWorkspaceFactory2 pWKF = new FileGDBWorkspaceFactoryClass();
IRasterWorkspaceEx pRasterWKEx = (IRasterWorkspaceEx)pWKF.OpenFromFile(m_pGDBHelper.GetResultsDBPath(), 0);
IRasterDataset3 pRasterDataset = pRasterWKEx.OpenRasterDataset(this.comboBoxOutputRaster.Text.ToString()) as IRasterDataset3;
Utilities.MapUtilites.AddRasterLayer(m_pMapCtrl.ActiveView, pRasterDataset, null);
}
MessageBox.Show("变化向量分析完成！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
}
}
this.DialogResult = DialogResult.OK;
this.Close();
}
else
{
MessageBox.Show("已存在同名数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
}
private ILayer GetILayerFromPath(string path)
{
ILayer pLayer = null;
string[] vectorExtension = { ".SHP", ".E00", ".DWG" };
string[] rasterExtension = { ".TIF", ".IMG", ".PIX", ".JPG", ".BMP" };
//if(System.IO.Path.GetExtension
if (System.IO.Path.GetDirectoryName(path).ToUpper().EndsWith(".GDB"))
{
IWorkspaceFactory2 pWKF = new FileGDBWorkspaceFactoryClass();
IRasterWorkspaceEx pRasterWKEx = (IRasterWorkspaceEx)pWKF.OpenFromFile(System.IO.Path.GetDirectoryName(path), 0);
IRasterDataset3 pRasterDataset = pRasterWKEx.OpenRasterDataset(System.IO.Path.GetFileName(path)) as IRasterDataset3;
}
else
{
if (rasterExtension.Contains(System.IO.Path.GetExtension(path).ToUpper()))
{
IWorkspaceFactory pWorkspaceFactory;
IRasterWorkspace pRasterWorkspace;
// 获取当前路径和文件名
if (path == "")
{
return(null);
}
int Index = path.LastIndexOf("//");
if (Index == -1)
{
Index = path.LastIndexOf("\\");
}
string fileName = path.Substring(Index + 1);
string filePath = path.Substring(0, Index);
pWorkspaceFactory = new RasterWorkspaceFactoryClass();
pRasterWorkspace = (IRasterWorkspace)pWorkspaceFactory.OpenFromFile(filePath, 0);
IRasterDataset pRasterDataset = (IRasterDataset)pRasterWorkspace.OpenRasterDataset(fileName);
IRasterLayer pRasterLayer = new RasterLayerClass();
pRasterLayer.CreateFromDataset(pRasterDataset);
return((ILayer)pRasterLayer);
}
else if (vectorExtension.Contains(System.IO.Path.GetExtension(path).ToUpper()))
{
IWorkspaceFactory pWorkspaceFactory;
IFeatureWorkspace pFeatureWorkspace;
IFeatureLayer pFeatureLayer;
if (path == "")
{
return(null);
}
int Index = path.LastIndexOf("\\");
string filePath = path.Substring(0, Index);
string fileName = path.Substring(Index + 1);
// 打开工作空间并添加shp文件
pWorkspaceFactory = new ShapefileWorkspaceFactoryClass();
pFeatureWorkspace = (IFeatureWorkspace)pWorkspaceFactory.OpenFromFile(filePath, 0);
pFeatureLayer = new FeatureLayerClass();
pFeatureLayer.FeatureClass = pFeatureWorkspace.OpenFeatureClass(fileName);
pFeatureLayer.Name = pFeatureLayer.FeatureClass.AliasName;
return((ILayer)pFeatureLayer);
}
}
return(pLayer);
}
