我正在尝试使用树莓派从卫星天气图像中提取数据，然后愿意通过单个树莓派运行数据挖掘。那可能吗？那很费时间吗？

我一直在使用 SMAP 数据卫星，专门用于湿度和土壤比例。

我遵循使用 GDAL 解决所有问题的想法，并制作了类似于Link to first approach to download SMAP data

修改代码和测试：

https://n5eil01u.ecs.nsidc.org/SMAP/SPL4SMAU.003/2017.08.01/SMAP_L4_SM_aup_20170801T030000_Vv3030_001.h5

因此，使用这种方法，结果是具有许多投影问题的全局地图，例如下面的图像，由上面的 python 代码生成。

为了与正确的数据进行比较，使用 HEG nasa 软件从 h5 中提取了相同的图像。

我目前正在研究一个使用 MODIS 数据集的机器学习项目。由于我的 PC 不符合项目的计算要求，所以我使用了 AWS 服务器。现在的问题是地球引擎正在将图像导出到 Google Drive 或 Google Cloud Storage，但我希望它们导出到我的 S3 存储桶。

我遇到了建议将数据下载到本地存储然后将它们上传到 S3 存储桶的答案。鉴于庞大的数据集和我糟糕的数据速度，我需要很长时间才能做到这一点。因此，我想直接使用 Earth Engine 将它们导出到我的 S3 存储桶。

我已经阅读了发生导出的文档（ee.batch.Export.image）。我正在考虑编写一个将 Geotiff 图像导出到 AWS S3 存储桶而不是 Google Drive 或 Cloud Storage 的函数。

附言

• 我已经使用 Amazon MODIS 公共数据集进行了验证，并且我想要的数据集（MOD09A1 和其他一些）不是由 Amazon 提供的。
• 我的电脑上安装了 Windows 10。

我是 Google 地球引擎的新手，正在尝试获取整个刚果民主共和国的 Landsat 图像。该项目涉及一些计算机视觉和图像分割，因此我需要尽可能获得最高分辨率。

我有一些从地球引擎开发的代码，用于从下面发布的 1 年陆地卫星图像合成。

问题是每次我运行脚本时都会询问我的scale值。所以我要求最高分辨率，但查询一直出错，因为它说我已经超过了图像导出的最大像素限制。

所以基本上我需要弄清楚如何将刚果分割成一组块，地球引擎允许我为这些块提取最大分辨率的合成图像。有谁知道我如何计算正确大小的多边形以适应账单？

我正在尝试做的事情：

我想要做的是有一个最近邻居的功能，用户提供纬度和经度数组和一个奇点。在函数中，我不介意应用来自不同库（例如 scipy kdtree）的函数。以前，我应用了以下函数，该函数可用于查找 1D lat 和 lon 数组的最近邻：

我知道有一种更有效的方法可以找到最近点的索引，但这种方法对于我之前正在做的应用程序来说很快。

更多细节：

我现在遇到的问题是我的函数只有在 lat 和 lon 是 1D 时才有效。使用我的新数据集，我的 lat 和 lon 数组都是 2D 的，每个都具有 (5500, 5500) 的形状，因此分解它并应用我的函数将非常低效，因为我在 3 个月内每 5 分钟执行一次此分析。这种方法的一个例子是：

我发现这个链接解释了如何使用纬度和经度坐标执行 kd 树，但纬度和经度也是一维的。

我想返回索引，因为我在后面的步骤中将它应用于我的 lat、lon 和 data 数组。

所以我的一般问题是任何人都可以解释/提供一个示例（或链接），说明如何为指定的纬度和经度坐标（例如 36.538537 N，126.330194 E）提取特定数量的邻居（在本例中为 4），前提是纬度和经度数组是二维的（即我需要从中提取邻居的数据集）。

如果有用，我正在处理卫星可见数据。

我已经尝试过：

以前我曾尝试解开数组（如上所示），这就是我发现它效率极低的原因。我也让它返回 True/False 值，然后查找我的数组在哪里是 True，但这比第一种方法效率更低。我尝试使用 heapq 库中的 nsmallest 没有运气，因为我不明白如何正确使用该库。

我进行了广泛的搜索，试图为我的问题找到解决方案，但没有成功，所以我没有选择，现在正在询问社区。

我正在尝试绘制海面温度数据并添加陆地的彩色图像，以免数据与NAs. 我尝试了多种方法来做到这一点，但正如您将在下图中看到的那样，地图相对于数据没有正确排列。

为了使这个问题可重现，这里是一个包含我正在使用的文件的保管箱的链接：https ://www.dropbox.com/s/e8pwgmnhvw4s0nf/sst.nc4?dl=0

这是我迄今为止开发的代码；

通过 ncdf4、rasterToPoints、map_data 和 ggplot2

通过 Raster、maptools 数据、SP 变换和基础绘图

- 我用于这个空间数据的投影是EPSG:4326

- 这是指示sst.nc4输出投影的 XML 片段

我也尝试使用map()带有mapproj'projection参数的函数，但它似乎没有一个伪墨卡托投影作为选项。

我想获取栅格（卫星图像）数据，并构建一个Datasetor DataArray，以加快我的图像处理速度（我必须经常处理多波段、多日期的卫星图像）。

数据作为每个图像日期的单独波段出现，我了解如何将每个波段日期转换为 xarray- DataArray。我认为为每个波段设置一个变量是最有意义的，并且在每个波段内都有空间（x，y）和时间维度。

但是，我不知道该怎么做。

我一直在与一些虚拟乐队合作以试图解决这个问题，因此将包括在内以澄清我的数据是什么样子以及我正在尝试做什么。

然后我想将它们组合成一个DataArray，带有两个变量（Band4 和 Band5），每个变量都包含两个图像日期......但不知道如何进行。

创建/导入数组时是否需要添加更多坐标或维度，然后concat沿着这些维度添加？

我正在使用随机森林分类器对 Sentinel 2 图像进行监督分类。我希望从图像中选择训练数据。谁能告诉我有效执行此操作的方法？

我正在尝试读取具有多个波段的图像（多光谱图像）。正常图像通常具有 3 个波段 (RGB)。我了解到某些类型的图像（例如：卫星图像）可以有 3 个以上的波段，如红外线等。

我知道如何Bitmap从类中读取普通图像，System.Drawing并以红色、绿色和蓝色值获取每个像素的颜色值。

使用 C# 如何读取其他频段？我发现很少有图书馆，但无法理解它们是如何工作的（甚至不确定它们是否相关）。LibTiff.net 就是其中之一（这听起来很有希望，因为多光谱图像扩展名是 .tiff）虽然我在图像中找不到任何关于“波段”的参考。

我正在使用谷歌地球引擎，我正在尝试使用另一个图像（波段）提取/过滤（剪辑）波段中的像素。我计算了 NDVI 并创建了一个阈值，以渲染 NDVI > 0.3 的图像，但我想提取可见 NIR 波段中的相应像素。
这是片段代码。

我现在想使用阈值变量（图像）在“B”、“G”、“R”和“NIR”波段中剪辑相应的像素。显然，我被困在这里。请让我知道是否有一种方法可以使用 GEE 中的另一个波段过滤/剪辑一个波段的像素。该任务类似于在 QGIS 中使用 Clipper，如果这不起作用，这是我剩下的选项。

谢谢你的帮助！