shp2json生成GeoJSON,其边界和坐标不在​​(-180,180)范围内

Question

我正在做@ mbostock 关于命令行制图的教程 1.以下是我所做的:

1)下载zip文件

wget https://www.arb.ca.gov/ei/gislib/gislib.htm/ca_air_basins.zip (来自加州空气资源委员会网站)

2)通过单击.zip文件手动解压缩(unzip命令不起作用,因为下载的文件没有'结束中心目录签名')

2)将shapefile转换为GeoJSON shp2json CaAirBasin.shp -o ca.json

3)使用California Albers(在CaAirBasin.prj文件中指定)进行GeoJSON投影:

geoproject 'd3.geoConicEqualArea().parallels([34, 40.5]).rotate([120, 0]).fitSize([960, 960], d)' < ca.json > ca-albers.json

结果:我得到ca-albers.json(223M)比ca.json(6.3M)大得多.在本教程中使用的人口普查区shapefile中,文件大小从10增加到14 M.

更新:似乎问题是shp2json步骤而不是geoProject步骤.看着CaAirBasin.shp.xml,有两组边界的,bounding并且lbounding.在我的情况,lbounding和坐标是在100,000s范围内使用,而不是在教程(-180,180)的值,使得投影失败.

有什么好的下一步 - 是否有替代方案shp2json可以适用于我的案例？或者如何将我的边界/坐标转换为适当的(-180,180)？

谢谢.

Answer 1

你的shapefile已经被预测了.通过d3.geoConicEqualArea传递它,你将它投射两次; d3.geoConicEqualArea期望WGS84(以度为单位的经度和纬度)作为输入,但是你给它投影的坐标.

您可以通过查看PRJ文件或通过查看shp2json生成的几何体的坐标来判断您的shapefile是否已投射:第一个点是[-298228.39936644124,437516.87775637675],它远远超出预期的范围[±180] °,±90°]经度和纬度.

不是将几何体重新投影到新投影,最简单的方法是使用它已经存在的投影.您需要做的就是平移和缩放以适应所需的960×960显示尺寸.您可以使用d3.geoIdentity.

geoproject 'd3.geoIdentity().reflectY(true).fitSize([960, 960], d)' < ca.json > ca-albers.json

设置标识 .reflectY也将修复几何体的垂直方向:GIS约定是+ y指向,但Canvas和SVG约定是+ y指向下方.