ValueError:零大小的数组，减少操作的最小值，没有标识。

#encoding=utf-8#

import os
import numpy as np
from matplotlib.mlab import griddata
import matplotlib.pyplot as plt

file = open('datalist.txt','r')
a = file.readlines()
file.close()
data = [lines.split('\t') for lines in a]
# read the rainfall data

file = open('lonlat.txt','r')
b = file.readlines()
file.close()
lonlat = [lines.split('\t') for lines in b]
# read the lon, lat data

main_data = []
for d_a in data:
    for l_b in lonlat:
        if d_a[2] == l_b[0]:
            main_data.append(d_a + l_b[1:])
# integrate them

numcols, numrows = len(main_data[0]), len(main_data)

Lon = [m_a[4] for m_a in main_data]
Lat = [m_a[5] for m_a in main_data]
z = [m_a[3] for m_a in main_data] # define the lon, lat and rainfall

xi = np.linspace(22, 23, 102)
yi = np.linspace(113, 115, 102)
xi, yi = np.meshgrid(xi, yi)

x, y, z = Lon, Lat, z
zi = griddata(x, y, z, xi, yi)

fig, ax = plt.subplots(figsize=(10, 5))
im = ax.contour(xi, yi, zi, colors = 'black', linewidth = 0.2)
im = ax.contourf(xi, yi, zi)
ax.scatter(Lon, Lat, c=z, s=50,
           vmin=zi.min(), vmax=zi.max())
fig.colorbar(im)

plt.xticks(())
plt.yticks(())
#encoding=utf-8#

import os
import numpy as np