Re: [R-sig-Geo] Spatial Regression: Moving window regression (MWR) or GWR adaptation...?

I find this topic very interesting as well. Juan, are you referring to R to perform your analysis? Clint, a Poisson point process takes into account the time that events occur, doesn't it? Juan did mention nothing about the time if I'm not wrong, does it matter? On Jul 31, 2013 11:19 PM, "Clin

Re: [R-sig-Geo] Spatial Regression: Moving window regression (MWR) or GWR adaptation...?

Isn't this amenable to a Poisson point process ? See Konisky, et al, 2013, "Strategic Placement of Air Polluters: An Application of Point Pattern Models". Clint BowmanINTERNET: cl...@ecy.wa.gov Air Quality Modeler

[R-sig-Geo] Spatial Regression: Moving window regression (MWR) or GWR adaptation...?

Hi I'm trying to make a model to estimate the number of roberys. I know the address of these people and the point on the map where they live. I have a sample of all the country, in different cities, different zones. The variable can be set to 0,1,2 .,and it is distributed as a poisson. Right n

Re: [R-sig-Geo] readOGR -- Dropping null geometries -- Possible issue with path as a point

Thanks for the quick response! Here is the proper link if you want it. http://www.spc.noaa.gov/gis/svrgis/zipped/tornado_tracks.zip Yes, the 2012 tracks have a similar behavior. One more question... I tried setting dropNULLGeometries=FALSE but I still seem to only get the non null geometries,

Re: [R-sig-Geo] readOGR -- Dropping null geometries -- Possible issue with path as a point

The data you link to are not those you refer to: > tornados.rg <- readOGR(".", "tornado") OGR data source with driver: ESRI Shapefile Source: ".", layer: "tornado" with 56221 features and 28 fields Feature type: wkbPoint with 2 dimensions that source does not provide any called "tornado_tracks".

[R-sig-Geo] readOGR -- Dropping null geometries -- Possible issue with path as a point

Greetings! I am attempting to read a shapefile with readOGR and I'm getting the message 'Dropping null geometries'. From what I have been able to discern the null geometries seem to be associated with lines where the begin and end lat/lon are the same. So effectively the line is really just a

Re: [R-sig-Geo] raster::stack() help

Thanks for the advice Robert. I'll give those a try if I continue to have problems. For now my code is running again. -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/raster-stack-help-tp7584298p7584313.html Sent from the R-sig-geo mailing list archive at Nabble.com. _

[R-sig-Geo] Script to get polygons area in square meters

Hello, I am really new to GIS and I am posting here to check if what I learned so far is right. My goal is to measure areas of polygons (representing species distributions) in square meters. I'm going to use data available as shapefiles on IUCN Red List website. To read the polygon into R: libra

[R-sig-Geo] cokriging for large data

Hello every one, I estimate the regional snow water equivalent using cokriging ( my data set has more than 70,000 data points). All the things run finely and quickly. However, once I run the predict.gstat function it takes much time, but it doesn't work and have no results. Would any one plea

[R-sig-Geo] cokriging for large data

Hello every one, I estimate the regional snow water equivalent using cokriging ( my data set has more than 70,000 data points). All the things run finely and quickly. However, once I run the predict.gstat function it takes much time, but it doesn't work and have no results. Would any one plea

Re: [R-sig-Geo] gstat: Cokriging with unusual neighbourhood choice - is this possible?

Dear Edzer Thank you greatly for your reply. I have attempted to implement the solution you suggested, but either it's not what I was looking for, or more likely, I'm missing something. Currently the relevant code passage looks like this: ## m.grid=data.frame(x=outgrid$V1,y=o

[R-sig-Geo] Walk around the block

Hi Everybody I have a SpatialPolygonsDataFrame with residential stands in a town. I have created a SpatialPolygonsDataFrame with the street blocks by using poly2nb in the spdep package. My sampling procedure (which was selected for ease of implementation by an interviewer) is to take a random star

[R-sig-Geo] projection issue - can one rotate a projection?

Hi all My colleague, a physical oceanographer, and I are having problems getting his geographic layers into the correct projection (after importing as a netcdf file), and we are hoping someone can help us. We are using the Raster package. My colleague obtains a Dutch Rijksdriehoek in the first

Re: [R-sig-Geo] how to crete a *.prj file from a Raster object

And even: require(raster) require(rgdal) r <- raster(list(x = 1:nrow(volcano), y = 1:ncol(volcano), z = volcano), crs = "+proj=laea +lat_0=-50 +datum=WGS84") writeRaster(r, "file.asc") showWKT(proj4string(r), file="file.prj") # which also shows the WKT version on the console Roger -- Roger Biva

Re: [R-sig-Geo] how to crete a *.prj file from a Raster object

One way is to write a dummy shapefile, and rename the .prj to suit. To do that you could: require(raster) require(rgdal) r <- raster(list(x = 1:nrow(volcano), y = 1:ncol(volcano), z = volcano), crs = "+proj=laea +lat_0=-50 +datum=WGS84") writeRaster(r, "file.asc") dummy <- as(extent(r), "SpatialPo

Re: [R-sig-Geo] how to crete a *.prj file from a Raster object

On 07/30/2013 03:53 PM, Emanuele Cordano wrote: Dear list, I have a RasterLayer object with CRS. I would like to write it in a *.asc ESRI ASCII file with writeRaster. I did it , I got the *.asc file but no related *.prj file was written. Is there a quick method to create directly the *.prj file