What is the best way to make directional h-scatterplots (lagged scatterplots)? There is hscat() in gstat, but I don't see a way to limit point pair selection to certain directions. I have a 3D dataset based on logs from 30 different wells, with vertical spacing between points and total extent both smaller than horizontal spacing and extent. Here is what I have figured out so far. For the horizontal direction, I can select points in a given x-y plane before calling hscat(). For the vertical direction, I can do a de facto selection by limiting the maximum lag distance to something smaller than the minimum distance between my wells. My main interest for this is in 2 directions, horizontal and vertical, but it would be nice to specify angles as in the anis argument for vgm(). Thanks, Scott Waichler Pacific Northwest National Laboratory scott.waichler at pnl.gov
How to make directional h-scatterplots
3 messages · Waichler, Scott R, Edzer Pebesma
you can pass the alpha, beta, tol.hor and tol.ver parameters to hscat, who will pass it on to variogram(). The problem is that the alpha parameter, when set, is also passed to xyplot() and messes up the plot there - I will choose a different name for that one.
On 09/10/2010 01:59 AM, Waichler, Scott R wrote:
What is the best way to make directional h-scatterplots (lagged scatterplots)? There is hscat() in gstat, but I don't see a way to limit point pair selection to certain directions. I have a 3D dataset based on logs from 30 different wells, with vertical spacing between points and total extent both smaller than horizontal spacing and extent. Here is what I have figured out so far. For the horizontal direction, I can select points in a given x-y plane before calling hscat(). For the vertical direction, I can do a de facto selection by limiting the maximum lag distance to something smaller than the minimum distance between my wells. My main interest for this is in 2 directions, horizontal and vertical, but it would be nice to specify angles as in the anis argument for vgm(). Thanks, Scott Waichler Pacific Northwest National Laboratory scott.waichler at pnl.gov
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebesma at wwu.de
In response to this response, I uploaded gstat_0.9-70 to CRAN.
On 09/10/2010 07:52 AM, Edzer Pebesma wrote:
you can pass the alpha, beta, tol.hor and tol.ver parameters to hscat, who will pass it on to variogram(). The problem is that the alpha parameter, when set, is also passed to xyplot() and messes up the plot there - I will choose a different name for that one. On 09/10/2010 01:59 AM, Waichler, Scott R wrote:
What is the best way to make directional h-scatterplots (lagged scatterplots)? There is hscat() in gstat, but I don't see a way to limit point pair selection to certain directions. I have a 3D dataset based on logs from 30 different wells, with vertical spacing between points and total extent both smaller than horizontal spacing and extent. Here is what I have figured out so far. For the horizontal direction, I can select points in a given x-y plane before calling hscat(). For the vertical direction, I can do a de facto selection by limiting the maximum lag distance to something smaller than the minimum distance between my wells. My main interest for this is in 2 directions, horizontal and vertical, but it would be nice to specify angles as in the anis argument for vgm(). Thanks, Scott Waichler Pacific Northwest National Laboratory scott.waichler at pnl.gov
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Edzer Pebesma Institute for Geoinformatics (ifgi), University of M?nster Weseler Stra?e 253, 48151 M?nster, Germany. Phone: +49 251 8333081, Fax: +49 251 8339763 http://ifgi.uni-muenster.de http://www.52north.org/geostatistics e.pebesma at wwu.de