从全球视角重新采样和重新投影气象卫星图像到区域范围

Question

在笔记本上阅读并绘制卫星图像后，我发现很难重新采样和重新投影（投影 = PlateCarre）到区域范围（范围 = -1、19、-1、17）。如果海岸线绘制在图像上，我将不胜感激。如果也可以绘制名为“DNKN”（长= 8.53，纬度= 12.05）的电台，我将不胜感激。读取和绘制卫星图像的代码如下。请在认为合适的情况下进行编辑。谢谢。

import matplotlib.pyplot as plt
import cartopy.crs as ccrs
ax = plt.axes()
sat_image1 = '2014_06_16_16Z.jpg'
df = plt.imread(sat_image1)
ax.set_axis_off()
plt.imshow(df)
plt.show()

卫星图像链接 ('2014_06_16_16Z.jpg') 是 https://drive.google.com/file/d/1fE4V9aHIklGdv7NizQommdWzAvR6DpUN/view?usp=sharing

Answer 1

你在这里有几个问题，几乎没有研究证据，在我看来，这对于 stackoverflow 问题来说还不够具体。我强烈建议您阅读 stackoverflow 杰出人物 Jon Skeet 的 article on asking a good question。

要解决的具体问题：

1. 您没有使用 cartopy 创建地理配准地图
2. 您尚未尝试定义图像的原点坐标参考系
3. 您尚未尝试转换为所需的坐标参考系。
4. 您尚未引用专为此目的设计的众多工具中的任何一个

将地理参考图像添加到地图是一项非常重要的任务。我认为你需要做更多的阅读和练习才能解决你的问题。那里有一些很棒的资源。您应该首先确保您了解坐标参考系是什么。这将帮助您理解为什么导入您包含的图像毫无意义。

如果您搜索 cartopy add satellite image，您会很快找到向您展示如何操作的文章：

cartopy 文档有一个带有 plotting an image on a map 的 source code 页面

内容广泛的blog post，其中包含在 cartopy 中绘制整个地球图像的完整指南

Another blog 展示了如何使用完整的 Python 代码在 Cartopy 地图上绘制来自 WMS 的卫星图像。

请花一些时间了解 cartopy 库和 GIS 的性质。从食谱中复制一些示例并根据您的需要调整它们，阅读 GIS 博客文章，处理一些实践数据集。当你有一个明确定义的特定问题时，回来提出一个新问题。