public override bool Draw(RenderContext11 renderContext, float opacity, bool flat)
{
if (shapefile == null)
{
return(false);
}
if (shapeFileVertex == null)
{
List <Vector3> vertList = new List <Vector3>();
List <UInt32> indexList = new List <UInt32>();
UInt32 firstItemIndex = 0;
Vector3 lastItem = new Vector3();
bool firstItem = true;
bool north = true;
double offsetX = 0;
double offsetY = 0;
double centralMeridian = 0;
double mapScale = 0;
double standardParallel = 70;
if (shapefile.Projection == ShapeFile.Projections.PolarStereo)
{
north = shapefile.FileHeader.ProjectionInfo.Name.ToLower().Contains("north");
standardParallel = shapefile.FileHeader.ProjectionInfo.GetParameter("standard_parallel_1");
centralMeridian = shapefile.FileHeader.ProjectionInfo.GetParameter("central_meridian");
mapScale = shapefile.FileHeader.ProjectionInfo.GetParameter("scale_factor");
offsetY = shapefile.FileHeader.ProjectionInfo.GetParameter("false_easting");
offsetX = shapefile.FileHeader.ProjectionInfo.GetParameter("false_northing");
}
UInt32 currentIndex = 0;
Color color = Color;
int count = 360;
for (int i = 0; i < shapefile.Shapes.Count; i++)
{
if (shapefile.Shapes[i].GetType() == typeof(Polygon))
{
Polygon p = (Polygon)shapefile.Shapes[i];
for (int z = 0; z < p.Rings.Length; z++)
{
count = (p.Rings[z].Points.Length);
// content from DBF
DataRow dr = p.Rings[z].Attributes;
for (int k = 0; k < p.Rings[z].Points.Length; k++)
{
// 2D Point coordinates. 3d also supported which would add a Z. There's also an optional measure (M) that can be used.
double Xcoord = p.Rings[z].Points[k].X;
double Ycoord = p.Rings[z].Points[k].Y;
if (shapefile.Projection == ShapeFile.Projections.Geo)
{
lastItem = Coordinates.GeoTo3d(Ycoord, Xcoord);
}
else if (shapefile.Projection == ShapeFile.Projections.PolarStereo)
{
lastItem = Coordinates.SterographicTo3d(Xcoord, Ycoord, 1, standardParallel, centralMeridian, offsetX, offsetY, mapScale, north).Vector3;
}
if (k == 0)
{
firstItemIndex = currentIndex;
firstItem = true;
}
vertList.Add(lastItem);
if (firstItem)
{
firstItem = false;
}
else
{
indexList.Add(currentIndex);
currentIndex++;
indexList.Add(currentIndex);
}
}
indexList.Add(currentIndex);
indexList.Add(firstItemIndex);
currentIndex++;
}
}
else if (shapefile.Shapes[i].GetType() == typeof(PolygonZ))
{
PolygonZ p = (PolygonZ)shapefile.Shapes[i];
for (int z = 0; z < p.Rings.Length; z++)
{
count = (p.Rings[z].Points.Length);
// content from DBF
DataRow dr = p.Rings[z].Attributes;
for (int k = 0; k < p.Rings[z].Points.Length; k++)
{
// 2D Point coordinates. 3d also supported which would add a Z. There's also an optional measure (M) that can be used.
double Xcoord = p.Rings[z].Points[k].X;
double Ycoord = p.Rings[z].Points[k].Y;
if (shapefile.Projection == ShapeFile.Projections.Geo)
{
lastItem = Coordinates.GeoTo3d(Ycoord, Xcoord);
}
else if (shapefile.Projection == ShapeFile.Projections.PolarStereo)
{
lastItem = Coordinates.SterographicTo3d(Xcoord, Ycoord, 1, standardParallel, centralMeridian, offsetX, offsetY, mapScale, north).Vector3;
}
if (k == 0)
{
firstItemIndex = currentIndex;
firstItem = true;
}
vertList.Add(lastItem);
if (firstItem)
{
firstItem = false;
}
else
{
indexList.Add(currentIndex);
currentIndex++;
indexList.Add(currentIndex);
}
}
indexList.Add(currentIndex);
indexList.Add(firstItemIndex);
currentIndex++;
}
}
else if (shapefile.Shapes[i].GetType() == typeof(PolyLine))
{
PolyLine p = (PolyLine)shapefile.Shapes[i];
for (int z = 0; z < p.Lines.Length; z++)
{
count = (p.Lines[z].Points.Length);
firstItem = true;
for (int k = 0; k < p.Lines[z].Points.Length; k++)
{
// 2D Point coordinates. 3d also supported which would add a Z. There's also an optional measure (M) that can be used.
double Xcoord = p.Lines[z].Points[k].X;
double Ycoord = p.Lines[z].Points[k].Y;
if (shapefile.Projection == ShapeFile.Projections.Geo)
{
lastItem = Coordinates.GeoTo3d(Ycoord, Xcoord);
}
else if (shapefile.Projection == ShapeFile.Projections.PolarStereo)
{
lastItem = Coordinates.SterographicTo3d(Xcoord, Ycoord, 1, standardParallel, centralMeridian, offsetX, offsetY, mapScale, north).Vector3;
}
if (k == 0)
{
firstItemIndex = currentIndex;
firstItem = true;
}
vertList.Add(lastItem);
if (firstItem)
{
firstItem = false;
}
else
{
indexList.Add(currentIndex);
currentIndex++;
indexList.Add(currentIndex);
}
}
currentIndex++;
}
}
else if (shapefile.Shapes[i].GetType() == typeof(PolyLineZ))
{
PolyLineZ p = (PolyLineZ)shapefile.Shapes[i];
for (int z = 0; z < p.Lines.Length; z++)
{
count = (p.Lines[z].Points.Length);
Vector3[] points = new Vector3[(count)];
firstItem = true;
for (int k = 0; k < p.Lines[z].Points.Length; k++)
{
// 2D Point coordinates. 3d also supported which would add a Z. There's also an optional measure (M) that can be used.
double Xcoord = p.Lines[z].Points[k].X;
double Ycoord = p.Lines[z].Points[k].Y;
if (shapefile.Projection == ShapeFile.Projections.Geo)
{
lastItem = Coordinates.GeoTo3d(Ycoord, Xcoord);
}
else if (shapefile.Projection == ShapeFile.Projections.PolarStereo)
{
lastItem = Coordinates.SterographicTo3d(Xcoord, Ycoord, 1, standardParallel, centralMeridian, offsetX, offsetY, mapScale, north).Vector3;
}
if (k == 0)
{
firstItemIndex = currentIndex;
firstItem = true;
}
vertList.Add(lastItem);
if (firstItem)
{
firstItem = false;
}
else
{
indexList.Add(currentIndex);
currentIndex++;
indexList.Add(currentIndex);
}
}
currentIndex++;
}
}
else if (shapefile.Shapes[i].GetType() == typeof(ShapefileTools.Point))
{
ShapefileTools.Point p = (ShapefileTools.Point)shapefile.Shapes[i];
// 2D Point coordinates. 3d also supported which would add a Z. There's also an optional measure (M) that can be used.
double Xcoord = p.X;
double Ycoord = p.Y;
if (shapefile.Projection == ShapeFile.Projections.Geo)
{
lastItem = Coordinates.GeoTo3d(Ycoord, Xcoord);
}
else if (shapefile.Projection == ShapeFile.Projections.PolarStereo)
{
lastItem = Coordinates.SterographicTo3d(Xcoord, Ycoord, 1, standardParallel, centralMeridian, offsetX, offsetY, mapScale, north).Vector3;
}
vertList.Add(lastItem);
currentIndex++;
lines = false;
}
}
shapeVertexCount = vertList.Count;
shapeFileVertex = new PositionVertexBuffer11(vertList.Count, RenderContext11.PrepDevice);
Vector3[] verts = (Vector3[])shapeFileVertex.Lock(0, 0); // Lock the buffer (which will return our structs)
int indexer = 0;
foreach (Vector3 vert in vertList)
{
verts[indexer++] = vert;
}
shapeFileVertex.Unlock();
shapeIndexCount = indexList.Count;
if (lines)
{
if (indexList.Count > 65500)
{
isLongIndex = true;
shapeFileIndex = new IndexBuffer11(typeof(UInt32), indexList.Count, RenderContext11.PrepDevice);
}
else
{
isLongIndex = false;
shapeFileIndex = new IndexBuffer11(typeof(short), indexList.Count, RenderContext11.PrepDevice);
}
if (isLongIndex)
{
indexer = 0;
UInt32[] indexes = (UInt32[])shapeFileIndex.Lock();
foreach (UInt32 indexVal in indexList)
{
indexes[indexer++] = indexVal;
}
shapeFileIndex.Unlock();
}
else
{
indexer = 0;
short[] indexes = (short[])shapeFileIndex.Lock();
foreach (UInt32 indexVal in indexList)
{
indexes[indexer++] = (short)indexVal;
}
shapeFileIndex.Unlock();
}
}
}
renderContext.DepthStencilMode = DepthStencilMode.Off;
renderContext.BlendMode = BlendMode.Alpha;
SimpleLineShader11.Color = Color.FromArgb((int)(opacity * 255), Color);
SharpDX.Matrix mat = (renderContext.World * renderContext.View * renderContext.Projection).Matrix11;
mat.Transpose();
SimpleLineShader11.WVPMatrix = mat;
SimpleLineShader11.CameraPosition = Vector3d.TransformCoordinate(renderContext.CameraPosition, Matrix3d.Invert(renderContext.World)).Vector3;
SimpleLineShader11.ShowFarSide = false;
SimpleLineShader11.Sky = false;
renderContext.SetVertexBuffer(shapeFileVertex);
SimpleLineShader11.Use(renderContext.devContext);
if (lines)
{
renderContext.devContext.InputAssembler.PrimitiveTopology = SharpDX.Direct3D.PrimitiveTopology.LineList;
renderContext.SetIndexBuffer(shapeFileIndex);
renderContext.devContext.DrawIndexed(shapeFileIndex.Count, 0, 0);
}
else
{
renderContext.devContext.InputAssembler.PrimitiveTopology = SharpDX.Direct3D.PrimitiveTopology.PointList;
renderContext.devContext.Draw(shapeVertexCount, 0);
}
renderContext.devContext.InputAssembler.PrimitiveTopology = SharpDX.Direct3D.PrimitiveTopology.TriangleList;
return(true);
}
public string ToWellKnownText()
{
if (shapefile == null)
{
return(null);
}
StringBuilder sb = new StringBuilder();
// Make Header
sb.Append("Geography");
if (shapefile.Table != null)
{
foreach (DataColumn col in shapefile.Table.Columns)
{
sb.Append("\t");
sb.Append(col.ColumnName);
}
}
sb.AppendLine();
bool firstItem = true;
bool firstShape = true;
bool north = true;
double offsetX = 0;
double offsetY = 0;
double centralMeridian = 0;
double mapScale = 0;
double standardParallel = 70;
if (shapefile.Projection == ShapeFile.Projections.PolarStereo)
{
north = shapefile.FileHeader.ProjectionInfo.Name.ToLower().Contains("north");
standardParallel = shapefile.FileHeader.ProjectionInfo.GetParameter("standard_parallel_1");
centralMeridian = shapefile.FileHeader.ProjectionInfo.GetParameter("central_meridian");
mapScale = shapefile.FileHeader.ProjectionInfo.GetParameter("scale_factor");
offsetY = shapefile.FileHeader.ProjectionInfo.GetParameter("false_easting");
offsetX = shapefile.FileHeader.ProjectionInfo.GetParameter("false_northing");
}
int color = Color.ToArgb();
int count = 360;
for (int i = 0; i < shapefile.Shapes.Count; i++)
{
if (shapefile.Shapes[i].GetType() == typeof(Polygon))
{
firstShape = true;
firstItem = true;
sb.Append("Polygon (");
Polygon p = (Polygon)shapefile.Shapes[i];
for (int z = 0; z < p.Rings.Length; z++)
{
if (firstShape)
{
firstShape = false;
}
else
{
sb.Append(",");
}
sb.Append("(");
count = (p.Rings[z].Points.Length);
// content from DBF
DataRow dr = p.Rings[z].Attributes;
firstItem = true;
for (int k = 0; k < p.Rings[z].Points.Length; k++)
{
// 2D Point coordinates. 3d also supported which would add a Z. There's also an optional measure (M) that can be used.
double Xcoord = p.Rings[z].Points[k].X;
double Ycoord = p.Rings[z].Points[k].Y;
if (firstItem)
{
firstItem = false;
}
else
{
sb.Append(",");
}
// sb.Append(Xcoord.ToString("F5") + " " + Ycoord.ToString("F5"));
sb.Append(Xcoord.ToString() + " " + Ycoord.ToString());
}
sb.Append(")");
}
sb.Append(")");
AddTableRow(sb, shapefile.Shapes[i].Attributes);
sb.AppendLine();
}
else if (shapefile.Shapes[i].GetType() == typeof(PolyLine))
{
PolyLine p = (PolyLine)shapefile.Shapes[i];
for (int z = 0; z < p.Lines.Length; z++)
{
firstShape = true;
firstItem = true;
sb.Append("LineString (");
count = (p.Lines[z].Points.Length);
PositionColoredTextured[] points = new PositionColoredTextured[(count)];
firstItem = true;
for (int k = 0; k < p.Lines[z].Points.Length; k++)
{
//if (firstShape)
//{
// firstShape = false;
//}
//else
//{
// sb.Append(",");
//}
// 2D Point coordinates. 3d also supported which would add a Z. There's also an optional measure (M) that can be used.
double Xcoord = p.Lines[z].Points[k].X;
double Ycoord = p.Lines[z].Points[k].Y;
if (firstItem)
{
firstItem = false;
}
else
{
sb.Append(",");
}
sb.Append(Xcoord.ToString() + " " + Ycoord.ToString());
}
sb.Append(")");
AddTableRow(sb, shapefile.Shapes[i].Attributes);
sb.AppendLine();
}
}
else if (shapefile.Shapes[i].GetType() == typeof(PolyLineZ))
{
PolyLineZ p = (PolyLineZ)shapefile.Shapes[i];
for (int z = 0; z < p.Lines.Length; z++)
{
firstShape = true;
firstItem = true;
sb.Append("LineString (");
count = (p.Lines[z].Points.Length);
PositionColoredTextured[] points = new PositionColoredTextured[(count)];
firstItem = true;
for (int k = 0; k < p.Lines[z].Points.Length; k++)
{
//if (firstShape)
//{
// firstShape = false;
//}
//else
//{
// sb.Append(",");
//}
// 2D Point coordinates. 3d also supported which would add a Z. There's also an optional measure (M) that can be used.
double Xcoord = p.Lines[z].Points[k].X;
double Ycoord = p.Lines[z].Points[k].Y;
double Zcoord = p.Lines[z].Points[k].Z;
if (firstItem)
{
firstItem = false;
}
else
{
sb.Append(",");
}
sb.Append(Xcoord.ToString() + " " + Ycoord.ToString() + " " + Zcoord.ToString());
}
sb.Append(")");
AddTableRow(sb, shapefile.Shapes[i].Attributes);
sb.AppendLine();
}
}
else if (shapefile.Shapes[i].GetType() == typeof(ShapefileTools.Point))
{
ShapefileTools.Point p = (ShapefileTools.Point)shapefile.Shapes[i];
// 2D Point coordinates. 3d also supported which would add a Z. There's also an optional measure (M) that can be used.
double Xcoord = p.X;
double Ycoord = p.Y;
sb.Append(string.Format("Point ({0} {1})"));
AddTableRow(sb, shapefile.Shapes[i].Attributes);
sb.AppendLine();
}
}
return(sb.ToString());
}
public override IPlace FindClosest(Coordinates target, float distance, IPlace closestPlace, bool astronomical)
{
bool pointFound = false;
int pointIndex = -1;
Vector3d target3d = Coordinates.GeoTo3dDouble(target.Lat, target.Lng, 1);
for (int i = 0; i < shapefile.Shapes.Count; i++)
{
if (shapefile.Shapes[i] is ComplexShape)
{
ComplexShape p = (ComplexShape)shapefile.Shapes[i];
if (PointInboundingBox(target, p.BoundingBox))
{
if (p.Attributes != null && p.Attributes.ItemArray.GetLength(0) > 0)
{
int nameCol = GetNameColumn();
if (nameCol == -1)
{
nameCol = 0;
}
TourPlace place = new TourPlace(p.Attributes.ItemArray[nameCol].ToString(), (p.BoundingBox[1] + p.BoundingBox[3]) / 2, (p.BoundingBox[0] + p.BoundingBox[2]) / 2, Classification.Unidentified, "", ImageSetType.Earth, -1);
Dictionary <String, String> rowData = new Dictionary <string, string>();
for (int r = 0; r < p.Attributes.ItemArray.GetLength(0); r++)
{
rowData.Add(p.Attributes.Table.Columns[r].ColumnName, p.Attributes.ItemArray[r].ToString());
}
place.Tag = rowData;
return(place);
}
return(closestPlace);
}
}
else if (shapefile.Shapes[i].GetType() == typeof(ShapefileTools.Point))
{
ShapefileTools.Point p = (ShapefileTools.Point)shapefile.Shapes[i];
Vector3d point = Coordinates.GeoTo3dDouble(p.Y, p.X, 1);
Vector3d dist = target3d - point;
if (dist.Length() < distance)
{
pointFound = true;
pointIndex = i;
distance = (float)dist.Length();
}
}
}
if (pointFound)
{
ShapefileTools.Point p = (ShapefileTools.Point)shapefile.Shapes[pointIndex];
if (p.Attributes.ItemArray.GetLength(0) > 0)
{
TourPlace place = new TourPlace(p.Attributes.ItemArray[0].ToString(), p.Y, p.X, Classification.Unidentified, "", ImageSetType.Earth, -1);
Dictionary <String, String> rowData = new Dictionary <string, string>();
for (int r = 0; r < p.Attributes.ItemArray.GetLength(0); r++)
{
rowData.Add(p.Attributes.Table.Columns[r].ColumnName, p.Attributes.ItemArray[r].ToString());
}
place.Tag = rowData;
return(place);
}
}
return(closestPlace);
}
