问题标签 [sentinel2]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
2214 浏览

python - 将 Sentinel-1 SAR 图像的像素位置转换为地理坐标(纬度、经度)

如何从 Sentinel-1 合成孔径雷达 (SAR) 卫星图像中的 x、y 位置获取地理坐标?

例如,我可以访问下载的图像sg信息

但是,在 Python中使用ESA 的 snap 引擎时,如何获得所需点 (x, y) 的纬度经度sg

0 投票
1 回答
2512 浏览

python - 将 Sentinel-1 SAR 图像的地理坐标(经纬度)转换为像素位置 (x,y)

如何从 Sentinel-1 合成孔径雷达 (SAR) 卫星图像中的地理坐标获取 (x, y) 像素位置?

我可以访问下载的图像信息 sg

但是,如何在 Python 中使用 ESA 的 snap 引擎获得所需纬度和经度的 (x, y) 像素位置?

0 投票
1 回答
1838 浏览

python - 从 ESA Sentinel-2 产品在 python 中加载 RGB 图像并使用 openCV 保存

从 ESA 快照来看,对于 RGB 图像,我们应该将波段 4 放入红色通道,将波段 3 放入绿色通道,将波段 2 放入蓝色通道。我们如何使用 python 将这些波段读取到一个numpy数组中,以便我们可以进行任何我们想要的图像处理,然后将 RGB 图像保存在磁盘上?

例如,这是上述变量之一的类型(其他变量相同):

如何从这些数据中获取 numpy 数组,然后将它们作为 jpg 图像保存到磁盘?

0 投票
1 回答
780 浏览

python-3.x - 合并卫星图像并保留坐标

感谢您光临这里。我目前正在做一个项目,但我对 python 还不是很擅长。所以我希望对这个问题有一些建设性的反馈。

我有一个包含核心样本的数据集,所有样本都存储有样本 id、纬度、经度、内容和与此问题无关的其他数据。

现在我已经导入了这个数据集并按照我想要的方式对其进行了切片。对于图像,我使用 rasterio 模块打开 2 个覆盖该地区的卫星图像。我正在使用 utm 模块在 latlong->UTM->Pixel 值之间来回转换(这似乎也在某些点给我带来了奇怪的坐标)。令人恼火的是,两张 Sentinel-2 图像正好穿过地图的中心。当我在采样的顶部做边界框时,这是一个问题,因为我需要从该区域中提取 10x10 像素的切口。这导致很多样品没有得到适当的切割。

所以我想为什么不将这两个图像合并成一个大的矩形位。但我仍然需要使用 UTM 坐标保留元数据。

你会建议我如何进行。能以简单的方式完成吗?我忽略了另一个角度吗?

感谢您的时间。

在此处输入图像描述 在此处输入图像描述

0 投票
1 回答
915 浏览

java - 使用 JAI 加载 JPEG2000 图像

我想使用 Java JAI 加载 jpeg2000 图像 (.jp2)。

我已经在使用仍然支持 jpeg2000 的旧 JAI 版本。我从 scihub.copernicus.eu 下载了这些图像。这些文件相当大(最大 100 MB)。我尝试了几种不同的方法:

  1. 使用 java 的标准方式是这样的:

    /li>

以下堆栈跟踪发生异常:

我希望得到 JPEP2000 文件的 BufferedImage。

  1. 这是 JAI 加载图片的默认方式。

    /li>

但是变量 image 不包含任何数据: 空对象 方法 getAsBufferedImage() 抛出以下异常:

任何其他图片格式都可以正常工作。您使用哪些 JAI 库来加载 JPEG2000 文件?Javax 或 github 或任何其他提供的库?你的代码是什么样的?

0 投票
0 回答
38 浏览

r - 如何从 R 中的 SAR 光学堆栈中提取反射率值

我在 QGIS 中堆叠了 sentinel-1 和 sentinel-2 图像。我在 SNAP 中进行了预处理,并使用 QGIS 中的“对齐栅格”工具将图像重新采样到 20m 并将它们注册到一个公共投影。我想在 R 中执行分类,但提取反射率/反向散射值需要很长时间。

在将我的 shapefile 和我的光栅读入 R 之后,我陷入了反射率值提取步骤。

我希望最后一行使用 shapefile 提取反射率值以生成数据帧

0 投票
1 回答
3337 浏览

python - Google Earth Engine Python API:使用波段列表在图像集合上映射函数

我一直在使用Sam Murphy的代码对 Google Earth Engine 中的 Sentinel-2 图像进行大气校正。一切顺利,它对于单个图像的运行速度非常快。我想做的是将以下代码映射到图像集合上:

我想我需要一个双映射函数(以避免在此处使用 for 循环),但尚未在 GEE 中看到任何 Python 示例。

到目前为止,这是我想出的:

但是,这不会返回所有波段的图像,所以我觉得我的嵌套函数不能正常工作。我是 python 新手,所以感谢所有帮助!

旁注:我已经尝试过 Sam 的第二个存储库,用于对图像集合进行大气校正,但它运行得太慢了,我更愿意使用建议的地图函数进行“服务器端”计算,因为我有大量图像要过程。

PS:下面是surface_reflectance函数的代码,从 Sam Murphy 的存储库中提取。它调用了他定制的一个名为“Atmospheric”的课程。用于大气校正的模型是 Py6S 模型。

0 投票
1 回答
2415 浏览

python-3.x - 如何在 tfrecord 中使用 python API 从谷歌地球引擎下载哨兵图像

在尝试下载特定位置的哨兵图像时,tif 文件默认在驱动器中生成,但 openCV 或 PIL.Image() 无法读取它。下面是相同的代码。如果我使用文件格式作为 tfrecord. 驱动器中没有下载图像。

我希望输出在 python 中是可读的,所以我可以转换为 numpy。如果文件格式为“tfrecord”,我希望该文件能够下载到我的驱动器中。

0 投票
1 回答
413 浏览

python-3.x - 如何堆叠哨兵乐队?

我使用许多 Sentinel-2 图像,我想将它们的 12 个条带堆叠到一个文件中。我的图像采用 envi 格式(img 和 hdr)。

我尝试通过应用以下代码与模块 rsgislib 进行连接:

但是无论我改变什么参数,我总是会得到同样的错误信息:

你对我有什么建议吗?

0 投票
1 回答
108 浏览

python - 使用 gdalwarp 重新采样会导致 IndentationError: unexpected indent

我使用 Sentinel2 图像,我正在尝试重新采样它们。

我尝试了以下代码:

但我仍然收到相同的错误消息:

我已经尝试过使用相同类型的代码来制作一个子集并且它有效。我不明白我的错误来自哪里。