*************第一 窗口初始化****************
*關(guān)閉已經(jīng)打開的窗口
dev_close_window () 

*關(guān)閉程序計(jì)數(shù)器,圖形變量更新,窗口圖形更新
dev_update_off () 

*讀取一張圖像
read_image (Image,'2.bmp') 

*獲取圖像大小
get_image_size (Image, Width, Height) 

*打開新的窗口
dev_open_window (0, 0, 640, 480, 'black', WindowHandle) 

*設(shè)置區(qū)域填充方式
dev_set_draw ('margin') 

*設(shè)置線寬度
dev_set_line_width (2) 

*設(shè)置輸出對(duì)象顯示顏色數(shù)目
dev_set_colored(12) 


*******************第二 圖像預(yù)處理**********************
*優(yōu)化FFT的速度
optimize_rft_speed (Width, Height, 'standard') 
Sigma1 := 10.0 
Sigma2 := 2.0 
*構(gòu)建高斯濾波器
gen_gauss_filter (GsFilter1, Sigma1, Sigma1, 0.0, 'none', 'rft', Width, Height) 
*構(gòu)建高斯濾波器
gen_gauss_filter (GsFilter2, Sigma2, Sigma2, 0.0, 'none', 'rft', Width, Height)

*兩個(gè)濾波器相減
sub_image (GsFilter1, GsFilter2, Filter, 1, 0) 


********************第三 開始處理圖像*************
* Image Acquisition 01: Code generated by Image Acquisition 01
ImageFiles := []
ImageFiles[0] := '1.bmp'
ImageFiles[1] := '2.bmp'



for Index := 0 to |ImageFiles| - 1 by 1
    
*讀取一張圖像
read_image (Image, ImageFiles[Index])

*將RGB彩色圖轉(zhuǎn)化為灰度圖
rgb1_to_gray (Image, Image) 

*對(duì)灰度圖進(jìn)行反選
invert_image(Image, ImageInvert) 
   
*將圖像轉(zhuǎn)化為頻域圖像
rft_generic (ImageInvert, ImageFFT, 'to_freq', 'sqrt', 'complex', Width) 

*將頻域圖像和濾波核進(jìn)行卷積運(yùn)算
convol_fft (ImageFFT, Filter, ImageConvol) 

*將頻域圖像轉(zhuǎn)化為空間域圖像
rft_generic (ImageConvol, ImageFiltered, 'from_freq', 'n', 'real', Width) 

*創(chuàng)建一個(gè)ROI
gen_rectangle1(Rectangle,30,25,450,610) 

*將ROI區(qū)域的圖像剪切
reduce_domain(ImageFiltered, Rectangle, ROI) 

*中值濾波
median_image(ROI, ImageMedian, 'circle', 20, 'mirrored') 

*使用高斯核進(jìn)行圖像平滑
smooth_image (ROI, ImageSmooth, 'gauss', 8) 

*通過(guò)閾值分割出缺陷區(qū)域
threshold (ImageSmooth, Regions, -0.0016783, -0.0006434) 

*對(duì)缺陷區(qū)域進(jìn)行連通處理
connection(Regions, ConnectedRegions) 

*對(duì)區(qū)域進(jìn)行膨脹操作
dilation_circle (ConnectedRegions, RegionDilation, 5) 

*過(guò)濾出指定面積大小的缺陷
select_shape (RegionDilation, SelectedRegions, 'area', 'and', 800, 99999) 

*根據(jù)區(qū)域生成XLD輪廓
gen_contour_region_xld (SelectedRegions, Contours, 'border') 

*將XLD擬合成圓
fit_circle_contour_xld (Contours, 'atukey', -1, 2, 0, 5, 2, Row, Column, Radius, StartPhi, EndPhi, PointOrder) 

*生成一個(gè)圓對(duì)象
gen_circle_contour_xld (ContCircle, Row, Column, Radius+20, 0, 6.28318, 'positive', 1)


*顯示圖像
dev_display (Image) 


*顯示圓標(biāo)記
dev_display(ContCircle) 

stop() 

endfor