Это команда r.seriesgrass, которую можно запустить в провайдере бесплатного хостинга OnWorks, используя одну из наших многочисленных бесплатных онлайн-рабочих станций, таких как Ubuntu Online, Fedora Online, онлайн-эмулятор Windows или онлайн-эмулятор MAC OS.
ПРОГРАММА:
ИМЯ
р.серия - Делает каждое значение выходной ячейки функцией значений, присвоенных
соответствующие ячейки во входных слоях растровой карты.
КЛЮЧЕВЫЕ СЛОВА
растр, агрегация, серия
СИНТАКСИС
р.серия
р.серия --Помогите
р.серия [-nz] [вход=имя[,имя, ...]] [файл=имя] выходной=имя[,имя, ...]
метод=string[,string, ...] [квантиль=плавать[,плавать, ...]] [весами=плавать[,плавать, ...]]
[ассортимент=вот, привет] [-затирать] [-помощь] [-подробный] [-тихий] [-ui]
Флаги:
-n
Распространять NULL
-z
Не оставляйте файлы открытыми
- перезаписать
Разрешить выходным файлам перезаписывать существующие файлы
--Помогите
Распечатать сводку использования
--подробный
Подробный вывод модуля
--тихий
Тихий выход модуля
--уи
Принудительный запуск диалогового окна GUI
Параметры:
вход=имя [, имя, ...]
Имя входной растровой карты (ов)
файл=имя
Входной файл с одним именем растровой карты и, необязательно, одним весом на строку, разделителем полей
между именем и весом |
выходной=имя [, имя, ...] [требуется]
Имя выходной растровой карты
метод=строка [, строка, ...] [требуется]
Агрегатная операция
Опции: средний, рассчитывать, медиана, Режим, минимум, мин_растр, максимум, макс_растр,
стандартное отклонение, диапазона, сумма дисперсия разнообразие, наклон компенсировать, Деткофф значение, кварт1,
кварт3, проц90, квантиль асимметрия эксцесс
квантиль=float [, float, ...]
Квантиль для расчета для метода = квантиль
Опции: 0.0-1.0
весами=float [, float, ...]
Весовой коэффициент для каждой входной карты, значение по умолчанию 1.0 для каждой входной карты.
ассортимент=вот, привет
Игнорировать значения вне этого диапазона
ОПИСАНИЕ
р.серия делает каждое значение выходной ячейки функцией значений, присвоенных
соответствующие ячейки во входных слоях растровой карты.
Доступны следующие методы:
· Средний: среднее значение
· Count: количество ненулевых ячеек
· Median: среднее значение
· Режим: наиболее часто встречающееся значение
· Минимум: наименьшее значение
· Максимум: максимальное значение
· Диапазон: диапазон значений (макс - мин)
· Stddev: стандартное отклонение
· Сумма: сумма значений
· Дисперсия: статистическая дисперсия
· Разнообразие: количество разных ценностей
· Наклон: наклон линейной регрессии
· Смещение: смещение линейной регрессии
· Detcoeff: коэффициент детерминации линейной регрессии
· Tvalue: значение t линейной регрессии
· Min_raster: номер растровой карты с минимальным значением временного ряда
· Max_raster: номер растровой карты с максимальным значением временного ряда
Обратите внимание, что большинство параметров принимают несколько ответов, что позволяет использовать несколько агрегатов.
вычисляется за один прогон, например:
r.series input = map1, ..., mapN \
output = map.mean, map.stddev \
метод = среднее, стандартное отклонение
или:
r.series input = map1, ..., mapN \
output = map.p10, map.p50, map.p90 \
метод = квантиль, квантиль, квантиль \
квантиль = 0.1,0.5,0.9
Для всех опций должно быть указано одинаковое количество значений.
ПРИМЕЧАНИЯ
Нет данных (НОЛЬ) обращение
В -n флаг, любая ячейка, для которой любая из соответствующих входных ячеек имеет значение NULL, является
автоматически устанавливается в NULL (распространение NULL). Агрегатная функция не вызывается, поэтому
все методы ведут себя так по отношению к -n флаг.
Без -n флаг, полный список входов для каждой ячейки (включая NULL) передается в
агрегатная функция. Отдельные агрегаты могут обрабатывать данные по своему усмотрению. В основном они
просто вычислить агрегат по значениям, отличным от NULL, и получить результат NULL, только если все
входы NULL.
Минимальный и максимальный анализ
Команда мин_растр и max_raster методы генерируют карту с номером растровой карты
который содержит минимальное / максимальное значение временного ряда. Нумерация начинается с 0 вплоть до n
для первого и последнего растра, указанного в input =, Соответственно.
Диапазон анализ
Если же линия индикатора диапазон = задана опция, любые значения, выходящие за пределы этого диапазона, будут рассматриваться как
если они были NULL. В ассортимент параметр может быть установлен на низко высоко пороги: значения за пределами
этого диапазона обрабатываются как NULL (т. е. они будут игнорироваться большинством агрегатов или будут
приведет к тому, что результат будет ПУСТО (NULL), если задано -n). В низко высоко пороги с плавающей точкой,
так что используйте -инф or инф для одного порога (например, диапазон = 0, инф игнорировать отрицательные значения или
диапазон = -inf, -200.4 игнорировать значения выше -200.4).
Линейные приводы регресс
Линейная регрессия (наклон, смещение, коэффициент детерминации, t-значение) предполагает равенство
временные интервалы. Если данные имеют нерегулярные временные интервалы, можно использовать NULL растровые карты.
вставляется во временной ряд, чтобы сделать временные интервалы равными (см. пример).
Квантили
р.серия может рассчитывать произвольные квантили.
Память потребление
Использование памяти не является проблемой, так как р.серия нужно только удерживать одну строку от каждой карты в
времени.
Руководство of открытый файл рамки
Количество обрабатываемых растровых карт задается лимитом операционной системы. Для
Например, как жесткие, так и мягкие пределы обычно равны 1024. Мягкое ограничение можно изменить.
например, ulimit -n 1500 (операционные системы на базе UNIX), но не выше жесткого
предел. Если он слишком низкий, вы можете как суперпользователь добавить запись в
/etc/security/limits.conf
#
your_username жесткий nofile 1500
Это повысит жесткое ограничение до 1500 файлов. Имейте в виду, что для открытия большего количества файлов требуется больше
БАРАН. См. Также Подсказки на странице Wiki для обработки больших растровых данных.
Для каждой карты можно указать весовой коэффициент с помощью весами вариант. Использование весов
может иметь значение при вычислении суммы или среднего числа карт с разной временной протяженностью.
Вес по умолчанию - 1.0. Количество гирь должно совпадать с количеством
входные карты и должны иметь такой же порядок. Во входном файле также можно указать вес.
Использовать файл возможность анализировать большое количество растровых карт, не открывая файлы
limit и предельный размер аргументов командной строки. Вычисление медленнее, чем
вход опционный метод. Для каждой отдельной строки в выходной карте (ах) открываются все входные карты.
и закрыто. Объем ОЗУ будет увеличиваться линейно с количеством указанных входных карт.
Параметры ввода и файла исключают друг друга. Ввод - текстовый файл с новой строкой
разделенный список имен растровых карт и дополнительных весов. Как разделитель между названием карты
и вес символа "|" должны быть использованы.
ПРИМЕРЫ
. р.серия с подстановочными знаками:
r.series input = "` g.list pattern = 'insitu_data. *' sep =, `" \
output = insitu_data.stddev метод = stddev
Обратите внимание г.список скрипт также поддерживает регулярные выражения для выбора названий карт.
. р.серия с NULL растровыми картами (чтобы рассматривать "полный" временной ряд):
r.mapcalc "dummy = null ()"
r.series in=map2001,map2002,dummy,dummy,map2005,map2006,dummy,map2008 \
out = res_slope, res_offset, res_coeff meth = slope, offset, detcoeff
Пример вычисления нескольких агрегатов за один прогон (3 результирующих агрегата из двух
входные карты):
r.series in = one, two out = result_avg, res_slope, result_count meth = sum, slope, count
Пример использования файловой опции r.series:
cat> input.txt << EOF
map1
map2
map3
EOF
r.series file = input.txt out = result_sum meth = sum
Пример использования опции файла для r.series, включая веса. Вес 0.75 должен быть
присвоено map2. Поскольку у других карт нет веса, мы можем не указывать его:
cat> input.txt << EOF
map1
map2 | 0.75
map3
EOF
r.series file = input.txt out = result_sum meth = sum
Пример подсчета количества дней выше определенной температуры с использованием среднесуточного значения
карты ('???' как подстановочный знак DOY):
# Подход для систем на основе оболочки
r.series input = `g.list rast pattern =" temp_2003 _ ??? _ avg "sep =,` \
output = temp_2003_days_over_25deg range = 25.0,100.0, XNUMX method = count
# Подход в два этапа (например, для систем Windows)
g.list rast pattern = "temp_2003 _ ??? _ avg" output = mapnames.txt
r.series file = mapnames.txt \
output = temp_2003_days_over_25deg range = 25.0,100.0, XNUMX method = count
Используйте r.seriesgrass в Интернете с помощью сервисов onworks.net