Shapefile到Topojson的转换

时间：2015-02-17 07:11:04

标签： geojson shapefile topojson

我正在尝试转换可以找到的加纳admin1 shapefile here。我的最终目标是获得这个question中描述的TopoJSON。 我在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...

我做错了什么？我应该安装其他＆＃34;驱动程序＆＃34;？但是怎么样？感谢。

2 个答案:

答案 0 :(得分：6)

下载GIS源

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

  

GADM.org＆gt; Download page：选择您的国家/地区＆＃34; shapefile＆＃34; ＆GT;确定。

或通过终端（用目标国家/地区的ISO代码替换新加纳代码GHA）：

# 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 to topojson

使用topojson command line，它更直接。如果您想保留所有属性：

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

您还可以从shapfiles中选择属性，并在go：

中重命名它们 
 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

答案 1 :(得分：1)

在这里工作得非常好：

  

GDAL 1.11.1，2014/09/24发布

也许你正在运行一个过时的GDAL版本？您可以使用以下网址检查您的版本：ogr2ogr --version

相关问题
最新问题