问题标签 [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 投票
0 回答
23 浏览

python - 为什么 sentinel-2 波段包含 25000 左右的值?

我正在尝试使用 Python 处理 .jp2 文件。数据来自 sentinel-2,因为它写入 MTD_MSIL1C.xml 文件。

当我使用 rasterio 库读取文件时,我收到值大于 20k 的 numpy 数组。

img = rasterio.open('some_band.jp2', driver='JP2OpenJPEG').read(1)

对于以下命令 img.max() ,我得到了这样的值23526

GDAL 版本 == 3.0.2 光栅版本 == 1.2.6

您每 90 分钟只能发布一次。

我想可视化 .jp2 文件。为此,我需要 0 到 255 范围内的值。简单的归一化和乘以 255 会导致图像变暗。

我试图将图像除以 3000 并将其传递给 matplotlib。生成的图像看起来不错,但为什么会这样?

我也尝试了这个命令,它产生的结果范围为 0 和 255 gdal_translate -scale 0 3000 0 255 -ot Byte

0 投票
0 回答
72 浏览

python - GEE 为图像集合中的图像添加波段,其中包含图像集合中其他图像的最接近信息

我有两个图像集。第一个带有 Sentinel 1 图像,第二个带有 ndvi 图像。

我想屏蔽所有 ndvi > 0.6 的 Sentinel 1 图像像素,以计算低到中等植被区域的平均值。

我面临的问题是两个集合中的图像不共享相同的日期。我写了一些代码,我得到了值,但我不确定我到底在做什么。

我想也许有一种方法可以将新乐队(ndvi)添加到 s1 集合中,然后屏蔽和减少,但我不知道如何在不同的日期完成这项工作?

0 投票
1 回答
20 浏览

python - Python 请求 - 由于 URI 和身份验证过程的性质,下载图像并写入文件不起作用

我正在编写一个使用sentinelsatPython API 下载 Sentinel 2 产品(卫星图像)的脚本。

产品的描述采用 JSON 结构并包含参数quicklook_url

例子:

https://apihub.copernicus.eu/apihub/odata/v1/Products('862619d6-9b82-4fe0-b2bf-4e1c78296990')/Products('Quicklook')/$value

任何 Sentinel API 调用都需要凭据。检索产品并打开存储在其中的链接也是如此quicklook_url。当我在浏览器中调用示例时,系统会要求我输入用户名和密码以获取

在此处输入图像描述

同名S2A_MSIL2A_20210625T065621_N0300_R063_T39NTJ_20210625T093748-ql.jpg

不用说我只是从 API 开始,所以我可能遗漏了一些东西,但是

产生 0KB 损坏的文件和

产生 1KB 损坏的文件。

我调查过requests.Session

问题是我找不到发布会话身份验证所需的 URL。我有点确定sentinelsatAPI 可以做到这一点,但我的外观并没有产生任何成功的结果。

我目前正在SentinelAPI上课。它具有download_quicklook()我现在正在使用的功能,但我仍然很好奇如何在没有该功能的情况下做到这一点。

0 投票
0 回答
25 浏览

javascript - 更快地从 Google 地球引擎导出结果

我有一个大约 68,000 个多边形的 shapefile,对于两个日期之间的所有 Sentinel 2 图像(例如,从 2021 年 7 月 1 日到 2021 年 8 月 1 日),对于所有 Sentinel 2 波段,我想计算所有 Sentinel 2 波段的平均值每个多边形中包含的像素。输出应该是每个波段的表格,其中每一行对应于每个多边形,每一列对应于每个 Sentinel 2 图像(+ 一些属性对应于 shapefile 中包含的多边形)。此输出将导出为 CSV 文件以在本地计算机中使用。

我试图只为一个波段 (B3) 执行此操作:https ://code.earthengine.google.com/476b39b8a93872506bb9bb17dd8705d3问题是完成这些约 68,000 个多边形的过程需要很长时间,而且只有一个月的数据. 这个想法是将多边形的数量增加到约 2,000,000 个,并将日期范围增加到 3-5 年。

有什么方法可以加快这个过程吗?

我想在服务器端计算每个多边形的每个波段的平均值的过程很快;但是,问题必须是从服务器端获取结果并将它们写入导出到驱动器的文件中。

0 投票
2 回答
33 浏览

python - 如何使用 python 从 XML 文件中检索特定信息?

我正在使用 Sentinel-2 图像,我想从 XML 文件中检索 Cloud_Coverage_Assessment。我需要用 Python 来做这件事。

有谁知道如何做到这一点?我想我必须使用 xml.etree.ElementTree 但我不确定如何使用?

XML 文件:

0 投票
0 回答
16 浏览

bash - gdal_translate 仅将前三个波段从 .vrt 翻译为 .tif

我正在尝试按照本教程将一组 .tif 图像转换为一个多波段 .tif 。我没有收到任何错误,但是每次我运行脚本时,当我有比这更多的输入 .tifs 时,它都会生成一个三波段图像。我已经尝试了 -r 标志,它没有改变任何东西。

我有:

"$outdir/btif_${roi}.txt"像这样的 GeoTIFF 文件路径的纺织品在哪里:

我正在处理 Sentinel-2 图像并使用 OSX 11.6。

0 投票
0 回答
12 浏览

tiff - 如何确定特定土地利用类型的两个图像的差异

我正在处理 2 个图像,image-1 是 xarray DataArray,image-2 是光栅 .tif 数据。我想覆盖这 2 个数据以查看落在 xarray (image-1) 中特定值内的土地利用类型 (image-2)。下面是我的代码:

现在我的问题是如何确定与两个图像之间的特定土地利用类型相对应的图像差异。