thr3ads.net - R help es - [R-es] Ayuda por favor [Jul 2014]

If this information is useful, please help other people find it:
Share via:

A. Robles R.

2014-Jul-29 06:24 UTC

[R-es] Ayuda por favor

Saludos, estoy intentando usar la función reprojectHDF() (
http://r-gis.net/?q=ModisDownload).

####
source('ModisDownload.R')
library(raster)
library(RCurl)
path<-"~/R/MODIS/Data/Test"
input<-list.files(path,"*.hdf",all.files=T,recursive=T,full.names=T)
input
outname<-paste(substr(input[1],40,45),".tif",sep='')
outname
reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin",
resample_type="NEAREST_NEIGHBOR",
proj_type="GEO",datum='WGS84')
####
Pero me resulta el siguiente mensaje de error:
####> reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin",+ resample_type="NEAREST_NEIGHBOR",
+ proj_type="GEO",datum='WGS84')
******************************************************************************
MODIS Reprojection Tool (v4.1 March 2009)
Start Time: Tue Jul 29 00:57:18 2014
------------------------------------------------------------------
Error: ReadHDFHeader : Opening Input Header File
: Unable to open
/home/rral/R/MODIS/Data/Test//home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h09v08.005.2008238073751.hdf
Fatal Error, Terminating...
[1] FALSE
####
Toda sugerencia es bienvenida.
####
PD: mis variables de entrada son:> input[1]
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h09v08.005.2008238073751.hdf"
[2]
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h09v09.005.2008238081536.hdf"
[3]
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h10v08.005.2008238053313.hdf"
[4]
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h10v09.005.2008238061408.hdf"
y> outname[1] "2000065.tif

-- 
*Atte:*   *A. Robles R.*

	[[alternative HTML version deleted]]

A. Robles R.

2014-Jul-29 06:47 UTC

head link

[R-es] Ayuda por favor

Cambie mi sintaxis de entrada a:
input<-list.files(path=".",".hdf",all.files=T,recursive=T)

y ahora mi error es:
> reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin",+
resample_type="NEAREST_NEIGHBOR",+
proj_type="GEO",datum='WGS84')******************************************************************************
MODIS Reprojection Tool (v4.1 March 2009)
Start Time:  Tue Jul 29 01:37:37 2014

------------------------------------------------------------------

Warning: gctp_call : Environmental Variable Not Found
     : MRT_DATA_DIR nor MRTDATADIR not defined
Error: GetInputGeoCorner : General Processing
     : Error converting input projection coordinates to lat/long coordinates.
Fatal Error, Terminating...
[1] FALSE



2014-07-29 1:24 GMT-05:00 A. Robles R. <albert.physik@gmail.com>:
> Saludos, estoy intentando usar la función reprojectHDF() (
> http://r-gis.net/?q=ModisDownload).
>
> ####
> source('ModisDownload.R')
> library(raster)
> library(RCurl)
> path<-"~/R/MODIS/Data/Test"
>
input<-list.files(path,"*.hdf",all.files=T,recursive=T,full.names=T)
> input
> outname<-paste(substr(input[1],40,45),".tif",sep='')
> outname
> reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin",
> resample_type="NEAREST_NEIGHBOR",
> proj_type="GEO",datum='WGS84')
> ####
> Pero me resulta el siguiente mensaje de error:
> ####
> >
reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin",
> + resample_type="NEAREST_NEIGHBOR",
> + proj_type="GEO",datum='WGS84')
>
>
******************************************************************************
> MODIS Reprojection Tool (v4.1 March 2009)
> Start Time: Tue Jul 29 00:57:18 2014
> ------------------------------------------------------------------
> Error: ReadHDFHeader : Opening Input Header File
> : Unable to open
>
/home/rral/R/MODIS/Data/Test//home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h09v08.005.2008238073751.hdf
> Fatal Error, Terminating...
> [1] FALSE
> ####
> Toda sugerencia es bienvenida.
> ####
> PD: mis variables de entrada son:
> > input
> [1]
>
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h09v08.005.2008238073751.hdf"
> [2]
>
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h09v09.005.2008238081536.hdf"
> [3]
>
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h10v08.005.2008238053313.hdf"
> [4]
>
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h10v09.005.2008238061408.hdf"
> y
> > outname
> [1] "2000065.tif
>
> --
> *Atte:*   *A. Robles R.*
>


