如何在qgis中编写统计特征属性的程序？

# 导入必要的模块
from qgis.core import QgsVectorLayer, QgsFeature, QgsField, QgsExpression, QgsFeatureRequest

# 加载矢量图层
layer = QgsVectorLayer('/path/to/your/layer.shp', 'layer_name', 'ogr')

# 检查图层是否加载成功
if not layer.isValid():
    print('图层加载失败！')

# 定义要统计的字段名
field_name = 'field_name'

# 获取字段索引
field_index = layer.fields().indexFromName(field_name)

# 定义统计特征属性
stats = QgsVectorLayerStatistics()
stats.calculate(layer, field_index, QgsStatisticalSummary.All)

# 打印统计结果
print('最大值:', stats.max())
print('最小值:', stats.min())
print('平均值:', stats.mean())
print('总和:', stats.sum())
print('计数:', stats.count())

# 通过表达式获取其他统计特征属性
expression = QgsExpression('median({})'.format(field_name))
expression.prepare(layer.pendingFields())
context = QgsExpressionContext()
context.setFeature(feature)
value = expression.evaluate(context)
print('中位数:', value)

# 通过特征请求获取其他统计特征属性
request = QgsFeatureRequest().setFilterExpression('"{}" > 100'.format(field_name))
count = 0
for feature in layer.getFeatures(request):
    count += 1
print('大于100的要素数量:', count)