基于缩放级别和中心点创建latlnbounds

我有一个特征/多边形列表，我需要基于相对于特征的地图边界框运行操作，而不移动地图。

为了进一步解释我的问题，我可以举一个例子，说明如果移动地图不是问题，我如何做到这一点:

features.forEach(function (feature) {
    var bbox,
        ne,
        sw,
        fBounds,
        zoom,
        mBounds;
    bbox = feature.geometry.bbox;
    sw = L.latLng(bbox[1], bbox[0]);
    ne = L.latLng(bbox[3], bbox[2]);
    fBounds = L.latLngBounds(sw, ne);
    map.fitBounds(bounds);
    mBounds = map.getBounds();
    zoom = map.getZoom();
    //Execute operation based on mBounds and zoom
}    

我已经测试了很多，这是我所拥有的最接近工作代码片段的东西:

        var self = this,
            bbox,
            sw,
            ne,
            bounds,
            zoom,
            swPoint,
            nePoint,
            center,
            factor,
            dw,
            dh,
            cpx;
        bbox = feature.geometry.bbox;
        sw = L.latLng(bbox[1], bbox[0]);
        ne = L.latLng(bbox[3], bbox[2]);
        bounds = L.latLngBounds(sw, ne);
        zoom = self.map.getBoundsZoom(bounds, false); //maxZoom?
        swPoint = self.map.project(bounds.getSouthWest(), zoom),
        nePoint = self.map.project(bounds.getNorthEast(), zoom),
        center = self.map.unproject(swPoint.add(nePoint).divideBy(2), zoom);
        factor = self.map.options.crs.scale(zoom) / 8388608;
        dw = self.map.getSize().x / 2*factor;
        dh = self.map.getSize().y / 2*factor;
        cpx = self.map.options.crs.latLngToPoint(center, zoom);            
        return {
            ne: self.map.options.crs.pointToLatLng(L.point(cpx.x + dw, cpx.y - dh, false), zoom),
            sw: self.map.options.crs.pointToLatLng(L.point(cpx.x - dw, cpx.y + dh, false), zoom),
            center: center,
            zoom: zoom
        }
        //Execute operation based on returned object, repeat for every feature

这个'工作'，但它没有给出相同的结果作为第一个代码片段(即结果是错误的)。