Recorte de imágenes Landsat mediante Extent en R

Question

Estoy siguiendo un tutorial sobre Teledetección dentro de R y estoy intentando recortar una imagen Landsat, sin embargo estoy obteniendo el error: Error in validityMethod(object) : invalid extent: ymin >= ymax

Mi código es:

# Find the current extent
extent(landsat)
class       : Extent 
xmin        : 269625 
xmax        : 337155 
ymin        : 3221685 
ymax        : 3289215 

# Create new extent object
e <- extent(250000, 325000, -3150000,-3250000)

# Crop Landsat
l_crop <- crop(landsat, extent(e))

¿Hay algún error en la sintaxis o me he perdido algo?

Answer 1

No estoy muy familiarizado con R pero parece que has confundido los valores ymin y ymax al crear tu objeto extent:

# Arguments are xmin,xmax,ymin,ymax
# Note: ymin > ymax which is invalid.
e <- extent(250000, 325000, -3150000, -3250000) 

Si lo cambias por:

e <- extent(250000, 325000, -3250000, -3150000)

Esto debería crear un objeto de extensión válido.