Arhn - архитектура программирования

R Как экспортировать мои результаты в несколько слоев растров одновременно?

В качестве первого шага я импортирую таблицу, содержащую ежедневные данные для нескольких станций отбора проб (см. Образец ниже). Впоследствии я захочу создать растровые слои для каждой даты (дня) и позже экспортировать эти растры с другими именами. Поскольку у меня есть длинный временной ряд (40 лет) с ежедневными данными, я хочу создать цикл (см. Мою программу ниже), но я не могу давать разные имена выходным растрам. Так каждый раз он давит мне предыдущий слой. Я новичок в программировании циклов, вы можете помочь мне решить эту проблему?

заранее спасибо

data=read.table(file="data.csv", sep=";", header=TRUE)
head(data)
LAMBX   LAMBY   DATE    S
600 24010   20180801    15.3
600 24010   20180802    12
600 24010   20180803    15
600 24010   20180804    14.8
600 24010   20180805    16.8
600 24010   20180806    15.1
601 24011   20180801    11
601 24011   20180802    14
601 24011   20180803    16.8
601 24011   20180804    15.1
601 24011   20180805    13.8
601 24011   20180806    12.7
602 24012   20180801    15.3
602 24012   20180802    14
602 24012   20180803    12
602 24012   20180804    16.8
602 24012   20180805    17.5
602 24012   20180806    15.1


library(dplyr)
library(sp)
library(raster)
# for each date I will do the following manipulations
for (i in data[,"DATE"]) {
# to delimit my study area
  table=filter(data, DATE == i & LAMBX %in% 601:602 & LAMBY %in% 24011:24012) 
# convert geographic coordinates
  table[,c("LAMBX", "LAMBY")]= 100*table[,c("LAMBX", "LAMBY")]

# spatialize the stations
  xy <- table[,c("LAMBX", "LAMBY")]
  sptable <- SpatialPointsDataFrame(coords = xy, data = table,
                                 proj4string = CRS("+proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs"))

# rasterize my SpatialPointsDataFrame
  rsptable <- rasterFromXYZ(as.data.frame(sptable)[, c("LAMBX", "LAMBY","S")],  crs="+proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs")
  # export my rasters layers
  writeRaster(rsptable, filename="S_date.tif", format="GTiff", overwrite=TRUE)
}
r loops dplyr lapply raster
28.11.2018

  • Я вижу T_Q в коде. Это в ваших данных? 29.11.2018

Ответы:


1

Я не вижу переменной T_Q в ваших данных, поэтому я собираюсь использовать вместо нее S. Их уловка состоит в том, чтобы создать пустой объект списка, который заполняется на каждой итерации цикла. Затем вы можете создать стек результатов. Каждый слой стопки - это отдельная дата.

library(dplyr)
library(sp)
library(raster)

data <- read.table(text = "
LAMBX   LAMBY   DATE    S
600 24010   20180801    15.3
600 24010   20180802    12
600 24010   20180803    15
600 24010   20180804    14.8
600 24010   20180805    16.8
600 24010   20180806    15.1
601 24011   20180801    11
601 24011   20180802    14
601 24011   20180803    16.8
601 24011   20180804    15.1
601 24011   20180805    13.8
601 24011   20180806    12.7
602 24012   20180801    15.3
602 24012   20180802    14
602 24012   20180803    12
602 24012   20180804    16.8
602 24012   20180805    17.5
602 24012   20180806    15.1", header = TRUE)

crs_string <- "+proj=lcc +lat_1=46.8 +lat_0=46.8 +lon_0=0 +k_0=0.99987742 +x_0=600000 +y_0=2200000 +a=6378249.2 +b=6356515 +towgs84=-168,-60,320,0,0,0,0 +pm=paris +units=m +no_defs"

res <- list()
for (i in seq_along(unique(data[,"DATE"]))){
  table <- filter(data,
                  DATE == unique(data[,"DATE"])[i] &
                  LAMBX %in% 601:602 &
                  LAMBY %in% 24011:24012)
  table[,c("LAMBX", "LAMBY")] <-  100 * table[,c("LAMBX", "LAMBY")]

  xy       <- table[,c("LAMBX", "LAMBY")]
  sptable  <- SpatialPointsDataFrame(coords = xy, data = table,
                                    proj4string = CRS(crs_string))
  rsptable <- rasterFromXYZ(
    as.data.frame(sptable)[, c("LAMBX", "LAMBY","S")],  crs = crs_string)
  res[[i]] <- rsptable
}

res <- stack(res)
writeRaster(res, filename="S_date.tif", format="GTiff", overwrite=TRUE)
29.11.2018
  • Большое спасибо jsta. Программа работает отлично 30.11.2018
    • Новые материалы
    Коллекции публикаций по глубокому обучению
    Последние пару месяцев я создавал коллекции последних академических публикаций по различным подполям глубокого обучения в моем блоге https://amundtveit.com - эта публикация дает обзор 25..
    Представляем: Pepita
    Фреймворк JavaScript с открытым исходным кодом Я знаю, что недостатка в фреймворках JavaScript нет. Но я просто не мог остановиться. Я хотел написать что-то сам, со своими собственными..
    Советы по коду Laravel #2
    1-) Найти // You can specify the columns you need // in when you use the find method on a model User::find(‘id’, [‘email’,’name’]); // You can increment or decrement // a field in..
    Работа с временными рядами спутниковых изображений, часть 3 (аналитика данных)
    Анализ временных рядов спутниковых изображений для данных наблюдений за большой Землей (arXiv) Автор: Рольф Симоэс , Жильберто Камара , Жильберто Кейрос , Фелипе Соуза , Педро Р. Андраде ,..
    3 способа решить квадратное уравнение (3-й мой любимый) -
    1. Методом факторизации — 2. Используя квадратичную формулу — 3. Заполнив квадрат — Давайте поймем это, решив это простое уравнение: Мы пытаемся сделать LHS,..
    Создание VR-миров с A-Frame
    Виртуальная реальность (и дополненная реальность) стали главными модными терминами в образовательных технологиях. С недорогими VR-гарнитурами, такими как Google Cardboard , и использованием..
    Демистификация рекурсии
    КОДЕКС Демистификация рекурсии Упрощенная концепция ошеломляющей О чем весь этот шум? Рекурсия, кажется, единственная тема, от которой у каждого начинающего студента-информатика..