Questions tagged [proj]

A library for performing conversions between cartographic projections.

106 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
5
votes
0answers
609 views

Fuller/Dymaxion in proj4/gdal/qgis?

I know this question has been asked a million times already, but I can't find many detailed responses on exactly what I want to do. As far as I can tell, the full C code for the fuller projection is ...
4
votes
0answers
67 views

Reproject WGS84 raster file to topocentric LTP-ENU using GDAL with a custom proj string

I have some raster files which are given in WGS84. As an example, here are some gdalinfo for one of those: Driver: GTiff/GeoTIFF Files: dem_r.tif Size is 2000, 911 Coordinate System is: GEOGCRS["...
4
votes
0answers
118 views

GIS proj pipeline

I would like to create a custom projection in QGIS 3.10. It is a local stereo and the reprojection between WGS84 and the local stereo uses proj pipeline. For example in cs2cs it is (which works well): ...
4
votes
1answer
1k views

Handling raster data in a “rotated pole projection” with R?

My question is about using R to transform raster data from a rotated pole projection to a regular latlon grid - and then plotting it. I am working with data from a numerical weather model, which is ...
4
votes
0answers
293 views

proj4js gives wrong result for EPSG:3375

Faced with very strange issue. Trying to do coordinate transformation from EPSG:4326 to EPSG:3375 using proj4js and getting wrong result. This is my code: var epsg4326 = '+proj=longlat +datum=WGS84 +...
4
votes
0answers
634 views

Creating Proj4 string for local coordinates

I am fairly new to coordinate system conversion, so am rather lost on the following problem. We have been using Proj4 to convert between various coordinate systems which is fine for known coordinate ...
4
votes
0answers
827 views

Creating custom CRS based on local project parameters in QGIS?

I'm having trouble creating custom proj4 CRS files to work with existing .shp files on local project coordinates. (Local in that the coordinates are relative only to themselves and not the 'real world'...
4
votes
0answers
296 views

Proj4 string for Hobo–Dyer projection?

I was fortunate to find an SVG of this projection at http://bl.ocks.org/mbostock/raw/4476487/ However, there appears to be little reference to this projection outside of wikipedia. Does anybody know ...
4
votes
0answers
2k views

How to set GDAL/cs2cs/proj to use NTv2 grid shift file when using EPSG codes?

I am using the GDAL API for Python for a script that takes only EPSG codes for inputs to specify source and target coordinate systems. While the cs2cs command seems to handle quite well ...
3
votes
0answers
54 views

Proj4 string missing datum information for SAGA outputs loaded into R

I've produced a set of topographic derivatives using both SAGA and ArcGIS. They are all in the same CRS. When I load them into R, however, the proj4 strings are different: >arcRaster class ...
3
votes
1answer
156 views

QGIS/proj errors with EPSG:2951

I've built various QGIS projects that use Canadian UTM and MTM CRS. They worked just fine in QGIS versions through 3.8 (based on proj4 and ?proj5). In QGIS windows version 3.10.0-3, things got borked ...
3
votes
0answers
519 views

Anaconda UnsatisfiableError installing QGIS

If I run conda install -c conda-forge qgis=2.18.2 from my Ubuntu bash environment, I obtain the following error: UnsatisfiableError: The following specifications were found to be in conflict: - ...
3
votes
0answers
249 views

Differences beetween proj4 string and WKT

What are the differences between proj4 strings and well-known text in terms of the ability to represent a projection, other than for WKT being more human readable?
2
votes
0answers
53 views

Visualizing resource dependencies and relationships underlying EPSG CRS

I'll go straight to the point; I remember long ago that I was quite shocked the first time I discovered that EPSG:4326 was kind of an "umbrella" definition on top of many others, which often ...
2
votes
0answers
25 views

Displaying (Correctly?) Projected PNG onto Leaflet Map

I am trying to overlay a PNG onto a Leaflet map. My GDAL commands: UL="-125.002000 52.013300" LL="-125.002000 22.013300" UR="-66.492200 52.013300" LR="-66.492200 22....
2
votes
0answers
30 views

Determine if a coordinate reference system is projected from Proj text

I'm looking for a reliable way to determine if a coordinate reference system is projected or not based on its Proj text. As examples I'm looking at EPSG:3857 and EPSG:4326 EPSG:3857: +proj=merc +a=...
2
votes
0answers
70 views

Did the syntax for proj4 change with QGIS 3.10?

In QGIS 2.8 I was able to successfully do a transformation using the following proj4 syntax : +proj=utm +zone=18 +ellps=clrk66 +datum=WGS84 +towgs84=-11.0584,163.8703,190.541,0,0,0,0 +units=m +...
2
votes
0answers
253 views

OpenLayer's graticule doesn't properly work

I am making a map of Krasnoyarskiy krai subject of Russia. I have defined a custom projection that should be used for this area using proj4 string. I need to show lat/long grid only for this territory....
2
votes
0answers
94 views

Why QGIS cannot georeference a Grib2 file with CRS (PROJ.4 and OGC WKT2)?

Situation: I have some .grib2 files that I downloaded from a (supposed to be) reliable source. When I do $ gdalsrsinfo file.grib2 I get: PROJ.4 : +proj=lcc +lat_0=38.5 +lon_0=262.5 +lat_1=38.5 +...
2
votes
0answers
29 views

Project point from one UTM zone to other

I am currently working on a project where we have to figure out the position of one point based on the known position of two (user-input) points. Now the problem is that the users could provide two ...
2
votes
0answers
39 views

How to run OSGeo4W shell for molodensky transform?

I'm a beginner with Proj, so I just downloaded proj, and I'm using the OSGeo4W Shell and when I try to run the example proj=molodensky a=6378160 rf=298.25 da=-23 df=-8.120449e-8 dx=-134 dy=-48 dz=...
2
votes
0answers
44 views

Create a Non-North/Rotated-Vertically User Coordinate System?

I need to plot 3d PostGIS data into vertical planes for symbolization and drawing vertical cross-section interpretations for mineral exploration models. I have written a short plpgsql function that ...
2
votes
0answers
39 views

Why is proj version dated in the future?

I compiled proj (also known as proj.4) from sources on Ubuntu (https://github.com/OSGeo/proj.4.git). I now wonder why is proj version dated in the future? E.g.: $ proj Rel. 6.0.0, March 1st,...
2
votes
0answers
593 views

No handlers could be found for logger “shapely.geos”

I have the following code to get area of polygon using shapely and pycrs. I was able to run this code successfully in a jupyter notebook but when I incorporate it into my main code, I had an error ...
2
votes
0answers
68 views

Can't open an image in CS63 in QGIS

I have an image in projection Pulkovo 1942 / CS63 zone X2. When I try to open it in QGIS, it fails to detect the projection and asks to select the projection manually. Does QGIS handle this family of ...
2
votes
0answers
175 views

Reproject vector layer in QGIS?

I want to warp a vector layer (set of polygons, but not a raster map) to fit a map. The problem is, the vector layer coordinates are from 0 to 360 and the basemap layer coordinates are from -180 to ...
2
votes
0answers
138 views

Countries cutting off at the edges of Mollweide projection

I have a custom CRS centering on the Pacific ocean with the Mollweide projection: +proj=moll +lon_0=-158 +datum=WGS84 +units=m +no_defs +lon_wrap=-158 My countries are split at the 22nd meridian, ...
2
votes
0answers
580 views

Converting Custom (local) .PRJ to .PROJ4

I have a Local PRJ file that I want to use in QGIS but how do I convert it to PROJ4? Details for the PRJ file as below (based of GDA94 / MGA Zone 53) JMG Authority: Custom Projection: ...
2
votes
0answers
724 views

Using PROJ.4 to convert from 3D Cartesian ECEF coordinates to Lat/Long coordinates with inverse Icosahedral Snyder Equal Area (ISEA) projection

I have a set of X, Y, Z coordinates in the Earth Centered Earth Fixed (ECEF) geocentric system. I need to convert them from ECEF coordinates into Latitude/Longitude coordinates using an inverse ...
2
votes
0answers
302 views

How to fix a shapefile reprojected from sinusoidal projection to lat long?

I have been work in convert the MODIS Tiles Grid from sinusoidal projection to lat long projection using python ogr2ogr. It doens't works well, the shapefile reprojected presents some lines that was ...
2
votes
0answers
265 views

Seeking example .PRJ file using WKT syntax for Concat_MT?

I've been digging into the WKT format definition found here:WKT Format with the aim of defining a PRJ with two chained affine transforms to go from a known projected system to a local grid. I've been ...
2
votes
0answers
59 views

mapnik proj4 tpers error

I'm working on a project where we need our tiles rendered with tilt perspective. I have installed and running openstreetmap server and it's working, but I need this projection to work: +proj=tpers +...
2
votes
0answers
372 views

“getting Unknown projection string error” serving EPSG:4326 tiles with open street maps mapnik renderd

The goal: I am attempting to generate map tiles with EPSG:4326 instead of web mercator by specifying EPSG:4326 in the mapnik style file srs attributes My current setup: I am running Ubuntu 14.14 LTS ...
2
votes
0answers
202 views

Projected Coordinate System issue using GDAL and qgis when defining x and y

I have a numpy.array which I would like to save as GeoTiff using GDAL in python. The array was created using the following PCS in mind: +proj=lcc +lat_1=48 +lat_2=48 +lat_0=48 +lon_0=0 +x_0=0 +y_0=...
2
votes
0answers
578 views

Converting from flat earth to WGS84 with python

I'm looking to move from basic 2d Cartesian coordinate system into WGS84. I know my starting coordinates, however all intermediate points are x/y/z points. of an origin 0,0,0 . Mathlab has a ...
2
votes
0answers
1k views

PyProj/Proj4 convert geographic to projected coordinates for Geostationary Data

I am trying to convert geographic coordinates (in degrees east/north lon lat) to projected coordinates in meters for geostationary full disc data (GVAR Goes East data from CLASS). To do this, I tried ...
2
votes
0answers
766 views

How can I fix badly specified geographic coordinate system with units in arc-seconds?

Check out this depth-to-rock soil data. It is linked from here. The spatial reference system is a latlon, but the units are arcseconds. I could not reason out how to override the SRS even after ...
2
votes
0answers
1k views

Projecting point in Leaflet from EPSG 900913 (same as 3857) to EPSG 2232 with proj4js returning unexpected results

I have a leaflet map from which I want to be able to query an API that accepts x y values in EPSG 2232 as inputs. Leaflet doesn't handle point transformations like OpenLayers does so I am trying to ...
2
votes
0answers
111 views

Proj4s.js Library Transformation Accuracy Documentation

I am using proj4s.js for transforming between coordinate systems for a GIS web application. Is there any documentation out there on transformations accuracy provided by the library?
1
vote
0answers
17 views

How to use a custom CRS in PROJ?

I want to change the CRS of a shapefile which uses EPSG:6257. Unfortunately, I am using Ubuntu 20.04 which comes with PROJ6. This particular CRS was only made available from PROJ7.2 onwards. So, I am ...
1
vote
0answers
45 views

SRID Transformation in PostGIS giving inaccurate transformation

I am trying to make a transformation in PostGIS from ESPG 4326 to ESPG 5247. The problem is that when the data is imported into ArcMap 10.7 it is clearly in the wrong place (by about 2.5km). I have ...
1
vote
0answers
62 views

Installing rgdal - cannot solve dependencies

I am trying to build a docker image based on osgeo/gdal:ubuntu-XXXXX. However it keeps failing when trying to install rgdal. I've tried this on gdal-3.1.1/proj-7.1.0 and gdal-3.0.3/proj-6.3.0, as well ...
1
vote
0answers
62 views

How correctly convert WGS84 to Pulkovo 1942 with height?

I have next python code: import pyproj import sys import logging from pyproj import Transformer from pyproj import CRS lat = sys.argv[2] lon = sys.argv[3] z = sys.argv[4] transg = Transformer....
1
vote
1answer
656 views

Steps for installing new GDAL version for Python in Ubuntu 18.04

If I'm installing the GDAL stock version from Ubuntu 18.04 then the installed version is 2.2 and I need at least 3.0.4 to better handle Sentinel 2 data. Is there a way to update GDAL to version 3.1.0?...
1
vote
0answers
67 views

Custom CRS: different results in PROJ6 and QGIS

I have created a custom CRS with PROJ6: echo 115.90428843055600 -1.92307681111111 126.014 | cct +proj=pipeline +step +proj=unitconvert +xy_in=deg +xy_out=rad +step +proj=tmerc +lat_0=-1....
1
vote
0answers
169 views

Updating PROJ GEOS and GDAL in R

I have been having trouble linking to GEOS 3.8.1, GDAL 3.0.4, PROJ 6.3.0. I have tried reinstalling sf, rgdal, raster, and rgeos packages. I still couldn't link to GEOS 3.8.1, GDAL 3.0.4, PROJ 6.3.0. ...
1
vote
0answers
91 views

Winkel II fails while reprojecting on the fly in QGIS 3.12

Currently, in QGIS 3.12 the Winkel Tripel projection working on the fly but Winkel II, which is also implemented, is not working. This is the fact or there is a bug? If I reproject dataset into ESRI:...
1
vote
0answers
53 views

Is there a piecewise-UTM coordinate system?

Is there an EPSG code, WKT or a proj4 string, for a single global piecewise-UTM coordinate system such as MGRS (the NATO military grid reference system)? I suppose this would work similar to the ...
1
vote
0answers
44 views

Creating custom PROJ projection from control points

I have a JPG raster, and I used QGis to georeference it with "control points", which are pixel coordinates linked to map coordinates, in EPSG:4326 (lon, lat). Here is an example of the control points ...
1
vote
0answers
94 views

gdalinfo returns PROJCS[“unnamed”] but QGIS returns the correct projection

I'm trying to understand why my gdalinfo call to a specific file returns "unnamed" for the PROJCS node. The call looks like this: ./gdalinfo /path/to/file/dem.tif. Output (notice the "unnamed"): ...