付丽叶变换图求水平线的斜率

付丽叶变换后的图，主要特点是中间有两条十字线，一横一竖，可能还有其他的线，研究不深，暂不做评论。十字线中间是亮线。如何获取该亮线的方程呢？假定中心点在图的中心，根据直线方程：

y=kx+b （1）

直线过点P(w/2,h/2),代入公式（1）得：

b=(h-kw)/2        （2）

因此直线方程可表示为：

y = k(x - w/2) + h/2        （3）

绕点P旋转，旋转步长为0.1度，最大倾角45度，因此总步长为450，求每个步长对应的直线，统计该线上的点的灰度累加值。越白的像素，灰度值越大，累加值最大的线为所求直线。

//倾斜校正
void CBmptestApp::OnLeanAdjust(){
    char* sMax = new char[100];
    
    int w = objCDib->GetWidth();
    int h = objCDib->GetHeight();
    DWORD size = w * h;

    BYTE* img = objCDib->GetData();
    BYTE* ftImg = new BYTE[size];
    memcpy(ftImg,img,size);

    CDC* pDC = m_pMainWnd->GetDC();
    int s = 27; /*工具条高度*/

    double pi = 3.1415926;
    int count = 450; //每次旋转0.1度
    long lines[450];
    
    for(int i=0;i<450;i++){
        lines[i] = 0;
    }

    AfxMessageBox("here0"); 

    for(i=0;i<count;i++){
        for(int j=0;j<w/2;j++){
            int x = w/2 + j;
            int y = tan(0.1*i*pi/180) * (x-w/2) + h/2;
            lines[i] += ftImg[w*y +x];
        }
    }

    AfxMessageBox("here1"); 

    int maxV = 0;
    int maxIdx = 0;
    for(i=0;i<count;i++){
        if(lines[i] > maxV){
            maxV = lines[i];
            maxIdx = i;
        }
    }


    sprintf(sMax,"maxV=%d,maxIdx=%d", maxV, maxIdx);
    AfxMessageBox(sMax); 
    delete[] sMax;


    //标识倾斜线
    int h0 = h/2 + w/2 * tan(0.1*maxIdx*pi/180) + s;
    CPen* pPen = new CPen;
    pPen->CreatePen(PS_SOLID,2,RGB(255,0,0));            
    // 选中新画笔
    CGdiObject* pOldPen = pDC->SelectObject(pPen);

    pDC->MoveTo(0,0);
    pDC->LineTo(w,0);

    pDC->MoveTo(0,s);
    pDC->LineTo(w,s);

    pDC->MoveTo(w/2,h/2+s);
    pDC->LineTo(w,h0);

    pDC->MoveTo(w/2,h/2+s);
    pDC->LineTo(w,h/2+s);

    // 选回以前的画笔
    pDC->SelectObject(pOldPen);
    delete pPen;

    //旋转图片
    BYTE* newImg =new BYTE[size];
    rotateImage(ftImg, w, h, newImg, w, h, -0.1*maxIdx);
    int rv = ShowGrayImg(pDC,3,s,newImg,w,h,NULL,NULL);
    delete[] ftImg;
}

原图：

标注水平线和亮线的图：

旋转后：

这只是一个简单的测试，仅仅根据图像的特征来求的直线，实际上是说：最亮的一条线怎么得到？