[R-sig-Geo] Converting raster to a matrix

[John Wasige]

​Dear all,
Greetings!

I want to convert a big raster (see attached) a matrix, however, I get the
error below from the script. Anybody with an Idea on what is going wrong in
my script?

###Errors at the following point in the script​

> mat <- as.matrix(rast)Error in .local(x, ...) : 'x' has no values

#####script



​library(raster)
library(rgdal)

setwd('D:/Mobile_testdata')
image <- raster('D:/Mobile_testdata/SOSug_mean2qgis.tif')

rast <-  raster(image)
mat <- as.matrix(rast)

Thanks for your help

Rgds John

<PAMDataset>
  <PAMRasterBand band="1">
    <Histograms>
      <HistItem>
        <HistMin>4</HistMin>
        <HistMax>362</HistMax>
        <BucketCount>256</BucketCount>
        <IncludeOutOfRange>1</IncludeOutOfRange>
        <Approximate>0</Approximate>
        <HistCounts>1|0|1|0|1|0|0|0|1|0|0|0|0|0|1|1|0|0|0|0|1|1|1|0|1|0|0|2|0|0|1|2|1|2|0|0|0|0|1|2|0|1|5|0|2|4|1|2|3|4|1|14|6|11|33|27|65|32|39|94|37|90|81|75|117|63|87|37|22|62|21|37|21|34|67|31|62|17|20|50|17|26|12|11|26|13|9|3|4|11|2|4|3|1|3|0|2|1|1|4|0|1|0|0|3|1|2|0|0|0|1|0|0|1|0|0|1|0|0|0|1|0|0|0|0|0|0|0|0|0|1|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|0|1|1|0|0|0|0|0|0|0|0|0|0|0|0|1|0|0|1|0|0|0|1|0|0|0|0|0|0|0|0|0|1|0|0|0|0|0|0|1|0|0|12|34|18|35|11|12|18|5|15|7|6|18|8|14|3|7|4|3|7|4|1|1|2|9|1|2|1|2|4|3|5|1|0|4|2|3|2|0|4|3|3|0|0|2|2|2|0|0|0|0|1|1|0|2|0|2|0|1|0|0|3|1|0|1</HistCounts>
      </HistItem>
    </Histograms>
    <Metadata>
      <MDI key="STATISTICS_MAXIMUM">362</MDI>
      <MDI key="STATISTICS_MEAN">131.04042553191</MDI>
      <MDI key="STATISTICS_MINIMUM">4</MDI>
      <MDI key="STATISTICS_SKIPFACTORX">1</MDI>
      <MDI key="STATISTICS_SKIPFACTORY">1</MDI>
      <MDI key="STATISTICS_STDDEV">74.594661020205</MDI>
    </Metadata>
  </PAMRasterBand>
</PAMDataset>

SOSug_mean2qgis.tfw
Description: Binary data

SOSug_mean2qgis.tif.vat.cpg
Description: Binary data

SOSug_mean2qgis.tif.vat.dbf
Description: Binary data

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo@r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo