private void AddRectangle(double LeftmostX, double LowestY, double RightmostX, double HighestY)
{
MapWinGIS.Shape newPolygon = new MapWinGIS.Shape();
int partIndex = 0;
newPolygon.Create(MapWinGIS.ShpfileType.SHP_POLYGON);
newPolygon.InsertPart(0, ref partIndex);
int shpIndex = 0;
MapWinGIS.Point newPt = new MapWinGIS.Point();
newPt.x = LeftmostX;
newPt.y = HighestY;
newPt.Z = 0;
newPolygon.InsertPoint(newPt, ref shpIndex);
newPt = new MapWinGIS.Point();
newPt.x = RightmostX;
newPt.y = HighestY;
newPt.Z = 0;
newPolygon.InsertPoint(newPt, ref shpIndex);
newPt = new MapWinGIS.Point();
newPt.x = RightmostX;
newPt.y = LowestY;
newPt.Z = 0;
newPolygon.InsertPoint(newPt, ref shpIndex);
newPt = new MapWinGIS.Point();
newPt.x = LeftmostX;
newPt.y = LowestY;
newPt.Z = 0;
newPolygon.InsertPoint(newPt, ref shpIndex);
// Finalize -- add it.
MapWinGIS.Shapefile sf = (MapWinGIS.Shapefile)g.MapWin.Layers[g.MapWin.Layers.CurrentLayer].GetObject();
if (!sf.EditingShapes || !sf.EditingTable)
{
sf.StartEditingShapes(true, null);
}
int addedShapes = sf.NumShapes;
sf.EditInsertShape(newPolygon, ref addedShapes);
g.CreateUndoPoint();
sf.StopEditingShapes(true, true, null);
// Release memory used by point reallocations above
GC.Collect();
// And show it:
g.UpdateView();
g.MapWin.Plugins.BroadcastMessage("ShapefileEditor: Layer " + g.MapWin.Layers.CurrentLayer.ToString() + ": New Shape Added");
}
private void AddEllipse(double LeftmostX, double LowestY, double RightmostX, double HighestY)
{
MapWinGIS.Shape newPolygon = new MapWinGIS.Shape();
int partIndex = 0;
newPolygon.Create(MapWinGIS.ShpfileType.SHP_POLYGON);
newPolygon.InsertPart(0, ref partIndex);
int shpIndex = 0;
double t, a, b, tinc, centx, centy;
a = RightmostX - LeftmostX;
b = HighestY - LowestY;
tinc = Math.PI * 2 / (a + b);
centx = (LeftmostX + RightmostX) * .5;
centy = (LowestY + HighestY) * .5;
MapWinGIS.Point newPt = new MapWinGIS.Point();
newPt.x = centx + a;
newPt.y = centy;
newPt.Z = 0;
newPolygon.InsertPoint(newPt, ref shpIndex);
for (t = 0; t < Math.PI * 2; t += tinc)
{
MapWinGIS.Point nextPt = new MapWinGIS.Point();
nextPt.x = centx + a * Math.Cos(t);
nextPt.y = centy - b * Math.Sin(t);
nextPt.Z = 0;
newPolygon.InsertPoint(nextPt, ref shpIndex);
}
// Finalize -- add it.
MapWinGIS.Shapefile sf = (MapWinGIS.Shapefile)g.MapWin.Layers[g.MapWin.Layers.CurrentLayer].GetObject();
if (!sf.EditingShapes || !sf.EditingTable)
{
sf.StartEditingShapes(true, null);
}
int addedShapes = sf.NumShapes;
sf.EditInsertShape(newPolygon, ref addedShapes);
g.CreateUndoPoint();
sf.StopEditingShapes(true, true, null);
// And show it:
g.UpdateView();
g.MapWin.Plugins.BroadcastMessage("ShapefileEditor: Layer " + g.MapWin.Layers.CurrentLayer.ToString() + ": New Shape Added");
}
/// <summary>
/// Returns a MapWinGIS Multi-part PolyLine representing this multiLineString
/// </summary>
/// <returns>MapWinGIS Shape</returns>
public static MapWinGIS.Shape mwShapeFromMultiLineString(MultiLineString MLS, MapWinGIS.ShpfileType sfType)
{
MapWinGIS.Shape mwShape = new MapWinGIS.Shape();
mwShape.Create(sfType);
int prt = 0; // Part Index
int pt = 0; // Point Index
int numParts = MLS.NumGeometries;
int maxpt;
while (prt < numParts)
{
LineString LS = MLS.GeometryN[prt] as LineString;
maxpt = LS.Coordinates.Count;
mwShape.InsertPart(pt, ref prt);
for (int I = 0; I < maxpt; I++)
{
MapWinGIS.Point mwPoint = new MapWinGIS.Point();
mwPoint.x = LS.Coordinates[I].X;
mwPoint.y = LS.Coordinates[I].Y;
bool result = mwShape.InsertPoint(mwPoint, ref pt);
if (result == false)
{
throw new ApplicationException(mwShape.get_ErrorMsg(mwShape.LastErrorCode));
}
pt++;
}
prt++;
}
return(mwShape);
}
/// <summary>
/// Converts a list of polylines into a single multi-part shape
/// </summary>
/// <param name="Polylines"></param>
/// <returns></returns>
public static MapWinGIS.Shape LineStrings_To_mwShape(List <LineString> Polylines)
{
MapWinGIS.Shape mwShape = new MapWinGIS.Shape();
bool res;
res = mwShape.Create(MapWinGIS.ShpfileType.SHP_POLYLINEZ);
if (res != true)
{
throw new ApplicationException(mwShape.get_ErrorMsg(mwShape.LastErrorCode));
}
int idx = 0;
for (int iPart = 0; iPart < Polylines.Count; iPart++)
{
LineString Part = Polylines[iPart];
// we only cal set_part if we have multiple polylines
if (Polylines.Count > 0)
{
mwShape.set_Part(iPart, idx);
}
for (int iPoint = 0; iPoint < Part.Points.Count; iPoint++)
{
mwShape.InsertPoint(Part.Points[iPoint].To_mwPoint(), ref idx);
idx++;
}
}
return(mwShape);
} // End LineStrings_To_mwShape
/// <summary>
/// Insers part to polygon based on given rectange
/// </summary>
private void InsertPart(MapWinGIS.Shape shp, double xMin, double xMax, double yMin, double yMax)
{
int numParts = shp.NumParts;
shp.InsertPart(shp.numPoints, ref numParts);
// to left
int index = shp.numPoints;
MapWinGIS.Point pnt = new MapWinGIS.Point();
pnt.x = xMin; pnt.y = yMax;
shp.InsertPoint(pnt, ref index); index++;
pnt = new MapWinGIS.Point();
pnt.x = xMax; pnt.y = yMax;
shp.InsertPoint(pnt, ref index); index++;
pnt = new MapWinGIS.Point();
pnt.x = xMax; pnt.y = yMin;
shp.InsertPoint(pnt, ref index); index++;
pnt = new MapWinGIS.Point();
pnt.x = xMin; pnt.y = yMin;
shp.InsertPoint(pnt, ref index); index++;
pnt = new MapWinGIS.Point();
pnt.x = xMin; pnt.y = yMax;
shp.InsertPoint(pnt, ref index); index++;
if (!shp.get_PartIsClockWise(0))
{
bool val = shp.ReversePointsOrder(0);
if (!val)
{
System.Diagnostics.Debug.Print("CCW");
}
}
if (!shp.get_PartIsClockWise(0))
{
System.Diagnostics.Debug.Print("CCW");
}
}
private void AddRegularPolygon(int sides, double sidelength, double centroidX, double centroidY)
{
double pi = 3.1415926535897932384626433832795;
double apothem = sidelength / (2 * Math.Tan(pi / sides));
double radius = sidelength / (2 * Math.Sin(pi / sides));
MapWinGIS.Shape newPolygon = new MapWinGIS.Shape();
int partIndex = 0;
newPolygon.Create(MapWinGIS.ShpfileType.SHP_POLYGON);
newPolygon.InsertPart(0, ref partIndex);
double i = 0;
int shpIndex = 0;
// Approximate the polygon with n line segments
double increment = ((2 * pi) / sides);
for (; i <= 2 * pi; i += increment, shpIndex++)
{
MapWinGIS.Point newPt = new MapWinGIS.Point();
newPt.x = (radius * Math.Cos(i) + centroidX);
newPt.y = (radius * Math.Sin(i) + centroidY);
newPt.Z = 0;
newPolygon.InsertPoint(newPt, ref shpIndex);
}
// Finalize -- add it.
MapWinGIS.Shapefile sf = (MapWinGIS.Shapefile)g.MapWin.Layers[g.MapWin.Layers.CurrentLayer].GetObject();
if (!sf.EditingShapes || !sf.EditingTable)
{
sf.StartEditingShapes(true, null);
}
int addedShapes = sf.NumShapes;
sf.EditInsertShape(newPolygon, ref addedShapes);
g.CreateUndoPoint();
sf.StopEditingShapes(true, true, null);
// And show it:
g.UpdateView();
g.MapWin.Plugins.BroadcastMessage("ShapefileEditor: Layer " + g.MapWin.Layers.CurrentLayer.ToString() + ": New Shape Added");
}
private bool AddPointToShapeFile(double x, double y, MapWinGIS.Shapefile shpFile, int shpIndex, int insertIndex)
{
try
{
MapWinGIS.Point p;
//start editing shapes
shpFile.StartEditingShapes(true, null);
MapWinGIS.Shape shp = shpFile.get_Shape(shpIndex);
p = new MapWinGIS.Point();
p.x = x;
p.y = y;
if (!shp.InsertPoint(p, ref insertIndex))
{
return(false);
}
else if (shp.NumParts > 1)
{
// Shift up part indices if needed (bugzilla 798)
bool shifting = false;
for (int i = 0; i < shp.NumParts; i++)
{
if (shp.get_Part(i) >= insertIndex)
{
// Shift up this one and all remaining
shifting = true;
}
if (shifting)
{
shp.set_Part(i, shp.get_Part(i) + 1);
}
}
}
m_Globals.CreateUndoPoint();
//stop editing and save changes
shpFile.StopEditingShapes(true, true, null);
return(true);
}
catch (System.Exception ex)
{
m_MapWin.ShowErrorDialog(ex);
}
return(false);
}
/// <summary>
/// Converts a list of 3d-points to a point shapefile with z-value field.
/// This function creates a new shapefile. The shapefile has two fields:
/// a 'MWShapeId' field and a field which contains the z-value.
/// </summary>
/// <param name="ShpFileName">Name of the resulting point shapefile</param>
/// <param name="ZFieldName">Name of the z-field in the shapefile</param>
public void ToShapefile(string ShpFileName, string ZFieldName)
{
MapWinGIS.Shapefile newSF = new MapWinGIS.Shapefile();
try
{
Hashtable FieldIndices = new Hashtable();
MapWinGIS.ShpfileType sftype;
sftype = MapWinGIS.ShpfileType.SHP_POINT;
int fldIdx = 0;
// if shapefile exists - open it and clear all shapes
if (System.IO.File.Exists(ShpFileName))
{
newSF.Open(ShpFileName, null);
newSF.StartEditingShapes(true, null);
newSF.EditClear();
}
else //else, create a new shapefile
{
if (!newSF.CreateNew(ShpFileName, sftype))
{
throw new InvalidOperationException
("Error creating shapefile " + newSF.get_ErrorMsg(newSF.LastErrorCode));
}
newSF.StartEditingShapes(true, null);
}
//check existing fields:
for (int i = 0; i < newSF.NumFields; ++i)
{
MapWinGIS.Field fl = newSF.get_Field(i);
if (fl.Name == "MWShapeID")
{
FieldIndices.Add("MWShapeID", i);
}
if (fl.Name == ZFieldName)
{
FieldIndices.Add(ZFieldName, i);
}
}
//Add the fields:
if (!FieldIndices.ContainsKey("MWShapeID"))
{
//First an ID field
MapWinGIS.Field idFld = new MapWinGIS.Field();
idFld.Name = "MWShapeID";
idFld.Type = MapWinGIS.FieldType.INTEGER_FIELD;
fldIdx = newSF.NumFields;
if (newSF.EditInsertField(idFld, ref fldIdx, null) == false)
{
throw new InvalidOperationException("error inserting field " +
newSF.get_ErrorMsg(newSF.LastErrorCode));
}
FieldIndices.Add("MWShapeID", fldIdx);
}
if (!FieldIndices.ContainsKey(ZFieldName))
{
//Second add a Z-field
MapWinGIS.Field zFld = new MapWinGIS.Field();
zFld.Name = "Z";
zFld.Type = MapWinGIS.FieldType.DOUBLE_FIELD;
fldIdx = newSF.NumFields;
if (newSF.EditInsertField(zFld, ref fldIdx, null) == false)
{
throw new InvalidOperationException("error inserting field " +
newSF.get_ErrorMsg(newSF.LastErrorCode));
}
FieldIndices.Add("Z", fldIdx);
}
foreach (ICoordinate pt in _points)
{
//first, add a point shape (geometry)
MapWinGIS.Shape newShp = new MapWinGIS.Shape();
newShp.Create(MapWinGIS.ShpfileType.SHP_POINT);
MapWinGIS.Point newPt = new MapWinGIS.Point();
newPt.x = pt.X;
newPt.y = pt.Y;
int ptIdx = 0;
newShp.InsertPoint(newPt, ref ptIdx);
int shpIdx = newSF.NumShapes;
newSF.EditInsertShape(newShp, ref shpIdx);
//second add the z-value
newSF.EditCellValue(fldIdx, shpIdx, pt.Z);
}
}
finally
{
//finally stop editing and close the shapefile
newSF.StopEditingShapes(true, true, null);
if (newSF.Close() == false)
{
throw new InvalidOperationException("error closing shapefile " +
newSF.get_ErrorMsg(newSF.LastErrorCode));
}
}
}
private void btnOK_Click(object sender, EventArgs e)
{
double tolerance = 0;
if (!double.TryParse(txtTolerance.Text, out tolerance))
{
MapWinUtility.Logger.Message("Please enter only numbers in the distance field.", "Enter Only Numbers", MessageBoxButtons.OK, MessageBoxIcon.Exclamation, DialogResult.OK);
return;
}
this.Cursor = Cursors.WaitCursor;
prgShape.Visible = true;
prgPoint.Visible = true;
sf.StartEditingShapes(true, null);
long totalRemoved = 0;
long fromShapes = 0;
int numPts;
System.Collections.ArrayList removeList = new System.Collections.ArrayList();
if (sf.NumShapes > 0)
{
prgShape.Maximum = sf.NumShapes;
for (int z = 0; z < sf.NumShapes; z++)
{
MapWinGIS.Shape shp = sf.get_Shape(z);
if (z < 100 || z % 5 == 0)
{
prgShape.Value = z;
this.Refresh();
}
Application.DoEvents();
numPts = shp.numPoints;
if (numPts > 0)
{
prgPoint.Maximum = numPts;
for (int i = numPts - 1; i > 0; i--) // 0 never needs to actually be hit, inner loop will compare against all zeros
{
if (i % 5 == 0)
{
prgPoint.Value = numPts - i;
this.Refresh();
}
for (int j = i - 1; j >= 0; j--)
{
if (Math.Sqrt(Math.Pow(shp.get_Point(i).x - shp.get_Point(j).x, 2) + Math.Pow(shp.get_Point(i).y - shp.get_Point(j).y, 2)) < tolerance)
{
// Make sure that !(i == numPts - 1 && j == 0) -- polygon completion point
if (!removeList.Contains(i) && !(i == numPts - 1 && j == 0))
{
removeList.Add(i);
}
}
}
}
}
if (removeList.Count > 0)
{
if (removeList.Count >= shp.numPoints - 1)
{
// Probably not a good thing.....
MapWinUtility.Logger.Message("Aborting: Proceeding will remove all points from one or more shapes. The distance may need to be smaller, particularly for unprojected (latitute and longitude) coordinate systems.", "Aborting -- All Points Would Be Removed", MessageBoxButtons.OK, MessageBoxIcon.Error, DialogResult.OK);
sf.StopEditingShapes(false, true, null);
this.Cursor = Cursors.Default;
prgPoint.Value = 0;
prgShape.Value = 0;
prgPoint.Visible = false;
prgShape.Visible = false;
this.Cursor = Cursors.Default;
return;
}
totalRemoved += removeList.Count;
fromShapes++;
while (removeList.Count > 0)
{
for (int part = 0; part < shp.NumParts; part++)
{
if (shp.get_Part(part) >= (int)removeList[0])
{
shp.set_Part(part, shp.get_Part(part) - 1);
}
}
shp.DeletePoint((int)removeList[0]);
removeList.RemoveAt(0);
}
// If this is a polygon and there are now less than 3 shapes, panic
if ((shp.ShapeType == MapWinGIS.ShpfileType.SHP_POLYGON || shp.ShapeType == MapWinGIS.ShpfileType.SHP_POLYGONZ || shp.ShapeType == MapWinGIS.ShpfileType.SHP_POLYGONM) && shp.numPoints < 3)
{
// Probably not a good thing.....
MapWinUtility.Logger.Message("Aborting: Proceeding will leave less than 3 points in a polygon. The distance may need to be smaller, particularly for unprojected (latitute and longitude) coordinate systems.", "Aborting -- Polygons Would Be Destroyed", MessageBoxButtons.OK, MessageBoxIcon.Error, DialogResult.OK);
sf.StopEditingShapes(false, true, null);
this.Cursor = Cursors.Default;
prgPoint.Value = 0;
prgShape.Value = 0;
prgPoint.Visible = false;
prgShape.Visible = false;
this.Cursor = Cursors.Default;
return;
}
// If the first and last points are not the same now, reclose it
if (shp.get_Point(0).x != shp.get_Point(shp.numPoints - 1).x || shp.get_Point(0).y != shp.get_Point(shp.numPoints - 1).y)
{
MapWinGIS.Point pnt = new MapWinGIS.Point();
pnt.x = shp.get_Point(0).x;
pnt.y = shp.get_Point(0).y;
pnt.Z = shp.get_Point(0).Z;
int ptidx = shp.numPoints;
shp.InsertPoint(pnt, ref ptidx);
}
}
}
}
prgPoint.Value = prgPoint.Maximum;
prgShape.Value = prgShape.Maximum;
g.CreateUndoPoint();
sf.StopEditingShapes(true, true, null);
this.Cursor = Cursors.Default;
if (totalRemoved > 0)
{
MapWinUtility.Logger.Message("There were " + totalRemoved.ToString() + " points removed from " + fromShapes.ToString() + " shapes.", "Finished", MessageBoxButtons.OK, MessageBoxIcon.Information, DialogResult.OK);
}
else
{
MapWinUtility.Logger.Message("Finished -- no extra points needed to be removed.", "Finished", MessageBoxButtons.OK, MessageBoxIcon.Information, DialogResult.OK);
}
this.DialogResult = DialogResult.OK;
this.Close();
}
/// <summary>
/// This creates a new shapefile that has Z values and follows along the same line segments.
/// The boundaries for grid cells are marked with vertices and the segment is given a Z value
/// that corresponds to the grid elevation it intersects.
/// </summary>
/// <param name="mwElevGrid">A MapWinGIS Grid that contains the elevations.</param>
/// <param name="mwPolyLine">A MapWinGIS Shapefile that shows the pathways of the cross sections in the X-Y direction.</param>
/// <param name="OutFileName">A string containing the full path of the desired output shapefile. The extension should be *.shp</param>
/// <param name="CrossSectionType">Clarifies the type of output. default = PolyLineWithZ</param>
/// <param name="ICallBack">A MapWinGIS.ICallback for progress messages. [Optional]</param>
/// <remarks>This function throws Argument or Application exceptions on errors, so it's recommended that coders enclose it in a try catch block.</remarks>
public static void GetCrossSection(MapWinGIS.Grid mwElevGrid, MapWinGIS.Shapefile mwPolyLine, string OutFileName, CrossSectionTypes CrossSectionType, MapWinGIS.ICallback ICallBack)
{
MapWinUtility.Logger.Dbg("GetCrossSection(mwElevGrid: " + Macro.ParamName(mwElevGrid) + ",\n" +
" mwPolyLine: " + Macro.ParamName(mwPolyLine) + ",\n" +
" OutFileName: " + OutFileName + ",\n" +
" CrossSectionType: " + CrossSectionType.ToString() + ",\n" +
" ICallback)");
bool res;
int Prog = 0;
int OldProg = 0;
double dS, dX, dY, XllCenter, YllCenter;
int NumRows, NumCols, ElevField, IDField;
ElevField = 1;
IDField = 0;
// Test to be sure that the elevation grid and shapefile are not null
if (mwElevGrid == null)
{
MapWinUtility.Logger.Dbg("Argument Exception: Elevation grid mwElevGrid can't be null.");
throw new ArgumentException("Elevation grid mwElevGrid can't be null.");
}
if (mwPolyLine == null)
{
MapWinUtility.Logger.Dbg("Argument Exception: The shapefile of input cross sections mwPolyLine can't be null.");
throw new ArgumentException("The shapefile of input cross sections mwPolyLine can't be null.");
}
// Clear any existing shapefile output filenames that might cause problems if they exist.
string fn;
if (System.IO.Directory.Exists(System.IO.Path.GetDirectoryName(OutFileName)) == false)
{
System.IO.Directory.CreateDirectory(System.IO.Path.GetDirectoryName(OutFileName));
}
if (System.IO.File.Exists(OutFileName))
{
System.IO.File.Delete(OutFileName);
}
fn = System.IO.Path.ChangeExtension(OutFileName, ".dbf");
if (System.IO.File.Exists(fn))
{
System.IO.File.Delete(fn);
}
fn = System.IO.Path.ChangeExtension(OutFileName, ".shx");
if (System.IO.File.Exists(fn))
{
System.IO.File.Delete(fn);
}
fn = System.IO.Path.ChangeExtension(OutFileName, ".prj");
if (System.IO.File.Exists(fn))
{
System.IO.File.Delete(fn);
}
// Since we are given the lower left coordinate, just make sure dX and dY are positive
dX = Math.Abs(mwElevGrid.Header.dX);
if (dX == 0)
{
throw new ApplicationException("mwElevGrid.Header.dX cannot be 0.");
}
dY = Math.Abs(mwElevGrid.Header.dY);
if (dY == 0)
{
throw new ApplicationException("mwElevGrid.Header.dY cannot be 0.");
}
// Determine the stepping distance from the grid coordintes
dS = dX / 2;
if (dY < dX)
{
dS = dY / 2;
}
XllCenter = mwElevGrid.Header.XllCenter;
YllCenter = mwElevGrid.Header.YllCenter;
NumRows = mwElevGrid.Header.NumberRows;
NumCols = mwElevGrid.Header.NumberCols;
// Test for intersection between the entire shapefile and the grid
double left, right, top, bottom;
left = XllCenter - dX / 2;
right = XllCenter + NumCols * dX - dX / 2;
bottom = YllCenter - dY / 2;
top = YllCenter + NumRows * dY - dY / 2;
MapWinGeoProc.Topology2D.Envelope gExt = new MapWinGeoProc.Topology2D.Envelope(left, right, bottom, top);
MapWinGeoProc.Topology2D.Envelope pExt = new MapWinGeoProc.Topology2D.Envelope(mwPolyLine.Extents);
if (gExt.Intersects(pExt) == false)
{
MapWinUtility.Logger.Dbg("Application Exception: The shapefile doesn't overlap the grid, so no cross sections were found.");
throw new ApplicationException("The shapefile doesn't overlap the grid, so no cross sections were found.");
}
// Setup the output shapefile and the basic shape objects
MapWinGIS.Shape mwShape;
MapWinGIS.Shapefile sfOut = new MapWinGIS.Shapefile();
sfOut.Projection = mwPolyLine.Projection;
if (CrossSectionType == CrossSectionTypes.PointsWithZAndElevField)
{
res = sfOut.CreateNew(OutFileName, MapWinGIS.ShpfileType.SHP_POINTZ);
if (res != true)
{
MapWinUtility.Logger.Dbg("Application Exception: " + sfOut.get_ErrorMsg(sfOut.LastErrorCode));
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
res = sfOut.StartEditingShapes(true, ICallBack);
if (res != true)
{
MapWinUtility.Logger.Dbg("Application Exception: " + sfOut.get_ErrorMsg(sfOut.LastErrorCode));
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
MapWinGIS.Field ID = new MapWinGIS.Field();
ID.Name = "ID";
ID.Type = MapWinGIS.FieldType.INTEGER_FIELD;
res = sfOut.EditInsertField(ID, ref IDField, ICallBack);
if (res != true)
{
MapWinUtility.Logger.Dbg("Application Exception: " + sfOut.get_ErrorMsg(sfOut.LastErrorCode));
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
MapWinGIS.Field Elev = new MapWinGIS.Field();
Elev.Name = "Elevation";
Elev.Type = MapWinGIS.FieldType.DOUBLE_FIELD;
res = sfOut.EditInsertField(Elev, ref ElevField, ICallBack);
if (res != true)
{
MapWinUtility.Logger.Dbg("Application Exception: " + sfOut.get_ErrorMsg(sfOut.LastErrorCode));
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
}
else
{
res = sfOut.CreateNew(OutFileName, MapWinGIS.ShpfileType.SHP_POLYLINEZ);
if (res != true)
{
MapWinUtility.Logger.Dbg("Application Exception: " + sfOut.get_ErrorMsg(sfOut.LastErrorCode));
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
res = sfOut.StartEditingShapes(true, ICallBack);
if (res != true)
{
MapWinUtility.Logger.Dbg("Application Exception: " + sfOut.get_ErrorMsg(sfOut.LastErrorCode));
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
MapWinGIS.Field ID = new MapWinGIS.Field();
ID.Name = "ID";
ID.Type = MapWinGIS.FieldType.INTEGER_FIELD;
res = sfOut.EditInsertField(ID, ref IDField, ICallBack);
if (res != true)
{
MapWinUtility.Logger.Dbg("Application Exception: " + sfOut.get_ErrorMsg(sfOut.LastErrorCode));
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
}
MapWinGIS.Shape mwOutShape;
// Loop through all the shapes in the polyline shapefile
int shIndx = -1;
for (int shp = 0; shp < mwPolyLine.NumShapes; shp++)
{
mwShape = mwPolyLine.get_Shape(shp);
// By turning multi-part polylines into multiple linestrings, we avoid the annoying multi-part logic
MultiLineString MLS = GeometryFactory.CreateMultiLineString(mwShape);
for (int ls = 0; ls < MLS.NumGeometries; ls++)
{
LineString lsSource = MLS.GeometryN[ls] as LineString;
for (int pt = 0; pt < lsSource.Coordinates.Count - 1; pt++)
{
Coordinate start = lsSource.Coordinates[pt];
Coordinate end = lsSource.Coordinates[pt + 1];
// Crop the parts of each segment that do not overlap with the grid.
if (start.X < left)
{
if (end.X < left)
{
// this segment is outside the grid
continue;
}
// crop this segment to only the portion on the grid
start.X = left;
}
if (end.X < left)
{
// crop this segment to only the portion on the grid
end.X = left;
}
if (start.X > right)
{
if (end.X > right)
{
// this segment is outside the grid
continue;
}
// crop to grid
start.X = right;
}
if (end.X > right)
{
// crop to the grid
end.X = right;
}
double length = Math.Sqrt((end.X - start.X) * (end.X - start.X) + (end.Y - start.Y) * (end.Y - start.Y));
int NumSteps = (int)Math.Floor(length / dS);
double segDx = (end.X - start.X) / NumSteps;
double segDy = (end.Y - start.Y) / NumSteps;
mwOutShape = new MapWinGIS.Shape();
if (CrossSectionType == CrossSectionTypes.PolyLineWithZ)
{
mwOutShape.Create(MapWinGIS.ShpfileType.SHP_POLYLINEZ);
}
// step by dS and get the grid value at that point at each step
int p = 0;
for (int I = 0; I < NumSteps; I++)
{
int row, col;
object Elev;
double X = start.X + segDx * I;
double Y = start.Y + segDy * I;
mwElevGrid.ProjToCell(X, Y, out col, out row);
Elev = mwElevGrid.get_Value(col, row);
MapWinGIS.Point pnt = new MapWinGIS.Point();
pnt.x = X;
pnt.y = Y;
pnt.Z = (double)Elev;
if (CrossSectionType == CrossSectionTypes.PointsWithZAndElevField)
{
p = 0;
mwOutShape = new MapWinGIS.Shape();
mwOutShape.Create(MapWinGIS.ShpfileType.SHP_POINTZ);
res = mwOutShape.InsertPoint(pnt, ref p);
if (res == false)
{
throw new ApplicationException(mwOutShape.get_ErrorMsg(mwOutShape.LastErrorCode));
}
res = sfOut.EditInsertShape(mwOutShape, ref shIndx);
if (res != true)
{
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
res = sfOut.EditCellValue(IDField, shIndx, shIndx);
if (res != true)
{
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
res = sfOut.EditCellValue(ElevField, shIndx, Elev);
if (res != true)
{
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
shIndx++;
}
else
{
res = mwOutShape.InsertPoint(pnt, ref p);
p++;
if (res == false)
{
throw new ApplicationException(mwOutShape.get_ErrorMsg(mwOutShape.LastErrorCode));
}
}
}
if (CrossSectionType == CrossSectionTypes.PolyLineWithZ)
{
if (mwOutShape.numPoints > 0)
{
res = sfOut.EditInsertShape(mwOutShape, ref shIndx);
if (res != true)
{
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
res = sfOut.EditCellValue(IDField, shIndx, shIndx);
if (res != true)
{
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
shIndx++;
}
}
}
}
Prog = Convert.ToInt32(shp * 100 / mwPolyLine.NumShapes);
if (Prog > OldProg)
{
MapWinUtility.Logger.Progress("Evaluating Cross Section..." + Prog.ToString() + "% Complete.", Prog, OldProg);
OldProg = Prog;
}
}
res = sfOut.StopEditingShapes(true, true, ICallBack);
if (res != true)
{
MapWinUtility.Logger.Dbg("Application Exception: " + sfOut.get_ErrorMsg(sfOut.LastErrorCode));
throw new ApplicationException(sfOut.get_ErrorMsg(sfOut.LastErrorCode));
}
}