-- 
*Atte:*   *A. Robles R.*

	[[alternative HTML version deleted]]

Carlos Ortega

2014-Jul-30 08:07 UTC

head link

[R-es] Ayuda por favor

Hola,

Creo que si sigues las recomendaciones que aparecen en esta referencia,
aumentarán mucho las posibilidades de que te podamos ayudar:

http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example

Saludos,
Carlos Ortega
www.qualityexcellence.es


2014-07-29 8:47 GMT+02:00 A. Robles R. <albert.physik@gmail.com>:
> Cambie mi sintaxis de entrada a:
>
input<-list.files(path=".",".hdf",all.files=T,recursive=T)
>
> y ahora mi error es:
>
> >
reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin",+
> resample_type="NEAREST_NEIGHBOR",+
>
proj_type="GEO",datum='WGS84')******************************************************************************
>
> MODIS Reprojection Tool (v4.1 March 2009)
> Start Time:  Tue Jul 29 01:37:37 2014
>
> ------------------------------------------------------------------
>
> Warning: gctp_call : Environmental Variable Not Found
>      : MRT_DATA_DIR nor MRTDATADIR not defined
> Error: GetInputGeoCorner : General Processing
>      : Error converting input projection coordinates to lat/long
> coordinates.
> Fatal Error, Terminating...
> [1] FALSE
>
>
>
> 2014-07-29 1:24 GMT-05:00 A. Robles R. <albert.physik@gmail.com>:
>
> > Saludos, estoy intentando usar la función reprojectHDF() (
> > http://r-gis.net/?q=ModisDownload).
> >
> > ####
> > source('ModisDownload.R')
> > library(raster)
> > library(RCurl)
> > path<-"~/R/MODIS/Data/Test"
> >
input<-list.files(path,"*.hdf",all.files=T,recursive=T,full.names=T)
> > input
> >
outname<-paste(substr(input[1],40,45),".tif",sep='')
> > outname
> >
reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin",
> > resample_type="NEAREST_NEIGHBOR",
> > proj_type="GEO",datum='WGS84')
> > ####
> > Pero me resulta el siguiente mensaje de error:
> > ####
> > >
reprojectHDF(hdfName=input,filename=outname,MRTpath="~/MRT/bin",
> > + resample_type="NEAREST_NEIGHBOR",
> > + proj_type="GEO",datum='WGS84')
> >
> >
>
******************************************************************************
> > MODIS Reprojection Tool (v4.1 March 2009)
> > Start Time: Tue Jul 29 00:57:18 2014
> > ------------------------------------------------------------------
> > Error: ReadHDFHeader : Opening Input Header File
> > : Unable to open
> >
>
/home/rral/R/MODIS/Data/Test//home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h09v08.005.2008238073751.hdf
> > Fatal Error, Terminating...
> > [1] FALSE
> > ####
> > Toda sugerencia es bienvenida.
> > ####
> > PD: mis variables de entrada son:
> > > input
> > [1]
> >
>
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h09v08.005.2008238073751.hdf"
> > [2]
> >
>
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h09v09.005.2008238081536.hdf"
> > [3]
> >
>
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h10v08.005.2008238053313.hdf"
> > [4]
> >
>
"/home/rral/R/MODIS/Data/Test/MOD13Q1.A2000065.h10v09.005.2008238061408.hdf"
> > y
> > > outname
> > [1] "2000065.tif
> >
> > --
> > *Atte:*   *A. Robles R.*
> >
>
>
>
> --
> *Atte:*   *A. Robles R.*
>
>         [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es@r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
>

-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]

R help es - Jul 2014 - Ayuda por favor

[R-es] Ayuda por favor

[R-es] Ayuda por favor

[R-es] Ayuda por favor

Possibly Parallel Threads