图像分割

阈值分割方法一般可以分为全局阈值法和局部阈值法。全局阈值法指利用全

局信息(例如整幅图像的灰度直方图)对整幅图像求出最优分割阈值,可以是单阈

值,也可以是多阈值:局部阈值法是把原始的整幅图像分为若干小的子图像,再

对每个子图像应用全局阈值法分别求出最优分割阈值。阈值法的最大特点是计算

简单,在重视运算效率的场合得到了广泛的应用。

imagesc是软件Matlab中的一个函数。imagesc(A) 将矩阵A中的元素数值按大小转化为不同颜色,并在座标轴对应位置处以这种颜色染色imagesc(x,y,A) x,y决定座标范围,x,y应是两个二维向量,即x=[x1 x2],y=[y1 y2],matlab会在[x1,x2]*[y1,y2]的范围内染色。 如果xy超过两维,则座标范围为[x(1),x(end)]*[y(1),y(end)]

B = sum(A)返回数组中各维的元素和。如果A是一个向量(即一个n1列的矩阵),sum(A)返回这个向量中所有元素的和;如果A是一个矩阵,sum(A)A的各列看做一个向量,并返回一个行向量(即一个1n列的矩阵),这个行向量的第n个元素是A中第n列元素的和;如果A是一个多维数组,sum(A)仅仅计算A中第一个非奇异维,并把它看成一个向量,计算后返回一个行向量组。B = sum(A,dim) 只对A中第dim维的元素进行计算。如果dim1,计算各列的元素之和;如果dim2,则计算A中各行元素之和。

最速下降法

下降法亦称极小化方法,是一类重要的迭代法。这类方法将方程组求解问题转化为求泛函极小问题。

设给出方程组 F(x)=0,其中 https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D122/sign=6e6c4ce5a3d3fd1f3209a638024c25ce/0824ab18972bd407d91d0c3670899e510eb30927.jpg ,令https://gss2.bdstatic.com/-fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D243/sign=c0168098fd246b607f0eb570d8f91a35/0df431adcbef76093529475425dda3cc7dd99edf.jpg F(x)=0 的充分必要条件是φ(x)=0. φ(x) 的极小点 x* 使φ(x*)=0,则 x* 也是方程组 F(x)=0 的解。只要构造迭代序列 {xk} 使https://gss0.bdstatic.com/94o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D125/sign=da2de00f10d5ad6eaef960e8b4ca39a3/42a98226cffc1e1782b457fa4190f603728de9f6.jpg且满足https://gss2.bdstatic.com/-fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D107/sign=3dfb12f89945d688a702b6a493c37dab/aa64034f78f0f73651c597ca0155b319eac413c6.jpg就可求得方程组的足够好的近似解。具体做法是选初始近似 x0,沿一个使φ(x) 下降的方向 p0,令 https://gss0.bdstatic.com/94o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D136/sign=0b1d5f85c18065387feaa010a1dca115/fcfaaf51f3deb48f145a7b15fb1f3a292cf578ce.jpg 然后选步长因子 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D45/sign=44014fbfa6345982c18ae4970ef4dc88/c9fcc3cec3fdfc0341993873df3f8794a5c2261c.jpg ,使 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D154/sign=2256ee123201213fcb334ad960e636f8/b2de9c82d158ccbfc3c76ff112d8bc3eb0354186.jpg ,得 https://gss2.bdstatic.com/-fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D104/sign=8e5008e28d18367aa9897bdd1a728b68/b3fb43166d224f4a1ea5b6d802f790529822d139.jpg 一般情况是从 xk 出发,沿φ(x) 下降方向 pk 求出https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D214/sign=d7ee686fafc27d1ea1263cc52fd7adaf/0b55b319ebc4b74540786dddc4fc1e178b821508.jpg https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D124/sign=3d00ed6e017b020808c93be356d8f25f/eaf81a4c510fd9f9842c90a22e2dd42a2934a4c6.jpg 直到

 https://gss1.bdstatic.com/9vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D71/sign=baf36c79da58ccbf1fbcb73b18d80154/1f178a82b9014a90ff2e8376a2773912b21beeed.jpg 为止,xm 就作为方程组解 x* 的近似,上述算法中也可选 uk 使https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D233/sign=4f956015fb1f3a295ec8d2cdaa27bce3/71cf3bc79f3df8dc12c14e5ac611728b46102879.jpg这是一个求一维的极小问题。以上算法即为下降法。如果 pk 选择不同,就可得到不同的下降法,特别地,若选 p φ 的负梯度,即 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D105/sign=9c13b6528926cffc6d2abbb28c004a7d/10dfa9ec8a136327498602f89a8fa0ec09fac796.jpg 则得梯度算法https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D159/sign=6a30e9cbe224b899da3c7d3d57071d59/0df3d7ca7bcb0a46b01748f26063f6246a60afc0.jpg其中https://gss1.bdstatic.com/9vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D232/sign=c26b0d09b9119313c343f8b357390c10/fcfaaf51f3deb48f1d7c6215fb1f3a292cf578e0.jpg此算法也成最速下降法,此法的优点是计算量少,程序简单,但收敛慢。在下降法中可去下降方向 pk 为牛顿方向,即https://gss1.bdstatic.com/-vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D155/sign=126b0525bf8f8c54e7d3c12a0f282dee/c8177f3e6709c93dcdf5cd3c943df8dcd00054f2.jpg特别地,当 https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D42/sign=a98fa300b9fb43161e1f7b7822a4938d/b58f8c5494eef01faecbdb0eebfe9925bd317d1a.jpg 时,就得到牛顿法,此外还可取沿座标方向下降的方法,实际上就是一步的 SOR 牛顿法。

