Python interface to PROJ.4 library. Performs cartographic transformations and geodetic computations.
Questions tagged [pyproj]
184 questions
0
votes
1 answer
Cartopy: Albers Equal Area projection not working
I am trying to plot a topographic raster in Cartopy. I have downloaded some sample GeoTIFF data from this database: https://zenodo.org/record/3940482. I then import the data and metadata using the Python GDAL library:
from osgeo import gdal …

Fab von Bellingshausen
- 1,387
- 1
- 8
- 16
0
votes
0 answers
calculate lat/long bounding box using Pyproj, Shapely, GeoPandas
I need to calculate the min/max latitude/longitude co-ordinates within a specified distance (in km) from a given latitude/longitude using Python.
From my googling, it seems like the best packages out there are GeoPandas, Shapely and Pyproj but I'm…

user3014653
- 735
- 2
- 8
- 13
0
votes
1 answer
Geopands to_crs() dropping z values
I've got a csv of metadata (metadata_df) for some images I'm trying to process. I extract the geometry in the metadata to make a Point geometry object, then create a geodataframe using this geometry & the original metadata.
df =…

E. Case
- 67
- 1
- 2
- 11
0
votes
0 answers
How to create a uniform continuous grid (i.e., a geopandas.GeoSeries) composed solely of uniform shapely.geometry.BBox in parallel?
I have a script that generates a continuous regular grid (i.e., a geopandas.GeoSeries) that is composed solely of shapely.geometry.BBox instances. This script runs correctly in a single-core. Given the fact that this kind of operation can take a…

Philipe Riskalla Leal
- 954
- 1
- 10
- 28
0
votes
1 answer
How to pick points inside grids?
I have one grid that consists of two boundary arrays: lon_bnds and lat_bnds. The goal is to pick the points inside the grid.
Here's an example:
import numpy as np
import matplotlib.pyplot as plt
import pandas as pd
lon_bnds = np.array([[-77.9645 …

zxdawn
- 825
- 1
- 9
- 19
0
votes
0 answers
How can I obtain pixel value from map co-ordinate [longitude, latitude] for image of band count 1?
so I am using Landsat-8 image of band-4. I want to obtain the n-array index for given longitude and latitude values. I have tried below…
0
votes
1 answer
Geopandas: dataframe to geodataframe with different espg code
I have a dataframe (df2): wherein x,y are specified in rd new epsg:28992 coordinates.
x y z batch_nr batch_description
0 117298.377 560406.392 0.612 5800 PRF Grasland (l)
1 117297.803 560411.756 1.015…

L Blauw
- 1
- 1
0
votes
1 answer
pyproj conversion UTM to lat/long out by 3 deg
I am reading UTM point data from a shape file. The geopandas CRS string associated with the shape file is:
PROJCS["WGS 84 / Falk",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS…

PetGriffin
- 495
- 1
- 4
- 13
0
votes
1 answer
Pyproj DLL load failed while importing _network
I'm trying to import geopandas in Python. The pyproj module fails to import, with the following error message:
File ~\.conda\envs\lizard\lib\site-packages\geopandas\_compat.py:9 in
import pyproj
File…

lizard
- 123
- 2
- 7
0
votes
1 answer
Using pyproj to transform shapely data is giving me strange results
I've got a polygon which looks something like this in WKT:
POLYGON ((-2.5079473598836624 51.34385834919997, -2.5081726654409133 51.34353499032948, -2.507909808957454 51.343441165566986, -2.507679138982173 51.34359530614682, -2.5079473598836624…

Tom
- 7,269
- 1
- 42
- 69
0
votes
1 answer
PyProj issue in working with OSMnx package
I have recently installed the OSMnx package using Anaconda environment (python 3.10.2) based on the installation instruction provided here. Once I try to import the OSMnx package into Python environment, it shows the following…
0
votes
1 answer
Anaconda-Error: cannot import name Transformer from pyproj
I try to install Transformer from pyproj and get the message "cannot import name".
I installed pyproy using Anaconda Navigator. There it tells me I have version 1.9.6. (newest update). I found a thread about same issue but it didn't help me (cannot…

Fabian Vogt
- 27
- 5
0
votes
1 answer
polygon porjection to another crs and backwards to original crs results in different polygon
I'm trying to place a 50km buffer on each country in the world. For this I transform the country shapes from epsg:4326 to Lambert azimuthal equal-area(laea), apply a buffer to it and then transform it back to espg:4326. This works well for the most…

seb2704
- 390
- 1
- 5
- 17
0
votes
0 answers
Error reading spatial files in Python "Invalid projection: EPSG:4326"
Im having problems reading any kind of spatial file in Python (such as GeoJson or Shapefile) the error that I get is the following:
CRSError: Invalid projection: EPSG:4326: (Internal Proj Error: proj_create: SQLite error on SELECT name, type,…

Orlando G
- 307
- 3
- 11
0
votes
0 answers
Calculation of the angle of a point from north using pyproj
Is there a way to get of angle of a point data with respect to the North in python? So far I have tried to calculate the forward azimuth using pyproj.
import pyproj
from pyproj import CRS
# lat long in WGS84
lat, long = (10.08079707822896,…

windcheater
- 25
- 1
- 4