问题标签 [satellite-image]

我正在尝试使用 matlab 从卫星数据创建一个真彩色 RBG 图像，但我不知道该怎么做。

假彩色RGB图像很简单，只要使用正确的红绿蓝通道就可以了

在这种情况下，我使用的是卫星频道的反射率，范围从 0 到 1，所以我不需要修改原始数据

但是当我尝试绘制真彩色图像时遇到了很多麻烦。

根据文献，以下配置文件应该从 MERIS Level-1b 数据产品（我正在使用的数据）中产生良好的 RGB 图像。红色、绿色和蓝色分量的线性组合基于 CIE 1931 颜色空间的颜色匹配函数。

Radiance 是从 0 到 400 的真实值（应用了比例因子），所以我想我必须标准化 RGB 数组（0-1 或 0-255）来创建图像。

但是自己进行标准化或仅使用 im2uint8 不会产生正确的图像。

很可能我做错了一切，因为我不熟悉颜色配置文件。matlab 中有没有一种方法可以直接使用 CIE rgb 组合（我认为我从上述公式中得到的那个）来创建图像？

有没有人熟悉使用 matlab 和卫星数据的图像？

谢谢！

我目前正在尝试使用 curl 从 esa 的 Copernicus / Sentinel 项目下载卫星图像。不幸的是，下载一直停止在 90% 左右，并且 php 脚本返回内部服务器错误 (500)。

因此，我想以特定的字节数恢复下载。似乎 esa 服务器只是忽略了 http-range-header (CURLOPT_RANGE) 并且 CURLOPT_RESUME_FROM 也没有改变任何东西。

如果我使用 Google Chrome 手动下载文件，下载也会中断，但会在一段时间后继续。

因此，如果谷歌浏览器可以恢复下载，curl 也应该能够做到这一点。对于如何做到这一点，我将不胜感激。

一些细节：

我要下载的文件在这里(420MB)，要访问它，您需要在 scihub.esa.int/dhus/ 注册。
内容类型是application/octet-stream

我的代码：

它非常适用于较小的文件（我已经用一些图像和 pdf 文件对其进行了测试），我还可以下载大部分卫星图像（下载了前 380MB）。我也尝试增加超时值，但脚本在 5 分钟之前就终止了。

我试过了curl_setopt($open_file, CURLOPT_RESUME_FROM, 1048576);，curl_setopt($open_file, CURLOPT_RANGE, "1048576-");但文件总是以相同的字节开头。

编辑：

我无法回答我的问题，但对于这种特定情况，我找到了一种解决方法。所以，如果有人读到这篇文章并且偶然想用 cURL 下载这些卫星图像，我就是这样做

的：和curl_setopt($open_file, CURLOPT_RESUME_FROM, $bytes_already_loaded);可以跳过先前加载的字节并继续下载（这对于图像文件是不可能的）。因此，请使用此链接而不是图像文件。

我正在尝试使用BEAM java API 导出 ALOS AVNIR-2 产品的 PNG quicklook。下图显示了产品的 RGB 预览，它出现在 beam 的 GUI 中。

如您所见，由于其地理编码，图像不是直立的。我开发了一个非常简单的 java 程序来导出产品的 quicklook。

该程序有效，但我从中获得的每个 PNG 图像都是一个直立的 PNG 图像，如下所示。

现在，我知道不可能在 PNG 图像中包含地理编码信息。我只需要重现图像的相同“旋转”。任何的想法？

我正在使用 Matlab 上的卫星图像，并正在编写用于计算净辐射的代码。输出包含带有NaN值的矩阵形式的图像。当我尝试使用以下方法绘制它时：

它将NaN值分配为与 0 相同的颜色。有什么方法可以将这些值分配为不同的颜色？

另外，我想将我的矩阵保存为图像，但其值不会改变的格式......可能是ASCII。是否有将矩阵转换为 ASCII 的工具？

我有来自 TLE 的地球观测卫星的星历，所以我可以获得卫星的位置（纬度/经度）。该文档为 ISS 提供了一个很好的示例，展示了它是如何完成的。

然而，我真正追求的是卫星传感器指向的地球表面点的位置。假设传感器指向垂直于飞行方向，所以我应该能够使用卫星的倾角计算方位角方向。然后我知道卫星传感器用来指向地面的角度。一些简单的几何应该给我我需要的高度。

现在的问题是：我如何将这两件事串在一起？计算（观察者）计算示例返回位置的高度和方位角（纬度/经度）。但是，我找不到任何相反的方法。有什么想法吗？任何线索都非常感谢。

我需要创建一个函数，以在时间跨度为 2（xt+1 和 xt-1）的平均值上对 landsat 图像时间序列上的噪声值进行插值。

我正在使用 fmask 产品来检测云和阴影，然后应用插值。

对于一个时间序列：

由于 c2 是 fmask 时间序列的向量（2 代表云，4 代表阴影），而 t2 是 evi 时间序列的向量：

但是使用 raster 包的 calc 函数是不可能的，因为它不适用于具有 2 个参数的函数。

关于如何处理此问题并将此插值应用于栅格时间序列的所有像素的任何建议？

我正在尝试这个，但它仍然无法正常工作：

我发现已经回答了一个类似的问题，但它对我不起作用。我需要调用 Google Satellite 作为基础层。可能我没有在正确的位置输入基础层代码以使其正常工作。这是我的页面的链接，如果我能获得有关在何处添加基础层代码以调用 Google Satellite 作为基础层而不是 OSM 的帮助，我将不胜感激。谢谢！

我在 GTiff 中有一张卫星图像，范围为 0 - 65535，我需要将其转换为字节格式（0-255）。

使用：gdal_translate [-a_nodata 0 和 -scale 0 65535 0 255] -ot Byte infile.if outfile.tif

这工作正常，但我得到很多像素被四舍五入（截断）为 0，这是我的 nodata 值，这意味着它们在可视化时变得透明。

我尝试过使用 -a_nodata 0 和 -scale 1 65535 0 255，但我找不到适合我的解决方案。

我正在寻找的是 0 作为 nodata 和 1-255 作为数据范围。

我正在编写代码以沿 numpy 数组的长度添加数据（用于组合卫星数据记录）。为了做到这一点，我的代码读取两个数组，然后使用函数

使一行中的两个组成一个数组

如果考虑中的数组宽度相等，则 swath_stack() 函数将简化为 np.vstack()。我的问题是我MemoryError在这个阶段不断遇到。我知道在算术运算符的情况下，在适当的位置进行算术（即array1 += array2，与array1 = array1 + array2

有人可以帮忙吗？

为了 100% 可重现，整个代码（和使用的图像）如下。

以下是代码中的步骤。

1. 我将轮廓 ( [[386, 330], [398, 320], [382, 300], [370, 310], [386, 330]]) 叠加在图像上。这很好用。

2. 使用 theta (= 90) 旋转图像warpAffine

3. 在完全相同大小和相同旋转功能的虚拟图像中旋转轮廓中的每个点。

4. 找到旋转点的坐标。

5. 在旋转图像上叠加旋转点。

我希望他们排队，但他们没有。

可能是什么原因？

这是具有对齐轮廓的原始图像

这是旋转后的图像和叠加。