search for: readgdal

Dear all, I have asked yesterday of how I can read a simple bitmap file in R cran. I was suggest to use either readGDAL or raster for loading my bitmap a. I have done it with readGDAL like     store<-readGDAL(fname='lena256.bmp')     and it works,... but it converts my matrix-like notion of a bitmap to a large vector b. Raster also returns a class that I can not understand. So 1, I want to load it an...

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] "40.702147,-74.015794...

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======== > library(rgdal) > > library(RgoogleMaps) > > GetMa...

I imported the attached tiff file and converted the coordinate system to long lat and graphed it: californiatiff<- readGDAL("california1.tif") proj4string(californiatiff) > rasterprojection <- spTransform(californiatiff, CRS("+proj=longlat") however, when using the plot command for rasterprojection, I get a blob. I can see the outline of the state of california and nevada, but rather than bei...

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'

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 should behave in the same way. From my initial research, the spplot function is the correct func...

...ows that rgdal library is loaded, also the rscproxy library is loaded, in addition to all others. I tried following: imshow(imread('testimg_p.tif')) displays the image properly, however calling this from R using rgdal from within Matlab gives error; >> [a b c] = evalR('img <- readGDAL("testimg_p.tif")') a = ?????[] b = ?????0 c = Invoke Error, Dispatch Exception: Object is static; operation not allowed The command works from R prompt: > img <- readGDAL("testimg_p.tif") testimg_p.tif has GDAL driver GTiff and has 280 rows and 272 columns It i...

...eally I would like to read in the png image and display in its original form as a plot background, hopefully with all its colors represented rather than being converted to a 256 color scale. A code example is below: png(file="file.png", 4000,4000) par(mar=c(0,0,0,0)) myTile <- readGDAL("basemap.png",silent=TRUE); #basemap.png is a PNG file returned from googlemaps myTile at data <- myTile at data[,1:3] col <- SGDF2PCT(myTile,ncolors=256) ## myTile is a spatialGridDataFrame with 3 bands myTile$ind <- col$idx ## add the colour index to the data frame myTi...

...pacetime but has not worked any suggestions? thanks setwd("C:/Documents and Settings/bob/Documenti/tuesday") library(rgdal) dts0 <- c("19980401", "19980701", "19981001", "19990101") fn <- "_Gambia__the_Extract.tif" ndv_gambia <- readGDAL(paste("NDV_", dts0[1], fn, sep="")) for (i in dts0[2:4]) ndv_gambia <- cbind(ndv_gambia, readGDAL(paste("NDV_", i, fn, sep=""))) names(ndv_gambia) <- paste("ndv", dts0, sep="_") dts <- as.Date(dts0, "%Y%m%d") #dal testo...

...ould like to reproject into UTM 18. I was wondering if it was possible in R. Furthermore, I looked into 'rgdal' package, but I can't really find out if I'm doing the right thing. So far, here is what I'm doing: library(rgdal) myTiff = 'L5014028_02819840617_B10.TIF' a = readGDAL(myTiff) aa = spTransform(a, CRS("+proj=utm +zone=18 +datum=WGS84")) However, after few seconds, the memory is full and R crashes. I used other proprietary tools to reproject and it takes only like 10 seconds to reproject. Therefore, I suspect I'm not using the 'rgdal' package...

...ast.m.+North.m. proj4string(data.grid) = CRS("+proj=tmerc +lat_0=0.00000000000 +lon_0=-81.00000000000 +k=0.99960000 +x_0=500000.000000 +y_0=0.0000000 +ellps=GRS80 +units=m") tr <- "e:\\JELA_veg\\test_gtiff\\test.tif" writeGDAL(data.grid["class.pca"], tr) mg3 <- readGDAL(tr) proj4string(mg3) [1] " +proj=utm +zone=17 +ellps=GRS80 +units=m +no_defs" I will really appreciate if anybody can point me in the right dirrection. Thanks, Monica _________________________________________________________________ [[alternative HTML version deleted]]

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) png(filename="hist"+tostring(i)+".png") } <li>I try to do some modyfication of this loop, but it stil...

Dear all, I have a so large image (43,000 x 18,000 pixels) and I need clip this image with a smallest one (1000x1000 pixels). I can read the second image using rgdal package. But the first image can´t be read on my system because if memory limitation (I have about 2GB availabe). So I would like hear from you if anyone have some suggestin in this regards. I have also ArcGis 8.6 ans Erdas 8.3

Is it possible to read jpeg files into R? If yes please guide, Thanks.. I tried to search many time but failed to do. Thankis in advance.. with Best Regards, Malik Shahzad Visiting Researcher National Institute of Informatics (NII) Tokyo, Japan Doctoral Student Asian Institute of Technology (AIT) Bangkok, Thailand +66-8-7676-5616 [[alternative HTML version deleted]]

...rary(rgdal) infile='elev' outfile_base='testout' outfile_ext='.bil' outfile=paste(outfile_base,outfile_ext,sep='') outcon <- file(outfile, "wb") infile_info=GDALinfo(infile) nl=infile_info[[1]] ns=infile_info[[2]] for (row in 1:nl) { templine <- readGDAL(infile,region.dim=c(1,ns),offset=c(row-1,0)) writeBin(templine[[1]], outcon,size=4) } close(outcon) # Below doesn't work # writeGDAL(templine,outfile_base,drivername='EHdr',type="Float32") ### The issue is this: I need to be able to effectively copy the geographic/heade...

I have this https://echange-fichiers.inra.fr/get?k=6TdTdezNeZwAqJtpwVm grib file which is a regular lat-lon 0.25x0.25 degrees and contains 4 fields . I tried to open it using : grib <- readGDAL("C:\\Users\\aalyaari\\Desktop\\Vol_025_H14_2010060700.grib") but I got this error: C:\Users\aalyaari\Desktop\Vol_025_H14_2010060700.grib has GDAL driver GRIB and has 721 rows and 1440 columns Error in validityMethod(as(object, superClass)) : Geographical CRS given to non-conformant...

I am not new to R, but I am new to .grib files. I am downloading some climate data and I would like to analyze it in R. R has a nice netcdf package, but I don?t see any package available to deal specifically with grib files. I see a few posts from other people using grib files in R. However, I was unclear if they used grib files in a different software program and then imported the data

Hello, The other day Justin Peter presented a mini program to plot a topographic map with an overlay of the worldHires. I seemed interesting so I checked the ETOPO5 site and find that there is a new file ETOPO1 with a 1 minute grid. I downloaded it and tried a similar procedure. Now the ETOPO1.gz is 1 Gb and the uncompressed file is 5 Gb. They do not fit into my laptop. I tried the following

Is there any R software that create an image from Yahoo maps together with points of known UTM coordinates (or lat/long marked? Note that my region of interest is not covered in sufficient detail by Google maps. It actually does not have to be Yahoo maps as long as it has sufficient coverage of my region but that's the one I have found with coverage of my region. The scale I am interested in

Um texto embutido e sem conjunto de caracteres especificado associado... Nome: n?o dispon?vel Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070312/112cd7f9/attachment.pl