Флаг "-c" используется только для типов данных Byte и UInt16, почитайте уже наконец справку по модулю.list писал(а): А вот после флага -c , все в порядке. Какой тогда применить флаг? -c?
Файл должен экспортироваться несмотря на ошибку.
Геоинформационные системы (ГИС) и Дистанционное зондирование Земли
Флаг "-c" используется только для типов данных Byte и UInt16, почитайте уже наконец справку по модулю.list писал(а): А вот после флага -c , все в порядке. Какой тогда применить флаг? -c?
На здоровьеlist писал(а):Спасибо
Т.е. считается всегда для какого-нибудь одного дня года.day=integer [required]
No. of day of the year (1-365)
Options: 1-365
Если я правильно понял то r.sun мне не сможет вычислить целый год в одном растре (четырех =beam, diffuse, global, insol_time). Но может вычислить их 365 х beam, diffuse, global, insol_time.Александр Мурый писал(а): Если вам нужно посчитать за все дни года, надо писать скрипт, как указано здесь.
Получается приказ, чтобы вычислить солярный потенциал для всех дней в году должен выглядеть так:Александр Мурый писал(а): Если вам нужно посчитать за все дни года, надо писать скрипт, как указано здесь.
Код: Выделить всё
r.sun -s elevin=elevation.dem day=$DAY lin=$LINKE step=0.05 \
beam_rad=rad_beam.$DAY_STR diff_rad=rad_diffuse.$DAY_STR \
refl_rad=rad_reflected.$DAY_STR glob_rad=rad_global.$DAY_STR \
insol_time=rad_insol_time.$DAY_STR
Код: Выделить всё
for DAY in `seq 1 365` ; do
DAY_STR=`echo $DAY | awk '{printf("%.03d", $1)}'`
echo "Processing day $DAY_STR at `date` ..."
Код: Выделить всё
DAY_STR=`echo $DAY | awk '{printf("%.03d", $1)}'`
Ну а как тогда написать скрипт с годовым вычислением. Меня конкретно интересует day=?. Все остальное знаю а day нет. Может так:KolesovDmitry писал(а):Например, если $DAY равен 45, то DAY_STR равен 045.
Код: Выделить всё
r.sun -s elevin=elevation.dem day=1-365 lin=$LINKE step=0.05 \
beam_rad=rad_beam.365 diff_rad=rad_diffuse.365 \
refl_rad=rad_reflected.365 glob_rad=rad_global.365 \
insol_time=rad_insol_time.365
Код: Выделить всё
r.sun.daily elevation=dtm@PERMANENT aspect=aspect slope=slope lin=4.2 alb=0.0 coefbh=beam@PERMANENT coefdh=dif@PERMANENT beam_rad=pri_r_smaz_year diff_rad=dif_r_smaz_year refl_rad=refl_r_smaz_year glob_rad=glo_r_smaz_year day=1-365 step=0.5
linke_value=float
A single value of the Linke atmospheric turbidity coefficient [-]
Options: 0.0-7.0
Default: 3.0
Хорошо, спасибо, буду пытаться. А что касается Linke atmospheric turbidity coefficient, так он у меня вроде бы указан lin=4.2 (это годовой)Александр Мурый писал(а):
Пример:Код: Выделить всё
r.sun.daily elevation=dtm@PERMANENT aspect=aspect slope=slope lin=4.2 alb=0.0 coefbh=beam@PERMANENT coefdh=dif@PERMANENT beam_rad=pri_r_smaz_year diff_rad=dif_r_smaz_year refl_rad=refl_r_smaz_year glob_rad=glo_r_smaz_year day=1-365 step=0.5
- Обратите внимание на параметр "linke_value" - важный коэффифиент, который у вас нигде не указан. Читайте в мануале про него.
linke_value=float
A single value of the Linke atmospheric turbidity coefficient [-]
Options: 0.0-7.0
Default: 3.0
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 8 гостей
© GIS-Lab и авторы, 2002-2017. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов (подробнее).