延长点污棉类国际分级实物标准使用期限的计算方法与流程

本发明涉及出入境棉花检测技术领域，特别涉及一种延长点污棉类国际分级实物标准使用期限的计算方法。

背景技术：

国际棉花分级通用实物标准有效期为一年，使用期限达到一年后全部更换，但并不是所有的实物标准均不能够再使用，因此一次性全部更换势必会造成资源浪费和经济损失，给国家带来严重的损失。

此外，分布于我国各地的棉花品级实验室温湿度存在很大差异，各实验室都是将标准物质置于自然环境中，不同的环境条件必然使标准物质产生不同程度的变化，标准物质发生不同程度的变化必将带来不同的实验室检测结果不一致的现象。多年来，随着我国棉花贸易范围的扩大，不同地区的检验检疫机构对同一批进口棉花的品级检验结果不一致的情况也非常普遍，引起国内外贸易双方的异议，从而导致贸易纠纷，检验检疫机构的公正性和准确性受到质疑，并带来很大的负面影响。

技术实现要素：

有鉴于此，本发明旨在提出一种延长点污棉类国际分级实物标准使用期限的计算方法，以可减少点污棉类国际通用棉花实物标准的更换次数，从而延长使用时间。

为达到上述目的，本发明的技术方案是这样实现的：

一种延长点污棉类国际分级实物标准使用期限的计算方法，该方法包括对点污棉类国际棉花分级实物标准各级标准单独定制使用期限，若该级标准达到使用期限则进行单独更换，具体步骤包括：

s1，将点污棉类国际棉花分级实物标准进行恒温恒湿保存，并对第一年的点污棉类国际棉花分级实物标准进行hvi测试，以确定点污棉类国际棉花分级实物标准第一年的反射率rd和和黄色深度+b，并分别记为rd初值和+b初值。

s2，将+b初值作为y值带入黄色深度公式y＝a1x⁴+a2x³+a3x²+a4x+m中，其中x为年份、a1、a2、a3、a4分别为与点污棉类中不同色征级棉花相对应的定值，计算得出m值，然后将m值带入公式以计算未来各年份所对应y值，该y值即为与点污棉类各色征级相对应的未来各年份的+b值；

s3，将rd初值作为y值带入公式y＝b1x⁴+b2x³+b3x²+b4x+n中，其中x为年份、b1、b2、b3、b4分别为与点污棉类中不同色征级棉花相对应的定值，计算得出n值，然后将n值带入公式以计算未来各年份所对应y值，该y值即为与点污棉类各色征级相对应的未来各年份的rd值；

s4，通过未来各年份的rd值和+b值判断点污棉类国际棉花分级实物标准继续使用或进行更换。

进一步的，步骤s3的具体步骤为：

根据得到的rd值或+b值，在棉花色征图的纵坐标和横坐标上分别找出与rd值和+b值相对应的位置点，若该位置点未超出其在棉花色征图上的规定区域，则该级标准继续使用，若该位置点超出其在棉花色征图上的规定区域，则判定该级标准失效，不再使用。

进一步的，在步骤s1中所述恒温恒湿环境为温度20±1℃，湿度65±2rh％。

进一步的，所述点污棉类国际分级实物标准使用期限≤7年。

相对于现有技术，本发明具有以下优势：

本发明所述的延长点污棉类国际分级实物标准使用期限的计算方法，通过将各级标准单独更换而不是统一规定一年全部更换，不仅延长了使用时间，而且避免因一次性全部更换而带来的资源浪费，从而减少资源浪费，节约开支；而采用rd值和+b值来判定是否达到使用期限的方法可减少人工判定误差，保证判定的准确性，且通过公式计算以直接得到未来年份的rd值和+b值，可对实物标准的有效性进行量化判断，保证检测结果准确，同时从而不需要每年对实物标准进行hvi测试，减少了繁琐的试验步骤，提高工作效率。

附图说明

构成本发明的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1为本发明实施例所述的美国陆地棉棉花色征图；

具体实施方式

需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。

下面将参考附图并结合实施例来详细说明本发明。

