[R-sig-Geo] writeRaster error: unknown file format
pecardoso
pauloeducardoso at gmail.com
Tue Jan 28 02:01:43 CET 2014
Robert,
Thank you for your message. First thing, raster package is simply amazing.
I've installed raster v2.2-16.
I've imported the TIF from the dropbox link with
i.l8 <- raster('D:\\Dropbox\\Public\\LC81810682013122LGN01_B2.TIF',
package = "raster")
> i.l8
class : RasterLayer
dimensions : 7311, 7521, 54986031 (nrow, ncol, ncell)
resolution : 30, 30 (x, y)
extent : 370485, 596115, -1388715, -1169385 (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=33 +datum=WGS84 +units=m +no_defs
+ellps=WGS84 +towgs84=0,0,0
data source : D:\Dropbox\Public\LC81810682013122LGN01_B2.TIF
names : LC81810682013122LGN01_B2
values : 0, 65535 (min, max)
I'm also able to generate a RST file from the rasterLayer i.l8 with
writeRaster(band, filename = 'D:\\l8_b2.rst', format = 'IDRISI', overwrite
= TRUE)
IDRISI Taiga can read the RST file but the file format is tagged as "IDRISI
x A1" and not "IDRISI Raster A1"
I still not being able to generate the RST file with
writeRaster(band, filename = 'D:\\l8_b2.rst', format = 'RST', overwrite =
TRUE)
Error in .local(.Object, ...) : Unable to create dataset
But it works fine with writeGDAL:
writeGDAL(as(i.l8, "SpatialGridDataFrame"), fname = "D:\l8_gdal.rst",
drivername = "RST")
IDRISI Taiga read it and tag the file format as "IDRISI Raster A1", just
like the ones generated internally with import functions from IDRISI.
writeRaster with format = 'IDRISI' generetes a 230,732 MB idrisi file while
a 214,790 MB file is generated with writeGDAL and drivername = 'RST'
I can't figure out why.
Paulo Eduardo Cardoso
2014-01-27 Robert Hijmans [via R-sig-geo] <
ml-node+s2731867n7585643h37 at n2.nabble.com>
> Paulo,
>
> This has been fixed in raster 2.2-16 (development version). You can
> install and try it:
> install.packages("raster", repos="http://R-Forge.R-project.org")
>
> Robert
>
> On Fri, Jan 24, 2014 at 7:53 AM, pecardoso <[hidden email]<http://user/SendEmail.jtp?type=node&node=7585643&i=0>>
> wrote:
>
> > Hello all. Thanks for your time.
> >
> > Robert, True about format, but 'IDRISI' file format is also there and
> I've
> > used format = 'IDRISI' ultil recently, and it was working just fine with
> > exactly the same procedure.
> >
> > About running scripts from RStudio and R console, I'm not skilled enough
> to
> > understand the differences.
> >
> > in my system, writeFormats() returns:
> >
> > name long_name
> > [1,] "raster" "R-raster"
> > [2,] "SAGA" "SAGA GIS"
> > [3,] "IDRISI" "IDRISI"
> > ...
> > [36,] "RST" "Idrisi Raster A.1"
> >
> > With older versions of raster package, the argument format='IDRISI' also
> > produced 'IDRISI Raster A.1' files.
> >
> > Latest version in fact is not accepting format = 'IDRISI' for the
> > RasterLayer created from reading the 16 bis Landsat TIF file and is
> > returning the error:
> > Error in .setFileExtensionHeader(filename, filetype) : unknown file
> format
> >
> > Changing to format = 'RST' isn't better. I get the error:
> > Error in .local(.Object, ...) : Unable to create dataset
> >
> > The behavior is the same running this from within RStudio and from the
> R
> > console.
> >
> > Maybe with the same data it would be better, or reproducible.
> >
> > Here
> > <https://www.dropbox.com/s/is2td4t6uwc44fq/LC81810682013122LGN01_B2.TIF>
>
> > the link to the landast image (107Mb).
> >
> > For the resulting RasterLayer
> > GDALinfo('...\\LC81810682013122LGN01_B2.tif')
> > rows 7311
> > columns 7521
> > bands 1
> > lower left origin.x 370485
> > lower left origin.y -1388715
> > res.x 30
> > res.y 30
> > ysign -1
> > oblique.x 0
> > oblique.y 0
> > driver GTiff
> > projection +proj=utm +zone=33 +datum=WGS84 +units=m +no_defs
> > file C:\Users\Strix\Downloads\LC81810682013122LGN01_B2.tif
> > apparent band summary:
> > GDType hasNoDataValue NoDataValue blockSize1 blockSize2
> > 1 UInt16 FALSE 0 1 7521
> > apparent band statistics:
> > Bmin Bmax Bmean Bsd
> > 1 0 65535 NA NA
> > Metadata:
> > AREA_OR_POINT=Point
> > Warning message:
> > statistics not supported by this driver
> >
> > Eddie, I generally do not use setwd when sharing scripts with
> colleagues. I
> > usually define paths to some shared folder. It's not working anyway.
> From my
> > sessionInfo string you will se that I'm running rgdal.
> >
> > Could it be something with rgdal (and gdal) installation after package
> > update? For a similar problem
> > <
> http://r-sig-geo.2731867.n2.nabble.com/problem-with-writeRaster-and-GeoTIFF-format-td7583541.html>
>
> > users suggests re-install of rgdal, gdal and its libraries.
> >
> >
> >
> > --
> > View this message in context:
> http://r-sig-geo.2731867.n2.nabble.com/writeRaster-error-unknown-file-format-tp7585627p7585632.html
>
> > Sent from the R-sig-geo mailing list archive at Nabble.com.
> >
> > _______________________________________________
> > R-sig-Geo mailing list
> > [hidden email] <http://user/SendEmail.jtp?type=node&node=7585643&i=1>
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
> _______________________________________________
> R-sig-Geo mailing list
> [hidden email] <http://user/SendEmail.jtp?type=node&node=7585643&i=2>
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
>
> ------------------------------
> If you reply to this email, your message will be added to the discussion
> below:
>
> http://r-sig-geo.2731867.n2.nabble.com/writeRaster-error-unknown-file-format-tp7585627p7585643.html
> To unsubscribe from writeRaster error: unknown file format, click here<http://r-sig-geo.2731867.n2.nabble.com/template/NamlServlet.jtp?macro=unsubscribe_by_code&node=7585627&code=cGF1bG9lZHVjYXJkb3NvQGdtYWlsLmNvbXw3NTg1NjI3fDIxMDEwOTczNDM=>
> .
> NAML<http://r-sig-geo.2731867.n2.nabble.com/template/NamlServlet.jtp?macro=macro_viewer&id=instant_html%21nabble%3Aemail.naml&base=nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.naml.namespaces.BasicNamespace-nabble.view.web.template.NabbleNamespace-nabble.view.web.template.NodeNamespace&breadcrumbs=notify_subscribers%21nabble%3Aemail.naml-instant_emails%21nabble%3Aemail.naml-send_instant_email%21nabble%3Aemail.naml>
>
--
View this message in context: http://r-sig-geo.2731867.n2.nabble.com/writeRaster-error-unknown-file-format-tp7585627p7585646.html
Sent from the R-sig-geo mailing list archive at Nabble.com.
More information about the R-sig-Geo
mailing list