WO2017156815A1 - Improved method for generating halftone projection and model for 3d printing - Google Patents

Improved method for generating halftone projection and model for 3d printing Download PDF

Info

Publication number
WO2017156815A1
WO2017156815A1 PCT/CN2016/079176 CN2016079176W WO2017156815A1 WO 2017156815 A1 WO2017156815 A1 WO 2017156815A1 CN 2016079176 W CN2016079176 W CN 2016079176W WO 2017156815 A1 WO2017156815 A1 WO 2017156815A1
Authority
WO
WIPO (PCT)
Prior art keywords
projection
printing
model
hole
circle
Prior art date
Application number
PCT/CN2016/079176
Other languages
French (fr)
Chinese (zh)
Inventor
陈宝权
科恩奥⋅丹尼尔
里奇斯⋅达尼
吕琳
沙夫⋅安德雷
魏源
赵海森
Original Assignee
山东大学
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 山东大学 filed Critical 山东大学
Publication of WO2017156815A1 publication Critical patent/WO2017156815A1/en

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/30—Polynomial surface description
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B29—WORKING OF PLASTICS; WORKING OF SUBSTANCES IN A PLASTIC STATE IN GENERAL
    • B29C—SHAPING OR JOINING OF PLASTICS; SHAPING OF MATERIAL IN A PLASTIC STATE, NOT OTHERWISE PROVIDED FOR; AFTER-TREATMENT OF THE SHAPED PRODUCTS, e.g. REPAIRING
    • B29C67/00—Shaping techniques not covered by groups B29C39/00 - B29C65/00, B29C70/00 or B29C73/00
    • B29C67/0007—Manufacturing coloured articles not otherwise provided for, e.g. by colour change
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B33—ADDITIVE MANUFACTURING TECHNOLOGY
    • B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
    • B33Y10/00—Processes of additive manufacturing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G06T15/50—Lighting effects
    • G06T15/506—Illumination models
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T15/00—3D [Three Dimensional] image rendering
    • G06T15/50—Lighting effects
    • G06T15/55—Radiosity

