A message from our CEO about the future of Stack Overflow and Stack Exchange. Read now.

Questions tagged [pyproj]

A Python interface to the PROJ.4 library to perform cartographic transformations and geodetic computations.

Filter by
Sorted by
Tagged with
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 ...