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

一大碗年糕

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

一大碗年糕

用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:}