4字节 经纬度_以一个经纬度为中心计算出四个顶点

public const double EARTH_RADIUS = 6378137.0;//地球半径(m)

///

/// 以一个经纬度为中心计算出四个顶点

///

/// 半径(米)

///

public static Degree[] GetDegreeCoordinates(Degree Degree1, double distance)

{

double dlng = 2 * Math.Asin(Math.Sin(distance / (2 * EARTH_RADIUS)) / Math.Cos(Degree1.X));

dlng = degrees(dlng);//一定转换成角度数 原PHP文章这个地方说的不清楚根本不正确 后来lz又查了很多资料终于搞定了

double dlat = distance / EARTH_RADIUS;

dlat = degrees(dlat);//一定转换成角度数

return new Degree[] { new Degree(Math.Round(Degree1.X + dlat,6), Math.Round(Degree1.Y - dlng,6)),//left-top

new Degree(Math.Round(Degree1.X - dlat,6), Math.Round(Degree1.Y - dlng,6)),//left-bottom

new Degree(Math.Round(Degree1.X + dlat,6), Math.Round(Degree1.Y + dlng,6)),//right-top

new Degree(Math.Round(Degree1.X - dlat,6), Math.Round(Degree1.Y + dlng,6)) //right-bottom

};

}

///

/// 角度数转换为弧度公式

///

///

///

private static double radians(double d)

{

return d * Math.PI / 180.0;

}

///

/// 弧度转换为角度数公式

///

///

///

private static double degrees(double d)

{

return d * (180 / Math.PI);

}

///

/// 代表经度, 纬度

///

public class Poin

{

/// 纬度 X

/// 经度 Y

public Poin(double lat, double lng)

{

this.lat = lat;

this.lng = lng;

}

// 纬度 X

private double lat;

// 经度 Y

private double lng;

///

/// 代表纬度 X轴

///

public double Lat { set; get; }

///

/// 代表经度 Y轴

///

public double Lng { get; set; }

public double RadLat { get { return lat * Math.PI / 180; } }

public double RadLng { get { return lng * Math.PI / 180; } }

}

///

/// 经纬度坐标

///

public class Degree

{

public Degree(double x, double y)

{

X = x;

Y = y;

}

private double x;

public double X

{

get { return x; }

set { x = value; }

}

private double y;

public double Y

{

get { return y; }

set { y = value; }

}

}

///

/// 计算坐标点的距离

///

public void GetDistance()

{

//x经度 y纬度

string enbid = "";

double distance = 0;

double x1 = ConverToRadian(31.499605);

double x2 = ConverToRadian(31.4996);

double y1 = ConverToRadian(104.73775);

double y2 = ConverToRadian(104.73783);

double d = 0;

if (Math.Sin(x1) * Math.Sin(x2) + Math.Cos(x1) * Math.Cos(x2) * Math.Cos(y1 - y2) > 1)

d = EARTH_RADIUS * Math.Acos(1);

else if (Math.Sin(x1) * Math.Sin(x2) + Math.Cos(x1) * Math.Cos(x2) * Math.Cos(y1 - y2) < -1)

d = EARTH_RADIUS * Math.Acos(-1);

else

d = EARTH_RADIUS * Math.Acos(Math.Sin(x1) * Math.Sin(x2) + Math.Cos(x1) * Math.Cos(x2) * Math.Cos(y1 - y2));

}