我在 QGIS-user forum 上发布了这个问题但这可能太技术性了,因此我也在这里尝试一下。
我的计划是有一个 QGIS 插件,它可以读取并绘制与矢量图层中选定点相对应的时间序列数据(这些点代表不同的测量站,有时我需要快速查看其中每个测量站的时间序列)同一地 block 内的多个站点)。因此,我需要 Python 从 Spatialite 表中读取时间序列数据,然后使用 Matplotlib.plot 绘制它:
from pyspatialite import dbapi2 as sqlite
import numpy as np
import matplotlib.pyplot as plt
import datetime
MyPath = r"C:\DATA\DBTEST\MyTestDb.sqlite"
conn = sqlite.connect(MyPath,detect_types=sqlite.PARSE_DECLTYPES|sqlite.PARSE_COLNAMES)
curs = conn.cursor()
sql = r"""SELECT Date as 'date [datetime]', Level_masl FROM MyTable"""
rs = curs.execute(sql)
recs = rs.fetchall()
My_format = [('date_time', datetime.datetime), ('values', float)]
table = np.array(recs, dtype=My_format)
table2=table.view(np.recarray)
fig = plt.figure()
ax = fig.add_subplot(111)
ax.plot(table2.date_time, table2.values, 'o-')
....
但似乎我没有获得正确的数据类型,因为 matplotlib 不接受我的“date_time”。下面显示一个打印表2,我猜小u表示它是unicode而不是datetime?
print table2
block 引用>recarray([(u'2011-04-20 09:42:00', 703.46000000000004), (u'2011-04-20 09:43:00', 705.35000000000002), ... dtype=[('date_time', '|O4'), ('values', '
我不是程序员,而且绝对是一个完整的 Python 和 Spatialite 初学者,因此我们将非常感谢任何帮助!
/JK
最佳答案
基本上,您正在尝试传递一个字符串。
Matplotlib 有几个方便的函数,用于将日期字符串转换为其(内部)日期时间表示形式。
尝试更多类似的东西:
...
from matplotlib.dates import datestr2num
...
...
time = datestr2num(table2.date_time)
ax.plot_date(time, table2.values, 'o-')
...
plot_date
与 plot
相同,它只是自动格式化 x 轴以具有时间(例如小时、天、月等)标签。
如果您想使用 plot
(或 bar
等),只需调用 ax.plot(time, y, ...)
或 ax.bar
像往常一样,然后调用 ax.xaxis_date()
设置日期标签。
关于Python:使用日期时间格式进行 numpy 重新排列的 SpatiaLite 表?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/7757461/