JP5684505B2 - Image processing apparatus, image processing method, and computer program - Google Patents
Image processing apparatus, image processing method, and computer program Download PDFInfo
- Publication number
- JP5684505B2 JP5684505B2 JP2010161125A JP2010161125A JP5684505B2 JP 5684505 B2 JP5684505 B2 JP 5684505B2 JP 2010161125 A JP2010161125 A JP 2010161125A JP 2010161125 A JP2010161125 A JP 2010161125A JP 5684505 B2 JP5684505 B2 JP 5684505B2
- Authority
- JP
- Japan
- Prior art keywords
- image
- pixel
- value
- distribution range
- distance
- Prior art date
- Legal status (The legal status 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 status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Investigating Materials By The Use Of Optical Means Adapted For Particular Applications (AREA)
- Image Processing (AREA)
- Image Analysis (AREA)
- Length Measuring Devices By Optical Means (AREA)
Description
本発明は、判定対象物を撮像して取得した多値画像及び撮像手段からの距離を濃淡値に変換した距離画像を、基準となる多値画像及び距離画像と比較する画像処理装置、画像処理方法、及びコンピュータプログラムに関する。 The present invention relates to an image processing apparatus and image processing for comparing a multi-value image obtained by imaging a determination target and a distance image obtained by converting a distance from an imaging unit into a gray value with a reference multi-value image and a distance image. The present invention relates to a method and a computer program.
従来、判定対象物を撮像して取得した多値画像を、基準となる多値画像と単純に比較することにより、判定対象物が良品であるか否かを判定する方法が開発されている。しかし、両画像を単純に比較するだけでは、良品自体の形状ばらつき、計算誤差、ノイズ等による画素値変動を考慮することができず、良品であるにもかかわらず良品ではないと誤って判定する等、正しく良品判定を行うことができないおそれがある等の問題があった。 Conventionally, a method has been developed for determining whether or not a determination target object is a non-defective product by simply comparing a multi-value image acquired by imaging the determination target object with a reference multi-value image. However, by simply comparing the two images, it is not possible to take into account pixel value fluctuations due to shape variations, calculation errors, noise, etc. of the non-defective product itself, and erroneously determine that it is non-defective even though it is non-defective For example, there is a possibility that correct product determination may not be performed correctly.
斯かる問題を解決するために、例えば特許文献1では、複数の良品に関する多値画像を準備しておき、それらの平均画像、標準偏差画像を求めることで、安定した良品判定を可能としている。具体的には、複数の良品に関する多値画像について位置合わせを行った後、同一座標の画素ごとに画素値の平均値及び標準偏差を算出する。判定対象物の多値画像について良品に関する多値(平均)画像との間で位置合わせを行った後、同一座標の画素ごとに、平均値との差分値を算出し、標準偏差に基づいて定義した閾値画像と画素ごとに比較することにより、良品判定を行っている。
In order to solve such a problem, for example, in
このようにすることで、良品自体の形状ばらつき、計算誤差、ノイズ等による画素値変動分を画素ごとの画素値のばらつき度合に応じて効果的に排除することができるとともに、画素ごとに良品であると判定する範囲が相違することから、良品の部位ごとに適切な良品判定の閾値を設定することができる。したがって、より高い精度で良品判定を行うことが可能となっている。 In this way, pixel value fluctuations due to non-defective product shape variations, calculation errors, noise, etc. can be effectively eliminated according to the degree of pixel value variation for each pixel, and non-defective products can be obtained for each pixel. Since the range to be determined is different, an appropriate non-defective threshold value can be set for each non-defective part. Therefore, the non-defective product can be determined with higher accuracy.
しかし、特許文献1に開示されている画像処理方法では、濃淡値の差異として表れにくい判定対象物表面の凹凸を検出することができない。したがって、色の相違、表面の凹凸の相違が良品の分布範囲に含まれているにもかかわらず、判定対象物を載置した位置のばらつき等により光の当たり方が変化し、濃淡値の差異が実際より大きく検出されることで、良品ではないと誤って判定されるおそれがあった。
However, the image processing method disclosed in
また、三次元形状の差異を判断するべく、判定対象物の三次元形状を測定した三次元データと、判定対象物の三次元CADデータとを比較することも考えられる。しかし、実測値として正確な三次元データを取得する必要があり、多くの方向から計測した複数の計測データで三次元データを構成する等、良品であるか否かを判定するまでに相当の時間を要し、しかも装置自体が高価となる。 It is also conceivable to compare the three-dimensional data obtained by measuring the three-dimensional shape of the determination target with the three-dimensional CAD data of the determination target in order to determine the difference in the three-dimensional shape. However, it is necessary to acquire accurate three-dimensional data as measured values, and it takes a considerable amount of time to determine whether or not the product is non-defective, such as configuring the three-dimensional data with a plurality of measurement data measured from many directions. And the device itself is expensive.
さらに、欠陥が輪郭線の近傍に生じている場合、輪郭線近傍では画素値の標準偏差が大きいことから欠陥として認識することが困難であり、欠陥の存在を検出する感度が低くなってしまうという問題点があった。輪郭線の近傍にて画素値の標準偏差が大きくなるのは、判定対象物を載置した位置のばらつき、撮像位置の相違、良品自体の形状ばらつき等に起因する。 Furthermore, when a defect occurs in the vicinity of the contour line, it is difficult to recognize the defect as the standard deviation of the pixel value is large in the vicinity of the contour line, and the sensitivity to detect the presence of the defect is low. There was a problem. The increase in the standard deviation of the pixel value in the vicinity of the contour line is caused by the variation in the position where the determination target is placed, the difference in the imaging position, the shape variation of the non-defective product itself, and the like.
本発明は上記事情に鑑みてなされたものであり、第二の多値画像又は第二の距離画像の性状に応じて適切な特徴量を選択することができ、欠陥の存在だけでなく、色の相違についても高い感度で検出することができる画像処理装置、画像処理方法、及びコンピュータプログラムを提供することを目的とする。 The present invention has been made in view of the above circumstances, and an appropriate feature amount can be selected according to the properties of the second multi-valued image or the second distance image. An object of the present invention is to provide an image processing apparatus, an image processing method, and a computer program that can detect the difference between the two with high sensitivity .
上記目的を達成するために第1発明に係る画像処理装置は、判定対象物を撮像した多値画像を、良品に関する多値画像群と比較して良品判定を行う画像処理装置において、撮像手段で撮像された、良品に関する複数の第一の多値画像を取得する良品画像取得手段と、取得した第一の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第一の距離画像を生成する良品距離画像生成手段と、前記第一の多値画像の画素ごとに、色成分に関する画素値を特徴量として抽出する第一の特徴量抽出手段と、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する第二の特徴量抽出手段とを少なくとも含む複数の特徴量抽出手段のうち、少なくとも1つの特徴量抽出手段の選択を受け付ける選択受付手段と、取得した第一の多値画像及び生成した前記第一の距離画像の少なくとも一方の画素ごとの濃淡値の、良品判定を行うための分布範囲を算出する分布範囲算出手段と、判定対象物に関する第二の多値画像を取得する多値画像取得手段と、取得した第二の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第二の距離画像を生成する対象距離画像生成手段と、取得した第二の多値画像及び生成した前記第二の距離画像の少なくとも一方の画素ごとの濃淡値が対応する前記分布範囲に含まれているか否かを判断する判断手段とを備え、前記分布範囲算出手段は、前記第一の特徴量抽出手段の選択を受け付けた場合には第一の分布範囲を、前記第二の特徴量抽出手段の選択を受け付けた場合には第二の分布範囲を、それぞれ選択を受け付けた特徴量抽出手段により抽出した特徴量に基づいて算出するようにしてあり、前記判断手段は、前記第二の多値画像又は前記第二の距離画像の画素ごとに、選択を受け付けた特徴量抽出手段により特徴量を抽出し、抽出された特徴量に対応する前記第一の分布範囲又は前記第二の分布範囲に含まれているか否かを判断するようにしてあることを特徴とする。 In order to achieve the above object, an image processing apparatus according to a first aspect of the present invention is an image processing apparatus for performing non-defective product determination by comparing a multi-value image obtained by capturing an object to be determined with a multi-value image group related to non-defective products. A first non-defective image acquisition unit that acquires a plurality of first multi-valued images related to the acquired non-defective product, and a first value obtained by converting the distance from the imaging unit into a gray value for each pixel of the acquired first multi-level image. A non-defective distance image generating unit that generates a distance image of the first multi-value image, a first feature amount extracting unit that extracts a pixel value related to a color component as a feature amount for each pixel of the first multi-value image, and the first multi-value image At least one feature among a plurality of feature quantity extraction means including at least a second feature quantity extraction means for extracting edge strengths in two different directions as feature quantities for each pixel of the value image or the first distance image Selection of quantity extraction means A selection accepting means for accepting the acquired first gray value of at least every one of the pixels of the multivalued image and the generated the first distance image, the distribution range calculation means for calculating a distribution range for performing good decision A second multi-valued image obtaining unit for obtaining a second multi-valued image relating to the determination target, and a second value obtained by converting the distance from the imaging unit into a gray value for each pixel of the obtained second multi-valued image. Whether or not the target distance image generating means for generating the distance image, the acquired second multi-value image, and the gray value for each pixel of the generated second distance image are included in the corresponding distribution range. Determination means for determining whether or not the distribution range calculation means selects the first distribution range when the selection of the first feature quantity extraction means is accepted, and selects the second feature quantity extraction means. If you accept the second The distribution range is calculated based on the feature amount extracted by the feature amount extraction unit that has received the selection, and the determination unit is configured for each pixel of the second multi-value image or the second distance image. In addition, a feature amount is extracted by the feature amount extraction unit that has received the selection, and it is determined whether the feature amount is included in the first distribution range or the second distribution range corresponding to the extracted feature amount. It is characterized by being.
また、第2発明に係る画像処理装置は、第1発明において、複数の前記撮像手段を備え、前記良品距離画像生成手段又は前記対象距離画像生成手段は、複数の前記撮像手段で撮像された複数の前記第一の多値画像又は前記第二の多値画像に基づいて前記第一の距離画像又は前記第二の距離画像を生成するようにしてあることを特徴とする。 An image processing apparatus according to a second aspect of the present invention includes the plurality of imaging units in the first invention, and the non-defective distance image generation unit or the target distance image generation unit includes a plurality of images captured by the plurality of imaging units. The first distance image or the second distance image is generated based on the first multi-value image or the second multi-value image.
また、第3発明に係る画像処理装置は、第1発明において、前記撮像手段と、判定対象物に対して光を照射する光照射手段とを備え、前記良品距離画像生成手段又は前記対象距離画像生成手段は、前記撮像手段で撮像された前記第一の多値画像又は前記第二の多値画像、及び照射した光に基づいて前記第一の距離画像又は前記第二の距離画像を生成するようにしてあることを特徴とする。 An image processing apparatus according to a third aspect of the present invention is the image processing apparatus according to the first aspect, further comprising the imaging unit and a light irradiation unit that irradiates a determination target with light, and the good product distance image generation unit or the target distance image. The generation unit generates the first distance image or the second distance image based on the first multi-value image or the second multi-value image captured by the imaging unit and the irradiated light. It is characterized by the above.
また、第4発明に係る画像処理装置は、第1乃至第3発明のいずれか1つにおいて、前記第一の特徴量抽出手段の選択を受け付けた場合、前記分布範囲算出手段は、前記第一の多値画像の画素ごとに、色成分に関する画素値の平均値を算出する画素平均値算出手段と、算出した色成分に関する画素値の平均値に基づいて、前記第一の多値画像の画素ごとに、色成分それぞれを座標軸とする多次元空間における画素値の前記第一の分布範囲を算出する画素値分布範囲算出手段とを有し、前記判断手段は、前記第二の多値画像の画素ごとに、色成分に関する画素値を算出して、色成分それぞれを座標軸とする多次元空間における画素値の前記第一の分布範囲に含まれているか否かを判断するようにしてあり、前記第二の特徴量抽出手段の選択を受け付けた場合、前記分布範囲算出手段は、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を算出する方向別エッジ強度算出手段と、算出したエッジ強度に基づいて、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度の平均値を算出するエッジ強度平均値算出手段と、算出した平均値を中心として、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度の相互相関分布範囲を前記第二の分布範囲として算出するエッジ強度分布範囲算出手段とを有し、前記判断手段は、前記第二の多値画像又は前記第二の距離画像の画素ごとに、異なる二方向におけるエッジ強度を算出して、算出したエッジ強度が対応する前記相互相関分布範囲に含まれているか否かを判断するようにしてあることを特徴とする。 The image processing apparatus according to a fourth aspect of the present invention is the image processing apparatus according to any one of the first to third aspects, wherein when the selection of the first feature amount extraction unit is received, the distribution range calculation unit Pixel average value calculating means for calculating an average value of pixel values relating to color components for each pixel of the multi-value image, and pixels of the first multi-value image based on the calculated average values of pixel values relating to color components Each having a pixel value distribution range calculation means for calculating the first distribution range of pixel values in a multi-dimensional space having each color component as a coordinate axis, and the determination means includes the second multi-value image. For each pixel, a pixel value related to the color component is calculated, and it is determined whether or not the pixel value is included in the first distribution range of the pixel value in the multidimensional space with each color component as a coordinate axis. Receives selection of second feature extraction means When attached, the distribution range calculation means includes edge strength calculation means for each direction for calculating edge strength in two different directions for each pixel of the first multi-value image or the first distance image, and the calculated edge Based on the intensity, edge intensity average value calculating means for calculating an average value of edge intensity in two different directions for each pixel of the first multi-value image or the first distance image, and centering the calculated average value An edge intensity distribution range calculation means for calculating a cross-correlation distribution range of edge strengths in two different directions as the second distribution range for each pixel of the first multi-valued image or the first distance image. And the determination means calculates edge strengths in two different directions for each pixel of the second multi-valued image or the second distance image, and the calculated edge strengths correspond to the corresponding phases. Characterized in that you have to determine whether or not included in the correlation distribution range.
また、第5発明に係る画像処理装置は、第4発明において、前記方向別エッジ強度算出手段は、前記エッジ強度を互いに直交する二方向にて算出するようにしてあることを特徴とする。 An image processing apparatus according to a fifth aspect of the present invention is the image processing apparatus according to the fourth aspect , wherein the direction-specific edge strength calculating means calculates the edge strength in two directions orthogonal to each other.
また、第6発明に係る画像処理装置は、第5発明において、前記方向別エッジ強度算出手段は、前記エッジ強度を二次元画像における行方向と列方向との二方向にて算出するようにしてあることを特徴とする。 The image processing apparatus according to a sixth aspect of the present invention is the image processing apparatus according to the fifth aspect , wherein the direction-specific edge strength calculating means calculates the edge strength in two directions of a row direction and a column direction in a two-dimensional image. It is characterized by being.
また、第7発明に係る画像処理装置は、第4乃至第6発明のいずれか1つにおいて、前記エッジ強度分布範囲算出手段は、仮想的な楕円領域として前記相互相関分布範囲を前記第二の分布範囲として算出するようにしてあることを特徴とする。 An image processing apparatus according to a seventh invention is the image processing apparatus according to any one of the fourth to sixth inventions, wherein the edge intensity distribution range calculation means sets the cross-correlation distribution range as the virtual ellipse region to the second It is characterized in that it is calculated as a distribution range .
また、第8発明に係る画像処理装置は、第1乃至第7発明のいずれか1つにおいて、前記第一の多値画像と前記第二の多値画像との位置合わせ又は前記第一の距離画像と前記第二の距離画像との位置合わせを行う位置調整手段を備えることを特徴とする。 According to an eighth aspect of the present invention, in any one of the first to seventh aspects, the first multi-value image and the second multi-value image are aligned or the first distance is set. Position adjustment means for aligning the image with the second distance image is provided.
また、第9発明に係る画像処理装置は、第1乃至第8発明のいずれか1つにおいて、前記選択受付手段が、複数の特徴量抽出手段の選択を受け付けた場合であって、前記判断手段が、選択を受け付けた特徴量抽出手段で抽出された特徴量のうち、いずれかの特徴量がそれぞれの特徴量に対応する前記第一の分布範囲又は前記第二の分布範囲に含まれていないと判断したとき、前記判定対象物は良品ではないと判定するようにしてあることを特徴とする。 An image processing apparatus according to a ninth invention is the image processing apparatus according to any one of the first to eighth inventions, wherein the selection accepting means accepts selection of a plurality of feature quantity extracting means. However, any one of the feature quantities extracted by the feature quantity extraction unit that has received the selection is not included in the first distribution range or the second distribution range corresponding to each feature quantity. When the determination is made, it is determined that the determination object is not a non-defective product.
次に、上記目的を達成するために第10発明に係る画像処理方法は、判定対象物を撮像した多値画像を、良品に関する多値画像群と比較して良品判定を行う画像処理装置で実行することが可能な画像処理方法において、撮像手段で撮像された、良品に関する複数の第一の多値画像を取得するステップと、取得した第一の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第一の距離画像を生成するステップと、前記第一の多値画像の画素ごとに、色成分に関する画素値を特徴量として抽出する第一の特徴量抽出手段と、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する第二の特徴量抽出手段とを少なくとも含む複数の特徴量抽出手段のうち、少なくとも1つの特徴量抽出手段の選択を受け付けるステップと、取得した第一の多値画像及び生成した前記第一の距離画像の少なくとも一方の画素ごとの濃淡値の、良品判定を行うための分布範囲を算出するステップと、判定対象物に関する第二の多値画像を取得するステップと、取得した第二の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第二の距離画像を生成するステップと、取得した第二の多値画像及び生成した前記第二の距離画像の少なくとも一方の画素ごとの濃淡値が対応する前記分布範囲に含まれているか否かを判断するステップとを含み、前記第一の特徴量抽出手段の選択を受け付けた場合には第一の分布範囲を、前記第二の特徴量抽出手段の選択を受け付けた場合には第二の分布範囲を、それぞれ選択を受け付けた特徴量抽出手段により抽出した特徴量に基づいて算出し、前記第二の多値画像又は前記第二の距離画像の画素ごとに、選択を受け付けた特徴量抽出手段により特徴量を抽出し、抽出された特徴量に対応する前記第一の分布範囲又は前記第二の分布範囲に含まれているか否かを判断することを特徴とする。 Next, in order to achieve the above object, the image processing method according to the tenth aspect of the present invention is executed by an image processing apparatus that performs a non-defective product determination by comparing a multi-value image obtained by imaging a determination target with a multi-value image group related to a good product In the image processing method that can be performed, a step of acquiring a plurality of first multi-value images related to non-defective products imaged by the imaging unit, and for each pixel of the acquired first multi-value image, from the imaging unit Generating a first distance image obtained by converting the distance into a gray value, and a first feature amount extraction means for extracting a pixel value related to a color component as a feature amount for each pixel of the first multi-value image; A plurality of feature quantity extraction means including at least a second feature quantity extraction means for extracting edge strengths in two different directions as feature quantities for each pixel of the first multi-value image or the first distance image. At least 1 Calculating a step of accepting a selection of a feature extraction means, the gray value of each at least one of the pixels of the first multivalued image, and generated the first distance image acquired, the distribution range for performing good determination of A second distance image obtained by converting a distance from the imaging unit into a gray value for each pixel of the acquired second multi-value image And determining whether or not a gray value for each pixel of at least one pixel of the acquired second multi-valued image and the generated second distance image is included in the corresponding distribution range; Including a first distribution range when the selection of the first feature amount extraction unit is received, and a second distribution range when the selection of the second feature amount extraction unit is received, respectively. Accept selection Calculated based on the feature amount extracted by the feature amount extraction means, and for each pixel of the second multi-valued image or the second distance image, the feature amount extraction means that has received a selection extracts the feature amount, It is characterized by determining whether it is included in said 1st distribution range or said 2nd distribution range corresponding to the extracted feature-value.
次に、上記目的を達成するために第11発明に係るコンピュータプログラムは、判定対象物を撮像した多値画像を、良品に関する多値画像群と比較して良品判定を行う画像処理装置で実行することが可能なコンピュータプログラムにおいて、前記画像処理装置を、撮像手段で撮像された、良品に関する複数の第一の多値画像を取得する良品画像取得手段、取得した第一の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第一の距離画像を生成する良品距離画像生成手段、前記第一の多値画像の画素ごとに、色成分に関する画素値を特徴量として抽出する第一の特徴量抽出手段と、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する第二の特徴量抽出手段とを少なくとも含む複数の特徴量抽出手段のうち、少なくとも1つの特徴量抽出手段の選択を受け付ける選択受付手段、取得した第一の多値画像及び生成した前記第一の距離画像の少なくとも一方の画素ごとの濃淡値の、良品判定を行うための分布範囲を算出する分布範囲算出手段、判定対象物に関する第二の多値画像を取得する多値画像取得手段、取得した第二の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第二の距離画像を生成する対象距離画像生成手段、及び取得した第二の多値画像及び生成した前記第二の距離画像の少なくとも一方の画素ごとの濃淡値が対応する前記分布範囲に含まれているか否かを判断する判断手段として機能させ、前記分布範囲算出手段を、前記第一の特徴量抽出手段の選択を受け付けた場合には第一の分布範囲を、前記第二の特徴量抽出手段の選択を受け付けた場合には第二の分布範囲を、それぞれ選択を受け付けた特徴量抽出手段により抽出した特徴量に基づいて算出する手段として機能させ、前記判断手段を、前記第二の多値画像又は前記第二の距離画像の画素ごとに、選択を受け付けた特徴量抽出手段により特徴量を抽出し、抽出された特徴量に対応する前記第一の分布範囲又は前記第二の分布範囲に含まれているか否かを判断する手段として機能させることを特徴とする。 Next, in order to achieve the above object, a computer program according to an eleventh aspect of the invention executes an image processing apparatus that performs a non-defective product determination by comparing a multi-value image obtained by imaging a determination target with a multi-value image group related to a non-defective product. A non-defective image acquiring unit that acquires a plurality of first multi-valued images related to a non-defective product that has been imaged by the imaging unit in the computer program capable of being acquired, for each pixel of the acquired first multi-valued image In addition, a non-defective distance image generation unit that generates a first distance image obtained by converting the distance from the imaging unit into a gray value, and extracts a pixel value related to a color component as a feature amount for each pixel of the first multi-value image First feature quantity extracting means, and second feature quantity extracting means for extracting edge strengths in two different directions as feature quantities for each pixel of the first multi-value image or the first distance image. Among the plurality of feature extraction means including at least a selection accepting means for accepting a selection of at least one feature extraction means, for each of at least one of the pixels of the first multivalued image, and generated the first distance image acquired Distribution range calculation means for calculating a distribution range for performing non-defective product determination of gray value, multi-value image acquisition means for acquiring a second multi-value image relating to the determination target, and for each pixel of the acquired second multi-value image In addition, target distance image generation means for generating a second distance image obtained by converting the distance from the imaging means into a gray value, and at least one of the acquired second multi-value image and the generated second distance image It functions as a determination unit that determines whether or not the gray value for each pixel is included in the corresponding distribution range, and the distribution range calculation unit accepts the selection of the first feature amount extraction unit If the selection of the second feature amount extraction unit is accepted, the second distribution range is selected based on the feature amount extracted by the feature amount extraction unit that has accepted the selection. It is made to function as a means for calculating, and the determination means extracts the feature amount by the feature amount extraction means that accepts selection for each pixel of the second multi-value image or the second distance image, and the extracted feature It is made to function as a means to judge whether it is contained in said 1st distribution range or said 2nd distribution range corresponding to quantity.
第1発明、第10発明及び第11発明では、撮像手段で撮像された、良品に関する複数の第一の多値画像を取得し、取得した第一の多値画像の画素ごとに、撮像手段からの距離を濃淡値に変換した第一の距離画像を生成する。第一の多値画像の画素ごとに、色成分に関する画素値を特徴量として抽出する第一の特徴量抽出手段と、第一の多値画像又は第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する第二の特徴量抽出手段とを少なくとも含む複数の特徴量抽出手段のうち、少なくとも1つの特徴量抽出手段の選択を受け付ける。取得した第一の多値画像及び生成した第一の距離画像の少なくとも一方の画素ごとの濃淡値の、良品判定を行うための分布範囲を算出する。判定対象物に関する第二の多値画像を取得し、取得した第二の多値画像の画素ごとに、撮像手段からの距離を濃淡値に変換した第二の距離画像を生成し、取得した第二の多値画像及び生成した第二の距離画像の少なくとも一方の画素ごとの濃淡値が対応する前記分布範囲に含まれているか否かを判断する。第二の距離画像に基づいて、良品であるか否かを判断するので、従来の方法では検出することが困難である判定対象物表面の凹凸を検出することができ、判定対象物を載置した位置のばらつき等により光の当たり方が変化した場合であっても、良品ではないと誤って判定される可能性を低減することが可能となる。また、第一の特徴量抽出手段の選択を受け付けた場合には第一の分布範囲を、第二の特徴量抽出手段の選択を受け付けた場合には第二の分布範囲を、それぞれ選択を受け付けた特徴量抽出手段により抽出した特徴量に基づいて算出し、第二の多値画像又は第二の距離画像の画素ごとに、選択を受け付けた特徴量抽出手段により特徴量を抽出し、抽出された特徴量に対応する第一の分布範囲又は第二の分布範囲に含まれているか否かを判断するので、第二の多値画像又は第二の距離画像の性状に応じて、特徴量抽出手段を選択することができ、適切な特徴量を選択することにより、欠陥が輪郭線近傍に生じている場合であっても欠陥の存在を検出することができ、色の相違についても高い感度で検出することが可能となる。例えば選択を受け付けた特徴量が色成分に関する画素値である場合、色の相違による濃淡のばらつきの存在を検出することができ、選択を受け付けた特徴量が異なる二方向におけるエッジ強度である場合、輪郭線近傍の欠陥の存在を検出することができる。 In the first invention, the tenth invention, and the eleventh invention, a plurality of first multi-value images related to non-defective products, which are imaged by the imaging means, are acquired, and the pixels of the acquired first multi-value image are obtained from the imaging means. A first distance image is generated by converting the distance to a gray value. The first feature value extraction means for extracting the pixel value related to the color component as a feature value for each pixel of the first multivalued image and two different pixels for each pixel of the first multivalued image or the first distance image. The selection of at least one feature quantity extraction means is accepted from among a plurality of feature quantity extraction means including at least a second feature quantity extraction means for extracting edge strength in the direction as a feature quantity. A distribution range for performing non-defective product determination of the gray value for each pixel of at least one of the acquired first multi-value image and the generated first distance image is calculated. A second multi-value image related to the determination target is acquired, and for each pixel of the acquired second multi-value image, a second distance image in which the distance from the imaging unit is converted into a gray value is generated and acquired. It is determined whether the gray value for each pixel of at least one of the second multi-valued image and the generated second distance image is included in the corresponding distribution range. Since it is determined based on the second distance image whether or not it is a non-defective product, it is possible to detect unevenness on the surface of the determination object that is difficult to detect by the conventional method, and to place the determination object Even if the way the light strikes changes due to variations in the position, etc., it is possible to reduce the possibility of erroneously determining that the light is not a good product. When the selection of the first feature quantity extraction unit is accepted, the first distribution range is accepted, and when the selection of the second feature quantity extraction unit is accepted, the second distribution range is accepted. Calculated based on the feature quantity extracted by the feature quantity extraction means, and extracted for each pixel of the second multi-valued image or the second distance image by extracting the feature quantity by the feature quantity extraction means that accepted the selection. Since it is determined whether or not it is included in the first distribution range or the second distribution range corresponding to the feature amount, the feature amount extraction is performed according to the property of the second multi-value image or the second distance image. By selecting an appropriate feature amount, it is possible to detect the presence of a defect even when the defect occurs in the vicinity of the contour line, and with high sensitivity for color differences. It becomes possible to detect. For example, when the feature amount that has received the selection is a pixel value related to a color component, it is possible to detect the presence or absence of shading variation due to a difference in color, and when the feature amount that has received the selection is an edge strength in two different directions, The presence of a defect near the contour line can be detected.
第2発明では、複数の撮像手段を備え、複数の撮像手段で撮像された複数の第一の多値画像又は第二の多値画像に基づいて第一の距離画像又は第二の距離画像を生成する。複数の撮像手段を備えることにより、複数の多値画像を取得することができ、ステレオ計測法等のいわゆるパッシブ方式にて距離画像を生成することができる。 In the second invention, the first distance image or the second distance image is provided on the basis of the plurality of first multi-value images or second multi-value images that are provided with a plurality of image pickup means and are imaged by the plurality of image pickup means. Generate. By providing a plurality of imaging means, a plurality of multi-value images can be acquired, and a distance image can be generated by a so-called passive method such as a stereo measurement method.
第3発明では、撮像手段と、判定対象物に対して光を照射する光照射手段とを備え、撮像手段で撮像された第一の多値画像又は第二の多値画像、及び照射した光に基づいて第一の距離画像又は第二の距離画像を生成する。これにより、高さ方向の距離を計測するために光を能動的に照射することができ、光切断法、パターン投影法等のいわゆるアクティブ方式にて距離画像を生成することができる。 In the third invention, the first multi-valued image or the second multi-valued image picked up by the image pickup means and the light irradiation means for irradiating the determination target with light, and the irradiated light The first distance image or the second distance image is generated based on the above. Accordingly, light can be actively irradiated to measure the distance in the height direction, and a distance image can be generated by a so-called active method such as a light cutting method or a pattern projection method.
第4発明では、第一の特徴量抽出手段の選択を受け付けた場合、第一の多値画像の画素ごとに、色成分に関する画素値の平均値を算出し、算出した色成分に関する画素値の平均値に基づいて、第一の多値画像の画素ごとに、色成分それぞれを座標軸とする多次元空間における画素値の第一の分布範囲を算出する。第二の多値画像の画素ごとに、色成分に関する画素値を算出して、色成分それぞれを座標軸とする多次元空間における画素値の第一の分布範囲に含まれているか否かを判断する。第二の特徴量抽出手段の選択を受け付けた場合、第一の多値画像又は第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を算出し、算出したエッジ強度に基づいて、第一の多値画像又は第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度の平均値を算出し、算出した平均値を中心として、第一の多値画像又は第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度の相互相関分布範囲を第二の分布範囲として算出する。第二の多値画像又は第二の距離画像の画素ごとに、異なる二方向におけるエッジ強度を算出して、算出したエッジ強度が対応する相互相関分布範囲に含まれているか否かを判断する。これにより、選択を受け付けた特徴量が色成分に関する画素値である場合、算出した画素値の第一の分布範囲に含まれているか否かにより色の相違を検出することができ、選択を受け付けた特徴量が異なる二方向におけるエッジ強度である場合、算出したエッジ強度の相互相関分布範囲に含まれているか否かにより、輪郭線近傍の欠陥の存在を検出することができる。 In the fourth aspect of the invention, when the selection of the first feature amount extraction unit is accepted, the average value of the pixel values related to the color component is calculated for each pixel of the first multi-value image, and the pixel value related to the calculated color component is calculated. Based on the average value, a first distribution range of pixel values in a multidimensional space with each color component as a coordinate axis is calculated for each pixel of the first multivalued image. For each pixel of the second multi-valued image, a pixel value related to the color component is calculated, and it is determined whether or not the pixel value is included in the first distribution range of the pixel value in the multi-dimensional space with each color component as a coordinate axis. . When the selection of the second feature amount extraction means is accepted, the edge strength in two different directions is calculated for each pixel of the first multi-value image or the first distance image, and the first edge value is calculated based on the calculated edge strength. For each pixel of one multi-value image or first distance image, an average value of edge strengths in two different directions is calculated, and the first multi-value image or first distance image is calculated around the calculated average value. For each pixel, a cross-correlation distribution range of edge strengths in two different directions is calculated as a second distribution range . For each pixel of the second multi-value image or the second distance image, the edge strength in two different directions is calculated, and it is determined whether or not the calculated edge strength is included in the corresponding cross-correlation distribution range. As a result, when the feature quantity accepted for selection is a pixel value related to a color component, it is possible to detect a color difference depending on whether the calculated pixel value is included in the first distribution range, and accept the selection. If the detected feature amounts are edge strengths in two different directions, it is possible to detect the presence of a defect near the contour line depending on whether or not it is included in the cross-correlation distribution range of the calculated edge strengths.
第5発明では、エッジ強度を互いに直交する二方向にて算出することにより、相互相関分布範囲の算出が容易となり、演算処理負荷を軽減することが可能となる。 In the fifth invention, by calculating the edge strength in two directions orthogonal to each other, the calculation of the cross-correlation distribution range is facilitated, and the processing load can be reduced.
第6発明では、エッジ強度を二次元画像における行方向と列方向との二方向にて算出することにより、さらに相互相関分布範囲の算出が容易となり、演算処理負荷をより軽減することが可能となる。 In the sixth invention, by calculating the edge strength in two directions of the row direction and the column direction in the two-dimensional image, the calculation of the cross-correlation distribution range is further facilitated, and the processing load can be further reduced. Become.
第7発明では、仮想的な楕円領域として相互相関分布範囲を第二の分布範囲として算出することにより、第二の多値画像又は第二の距離画像の画素ごとに算出した異なる二方向におけるエッジ強度が、相互相関分布範囲に含まれているか否かをより容易に判断することが可能となる。 In the seventh invention, by calculating the cross-correlation distribution range as the second distribution range as a virtual elliptical region, edges in different two directions calculated for each pixel of the second multi-valued image or the second distance image It is possible to more easily determine whether or not the intensity is included in the cross-correlation distribution range.
第8発明では、第一の多値画像と第二の多値画像との位置合わせ又は第一の距離画像と第二の距離画像との位置合わせを行うことにより、第二の多値画像又は第二の距離画像の画素ごとに算出した異なる二方向におけるエッジ強度が、相互相関分布範囲に含まれているか否かをより正確に判断することが可能となる。 In the eighth invention, the second multi-value image or the second multi-value image is aligned by aligning the first multi-value image and the second multi-value image, or the first distance image and the second distance image. It is possible to more accurately determine whether or not the edge intensity in two different directions calculated for each pixel of the second distance image is included in the cross-correlation distribution range.
第9発明では、複数の特徴量抽出手段の選択を受け付けた場合であって、選択を受け付けた特徴量抽出手段で抽出された特徴量のうち、いずれかの特徴量がそれぞれの特徴量に対応する第一の分布範囲又は第二の分布範囲に含まれていないと判断したとき、判定対象物は良品ではないと判定する。これにより、従来は不適切な特徴量に基づいて、本来良品ではない判定対象物を誤って良品であると判定していた場合であっても、適切な特徴量に基づいて良品ではないと正しく判定することが可能となる。 In the ninth invention, when selection of a plurality of feature quantity extraction means is accepted, one of the feature quantities extracted by the feature quantity extraction means that has accepted the selection corresponds to each feature quantity. When it is determined that it is not included in the first distribution range or the second distribution range , it is determined that the determination target object is not a non-defective product. As a result, even if it was determined that a non-defective product was determined to be a non-defective product based on an inappropriate feature value, it was correctly determined that it was not a non-defective product based on an appropriate feature value. It becomes possible to judge.
本発明によれば、第二の距離画像に基づいて、良品であるか否かを判断するので、従来の方法では検出することが困難である判定対象物表面の凹凸を検出することができ、判定対象物を載置した位置のばらつき等により光の当たり方が変化した場合であっても、良品ではないと誤って判定される可能性を低減することが可能となる。また、第二の多値画像又は第二の距離画像の性状に応じて、特徴量抽出手段を選択することができ、適切な特徴量を選択することにより、欠陥が輪郭線近傍に生じている場合であっても欠陥の存在を検出することができ、色の相違についても高い感度で検出することが可能となる。例えば選択を受け付けた特徴量が色成分に関する画素値である場合、色の相違による濃淡のばらつきの存在を検出することができ、選択を受け付けた特徴量が異なる二方向におけるエッジ強度である場合、輪郭線近傍の欠陥の存在を検出することができる。 According to the present invention, since it is determined whether or not it is a non-defective product based on the second distance image, it is possible to detect unevenness on the surface of the determination target that is difficult to detect by the conventional method, Even if the way the light hits changes due to variations in the position where the determination target is placed, it is possible to reduce the possibility that it is erroneously determined that it is not a non-defective product. In addition, the feature quantity extraction unit can be selected according to the properties of the second multi-value image or the second distance image, and defects are generated in the vicinity of the contour line by selecting an appropriate feature quantity. Even in this case, the presence of a defect can be detected, and a difference in color can be detected with high sensitivity. For example, when the feature amount that has received the selection is a pixel value related to a color component, it is possible to detect the presence or absence of shading variation due to a difference in color, and when the feature amount that has received the selection is an edge strength in two different directions, The presence of a defect near the contour line can be detected.
以下、本発明の実施の形態に係る画像処理装置について、図面を参照して説明する。なお、参照する図面を通じて、同一又は同様の構成又は機能を有する要素については、同一又は同様の符号を付して、詳細な説明を省略する。 Hereinafter, an image processing apparatus according to an embodiment of the present invention will be described with reference to the drawings. Throughout the drawings to be referred, elements having the same or similar configuration or function are denoted by the same or similar reference numerals, and detailed description thereof is omitted.
(実施の形態1)
図1は、本発明の実施の形態1に係る画像処理装置の構成を模式的に示すブロック図である。図1に示すように本実施の形態1に係る画像処理装置2は、判定対象物又は良品の多値画像を撮像する撮像手段であるカメラ1、判定対象物又は良品にレーザ光を照射するライト9及び撮像された多値画像、生成された距離画像又は演算処理の途上で生成された画像を表示する画像表示手段である表示装置3に接続されている。
(Embodiment 1)
FIG. 1 is a block diagram schematically showing a configuration of an image processing apparatus according to
画像処理装置2は、少なくともCPU(中央演算装置)、LSI等で構成された主制御部21、メモリ22、記憶手段23、入力手段24、出力手段25、通信手段26、補助記憶手段27及び上述したハードウェアを接続する内部バス28で構成されている。主制御部21は、内部バス28を介して画像処理装置2の上述したようなハードウェア各部と接続されており、上述したハードウェア各部の動作を制御するとともに、記憶手段23に記憶されているコンピュータプログラム5に従って、種々のソフトウェア的機能を実行する。メモリ22は、SRAM、SDRAM等の揮発性メモリで構成され、コンピュータプログラム5の実行時にロードモジュールが展開され、コンピュータプログラム5の実行時に発生する一時的なデータ等を記憶する。
The
記憶手段23は、内蔵される固定型記憶装置(ハードディスク、フラッシュメモリ)、ROM等で構成されている。記憶手段23に記憶されているコンピュータプログラム5は、プログラム及びデータ等の情報を記録したDVD、CD−ROM、フラッシュメモリ等の可搬型記録媒体4から、補助記憶手段27によりダウンロードされ、実行時には記憶手段23からメモリ22へ展開して実行される。もちろん、通信手段26を介して外部のコンピュータからダウンロードされたコンピュータプログラムであっても良い。
The storage means 23 is composed of a built-in fixed storage device (hard disk, flash memory), ROM or the like. The
記憶手段23は、取得した良品に関する複数の多値画像に関する多値画像データを記憶しておく良品画像データ記憶部231と、良品画像データ記憶部231に記憶してある良品に関する複数の多値画像データに基づいて、画素ごとに濃淡値の分布範囲に関する情報(分布範囲情報)を記憶する分布範囲情報記憶部232とを備えている。判定対象物の特徴量が、分布範囲情報記憶部232に記憶してある分布範囲情報に基づいて特定される分布範囲に含まれているか否かの判断に応じて、判定対象物が良品であるか否かを判定する。
The
通信手段26は内部バス28に接続されており、インターネット、LAN、WAN等の外部のネットワークに接続されることにより、外部のコンピュータ等とデータ送受信を行うことが可能となっている。すなわち、上述した記憶手段23は、画像処理装置2に内蔵される構成に限定されるものではなく、通信手段26を介して接続されている外部のサーバコンピュータ等に設置されているハードディスク等の外部記録媒体であっても良い。
The communication means 26 is connected to an
入力手段24は、キーボード及びマウス等のデータ入力媒体の他、液晶パネル等と一体となったタッチパネル等の入力情報を取得する装置全般を含む広い概念である。出力手段25は、レーザプリンタ、ドットプリンタ等の印刷装置等を意味する。 The input means 24 is a broad concept including all devices for acquiring input information such as a touch panel integrated with a liquid crystal panel or the like in addition to a data input medium such as a keyboard and a mouse. The output means 25 means a printing device such as a laser printer or a dot printer.
カメラ(撮像手段)1は、CCD撮像素子を備えたCCDカメラ等である。表示装置3は、CRT、液晶パネル等を有する表示装置である。また、ライト(光照射手段)9は、判定対象物及び良品に対してレーザ光を照射し、例えば正弦波縞パターンの位相をずらして撮像した複数の多値画像を取得することができる。
The camera (imaging means) 1 is a CCD camera or the like provided with a CCD image sensor. The
なお、一基のカメラ1と一基のライト9とを用いる構成に限定されるものではなく、カメラ1を複数用いても良い。複数のカメラ1、1、・・・により取得した複数の多値画像に基づいてステレオ計測法等を用いることにより、カメラ1から判定対象物までの距離を求め、距離画像を生成することができるからである。
The configuration is not limited to using one
また、カメラ1、表示装置3、ライト9等は、画像処理装置2と一体化されていても良いし、分離されていても良い。
The
外部制御機器6は、通信手段26を介して接続されている制御機器であり、例えばPLC(プログラマブルロジックコントローラ)等が相当する。ここで外部制御機器6とは、画像処理装置2による画像処理結果に応じて後処理を実行する機器全般を意味している。
The
図2は、本発明の実施の形態1に係る画像処理装置2の一構成例を示す機能ブロック図である。図2において、本実施の形態1に係る画像処理装置2は、カメラ1と、ライト9と、画像処理装置2の処理を実行する画像処理部7と、記憶手段23と、画像表示部8とから構成される。
FIG. 2 is a functional block diagram showing a configuration example of the
カメラ1は、例えばディジタルカメラであり、判定対象物として例えばフィルム表面を撮像し多値画像を取得して画像処理部7へ出力する。また、ライト9から光を照射することにより、判定対象物又は良品の表面に縞パターンを投影する。カメラ1は、投影された縞パターンを撮像する。画像処理装置2は、カメラ1により撮像された縞パターンを縞パターン画像として取得して、縞パターンの位相解析を実行することにより距離画像を生成し、記憶装置23に記憶する。なお、ライト9の光照射制御は、例えば主制御部21により行うことができる。
The
画像処理部7は、良品画像取得手段71と、良品距離画像生成手段72と、分布範囲算出手段73と、多値画像取得手段74と、対象距離画像生成手段75と、位置調整手段76と、判断手段77とを含む。また、画像処理部7は、主制御部21、メモリ22、外部I/F等を含んで構成され、良品画像取得手段71、良品距離画像生成手段72、分布範囲算出手段73、多値画像取得手段74、対象距離画像生成手段75、位置調整手段76及び判断手段77の処理動作を制御する。
The
記憶手段23は、画像メモリとして機能し、カメラ1で撮像した多値画像に関する多値画像データ、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像に基づいて生成した距離画像に関する距離画像データ、及び画像処理部7において位置合わせ、平均値算出等の各種処理を行った後の画像データを随時記憶する。画像データとしてではなく、画素ごとの画素値データとして記憶しても良い。
The
画像表示部8は、コンピュータ用のモニタ等の表示装置3で構成される。画像表示部8は、良品判定の対象となる判定対象物を撮像した多値画像、良品であるか否かの判定結果を、表示装置3の表示画面上に表示させる。すなわち多値画像表示手段81は、画像処理部7の指示に応じた多値画像を表示装置3の表示画面上に表示させ、判定結果表示手段82は、判定対象物が良品であるか否かの判定結果を表示装置3の表示画面上に表示させる。
The image display unit 8 includes a
次に、画像処理部7の各構成について説明する。
Next, each configuration of the
良品画像取得手段71は、カメラ1で撮像した、良品に関する複数の多値画像(第一の多値画像)を取得する。すなわち、良品に関する多値画像を複数取得しておくことにより、取得した多値画像の画素ごとに複数の特徴量を抽出することができる。取得した多値画像に関する多値画像データは、記憶手段23の良品画像データ記憶部231に記憶される。
The non-defective product
良品距離画像生成手段72は、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像に基づいて、縞パターンの位相解析を実行することにより、良品の表面からカメラ1までの距離に応じて画素ごとに濃淡値を特定した良品距離画像(第一の距離画像)を生成する。図3は、本発明の実施の形態1に係る距離画像を説明するための斜視図及び距離画像の例示図である。
The non-defective product distance image generation means 72 performs the phase analysis of the fringe pattern based on the fringe pattern image obtained by imaging the fringe pattern projected by irradiating light from the light 9 with the
図3(a)は、本発明の実施の形態1に係る距離画像を説明するための斜視図である。図3(a)の例では、半径の異なる2つの円柱を重ね合わせた形状を有するワーク301の上方にカメラ1を設置してある。ワーク301の最上面S1とカメラ1との距離をL1、ワーク301の中段面S2とカメラ1との距離をL2、ワーク301の載置面S3とカメラ1との距離をL3とし、ワークの高さHを(L3−L1)とする。
FIG. 3A is a perspective view for explaining a distance image according to
図3(b)は、本発明の実施の形態1に係る距離画像の例示図である。図3(b)に示すように、ワーク301の最上面S1が最も濃く(例えば黒色)、ワーク301の載置面S3が最も薄く(例えば白色)、ワーク301の中段面S2が中間色(例えば灰色)となっている。すなわち、ワーク301とカメラ1との間の距離L1、L2、L3に応じて、最上面S1、中段面S2、載置面S3の順に各画素の濃淡値が小さくなる距離画像を生成することができる。
FIG.3 (b) is an illustration figure of the distance image which concerns on
なお、距離画像の生成方法は、大別すると2つの方式に分類される。一方は、画像を取得するために特段の光照射手段を用いることなく撮像した画像を用いて生成するパッシブ方式であり、他方は、高さ方向の距離を計測するために光を能動的に照射して生成するアクティブ方式である。本実施の形態1では、ライト9から光を照射して生成するアクティブ方式を採用している。 The generation method of the distance image is roughly classified into two methods. One is a passive method that uses images taken without using special light irradiation means to acquire images, and the other is active irradiation with light to measure the distance in the height direction. It is an active method to generate. In the first embodiment, an active method is used in which light is emitted from the light 9 and generated.
二基のカメラを用いるステレオ計測法は、パッシブ方式の代表例である。ステレオ計測法は、二基のカメラを所定の位置に配置するだけで距離画像を生成することができる。 The stereo measurement method using two cameras is a representative example of the passive method. In the stereo measurement method, a distance image can be generated simply by placing two cameras at a predetermined position.
一方、アクティブ方式では、光切断法、パターン投影法が代表例として挙げられる。光切断法は、一基のカメラと一基のライト(光照射手段)とを用いて、ワークに対して線状のレーザ光を照射し、ワーク表面の形状に応じた光の像の歪具合に基づいて、ワークの三次元形状を復元する。 On the other hand, in the active method, a light cutting method and a pattern projection method are given as typical examples. The light cutting method uses a single camera and a single light (light irradiation means) to irradiate a workpiece with a linear laser beam and distort the light image according to the shape of the workpiece surface. To restore the three-dimensional shape of the workpiece.
パターン投影法は、ワークに投影された所定パターン(例えば縞パターン)の形状、位相等をずらして撮像した複数の多値画像を解析することで、ワークの三次元形状を復元する。本実施の形態1では、正弦波縞パターンの位相をずらして少なくとも三枚の多値画像を撮像し、画素ごとに正弦波の位相を求め、ワークの三次元形状を求める位相シフト法を用いている。 The pattern projection method restores the three-dimensional shape of the workpiece by analyzing a plurality of multi-valued images captured by shifting the shape, phase, and the like of a predetermined pattern (for example, a fringe pattern) projected onto the workpiece. In the first embodiment, the phase of the sine wave fringe pattern is shifted to capture at least three multivalued images, the phase of the sine wave is obtained for each pixel, and the phase shift method for obtaining the three-dimensional shape of the workpiece is used. Yes.
もちろん、これに限定されるものではなく、2つの規則的なパターンが合成される場合に生じる一種の空間周波数のうなり現象を利用して、ワークの三次元形状を復元するモアレポトグラフィ法を用いても良いし、光レーダ法、合焦点法、共焦点法、白色光干渉法等であっても良い。 Of course, the present invention is not limited to this, and a moire-potography method is used that restores the three-dimensional shape of the workpiece by using a kind of spatial frequency beat phenomenon that occurs when two regular patterns are synthesized. Alternatively, an optical radar method, a focusing method, a confocal method, a white light interference method, or the like may be used.
図2に戻って、分布範囲算出手段73は、生成した距離画像の画素ごとの濃淡値の、良品判定を行うための分布範囲を算出する。算出した分布範囲に関する情報は、分布範囲の境界を示す関数式、座標値、閾値等として、記憶手段23の分布範囲情報記憶部232に記憶される。
Returning to FIG. 2, the distribution range calculation means 73 calculates a distribution range for performing non-defective product determination of the gray value for each pixel of the generated distance image. Information on the calculated distribution range is stored in the distribution range
多値画像取得手段74は、カメラ1で撮像した、判定対象物に関する多値画像(第二の多値画像)を取得する。対象距離画像生成手段75は、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像に基づいて、縞パターンの位相解析を実行することにより、判定対象物の表面からカメラ1までの距離に応じて画素ごとの濃淡値を特定した対象距離画像(第二の距離画像)を生成する。良品距離画像の生成方法は、良品距離画像生成手段72と同様である。
The multi-value
位置調整手段76は、生成した判定対象物に関する対象距離画像と、良品に関する良品距離画像との位置合わせを行う。具体的には、良品に関する複数の良品距離画像の平均画像を算出しておき、平均画像との間で位置合わせを行う。距離画像間で位置合わせを行う手段は、周知の技術であれば特に限定されるものではない。例えば、パターンマッチングにより双方の距離画像の特徴点の位置を検出し、一致していると判断することが可能な上限値及び下限値を用いて位置合わせを行っても良い。また、正規化相関等を算出することにより双方の距離画像の一致度を算出し、一致度が所定値より大きくなるよう位置合わせを行っても良い。もちろん、双方の距離画像の輪郭線、面積、重心等を一致させるよう位置合わせを行っても良い。
The
なお、本実施の形態では、位置調整手段76は、生成した判定対象物に関する対象距離画像と、良品に関する良品距離画像との位置合わせを行っているが、特にこれに限定されるものではなく、例えば良品に関する複数の良品距離画像の平均画像を算出する際に複数の良品距離画像間で位置合わせを行っても良い。良品に関する複数の良品距離画像の平均画像を算出する際に複数の良品距離画像間で位置合わせを行う場合、位置調整手段76により位置合わせを行っても良いし、別途、良品距離画像間位置調整手段等を設けて行っても良い。
In the present embodiment, the
判断手段77は、生成した判定対象物に関する対象距離画像の画素ごとの濃淡値が、算出してある良品に関する良品距離画像の画素ごとの濃淡値の分布範囲に含まれているか否かを判断する。判定対象物に関する対象距離画像の画素ごとの濃淡値が分布範囲に含まれていると判断した場合には、判定対象物が良品であると判定し、含まれていないと判断した場合には、判定対象物が良品ではないと判定する。
The
図4は、本発明の実施の形態1に係る画像処理装置2の主制御部21の良品情報収集処理の手順を示すフローチャートである。画像処理装置2の主制御部21は、カメラ1で撮像した、良品に関する複数の多値画像(第一の多値画像)を取得する(ステップS401)。主制御部21は、取得した多値画像群の多値画像データを、記憶手段23の良品画像データ記憶部231に記憶する(ステップS402)。
FIG. 4 is a flowchart showing a non-defective information collection process procedure of the
主制御部21は、複数の良品に関して、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像を取得し(ステップS403)、縞パターンの位相解析を実行することにより、良品の表面からカメラ1までの距離に応じて画素ごとに濃淡値を特定した複数の良品距離画像を生成する(ステップS404)。主制御部21は、生成した良品距離画像群の距離画像データを、記憶手段23に記憶する(ステップS405)。
The
主制御部21は、生成した良品距離画像に基づいて、良品判定を行うための分布範囲を算出する(ステップS406)。主制御部21は、算出した分布範囲に関する情報を、分布範囲の境界を示す関数式、座標値、閾値等として、記憶手段23の分布範囲情報記憶部232に記憶する(ステップS407)。なお、良品判定を行うための分布範囲を算出する場合、複数の良品距離画像間で位置合わせを行っても良い。
The
図5は、本発明の実施の形態1に係る画像処理装置2の主制御部21の良品判定処理の手順を示すフローチャートである。画像処理装置2の主制御部21は、カメラ1で撮像した、判定対象物に関する多値画像(第二の多値画像)を取得する(ステップS501)。主制御部21は、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像を取得し(ステップS502)、縞パターンの位相解析を実行することにより、判定対象物の表面からカメラ1までの距離に応じて画素ごとに濃淡値を特定した対象距離画像を生成する(ステップS503)。
FIG. 5 is a flowchart showing a non-defective product determination process performed by the
主制御部21は、生成した判定対象物に関する対象距離画像と、良品に関する良品距離画像との位置合わせを行う(ステップS504)。具体的には、良品に関する複数の良品距離画像の平均距離画像を算出しておき、平均距離画像との間で位置合わせを行う。距離画像間で位置合わせを行う手段は、周知の技術であれば特に限定されるものではない。例えば、パターンマッチングにより双方の距離画像の特徴点の位置を検出し、一致していると判断することが可能な上限値及び下限値を用いて位置合わせを行っても良い。また、正規化相関等を算出することにより双方の距離画像の一致度を算出し、一致度が所定値より大きくなるよう位置合わせを行っても良い。もちろん、双方の距離画像の輪郭線、面積、重心等を一致させるよう位置合わせを行っても良い。
The
主制御部21は、生成した判定対象物に関する対象距離画像の画素ごとの濃淡値が、算出してある良品に関する良品距離画像の画素ごとの濃淡値の分布範囲に含まれているか否かを判断する(ステップS505)。主制御部21が、判定対象物に関する対象距離画像の画素ごとの濃淡値が分布範囲に含まれていると判断した場合(ステップS505:YES)、主制御部21は、判定対象物が良品であると判定し、判定対象物が良品である旨を示す情報を判定結果として、表示装置3の表示画面上に表示する(ステップS506)。
The
主制御部21が、判定対象物に関する対象距離画像の画素ごとの濃淡値が分布範囲に含まれていないと判断した場合(ステップS505:NO)、主制御部21は、判定対象物が良品ではないと判定し、判定対象物が良品ではない旨を示す情報を判定結果として、表示装置3の表示画面上に表示する(ステップS507)。
When the
以上のように本実施の形態1によれば、第二の距離画像に基づいて、良品であるか否かを判断するので、従来の方法では検出することが困難である判定対象物表面の凹凸を検出することができ、判定対象物を載置した位置のばらつき等により光の当たり方が変化した場合であっても、良品ではないと誤って判定される可能性を低減することが可能となる。 As described above, according to the first embodiment, since it is determined whether or not it is a non-defective product based on the second distance image, the unevenness on the surface of the determination target that is difficult to detect by the conventional method. It is possible to reduce the possibility of erroneously determining that the product is not a good product even if the way the light hits changes due to variations in the position where the determination target is placed, etc. Become.
(実施の形態2)
本発明の実施の形態2に係る画像処理装置の構成は、実施の形態1と同様であることから、同一の符号を付することにより詳細な説明を省略する。本実施の形態2は、複数の特徴量抽出手段を有し、取得した多値画像又は生成した距離画像それぞれの特徴量を抽出することにより、良品であるか否かを判断する点で、実施の形態1とは相違する。
(Embodiment 2)
Since the configuration of the image processing apparatus according to the second embodiment of the present invention is the same as that of the first embodiment, detailed description thereof is omitted by attaching the same reference numerals. The second embodiment has a plurality of feature amount extraction means, and it is implemented by determining whether or not it is a non-defective product by extracting the feature amounts of each of the acquired multi-value image or generated distance image. This is different from the first embodiment.
図6は、本発明の実施の形態2に係る画像処理装置2の一構成例を示す機能ブロック図である。図6において、本実施の形態2に係る画像処理装置2は、実施の形態1と同様、カメラ1と、ライト9と、画像処理装置2の処理を実行する画像処理部7と、記憶手段23と、画像表示部8とから構成される。
FIG. 6 is a functional block diagram showing a configuration example of the
カメラ1は、例えばディジタルカメラであり、判定対象物として例えばフィルム表面を撮像し多値画像を取得して画像処理部7へ出力する。また、ライト9から光を照射することにより、判定対象物又は良品の表面に縞パターンを投影する。カメラ1は、投影された縞パターンを撮像する。画像処理装置2は、カメラ1により撮像された縞パターンを縞パターン画像として取得して、縞パターンの位相解析を実行することにより距離画像を生成し、記憶装置23に記憶する。
The
画像処理部7は、良品画像取得手段71と、良品距離画像生成手段72と、選択受付手段78と、特徴量抽出手段79(選択特徴量抽出手段79c)と、分布範囲算出手段73と、多値画像取得手段74と、対象距離画像生成手段75と、位置調整手段76と、判断手段77とを含む。また、画像処理部7は、主制御部21、メモリ22、外部I/F等を含んで構成され、良品画像取得手段71、良品距離画像生成手段72、選択受付手段78、特徴量抽出手段79(選択特徴量抽出手段79c)、分布範囲算出手段73、多値画像取得手段74、対象距離画像生成手段75、位置調整手段76、及び判断手段77の処理動作を制御する。
The
記憶手段23は、画像メモリとして機能し、カメラ1により撮像した多値画像に関する多値画像データ、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像に基づいて生成した距離画像に関する距離画像データ、及び画像処理部7において位置合わせ、平均値算出等の各種処理を行った後の画像データを随時記憶する。画像データとしてではなく、画素ごとの画素値データとして記憶しても良い。
The
画像表示部8は、コンピュータ用のモニタ等の表示装置3で構成される。画像表示部8は、良品判定の対象となる判定対象物を撮像した多値画像、良品であるか否かの判定結果を、表示装置3の表示画面上に表示させる。すなわち多値画像表示手段81は、画像処理部7の指示に応じた多値画像を表示装置3の表示画面上に表示させ、判定結果表示手段82は、判定対象物が良品であるか否かの判定結果を表示装置3の表示画面上に表示させる。
The image display unit 8 includes a
次に、画像処理部7の各構成について説明する。
Next, each configuration of the
良品画像取得手段71は、カメラ1で撮像した、良品に関する複数の多値画像(第一の多値画像)を取得する。すなわち、良品に関する多値画像を複数取得しておくことにより、取得した多値画像の画素ごとに複数の特徴量を抽出することができる。取得した多値画像に関する多値画像データは、記憶手段23の良品画像データ記憶部231に記憶される。
The non-defective product
良品距離画像生成手段72は、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像に基づいて、縞パターンの位相解析を実行することにより、良品の表面からカメラ1までの距離に応じて画素ごとに濃淡値を特定した良品距離画像(第一の距離画像)を生成する。
The non-defective product distance image generation means 72 performs the phase analysis of the fringe pattern based on the fringe pattern image obtained by imaging the fringe pattern projected by irradiating light from the light 9 with the
選択受付手段78は、特徴量を抽出する複数の特徴量抽出手段のうち、少なくとも1つの特徴量抽出手段の選択を受け付ける。すなわち、本実施の形態2では、後述する第一の特徴量抽出手段79a又は第二の特徴量抽出手段79bの選択を受け付ける。なお、第一の特徴量抽出手段79a又は第二の特徴量抽出手段79bのいずれか1つの選択を受け付けても良いし、両方の選択を受け付けても良い。また、その他の特徴量抽出手段を含めた複数の特徴量抽出手段の中から選択を受け付けても良いことは言うまでもない。
The
特徴量抽出手段79は、第一の特徴量抽出手段79aと第二の特徴量抽出手段79bとを備えている。第一の特徴量抽出手段79aは、取得した多値画像又は生成した距離画像の画素ごとに、色成分に関する画素値を特徴量として抽出する。第二の特徴量抽出手段79bは、取得した多値画像又は生成した距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する。特徴量が色成分に関する画素値である場合、色の相違を検出することができ、特徴量が異なる二方向におけるエッジ強度である場合、輪郭線近傍の欠陥の存在を検出することができる。
The feature
分布範囲算出手段73は、選択を受け付けた特徴量抽出手段79により抽出した特徴量に基づいて、良品判定を行うための分布範囲を算出する。第一の特徴量抽出手段79aの選択を受け付けた場合には、分布範囲算出手段73aにより、第二の特徴量抽出手段79bの選択を受け付けた場合には、分布範囲算出手段73bにより、それぞれ分布範囲を算出する。算出した分布範囲に関する情報は、分布範囲の境界を示す関数式、座標値、閾値等として、記憶手段23の分布範囲情報記憶部232に記憶される。
The distribution
多値画像取得手段74は、カメラ1で撮像した、判定対象物に関する多値画像(第二の多値画像)を取得する。対象距離画像生成手段75は、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像に基づいて、縞パターンの位相解析を実行することにより、判定対象物の表面からカメラ1までの距離に応じて画素ごとに濃淡値を特定した対象距離画像(第二の距離画像)を生成する。対象距離画像の生成方法は、良品距離画像生成手段72と同様である。
The multi-value
位置調整手段76は、判定対象物に関する多値画像と、良品に関する多値画像との位置合わせ、及び生成した判定対象物に関する対象距離画像と、良品に関する良品距離画像との位置合わせを行う。具体的には、良品に関する複数の多値画像又は良品距離画像の平均画像を算出しておき、平均画像との間で位置合わせを行う。多値画像間又は距離画像間で位置合わせを行う手段は、周知の技術であれば特に限定されるものではない。例えば、パターンマッチングにより双方の多値画像又は距離画像の特徴点の位置を検出し、一致していると判断することが可能な上限値及び下限値を用いて位置合わせを行っても良い。また、正規化相関等を算出することにより双方の多値画像又は距離画像の一致度を算出し、一致度が所定値より大きくなるよう位置合わせを行っても良い。もちろん、双方の多値画像又は距離画像の輪郭線、面積、重心等を一致させるよう位置合わせを行っても良い。
The
選択特徴量抽出手段79cは、取得した判定対象物に関する多値画像又は生成した対象距離画像の画素ごとに、選択受付手段78にて選択を受け付けた特徴量抽出手段79と同じ特徴量抽出手段により特徴量を抽出する。例えば、第一の特徴量抽出手段79aの選択を受け付けた場合、判定対象物に関する多値画像又は対象距離画像の画素ごとに、色成分に関する画素値を特徴量として抽出し、第二の特徴量抽出手段79bの選択を受け付けた場合、判定対象物に関する多値画像又は対象距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する。
The selected feature
判断手段77は、取得した判定対象物に関する多値画像又は生成した対象距離画像の画素ごとに、選択を受け付けた特徴量抽出手段79により対応する特徴量を抽出し、抽出した特徴量が、算出してある良品に関する多値画像又は良品距離画像の、同じ特徴量の分布範囲に含まれているか否かを判断する。抽出した特徴量がすべて分布範囲に含まれていると判断した場合には、判定対象物が良品であると判定し、1つでも含まれていないと判断した場合には、判定対象物が良品ではないと判定する。
The
図7は、本発明の実施の形態2に係る画像処理装置2の主制御部21の良品情報収集処理の手順を示すフローチャートである。画像処理装置2の主制御部21は、カメラ1で撮像した、良品に関する複数の多値画像(第一の多値画像)を取得する(ステップS701)。主制御部21は、取得した多値画像群の多値画像データを、記憶手段23の良品画像データ記憶部231に記憶する(ステップS702)。
FIG. 7 is a flowchart showing a non-defective product information collection process performed by the
主制御部21は、複数の良品に関して、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像を取得し(ステップS703)、縞パターンの位相解析を実行することにより、良品の表面からカメラ1までの距離に応じて画素ごとの濃淡値を特定した複数の良品距離画像を生成する(ステップS704)。主制御部21は、生成した良品距離画像群の距離画像データを、記憶手段23に記憶する(ステップS705)。
The
主制御部21は、多値画像又は良品距離画像から特徴量を抽出する特徴量抽出手段79の選択を受け付ける(ステップS706)。選択対象となる特徴量抽出手段79は、取得した多値画像又は生成した良品距離画像の画素ごとに、色成分に関する画素値を特徴量として抽出する第一の特徴量抽出手段79a、及び取得した多値画像又は生成した良品距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する第二の特徴量抽出手段79bである。なお、選択を受け付ける特徴量抽出手段79は一つであっても良いし、複数であっても良い。
The
主制御部21は、選択を受け付けた特徴量抽出手段79により、多値画像又は良品距離画像の画素ごとに特徴量を抽出し(ステップS707)、良品判定を行うための分布範囲を抽出した特徴量ごとに算出する(ステップS708a、S708b)。選択を受け付けた特徴量抽出手段79が一つである場合には、例えばステップS708aのみにより、選択を受け付けた特徴量抽出手段79が二つである場合には、ステップS708a、S708bにより、それぞれに対応する分布範囲を、抽出した特徴量ごとに算出する。主制御部21は、算出した分布範囲に関する情報を、分布範囲の境界を示す関数式、座標値、閾値等として、記憶手段23の分布範囲情報記憶部232に記憶する(ステップS709)。
The
図8は、本発明の実施の形態2に係る画像処理装置2の主制御部21の良品判定処理の手順を示すフローチャートである。画像処理装置2の主制御部21は、カメラ1で撮像した、判定対象物に関する多値画像(第二の多値画像)を取得する(ステップS801)。主制御部21は、取得した判定対象物に関する多値画像と、良品に関する多値画像(第一の多値画像)との位置合わせを行う(ステップS802)。具体的には、良品に関する複数の多値画像の平均画像を算出しておき、平均画像との間で位置合わせを行う。
FIG. 8 is a flowchart showing a non-defective product determination process performed by the
主制御部21は、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像を取得し(ステップS803)、縞パターンの位相解析を実行することにより、判定対象物の表面からカメラ1までの距離に応じて画素ごとに濃淡値を特定した対象距離画像を生成する(ステップS804)。
The
主制御部21は、生成した判定対象物に関する対象距離画像と、良品に関する良品距離画像との位置合わせを行う(ステップS805)。具体的には、良品に関する複数の良品距離画像の平均距離画像を算出しておき、平均距離画像との間で位置合わせを行う。
The
主制御部21は、取得した判定対象物に関する多値画像又は生成した対象距離画像の画素ごとに、選択を受け付けた特徴量抽出手段79により特徴量を抽出する(ステップS806)。例えば、特徴量抽出手段79として色成分に関する画素値を抽出する第一の特徴量抽出手段79aの選択を受け付けた場合には、判定対象物に関する多値画像又は対象距離画像の画素ごとに、色成分に関する画素値を特徴量として抽出し、異なる二方向におけるエッジ強度を抽出する第二の特徴量抽出手段79bの選択を受け付けた場合には、判定対象物に関する多値画像又は対象距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する。
The
主制御部21は、多値画像又は対象距離画像の画素ごとに抽出した特徴量が、算出してある良品に関する多値画像又は良品距離画像の同じ特徴量に対応する分布範囲に含まれているか否かを判断する(ステップS807)。主制御部21が、抽出した多値画像又は対象距離画像の特徴量がすべて分布範囲に含まれていると判断した場合(ステップS807:YES)、主制御部21は、判定対象物が良品であると判定し、判定対象物が良品である旨を示す情報を判定結果として、表示装置3の表示画面上に表示する(ステップS808)。
Whether the
主制御部21が、抽出した多値画像又は対象距離画像の特徴量が1つでも分布範囲に含まれていないと判断した場合(ステップS807:NO)、主制御部21は、判定対象物が良品ではないと判定し、判定対象物が良品ではない旨を示す情報を判定結果として、表示装置3の表示画面上に表示する(ステップS809)。なお、特徴量抽出手段79の選択を複数受け付けた場合には、選択を受け付けた特徴量抽出手段79により抽出された特徴量を用いた良品判定の判定結果のうち、一つでも良品ではないと判定された場合には、判定対象物が良品ではないと判定すれば良い。
When the
また、上述したステップS806では、取得した判定対象物に関する多値画像又は生成した対象距離画像の画素ごとに、選択を受け付けた特徴量抽出手段79により特徴量を抽出しているが、例えば取得した判定対象物に関する多値画像と生成した対象距離画像との双方について、画素ごとに選択を受け付けた特徴量抽出手段79により特徴量を抽出することもできる。つまり、取得した判定対象物に関する多値画像の画素ごとに、選択を受け付けた特徴量抽出手段79により特徴量を抽出するとともに、生成した対象距離画像の画素ごとに、選択を受け付けた特徴量抽出手段79により特徴量を抽出することも可能である(特徴量は同じであっても異なっていても良い)。そして、取得した判定対象物に関する多値画像と生成した対象距離画像とのそれぞれについて、特徴量が分布範囲に含まれているか否かを判断し、取得した判定対象物に関する多値画像に着目したときの第一の良品判定及び生成した対象距離画像に着目したときの第二の良品判定を行う。例えば、いずれの良品判定においても特徴量が分布範囲に含まれていると判断した場合には、最終的に良品であると判定し、いずれか一方の良品判定において特徴量が分布範囲に含まれていないと判断した場合には、最終的に不良品であると判定する。このように、分布範囲算出手段73は、生成した距離画像だけでなく、多値画像の画素ごとの濃淡値の、良品判定を行うための分布範囲を算出しても良く、多値画像の画素ごとの濃淡値が、算出してある良品に関する良品多値画像の画素ごとの濃淡値の分布範囲に含まれているか否かを判断しても良い。なお、最終的に良品であると判定する条件のユーザ設定を受け付ける受付手段(マウス、モニタ等)を設けても良い。
Further, in step S806 described above, the feature amount is extracted by the feature
ここで、第一の特徴量抽出手段79aの選択を受け付けた場合、すなわち抽出される特徴量が、色成分に関する画素値である場合には、良品に関する多値画像又は良品距離画像の画素ごとの色成分に関する画素値の平均値を算出し、算出した色成分に関する画素値の平均値に基づいて、良品に関する多値画像又は良品距離画像の画素ごとに画素値の分布範囲を算出する。図9は、本発明の実施の形態2に係る画像処理装置2の第一の特徴量抽出手段79aの選択を受け付けた場合の分布範囲算出手段73aの機能ブロック図である。
Here, when the selection of the first feature
図9に示すように、分布範囲算出手段73aは、画素平均値算出手段731、画素値分布範囲算出手段732を含む。画像処理部7の主制御部21は、画素平均値算出手段731、画素値分布範囲算出手段732の処理動作を制御する。
As shown in FIG. 9, the distribution
画素平均値算出手段731は、取得した良品に関する多値画像又は生成した良品距離画像の画素ごとに、色成分に関する画素値、例えば濃淡値の平均値を算出する。画素値分布範囲算出手段732は、算出した色成分に関する画素値、例えば濃淡値の平均値に基づいて、良品に関する多値画像又は良品距離画像の画素ごとに色成分それぞれを座標軸とする多次元空間における画素値の分布範囲を算出する。
The pixel average
例えば、R成分、G成分、B成分それぞれを座標軸とした三次元座標系を想定する。単純にR成分、G成分、B成分に関する画素値の平均値をそれぞれ算出して、各成分の画素値の分布範囲を算出して良品判定を行う方法では、色成分に基づいて正しく判定することができない場合が生じうる。 For example, assume a three-dimensional coordinate system having R, G, and B components as coordinate axes. In the method of simply calculating the average value of the pixel values related to the R component, G component, and B component, and calculating the distribution range of the pixel values of each component to perform the non-defective product determination, the correct determination is made based on the color components. There are cases where it is not possible.
例えば良品が白に近似する色から濃いグレイまでを含む色、すなわちR成分の画素値が30〜240、G成分の画素値が30〜240、B成分の画素値が30〜240に分布する場合、赤に近い色の判定対象物、例えばR成分の画素値が200、G成分の画素値が50、B成分の画素値が100である判定対象物については、明らかに良品ではないにもかかわらず、各成分に関する画素値の分布範囲に含まれていると判断されるおそれがある。したがって、R成分、G成分、B成分について個別に判断するのではなく、R成分、G成分、B成分を座標軸とした仮想三次元空間における画素値の分布範囲に含まれているか否かに応じて良品であるか否かを判定する。 For example, when the non-defective product is distributed in a color including a color similar to white to dark gray, that is, the R component pixel values are 30 to 240, the G component pixel values are 30 to 240, and the B component pixel values are 30 to 240. A determination object having a color close to red, such as a determination object having an R component pixel value of 200, a G component pixel value of 50, and a B component pixel value of 100, is clearly not a good product. However, it may be determined that the pixel value is included in the distribution range of the pixel values for each component. Therefore, R component, G component, and B component are not individually determined, but depending on whether they are included in the pixel value distribution range in the virtual three-dimensional space with R component, G component, and B component as coordinate axes. Judge whether it is a good product.
図10は、本発明の実施の形態2に係る画像処理装置2の第一の特徴量抽出手段79aの選択を受け付けた場合の、画素値の分布範囲を示す例示図である。図10では、複数の良品に関する多値画像又は良品距離画像に基づいて、各色成分の画素値をR成分、G成分、B成分を互いに直交する座標軸とした三次元空間に点60、60、・・・としてプロットしている。そして、各色成分の画素値の平均値Avを算出して、算出した平均値Avを中心として、三次元空間に良品判定を行うための分布範囲61を算出している。
FIG. 10 is an exemplary diagram showing a distribution range of pixel values when selection by the first feature
例えば判定対象物の多値画像又は対象距離画像の画素ごとの各色成分に関する画素値をプロットした点62、63では、点63の方が平均値Avから離れている。しかし、点63は分布範囲61内であるのに対し、点62は分布範囲61の外側にプロットされている。したがって、点63としてプロットされた判定対象物は良品であると判定され、点62としてプロットされた判定対象物は良品ではないと判定される。
For example, at
分布範囲61の算出方法は特に限定されるものではない。例えば空間距離であるユークリッド距離を用いても良いし、マハラノビス距離を用いても良い。以下、マハラノビス距離を用いる場合について説明する。
The calculation method of the
まず、画素ごとの色成分に関する画素値のマハラノビス距離DM は、三次元座標軸、すなわち図10の例ではR軸方向、G軸方向及びB軸方向における画素ごとの画素ベクトルx(r、g、b)の平均値をそれぞれ平均値ベクトルμi (rバー、gバー、bバー)、分散共分散行列の逆行列をΣi -1 とした場合、ベクトル行列式を用いて(式1)のように算出することができる。 First, the Mahalanobis distance D M of the pixel value related to the color component for each pixel is the pixel vector x (r, g, The average value of b) is the average value vector μ i (r bar, g bar, b bar), and the inverse matrix of the covariance matrix is Σ i −1. In this case, it can be calculated as (Equation 1) using a vector determinant.
(式1)において、λj は固有値であり、ベクトルφj は固有値λj に対応した固有ベクトルである。すなわち(式1)で算出するマハラノビス距離DM は、点x(r、g、b)と分布範囲61の重心(平均値)との距離を、固有ベクトルφj 方向成分に分解し、固有ベクトルφj 方向成分を分散λj で正規化した距離と言える。なお、(式1)におけるiは分布の区別を示しており、上述した計算では分布範囲61が1つであることから特段の意味を有してはいない。
In (Expression 1), λ j is an eigenvalue, and the vector φ j is an eigenvector corresponding to the eigenvalue λ j . That is, the Mahalanobis distance D M calculated by (Equation 1) is obtained by decomposing the distance between the point x (r, g, b) and the center of gravity (average value) of the
仮想三次元空間であることから、(式1)において、n=3であり、固有値λj 、固有ベクトルφj はそれぞれ3個であり、固有ベクトルφj も三次元の行列となる。すなわち、分散共分散行列Σi は、(式2)のように表すことができる。 Since it is a virtual three-dimensional space, in (Expression 1), n = 3, there are three eigenvalues λ j and eigenvectors φ j , and the eigenvector φ j is also a three-dimensional matrix. That is, the variance-covariance matrix Σ i can be expressed as (Equation 2).
(式2)において、Vは分散を表しており、R成分、G成分、B成分それぞれの分散は(式3)により求めることができる。 In (Expression 2), V represents the dispersion, and the dispersion of each of the R component, the G component, and the B component can be obtained by (Expression 3).
また、(式2)において、Covは共分散を表しており、R成分とG成分との共分散、R成分とB成分との共分散、G成分とB成分との共分散は、それぞれ(式4)により求めることができる。 In (Expression 2), Cov represents covariance, and the covariance between the R component and the G component, the covariance between the R component and the B component, and the covariance between the G component and the B component are ( It can be obtained from equation 4).
(式4)の分母及び分子にnを乗算して変形することで、(式5)を得る。 (Formula 5) is obtained by multiplying n by the denominator and numerator of (Formula 4).
(式5)から、良品に関する多値画像又は良品距離画像に基づいて、画素ごとに、色成分に関する画素値の総和(平均値)、色成分に関する画素値の二乗和、及び二組ずつの色成分の画素値を乗算した値の総和を算出することにより、分散共分散行列Σi を求めることができることがわかる。色成分に関する画素値の総和(平均値)、色成分に関する画素値の二乗和、及び二組ずつの色成分の画素値を乗算した値の総和は、良品画像データ記憶部231に記憶しておき、良品であると判定する分布範囲61を算出する。
From (Equation 5), based on the multi-value image or the non-defective distance image regarding the non-defective product, the sum (average value) of the pixel values regarding the color component, the square sum of the pixel values regarding the color component, and two sets of colors for each pixel. It can be seen that the variance-covariance matrix Σ i can be obtained by calculating the sum of the values obtained by multiplying the component pixel values. The sum total of pixel values related to color components (average value), the sum of squares of pixel values related to color components, and the sum of values obtained by multiplying pixel values of two sets of color components are stored in the non-defective image
次に、判定対象物に関する多値画像又は距離画像の画素ごとの色成分に関する画素値について、(式1)で算出されるマハラノビス距離DM 、すなわち分布範囲61の重心から色成分に関する画素値までの距離を、分布範囲61の重心から色成分に関する画素値に向かう方向における分布範囲61の境界と重心との距離で正規化した値が、所定値より大きい場合には、画素値は分布範囲61に含まれていない。したがって、色が良品範囲を超えて変化したことを示しているので、良品ではないと判定することができる。
Next, regarding the pixel value related to the color component for each pixel of the multivalued image or the distance image related to the determination object, the Mahalanobis distance D M calculated by (Equation 1), that is, from the center of gravity of the
別の判定方法として、例えば仮想的な楕円球のような形状の範囲(以下、楕円球範囲)の重心から判定対象物に関する多値画像又は対象距離画像の画素ごとの画素値までの距離であるユークリッド距離を算出し、算出したユークリッド距離が、楕円球範囲の重心から判定対象物に関する多値画像又は対象距離画像の画素ごとの画素値に向かう方向における楕円球範囲の境界と重心との距離に基づく所定値より小さいか否かで相互相関分布範囲に含まれているか否かを判断する。ユークリッド距離を用いることにより、分布広がりの大小により感度が大きく変化することなく良品として判定される相互相関分布範囲から外れている度合いを算出することが可能となる。 As another determination method, for example, a distance from the center of gravity of a range of a shape like a virtual ellipsoid (hereinafter referred to as an elliptical sphere range) to a pixel value for each pixel of the multivalued image or target distance image related to the determination target. The Euclidean distance is calculated, and the calculated Euclidean distance is the distance between the center of gravity of the elliptic sphere range and the center of gravity in the direction from the centroid of the elliptic sphere range to the pixel value for each pixel of the multi-valued image or target distance image regarding the determination target. It is determined whether or not it is included in the cross-correlation distribution range based on whether or not it is smaller than a predetermined value. By using the Euclidean distance, it is possible to calculate the degree of deviating from the cross-correlation distribution range determined as non-defective without the sensitivity greatly changing due to the size of the distribution spread.
また、分布範囲61を楕円球範囲に限定する必要はなく、分布範囲61を内接する直方体形状の範囲を仮想的に算出しても良い。
Further, it is not necessary to limit the
図11は、本発明の実施の形態2に係る画像処理装置2の主制御部21の、第一の特徴量抽出手段79aの選択を受け付けた場合の、分布範囲算出処理の手順を示すフローチャートである。斯かる処理は、図7のステップS708aに相当する。
FIG. 11 is a flowchart showing the procedure of the distribution range calculation process when the
図11において、画像処理装置2の主制御部21は、良品に関する複数の多値画像又は良品距離画像に基づいて算出して、良品画像データ記憶部231に記憶してある、多値画像又は良品距離画像の画素ごとの色成分に関する画素値の総和(平均値)、色成分に関する画素値の二乗和、及び二組ずつの色成分に関する画素値を乗算した値の総和を読み出し(ステップS1101)、画素ごとの色成分に関する画素値の分散共分散行列Σiを作成する(ステップS1102)。すなわち、(式2)に示す分散共分散行列Σi を作成する。
In FIG. 11, the
主制御部21は、分散共分散行列Σi の固有値λj 、固有ベクトルφj を算出し(ステップS1103)、分布範囲に関する分布範囲情報として分布範囲情報記憶部232へ記憶する(ステップS1104)。
The
そして、図8のステップS807における分布範囲に含まれているか否かの判断は、取得した判定対象物に関する多値画像又は生成した対象距離画像の画素ごとの特徴量として、色成分に関する画素値を算出し、(式1)にて算出することができるマハラノビス距離DM が所定値よりも大きい場合には分布範囲61に含まれていないと判断し、良品ではないと判定することができる。
The determination as to whether or not the image is included in the distribution range in step S807 in FIG. 8 is performed by using a pixel value related to a color component as a feature value for each pixel of the acquired multi-value image or the generated target distance image. When the Mahalanobis distance D M that can be calculated and calculated by (Equation 1) is larger than a predetermined value, it is determined that it is not included in the
次に、第二の特徴量抽出手段79bの選択を受け付けた場合、すなわち抽出される特徴量が、異なる二方向におけるエッジ強度である場合には、良品に関する多値画像又は距離画像の画素ごとの異なる二方向におけるエッジ強度の平均値を算出し、算出した平均値を中心として、良品に関する多値画像又は良品距離画像の画素ごとに、異なる二方向にエッジ強度の相互相関分布範囲を算出する。図12は、本発明の実施の形態2に係る画像処理装置2の第二の特徴量抽出手段79bの選択を受け付けた場合の分布範囲算出手段73bの機能ブロック図である。
Next, when the selection of the second feature
図12に示すように、分布範囲算出手段73bは、方向別エッジ強度算出手段734、エッジ強度平均値算出手段735、及びエッジ強度分布範囲算出手段736を含む。画像処理部7の主制御部21は、方向別エッジ強度算出手段734、エッジ強度平均値算出手段735、及びエッジ強度分布範囲算出手段736の処理動作を制御する。
As shown in FIG. 12, the distribution
方向別エッジ強度算出手段734は、取得した良品に関する多値画像又は生成した良品距離画像の画素ごとの異なる二方向におけるエッジ強度を算出する。エッジ強度を算出する方向は特に限定されるものではないが、本実施の形態2では、異なる二方向である二次元画像の行方向及び列方向について、それぞれエッジ強度を算出する。具体的にはソーベルフィルタで用いる3×3のウインドウテンプレートを用いて、行方向及び列方向のエッジ強度を算出すれば良い。 The direction-specific edge strength calculation means 734 calculates edge strengths in two different directions for each pixel of the acquired multi-value image or the generated good product distance image. Although the direction in which the edge strength is calculated is not particularly limited, in the second embodiment, the edge strength is calculated for each of the row direction and the column direction of the two-dimensional image which are two different directions. Specifically, the edge intensity in the row direction and the column direction may be calculated using a 3 × 3 window template used in the Sobel filter.
エッジ強度平均値算出手段735は、算出した画素ごとの異なる二方向におけるエッジ強度に基づいて、良品に関する多値画像又は良品距離画像の画素ごとに、異なる二方向におけるエッジ強度の平均値を算出する。 The edge strength average value calculation means 735 calculates an average value of edge strengths in two different directions for each pixel of the multi-valued image or the non-defective distance image regarding the non-defective product based on the calculated edge strengths in two different directions for each pixel. .
エッジ強度分布範囲算出手段736は、算出した平均値を中心として、良品に関する多値画像又は良品距離画像の画素ごとに、異なる二方向におけるエッジ強度の相互相関分布範囲を算出する。 The edge intensity distribution range calculation means 736 calculates the cross-correlation distribution range of the edge intensity in two different directions for each pixel of the multi-valued image or the non-defective distance image regarding the non-defective product with the calculated average value as the center.
図13は、本発明の実施の形態2に係る画像処理装置2の所定の画素ごとの異なる二方向におけるエッジ強度の相互相関分布範囲の例示図である。図13では、二次元画像における行方向のエッジ強度exをX軸、列方向のエッジ強度eyをY軸としている。
FIG. 13 is an exemplary diagram of a cross-correlation distribution range of edge strengths in two different directions for each predetermined pixel of the
エッジ強度平均値算出手段735により、良品に関する多値画像又は良品距離画像の画素ごとに、二次元画像における行方向のエッジ強度exの平均値及び列方向のエッジ強度eyの平均値が算出されているので、良品に関する多値画像又は良品距離画像の画素ごとのエッジ強度の平均値ベクトルEAVE は、エッジ角度θを傾斜角とする直線上に位置する。そして、例えば取得した複数の良品に関する多値画像(第一の多値画像)又は生成した複数の良品に関する良品距離画像(第一の距離画像)の画素ごとのエッジ強度のばらつきが大きい場合であって、エッジ方向がほぼ揃っているときには、良品に関する多値画像又は良品距離画像のエッジ強度の平均値ベクトルEAVE を長軸及び短軸が交差する点を中心とした仮想的な楕円領域として相互相関分布領域91を算出する。したがって、判定対象物に関する多値画像又は対象距離画像の画素ごとのエッジ強度を算出して、同一の画素について、二次元画像における行方向のエッジ強度ex及び列方向のエッジ強度eyをプロットしたエッジ点92が、相互相関分布範囲91に含まれているか否かを判断することにより、判定対象物が良品であるか否かを判定することができる。
The edge strength average value calculation means 735 calculates the average value of the edge strength ex in the row direction and the average value of the edge strength ey in the column direction in the two-dimensional image for each pixel of the multi-value image or the good product distance image related to the good product. Therefore, the average value vector E AVE of the edge strength for each pixel of the multi-value image or the non-defective distance image regarding the non-defective product is located on a straight line having the edge angle θ as the inclination angle. For example, this is a case where there is a large variation in edge intensity for each pixel of the acquired multi-value image (first multi-value image) related to a plurality of non-defective products or the generated good-quality distance image (first distance image) related to a plurality of non-defective products. When the edge directions are almost aligned, the average value vector E AVE of the edge strengths of the multi-valued image or the good-quality distance image for the non-defective product is mutually displayed as a virtual elliptical area centered on the point where the major axis and the minor axis intersect. A
相互相関分布範囲91の算出方法は特に限定されるものではない。例えば空間距離であるユークリッド距離を用いても良いし、マハラノビス距離を用いても良い。以下、マハラノビス距離を用いる場合について説明する。
The calculation method of the
この場合、図6に示す判断手段77の距離算出手段70は、分布範囲の重心、例えば図13の例では仮想的な楕円領域である相互相関分布範囲91の重心(中心点)から判定対象物に関する多値画像又は対象距離画像の画素ごとのエッジ点(特徴量)までの距離を、重心からエッジ点(特徴量)に向かう方向における相互相関分布範囲91の境界と重心との距離で正規化したマハラノビス距離を算出する。
In this case, the distance calculation means 70 of the determination means 77 shown in FIG. 6 determines the object to be determined from the centroid of the distribution range, for example, the centroid (center point) of the
画素ごとのエッジ強度のマハラノビス距離DM
は、異なる二方向、すなわち図13の例ではX軸方向及びY軸方向における画素ごとのエッジ強度ベクトルx(ex 、ey )の平均値をそれぞれ平均値ベクトルμi (ex バー、ey バー)、分散共分散行列の逆行列をΣi -1 とした場合、ベクトル行列式を用いて(式6)のように算出することができる。
Mahalanobis distance D M of edge intensity for each pixel
Is two different directions, i.e. the average value of each average value vector μ i (e x bar in the example of FIG. 13 X-axis direction and Y each pixel in the direction of the edge intensity vector x (e x, e y) , e y bar), and the inverse of the variance-covariance matrix Σ i -1 In this case, it can be calculated as (Expression 6) using a vector determinant.
(式6)において、λj は固有値であり、ベクトルφj は固有値λj に対応した固有ベクトルである。すなわち(式6)で算出するマハラノビス距離DM は、点x(ex 、ey )と相互相関分布範囲91の重心(平均値)との距離を、固有ベクトルφj 方向成分に分解し、固有ベクトルφj 方向成分を分散λj で正規化した距離と言える。なお、(式6)におけるiは分布の区別を示しており、上述した計算では相互相関分布範囲91が1つであることから特段の意味を有してはいない。
In (Expression 6), λ j is an eigenvalue, and vector φ j is an eigenvector corresponding to the eigenvalue λ j . That Mahalanobis distance D M is calculated by (Equation 6), the point x (e x, e y) the distance between the center of gravity of the cross-correlation distribution range 91 (average value), and degraded eigenvector phi j direction component, eigenvector It can be said that the φ j direction component is normalized by the variance λ j . Note that i in (Equation 6) indicates the distinction of the distribution, and since the cross
図10とは異なり、異なる二方向の二次元空間であることから、(式6)において、n=2であり、固有値λj 、固有ベクトルφj はそれぞれ2個であり、固有ベクトルφj も二次元の行列となる。すなわち、分散共分散行列Σi は、(式7)のように表すことができる。 Unlike FIG. 10, since it is a two-dimensional space in two different directions, in (Equation 6), n = 2, eigenvalue λ j and eigenvector φ j are two, and eigenvector φ j is also two-dimensional. It becomes the matrix of. That is, the variance-covariance matrix Σ i can be expressed as (Equation 7).
(式7)において、Vは分散を表しており、X軸方向のエッジ強度、Y軸方向のエッジ強度それぞれの分散は(式8)により求めることができる。 In (Expression 7), V represents the dispersion, and the dispersion of the edge strength in the X-axis direction and the edge strength in the Y-axis direction can be obtained by (Expression 8).
また、(式7)において、Covは共分散を表しており、X軸方向のエッジ強度とY軸方向のエッジ強度との共分散は、(式9)により求めることができる。 In (Expression 7), Cov represents covariance, and the covariance between the edge strength in the X-axis direction and the edge strength in the Y-axis direction can be obtained by (Expression 9).
(式9)から、良品に関する多値画像又は良品距離画像に基づいて、画素ごとに、異なる二方向におけるエッジ強度の総和(平均値)、異なる二方向におけるエッジ強度の二乗和、及び異なる二方向におけるエッジ強度を互いに乗算した値の総和を算出することにより、分散共分散行列Σi を求めることができることがわかる。異なる二方向におけるエッジ強度の総和(平均値)、異なる二方向におけるエッジ強度の二乗和、及び異なる二方向におけるエッジ強度を互いに乗算した値の総和は、良品画像データ記憶部231に記憶しておき、良品であると判定する相互相関分布範囲91を算出する。
From (Equation 9), based on a multi-value image or a non-defective distance image regarding a non-defective product, for each pixel, the sum (average value) of edge strengths in two different directions, the square sum of edge strengths in two different directions, and two different directions It can be seen that the variance-covariance matrix Σ i can be obtained by calculating the total sum of the values obtained by multiplying the edge intensities at. The sum (average value) of edge strengths in two different directions, the square sum of edge strengths in two different directions, and the sum of values obtained by multiplying the edge strengths in two different directions are stored in the non-defective image
次に、判定対象物に関する多値画像又は対象距離画像の画素ごとの異なる二方向(X軸方向及びY軸方向)におけるエッジ強度について、(式6)で算出されるマハラノビス距離DM 、すなわち相互相関分布範囲91の重心からエッジ点(特徴量)までの距離を、相互相関分布範囲91の重心からエッジ点(特徴量)に向かう方向における相互相関分布範囲91の境界と重心との距離で正規化した値が、所定値より大きい場合には、エッジ点は相互相関分布範囲91に含まれていない。したがって、エッジ角度が良品範囲を超えて変動したことを示しているので輪郭線の形状が大きく変化したと判断することができ、良品ではないと判定することができる。
Next, regarding the edge strength in two different directions (X-axis direction and Y-axis direction) for each pixel of the multivalued image or the target distance image related to the determination target, the Mahalanobis distance D M calculated by (Expression 6), that is, mutual The distance from the centroid of the
図14は、本発明の実施の形態2に係る画像処理装置2の主制御部21の、第二の特徴量抽出手段79bの選択を受け付けた場合の、分布範囲算出処理の手順を示すフローチャートである。斯かる処理は、図7のステップS708bに相当する。
FIG. 14 is a flowchart showing the procedure of the distribution range calculation process when the
図14において、画像処理装置2の主制御部21は、良品に関する多値画像又は良品距離画像に基づいて算出して記憶してある、多値画像又は良品距離画像の画素ごとの異なる二方向におけるエッジ強度の総和(平均値)、異なる二方向におけるエッジ強度の二乗和、及び異なる二方向におけるエッジ強度を互いに乗算した値の総和を読み出し(ステップS1401)、画素ごとの異なる二方向におけるエッジ強度の分散共分散行列Σiを作成する(ステップS1402)。
In FIG. 14, the
主制御部21は、分散共分散行列Σi の固有値λj 、固有ベクトルφj を算出し(ステップS1403)、相互相関分布範囲に関する分布範囲情報として分布範囲情報記憶部232へ記憶する(ステップS1404)。
The
そして、図8のステップS807における分布範囲に含まれているか否かの判断は、取得した判定対象物に関する多値画像又は生成した対象距離画像の画素ごとの特徴量として、異なる二方向におけるエッジ強度を算出し、(式6)にて算出することができるマハラノビス距離DM が所定値よりも大きい場合には相互相関分布範囲91に含まれていないと判断し、良品ではないと判定することができる。
Then, the determination as to whether or not the image is included in the distribution range in step S807 in FIG. 8 is based on the edge strength in two different directions as the feature value for each pixel of the obtained multi-value image or the generated target distance image. When the Mahalanobis distance D M that can be calculated by (Equation 6) is larger than a predetermined value, it is determined that it is not included in the
図15は、本発明の実施の形態2に係る画像処理装置2の主制御部21の、第二の特徴量抽出手段79bの選択を受け付けた場合の、良品判定処理の手順を示すフローチャートである。画像処理装置2の主制御部21は、カメラ1で撮像した、判定対象物に関する多値画像を取得する(ステップS1501)。主制御部21は、取得した判定対象物に関する多値画像と、良品に関する多値画像との位置合わせを行う(ステップS1502)。具体的には、良品に関する複数の多値画像の平均画像を算出しておき、平均画像との間で位置合わせを行う。多値画像間で位置合わせを行う手段は、周知の技術であれば特に限定されるものではない。
FIG. 15 is a flowchart showing a non-defective product determination process when the
主制御部21は、ライト9から光を照射することにより投影された縞パターンをカメラ1で撮像した縞パターン画像を取得し(ステップS1503)、縞パターンの位相解析を実行することにより、判定対象物の表面からカメラ1までの距離に応じて画素ごとに濃淡値を特定した対象距離画像を生成する(ステップS1504)。
The
主制御部21は、生成した判定対象物に関する対象距離画像と、良品に関する良品距離画像との位置合わせを行う(ステップS1505)。具体的には、良品に関する複数の良品距離画像の平均距離画像を算出しておき、平均距離画像との間で位置合わせを行う。
The
主制御部21は、取得した判定対象物に関する多値画像又は生成した対象距離画像の画素ごとに、異なる二方向におけるエッジ強度を算出する(ステップS1506)。エッジ強度を算出する方向は特に限定されるものではないが、上述した方法と同様、異なる二方向である二次元画像における列方向及び行方向について、それぞれエッジ強度を算出する。
The
主制御部21は、算出した判定対象物に関する多値画像又は対象距離画像の画素ごとのエッジ強度が、算出した良品に関する多値画像又は良品距離画像の画素ごとのエッジ強度の相互相関分布範囲に含まれているか否かを判断する(ステップS1507)。主制御部21が、算出したエッジ強度がすべて相互相関分布範囲に含まれていると判断した場合(ステップS1507:YES)、主制御部21は、判定対象物が良品であると判定し、判定対象物が良品である旨を示す情報を判定結果として、表示装置3の表示画面上に表示する(ステップS1508)。主制御部21が、算出したエッジ強度が1つでも相互相関分布範囲に含まれていないと判断した場合(ステップS1507:NO)、主制御部21は、判定対象物が良品ではないと判定し、判定対象物が良品ではない旨を示す情報を判定結果として、表示装置3の表示画面上に表示する(ステップS1509)。
The
なお、上述した実施の形態2では、相互相関分布範囲が、仮想的な楕円領域である場合について説明しているが、特に楕円領域に限定されるものではない。 In the second embodiment described above, the case where the cross-correlation distribution range is a virtual elliptical area has been described. However, the present invention is not particularly limited to the elliptical area.
以上のように本実施の形態2によれば、取得した判定対象物に関する多値画像(第二の多値画像)又は対象距離画像(第二の距離画像)の性状に応じて、良品であるか否かを判定する基礎となる特徴量を抽出する手段を少なくとも一つ選択することができるので、適切な特徴量抽出手段79を選択することにより、欠陥が輪郭線近傍に生じている場合であっても欠陥の存在を検出することができ、色の相違についても高い感度で検出することが可能となる。例えば選択を受け付けた特徴量抽出手段79が、色成分に関する画素値を特徴量として抽出する場合、色の相違による濃淡のばらつきの存在を検出することができ、選択を受け付けた特徴量抽出手段79が異なる二方向におけるエッジ強度を特徴量として抽出する場合、輪郭線近傍の欠陥の存在を検出することができる。
As described above, according to the second embodiment, the multi-value image (second multi-value image) or the target distance image (second distance image) related to the acquired determination target object is a non-defective product. Since at least one means for extracting a feature quantity as a basis for determining whether or not a defect has occurred in the vicinity of the contour line by selecting an appropriate feature quantity extraction means 79 can be selected. Even if it exists, the presence of a defect can be detected, and a difference in color can be detected with high sensitivity. For example, when the feature
なお、本発明は上記実施例に限定されるものではなく、本発明の趣旨の範囲内であれば多種の変更、改良等が可能である。異なる二方向を二次元画像の行方向及び列方向としているが、特に限定されるものではなく、互いに直交する二方向であれば足り、色成分についてもR成分、G成分、B成分だけではなく、シアン成分、マゼンタ成分、イエロー成分であっても良い。 The present invention is not limited to the above-described embodiments, and various changes and improvements can be made within the scope of the present invention. Two different directions are used as the row direction and the column direction of the two-dimensional image. However, the two directions are not particularly limited, and the two directions are orthogonal to each other. The color components are not limited to the R component, the G component, and the B component. , A cyan component, a magenta component, and a yellow component.
1 カメラ(撮像手段)
2 画像処理装置
3 表示装置(画像表示手段)
4 可搬型記録媒体
5 コンピュータプログラム
6 外部制御機器
7 画像処理部
8 画像表示部
9 ライト(光照射手段)
21 主制御部
22 メモリ
23 記憶手段
24 入力手段
25 出力手段
26 通信手段
27 補助記憶手段
28 内部バス
1 Camera (imaging means)
2
4
21
Claims (11)
撮像手段で撮像された、良品に関する複数の第一の多値画像を取得する良品画像取得手段と、
取得した第一の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第一の距離画像を生成する良品距離画像生成手段と、
前記第一の多値画像の画素ごとに、色成分に関する画素値を特徴量として抽出する第一の特徴量抽出手段と、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する第二の特徴量抽出手段とを少なくとも含む複数の特徴量抽出手段のうち、少なくとも1つの特徴量抽出手段の選択を受け付ける選択受付手段と、
取得した第一の多値画像及び生成した前記第一の距離画像の少なくとも一方の画素ごとの濃淡値の、良品判定を行うための分布範囲を算出する分布範囲算出手段と、
判定対象物に関する第二の多値画像を取得する多値画像取得手段と、
取得した第二の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第二の距離画像を生成する対象距離画像生成手段と、
取得した第二の多値画像及び生成した前記第二の距離画像の少なくとも一方の画素ごとの濃淡値が対応する前記分布範囲に含まれているか否かを判断する判断手段と
を備え、
前記分布範囲算出手段は、前記第一の特徴量抽出手段の選択を受け付けた場合には第一の分布範囲を、前記第二の特徴量抽出手段の選択を受け付けた場合には第二の分布範囲を、それぞれ選択を受け付けた特徴量抽出手段により抽出した特徴量に基づいて算出するようにしてあり、
前記判断手段は、前記第二の多値画像又は前記第二の距離画像の画素ごとに、選択を受け付けた特徴量抽出手段により特徴量を抽出し、抽出された特徴量に対応する前記第一の分布範囲又は前記第二の分布範囲に含まれているか否かを判断するようにしてあることを特徴とする画像処理装置。 In an image processing apparatus that performs a non-defective product determination by comparing a multi-value image obtained by imaging a determination target with a multi-value image group related to a good product,
A non-defective image acquisition means for acquiring a plurality of first multi-value images related to the non-defective product imaged by the imaging means;
A non-defective distance image generating means for generating a first distance image obtained by converting the distance from the imaging means into a gray value for each pixel of the acquired first multi-valued image;
For each pixel of the first multi-valued image, for each pixel of the first multi-valued image or the first distance image, a first feature amount extracting means for extracting a pixel value related to a color component as a feature amount. Selection accepting means for accepting selection of at least one feature quantity extracting means among a plurality of feature quantity extracting means including at least second feature quantity extracting means for extracting edge strengths in two different directions as feature quantities;
The gray value of each at least one pixel of the acquired first multivalued image and generated the first distance image, the distribution range calculation means for calculating a distribution range for performing good judgment,
Multi-valued image acquisition means for acquiring a second multi-valued image related to the determination target;
Target distance image generation means for generating a second distance image obtained by converting the distance from the imaging means into a gray value for each pixel of the acquired second multi-valued image;
Determination means for determining whether or not the gray value for each pixel of at least one pixel of the acquired second multi-valued image and the generated second distance image is included in the corresponding distribution range, and
The distribution range calculation means displays the first distribution range when the selection of the first feature quantity extraction means is accepted, and the second distribution when the selection of the second feature quantity extraction means is accepted. The range is calculated based on the feature amount extracted by the feature amount extraction unit that has accepted the selection,
The determination unit extracts a feature amount by a feature amount extraction unit that has received a selection for each pixel of the second multi-valued image or the second distance image, and the first corresponding to the extracted feature amount An image processing apparatus is characterized in that it is determined whether or not the distribution range is included in the distribution range or the second distribution range.
前記第一の多値画像の画素ごとに、色成分に関する画素値の平均値を算出する画素平均値算出手段と、
算出した色成分に関する画素値の平均値に基づいて、前記第一の多値画像の画素ごとに、色成分それぞれを座標軸とする多次元空間における画素値の前記第一の分布範囲を算出する画素値分布範囲算出手段と
を有し、
前記判断手段は、前記第二の多値画像の画素ごとに、色成分に関する画素値を算出して、色成分それぞれを座標軸とする多次元空間における画素値の前記第一の分布範囲に含まれているか否かを判断するようにしてあり、
前記第二の特徴量抽出手段の選択を受け付けた場合、前記分布範囲算出手段は、
前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を算出する方向別エッジ強度算出手段と、
算出したエッジ強度に基づいて、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度の平均値を算出するエッジ強度平均値算出手段と、
算出した平均値を中心として、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度の相互相関分布範囲を前記第二の分布範囲として算出するエッジ強度分布範囲算出手段と
を有し、
前記判断手段は、前記第二の多値画像又は前記第二の距離画像の画素ごとに、異なる二方向におけるエッジ強度を算出して、算出したエッジ強度が対応する前記相互相関分布範囲に含まれているか否かを判断するようにしてあることを特徴とする請求項1乃至3のいずれか一項に記載の画像処理装置。 When receiving the selection of the first feature amount extraction unit, the distribution range calculation unit,
Pixel average value calculating means for calculating an average value of pixel values related to color components for each pixel of the first multi-valued image;
Pixels for calculating the first distribution range of pixel values in a multidimensional space with each color component as a coordinate axis for each pixel of the first multi-valued image, based on the average value of the pixel values relating to the calculated color component A value distribution range calculation means, and
The determination means calculates a pixel value related to a color component for each pixel of the second multi-value image, and is included in the first distribution range of pixel values in a multi-dimensional space with each color component as a coordinate axis. To judge whether or not
When receiving the selection of the second feature amount extraction means, the distribution range calculation means,
Direction-specific edge strength calculating means for calculating edge strength in two different directions for each pixel of the first multi-valued image or the first distance image;
Based on the calculated edge strength, for each pixel of the first multi-value image or the first distance image, an edge strength average value calculating means for calculating an average value of edge strengths in two different directions;
Edge strength for calculating the cross-correlation distribution range of edge strengths in two different directions as the second distribution range for each pixel of the first multi-valued image or the first distance image with the calculated average value as the center A distribution range calculation means, and
The determining means calculates edge strengths in two different directions for each pixel of the second multi-valued image or the second distance image, and the calculated edge strengths are included in the corresponding cross-correlation distribution range. 4. The image processing apparatus according to claim 1, wherein the image processing apparatus is configured to determine whether or not the image processing apparatus is present.
撮像手段で撮像された、良品に関する複数の第一の多値画像を取得するステップと、
取得した第一の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第一の距離画像を生成するステップと、
前記第一の多値画像の画素ごとに、色成分に関する画素値を特徴量として抽出する第一の特徴量抽出手段と、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する第二の特徴量抽出手段とを少なくとも含む複数の特徴量抽出手段のうち、少なくとも1つの特徴量抽出手段の選択を受け付けるステップと、
取得した第一の多値画像及び生成した前記第一の距離画像の少なくとも一方の画素ごとの濃淡値の、良品判定を行うための分布範囲を算出するステップと、
判定対象物に関する第二の多値画像を取得するステップと、
取得した第二の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第二の距離画像を生成するステップと、
取得した第二の多値画像及び生成した前記第二の距離画像の少なくとも一方の画素ごとの濃淡値が対応する前記分布範囲に含まれているか否かを判断するステップと
を含み、
前記第一の特徴量抽出手段の選択を受け付けた場合には第一の分布範囲を、前記第二の特徴量抽出手段の選択を受け付けた場合には第二の分布範囲を、それぞれ選択を受け付けた特徴量抽出手段により抽出した特徴量に基づいて算出し、
前記第二の多値画像又は前記第二の距離画像の画素ごとに、選択を受け付けた特徴量抽出手段により特徴量を抽出し、抽出された特徴量に対応する前記第一の分布範囲又は前記第二の分布範囲に含まれているか否かを判断することを特徴とする画像処理方法。 In an image processing method that can be executed by an image processing apparatus that performs a non-defective product determination by comparing a multi-value image obtained by imaging a determination target with a multi-value image group related to a non-defective product,
Acquiring a plurality of first multi-valued images related to non-defective products imaged by the imaging means;
Generating a first distance image obtained by converting the distance from the imaging unit into a gray value for each pixel of the acquired first multi-valued image;
For each pixel of the first multi-valued image, for each pixel of the first multi-valued image or the first distance image, a first feature amount extracting means for extracting a pixel value related to a color component as a feature amount. Receiving a selection of at least one feature quantity extraction means among a plurality of feature quantity extraction means including at least a second feature quantity extraction means for extracting edge strengths in two different directions as feature quantities;
Calculating the gray value of each at least one pixel of the acquired first multivalued image and generated the first distance image, the distribution range for performing good judgment,
Obtaining a second multi-valued image relating to the determination object;
Generating a second distance image obtained by converting the distance from the imaging unit into a gray value for each pixel of the acquired second multi-valued image;
Determining whether a gray value for each pixel of at least one pixel of the acquired second multi-valued image and the generated second distance image is included in the corresponding distribution range, and
When the selection of the first feature quantity extraction unit is accepted, the selection is accepted, and when the selection of the second feature quantity extraction unit is accepted, the second distribution range is accepted. Calculated based on the feature amount extracted by the feature amount extraction means,
For each pixel of the second multi-value image or the second distance image, a feature amount is extracted by a feature amount extraction unit that has received a selection, and the first distribution range corresponding to the extracted feature amount or the An image processing method comprising determining whether or not the second distribution range is included.
前記画像処理装置を、
撮像手段で撮像された、良品に関する複数の第一の多値画像を取得する良品画像取得手段、
取得した第一の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第一の距離画像を生成する良品距離画像生成手段、
前記第一の多値画像の画素ごとに、色成分に関する画素値を特徴量として抽出する第一の特徴量抽出手段と、前記第一の多値画像又は前記第一の距離画像の画素ごとに、異なる二方向におけるエッジ強度を特徴量として抽出する第二の特徴量抽出手段とを少なくとも含む複数の特徴量抽出手段のうち、少なくとも1つの特徴量抽出手段の選択を受け付ける選択受付手段、
取得した第一の多値画像及び生成した前記第一の距離画像の少なくとも一方の画素ごとの濃淡値の、良品判定を行うための分布範囲を算出する分布範囲算出手段、
判定対象物に関する第二の多値画像を取得する多値画像取得手段、
取得した第二の多値画像の画素ごとに、前記撮像手段からの距離を濃淡値に変換した第二の距離画像を生成する対象距離画像生成手段、及び
取得した第二の多値画像及び生成した前記第二の距離画像の少なくとも一方の画素ごとの濃淡値が対応する前記分布範囲に含まれているか否かを判断する判断手段
として機能させ、
前記分布範囲算出手段を、前記第一の特徴量抽出手段の選択を受け付けた場合には第一の分布範囲を、前記第二の特徴量抽出手段の選択を受け付けた場合には第二の分布範囲を、それぞれ選択を受け付けた特徴量抽出手段により抽出した特徴量に基づいて算出する手段として機能させ、
前記判断手段を、前記第二の多値画像又は前記第二の距離画像の画素ごとに、選択を受け付けた特徴量抽出手段により特徴量を抽出し、抽出された特徴量に対応する前記第一の分布範囲又は前記第二の分布範囲に含まれているか否かを判断する手段として機能させることを特徴とするコンピュータプログラム。 In a computer program that can be executed by an image processing apparatus that performs a non-defective product determination by comparing a multi-value image obtained by imaging a determination target with a multi-value image group related to a non-defective product,
The image processing apparatus;
Non-defective image acquisition means for acquiring a plurality of first multi-value images related to non-defective products imaged by the imaging means,
Non-defective distance image generation means for generating a first distance image obtained by converting the distance from the imaging means into a gray value for each pixel of the acquired first multi-valued image,
For each pixel of the first multi-valued image, for each pixel of the first multi-valued image or the first distance image, a first feature amount extracting means for extracting a pixel value related to a color component as a feature amount. Selection accepting means for accepting selection of at least one feature quantity extracting means among a plurality of feature quantity extracting means including at least second feature quantity extracting means for extracting edge strengths in two different directions as feature quantities;
Acquired first gray value of at least every one of the pixels of the multivalued image and the generated the first distance image, the distribution range calculation means for calculating a distribution range for performing good judgment,
Multi-valued image obtaining means for obtaining a second multi-valued image relating to the determination object;
Target distance image generation means for generating a second distance image obtained by converting the distance from the imaging means into a gray value for each pixel of the acquired second multi-valued image; and
Function as a determination means for determining whether or not the gray value for each pixel of at least one pixel of the acquired second multi-valued image and the generated second distance image is included in the corresponding distribution range;
When the selection of the first feature amount extraction unit is accepted, the distribution range calculation unit displays the first distribution range, and when the selection of the second feature amount extraction unit is accepted, the second distribution The range functions as a means for calculating based on the feature amount extracted by the feature amount extraction means that accepts the selection,
The determination unit extracts the feature amount by the feature amount extraction unit that has received selection for each pixel of the second multi-valued image or the second distance image, and the first feature corresponding to the extracted feature amount A computer program that functions as means for determining whether or not the distribution range is included in the second distribution range.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010161125A JP5684505B2 (en) | 2010-07-16 | 2010-07-16 | Image processing apparatus, image processing method, and computer program |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010161125A JP5684505B2 (en) | 2010-07-16 | 2010-07-16 | Image processing apparatus, image processing method, and computer program |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2012021919A JP2012021919A (en) | 2012-02-02 |
JP5684505B2 true JP5684505B2 (en) | 2015-03-11 |
Family
ID=45776311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010161125A Expired - Fee Related JP5684505B2 (en) | 2010-07-16 | 2010-07-16 | Image processing apparatus, image processing method, and computer program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5684505B2 (en) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4054188B2 (en) * | 2001-11-30 | 2008-02-27 | 富士通株式会社 | Semiconductor device |
JP2014106164A (en) * | 2012-11-29 | 2014-06-09 | Mitsubishi Motors Corp | Method and apparatus for inspecting fitting state of connector |
JP6405124B2 (en) * | 2014-06-09 | 2018-10-17 | 株式会社キーエンス | Inspection device, inspection method, and program |
JP6314798B2 (en) * | 2014-11-12 | 2018-04-25 | Jfeスチール株式会社 | Surface defect detection method and surface defect detection apparatus |
US11158039B2 (en) | 2015-06-26 | 2021-10-26 | Cognex Corporation | Using 3D vision for automated industrial inspection |
EP3109826B1 (en) * | 2015-06-26 | 2019-02-13 | Cognex Corporation | Using 3d vision for automated industrial inspection |
JP2018112471A (en) * | 2017-01-11 | 2018-07-19 | リコーエレメックス株式会社 | Inspection system and inspection method |
JP2018112470A (en) * | 2017-01-11 | 2018-07-19 | リコーエレメックス株式会社 | Inspection system and inspection method |
JP6530536B2 (en) * | 2018-06-20 | 2019-06-12 | 株式会社キーエンス | Three-dimensional shape measuring apparatus, measurement data processing unit, measurement data processing method, and computer program |
JP6585793B2 (en) * | 2018-09-18 | 2019-10-02 | 株式会社キーエンス | Inspection device, inspection method, and program |
JP7416071B2 (en) * | 2019-08-23 | 2024-01-17 | コニカミノルタ株式会社 | Judgment device and judgment program |
CN116297524B (en) * | 2023-04-07 | 2023-11-14 | 苏州矽行半导体技术有限公司 | Multi-mode detection method for wafer with image |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2006282725A (en) * | 2005-03-31 | 2006-10-19 | Kri Inc | Silicon-containing new optically active compound |
JP4928862B2 (en) * | 2006-08-04 | 2012-05-09 | 株式会社日立ハイテクノロジーズ | Defect inspection method and apparatus |
JP2008139074A (en) * | 2006-11-30 | 2008-06-19 | Rozefu Technol:Kk | Method for detecting defect in image |
JP2010146357A (en) * | 2008-12-19 | 2010-07-01 | Panasonic Electric Works Co Ltd | Method and apparatus for three-dimensional image processing |
-
2010
- 2010-07-16 JP JP2010161125A patent/JP5684505B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2012021919A (en) | 2012-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP5684505B2 (en) | Image processing apparatus, image processing method, and computer program | |
US10775160B2 (en) | System and method for efficient surface measurement using a laser displacement sensor | |
WO2022052313A1 (en) | Calibration method for 3d structured light system, and electronic device and storage medium | |
EP1788345B1 (en) | Image processing device and image processing method performing 3d measurement | |
JP4885584B2 (en) | Rangefinder calibration method and apparatus | |
JP6506731B2 (en) | System and method for scoring clutter used for 3D point cloud matching in vision system | |
US8581162B2 (en) | Weighting surface fit points based on focus peak uncertainty | |
US10497111B2 (en) | Information processing apparatus and method of selecting viewpoint for measuring object, and measuring system | |
JP6899189B2 (en) | Systems and methods for efficiently scoring probes in images with a vision system | |
AU2019203928B2 (en) | Face location detection | |
US10508995B2 (en) | Three-dimensional image inspection device, three-dimensional image inspection method, three-dimensional image inspection program, and computer-readable recording medium and recording equipment | |
JP5956814B2 (en) | Appearance inspection apparatus, appearance inspection method, and computer program | |
JP6256475B2 (en) | Image processing apparatus, image processing method, and image processing program | |
US10771776B2 (en) | Apparatus and method for generating a camera model for an imaging system | |
KR20150128300A (en) | method of making three dimension model and defect analysis using camera and laser scanning | |
JP5469532B2 (en) | Image processing apparatus, image processing method, and computer program | |
WO2013061976A1 (en) | Shape inspection method and device | |
KR102073468B1 (en) | System and method for scoring color candidate poses against a color image in a vision system | |
JP7064404B2 (en) | Optical displacement meter | |
KR101510206B1 (en) | Urban Change Detection Method Using the Aerial Hyper Spectral images for Digital Map modify Drawing | |
WO2021024499A1 (en) | Reinforcing bar determination device and reinforcing bar determination method | |
JP2010133744A (en) | Defect detection method, and visual inspection device using the same | |
JP5486403B2 (en) | Image processing apparatus, image processing method, and computer program | |
JP4001162B2 (en) | Image processing method, image processing program and storage medium therefor, and image processing apparatus | |
JP6867766B2 (en) | Information processing device and its control method, program |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130409 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131126 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131129 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140114 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140722 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140917 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20150113 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150115 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5684505 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |