Openlayers利用kriging.js实现纯前端插值

最新推荐文章于 2024-03-15 09:36:10 发布
最新推荐文章于 2024-03-15 09:36:10 发布
阅读量1.5w 收藏 42
点赞数 5
分类专栏： OpenLayer
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_40184249/article/details/90521127
版权

关于空间插值这一块，也是GIS一个重要的分析。前端可以通过turf.js后端可以采用contour，如果仅有少量数据需要进行插值分析并进行展示，那应当如何做？可以采用turf.js但是turf.js插值效果是真的烂。通过检索发现可以通过kriging.js去做。kriging.js是基于canvas实现的在OL中可以通过imagecanvas来实现。

一、关于kriging.js
1、GitHub地址（https://github.com/oeo4b/kriging.js）
2、核心功能

  • kriging.train(t, x, y, model, sigma2, alpha):使用gaussian、exponential或spherical模型对数据集进行训练，返回的是一个variogram对象；
  • kriging.grid(polygons,variogram,width);使用刚才的variogram对象使polygons描述的地理位置内的格网元素具备不一样的预测值；
  • kriging.plot(canvas,grid,xlim,ylim,colors);将得到的格网grid渲染至canvas上。
    二、效果图

在这里插入图片描述
三、demo示例

<!DOCTYPE html>
<html>
<head>
    <meta charset="utf-8" />
    <title>前端空间插值</title>
    <style>
        #map {
            height: 100%;
            width: 100%;
        }
    </style>
    <link href="../script/ol4/ol.css" rel="stylesheet" />
    <script src="../script/ol4/ol.js"></script>
    <script src="../script/kriging.js-master/data.js"></script>
    <script src="../script/kriging.js-master/kriging.js"></script>
