我正在尝试使用python ppygis包设置基本的postgis设置。
>>> import psycopg2
>>> import ppygis
>>> connection = psycopg2.connect(database='spre', user='postgres')
>>> cursor = connection.cursor()
>>> cursor.execute('CREATE TABLE test (geometry GEOMETRY)')
>>> cursor.execute('INSERT INTO test VALUES(%s)', (ppygis.Point(1.0, 2.0),))
>>> cursor.execute('SELECT * from test')
>>> point = cursor.fetchone()[0]
>>> print point
0101000000000000000000F03F0000000000000040
>>>
我应该有一个带有单独的X和Y坐标的python对象。像
这样的东西>>> Point(X: 1.0, Y: 2.0)
我做错了什么?
答案 0 :(得分:3)
你没有做错任何事。在the same steps as the PPyGIS basic example之后,我得到了问题(010100 ...)中显示的hex-encoded EWKB,这通常是预期的。也许这适用于旧版本的PPyGIS / Psycopg?今天没有。
该软件包似乎没有正确地将自己注册为PostGIS类型的适配器,因此我的建议是不使用该软件包。此外,您不需要额外的包来从Psycopg2使用PostGIS。
这是读取/写入点的常规方法,没有任何额外的包:
# Assuming PostGIS 2.x, use a typmod
cursor.execute('CREATE TEMP TABLE test (geom geometry(PointZ,4326));')
# Longyearbyen, 78.22°N 15.65°E, altitude 10 m
cursor.execute('''\
INSERT INTO test (geom)
VALUES(ST_SetSRID(ST_MakePoint(%s, %s, %s), 4326));
''', (15.65, 78.22, 10.0))
cursor.execute('''\
SELECT ST_Y(geom) AS latitude, ST_X(geom) AS longitude, ST_Z(geom) AS altitude
FROM test;
''')
print(cursor.fetchone()) # (78.22, 15.65, 10.0)
cursor.execute('SELECT ST_AsText(geom) FROM test;')
print(cursor.fetchone()[0]) # POINT Z (15.65 78.22 10)
cursor.execute('SELECT ST_AsLatLonText(geom) FROM test;')
print(cursor.fetchone()[0]) # 78°13'12.000"N 15°39'0.000"E
如果您希望客户端的几何对象更多地使用实际几何体,请考虑使用Shapely,它可以使用WKB数据进行接口:
from shapely.wkb import loads
cursor.execute('SELECT geom FROM test;')
pt = loads(cursor.fetchone()[0], hex=True)
print(pt) # POINT Z (15.65 78.22 10)