[R] MODISTools help - with reproducible examples
David Winsemius
dwinsemius at comcast.net
Fri Jun 23 05:25:05 CEST 2017
> On Jun 22, 2017, at 4:44 PM, Caroline <gliddeca at science.oregonstate.edu> wrote:
>
> I am using the R-package MODISTools (different than MODIS) and am having a lot of difficulty troubleshooting my code. I have spent awhile going through MODISTools tutorials, 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:
>
> Warning in MODISSummaries(LoadDat = period, FileSep = ",", Product = "MOD13Q1", :
> Only single data point that passed the quality screen: cannot summarise
>
> When using ‘250m_16_days_VI_Quality’ as my quality control band I keep getting the warning message:
>
> Error in QualityCheck(Data = band.time.series, QualityScores = QA.time.series, :
> QualityScores not all in range of MOD13Q1's QC: 0-3
>
> I seem to get this message with all subsets of my data (I have tried running all of my data at once and then just one data point at a time). I have also tried using wider date ranges as well as wider size ranges (in case the pixel reliability is poor within a certain area or time frame) but still get the same messages.
>
> I have attached the data file I have been using as well as my code (subsetted to just one animal so the run time is faster)
>
> Has anyone seen this error message before/have any suggestions for why I may be getting these errors?
I don't and you have not provided the accepted means to investigate which would be code and data . Have you contacted the package maintainer?
--
David Winsemius
Alameda, CA, USA
More information about the R-help
mailing list