文件后綴為.asc
文件頭部分格式：
nrows 1232
xllcorner 484314.627259768022
yllcorner 4327241.646642711945
cellsize 9.994747084015
NODATA_value -9999
如果QGIS或者GDAL命令中導(dǎo)出來的asc文件格式是：
nrows 1232
xllcorner 484314.627259768022
yllcorner 4327241.646642711945
dx 9.994747084015
dy 9.994747084015
NODATA_value -9999
則需要手動(dòng)進(jìn)行修改一下。
如果要自動(dòng)化執(zhí)行填充泡嘴，注意文件格式的不同量愧。或者修改anuga的源代碼
anuga.asc2dem()方法辑畦，但修改的話有風(fēng)險(xiǎn)，因?yàn)椴恢肋@個(gè)方法是否也用在別處腿倚。
2)很奇怪的是第二種格式的asc纯出，也就是anuga使用的asc格式，QGIS也能打開敷燎。
所以最佳的修改方式是anuga源碼中修改的都兼容暂筝。（尚未修改20200222）
相關(guān)asc的投影prj僅僅支持下面的格式（源代碼中默認(rèn)的）
Projection UTM
Zone 55
Datum WGS84
Zunits NO
Units m
Spheroid WGS84
Xshift 500000
Yshift 10000000
Parameters
但在asc2dem中沒有看到投影轉(zhuǎn)換的地方。數(shù)據(jù)是直接從asc讀出來放到dem（一個(gè)netcdf文件懈叹，不是地理的dem文件）中的乖杠。
這個(gè)投影文件的使用地方在set_quantity('elevation')的時(shí)候分扎〕纬桑看看到底有什么用。（20200222）
沒有仔細(xì)研究畏吓。投影保持源代碼默認(rèn)的 或者改為如下墨状，
Projection NO
Zone 50
Datum WGS84
Zunits NO
Units m
Spheroid WGS84
Xshift 500000
Yshift 0
Parameters
運(yùn)行結(jié)果：sww數(shù)據(jù)投影將和輸入的原DEM數(shù)據(jù)坐標(biāo)一樣。所以懷疑投影文件沒有發(fā)揮作用菲饼。
2020-02-27日：用來計(jì)算的prj文件保持和源anuga一樣的投影肾砂，并不影響計(jì)算結(jié)果。但結(jié)果計(jì)算出來后宏悦，在QGIS顯示是镐确，依然使用源tif文件的投影包吝。jumahe案例的投影文件為：
+proj=tmerc +lat_0=0 +lon_0=116 +k=1 +x_0=500000 +y_0=0 +ellps=GRS80 +units=m +no_defs
