小编ham*_*mmu的帖子

我想创建一个矩阵子图并显示每个BMP文件,从一个目录,在不同的子图中,但我找不到适合我的问题的解决方案,有人可以帮助我吗？

这是我的代码:

import os, sys
from PIL import Image
import matplotlib.pyplot as plt
from glob import glob

bmps = glob('*trace*.bmp')

fig, axes = plt.subplots(3, 3)

for arch in bmps:
    i = Image.open(arch)
    iar = np.array(i)
    for i in range(3):
        for j in range(3):
            axes[i, j].plot(iar)
            plt.subplots_adjust(wspace=0, hspace=0)
plt.show()

执行后我遇到以下错误:

我有一个具有以下结果的 DataFrame：

我想将这些坐标列从 WGS84 转换为经度和纬度，最后在我的数据框中添加这些新列：

对于转换，我使用以下代码，但我认为应该有一种更好的方法，无需将坐标列转换为列表并创建一个新的 DataFrame。

import pyproj as pp
from mpl_toolkits.basemap import Basemap
import pandas as pd

cx =dfb.COORDENADA_X.tolist()
cy =dfb.COORDENADA_Y.tolist()

utm15_wgs84 = pp.Proj(init='epsg:32615')

for ix, iy in zip(cx, cy):
    lon, lat = utm15_wgs84(ix, iy, inverse=True)
    print(lon, lat)

这样做有什么建议吗？