问题标签 [geotiff]

我有一堆大的 GeoTiff 文件（1.4GB，4 纬度 x 8 经度）。

我需要将每个切片切成 1 度纬度 x 1 度长的瓷砖（每个瓷砖只有一点点超龄）。

这是 GeoTiff 文件之一的 gdalinfo 输出：

我希望一切都一样，除了大小。所以每一个都需要被切成 32 块瓷砖，每块 1 度 x 1 度。

例如
，一个将涵盖：-112.0001042、56.0001042 到 -110.9998958、54.9998958（基本上是 -112、56 到 -111、55）。

我看到了像 gdal_retile.py 和 gdal_grid 这样的实用程序，但我没有取得进展。什么是正确的工具/命令行选项？

如何从将每个块geotiff分成常规块（例如，50 x 50 像素）的图像中计算每个块的平均纬度和平均经度。

输入数据只是图像，例如从以下网址下载：http ://eoimages.gsfc.nasa.gov/images/imagerecords/57000/57752/land_shallow_topo_2048.tif

可以使用 gdal 将其打开到 python 中，如下所示：

我有一个 jpg 图像，我想将其转换为 GeoTiff。我用GDAL做了以下事情：

这为 GeoTiff 提供了以下地理数据（使用 gdalinfo）：

我现在想使用 gdalwarp 向文件添加投影，所以我这样做：

gdal info 的输出是：

我期望输出包含 UTM 中的投影信息，如下所示：
坐标系是：
PROJCS["UTM",
GEOGCS["WGS 84",

我究竟做错了什么？

提前致谢。

JC

我有两个 geotiff 图像（说“A”和“B”）在 Matlab 中作为矩阵导入Geotiffread。一个有不同的值，而第二个只有 0 和 255s。我想做的是根据它们的位置将所有 255 替换为另一个图像（或矩阵）中的值。A 和 B 的大小不同，但它们具有相同的投影。

我试过这个：

输出是错误：

否则，我也尝试了合乎逻辑的方法：

什么也没有发生。

有没有办法根据特定值和引用空间中的位置将 A 的值替换为 B 的值？

我正在尝试gdal 1.9.2从源代码构建。该系统是一个Windows 7 64 bits，我在msys环境中工作。

我已经成功构建zlib 1.2.7, jpeg9a, proj 4.8.0,tiff和.geotiff 1.4.0geos 3.4.2

现在，当gdal使用这个进行配置时：

我得到了这个（修剪了一下以显示相关部分）：

所以我知道所有配置都正确，因为我从源代码构建的库被检测到并显示为external。

现在，如果我尝试gdal构建geotiff：

我结束了这条消息：

我已经为此苦苦挣扎了一段时间，但我没有找到解决它的方法。您对如何解决这个问题有任何想法吗？

注意： 一些库是仅静态的：zlib, libgeotiff, tiff.

其他的可用作静态和动态库。

静态库位于：path_to_the_installation/lib.

动态库位于：path_to_the_installation/bin.

以下 MATLAB 脚本在 300x400 数组中生成随机位置，并使用 1-12 的值对这些位置进行编码。 如何将此非空间数组转换为 geotiff？我希望使用 geotiff 输出来进行一些试验分析。任何投影坐标系（例如 UTM）都可以用于此分析。

我尝试使用geotiffwrite()以下实现但没有成功：

这会产生以下错误：

编辑：

我遇到的主要问题是缺乏对该geotiffwrite()函数的输入。我不确定如何处理这个问题。例如，我没有A或R变量，因为数组没有空间参考。只要每个像素都在某处进行地理参考，我不在乎空间参考是什么。这样做的目的是创建一个示例数据集，我可以使用 MATLAB 空间函数进行试验。

我想将 2-D Numpy 数组转换为 3 波段 (RGB) GeoTIFF 文件。你会在下面找到我的一段代码，这就是我打算做的：

1. 获取数组“topo”（由 nan 值组成）
2. 将 numpy 数组转换为 RGB（如此处）
3. 项目清单

我真的不知道如何处理第 2 步和第 3 步中的 nan 值。你能帮我吗？谢谢

我有许多使用朗伯变换的图表。其中一些使用 Tie Point 标签。为此，我已经弄清楚了Lat/Long<=>X,Y转换。

我也有一些使用模型转换标签的图表。到目前为止，我无法弄清楚如何进行点转换。

Geotiff 文档说Model-Coordinates = Matrix * Image Coordinates

一张图表有这个矩阵。

通过 Photoshop 我发现

当我将矩阵（转置的向量）相乘时，我得到

我没有得到我期望的度数（弧度）。也不会反转 X/Y。这让我相信我错过了一些东西。

我从像素到地理坐标转换中遗漏了什么？

以下是其他一些 Geotiff 数据

我有两个 geotiff 格式的栅格，我想在其中执行矩阵乘法以“屏蔽”黑白图像中与二进制图像中的“0”值相对应的区域。下面介绍两个栅格数据集：

1. 1m 空间分辨率的黑白图像（8 位无符号整数；6190 列，7750 行）
2. 二进制图像（0 和 1）（2 位无符号整数；8573 列，11077 行）

此屏幕截图显示图像 2 覆盖图像 1 并突出显示像素大小的差异

这是我的掩码操作的错误结果。您可以看到遮罩图像被缩小到黑白图像的范围，而不是保持它们的空间尺寸。

如何使用不同大小和不同空间分辨率的 geotiff 执行遮罩操作？

这是我用来应用蒙版的方法：

我有一个 10° x 10° 未地理参考的 HDF5，我想将其转换为 geotiff。我尝试了 2 种方法（1）gdalwarp和（2）gdal_translate，但无法找到如何转换文件。

1）gdalwarp

我有以下错误：

2) gdal_translate

我有以下错误：

知道什么是错的或者我应该如何转换这些数据？