search for: utm

Hi, is there a function in R to convert data read with read.shape and which is originally in UTM coordinates into longitude / latitude coordinates? I found the convUL() function from the PBSmapping package but I have no idea how I could apply that to the read.shape object. Many thanks, Werner __________________________________________________ Do sragenden Schutz gegen Massenmails. http:/...

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")) coor...

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...

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 <- spTransform(study_area, CRS("+proj=utm +zone=19...

Hi R helpers, I tried to convert a list of LatLong coordinates (in DD format) into UTM zone 11 NAD 27. I first tried this from PBSmapping: library(PBSmapping) LatLong<-cbind(c(56.85359, 56.85478),c(-118.4109, -118.4035)) colnames(LatLong)<-c("X","Y") attr(LatLong, "projection") <- "UTM" attr(LatLong, "zone") <- 11 UTM&lt...

Dear all: I have a shapefile in UTM coordinate system and I would like to transform it to Lat/Log coordinates (WSG84). The package PBSmapping contains function convUL to transform between the two coordinate systems when data is in the form of a data frame with attributes specifying the coordinate system. However, when shapefiles are...

...d().? I have increased the memory to 4GB, but am being told I can't allocate a vector of size 240 Kb.? I am running R 2.11.1 on a Dell Optiplex 760 with Windows XP.? I have pasted the error message and summaries of the objects below.? Thanks for your help.? Tim >???? xyz<-cbind(hi.to.utm,z=b.depth$z) Error: cannot allocate vector of size 240 Kb > memory.limit() [1] 4000 > memory.size() [1] 1971.68 > summary(hi.to.utm) Object of class SpatialPoints Coordinates: ??????? min?????? max x? 708745.5? 923406.7 y 2046153.1 2327910.9 Is projected: TRUE proj4string : [+proj=utm +...

Hi: Is there any function in R which can convert the long/lat to UTM(Universal Transverse Mercator)? There are quite a few converters on Internet. However, the interface is designed as input->output which I can not convert lots of locations at the same time. Another question is whether there is a function in R which can tell the time zone from the location...

Hi, I'm trying to use transform my mercator locations into utm but I'm doing something wrong because only x is transformed ... see > xy.sp SpatialPoints: loni lati [1,] -8.85 38.16 [2,] -9.19 37.99 [3,] -9.11 37.97 [4,] -9.06 38.15 [5,] -9.03 37.87 [6,] -9.14 37.81 [7,] -9.09 37.70 [8,] -8.95 37.45 [9,] -9.17 37.37 [10,] -9...

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

...w if there are in R already loaded 2 different line shapefiles with projection then i will see if i can replicate this problem. So, first file: buff class : SpatialLinesDataFrame features : 1 extent : 560525.5, 561302.4, 4698908, 4701110 (xmin, xmax, ymin, ymax) crs : +proj=utm +zone=16 +datum=WGS84 +units=m +no_defs variables : 2 names : ID, buff_dist value : 2, 135 second file: pln class : SpatialLinesDataFrame features : 1 extent : 560615.6, 560705.4, 4698905, 4699180 (xmin, xmax, ymin, ymax) crs : +proj=utm +zone=16 +datum...

...you should not use latitude and longitude...can anyone expand on this a little for me, and what would be the best coordinate system to use? I have my data in a geographic coordinate system, WGS84, decimal degrees....is this the wrong format for such analyses? I have also converted my data in the UTM projection and so have it in metres....(ranging from 480,000 to 550,000 E and 170,000 to 230,000 N). If I was to use the UTM coordinates, should I be using the actual coordinates in metres, or should I convert this into an arbitrary coordinate system (i.e. from 0 - 1) somehow? I have noticed...

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) myTiff = 'L5014028_02819840617...

Hello I'm working with a large hydrological data set stored in a netCDF format. The file stores x and y coordinates in the UTM projected coordinate system, yet when I use image to graphically display the z variable, the image is distorted in the sense that it does not plot the map in the correct spatial organization. I'm wondering if I need to define the projection of the netCDF file with rgdal or proj4 routines first...

...Consider the following toy example: ### Begin Example ### library( PBSmapping ) # Define some EventData events <- as.EventData( read.table( textConnection( 'EID X Y Color 1 494 1494 red 2 497 1497 blue 3 500 1500 green 4 503 1503 yellow' ), header=TRUE, strings=FALSE ), proj='UTM', zone=10 ) par( mfrow=c(3,1) ) # Plot the events with plot limits large enough to show # the full extent of all the symbols plotPoints( events, pch=16, cex=5, col=events$Color, xlim=c(490,508), ylim=c(1490,1508), proj=TRUE ) with( events, text( X, Y, toupper( substr( Color, 1, 1 ) ), f...

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 interes...

I'm pretty new to KVM and have a KVM CentOS 7.1 hypervisor running a few VMs. I'm moving all my VMs from an ESXi host as I want to use KVM in future. Most of my VMs are working except for one which is running a Sophos UTM router (Sophos UTM is similar to products like pfSense <https://www.pfsense.org/>, Smoothwall <http://www.smoothwall.org/> etc). The host has 3 physical NICs which are configured on the Sophos VM as: 1) LAN (fixed IP) 2) DMZ (fixed IP) 3) WAN (which is directly plugged into a cable mod...

...about spatial join on the minimum distance between 2 sets of points with coordinates and attributes in 2 different data frames. Simon Knapp sent code to do it when calculating distance on a sphere using lat, long coordinates and I've change his code to use Euclidian distances since my data had UTM coordinates. Typically one data frame has around 30 000 points and the classification data frame has around 4000 points, and the aim is to add to each point from the first data frame all the attributes from the second data frame of the point that is closest to it. On my PC (Dell, OptiPlex GX620...

Hola Ángel: Yo creo que tendrías que establecer el sistema de coordenadas de referencia de tu objeto raster antes de salvarlo como GTiff. Algo así: crs(r1) <-"+proj=utm +zone=14 +datum=WGS84" Saludos, Marcelino. El 20/02/2020 a las 1:42, Angel Cervantes escribió: > Hola a todos, quisiera pedirles su ayuda. Estoy tratando de crear un raster a partir de una tabla de datos que contiene coordenadas UTM. Determino los renglones y columnas que debe de tener...

...om_polygon(data = base_shape, aes(x = long, y = lat, group = group, fill = valor), color = "black", size = 0.25) + scale_fill_distiller(palette = "Reds", guide="colourbar") + theme_nothing(legend = TRUE) + geom_point(data = base_puntos, aes(x = long.agencia.utm, y =lat.agencia.utm, colour="Agencia"), size = 3.05, pch = 20, fill = "Blue")+ labs(title=paste(base_titulos,"\n(por provincias)"),fill="Temperatura",colour=NULL)+ theme(title = element_text(colour="black", size = 16, face = "bold"...