Questions tagged [proj]
A library for performing conversions between cartographic projections.
533
questions
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-...
- 41
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 ...
- 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> ...
- 1
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 ...
- 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 ...
- 1
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/...
- 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 ...
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=...
- 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, ...
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 +...
- 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 +...
- 160
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 ...
- 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=&...
- 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 ...
- 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 ...
- 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 ...
- 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 ...
- 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 ...
- 1
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 ...
- 161
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 ...
- 328
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 ...
- 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 ...
- 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 ...
- 11
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 ...
- 31
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, ...
- 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#...
- 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....
- 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 ...
- 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-...
- 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 ...
- 202
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 ...
- 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:\...
- 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 +...
- 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 ...
- 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 ...
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 ...
- 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 ...
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 ...
- 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(...
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 ...
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: ...
- 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 ...
- 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 ...
- 31
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 ...
- 121
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 ...
- 33
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 ...
- 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 ...
- 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' ...
- 141
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 ...
- 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 ...
- 11