问题标签 [topography]

For questions regarding programming in ECMAScript (JavaScript/JS) and its various dialects/implementations (excluding ActionScript). Note JavaScript is NOT the same as Java! Please include all relevant tags on your question; e.g., [node.js], [jquery], [json], [reactjs], [angular], [ember.js], [vue.js], [typescript], [svelte], etc.

0 投票
1 回答
62 浏览

r - 基于已知重叠将地形剖面拼接在一起

我有一系列地形剖面扫描,我想将它们组合起来创建一个连续的剖面。唯一的问题是每次扫描可能或可能不是从不同的高度进行的,因此虽然不同的文件在覆盖区域方面有相当多的重叠,但不同的数据可能没有共同的参考点绝对高度。

以下是 4 种不同的扫描。每次扫描包含大约 30 次测量,最后几次测量代表新数据,其余的与前一次扫描重叠。第一次扫描包含唯一已知的绝对值,因此第一次扫描是“黄金标准”。第二次扫描恰好是从相同的高度进行的,因此重叠(几乎)完美匹配,并且仅在前一次扫描中增加了 4 个新点。第三次和第四次扫描是从不同的高度拍摄的,所以虽然重叠覆盖了相同的区域(相对),但我不能简单地将它拼接到前两次扫描上。

使用 R,有没有办法将这 4 次扫描拼接在一起以制作连续的地形剖面?绝对高度需要基于第一次扫描,每次连续扫描都被缝合到前一次扫描上。IE-Scan2 被拼接到 Scan 1 上添加 4 个数据点,然后来自 Scan 3 的新数据被添加到 Scan1 和 Scan2 的组合中,然后来自 Scan4 的新数据被添加到 Scans 1,2 和 3 的组合中, 等等....

我假设有一种方法可以通过匹配扫描之间的大重叠来规范化所有数据,使用某种模式识别来确定 Scan3 与 Scan1 大约相差 8 个单位,而 Scan4 大约相差 11 个单位。但请注意,我的数据中有一些“噪音”,重叠的模式并不完美。

最终结果应包含包含所有 4 次扫描的完整地形剖面,并在实际数字不同时进行某种调整。类似于以下内容:

0 投票
1 回答
1035 浏览

r - 在 R 中使用 plot3D {rasterVis} 在同一窗口中显示多个 3d 图

我有兴趣制作两个并排显示在同一个 Xquartz 设备窗口中的 3d 地形图。使用 rgl 包在同一个窗口中显示两个 3d 图很简单——文档中有很多示例 usingmfrow3d()和其他方法。

但是,我需要使用plot3D()rasterVis 包中的函数。我使用plot3D()rasterVis 而不是 rgl的原因surface3d()是我需要使用drape参数 inplot3D()将栅格中的值显示为 3d 地形图上的颜色(并且该栅格的值与创建该栅格的值不同)图上的 z 轴)。如果有人对 rgl 函数中类似于悬垂的东西有提示,我也会对此感兴趣!

当我尝试mfrow3d()与该plot3D()功能一起使用时,它会显示一系列空白设备窗口,而不是在同一个窗口中并排显示两个图。

以下是 rasterVis 文档中使用 plot3D 制作地形图的一些代码:

在这里,我尝试使用 mfrow3d 并排绘制 2 个相同的火山图,一个蓝色,一个红色,我改编自 rgl 文档:

我正在尝试使用 rasterVis 的 plot3D 做些什么?

0 投票
0 回答
138 浏览

r - 如何将地形图添加到 shapefile?

我想知道如何将这种地形信息(海拔)添加到我在这里拥有的圣克鲁斯岛和所有加拉帕戈斯岛的地图上。

你有什么主意吗?你在 R 中有一个包可以做到这一点吗?

谢谢,

0 投票
1 回答
776 浏览

gpgpu - Delaunay三角剖分内整数坐标的双线性插值

我有一个由大约 100 万个三角形组成的平面 Delaunay 三角剖分。每个顶点都标有几个标量指标 [1],我希望在同一个规则网格上看到每个指标的快速、简单插值。作为参考,我的三角形的并集覆盖了大约 1000 万个具有(整数)坐标的网格单元。[2]

当我说简单时,我的意思是简单。双线性就好了!我的理解是,这(a)基本上是 GPU 为生而做的,(b)可能是无数家庭作业的主题。我自己是公共卫生领域的政府研究员,所以这不是我的功课。:-)

在我缓慢但正确的参考实现中,我可以在大约 10 分钟内计算出以下内容:

对于每个三角形 T:

  1. T的边界框内所有(整数)笛卡尔坐标的集合G;
  2. G中每个(x,y)的重心坐标(u,v,w);
  3. 拒绝不全为正的(u, v, w)——即在T内部;
  4. T 中每个剩余坐标的加权和 (u z_1 + v z_2 + w*z_3),其中 z_1、z_2 和 z_3 对于给定度量 [1] 是 T 顶点处的标量值。

