Question

我试图将Open Street Maps纬度/经度转换为像素位置。我在java中发现了这个example。当我将代码转换为Delphi时，X结果是正确的，但Y结果稍微有点。

Java代码：

double lon = -89.08805;
double lat = 30.6393;
double zoom = 6; // 6.5156731549786215 would be possible too

double lon_rad = Math.toRadians(lon);
double lat_rad = Math.toRadians(lat);
double n = Math.pow(2.0, zoom);

double tileX = ((lon + 180) / 360) * n;
double tileY = (1 - (Math.log(Math.tan(lat_rad) + 1.0/Math.cos(lat_rad)) / Math.PI)) * n / 2.0;

System.out.println("tileX = "+tileX+" tileY = "+tileY);

德尔福代码：

function LatLngToTilePixels(Lat, Lng: Double; ZoomLevel: Integer): TPointF;
var
  lon_rad, lat_rad, n: Double;
  TileX, TileY: Double;
begin
  lon_rad := DegToRad(Lng);
  lat_rad := DegToRad(Lat);

  n := Power(2.0, ZoomLevel);

  Result.X := ((Lng + 180) / 360) * n;
  Result.Y := (1 - (Math.Log10(Math.tan(lat_rad) + 1.0/cos(lat_rad)) / PI)) * n / 2.0;
end;

结果，使用示例中给出的输入（lat = 30.6393，lng = -89.08805，zoom = 6）：

JAVA：

tileX = 16.162124444444444

tileY = 26.273150713795616

的Delphi：

tileX = 16.1621244444444

tileY = 29.5128609563099＆lt; ---错误的结果

我觉得问题可能是math.log调用。 JAVA使用日志，但Delphi使用log10。我试过log2，但结果更糟。

Answer 1

您使用不同基数的对数：

Java：基础e，即自然对数，请参阅Math.log(double)

德尔福：显然基础10

您可以使用以下公式计算不受支持的基数的对数：

log_a(b) = log_x(b) / log_x(a)

AFAIK在delphi中有一个Ln函数，它计算自然对数，至少列出here。

JAVA和Delphi中的数学函数给出了不同的结果

1 个答案: