问题标签 [satellite-image]

你知道我可以免费获得地点卫星图像的任何图像吗？

我有一个地形编辑器，它使用 GIS 数据来可视化高度图，但我想用卫星图像对它们进行纹理化。

我希望你们中的一些人可以提供一些来源？

我必须开始应用分析卫星图像来识别一些人造结构。我想为此使用 C 或 Java。

对于卫星，我计划使用谷歌地图数据。

我在这里有三个问题：

1. 除了谷歌地图/地球之外，什么是 GIS 数据的最佳来源。
2. 考虑到我将不得不使用第三方 API，编写此类应用程序的最佳语言
3. 是否有可识别人造结构的开放式图像处理引擎？

那是很多问题，但我希望这里的聪明人可以在这里帮助我。

当卫星图像在某个缩放级别不可用时，是否可以使用谷歌地图 api（javascript，v2）以编程方式检测？我正在自动创建一些地图并将它们设置为某个缩放级别，但在某些情况下，该级别没有可用的卫星图像，在这种情况下，我想自动备份一个缩放级别。api 是否公开了某种确定这种缺乏图像的方法？

到目前为止，我唯一能想到的就是尝试使用 javascript 找出 api 发出的图像请求是否失败，然后根据这些失败的图像请求做出反应，但我真的无法做到工作......这似乎有点不雅。谢谢你的帮助！

背景

我使用来自合成孔径雷达卫星的非常大的数据集。这些可以被认为是一侧大约 10k 像素的高动态范围灰度图像。

最近，我一直在开发Lindeberg 的尺度空间脊检测算法的单尺度变体的应用程序，用于检测 SAR 图像中的线性特征。这是对使用方向滤波器或使用霍夫变换的改进，这两种方法以前都使用过，因为它的计算成本比任何一种都低。（我将在 4 月的JURSE 2011上展示一些最近的结果，如果有帮助，我可以上传预印本）。

我目前使用的代码会生成一个记录数组，每个像素一个，每个记录都描述了像素右下角矩形中的一个脊段，并以相邻像素为界。

一个条目ridges包含一个山脊段，如果恰好两个,top和具有在 0 - 128 范围内的值。假设我有：leftrightbottom

然后我可以找到山脊段的起点 (x1,y1) 和终点 (x2,y2)：

当这些单独的山脊段被渲染时，我得到一个像这样的图像（一个大得多的图像的一个非常小的角落）：

这些长曲线中的每一条都是由一系列微小的山脊段渲染而成的。

确定两个包含山脊段的相邻位置是否相连是很简单的。如果我有ridge1at (x, y) 和at (x+1, y)，那么如果 0 <= <= 128和= ridge2，它们是同一行的一部分。ridge1.right ridge2.leftridge1.right

问题

理想情况下，我想将所有山脊段拼接成线，这样我就可以遍历图像中找到的每条线以应用进一步的计算。不幸的是，我发现很难找到一种既低复杂度又节省内存并且适合多处理的算法（在处理非常大的图像时所有重要的考虑因素！）

我考虑过的一种方法是扫描图像，直到找到一个只有一个链接的山脊段的山脊，然后遍历结果线，将线中的任何山脊标记为已访问。但是，这不适合多处理，因为如果没有昂贵的锁定，就无法判断是否没有另一个线程从另一个方向（例如）走同一条线。

读者建议作为一种可能的方法是什么？过去似乎有人会想出一种有效的方法来做这种事情……

我已经下载了 USDA 作物数据层（数据可在此处获得，请注意大文件： http: //www.nass.usda.gov/research/Cropland/SARS1a.htm）并希望将其导入 R 进行分析。我导入了其他大型 GTIFF 文件（例如 LANDSAT 数据），但在调用该函数时出现以下错误：

其中“文件”是美国农业部提供的任何 .tif 文件（我使用的是中西部地区汇编）。

我怀疑问题是 GTiff 是分类数据（土地利用类别）而不是真正的序数/光谱数据，但我不明白为什么它会在图像编辑器（例如预览）中呈现良好并在 R / GDAL 中退缩.