</head>
<body>
    <div id="map"></div>
    <script type="text/javascript">
        let params = {
            mapCenter: [116.40, 39.90],
            krigingModel: 'exponential',//model还可选'gaussian','spherical'
            krigingSigma2: 0,
            krigingAlpha: 226,
            canvasAlpha: 0.75,//canvas图层透明度
            colors: ["#00A600", "#01A600", "#03A700", "#04A700", "#05A800", "#07A800", "#08A900", "#09A900", "#0BAA00", "#0CAA00", "#0DAB00", "#0FAB00", "#10AC00", "#12AC00", "#13AD00", "#14AD00", "#16AE00", "#17AE00", "#19AF00", "#1AAF00", "#1CB000", "#1DB000", "#1FB100", "#20B100", "#22B200", "#23B200", "#25B300", "#26B300", "#28B400", "#29B400", "#2BB500", "#2CB500", "#2EB600", "#2FB600", "#31B700", "#33B700", "#34B800", "#36B800", "#37B900", "#39B900", "#3BBA00", "#3CBA00", "#3EBB00", "#3FBB00", "#41BC00", "#43BC00", "#44BD00", "#46BD00", "#48BE00", "#49BE00", "#4BBF00", "#4DBF00", "#4FC000", "#50C000", "#52C100", "#54C100", "#55C200", "#57C200", "#59C300", "#5BC300", "#5DC400", "#5EC400", "#60C500", "#62C500", "#64C600", "#66C600", "#67C700", "#69C700", "#6BC800", "#6DC800", "#6FC900", "#71C900", "#72CA00", "#74CA00", "#76CB00", "#78CB00", "#7ACC00", "#7CCC00", "#7ECD00", "#80CD00", "#82CE00", "#84CE00", "#86CF00", "#88CF00", "#8AD000", "#8BD000", "#8DD100", "#8FD100", "#91D200", "#93D200", "#95D300", "#97D300", "#9AD400", "#9CD400", "#9ED500", "#A0D500", "#A2D600", "#A4D600", "#A6D700", "#A8D700", "#AAD800", "#ACD800", "#AED900", "#B0D900", "#B2DA00", "#B5DA00", "#B7DB00", "#B9DB00", "#BBDC00", "#BDDC00", "#BFDD00", "#C2DD00", "#C4DE00", "#C6DE00", "#C8DF00", "#CADF00", "#CDE000", "#CFE000", "#D1E100", "#D3E100", "#D6E200", "#D8E200", "#DAE300", "#DCE300", "#DFE400", "#E1E400", "#E3E500", "#E6E600", "#E6E402", "#E6E204", "#E6E105", "#E6DF07", "#E6DD09", "#E6DC0B", "#E6DA0D", "#E6D90E", "#E6D710", "#E6D612", "#E7D414", "#E7D316", "#E7D217", "#E7D019", "#E7CF1B", "#E7CE1D", "#E7CD1F", "#E7CB21", "#E7CA22", "#E7C924", "#E8C826", "#E8C728", "#E8C62A", "#E8C52B", "#E8C42D", "#E8C32F", "#E8C231", "#E8C133", "#E8C035", "#E8BF36", "#E9BE38", "#E9BD3A", "#E9BC3C", "#E9BB3E", "#E9BB40", "#E9BA42", "#E9B943", "#E9B945", "#E9B847", "#E9B749", "#EAB74B", "#EAB64D", "#EAB64F", "#EAB550", "#EAB552", "#EAB454", "#EAB456", "#EAB358", "#EAB35A", "#EAB35C", "#EBB25D", "#EBB25F", "#EBB261", "#EBB263", "#EBB165", "#EBB167", "#EBB169", "#EBB16B", "#EBB16C", "#EBB16E", "#ECB170", "#ECB172", "#ECB174", "#ECB176", "#ECB178", "#ECB17A", "#ECB17C", "#ECB17E", "#ECB27F", "#ECB281", "#EDB283", "#EDB285", "#EDB387", "#EDB389", "#EDB38B", "#EDB48D", "#EDB48F", "#EDB591", "#EDB593", "#EDB694", "#EEB696", "#EEB798", "#EEB89A", "#EEB89C", "#EEB99E", "#EEBAA0", "#EEBAA2", "#EEBBA4", "#EEBCA6", "#EEBDA8", "#EFBEAA", "#EFBEAC", "#EFBFAD", "#EFC0AF", "#EFC1B1", "#EFC2B3", "#EFC3B5", "#EFC4B7", "#EFC5B9", "#EFC7BB", "#F0C8BD", "#F0C9BF", "#F0CAC1", "#F0CBC3", "#F0CDC5", "#F0CEC7", "#F0CFC9", "#F0D1CB", "#F0D2CD", "#F0D3CF", "#F1D5D1", "#F1D6D3", "#F1D8D5", "#F1D9D7", "#F1DBD8", "#F1DDDA", "#F1DEDC", "#F1E0DE", "#F1E2E0", "#F1E3E2", "#F2E5E4", "#F2E7E6", "#F2E9E8", "#F2EBEA", "#F2ECEC", "#F2EEEE", "#F2F0F0", "#F2F2F2"]
        };
        let baseLayer = new ol.layer.Tile({
            title: "base",
            source: new ol.source.OSM()
        });
        let map = new ol.Map({
            target: 'map',
            layers: [baseLayer],
            view: new ol.View({
                center: params.mapCenter,
                projection: 'EPSG:4326',
                zoom: 8
            })
        });
        let WFSVectorSource = new ol.source.Vector();
        let WFSVectorLayer = new ol.layer.Vector(
            {
                source: WFSVectorSource,
            });
        map.addLayer(WFSVectorLayer);
        //创建要素
        for (let i = 0; i < data.features.length; i++) {
            let feature = new ol.Feature({
                geometry: new ol.geom.Point([data.features[i].attributes.x, data.features[i].attributes.y]),
                value: data.features[i].attributes.z
            });
            feature.setStyle(new ol.style.Style({
                image: new ol.style.Circle({
                    radius: 6,
                    fill: new ol.style.Fill({ color: "#00F" })
                })
            }));
            WFSVectorSource.addFeature(feature);
        }
        //定义裁剪边界
        let coord = [[
            [117.315375, 40.181212],
            [117.367916, 40.135762],
            [116.751758, 40.002595],
            [116.754136, 39.870341],
            [116.913383, 39.804999],
            [116.901858, 39.680614],
            [116.788145, 39.56255],
            [116.535646, 39.590346],
            [116.392103, 39.491124],
            [116.4293, 39.433841],
            [116.387072, 39.417336],
            [116.237232, 39.476253],
            [116.172242, 39.578854],
            [115.728745, 39.479123],
            [115.610225, 39.588658],
            [115.508537, 39.59137],
            [115.416399, 39.733407],
            [115.416624, 39.776734],
            [115.550565, 39.772964],
            [115.408433, 40.015829],
            [115.85422, 40.144999],
            [115.922315, 40.216777],
            [115.708758, 40.499032],
            [115.89819, 40.585919],
            [116.03778, 40.577909],
            [116.208725, 40.741562],
            [116.454984, 40.739689],
            [116.297615, 40.910402],
            [116.43816, 40.870116],
            [116.405424, 40.94854],
            [116.537137, 40.9661],
            [116.621495, 41.057333],
            [116.703349, 40.853574],
            [116.93405, 40.675155],
            [117.454502, 40.647216],
            [117.387854, 40.533274],
            [117.166811, 40.503979],
            [117.164198, 40.373887],
            [117.315375, 40.181212]]];       
        let clipgeom = new ol.geom.Polygon(coord);            
        //绘制kriging插值图
        let canvasLayer = null;
        let drawKriging =function (extent){
            let values = [], lngs = [], lats = [];
            WFSVectorSource.forEachFeature(function (feature) {
                values.push(feature.getProperties().value);
                lngs.push(feature.getGeometry().getCoordinates()[0]);
                lats.push(feature.getGeometry().getCoordinates()[1]);
            })
            console.log(values.length)
            if (values.length > 3) {
                let letiogram = kriging.train(values, lngs, lats,
                    params.krigingModel, params.krigingSigma2, params.krigingAlpha);
                let ex = clipgeom.getExtent();
                let grid = kriging.grid(coord, letiogram, (ex[2] - ex[0]) / 200);               
                //移除已有图层
                if (canvasLayer !== null) {
                    map.removeLayer(canvasLayer);
                }
                //创建新图层
                canvasLayer = new ol.layer.Image({
                    source: new ol.source.ImageCanvas({
                        canvasFunction: (extent, resolution, pixelRatio, size, projection) => {
                            console.log(extent);
                            let canvas = document.createElement('canvas');
                            canvas.width = size[0];
                            canvas.height = size[1];
                            canvas.style.display = 'block';
                            //设置canvas透明度
                            canvas.getContext('2d').globalAlpha = params.canvasAlpha;
                            //使用分层设色渲染
                            kriging.plot(canvas, grid,
                                [extent[0], extent[2]], [extent[1], extent[3]], params.colors);                            
                            return canvas;
                        },
                        projection: 'EPSG:4326'
                    })
                })             
                //向map添加图层
                map.addLayer(canvasLayer);
            } else {
                alert("有效样点个数不足，无法插值");
            }
        }      
        //首次加载，自动渲染一次差值图
        let extent = clipgeom.getExtent();
        drawKriging(extent);
        
    </script>
