测试 halcon measure_projection 算子

期望结果完全相同，但是下面的测试结果和halcon的差值如下:

[0.132838, 0.231991, 0.265157, 0.296903, 0.0998573, 0.165907, 0.230686, 0.130266, 0.0977104, 0.197109, 0.198173, 0.197086, 0.190943, 0.177665, 0.163521, 0.146541, 0.161362, 0.166666, 0.22818, 0.186114, 0.196333, 0.255945, 0.14496, 0.205199, 0.262202, 0.201859, 0.173616, 0.173016, 0.173485, 0.173227, 0.171494, 0.257541, 0.171186, 0.141769, 0.0852231, 0.257203, 0.198859, 0.170551, 0.25673, 0.171491, 0.170175, 0.169464, 0.198114, 0.169922, 0.169473, 0.1695, 0.168997, 0.139408, 0.195358, 0.193962, 0.221091, 0.19246, 0.246518, 0.138458, 0.194492, 0.193142, 0.217299, 0.188277, 0.16121, 0.106746, 0.106695, 0.104717, 0.154996, 0.17956, 0.199884, 0.174948, 0.174217, 0.125464, 0.234857, 0.197361, 0.239242, 0.217724, 0.220922, 0.160648, 0.19362, 0.1959, 0.196964, 0.197266, 0.229685, 0.193291, 0.191924, 0.191487, 0.288308, 0.191503, 0.223271, 0.285233, 0.0954785, 0.157538, 0.187701, 0.284095, 0.190475, 0.18948, 0.188329, 0.187838, 0.186164, 0.217482, 0.277055, 0.215032, 0.15109, 0.27215, 0.211512, 0.208731, 0.286509, 0.262083, 0.31329, 0.315317, 0.358043, 0.396789, 0.403176, 0.397609, 0.39305, 0.392336, 0.19589, 0.260681, 0.459158, 0.330526, 0.19935, 0.587357, 0.518498, 0.44198, 0.24967]

结果并不小,不知道问题出在哪里?

get_grayval_interpolated 和自己实现的双线性差值的结果是完全相同的

dev_close_window()
dev_open_window (0, 0, 512, 512, 'black', WindowHandle)
dev_set_window (WindowHandle)
read_image (Image, 'fabrik')
get_image_size (Image, Width, Height)
CenterRow:=Height/2-50
CenterCol:=Width/2
degPhi:=45
Phi:=rad(degPhi)
Length1:=60
Length1:=int(floor(Length1))
Length2:=1
gen_measure_rectangle2 (CenterRow, CenterCol, -Phi, Length1, Length2, Width, Height, 'bilinear', MeasureHandle)
measure_projection (Image, MeasureHandle, GrayValues)
disp_rectangle2(WindowHandle,CenterRow, CenterCol, -Phi, Length1, Length2)

RowAry:=[]
ColAry:=[]
RowAry[Length1]:=CenterRow
ColAry[Length1]:=CenterCol

RowOffsetAry[Length2]:=0
ColOffsetAry[Length2]:=0
clockWise:=true
if(clockWise)
  Phi2:=rad(degPhi-90)  
else
   Phi2:=rad(degPhi+90)    
endif
RowOffsetStep:=sin(Phi2)
ColOffsetStep:=cos(Phi2)
for Index := 1 to Length2 by 1
    RowOffsetAry[Length2+Index]:=RowOffsetStep*Index
    ColOffsetAry[Length2+Index]:=ColOffsetStep*Index
    RowOffsetAry[Length2-Index]:=-RowOffsetAry[Length2+Index]
    ColOffsetAry[Length2-Index]:=-ColOffsetAry[Length2+Index]
endfor
gen_image_const (Image1, 'real', Length1*2+1, Length2*2+1)
get_image_size (Image1, Width1, Height1)
RowOffsetStep:=sin(Phi)
ColOffsetStep:=cos(Phi)
for Index := 1 to Length1 by 1
    RowAry[Length1+Index]:=CenterRow+Index*RowOffsetStep
    ColAry[Length1+Index]:=CenterCol+Index*ColOffsetStep
    RowAry[Length1-Index]:=CenterRow-Index*RowOffsetStep
    ColAry[Length1-Index]:=CenterCol-Index*ColOffsetStep
