3

抱歉，如果我错过了明显的...

我正在用 rgl 绘制一个 3d 曲面。我的代码是

library(rgl)
dem1 = read.table(file="file.txt",skip=5,header=F,na.strings="0.")
dem = dem1[order(dem1$V1,dem1$V2),] 
z = matrix(dem$V3,nrow=1250,ncol=1250)
is.na(z) = (z < 200)

#create x y dimensions
x=4*(1:nrow(z))
y=4*(1:ncol(z))

open3d()
bg3d("white")
persp3d(x,y,z)

这给出了这张地图（添加了颜色以更好地查看功能，即使我没有将代码放在上面）

在此处输入图像描述

问题是，无论我对这张地图做什么，它都是颠倒的，即 x 应该是 y，而当前的 y 从西 (0) 到东 (5000) 应该是相反的，这样高架特征实际上应该是左下而不是右下。

我使用与此脚本相同的文件绘制了一个非常简单的等高线图

dem=read.table("file.txt",header=F,skip=5,na.strings="0.")
library(lattice)
contourplot(dem$V3 ~ dem$V1+dem$V2)

这使

在此处输入图像描述

并得到右轴和左下角最高的区域，正是它应该在的位置，所以数据没有问题。

我解释了数据在这里的样子以及为什么我觉得需要重新排序

dem = dem1[order(dem1$V1,dem1$V2),]

奇怪的是我是否使用了上面的命令，3d 表面图看起来完全一样，这让我想知道代码是否真的使用了使用order命令创建的“dem”数据集，或者它是否仍在使用原始的“dem1” " 它从文件中读取并且顺序错误的数据。

我很高兴根据请求发送数据或将其放在可以看到的地方，但我无法在此处复制它，因为它是 1250 行 x1250 列。

提前致谢！

4

1 回答 1

2

问题在于z高程矩阵的创建。R 在创建矩阵时按列填充矩阵。正是这种由柱子填充的情况正在重新排列相对于彼此的高程。矩阵是正方形的事实使情况更加复杂。x如果矩阵不是方阵，则y和之间的关系z会发生更显着的变化，而不仅仅是被翻转。

解决方案是让 R 按行填充矩阵，例如定义z使用：

z <- matrix(dem$V3, nrow=1250, ncol=1250, byrow = TRUE)
于 2011-08-24T16:06:16.327 回答