similar to: Re-projecting geotiff

Dear list, I have to compare two digital elevation models in raster format (geotiff). I then have to calculate the differences in altitude for each cell and make some statistics (basic as mean, median, std, range but also more "advanced" as RMSE) on that. I do not know very much how to proceed: 1) is it possible to import the geotiff in R? If so with which package? if not which is the

Hi, Doing more search i've discovered package RGDAL that can write a geotiff file with projection. I saved a geotiff file in UTM projection and if i read the file back in R and check the projection seems that everything is OK. But if i load the file in ArcGIS (ESRI product) i get the warning that the file is missing spatial reference so it cannot be projected, but it is displayed correctly.

Hi I have a netCDF file of volumetric soil water content at four different soil layers and want to convert each soil layer in the netCDF file to a GeoTIFF layer. This code converts the netCDF file to one GeoTIFF layer, i.e. unclear which soil depth. file.nc <- "C:/Soil_Weather_data/Agro/VMC21/VMC21.nc" file.tiff <- "C:/Soil_Weather_data/Agro /VMC21/VMC21" importnetcdf

Hi Ahmed, When reading from a ncdf file you can use the 'varname', 'lvar' and 'level' arguments - see the 'Details' section in the docs https://www.rdocumentation.org/packages/raster/versions/2.6-7/topics/raster <https://www.rdocumentation.org/packages/raster/versions/2.6-7/topics/raster> We can't tell what is in the ncdf file from what you report other

Hi, I have a matrix of data which i can vizualize as an image - for example. I would like to save this image as a geotiff file or at a tiff file with a world file which holds the projection of my data (ultimately the data represent a map of some sort). I know i can save the data as an ESRI grid, but i am not interested in that. I wonder if anybody knows about any code which will help me do

Hi all, We followed some books and sample codes and did some EMM estimation, only to find it won't be able to generate forecast. This is because in the stochastic volatility models we are estimating, the volatilities are latent variables, and we want to forecast 1-step ahead or h-step ahead volatilities. So it is nice to have the system estimated, but we couldn't get it to forecast at

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

Hi everyone, I wonder if anyone knows the reason why the outputs of the same reprojection in r and arcgis are different?. The magnitude of the change can be up to 40 km in the poles. Basically, I have a database of points equally separated by one degree over the globe. In ARCGIS, I am projecting the data in GCS-WGS-1984 and then reprojected it to Berhmann to ensure equal area

GDAL Package for R http://sourceforge.net/projects/rgdal The R GDAL package is an interface for accessing Frank Warmerdam's Geographic Data Abstraction Library from within R. GDAL is capable of reading and writing a wide range of geographic data formats including ESRI grid format and geotiff. On Wed, 2004-11-17 at 09:09, Witold Eryk Wolski wrote: > <SNIP> > SourceForge.net

I know there is some integration between Grass and R available, but has anyone produced a link between R and MFWorks. I am considering Buying MFWorks (Windows) and it would be useful to know if anyone has written such a link. Many thanks, Graham S -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read

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

In the manual page for prcomp(), it says that sdev is "the standard deviations of the principal components (i.e., the square roots of the eigenvalues of the covariance/correlation matrix, though the calculation is actually done with the singular values of the data matrix)." ?However, this is not what I'm finding. ?The values appear to be the standard deviations of a reprojection of

Saludos, estoy intentando usar la función reprojectHDF() ( http://r-gis.net/?q=ModisDownload). #### source('ModisDownload.R') library(raster) library(RCurl) path<-"~/R/MODIS/Data/Test" input<-list.files(path,"*.hdf",all.files=T,recursive=T,full.names=T) input outname<-paste(substr(input[1],40,45),".tif",sep='') outname

Hi I am trying to install the rgdal package through r-cran-rgdal (I had added the ppa), and run into an error in Ubuntu 12.04: r-cran-rgdal : D?pend: libgdal1 (>= 1.9.0) mais il n'est pas installable in English, libgdal1 1.9 does not seem available in Ubuntu 12.04 (I think version 1.7) comes. And indeed, a: $sudo apt-get install libgdal1 Gives also an error message, saying libgdal1

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)

Dear all, I'm working with satellite images in R and plotting them via the code below. I was wondering whether coordinates (spatial["x"], spatial["y"]) are used as centre coordinates of the pixels in the GRID? In this script; spatial["x"] & spatial["y"] are the centre coordinates of the satellite image pixels. I'm asking this because some

> > Hi, > > I used > GDALinfo("MOD13Q1.A2001049.h13v11.005.2007002215512.250m_16_days_EVI.tif") and > got the results: > > rows 10 > columns 11 > bands 1 > origin.x 150701.4 > origin.y 7744897 > res.x 250 > res.y 250 > ysign -1 > oblique.x 0 > oblique.y 0 > driver GTiff >

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

Hi, I have a dataset (CSV) with some counts of firms located around the globe. Each count is assigned to the longitude and latitude of the specific location. Now I want to plot these counts on a world map using dots (size of dots represent the count). I have been unable to find any info on whether this is possible and if so, how? Can you please help me? Thanks! Mathijs -- View this message

Hello all and thanks for the help. I am analyzing a dataset using MetaMDS and I would like to project some extra samples onto the plot such that the extra samples do not play a role in defining the axes. I have been thinking of different ways of doing this and I was wondering if anyone had a suggestion for an easy way to do this (possibly a way to weight samples in the analysis?) Thanks --