Questions tagged [proj]

A library for performing conversions between cartographic projections.

Filter by
Sorted by
Tagged with
1
vote
0answers
15 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 ...
0
votes
0answers
7 views

CDO without proj and hdf libraries (Regrid HDF file AMSR-2)

I am trying to use CDO for AMSR-2 data analysis, but when I installed CDO 1.9.9, it does not include Proj and HDF library. Even when I installed Proj and HDF libraries in Ubuntu, the CDO still unable ...
2
votes
1answer
31 views

Howto list EPSG codes as proj.4 strings from the new Proj 6 sqlite3 proj.db

I want to list all EPSG codes from the new PROJ.6 tool in form of proj.4 strings. Is there a tool / SQL script to replicate the former /usr/share/proj/epsg file from the SQLite3 proj.db?
0
votes
0answers
74 views

Georeferencing map with unknown projection and grid

It is another question about georeferencing a map with unknown projection and grid. You can find map here. As you can see it is already georeferenced, but poorly, I think. Projection used is ESRI:...
0
votes
1answer
38 views

Use custom projection in MapProxy

I'm trying to use the proj_data_dir config setting in MapProxy to add support for a custom projection. So far I've been unsuccessful and suspect I'm almost there, but not quite... I've created a file ...
0
votes
1answer
31 views

Using ArcMap georeferencing in AutoCAD

I am working on a project in ArcMap. I added an auto cad .dwg to my map, and in ArcMap, and I Georeferenced this .dwg file. So now the DWG drawing is moved to be exactly on my drawing in ArcMap. I ...
0
votes
0answers
24 views

Use American Polyconic Projection North American Datum 1927 with GDAL and QGIS

I am trying to use modern USGS DEM data with older USGS maps. The first step, as far as I am aware, is to reproject the contemporary DEM into the same projection as the old USGS map. The old USGS map ...
1
vote
1answer
295 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
66 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["...
2
votes
1answer
78 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
61 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
174 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
24 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
82 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). ...
1
vote
1answer
56 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
58 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
40 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
24 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
43 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
77 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
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 ...
0
votes
1answer
43 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
44 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
34 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
45 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
38 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
147 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
24 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
108 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
122 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
112 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
2answers
104 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
43 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
381 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
100 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
101 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
1k 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
325 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
76 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
119 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
86 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
165 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
141 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
61 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
25 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
34 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
94 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
112 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
2 3 4 5
10