search for: ndvi

Hello, I have conducted an SEM in which the resultant standardized path coefficients are much higher than would be expected from the raw correlation matrix. To explore further, I stripped the model down to a simple bivariate relationship between two variables (NDVI, and species richness), where it's my understanding that the SEM's standardized path coefficient should equal the correlation coefficient. The datafile can be downloaded at http://labs.bio.unc.edu/Hurlbert/pubs/NDVI_lep_data.csv The model specification file can be downloaded at http://labs...

Hi People! I have following concern consisting of some steps to do in R: I have an ascii file (table) consisting of many columns and rows. 1. I would like to order all values of the columns one under each other. It will begin with column 1, then column 2 under column 1, column 3 under column 2 etc. until at the end there is only 1 column. How do I do it? 2. Second problem is to make a

Hi All, I'm running into a problem with GAM (in the MGCV package). When I try to estimate the model, I get the following error message: 1> fit <- gam(ndvi~s(rain)+s(temp)+s(rainl1)+s(rainl2)+s(rainxY)+s(rainl1xY)+s(rainl2xY)+s(tempxY), data=dsub, weights=wvec) Error in while (mean(ldxx/(ldxx + ldss)) > 0.4) { : missing value where TRUE/FALSE needed Using traceback, I get the following output (which I don't understand) 1> tr...

...order legend.txt <- c("NDII", "Inverse KBDI") legend(2001.8,0.29, legend.txt,col=c(1,2), lty=1) par(new=T) plot(ts.Origi[,2], yaxt="n", type="l", col=2, las=1, ylab="") axis(4) plot(ts.NDVIKB[,1],ylab='NDVI', main=name) # Add the legend text in the right order legend.txt <- c("NDVI", "Inverse KBDI") legend(2001.8,0.42, legend.txt,col=c(1,2), lty=1) par(new=T) plot(ts.NDVIKB[,2], yaxt="n"...

...t it would help much with my problem: ######## # It made of 4 functions (sorry, it's french): ########################################################################## ########################################################################## ### Ensemble des fonctions pour faire les images NDVI rouge et verte ### ########################################################################## ###### Bastien Ferland-Raymond, 5 oct 2010 ############################# ########################################################################## ######## ## Simplement rouler le script au complet ##...

...utorials, searching for my error code, looking at the source code. However, I have not been able to find any relevant information related to my error message. The study I am currently working on involves a herd of 200 African buffalo caught every six months for 4 years. I am trying to use EVI and NDVI to assess seasonal variation thus I would like mean EVI and NDVI for each observation (each time each buffalo was captured). I have capture date, lat and long for each observation. However, when using ?250m_16_days_pixel_reliability? as my quality control band I keep getting the warning message:...

...g for my error code, looking at the source code. However, I have not been able to find any relevant information related to my error message. >> >> The study I am currently working on involves a herd of 200 African buffalo caught every six months for 4 years. I am trying to use EVI and NDVI to assess seasonal variation thus I would like mean EVI and NDVI for each observation (each time each buffalo was captured). I have capture date, lat and long for each observation. >> >> However, when using ?250m_16_days_pixel_reliability? as my quality control band I keep getting the...

...stion): #---> I first define the functions image.data<-function(...){ .... #-->Plot the data par(fig=c(0,1,0.25,1),new=TRUE) par(pin=c(1.5,(num.y/num.x*1.5))) image(c(1:(num.x+1)),c(1:(num.y+1)),output.matrix,xlab="X coord",ylab="Y coord",main='NDVI',las =1,cex.axis=0.7, cex.lab=0.7, col=topo.colors(255), asp=1) #---> Add a legend (colored rectangles) + legend axis (=values) par(fig=c(0,1,0.3,0.4),new=TRUE) plot.window(xlim = c(0,1), ylim = c(0,1), yaxs="i") rect((1:254)/255,0,(2:255)/255,0.1,col=topo.colors(...

...sing MODIS Tools and am having a lot of difficulty troubleshooting my code. I am a PhD student studying African buffalo in Kruger National Park, South Africa. The study I am currently working on involves a herd of 200 African buffalo caught every six months for 4 years. I am trying to use EVI and NDVI to assess seasonal variation thus I would like mean EVI and NDVI for each observation (each time each buffalo was captured). I have capture date, lat and long for each observation. However, when using ?250m_16_days_pixel_reliability? as my quality control band I keep getting the warning message:...

...having a lot of difficulty troubleshooting my code. >>> >>> I am a PhD student studying African buffalo in Kruger National Park, South Africa. The study I am currently working on involves a herd of 200 African buffalo caught every six months for 4 years. I am trying to use EVI and NDVI to assess seasonal variation thus I would like mean EVI and NDVI for each observation (each time each buffalo was captured). I have capture date, lat and long for each observation. >>> >>> However, when using ?250m_16_days_pixel_reliability? as my quality control band I keep getti...

