.nc地形數(shù)據(jù)的python轉(zhuǎn)換實(shí)現(xiàn)tiff

之前做的一個(gè)海底地形的數(shù)據(jù)導(dǎo)入arcgis的問題，由于是全球的地形數(shù)據(jù)（精度為15弧秒），從GEBCO下載獲取，數(shù)據(jù)達(dá)到12G，在arcgis中轉(zhuǎn)換也可以，但是實(shí)在是太慢，我的電腦i7-7700的轉(zhuǎn)了一晚上沒成功，后來看了好多arcgis的官方文檔和網(wǎng)上的轉(zhuǎn)化教程，改寫成下面這一小段，速度很快，大約兩小時(shí)即可完成。
我印象中有兩篇參考的非官方文檔給我提供了很大幫助，可惜找不到了，在此不能引用，但表示感謝！

#coding=utf8 # Import system modules import arcpy from arcpy import env from arcpy.sa import * # Input data source arcpy.env.workspace = "E:/ALL_change" arcpy.env.overwriteOutput = True # Set output folder OutputFolder = "E:/ALL_change" # Loop through a list of files in the workspace filename = "GEBCO_2019.nc" #for filename in NCfiles: #print("Processing: " + filename) inNCfiles = arcpy.env.workspace + "/" + filename fileroot = "npr"+filename variable = "elevation" outRaster = OutputFolder + "/" + fileroot # Process: Make NetCDF Raster Layer arcpy.MakeNetCDFRasterLayer_md(inNCfiles, "elevation", "lon", "lat", variable, "", "", "BY_VALUE") # Process: Copy Raster arcpy.CopyRaster_management(variable, outRaster + ".tif", "", "", "", "NONE", "NONE", "")

總結(jié)

以上是生活随笔為你收集整理的.nc地形数据的python转换实现tiff的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。