问题标签 [satellite-image]

我正在使用 CNN 研究卫星分割问题。但在达到 50 纪元时，骰子系数达到 0.995。网络代码如下。

所以我想知道它是网络训练的问题还是与训练/验证数据集有关？模型有问题吗？

使用时，rasterio我可以执行以下任一方法来获取栅格的单个波段：

它会打印出来：

但我很好奇哪个“更好”？我假设索引到 numpy 数组比从文件中读取要快得多，但是打开其中一些大型卫星图像会占用大量内存。有什么好的理由去做其中一个吗？

以前，我能够在卫星地图上叠加疾病地图。现在看来需要一个 API 密钥。我无法将其集成到 R 代码中。

到目前为止，我已经尝试访问谷歌地图上的代码。我得到了一个代码，但我无法将它集成到我的 R 代码中。

我认为我需要在此处包含一些内容，其中包括相关地图的 api 密钥

我正在尝试使用 OpenLayers 5 在我的地图上添加卫星图像。

问题是我无法做到这一点，因为我刚刚找到了一个选项，可以在地图上添加图像，传递图像范围（xmin、ymin、xmax、ymax）而不是边界框。图像应适合边界框。由于这个原因，图像被扭曲了。

图像为 JPG 文件（属性 feature.properties.icon）。示例：http ://exampleserver.com/220/063/353LGN00/353LGN00_thumb_large.jpg

我想要的结果是这样的：

我得到的结果是：

我在地图上添加此图像的代码如下：

有人知道如何传递图像边界框而不仅仅是图像范围吗？

先感谢您。

编辑 1：迈克的解决方案

通过 Mike 的解决方案，我能够修复一些图像存在的错误（靠近赤道线）。出于这个原因，他的回答解决了我的问题，并将图像插入到我在创建问题时所期望的更好位置。

然而，这个解决方案对我来说适用于赤道线附近的图像。两极旁边的图像保持扭曲（编辑 2）。

我在下面发送了一张说明最终结果的图片：

编辑2：新问题？

我正在测试一些图像，我发现了一个新错误。现在我发现图像应该适合边界框。如果图像不适合 bbox，它会保持失真，例如我在下面发送的打印说明。

该图像应适合 bbox 内，如下图[PS 1]所示：

我相信这可能是重投影的问题，但我不知道，因为视图投影和图像投影都是EPSG：4326。

我试图按照 Openlayers 网站上关于 Raster Reprojection[1.] 的解释进行操作，但是我无法重现它，因为正如我所说，两个投影（视图和图像）是相同的（或者它们应该是相同的）。

我在包含与上图相关信息的 GeoJSON 下方发送。该图像可以在“properties.icon”（http://www.dpi.inpe.br/newcatalog/tmp/MOD13Q1/2018/MOD13Q1.A2018017.h13v14.jpg）中找到。bbox 坐标可以在“geometry.coordinates”或“properties.bl_latitude”、“properties.bl_longitude”、“properties.br_latitude”等中找到。“bl”表示“左下”，“br”表示“右下”，“tl”表示“左上”，“tr”表示“右上”。“properties”中的这些坐标在“geometry.coordinates”中是相同的。

有人会有新想法吗？

[PS 1]：使图像适合 bbox 的原始代码是传单代码 [2.]，我将其发送到下面：

[1.] https://openlayers.org/en/latest/doc/tutorials/raster-reprojection.html

[2.] https://github.com/ScanEx/Leaflet.imageTransform

我试图总结 NDVI 值的卫星图像栅格的值。我已经提取了数据，它似乎已经填充了一个我可以在控制台窗口中看到的列表，但是我无法使用“查看”功能看到它，也无法将其转换为 .csv。

示例数据

提取值

写入文件

我希望能够在某个时候使用 NDVI 图像的单元格值的分布将其与另一组数据相关联，所以为了做到这一点，我需要有一个所有值的工作 csv细胞。先感谢您。

通过非 matlab/非本地工具 (GDAL) 以及本地工具 (geoimread) 的组合，我可以将 Sentinel-2A 数据提取为单个波段或作为采用 gdal 合并的 RGB 图像。我被困在使用

产生黑色图像，显然没有信号。图像中的强度值范围是 271 - 4349。我知道图像中有一个很好的信号，因为当我这样做时：

即索引图像，收集颜色图，设置颜色图然后调用image函数，我得到了我正在探索的区域的相似之处（尽管颜色很奇怪）

