Re: [R-sig-Geo] ENVI data and R

2008-07-29 Thread PUJAN RAJ REGMI
This might help to mange the orientation of image: # To read ENVI formatcir.image<-("YOUR_ENVI_FILE")CIR.envi<-read.ENVI(cir.image,headerfile=paste(cir.image,".hdr",sep=""))# To Show imageCIR.envi.band1<-CIR.envi[,,1]CIR.envi.band1.s<-CIR.envi.band1[order(nrow(CIR.envi.band1):1),]CIR.envi.ban

Re: [R-sig-Geo] Spatial Resampling

2008-07-25 Thread PUJAN RAJ REGMI
ial > Resampling> From: [EMAIL PROTECTED]> To: r-sig-geo@stat.math.ethz.ch> CC: > [EMAIL PROTECTED]> > Pujan Raj Regmi <[EMAIL PROTECTED]> writes:> > > I have > a satellite image of 12m spatial resolution (4 bands). I want to> > resample > to coarser

Re: [R-sig-Geo] Spatial Resampling

2008-07-25 Thread PUJAN RAJ REGMI
Dear sir, Thank you very much for help. Pujan > Date: Fri, 25 Jul 2008 09:59:50 +0800> Subject: [R-sig-Geo] Spatial > Resampling> From: [EMAIL PROTECTED]> To: r-sig-geo@stat.math.ethz.ch> CC: > [EMAIL PROTECTED]> > Pujan Raj Regmi <[EMAIL PROTECTED]> writes:

Re: [R-sig-Geo] To fix the axis increment

2008-07-25 Thread PUJAN RAJ REGMI
Dear, Thank you very much. Pujan CC: [EMAIL PROTECTED]: [EMAIL PROTECTED]: [EMAIL PROTECTED]: Re: [R-sig-Geo] To fix the axis incrementDate: Fri, 25 Jul 2008 00:04:13 -0500Take a look at eqscplot(), a function in the MASS library. On Jul 22, 2008, at 8:00 AM, PUJAN RAJ REGMI wrote: Dear

[R-sig-Geo] Spatial Resampling

2008-07-23 Thread PUJAN RAJ REGMI
Dear all, I have a satellite image of 12m spatial resolution (4 bands). I want to resample to coarser level. i.e. to 18m, 24m,30m,36m and so on. The result I got from ENVI are not realistic so I want to do it in R using weighted average from neighboring pixels. I would appreciated any help. Thank

[R-sig-Geo] To fix the axis increment

2008-07-22 Thread PUJAN RAJ REGMI
Dear all, I am trying to plot 2-D scatter data in R, but i am having problem in fixing the increment in x and y axes. I am referring the function "axis equal" in Matlab by which the axes have same increment. Thanking you, Pujan _ e

[R-sig-Geo] Correlogram

