在shapefile的字段中查找最大值

Question

我有一个shapefile（mich_co.shp），我试图找到最大人口的县。我的想法是使用max（）函数，这是不可能的。到目前为止，这是我的代码：

from osgeo import ogr
import os

shapefile = "C:/Users/root/Python/mich_co.shp"
driver = ogr.GetDriverByName("ESRI Shapefile")
dataSource = driver.Open(shapefile, 0)
layer = dataSource.GetLayer()

for feature in layer:
    print(feature.GetField("pop"))
layer.ResetReading()

但是上面的代码只打印“pop”字段的所有值，如下所示：

10635.0
9541.0
112039.0
29234.0
23406.0
15477.0
8683.0
58990.0
106935.0
17465.0
156067.0
43868.0
135099.0

我试过了：

print(max(feature.GetField("pop")))

但它返回TypeError：'float'对象不可迭代。为此，我也尝试过：

for feature in range(layer):

它返回TypeError：'Layer'对象不能解释为整数。

任何提示的帮助将非常感激。

谢谢！

Answer 1

max()需要一个可迭代的，例如列表。尝试建立一个列表：

pops = [ feature.GetField("pop") for feature in layer ]
print(max(pops))