怎样将Esri ASCII Grid转bin文件
发布网友
发布时间：2023-06-11 02:20
共1个回答
热心网友
时间：2023-10-09 04:43
问题
ASCII grid由matlab代码生成，存储于 LDAS.data 目录。ASCII grid内容如下：
ncols 32
nrows 32
xllcorner -88
yllcorner 37.75
cellsize 0.125
NODATA_value -9999
0.000000 0.000000 ...
各ASCII grids命名有以下规律：
nbwi3.[yyyy].[mm].[dd].[hh]h.txt
其中[yyyy]是4位年份，[mm]是月份，[dd]是2位日期，[hh]是24制的小时。比如，nbwi3.2001.07.08.15h.txt等。
生成的ESRI Binary grid有以下命名规律：n[yy][mm][dd][hh]，其中[yy]是2位年份，[mm]是月份，[dd]是2位日期，[hh]是24制的小时。如n01070814。注意，ESRI grid的命名和长度是有一定要求*的，所以我们必须要选择合适的，并能区别数据的命名方式。
我们需要做的是将数百个ASCII grids转换成对应的ESRI binary grids。
实现
脚本编写环境我们这里采用免费的PyDev/Eclipse。Python脚本存储在与LDAS.data并列的src目录下。事实上，可以使用任何文本编辑软件来用脚本，但Eclipse提供很多编程增强功能，如智能提示等。
import arcgisscripting
import os
datadir=os.path.abspath(os.path.join(os.getcwd(),"../LDAS.data"))
#print datadir
children=os.listdir(datadir)
#only ".txt" files
txtfiles=[]
grdnames=[]
for child in children:
root,ext=os.path.splitext(child)
if (ext=='.txt'):
txtfiles.append(child);
grdname='n'+root[8:19]
grdname=grdname.replace('.','')
grdnames.append(grdname)
# grid names now ready
gp=arcgisscripting.create()
gp.workspace=datadir
# begin iteration
for i in range(len(grdnames)):
gp.ASCIIToRaster(txtfiles[i], grdnames[i], "FLOAT")
print txtfiles[i]
首先，需要通过import arcgisscripting，将arcgisscripting导入Python。arcgisscripting事实上是位于C:\Program Files\arcgis\bin目录下的arcgisscripting.dll。如果在运行时，提示找不到arcgisscripting库，则需要配置Python库路径。
库导入成功后，gp=arcgisscripting.create() 创建script对象，此后通过“gp.命令”的方式调用script支持的命令。一般来讲，可以在arcgis command window里使用的命令，都会有对象的script调用方式。比如，
ASCIIToRaster的命令语法是 ASCIIToRaster_conversion <in_ascii_file> <out_raster> {INTEGER | FLOAT}，对应的脚本语法是ASCIIToRaster_conversion (in_ascii_file, out_raster, data_type)。我们可以在ArcGIS Desktop帮助文档里找到。
此外，gp还有一些专用的script命令，比如ListRasters、GetMessage等，大家可以在ArcGIS Desktop帮助文档查询Geoprocessor object。这些命令对于进行批处理迭代时十分有用。
在以上代码中，与geoprocessor有关的代码有，
gp=arcgisscripting.create()
gp.workspace=datadir
gp.ASCIIToRaster(txtfiles[i], grdnames[i], "FLOAT")
workspace设置了当前的工作目录，相当于在command window里执行 “workspace 目录”。
ASCIIToRaster执行实质的转换工作。如果此处我们需要执行位于ArcGIS扩展模块的命令，则需要在调用命令前，先对扩展模块（如spatial analysis）进行checkout。使用以下命令：
gp.CheckOutExtension("Spatial")
gp.CellStatistics_sa(...)
否则CellStatistic_sa会抱怨说没有License可用。
代码其余部分是常规Python代码。以下取得指定目录下的txt扩展名的全部文件名，并生成对应的grid文件名。
children=os.listdir(datadir)
#only ".txt" files
txtfiles=[]
grdnames=[]
for child in children:
root,ext=os.path.splitext(child)
if (ext=='.txt'):
txtfiles.append(child);
grdname='n'+root[8:19]
grdname=grdname.replace('.','')
grdnames.append(grdname)
Python 十分有趣，如果有其它语言基础，学习起来并不难。对我们写ArcGIS脚本，有一般的Python知识就足够了。尽管效率不如传统的arcinfo AML高，但在功能上则不逊色。如果需要实现十分复杂的GIS功能，则建议使用c++, vb或.net语言操作ArcObjects。
