similar to: readGDAL or raster for reading bit map files

Hello: I would like to superimpose vector polygons (state outlines) from a Shape file on top of a satellite image, imported into a SpatialGridDataFrame from GEOTIFF via gdal_translate and readGDAL. When I plot polygon and point shape files in R, into SpatialPointDataFrame and SpatialPolygonDataFrame, the two feature sets line up geographically, so it seems logical that a SpatialGridDataFrame

Dear R users, I was wondering how I can convert a raster image (that made R through interpolation) into an ascii or csv format? this is the last line of my command p <- interpolate(r, tpsfit) So p is my raster file which I want to convert into ascii or csv Many thanks Regards Mintewab

library(rgdal) my_asc=dir("~/Pulpit/dods/karol/TVDI 113_121",pattern=".asc",recursive=T,full.names=T) for (i in 1:length(my_asc)) { r <- readGDAL(my_asc[i]) z <- as.matrix(r) vectordata[i] <- mean(z) vectordatamax[i] <- max(z) vectordatamin[i] <- min(z) vectordev[i] <- sd(z,na.rm=True) hist(z)

Dear All, I'v tried to read in data in numerous ways including:- Spain10km<-data.frame(readAsciiGrid("F://RMap//sp10kpointid1.aux")) Error in readAsciiGrid("F://RMap//sp10kpointid1.aux") : object 'cellsize' not found In addition: Warning message: In readLines(t, n = 6) : incomplete final line found on 'F://RMap//sp10kpointid1.aux'

Hi all, I am trying to plot a raster object (I can explain why but the point is that it would be a raster objeçt).. I have selected a small code to show you exactly the problem require(raster) test<-matrix(data=runif(10000),nrow=100) m<-raster(test) plot(m,axes="FALSE") axis(1,at=c(0,1),labels=c("a","b")) # THIS DOES NOT CHANGE THE INVISIBLE AXIS WHILE

Dear all, I am looking for a function that can read a simple bmp file and return me a matrix with their corresponding values. I have already looked in the packages require(ReadImages) require(pixmap) and also did apropos("image") but I could not find something that can do something like that. I have also tried the read.ftable but it also failed. B.R Alex [[alternative HTML

Hello list I am trying to use a Googlemaps tile (png file, 640 X 640 px tile) as a background to a plot and have been using the rgdal library to read in the PNG file (modified from code in the RGoogleMaps package). This works OK. My problem is is that the SGDF2PCT function in rgdal seems to be limited to 256 colors and this is introducing (I think) some degradation to the image which appears on

Dear all, Several days ago, I posted How to write a bmp file pixel by pixel. Instead of bmp, I succeeded in writing a PPM file by using the pixmap package. Thanks for the hint generously provided by Uwe Ligges. Now I have a new question. How to convert a PPM file to BMP file in R? I know I can do this in photoshop or by some other softwares, but I think if I Can do this in R, that

Dear listers, I have a collection of tif images that I would like to convert, in R, to a matrix containing the values of the 8bit colour. Ideally, I would like a matrix for each of the colour channels (red, blue and green). I have 'googled' and searched the help list but have yet to find a solution and hope that someone can point me in the right direction. I currently use

Hi, everybody! i created a imagem by kriging using geoR package. I imported points from GRASS("zn", after converted to geodata "zn_geo"), the border "zn_border" and a raster mask. Then i interpolated the points by kriging and created a raster image. Now, i need export this image back to GRASS to use it in the module r.mapcalc. I can't do it. I tried use

Dear All, I have a raster map of the class 'SpatialPointsDataFrame' and coordinates of the class 'SpatialPoints'. I would like to retrieve the values that are contained in the raster map at the specific locations given by the coordinates. Can anyone help me out? Kind regards, Katrin Fleischer

Hello everyone. I am trying to build up understanding in R by trying to develop just some simple scenarios. I would like to explain you what I am trying to do and what I did so far. I would like to put inside a RasterLayer (raster package) a Gaussian field (for given covariance) using grf function (geoR package) 1. First I created a Raster Layer object r <- raster() # Default values are

Hello, I have been exploring the possibility to transition some code that currently uses image() to use the new rasterImage(). To date, I haven't been able to specify a color look-up strategy that works. For example... nx <- 100 ny <- 100 m <- matrix(data = rep(seq(0,1, length = nx), ny), ncol = nx, nrow = ny, byrow = TRUE) plot(1:nx, 1:ny, type = "n") my.color

Hello everyone The following two commands plot.default(seq(1,5),seq(2,6)) plot(seq(1,5),seq(2,6)) plot nothing. One day ago this would create a simple plot diagram but unfortunately right now no plot appears. ?plot returns Help on topic 'plot' was found in the following packages: Plot a Raster* object (in package raster in library

Hi, I just imported two raster maps into R using the SPGRASS6 package, one containing elevation data and the other containing an erosion index: Kar_inc <-readRAST6("Incis_Kar", plugin=FALSE) Kar_dem <- readRAST6("DEM_Kar", plugin=FALSE) I just wanted to make a xy plot of erosion parameter vs elevation. How does this work? I don't get how to handle

Dear Ruser, I am trying to use the map library. But I have been problem the follow message. Warning message: *In readGDAL(destfile, silent = TRUE) : GeoTransform values not available * library(rgdal) library(RgoogleMaps) > GetMap(markers = '40.702147,-74.015794,blues%7C40.711614,-74.012318,greeng%7C40.718217,-73.998284,redc', destfile = "MyTile1.png"); Read 1 item [1]

This code will read binary file and display it as a map. may problem is that this code is using a continuous colour scheme, even though I have discrete data (which is a classification scheme). How can I map numbers to colours with raster? Please require(raster) conne <- file("C:\\lai.bin", "rb") sd<- readBin(conne, integer(), size=1, n=360*720, signed=F)

R-geographers... I'm trying to solve a problem to implement a line-by-line tiled processing using RGDAL (read 1 line of an image, process the one line, write one line of the image to a binary file). Everything except for the final step I'm able to do using a combination of RGDAL and r-base commands. Below is the basic structure, the input file "elev" can be any image file

Hi I want to overlay my raster and its boundary which is a shapefile. When I used thise code separately, all is ok: # Open raster >Image<-read.table("C:\\Users\\Documents\\Nouveau\\Frequence.txt",sep="",dec=",",header=TRUE) >testo<-rasterFromXYZ(Image) >plot(testo) >testo2 <- aggregate(testo,fact=10, fun=mean) >plot(testo2) # open

Dear Ruser, > > > I am trying to use the map library. But unfortunately I am facing the > warning message above: > Warning message: > *In readGDAL(destfile, silent = TRUE) : GeoTransform values not available > * I want to use that library because I intend to plot some station in a Brazilian map. I really appreciate any help. > > ===here what I was testing========