谢谢您的帮助。

编辑：如果这应该在 GIS 论坛上进行，请告诉我。一般来说，空间问题的编程方法似乎就在这里，不是吗？

我有一个以波段交错像素 (BIP) 格式存储的多波段卫星图像以及一个单独的头文件。头文件提供了图像中的行数和列数，以及波段数（可以超过标准3）等详细信息。

图像本身是这样存储的（假设一个 5 波段图像）：

[B1][B2][B3][B4][B5][B1][B2][B3][B4][B5] ... 等等（基本上 5 个字节 - 每个波段一个 - 每个像素开始从图像的左上角）。

我需要在 Python 3.2（在 Windows 7 64 位上）中将这些波段中的每一个分离为 PIL 图像，目前我认为我正在错误地解决问题。我目前的代码如下：

这段代码打开 BIP 文件的时间太长了，肯定有更好的方法来做到这一点。我也有 numpy 和 scipy 库，但我不确定如何使用它们，或者它们是否会以任何方式提供帮助。

由于图像中的波段数量也是可变的，我发现很难找到一种快速读取文件并将图像分成其组成波段的方法。

只是为了记录，我尝试在循环中弄乱列表方法（使用切片，不使用切片，仅使用追加，仅使用扩展等），它并没有特别重要，因为主要时间丢失了，因为涉及的迭代次数 - (width * height * numberOfBands)。

任何建议或建议都会非常有帮助。谢谢。

我使用 python 读取了一个 hdf 文件（分层数据格式），其中包含来自卫星的图像数据（名称应该是 MODIS）。python代码能够获取图像的经度、纬度和辐射度。然后使用数据生成图像。

我的老师告诉我目标是把图像放在地图上，我不明白这是什么意思。但她提到要将图像“放置”在地图上，我很可能必须使用“近似最近邻”库（ANN 库）。

图像“放置”在地图上后，应如下所示：

那么在地图上“放置”图像是什么意思呢？谢谢你。

我试图在双时相 RapidEye 多光谱图像上运行 Mort Canty 的http://mcanty.homepage.t-online.de/ Python iMAD 实现。它基本上计算了两个图像的典型相关，然后减去它们。我遇到的问题是图像的大小为 5000 x 5000 x 5（条带）像素。如果我尝试在整个图像上运行它，我会收到内存错误。

使用 pyTables 之类的东西会帮助我吗？

Mort Canty 的代码试图做的是它使用 gdal 加载图像，然后将它们存储在一个 10 x 25,000,000 的数组中。

即使只是创建一个 10 x 25,000,000 的 numpy 浮点数组也会引发内存错误。任何人都知道如何解决这个问题？这是我的第一篇文章，所以任何关于如何发布的建议也将受到欢迎。

问候

我想使用 MRT-MODIS Reprojection Tools 批量拼接 MODIS MCD15A2 数据，这是我的脚本：

它没有用，我收到“警告：CheckMosaicArgs：常规处理：处理镶嵌工具的 (MOSAICINPUT.TXT) 中的 input_filenames 时出错”

我怎样才能解决这个问题？我正在使用windows7。谢谢！

可能的重复：
大卫星图像处理

我正在尝试在双时相 RapidEye 多光谱图像上运行 Mort Canty 的 Python iMAD 实现。它基本上计算了两个图像的典型相关，然后减去它们。我遇到的问题是图像的大小为 5000 x 5000 x 5（条带）像素。当我在整个图像上运行此程序时，我的计算机严重崩溃，我必须将其关闭。

有谁知道什么可以让 python像这样使计算机崩溃？例如，如果我选择每个波段 2999x2999 像素，一切运行正常。

8 GB 内存，I7-2617M 1.5 1.5 ghz，Windows7 64 位。我使用 64 位版本的一切：python(2.7)、numpy、scipy 和 gdal。

先感谢您！

如果名称== '主要': main()