[GRASS-user] FW: Temporal + export

Matias_Alejandro_Cas · May 5, 2016, 3:41pm

Hi everyone!

I never used a list… I think I made a mistake when sent the message yesterday…

I was construct a strds in grass 7.0.3, and I want to export it for use in R to do some statistics not included with tgrass… But when I do the export operation, I receive this message:

Usising the GUI:

t.rast.export --overwrite --verbose input=LST2@MOD11A2 output=C:\SIG\Base de datos area de estudio\Datasets_R\LSR_R compression=gzip
Traceback (most recent call last):
File “C:\OSGEO4~1\apps\grass\grass-7.0.3/scripts/t.rast.ex
port.py”, line 87, in
main()
File “C:\OSGEO4~1\apps\grass\grass-7.0.3/scripts/t.rast.ex
port.py”, line 82, in main
_input, output, compression, directory, where, _format,
“strds”)
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass
temporal\stds_export.py”, line 268, in export_stds
new_cwd = tempfile.mkdtemp(dir=directory)
File “C:\OSGEO4~1\apps\Python27\lib\tempfile.py”, line
329, in mkdtemp
_os.mkdir(file, 0700)
WindowsError: [Error 3] El sistema no puede encontrar la
ruta especificada: ‘/tmp\tmppi2oh1’
(Wed May 04 16:09:44 2016) Comando finalizado. (1 segundos)

In python (very difficult):

input=“LST2@MOD11A2”
output=“LST_R.tar.gz”
compression=“gzip”
directory=“C:/SIG/Base de datos area de estudio/Datasets_R/”
where=None
format_=“GTiff”
type_=“strds”
tgis.export_stds(input, output, compression, directory, where, format_, type_)

Traceback (most recent call last):
File “”, line 1, in
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py”, line 298, in export_stds
rows, tar, list_file, new_cwd, fs, format_)
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py”, line 74, in _export_raster_maps_as_gdal
nodata = max_val + 1
TypeError: unsupported operand type(s) for +: ‘NoneType’ and ‘int’

Some idea? I was try to set the nodata values to 0 using gdalwarp but this doesn’t work apparently…

Thanks you!!

Matías

neteler · May 5, 2016, 4:39pm

Hi Matías,

On Thu, May 5, 2016 at 5:41 PM, Matías Alejandro Castillo Moine
<matialecastillo@hotmail.com> wrote:

Hi everyone!

I never used a list... I think I made a mistake when sent the message
yesterday...

now it reached the list

I was construct a strds in grass 7.0.3,

(Would it be possible for you to update to 7.0.4? Just to be sure that
the problem is yes/no solved?

and I want to export it for use in R
to do some statistics not included with tgrass.. But when I do the export
operation, I receive this message:

Usising the GUI:
t.rast.export --overwrite --verbose input=LST2@MOD11A2 output=C:\SIG\Base de
datos area de estudio\Datasets_R\LSR_R compression=gzip

Hint: You are using white space in the path ("Base de datos area de
estudio\"). Maybe it doesn't matter but just to exclude this potential
issue, please try with a different path name.

...

    new_cwd = tempfile.mkdtemp(dir=directory)
  File "C:\OSGEO4~1\apps\Python27\lib\tempfile.py", line
329, in mkdtemp
    _os.mkdir(file, 0700)
WindowsError: [Error 3] El sistema no puede encontrar la
ruta especificada: '/tmp\\tmppi2oh1'

This looks like a bug: it uses /tmp/ which does not exist on Windows.
The current code is:

/lib/python/temporal/stds_export.py, line 274
    # Create the temporary directory and jump into it
    new_cwd = tempfile.mkdtemp(dir=directory)
    os.chdir(new_cwd)

but AFAIK it should use this portable function:
lib/python/script/core.py:def tempdir():

@devs: opinions?

Concerning the other issue:

In python (very difficult):

input="LST2@MOD11A2"
output="LST_R.tar.gz"
compression="gzip"
directory="C:/SIG/Base de datos area de estudio/Datasets_R/"
where=None
format_="GTiff"
type_="strds"
tgis.export_stds(input, output, compression, directory, where, format_,
type_)

Traceback (most recent call last):

  File "<input>", line 1, in <module>
  File
"C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py",
line 298, in export_stds
    rows, tar, list_file, new_cwd, fs, format_)
  File
"C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py",
line 74, in _export_raster_maps_as_gdal
    nodata = max_val + 1
TypeError: unsupported operand type(s) for +: 'NoneType' and 'int'

...here I have no suggestion. Anyone?

Markus

veroandreo · May 5, 2016, 5:07pm

Hola Mati!!!

I’m so glad you are using grass and writing to the list!

Markus won me, but I had the same pieces of advice.

First idea/advice: do not use folder names (nor file names) with spaces, second: maybe you wanna update to 7.0.4, the new stable release.

I’ve tried here (grass71, fedora23 box) with:

t.rast.export input=cla@clorofila output=cla.gzip compression=gzip directory=/home/veroandreo where=“start_time >= ‘2013-01-01’”

Checking GDAL data type and nodata value…
100%
Using GDAL data type
Input raster map contains cells with NULL-value (no-data). The value -nan
will be used to represent no-data values in the input map. You can specify
a nodata value with the nodata option.
Exporting raster data to GTiff format…
100%
r.out.gdal complete. File <A20130012013008.L3m_8D_CHL_chlor_a_4km_arg.tif>
created.

…
tar xvf cla.gzip

all is there and looks fine…

The nodata option comes from r.out.gdal [1]. t.rast.export uses it to export maps in the strds. Don’t know about gdalwarp, but in any case, if you will then import data into R for other analyses, you can change the value there.

Un abrazo,
Vero

[1] https://grass.osgeo.org/grass70/manuals/r.out.gdal.html

···

2016-05-05 12:41 GMT-03:00 Matías Alejandro Castillo Moine <matialecastillo@hotmail.com>:

Hi everyone!

I never used a list… I think I made a mistake when sent the message yesterday…

I was construct a strds in grass 7.0.3, and I want to export it for use in R to do some statistics not included with tgrass… But when I do the export operation, I receive this message:

Usising the GUI:

t.rast.export --overwrite --verbose input=LST2@MOD11A2 output=C:\SIG\Base de datos area de estudio\Datasets_R\LSR_R compression=gzip
Traceback (most recent call last):
File “C:\OSGEO4~1\apps\grass\grass-7.0.3/scripts/t.rast.ex
port.py”, line 87, in
main()
File “C:\OSGEO4~1\apps\grass\grass-7.0.3/scripts/t.rast.ex
port.py”, line 82, in main
_input, output, compression, directory, where, _format,
“strds”)
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass
temporal\stds_export.py”, line 268, in export_stds
new_cwd = tempfile.mkdtemp(dir=directory)
File “C:\OSGEO4~1\apps\Python27\lib\tempfile.py”, line
329, in mkdtemp
_os.mkdir(file, 0700)
WindowsError: [Error 3] El sistema no puede encontrar la
ruta especificada: ‘/tmp\tmppi2oh1’
(Wed May 04 16:09:44 2016) Comando finalizado. (1 segundos)

In python (very difficult):

input=“LST2@MOD11A2”
output=“LST_R.tar.gz”
compression=“gzip”
directory=“C:/SIG/Base de datos area de estudio/Datasets_R/”
where=None
format_=“GTiff”
type_=“strds”
tgis.export_stds(input, output, compression, directory, where, format_, type_)

Traceback (most recent call last):
File “”, line 1, in
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py”, line 298, in export_stds
rows, tar, list_file, new_cwd, fs, format_)
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py”, line 74, in _export_raster_maps_as_gdal
nodata = max_val + 1
TypeError: unsupported operand type(s) for +: ‘NoneType’ and ‘int’

