8

我有一点问题。我正在尝试使用 Javascript 和 Google Maps API v2 执行以下操作:

myhurricane.net - Wind Radii Profile

我可以使用在互联网上找到的公式很好地绘制单个圆圈。我面临的问题是圈子必须:

A. 同心,并且 B. 每个“象限”必须有不同的半径,即 NE、NW、SE 和 SW

我几乎在互联网上能想到的所有地方都进行了搜索,但没有想出如何做到这一点。很明显,以前有人这样做过,因此我为什么要在程序员论坛上提问。:)

谢谢!

更新:我已经使用以下代码绘制了我认为每个点的坐标。对于下图:

快照 1252125257.781397

这是使用以下 JS 获得的:

http://gist.github.com/181290

注意:此 javascript 来自(稍作修改)以下站点,该站点可能包含更多关于算法最终可能的答案:http ://www.movable-type.co.uk/scripts/latlong.html

更新2:我能够在谷歌地图中得到这个:

同心圆进展

使用以下代码创建:

var NEQ = [0, 90];
var SEQ = [90, 180];
var SWQ = [180, 270];
var NWQ = [270, 0];

// var centrePoint = new LatLon(25.0, -83.1);
// pointsForWindQuadrant(NEQ, centrePoint, 50);
function pointsForWindQuadrant(quadrantDegrees, centrePoint, radius){
  var points = [];

  // Points must be pushed into the array in order
  points.push(new google.maps.LatLng(centrePoint.lat, centrePoint.lon));

  for(i = quadrantDegrees[0]; i <= quadrantDegrees[1]; i++){
    var point = centrePoint.destPoint(i, radius * 1.85);
    points.push(new google.maps.LatLng(point.lat, point.lon)); // Radius should be in nautical miles from NHC
  }

  points.push(new google.maps.LatLng(centrePoint.lat, centrePoint.lon));

  return points;
}

更新 3:我可能还应该指出,这是一个地理坐标系(因为整个事情都是针对热带气旋风半径),而不是笛卡尔坐标系。谢谢!

4

2 回答 2

6

基本上,将圆计算为 x,y = (cos(a), sin(a)),然后将这个(两项)乘以作为角度的适当函数的半径。我不太了解 Javascript 或 Google 地图,所以我将在 Python 中执行此操作,希望从中可以清楚地了解。

from pylab import *

def Rscale(a):
    if a>3*pi/2:  # lower right, and then work CW around the circle
        return 1.
    elif a>pi:  # lower left
        return .9
    elif a>pi/2:   # upper left
        return .8
    else:       # upper right
        return 1.

def step_circle(R):
    return array([(R*Rscale(a))*array([cos(a), sin(a)]) for a in arange(0, 2*pi, .001)])

for R in (.5, .7, .9):  # make three concentric circles
    c = step_circle(R)
    plot(c[:,0], c[:,1])

show()

这使 替代文字

我无法真正按照您的草图进行操作,所以我只是猜测了数字。另外,我将最右边的两个象限设置为相同,因为这就是您的情节,但这当然是可选的。

于 2009-09-05T05:09:45.773 回答
4

我想到了。这是最终的代码。也许它可以重构一下?

// Returns points for a wind field for a cyclone. Requires
// a LatLon centre point, and an array of wind radii, starting
// from the northeast quadrant (NEQ), i.e., [200, 200, 150, 175]
//
// Returns points to be used in a GPolyline object.
function pointsForWindQuadrant(centrePoint, radii){
  if(radii.length != 4){ return false; }

  var points = [];
  var angles = [0, 90, 180, 270];

  // For each angle 0, 90, 180, 270...
  for(a = 0; a < angles.length; a++){
    // For each individual angle within the range, create a point...
    for(i = angles[a]; i <= angles[a] + 90; i++){
      var point = centrePoint.destPoint(i, radii[a] * 1.85); // Radius should be in nautical miles from NHC
      points.push(new google.maps.LatLng(point.lat, point.lon));
    }
  }

  // Add the first point again, to be able to close the GPolyline
  var point = centrePoint.destPoint(0, radii[0] * 1.85);
  points.push(new google.maps.LatLng(point.lat, point.lon));

  return points;
}

这导致以下结果:

新 myhurricane.net - 风半径(地图视图) 新 myhurricane.net - 风半径(卫星视图)

于 2009-09-05T07:42:01.270 回答