问题标签 [rasterio]

我有一个在使用AWS Deep Learning AMI (DLAMI)创建的 Ubuntu 16.04.4 LTS 上运行的 AWS EC2 p2.xlarge 实例。我正在使用 keras/Tensor Flow conda 环境：

我正在尝试使用以下命令（源）在 Keras - Tensor Flow AMI 安装之上安装 Rasterio 和 GDAL：

GDAL 安装似乎有效：

但是，当我尝试在 Python 中导入 rasterio 时，会产生以下错误：

如何清除此错误以便在 Python 中导入和使用 rasterio？

我想以特定格式堆叠大约 700 个 Landsat 图像作为YATSM连续变化检测和分类 (CCDC) 算法的输入。我尝试使用这个 bash 脚本4_stack_TIFF.sh。在这个脚本中，堆叠是使用这个文件完成的landsat_stack.py。我收到两个错误：

• 首先与缺少一个模块有关module load batch_landsat。我可以找到共享存储库

• 第二个是相关的landsat_stack.py，我得到以下信息： Stacking to extent: Traceback (most recent call last): File "./landsat_stack.py", line 936, in <module> sys.exit(main()) File "./landsat_stack.py", line 920, in main utm, resume, fformat, creation_opts)) File "./landsat_stack.py", line 757, in landsat_stack print('\tUpper Left: {ulx},{uly}'.format(ulx=extent[0], uly=extent[1])) TypeError: 'NoneType' object is not subscriptable Done!

问题：

1. 是否有人使用过 YATSM 并在陆地卫星数据预处理中面临同样的问题？
2. 如果不; 你能建议我使用 Rasterio 或 Rio 的解决方案吗？

格式：

我正在尝试使用光栅从 .TIF 中提取坐标和像素值数据。 在最近更新的文档（2018 年 4 月）中，似乎有一个 .xy() 函数来获取像素的空间坐标，例如：

但是，当我运行相同的示例时：

我收到错误：“AttributeError：'rasterio._io.RasterReader' 对象没有属性 'xy'”。

我是在误解文档还是这类似于未来版本的预文档？在当前版本/官方文档中是否有另一种我没见过的方法？

谢谢！

我正在尝试从使用 python3.4 在 Windows 7 - 32 位上开发的 python 应用程序制作可移植应用程序。我的第一个想法是将 py2exe 与 ThinApp 等效项结合使用，但由于我使用的不是纯 python 库的 rasterio，因此 py2exe 无法找到一些 DLL。如果我尝试运行生成的可执行文件，我会收到此错误：

我读到我必须在 dist 目录中手动添加这些 DLL，但我想知道如何获得所需 DLL 的列表？

我正在使用以下软件包：

我有以下存储数据的对象：

这看起来如下：

我的 shapefilegeopandas.GeoDataFrame代表一个多边形。

如下所示：

在另一个之上绘制一个，它们看起来像这样：

我的问题是，如何xarray.DataArray通过检查经纬度点是否位于存储为 a 的 shapefile 内部来掩盖geopandas.GeoDataFrame？

 所以我只想要那个形状文件中的降水值（毫米/天）。

我想做如下的事情：

或者

编辑 1

我曾考虑过使用该rasterio.mask 模块，但我不知道输入数据需要采用什么格式。听起来好像它做了正确的事情：

“使用输入形状创建一个蒙版或填充的数组。像素被蒙版或设置为输入形状之外的 nodata ”

此处从 GIS Stack Exchange 转贴

我正在使用 rasterio 从 geotiff 中切出我的 shapefile 中勾勒出的湖的形状。当使用带有 geotiff 的单一湖泊的 shapefile 时，这可以完美地工作，但在使用包含许多不同湖泊轮廓的不同 shapefile 切割单个湖泊的 geotiff 时就不行了。我想知道我的代码编写方式或我正在使用的文件是否有问题？

这是我的代码：

我正在尝试重新创建这个光栅示例：

我检查了rasterize( print(result)) 的结果，这似乎是正常的，并且结果是 2x2 像素正方形：

但我得到一个完全黑色的图像（nodata图像）。

我在 Windows 7 上安装了 Anaconda Python 3.6。我经常使用 GeoPandas，它运行良好。昨天我尝试安装 rasterio 和 georasters，他们似乎使用建议的方法安装得很好

然后使用安装georasters

并收到安装成功的消息。当尝试导入它们时，它们不起作用 -import rasterio或上有一条错误消息import georasters。深入研究，我发现 geopandas、fiona 和 GDAL 现在也坏了，在尝试使用 conda 卸载/安装软件包 1 小时后，我不确定问题和解决方案是什么。根本问题似乎是GDAL，但我不确定。

导入 GDAL 会给出以下错误消息

有谁知道修复或如何干净地重新安装这些软件包？

我尝试将 tiff 转换为 netcdf 文件。错误是说索引错误：

错误是索引错误：数据数组的大小不符合切片。有人可以看看并帮助我在哪里做错了。非常感激！

感谢您光临这里。我目前正在做一个项目，但我对 python 还不是很擅长。所以我希望对这个问题有一些建设性的反馈。

我有一个包含核心样本的数据集，所有样本都存储有样本 id、纬度、经度、内容和与此问题无关的其他数据。

现在我已经导入了这个数据集并按照我想要的方式对其进行了切片。对于图像，我使用 rasterio 模块打开 2 个覆盖该地区的卫星图像。我正在使用 utm 模块在 latlong->UTM->Pixel 值之间来回转换（这似乎也在某些点给我带来了奇怪的坐标）。令人恼火的是，两张 Sentinel-2 图像正好穿过地图的中心。当我在采样的顶部做边界框时，这是一个问题，因为我需要从该区域中提取 10x10 像素的切口。这导致很多样品没有得到适当的切割。

所以我想为什么不将这两个图像合并成一个大的矩形位。但我仍然需要使用 UTM 坐标保留元数据。

你会建议我如何进行。能以简单的方式完成吗？我忽略了另一个角度吗？

感谢您的时间。