Questions tagged [pyproj]
A Python interface to the PROJ.4 library to perform cartographic transformations and geodetic computations.
0
votes
1answer
43 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
0answers
26 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
0answers
29 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
1answer
50 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
0answers
115 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
1answer
160 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
322 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
43 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
542 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
496 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
57 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
815 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
305 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
884 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
372 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
2k 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
314 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
190 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
591 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
141 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
198 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
360 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
128 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
0answers
106 views
How can I transform 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
1k 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
305 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
427 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
392 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
98 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
8k 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
1answer
132 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
0answers
402 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
269 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
572 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
7k 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/...
1
vote
0answers
709 views
Projecting data from 'EPSG:4326' to Lambert Azimuthal Equal Area with PyProj [closed]
I have a list of data points given in lon/lat coordinates, which I want to convert to cartesian ones using pyproj. My dataset refers to Italy, so I was considering using either the Lambert Azimuthal ...
1
vote
2answers
2k views
Latitude Longtitude to X, Y, Z coordinate system using python
does anyone know how to go about converting lat, long, alt to x, y, z where you have an origin for x/y coordinate (specific place in London)? In other words, I have a lat long and alt data and would ...
3
votes
1answer
440 views
Incorrect transformation when using pyproj
I am using pyproj to work with shp files and am relatively new to GIS. I am tryign to visualize the shp data, but am having a problem as I am incorrectly transforming the coordionates.
Here is an ...
1
vote
1answer
224 views
What does the “zone” stand for in proj?
I am using proj, in particular pyproj, to convert from UTM coordinates to the (x,y) coordinates used, for example, by Google Maps. There is a parameter called zone which I could not understand how I ...
1
vote
1answer
412 views
Blank crs attributes in fiona: Can't get proj4 parameters from OpenFileGeodatabase
I’m having trouble obtaining crs values for geodatabase feature classes in fiona (version 1.7.5). Specifically, the crs attribute returns empty for many of the geodatabases I try to read. This is not ...
2
votes
1answer
1k views
Why does Pyproj fail converting lon lat alt to ECEF?
I'm trying to convert a WGS84 lon, lat and alt to ECEF and Pyproj is giving me really weird results. My code is simply:
lon, lat, alt = 49.74761271, -113.2179781, 0
ecef = pyproj.Proj(proj='geocent', ...
2
votes
1answer
672 views
Are there EPSG codes for solving China GPS shift problem?
I am currently using PyProj for a project, and would like to transform from WGS-84 (or other systems) to the China-specific datums for writing into kml files. I have seen other packages specifically ...
1
vote
2answers
5k views
Converting WGS84 to ECEF in Python?
I am trying to convert some coordinates geodetic coordinates (WGS84) to ECEF. I've tried my own implementation in Python and using pyproj and am getting difference results. Could someone please help ...
6
votes
2answers
4k views
Converting EPSG:2284 to EPSG:4326 with pyproj
I have some points in EPSG:2284 and I'm trying to convert them to EPSG:4326 using pyproj, however I get weird results. When I go to the spatial reference website
I see that there is an endpoint that ...
2
votes
1answer
257 views
Difficulty converting State Plane to WGS 84
I have some parcel coordinates that I cannot convert to WGS 84 (cf. http://www.bcad.org/mapSearch/?cid=1&p=448116)
What's the coordinate system of these parcel coordinates, you ask? Great ...
1
vote
0answers
210 views
Transforming azimuthal to wgs projection using Python?
I want to plot a circle around a given lat/long with a certain radius on the Earth. The shapefile of the earth has WGS84 coordinate system.
First, I projected the circle to azimuthal equidistant map ...
9
votes
1answer
3k views
Convert shapefile from projected coordinates using Python
Newbie here struggling with GIS. I am trying to map out wards for the city of Milwuakee using shapefiles found on their county website county website. I am following the thread here with some success. ...
2
votes
1answer
461 views
What projection to use to calculate polygon areas of city to global size from lat/lon data?
I have a set of bounding boxes in WGS:84, and I would like to calculate their areas using Python.
The bounding boxes are of all sizes, from city sized to covering the entire Earth.
The bounding boxes ...
3
votes
1answer
466 views
Python pyproj converting from wgs84 to Israel Grid IG 05/12 EPSG:6991
I am trying to convert from lat/long to xy in Israel Grid 05/12:
links https://epsg.io/6991 (https://epsg.io/6991.gml)
I tried this :
p = pyproj.Proj("+init=ESPG:6991") using python package pyproj ...
