An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110901/9ebd912a/attachment.pl>
convert to grid file
6 messages · azam jaafari, PIKAL Petr
Hi
Hi
I computed probability in each cell.
I have:
[99883,] -0.0062412957690
[99884,] -0.0062412957690
[99885,] -0.0062412957690
[99886,] -0.0062412957690
[99887,] -0.0062412957690
[99888,] -0.0062412957690
[99889,] 0.9909126638948
[99890,] 0.9909126638948
[99891,] 0.9909126638948
[99892,] 0.9909126638948
[99893,] 0.9909126638948
[99894,] 0.9909126638948
[99895,] 0.9909126638948
[99896,] 0.9909126638948
[99897,] 0.9909126638948
[99898,] 0.9909126638948
[99899,] 0.9909126638948
[99900,] 0.9909126638948
[99901,] 0.9909126638948
[99902,] 0.9909126638948
[99903,] 0.9909126638948
[99999,] -0.0062412957690
[ reached getOption("max.print") -- omitted 839931 rows ]]
How want to convert this matrix to a grid file with 970*960 pixel.
Assuming your object is one column matrix called mat1 You can simply change its dimension dim(mat) <- c(970, 960) regards Petr
Thanks alot [[alternative HTML version deleted]]
______________________________________________ R-help at 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.
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110901/36104474/attachment.pl>
Hi
Thank you Petr It work. Now I have a matrix 970*960. If I want to convert to spatial grid (each pixel has x and y coordinate). How can I do?
I do not understand. What do you want to do with your data? Maybe you could consult spatial package or CRAN Task views. One option could be to make vectors of row and columns coordinates. But it depends on what you want to do with your data. Regards Petr
Thanks From: Petr PIKAL <petr.pikal at precheza.cz> To: azam jaafari <azamjaafari at yahoo.com> Cc: R-help <r-help at r-project.org> Sent: Thursday, September 1, 2011 8:59 AM Subject: Re: [R] convert to grid file Hi
Hi
I computed probability in each cell.
I have:
[99883,] -0.0062412957690
[99884,] -0.0062412957690
[99885,] -0.0062412957690
[99886,] -0.0062412957690
[99887,] -0.0062412957690
[99888,] -0.0062412957690
[99889,] 0.9909126638948
[99890,] 0.9909126638948
[99891,] 0.9909126638948
[99892,] 0.9909126638948
[99893,] 0.9909126638948
[99894,] 0.9909126638948
[99895,] 0.9909126638948
[99896,] 0.9909126638948
[99897,] 0.9909126638948
[99898,] 0.9909126638948
[99899,] 0.9909126638948
[99900,] 0.9909126638948
[99901,] 0.9909126638948
[99902,] 0.9909126638948
[99903,] 0.9909126638948
[99999,] -0.0062412957690
[ reached getOption("max.print") -- omitted 839931 rows ]]
How want to convert this matrix to a grid file with 970*960 pixel.
Assuming your object is one column matrix called mat1 You can simply change its dimension dim(mat) <- c(970, 960) regards Petr
Thanks alot [[alternative HTML version deleted]]
______________________________________________ R-help at 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.
An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20110901/26a865be/attachment.pl>
Hi I do not know sp package therefore I can not help you. Anyway without knowing "what is wrong" and without some reproducible example you hardly get any reasonable answer. It seems to me that your d list especially x and y is not what you expect to be. see str(d) 430071.4887:460006.8067 will resilt in row of integers from 430071 to 460006. I believe it is not what you want. Regards Petr
Thank you I want to make a map from my spatial data that I can show it in GIS. I used the package sp and d<- list(x=430071.4887:460006.8067, y =3390040.0591:3420006.2701, z = matrix(mat, 970, 960)) gt<- GridTopology(cellcentre.offset = c(d$x[1], d$y[1]), cellsize=c(diff(d$x[1:2]), diff(d$y[1:2])), cells.dim = dim(d$z)) grdatts<-SpatialGridDataFrame(gt, data.frame(depth = as.vector(d$z[,ncol(d$z):1]))) image(grd.atts, axes=TRUE) grd<-SpatialGridDataFrame(gt, data.frame(depth = as.vector(d$z[,ncol(d$z):1]))) It doesn't work and I think it is wrong. Thanks From: Petr PIKAL <petr.pikal at precheza.cz> To: azam jaafari <azamjaafari at yahoo.com> Cc: R-help <r-help at r-project.org> Sent: Thursday, September 1, 2011 10:22 AM Subject: Re: [R] convert to grid file Hi
Thank you Petr It work. Now I have a matrix 970*960. If I want to convert to spatial grid
(each
pixel has x and y coordinate). How can I do?
I do not understand. What do you want to do with your data? Maybe you could consult spatial package or CRAN Task views. One option could be to
make vectors of row and columns coordinates. But it depends on what you want to do with your data. Regards Petr
Thanks From: Petr PIKAL <petr.pikal at precheza.cz> To: azam jaafari <azamjaafari at yahoo.com> Cc: R-help <r-help at r-project.org> Sent: Thursday, September 1, 2011 8:59 AM Subject: Re: [R] convert to grid file Hi
Hi
I computed probability in each cell.
I have:
[99883,] -0.0062412957690
[99884,] -0.0062412957690
[99885,] -0.0062412957690
[99886,] -0.0062412957690
[99887,] -0.0062412957690
[99888,] -0.0062412957690
[99889,] 0.9909126638948
[99890,] 0.9909126638948
[99891,] 0.9909126638948
[99892,] 0.9909126638948
[99893,] 0.9909126638948
[99894,] 0.9909126638948
[99895,] 0.9909126638948
[99896,] 0.9909126638948
[99897,] 0.9909126638948
[99898,] 0.9909126638948
[99899,] 0.9909126638948
[99900,] 0.9909126638948
[99901,] 0.9909126638948
[99902,] 0.9909126638948
[99903,] 0.9909126638948
[99999,] -0.0062412957690
[ reached getOption("max.print") -- omitted 839931 rows ]]
How want to convert this matrix to a grid file with 970*960 pixel.
Assuming your object is one column matrix called mat1 You can simply change its dimension dim(mat) <- c(970, 960) regards Petr
Thanks alot [[alternative HTML version deleted]]
______________________________________________ R-help at 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.