affineCoords: Affine Transformation (rotation and shift) on Spatial...

In hyperbrick: Accessory Tools for Preprocessing Hyper-Spectral Images

Affine Transformation (rotation and shift) on Spatial Coordinates

Description

Affine transformations are of type f(x) = Ax + b, where x is the spatial coordinates (2D in this case), A is a rotation matrix (it can also include scale/shear parameters, but only rotation is considered here), and b is the translation (xy shift) parameters.

Usage

affineCoords(s, angle = 0, xy_shift = c(0, 0))

Arguments

s	A spatial object where raster::coordinates() can be extracted from. Example classes: Extent, RasterLayer, RasterBrick, RasterStack (from package raster), SpatialPolygon (from package sp).
angle	A numeric value of the angle (in degrees, 0-360) to rotate s. A negative value will change the direction of rotation to clockwise.
xy_shift	A numeric vector of length two with the x and y shift (the translation parameters).

Value

A two-column matrix with the transformed coordinates (xy).

See Also

affineBrick()

Examples

p <- system.file('exdata', 'soybean.tif', package = 'hyperbrick')
im <- brick(p)
print(im)

# view band-3
plot(im[[3]], col = gray.colors(20), asp = 0)

# draw a spatial polygon on image
pol <- Polygon(extent(c(40, 85, 50, 150)))
lines(pol)

# rotate and shift the spatial polygon
new_pol <- affineCoords(pol, angle = -3, xy_shift = c(-11, 0))
plot(im[[3]], col = gray.colors(20), asp = 0)
lines(new_pol)

# do some analysis within it, like:
new_pols <- SpatialPolygons(list(Polygons(list(Polygon(new_pol)), "id0")))
plot(mask(im[[3]], new_pols))
mean(extract(im[[3]], new_pols)[[1]])

hyperbrick documentation built on April 1, 2022, 9:07 a.m.