gpt4 book ai didi

geojson - Shapefile 到 Topojson 的转换

转载 作者:行者123 更新时间:2023-12-02 02:57:03 27 4
gpt4 key购买 nike

我正在尝试转换可在 here 找到的加纳 admin1 shapefile 。我的最终目标是获得 TopoJSON,如本文 question 中所述。 。我在Linux机器上使用过这个命令:

ogr2ogr -f GeoJSON GHA_adm1.json GHA_adm1.shp

它返回这个:

Unable to open datasource `GHA_adm1.shp' with the following drivers.
... here goes a long list of drivers...

我是不是做错了什么?我应该安装其他“驱动程序”吗?但如何呢?谢谢。

最佳答案

下载 GIS 源

GADM 是管理 GIS 文件的完美来源。

GADM.org > Download page : select your country and format "shapefile" > ok.

或通过终端(将新加纳代码 GHA 替换为目标国家/地区的 ISO 代码):

# prepare folder
mkdir -p ./map; cd ./map
curl http://d3js.org/d3.v3.min.js -O -J
curl http://d3js.org/topojson.v1.min.js -O -J
# download data
curl \
-L -C - 'http://biogeo.ucdavis.edu/data/gadm2/shp/GHA_adm.zip' \
-o ./GHA_adm.zip
unzip -n ./GHA_adm.zip -d ./

shp 到 topojson

使用topojson command line ,这样更直接。如果您想保留所有属性:

 topojson -q 1e4 \
-o out.json \
-- in1.shp in2.shp

您还可以从 shapfile 中选择属性,并随时重命名它们:

 topojson \
--bbox \
--id-property none \
-p name=NAME_1 \
-p code=ID_1 \
-p L0=NAME_0 \
-q 1e4 \
--filter=small \
-o GHA_adm_w.topo.json \
-- admin_1=GHA_adm1.shp admin_2=GHA_adm2.shp

对于 GHA,.shp 中没有适合良好 ID 的属性。 NAME_1 包含空格,这会在您的 HTML 中提供无效的 id

enter image description here

检查 json

使用http://jsoneditoronline.org 。检查 json 将为您提供什么可用数据以及何处(点符号路径)的线索。 topojson distillery帮助预览任何 topojson 以及代码是否正确。

D3js 调用

<!DOCTYPE html>
<meta http-equiv="content-type" content="text/html; charset=utf-8">
<style>
.L1 {
fill: #E0E0E0;
stroke: #FFF;
stroke-width:1px;
}
</style>
<body>
<script src="./d3.v3.min.js"></script>
<script src="./topojson.v1.min.js"></script>
<script>
var mapIt = function(width, url){
console.log("mapIt(): start");
var height = width/960*500;
var svg = d3.select('body').append('svg')
.attr('width', width)
.attr('height', height);
var projection = d3.geo.mercator()
.scale(1)
.translate([0, 0]);

var path = d3.geo.path()
.projection(projection);

d3.json(url, function (error, json) {
var admin_1 = topojson.feature(json, json.objects.admin_1);

/* Autofocus code comes here ! */

svg.selectAll("path")
.data(admin_1.features)
.enter().append("path")
.attr('d', path)
.attr('class', 'L1');
});
};
mapIt(960,"http://somesite.org/data/NZL_adm.topo.json");
</script>
</body>
</html>

焦点

正确的自动对焦需要 Mike Bostocks 中的一小段代码,例如here :

// Compute the bounds of a feature of interest, then derive scale & translate.
var b = path.bounds(admin_1),
s = .95 / Math.max((b[1][0] - b[0][0]) / width, (b[1][1] - b[0][1]) / height),
t = [(width - s * (b[1][0] + b[0][0])) / 2, (height - s * (b[1][1] + b[0][1])) / 2];

// Update the projection to use computed scale & translate.
projection
.scale(s)
.translate(t);
<小时/>

编辑:现在应该可以工作了。现场演示:bl.ocks.org

关于geojson - Shapefile 到 Topojson 的转换,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/28556524/

27 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com