我有两个相同大小的矢量,一个用于波高,一个用于与测量时间相同的时间段相对应的周期 . 我想知道重复两个特定数据的次数,例如:
Hs = [0.5 1.0 2.3 0.5 0.5]
Tm = [2.0 2.5 2.0 2.0 3.0]
所以你可以看到:
Hs Tm Count
0.5 2.0 2
0.5 2.5 0
0.5 3.0 1
1.0 2.0 0
1.0 2.5 1 ...
我尝试但发生了以下错误,因为我看到整个行和列没有数据,当我看到我的信息值 .
from numpy import *
from matplotlib.pyplot import *
import matplotlib.pyplot as plt
from time import *
clf; cla; close
dat = loadtxt("ecmwf.dat", unpack=True)
HSf = dat[0,:]
HSf = around(HSf,decimals=1)
TMf = dat[1,:]
TMf = around(TMf,decimals=1)
mmat = zeros((31,141))
vhs = linspace(0.0,3.0,31)
vtm = linspace(0.0,14.0,141)
for i in xrange(0, vtm.size):
for k in xrange(0, vhs.size):
if all((k <= vhs.size) & (i <= vtm.size)):
lg1 = (TMf == vtm[i]) & (HSf == vhs[k])
lg2 = sum(lg1)
if lg2>=1:
fg1 = text(i,k, str(lg2),horizontalalignment='center', verticalalignment='center',fontsize=6)
mmat[k,i] = lg2