2008-06-26 Thread PUJAN RAJ REGMI
Dear list, I am looking for R- package to compute correlogram and fit a model in it. Thanking you in advance Regards Pujan _ The other season of giving begins 6/24/08. Check out the i’m Talkathon. [[alternative HTML version de

[R-sig-Geo] Normality Test for Large sample

2008-06-20 Thread PUJAN RAJ REGMI
Dear list, I am trying to perform normality test for my data set. my data set is quite large N>13000 with 4 valriables under study , so When I tried to do normality test for each variable using lillie.test of nortest package the p values is almost zero "2.16e-16". So for all variable the test r

[R-sig-Geo] Semivariogram Plot

2008-06-09 Thread PUJAN RAJ REGMI
Dear list, I am trying to use semivariogram for texture classification for RS imagery. For this I am using the techniques of minimum distance algorithm to classy the image. To calculate minimum distance which is actually based on the gamma value of training site for each lag and those of pixel

[R-sig-Geo] Replacing the Matrix vlaues

2008-05-25 Thread PUJAN RAJ REGMI
Dear all, I am trying to replace the value in a matrix by specific value form noter matrix. The problem I am facing is: The size of one matrix is 30294 by 23(say matrix A) while the other matrix from which the values are intented to be extrated has the dimensions i.e. 14813 by 23(say matrix B).

Re: [R-sig-Geo] fit.variogram

2008-05-21 Thread PUJAN RAJ REGMI
Dear Sir, When I looked at the sent file in my mail box , the mail seems ok but the one i get it back from r-sig-geo@stat.math.ethz.ch it was completely meshed up. I am sorry for that. Now I am sending a copy of my code again: A<-read.table(file.choose(new=FALSE), sep = ",",skip=1,header=FALSE)

[R-sig-Geo] To fit variogram model

2008-05-21 Thread PUJAN RAJ REGMI
Dear list, I tried to fit the model variogram to the experimental semivariogram. I used the following code: A<-read.table(file.choose(new=FALSE), sep = ",",skip=1,header=FALSE) names(A) <- c("S.N","x", "y", "B1","B2","B3","B4") library(sp) library(gstat) coordinates(A) = ~x+y C<-vari

[R-sig-Geo] fit.variogram

2008-05-20 Thread PUJAN RAJ REGMI
Dear list, I tried to fit the model variogram to the experimental semivariogram. I used the following code: A<-read.table(file.choose(new=FALSE), sep = "",skip=8,header=FALSE) names(A) <- c("x", "y", "B1","B2","B3","B4") library(sp)library(gstat)coordinates(A) = ~x+yC<-variogram(B1~x+y+x*y+x*x+y

Re: [R-sig-Geo] variogram using gstat

2008-05-08 Thread PUJAN RAJ REGMI
> and width for omnidirectional and directional semivariogram > > (especially in 45 and 135 directions)?? as the lag in 0 and 90 > > directions is 12 where as in 45 and 135 is 16.97m.> > Though my question seems very basic, hope you find sometime to answer them> > Thanking you> > Pujan> >> > -

Re: [R-sig-Geo] variogram using gstat

2008-05-07 Thread PUJAN RAJ REGMI
hz.ch > Subject: Re: [R-sig-Geo] variogram using gstat > > PUJAN RAJ REGMI wrote: > > > > Dear sir, > > Ya, Now I am comparing the results from gslib and R/gsat. The shape of > > variograms exactly overlaps but the position of points varies, I don

Re: [R-sig-Geo] variogram using gstat

2008-05-07 Thread PUJAN RAJ REGMI
Dear sir, Ya, Now I am comparing the results from gslib and R/gsat. The shape of variograms exactly overlaps but the position of points varies, I don't know why? or it is normal! I first checked the result of omnidirectional variogram but still to have checked for directional cases. I hope the

Re: [R-sig-Geo] variogram using gstat

2008-05-07 Thread PUJAN RAJ REGMI
Dear Thierry, Thank you very much. Ya it worked after waiting for so long time. But i am bit confused with the result! as my data set have 12m interval both x and y direction but the result confused me. Please have a look at: A<-read.table(file.choose(new=FALSE), sep = "",skip=8,header=FALSE) n

[R-sig-Geo] variogram using gstat

2008-05-07 Thread PUJAN RAJ REGMI
Dear list, I am trying to plot semivariogram suing gstat packages. But as soon as i run the command variogram the programs hangs!! I don't know what happened! Actually in my data set i have many rows (96000). My data has 5 columns out of these first 2 are x and y coordinates and the rest 3 are t

Re: [R-sig-Geo] Scale Bar

2008-05-02 Thread PUJAN RAJ REGMI
le > Bar> > On Fri, 2 May 2008, PUJAN RAJ REGMI wrote:> > >> >> > Dear sir,> > My > data has been acquired from ENVI 4.0 and it has the following information > regarding pixel size and projection system:> > samples = 596lines = 440bands >

Re: [R-sig-Geo] Scale Bar

2008-05-02 Thread PUJAN RAJ REGMI
n. This is my source of error. Thanking you in advance Pujan > Date: Fri, 2 May 2008 19:49:07 +0200> From: [EMAIL PROTECTED]> To: [EMAIL > PROTECTED]> CC: r-sig-geo@stat.math.ethz.ch> Subject: RE: [R-sig-Geo] Scale > Bar> > On Fri, 2 May 2008, PUJAN RAJ REGMI wrote:&

Re: [R-sig-Geo] Scale Bar

2008-05-02 Thread PUJAN RAJ REGMI
; To: [EMAIL > PROTECTED]> CC: r-sig-geo@stat.math.ethz.ch> Subject: Re: [R-sig-Geo] Scale > Bar> > On Thu, 1 May 2008, PUJAN RAJ REGMI wrote:> > > Dear,> > > I followed > the instruction given by Roger. But i got the following error > > mess

[R-sig-Geo] Scale Bar

2008-05-01 Thread PUJAN RAJ REGMI
Dear, I followed the instruction given by Roger. But i got the following error message; coordinates(A) <- c("x", "y")> gridded(A) <- TRUE> greenchannel<-matrix(A[,3],nrow=155,ncol=191,byrow=TRUE)Error in as.vector(x, mode) : cannot coerce to vectorredchannel<-matrix(A[,4],nrow=155,ncol=191,byro

[R-sig-Geo] Adding Scale bar in the Image

2008-05-01 Thread PUJAN RAJ REGMI
Dear List, I would be thankful if anybody give me an way-out that How should can i add scale bar in the RGB color image. Actually i want to use this image for other purpose where i should mention the scale of the image. Thanking you in advance Pujan ___

Re: [R-sig-Geo] Trend

2008-03-29 Thread PUJAN RAJ REGMI
Dear all, I am trying to plot semivariogram for my data set consisting of 1373 values. As I plotted semivariogram, trend was seen. Then I decided to remove the trend so I tried to find the trend's coefficients for Linear, quadratic and cubic using geoR package. Fortunately its worked thanks to a

Re: [R-sig-Geo] Variogram Plot

2008-03-22 Thread PUJAN RAJ REGMI
Dear all, I think earlier i didn't clear myself. I ahve a data set of 3 colun and 14000 rows in which first and second column is X and Y coordinates and the third column is the value of samples of each coordinates points. Now I am trying to plot semi variogram for the given sample with respect

Re: [R-sig-Geo] Variogram Plot

2008-03-19 Thread PUJAN RAJ REGMI
Dear all, I think earlier i didn't clear myself. Now i am sending my data file and the code. My data set is attached with the mail where first and second column is X and Y coordinates and the column is the value of samples. I hope this time i am able to communicate. Thanking you Pujan Cod

Re: [R-sig-Geo] Variogram Plot

2008-03-18 Thread PUJAN RAJ REGMI
Thank you very much for response. Here i have attached the my code below . Pls have a look and suggest me the solution. Thanking you in advance Yours Pujan mol.gr <- read.table(“D:/Pujan_IUPWARE/Second_Year/Thesis/mol.dat”) mol.gr[,"x"] <-mol.gr[,1] mol.gr[,"y"] <-mol.gr[,2] mol.gr[,"Green Ban

[R-sig-Geo] Variogram Plot

2008-03-17 Thread PUJAN RAJ REGMI
Dear all, I am trying to plot semivarigram using the function "variog" but i am getting following error message "cannot allocate vector of size 757.6 Mb" I suppose my data set seems large?? it has 14090 samples. I will be thankful if anybody suggest me the alternative ways to overcome the probl

Re: [R-sig-Geo] Rimage Package

2008-02-26 Thread PUJAN RAJ REGMI
s fairly simple and matches with what you already tried).> > Hope this > helps,> > -- > Gregory (Greg) L. Snow Ph.D.> Statistical Data Center> > Intermountain Healthcare> [EMAIL PROTECTED]> (801) 408-8111> > > > > > -Original Message-> &g

[R-sig-Geo] Rimage Package

2008-02-25 Thread PUJAN RAJ REGMI
Thanks a lot it gives me a wayout but I still want favour from you. Could you please explain how can i prepared three mat file required for the following code? imagematrix(mat, type=NULL, ncol=dim(mat)[1], nrow=dim(mat)[2], noclipping=FALSE) I have three R files named as Redchannel, Greenchann

[R-sig-Geo] RGB Image Plotting!

2008-02-25 Thread PUJAN RAJ REGMI
Dear Madam/ Sir I am New user of R. I am looking for the way to plot true color image from matrix data of of order (m,n,3) as in MAT LAB where 3 represent the three column matrix where red, green and blue values of pixels are stored. I hope you understand my question. Looking for your help for