diff --git a/package.json b/package.json
index 4dc7f994f..1cdbaab35 100644
--- a/package.json
+++ b/package.json
@@ -54,8 +54,7 @@
         "d3-geo": "^3.0.1",
         "d3-random": "^3.0.1",
         "d3-selection": "^3.0.0",
-        "d3-zoom": "^3.0.0",
-        "fast-kde": "^0.2.1"
+        "d3-zoom": "^3.0.0"
     },
     "devDependencies": {
         "@babel/core": "^7.9.0",
diff --git a/src/index.js b/src/index.js
index 958b61dce..71c2a395f 100644
--- a/src/index.js
+++ b/src/index.js
@@ -22,9 +22,6 @@ export { JoyPlotStyle } from './style/JoyPlotStyle.js'
 export { CompositionStyle } from './style/CompositionStyle.js'
 export { SegmentStyle } from './style/SegmentStyle.js'
 export { TextStyle } from './style/TextStyle.js'
-export { KernelSmoothingStyle } from './style/KernelSmoothingStyle.js'
-//export { KernelSmoothingOldStyle } from "./style/KernelSmoothingStyle"
-//export { KernelSmoothingWGLStyle } from "./style/KernelSmoothingStyle"
 export { PillarStyle } from './style/PillarStyle.js'
 export { SideStyle } from './style/SideStyle.js'
 export { ContourStyle } from './style/ContourStyle.js'
diff --git a/src/style/KernelSmoothingOldStyle.js b/src/style/KernelSmoothingOldStyle.js
deleted file mode 100644
index 1de177e4f..000000000
--- a/src/style/KernelSmoothingOldStyle.js
+++ /dev/null
@@ -1,221 +0,0 @@
-//@ts-check
-'use strict'
-
-import { Style } from '../Style.js'
-
-/**
- * A style representing the cell as a smoothed layer, to smoothing local variations and show main trends across space.
- *
- * @author Julien Gaffuri
- */
-export class KernelSmoothingStyle extends Style {
-    //bandwidth
-    // https://observablehq.com/@uwdata/fast-kde
-    // https://observablehq.com/d/3127b6d89ada959f
-    //TODO use https://github.com/uwdata/fast-kde ?
-    //TODO https://observablehq.com/@sahilchinoy/areal-interpolation-iii ?
-
-    /** @param {object} opts */
-    constructor(opts) {
-        super(opts)
-        opts = opts || {}
-
-        /** A function specifying the value to consider for each cell. This is the value to smooth.
-         * @type {function(import("../Dataset").Cell):number} */
-        this.value = opts.value
-
-        /** The smoothing parameter, in geo unit. The larger, the more smoothed.
-         * @type {function(number,number):number}
-         */
-        this.sigma = opts.sigma // (r, zf)=>...
-
-        /** The styles to represent the smoothed cells.
-         * @type {Array.<Style>}
-         */
-        this.styles = opts.styles
-    }
-
-    /**
-     * Compute the kernel matrix, that is the matrix of the weights.
-     * One quadrant is necessary only, since it is symetrical (along both x and y axes).
-     * @param {number} s Sigma (in grid pixel !)
-     * @returns {Array.<Array<number>>}
-     * @private
-     */
-    getKernelMatrix(s) {
-        //the size of the kernel: lets limit that to ~3 times the standard deviation, as an approximation.
-        const kernelSize = Math.floor(3.25 * s)
-
-        //pre-compute coefficient for gaussian computation, to avoid computing them every time.
-        const c2 = 2 * s * s
-
-        /**
-         * The gaussian function. TODO expose that function as a parameter, to use other kernels ?
-         * @param {number} x
-         * @param {number} y
-         * @returns {number}
-         */
-        const gaussian = (x, y) => (c2 != 0 ? Math.exp(-(x * x + y * y) / c2) : x == 0 && y == 0 ? 1 : 0) //dirac case
-
-        /** @type {Array.<Array<number>>} */
-        const kw = []
-        for (let wi = 0; wi <= kernelSize; wi++) {
-            /** @type {Array<number>} */
-            const col = []
-            for (let wj = 0; wj <= kernelSize; wj++) {
-                //compute weight at wi,wj
-                col.push(gaussian(wi, wj))
-            }
-            kw.push(col)
-        }
-
-        /*
-        //debug: show values
-        for (let wj = kernelSize-1; wj >=0; wj--) {
-            let st = "";
-            for (let wi = 0; wi <= kernelSize; wi++) {
-                st += "   " + Math.floor(kw[wi][wj].w *1000) /1000
-            }
-            console.log(st)
-        }*/
-
-        return kw
-    }
-
-    /**
-     * Compute kernel smoothing of some cells.
-     *
-     * @private
-     * @param {Array.<import("../Dataset").Cell>} cells The cells to be smoothed.
-     * @param {import("../Dataset").Envelope} e Geo envelope to consider.
-     * @param {number} r Resolution, in geo unit.
-     * @param {number} s Sigma (in grid pixel !)
-     * @returns {Array.<import("../Dataset").Cell>} The list of cells, including the initial ones and the ones with smoothed values, in "ksmval" property.
-     */
-    kernelSmoothing(cells, e, r, s) {
-        //filter
-        if (this.filter) cells = cells.filter(this.filter)
-
-        const wThr = 0.01 //gaussian weights below this value will be ignored
-
-        //compute extent, in grid position
-        const xMin = Math.floor(e.xMin / r) * r
-        const xMax = Math.floor(e.xMax / r) * r
-        const yMin = Math.floor(e.yMin / r) * r
-        const yMax = Math.floor(e.yMax / r) * r
-
-        //compute matrix dimensions
-        const nbX = (xMax - xMin) / r + 1
-        const nbY = (yMax - yMin) / r + 1
-
-        //initialise smoothed value to 0
-        for (const c of cells) c['ksmval'] = 0
-
-        //make matrix. add input cells by i/j (grid position)
-        const index = {}
-        for (const c of cells) {
-            // i,j of the cell
-            const i = Math.floor((c.x - xMin) / r)
-            const j = Math.floor((c.y - yMin) / r)
-            if (!index[i]) index[i] = {}
-            index[i][j] = c
-        }
-
-        /** Check that the cell i,j is within the frame */
-        const isWithinFrame = (i, j) => i >= 0 && i < nbX && j >= 0 && j < nbY
-
-        //get kernel matrix
-        /** @type {Array.<Array.<number>>} */
-        const kernelMatrix = this.getKernelMatrix(s)
-        const kernelSize = kernelMatrix.length
-
-        //compute summ of the weights over entire kernel window
-        /** @type {number} */
-        let sumWeights = 0
-        for (let ki = 1 - kernelSize; ki < kernelSize; ki++) {
-            for (let kj = 1 - kernelSize; kj < kernelSize; kj++) {
-                let w = kernelMatrix[Math.abs(ki)][Math.abs(kj)]
-                if (w < wThr) continue
-                sumWeights += w
-            }
-        }
-
-        /** Add v as a contribution to the cell i,j */
-        const addContributionTo = (i, j, v) => {
-            //get cell at (i,j)
-            const c_ = index[i] ? index[i][j] : undefined
-
-            if (c_) {
-                //cell exists: add contribution
-                if (c_.ksmval) c_.ksmval += v
-                else c_.ksmval = v
-            } else {
-                //cell does not exist: create a new one with the smoothed value
-                if (!index[i]) index[i] = {}
-                index[i][j] = { x: xMin + i * r, y: yMin + j * r, ksmval: v }
-            }
-        }
-
-        //compute smoothing, input cell by input cell
-        for (const c of cells) {
-            /** get value of cell c
-             * @type {number} */
-            const val = this.value(c)
-            if (!val) continue
-
-            //cell matrix coordinates
-            const i = Math.floor((c.x - xMin) / r)
-            const j = Math.floor((c.y - yMin) / r)
-
-            //add contributions to smoothed values
-            for (let ki = 1 - kernelSize; ki < kernelSize; ki++) {
-                for (let kj = 1 - kernelSize; kj < kernelSize; kj++) {
-                    //check cell is within the frame
-                    if (!isWithinFrame(i + ki, j + kj)) continue
-
-                    //get contribution (ki,kj)
-                    let w = kernelMatrix[Math.abs(ki)][Math.abs(kj)]
-                    if (!w || w < wThr) continue
-                    let v = w * val
-                    if (!v) continue
-                    v /= sumWeights
-                    if (!v) continue
-
-                    //add contribution
-                    addContributionTo(i + ki, j + kj, v)
-                }
-            }
-        }
-
-        //make output list
-        const out = []
-        for (let i of Object.keys(index)) for (const j of Object.keys(index[i])) out.push(index[i][j])
-
-        return out
-    }
-
-    /**
-     * Draw the smoothed cells depending on the list of styles specified.
-     *
-     * @param {Array.<import("../Dataset").Cell>} cells
-     * @param {number} r
-     * @param {import("../GeoCanvas").GeoCanvas} cg
-     */
-    draw(cells, r, cg) {
-        if (!cells || cells.length == 0) return
-
-        //get smoothing param in geo unit
-        /** @type {number} */
-        const sG = this.sigma(r, cg.zf)
-
-        //apply kernel smoothing
-        cells = this.kernelSmoothing(cells, cg.extGeo, r, sG / r)
-
-        //draw smoothed cells from styles
-        for (let s of this.styles) s.draw(cells, r, cg)
-
-        //update legends
-        //for (let s of this.styles)
-        //    s.updateLegends({ style: s, r: r, zf: cg.getZf() });
-    }
-}
diff --git a/src/style/KernelSmoothingStyle.js b/src/style/KernelSmoothingStyle.js
deleted file mode 100644
index e3d7cfbe2..000000000
--- a/src/style/KernelSmoothingStyle.js
+++ /dev/null
@@ -1,115 +0,0 @@
-//@ts-check
-'use strict'
-
-import { Style } from '../Style.js'
-import { density2d } from 'fast-kde'
-
-/**
- * A style representing the cell as a smoothed layer, to smoothing local variations and show main trends across space.
- *
- * @author Julien Gaffuri
- */
-export class KernelSmoothingStyle extends Style {
-    // https://observablehq.com/d/5dd1cb5e4d21c021
-    // https://observablehq.com/@uwdata/fast-kde
-    // https://observablehq.com/d/3127b6d89ada959f
-    //TODO https://observablehq.com/@sahilchinoy/areal-interpolation-iii ?
-
-    /** @param {object} opts */
-    constructor(opts) {
-        super(opts)
-        opts = opts || {}
-
-        /** A function returning the value to consider for each cell. This is the value to be smoothed.
-         * @type {function(import("../Dataset").Cell):number} */
-        this.value = opts.value
-
-        /** The smoothing parameter, in geo unit. The larger, the more smoothed.
-         * @type {function(number,number):number}
-         */
-        this.sigma = opts.sigma // (r, zf)=>...
-
-        /** A factor to adjust the smoothed grid resolution.
-         * When set to 1, the smoothed grid is exactly the screen resolution.
-         * Set to 2 to degrade the resolution to a factor 2.
-         * The higher, the more pixelised and the faster to compute.
-         * @type { number }
-         */
-        this.factor = opts.factor || 2
-
-        /** A filter function to filter the smoothed cells based on their smoothed value.
-         *  Return true to keep the cell, false otherwise.
-         * @type { function(number):boolean }
-         */
-        this.filterSm = opts.filterSm
-
-        /** The name of the attribute where the smoothed value is stored in the output smoothed grid.
-         * @type { string }
-         */
-        this.sCol = opts.sCol || 'ksmval'
-
-        /** The styles to represent the smoothed grid.
-         * @type {Array.<Style>}
-         */
-        this.styles = opts.styles || []
-    }
-
-    /**
-     * Draw the smoothed cells depending on the list of styles specified.
-     *
-     * @param {Array.<import("../Dataset").Cell>} cells
-     * @param {number} r
-     * @param {import("../GeoCanvas").GeoCanvas} cg
-     */
-    draw(cells, r, cg) {
-        //filter
-        if (this.filter) cells = cells.filter(this.filter)
-
-        if (!cells || cells.length == 0) return
-
-        //get smoothing param in geo unit
-        /** @type {number} */
-        const sG = this.sigma(r, cg.zf)
-
-        //compute smoothed grid dimensions
-        const nbX = Math.ceil(cg.w / this.factor)
-        const nbY = Math.ceil(cg.h / this.factor)
-        //compute smoothed grid geo extent
-        const e_ = [
-            [cg.pixToGeoX(0), cg.pixToGeoX(nbX * this.factor)],
-            [cg.pixToGeoY(nbY * this.factor), cg.pixToGeoY(0)],
-        ]
-
-        //compute smoothed grid
-        let g = density2d(cells, {
-            x: (c) => c.x + r / 2,
-            y: (c) => c.y + r / 2,
-            weight: (c) => this.value(c),
-            bins: [nbX, nbY],
-            bandwidth: sG,
-            extent: e_,
-        }).grid()
-
-        //compute the resolution of the smoothed grid
-        const resSmoothed = (e_[0][1] - e_[0][0]) / nbX
-
-        //make smoothed cells
-        cells = []
-        for (let ind = 0; ind < g.length; ind++) {
-            const v = g[ind]
-            if (this.filterSm && !this.filterSm(v)) continue
-            const row = Math.floor(ind / nbX)
-            const col = ind % nbX
-            const c = { x: e_[0][0] + col * resSmoothed, y: e_[1][0] + row * resSmoothed }
-            c[this.sCol] = v
-            cells.push(c)
-        }
-
-        //draw smoothed cells from styles
-        for (let s of this.styles) s.draw(cells, resSmoothed, cg)
-
-        //update legends
-        //for (let s of this.styles)
-        //    s.updateLegends({ style: s, r: r, zf: cg.getZf() });
-    }
-}
diff --git a/src/style/KernelSmoothingWGLStyle.js b/src/style/KernelSmoothingWGLStyle.js
deleted file mode 100644
index d4bf8c837..000000000
--- a/src/style/KernelSmoothingWGLStyle.js
+++ /dev/null
@@ -1,130 +0,0 @@
-//@ts-check
-'use strict'
-
-import { Style } from '../Style.js'
-import { makeWebGLCanvas } from '../utils/webGLUtils'
-import { WebGLSquareColoringKS } from '../utils/WebGLSquareColoringKS'
-
-/**
- * @author Julien Gaffuri
- */
-export class KernelSmoothingStyle extends Style {
-    //see
-    //https://stackoverflow.com/questions/8099979/creating-a-glsl-arrays-of-uniforms
-    //see https://gist.github.com/jasonkit/c5b4fd62e8cbfe2780cc
-
-    /** @param {object} opts */
-    constructor(opts) {
-        super(opts)
-        opts = opts || {}
-
-        /**
-         * The name of the column/attribute of the tabular data where to retrieve the variable for color.
-         * @type {string} */
-        this.colorCol = opts.colorCol
-
-        /**
-         * A function returning the t value (within [0,1]) of the cell.
-         * @type {function(number,number,import("../Style").Stat):number} */
-        this.tFun = opts.tFun || ((v, r, s) => v / s.max)
-
-        /**
-         * Distribution stretching method.
-         * The stretching is performed on GPU side (fragment shader).
-         * @type {{ fun:string, alpha:number }} */
-        this.stretching = opts.stretching
-
-        /**
-         * The sample of the color ramp.
-         * The color is computed on GPU side (fragment shader) based on those values (linear interpolation).
-         * @type {Array.<string>} */
-        this.colors =
-            opts.colors ||
-            [
-                'rgb(158, 1, 66)',
-                'rgb(248, 142, 83)',
-                'rgb(251, 248, 176)',
-                'rgb(137, 207, 165)',
-                'rgb(94, 79, 162)',
-            ].reverse()
-        if (opts.color)
-            this.colors = [
-                opts.color(0),
-                opts.color(0.2),
-                opts.color(0.4),
-                opts.color(0.6),
-                opts.color(0.8),
-                opts.color(1),
-            ]
-
-        /**
-         * A function returning the size of the cells, in geographical unit. All cells have the same size.
-         * @type {function(number,number):number} */
-        this.size = opts.size // (resolution, zf) => ...
-    }
-
-    /**
-     * @param {Array.<import("../Dataset").Cell>} cells
-     * @param {number} r
-     * @param {import("../GeoCanvas").GeoCanvas} cg
-     */
-    draw(cells, r, cg) {
-        //filter
-        if (this.filter) cells = cells.filter(this.filter)
-
-        //zoom factor
-        const zf = cg.getZf()
-
-        //compute color variable statistics
-        const statColor = Style.getStatistics(cells, (c) => c[this.colorCol], true)
-
-        if (!statColor) return
-
-        //create canvas and webgl renderer
-        const cvWGL = makeWebGLCanvas(cg.w + '', cg.h + '')
-        if (!cvWGL) {
-            console.error('No webGL')
-            return
-        }
-
-        //add vertice and fragment data
-        const r2 = r / 2
-        let c,
-            nb = cells.length
-        const verticesBuffer = []
-        //const tBuffer = []
-        for (let i = 0; i < nb; i++) {
-            c = cells[i]
-            const t = this.tFun(c[this.colorCol], r, statColor)
-            if (t == null || t == undefined) continue
-            verticesBuffer.push(c.x + r2, c.y + r2)
-            //tBuffer.push(t > 1 ? 1 : t < 0 ? 0 : t)
-        }
-
-        const sizeGeo = this.size ? this.size(r, zf) : r + 0.2 * zf
-        const wgp = new WebGLSquareColoringKS(cvWGL.gl, this.colors, this.stretching, sizeGeo / zf)
-
-        //TODO - [i,j,t]
-        let data = []
-        for (const c of cells) {
-            const xGL = cg.geoToPixX(c.x + r2) //TODO within [-1,1] ?
-            data.push(xGL)
-            const yGL = cg.geoToPixY(c.y + r2) //TODO within [-1,1] ?
-            data.push(yGL)
-            const t = this.tFun(c[this.colorCol], r, statColor)
-            data.push(t)
-        }
-
-        data = [0.0, 2.0, 1.0, 0.0, 0.0, 0.0]
-
-        //draw
-        wgp.draw(verticesBuffer, data, cg.getWebGLTransform())
-
-        //draw in canvas geo
-        cg.initCanvasTransform()
-        cg.ctx.drawImage(cvWGL.canvas, 0, 0)
-
-        //update legends
-        this.updateLegends({ style: this, r: r, zf: zf, sColor: statColor })
-    }
-}