/// <summary>
/// Exports the selected database to fieldAttributes filegeodatabase
/// </summary>
/// <param name="progressBar">A progress bar object</param>
/// <param name="logFunction">A log function</param>
/// <param name="mdbFile">The outputShapefileName of an addressing database</param>
/// <param name="approvedOnly">1 to include only approved, 0 to include all</param>
/// <returns>Feature set</returns>
public static FeatureSet GetRoadFeatureSetFromAdmAdrMdb(ref ToolStripProgressBar progressBar, Action <string, bool> logFunction, string mdbFile, int approvedOnly = 0)
{
// Setup SQL string to select all or only some streets
string sqlStatement;
if (approvedOnly == 1)
{
sqlStatement = "SELECT ADRROADID, NAMEENGLISH, NAMEARABIC, NAMEPOPULARENGLISH, NAMEPOPULARARABIC, ROADTYPE, ADRDISTRICTID, APPROVED FROM ADM_ADRROAD WHERE APPROVED = 1";
}
else
{
sqlStatement = "SELECT ADRROADID, NAMEENGLISH, NAMEARABIC, NAMEPOPULARENGLISH, NAMEPOPULARARABIC, ROADTYPE, ADRDISTRICTID, APPROVED FROM ADM_ADRROAD";
}
logFunction(sqlStatement, true);
// Setup fieldAttributes dictionary to hold road names
var roadNames = new Dictionary <int, DataRow>();
// Connect to database and load names
using (var database = new Database(mdbFile))
{
var dataTable = database.Query(sqlStatement);
progressBar.Value = 0;
foreach (DataRow mRow in dataTable.Rows)
{
int mRoadId;
if (int.TryParse(mRow["ADRROADID"].ToString(), out mRoadId))
{
roadNames.Add(mRoadId, mRow);
}
}
// Connect to database using OGR to load geometries
var ogrDriver = Ogr.GetDriverByName("PGeo");
if (ogrDriver == null)
{
return(null);
}
logFunction("Loaded driver", true);
var ogrDataSource = ogrDriver.Open(mdbFile, 0);
if (ogrDataSource == null)
{
return(null);
}
logFunction("Loaded datasource", true);
var roadSegmentLayer = ogrDataSource.GetLayerByName("ADRROADSEGMENT");
if (roadSegmentLayer == null)
{
logFunction("Could not load layer", true);
return(null);
}
logFunction("Loaded layer", true);
// Create simplified roads
var simplifiedRoadsFeatureSet = new SimplifiedRoads();
int numberOfRoads = roadNames.Count();
progressBar.Value = 0;
var roadCounter = 0;
var segmentCounter = 0;
foreach (var roadName in roadNames)
{
roadCounter++;
int roadIdentifier = roadName.Key;
DataRow roadData = roadName.Value;
roadSegmentLayer.ResetReading();
roadSegmentLayer.SetAttributeFilter("ADRROADID=" + roadIdentifier);
OSGeo.OGR.Feature roadSegment;
var mWkbReader = new WkbReader();
while (null != (roadSegment = roadSegmentLayer.GetNextFeature()))
{
segmentCounter++;
OSGeo.OGR.Geometry mGeometry = roadSegment.GetGeometryRef();
if (mGeometry != null)
{
mGeometry.FlattenTo2D();
byte[] mWkb = new Byte[mGeometry.WkbSize()];
mGeometry.ExportToWkb(mWkb);
IMultiLineString multiLine;
IGeometry roadGeometry = mWkbReader.Read(mWkb);
if (roadGeometry.GetType() == typeof(DotSpatial.Topology.LineString))
{
var mLineStrings = new List <LineString>();
mLineStrings.Add(roadGeometry as LineString);
multiLine = new DotSpatial.Topology.MultiLineString(mLineStrings);
}
else
{
multiLine = roadGeometry as IMultiLineString;
}
int districtIdentifier, roadClass, roadApproved;
int.TryParse(roadData["ADRDISTRICTID"].ToString(), out districtIdentifier);
var nameArabic = roadData["NAMEARABIC"].ToString();
var nameLatin = roadData["NAMEENGLISH"].ToString();
var namePopularArabic = roadData["NAMEPOPULARARABIC"].ToString();
var namePopularLatin = roadData["NAMEPOPULARENGLISH"].ToString();
int.TryParse(roadData["ROADTYPE"].ToString(), out roadClass);
int.TryParse(roadData["APPROVED"].ToString(), out roadApproved);
simplifiedRoadsFeatureSet.AddNewRow(
roadGeometry,
districtIdentifier,
roadIdentifier,
nameArabic,
nameLatin,
namePopularArabic,
namePopularLatin,
roadClass,
roadApproved);
}
else
{
logFunction("No geometry", true);
}
if (roadCounter % 100 == 0 || roadCounter == numberOfRoads)
{
progressBar.ProgressBar.Value = (int)Math.Round((double)(roadCounter * 100 / numberOfRoads));
Application.DoEvents();
}
}
}
simplifiedRoadsFeatureSet.UpdateExtent();
logFunction("Done creating featureset", true);
return(simplifiedRoadsFeatureSet);
}
}
public static void AnalyzeRoadBoundingBoxes(string pMdbFile, string pLogFile, Action <Object, bool> aLog, int pOnlyApproved = 0)
{
var roadSuspects = new List <RoadSuspect>();
// Setup SQL string to select all or only some streets
string sql;
if (pOnlyApproved == 1)
{
sql = "SELECT ADRROADID, NAMEENGLISH, NAMEARABIC, NAMEPOPULARENGLISH, NAMEPOPULARARABIC, ROADTYPE, ADRDISTRICTID, APPROVED FROM ADM_ADRROAD WHERE APPROVED = 1 ORDER BY ADRROADID";
}
else
{
sql = "SELECT ADRROADID, NAMEENGLISH, NAMEARABIC, NAMEPOPULARENGLISH, NAMEPOPULARARABIC, ROADTYPE, ADRDISTRICTID, APPROVED FROM ADM_ADRROAD ORDER BY ADRROADID";
}
// Setup fieldAttributes dictionary to hold road names
var mRoads = new Dictionary <int, DataRow>();
// Connect to database and load names
var mOdbcConn = new OdbcConnection("Driver={Microsoft Access Driver (*.mdb, *.accdb)};Dbq=" + pMdbFile + ";Uid=Admin;Pwd=;");
var mDataAdapter = new OdbcDataAdapter(sql, mOdbcConn);
var mDataTable = new DataTable();
mDataAdapter.Fill(mDataTable);
int ctr1 = 0;
int total1 = mDataTable.Rows.Count;
foreach (DataRow mRow in mDataTable.Rows)
{
int mRoadId;
if (int.TryParse(mRow["ADRROADID"].ToString(), out mRoadId))
{
mRoads.Add(mRoadId, mRow);
}
ctr1++;
if (ctr1 % 100 == 0)
{
double mFraction1 = (((double)ctr1 / (double)total1) * 100);
Application.DoEvents();
}
}
Debug.WriteLine("Done loading road names");
mOdbcConn.Close();
mDataAdapter = null;
mOdbcConn = null;
// Connect to database using OGR to load geometries
var mDriver = Ogr.GetDriverByName("PGeo");
if (mDriver == null)
{
return;
}
Debug.WriteLine("Loaded driver");
var mSource = mDriver.Open(pMdbFile, 0);
if (mSource == null)
{
return;
}
Debug.WriteLine("Loaded datasource");
var mLayer = mSource.GetLayerByName("ADRROADSEGMENT");
if (mLayer == null)
{
return;
}
Debug.WriteLine("Loaded layer");
// Create simplified roads
var mSimplifiedRoadsFeatureSet = new SimplifiedRoads();
int mTotal2 = mLayer.GetFeatureCount(0);
foreach (var mRoadID in mRoads)
{
mLayer.SetAttributeFilter("ADRROADID = " + mRoadID.Key);
OSGeo.OGR.Feature mFeature;
int mNumShapes = 0;
double mLength = 0;
double mLeft = 0;
double mBottom = 0;
double mTop = 0;
double mRight = 0;
int mDuplicates = 0;
List <string> mGeoms = new List <string>();
while (null != (mFeature = mLayer.GetNextFeature()))
{
OSGeo.OGR.Geometry mGeom = mFeature.GetGeometryRef();
if (mGeom != null)
{
string wkt;
mGeom.ExportToWkt(out wkt);
OSGeo.OGR.Envelope mExtent = new OSGeo.OGR.Envelope();
mGeom.GetEnvelope(mExtent);
if (mTop == 0 || mExtent.MaxY > mTop)
{
mTop = mExtent.MaxY;
}
if (mBottom == 0 || mExtent.MinY < mBottom)
{
mBottom = mExtent.MinY;
}
if (mLeft == 0 || mExtent.MinX < mLeft)
{
mLeft = mExtent.MinX;
}
if (mRight == 0 || mExtent.MaxX > mRight)
{
mRight = mExtent.MaxX;
}
if (!mGeoms.Contains(wkt))
{
mLength = mLength + mGeom.Length();
mGeoms.Add(wkt);
}
else
{
mDuplicates++;
}
}
mNumShapes++;
}
var mBoxWidth = mRight - mLeft;
var mBoxHeight = mTop - mBottom;
var mBoxDiagonal = Math.Sqrt(mBoxWidth * mBoxWidth + mBoxHeight * mBoxHeight);
roadSuspects.Add(new RoadSuspect(mRoadID.Key,
mNumShapes,
mDuplicates,
Math.Round(mBoxDiagonal, 1),
Math.Round(mLength, 1)));
if (mLength > (2.2 * mBoxDiagonal) || mDuplicates > 0)
{
aLog(String.Format("r: {0}, s: {1}, l: {2}, d: {3}, bb: {4}",
mRoadID.Key,
mNumShapes,
Math.Round(mLength, 1),
mDuplicates,
Math.Round(mBoxDiagonal, 1)),
true);
}
}
using (var mCsvWriter = new CsvHelper.CsvWriter(new ExcelSerializer(pLogFile)))
{
mCsvWriter.WriteRecords(roadSuspects);
aLog("Wrote output to " + pLogFile, true);
aLog("Operation completed...", true);
}
}