Some idea? I was try to set the nodata values to 0 using gdalwarp but this doesn’t work apparently…

Thanks you!!

Matías

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Matias_Alejandro_Cas · May 6, 2016, 9:24pm

Hola Verooooo!
Thanks you for answer!
I’m working to solve it. Apparently all my imported rasters are zero cell values. There is something wrong with my gdalwarp command… I’ll tell you later…
Thanks!!!

Date: Thu, 5 May 2016 14:07:56 -0300
Subject: Re: [GRASS-user] FW: Temporal + export
From: veroandreo@gmail.com
To: matialecastillo@hotmail.com
CC: grass-user@lists.osgeo.org

Hola Mati!!!

I’m so glad you are using grass and writing to the list!

Markus won me, but I had the same pieces of advice.

First idea/advice: do not use folder names (nor file names) with spaces, second: maybe you wanna update to 7.0.4, the new stable release.

I’ve tried here (grass71, fedora23 box) with:

t.rast.export input=cla@clorofila output=cla.gzip compression=gzip directory=/home/veroandreo where=“start_time >= ‘2013-01-01’”

Checking GDAL data type and nodata value…
100%
Using GDAL data type
Input raster map contains cells with NULL-value (no-data). The value -nan
will be used to represent no-data values in the input map. You can specify
a nodata value with the nodata option.
Exporting raster data to GTiff format…
100%
r.out.gdal complete. File <A20130012013008.L3m_8D_CHL_chlor_a_4km_arg.tif>
created.

…
tar xvf cla.gzip

all is there and looks fine…

The nodata option comes from r.out.gdal [1]. t.rast.export uses it to export maps in the strds. Don’t know about gdalwarp, but in any case, if you will then import data into R for other analyses, you can change the value there.

Un abrazo,
Vero

[1] https://grass.osgeo.org/grass70/manuals/r.out.gdal.html

2016-05-05 12:41 GMT-03:00 Matías Alejandro Castillo Moine <matialecastillo@hotmail.com>:

Hi everyone!

I never used a list… I think I made a mistake when sent the message yesterday…

I was construct a strds in grass 7.0.3, and I want to export it for use in R to do some statistics not included with tgrass… But when I do the export operation, I receive this message:

Usising the GUI:

t.rast.export --overwrite --verbose input=LST2@MOD11A2 output=C:\SIG\Base de datos area de estudio\Datasets_R\LSR_R compression=gzip
Traceback (most recent call last):
File “C:\OSGEO4~1\apps\grass\grass-7.0.3/scripts/t.rast.ex
port.py”, line 87, in
main()
File “C:\OSGEO4~1\apps\grass\grass-7.0.3/scripts/t.rast.ex
port.py”, line 82, in main
_input, output, compression, directory, where, _format,
“strds”)
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass
temporal\stds_export.py”, line 268, in export_stds
new_cwd = tempfile.mkdtemp(dir=directory)
File “C:\OSGEO4~1\apps\Python27\lib\tempfile.py”, line
329, in mkdtemp
_os.mkdir(file, 0700)
WindowsError: [Error 3] El sistema no puede encontrar la
ruta especificada: ‘/tmp\tmppi2oh1’
(Wed May 04 16:09:44 2016) Comando finalizado. (1 segundos)

In python (very difficult):

input=“LST2@MOD11A2”
output=“LST_R.tar.gz”
compression=“gzip”
directory=“C:/SIG/Base de datos area de estudio/Datasets_R/”
where=None
format_=“GTiff”
type_=“strds”
tgis.export_stds(input, output, compression, directory, where, format_, type_)

Traceback (most recent call last):
File “”, line 1, in
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py”, line 298, in export_stds
rows, tar, list_file, new_cwd, fs, format_)
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py”, line 74, in _export_raster_maps_as_gdal
nodata = max_val + 1
TypeError: unsupported operand type(s) for +: ‘NoneType’ and ‘int’

Some idea? I was try to set the nodata values to 0 using gdalwarp but this doesn’t work apparently…

