如何从.xyz文件建立一个特定分辨率的数字高程模型?

21次浏览(过去30天)

显示旧的注释

Torcuato点 2020年8月30日

0
链接

直接链接到这个问题

https://es.mathworks.com/matlabcentral/answers/586475-how-i-can-build-a-digital-elevation-model-at-a-specific-resolution-from-xyz-file

评论道: Torcuato点2020年9月8日

答:接受 KSSV

我需要从点云数据集(.xyz格式)生成不同分辨率(0.01 m, 0.05 m, 0.1 m)的数字高程模型。

然后我想绘制得到的dem，以便可视化不同分辨率下的差异。

最后，我想将不同分辨率的dem保存为.xyz文件。

我有xyz数据在本地度量坐标(不是地理)，我加载它们如下:

                         数据=负载(“pointcloud_3D.xyz”）;
                        
                         X = data(:，1);
                        
                         Y = data(:，2);
                        
                         Z = data(:，3);
                        
                         清晰的数据；
                        
                         zzero = z-min (z);”

提前感谢您的帮助。

2的评论
显示隐藏旧的注释

Torcuato点 2020年9月1日

谢谢你，大卫!

这执行得很好，但我需要插值值与IDW ( 逆距离加权 )方法。

你有什么想法吗?!

提前谢谢你

登录评论。

登录回答这个问题。

接受的答案

KSSV 2020年8月31日

1
链接

直接链接到这个答案

https://es.mathworks.com/matlabcentral/answers/586475-how-i-can-build-a-digital-elevation-model-at-a-specific-resolution-from-xyz-file#answer_487394

这取决于您的数据是否被网格化。设A为m*3数据。检查以下选项继续:

                             x = A(:，1);y = A(:，2);z = A(:，3);
                            
                             % %结构化
                            
                             Yi =唯一的(y);
                            
                             [X,Y] = meshgrid(xi,yi);
                            
                             Z =重塑(Z，大小(X));
                            
                             数字
                            
                             冲浪(X, Y, Z)
                            
                             并使用interp2进行插值。
                            
                             % %非结构化
                            
                             dt = delaunayTriangulation(x,y);
                            
                             Tri = dt。ConnectivityList;
                            
                             数字
                            
                             trisurf(三,x, y, z)
                            
                             并使用分散插值。

14日的评论
显示隐藏旧的评论

Torcuato点 2020年9月7日

编辑:KSSV 2020年9月7日

下面的代码和错误消息。无论如何，我只需要分析度量局部坐标中的粗糙高度值(Z)。因此我可以得到- X和Y。

                                  文件=“prova_Geotiff.tif”；
                                 
                                  [path,name,ext] = fileparts(file);
                                 
                                  I = imread(文件);
                                 
                                  Lonmin = -0.8930;Lonmax = 5.5936;Latmin = -1.3711;Latmax = -1.3711;
                                 
                                  x = linspace(lonmin,lonmax,size(I,2));
                                 
                                  y = linspace(latmin,latmax,size(I,1));
                                 
                                  写入geotiff
                                 
                                  R = georasterref(“RasterSize”、大小(我),“LatitudeLimits”(latmin latmax),“LongitudeLimits”[lonmax铂金生产商lonmin]);
                                 
                                  tiffile = strcat(Geotiff，“.tif”);
                                 
                                  geotiffwrite (tiffile, R)

使用map.rasterref.GeographicRasterReference/set时出错。纬度限制(第915行)

latlim元素的值必须是递增的。

setSpatialReferencingProperties错误(第38行)

R.(name) = value;

georasterref错误(第111行)

R = setSpatialReferencingProperties(R，…

Torcuato点 2020年9月8日

Geotiff transpose.jpg

转置I总是有变形，但geotiff被转了90度。

从下面的代码(meshgrid和griddata)开始，是否有可能提取一个带有xyz坐标的.txt文件，以便我能够直接用GIS软件在公制坐标中构建DEM ?

数据= importdata( "input_1cm - Copia.txt" );

X = data(:，1);Y = data(:，2);Z = data(:，3);

Dr = 0.01; %的决议

Yi = min(y):dr:max(y);

[X,Y] = meshgrid(xi,yi);

Z = griddata(x,y, Z, x,y);

pcolor (X, Y, Z)

阴影插值函数

colorbar

登录评论。

类别

数学与优化映射的工具箱数据分析

雷达映射的工具箱数据分析

了解更多数据分析在帮助中心和文件交换

社区寻宝

在MATLAB Central中找到宝藏，并发现社区如何帮助您!

开始狩猎!

翻译的

如何从.xyz文件建立一个特定分辨率的数字高程模型?

2的评论
显示隐藏旧的注释

接受的答案

14日的评论
显示隐藏旧的评论

更多答案(1)

0评论
显示隐藏-1旧的注释

另请参阅

类别

标签

社区寻宝

如何从.xyz文件建立一个特定分辨率的数字高程模型?

2的评论 显示隐藏旧的注释

接受的答案

14日的评论 显示隐藏旧的评论

更多答案(1)

0评论 显示隐藏-1旧的注释

另请参阅

类别

标签

社区寻宝

2的评论
显示隐藏旧的注释

14日的评论
显示隐藏旧的评论

0评论
显示隐藏-1旧的注释