similar to: Memory limit problem

I have a simple problem: I need to load a ERSI shapefile of US states and check whether or not a set of points are within the boundary of these states. I have the shapefile, I have the coordinates but I'm having a great deal of difficulty bringing the two together. The problem is the various GIS packages for R do not play well with each other. sp, shapefiles, maptools, etc all use different

Hi I'm trying to analyse some data and need to set the geographic coordinate system before I can do the analysis. I've been trying to use the project command in rgdal but keep getting an error message saying: Error in project(locationsMatrix, PROJECTION.OUT) : latitude or longitude exceeded limits ( PROJECTION.OUT <- "+proj=aea +lat_1=-18 +lat_2=-36 +lat_0=0 +lon_0=132

Hola a todos Esta es mi primera pregunta en el grupo, y es sencilla pero me tiene atascado. Estoy tratando de cortar mi mapa de (poner limites en UTM) en un lugar definido como mi area de estudio (en este caso el sur de chile). Pero creo no estar usando bien la función CRS ponendo bien los limites requeridos. > study_area <- readRDS("CHL_adm0.rds") > study_area_UTM <-

Hello R Community, I have imported a dataset which contain X Y coordinates and would like to recreate a shape file after some data analysis. What i have done is to import some taxlot data and join them based on some criteria. I want to check to see how well the joining went by reviewing the results in GIS. A couple things. I cant seem to import a shape file correctly using the maptools

I have imported a CSV file: rfishR <- read.csv(file="rfishR.csv",stringsAsFactors = FALSE, strip.white = TRUE, na.strings = c("NA","") ) attach(rfishR) When I call it up in R, it starts with line 2066 rather than 1 and some of the headers (used Headers = TRUE, too) are "masked"? Sample data loc lat lon datum water date obs net

Hi, ? I am trying to write a shapefile with projection. I have my data in a data.frame called try and consists in xy coordinates and a numerical attribute value z1. ? Libraries loaded are: sp, rgdal, raster, maptools ? head(try) ???????? x?????? ??? y??? ? ? ? ? ? ? ??? z1 1 610237.1 ???????? 3375751 ???????????? 8.221 2 610236.1 ???????? 3375750 ???????????? 8.153 3 610236.1 ????????

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 seem to be having a problem adding the following two raster objects together - one is a rasterLayer, the other is a rasterBrick. The extent, resolution, and origin are the same, so according to my understand it should work. The objects look like so: > obs.clim class : RasterLayer dimensions : 60, 200, 12000 (nrow, ncol, ncell) resolution : 0.5, 0.5 (x, y) extent : -70,

I am trying to set up a Grass project and need to set up the region so that I can view the map. I can look at a map and find the lat/lon, but the map projection is in UTM NAD38 WGS84 and I need to set the eastings and northings. Is there a package that will help me calculate this in R. thanks -- Stephen Sefick Let's not spend our time and resources thinking about things that are so little

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

Salutations, Long story short: I need to create an SPDF to superimpose on a Google map. I have .shp file which has a long list of locations listed in (what I guess must be) UTM. I need to transform these into lon, lat. My use of spTransform, however, keeps giving me the wrong results. For what it's worth, it only seems to be messing up the latitude. My final line of code looks like this:

Dear all, I am trying to project my LongLat-maps to a plane. The ultimate purpose is to do a search of points in vicinity of other points using overlay-commands (sp) with radius in km. I am applying spTransform (package rgdal) and it gives my some curious results. An example. Let's take a point lying somewhere in Germany, zone=32U x <- 8.968735 y <- 49.454735 After conversion I

I have population data (?totpopE?) at the census tract level (?GEOID?), which are nested within Precincts (?Precinct?). Please see below my data structure. I used the code to sum population data per precinct: inters <- inters %>% group_by(Precinct) %>% mutate(TotalPop = sum(totpopE) ) However, said code produced too large sums because each census tract (?GEOID?) has multiple

I have several rasters that I want to do some calculations ,basically calculating the moving average. dir2 <- list.files("D:\\2010+2011", "*.bin", full.names = TRUE) saf=stack(dir2) movi <- overlay(stack(saf),fun=function(x) movingFun(x, fun=mean, n=3, na.rm=TRUE)) Error in .overlayList(x, fun = fun, filename = filename,

Hi R-Users, I have plotted a region whose polygon coordinates are given in shp format ED50 UTM (zone=30) ) using "readShapePoly" in library(maptools). Now I need to plot a set of points in that region (my.dataframe, with X and Y geographic coordinates), which have been read using GPS in Longitud-Latitud form (using WGS84 system), so I first need to convert these Longitud-Latitud data

Dear all, I have been trying to convert coordinates from longitude/latitude to utm but I got an error. As soon as the longitude coordinate is greater than 90, I get the folloowing error message: "error in pj_transform: latitude or longitude exceeded limits" Here is what I did: SP<-SpatialPoints(cbind(126.59,-14.30),proj4string=CRS("+proj=longlat"))

Hi, I am stuck in a looping problem. It might be an easy problem for experienced R users but I have been unable to do it. Any kind of help or advice will be great appreciated. I am creating a spatial grid and have a list of spatial points in a folder. I can read the spatial points separately and get a count of the points in the spatial grid in R. The output should be a table showing the counts of

What fun. I don't often stoke up the boiler and actually compile open-source code, but here I am at the beginning of just such a joyous and wondrous adventure. I'm running the latest level of AIX (V5.3...V4.3 has been obsolete for at least, what, 5 years) and using the IBM compiler (xlc is similarly obsolete, from back when the compiler came with the O/S...oh, for the halcyon days of

Hello, I have a problem writing a raster with the raster package. I have the raster object "mask" which has the following geoinformation: ---- R output > mask class : RasterLayer filename : E:/Daten/FE/HyMAP/Luxembourg_2010/Kehlen_Useldange/mosaik/LUX_LC_noOverlap_mask nrow : 5198 ncol : 2813 ncell : 14621974 min value : 0 max value : 255