高分求等值线算法

在一幅灰度图中，根据制定灰度值绘制等值线
分不够我还会加的
效率越高分数越高

18 个解决方案

#1

我多次发贴未成功!

#2

成功了! 再发发看:

按灰度级大小,由小到大,用轮廓跟踪法,将小于,和大于等于该灰度级的区域提取边界线即可.

#3

给我一幅实际的图像,我编个程序试一试,把结果发给您,看行不行.
我的Email: zzwu@citiz.net

#4

关注

#5

因为我需要对一个热成像仪传来的图像绘制等值线，而图像就是一幅灰度图，所以你只要随便找一幅320*240的256级灰度图就可以了
只是有可能会需要绘制不止一条等值线

#6

等值线通常都不会只有一条的,像显示地形的地图那样,除非图像本来只有2种颜色.

按灰度等级,用轮廓跟踪法来画等值线时,一种灰度就可产生一组等值线,一组可以是0根线,1根线,或很多根线. 所以256灰级的图像的等值线特殊情况说不定会有成千上万条的通常都是封闭的等值线.

但担心等值线太多也不必,因只要把灰度的等级数少分一些就可以了.

#7

                    333333333
                   3         333     333333333
                 33  4444444    33333          333
                3   4       444       44444444    33
               3  4    555     4     4        44    33
              3  4    5   5    4     4   55555  4     3
              3  4    5   5   4     4   5 666 5  4    3
              3   4    555    4     4   5 6  6 5  4   3
               3   4         4  33   4  5 666 5  4   3
               3    444444444  3  3   4  55555  4    3
                3             3    3    44444444    3
                  333333333333      3              3
                                     33333333333333

#8

这不是好办法,因为要一次一次地扫描,等级少还可以,多了就慢了.

一次扫描,就提取所有等高线,也是能办得到的,我已想出了办法.

#9

zzwu(未名) :
办法是。。。？？

#10

一次扫描,就提取所有等高线的办法，告诉我吧

#11

An_Apple (我爱读书) ：

一次扫描提取所有等高线的办法，实际就是一次扫描提取图像所有边界线的办法一样，仍然是轮廓跟踪法。只是把所有有不同高度（亦即灰度值）的连通区域作为区域就行了。

所以，需要弄清轮廓跟踪法的基本原理，以及其应用于二值图像时的基本算法。

下面是轮廓跟踪法的原理及一个Ｃ程序，供您参考。

＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊＊

轮廓跟踪法顾名思义就是通过顺序找出边缘点来跟踪出边界。

一个简单二值图象闭合边界的轮廓跟踪算法很简单：首先按从上到下，从左到右的顺序搜索，找到的第一个黑点一定是最左上方的边界点，记为A。它的右，右下，下，左下四个邻点中至少有一个是边界点，记为B。从开始B找起，按右，右下，下，左下，左，左上，上，右上的顺序找相邻点中的边界点C。如果C就是A点，则表明已经转了一圈，程序结束；否则从C点继续找，直到找到A为止。判断是不是边界点很容易：如果它的上下左右四个邻居都是黑点则不是边界点，否则是边界点。源程序如下，其中函数IsContourP用来判断某点是不是边界点。

BOOL Contour(HWND hWnd)

