使用 Python 图像库抗锯齿时，通过透明消除边缘的细边框

Question

我有一张高分辨率图像，想使用 Google Maps API v3 将其用作平铺地图叠加层。

我使用 MapTiler 将其分解为所需缩放级别的适当图块。效果很好，只是在作为原始图像边缘的图块上有一条细细的灰黑色边框。

遵循 http://groups.google.com 上的第二篇文章的建议/group/maptiler/browse_thread/thread/70a4c5610538332a/42fefedb4a0bc6d2，我尝试使用 gdal2tiles.py 代替，向其传递 -r antialias 选项，但细边框仍然存在。

如果我打开实际生成的图像图块，细边框确实看起来是生成图块的一部分，但它不是原始图像的一部分。

我怀疑正在发生的事情是，当程序生成图块图像文件时，我没有数据的 Google 地图图块部分被视为黑色像素，结果是灰色边框。

以下是我认为来自 gdal2tiles.py 的相关代码：

                    # Scaling by PIL (Python Imaging Library) - improved Lanczos
                    array = numpy.zeros((querysize, querysize, tilebands), numpy.uint8)
                    for i in range(tilebands):
                            array[:,:,i] = gdalarray.BandReadAsArray(dsquery.GetRasterBand(i+1), 0, 0, querysize, querysize)
                    im = Image.fromarray(array, 'RGBA') # Always four bands
                    im1 = im.resize((tilesize,tilesize), Image.ANTIALIAS)
                    if os.path.exists(tilefilename):
                            im0 = Image.open(tilefilename)
                            im1 = Image.composite(im1, im0, im1) 
                    im1.save(tilefilename,self.tiledriver)

知道如何制作它以使边框不存在，而不是在图像编辑器中打开相关生成的图块图像文件并将相关像素设置为透明吗？

我怀疑答案涉及找到某种方法来表示透明像素，以便抗锯齿功能在采样时忽略它们。

更新：它可能不会赢得任何优雅或性能奖项，但我认为我已经做到了这一点。当然，最后几码也是最艰难的。

如果我将 ANTIALIAS 更改为 BICUBIC，然后处理 alpha 通道，使任何半透明像素呈现完全透明，我就消除了大多数图块上的大部分边框。然而，一些浅色边界仍然存在。我不知道该怎么办。还值得注意的是，我猜想如果图像中存在不在实际图像区域边缘之外的透明或半透明像素，则此策略可能不会那么有效。

以下是经过这些修改的代码：

                # Scaling by PIL (Python Imaging Library) - improved Lanczos
                array = numpy.zeros((querysize, querysize, tilebands), numpy.uint8)
                for i in range(tilebands):
                        array[:,:,i] = gdalarray.BandReadAsArray(dsquery.GetRasterBand(i+1), 0, 0, querysize, querysize)
                im = Image.fromarray(array, 'RGBA') # Always four bands
                im1 = im.resize((tilesize,tilesize), Image.BICUBIC)
                if os.path.exists(tilefilename):
                        im0 = Image.open(tilefilename)
                        im1 = Image.composite(im1, im0, im1)
                im1AsArray = numpy.array(im1)
                alpha = im1AsArray[:,:,3]
                semiTransparentIndices = alpha < 255
                alpha[semiTransparentIndices] = 0
                im1AsArray[:,:,3] = alpha
                im1 = Image.fromarray(im1AsArray, 'RGBA')
                im1.save(tilefilename,self.tiledriver)

Answer 1

答案是将重采样更改为 BILINEAR（而不是 BICUBIC，这是我在发布到问题的更新中尝试的），然后确保将任何半透明像素更改为完全透明像素。

正如我在更新中所说，我在这里所做的代码修改可能不会赢得任何优雅或性能奖项，但它确实有效。以下是 gdal2tiles.py 中原始发布的代码片段需要更改为：

            # Scaling by PIL (Python Imaging Library) - improved Lanczos
            array = numpy.zeros((querysize, querysize, tilebands), numpy.uint8)
            for i in range(tilebands):
                    array[:,:,i] = gdalarray.BandReadAsArray(dsquery.GetRasterBand(i+1), 0, 0, querysize, querysize)
            im = Image.fromarray(array, 'RGBA') # Always four bands
            im1 = im.resize((tilesize,tilesize), Image.BILINEAR)
            if os.path.exists(tilefilename):
                    im0 = Image.open(tilefilename)
                    im1 = Image.composite(im1, im0, im1)
            im1AsArray = numpy.array(im1)
            alpha = im1AsArray[:,:,3]
            semiTransparentIndices = alpha < 255
            alpha[semiTransparentIndices] = 0
            im1AsArray[:,:,3] = alpha
            im1 = Image.fromarray(im1AsArray, 'RGBA')
            im1.save(tilefilename,self.tiledriver)

另请注意，仅当您向 gdal2tiles.py 传递 -r antialias 标志时，上述代码才会执行​​。是的，没错：我们更改了 -r antialias 代码，使其不进行抗锯齿。但如果您遇到了我所遇到的问题并且只是想要一个解决方案，那就是。