使用 d3-geo-projection 命令行工具设置自定义投影

Question

我正在尝试使用 geoproject 在 geojson 文件上设置投影。具体来说，我正在尝试将投影设置为 BCalbers (http://spatialreference.org/ref/epsg/3005/)

我看到 geoproject 有许多投影选项，即

geoproject 'd3.geoAlbersUsa()' us.json \
  > us-albers.json

但是可以使用命令行工具设置自定义投影吗？我希望这样的事情是可能的：

geoproject '+proj=aea +lat_1=50 +lat_2=58.5 +lat_0=45 +lon_0=-126 +x_0=1000000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs' build/airzones.geojson \
  > bc-albers.json

但没有骰子。这适用于 ogr2ogr

ogr2ogr -f GeoJSON -t_srs "+proj=aea +lat_1=50 +lat_2=58.5 +lat_0=45 +lon_0=-126 +x_0=1000000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs  " \

Answer 1

D3 地理投影支持一系列通用投影（和自定义投影），通常可以重新创建任何给定的特定投影。但是，在重新创建投影时，d3 投影不会复制地图单位*。 D3 投影在 svg 坐标 space 中创建单位（投影坐标将从左上角 [0,0] 开始）。这允许跳过投影数据的步骤（通常使用 d3 在运行中），然后重新缩放和平移它以显示地图。

比较：使用引用的 d3.geoAlbersUsa()，如果使用默认比例，d3 将沿 x 轴投影大约 960 个像素的数据。非常适合在 SVG 中未经修改的使用——坐标可以用作直线 SVG 坐标。在 ArcGIS 或 QGIS 中使用 BC Albers 将投影数百万米的数据。

但是，除了单位之外，您可以重新创建按比例缩放的 BC Albers，以保留常规 BC Albers 的形状、距离、方向和面积。但由于 d3 不采用 .prj 文件或其他定义等投影定义，因此您需要使用 d3 的投影方法和适当的参数。

对于 BC Albers，您的纬线是：50 和 58.5，您的中心经度是 -126，您的投影类型是 Albers。这就是您所需要的 - (false) easting/northing 参考是重新创建地图单元 - 在 web 场景中通常不需要它（如果需要，使用更完整的 GIS 平台会更合适）。

因此，要设置投影，您可以使用：

d3.geoAlbers()
  .center([0,y])
  .rotate([-x,0])
  .parallels([a,b])
  .scale(k)

哪里

x = center longitude (negative because we rotate the globe under the map)

a, b = standard parallels

k = scale factor(for a d3.geoAlbers(): whole world is 960 px across with a scale of 1070, the default scale, larger numbers expand this/zoom in)

y = centering latitude.

注意：y 不会改变投影，它只是对其进行平移 - BC Albers 的 y 参考在 BC 以南，它只是北向的参考，因为它对地图形状、面积、距离，或方向。对于 BC Albers，我可能会选择 50.5，因为它位于育空边界和华盛顿边界之间，这是 BC 的北部和南部边界（嗯，除了温哥华岛和一些海湾群岛，所以假设 50 度北方，对不起，我忘了你维多利亚）。

另请记住，d3 的投影功能假设数据是未投影的（或 WGS84 中的 "projected"），由长纬度对组成。

您可以在他的命令行制图文章中看到 Mike Bostock 如何使用这些方法 here:

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

fitSize 在这个例子中将特征缩放和平移到指定的边界框 - 这种平移和缩放不会改变投影参数，就像 [=16 中的 y 坐标=] 方法，不会改变距离、面积、形状或角度（嗯，距离和面积与适当的 BC Albers 保持成比例）。

*您可以重新创建地图单位（false eastings/northings 可能需要一些自定义投影工作），但这并不是真正适合它的平台，使用许多其他平台会更容易。

另请参阅此问答：。

使用 d3-geo-projection 命令行工具设置自定义投影

Setting a custom projection using d3-geo-projection command line tool

map-projections

d3.js