LeafletでL.circleを使わずに円を描く

/*
  引数はL.Circle()に準拠
  なお、options.radius は "m"単位
*/
function pseudoCircle(latlng, options) {

  // どんな形で来るのかわからんので
  var pos = L.latLng(latlng);

  // 360角形で近似
  var vertex = 360;

  // 赤道半径(m) (WGS-84)
  var EquatorialRadius = 6378137;

  // 扁平率の逆数 : 1/f (WGS-84)
  var F = 298.257223563;

  // 離心率の２乗
  var E = ((2 * F) -1) / Math.pow(F, 2);

  // 赤道半径 × π
  var PI_ER = Math.PI * EquatorialRadius;

  // 1 - e^2 sin^2 (θ)
  var TMP = 1 - E * Math.pow(pos.lat * Math.PI / 180, 2);

  // 経度１度あたりの長さ(m)
  var arc_lat = (PI_ER * (1 - E)) / (180 * Math.pow(TMP, 3/2)); 

  // 緯度１度あたりの長さ(m)
  var arc_lng = (PI_ER * Math.cos(pos.lat * Math.PI / 180)) / (180 * Math.pow(TMP, 1/2)); 

  var points = new Array();
  for (var i = 0; i <= vertex; i++) {
    var rad = (i / (vertex / 2)) * Math.PI;
    var lat = (options.radius / arc_lat) * Math.sin(rad) + pos.lat;
    var lng = (options.radius / arc_lng) * Math.cos(rad) + pos.lng;
    points[i] = [lat, lng];
  }
  return L.polygon(points, options);
}