anglesDif: Angles differences of two SAR records

In ragram: Radargrammetry

Description

Returns the differences in angles within the intersection area of two SAR records. Uses the angles provided by GCPs/TPs of both records. The angles for the slave image are interpolated using Kriging (krige). Returns the results as a SpatialPointsDataFrame object.

Usage

anglesDif(object, slave = NULL, z = "thetaIn", variogram.fit = T,
  plot.fit = F, interpolate = F, aggregate.fact = 1000, ...)

## S4 method for signature 'SAR,ANY'
anglesDif(object, slave = NULL, z = "thetaIn",
  variogram.fit = T, plot.fit = F, interpolate = F,
  aggregate.fact = 1000, ...)

## S4 method for signature 'SAR,SAR'
anglesDif(object, slave = NULL, z = "thetaIn",
  variogram.fit = T, plot.fit = F, interpolate = F,
  aggregate.fact = 1000, ...)

## S4 method for signature 'SARSet,ANY'
anglesDif(object, slave = NULL, z = "thetaIn",
  variogram.fit = T, plot.fit = F, interpolate = F,
  aggregate.fact = 1000, ...)

Arguments

object	Object of the SAR-class or a subclass (e.g. Sentinel-class or TSX-class).
slave	Object of the SAR-class or a subclass (e.g. Sentinel-class or TSX-class).
z	Character. Either incidence ('thetaIn', default) or elevation angles 'thetaEl'.
variogram.fit	Logical. Fit a Gaussian variogram?
plot.fit	Logical. Plot the fitted variogram?
interpolate	Logical. Interpolate angles?
aggregate.fact	Integer. Aggregation factor expressed as number of cells in each direction (horizontally and vertically). Or two integers (horizontal and vertical aggregation factor) or three integers (when also aggregating over layers). See aggregate method of raster package.
...	Further arguments (ignored).

Value

SpatialPointsDataFrame

See Also

angles, plotAnglesDif, GeolocationPoints-class

Examples

data(kili)

anglesDif(kili[[1]], kili[[5]])
plotAnglesDif(kili[[1]], kili[[5]])

anglesDif(master, slave, interpolate=TRUE, aggregate.fact=100)
plotAnglesDif(master, slave, interpolate=TRUE, aggregate.fact=100)

anglesDif.sp <- anglesDif(master, slave, interpolate = TRUE, aggregate.fact=10)
anglesDif.raster <- rasterFromXYZ(anglesDif.sp)
plot(anglesDif.raster)

ragram documentation built on May 2, 2019, 4:42 p.m.