Displaying 1 result from an estimated 1 matches for "your_lon_low_limit".
Did you mean:
your_lat_low_limit
2010 Jul 27
0
3d topographic map [SEC=UNCLASSIFIED]
...- attr(ETOPO5,'lon') # all lon
yf <- attr(ETOPO5,'lat') # all lat
zf <- ETOPO5 # elevation
(3) The topography data has to be flipped horizontally and vertically and you also require the data to be placed on a regular grid and you need to specify limits,
xn <- which(xf>=your_lon_low_limit & xf<=your_lon_high_limit)
yn <- which(yf>=your_lat_low_limit & yf<=your_lat_high_limit )
x <- xf[xn]
y <- yf[yn]
z <- zf[xn,yn]
za <- apply(z,1,rev)
zb <- t(za)
yb <- -rev(y)
xb <- x
nx <- length(xb)
dx <- (xb[nx]-xb[1])/(nx-1)
ny <- length(yb)...