private static void _标准差(string dzx, out double aue, out double bzc)
{
OSGeo.OGR.Ogr.RegisterAll();
OSGeo.OGR.Driver dr = OSGeo.OGR.Ogr.GetDriverByName("ESRI shapefile");
OSGeo.OGR.DataSource dzxDS = dr.Open(dzx, 0);
OSGeo.OGR.Layer dzxLayer = dzxDS.GetLayerByIndex(0);
//获取Featuer数
int featCount = dzxLayer.GetFeatureCount(0);
// 1 拿到每个Featuer的Value
double[] values = new double[featCount];
for (int i = 0; i < featCount; i++)
{
OSGeo.OGR.Feature fileFeat = dzxLayer.GetFeature(i);
OSGeo.OGR.Geometry fileGeom = fileFeat.GetGeometryRef();
values[i] = fileFeat.GetFieldAsDouble("EVE");
fileGeom.Dispose();
fileFeat.Dispose();
}
dzxDS.Dispose();
// 2 求Values的平均值
aue = values.Average();
// 3 求values与平均值差的平方和
double pingFangHe = 0;
for (int i = 0; i < featCount; i++)
{
pingFangHe += (values[i] - aue) * (values[i] - aue);
}
// 4 每个值与平均值的差相加,除Featuer数.再开方,得到标准差
bzc = Math.Sqrt(pingFangHe / featCount);
}
/// <summary>
/// 标准差
/// </summary>
/// <param name="dzxLayer"></param>
/// <param name="aue"></param>
/// <param name="bzc"></param>
public static void getBZC(OSGeo.OGR.Layer dzxLayer, out double aue, out double bzc)
{
//获取Featuer数
int featCount = dzxLayer.GetFeatureCount(0);
// 1 拿到每个Featuer的Value
double[] values = new double[featCount];
for (int i = 0; i < featCount; i++)
{
OSGeo.OGR.Feature fileFeat = dzxLayer.GetFeature(i);
values[i] = fileFeat.GetFieldAsDouble("EVE");
fileFeat.Dispose();
}
// 2 求Values的平均值
aue = values.Average();
// 3 求values与平均值差的平方和
double pingFangHe = 0;
for (int i = 0; i < featCount; i++)
{
pingFangHe += (values[i] - aue) * (values[i] - aue);
}
// 4 每个值与平均值的差相加,除Featuer数.再开方,得到标准差
bzc = Math.Sqrt(pingFangHe / featCount);
}
private static void getH(OSGeo.OGR.Layer oriLayer)
{
//创建一个BUFFER，BUFFER距离为1米
OSGeo.OGR.Layer bufLayer = bufferFile(oriLayer, 1);
//判断原文件中是否有以下字段，没有就创建
if (oriLayer.FindFieldIndex("MIN", 1) == -1)
{
OSGeo.OGR.FieldDefn min = new OSGeo.OGR.FieldDefn("MIN", OSGeo.OGR.FieldType.OFTReal);
oriLayer.CreateField(min, 1);
}
if (oriLayer.FindFieldIndex("MAX", 1) == -1)
{
OSGeo.OGR.FieldDefn max = new OSGeo.OGR.FieldDefn("MAX", OSGeo.OGR.FieldType.OFTReal);
oriLayer.CreateField(max, 1);
}
if (oriLayer.FindFieldIndex("HIGHT", 1) == -1)
{
OSGeo.OGR.FieldDefn hight = new OSGeo.OGR.FieldDefn("HIGHT", OSGeo.OGR.FieldType.OFTReal);
oriLayer.CreateField(hight, 1);
}
if (oriLayer.FindFieldIndex("TEMP", 1) == -1)
{
OSGeo.OGR.FieldDefn hight = new OSGeo.OGR.FieldDefn("TEMP", OSGeo.OGR.FieldType.OFTReal);
oriLayer.CreateField(hight, 1);
}
double[] transfrom = new double[6];
dsmDataset.GetGeoTransform(transfrom);
int allX = dsmDataset.RasterXSize;
int allY = dsmDataset.RasterYSize;
StaticTools.msgLine("getH...");
//开始计算每个Feature需要读取的栅格参数
int featCount = oriLayer.GetFeatureCount(0);
for (int i = 0; i < featCount; i++)
{
OSGeo.OGR.Feature oriFeat = oriLayer.GetFeature(i);
OSGeo.OGR.Feature bufFeat = bufLayer.GetFeature(i);
int[] subRasterOff_Size = subRasterInfo(transfrom, allX, allY, bufFeat);
// getMaxMinValue(dsmDataset, oriFeat, bufFeat, subRasterOff_Size);
oriLayer.SetFeature(oriFeat);
oriFeat.Dispose();
bufFeat.Dispose();
StaticTools.progress((i + 1) * 100 / featCount, $"{i} / {featCount}");
}
if (IsDelete)
{
shpDataSet.deleteLayerByName(bufLayer.GetName());
}
}
public static void claenPoint(this OSGeo.OGR.Layer resLayer, double jiaodu, int cishu)
{
int featCount = resLayer.GetFeatureCount(0);
for (int i = 0; i < featCount; i++)
{
OSGeo.OGR.Feature oriFeat = resLayer.GetFeature(i);
OSGeo.OGR.Geometry oriGeom = oriFeat.GetGeometryRef();
OSGeo.OGR.Geometry subGeom = oriGeom.GetGeometryRef(0);
int pointCount = subGeom.GetPointCount();
Point[] aFeat = new Point[pointCount];
for (int c = 0; c < pointCount; c++)
{
aFeat[c].X = subGeom.GetX(c);
aFeat[c].Y = subGeom.GetY(c);
aFeat[c].Z = subGeom.GetZ(c);
}
OSGeo.OGR.Geometry newGeom = null;
if (aFeat.Length > cishu * 3)
{
newGeom = JID(aFeat, jiaodu, cishu);
}
else
{
oriFeat.Dispose();
continue;
}
if (newGeom != null)
{
oriFeat.SetGeometry(newGeom);
resLayer.SetFeature(oriFeat);
}
oriFeat.Dispose();
}
}
/******************************** 等值线转换为多边形 ***********************************************/
/// <summary>
/// 等值线转为POLYGON
/// </summary>
/// <param name="filePath"></param>
/// <returns></returns>
public static string dzPoly(string filePath)
{
OSGeo.OGR.Ogr.RegisterAll();
OSGeo.OGR.Driver dr = OSGeo.OGR.Ogr.GetDriverByName("ESRI shapefile");
string a = StaticTools.tempFilePath("shp", "等值线POLY");
OSGeo.OGR.DataSource newDS = dr.CreateDataSource(a, null);
OSGeo.OGR.Layer polyLayer = newDS.CreateLayer("", null, OSGeo.OGR.wkbGeometryType.wkbPolygon, null);
OSGeo.OGR.FieldDefn fieldDf0 = new OSGeo.OGR.FieldDefn("LID", OSGeo.OGR.FieldType.OFTInteger);
OSGeo.OGR.FieldDefn fieldDf1 = new OSGeo.OGR.FieldDefn("EVE", OSGeo.OGR.FieldType.OFTReal);
polyLayer.CreateField(fieldDf0, 1); //ID
polyLayer.CreateField(fieldDf1, 1); //Value
OSGeo.OGR.FeatureDefn featDF = new OSGeo.OGR.FeatureDefn("");
Console.WriteLine("开始等值线转POLY！");
OSGeo.OGR.DataSource cleanDS = dr.Open(filePath, 0);
OSGeo.OGR.Layer cleanLayer = cleanDS.GetLayerByIndex(0);
for (int i = 0; i < cleanLayer.GetFeatureCount(0); i++)
{
OSGeo.OGR.Feature lineFeat = cleanLayer.GetFeature(i);
OSGeo.OGR.Geometry lineGeom = lineFeat.GetGeometryRef();
OSGeo.OGR.Feature polyFeat = new OSGeo.OGR.Feature(featDF);
OSGeo.OGR.Geometry polyGeom = new OSGeo.OGR.Geometry(OSGeo.OGR.wkbGeometryType.wkbPolygon);
OSGeo.OGR.Geometry subGeom = new OSGeo.OGR.Geometry(OSGeo.OGR.wkbGeometryType.wkbLinearRing);
int u = lineGeom.GetPointCount();
for (int s = 0; s < u; s++)
{
double x = lineGeom.GetX(s);
double y = lineGeom.GetY(s);
double z = lineGeom.GetZ(s);
subGeom.AddPoint(x, y, z);
}
polyGeom.AddGeometry(subGeom);
polyFeat.SetGeometry(polyGeom);
polyLayer.CreateFeature(polyFeat);
lineGeom.Dispose();
polyGeom.Dispose();
subGeom.Dispose();
lineFeat.Dispose();
polyFeat.Dispose();
}
cleanLayer.Dispose();
polyLayer.Dispose();
cleanDS.Dispose();
newDS.Dispose();
Console.WriteLine("等值线转POLY完成！");
return(a);
}
/// <summary>
/// 清理Layer中重复的，面积过小的featuer,返回一个新的Layer
/// </summary>
/// <param name="inLayer">源layer</param>
/// <param name="andArea">是否判断最小面积</param>
/// <param name="minArea">最小面积</param>
/// <returns></returns>
public static OSGeo.OGR.Layer cleanLayer_FF(OSGeo.OGR.Layer inLayer)
{
string oldLayerName = inLayer.GetName();
string newLayerName = oldLayerName + "Clear";
shpDataSet.deleteLayerByName(newLayerName);
OSGeo.OGR.Layer outLayer = shpDataSet.CreateLayer(newLayerName, srs, inLayer.GetGeomType(), null);
int featCount = inLayer.GetFeatureCount(0);
StaticTools.msgLine($"clean resLine... before {featCount}");
for (int i = 0; i < featCount - 1; i++)
{
bool isOnly = true;
OSGeo.OGR.Feature ori = inLayer.GetFeature(i);
OSGeo.OGR.Envelope oriEnve = new OSGeo.OGR.Envelope();
ori.GetGeometryRef().GetEnvelope(oriEnve);
double maxCha = 1;
for (int j = i + 1; j < featCount; j++)
{
OSGeo.OGR.Feature next = inLayer.GetFeature(j);
OSGeo.OGR.Envelope nextEnve = new OSGeo.OGR.Envelope();
next.GetGeometryRef().GetEnvelope(nextEnve);
if (Math.Abs(oriEnve.MaxX - nextEnve.MaxX) < maxCha &&
Math.Abs(oriEnve.MaxY - nextEnve.MaxY) < maxCha &&
Math.Abs(oriEnve.MinX - nextEnve.MinX) < maxCha &&
Math.Abs(oriEnve.MinY - nextEnve.MinY) < maxCha)
{
isOnly = false;
break;
}
}
if (isOnly)
{
outLayer.CreateFeature(ori);
}
ori.Dispose();
StaticTools.progress((i + 2) * 100 / featCount, $"{i + 1} / {featCount}");
}
if (IsDelete)
{
shpDataSet.deleteLayerByName(inLayer.GetName());
}
StaticTools.msgLine($"clean resLine... after {outLayer.GetFeatureCount(0)}");
return(outLayer);
}
/************************************** 清理等值线 *****************************************/
/// <summary>
/// 通过线长,值 清理等值线
/// </summary>
/// <param name="filePath"></param>
public static string cleanDS(string dzx)
{
Console.WriteLine("开始清理等值线！");
double aue, bzc; _标准差(dzx, out aue, out bzc);
double minLength = 50;
double maxLength = 2600;
double minValue = aue - bzc * 2;
double maxValue = aue + bzc * 2;
//open dzx
OSGeo.OGR.Ogr.RegisterAll();
OSGeo.OGR.Driver dr = OSGeo.OGR.Ogr.GetDriverByName("ESRI shapefile");
OSGeo.OGR.DataSource dzxDS = dr.Open(dzx, 0);
OSGeo.OGR.Layer dzxLayer = dzxDS.GetLayerByIndex(0);
//new a shp
string cleanline = StaticTools.tempFilePath("shp", "清理后的等值线");
OSGeo.OGR.DataSource newdzxDS = dr.CreateDataSource(cleanline, null);
OSGeo.OGR.Layer newdzxLayer = newdzxDS.CreateLayer(dzxLayer.GetName(), dzxLayer.GetSpatialRef(), dzxLayer.GetGeomType(), null);
for (int i = 0; i < dzxLayer.GetFeatureCount(0); i++)
{
OSGeo.OGR.Feature fileFeat = dzxLayer.GetFeature(i);
OSGeo.OGR.Geometry fileGeom = fileFeat.GetGeometryRef();
double FeatLength = fileGeom.Length();
bool s1 = FeatLength > minLength && FeatLength < maxLength;
double featValue = fileFeat.GetFieldAsDouble("EVE");
bool s2 = featValue > minValue && featValue < maxValue;
bool isR = fileGeom.IsRing();
if (s1 && s2 && isR)
{
newdzxLayer.CreateFeature(fileFeat);
}
fileFeat.Dispose();
}
newdzxDS.Dispose();
dzxDS.Dispose();
Console.WriteLine("清理等值线完成！");
return(cleanline);
}
/// <summary>
/// 创建一个Buffer file,用于获取最小高度,juli是buffer的距离
/// </summary>
/// <param name="infile"></param>
/// <param name="juli"></param>
/// <returns></returns>
private static OSGeo.OGR.Layer bufferFile(OSGeo.OGR.Layer inLayer, int juli)
{
string buf = inLayer.GetName() + "buf";
shpDataSet.deleteLayerByName(buf);
OSGeo.OGR.Layer bufferLayer = shpDataSet.CreateLayer(buf, inLayer.GetSpatialRef(), inLayer.GetGeomType(), null);
int featCount = inLayer.GetFeatureCount(0);
for (int i = 0; i < featCount; i++)
{
OSGeo.OGR.Feature inFeat = inLayer.GetFeature(i);
OSGeo.OGR.Geometry inGeom = inFeat.GetGeometryRef();
OSGeo.OGR.Geometry outGeom = inGeom.Buffer(juli, 0);
OSGeo.OGR.Feature outFeat = new OSGeo.OGR.Feature(new OSGeo.OGR.FeatureDefn(""));
outFeat.SetGeometry(outGeom);
bufferLayer.CreateFeature(outFeat);
inFeat.Dispose();
outFeat.Dispose();
}
return(bufferLayer);
}
private void pdfToolStripMenuItem_Click(object sender, EventArgs e)
{
for (int i = 0; i < mapBox1.Map.Layers.Count; i++)
{
if (mapBox1.Map.Layers[i].LayerName == SelectedLayer)
{
var layer = mapBox1.Map.Layers[i] as VectorLayer;
var dataSource = layer.DataSource as GeometryFeatureProvider;
OSGeo.OGR.Driver driverSH = OSGeo.OGR.Ogr.GetDriverByName("PDF");
var referance = new SpatialReference("PROJCS[\"WGS 84 / Pseudo-Mercator\",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\"]],PROJECTION[\"Mercator_1SP\"],PARAMETER[\"central_meridian\",0],PARAMETER[\"scale_factor\",1],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH],EXTENSION[\"PROJ4\",\"+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs\"],AUTHORITY[\"EPSG\",\"3857\"]]");
OSGeo.OGR.DataSource dataSourceSH = driverSH.CreateDataSource(Application.StartupPath + "\\" + SelectedLayer + ".pdf", new string[] { "ENCODING=UTF-16" });
var layerSH = dataSourceSH.CreateLayer("PointLayer", referance, OSGeo.OGR.wkbGeometryType.wkbPoint, new string[] { "ENCODING=UTF-16" });
layerSH.StartTransaction();
for (int k = 0; k < dataSource.Features.Columns.Count; k++)
{
OSGeo.OGR.FieldDefn field = new OSGeo.OGR.FieldDefn(dataSource.Features.Columns[k].Caption, OSGeo.OGR.FieldType.OFTString);
layerSH.CreateField(field, 1);
}
for (int k = 0; k < dataSource.Features.Rows.Count; k++)
{
OSGeo.OGR.Feature feature = new OSGeo.OGR.Feature(layerSH.GetLayerDefn());
if (!string.IsNullOrEmpty(dataSource.Features.Rows[k]["geom_wkt"].ToString()))
{
OSGeo.OGR.Geometry geom = OSGeo.OGR.Geometry.CreateFromWkt(dataSource.Features.Rows[k]["geom_wkt"].ToString());
feature.SetGeometry(geom);
}
for (int x = 0; x < dataSource.Features.Columns.Count; x++)
{
var data = dataSource.Features.Rows[k][x].ToString();
feature.SetField(dataSource.Features.Columns[x].Caption, data);
}
try
{
layerSH.CreateFeature(feature);
feature.Dispose();
}
catch (Exception)
{
}
}
layerSH.CommitTransaction();
layerSH.SyncToDisk();
layerSH.Dispose();
dataSourceSH.FlushCache();
dataSourceSH.Dispose();
driverSH.Dispose();
}
}
}
/// <summary>
/// Creates an OGR data source from a FeatureDataTable
/// </summary>
/// <param name="table">The name of the table</param>
/// <param name="geometryType">The geometry type</param>
/// <param name="driver">The driver</param>
/// <param name="connection">The connection string</param>
/// <param name="driverOptions">The options for the driver</param>
/// <param name="layerOptions">The options for the layer</param>
public static void CreateFromFeatureDataTable(FeatureDataTable table,
OgcGeometryType geometryType, int srid, string driver, string connection, string[] driverOptions = null, string[] layerOptions = null)
{
if (table == null)
{
throw new ArgumentNullException("table");
}
if (table.Rows.Count == 0)
{
throw new ArgumentException("The table contains no rows", "table");
}
if (geometryType < OgcGeometryType.Point || geometryType > OgcGeometryType.MultiPolygon)
{
throw new ArgumentException("Invalid geometry type", "geometryType");
}
if (string.IsNullOrWhiteSpace(driver))
{
throw new ArgumentException("No driver specified", "driver");
}
var dr = OSGeo.OGR.Ogr.GetDriverByName(driver);
if (dr == null)
{
throw new Exception(string.Format("Cannot load driver '{0}'!", driver));
}
//if (!dr.TestCapability("ODrCCreateDataSource"))
// throw new Exception(string.Format("Driver '{0}' cannot create a data source!", driver));
// Create the data source
var ds = dr.CreateDataSource(connection, driverOptions);
//if (!ds.TestCapability("ODsCCreateLayer"))
// throw new Exception(string.Format("Driver '{0}' cannot create a layer!", driver));
// Create the spatial reference
var sr = new OSGeo.OSR.SpatialReference(string.Empty);
sr.ImportFromEPSG(srid);
// Create the layer
var lyr = ds.CreateLayer(table.TableName, sr, (OgrGeometryType)geometryType, layerOptions);
sr.Dispose();
//lyr.GetSpatialRef();
foreach (System.Data.DataColumn dc in table.Columns)
{
using (var fldDef = GetFieldDefinition(dc))
lyr.CreateField(fldDef, 0);
}
using (var ld = lyr.GetLayerDefn())
{
foreach (FeatureDataRow fdr in table.Rows)
{
if ((int)fdr.Geometry.OgcGeometryType != (int)geometryType)
{
continue;
}
using (var feature = new OgrFeature(ld))
{
feature.SetGeometry(OgrGeometry.CreateFromWkb(fdr.Geometry.AsBinary()));
var idx = -1;
foreach (System.Data.DataColumn dc in table.Columns)
{
idx++;
var fd = ld.GetFieldDefn(idx);
DateTime dt;
switch (fd.GetFieldType())
{
case OgrFieldType.OFTBinary:
//Nothing
break;
case OgrFieldType.OFTDate:
dt = ((DateTime)fdr[dc]).Date;
feature.SetField(idx, dt.Year, dt.Month, dt.Day, dt.Hour, dt.Minute, dt.Second, 0);
break;
case OgrFieldType.OFTDateTime:
dt = (DateTime)fdr[dc];
feature.SetField(idx, dt.Year, dt.Month, dt.Day, dt.Hour, dt.Minute, dt.Second, 0);
break;
case OgrFieldType.OFTTime:
var tod = ((DateTime)fdr[dc]).TimeOfDay;
feature.SetField(idx, 0, 0, 0, tod.Hours, tod.Minutes, tod.Seconds, 0);
break;
case OgrFieldType.OFTInteger:
feature.SetField(idx, Convert.ToInt32(fdr[dc]));
break;
case OgrFieldType.OFTIntegerList:
var il = GetIntegerList(fdr[dc], dc.DataType);
feature.SetFieldIntegerList(idx, il.Length, il);
break;
case OgrFieldType.OFTReal:
feature.SetField(idx, Convert.ToDouble(fdr[dc]));
break;
case OgrFieldType.OFTRealList:
var dl = GetDoubleList(fdr[dc], dc.DataType);
feature.SetFieldDoubleList(idx, dl.Length, dl);
break;
case OgrFieldType.OFTString:
feature.SetField(idx, Convert.ToString(fdr[dc]));
break;
case OgrFieldType.OFTStringList:
var sl = (string[])fdr[dc];
feature.SetFieldStringList(idx, sl);
break;
}
fd.Dispose();
}
lyr.CreateFeature(feature);
feature.Dispose();
}
//ld.Dispose();
}
}
lyr.Dispose();
ds.Dispose();
dr.Dispose();
}
/// <summary>
/// Creates an OGR data source from a FeatureDataTable
/// </summary>
/// <param name="table">The name of the table</param>
/// <param name="geometryType">The geometry type</param>
/// <param name="driver">The driver</param>
/// <param name="connection">The connection string</param>
/// <param name="driverOptions">The options for the driver</param>
/// <param name="layerOptions">The options for the layer</param>
public static void CreateFromFeatureDataTable(FeatureDataTable table,
OgcGeometryType geometryType, int srid, string driver, string connection, string[] driverOptions = null, string[] layerOptions = null)
{
if (table == null)
throw new ArgumentNullException("table");
if (table.Rows.Count == 0)
throw new ArgumentException("The table contains no rows", "table");
if (geometryType < OgcGeometryType.Point || geometryType > OgcGeometryType.MultiPolygon)
throw new ArgumentException("Invalid geometry type", "geometryType");
if (string.IsNullOrWhiteSpace(driver))
throw new ArgumentException("No driver specified", "driver");
var dr = OSGeo.OGR.Ogr.GetDriverByName(driver);
if (dr == null)
throw new Exception(string.Format("Cannot load driver '{0}'!", driver));
//if (!dr.TestCapability("ODrCCreateDataSource"))
// throw new Exception(string.Format("Driver '{0}' cannot create a data source!", driver));
// Create the data source
var ds = dr.CreateDataSource(connection, driverOptions);
//if (!ds.TestCapability("ODsCCreateLayer"))
// throw new Exception(string.Format("Driver '{0}' cannot create a layer!", driver));
// Create the spatial reference
var sr = new OSGeo.OSR.SpatialReference(string.Empty);
sr.ImportFromEPSG(srid);
// Create the layer
var lyr = ds.CreateLayer(table.TableName, sr, (OgrGeometryType)geometryType, layerOptions);
sr.Dispose();
//lyr.GetSpatialRef();
foreach (System.Data.DataColumn dc in table.Columns)
{
using (var fldDef = GetFieldDefinition(dc))
lyr.CreateField(fldDef, 0);
}
using (var ld = lyr.GetLayerDefn())
{
foreach (FeatureDataRow fdr in table.Rows)
{
if ((int)fdr.Geometry.OgcGeometryType != (int)geometryType)
continue;
using (var feature = new OgrFeature(ld))
{
feature.SetGeometry(OgrGeometry.CreateFromWkb(fdr.Geometry.AsBinary()));
var idx = -1;
foreach (System.Data.DataColumn dc in table.Columns)
{
idx++;
var fd = ld.GetFieldDefn(idx);
DateTime dt;
switch (fd.GetFieldType())
{
case OgrFieldType.OFTBinary:
//Nothing
break;
case OgrFieldType.OFTDate:
dt = ((DateTime)fdr[dc]).Date;
feature.SetField(idx, dt.Year, dt.Month, dt.Day, dt.Hour, dt.Minute, dt.Second, 0);
break;
case OgrFieldType.OFTDateTime:
dt = (DateTime)fdr[dc];
feature.SetField(idx, dt.Year, dt.Month, dt.Day, dt.Hour, dt.Minute, dt.Second, 0);
break;
case OgrFieldType.OFTTime:
var tod = ((DateTime)fdr[dc]).TimeOfDay;
feature.SetField(idx, 0, 0, 0, tod.Hours, tod.Minutes, tod.Seconds, 0);
break;
case OgrFieldType.OFTInteger:
feature.SetField(idx, Convert.ToInt32(fdr[dc]));
break;
case OgrFieldType.OFTIntegerList:
var il = GetIntegerList(fdr[dc], dc.DataType);
feature.SetFieldIntegerList(idx, il.Length, il);
break;
case OgrFieldType.OFTReal:
feature.SetField(idx, Convert.ToDouble(fdr[dc]));
break;
case OgrFieldType.OFTRealList:
var dl = GetDoubleList(fdr[dc], dc.DataType);
feature.SetFieldDoubleList(idx, dl.Length, dl);
break;
case OgrFieldType.OFTString:
feature.SetField(idx, Convert.ToString(fdr[dc]));
break;
case OgrFieldType.OFTStringList:
var sl = (string[])fdr[dc];
feature.SetFieldStringList(idx, sl);
break;
}
fd.Dispose();
}
lyr.CreateFeature(feature);
feature.Dispose();
}
//ld.Dispose();
}
}
lyr.Dispose();
ds.Dispose();
dr.Dispose();
}
