113 questions
0
votes
0
answers
31
views
Maps used as context using contextily doesn't match with basemap
Both sources use the the same WGS84 Web Mercator projection, but the maps doesn't match:
import contextily as cx
import matplotlib.pyplot as plt
from mpl_toolkits.basemap import Basemap
# Map source
...
- 159
1
vote
1
answer
34
views
problem to transform from epsg:3116 to lat, long (epsg:4326) in dart-flutter
I was trying to convert from epsg:3116(gcm_magna) to epsg:4326(lat and long) using the lib but the outputs were inaccurate. Can some one explain me or correct me?
var destProjection = Projection....
1
vote
0
answers
102
views
MapProxy: The requested tile is outside the bounding box of the tile map
Some context:
I have a Map Proxy server with the following configuration:
services:
demo:
tms:
use_grid_names: true
layers:
- name: my_layer
title: My example layer
sources: [...
- 53
0
votes
1
answer
39
views
Which EPSG does SUMO use?
I want to connect sumo with unreal engine use real world map.
According to this page (https://sumo.dlr.de/docs/Simulation/Distances.html), Sumo simulation network uses multiple coordination system:
x,...
- 1,112
0
votes
1
answer
54
views
Why can't I cut a buffer, both in R and QGIS?
My question is very straightforward.
I am working with a large coordinates dataset (animal GPS tracks). I would like to create a buffer of a 50km radius around the point where all my tracks start.
No ...
- 29
1
vote
1
answer
298
views
Wrong coordinate transformation with osr.CoordinateTransformation
I'm running a Python script (Python 3.9.16) in ubuntu 22.04.3 LTS. I'm doing it through Docker, althouh I've doubble-checked in Jupyter Notebook and I'm obtaining the same result.
I want to transform ...
0
votes
0
answers
116
views
Creating openlayers feature in the backend
I have a website where I render a map using openlayers, I save the features (polygons) on the map to the database so I can load them later by doing this:
const format = new GeoJSON({ featureProjection:...
- 31
1
vote
0
answers
363
views
EPSG:4326 to EPSG:27700 conversion
I am projecting some EPSG:4326 coordinates into EPSG:27700. The results are not what I am expecting, so I am wondering if there my code is wrong. Do I need to specify some additional datum in the ...
1
vote
1
answer
77
views
How can I choose to not change the image size when reprojecting with Python rioxarray?
I have an image in tif-format of size 18346x10218 in coordinate system EPSG:4326 that I want to convert to a new coordinate system, which I do with the following method with rioxarray in Python.
...
- 403
0
votes
1
answer
251
views
OpenLayers transform EPSG
I'm preparing a project that visualizes MULTIPOLYGON data on an OpenStreetMap. I'm testing the solution on a sample area defined as MULTIPOLYGON(((497700.486018679 520452.398175831,497716.241129864 ...
- 3
1
vote
0
answers
115
views
QGIS - Incorrect distance measurement using EPSG: 25832
I have an issue regarding distance and area measurements in QGIS (includes buffer etc.). My CRS is set to EPSG: 25832 (Germany). Measurements differ substantially from those of my colleagues (using ...
- 11
0
votes
1
answer
232
views
How can I use lon/lat in a netCDF with rlon/rlat dimensions?
I have a netCDF files with 5 coordinates but 3 dimensions :
xarray.Dataset {
dimensions:
rlat = 183 ;
rlon = 182 ;
time = 1 ;
coordinates:
float32 lat(rlat, rlon) ;
float32 lon(...
- 3
1
vote
0
answers
256
views
Transform coordinates with PostGIS based on Proj4 transformation text
I am working on a project where I have to transform coordinates expressed in LUREF (LUxemburg REference System/Frame) to WGS84 (Lat/long) by using the bundle PostGIS of Postgres database.
The Proj4 ...
- 13
0
votes
1
answer
424
views
Fundamental misunderstanding with crs and map projections in R
I think I have a fundamental misunderstanding with how using crs and map projections work in R. First let me show what my final goal is here. I would like to create something like the following:
...
- 111
0
votes
1
answer
154
views
geopandas to_crs epsg:3347
import geopandas as gpd
import matplotlib.pyplot as plt
world = gpd.read_file(gpd.datasets.get_path('naturalearth_lowres'))
world.to_crs('epsg:3035').plot()
plt.show()
changing to epsg:3347 gives ...
- 154
0
votes
0
answers
264
views
OpenLayers not showing Vector from WFS
I'm working on creating a web GIS app with OpenLayers. I'm having a difficulty with fetching data from WFS GeoServer.
The data is correctly downloaded, however it is not being drawn in any way. It's ...
- 21
4
votes
2
answers
4k
views
geopandas to_crs() return ProjError: x, y, z, and time must be same size
How can I convert this geodataframe with polygons in it from epsg:32748 to epsg:4326?
I have a geodataframe called poi_shp_utm_buffered which looks like this
| | ...
- 194
0
votes
1
answer
448
views
QGIS: 2 layers using different CRS
I am trying to draw a map of tree protection orders locally. I have a .shp file which uses CRS: EPSG 4326. I can add this to a project and see something that looks sensible, but with no context of ...
0
votes
0
answers
224
views
How to create Geopandas buffers for locations around the World
I have a location dataset (points long & lat) for events that happened around the World. I want to create 5km rings around these locations using Geopandas.buffer function. As explained in a ...
- 1
1
vote
2
answers
708
views
Transform laltitude and longitude in python pandas using pyproj
I have a dataframe as below:
df = pd.DataFrame(
{
'epsg': [4326, 4326, 4326, 4203, 7844],
'latitude': [-34.58, -22.78, -33.45, -33.60, -30.48],
'longitude': [122.31, 120.2, 118.55, ...
- 1,043
0
votes
1
answer
695
views
Proj4: How to convert coordinates from EPSG 4326 to EPSG 31370
Need to convert coordinates of Leaflet.js Map from EPSG 4326 to EPSG 31370.
What I did till now: Defining EPSG 31370
proj4.defs["EPSG:31370"] = "+proj=lcc +lat_1=51.16666723333333 +...
0
votes
1
answer
1k
views
proj_create_from_database when reading a raster with terra package
I am trying to read several rasters using the terra package; first I get a vector with the file names and then I use the rast() function to read the files.
Land_variables_files <- list.files("...
0
votes
1
answer
243
views
GeoPandas: Changing the CRS in a GeoDataFrame created through iteration
I'm having a bit of trouble changing the CRS in a GeoDataFrame created by iterating over a larger dataset. The epsg is set to the correct value (5070, an Albers Equal Area projection using meters) but ...
- 1
0
votes
0
answers
389
views
GDAL/PROJ fails to get SRS for gdalwarp
I'm having an ongoing issue with latest versions of GDAL (3.5.2, 3.6.0). I have an application that depends on the GDAL library for a number of functions, in this case merging multiple rasters and ...
0
votes
0
answers
265
views
Combining WMTS and WMS in OpenLayers for custom projection
I want to combine two layers for my map, ortophoto map and on top of that, polygons showing parcel boundaries. For ortophoto I have available WMTS and for polygons WMS.
Projection I need to use is ...
- 21
0
votes
1
answer
759
views
Raster and Shapefiles not lining up using Geopandas, Rasterio, and Contextily
I am trying to get a DEM raster to line up with a shapefile in Python, but it will not show up no matter what I do. This is for lab exercise, the entire rest of the exercise relies on these lining up, ...
0
votes
1
answer
271
views
datum conversion from NAD27 + NGVD29 height (ftUS) to NAD27 + NAVD88 height (ftUS) using python
I am using python pyproj to convert xyz values of Louisiana, USA region from NAD27 + NGVD29 height (ftUS) to NAD27 + NAVD88 height (ftUS). So, I am not changing the horizontal projection (because it'...
- 404
0
votes
0
answers
300
views
Offset of 170m when plotting a map with Cartopy add_wms in EPSG 2056 (LV95) from SwissTopo
Using Matplotlib and Cartopy, plotting a map from Swisstopo is shifted by about 170m. The official coordinate system of Switzerland is EPSG 2056 ("LV95"). I create an axes in matplotlib with ...
- 101
0
votes
0
answers
85
views
features at "wrong" position / flipped after EPSG:31287 -> EPSG:4326 transformation
I am currently trying to convert WFS2/GML3.2 data in EPSG:31287 to geojson using openlayers 6.15.1. For geojson I am trying to convert it to EPSG:4326/WGS84 using the following Typescript snippit:
...
- 315
-1
votes
1
answer
46
views
Is there any python method that is able to extract value of key between [ ] sign in a long string?
I am using geodjango ,and i would extract the value of "EPSG" from the long string returned from a code as below:
PROJCS["WGS_1984_UTM_Zone_42N",GEOGCS["WGS 84",DATUM[&...
-1
votes
1
answer
285
views
Follium map : my European locations are plotted in Africa
I would like to use Folium map to plot markers. My locations are in France. I have latitude and longitude information. So I create POINT geometry in order to implement them in Folium map.
df = pd....
user16313671
1
vote
2
answers
195
views
Does sgo R package give accurate UK latitude and longitude?
Newbie here, so please be gentle with me!
My council has released some data about tree protection orders (TPOs) in EPSG 27700 format. So an example would be 448098.4741, 271753.1604
I want to convert ...
0
votes
1
answer
2k
views
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 &...
- 67
0
votes
2
answers
1k
views
How do i create 1mile buffer for each point geometry and project it to a map?
When ever I try to create buffer, it covers the whole area and a big dot is seen. Here, I am trying to make buffer on points which are located in Travis county, Texas (Austin)
gdf_2017 = gpd....
- 1
-1
votes
1
answer
397
views
Python: correcting offset of points on google maps (due to satellite image tilt)
I'm using the google maps static api to get top view satellite images of objects of which I have the surface coordinates (LoD1 / LoD2).
the points are always slightly off, I think this is due to a ...
- 484
0
votes
1
answer
1k
views
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....
- 7,850
0
votes
0
answers
128
views
Openlayers: Why is (0, 0) in the bottom right for me?
I'm displaying the coords of where the mouse is on the screen, only it appears the coordinates are centered such that the x-axis is negative and the y-axis is positive (0, 0 is at the bottom right of ...
0
votes
1
answer
647
views
Losing the crs when writing to .gpkg with geopandas
When I write my .gpkg I am losing the CRS. I have tried setting the CRS with .set_crs, or adding the CRS when writing the .gpkg (which creates a fault - "fiona._env - WARNING - dataset filename....
- 69
0
votes
1
answer
1k
views
How to convert EPSG:28992 to EPSG:4326?
I have a list of coordinates which i need to convert to longitude and latitude in the folowing way:
What are the math transformations and formulas used to achieve this?
0
votes
1
answer
1k
views
GeoTools transform from CRS in meters to CRS in degree (WGS84)
I am trying to use GeoTools to transform between two coordinate systems, one in meter and the other in degree, but no matter what i try the converted values are not correct. I have tried to follow ...
- 5,576
0
votes
1
answer
1k
views
Transforming coordinates from EPSG:3857 to EPSG:32633 using Proj4js
I'm trying to transform coordinates from EPSG:3857 to EPSG:32633 using Proj4js. I added this script and the code is as follows:
<script src="https://cdnjs.cloudflare.com/ajax/libs/...
0
votes
1
answer
417
views
EPSG:3857 coordinates system compatibility with Azure Maps
Is EPSG:3857 coordinates system compatibility with Azure Maps? I am working on some project which require a geoJson to be loaded over the Azure maps but its not able to identify the coordinates ...
0
votes
0
answers
658
views
QGIS does not change y and x coordinates to decimal degrees
I'm not quite familiar with QGIS, but now I have used too much time trying to find a solution of my probably trivial question.
So I have data for Danish administrative areas. So far I have calculated ...
3
votes
1
answer
1k
views
How can I display the global UTM grid on a map in openlayers?
I want to display the global UTM Grid on a map in openlayers. There is no proj definition for the whole system, just the individual zones. For the projection of my map I want to use EPSG:3857.Im using ...
- 31
0
votes
1
answer
368
views
adehabitatHR home range estimation is too small
I have lat/long data of two animals tracked in Western Australia and I'd like to find their home ranges using adehabitatHR.
library(sp)
library(rgdal)
library(raster)
library(adehabitatHR)
library(sf)
...
- 25
-1
votes
1
answer
869
views
Translate EPSG 3035 to EPSG 4326(lon lat) from shapefile
I have a dataset in a shapefile with CRS: EPSG 3035, which is the extent is xmax,xmin,ymax,ymin. I'm struggling to translate it into long and lat format. Do you know how to solve this?
Thanks
- 3
0
votes
0
answers
772
views
How to convert epsg:28992 to epsg:3857 in python for a shapely polygon
The title says it, how do i change the coordinates of an shapely geom from epsg:28992 to epsg: 3857 in python.
- 77
1
vote
0
answers
471
views
Find lat/lon coordinates using azimuth,inclination and measured_depth(MD) for a wellbore/pipeline
I'm trying to find out the missing lat and lon points for a oil pipeline so that it can be 3D-plotted for our project.I am using a third-party api which provides me with the starting lat and long of ...
- 11
6
votes
8
answers
21k
views
CRSError: Invalid projection: epsg:4326: for geopandas
I am using anaconda for geopandas.
However, everytime I try to use epsg:4326:, it gives an error.
CRSError: Invalid projection: epsg:4326: (Internal Proj Error: proj_create: SQLite error on SELECT ...
- 125
0
votes
1
answer
326
views
How to transform EPSG:6249 coordinates to "standard" lat/lon coordinates?
I have a set of coordinate points in a format refer to as "MAGNA-SIRGAS / Cali urban grid" or "EPSG:6249". Each point coordinates are labeled as North and East. I would like to ...
- 105