经纬度5位数和6位数差多少_经纬度坐标小数位与精度的对应关系

本文通过测试展示了经纬度坐标的小数位数与精度的对应关系,例如:经度相差1位小数点后6位对应0.1米,5位对应1米。还探讨了全国经纬度数据保存时,3位小数点精度约为100米,可能需要保存约10亿条记录。
摘要由CSDN通过智能技术生成

纬度不变 经度相差1的情况

测试结果如下:

小数点后6位,精度:0.1米

0.10139558886691133

小数点后5位,精度:1米

1.0139558860891593

小数点后4位,精度:10米

10.139558845411655

小数点后3位,精度:101米

101.39558846680167

小数点后3位(经纬度同时改变1),精度:150米

150.4837860345681

小数点后2位,精度:1014米

1013.9558844430293

测试代码:

public class DistanceCalcUtil {

// 地球半径

private static final double EARTH_RADIUS = 6370996.81;

// 弧度

private static double radian(double d) {

return d * Math.PI / 180.0;

}

/**

* 计算两点之间的距离

*

* @param lat1

* @param lng1

* @param lat2

* @param lng2

* @return

*/

public static double distanceOfTwoPoints(double l

  • 1
    点赞
  • 8
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
要根据经纬度获取h5文件中的像元坐标,你可以使用一些地理信息处理库来实现。下面是一个示例代码,使用了`pyproj`库来进行坐标转换和投影操作: ```python import pandas as pd import h5py from pyproj import Proj, transform # 读取Excel表格中的经纬度数据 df = pd.read_excel('data.xlsx') # 假设Excel文件名为data.xlsx,需要与代码在同一目录下 latitudes = df['纬度列名'] longitudes = df['经度列名'] # 定义投影坐标系 in_proj = Proj(init='epsg:4326') # 输入经纬度坐标系,这里使用WGS84坐标系 out_proj = Proj(init='epsg:XXXX') # 输出坐标系,根据你的具体数据情况指定,例如UTM投影坐标系 # 打开h5文件 h5_file = h5py.File('data.h5', 'r') # 假设h5文件名为data.h5,需要与代码在同一目录下 # 提取像元值 pixel_values = [] for lat, lon in zip(latitudes, longitudes): # 经纬度转换为目标坐标系的坐标 x, y = transform(in_proj, out_proj, lon, lat) # 将坐标四舍五入为整数,表示像元坐标 pixel_coords = int(round(x)), int(round(y)) # 从h5文件中读取像元值 pixel_value = h5_file['dataset_name'][pixel_coords[1], pixel_coords[0]] # dataset_name是h5文件中的数据集名字 pixel_values.append(pixel_value) # 关闭h5文件 h5_file.close() # 打印提取到的像元值 print(pixel_values) ``` 上述代码中的`'epsg:XXXX'`需要根据你的实际数据情况指定输出坐标系,例如使用UTM投影坐标系。你可以根据需要修改该参数。 希望这能帮助到你!如果还有其他问题,请随时提问。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值