266 questions
1
vote
1
answer
64
views
OGRCoordinateTransformation.Transform() gives wrong coordinates
I'm trying to do some local to global (EPSG:3042 to EPSG:4326) transformations in a ROS2 c++ script using GDAL.
OGRCoordinateTransformation *Transformation;
OGRSpatialReference poLocal, poGlobal;
...
2
votes
0
answers
182
views
"PROJ: proj_create_from_database" error when transferring a KML file to PostgreSQL using ogr2ogr
I am trying to import a KML file into a PostgreSQL database using the ogr2ogr tool, but I am encountering the following error:
ERROR 1: PROJ: proj_create_from_database: C:\Program Files\PostgreSQL\17\...
3
votes
2
answers
168
views
C++ libproj, swap coordinates
I have a C++ program running on linux for years. It uses some coordinate transformation, we use the libproj (version 6.3.1) for this purpose.
Interestingly, one client gave the input parameter epsg:...
- 1,675
1
vote
1
answer
69
views
Differences in CRS when reading in NetCDF using terra::spatRaster
I've noticed a difference in behaviour in terra::rast on two different computers running:
Windows; R4.4; terra 1.71-78; GEOS 3.12.1, GDAL 3.8.4, PROJ 9.3.1
Linux; R4.4; terra 1.71-71; 3.10.1, GDAL 3....
- 13
1
vote
0
answers
200
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
1
vote
0
answers
76
views
I need assistance with the error Proj_db in QGIS
I'm trying to merge 3 rasters to complete the map of a country but when i run the merge function on QGIS i obtain this error:
GDAL command:gdal_merge.py -separate -ot Float32 -of GTiff -co COMPRESS=...
1
vote
0
answers
307
views
How to crop an xarray Dataset using a Shapefile with non-standard grid coordinates
I have a netCDF dataset with dimensions represented by 'xc' and 'yc', which are in meters, and the CRS of the dataset is Lambert Azimuthal Equal Area projection with the following parameters:
Grid ...
- 133
1
vote
0
answers
53
views
How do I use inverse projection in gdalwarp?
I have a small orthographic image I want in lat-long (equirectangular). Using gdalwarp to reproject causes a truncation at the poles above 80 degrees (feature or bug?). Does unprojecting using the -...
- 153
1
vote
3
answers
325
views
Highcharts US Bubble Map: Tooltip Issue When Hovering Over Bubbles
I'm using Highcharts in a React application to create a US Bubble Map. The version details for the relevant libraries are as follows:
"@highcharts/map-collection": "^2.1.0"
"...
1
vote
1
answer
4k
views
Warning: GDAL Error '1' and Error: Node stack overflow
OS: Windows
I am plotting a raster file after changing the CRS of the same with st_transform. Here, is a warning of GDAL Error 1 being generated.
Later, I changed the CRS of the vector/feature file ...
- 129
1
vote
0
answers
158
views
Unable to install Pyproj on a raspberry pi zero 2W (debian bullseye)
I tried to install Pyproj=3.0.0 and can't get it to work. The most weird thing is, that I was able to install it on a Raspi 1 B in the past. Sadly I can't remember how I did it and now I'm not able to ...
- 11
0
votes
1
answer
1k
views
How does one build GDAL 3.7.2 binaries
I'm trying to build the GDAL DLL and .lib files, so GDAL can be used for linker and includes in an other C++ project. We're needing to upgrade from GDAL 2.0.2 to 3.7.2.
Has anyone built GDAL from ...
- 3
0
votes
0
answers
57
views
I am getting incorrect north and east values after carthographic projection using PROJ library in c++
I am getting incorrect north and east values after carthographic projection from global wgs84 crs(i.e., EPSG:4326) to crs Korean 1985/west belt using proj c++ library. I got an error of approx. 1.5 cm ...
0
votes
1
answer
1k
views
Resample and plot satpy scene data from geos to mercator using pyresample
I have a scene object from eumetsat satellite.
My area of interest:
(xmin_corner, ymin_corner, xmax_corner, ymax_corner)
[40,-10,105,40] degrees
looks like this
area definition looks like this:
Area ...
- 117
0
votes
0
answers
100
views
Issues while interfacing R with GRASS GIS using rgrass package
I want to use some GRASS-GIS tools in R and I am using the rgrass package to interface R to GRASS GIS installed on my windows machine.
When I run initGRASS(gisBase = "C:/Program Files/GRASS GIS 8....
1
vote
1
answer
83
views
Difference in Z result of ST_Transform
Query:
with geom as (
select st_transform(
st_geometryfromtext(
'POINT Z(114.19515060316667 22.3979100265 81.8574198829997)',4326) ,2326)as a)
select st_x(a),st_y(a),st_z(a) ...
- 11
1
vote
0
answers
986
views
Error: PROJ: proj_create: unrecognized format / unknown name
I have installed Terria JS and Cartoview. when i upload layers in Cartoview then they are uploaded but the following errors come:
Error Message: PROJ: proj_create_from_database: SQLite error on ...
0
votes
1
answer
1k
views
Using PROJ C++ API for CRS-to-CRS transformations
pyproj is a language binding for PROJ. The pyproj module provides easy-to-use methods for CRS-to-CRS transformations, For instance, use it to convert global latitude/longitude (degrees) to local ...
- 45
4
votes
0
answers
1k
views
Trying to build GDAL from github using Visual Studio 2019 on windows 64
I want to use GDAL for a project that I am working on which is written in C++ using Visual Studio 16 2019. I am getting errors and I am finding difficulties in building the library. The steps I ...
- 51
0
votes
1
answer
118
views
pyproj import fails with array.array error when running from command line, macOS, 13.0.1, Python 3.8.10, pyproj 3.4.0
System is configured as macOS 13.0.1 (Ventura), python 3.8.10, pyproj-3.4.0.
python is installed from python.org
>which python3.8
/Library/Frameworks/Python.framework/Versions/3.8/bin/python3.8
&...
1
vote
1
answer
709
views
Using World Files (jgw, pgw) in Leaflet
I have a list of image files along with their corresponding world files in jgw, pgw format etc. how can I use the world files to get the real world coordinates of these images so I can generate them ...
- 43
0
votes
0
answers
361
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
162
views
Creating a spherical CRS in Pyproj with "R_A"
The Proj documentation states that the R_A parameter indicates "A sphere with the same surface area as the ellipsoid." (source)
So I tried this, believing I would get a spherical CRS with ...
- 14.1k
2
votes
1
answer
270
views
Recreate sumo projection
After converting osm file to sumo network I got this info
<location netOffset="-395885.58,-6185860.67" convBoundary="0.00,0.00,3442.30,2845.95" origBoundary="37.338524,55....
- 23
5
votes
4
answers
4k
views
Can't install Proj 8.0.0 for cartopy linux
I am trying to install Cartopy on Ubuntu and need to install proj v8.0.0 binaries for Cartopy. However when I try to apt-get install proj-bin I can only get proj v6.3.1. How do I install the latest (...
- 1,979