private IVector GetClosingVector(List <IVector3D> TraverseCourses, IPoint StartPoint, IPoint EndPoint, out double SUMofLengths)
{
IVector SumVec = null;
SUMofLengths = 0;
for (int i = 0; i < TraverseCourses.Count - 1; i++)
{
if (i == 0)
{
SUMofLengths = TraverseCourses[0].Magnitude + TraverseCourses[1].Magnitude;
SumVec = TraverseCourses[0].AddVector(TraverseCourses[1]);
}
else
{
IVector3D SumVec3D = SumVec as IVector3D;
SUMofLengths += TraverseCourses[i + 1].Magnitude;
SumVec = SumVec3D.AddVector(TraverseCourses[i + 1]);
}
}
double dCalcedEndX = StartPoint.X + SumVec.ComponentByIndex[0];
double dCalcedEndY = StartPoint.Y + SumVec.ComponentByIndex[1];
IVector3D CloseVector3D = new Vector3DClass();
CloseVector3D.SetComponents(dCalcedEndX - EndPoint.X, dCalcedEndY - EndPoint.Y, 0);
IVector CloseVector = CloseVector3D as IVector;
return(CloseVector);
}
/// <summary>
/// 更新半径
/// </summary>
/// <param name="radius"></param>
public void UpdatePosition(double radius)
{
this.radius = radius;
IZAware zAware = (IGeometry)centerPoint as IZAware;
zAware.ZAware = true;
upperAxisVector3D = new Vector3DClass();
upperAxisVector3D.SetComponents(0, 0, 10);
lowerAxisVector3D = new Vector3DClass();
lowerAxisVector3D.SetComponents(0, 0, -10);
lowerAxisVector3D.XComponent -= vectorComponentOffset;
normalVector3D = upperAxisVector3D.CrossProduct(lowerAxisVector3D) as IVector3D;
normalVector3D.Magnitude = this.radius;
double rotationAngleInRadians = 2 * (Math.PI / 180);
geometryCollection = new MultiPatchClass();
pointCollection.RemovePoints(0, pointCollection.PointCount);
for (int i = 0; i < 180; i++)
{
normalVector3D.Rotate(-1 * rotationAngleInRadians, upperAxisVector3D);
IPoint vertexPoint = new PointClass();
vertexPoint.X = centerPoint.X + normalVector3D.XComponent;
vertexPoint.Y = centerPoint.Y + normalVector3D.YComponent;
vertexPoint.Z = centerPoint.Z;
pointCollection.AddPoint(vertexPoint, missing, missing);
}
base.Geometry = pointCollection as IGeometry;
this.Update();
}
public IGeometry CreateGeometry(double radius, IPolyline polyline, double qdgc, double zdgc)
{
IPointCollection pointCollection = CreatePointCollectionForCircle(radius);
IVector3D pVectorZ = new Vector3DClass();
pVectorZ.SetComponents(0, 0, 1);
IConstructMultiPatch patch = new MultiPatchClass();
IZAware zAware = pointCollection as IZAware;
if (zAware == null)
{
return(null);
}
zAware.ZAware = true;
// 依据管线长度拉伸
patch.ConstructExtrude(polyline.Length, pointCollection as IGeometry);
// 依据管线角度旋转
IVector3D pVector3D = new Vector3DClass();
pVector3D.SetComponents(polyline.ToPoint.X - polyline.FromPoint.X, polyline.ToPoint.Y - polyline.FromPoint.Y, zdgc - qdgc);
double rotateAngle = Math.Acos(pVector3D.ZComponent / pVector3D.Magnitude);
IVector3D vectorAxis = pVectorZ.CrossProduct(pVector3D) as IVector3D;
ITransform3D transform3D = patch as ITransform3D;
transform3D.RotateVector3D(vectorAxis, rotateAngle);
// 平移到指定位置
transform3D.Move3D(polyline.FromPoint.X, polyline.FromPoint.Y, qdgc);
return(patch as IGeometry);
}
public static IVector3D ConstructVector3D(double xComponent, double yComponent, double zComponent)
{
IVector3D vector3D = new Vector3DClass();
vector3D.SetComponents(xComponent, yComponent, zComponent);
return(vector3D);
}
public static IVector3D ConstructVector3D(double xComponent, double yComponent, double zComponent)
{
IVector3D vector3D = new Vector3DClass();
vector3D.SetComponents(xComponent, yComponent, zComponent);
return vector3D;
}
/// <summary>
/// create a Vector3D from components
/// </summary>
/// <param name="componentX">component X</param>
/// <param name="componentY">component Y</param>
/// <param name="componentZ">component Z</param>
/// <returns>object Vector 3D</returns>
public static IVector3D ConstructVector3D(double componentX, double componentY, double componentZ)
{
IVector3D vector3D = new Vector3DClass();
vector3D.SetComponents(componentX, componentY, componentZ);
return(vector3D);
}
public override IGeometry CreateGeometry()
{
IPointCollection pointCollection = new PolygonClass();
IZAware zAware = pointCollection as IZAware;
zAware.ZAware = true;
IPoint point = new PointClass();
point.X = -_width / 2;
point.Y = -_height / 2;
point.Z = 0;
pointCollection.AddPoint(point);
point = new PointClass();
point.X = -_width / 2;
point.Y = _height / 2;
point.Z = 0;
pointCollection.AddPoint(point);
point = new PointClass();
point.X = _width / 2;
point.Y = _height / 2;
point.Z = 0;
pointCollection.AddPoint(point);
point = new PointClass();
point.X = _width / 2;
point.Y = -_height / 2;
point.Z = 0;
pointCollection.AddPoint(point);
((IPolygon)pointCollection).Close();
IConstructMultiPatch patch = new MultiPatchClass();
patch.ConstructExtrude(_polyline.Length, pointCollection as IGeometry);
IVector3D vectorZ = new Vector3DClass();
vectorZ.SetComponents(0, 0, 1);
IVector3D vector3D = new Vector3DClass();
vector3D.SetComponents(_polyline.ToPoint.X - _polyline.FromPoint.X, _polyline.ToPoint.Y - _polyline.FromPoint.Y, _zdgc - _qdgc);
double rotateAngle = Math.Acos(vector3D.ZComponent / vector3D.Magnitude);
IVector3D vectorAxis = vectorZ.CrossProduct(vector3D) as IVector3D;
ITransform3D transform3D = patch as ITransform3D;
transform3D.RotateVector3D(vectorAxis, rotateAngle);
transform3D.Move3D(_polyline.FromPoint.X, _polyline.FromPoint.Y, _qdgc);
return(patch as IGeometry);
}
public void LoadModel(string sFile)
{
IImport3DFile file = new Import3DFileClass();
file.CreateFromFile(sFile);
IGeometry geometry = file.Geometry;
this.m_pMarker3DSymbol.Shape = geometry;
if (System.IO.Path.GetExtension(sFile).ToLower() == ".wrl")
{
IVector3D pAxis = new Vector3DClass();
pAxis.SetComponents(1.0, 0.0, 0.0);
this.RotateGeometry(this.m_pMarker3DSymbol, pAxis, 90.0);
}
this.m_pMarker3DSymbol.UseMaterialDraping = true;
this.m_pMarker3DSymbol.RestrictAccessToShape();
}
protected override void OnClick()
{
ICadastralEditor pCadEd = (ICadastralEditor)ArcMap.Application.FindExtensionByName("esriCadastralUI.CadastralEditorExtension");
IParcelEditManager pParcEditorMan = (IParcelEditManager)pCadEd;
ICadastralPacketManager pCadPacketMan = (ICadastralPacketManager)pCadEd;
//bool bStartedWithPacketOpen = pCadPacketMan.PacketOpen;
if (pParcEditorMan == null)
{
return;
}
IEditor pEd = (IEditor)ArcMap.Application.FindExtensionByName("esri object editor");
if (pEd.EditState == esriEditState.esriStateNotEditing)
{
MessageBox.Show("Please start editing and try again.");
return;
}
IParcelConstruction pConstr = pParcEditorMan.ParcelConstruction;
IParcelConstruction4 pConstr4 = pConstr as IParcelConstruction4;
ICadastralPoints pCadastralPts = pConstr4 as ICadastralPoints;
ICadastralEditorSettings2 pCadastralEditorSettings2 = pCadEd as ICadastralEditorSettings2;
ICadastralFixedPoints pFixedPoints = pCadastralPts as ICadastralFixedPoints;
IPointCalculation pPointCalc = new PointCalculationClass();
if (pConstr == null)
{
return;
}
IGSLine pParcelLine = null;
IMetricUnitConverter pMetricUnitConv = (IMetricUnitConverter)pCadEd;
IGSPoint pStartPoint = null;
List <int> lstPointIds = new List <int>();
List <IVector3D> Traverse = new List <IVector3D>();
//get rotation here
IParcelConstructionData pConstrData = pConstr4.ConstructionData;
IConstructionParentParcels pConstructionParentParcels = pConstrData as IConstructionParentParcels;
ICadastralUndoRedo pCadUndoRedo = pConstr as ICadastralUndoRedo;
try
{
int iParcelID = -1;
if (pConstructionParentParcels.ParentParcelCount > 0)
{
pConstructionParentParcels.GetParentParcel(0, ref iParcelID);
}
ICadastralParcel pCadaParcel = pCadPacketMan.JobPacket as ICadastralParcel;
IGSParcel pGSParcel = null;
if (pCadaParcel != null)
{
pGSParcel = pCadaParcel.GetParcel(iParcelID);
}
//if in measurement view then rotation is 0
double TheRotation = 0;
if (pGSParcel == null)
{
pGSParcel = pConstr.Parcel;
}
if (!pCadastralEditorSettings2.MeasurementView)
{
TheRotation = pGSParcel.Rotation;//radians
}
if (TheRotation == 123456789)
{
TheRotation = 0;
}
pPointCalc.Rotation = TheRotation;
IGSPoint pClosingPoint = null;
#region simple method as fall-back
bool bUseSimpleStackSelection = false;
if (bUseSimpleStackSelection)
{
//bool bLineSelectionSequence = false;
//IGSLine pLastSelectedGSLineInGrid = null;
//for (int i = 0; i < pConstr.LineCount; i++)
//{
// if (pConstr.GetLineSelection(i))
// {
// if (pConstr.GetLine(i, ref pParcelLine))
// {
// if (!bLineSelectionSequence) //first line
// {
// pStartPoint = pCadastralPts.GetPoint(pParcelLine.FromPoint);
// pToPoint = pCadastralPts.GetPoint(pParcelLine.ToPoint);
// }
// pPointCalc.AddLine(pParcelLine);
// pLastSelectedGSLineInGrid = pParcelLine;
// pToPoint = pCadastralPts.GetPoint(pParcelLine.ToPoint);
// lstPointIds.Add(pToPoint.Id);
// }
// bLineSelectionSequence = true;
// double dBear = pParcelLine.Bearing; //Azimuth of IVector3D is north azimuth radians zero degrees north
// double dDist = pParcelLine.Distance;
// IVector3D vec = new Vector3DClass();
// vec.PolarSet(dBear, 0, dDist); ////Azimuth of IVector3D is north azimuth radians zero degrees north
// Traverse.Add(vec);
// }
// else
// {
// if (bLineSelectionSequence && pConstr.GetLine(i, ref pParcelLine) && HasLineSelectionAfter(pConstr, i))
// //if there was a prior selection and this line is a complete line, and there is no later selection
// {
// MessageBox.Show("Please select a continuous set of lines for closure.");
// return;
// }
// }
//}
//pClosingPoint = pCadastralPts.GetPoint(pLastSelectedGSLineInGrid.ToPoint);
}
else
#endregion
{//build a forward star for the selected lines
IEnumCELines pCELines = new EnumCELinesClass();
IEnumGSLines pEnumGSLines = (IEnumGSLines)pCELines;
ILongArray pLongArray = new LongArrayClass();
int iFirstToNode = -1;
for (int i = 0; i < pConstr.LineCount; i++)
{
if (pConstr.GetLineSelection(i))
{
if (pConstr.GetLine(i, ref pParcelLine))
{
if (iFirstToNode < 0)
{
iFirstToNode = pParcelLine.ToPoint;
}
pLongArray.Add(i);
pCELines.Add(pParcelLine);
}
}
}
if (pCELines.Count == 0)
{
MessageBox.Show("No lines selected. Please select a continuous set of lines for closure." + Environment.NewLine +
"Line selection should not have branches.", "Traverse");
return;
}
IParcelLineFunctions3 ParcelLineFx = new ParcelFunctionsClass();
IGSForwardStar pFwdStar = ParcelLineFx.CreateForwardStar(pEnumGSLines);
//forward star object is now created for all the selected lines,
//need to first re-sequence the lines, and test for branching and discontinuity
int iBranches = 0; int iTracedLines = 0;
int iLoops = 0; int iTerminals = 0;
List <int> LineIDList = new List <int>();
List <int> FromList = new List <int>();
List <int> ToList = new List <int>();
List <string> FromToLine = new List <string>();
bool bTraceSucceeded = TraceLines(ref pFwdStar, iFirstToNode, ref iBranches, ref iTracedLines,
ref iLoops, ref iTerminals, ref FromToLine, ref FromList, ref ToList, 0);
if (iBranches > 0)
{
MessageBox.Show("Please select a continuous set of lines for closure." + Environment.NewLine +
"Line selection should not have branches.", "Traverse");
return;
}
if (iTracedLines < pLongArray.Count)
{
MessageBox.Show("Please select a continuous set of lines for closure." + Environment.NewLine +
"Selected Lines should be connected in a single sequence without branches.", "Traverse");
return;
}
//if it's a single loop check to see if the sequence needs to be reversed
//CW or CCW based on bearings
if (iLoops == 1)
{
bool bIsReversed = false;
foreach (int i in FromList)
{
if (i < 0)
{
bIsReversed = true;
}
else
{
bIsReversed = false;
break;
}
}
if (bIsReversed)
{//all courses are running reversed, so reverse the whole sequence
FromToLine.Clear();
FromList.Reverse();
ToList.Reverse();
int iNewFrom = -ToList[ToList.Count - 1];
int iNewTo = -FromList[ToList.Count - 1];
string sNewFromTo = iNewFrom.ToString() + "," + iNewTo.ToString();
FromToLine.Add(sNewFromTo);
for (int i = 1; i < ToList.Count; i++)
{
iNewFrom = -ToList[i - 1];
iNewTo = -FromList[i - 1];
sNewFromTo = iNewFrom.ToString() + "," + iNewTo.ToString();
FromToLine.Add(sNewFromTo);
}
}
}
LineIDList.Clear();
FromList.Clear();
ToList.Clear();
pLongArray.RemoveAll();
pCadUndoRedo.StartUndoRedoSession("Adjust Traverse");
if (iLoops == 0)
{
//re-sequence using TraceLines function based on either end point, because the order of
//selected construction lines in grid don't control start or end point
FromToLine.Clear();
int iTerminus = -1;
iTracedLines = 0;
iBranches = 0;
iLoops = 0; iTerminals = 0;
FindTerminusForSequence(ref pFwdStar, iFirstToNode, ref iTerminus, 0);
if (iTerminus == -1)
{
pCadUndoRedo.WriteUndoRedoSession(false);
return;
}
TraceLines(ref pFwdStar, iTerminus, ref iBranches, ref iTracedLines, ref iLoops, ref iTerminals,
ref FromToLine, ref FromList, ref ToList, 0);
}
List <IVector3D> SequencedTraverse = new List <IVector3D>();
IGSLine pGSLineInPath = null;
foreach (string s in FromToLine)
{
string[] sFromTo = s.Split(',');
int iFrom = Convert.ToInt32(sFromTo[0]);
int iTo = Convert.ToInt32(sFromTo[1]);
bool bReversed = pFwdStar.GetLine(iFrom, iTo, ref pGSLineInPath);
if (bReversed)
{
IGSLine pGSLine180 = new GSLineClass();
pGSLine180.FromPoint = pGSLineInPath.ToPoint;
pGSLine180.ToPoint = pGSLineInPath.FromPoint;
pGSLine180.Bearing = pGSLineInPath.Bearing + Math.PI;
pGSLine180.Distance = pGSLineInPath.Distance;
IVector3D vec180 = new Vector3DClass();
vec180.PolarSet(pGSLine180.Bearing, 0, pGSLine180.Distance); //Azimuth of IVector3D is north azimuth radians zero degrees north
Traverse.Add(vec180);
lstPointIds.Add(pGSLine180.ToPoint);
pPointCalc.AddLine(pGSLine180);
}
else
{
double dBear = pGSLineInPath.Bearing;
double dDist = pGSLineInPath.Distance;
IVector3D vec = new Vector3DClass();
vec.PolarSet(dBear, 0, dDist); //Azimuth of IVector3D is north azimuth radians zero degrees north
Traverse.Add(vec);
lstPointIds.Add(pGSLineInPath.ToPoint);
pPointCalc.AddLine(pGSLineInPath);
}
if (pStartPoint == null)
{
if (bReversed)
{
pStartPoint = pCadastralPts.GetPoint(pGSLineInPath.ToPoint);
}
else
{
pStartPoint = pCadastralPts.GetPoint(pGSLineInPath.FromPoint);
}
}
if (bReversed)
{
pClosingPoint = pCadastralPts.GetPoint(pGSLineInPath.FromPoint);
}
else
{
pClosingPoint = pCadastralPts.GetPoint(pGSLineInPath.ToPoint);
}
}
}
if (pStartPoint == null)
{
pCadUndoRedo.WriteUndoRedoSession(false);
return;
}
IPoint pStart = new PointClass();
pStart.X = pStartPoint.X;
pStart.Y = pStartPoint.Y;
string sAdjustMethod = "Compass";
esriParcelAdjustmentType eAdjMethod = esriParcelAdjustmentType.esriParcelAdjustmentCompass;
if (pCadastralEditorSettings2.ParcelAdjustment == esriParcelAdjustmentType.esriParcelAdjustmentNone ||
pCadastralEditorSettings2.ParcelAdjustment == esriParcelAdjustmentType.esriParcelAdjustmentCompass)
{
eAdjMethod = esriParcelAdjustmentType.esriParcelAdjustmentCompass;
}
else if (pCadastralEditorSettings2.ParcelAdjustment == esriParcelAdjustmentType.esriParcelAdjustmentCrandall)
{
sAdjustMethod = "Crandall";
eAdjMethod = pCadastralEditorSettings2.ParcelAdjustment;
}
else if (pCadastralEditorSettings2.ParcelAdjustment == esriParcelAdjustmentType.esriParcelAdjustmentTransit)
{
sAdjustMethod = "Transit";
eAdjMethod = pCadastralEditorSettings2.ParcelAdjustment;
}
pPointCalc.CalculatePoints(eAdjMethod, pStartPoint.Id, pStartPoint, pClosingPoint.Id, pClosingPoint, true);
ITraverseClosure pClose = pPointCalc.Closure;
List <string> lstCoursesFromTo = new List <string>();
List <IVector3D> AdjustedTraverse = new List <IVector3D>();
double dAdjustedPointX = 0; double dAdjustedPointY = 0;
double dPreviousPointX = 0; double dPreviousPointY = 0;
for (int i = 0; i < pClose.CourseCount; i++)
{
IGSPoint pPt = pCadastralPts.GetPoint(lstPointIds[i]);
dAdjustedPointY = pPointCalc.GetCalculatedPoint(lstPointIds[i], ref dAdjustedPointX);
string sFromTo = "";
IVector3D pAdjustedLine = new Vector3DClass();
if (i == 0)
{
sFromTo = pStartPoint.Id.ToString() + "-" + lstPointIds[i].ToString();
pAdjustedLine.SetComponents(dAdjustedPointX - pStartPoint.X, dAdjustedPointY - pStartPoint.Y, 0);
}
else
{
sFromTo = lstPointIds[i - 1].ToString() + "-" + lstPointIds[i].ToString();
pAdjustedLine.SetComponents(dAdjustedPointX - dPreviousPointX, dAdjustedPointY - dPreviousPointY, 0);
}
lstCoursesFromTo.Add(sFromTo);
IVector3D Z_Axis = new Vector3DClass();
Z_Axis.SetComponents(0, 0, 100);
pAdjustedLine.Rotate(TheRotation, Z_Axis);
AdjustedTraverse.Add(pAdjustedLine);
dPreviousPointX = dAdjustedPointX;
dPreviousPointY = dAdjustedPointY;
pPt.X = dAdjustedPointX;
pPt.Y = dAdjustedPointY;
if (!pCadastralEditorSettings2.MeasurementView)
{
pFixedPoints.SetFixedPoint(lstPointIds[i], true);
}
}
double dMisclosureDistance = pClose.MisclosureDistance; double dMisclosureBearing = pClose.MisclosureDirection;
IVector MiscloseVector = new Vector3DClass();
IEditProperties2 pEdProps = pEd as IEditProperties2;
IAngularConverter pAngConv = new AngularConverterClass();
pAngConv.SetAngle(dMisclosureBearing, esriDirectionType.esriDTNorthAzimuth, esriDirectionUnits.esriDURadians);
//int iPrec = 7;
//if (pConstr.Parcel.Plan.AngleUnits == esriDirectionUnits.esriDUDegreesMinutesSeconds)
// iPrec = 0;
string sMiscloseBearing = pAngConv.GetString(pEdProps.DirectionType, pEdProps.DirectionUnits, pEdProps.AngularUnitPrecision);
Utilities UTIL = new Utilities();
string sRatio = "High Accuracy";
if (pClose.RelativeErrorRatio < 10000)
{
sRatio = "1:" + pClose.RelativeErrorRatio.ToString("0");
}
if (dMisclosureDistance >= 0.001)
{
sMiscloseBearing = UTIL.FormatDirectionDashesToDegMinSecSymbols(sMiscloseBearing);
}
else
{
sMiscloseBearing = "----";
}
ICadastralUnitConversion pCadUnitConverter = new CadastralUnitConversionClass();
double dMetersPerUnit = pCadUnitConverter.ConvertDouble(1, pConstr.Parcel.Plan.DistanceUnits, esriCadastralDistanceUnits.esriCDUMeter);
string sReport = "Closure:" + Environment.NewLine +
" error: " + sRatio + Environment.NewLine +
" distance: " + (dMisclosureDistance / dMetersPerUnit).ToString("0.000") + Environment.NewLine +
" bearing: " + sMiscloseBearing + Environment.NewLine +
" xdist: " + (pClose.MisclosureX / dMetersPerUnit).ToString("0.000") + Environment.NewLine +
" ydist: " + (pClose.MisclosureY / dMetersPerUnit).ToString("0.000") + Environment.NewLine +
" courses: " + (pClose.CourseCount) + Environment.NewLine +
Environment.NewLine + "Adjustment:" + Environment.NewLine +
" method: " + sAdjustMethod;
dlgTraverseResults dlgTraverseResults = new dlgTraverseResults();
AddTraverseInfoToGrid(pConstr.Parcel.Plan, dlgTraverseResults.dataGridView1, Traverse, AdjustedTraverse, lstCoursesFromTo);
dlgTraverseResults.txtMiscloseReport.Text = sReport;
DialogResult dRes = dlgTraverseResults.ShowDialog();
if (dRes == DialogResult.Cancel)
{
//since we cancelled, set the points back
foreach (int i in lstPointIds)
{
pFixedPoints.SetFixedPoint(i, false);
}
pCadUndoRedo.WriteUndoRedoSession(false);
}
else
{
pCadUndoRedo.WriteUndoRedoSession(true);
}
}
catch (Exception ex)
{
MessageBox.Show(ex.Message + Environment.NewLine + "Line number:" + ex.LineNumber().ToString()
+ " in " + ex.TargetSite.Name, "Traverse");
pCadUndoRedo.WriteUndoRedoSession(false);
}
}
/// <summary>
/// 构造函数
/// </summary>
/// <param name="graphicsLayer">图层</param>
/// <param name="kmlCircle">圆的kml</param>
public Circle_ArcGlobe(IGlobeGraphicsLayer _graphicsLayer, KmlCircle kmlCircle)
{
this.ElementType = Core.Model.ElementTypeEnum.Circle;
graphicsLayer = _graphicsLayer;
lineSymbol = new SimpleLineSymbolClass();
lineSymbol.Color = new RgbColorClass()
{
Red = kmlCircle.StrokeColor.R,
Green = kmlCircle.StrokeColor.G,
Blue = kmlCircle.StrokeColor.B
};
if (kmlCircle.StrokeWidth == 0)
{
kmlCircle.StrokeWidth = 2;
}
lineSymbol.Width = kmlCircle.StrokeWidth;
fillSymbol = new SimpleFillSymbolClass();
fillSymbol.Outline = lineSymbol;
fillSymbol.Color = new RgbColorClass()
{
Red = kmlCircle.FillColor.R,
Green = kmlCircle.FillColor.G,
Blue = kmlCircle.FillColor.B
};
radius = kmlCircle.Radius;
outlineColor = kmlCircle.StrokeColor;
fillColor = kmlCircle.FillColor;
centerPoint = new PointClass();//圆心坐标
centerPoint.PutCoords(kmlCircle.Position.Lng, kmlCircle.Position.Lat);
centerPoint.Z = kmlCircle.Position.Alt;
missing = System.Type.Missing;
IZAware zAware = (IGeometry)centerPoint as IZAware;
zAware.ZAware = true;
upperAxisVector3D = new Vector3DClass();
upperAxisVector3D.SetComponents(0, 0, 2);
lowerAxisVector3D = new Vector3DClass();
lowerAxisVector3D.SetComponents(0, 0, -2);
lowerAxisVector3D.XComponent -= vectorComponentOffset;
lowerAxisVector3D.YComponent -= vectorComponentOffset;//TODO
normalVector3D = upperAxisVector3D.CrossProduct(lowerAxisVector3D) as IVector3D;
normalVector3D.Magnitude = kmlCircle.Radius;
double rotationAngleInRadians = 2 * (Math.PI / 180);
//geometryCollection = new MultiPatchClass();
pointCollection = new PolygonClass();
for (int i = 0; i < 180; i++)
{
normalVector3D.Rotate(-1 * rotationAngleInRadians, upperAxisVector3D);
IPoint vertexPoint = new PointClass();
vertexPoint.X = centerPoint.X + normalVector3D.XComponent;
vertexPoint.Y = centerPoint.Y + normalVector3D.YComponent;
vertexPoint.Z = centerPoint.Z;
pointCollection.AddPoint(vertexPoint, missing, missing);
}
base.Symbol = fillSymbol;
base.Geometry = pointCollection as IGeometry;
flashTimer = new System.Timers.Timer();
flashTimer.Elapsed += new System.Timers.ElapsedEventHandler(flashTimer_Elapsed);
flashTimer.Interval = 500;
}
private IVector GetClosingVector(List<IVector3D> TraverseCourses, IPoint StartPoint, IPoint EndPoint, out double SUMofLengths)
{
IVector SumVec = null;
SUMofLengths = 0;
for (int i = 0; i < TraverseCourses.Count - 1; i++)
{
if (i == 0)
{
SUMofLengths = TraverseCourses[0].Magnitude + TraverseCourses[1].Magnitude;
SumVec = TraverseCourses[0].AddVector(TraverseCourses[1]);
}
else
{
IVector3D SumVec3D = SumVec as IVector3D;
SUMofLengths += TraverseCourses[i + 1].Magnitude;
SumVec = SumVec3D.AddVector(TraverseCourses[i + 1]);
}
}
double dCalcedEndX = StartPoint.X + SumVec.ComponentByIndex[0];
double dCalcedEndY = StartPoint.Y + SumVec.ComponentByIndex[1];
IVector3D CloseVector3D = new Vector3DClass();
CloseVector3D.SetComponents(dCalcedEndX - EndPoint.X, dCalcedEndY - EndPoint.Y, 0);
IVector CloseVector = CloseVector3D as IVector;
return CloseVector;
}
protected override void OnClick()
{
ICadastralEditor pCadEd = (ICadastralEditor)ArcMap.Application.FindExtensionByName("esriCadastralUI.CadastralEditorExtension");
IParcelEditManager pParcEditorMan = (IParcelEditManager)pCadEd;
ICadastralPacketManager pCadPacketMan = (ICadastralPacketManager)pCadEd;
//bool bStartedWithPacketOpen = pCadPacketMan.PacketOpen;
if (pParcEditorMan == null)
return;
IEditor pEd = (IEditor)ArcMap.Application.FindExtensionByName("esri object editor");
if (pEd.EditState == esriEditState.esriStateNotEditing)
{
MessageBox.Show("Please start editing and try again.");
return;
}
IParcelConstruction pConstr = pParcEditorMan.ParcelConstruction;
IParcelConstruction4 pConstr4 = pConstr as IParcelConstruction4;
ICadastralPoints pCadastralPts = pConstr4 as ICadastralPoints;
ICadastralEditorSettings2 pCadastralEditorSettings2 = pCadEd as ICadastralEditorSettings2;
ICadastralFixedPoints pFixedPoints = pCadastralPts as ICadastralFixedPoints;
IPointCalculation pPointCalc = new PointCalculationClass();
if (pConstr == null)
return;
IGSLine pParcelLine = null;
IMetricUnitConverter pMetricUnitConv = (IMetricUnitConverter)pCadEd;
IGSPoint pStartPoint = null;
List<int> lstPointIds = new List<int>();
List<IVector3D> Traverse = new List<IVector3D>();
//get rotation here
IParcelConstructionData pConstrData = pConstr4.ConstructionData;
IConstructionParentParcels pConstructionParentParcels = pConstrData as IConstructionParentParcels;
ICadastralUndoRedo pCadUndoRedo = pConstr as ICadastralUndoRedo;
try
{
int iParcelID = -1;
if (pConstructionParentParcels.ParentParcelCount>0)
pConstructionParentParcels.GetParentParcel(0, ref iParcelID);
ICadastralParcel pCadaParcel = pCadPacketMan.JobPacket as ICadastralParcel;
IGSParcel pGSParcel = null;
if (pCadaParcel != null)
pGSParcel = pCadaParcel.GetParcel(iParcelID);
//if in measurement view then rotation is 0
double TheRotation = 0;
if (pGSParcel == null)
pGSParcel = pConstr.Parcel;
if (!pCadastralEditorSettings2.MeasurementView)
TheRotation = pGSParcel.Rotation;//radians
if (TheRotation == 123456789)
TheRotation = 0;
pPointCalc.Rotation = TheRotation;
IGSPoint pClosingPoint = null;
#region simple method as fall-back
bool bUseSimpleStackSelection = false;
if (bUseSimpleStackSelection)
{
//bool bLineSelectionSequence = false;
//IGSLine pLastSelectedGSLineInGrid = null;
//for (int i = 0; i < pConstr.LineCount; i++)
//{
// if (pConstr.GetLineSelection(i))
// {
// if (pConstr.GetLine(i, ref pParcelLine))
// {
// if (!bLineSelectionSequence) //first line
// {
// pStartPoint = pCadastralPts.GetPoint(pParcelLine.FromPoint);
// pToPoint = pCadastralPts.GetPoint(pParcelLine.ToPoint);
// }
// pPointCalc.AddLine(pParcelLine);
// pLastSelectedGSLineInGrid = pParcelLine;
// pToPoint = pCadastralPts.GetPoint(pParcelLine.ToPoint);
// lstPointIds.Add(pToPoint.Id);
// }
// bLineSelectionSequence = true;
// double dBear = pParcelLine.Bearing; //Azimuth of IVector3D is north azimuth radians zero degrees north
// double dDist = pParcelLine.Distance;
// IVector3D vec = new Vector3DClass();
// vec.PolarSet(dBear, 0, dDist); ////Azimuth of IVector3D is north azimuth radians zero degrees north
// Traverse.Add(vec);
// }
// else
// {
// if (bLineSelectionSequence && pConstr.GetLine(i, ref pParcelLine) && HasLineSelectionAfter(pConstr, i))
// //if there was a prior selection and this line is a complete line, and there is no later selection
// {
// MessageBox.Show("Please select a continuous set of lines for closure.");
// return;
// }
// }
//}
//pClosingPoint = pCadastralPts.GetPoint(pLastSelectedGSLineInGrid.ToPoint);
}
else
#endregion
{//build a forward star for the selected lines
IEnumCELines pCELines = new EnumCELinesClass();
IEnumGSLines pEnumGSLines = (IEnumGSLines)pCELines;
ILongArray pLongArray = new LongArrayClass();
int iFirstToNode = -1;
for (int i = 0; i < pConstr.LineCount; i++)
{
if (pConstr.GetLineSelection(i))
{
if (pConstr.GetLine(i, ref pParcelLine))
{
if (iFirstToNode < 0)
iFirstToNode = pParcelLine.ToPoint;
pLongArray.Add(i);
pCELines.Add(pParcelLine);
}
}
}
if (pCELines.Count==0)
{
MessageBox.Show("No lines selected. Please select a continuous set of lines for closure." + Environment.NewLine +
"Line selection should not have branches.", "Traverse");
return;
}
IParcelLineFunctions3 ParcelLineFx = new ParcelFunctionsClass();
IGSForwardStar pFwdStar = ParcelLineFx.CreateForwardStar(pEnumGSLines);
//forward star object is now created for all the selected lines,
//need to first re-sequence the lines, and test for branching and discontinuity
int iBranches = 0; int iTracedLines = 0;
int iLoops = 0; int iTerminals = 0;
List<int> LineIDList = new List<int>();
List<int> FromList = new List<int>();
List<int> ToList = new List<int>();
List<string> FromToLine = new List<string>();
bool bTraceSucceeded = TraceLines(ref pFwdStar, iFirstToNode, ref iBranches, ref iTracedLines,
ref iLoops, ref iTerminals, ref FromToLine, ref FromList, ref ToList, 0);
if (iBranches > 0)
{
MessageBox.Show("Please select a continuous set of lines for closure." + Environment.NewLine +
"Line selection should not have branches.", "Traverse");
return;
}
if (iTracedLines < pLongArray.Count)
{
MessageBox.Show("Please select a continuous set of lines for closure." + Environment.NewLine +
"Selected Lines should be connected in a single sequence without branches.", "Traverse");
return;
}
//if it's a single loop check to see if the sequence needs to be reversed
//CW or CCW based on bearings
if (iLoops == 1)
{
bool bIsReversed = false;
foreach (int i in FromList)
{
if (i < 0)
bIsReversed = true;
else
{
bIsReversed = false;
break;
}
}
if (bIsReversed)
{//all courses are running reversed, so reverse the whole sequence
FromToLine.Clear();
FromList.Reverse();
ToList.Reverse();
int iNewFrom = -ToList[ToList.Count-1];
int iNewTo = -FromList[ToList.Count-1];
string sNewFromTo = iNewFrom.ToString() + "," + iNewTo.ToString();
FromToLine.Add(sNewFromTo);
for (int i =1; i < ToList.Count ;i++)
{
iNewFrom = -ToList[i-1];
iNewTo = -FromList[i-1];
sNewFromTo = iNewFrom.ToString() + "," + iNewTo.ToString();
FromToLine.Add(sNewFromTo);
}
}
}
LineIDList.Clear();
FromList.Clear();
ToList.Clear();
pLongArray.RemoveAll();
pCadUndoRedo.StartUndoRedoSession("Adjust Traverse");
if (iLoops == 0)
{
//re-sequence using TraceLines function based on either end point, because the order of
//selected construction lines in grid don't control start or end point
FromToLine.Clear();
int iTerminus = -1;
iTracedLines = 0;
iBranches = 0;
iLoops = 0; iTerminals = 0;
FindTerminusForSequence(ref pFwdStar, iFirstToNode, ref iTerminus, 0);
if (iTerminus == -1)
{
pCadUndoRedo.WriteUndoRedoSession(false);
return;
}
TraceLines(ref pFwdStar, iTerminus, ref iBranches, ref iTracedLines, ref iLoops, ref iTerminals,
ref FromToLine, ref FromList, ref ToList, 0);
}
List<IVector3D> SequencedTraverse = new List<IVector3D>();
IGSLine pGSLineInPath = null;
foreach (string s in FromToLine)
{
string[] sFromTo = s.Split(',');
int iFrom = Convert.ToInt32(sFromTo[0]);
int iTo = Convert.ToInt32(sFromTo[1]);
bool bReversed = pFwdStar.GetLine(iFrom, iTo, ref pGSLineInPath);
if (bReversed)
{
IGSLine pGSLine180 = new GSLineClass();
pGSLine180.FromPoint = pGSLineInPath.ToPoint;
pGSLine180.ToPoint = pGSLineInPath.FromPoint;
pGSLine180.Bearing = pGSLineInPath.Bearing + Math.PI;
pGSLine180.Distance = pGSLineInPath.Distance;
IVector3D vec180 = new Vector3DClass();
vec180.PolarSet(pGSLine180.Bearing, 0, pGSLine180.Distance); //Azimuth of IVector3D is north azimuth radians zero degrees north
Traverse.Add(vec180);
lstPointIds.Add(pGSLine180.ToPoint);
pPointCalc.AddLine(pGSLine180);
}
else
{
double dBear = pGSLineInPath.Bearing;
double dDist = pGSLineInPath.Distance;
IVector3D vec = new Vector3DClass();
vec.PolarSet(dBear, 0, dDist); //Azimuth of IVector3D is north azimuth radians zero degrees north
Traverse.Add(vec);
lstPointIds.Add(pGSLineInPath.ToPoint);
pPointCalc.AddLine(pGSLineInPath);
}
if (pStartPoint == null)
{
if (bReversed)
pStartPoint = pCadastralPts.GetPoint(pGSLineInPath.ToPoint);
else
pStartPoint = pCadastralPts.GetPoint(pGSLineInPath.FromPoint);
}
if (bReversed)
pClosingPoint = pCadastralPts.GetPoint(pGSLineInPath.FromPoint);
else
pClosingPoint = pCadastralPts.GetPoint(pGSLineInPath.ToPoint);
}
}
if (pStartPoint == null)
{
pCadUndoRedo.WriteUndoRedoSession(false);
return;
}
IPoint pStart = new PointClass();
pStart.X = pStartPoint.X;
pStart.Y = pStartPoint.Y;
string sAdjustMethod = "Compass";
esriParcelAdjustmentType eAdjMethod = esriParcelAdjustmentType.esriParcelAdjustmentCompass;
if (pCadastralEditorSettings2.ParcelAdjustment == esriParcelAdjustmentType.esriParcelAdjustmentNone ||
pCadastralEditorSettings2.ParcelAdjustment == esriParcelAdjustmentType.esriParcelAdjustmentCompass)
eAdjMethod = esriParcelAdjustmentType.esriParcelAdjustmentCompass;
else if (pCadastralEditorSettings2.ParcelAdjustment == esriParcelAdjustmentType.esriParcelAdjustmentCrandall)
{
sAdjustMethod = "Crandall";
eAdjMethod = pCadastralEditorSettings2.ParcelAdjustment;
}
else if (pCadastralEditorSettings2.ParcelAdjustment == esriParcelAdjustmentType.esriParcelAdjustmentTransit)
{
sAdjustMethod = "Transit";
eAdjMethod = pCadastralEditorSettings2.ParcelAdjustment;
}
pPointCalc.CalculatePoints(eAdjMethod, pStartPoint.Id, pStartPoint, pClosingPoint.Id, pClosingPoint, true);
ITraverseClosure pClose = pPointCalc.Closure;
List<string> lstCoursesFromTo = new List<string>();
List<IVector3D> AdjustedTraverse = new List<IVector3D>();
double dAdjustedPointX = 0; double dAdjustedPointY = 0;
double dPreviousPointX = 0; double dPreviousPointY = 0;
for (int i = 0; i < pClose.CourseCount; i++)
{
IGSPoint pPt = pCadastralPts.GetPoint(lstPointIds[i]);
dAdjustedPointY = pPointCalc.GetCalculatedPoint(lstPointIds[i], ref dAdjustedPointX);
string sFromTo = "";
IVector3D pAdjustedLine = new Vector3DClass();
if (i == 0)
{
sFromTo = pStartPoint.Id.ToString() + "-" + lstPointIds[i].ToString();
pAdjustedLine.SetComponents(dAdjustedPointX - pStartPoint.X, dAdjustedPointY - pStartPoint.Y, 0);
}
else
{
sFromTo = lstPointIds[i - 1].ToString() + "-" + lstPointIds[i].ToString();
pAdjustedLine.SetComponents(dAdjustedPointX - dPreviousPointX, dAdjustedPointY - dPreviousPointY, 0);
}
lstCoursesFromTo.Add(sFromTo);
IVector3D Z_Axis = new Vector3DClass();
Z_Axis.SetComponents(0, 0, 100);
pAdjustedLine.Rotate(TheRotation, Z_Axis);
AdjustedTraverse.Add(pAdjustedLine);
dPreviousPointX = dAdjustedPointX;
dPreviousPointY = dAdjustedPointY;
pPt.X = dAdjustedPointX;
pPt.Y = dAdjustedPointY;
if (!pCadastralEditorSettings2.MeasurementView)
pFixedPoints.SetFixedPoint(lstPointIds[i], true);
}
double dMisclosureDistance = pClose.MisclosureDistance; double dMisclosureBearing = pClose.MisclosureDirection;
IVector MiscloseVector = new Vector3DClass();
IEditProperties2 pEdProps = pEd as IEditProperties2;
IAngularConverter pAngConv = new AngularConverterClass();
pAngConv.SetAngle(dMisclosureBearing, esriDirectionType.esriDTNorthAzimuth, esriDirectionUnits.esriDURadians);
//int iPrec = 7;
//if (pConstr.Parcel.Plan.AngleUnits == esriDirectionUnits.esriDUDegreesMinutesSeconds)
// iPrec = 0;
string sMiscloseBearing = pAngConv.GetString(pEdProps.DirectionType, pEdProps.DirectionUnits, pEdProps.AngularUnitPrecision);
Utilities UTIL = new Utilities();
string sRatio = "High Accuracy";
if (pClose.RelativeErrorRatio < 10000)
sRatio = "1:" + pClose.RelativeErrorRatio.ToString("0");
if (dMisclosureDistance >= 0.001)
sMiscloseBearing = UTIL.FormatDirectionDashesToDegMinSecSymbols(sMiscloseBearing);
else
sMiscloseBearing = "----";
ICadastralUnitConversion pCadUnitConverter = new CadastralUnitConversionClass();
double dMetersPerUnit = pCadUnitConverter.ConvertDouble(1, pConstr.Parcel.Plan.DistanceUnits, esriCadastralDistanceUnits.esriCDUMeter);
string sReport = "Closure:" + Environment.NewLine +
" error: " + sRatio + Environment.NewLine +
" distance: " + (dMisclosureDistance / dMetersPerUnit).ToString("0.000") + Environment.NewLine +
" bearing: " + sMiscloseBearing + Environment.NewLine +
" xdist: " + (pClose.MisclosureX / dMetersPerUnit).ToString("0.000") + Environment.NewLine +
" ydist: " + (pClose.MisclosureY / dMetersPerUnit).ToString("0.000") + Environment.NewLine +
" courses: " + (pClose.CourseCount) + Environment.NewLine +
Environment.NewLine + "Adjustment:" + Environment.NewLine +
" method: " + sAdjustMethod;
dlgTraverseResults dlgTraverseResults = new dlgTraverseResults();
AddTraverseInfoToGrid(pConstr.Parcel.Plan, dlgTraverseResults.dataGridView1, Traverse, AdjustedTraverse, lstCoursesFromTo);
dlgTraverseResults.txtMiscloseReport.Text = sReport;
DialogResult dRes = dlgTraverseResults.ShowDialog();
if (dRes == DialogResult.Cancel)
{
//since we cancelled, set the points back
foreach (int i in lstPointIds)
pFixedPoints.SetFixedPoint(i, false);
pCadUndoRedo.WriteUndoRedoSession(false);
}
else
pCadUndoRedo.WriteUndoRedoSession(true);
}
catch (Exception ex)
{
MessageBox.Show(ex.Message + Environment.NewLine + "Line number:" + ex.LineNumber().ToString()
+ " in " + ex.TargetSite.Name, "Traverse");
pCadUndoRedo.WriteUndoRedoSession(false);
}
}