Re: [R-sig-Geo] Broken conversion from sp to spatstat ppp

David, iirc these conversion methods are in package maptools. On 06/18/2013 08:40 AM, d.g.rossi...@utwente.nl wrote: > Since upgrading to R 3.0.1 and the latest version of packages sp 1.0.9, gstat > 1.0.16, spatstat 1.31-13 the following conversion no longer works: > >> require(gstat) >> data(me

[R-sig-Geo] Broken conversion from sp to spatstat ppp

Since upgrading to R 3.0.1 and the latest version of packages sp 1.0.9, gstat 1.0.16, spatstat 1.31-13 the following conversion no longer works: > require(gstat) > data(meuse) > coordinates(meuse) <- ~x + y > meuse.ppp <- as(meuse, "ppp") Error message: Error in as(ds, "ppp") : no method or

Re: [R-sig-Geo] Create a mask with 1 deg latitudinal bands

How about this? library(raster) ## quick way to get 1-degree raster with your extents r <- crop(raster(), extent(-58,-48,-17,-2)) r2 <- setValues(r, rowFromCell(r, 1:ncell(r))) plot(r2) On Tue, Jun 18, 2013 at 1:10 AM, ppanday wrote: > Hi, I am trying to create a mask within a specified exten

Re: [R-sig-Geo] Lat Long --> UTM again

Hello, it was just a silly mistake: X and Y flipped. Thanks for the person who noticed it and for the time you took to give it a look. here goes the corrected fragment: #convert to sp object: coordinates(LatLong) <- ~ X + Y #longitude first > UTM <- spTransform(LatLong, CRS=CRS("+proj=utm +zone

Re: [R-sig-Geo] Lat Long --> UTM again

Hi Moshood Agba, thank you for the quick reply. i get the following message: Error: unexpected '=' in "UTM <- spTransform(LatLong, CRS("+proj=utm +zone=29 +ellps=WGS84 +datum=WGS84" +units=" i noticed that difference in the result when using the espg information in the code but i don't know wh

Re: [R-sig-Geo] Lat Long --> UTM again

Hi Joana, Try to modify your code as shown below and re run the script. #project using sp transform: UTM <- spTransform(LatLong, CRS("+proj=utm +zone=29 +ellps=WGS84 +datum=WGS84" +units=m +no_defs +towgs84=0,0,0")) Let me know if it works for you. Moshood On Mon, Jun 17, 2013 at 2:53 PM,

Re: [R-sig-Geo] Lat Long --> UTM again

Hi Joana: You may want to check out this website: http://georepository.com/crs_32729/WGS-84-UTM-zone-29S.html I think that you're getting the false northing and easting. Thanks, Erin -Original Message- From: r-sig-geo-boun...@r-project.org [mailto:r-sig-geo-boun...@r-project.org] On B

[R-sig-Geo] Lat Long --> UTM again

Hello, perhaps someone can help. i'm trying to convert coordinates from decimal degrees (WGS84) to UTM zone 29S (Portugal), but i have two problems: 1) i always get a wrong northing. I've read all messages regarding this subject and online documentation, which greatly helped, but i still couldn't

[R-sig-Geo] Anova between GLM and Spatial Eigenvector GLM

Dear list, Recently I used a Spatial Eigenvector Generalized Linear Model (SEVM-GLM) to understand whales distribution in Rio de Janeiro in spdep package and when I compared the SEVM-GLM with a "normal" GLM (as the example in page118 indicated) I had no differences (p>0.05) between them. What d

[R-sig-Geo] Create a mask with 1 deg latitudinal bands

Hi, I am trying to create a mask within a specified extent [ext <- extent(-58,-48,-17,-2)] that has 1 degree latitudinal band masks. So every 1 degree is a different number category. Eventually, I want to use it in zonal (as a mask) to obtain zonal statistics. I have a way that is a bit roundabout

Re: [R-sig-Geo] plotting counties of two states

Hi, Since you sent me your data (I'd still prefer dput() and you should use that in the future so you can send your data to the list): I don't see how your problem has anything to do with YNP. See below. On Mon, Jun 17, 2013 at 9:18 AM, Milan Sharma wrote: > Could you please run my code and see

Re: [R-sig-Geo] plotting counties of two states

Hi, On Sun, Jun 16, 2013 at 10:46 PM, Milan Sharma wrote: > Hi all, > As I already mentioned, this is the problem due to one extra county > "yellowstone national" in montana. I am trying to plot Montana and wyoming > together for further analysis.. At the last line, index does not match > betw

[R-sig-Geo] band/wavelength names

Dear friends, Is there a possibility to write/keep band names or wavelength names within an Envi *.hdr file when writing a RasterStack to a file? I only get a standard name output (band1, band2...) in the header file, althouhgh the original file contains wavelength numbers. I would like to avoid r

Re: [R-sig-Geo] STIDF subsetting problem

I don't know what you would like to modify: > library(spacetime) Loading required package: sp > data(air) > class(rural) [1] "STFDF" attr(,"package") [1] "spacetime" > dim(rural) space time variables 70 4383 1 > dim(rural[c(1,1:10),]) space time variables

Re: [R-sig-Geo] STIDF subsetting problem

Thanks for your interest, and timely fix, Edzer. Do you also think that some modifications can be done in the way STFDF and STSDF are subsetted? I understand that they are different than STIDF as length(rural@sp)!=nrow(rural@data) Servet -- View this message in context: http://r-sig-geo.273

Re: [R-sig-Geo] STIDF subsetting problem

I think this is a fair request, and changed the code accordingly. The updated code is on r-forge; the next CRAN version will support this. Also, this version has an stplot method for STIDF that shows attributes per time (slice), e.g. point pattern marks, the same way spplot does for multiple attri

Re: [R-sig-Geo] Help with creating a weight matrix for point dataset of 78, 000 observations

On Sun, 16 Jun 2013, Susan Gorelick wrote: Hi, I need to create a weight matrix in R for point dataset of 78,000 observations from a shapefile( from ArcGIS) which has both sets of longitudes and latitudes and X-Y coordinates (based on state plane feet). So far I projected the dataset in ESRI