mapteksdk.geologycore.drillholes module
Drillhole data types.
- class Drillhole(object_id=None, lock_type=LockType.READ)
Bases:
Topology
,TablesMixin
Class representing a single drillhole.
Drillholes cannot be created with Project.new() like other objects. Instead they are created through DrillholeDatabase.new_drillhole(). An error will be raised if you open the new drillhole before closing the drillhole database.
Notes
Though the values representing a Drillhole are accessed through this object, they are actually stored in the DrillholeDatabase which the Drillhole is inside.
Drillholes in the same database will have the same tables and fields.
- Raises
TypeError – If passed to Project.new().
OrphanDrillholeError – If the drillhole is not inside of a drillhole database container. This error will also be raised if a new drillhole is opened before the drillhole database is closed.
- Parameters
object_id (ObjectID) –
lock_type (LockType) –
See also
mapteksdk.geologycore.database.DrillholeDatabase.new_drillhole
Create new drillholes.
- drillholes
Help page for this class.
- classmethod static_type()
Return the type of a topology as stored in a Project.
This can be used for determining if the type of an object is topology.
- property id: ObjectID[Drillhole]
Object ID that uniquely references this object in the project.
- Returns
The unique id of this object.
- Return type
- close()
Closes the object and saves the changes to the Project.
Attempting to read or edit properties of an object after closing it will raise a ReadOnlyError.
- property name: str
The name of the drillhole.
- property converted_collar: ndarray
The collar point of the drillhole adjusted for the coordinate system.
This is the collar point adjusted for the drillhole’s coordinate system and thus where the collar point appears when the drillhole is viewed in the application.
This will not match the ordinates read from the collar table (and thus the raw_collar property) if the drillhole has a coordinate system.
- Raises
TableNotFoundError – If there is no collar table.
TooManyTablesError – If there are multiple collar tables in the database.
Warning
If the collar table or the raw_collar property is edited, changes will only be reflected in this property after save() is called.
Notes
If there is no elevation field in the collar table, the Z ordinate will always be zero.
- property raw_collar: ndarray
The collar point as it appears in the collar table.
This does not take into account the drillhole’s coordinate system.
This will not match the converted_collar property if the drillhole has a coordinate system.
- Warns
RuntimeWarning – If the table contains no elevation field and the collar is set to a point containing a non-zero Z ordinate.
Notes
If there is no elevation field in the collar table, the Z ordinate will always be zero.
Examples
If the table has an elevation field, when assigning to this property the caller should provide a northing, easting and elevation value:
>>> drillhole.raw_collar = [1.1, 2.2, 3.3]
If the table does not have an elevation field, when assigning to this property the elevation may be omitted:
>>> drillhole.raw_collar = [1.1, 2.2]
- property displayed_table: BaseDrillholeTable
Returns the displayed table for the drillhole.
This is the table which contains the displayed field used to colour the intervals of the drillhole in views.
- Returns
The displayed table.
- Return type
- Raises
TableNotFoundError – If the displayed table cannot be found.
- property displayed_field: DrillholeDatabaseField
Returns the displayed field for the drillhole.
When the drillhole is displayed in the application, this is the field which is used to determine the intervals used to colour the drillhole. This is also the field used to determine the points and edges returned by the drillhole’s point and edge properties.
- Returns
The displayed field.
- Return type
- Raises
TableNotFoundError – If the displayed table cannot be found.
FieldNotFoundError – If the displayed field cannot be found.
- get_colour_map()
Return the ID of the colour map object associated with this object.
- Returns
The ID of the colour map object or null object ID if there is no colour map.
- Return type
- set_visualisation(field, colour_map)
Set the field and colour map used to display the Drillhole.
- Parameters
field (DrillholeDatabaseField) – Field to use to display the Drillhole.
colour_map (Union[ObjectID[Union[StringColourMap, NumericColourMap]], StringColourMap, NumericColourMap]) – ObjectID of the colour map to use to display the Drillhole, or the colour map object itself.
- Raises
ReadOnlyError – If the drillhole is open for read-only.
TypeError – If colour_map is not an ObjectID.
TypeError – If field is not a DrillholeDatabaseField.
TypeError – If field contains numeric values and colour_map is a StringColourMap.
TypeError – If field contains string values and colour_map is a NumericColourMap.
ValueError – If field stores boolean values.
ValueError – If field is not a field read from this drillhole.
- property points: ndarray
The points used to visualise the drillhole.
These points mark the boundaries of intervals for the displayed field in the displayed table. Thus these are a property of the visualisation of the drillhole and not the drillhole itself. This array cannot be edited directly, however edits to the displayed field will propagate to this property when the drillhole is saved.
Notes
The points of a drillhole are derived from the displayed table and cannot be edited directly from Python.
- property point_selection: ndarray
Point selection array for the drillhole.
If point_selection[i] is True then the point located at points[i] is selected.
Notes
The point selection of a drillhole cannot be edited from Python.
- property edges: ndarray
The edges used to visualise the drillhole.
The edges which represent the intervals of the displayed field in the displayed table. Thus these are a property of the visualisation of the drillhole and not the drillhole itself. This array cannot be edited directly, however edits to the displayed field will propagate to this property when the drillhole is saved.
Notes
The edges of a drillhole are derived from the displayed table and cannot be edited directly from Python.
- property edge_selection: ndarray
Edge selection array for the drillhole.
If edge_selection[i] is True then the edge located at edges[i] is selected.
Notes
The edge selection of a drillhole cannot be edited from Python.
- property assay_table: AssayTable
Returns the assay table if it exists.
This property should only be used if the caller is certain that the database only contains one assay table.
- Raises
TableNotFoundError – If there is no assay table.
- Warns
TooManyTablesWarning – If the database contains multiple assay tables. The first table was returned.
- property collar_table: CollarTable
Returns the collar table if it exists.
The collar table represents the location on the surface which the drillhole was taken from.
A database can only have one collar table.
- Raises
TableNotFoundError – If there is no collar table.
TooManyTablesError – If there are multiple collar tables in the database. This should not be possible.
- property survey_table: SurveyTable
Returns the survey table if it exists.
A database can only contain one survey table.
- Raises
TableNotFoundError – If there is no survey table.
TooManyTablesError – If there are multiple survey tables in the database. This should not be possible.
- property geology_table: GeologyTable
Returns the geology table if it exists.
A database may contain multiple geology tables. This property should only be used if the caller is certain that the database only contains one geology table.
- Raises
TableNotFoundError – If there is no geology table.
- Warns
TooManyTablesWarning – If the database contains multiple geology tables. The first table was returned.
- property downhole_table: DownholeTable
Returns the downhole table if it exists.
A database may contain multiple downhole tables. This property should only be used if the caller is certain that the database only contains one downhole table.
- Raises
TableNotFoundError – If there is no geology table.
- Warns
TooManyTablesWarning – If the database contains multiple geology tables. The first table was returned.
- property quality_table: QualityTable
Returns the quality table if it exists.
A database may contain multiple quality tables. This property should only be used if the caller is certain that the database only contains one quality table.
- Raises
TableNotFoundError – If there is no geology table.
- Warns
TooManyTablesWarning – If the database contains multiple geology tables. The first table was returned.
- property tables: list[mapteksdk.geologycore.tables.BaseDrillholeTable]
The tables representing the drillhole.
- Returns
List of BaseDrillholeTable for the drillhole.
- Return type
list
- point_at_depth(depth)
Get the point in the drillhole at the specified depth.
The returned point is relative to the collar point. To get the point at the specified depth in world coordinates, you must add the collar point to this point.
This takes into account the desurvey method assigned to the drillhole database and any values in the survey table.
- Parameters
depth (float) – The depth down the hole for which the point should be returned.
- Returns
Numpy array of shape (3,) containing the X, Y and Z coordinate of the point at the specified depth down the hole. This may not correspond to any point in the points array.
- Return type
np.ndarray
- Raises
SurveyTableLoadedError – If the drillhole is open for editing and may have unsaved changes.
- tables_by_type(table_type)
Returns a list of tables with the specified type.
- Parameters
table_type (DrillholeTableType) – The type of table to include in the list.
- Returns
List of BaseDrillholeTable objects with the specified table type.
- Return type
list
- Raises
KeyError – If table_type is not a DrillholeTableType.
- table_by_name(name)
Returns the table with the specified name.
- Parameters
name (str) – The name of the table to return.
- Returns
The table with the specified name.
- Return type
- Raises
TableNotFoundError – If there is no table with the specified name.
- attribute_names()
Returns a list containing the names of all object-level attributes.
Use this to iterate over the object attributes.
- Returns
List containing the attribute names.
- Return type
list
Examples
Iterate over all object attributes of the object stared at “target” and print their values.
>>> from mapteksdk.project import Project >>> project = Project() >>> with project.read("target") as read_object: ... for name in read_object.attribute_names(): ... print(name, ":", read_object.get_attribute(name))
- property closed: bool
If this object has been closed.
Attempting to read or edit a closed object will raise an ObjectClosedError. Such an error typically indicates an error in the script and should not be caught.
Examples
If the object was opened with the Project.new(), Project.edit() or Project.read() in a “with” block, this will be True until the with block is closed and False afterwards.
>>> with self.project.new("cad/point_set", PointSet) as point_set: >>> point_set.points = [[1, 2, 3], [4, 5, 6]] >>> print("closed?", point_set.closed) >>> print("closed?", point_set.closed) closed? False closed? True
- property coordinate_system: mapteksdk.data.coordinate_systems.CoordinateSystem | None
The coordinate system the points of this object are in.
- Raises
ReadOnlyError – If set on an object open for read-only.
Warning
Setting this property does not change the points. This is only a label stating the coordinate system the points are in.
Notes
If the object has no coordinate system, this will be None.
Changes are done directly in the project and will not be undone if an error occurs.
Examples
Creating an edge network and setting the coordinate system to be WGS84. Note that setting the coordinate system does not change the points. It is only stating which coordinate system the points are in.
>>> from pyproj import CRS >>> from mapteksdk.project import Project >>> from mapteksdk.data import Polygon >>> project = Project() >>> with project.new("cad/rectangle", Polygon) as new_edges: ... # Coordinates are in the form [longitude, latitude] ... new_edges.points = [[112, 9], [112, 44], [154, 44], [154, 9]] ... new_edges.coordinate_system = CRS.from_epsg(4326)
Often a standard map projection is not convenient or accurate for a given application. In such cases a local transform can be provided to allow coordinates to be specified in a more convenient system. The below example defines a local transform where the origin is translated 1.2 degrees north and 2.1 degree east, points are scaled to be twice as far from the horizontal origin and the coordinates are rotated 45 degrees clockwise about the horizontal_origin. Note that the points of the polygon are specified in the coordinate system after the local transform has been applied.
>>> import math >>> from pyproj import CRS >>> from mapteksdk.project import Project >>> from mapteksdk.data import Polygon, CoordinateSystem, LocalTransform >>> project = Project() >>> transform = LocalTransform( ... horizontal_origin = [1.2, 2.1], ... horizontal_scale_factor = 2, ... horizontal_rotation = math.pi / 4) >>> system = CoordinateSystem(CRS.from_epsg(20249), transform) >>> with project.new("cad/rectangle_transform", Polygon) as new_edges: ... new_edges.points = [[112, 9], [112, 44], [154, 44], [154, 9]] ... new_edges.coordinate_system = system
See also
mapteksdk.data.coordinate_systems.CoordinateSystem
Allows for a coordinate system to be defined with an optional local transform.
- property created_date: datetime
The date and time (in UTC) of when this object was created.
- Returns
The date and time the object was created. 0:0:0 1/1/1970 if the operation failed.
- Return type
datetime.datetime
- delete_all_attributes()
Delete all object attributes attached to an object.
This only deletes object attributes and has no effect on PrimitiveAttributes.
- Raises
RuntimeError – If all attributes cannot be deleted.
- delete_attribute(attribute)
Deletes a single object-level attribute.
Deleting a non-existent object attribute will not raise an error.
- Parameters
attribute (str) – Name of attribute to delete.
- Returns
True if the object attribute existed and was deleted; False if the object attribute did not exist.
- Return type
bool
- Raises
RuntimeError – If the attribute cannot be deleted.
- dissociate_raster(raster)
Removes the raster from the object.
This is done directly on the Project and will not be undone if an error occurs.
- Parameters
raster (Raster | ObjectID[Raster]) – The raster to dissociate.
- Returns
True if the raster was successfully dissociated from the object, False if the raster was not associated with the object.
- Return type
bool
- Raises
TypeError – If raster is not a Raster.
ReadOnlyError – If this object is open for read-only.
Notes
This only removes the association between the Raster and the object, it does not clear the registration information from the Raster.
Examples
Dissociate the first raster found on a picked object.
>>> from mapteksdk.project import Project >>> from mapteksdk import operations >>> project = Project() >>> oid = operations.object_pick( ... support_label="Pick an object to remove a raster from.") ... with project.edit(oid) as data_object: ... report = f"There were no raster to remove from {oid.path}" ... for index in data_object.rasters: ... data_object.dissociate_raster(data_object.rasters[index]) ... report = f"Removed raster {index} from {oid.path}" ... break ... # Now that the raster is dissociated and the object is closed, ... # the raster can be associated with a different object. ... operations.write_report("Remove Raster", report)
- get_attribute(name)
Returns the value for the attribute with the specified name.
- Parameters
name (str) – The name of the object attribute to get the value for.
- Returns
The value of the object attribute name. For dtype = datetime.datetime this is an integer representing the number of milliseconds since 1st Jan 1970. For dtype = datetime.date this is a tuple of the form: (year, month, day).
- Return type
ObjectAttributeTypes
- Raises
KeyError – If there is no object attribute called name.
Warning
In the future this function may be changed to return datetime.datetime and datetime.date objects instead of the current representation for object attributes of type datetime.datetime or datetime.date.
- get_attribute_type(name)
Returns the type of the attribute with the specified name.
- Parameters
name (str) – Name of the attribute whose type should be returned.
- Returns
The type of the object attribute name.
- Return type
ObjectAttributeDataTypes
- Raises
KeyError – If there is no object attribute called name.
- property is_read_only: bool
If this object is read-only.
This will return True if the object was open with Project.read() and False if it was open with Project.edit() or Project.new(). Attempting to edit a read-only object will raise an error.
- property lock_type: LockType
Indicates whether operating in read-only or read-write mode.
Use the is_read_only property instead for checking if an object is open for reading or editing.
- Returns
The type of lock on this object. This will be LockType.ReadWrite if the object is open for editing and LockType.Read if the object is open for reading.
- Return type
LockType
- property modified_date: datetime
The date and time (in UTC) of when this object was last modified.
- Returns
The date and time this object was last modified. 0:0:0 1/1/1970 if the operation failed.
- Return type
datetime.datetime
- property rasters: dict[int, ObjectID[Raster]]
The raster associated with this object.
This is a dictionary of raster indices and Object IDs of the raster images currently associated with this object.
The keys are the raster ids and the values are the Object IDs of the associated rasters. Note that all raster ids are integers however they may not be consecutive - for example, an object may have raster ids 0, 1, 5 and 200.
Notes
Rasters with higher indices appear on top of rasters with lower indices. The maximum possible raster id is 255.
Removing a raster from this dictionary will not remove the raster association from the object. Use dissociate_raster to do this.
Examples
Iterate over all rasters on an object and invert the colours. Note that this will fail if there is no object at the path “target” and it will do nothing if no rasters are associated with the target.
>>> from mapteksdk.project import Project >>> project = Project() >>> with project.read("target") as read_object: ... for raster in read_object.rasters.values(): ... with project.edit(raster) as edit_raster: ... edit_raster.pixels[:, :3] = 255 - edit_raster.pixels[:, :3]
- save()
Save the changes made to the object.
Generally a user does not need to call this function, because it is called automatically at the end of a with block using Project.new() or Project.edit().
- set_attribute(name, dtype, data)
Sets the value for the object attribute with the specified name.
This will overwrite any existing attribute with the specified name.
- Parameters
name (str) – The name of the object attribute for which the value should be set.
dtype (type[Union[NoneType, Type[NoneType], ctypes.c_bool, ctypes.c_byte, ctypes.c_ubyte, ctypes.c_short, ctypes.c_ushort, ctypes.c_long, ctypes.c_ulong, ctypes.c_longlong, ctypes.c_ulonglong, ctypes.c_float, ctypes.c_double, ctypes.c_char_p, datetime.datetime, datetime.date]] | None) – The type of data to assign to the attribute. This should be a type from the ctypes module or datetime.datetime or datetime.date. Passing bool is equivalent to passing ctypes.c_bool. Passing str is equivalent to passing ctypes.c_char_p. Passing int is equivalent to passing ctypes.c_int16. Passing float is equivalent to passing ctypes.c_double.
data (Any) – The value to assign to object attribute name. For dtype = datetime.datetime this can either be a datetime object or timestamp which will be passed directly to datetime.utcfromtimestamp(). For dtype = datetime.date this can either be a date object or a tuple of the form: (year, month, day).
- Raises
ValueError – If dtype is an unsupported type.
TypeError – If value is an inappropriate type for object attribute name.
ValueError – If name starts or ends with whitespace or is empty.
RuntimeError – If a different error occurs.
Warning
Object attributes are saved separately from the object itself - any changes made by this function (assuming it does not raise an error) will be saved even if save() is not called (for example, due to an error being raised by another function).
Examples
Create an object attribute on an object at “target” and then read its value.
>>> import ctypes >>> from mapteksdk.project import Project >>> project = Project() >>> with project.edit("target") as edit_object: ... edit_object.set_attribute("count", ctypes.c_int16, 0) ... with project.read("target") as read_object: ... print(read_object.get_attribute("count")) 0
- property table_count: int
The number of tables in the database.