這一部分主要講的是怎么去 除區域，，如果我們知道亮點周圍的像素值，，那么我們可以用插值的方法把亮點部分的像素插值成跟亮斑周圍像素值

這里如果想要精確的去除亮斑我們可以用徑向基函數插值法

這里我們用線性插值，簡單操作

假設我們已知坐標(x0,y0)與(x1,y1),要得到[x0,x1]區間內某一位置x在直線上的y值。 根據圖中所示，假設AB上有一點（x,y），可作出兩個相似三角形，我們得到（y-y0)/(x-x0)=(y1-y0)/(x1-x0) 假設方程兩邊的值為α，那么這個值就是插值系數—從x0到x的距離與從x0到x1距離的比值。由于x值已知，所以可以從公式得到α的值 α=(x-x0)/(x1-x0) 同樣，α=(y-y0)/(y1-y0) 這樣，在代數上就可以表示成為： y?= (1- α)y0 + αy1 或者， y?=?y0 + α(y1 -?y0) 這樣通過α就可以直接得到 y。實際上，即使x不在x0到x1之間并且α也不是介于0到1之間，這個公式也是成立的。在這種情況下，這種方法叫作線性外插—參見 外插值。

% Input: % I = input image % angle_delta = angle step size % [crx cry] = corneal reflection center % crr = corneal reflection radius % Output: % I = output image with the corneal reflection removedif crx==0 | cry==0 | crr<=0return; end[height width] = size(I);if crx-crr < 1 || crx+crr > width || cry-crr < 1 || cry+crr > heightfprintf(1, 'Error! Corneal reflection is too near the image border\n');return; endtheta=[0:pi/360:2*pi]; tmat = repmat(theta',[1 crr]);%拓展角度0----360度 rmat = repmat([1:crr],[length(theta) 1]);%一角度的長度為行數，半徑的增幅是1， [xmat,ymat]=pol2cart(tmat,rmat); %在各個角度的上各個半徑的xy坐標值 就是在亮斑區域的坐標值 xv=reshape(xmat,[1 prod(size(xmat))]); yv=reshape(ymat,[1 prod(size(ymat))]);avgmat=ones(size(rmat))*mean(I(round(ymat(:,end)+cry)+(round(xmat(:,end)+crx)-1)*height));%取亮斑邊界的值 permat=repmat(I(round(ymat(:,end)+cry)+(round(xmat(:,end)+crx)-1)*height),[1 crr]); wmat=repmat((1:crr)/crr,[length(theta) 1]); imat=avgmat.*(1-wmat) + permat.*(wmat); %線性插值 在代數上就可以表示成為：y = (1- α)y0 + αy1 I(round(yv+cry)+(round(xv+crx-1).*height))=imat; %亮斑區域的像素值設定

總結

以上是生活随笔為你收集整理的人眼中亮斑的检测、定位和去除（3）的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。