plotAnglesDif: Plot differences of two SAR records

In ragram: Radargrammetry

Description

Plots the difference in angles of two provided by GCPs/TPs as points.

Usage

plotAnglesDif(object, slave, z = "thetaIn", xlab = "Longitude",
  ylab = "Latitude", breaks = 255,
  col.regions = colorRampPalette(colors)(breaks), cex = 1, cex.axis = 1,
  grid = T, variogram.fit = T, plot.fit = F, plot.legend = T,
  pch = 20, legend.lab = expression(paste(Delta, theta[i], " [deg]")),
  interpolate = F, aggregate.fact = 100, disparity = F, h = 100, ...)

## S4 method for signature 'SAR,ANY'
plotAnglesDif(object, slave, z = "thetaIn",
  xlab = "Longitude", ylab = "Latitude", breaks = 255,
  col.regions = colorRampPalette(colors)(breaks), cex = 1, cex.axis = 1,
  grid = T, variogram.fit = T, plot.fit = F, plot.legend = T,
  pch = 20, legend.lab = expression(paste(Delta, theta[i], " [deg]")),
  interpolate = F, aggregate.fact = 100, disparity = F, h = 100, ...)

## S4 method for signature 'SAR,SAR'
plotAnglesDif(object, slave, z, xlab, ylab, breaks,
  col.regions, cex, cex.axis, grid, variogram.fit, plot.fit, plot.legend, pch,
  legend.lab, interpolate, aggregate.fact, disparity, h,
  xlim = c(min(object@extent@xmin, slave@extent@xmin), max(object@extent@xmax,
  slave@extent@xmax)), ylim = c(min(object@extent@ymin, slave@extent@ymin),
  max(object@extent@ymax, slave@extent@ymax)), asp = 1/cos((mean(range(ylim))
  * pi)/180), zlim = c(min(angleDif@data[, 1]), max(angleDif@data[, 1])), ...)

## S4 method for signature 'SARSet,ANY'
plotAnglesDif(object, slave, z = "thetaIn",
  xlab = "Longitude", ylab = "Latitude", breaks = 255,
  col.regions = colorRampPalette(colors)(breaks), cex = 1, cex.axis = 1,
  grid = T, variogram.fit = T, plot.fit = F, plot.legend = T,
  pch = 20, legend.lab = expression(paste(Delta, theta[i], " [deg]")),
  interpolate = F, aggregate.fact = 100, disparity = F, h = 100, ...)

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'.
xlab	Character. Title of the x axis.
ylab	Character. Title of the y axis.
breaks	Integer. Number of breaks of the color ramp.
col.regions	Color ramp.
cex	Numeric. Size of label text.
cex.axis	Numeric. Size of axes.
grid	Logical. Plot a grid?
variogram.fit	Logical. Fit a Gaussian variogram?
plot.fit	Logical. Plot the fitted variogram?
plot.legend	Logical. Plot the legend?
pch	Integer. Point type.
legend.lab	Character. Title of the legend.
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. If no aggregation is done, the interpolation might take a very long time...
disparity	Logical. Should the expected disparity be plotted?
h	Numeric. Relative height for plotting the expected disparity.
...	graphical parameters. Any argument that can be passed to plot, such as axes=FALSE and main='title'.
xlim	Numeric vector of length 2, giving the x coordinate range.
ylim	Numeric vector of length 2, giving the y coordinate range.
asp	Numeric vector. Aspect, default is 1/cos((mean(range(ylim)) * pi)/180).
zlim	Numeric vector of length 2, giving the z coordinate range. Range of values that are plotted.

See Also

anglesDif, angles, plotAngles, GeolocationPoints-class

Examples

data(kili)

plotAnglesDif(kili[[1]], kili[[3]])

plotAnglesDif(kili)

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

plotAnglesDif(kili[[1]], kili[[5]], disparity = TRUE,
legend.lab = 'Disparity [m]', h=100)

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