多条水平线和垂直线未显示

Question

在这里，我正在尝试使用opencv，霍夫变换来查找水的高度。该程序包括Canny边缘检测，背景减法，霍夫变换和高度估计。然后使用高度绘制图表，但这只会检测垂直线。没有显示多条水平和垂直线。代码有问题吗？

程序输出：

using Distances
using RCall
using Distributions
using BSON: @save, @load
using StatsBase
using LinearAlgebra

R"simul<-function(m){
  comb<-expand.grid(c(0.01,0.2,0.4),
                c(sample(2:7,1),sample(8:12,1),sample(13:20,1)),
                c(sample(2:5,1),sample(6:10,1),sample(11:20,1)),
                c(150,500,1500))
gener<-function(i){
maxoverlap<-comb[i,1]
nbvar<-comb[i,2]
nbclass<-comb[i,3]
propmix<-runif(1,0.001,1/nbclass)
Q<-MixSim(MaxOmega = maxoverlap, K = nbclass, p = nbvar,PiLow = propmix,resN = 1000)
A <- simdataset(n = comb[i,4], Pi = Q$Pi, Mu = Q$Mu, S = Q$S)
results<-list(Q,A)
return(results)
}
donnees<-sapply(1:nrow(comb),gener)
}
library(MixSim)
donneesimul=simul(1)"
@rget donneesimul

function pointsdpp(t)
 datasim=donneesimul[2,t][:X]

 Eucldist=pairwise(Euclidean(),transpose(datasim))
 D=maximum(Eucldist.^2)
 sigma2hat=mean(((Eucldist.^2)./D)[tril!(trues(size((Eucldist.^2)./D)),-1)])
 L=exp.(-(Eucldist.^2/D)/(2*sigma2hat))

 eigenv=eigvals(L)
 prob=eigenv./(eigenv.+1)
 eigenvectors=eigvecs(L)
function sampledpp(m)
u=rand(size(L,1))
V=eigenvectors[:,findall(u.<=prob)]
k=size(V,2)
Y=zeros(Int64,k)
for i=k:-1:1
P=sum(V.^2,dims=2)
Pri=P / sum(P)
Cumpri=cumsum(Pri,dims=1)
u=rand()
Y[i]=findfirst(u.<=Cumpri)[1]
if i==1 break end
j=findfirst(V[Y[i],:].!=0)
Vj=V[:,j]
V=V[:,deleteat!(collect(1:1:size(V,2)),j)]
V=V-repeat(Vj,1,size(V,2)).*repeat(transpose(V[Y[i],:]/Vj[Y[i]]),size(V,1))

for a = 1:i-1
    for b = 1:a-1
      V[:,a] = V[:,a] - transpose(V[:,a])*V[:,b]*V[:,b]
    end
    V[:,a] = V[:,a] / norm(V[:,a])
end
end
Y=sort(Y)
return(Y)
end

m=collect(1:1000)
sampleY_repet=map(sampledpp,m)
end
w=collect(1:1:81)
echantdpp=map(pointsdpp,w)
@save "echantdppdatasim1.bson" echantdpp

Answer 1

您只画一条线，因为您的for循环行仅循环单行。

if lines is not None:
    for line in lines: #for each line...not just one of them
        for x1,y1,x2,y2 in line:
            dot1 = (x1,y1)
            dot2 = (x2,y2)
            cv2.line(frame, dot1, dot2, (255,0,0), 3)
            length = y1 - y2
            print(length)
            height.append(length)