In [1]:
Copied!
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 Lines¶
Lines can be specified using a list of tuples. A spatial reference system should also be applied, as the coordinates cannot be interpreted without it.
In [2]:
Copied!
# Create a LineString geometry from a list of coordinate tuples.
# Define a projection for Aachen area
aachen_centered_srs = geokit.core.srs.centeredLAEA(6.083, 50.775)
# Create a linestring geometry using a list of xy-coordinate tuples
house_line = geokit.core.geom.line(
[
(0, 0),
(40000, 0),
(40000, 30000),
(20000, 50000),
(0, 30000),
(0, 0),
],
srs=aachen_centered_srs,
)
print(type(house_line))
# Create a LineString geometry from a list of coordinate tuples.
# Define a projection for Aachen area
aachen_centered_srs = geokit.core.srs.centeredLAEA(6.083, 50.775)
# Create a linestring geometry using a list of xy-coordinate tuples
house_line = geokit.core.geom.line(
[
(0, 0),
(40000, 0),
(40000, 30000),
(20000, 50000),
(0, 30000),
(0, 0),
],
srs=aachen_centered_srs,
)
print(type(house_line))
<class 'osgeo.ogr.Geometry'>
Visualize Line¶
In [3]:
Copied!
ax_handle = geokit.core.geom.drawGeoms(
geoms=house_line,
figsize=(4, 4),
srs=aachen_centered_srs,
)
ax_handle = geokit.core.geom.drawGeoms(
geoms=house_line,
figsize=(4, 4),
srs=aachen_centered_srs,
)
Store Line as Vector File¶
In [4]:
Copied!
# Simple vector creation
geokit.core.vector.createVector(geoms=house_line, output="example_lines.shp")
# Simple vector creation
geokit.core.vector.createVector(geoms=house_line, output="example_lines.shp")
Out[4]:
'/home/docs/checkouts/readthedocs.org/user_builds/geokit/checkouts/latest/docs/Examples/_03_geometries/_01_create_geometries/example_lines.shp'