0

我想使用品脱将度数（地理 CRS 中的距离）转换为海里。

https://geopandas.org/docs/reference/api/geopandas.GeoDataFrame.sjoin_nearest.html以 epsg:4326 为单位输出距离。

给定距离（以纳米为单位）因赤道到两极而异，我不确定这是否可能。

我可以使用 1 度 ~= 111 公里 ~= 60 海里的经验法则。

也许它可以使用起点和距离来计算，例如：https ://github.com/anitagraser/movingpandas/blob/master/movingpandas/geometry_utils.py#L38

此代码也很有用：https ://geopy.readthedocs.io/en/stable/#module-geopy.distance

这是一些要测试的代码：

import pandas as pd
import geopandas as gpd

df = pd.DataFrame({"lon": [0], "lat": [0]})
gdf_pt = gpd.GeoDataFrame(df, geometry=gpd.points_from_xy(df["lon"], df["lat"]), crs="epsg:4326")

df2 = pd.DataFrame({"lon": [1, 2], "lat": [0, 0]})
gdf_pts = gpd.GeoDataFrame(df2, geometry=gpd.points_from_xy(df2["lon"], df2["lat"]), crs="epsg:4326")

value = gdf_pt.sjoin_nearest(gdf_pts, distance_col="distances")["distances"].values[0]


import pint

l = value * ureg.arcdegree
4

2 回答 2

0

可能最好把它扔给墨卡托，如果可以的话就用它

import pint_pandas

gdf = gdf_pt.to_crs("EPSG:3395").sjoin_nearest(gdf_pts.to_crs("EPSG:3395"), distance_col="distances")
gdf["distance"] = gdf["distance"].astype("pint[meter]").pint.to("nautical_mile")
于 2021-10-28T02:06:37.567 回答
0

我从现有代码中提取的这个函数计算两个纬度/经度集之间的距离（以米为单位）。“rlat”和“rlong”以弧度表示；您必须从度数进行转换。要获取 nm 而不是米，只需将 R 设置为 3440。

from math import *

# Radius of the earth, in meters.

R = 6371000

# Return distance between two lat/longs.

def distance( pt1, pt2 ):
    rlat1 = pt1.rlat
    rlat2 = pt2.rlat
    dlat = pt2.rlat - pt1.rlat
    dlong = pt2.rlong - pt1.rlong

    a = sin(dlat/2) * sin(dlat/2) + cos(rlat1) * cos(rlat2) * sin(dlong/2) * sin(dlong/2)
    c = 2 * atan2(sqrt(a), sqrt(1-a))
    return R * c
于 2021-10-26T19:06:49.753 回答