Create Points¶

This example shows to create, visualize and store points using GeoKit.

Import required packages¶

In [1]:

import geokit.core.geom
import geokit.core.vector
import geokit.core.srs

import geokit.core.geom import geokit.core.vector import geokit.core.srs

Create Two Dimensional Points¶

Points can be specified either from tuples or using two coordinates. A spatial reference system should also be applied, as the coordinates cannot be interpreted without it.

In [2]:

# Create individual point geometries using either a coordinate tuple or two separate coordinate arguments.

aachen_centered_srs = geokit.core.srs.centeredLAEA(6.083, 50.775)
# Tuple representing the coordinates of your points (in meters)
point_1 = (0, 0)
point_2 = (40000, 0)
point_3 = (40000, 30000)

# Create point geometry from tuple
aachen_pt_1_from_tuple = geokit.core.geom.point(point_1, srs=aachen_centered_srs)
aachen_pt_2_from_tuple = geokit.core.geom.point(point_2, srs=aachen_centered_srs)
aachen_pt_3_from_tuple = geokit.core.geom.point(point_3, srs=aachen_centered_srs)

print(type(aachen_pt_1_from_tuple))

# Create point geometry from separate x and y values (in longitude and latitude)
aachen_pt_4_from_two_arguments = geokit.core.geom.point(6.373, 51.225, srs=4326)
aachen_pt_5_from_two_arguments = geokit.core.geom.point(6.083, 51.045, srs=4326)


# Define a projection for the Aachen area
aachen_centered_srs = geokit.core.srs.centeredLAEA(6.083, 50.775)

# Transform the points to the defined projection
aachen_pt_4_from_two_arguments = geokit.core.geom.transform(aachen_pt_4_from_two_arguments, aachen_centered_srs)
aachen_pt_5_from_two_arguments = geokit.core.geom.transform(aachen_pt_5_from_two_arguments, aachen_centered_srs)

list_of_geometry_points = [
    aachen_pt_1_from_tuple,
    aachen_pt_2_from_tuple,
    aachen_pt_3_from_tuple,
    aachen_pt_4_from_two_arguments,
    aachen_pt_5_from_two_arguments,
]

# Create individual point geometries using either a coordinate tuple or two separate coordinate arguments. aachen_centered_srs = geokit.core.srs.centeredLAEA(6.083, 50.775) # Tuple representing the coordinates of your points (in meters) point_1 = (0, 0) point_2 = (40000, 0) point_3 = (40000, 30000) # Create point geometry from tuple aachen_pt_1_from_tuple = geokit.core.geom.point(point_1, srs=aachen_centered_srs) aachen_pt_2_from_tuple = geokit.core.geom.point(point_2, srs=aachen_centered_srs) aachen_pt_3_from_tuple = geokit.core.geom.point(point_3, srs=aachen_centered_srs) print(type(aachen_pt_1_from_tuple)) # Create point geometry from separate x and y values (in longitude and latitude) aachen_pt_4_from_two_arguments = geokit.core.geom.point(6.373, 51.225, srs=4326) aachen_pt_5_from_two_arguments = geokit.core.geom.point(6.083, 51.045, srs=4326) # Define a projection for the Aachen area aachen_centered_srs = geokit.core.srs.centeredLAEA(6.083, 50.775) # Transform the points to the defined projection aachen_pt_4_from_two_arguments = geokit.core.geom.transform(aachen_pt_4_from_two_arguments, aachen_centered_srs) aachen_pt_5_from_two_arguments = geokit.core.geom.transform(aachen_pt_5_from_two_arguments, aachen_centered_srs) list_of_geometry_points = [ aachen_pt_1_from_tuple, aachen_pt_2_from_tuple, aachen_pt_3_from_tuple, aachen_pt_4_from_two_arguments, aachen_pt_5_from_two_arguments, ]

<class 'osgeo.ogr.Geometry'>

Visualize Points¶

In [3]:

ax_handle = geokit.core.geom.drawGeoms(
    geoms=list_of_geometry_points,
    figsize=(4, 4),
    srs=aachen_centered_srs,
)

ax_handle = geokit.core.geom.drawGeoms( geoms=list_of_geometry_points, figsize=(4, 4), srs=aachen_centered_srs, )

Store Points as Vector File¶

In [4]:

# Simple vector creation
geokit.core.vector.createVector(geoms=list_of_geometry_points, output="example_points.shp")

# Simple vector creation geokit.core.vector.createVector(geoms=list_of_geometry_points, output="example_points.shp")

Out[4]:

'/home/docs/checkouts/readthedocs.org/user_builds/geokit/checkouts/latest/docs/Examples/_03_geometries/_01_create_geometries/example_points.shp'