Create Polygon¶

This example shows to create, visualize and store a polygon 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 Polygon¶

Similar to the line object polygons are specified using a list of tuples. A spatial reference system should also be applied, as the coordinates cannot be interpreted without it.

In [2]:

# Create a Polygon geometry from a list of coordinate tuples.

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

aachen_polygon = geokit.core.geom.polygon(
    [
        (0, 0),
        (40000, 0),
        (40000, 30000),
        (20000, 50000),
        (0, 30000),
        (0, 0),
    ],
    srs=4326,
)

print(type(aachen_polygon))

# Create a Polygon geometry from a list of coordinate tuples. # Define a projection for the Aachen area aachen_centered_srs = geokit.core.srs.centeredLAEA(6.083, 50.775) aachen_polygon = geokit.core.geom.polygon( [ (0, 0), (40000, 0), (40000, 30000), (20000, 50000), (0, 30000), (0, 0), ], srs=4326, ) print(type(aachen_polygon))

<class 'osgeo.ogr.Geometry'>

Visualize Polygon¶

In [3]:

ax_handle = geokit.core.geom.drawGeoms(
    geoms=aachen_polygon,
    figsize=(4, 4),
)

ax_handle = geokit.core.geom.drawGeoms( geoms=aachen_polygon, figsize=(4, 4), )

Store Polygon as Vector File¶

In [4]:

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

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

Out[4]:

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