Hi R Users, Is it possible to visualize river flow in 3D (latitude, longitude with respect to depth)? The example of my data looks like. Any suggestions?> dat1long lat depth flow 1 1015.9 857 1.00 1.50 2 1015.9 857 1.25 1.23 3 1015.9 857 0.50 2.00 4 1015.9 858 0.10 1.95 5 1015.9 858 0.20 1.50 6 1025.0 858 0.30 1.20 7 1025.0 858 0.40 0.50 8 1025.0 858 0.35 0.70 9 1025.0 858 0.24 1.20 Thanks for your help. thanks [[alternative HTML version deleted]]
Dear Marna, This might come handy: http://www.geo.ut.ee/aasa/LOOM02331/heatmap_in_R.html Best regards, Michal *Michal Kubi?ta* Retail Analytics Support Europe Central/East Nielsen tel (+420) 242 438 737 www.nielsen.com [image: nielsen] <http://www.nielsen.com/> 2016-10-12 10:49 GMT+02:00 Marna Wagley <marna.wagley at gmail.com>:> Hi R Users, > Is it possible to visualize river flow in 3D (latitude, longitude with > respect to depth)? > The example of my data looks like. Any suggestions? > > > dat1 > long lat depth flow > 1 1015.9 857 1.00 1.50 > 2 1015.9 857 1.25 1.23 > 3 1015.9 857 0.50 2.00 > 4 1015.9 858 0.10 1.95 > 5 1015.9 858 0.20 1.50 > 6 1025.0 858 0.30 1.20 > 7 1025.0 858 0.40 0.50 > 8 1025.0 858 0.35 0.70 > 9 1025.0 858 0.24 1.20 > > Thanks for your help. > thanks > > [[alternative HTML version deleted]] > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/ > posting-guide.html > and provide commented, minimal, self-contained, reproducible code. >[[alternative HTML version deleted]]
Hi Marna,
Isn't the conventional way to visualize depth as shades of blue?
library(plotrix)
depth.col<-color.scale(dat1$depth,extremes=c("lightblue",blue"))
Then color the lon/lat rectangles with depth.col
Jim
On Wed, Oct 12, 2016 at 7:49 PM, Marna Wagley <marna.wagley at gmail.com>
wrote:> Hi R Users,
> Is it possible to visualize river flow in 3D (latitude, longitude with
> respect to depth)?
> The example of my data looks like. Any suggestions?
>
>> dat1
> long lat depth flow
> 1 1015.9 857 1.00 1.50
> 2 1015.9 857 1.25 1.23
> 3 1015.9 857 0.50 2.00
> 4 1015.9 858 0.10 1.95
> 5 1015.9 858 0.20 1.50
> 6 1025.0 858 0.30 1.20
> 7 1025.0 858 0.40 0.50
> 8 1025.0 858 0.35 0.70
> 9 1025.0 858 0.24 1.20
>
> Thanks for your help.
> thanks
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
On 12/10/2016 4:49 AM, Marna Wagley wrote:> Hi R Users, > Is it possible to visualize river flow in 3D (latitude, longitude with > respect to depth)? > The example of my data looks like. Any suggestions? > >> dat1 > long lat depth flow > 1 1015.9 857 1.00 1.50 > 2 1015.9 857 1.25 1.23 > 3 1015.9 857 0.50 2.00 > 4 1015.9 858 0.10 1.95 > 5 1015.9 858 0.20 1.50 > 6 1025.0 858 0.30 1.20 > 7 1025.0 858 0.40 0.50 > 8 1025.0 858 0.35 0.70 > 9 1025.0 858 0.24 1.20 > > Thanks for your help.It may be, but it's hard to give a nice looking graphic of that small dataset. You could try the rgl package and use plot3d to show spheres with radius depending on the flow rate, for example plot3d(cbind(long, lat, depth), type="s", col="blue", radius=flow/5) Duncan Murdoch
> On Oct 12, 2016, at 4:28 AM, Duncan Murdoch <murdoch.duncan at gmail.com> wrote: > > On 12/10/2016 4:49 AM, Marna Wagley wrote: >> Hi R Users, >> Is it possible to visualize river flow in 3D (latitude, longitude with >> respect to depth)? >> The example of my data looks like. Any suggestions? >> >>> dat1 >> long lat depth flow >> 1 1015.9 857 1.00 1.50 >> 2 1015.9 857 1.25 1.23 >> 3 1015.9 857 0.50 2.00 >> 4 1015.9 858 0.10 1.95 >> 5 1015.9 858 0.20 1.50 >> 6 1025.0 858 0.30 1.20 >> 7 1025.0 858 0.40 0.50 >> 8 1025.0 858 0.35 0.70 >> 9 1025.0 858 0.24 1.20 >> >> Thanks for your help. > > It may be, but it's hard to give a nice looking graphic of that small dataset. You could try the rgl package and use plot3d to show spheres with radius depending on the flow rate, for example > > plot3d(cbind(long, lat, depth), type="s", col="blue", radius=flow/5)A complementary option is to install the plot3D package which I see also has a plot3Drgl "co-package". The advantage to this option is the association with beautiful modeling packages that Karline Soetaert, Peter M. J. Herman, and Thomas Petzoldt have been offering to ecologists for the last decade. (Packages: deSolve, marelac, seacarb, AquaEnv) A lot of her work has been on flows within systems. I usually think of "flows" in rivers as being vector fields in an incompressible fluid (water) with 6 components per point, but you can also think of them as being scalar state variables. So I suppose you could be modeling something other than mass flows. (See Package::ReacTran for the R portal to that mathematical world.) Best; David Winsemius> > Duncan Murdoch > > ______________________________________________ > R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code.David Winsemius Alameda, CA, USA
Hi
With a small data set 3D is not really an option; reduce the number of
dimensions-- 2D conditional
library(lattice)
xyplot(flow ~ depth|factor(long), dat1, groups = lat, type =
c("p","r"), pch
= 16)
I had started with lat and long reversed doing EDA gave the above
? latitude effect
Regards
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mackay at northnet.com.au
-----Original Message-----
From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Marna Wagley
Sent: Wednesday, 12 October 2016 19:49
To: r-help mailing list
Subject: [R] can we visualize water flows with 3d in R?
Hi R Users,
Is it possible to visualize river flow in 3D (latitude, longitude with
respect to depth)?
The example of my data looks like. Any suggestions?
> dat1
long lat depth flow
1 1015.9 857 1.00 1.50
2 1015.9 857 1.25 1.23
3 1015.9 857 0.50 2.00
4 1015.9 858 0.10 1.95
5 1015.9 858 0.20 1.50
6 1025.0 858 0.30 1.20
7 1025.0 858 0.40 0.50
8 1025.0 858 0.35 0.70
9 1025.0 858 0.24 1.20
Thanks for your help.
thanks
[[alternative HTML version deleted]]
______________________________________________
R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Hi Marna
These are a few suggestions for EDA of your data
Firstly check it ? it is laborious work but it needs to be done ? use a proper
database if R is not suitable (memory,time)
dat <- read.csv("G:/1/example_subsetdata.csv")
dat <- subset(dat, !is.na(depth))
apply(dat,2,range,na.rm=T)
table(cut(dat$depth, breaks = seq(496.5, 498.5, 0.25), include.lowest = TRUE,
labels = F)) # needs more groups to even out numbers
# use hand drafted breaks
table(dat$Y)
dat$ldepth <- cut(dat$depth, breaks = seq(496.5, 498.5, 0.25), include.lowest
= TRUE, labels = F)
xyplot(Y ~ X|factor(ldepth), subset(dat, Y > 400), pch = ".",
as.table = TRUE)
xyplot(flow ~ X|factor(ldepth), subset(dat, Y > 400), pch = ".",
as.table = TRUE, col =1)
xyplot(flow ~ Y|factor(ldepth), subset(dat, Y > 400), pch = ".",
as.table = TRUE, col =1)
I have made pch = ?.? so that points are not overtopping. increase by by units
ie cex = 2 , cex =3 if you want to make them larger
Have a look a lattice::par.settings and names(trellis.par.get()) these will
make things easier if you you are using panel functions.
I did not go further but the next step in this would be lattice 3d plots
contourplot wireframe etc.
You will get more information in working with the full dataset
I have not really gone into spatials but you should look at the task view on
spatial analysis.
and the spatial SIG
How you proceed will be partly determined by your future requirements of which I
am not qualified.
Regards
Duncan
From: Marna Wagley [mailto:marna.wagley at gmail.com]
Sent: Friday, 14 October 2016 03:13
To: Duncan Mackay
Subject: Re: [R] can we visualize water flows with 3d in R?
Hi Duncan,
Thank you very much for the message. Indeed I do have a very big data set with a
grid size of 10cm * 10 cm. Every 10 cm*10cm grid has the value of X,Y, depth and
flow measurement. I have attached a subset of the data.
Thanks
On Wed, Oct 12, 2016 at 11:05 PM, Duncan Mackay <dulcalma at bigpond.com>
wrote:
Hi
With a small data set 3D is not really an option; reduce the number of
dimensions-- 2D conditional
library(lattice)
xyplot(flow ~ depth|factor(long), dat1, groups = lat, type =
c("p","r"), pch
= 16)
I had started with lat and long reversed doing EDA gave the above
? latitude effect
Regards
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mackay at northnet.com.au
-----Original Message-----
From: R-help [mailto:r-help-bounces at r-project.org] On Behalf Of Marna Wagley
Sent: Wednesday, 12 October 2016 19:49
To: r-help mailing list
Subject: [R] can we visualize water flows with 3d in R?
Hi R Users,
Is it possible to visualize river flow in 3D (latitude, longitude with
respect to depth)?
The example of my data looks like. Any suggestions?
> dat1
long lat depth flow
1 1015.9 857 1.00 1.50
2 1015.9 857 1.25 1.23
3 1015.9 857 0.50 2.00
4 1015.9 858 0.10 1.95
5 1015.9 858 0.20 1.50
6 1025.0 858 0.30 1.20
7 1025.0 858 0.40 0.50
8 1025.0 858 0.35 0.70
9 1025.0 858 0.24 1.20
Thanks for your help.
thanks
[[alternative HTML version deleted]]
______________________________________________
R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
[[alternative HTML version deleted]]