平行四边形区域的垂直剖面该用什么画呢？

一大碗年糕

想请教一下大家，Python画左图里绿色平行四边形的垂直剖面环流该用什么比较方便呀？😂metpy的cross section好像只能画任意两点之间的剖面2025.9.20更新，添加平行四边形

1. from matplotlib.patches import Polygon as MplPolygon
   
2. import numpy as np
   
3. import math
   
4. import cartopy.crs as ccrs
   
5. 
   
6. def add_parallelogram(ax, center_lon, center_lat, width_deg, height_deg, skew_deg=0.3,
   
7.                       angle_deg=0, edgecolor='blue', linewidth=2, linestyle='-',
   
8.                       zorder=2, label=None):
   
9.     """
   
10.     在 cartopy 地理坐标中添加一个平行四边形。
    
11.    
    
12.     参数：
    
13.     - ax: cartopy axes
    
14.     - center_lon, center_lat: 中心经纬度
    
15.     - width_deg, height_deg: 平行四边形宽高（单位：度）
    
16.     - skew_deg: 水平方向的偏移量（控制平行四边形形状）
    
17.     - angle_deg: 绕中心顺时针旋转角度（度）
    
18.     - edgecolor, linewidth, linestyle: 样式
    
19.     """
    
20.     angle_rad = math.radians(angle_deg)
    
21.     dx, dy = width_deg / 2, height_deg / 2
    
22. 
    
23.     # 定义平行四边形的四个角
    
24.     corners = [(-dx, -dy), (dx, -dy), (dx + skew_deg, dy), (-dx + skew_deg, dy)]
    
25.    
    
26.     # 旋转 + 平移
    
27.     rotated_corners = []
    
28.     for x, y in corners:
    
29.         x_rot = x * math.cos(angle_rad) - y * math.sin(angle_rad) + center_lon
    
30.         y_rot = x * math.sin(angle_rad) + y * math.cos(angle_rad) + center_lat
    
31.         rotated_corners.append((x_rot, y_rot))
    
32.    
    
33.     # 绘制
    
34.     poly = MplPolygon(rotated_corners, closed=True, edgecolor=edgecolor,
    
35.                       facecolor='none', linewidth=linewidth, linestyle=linestyle,
    
36.                       transform=ccrs.PlateCarree(), zorder=zorder, label=label)
    
37.     ax.add_patch(poly)
    
38.     return poly
    

复制代码

一大碗年糕

用rioxarray可以实现

tosaka

一大碗年糕 发表于 2025-6-23 12:58
用rioxarray可以实现

大佬用的哪个函数实现的

一大碗年糕

tosaka 发表于 2025-7-31 09:45
大佬用的哪个函数实现的

自己定义区域，然后用rixarray的clip截取

1. import rioxarray as rio
   
2. def get_rotated_rectangle_polygon(center_lon, center_lat, width_deg, height_deg, angle_deg):
   
3.     angle_rad = math.radians(angle_deg)
   
4.     dx, dy = width_deg / 2, height_deg / 2
   
5.     offsets = [(-dx, -dy), (dx, -dy), (dx, dy), (-dx, dy)]
   
6.     corners = []
   
7.     for ox, oy in offsets:
   
8.         x_rot = ox * math.cos(angle_rad) - oy * math.sin(angle_rad) + center_lon
   
9.         y_rot = ox * math.sin(angle_rad) + oy * math.cos(angle_rad) + center_lat
   
10.         corners.append((x_rot, y_rot))
    
11.     return Polygon(corners)
    
12. 
    
13. AR = AR.rio.write_crs("EPSG:4326")  # 明确是地理坐标
    
14. 
    
15. polygon_B = get_rotated_rectangle_polygon(127, 38, 4, 8, 25)  # 朝鲜半岛区域
    
16. gdf_B = gpd.GeoDataFrame({"geometry": [polygon_B]}, crs="EPSG:4326")
    
17. AR_B = AR.rio.clip(gdf_B.geometry, gdf_B.crs, drop=True)
    
18. AR_B = mo.AreaWeightMean(AR_B)
    
19. AR_B = mo.Standardization(AR_B)
    
20. AR_B.plot(label=region[re])
    

复制代码

tosaka

{:5_213:}

一大碗年糕

更新一下