1 Introduction
The emergence of 3D printing technology opens new interesting opportunities for computerized halftoning. In this paper, we introduce the generation of 3Dprinted perforated lampshades designed to project grayscale images onto surrounding surfaces. The spatial distribution of tone across the projection is generated by carefully controlling the amount of light shining through the printed surface. Similarly to the dots used in halftoning, our basic idea is to perforate the physical 3D surface of the lamp with small tubes, whose parameters are used to control the amount and the spatial distribution of light that emanates from the lamp and reaches a receiving surface. Figure 1 shows an example of such a lampshade and its projected image. Projecting continuous tone imagery via a perforated surface can thus be regarded as halftoning with light or 3D halftoning.
Halftoning is a technique used to represent continuous shades of gray through the use of discrete dots of the same color, varying either in size, in shape, or in spacing. The method relies on the natural lowpass spatial filtering of the human visual system that blends the discrete dot pattern into a continuous tone. Traditional halftoning uses simpleshaped dots, typically circular, of continuously varying size. However, when the dots all have the same size, spatial resolution can be traded for perceived tone resolution [Ulichney (1988), Lau and Arce (2008)].
The 3D halftoning technique that we present differs from ordinary halftoning in a number of ways. First, unlike other common digital media, here one can generate dots of continuous sizes. In that sense, the technique is closer to analog halftoning, where the dot sizes are continuous. Second, dots are holes in a surface, realized as tiny 3D tubes, hence having both radius and direction with respect to the light source. Third, the resulting image is not printed but projected, which requires to consider the geometry of both the projecting surface (the lampshade) and the receiving surface (the wall). Lastly, since the projected light footprint of each tube is slightly blurred, and multiple footprints add up in areas of overlap, there is already some lowpass filtering inherent in the image formation process. We explicitly account for and take advantage of the effect of overlapping footprints of adjacent tubes, which is inhibited in digital halftoning.
Moreover, it should be stressed that in our setting, the major challenge is to ensure the printability and structural integrity of the perforated lampshade. Specifically, there are strict lower bounds on the radii of the tubes, and on the intertube spacing. Violating the former constraint would result in clogged tubes, and the latter would render the shell fragile and prone to breakage. In particular, having a lower bound on the tube radius means that darker tones cannot be achieved by simply using smaller holes; instead, we reduce the amount of light passing through a tube by tilting it away from the light source center.
We present a 3D halftoning technique that, given the geometry of the lampshade surface, a target grayscale image, and a receiving surface, produces a spatial distribution of tubes, along with their radii and directions (tilt angles), such that the resulting projected image faithfully reproduces the input image. The goal is to produce a continuous tone projected image, which strives to match the distribution of tones and the spatial detail of the original.
Our approach attempts to match the target image intensity at each projected location by placing tubes with suitable radii and tilt angles around the corresponding location on the lampshade. Brighter areas are reproduced using tubes with a larger radius, while in darker areas we place minimal radius tubes, tilted away from the light source center. In order to maximize the spatial resolution, the tubes must be placed as densely as possible, but they must not violate the intertube distance constraint. To achieve this, we embed each tube inside a disk that incorporates a safety margin around the tube. Note that differently from ordinary halftoning or stippling, in our case both brighter and darker areas require larger disks, with the maximal density of disks corresponding to the middle of the grayscale range.
Having reduced the problem to one of finding a dense packing of disks with spatially varying radii, we solve it by computing a capacityconstrained Voronoi diagram over a suitably defined density function. The density function for a particular target image is determined via lookup from a series of precomputed simulated projected images, each corresponding to a different uniform density tube pattern on the lampshade. The entire process is illustrated in Figure 2.
In summary, the contributions of this paper are as follows:

We tackle the novel problem of generating projected imagery by shining light through a 3Dprinted perforated surface.

We introduce a novel 3D halftoning approach, where the halftoning dots are realized as a distribution of tubes with varying radii and orientation, passing through a solid surface.

