测定部11在印刷介质通过印刷头33的下侧的状态下开始测定。测定部11使光源部42亮灯。另外，测定部11使滑架32动作而使受光传感器41向头扫描方向移动。并且，测定部11分别在头扫描方向的各个头位置，控制受光传感器41，使多个波长光由分光部411进行透射，同时从受光部412输出该多个波长的光的电信号。并且，测定部11在头扫描方向的多个头位置中的每个头位置，取得从受光传感器41输出的多个波长的光的电信号(光强度)。并且，测定部11在头扫描方向的多个头位置中的每个头位置，通过存储在存储部20等的基准板43的各波长的光的强度对取得的各波长的光强度进行除算，算出相当于各波长的光的反射率的值，并且存储于存储部20等。由此，在多个头位置，都能够测定各波长的光的反射率(分光反射率)。
[0121]在步骤S21中，针对一条线，进行多个头位置的反射率的测定。当然，也可以针对多条线，进行多个头位置的反射率的测定。
[0122]接着，印刷介质种类判定部13基于在步骤S21测定的分光反射率、以及保存于印刷介质种类数据存储部21的印刷介质种类数据，判别印刷介质的印刷介质种类(步骤S22)。下面，进行详细的说明。
[0123]印刷介质种类判定部13算出步骤S21测定的分光反射率与各样品(各印刷介质种类)之间的马氏距离。在本实施方式中，由于在印刷介质种类数据存储部21存储有各样品的印刷介质种类数据(平均分光反射率及协方差矩阵)，所以可以算出这些印刷介质种类数据的几分的马氏距离。
[0124]图11是马氏距离的说明图。假设以已知某组(此处设为Gr.A)有关的平均值a与另一组(此处设为Gr.B)有关的平均值b的情况下，得到了有关新的样本的测定值X的情况为例。如图11的(a)所示，如果测定值X相比平均值b更接近于平均值a的话(偏差很小的情况下)，则该样本通常认为是属于Gr.A0然而，这是以在Gr.A与Gr.B测定值的偏差是相同程度为前提的基础上成立的。因此，考虑到Gr.A及Gr.B的测定值的偏差，也可能引起结论相反的情况。
[0125]例如，如图11的(b)所示的情况下，即使样本的测定值X相比于平均值b更接近于平均值a，但是也考虑该样本属于Gr.B。S卩，的确是测定值x与平均值a的偏差是比测定值X与平均值b的偏差小，但认为该样本属于Gr.A的话，测定值X与平均值a的偏差与Gr.A的测定值的偏差相比过大。对此，由于测定值X与平均值b的偏差比Gr.B的测定值的偏差小，所以自然认为该样本属于Gr.B。由此，不仅考虑测定值与平均值的偏差，也可通过考虑相对于测定值的偏差(方差)的偏差的比例，从而能够更加正确地判断。马氏距离是指，在考虑如上测定值的偏差的基础上，示出认为样本是属于哪组的指标。属于某组的可能性越高，则相对于该组的马氏距离就越小。
[0126]如图11的(b)所示，测定值X是一维时，测定值与平均值的偏差进行平方，通过方差对该值进行除算，能够算出马氏距离(正确为马氏距离的平方值)。在图11的(c)中，示出求得马氏距离的计算式。计算式中的X表示测定值，av表示为平均值，s表示为方差。
[0127]另外，马氏距离也能够扩展为多维。对作为最简单的情况的二维马氏距离进行说明。在二维时，每次测定得到xl、x2这两个测定值。另外，有关方差，不仅有关Xl的方差、有关x2的方差，也能够考虑xl与x2之间的协方差。在此，将如图11的(c)所示的计算式的测定值X通过向量(xl，x2)置换，根据由协方差矩阵置换的方差S，从而能够得到求得如图11的(d)所示的二维的马氏距离的计算式。另外，计算式中的avl、av2分别表示有关xl、x2的平均值。另外，协方差矩阵的右边带有的“-1”表示为逆矩阵。并且，带有向量(xl，x2)的“T”表示换位向量。
[0128]如图11的(d)所示的算出式，将向量(xl-avl.x2-av2)以大写字母的“X”表示，如果协方差矩阵以大写字母“R”表示，则二维的马氏距离的算出式能够以图11的(e)表示。该图11的(e)的算出式能够直接作为三维以上的多维马氏距离的算出式使用。S卩，当算出η维的马氏距离的情况下，成为“X”具有η各成分的向量，“R”成为η行η列的协方差矩阵。
[0129]在图7的步骤S21得到的分光反射率是通过各头位置P (I?P的范围、P为I以上的自然数)取得η个波长的、ηΧΡ个波长的反射率，所以是ηΧρ维的测定值(参考图8)。另外，从这种的分光反射率求得的协方差矩阵成为ηΧΡ行ηΧΡ列的矩阵(参照图10)。并且，在印刷介质种类数据存储部21保存有与多个样本相关的平均分光反射率及协方差矩阵。在此，在步骤S22中，使用图11的(e)的算出式算出由步骤S21测定的分光反射率开始的各样本的马氏距离。并且，选择马氏距离为最小值的样品(印刷介质种类)，并判别该印刷介质种类是进行印刷的印刷介质的印刷介质种类。
[0130]另外，在步骤S22中，虽然选择马氏距离是最小的印刷介质种类，但是也可以例如，印刷介质种类判定部13判定该最小的马氏距离是否超过规定的阈值，当超过时，经由通信部50在能通信的PC等终端、操作面板60上输出并显示未登记合适的印刷介质种类、催促注册印刷介质种类数据及印刷参数的信息的内容。
[0131]接着，印刷控制部10从印刷参数存储部22取得与通过步骤S22判别的印刷介质种类对应的印刷参数(步骤S23)。并且，印刷控制部10取得存储于存储部20等的印刷对象数据(步骤S24)。印刷对象数据是具有例如R(红色)、G(绿色)、B(蓝色)各色的数据的图像数据。
[0132]然后，印刷控制部10对由步骤S24取得的印刷对象数据，基于由步骤S23取得的印刷参数，执行各种的图像处理(步骤S25)。具体而言，图像处理部14基于由步骤S23取得的印刷参数(例如，有关颜色转换处理、半色调处理、墨水喷出量、印刷介质输送处理等的参数)，对印刷对象数据进行图像处理，并生成印刷数据。印刷数据是例如标示打开及关闭墨水喷出的点数据。
[0133]并且，印刷控制部10基于由步骤S25生成的印刷数据，控制印刷发动机部30，执行向印刷介质的印刷(步骤S26)。另外，步骤S6的处理也可以基于印刷参数(例如，与印刷介质输送处理等相关的参数)而执行。另外，印刷控制部10进行干燥处理(步骤S27)。具体而言，干燥处理部15基于由步骤S23取得的印刷参数(与干燥处理相关的参数)，控制加热器34，并在步骤S26使印刷结束的地方干燥。
[0134]接着，印刷控制部10判定规定范围(例如一页)的印刷是否结束(步骤S28)。当未结束印刷时(步骤S28的否)，印刷控制部10再一次进行S25的处理。
[0135]另外，在本实施方式中，印刷控制部10控制输送部31输送印刷介质，并进行步骤S25?步骤S28的处理。
[0136]当印刷结束时(步骤S28的是)，印刷控制部10判定是否存在其他的印刷(例如下一页)(步骤S29)。在存在其他的印刷时(步骤S29的是)，印刷控制部10返回处理至步骤S20。没有其他的印刷时(步骤S29的否)，印刷控制部10结束如图7所示的流程图的处理。
[0137]另外，为了容易理解印刷装置I的处理，根据主要的处理内容，将如上述的图6及图7所示的流程图的各处理单位进行了分割。处理单位的分割的方法、名称并不限制本发明。印刷装置I的处理根据处理内容能够进一步分割为更多的处理单位。另外，能够分隔使一个处理单位进一步包括更多的处理。
[0138]以上，对本发明的一实施方式进行了说明。在本实施方式中，在图像的印刷之前判别印刷介质的种类，并根据判别结果印刷图像。因此，印刷装置的用户即使未设定印刷介质的种类，也能够自动判别印刷介质的种类而适当印刷图像。并且，可以自动进行印刷介质的印刷参数的设定等。
[0139]另外，当判别印刷介质的种类时，预先存储有关多个样本的平均分光反射率及协方差矩阵，通过相对于各样本算出马氏距离，从而判别印刷介质的种类。如果使用协方差，不仅可以考虑各波长的测定值的偏差(方差)，也可以考虑多个波长间的关系(分布状况等)。因此，如果基于使用平均分光反射率及协方差矩阵算出的马氏距离来判别印刷介质的种类，则在不仅考虑各波长的反射率也考虑各波长间的反射率的关系的基础上，能够选择最接近测定的分光反射率的样本。由此，可以通过非常高的概率正确地判别印刷介质的种类。其结果，由于印刷装置的用户弄错印刷介质的种类的设定、或由于忘记设定而成为以前的设定等因素，在弄错印刷介质的种类并设定的情况下不印刷图像，从而能印刷适当地图像。
[0140]当然，测定了分光反射率的印刷介质也可以使用没有存储数据作为样本的种类的印刷介质(未知的印刷介质)。此种情况时，可以判别与实际不同种类的印刷介质。然而，判断与马氏距离最小的样品为相同种类是指判断与该印刷介质最接近的样品。因此，即使是未知的印刷介质，由于印刷装置也在可以的范围内进行认为是最适当地图像处理，所以可以适当地印刷图像。
[0141]另外，在实施方式中，对于多个头位置进行各波长的反射率的测定，基于该测定结果，算出平均反射率及协方差矩阵，从而作为印刷介质种类数据进行保存。并且，在实际的印刷中，使用对应多个头位置测定的各波长反射率、被保存的各样品的平均反射率及协方差矩阵，从而进行印刷介质种类的判别。由此，通过使用多个头位置的测定结果，从而由于考虑多个不同头位置的反射光(正反射光及漫反射光)的马氏距离被求得，所以能够进行更高精度的印刷介质种类判别。即，由于不仅印刷介质的颜色，且厚度、材质、光泽(表面状态)等的性质的辨别精度得以提高，所以即使性质具有微妙的不同的情况下，也能够更准确地判别印刷介质种类，并对该印刷介质种类选择适当地印刷参数。
[0142]另外，在本实施方式中，在作为可动部分的印刷头设置受光传感器，通过在印刷装置的非可动部分设置光源，从而利用印刷头的扫描能够简单测定多个头位置的光强度。另夕卜，由于可以向印刷装置追加受光传感器及光源，所以能够抑制硬件的制造成本。另外，通过使用法布里珀罗光谱仪，从而能够使受光传感器更加小型化。
[0143]另外，在本实施方式中，作为印刷参数，准备与颜色转换处理、半色调处理、墨水喷出量、印刷介质输送处理、干燥处理等相关的各种参数，基于这些参数进行印刷。如此，能够减少根据印刷介质的种类设定多个印刷参数的用户的负担。
[0144]〈变形例〉
[0145]在上述的实施方式中，直接使用测定的分光反射率而算出马氏距离。测定出的分光反射率是多个波长的反射率，所以通过多维算出马氏距离。如果减少测定的波长的数量，则测定变快，虽然马氏距离的计算也变得快速，但是一旦测定波长的数量减少，则担心印刷介质的判别精度将恶化。在此，不仅在直接使用测定的分光反射率而算出马氏距离，还使用主成分分析减少分光反射率的维数的基础上，通过少量维数算出马氏距离来判别印刷介质的种类。
[0146]在上述实施方式中，虽然已经对将波长的数量设为nXP个的情况进行说明，在下面，以n = n X P为例进行说明。即，测定值为η维，协方差矩阵为η行η列。
[0147]首先，开始对使用主成分分析减少维数的方法进行说明。另外，由于主成