本实施例涉及一种延长点污棉类国际分级实物标准使用期限的计算方法，该方法包括：对国际通用棉花实物标准各级标准单独定制使用期限，若该级标准达到使用期限则进行单独更换。本发明所述的延长点污棉类国际分级实物标准使用期限的计算方法，通过将各级标准单独更换而不是统一规定一年全部更换，不仅延长了使用时间，而且避免因一次性全部更换而带来的资源浪费，从而减少资源浪费，节约开支。

上述的描述中，具体步骤包括：

s1，将点污棉类国际棉花分级实物标准在温度20±1℃、湿度65±2rh％的恒温恒湿保存，并对第一年的点污棉类国际棉花分级实物标准进行hvi测试，以确定点污棉类国际棉花分级实物标准第一年的反射率rd和黄色深度+b，并分别记为rd初值和+b初值。

s2，将+b初值作为y值带入黄色深度公式y＝a1x⁴+a2x³+a3x²+a4x+m中，其中x为年份、a1、a2、a3、a4分别为与点污棉类中不同色征级棉花相对应的定值，计算得出m值，然后将m值带入公式以计算未来各年份所对应y值，该y值即为与点污棉类各色征级相对应的未来各年份的+b值；

s3，将rd初值作为y值带入公式y＝b1x⁴+b2x³+b3x²+b4x+n中，其中x为年份、b1、b2、b3、b4分别为与点污棉类中不同色征级棉花相对应的定值，计算得出n值，然后将n值带入公式以计算未来各年份所对应y值，该y值即为与点污棉类各色征级相对应的未来各年份的rd值。

s4，通过未来各年份的rd值和+b值判断点污棉类国际棉花分级实物标准继续使用或进行更换，该步骤具体包括：根据得到的rd值或+b值，在棉花色征图的纵坐标和横坐标上分别找出与rd值和+b值相对应的位置点，若该位置点未超出其在棉花色征图上的规定区域，则该级标准继续使用，若该位置点超出其在棉花色征图上的规定区域，则判定该级标准失效，不再使用。

国际通用棉花分级实物标准主要包括陆地棉和皮马棉，由于进口皮马棉很少，因此皮马棉标准使用率也很低，最常用的是陆地棉国际通用棉花分级实物标准。陆地棉分级实物标准主要包括白棉、点污棉、淡黄染棉3类共15盒实物标准，其余都是文字标准。其中白棉实物标准包括gm、sm、m、slm、lm、sgo、go共7盒色征级，点污棉类实物标准包括sm.sp、m.sp、slm.sp、lm.sp、sgo.sp共5盒色征级，淡黄染类实物标准包括sm.sp.tg、m.sp.tg、lm.sp.tg共3盒色征级。基于如上的描述，本实施例中以点污棉(spotted)为例进行详细说明，上述样品的选取具备代表性且能覆盖主要的点污棉类棉花品级实物标准类型。

国际通用棉花实物标准中，在点污棉类实物标准七盒色征级sm.sp、m.sp、slm.sp、lm.sp、sgo.sp中，其中每盒色征级中包含六块不同产地的样品，依次编号为：

sm.sp盒色征级中六块样品编号分别为23-1、23-2、23-3、23-4、23-5、23-6；

m.sp盒色征级中六块样品编号分别为33-1、33-2、33-3、33-4、33-5、33-6；

slm.sp盒色征级中六块样品编号分别为43-1、43-2、43-3、43-4、43-5、43-6；

lm.sp盒色征级中六块样品编号分别为53-1、53-2、53-3、53-4、53-5、53-6；

sgo.sp盒色征级中六块样品编号分别为63-1、63-2、63-3、63-4、63-5、63-6。

表1点污棉标准样品黄色深度测试数据

表2点污棉标准样品反射率测试数据

表3点污棉标准样品黄色深度和反射率计算公式

表1和表2分别给出了点污棉标准样品黄色深度测试数据和点污棉标准样品反射率测试数据，表3给出了点污棉各盒色征级的六个样品的黄色深度和反射率计算公式，下文以2011年购买的点污棉类国际通用棉花实物标准并结合表1至表3对步骤s2进行详细说明。

