private SdoGeometry WritePoint(IGeometry geometry)
{
int dimension = GetGeometryDimension(geometry);
double[] coord = ConvertCoordinates(geometry.Coordinates, dimension);
SdoGeometry sdoGeometry = new SdoGeometry();
sdoGeometry.GeometryType = (int)SdoGeometryTypes.GTYPE.POINT;
sdoGeometry.Dimensionality = dimension;
sdoGeometry.LRS = 0;
sdoGeometry.Sdo_Srid = geometry.SRID;
sdoGeometry.ElemArrayOfInts = new [] { 1, (int)SdoGeometryTypes.ETYPE_SIMPLE.POINT, 1 };
sdoGeometry.OrdinatesArrayOfDoubles = coord;
sdoGeometry.PropertiesToGTYPE();
return sdoGeometry;
}
public IGeometry Read(SdoGeometry geometry)
{
geometry.PropertiesToGTYPE();
return this.ReadGeometry(geometry);
}
public BIMRLProcessModel(IModel model, bool update)
{
//IfcProject proj = model.IfcProject;
IfcStore modelStore = model as IfcStore;
string currStep = String.Empty;
_bimrlCommon.resetAll();
// Connect to Oracle DB
DBOperation.refBIMRLCommon = _bimrlCommon; // important to ensure DBoperation has reference to this object!!
try
{
DBOperation.ExistingOrDefaultConnection();
}
catch
{
if (DBOperation.UIMode)
{
BIMRLErrorDialog erroDlg = new BIMRLErrorDialog(_bimrlCommon);
erroDlg.ShowDialog();
}
else
{
Console.Write(_bimrlCommon.ErrorMessages);
}
return;
}
DBOperation.commitInterval = 5000;
// Initial Spatial index for later use
BIMRLSpatialIndex spIdx = new BIMRLSpatialIndex(_bimrlCommon);
try
{
DBOperation.beginTransaction();
IIfcProject firstProject;
if (modelStore.IsFederation)
{
IfcStore firstModel = modelStore.ReferencedModels.FirstOrDefault().Model as IfcStore;
firstProject = firstModel.Instances.OfType <IIfcProject>().FirstOrDefault();
}
else
{
firstProject = modelStore.Instances.OfType <IIfcProject>().FirstOrDefault();
}
string projLName;
// Check whether Model has been defined before
if (string.IsNullOrEmpty(firstProject.LongName))
{
projLName = firstProject.Name + " - Federated";
}
else
{
projLName = firstProject.LongName;
}
string modelNameFromFile;
if (!string.IsNullOrEmpty(modelStore.FileName))
{
modelNameFromFile = Path.GetFileNameWithoutExtension(modelStore.FileName);
}
else
{
modelNameFromFile = firstProject.Name + " - " + firstProject.LongName;
}
currStep = "Getting Federated ID from BIMRL_FEDERATEDMODEL - Model name, Project name and longname: " + modelNameFromFile + "; " + firstProject.Name + "; " + firstProject.LongName;
FederatedModelInfo fedModel;
FedIDStatus stat = DBOperation.getFederatedModel(modelNameFromFile, projLName, firstProject.Name, out fedModel);
if (stat == FedIDStatus.FedIDNew)
{
DBOperation.currFedModel = fedModel;
// Create new set of tables using the fedID as suffix
int retStat = DBOperation.createModelTables(DBOperation.currFedModel.FederatedID);
}
else
{
DBOperation.currFedModel = fedModel;
if (!fedModel.Owner.Equals(DBOperation.DBUserID))
{
_bimrlCommon.StackPushError("%Error: Only the Owner (" + fedModel.Owner + ") can delete or override existing model (!" + fedModel.ModelName + ")");
throw new Exception("%Error: Unable to overwrite exisitng model");
}
// Drop and recreate tables
currStep = "Dropping existing model tables (ID: " + DBOperation.currFedModel.FederatedID.ToString("X4") + ")";
int retStat = DBOperation.dropModelTables(DBOperation.currFedModel.FederatedID);
currStep = "Creating model tables (ID: " + DBOperation.currFedModel.FederatedID.ToString("X4") + ")";
retStat = DBOperation.createModelTables(DBOperation.currFedModel.FederatedID);
}
DBOperation.currSelFedID = DBOperation.currFedModel.FederatedID; // set the static variable keeping the selected Fed Id
if (modelStore.IsFederation)
{
// get all models
foreach (IReferencedModel refModel in modelStore.ReferencedModels)
{
IfcStore m = refModel.Model as IfcStore;
currStep = "Getting Model ID for Federated model ID:" + DBOperation.currFedModel.FederatedID.ToString("X4");
_bimrlCommon.ClearDicts();
_ModelID = DBOperation.getModelID(DBOperation.currFedModel.FederatedID);
doModel(m);
BIMRLUtils.ResetIfcUnitDicts();
}
}
else
{
currStep = "Getting Model ID for Federated model ID:" + DBOperation.currFedModel.FederatedID.ToString("X4");
_ModelID = DBOperation.getModelID(DBOperation.currFedModel.FederatedID);
doModel(modelStore);
BIMRLUtils.ResetIfcUnitDicts();
}
}
catch (Exception e)
{
string excStr = "%%Error - " + e.Message + "\n\t" + currStep;
_bimrlCommon.StackPushError(excStr);
DBOperation.endTransaction(false); // rollback
}
DBOperation.endTransaction(true); // Commit = true
try
{
DBOperation.beginTransaction();
OracleCommand cmd = new OracleCommand("", DBOperation.DBConn);
// Define the spatial index metadata
double marginX = (_bimrlCommon.URT_X - _bimrlCommon.LLB_X) * 0.2; // 20% margin
double marginY = (_bimrlCommon.URT_Y - _bimrlCommon.LLB_Y) * 0.2; // 20% margin
double marginZ = (_bimrlCommon.URT_Z - _bimrlCommon.LLB_Z) * 0.2; // 20% margin
double lowerX = _bimrlCommon.LLB_X - marginX;
double upperX = _bimrlCommon.URT_X + marginX;
double lowerY = _bimrlCommon.LLB_Y - marginY;
double upperY = _bimrlCommon.URT_Y + marginY;
double lowerZ = _bimrlCommon.LLB_Z - marginZ;
double upperZ = _bimrlCommon.URT_Z + marginZ;
string sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_ELEMENT_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','GEOMETRYBODY',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', " + lowerX.ToString() + ", " + upperX.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Y', " + lowerY.ToString() + ", " + upperY.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Z', " + lowerZ.ToString() + ", " + upperZ.ToString() + ", 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_ELEMENT_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','GEOMETRYBODY_BBOX',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', " + lowerX.ToString() + ", " + upperX.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Y', " + lowerY.ToString() + ", " + upperY.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Z', " + lowerZ.ToString() + ", " + upperZ.ToString() + ", 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_ELEMENT_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','GEOMETRYBODY_BBOX_CENTROID',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', " + lowerX.ToString() + ", " + upperX.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Y', " + lowerY.ToString() + ", " + upperY.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Z', " + lowerZ.ToString() + ", " + upperZ.ToString() + ", 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_ELEMENT_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','GEOMETRYFOOTPRINT',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', " + lowerX.ToString() + ", " + upperX.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Y', " + lowerY.ToString() + ", " + upperY.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Z', " + lowerZ.ToString() + ", " + upperZ.ToString() + ", 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_ELEMENT_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','GEOMETRYAXIS',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', " + lowerX.ToString() + ", " + upperX.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Y', " + lowerY.ToString() + ", " + upperY.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Z', " + lowerZ.ToString() + ", " + upperZ.ToString() + ", 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_ELEMENT_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','BODY_MAJOR_AXIS1',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', -1.01, 1.01, 0.000001),"
+ "SDO_DIM_ELEMENT('Y', -1.01, 1.01, 0.000001),"
+ "SDO_DIM_ELEMENT('Z', -1.01, 1.01, 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_ELEMENT_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','BODY_MAJOR_AXIS2',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', -1.01, 1.01, 0.000001),"
+ "SDO_DIM_ELEMENT('Y', -1.01, 1.01, 0.000001),"
+ "SDO_DIM_ELEMENT('Z', -1.01, 1.01, 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_ELEMENT_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','BODY_MAJOR_AXIS3',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', -1.01, 1.01, 0.000001),"
+ "SDO_DIM_ELEMENT('Y', -1.01, 1.01, 0.000001),"
+ "SDO_DIM_ELEMENT('Z', -1.01, 1.01, 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_ELEMENT_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','BODY_MAJOR_AXIS_CENTROID',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', " + lowerX.ToString() + ", " + upperX.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Y', " + lowerY.ToString() + ", " + upperY.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Z', " + lowerZ.ToString() + ", " + upperZ.ToString() + ", 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_TOPO_FACE_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','POLYGON',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', " + lowerX.ToString() + ", " + upperX.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Y', " + lowerY.ToString() + ", " + upperY.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Z', " + lowerZ.ToString() + ", " + upperZ.ToString() + ", 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_TOPO_FACE_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','NORMAL',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', -1.01, 1.01, 0.000001),"
+ "SDO_DIM_ELEMENT('Y', -1.01, 1.01, 0.000001),"
+ "SDO_DIM_ELEMENT('Z', -1.01, 1.01, 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "insert into USER_SDO_GEOM_METADATA (TABLE_NAME, COLUMN_NAME, DIMINFO, SRID) VALUES "
+ "('BIMRL_TOPO_FACE_" + DBOperation.currFedModel.FederatedID.ToString("X4") + "','CENTROID',"
+ "SDO_DIM_ARRAY(SDO_DIM_ELEMENT('X', " + lowerX.ToString() + ", " + upperX.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Y', " + lowerY.ToString() + ", " + upperY.ToString() + ", 0.000001),"
+ "SDO_DIM_ELEMENT('Z', " + lowerZ.ToString() + ", " + upperZ.ToString() + ", 0.000001)),"
+ "NULL)";
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
sqlStmt = "Update BIMRL_FEDERATEDMODEL SET LastUpdateDate=sysdate WHERE FederatedID=" + DBOperation.currFedModel.FederatedID.ToString();
currStep = sqlStmt;
cmd.CommandText = sqlStmt;
cmd.ExecuteNonQuery();
// get the world BBox coordinates and update the BIMRL_FEDERATEDMODEL Table
SdoGeometry Bbox = new SdoGeometry();
Bbox.Dimensionality = 3;
Bbox.LRS = 0;
Bbox.GeometryType = (int)SdoGeometryTypes.GTYPE.POINT;
int gType = Bbox.PropertiesToGTYPE();
int[] elemInfoArr = { 1, (int)SdoGeometryTypes.ETYPE_SIMPLE.POLYGON_EXTERIOR, 1 };
Bbox.ElemArrayOfInts = elemInfoArr;
double[] coordArr = new double[6];
coordArr[0] = _bimrlCommon.LLB_X;
coordArr[1] = _bimrlCommon.LLB_Y;
coordArr[2] = _bimrlCommon.LLB_Z;
coordArr[3] = _bimrlCommon.URT_X;
coordArr[4] = _bimrlCommon.URT_Y;
coordArr[5] = _bimrlCommon.URT_Z;
Bbox.OrdinatesArrayOfDoubles = coordArr;
// Create spatial index from the new model (list of triangles are accummulated during processing of geometries)
sqlStmt = "update BIMRL_FEDERATEDMODEL SET WORLDBBOX=:1 , MAXOCTREELEVEL=:2 WHERE FEDERATEDID=" + DBOperation.currFedModel.FederatedID.ToString();
currStep = sqlStmt;
OracleCommand command = new OracleCommand(" ", DBOperation.DBConn);
command.CommandText = sqlStmt;
OracleParameter[] sdoGeom2 = new OracleParameter[3];
sdoGeom2[0] = command.Parameters.Add("1", OracleDbType.Object);
sdoGeom2[0].Direction = ParameterDirection.Input;
sdoGeom2[0].UdtTypeName = "MDSYS.SDO_GEOMETRY";
sdoGeom2[0].Value = Bbox;
sdoGeom2[0].Size = 1;
sdoGeom2[1] = command.Parameters.Add("2", OracleDbType.Int16);
sdoGeom2[1].Direction = ParameterDirection.Input;
sdoGeom2[1].Value = DBOperation.OctreeSubdivLevel;
sdoGeom2[1].Size = 1;
int commandStatus = command.ExecuteNonQuery();
if (DBOperation.OnepushETL)
{
DBOperation.commitInterval = 10000;
int octreeLevel = DBOperation.computeRecomOctreeLevel(DBOperation.currFedModel.FederatedID);
// 1. Create Octree spatial indexes and the Brep Topology Faces
spIdx.createSpatialIndexFromBIMRLElement(DBOperation.currFedModel.FederatedID, null, true, true);
// 2. Update major Axes and OBB
BIMRLUtils.updateMajorAxesAndOBB(DBOperation.currFedModel.FederatedID, null);
// 3. Enhance Space Boundary
EnhanceBRep eBrep = new EnhanceBRep();
eBrep.enhanceSpaceBoundary(null);
// 4. Process Face orientations. We will procees the normal face first and then after that the spacial ones (OBB, PROJOBB)
string whereCond2 = "";
BIMRLCommon.appendToString(" TYPE NOT IN ('OBB','PROJOBB')", " AND ", ref whereCond2);
eBrep.ProcessOrientation(whereCond2);
whereCond2 = "";
BIMRLCommon.appendToString(" TYPE='OBB'", " AND ", ref whereCond2);
eBrep.ProcessOrientation(whereCond2);
whereCond2 = "";
BIMRLCommon.appendToString(" TYPE='PROJOBB'", " AND ", ref whereCond2);
eBrep.ProcessOrientation(whereCond2);
// 5. Create Graph Data
BIMRLGraph.GraphData graphData = new BIMRLGraph.GraphData();
graphData.createCirculationGraph(DBOperation.currFedModel.FederatedID);
graphData.createSpaceAdjacencyGraph(DBOperation.currFedModel.FederatedID);
sqlStmt = "UPDATE BIMRL_FEDERATEDMODEL SET LASTUPDATEDATE=sysdate";
BIMRLCommon.appendToString("MAXOCTREELEVEL=" + octreeLevel.ToString(), ", ", ref sqlStmt);
BIMRLCommon.appendToString("WHERE FEDERATEDID=" + DBOperation.currFedModel.FederatedID.ToString(), " ", ref sqlStmt);
DBOperation.executeSingleStmt(sqlStmt);
}
else
{
// Minimum (without One-push ETL, create the bounding boxes (AABB)
spIdx.createSpatialIndexFromBIMRLElement(DBOperation.currFedModel.FederatedID, null, false, false);
}
var location = new Uri(Assembly.GetEntryAssembly().GetName().CodeBase);
string exePath = new FileInfo(location.AbsolutePath).Directory.FullName;
// (Re)-Create the spatial indexes
DBOperation.executeScript(Path.Combine(exePath, "script", "BIMRL_Idx_SpatialIndexes.sql"), DBOperation.currFedModel.FederatedID);
DBOperation.executeScript(Path.Combine(exePath, "script", "BIMRL_Idx_TopoFace.sql"), DBOperation.currFedModel.FederatedID);
DBOperation.executeScript(Path.Combine(exePath, "script", "BIMRL_Idx_MajorAxes.sql"), DBOperation.currFedModel.FederatedID);
//sqlStmt = "Create Index IDX_BIMRLELEM_GEOM_" + currFedID.ToString("X4") + " on BIMRL_ELEMENT_" + currFedID.ToString("X4")
// + " (GEOMETRYBODY) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_BIMRLELEM_GEOMBB_" + currFedID.ToString("X4") + " on BIMRL_ELEMENT_" + currFedID.ToString("X4")
// + " (GeometryBody_BBOX) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_BIMRLELEM_GEOMBBC_" + currFedID.ToString("X4") + " on BIMRL_ELEMENT_" + currFedID.ToString("X4")
// + " (GeometryBody_BBOX_CENTROID) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_BIMRLELEM_GEOMFP_" + currFedID.ToString("X4") + " on BIMRL_ELEMENT_" + currFedID.ToString("X4")
// + " (GeometryFootprint) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_BIMRLELEM_GEOMAX_" + currFedID.ToString("X4") + " on BIMRL_ELEMENT_" + currFedID.ToString("X4")
// + " (GeometryAxis) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_BIMRLELEM_GEOMMJ1_" + currFedID.ToString("X4") + " on BIMRL_ELEMENT_" + currFedID.ToString("X4")
// + " (BODY_MAJOR_AXIS1) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_BIMRLELEM_GEOMMJ2_" + currFedID.ToString("X4") + " on BIMRL_ELEMENT_" + currFedID.ToString("X4")
// + " (BODY_MAJOR_AXIS2) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_BIMRLELEM_GEOMMJ3_" + currFedID.ToString("X4") + " on BIMRL_ELEMENT_" + currFedID.ToString("X4")
// + " (BODY_MAJOR_AXIS3) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_BIMRLELEM_GEOMMJC_" + currFedID.ToString("X4") + " on BIMRL_ELEMENT_" + currFedID.ToString("X4")
// + " (BODY_MAJOR_AXIS_CENTROID) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_TOPOF_POLYGON_" + currFedID.ToString("X4") + " on BIMRL_TOPO_FACE_" + currFedID.ToString("X4")
// + " (POLYGON) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_TOPOF_NORMAL_" + currFedID.ToString("X4") + " on BIMRL_TOPO_FACE_" + currFedID.ToString("X4")
// + " (NORMAL) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
//sqlStmt = "Create Index IDX_TOPOF_CENTROID_" + currFedID.ToString("X4") + " on BIMRL_TOPO_FACE_" + currFedID.ToString("X4")
// + " (CENTROID) INDEXTYPE is MDSYS.SPATIAL_INDEX PARAMETERS ('sdo_indx_dims=3')";
//currStep = sqlStmt;
//cmd.CommandText = sqlStmt;
//cmd.ExecuteNonQuery();
DBOperation.commitTransaction();
}
catch (Exception e)
{
string excStr = "%%Error - " + e.Message + "\n\t" + currStep;
_bimrlCommon.StackPushError(excStr);
if (DBOperation.UIMode)
{
BIMRLErrorDialog errorDlg = new BIMRLErrorDialog(_bimrlCommon);
errorDlg.ShowDialog();
}
else
{
Console.Write(_bimrlCommon.ErrorMessages);
}
}
// There are entries in the error stack, show them at the end
if (_bimrlCommon.BIMRLErrorStackCount > 0)
{
if (DBOperation.UIMode)
{
BIMRLErrorDialog errorDlg = new BIMRLErrorDialog(_bimrlCommon);
errorDlg.ShowDialog();
}
else
{
Console.Write(_bimrlCommon.ErrorMessages);
}
}
}
public void createSpatialIndexFromBIMRLElement(int federatedId, string whereCond, bool createFaces = true, bool createSpIdx = true)
{
DBOperation.beginTransaction();
string currStep = string.Empty;
OracleCommand command = new OracleCommand(" ", DBOperation.DBConn);
OracleDataReader reader;
bool selectedRegen = false;
Point3D llb;
Point3D urt;
DBOperation.getWorldBB(federatedId, out llb, out urt);
Octree.WorldBB = new BoundingBox3D(llb, urt);
Octree.MaxDepth = DBOperation.OctreeSubdivLevel;
try
{
command.CommandText = "SELECT COUNT(*) FROM " + DBOperation.formatTabName("BIMRL_ELEMENT", federatedId) + " where geometrybody is not null ";
object rC = command.ExecuteScalar();
int totalRowCount = Convert.ToInt32(rC.ToString()) * (int)Math.Pow(8, 2);
string sqlStmt = "select elementid, elementtype, geometrybody from " + DBOperation.formatTabName("BIMRL_ELEMENT", federatedId) + " where geometrybody is not null ";
if (!string.IsNullOrEmpty(whereCond))
{
sqlStmt += " and " + whereCond;
selectedRegen = true;
}
currStep = sqlStmt;
command.CommandText = sqlStmt;
command.FetchSize = 20;
// The following is needed to update the element table with Bbox information
string sqlStmt3 = "UPDATE " + DBOperation.formatTabName("BIMRL_ELEMENT", federatedId) + " SET GeometryBody_BBOX = :bbox, "
+ "GeometryBody_BBOX_CENTROID = :cent WHERE ELEMENTID = :eid";
OracleCommand commandUpdBbox = new OracleCommand(" ", DBOperation.DBConn);
commandUpdBbox.CommandText = sqlStmt3;
commandUpdBbox.Parameters.Clear();
OracleParameter[] Bbox = new OracleParameter[3];
Bbox[0] = commandUpdBbox.Parameters.Add("bbox", OracleDbType.Object);
Bbox[0].Direction = ParameterDirection.Input;
Bbox[0].UdtTypeName = "MDSYS.SDO_GEOMETRY";
List <List <SdoGeometry> > bboxListList = new List <List <SdoGeometry> >();
List <SdoGeometry> bboxList = new List <SdoGeometry>();
Bbox[1] = commandUpdBbox.Parameters.Add("cent", OracleDbType.Object);
Bbox[1].Direction = ParameterDirection.Input;
Bbox[1].UdtTypeName = "MDSYS.SDO_GEOMETRY";
List <List <SdoGeometry> > centListList = new List <List <SdoGeometry> >();
List <SdoGeometry> centList = new List <SdoGeometry>();
Bbox[2] = commandUpdBbox.Parameters.Add("eid", OracleDbType.Varchar2);
Bbox[2].Direction = ParameterDirection.Input;
List <List <string> > eidUpdListList = new List <List <string> >();
List <string> eidUpdList = new List <string>();
int sublistCnt = 0;
// end for Bbox
Octree octreeInstance = null;
if (selectedRegen)
{
octreeInstance = new Octree(federatedId, totalRowCount, DBOperation.OctreeSubdivLevel);
}
else
{
octreeInstance = new Octree(federatedId, totalRowCount, DBOperation.OctreeSubdivLevel, false, true); // Since it is not selectedRegen, we will rebuild the entire tree, skip Dict regen for this case
}
reader = command.ExecuteReader();
while (reader.Read())
{
string elemID = reader.GetString(0);
string elemTyp = reader.GetString(1);
SdoGeometry sdoGeomData = reader.GetValue(2) as SdoGeometry;
Polyhedron geom;
if (!SDOGeomUtils.generate_Polyhedron(sdoGeomData, out geom))
{
continue; // if there is something not right, skip the geometry
}
// - Update geometry info with BBox information
SdoGeometry bbox = new SdoGeometry();
bbox.Dimensionality = 3;
bbox.LRS = 0;
bbox.GeometryType = (int)SdoGeometryTypes.GTYPE.POLYGON;
int gType = bbox.PropertiesToGTYPE();
double[] arrCoord = new double[6];
int[] elemInfoArr = { 1, (int)SdoGeometryTypes.ETYPE_SIMPLE.POLYGON_EXTERIOR, 1 };
arrCoord[0] = geom.boundingBox.LLB.X;
arrCoord[1] = geom.boundingBox.LLB.Y;
arrCoord[2] = geom.boundingBox.LLB.Z;
arrCoord[3] = geom.boundingBox.URT.X;
arrCoord[4] = geom.boundingBox.URT.Y;
arrCoord[5] = geom.boundingBox.URT.Z;
bbox.ElemArrayOfInts = elemInfoArr;
bbox.OrdinatesArrayOfDoubles = arrCoord;
bboxList.Add(bbox);
SdoGeometry centroid = new SdoGeometry();
centroid.Dimensionality = 3;
centroid.LRS = 0;
centroid.GeometryType = (int)SdoGeometryTypes.GTYPE.POINT;
gType = centroid.PropertiesToGTYPE();
SdoPoint cent = new SdoPoint();
cent.XD = geom.boundingBox.Center.X;
cent.YD = geom.boundingBox.Center.Y;
cent.ZD = geom.boundingBox.Center.Z;
centroid.SdoPoint = cent;
centList.Add(centroid);
eidUpdList.Add(elemID);
sublistCnt++;
// Set 1000 records as a threshold for interval commit later on
if (sublistCnt >= 500)
{
bboxListList.Add(bboxList);
centListList.Add(centList);
eidUpdListList.Add(eidUpdList);
sublistCnt = 0;
bboxList = new List <SdoGeometry>();
centList = new List <SdoGeometry>();
eidUpdList = new List <string>();
}
// We will skip large buildinglementproxy that has more than 5000 vertices
bool largeMesh = (string.Compare(elemTyp, "IFCBUILDINGELEMENTPROXY", true) == 0) && geom.Vertices.Count > 5000;
if ((createFaces && !largeMesh) || (createFaces && selectedRegen))
{
// - Process face information and create consolidated faces and store them into BIMRL_TOPO_FACE table
BIMRLGeometryPostProcess processFaces = new BIMRLGeometryPostProcess(elemID, geom, _refBIMRLCommon, federatedId, null);
processFaces.simplifyAndMergeFaces();
processFaces.insertIntoDB(false);
}
if (createSpIdx)
{
octreeInstance.ComputeOctree(elemID, geom);
}
}
reader.Dispose();
if (createSpIdx)
{
// Truncate the table first before reinserting the records
FederatedModelInfo fedModel = DBOperation.getFederatedModelByID(federatedId);
if (DBOperation.DBUserID.Equals(fedModel.FederatedID))
{
DBOperation.executeSingleStmt("TRUNCATE TABLE " + DBOperation.formatTabName("BIMRL_SPATIALINDEX", federatedId));
}
else
{
DBOperation.executeSingleStmt("DELETE FROM " + DBOperation.formatTabName("BIMRL_SPATIALINDEX"));
}
collectSpatialIndexAndInsert(octreeInstance, federatedId);
}
if (sublistCnt > 0)
{
bboxListList.Add(bboxList);
centListList.Add(centList);
eidUpdListList.Add(eidUpdList);
}
for (int i = 0; i < eidUpdListList.Count; i++)
{
Bbox[0].Value = bboxListList[i].ToArray();
Bbox[0].Size = bboxListList[i].Count;
Bbox[1].Value = centListList[i].ToArray();
Bbox[1].Size = centListList[i].Count;
Bbox[2].Value = eidUpdListList[i].ToArray();
Bbox[2].Size = eidUpdListList[i].Count;
commandUpdBbox.ArrayBindCount = eidUpdListList[i].Count; // No of values in the array to be inserted
int commandStatus = commandUpdBbox.ExecuteNonQuery();
DBOperation.commitTransaction();
}
if (!string.IsNullOrEmpty(whereCond) && createSpIdx)
{
command.CommandText = "UPDATE BIMRL_FEDERATEDMODEL SET MAXOCTREELEVEL=" + Octree.MaxDepth.ToString() + " WHERE FEDERATEDID=" + federatedId.ToString();
command.ExecuteNonQuery();
DBOperation.commitTransaction();
}
}
catch (OracleException e)
{
string excStr = "%%Read Error - " + e.Message + "\n\t" + currStep;
_refBIMRLCommon.StackPushError(excStr);
}
catch (SystemException e)
{
string excStr = "%%Read Error - " + e.Message + "\n\t" + currStep;
_refBIMRLCommon.StackPushError(excStr);
throw;
}
command.Dispose();
}
private SdoGeometry WriteMultiLineString(IGeometry geometry)
{
int dimension = GetGeometryDimension(geometry);
SdoGeometry sdoGeometry = new SdoGeometry();
sdoGeometry.GeometryType = (int)SdoGeometryTypes.GTYPE.MULTILINE;
sdoGeometry.Dimensionality = dimension;
sdoGeometry.LRS = 0;
sdoGeometry.Sdo_Srid = geometry.SRID;
int[] elements = new int[geometry.NumGeometries * ElementTupleSize];
int oordinatesOffset = 1;
double[] ordinates = new double[] { };
for (int i = 0; i < geometry.NumGeometries; i++)
{
elements[i * ElementTupleSize + 0] = oordinatesOffset;
elements[i * ElementTupleSize + 1] = (int)SdoGeometryTypes.ETYPE_SIMPLE.LINE;
elements[i * ElementTupleSize + 2] = 1;
ordinates = AppendCoordinates(ordinates, ConvertCoordinates(geometry.GetGeometryN(i).Coordinates, dimension));
oordinatesOffset = ordinates.Length + 1;
}
sdoGeometry.ElemArrayOfInts = elements;
sdoGeometry.OrdinatesArrayOfDoubles = ordinates;
sdoGeometry.PropertiesToGTYPE();
return sdoGeometry;
}
private SdoGeometry WriteMultiPolygon(IGeometry geometry)
{
int dimension = GetGeometryDimension(geometry);
SdoGeometry sdoGeometry = new SdoGeometry();
sdoGeometry.GeometryType = (int)SdoGeometryTypes.GTYPE.MULTIPOLYGON;
sdoGeometry.Dimensionality = dimension;
sdoGeometry.LRS = 0;
sdoGeometry.Sdo_Srid = geometry.SRID;
for (int i = 0; i < geometry.NumGeometries; i++)
{
try
{
IPolygon pg = (IPolygon)geometry.GetGeometryN(i);
AddPolygon(sdoGeometry, pg, dimension);
}
catch (Exception e)
{
throw new ApplicationException(
"Found geometry that was not a Polygon in MultiPolygon", e);
}
}
sdoGeometry.PropertiesToGTYPE();
return sdoGeometry;
}
private SdoGeometry WritePolygon(IGeometry geometry)
{
int dimension = GetGeometryDimension(geometry);
SdoGeometry sdoGeometry = new SdoGeometry();
sdoGeometry.GeometryType = (int)SdoGeometryTypes.GTYPE.POLYGON;
sdoGeometry.Dimensionality = dimension;
sdoGeometry.LRS = 0;
sdoGeometry.Sdo_Srid = geometry.SRID;
AddPolygon(sdoGeometry, geometry as IPolygon, dimension);
sdoGeometry.PropertiesToGTYPE();
return sdoGeometry;
}
public static void UpdateElementTransform(IfcStore _model, string projectNumber, string projectName)
{
DBOperation.beginTransaction();
DBOperation.commitInterval = 5000;
string currStep = string.Empty;
BIMRLCommon _refBIMRLCommon = new BIMRLCommon();
int commandStatus = -1;
int currInsertCount = 0;
OracleCommand command = new OracleCommand(" ", DBOperation.DBConn);
XbimMatrix3D m3D = new XbimMatrix3D();
if (string.IsNullOrEmpty(projectName))
{
projectName = projectNumber + " - Federated";
}
string modelName;
if (!string.IsNullOrEmpty(_model.FileName))
{
modelName = Path.GetFileNameWithoutExtension(_model.FileName);
}
else
{
modelName = projectNumber + " - " + projectName;
}
command.CommandText = "SELECT FEDERATEDID FROM bimrl_federatedmodel WHERE MODELNAME = '" + modelName + "' AND PROJECTNUMBER='" + projectNumber + "' AND PROJECTNAME='" + projectName + "'";
object oFedID = command.ExecuteScalar();
if (oFedID == null)
{
return;
}
int fedID = int.Parse(oFedID.ToString());
command.CommandText = "SELECT ELEMENTID, LINENO FROM " + DBOperation.formatTabName("bimrl_element", fedID) + " WHERE GEOMETRYBODY IS NOT NULL";
OracleDataReader reader = command.ExecuteReader();
SortedDictionary <int, string> elemList = new SortedDictionary <int, string>();
while (reader.Read())
{
string elemid = reader.GetString(0);
int lineNo = reader.GetInt32(1);
elemList.Add(lineNo, elemid);
}
reader.Close();
Xbim3DModelContext context = new Xbim3DModelContext(_model);
foreach (KeyValuePair <int, string> elemListItem in elemList)
{
//IEnumerable<XbimGeometryData> geomDataList = _model.GetGeometryData(elemListItem.Key, XbimGeometryType.TriangulatedMesh);
IIfcProduct product = _model.Instances[elemListItem.Key] as IIfcProduct;
IEnumerable <XbimShapeInstance> shapeInstances = context.ShapeInstancesOf(product).Where(x => x.RepresentationType == XbimGeometryRepresentationType.OpeningsAndAdditionsIncluded);
if (shapeInstances.Count() == 0)
{
continue; // SKip if the product has no geometry
}
XbimMeshGeometry3D prodGeom = new XbimMeshGeometry3D();
IXbimShapeGeometryData shapeGeom = context.ShapeGeometry(shapeInstances.FirstOrDefault().ShapeGeometryLabel);
//XbimModelExtensions.Read(prodGeom, shapeGeom.ShapeData, shapeInstances.FirstOrDefault().Transformation);
//XbimGeometryData sdoGeomData = geomDataList.First();
//m3D = sdoGeomData.Transform;
//m3D = XbimMatrix3D.FromArray(sdoGeomData.DataArray2); // Xbim 3.0 removes Tranform property
m3D = shapeInstances.FirstOrDefault().Transformation;
string sqlStmt = "update " + DBOperation.formatTabName("BIMRL_ELEMENT", fedID) + " set TRANSFORM_COL1=:1, TRANSFORM_COL2=:2, TRANSFORM_COL3=:3, TRANSFORM_COL4=:4"
+ " Where elementid = '" + elemListItem.Value + "'";
// int status = DBOperation.updateGeometry(sqlStmt, sdoGeomData);
currStep = sqlStmt;
command.CommandText = sqlStmt;
try
{
OracleParameter[] sdoGeom = new OracleParameter[4];
for (int i = 0; i < sdoGeom.Count(); ++i)
{
sdoGeom[i] = command.Parameters.Add((i + 1).ToString(), OracleDbType.Object);
sdoGeom[i].Direction = ParameterDirection.Input;
sdoGeom[i].UdtTypeName = "MDSYS.SDO_GEOMETRY";
sdoGeom[i].Size = 1;
}
SdoGeometry trcol1 = new SdoGeometry();
trcol1.Dimensionality = 3;
trcol1.LRS = 0;
trcol1.GeometryType = (int)SdoGeometryTypes.GTYPE.POINT;
int gType = trcol1.PropertiesToGTYPE();
SdoPoint trcol1V = new SdoPoint();
trcol1V.XD = m3D.M11;
trcol1V.YD = m3D.M12;
trcol1V.ZD = m3D.M13;
trcol1.SdoPoint = trcol1V;
sdoGeom[1].Value = trcol1;
SdoGeometry trcol2 = new SdoGeometry();
trcol2.Dimensionality = 3;
trcol2.LRS = 0;
trcol2.GeometryType = (int)SdoGeometryTypes.GTYPE.POINT;
gType = trcol2.PropertiesToGTYPE();
SdoPoint trcol2V = new SdoPoint();
trcol2V.XD = m3D.M21;
trcol2V.YD = m3D.M22;
trcol2V.ZD = m3D.M23;
trcol2.SdoPoint = trcol2V;
sdoGeom[2].Value = trcol2;
SdoGeometry trcol3 = new SdoGeometry();
trcol3.Dimensionality = 3;
trcol3.LRS = 0;
trcol3.GeometryType = (int)SdoGeometryTypes.GTYPE.POINT;
gType = trcol3.PropertiesToGTYPE();
SdoPoint trcol3V = new SdoPoint();
trcol3V.XD = m3D.M31;
trcol3V.YD = m3D.M32;
trcol3V.ZD = m3D.M33;
trcol3.SdoPoint = trcol3V;
sdoGeom[3].Value = trcol3;
SdoGeometry trcol4 = new SdoGeometry();
trcol4.Dimensionality = 3;
trcol4.LRS = 0;
trcol4.GeometryType = (int)SdoGeometryTypes.GTYPE.POINT;
gType = trcol4.PropertiesToGTYPE();
SdoPoint trcol4V = new SdoPoint();
trcol4V.XD = m3D.OffsetX;
trcol4V.YD = m3D.OffsetY;
trcol4V.ZD = m3D.OffsetZ;
trcol4.SdoPoint = trcol4V;
sdoGeom[4].Value = trcol4;
commandStatus = command.ExecuteNonQuery();
command.Parameters.Clear();
currInsertCount++;
if (currInsertCount % DBOperation.commitInterval == 0)
{
//Do commit at interval but keep the long transaction (reopen)
DBOperation.commitTransaction();
}
}
catch (OracleException e)
{
string excStr = "%%Error - " + e.Message + "\n\t" + currStep;
_refBIMRLCommon.StackPushError(excStr);
//command.Dispose(); // Log Oracle error and continue
command = new OracleCommand(" ", DBOperation.DBConn);
// throw;
}
catch (SystemException e)
{
string excStr = "%%Insert Error - " + e.Message + "\n\t" + currStep;
_refBIMRLCommon.StackPushError(excStr);
throw;
}
DBOperation.commitTransaction();
command.Dispose();
}
}
