.nc地形数据的python转换实现tiff
之前做的一个海底地形的数据导入arcgis的问题,由于是全球的地形数据(精度为15弧秒),从 GEBCO下载获取,数据达到12G,在arcgis中转换也可以,但是实在是太慢,我的电脑i7-7700的转了一晚上没成功,后来看了好多arcgis的官方文档和网上的转化教程,改写成下面这一小段,速度很快,大约两小时即可完成。
我印象中有两篇参考的非官方文档给我提供了很大帮助,可惜找不到了,在此不能引用,但表示感谢!
#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", "")