Georeferencing utility to convert an angle in degrees to X Axis Abscissa and Ordinate

src/blenderbim/blenderbim/bim/module/georeference/__init__.py

@@ -33,6 +33,7 @@ classes = (
     operator.ConvertGlobalToLocal,
     operator.GetCursorLocation,
     operator.SetCursorLocation,
+    operator.ConvertAngleToCoordinates,
     prop.BIMGeoreferenceProperties,
     ui.BIM_PT_gis,
     ui.BIM_PT_gis_utilities,

src/blenderbim/blenderbim/bim/module/georeference/operator.py

@@ -171,3 +171,19 @@ class ConvertGlobalToLocal(bpy.types.Operator, tool.Ifc.Operator):
 
     def _execute(self, context):
         core.convert_global_to_local(tool.Georeference)
+
+class ConvertAngleToCoordinates(bpy.types.Operator, tool.Ifc.Operator):
+    bl_idname = "bim.convert_angle_to_coord"
+    bl_label = "Convert Angle To Y Axis"
+    bl_options = {"REGISTER", "UNDO"}
+    bl_description = "Convert angle to Y axis"
+    type: bpy.props.StringProperty()
+
+    @classmethod
+    def poll(cls, context):
+        file = tool.Ifc.get()
+        props = context.scene.BIMGeoreferenceProperties
+        return file and (props.angle_degree_input_x or props.angle_degree_input_y)
+
+    def _execute(self, context):
+        core.convert_angle_to_coord(tool.Georeference, type=self.type)

src/blenderbim/blenderbim/bim/module/georeference/prop.py

@@ -37,6 +37,12 @@ class BIMGeoreferenceProperties(PropertyGroup):
     projected_crs: CollectionProperty(name="Projected CRS", type=Attribute)
     coordinate_input: StringProperty(name="Coordinate Input", description='Formatted "x,y,z" (without quotes)')
     coordinate_output: StringProperty(name="Coordinate Output", description='Formatted "x,y,z" (without quotes)')
+    angle_degree_input_x: FloatProperty(name="Angle Degree Input", description="Angle (in degrees) rel to Easting")
+    angle_degree_input_y: FloatProperty(name="Angle Degree Input", description="Angle (in degrees) rel to +Y")
+    x_axis_abscissa_output: StringProperty(name="X Axis Abscissa Ordinate Output", description="X axis abscissa and ordinate", )
+    x_axis_ordinate_output: StringProperty(name="X Axis Abscissa Ordinate Output", description="X axis abscissa and ordinate", )
+    y_axis_abscissa_output: StringProperty(name="Y Axis Abscissa Ordinate Output", description="Y axis abscissa and ordinate", )
+    y_axis_ordinate_output: StringProperty(name="Y Axis Abscissa Ordinate Output", description="Y axis abscissa and ordinate", )
     has_blender_offset: BoolProperty(name="Has Blender Offset")
     blender_eastings: StringProperty(name="Blender Eastings", default="0")
     blender_northings: StringProperty(name="Blender Northings", default="0")

src/blenderbim/blenderbim/bim/module/georeference/ui.py

@@ -167,3 +167,21 @@ class BIM_PT_gis_utilities(Panel):
         row = self.layout.row(align=True)
         row.operator("bim.convert_local_to_global", text="Local to Global")
         row.operator("bim.convert_global_to_local", text="Global to Local")
+
+        row = self.layout.row(align=True)
+        row.label(text="Orientation Calculator", icon="TRACKING_REFINE_FORWARDS")
+        row = self.layout.row(align=True)
+        row.prop(props, "angle_degree_input_x", text="Angle (°) rel to Easting")
+        row.operator("bim.convert_angle_to_coord", text="", icon="FILE_REFRESH").type = "rel_x"
+        row = self.layout.row(align=True)
+        row.prop(props, "x_axis_abscissa_output", text="XAxis Abscissa")
+        row = self.layout.row(align=True)
+        row.prop(props, "x_axis_ordinate_output", text="XAxis Ordinate")
+
+        row = self.layout.row(align=True)
+        row.prop(props, "angle_degree_input_y", text="Angle (°) rel to to +Y")
+        row.operator("bim.convert_angle_to_coord", text="", icon="FILE_REFRESH").type = "rel_y"
+        row = self.layout.row(align=True)
+        row.prop(props, "y_axis_abscissa_output", text="North Abscissa")
+        row = self.layout.row(align=True)
+        row.prop(props, "y_axis_ordinate_output", text="North Ordinate")

src/blenderbim/blenderbim/core/georeference.py

@@ -80,3 +80,7 @@ def convert_global_to_local(georeference):
     coordinates = georeference.enh2xyz(georeference.get_coordinates("input"), georeference.get_map_conversion())
     georeference.set_coordinates("output", coordinates)
     georeference.set_cursor_location(coordinates)
+
+def convert_angle_to_coord(georeference, type):
+    vector_coordinates = georeference.angle2coords(georeference.get_angle(type), type)
+    georeference.set_vector_coordinates(vector_coordinates,type)

src/blenderbim/blenderbim/core/tool.py

@@ -370,13 +370,15 @@ class Geometry:
 
 @interface
 class Georeference:
+    def angle2coords(cls, angle, type): pass
     def disable_editing(cls): pass
     def enable_editing(cls): pass
     def enh2xyz(cls, map_conversion, coordinates): pass
+    def get_angle(cls, type): pass
     def get_coordinates(cls, io): pass
     def get_cursor_location(cls): pass
-    def get_map_conversion(cls): pass
     def get_map_conversion_attributes(cls): pass
+    def get_map_conversion(cls): pass
     def get_projected_crs_attributes(cls): pass
     def get_true_north_attributes(cls): pass
     def import_map_conversion(cls): pass
@@ -388,6 +390,7 @@ class Georeference:
     def set_cursor_location(cls, coordinates): pass
     def set_ifc_grid_north(cls): pass
     def set_ifc_true_north(cls): pass
+    def set_vector_coordinates(cls, vector_coordinates, type): pass
     def xyz2enh(cls, map_conversion, coordinates): pass
 
 

src/blenderbim/blenderbim/tool/georeference.py

@@ -272,3 +272,27 @@ class Georeference(blenderbim.core.tool.Georeference):
             float(props.map_conversion.get("XAxisOrdinate").string_value),
         )
         bpy.context.scene.sun_pos_properties.north_offset = -radians(angle)
