问题标签 [pyresample]

我编写了一个 python 脚本来投影和叠加来自邓迪大学的地球静止卫星图像，因此生成的图像可用于 xplanet 渲染地球表面。该工具的源代码可以在https://github.com/jmozmoz/cloudmap/tree/cartopy找到（这是支持 cartopy 的分支）

该工具支持两个不同的 python 库来在平面地图上投影地球静止图像：pyresample 和 cartopy。

我发现了以下差异/问题：

1. pyresample 比 cartopy 快得多（取决于输出图像的大小高达 10 倍）
2. 输出图像不同：使用 pyresample 的结果显示出更强的对比度。有关示例，请参见https://github.com/jmozmoz/cloudmap/tree/cartopy/debug的调试目录

3. 如果使用多处理库并行进行投影，cartopy 版本会崩溃并显示以下错误消息：

   /li>

那么为什么 cartopy 这么慢呢？pyresample 是否在 C 代码中完成工作？cartopy 应该支持多处理吗？以及如何解决对比度问题？

谢谢您的帮助

我正在尝试使用 pyresample 模块绘制 ASCAT 海洋风矢量及其 wvc 质量标志的完整条带轨道。可以在此处找到指向 ASCAT netcdf 文件的链接：ftp: //podaac-ftp.jpl.nasa.gov/allData/ascat/preview/L2/metop_a/12km/ 从我读过的关于模块的文档中，确实如此没有很好地描述我们如何在文件中找到满足几何区域定义的信息。下面是一个示例代码，用于绘制大量卫星数据

我从 netcdf 文件中提取了 lons、lats 和 wvc_quality_flag

首先，我收到一个错误，即 pyresample 模块没有该属性图，而文档说明它有，并且“area_def”从未在 ASCAT netcdf 文件中的任何位置定义。pyresample 对于这些类型的文件是不合理的，还是我没有在 ASCAT 文件的元数据中寻找这些定义的正确位置？关于这个模块的一些澄清真的很有帮助！再次感谢！

我需要将不规则网格（朗伯圆锥形）上的数据重新网格化为规则网格。我认为 pyresample 是我最好的选择。事实上，我原来的纬度，经度不是一维的（似乎需要使用 basemap.interp 或 scipy.interpolate.griddata）。

我发现这个 SO 的答案很有帮助。但是我得到空的插值数据。我认为这与我的影响半径的选择以及我的数据被包装的事实有关（？？）。

这是我的代码：

lat_n 和 lon_n 是不规则的，经纬度对应投影坐标 x,y。

由于 lon_n 的方式，我补充说：

所以现在如果我绘制它们，它们看起来不错并且还可以：

然后我创建了一组新的常规坐标：

按照上面的答案，我编写了以下代码：

从 pyresample.geometry 导入 SwathDefinition 从 pyresample.kd_tree 导入 resample_nearest

数据的分辨率大约是30km，所以我放了70km，我放的fill_value是数据中的那个，当然我也可以放0或者nan。

但是我得到一个空数组。

我做错了什么？另外-如果有另一种方法，我有兴趣知道它。Pyresample 文档有点薄，我需要更多帮助。

我确实发现这个答案建议使用另一个 griddata 函数：

似乎没问题：

但我想更多地了解 pyresample，因为它看起来很强大。

我正在按照Pyresample 示例绘制 Sentinel-3 场景和海岸线。但是，似乎pyresample.kd_tree.resample_nearest并pyresample.save_quicklook期望将 2D 数组作为数据，而不是我试图绘制的 3D（3 波段）Sentinel 场景：

伪代码和错误是：

如何绘制这些数据？