writeRaster to ascii file (.asc)

classic Classic list List threaded Threaded
6 messages Options
Reply | Threaded
Open this post in threaded view
|

writeRaster to ascii file (.asc)

Manuel Spínola
Dear list members,

I am trying to export a raster file to a ascii file (.asc) but the min and
max values are not the same. Any idea why?

> bio1res
class       : RasterLayer
dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
resolution  : 0.008333333, 0.008333333  (x, y)
extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin, ymax)
coord. ref. : NA
data source : in memory
names       : bio1_23
values      : 47, 272  (min, max)

> writeRaster(bio1res, filename="temp.asc", format="ascii",overwrite=TRUE)
class       : RasterLayer
dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
resolution  : 0.008333333, 0.008333333  (x, y)
extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin, ymax)
coord. ref. : NA
data source : /Users/manuelspinola/Dropbox/Distribución de
especies/World_Clim_Costa_Rica/temp.asc
names       : temp
values      : -2147483648, 2147483647  (min, max)


Best,

Manuel



--
*Manuel Spínola, Ph.D.*
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
[hidden email]
[hidden email]
Teléfono: (506) 2277-3598
Fax: (506) 2237-7036
Personal website: Lobito de río <https://sites.google.com/site/lobitoderio/>
Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>

        [[alternative HTML version deleted]]


_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

Re: writeRaster to ascii file (.asc)

Mauricio Zambrano-Bigiarini
On 10/07/13 18:38, Manuel Spínola wrote:
> Dear list members,
>
> I am trying to export a raster file to a ascii file (.asc) but the min and
> max values are not the same. Any idea why?

Hi Manuel,

Could you provide your sessionInfo() ?

>
>> bio1res
> class       : RasterLayer
> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
> resolution  : 0.008333333, 0.008333333  (x, y)
> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin, ymax)
> coord. ref. : NA
> data source : in memory
> names       : bio1_23
> values      : 47, 272  (min, max)
>
>> writeRaster(bio1res, filename="temp.asc", format="ascii",overwrite=TRUE)
> class       : RasterLayer
> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
> resolution  : 0.008333333, 0.008333333  (x, y)
> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin, ymax)
> coord. ref. : NA
> data source : /Users/manuelspinola/Dropbox/Distribucio�n de
> especies/World_Clim_Costa_Rica/temp.asc
> names       : temp
> values      : -2147483648, 2147483647  (min, max)
Unfortunately, your example is not reproducible, so I tried with a
different example and the values in the original and the written ascii
file were the same (rgdal_0.8-10, raster_2.1-48, sp_1.0-11):

----------  START --------------
library(raster)
Loading required package: sp

# Creation of a raster file for testing
r <- raster(system.file("external/test.grd", package="raster"))

r
#class       : RasterLayer
#dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
#resolution  : 40, 40  (x, y)
#extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
#coord. ref. : +init=epsg:28992
+towgs84=565.237,50.0087,465.658,-0.406857,0.350733,-1.87035,4.0812
+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889
+k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs
#data source : /usr/lib64/R/library/raster/external/test.grd
#names       : test
#values      : 128.434, 1805.78  (min, max)

writeRaster(r, filename="temp.asc", format="ascii",overwrite=TRUE)
#rgdal: version: 0.8-10, (SVN revision 478)
#Geospatial Data Abstraction Library extensions to R successfully loaded
#Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
#Path to GDAL shared files: /usr/share/gdal
#Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
#Path to PROJ.4 shared files: (autodetected)
#class       : RasterLayer
#dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
#resolution  : 40, 40  (x, y)
#extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
#coord. ref. : NA
#data source : /home/zambrhe/temp.asc
#names       : temp

# reading the .asc file form disk
r2 <- raster("temp.asc")

r2
#class       : RasterLayer
#dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
#resolution  : 40, 40  (x, y)
#extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
#coord. ref. : NA
#data source : /home/zambrhe/temp.asc
#names       : temp

# Summary of the ascii file in the hard disk
#summary(r2)
#             temp
#Min.     128.4340
#1st Qu.  293.2325
#Median   371.4120
#3rd Qu.  499.8195
#Max.    1805.7800
#NA's    6097.0000


# Summary of the original raster file
summary(r)
#             test
#Min.     128.4340
#1st Qu.  293.2325
#Median   371.4120
#3rd Qu.  499.8195
#Max.    1805.7800
#NA's    6097.0000
----------  END --------------

Kind regards,

Mauricio Zambrano-Bigiarini, Ph.D

--
=================================================
Water Resources Unit
Institute for Environment and Sustainability (IES)
Joint Research Centre (JRC), European Commission
TP 261, Via Enrico Fermi 2749, 21027 Ispra (VA), IT
Work Phone : +39 0332 789588
webinfo    : http://floods.jrc.ec.europa.eu/
=================================================
DISCLAIMER:
"The views expressed are purely those of the writer
and may not in any circumstances be regarded as sta-
ting an official position of the European Commission"
=================================================
"The journey is the reward"  (Steve Jobs)


