问题标签 [reprojection-error]

我尝试从以下位置重新投影 OpenLayers 3 中的 Geojson 文件：

• Geojson 链接
• 有了这些 JS 在头脑中：

- 使用此代码使用 newProjCode 设置新的 ol.View：

• 但是在设置了新的 ol.View 之后，我只看到了白屏。

问题出在哪里？

下面的代码是 OpenCV 的相机校准部分。我只想改变一点：变量err是使用欧几里得距离公式计算的sqrt(x-coordinate diff^2+y-coordinate diff^2)，但我只想考虑 x 坐标差异。所以我重写了它

imagePoints[i]但是，我收到错误，因为我无法访问&的向量的 x 坐标值imagepoints2。

err有人可以帮助仅考虑 x 坐标增量来重写变量吗？

从给定的世界点（原始坐标）、内在矩阵、旋转矩阵和平移向量计算重投影点、重投影误差和平均重投影误差的过程是什么？

是否有任何内置的opencv函数或者我们应该手动计算？

如果我们必须手动计算，那么获得重投影点的最佳方法是什么？

我一直在使用 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 中提取了相同的图像。

我从http://www1.kew.org/gis/tdwg/index.html下载了 TDWG level4 形状。此形状具有 WGS84 投影。我需要将此层重新投影到欧洲 LAEA 或 CRS("+init=epsg:3035")。

我使用了以下代码：

并得到以下错误：

.spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, : Polygons 中的错误 37 Polygon 1 分另外：警告消息：在 .spTransform_Polygon(input[[i]], to_args = to_args, from_args = from_args, : 361 个投影点不是有限的

有什么建议可以解决这个问题吗？

Step我正在尝试对使用同一台相机拍摄的一系列立体图像（类）执行捆绑调整（BA） 。

每个Step都有左右图像（校正和同步）、生成的深度图、左图像的关键点+描述符和 2 个 4x4 矩阵 - 1 个用于局部（图像平面）到全局（3D 世界），以及它的逆矩阵（T_L2G 和 T_G2L分别）。

步骤相对于第一张图像注册。

我正在尝试对结果运行 BA 以改进转换，并且我正在尝试使用 PBA ( https://grail.cs.washington.edu/projects/mcba/ )

设置摄像头的代码：

然后，我通过在所有图像对上运行并匹配（当前使用 SURF）来生成全局关键点。

然后，生成 BA 点数据：

然后，运行 BA：

然后，在我遇到问题的地方，从 PBA 中提取数据：

一切都按照我期望的方式运行。PBA 报告的初始错误相对较小（目前使用少量成对注册图像进行测试，因此错误不应太大），运行后报告的错误较小。（快速收敛，通常少于 3 次迭代）

然而，当我使用新发现的转换来转储关键点时，云似乎彼此之间的距离更远了。

（我也尝试在 T_G2L 和 T_L2G 之间切换以“让它们更接近”。不起作用）。

我想知道使用它是否缺少一些东西。

我使用 Matlab Stereo camera app 校准了立体相机，我的重投影误差接近 0.19 像素。但是当我使用相同的校准参数来找出对象的 3D 坐标时，重投影误差在较高的一侧接近 7-8 个像素。我正在做的是找出左右图像中对象的中心分别表示（xl，yl）和（xr，yr），然后使用这些点来找出重投影误差。我做错了什么？为什么重投影误差与校准重投影误差相比如此之高？

谢谢。

我有一个来自 sentinel1 数据的投影 tiff 文件（s1a-ew-grd-hh-20180208t071218-20180208t071323-020512-023164-001.tiff）（我们称之为 test.tiff），其投影为 WKT 格式：

这应该相当于北极极地立体（epsg 3995）。我的目标是通过插入一些 GCP 并扭曲它来稍微调整图像。为了仔细检查这个方法，我从一个简单的案例开始，我从投影图像 (test.tiff) 中获取角坐标：

我据此设置GCP。因此：

然后我使用 -tps 选项扭曲图像以使用 GCP 信息强制扭曲。

人们会期望 test.tiff 和 gcp_warped.tiff 文件应该重合，但不幸的是情况并非如此：

在上述过程中有什么我想念的吗？还是有其他解决方法？

我有每日温度和降水数据的 netcdf 文件 如何更改 netcdf 的投影？

我曾尝试使用 raster 函数来读取文件并使用 projectraster 重新投影它们？

代码_使用

错误_

（函数（类，fdef，mtable）中的错误：无法找到签名“字符”的函数“res”的继承方法另外：警告消息：在min（dim（to）[1：2]）中： min 没有非缺失参数；返回 Inf

我正在尝试在 SQL Server 中进行一些重新投影，但不幸的是我找不到任何解决方案。我知道这在 SQL Server 中没有实现，所以想知道是否有任何快捷方式。

有些人建议使用诸如 ogr2ogr 之类的外部库，而这对我来说不是一个好选择，因为我需要即时进行。我在这里找到了一篇文章 Transform/Project a geometry from one SRID to another，但是这个转换为 lat 和 long 并且具有一些有限的能力。

我正在寻找一种解决方案，允许我在任何 srid 中进行重新投影，这与ST_Transformpostgresql (postgis) 中的非常相似