Erstellen Sie Spatia Lite in Python

Über diesen Artikel

Hinweise zum Erstellen von Spatialite in Python

Ausführungsumgebung

OS Windows7-64bit Python 2.7.10

DLL herunterladen

Laden Sie die Windows-Binärdatei von [hier] [Link-1] herunter. [link-1]:https://www.gaia-gis.it/spatialite-2.3.1/binaries.html

• libspatialite-win-x86-2.3.1.zip
• libiconv-win-x86-1.9.2.zip
• proj-win-x86-4.6.1.zip

DLL-Kopie

Wenn Sie es entpacken, befindet sich in jedem bin-Ordner eine DLL. Kopieren Sie sie daher in einen Ordner, der den Pfad übergibt

• libspatialite-1.dll
• libiconv2.dll
• libproj-0.dll
• Es scheint, dass geos * .dll nicht erforderlich ist, wenn Sie nur eine Datenbank erstellen.

Skript

Ein einfaches Skript zum Erstellen eines Spatialites mit POINT, LINESTRING, POLYGON Auf die erstellte sample.spatialite kann mit GIS-Tools wie QGIS verwiesen werden.

sample.py

# -*- coding: utf-8 -*- 
import sys
import os
import sqlite3

if __name__ == "__main__":
  conn = sqlite3.connect("sample.spatialite")
  if conn:
    print 'connect success'
  else:
    print 'connect failes'
    sys.exit()
  conn.enable_load_extension(True)
  conn.execute('SELECT load_extension("libspatialite-1.dll")')
  conn.execute('SELECT InitSpatialMetaData()')
  conn.execute("INSERT INTO spatial_ref_sys (srid, auth_name, auth_srid, ref_sys_name, proj4text) VALUES (4326, 'epsg', 4326, 'WGS 84', '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs')")

  #POINT
  conn.execute('CREATE TABLE "point" ("OBJECTID" INTEGER PRIMARY KEY AUTOINCREMENT)')
  conn.execute('Select AddGeometryColumn ("point", "Geometry", 4326, "POINT", 2)')
  conn.execute('INSERT INTO point (Geometry) VALUES(GeomFromText("POINT(139.69 35.679)",4326))')

  #LINESTRING
  conn.execute('CREATE TABLE "line" ("OBJECTID" INTEGER PRIMARY KEY AUTOINCREMENT)')
  conn.execute('Select AddGeometryColumn ("line", "Geometry", 4326, "LINESTRING", 2)')
  conn.execute('INSERT INTO line (Geometry) VALUES(GeomFromText("LINESTRING(139.69 35.68, 139.691 35.681, 139.692 35.68)",4326))')

  #POLYGON
  conn.execute('CREATE TABLE "polygon" ("OBJECTID" INTEGER PRIMARY KEY AUTOINCREMENT)')
  conn.execute('Select AddGeometryColumn ("polygon", "Geometry", 4326, "POLYGON", 2)')
  #Exterior CCW - Interior:CW
  conn.execute('INSERT INTO polygon (Geometry) VALUES(GeomFromText("POLYGON((139.69 35.682, 139.69 35.681, 139.691 35.681, 139.691 35.682, 139.69 35.682),(139.6902 35.6818, 139.6908 35.6818, 139.6908 35.6812, 139.6902 35.6812, 139.6902 35.6818))",4326))')

  conn.commit()
  conn.close()

Andere

• OGC:Simple Feature Access - Part 1: Common Architecture http://www.opengeospatial.org/standards/sfa

• Spatialite SQL-Liste https://github.com/azavea/acs-alchemist/blob/master/csharp/Azavea.NijPredictivePolicing.Common/init_spatialite.sql

• Die Geometrie ist nur vom Typ XY Spatialite (2.3.1) scheint XYZ nicht zu unterstützen Wenn definiert, tritt der folgende Fehler auf AddGeometryColumn() error: argument 5 [dimension] current version only accepts dimension=2