</body>
</html>

四、参考文章
1、https://www.jianshu.com/p/c8473ac0b08a
2、https://www.cnblogs.com/naaoveGIS/p/9444480.html
3、https://www.cnblogs.com/giserhome/p/10851799.html
4、https://mathiasleroy.com/georeferenced-dynamic-kriging/

五、data.js数据


  var data={
        "features": [
            {
                "attributes": {
                    "FID": 0,
                    "id": 1,
                    "name": "Beijing US Embassy",
                    "x": 116.468,
                    "y": 39.954999999999998,
                    "z": 46,
                    "xu": 454558.72859999997,
                    "yu": 4422898.159
                },
                "geometry": {
                    "x": 116.468,
                    "y": 39.954999999999998
                }
            },
            {
                "attributes": {
                    "FID": 1,
                    "id": 2,
                    "name": "Temple of Heaven",
                    "x": 116.407,
                    "y": 39.886000000000003,
                    "z": 36,
                    "xu": 449297.45990000002,
                    "yu": 4415272.716
                },
                "geometry": {
                    "x": 116.407,
                    "y": 39.886000000000003
                }
            },
            {
                "attributes": {
                    "FID": 2,
                    "id": 3,
                    "name": "Haidian Beijing Botanical Garden",
                    "x": 116.20699999999999,
                    "y": 40.002000000000002,
                    "z": 40,
                    "xu": 432311.35320000001,
                    "yu": 4428280.3169999998
                },
                "geometry": {
                    "x": 116.20699999999999,
                    "y": 40.002000000000002
                }
            },
            {
                "attributes": {
                    "FID": 3,
                    "id": 4,
                    "name": "Chaoyang Olympic Sports Center116.397",
                    "x": 116.39700000000001,
                    "y": 39.981999999999999,
                    "z": 43,
                    "xu": 448514.42090000003,
                    "yu": 4425933.4840000002
                },
                "geometry": {
                    "x": 116.39700000000001,
                    "y": 39.981999999999999
                }
            },
            {
                "attributes": {
                    "FID": 4,
                    "id": 5,
                    "name": "Chaoyang Agricultural Exhibition Hall",
                    "x": 116.461,
                    "y": 39.936999999999998,
                    "z": 52,
                    "xu": 453948.74660000001,
                    "yu": 4420903.926
                },
                "geometry": {
                    "x": 116.461,
                    "y": 39.936999999999998
                }
            },
            {
                "attributes": {
                    "FID": 5,
                    "id": 6,
                    "name": "Liangxiang",
                    "x": 116.136,
                    "y": 39.741999999999997,
                    "z": 123,
                    "xu": 425971.88909999997,
                    "yu": 4399479.2439999999
                },
                "geometry": {
                    "x": 116.136,
                    "y": 39.741999999999997
                }
            },
            {
                "attributes": {
                    "FID": 6,
                    "id": 7,
                    "name": "Fengtai Yungang",
                    "x": 116.146,
                    "y": 39.823999999999998,
                    "z": 85,
                    "xu": 426915.51539999997,
                    "yu": 4408572.0729999999
                },
                "geometry": {
                    "x": 116.146,
                    "y": 39.823999999999998
                }
            },
            {
                "attributes": {
                    "FID": 7,
                    "id": 8,
                    "name": "Shunyi New Town",
                    "x": 116.655,
                    "y": 40.127000000000002,
                    "z": 154,
                    "xu": 470605.50309999997,
                    "yu": 4441910.1670000004
                },
                "geometry": {
                    "x": 116.655,
                    "y": 40.127000000000002
                }
            },
            {
                "attributes": {
                    "FID": 8,
                    "id": 12,
                    "name": "Donggaocun Zhen",
                    "x": 117.12,
                    "y": 40.100000000000001,
                    "z": 172,
                    "xu": 510228.20750000002,
                    "yu": 4438863.2359999996
                },
                "geometry": {
                    "x": 117.12,
                    "y": 40.100000000000001
                }
            },
            {
                "attributes": {
                    "FID": 9,
                    "id": 13,
                    "name": "Tongzhou New Town",
                    "x": 116.663,
                    "y": 39.886000000000003,
                    "z": 216,
                    "xu": 471185.97159999999,
                    "yu": 4415158.7980000004
                },
                "geometry": {
                    "x": 116.663,
                    "y": 39.886000000000003
                }
            },
            {
                "attributes": {
                    "FID": 10,
                    "id": 14,
                    "name": "Huangcunzhen",
                    "x": 116.404,
                    "y": 39.718000000000004,
                    "z": 226,
                    "xu": 448916.78909999999,
                    "yu": 4396628.5209999997
                },
                "geometry": {
                    "x": 116.404,
                    "y": 39.718000000000004
                }
            },
            {
                "attributes": {
                    "FID": 11,
                    "id": 15,
                    "name": "BDA",
                    "x": 116.506,
                    "y": 39.795000000000002,
                    "z": 216,
                    "xu": 457706.38530000002,
                    "yu": 4405121.3250000002
                },
                "geometry": {
                    "x": 116.506,
                    "y": 39.795000000000002
                }
            },
            {
                "attributes": {
                    "FID": 12,
                    "id": 16,
                    "name": "Yufazhen",
                    "x": 116.3,
                    "y": 39.520000000000003,
                    "z": 204,
                    "xu": 439831.65490000002,
                    "yu": 4374718.0180000002
                },
                "geometry": {
                    "x": 116.3,
                    "y": 39.520000000000003
                }
            },
            {
                "attributes": {
                    "FID": 13,
                    "id": 17,
                    "name": "Yongledianzhen",
                    "x": 116.783,
                    "y": 39.712000000000003,
                    "z": 198,
                    "xu": 481399.35399999999,
                    "yu": 4395815.3370000003
                },
                "geometry": {
                    "x": 116.783,
                    "y": 39.712000000000003
                }
            },
            {
                "attributes": {
                    "FID": 14,
                    "id": 18,
                    "name": "Xianghe EPA",
                    "x": 117.009,
                    "y": 39.765999999999998,
                    "z": 204,
                    "xu": 500770.85320000001,
                    "yu": 4401786.0719999997
                },
                "geometry": {
                    "x": 117.009,
                    "y": 39.765999999999998
                }
            },
            {
                "attributes": {
                    "FID": 15,
                    "id": 19,
                    "name": "Badaling Northwest",
                    "x": 115.988,
                    "y": 40.365000000000002,
                    "z": 46,
                    "xu": 414076.95390000002,
                    "yu": 4468761.409
                },
                "geometry": {
                    "x": 115.988,
                    "y": 40.365000000000002
                }
            },
            {
                "attributes": {
                    "FID": 16,
                    "id": 20,
                    "name": "East Fourth Ring Road",
                    "x": 116.483,
                    "y": 39.939,
                    "z": 181,
                    "xu": 455829.68479999999,
                    "yu": 4421114.7860000003
                },
                "geometry": {
                    "x": 116.483,
                    "y": 39.939
                }
            },
            {
                "attributes": {
                    "FID": 17,
                    "id": 21,
                    "name": "Yanqing town",
                    "x": 115.97199999999999,
                    "y": 40.453000000000003,
                    "z": 59,
                    "xu": 412832.08510000003,
                    "yu": 4478545.159
                },
                "geometry": {
                    "x": 115.97199999999999,
                    "y": 40.453000000000003
                }
            },
            {
                "attributes": {
                    "FID": 18,
                    "id": 22,
                    "name": "Miyun Reservoir",
                    "x": 116.911,
                    "y": 40.499000000000002,
                    "z": 63,
                    "xu": 492458.57429999998,
                    "yu": 4483147.5360000003
                },
                "geometry": {
                    "x": 116.911,
                    "y": 40.499000000000002
                }
            },
            {
                "attributes": {
                    "FID": 19,
                    "id": 23,
                    "name": "Haidian Wanliu",
                    "x": 116.28700000000001,
                    "y": 39.987000000000002,
                    "z": 180,
                    "xu": 439126.71720000001,
                    "yu": 4426557.75
                },
                "geometry": {
                    "x": 116.28700000000001,
                    "y": 39.987000000000002
                }
            },
            {
                "attributes": {
                    "FID": 20,
                    "id": 24,
                    "name": "Yongdingmen Inner St",
                    "x": 116.39400000000001,
                    "y": 39.875999999999998,
                    "z": 196,
                    "xu": 448178.40250000003,
                    "yu": 4414170.284
                },
                "geometry": {
                    "x": 116.39400000000001,
                    "y": 39.875999999999998
                }
            },
            {
                "attributes": {
                    "FID": 21,
                    "id": 25,
                    "name": "Jianshe Road",
                    "x": 117.304,
                    "y": 39.719000000000001,
                    "z": 55,
                    "xu": 526055.42509999999,
                    "yu": 4396613.8959999997
                },
                "geometry": {
                    "x": 117.304,
                    "y": 39.719000000000001
                }
            },
            {
                "attributes": {
                    "FID": 22,
                    "id": 26,
                    "name": "Ligong, Chengde",
                    "x": 117.938,
                    "y": 41.011000000000003,
                    "z": 70,
                    "xu": 578874.61600000004,
                    "yu": 4540401.8559999997
                },
                "geometry": {
                    "x": 117.938,
                    "y": 41.011000000000003
                }
            },
            {
                "attributes": {
                    "FID": 23,
                    "id": 27,
                    "name": "Fengning County City Hall",
                    "x": 116.652,
                    "y": 41.215000000000003,
                    "z": 59,
                    "xu": 470827.95390000002,
                    "yu": 4562682.9199999999
                },
                "geometry": {
                    "x": 116.652,
                    "y": 41.215000000000003
                }
            },
            {
                "attributes": {
                    "FID": 24,
                    "id": 28,
                    "name": "Xiahuayuan EPA",
                    "x": 115.29600000000001,
                    "y": 40.508000000000003,
                    "z": 61,
                    "xu": 355627.29269999999,
                    "yu": 4485537.4510000004
                },
                "geometry": {
                    "x": 115.29600000000001,
                    "y": 40.508000000000003
                }
            },
            {
                "attributes": {
                    "FID": 25,
                    "id": 29,
                    "name": "Yuxian Secondary School",
                    "x": 114.596,
                    "y": 39.845999999999997,
                    "z": 61,
                    "xu": 294324.1188,
                    "yu": 4413430.4780000001
                },
                "geometry": {
                    "x": 114.596,
                    "y": 39.845999999999997
                }
            },
            {
                "attributes": {
                    "FID": 26,
                    "id": 30,
                    "name": "Zhouzhou Monitoring Stations",
                    "x": 116.03400000000001,
                    "y": 39.491999999999997,
                    "z": 80,
                    "xu": 416933.99089999998,
                    "yu": 4371822.0190000003
                },
                "geometry": {
                    "x": 116.03400000000001,
                    "y": 39.491999999999997
                }
            },
            {
                "attributes": {
                    "FID": 27,
                    "id": 31,
                    "name": "Guan Party School, Langfang",
                    "x": 116.30500000000001,
                    "y": 39.445,
                    "z": 70,
                    "xu": 440197.23119999998,
                    "yu": 4366391.0480000004
                },
                "geometry": {
                    "x": 116.30500000000001,
                    "y": 39.445
                }
            },
            {
                "attributes": {
                    "FID": 28,
                    "id": 32,
                    "name": "Guyuan County welfare",
                    "x": 115.681,
                    "y": 41.667999999999999,
                    "z": 63,
                    "xu": 390196.44390000001,
                    "yu": 4613756.0429999996
                },
                "geometry": {
                    "x": 115.681,
                    "y": 41.667999999999999
                }
            }
        ]
}
hpugisers
关注
  • 5
    点赞
  • 踩
  • 42
    收藏
    觉得还不错? 一键收藏
  • 7
    评论