{

       DWORD                             OffBits，BufSize;

LPBITMAPINFOHEADER    lpImgData;

       LPSTR                   lpPtr;

       HLOCAL                  hTempImgData;

       LPBITMAPINFOHEADER    lpTempImgData;

       LPSTR                   lpTempPtr;

       HDC                      hDc;

       HFILE                    hf;

       LONG                    x,y;

       POINT                   StartP,CurP;

       BOOL                     found;

       int                        i;

int       direct[8][2]={{1,0},{1,1},{0,1},{-1,1},{-1,0},{-1,-1},{0,-1},{1,-1}};

//我们处理的实际上是256级灰度图，不过只用到了0和255两种颜色。

       if( NumColors!=256){

MessageBox(hWnd,"Must be a mono bitmap with grayscale palette!",

"Error Message",MB_OK|MB_ICONEXCLAMATION);

return FALSE;

}

//到位图数据的偏移值

       OffBits=bf.bfOffBits-sizeof(BITMAPFILEHEADER);

       //缓冲区大小

BufSize=OffBits+bi.biHeight*LineBytes;

//为新图缓冲区分配内存

       if((hTempImgData=LocalAlloc(LHND,BufSize))==NULL)

     {

       MessageBox(hWnd,"Error alloc memory!","Error Message",

MB_OK|MB_ICONEXCLAMATION);

return FALSE;

}

      lpImgData=(LPBITMAPINFOHEADER)GlobalLock(hImgData);

       lpTempImgData=(LPBITMAPINFOHEADER)LocalLock(hTempImgData);

       //新图缓冲区初始化为255

       memset(lpTempImgData,(BYTE)255,BufSize);

       //拷贝头信息

       memcpy(lpTempImgData,lpImgData,OffBits);

       //找到标志置为假

       found=FALSE;

       for (y=0;y<bi.biHeight && !found; y++){

              lpPtr=(char *)lpImgData+(BufSize-LineBytes-y*LineBytes);

              for (x=0;x<bi.biWidth && !found; x++)

                     if (*(lpPtr++) ==0) found=TRUE;

//找到了最左上的黑点，一定是个边界点

       }

       if(found){ //如果找到了，才做处理

//从循环退出时，x，y坐标都做了加1的操作。在这里把它们减1，得到

//起始点坐标StartP

              StartP.x=x-1;

              StartP.y=y-1;

              lpTempPtr=(char*)lpTempImgData+

(BufSize-LineBytes-StartP.y*LineBytes)+StartP.x;

              *lpTempPtr=(unsigned char)0; //起始点涂黑

              //右邻点

            CurP.x=StartP.x+1;

              CurP.y=StartP.y;

              lpPtr=(char *)lpImgData+(BufSize-LineBytes-CurP.y*LineBytes)+CurP.x;

              if(*lpPtr!=0){ //若右邻点为白，则找右下邻点

                   CurP.x=StartP.x+1;

                     CurP.y=StartP.y+1;

                     lpPtr=(char*)lpImgData+

(BufSize-LineBytes-CurP.y*LineBytes)+CurP.x;

                     if(*lpPtr!=0){ //若仍为白，则找下邻点

                          CurP.x=StartP.x;

                            CurP.y=StartP.y+1;

                     }

                     else{ //若仍为白，则找左下邻点

                          CurP.x=StartP.x-1;

                            CurP.y=StartP.y+1;

                     }

              }

              while (! ( (CurP.x==StartP.x) &&(CurP.y==StartP.y))){ //知道找到起始点，

//循环才结束

                     lpTempPtr=(char*)lpTempImgData+

(BufSize-LineBytes-CurP.y*LineBytes)+CurP.x;

                     *lpTempPtr=(unsigned char)0;

                     for(i=0;i<8;i++){

//按右，右上，上，左上，左，左下，下，右下的顺序找相邻点

//direct[i]中存放的是该方向x,y的偏移值

                            x=CurP.x+direct[i][0];

                            y=CurP.y+direct[i][1];

              //lpPtr指向原图数据，lpTempPtr指向新图数据

                            lpTempPtr=(char*)lpTempImgData+

(BufSize-LineBytes-y*LineBytes)+x;

                            lpPtr=(char *)lpImgData+(BufSize-LineBytes-y*LineBytes)+x;

                            if(((*lpPtr==0)&&(*lpTempPtr!=0))||

((x==StartP.x)&&(y==StartP.y)))

                            //原图中为黑点，且新图中为白点(表示还没搜索过)时才处理

                            //另一种可能是找到了起始点

                                   if(IsContourP(x,y,lpPtr)){ //若是个边界点

                                          //记住当前点的位置

                            CurP.x=x;

                                          CurP.y=y;

                                          break;

                                   }

                     }

              }

       }

    if(hBitmap!=NULL)

           DeleteObject(hBitmap);

       hDc=GetDC(hWnd);

       //创立一个新的位图

       hBitmap=CreateDIBitmap(hDc,(LPBITMAPINFOHEADER)lpTempImgData,

(LONG)CBM_INIT,

(LPSTR)lpTempImgData+

sizeof(BITMAPINFOHEADER)+

NumColors*sizeof(RGBQUAD),

                                  (LPBITMAPINFO)lpTempImgData,

DIB_RGB_COLORS);

       hf=_lcreat("c:\\contour.bmp",0);

       _lwrite(hf,(LPSTR)&bf,sizeof(BITMAPFILEHEADER));

       _lwrite(hf,(LPSTR)lpTempImgData,BufSize);

       _lclose(hf);

       //释放内存和资源

      ReleaseDC(hWnd,hDc);

       LocalUnlock(hTempImgData);

       LocalFree(hTempImgData);

       GlobalUnlock(hImgData);

       return TRUE;

}

//判断某点是不是边界点，参数x,y 为该点坐标，lpPtr为指向原数据的指针

BOOL IsContourP(LONG x,LONG y, char *lpPtr)

{

       int    num,n,w,e,s;

       n=(unsigned char)*(lpPtr+LineBytes); //上邻点

       w=(unsigned char)*(lpPtr-1); //左邻点

       e=(unsigned char)*(lpPtr+1); //右邻点

       s=(unsigned char)*(lpPtr-LineBytes); //下邻点

       num=n+w+e+s;

       if(num==0) //全是黑点，说明是个内部点而不是边界点

              return FALSE;

       return TRUE;

}

#12

to zzwu
这是数字图像处理编程入门里的轮廓跟踪，不能实现一次扫描画多个值的等值线，不过还是给你一些分吧，到下边去领
http://www.csdn.net/expert/topic/830/830441.xml?temp=.9140131
谢谢

还有没有高手能搞定啊，救救我吧

#13

1。您说的对，这确实不是实现一次扫描画多个值的等值线的算法。这只是为弄清轮廓跟踪法的基本原理所需要的求单根边界的轮廓跟踪。
2。如果轮廓跟踪原理您已熟悉，下面就可以来考察我的想法了。
3。首先要说明，下面谈的是我看了您的问题后才生成的新的想法，不是我已做工作的归纳。我确实编过类似问题的程序，但不是和您所要的完全相同。所以，所说的想法完全可能会不全面，如有问题，希望您能毫不客气地批评指正，让我们一起来完善算法的思路。
4。在谈我的想法之前，还想了解一下：等高线只要画出来就可以了吗？还是要把它们作为链码保存起来？

#14

画出来就可以了
还有千万别再用您这个字眼了，这辈子还没几个人这么称呼过我，别吓坏我了

#15

to An_Apple(我爱读书)
你是否已经得到了一次扫描得出全部扫描线的算法？还需要兄弟帮忙吗？

#16

还没有啊!!!
需要帮忙啊!!

#17

而且要高效的，因为需要在实时采集的图像上处理

#18

#1