Questions tagged [proj]

A library for performing conversions between cartographic projections.

Filter by
Sorted by
Tagged with
0 votes
0 answers
25 views

Updating GDAL version in dockerfile giving error PROJ 6 symbols not found

I'm very new ubuntu, but I am using it in my dockerfile like below FROM ubuntu:bionic ENV ROOTDIR /usr/local/ ENV GDAL_VERSION 3.4.0 ENV OPENJPEG_VERSION 2.3.0 ENV CURL_CA_BUNDLE /etc/ssl/certs/ca-...
user avatar
1 vote
1 answer
18 views

Does GDAL 3.x/PROJ 7.x do some kind of bounds checking when doing a gdalwarp operation?

This is more of a question about how GDAL 3.x and PROJ 7.x work more than a question about how to do something. I have been experimenting with GDAL 3.4.2 and have noticed some things that are ...
user avatar
  • 431
0 votes
2 answers
27 views

why error msProcessProjection(): Projection library error?

I'm using mapserver in Centos 7. I'm using fgs from this (http://www.maptools.org/fgs/index.phtml?page=downloads.html). I'm using mapfile with epsg:4326, but I get error <ServiceException> ...
user avatar
0 votes
0 answers
60 views

How to remove banding artifact in R after reprojecting sf with st_transform

I have several global map layers I am plotting using ggplot() + geom_sf() + etc. All data is "+proj=longlat +datum=WGS84 +no_defs". My goal is to have Eurasia be a contiguous land mass. This ...
user avatar
  • 1
0 votes
0 answers
21 views

Append original shapefile projection definition to GeoJSON using Javascript

In our web application we upload shapefiles and convert them into geoJSON format. These shapefiles will be in different projections and I want to be able to add the original projection of the ...
user avatar
1 vote
0 answers
44 views

Reprojecting raster or vector CRS and extracting raster values based on vector boundaries

I am working with a EUR-11 CORDEX dataset of daily temperature values (nc file) from ESGF database and EUROSTAT NUTS data (from here: https://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/...
user avatar
  • 11
0 votes
1 answer
48 views

GeoTools: Geo Transform incorrect

I am transforming a coordinate between Lambert Conformal and WGS84. The values returned by GeoTools do not seem correct. From: PROJCS["unnamed", GEOGCS["Coordinate System imported ...
user avatar
2 votes
1 answer
64 views

Running osm2pgsql 1.4+ with custom projection, proj 8.2

first some context... I have a custom projection, here is the projection string: <5514> +proj=krovak +lat_0=49.5 +lon_0=24.83333333333333 +alpha=30.28813972222222 +k=0.9999 +x_0=0 +y_0=0 +ellps=...
user avatar
  • 21
0 votes
1 answer
31 views

QGIS PROJ custom CRS for only shifting X and Y

I'm not very fluent with CRS and custom CRS but I'm only want to shift X and Y WGS84/32N coordinates. For example: with point A(115633.49, 57104.49) of my maps, in EPGS:32632 would be A(1112126.15, ...
user avatar
1 vote
1 answer
157 views

Convert local coordinates to WGS84 in C#

Some people working on a the construction site use a local custom coordinate system. I have been given a proj4 string like this: +proj=omerc +lat_0=46.325454996 +lonc=7.99050231 +alpha=60.4026421358 +...
user avatar
  • 133
2 votes
1 answer
78 views

Proj pipeline support in PostGIS

I would like to use proj pipeline in PostGIS ST_Transform function. My sample SRS is as follows: +proj=pipeline +step +inv +proj=unitconvert +xy_in=m +xy_out=us-ft +z_in=m +z_out=us-ft +step +inv +...
user avatar
0 votes
1 answer
34 views

Why does this suggested change in a Python pyproj initalisation future warning result in an error?

I Recently updated to pyproj v2.6. and started getting future warnings when initalisating a Proj class like so: Proj("+init=world:bng") Proj("+init=EPSG:22770") Both of these ...
user avatar
  • 143
2 votes
0 answers
69 views

Wrong transformation between British National Grid and WGS84

I have access to a remote DB with the following specifications: PostgreSQL 11.12 POSTGIS="2.5.2 r17328" [EXTENSION] PGSQL="110" GEOS="3.7.0-CAPI-1.11.0 673b9939" PROJ=&...
user avatar
  • 277
1 vote
0 answers
68 views

Safe upgrade of PROJ for PostGIS?

Is there a safe way to install fresh PROJ, for example v8.2, in a Debian Linux server, for use with (for example) PostGIS v3? The installation of PROJ was implicit in the PostGIS's, I never do ...
user avatar
  • 2,084
0 votes
0 answers
54 views

gdal.translate via os.system/subprocess

I want to run a gdal.translate command via os.system() or subprocess.run() inside a script. I tested the following ways of running the command: Directly in the command-line: Starting python from the ...
user avatar
  • 127
1 vote
0 answers
50 views

QGIS cannot find transformation grid

When I open an existing QGIS project I get an error bar at the top of the mapping window stating that some data cannot be projected (see below). I can click on the details button and I get the ...
user avatar
  • 7,172
3 votes
2 answers
93 views

Convert geodetic to topocentric coordinates in pyproj

I'm trying to convert between topocentric coordinates and geographic coordinates with pyproj. First, I made a topocentric CRS using the proj4 string from the example here. Then, I made a geodetic CRS ...
user avatar
  • 31
0 votes
1 answer
49 views

Does repeated coordinate re-projection cause a loss in accuracy?

If I take a set of GIS data (shapes, points, lines, etc) and convert the projection many times, will there be a loss in accuracy in the data? The reason I ask is that I have a dataset from OSM which ...
user avatar
1 vote
0 answers
41 views

Creating custom CRS by rotating EPSG:27700 OSGB (for mine grid)?

I am trying to georeference a mine plan, but it uses a custom CRS based on the following: I gather I just need to rotate the OSGB grid system by 25 17'59" ? is this just a case of modifying the ...
user avatar
2 votes
0 answers
63 views

Use Mollweide projection in Leaflet, for a Geoserver WMS layer

I did a world map in Geoserver, showing in Leaflet in the default projection, which is EPSG:900913 (aka GOOGLE Mercator) but I would like to use the Mollweide projection instead. I successfully added ...
user avatar
2 votes
1 answer
386 views

Fixing mismatch between PROJ version and proj.db?

When transforming the CRS in QGIS, I get the following error: No transform available between EPSG:32722 and Custom CRS (or any CRS). When I click on Details, I read the following message: No ...
user avatar
  • 431
1 vote
1 answer
148 views

Can we construct WKT2 from PROJ4String correctly?

There is the new WKT2, with great features, but there is a lot of data saved in WKT/PROJ4String too. The docs I read say the PROJ4String doesn't have enough info for a complete reference for ...
user avatar
  • 21
1 vote
0 answers
37 views

Reprojecting WGS84 to Winkel-Tripel with new prime meridian breaks polygons

I've been trying to re-project a rather large shapefile of the ocean's polygons (found here) to a Winkel-Tripel using a prime meridian of -96 (US) rather than the standard 0. However, I'm noticing ...
user avatar
3 votes
1 answer
92 views

Coordinate transformation with NTv2 and cs2cs works in console but not in GDAL

Our Java application uses GDAL (v3.2.0) and I am trying to convert some coordinates from EPSG:31467 (DHDN GK3) to EPSG:25832 (ETRS89 UTM32). For this I downloaded the grid shift file BETA2007.gsb and ...
user avatar
5 votes
0 answers
74 views

How to use a Spilhaus/ESRI:54099 projection in QGIS

I am trying to make a map using a Spilhaus projection, similar to the one implemented by Esri here. QGIS has ESRI:54099 WGS_1984_Spilhaus_Ocean_Map_in_Square in its list of available projections, ...
user avatar
  • 333
0 votes
1 answer
55 views

Proj4 string for USGS Bulletin 1532

I'm using an application that uses equations from USGS Bulletin 1532 (see source code - line 295: http://docs.ros.org/en/kinetic/api/robot_localization/html/api/navsat__conversions_8h_source.html#...
user avatar
  • 81
0 votes
0 answers
32 views

How to find scale error for map projection?

I've just upgraded from proj version 7.1 to 8.2, and some code has stopped working. I used to use proj to determine the map projection scale error at a point with the following code: echo "-55....
user avatar
  • 1,588
4 votes
2 answers
146 views

How to use ISEA projection?

The projection named ISEA - Icosahedral Snyder Equal Area has a (supposedly) good implementation in the PROJ Library, proj.org/isea, that is used by PostGIS, QGIS, etc. The example is agnostic - was ...
user avatar
  • 2,084
2 votes
1 answer
51 views

Peirce Quincuncial Projection in QGIS

I'd like to use a Peirce quincuncial projection in QGIS to produce a whole sphere map, but struggling. (To see some nice examples of the projection:- https://github.com/cspersonal/peirce-quincuncial-...
user avatar
  • 31
0 votes
0 answers
140 views

Defining projection of Polar stereographic data in R via raster package

I'm having a hard time reading the following data into R: https://resources.marine.copernicus.eu/product-detail/ARCTIC_MULTIYEAR_BGC_002_005/INFORMATION It uses a polar stereographic projection, and ...
user avatar
0 votes
0 answers
108 views

Removing all spatially related software from Ubuntu

I managed to mess up my installation of gdal, proj, qgis and sf. I am not a hundred percent sure what I did wrong but I earlier today tried to update my system gdal installation by using the ubuntugis ...
user avatar
  • 707
0 votes
1 answer
47 views

Why is PROJ not figuring out a CRS?

I am trying to use gdal_translate to convert a WCS source to a tif target. I am trying to achieve this using a WCS XML service description file. The gdal_translate call looks like this "C:\...
user avatar
  • 608
2 votes
1 answer
99 views

Running proj/cs2cs command in Powershell

I'm reading a book where it says to run the following in the commandline, I think he used Linux. Therefore the command doesn't run in Windows Powershell cs2cs +proj=latlong +datum=NAD83 +to +proj=utm +...
user avatar
  • 1,159
1 vote
0 answers
137 views

Undefined reference to `proj_context_create' while run PROJ example

I have installed Proj 6.3.1 and I tried to run a Proj example from https://proj.org/development/quickstart.html. I implemented it with CLion on Ubuntu 20.04. And I got many undefined reference errors ...
user avatar
  • 11
0 votes
0 answers
102 views

Proj Error when trying to save data in PostGIS database from geodjango

I want to store my raster image(tiff) in my PostGIS database I have a Django model ndvitime and I want to store my tiff image from earth engine to my database. The error arises when I want to save the ...
user avatar
3 votes
0 answers
889 views

QGIS: No transform available error

I installed QGIS 3.16 on Ubuntu-20.04. Details on 'About QGIS': QGIS version 3.16.12-Hannover QGIS code revision be6cb030ac Compiled against Qt 5.12.8 Running against Qt 5.12.8 Compiled against ...
user avatar
  • 31
5 votes
0 answers
63 views

Function and meaning of parameter lat_0 in Transverse Mercator Projection in PROJ String

One may specify a custom Transverse Mercator in PROJ library string with something like: proj +proj=tmerc +lat_0=-25 +lon_0=-49 For me the function of the parameter lon_0 is clear, is the Central ...
user avatar
2 votes
0 answers
68 views

Change vertical unit in a LAS/Z file with PDAL

I am using pdal translate to filter and transform some lidar files. One issue we have are the vertical units. Sometimes they are correctly specified, sometimes they are not, and sometimes they are ...
user avatar
  • 1,214
0 votes
0 answers
78 views

Basemap npstere to cartopy NorthPolarStereo CRS

Somebody has given me the following code to convert north-polar stereographic (x,y) coordinates to WGS84 (lon,lats): m = Basemap(projection='npstere',boundinglat=60,lon_0=0,resolution='l') lon,lat = m(...
user avatar
2 votes
1 answer
107 views

Checking whether CRS is geographic coordinate system or projected coordinate system using R

I am writing a generic script that would calculate areas covered by different classes in a categorical raster in R. I would like to provide two functions; one for use when the CRS of the raster is ...
user avatar
3 votes
0 answers
923 views

"Failed to process SRS definition" error when reprojecting vector layer via ogr2ogr

I'm using GDAL 3.3.1 with all the necessary variables already setup. I'm reprojecting a vector layer using ogr2ogr ogr2ogr -t_srs EPSG:32651 test.shp new_layer.shp and an error says ERROR 1: PROJ: ...
user avatar
  • 1,159
1 vote
0 answers
126 views

Converting from Lat Long Altitude to Geocentric to Topocentric using modern PROJ C++

I'm trying to use Proj to perform lat long alt to geocentric to topocentric transforms. I've been able to do this fine in python, but the API pyproj uses underneath is apparently very old. I want to ...
user avatar
  • 111
2 votes
1 answer
47 views

OGRSpatialReference works with proj4 string in GDAL 2 but not GDAL 3

I have the following proj string: +pm=greenwich +ellps=GRS80 +a=6.37814e+06 +rf=298.257 +proj=tmerc +units=m +to_meter=1 +lat_0=42.5 +x_0=500000 +y_0=0 +lon_0=-72.5 +k_0=0.999964 +vunits=m and I have ...
user avatar
1 vote
0 answers
37 views

Custom Projection Definition - Excel Formula to ArcGIS Pro

I have been given an Excel formula for converting a custom projection coordinate point to a known defined system (state plane AZ East) I'm trying to get create a custom projection definition in ArcGIS ...
user avatar
3 votes
1 answer
91 views

Are some PROJ.4 projections missing in QGIS?

I am using QGIS 3.16.8 on Windows 10 (64). I have tried to apply the Apian II (aka Arago) projection with: +proj=apian2 +lon_0=0 +lat_0=0 +x_0=0 +y_0=0 +a=6371000 +b=6371000 +units=m +no_defs but the ...
user avatar
1 vote
1 answer
63 views

Reproject shapefile between two user projections

I need to tranform a set of geometries that have been drawn in the wrong projection system in an other projection system. I tried to do this using this command : ogr2ogr -s_srs bad.prj -t_srs good.prj ...
user avatar
  • 1,259
2 votes
1 answer
197 views

Converting to a CRS and the back to the original CRS, not the same original coordinates

I am using Proj4Net to convert between coordinate systems. I've noticed that when I convert from one coordinate systems (CRS) and them back to the original CRS, the coordinate Y values have changed ...
user avatar
  • 163
1 vote
0 answers
100 views

Helmert transformation in pyproj

I have a set of coordinate points projected in Web Mercator (EPSG:3857) and a raster specific coordinate system and have calculated the translation, rotation and scale for a Helmert transformation (P' ...
user avatar
2 votes
0 answers
111 views

How do I dump towgs84 conversions for PROJ strings?

Building off of the awesome answer by mikewatt here: How do I dump PROJ EPSG database to PROJ4JS-compatible PROJ strings? How would one dump the towgs84 conversion from the PROJ database for a given ...
user avatar
  • 2,609
1 vote
3 answers
263 views

Transforming EPSG:26918 to EPSG:4236 in GeoPandas

I am trying to transform a shapefile of polygons projected on CRS: EPSG:26918 to EPSG:4326 so that I can perform a spatial join in GeoPandas over two shapefiles. My previous code allowed me to simply ...
user avatar
  • 11

1
2 3 4 5
11