-
-
Save Flix01/254005101903fe1bc22fe0b9f3ed3931 to your computer and use it in GitHub Desktop.
| /* | |
| * edtaa3() | |
| * | |
| * Sweep-and-update Euclidean distance transform of an | |
| * image. Positive pixels are treated as object pixels, | |
| * zero or negative pixels are treated as background. | |
| * An attempt is made to treat antialiased edges correctly. | |
| * The input image must have pixels in the range [0,1], | |
| * and the antialiased image should be a box-filter | |
| * sampling of the ideal, crisp edge. | |
| * If the antialias region is more than 1 pixel wide, | |
| * the result from this transform will be inaccurate. | |
| * | |
| * By Stefan Gustavson ([email protected]). | |
| * | |
| * Originally written in 1994, based on a verbal | |
| * description of the SSED8 algorithm published in the | |
| * PhD dissertation of Ingemar Ragnemalm. This is his | |
| * algorithm, I only implemented it in C. | |
| * | |
| * Updated in 2004 to treat border pixels correctly, | |
| * and cleaned up the code to improve readability. | |
| * | |
| * Updated in 2009 to handle anti-aliased edges. | |
| * | |
| * Updated in 2011 to avoid a corner case infinite loop. | |
| * | |
| * Updated 2012 to change license from LGPL to MIT. | |
| * | |
| * Updated 2014 to fix a bug in the gradient computations. Ahem. | |
| * | |
| */ | |
| /* | |
| Copyright (C) 2009-2012 Stefan Gustavson ([email protected]) | |
| The code in this file is distributed under the MIT license: | |
| Permission is hereby granted, free of charge, to any person obtaining a copy | |
| of this software and associated documentation files (the "Software"), to deal | |
| in the Software without restriction, including without limitation the rights | |
| to use, copy, modify, merge, publish, distribute, sublicense, and/or sell | |
| copies of the Software, and to permit persons to whom the Software is | |
| furnished to do so, subject to the following conditions: | |
| The above copyright notice and this permission notice shall be included in | |
| all copies or substantial portions of the Software. | |
| THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR | |
| IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, | |
| FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE | |
| AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER | |
| LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM, | |
| OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN | |
| THE SOFTWARE. | |
| */ | |
| // Flix01: Added the EDTAA3FUNC_HAS_IMGUI_SUPPORT definition that can be used | |
| // to provide Dear Imgui (at https://github.com/ocornut/imgui) the method: | |
| // ImGui::PostBuildForSignedDistanceFontEffect(ImFontAtlas* atlas). | |
| // USAGE: | |
| /* | |
| // In imgui_impl_glfw.cpp (or similiar): | |
| #define IMIMPL_BUILD_SDF // Twickable | |
| #ifdef IMIMPL_BUILD_SDF | |
| #define EDTAA3FUNC_HAS_IMGUI_SUPPORT | |
| #include "edtaa3func.h" | |
| #endif //IMIMPL_BUILD_SDF | |
| bool ImGui_ImplGlfw_CreateFontsTexture() | |
| { | |
| // Build texture atlas | |
| ImGuiIO& io = ImGui::GetIO(); | |
| unsigned char* pixels; | |
| int width, height; | |
| // New code ------------------------------------------------- | |
| if (io.Fonts->ConfigData.empty()) io.Fonts->AddFontDefault(); | |
| io.Fonts->Build(); // (change this line if you use [imgui_freetype](https://github.com/Vuhdo/imgui) appropriately) | |
| # ifdef IMIMPL_BUILD_SDF | |
| ImGui::PostBuildForSignedDistanceFontEffect(io.Fonts); | |
| # endif //IMIMPL_BUILD_SDF | |
| // End new code ---------------------------------------------- | |
| ... | |
| } | |
| */ | |
| // The methods ImGuiEdtaaHelper::make_distance_mapd(...) and ImGuiEdtaaHelper::make_distance_mapb(...) | |
| // are based on https://github.com/rougier/freetype-gl/blob/master/distance-field.c [Copyright 2011,2012 Nicolas P. Rougier. All rights reserved.] | |
| // freetype-gl license: | |
| /* | |
| /* ========================================================================= | |
| * Freetype GL - A C OpenGL Freetype engine | |
| * Platform: Any | |
| * WWW: https://github.com/rougier/freetype-gl | |
| * ------------------------------------------------------------------------- | |
| * Copyright 2011,2012 Nicolas P. Rougier. All rights reserved. | |
| * | |
| * Redistribution and use in source and binary forms, with or without | |
| * modification, are permitted provided that the following conditions are met: | |
| * | |
| * 1. Redistributions of source code must retain the above copyright notice, | |
| * this list of conditions and the following disclaimer. | |
| * | |
| * 2. Redistributions in binary form must reproduce the above copyright | |
| * notice, this list of conditions and the following disclaimer in the | |
| * documentation and/or other materials provided with the distribution. | |
| * | |
| * THIS SOFTWARE IS PROVIDED BY NICOLAS P. ROUGIER ''AS IS'' AND ANY EXPRESS OR | |
| * IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF | |
| * MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO | |
| * EVENT SHALL NICOLAS P. ROUGIER OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, | |
| * INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES | |
| * (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES; | |
| * LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND | |
| * ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT | |
| * (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF | |
| * THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. | |
| * | |
| * The views and conclusions contained in the software and documentation are | |
| * those of the authors and should not be interpreted as representing official | |
| * policies, either expressed or implied, of Nicolas P. Rougier. | |
| * ========================================================================= | |
| */ | |
| #ifndef EDTAA3FUNC_H_ | |
| #define EDTAA3FUNC_H_ | |
| #pragma once | |
| #include <math.h> | |
| #ifdef __cplusplus | |
| extern "C" { | |
| #endif //__cplusplus | |
| /* | |
| * Compute the local gradient at edge pixels using convolution filters. | |
| * The gradient is computed only at edge pixels. At other places in the | |
| * image, it is never used, and it's mostly zero anyway. | |
| */ | |
| void computegradient(double *img, int w, int h, double *gx, double *gy) | |
| { | |
| int i,j,k,p,q; | |
| double glength, phi, phiscaled, ascaled, errsign, pfrac, qfrac, err0, err1, err; | |
| #define SQRT2 1.4142136 | |
| for(i = 1; i < h-1; i++) { // Avoid edges where the kernels would spill over | |
| for(j = 1; j < w-1; j++) { | |
| k = i*w + j; | |
| if((img[k]>0.0) && (img[k]<1.0)) { // Compute gradient for edge pixels only | |
| gx[k] = -img[k-w-1] - SQRT2*img[k-1] - img[k+w-1] + img[k-w+1] + SQRT2*img[k+1] + img[k+w+1]; | |
| gy[k] = -img[k-w-1] - SQRT2*img[k-w] - img[k-w+1] + img[k+w-1] + SQRT2*img[k+w] + img[k+w+1]; | |
| glength = gx[k]*gx[k] + gy[k]*gy[k]; | |
| if(glength > 0.0) { // Avoid division by zero | |
| glength = sqrt(glength); | |
| gx[k]=gx[k]/glength; | |
| gy[k]=gy[k]/glength; | |
| } | |
| } | |
| } | |
| } | |
| // TODO: Compute reasonable values for gx, gy also around the image edges. | |
| // (These are zero now, which reduces the accuracy for a 1-pixel wide region | |
| // around the image edge.) 2x2 kernels would be suitable for this. | |
| } | |
| /* | |
| * A somewhat tricky function to approximate the distance to an edge in a | |
| * certain pixel, with consideration to either the local gradient (gx,gy) | |
| * or the direction to the pixel (dx,dy) and the pixel greyscale value a. | |
| * The latter alternative, using (dx,dy), is the metric used by edtaa2(). | |
| * Using a local estimate of the edge gradient (gx,gy) yields much better | |
| * accuracy at and near edges, and reduces the error even at distant pixels | |
| * provided that the gradient direction is accurately estimated. | |
| */ | |
| double edgedf(double gx, double gy, double a) | |
| { | |
| double df, glength, temp, a1; | |
| if ((gx == 0) || (gy == 0)) { // Either A) gu or gv are zero, or B) both | |
| df = 0.5-a; // Linear approximation is A) correct or B) a fair guess | |
| } else { | |
| glength = sqrt(gx*gx + gy*gy); | |
| if(glength>0) { | |
| gx = gx/glength; | |
| gy = gy/glength; | |
| } | |
| /* Everything is symmetric wrt sign and transposition, | |
| * so move to first octant (gx>=0, gy>=0, gx>=gy) to | |
| * avoid handling all possible edge directions. | |
| */ | |
| gx = fabs(gx); | |
| gy = fabs(gy); | |
| if(gx<gy) { | |
| temp = gx; | |
| gx = gy; | |
| gy = temp; | |
| } | |
| a1 = 0.5*gy/gx; | |
| if (a < a1) { // 0 <= a < a1 | |
| df = 0.5*(gx + gy) - sqrt(2.0*gx*gy*a); | |
| } else if (a < (1.0-a1)) { // a1 <= a <= 1-a1 | |
| df = (0.5-a)*gx; | |
| } else { // 1-a1 < a <= 1 | |
| df = -0.5*(gx + gy) + sqrt(2.0*gx*gy*(1.0-a)); | |
| } | |
| } | |
| return df; | |
| } | |
| double distaa3(double *img, double *gximg, double *gyimg, int w, int c, int xc, int yc, int xi, int yi) | |
| { | |
| double di, df, dx, dy, gx, gy, a; | |
| int closest; | |
| closest = c-xc-yc*w; // Index to the edge pixel pointed to from c | |
| a = img[closest]; // Grayscale value at the edge pixel | |
| gx = gximg[closest]; // X gradient component at the edge pixel | |
| gy = gyimg[closest]; // Y gradient component at the edge pixel | |
| if(a > 1.0) a = 1.0; | |
| if(a < 0.0) a = 0.0; // Clip grayscale values outside the range [0,1] | |
| if(a == 0.0) return 1000000.0; // Not an object pixel, return "very far" ("don't know yet") | |
| dx = (double)xi; | |
| dy = (double)yi; | |
| di = sqrt(dx*dx + dy*dy); // Length of integer vector, like a traditional EDT | |
| if(di==0) { // Use local gradient only at edges | |
| // Estimate based on local gradient only | |
| df = edgedf(gx, gy, a); | |
| } else { | |
| // Estimate gradient based on direction to edge (accurate for large di) | |
| df = edgedf(dx, dy, a); | |
| } | |
| return di + df; // Same metric as edtaa2, except at edges (where di=0) | |
| } | |
| // Shorthand macro: add ubiquitous parameters dist, gx, gy, img and w and call distaa3() | |
| #define DISTAA(c,xc,yc,xi,yi) (distaa3(img, gx, gy, w, c, xc, yc, xi, yi)) | |
| void edtaa3(double *img, double *gx, double *gy, int w, int h, short *distx, short *disty, double *dist) | |
| { | |
| int x, y, i, c; | |
| int offset_u, offset_ur, offset_r, offset_rd, | |
| offset_d, offset_dl, offset_l, offset_lu; | |
| double olddist, newdist; | |
| int cdistx, cdisty, newdistx, newdisty; | |
| int changed; | |
| double epsilon = 1e-3; | |
| /* Initialize index offsets for the current image width */ | |
| offset_u = -w; | |
| offset_ur = -w+1; | |
| offset_r = 1; | |
| offset_rd = w+1; | |
| offset_d = w; | |
| offset_dl = w-1; | |
| offset_l = -1; | |
| offset_lu = -w-1; | |
| /* Initialize the distance images */ | |
| for(i=0; i<w*h; i++) { | |
| distx[i] = 0; // At first, all pixels point to | |
| disty[i] = 0; // themselves as the closest known. | |
| if(img[i] <= 0.0) | |
| { | |
| dist[i]= 1000000.0; // Big value, means "not set yet" | |
| } | |
| else if (img[i]<1.0) { | |
| dist[i] = edgedf(gx[i], gy[i], img[i]); // Gradient-assisted estimate | |
| } | |
| else { | |
| dist[i]= 0.0; // Inside the object | |
| } | |
| } | |
| /* Perform the transformation */ | |
| do | |
| { | |
| changed = 0; | |
| /* Scan rows, except first row */ | |
| for(y=1; y<h; y++) | |
| { | |
| /* move index to leftmost pixel of current row */ | |
| i = y*w; | |
| /* scan right, propagate distances from above & left */ | |
| /* Leftmost pixel is special, has no left neighbors */ | |
| olddist = dist[i]; | |
| if(olddist > 0) // If non-zero distance or not set yet | |
| { | |
| c = i + offset_u; // Index of candidate for testing | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx; | |
| newdisty = cdisty+1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_ur; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx-1; | |
| newdisty = cdisty+1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| changed = 1; | |
| } | |
| } | |
| i++; | |
| /* Middle pixels have all neighbors */ | |
| for(x=1; x<w-1; x++, i++) | |
| { | |
| olddist = dist[i]; | |
| if(olddist <= 0) continue; // No need to update further | |
| c = i+offset_l; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx+1; | |
| newdisty = cdisty; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_lu; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx+1; | |
| newdisty = cdisty+1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_u; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx; | |
| newdisty = cdisty+1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_ur; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx-1; | |
| newdisty = cdisty+1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| changed = 1; | |
| } | |
| } | |
| /* Rightmost pixel of row is special, has no right neighbors */ | |
| olddist = dist[i]; | |
| if(olddist > 0) // If not already zero distance | |
| { | |
| c = i+offset_l; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx+1; | |
| newdisty = cdisty; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_lu; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx+1; | |
| newdisty = cdisty+1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_u; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx; | |
| newdisty = cdisty+1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| changed = 1; | |
| } | |
| } | |
| /* Move index to second rightmost pixel of current row. */ | |
| /* Rightmost pixel is skipped, it has no right neighbor. */ | |
| i = y*w + w-2; | |
| /* scan left, propagate distance from right */ | |
| for(x=w-2; x>=0; x--, i--) | |
| { | |
| olddist = dist[i]; | |
| if(olddist <= 0) continue; // Already zero distance | |
| c = i+offset_r; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx-1; | |
| newdisty = cdisty; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| changed = 1; | |
| } | |
| } | |
| } | |
| /* Scan rows in reverse order, except last row */ | |
| for(y=h-2; y>=0; y--) | |
| { | |
| /* move index to rightmost pixel of current row */ | |
| i = y*w + w-1; | |
| /* Scan left, propagate distances from below & right */ | |
| /* Rightmost pixel is special, has no right neighbors */ | |
| olddist = dist[i]; | |
| if(olddist > 0) // If not already zero distance | |
| { | |
| c = i+offset_d; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx; | |
| newdisty = cdisty-1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_dl; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx+1; | |
| newdisty = cdisty-1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| changed = 1; | |
| } | |
| } | |
| i--; | |
| /* Middle pixels have all neighbors */ | |
| for(x=w-2; x>0; x--, i--) | |
| { | |
| olddist = dist[i]; | |
| if(olddist <= 0) continue; // Already zero distance | |
| c = i+offset_r; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx-1; | |
| newdisty = cdisty; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_rd; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx-1; | |
| newdisty = cdisty-1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_d; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx; | |
| newdisty = cdisty-1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_dl; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx+1; | |
| newdisty = cdisty-1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| changed = 1; | |
| } | |
| } | |
| /* Leftmost pixel is special, has no left neighbors */ | |
| olddist = dist[i]; | |
| if(olddist > 0) // If not already zero distance | |
| { | |
| c = i+offset_r; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx-1; | |
| newdisty = cdisty; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_rd; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx-1; | |
| newdisty = cdisty-1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| olddist=newdist; | |
| changed = 1; | |
| } | |
| c = i+offset_d; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx; | |
| newdisty = cdisty-1; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| changed = 1; | |
| } | |
| } | |
| /* Move index to second leftmost pixel of current row. */ | |
| /* Leftmost pixel is skipped, it has no left neighbor. */ | |
| i = y*w + 1; | |
| for(x=1; x<w; x++, i++) | |
| { | |
| /* scan right, propagate distance from left */ | |
| olddist = dist[i]; | |
| if(olddist <= 0) continue; // Already zero distance | |
| c = i+offset_l; | |
| cdistx = distx[c]; | |
| cdisty = disty[c]; | |
| newdistx = cdistx+1; | |
| newdisty = cdisty; | |
| newdist = DISTAA(c, cdistx, cdisty, newdistx, newdisty); | |
| if(newdist < olddist-epsilon) | |
| { | |
| distx[i]=newdistx; | |
| disty[i]=newdisty; | |
| dist[i]=newdist; | |
| changed = 1; | |
| } | |
| } | |
| } | |
| } | |
| while(changed); // Sweep until no more updates are made | |
| /* The transformation is completed. */ | |
| } | |
| #ifdef __cplusplus | |
| } // extern "C" | |
| #endif //__cplusplus | |
| #ifdef EDTAA3FUNC_HAS_IMGUI_SUPPORT | |
| #include <imgui.h> | |
| namespace ImGuiEdtaaHelper { | |
| // The purpose of TmpData is to allow releasing the memory all togther | |
| // at the end of a multiple call to make_distance_mapb(...)/make_distance_mapd(...) | |
| // [However we don't use this feature here] | |
| struct TmpData { | |
| ImVector<double> v_data; | |
| ImVector<short> v_xdist; | |
| ImVector<short> v_ydist; | |
| ImVector<double> v_gx; | |
| ImVector<double> v_gy; | |
| ImVector<double> v_outside; | |
| ImVector<double> v_inside; | |
| void freeMemory() { | |
| TmpData tmp; | |
| v_data.swap(tmp.v_data); | |
| v_xdist.swap(tmp.v_xdist); | |
| v_ydist.swap(tmp.v_ydist); | |
| v_gx.swap(tmp.v_gx); | |
| v_outside.swap(tmp.v_outside); | |
| v_inside.swap(tmp.v_inside); | |
| } | |
| }; | |
| static TmpData tmpData; | |
| static void FreeTmpMemory() {tmpData.freeMemory();} | |
| // These two methods are adapted from: https://github.com/rougier/freetype-gl (Copyright 2011,2012 Nicolas P. Rougier), | |
| // so that the memory they use can optionally be freed once, after all the glyphs are processed. | |
| static double * make_distance_mapd( double *data, unsigned int width, unsigned int height, TmpData* pOptionalTmpData=NULL) | |
| { | |
| if (!pOptionalTmpData) pOptionalTmpData=&tmpData; | |
| const unsigned int area = width * height; | |
| IM_ASSERT(area>0); | |
| pOptionalTmpData->v_xdist.resize(area); short* xdist = &pOptionalTmpData->v_xdist[0]; | |
| pOptionalTmpData->v_ydist.resize(area); short* ydist = &pOptionalTmpData->v_ydist[0]; | |
| pOptionalTmpData->v_gx.resize(area); double* gx = &pOptionalTmpData->v_gx[0]; | |
| pOptionalTmpData->v_gy.resize(area); double* gy = &pOptionalTmpData->v_gy[0]; | |
| pOptionalTmpData->v_outside.resize(area); double* outside=&pOptionalTmpData->v_outside[0]; | |
| pOptionalTmpData->v_inside.resize(area); double* inside= &pOptionalTmpData->v_inside[0]; | |
| double vmin = DBL_MAX; | |
| unsigned int i; | |
| // Compute outside = edtaa3(bitmap); % Transform background (0's) | |
| computegradient( data, width, height, gx, gy); | |
| edtaa3(data, gx, gy, width, height, xdist, ydist, outside); | |
| for( i=0; i<area; ++i) { | |
| if( outside[i] < 0.0 ) outside[i] = 0.0; | |
| } | |
| // Compute inside = edtaa3(1-bitmap); % Transform foreground (1's) | |
| memset( gx, 0, sizeof(double)*area ); | |
| memset( gy, 0, sizeof(double)*area ); | |
| for( i=0; i<area; ++i) data[i] = 1 - data[i]; | |
| computegradient( data, width, height, gx, gy ); | |
| edtaa3( data, gx, gy, width, height, xdist, ydist, inside ); | |
| for( i=0; i<area; ++i ) { | |
| if( inside[i] < 0 ) inside[i] = 0.0; | |
| } | |
| // distmap = outside - inside; % Bipolar distance field | |
| for( i=0; i<area; ++i) { | |
| outside[i] -= inside[i]; | |
| if( outside[i] < vmin ) vmin = outside[i]; | |
| } | |
| vmin = fabs(vmin); | |
| for( i=0; i<area; ++i) { | |
| double v = outside[i]; | |
| if ( v < -vmin) outside[i] = -vmin; | |
| else if( v > +vmin) outside[i] = +vmin; | |
| data[i] = (outside[i]+vmin)/(2*vmin); | |
| } | |
| return data; | |
| } | |
| static unsigned char * make_distance_mapb(const unsigned char *img,unsigned int width, unsigned int height, unsigned char* imgOutSameSize, TmpData* pOptionalTmpData=NULL) | |
| { | |
| IM_ASSERT(imgOutSameSize); | |
| if (!pOptionalTmpData) pOptionalTmpData=&tmpData; | |
| const unsigned int area = width * height; | |
| IM_ASSERT(area>0); | |
| pOptionalTmpData->v_data.resize(area); | |
| double* data = &pOptionalTmpData->v_data[0]; | |
| unsigned int i; | |
| // find minimimum and maximum values | |
| double img_min = DBL_MAX,img_max = DBL_MIN,v; | |
| for( i=0; i<area; ++i) { | |
| data[i] = v = img[i]; | |
| if (v > img_max) img_max = v; | |
| if (v < img_min) img_min = v; | |
| } | |
| // Map values from 0 - 255 to 0.0 - 1.0 | |
| for( i=0; i<area; ++i) data[i] = (img[i]-img_min)/img_max; | |
| data = make_distance_mapd(data, width, height,pOptionalTmpData); | |
| // map values from 0.0 - 1.0 to 0 - 255 | |
| for( i=0; i<area; ++i) imgOutSameSize[i] = (unsigned char)(255*(1-data[i])); | |
| return imgOutSameSize; | |
| } | |
| } //ImGuiEdtaaHelper | |
| namespace ImGui { | |
| void ProcessAlpha8ImageForSignedDistanceFontEffect(unsigned char* pixels,int w,int h) { | |
| IM_ASSERT(pixels && w>0 && h>0); | |
| ImVector<unsigned char> tmpData;tmpData.resize(w*h); | |
| ImGuiEdtaaHelper::make_distance_mapb(pixels,w,h,&tmpData[0]); | |
| ImGuiEdtaaHelper::FreeTmpMemory(); | |
| memcpy(pixels,&tmpData[0],w*h); | |
| } | |
| void PostBuildForSignedDistanceFontEffect(ImFontAtlas* atlas) { | |
| if (!atlas->TexPixelsAlpha8) atlas->GetTexDataAsAlpha8(&atlas->TexPixelsAlpha8,NULL,NULL); | |
| ProcessAlpha8ImageForSignedDistanceFontEffect(atlas->TexPixelsAlpha8,atlas->TexWidth,atlas->TexHeight); | |
| } | |
| } | |
| #endif //EDTAA3FUNC_HAS_IMGUI_SUPPORT | |
| #endif //EDTAA3FUNC_H_ |
Making further experiments, it seems that the edtaa approach used by freetype_gl is not only much easier, but much better: see some results:
Ubuntu-B.ttf 13 pxls cfg.OversampleH=cfg.OversampleV=1; (16 times less texture memory than in the pic posted above!)

