0

我正在使用一个从激光雷达数据生成数字高程模型的程序。代码是使用 PDAL 2.3.0 库在 python 3.8 中编写的。我需要将重投影从 EPSG 3395 更新到 EPSG 4326。当我使用 EPSG 4326 运行管道时，会生成一个全黑的 geotiff，但没有数据。不知道我在这里做错了什么。我已经附加了传递给 pdal 的 JSON 数据。

 dtm_json = json.dumps([file_path,
                {
                    "type": "filters.reprojection",
                   # "out_srs": epsg
                    "out_srs": "EPSG:4326"
                },
                {
                    "type": "filters.smrf",
                    "window": self.window_size,
                    "slope": self.slope,
                    "threshold": self.threshold,
                    "cell": self.cell_size,
                    "where": "(NumberOfReturns > 0 && ReturnNumber > 0)"
                },
                {
                    "type": "filters.range",
                    "limits": "Classification[2:2]"
                },
                {
                    "type": "writers.gdal",
                    "filename": self.output_file,
                    "output_type": "min",
                    "gdaldriver": "GTiff",
                    "window_size": 3,
                    "resolution": 1.0,
                    "data_type": "float32",
                    "nodata": "nan"
                },
                {
                    "type": self.FILTERS_INFO,
                }
            ])
4

0 回答 0