//解决路径
public void SolvePath(string WeightName)
{
try
{
int intEdgeUserClassID;
int intEdgeUserID;
int intEdgeUserSubID;
int intEdgeID;
IPoint ipFoundEdgePoint;
double dblEdgePercent;
ITraceFlowSolverGEN ipTraceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
INetSolver ipNetSolver = ipTraceFlowSolver as INetSolver;
INetwork ipNetwork = m_ipGeometricNetwork.Network;
ipNetSolver.SourceNetwork = ipNetwork;
INetElements ipNetElements = ipNetwork as INetElements;
int intCount = m_ipPoints.PointCount;//这里的points有值吗？
//定义一个边线旗数组
IEdgeFlag[] pEdgeFlagList = new EdgeFlagClass[intCount];
for (int i = 0; i < intCount; i++)
{
INetFlag ipNetFlag = new EdgeFlagClass() as INetFlag;
IPoint ipEdgePoint = m_ipPoints.get_Point(i);
//查找输入点的最近的边线
m_ipPointToEID.GetNearestEdge(ipEdgePoint, out intEdgeID, out ipFoundEdgePoint, out dblEdgePercent);
ipNetElements.QueryIDs(intEdgeID, esriElementType.esriETEdge, out intEdgeUserClassID, out intEdgeUserID, out intEdgeUserSubID);
ipNetFlag.UserClassID = intEdgeUserClassID;
ipNetFlag.UserID = intEdgeUserID;
ipNetFlag.UserSubID = intEdgeUserSubID;
IEdgeFlag pTemp = (IEdgeFlag)(ipNetFlag as IEdgeFlag);
pEdgeFlagList[i] = pTemp;
}
ipTraceFlowSolver.PutEdgeOrigins(ref pEdgeFlagList);
INetSchema ipNetSchema = ipNetwork as INetSchema;
INetWeight ipNetWeight = ipNetSchema.get_WeightByName(WeightName);
INetSolverWeights ipNetSolverWeights = ipTraceFlowSolver as INetSolverWeights;
ipNetSolverWeights.FromToEdgeWeight = ipNetWeight; //开始边线的权重
ipNetSolverWeights.ToFromEdgeWeight = ipNetWeight; //终止边线的权重
object[] vaRes = new object[intCount - 1];
//通过findpath得到边线和交汇点的集合
ipTraceFlowSolver.FindPath(esriFlowMethod.esriFMConnected,
esriShortestPathObjFn.esriSPObjFnMinSum,
out m_ipEnumNetEID_Junctions, out m_ipEnumNetEID_Edges, intCount - 1, ref vaRes);
//计算元素成本
m_dblPathCost = 0;
for (int i = 0; i < vaRes.Length; i++)
{
double m_Va = (double)vaRes[i];//我修改过
m_dblPathCost = m_dblPathCost + m_Va;
}
m_ipPolyline = null;
}
catch (Exception ex)
{
Console.WriteLine(ex.Message);
}
}
public int [] GetRouteFromEdgeFlags(ArrayList edgeFlags)
{
IEdgeFlag [] edgeflags = new IEdgeFlag[edgeFlags.Count];
for (int i = 0; i < edgeFlags.Count; i++)
{
edgeflags[i] = edgeFlags[i] as IEdgeFlag;
}
double length = -1;
return(FindPath(edgeflags, out length));
}
public void SolvePath(string WeightName)
{
int num;
ITraceFlowSolver solver = (ITraceFlowSolver) new TraceFlowSolverClass();
INetSolver solver2 = (INetSolver)solver;
INetwork network = this.m_ipGeometricNetwork.Network;
solver2.SourceNetwork = network;
INetElements elements = (INetElements)network;
IEdgeFlag[] edgeOrigins = new IEdgeFlag[this.m_ipPoints.PointCount + 1];
for (num = 0; num < this.m_ipPoints.PointCount; num++)
{
IPoint point;
int num2;
double num3;
int num4;
int num5;
int num6;
INetFlag flag = new EdgeFlagClass();
IPoint inputPoint = this.m_ipPoints.get_Point(num);
this.m_ipPointToEID.GetNearestEdge(inputPoint, out num2, out point, out num3);
elements.QueryIDs(num2, esriElementType.esriETEdge, out num4, out num5, out num6);
flag.UserClassID = num4;
flag.UserID = num5;
flag.UserSubID = num6;
edgeOrigins[num] = (IEdgeFlag)flag;
}
(solver as ITraceFlowSolverGEN).PutEdgeOrigins(ref edgeOrigins);
try
{
INetSolverWeights weights = (INetSolverWeights)solver;
weights.JunctionWeight = this.JunctionWeight;
weights.FromToEdgeWeight = this.FromToEdgeWeight;
weights.ToFromEdgeWeight = this.ToFromEdgeWeight;
}
catch
{
}
object[] segmentCosts = new object[this.m_ipPoints.PointCount];
(solver as ITraceFlowSolverGEN).FindPath(esriFlowMethod.esriFMConnected,
esriShortestPathObjFn.esriSPObjFnMinSum, out this.m_ipEnumNetEID_Junctions,
out this.m_ipEnumNetEID_Edges, this.m_ipPoints.PointCount - 1, ref segmentCosts);
this.m_dblPathCost = 0.0;
for (num = 0; num < segmentCosts.Length; num++)
{
this.m_dblPathCost += (double)segmentCosts[num];
}
this.m_ipPolyline = null;
}
private void ResultBtn_Click(object sender, EventArgs e)
{
try
{
#region 设置管点分析分析条件
//为追踪任务分析器设置管点
IJunctionFlag[] arrayJunctionFlag = new IJunctionFlag[listJunctionsFlag.Count];
for (int i = 0; i < arrayJunctionFlag.Length; i++)
{
arrayJunctionFlag[i] = listJunctionsFlag[i];
}
traceFlowsolverGen.PutJunctionOrigins(ref arrayJunctionFlag);
//为追踪任务分析器设置管线
IEdgeFlag[] arrayEdgeFlag = new IEdgeFlag[listEdgeFlag.Count];
for (int i = 0; i < arrayEdgeFlag.Length; i++)
{
traceFlowsolverGen.PutEdgeOrigins(ref arrayEdgeFlag);
}
//为管点分析器设置障碍点
INetElementBarriersGEN netElementBarriersGEN = new NetElementBarriersClass();
netElementBarriersGEN.Network = m_GeometryNetwork.Network;
if (JunctionBarrierEIDs.Count > 0)
{
int[] junctionBarrierEIDs = new int[JunctionBarrierEIDs.Count];
for (int j = 0; j < junctionBarrierEIDs.Length; j++)
{
junctionBarrierEIDs[j] = JunctionBarrierEIDs[j];
}
netElementBarriersGEN.ElementType = esriElementType.esriETJunction;
netElementBarriersGEN.SetBarriersByEID(ref junctionBarrierEIDs);
netsolver.set_ElementBarriers(esriElementType.esriETJunction, netElementBarriersGEN as INetElementBarriers);
}
else //否则将管点设置为空
{
netsolver.set_ElementBarriers(esriElementType.esriETJunction, null);
}
//未管点分析器设置障碍线
if (EdgeBarrierEIDs.Count > 0)
{
int[] edgeBarrierEIDs = new int[EdgeBarrierEIDs.Count];
for (int j = 0; j < EdgeBarrierEIDs.Count; j++)
{
edgeBarrierEIDs[j] = EdgeBarrierEIDs[j];
}
netElementBarriersGEN.ElementType = esriElementType.esriETEdge;
netElementBarriersGEN.SetBarriersByEID(ref edgeBarrierEIDs);
netsolver.set_ElementBarriers(esriElementType.esriETEdge, netElementBarriersGEN as INetElementBarriers);
}
else //否则将管线设置为空
{
netsolver.set_ElementBarriers(esriElementType.esriETEdge, null);
}
#endregion
//定义相关变量
IEnumNetEID junctionEIDs = new EnumNetEIDArrayClass();
IEnumNetEID edgeEIDs = new EnumNetEIDArrayClass();
object[] segmentsCosts = null;
object totalCost = null;
int Counts = -1;
#region 各种几何网络分析的结果
switch (this.AnalysisCategoryCmbx.SelectedIndex)
{
//查询共同祖先
case 0:
traceFlowsolverGen.FindCommonAncestors(esriFlowElements.esriFEJunctionsAndEdges, out junctionEIDs, out edgeEIDs);
break;
//查找相连接的网络要素
case 1:
traceFlowsolverGen.FindFlowElements(esriFlowMethod.esriFMConnected, esriFlowElements.esriFEJunctionsAndEdges, out junctionEIDs, out edgeEIDs);
break;
//查找网络中连接的环
case 2:
traceFlowsolverGen.FindCircuits(esriFlowElements.esriFEJunctionsAndEdges, out junctionEIDs, out edgeEIDs);
break;
//查找未连接的网络要素
case 3:
traceFlowsolverGen.FindFlowUnreachedElements(esriFlowMethod.esriFMConnected, esriFlowElements.esriFEJunctionsAndEdges, out junctionEIDs, out edgeEIDs);
break;
//查找上游路径，同时跟踪耗费
case 4:
Counts = GetSegmentCounts();
segmentsCosts = new object[Counts];
traceFlowsolverGen.FindSource(esriFlowMethod.esriFMUpstream, esriShortestPathObjFn.esriSPObjFnMinSum, out junctionEIDs, out edgeEIDs, Counts, ref segmentsCosts);
break;
//获取网络路径，追踪网络耗费。count比所有的的标识总数少一个，但如果管点和管线标识同时存在的时候该功能不可用
case 5:
if (listEdgeFlag.Count > 0 && listJunctionsFlag.Count > 0)
{
break;
}
else if (listJunctionsFlag.Count > 0)
{
Counts = listJunctionsFlag.Count - 1;
}
else if (listEdgeFlag.Count > 0)
{
Counts = listEdgeFlag.Count - 1;
}
else
{
break;
}
segmentsCosts = new object[Counts];
traceFlowsolverGen.FindPath(esriFlowMethod.esriFMConnected, esriShortestPathObjFn.esriSPObjFnMinSum, out junctionEIDs, out edgeEIDs, Counts, ref segmentsCosts);
break;
//下游追踪
case 6:
traceFlowsolverGen.FindFlowElements(esriFlowMethod.esriFMDownstream, esriFlowElements.esriFEJunctionsAndEdges, out junctionEIDs, out edgeEIDs);
break;
//查找上游路径消耗，同时获取网络追踪的耗费
case 7:
totalCost = new object();
traceFlowsolverGen.FindAccumulation(esriFlowMethod.esriFMUpstream, esriFlowElements.esriFEJunctionsAndEdges, out junctionEIDs, out edgeEIDs, out totalCost);
break;
//上游追踪
case 8:
Counts = GetSegmentCounts();
segmentsCosts = new object[Counts];
traceFlowsolverGen.FindSource(esriFlowMethod.esriFMUpstream, esriShortestPathObjFn.esriSPObjFnMinSum, out junctionEIDs, out edgeEIDs, Counts, ref segmentsCosts);
break;
default:
break;
}
#endregion
//绘制结果图像
DrawTraceRsult(junctionEIDs, edgeEIDs);
this.Mapcontrol.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGraphics, null, this.Mapcontrol.ActiveView.Extent);
}
catch (Exception ex)
{
MessageBox.Show(ex.Message);
return;
}
}
public ITraceFlowSolverGEN UTIL_coreTraceSetup()
{
// get the current network's logical network
INetworkAnalysisExt nax = m_utilNetExt as INetworkAnalysisExt;
INetwork net = nax.CurrentNetwork.Network;
// create a new TraceFlowSolver object and
// set the source network for the solve
ITraceFlowSolverGEN tfs = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
INetSolver netSolver = tfs as INetSolver;
netSolver.SourceNetwork = net;
// get the barriers for the network, using the element barriers and
// selection barriers that have been added using the user interface
INetworkAnalysisExtBarriers naxBarriers = m_utilNetExt as INetworkAnalysisExtBarriers;
INetElementBarriers juncElemBarriers, edgeElemBarriers;
naxBarriers.CreateElementBarriers(out juncElemBarriers, out edgeElemBarriers);
ISelectionSetBarriers selSetBarriers;
naxBarriers.CreateSelectionBarriers(out selSetBarriers);
netSolver.set_ElementBarriers(esriElementType.esriETJunction, juncElemBarriers);
netSolver.set_ElementBarriers(esriElementType.esriETEdge, edgeElemBarriers);
netSolver.SelectionSetBarriers = selSetBarriers;
// set up the disabled layers for the network solver
// for each feature layer belonging to this network, determine if it is
// enabled or disabled; if it's disabled, then notify the network solver
for (int i = 0; i < nax.FeatureLayerCount; i++)
{
IFeatureLayer featureLayer = nax.get_FeatureLayer(i);
if (naxBarriers.GetDisabledLayer(featureLayer))
{
netSolver.DisableElementClass(featureLayer.FeatureClass.FeatureClassID);
}
}
INetworkAnalysisExtWeightFilter naxWeightFilter = m_utilNetExt as INetworkAnalysisExtWeightFilter;
INetSolverWeightsGEN netSolverWeights = netSolver as INetSolverWeightsGEN;
INetSchema netSchema = net as INetSchema;
// create the junction weight filter
int juncFilterRangeCount = naxWeightFilter.get_FilterRangeCount(esriElementType.esriETJunction);
if (juncFilterRangeCount > 0)
{
INetWeight netWeight = netSchema.get_WeightByName(naxWeightFilter.JunctionWeightFilterName);
netSolverWeights.JunctionFilterWeight = netWeight;
esriWeightFilterType juncWeightFilterType;
bool juncApplyNotOperator;
naxWeightFilter.GetFilterType(esriElementType.esriETJunction,
out juncWeightFilterType, out juncApplyNotOperator);
netSolverWeights.SetFilterType(esriElementType.esriETJunction,
juncWeightFilterType, juncApplyNotOperator);
object[] juncFromValues = new object[juncFilterRangeCount];
object[] juncToValues = new object[juncFilterRangeCount];
for (int i = 0; i < juncFilterRangeCount; i++)
{
naxWeightFilter.GetFilterRange(esriElementType.esriETJunction, i,
out juncFromValues[i], out juncToValues[i]);
}
netSolverWeights.SetFilterRanges(esriElementType.esriETJunction,
ref juncFromValues, ref juncToValues);
}
// create the edge weight filters
int edgeFilterRangeCount = naxWeightFilter.get_FilterRangeCount(esriElementType.esriETEdge);
if (edgeFilterRangeCount > 0)
{
INetWeight fromToNetWeight = netSchema.get_WeightByName(naxWeightFilter.FromToEdgeWeightFilterName);
netSolverWeights.FromToEdgeFilterWeight = fromToNetWeight;
INetWeight toFromNetWeight = netSchema.get_WeightByName(naxWeightFilter.ToFromEdgeWeightFilterName);
netSolverWeights.ToFromEdgeFilterWeight = toFromNetWeight;
esriWeightFilterType edgeWeightFilterType;
bool edgeApplyNotOperator;
naxWeightFilter.GetFilterType(esriElementType.esriETEdge,
out edgeWeightFilterType, out edgeApplyNotOperator);
netSolverWeights.SetFilterType(esriElementType.esriETEdge,
edgeWeightFilterType, edgeApplyNotOperator);
object[] edgeFromValues = new object[edgeFilterRangeCount];
object[] edgeToValues = new object[edgeFilterRangeCount];
for (int i = 0; i < edgeFilterRangeCount; i++)
{
naxWeightFilter.GetFilterRange(esriElementType.esriETEdge, i,
out edgeFromValues[i], out edgeToValues[i]);
}
netSolverWeights.SetFilterRanges(esriElementType.esriETEdge,
ref edgeFromValues, ref edgeToValues);
}
INetworkAnalysisExtFlags naxFlags = m_utilNetExt as INetworkAnalysisExtFlags;
// assign the edge flags to the network solver
int edgeFlagCount = naxFlags.EdgeFlagCount;
if (edgeFlagCount > 0)
{
IEdgeFlag[] edgeFlags = new IEdgeFlag[edgeFlagCount];
for (int i = 0; i < edgeFlagCount; i++)
{
IEdgeFlagDisplay edgeFlagDisplay = naxFlags.get_EdgeFlag(i);
IFlagDisplay flagDisplay = edgeFlagDisplay as IFlagDisplay;
IEdgeFlag edgeFlag = new EdgeFlagClass();
edgeFlag.Position = Convert.ToSingle(edgeFlagDisplay.Percentage);
INetFlag netFlag = edgeFlag as INetFlag;
netFlag.UserClassID = flagDisplay.FeatureClassID;
netFlag.UserID = flagDisplay.FID;
netFlag.UserSubID = flagDisplay.SubID;
edgeFlags[i] = edgeFlag;
}
tfs.PutEdgeOrigins(ref edgeFlags);
}
// assign the junction flags to the network solver
int juncFlagCount = naxFlags.JunctionFlagCount;
if (juncFlagCount > 0)
{
IJunctionFlag[] juncFlags = new IJunctionFlag[juncFlagCount];
for (int i = 0; i < juncFlagCount; i++)
{
IJunctionFlagDisplay juncFlagDisplay = naxFlags.get_JunctionFlag(i);
IFlagDisplay flagDisplay = juncFlagDisplay as IFlagDisplay;
IJunctionFlag juncFlag = new JunctionFlagClass();
INetFlag netFlag = juncFlag as INetFlag;
netFlag.UserClassID = flagDisplay.FeatureClassID;
netFlag.UserID = flagDisplay.FID;
netFlag.UserSubID = flagDisplay.SubID;
juncFlags[i] = juncFlag;
}
tfs.PutJunctionOrigins(ref juncFlags);
}
// set the option for tracing on indeterminate flow
ITraceTasks traceTasks = m_utilNetExt as ITraceTasks;
tfs.TraceIndeterminateFlow = traceTasks.TraceIndeterminateFlow;
return(tfs);
}
public ITraceFlowSolverGEN UTIL_coreTraceSetup()
{
// get the current network's logical network
INetworkAnalysisExt nax = m_utilNetExt as INetworkAnalysisExt;
INetwork net = nax.CurrentNetwork.Network;
// create a new TraceFlowSolver object and
// set the source network for the solve
ITraceFlowSolverGEN tfs = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
INetSolver netSolver = tfs as INetSolver;
netSolver.SourceNetwork = net;
// get the barriers for the network, using the element barriers and
// selection barriers that have been added using the user interface
INetworkAnalysisExtBarriers naxBarriers = m_utilNetExt as INetworkAnalysisExtBarriers;
INetElementBarriers juncElemBarriers, edgeElemBarriers;
naxBarriers.CreateElementBarriers(out juncElemBarriers, out edgeElemBarriers);
ISelectionSetBarriers selSetBarriers;
naxBarriers.CreateSelectionBarriers(out selSetBarriers);
netSolver.set_ElementBarriers(esriElementType.esriETJunction, juncElemBarriers);
netSolver.set_ElementBarriers(esriElementType.esriETEdge, edgeElemBarriers);
netSolver.SelectionSetBarriers = selSetBarriers;
// set up the disabled layers for the network solver
// for each feature layer belonging to this network, determine if it is
// enabled or disabled; if it's disabled, then notify the network solver
for (int i = 0; i < nax.FeatureLayerCount; i++)
{
IFeatureLayer featureLayer = nax.get_FeatureLayer(i);
if (naxBarriers.GetDisabledLayer(featureLayer))
netSolver.DisableElementClass(featureLayer.FeatureClass.FeatureClassID);
}
INetworkAnalysisExtWeightFilter naxWeightFilter = m_utilNetExt as INetworkAnalysisExtWeightFilter;
INetSolverWeightsGEN netSolverWeights = netSolver as INetSolverWeightsGEN;
INetSchema netSchema = net as INetSchema;
// create the junction weight filter
int juncFilterRangeCount = naxWeightFilter.get_FilterRangeCount(esriElementType.esriETJunction);
if (juncFilterRangeCount > 0)
{
INetWeight netWeight = netSchema.get_WeightByName(naxWeightFilter.JunctionWeightFilterName);
netSolverWeights.JunctionFilterWeight = netWeight;
esriWeightFilterType juncWeightFilterType;
bool juncApplyNotOperator;
naxWeightFilter.GetFilterType(esriElementType.esriETJunction,
out juncWeightFilterType, out juncApplyNotOperator);
netSolverWeights.SetFilterType(esriElementType.esriETJunction,
juncWeightFilterType, juncApplyNotOperator);
object[] juncFromValues = new object[juncFilterRangeCount];
object[] juncToValues = new object[juncFilterRangeCount];
for (int i = 0; i < juncFilterRangeCount; i++)
{
naxWeightFilter.GetFilterRange(esriElementType.esriETJunction, i,
out juncFromValues[i], out juncToValues[i]);
}
netSolverWeights.SetFilterRanges(esriElementType.esriETJunction,
ref juncFromValues, ref juncToValues);
}
// create the edge weight filters
int edgeFilterRangeCount = naxWeightFilter.get_FilterRangeCount(esriElementType.esriETEdge);
if (edgeFilterRangeCount > 0)
{
INetWeight fromToNetWeight = netSchema.get_WeightByName(naxWeightFilter.FromToEdgeWeightFilterName);
netSolverWeights.FromToEdgeFilterWeight = fromToNetWeight;
INetWeight toFromNetWeight = netSchema.get_WeightByName(naxWeightFilter.ToFromEdgeWeightFilterName);
netSolverWeights.ToFromEdgeFilterWeight = toFromNetWeight;
esriWeightFilterType edgeWeightFilterType;
bool edgeApplyNotOperator;
naxWeightFilter.GetFilterType(esriElementType.esriETEdge,
out edgeWeightFilterType, out edgeApplyNotOperator);
netSolverWeights.SetFilterType(esriElementType.esriETEdge,
edgeWeightFilterType, edgeApplyNotOperator);
object[] edgeFromValues = new object[edgeFilterRangeCount];
object[] edgeToValues = new object[edgeFilterRangeCount];
for (int i = 0; i < edgeFilterRangeCount; i++)
{
naxWeightFilter.GetFilterRange(esriElementType.esriETEdge, i,
out edgeFromValues[i], out edgeToValues[i]);
}
netSolverWeights.SetFilterRanges(esriElementType.esriETEdge,
ref edgeFromValues, ref edgeToValues);
}
INetworkAnalysisExtFlags naxFlags = m_utilNetExt as INetworkAnalysisExtFlags;
// assign the edge flags to the network solver
int edgeFlagCount = naxFlags.EdgeFlagCount;
if (edgeFlagCount > 0)
{
IEdgeFlag[] edgeFlags = new IEdgeFlag[edgeFlagCount];
for (int i = 0; i < edgeFlagCount; i++)
{
IEdgeFlagDisplay edgeFlagDisplay = naxFlags.get_EdgeFlag(i);
IFlagDisplay flagDisplay = edgeFlagDisplay as IFlagDisplay;
IEdgeFlag edgeFlag = new EdgeFlagClass();
edgeFlag.Position = Convert.ToSingle(edgeFlagDisplay.Percentage);
INetFlag netFlag = edgeFlag as INetFlag;
netFlag.UserClassID = flagDisplay.FeatureClassID;
netFlag.UserID = flagDisplay.FID;
netFlag.UserSubID = flagDisplay.SubID;
edgeFlags[i] = edgeFlag;
}
tfs.PutEdgeOrigins(ref edgeFlags);
}
// assign the junction flags to the network solver
int juncFlagCount = naxFlags.JunctionFlagCount;
if (juncFlagCount > 0)
{
IJunctionFlag[] juncFlags = new IJunctionFlag[juncFlagCount];
for (int i = 0; i < juncFlagCount; i++)
{
IJunctionFlagDisplay juncFlagDisplay = naxFlags.get_JunctionFlag(i);
IFlagDisplay flagDisplay = juncFlagDisplay as IFlagDisplay;
IJunctionFlag juncFlag = new JunctionFlagClass();
INetFlag netFlag = juncFlag as INetFlag;
netFlag.UserClassID = flagDisplay.FeatureClassID;
netFlag.UserID = flagDisplay.FID;
netFlag.UserSubID = flagDisplay.SubID;
juncFlags[i] = juncFlag;
}
tfs.PutJunctionOrigins(ref juncFlags);
}
// set the option for tracing on indeterminate flow
ITraceTasks traceTasks = m_utilNetExt as ITraceTasks;
tfs.TraceIndeterminateFlow = traceTasks.TraceIndeterminateFlow;
return tfs;
}
private void SolvePath(string weightName)
{
//创建ITraceFlowSolverGEN
ITraceFlowSolverGEN pTraceFlowSolverGEN = new TraceFlowSolverClass();
INetSolver pNetSolver = pTraceFlowSolverGEN as INetSolver;
//初始化用于路径计算的Network
INetwork pNetWork = mGeometricNetwork.Network;
pNetSolver.SourceNetwork = pNetWork;
//获取分析经过的点的个数
int intCount = mPointCollectionNet.PointCount;
if (intCount < 1)
{
return;
}
INetFlag pNetFlag;
//用于存储路径计算得到的边
IEdgeFlag[] pEdgeFlags = new IEdgeFlag[intCount];
IPoint pEdgePoint = new PointClass();
int intEdgeEID;
IPoint pFoundEdgePoint;
double dblEdgePercent;
//用于获取几何网络元素的UserID, UserClassID,UserSubID
INetElements pNetElements = pNetWork as INetElements;
int intEdgeUserClassID;
int intEdgeUserID;
int intEdgeUserSubID;
for (int i = 0; i < intCount; i++)
{
pNetFlag = new EdgeFlagClass();
//获取用户点击点
pEdgePoint = mPointCollectionNet.get_Point(i);
//获取距离用户点击点最近的边
mPointToEID.GetNearestEdge(pEdgePoint, out intEdgeEID, out pFoundEdgePoint, out dblEdgePercent);
if (intEdgeEID <= 0)
{
continue;
}
//根据得到的边查询对应的几何网络中的元素UserID, UserClassID,UserSubID
pNetElements.QueryIDs(intEdgeEID, esriElementType.esriETEdge,
out intEdgeUserClassID, out intEdgeUserID, out intEdgeUserSubID);
if (intEdgeUserClassID <= 0 || intEdgeUserID <= 0)
{
continue;
}
pNetFlag.UserClassID = intEdgeUserClassID;
pNetFlag.UserID = intEdgeUserID;
pNetFlag.UserSubID = intEdgeUserSubID;
pEdgeFlags[i] = pNetFlag as IEdgeFlag;
}
//设置路径求解的边
pTraceFlowSolverGEN.PutEdgeOrigins(ref pEdgeFlags);
//路径计算权重
INetSchema pNetSchema = pNetWork as INetSchema;
INetWeight pNetWeight = pNetSchema.get_WeightByName(weightName);
if (pNetWeight == null)
{
return;
}
//设置权重，这里双向的权重设为一致
INetSolverWeights pNetSolverWeights = pTraceFlowSolverGEN as INetSolverWeights;
pNetSolverWeights.ToFromEdgeWeight = pNetWeight;
pNetSolverWeights.FromToEdgeWeight = pNetWeight;
object[] arrResults = new object[intCount - 1];
//执行路径计算
pTraceFlowSolverGEN.FindPath(esriFlowMethod.esriFMConnected, esriShortestPathObjFn.esriSPObjFnMinSum,
out mEnumNetEID_Junctions, out mEnumNetEID_Edges, intCount - 1, ref arrResults);
//获取路径计算总代价（cost）
mdblPathCost = 0;
for (int i = 0; i < intCount - 1; i++)
{
mdblPathCost += (double)arrResults[i];
}
}
private static ITraceFlowSolver CreateTraceFlowSolver()
{
int num;
INetFlag netFlag;
ITraceFlowSolver solver = (ITraceFlowSolver) new TraceFlowSolverClass();
INetSolver solver2 = (INetSolver)solver;
INetwork network = m_pAnalystGN.Network;
solver2.SourceNetwork = network;
INetElements elements = (INetElements)network;
if (m_JunFlags.Count > 0)
{
IJunctionFlag[] junctionOrigins = new IJunctionFlag[m_JunFlags.Count];
for (num = 0; num < m_JunFlags.Count; num++)
{
netFlag = (m_JunFlags[num] as NetFlagWrap).NetFlag;
junctionOrigins[num] = (IJunctionFlag)netFlag;
}
(solver as ITraceFlowSolverGEN).PutJunctionOrigins(ref junctionOrigins);
}
if (m_EdgeFlags.Count > 0)
{
IEdgeFlag[] edgeOrigins = new IEdgeFlag[m_EdgeFlags.Count];
for (num = 0; num < m_EdgeFlags.Count; num++)
{
netFlag = (m_EdgeFlags[num] as NetFlagWrap).NetFlag;
edgeOrigins[num] = (IEdgeFlag)netFlag;
}
(solver as ITraceFlowSolverGEN).PutEdgeOrigins(ref edgeOrigins);
}
INetSolverWeightsGEN sgen = (INetSolverWeightsGEN)solver;
if (JunctionWeight != null)
{
sgen.JunctionWeight = JunctionWeight;
}
if (FromToEdgeWeight != null)
{
sgen.JunctionWeight = FromToEdgeWeight;
}
if (ToFromEdgeWeight != null)
{
sgen.JunctionWeight = ToFromEdgeWeight;
}
if (JunctionFilterWeight != null)
{
sgen.JunctionFilterWeight = JunctionFilterWeight;
}
if (FromToEdgeFilterWeight != null)
{
sgen.FromToEdgeFilterWeight = FromToEdgeFilterWeight;
}
if (ToFromEdgeFilterWeight != null)
{
sgen.ToFromEdgeFilterWeight = ToFromEdgeFilterWeight;
}
if (JuncfromValues != null)
{
sgen.SetFilterRanges(esriElementType.esriETJunction, ref JuncfromValues, ref JunctoValues);
sgen.SetFilterType(esriElementType.esriETJunction, esriWeightFilterType.esriWFRange,
ApplyJuncFilterWeight);
}
else
{
sgen.SetFilterType(esriElementType.esriETJunction, esriWeightFilterType.esriWFNone,
ApplyJuncFilterWeight);
}
if (EdgefromValues != null)
{
sgen.SetFilterRanges(esriElementType.esriETEdge, ref EdgefromValues, ref EdgetoValues);
sgen.SetFilterType(esriElementType.esriETEdge, esriWeightFilterType.esriWFRange, ApplyJuncFilterWeight);
}
else
{
sgen.SetFilterType(esriElementType.esriETEdge, esriWeightFilterType.esriWFNone, ApplyEdgeFilterWeight);
}
INetElementBarriers netElementBarriers = GetNetElementBarriers(esriElementType.esriETJunction);
if (netElementBarriers != null)
{
(solver as INetSolver).set_ElementBarriers(esriElementType.esriETJunction, netElementBarriers);
}
netElementBarriers = GetNetElementBarriers(esriElementType.esriETEdge);
if (netElementBarriers != null)
{
(solver as INetSolver).set_ElementBarriers(esriElementType.esriETEdge, netElementBarriers);
}
return(solver);
}
private void SetButton_Click(object obj, EventArgs eventArgs)
{
int count = this.listJunctionFlag.Count;
int count2 = this.listEdgeFlag.Count;
int count3 = this.listJunctionBarrier.Count;
int count4 = this.listEdgeBarrier.Count;
if (count < 1 && count2 < 1)
{
MessageBox.Show("请用户选择要分析的点或线!");
}
else
{
IEdgeFlag[] array = new EdgeFlag[count2];
IJunctionFlag[] array2 = new JunctionFlag[count];
INetwork network = null;
INetworkClass networkClass = null;
if (count > 0)
{
for (int i = 0; i < count; i++)
{
IFeature feature = this.listJunctionFlag[i];
networkClass = (feature.Class as INetworkClass);
network = networkClass.GeometricNetwork.Network;
INetElements netElements = network as INetElements;
INetFlag netFlag = new JunctionFlag() as INetFlag;
ISimpleJunctionFeature simpleJunctionFeature = feature as ISimpleJunctionFeature;
int userClassID;
int userID;
int userSubID;
netElements.QueryIDs(simpleJunctionFeature.EID, (esriElementType)1, out userClassID, out userID,
out userSubID);
netFlag.UserClassID = (userClassID);
netFlag.UserID = (userID);
netFlag.UserSubID = (userSubID);
IJunctionFlag junctionFlag = netFlag as IJunctionFlag;
array2[i] = junctionFlag;
}
}
if (count2 > 0)
{
for (int j = 0; j < count2; j++)
{
IFeature feature2 = this.listEdgeFlag[j];
networkClass = (feature2.Class as INetworkClass);
network = networkClass.GeometricNetwork.Network;
INetElements netElements2 = network as INetElements;
INetFlag netFlag2 = new EdgeFlag() as INetFlag;
ISimpleEdgeFeature simpleEdgeFeature = feature2 as ISimpleEdgeFeature;
int userClassID2;
int userID2;
int userSubID2;
netElements2.QueryIDs(simpleEdgeFeature.EID, (esriElementType)2, out userClassID2, out userID2,
out userSubID2);
netFlag2.UserClassID = (userClassID2);
netFlag2.UserID = (userID2);
netFlag2.UserSubID = (userSubID2);
IEdgeFlag edgeFlag = netFlag2 as IEdgeFlag;
array[j] = edgeFlag;
}
}
ITraceFlowSolverGEN traceFlowSolverGEN = new TraceFlowSolver() as ITraceFlowSolverGEN;
INetSolver netSolver = traceFlowSolverGEN as INetSolver;
INetElementBarriersGEN netElementBarriersGEN = new NetElementBarriers();
netElementBarriersGEN.Network = (network);
netElementBarriersGEN.ElementType = (esriElementType)(1);
int[] array3 = new int[count3];
int num = 0;
if (count3 > 0)
{
for (int k = 0; k < count3; k++)
{
IFeature feature3 = this.listJunctionBarrier[k];
networkClass = (feature3.Class as INetworkClass);
network = networkClass.GeometricNetwork.Network;
INetElements netElements3 = network as INetElements;
new EdgeFlag();
ISimpleJunctionFeature simpleJunctionFeature2 = feature3 as ISimpleJunctionFeature;
int num2;
int num3;
netElements3.QueryIDs(simpleJunctionFeature2.EID, (esriElementType)1, out num, out num2,
out num3);
array3[k] = num2;
}
netElementBarriersGEN.SetBarriers(num, ref array3);
netSolver.set_ElementBarriers((esriElementType)1, (INetElementBarriers)netElementBarriersGEN);
}
INetElementBarriersGEN netElementBarriersGEN2 = new NetElementBarriers();
netElementBarriersGEN2.Network = (network);
netElementBarriersGEN2.ElementType = (esriElementType)(2);
int[] array4 = new int[count4];
if (count4 > 0)
{
for (int l = 0; l < count4; l++)
{
IFeature feature4 = this.listEdgeBarrier[l];
networkClass = (feature4.Class as INetworkClass);
network = networkClass.GeometricNetwork.Network;
INetElements netElements4 = network as INetElements;
new EdgeFlag();
ISimpleEdgeFeature simpleEdgeFeature2 = feature4 as ISimpleEdgeFeature;
int num4;
int num5;
netElements4.QueryIDs(simpleEdgeFeature2.EID, (esriElementType)2, out num, out num4, out num5);
array4[l] = num4;
}
netElementBarriersGEN2.SetBarriers(num, ref array4);
netSolver.set_ElementBarriers((esriElementType)2, (INetElementBarriers)netElementBarriersGEN2);
}
netSolver.SourceNetwork = (network);
if (count > 0)
{
traceFlowSolverGEN.PutJunctionOrigins(ref array2);
}
if (count2 > 0)
{
traceFlowSolverGEN.PutEdgeOrigins(ref array);
}
IEnumNetEID enumNetEID = null;
IEnumNetEID enumNetEID2 = null;
object[] array5 = new object[1];
if (this.WayCom.SelectedIndex == 0)
{
traceFlowSolverGEN.FindSource(0, (esriShortestPathObjFn)1, out enumNetEID, out enumNetEID2,
count + count2, ref array5);
}
if (this.WayCom.SelectedIndex == 1)
{
traceFlowSolverGEN.FindSource((esriFlowMethod)1, (esriShortestPathObjFn)1, out enumNetEID,
out enumNetEID2, count + count2, ref array5);
}
IPolyline polyline = new Polyline() as IPolyline;
IGeometryCollection geometryCollection = polyline as IGeometryCollection;
ISpatialReference spatialReference = this.m_iApp.FocusMap.SpatialReference;
IEIDHelper iEIDHelper = new EIDHelper();
iEIDHelper.GeometricNetwork = (networkClass.GeometricNetwork);
iEIDHelper.OutputSpatialReference = (spatialReference);
iEIDHelper.ReturnGeometries = (true);
iEIDHelper.ReturnFeatures = (true);
int selectedIndex = this.LayerCom.SelectedIndex;
if (selectedIndex >= 0 && this.MapControl != null)
{
this.LayerCom.SelectedItem.ToString();
IFeatureLayer ifeatureLayer_ = ((TrackingAnalyForm.LayerInfo) this.LayerCom.SelectedItem)._layer;
if (ifeatureLayer_ != null)
{
IFeatureSelection featureSelection = (IFeatureSelection)ifeatureLayer_;
featureSelection.Clear();
if (enumNetEID2 != null)
{
IEnumEIDInfo enumEIDInfo = iEIDHelper.CreateEnumEIDInfo(enumNetEID2);
int count5 = enumEIDInfo.Count;
enumEIDInfo.Reset();
for (int m = 0; m < count5; m++)
{
IEIDInfo iEIDInfo = enumEIDInfo.Next();
featureSelection.Add(iEIDInfo.Feature);
IGeometry geometry = iEIDInfo.Geometry;
geometryCollection.AddGeometryCollection(geometry as IGeometryCollection);
}
}
featureSelection.SelectionSet.Refresh();
IActiveView activeView = this.m_iApp.ActiveView;
activeView.Refresh();
CMapOperator.ShowFeatureWithWink(this.m_iApp.ActiveView.ScreenDisplay, polyline);
}
}
}
}
public int [] FindPath(IEdgeFlag [] edgeFlags, out double length)
{
int [] roadIDs = {};
int [] roadIDList = {};
length = -1;
try
{
ITraceFlowSolverGEN traceSolver = new TraceFlowSolverClass();
INetSolver solver = traceSolver as INetSolver;
solver.SourceNetwork = roadNet;
traceSolver.PutEdgeOrigins(ref edgeFlags);
INetSchema schema = roadNet as INetSchema;
INetSolverWeights netWeights = traceSolver as INetSolverWeights;
INetWeight weight = schema.get_WeightByName("RoadTypeWeight");
netWeights.FromToEdgeWeight = weight;
netWeights.ToFromEdgeWeight = weight;
IEnumNetEID roadJunctions = new EnumNetEIDArrayClass();
IEnumNetEID roadEdges = new EnumNetEIDArrayClass();
object [] segmentCosts = new object [edgeFlags.Length];
traceSolver.FindPath(esriFlowMethod.esriFMConnected,
esriShortestPathObjFn.esriSPObjFnMinSum,
out roadJunctions, out roadEdges, edgeFlags.Length-1, ref segmentCosts);
// arrays for storing the featureclass id, object id, and subclass id
roadIDs = new int[roadEdges.Count];
int [] FCIDList = new int[roadEdges.Count];
int [] SCIDList = new int[roadEdges.Count];
//index of feature fields of interest
int lengthField = RoadEdge.FindField("SHAPE_Length");
//int roadIDField = RoadEdge.FindField("roadid");
//collect feature ids, roadids, and length by iterating through uid enumeration
length = 0;
roadIDList = new int[roadEdges.Count];
for(int i=0; i < roadEdges.Count; i++)
{
RoadNetElements.QueryIDs(roadEdges.Next(),ESRI.ArcGIS.Geodatabase.esriElementType.esriETEdge,out FCIDList[i], out roadIDs[i],out SCIDList[i]);
roadIDList[i] = RoadEdge.GetFeature(roadIDs[i]).OID;//get_Value(roadIDField).ToString());
length += double.Parse(RoadEdge.GetFeature(roadIDs[i]).get_Value(lengthField).ToString());
}
}
catch(Exception err)
{
throw err;
}
return roadIDList;//roadIDs;
}
public int [] GetRouteFromEdgeFlags(ArrayList edgeFlags)
{
IEdgeFlag [] edgeflags = new IEdgeFlag[edgeFlags.Count];
for(int i = 0; i < edgeFlags.Count; i++)
{
edgeflags[i] = edgeFlags[i] as IEdgeFlag;
}
double length = -1;
return FindPath(edgeflags, out length);
}
private bool slovePath()
{
try
{
double costLength = 0;
m_ipPointToEID = new PointToEIDClass();
m_ipPointToEID.SourceMap = m_HookHelper.ActiveView.FocusMap;
m_ipPointToEID.GeometricNetwork = m_ipGeometricNetwork;
m_ipPointToEID.SnapTolerance = m_HookHelper.ActiveView.Extent.Envelope.Width / 100;
INetwork ipNetwork;
INetElements ipNetElements;
INetFlag ipNetFlag;
//IEdgeFlag ipEdgeFlag;
IEdgeFlag [] EdgeFlagArray;
int pEdgeID;
IPoint pLocation;
double dblEdgePercent;
int iUserClassID, iUserID, iUserSubID;
// look up the EID for the current point (this will populate intEdgeID and dblEdgePercent)
ipTraceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolver;
//get the inetsolver interface
ipNetSolver = ipTraceFlowSolver as INetSolver;
ipNetwork = m_ipGeometricNetwork.Network;
//set the net solver's source network
ipNetSolver.SourceNetwork = ipNetwork;
//make edge flags from the points
//the INetElements interface is needed to get UserID, UserClassID,
//and UserSubID from an element id
EdgeFlagArray = new IEdgeFlag[mPointArray.PointCount];
ipNetElements = ipNetwork as INetElements;
for (int ii = 0; ii < mPointArray.PointCount; ii++)
{
ipNetFlag = new EdgeFlagClass();
point = mPointArray.get_Point(ii);
m_ipPointToEID.GetNearestEdge(point, out pEdgeID, out pLocation, out dblEdgePercent);
//iUserClassID是featureclass ID,userID是objectID,userSubID
ipNetElements.QueryIDs(pEdgeID, esriElementType.esriETEdge, out iUserClassID, out iUserID, out iUserSubID);
ipNetFlag.UserClassID = iUserClassID;
ipNetFlag.UserID = iUserID;
ipNetFlag.UserSubID = iUserSubID;
EdgeFlagArray[ii] = ipNetFlag as IEdgeFlag;
//MessageBox.Show(iUserID.ToString());
}
//添加到TraceFlowSolver中
ITraceFlowSolverGEN pTFSolverGEN = ipTraceFlowSolver as ITraceFlowSolverGEN;
pTFSolverGEN.PutEdgeOrigins(ref EdgeFlagArray);
//get the INetSchema interface
INetSchema ipNetSchema = ipNetwork as INetSchema;
INetWeight ipNetWeight = ipNetSchema.get_Weight(0);
//MessageBox.Show(ipNetWeight.WeightName);
INetSolverWeights ipNetSolverWeights = ipTraceFlowSolver as INetSolverWeights;
//来回使用同一个权重
ipNetSolverWeights.FromToEdgeWeight = ipNetWeight;
ipNetSolverWeights.ToFromEdgeWeight = ipNetWeight;
object [] totalCost = new object[mPointArray.PointCount - 1];
pTFSolverGEN.FindPath(esriFlowMethod.esriFMConnected,
esriShortestPathObjFn.esriSPObjFnMinSum,
out m_ipEnumNetEID_Junctions,
out m_ipEnumNetEID_Edges,
mPointArray.PointCount - 1, ref totalCost);
//MessageBox.Show(totalCost[0].ToString());
for (int ii = 0; ii < totalCost.Length; ii++)
{
costLength += (double)(totalCost[ii]);
}
pStatusBar.Panels[0].Text = "水管长度:" + costLength.ToString();
return(true);
}
catch (Exception e) {
MessageBox.Show(e.Message);
return(false);
}
}
private bool slovePath()
{
try
{
double costLength=0;
m_ipPointToEID = new PointToEIDClass();
m_ipPointToEID.SourceMap=m_HookHelper.ActiveView.FocusMap;
m_ipPointToEID.GeometricNetwork = m_ipGeometricNetwork;
m_ipPointToEID.SnapTolerance=m_HookHelper.ActiveView.Extent.Envelope.Width/100;
INetwork ipNetwork;
INetElements ipNetElements;
INetFlag ipNetFlag;
//IEdgeFlag ipEdgeFlag;
IEdgeFlag [] EdgeFlagArray;
int pEdgeID;
IPoint pLocation;
double dblEdgePercent;
int iUserClassID,iUserID,iUserSubID;
// look up the EID for the current point (this will populate intEdgeID and dblEdgePercent)
ipTraceFlowSolver=new TraceFlowSolverClass() as ITraceFlowSolver;
//get the inetsolver interface
ipNetSolver=ipTraceFlowSolver as INetSolver;
ipNetwork=m_ipGeometricNetwork.Network;
//set the net solver's source network
ipNetSolver.SourceNetwork=ipNetwork;
//make edge flags from the points
//the INetElements interface is needed to get UserID, UserClassID,
//and UserSubID from an element id
EdgeFlagArray=new IEdgeFlag[mPointArray.PointCount];
ipNetElements = ipNetwork as INetElements;
for(int ii=0;ii<mPointArray.PointCount;ii++)
{
ipNetFlag = new EdgeFlagClass();
point=mPointArray.get_Point(ii);
m_ipPointToEID.GetNearestEdge(point,out pEdgeID,out pLocation,out dblEdgePercent);
//iUserClassID��featureclass ID,userID��objectID,userSubID
ipNetElements.QueryIDs(pEdgeID,esriElementType.esriETEdge,out iUserClassID,out iUserID,out iUserSubID);
ipNetFlag.UserClassID=iUserClassID;
ipNetFlag.UserID=iUserID;
ipNetFlag.UserSubID=iUserSubID;
EdgeFlagArray[ii]=ipNetFlag as IEdgeFlag;
//MessageBox.Show(iUserID.ToString());
}
//��ӵ�TraceFlowSolver��
ITraceFlowSolverGEN pTFSolverGEN=ipTraceFlowSolver as ITraceFlowSolverGEN;
pTFSolverGEN.PutEdgeOrigins(ref EdgeFlagArray);
//get the INetSchema interface
INetSchema ipNetSchema = ipNetwork as INetSchema;
INetWeight ipNetWeight = ipNetSchema.get_Weight(0);
//MessageBox.Show(ipNetWeight.WeightName);
INetSolverWeights ipNetSolverWeights=ipTraceFlowSolver as INetSolverWeights;
//����ʹ��ͬһ��Ȩ��
ipNetSolverWeights.FromToEdgeWeight=ipNetWeight;
ipNetSolverWeights.ToFromEdgeWeight=ipNetWeight;
object [] totalCost=new object[mPointArray.PointCount-1];
pTFSolverGEN.FindPath(esriFlowMethod.esriFMConnected,
esriShortestPathObjFn.esriSPObjFnMinSum,
out m_ipEnumNetEID_Junctions,
out m_ipEnumNetEID_Edges,
mPointArray.PointCount-1,ref totalCost);
//MessageBox.Show(totalCost[0].ToString());
for(int ii=0;ii<totalCost.Length;ii++){
costLength+=(double)(totalCost[ii]);
}
pStatusBar.Panels[0].Text="ˮ�ܳ���:"+costLength.ToString();
return true;
}
catch(Exception e){
MessageBox.Show(e.Message);
return false;
}
}
