search for: lat_ts

I have two different datasets 1) is in monthly format (obs) 2) yearly format (model) in obs I have 84 files ( 2003:2009)for different months & in model I have 4 different files which has yearly data (2005:2008) So for calculating my requirement I need these both data sets. The sample calculations are as follows file_o<-list.files(path=' ', pattern="0.2.text") #

...t to Berhmann to ensure equal area distribution of the points. In R, I am using: spPoint <- SpatialPoints(coords=coordinates(Data),proj4string=CRS("+proj=longlat +datum=WGS84")) and then reprojecting it to Berhmann with: spPointReprj=spTransform(Data,CRS("+proj=cea +lon_0=0 +lat_ts=30 +x_0=0 +y_0=0 +ellps=WGS84")) If I put the two outputs of the reprojections in the same map, they are off by few meters in the tropics by up to 40km in the poles. I decided to investigate how the reprojections are done and my calculations are different from both R and ArcGis: First...

Dear all, We are trying the code from the paper (http://www.jstatsoft.org/v40/i01/) to build stars for each point on the map (plyr package). It uses ggplot2 package and works fine for ozone data but we can't get it work for our data (probably due to misunderstanding of some issues). Our data is SO2 time series for each point on the map with 0.25 degrees step. It has both NA, negative and