Thanks you!!

Matías

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

veroandreo · May 7, 2016, 2:24am

Hola Mati!

Maybe if you tell us more about your data, workflow and the commands you are using, we could help you. What do you use gdalwarp for?

Cheers,
Vero

···

2016-05-06 18:24 GMT-03:00 Matías Alejandro Castillo Moine <matialecastillo@hotmail.com>:

Hola Verooooo!
Thanks you for answer!
I’m working to solve it. Apparently all my imported rasters are zero cell values. There is something wrong with my gdalwarp command… I’ll tell you later…
Thanks!!!

Date: Thu, 5 May 2016 14:07:56 -0300
Subject: Re: [GRASS-user] FW: Temporal + export
From: veroandreo@gmail.com
To: matialecastillo@hotmail.com
CC: grass-user@lists.osgeo.org

Hola Mati!!!

I’m so glad you are using grass and writing to the list!

Markus won me, but I had the same pieces of advice.

First idea/advice: do not use folder names (nor file names) with spaces, second: maybe you wanna update to 7.0.4, the new stable release.

I’ve tried here (grass71, fedora23 box) with:

t.rast.export input=cla@clorofila output=cla.gzip compression=gzip directory=/home/veroandreo where=“start_time >= ‘2013-01-01’”

Checking GDAL data type and nodata value…
100%
Using GDAL data type
Input raster map contains cells with NULL-value (no-data). The value -nan
will be used to represent no-data values in the input map. You can specify
a nodata value with the nodata option.
Exporting raster data to GTiff format…
100%
r.out.gdal complete. File <A20130012013008.L3m_8D_CHL_chlor_a_4km_arg.tif>
created.

…
tar xvf cla.gzip

all is there and looks fine…

The nodata option comes from r.out.gdal [1]. t.rast.export uses it to export maps in the strds. Don’t know about gdalwarp, but in any case, if you will then import data into R for other analyses, you can change the value there.

Un abrazo,
Vero

[1] https://grass.osgeo.org/grass70/manuals/r.out.gdal.html

2016-05-05 12:41 GMT-03:00 Matías Alejandro Castillo Moine <matialecastillo@hotmail.com>:

Hi everyone!

I never used a list… I think I made a mistake when sent the message yesterday…

I was construct a strds in grass 7.0.3, and I want to export it for use in R to do some statistics not included with tgrass… But when I do the export operation, I receive this message:

Usising the GUI:

t.rast.export --overwrite --verbose input=LST2@MOD11A2 output=C:\SIG\Base de datos area de estudio\Datasets_R\LSR_R compression=gzip
Traceback (most recent call last):
File “C:\OSGEO4~1\apps\grass\grass-7.0.3/scripts/t.rast.ex
port.py”, line 87, in
main()
File “C:\OSGEO4~1\apps\grass\grass-7.0.3/scripts/t.rast.ex
port.py”, line 82, in main
_input, output, compression, directory, where, _format,
“strds”)
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass
temporal\stds_export.py”, line 268, in export_stds
new_cwd = tempfile.mkdtemp(dir=directory)
File “C:\OSGEO4~1\apps\Python27\lib\tempfile.py”, line
329, in mkdtemp
_os.mkdir(file, 0700)
WindowsError: [Error 3] El sistema no puede encontrar la
ruta especificada: ‘/tmp\tmppi2oh1’
(Wed May 04 16:09:44 2016) Comando finalizado. (1 segundos)

In python (very difficult):

input=“LST2@MOD11A2”
output=“LST_R.tar.gz”
compression=“gzip”
directory=“C:/SIG/Base de datos area de estudio/Datasets_R/”
where=None
format_=“GTiff”
type_=“strds”
tgis.export_stds(input, output, compression, directory, where, format_, type_)