endfor
GrayValues2:=[]
copy_image(Image, DupImage)

for Index := 0 to Width1-1 by 1
   
    get_grayval_interpolated (Image, RowAry[Index], ColAry[Index], 'bilinear', Grayval)
    
    left:= int(floor(ColAry[Index]))
    right:=left+1
    top:= int(floor(RowAry[Index]))
    bottom:=top+1
    a:=ColAry[Index]-left
    b:=RowAry[Index]-top
    get_grayval(Image,top,left,topLeftGrayVal)
    get_grayval(Image,top,right,topRightGrayVal)
    get_grayval(Image,bottom,left,bottomLeftGrayVal)
    get_grayval(Image,bottom,right,bottomRightGrayVal)
    Grayval:=(1 - a)*(1 - b)*topLeftGrayVal+ a*(1 - b)*topRightGrayVal +(1 - a)*b*bottomLeftGrayVal + a*b*bottomRightGrayVal
    
    
    set_grayval (Image1, Length2, Index, Grayval)
    set_grayval(DupImage,int(RowAry[Index]),int(ColAry[Index]),0)
    AVGGrayVal:=Grayval
    for Index2 := 1 to Length2 by 1
        Row:=RowAry[Index]+RowOffsetAry[Length2+Index2]
        Col:=ColAry[Index]+ColOffsetAry[Length2+Index2]
        
        get_grayval_interpolated (Image, Row, Col, 'bilinear', Grayval)
        
        
        left:= int(floor(Col))
        right:=left+1
        top:= int(floor(Row))
        bottom:=top+1
        a:=Col-left
        b:=Row-top
        get_grayval(Image,top,left,topLeftGrayVal)
        get_grayval(Image,top,right,topRightGrayVal)
        get_grayval(Image,bottom,left,bottomLeftGrayVal)
        get_grayval(Image,bottom,right,bottomRightGrayVal)
        Grayval:=(1 - a)*(1 - b)*topLeftGrayVal+ a*(1 - b)*topRightGrayVal +(1 - a)*b*bottomLeftGrayVal + a*b*bottomRightGrayVal
        
        
        
        set_grayval (Image1, Length2+Index2, Index, Grayval)
        set_grayval(DupImage,int(Row),int(Col),0)
        AVGGrayVal:=AVGGrayVal+Grayval
        Row:=RowAry[Index]+RowOffsetAry[Length2-Index2]
        Col:=ColAry[Index]+ColOffsetAry[Length2-Index2]
        get_grayval_interpolated (Image, Row, Col, 'bilinear', Grayval)
        
        
        
        left:= int(floor(Col))
        right:=left+1
        top:= int(floor(Row))
        bottom:=top+1
        a:=Col-left
        b:=Row-top
        get_grayval(Image,top,left,topLeftGrayVal)
        get_grayval(Image,top,right,topRightGrayVal)
        get_grayval(Image,bottom,left,bottomLeftGrayVal)
        get_grayval(Image,bottom,right,bottomRightGrayVal)
        Grayval:=(1 - a)*(1 - b)*topLeftGrayVal+ a*(1 - b)*topRightGrayVal +(1 - a)*b*bottomLeftGrayVal + a*b*bottomRightGrayVal
        
        
        
        
        set_grayval (Image1, Length2-Index2, Index, Grayval)
        set_grayval(DupImage,int(Row),int(Col),0)
        AVGGrayVal:=AVGGrayVal+Grayval
    endfor
    AVGGrayVal:=AVGGrayVal/Height1
    GrayValues2[Index]:=AVGGrayVal
endfor
for Index := 0 to |GrayValues2|-1  by 1
    GrayValuesDidff[Index]:=GrayValues2[Index]-GrayValues[Index]
endfor