Questions tagged [pyproj]
A Python interface to the PROJ.4 library to perform cartographic transformations and geodetic computations.
107
questions
0
votes
1answer
17 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
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 ...
1
vote
2answers
88 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
59 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
0answers
21 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
56 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
51 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 ...
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 ...
1
vote
1answer
48 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
61 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
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
105 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
122 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
44 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
174 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
46 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
63 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
118 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
254 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 ...
0
votes
2answers
347 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....
1
vote
1answer
576 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
56 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
885 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
995 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
77 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
1k 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
444 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
1k 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
490 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
3k 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
470 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
257 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
861 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 ...
0
votes
1answer
152 views
How to align Aerial images from MassGIS with OSM using python?
MassGIS Dataset: I'm using the MassGIS Data: USGS Color Ortho Imagery (2013/2014) aerial images of Boston. Basically, the dataset is divided into small patches of jpeg2000 images, where each of the ...
1
vote
1answer
292 views
pyproj does not properly convert EPSG:3760 to WGS 84 [duplicate]
I have coordinates stored in EPSG:3760. I can use https://epsg.io/transform to transform these properly to WGS 84 GPS coordinates, like so:
As seen in the photo, my coordinates (1728779.950564162, ...
2
votes
2answers
459 views
What is the difference between a datum and an ellipsoid? [duplicate]
For example I am converting from ECEF to lla using
def ecef_to_lla(x, y, z):
lla = pyproj.Proj(proj='latlong', ellps='WGS84', datum='WGS84')
ecef = pyproj.Proj(proj='geocent', ellps='WGS84', ...
1
vote
1answer
171 views
Transform Box from SR-ORG:6630 to WGS-84 in Python
I have an image file that I'm opening with gdal and using GetGeoTransform() to determine the minimum and maximum lats and lons:
Lon_Min=-2493045.0
Lon_Max=2342655.0
Lat_Max=3310005.0
Lat_Min=177285.0
...
2
votes
1answer
136 views
Transforming coordiantes to a topocentric system using PROJ.4?
I would like to transform coordinates (e.g. WGS84 or Bejing1954) to a topocentric system.
There seem to be several EPSG codes for this (5819, 5820, 5821, 15594, found here or here) but none of them ...
1
vote
2answers
2k views
geopandas can't import
I've installed geopandas on ubuntu 17.10 and can't import geopandas. I can't find this particular error anywhere.
from pyproj import _proj
ImportError: cannot import name '_proj'
I'm using a virtual ...
1
vote
0answers
379 views
Shapely LineString contains method not showing contained points
I have some GeoJSON files that I converted from some KML someone gave me to analyze.
I have used GPSBabel to convert them from KML to GeoJSON. The reason for doing so is that I am able to easily ...
1
vote
1answer
511 views
Converting coordinates from shapefile using python?
I am really new to working with this kind of data.
I have been trying to work with a shapefile that I wish to convert to geojson. The source has both shapefiles and geojson downloads available, but ...
2
votes
1answer
440 views
Bessel_1841 to Mercator in Python
Despite similar Q&As I can't seem to get my answer.
So I have this set of coordinate (124616.1535508, 486555.4631204) are mapped to the Bessel_1841 coordinate space. I'd like to transform these ...
1
vote
0answers
131 views
Cartesian / Meters issue while making ellipses
I am trying to create an ellipse to be stored in my django project. I'm provided a lat, lon, major/minor in Nautical Miles, and an orientation in degrees. Between this post and the docs I've been able ...
6
votes
3answers
11k views
Converting projected coordinates to lat/lon using pyproj fails
I have been trying to work on this issue on myself mainly referring to this post (Converting projected coordinates to lat/lon using Python?) but still didn't succeed.
I have some 2-D data function ...
1
vote
2answers
168 views
Something going wrong trying to transform longitude, latitude to NAD 83 using pyproj
I'd like to convert a latitude longitude pair to NAD 83. I've read NAD 83 is equivalent to epsg:4269 and the shapefile I am working with is a US Census TIGER shapefile, all of which are supposed to ...
1
vote
1answer
487 views
PyProj - Transform coordinate in “HK80 grid” to “WGS 84 / UTM zone 49N or 49S”
Am trying to use python pyproj to do the above conversion (convert/transform from Hong Kong GRID to Hong Kong UTM), but am getting different values from the expected result.
My expected conversion ...
1
vote
1answer
312 views
Python Egg Extraction Error When Trying to Use pyproj in Pyspark
I have pyproj installed in my conda environment on my local machine (MacBook Pro). Everything works properly. I am able to use pyproj within the pyspark installation on my local machine. However, I ...
2
votes
2answers
687 views
Buffer does not intersect or contain the point used for the center of the buffer when creating the buffer
I made the buffer below (see code). I then tried to see if the buffer contains or intersects the center of the buffer, but False is returned. Why?
from shapely.geometry import Point, Polygon
...
2
votes
3answers
2k views
Projection of all the polygons contained in a geojson at once
I would like to project all the features in a geojson from WGS84 (EPSG:4326) to Web Mercator (EPSG:3857) using python.
This geojson contains only polygons.
I found out the following shapely code ...
3
votes
1answer
8k views
Convert GPS coordinates to Web Mercator EPSG:3857 using python/pyproj
I've created a couple of functions using the math package from python 3.5 to transform GPS coordinates to x/y web mercator (pseudo mercator) as explained here: http://wiki.openstreetmap.org/wiki/...
