public override void Close() { try { if (this.intptr_0 != IntPtr.Zero) { ShapeLib.SHPClose(this.intptr_0); this.intptr_0 = IntPtr.Zero; } } catch { } try { if (this.intptr_1 != IntPtr.Zero) { ShapeLib.DBFClose(this.intptr_1); this.intptr_1 = IntPtr.Zero; } } catch { } }
//从源文件获取字段信息 public List <R_FieldInf> GetFieldInfs(string path) { List <R_FieldInf> FieldInfs = new List <R_FieldInf>(); string dbfpath = path; #region 该部分使用FastDBF获取字段名称 返回List<string> fieldNames DbfFile dbf = new DbfFile(Encoding.Default); dbf.Open(dbfpath, FileMode.Open); DbfHeader dh = dbf.Header; List <string> fieldNames = new List <string>(); int fieldCount = dh.ColumnCount; for (int index = 0; index < fieldCount; index++) { fieldNames.Add(dh[index].Name); } dbf.Close(); #endregion #region 该部分使用Shapelib获取字段类型 返回List<string> fieldTypes //获取字段类型 IntPtr hDbf = ShapeLib.DBFOpen(dbfpath, "rb+");//"rb"(只读)"rb+"(读/写) int pointCount = ShapeLib.DBFGetRecordCount(hDbf); List <string> fieldTypes = new List <string>(); StringBuilder stringBuilder = new StringBuilder(20); int pnWidth = 10; int pnDecimals = 10; for (int index = 0; index < fieldCount; index++) { string type = TypeConvert(ShapeLib.DBFGetFieldInfo(hDbf, index, stringBuilder, ref pnWidth, ref pnDecimals).ToString()); fieldTypes.Add(type); } ShapeLib.DBFClose(hDbf); #endregion //实例化类型 for (int index = 0; index < fieldCount; index++) { FieldInfs.Add(new R_FieldInf(fieldNames[index], fieldTypes[index])); } return(FieldInfs); }
public static void ConvertToOra(string[] args) { string pkValue = string.Empty; if (args == null || args.Length < 4) { Console.Write("\nUsage: shape2ora <username/password>@dbalias> <spatial_table_name> <shape_col> <shapefile> <srid> \n\nPress any key to exit"); Console.ReadLine(); return; } string connectionstring = Utils.ParseConnectionString(args[0].ToString()); string spatial_table = args[1].ToString(); string shape_col = args[2].ToString(); string inShapeFile = args[3].ToString(); string strSRID = "NULL"; if (args.Length == 5) { strSRID = args[4].ToString(); } if (inShapeFile.ToUpper().EndsWith(".SHP")) { inShapeFile = inShapeFile.Substring(0, inShapeFile.Length - 4); } IntPtr hShp = ShapeLib.SHPOpen(inShapeFile, "rb"); IntPtr hDbf = ShapeLib.DBFOpen(inShapeFile, "rb"); if (hDbf.Equals(IntPtr.Zero)) { Console.WriteLine("\nCould not open {0}.dbf\nProbable cause: You do not have permissions or filename is incorrect\n\nPress any key to exit", inShapeFile); Console.ReadLine(); return; } OracleConnection oracon = new OracleConnection(connectionstring); oracon.Open(); try { //Check if table exists bool tabExists = TabExists(oracon, spatial_table); //Get dbf info int recCount = ShapeLib.DBFGetRecordCount(hDbf); // get shape info double[] minXY = new double[4]; double[] maxXY = new double[4]; int nEntities = 0; ShapeLib.ShapeType shapeType = 0; ShapeLib.SHPGetInfo(hShp, ref nEntities, ref shapeType, minXY, maxXY); string sqlCreateTab = @"CREATE TABLE " + spatial_table.ToUpper() + "("; bool flag = true; for (int i = 0; i < nEntities; i++) { if (flag) { sqlCreateTab = sqlCreateTab + GetShapeColumnNames(ref hDbf); sqlCreateTab = sqlCreateTab + shape_col + "\tMDSYS.SDO_GEOMETRY\n)"; if (!tabExists) { //EXECUTE table creation sql ExecuteStatement(oracon, sqlCreateTab); } flag = false; } #region Deal With Geometry IntPtr pshpObj = ShapeLib.SHPReadObject(hShp, i); // Get the SHPObject associated with our IntPtr pshpObj // We create a new SHPObject in managed code, then use Marshal.PtrToStructure // to copy the unmanaged memory pointed to by pshpObj into our managed copy. ShapeLib.SHPObject shpObj = new ShapeLib.SHPObject(); Marshal.PtrToStructure(pshpObj, shpObj); //Number of parts int nParts = shpObj.nParts; //Part starts List int[] partStarts = new int[nParts]; //Coordinate arrays double[] xCoord = new double[shpObj.nVertices]; double[] yCoord = new double[shpObj.nVertices]; double[] zCoord = new double[shpObj.nVertices]; double[] MCoord = new double[shpObj.nVertices]; // Use Marshal.Copy to copy the memory pointed // to by shpObj.padfX and shpObj.padfX (each an IntPtr) to an actual array. Marshal.Copy(shpObj.padfX, xCoord, 0, shpObj.nVertices); Marshal.Copy(shpObj.padfY, yCoord, 0, shpObj.nVertices); if (shapeType == ShapeLib.ShapeType.MultiPointM || shapeType == ShapeLib.ShapeType.PointM || shapeType == ShapeLib.ShapeType.PolygonM || shapeType == ShapeLib.ShapeType.PolyLineM) { Marshal.Copy(shpObj.padfM, zCoord, 0, shpObj.nVertices); } if (shapeType == ShapeLib.ShapeType.MultiPointZ || shapeType == ShapeLib.ShapeType.PointZ || shapeType == ShapeLib.ShapeType.PolygonZ || shapeType == ShapeLib.ShapeType.PolyLineZ) { Marshal.Copy(shpObj.padfZ, zCoord, 0, shpObj.nVertices); } if (nParts > 0) { Marshal.Copy(shpObj.paPartStart, partStarts, 0, nParts); } string sqlInsertShape = " MDSYS.SDO_GEOMETRY("; string gType = GetGTYPE(shapeType, nParts); string elem_info = "NULL"; string sdo_point = "NULL"; if (shapeType == ShapeLib.ShapeType.Point || shapeType == ShapeLib.ShapeType.PointM || shapeType == ShapeLib.ShapeType.PointZ) { sdo_point = "MDSYS.SDO_POINT_TYPE(" + xCoord[0].ToString() + "," + yCoord[0].ToString(); if (shapeType == ShapeLib.ShapeType.PointZ) { sdo_point = sdo_point + "," + zCoord[0].ToString(); } else if (shapeType == ShapeLib.ShapeType.PointM) { sdo_point = sdo_point + "," + MCoord[0].ToString(); } else { sdo_point = sdo_point + ", NULL"; } sdo_point = sdo_point + ")"; } else { elem_info = GetElemInfoString(shapeType, nParts, partStarts, shpObj.nVertices); } string vert_String = GetVerticesString(shapeType, xCoord, yCoord, zCoord, MCoord); //Construct the geometry statement sqlInsertShape = sqlInsertShape + gType + "," + strSRID + "," + sdo_point + "," + elem_info + "," + vert_String + ")"; # endregion #region Deal with Attributes string[] attrs = InsAttrSQL(ref hDbf, i); string insStatement = "INSERT INTO " + spatial_table.ToUpper() + "\n" + "(" + attrs[0] + "," + shape_col + ")\n" + " VALUES (" + attrs[1] + "," + sqlInsertShape + ")"; //Do the insert ExecuteStatement(oracon, insStatement); #endregion } //Create user_sdo_geom_metadata and spatial index for a new table if (!tabExists) { string usgm = GetUSGMString(shapeType, spatial_table, shape_col, minXY, maxXY); ExecuteStatement(oracon, usgm); string spidx = GetSPIDXString(spatial_table, shape_col); ExecuteStatement(oracon, spidx); } // free resources ShapeLib.SHPClose(hShp); ShapeLib.DBFClose(hDbf); Console.Write("Done.\nPress any key to exit"); Console.ReadLine(); }
public static void CreateShapeFile(string[] args) { string pkValue = string.Empty; bool pkValIsString = false; bool isSdoPoint = false; if (args == null || args.Length < 5) { Console.Write("\nUsage: sdo2shp <username/password>@dbalias> <spatial_table_name> <PK_col> <shape_col> <shapefile> [\"optional_where_clause\"]\n\nPress any key to exit"); Console.ReadLine(); return; } System.Reflection.Assembly assembly = System.Reflection.Assembly.GetExecutingAssembly(); string version = assembly.GetName().Version.ToString(); Console.Write("\nORA2SHP Oracle SDO to Shapefile Utility. Version: " + version + "\n"); Console.Write("====================================================================\n"); string connectionstring = Utils.ParseConnectionString(args[0].ToString()); string where_clause = string.Empty; try { where_clause = args[5].ToString().Trim('"'); } catch (Exception) { } if (where_clause != string.Empty) { if (!where_clause.ToUpper().StartsWith("WHERE")) { where_clause = " where " + where_clause; } } string outShpFile = args[4].ToString(); string tab_name = args[1].ToString(); string pk_column = args[2].ToString(); string sdo_column = args[3].ToString(); if (outShpFile.ToUpper().EndsWith(".SHP")) { outShpFile = outShpFile.Substring(0, outShpFile.Length - 4); } OracleConnection oracon = new OracleConnection(connectionstring); List <GeoInfo> lShpInfo = new List <GeoInfo>(); try { oracon.Open(); string selFirstShp = string.Empty; if (where_clause == string.Empty) { selFirstShp = "select " + args[2].ToString() + "," + args[3].ToString() + " from " + args[1] + " where rownum = 1"; } else { selFirstShp = "select " + pk_column + "," + sdo_column + " from " + tab_name + " " + where_clause + " and rownum = 1"; } try { lShpInfo = getGeomInfo(oracon, selFirstShp, sdo_column); if (lShpInfo.Count == 0) { Console.WriteLine("\nNo rows found it table {0}", args[1]); Console.ReadLine(); oracon.Close(); oracon.Dispose(); return; } } catch (Exception ex) { Console.WriteLine("\nORA2SHP ERROR:\n" + ex.Message + "\n" + ex.StackTrace); } finally { } //Check if its an SDO point geometry GeoInfo geomInfo = lShpInfo[0]; isSdoPoint = isSDOPointType(geomInfo); //Initialize shape type ShapeLib.ShapeType shpType = ShapeLib.ShapeType.Point; try { shpType = GetShapeType(geomInfo); } catch (Exception ex) { Console.Write(ex.Message); Console.ReadLine(); return; } IntPtr hShp = ShapeLib.SHPCreate(outShpFile, shpType); if (hShp.Equals(IntPtr.Zero)) { Console.WriteLine("\nCould not create {0}.shp\nProbable cause: You do not have permissions or file is in use by another process\n\nPress any key to exit", outShpFile); Console.ReadLine(); return; } string sqlselect = string.Empty; if (where_clause == string.Empty) { sqlselect = @"SELECT " + GetColumnNames(oracon, tab_name) + ", " + sdo_column + " from " + tab_name; } else { sqlselect = @"SELECT " + GetColumnNames(oracon, tab_name) + ", " + sdo_column + " from " + tab_name + " " + where_clause; } OracleCommand command1 = oracon.CreateCommand(); command1.CommandText = sqlselect; OracleDataReader dr = command1.ExecuteReader(); int iShape = 0; IntPtr hDbf = IntPtr.Zero; DataTable dt = new DataTable(); System.Collections.Hashtable ht = new System.Collections.Hashtable(); //#if (DEBUG==true) // Console.WriteLine("Before Loop Elapsed Time: " + (DateTime.Now - tstart).ToString()); tstart = DateTime.Now; //#endif #region create dbf // create dbase file hDbf = ShapeLib.DBFCreate(outShpFile); if (!hDbf.Equals(IntPtr.Zero)) { dt = dr.GetSchemaTable(); //Get table metadata dt.Columns.Add("dBaseName"); } int ordinal1 = 0; foreach (DataRow row in dt.Rows) { string name = row["ColumnName"].ToString().ToUpper(); if (name.Length > 10) //Truncate column name to 10 chars { name = name.Substring(0, 10); } int i = 0; while (ht.ContainsKey(name)) //Check if column name exists after truncation { string iVal = (i++).ToString(); if (name.Length + iVal.Length > 10) { name = name.Substring(0, 10 - iVal.Length) + iVal; } else { name = name + iVal; } } ht.Add(name, ordinal1++); row["dBaseName"] = name; string type = row["DataType"].ToString(); switch (type) { case "System.Int32": case "System.Int16": ShapeLib.DBFAddField(hDbf, name, ShapeLib.DBFFieldType.FTInteger, 16, 0); break; case "System.String": int len = Math.Min(255, int.Parse(row["ColumnSize"].ToString())); ShapeLib.DBFAddField(hDbf, name, ShapeLib.DBFFieldType.FTString, len, 0); if (name == args[1].ToString()) { pkValIsString = true; } break; case "System.Boolean": ShapeLib.DBFAddField(hDbf, name, ShapeLib.DBFFieldType.FTLogical, 5, 0); break; case "System.Double": case "System.Float": case "System.Decimal": int prec = int.Parse(row["NumericPrecision"].ToString()); int scale = int.Parse(row["NumericScale"].ToString()); ShapeLib.DBFAddField(hDbf, name, ShapeLib.DBFFieldType.FTDouble, prec, scale); break; case "System.DateTime": ShapeLib.DBFAddField(hDbf, name, ShapeLib.DBFFieldType.FTDate, 8, 0); break; default: ht.Remove(name); row["dBaseName"] = null; ordinal1--; break; } } #endregion create dbf ShapeCreator sc = new ShapeCreator(); IntPtr pShp = new IntPtr(); while (dr.Read()) { //Console.Write("Processing record ..." + iShape.ToString() +"."); pkValue = dr[args[2]].ToString(); OraShape shp = new OraShape(); ShapePartInfo shpInfo = new ShapePartInfo(); try { GeoInfo geoInfo = new GeoInfo(); geoInfo.Geo = (SdoGeometry)dr[sdo_column]; shp = sc.getOraShape(geoInfo, shpType, Convert.ToInt32(dr[pk_column])); } catch (Exception) { continue; } if (!isSdoPoint) { try { shpInfo = sc.getShapePartInfo(geomInfo); } catch (Exception) { continue; } shp.nParts = shpInfo.nParts; shp.PartType = null; shp.PartStarts = shpInfo.PartStarts; if (shp.MList == null || shp.MList.Length == 0) { if (shp.ZList == null || shp.ZList.Length == 0) { pShp = ShapeLib.SHPCreateObject(shpType, -1, shp.nParts, shp.PartStarts, null, shp.nVertices, shp.XList, shp.YList, null, null); } else { pShp = ShapeLib.SHPCreateObject(shpType, -1, shp.nParts, shp.PartStarts, null, shp.nVertices, shp.XList, shp.YList, shp.ZList, null); } } else { pShp = ShapeLib.SHPCreateObject(shpType, -1, shp.nParts, shp.PartStarts, null, shp.nVertices, shp.XList, shp.YList, null, shp.MList); } } else { shp = sc.CreatePointShape(geomInfo, shpType, pkValue); if (shpType == ShapeLib.ShapeType.PointM) { pShp = ShapeLib.SHPCreateObject(shpType, -1, 0, null, null, 1, shp.XList, shp.YList, null, shp.MList); } else if (shpType == ShapeLib.ShapeType.PointZ) { pShp = ShapeLib.SHPCreateObject(shpType, -1, 0, null, null, 1, shp.XList, shp.YList, shp.ZList, null); } else { pShp = ShapeLib.SHPCreateObject(shpType, -1, 0, null, null, 1, shp.XList, shp.YList, null, null); } } try //In case of invalid shapes { ShapeLib.SHPWriteObject(hShp, -1, pShp); } catch (Exception) { Console.WriteLine("Shape with " + args[2].ToString().ToUpper() + "= " + pkValue + " is invalid. Number of vertices: " + shp.nVertices.ToString() + " Number of parts: " + shp.nParts.ToString()); continue; } ShapeLib.SHPDestroyObject(pShp); foreach (DataRow row in dt.Rows) { if (row["dBaseName"] == null || row["dBaseName"].ToString() == string.Empty) { continue; } int ordinal = (int)ht[row["dBaseName"].ToString()]; string fieldName = row["ColumnName"].ToString(); if (dr[fieldName] is DBNull) { continue; } switch (row["DataType"].ToString()) { case "System.Int32": case "System.Int16": ShapeLib.DBFWriteIntegerAttribute(hDbf, iShape, ordinal, int.Parse(dr[fieldName].ToString())); break; case "System.String": ShapeLib.DBFWriteStringAttribute(hDbf, iShape, ordinal, dr[fieldName].ToString()); break; case "System.Boolean": ShapeLib.DBFWriteLogicalAttribute(hDbf, iShape, ordinal, bool.Parse(dr[fieldName].ToString())); break; case "System.Double": case "System.Float": case "System.Decimal": ShapeLib.DBFWriteDoubleAttribute(hDbf, iShape, ordinal, double.Parse(dr[fieldName].ToString())); break; case "System.DateTime": DateTime date = DateTime.Parse(dr[fieldName].ToString()); ShapeLib.DBFWriteDateAttribute(hDbf, iShape, ordinal, date); break; } } iShape++; //dr.NextResult(); } Console.WriteLine("Converted " + iShape.ToString() + " shapes in: " + (DateTime.Now - tstart).ToString()); // free resources ShapeLib.SHPClose(hShp); ShapeLib.DBFClose(hDbf); //Create projection file if needed CreatePrjFile(oracon, args[1].ToString(), args[3].ToString(), outShpFile); //Console.Write("Done.\n"); Console.Write("\nCreated shapefile " + args[4].ToString() + " with " + iShape.ToString() + " records.\nTotal Elapsed Time: " + (DateTime.Now - tProgramstart).ToString() + "\nPress any key to exit."); Console.ReadLine(); } catch (Exception ex) { Console.Write("ERROR: " + ex.Message + "\nPress any key to exit"); Utils.WriteErrLog(ex); Console.ReadLine(); } finally { //Close and free connection oracon.Close(); oracon.Dispose(); } }
/// <summary> /// Creates the database file that contains meta data about the shapes. /// </summary> /// <param name="shapetype"></param> /// <param name="SHPFile"></param> private void WriteDBF(ShapeLib.ShapeType shapetype, string SHPFile, bool convertTags) { string filename; int shapes = 0; List <MetaData> elementData; List <string> fields = new List <string>(); switch (shapetype) { case ShapeLib.ShapeType.Polygon: filename = SHPFile + "-polygons"; shapes = areas; elementData = polygonsData; break; case ShapeLib.ShapeType.PolyLine: filename = SHPFile + "-lines"; shapes = ways; elementData = linesData; break; case ShapeLib.ShapeType.Point: filename = SHPFile + "-points"; shapes = points; elementData = pointsData; break; default: return; } IntPtr hDbf = ShapeLib.DBFCreate(filename); if (hDbf.Equals(IntPtr.Zero)) { Console.WriteLine("Error: Unable to create {0}.dbf!", filename); return; } // add some fields. Fields have to be initialized before data is added int iRet = ShapeLib.DBFAddField(hDbf, "shapeID", ShapeLib.DBFFieldType.FTInteger, 4, 0); iRet = ShapeLib.DBFAddField(hDbf, "shapeName", ShapeLib.DBFFieldType.FTString, 50, 0); // if tags have to be converted, collect the names of all tags if (convertTags) { for (int iShape = 0; iShape < shapes; iShape++) { foreach (KeyValuePair <string, string> entry in elementData[iShape].Tags) { // Make sure that fields are only added once if (!fields.Contains(entry.Key)) { fields.Add(entry.Key); } } } // create fields in the DB for all tags foreach (string field in fields) { iRet = ShapeLib.DBFAddField(hDbf, field, ShapeLib.DBFFieldType.FTString, 128, 0); } } // populate for (int iShape = 0; iShape < shapes; iShape++) { iRet = (ShapeLib.DBFWriteIntegerAttribute(hDbf, iShape, 0, iShape)); iRet = (ShapeLib.DBFWriteStringAttribute(hDbf, iShape, 1, elementData[iShape].Name)); // If tags should be converted, write their values to the appropriate fields if (convertTags) { foreach (KeyValuePair <string, string> entry in elementData[iShape].Tags) { // Cut the entry after 128 characters int stringLength = entry.Value.Length; if (stringLength > 127) { stringLength = 127; } iRet = ShapeLib.DBFWriteStringAttribute(hDbf, iShape, fields.IndexOf(entry.Key) + 2, entry.Value.Substring(0, stringLength)); } } } ShapeLib.DBFClose(hDbf); }
public List <InputFields> GetFieldValues(List <R_FieldInf> fieldInfs, FieldMapping vitalFields, string path) { List <InputFields> result = new List <InputFields>(); string dbfpath = path; IntPtr hDbf = ShapeLib.DBFOpen(dbfpath, "rb+"); int FieldCount = fieldInfs.Count(); int PointCount = ShapeLib.DBFGetRecordCount(hDbf); for (int pi = 0; pi < PointCount; pi++) { InputFields point = new InputFields(); string other = "{"; for (int i = 0; i < FieldCount; i++) { if (fieldInfs[i].FieldName == vitalFields.ID) { IntPtr FieldvaluePtr = ShapeLib.DBFReadStringAttribute(hDbf, pi, i); string FieldValue = Marshal.PtrToStringAnsi(FieldvaluePtr, 255).Replace("\0", ""); point.ID = FieldValue; } else { if (fieldInfs[i].FieldName == vitalFields.Name) { IntPtr FieldvaluePtr = ShapeLib.DBFReadStringAttribute(hDbf, pi, i); string FieldValue = Marshal.PtrToStringAnsi(FieldvaluePtr, 255).Replace("\0", ""); point.Name = FieldValue; } else { if (fieldInfs[i].FieldName == vitalFields.Address) { IntPtr FieldvaluePtr = ShapeLib.DBFReadStringAttribute(hDbf, pi, i); string FieldValue = Marshal.PtrToStringAnsi(FieldvaluePtr, 255).Replace("\0", ""); point.Addresss = FieldValue; } else { if (fieldInfs[i].FieldName == vitalFields.Xcoordinate) { IntPtr FieldvaluePtr = ShapeLib.DBFReadStringAttribute(hDbf, pi, i); string FieldValue = Marshal.PtrToStringAnsi(FieldvaluePtr, 255).Replace("\0", ""); point.X = FieldValue; } else { if (fieldInfs[i].FieldName == vitalFields.Ycoordinate) { IntPtr FieldvaluePtr = ShapeLib.DBFReadStringAttribute(hDbf, pi, i); string FieldValue = Marshal.PtrToStringAnsi(FieldvaluePtr, 255).Replace("\0", ""); point.Y = FieldValue; } else { if (fieldInfs[i].FieldType == "text") { IntPtr FieldvaluePtr = ShapeLib.DBFReadStringAttribute(hDbf, pi, i); string FieldValue = Marshal.PtrToStringAnsi(FieldvaluePtr, 255).Replace("\0", ""); other += "\"" + fieldInfs[i].FieldName + "\":\"" + FieldValue + "\","; } else { if (fieldInfs[i].FieldType == "integer") { int FieldValue = ShapeLib.DBFReadIntegerAttribute(hDbf, pi, i); other += "\"" + fieldInfs[i].FieldName + "\":" + FieldValue + ","; } else { double FieldValue = ShapeLib.DBFReadDoubleAttribute(hDbf, pi, i); other += "\"" + fieldInfs[i].FieldName + "\":" + FieldValue + ","; } } } } } } } } other = other.Substring(0, other.Length - 1) + "}"; point.Other = other; result.Add(point); } ShapeLib.DBFClose(hDbf); return(result); }