And the default font now works much better too:

All these pics have been taken using the standard stb_truetype.h. Chances are that using the freetype lib we can get even better results.
ROADMAP:
- Although it's trivial, flush all this code and post another one using edtaa.
Some reference code notes for the impatients:
// At the end of a standard Build(ImFontAtlas* atlas) we just need to append something like:
int w,h;
unsigned char* texData;
atlas->GetTexDataAsAlpha8(&texData,&w,&h);
ImVector<unsigned char> tmpData;tmpData.resize(w*h);
ImGuiEdtaaHelper::make_distance_mapb(texData,w,h,&tmpData[0]);
memcpy(texData,&tmpData[0],w*h);
// where make_distance_mapb(...) is a trivial adaptation of the code that uses edtaa3func.c and that can be found here:
// https://github.com/rougier/freetype-gl/blob/master/distance-field.cDONE in Revision 2 😃!
And now it's independent on the "font library" stb_truetype.h (and can be used in imgui_freetype too).
Revision 3
Same as Rev.2, but all embedded in the single file: edtaa3func.h.
See the usage documentation there (however the fragment shader code has been stripped from the comment: you may need to download Rev.2 if you need it).
How do you implement font stroke?
@NewHand007: the fragment shader code has been stripped from the comment: you may need to download Rev.2 if you need it.
sorry ,where can i download "Rev.2" ?
I've been working on font strokes for a few days. Can you share a copy of the code that works my bro?
@NewHand007: the fragment shader code has been stripped from the comment: you may need to download Rev.2 if you need it.
bro
The Rev.2 file with the comment is here
The Rev.2 file with the comment is here
thank you bro !
Docs are at the top of imguiSdfGenerator.h.

Ubuntu-B.ttf 13 pxls cfg.OversampleH=cfg.OversampleV=4;
[Edit:]: In Rev.1 I used direct calculations on the vectorial shape of the glyphs.
Another (more common) approach, easier to implement, is to simply post-process the normal glyph texture using the file edtaa by Stefan Gustavson. This is used in a lot of projects, such as the popular freetype_gl.