+
+    @classmethod
+    def angle2coords(cls, angle, type):
+        if type == "rel_x":
+            return ifcopenshell.util.geolocation.angle2xaxis(angle)
+        elif type == "rel_y":
+            return ifcopenshell.util.geolocation.angle2yaxis(angle)
+
+    @classmethod
+    def get_angle(cls, type):
+        if type == "rel_x":
+            return bpy.context.scene.BIMGeoreferenceProperties.angle_degree_input_x
+        elif type == "rel_y":
+            return bpy.context.scene.BIMGeoreferenceProperties.angle_degree_input_y
+
+    @classmethod
+    def set_vector_coordinates(cls, vector_coordinates, type):
+        x, y = vector_coordinates
+        if type == "rel_x":
+            bpy.context.scene.BIMGeoreferenceProperties.x_axis_abscissa_output = str(x)
+            bpy.context.scene.BIMGeoreferenceProperties.x_axis_ordinate_output = str(y)
+        elif type == "rel_y":
+            bpy.context.scene.BIMGeoreferenceProperties.y_axis_abscissa_output = str(x)
+            bpy.context.scene.BIMGeoreferenceProperties.y_axis_ordinate_output = str(y)

src/ifcopenshell-python/ifcopenshell/util/geolocation.py

@@ -187,3 +187,18 @@ def get_true_north(ifc_file):
     except:
         return 0
     return 0
+
+
+# Used  for converting an angle in degrees to return the X and Y vectors of the X Axis in IFC grid north geolocation:
+def angle2xaxis(angle):
+    angle_rad = math.radians(angle)
+    x = math.cos(angle_rad)
+    y = - math.sin(angle_rad)
+    return x, y
+
+# Used for converting True North angle as seen in CAD (relative to +Y)
+def angle2yaxis(angle):
+    angle_rad = math.radians(angle)
+    x = - math.sin(angle_rad)
+    y = math.cos(angle_rad)
+    return x, y