...74, 175, 172, 171, 171, 174, 160, 175, 160, 171, 170, 175, 170, 170, 160, 160, 160, 171, 171, 171, 171, 160, 160, 160, 171, 171, 176, 171, 176, 176, 171, 176, 171, 176, 176, 176, 176, 159, 166, 159, 159, 166, 168, 169, 159, 168, 169, 166, 163, 180, 163, 165, 164, 180, 166, 166, 164, 164, 177, 166), NDVI = c(0.187, 0.2, 0.379, 0.253, 0.356, 0.341, 0.268, 0.431, 0.282, 0.181, 0.243, 0.327, 0.26, 0.232, 0.438, 0.275, 0.169, 0.288, 0.138, 0.404, 0.386, 0.194, 0.266, 0.23, 0.333, 0.234, 0.258, 0.333, 0.234, 0.096, 0.354, 0.394, 0.304, 0.162, 0.565, 0.348, 0.345, 0.226, 0.316, 0.312, 0.333, 0.28, 0.325,...

...searching for my error code, looking at the source code. However, I have not been able to find any relevant information related to my error message. > > The study I am currently working on involves a herd of 200 African buffalo caught every six months for 4 years. I am trying to use EVI and NDVI to assess seasonal variation thus I would like mean EVI and NDVI for each observation (each time each buffalo was captured). I have capture date, lat and long for each observation. > > However, when using ?250m_16_days_pixel_reliability? as my quality control band I keep getting the warning...

...5 PM, Caroline wrote: > ##MODISTools example > library(MODISTools) > library(lubridate) > setwd('~/Documents/Modis data') > > #####MODISTools with buffalo data > > ###Read in data rename for easier coding > tbdata <- read.csv('~/Desktop/All TB data for EVI, NDVI.csv') Since this dataset is only on your desktop it cannot help us reproduce your error. Can you supply a small dataset that cause the error you are talking about? One way to do this is to use supply the results of dput(tbdata) if it is small enough. If not, maybe create a subset of the dat...

...;looking at the source code. However, I have not been able to find any >relevant information related to my error message. >>> >>> The study I am currently working on involves a herd of 200 African >buffalo caught every six months for 4 years. I am trying to use EVI and >NDVI to assess seasonal variation thus I would like mean EVI and NDVI >for each observation (each time each buffalo was captured). I have >capture date, lat and long for each observation. >>> >>> However, when using ?250m_16_days_pixel_reliability? as my quality >control ban...

##MODISTools example library(MODISTools) library(lubridate) setwd('~/Documents/Modis data') #####MODISTools with buffalo data ###Read in data rename for easier coding tbdata <- read.csv('~/Desktop/All TB data for EVI, NDVI.csv') firstobs <- subset(tbdata, capture.ID == 'B1-1108') firstobs <- firstobs[,c(1,2,2,3,4)] colnames(firstobs) <- c('id', 'start.date','end.date','lat','long') ###change date format and change start date to previous 14 days firstobs$st...

...geotiffs contained within them. I tried using the untar() function to simplify this task but I am stumped by an error. I've combed the Internet for a solution without luck. The details are below, and any help in solving this matter is appreciated. > files = list.files(path = "J:/GIMMS/NDVI", pattern = "data.tif.gz", all.files = TRUE, full.names = TRUE, recursive = TRUE, ignore.case = TRUE, include.dirs = TRUE) > lapply(files, untar) Error in rawToChar(block[seq_len(ns)]) : embedded nul in string: 'II*\0Ì <\001´ \0\0`G\0\0\fn\0\0¸”\0\0d»\0\0\020â\0\0¼\b\001...

...0, + 1, 1, + 0, 0, + 0, 0, + 1, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 1, 1, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0, + 0, 0), nrow=R, ncol=J, byrow=TRUE) > > ndvi<-data.frame(x1=factor(c(4, 4, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, > 3, 4, 4, 4, 4, 3, 3, 4, 4, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, > 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, > 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3,...

...statsoft.org/v15/i05/paper). I have divided my variable in intervals: DATA$s_Shape_2 <- cut(x=DATA$Shape_inde,breaks=c(0.9999999,1.2,1.4,1.7,2.1,5)) I will use the index of the categories as the radius of the circles. Now, following the aforementioned figure 5.6, I have the following code: NDVI.breaks <- do.breaks(range(DATA.ord$NDVI), 50) xyplot(Y_Center_P~X_Center_P|Level,data=DATA.ord,col = "black",aspect = "iso", ??????????????????????????????????? fill.color = DATA.ord$c_RNDVI, cex = DATA$s_Shape_2, ??????????????????????????????????? panel = function(x, y, fi...

Dear R Helpers, Currently I'm working with the ts package of R and created a TimeSerie from pixels extracted from satellite imagery(S10 NDVI data, 10 daily composites). I'm trying to decompose this signal in different signals (seasonal and trend). When testing out the STL method is says => Only univariate timeseries are allowed, but the current Timeserie I'm using is univariate! => The problem is probably that this time s...

...ave been using plot(x,y,"") to visualise my variables. I would now like to label my points using data stored in the data.frame used for the regression analysis. For example each of my data points is made up of a field measured forest volume value and a remotely sensed vegetation estimate (NDVI). Each point is an enumeration plot and I would like to label each the points in the xy-scatterplot with their respective plot numbers. Is this possible in R, if so how do I go about doing it? Many thanks for your help Wesley Wesley Roberts MSc. Researcher: Forest Assessment (Remote Sensing &amp...