Tag Info

Hot answers tagged

3

So, based on AndreJ's answer above, I whipped this up. cleanSpLinesForProjection <- function(w){ #We'll have to go through this one, line by line... slot(w, "lines") <- lapply(slot(worldmapLines, "lines"), function(x) { coords <- slot(x, "Lines")[[1]]@coords #get the rows with too large longitude values rIDX <- which(coords ...


3

Your data exceeds the 180°E limit in eastern russia and some pacific island. If you want a good looking map in degrees, you have to cut your source data at 179.9°E/W. See my answer here (though it deals with a pacific centered view) QGIS display world country shape files centered on pacific ocean using Robinson, Miller Cylindrical or other projection


2

If the degree coordinates should be in Helsinki (for which UTM zone 35 is correct), you should enter 24.9993393, 60.184716. Note that proj expects East - North order, with negative singns for Southern and western coodinates. Converted to EPSG:3067, the coordinates are 389042, 6673664.


1

You could use the GDALAutoCreateWarpedVRT function. This function creates a virtual raster dataset that represents the projected version of your raster. If you use this function, then you can read the properties of your image in the new coordinate system from the virtual raster's geotransform parameters.


1

Thanks for all of the suggestions. Over the weekend my boss re-imported the shapefiles that were giving me trouble. QGIS assigned them custom "user:100000, *generated CRS" to match the file CRS (EPSG: 3857). It looks like the shapefiles were arbitrarily assigned to WGS84 by one of us at some point last week. So I suppose the lesson learned here is: know ...


1

I was able to figure it out. I had to create a custom project and add it to geoserver. At first, I was trying to find pre-defined projection thinking they were complicated to add to the system but it was much easier than what I though. I took a previously defined user projection in geoserver and just copied and pasted parts of the projection of the layer to ...


1

You will not get lucky with predefined EPSG codes, because these are not defined for all local CRS in the world. gdalsrsinfo on the WKT definition returns: +proj=lcc +lat_1=0 +lat_2=0 +lat_0=46.66964837722222 +lon_0=-91.15277777777777 +x_0=228600.4575057151 +y_0=148551.4835661671 +ellps=GRS80 +units=us-ft +no_defs So you have to create a custom CRS to ...



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