Questions tagged [proj]

A library for performing conversions between cartographic projections.

Filter by
Sorted by
Tagged with
1
vote
1answer
22 views

How to update CRS r object from PROJ4 to PROJ6

I have a script which contains several r objects which are CRSs that I use for transforming spatial objects. After updating R and many packages, there is a warning which wasn't there before. LCC <- ...
4
votes
0answers
51 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["...
1
vote
1answer
41 views

Adding HS2 SnakeGrid Projection to QGIS

How do I add the HS2 Snake Grid Projection to QGIS? I have the projection file HS2TN15_NTv2.GSB and I have added to the shared projections folder C:\Program Files\QGIS 3.10\share\proj. I also know ...
0
votes
0answers
58 views

Converting LCC crs to another format

I collected metadata about POIs which I want to visualize. The original data however makes use of a projected CRS ("LCC projection") which I want to transform to lat/lon (WGS84): "crs&...
1
vote
1answer
69 views

Unable to install proj related Python packages in a Linux alpine environment

I'm unable to install pyproj or any proj related Python packages such as geopandas, in an alpine dockerized environment, e.g. [python:3.9.1-alpine3.12][1]. Error is, e.g. for pyproj itself here: # pip ...
0
votes
0answers
22 views

Difference between transformation values in Proj4 npm and epsg.io

The values are not matching with other online tools for the conversion using Proj string and in the case of wkt format, the conversion using -ive value doesn't match. While using WKT string in epsg ...
1
vote
1answer
49 views

Why is this coordinate system assigned to this raster (R)?

I'm trying to verify that my raster has the correct coordinate reference system represented in its WKT, as I had to manually re-assign the datum after importing it (due to the PROJ6 migration). ...
2
votes
1answer
43 views

How to use Leaflet to display a map in Lambert Conformal Conic without warping to EPSG:3857

I would like to display a map in leaflet which is in Lambert Conformal Conic. The CRS of this map is: '+proj=lcc +lat_0=38 +lon_0=-100 +lat_1=45 +lat_2=33 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs' ...
0
votes
0answers
40 views

Understanding difference between osgeo4mac libraries and basic libraries (PROJ, GDAL, etc)

Is there a meaningful difference between osgeo-gdal and gdal that would make me choose the osgeo version over the basic version? The context here is that I'm trying to install grass on my mac (...
0
votes
1answer
32 views

build gdal3.2 encounter proj_api.h error: conflicting declaration ‘typedef struct projUV projUV’

How do I fix this error when making GDAL 3.2.0 (have also tried 3.0.4, 3.1.0) from source. The messages are: In file included from /usr/include/projects.h:284:0, from /usr/include/...
0
votes
0answers
22 views

Get projection parameters from Leaflet

I have a Leaflet map that is customized by an end-user by panning and zooming to an area of interest. I'd like to then take all of the relevant parameters (projection, zoom level, bounds, etc) and ...
0
votes
1answer
40 views

Data from EPSG dataset in proj format

In the EPSG dataset obtained from https://epsg.org/home.html where to find the values to fill up the field of +proj =" or +zone= " (to be in proj.4 string)from the database itself? Will the ...
2
votes
2answers
70 views

Test whether coordinate reference system has units of meters

Is there a simple way to test whether a spatial object's coordinate reference system expresses coordinates in units of meters? The snippet below shows my initial attempt at a function that provides ...
2
votes
0answers
51 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 ...
0
votes
1answer
40 views

Creating an oblique Mercator projection using proj and QGIS

I am trying to create a coordinate system within QGIS that's aligned with a great circle between two points. The points defining the great circle are (-44.0597,168.7192) and (-42.7365,171.4428) and, ...
1
vote
0answers
40 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
1answer
25 views

PostGIS unknown projection id in ST_Transform proj4 string

Using PostGIS I'm trying to transform a point coordinate from EPSG:31468 to EPSG:4326 using a proj4-string that I obtained from QGIS "Select Datum Transform" (is it one?): SQL is of the ...
3
votes
1answer
41 views

Having problem with PROJ string for custom Coordinate System

I was looking for an exhaustive reference on instructions supported by PROJ string, but couldn't find anything yet. Only this page in PROJ library docs - I guess that's all there is. I have quite a ...
1
vote
1answer
37 views

How to get rid of the transformation grid in st_transform?

I use the last versions of R -4.0.3- and sf -0.9.6- (thus relying on PROJ version 6.3.1). I want to transform a layer from one system of coordinates (EPSG:27572) to another (EPSG:2154). This can be ...
0
votes
0answers
99 views

Why does my rgdal discard datums?

I'm using rgdal to process rasters and shapefiles that contain CRS information [CH1903+ LV95 (EPSG2056) or CH1903 LV03 (EPSG:21781)]. For some reason part of the CRS information is always discarded (...
0
votes
0answers
19 views

Change Projection of a Web Mercator PCS to rotate it around center

I am using the GIS tool SharpMap together with ProjNet(4GeoAPI) as a tool to switch between CRS and transform coordinates (This information are just for context; my question is more abstract). The app ...
2
votes
1answer
91 views

Deprecated CRS (EPSG:2163) gets reinterpreted as another CRS (EPSG:9311) by GDAL

I recently noticed that all vector and raster files with CRS EPSG:2163 (US National Atlas Equal Area) are being automatically reinterpreted as EPSG:9311. When I load the layers into QGIS initially ...
1
vote
1answer
98 views

Clean way to convert from US feet to meters with lidR - LiDAR point cloud

I would like to convert this point cloud from US feet to meters. Metadata of the point cloud here. I firstly though to convert it properly by re-projecting to a metric CRS but this caused me a lot of ...
2
votes
1answer
91 views

Reproject using WKT for just the horizontal CRS

I'm trying to reproject LiDAR point clouds from California zone 2 (ftUS) using spTransform from sp R package to a metric CRS like EPSG 3857 (initial post here) and got non finite transformation ...
1
vote
1answer
76 views

“non finite transformation detected” using sp R

This is a part of a broader question already posted (Non finite transformation detected - lidR spTransform) but I am creating a new question for the sp related part of that. I'm trying to reproject ...
0
votes
1answer
32 views

zonal_stats returns PROJ error

I have been using zonal_stats from the rasterstats module for quite some time in Python. And I have never come across this error, which I tried to Google to find some solution with no succes. ...
7
votes
3answers
295 views

sp::proj4string equivalent in the PROJ6 framework?

I am confused about the shift from PROJ4 to PROJ6 in the (relatively) recent versions of the spatial packages of R (and gdal for that matter). I used to use sp::proj4string() to get the projection ...
2
votes
1answer
61 views

Using specific transformation in R sf::st_transform() to reproduce ESRI’s Project tool

I am reading in a CSV with latitude/longitude points (WGS 84, EPSG 4326), transforming into the NAD_1983_HARN_StatePlane_Maryland_FIPS_1900_Feet projection (WKID 2893) and saving as a feature class. I ...
2
votes
1answer
85 views

First call to transform() fails with inf, all subsequent calls are OK - what could be the reason?

I'm trying to get to the bottom of the following (weird) problem. When applying pyproj.Transformer transformation to a certain polygon, it fails the first time with (inf,inf) throughout. Second and ...
6
votes
1answer
828 views

How should I handle CRS properly after the major change in PROJ library?

As explained here, the PROJ library has changed a lot and one should stop using proj.4 strings. To me it is not clear, what I should do instead. This answer explains, how I could deal with the problem ...
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....
3
votes
2answers
257 views

What is the order of parameters in +towgs84?

What is the order of parameters in +towgs84 attribute within a proj string. I assume it is dx,dy,dz,rx,ry,rz,sf dx,dy,dz are expressed in meters rx,ry,rz are these in degrees/radians/arc seconds? sf ...
1
vote
1answer
73 views

Coordinate Transformation Difference between QGIS and Proj

I transformed coordinates on PyQgis and Proj command line and there is 3.8 meters difference between the transformations. I used the same EPSG codes. For Proj 6.3: echo 35.53 38.75 | cs2cs +init=epsg:...
0
votes
0answers
90 views

Unable to run Proj 'C' example

I use Proj 6.3 and I try to run Proj 'C' example from here. I tried it using Qt 5.14.2 on Ubuntu 20.04. I created a 'C' project from Qt Creator and tried the example. I get the errors: g++ -o proj_c ...
1
vote
1answer
65 views

Attached image is missing unit support in GDAL 3+ (and associated PROJ4)

This image below worked with the GeoTIFF tools in GDAL v2. When loading it in v3 I get "proj_uom_get_info_from_database: unit of measure not found". What unit of measure isn't found, and is ...
1
vote
1answer
121 views

get “Cannot find coordinate operations” in GDAL 3; works in GDAL 2

I get the following error when trying to use GDAL 3.1.2 for computing the LLA from my custom gnomonic projections. Am I setting it up wrong? Missing some information/files? Would it be a misconfigured ...
2
votes
2answers
136 views

Trying to set custom projection for shapefile in QGIS

I'm trying to open and use a dataset (in the form of a shapefile) (https://webarchive.iiasa.ac.at/Research/FOR/forest_cdrom/data/landscapes/landscapes.zip) with QGis. The readme of the dataset says ...
0
votes
1answer
47 views

Converting from Local Coordinate to lat/long and vice verse gives out different results

I am trying to convert a set of local co ordinates into lat/long and vice versa and I have the code for it taken from this site itself , the issue is When I do these conversions they are not resulting ...
1
vote
0answers
43 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 ...
0
votes
0answers
24 views

Is it possible to create a geometry from a proj4 string and Cartesian coordinates in PostGIS?

Considering I have a proj4 string such as: +proj=tmerc +lat_0=12 +lon_0=-12 +k=1 +x_0=0 +y_0=0 +datum=WGS84 +units=m Is it possible to create a geometry given X and Y offset (in meters)? For example ...
0
votes
1answer
33 views

Difference in area size calculation since update packages

On my mac I use brew to install PostGIS and all dependencies, including geos. Since month or so, I have a failing test in my Python code. A polygon that is entered in WGS84 (srid 4326), the converted ...
0
votes
1answer
73 views

Proj/PyProj Convert between ocentric and ographic latitudes

I am trying to use Pyproj to perform the conversion between ocentric and ographic latitudes. I can do this using the cct command line tool: echo 0 45 | cct -z0 +proj=geoc +a=3396190 +b=3376200 0....
4
votes
1answer
92 views

Custom CRS for historic “site” grid

I have a set of historic drawings which show a local site grid. I'd like to setup a custom CRS for this grid in QGIS so i can convert layers to/from this grid and MGA94 Zone 56 (EPSG:28356). ...
1
vote
3answers
82 views

Finding applied transformation method and accuracy information when reprojecting layer in QGIS

When I load two layers with different CRS I'm prompted to choose a transformation method and there is a column indicating the accuracy (see image 1). I would like to know where the information as part ...
0
votes
0answers
13 views

How use horner transformation in proj?

I try make horner transformation, but i get error # Record 0 TRANSFORMATION ERROR: x.xx y.yy Such an error even gives an example from the official website(https://proj.org/operations/transformations/...
1
vote
0answers
43 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
300 views

Transformation from CRS EPSG:3857 to ESRI:54009

I am trying to convert from EPSG:3857 to ESRI:54009 and facing some issues (pyproj version -2.6.1). ERROR 1: PROJ: proj_create_from_database: crs not found ERROR 1: PROJ: proj_create: unrecognized ...
0
votes
1answer
291 views

Reprojecting unequal number of latitude and longitude: pyproj.exceptions.ProjError: x, y, z, and time must be same size

I have a rectangular matrix of geographic climate data with a 1d row of latitude coordinates (size=834) and a 1d column with longitude coordinates (size=1115). They are both in unit of m (EPSG:32630) ...
1
vote
1answer
475 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?...
2
votes
2answers
85 views

Geodetic transformation with Mars data: 'ographic to 'ocentric latitude

Introduction I am working with a few raster datasets covering the south polar region of Mars. Some use a coordinate system with aerographic latitudes, which are equivalent in concept to geodetic ...

1
2 3 4 5
10