Yahoo!地図APIで円を描きたい(JavaScript版地図「旧地図Web API」対応版)

// 測地系どうしよう？
var map = new YahooMapsCtrl("map", new YLLPoint(35, 135), 4, YMapMode.MAP, YDatumType.TOKYO97);
// var map = new YahooMapsCtrl("map", new YLLPoint(35, 135), 4, YMapMode.MAP, YDatumType.WGS84);

// 中心地
var point = new YLLPoint(35, 135);

// 描く円の半径(km)
var radius = 1;

// 赤道半径(m)
var EquatorialRadius;

// 扁平率の逆数 : 1/f
var F;

if (map.getDatum() == YDatumType.TOKYO97) {
  EquatorialRadius = 6377397.155;
  F = 299.152813;
} else if (map.getDatum() == YDatumType.WGS84){
  EquatorialRadius = 6378137;
  F = 298.257223;
} else {
  alert('unknown datum');
}

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

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

// 緯度のラジアン表記
var pLatRadians = (point.lat * Math.PI) / 180;

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

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

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

// 半径をｍ単位に
var R = radius * 1000;

// 0～360°で1°ずつ頂点を求めることで360角形を生成
var points = new Array(360);
for (i = 0; i <= 360; i++) {
  var rad = i / 180 * Math.PI;
  var lat = (R / arc_lat) * Math.sin(rad) + point.lat;
  var lon = (R / arc_lng) * Math.cos(rad) + point.lon;
  points[i] = new YLLPoint(lat, lon);
}

// id
var circle_id = 'circle-id';
// popup
var circle_popup ='';
// style
var circle_style = { strokeWidth: 1, strokeColor: 'FF0000', strokeOpacity: 1 };

map.addPolyline(circle_id, points, circle_popup, circle_style);