first working map

dwtkns · Nov 1, 2012 · 4fa5889 · 4fa5889
1 parent 7ff55c3
commit 4fa5889
Show file tree

Hide file tree

Showing 8 changed files with 1,088 additions and 21 deletions.
diff --git a/.DS_Store b/.DS_Store
diff --git a/README.md b/README.md
@@ -1,4 +1,8 @@
 d3_dymaxion
 ===========
 
-d3 dymaxion fiddling
+Fiddling with implementing Protovis' dymaxion projection in D3.
+
+Code hacked together from:
+http://mbostock.github.com/protovis/ex/dymax.html
+http://bl.ocks.org/3712397
diff --git a/boundaries.json b/boundaries.json
diff --git a/d3_dymaxion.js b/d3_dymaxion.js
@@ -0,0 +1,107 @@
+(function() {
+  var ε = 1e-6,
+      π = Math.PI,
+      sqrtπ = Math.sqrt(π);
+
+  var robinsonConstants = [
+    [1.0000, 0.0000],
+    [0.9986, 0.0620],
+    [0.9954, 0.1240],
+    [0.9900, 0.1860],
+    [0.9822, 0.2480],
+    [0.9730, 0.3100],
+    [0.9600, 0.3720],
+    [0.9427, 0.4340],
+    [0.9216, 0.4958],
+    [0.8962, 0.5571],
+    [0.8679, 0.6176],
+    [0.8350, 0.6769],
+    [0.7986, 0.7346],
+    [0.7597, 0.7903],
+    [0.7186, 0.8435],
+    [0.6732, 0.8936],
+    [0.6213, 0.9394],
+    [0.5722, 0.9761],
+    [0.5322, 1.0000]
+  ];
+
+  function sinci(x) {
+    return x ? x / Math.sin(x) : 1;
+  }
+
+  function sgn(x) {
+    return x > 0 ? 1 : x < 0 ? -1 : 0;
+  }
+
+  function asqrt(x) {
+    return x > 0 ? Math.sqrt(x) : 0;
+  }
+
+  // Calculate F(φ+iψ|m).
+  // See Abramowitz and Stegun, 17.4.11.
+  function ellipticFi(φ, ψ, m) {
+    var r = Math.abs(φ),
+        i = Math.abs(ψ),
+        sinhψ = .5 * ((sinhψ = Math.exp(i)) - 1 / sinhψ);
+    if (r) {
+      var cscφ = 1 / Math.sin(r),
+          cotφ2 = (cotφ2 = Math.cos(r) * cscφ) * cotφ2,
+          b = -(cotφ2 + m * (sinhψ * sinhψ * cscφ * cscφ + 1) - 1),
+          cotλ2 = .5 * (-b + Math.sqrt(b * b - 4 * (m - 1) * cotφ2));
+      return [
+        ellipticF(Math.atan(1 / Math.sqrt(cotλ2)), m) * sgn(φ),
+        ellipticF(Math.atan(Math.sqrt(Math.max(0, cotλ2 / cotφ2 - 1) / m)), 1 - m) * sgn(ψ)
+      ];
+    }
+    return [
+      0,
+      ellipticF(Math.atan(sinhψ), 1 - m) * sgn(ψ)
+    ];
+  }
+
+  // Calculate F(φ|m) where m = k² = sin²α.
+  // See Abramowitz and Stegun, 17.6.7.
+  function ellipticF(φ, m) {
+    var a = 1,
+        b = Math.sqrt(1 - m),
+        c = Math.sqrt(m);
+    for (var i = 0; Math.abs(c) > ε; i++) {
+      if (φ % π) {
+        var dφ = Math.atan(b * Math.tan(φ) / a);
+        if (dφ < 0) dφ += π;
+        φ += dφ + ~~(φ / π) * π;
+      } else φ += φ;
+      c = (a + b) / 2;
+      b = Math.sqrt(a * b);
+      c = ((a = c) - b) / 2;
+    }
+    return φ / (Math.pow(2, i) * a);
+  }
+
+//---  
+
+
+function dymaxion(λ, φ) {
+    // var α = Math.sqrt(4 - 3 * Math.sin(Math.abs(φ)));
+    // return [
+    //   2 / Math.sqrt(6 * π) * λ * α,
+    //   sgn(φ) * Math.sqrt(2 * π / 3) * (2 - α)
+    // ];
+    var cartesian = convert_s_t_p(λ,φ);
+    return [
+    cartesian.x
+    ,
+    cartesian.y
+    ]
+  }
+
+
+//---
+
+
+  var projection = d3.geo.projection,
+      projectionMutator = d3.geo.projectionMutator;
+
+  d3.geo.dymaxion = function() { return projection(dymaxion); };
+
+})();
diff --git a/index.html b/index.html
@@ -12,15 +12,15 @@
   stroke-width: 1.5px;
 }
 
-.graticule {
+/*.graticule {
   fill: none;
-  stroke: #fff;
+  stroke: #000;
   stroke-width: .5px;
 }
 
 .graticule:nth-child(2n) {
   stroke-dasharray: 2,2;
-}
+}*/
 
 .land {
   fill: #d7c7ad;
@@ -36,13 +36,16 @@
 <body>
 <script src="http://d3js.org/d3.v3.min.js"></script>
 <script src="https://raw.github.com/d3/d3-plugins/master/geo/projection/projection.js"></script>
+<script src="/protovis_dymaxion.js"></script>
+<script src="/d3_dymaxion.js"></script>
+
 <script>
 
 var width = 960,
     height = 500;
 
-var projection = d3.geo.hammer()
+var projection = d3.geo.dymaxion()
-    .translate([width / 2 - .5, height / 2 - .5]);
+    .translate([width/3, height/6])
 
 var path = d3.geo.path()
     .projection(projection);
@@ -53,30 +56,30 @@
     .attr("width", width)
     .attr("height", height);
 
-svg.append("path")
+// svg.append("path")
-    .datum(graticule.outline)
+//     .datum(graticule.outline)
-    .attr("class", "background")
+//     .attr("class", "background")
-    .attr("d", path);
+//     .attr("d", path);
 
-svg.selectAll(".graticule")
+// svg.selectAll(".graticule")
-    .data(graticule.lines)
+//     .data(graticule.lines)
-  .enter().append("path")
+//   .enter().append("path")
-    .attr("class", "graticule")
+//     .attr("class", "graticule")
-    .attr("d", path);
+//     .attr("d", path);
 
-svg.append("path")
+// svg.append("path")
-    .datum(graticule.outline)
+//     .datum(graticule.outline)
-    .attr("class", "foreground")
+//     .attr("class", "foreground")
-    .attr("d", path);
+//     .attr("d", path);
 
-d3.json("/d/3682676/readme-boundaries.json", function(error, collection) {
+d3.json("boundaries.json", function(error, collection) {
   svg.insert("path", ".graticule")
       .datum(collection)
       .attr("class", "boundary")
       .attr("d", path);
 });
 
-d3.json("/d/3682676/readme-land.json", function(error, collection) {
+d3.json("/land.json", function(error, collection) {
   svg.insert("path", ".graticule,.boundary")
       .datum(collection)
       .attr("class", "land")