标签: reprojection-error

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

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

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

错误：

---------------------------------------------------------------------------
CPLE_BaseError                            Traceback (most recent call last)
rasterio/_crs.pyx in rasterio._crs._CRS.from_user_input()

rasterio/_err.pyx in rasterio._err.exc_wrap_ogrerr()

CPLE_BaseError: OGR Error code 6

During handling of the above exception, another exception occurred:

CRSError                                  Traceback (most recent call last)
<ipython-input-35-aa94757a1fc2> in <module>
      7 with rasterio.open('/home/DATA/Test.tif') as src:
      8     transform, width, height = calculate_default_transform(
----> 9         src.crs, dst_crs, src.width, src.height, *src.bounds)
     10     kwargs = src.meta.copy()
     11     kwargs.update({

~/anaconda3/envs/elsa/lib/python3.7/site-packages/rasterio/env.py in wrapper(*args, **kwds)
    380     def wrapper(*args, …

我有一个平面标记,我运行SIFT算法来提取特征.然后我运行一个探测器在场景中找到这个标记并再次提取特征.我匹配点并使用OpenCV从匹配对中提取单应性findHomography().

现在,我想用计算的单应性投影在标记中检测到的2D点,以将位置与从场景测量的3D点进行比较,并计算重投影误差.我对像素坐标,厘米,校准矩阵感到困惑,我不知道我应该先做哪些转换.

有人知道这方面的链接还是可以解释方法？

我使用Nister中的5点方法来计算基本矩阵.使用RANSAC和Sampson误差阈值进一步改进了异常值抑制.我随机选择5个点集,估计基本矩阵并评估匹配向量的Sampson误差.Sampson误差低于阈值t(在我所拥有的示例中设置为0.01)的点坐标被设置为内点.对所有基本矩阵重复该过程,并保留具有最佳内点的分数.

我注意到d的大多数值,即sampson误差的向量太大了:例如,如果d的大小是(1x1437),如果我做g = find(abs(d)> 0.01); 长度(G)

然后长度(g)= 1425这意味着只有7个值是这个阈值的内部值,这是不正确的!

如何设置门槛？如何解释Sampson错误值？

请帮帮我.谢谢

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

有些人建议使用外部库，例如 ogr2ogr，但这对我来说不是一个好的选择，因为我需要即时完成。我在这里找到了一篇文章 Transform/Project a geometry from one SRID to another，但这个转换为纬度和经度并且具有一些有限的功能。

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