A selected physical attribute of a material can be analyzed to determine the uniformity of the material, which in turn can provide useful information regarding the appearance and functionality of the material in a particular product application. Methods for analyzing and determining uniformity have relied on pictorial standards and the judgment of human experts, but such qualitative methods lack precision and cannot be utilized in real-time as a product is manufactured.
Optical methods have been used to measure physical properties of materials in real-time. However, rapidly evaluating the overall uniformity of a material based on these measurements has proven to be difficult, as some non-uniformities are present at small size scales, while others are apparent only at larger size scales.
A. A method for characterizing the uniformity of a material, including selecting a set of size scales at which to measure uniformity within an area of interest in an image of the material; suppressing features in the image smaller than a selected size scale of interest within the set of size scales; dividing the image into patches equal to the size scale of interest; and calculating a uniformity value within each patch.
B. The method of embodiment A, wherein suppressing the features comprises processing the image with a low-pass filter.
C. The method of embodiments A-B, wherein the low pass filter comprises a box filter with a cutoff frequency equal to a predetermined fraction of the size scale of interest.
D. The method of embodiments A-B, wherein the low-pass filter comprises a two-dimensional Gaussian kernel.
E. The method of any of embodiments A-D, wherein the uniformity value is calculated by determining at least one of a standard deviation, an inter-quartile range (IQR), a median absolute deviation, or an information entropy of a selected characteristic of the patch.
F. The method of embodiment E, wherein the selected characteristic of the patch comprises an intensity of light transmitted through the patch or reflected off a surface of the material comprising the patch.
G. The method of embodiment E, wherein the uniformity value is calculated by determining the inter-quartile range (IQR) of an intensity of the light transmitted through the patch.
H. The method of any of embodiments A-G, further comprising scaling the area of interest to a predetermined size prior to removing the features.
I. The method of any of embodiments A-H, further comprising calibrating the area of interest prior to removing the features.
J. The method of any of embodiments A-G, further comprising aggregating the uniformity values of the patches to determine a uniformity value for the area of interest.
K. The method of any of embodiments A-J, further comprising aggregating the uniformity values of a selected array of patches within the area of interest to provide an uniformity value for the area of interest.
L. The method of any of embodiments A-L, wherein the material is selected from wovens, non-wovens, paper, coatings, polymeric films and combinations thereof.
M. The method of embodiment L, wherein the material is a non-woven.
N. The method of any of embodiments A-M, wherein the image is obtained by transmitting light through or reflecting light off a surface of the material comprising the area of interest.
O. The method of embodiment N, wherein the image is obtained by transmitting light through the material to an optical receiving device.
P. A method for characterizing the uniformity of a material, including obtaining an image of an area of interest of the material by transmitting light through the material to an optical receiving device; selecting a graduated set of size scales at which to measure uniformity within the area of interest; convolving a low-pass filter with the image to suppress features in the image smaller than a selected size scale of interest within the graduated set of size scales; dividing the image into patches equal to the size scale of interest, wherein the patches each comprise an array of pixels; and determining a standard deviation of the light intensity in the pixels in the array to calculate a uniformity value within each patch.
Q. The method of embodiment P, wherein the low pass filter comprises a box filter with a width equal to a predetermined fraction of the pixels within the array.
R. The method of embodiment Q, wherein the low-pass filter replaces a selected pixel in the array with a weighted average of the light intensities of the pixels surrounding the selected pixel, and wherein the weighted average is determined by a two-dimensional Gaussian kernel.
S. The method of any of embodiments P-R, further comprising determining an ideal pixel size for analyzing a selected non-uniformity, and scaling the area of interest to the ideal pixel size prior to removing the features.
T. The method of any of embodiments P-S, further comprising calibrating the area of interest prior to removing the features.
U. The method of any of embodiments P-T, further comprising aggregating the uniformity values of the patches to determine a uniformity value for the area of interest.
V. The method of any of embodiments P-U, further comprising aggregating the uniformity values of a selected array of patches within the area of interest to provide an uniformity value for the area of interest.
W. The method of any of embodiments P-V, wherein the material is selected from wovens, non-wovens, paper, coatings, polymeric films and combinations thereof.
X. The method of embodiment W, wherein the material is a non-woven.
Y. An apparatus, comprising: at least one light source illuminating a web of a material; a camera that captures light transmitted through or reflected from an area of interest on the material to generate an image of the area of interest; and a processor which, in response to an input of a set of size scales at which to measure uniformity within the area of interest: convolves a low-pass filter with the image to suppress features in the image smaller than a selected size scale of interest within the set of size scales; divides the image into patches equal to the size scale of interest, wherein the patches each comprise an array of pixels; and calculates a uniformity value within each patch.
Z. The apparatus of embodiment Y, wherein the processor calculates the uniformity value by determining at least one of a standard deviation, an inter-quartile range (IQR), a median absolute deviation, or an information entropy of a light intensity in the pixels in the array.
AA. The apparatus of any of embodiments Y-Z, wherein the processor calculates the uniformity value by determining the inter-quartile range (IQR).
BB. The apparatus of any of embodiments Y-AA, wherein the low pass filter comprises a box filter with a width equal to a predetermined fraction of the pixels within the array.
CC. The apparatus of any of embodiments Y-BB, wherein the low-pass filter replaces a selected pixel in the array with a weighted average of the light intensities of the pixels surrounding the selected pixel, and wherein the weighted average is determined by a two-dimensional Gaussian kernel.
DD. The apparatus of any of embodiments Y-CC, wherein the processor further determines an ideal pixel size for analyzing a selected non-uniformity in the material, and scales the area of interest to the ideal pixel size prior to removing the features.
EE. The apparatus of any of embodiments Y-DD, wherein the processor calibrates the area of interest prior to removing the features.
FF. The apparatus of any of embodiments Y-EE, wherein the processor aggregates the uniformity values of the patches to determine a uniformity value for the area of interest.
GG. The apparatus of any of embodiments Y-FF, wherein the processor aggregates the uniformity values of a selected array of patches within the area of interest to provide a uniformity value for the area of interest.
HH. The apparatus of any of embodiments Y-GG, wherein the material is selected from non-wovens and polymeric films.
II. The apparatus of embodiment HH, wherein the material is a non-woven.
JJ. The apparatus of any of embodiments Y-HH, wherein the camera captures light transmitted through the area of interest.
KK. The apparatus of embodiment JJ, wherein only scattered light is captured by the camera to form the image.
LL. The apparatus of embodiment KK, wherein a dark stripe is placed across the light source, and the camera is aimed directly at the dark stripe.
MM. A method for inspecting web material in real time and computing a level of uniformity of an area of interest in a surface of the web material as the web material is manufactured, comprising: obtaining an image of the area of interest by transmitting light through the material to an optical receiving device; selecting a set of size scales at which to measure uniformity within the area of interest; convolving a low-pass filter with the image to suppress features in the image smaller than a selected size scale of interest within the set of size scales; dividing the image into patches equal to the size scale of interest, wherein the patches each comprise an array of pixels; and calculating a uniformity value within each patch.
NN. The method of embodiment MM, wherein the uniformity value is calculated by determining the inter-quartile range (IQR) of the light intensity in the pixels in the array.
OO. The method of any of embodiments MM-NN, wherein the low pass filter comprises a box filter with a width equal to a predetermined fraction of the pixels within the array.
PP. The method of any of embodiments MM-OO, wherein the low-pass filter replaces a selected pixel in the array with a weighted average of the light intensities of the pixels surrounding the selected pixel, and wherein the weighted average is determined by a two-dimensional Gaussian kernel.
QQ. The method of any of embodiments MM-PP, further comprising determining an ideal pixel size for analyzing a selected non-uniformity, and scaling the area of interest to the ideal pixel size prior to removing the features.
RR. The method of any of embodiments MM-QQ, further comprising calibrating the area of interest prior to removing the features.
SS. The method of any of embodiments MM-RR, further comprising aggregating the uniformity values of the patches to determine a uniformity value for the area of interest.
TT. The method of any of embodiments MM-SS, further comprising aggregating the uniformity values of a selected array of patches within the area of interest to provide an uniformity value for the area of interest.
UU. The method of any of embodiments MM-TT, wherein the material is a non-woven.
VV. An online computerized inspection system for inspecting web material in real time, the system comprising: at least one light source illuminating a web of a material; a camera that captures light transmitted through or reflected from an area of interest on the material to generate an image of the area of interest; and a computer executing software to characterize the uniformity of the material in the area of interest, wherein the computer comprises a processor which, in response to an input of a set of size scales at which to measure uniformity within the area of interest: convolves a low-pass filter with the image to suppress features in the image smaller than a selected size scale of interest within the set of size scales; divides the image into patches equal to the size scale of interest, wherein the patches each comprise an array of pixels; and calculates a uniformity value within each patch.
WW. The system of embodiment VV, further comprising a memory to store a web inspection model, wherein the computer executes software to compare the uniformity in the area of interest to the model and compute a severity of a non-uniformity defect in the material.
XX. The system of any of embodiments VV-WW, further comprising a user interface to output the severity of the defect to a user.
YY. The system of any of embodiments VV-XX, wherein the material is a non-woven.
ZZ. A non-transitory computer readable medium comprising software instructions to cause a computer processor to: receive, with an online computerized inspection system, an image of one or more areas of interest on a surface of a web material during the manufacture thereof; select a set of size scales at which to measure uniformity within the area of interest; convolve a low-pass filter with the image to suppress features in the image smaller than a selected size scale of interest within the set of size scales; divide the image into patches equal to the size scale of interest, wherein the patches each comprise an array of pixels; and calculate a uniformity value within each patch.
AAA. The computer readable medium of embodiment ZZ, wherein the uniformity value is calculated by determining the inter-quartile range (IQR) of the light intensity in the pixels in the array.
BBB. The computer readable medium of any of embodiments ZZ-AAA, wherein the low pass filter comprises a box filter with a width equal to a predetermined fraction of the pixels within the array.
CCC. The computer readable medium of any of embodiments ZZ-BBB, wherein the low-pass filter replaces a selected pixel in the array with a weighted average of the light intensities of the pixels surrounding the selected pixel, and wherein the weighted average is determined by a two-dimensional Gaussian kernel.
DDD. The computer readable medium of any of embodiments ZZ-CCC, wherein the processor further determines an ideal pixel size for analyzing a selected non-uniformity, and scaling the area of interest to the ideal pixel size prior to removing the features.
EEE. The computer readable medium of any of embodiments ZZ-DDD, wherein the processor calibrates the area of interest prior to removing the features.
FFF. The computer readable medium of any of embodiments ZZ-EEE, wherein the processor aggregates the uniformity values of the patches to determine a uniformity value for the area of interest.
GGG. The computer readable medium of any of embodiments ZZ-FFF, wherein the processor aggregates the uniformity values of a selected array of patches within the area of interest to provide a uniformity value for the area of interest.
HHH. The computer readable medium of any of embodiments ZZ-GGG, wherein the material is a non-woven.
The terms “about” or “approximately” with reference to a numerical value, property, or characteristic, means+/−five percent of the numerical value, property, characteristic, but also expressly includes any narrow range within the +/−five percent of the numerical value or property or characteristic as well as the exact numerical value. For example, a temperature of “about” 100° C. refers to a temperature from 95° C. to 105° C., inclusive, but also expressly includes any narrower range of temperature or even a single temperature within that range, including, for example, a temperature of exactly 100° C.
The term “substantially” with reference to a property or characteristic means that the property or characteristic is exhibited to within 98% of that property or characteristic, but also expressly includes any narrow range within the two percent of the property or characteristic, as well as the exact value of the property or characteristic. For example, a substrate that is “substantially” transparent refers to a substrate that transmits 98-100%, inclusive, of the incident light.
The details of one or more embodiments of the invention are set forth in the accompanying drawings and the description below. Other features, objects, and advantages of the invention will be apparent from the description and drawings, and from the claims.
Depending on the needs of a particular application, the uniformity of a material at some size scales might be deemed as more important than at others. For example, in applications where the material is to be converted into small pieces, a non-uniformity at a size scale much larger than these pieces may not have any cosmetic impact, since it will not be visible within the extent of a single small piece. On the other hand, larger-scale non-uniformities may cause differences in functional properties between samples.
In general, the present disclosure is directed to an inspection technique to characterize the uniformity of an area of interest in a sample of a material. Images of the area of interest are captured using an optical inspection system and processed using image processing algorithms to determine the uniformity (or conversely, non-uniformity) over a set of size scales (for example, from 1 mm to 10 cm in increments of 1 mm) selected to evaluate a particular appearance or performance metric of the material.
At each size scale of interest within the set of size scales, the captured image is processed to remove image features significantly smaller than the size scale of interest. The image is then divided into areas equal to the size scale of interest, referred to herein as patches, and the uniformity of each patch is evaluated using a metric selected to provide a quantitative and repeatable measurement of the overall uniformity for that patch size.
The measurements obtained for each patch may then be utilized for further analysis of the area of interest or of the material. For example, the uniformity measurements for all or a selected group of patches within the area of interest may be aggregated to calculate a single uniformity value for the area of interest. The calculated uniformity values may be utilized in a wide variety of ways. In one non-limiting example, the uniformity values may be monitored in real-time during manufacturing operations to assess the appearance of the material or the functionality of the product of which the material is a part.
In
For applications where the web material is more dense or sheet-like, a top-light configuration may be used (not shown in
The dark field light sources 118B and 118C also can be conveniently provided as fiber light lines, although in some embodiments a laser source or other source may be employed. The dark field sources 118B and 118C conveniently illuminate a strip across the entire width of the web of material 112, oriented along the cross web direction. However, in some embodiments, they are mounted at an angle to the direction normal to the plane of the web of material 112
In some embodiments, lenses can be used to focus the light emerging from the direct light source 118A and the two dark field sources 118B, 118C. For example, when fiber light lines are used as the sources, cylindrical lenses oriented parallel to the fiber light lines may be used. The cylindrical lenses 132 and 134 preferably focus the light from far dark field sources 118B and 118C onto the underside of the web 112 on a line directly under the camera 120. The cylindrical lens 136 that focuses light from direct light source 118A can have the same focal length as cylindrical lenses 132 and 134, but the light from the direct source 11A is directed onto a diffuser 138. In one embodiment, the diffuser 138 is a diffuser film. A dark line on the diffuser film, or another solid object, for example a taut cable or narrow piece of metal 140 mounted just above (or on) the diffuser 138, can provide the dark stripe as shown in
Referring to
In some embodiments, after the image is acquired by the camera in step 202, prior to application of further image processing algorithms, the image may optionally be calibrated in step 204, and the image intensities mapped according to the calibration. The images obtained by the camera in
For example, for a dark-stripe illumination inspection system such as described in
One of the intensity profiles of
IB(L),L=1,2, . . . ,N, wherein N is the number of levels,
and some other intensity profile, from a system to be calibrated, is given by
I(L),L=1,2, . . . ,N, wherein N is the number of levels.
Then we assume that there is a linear mapping from I(L)→Ig(L), parameterized by m and b, such that
IB(L)=m I(L)+b.
The linear mapping is acquired via a linear least-squares regression, as illustrated in the example of
After the optional image calibration step 204, an image processing technique may be applied to measure the overall uniformity of an area of interest in a material at various size scales s1, s2, . . . , sn. Computation of uniformity across a range of sizes may capture some non-uniformities present at small size scales, while others are apparent only on larger size scales. Depending on the needs of a particular application, the uniformity at some size scales might be deemed as more important than at others. For example, in applications where a material is to be converted into small patches, a non-uniformity at a size scale much larger than these patches may not have any cosmetic or functional impact, since it will not be visible within the extend of and single small patch. On the other hand, larger-scale non-uniformities may cause differences in functional properties between samples. These are but two examples of the types of application-specific considerations that can be taken into account when choosing the range of size scales over which to estimate uniformity.
Referring again to
In step 208 of
In some embodiments, the low pass filter is a “box filter,” which consists of a two-dimensional kernel consisting of identical values. When convolved with an image, the box filter replaces each pixel in the size scale under consideration with the average of all neighboring pixel values. In other embodiments, a two-dimensional Gaussian kernel low-pass filter may be used, which can have more favorable characteristics in the frequency domain. When convolved with an image, the two-dimensional Gaussian kernel replaces each pixel with a weighted average of the intensities of the surrounding pixels, where the weights are given by the Gaussian kernel.
Regardless of the type of low-pass filter selected for a particular application, the algorithm suppresses high-frequency components of the image, which consist of image features much smaller than the size scale of interest. The low-pass filter allows measurement of only non-uniformities that are roughly near the size scale of interest, which removes the effect in a given patch caused by non-uniformities at much lower size scales. The smaller non-uniformities are captured at smaller size scales in the multiscale processing algorithms.
The application of a low-pass filter can be thought of in terms of how an observer visually perceives non-uniformities when physically looking at a sample. That is, when the observer stands close to the sample, very fine details of the surface are apparent, but not the overall uniformity on a large scale. On the other hand, when the observer stands far away from the sample, the overall uniformity and variations dominate the image, but the observer can no longer detect the fine level of detail that may exist at smaller size scales.
For example, in each iteration of the low-pass filtering algorithm described above, the low-pass filter can be selected to have a cutoff frequency equal to a predefined fraction of the current size scale at which to measure uniformity. In one specific example, if the size scale under consideration corresponds to 100 pixels, a box filter with a width of 20 pixels might be selected to suppress non-uniformities that are outside the size scale of interest.
Once the image is filtered to remove or reduce the impact of image features that are non-essential to the uniformity analysis at the selected size scale, in step 210 of
Referring to step 212 in
where μ(I) is the mean intensity in the patch, and N is the total number of pixels in it.
Other possible uniformity metrics include inter-quartile range (IQR), median absolute deviation (MAD), and the information entropy, among others. In some embodiments, the IQR, which is defined as the difference between the 75th and 25th percentile intensity values in the sample area, is more robust to outliers.
This uniformity analysis is computed for each patch using the metrics each time a new image is acquired by the camera 20 and the acquisition module 22 (
In step 215, the image processing steps 208, 210, 212, 214 are repeated for each size scale s1, s2, . . . , sn.
In some embodiments, as shown in step 216 of
The optical inspection system shown in
The analysis computer 24 (
The database 25 may be implemented in any of a number of different forms including a data storage file or one or more database management systems (DBMS) executing on one or more database servers. The database management systems may be, for example, a relational (RDBMS), hierarchical (HDBMS), multidimensional (MDBMS), object oriented (ODBMS or OODBMS) or object relational (ORDBMS) database management system. As one example, the database 25 is implemented as a relational database available under the trade designation SQL Server from Microsoft Corporation, Redmond, Wash.
Once the process has ended, the analysis computer 24 may transmit the data collected in the database 25 to a conversion control system 60 via a network 50. For example, the analysis computer 24 may communicate the roll information as well as the uniformity information and respective sub-images for each uniformity measurement to the conversion control system 60 for subsequent, offline, detailed analysis. For example, the uniformity information may be communicated by way of database synchronization between the database 25 and the conversion control system 60.
In some embodiments, the conversion control system 60 may determine those products of products for which each anomaly may cause a defect, rather than the analysis computer 24. Once data for the finished web roll has been collected in the database 25, the data may be communicated to converting sites and/or used to mark anomalies on the web roll, either directly on the surface of the web with a removable or washable mark, or on a cover sheet that may be applied to the web before or during marking of anomalies on the web.
The components of the analysis computer 24 may be implemented, at least in part, as software instructions executed by one or more processors of the analysis computer 24, including one or more hardware microprocessors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field programmable gate arrays (FPGAs), or any other equivalent integrated or discrete logic circuitry, as well as any combinations of such components. The software instructions may be stored within in a non-transitory computer readable medium, such as random access memory (RAM), read only memory (ROM), programmable read only memory (PROM), erasable programmable read only memory (EPROM), electronically erasable programmable read only memory (EEPROM), flash memory, a hard disk, a CD-ROM, a floppy disk, a cassette, magnetic media, optical media, or other computer-readable storage media.
Although shown for purposes of example as positioned within a manufacturing plant, the analysis computer 24 may be located external to the manufacturing plant, e.g., at a central location or at a converting site. For example, the analysis computer 24 may operate within the conversion control system 60. In another example, the described components execute on a single computing platform and may be integrated into the same software system.
The subject matter of the present disclosure will now be described with reference to the following non-limiting examples.
A dark field imaging system generally as depicted in
The nonwoven web material was formed from polymeric strands using a conventional air-laid process, having a mean thickness of approximately 1.25 cm, but also having significant variations in thickness, including relatively thin and thick portions in close proximity.
The light source was placed approximately 10 cm below the surface of the web, and the camera was located approximately 102 cm above it. Dark-stripe illumination was used, as described in
While the production line was running and the web was being conveyed, crossweb image lines at a fixed spacing along the downweb direction were acquired by the imaging system. These lines, when stacked together, formed a 2-dimensional image. These images were then processed according to the multi-scale uniformity algorithms described above.
In this case, the non-uniformity was measured at 5 different size scales of interest. These were 10 mm, 25 mm, 50 mm, 75 mm, and 100 mm. For each image, a single metric of non-uniformity was computed at each size scale. These were displayed as a time-series plot, with non-uniformity vs. time displayed at each of the 5 size scales, resulting in 5 different curves on the plot. An example of a plot obtained in this manner is depicted in
Various embodiments of the invention have been described. These and other embodiments are within the scope of the following claims.
This application is a national stage filing under 35 U.S.C. 371 of PCT/US2014/060898, filed Oct. 16, 2014, which claims the benefit of U.S. Application No. 61/898,304 filed Oct. 31, 2013, the disclosure of which is incorporated by reference in its/their entirety herein.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2014/060898 | 10/16/2014 | WO | 00 |
Publishing Document | Publishing Date | Country | Kind |
---|---|---|---|
WO2015/065726 | 5/7/2015 | WO | A |
Number | Name | Date | Kind |
---|---|---|---|
5774574 | Hoki | Jun 1998 | A |
6239554 | Tessadro | May 2001 | B1 |
6301373 | Bernie | Oct 2001 | B1 |
20050004956 | Pourdeyhimi | Jan 2005 | A1 |
20090051907 | Li | Feb 2009 | A1 |
20090310864 | Takagi et al. | Dec 2009 | A1 |
20110157584 | Lin | Jun 2011 | A1 |
20120327214 | McEntyre | Dec 2012 | A1 |
20130033620 | Polidor | Feb 2013 | A1 |
20130083324 | Wilhelm | Apr 2013 | A1 |
Number | Date | Country |
---|---|---|
11-94767 | Apr 1999 | JP |
2001-266122 | Sep 2001 | JP |
2002-228594 | Aug 2002 | JP |
200951878 | Dec 2009 | TW |
2012115819 | Aug 2012 | WO |
WO 2013-015699 | Jan 2013 | WO |
WO 2013-051786 | Apr 2013 | WO |
Entry |
---|
International Search Report for PCT International Application No. PCT/US2014/060898 dated Jan. 27, 2015, 4 pages. |
Kenneth Ivan Laws A: “Textured image segmentation”, Uscript Report 940, Universtiy of Southern California 1980, US, val. 84, Jan. 1, 1980 (Jan. 1, 1980), pp. 1-195, XP008184787, Retrieved from the Internet: URL:http://www.dlic.mil/docs/citations/ADA083283. |
Bresee R R et al: “Characterizing Nonwoven Web Structure Using Image Analysis Techniques”, Tappi Journal, Technical Association of the Pulp & Paper Industry. Atlanta, US, vol. 80, No. 7, Jul. 1, 1997 (Jul. 1, 1997), pp. 133-138, XP000751858, ISSN: 0734-1415. |
Nickolay B et al: “Automatisierte Visuelle Inspektion Auf Basis Einer Hierarchischen Bildaufloesung”, Zvvf Zeitschrift Fur Wirtschaftliche Fertigung Undautomatisierung, Carl Hanser Verlag. Munchen, DE, vol. 87, No. 1, Jan. 1, 1992 (Jan. 1, 1992), pp. 34-37, XP000271912, ISSN: 0947-0085. |
Nickolay B: “Uberwacht lernendes Bildauswertungssystem zur Erkennung von Oberflachenfehlern”, Aug. 1, 1990 (Aug. 1, 1990), Carl Hanser Verlag Munchen, Berlin, ISBN: 3-446-16176-7, val. 1, pp. 32-39-116-137. |
Nickolay B et al: “Texturanalyse mittels Co-OccuranceMatrizen zur autmatisierten Oberflaechenpruefung”, Jan. 1, 1988 (Jan. 1, 1988), Conference Info: ZWF83 (1988), Carl Hanser Verlag, Miinchen 1988, DE, pp. 90 94, XP008184789. |
Number | Date | Country | |
---|---|---|---|
20160282279 A1 | Sep 2016 | US |
Number | Date | Country | |
---|---|---|---|
61898304 | Oct 2013 | US |