Questions tagged [pyproj]
A Python interface to the PROJ.4 library to perform cartographic transformations and geodetic computations.
124
questions
4
votes
1answer
65 views
Setting CRS to ESRI:102001 using GeoPandas
I've been passed some GeoJSON shape data that I've been told is in the CRS ESRI:102001, although that property has not been written into the GeoJSON file. I've been able to confirm that is the correct ...
1
vote
0answers
26 views
Manipulating Transformation parameters with PyQGIS
Is it possible to add a layer with a local coordinate system, transformed on the fly by a custom transformation entry from this local CRS -towgs84 and furthermore to manipulate the -towgs84 Parameters?...
1
vote
1answer
35 views
Updating Pyproj syntax results in other values
I want to know the length in meters of a LineString with coordinates in EPSG4326.
Based on this post, this is code that works for me:
from shapely.geometry import LineString
from shapely.ops import ...
0
votes
0answers
21 views
Convert heights from DHDN Gauß Krüger Zone 4 to UTM 32 ETRS 89 in pyproj
I am using pyproj.
X/Y conversion is pretty straightforward through
transformer = Transformer.from_crs(31468, 25832,always_xy=True)
However, this does not transform the heights.
From
Vertical ...
3
votes
2answers
134 views
Using QGIS Python interpreter outside QGIS
I'm writing a simple Python Geoprocessing script which leverages two specific modules:
from osgeo import gdal
import pyproj
The script must be shared with colleagues in my organisation.To avoid ...
1
vote
2answers
44 views
Project large collection of shapely object to different CRS [closed]
What is the most efficient way of projecting several thousands of polygons and multipolygons to a different CRS in Python? I have found that looping over the polygons and reprojecting them one by one ...
6
votes
1answer
110 views
How can I convert the Gaussian coordinates( Gauß-Krüger) to the WGS84 coordinates
I have a list of WGS84 and Gauss coordinates, but how to convert another WGS84 coordinates.
I want to use PROJ for conversion.
What is PROJ string like?
| Name | x | y | H | ...
1
vote
1answer
33 views
Creating custom proj4 string from components
I'm trying to work with NOAA's NDFD weather data, which are released as rasters in GRIB format. Specifically, I'm currently trying to get the pixels that intersect an arbitrary geometry.
Information ...
1
vote
1answer
103 views
MapProxy Installation problem with PyProj
I'm testing MapProxy (1.12.0) on a Windows Machine and followed each step from the "Installation on Windows" docs.
I have a running virtualenv, used pip to install MapProxy and PyProj.
The command ...
1
vote
1answer
30 views
Conversion of NAM Lambert Conformal coordinate system to NAD 1983 UTM Zone 11N
I am attempting to convert the NAM wind data coordinate system using data from the thredds server:
Catalog = TDSCatalog('http://thredds.ucar.edu/thredds/catalog/grib/NCEP/GFS/'
...
0
votes
1answer
42 views
Get type of projection from pyproj CRS object
I'm used to using PROJ strings/dictionaries in python where I have access to a proj_dict['proj'] element. I have some calculations that I can only do with the PROJ.4 geos projection and otherwise ...
2
votes
1answer
41 views
Gravity related height from ellipsoidal height using pyproj
I'm using the pyproj library (http://pyproj4.github.io/pyproj/stable/index.html) to get translate ellipsoidal heights into gravity related height.
Code:
from pyproj import CRS
from pyproj import ...
2
votes
1answer
89 views
How to create an hexagonal grid of regular hexagons of definite area anywhere on the globe, using shapely?
I have a geojson file describing a polygon. I want to create an hexagonal grid on top of this polygon, with regular hexagons of area 90000 square meters.
Right now, I can either guarantee the area, ...
0
votes
1answer
131 views
Vertical Datum transformation using Pyproj
I have been using Pyproj for some time now doing only simple projection and transformation. I am totally new to vertical datum and differences in GEOID among Coordinate systems.
I have been trying to ...
1
vote
1answer
42 views
Pyproj projection 0.5 meter differnce from official results (IGD05/12,IG05/12) Israel Grid
I'm trying to use those parameters and I receive a difference of ~0.5m
ALON and BSHM APN Points:
coords = [
(dms2rad([31,42,28.5694]), dms2rad([34,36,23.92196]), 55.671),
(dms2rad([32,46,44....
0
votes
1answer
910 views
Why do I get error “Runtime: b'no arguments in initialization list'” when running to_crs?
I have a geopandas dataframe that has POINT geometry in NAD83. I want to use to_crs to change them to espg: 4326. When I first bring in the feature class as a gpd dataframe, it shows no crs, so I use ...
2
votes
1answer
92 views
Transforming geometry points coordinates to WGS84 in geodataframe
I have shapefile that has been brought into my script as a geodataframe. Initially, it shows no crs although the source originally had a datum = NAD83.
I want the coordinates for the geometry (the ...
0
votes
1answer
103 views
`pyproj` initialization error in PyCharm but not in cmd?
I am currently struggling to get pyproj to initialize when running in PyCharm. Oddly enough, I am able to get it to work in cmd.
The version of pyproj I am using is 1.9.6, and the version of proj4 I ...
1
vote
1answer
32 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 ...
2
votes
2answers
454 views
Converting large data with lat and long into X and Y?
I have 9888562 records in dataframe and I would like to convert my lat, long to UTM x,y. according to my code, I have used pyproj package but because my data are too much it takes a long time and ...
0
votes
1answer
122 views
Converting XY coordinates in the Iberian peninsulate to lat, lon?
I have been looking for a way to convert XY coordinates (from a coordinate system based on ETRS89) to latitude and longitude coordinates.
So far I have this function:
from pyproj import Proj, ...
1
vote
1answer
37 views
Is height conversion between ellipsoids necessary?
I have measurement data that uses X/Y in different zones of JGD2011 and an ellipsoid height above GRS1980 (but it also includes a geoid height as well). My task is to transform this data to "WGS84 LLA"...
2
votes
1answer
74 views
Polygon area from Google Earth so different to QGIS & Shapely/PyProj?
I have a set of points defined by their longitudes and latitudes in WGS84 that define a polygon that I need to calculate the area of.
I originally tried to approach this using pyproj and Shapely as ...
1
vote
0answers
104 views
Is the Albers Equal Area projection only useful for calculating the area over small differences in longitude?
I've been trying to calculate the area of a map, given its WGS 84 coordinates, with the following code:
import pyproj
import shapely
import shapely.ops as ops
from shapely.geometry.polygon import ...
3
votes
2answers
154 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 ...
1
vote
1answer
134 views
Generating custom 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. In aerospace ...
2
votes
1answer
63 views
Why is a projection conversion off by factor of SQRT(2)?
I am attempting to draw a 3 km ring around downtown Minneapolis. Let downtown be defined as Government Plaza (-93.2661, 44.9765). I am using this ring to calculate census tracts within a distance of a ...
5
votes
1answer
182 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
210 views
Getting EPSG code in pyproj with version==2.1.3?
I'm using pyproj in python and I want to get epsg code based on an existing Proj class. An example in version=1.9.6:
from pyproj import Proj
CGCS2000metre = Proj(init='EPSG:4544')
print(CGCS2000metre....
2
votes
2answers
262 views
pyproj EPSG:4326 converts degrees to radians
Case
My input is a point given as x, y and an SRS name which specifies how to interpret the numbers. My goal is to convert these into lon, lat in degrees.
I convert this by using this piece of code:...
1
vote
1answer
66 views
Custom projection in GeoServer
I am trying to create a custom projection in GeoServer, but I see nothing (from the logs, neither in the SRS_list) coming out.
The projection I'm trying to create is a Mercator for Mercury.
I start ...
0
votes
1answer
263 views
Change CRS of large GeoTIFF and extract windows from large raster in Python
This is the first time I am working with GeoTIFF formats and I want to do this in Python. I am trying to use rasterio and pyproj for the following tasks.
I have already made some first steps based on ...
0
votes
1answer
48 views
Retrieving lat-lon coordinate from projdef?
I am trying to plot an hdf5 file where geographical domain is defined like :
group
LL_lat : 23.55908203125102
LL_lon : -138.4263916015652
LR_lat : 16....
0
votes
1answer
85 views
Get central meridian, PyProj?
I am using PyProj and I am trying to find a way how to get Central Meridian and Latitude of the center of the zone for my State Plane coordinate system.
1
vote
2answers
162 views
PyProj giving RuntimeError: b'argument not numerical or out of range'?
I have this string:
string = '+proj=lcc +a=6371200.0 +f=0.0 +pm=0.0 +x_0=5632642.22547 +y_0=4612545.65137 +lon_0=-107.0 +lat_1=50.0 +lat_2=50.0 +lat_0=50.0 +units=m +axis=enu +no_defs'
Then I do:
...
2
votes
1answer
421 views
Transform coordinate from rotated lat/lon to normal lat/lon WGS84
I have access 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 ...
0
votes
2answers
495 views
Projection conversion with Python/Pyproj: Canada Albers Equal Area to WGS
How do I get the geographical coordinates (lat and lon) from x,y values in Canada Albers Equal Area projection?
The range of x and y values as follows:
min(x) = -5171461.906673405
max(x) = 3372679....
2
votes
1answer
925 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
58 views
Trying to implement Albers projection in Python, what is going wrong?
I have a platform where I can't import Proj4. I only need to convert points to Albers projection, so it seemed simple, but my conversion has some errors:
Projection definition: +proj=aea +lat_1=52 +...
2
votes
1answer
1k views
Transform coordinates from custom projection to lat long PyProj
I am trying to convert a pair of coordinates specified in US survey feet on a custom transverse mercator projection to lat/long using PyProj. I must be doing something wrong because I cannot get the ...
2
votes
1answer
2k views
Converting lat,lon to x,y coordinates with pyproj
I am trying to convert lat,lon to x,y coordinates in France, using the pyproj library and what I am doing is obviously wrong.
Here is some of my code:
import pyproj
import math
PROJ='+proj=utm +...
0
votes
1answer
88 views
Why is pyproj shifted by 33 feet from QGIS?
I have written a script to take WGS84 coordinates and convert them into a projected system for importing into a CAD program. I thought I had everything working when I decided to use QGIS to confirm ...
1
vote
1answer
2k views
Easting and northing conversion into lat and long gives wrong information python
This is my easting and northing
315846.356,5808849.713
p2 = Proj(proj="utm",zone='55H')
x,y=315846.356 ,5808849.713
lat,long=p2(x,y,inverse=True)
I got Longitude 144.29314340365048 and Lattitude 52....
1
vote
1answer
541 views
Reprojecting coordinate system with Python and lat/lon is only slightly off [duplicate]
I am using Python to reproject a shapefile (called "SCL" in my project) into a new coordinate system that will allow me to plot it on the same map as another shapefile (called "PUMA" in my project). ...
0
votes
3answers
2k views
Translate Cartesian coordinates to WGS84, preferably in pyproj?
I have some Cartesian coordinates:
x = 3862525.511
y = -109738.522
I need to translate them into WGS84 coordinates, preferably using pyproj. Is this possible, or have I misunderstood the nature of ...
0
votes
2answers
625 views
Rinex Header GPS Position Coordinate System
I am extracting data from rinex observation files and would like to convert the approximate position to latitude and longitude using Proj.4. The rinex standard states that these coordinates are WGS84, ...
2
votes
1answer
4k views
Converting EPSG:2263 to WGS84 using Python Pyproj
I'm trying to convert Projected bounds to WGS84 using Pyproj but getting strange results.
The Projected Bounds are : 909126.0155, 110626.2880, 1610215.3590, 424498.0529
The corresponding WGS84 ...
0
votes
1answer
565 views
Convert NAM Grid (Lambert Conformal Conic) to MERRA2 Grid (lat/lon) using pyproj
I am attempting to convert NAM grid data into a simple lat/lon grid that my current MERRA2 data is on (NAM Grid Info). The NAM data is currently in lat/lon but on a Lambert Conformal Projection. I've ...
0
votes
1answer
312 views
PyProj: Transform from LLA to ORTHO with WGS84 gives sphere instead of ellipsoid
I want to convert my coodrinates from lon, lat to x, y of orthographic projection. I tried to use PyProj to do it, but it uses sphere, not ellipsoid. Here it is:
lonlat_proj = pyproj.Proj(proj="...
7
votes
1answer
1k views
Different origins of UTM projection
I have latitudes and longitudes from a GPS unit which uses WGS84, and would like to project those to some kind of plane.
I'm totally new to this, but after a lot of reading it sounds like UTM is a ...