共轭梯度法

另一类较重要的下降法为共轭梯度法。共轭梯度法是最简单的下降法,早在 1847 年就由法国数学家、力学家柯西 (Cauchy,A.-L.)提出,以后坦普尔 (Temple,G.)、柯里 (Curry,H.) 等人也进行过研究并证明了方法的收敛性。20世纪5060年代,又有不少学者对下降法做了很多研究,提出不少具体算法并建立了收敛性理论,使这类算法在解非线性方程组和最优化计算中得到广泛应用。 [1] 

参考资料

  • 1.  《数学辞海》总编辑委员会.《数学辞海》第3卷.南京:东南大学出版社,2002.8

 

符号距离函数sign distance function),简称SDF,又可以称为定向距离函数(oriented distance function),在空间中的一个有限区域上确定一个点到区域边界的距离并同时对距离的符号进行定义:点在区域边界内部为正,外部为负,位于边界上时为0

在欧几里得空间属性

如果 https://gss0.bdstatic.com/94o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=70d69ce2a718972ba73a04cee7cd4a3b/9922720e0cf3d7cadb0b9f80fe1fbe096a63a943.jpg 是具有分段平滑边界的欧几里德空间 https://gss1.bdstatic.com/9vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D18/sign=686a2771dca20cf44290fad77409d9e7/d53f8794a4c27d1ee3f39b3217d5ad6edcc4386a.jpg 一个子集,那么带符号的距离函数几乎在任何地方都是可微的,并且它的梯度满足方程

https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D58/sign=95e16964c111728b342d8c2acafc64df/1c950a7b02087bf44e30b11cfed3572c10dfcf63.jpg

