similar to: geotiff calculations

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.

Dear R users, I'm currently trying to re-project a geotiff in another coordinate system. For instance, I have a tif image in UTM 19 zone which I would 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)

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

Dear all, sorry for this very basic question, but I did not find any good example yet. I would like to set up a variable that can be recall later to substitute a part of a file name. As example: var_filename = as.name("aaa") jpeg("var_filename.jpg") plot() dev.off() It would be very useful in order to avoid to change manually many strings. I guess it is possible, but I did

Estimados Escribo para consultar sobre el uso de modelos mixtos anidados. Los datos que estoy analizando provienen de censos de malezas en cuatro tipos de paisajes de la región pampeana, en los que seleccioné al azar igual número de lotes agrícolas cultivados con tres cultivos (maíz, soja y trigo-soja). En cada lote censé el número de especies de malezas en tres posiciones: el alambrado, el borde

Dear all, I have been trying to speed up a process we have been done in ArcGIS. We have to read a single layer TIFF (monochrome image) in . For this, I have used the "rtiff" package. After reading the TIFF file, I compared the raw values for each pixel that I have in ArcGIS to the ones obtained in R. In ArcGIS I have discrete values in the range 0..255, while in R I have continuous

Dear members, How can I read and plot an ArcGIS raster file into R ? The file has extension .aux and contains floating point bathymetry data. The purpose is to create a spatial model in R that uses ArcGIS map data. I have managed to read and plot various shape files into my R project, but I am stuck with this now. I am new to this list and also to R, so any help would be much appreciated. Many

Hello every body, i?m testing some w32 apps with wine and I found some problems. Trying to install ArcGIS (My first installation), I had the following error: Need internet explorer 6. Well I installed IE 6, and continue installing ArcGIS, it was installed succesfully, but it installed Phyton, and throw some problems. Well finally installed I run ArcGIS and the fonts were crazy, I didnt see

Can anyone tell me why I would get different average nearest neighbor values for the same set of coordinates between ArcGIS 10 and R? Sometimes the difference in distance is over 1.3 km. Alexis -- View this message in context: http://r.789695.n4.nabble.com/nndist-R-vs-ArcGIS-tp3442375p3442375.html Sent from the R help mailing list archive at Nabble.com.

On R 2.15.2 and ArcGIS 9.3.1, it works for me in ArcCatalog but you have to follow the particulars here: http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=Accessing_delimited_text_file_data For example: write.table(test, '***.tab', sep = '\t', row.names = F) The extension .tab and sep = '\t' are required for text files. Didn't test row.names=T but I

Hi everybody, I heard about "R(D)-COM Stat connector" for ArcGIS, but i am not sure what that is. I did a search in the archive but it seems i am not getting anything back. can anybody explain me what that is, and where i can find more info about it? There is any possibility to run R from inside ArcGIS? there is more than RArcInfo and Shapefile which can

Dear list, I just tried to use the function cluster.stat in the package fpc. I just have a couple of questions about the syntax: cluster.stats(d,clustering,alt.clustering=NULL, silhouette=TRUE,G2=FALSE,G3=FALSE) 1) the distance object (d) is an object obtained by the function dist() on my own original matrix? 2) clustering is the clusters vector as result of one of the many clustering methods?

It seems within the last ~3 months Ive been having issues with writing text or csv files from a R data frame.  The problem is multifold and it is hard to filter  out what is going on and where the problem is.  So, Im hoping someone else has come across this and may provide insight.   My current settings for R: R version 2.15.2 (2012-10-26) Platform: x86_64-w64-mingw32/x64 (64-bit) locale: [1]

Can someone help me out with a "small" problem? I've started using netcdf files recently, and I want to extract the grid id and also the coordinates from a HIRHAM netcdf file. I know how to extract a slice of dataset both in space and in time and I also know the area that this file should cover, however I have no idea regarding the reference for both LAT/LON and RLAT/RLON. I tried

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

Fwiw the engine behind geom_raster needs explicit observation-per-row form for input (with no structural normalization), so conversion to points is perfectly proper here, albeit confusing in context. (It's closer to what graphics devices actually use ultimately, but the expansion is laid out very early in ggplot2 because there's no standard for intermediate forms.) Cheers, Mike On Wed,

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

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