如点污棉sm.sp23-1的黄色深度+b的计算公式中a1＝-0.0006、a2＝0.0184、a3＝-0.0804、a4＝0.1292，带入公式得其黄色深度计算公式为y＝-0.0006x⁴+0.0184x³-0.0804x²+0.1292x+m，经测2011年刚购买的点污棉sm.sp23-1黄色深度+b初值为11.03，x＝1，带入公式得m＝10.9634，即sm.sp23-1未来年份x黄色深度的计算公式为y＝-0.0006x⁴+0.0184x³-0.0804x²+0.1292x+10.9634，2≤x。由该公式计算得点污棉sm.sp23-1未来各年份的黄色深度值见表4(由公式计算的数值均采用四舍五入的方法保留两位小数点，下同)，由表4可知经由公式计算得出的黄色深度值与hvi的测试值相等，且在下文判断能否继续使用的结果与hvi测试数据的判断结果完全相同，可作等同处理。

表4点污棉sm.sp23-1由公式计算得到的未来年份黄色深度数据与hvi测试值对比

点污棉sm.sp23-1反射率的计算公式中，b1＝0.0007、b2＝-0.0059、b3＝-0.01、b4＝-0.0363，带入公式得反射率计算公式为y＝0.0007x⁴-0.0059x³-0.01x²-0.0363x+n，经测2011年刚购买的点污棉sm.sp23-1的反射率rb初值为74.44，x＝1，带入公式得n＝74.4915，即点污棉sm.sp23-1未来年份x反射率的计算公式为y＝0.0007x⁴-0.0059x³-0.01x²-0.0363x+74.4915，2≤x。该公式计算得点污棉sm.sp23-1未来各年份的黄色深度值见表5，由表5可知经由公式计算得出的反射率与hvi的测试值相等或者近似相等，且在下文判断能否继续使用的结果与hvi测试数据的判断结果完全相同，可作等同处理。

表5点污棉sm.sp23-1由公式计算得到的未来年份反射率数据与hvi测试值对比

如点污棉lm.sp53-1的黄色深度+b的计算公式中，a1＝0.0006、a2＝0.0084、a3＝-0.0062、a4＝0.016，带入公式得其黄色深度计算公式为y＝-0.0006x⁴+0.0084x³-0.0062x²-0.016x+m，经测2011年刚购买的点污棉lm.sp53-1黄色深度+b初值为10.87，x＝1，带入公式得m＝10.8844，即点污棉lm.sp53-1未来年份x黄色深度的计算公式为y＝-0.0006x⁴+0.0084x³-0.0062x²-0.016x+10.8844，2≤x。由该公式计算得点污棉lm.sp53-1未来各年份的黄色深度值见表6，由表6可知经由公式计算得出的黄色深度值与hvi的测试值近似相等，且在下文判断能否继续使用的结果与hvi测试数据的判断结果完全相同，可作等同处理。

表6点污棉lm.sp53-1由公式计算得到的未来年份黄色深度数据与hvi测试值对比

如点污棉lm.sp53-1的反射率rd的计算公式中，b1＝-0.001、b2＝0.0216、b3＝-0.1546、b4＝0.3249，带入公式得其黄色深度计算公式为y＝-0.001x⁴+0.0216x³-0.1546x²+0.3249x+n，经测2011年刚购买的点污棉lm.sp53-1反射率rd初值为60.65，x＝1，带入公式得n＝60.4591，即点污棉lm.sp53-1未来年份x反射率的计算公式为y＝-0.001x⁴+0.0216x³-0.1546x²+0.3249x+60.4591，2≤x。由该公式计算得点污棉lm.sp53-1未来各年份的黄色深度值见表7，由表7可知经由公式计算得出的黄色深度值与hvi的测试值相等或者近似相等，且在下文判断能否继续使用的结果与hvi测试数据的判断结果完全相同，可作等同处理。

表7点污棉lm.sp53-1由公式计算得到的未来年份反射率数据与hvi测试值对比

上述描述通过2011刚购买的点污棉为例并以点污棉sm.sp23-1、lm.sp53-1进行举例说明，描述了对点污棉各色征级黄色深度和反射率的m.sp值和n值的确定步骤，点污棉其它色征级的m.sp值和n值的确定方法与上述相同，在此则不再一一详细描述，表8给出了以2011年购买的点污棉为基准以表3给出的公式通过上述方法得出的计算公式。

表82011年点污棉国际棉花实物标准各色征级黄色深度和反射率计算公式

