- 积分
- 11
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2023-5-31
- 最后登录
- 1970-1-1

|
发表于 2023-6-1 12:40:26
|
显示全部楼层
#!/bin/bash
j=1
# 维度变化根据实际情况修改
for((lat=26; lat<=41; lat++))
do
lat1=$[($j-1)*3601+1]
lat2=$[$j*3601]
aster_name="ASTGTMV003_N"
aster_name=$aster_name$lat"*"
i=1
for filename in $aster_name
do
#echo $filename
lon1=$[($i-1)*3601+1]
lon2=$[$i*3601]
new=$(printf "%06d-" ${lon1})
new=$new$(printf "%06d." ${lon2})
new=$new$(printf "%06d-" ${lat1})
new=$new$(printf "%06d" ${lat2})
#echo $new
gdal_translate -of ENVI -co INTERLEAVE=BSQ ${filename} ${new}
let i=i+1
done
let j=j+1
done
改了一下脚本,解决多行多列,保存成bash脚本即可 |
|