[go: up one dir, main page]

TWI461650B - An image processing system and a method thereof - Google Patents

An image processing system and a method thereof Download PDF

Info

Publication number
TWI461650B
TWI461650B TW101113556A TW101113556A TWI461650B TW I461650 B TWI461650 B TW I461650B TW 101113556 A TW101113556 A TW 101113556A TW 101113556 A TW101113556 A TW 101113556A TW I461650 B TWI461650 B TW I461650B
Authority
TW
Taiwan
Prior art keywords
image processing
light intensity
interference pattern
processing method
tested
Prior art date
Application number
TW101113556A
Other languages
Chinese (zh)
Other versions
TW201344147A (en
Inventor
Chien Ching Ma
Ching Yuan Chang
Original Assignee
Univ Nat Taiwan
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Univ Nat Taiwan filed Critical Univ Nat Taiwan
Priority to TW101113556A priority Critical patent/TWI461650B/en
Publication of TW201344147A publication Critical patent/TW201344147A/en
Application granted granted Critical
Publication of TWI461650B publication Critical patent/TWI461650B/en

Links

Landscapes

  • Image Processing (AREA)
  • Length Measuring Devices By Optical Means (AREA)

Description

影像處理系統及其方法Image processing system and method thereof

本發明係有關一種影像處理系統及其方法,詳而言之,係涉及一種應用於干涉圖樣之影像處理系統及其方法。The present invention relates to an image processing system and method thereof, and more particularly to an image processing system and method thereof for use in an interference pattern.

近年來,光學量測已成功發展出一些高效率的量測技術並廣泛運用到各種變形量測上,如骨頭、IC封裝、或薄膜變形等,其中,電子斑點干涉術(Electronic Speckle Pattern Interferometry;簡稱ESPI)可量測,待測物因壓力、溫度或應力所引起的面內與面外振動、待測物表面的起伏變化情形、以及待測物振動模態等。In recent years, optical metrology has successfully developed some high-efficiency measurement techniques and is widely used in various deformation measurements, such as bone, IC packaging, or film deformation, among which, Electronic Speckle Pattern Interferometry (Electronic Speckle Pattern Interferometry; Referring to ESPI, it can measure the in-plane and out-of-plane vibration caused by pressure, temperature or stress, the fluctuation of the surface of the object to be tested, and the vibration mode of the object to be tested.

當雷射光照射表面粗糙的待測物時,由於其凹凸不平的表面導致反射光因散射而在空間互相干涉進而在反射處形成複雜的光強分布,此即為斑點效應。此外,當待測物表面產生位移時,斑點之分佈情形亦隨之變化。因此,電子斑點干涉術係根據斑點的干涉原理,於待測物變形前後利用CCD擷取影像,再透過影像軟體疊像產生干涉現象,此干涉現象不僅記錄待測物反射光的振幅,同時也記錄了反射光的相位資料等。When the laser light illuminates the object to be tested with a rough surface, the uneven light-emitting surface causes the reflected light to interfere with each other due to scattering and form a complex light intensity distribution at the reflection, which is a speckle effect. In addition, when the surface of the object to be tested is displaced, the distribution of the spots also changes. Therefore, according to the interference principle of the speckle, the electron speckle interferometry uses the CCD to capture the image before and after the deformation of the object to be detected, and then interferes with the image superficial image to generate an interference phenomenon, which not only records the amplitude of the reflected light of the object to be tested, but also The phase data of the reflected light and the like are recorded.

再者,電子斑點干涉術為全域式量測技術,但量測結果往往由於雜訊,例如電子儀器、光學設備或機械擾動所產生之雜訊等,而導致影像品質降低,其中電子儀器所產生的為背景雜訊,此外,藉由雷射光和透鏡能減少光學設備所產生的雜訊,但關於機械擾動(一般是因為空氣擾動)則難以移除。Furthermore, electronic speckle interferometry is a global measurement technique, but the measurement results are often caused by noise, such as noise generated by electronic instruments, optical equipment or mechanical disturbances, resulting in reduced image quality, which is generated by electronic instruments. For background noise, in addition, lasers and lenses can reduce the noise generated by optical devices, but mechanical disturbances (generally due to air disturbances) are difficult to remove.

於習知技術中多利用相減法或平均法濾除雜訊以提升影像品質,然而,相減法可減少背景雜訊的影響但難以消除擾動雜訊,而平均法可消除擾動雜訊但仍殘留有背景雜訊。目前,尚未能提出一種能消除或降低背景雜訊和擾動雜訊之量測技術。In the prior art, the subtraction method or the averaging method is used to filter out noise to improve image quality. However, the subtraction method can reduce the influence of background noise but it is difficult to eliminate the disturbance noise, and the averaging method can eliminate the disturbance noise but still remains. There is background noise. At present, a measurement technique that can eliminate or reduce background noise and disturbance noise has not been proposed.

為解決前述習知技術之種種問題,本發明提供一種影像處理系統及其方法,用以降低影像的背景雜訊和擾動雜訊。To solve the above problems of the prior art, the present invention provides an image processing system and method thereof for reducing background noise and disturbance noise of an image.

本發明之影像處理方法,係應用於干涉圖樣,該影像處理方法包括以下步驟:(1)利用感光元件對該干涉圖樣擷取不同曝光時間的影像序列,以使該感光元件之複數個像素點取得該干涉圖樣之光強度;(2)基於該曝光時間,計算各該像素點上的光強度平均值;以及(3)基於該影像序列的張數,對各該像素點上的光強度平均值執行標準差之運算,俾利用該標準差之運算結果消除該干涉圖樣的背景雜訊及擾動雜訊。此外,還可包括步驟(4)對各該像素點上之該標準差之運算結果執行頻率掃描,以取得該待測物之共振頻及對應該共振頻之干涉圖樣的影像。The image processing method of the present invention is applied to an interference pattern, and the image processing method comprises the following steps: (1) using a photosensitive element to capture an image sequence of different exposure time for the interference pattern, so as to make a plurality of pixels of the photosensitive element Obtaining the light intensity of the interference pattern; (2) calculating an average value of the light intensity at each of the pixel points based on the exposure time; and (3) averaging the light intensity at each of the pixel points based on the number of sheets of the image sequence The value performs the operation of the standard deviation, and the background noise and the disturbance noise of the interference pattern are eliminated by using the result of the standard deviation. In addition, the step (4) may be performed to perform frequency scanning on the operation result of the standard deviation on each of the pixel points to obtain an image of the resonance frequency of the object to be tested and the interference pattern corresponding to the resonance frequency.

本發明之影像處理系統,係應用於待測物之干涉圖樣,該影像處理系統係包括:感光元件,用以對該干涉圖樣擷取不同曝光時間的影像序列,以使所具有之複數個像素點取得該干涉圖樣之光強度;以及處理器,係連接該感光元件並取得該干涉圖樣之光強度,該處理器包括:時域平均模組,用以基於該曝光時間,計算各該像素點上的光強度平均值;及雜訊消除模組,用以基於該影像序列的影像張數,對各該像素點上的光強度平均值執行標準差之運算,俾利用該標準差之運算結果消除該干涉圖樣的背景雜訊及擾動雜訊。The image processing system of the present invention is applied to an interference pattern of a sample to be tested, the image processing system comprising: a photosensitive element for capturing an image sequence of different exposure time for the interference pattern, so as to have a plurality of pixels Obtaining a light intensity of the interference pattern; and a processor connecting the photosensitive element and obtaining a light intensity of the interference pattern, the processor comprising: a time domain averaging module, configured to calculate each pixel point based on the exposure time The average value of the light intensity; and the noise cancellation module is configured to perform a standard deviation operation on the average value of the light intensity at each pixel point based on the number of images of the image sequence, and use the standard deviation operation result Eliminate background noise and disturbance noise of the interference pattern.

於本發明之影像處理系統及其方法中,所述之干涉圖樣可由兩同調光同時照射一面內振動的待測物,該兩同調光所形成之斑點相互干涉而產生;或可由兩同調光分別照射一參考面和一面外振動的待測物,該兩同調光所形成之斑點相互干涉而產生。In the image processing system and method of the present invention, the interference pattern may be irradiated by two dimming lights simultaneously to illuminate an object to be tested vibrating in one side, and the spots formed by the two dimming lights interfere with each other; or may be separately dimmed by two The object to be tested is irradiated with a reference surface and an external vibration, and the spots formed by the two dimming light interfere with each other to generate.

此外,於該雜訊消除模組尚未執行標準差運算之前,各該像素點上的光強度平均值仍包括該擾動雜訊,且該擾動雜訊的發生機率可利用高斯分佈函數來估計。In addition, before the noise cancellation module has performed the standard deviation operation, the average light intensity at each pixel point still includes the disturbance noise, and the probability of occurrence of the disturbance noise can be estimated by using a Gaussian distribution function.

相較於習知技術,本發明將時間域標準差(temporal standard deviation,TSTD)演算法應用在電子斑點干涉術上,可同時兼顧相減法與平均法之優點,更可降低空氣擾動雜訊和背景雜訊對電子斑點干涉術的影響,以獲得高解析度的干涉圖樣,提升影像品質。Compared with the prior art, the present invention applies the temporal standard deviation (TSTD) algorithm to the electronic speckle interferometry, and can simultaneously take advantage of the subtraction method and the averaging method, and can reduce the air disturbance noise and The effect of background noise on electronic speckle interferometry to obtain high-resolution interference patterns and improve image quality.

以下藉由特定的具體實施形態說明本發明之實施方式,熟悉此技術之人士可由本說明書所揭示之內容輕易地了解本發明之其他優點與功效,亦可藉由其他不同的具體實施形態加以施行或應用。The embodiments of the present invention are described in the following specific embodiments, and those skilled in the art can easily understand other advantages and functions of the present invention by the disclosure of the present disclosure, and can also be implemented by other different embodiments. Or application.

第1A圖係概略繪示量測面內振動的待測物之光學架構圖,第1B圖係用以表示待測物的面內振動之示意圖,得藉由第1A圖所示之光學架構對第1B圖所示之待測物的干涉圖樣擷取不同曝光時間的影像序列。FIG. 1A is a schematic diagram showing an optical architecture of a test object in a measurement surface vibration, and FIG. 1B is a schematic diagram showing an in-plane vibration of the object to be tested, which is obtained by the optical architecture shown in FIG. 1A. The interference pattern of the object to be tested shown in Fig. 1B captures an image sequence of different exposure times.

待測物10,例如壓電片,可藉由函數波產生器11和功率放大器12而呈面內振動狀態,如第1B圖所示,面內振動包括U方向和V方向的位移。The object to be tested 10, for example, a piezoelectric piece, can be in an in-plane vibration state by the function wave generator 11 and the power amplifier 12, and as shown in FIG. 1B, the in-plane vibration includes displacements in the U direction and the V direction.

光束發射器13發射同調光,例如氦-氖雷射,該同調光經分光器14分成兩束光,此兩束光分別由空間濾波器16a和16b濾除空間雜訊後,以等角度入射在待測物10的法線方向。而待測物10將該兩束光反射,以使該兩束光所形成的斑點相互干涉而產生干涉圖樣。The beam emitter 13 emits the same dimming light, such as a 氦-氖 laser, which is split into two beams by the beam splitter 14 which are spatially filtered by the spatial filters 16a and 16b, respectively, and incident at equal angles. In the normal direction of the object 10 to be tested. The object to be tested 10 reflects the two beams of light so that the spots formed by the two beams interfere with each other to generate an interference pattern.

感光元件17放置於待測物10的法線方向上,以對該干涉圖樣擷取不同曝光時間的影像序列,使感光元件17的複數個像素點取得該干涉圖樣之光強度,並將該些光強度傳輸至處理器18及將所擷取之影像顯示在顯示器19上。於本光學架構中,圖式中處理器18為電腦,然亦可以其他任何具有運算功能之電子裝置代替,處理器18對該影像之影像處理係於後文敘述。The photosensitive element 17 is placed in the normal direction of the object to be tested 10, and the image sequence of the different exposure time is captured for the interference pattern, so that the plurality of pixels of the photosensitive element 17 obtain the light intensity of the interference pattern, and the light intensity is obtained. The light intensity is transmitted to the processor 18 and the captured image is displayed on the display 19. In the present optical architecture, the processor 18 in the figure is a computer, but it can be replaced by any other electronic device having a computing function, and the image processing of the image by the processor 18 will be described later.

第2A圖係概略繪示量測面外振動的待測物之光學架構圖,第2B圖係用以表示待測物的面外振動之示意圖,得藉由第2A圖所示之光學架構對第2B圖所示之待測物的干涉圖樣擷取不同曝光時間的影像序列。2A is a schematic diagram showing the optical architecture of the object to be tested for vibration outside the surface, and FIG. 2B is a schematic diagram showing the vibration of the surface of the object to be tested, which is obtained by the optical architecture shown in FIG. 2A. The interference pattern of the object to be tested shown in Fig. 2B captures an image sequence of different exposure times.

待測物的面外振動係如第2B圖所示,係在W方向的 位移,第2A圖與第1A圖之差異在於,同調光經反射器15反射及空間濾波器16消除空間雜訊後,分光器14’將該同調光分成兩束,一束光照射於待測物10上再反射回感光元件17(物光),另一束光照射在參考板20上當作參考光,物光及參考光所形成的斑點相互干涉而產生干涉圖樣。感光元件17放置於待測物10的法線方向上,以對該干涉圖樣擷取不同曝光時間的影像序列,使感光元件17的複數個像素點取得該干涉圖樣之光強度,並將該些光強度傳輸至處理器18及將所擷取之影像顯示在顯示器19上。於本光學架構中,圖式中處理器18為電腦,然亦可以其他任何具有運算功能之電子裝置代替,處理器18對該影像之影像處理係於後文敘述。The out-of-plane vibration of the object to be tested is shown in Figure 2B and is in the W direction. The difference between the 2A and 1A is that after the dimming is reflected by the reflector 15 and the spatial filter 16 cancels the spatial noise, the beam splitter 14' splits the dimming into two beams, and one beam is irradiated to be tested. The object 10 is reflected back to the photosensitive element 17 (object light), and the other light is irradiated on the reference plate 20 as a reference light, and the spots formed by the object light and the reference light interfere with each other to generate an interference pattern. The photosensitive element 17 is placed in the normal direction of the object to be tested 10, and the image sequence of the different exposure time is captured for the interference pattern, so that the plurality of pixels of the photosensitive element 17 obtain the light intensity of the interference pattern, and the light intensity is obtained. The light intensity is transmitted to the processor 18 and the captured image is displayed on the display 19. In the present optical architecture, the processor 18 in the figure is a computer, but it can be replaced by any other electronic device having a computing function, and the image processing of the image by the processor 18 will be described later.

需說明的是,本發明之影像處理系統及方法可用以針對面內振對、面外振動及整體振動等振動模式之待測物之干涉圖像進行影像處理,其中,面內振動為待測物平行於其水平面振動、面外振動為待測物垂直於其水平面振動。上述第1A和1B僅為舉例說明之光學架構,所屬技術中具有通常知識者應可針對所需量測的振動模式而改變光學架構。It should be noted that the image processing system and method of the present invention can perform image processing on an interference image of a test object of a vibration mode such as an in-plane vibration, an out-of-plane vibration, and an overall vibration, wherein the in-plane vibration is to be tested. The object vibrates parallel to its horizontal plane, and the out-of-plane vibration causes the object to be tested to vibrate perpendicular to its horizontal plane. The above 1A and 1B are merely illustrative optical architectures, and those of ordinary skill in the art should be able to change the optical architecture for the desired measured vibration modes.

第3圖係為本發明之影像處理系統之方塊圖。本發明之影像處理系統基本包括感光元件31和處理器32。Figure 3 is a block diagram of the image processing system of the present invention. The image processing system of the present invention basically includes a photosensitive element 31 and a processor 32.

感光元件31可例如感光耦合元件(Charge-coupled Device,CCD)或互補式金氧半導體(Complementary Metal-Oxide-Semiconductor,CMOS),用以對該干涉圖樣擷取不同曝光時間的影像序列,以使感光元件31的複數個像素點取得該干涉圖樣之光強度。換言之,感光元件31每隔一段時間(固定或非固定時段皆可)對待測物拍照,其中,該待測物可呈面內振動、面外振動或整體振動狀態, 則感光元件31的複數個像素點便取得該干涉圖樣的光強度,再將每一個像素點所取得之光強度的電荷信號傳輸至處理器32。The photosensitive element 31 can be, for example, a Charge-coupled Device (CCD) or a Complementary Metal-Oxide-Semiconductor (CMOS), for capturing an image sequence of different exposure time for the interference pattern, so that The plurality of pixels of the photosensitive element 31 acquire the light intensity of the interference pattern. In other words, the photosensitive element 31 takes a picture of the object to be measured at intervals (either fixed or non-fixed time period), wherein the object to be tested may be in-plane vibration, out-of-plane vibration or overall vibration state. Then, the plurality of pixels of the photosensitive element 31 acquire the light intensity of the interference pattern, and then transmit the charge signal of the light intensity obtained by each pixel to the processor 32.

處理器32,係連接感光元件31並接收各該像素點所取得之干涉圖樣之光強度的電荷訊號。處理器32包括時域平均模組321和雜訊消除模組322。時域平均模組321,用以基於該曝光時間計算各該像素點上的光強度平均值。雜訊消除模組322,用以基於影像序列的影像張數對在各該像素點上的光強度平均值執行標準差之運算。因此,本發明之影像處理系統能利用標準差消除之干涉圖樣的背景雜訊及擾動雜訊,以增進影像品質。The processor 32 is connected to the photosensitive element 31 and receives a charge signal of the intensity of the interference pattern obtained by each pixel. The processor 32 includes a time domain averaging module 321 and a noise cancellation module 322. The time domain averaging module 321 is configured to calculate an average value of the light intensity at each of the pixel points based on the exposure time. The noise cancellation module 322 is configured to perform a standard deviation operation on the average value of the light intensity at each pixel point based on the number of images of the image sequence. Therefore, the image processing system of the present invention can utilize the standard noise to eliminate the background noise and disturbance noise of the interference pattern to enhance the image quality.

以下係以第4圖所示之流程圖配合演算法來說明本發明之影像處理方法。Hereinafter, the image processing method of the present invention will be described with the flowchart shown in FIG. 4 in conjunction with an algorithm.

於步驟S301中,利用感光元件對待測物所產生之干涉圖樣擷取不同曝光時間的影像序列,則使感光元件之複數個像素點取得該干涉圖樣的光強度。In step S301, the image sequence of the different exposure time is captured by the interference pattern generated by the photosensitive element to be measured, and the light intensity of the interference pattern is obtained by a plurality of pixels of the photosensitive element.

於步驟S302中,基於選定曝光時間(i τ),計算各該像素點上的光強度平均值(I i (x ,y ))。詳言之,兩同調光發生干涉後,干涉圖樣的光強度公式為: 其中,若考慮擾動所造成的雜訊,則於公式中增加ψ (x ,y ,t ),則干涉圖樣的光強度公式為: 則,各像素點的光強度時間平均值為: 其中,i 表示影像序列中的第i張影像,(x,y)表示一像素點,J 0 為零階貝索函數,Γ為靈敏度因子,(x ,y )表示相位,A (x ,y )表示振幅。又,擾動雜訊ψ (x ,y ,t )發生的機率可以高斯分布函數來估計,因此,各像素點的光強度時間平均值為: 其中,I0 (x ,y )+I R (x ,y )表示背景雜訊,ψ i sinφ (x ,y )表示擾動雜訊。In step S302, an average light intensity ( I i ( x , y )) at each of the pixel points is calculated based on the selected exposure time (i τ). In detail, after the interference of the two dimming lights, the light intensity formula of the interference pattern is: If we consider the noise caused by the disturbance, add ψ ( x , y , t ) to the formula, then the light intensity formula of the interference pattern is: Then, the light intensity time average of each pixel is: Where i denotes the ith image in the image sequence, (x, y) denotes a pixel point, J 0 is a zero-order Bessel function, and Γ is a sensitivity factor. ( x , y ) represents the phase, and A ( x , y ) represents the amplitude. Moreover, the probability of disturbing the noise ψ ( x , y , t ) can be estimated by a Gaussian distribution function. Therefore, the average value of the light intensity time of each pixel is: Where I 0 ( x , y )+ I R ( x , y ) represents background noise, and ψ i sin φ ( x , y ) represents disturbance noise.

於步驟S303中,基於影像序列的影像張數(n ),對各該像素點上的光強度平均值執行標準差運算。光強度的標準差演算式為: 其中,為不包含時間的常數,亦即,該擾動雜訊於大量的影像序列時會收斂至常數,故可消除背景雜訊和擾動雜訊。In step S303, a standard deviation operation is performed on the average value of the light intensity at each of the pixel points based on the number of images ( n ) of the image sequence. The standard deviation formula of light intensity is: among them, It is a constant that does not contain time, that is, the disturbance noise converges to a constant when a large number of image sequences are used, thereby eliminating background noise and disturbing noise.

另外,在消除背景雜訊和擾動雜訊之後,計算在一特定振動頻率下該待測物所接收到的光強度,以對該待測物執行頻率掃描(frequency-sweeping)。可以下列演算法來表達: 其中LW 分別為該待測物的長和寬,ω 為特定頻率。In addition, after the background noise and the disturbance noise are eliminated, the light intensity received by the object to be tested is calculated at a specific vibration frequency to perform frequency-sweeping on the object to be tested. It can be expressed by the following algorithms: Where L and W are the length and width of the object to be tested, respectively, and ω is a specific frequency.

具體實施時,首先分別於每一特定頻率擷取例如50張影像序列,接著執行時域平均和標準差運算,並執行E (ω )運算以建立頻率掃描曲線,如第5圖所示,接著選擇相對應的影像。參閱第5圖,其顯示在各個特定頻率下之待測物所接收到的光強度之標準化,並繪示最適曲線(fitting curve),由該最適曲線的轉折得知該待測物之共振頻率,例如約852Hz,且該共振頻之下,感光元件所取得之該待測物的干涉圖樣的影像品質係為較佳。In a specific implementation, first, for example, 50 image sequences are respectively taken at each specific frequency, then time domain average and standard deviation operations are performed, and an E ( ω ) operation is performed to establish a frequency scan curve, as shown in FIG. 5, and then Select the corresponding image. Referring to FIG. 5, it shows the normalization of the light intensity received by the object under test at each specific frequency, and shows the fitting curve, and the resonance frequency of the object to be tested is obtained from the turning of the optimum curve. For example, about 852 Hz, and below the resonance frequency, the image quality of the interference pattern of the object to be tested obtained by the photosensitive element is preferable.

參閱第6A、6B及6C圖,其分別為利用平均法、相減法、標準差法之影像處理方法所擷取之影像。第6C圖相較於第6A和6B圖明顯地提升影像品質,平均法仍殘留有背景雜訊。相減法的演算式為: 其中相減法之結果仍殘留有擾動雜訊|ψ j |-|ψ i |。Refer to Figures 6A, 6B, and 6C, which are images captured by the image processing method of the averaging method, the subtraction method, and the standard deviation method, respectively. The 6C picture significantly improves the image quality compared to the 6A and 6B pictures, and the background method still has background noise. The calculation formula of the subtraction method is: The result of the subtraction method still has disturbing noise | ψ j |-| ψ i |.

綜上所述,本發明之影像處理系統及其方法,主要對干涉圖樣之不同曝光時間的影像序列執行時間平均運算和標準差運算,亦即時間域標準差演算法,以消除背景雜訊和空氣擾動,進而提升影像品質。此外,本發明除了量測壓電片的微小形變外,更可應用於壓電片之電極、壓電換能器、高頻振動元件之設計與開發。In summary, the image processing system and method thereof of the present invention mainly perform time averaging operation and standard deviation operation, that is, time domain standard deviation algorithm, on image sequences of different exposure times of interference patterns, to eliminate background noise and Air disturbances, which in turn improve image quality. In addition, the present invention can be applied to the design and development of electrodes, piezoelectric transducers, and high-frequency vibration elements of piezoelectric sheets in addition to measuring minute deformation of piezoelectric sheets.

上述實施形態僅例示性說明本發明之原理、特點及其功效,並非用以限制本發明之可實施範疇,任何熟習此項技藝之人士均可在不違背本發明之精神及範疇下,對上述實施形態進行修飾與改變。任何運用本發明所揭示內容而完成之等效改變及修飾,均仍應為下述之申請專利範圍所涵蓋。因此,本發明之權利保護範圍,應如後述之申請專利範圍所列。The above-described embodiments are merely illustrative of the principles, features, and effects of the present invention, and are not intended to limit the scope of the present invention. Any person skilled in the art can recite the above without departing from the spirit and scope of the present invention. The embodiment is modified and changed. Any equivalent changes and modifications made by the disclosure of the present invention should still be covered by the following claims. Therefore, the scope of protection of the present invention should be as set forth in the scope of the claims described below.

10...待測物10. . . Analyte

11...函數波產生器11. . . Function wave generator

12...功率放大器12. . . Power amplifier

13...光束發射器13. . . Beam emitter

14、14’...分光器14, 14’. . . Splitter

15...反射器15. . . reflector

16、16a、16b...空間濾波器16, 16a, 16b. . . Spatial filter

17...感光元件17. . . Photosensitive element

18...處理器18. . . processor

19...顯示器19. . . monitor

20...參考板20. . . Reference board

31...感光元件31. . . Photosensitive element

32...處理器32. . . processor

321...時域平均模組321. . . Time domain average module

322...雜訊消除模組322. . . Noise cancellation module

U、V、W...方向U, V, W. . . direction

S301~S303...步驟S301~S303. . . step

第1A圖係概略繪示量測面內振動的待測物之光學架構圖;第1B圖係用以表示待測物的面內振動之示意圖:第2A圖係概略繪示量測面外振動的待測物之光學架構圖;第2B圖係用以表示待測物的面外振動之示意圖;第3圖係為本發明之影像處理系統之方塊圖;第4圖係為本發明之影像處理方法之流程圖;第5圖係為本發明之影像處理方法之執行頻率掃描之說明圖;以及Figure 1A is a schematic diagram showing the optical architecture of the object under test in the measurement plane; Figure 1B is a schematic diagram showing the in-plane vibration of the object to be tested: Figure 2A is a schematic diagram showing the measurement of the out-of-plane vibration FIG. 2B is a schematic diagram showing the out-of-plane vibration of the object to be tested; FIG. 3 is a block diagram of the image processing system of the present invention; FIG. 4 is an image of the present invention; A flowchart of a processing method; FIG. 5 is an explanatory diagram of an execution frequency scan of the image processing method of the present invention;

第6A至6C圖分別為利用平均演算法、相減演算法、標準差演算法之影像處理方法所擷取之影像。Figures 6A to 6C are images captured by an image processing method using an average algorithm, a subtraction algorithm, and a standard difference algorithm, respectively.

S301~S303...步驟S301~S303. . . step

Claims (10)

一種影像處理方法,係應用於干涉圖樣,該影像處理方法包括以下步驟:(1)利用感光元件對該干涉圖樣擷取不同曝光時間的影像序列,以使該感光元件之複數個像素點取得該干涉圖樣之光強度;(2)基於該曝光時間,計算各該像素點上的光強度平均值;以及(3)基於該影像序列的影像張數,對各該像素點上的光強度平均值執行標準差之運算,俾利用該標準差之運算結果消除該干涉圖樣的背景雜訊及擾動雜訊。An image processing method is applied to an interference pattern, the image processing method comprising the following steps: (1) using a photosensitive element to capture an image sequence of different exposure time for the interference pattern, so that the plurality of pixels of the photosensitive element acquire the image Interfering with the light intensity of the pattern; (2) calculating an average value of the light intensity at each pixel point based on the exposure time; and (3) averaging the light intensity at each of the pixel points based on the number of images of the image sequence The standard deviation operation is performed, and the background noise and the disturbance noise of the interference pattern are eliminated by using the operation result of the standard deviation. 如申請專利範圍第1項所述之影像處理方法,其中,於步驟(1)之前復包括:令兩同調光照射一待測物,以使該兩同調光所形成的斑點相互干涉而產生該干涉圖樣之步驟。The image processing method according to claim 1, wherein before the step (1), the method further comprises: dimming the same object to be irradiated, so that the spots formed by the two dimming light interfere with each other to generate the The step of interfering with the pattern. 如申請專利範圍第2項所述之影像處理方法,其中,該待測物係呈面內振動狀態。The image processing method according to claim 2, wherein the object to be tested is in an in-plane vibration state. 如申請專利範圍第1項所述之影像處理方法,其中,於步驟(1)之前復包括:令兩同調光分別照射一參考物和一待測物,以使該兩同調光所形成的斑點相互干涉而產生該干涉圖樣之步驟。The image processing method of claim 1, wherein before the step (1), the method further comprises: dimming the two dimming lights to respectively irradiate a reference object and a test object to make the spot formed by the dimming The step of interfering with each other to produce the interference pattern. 如申請專利範圍第4項所述之影像處理方法,其中,該待測物係呈面外振動狀態。The image processing method according to claim 4, wherein the object to be tested is in an out-of-plane vibration state. 如申請專利範圍第1項所述之影像處理方法,其中,該待測物係呈整體振動狀態。The image processing method according to claim 1, wherein the object to be tested is in an overall vibration state. 如申請專利範圍第1項所述之影像處理方法,其中,於步驟(2)中,各該像素點上的光強度平均值係包括該擾動雜訊。The image processing method according to claim 1, wherein in step (2), the average value of the light intensity at each of the pixels includes the disturbance noise. 如申請專利範圍第1項所述之影像處理方法,其中,於步驟(2)中,該擾動雜訊的發生機率係利用機率分佈函數來估計。The image processing method according to claim 1, wherein in step (2), the probability of occurrence of the disturbance noise is estimated by using a probability distribution function. 如申請專利範圍第1項所述之影像處理方法,復包括步驟(4),對各該像素點上之該標準差之運算結果執行頻率掃描,以取得該待測物之共振頻及對應該共振頻之干涉圖樣的影像。The image processing method of claim 1, further comprising the step (4) of performing a frequency sweep on the operation result of the standard deviation at each pixel to obtain a resonance frequency of the object to be tested and correspondingly An image of the interference pattern of the resonant frequency. 一種影像處理系統,係應用於待測物之干涉圖樣,該影像處理系統係包括:感光元件,用以對該干涉圖樣擷取不同曝光時間的影像序列,以使所具有之複數個像素點取得該干涉圖樣之光強度;以及處理器,係連接該感光元件以取得該干涉圖樣之光強度,該處理器包括:時域平均模組,用以基於該曝光時間,計算各該像素點上的光強度平均值;及雜訊消除模組,用以基於該影像序列的影像張數,對各該像素點上的光強度平均值執行標準差之運算,俾利用該標準差之運算結果消除該干涉圖樣的背景雜訊及擾動雜訊。An image processing system is applied to an interference pattern of a sample to be tested, the image processing system comprising: a photosensitive element for capturing an image sequence of different exposure time for the interference pattern, so as to obtain a plurality of pixel points a light intensity of the interference pattern; and a processor connecting the photosensitive element to obtain a light intensity of the interference pattern, the processor comprising: a time domain averaging module, configured to calculate each pixel point based on the exposure time An average value of the light intensity; and a noise cancellation module for performing a standard deviation operation on the average value of the light intensity at each pixel point based on the number of images of the image sequence, and eliminating the result by using the standard deviation Interference noise background noise and disturbance noise.
TW101113556A 2012-04-17 2012-04-17 An image processing system and a method thereof TWI461650B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
TW101113556A TWI461650B (en) 2012-04-17 2012-04-17 An image processing system and a method thereof

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
TW101113556A TWI461650B (en) 2012-04-17 2012-04-17 An image processing system and a method thereof

Publications (2)

Publication Number Publication Date
TW201344147A TW201344147A (en) 2013-11-01
TWI461650B true TWI461650B (en) 2014-11-21

Family

ID=49990151

Family Applications (1)

Application Number Title Priority Date Filing Date
TW101113556A TWI461650B (en) 2012-04-17 2012-04-17 An image processing system and a method thereof

Country Status (1)

Country Link
TW (1) TWI461650B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MX2017008537A (en) 2014-12-27 2018-02-21 Guardian Optical Tech Ltd System and method for detecting surface vibrations.
TWI727497B (en) * 2019-11-11 2021-05-11 瑞昱半導體股份有限公司 Image processing method based on sensor characteristics

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6128082A (en) * 1998-09-18 2000-10-03 Board Of Trustees Operating Michigan State University Technique and apparatus for performing electronic speckle pattern interferometry
US6943870B2 (en) * 2002-09-20 2005-09-13 President Of Saitama University Deformation measuring method and apparatus using electronic speckle pattern interferometry
TWI292471B (en) * 2006-10-19 2008-01-11 Metal Ind Res & Dev Ct
TW200844395A (en) * 2006-11-30 2008-11-16 Corning Inc Method and apparatus for distortion measurement imaging

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6128082A (en) * 1998-09-18 2000-10-03 Board Of Trustees Operating Michigan State University Technique and apparatus for performing electronic speckle pattern interferometry
US6943870B2 (en) * 2002-09-20 2005-09-13 President Of Saitama University Deformation measuring method and apparatus using electronic speckle pattern interferometry
TWI292471B (en) * 2006-10-19 2008-01-11 Metal Ind Res & Dev Ct
TW200844395A (en) * 2006-11-30 2008-11-16 Corning Inc Method and apparatus for distortion measurement imaging

Also Published As

Publication number Publication date
TW201344147A (en) 2013-11-01

Similar Documents

Publication Publication Date Title
US10254252B2 (en) Surface and subsurface detection sensor
JP6574989B2 (en) Stress measuring apparatus, stress measuring system and stress measuring method
US11692810B2 (en) Photoacoustic excitation sensing enhanced by cross-correlated unfocused speckle images
US20090180124A1 (en) Method and apparatus for identifying dynamic characteristics of a vibratory object
Zhong et al. Quasi-optical coherence vibration tomography technique for damage detection in beam-like structures based on auxiliary mass induced frequency shift
US11867611B2 (en) Optical-interference analysis
Zaletelj et al. Laser-light speckle formation for deflection-shape identification using digital image correlation
US11187648B2 (en) Two-stage photo-acoustic excitation sensing systems and methods
JP7126670B2 (en) Defect detection method and apparatus using speckle image
TWI461650B (en) An image processing system and a method thereof
JP2000111490A (en) Painted surface detection device
Hermawanto et al. Determination of microphone acoustic center from sound field projection measured by optical interferometry
CN113167767B (en) Defect inspection device and defect inspection method
JP2017101976A (en) Inspection system and inspection method
JP6047764B2 (en) White interferometer, image processing method, and image processing program
TW201814239A (en) Interferometer and imaging method therefor
JP5518187B2 (en) Deformation measurement method
JP7480915B2 (en) Defect inspection device and defect inspection method
CN108254379B (en) Defect detection device and method
Zhong et al. Quasi‐OCVT technique for response‐only experimental modal analysis of beam‐like structures
JP6457846B2 (en) Method and apparatus for measuring shape of transparent plate
Yano et al. Real-time three-dimensional vibration mode visualization system based on laser speckle interference
Bointon et al. The effects of vibration on fringe projection systems
JP4210236B2 (en) Optical shape measurement method
Wang et al. High precision full-field vibration measurement by LDV-induced stroboscopic fringe projection