public static OSGeo.OGR.Geometry ToGeometry(this Envelope envelope)
{
var ring = new OSGeo.OGR.Geometry(wkbGeometryType.wkbLinearRing);
ring.AddPoint_2D(envelope.MinX, envelope.MinY);
ring.AddPoint_2D(envelope.MaxX, envelope.MinY);
ring.AddPoint_2D(envelope.MaxX, envelope.MaxY);
ring.AddPoint_2D(envelope.MinX, envelope.MaxY);
ring.CloseRings();
var polygon = new OSGeo.OGR.Geometry(wkbGeometryType.wkbPolygon);
polygon.AddGeometry(ring);
return(polygon);
}
private static OSGeo.OGR.Geometry GetOgrGeometry(wkbGeometryType geometryType, IEnumerable <ICoordinate> coordinates)
{
OSGeo.OGR.Geometry ogrGeometry = null;
var firstCoord = coordinates?.FirstOrDefault();
if (firstCoord != null)
{
ogrGeometry = new OSGeo.OGR.Geometry(geometryType);
Action <ICoordinate> addCoordAction;
switch (firstCoord.Dimension)
{
case 2:
addCoordAction = (coordinate) => ogrGeometry.AddPoint_2D(coordinate.X, coordinate.Y);
break;
case 3:
addCoordAction = (coordinate) => ogrGeometry.AddPoint(coordinate.X, coordinate.Y, coordinate.Z);
break;
case 4:
addCoordAction = (coordinate) => ogrGeometry.AddPointZM(coordinate.X, coordinate.Y, coordinate.Z, coordinate.M);
break;
default:
return(ogrGeometry);
}
foreach (var item in coordinates)
{
addCoordAction(item);
}
}
return(ogrGeometry);
}
private static void AddPoint(this OSGeo.OGR.Geometry geometry, double x, double y, double z = double.NaN, double m = double.NaN)
{
if (geometry != null)
{
if (double.IsNaN(z))
{
if (double.IsNaN(m))
{
geometry.AddPoint_2D(x, y);
}
else
{
geometry.AddPointM(x, y, m);
}
}
else
{
if (double.IsNaN(m))
{
geometry.AddPoint(x, y, z);
}
else
{
geometry.AddPointZM(x, y, z, m);
}
}
}
}
public static OSGeo.OGR.Geometry ToOgrPoint(double x, double y, double z = double.NaN, double m = double.NaN)
{
OSGeo.OGR.Geometry ogrGeometry = new OSGeo.OGR.Geometry(wkbGeometryType.wkbPoint);
if (double.IsNaN(z))
{
ogrGeometry.AddPoint_2D(x, y);
}
else
{
if (double.IsNaN(m))
{
ogrGeometry.AddPoint(x, y, z);
}
else
{
ogrGeometry.AddPointZM(x, y, z, m);
}
}
return(ogrGeometry);
}
/// <summary>
/// 转换是否纠偏，只使用与火星坐标纠偏
/// </summary>
/// <param name="dataTable">数据表，如果是点表中含X,Y字段，如果是线和面，表中含有PATH字段</param>
/// <param name="filePath">存储路径</param>
/// <param name="geoType">数据类型，暂时支持简单的点线面</param>
/// <param name="jiupian">是否纠偏[火星坐标转84坐标]</param>
public static void SaveShpFile(DataTable dataTable, string filePath, wkbGeometryType geoType, ProjectConvert jiupian)
{
if (geoType == wkbGeometryType.wkbLineString || geoType == wkbGeometryType.wkbPolygon)
{
if (!dataTable.Columns.Contains("PATH"))
{
return;
}
}
else if (geoType == wkbGeometryType.wkbPoint)
{
if (!dataTable.Columns.Contains("X") || !dataTable.Columns.Contains("Y"))
{
return;
}
}
else
{
return;
}
// 为了支持中文路径，请添加下面这句代码
OSGeo.GDAL.Gdal.SetConfigOption("GDAL_FILENAME_IS_UTF8", "NO");
// 为了使属性表字段支持中文，请添加下面这句
OSGeo.GDAL.Gdal.SetConfigOption("SHAPE_ENCODING", "");
string strVectorFile = filePath;
// 注册所有的驱动
Ogr.RegisterAll();
//创建数据，这里以创建ESRI的shp文件为例
string strDriverName = "ESRI Shapefile";
int count = Ogr.GetDriverCount();
Driver oDriver = Ogr.GetDriverByName(strDriverName);
if (oDriver == null)
{
log.ErrorFormat("{0}驱动不可用！\n", strVectorFile);
Console.WriteLine("{0}驱动不可用！\n", strVectorFile);
return;
}
// 创建数据源
DataSource oDS = oDriver.CreateDataSource(strVectorFile, null);
if (oDS == null)
{
log.ErrorFormat("创建矢量文件{0}失败！\n", strVectorFile);
Console.WriteLine("创建矢量文件【%s】失败！\n", strVectorFile);
return;
}
string strwkt = "GEOGCS[\"WGS 84\",DATUM[\"WGS_1984\",SPHEROID[\"WGS 84\",6378137,298.257223563,AUTHORITY[\"EPSG\",\"7030\"]],AUTHORITY[\"EPSG\",\"6326\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.0174532925199433,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4326\"]]";
OSGeo.OSR.SpatialReference srs = new OSGeo.OSR.SpatialReference(strwkt);
// 创建图层，创建一个多边形图层，这里没有指定空间参考，如果需要的话，需要在这里进行指定
Layer oLayer = oDS.CreateLayer("test", null, geoType, null);
if (oLayer == null)
{
Console.WriteLine("图层创建失败！\n");
return;
}
for (int i = 0; i < dataTable.Columns.Count; i++)
{
if (dataTable.Columns[i].ColumnName == "PATH")
{
continue;
}
string columnName = dataTable.Columns[i].ColumnName;
FieldDefn oFieldName = new FieldDefn(columnName, FieldType.OFTString);
oFieldName.SetWidth(255);
oLayer.CreateField(oFieldName, 1);
}
FeatureDefn oDefn = oLayer.GetLayerDefn();
foreach (DataRow row in dataTable.Rows)
{
if ((geoType == wkbGeometryType.wkbPolygon || geoType == wkbGeometryType.wkbLineString) && string.IsNullOrEmpty(row["PATH"].ToString()))
{
continue;
}
// 创建Multi要素
Feature oFeatureMulty = new Feature(oDefn);
for (int i = 0; i < dataTable.Columns.Count; i++)
{
if (geoType == wkbGeometryType.wkbPolygon || geoType == wkbGeometryType.wkbLineString)
{
if (dataTable.Columns[i].ColumnName == "PATH")
{
continue;
}
}
string name = row[i] != null ? row[i].ToString() : "";
if (!string.IsNullOrEmpty(name))
{
try
{
//oFeatureMulty.SetField(dataTable.Columns[i].ColumnName,name);
var value = name.Trim().Replace("'", "").Replace(",", "");
if (value.Length > 255)
{
value = value.Substring(0, 250);
}
oFeatureMulty.SetField(i, value);
}
catch (Exception ex)
{
Console.WriteLine(ex);
}
}
}
if (geoType == wkbGeometryType.wkbPolygon || geoType == wkbGeometryType.wkbLineString)
{
string path = row["PATH"] != null ? row["PATH"].ToString() : "";
if (string.IsNullOrEmpty(path))
{
continue;
}
var geoCount = path.Split('|').Length;
wkbGeometryType newGeometryType = geoType;
OSGeo.OGR.Geometry geo1 = new OSGeo.OGR.Geometry(geoType);
if (geoType == wkbGeometryType.wkbPolygon)
{
if (geoCount > 1)
{
geo1 = new OSGeo.OGR.Geometry(wkbGeometryType.wkbMultiPolygon);
for (int i = 0; i < geoCount; i++)
{
var temp = path.Split('|')[i];
var points = temp.Split(';');
OSGeo.OGR.Geometry innGeometry = new Geometry(wkbGeometryType.wkbLinearRing);
for (int j = 0; j < points.Length; j++)
{
var c = ProjectChange(points[j], jiupian);
innGeometry.AddPoint_2D(c.lon, c.lat);
}
var polygon = new Geometry(wkbGeometryType.wkbPolygon);
polygon.AddGeometry(innGeometry);
geo1.AddGeometry(polygon);
}
}
else
{
geo1 = new OSGeo.OGR.Geometry(wkbGeometryType.wkbPolygon);
for (int i = 0; i < geoCount; i++)
{
var temp = path.Split('|')[i];
var points = temp.Split(';');
OSGeo.OGR.Geometry innGeometry = new Geometry(wkbGeometryType.wkbLinearRing);
for (int j = 0; j < points.Length; j++)
{
var c = ProjectChange(points[j], jiupian);
innGeometry.AddPoint_2D(c.lon, c.lat);
}
geo1.AddGeometry(innGeometry);
}
}
}
else
{
if (geoCount > 1)
{
geo1 = new OSGeo.OGR.Geometry(wkbGeometryType.wkbMultiLineString);
for (int i = 0; i < geoCount; i++)
{
var temp = path.Split('|')[i];
var points = temp.Split(';');
OSGeo.OGR.Geometry innGeometry = new Geometry(wkbGeometryType.wkbLinearRing);
for (int j = 0; j < points.Length; j++)
{
var c = ProjectChange(points[j], jiupian);
innGeometry.AddPoint_2D(c.lon, c.lat);
}
geo1.AddGeometry(innGeometry);
}
}
else
{
for (int i = 0; i < geoCount; i++)
{
var temp = path.Split('|')[i];
var points = temp.Split(';');
geo1 = new OSGeo.OGR.Geometry(wkbGeometryType.wkbLinearRing);
for (int j = 0; j < points.Length; j++)
{
var c = ProjectChange(points[j], jiupian);
geo1.AddPoint_2D(c.lon, c.lat);
}
}
}
}
//if (geoType == wkbGeometryType.wkbLineString)
//{
// for (int i = 0; i < ps.Length; i++)
// {
// var c = ProjectChange(ps[i], jiupian);
// geo1.AddPoint_2D(c.lon, c.lat);
// }
//}
//else if (geoType == wkbGeometryType.wkbPolygon)
//{
// List<int> indexs = new List<int>();
// indexs.Add(0);
// for (int i = 0; i < ps.Length; i++)
// {
// if (ps[i].Split('|', ',').Length > 2)
// {
// indexs.Add(i);
// }
// }
// indexs.Add(ps.Length - 1);
// if (indexs.Count == 2)
// {
// OSGeo.OGR.Geometry outGeo = new OSGeo.OGR.Geometry(OSGeo.OGR.wkbGeometryType.wkbLinearRing);
// for (int i = 0; i < ps.Length; i++)
// {
// var c = ProjectChange(ps[i], jiupian);
// outGeo.AddPoint_2D(c.lon, c.lat);
// }
// var c1 = ProjectChange(ps[0], jiupian);
// outGeo.AddPoint_2D(c1.lon, c1.lat);
// geo1.AddGeometryDirectly(outGeo);
// }
// else
// {
// for (int j = 0; j < indexs.Count - 1; j++)
// {
// OSGeo.OGR.Geometry outGeo = new OSGeo.OGR.Geometry(OSGeo.OGR.wkbGeometryType.wkbLinearRing);
// Coord sc = new Coord();
// int first = indexs[j];
// int end = indexs[j + 1];
// for (int i = first; i < end + 1; i++)
// {
// if (i == first && first != 0)
// {
// var c = ProjectChange(ps[i], jiupian);
// outGeo.AddPoint(c.lon, c.lat,0);
// }
// else
// {
// if (jiupian != ProjectConvert.NONE)
// {
// var c = ProjectChange(ps[i], jiupian);
// outGeo.AddPoint(c.lon, c.lat,0);
// }
// else
// {
// double x = double.Parse(ps[i].Split('|', ',')[0]);
// double y = double.Parse(ps[i].Split('|', ',')[1]);
// if (i == first)
// {
// sc.lon = x;
// sc.lat = y;
// }
// outGeo.AddPoint(x, y, 0);
// }
// }
// }
// outGeo.AddPoint(sc.lon, sc.lat, 0);
// geo1.AddGeometryDirectly(outGeo);
// }
// }
//}
oFeatureMulty.SetGeometry(geo1);
}
else if (geoType == wkbGeometryType.wkbPoint)
{
OSGeo.OGR.Geometry geo1 = new OSGeo.OGR.Geometry(wkbGeometryType.wkbPoint);
string xstring = row["X"] != null ? row["X"].ToString() : "";
string ystring = row["Y"] != null ? row["Y"].ToString() : "";
double x = double.Parse(xstring);
double y = double.Parse(ystring);
Coord c = CoordHelper.Gcj2Wgs(x, y);
geo1.AddPoint(c.lon, c.lat, 0);
oFeatureMulty.SetGeometry(geo1);
}
oLayer.CreateFeature(oFeatureMulty);
}
oLayer.Dispose();
oDS.Dispose();
}