We present a method that determines the spatial distribution and the parameters of the individual tubes, while striving to match the target image, subject to fabrication constraints.
Our results demonstrate the effectiveness of the proposed technique using a variety of target images. The limitations of the process are also demonstrated and discussed.
2 Related Work
Halftoning is a classical technique which played a major role in traditional paper printing and in digital displays [Kipphan (2001)]. Our work leverages this classical technique in the new domain of digital fabrication and 3D printing. Throughout the evolution of halftoning, the dominating issues have persistently been ones of structure preservation, tone reproduction, point density and spatial resolution.
2.1 Halftoning and Stippling
As a mature research problem, digital halftoning is well studied with numerous variants and solutions [Pang et al. (2008), Schmaltz et al. (2010)]. However, fewer researchers have addressed the problem of structure preserving halftoning and stippling. Kim et al. [Kim et al. (2008)] present a dot placement algorithm that guides the dot placement along feature flow extracted from the feature lines, resulting in a dot distribution that conforms to feature shapes. Pang et al. [Pang et al. (2008)] introduce an optimizationbased halftoning technique that preserves both the structure and tone metrics by optimizing for both features. Chang et al. [Chang et al. (2009)] further extend the iterative optimization to an error diffusion method, which is more efficient.
Li et al. [Li et al. (2010)] show an efficient application for anisotropic blue noise sampling to stippling of images. Li and Mould [Li and Mould (2011)] present an automatic structureaware stippling algorithm. The method starts from a contrastaware halftoning result and reduces the number of dots using a nonlinear prioritybased error adjustment function.
Balzer et al. [Balzer et al. (2009)] propose using the capacityconstrained Voronoi tessellation (CCVT) for enforcing the constraint of equal weighted area for the region around each point in a stippled image. However, the method they propose suffers from high computational complexity. Consequently, considerable attention was given to developing faster alternatives. In particular, the remarkable work of de Goes et al. [de Goes et al. (2012)] shows that CCVT can be formulated as a constrained optimal transport problem. This results in a fast continuous constrained minimization method, which is able to enforce the capacity constraints exactly. In this work, we show how to encode the objectives and the requirements of our 3D halftoning method as a density function, which enables solving the problem by constructing a CCVT, and use the method of de Goes et al. for this purpose.
The works mentioned earlier are concerned with digital halftoning and stippling of 2D images, while our work focuses on halftoning of 3D surfaces (with finite thickness) in the context of 3D printing. In this context, Stucki [Stucki (1997)] first introduced the idea of 3D digital halftoning for transforming continuousdensity objects into binary representations for rendition with additive fabrication technologies. Lou and Stucki [Lou and Stucki (1998)] further adapt the ordereddither and error diffusion algorithms to the 3D case. Zhou and Chen [Zhou and Chen (2009)] utilize 3D digital halftoning, to reduce the fabricating time in layered manufacturing. In contrast to these works, we focus on 3D halftoning for the purpose of forming projected imagery that reproduces the tones of a target image, subject to fabrication constraints.
2.2 Opticsrelated Fabrication
Recently, design or generation of illumination effects via geometric modulation has been drawing increasing attention. Mitra and Pauly [Mitra and Pauly (2009)] introduce shadow art, an algorithm for computing a 3D volume, whose shadows best approximate multiple binary images. Subsequent approaches generalize this to colored shadows using volumetric objects manufactured by transparent acrylic [Wetzstein et al. (2011), Baran et al. (2012)]. In general, they build transmittance functions to simulate light attenuation through multilayered attenuators.
Our work is inspired by the technique of Alexa and Matusik [Alexa and Matusik (2012)], who drill holes with varying depths on a surface to induce the given image based on the occlusion of small holes.
Research in this area has also been addressing light reflections and refractions using surface modeling methods. These works produce a desired caustic image by casting a caustic, modulating the geometry using microfacets [Weyrich et al. (2009)], micropatches [Papas et al. (2011)], Bspline surfaces [Finckh et al. (2010)], continuous surfaces [Kiser et al. (2013), Yue et al. (2014)] or normal fields [Schwartzburg et al. (2014)], and then milling or engraving the surfaces.
Researchers have also shown light effects in a wide range of applications, such as steganography and appearance design. Papas et al. [Papas et al. (2012)] develop passive display devices to hide images that can be decoded with fabricated lens. Malzbender et al. [Malzbender et al. (2012)] present a method for printing the 4D reflectance function of a surface for a fixed viewing direction. Levin et al. [Levin et al. (2013)] use wave optics instead of geometric optics to produce custom designed reflectance functions with high spatial resolution. Lan et al. [Lan et al. (2013)] design the appearance fabrication by spatiallyvarying changes to both local shading frames and reflectance.
Printing optical fibers that control light propagation through total internal reflection between two surfaces has drawn some research attention thanks to the cutting edge multimaterial 3D printers. Willis et al. [Willis et al. (2012)] design optical fibers for customizing interactive devices in display and illumination. Pereira et al. [Pereira et al. (2014)] introduce automatic fiber routing algorithm that minimizes curvature and compression for optimal light transmission.
All of the above are highend techniques, which rely on expensive manufacturing equipment and specific materials. In contrast, our method is designed for standard 3D printing, which is rapidly becoming an offtheshelf, popular technology. Furthermore, for widely used materials, such as plastic and powder, the transmittance of light is very weak, as light mostly scatters due to the surface roughness. Therefore, light effects based on surface reflection and refraction cannot be easily adapted to such materials.
3 Perforated lampshade design
Recall that our goal is to design a 3Dprintable perforated lampshade, such that the light emanating from the lamp forms a continuous tone image, which is as close as possible to a target grayscale input. Light emanates from the lamp through holes in the lamp’s surface. Due to the finite thickness of the lamp’s shell, the holes may be thought of as tiny tubes, whose density, radii, and orientations vary across the surface, as may be seen in Figure 3. Since the total distribution of light across the projected image is determined by the combined effect of light passing through the tubes, our challenge is to determine the tubes’ parameters, while respecting fabrication constraints imposed by the need to obtain a printable and structurally sound surface.
Specifically, given a target grayscale image our task is to configure a set of tubes perforating the lampshade shell, such that the following objectives and requirements are satisfied:

