我正在使用LandScan Raster数据,并尝试将数据放入csv文件中。我这样做是为了对数据进行统计分析。我需要能够获得纬度,经度和人口。将来,我将把其他栅格格式的数据添加到CSV文件中,因此每个纬度和经度都必须与总体相关联。
最初,我尝试使用我感兴趣的坐标创建一个新的数据numpy对象:
dataBand = dataBand[minLong:maxLong, minLat:maxLat]
但是,一旦执行此操作,纬度和经度将不再对应于数组索引。因此,我认为必须进行循环,但这显然会花费很长时间:
for latitudeIndex in range(minLat, maxLat):
for longitudeIndex in range(minLong, maxLong):
x,y = data.xy(longitudeIndex, latitudeIndex)
dataDF = dataDF.append({'long':x, 'lat':y, 'population':dataBand[longitudeIndex, latitudeIndex]}, ignore_index = True)```
是否有任何方法可以优化此代码,或者有更好的方法吗?
谢谢!