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

2024-05-23

我正在尝试使用 geoproject 在 geojson 文件上设置投影。具体来说,我正在尝试将投影设置为 BCalbers (http://spatialreference.org/ref/epsg/3005/ 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  " \

D3 地理投影支持一系列通用投影(和自定义投影),通常可以重新创建任何给定的特定投影。但是,重新创建投影时,d3 投影不会复制地图单位*。 D3 投影在 svg 坐标空间中创建单位(投影坐标将从[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。这就是您所需要的 - (假)东移/北移参考是为了重新创建地图单位 - 在 Web 场景中通常不需要这(如果需要,使用更完整的 GIS 平台会更合适)。

因此,要设置投影,您可以使用:

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

Where

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 阿尔伯斯,我可能会选择 50.5,因为它位于育空地区边界和华盛顿边界之间的中间位置,这是 BC 的北部和南部边界(嗯,温哥华岛和一些海湾群岛除外,所以假设 50 度北,抱歉我忘了你维多利亚)。

另请记住,d3 的投影功能假定数据未投影(或在 WGS84 中“投影”),由长纬度对组成。

您可以在 Mike Bostock 的命令行制图文章中了解如何使用这些方法here https://medium.com/@mbostock/command-line-cartography-part-1-897aa8f8ca2c:

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

fitSize在这个例子中,缩放和平移特征到指定的边界框 - 这个平移和缩放不会改变投影参数,就像center方法,不会改变距离、面积、形状或角度(距离和面积仍然与适当的 BC Albers 成比例)。

*您可以重新创建地图单位(东东/北东错误可能需要一些自定义投影工作),但这实际上并不是它的平台,使用许多其他平台会更容易。

另请参阅此问题和答案:将 EPSG 投影范围转换为 D3.js https://stackoverflow.com/q/42259132/7106086.

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

使用 d3-geo-projection 命令行工具设置自定义投影 的相关文章

随机推荐