卫星数据Python中u、v风分量的卷曲计算问题的回答

卫星数据Python中u、v风分量的卷曲计算

回答此问题可获得 20 贡献值，回答如果被采纳可获得 50 分。

<p>我不知道如何获取卫星数据中风的u和v分量的导数。我想我可以用数字梯度以这种方式：</p> <pre><code> from netCDF4 import <a href="https://www.cnpython.com/pypi/dataset" class="inner-link">dataset</a> import numpy as np import matplotlib.pyplot as plt GridSat = Dataset('analysis_20040713_v11l30flk.nc4','r',format='NETCDF4') missing_data = -9999.0 lat = GridSat.variables['lat'] lat = lat[:] lat[np.where(lat==missing_data)] = np.nan lat[np.where(lat > 90.0)] = np.nan lon = GridSat.variables['lon'] lon = lon[:] lon[np.where(lon==missing_data)] = np.nan uwind_data = GridSat.variables['uwnd'] uwind = GridSat.variables['uwnd'][:] uwind_sf = uwind_data.scale_factor uwind_ao = uwind_data.add_offset miss_uwind = uwind_data.missing_value uwind[np.where(uwind==miss_uwind)] = np.nan vwind_data = GridSat.variables['vwnd'] vwind = GridSat.variables['vwnd'][:] vwind_sf = vwind_data.scale_factor vwind_ao = vwind_data.add_offset miss_vwind = vwind_data.missing_value vwind[np.where(vwind==miss_vwind)] = np.nan uwind = uwind[2,:,:] vwind = vwind[2,:,:] dx = 28400.0 # meters calculated from the 0.25 degree spatial gridding dy = 28400.0 # meters calculated from the 0.25 degree spatial gridding dv_dx, dv_dy = np.gradient(vwind, [dx,dy]) du_dx, du_dy = np.gradient(uwind, [dx,dy]) File "<ipython-input-229-c6a5d5b09224>", line 1, in <module> np.gradient(vwind, [dx,dy]) File "/Users/anaconda/lib/python2.7/site-packages/nump/lib/function_base.py", line 1040, in gradient out /= dx[axis] ValueError: operands could not be broadcast together with shapes (628,1440) (2,) (628,1440) </code></pre> <p>老实说，我不知道如何计算（0.25x0.25）度间距卫星数据的中心差。我认为我的dx和dy也不正确。如果有人能在卫星数据中处理这些类型的计算，我将非常感激。谢谢您！！在</p>

0 条评论
分类：Python问答

默认排序时间排序

1 个回答

匿名 1天前

　擅长：python、mysql、java

卫星数据Python中u、v风分量的卷曲计算

1 个回答

相关Python问题