Traceback (most recent call last):
File “”, line 1, in
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py”, line 298, in export_stds
rows, tar, list_file, new_cwd, fs, format_)
File “C:\OSGEO4~1\apps\grass\grass-7.0.3\etc\python\grass\temporal\stds_export.py”, line 74, in _export_raster_maps_as_gdal
nodata = max_val + 1
TypeError: unsupported operand type(s) for +: ‘NoneType’ and ‘int’

Some idea? I was try to set the nodata values to 0 using gdalwarp but this doesn’t work apparently…

Thanks you!!

Matías

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

neteler · May 7, 2016, 12:09pm

On May 6, 2016 11:24 PM, “Matías Alejandro Castillo Moine” <matialecastillo@hotmail.com> wrote:

Hola Verooooo!
Thanks you for answer!
I’m working to solve it. Apparently all my imported rasters are zero cell values.
… maybe the data are there but the computational region wrong?

Markus

Matias_Alejandro_Cas · May 13, 2016, 4:17pm

Hi! I´m again!
The problem whith galwarp is solved (spatial reference system was wrong).
But now, I have my time serie registered in Grass and when I export it, I can only export this using the AAgrid or pack format, and R needs the GTiff formatt into the tar.gz archive. I believe that I´m using correctly the command

t.rast.export --overwrite input=LST2@MOD11A2 output=MOD11A2.tar.gz directory=C:\SIG\Grass compression=gzip
WARNING: El tipo de datos GDAL seleccionado no cubre el rango de datos
WARNING: Tipo de datos GDAL: Byte, rango: 0 - 255
WARNING: Rango a exportar: 14837 - 15732
ERROR: La exportación raster puede resultar en una pérdida completa de datos, cancelando.
Traceback (most recent call last):
File “C:\OSGEO4~1\apps\grass\grass-7.0.4/scripts/t.rast.ex
port.py”, line 87, in
main()
File “C:\OSGEO4~1\apps\grass\grass-7.0.4/scripts/t.rast.ex
port.py”, line 82, in main
_input, output, compression, directory, where, format,
“strds”)
File “C:\OSGEO4~1\apps\grass\grass-7.0.4\etc\python\grass
temporal\stds_export.py”, line 298, in export_stds
rows, tar, list_file, new_cwd, fs, format)
File “C:\OSGEO4~1\apps\grass\grass-7.0.4\etc\python\grass
temporal\stds_export.py”, line 97, in
_export_raster_maps_as_gdal
shutil.rmtree(new_cwd)
File “C:\OSGEO4~1\apps\Python27\lib\shutil.py”, line 252,
in rmtree
onerror(os.remove, fullname, sys.exc_info())
File “C:\OSGEO4~1\apps\Python27\lib\shutil.py”, line 250,
in rmtree
os.remove(fullname)
WindowsError: [Error 32] El proceso no tiene acceso al
archivo porque está siendo utilizado por otro proceso:
‘C:\SIG\Grass\tmp883adw\archive’
(Fri May 13 13:09:21 2016) Comando finalizado. (2 segundos)
(Fri May 13 13:14:02 2016)

Some new idea?

Thanks!!!

Date: Sat, 7 May 2016 14:09:33 +0200
Subject: Re: [GRASS-user] FW: Temporal + export
From: neteler@osgeo.org
To: matialecastillo@hotmail.com
CC: veroandreo@gmail.com; grass-user@lists.osgeo.org

On May 6, 2016 11:24 PM, “Matías Alejandro Castillo Moine” <matialecastillo@hotmail.com> wrote:

Hola Verooooo!
Thanks you for answer!
I’m working to solve it. Apparently all my imported rasters are zero cell values.
… maybe the data are there but the computational region wrong?

Markus

veroandreo · May 13, 2016, 6:07pm

Hey, Mati! Hola!

(Hint: you may want to set language to English so more people can understand the error messages and help)

El may 13, 2016 1:17 PM, “Matías Alejandro Castillo Moine” <matialecastillo@hotmail.com> escribió:

Hi! I´m again!
The problem whith galwarp is solved (spatial reference system was wrong).
But now, I have my time serie registered in Grass and when I export it, I can only export this using the AAgrid or pack format, and R needs the GTiff formatt into the tar.gz archive. I believe that I´m using correctly the command

t.rast.export --overwrite input=LST2@MOD11A2 output=MOD11A2.tar.gz directory=C:\SIG\Grass compression=gzip

If this is the command, you are not setting format there, so by default you should be exporting GTiff.

In the manual page for 7.1 also says that if you use gzip compression, you should add .tar.gzip to the output file name. This latter is not specified in 70 manual but I guess it applies here, too [0].

WARNING: El tipo de datos GDAL seleccionado no cubre el rango de datos
WARNING: Tipo de datos GDAL: Byte, rango: 0 - 255
WARNING: Rango a exportar: 14837 - 15732

What I can see there is that you are exporting your maps as byte (0-255) when your data is not in that range (WARNING: Rango a exportar: 14837 - 15732). In t.rast.export for grasss 7.1 there’s a new parameter called type which allows to set data type and would allow you to change that [1]. Dunno which is default type for AAgrid or pack (not mentioned in the manual page), but for GTiff you have several options available (in that case you could try forcing type to UInt16, which would cover your data range). Anyway, this seems odd since you are using 7.0.4 and not 7.1, right?

@Luca, might the error and/or message be related to [2]?

ERROR: La exportación raster puede resultar en una pérdida completa de datos, cancelando.
Traceback (most recent call last):
File “C:\OSGEO4~1\apps\grass\grass-7.0.4/scripts/t.rast.ex
port.py”, line 87, in
main()
File “C:\OSGEO4~1\apps\grass\grass-7.0.4/scripts/t.rast.ex
port.py”, line 82, in main
_input, output, compression, directory, where, format,
“strds”)
File “C:\OSGEO4~1\apps\grass\grass-7.0.4\etc\python\grass
temporal\stds_export.py”, line 298, in export_stds
rows, tar, list_file, new_cwd, fs, format)
File “C:\OSGEO4~1\apps\grass\grass-7.0.4\etc\python\grass
temporal\stds_export.py”, line 97, in
_export_raster_maps_as_gdal
shutil.rmtree(new_cwd)
File “C:\OSGEO4~1\apps\Python27\lib\shutil.py”, line 252,
in rmtree
onerror(os.remove, fullname, sys.exc_info())
File “C:\OSGEO4~1\apps\Python27\lib\shutil.py”, line 250,
in rmtree
os.remove(fullname)
WindowsError: [Error 32] El proceso no tiene acceso al
archivo porque está siendo utilizado por otro proceso:
‘C:\SIG\Grass\tmp883adw\archive’
(Fri May 13 13:09:21 2016) Comando finalizado. (2 segundos)
(Fri May 13 13:14:02 2016)

Some new idea?

Thanks!!!

Does r.out.gdal for only one map works? Something like

r.out.gdal input=LST output=LST.tif format=GTiff type=UInt16

Let us know

Vero

[0] https://trac.osgeo.org/grass/ticket/2416
[1] https://grass.osgeo.org/grass71/manuals/t.rast.export.html
[2] https://trac.osgeo.org/grass/changeset/67961

lucadelu · June 9, 2016, 1:19pm

On 13 May 2016 at 20:07, Veronica Andreo <veroandreo@gmail.com> wrote:

Hey, Mati! Hola!

Hi Vero, Matias

@Luca, might the error and/or message be related to [2]?

I don't think so

ERROR: La exportación raster puede resultar en una pérdida completa de
datos, cancelando.
Traceback (most recent call last):
  File "C:\OSGEO4~1\apps\grass\grass-7.0.4/scripts/t.rast.ex
port.py", line 87, in <module>
    main()
  File "C:\OSGEO4~1\apps\grass\grass-7.0.4/scripts/t.rast.ex
port.py", line 82, in main
    _input, output, compression, directory, where, _format,
