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

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

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; }

原图：

标注水平线和亮线的图：

旋转后：

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