The projected image on a given diffuse surface, henceforth referred to as the wall, closely approximates the target grayscale values across the projection. Note that the approximation is up to some scaling factor, since we have control over the total amount of luminous flux emitted by the light source, as well as over the exposure time, when capturing the projection with a camera.

The projected image should exhibit continuous tones, while resolving fine spatial detail, as much as possible. In order to achieve this objective the density of the tubes should be maximized, while their radii should be minimized.

Fabrication constraints: (a) the radius of a tube is bounded below by ; and (b) any two adjacent tubes must have a gap of width greater than of solid material between them. In our current setting, and .
The densest arrangement of tubes that satisfies the fabrication constraints above is shown in Figure 4(a). This arrangement corresponds to a hexagonal packing of disks of radius , with a tube of radius embedded at the center of each disk. The resulting projected pattern is shown in Figure 4(c). In order to match a target image , the tube arrangement must be adjusted in a spatially varying fashion. In order to achieve brighter tones we must increase the tube radii, but darker tones cannot be achieved by decreasing the radii, since this would violate the first fabrication constraint above.
One alternative is to employ a stippling approach to generate a variable density set of tube positions that would satisfy the fabrication constraints above. However, such an approach necessitates spacing tubes much farther apart from each other in the darker areas of the image, resulting in projection of isolated light dots inside such areas. Thus, continuity and spatial detail resolution has to be severely compromised. Instead, we propose decreasing the amount of light passing through tubes in darker areas by tilting their directions away from the light source center. The tilting direction is random, in order to avoid undesirable structured artifacts. This solution avoids unnecessarily isolating tubes, and plays an essential role in achieving continuous gray tones in the projected imagery. It is worth mentioning that the tilting angle is not restricted by fabrication capabilities, which offers enough variable tuning range.
The inadequacy of a stipplingbased approach is demonstrated in Figure 5. Here, starting from an input image (a), we use a stateoftheart stippling method [de Goes et al. (2012)] to generate the pattern in (b). In order to satisfy the fabrication constraints, the number of dots in this pattern is limited to roughly 3K. This results in a (simulated) projected image with poor continuity and isolated light spots, especially inside the darker tone areas around the dog’s eyes and nose. Using our proposed solution, we are able to distribute roughly 6K tubes across the lampshade, thereby achieving the same apparent grayscale range, but in a much more continuous fashion.
Note that in order to ensure a minimal intertube distance, both widened tubes and tilted tubes must be placed inside larger disks on the lampshade surface, as illustrated in Figure 6. The radius of each disk is chosen so as to accommodate the widened or tilted tube, in addition to a safety margin of width . In summary, each target intensity at a particular location corresponds to a certain disk radius, and our goal is to distribute a set of tightly packed nonoverlapping disks with the desired radii across the relevant portion of the lampshade surface.
We attempt to achieve the required disk arrangement by computing a capacityconstrained Voronoi tessellation over a suitably defined density function, and inscribing a maximal circle inside each tessellation cell. The density function for a particular target image is determined via lookup from a series of precomputed simulated projected images, each corresponding to a different uniform density tube pattern on the lampshade; one such image is shown in Figure 4(b). The entire design process is visualized in Figure 2 and summarized in Algorithm 1. The main stages of the process are described in more detail in the remainder of this section.
3.1 Simulating reference images
We compute a series of simulated reference images , each of which predicts the spatially varying reflected light intensity across the projection domain that results from a dense uniform arrangement of tubes with a given fixed radius and tilting angle. The reference images are precomputed and stored once for each setup (lampshade and wall geometry). The simulation process is described in detail in Section 3.5.
Reference image corresponds to the densest arrangement of the smallest size tubes (radius ) through the lampshade shell. The arrangement, the result of the simulation, and a photograph of the actual projected image are shown in Figure 4. Images with correspond to a similarly arranged grid of wider tubes, thus resulting in brighter images. Images with correspond to grids of tilted tubes, resulting in darker images.
Given a target image , these reference images enable us to determine the desired tube parameters at each position on the lampshade, in turn defining a density function on the wall. This density function is then used to compute a capacityconstrained Voronoi tessellation.
3.2 Computing the density function
We use the reference images to determine the spatially varying density on the wall. Since the projected image is formed on the wall by a linear combination of light transmitted through the tubes, we first linearize the given target image by applying an inverse gamma correction. We assume that the target image is gamma corrected with a standard gamma value of 2.2.
For each projected image pixel , whose target grayscale value is given by , we determine , such that . Let and denote the radii of the disks in which the tubes were embedded when generating and . This means that the desired target value at position could be achieved by embedding a tube inside a disk of radius , , and we determine the value of using linear interpolation.
In our current implementation, we compute the capacityconstrained Voronoi tessellation on the planar wall domain. Thus, we need to convert the desired radii on the lampshade to ones on the wall. Specifically, given a radius of a disk on the lampshade, we compute the radius of a disk on the wall, such that its projection back onto the lampshade (ellipse in the spherical case) has the same area as the original disk of radius .
Given the target disk radii across the projection area, we set the density to be proportional to the inverse of the disk area, . The density values are normalized, such that . In practice, we impose an upper bound of 1.3mm on the tube radius. Thus, the density function is bounded from below, with the actual bound varying across the projection surface.
3.3 Computing the disk distribution
Given the density function , and a desired number of disks , we use the method of de Goes et al. [de Goes et al. (2012)] to compute a capacityconstrained Voronoi tessellation, and fit a maximal inscribed disk inside each of the resulting tessellation cells.
Although the above method is able to enforce the capacity constraints exactly, the resulting Voronoi cells are not necessarily hexagonal, and thus some of the maximal inscribed disks may deviate from their intended size. We estimate an initial value for
by computing , the amount of density per disk in the dense reference pattern , and dividing the integral of the target density function by . We found that this estimate typically results in a large percentage of disks being too small. We thus employ a binary search to find the number of disks for which the percentage of disks which achieve their intended radius (within 0.05mm) is greatest.While most disks are acceptably close to their desired radius, there are also disks which are too large or too small. Consequently, the amount of light transmitted through the embedded tubes could be smaller or larger than the desired intensity. Larger disks are handled by reducing the embedded tube radius or tilt angle, as needed to match the originally intended ones. Smaller disks are resolved by removing a fraction of them. For each group of such disks, whose Voronoi cells are connected, we remove disks randomly one by one. After each disk removal, we enlarge the group of cells along the onering neighborhood to include their immediately neighboring cells, and locally relax the remaining disks in this area via centroidal Voronoi tessellation using Lloyd’s method.
3.4 Tube generation
For each disk centered at with radius on , we first subtract the safety margin to obtain a smaller concentric disk . In the case where , we generate a tube directed at the center of the light source by intersecting the oblique circular cone defined by and the light source center with the lampshade shell, as shown in Figure 6.
If , a minimal radius tilted tube must be generated. We select a random diameter of and place two random minimal radius disks and near the ends of this diameter. and are then projected onto the outer and the inner surfaces of the lampshade shell, respectively (see Figure 6). The tilted tube is formed by connecting the two resulting contours.
At this point, a printable 3D model of the perforated lampshade may be generated and printed by a 3D printer.
3.5 Projected image simulation
The projected image is formed by the light emanating from the light source inside the lamp, passing through the tubes in the lampshade and reaching a receiving surface (wall). Thus, the amount of reflected light at each point on the wall depends on the emission characteristics of the light source, the geometric relationship between the light source and the receiving point, the visibility between them, and the reflectance characteristics of the wall.
The light source that we use in this work is a small and powerful COB (Chips on Board) LED, which is roughly diskshaped with a diameter of 9mm. In our simulations, we represent this area light source as a collection of point light sources , each emitting a luminous flux of , where is the luminous flux for the entire light source, measured in lumens. We use in our results.
For a Lambertian light source and a Lambertian reflector, the total illuminance (in lux) reaching a point on the receiving surface is given by
(1) 
where is the distance from to , and and are the angles between the line connecting the two points and the two normals, as shown in Figure 7. is the visibility between and .
In practice, we found that our light source is not Lambertain, and its emitted luminance diminishes as the angle increases. Furthermore, we found that this behavior is slightly anisotropic, and the luminance diminishes more quickly in the vertical direction. Empirically, we found that the nonLambertian behavior and the anisotropy are still well modeled by Equation (1), provided that the horizontal and vertical coordinates of the point (with respect to the center of the image) are each scaled by an appropriate factor. Specifically, the scaling factors that we use are 1.7 and 1.9 for the horizontal and the vertical coordinate, respectively.
The accuracy of our empirical model is demonstrated in Figure 8. The red curves show a horizontal (left) and a vertical (right) cross section through the projected image shown in Figure 4(c), and the blue curves are the values predicted by our model.
We evaluate the visibility of points on the light source using ray casting (see Figure 7). Note that it is not necessary to intersect rays with the complete geometric model of the perforated lampshade shell. Instead, we quickly determine a small set of tubes that are relevant for a given ray, and check whether the ray passes through one of these tubes, without intersecting the tube’s surface. Note that the above approach accounts only for direct illumination, without taking into account reflections, subsurface scattering, etc. In order to assess the magnitude of these effects, we have sprayed the interior of lamps with black paint, and found the differences in the measured intensity very small. The main difference is that this reduces the amount of ambient light in the room, leading to projected images with slightly better contrast.
Accurate light source visibility estimation is particularly important for correctly simulating tilted tubes. Figure 9 shows the simulated and the actual light footprints of a single tube, for a variety of different tilt angles. The sum of the intensity across the area of the footprint is plotted below (normalized such that the maximal sum is 1.0). We can observe that the simulation predicts the illuminance reasonably well, with some deviations between the two curves due to imprecisions in the manufacturing process.
Finally, in order to visualize the simulated results we convert them into grayscale values by scaling the result into a range, and applying a gamma correction (we use gamma 2.2 in all our results).
4 Results
We have implemented the method described in the previous section, and have used the resulting system to design a number of lampshades for projecting a variety of target images. The results may be seen in Figures 1 and 11. The lampshades in all of these results have spherical geometry, with a diameter of 22cm. To save material we typically perforate each spherical lamp with two images on two opposing sides. As mentioned earlier, the light source we use is a Cree XLamp CXA1507 LED with 3000K color temperature, which is roughly diskshaped with a radius of 9mm. The perforated lamp projects its image onto a planar rear projection screen that offers high light transmission, located at a distance of 40cm from the light source. The physical size of the projected image is 100 100 cm. Figure 10 shows a photograph of our setup. The photographs of the projected images shown in this paper are shot from the opposite side of the screen for eliminating distortion, using a Canon EOS 5D Mark II, with exposure 1/80 sec, f/4.0, and ISO 400 settings, in a room without any additional light sources.
Model  Initial #disks  Final #disks  % correct size 

