需要wget脚本软件具体见GNSS观测数据o文件批量下载(Linux Shell)
具体实现代码如下:
#!/bin/bash
# function: Download the navigation file
# Usage : wgetbrd YYYY DOY SITE outDir
###########################################################################################
if [ $# != 3 ]; then #判断输入格式是否正确
echo "Usage: wgetbrd YYYY DOY outDir"
echo "Usage: wgetbrd 2019 001 /date/gnss/"
exit
fi
yyyy=$1
day=$2
ObsPath=$3
doy=`echo ${day} | awk '{printf("%03d",$0)}'`; #年积日
yy=`echo ${yyyy} | awk '{printf("%2.2d",$1-int($1/100)*100)}'`; #年份
#Check the file existence #判断文件是否存在
if [ -f ${ObsPath}brdm${doy}0.${yy}p ] ; then
echo "Obsfile : ${ObsPath}brdm${doy}0.${yy}p exist"
exit
fi
#数据下载地址
URL_EPNCB="ftp://epncb.oma.be/pub/obs/BRDC/${yyyy}/" #brdm文件下载地址
RINEX3L="BRDC*${yyyy}${doy}0000*MN.rnx.gz" #Rinex3 LONG NAME 用于匹配BRDC00GOP_R_20200010000_01D_MN.rnx.gz
wget -t 10 -T 120 -q ${URL_EPNCB}${RINEX3L} -P ${ObsPath}
# Rename the 3.X observations
if [ -f ${ObsPath}${RINEX3L} ]; then
gzip -d ${ObsPath}BRDC*${yyyy}${doy}*.rnx.gz
file=`ls ${ObsPath}BRDC*${yyyy}${doy}*.rnx`
mv ${file} ${ObsPath}brdm${doy}0.${yy}p #将文件改名为.p
fi
# If the 3.x file donoe exit, download the 2.x file
if [ -f ${ObsPath}brdm${doy}0.${yy}p ] ; then
echo "Download completed: brdm${doy}0.${yy}p"
else
echo "***WARNING***:Download failed: brdm${doy}0.${yy}p"
exit
fi
运行命令为:
wgetbrd 2021 001 /home/zh/gnss/
此处表示下载2021年001日的广播星历,输出路径为/home/zh/gnss/。
若成功下载会有如下提示:
Download completed: brdm0010.21p
如有问题地方欢迎咨询与交流。