Re: [R] How to plot wind direction and strength field
Thanks Greg, I'll give it a whirl and let you know how it goes (although it might not be till after the bank holiday weekend now sorry for delay!) Thanks for replying and spending time on this problem. All the best, Jenny On Thu, 1 May 2008, Greg Snow wrote: Does the following do what you want (or at least move in that direction)? u - array(NA,c(5,8)) t - seq(from=0.5, to=0.11,length=15) t2 - seq(from=(-0.7),to=(-0.1),length=25) u[1:15] - t u[16:40] - t2 v - array(NA,c(5,8)) y - seq(from=(-0.9), to=(-0.01),length=40) v[1:40] - y library(TeachingDemos) library(maps) map('state',col='grey') lats - seq( 27.5, 48.7, length=5 ) longs - seq( -123, -73, length=8 ) x - longs[ col(u) ] y - lats[ row(u) ] speed - sqrt( u*u+v*v ) dir - atan2(v,u) my.symbols(x,y,ms.arrows, angle=dir, r=speed, add=TRUE, length=.05) Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jenny Barnes Sent: Tuesday, April 29, 2008 6:26 AM To: Jim Lemon Cc: r-help@r-project.org Subject: Re: [R] How to plot wind direction and strength field Thanks Jim - here is some data: u - array(NA,c(5,8)) t - seq(from=0.5, to=0.11,length=15) t2 - seq(from=(-0.7),to=(-0.1),length=25) u[1:15] - t u[16:40] - t2 v - array(NA,c(5,8)) y - seq(from=(-0.9), to=(-0.01),length=40) v[1:40] - y I've made up the data but it's similar magnitude and the longitude direction is u and latitude direction is v as in the real data. I really appreciate your help! Jenny On Tue, 29 Apr 2008, Jim Lemon wrote: Jenny Barnes wrote: Hi Jim, I would like to plot something like figure 2 on this webpage: http://www.cnrfc.noaa.gov/storm_summaries/jan1997storms.php My data is very large - covering the whole globe at 2.5deg resolution so longitude=144 girds, latitude=73 grids and time=32 years - hard to give you that data..Would it help to give you a couple of grid squares worth of data for one year? Okay, Figures 2 and 3 look like 5 degree squares, and given a larger plot, 2.5 degree should be okay. Do you want curved arrows? That would take some programming as the standard arrows are straight. Also, the arrows in the figures all seem to be the same length. I can make the lengths proportional to wind speed. Yes, some data would be helpful, as getting the right plot usually involves trying out real data. Doesn't matter if it's faked as long as it has the same format and similar numbers to the real thing. Jim __ R-help@r-project.org mailing list 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@r-project.org mailing list 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.
Re: [R] How to plot wind direction and strength field
Does the following do what you want (or at least move in that direction)? u - array(NA,c(5,8)) t - seq(from=0.5, to=0.11,length=15) t2 - seq(from=(-0.7),to=(-0.1),length=25) u[1:15] - t u[16:40] - t2 v - array(NA,c(5,8)) y - seq(from=(-0.9), to=(-0.01),length=40) v[1:40] - y library(TeachingDemos) library(maps) map('state',col='grey') lats - seq( 27.5, 48.7, length=5 ) longs - seq( -123, -73, length=8 ) x - longs[ col(u) ] y - lats[ row(u) ] speed - sqrt( u*u+v*v ) dir - atan2(v,u) my.symbols(x,y,ms.arrows, angle=dir, r=speed, add=TRUE, length=.05) Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare [EMAIL PROTECTED] (801) 408-8111 -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Jenny Barnes Sent: Tuesday, April 29, 2008 6:26 AM To: Jim Lemon Cc: r-help@r-project.org Subject: Re: [R] How to plot wind direction and strength field Thanks Jim - here is some data: u - array(NA,c(5,8)) t - seq(from=0.5, to=0.11,length=15) t2 - seq(from=(-0.7),to=(-0.1),length=25) u[1:15] - t u[16:40] - t2 v - array(NA,c(5,8)) y - seq(from=(-0.9), to=(-0.01),length=40) v[1:40] - y I've made up the data but it's similar magnitude and the longitude direction is u and latitude direction is v as in the real data. I really appreciate your help! Jenny On Tue, 29 Apr 2008, Jim Lemon wrote: Jenny Barnes wrote: Hi Jim, I would like to plot something like figure 2 on this webpage: http://www.cnrfc.noaa.gov/storm_summaries/jan1997storms.php My data is very large - covering the whole globe at 2.5deg resolution so longitude=144 girds, latitude=73 grids and time=32 years - hard to give you that data..Would it help to give you a couple of grid squares worth of data for one year? Okay, Figures 2 and 3 look like 5 degree squares, and given a larger plot, 2.5 degree should be okay. Do you want curved arrows? That would take some programming as the standard arrows are straight. Also, the arrows in the figures all seem to be the same length. I can make the lengths proportional to wind speed. Yes, some data would be helpful, as getting the right plot usually involves trying out real data. Doesn't matter if it's faked as long as it has the same format and similar numbers to the real thing. Jim __ R-help@r-project.org mailing list 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@r-project.org mailing list 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] How to plot wind direction and strength field
Dear R-help community, I have searched through the archives and not been able ot find any advice on how to plot a wind field with one arrow per grid square with the arrow pointing in the direction of the wind and it's size proportional to the wind strength. I have the wind speed data in arrays of [lon,lat,uwind] and [lon,lat,vwind] so it is broken down into u and v components. How do I plot it though?!?! Any suggestions very wecome indeed - I seem to have hit a brick wall. All the best, Jenny ~~ Jennifer Barnes PhD student: long range drought prediction Climate Extremes Group Department of Space and Climate Physics University College London Holmbury St Mary Dorking, Surrey, RH5 6NT Web: http://climate.mssl.ucl.ac.uk __ R-help@r-project.org mailing list 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.
Re: [R] How to plot wind direction and strength field
Jenny Barnes wrote: Dear R-help community, I have searched through the archives and not been able ot find any advice on how to plot a wind field with one arrow per grid square with the arrow pointing in the direction of the wind and it's size proportional to the wind strength. I have the wind speed data in arrays of [lon,lat,uwind] and [lon,lat,vwind] so it is broken down into u and v components. How do I plot it though?!?! Any suggestions very wecome indeed - I seem to have hit a brick wall. Hi Jenny, Some time ago, there was a request for a direction field plot. I programmed a basic function, but perhaps the person found another solution, for I never heard any more of it. However, this is not too hard to do in R, even adding the arrows to a geographical plot. Could you post some data and perhaps a link to an example of what output you would like? Jim __ R-help@r-project.org mailing list 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.
Re: [R] How to plot wind direction and strength field
Hi Jim, I would like to plot something like figure 2 on this webpage: http://www.cnrfc.noaa.gov/storm_summaries/jan1997storms.php My data is very large - covering the whole globe at 2.5deg resolution so longitude=144 girds, latitude=73 grids and time=32 years - hard to give you that data..Would it help to give you a couple of grid squares worth of data for one year? Thanks, Jenny On Tue, 29 Apr 2008, Jim Lemon wrote: Jenny Barnes wrote: Dear R-help community, I have searched through the archives and not been able ot find any advice on how to plot a wind field with one arrow per grid square with the arrow pointing in the direction of the wind and it's size proportional to the wind strength. I have the wind speed data in arrays of [lon,lat,uwind] and [lon,lat,vwind] so it is broken down into u and v components. How do I plot it though?!?! Any suggestions very wecome indeed - I seem to have hit a brick wall. Hi Jenny, Some time ago, there was a request for a direction field plot. I programmed a basic function, but perhaps the person found another solution, for I never heard any more of it. However, this is not too hard to do in R, even adding the arrows to a geographical plot. Could you post some data and perhaps a link to an example of what output you would like? Jim __ R-help@r-project.org mailing list 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.
Re: [R] How to plot wind direction and strength field
Jenny, Have a look at the R Newsletter Volume 3/2, October 2003 Regards, Sean Jenny Barnes wrote: Dear R-help community, I have searched through the archives and not been able ot find any advice on how to plot a wind field with one arrow per grid square with the arrow pointing in the direction of the wind and it's size proportional to the wind strength. I have the wind speed data in arrays of [lon,lat,uwind] and [lon,lat,vwind] so it is broken down into u and v components. How do I plot it though?!?! Any suggestions very wecome indeed - I seem to have hit a brick wall. All the best, Jenny ~~ Jennifer Barnes PhD student: long range drought prediction Climate Extremes Group Department of Space and Climate Physics University College London Holmbury St Mary Dorking, Surrey, RH5 6NT Web: http://climate.mssl.ucl.ac.uk __ R-help@r-project.org mailing list 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. -- View this message in context: http://www.nabble.com/generic-question--%3E-Genomics-with-R-tp16954827p16958167.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org mailing list 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.
Re: [R] How to plot wind direction and strength field
Check out: http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=80 and RSiteSearch(quiver) On Tue, Apr 29, 2008 at 6:08 AM, Jenny Barnes [EMAIL PROTECTED] wrote: Dear R-help community, I have searched through the archives and not been able ot find any advice on how to plot a wind field with one arrow per grid square with the arrow pointing in the direction of the wind and it's size proportional to the wind strength. I have the wind speed data in arrays of [lon,lat,uwind] and [lon,lat,vwind] so it is broken down into u and v components. How do I plot it though?!?! Any suggestions very wecome indeed - I seem to have hit a brick wall. All the best, Jenny ~~ Jennifer Barnes PhD student: long range drought prediction Climate Extremes Group Department of Space and Climate Physics University College London Holmbury St Mary Dorking, Surrey, RH5 6NT Web: http://climate.mssl.ucl.ac.uk __ R-help@r-project.org mailing list 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@r-project.org mailing list 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.
Re: [R] How to plot wind direction and strength field
Thanks Jim - here is some data: u - array(NA,c(5,8)) t - seq(from=0.5, to=0.11,length=15) t2 - seq(from=(-0.7),to=(-0.1),length=25) u[1:15] - t u[16:40] - t2 v - array(NA,c(5,8)) y - seq(from=(-0.9), to=(-0.01),length=40) v[1:40] - y I've made up the data but it's similar magnitude and the longitude direction is u and latitude direction is v as in the real data. I really appreciate your help! Jenny On Tue, 29 Apr 2008, Jim Lemon wrote: Jenny Barnes wrote: Hi Jim, I would like to plot something like figure 2 on this webpage: http://www.cnrfc.noaa.gov/storm_summaries/jan1997storms.php My data is very large - covering the whole globe at 2.5deg resolution so longitude=144 girds, latitude=73 grids and time=32 years - hard to give you that data..Would it help to give you a couple of grid squares worth of data for one year? Okay, Figures 2 and 3 look like 5 degree squares, and given a larger plot, 2.5 degree should be okay. Do you want curved arrows? That would take some programming as the standard arrows are straight. Also, the arrows in the figures all seem to be the same length. I can make the lengths proportional to wind speed. Yes, some data would be helpful, as getting the right plot usually involves trying out real data. Doesn't matter if it's faked as long as it has the same format and similar numbers to the real thing. Jim __ R-help@r-project.org mailing list 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.
Re: [R] How to plot wind direction and strength field
Hi Jim, I forgot to say that I don't want to create curved arrows, just straight - and I would like them to be proportional to the magnitude :o) Thanks again for your time! Jenny On Tue, 29 Apr 2008, Jenny Barnes wrote: Hi Jim, I would like to plot something like figure 2 on this webpage: http://www.cnrfc.noaa.gov/storm_summaries/jan1997storms.php My data is very large - covering the whole globe at 2.5deg resolution so longitude=144 girds, latitude=73 grids and time=32 years - hard to give you that data..Would it help to give you a couple of grid squares worth of data for one year? Thanks, Jenny On Tue, 29 Apr 2008, Jim Lemon wrote: Jenny Barnes wrote: Dear R-help community, I have searched through the archives and not been able ot find any advice on how to plot a wind field with one arrow per grid square with the arrow pointing in the direction of the wind and it's size proportional to the wind strength. I have the wind speed data in arrays of [lon,lat,uwind] and [lon,lat,vwind] so it is broken down into u and v components. How do I plot it though?!?! Any suggestions very wecome indeed - I seem to have hit a brick wall. Hi Jenny, Some time ago, there was a request for a direction field plot. I programmed a basic function, but perhaps the person found another solution, for I never heard any more of it. However, this is not too hard to do in R, even adding the arrows to a geographical plot. Could you post some data and perhaps a link to an example of what output you would like? Jim __ R-help@r-project.org mailing list 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.