.nc地形数据的python转换实现tiff

.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", "")  
  • 1
    点赞
  • 3
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值