颜色分解装置及颜色分解方法
【专利摘要】本发明提供颜色分解装置及颜色分解方法，能够容易地进行用于防止莫尔条纹、罗塞塔图案等图像结构产生的颜色分解。设备信号值备选计算部(106)基于取得的印刷轮廓，算出使在记录纸上重现的目标颜色所对应的表色系的值在容许色差的范围内能够重现的设备信号值的备选。颜色评价值保持部(110)保持设备信号值的各备选的色差作为颜色评价值。画质评价值计算部(116)使用网点重叠模拟的结果来算出各设备信号值备选的画质评价值。综合评价值计算部(118)基于上述颜色评价值及画质评价值，算出各设备信号值备选的综合评价值。最佳信号值决定部(120)基于该综合评价值从设备信号值备选之中决定最佳的设备信号值。
【专利说明】颜色分解装置及颜色分解方法
【技术领域】
[0001] 本发明涉及颜色分解装置、颜色分解方法及颜色分解程序，尤其是涉及从印刷用 的图像数据颜色分解成印刷装置所使用的颜色的颜色分解装置、颜色分解方法及颜色分解 程序。
【背景技术】
[0002] 在对图像数据进行印刷时，将图像数据颜色分解成各色材（例如，印刷色（CMYK)) 的网点图像，使用形成有与各色材对应的网点的印刷版来进行重叠印刷，由此来表现图像 数据的灰度。如此在使与多个色材对应的网点图像重合的情况下，根据网点的位置、角度、 网线数、形状，有时会产生具有莫尔条纹、罗塞塔图案等图像结构的花纹。
[0003] 在专利文献1中公开了一种在对图像数据进行颜色分解时对网点图像中的具有 强周期性的噪声信息进行评价的方法。在专利文献1中，读入印刷的块图像，将该图像数据 的各像素数据分解成RGB数据，变换成均等颜色空间的颜色信息（例如，L*、a*、b*)而向最 佳化部12转送。并且，最佳化部12基于由n部11选择的粒状度、油墨浓度的参数，并基 于L*，算出粒状性评价值（段落[0019]?[0023])。
[0004] 专利文献
[0005] 专利文献1 :日本特开2007-116287号公报
【发明内容】
[0006] 发明要解决的课题
[0007] 根据专利文献1，在进行颜色分解时，对与全部设备信号值（每一种颜色为256 色）的组合对应的块进行印刷，从该块的图像获得颜色信息来算出粒状度的评价值。因此， 在进行向两种颜色的颜色分解的情况下，需要对256X256 ( = 65536)这样的多个块进行印 刷而测定。在进行向CMYK的四种颜色的颜色分解的情况下，对庞大的数目的块进行印刷而 测定，因此耗费劳力。而且，粒状度等的评价值也根据网点的种类而不同，因此需要对应各 网点的种类进行印刷、测定的作业，因此负荷大。
[0008] 本发明鉴于这样的情况而作出，目的在于提供能够容易地进行用于防止莫尔条 纹、罗塞塔图案等图像结构的产生的颜色分解的颜色分解装置、颜色分解方法及颜色分解 程序。
[0009] 用于解决课题的方案
[0010] 为了解决上述课题，本发明的第一方案的颜色分解装置具备：目标值取得单元，在 根据连续色调图像数据而生成在印刷装置中表示构成图像的网点的形状及配置的各色材 的二值图像数据时，取得在印刷装置中重现的颜色的目标值；网点阈值数据取得单元，取得 包含用于将连续色调图像数据变换成各色材的二值图像数据的各网点的阈值的信息的网 点阈值数据；印刷轮廓取得单元，取得表示印刷装置中的设备信号值与表色系的值之间的 对应关系的印刷轮廓；及颜色分解单元，基于由目标值取得单元取得的颜色的目标值和印 刷轮廓，在印刷装置中算出设备信号值的备选，基于网点阈值数据及印刷轮廓，从设备信号 值的备选之中决定用于使目标值所对应的颜色重现的设备信号值。
[0011] 根据本方案，无需进行各设备信号值的块的印刷及测定，通过更简易的运算，能够 实现可抑制色差及抑制图像结构的产生的最佳的颜色分解。
[0012] 本发明的第二方案的颜色分解装置以第一方案为基础，颜色分解装置还具备：颜 色评价值保持单元，对应设备信号值的各备选，保持表示与目标值所对应的颜色之间的色 差的颜色的评价值；及画质评价值计算单元，基于网点阈值数据及印刷轮廓，对应设备信号 值的各备选而算出画质的评价值，颜色分解单元基于颜色的评价值及画质的评价值来决定 用于使目标值所对应的颜色重现的设备信号值。
[0013] 本发明的第三方案的颜色分解装置以第二方案为基础，颜色分解装置还具备模拟 单元，该模拟单元使用网点阈值数据对设备信号值的备选实施二值化处理，对于各备选生 成各色材的二值图像数据，使各色材的二值图像数据重叠，由此进行由印刷装置形成于印 刷物上的图像的模拟，画质评价值计算单元基于形成于印刷物上的图像的模拟的结果，对 应各备选而算出画质的评价值。
[0014] 本发明的第四方案的颜色分解装置以第三方案为基础，模拟单元在印刷装置中， 基于形成了二值图像数据所对应的网点时的响应特性，考虑由印刷装置印刷网点时产生的 网点的形状的变化来生成二值图像数据。
[0015] 本发明的第五方案的颜色分解装置以第二至第四方案中任一方案为基础，颜色分 解装置还具备油墨成本评价值计算单元，该油墨成本评价值计算单元对于设备信号值的各 备选，基于各色材的网点面积率来算出油墨成本的评价值，颜色分解单元基于颜色的评价 值、画质的评价值及油墨成本的评价值来决定用于使目标值所对应的颜色重现的设备信号 值。
[0016] 本发明的第六方案的颜色分解装置以第二至第五方案中任一方案为基础，颜色分 解装置还具备：加权系数设定单元，对各评价值分别设定加权系数；及综合评价值计算单 元，通过使用由加权系数设定单元设定的加权系数进行各评价值的加权相加，对应设备信 号值的各备选而算出综合评价值，颜色分解单元基于综合评价值来决定用于使目标值所对 应的颜色重现的设备信号值。
[0017] 本发明的第七方案的颜色分解装置以第一至第六方案中任一方案为基础，目标值 取得单元通过接受用于确定颜色的标识符的输入或测定颜色，取得在印刷装置中重现的颜 色的目标值。
[0018] 本发明的第八方案涉及一种颜色分解方法，由颜色分解装置执行，所述颜色分解 方法具备：目标值取得工序，在根据连续色调图像数据而生成在印刷装置中表示构成图像 的网点的形状及配置的各色材的二值图像数据时，取得在印刷装置中重现的颜色的目标 值；网点阈值数据取得工序，取得包含用于将连续色调图像数据变换成各色材的二值图像 数据的各网点的阈值的信息的网点阈值数据；印刷轮廓取得工序，取得表示印刷装置中的 设备信号值与表色系的值之间的对应关系的印刷轮廓；及颜色分解工序，基于由目标值取 得工序取得的颜色的目标值和印刷轮廓，在印刷装置中算出设备信号值的备选，基于网点 阈值数据及印刷轮廓，从设备信号值的备选之中决定用于使目标值所对应的颜色重现的设 备信号值。
[0019] 本发明的第九方案的颜色分解程序使计算机实现如下功能：目标值取得功能，在 根据连续色调图像数据而生成在印刷装置中表示构成图像的网点的形状及配置的各色材 的二值图像数据时，取得在印刷装置中重现的颜色的目标值；网点阈值数据取得功能，取得 包含用于将连续色调图像数据变换成各色材的二值图像数据的各网点的阈值的信息的网 点阈值数据；印刷轮廓取得功能，取得表示印刷装置中的设备信号值与表色系的值之间的 对应关系的印刷轮廓；及颜色分解功能，基于由目标值取得功能取得的颜色的目标值和印 刷轮廓，在印刷装置中算出设备信号值的备选，基于网点阈值数据及印刷轮廓，从设备信号 值的备选之中决定用于使目标值所对应的颜色重现的设备信号值。
[0020] 发明效果
[0021] 根据本发明，除了表示色差的大小的色差评价值之外，还使用画质评价值进行颜 色分解，该画质评价值使用网点重叠连续调明度图像的功率谱和视觉特性的功率谱而算 出，因此能够抑制与特殊颜色之间产生的色差，并抑制莫尔条纹、罗塞塔图案那样的图像结 构的产生。而且，通过进行网点重叠模拟，无需进行各设备信号值的块的印刷及测定，通过 更简易的运算，能够实现可抑制色差及抑制图像结构的产生的最佳的颜色分解。
【专利附图】
【附图说明】
[0022] 图1是表示本发明的一实施方式的图像形成系统的结构的框图。
[0023] 图2是表示本发明的一实施方式的图像处理装置的结构的框图。
[0024] 图3是用于说明网点重叠模拟的基本处理的图。
[0025] 图4是用于说明根据设备信号值（CMYK)的备选来生成单色的网点图像的处理的 图。
[0026] 图5是表示进行网点形状变化的推定处理时的网点重叠模拟处理部的结构的框 图。
[0027] 图6是表示版上网点形状推定处理部的结构的框图。
[0028] 图7是表示印刷物上网点形状推定处理部的结构的框图。
[0029] 图8是用于说明系统响应特性的推定处理的框图。
[0030] 图9(a)是表不响应特性取得用网点图像数据的俯视图，图9 (b)是表不印刷版 (版实物）中的网点的形状的俯视图，图9(c)是表示印刷物（印刷物实物）的网点的形状 的俯视图。
[0031] 图10(a)至图10(f)是表示曝光响应特性的例子的坐标图。
[0032] 图11是表示曝光响应特性推定部的结构的框图。
[0033] 图12是表示印刷响应特性推定部的结构的框图。
[0034] 图13是表示画质评价值计算部的运算处理的框图。
[0035] 图14是表示将观察距离（L)设为300 (mm)时的视觉特性（Dooley-Shaw函数）的 坐标图。
[0036] 图15(a)至图15(d)是表示网点重叠连续调明度图像的功率谱的例子的坐标图。
[0037] 图16是表示本发明的一实施方式的颜色分解处理的流程图。
[0038] 图17是表示用于选择颜色分解结果的⑶I的图。
[0039] 附图标记说明
[0040] 10…图像形成系统，100…图像处理装置，102…特殊颜色数据库（DB)，104…测色 部，106…设备信号值备选计算部，108…印刷轮廓存储部，110…颜色评价值保持部，112"-网点阈值数据存储部，114…网点重叠模拟部，116…画质评价值计算部，118…综合评价值 计算部，120…最佳值决定部，200…制版装置，300…印刷装置
【具体实施方式】
[0041] 以下，按照附图，说明本发明的颜色分解装置、颜色分解方法及颜色分解程序的实 施方式。
[0042][图像形成系统的结构]
[0043]图1是表示本发明的一实施方式的图像形成系统的结构的框图。
[0044] 如图1所示，本实施方式的图像形成系统10具备图像处理装置100、制版装置200 及印刷装置300。
[0045] 印刷用图像数据D10包含连续色调的图像信号。图像处理装置100对输入的印刷 用图像数据D10进行颜色分解处理，生成表示印刷装置300中使用的各色材（在本实施方 式中，为青色（C)、品红色（M)、黄色（Y)、黑色（K)这四种颜色）的网点图像的制版用数据 D20而向制版装置200输出。
[0046] 制版装置200是基于从图像处理装置100供给的制版用数据D20而生成具备与多 个网点对应的突起部的印刷版250的CTP(ComputerToPlate)描绘装置。制版装置200 具备激光雕刻机，在印刷版250的主面上雕刻与多个网点对应的突起部，由此生成印刷版 250。作为印刷版250的材料，在苯胺印刷的情况下，可使用例如橡胶片、光固化性树脂片等 高弹性材料。
[0047] 印刷装置300向由制版装置200生成的印刷版250赋予色材（油墨），使油墨向记 录介质（记录纸）的版面转移而形成印刷物400。在苯胺印刷的情况下，作为记录纸，除了 印刷用途的单张纸、卷筒纸之外，还使用表面的凹凸形状大的印刷介质（例如，瓦楞纸、薄 膜、布）等。印刷装置300与制版装置200 -起构成向记录纸的版面形成图像的图像形成 装直。
[0048][图像处理装置的结构]
[0049] 图2是表示本发明的一实施方式的图像处理装置（颜色分解装置）的结构的框 图。
[0050] 图2所示的图像处理装置100取得在印刷装置300中重现的颜色的目标值。并 且，图像处理装置100进行颜色分解处理，基于上述取得的颜色的目标值和印刷轮廓，取得 在印刷装置300中最佳的设备信号值。
[0051] 在图2中，特殊颜色数据库（DB) 102及测色部104构成目标值取得单元，网点阈值 数据存储部112及网点重叠模拟部114构成网点阈值数据取得单元，设备信号值备选计算 部106及印刷轮廓存储部108构成印刷轮廓取得单元，设备信号值备选计算部106、颜色评 价值保持部110、网点重叠模拟部114、画质评价值计算部116、综合评价值计算部118及最 佳值决定部120构成颜色分解单元。
[0052] 另外，在本实施方式中，说明使用L*a*b*表色系（JIS(JapaneseIndustrial StandardsCommittee)Z8729)作为表示颜色的目标值的表色系及印刷轮廓的表色系的例 子，但本发明的表色系没有限定于此。例如，除了可以使用国际照明委员会（CIE)规定的XYZ表色系（包含亮度（明亮度）的刺激值Y、颜色的刺激值X、Z)、Yxy表色系（亮度Y、 色度坐标x、y)、表色系（JISZ8518)之外，还可以使用HSV表色系（色相H(hue)、 彩度S(saturation)、明度V(value)或B(brightness))、HLS表色系（色相H(hue)、彩度 S(saturation)、亮度L(luminance))、YCbCr表色系（亮度Y、色差Cb、Cr)。
[0053]〈特殊颜色的取得〉
[0054] 在特殊颜色数据库（特殊颜色DB) 102中存储有表示PANT0NE、DIC、T0Y0等油墨厂 家的颜色样本中的特殊颜色的编号（标识符、特殊颜色ID)与特殊颜色ID所对应的表色系 的值（在本实施方式中，为Lab值（明度值L*、色度值a*、b*))的对应关系的一览表（LUT)。 图像处理装置100从使用者接受用于确定在印刷物400上重现的成为目标的颜色的特殊颜 色ID的输入，从特殊颜色DB102取得输入的特殊颜色ID所对应的表色系的值（Lab值）。
[0055] 测色部104具备测色机（例如，分光测色机等光学测色机）。测色部104使用测 色机来测定特殊颜色实物（例如，取色条、印刷物）或能够测定颜色的物品（成为颜色标本 的物品）的颜色所对应的表色系的值（Lab值），该特殊颜色实物涂敷有用于确定在印刷物 400上重现的目标颜色的特殊颜色。
[0056] 由上述特殊颜色ID确定出的特殊颜色的Lab值或由测色机测定出的特殊颜色的 Lab值（表示在印刷物400上重现的目标颜色的目标值）向设备信号值备选计算部106输 出。
[0057] 另外，在本实施方式中，说明颜色分解成CMYK这四种颜色的例子，但是印刷装置 300中的色材的数目及种类没有限定于此。即使在例如色材的数目为三种颜色（例如，CMY) 的情况下，也能够适用本发明的颜色分解方法。而且，即使在色材的数目为五种颜色以上 (例如，CMYK+追加的特殊颜色一种颜色（例如，使用者指定的颜色、或基于使用者指定的颜 色的色相而图像处理装置100所选择的颜色）这五种颜色、CMYK+追加的特殊颜色两种颜 色（例如，"橙色（0)、绿色（G)"、"橙色（0)、紫色（V)"）这六种颜色、CMYK+追加的特殊颜 色三种颜色（例如，"红色（R)、绿色（G)、蓝色（B)"、"0、G、V"）这七种颜色）的情况下，也 能够适用本发明的颜色分解方法。
[0058]〈设备信号值备选及颜色评价值的计算〉
[0059] 在印刷轮廓存储部108中存储有表示色材（例如，印刷色（CMYK))所对应的设 备信号值与表色系的值（在本实施方式中，为Lab值）的对应关系的一览表（LUT;用于 将设备信号值变换成表色系的值（Lab值）的印刷轮廓）。在此，作为印刷轮廓，可以使 用例如国际色彩委员会（InternationalColorConsortium:ICC)的轮廓的A2B1表。另 夕卜，作为颜色空间变换的方法，除了A2Bl(RelativeColormetric(测色一致））以外，还有 A2B0 (Perc印tual(感觉上的一致））、A2B2(Saturation(彩度增强）），但是在采用A2B1表 的情况下，能维持相对的色域，因此能够使色域（Gamut)变换（压缩）的影响为最小限度。
[0060] 设备信号值备选计算部106从印刷轮廓存储部108取得用于将色材（例如，印刷 色（CMYK))所对应的设备信号值变换成表色系的值（Lab值）的印刷轮廓。并且，设备信号 值备选计算部106使用该印刷轮廓进行逆运算，由此算出使在记录纸上重现的目标的Lab 值（目标Lab值）在容许色差的范围内能够重现的设备信号值（印刷装置300的CMYK值） 的备选（(cl、ml、yl、kl)、（c2、m2、y2、k2)、（c3、m3、y3、k3)、…）。
[0061] 具体而言，在将容许色差设为dEa、将目标Lab值设为（Lt、at、bt)且将向印刷轮 廓（LUT)输入了（C、M、Y、K)时的输出值（Lab值）设为（Lp、ap、bp) =LUT(C、M、Y、K)时， 设备信号值备选计算部106检索满足下述的[数学式1]的条件的CMYK值，求出CMYK值的 备选。
[0062][数学式1]
【权利要求】
1. 一种颜色分解装置，具备： 目标值取得单元，在根据连续色调图像数据而生成在印刷装置中表示构成图像的网点 的形状及配置的各色材的二值图像数据时，取得在所述印刷装置中重现的颜色的目标值； 网点阈值数据取得单元，取得包含用于将所述连续色调图像数据变换成各所述色材的 二值图像数据的各所述网点的阈值的信息的网点阈值数据； 印刷轮廓取得单元，取得表示所述印刷装置中的设备信号值与表色系的值之间的对应 关系的印刷轮廓；及 颜色分解单元，基于由所述目标值取得单元取得的所述颜色的目标值和所述印刷轮 廓，在所述印刷装置中算出所述设备信号值的备选，基于所述网点阈值数据及所述印刷轮 廓，从所述设备信号值的备选之中决定用于使所述目标值所对应的颜色重现的设备信号 值。
2. 根据权利要求1所述的颜色分解装置，其中，还具备： 颜色评价值保持单元，对应所述设备信号值的各备选，保持表示与所述目标值所对应 的颜色之间的色差的颜色的评价值；及 画质评价值计算单元，基于所述网点阈值数据及所述印刷轮廓，对应所述设备信号值 的各备选而算出画质的评价值， 所述颜色分解单元基于所述颜色的评价值及所述画质的评价值来决定用于使所述目 标值所对应的颜色重现的设备信号值。
3. 根据权利要求2所述的颜色分解装置，其中， 还具备模拟单元，该模拟单元使用所述网点阈值数据对所述设备信号值的备选实施二 值化处理，对于所述各备选生成各所述色材的二值图像数据，使各所述色材的二值图像数 据重叠，由此进行由所述印刷装置形成于印刷物上的图像的模拟， 所述画质评价值计算单元基于形成于所述印刷物上的图像的模拟的结果，对应各所述 备选而算出所述画质的评价值。
4. 根据权利要求3所述的颜色分解装置，其中， 所述模拟单元在所述印刷装置中，基于形成了所述二值图像数据所对应的网点时的响 应特性，考虑由所述印刷装置印刷所述网点时产生的所述网点的形状的变化来生成所述二 值图像数据。
5. 根据权利要求2所述的颜色分解装置，其中， 还具备油墨成本评价值计算单元，该油墨成本评价值计算单元对于所述设备信号值的 各备选，基于各色材的网点面积率来算出油墨成本的评价值， 所述颜色分解单元基于所述颜色的评价值、所述画质的评价值及所述油墨成本的评价 值来决定用于使所述目标值所对应的颜色重现的设备信号值。
6. 根据权利要求3所述的颜色分解装置，其中， 还具备油墨成本评价值计算单元，该油墨成本评价值计算单元对于所述设备信号值的 各备选，基于各色材的网点面积率来算出油墨成本的评价值， 所述颜色分解单元基于所述颜色的评价值、所述画质的评价值及所述油墨成本的评价 值来决定用于使所述目标值所对应的颜色重现的设备信号值。
7. 根据权利要求4所述的颜色分解装置，其中， 还具备油墨成本评价值计算单元，该油墨成本评价值计算单元对于所述设备信号值的 各备选，基于各色材的网点面积率来算出油墨成本的评价值， 所述颜色分解单元基于所述颜色的评价值、所述画质的评价值及所述油墨成本的评价 值来决定用于使所述目标值所对应的颜色重现的设备信号值。
8. 根据权利要求2所述的颜色分解装置，其中，还具备： 加权系数设定单元，对所述各评价值分别设定加权系数；及 综合评价值计算单元，通过使用由所述加权系数设定单元设定的加权系数进行所述各 评价值的加权相加，对应所述设备信号值的各备选而算出综合评价值， 所述颜色分解单元基于所述综合评价值来决定用于使所述目标值所对应的颜色重现 的设备信号值。
9. 根据权利要求3所述的颜色分解装置，其中，还具备： 加权系数设定单元，对所述各评价值分别设定加权系数；及 综合评价值计算单元，通过使用由所述加权系数设定单元设定的加权系数进行所述各 评价值的加权相加，对应所述设备信号值的各备选而算出综合评价值， 所述颜色分解单元基于所述综合评价值来决定用于使所述目标值所对应的颜色重现 的设备信号值。
10. 根据权利要求4所述的颜色分解装置，其中，还具备： 加权系数设定单元，对所述各评价值分别设定加权系数；及 综合评价值计算单元，通过使用由所述加权系数设定单元设定的加权系数进行所述各 评价值的加权相加，对应所述设备信号值的各备选而算出综合评价值， 所述颜色分解单元基于所述综合评价值来决定用于使所述目标值所对应的颜色重现 的设备信号值。
11. 根据权利要求5所述的颜色分解装置，其中，还具备： 加权系数设定单元，对所述各评价值分别设定加权系数；及 综合评价值计算单元，通过使用由所述加权系数设定单元设定的加权系数进行所述各 评价值的加权相加，对应所述设备信号值的各备选而算出综合评价值， 所述颜色分解单元基于所述综合评价值来决定用于使所述目标值所对应的颜色重现 的设备信号值。
12. 根据权利要求6所述的颜色分解装置，其中，还具备： 加权系数设定单元，对所述各评价值分别设定加权系数；及 综合评价值计算单元，通过使用由所述加权系数设定单元设定的加权系数进行所述各 评价值的加权相加，对应所述设备信号值的各备选而算出综合评价值， 所述颜色分解单元基于所述综合评价值来决定用于使所述目标值所对应的颜色重现 的设备信号值。
13. 根据权利要求7所述的颜色分解装置，其中，还具备： 加权系数设定单元，对所述各评价值分别设定加权系数；及 综合评价值计算单元，通过使用由所述加权系数设定单元设定的加权系数进行所述各 评价值的加权相加，对应所述设备信号值的各备选而算出综合评价值， 所述颜色分解单元基于所述综合评价值来决定用于使所述目标值所对应的颜色重现 的设备信号值。
14. 根据权利要求1?13中任一项所述的颜色分解装置，其中， 所述目标值取得单元通过接受用于确定颜色的标识符的输入或测定颜色，取得在所述 印刷装置中重现的所述颜色的目标值。
15. -种颜色分解方法，由颜色分解装置执行，所述颜色分解方法具备： 目标值取得工序，在根据连续色调图像数据而生成在印刷装置中表示构成图像的网点 的形状及配置的各色材的二值图像数据时，取得在所述印刷装置中重现的颜色的目标值； 网点阈值数据取得工序，取得包含用于将所述连续色调图像数据变换成各所述色材的 二值图像数据的各所述网点的阈值的信息的网点阈值数据； 印刷轮廓取得工序，取得表示所述印刷装置中的设备信号值与表色系的值之间的对应 关系的印刷轮廓；及 颜色分解工序，基于由所述目标值取得工序取得的所述颜色的目标值和所述印刷轮 廓，在所述印刷装置中算出所述设备信号值的备选，基于所述网点阈值数据及所述印刷轮 廓，从所述设备信号值的备选之中决定用于使所述目标值所对应的颜色重现的设备信号 值。
【文档编号】H04N1/46GK104427204SQ201410437695
【公开日】2015年3月18日 申请日期:2014年8月29日 优先权日:2013年8月30日
【发明者】堀田修平 申请人:富士胶片株式会社