Python地理位置信息库geopy的使用:根据中心点坐标,方向,距离计算坐标; 利用两点经纬度计算地理空间距离

Python地理位置信息库geopy的使用:根据中心点坐标,方向,距离计算坐标

import geopy.distance

def get_distance_point(lat, lon, distance, direction):
    """
    根据经纬度,距离,方向获得一个地点
    :param lat: 纬度
    :param lon: 经度
    :param distance: 距离(千米)
    :param direction: 方向(北:0,东:90,南:180,西:360)
    :return:
    """
    start = geopy.Point(lat, lon)
    d = geopy.distance.VincentyDistance(kilometers=distance)
    return d.destination(point=start, bearing=direction)


p = get_distance_point(34.490194, 117.281568, 0.5, 0)
print(p.latitude, p.longitude)

https://www.cnblogs.com/lloyd-zh/p/9805828.html

利用两点经纬度计算地理空间距离

from math import sin, asin, cos, radians, fabs, sqrt
 
EARTH_RADIUS = 6371.137      # 地球平均半径大约6371km
 
 
def hav(theta):
    s = sin(theta / 2)
    return s * s
 
 
def get_distance_hav(lat0, lng0, lat1, lng1):
    # 用haversine公式计算球面两点间的距离
    # 经纬度转换成弧度
    lat0 = radians(lat0)
    lat1 = radians(lat1)
    lng0 = radians(lng0)
    lng1 = radians(lng1)
    dlng = fabs(lng0 - lng1)
    dlat = fabs(lat0 - lat1)
    h = hav(dlat) + cos(lat0) * cos(lat1) * hav(dlng)
    distance = 2 * EARTH_RADIUS * asin(sqrt(h))      # km
    return distance
 
if __name__ == "__main__":
    result = get_distance_hav(39.995304, 116.308264, 40.003304, 116.326759)
    print(result)

https://blog.csdn.net/chinawangfei/article/details/125458095

  • 0
    点赞
  • 5
    收藏
    觉得还不错? 一键收藏
  • 1
    评论
获取高分一号影像中心点经纬度坐标可以通过以下步骤实现: 1. 首先,需要获取高分一号影像的元数据信息,其中包括影像的行列数、像元分辨率、中心点坐标等。可以使用GDAL库读取影像元数据信息。 2. 然后,根据影像的行列数、像元分辨率和中心点坐标计算出影像的四个角点坐标。 3. 最后,利用GIS软件或者在线地图服务,将四个角点坐标输入到地图中,即可获取高分一号影像中心点经纬度坐标。 以下是示例代码: ```python import gdal import numpy as np # 读取影像元数据信息 dataset = gdal.Open('GF1.tif') cols = dataset.RasterXSize rows = dataset.RasterYSize transform = dataset.GetGeoTransform() # 计算影像四个角点坐标 ulx = transform[0] uly = transform[3] lrx = ulx + cols * transform[1] lry = uly + rows * transform[5] # 计算影像中心点坐标 center_x = (ulx + lrx) / 2 center_y = (uly + lry) / 2 # 将中心点坐标转换为经纬度坐标 driver = gdal.GetDriverByName('MEM') ds = driver.Create('', 1, 1, 1, gdal.GDT_Float32) ds.SetGeoTransform((center_x, transform[1], transform[2], center_y, transform[4], transform[5])) srs = gdal.osr.SpatialReference() srs.ImportFromWkt(dataset.GetProjection()) ds.SetProjection(srs.ExportToWkt()) res = gdal.ReprojectImage(dataset, ds, dataset.GetProjection(), srs.ExportToWkt(), gdal.GRA_NearestNeighbour) center_lon, center_lat, _ = np.array(res.GetRasterBand(1).ReadAsArray()).ravel() print('高分一号影像中心点经度为:', center_lon) print('高分一号影像中心点纬度为:', center_lat) ``` 需要注意的是,以上代码中的`GF1.tif`是高分一号影像的文件路径,需要根据实际情况进行修改。此外,代码中使用了GDAL库进行影像处理。如果没有安装GDAL库,可以使用以下命令进行安装: ```bash pip install gdal ```
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值