"strds")
  File "C:\OSGEO4~1\apps\grass\grass-7.0.4\etc\python\grass\
temporal\stds_export.py", line 298, in export_stds
    rows, tar, list_file, new_cwd, fs, format_)
  File "C:\OSGEO4~1\apps\grass\grass-7.0.4\etc\python\grass\
temporal\stds_export.py", line 97, in
_export_raster_maps_as_gdal
    shutil.rmtree(new_cwd)
  File "C:\OSGEO4~1\apps\Python27\lib\shutil.py", line 252,
in rmtree
    onerror(os.remove, fullname, sys.exc_info())
  File "C:\OSGEO4~1\apps\Python27\lib\shutil.py", line 250,
in rmtree
    os.remove(fullname)
WindowsError: [Error 32] El proceso no tiene acceso al
archivo porque está siendo utilizado por otro proceso:
'C:\\SIG\\Grass\\tmp883adw\\archive'
(Fri May 13 13:09:21 2016) Comando finalizado. (2 segundos)
(Fri May 13 13:14:02 2016)

it seems related to r.out.gdal, Matias could you try again the command?

Some new idea?

Thanks!!!

Does r.out.gdal for only one map works? Something like

r.out.gdal input=LST output=LST.tif format=GTiff type=UInt16

Matias is this working for you?

Let us know

Vero

--
ciao
Luca

www.lucadelu.org

Matias_Alejandro_Cas · July 5, 2016, 9:00pm

Hi everyone!
I solved this generating a .tif file with the target SRC, and then reading their SCR from .tif file when the grass location is created. Then I import all the .tif files (-o flag is not activated) and finally t.rast.export works correctly.

Thanks you!!!

From: lucadeluge@gmail.com
Date: Thu, 9 Jun 2016 15:19:20 +0200
Subject: Re: [GRASS-user] FW: Temporal + export
To: veroandreo@gmail.com
CC: matialecastillo@hotmail.com; grass-user@lists.osgeo.org

On 13 May 2016 at 20:07, Veronica Andreo veroandreo@gmail.com wrote:

Hey, Mati! Hola!

Hi Vero, Matias

@Luca, might the error and/or message be related to [2]?

I don’t think so

ERROR: La exportación raster puede resultar en una pérdida completa de
datos, cancelando.
Traceback (most recent call last):
File “C:\OSGEO4~1\apps\grass\grass-7.0.4/scripts/t.rast.ex
port.py”, line 87, in
main()
File “C:\OSGEO4~1\apps\grass\grass-7.0.4/scripts/t.rast.ex
port.py”, line 82, in main
_input, output, compression, directory, where, format,
“strds”)
File “C:\OSGEO4~1\apps\grass\grass-7.0.4\etc\python\grass
temporal\stds_export.py”, line 298, in export_stds
rows, tar, list_file, new_cwd, fs, format)
File “C:\OSGEO4~1\apps\grass\grass-7.0.4\etc\python\grass
temporal\stds_export.py”, line 97, in
_export_raster_maps_as_gdal
shutil.rmtree(new_cwd)
File “C:\OSGEO4~1\apps\Python27\lib\shutil.py”, line 252,
in rmtree
onerror(os.remove, fullname, sys.exc_info())
File “C:\OSGEO4~1\apps\Python27\lib\shutil.py”, line 250,
in rmtree
os.remove(fullname)
WindowsError: [Error 32] El proceso no tiene acceso al
archivo porque está siendo utilizado por otro proceso:
‘C:\SIG\Grass\tmp883adw\archive’
(Fri May 13 13:09:21 2016) Comando finalizado. (2 segundos)
(Fri May 13 13:14:02 2016)

it seems related to r.out.gdal, Matias could you try again the command?

Some new idea?

Thanks!!!

Does r.out.gdal for only one map works? Something like

r.out.gdal input=LST output=LST.tif format=GTiff type=UInt16

Matias is this working for you?

Let us know

Vero

–
ciao
Luca

www.lucadelu.org