如果 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=ac1206297ac6a7efbd26ac22fcfa3ef8/c995d143ad4bd11338f314c456afa40f4bfb050d.jpg 的边界对于 https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D36/sign=4245ee4f92510fb37c197191d8334401/7acb0a46f21fbe09cfeb63e267600c338644ad8b.jpg  https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D17/sign=f62e05949422720e7fcee6fd79cbebd6/5fdf8db1cb134954f6fa48925a4e9258d0094a64.jpg ,那么在充分接近 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=ac1206297ac6a7efbd26ac22fcfa3ef8/c995d143ad4bd11338f314c456afa40f4bfb050d.jpg 边界的点上, https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D13/sign=fd706518b88f8c54e7d3c12c3b29379b/a50f4bfbfbedab64516eee36fb36afc379311e23.jpg  https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D17/sign=f62e05949422720e7fcee6fd79cbebd6/5fdf8db1cb134954f6fa48925a4e9258d0094a64.jpg 。尤其在边界 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D16/sign=764c2571dca20cf44290fad97409d9cb/d53f8794a4c27d1efdd5993217d5ad6edcc43850.jpg 满足https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D98/sign=181d51c756afa40f38c6c2d5a96459eb/f7246b600c3387448015c9895d0fd9f9d62aa077.jpg其中 https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=3623ead8de1b0ef468e89c5adfc476da/279759ee3d6d55fb201ad4d061224f4a21a4dd6c.jpg 是内部法向矢量场。有符号距离函数因此是法向量场的可微分延伸。特别地, https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=8f98f516fb03918fd3d139ce503d9c3b/a6efce1b9d16fdfa43ee6418b88f8c5495ee7b42.jpg 边界上的带符号距离函数的Hessian给出了Weingarten映射。进一步,如果是一个足够接近Ω边界的区域,那么 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D16/sign=fff895323dadcbef05347a00adafd757/a5c27d1ed21b0ef4348f34a8d1c451da80cb3ee4.jpg 是两次连续可微分的,那么就存在一个明确的公式,该公式涉及Weingarten映射 https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D21/sign=91b21cf3e9dde711e3d244f7a5ef21c6/0eb30f2442a7d9338c454885a14bd11372f00152.jpg 用于根据有符号距离函数改变变量的雅可比矩阵和最近的边界点。具体地,如果 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D12/sign=180c9a8974f0f736dcfe48030b55535d/023b5bb5c9ea15ce85bb4cffba003af33b87b2eb.jpg是距离内的点的集合 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D8/sign=eaec1cf3e9dde711e3d24ec7a778a4/d50735fae6cd7b89a1df9675032442a7d8330e90.jpg 的边界的 https://gss3.bdstatic.com/7Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D14/sign=f5d33dc947fbfbedd859327b79f05777/242dd42a2834349bac47c273c5ea15ce37d3be87.jpg (即管状邻域半径的 https://gss3.bdstatic.com/-Po3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D8/sign=40ad5a3bf0039245a5b5ec3e872451/63d9f2d3572c11df2a8f13d56f2762d0f603c2ee.jpg ),并 https://gss1.bdstatic.com/9vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D12/sign=72616234be119313c343fbb26438c59b/d000baa1cd11728b663220d4c4fcc3cec3fd2c2c.jpg 是上一个绝对积函数Γ,然后

https://gss0.bdstatic.com/-4o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D432/sign=2058bac311950a7b71354fc738d0625c/0d338744ebf81a4c3d5f0d3cdb2a6059242da6d2.jpg

其中 https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D23/sign=e80afb1ffed3572c62e29bdf881383d0/9a504fc2d562853507067ed79cef76c6a6ef6366.jpg 表示行列式 https://gss2.bdstatic.com/-fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D25/sign=090e3d04711ed21b7dc929e0ac6e7710/09fa513d269759eeb3f7cc3dbefb43166c22dfa5.jpg 表示表面积分 [2]  

应用

例如,在计算机视觉中使用签名的距离函数。如使用了一种方法(由Valve Corporation开发),使用GPU加速来渲染大尺寸(或者高DPI)的平滑字体。Valve的方法计算了光栅空间中的有符号距离场,以避免解决(连续)向量空间中问题的计算复杂性。最近已经提出了分段逼近解决方案(例如,近似具有弧形样条的Bézier),但即使如此,计算对于实时渲染而言可能太慢,并且必须通过基于网格的离散化技术来近似(并从计算中剔除)到距离太远的点的距离。

参考资料

  • 1.  [1]董吉文,杨海英,张冰. 水平集方法中符号距离函数的快速生成[J]. 计算机工程与应用,2009,45(34):180-182.
  • 2.  [2] Stanley J. Osher and Ronald P. Fedkiw (2003). Level Set Methods and Dynamic Implicit Surfaces. Springer.

 

發表評論
所有評論
還沒有人評論,想成為第一個評論的人麼? 請在上方評論欄輸入並且點擊發布.
相關文章