Autumn Mannsfeld's questions

Criei um SpatRaster com o pacote terra em R e tenho dois agrupamentos de pontos dentro do SpatRaster maior, que então separo em grupos para observar os agrupamentos com mais detalhes. Uso terra::plot para plotar o SpatRaster maior e os grupos individuais, mas quero manter a escala de gradiente igual para todos os gráficos. O problema é que o intervalo de valores para cada um dos gráficos é diferente, e não consegui descobrir como definir os limites inferior e superior para o gradiente de cor de cada terra::plot para que sejam iguais, de modo que as cores sejam consistentes.

O pacote ggplot2 não funcionou para esses dados porque são dados SpatRaster, portanto scale_color_gradientnão funciona, tentei uma combinação do colargumento em terra::plot, como

plot(DTM_50570, main = "DTM", col = data.frame(from = 0, to = 2000, (terra.pal = terrain.colors(255))))

plot(DTM_50570, main = "DTM", col = data.frame(from = c(0, 1000), to = c(1000, 2000), col = terrain.colors(255)))

plot(DTM_50570, main = "DTM", col = data.frame(from = c(0:10:2550), to = c(0:10:2550), col = terrain.colors(255)))

sem sucesso.

Os dados que estou usando estão neste link: https://tinitaly.pi.ingv.it/Download_Area1_1.html , seção W50570.

Aqui está o código que uso para plotar e separar os grupos:

library(terra)

DTM_50570_32632 <- rast(paste0("Downloads/w50570_s10/w50570_s10/w50570_s10.tif"))
DTM_50570 <- project(x = DTM_50570_32632, "EPSG:4326")

DTM_longitude <- c(11.80111, 11.80222, 11.80333, 11.95111, 11.95222, 11.95333)
DTM_latitude <-c(45.98333, 45.98777, 45.98555, 45.94111, 45.94888, 45.94333)

group_1_extent <- ext(11.80, 11.81, 45.98, 45.99)
DTM_group_1 <- crop(DTM_50570, group_1_extent)

group_2_extent <- ext(11.95, 11.96, 45.94, 45.95)
DTM_group_2 <- crop(DTM_50570, group_2_extent)

DTM_longitude_group_1 <- c(11.80111, 11.80222, 11.80333)
DTM_latitude_group_1 <-c(45.98333, 45.98777, 45.98555)

DTM_longitude_group_2 <- c(11.95111, 11.95222, 11.95333)
DTM_latitude_group_2 <-c(45.94111, 45.94888, 45.94333)

par(mfrow = c(1, 3))
plot(DTM_50570, main = "DTM")
points(DTM_longitude, DTM_latitude, pch = 20)

plot(DTM_group_1, main = "DTM group 1")
points(DTM_longitude_group_1, DTM_latitude_group_1, pch = 20)

plot(DTM_group_2, main = "DTM group 2")
points(DTM_longitude_group_2, DTM_latitude_group_2, pch = 20)

E os gráficos resultantes (todos com escalas de gradiente de cores diferentes)

Tenho um conjunto de dados brutos que usa um decimal como separador de tempo e quero que ele tenha dois pontos, mas não tenho ideia de como alterá-lo. Meu tempo bruto (com decimais) está na coluna "df$Time", e meu tempo correto (com os dois pontos) estará em "df$time" (para manter a consistência com os outros dataframes no meu código). df$Date é uma coluna de caracteres, df$Time é uma coluna numérica. Estou usando uma infinidade de bibliotecas, algumas incluindo dplyre stringr.

Não sei bem como criar uma tabela bonita com dados recriáveis. Escrever nesta plataforma é novidade para mim. Mas abaixo, incluo um trecho dos dados de tempo e alguns métodos que tentei. Qualquer ajuda é bem-vinda! Parece tão simples, mas realmente não consigo entender.

(Também peço desculpas pelo formato dos meus dados abaixo, tentei fazer uma tabela, mas uma que parecesse melhor e, novamente, nada funcionou. Acho que preciso dar uma pausa na solução de problemas, rs)

Aqui está o resultado após a execução dput(head(df)):

structure(list(Date = c("01/03/2024", "01/03/2024", "01/03/2024", "01/03/2024", "01/03/2024", "01/03/2024"), Time = c(0, 0.05, 0.1, 0.15, 0.2, 0.25)), row.names = c(NA, 6L), class = "data.frame")

Aqui está um trecho de duas horas de dados simplificados (incluindo apenas colunas relevantes) que mostra o problema:

          Date  Time
181 01/03/2024 15.00  
182 01/03/2024 15.05
183 01/03/2024 15.10
184 01/03/2024 15.15
185 01/03/2024 15.20
186 01/03/2024 15.25
187 01/03/2024 15.30
188 01/03/2024 15.35
189 01/03/2024 15.40
190 01/03/2024 15.45
191 01/03/2024 15.50
192 01/03/2024 15.55
193 01/03/2024 16.00
194 01/03/2024 16.05
195 01/03/2024 16.10
196 01/03/2024 16.15
197 01/03/2024 16.20
198 01/03/2024 16.25
199 01/03/2024 16.30
200 01/03/2024 16.35
201 01/03/2024 16.40
202 01/03/2024 16.45
203 01/03/2024 16.50
204 01/03/2024 16.55
205 01/03/2024 17.00

Aqui estão alguns dos métodos que tentei:

Método 1:

df$time <- df %>% mutate(across(everything(), ~str_replace(., ",", ".")))

Método 2:

df$time <- str_replace(df$Time, ".", ":")

Método 3:

df$Time <- as.character(df$Time)

df$time <- gsub(".", ":", df$Time, fixed = T)

Método 4:

df$time <- strptime(as.character(df$Time), "%H.%M")

df$time <- format(df$time, "%H:%M")

Método 5:

df$time <- as.POSIXlt(df$Time)

Método 6:

df$time <- df$Time df$time <- gsub(".", ":", df$time)

Método 7:

df$time <- df$Time df$time <- str_replace_all(df$time, ".", ":")