Definitions

  • the present invention relates to an improved halftone projection and model generation method for 3D printing.
  • 3D Printing also known as Additive Manufacturing (AM), Rapid Prototyping, etc.
  • AM Additive Manufacturing
  • Rapid Prototyping etc.
  • 3D printing can be made of metal, photosensitive resin, plastic and other materials, directly using digital model files as input to create 3D solids of any complex shape, suitable for customizable manufacturing.
  • Continuously tuned image refers to a true color image in which the hue, brightness, and saturation continuously change on an image, and the continuous change is composed of the density of the image material particles per unit area, such as a CRT display.
  • the depth change of the continuous tone image is stepless.
  • the corresponding halftone image also known as the mesh image, exhibits a relatively small hue, and the image detail changes discontinuously, such as an inkjet type printer, by the size of the dots or the level of sparsely expressed images.
  • a method of generating a halftone image from a continuous tone image may be referred to as a halftone image generation technique, or a halftone technique.
  • Halftone technology has been widely used in traditional paper printing and digital display. At its core are problems such as tone reproduction, structure retention, dot density and image resolution matching. After decades of exploration and research, many halftone generation methods have emerged. In order to maintain the relative color of the original image, domestic and foreign research researchers have proposed a number of corresponding halftone techniques.
  • Patent No. CN201410420912.4 proposes a halftone projection and model generation method for 3D printing. By controlling the position, size and length of each hole, a macroscopically complete projection image is obtained, and feature layers and tones are introduced. The layers respectively generate corresponding holes, and finally the two layers are merged and printed to generate a model.
  • its projection effect may be discontinuous with discrete light halves.
  • the model Compared with this patent, the model generates calculations using completely different techniques, and can obtain more continuous grayscale projection images without visible visible light spots.
  • the final projection effect is significantly better than the patent CN201410420912.4.
  • the technology proposed in this patent is directed to the field of 3D printing, in which a spot formed by direct projection of light is used as a basic unit constituting a halftone image, and a continuous projection image having a progressive gradation change is projected in a three-dimensional space.
  • Improved existing adoption The above problem is solved by optimizing the geometry of the object to obtain various related techniques of light and shadow effects, and using a technique of refracting light to obtain a specified projected image.
  • the present invention proposes an improved halftone projection and model generation method for 3D printing.
  • the method sets a surface on the surface of the three-dimensional solid model according to any grayscale image and three-dimensional solid model given by the user.
  • An improved halftone projection and model generation method for 3D printing comprising the following steps:
  • the specific method of the projection simulation method includes:
  • the correction method in the step (1-4) is corrected by Gamma.
  • a specific method for calculating a projected reference image includes:
  • satisfying the printability condition means that the radius of the hole in the three-dimensional model is not less than the minimum radius of the hole that can be printed in the selected printing technique, and the distance between the holes is not less than the Printing technology prints the minimum distance between two holes.
  • the density map is calculated, and specifically includes the following steps:
  • the shading rate of the corresponding three-dimensional model is set according to the irradiance of the illuminance, and the corresponding target circle radius is calculated according to the shading rate;
  • the target circle radius corresponding to each point of the illuminance map is mapped to the corresponding density value, and normalized to obtain a density map.
  • the specific method for calculating the optimal circle arrangement includes the following steps:
  • (3-1) calculating the optimal target circle number according to the accumulated density value of the density map, the number of circles of the projected reference image, and the accumulated density value;
  • the specific method includes:
  • the desired hole type includes an enlarged hole and a slanted hole.
  • the method is: nesting an inscribed circle that reduces a safe distance in a maximum inscribed circle of the Voronoi region.
  • the inscribed circle is projected on the outer surface and the inner surface of the three-dimensional shell model by central projection, respectively forming two intersecting ellipse, and the enlarged hole is formed by connecting an ellipse of the inner and outer surfaces with a cylindrical structure.
  • the method is: inserting an inscribed circle that reduces a safe distance in a maximum inscribed circle of the Voronoi region. Selecting a random direction to place two minimum circles D 1 and D 2 satisfying the printable condition in the inscribed circle, respectively projecting D 1 and D 2 on the inner surface of the three-dimensional shell model by central projection and The outer surface forms two intersecting ellipse, respectively, and the elliptical shape connecting the inner and outer surfaces is formed by a cylindrical structure to form the inclined hole.
  • the present invention proposes a halftone projection and model generation method for extending a medium in a halftone technique from digital printing ink to projection light;
  • the model generated by the present invention causes the light source to form a continuous grayscale image closest to a given grayscale image on the projection receiving surface through the hole in the model, and no visible discrete spot is visible;
  • the invention directly generates a physical model satisfying 3D printing constraints for 3D printing, and supports the user to customize any model shape and target projection image;
  • Figure 1 is a flow chart of the present invention
  • Figure 3 (a) is an arrangement pattern of holes arranged to satisfy the most closely arranged printable conditions
  • Figure 3 (b) is a projection projection image generated by calling a projection simulation method
  • Figure 3 (c) is a corresponding projected image obtained by actual shooting
  • Figure 5 is a view showing an application of the present invention.
  • Fig. 6 is an analytical diagram showing the generation of enlarged holes and inclined holes.
  • a halftone projection and model generation method for 3D printing includes the following steps:
  • the user specifies the three-dimensional model, the projection receiving surface and the position and size range of the light source, and the projection reference image B 0 is calculated by the projection simulation method;
  • the projection simulation method specifically includes the following steps:
  • the total irradiance E v (p) of the projection receiving point p is the sum of the illuminances of all the point sources to the point.
  • r i is the Euclidean distance from the point p to the light source l i
  • ⁇ i and ⁇ p are the straight lines connecting the points p and l i and the normals at the points p and l i with
  • the angle of the intersection, V(p, l i ) is the visible relationship between points p and l i , the value of 0 means invisible, and the value of 1 means visible;
  • the projection reference image is calculated, and specifically includes the following steps:
  • the holes having the smallest radius of the printability condition are closely arranged on the three-dimensional model, as shown in Fig. 3(a). Satisfying the printability condition here means that the radius of the hole in the three-dimensional model is not less than r min , and the distance between the holes is not less than d min .
  • r min is the minimum radius at which a particular printing technique can print holes
  • d min is the minimum distance between two holes that the printing technique can print.
  • the holes that meet the printable conditions are arranged closely.
  • the holes with a radius of r min are arranged in the most closely arranged arrangement with the spacing d min , as shown in Fig. 3(a);
  • Fig. 3(b) The projection simulation method is called to generate a projection reference image B 0 as shown in Fig. 3(b).
  • Fig. 3(c) is a corresponding projected image obtained by actual shooting.
  • the projected reference image B 0 is the lowest gray value that the enlarged hole can reach. If the brightness value of the specific region of the input image I t is not lower than B 0 , it is necessary to arrange the enlarged holes; if the brightness value of the specific region below the input image I t is lower than B 0 , it is necessary to arrange the inclined holes.
  • the density map M is calculated, which specifically includes the following steps:
  • the illuminance map is obtained by inverse Gamma correction for a given grayscale image I t .
  • the corresponding irradiance is E v (p).
  • K can be expressed as a functional form of r. As shown in FIG 4, r of the inscribed circle radius of the largest regular hexagonal modification, r min is the minimum radius of the holes may be printed, d min is the minimum printable distance between the two holes.
  • the radius of the largest inscribed circle is expected to be r.
  • r is the radius of the maximum inscribed circle corresponding to the upper surface of the three-dimensional model, which receives the projection facing r is the radius of the circle to be W, the respective position of the projection-receiving surface of the circle of radius r W at point p (x, y)
  • the projected area at the corresponding three-dimensional model is equal to the area of the circle having the radius r;
  • the optimal circle arrangement is calculated, which specifically includes the following steps:
  • the CCVT is calculated by calling the method of de Goes. Calculating its maximum inscribed circle in each Voronoi region of the resulting CCVT results in a set of circular arrangements. For a certain Voronoi region, it corresponds to a certain block region of the density map M, and the mean value of the radius r w corresponding to all the points in the region is taken as the desired circle size of the region. It is thus determined whether the maximum inscribed circle of the region has reached its desired circle size. With the target circle number N as the upper bound, the binary search is used to find the CCVT that achieves the optimal circle correct rate. As shown in Fig. 5, 5(a) is the user input grayscale image, 5(b) is the corresponding density map M, 5(c) is the CCVT calculated by the method of de Goes, and 5(d) is the calculated optimal circle. arrangement.
  • step (4) determining the size, position and relative light source orientation of the corresponding holes, specifically including the following steps:
  • step 3 From step 3, the CCVT and its corresponding circle arrangement can be obtained to achieve the optimal correct rate.
  • Any point in the density map M records its corresponding density value ⁇ w (x, y) and its desired hole type, enlarged hole or inclined hole.
  • the hole type of the surface of most points in the region is taken as the type of the corresponding hole of the Voronoi region.
  • 3D printing can be directly performed according to the model.
  • the printed model allows the light source to be projected through the holes in the model on the projection receiving surface to form a continuous grayscale image that is closest to the given grayscale image of the user.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Graphics (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Analysis (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Algebra (AREA)
  • Manufacturing & Machinery (AREA)
  • Mechanical Engineering (AREA)
  • Chemical & Material Sciences (AREA)
  • Materials Engineering (AREA)

Abstract

An improved method for generating halftone projection and a model for 3D printing, the method comprising: calculating an reference projection image by using a project simulation method according to a given three-dimensional model, a projection receiving surface and basic parameters of a light source; correcting a given grayscale image to obtain an illuminance diagram; corresponding points in the illuminance diagram to the reference projection image; normalizing to obtain a density graph by means of calculation; using a Voronoi division method having capacity constraints to provide an optimal circle permutation according to the density graph; and according to the type of a hole as determined by the optimal circular permutation, determining the corresponding size and position of the hole as well as the direction and angle of the light source, and generating a physical model for 3D printing. By means of the present method, the light source transmits through the hole on the model to form a continuous grayscale image which is closest to the given grayscale image, and no obviously visible discrete spot is present.

Description

一种改进的面向3D打印的半色调投影与模型生成方法An improved halftone projection and model generation method for 3D printing 技术领域Technical field
本发明涉及一种改进的面向3D打印的半色调投影与模型生成方法。The present invention relates to an improved halftone projection and model generation method for 3D printing.
背景技术Background technique
3D打印(3D Printing)，又称增材制造(Additive Manufacturing,AM)，快速原型制造(Rapid Prototyping)等，是一种基于离散-堆积原理，采用材料逐层累加的方法制造实体零件的技术。3D打印可采用金属、光敏树脂、塑料等多种材料，直接以数字模型文件为输入制造任意复杂形状的三维实体，适用于可定制化的制造。3D Printing, also known as Additive Manufacturing (AM), Rapid Prototyping, etc., is a technique for manufacturing solid parts based on the principle of discrete-stacking. 3D printing can be made of metal, photosensitive resin, plastic and other materials, directly using digital model files as input to create 3D solids of any complex shape, suitable for customizable manufacturing.
连续调图像是指在一幅图像上，存在着色调、亮度与饱和度连续变化的真彩色图像，其连续变化是以单位面积成像物质颗粒的密度构成的，如CRT显示器。连续调图像的深浅变化是无级的。与之对应的半色调图像，又称为网目图像，表现的色调则相对少一些，通过网点的大小或稀疏表达图像的层次，图像细节的变化不连续，如喷墨类型的打印机。由连续调图像生成半色调图像的方法，可称之为半色调图像生成技术，或半色调技术。Continuously tuned image refers to a true color image in which the hue, brightness, and saturation continuously change on an image, and the continuous change is composed of the density of the image material particles per unit area, such as a CRT display. The depth change of the continuous tone image is stepless. The corresponding halftone image, also known as the mesh image, exhibits a relatively small hue, and the image detail changes discontinuously, such as an inkjet type printer, by the size of the dots or the level of sparsely expressed images. A method of generating a halftone image from a continuous tone image may be referred to as a halftone image generation technique, or a halftone technique.
半色调技术已经广泛应用于传统的纸面印刷和数字显示等领域。其核心在于色调再现，结构保持，点密度和图像分辨率的匹配等问题。经过几十年的探索研究，出现了很多半色调生成方法。以保持原始图像的相对色调为目的，国内外的研究学者们提出了很多相应的半色调技术。Halftone technology has been widely used in traditional paper printing and digital display. At its core are problems such as tone reproduction, structure retention, dot density and image resolution matching. After decades of exploration and research, many halftone generation methods have emerged. In order to maintain the relative color of the original image, domestic and foreign research scholars have proposed a number of corresponding halftone techniques.
但是，已有的半色调图像生成的技术，所面向的是数字半色调图像生成或者给定图像的点刻画表达。However, existing techniques for generating halftone images are directed to digital halftone image generation or point characterization of a given image.
申请号为CN201410420912.4的专利，提出一种面向3D打印的半色调投影与模型生成方法，通过控制每个孔的位置，大小和长度，得到宏观完整的投影图像，并引入了特征层和色调层分别生成对应孔洞，最终将两层融合并打印生成模型。但是，其投影效果可能存在不连续、有离散光半点。Patent No. CN201410420912.4 proposes a halftone projection and model generation method for 3D printing. By controlling the position, size and length of each hole, a macroscopically complete projection image is obtained, and feature layers and tones are introduced. The layers respectively generate corresponding holes, and finally the two layers are merged and printed to generate a model. However, its projection effect may be discontinuous with discrete light halves.
本专利与之相比采用完全不同的技术进行模型生成计算，并且能得到更连续的灰度投影图像，无明显可见的离散光斑点，最终投影效果显著优于专利CN201410420912.4。Compared with this patent, the model generates calculations using completely different techniques, and can obtain more continuous grayscale projection images without visible visible light spots. The final projection effect is significantly better than the patent CN201410420912.4.
本专利提出的技术面向3D打印领域，以光通过直射投影所形成的光斑作为组成半色调图像的基本单元，在三维空间中投影出有渐进灰度变化的连续投影图像。改进了现有的通过 优化调整物体的几何结构得到各种不同的光影效果的相关技术、利用光线折射得到指定的投影图像的技术，解决了上述问题。The technology proposed in this patent is directed to the field of 3D printing, in which a spot formed by direct projection of light is used as a basic unit constituting a halftone image, and a continuous projection image having a progressive gradation change is projected in a three-dimensional space. Improved existing adoption The above problem is solved by optimizing the geometry of the object to obtain various related techniques of light and shadow effects, and using a technique of refracting light to obtain a specified projected image.
发明内容Summary of the invention
本发明为了解决上述问题，提出了一种改进的面向3D打印的半色调投影与模型生成方法，本方法根据用户给定的任一灰度图片和三维实体模型，通过在三维实体模型表面设置一组不同大小，位置和相对光源朝向角度的孔洞，分别使用扩大型孔洞和倾斜型孔洞投影表示灰度图像中较高亮度和较低亮度的区域，使光源透过这些孔洞在投影接收面上形成一幅与给定灰度图像最相近的连续灰度图像。In order to solve the above problems, the present invention proposes an improved halftone projection and model generation method for 3D printing. The method sets a surface on the surface of the three-dimensional solid model according to any grayscale image and three-dimensional solid model given by the user. Grouping holes of different sizes, positions and angles with respect to the light source, respectively, using enlarged holes and oblique holes to respectively represent areas of higher brightness and lower brightness in the gray image, so that the light source is formed through the holes on the projection receiving surface A continuous grayscale image that is closest to a given grayscale image.
为了实现上述目的，本发明采用如下技术方案：In order to achieve the above object, the present invention adopts the following technical solutions:
一种改进的面向3D打印的半色调投影与模型生成方法，包括以下步骤：An improved halftone projection and model generation method for 3D printing, comprising the following steps:
(1)根据给定的三维模型、投影接收面以及光源的基本参数，利用投影模拟方法，计算投影参考图像；(1) calculating a projection reference image by using a projection simulation method according to a given three-dimensional model, a projection receiving surface, and basic parameters of the light source;
(2)将给定的灰度图像进行校正，得到照度图，将照度图的点与投影参考图像相对应，进行归一化处理后，计算得到密度图；(2) correcting the given grayscale image to obtain an illuminance map, corresponding to the projected reference image, and performing normalization processing to calculate the density map;
(3)利用带容量约束的Voronoi划分的方法，根据密度图进行最优圆排列；(3) Using the Voronoi partitioning method with capacity constraints, the optimal circle arrangement is performed according to the density map;
(4)根据最优圆排列确定的孔洞类型，确定相应孔洞的大小、位置和相对光源朝向角度，生成物理模型，进行3D打印。(4) According to the hole type determined by the optimal circle arrangement, determine the size, position and relative light source orientation angle of the corresponding hole, generate a physical model, and perform 3D printing.
所述步骤(1)中，投影模拟方法的具体方法包括：In the step (1), the specific method of the projection simulation method includes:
(1-1)将光源离散化为若干个点光源；(1-1) discretizing the light source into a plurality of point sources;
(1-2)在投影接收面的投影区域离散采样为若干个的投影接收点；(1-2) discretely sampling the projection area of the projection receiving surface into a plurality of projection receiving points;
(1-3)在多孔模型的遮挡作用下，计算投影接收点的所有点光源的总辐射照度；(1-3) Calculating the total irradiance of all point sources of the projection receiving point under the occlusion of the porous model;
(1-4)对投影接收点的总辐射照度进行校正，得到投影模拟图像灰度值。(1-4) Correcting the total irradiance of the projection receiving point to obtain the gradation value of the projected analog image.
进一步的，所述步骤(1-4)中校正方法为通过Gamma校正。Further, the correction method in the step (1-4) is corrected by Gamma.
所述步骤(1)中，计算投影参考图像的具体方法包括：In the step (1), a specific method for calculating a projected reference image includes:
(1-a)在三维模型上紧密排列满足可打印性条件半径最小的孔洞；(1-a) closely aligning the holes satisfying the minimum radius of printability conditions on the three-dimensional model;
(1-b)使用投影模拟方法，在孔洞最紧密排列的约束条件下，生成投影参考图像。(1-b) Using the projection simulation method, a projection reference image is generated under the constraint that the holes are most closely arranged.
所述步骤(1-a)中，满足可打印性条件指的是三维模型上的孔洞的半径不小于选择的打印技术中能够打印的孔洞的最小半径，同时，孔洞之间的距离不小于该打印技术可打印两孔洞间的最小距离。 In the step (1-a), satisfying the printability condition means that the radius of the hole in the three-dimensional model is not less than the minimum radius of the hole that can be printed in the selected printing technique, and the distance between the holes is not less than the Printing technology prints the minimum distance between two holes.
所述步骤(2)中，对于输入的任一灰度图像，计算密度图，具体包括以下步骤：In the step (2), for any grayscale image input, the density map is calculated, and specifically includes the following steps:
(2-1)将给定的灰度图像，通过逆校正得到照度图；(2-1) obtaining a illuminance map by inverse correction of a given grayscale image;
(2-2)对于照度图的每一点，根据其辐射照度设定对应的三维模型处遮光率，根据遮光率计算其对应的目标圆半径；(2-2) For each point of the illuminance map, the shading rate of the corresponding three-dimensional model is set according to the irradiance of the illuminance, and the corresponding target circle radius is calculated according to the shading rate;
(2-3)将照度图每点对应的目标圆半径映射为相应密度值，进行归一化处理，得到密度图。(2-3) The target circle radius corresponding to each point of the illuminance map is mapped to the corresponding density value, and normalized to obtain a density map.
所述步骤(3)中，计算最优圆排列的具体方法包括以下步骤：In the step (3), the specific method for calculating the optimal circle arrangement includes the following steps:
(3-1)根据密度图的累加密度值、投影参考图像的圆个数和累加密度值，计算最优目标圆个数；(3-1) calculating the optimal target circle number according to the accumulated density value of the density map, the number of circles of the projected reference image, and the accumulated density value;
(3-2)根据密度图和目标圆个数，通过二分搜索查找达到最优圆正确率的带容量约束的Voronoi划分。(3-2) According to the density map and the number of target circles, the Voronoi partition with capacity constraint that achieves the optimal circle correct rate is searched by binary search.
所述步骤(4)中，具体方法包括：In the step (4), the specific method includes:
(4-1)根据得到的最优正确率的Voronoi区域和其对应圆排列，将其对应的密度图的区域中多数点标明的孔洞类型作为该Voronoi区域期望的孔洞类型；(4-1) according to the obtained optimal correct rate of the Voronoi region and its corresponding circle arrangement, the hole type indicated by the majority of the points in the corresponding density map region is taken as the desired hole type of the Voronoi region;
(4-2)根据任意Voronoi区域期望的孔洞类型，在其三维模型的对应位置生成相应的孔洞。(4-2) According to the type of hole desired in any Voronoi region, a corresponding hole is generated at a corresponding position of the three-dimensional model.
所述步骤(4-1)中，期望的孔洞类型包括扩大型孔洞和倾斜型孔洞。In the step (4-1), the desired hole type includes an enlarged hole and a slanted hole.
进一步的，所述步骤(4-2)中，若某Voronoi区域对应的孔洞类型为扩大型孔洞，生成方法为：在该Voronoi区域的最大内切圆嵌套一个缩小一个安全距离的内切圆，将该内切圆通过中心投影的方式投影在三维壳状模型的外表面和内表面，分别形成两个相交的椭圆，使用一个圆柱形结构连接内外表面的椭圆形成该扩大型孔洞。Further, in the step (4-2), if the hole type corresponding to a Voronoi region is an enlarged hole, the method is: nesting an inscribed circle that reduces a safe distance in a maximum inscribed circle of the Voronoi region. The inscribed circle is projected on the outer surface and the inner surface of the three-dimensional shell model by central projection, respectively forming two intersecting ellipse, and the enlarged hole is formed by connecting an ellipse of the inner and outer surfaces with a cylindrical structure.
进一步的，所述步骤(4-2)中，若某Voronoi区域对应的孔洞类型为倾斜型孔洞，生成方法为：在该Voronoi区域的最大内切圆嵌套一个缩小一个安全距离的内切圆，在该内切圆内选一个随机方向放置两个半径满足可打印条件的最小圆D1和D2，分别将D1和D2通过中心投影的方式投影在三维壳状模型的内表面和外表面，分别形成两个相交的椭圆，使用一个圆柱形结构连接内外表面的椭圆形成该倾斜型孔洞。Further, in the step (4-2), if a hole type corresponding to a Voronoi region is an inclined hole, the method is: inserting an inscribed circle that reduces a safe distance in a maximum inscribed circle of the Voronoi region. Selecting a random direction to place two minimum circles D 1 and D 2 satisfying the printable condition in the inscribed circle, respectively projecting D 1 and D 2 on the inner surface of the three-dimensional shell model by central projection and The outer surface forms two intersecting ellipse, respectively, and the elliptical shape connecting the inner and outer surfaces is formed by a cylindrical structure to form the inclined hole.
本发明的有益效果为： The beneficial effects of the invention are:
(1)本发明提出半色调投影与模型生成方法，将半色调技术中的介质由数字印刷墨水推广到投影光线；(1) The present invention proposes a halftone projection and model generation method for extending a medium in a halftone technique from digital printing ink to projection light;
(2)本发明生成的模型，使光源透过模型上的孔洞在投影接收面上形成一幅与给定灰度图像最相近的连续灰度图像，无明显的离散光斑可见；(2) The model generated by the present invention causes the light source to form a continuous grayscale image closest to a given grayscale image on the projection receiving surface through the hole in the model, and no visible discrete spot is visible;
(3)本发明直接面向3D打印生成满足3D打印约束的物理模型，支持用户个性化定制任意的模型外形和目标投影图像；(3) The invention directly generates a physical model satisfying 3D printing constraints for 3D printing, and supports the user to customize any model shape and target projection image;
(4)应用范围广泛，适用于灯光艺术造型，灯具定制等多种场合。(4) A wide range of applications, suitable for a variety of occasions such as lighting art modeling, lighting customization.
附图说明DRAWINGS
图1为本发明的流程图；Figure 1 is a flow chart of the present invention;
图2为投影模拟方法解析图；2 is an analytical diagram of a projection simulation method;
图3(a)为满足可打印条件最紧密排列的孔洞的排布方式排列图像；Figure 3 (a) is an arrangement pattern of holes arranged to satisfy the most closely arranged printable conditions;
图3(b)为调用投影模拟方法生成投影参考图像；Figure 3 (b) is a projection projection image generated by calling a projection simulation method;
图3(c)为对应的实际拍摄得到的投影图像；Figure 3 (c) is a corresponding projected image obtained by actual shooting;
图4为扩大型孔洞和倾斜型孔洞遮光率解析图；4 is an analytical diagram of the aperture ratio of the enlarged hole and the inclined hole;
图5为本发明的适用实例图；Figure 5 is a view showing an application of the present invention;
图6为扩大型孔洞和倾斜型孔洞生成解析图。Fig. 6 is an analytical diagram showing the generation of enlarged holes and inclined holes.
具体实施方式：detailed description:
下面结合附图与实施例对本发明作进一步说明。The invention will be further described below in conjunction with the drawings and embodiments.
如图1所示，一种面向3D打印的半色调投影与模型生成方法，包括以下步骤：As shown in FIG. 1, a halftone projection and model generation method for 3D printing includes the following steps:
(1)用户指定三维模型，投影接收面及光源的位置大小范围等参数，通过投影模拟方法，计算投影参考图像B0；(1) The user specifies the three-dimensional model, the projection receiving surface and the position and size range of the light source, and the projection reference image B 0 is calculated by the projection simulation method;
(2)对于输入的任一灰度图像It，结合投影参考图像B0，计算一幅密度图M；(2) For any grayscale image I t input, combined with the projected reference image B 0 , calculate a density map M;
(3)根据密度图M，应用带容量约束的Voronoi划分(CCVT)的方法计算最优圆排列；(3) According to the density map M, the method of calculating the optimal circle by applying the Voronoi partition (CCVT) with capacity constraint;
(4)根据这组圆排列，生成相应的扩大型孔洞和倾斜型孔洞。(4) According to the set of circular arrays, corresponding enlarged holes and inclined holes are generated.
所述步骤(1)中，投影模拟方法，具体包括以下步骤：In the step (1), the projection simulation method specifically includes the following steps:
(1-1)将光通量为Φ的光源L离散化为n个点光源
Figure PCTCN2016079176-appb-000001
每个点光源li的光通量为
Figure PCTCN2016079176-appb-000002
(1-1) Discretizing the light source L having the luminous flux Φ into n point light sources
Figure PCTCN2016079176-appb-000001
The luminous flux of each point source l i is
Figure PCTCN2016079176-appb-000002
(1-2)在投影接收面的投影区域离散采样有限个数的投影接收点，如图2中的点p； (1-2) discretely sampling a limited number of projection receiving points in the projection area of the projection receiving surface, as shown by point p in FIG. 2;
(1-3)如图2，在多孔模型(lampshade)的遮挡作用下，计算投影接收点p的总辐射照度Ev(p)为:(1-3) As shown in Fig. 2, under the occlusion of the porous model, the total irradiance E v (p) of the projection receiving point p is calculated as:
Figure PCTCN2016079176-appb-000003
Figure PCTCN2016079176-appb-000003
投影接收点p的总辐射照度Ev(p)为所有点光源到该点的辐射照度累加和。如图2中，ri为点p到光源li的欧式距离，θi和θp为连接点p和li的直线和点p和li处法线
Figure PCTCN2016079176-appb-000004
和
Figure PCTCN2016079176-appb-000005
的夹角，V(p,li)为点p和li的可见关系，取值为0代表不可见，取值为1代表可见；The total irradiance E v (p) of the projection receiving point p is the sum of the illuminances of all the point sources to the point. As shown in Fig. 2, r i is the Euclidean distance from the point p to the light source l i , and θ i and θ p are the straight lines connecting the points p and l i and the normals at the points p and l i
Figure PCTCN2016079176-appb-000004
with
Figure PCTCN2016079176-appb-000005
The angle of the intersection, V(p, l i ) is the visible relationship between points p and l i , the value of 0 means invisible, and the value of 1 means visible;
(1-4)将投影接收点的总辐射照度Ev(p)，通过Gamma校正得到投影模拟图像灰度值(1-4) The total irradiance E v (p) of the projection receiving point is corrected by Gamma to obtain the gray value of the projected analog image.
It(p)＝g(Ev(p))＝(Ev(p))1/γ，其中g(.)表示Gamma校正过程，通常Gamma函数γ取值为2.2。I t (p)=g(E v (p))=(E v (p)) 1/γ , where g(.) represents the gamma correction process, and usually the Gamma function γ takes a value of 2.2.
所述步骤(1)中，计算投影参考图像，具体包括以下步骤：In the step (1), the projection reference image is calculated, and specifically includes the following steps:
(1-1)在三维模型上紧密排列满足可打印性条件半径最小的孔洞，如图3(a)。在此满足可打印性条件指的是三维模型上的孔洞的半径不小于rmin，孔洞之间的距离不小于dmin。其中，rmin为某特定打印技术可打印孔洞的最小半径，dmin为该打印技术可打印两孔洞间的最小距离。满足可打印条件最紧密排列的孔洞为一组半径为rmin的孔洞按照间距为dmin的最紧密的排布方式排列，如图3(a)所示；(1-1) The holes having the smallest radius of the printability condition are closely arranged on the three-dimensional model, as shown in Fig. 3(a). Satisfying the printability condition here means that the radius of the hole in the three-dimensional model is not less than r min , and the distance between the holes is not less than d min . Where r min is the minimum radius at which a particular printing technique can print holes, and d min is the minimum distance between two holes that the printing technique can print. The holes that meet the printable conditions are arranged closely. The holes with a radius of r min are arranged in the most closely arranged arrangement with the spacing d min , as shown in Fig. 3(a);
(1-2)调用投影模拟方法生成投影参考图像B0，如图3(b)。图3(c)为其对应的实际拍摄得到的投影图像。在孔洞最紧密排列的约束条件下，投影参考图像B0为扩大型孔洞能到达的最低灰度值。若输入图像It特定区域的亮度值不低于B0，需要排列扩大型孔洞；若低于输入图像It特定区域的亮度值低于B0，则需要排列倾斜型孔洞。(1-2) The projection simulation method is called to generate a projection reference image B 0 as shown in Fig. 3(b). Fig. 3(c) is a corresponding projected image obtained by actual shooting. Under the constraint that the holes are most closely arranged, the projected reference image B 0 is the lowest gray value that the enlarged hole can reach. If the brightness value of the specific region of the input image I t is not lower than B 0 , it is necessary to arrange the enlarged holes; if the brightness value of the specific region below the input image I t is lower than B 0 , it is necessary to arrange the inclined holes.
所述步骤(2)中，对于输入的任一灰度图像It，计算密度图M，具体包括以下步骤：In the step (2), for any grayscale image I t input, the density map M is calculated, which specifically includes the following steps:
(2-1)将给定的灰度图像It，通过逆Gamma校正得到照度图。对于灰度图像It上某一点p(x,y)的灰度It(p(x,y))，辐射照度Ev(p)＝g-1(It(p))，其中g-1(.)为逆Gamma校正过程。(2-1) The illuminance map is obtained by inverse Gamma correction for a given grayscale image I t . For grayscale images I t on a point p (x, y) is the gray I t (p (x, y )), irradiance E v (p) = g -1 (I t (p)), where g -1 (.) is the inverse gamma correction process.
(2-2)对于照度图的每一点p(x,y)，其对应辐射照度为Ev(p)。为达到辐射照度Ev(p)，在点p(x,y)对应三维模型处遮光率设为K，有
Figure PCTCN2016079176-appb-000006
其中，
Figure PCTCN2016079176-appb-000007
为点p(x,y)在没有任何三维模型遮挡的情况下的总辐射照度
Figure PCTCN2016079176-appb-000008
K表示点p(x,y)对应三维模型处对应的遮光率，K＝Area(unoccluded)/Area(Cell)。可将K表示为r的 函数形式。如图4，r为正六变形内最大内切圆的半径，rmin为可打印孔洞的最小半径，dmin为可打印两孔洞间的最小距离。(2-2) For each point p(x, y) of the illuminance map, the corresponding irradiance is E v (p). In order to achieve the irradiance E v (p), the shading rate is set to K at the point p(x, y) corresponding to the three-dimensional model,
Figure PCTCN2016079176-appb-000006
among them,
Figure PCTCN2016079176-appb-000007
Total irradiance for point p(x, y) without any 3D model occlusion
Figure PCTCN2016079176-appb-000008
K represents the point p(x, y) corresponding to the corresponding shading rate at the three-dimensional model, K = Area (unoccluded) / Area (Cell). K can be expressed as a functional form of r. As shown in FIG 4, r of the inscribed circle radius of the largest regular hexagonal modification, r min is the minimum radius of the holes may be printed, d min is the minimum printable distance between the two holes.
若It(p(x,y))≥B0(p(x,y))，需要排列扩大型孔洞，有：If I t (p(x, y)) ≥ B 0 (p(x, y)), it is necessary to arrange the enlarged holes, which are:
Figure PCTCN2016079176-appb-000009
Figure PCTCN2016079176-appb-000009
若It(p(x,y))<B0(p(x,y))，需要排列倾斜型孔洞，有：If I t (p(x, y)) < B 0 (p(x, y)), it is necessary to arrange the inclined holes, which are:
Figure PCTCN2016079176-appb-000010
其中d＝r-rmin-0.5dmin.
Figure PCTCN2016079176-appb-000010
Where d = rr min -0.5d min .
由上述公式可得，对于照度图上点p(x,y)，为达到Ev(p)，期望最大内切圆的半径为r。r为三维模型表面上相应最大内切圆的半径，其在投影接收面对应圆的半径为rw，使投影接收面上相应位置上半径为rw的圆在点p(x,y)对应三维模型处的投影面积与半径为r的圆的面积相等；It can be obtained from the above formula that for the point p(x, y) on the illuminance map, to reach E v (p), the radius of the largest inscribed circle is expected to be r. r is the radius of the maximum inscribed circle corresponding to the upper surface of the three-dimensional model, which receives the projection facing r is the radius of the circle to be W, the respective position of the projection-receiving surface of the circle of radius r W at point p (x, y) The projected area at the corresponding three-dimensional model is equal to the area of the circle having the radius r;
(2-3)对于照度图上点p(x,y)，其相应的投影接收面上圆半径为rw，定义点p(x,y)处密度
Figure PCTCN2016079176-appb-000011
经过归一化，得到密度图M。(2-3) For the point p(x, y) on the illuminance map, the radius of the corresponding projection receiving surface is r w , and the density at the point p (x, y) is defined.
Figure PCTCN2016079176-appb-000011
After normalization, a density map M is obtained.
所述步骤(3)中，计算最优圆排列，具体包括以下步骤：In the step (3), the optimal circle arrangement is calculated, which specifically includes the following steps:
(3-1)计算最优目标圆个数N＝ρmN0/ρ0，其中ρm为密度图M的累加密度值，ρ0为B0的累加密度值，N0为B0的圆个数。(3-1) calculating an optimum target round number N = ρ m N 0 / ρ 0, wherein ρ m is the density of map M cumulative density values, ρ 0 is the cumulative density value of B 0, N 0 is of B 0 Round number.
(3-2)给定密度图M和目标圆个数N，通过调用de Goes的方法计算CCVT。在结果CCVT的每个Voronoi区域中计算其最大内切圆得到一组圆排列。对于某个Voronoi区域，其对应密度图M的某一块区域，将该区域中所有点对应的半径rw的均值作为该区域的期望圆大小。由此判定该区域的最大内切圆是否达到其期望圆大小。以目标圆个数N为上界，通过二分搜索查找达到最优圆正确率的CCVT。如图5，5(a)为用户输入灰度图像，5(b)为对应密度图M，5(c)为de Goes的方法计算得到的CCVT，5(d)为计算得到的最优圆排列。(3-2) Given the density map M and the number of target circles N, the CCVT is calculated by calling the method of de Goes. Calculating its maximum inscribed circle in each Voronoi region of the resulting CCVT results in a set of circular arrangements. For a certain Voronoi region, it corresponds to a certain block region of the density map M, and the mean value of the radius r w corresponding to all the points in the region is taken as the desired circle size of the region. It is thus determined whether the maximum inscribed circle of the region has reached its desired circle size. With the target circle number N as the upper bound, the binary search is used to find the CCVT that achieves the optimal circle correct rate. As shown in Fig. 5, 5(a) is the user input grayscale image, 5(b) is the corresponding density map M, 5(c) is the CCVT calculated by the method of de Goes, and 5(d) is the calculated optimal circle. arrangement.
所述步骤(4)中，确定相应孔洞的大小，位置和相对光源朝向角度，具体包括以下步骤：In the step (4), determining the size, position and relative light source orientation of the corresponding holes, specifically including the following steps:
(4-1)由步骤3可得达到最优正确率的CCVT及其对应圆排列。密度图M中任一点记录了其对应的密度值ρw(x,y)及其期望的孔洞类型，扩大型孔洞或者倾斜型孔洞。对于某个Voronoi区域，其对应的密度图的某块区域，将该区域中多数点表面的孔洞类型作为该Voronoi区域对应孔洞的类型。 (4-1) From step 3, the CCVT and its corresponding circle arrangement can be obtained to achieve the optimal correct rate. Any point in the density map M records its corresponding density value ρ w (x, y) and its desired hole type, enlarged hole or inclined hole. For a certain Voronoi region, a certain region of the corresponding density map, the hole type of the surface of most points in the region is taken as the type of the corresponding hole of the Voronoi region.
(4-2)如图6所示，若某Voronoi区域对应的孔洞类型为扩大型孔洞，该Voronoi区域的最大内切圆为D，D中嵌套一个缩小一个安全距离的D′，安全距离为0.5dmin。将D′通过中心投影的方式投影在三维壳状模型的外表面和内表面，分别形成两个相交的椭圆，使用一个圆柱形结构连接内外表面的椭圆形成该扩大型孔洞。(4-2) As shown in Fig. 6, if the hole type corresponding to a Voronoi region is an enlarged hole, the maximum inscribed circle of the Voronoi region is D, and a D' in which a safety distance is reduced is nested in D, the safety distance is It is 0.5d min . D' is projected on the outer and inner surfaces of the three-dimensional shell-like model by means of a central projection, respectively forming two intersecting ellipses, and the enlarged-type holes are formed by connecting an ellipse of the inner and outer surfaces with a cylindrical structure.
(4-3)如图6所示，若某Voronoi区域对应的孔洞类型为倾斜型孔洞，在D′选一个随机方向放置两个半径满足可打印条件的最小圆D1和D2，分别将D1和D2通过中心投影的方式投影在三维壳状模型的内表面和外表面，分别形成两个相交的椭圆，使用一个圆柱形结构连接内外表面的椭圆形成该倾斜型孔洞。(4-3) As shown in Fig. 6, if the hole type corresponding to a Voronoi region is an inclined hole, a minimum circle D 1 and D 2 whose radius meets the printable condition is placed in a random direction in D', respectively D 1 and D 2 are projected on the inner and outer surfaces of the three-dimensional shell-like model by means of a central projection, respectively forming two intersecting ellipses, and the slanted holes are formed by connecting an ellipse of the inner and outer surfaces with a cylindrical structure.
生成半色调投影模型后，可依照该模型直接进行3D打印。打印得到的模型可使光源透过模型上的孔洞在投影接收面上投影形成一幅与用户给定灰度图像最相近的连续灰度图像。After the halftone projection model is generated, 3D printing can be directly performed according to the model. The printed model allows the light source to be projected through the holes in the model on the projection receiving surface to form a continuous grayscale image that is closest to the given grayscale image of the user.
上述虽然结合附图对本发明的具体实施方式进行了描述，但并非对本发明保护范围的限制，所属领域技术人员应该明白，在本发明的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。 The above description of the specific embodiments of the present invention has been described with reference to the accompanying drawings, but it is not intended to limit the scope of the present invention. Those skilled in the art should understand that the skilled in the art does not require the creative work on the basis of the technical solutions of the present invention. Various modifications or variations that can be made are still within the scope of the invention.

Claims (10)

  1. 一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：包括以下步骤：An improved halftone projection and model generation method for 3D printing, comprising: the following steps:
    (1)根据给定的三维模型、投影接收面以及光源的基本参数，利用投影模拟方法，计算投影参考图像；(1) calculating a projection reference image by using a projection simulation method according to a given three-dimensional model, a projection receiving surface, and basic parameters of the light source;
    (2)将给定的灰度图像进行校正，得到照度图，将照度图的点与投影参考图像相对应，进行归一化处理后，计算得到密度图；(2) correcting the given grayscale image to obtain an illuminance map, corresponding to the projected reference image, and performing normalization processing to calculate the density map;
    (3)利用带容量约束的Voronoi划分的方法，根据密度图进行最优圆排列；(3) Using the Voronoi partitioning method with capacity constraints, the optimal circle arrangement is performed according to the density map;
    (4)根据最优圆排列确定的孔洞类型，确定相应孔洞的大小、位置和相对光源朝向角度，生成物理模型，进行3D打印。(4) According to the hole type determined by the optimal circle arrangement, determine the size, position and relative light source orientation angle of the corresponding hole, generate a physical model, and perform 3D printing.
  2. 如权利要求1所述的一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：所述步骤(1)中，投影模拟方法的具体方法包括：An improved 3D printing-oriented halftone projection and model generation method according to claim 1, wherein in the step (1), the specific method of the projection simulation method comprises:
    (1-1)将光源离散化为若干个点光源；(1-1) discretizing the light source into a plurality of point sources;
    (1-2)在投影接收面的投影区域离散采样为若干个的投影接收点；(1-2) discretely sampling the projection area of the projection receiving surface into a plurality of projection receiving points;
    (1-3)在多孔模型的遮挡作用下，计算投影接收点的所有点光源的总辐射照度；(1-3) Calculating the total irradiance of all point sources of the projection receiving point under the occlusion of the porous model;
    (1-4)对投影接收点的总辐射照度进行校正，得到投影模拟图像灰度值。(1-4) Correcting the total irradiance of the projection receiving point to obtain the gradation value of the projected analog image.
  3. 如权利要求2所述的一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：所述步骤(1-4)中校正方法为通过Gamma校正。An improved 3D printing-oriented halftone projection and model generation method according to claim 2, wherein the correction method in the step (1-4) is corrected by Gamma.
  4. 如权利要求1所述的一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：所述步骤(1)中，计算投影参考图像的具体方法包括：The improved 3D printing-oriented halftone projection and model generation method according to claim 1, wherein in the step (1), the specific method for calculating the projection reference image comprises:
    (1-a)在三维模型上紧密排列满足可打印性条件半径最小的孔洞；(1-a) closely aligning the holes satisfying the minimum radius of printability conditions on the three-dimensional model;
    (1-b)使用投影模拟方法，在孔洞最紧密排列的约束条件下，生成投影参考图像。(1-b) Using the projection simulation method, a projection reference image is generated under the constraint that the holes are most closely arranged.
  5. 如权利要求1所述的一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：所述步骤(1-a)中，满足可打印性条件指的是三维模型上的孔洞的半径不小于选择的打印技术中能够打印的孔洞的最小半径，同时，孔洞之间的距离不小于该打印技术可打印两孔洞间的最小距离。An improved 3D printing-oriented halftone projection and model generation method according to claim 1, wherein in the step (1-a), the printability condition is satisfied to refer to a hole in the three-dimensional model. The radius is not less than the minimum radius of the holes that can be printed in the selected printing technique, and the distance between the holes is not less than the minimum distance between the two holes that the printing technique can print.
  6. 如权利要求1所述的一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：所述步骤(2)中，对于输入的任一灰度图像，计算密度图，具体包括以下步骤：An improved 3D printing-oriented halftone projection and model generation method according to claim 1, wherein in the step (2), for any grayscale image input, a density map is calculated, specifically including The following steps:
    (2-1)将给定的灰度图像，通过逆校正得到照度图；(2-1) obtaining a illuminance map by inverse correction of a given grayscale image;
    (2-2)对于照度图的每一点，根据其辐射照度设定对应的三维模型处遮光率，根据遮光率计算其对应的目标圆半径； (2-2) For each point of the illuminance map, the shading rate of the corresponding three-dimensional model is set according to the irradiance of the illuminance, and the corresponding target circle radius is calculated according to the shading rate;
    (2-3)将照度图每点对应的目标圆半径映射为相应密度值，进行归一化处理，得到密度图。(2-3) The target circle radius corresponding to each point of the illuminance map is mapped to the corresponding density value, and normalized to obtain a density map.
  7. 如权利要求1所述的一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：所述步骤(3)中，计算最优圆排列的具体方法包括以下步骤：An improved 3D printing-oriented halftone projection and model generation method according to claim 1, wherein in the step (3), the specific method for calculating the optimal circle arrangement comprises the following steps:
    (3-1)根据密度图的累加密度值、投影参考图像的圆个数和累加密度值，计算最优目标圆个数；(3-1) calculating the optimal target circle number according to the accumulated density value of the density map, the number of circles of the projected reference image, and the accumulated density value;
    (3-2)根据密度图和目标圆个数，通过二分搜索查找达到最优圆正确率的带容量约束的Voronoi划分。(3-2) According to the density map and the number of target circles, the Voronoi partition with capacity constraint that achieves the optimal circle correct rate is searched by binary search.
  8. 如权利要求1所述的一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：所述步骤(4)中，具体方法包括：An improved 3D printing-oriented halftone projection and model generation method according to claim 1, wherein in the step (4), the specific method comprises:
    (4-1)根据得到的最优正确率的Voronoi区域和其对应圆排列，将其对应的密度图的区域中多数点标明的孔洞类型作为该Voronoi区域期望的孔洞类型；(4-1) according to the obtained optimal correct rate of the Voronoi region and its corresponding circle arrangement, the hole type indicated by the majority of the points in the corresponding density map region is taken as the desired hole type of the Voronoi region;
    (4-2)根据任意Voronoi区域期望的孔洞类型，在其三维模型的对应位置生成相应的孔洞。(4-2) According to the type of hole desired in any Voronoi region, a corresponding hole is generated at a corresponding position of the three-dimensional model.
  9. 如权利要求8所述的一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：所述步骤(4-2)中，若某Voronoi区域对应的孔洞类型为扩大型孔洞，生成方法为：在该Voronoi区域的最大内切圆嵌套一个缩小一个安全距离的内切圆，将该内切圆通过中心投影的方式投影在三维壳状模型的外表面和内表面，分别形成两个相交的椭圆，使用一个圆柱形结构连接内外表面的椭圆形成该扩大型孔洞。An improved 3D printing-oriented halftone projection and model generation method according to claim 8, wherein in the step (4-2), if a hole type corresponding to a Voronoi region is an enlarged hole, The generating method is: nesting an inscribed circle which is reduced by a safe distance in a maximum inscribed circle of the Voronoi region, and projecting the inscribed circle by a central projection on the outer surface and the inner surface of the three-dimensional shell model, respectively forming Two intersecting ellipses are formed by joining an ellipse of the inner and outer surfaces using a cylindrical structure to form the enlarged holes.
  10. 如权利要求8所述的一种改进的面向3D打印的半色调投影与模型生成方法，其特征是：所述步骤(4-2)中，若某Voronoi区域对应的孔洞类型为倾斜型孔洞，生成方法为：在该Voronoi区域的最大内切圆嵌套一个缩小一个安全距离的内切圆，在该内切圆选一个随机方向放置两个半径满足可打印条件的最小圆D1和D2，分别将D1和D2通过中心投影的方式投影在三维壳状模型的内表面和外表面，分别形成两个相交的椭圆，使用一个圆柱形结构连接内外表面的椭圆形成该倾斜型孔洞。 An improved 3D printing-oriented halftone projection and model generation method according to claim 8, wherein in the step (4-2), if a hole type corresponding to a Voronoi region is an inclined hole, The generating method is: nesting an inscribed circle which is reduced by a safe distance in the largest inscribed circle of the Voronoi region, and selecting a minimum circle D 1 and D 2 in the inscribed circle to select two radii satisfying the printable condition. D 1 and D 2 are respectively projected on the inner surface and the outer surface of the three-dimensional shell model by central projection, respectively forming two intersecting ellipse, and the inclined holes are formed by connecting an ellipse of the inner and outer surfaces with a cylindrical structure.
PCT/CN2016/079176 2016-03-16 2016-04-13 Improved method for generating halftone projection and model for 3d printing WO2017156815A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201610150875.9 2016-03-16
CN201610150875.9A CN105844712B (en) 2016-03-16 2016-03-16 A kind of improved halftoning projection and model generating method towards 3D printing

Publications (1)

Publication Number Publication Date
WO2017156815A1 true WO2017156815A1 (en) 2017-09-21

Family

ID=56587933

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2016/079176 WO2017156815A1 (en) 2016-03-16 2016-04-13 Improved method for generating halftone projection and model for 3d printing

Country Status (2)

Country Link
CN (1) CN105844712B (en)
WO (1) WO2017156815A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111757790A (en) * 2018-02-20 2020-10-09 本地汽车知识产权有限责任公司 Method and apparatus for additive manufacturing
CN114290683A (en) * 2022-01-07 2022-04-08 宁波智造数字科技有限公司 Photocuring 3D printing uniformity adjusting process method utilizing cylinder

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108257112B (en) * 2017-12-27 2020-08-18 北京七鑫易维信息技术有限公司 Method and device for filtering light spots
CN111914974B (en) * 2020-07-07 2022-02-18 山东大学 Three-dimensional two-dimensional code generation method and system for 3D printing
CN113261973B (en) * 2021-04-13 2022-09-02 山东师范大学 Method and system for identifying atrial fibrillation electrocardiosignals through weighted multi-scale finite-crossing visual image
CN113946987B (en) * 2021-12-20 2022-03-22 浙江闪铸三维科技有限公司 3D model nesting and placing method based on model similarity

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110075944A1 (en) * 2009-09-28 2011-03-31 Yuan Xiaoru Sample based tone mapping method for high dynamic range images
CN103978690A (en) * 2014-05-28 2014-08-13 山东大学 Method for optimizing internal structure of 3D (Three-Dimensional) printed object
CN105023296A (en) * 2014-08-22 2015-11-04 山东大学 3D-printing-oriented halftone projection and model generation method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140032448A1 (en) * 2012-07-24 2014-01-30 Telefonica, S.A. Method, computer programs and a use for the prediction of the socioeconomic level of a region

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110075944A1 (en) * 2009-09-28 2011-03-31 Yuan Xiaoru Sample based tone mapping method for high dynamic range images
CN103978690A (en) * 2014-05-28 2014-08-13 山东大学 Method for optimizing internal structure of 3D (Three-Dimensional) printed object
CN105023296A (en) * 2014-08-22 2015-11-04 山东大学 3D-printing-oriented halftone projection and model generation method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
LIANG, GUANGHUI ET AL.: "Poisson disk sampling through disk packing", COMPUTATIONAL VISUAL MEDIA, vol. 1, no. 1, 31 March 2015 (2015-03-31), XP055424061 *
LU , LIN ET AL.: "Build-to-Last:Strength to Weight 3D printed Objects", ACM TRANSACTION ON GRAPHICS, vol. 33, no. 4, 31 July 2014 (2014-07-31), XP058051982 *
ZHAO, - HAISEN ET AL.: "Printed Perforated Lampshades for Continuous Projective Images", ARXIV E-PRINTS ABSTRACT SERVICE, 31 October 2015 (2015-10-31), XP058259728 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111757790A (en) * 2018-02-20 2020-10-09 本地汽车知识产权有限责任公司 Method and apparatus for additive manufacturing
CN114290683A (en) * 2022-01-07 2022-04-08 宁波智造数字科技有限公司 Photocuring 3D printing uniformity adjusting process method utilizing cylinder

Also Published As

Publication number Publication date
CN105844712A (en) 2016-08-10
CN105844712B (en) 2018-11-13

Similar Documents

Publication Publication Date Title
WO2017156815A1 (en) Improved method for generating halftone projection and model for 3d printing
US11657470B2 (en) Generating signal bearing art using Stipple, Voronoi and Delaunay methods and reading same
CN107206696B (en) Three-dimensional object substructure
CN101322153B (en) Adjusting digital image exposure and tone scale
US8681383B2 (en) 2nd generation dot-off-dot stochastic halftone
US8547598B2 (en) Polygonal-boundary-based halftone method
US6710778B2 (en) Method for halftoning using a difference weighting function
CN105023296B (en) A kind of halftoning projection and model generating method towards 3D printing
CN1734503A (en) Stretch-driven mesh parameterization method using spectral analysis
JP5807489B2 (en) Image processing apparatus, image processing method, and program
EP3458893B1 (en) Optic, luminaire, method for fabricating optic and non- transitory computer readable storage medium
Zhao et al. Printed perforated lampshades for continuous projective images
US20110157654A1 (en) System and method for halftoning using a parametrically controlled hexagonal halftone dot shape threshold function
CN110001066B (en) Method for determining filling direction of minimum partition in three-dimensional printing
US20130148168A1 (en) Visible and invisible watermarking of printed images via 2nd generation stochastic seed frequency modulation
JP6486096B2 (en) Print data generation apparatus, print data generation method and program
CN108389259B (en) Random center aggregation image mesh tone method and system
CN105335978A (en) Image edge detection method and system
Gooran et al. 3D surface structures and 3D halftoning
Stoppel et al. LinesLab: A Flexible Low‐Cost Approach for the Generation of Physical Monochrome Art
CN109154991B (en) Three-dimensional object substructure
US8576448B2 (en) Clustered halftone generation
Abedini et al. 3D halftoning based on iterative method controlling dot placement
US10587774B2 (en) 3D printed object halftone image generation containing updated voxel data
WO2023103548A1 (en) Automobile parametric texture generation method and system, and storage medium

Legal Events

Date Code Title Description
NENP Non-entry into the national phase

Ref country code: DE

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 16893978

Country of ref document: EP

Kind code of ref document: A1

122 Ep: pct application non-entry in european phase

Ref document number: 16893978

Country of ref document: EP

Kind code of ref document: A1