我正在尝试使用GDAL将netCDF文件转换为EPSG:3857以便与Mapbox一起使用,这将是.nc到.nc的转换,而不是栅格化。我可以使用GDAL或其他方法来完成此操作。在将此数据发送到控制台应用程序之前,必须重新投影此数据,而这个过程需要几周时间来找到解决方案 - 我认为这很简单。
我正在处理着色卫星数据。当组合和处理三个.nc文件(蓝色、红色和红外线)时,会创建彩色图像。在下载这3个文件(从Amazon AWS)后,一个Python控制台应用程序进行处理,并将.jpg倒入同一文件夹中。该应用程序的源代码位于这里,以验证数据。(由于文件分辨率超高,因此速度很慢)。
我尝试的代码是:
gdalwarp -t_srs EPSG:3857 test.nc test-projected.nc
然而,已经尝试了其他几种变化,但都无效。
我不是专业人士,但我是否应该使用gdalwarp来做到这一点?我只想更改投影 - 仅此而已,以便Python应用程序仍然可以使用数据。它必须能够使用重新投影的文件创建.jpg。
以下链接是需要转换的数据样本:
AWS上的 .nc 文件 > 色彩通道 1(蓝色 1km 分辨率)
AWS上的 .nc 文件 > 色彩通道 2(红色,0.5km 分辨率更高且文件大小更大)
AWS上的 .nc 文件 > 色彩通道 3(红外线 - 用作绿色)
此外,其他人使用类似的投影通过pyproj模块在https://github.com/blaylockbk/pyBKB_v2/tree/master/BB_GOES16中完成了这项工作。(我的必须是EPSG:3857,以便与Mapbox一起使用)。如果Python代码被修改为一次完成所有操作,那就太好了。我正在发放赏金,这是最后的希望。
我不懂Python,所以大部分时间都在尝试GDAL - 然而,将可行的Python代码添加到我的源代码中以实现预期结果(或有效的GDAL脚本)将赢得赏金。