/// <summary>
/// 设置障碍点
/// </summary>
/// <param name="pNetSolver"></param>
/// <param name="breakPoints"></param>
/// <param name="pGeometricNet"></param>
/// <param name="pDisc"></param>
private void SetBarries(INetSolver pNetSolver, List <IPoint> breakPoints, IGeometricNetwork pGeometricNet, double pDisc)
{
INetwork pNetwork = pGeometricNet.Network;
INetElements pNetElements = pNetwork as INetElements;
IPointToEID pPointToEID = new PointToEIDClass();
pPointToEID.GeometricNetwork = pGeometricNet;
pPointToEID.SnapTolerance = pDisc;
ISelectionSetBarriers barriers = new SelectionSetBarriersClass();
foreach (var point in breakPoints)
{
//寻找障碍点最近的线要素，从而添加爆管线
int nearEdgeID;
IPoint outIpoint;
double precent;
pPointToEID.GetNearestEdge(point, out nearEdgeID, out outIpoint, out precent);
int userClassID;
int userID;
int userSubID;
//查询相关ID
pNetElements.QueryIDs(nearEdgeID, esriElementType.esriETJunction, out userClassID, out userID, out userSubID);
//添加障碍点
barriers.Add(userClassID, userID);
}
pNetSolver.SelectionSetBarriers = barriers;
}
//解决路径
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 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;
}
/// <summary>
/// 最短路径分析
/// </summary>
/// <param name="pMap"></param>
/// <param name="pGeometricNet"></param>
/// <param name="pCollection"></param>
/// <param name="weightName"></param>
/// <param name="pDisc"></param>
/// <returns></returns>
public static IPolyline DistanceFun(IMap pMap, IGeometricNetwork pGeometricNet, IPointCollection pCollection, string weightName, double pDisc)
{
//几何网络分析接口
ITraceFlowSolverGEN pTraceFlowGen = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
INetSolver pNetSolver = pTraceFlowGen as INetSolver;
//获取网络
INetwork pNetwork = pGeometricNet.Network;
//设置该网络为几何分析处理网络
pNetSolver.SourceNetwork = pNetwork;
//网络元素
INetElements pNetElements = pNetwork as INetElements;
//根据输入点建立旗帜数组
IJunctionFlag[] pJunctionFlags = GetJunctionFlags(pMap, pGeometricNet, pCollection, pDisc);
//将旗帜数组添加到处理类中
pTraceFlowGen.PutJunctionOrigins(ref pJunctionFlags);
//设置边线的权重 可以考虑使用点权重
SetWeight(weightName, pTraceFlowGen, pNetwork);
//获取边线和交汇点的集合
IEnumNetEID junctionEIDs;
IEnumNetEID netEIDS;
object[] pRec = new object[pCollection.PointCount - 1];
pTraceFlowGen.FindPath(esriFlowMethod.esriFMConnected, esriShortestPathObjFn.esriSPObjFnMinSum, out junctionEIDs, out netEIDS, pCollection.PointCount - 1, ref pRec);
//获取最短路径
IPolyline pPolyline = new PolylineClass();
IGeometryCollection pGeometryCollection = pPolyline as IGeometryCollection;
ISpatialReference pSpatialReference = pMap.SpatialReference;
IEIDHelper pEIDHelper = new EIDHelperClass();
pEIDHelper.GeometricNetwork = pGeometricNet;
pEIDHelper.OutputSpatialReference = pSpatialReference;
pEIDHelper.ReturnGeometries = true;
IEnumEIDInfo pEnumEIDInfo = pEIDHelper.CreateEnumEIDInfo(netEIDS);
int Count = pEnumEIDInfo.Count;
pEnumEIDInfo.Reset();
for (int i = 0; i < Count; i++)
{
IEIDInfo pEIDInfo = pEnumEIDInfo.Next();
IGeometry pGeometry = pEIDInfo.Geometry;
pGeometryCollection.AddGeometryCollection(pGeometry as IGeometryCollection);
}
return(pPolyline);
}
/// <summary>
/// set barrier for traceFlowSolver
/// </summary>
/// <param name="traceFlowSolver">object ITraceFlowSolverGEN</param>
protected void SetBarriers(ref ITraceFlowSolverGEN traceFlowSolver)
{
INetSolver netSolver = traceFlowSolver as INetSolver;
if (this.edgeBarriers.Length > 0)
{
INetElementBarriersGEN netElementBarriersGEN = new NetElementBarriersClass() as INetElementBarriersGEN;
netElementBarriersGEN.ElementType = esriElementType.esriETEdge;
netElementBarriersGEN.Network = this.geometricNetwork.Network;
netElementBarriersGEN.SetBarriersByEID(ref this.edgeBarriers);
INetElementBarriers netElementBarriers = netElementBarriersGEN as INetElementBarriers;
netSolver.set_ElementBarriers(esriElementType.esriETEdge, netElementBarriers);
}
if (this.junctionBarriers.Length > 0)
{
INetElementBarriersGEN netElementBarriersGEN = new NetElementBarriersClass() as INetElementBarriersGEN;
netElementBarriersGEN.ElementType = esriElementType.esriETJunction;
netElementBarriersGEN.Network = this.geometricNetwork.Network;
netElementBarriersGEN.SetBarriersByEID(ref this.junctionBarriers);
INetElementBarriers netElementBarriers = netElementBarriersGEN as INetElementBarriers;
netSolver.set_ElementBarriers(esriElementType.esriETJunction, netElementBarriers);
}
}
private void StartAnalysis(IFeature feature)
{
if (feature.FeatureType != esriFeatureType.esriFTSimpleJunction)
{
MessageService.Current.Warn("请选择管线点");
return;
}
if (!_pipelineConfig.IsPipelineLayer(feature.Class.AliasName, enumPipelineDataType.Point))
{
MessageService.Current.Warn("请选择管线点");
return;
}
double snapDist = CommonUtils.ConvertPixelsToMapUnits(_context.ActiveView,
_context.Config.SnapTolerance);
IBasicLayerInfo lineConfig =
_plugin.PipeConfig.GetBasicLayerInfo(feature.Class as IFeatureClass) as IBasicLayerInfo;
if (this._startPoint == null && _startEid == 0)
{
//开始记录起始点
IPipelineLayer oldLayer = _pipelineConfig.GetPipelineLayer(feature.Class.AliasName,
enumPipelineDataType.Point);
if (oldLayer == null)
{
MessageService.Current.Warn("你选择的图层不是合法的管线图层!");
return;
}
List <IBasicLayerInfo> basicInfos = oldLayer.GetLayers(enumPipelineDataType.Junction);
IFeatureClass featureClass = basicInfos.Count > 0 ? basicInfos[0].FeatureClass : null;
if (featureClass == null)
{
MessageService.Current.Warn("管线图层没有构建网络图层!");
return;
}
INetworkClass networkClass = featureClass as INetworkClass;
_geometricNetwork = networkClass.GeometricNetwork;
IPointToEID pnToEid = new PointToEIDClass();
pnToEid.GeometricNetwork = _geometricNetwork;
pnToEid.SnapTolerance = snapDist;
pnToEid.SourceMap = _context.FocusMap;
pnToEid.GetNearestJunction(feature.Shape as IPoint, out _startEid, out _startPoint);
return;
}
IPipelineLayer newLayer = _pipelineConfig.GetPipelineLayer(feature.Class.AliasName,
enumPipelineDataType.Point);
if (newLayer == null)
{
MessageService.Current.Warn("你选择的图层不是合法的管线图层!");
return;
}
List <IBasicLayerInfo> basicInfos1 = newLayer.GetLayers(enumPipelineDataType.Junction);
IFeatureClass featureClass2 = basicInfos1.Count > 0 ? basicInfos1[0].FeatureClass : null;
if (featureClass2 == null)
{
MessageService.Current.Warn("第二个管线图层没有构建网络图层!");
return;
}
INetworkClass networkClass2 = featureClass2 as INetworkClass;
if (networkClass2.GeometricNetwork != _geometricNetwork)
{
if (MessageService.Current.Ask("两个点位属于不同的网络图层，使用第二个网络图层作为分析图层吗?") == false)
{
return;
}
_geometricNetwork = networkClass2.GeometricNetwork;
IPointToEID pnToEid = new PointToEIDClass();
pnToEid.GeometricNetwork = _geometricNetwork;
pnToEid.SnapTolerance = snapDist;
pnToEid.SourceMap = _context.FocusMap;
pnToEid.GetNearestJunction(feature.Shape as IPoint, out _startEid, out _startPoint);
return;
}
try
{
IPointToEID pntEid = new PointToEIDClass();
pntEid.GeometricNetwork = _geometricNetwork;
pntEid.SourceMap = _context.FocusMap;
pntEid.SnapTolerance = snapDist;
pntEid.GetNearestJunction(feature.Shape as IPoint, out _endEid, out _endPoint);
if (_endEid < 1)
{
MessageService.Current.Warn("未能找到第二个分析点!");
return;
}
if (_startEid == _endEid)
{
MessageService.Current.Warn("起点终点为同一个点!");
return;
}
INetElements netElements = _geometricNetwork.Network as INetElements;
INetworkClass networkClass = feature.Class as INetworkClass;
IJunctionFlag[] array = new JunctionFlag[2];
INetFlag netFlag = new JunctionFlag() as INetFlag;
int userClassID;
int userID;
int userSubID;
netElements.QueryIDs(_endEid, esriElementType.esriETJunction, out userClassID, out userID, out userSubID);
netFlag.UserClassID = (userClassID);
netFlag.UserID = (userID);
netFlag.UserSubID = (userSubID);
IJunctionFlag value = netFlag as IJunctionFlag;
array.SetValue(value, 0);
INetFlag netFlag2 = new JunctionFlag() as INetFlag;
netElements.QueryIDs(_startEid, esriElementType.esriETJunction, out userClassID,
out userID, out userSubID);
netFlag2.UserClassID = (userClassID);
netFlag2.UserID = (userID);
netFlag2.UserSubID = (userSubID);
value = (netFlag2 as IJunctionFlag);
array.SetValue(value, 1);
ITraceFlowSolverGEN traceFlowSolverGEN = new TraceFlowSolver() as ITraceFlowSolverGEN;
INetSolver netSolver = traceFlowSolverGEN as INetSolver;
netSolver.SourceNetwork = _geometricNetwork.Network;
traceFlowSolverGEN.PutJunctionOrigins(ref array);
object[] array2 = new object[1];
IEnumNetEID enumNetEID;
IEnumNetEID enumNetEID2;
traceFlowSolverGEN.FindPath((esriFlowMethod)2, (esriShortestPathObjFn)1, out enumNetEID,
out enumNetEID2, 1, ref array2);
if (this.ipolyline_0 == null)
{
this.ipolyline_0 = new Polyline() as IPolyline;
}
IGeometryCollection geometryCollection = this.ipolyline_0 as IGeometryCollection;
geometryCollection.RemoveGeometries(0, geometryCollection.GeometryCount);
if (enumNetEID2.Count <= 0)
{
this.ifeature_0 = null;
MessageBox.Show("两点之间不存在路径可以连通！");
}
else
{
ShowShortObjectForm showShortObjectForm = new ShowShortObjectForm(_context);
showShortObjectForm.pApp = _context;
ISpatialReference spatialReference = _context.FocusMap.SpatialReference;
IEIDHelper eIDHelperClass = new EIDHelper();
eIDHelperClass.GeometricNetwork = (networkClass.GeometricNetwork);
eIDHelperClass.OutputSpatialReference = (spatialReference);
eIDHelperClass.ReturnGeometries = (true);
eIDHelperClass.ReturnFeatures = (true);
IEnumEIDInfo enumEIDInfo = eIDHelperClass.CreateEnumEIDInfo(enumNetEID2);
int count = enumEIDInfo.Count;
enumEIDInfo.Reset();
for (int i = 0; i < count; i++)
{
IEIDInfo iEIDInfo = enumEIDInfo.Next();
IGeometry geometry = iEIDInfo.Geometry;
if (i == 0)
{
showShortObjectForm.AddPipeName(this.string_0);
}
showShortObjectForm.AddFeature(iEIDInfo.Feature);
geometryCollection.AddGeometryCollection(geometry as IGeometryCollection);
}
showShortObjectForm.AddShortPath(this.ipolyline_0);
showShortObjectForm.AddLenght(this.ipolyline_0.Length);
this.ifeature_2 = feature;
EsriUtils.ZoomToGeometry(this.ipolyline_0, _context.MapControl.Map, 1.3);
FlashUtility.FlashGeometry(this.ipolyline_0, _context.MapControl);
this.ifeature_0 = null;
_startEid = 0;
_startPoint = null;
_geometricNetwork = null;
showShortObjectForm.Show();
}
}
catch (Exception ex)
{
this.ifeature_0 = null;
MessageBox.Show(ex.Message);
}
}
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);
}
}
}
}
/// <summary>
/// solve isolation
/// </summary>
/// <returns>object JsonObject</returns>
public override JsonObject Solve()
{
try
{
if ((this.edgeFlags.Length == 0) && (this.junctionFlags.Length == 0))
{
string message = "No input valid flags found";
if (this.flagNotFound.Count == 0)
{
throw new GeometricNetworkException(message);
}
else
{
JsonObject error = (new ObjectError(message)).ToJsonObject();
error.AddArray("flagsNotFound", Helper.GetListJsonObjects(this.flagNotFound));
return(error);
}
}
INetwork network = this.geometricNetwork.Network;
ITraceFlowSolverGEN traceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
INetSolver netSolver = traceFlowSolver as INetSolver;
netSolver.SourceNetwork = network;
this.SetFlagsOrigin(ref traceFlowSolver);
////now create barries based on valves
INetElementBarriersGEN netElementBarriersGEN = new NetElementBarriersClass() as INetElementBarriersGEN;
netElementBarriersGEN.ElementType = esriElementType.esriETJunction;
netElementBarriersGEN.Network = network;
int[] userIds = Helper.GetOIDs(this.Valve);
netElementBarriersGEN.SetBarriers(this.Valve.FeatureClassID, ref userIds);
INetElementBarriers netElementBarriers = netElementBarriersGEN as INetElementBarriers;
netSolver.set_ElementBarriers(esriElementType.esriETJunction, netElementBarriers);
IEnumNetEID junctionEIDs;
IEnumNetEID edgeEIDs;
traceFlowSolver.FindFlowEndElements(esriFlowMethod.esriFMConnected, esriFlowElements.esriFEJunctionsAndEdges, out junctionEIDs, out edgeEIDs);
IEIDHelper eidHelper = new EIDHelperClass();
eidHelper.GeometricNetwork = this.geometricNetwork;
eidHelper.ReturnFeatures = true;
List <IEIDInfo> valveEIDInfoHT = Helper.GetEIDInfoListByFeatureClass(this.Valve.FeatureClassID, junctionEIDs, eidHelper);
traceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
netSolver = traceFlowSolver as INetSolver;
netSolver.SourceNetwork = network;
this.SetFlagsOrigin(ref traceFlowSolver);
netSolver.DisableElementClass(this.Station.FeatureClassID);
traceFlowSolver.FindFlowEndElements(esriFlowMethod.esriFMConnected, esriFlowElements.esriFEJunctions, out junctionEIDs, out edgeEIDs);
List <IEIDInfo> sourceEIDInfoHT = Helper.GetEIDInfoListByFeatureClass(this.Station.FeatureClassID, junctionEIDs, eidHelper);
////set up trace to find directly reachable sources
traceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
netSolver = traceFlowSolver as INetSolver;
netSolver.SourceNetwork = network;
this.SetFlagsOrigin(ref traceFlowSolver);
////set barriers in the network based on the saved valves
ISelectionSetBarriers netElementBarrier = new SelectionSetBarriersClass();
foreach (IEIDInfo eidInfo in valveEIDInfoHT)
{
netElementBarrier.Add(this.Valve.FeatureClassID, eidInfo.Feature.OID);
}
netSolver.SelectionSetBarriers = netElementBarrier;
////run trace to find directly reachable sources(without passing valve)
traceFlowSolver.FindFlowElements(esriFlowMethod.esriFMConnected, esriFlowElements.esriFEJunctionsAndEdges, out junctionEIDs, out edgeEIDs);
List <IEIDInfo> sourceDirectEIDInfoHT = Helper.GetEIDInfoListByFeatureClass(this.Station.FeatureClassID, junctionEIDs, eidHelper);
////remove directly reachable sources from source array
foreach (IEIDInfo eidInfo in sourceDirectEIDInfoHT)
{
sourceEIDInfoHT.RemoveAll(x => x.Feature.OID == eidInfo.Feature.OID);
}
List <IEIDInfo> noSourceValveHT = new List <IEIDInfo>();
List <IEIDInfo> hasSourceValveHT = new List <IEIDInfo>();
List <int> listBarrierIds;
bool foundSource;
foreach (IEIDInfo valve in valveEIDInfoHT)
{
foundSource = false;
////create array of all isolation valve excluding the current one
listBarrierIds = new List <int>();
if (valveEIDInfoHT.Count > 1)
{
listBarrierIds.AddRange(valveEIDInfoHT.ConvertAll <int>(i => i.Feature.OID));
if (listBarrierIds.Contains(valve.Feature.OID))
{
listBarrierIds.Remove(valve.Feature.OID);
}
}
////for each source attempt to trace
INetFlag netFlag1;
INetFlag netFlag2;
foreach (IEIDInfo source in sourceEIDInfoHT)
{
////setup trace to test each valve
traceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
netSolver = traceFlowSolver as INetSolver;
netSolver.SourceNetwork = network;
////set the first junction flag for path finding based this current valve
netFlag1 = new JunctionFlagClass();
netFlag1.UserClassID = valve.Feature.Class.ObjectClassID;
netFlag1.UserID = valve.Feature.OID;
netFlag1.UserSubID = 0;
netFlag1.Label = "Origin";
////set the second and last trace flag at this source
netFlag2 = new JunctionFlagClass();
netFlag2.UserClassID = source.Feature.Class.ObjectClassID;
netFlag2.UserID = source.Feature.OID;
netFlag2.UserSubID = 0;
netFlag2.Label = "Destination";
Helper.AddTwoJunctionFlagsToTraceSolver(ref traceFlowSolver, netFlag1, netFlag2);
if (listBarrierIds.Count > 0)
{
netElementBarriersGEN = new NetElementBarriersClass() as INetElementBarriersGEN;
netElementBarriersGEN.ElementType = esriElementType.esriETJunction;
netElementBarriersGEN.Network = network;
int[] barrierIds = listBarrierIds.ToArray();
netElementBarriersGEN.SetBarriers(this.Valve.FeatureClassID, ref barrierIds);
netElementBarriers = netElementBarriersGEN as INetElementBarriers;
netSolver.set_ElementBarriers(esriElementType.esriETJunction, netElementBarriers);
}
////run trace
object[] segCosts = new object[1];
segCosts[0] = new object();
edgeEIDs = null;
traceFlowSolver.FindPath(esriFlowMethod.esriFMConnected, esriShortestPathObjFn.esriSPObjFnMinSum, out junctionEIDs, out edgeEIDs, 1, ref segCosts);
if (edgeEIDs != null && edgeEIDs.Count > 0)
{
foundSource = true;
break;
}
}
if (foundSource)
{
hasSourceValveHT.Add(valve);
}
else
{
noSourceValveHT.Add(valve);
}
}
traceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
netSolver = traceFlowSolver as INetSolver;
netSolver.SourceNetwork = network;
this.SetFlagsOrigin(ref traceFlowSolver);
////set the barrier in the network based on the saved valves
netElementBarrier = new SelectionSetBarriersClass();
foreach (IEIDInfo eidInfo in hasSourceValveHT)
{
netElementBarrier.Add(this.Valve.FeatureClassID, eidInfo.Feature.OID);
}
netSolver.SelectionSetBarriers = netElementBarrier;
////run last trace
traceFlowSolver.FindFlowElements(esriFlowMethod.esriFMConnected, this.FlowElements, out junctionEIDs, out edgeEIDs);
////deal with out put juncEIDs and edgeEIDs
JsonObject objectJson = new JsonObject();
this.SetResults(ref objectJson, edgeEIDs, junctionEIDs);
objectJson.AddArray("flagsNotFound", Helper.GetListJsonObjects(this.flagNotFound));
////valves
if (hasSourceValveHT.Count > 0)
{
JsonObject[] featureSetValves;
this.GetTrace(Helper.GetEnumNetEID(hasSourceValveHT, esriElementType.esriETJunction), out featureSetValves);
objectJson.AddJsonObject("valves", featureSetValves[0]);
}
return(objectJson);
}
catch (Exception e)
{
return((new ObjectError(e.Message)).ToJsonObject());
}
}
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 void setUpstreamJunctionFlag(IPoint pPoint)
{
try
{
if (pGeometricNetwork == null)
{
openGeoNetwork();
}
IEnumNetEID m_ipEnumNetEID_Junctions, m_ipEnumNetEID_Edges; //用来存放搜索结果的Junctions和Edges
ipNetwork = pGeometricNetwork.Network;
/*
* IMarkerSymbol pMarkerSym=SpatialHelperFunction.FindMarkerSym("Civic.ServerStyle", "Default", "Pin Flag Square");
* if(pMarkerSym==null){
* MessageBox.Show("无法取得Symbol");
* return;
* }
*
* IRgbColor pRGBColor=Utility.getRGBColor(255, 0, 0);
*
* pMarkerSym.Color=pRGBColor;
* pMarkerSym.Size=26;
* pMarkerSym.XOffset=10;
* pMarkerSym.YOffset=10;
* //无法使用 IFlagDisplay，只有 INetworkAnalysisExtFlags才可以使用IFlagDisplay
*
* Utility.drawPointToGraphicLayer(m_HookHelper.ActiveView,pPoint,pMarkerSym);
*/
int pJunctionID;
IPoint pLocation;
int iUserClassID, iUserID, iUserSubID;
INetElements ipNetElements;
ipNetElements = ipNetwork as INetElements; //得到network中所有的Elements
mPointToEID = new PointToEIDClass();
//设置它的GeometricNetwork
mPointToEID.GeometricNetwork = pGeometricNetwork;
//设置容限值
mPointToEID.SnapTolerance = m_HookHelper.ActiveView.Extent.Envelope.Width / 100;
//设置Source Map
mPointToEID.SourceMap = m_HookHelper.ActiveView.FocusMap;
//得到最近的JunctionID
mPointToEID.GetNearestJunction(pPoint, out pJunctionID, out pLocation);
//Each element in a logical network has a unique element ID (EID).
ipNetElements.QueryIDs(pJunctionID, esriElementType.esriETJunction, out iUserClassID, out iUserID, out iUserSubID);
//得到所在Junction的ID
//无法使用 IFlagDisplay，只有 INetworkAnalysisExtFlags才可以使用IFlagDisplay
INetFlag ipNetFlag = new JunctionFlagClass();
ipNetFlag.UserClassID = iUserClassID;
ipNetFlag.UserID = iUserID;
ipNetFlag.UserSubID = iUserSubID;
IJunctionFlag [] pArrayJFlag = new IJunctionFlag[1];
pArrayJFlag[0] = ipNetFlag as IJunctionFlag;
ipTraceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
//get the inetsolver interface
ipNetSolver = ipTraceFlowSolver as INetSolver;
ipNetSolver.SourceNetwork = pGeometricNetwork.Network;
ipTraceFlowSolver.PutJunctionOrigins(ref pArrayJFlag);
ipTraceFlowSolver.TraceIndeterminateFlow = false;
object[] totalCost = new object[1];
//ipTraceFlowSolver.FindSource(esriFlowMethod.esriFMUpstream,esriShortestPathObjFn.esriSPObjFnMinSum,out m_ipEnumNetEID_Junctions,out m_ipEnumNetEID_Edges,1,ref totalCost);
ipTraceFlowSolver.FindFlowElements(esriFlowMethod.esriFMUpstream,
esriFlowElements.esriFEEdges,
out m_ipEnumNetEID_Junctions,
out m_ipEnumNetEID_Edges);
SpatialHelperFunction.pathToPolyline(pGeometricNetwork, m_HookHelper.ActiveView, m_ipEnumNetEID_Edges);
//ipTraceFlowSolver.FindFlowEndElements(esriFlowMethod.esriFMUpstream,esriFlowElements.esriFEJunctionsAndEdges,out m_ipEnumNetEID_Junctions,out m_ipEnumNetEID_Edges);
m_HookHelper.ActiveView.Refresh();
}
catch (Exception e) {
MessageBox.Show(e.Message);
}
}
private void setUpstreamJunctionFlag(IPoint pPoint)
{
try
{
if(pGeometricNetwork==null){
openGeoNetwork();
}
IEnumNetEID m_ipEnumNetEID_Junctions,m_ipEnumNetEID_Edges;//����������������Junctions��Edges
ipNetwork=pGeometricNetwork.Network;
/*
IMarkerSymbol pMarkerSym=SpatialHelperFunction.FindMarkerSym("Civic.ServerStyle", "Default", "Pin Flag Square");
if(pMarkerSym==null){
MessageBox.Show("��ȡ��Symbol");
return;
}
IRgbColor pRGBColor=Utility.getRGBColor(255, 0, 0);
pMarkerSym.Color=pRGBColor;
pMarkerSym.Size=26;
pMarkerSym.XOffset=10;
pMarkerSym.YOffset=10;
//��ʹ�� IFlagDisplay��ֻ�� INetworkAnalysisExtFlags�ſ���ʹ��IFlagDisplay
Utility.drawPointToGraphicLayer(m_HookHelper.ActiveView,pPoint,pMarkerSym);
*/
int pJunctionID;
IPoint pLocation;
int iUserClassID,iUserID,iUserSubID;
INetElements ipNetElements;
ipNetElements = ipNetwork as INetElements; //�õ�network�����е�Elements
mPointToEID=new PointToEIDClass();
//��������GeometricNetwork
mPointToEID.GeometricNetwork=pGeometricNetwork;
//��������ֵ
mPointToEID.SnapTolerance=m_HookHelper.ActiveView.Extent.Envelope.Width/100;
//����Source Map
mPointToEID.SourceMap=m_HookHelper.ActiveView.FocusMap;
//�õ������JunctionID
mPointToEID.GetNearestJunction(pPoint,out pJunctionID,out pLocation);
//Each element in a logical network has a unique element ID (EID).
ipNetElements.QueryIDs(pJunctionID,esriElementType.esriETJunction,out iUserClassID,out iUserID,out iUserSubID);
//�õ�����Junction��ID
//��ʹ�� IFlagDisplay��ֻ�� INetworkAnalysisExtFlags�ſ���ʹ��IFlagDisplay
INetFlag ipNetFlag=new JunctionFlagClass();
ipNetFlag.UserClassID=iUserClassID;
ipNetFlag.UserID=iUserID;
ipNetFlag.UserSubID=iUserSubID;
IJunctionFlag [] pArrayJFlag=new IJunctionFlag[1];
pArrayJFlag[0]=ipNetFlag as IJunctionFlag;
ipTraceFlowSolver=new TraceFlowSolverClass() as ITraceFlowSolverGEN;
//get the inetsolver interface
ipNetSolver=ipTraceFlowSolver as INetSolver;
ipNetSolver.SourceNetwork=pGeometricNetwork.Network;
ipTraceFlowSolver.PutJunctionOrigins(ref pArrayJFlag);
ipTraceFlowSolver.TraceIndeterminateFlow=false;
object[] totalCost=new object[1];
//ipTraceFlowSolver.FindSource(esriFlowMethod.esriFMUpstream,esriShortestPathObjFn.esriSPObjFnMinSum,out m_ipEnumNetEID_Junctions,out m_ipEnumNetEID_Edges,1,ref totalCost);
ipTraceFlowSolver.FindFlowElements(esriFlowMethod.esriFMUpstream,
esriFlowElements.esriFEEdges,
out m_ipEnumNetEID_Junctions,
out m_ipEnumNetEID_Edges);
SpatialHelperFunction.pathToPolyline(pGeometricNetwork,m_HookHelper.ActiveView,m_ipEnumNetEID_Edges);
//ipTraceFlowSolver.FindFlowEndElements(esriFlowMethod.esriFMUpstream,esriFlowElements.esriFEJunctionsAndEdges,out m_ipEnumNetEID_Junctions,out m_ipEnumNetEID_Edges);
m_HookHelper.ActiveView.Refresh();
}
catch(Exception e){
MessageBox.Show(e.Message);
}
}
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 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];
}
}
public static List <int> StreamTrace(IGeometricNetwork geometricNetwork, StartFlagEdge edge, List <int> disabledFeatureClassIds, StopperJunctions stoppers, bool isUpStream, int maxFeatureCount, ServerLogger logger)
{
esriFlowMethod direction = isUpStream ? esriFlowMethod.esriFMUpstream : esriFlowMethod.esriFMDownstream;
if (null == geometricNetwork || null == edge || maxFeatureCount <= 0)
{
return(null);
}
ITraceFlowSolverGEN traceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
INetSolver netSolver = traceFlowSolver as INetSolver;
netSolver.SourceNetwork = geometricNetwork.Network;
INetFlag netFlag = new EdgeFlagClass();
netFlag.UserClassID = edge.FeatureClassID;
netFlag.UserID = edge.FeatureID;
//no idea when to assign -1, when to do 0
netFlag.UserSubID = -1;
traceFlowSolver.PutEdgeOrigins(new IEdgeFlag[1] {
netFlag as IEdgeFlag
});
if (null != disabledFeatureClassIds)
{
foreach (int il in disabledFeatureClassIds)
{
if (il > 0)
{
netSolver.DisableElementClass(il);
}
}
}
if (null != stoppers && null != stoppers.Stoppers && stoppers.Stoppers.Length > 0)
{
INetElementBarriersGEN netBarriersGEN = null;
netBarriersGEN = new NetElementBarriersClass() as INetElementBarriersGEN;
netBarriersGEN.ElementType = esriElementType.esriETJunction;
netBarriersGEN.Network = geometricNetwork.Network;
netBarriersGEN.SetBarriers(stoppers.FeatureClassID, stoppers.Stoppers);
netSolver.set_ElementBarriers(esriElementType.esriETJunction, netBarriersGEN as INetElementBarriers);
}
IEnumNetEID junctionEIDs = null;
IEnumNetEID edgeEIDs = null;
traceFlowSolver.TraceIndeterminateFlow = false;
try
{
traceFlowSolver.FindFlowElements(direction, esriFlowElements.esriFEEdges, out junctionEIDs, out edgeEIDs);
if (null != edgeEIDs)
{
if (edgeEIDs.Count <= maxFeatureCount)
{
IEIDHelper eidHelper = new EIDHelperClass();
eidHelper.GeometricNetwork = geometricNetwork;
eidHelper.ReturnGeometries = false;
eidHelper.ReturnFeatures = true;
//eidHelper.AddField("FType");
IEnumEIDInfo eidInfos = eidHelper.CreateEnumEIDInfo(edgeEIDs);
eidInfos.Reset();
IEIDInfo eidInfo = null;
List <int> ftrs = new List <int>();
//IFeature cadFtr = null;
//int ftype;
while ((eidInfo = eidInfos.Next()) != null)
{
ftrs.Add(eidInfo.Feature.OID);
/*cadFtr = eidInfo.Feature;
* if (null != cadFtr.get_Value(edgeTypeId) && int.TryParse(cadFtr.get_Value(edgeTypeId).ToString(), out ftype))
* {
* if(460 == ftype || 558 == ftype)
* ftrs.Add(cadFtr);
* }*/
}
return(ftrs);
}
}
}
catch (Exception e)
{
if (null != logger)
{
logger.LogMessage(ServerLogger.msgType.error, typeof(AOUtilities).Name + "." + System.Reflection.MethodBase.GetCurrentMethod().Name, ErrorCode, e.Message);
}
}
finally
{
ReleaseCOMObj(traceFlowSolver);
}
return(null);
}
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);
}
/// <summary>
/// operation solve trace
/// </summary>
/// <returns>solve trace</returns>
public override JsonObject Solve()
{
try
{
if ((this.edgeFlags.Length == 0) && (this.junctionFlags.Length == 0))
{
string message = "No input valid flags found";
if (this.flagNotFound.Count == 0)
{
throw new GeometricNetworkException(message);
}
else
{
JsonObject error = (new ObjectError(message)).ToJsonObject();
error.AddArray("flagsNotFound", Helper.GetListJsonObjects(this.flagNotFound));
if (this.barrierNotFound.Count > 0)
{
error.AddArray("barriersNotFound", Helper.GetListJsonObjects(this.barrierNotFound));
}
return(error);
}
}
ITraceFlowSolverGEN traceFlowSolver = new TraceFlowSolverClass() as ITraceFlowSolverGEN;
INetSolver netSolver = traceFlowSolver as INetSolver;
INetwork network = this.geometricNetwork.Network;
netSolver.SourceNetwork = network;
// flag origin
this.SetFlagsOrigin(ref traceFlowSolver);
// barrier
this.SetBarriers(ref traceFlowSolver);
// set disabled layers
foreach (int i in this.DisableLayers)
{
netSolver.DisableElementClass(i);
}
// set weight
this.SetWeights(ref traceFlowSolver);
if ((this.TraceSolverType != TraceSolverType.FindCircuits) && (this.TraceSolverType != TraceSolverType.FindCommonAncestors))
{
// The TraceIndeterminateFlow property affects the trace results of trace solvers whose FlowMethod parameter is esriFMUpstream or esriFMDownstream
if ((this.FlowMethod == esriFlowMethod.esriFMDownstream) || (this.FlowMethod == esriFlowMethod.esriFMUpstream))
{
traceFlowSolver.TraceIndeterminateFlow = this.TraceIndeterminateFlow;
}
}
IEnumNetEID junctionEIDs = null;
IEnumNetEID edgeEIDs = null;
object totalCost;
int count;
object[] segmentCosts;
JsonObject objectJson = new JsonObject();
switch (this.TraceSolverType)
{
case TraceSolverType.FindAccumulation:
traceFlowSolver.FindAccumulation(this.FlowMethod, this.FlowElements, out junctionEIDs, out edgeEIDs, out totalCost);
objectJson.AddDouble("totalCost", (double)totalCost);
break;
case TraceSolverType.FindCircuits:
traceFlowSolver.FindCircuits(this.FlowElements, out junctionEIDs, out edgeEIDs);
break;
case TraceSolverType.FindCommonAncestors:
traceFlowSolver.FindCommonAncestors(this.FlowElements, out junctionEIDs, out edgeEIDs);
break;
case TraceSolverType.FindFlowElements:
traceFlowSolver.FindFlowElements(this.FlowMethod, this.FlowElements, out junctionEIDs, out edgeEIDs);
break;
case TraceSolverType.FindFlowEndElements:
traceFlowSolver.FindFlowEndElements(this.FlowMethod, this.FlowElements, out junctionEIDs, out edgeEIDs);
break;
case TraceSolverType.FindFlowUnreachedElements:
traceFlowSolver.FindFlowUnreachedElements(this.FlowMethod, this.FlowElements, out junctionEIDs, out edgeEIDs);
break;
case TraceSolverType.FindPath:
count = Math.Max(this.junctionFlags.Length, this.edgeFlags.Length);
if (count < 2)
{
throw new GeometricNetworkException("Edge or junction found < 2!");
}
--count;
segmentCosts = new object[count];
traceFlowSolver.FindPath(this.FlowMethod, this.ShortestPathObjFn, out junctionEIDs, out edgeEIDs, count, ref segmentCosts);
objectJson.AddArray("segmentCosts", segmentCosts);
break;
case TraceSolverType.FindSource:
count = this.junctionFlags.Length + this.edgeFlags.Length;
segmentCosts = new object[count];
traceFlowSolver.FindSource(this.FlowMethod, this.ShortestPathObjFn, out junctionEIDs, out edgeEIDs, count, ref segmentCosts);
objectJson.AddArray("segmentCosts", segmentCosts);
break;
case TraceSolverType.FindLongest:
// get end junctions in upstream
IEnumNetEID junctionEIDsFindLongest = null;
IEnumNetEID edgeEIDsFindLongest = null;
traceFlowSolver.FindFlowEndElements(esriFlowMethod.esriFMUpstream, esriFlowElements.esriFEJunctions, out junctionEIDsFindLongest, out edgeEIDsFindLongest);
int? eidLongest = null;
double eidLongestLenght = double.MinValue;
if (junctionEIDsFindLongest.Count > 0)
{
double eidLongestLenghtCurrent = double.NaN;
for (int i = 0; i < junctionEIDsFindLongest.Count; i++)
{
int netEIDCurrent = junctionEIDsFindLongest.Next();
object[] segmentLenghts = new object[1];
IEnumNetEID junctionEIDsLongest = null;
IEnumNetEID edgeEIDsLongest = null;
INetFlag netFlag = new JunctionFlagClass();
INetElements netElements = this.geometricNetwork.Network as INetElements;
int featureClassID, featureID, subID;
netElements.QueryIDs(netEIDCurrent, esriElementType.esriETJunction, out featureClassID, out featureID, out subID);
netFlag.UserClassID = featureClassID;
netFlag.UserID = featureID;
netFlag.UserSubID = subID;
IJunctionFlag[] junctionFlags = new IJunctionFlag[] { this.junctionFlags[0], netFlag as IJunctionFlag };
traceFlowSolver.PutJunctionOrigins(ref junctionFlags);
traceFlowSolver.FindPath(esriFlowMethod.esriFMUpstream, esriShortestPathObjFn.esriSPObjFnMinMax, out junctionEIDsLongest, out edgeEIDsLongest, 1, ref segmentLenghts);
if (segmentLenghts[0] != null)
{
eidLongestLenghtCurrent = Convert.ToDouble(segmentLenghts[0]);
if (eidLongestLenghtCurrent > eidLongestLenght)
{
eidLongestLenght = eidLongestLenghtCurrent;
eidLongest = netEIDCurrent;
edgeEIDs = edgeEIDsLongest;
junctionEIDs = junctionEIDsLongest;
}
}
}
}
else
{
throw new GeometricNetworkException("Junction end not found!");
}
if (eidLongest.HasValue)
{
objectJson.AddDouble("totalCost", eidLongestLenght);
}
else
{
throw new GeometricNetworkException("EID longest not found!");
}
break;
default:
throw new GeometricNetworkException("Trace solver type not found");
}
this.SetResults(ref objectJson, edgeEIDs, junctionEIDs);
objectJson.AddArray("flagsNotFound", Helper.GetListJsonObjects(this.flagNotFound));
objectJson.AddArray("barriersNotFound", Helper.GetListJsonObjects(this.barrierNotFound));
return(objectJson);
}
catch (Exception e)
{
return((new ObjectError(e.Message)).ToJsonObject());
}
}
