Questions tagged [proj]

A library for performing conversions between cartographic projections.

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

Which temporary projection is most suitable for calculating a buffer around a route

To create a buffer around route (for purpose of querying around a route), I figured I need to: convert route points from geographic to planar coordinates calculate a buffer convert result back to ...
  • 111
0 votes
1 answer
25 views

Reproducing SWEREF99 to ITRF2014 Transformation using proj.org cct

The Swedish Geodetic Agency Lantmateriet provides a document (see references) with a procedure for transforming coordinates between ITRF2014 and SWEREF99. The document also provides and example of the ...
  • 111
1 vote
0 answers
12 views

PDAL +init use deprecated but what is the alternative

With PDAL I was trying to using a proj4 string like +init=epsg:32615 +geoidgrids=C:\vdatum_GEOID12A\vdatum\core\geoid12a\g2012au7.gtx based on an example given on the PDAL website https://pdal.io/en/...
1 vote
0 answers
38 views

Transform latitude and longitude to Transverse Mercator

I would like to transform point in WGS84 lat lon to Transverse Mercator: +proj=tmerc +lat_0=0 +lon_0=114 +k=1 +x_0=-49099.65399999998 +y_0=-2544749.4072 +ellps=WGS84 +units=m +no_defs I could easily ...
  • 103
1 vote
1 answer
33 views

Proj R_A discrepancy with computed authalic radius

Using Pyproj, I attempted to create a spherical CRS using the authalic radius from WGS84: import pyproj wgs84_spherical = pyproj.CRS.from_proj4('+proj=latlon +ellps=sphere +R_A=WGS84') print(...
1 vote
1 answer
33 views

Vertical grid shift with PROJ 9

I want to perform a vertical grid shift using PROJ 9.I have never used PROJ syntax and I am having some issues to have it right. I have checked documentation and this previous post from which I got a ...
0 votes
0 answers
23 views

GDALWarp Reprojecting raster extent from (0, 360) to (-180, 180) filling left hemisphere with blank data

Understanding that this question has been asked a few times (Reprojecting raster from 0 360 to -180 180 with cutting 180 meridian using gdalwarp seems to be the canonical one), I've found that any ...
1 vote
0 answers
13 views

How to get CS2CS to recognize a 7-parameter Helmert shift (the +towgs parameter)

I have defined a transformation from WGS84 to ETRS89 Poland CS92 (Epoch 2019.5) in CS2CS. I can do the basic map projection, but the datum shift is being ignored. It doesn't matter what numbers I type ...
  • 111
5 votes
1 answer
59 views

What is a source of Proj projections definitions?

I'm trying to create my own implementation of most popular projections (for reasons unrelated with GIS). But I want my implementation to be as compatible with Proj as possible (because it is an ...
2 votes
0 answers
40 views

ST_TRANSFORM NAD27 to WGS84

I've encountered an issue when transforming NAD27 coordinates to WGS84. The resulting WGS84 coordinates are off by quite a bit (10-12 feet in most cases) and it lead me to digging into why that is so. ...
2 votes
1 answer
47 views

gdalwarp to change my raster extents

I'm trying the gdalwrap command to set center and change the extent to -180 to 180. But for some reason, the bounding value is set but the pixel values are not transferred. This solution not working: ...
0 votes
0 answers
71 views

Get geoid height at a particular lat/lon (Python/Proj)

I am trying to calculate the geoid height at a particular lat/lon. Below is the sample code. I see that z3 = 10 and I was expecting a different number than 10 (input height was 10). What would need to ...
  • 103
0 votes
1 answer
33 views

Convert lat/lon to data in Georgia state plane

I am trying to convert some data in x,y,z (ECEF) to Georgia state plane. First I convert x,y,z to lat,lon,alt and then try to convert it into Georgia state plane. When I look at the data, the numbers ...
  • 103
0 votes
0 answers
52 views

Set Mapbox projection to be compatible with EPSG:27700

I have 5m rasters that I am adding as layers to mapbox. The rasters themselves are projected in EPSG:27700 (OSGB36), whereas mapbox is projected in WGS84. Because of this I am getting certain ...
1 vote
0 answers
21 views

PROJ pipeline not working in python

I am using gdalwarp in the commandline as so gdalwarp -ct "+proj=pipeline +step +inv +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +step +proj=hgridshift +...
0 votes
0 answers
49 views

Defining exact projection in GDAL command

I am using gdalwarp to convert a tiff from one projection to another. I want to specify the t_srs as the same transformation/projection type I've been using in QGIS which reads as below +proj=pipeline ...
0 votes
0 answers
103 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-...
1 vote
1 answer
24 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 ...
  • 431
0 votes
2 answers
52 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> ...
0 votes
0 answers
79 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 ...
  • 1
0 votes
0 answers
47 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 ...
1 vote
0 answers
78 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/...
  • 11
0 votes
1 answer
63 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 ...
2 votes
1 answer
80 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=...
  • 21
0 votes
1 answer
34 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, ...
1 vote
1 answer
350 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 +...
  • 143
2 votes
1 answer
115 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 +...
0 votes
1 answer
44 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 ...
  • 163
2 votes
0 answers
82 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=&...
  • 277
1 vote
0 answers
104 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 ...
  • 2,154
0 votes
0 answers
78 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 ...
  • 157
1 vote
0 answers
63 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 ...
  • 7,527
3 votes
2 answers
203 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 ...
0 votes
1 answer
60 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 ...
1 vote
0 answers
44 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 ...
2 votes
0 answers
89 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 ...
3 votes
1 answer
605 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 ...
  • 451
1 vote
1 answer
185 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 ...
  • 21
1 vote
0 answers
50 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 ...
3 votes
1 answer
114 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 ...
6 votes
0 answers
121 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, ...
  • 353
0 votes
1 answer
67 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#...
  • 81
0 votes
1 answer
45 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 (source: https://...
  • 1,600
4 votes
2 answers
166 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 ...
  • 2,154
2 votes
1 answer
68 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-...
  • 31
0 votes
0 answers
182 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 ...
0 votes
0 answers
171 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 ...
  • 717
1 vote
2 answers
113 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:\...
  • 676
2 votes
1 answer
135 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 +...
  • 1,169
1 vote
0 answers
283 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 ...
  • 11

1
2 3 4 5
11