数组相加出错

有事秦天

@MeteoInfo 王老师您好，我在处理lisotd闪电数据时遇到问题不知道怎么解决？对区域范围一年的闪电数叠加，但因为有个别月份一部分数为NaN，造成这个点上的数相加后数值为0，在图上显示有一片变成空的了，请问这种情况该怎么处理？
LISOTD_HRMC_V2.3.2014.hdf (14.02 MB, 下载次数: 11)

2017-5-18 17:22 上传

点击文件名下载附件

f = addfile('E:/lightning/LISOTD_HRMC_V2.3.2014.hdf')
n = 0
i = 0
for i in range(0,11):
    lon = '101:109'
    lat = '25:34'
    flash = f['HRMC_COM_FR'][lat,lon,i]
    n = flash + n
    i = i + 1
lon = f['Longitude'][:]
lat = f['Latitude'][:]
c = n.ave()
print(c)
#Plot
axesm()
mlayer = shaperead('D:/MeteoInfo/map/province_2004.shp')
geoshow(mlayer)
player = shaperead('D:/MeteoInfo/map/rivers.shp')
geoshow(player)
#levs = [0,0.005,0.007,0.010,0.015,0.020,0.025,0.030,0.035,0.040,0.050,0.060,0.070,0.080,0.090,0.100]
#cols = [(255,250,250),(220,220,220),(255,228,196),(190,190,190),(119,136,153),(105,105,105),(100,149,237),(132,112,255),(0,0,255),(0,225,225),(32,178,170),(0,100,0),(255,215,0),(255,255,0),(240,128,128),(255,69,0),(255,0,0)]
layer = contourfm(n)
colorbar(layer)
xlim(102, 108)
ylim(26, 33)

MeteoInfo

可以在相加之前将nan转为0：

flash = f['HRMC_COM_FR'][lat,lon,i]
flash[flash==nan] = 0

有事秦天

MeteoInfo 发表于 2017-5-18 22:46
可以在相加之前将nan转为0：

flash = f['HRMC_COM_FR'][lat,lon,i]

谢谢王老师，图画出来了