Marilyn  7323  5789  79.47 
CircularRamp  7914  6321  87.73 
Hepburn  7617  5945  82.50 
Kelly  7563  5914  80.74 
Marlon  7167  5607  81.34 
Toucan  6926  5352  75.76 
Lion  7243  5650  74.81 
Dog  7891  6055  82.12 
Statistics for our results. Initial disks is our initial estimate of the required number of disks, which serves as a starting point for a binary search, after which we achieve the final number of disks. The percentage of correctly sized disks among these are reported in the right column.
The resulting lampshade designs have been 3Dprinted to allow a qualitative evaluation of the projected imagery. The lampshades were printed on a powderbased binderjet printer (ProJet 660Pro). The thickness of the lampshade shells was set to 3mm. The manufacturing time for such a spherical lampshade includes 16.5 hours of printing and 1 hour of drying. The inside surface of the lamps were sprayed with black paint, which somewhat improves the contrast of the projected images.
Table 4 reports several statistics for each of these results. The computation time of our method consists of two main parts: computing the disk distribution and generating the 3D model of the lampshade perforated with tubes. For the lamps shown in the paper it takes less than 5 minutes to run the method of de Goes et al. [de Goes et al. (2012)] up to 10 times, while searching for the optimal number of disks. The resulting lamps have around 6000 tubes, and the 3D model generation takes roughly 3 minutes. For verification, we typically compute a simulated result before sending a lamp to be printed, which takes roughly 1 minute. All of the above times are measured on an Intel Core^{TM} i5 CPU 3.3GHz with 8GB RAM.
Figure 12 shows two qualitative comparisons between the target image, the simulated result, and the physical projected image. This is done on two examples: a concentric circles pattern featuring four regions of a constant grayscale value, and a more complex natural image. Below each image, we show the corresponding histogram, in order to better examine the differences between the corresponding tone distributions.
In the concentric ramp case (12a), it obviously infeasible to reproduce the ideal fourway distribution of grayscale values with our process. Our design process approximates this distribution with four wider peaks (the black peak coincides with the left edge of the histogram). A photograph of the actual projection features four similar peaks, although there are slight differences in their positions and spread. The target image of Marilyn Monroe (12h) has a nearly uniform histogram. This histogram is reproduced reasonably well by the our process, however, both the simulated and the actual projected image have more mass in the lower part of the grayscale range.
These two examples illustrate well some of limitations of this new medium. Achieving bright tones in the middle of the image is relatively easy (since the distance between the receiver and the light source is minimal there, and both cosines in Eq. (1) are near 1.0). However, reproducing dark tones in the central region is challenging because of the necessity to use large tilt angles. As explained earlier, extensive use of tilting requires reducing the tube density, thereby reducing spatial resolution. Conversely, achieving light tones at the periphery of the image is difficult, because of the falloff in the illuminance, and increasingly wider tubes are needed to produce a light tone in the periphery. Thus, there’s a tradeoff between our ability to reproduce the full tonal range, and our desire for better spatial resolution. In particular, by imposing an upper bound (1.3mm) on the embedding disk radii across the lamp surface, we are not able to match the lighter tones of the target image beyond a certain radius.
The above can be observed in a more quantitative manner in Figure 13, which shows the response of our method to several radial cosine waves with the frequency increasing from top to bottom row. Next to each input pattern we show the simulated result produced by our method and plot the amplitudes of the input wave (in red) and the result (blue), as a function of the radial distance. It may be observed that for all frequencies our method fails to reproduce the waveform beyond a certain radial distance, with the distance slightly decreasing as the frequency increases. It may also be seen that in the areas where the waveform is reproduced, the contrast response also decreases with frequency: while in the top row the ratio between the reproduced amplitude to the input one is roughly 0.89, this ratio drops to 0.40 in the bottom row.
5 Concluding Remarks
We have presented a technique to design printable perforated lampshades. The main challenge was to control the light emanating through the lampshade under the unique printability constraints and engineering setting. The technique that we developed extends classical halftoning to a novel domain, where the “dots” are not printed, or physically tangible, but projected, effectively halftoning with light.
Although we have not yet adapted our process to nonspherical lampshades, note that the process is general enough to accommodate such shapes. Specifically, given a nonspherical lampshade geometry, one would need to (a) distribute dense uniform patterns of tubes across the lampshade, in order to compute the simulated reference images ; and (b) determine an appropriate transformation between disks on the lampshade surface and corresponding disks on the projection surface. The rest of our process would remain unchanged.
The 3D printing technology is emerging and growing quickly, and many new intriguing applications are introduced. We believe that the combination of light and 3D printing has much more to offer. One direction which we now consider is not to use the light directly, but indirectly, by using the printed surface as reflector. The printed surface can be customized to the given environment, to its geometry and photometric properties, so as to optimize the distribution of light by a custommade reflector.
Another direction for future work is to go beyond the printability constraints, by assembling a large scale lampshade from surface pieces printed separately. This will allow using bigger and stronger light sources, and possibly compound arrays of light sources. A larger surface area will allow increasing the relative resolution of printed holes, and project an image on larger and more distant receiving surfaces. We are also considering focusing on the design of aesthetic imagery on the lampshade itself. Halftoned and stippled images have their own aesthetic virtue. In Figure 14, we show our lamp in a bedroom. Light effects combined with halftoning can lead to creative artistic media.
This work is supported in part by grants from National 973 Program (2015CB352501), by NSFC (61232011, 61202147, 61332015), and by the Israel Science Foundation (ISF).
References
 Alexa and Matusik (2012) Alexa, M. and Matusik, W. 2012. Irregular pit placement for dithering images by selfocclusion. Computers & Graphics 36, 6, 635 – 641. 2011 Joint Symposium on Computational Aesthetics (CAe), NonPhotorealistic Animation and Rendering (NPAR), and SketchBased Interfaces and Modeling (SBIM).
 Balzer et al. (2009) Balzer, M., Schlömer, T., and Deussen, O. 2009. Capacityconstrained point distributions: A variant of Lloyd’s method. ACM Trans. Graph. 28, 3 (July), 86:1–86:8.
 Baran et al. (2012) Baran, I., Keller, P., Bradley, D., Coros, S., Jarosz, W., Nowrouzezahrai, D., and Gross, M. 2012. Manufacturing layered attenuators for multiple prescribed shadow images. Comp. Graph. Forum 31, 2pt3 (May), 603–610.
 Chang et al. (2009) Chang, J., Alain, B., and Ostromoukhov, V. 2009. Structureaware error diffusion. ACM Trans. Graph. 28, 5 (Dec.), 162:1–162:8.
 de Goes et al. (2012) de Goes, F., Breeden, K., Ostromoukhov, V., and Desbrun, M. 2012. Blue noise through optimal transport. ACM Trans. Graph. 31, 6 (Nov.), 171:1–171:11.
 Finckh et al. (2010) Finckh, M., Dammertz, H., and Lensch, H. 2010. Geometry construction from caustic images. In Computer Vision  ECCV 2010, K. Daniilidis, P. Maragos, and N. Paragios, Eds. Lecture Notes in Computer Science, vol. 6315. Springer Berlin Heidelberg, 464–477.
 Kim et al. (2008) Kim, D., Son, M., Lee, Y., Kang, H., and Lee, S. 2008. Featureguided image stippling. In Proceedings of the Nineteenth Eurographics Conference on Rendering. EGSR’08. Eurographics Association, AirelaVille, Switzerland, Switzerland, 1209–1216.
 Kipphan (2001) Kipphan, H. 2001. Handbook of Print Media: Technologies and Production Methods. SpringerVerlag.
 Kiser et al. (2013) Kiser, T., Eigensatz, M., Nguyen, M., Bompas, P., and Pauly, M. 2013. Architectural caustics – controlling light with geometry. In Advances in Architectural Geometry 2012, L. Hesselgren, S. Sharma, J. Wallner, N. Baldassini, P. Bompas, and J. Raynaud, Eds. Springer Vienna, 91–106.
 Lan et al. (2013) Lan, Y., Dong, Y., Pellacini, F., and Tong, X. 2013. Biscale appearance fabrication. ACM Trans. Graph. 32, 4 (July), 145:1–145:12.
 Lau and Arce (2008) Lau, D. L. and Arce, G. R. 2008. Modern Digital Halftoning, Second Edition (Signal Processing and Communications), 2 ed. CRC Press.
 Levin et al. (2013) Levin, A., Glasner, D., Xiong, Y., Durand, F., Freeman, W., Matusik, W., and Zickler, T. 2013. Fabricating BRDFs at high spatial resolution using wave optics. ACM Trans. Graph. 32, 4 (July), 144:1–144:14.
 Li and Mould (2011) Li, H. and Mould, D. 2011. Structurepreserving stippling by prioritybased error diffusion. In Proceedings of Graphics Interface 2011. GI ’11. Canadian HumanComputer Communications Society, School of Computer Science, University of Waterloo, Waterloo, Ontario, Canada, 127–134.
 Li et al. (2010) Li, H., Wei, L.Y., Sander, P. V., and Fu, C.W. 2010. Anisotropic blue noise sampling. ACM Trans. Graph. 29, 6 (Dec.), 167:1–167:12.
 Lou and Stucki (1998) Lou, Q. and Stucki, P. 1998. Funfamentals of 3D halftoning. In Electronic Publishing, Artistic Imaging, and Digital Typography, 7th International Conference on Electronic Publishing, EP ’98, Held Jointly with the 4th International Conference on Raster Imaging and Digital Typography, RIDT ’98, St. Malo, France, March 30  April 3, 1998, Proceedings. 224–239.
 Malzbender et al. (2012) Malzbender, T., Samadani, R., Scher, S., Crume, A., Dunn, D., and Davis, J. 2012. Printing reflectance functions. ACM Trans. Graph. 31, 3 (June), 20:1–20:11.
 Mitra and Pauly (2009) Mitra, N. J. and Pauly, M. 2009. Shadow art. ACM Trans. Graph. 28, 5 (Dec.), 156:1–156:7.
 Pang et al. (2008) Pang, W.M., Qu, Y., Wong, T.T., CohenOr, D., and Heng, P.A. 2008. Structureaware halftoning. ACM Trans. Graph. 27, 3 (Aug.), 89:1–89:8.
 Papas et al. (2012) Papas, M., Houit, T., Nowrouzezahrai, D., Gross, M., and Jarosz, W. 2012. The magic lens: Refractive steganography. ACM Trans. Graph. 31, 6 (Nov.), 186:1–186:10.
 Papas et al. (2011) Papas, M., Jarosz, W., Jakob, W., Rusinkiewicz, S., Matusik, W., and Weyrich, T. 2011. Goalbased caustics. Comp. Graph. Forum 30, 2 (June), 503–511.
 Pereira et al. (2014) Pereira, T., Rusinkiewicz, S., and Matusik, W. 2014. Computational light routing: 3D printed optical fibers for sensing and display. ACM Trans. Graph. 33, 3 (June), 24:1–24:13.
 Schmaltz et al. (2010) Schmaltz, C., Gwosdek, P., Bruhn, A., and Weickert, J. 2010. Electrostatic halftoning. Comp. Graph. Forum 29, 8, 2313–2327.
 Schwartzburg et al. (2014) Schwartzburg, Y., Testuz, R., Tagliasacchi, A., and Pauly, M. 2014. Highcontrast computational caustic design. ACM Trans. Graph. 33, 4 (July), 74:1–74:11.
 Stucki (1997) Stucki, P. 1997. 3D halftoning. In Proc. SPIE 2949, Imaging Sciences and Display Technologies. Proc. SPIE 2949, 314–317.
 Ulichney (1988) Ulichney, R. 1988. Dithering with blue noise. Proceedings of the IEEE 76, 1 (Jan), 56–79.
 Wetzstein et al. (2011) Wetzstein, G., Lanman, D., Heidrich, W., and Raskar, R. 2011. Layered 3D: Tomographic image synthesis for attenuationbased light field and high dynamic range displays. ACM Trans. Graph. 30, 4 (July), 95:1–95:12.
 Weyrich et al. (2009) Weyrich, T., Peers, P., Matusik, W., and Rusinkiewicz, S. 2009. Fabricating microgeometry for custom surface reflectance. ACM Trans. Graph. 28, 3 (July), 32:1–32:6.
 Willis et al. (2012) Willis, K., Brockmeyer, E., Hudson, S., and Poupyrev, I. 2012. Printed optics: 3d printing of embedded optical elements for interactive devices. In Proceedings of the 25th Annual ACM Symposium on User Interface Software and Technology. UIST ’12. ACM, New York, NY, USA, 589–598.
 Yue et al. (2014) Yue, Y., Iwasaki, K., Chen, B.Y., Dobashi, Y., and Nishita, T. 2014. Poissonbased continuous surface generation for goalbased caustics. ACM Trans. Graph. 33, 3 (June), 31:1–31:7.
 Zhou and Chen (2009) Zhou, C. and Chen, Y. 2009. Threedimensional digital halftoning for layered manufacturing based on droplets. Transactions of North American Manufacturing Research Institute of SME 37, 175–182.
Comments
There are no comments yet.