我目前正在考虑是否应该尝试：

1. 查找 Sentinel-2A 数据的低级描述，实现缩放/校正

2. 使用工具箱，可能是这个。

3. 可能在涉及 GDAL 的早期步骤之一中调整输出设置 非常感谢您的评论或建议。

一个基本的缩放方案是：

uint16注意到从doublefor转换为的重要性imshow。

我以前见过这个问题，但似乎专注于导出到 GDrive，我想将图像导出/下载到我的计算机而不是谷歌驱动器。

我正在尝试使用地球引擎的 python api 从图像集合中下载每个图像。

按照他们的示例，抓取单个图像并获取下载 URL 很容易：

这给了我一个 URL，我可以从中下载图像。

同样对于图像集合，我可以使用以下方法获取图像集合：

是否可以遍历每个图像并获取每个下载 URL？即想在不使用谷歌驱动器的情况下每月抓取一次相同区域的图像或类似的东西。现在我把问题简单化了，但会在下载之前将图像剪辑到我的区域并根据需要执行其他操作。

我在这里找到了这个问题：How to iterate over and download each image in an image collection from the Google Earth Engine python api

但无法提出最终解决方案...

帮助将不胜感激！

I have an array that is 1648 x 2608 (from a GOES-16 image, part of a full disk image).

When using cartopy, I can't create a 1648 x 2608 PNG, with the size of my array.

Instead, it creates a 1648 x 2608 white box with the satellite image inside (smaller than 1648 x 2608).

Without using cartopy, I can easily get a full resolution 1648 x 2608 pixels PNG with the code below:

This is the result I get: Full Resolution (1648 x 2608) PNG Without Cartopy

However, when using cartopy, I do not get a full resolution image as before:

This is the result I get with cartopy: Image with 1648 x 2608 pixels, however, white bounding box

How can I save a PNG with cartopy with the size of my array, in order to get a full resolution plot?

Note: I could create it with Basemap without any problems, however, I'm trying to migrate to Cartopy.

Thanks!

我正在尝试使用 Sentinal 下载某个地区的地图。我正在使用“BANDS-S2-L1C”进行身份验证。当我尝试运行算法来获取该区域的数据时，它会抛出“服务器响应：“未找到层 BANDS-S2-L1C”

“谁能告诉我如何将我的哨兵 ID 与程序链接。如何将图层链接到 BANDS-S2-L1C” layer = 'BANDS-S2-L1C'

我得到的错误是

DownloadFailedException: 在执行任务 S2L1CWCSInput: 下载失败: https://services.sentinel-hub.com/ogc/wcs/ecc7a293-1882-4cff-8bf8-918a2e74baff?SERVICE=wcs&MAXCC=30.0&ShowLogo=False&Transparent=True&BBOX= 44.97%2C27.67%2C45.26%2C28.03&FORMAT=image%2Ftiff%3Bdepth%3D32f&CRS=EPSG%3A4326&TIME=2019-04-01T07%3A18%3A16%2F2019-04-01T11%3A18%3A16&RESX&RESX=20m&RESY=2 BANDS-S2-L1C&REQUEST=GetCoverage&VERSION=1.1.2 with HTTPError: 400 Client Error: Bad Request for url:https://services.sentinel-hub.com/ogc/wcs/ecc7a293-1882-4cff-8bf8-918a2e74baff?SERVICE=wcs&MAXCC=30.0&ShowLogo=False&Transparent=True&BBOX=44.97%2C27.67%2C45.26%2C28.03&FORMAT =image%2Ftiff%3Bdepth%3D32f&CRS=EPSG%3A4326&TIME=2019-04-01T07%3A18%3A16%2F2019-04-01T11%3A18%3A16&RESX=20m&RESY=20m&COVERAGE=BANDS-S2-L1C&REQUEST=2 响应。 ：“未找到层 BANDS-S2-L1C”

给定一个栅格时间序列（相同空间范围内的卫星图像，每个图像都在不同的日期拍摄），如何将日期添加为（地理）Tiff 中每个图层/波段的名称？它可能未指定吗？

以 R 为例：

gdalinfo输出如下（无波段名称信息）：

有一个Tiff 标签库，但我找不到任何（语义上合适的）标签可以满足我的需要。

是否有一个好的做法如何将图层/波段名称添加到（Geo）Tiff - 最好可以由 QGis 等其他软件解释/读取。