Hot answers tagged

2

Perhaps you need to first look at how ifelse works. I get the same results when I use it "stand-alone" and within a call to raster::overlay. a <- rep(2, 5) b <- rep(1, 5) d <- c(2, NA, 2, NA, 2) library(raster) r <- raster(nrow=1, ncol=5) A <- setValues(r, a) B <- setValues(r, b) D <- setValues(r, d) s <- stack(A,B,D) ifelse(a==2 ...


1

Tx for the feedback, this is important for others who find this question in the future. Are you using WGS84 in your tilestache configuration? There is also a 'projected' attribute you might have to set. You are using VecTile class, I haven't played with that one yet, but you might be able to get it work (or atleast it will give you some clue/ideas to ...


1

I got the tiles to overlay correctly. The problem was in the re-projection done by both ArcMap and QGIS. When I was checking the reprojected shapefiles in ArcMap and QGIS, they were overlaying correctly and had the correct SRIDs. So I imported the shapefiles in WGS84 in PostgreSQL using the SRID4326 with shp2pgsql then used ST_Transfrom to reproject the ...


1

I guess your dataset about population density is the GEOSTAT grid dataset from Eurostat. I have already worked with it so I can give you the main steps. First, your GEOSTAT dataset should come with a shapefile (for me, it's called Grid_ETRS89_...). You have to import in QGIS both this shapefile and the .csv file that contains the population data. Then, you ...


1

Perhaps a bit late to the party, but for future reference: a simple Spatial Join will do the trick. Define a Merge rule to the concatenate the IDs of multiple features into a single field. Target features: Polyline layer Join features: Polygon layer Join type: One-to-One Field map: right-click a textual ID field of your Polygon layer and set Merge rule ...



Only top voted, non community-wiki answers of a minimum length are eligible