Help for STFDF Object Creation

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
3 messages Options
Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Help for STFDF Object Creation

R-sig-geo mailing list
Hi,
I am trying to create a STFDF object and draw variogram. However, I am getting several errors. Here is my r code:
#create count data
c0157ts <- as.numeric(c(t(c0157)))
c0001ts <- as.numeric(c(t(c0001)))

#create coordinates
  lat <- 0
  lon <- 0
  lat[1] <- 34.0793684
  lon[1] <- -81.1301722
  lat[2] <- 33.979127
  lon[2] <- -81.321166398
  sp <- data.frame(lat,lon)
  coordinates(sp)=~lon+lat
  projection(sp)=CRS("+init=epsg:4326")

#create time object
time_index <- seq(from = as.POSIXct("2017-01-02 01:00", tz = 'UTC'), to = as.POSIXct("2017-03-22 00:00", tz = 'UTC'), by = "hour")

#combine the data
mydata <- data.frame("count"=c(c0157ts,c0001ts))

#create STFDF object
stfdf = STFDF(ozoneSP, time_index, mydata)
#trying to create variogram
    var <- variogramST(count~1,data=stfdf,tunit="hours",assumeRegular=F,na.omit=T)which gives me following error:0%Error in apply(do.call(cbind, lapply(ret, function(x) x$np)), 1, sum,  :
      dim(X) must have a positive length
I guess my STFDF object creation is not correct. Can anyone give us any hint? Which elements in STFDF object is not correct.Here is the description of STFDF object:> str(stfdf)
Formal class 'STFDF' [package "spacetime"] with 4 slots
  ..@ data   :'data.frame':     3792 obs. of  1 variable:
  .. ..$ count: num [1:3792] 40 32 64 40 55 89 71 43 69 44 ...
  ..@ sp     :Formal class 'SpatialPoints' [package "sp"] with 3 slots
  .. .. ..@ coords     : num [1:2, 1:2] -9031369 -9052631 4015522 4002120
  .. .. .. ..- attr(*, "dimnames")=List of 2
  .. .. .. .. ..$ : NULL
  .. .. .. .. ..$ : chr [1:2] "lon" "lat"
  .. .. ..@ bbox       : num [1:2, 1:2] -9052631 4002120 -9031369 4015522
  .. .. .. ..- attr(*, "dimnames")=List of 2
  .. .. .. .. ..$ : chr [1:2] "lon" "lat"
  .. .. .. .. ..$ : chr [1:2] "min" "max"
  .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
  .. .. .. .. ..@ projargs: chr "+init=epsg:3395 +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
  ..@ time   :An ‘xts’ object on 2017-01-02 01:00:00/2017-03-22 containing:
  Data: int [1:1896, 1] 1 2 3 4 5 6 7 8 9 10 ...
 - attr(*, "dimnames")=List of 2
  ..$ : NULL
  ..$ : chr "timeIndex"
  Indexed by objects of class: [POSIXct,POSIXt] TZ: UTC
  xts Attributes:  
 NULL
Thanks for your help!
With Regards,Rukna
        [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: Help for STFDF Object Creation

edzer
Please provide a script of an example that we can reproduce, preferably
using data that is loaded from a package (to safe other people's time),
otherwise it is impossible for anyone to help you, or unlikely somebody
will. Your script does not create count data but assumes it is there;
you don't show which package you loaded; it is unclear what your problem is.

Your second post on this list came through completely scrambled, and
within 24 hours. Don't expect that reposting, or scrambling, helps.

On 23/04/17 06:02, Ankur Sarker via R-sig-Geo wrote:

> Hi,
> I am trying to create a STFDF object and draw variogram. However, I am getting several errors. Here is my r code:
> #create count data
> c0157ts <- as.numeric(c(t(c0157)))
> c0001ts <- as.numeric(c(t(c0001)))
>
> #create coordinates
>   lat <- 0
>   lon <- 0
>   lat[1] <- 34.0793684
>   lon[1] <- -81.1301722
>   lat[2] <- 33.979127
>   lon[2] <- -81.321166398
>   sp <- data.frame(lat,lon)
>   coordinates(sp)=~lon+lat
>   projection(sp)=CRS("+init=epsg:4326")
>
> #create time object
> time_index <- seq(from = as.POSIXct("2017-01-02 01:00", tz = 'UTC'), to = as.POSIXct("2017-03-22 00:00", tz = 'UTC'), by = "hour")
>
> #combine the data
> mydata <- data.frame("count"=c(c0157ts,c0001ts))
>
> #create STFDF object
> stfdf = STFDF(ozoneSP, time_index, mydata)
> #trying to create variogram
>     var <- variogramST(count~1,data=stfdf,tunit="hours",assumeRegular=F,na.omit=T)which gives me following error:0%Error in apply(do.call(cbind, lapply(ret, function(x) x$np)), 1, sum,  :
>       dim(X) must have a positive length
> I guess my STFDF object creation is not correct. Can anyone give us any hint? Which elements in STFDF object is not correct.Here is the description of STFDF object:> str(stfdf)
> Formal class 'STFDF' [package "spacetime"] with 4 slots
>   ..@ data   :'data.frame':     3792 obs. of  1 variable:
>   .. ..$ count: num [1:3792] 40 32 64 40 55 89 71 43 69 44 ...
>   ..@ sp     :Formal class 'SpatialPoints' [package "sp"] with 3 slots
>   .. .. ..@ coords     : num [1:2, 1:2] -9031369 -9052631 4015522 4002120
>   .. .. .. ..- attr(*, "dimnames")=List of 2
>   .. .. .. .. ..$ : NULL
>   .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>   .. .. ..@ bbox       : num [1:2, 1:2] -9052631 4002120 -9031369 4015522
>   .. .. .. ..- attr(*, "dimnames")=List of 2
>   .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>   .. .. .. .. ..$ : chr [1:2] "min" "max"
>   .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
>   .. .. .. .. ..@ projargs: chr "+init=epsg:3395 +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
>   ..@ time   :An ‘xts’ object on 2017-01-02 01:00:00/2017-03-22 containing:
>   Data: int [1:1896, 1] 1 2 3 4 5 6 7 8 9 10 ...
>  - attr(*, "dimnames")=List of 2
>   ..$ : NULL
>   ..$ : chr "timeIndex"
>   Indexed by objects of class: [POSIXct,POSIXt] TZ: UTC
>   xts Attributes:  
>  NULL
> Thanks for your help!
> With Regards,Rukna
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> [hidden email]
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
--
Edzer Pebesma
Institute for Geoinformatics  (ifgi),  University of Münster
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software:   http://www.jstatsoft.org/
Computers & Geosciences:   http://elsevier.com/locate/cageo/


_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

signature.asc (484 bytes) Download Attachment
Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

Re: Help for STFDF Object Creation

R-sig-geo mailing list
Thanks for your reply. Sorry for my second post, I tried to be more clear.

Here is link of my complete source code, data and R instance:https://drive.google.com/drive/folders/0B3tdiSZ9hBj4dVZXOXItUjdrRlE?usp=sharing


Hope, it would be easier to understand my problem.

    On Sunday, April 23, 2017 5:25 PM, Edzer Pebesma <[hidden email]> wrote:
 

 Please provide a script of an example that we can reproduce, preferably
using data that is loaded from a package (to safe other people's time),
otherwise it is impossible for anyone to help you, or unlikely somebody
will. Your script does not create count data but assumes it is there;
you don't show which package you loaded; it is unclear what your problem is.

Your second post on this list came through completely scrambled, and
within 24 hours. Don't expect that reposting, or scrambling, helps.

On 23/04/17 06:02, Ankur Sarker via R-sig-Geo wrote:

> Hi,
> I am trying to create a STFDF object and draw variogram. However, I am getting several errors. Here is my r code:
> #create count data
> c0157ts <- as.numeric(c(t(c0157)))
> c0001ts <- as.numeric(c(t(c0001)))
>
> #create coordinates
>  lat <- 0
>  lon <- 0
>  lat[1] <- 34.0793684
>  lon[1] <- -81.1301722
>  lat[2] <- 33.979127
>  lon[2] <- -81.321166398
>  sp <- data.frame(lat,lon)
>  coordinates(sp)=~lon+lat
>  projection(sp)=CRS("+init=epsg:4326")
>
> #create time object
> time_index <- seq(from = as.POSIXct("2017-01-02 01:00", tz = 'UTC'), to = as.POSIXct("2017-03-22 00:00", tz = 'UTC'), by = "hour")
>
> #combine the data
> mydata <- data.frame("count"=c(c0157ts,c0001ts))
>
> #create STFDF object
> stfdf = STFDF(ozoneSP, time_index, mydata)
> #trying to create variogram
>    var <- variogramST(count~1,data=stfdf,tunit="hours",assumeRegular=F,na.omit=T)which gives me following error:0%Error in apply(do.call(cbind, lapply(ret, function(x) x$np)), 1, sum,  :
>      dim(X) must have a positive length
> I guess my STFDF object creation is not correct. Can anyone give us any hint? Which elements in STFDF object is not correct.Here is the description of STFDF object:> str(stfdf)
> Formal class 'STFDF' [package "spacetime"] with 4 slots
>  ..@ data  :'data.frame':    3792 obs. of  1 variable:
>  .. ..$ count: num [1:3792] 40 32 64 40 55 89 71 43 69 44 ...
>  ..@ sp    :Formal class 'SpatialPoints' [package "sp"] with 3 slots
>  .. .. ..@ coords    : num [1:2, 1:2] -9031369 -9052631 4015522 4002120
>  .. .. .. ..- attr(*, "dimnames")=List of 2
>  .. .. .. .. ..$ : NULL
>  .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>  .. .. ..@ bbox      : num [1:2, 1:2] -9052631 4002120 -9031369 4015522
>  .. .. .. ..- attr(*, "dimnames")=List of 2
>  .. .. .. .. ..$ : chr [1:2] "lon" "lat"
>  .. .. .. .. ..$ : chr [1:2] "min" "max"
>  .. .. ..@ proj4string:Formal class 'CRS' [package "sp"] with 1 slot
>  .. .. .. .. ..@ projargs: chr "+init=epsg:3395 +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0"
>  ..@ time  :An ‘xts’ object on 2017-01-02 01:00:00/2017-03-22 containing:
>  Data: int [1:1896, 1] 1 2 3 4 5 6 7 8 9 10 ...
>  - attr(*, "dimnames")=List of 2
>  ..$ : NULL
>  ..$ : chr "timeIndex"
>  Indexed by objects of class: [POSIXct,POSIXt] TZ: UTC
>  xts Attributes: 
>  NULL
> Thanks for your help!
> With Regards,Rukna
>     [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> [hidden email]
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

--
Edzer Pebesma
Institute for Geoinformatics  (ifgi),  University of Münster
Heisenbergstraße 2, 48149 Münster, Germany; +49 251 83 33081
Journal of Statistical Software:  http://www.jstatsoft.org/
Computers & Geosciences:  http://elsevier.com/locate/cageo/
_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo

   
        [[alternative HTML version deleted]]

_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Loading...