Re: [R] Errors in Raster to Point

2016-09-20 Thread David Remotti
First answer is that R is not the proper environment for such a 
question. There a re many free package for image analysis or even GIS. 
Try for example Q-GIS


David


Il 19/09/2016 22:34, GwanSeon Kim ha scritto:

Hi, all
I am just beginner to use R.
I am working with TIF image file, and the information about the raster is
following:

class   : RasterLayer
dimensions  : 11150, 21808, 243159200  (nrow, ncol, ncell)
resolution  : 30, 30  (x, y)
extent  : 569685, 1223925, 1513995, 1848495  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0
+y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0
data source :
C:\Users\Gwan\AppData\Local\Temp\Rtmpg506Ee\raster\r_tmp_2016-09-14_122409_6260_09589.grd
names   : test_map
values  : 1, 225  (min, max)
attributes  :
ID OBJECTID Value Red Green Blue   Count   Class_Name
Opacity
  from:  02 1   1 00 5982503 Corn
 1
  to  : 48  255   254   0 00   10336 Dbl Crop Barley/Soybeans
 1



>From this Rasterlayer, I want to convert raster to point for each pixel
based on "Value (one of column name)" and create a raster with
georeferenced information.
I used code as following: RP <- rasterToPoints(KY_raster)
However, I could not get the points and have an error message "cannot
allocate vector of size 5.4 Gb" and "Your computer is low on memory. Save
your files and close these programs".
Could someone please help me how I can convert to raster to points??
Best,

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.



__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


Re: [R] Errors in Raster to Point

2016-09-20 Thread Michael Sumner
On Tue, 20 Sep 2016, 15:55 GwanSeon Kim  wrote:

> Hi, all
> I am just beginner to use R.
> I am working with TIF image file, and the information about the raster is
> following:
>
> class   : RasterLayer
> dimensions  : 11150, 21808, 243159200  (nrow, ncol, ncell)
> resolution  : 30, 30  (x, y)
> extent  : 569685, 1223925, 1513995, 1848495  (xmin, xmax, ymin, ymax)
> coord. ref. : +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0
> +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0
> data source :
>
> C:\Users\Gwan\AppData\Local\Temp\Rtmpg506Ee\raster\r_tmp_2016-09-14_122409_6260_09589.grd
> names   : test_map
> values  : 1, 225  (min, max)
> attributes  :
>ID OBJECTID Value Red Green Blue   Count   Class_Name
> Opacity
>  from:  02 1   1 00 5982503 Corn
> 1
>  to  : 48  255   254   0 00   10336 Dbl Crop Barley/Soybeans
> 1
>
>
>
> >From this Rasterlayer, I want to convert raster to point for each pixel
> based on "Value (one of column name)" and create a raster with
> georeferenced information.
> I used code as following: RP <- rasterToPoints(KY_raster)
> However, I could not get the points and have an error message "cannot
> allocate vector of size 5.4 Gb" and "Your computer is low on memory. Save
> your files and close these programs".
> Could someone please help me how I can convert to raster to points??
> Best,
>


The First question is why? The point (centre)coordinates of every pixel are
massively redundant since they are a simple function of cell index and the
raster's extent.

You might try as.data.frame with xy =TRUE to avoid any overhead in casting
to Spatial, but still it's very likely that this is just a step towards
your actual goal. Tell us what you want to do and I am sure there is a
better way.

Cheers, Mike

>
> [[alternative HTML version deleted]]
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
-- 
Dr. Michael Sumner
Software and Database Engineer
Australian Antarctic Division
203 Channel Highway
Kingston Tasmania 7050 Australia

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.


[R] Errors in Raster to Point

2016-09-19 Thread GwanSeon Kim
Hi, all
I am just beginner to use R.
I am working with TIF image file, and the information about the raster is
following:

class   : RasterLayer
dimensions  : 11150, 21808, 243159200  (nrow, ncol, ncell)
resolution  : 30, 30  (x, y)
extent  : 569685, 1223925, 1513995, 1848495  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=aea +lat_1=29.5 +lat_2=45.5 +lat_0=23 +lon_0=-96 +x_0=0
+y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80 +towgs84=0,0,0
data source :
C:\Users\Gwan\AppData\Local\Temp\Rtmpg506Ee\raster\r_tmp_2016-09-14_122409_6260_09589.grd
names   : test_map
values  : 1, 225  (min, max)
attributes  :
   ID OBJECTID Value Red Green Blue   Count   Class_Name
Opacity
 from:  02 1   1 00 5982503 Corn
1
 to  : 48  255   254   0 00   10336 Dbl Crop Barley/Soybeans
1



>From this Rasterlayer, I want to convert raster to point for each pixel
based on "Value (one of column name)" and create a raster with
georeferenced information.
I used code as following: RP <- rasterToPoints(KY_raster)
However, I could not get the points and have an error message "cannot
allocate vector of size 5.4 Gb" and "Your computer is low on memory. Save
your files and close these programs".
Could someone please help me how I can convert to raster to points??
Best,

[[alternative HTML version deleted]]

__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.