我有两个大小相同的矢量,一个是波高,一个是周期,分别对应于进行测量的同一时间点。我想知道两个特定数据重复多少次,例如:
Hs = [0.5 1.0 2.3 0.5 0.5]
Tm = [2.0 2.5 2.0 2.0 3.0]
这样您就可以看到:
Hs Tm计数
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
我建议使用Counter
来计算你的双。
from collections import Counter
Hs = [0.5, 1.0, 2.3, 0.5, 0.5]
Tm = [2.0, 2.5, 2.0, 2.0, 3.0]
occurrences = Counter(zip(Hs, Tm))
for h in sorted(set(Hs)):
for t in sorted(set(Tm)):
print h, t, occurrences[(h,t)]
结果:
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
1.0 3.0 0
2.3 2.0 1
2.3 2.5 0
2.3 3.0 0
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句