问题标签 [reprojection-error]

我正在尝试按照这个 youtube 教程在 QGIS 中制作热图。但是，我的数据是度数而不是米数。我正在使用的数据的文档说 CRS 应该是 EPSG:4326 - WGS 84，这是我导入的。我还需要重新投影吗？

，但我试图获得米而不是度数。我使用的数据来自本教程。

I have two corresponding image points (2D) visualized by the same camera with intrinsic matrix K each coming from different camera poses (R1, t1, R2, t2). If I triangulate the corresponding image points to a 3D point and then reproject it back to the original cameras it only closely matches the original image point in the first camera. Can someone help me understand why? Here is a minimal example showing the issue:

The reprojection error in the first camera is always good (< 1px) but the second one is always large.

我有大量的栅格（ASCII 文件），我想： 1. 从 Lamberts Equal area 重新投影到 WGS 84 2. 将生成的 WGS 84 栅格裁剪到一定程度 3. 将生成的栅格写入目录

我知道堆叠栅格和重新投影会遇到内存问题。因此，我尝试了一个 for 循环，该循环遍历每个光栅，重新投影作物并将其保存到输出目录。即便如此，我也遇到了内存问题。如果我只是裁剪和掩盖栅格，代码工作得非常快，当我想重新投影时会出现问题。这是我的代码

在geopandas 文档中它说

AGeoDataFrame还可能包含其他具有几何（形状）对象的列，但一次只能有一个列是活动几何。要更改哪一列是活动几何列，请使用该set_geometry方法。

如果目标是将这些不同列中的几何数据灵活地重新投影到一个或多个其他坐标参考系统，我想知道如何使用这样的 GeoDataFrame 。这是我尝试过的。

第一次尝试

到目前为止，一切都很好。如果我想设置geom2为几何列，事情就会出轨，并重新投影该列：

好的，所以，显然，在更改用作几何图形的列时，.crs属性并没有重置为其原始值 - 看起来，crs 没有为各个列存储。gdf如果是这种情况，我看到对此数据框使用重新投影的唯一方法是回溯：开始->选择列作为几何图形->将gdf重新投影到crs_new->使用/可视化/...->将 gdf ​​重新投影回 crs_lonlat --> goto start。如果我想在一个图中可视化两列，这将不可用。

第二次尝试

我的第二次尝试是，crs通过将上面脚本中的相应行更改为：

然而，很快就很清楚，虽然初始化为 a GeoSeries，但这些列是 normal pandas Series，并且没有.crs相同的属性GeoSeries：

我在这里缺少什么吗？

唯一的解决方案是事先决定“最终”crs - 并在添加列之前进行所有重新投影？像这样...

...然后，当需要另一个 crs 时，gdf从头开始重建整个？这不可能是打算使用的方式。

我正在尝试使用 Google Earth Engine (GEE) 中的 reduceResoluton().reproject() 函数重新缩放图像（30 到 100m 空间分辨率）。我想要的是将 30 m 图像的平均值计算到新的 100 m 分辨率图像中。我试过以下

我使用了一个 Landsat 场景，它提示瓦片错误“瓦片错误：图像计算的输出太大”： 错误 我错过了什么吗？是否有另一种方法来执行此任务？非常感谢大家

显然，NOAA 和 NWS 对他们的一些降雨数据使用了非传统的投影方式，并且在将其投影为其他用户的传统格式方面没有提供很多帮助。我在让栅格覆盖美国部分地区方面取得了一些成功，但仍然不太正确。

我希望有人可以帮助我破译我所缺少的内容并纠正这些数据的投影。

您可以在此处找到有关此数据的更多信息：https ://polyploid.net/blog/?p=216 https://water.weather.gov/precip/download.php

我设法从 HDF 文件中提取 MODIS 土地覆盖数据并将其放入栅格中。

我想把它放到一个数据框中，但是从原来的正弦曲线中重新投影+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs，我认为。

对于与我正在分析的其他数据集兼容的普通椭圆体/WGS84。

这是我尝试过并且似乎有效的方法：

但是，当我随后将我的土地覆盖数据放入这个新投影中的数据框时，所有土地覆盖值都变为 NA。

这就是数据框在正弦投影中的样子（4 是土地覆盖分类）

通过我的重新投影，它看起来像这样：

关于我做错了什么或如何让土地覆盖坐标正确重新投影的任何建议？

干杯！！！！

我还读到有一个 NASA MODIS 重投影工具，但它不再存在/可用。有人知道吗？

我需要将 .tiff 文件从 EPSG:4326 - WGS 84 - Geographic 重新投影到 EPSG: 32638 (UTM Zone 38N)。当我运行下面的代码时，输​​出文件是一个重新投影的 .tiff 大小的原始文件的一半。

我正在尝试使用以下代码将 .tiff 文件重新投影到 EPSG:32638。

我安装的版本：Rasterio 版本 1.1.5 Numpy 版本 1.18.1

这是我正在使用的代码：https ://rasterio.readthedocs.io/en/latest/topics/reproject.html（重新投影 GeoTIFF 数据集）

错误：

我一直在处理栅格数据，代码运行良好。意外地使用以下命令重新投影时：
reproject()
我得到“wkt(projfrom) 中的错误：找不到函数“wkt”。我正在使用... CRS 参数：+proj=utm +zone=5 +datum=WGS84 +units =m +no_defs
我也试过 projectRaster() 并得到类似的错误可能是 rgdal 包的问题？还有其他人得到这个吗？谢谢！