python地球重量计算_python 计算曼哈顿距离

def clip(self, n, min, max):

if n < min:

n = min

if n > max:

n = max

return n

def longlat2tile(self, x1, y1, level):

M_PI = 3.1415926

EarthRadiusInMeters = 6378137

PixelsPerTile = 256

EarthCircumferenceInMeters = 2 * M_PI * EarthRadiusInMeters

rad_long = x1 * M_PI / 180

rad_lat = y1 * M_PI / 180

x_meters = EarthRadiusInMeters * rad_long

sin_latitude = math.sin(rad_lat)

i_log = math.log((1 + sin_latitude) / (1 - sin_latitude))

y_meters = EarthRadiusInMeters / 2 * i_log

num_pixels = long(PixelsPerTile << level)

meters_per_pixel = EarthCircumferenceInMeters / num_pixels

tile_long = int(self.clip((EarthCircumferenceInMeters / 2 + x_meters) / meters_per_pixel + 0.5, 0, num_pixels - 1))

tmp = long((EarthCircumferenceInMeters / 2 - y_meters))

tile_lat = int(self.clip(tmp / meters_per_pixel + 0.5, 0, num_pixels - 1))

return tile_long, tile_lat

计算直线距离

def getDistance(self, x1,y1,x2,y2):

dx = x2 - x1

dy = y2 - y1

sx = math.cos(x1* 0.01745329252)

dis = math.sqrt(dx*dx*sx*sx + dy*dy) * 111195.0

return dis