如何在GeoViews中添加底图瓦片时指定地图范围

问题核心是投影不匹配：你使用的瓦片图层（如OSM）采用Web Mercator投影（EPSG:3857，单位为米），而你的点数据是WGS84（EPSG:4326，单位为十进制度数）。直接用度数设置xlim会被解析为米单位，导致显示范围完全错误。

以下是两种可行的解决方式：

方式一：转换点数据到Web Mercator投影 #

先将GeoDataFrame的坐标转换为瓦片使用的EPSG:3857，再设置对应范围：

import geopandas as gpd
import geoviews as gv
import numpy as np
from pyproj import Transformer
gv.extension('bokeh')

# 生成原始WGS84坐标的点数据
gdf = gpd.GeoDataFrame(
    geometry=gpd.points_from_xy(np.linspace(-88, -82, 100), np.linspace(32,38, 100)), 
    crs='WGS84'
)

# 转换到Web Mercator投影（EPSG:3857）
gdf_3857 = gdf.to_crs('EPSG:3857')

# 创建Points并叠加瓦片
points = gv.Points(data=gdf_3857)
points_map = gv.tile_sources.OSM * points

# 将目标经度范围转换为3857坐标系的米单位
transformer = Transformer.from_crs('EPSG:4326', 'EPSG:3857')
x_min, _ = transformer.transform(0, -86)
x_max, _ = transformer.transform(0, -84)

points_map = points_map.opts(xlim=(x_min, x_max))
points_map

方式二：让GeoViews自动处理投影，用地理坐标设置范围 #

通过指定projection和x_range/y_range，让GeoViews自动完成投影转换，直接用十进制度数设置范围：

import geopandas as gpd
import geoviews as gv
import numpy as np
gv.extension('bokeh')

gdf = gpd.GeoDataFrame(
    geometry=gpd.points_from_xy(np.linspace(-88, -82, 100), np.linspace(32,38, 100)), 
    crs='WGS84'
)

# 创建Points时指定原始CRS
points = gv.Points(data=gdf).opts(crs='EPSG:4326')
points_map = gv.tile_sources.OSM * points

# 设置投影和地理坐标范围
points_map = points_map.opts(
    projection=gv.CRS('EPSG:3857'),
    x_range=(-86, -84),
    y_range=(33, 37)  # 可根据需求调整纬度范围
)
points_map

第二种方式更直观，无需手动转换坐标，GeoViews会自动处理WGS84到Web Mercator的投影转换，直接用度数设置范围即可。

内容的提问来源于stack exchange，提问作者hm8