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.

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: ObjectID

close(): Closes the object and saves the changes to the Project, preventing any further changes.

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().

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: BaseDrillholeTable
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

DrillholeDatabaseField

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 currently associated with this object.

Returns: The ID of the colour map object or null object ID if there is no colour map.
Return type: ObjectID

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 (ObjectID[ColourMap] | ColourMap) – 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 displayed field in the displayed table. Thus these are a property of the visualisation of the drillhole and not the drillhole itself. The points will change if the visualised field is changed.

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.

property edges: ndarray

The edges used to visualise the drillhole.

These edges mark the boundaries of intervals for displayed field in the displayed table. Thus these are a property of the visualisation of the drillhole and not the drillhole itself. The edges will change if the visualised field is changed.

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.

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[BaseDrillholeTable]

The tables representing the drillhole.

Returns: List of BaseDrillholeTable for the drillhole.
Return type: list

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: BaseDrillholeTable
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 coordinate_system: CoordinateSystem | None

The coordinate system the points of this object are in.

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.

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)

property extent: Extent: The axes aligned bounding extent of the object.

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 lock_type: LockType

Indicates whether operating in read-only or read-write mode.

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]]

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]

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 – The name of the object attribute for which the value should be set.
dtype – 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 – 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.
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.