专栏目录
掌握 OpenLayers：集成 Canvas、ECharts 和 Turf.js 的终极指南
qq_44103359的博客
05-23 396
OpenLayers 是一个强大的开源地图库，用于在网页上展示地图数据。它支持多种地图源，包括 Google Maps、Bing Maps、OpenStreetMap 等，并且可以轻松与各种前端技术如 Canvas、ECharts 和 Turf.js 集成。本教程将指导您如何使用 OpenLayers，以及如何将其与 Canvas、ECharts 和 Turf.js 结合使用，以创建丰富、交互式的地图应用。
kriging克里金插值以及前端渲染jS代码部分解释
振云博客
05-14 2761
.html中的方法调用 //训练使用高斯过程与贝叶斯先验 let variogram=kriging.train(positions.map(pos=>pos[2]),positions.map(pos=>pos[0]),positions.map(pos=>pos[1]),params.krigingModel,params.krigingSigma2,params.krigingAlpha); //网格矩阵或轮廓路径 let grid=kriging.grid(po
Kriging 算法实现 2维和3维地图等高线
02-23
A year and a half year ago, I published this article to the Codeguru site and got a number of requests about the Kriging algorithm contour map. Unfortunately, my project was changed shortly after that article and later I quit the company so I couldn‘t find time to finish this Contour business. A week ago, I happened to need a contour map again so I decided to solve the Kriging algorithm. I searched the Internet for a commercial library but they all look ugly and hard to use. So, I made up my mind to make my own algorithm. The Kriging algorithm is easy to find, but this algorithm needs a Matrix and solver (LU-Decomposition). Again, I couldn‘t find suitable code for this. I tried to use GSL first but this made my code too big and was slower. Finally, I went back to "Numerical Recipe in C"—yes, that horrible-looking C code—and changed the code there to my taste.If you read this article before, the rendering part hasn‘t been changed much. I added the Kriging algorithm and revised the codes a little bit. Following is the Kriging Algorithm:templatedouble GetDistance(const ForwardIterator start, int i, int j){ return ::sqrt(::pow(((*(start+i)).x - (*(start+j)).x), 2) + ::pow(((*(start+i)).y - (*(start+j)).y), 2));}templatedouble GetDistance(double xpos, double ypos, const ForwardIterator start, int i){ return ::sqrt(::pow(((*(start+i)).x - xpos), 2) + ::pow(((*(start+i)).y - ypos), 2));}templateclass TKriging : public TInterpolater{public: TKriging(const ForwardIterator first, const ForwardIterator last, double dSemivariance) : m_dSemivariance(dSemivariance) { m_nSize = 0; ForwardIterator start = first; while(start != last) { ++m_nSize; ++start; } m_matA.SetDimension(m_nSize, m_nSize); for(int j=0; j<m_nSize; j++) { for(int i=0; i<m_nSize; i++) { if(i == m_nSize-1 || j == m_nSize-1) { m_matA(i, j) = 1; if(i == m_nSize-1 && j == m_nSize-1) m_matA(i, j) = 0; continue; } m
Kriging.js - 地统计插值库
最新发布
gitblog_00069的博客
03-15 702
Kriging.js - 地统计插值库 Kriging.js 是一个用于地统计插值的 JavaScript 库。它提供了多种插值方法，并且可以在 Web 浏览器中运行。 什么是 Kriging？ Kriging 是一种用于估计未知点的变量值的技术，基于已知点的数据。这种方法广泛应用于地球科学、环境学、地理信息系统等领域。Kriging 的主要思想是利用空间相关性来预测未知点的变量值。 Krigin...
Cesium+kriging.js实现雨量插值
weixin_44265800的博客
11-17 9884
0.前言 网上有很多关于openlayers的克里金插值，但是最近在学习cesium，也想在cesium中做，但是好像网上相关的例子一个都没有。所以我就自己尝试去做。 1.克里金插值 克里金插值也称作空间局部插值法，或空间自协方差最佳插值法，它以变异函数理论和结构分析为基础， 在特定区域内对区域化变量进行无偏最优估计，被广泛应用于土壤制图领域，是一种非常重要的地统计学方法。 主要是用来预测一些地理...
kriging.zip
08-24
克里金差值渲染演示代码，基于vue+threejs
OpenLayers6（5）：基于kriging.js实现插值渲染图
碰碰的博客
04-27 2043
1 版本 OpenLayers：6.4.3 2 相关配置 import kriging from '@sakitam-gis/kriging'; 3 使用kriging.js 3.1 数据准备 待准备数据以通用的geojson数据格式即可 待插值得矢量点数据：本例使用河南省XXX点数据 待插值的范围边界数据：本例使用河南省边界 3.2 样本训练 The train method in the kriging object fits your input to whatever va
openlayers : ol.js、ol.css 下载
07-23
`ol.js`是OpenLayers的核心JavaScript库，包含了所有的地图渲染、图层管理、交互操作等功能；而`ol.css`则是配套的样式表文件，提供了默认的用户界面样式和地图元素的布局。 1. **OpenLayers 概述** OpenLayers ...
openlayers : ol.css ol.js
04-09
`ol.js` 文件则是OpenLayers的核心JavaScript库，包含了所有地图操作、图层管理、事件处理等功能的实现。通过这个文件，开发者可以创建地图实例、加载不同类型的图层、添加交互功能、处理用户事件等。它提供了一套...
wegue:使用OpenLayers和Vue.js进行Webmapping应用程序的模板和组件
04-28
Wegue（带有OpenLayers和Vue.js的WebGIS） 使用OpenLayers和Vue.js进行Webmapping应用程序的模板和可重用组件 转到上的在线演示关于Wegue（的We b摹用的OpenLayers和V UE被）结合的功率和地理空间精明使轻量级电子...
克里金算法kriging.js
03-24
克里金算法
openlayers中文api.zip
07-29
OpenLayers 是一个强大的开源JavaScript库，专用于在Web上创建交互式地图应用。这个"openlayers中文api.zip"文件包含了OpenLayers中文版的API文档，为开发者提供了详细且易懂的指南，帮助他们更好地理解和使用这个库...
oeo4b.github.io:kriging.jsHTML页面
05-15
4. **kriging.js库**：此项目可能包含一个名为"kriging.js"的自定义JavaScript库，用于实现kriging算法。理解kriging的统计原理，包括变程分析、协方差函数、最佳线性无偏预测（OLUP）和误差变异函数等概念，对于...
根据经纬度在layui弹出层中显示百度地图
奔跑的菜鸡
12-22 1204
先上效果图 实现js代码 首选要引入layui，这个大家可以到layui官网看看怎么引入layer 附上layui官网url https://www.layui.com/doc/ 注意：这里要提示一点content中内容里的div一定要给一个宽高，不然地图是加载不出来的 return layer.open({ type: 1, title: false, closeBtn: false, area: '500px',
leaflet使用百度地图
Kevinvcc200的博客
05-16 6984
参考：https://www.jianshu.com/p/b48ec99011c2 leaflet使用百度地图，叠加层使用硬件/谷歌坐标。 leaflet默认使用硬件/谷歌坐标，使用百度地图需要进行坐标偏移。就是把硬件/谷歌坐标转换成百度坐标。 网上搜到一些解决方法，但是总是不奏效。最后还是拿网上的改了下 参考了以下文章：https://blog.csdn.net/u012087400/ar...
Leaflet自定义图片作为地图
VaeEcho的博客
05-04 1201
vue中使用leaflet展示室内地图，并解决图片初始化问题；
Leaflet使用GP服务实现在线插值
孙霸天的专栏
07-31 1574
前期准备： 使用ArcMap发布GP工具 使用leaflet官方的js和css 点数据准备   正式开始： 发布gp工具的流程不说了直接上GP截图 在发布GP服务的时候需要注意几点 为了保证GP的插值是在我们所选的范围内执行，需要设置它的插值范围，右击IDW选择属性-&gt;处理范围-&gt;选择你要的范围 如果需要给插值服务设置颜色，就必须给输出的结果设置图层符号系统，右击属...
bar图设置距离 python_Python可视化 | 反距离加权法插值及可视化绘制
weixin_39715926的博客
12-18 563
本期推文，我们将介绍IDW(反距离加权法(Inverse Distance Weighted))插值的Python计算方法及插值结果的可视化绘制过程。主要涉及的知识点如下：IDW简介自定义Python代码计算空间IDW分别使用plotnine、Basemap进行IDW插值结果可视化绘制IDW简介反距离权重 (IDW)插值假设：彼此距离较近的事物要比彼此距离较远的事物更相似。当为任何未测量的位置预...
openlayers 报错VectorLayer.js:179 Uncaught TypeError: Cannot read properties of undefined (reading 'opacity')
11-29
这个错误通常是由于openlayers中的矢量图层（VectorLayer）没有正确设置透明度（opacity）属性导致的。解决此问题的方法是在创建矢量图层时设置透明度属性。以下是一个示例代码： ```javascript var vectorLayer = new ol.layer.Vector({ opacity: 0.5, // 设置透明度为0.5 source: new ol.source.Vector({ url: 'data/geojson/countries.geojson', format: new ol.format.GeoJSON() }) }); ```

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 7
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值