_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

Re: writeRaster to ascii file (.asc)

Manuel Spínola
Thank you Mauricio.

> sessionInfo()
R version 3.0.1 (2013-05-16)
Platform: x86_64-apple-darwin10.8.0 (64-bit)

locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8

attached base packages:
[1] grid      stats     graphics  grDevices utils     datasets  methods
base

other attached packages:
[1] SDMTools_1.1-13 dismo_0.8-11    raster_2.1-37   rgdal_0.8-10
 maptools_0.8-25
[6] lattice_0.20-15 sp_1.0-11       foreign_0.8-54

loaded via a namespace (and not attached):
[1] tools_3.0.1


2013/7/11 Mauricio Zambrano-Bigiarini <[hidden email]>

> On 10/07/13 18:38, Manuel Spínola wrote:
>
>> Dear list members,
>>
>> I am trying to export a raster file to a ascii file (.asc) but the min and
>> max values are not the same. Any idea why?
>>
>
> Hi Manuel,
>
> Could you provide your sessionInfo() ?
>
>>
>>  bio1res
>>>
>> class       : RasterLayer
>> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
>> resolution  : 0.008333333, 0.008333333  (x, y)
>> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin,
>> ymax)
>> coord. ref. : NA
>> data source : in memory
>> names       : bio1_23
>> values      : 47, 272  (min, max)
>>
>>  writeRaster(bio1res, filename="temp.asc", format="ascii",overwrite=TRUE)
>>>
>> class       : RasterLayer
>> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
>> resolution  : 0.008333333, 0.008333333  (x, y)
>> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin,
>> ymax)
>> coord. ref. : NA
>> data source : /Users/manuelspinola/Dropbox/**DistribucioÌ n de
>>
>> especies/World_Clim_Costa_**Rica/temp.asc
>> names       : temp
>> values      : -2147483648, 2147483647  (min, max)
>>
>
> Unfortunately, your example is not reproducible, so I tried with a
> different example and the values in the original and the written ascii file
> were the same (rgdal_0.8-10, raster_2.1-48, sp_1.0-11):
>
> ----------  START --------------
> library(raster)
> Loading required package: sp
>
> # Creation of a raster file for testing
> r <- raster(system.file("external/**test.grd", package="raster"))
>
> r
> #class       : RasterLayer
> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
> #resolution  : 40, 40  (x, y)
> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
> #coord. ref. : +init=epsg:28992 +towgs84=565.237,50.0087,465.**
> 658,-0.406857,0.350733,-1.**87035,4.0812 +proj=sterea
> +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000
> +y_0=463000 +ellps=bessel +units=m +no_defs
> #data source : /usr/lib64/R/library/raster/**external/test.grd
> #names       : test
> #values      : 128.434, 1805.78  (min, max)
>
> writeRaster(r, filename="temp.asc", format="ascii",overwrite=TRUE)
> #rgdal: version: 0.8-10, (SVN revision 478)
> #Geospatial Data Abstraction Library extensions to R successfully loaded
> #Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
> #Path to GDAL shared files: /usr/share/gdal
> #Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
> #Path to PROJ.4 shared files: (autodetected)
> #class       : RasterLayer
> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
> #resolution  : 40, 40  (x, y)
> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
> #coord. ref. : NA
> #data source : /home/zambrhe/temp.asc
> #names       : temp
>
> # reading the .asc file form disk
> r2 <- raster("temp.asc")
>
> r2
> #class       : RasterLayer
> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
> #resolution  : 40, 40  (x, y)
> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
> #coord. ref. : NA
> #data source : /home/zambrhe/temp.asc
> #names       : temp
>
> # Summary of the ascii file in the hard disk
> #summary(r2)
> #             temp
> #Min.     128.4340
> #1st Qu.  293.2325
> #Median   371.4120
> #3rd Qu.  499.8195
> #Max.    1805.7800
> #NA's    6097.0000
>
>
> # Summary of the original raster file
> summary(r)
> #             test
> #Min.     128.4340
> #1st Qu.  293.2325
> #Median   371.4120
> #3rd Qu.  499.8195
> #Max.    1805.7800
> #NA's    6097.0000
> ----------  END --------------
>
> Kind regards,
>
> Mauricio Zambrano-Bigiarini, Ph.D
>
> --
> ==============================**===================
> Water Resources Unit
> Institute for Environment and Sustainability (IES)
> Joint Research Centre (JRC), European Commission
> TP 261, Via Enrico Fermi 2749, 21027 Ispra (VA), IT
> Work Phone : +39 0332 789588
> webinfo    : http://floods.jrc.ec.europa.**eu/<http://floods.jrc.ec.europa.eu/>
> ==============================**===================
> DISCLAIMER:
> "The views expressed are purely those of the writer
> and may not in any circumstances be regarded as sta-
> ting an official position of the European Commission"
> ==============================**===================
> "The journey is the reward"  (Steve Jobs)
>
>
> _______________________________________________
> R-sig-Geo mailing list
> [hidden email]
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>

--
*Manuel Spínola, Ph.D.*
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
[hidden email]
[hidden email]
Teléfono: (506) 2277-3598
Fax: (506) 2237-7036
Personal website: Lobito de río <https://sites.google.com/site/lobitoderio/>
Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>

        [[alternative HTML version deleted]]


_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

Re: writeRaster to ascii file (.asc)

Robert Hijmans
I do not the data have changed. Wrong min and max values are reported,
however. Seems to me that in this case gdal does not know the min and
max values but reports  -2147483648, 2147483647. You could check with
GDALinfo. I guess that this is a  shortcoming of (your version of?)
gdal, not rgdal. Robert

On Thu, Jul 11, 2013 at 5:04 AM, Manuel Spínola <[hidden email]> wrote:

> Thank you Mauricio.
>
>> sessionInfo()
> R version 3.0.1 (2013-05-16)
> Platform: x86_64-apple-darwin10.8.0 (64-bit)
>
> locale:
> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
>
> attached base packages:
> [1] grid      stats     graphics  grDevices utils     datasets  methods
> base
>
> other attached packages:
> [1] SDMTools_1.1-13 dismo_0.8-11    raster_2.1-37   rgdal_0.8-10
>  maptools_0.8-25
> [6] lattice_0.20-15 sp_1.0-11       foreign_0.8-54
>
> loaded via a namespace (and not attached):
> [1] tools_3.0.1
>
>
> 2013/7/11 Mauricio Zambrano-Bigiarini <[hidden email]>
>
>> On 10/07/13 18:38, Manuel Spínola wrote:
>>
>>> Dear list members,
>>>
>>> I am trying to export a raster file to a ascii file (.asc) but the min and
>>> max values are not the same. Any idea why?
>>>
>>
>> Hi Manuel,
>>
>> Could you provide your sessionInfo() ?
>>
>>>
>>>  bio1res
>>>>
>>> class       : RasterLayer
>>> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
>>> resolution  : 0.008333333, 0.008333333  (x, y)
>>> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin,
>>> ymax)
>>> coord. ref. : NA
>>> data source : in memory
>>> names       : bio1_23
>>> values      : 47, 272  (min, max)
>>>
>>>  writeRaster(bio1res, filename="temp.asc", format="ascii",overwrite=TRUE)
>>>>
>>> class       : RasterLayer
>>> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
>>> resolution  : 0.008333333, 0.008333333  (x, y)
>>> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin,
>>> ymax)
>>> coord. ref. : NA
>>> data source : /Users/manuelspinola/Dropbox/**DistribucioÌ n de
>>>
>>> especies/World_Clim_Costa_**Rica/temp.asc
>>> names       : temp
>>> values      : -2147483648, 2147483647  (min, max)
>>>
>>
>> Unfortunately, your example is not reproducible, so I tried with a
>> different example and the values in the original and the written ascii file
>> were the same (rgdal_0.8-10, raster_2.1-48, sp_1.0-11):
>>
>> ----------  START --------------
>> library(raster)
>> Loading required package: sp
>>
>> # Creation of a raster file for testing
>> r <- raster(system.file("external/**test.grd", package="raster"))
>>
>> r
>> #class       : RasterLayer
>> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
>> #resolution  : 40, 40  (x, y)
>> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
>> #coord. ref. : +init=epsg:28992 +towgs84=565.237,50.0087,465.**
>> 658,-0.406857,0.350733,-1.**87035,4.0812 +proj=sterea
>> +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000
>> +y_0=463000 +ellps=bessel +units=m +no_defs
>> #data source : /usr/lib64/R/library/raster/**external/test.grd
>> #names       : test
>> #values      : 128.434, 1805.78  (min, max)
>>
>> writeRaster(r, filename="temp.asc", format="ascii",overwrite=TRUE)
>> #rgdal: version: 0.8-10, (SVN revision 478)
>> #Geospatial Data Abstraction Library extensions to R successfully loaded
>> #Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
>> #Path to GDAL shared files: /usr/share/gdal
>> #Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
>> #Path to PROJ.4 shared files: (autodetected)
>> #class       : RasterLayer
>> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
>> #resolution  : 40, 40  (x, y)
>> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
>> #coord. ref. : NA
>> #data source : /home/zambrhe/temp.asc
>> #names       : temp
>>
>> # reading the .asc file form disk
>> r2 <- raster("temp.asc")
>>
>> r2
>> #class       : RasterLayer
>> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
>> #resolution  : 40, 40  (x, y)
>> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
>> #coord. ref. : NA
>> #data source : /home/zambrhe/temp.asc
>> #names       : temp
>>
>> # Summary of the ascii file in the hard disk
>> #summary(r2)
>> #             temp
>> #Min.     128.4340
>> #1st Qu.  293.2325
>> #Median   371.4120
>> #3rd Qu.  499.8195
>> #Max.    1805.7800
>> #NA's    6097.0000
>>
>>
>> # Summary of the original raster file
>> summary(r)
>> #             test
>> #Min.     128.4340
>> #1st Qu.  293.2325
>> #Median   371.4120
>> #3rd Qu.  499.8195
>> #Max.    1805.7800
>> #NA's    6097.0000
>> ----------  END --------------
>>
>> Kind regards,
>>
>> Mauricio Zambrano-Bigiarini, Ph.D
>>
>> --
>> ==============================**===================
>> Water Resources Unit
>> Institute for Environment and Sustainability (IES)
>> Joint Research Centre (JRC), European Commission
>> TP 261, Via Enrico Fermi 2749, 21027 Ispra (VA), IT
>> Work Phone : +39 0332 789588
>> webinfo    : http://floods.jrc.ec.europa.**eu/<http://floods.jrc.ec.europa.eu/>
>> ==============================**===================
>> DISCLAIMER:
>> "The views expressed are purely those of the writer
>> and may not in any circumstances be regarded as sta-
>> ting an official position of the European Commission"
>> ==============================**===================
>> "The journey is the reward"  (Steve Jobs)
>>
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> [hidden email]
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>>
>
>
> --
> *Manuel Spínola, Ph.D.*
> Instituto Internacional en Conservación y Manejo de Vida Silvestre
> Universidad Nacional
> Apartado 1350-3000
> Heredia
> COSTA RICA
> [hidden email]
> [hidden email]
> Teléfono: (506) 2277-3598
> Fax: (506) 2237-7036
> Personal website: Lobito de río <https://sites.google.com/site/lobitoderio/>
> Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>
>
>         [[alternative HTML version deleted]]
>
>
> _______________________________________________
> R-sig-Geo mailing list
> [hidden email]
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>

_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

Re: writeRaster to ascii file (.asc)

Mauricio Zambrano-Bigiarini
On 07/13/2013 01:44 AM, Robert J. Hijmans wrote:
> I do not the data have changed. Wrong min and max values are reported,
> however. Seems to me that in this case gdal does not know the min and
> max values but reports  -2147483648, 2147483647. You could check with
> GDALinfo. I guess that this is a  shortcoming of (your version of?)
> gdal, not rgdal. Robert

Just in case still be useful.

I tried to reproduce the wrongly reported values for min and max using
an .asc file I'm using, but I couldn't:

------------ START -----------
library(raster)

# entering to the directory where the .asc file is stored
# setwd("~")

# Creation of a raster file for testing
r <- raster("temp2.asc")

# rgdal: version: 0.8-10, (SVN revision 478)
# Geospatial Data Abstraction Library extensions to R successfully loaded
# Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
# Path to GDAL shared files: /usr/share/gdal
# Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
# Path to PROJ.4 shared files: (autodetected)

# Summary of the original raster file
summary(r)

#                temp2
# Min.    2.668977e-02
# 1st Qu. 1.474479e-01
# Median  2.109602e-01
# 3rd Qu. 4.575650e-01
# Max.    7.108980e+00
# NA's    1.690500e+04

# Writing the raster from memory to disk (as a new file)
writeRaster(r, filename="temp2_new.asc", format="ascii",overwrite=TRUE)

# class       : RasterLayer
# dimensions  : 254, 262, 66548  (nrow, ncol, ncell)
# resolution  : 5000, 5000  (x, y)
# extent      : -1045000, 265000, -620000, 650000  (xmin, xmax, ymin, ymax)
# coord. ref. : NA
# data source : /home/hzambran/temp2_new.asc
# names       : temp2_new

# reading the .asc file form disk
r2 <- raster("temp2_new.asc")

# Summary of the .asc file in the hard disk
summary(r2)

#            temp2_new
# Min.    2.668977e-02
# 1st Qu. 1.474479e-01
# Median  2.109602e-01
# 3rd Qu. 4.575650e-01
# Max.    7.108980e+00
# NA's    1.690500e+04

# the same summary values are the same as for the original raster (r)

------------ END -----------

I'm using the same rgdal version of Manuel (0.8-10), so as Roger
mentioned before, the error should come from your GDAL version. I'm
using GDAL 1.9.2, released 2012/10/08 (not the latest...)

Kind regards,

Mauricio

--
=================================================
Water Resources Unit
Institute for Environment and Sustainability (IES)
Joint Research Centre (JRC), European Commission
webinfo    : http://floods.jrc.ec.europa.eu/
=================================================
DISCLAIMER:
"The views expressed are purely those of the writer
and may not in any circumstances be regarded as sta-
ting an official position of the European Commission"
=================================================
"The journey is the reward"  (Steve Jobs)


>
> On Thu, Jul 11, 2013 at 5:04 AM, Manuel Spínola <[hidden email]> wrote:
>> Thank you Mauricio.
>>
>>> sessionInfo()
>> R version 3.0.1 (2013-05-16)
>> Platform: x86_64-apple-darwin10.8.0 (64-bit)
>>
>> locale:
>> [1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
>>
>> attached base packages:
>> [1] grid      stats     graphics  grDevices utils     datasets  methods
>> base
>>
>> other attached packages:
>> [1] SDMTools_1.1-13 dismo_0.8-11    raster_2.1-37   rgdal_0.8-10
>>   maptools_0.8-25
>> [6] lattice_0.20-15 sp_1.0-11       foreign_0.8-54
>>
>> loaded via a namespace (and not attached):
>> [1] tools_3.0.1
>>
>>
>> 2013/7/11 Mauricio Zambrano-Bigiarini <[hidden email]>
>>
>>> On 10/07/13 18:38, Manuel Spínola wrote:
>>>
>>>> Dear list members,
>>>>
>>>> I am trying to export a raster file to a ascii file (.asc) but the min and
>>>> max values are not the same. Any idea why?
>>>>
>>>
>>> Hi Manuel,
>>>
>>> Could you provide your sessionInfo() ?
>>>
>>>>
>>>>   bio1res
>>>>>
>>>> class       : RasterLayer
>>>> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
>>>> resolution  : 0.008333333, 0.008333333  (x, y)
>>>> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin,
>>>> ymax)
>>>> coord. ref. : NA
>>>> data source : in memory
>>>> names       : bio1_23
>>>> values      : 47, 272  (min, max)
>>>>
>>>>   writeRaster(bio1res, filename="temp.asc", format="ascii",overwrite=TRUE)
>>>>>
>>>> class       : RasterLayer
>>>> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
>>>> resolution  : 0.008333333, 0.008333333  (x, y)
>>>> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin,
>>>> ymax)
>>>> coord. ref. : NA
>>>> data source : /Users/manuelspinola/Dropbox/**DistribucioÌ n de
>>>>
>>>> especies/World_Clim_Costa_**Rica/temp.asc
>>>> names       : temp
>>>> values      : -2147483648, 2147483647  (min, max)
>>>>
>>>
>>> Unfortunately, your example is not reproducible, so I tried with a
>>> different example and the values in the original and the written ascii file
>>> were the same (rgdal_0.8-10, raster_2.1-48, sp_1.0-11):
>>>
>>> ----------  START --------------
>>> library(raster)
>>> Loading required package: sp
>>>
>>> # Creation of a raster file for testing
>>> r <- raster(system.file("external/**test.grd", package="raster"))
>>>
>>> r
>>> #class       : RasterLayer
>>> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
>>> #resolution  : 40, 40  (x, y)
>>> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
>>> #coord. ref. : +init=epsg:28992 +towgs84=565.237,50.0087,465.**
>>> 658,-0.406857,0.350733,-1.**87035,4.0812 +proj=sterea
>>> +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079 +x_0=155000
>>> +y_0=463000 +ellps=bessel +units=m +no_defs
>>> #data source : /usr/lib64/R/library/raster/**external/test.grd
>>> #names       : test
>>> #values      : 128.434, 1805.78  (min, max)
>>>
>>> writeRaster(r, filename="temp.asc", format="ascii",overwrite=TRUE)
>>> #rgdal: version: 0.8-10, (SVN revision 478)
>>> #Geospatial Data Abstraction Library extensions to R successfully loaded
>>> #Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
>>> #Path to GDAL shared files: /usr/share/gdal
>>> #Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
>>> #Path to PROJ.4 shared files: (autodetected)
>>> #class       : RasterLayer
>>> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
>>> #resolution  : 40, 40  (x, y)
>>> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
>>> #coord. ref. : NA
>>> #data source : /home/zambrhe/temp.asc
>>> #names       : temp
>>>
>>> # reading the .asc file form disk
>>> r2 <- raster("temp.asc")
>>>
>>> r2
>>> #class       : RasterLayer
>>> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
>>> #resolution  : 40, 40  (x, y)
>>> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
>>> #coord. ref. : NA
>>> #data source : /home/zambrhe/temp.asc
>>> #names       : temp
>>>
>>> # Summary of the ascii file in the hard disk
>>> #summary(r2)
>>> #             temp
>>> #Min.     128.4340
>>> #1st Qu.  293.2325
>>> #Median   371.4120
>>> #3rd Qu.  499.8195
>>> #Max.    1805.7800
>>> #NA's    6097.0000
>>>
>>>
>>> # Summary of the original raster file
>>> summary(r)
>>> #             test
>>> #Min.     128.4340
>>> #1st Qu.  293.2325
>>> #Median   371.4120
>>> #3rd Qu.  499.8195
>>> #Max.    1805.7800
>>> #NA's    6097.0000
>>> ----------  END --------------
>>>
>>> Kind regards,
>>>
>>> Mauricio Zambrano-Bigiarini, Ph.D
>>>
>>> --
>>> ==============================**===================
>>> Water Resources Unit
>>> Institute for Environment and Sustainability (IES)
>>> Joint Research Centre (JRC), European Commission
>>> TP 261, Via Enrico Fermi 2749, 21027 Ispra (VA), IT
>>> Work Phone : +39 0332 789588
>>> webinfo    : http://floods.jrc.ec.europa.**eu/<http://floods.jrc.ec.europa.eu/>
>>> ==============================**===================
>>> DISCLAIMER:
>>> "The views expressed are purely those of the writer
>>> and may not in any circumstances be regarded as sta-
>>> ting an official position of the European Commission"
>>> ==============================**===================
>>> "The journey is the reward"  (Steve Jobs)
>>>
>>>
>>> _______________________________________________
>>> R-sig-Geo mailing list
>>> [hidden email]
>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>>
>>>
>>
>>
>> --
>> *Manuel Spínola, Ph.D.*
>> Instituto Internacional en Conservación y Manejo de Vida Silvestre
>> Universidad Nacional
>> Apartado 1350-3000
>> Heredia
>> COSTA RICA
>> [hidden email]
>> [hidden email]
>> Teléfono: (506) 2277-3598
>> Fax: (506) 2237-7036
>> Personal website: Lobito de río <https://sites.google.com/site/lobitoderio/>
>> Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>
>>
>>          [[alternative HTML version deleted]]
>>
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> [hidden email]
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>

_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo
Reply | Threaded
Open this post in threaded view
|

Re: writeRaster to ascii file (.asc)

Manuel Spínola
Thank you Mauricio,

I solved the problem working with tif files instead of grd files.  I don't
know why but that solved the problem, now the ascii file keep the values of
the tif file.

Manuel


2013/7/18 Mauricio Zambrano-Bigiarini <[hidden email]>

> On 07/13/2013 01:44 AM, Robert J. Hijmans wrote:
>
>> I do not the data have changed. Wrong min and max values are reported,
>> however. Seems to me that in this case gdal does not know the min and
>> max values but reports  -2147483648, 2147483647. You could check with
>> GDALinfo. I guess that this is a  shortcoming of (your version of?)
>> gdal, not rgdal. Robert
>>
>
> Just in case still be useful.
>
> I tried to reproduce the wrongly reported values for min and max using an
> .asc file I'm using, but I couldn't:
>
> ------------ START -----------
> library(raster)
>
> # entering to the directory where the .asc file is stored
> # setwd("~")
>
>
> # Creation of a raster file for testing
> r <- raster("temp2.asc")
>
> # rgdal: version: 0.8-10, (SVN revision 478)
> # Geospatial Data Abstraction Library extensions to R successfully loaded
> # Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
> # Path to GDAL shared files: /usr/share/gdal
> # Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
> # Path to PROJ.4 shared files: (autodetected)
>
>
> # Summary of the original raster file
> summary(r)
>
> #                temp2
> # Min.    2.668977e-02
> # 1st Qu. 1.474479e-01
> # Median  2.109602e-01
> # 3rd Qu. 4.575650e-01
> # Max.    7.108980e+00
> # NA's    1.690500e+04
>
> # Writing the raster from memory to disk (as a new file)
> writeRaster(r, filename="temp2_new.asc", format="ascii",overwrite=TRUE)
>
> # class       : RasterLayer
> # dimensions  : 254, 262, 66548  (nrow, ncol, ncell)
> # resolution  : 5000, 5000  (x, y)
> # extent      : -1045000, 265000, -620000, 650000  (xmin, xmax, ymin, ymax)
> # coord. ref. : NA
> # data source : /home/hzambran/temp2_new.asc
> # names       : temp2_new
>
>
> # reading the .asc file form disk
> r2 <- raster("temp2_new.asc")
>
> # Summary of the .asc file in the hard disk
> summary(r2)
>
> #            temp2_new
> # Min.    2.668977e-02
> # 1st Qu. 1.474479e-01
> # Median  2.109602e-01
> # 3rd Qu. 4.575650e-01
> # Max.    7.108980e+00
> # NA's    1.690500e+04
>
> # the same summary values are the same as for the original raster (r)
>
> ------------ END -----------
>
> I'm using the same rgdal version of Manuel (0.8-10), so as Roger mentioned
> before, the error should come from your GDAL version. I'm using GDAL 1.9.2,
> released 2012/10/08 (not the latest...)
>
> Kind regards,
>
> Mauricio
>
> --
> ==============================**===================
>
> Water Resources Unit
> Institute for Environment and Sustainability (IES)
> Joint Research Centre (JRC), European Commission
> webinfo    : http://floods.jrc.ec.europa.**eu/<http://floods.jrc.ec.europa.eu/>
> ==============================**===================
>
> DISCLAIMER:
> "The views expressed are purely those of the writer
> and may not in any circumstances be regarded as sta-
> ting an official position of the European Commission"
> ==============================**===================
>
> "The journey is the reward"  (Steve Jobs)
>
>
>
>> On Thu, Jul 11, 2013 at 5:04 AM, Manuel Spínola <[hidden email]>
>> wrote:
>>
>>> Thank you Mauricio.
>>>
>>>  sessionInfo()
>>>>
>>> R version 3.0.1 (2013-05-16)
>>> Platform: x86_64-apple-darwin10.8.0 (64-bit)
>>>
>>> locale:
>>> [1] en_US.UTF-8/en_US.UTF-8/en_US.**UTF-8/C/en_US.UTF-8/en_US.UTF-**8
>>>
>>> attached base packages:
>>> [1] grid      stats     graphics  grDevices utils     datasets  methods
>>> base
>>>
>>> other attached packages:
>>> [1] SDMTools_1.1-13 dismo_0.8-11    raster_2.1-37   rgdal_0.8-10
>>>   maptools_0.8-25
>>> [6] lattice_0.20-15 sp_1.0-11       foreign_0.8-54
>>>
>>> loaded via a namespace (and not attached):
>>> [1] tools_3.0.1
>>>
>>>
>>> 2013/7/11 Mauricio Zambrano-Bigiarini <[hidden email].**
>>> europa.eu <[hidden email]>>
>>>
>>>  On 10/07/13 18:38, Manuel Spínola wrote:
>>>>
>>>>  Dear list members,
>>>>>
>>>>> I am trying to export a raster file to a ascii file (.asc) but the min
>>>>> and
>>>>> max values are not the same. Any idea why?
>>>>>
>>>>>
>>>> Hi Manuel,
>>>>
>>>> Could you provide your sessionInfo() ?
>>>>
>>>>
>>>>>   bio1res
>>>>>
>>>>>>
>>>>>>  class       : RasterLayer
>>>>> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
>>>>> resolution  : 0.008333333, 0.008333333  (x, y)
>>>>> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin,
>>>>> ymax)
>>>>> coord. ref. : NA
>>>>> data source : in memory
>>>>> names       : bio1_23
>>>>> values      : 47, 272  (min, max)
>>>>>
>>>>>   writeRaster(bio1res, filename="temp.asc",
>>>>> format="ascii",overwrite=TRUE)
>>>>>
>>>>>>
>>>>>>  class       : RasterLayer
>>>>> dimensions  : 382, 407, 155474  (nrow, ncol, ncell)
>>>>> resolution  : 0.008333333, 0.008333333  (x, y)
>>>>> extent      : -85.95, -82.55833, 8.041667, 11.225  (xmin, xmax, ymin,
>>>>> ymax)
>>>>> coord. ref. : NA
>>>>> data source : /Users/manuelspinola/Dropbox/****DistribucioÌ n de
>>>>>
>>>>> especies/World_Clim_Costa_****Rica/temp.asc
>>>>> names       : temp
>>>>> values      : -2147483648, 2147483647  (min, max)
>>>>>
>>>>>
>>>> Unfortunately, your example is not reproducible, so I tried with a
>>>> different example and the values in the original and the written ascii
>>>> file
>>>> were the same (rgdal_0.8-10, raster_2.1-48, sp_1.0-11):
>>>>
>>>> ----------  START --------------
>>>> library(raster)
>>>> Loading required package: sp
>>>>
>>>> # Creation of a raster file for testing
>>>> r <- raster(system.file("external/****test.grd", package="raster"))
>>>>
>>>> r
>>>> #class       : RasterLayer
>>>> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
>>>> #resolution  : 40, 40  (x, y)
>>>> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
>>>> #coord. ref. : +init=epsg:28992 +towgs84=565.237,50.0087,465.****
>>>> 658,-0.406857,0.350733,-1.****87035,4.0812 +proj=sterea
>>>> +lat_0=52.15616055555555 +lon_0=5.38763888888889 +k=0.9999079
>>>> +x_0=155000
>>>> +y_0=463000 +ellps=bessel +units=m +no_defs
>>>> #data source : /usr/lib64/R/library/raster/****external/test.grd
>>>> #names       : test
>>>> #values      : 128.434, 1805.78  (min, max)
>>>>
>>>> writeRaster(r, filename="temp.asc", format="ascii",overwrite=TRUE)
>>>> #rgdal: version: 0.8-10, (SVN revision 478)
>>>> #Geospatial Data Abstraction Library extensions to R successfully loaded
>>>> #Loaded GDAL runtime: GDAL 1.9.2, released 2012/10/08
>>>> #Path to GDAL shared files: /usr/share/gdal
>>>> #Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
>>>> #Path to PROJ.4 shared files: (autodetected)
>>>> #class       : RasterLayer
>>>> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
>>>> #resolution  : 40, 40  (x, y)
>>>> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
>>>> #coord. ref. : NA
>>>> #data source : /home/zambrhe/temp.asc
>>>> #names       : temp
>>>>
>>>> # reading the .asc file form disk
>>>> r2 <- raster("temp.asc")
>>>>
>>>> r2
>>>> #class       : RasterLayer
>>>> #dimensions  : 115, 80, 9200  (nrow, ncol, ncell)
>>>> #resolution  : 40, 40  (x, y)
>>>> #extent      : 178400, 181600, 329400, 334000  (xmin, xmax, ymin, ymax)
>>>> #coord. ref. : NA
>>>> #data source : /home/zambrhe/temp.asc
>>>> #names       : temp
>>>>
>>>> # Summary of the ascii file in the hard disk
>>>> #summary(r2)
>>>> #             temp
>>>> #Min.     128.4340
>>>> #1st Qu.  293.2325
>>>> #Median   371.4120
>>>> #3rd Qu.  499.8195
>>>> #Max.    1805.7800
>>>> #NA's    6097.0000
>>>>
>>>>
>>>> # Summary of the original raster file
>>>> summary(r)
>>>> #             test
>>>> #Min.     128.4340
>>>> #1st Qu.  293.2325
>>>> #Median   371.4120
>>>> #3rd Qu.  499.8195
>>>> #Max.    1805.7800
>>>> #NA's    6097.0000
>>>> ----------  END --------------
>>>>
>>>> Kind regards,
>>>>
>>>> Mauricio Zambrano-Bigiarini, Ph.D
>>>>
>>>> --
>>>> ==============================****===================
>>>> Water Resources Unit
>>>> Institute for Environment and Sustainability (IES)
>>>> Joint Research Centre (JRC), European Commission
>>>> TP 261, Via Enrico Fermi 2749, 21027 Ispra (VA), IT
>>>> Work Phone : +39 0332 789588
>>>> webinfo    : http://floods.jrc.ec.europa.****eu/<http://floods.jrc.ec.*
>>>> *europa.eu/ <http://floods.jrc.ec.europa.eu/>>
>>>> ==============================****===================
>>>> DISCLAIMER:
>>>> "The views expressed are purely those of the writer
>>>> and may not in any circumstances be regarded as sta-
>>>> ting an official position of the European Commission"
>>>> ==============================****===================
>>>> "The journey is the reward"  (Steve Jobs)
>>>>
>>>>
>>>> ______________________________**_________________
>>>> R-sig-Geo mailing list
>>>> [hidden email]
>>>> https://stat.ethz.ch/mailman/**listinfo/r-sig-geo<https://stat.ethz.ch/mailman/listinfo/r-sig-geo>
>>>>
>>>>
>>>>
>>>
>>> --
>>> *Manuel Spínola, Ph.D.*
>>> Instituto Internacional en Conservación y Manejo de Vida Silvestre
>>> Universidad Nacional
>>> Apartado 1350-3000
>>> Heredia
>>> COSTA RICA
>>> [hidden email]
>>> [hidden email]
>>> Teléfono: (506) 2277-3598
>>> Fax: (506) 2237-7036
>>> Personal website: Lobito de río <https://sites.google.com/**
>>> site/lobitoderio/ <https://sites.google.com/site/lobitoderio/>>
>>> Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/**>
>>>
>>>          [[alternative HTML version deleted]]
>>>
>>>
>>> ______________________________**_________________
>>> R-sig-Geo mailing list
>>> [hidden email]
>>> https://stat.ethz.ch/mailman/**listinfo/r-sig-geo<https://stat.ethz.ch/mailman/listinfo/r-sig-geo>
>>>
>>>
>>
>
>
>

--
*Manuel Spínola, Ph.D.*
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
[hidden email]
[hidden email]
Teléfono: (506) 2277-3598
Fax: (506) 2237-7036
Personal website: Lobito de río <https://sites.google.com/site/lobitoderio/>
Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/>

        [[alternative HTML version deleted]]


_______________________________________________
R-sig-Geo mailing list
[hidden email]
https://stat.ethz.ch/mailman/listinfo/r-sig-geo