Questions tagged [proj]
A library for performing conversions between cartographic projections.
366
questions
1
vote
1answer
12 views
Incorrect false northing in pyproj?
I am using version Rel. 5.2.0, September 15th, 2018 on OSX but get the same answers on other platforms with other versions.
prompt> echo 3 0 | proj +proj=utm +zone=1 +ellips=wgs84
500000.00 ...
3
votes
1answer
35 views
Geos projection not centered on the equator?
I'm working with QGIS and satellite view projections, in order to achieve something similar to the maps created by Richard Edes Harrison for Fortune magazine in the '40s (reference). I have a question ...
2
votes
2answers
40 views
Generating a custom flat projection in pyproj
I have little experience with pyproj, but I want to achieve the following task:
I have a trajectory of a robot, given in lat/long coordinates which I want to convert to x/y Coordinates. And I want to ...
2
votes
0answers
14 views
How to run OSGeo4W shell for molodensky transform?
I'm a beginner with Proj, so I just downloaded proj, and I'm using the OSGeo4W Shell and when I try to run the example
proj=molodensky a=6378160 rf=298.25 da=-23 df=-8.120449e-8 dx=-134 dy=-48 dz=...
0
votes
0answers
20 views
Custom metric projection definition in QGIS with reference point and azimuth [duplicate]
I need to create a custom CRS.
I have a reference point in UTM NAD83 18N, an azimuth (30) and the extent. The projection must be metric, I know I have to use PROJ4 but the language is not very ...
5
votes
3answers
163 views
Making a GeoTIFF using Google Earth Data
I'm trying to make a GeoTIFF using geotifcp, which I'll then upload to GeoServer.
I've got an image overlay (part of a .kml) in Google Earth, so I tried to use its bounding box to create a world file,...
2
votes
1answer
50 views
CRS validity and PROJ versions
I'm not sure where the right place to raise this is, so:
The PROJ strings for EPSG codes 4151 and 4283 are identical:
+proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs
which is a hassle ...
2
votes
1answer
102 views
How to use a specific transformation for a coordinate projection with Python
I try to use Python to project single coordinates from Gauß-Krueger-coordinate system (EPSG:31467) to ETRS89-UTM (EPSG:25832).
I'm using the way Antonio Falciano showed here.
So my script looks like ...
5
votes
1answer
92 views
Explaining pyproj to_epsg min_confidence parameter?
Following on from this question How can I get epsg code in pyproj with version==2.1.3, I discovered that sometimes when intialising a pyproj Proj instance with a string ESPG code eventually results in ...
2
votes
1answer
145 views
How to create a grid of pixel coordinates from the corners of a MODIS tile?
I have an HDF4 file whose StructMetadata.0 contains the following attributes:
UpperLeftPointMtrs = (-20015109.354000,1111950.519667)
LowerRightMtrs = (-18903158.834333,0.000000)
These are X and ...
-3
votes
1answer
35 views
Proj4 give always a false answer [closed]
I'm trying to use Proj4 (dll) to convert from Lambert93 to WGS84 but it always give a wrong answer:
// EPSG:2154 : Lambert 93 definition: https://spatialreference.org/ref/epsg/2154/proj4/
const char*...
1
vote
0answers
39 views
Applying 3D-Helmert Transformation Parameters in QGIS?
I have received a set of local transformation parameters from my surveyor to go from WGS84/UTM36S to ARC1950/UTM36S (See picture below)
The way QGIS does the transformation is not the best. I am ...
0
votes
0answers
15 views
NDFD Grib2 geotransform
I am trying to transform the projection for my data in grib2 format to long/lat, but I'm not sure I'm using the correct EPSG. The data says it uses a "Lambert_Conformal" projection with a WGS84 ...
2
votes
0answers
57 views
WGS84 to local coordinate system [closed]
I have a point cloud and would like to move this to the local coordinate system for example from Revit facing toward the first quadrant. We have a file with 3 ground reference points and on-site ...
0
votes
0answers
18 views
proj: Oblique Mercator specifying two points on central line
I'm trying to use Oblique Mercator in QGIS 3.4, which uses proj to calculate coordinates conversion. I specify two points in the central line, rather than specifying the center and alpha.
Oblique ...
3
votes
1answer
40 views
Creating custom CRS from existing CRS EPSG:22780?
I am currently trying to figure out how to create a custom CRS in QGIS based on Deir ez Zor / Levant Stereographic (EPSG:22780). For this I am working with a map which shows both the EPSG:22780 ...
1
vote
1answer
105 views
Transforming coordinates using R - SWEREF to WGS84
I have some coordinated in SWEREF 99 TM which I cannot convert into WGS84 using the proj4 package in R.
temp <- data.frame (x = c(598223, 598812, 598824, 598232, 597614, 597629),
y = c(...
0
votes
0answers
30 views
Are Azimuthal Equidistant and ESRI:54032 the same thing?
I am looking at
(A): https://proj4.org/operations/projections/aeqd.html#cmdoption-arg-ellps
and
(B): http://epsg.io/54032 Description: +proj=aeqd +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +...
0
votes
0answers
48 views
Calculating tile coordinates based on location and crs for ArcGIS maps?
I am working on a solution to download map tiles based on given location:(latitude, longitude) and zoom level
Followed https://wiki.openstreetmap.org/wiki/Slippy_map_tilenames
The key concept is ...
6
votes
3answers
3k views
Build gdal with proj version 6
Description of the issue:
Building gdal from sources after having build proj in its latest version is failing at the configure step because it's not finding proj 6.
Details to reproduce:
Work ...
0
votes
1answer
82 views
Get proj4 parameters from EPSG in C++
I have the EPSG number of the projection and I want to use proj4 to get the projection descriptive string (like +proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 ...
0
votes
1answer
42 views
MapProxy doesn't start with EPSG:28992 while it starts and works fine with EPSG:4326
We are moving from Debian to CentOS and in that process we discover that MapProxy installation is turning out tougher in CentOS than it was in Debian.
One of the dependencies 'libproj' is not ...
0
votes
0answers
39 views
Creating a custom “flat” projection definition in proj4
I want to create a simple, custom projection in proj4. I need it to be completly flat ( rectangular coordinates ). Starts in point (0,0), meters as unit, everything else - as simple as possible. What ...
2
votes
1answer
32 views
Managing custom projections
I am developing a mining sampling database. There will be many project specific, user created projections. Typically they will be Oblique Mercator defining a local coordinate projection system.
...
1
vote
0answers
33 views
Create a Non-North/Rotated-Vertically User Coordinate System?
I need to plot 3d PostGIS data into vertical planes for symbolization and drawing vertical cross-section interpretations for mineral exploration models. I have written a short plpgsql function that ...
2
votes
1answer
102 views
Lambert conformal conic alternative projection in PROJ4
What's the difference between Lambert conformal conic projection and Lambert conformal conic alternative projection in PROJ4?
0
votes
0answers
68 views
UTM to latlong in R - with proj4
I have received a shapefile with measurement points to be analyzed, but no matter how I try, I can't find how to convert to lat lon. I suppose it would be good to know the projection but when I look ...
0
votes
0answers
19 views
Coordinate system transformation deviation [duplicate]
I'm having the same problem as described in this Question from September 2018.
Has it ever been answered satisfactorily?
I'm attempting to transform from EPSG 3785 to 4326. The X coordinate output is ...
2
votes
1answer
258 views
Transform coordinate from rotated lat/lon to normal lat/lon WGS84
I have acces to an OPeNDAP service that allows me to peek into a humongous ocean wave dataset. It takes latitude and longitude as inputs (in fact even a range of the two) - but the coordinate system ...
6
votes
1answer
55 views
Are proj4 strings case sensitive? [closed]
Is there any reason that case of characters in Proj4 strings would affect their interpretation? Before parsing a Proj4 string, can I lower case all the characters first?
1
vote
0answers
94 views
Understanding proj4, epsg and esri *.prj in depth? [closed]
Where can I find a good resource that teaches me
1) about the difference between proj4, EPSG and ESRI *.prj? (Why do several coordinate system definitions exist? When should I use which? How are they ...
1
vote
0answers
22 views
WFS how to consume with different coordinate system? [duplicate]
I have a WFS in EPSG:4167. My map view has EPSG:3857. How to make the WFS work on a map view with different coordinates?
I used this to define my coordinate system from wfs:
import { register } ...
1
vote
2answers
119 views
Creating generic coordinate system in QGIS
So I'm new to QGIS (3.4) and haven't touched any mapping software in 15 or more years. I'm slowly finding my way around but I'm currently struggling to create a generic coordinate system for a map ...
0
votes
1answer
159 views
Is MODIS sinusoidal tile boundary in Eathdata website real?
I was using MOIDS data MCD12Q1 (MODIS/Terra+Aqua Land Cover Type Yearly L3 Global 500m SIN Grid V00). I try to use PyProj in Python to transform the Sinusoidal projection to lat/lon. It actually ...
5
votes
2answers
256 views
Convert WKT coordinate reference system to PROJ4
I'm currently working on a better support of LAS 1.4 file format in the rlas R package. In LAS 1.4 format the CRS of the point cloud can be stored as a WKT string. For example this is what I read from ...
2
votes
1answer
102 views
RSAGA error using rsaga.pisr2
I just downloaded SAGA and the RSAGA package for R. I would like to use the rsaga.pisr2 function to compute solar radiation for my region of interest. I am currently getting errors related to proj4 (...
0
votes
0answers
38 views
Gdalwarp producing odd results around 0 degree longitude
I am warping ECMWF grib data (full globe, 0.25 degree resolution) to match the SEVIRI geostationary satellite projection. I use this command:
gdalwarp -overwrite -ot Float32 -t_srs "+proj=geos +a=...
0
votes
0answers
48 views
How to import a data with unknown coordinate system to a data with known coordinate system [duplicate]
I have recieved a Group of shape files with unknown coordinate System. The coordinates are like: 536605 , 5658923 and I can not find the corresponding projection/coordinate System. It looks like that ...
0
votes
1answer
67 views
What parameters do azimuthal and transverse cylindrical projection in secant case have as proj4 string
conical equal-area projection, secant case has e.g.
+proj=aea +lat_1=... +lat_2=... +lon_0=... +lat_0=... +x_0=... +y_0=... etc...
I would ask if there is such thing, for transverse cylindrical ...
2
votes
0answers
33 views
Why is proj version dated in the future?
I compiled proj (also known as proj.4) from sources on Ubuntu (https://github.com/OSGeo/proj.4.git).
I now wonder why is proj version dated in the future?
E.g.:
$ proj
Rel. 6.0.0, March 1st,...
3
votes
0answers
48 views
Accurate area calculation for countries worldwide individually [closed]
I have polygons of countries with its administrative divisions, I would calculate the areas of each division in a country accurately.
I got the bounding boxes for each country, then I divided its ...
2
votes
2answers
116 views
Hammer projection in GDAL
I am working on Ubuntu 18.04 with PROJ 4.9.3 and GDAL 2.2.3 (installed from the ubuntu-gis PPA).
In a previous question on Hammer's projection in QGIS it was concluded that with older PROJ versions ...
1
vote
0answers
71 views
Create custom projection for Shapefile, given specific target coordinates [duplicate]
I have a somewhat odd situtaion. A client has provided me with a Shapefile that has no projection information, along with an image of the map that it corresponds to. I have run some coordinates ...
0
votes
0answers
46 views
Fishnet programmed in python does not display correctly in ArcMap but displayed correctly in QGIS
I programmed a fishnet in python to use as a reference map for georeferencing. I am facing the following issue:
The fishnet is not displayed correctly in ArcMap 10.5.1. The fishnet has parameters lat-...
3
votes
2answers
270 views
Handling raster data in a “rotated pole projection” with R?
My question is about using R to transform raster data from a rotated pole projection to a regular latlon grid - and then plotting it.
I am working with data from a numerical weather model, which is ...
1
vote
1answer
581 views
Converting NAD83 (epsg:4269) to WGS84 (epsg:4326) using pyproj
I'm trying to convert US census data (NAD83) to WGS84.
I know the difference between the two is quite minor, but trying to apply the transformation using pyproj does not do anything. e.g.:
from ...
2
votes
1answer
194 views
Can't perform Helmert transformation in Proj.4
I'm trying to convert a point from one SRS to another using proj4 strings and GDAL. My origin SRS uses Helmert transformation to convert point to EPSG:2229 so I tried to use pipeline operator (https://...
0
votes
1answer
124 views
Moving MS4W from Apache to IIS
I'm trying to switch to IIS on my MS4W MapServer install. It works fine in Apache but on IIS I get:
msProcessProjection(): Projection library error. proj error "no arguments in initialization list"...
0
votes
0answers
132 views
Different ways of NTv2-Grid integration
In QGIS there is an entry in the tbl_datum_transform for a NTv2-grid called BWTA2017.gsb ( https://github.com/qgis/QGIS/pull/6980 ).
So theoretically it would be enough just to activate the option "...
1
vote
0answers
32 views
Rotating North in Custom CRS using proj4?
I'm trying to create a CRS for a local mine grid that involves the translation and rotation of GDA94 / MGA zone 55. The origin is translated 2195129.073 in the +x direction and by 5489615.789 in the +...
