d3.js - GPX は d3.js で RGF93 を追跡します

Question

このブログで説明されているように、GPX トラックが挿入された 3D サーフェスを実装しています: http://blog.mastermaps.com/2013/11/showing-gps-tracks-in-3d-with-threejs.html

私の 3D マップは、RGF93 投影法を使用して QGIS で作成されたカスタム マップであり、d3.js で適切な投影法を見つけることができません。

私の地図の gdalinfo は

["RGF93 / Lambert-93",
GEOGCS["RGF93",
    DATUM["Reseau_Geodesique_Francais_1993",
    SPHEROID["GRS 1980",6378137,298.257222101,
            AUTHORITY["EPSG","7019"]],
        TOWGS84[0,0,0,0,0,0,0],
        AUTHORITY["EPSG","6171"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9122"]],
    AUTHORITY["EPSG","4171"]],
PROJECTION["Lambert_Conformal_Conic_2SP"],
PARAMETER["standard_parallel_1",49],
PARAMETER["standard_parallel_2",44],
PARAMETER["latitude_of_origin",46.5],
PARAMETER["central_meridian",3],
PARAMETER["false_easting",700000],
PARAMETER["false_northing",6600000],
UNIT["metre",1,
    AUTHORITY["EPSG","9001"]],
AXIS["X",EAST],
AXIS["Y",NORTH],
AUTHORITY["EPSG","2154"]]
Origin = (585097.372746756765991,6171321.216855648905039)
Pixel Size = (15.149120309400001,-14.887674158699999)
Metadata:
    AREA_OR_POINT=Area
Image Structure Metadata:
    INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  585097.373, 6171321.217) (  1d36' 3.31"E, 42d37'56.65"N)
Lower Left  (  585097.373, 6165946.766) (  1d36' 7.48"E, 42d35' 2.74"N)
Upper Right (  590475.310, 6171321.217) (  1d39'59.00"E, 42d37'59.67"N)
Lower Right (  590475.310, 6165946.766) (  1d40' 2.98"E, 42d35' 5.75"N)
Center      (  587786.342, 6168633.992) (  1d38' 3.19"E, 42d36'31.22"N)
Band 1 Block=355x5 Type=Float32, ColorInterp=Gray
    Computed Min/Max=1735.500,2908.485
NoData Value=-3.4028234663852886e+38

そして、私が（他の多くの組み合わせの中で）試みている投影ですが、成功していません。scale(1) のときはポイントがサーフェスの中心に配置されていますが、スケール値を大きくすると遠くに移動することがわかります。

var projection = d3.geoConicConformal()
    .rotate([-3, 0])
    .center([1.6342194, 42.6086722])
    .parallels([49, 44])
    .translate([100, 100])
    .scale(terrain.map);

私は何が間違っているのか、私の投影に適切な値は何かを知りたいと思います。

---

編集

私はほとんど解決策を見つけました：

var projection = d3.geoConicConformal()
    .rotate([-terrain.center[1], 0])
    .center([0, terrain.center[0]])
    .parallels([44, 49])
    .translate([100, 100])
    .scale(terrain.map);

ここで、terrain.center[lat, lon] は gdalinfo のコーナー座標: Center からの値です。

完全に正確でなくても、ポイントは予想される位置とほぼ一致します。少し反時計回りに回転しているように見えます。回転値で遊んでみましたが、すべての投影を左または右に変換するだけです。