public string RasterCalculatorNKS(string nN, string nP, string Y, string S, string gS)
{
try
{
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator rasterCalculator = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
if (nP == "12")
{
rasterCalculator.expression = "('" + Y + "' * '" + nN + "' * 12) / ('" + S + "' * '" + gS + "' + 15 )";
}
else
{
rasterCalculator.expression = "('" + Y + "' * '" + nN + "' * '" + nP + "') / ('" + S + "' * '" + gS + "' + 15)";
}
ESRI.ArcGIS.Geoprocessor.Geoprocessor gp = new ESRI.ArcGIS.Geoprocessor.Geoprocessor();
rasterCalculator.output_raster = AppSingleton.Instance().WorkspacePath + "\\nks";
gp.AddOutputsToMap = AppSingleton.Instance().AralariEkle;
gp.OverwriteOutput = true;
gp.Execute(rasterCalculator, null);
return(rasterCalculator.output_raster.ToString());
}
catch (Exception ex)
{
return(string.Empty);
}
}
public string RasterCalculatorYillik(string P, string T, string layerName)
{
try
{
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator rasterCalculator = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
rasterCalculator.expression = "'" + P + "'/('" + T + "' + 10)";
ESRI.ArcGIS.Geoprocessor.Geoprocessor gp = new ESRI.ArcGIS.Geoprocessor.Geoprocessor();
rasterCalculator.output_raster = AppSingleton.Instance().WorkspacePath + "\\" + layerName;
gp.AddOutputsToMap = true;//AppSingleton.Instance().AralariEkle;
gp.OverwriteOutput = true;
gp.Execute(rasterCalculator, null);
return(rasterCalculator.output_raster.ToString());
}
catch (Exception ex)
{
return(string.Empty);
}
}
public string RasterCalculatorKi(string d, string etp, string layerName)
{
try
{
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator rasterCalculator = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
rasterCalculator.expression = "((100 * '" + d + "') / '" + etp + "')*10";
ESRI.ArcGIS.Geoprocessor.Geoprocessor gp = new ESRI.ArcGIS.Geoprocessor.Geoprocessor();
rasterCalculator.output_raster = AppSingleton.Instance().WorkspacePath + "\\" + layerName;
gp.AddOutputsToMap = AppSingleton.Instance().AralariEkle;
gp.OverwriteOutput = true;
gp.Execute(rasterCalculator, null);
return(rasterCalculator.output_raster.ToString());
}
catch (Exception ex)
{
return(string.Empty);
}
}
public string RasterCalculatorKKS(string nKS)
{
try
{
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator rasterCalculator = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
rasterCalculator.expression = "1 / '" + nKS + "'";
ESRI.ArcGIS.Geoprocessor.Geoprocessor gp = new ESRI.ArcGIS.Geoprocessor.Geoprocessor();
rasterCalculator.output_raster = AppSingleton.Instance().WorkspacePath + "\\kks";
gp.AddOutputsToMap = AppSingleton.Instance().AralariEkle;
gp.OverwriteOutput = true;
gp.Execute(rasterCalculator, null);
return(rasterCalculator.output_raster.ToString());
}
catch (Exception ex)
{
return(string.Empty);
}
}
private void bt_ok_Click(object sender, EventArgs e)
{
//判断输入路径是否正确
#region
try
{
if (!File.Exists(tb_demnput.Text))
{
MessageBox.Show("平均降水量数据输入路径不正确或文件被占用，请重新输入路径！", "提示", MessageBoxButtons.OK);
return;
}
if (Convert.ToDouble(tb_highinput.Text) < 0)
{
MessageBox.Show("请输入正确的水位高程值！", "提示", MessageBoxButtons.OK);
return;
}
}
catch (Exception ex)
{
MessageBox.Show("设置路径或水位高程值不合法，请检查！");
return;
}
#endregion
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator rc = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
Geoprocessor gp = new Geoprocessor();
gp.OverwriteOutput = true;
object sev = null;
try
{
tb_state.Text = "正在处理……";
//判断输入水位数据的单位
double waterunit = 0;
if (cb_waterunits.Text == "米（m）")
{
waterunit = 1;
}
else
{
waterunit = 0.001;
}
//判断输入缓冲距离及缓冲单位
double bufferunit = 0;
if (cb_bufferunits.Text == "米（m）")
{
bufferunit = 1;
}
else
{
bufferunit = 1000;
}
double[] bufferdistance = { Convert.ToDouble(nud_one.Value) * bufferunit, Convert.ToDouble(nud_two.Value) * bufferunit, Convert.ToDouble(nud_three.Value) * bufferunit, Convert.ToDouble(nud_four.Value) * bufferunit };
//提取高程高于水位的区域
tb_state.Text = "正在提取未淹没区域……";
string flood_high = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_high" + System.IO.Path.GetExtension(tb_flooutput.Text);
rc.output_raster = flood_high;
rc.expression = "Con(\"" + tb_demnput.Text + "\" >= (" + tb_highinput.Text + " * " + waterunit.ToString() + "),\"" + tb_demnput.Text + "\")";
gp.Execute(rc, null);
tb_state.Text = "未淹没区域提取完成！";
//水文分析填洼
tb_state.Text = "正在进行水文分析计算填洼……";
ESRI.ArcGIS.SpatialAnalystTools.Fill fill = new ESRI.ArcGIS.SpatialAnalystTools.Fill();
string fill_input, fill_output;
fill_input = flood_high;
fill_output = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_fill" + System.IO.Path.GetExtension(tb_flooutput.Text);
fill.in_surface_raster = fill_input;
fill.out_surface_raster = fill_output;
gp.Execute(fill, null);
tb_state.Text = "水文分析填洼计算完成！";
//水文分析计算流向
tb_state.Text = "正在进行水文分析计算流向……";
ESRI.ArcGIS.SpatialAnalystTools.FlowDirection flowdirection = new ESRI.ArcGIS.SpatialAnalystTools.FlowDirection();
string flowdirection_input, flowdirection_output;
flowdirection_input = fill_output;
flowdirection_output = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_FlowDir" + System.IO.Path.GetExtension(tb_flooutput.Text);
flowdirection.in_surface_raster = flowdirection_input;
flowdirection.out_flow_direction_raster = flowdirection_output;
gp.Execute(flowdirection, null);
tb_state.Text = "水文分析流向计算完成！";
//水文分析计算流量
tb_state.Text = "正在进行水文分析计算流量……";
ESRI.ArcGIS.SpatialAnalystTools.FlowAccumulation flowaccumulation = new ESRI.ArcGIS.SpatialAnalystTools.FlowAccumulation();
string flowaccumulation_input, flowaccumulation_output;
flowaccumulation_input = flowdirection_output;
flowaccumulation_output = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_FlowAcc" + System.IO.Path.GetExtension(tb_flooutput.Text);
flowaccumulation.in_flow_direction_raster = flowaccumulation_input;
flowaccumulation.out_accumulation_raster = flowaccumulation_output;
gp.Execute(flowaccumulation, null);
tb_state.Text = "水文分析流量计算完成！";
//盆域分析
tb_state.Text = "正在进行水文分析盆域分析……";
ESRI.ArcGIS.SpatialAnalystTools.Basin basin = new ESRI.ArcGIS.SpatialAnalystTools.Basin();
string basin_input, basin_output;
basin_input = flowdirection_output;
basin_output = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_Basin" + System.IO.Path.GetExtension(tb_flooutput.Text);
basin.in_flow_direction_raster = basin_input;
basin.out_raster = basin_output;
gp.Execute(basin, null);
tb_state.Text = "水文分析盆域分析计算完成！";
//阈值提取河流
tb_state.Text = "正在阈值提取河流……";
string flood_rasterriver = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_river" + System.IO.Path.GetExtension(tb_flooutput.Text);
rc.expression = "Con(\"" + flowaccumulation_output + "\" >= 800,1)";
rc.output_raster = flood_rasterriver;
gp.Execute(rc, null);
tb_state.Text = "阈值提取河流计算完成！";
//水文分析栅格河网矢量化
tb_state.Text = "正在进行水文分析栅格河网矢量化……";
ESRI.ArcGIS.SpatialAnalystTools.StreamToFeature streamtofeature = new ESRI.ArcGIS.SpatialAnalystTools.StreamToFeature();
string streamtofeature_output;
streamtofeature_output = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_river.shp";
streamtofeature.in_flow_direction_raster = flowdirection_output;
streamtofeature.in_stream_raster = flood_rasterriver;
streamtofeature.out_polyline_features = streamtofeature_output;
gp.Execute(streamtofeature, null);
tb_state.Text = "水文分析栅格河网矢量化计算完成！";
//多环缓冲区
tb_state.Text = "正在计算多环缓冲区……";
ESRI.ArcGIS.AnalysisTools.MultipleRingBuffer multipleringbuffer = new ESRI.ArcGIS.AnalysisTools.MultipleRingBuffer();
string multipleringbuffer_output = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_multbuffer.shp";
multipleringbuffer.Input_Features = streamtofeature_output;
multipleringbuffer.Distances = bufferdistance[0].ToString() + ";" + bufferdistance[1].ToString() + ";" + bufferdistance[2].ToString() + ";" + bufferdistance[3].ToString() + ";";
;
multipleringbuffer.Buffer_Unit = "meters";
multipleringbuffer.Output_Feature_class = multipleringbuffer_output;
gp.Execute(multipleringbuffer, null);
tb_state.Text = "多环缓冲区计算完成！";
//多环缓冲区矢量转栅格
tb_state.Text = "正在进行多环缓冲区矢量转栅格……";
ESRI.ArcGIS.ConversionTools.FeatureToRaster featuretoraster = new ESRI.ArcGIS.ConversionTools.FeatureToRaster();
string featuretoraster_output = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_multbuffer" + System.IO.Path.GetExtension(tb_flooutput.Text);
featuretoraster.in_features = System.IO.Path.GetDirectoryName(tb_flooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "_multbuffer.shp";
featuretoraster.field = "distance";
featuretoraster.out_raster = featuretoraster_output;
gp.Execute(featuretoraster, null);
tb_state.Text = "多环缓冲区矢量转栅格计算完成！";
//删除临时文件
if (tb_state.Text == "多环缓冲区矢量转栅格计算完成！")
{
string[] files = Directory.GetFiles(System.IO.Path.GetDirectoryName(tb_flooutput.Text), System.IO.Path.GetFileNameWithoutExtension(tb_flooutput.Text) + "*");
foreach (string file in files)
{
if (System.IO.Path.GetFileNameWithoutExtension(file) != System.IO.Path.GetFileNameWithoutExtension(basin_output) && System.IO.Path.GetFileNameWithoutExtension(file) != System.IO.Path.GetFileNameWithoutExtension(featuretoraster_output))
{
File.Delete(file);
}
}
}
//将结果加载显示
#region
string mxfile = System.Environment.CurrentDirectory + "\\Data\\EcoRedLine\\洪水淹没生态红线划分.mxd";
IMapDocument pMapDocument = new MapDocumentClass();
pMapDocument.Open(mxfile); //打开本地的地图文档，用来操作改mxd文件
IMap pMap = pMapDocument.get_Map(0);
IMapLayers pMapLayer = pMap as IMapLayers;
IRasterLayer pRasterLayer = new RasterLayerClass();
IWorkspaceFactory rasterWorkspaceFactory = new RasterWorkspaceFactoryClass();
IRasterWorkspace rasterWorkspace = (IRasterWorkspace)rasterWorkspaceFactory.OpenFromFile(System.IO.Path.GetDirectoryName(tb_flooutput.Text), 0);
IRasterDataset pRasterDataset1 = rasterWorkspace.OpenRasterDataset("Flood_multbuffer.tif"); //打开栅格图的文件名
pRasterLayer.CreateFromDataset(pRasterDataset1); //创建
pRasterLayer = pMapLayer.get_Layer(0) as IRasterLayer;
pRasterLayer.CreateFromDataset(pRasterDataset1);
IRasterLayer pRasterLayer2 = new RasterLayerClass();
IRasterDataset pRasterDataset2 = rasterWorkspace.OpenRasterDataset("Flood_Basin.tif"); //打开栅格图的文件名
pRasterLayer2.CreateFromDataset(pRasterDataset2); //创建
pRasterLayer2 = pMapLayer.get_Layer(1) as IRasterLayer;
pRasterLayer2.CreateFromDataset(pRasterDataset2);
pMapDocument.Save(true, true);//保存更改完路径后的mxd文件
_PageLayoutControl.LoadMxFile(mxfile);
_PageLayoutControl.Extent = _PageLayoutControl.FullExtent;
_PageLayoutControl.ZoomToWholePage();
#endregion
}
catch (Exception ex)
{
MessageBox.Show(gp.GetMessages(ref sev), "提示", MessageBoxButtons.OK);
}
}
private void bt_ok_Click(object sender, EventArgs e)
{
//判断输入路径是否正确
#region
try
{
if (!File.Exists(tb_nppinput.Text))
{
MessageBox.Show("生态系统净初级生产力平均值数据输入路径不正确或文件被占用，请重新输入路径！", "提示", MessageBoxButtons.OK);
return;
}
if (!File.Exists(tb_preinput.Text))
{
MessageBox.Show("平均年降水量数据输入路径不正确或文件被占用，请重新输入路径！", "提示", MessageBoxButtons.OK);
return;
}
if (!File.Exists(tb_tminput.Text))
{
MessageBox.Show("平均温度数据输入路径不正确或文件被占用，请重新输入路径！", "提示", MessageBoxButtons.OK);
return;
}
if (!File.Exists(tb_altinput.Text))
{
MessageBox.Show("海拔参数数据输入路径不正确或文件被占用，请重新输入路径！", "提示", MessageBoxButtons.OK);
return;
}
}
catch (Exception ex)
{
MessageBox.Show("设置路径不合法，请检查！");
return;
}
#endregion
object sev = null;
tb_state.Text = "正在处理……";
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator rc = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.Float tofloat = new ESRI.ArcGIS.SpatialAnalystTools.Float();
Geoprocessor gp = new Geoprocessor();
gp.OverwriteOutput = true;
try
{
//输入数据类型转换
#region
//npp输入数据转化为浮点型
string floatfilepath_npp = System.IO.Path.GetDirectoryName(tb_biooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) + "_float_npp" + System.IO.Path.GetExtension(tb_biooutput.Text);
tofloat.in_raster_or_constant = tb_nppinput.Text;
tofloat.out_raster = floatfilepath_npp;
gp.Execute(tofloat, null);
//pre输入数据转化为浮点型
string floatfilepath_pre = System.IO.Path.GetDirectoryName(tb_biooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) + "_float_pre" + System.IO.Path.GetExtension(tb_biooutput.Text);
tofloat.in_raster_or_constant = tb_preinput.Text;
tofloat.out_raster = floatfilepath_pre;
gp.Execute(tofloat, null);
//tm输入数据转化为浮点型
string floatfilepath_tm = System.IO.Path.GetDirectoryName(tb_biooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) + "_float_tm" + System.IO.Path.GetExtension(tb_biooutput.Text);
tofloat.in_raster_or_constant = tb_tminput.Text;
tofloat.out_raster = floatfilepath_tm;
gp.Execute(tofloat, null);
//alt输入数据转化为浮点型
string floatfilepath_alt = System.IO.Path.GetDirectoryName(tb_biooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) + "_float_alt" + System.IO.Path.GetExtension(tb_biooutput.Text);
tofloat.in_raster_or_constant = tb_altinput.Text;
tofloat.out_raster = floatfilepath_alt;
gp.Execute(tofloat, null);
#endregion
//输入数据进行归一化处理
#region
//npp数据归一化处理
string npp;
npp = System.IO.Path.GetDirectoryName(tb_biooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) + "_npp" + System.IO.Path.GetExtension(tb_biooutput.Text);
rc.expression = "(\"" + floatfilepath_npp + "\" - " + GetMinPixelValue(floatfilepath_npp).ToString() + ") / (" + (GetMaxPixelValue(floatfilepath_npp) - GetMinPixelValue(floatfilepath_npp)).ToString() + ")";
rc.output_raster = npp;
tb_state.Text = "正在对生态系统净初级生产力平均值进行归一化……";
gp.Execute(rc, null);
tb_state.Text = "生态系统净初级生产力平均值归一化处理完成！";
//pre数据归一化处理
string pre;
pre = System.IO.Path.GetDirectoryName(tb_biooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) + "_pre" + System.IO.Path.GetExtension(tb_biooutput.Text);
rc.expression = "(\"" + floatfilepath_pre + "\" - " + GetMinPixelValue(floatfilepath_pre).ToString() + ") / (" + (GetMaxPixelValue(floatfilepath_pre) - GetMinPixelValue(floatfilepath_pre)).ToString() + ")";
rc.output_raster = pre;
tb_state.Text = "正在对平均年降水量数据进行归一化……";
gp.Execute(rc, null);
tb_state.Text = "平均年降水量数据归一化处理完成！";
//tm数据归一化处理
string tm;
tm = System.IO.Path.GetDirectoryName(tb_biooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) + "_tm" + System.IO.Path.GetExtension(tb_biooutput.Text);
rc.expression = "(\"" + floatfilepath_tm + "\" - " + GetMinPixelValue(floatfilepath_tm).ToString() + ") / (" + (GetMaxPixelValue(floatfilepath_tm) - GetMinPixelValue(floatfilepath_tm)).ToString() + ")";
rc.output_raster = tm;
tb_state.Text = "正在对平均温度进行归一化……";
gp.Execute(rc, null);
tb_state.Text = "平均温度归一化处理完成！";
//alt数据归一化处理
string alt;
alt = System.IO.Path.GetDirectoryName(tb_biooutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) + "_alt" + System.IO.Path.GetExtension(tb_biooutput.Text);
rc.expression = "(\"" + floatfilepath_alt + "\" - " + GetMinPixelValue(floatfilepath_alt).ToString() + ") / (" + (GetMaxPixelValue(floatfilepath_alt) - GetMinPixelValue(floatfilepath_alt)).ToString() + ")";
rc.output_raster = alt;
tb_state.Text = "正在对海拔参数数据进行归一化……";
gp.Execute(rc, null);
tb_state.Text = "海拔参数数据归一化处理完成！";
#endregion
//计算生物多样性保护服务能力指数
rc.expression = "\"" + npp + "\" * \"" + pre + "\" * \"" + tm + "\" * ( 1 - \"" + alt + "\" )";
rc.output_raster = tb_biooutput.Text;
tb_state.Text = "正在计算计算生物多样性保护服务能力指数……";
gp.Execute(rc, null);
tb_state.Text = "计算生物多样性保护服务能力指数计算完成！";
//删除临时文件
if (tb_state.Text == "计算生物多样性保护服务能力指数计算完成！")
{
string[] files = Directory.GetFiles(System.IO.Path.GetDirectoryName(tb_biooutput.Text), System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) + "*");
foreach (string file in files)
{
if (System.IO.Path.GetFileNameWithoutExtension(tb_biooutput.Text) != System.IO.Path.GetFileNameWithoutExtension(file))
{
File.Delete(file);
}
}
}
//将结果加载显示
#region
string mxfile = System.Environment.CurrentDirectory + "\\Data\\EcoRedLine\\生物多样性生态红线划分.mxd";
IMapDocument pMapDocument = new MapDocumentClass();
pMapDocument.Open(mxfile); //打开本地的地图文档，用来操作改mxd文件
IMap pMap = pMapDocument.get_Map(0);
IMapLayers pMapLayer = pMap as IMapLayers;
IRasterLayer pRasterLayer = new RasterLayerClass();
IWorkspaceFactory rasterWorkspaceFactory = new RasterWorkspaceFactoryClass();
IRasterWorkspace rasterWorkspace = (IRasterWorkspace)rasterWorkspaceFactory.OpenFromFile(System.IO.Path.GetDirectoryName(tb_biooutput.Text), 0);
IRasterDataset pRasterDataset1 = rasterWorkspace.OpenRasterDataset("Biodiversity.tif"); //打开栅格图的文件名
pRasterLayer.CreateFromDataset(pRasterDataset1); //创建
pRasterLayer = pMapLayer.get_Layer(0) as IRasterLayer;
pRasterLayer.CreateFromDataset(pRasterDataset1);
pMapDocument.Save(true, true);//保存更改完路径后的mxd文件
_PageLayoutControl.LoadMxFile(mxfile);
_PageLayoutControl.Extent = _PageLayoutControl.FullExtent;
_PageLayoutControl.ZoomToWholePage();
#endregion
}
catch (Exception ex)
{
MessageBox.Show(gp.GetMessages(ref sev), "提示", MessageBoxButtons.OK);
}
}
private void btOk_Click(object sender, EventArgs e)
{
int sExcute = 0;
int sIsNul = 0;
int sIsThan = 0;
double sAll = 0;
for (int i = 0; i < this.dataGridView1.Rows.Count - 1; i++)
{
if (this.dataGridView1.Rows[i].Cells[1].Value == null)
{
sIsNul = 1; //如果权重存在空值，sIsNul=1
}
else
{
double ssRow = double.Parse(this.dataGridView1.Rows[i].Cells[1].Value.ToString());
sAll = sAll + ssRow;
}
}
if ((sAll != 1.0) && (sIsNul == 0))
{
sIsThan = 1; //如果权重之和不为1，sIsThan=1
}
//判断权重是否为空
if (sIsNul == 1)
{
MessageBox.Show("请您输入数据权重！");
}
//判断权重之和是否等于1
if (sIsThan == 1)
{
MessageBox.Show("请您确保输入的权重之和等于1");
}
if (txtSavePath.Equals(""))
{
MessageBox.Show("请选择输出路径！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
//准备后续执行
if ((sIsNul == 0) && (sIsThan == 0) && (Onebt == 1))
{
//现创建一个存放临时文件的临时文件夹
string newPath = System.IO.Path.Combine(Temfile, "");
System.IO.Directory.CreateDirectory(newPath);
this.rtxtState.AppendText("正在执行，请您耐心等待...\n");
this.rtxtState.ScrollToCaret();
this.rtxtState.AppendText("准备调用GP工具箱...\n");
this.rtxtState.ScrollToCaret();
IVariantArray parameters = new VarArrayClass();
Geoprocessor GP = new Geoprocessor();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCalR = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
this.rtxtState.AppendText("开始遍历读取表中数据...\n");
this.rtxtState.ScrollToCaret();
this.rtxtState.AppendText("开始对栅格数据重分类...\n");
this.rtxtState.ScrollToCaret();
for (int m = 0; m < this.dataGridView1.Rows.Count - 1; m++)
{
string sFileName;
sFileName = Temfile + "\\重分类" + (m + 1).ToString() + ".tif";
//计算栅格最小值
double dMin0 = 0;
IGeoProcessor2 gp = new GeoProcessorClass();
gp.OverwriteOutput = true;
// Create a variant array to hold the parameter values.
IVariantArray parameters2 = new VarArrayClass();
IGeoProcessorResult result = new GeoProcessorResultClass();
// Set parameter values.
parameters2.Add(this.dataGridView1.Rows[m].Cells[0].Value.ToString());
parameters2.Add("MINIMUM");
result = gp.Execute("GetRasterProperties_management", parameters2, null);
dMin0 = (double)result.ReturnValue;
int dMin = (int)dMin0;
//计算栅格最大值
double dMax0 = 0;
IGeoProcessor2 gp2 = new GeoProcessorClass();
gp2.OverwriteOutput = true;
// Create a variant array to hold the parameter values.
IVariantArray parameters3 = new VarArrayClass();
IGeoProcessorResult result3 = new GeoProcessorResultClass();
// Set parameter values.
parameters3.Add(this.dataGridView1.Rows[m].Cells[0].Value.ToString());
parameters3.Add("MAXIMUM");
result3 = gp2.Execute("GetRasterProperties_management", parameters3, null);
dMax0 = (double)result3.ReturnValue;
int dMax = (int)dMax0;
//计算分类区间
int Avera = (dMax - dMin) / 4;
int interval1 = dMin + Avera;
int interval2 = dMin + 2 * Avera;
int interval3 = dMin + 3 * Avera;
//开始对数据进行重分类
ESRI.ArcGIS.SpatialAnalystTools.Reclassify tReclass = new ESRI.ArcGIS.SpatialAnalystTools.Reclassify();
tReclass.in_raster = this.dataGridView1.Rows[m].Cells[0].Value.ToString();
tReclass.missing_values = "NODATA";
tReclass.out_raster = sFileName;
tReclass.reclass_field = "VALUE";
//tReclass.remap = "1400 2176 1;2176 2538 2;2538 3040 3;3040 4073 4";
tReclass.remap = dMin.ToString() + " " + interval1.ToString() + " " + "1;" + interval1.ToString() + " " + interval2.ToString() + " " + "2;" + interval2.ToString() + " " + interval3.ToString() + " " + "3;" + interval3.ToString() + " " + dMax.ToString() + " " + "4;";
// ScrollToBottom("Reclassify");
//tGeoResult = (IGeoProcessorResult)tGp.Execute(tReclass, null);
GP.Execute(tReclass, null);
}
this.rtxtState.AppendText("输入数据重分类完成...\n");
this.rtxtState.ScrollToCaret();
//开始进行栅格计算
this.rtxtState.AppendText("开始准备进行栅格加权运算...\n");
this.rtxtState.ScrollToCaret();
string sFileName2 = "重分类";
string sRoad = Temfile + "\\重分类";
for (int n = 1; n < this.dataGridView1.Rows.Count; n++)
{
sFileName2 = sRoad + n.ToString() + "." + "tif";
strExp = "\"" + sFileName2 + "\"" + "*" + double.Parse(this.dataGridView1.Rows[n - 1].Cells[1].Value.ToString());
if (n < this.dataGridView1.Rows.Count - 1)
{
sR2 = sR2 + strExp + "+";
}
else
{
sR2 = sR2 + strExp;
}
}
sR3 = sR2;
sCalR.expression = sR3;
sR = Temfile + "\\地质灾害.tif";
sCalR.output_raster = sR;
GP.Execute(sCalR, null);
this.rtxtState.AppendText("栅格计算成功，得到地质灾害分布栅格影像...\n");
this.rtxtState.ScrollToCaret();
this.rtxtState.AppendText("开始对地质灾害影像进行重分类...\n");
this.rtxtState.ScrollToCaret();
//开始准备对生成的地质灾害.tif进行分类
//计算栅格最小值
double dMin02 = 0;
IGeoProcessor2 gp3 = new GeoProcessorClass();
gp3.OverwriteOutput = true;
// Create a variant array to hold the parameter values.
IVariantArray parameters22 = new VarArrayClass();
IGeoProcessorResult result22 = new GeoProcessorResultClass();
// Set parameter values.
parameters22.Add(sR);
parameters22.Add("MINIMUM");
result22 = gp3.Execute("GetRasterProperties_management", parameters22, null);
dMin02 = (double)result22.ReturnValue;
//计算栅格最大值
double dMax02 = 0;
IGeoProcessor2 gp4 = new GeoProcessorClass();
gp4.OverwriteOutput = true;
// Create a variant array to hold the parameter values.
IVariantArray parameters33 = new VarArrayClass();
IGeoProcessorResult result33 = new GeoProcessorResultClass();
// Set parameter values.
parameters33.Add(sR);
parameters33.Add("MAXIMUM");
result33 = gp4.Execute("GetRasterProperties_management", parameters33, null);
dMax02 = (double)result33.ReturnValue;
//计算分类区间
double Avera2 = (dMax02 - dMin02) / 4;
double interval12 = dMin02 + Avera2;
double interval22 = dMin02 + 2 * Avera2;
double interval32 = dMin02 + 3 * Avera2;
//开始对地质灾害.tif重分类
ESRI.ArcGIS.SpatialAnalystTools.Reclassify tReclass2 = new ESRI.ArcGIS.SpatialAnalystTools.Reclassify();
tReclass2.in_raster = Temfile + "\\地质灾害.tif";
tReclass2.missing_values = "NODATA";
tReclass2.out_raster = txtSavePath.Text;
tReclass2.reclass_field = "VALUE";
//tReclass.remap = "1400 2176 1;2176 2538 2;2538 3040 3;3040 4073 4";
tReclass2.remap = dMin02.ToString() + " " + interval12.ToString() + " " + "1;" + interval12.ToString() + " " + interval22.ToString() + " " + "2;" + interval22.ToString() + " " + interval32.ToString() + " " + "3;" + interval32.ToString() + " " + dMax02.ToString() + " " + "4;";
// ScrollToBottom("Reclassify");
//tGeoResult = (IGeoProcessorResult)tGp.Execute(tReclass, null);
GP.Execute(tReclass2, null);
//删除临时文件夹
string deleteFile = Temfile;
DeleteFolder(deleteFile);
this.rtxtState.AppendText("完成地质灾害生态红线的划分，已将结果成功保存...\n");
this.rtxtState.ScrollToCaret();
}
}
private void bt_ok_Click(object sender, EventArgs e)
{
//判断输入路径是否正确
#region
try
{
if (!File.Exists(tb_raininput.Text))
{
MessageBox.Show("平均降水量数据输入路径不正确或文件被占用，请重新输入路径！", "提示", MessageBoxButtons.OK);
return;
}
if (!File.Exists(tb_suninput.Text))
{
MessageBox.Show("平均日照时长数据输入路径不正确或文件被占用，请重新输入路径！", "提示", MessageBoxButtons.OK);
return;
}
if (!File.Exists(tb_tminput.Text))
{
MessageBox.Show("平均温度数据输入路径不正确或文件被占用，请重新输入路径！", "提示", MessageBoxButtons.OK);
return;
}
if (!File.Exists(tb_luinput.Text))
{
MessageBox.Show("土地利用类型数据输入路径不正确或文件被占用，请重新输入路径！", "提示", MessageBoxButtons.OK);
return;
}
}
catch (Exception ex)
{
MessageBox.Show("设置路径不合法，请检查！");
return;
}
#endregion
//水源涵养生态红线
#region
object sev = null;
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator rc = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
Geoprocessor gp = new Geoprocessor();
gp.OverwriteOutput = true;
try
{
//计算地表净辐射量（地表净辐射量=日照时长/20）
string r;
r = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_r" + System.IO.Path.GetExtension(tb_wcoutput.Text);
rc.expression = "(\"" + tb_suninput.Text + "\")" + " / 20";
rc.output_raster = r;
tb_state.Text = "正在计算地表净辐射量……";
gp.Execute(rc, null);
tb_state.Text = "地表净辐射量计算完成！";
//计算多年平均潜在蒸发量
string pet, pet1, pet2, pet3, pet4, pet5;
pet1 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_pet1" + System.IO.Path.GetExtension(tb_wcoutput.Text);
pet2 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_pet2" + System.IO.Path.GetExtension(tb_wcoutput.Text);
pet3 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_pet3" + System.IO.Path.GetExtension(tb_wcoutput.Text);
pet4 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_pet4" + System.IO.Path.GetExtension(tb_wcoutput.Text);
pet5 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_pet5" + System.IO.Path.GetExtension(tb_wcoutput.Text);
pet = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_pet" + System.IO.Path.GetExtension(tb_wcoutput.Text);
tb_state.Text = "正在计算多年平均潜在蒸发量……";
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 1";
rc.output_raster = pet1;
gp.Execute(rc, null);
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 2";
rc.output_raster = pet2;
gp.Execute(rc, null);
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 3";
rc.output_raster = pet3;
gp.Execute(rc, null);
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 4";
rc.output_raster = pet4;
gp.Execute(rc, null);
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 5";
rc.output_raster = pet5;
gp.Execute(rc, null);
rc.expression = "(\"" + pet1 + "\"* 336.30)" + " + (\"" + pet2 + "\"* 314.70)" + " + (\"" + pet3 + "\"* 384.58)" + " + (\"" + pet4 + "\"* 182.81)" + " + (\"" + pet5 + "\"* 219.87)";
rc.output_raster = pet;
gp.Execute(rc, null);
tb_state.Text = "多年平均潜在蒸发量计算完成！";
//计算下垫面（土地覆盖）影响系数
string w, w1, w2, w3, w4, w5;
w1 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_w1" + System.IO.Path.GetExtension(tb_wcoutput.Text);
w2 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_w2" + System.IO.Path.GetExtension(tb_wcoutput.Text);
w3 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_w3" + System.IO.Path.GetExtension(tb_wcoutput.Text);
w4 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_w4" + System.IO.Path.GetExtension(tb_wcoutput.Text);
w5 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_w5" + System.IO.Path.GetExtension(tb_wcoutput.Text);
w = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_w" + System.IO.Path.GetExtension(tb_wcoutput.Text);
tb_state.Text = "正在计算下垫面（土地覆盖）影响系数……";
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 1";
rc.output_raster = w1;
gp.Execute(rc, null);
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 2";
rc.output_raster = w2;
gp.Execute(rc, null);
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 3";
rc.output_raster = w3;
gp.Execute(rc, null);
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 4";
rc.output_raster = w4;
gp.Execute(rc, null);
rc.expression = "\"" + tb_luinput.Text + "\"" + " == 5";
rc.output_raster = w5;
gp.Execute(rc, null);
rc.expression = "(\"" + pet1 + "\"* 0.5)" + " + (\"" + pet2 + "\"* 0.1)" + " + (\"" + pet3 + "\"* 1.5)" + " + (\"" + pet4 + "\"* 0.1)" + " + (\"" + pet5 + "\"* 0.5)";
rc.output_raster = w;
gp.Execute(rc, null);
tb_state.Text = "下垫面（土地覆盖）影响系数计算完成！";
//计算蒸散量（蒸散量=0.489+0.289*表净辐射量+0.023*平均温度）
string et, et1, et2, et2_t;
et1 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_et1" + System.IO.Path.GetExtension(tb_wcoutput.Text);
et2 = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_et2" + System.IO.Path.GetExtension(tb_wcoutput.Text);
et2_t = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_et2_t" + System.IO.Path.GetExtension(tb_wcoutput.Text);
et = System.IO.Path.GetDirectoryName(tb_wcoutput.Text) + "\\" + System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "_et" + System.IO.Path.GetExtension(tb_wcoutput.Text);
tb_state.Text = "正在计算蒸散量……";
rc.expression = "\"" + tb_raininput.Text + "\" + " + "\"" + w + "\" * \"" + pet + "\"";
rc.output_raster = et1;
gp.Execute(rc, null);
rc.expression = " 1 + " + "(\"" + w + "\" * \"" + pet + "\" / (\"" + tb_raininput.Text + "\") + \"" + tb_raininput.Text + "\" / \"" + pet + "\")";
rc.output_raster = et2_t;
gp.Execute(rc, null);
rc.expression = "Con(\"" + et2_t + "\"" + " != 0 ,\"" + et2_t + "\")";
rc.output_raster = et2;
gp.Execute(rc, null);
rc.expression = "\"" + et1 + "\" / \"" + et2 + "\"";
rc.output_raster = et;
gp.Execute(rc, null);
tb_state.Text = "蒸散量计算完成！";
//计算水源涵养量(水源涵养量=多年平均降水量-蒸散量)
rc.expression = "\"" + tb_raininput.Text + "\" - \"" + et + "\"";
rc.output_raster = tb_wcoutput.Text;
tb_state.Text = "正在计算水源涵养量……";
gp.Execute(rc, null);
tb_state.Text = "水源涵养量计算完成！";
//删除临时文件
if (tb_state.Text == "水源涵养量计算完成！")
{
string[] files = Directory.GetFiles(System.IO.Path.GetDirectoryName(tb_wcoutput.Text), System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) + "*");
foreach (string file in files)
{
if (System.IO.Path.GetFileNameWithoutExtension(tb_wcoutput.Text) != System.IO.Path.GetFileNameWithoutExtension(file))
{
File.Delete(file);
}
}
}
#endregion
//将结果加载显示
#region
string mxfile = System.Environment.CurrentDirectory + "\\Data\\EcoRedLine\\水源涵养生态红线划分.mxd";
IMapDocument pMapDocument = new MapDocumentClass();
pMapDocument.Open(mxfile); //打开本地的地图文档，用来操作改mxd文件
IMap pMap = pMapDocument.get_Map(0);
IMapLayers pMapLayer = pMap as IMapLayers;
IRasterLayer pRasterLayer = new RasterLayerClass();
IWorkspaceFactory rasterWorkspaceFactory = new RasterWorkspaceFactoryClass();
IRasterWorkspace rasterWorkspace = (IRasterWorkspace)rasterWorkspaceFactory.OpenFromFile(System.IO.Path.GetDirectoryName(tb_wcoutput.Text), 0);
IRasterDataset pRasterDataset1 = rasterWorkspace.OpenRasterDataset("waterConervations.tif"); //打开栅格图的文件名
pRasterLayer.CreateFromDataset(pRasterDataset1); //创建
pRasterLayer = pMapLayer.get_Layer(0) as IRasterLayer;
pRasterLayer.CreateFromDataset(pRasterDataset1);
pMapDocument.Save(true, true);//保存更改完路径后的mxd文件
_PageLayoutControl.LoadMxFile(mxfile);
_PageLayoutControl.Extent = _PageLayoutControl.FullExtent;
_PageLayoutControl.ZoomToWholePage();
#endregion
}
catch (Exception ex)
{
MessageBox.Show(gp.GetMessages(ref sev), "提示", MessageBoxButtons.OK);
}
}
private void btOk_Click(object sender, EventArgs e)
{
this.rtxtState.AppendText("正在执行，请您耐心等待...\n");
this.rtxtState.ScrollToCaret();
#region 步骤一：接收输入数据的路径
string sInputFile = this.txtInput.Text;
if (sInputFile.Equals(""))
{
MessageBox.Show("请选择输入DEM数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
this.rtxtState.AppendText("读取DEM数据...\n");
this.rtxtState.ScrollToCaret();
string sOutputFile = this.txtOutput.Text;
if (sOutputFile.Equals(""))
{
MessageBox.Show("请选择输出路径！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
this.rtxtState.AppendText("读取输出路径...\n");
this.rtxtState.ScrollToCaret();
#endregion
//string fileName = "地形地貌.mxd";
//axMapControl1.LoadMxFile(fileName);
////将地图全屏最大化
//axMapControl1.Extent = axMapControl1.FullExtent;
IVariantArray parameters = new VarArrayClass();
Geoprocessor GP = new Geoprocessor();
object sev = null;
this.rtxtState.AppendText("调用GP工具箱功能...\n");
this.rtxtState.ScrollToCaret();
ESRI.ArcGIS.SpatialAnalystTools.Slope slo = new ESRI.ArcGIS.SpatialAnalystTools.Slope();//计算坡度
slo.in_raster = sInputFile;
slo.output_measurement = smeasurement;
slo.z_factor = sfactor;
string sRoad1 = sOutputFile + sslope;
slo.out_raster = sRoad1;
this.rtxtState.AppendText("利用DEM数据计算坡度信息...\n");
this.rtxtState.ScrollToCaret();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCal1 = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCal2 = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCal3 = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCal4 = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
//try
//{
GP.Execute(slo, null);//坡度计算
this.rtxtState.AppendText("计算坡度完成...\n");
this.rtxtState.ScrollToCaret();
//坡度<=8°&高程<=2500m
sOrder1 = "(\"" + sOutputFile + "/slope.tif" + "\"" + "<=8" + ")" + "&" + "(\"" + sInputFile + "\"" + "<=2500" + ")";
sCal1.expression = sOrder1;
sCal1.output_raster = sOutputFile + sSlope08;
this.rtxtState.AppendText("开始进行建设区分级计算...\n");
this.rtxtState.ScrollToCaret();
GP.Execute(sCal1, null);
this.rtxtState.AppendText("计算出“优化建设区”...\n");
this.rtxtState.ScrollToCaret();
//8°<坡度<=15°& 高程<=2500m
sOrder2 = "(" + "8<" + "\"" + sOutputFile + "/slope.tif" + "\")" + "&" + "(\"" + sOutputFile + "/slope.tif" + "\"" + "<=15" + ")" + "&" + "(\"" + sInputFile + "\"" + "<=2500" + ")";
sCal2.expression = sOrder2;
sCal2.output_raster = sOutputFile + sSlope815;
GP.Execute(sCal2, null);
this.rtxtState.AppendText("计算出“允许建设区”...\n");
this.rtxtState.ScrollToCaret();
//15°<坡度<=25°&高程<=2500m
sOrder3 = "(" + "15<" + "\"" + sOutputFile + "/slope.tif" + "\")" + "&" + "(\"" + sOutputFile + "/slope.tif" + "\"" + "<=25" + ")" + "&" + "(\"" + sInputFile + "\"" + "<=2500" + ")";
sCal3.expression = sOrder3;
sCal3.output_raster = sOutputFile + sSlope1525;
GP.Execute(sCal3, null);
this.rtxtState.AppendText("计算出“有条件建设区”...\n");
this.rtxtState.ScrollToCaret();
//25°<坡度 | 高程<=2500m
sOrder4 = "(\"" + sOutputFile + "/slope.tif" + "\"" + ">25" + ")" + "|" + "(\"" + sInputFile + "\"" + ">2500" + ")";
sCal4.expression = sOrder4;
sCal4.output_raster = sOutputFile + sSlope25;
GP.Execute(sCal4, null);
this.rtxtState.AppendText("计算出禁止建设区...\n");
this.rtxtState.ScrollToCaret();
this.rtxtState.AppendText("结果计算完成...\n");
this.rtxtState.ScrollToCaret();
//MessageBox.Show(GP.GetMessages(ref sev));
//修改mxd数据源路径
this.rtxtState.AppendText("开始将结果用地图输出展示，请稍等...\n");
this.rtxtState.ScrollToCaret();
IMapDocument pMapDocument = new MapDocumentClass();
pMapDocument.Open("模板地形地貌.mxd"); //打开本地的“模板地形地貌.mxd”地图文档，用来操作改mxd文件
IMap pMap = pMapDocument.get_Map(0);
IMapLayers pMapLayer = pMap as IMapLayers;
IRasterLayer pRasterLayer = new RasterLayerClass();
IWorkspaceFactory rasterWorkspaceFactory = new RasterWorkspaceFactoryClass();
IRasterWorkspace rasterWorkspace = (IRasterWorkspace)rasterWorkspaceFactory.OpenFromFile(txtOutput.Text, 0);
IRasterDataset pRasterDataset1 = rasterWorkspace.OpenRasterDataset("优化建设区.tif"); //打开"优化建设区.tif"为栅格图的文件名
pRasterLayer.CreateFromDataset(pRasterDataset1); //创建
pRasterLayer = pMapLayer.get_Layer(1) as IRasterLayer;
pRasterLayer.CreateFromDataset(pRasterDataset1);
//更改“允许建设区.tif”数据源路径
IRasterLayer pRasterLayer2 = new RasterLayerClass();
IRasterDataset pRasterDataset2 = rasterWorkspace.OpenRasterDataset("允许建设区.tif"); //打开"允许建设区.tif"为栅格图的文件名
pRasterLayer2.CreateFromDataset(pRasterDataset2); //创建
pRasterLayer2 = pMapLayer.get_Layer(2) as IRasterLayer;
pRasterLayer2.CreateFromDataset(pRasterDataset2);
//更改"有条件建设区.tif"数据源路径
IRasterLayer pRasterLayer3 = new RasterLayerClass();
IRasterDataset pRasterDataset3 = rasterWorkspace.OpenRasterDataset("有条件建设区.tif"); //打开"有条件建设区.tif"为栅格图的文件名
pRasterLayer3.CreateFromDataset(pRasterDataset3); //创建
pRasterLayer3 = pMapLayer.get_Layer(3) as IRasterLayer;
pRasterLayer3.CreateFromDataset(pRasterDataset3);
//更改"禁止建设区.tif"数据源路径
IRasterLayer pRasterLayer4 = new RasterLayerClass();
IRasterDataset pRasterDataset4 = rasterWorkspace.OpenRasterDataset("禁止建设区.tif"); //打开"禁止建设区.tif"为栅格图的文件名
pRasterLayer4.CreateFromDataset(pRasterDataset4); //创建
pRasterLayer4 = pMapLayer.get_Layer(4) as IRasterLayer;
pRasterLayer4.CreateFromDataset(pRasterDataset4);
pMapDocument.Save(true, true);//保存更改完路径后的mxd文件模板地形地貌2
//string fileName = "南湖村.mxd";// @"C:\Users\tangmeng\Desktop\农用地籍数据库\实验四\南湖数据\南湖村.mxd";
//pMap.LoadMxFile(fileName);
////将地图全屏最大化
//pMap.Extent = pMap.FullExtent;
string MdbFile = "模板地形地貌.mxd";
_PageLayoutControl.LoadMxFile(MdbFile);
//将mxd全屏最大化
_PageLayoutControl.Extent = _PageLayoutControl.FullExtent;
// tabControl1.SelectedTab = tabPage1;
this.rtxtState.AppendText("地图输出完成...\n");
this.rtxtState.ScrollToCaret();
}
private void btnOK_Click(object sender, EventArgs e)
{
if (txtSoilClay.Text.Equals(""))
{
MessageBox.Show("请选择输入土壤黏粒含量(%)数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if (txtSoilSlit.Text.Equals(""))
{
MessageBox.Show("请选择输入土壤粉粒含量(%)数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if (txtSoilSand.Text.Equals(""))
{
MessageBox.Show("请选择输入土壤砂粒含量(%)数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if (txtSoilOrganic.Text.Equals(""))
{
MessageBox.Show("请选择输入土壤有机物含量(%)数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if (txtDem.Text.Equals(""))
{
MessageBox.Show("请选择输入DEM数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if (txtC.Text.Equals(""))
{
MessageBox.Show("请选择输入地表覆被因子相关数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if (txtP.Text.Equals(""))
{
MessageBox.Show("请选择输入水土保持措施因子数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if (txtSavePath.Text.Equals(""))
{
MessageBox.Show("请选择结果输出路径！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if (rbtnR.Checked && txtR.Text.Equals(""))
{
MessageBox.Show("请选择输入降雨侵蚀力因子R相关数据！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if ((!rbtnR.Checked) && txtPcpPath.Text.Equals(""))
{
MessageBox.Show("请选择输入多年平均各月降雨量数据所在路径！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if ((!rbtnR.Checked) && txtPcpPrefix.Text.Equals(""))
{
MessageBox.Show("请输入数据前缀(如：pcp_*)！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if ((!rbtnR.Checked) && txtPcpSuffix.Text.Equals(""))
{
MessageBox.Show("请输入数据后缀(如：tif)！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
if ((!rbtnR.Checked) && txtYear.Text.Equals(""))
{
MessageBox.Show("请输入年平均降水量！", "提示", MessageBoxButtons.OK, MessageBoxIcon.Information);
return;
}
//现创建一个存放临时文件的临时文件夹
string newPath = System.IO.Path.Combine(Temfile, "");
System.IO.Directory.CreateDirectory(newPath);
this.rtxtState.AppendText("正在执行，请您耐心等待...\n");
this.rtxtState.ScrollToCaret();
this.rtxtState.AppendText("开始准备配置文件...\n");
this.rtxtState.ScrollToCaret();
IVariantArray parameters = new VarArrayClass();
Geoprocessor GP = new Geoprocessor();
this.rtxtState.AppendText("准备调用GP工具箱...\n");
this.rtxtState.ScrollToCaret();
//ESRI.ArcGIS.DataManagementTools.GetRasterProperties NDVIMin = new ESRI.ArcGIS.DataManagementTools.GetRasterProperties();
//ESRI.ArcGIS.DataManagementTools.GetRasterProperties NDVIMax = new ESRI.ArcGIS.DataManagementTools.GetRasterProperties();
ESRI.ArcGIS.SpatialAnalystTools.Slope slo = new ESRI.ArcGIS.SpatialAnalystTools.Slope();//计算坡度
ESRI.ArcGIS.SpatialAnalystTools.Fill demFill = new ESRI.ArcGIS.SpatialAnalystTools.Fill();
ESRI.ArcGIS.SpatialAnalystTools.FlowDirection Filldec = new ESRI.ArcGIS.SpatialAnalystTools.FlowDirection();
ESRI.ArcGIS.SpatialAnalystTools.FlowAccumulation DecAcc = new ESRI.ArcGIS.SpatialAnalystTools.FlowAccumulation();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCals = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCalm = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCalLS = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCalR = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCalK = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCalA = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator sCalC = new ESRI.ArcGIS.SpatialAnalystTools.RasterCalculator();
//1计算参数R
if (rbtnR.Checked)
{
sR = txtR.Text;
this.rtxtState.AppendText("降水侵蚀因子R读取成功，准备计算参数K...\n");
this.rtxtState.ScrollToCaret();
}
else
{
this.rtxtState.AppendText("开始计算参数R...\n");
this.rtxtState.ScrollToCaret();
string sFileName = "";
string sYear = txtYear.Text;
for (int i = 1; i < 13; i++)
{
sFileName = txtPcpPath.Text + "\\" + this.txtPcpPrefix.Text + i.ToString() + "." + this.txtPcpSuffix.Text;
strExp = "(1.735 * Power(10,1.5 * Log10((" + "\"" + sFileName + "\"" + " * " + "\"" + sFileName + "\"" + ") /" + "\"" + txtYear.Text + "\")" + " - 0.08188))";
if (i < 12)
{
sR2 = sR2 + strExp + "+";
}
else
{
sR2 = sR2 + strExp;
}
}
sR3 = sR2;
sCalR.expression = sR3;
sR = Temfile + "\\CalR.tif";
sCalR.output_raster = sR;
GP.Execute(sCalR, null);
this.rtxtState.AppendText("降水侵蚀因子R计算成功，准备计算参数K...\n");
this.rtxtState.ScrollToCaret();
}
//2计算参数K
sK = "(0.2 + 0.3 * Exp(-0.0256 *" + "\"" + txtSoilSand.Text + "\"" + "* (1.0 - " + "\"" + txtSoilSlit.Text + "\"" + " / 100.0))) * Power((" + "\"" + txtSoilSlit.Text + "\"" + " * 1.0 / (" + "\"" + txtSoilClay.Text + "\"" + " * 1.0 + " + "\"" + txtSoilSlit.Text + "\"" + " * 1.0)), 0.3) * (1.0 - 0.25 * " + "\"" + txtSoilOrganic.Text + "\"" + " * 0.58 / (" + "\"" + txtSoilOrganic.Text + "\"" + " * 0.58 + Exp(3.72 - 2.95 * " + "\"" + txtSoilOrganic.Text + "\"" + " * 0.58))) * (1.0 - (0.7 * (1.0 - " + "\"" + txtSoilSand.Text + "\"" + " / 100.0)) / ((1.0 - " + "\"" + txtSoilSand.Text + "\"" + " / 100.0) + Exp(-5.51 + 22.9 * (1.0 - " + "\"" + txtSoilSand.Text + "\"" + " / 100.0))))";
sCalK.expression = sK;
CalKpath = Temfile + "\\CalK.tif";
sCalK.output_raster = CalKpath;
GP.Execute(sCalK, null);
this.rtxtState.AppendText("土壤可蚀性因子K计算成功...\n");
this.rtxtState.ScrollToCaret();
this.rtxtState.AppendText("准备计算地形因子LS...\n");
this.rtxtState.ScrollToCaret();
//3计算参数LS
//Fill Dem
this.rtxtState.AppendText("开始填充洼地...\n");
this.rtxtState.ScrollToCaret();
demFill.in_surface_raster = txtDem.Text;
demFill.out_surface_raster = Temfile + "\\demFill.tif";
GP.Execute(demFill, null);
// cal FlowDirection
this.rtxtState.AppendText("开始计算流向...\n");
this.rtxtState.ScrollToCaret();
Filldec.in_surface_raster = Temfile + "\\demFill.tif";
Filldec.out_flow_direction_raster = Temfile + "\\FillDec.tif";
GP.Execute(Filldec, null);
//cal FlowAccumulation
this.rtxtState.AppendText("开始计算流量...\n");
this.rtxtState.ScrollToCaret();
Filldec.in_surface_raster = Temfile + "\\FillDec.tif";
FlowAcc = Temfile + "\\FlowAcc.tif";
Filldec.out_flow_direction_raster = FlowAcc;
GP.Execute(Filldec, null);
//先计算坡度
this.rtxtState.AppendText("开始计算坡度...\n");
this.rtxtState.ScrollToCaret();
slo.in_raster = txtDem.Text;
slo.output_measurement = "DEGREE";
slo.z_factor = 1;
string sRoad1 = Temfile + "\\Slope.tif";
slo.out_raster = sRoad1;
GP.Execute(slo, null);//坡度计算
//cal S
sS = "Con(" + "\"" + sRoad1 + "\"" + " < 5,10.8 * Sin(" + "\"" + sRoad1 + "\"" + " * 3.14 / 180) + 0.03,Con(" + "\"" + sRoad1 + "\"" + " >= 10,21.9 * Sin(" + "\"" + sRoad1 + "\"" + " * 3.14 / 180) - 0.96,16.8 * Sin(" + "\"" + sRoad1 + "\"" + " * 3.14 / 180) - 0.5))";
sCals.expression = sS;
Calspath = Temfile + "\\CalS.tif";
sCals.output_raster = Calspath;
GP.Execute(sCals, null);
//cal m
sM = "Con(" + "\"" + sRoad1 + "\"" + " <= 1,0.2,Con(" + "\"" + sRoad1 + "\"" + " <= 3,0.3,Con(" + "\"" + sRoad1 + "\"" + " <= 5,0.4,0.5)))";
sCalm.expression = sM;
Calmpath = Temfile + "\\CalM.tif";
sCalm.output_raster = Calmpath;
GP.Execute(sCalm, null);
//cal ls
sLs = "\"" + Calspath + "\"" + " * Power((" + "\"" + CalKpath + "\"" + " * " + this.txtCellSize.Text + " / 22.1)," + "\"" + Calmpath + "\"" + ")";
sCalLS.expression = sLs;
Callspath = Temfile + "\\CalLS.tif";
sCalm.output_raster = Calmpath;
sCalLS.output_raster = Callspath;
GP.Execute(sCalLS, null);
this.rtxtState.AppendText("地形因子LS计算成功...\n");
this.rtxtState.ScrollToCaret();
//4计算参数C
if (rbtnVegCover.Checked)
{
CalCpath = txtC.Text;
this.rtxtState.AppendText("地表覆盖因子C读取成功...\n");
this.rtxtState.ScrollToCaret();
}
else
{
this.rtxtState.AppendText("准备计算地表覆盖因子C...\n");
this.rtxtState.ScrollToCaret();
//计算NDVI最小值
//NDVIMin.in_raster = txtC.Text;
//CalKpath = txtDem.Text + "/CalK.tif";
//NDVIMin.property_type = "MINIMUM";
//GP.Execute(sCalK, null);
double dMin = 0;
IGeoProcessor2 gp = new GeoProcessorClass();
gp.OverwriteOutput = true;
// Create a variant array to hold the parameter values.
IVariantArray parameters2 = new VarArrayClass();
IGeoProcessorResult result = new GeoProcessorResultClass();
// Set parameter values.
parameters2.Add(txtC.Text);
parameters2.Add("MINIMUM");
result = gp.Execute("GetRasterProperties_management", parameters2, null);
dMin = (double)result.ReturnValue;
//计算NDVI最大值
double dMax = 0;
IGeoProcessor2 gp2 = new GeoProcessorClass();
gp2.OverwriteOutput = true;
// Create a variant array to hold the parameter values.
IVariantArray parameters3 = new VarArrayClass();
IGeoProcessorResult result3 = new GeoProcessorResultClass();
// Set parameter values.
parameters3.Add(txtC.Text);
parameters3.Add("MAXIMUM");
result3 = gp2.Execute("GetRasterProperties_management", parameters3, null);
dMax = (double)result3.ReturnValue;
//最后计算C
sC = "(" + "\"" + txtC.Text + "\"" + " - " + dMin + ") / (" + dMax + " - " + dMin + ")";
sCalC.expression = sC;
CalCpath = Temfile + "\\CalC.tif";
sCalC.output_raster = CalCpath;
GP.Execute(sCalC, null);
this.rtxtState.AppendText("地表覆盖因子C计算成功...\n");
this.rtxtState.ScrollToCaret();
}
//5计算P
sP = txtP.Text;
this.rtxtState.AppendText("读取水土保持措施因子P...\n");
this.rtxtState.ScrollToCaret();
//最后开始计算A=R*K*LS*C*P
this.rtxtState.AppendText("准备计算水土流失方程...\n");
this.rtxtState.ScrollToCaret();
sA = "\"" + sR + "\"" + " * " + "\"" + CalKpath + "\"" + " * " + "\"" + Callspath + "\"" + " * (1 - " + "\"" + CalCpath + "\"" + ") * " + "\"" + sP + "\"";
sCalA.expression = sA;
sCalA.output_raster = txtSavePath.Text;
GP.Execute(sCalA, null);
//删除临时文件夹
string deleteFile = Temfile;
DeleteFolder(deleteFile);
this.rtxtState.AppendText("计算成功，已将结果成功保存...\n");
this.rtxtState.ScrollToCaret();
}
