如何在世界 scatter_geo 上仅显示美国的国家边界？

Question

我正在使用plotly express scattergeo 在世界地图上绘制一些数据点，但我也希望在州一级显示美国。当我使用时，“世界”范围不会对 showsubunits 参数做出反应：

fig.update_geos(
   visible=True, resolution=50, scope='world',
   showcountires=True, countrycolor="Black",
   showsubunits=True, subunitcolor='Brown')

但是，当我将范围仅限于美国时，“美国”范围会做出反应并显示州界线。 有没有一种方法可以以某种方式覆盖仅包含国家边界的世界地图和包含国家边界的美国地图，以获取同一组散布纬度/经度数据？作为一个更普遍的问题，我可以有选择地定义在州一级显示哪些国家/地区吗？

Answer 1

• 您可以使用 https://github.com/nvkelso/natural-earth-vector/blob/master/geojson/ne_110m_admin_1_states_provinces_lines.geojson
• 然后是添加额外跟踪以从 geojson 中提取州边界详细信息的情况根据需要纬度和经度

import requests
import numpy as np
import plotly.express as px
import plotly.graph_objects as go

states_geojson = requests.get(
    "https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_110m_admin_1_states_provinces_lines.geojson"
).json()


df = px.data.gapminder().query("year == 2007")
fig = px.scatter_geo(
    df,
    locations="iso_alpha",
    size="pop",  # size of markers, "pop" is one of the columns of gapminder
)


fig = fig.add_trace(
    go.Scattergeo(
        lat=[
            v
            for sub in [
                np.array(f["geometry"]["coordinates"])[:, 1].tolist() + [None]
                for f in states_geojson["features"]
            ]
            for v in sub
        ],
        lon=[
            v
            for sub in [
                np.array(f["geometry"]["coordinates"])[:, 0].tolist() + [None]
                for f in states_geojson["features"]
            ]
            for v in sub
        ],
        line_color="brown",
        line_width=1,
        mode="lines",
        showlegend=False,
    )
)

fig.update_geos(
    visible=True, resolution=50, scope="world", showcountries=True, countrycolor="Black"
)