我真的需要步骤 1-3 快速;第 4 步是微不足道的,但这是我的最终目标。理想情况下,解决方案将采用以下任一形式:

  • 一个经过适当许可(GPL 可以)的库,带有一个非常简单的 API;或者
  • 一个足够清楚的解释,很明显中级程序员如何用 Fortran、R、Python 或 C 对其进行编码。

此任务的经典表述是“TIN 到 DEM”地形建模作业。但现在似乎更需要相反的情况(?)

一些基本的清理,比如当一个点恰好落在由 2 个以上三角形共享的边或顶点上时删除重复项,也是可以的。

非常感谢您的时间和关注。下火车后,我将根据建议清理格式和编辑!

脚注:

[1] 海拔、温度和湿度。[2] 整数,即它们在 UTM 网格上的间距为 20x20m。所以只需按 20 缩放。

0 投票
1 回答
4504 浏览

python - Matlab 的 Demcmap 的 Python 等效项(海拔 +/- 适当的颜色图)

我正在寻找一种方法来为matplotlib.

cmap ' terrain' 看起来不错,但颜色缩放不是基于零(即,如果比例为 0->5000m,则 0->1000m 范围可能是蓝色阴影,您会假设它低于海平面)

例如: 在此处输入图像描述

Matlab 的等效函数是: demcmap

让 matplotlib 在零高程标记周围移动地形颜色图的绿色/棕色和蓝色的最佳方法是什么?

0 投票
0 回答
69 浏览

image - 从等高线函数中分离一条等值线

我在 RGB 图像上使用了轮廓函数(在分离了三个通道之后),现在我想分别处理每个等值线。

如何在单独的矩阵中获取每个等值线信息?

我有 9 个等值线,我想要 9 个矩阵。

我试过这个(从 Matlab 轮廓函数中选择等值线),但它没有用。

0 投票
1 回答
341 浏览

python - LandXML to Topography.ByPoints 不能按预期工作

我是新来的,所以我希望这个线程不存在。

我从 Civil3D 导出​​了一个 LandXML 文件。在这个文件中列出了 3 个表面。现在我的目标是将其导入 revit 并将其转换为地形。提取所有点和面并生成网格是没有问题的。但是使用 Topography.ByPoint 节点,我得到了一个看起来不正确的地形。

[网格和地形][1]

我已经尝试了很多。但现在我不知道出了什么问题。有人有想法吗?

0 投票
1 回答
426 浏览

python - 如何编写脑电图的地形函数

我正在编写一个硬件神经反馈应用程序。哪个 wotks 像这样: https ://www.youtube.com/watch?v=pjCghiq5FoQ

在这种情况下,它是一个实际上不起作用的演示。我有来自 EEG 的数据,我想通过 LED 将它们绘制在表面上。该技术类似于此输出: http ://www.fieldtriptoolbox.org/_media/tutorial/natmeg_temp/natmeg_freq11.png?w=400&tok=6e5a3c

但我需要自己编写,因为我需要点亮 LED 来显示 2D 图像。基本上我不知道从哪里开始。

我的目标是通过每个 LED 可视化每个 EEG 通道的频谱密度,如您在 youtube 上看到的演示。我愿意提供任何帮助,即使是理论上的帮助。

我知道我需要 x、y、z 中的信号和电极位置,例如雕塑上 LED 编号的标识。¨

谢谢,

迈克尔

0 投票
1 回答
61 浏览

matlab - 计算每个三角形的地形与矩阵的角度

我有一个包含地形数据的矩阵,比如说几座山。我想知道每个数据点与垂直线的角度信息。这里有两个例子:

  1. 如果我考虑一个完全平坦的山脚附近的地方,我的角度是 90°(与垂直线成 90°)。
  2. 如果我在山的最陡点,我有一个较低的角度,比如说 50°。

为了计算这个,我想我必须连接所有的地形数据,以便(至少)三个近像素形成三角形。之后我必须计算这个三角形的角度。

我可以使用现有的算法吗?

0 投票
1 回答
273 浏览

c++ - 窗口中出现 glTexImage2D 切片和对齐问题

我正在制作自己的地形图,并且一直在使用来自 NASA 的 .hgt 文件。

我正在使用加载文件

然后使用以下命令将它们添加到内存位图中:

然而,这会导致图像的角被切开,就像这样 切片图像

在 loadTextureImage() 中使用带注释的版本看起来略有不同,但是切片的角相同。 切片图像 2

有没有人暗示发生了什么?我尝试使用图像作为纹理,使用 stbi 库加载,效果很好,所以我不确定哪里出了问题。

(图像的坐标是 N10E099)