我有一些包含曲面的og geojson文件,我需要为每个文件制作一个DEM。这意味着对于每个表面,我想要一个geotiff,其中每个像素是表面的高度。
我尝试使用gdal_grid:
gdal_grid -a linear:radius=0 -co "COMPRESS=LZW" -outsize 15457 27059 surface.geojson DEM.tif
根据文档linear:radius=0
应该使它在点之间使用Delaunay三角剖分来找到每个像素的高度。
这适用于我的很多文件,但是对于其中一些文件,gdal_grid崩溃并显示" gdal_grid.exe已停止工作" -window。根据Windows调试器,它是一个访问冲突异常。我只是不明白为什么。
我找不到文件导致程序崩溃的模式。它不是最大的文件或最小的文件,它们都是非常简单的三维多边形。
所有geojson文件都在UTM32(丹麦),但它并不总是适用于WGS84。
我也在python中尝试过相同的结果,但结果相同:
raster = gdal.Grid(rasterName, featureJson, outputSRS = outputCrs, width = rasterWidth, height = rasterHeight, algorithm = 'linear:radius=0')
异常并不总是在代码中同时抛出,有时其中一个问题文件实际上完成到最后,但DEM在多边形的中间有奇怪的区域,零。
请帮我找出为什么gdal_grid有时会抛出访问冲突异常。