Questions tagged [pyproj]

Python interface to PROJ.4 library. Performs cartographic transformations and geodetic computations.

See also:

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…
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…
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