下文将结合表1至表3以及图1对步骤s3如何通过rd值和+b值判断点污棉国际棉花分级实物标准继续使用或进行更换进行详细说明。

点污棉类其在棉花色征图上对应的区域为spotted区域，sm.sp、m.sp、slm.sp、lm.sp、sgo分别对应于spotted-sm.sp、spotted-m.sp、spotted-slm.sp、spotted-lm.sp、spotted-sgo。如点污棉sm.sp23-1在2011年的反射率和黄色深度hvi测试数据分别为(74.44，11.03)，其在图1(图1中纵坐标为反射率rd，横坐标为黄色深度+b)中的交叉点位于spotted-sm.sp区域内，满足要求继续使用；由公式计算得到的点污棉sm.sp23-1在2012至2016年的反射率和黄色深度分别为(74.34，11.04)、(74.19，11.08)、(73.99，11.22)、(73.76，11.52)、(73.55，12.04)，其在图1中的交叉点均位于spotted-sm.sp区域内，满足要求继续使用；由公式计算得到的点污棉sm.sp23-1在2017年的反射率和黄色深度分别为(73.40，12.80)，其在图1中的交叉点则不再位于spotted-sm.sp区域内，不能满足要求，因此不再使用。

如点染棉lm.sp53-1在2011年的反射率和黄色深度hvi测试数据分别为(60.65，10.87)，其在图1中的交叉点位于spotted-lm.sp区域内，满足要求继续使用；由公式计算得到的点染棉lm.sp53-1在2012至2015年的反射率和黄色深度hvi测试数据分别为(60.65，10.89)(60.54,10.96)、(60.41,11.11)、(60.29,11.32)，其在图1中的交叉点均位于spotted-sm.sp区域内；由公式计算得到的点污棉lm.sp53-1在2016年的反射率和黄色深度分别为(60.21，11.60)，其在图1中的交叉点则不再位于spotted-lm.sp区域内，不能满足要求，因此不再使用。

在对步骤s3如何通过rd值和+b值判断点污棉国际棉花分级实物标准继续使用或进行更换时，由于部分黄色深度和反射率由公式得出的计算数据与hvi测试数据部分近似相等，为保证公式的准确性，下文将通过hvi测试数据进行再次验证。

由hvi测试数据得到的点污棉sm.sp23-1在2012至2016年的反射率和黄色深度分别为(74.34，11.04)、(74.20，11.07)、(73.96，11.20)、(73.77，11.50)、(73.51，11.97)，其在图1中的交叉点均位于spotted-sm.sp区域内，满足要求继续使用；点污棉sm.sp23-1在2017年的反射率和黄色深度hvi测试数据分别为(73.36，12.68)，其在图1中的交叉点则不再位于spotted-sm.sp区域内，与公式计算得出的点污棉sm.sp23-1使用期限结论完全相同，由此可见公式满足使用要求。

由hvi测试数据得到的点染棉lm.sp53-1在2012至2015年的反射率和黄色深度hvi测试数据分别为(60.64，10.89)(60.57,10.95)、(60.41,11.09)、(60.28,11.30)，其在图1中的交叉点均位于spotted-lm.sp区域内，满足要求继续使用；点染棉lm.sp53-1在2016年的反射率和黄色深度hvi测试数据分别为(60.21，11.53)，其在图1中的交叉点则不再位于spotted-lm.sp区域内，与公式计算得出的点染棉lm.sp53-1使用期限结论完全相同。

同时，由上述描述可知，点污棉类国际分级实物各色征级标准使用期限不等，因七年后点污棉各色征级反射率和黄色深度变化规律不确定性较高，因此各点污棉类各国际分级实物标准使用期限≤7年。

本发明所述的延长点污棉类国际分级实物标准使用期限的计算方法，通过将各级标准单独更换而不是统一规定一年全部更换，不仅延长了使用时间，而且避免因一次性全部更换而带来的资源浪费，从而减少资源浪费，节约开支；而采用rd值和+b值来判定是否达到使用期限的方法可减少人工判定误差，保证判定的准确性，且通过公式计算以直接得到未来年份的rd值和+b值，可对实物标准的有效性进行量化判断，保证检测结果准确，同时从而不需要每年对实物标准进行hvi测试，减少了繁琐的试验步骤，提高工作效率。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。