CN104849278A - Method for measuring cut tobacco filament rate on basis of similarity coefficient - Google Patents
Method for measuring cut tobacco filament rate on basis of similarity coefficient Download PDFInfo
- Publication number
- CN104849278A CN104849278A CN201510260162.3A CN201510260162A CN104849278A CN 104849278 A CN104849278 A CN 104849278A CN 201510260162 A CN201510260162 A CN 201510260162A CN 104849278 A CN104849278 A CN 104849278A
- Authority
- CN
- China
- Prior art keywords
- pipe tobacco
- long filament
- likeness coefficient
- filament rate
- broken silk
- 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.)
- Granted
Links
Landscapes
- Manufacture Of Tobacco Products (AREA)
Abstract
The invention discloses a method for measuring the cut tobacco filament rate on the basis of a similarity coefficient. The method is applicable to the measurement of the cut tobacco filament rate of the same recipe ingredient, and has the advantages that the goal of measuring the filament rate at high precision, high accuracy and lower consumption is achieved, the important significance is realized on evaluating the reasonableness of the cigarette processing process and controlling the tobacco damage ratio in the cut tobacco making process, and good industrialized popularization and application prospects are realized.
Description
Technical field
The present invention relates to a kind of method measuring pipe tobacco long filament rate based on likeness coefficient, can in pipe tobacco process effective monitoring product quality; Belong to cigarette product quality monitoring technology field.
Background technology
In the production of cigarettes process of tobacco business and the quality testing at cigarette product, the detection of tobacco structure is all important link, wherein, the long filament rate of pipe tobacco clearly can reflect the quality control situation of pipe tobacco in process and make broken situation, also directly affects the quality index of seaming operation.
The method that tobacco structure is employing grab sample 4000g pipe tobacco, inquartation is reduced to the sorting of 1000g, YQ-2 type vibratory screening apparatus, sensibility reciprocal 1/10g balance is weighed is measured in current China tobacco business standard, the indivedual sample of minority is extracted by the every order of classes or grades at school of Quality Inspector, static state is carried out detection and is analyzed, and its result is as decision-making foundation afterwards.Certain deviation is there is between the tobacco structure data that this detection method obtains and the process that tobacco business requires at present accurately control, also there is following defect in this detection method: (1), measuring accuracy are lower, the balance system error of sensibility reciprocal 1/10g is larger, and the data obtained only retain a decimal, and the result of acquisition can not react the difference between different batches pipe tobacco; (2), sampling consumes greatly, and every batch of pipe tobacco normally samples 3-5 group, consumes pipe tobacco and is worth about 2000 yuan; (3), the accuracy of the method is lower, shake screening from time, exist and leak silk, carry the situations such as pipe tobacco secretly, cause pipe tobacco to be separated and there is deviation; (4), sampling risk is larger, be limited to sampling to consume, only getting in 5 minutes 3 parallel sample during each sampling Detection, for the whole production run production time of nearly 90 minutes, the sample size extracted and time period of covering less, the risk obtaining data is larger.
As can be seen here, set up that a kind of accuracy is high, the detection method of the pipe tobacco long filament rate of the high and low consumption of precision has very important significance.
Summary of the invention
For solving the deficiencies in the prior art, the object of the present invention is to provide a kind of method measuring pipe tobacco long filament rate based on likeness coefficient, there is high precision, high accuracy and low consumed advantage.
In order to realize above-mentioned target, the present invention adopts following technical scheme:
The method utilizing fractal theory, measure pipe tobacco long filament rate based on likeness coefficient of the present invention, comprises the steps:
S1, detection many groups tobacco structure data, obtain pipe tobacco long filament rate data;
S2, gather the broken silk image after the screening of each group of pipe tobacco;
S3, broken silk image to be analyzed, determine the pixel distribution scope of broken silk in image;
S4, carry out pixel range division according to pixel distribution scope, area percentage statistics is carried out to the broken silk component in the different pixels interval often organized in broken silk, and calculates cumulative percentage;
S5, linear fit is carried out to area logarithm and cumulative percentage logarithm, calculate the likeness coefficient of each group of broken silk;
S6, set up the fit equation of likeness coefficient and pipe tobacco long filament rate;
S7, get pipe tobacco to be measured and sieve, gather broken silk image and repeat step S4 and S5, obtaining the likeness coefficient of broken silk to be measured;
S8, the likeness coefficient obtained according to step S7, the fit equation of contrast step S6, obtains the long filament rate of pipe tobacco to be measured.
Preferably, in step sl, first utilize the sorting of YQ-2 type vibratory screening apparatus, then adopt sensibility reciprocal 1/10g balance to detect tobacco structure, tobacco structure can comprise long filament rate, short silk rate, filament broken rate etc., this step focus on the numerical value obtaining long filament rate.
Particularly, in step s 2, broken silk zero lap is laid on white background paper, the sampling amount of broken silk is preferably 1g, tiling area is 5cm*5cm, adopt digital camera (as Nikon D90 camera) to carry out image acquisition, pixel coverage is 6,000,000-1,800 ten thousand, and gathering distance is 10cm-100cm.
Preferably, in step s3, adopt the elemental area scope of MATLAB or IPP software to broken silk to add up, broken silk size represents with pixel, and the pixel coverage planted agent of division comprises the broken silk of more than 99%.
Usefulness of the present invention is: the assay method of pipe tobacco long filament rate of the present invention is applicable to the mensuration of same formula constituent pipe tobacco long filament rate, achieve high precision, high accuracy and low consumption low mensuration long filament rate, for evaluating cigarette process rationality, to control the broken rate of Primary Processing significant, has the good prospect of industrial application.
Accompanying drawing explanation
Fig. 1 is the area logarithm of one of them sample of the preferred embodiments of the present invention and the linear fit result of cumulative percentage logarithm;
Fig. 2 is the likeness coefficient of one of them sample of the preferred embodiments of the present invention and the linear fit result of long filament rate.
Embodiment
Below in conjunction with the drawings and specific embodiments, concrete introduction is done to the present invention.
The assay method of pipe tobacco long filament rate of the present invention realizes based on likeness coefficient, specific as follows:
(1), get certain batch of 24 groups, pipe tobacco before sieve is shaken in the outlet of product cigarette drying process, often organize 1000g, every 4 groups adopt inquartation to be reduced to 1000g, totally 6 parts;
(2), utilize the sorting of YQ-2 type vibratory screening apparatus, structure that sensibility reciprocal 1/10 balance detects 6 increment product pipe tobaccos, testing result is in table 1:
| Sample number into spectrum | Long filament rate (%) | Middle silk rate (%) | Whole cut rate (%) |
| 1 | 66.4 | 16.0 | 82.4 |
| 2 | 68.1 | 14.4 | 82.5 |
| 3 | 67.9 | 14.9 | 82.8 |
| 4 | 66.6 | 15.8 | 82.4 |
| 5 | 66.7 | 15.9 | 82.6 |
| 6 | 67.7 | 14.7 | 82.4 |
The tobacco structure testing result of table 1 sample 1-6
YQ-2 type vibratory screening apparatus has three-layer network, is sieved successively from top to bottom by pipe tobacco to be: long filament, middle silk, broken silk and broken end, in the present invention, long filament and middle silk are referred to as whole silk.
(3), get every increment product pipe tobacco screening after broken silk 1g, be laid on white background paper, tiling area be about 5cm*5cm, zero lap, carry out image acquisition with Nikon D90 camera, altogether gather 6 broken silk images;
(4), utilize the broken silk image of image analysis software to sample 1 to analyze, analyze the pixel coverage that broken silk comprises, determine afterwards by analysis to get comprise more than 99% broken silk pixel coverage 2-300 as the division scope of pixel range, in table 2:
| Minimum pixel | 2 | Pixel coverage | 2-30 | 2-60 | 2-90 | 2-120 | 2-150 | 2-180 | 2-210 | 2-240 | 2-300 |
| Maximum pixel | 327 | Area ratio | 10.23 | 29.52 | 43.59 | 57.98 | 70.59 | 81.02 | 86.43 | 93.79 | 99.12 |
The pixel range analysis result of table 2 sample 1
(5), in the present invention, the pixel coverage demarcation interval of same recipe pipe tobacco is identical, therefore, pixel range is divided according to the pixel coverage in table 2,6 broken silk images are analyzed, from analysis result, obtain the total pixel size shared by broken silk of respective area scope, analysis result is in table 3:
The broken silk image analysis result of pipe tobacco of table 3 sample 1-6
(6), according to the result of table 3, obtain area-area percentage result after conversion, in table, each numerical value is number percent, and for convenience of description, eliminate percentage sign " % ", result is as table 4:
Area-cumulative percentage the result of table 4 sample 1-6
Then, reference area logarithm and cumulative percentage logarithm, table 5 calculates for sample 1:
Lgy and the lgsi result of calculation of table 5 sample 1
Then, carry out linear fit to area logarithm and cumulative percentage logarithm, the fitting result of sample 1 is shown in Fig. 1, according to the slope of fit equation, and application of formula
calculate the likeness coefficient D of the broken silk of respective sample.
The computing method of sample 2-sample 6 are the same, are also application of formula
calculate the likeness coefficient D of the broken silk of respective sample, obtain the likeness coefficient D of 6 samples thus, the results are shown in Table 6:
The likeness coefficient D that table 6 sample 1-6 calculates
(7), by the long filament rate in each sample tobacco structure data in each sample likeness coefficient in table 6 and table 1 correspond to table 7, then carry out linear fit, obtain fit equation, fitting result is shown in Fig. 2,
| Sample number into spectrum | Likeness coefficient | Long filament rate (%) |
| 1 | 0.1052 | 66.4 |
| 2 | 0.083 | 68.1 |
| 3 | 0.086 | 67.9 |
| 4 | 0.1018 | 66.6 |
| 5 | 0.1012 | 66.7 |
| 6 | 0.0884 | 67.7 |
The likeness coefficient result corresponding to long filament rate of table 7 sample 1-6
(8), with 5 minutes intervals, get drying process and export the rear broken silk 10 groups of sieve screening that shakes, often organize about 1g, repeat (5) ~ (6) above, obtain the likeness coefficient organizing broken silk sample more; Finally, in the fit equation that the likeness coefficient obtained substitution step (7) is obtained, calculate the long filament rate often organizing the corresponding pipe tobacco of broken silk, the results are shown in Table 8.
| Batch | Likeness coefficient average | Long filament rate measured value (%) |
| 1 | 0.0952 | 67.2 |
| 2 | 0.0932 | 67.3 |
| 3 | 0.0783 | 68.5 |
| 4 | 0.1015 | 66.7 |
| 5 | 0.0899 | 67.6 |
| 6 | 0.1026 | 66.6 |
| 7 | 0.0898 | 67.6 |
| 8 | 0.1019 | 66.7 |
| 9 | 0.0944 | 67.2 |
| 10 | 0.0916 | 67.5 |
The likeness coefficient of 10 batches, table 8 and long filament rate result
To sum up, the assay method of pipe tobacco long filament rate of the present invention utilizes fractal theory, based on likeness coefficient, achieve high precision, high accuracy and low consumption low mensuration long filament rate, the pipe tobacco long filament rate being applicable to same formula constituent measures, for evaluating cigarette process rationality, to control the broken rate of Primary Processing significant, has the good prospect of industrial application.
More than show and describe ultimate principle of the present invention, principal character and advantage.The technician of the industry should understand, and above-described embodiment does not limit the present invention in any form, the technical scheme that the mode that all employings are equal to replacement or equivalent transformation obtains, and all drops in protection scope of the present invention.
Claims (10)
1. measure a method for pipe tobacco long filament rate based on likeness coefficient, it is characterized in that, comprise the steps:
S1, detection many groups tobacco structure data, obtain pipe tobacco long filament rate data;
S2, gather the broken silk image after the screening of each group of pipe tobacco;
S3, broken silk image to be analyzed, determine the pixel distribution scope of broken silk in image;
S4, carry out pixel range division according to pixel distribution scope, area percentage statistics is carried out to the broken silk component in the different pixels interval often organized in broken silk, and calculates cumulative percentage;
S5, linear fit is carried out to area logarithm and cumulative percentage logarithm, calculate the likeness coefficient of each group of broken silk;
S6, set up the fit equation of likeness coefficient and pipe tobacco long filament rate;
S7, get pipe tobacco to be measured and sieve, gather broken silk image and repeat step S4 and S5, obtaining the likeness coefficient of broken silk to be measured;
S8, the likeness coefficient obtained according to step S7, the fit equation of contrast step S6, obtains the long filament rate of pipe tobacco to be measured.
2. a kind of method measuring pipe tobacco long filament rate based on likeness coefficient according to claim 1, is characterized in that, in step S1, first utilize the sorting of YQ-2 type vibratory screening apparatus, then adopts sensibility reciprocal 1/10g balance to detect tobacco structure.
3. a kind of method measuring pipe tobacco long filament rate based on likeness coefficient according to claim 1, it is characterized in that, in step S2, the sampling amount of broken silk is 1g, is laid on white background paper by broken silk zero lap, and tiling area is 5cm*5cm.
4. a kind of method measuring pipe tobacco long filament rate based on likeness coefficient according to claim 1, is characterized in that, in step S2, the instrument gathering image is digital camera, and pixel coverage is 6,000,000-1,800 ten thousand, and gathering distance is 10cm-100cm.
5. a kind of method measuring pipe tobacco long filament rate based on likeness coefficient according to claim 4, it is characterized in that, in step S2, described digital camera is Nikon D90 camera.
6. a kind of method measuring pipe tobacco long filament rate based on likeness coefficient according to claim 1, is characterized in that, in step S3, adopt the elemental area scope of MATLAB or IPP software to broken silk to add up, broken silk area represents with pixel value.
7. a kind of method measuring pipe tobacco long filament rate based on likeness coefficient according to claim 6, it is characterized in that, in step S3, pixel distribution scope planted agent comprises the broken silk of more than 99%.
8. a kind of method measuring pipe tobacco long filament rate based on likeness coefficient according to claim 7, it is characterized in that, in step S4, same recipe pipe tobacco is suitable for same pixel interval.
9. a kind of method measuring pipe tobacco long filament rate based on likeness coefficient according to any one of claim 1-8, it is characterized in that, in step S4, likeness coefficient computing formula is:
wherein lgy is the logarithm of cumulative percentage, lgs
ifor the logarithm of area, lgs
maxfor the logarithm of the total area, D is likeness coefficient.
10. a kind of method measuring pipe tobacco long filament rate based on likeness coefficient according to claim 9, it is characterized in that, in step S6, the mensuration of same recipe pipe tobacco long filament rate is suitable for same fit equation.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510260162.3A CN104849278B (en) | 2015-05-20 | 2015-05-20 | A kind of method that pipe tobacco long filament rate is determined based on likeness coefficient |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510260162.3A CN104849278B (en) | 2015-05-20 | 2015-05-20 | A kind of method that pipe tobacco long filament rate is determined based on likeness coefficient |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104849278A true CN104849278A (en) | 2015-08-19 |
| CN104849278B CN104849278B (en) | 2017-06-09 |
Family
ID=53849084
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510260162.3A Active CN104849278B (en) | 2015-05-20 | 2015-05-20 | A kind of method that pipe tobacco long filament rate is determined based on likeness coefficient |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104849278B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106770303A (en) * | 2017-03-31 | 2017-05-31 | 河南农业大学 | Cigarette shreds structure characterization methods based on graphical analysis |
| CN109239082A (en) * | 2018-09-21 | 2019-01-18 | 杭州安脉盛智能技术有限公司 | Tobacco structure quality online test method and system based on machine vision technique |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102628813A (en) * | 2012-03-30 | 2012-08-08 | 中国烟草总公司郑州烟草研究院 | Measuring and representing method of slice tobacco raw material filamentation structure |
| CN102829730A (en) * | 2012-09-05 | 2012-12-19 | 江苏中烟工业有限责任公司 | Measuring method for shredding width |
| CN102879535A (en) * | 2012-10-25 | 2013-01-16 | 山东中烟工业有限责任公司 | Method for determining required tobacco length for improving cigarette quality stability |
| CN104237238A (en) * | 2014-10-11 | 2014-12-24 | 中国烟草总公司郑州烟草研究院 | Shredded tobacco structure prediction method |
-
2015
- 2015-05-20 CN CN201510260162.3A patent/CN104849278B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102628813A (en) * | 2012-03-30 | 2012-08-08 | 中国烟草总公司郑州烟草研究院 | Measuring and representing method of slice tobacco raw material filamentation structure |
| CN102829730A (en) * | 2012-09-05 | 2012-12-19 | 江苏中烟工业有限责任公司 | Measuring method for shredding width |
| CN102879535A (en) * | 2012-10-25 | 2013-01-16 | 山东中烟工业有限责任公司 | Method for determining required tobacco length for improving cigarette quality stability |
| CN104237238A (en) * | 2014-10-11 | 2014-12-24 | 中国烟草总公司郑州烟草研究院 | Shredded tobacco structure prediction method |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106770303A (en) * | 2017-03-31 | 2017-05-31 | 河南农业大学 | Cigarette shreds structure characterization methods based on graphical analysis |
| CN109239082A (en) * | 2018-09-21 | 2019-01-18 | 杭州安脉盛智能技术有限公司 | Tobacco structure quality online test method and system based on machine vision technique |
| CN109239082B (en) * | 2018-09-21 | 2021-01-26 | 杭州安脉盛智能技术有限公司 | Tobacco shred structure quality online detection method and system based on machine vision technology |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104849278B (en) | 2017-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104502245B (en) | Method for measuring fineness modulus of sand by utilizing image analysis technology | |
| CN110555596B (en) | Sampling inspection strategy making method and system based on power distribution material quality evaluation | |
| CN104330675B (en) | Multivariate time series based power transformation equipment online monitoring and analysis system and method thereof | |
| CN102095831A (en) | Moisture meter management method and system | |
| CN104425300B (en) | The method of sampling and device are measured in product | |
| CN109309022A (en) | A kind of defect sampling observation method | |
| CN106844901B (en) | Structural part residual strength evaluation method based on multi-factor fusion correction | |
| CN107358339A (en) | A kind of track quality state evaluating method based on big data | |
| CN112462312B (en) | Automatic calibration method for chip testing machine and application thereof | |
| CN104849278A (en) | Method for measuring cut tobacco filament rate on basis of similarity coefficient | |
| CN110261272A (en) | Based on geographical detection with PCA to the Key Influential Factors screening technique of PM2.5 concentration distribution | |
| CN108921305B (en) | Component life period monitoring method | |
| CN118051863B (en) | Health data acquisition system and method based on digital metering technology | |
| CN110458809B (en) | Yarn evenness detection method based on sub-pixel edge detection | |
| CN110186962B (en) | Incomplete measurement data imaging method for capacitance tomography | |
| CN107368303A (en) | A kind of evolution measure based on software architecture | |
| CN103902798B (en) | Data preprocessing method | |
| CN111122775A (en) | Pollution concentration monitoring equipment-oriented segmentation data calibration method and system | |
| CN105628646B (en) | A kind of online tar prediction of cigarette and method for early warning | |
| CN105844646A (en) | Mixture spectrum consistency determination method | |
| CN204795439U (en) | A predistortion test card molding system for image acquisition system that distorts greatly | |
| CN111337452A (en) | Method for verifying feasibility of spectral data model transfer algorithm | |
| CN107271780A (en) | A kind of abnormal detection method of biological tissue's electrical impedance characteristic parameter | |
| CN101592692B (en) | Evaluation method of measuring machines | |
| CN107917914B (en) | Method for evaluating mixing uniformity of formula components in formula tobacco shreds in cigarette processing industry |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| EXSB | Decision made by sipo to initiate substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |