D3 制图：lon/lat 圆圈在地图上的错误位置（投影）

Question

我有一个关于 D3 制图的问题。我正在做一个小项目，我是 D3 的新手。

我从这个例子开始：http://bl.ocks.org/mbostock/5914438 我不想显示状态网格，而是想在地图上的某些位置 (lon/lat) 显示圆圈。我目前面临的问题是圆圈不在地图上的正确位置。我怀疑问题出在迈克使用的特殊投影上。他使用 1x1 正方形投影。这可能是显示图块所必需的。当我投影坐标时，值都在 -1 和 1 之间。我想我可以通过将它的宽度乘以高度和宽度来修复它，但它没有用。下面是我的代码（片段没有运行因为它缺少一个文件）。感谢您的协助！

<!DOCTYPE html>
<meta charset="utf-8">
<style>

body {
  margin: 0;
}

path {
  fill: none;
  stroke: red;
  stroke-linejoin: round;
  stroke-width: 1.5px;
}

circle {
  fill: #fff;
  fill-opacity: 0.4;
  stroke: #111;
}

</style>
<svg>
</svg>
<script src="//d3js.org/d3.v4.min.js"></script>
<script src="//d3js.org/d3-tile.v0.0.min.js"></script>
<script src="//d3js.org/topojson.v1.min.js"></script>
<script>
var pi = Math.PI,
    tau = 2 * pi;

var width = Math.max(960, window.innerWidth),
    height = Math.max(500, window.innerHeight);

// Initialize the projection to fit the world in a 1×1 square centered at the origin.
var projection = d3.geoMercator()
    .scale(1 / tau)
    .translate([0, 0]);

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

var tile = d3.tile()
    .size([width, height]);

var zoom = d3.zoom()
    .scaleExtent([1 << 9, 1 << 23])
    .on("zoom", zoomed);

var svg = d3.select("svg")
    .attr("width", width)
    .attr("height", height);

var raster = svg.append("g");

var vector = svg.append("path");

var circle = svg.append("g")
      
d3.json("/data/flyingsites/AD.json", function(error, flyingsites) {
  if (error) console.log(error); 
  
  // Compute the projected initial center.
  var center = projection([6.2, 45.8]);//45,809718, 6,252314

  // Apply a zoom transform equivalent to projection.{scale,translate,center}.
  svg
      .call(zoom)
      .call(zoom.transform, d3.zoomIdentity
          .translate(width / 2, height / 2)
          .scale(1 << 12)
          .translate(-center[0], -center[1]));

//add flying sites
  circle.selectAll("circle")
   .data(flyingsites.features)
  .enter().append("circle")
   .attr('r',5)
    .attr('cx',function(d) { return projection(d.geometry.coordinates)[0]*width})
    .attr('cy',function(d) { return projection(d.geometry.coordinates)[1]*height})
    .style('fill','red')
    
//console.log(flyingsites.features);
//console.log(circle);
});


function zoomed() {
  var transform = d3.event.transform;

  var tiles = tile
      .scale(transform.k)
      .translate([transform.x, transform.y])
      ();

  vector
      .attr("transform", transform)
      .style("stroke-width", 1 / transform.k);
     
  circle
      .attr("transform", "translate(" + transform.x + "," + transform.y + ")");
      
  var image = raster
      .attr("transform", stringify(tiles.scale, tiles.translate))
    .selectAll("image")
    .data(tiles, function(d) { return d; });

  image.exit().remove();

  image.enter().append("image")
      .attr("xlink:href", function(d) { return "http://" + "abc"[d[1] % 3] + ".tile.openstreetmap.org/" + d[2] + "/" + d[0] + "/" + d[1] + ".png"; })
      .attr("x", function(d) { return d[0] * 256; })
      .attr("y", function(d) { return d[1] * 256; })
      .attr("width", 256)
      .attr("height", 256);
}

function stringify(scale, translate) {
  var k = scale / 256, r = scale % 1 ? Number : Math.round;
  return "translate(" + r(translate[0] * scale) + "," + r(translate[1] * scale) + ") scale(" + k + ")";
}

</script>

Answer 1

您采用的方法行不通，一方面，它没有考虑规模（只是翻译）。这一点很关键，因为 d3-tile 使用几何缩放 - 它对所有矢量元素（不是图块）应用缩放变换（缩放和平移），这就是为什么投影将所有内容投影到一个像素的正方形区域并且永远不会随着缩放。

要解决此问题，请将您的圆圈放置在示例中放置（和大小）多边形的位置：

vector
  .attr("d", path(topojson.mesh(us, us.objects.counties)));

circle
  .attr("cx", projection(coord)[0])
  .attr("cy", projection(coord)[1])
  .attr("r", 5/(1<<12));

这两个位置的特点是一样的：只有投影，投影到一个像素的正方形。缩放应用变换以覆盖整个 svg。 此外，由于我们正在缩放那个像素以适合 svg，我们希望半径也能适当缩放。

现在我们可以像对多边形一样对圆应用变换：

circle
  .attr("transform", transform);

当然，我们也可以缩小每次缩放的半径，使用缩放 k 来修改圆的大小：

var pi = Math.PI,
    tau = 2 * pi;

var width = Math.max(960, window.innerWidth),
    height = Math.max(500, window.innerHeight);

// Initialize the projection to fit the world in a 1×1 square centered at the origin.
var projection = d3.geoMercator()
    .scale(1 / tau)
    .translate([0, 0]);

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

var tile = d3.tile()
    .size([width, height]);

var zoom = d3.zoom()
    .scaleExtent([1 << 11, 1 << 14])
    .on("zoom", zoomed);

var svg = d3.select("svg")
    .attr("width", width)
    .attr("height", height);

var raster = svg.append("g");

var vector = svg.append("circle");

  // Compute the projected initial center.
  var center = projection([-98.5, 39.5]);

  // Apply a zoom transform equivalent to projection.{scale,translate,center}.
  svg
      .call(zoom)
      .call(zoom.transform, d3.zoomIdentity
          .translate(width / 2, height / 2)
          .scale(1 << 12)
          .translate(-center[0], -center[1]));

  var coord = [-100,40]
    
  vector
      .attr("cx", projection(coord)[0])
   .attr("cy", projection(coord)[1])
   .attr("r", 5/(1<<12));

function zoomed() {
  var transform = d3.event.transform;

  var tiles = tile
      .scale(transform.k)
      .translate([transform.x, transform.y])
      ();

  vector
      .attr("transform", transform)
      .attr("r", 5/transform.k);

  var image = raster
      .attr("transform", stringify(tiles.scale, tiles.translate))
    .selectAll("image")
    .data(tiles, function(d) { return d; });

  image.exit().remove();

  image.enter().append("image")
      .attr("xlink:href", function(d) { return "http://" + "abc"[d[1] % 3] + ".tile.openstreetmap.org/" + d[2] + "/" + d[0] + "/" + d[1] + ".png"; })
      .attr("x", function(d) { return d[0] * 256; })
      .attr("y", function(d) { return d[1] * 256; })
      .attr("width", 256)
      .attr("height", 256);
}

function stringify(scale, translate) {
  var k = scale / 256, r = scale % 1 ? Number : Math.round;
  return "translate(" + r(translate[0] * scale) + "," + r(translate[1] * scale) + ") scale(" + k + ")";
}

<svg></svg>
<script src="https://d3js.org/d3.v4.min.js"></script>
<script src="https://d3js.org/d3-tile.v0.0.min.js"></script>
<script src="https://d3js.org/topojson.v1.min.js"></script>

最终，d3-tile 开始时可能会有点混乱，因为您使用了相当多的坐标系（平铺、缩放、像素、投影、地理），并且通常不会将整个地图投影到 1像素方块。

D3 制图：lon/lat 圆圈在地图上的错误位置（投影）

D3 cartography: lon/lat circles in wrong place on map (projection)

cartography

map-projections

d3.js