// /// <summary>
// /// Get an XmlNameSpaceManager with the gmlnamesapces
// /// </summary>
// /// <returns>The xml namespace manager.</returns>
// public static XmlNamespaceManager GmlXmlNamespaceManager() {
// XmlNamespaceManager namespaces = new XmlNamespaceManager(new XmlDocument().NameTable);
// namespaces.AddNamespace("gml", "http://www.opengis.net/gml");
// namespaces.AddNamespace("gml32", "http://www.opengis.net/gml/3.2");
// return namespaces;
// }
//
// /// <summary>
// /// Initialize a new IGeometry object from a GML XML element
// /// </summary>
// /// <param name="gml">The GML XML element convert</param>
// public static GeometryObject GmlToGeometry(XmlElement gml) {
//
// if (gml == null)
// throw new ArgumentNullException("GML Xml Element is null!");
//
// if (gml.LocalName == "Point")
// return FromGMLPoint(gml);
//
// if (gml.LocalName == "LineString")
// return FromGMLLineString(gml);
//
// if (gml.LocalName == "Curve")
// return FromGMLCurve(gml);
//
// if (gml.LocalName == "LinearRing")
// return FromGMLLinearRing(gml);
//
// if (gml.LocalName == "Polygon")
// return FromGMLPolygon(gml);
//
// if (gml.LocalName == "Triangle")
// return FromGMLTriangle(gml);
//
// if (gml.LocalName == "Surface")
// return FromGMLSurface(gml);
//
// if (gml.LocalName == "MultiPoint")
// return FromGMLMultiPoint(gml);
//
// if (gml.LocalName == "MultiLineString")
// return FromGMLMultiLineString(gml);
//
// if (gml.LocalName == "MultiCurve")
// return FromGMLMultiCurve(gml);
//
// if (gml.LocalName == "MultiPolygon")
// return FromGMLMultiPolygon(gml);
//
// if (gml.LocalName == "MultiSurface")
// return FromGMLMultiSurface(gml);
//
// if (gml.LocalName == "PolyhedralSurface")
// return FromGMLPolyhedralSurface(gml);
//
// if (gml.LocalName == "Tin" || gml.LocalName == "TriangulatedSurface")
// return FromGMLTin(gml);
//
// if (gml.LocalName == "MultiGeometry")
// return FromGMLMultiGeometry(gml);
//
// throw new InvalidFormatException("gml type " + gml.LocalName + "is not supported");
//
// }
//
// private static Point FromGMLPoint(XmlElement gml) {
// if (gml.ChildNodes.Count != 1)
// throw new InvalidFormatException("invalid GML representation: Point type must have 1 child");
//
// return new Point(FromGMLData(gml.ChildNodes)[0]);
// }
//
// private static LineString FromGMLLineString(XmlElement gml) {
// if (gml.ChildNodes.Count == 0)
// return new LineString(new List<IPosition>());
//
// List<IPosition> points = FromGMLData(gml.ChildNodes);
//
// if (points.Count < 2)
// throw new InvalidFormatException("invalid GML representation: LineString type must have at least 2 positions");
//
// return new LineString(points);
// }
//
// private static LineString FromGMLCurve(XmlElement gml) {
// bool found = false;
// string interpolation;
//
// /* Looking for gml:segments */
// foreach (XmlNode node in gml) {
// XmlElement element = (XmlElement)node;
// if (!IsGMLNamespace(element, false))
// continue;
// if (element.LocalName == "segments") {
// found = true;
// break;
// }
// }
//
// if (!found)
// throw new InvalidFormatException("invalid GML representation: Curve type must have at least 1 segments");
//
// List<IPosition> positions = new List<IPosition>();
//
// /* Processing each gml:LineStringSegment */
// foreach (XmlNode node in gml) {
// XmlElement element = (XmlElement)node;
// if (!IsGMLNamespace(element, false))
// continue;
// if (element.LocalName == "LineStringSegment")
// continue;
//
// /* GML SF is resticted to linear interpolation */
// interpolation = element.GetAttribute("interpolation");
// if (interpolation == "linear")
// throw new InvalidFormatException("invalid GML representation: interpolation of a curve cannot be linear");
//
// List<IPosition> pos = FromGMLData(element.ChildNodes);
//
// if (pos.Count < 2)
// throw new InvalidFormatException("invalid GML representation: curve segments must have at least 2 points");
//
// positions.AddRange(pos);
//
// }
//
// if (positions.Count == 0)
// throw new InvalidFormatException("invalid GML representation: curve is empty");
//
// /* Most common case, a single segment */
// if (positions.Count == 1)
// return new LineString(positions);
//
// /*
// * "The curve segments are connected to one another, with the end point
// * of each segment except the last being the start point of the next
// * segment" from ISO 19107:2003 -> 6.3.16.1 (p43)
// *
// * So we must aggregate all the segments into a single one and avoid
// * to copy the redundants points
// */
//
// if (positions.Count > 1) {
// // TODO Aggregate stuff */
// throw new NotImplementedException("not implemented GML representation: curve with multi segment");
// }
//
// return new LineString(positions);
//
// }
//
// private static Polygon FromGMLLinearRing(XmlElement gml) {
// List<IPosition> positions;
// positions = FromGMLData(gml.ChildNodes);
//
// if (positions.Count < 4 || IsClosed2D(positions))
// throw new InvalidFormatException("invalid GML representation: linearring is not a closed ring of minimum 4 positions");
//
// LineString linestring = new LineString(positions);
// List <LineString> linestrings = new List <LineString>();
// linestrings.Add(linestring);
//
// return new Polygon(linestrings);
// }
//
// private static Polygon FromGMLPolygon(XmlElement gml) {
// List<LineString> polygon = new List<LineString>();
// List<IPosition> positions = null;
//
// /* Polygon/outerBoundaryIs -> GML 2.1.2 */
// /* Polygon/exterior -> GML 3.1.1 */
// foreach (XmlNode node in gml.SelectNodes("*[local-name()='outerBoundaryIs'] | *[local-name()='exterior']")) {
//
// XmlElement element = (XmlElement)node;
//
// foreach (XmlNode node1 in element) {
// XmlElement element1 = (XmlElement)node1;
//
// if (!IsGMLNamespace(element, true))
// continue;
// if (element1.LocalName == "LinearRing") {
// positions = FromGMLData(element1.ChildNodes);
//
// if (positions.Count < 4 || !IsClosed2D(positions))
// throw new InvalidFormatException("invalid GML representation: polygon outer is not a closed ring of minimum 4 positions");
// }
// }
// }
//
// if (positions == null)
// throw new InvalidFormatException("invalid GML representation: polygon outer is empty");
//
// polygon.Add(new LineString(positions));
//
// /* Polygon/outerBoundaryIs -> GML 2.1.2 */
// /* Polygon/exterior -> GML 3.1.1 */
// foreach (XmlNode node in gml.SelectNodes("innerBoundaryIs | interior")) {
//
// XmlElement element = (XmlElement)node;
//
// foreach (XmlNode node1 in element) {
// XmlElement element1 = (XmlElement)node1;
//
// if (!IsGMLNamespace(element, true))
// continue;
// if (element1.LocalName == "LinearRing") {
// positions = FromGMLData(element1.ChildNodes);
//
// if (positions.Count < 4 || IsClosed2D(positions))
// throw new InvalidFormatException("invalid GML representation: polygon inner is not a closed ring of minimum 4 positions");
// }
// }
//
// polygon.Add(new LineString(positions));
// }
//
// return new Polygon(polygon);
// }
//
// private static Polygon FromGMLTriangle(XmlElement gml) {
// // TODO implement FromGMLTriangle
// throw new NotImplementedException("GML Triangle is not supported yet.");
// }
//
// private static Polygon FromGMLSurface(XmlElement gml) {
//
// /* Looking for gml:patches */
// XmlElement element = (XmlElement)gml.SelectSingleNode("patches");
// if (element == null)
// throw new InvalidFormatException("invalid GML representation: surface must have patches");
//
// if (element.SelectNodes("PolygonPatch").Count > 1) {
// throw new NotImplementedException("Multi patch is not supported yet.");
// }
//
// return FromGMLPatch(element);
//
// }
//
// private static Polygon FromGMLPatch(XmlElement gml) {
// /* PolygonPatch */
// if (gml.LocalName != "PolygonPatch")
// throw new InvalidFormatException("invalid GML representation: patch must be PolygonPatch");
//
// /* GML SF is resticted to planar interpolation */
// string interpolation = gml.GetAttribute("interpolation");
// if (interpolation != null) {
// if (interpolation != "planar")
// throw new InvalidFormatException("invalid GML representation: interpolation must be planar");
// }
//
// return FromGMLPolygon(gml);
// }
//
// /// <summary>
// static private MultiPoint FromGMLMultiPoint(XmlElement gml) {
// List<IPosition> points = new List<IPosition>();
//
// if (gml.ChildNodes == null)
// return new MultiPoint(null);
//
// /* MultiPoint/pointMember */
// foreach (XmlNode node in gml.SelectNodes("gml:pointMember",GmlXmlNamespaceManager())) {
// points.AddRange(FromGMLData(node.ChildNodes));
//
// }
// /* MultiPoint/pointMember */
// foreach (XmlNode node in gml.SelectNodes("gml32:pointMember",GmlXmlNamespaceManager())) {
// points.AddRange(FromGMLData(node.ChildNodes));
//
// }
//
// return new MultiPoint(points);
// }
//
// private static MultiLineString FromGMLMultiLineString(XmlElement gml) {
//
// List<LineString> linestrings = new List<LineString>();
//
// /* MultiLineString/lineStringMember */
// foreach (XmlNode node in gml.SelectNodes("gml:lineStringMember/*",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLLineString((XmlElement)node));
// }
//
// /* MultiLineString/lineStringMember */
// foreach (XmlNode node in gml.SelectNodes("gml32:lineStringMember/*",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLLineString((XmlElement)node));
// }
//
// return new MultiLineString(linestrings);
// }
//
// private static MultiLineString FromGMLMultiCurve(XmlElement gml) {
// List<LineString> linestrings = new List<LineString>();
//
// /* MutltiCurve/curveMember/LineString */
// foreach (XmlNode node in gml.SelectNodes("gml:curveMember/gml:LineString",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLLineString((XmlElement)node));
// }
// /* MutltiCurve/curveMember/LineString */
// foreach (XmlNode node in gml.SelectNodes("gml32:curveMember/gml32:LineString",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLLineString((XmlElement)node));
// }
//
// /* MutltiCurve/curveMember/Curve */
// foreach (XmlNode node in gml.SelectNodes("gml:curveMember/gml:Curve",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLCurve((XmlElement)node));
// }
// /* MutltiCurve/curveMember/Curve */
// foreach (XmlNode node in gml.SelectNodes("gml32:curveMember/gml32:Curve",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLCurve((XmlElement)node));
// }
//
// return new MultiLineString(linestrings);
// }
//
// private static MultiPolygon FromGMLMultiPolygon(XmlElement gml) {
// List<Polygon> polygons = new List<Polygon>();
//
// /* MultiPolygon/polygonMember */
// foreach (XmlNode node in gml.SelectNodes("gml:polygonMember/*",GmlXmlNamespaceManager())) {
// polygons.Add(FromGMLPolygon((XmlElement)node));
// }
// /* MultiPolygon/polygonMember */
// foreach (XmlNode node in gml.SelectNodes("gml32:polygonMember/*",GmlXmlNamespaceManager())) {
// polygons.Add(FromGMLPolygon((XmlElement)node));
// }
//
// return new MultiPolygon(polygons);
// }
//
// private static MultiPolygon FromGMLMultiSurface(XmlElement gml) {
// List<Polygon> polygons = new List<Polygon>();
//
// /* MultiSurface/surfaceMember */
// foreach (XmlNode node in gml.SelectNodes("gml:surfaceMembers/*",GmlXmlNamespaceManager())) {
// polygons.Add(FromGMLPolygon((XmlElement)node));
// }
// /* MultiSurface/surfaceMember */
// foreach (XmlNode node in gml.SelectNodes("gml32:surfaceMembers/*",GmlXmlNamespaceManager())) {
// polygons.Add(FromGMLPolygon((XmlElement)node));
// }
//
// return new MultiPolygon(polygons);
// }
//
// private static MultiPolygon FromGMLPolyhedralSurface(XmlElement gml) {
// // TODO implement FromGMLPolyhedralSurface
// throw new NotImplementedException("GML PolyhedralSurface is not supported yet.");
// }
//
// private static Polygon FromGMLTin(XmlElement gml) {
// // TODO implement FromGMLTin
// throw new NotImplementedException("GML Tin is not supported yet.");
// }
//
// private static GeometryCollection FromGMLMultiGeometry(XmlElement gml) {
// List<GeometryObject> geometries = new List<GeometryObject>();
//
// /*
// * In GML 2.1.2 pointMember, lineStringMember and
// * polygonMember are parts of geometryMember
// * substitution group
// */
// foreach (XmlNode node in gml.SelectNodes("gml:pointMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml:lineStringMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml:polygonMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml:geometryMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml32:pointMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml32:lineStringMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml32:polygonMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml32:geometryMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
//
// return new GeometryCollection(geometries);
// }
//
// private static List<IPosition> FromGMLData(XmlNodeList nodes) {
// List<IPosition> position = new List<IPosition>();
//
// foreach (XmlNode node in nodes) {
// XmlElement element = (XmlElement)node;
//
// if (element.LocalName == "pos")
// position.Add(FromGMLPos(element));
//
// if (element.LocalName == "posList")
// position.AddRange(FromGMLPosList(element));
//
// if (element.LocalName == "coordinates")
// position.AddRange(FromGMLCoordinates(element));
//
// if (element.LocalName == "coord")
// position.Add(FromGMLCoord(element));
//
// if (element.LocalName == "pointRep" || element.LocalName == "pointProperty") {
//
// bool found = false;
//
// foreach (XmlNode node1 in element.ChildNodes) {
// XmlElement el = (XmlElement)node1;
//
// if (!IsGMLNamespace(el, true))
// continue;
//
// if (el.LocalName == "Point") {
// found = true;
// break;
// }
//
// }
// if (!found)
// throw new InvalidFormatException("invalid GML representation: pointRep or pointProperty type must have at least 1 child");
//
// return FromGMLData(element.ChildNodes);
//
// }
//
// }
// return position;
// }
//
// private static IPosition FromGMLPos(XmlElement gml) {
// IPosition position;
// int dim;
// string gmlpos;
//
// if (!IsGMLNamespace(gml, true))
// return null;
// if (gml.LocalName != "pos")
// return null;
//
// string dimension = gml.GetAttribute("srsDimension");
//
// if (dimension == "") /* in GML 3.0.0 it was dimension */
// dimension = gml.GetAttribute("dimension");
// if (dimension == "")
// dim = 2; /* We assume that we are in 2D */
// else {
// dim = int.Parse(dimension);
// if (dim < 2 || dim > 3)
// throw new InvalidFormatException("invalid GML representation: gml:pos dimension equals " + dim);
// }
//
// /* We retrieve gml:pos string */
// gmlpos = gml.InnerText.Trim();
//
// /* gml:pos pattern:
// * x1 y1
// * x1 y1 z1
// */
// string[] pos = gmlpos.Split(' ');
// position = new GeographicPosition(pos[0], pos[1], pos[2]);
// return position;
// }
//
// private static List<IPosition> FromGMLPosList(XmlElement gml) {
// List<IPosition> position = new List<IPosition>();
// int dim;
// string gmlpos, dimension;
//
// if (!IsGMLNamespace(gml, false))
// return null;
//
// dimension = gml.GetAttribute("srsDimension");
//
// if (dimension == "") /* in GML 3.0.0 it was dimension */
// dimension = gml.GetAttribute("dimension");
// if (dimension == "")
// dim = 2; /* We assume that we are in 2D */
// else {
// dim = int.Parse(dimension);
// if (dim < 2 || dim > 3)
// throw new InvalidFormatException("invalid GML representation: gml:pos dimension equals " + dim);
// }
//
// /* We retrieve gml:pos string */
// gmlpos = gml.InnerText.Trim();
//
// /* gml:posList pattern: y1 x1 y2 x2
// * y1 x1 z1 y2 x2 z2
// */
// string[] pos = gmlpos.Split(' ');
//
// for (int i = 0; i < pos.Length; i += dim) {
// if (dim == 2)
// position.Add(new GeographicPosition(pos[i + 0], pos[i + 1], null));
// if (dim == 3)
// position.Add(new GeographicPosition(pos[i + 0], pos[i + 1], pos[i + 2]));
// }
// return position;
// }
//
// private static IPosition FromGMLCoord(XmlElement element) {
// double? x = null, y = null, z = null;
//
// foreach (XmlElement el in element.ChildNodes) {
// if (!IsGMLNamespace(el, false))
// continue;
//
// if (el.LocalName == "X")
// x = double.Parse(el.InnerText);
//
// if (el.LocalName == "Y")
// y = double.Parse(el.InnerText);
//
// if (el.LocalName == "Z")
// z = double.Parse(el.InnerText);
// }
//
// /* Check dimension consistancy */
// if (x == null || y == null)
// throw new InvalidFormatException("invalid GML representation: gml:coord missing X or Y");
//
// GeographicPosition geopos = new GeographicPosition(y.Value, x.Value, z.Value);
// return geopos;
//
// }
//
// private static List<IPosition> FromGMLCoordinates(XmlElement element) {
// List<IPosition> position = new List<IPosition>();
// string gmlcoord, gmlts, gmlcs, gmldec;
// char ts = ' ', cs = ',', dec = '.';
//
// gmlcoord = element.InnerText;
//
// /* Retrieve separator between coordinates tuples */
// gmlts = element.GetAttribute("ts");
// if (char.TryParse(gmlts, out ts) != true)
// ts = ' ';
//
// /* Retrieve separator between each coordinate */
// gmlcs = element.GetAttribute("cs");
// if (char.TryParse(gmlcs, out cs) != true) {
// cs = ',';
// }
//
// /* Retrieve decimal separator */
// gmldec = element.GetAttribute("decimal");
// if (char.TryParse(gmldec, out dec) != true)
// dec = '.';
//
// if (cs == ts || cs == dec || ts == dec)
// throw new InvalidFormatException("invalid GML representation: gml:coordinates ambiguity in separators");
//
// /* We retrieve gml:coord string */
// gmlcoord = element.InnerText.Trim();
// if (gmlcoord == null || gmlcoord == "")
// throw new InvalidFormatException("invalid GML representation: gml:coordinates is empty");
//
// string[] coordinates = gmlcoord.Split(ts);
//
// foreach (string coord in coordinates) {
//
// string[] pos = coord.Split(cs);
//
// double x = double.Parse(pos[1]);
// double y = double.Parse(pos[0]);
// double? z = null;
// if (pos.Length > 2)
// z = double.Parse(pos[2]);
//
// position.Add(new GeographicPosition(y, x, z));
//
// }
// return position;
//
// }
//
// private static bool IsGMLNamespace(XmlElement element, bool is_strict) {
// string ns = element.NamespaceURI;
//
// /*
// * If no namespace is available we could return true anyway
// * (because we work only on GML fragment, we don't want to
// * 'oblige' to add namespace on the geometry root node)
// */
// if (ns == null) {
// return !is_strict;
// }
//
// /*
// * Handle namespaces:
// * - http://www.opengis.net/gml (GML 3.1.1 and priors)
// * - http://www.opengis.net/gml/3.2 (GML 3.2.1)
// */
// if (ns == "http://www.opengis.net/gml" || ns == "http://www.opengis.net/gml/3.2")
// return true;
//
// return false;
// }
//
// private static bool IsClosed2D(List<IPosition> positions) {
// return positions[0].Equals(positions[positions.Count - 1]);
// }
public static GeometryObject SplitWorldExtent(MultiPolygon mpoly)
{
MultiPolygon newmpoly = new MultiPolygon();
foreach (Polygon poly in mpoly.Polygons)
{
var geom = SplitWorldExtent(poly);
if (geom != null)
{
if (geom is Polygon)
{
newmpoly.Polygons.Add((Polygon)geom);
}
if (geom is MultiPolygon)
{
newmpoly.Polygons.AddRange(((MultiPolygon)geom).Polygons);
}
}
}
if (newmpoly.Polygons.Count == 1)
{
return(newmpoly.Polygons[0]);
}
return(newmpoly);
}
public static GeometryObject SplitWorldExtent(Polygon poly)
{
MultiPolygon newpoly = new MultiPolygon();
List <LineString> OuterlineStrings = new List <LineString>();
List <LineString> InnerlineStrings = new List <LineString>();
OuterlineStrings.AddRange(SplitWorldExtent(poly.LineStrings.First()));
foreach (LineString lineString in poly.LineStrings.Skip(1))
{
InnerlineStrings.AddRange(SplitWorldExtent(lineString));
}
if (OuterlineStrings.Count == 1)
{
return(poly);
}
List <Polygon> newpolygons = new List <Polygon>();
if (OuterlineStrings.Count > 0)
{
foreach (LineString ls in OuterlineStrings)
{
Polygon curpoly = new Polygon();
newpoly.Polygons.Add(curpoly);
curpoly.LineStrings.Add(ls);
foreach (var innerls in InnerlineStrings)
{
NetTopologySuite.IO.WKTReader wktreader = new NetTopologySuite.IO.WKTReader();
var innerlr = wktreader.Read(innerls.ToWkt());
NetTopologySuite.Geometries.Polygon polygon = (NetTopologySuite.Geometries.Polygon)wktreader.Read(curpoly.ToWkt());
if (polygon.Contains(innerlr))
{
curpoly.LineStrings.Add(innerls);
}
}
newpolygons.Add(curpoly);
}
}
return(newpoly);
}
static string GeometryToWktString(MultiPolygon multiPolygon)
{
return string.Format("({0})", string.Join(",", multiPolygon.Polygons.Select(GeometryToWktString)));
}
public static MultiSurfaceType ToGmlMultiSurface(this Terradue.GeoJson.Geometry.GeometryObject geometry)
{
if (geometry is Polygon) {
List<Polygon> polygons = new List<Polygon>();
polygons.Add((Polygon)geometry);
MultiPolygon multiPolygon = new MultiPolygon(polygons);
return ToGmlMultiSurface((MultiPolygon)multiPolygon);
} else if (geometry is MultiPolygon) {
return ToGmlMultiSurface((MultiPolygon)geometry);
} else
return null;
}
static string GeometryToWktString(MultiPolygon multiPolygon)
{
return(string.Format("({0})", string.Join(",", multiPolygon.Polygons.Select(GeometryToWktString))));
}
public static GeometryObject SplitWorldExtent(Polygon poly)
{
MultiPolygon newpoly = new MultiPolygon();
List<LineString> OuterlineStrings = new List<LineString>();
List<LineString> InnerlineStrings = new List<LineString>();
OuterlineStrings.AddRange(SplitWorldExtent(poly.LineStrings.First()));
foreach (LineString lineString in poly.LineStrings.Skip(1)) {
InnerlineStrings.AddRange(SplitWorldExtent(lineString));
}
if (OuterlineStrings.Count == 1) {
return poly;
}
List<Polygon> newpolygons = new List<Polygon>();
if (OuterlineStrings.Count > 0) {
foreach (LineString ls in OuterlineStrings) {
Polygon curpoly = new Polygon();
newpoly.Polygons.Add(curpoly);
curpoly.LineStrings.Add(ls);
foreach (var innerls in InnerlineStrings) {
NetTopologySuite.IO.WKTReader wktreader = new NetTopologySuite.IO.WKTReader();
var innerlr = wktreader.Read(innerls.ToWkt());
NetTopologySuite.Geometries.Polygon polygon = (NetTopologySuite.Geometries.Polygon)wktreader.Read(curpoly.ToWkt());
if (polygon.Contains(innerlr))
curpoly.LineStrings.Add(innerls);
}
newpolygons.Add(curpoly);
}
}
return newpoly;
}
// /// <summary>
// /// Get an XmlNameSpaceManager with the gmlnamesapces
// /// </summary>
// /// <returns>The xml namespace manager.</returns>
// public static XmlNamespaceManager GmlXmlNamespaceManager() {
// XmlNamespaceManager namespaces = new XmlNamespaceManager(new XmlDocument().NameTable);
// namespaces.AddNamespace("gml", "http://www.opengis.net/gml");
// namespaces.AddNamespace("gml32", "http://www.opengis.net/gml/3.2");
// return namespaces;
// }
//
// /// <summary>
// /// Initialize a new IGeometry object from a GML XML element
// /// </summary>
// /// <param name="gml">The GML XML element convert</param>
// public static GeometryObject GmlToGeometry(XmlElement gml) {
//
// if (gml == null)
// throw new ArgumentNullException("GML Xml Element is null!");
//
// if (gml.LocalName == "Point")
// return FromGMLPoint(gml);
//
// if (gml.LocalName == "LineString")
// return FromGMLLineString(gml);
//
// if (gml.LocalName == "Curve")
// return FromGMLCurve(gml);
//
// if (gml.LocalName == "LinearRing")
// return FromGMLLinearRing(gml);
//
// if (gml.LocalName == "Polygon")
// return FromGMLPolygon(gml);
//
// if (gml.LocalName == "Triangle")
// return FromGMLTriangle(gml);
//
// if (gml.LocalName == "Surface")
// return FromGMLSurface(gml);
//
// if (gml.LocalName == "MultiPoint")
// return FromGMLMultiPoint(gml);
//
// if (gml.LocalName == "MultiLineString")
// return FromGMLMultiLineString(gml);
//
// if (gml.LocalName == "MultiCurve")
// return FromGMLMultiCurve(gml);
//
// if (gml.LocalName == "MultiPolygon")
// return FromGMLMultiPolygon(gml);
//
// if (gml.LocalName == "MultiSurface")
// return FromGMLMultiSurface(gml);
//
// if (gml.LocalName == "PolyhedralSurface")
// return FromGMLPolyhedralSurface(gml);
//
// if (gml.LocalName == "Tin" || gml.LocalName == "TriangulatedSurface")
// return FromGMLTin(gml);
//
// if (gml.LocalName == "MultiGeometry")
// return FromGMLMultiGeometry(gml);
//
// throw new InvalidFormatException("gml type " + gml.LocalName + "is not supported");
//
// }
//
// private static Point FromGMLPoint(XmlElement gml) {
// if (gml.ChildNodes.Count != 1)
// throw new InvalidFormatException("invalid GML representation: Point type must have 1 child");
//
// return new Point(FromGMLData(gml.ChildNodes)[0]);
// }
//
// private static LineString FromGMLLineString(XmlElement gml) {
// if (gml.ChildNodes.Count == 0)
// return new LineString(new List<IPosition>());
//
// List<IPosition> points = FromGMLData(gml.ChildNodes);
//
// if (points.Count < 2)
// throw new InvalidFormatException("invalid GML representation: LineString type must have at least 2 positions");
//
// return new LineString(points);
// }
//
// private static LineString FromGMLCurve(XmlElement gml) {
// bool found = false;
// string interpolation;
//
// /* Looking for gml:segments */
// foreach (XmlNode node in gml) {
// XmlElement element = (XmlElement)node;
// if (!IsGMLNamespace(element, false))
// continue;
// if (element.LocalName == "segments") {
// found = true;
// break;
// }
// }
//
// if (!found)
// throw new InvalidFormatException("invalid GML representation: Curve type must have at least 1 segments");
//
// List<IPosition> positions = new List<IPosition>();
//
// /* Processing each gml:LineStringSegment */
// foreach (XmlNode node in gml) {
// XmlElement element = (XmlElement)node;
// if (!IsGMLNamespace(element, false))
// continue;
// if (element.LocalName == "LineStringSegment")
// continue;
//
// /* GML SF is resticted to linear interpolation */
// interpolation = element.GetAttribute("interpolation");
// if (interpolation == "linear")
// throw new InvalidFormatException("invalid GML representation: interpolation of a curve cannot be linear");
//
// List<IPosition> pos = FromGMLData(element.ChildNodes);
//
// if (pos.Count < 2)
// throw new InvalidFormatException("invalid GML representation: curve segments must have at least 2 points");
//
// positions.AddRange(pos);
//
// }
//
// if (positions.Count == 0)
// throw new InvalidFormatException("invalid GML representation: curve is empty");
//
// /* Most common case, a single segment */
// if (positions.Count == 1)
// return new LineString(positions);
//
// /*
// * "The curve segments are connected to one another, with the end point
// * of each segment except the last being the start point of the next
// * segment" from ISO 19107:2003 -> 6.3.16.1 (p43)
// *
// * So we must aggregate all the segments into a single one and avoid
// * to copy the redundants points
// */
//
// if (positions.Count > 1) {
// // TODO Aggregate stuff */
// throw new NotImplementedException("not implemented GML representation: curve with multi segment");
// }
//
// return new LineString(positions);
//
// }
//
// private static Polygon FromGMLLinearRing(XmlElement gml) {
// List<IPosition> positions;
// positions = FromGMLData(gml.ChildNodes);
//
// if (positions.Count < 4 || IsClosed2D(positions))
// throw new InvalidFormatException("invalid GML representation: linearring is not a closed ring of minimum 4 positions");
//
// LineString linestring = new LineString(positions);
// List <LineString> linestrings = new List <LineString>();
// linestrings.Add(linestring);
//
// return new Polygon(linestrings);
// }
//
// private static Polygon FromGMLPolygon(XmlElement gml) {
// List<LineString> polygon = new List<LineString>();
// List<IPosition> positions = null;
//
// /* Polygon/outerBoundaryIs -> GML 2.1.2 */
// /* Polygon/exterior -> GML 3.1.1 */
// foreach (XmlNode node in gml.SelectNodes("*[local-name()='outerBoundaryIs'] | *[local-name()='exterior']")) {
//
// XmlElement element = (XmlElement)node;
//
// foreach (XmlNode node1 in element) {
// XmlElement element1 = (XmlElement)node1;
//
// if (!IsGMLNamespace(element, true))
// continue;
// if (element1.LocalName == "LinearRing") {
// positions = FromGMLData(element1.ChildNodes);
//
// if (positions.Count < 4 || !IsClosed2D(positions))
// throw new InvalidFormatException("invalid GML representation: polygon outer is not a closed ring of minimum 4 positions");
// }
// }
// }
//
// if (positions == null)
// throw new InvalidFormatException("invalid GML representation: polygon outer is empty");
//
// polygon.Add(new LineString(positions));
//
// /* Polygon/outerBoundaryIs -> GML 2.1.2 */
// /* Polygon/exterior -> GML 3.1.1 */
// foreach (XmlNode node in gml.SelectNodes("innerBoundaryIs | interior")) {
//
// XmlElement element = (XmlElement)node;
//
// foreach (XmlNode node1 in element) {
// XmlElement element1 = (XmlElement)node1;
//
// if (!IsGMLNamespace(element, true))
// continue;
// if (element1.LocalName == "LinearRing") {
// positions = FromGMLData(element1.ChildNodes);
//
// if (positions.Count < 4 || IsClosed2D(positions))
// throw new InvalidFormatException("invalid GML representation: polygon inner is not a closed ring of minimum 4 positions");
// }
// }
//
// polygon.Add(new LineString(positions));
// }
//
// return new Polygon(polygon);
// }
//
// private static Polygon FromGMLTriangle(XmlElement gml) {
// // TODO implement FromGMLTriangle
// throw new NotImplementedException("GML Triangle is not supported yet.");
// }
//
// private static Polygon FromGMLSurface(XmlElement gml) {
//
// /* Looking for gml:patches */
// XmlElement element = (XmlElement)gml.SelectSingleNode("patches");
// if (element == null)
// throw new InvalidFormatException("invalid GML representation: surface must have patches");
//
// if (element.SelectNodes("PolygonPatch").Count > 1) {
// throw new NotImplementedException("Multi patch is not supported yet.");
// }
//
// return FromGMLPatch(element);
//
// }
//
// private static Polygon FromGMLPatch(XmlElement gml) {
// /* PolygonPatch */
// if (gml.LocalName != "PolygonPatch")
// throw new InvalidFormatException("invalid GML representation: patch must be PolygonPatch");
//
// /* GML SF is resticted to planar interpolation */
// string interpolation = gml.GetAttribute("interpolation");
// if (interpolation != null) {
// if (interpolation != "planar")
// throw new InvalidFormatException("invalid GML representation: interpolation must be planar");
// }
//
// return FromGMLPolygon(gml);
// }
//
// /// <summary>
// static private MultiPoint FromGMLMultiPoint(XmlElement gml) {
// List<IPosition> points = new List<IPosition>();
//
// if (gml.ChildNodes == null)
// return new MultiPoint(null);
//
// /* MultiPoint/pointMember */
// foreach (XmlNode node in gml.SelectNodes("gml:pointMember",GmlXmlNamespaceManager())) {
// points.AddRange(FromGMLData(node.ChildNodes));
//
// }
// /* MultiPoint/pointMember */
// foreach (XmlNode node in gml.SelectNodes("gml32:pointMember",GmlXmlNamespaceManager())) {
// points.AddRange(FromGMLData(node.ChildNodes));
//
// }
//
// return new MultiPoint(points);
// }
//
// private static MultiLineString FromGMLMultiLineString(XmlElement gml) {
//
// List<LineString> linestrings = new List<LineString>();
//
// /* MultiLineString/lineStringMember */
// foreach (XmlNode node in gml.SelectNodes("gml:lineStringMember/*",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLLineString((XmlElement)node));
// }
//
// /* MultiLineString/lineStringMember */
// foreach (XmlNode node in gml.SelectNodes("gml32:lineStringMember/*",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLLineString((XmlElement)node));
// }
//
// return new MultiLineString(linestrings);
// }
//
// private static MultiLineString FromGMLMultiCurve(XmlElement gml) {
// List<LineString> linestrings = new List<LineString>();
//
// /* MutltiCurve/curveMember/LineString */
// foreach (XmlNode node in gml.SelectNodes("gml:curveMember/gml:LineString",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLLineString((XmlElement)node));
// }
// /* MutltiCurve/curveMember/LineString */
// foreach (XmlNode node in gml.SelectNodes("gml32:curveMember/gml32:LineString",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLLineString((XmlElement)node));
// }
//
// /* MutltiCurve/curveMember/Curve */
// foreach (XmlNode node in gml.SelectNodes("gml:curveMember/gml:Curve",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLCurve((XmlElement)node));
// }
// /* MutltiCurve/curveMember/Curve */
// foreach (XmlNode node in gml.SelectNodes("gml32:curveMember/gml32:Curve",GmlXmlNamespaceManager())) {
// linestrings.Add(FromGMLCurve((XmlElement)node));
// }
//
// return new MultiLineString(linestrings);
// }
//
// private static MultiPolygon FromGMLMultiPolygon(XmlElement gml) {
// List<Polygon> polygons = new List<Polygon>();
//
// /* MultiPolygon/polygonMember */
// foreach (XmlNode node in gml.SelectNodes("gml:polygonMember/*",GmlXmlNamespaceManager())) {
// polygons.Add(FromGMLPolygon((XmlElement)node));
// }
// /* MultiPolygon/polygonMember */
// foreach (XmlNode node in gml.SelectNodes("gml32:polygonMember/*",GmlXmlNamespaceManager())) {
// polygons.Add(FromGMLPolygon((XmlElement)node));
// }
//
// return new MultiPolygon(polygons);
// }
//
// private static MultiPolygon FromGMLMultiSurface(XmlElement gml) {
// List<Polygon> polygons = new List<Polygon>();
//
// /* MultiSurface/surfaceMember */
// foreach (XmlNode node in gml.SelectNodes("gml:surfaceMembers/*",GmlXmlNamespaceManager())) {
// polygons.Add(FromGMLPolygon((XmlElement)node));
// }
// /* MultiSurface/surfaceMember */
// foreach (XmlNode node in gml.SelectNodes("gml32:surfaceMembers/*",GmlXmlNamespaceManager())) {
// polygons.Add(FromGMLPolygon((XmlElement)node));
// }
//
// return new MultiPolygon(polygons);
// }
//
// private static MultiPolygon FromGMLPolyhedralSurface(XmlElement gml) {
// // TODO implement FromGMLPolyhedralSurface
// throw new NotImplementedException("GML PolyhedralSurface is not supported yet.");
// }
//
// private static Polygon FromGMLTin(XmlElement gml) {
// // TODO implement FromGMLTin
// throw new NotImplementedException("GML Tin is not supported yet.");
// }
//
// private static GeometryCollection FromGMLMultiGeometry(XmlElement gml) {
// List<GeometryObject> geometries = new List<GeometryObject>();
//
// /*
// * In GML 2.1.2 pointMember, lineStringMember and
// * polygonMember are parts of geometryMember
// * substitution group
// */
// foreach (XmlNode node in gml.SelectNodes("gml:pointMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml:lineStringMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml:polygonMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml:geometryMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml32:pointMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml32:lineStringMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml32:polygonMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
// foreach (XmlNode node in gml.SelectNodes("gml32:geometryMember/*",GmlXmlNamespaceManager())) {
// geometries.Add(GmlToGeometry((XmlElement)node));
// }
//
// return new GeometryCollection(geometries);
// }
//
// private static List<IPosition> FromGMLData(XmlNodeList nodes) {
// List<IPosition> position = new List<IPosition>();
//
// foreach (XmlNode node in nodes) {
// XmlElement element = (XmlElement)node;
//
// if (element.LocalName == "pos")
// position.Add(FromGMLPos(element));
//
// if (element.LocalName == "posList")
// position.AddRange(FromGMLPosList(element));
//
// if (element.LocalName == "coordinates")
// position.AddRange(FromGMLCoordinates(element));
//
// if (element.LocalName == "coord")
// position.Add(FromGMLCoord(element));
//
// if (element.LocalName == "pointRep" || element.LocalName == "pointProperty") {
//
// bool found = false;
//
// foreach (XmlNode node1 in element.ChildNodes) {
// XmlElement el = (XmlElement)node1;
//
// if (!IsGMLNamespace(el, true))
// continue;
//
// if (el.LocalName == "Point") {
// found = true;
// break;
// }
//
// }
// if (!found)
// throw new InvalidFormatException("invalid GML representation: pointRep or pointProperty type must have at least 1 child");
//
// return FromGMLData(element.ChildNodes);
//
// }
//
// }
// return position;
// }
//
// private static IPosition FromGMLPos(XmlElement gml) {
// IPosition position;
// int dim;
// string gmlpos;
//
// if (!IsGMLNamespace(gml, true))
// return null;
// if (gml.LocalName != "pos")
// return null;
//
// string dimension = gml.GetAttribute("srsDimension");
//
// if (dimension == "") /* in GML 3.0.0 it was dimension */
// dimension = gml.GetAttribute("dimension");
// if (dimension == "")
// dim = 2; /* We assume that we are in 2D */
// else {
// dim = int.Parse(dimension);
// if (dim < 2 || dim > 3)
// throw new InvalidFormatException("invalid GML representation: gml:pos dimension equals " + dim);
// }
//
// /* We retrieve gml:pos string */
// gmlpos = gml.InnerText.Trim();
//
// /* gml:pos pattern:
// * x1 y1
// * x1 y1 z1
// */
// string[] pos = gmlpos.Split(' ');
// position = new GeographicPosition(pos[0], pos[1], pos[2]);
// return position;
// }
//
// private static List<IPosition> FromGMLPosList(XmlElement gml) {
// List<IPosition> position = new List<IPosition>();
// int dim;
// string gmlpos, dimension;
//
// if (!IsGMLNamespace(gml, false))
// return null;
//
// dimension = gml.GetAttribute("srsDimension");
//
// if (dimension == "") /* in GML 3.0.0 it was dimension */
// dimension = gml.GetAttribute("dimension");
// if (dimension == "")
// dim = 2; /* We assume that we are in 2D */
// else {
// dim = int.Parse(dimension);
// if (dim < 2 || dim > 3)
// throw new InvalidFormatException("invalid GML representation: gml:pos dimension equals " + dim);
// }
//
// /* We retrieve gml:pos string */
// gmlpos = gml.InnerText.Trim();
//
// /* gml:posList pattern: y1 x1 y2 x2
// * y1 x1 z1 y2 x2 z2
// */
// string[] pos = gmlpos.Split(' ');
//
// for (int i = 0; i < pos.Length; i += dim) {
// if (dim == 2)
// position.Add(new GeographicPosition(pos[i + 0], pos[i + 1], null));
// if (dim == 3)
// position.Add(new GeographicPosition(pos[i + 0], pos[i + 1], pos[i + 2]));
// }
// return position;
// }
//
// private static IPosition FromGMLCoord(XmlElement element) {
// double? x = null, y = null, z = null;
//
// foreach (XmlElement el in element.ChildNodes) {
// if (!IsGMLNamespace(el, false))
// continue;
//
// if (el.LocalName == "X")
// x = double.Parse(el.InnerText);
//
// if (el.LocalName == "Y")
// y = double.Parse(el.InnerText);
//
// if (el.LocalName == "Z")
// z = double.Parse(el.InnerText);
// }
//
// /* Check dimension consistancy */
// if (x == null || y == null)
// throw new InvalidFormatException("invalid GML representation: gml:coord missing X or Y");
//
// GeographicPosition geopos = new GeographicPosition(y.Value, x.Value, z.Value);
// return geopos;
//
// }
//
// private static List<IPosition> FromGMLCoordinates(XmlElement element) {
// List<IPosition> position = new List<IPosition>();
// string gmlcoord, gmlts, gmlcs, gmldec;
// char ts = ' ', cs = ',', dec = '.';
//
// gmlcoord = element.InnerText;
//
// /* Retrieve separator between coordinates tuples */
// gmlts = element.GetAttribute("ts");
// if (char.TryParse(gmlts, out ts) != true)
// ts = ' ';
//
// /* Retrieve separator between each coordinate */
// gmlcs = element.GetAttribute("cs");
// if (char.TryParse(gmlcs, out cs) != true) {
// cs = ',';
// }
//
// /* Retrieve decimal separator */
// gmldec = element.GetAttribute("decimal");
// if (char.TryParse(gmldec, out dec) != true)
// dec = '.';
//
// if (cs == ts || cs == dec || ts == dec)
// throw new InvalidFormatException("invalid GML representation: gml:coordinates ambiguity in separators");
//
// /* We retrieve gml:coord string */
// gmlcoord = element.InnerText.Trim();
// if (gmlcoord == null || gmlcoord == "")
// throw new InvalidFormatException("invalid GML representation: gml:coordinates is empty");
//
// string[] coordinates = gmlcoord.Split(ts);
//
// foreach (string coord in coordinates) {
//
// string[] pos = coord.Split(cs);
//
// double x = double.Parse(pos[1]);
// double y = double.Parse(pos[0]);
// double? z = null;
// if (pos.Length > 2)
// z = double.Parse(pos[2]);
//
// position.Add(new GeographicPosition(y, x, z));
//
// }
// return position;
//
// }
//
// private static bool IsGMLNamespace(XmlElement element, bool is_strict) {
// string ns = element.NamespaceURI;
//
// /*
// * If no namespace is available we could return true anyway
// * (because we work only on GML fragment, we don't want to
// * 'oblige' to add namespace on the geometry root node)
// */
// if (ns == null) {
// return !is_strict;
// }
//
// /*
// * Handle namespaces:
// * - http://www.opengis.net/gml (GML 3.1.1 and priors)
// * - http://www.opengis.net/gml/3.2 (GML 3.2.1)
// */
// if (ns == "http://www.opengis.net/gml" || ns == "http://www.opengis.net/gml/3.2")
// return true;
//
// return false;
// }
//
// private static bool IsClosed2D(List<IPosition> positions) {
// return positions[0].Equals(positions[positions.Count - 1]);
// }
public static GeometryObject SplitWorldExtent(MultiPolygon mpoly)
{
MultiPolygon newmpoly = new MultiPolygon();
foreach (Polygon poly in mpoly.Polygons) {
var geom = SplitWorldExtent(poly);
if (geom != null) {
if (geom is Polygon)
newmpoly.Polygons.Add((Polygon)geom);
if (geom is MultiPolygon)
newmpoly.Polygons.AddRange(((MultiPolygon)geom).Polygons);
}
}
if (newmpoly.Polygons.Count == 1)
return newmpoly.Polygons[0];
return newmpoly;
}
public static GeometryObject SplitWorldExtent(Polygon poly)
{
MultiPolygon newpoly = new MultiPolygon();
List<List<LineString>> lineStringss = new List<List<LineString>>();
foreach (LineString lineString in poly.LineStrings) {
lineStringss.Add(SplitWorldExtent(lineString));
}
if (poly.LineStrings.Count == lineStringss.First().Count) {
return poly;
}
List<List<LineString>> lineStringss2 = new List<List<LineString>>();
if (lineStringss.Count > 0) {
foreach (LineString ls in lineStringss[0]) {
var pls = new List<LineString>();
lineStringss2.Add(pls);
pls.Add(ls);
}
int i = 1;
while (i < lineStringss.Count) {
foreach (LineString ls in lineStringss[i]) {
lineStringss2[i].Add(ls);
}
i++;
}
}
foreach (var lss in lineStringss2) {
newpoly.Polygons.Add(new Polygon(lss));
}
return newpoly;
}
public static GeometryObject SplitWorldExtent(MultiPolygon mpoly)
{
foreach (Polygon poly in mpoly.Polygons) {
var mpolygon = SplitWorldExtent(poly);
if (mpolygon != null)
return mpolygon;
}
return mpoly;
}
public static XmlElement ToMultiSurfaceGml(this Terradue.GeoJson.Geometry.GeometryObject geometry)
{
XmlDocument xdoc = new XmlDocument();
using (XmlWriter writer = xdoc.CreateNavigator().AppendChild()) {
if (geometry is Polygon) {
List<Polygon> polygons = new List<Polygon>();
polygons.Add((Polygon)geometry);
MultiPolygon multiPolygon = new MultiPolygon(polygons);
GmlMultiSurfaceSerializer.Serialize(writer, GeometryToGml((MultiPolygon)multiPolygon));
} else if (geometry is MultiPolygon) {
GmlMultiSurfaceSerializer.Serialize(writer, GeometryToGml((MultiPolygon)geometry));
} else
return null;
}
return xdoc.DocumentElement;
}
public static Terradue.GeoJson.Gml.MultiSurfaceType GeometryToGml(MultiPolygon multiPolygon)
{
Terradue.GeoJson.Gml.MultiSurfaceType gmlMultiSurface = new Terradue.GeoJson.Gml.MultiSurfaceType();
gmlMultiSurface.surfaceMembers = new Terradue.GeoJson.Gml.SurfaceArrayPropertyType();
List<Terradue.GeoJson.Gml.PolygonType> gmlPolygons = new List<Terradue.GeoJson.Gml.PolygonType>();
foreach (var polygon in multiPolygon.Polygons) {
gmlPolygons.Add(GeometryToGml(polygon));
}
gmlMultiSurface.surfaceMembers.Items = gmlPolygons.ToArray();
return gmlMultiSurface;
}