New answers tagged

3 votes
Accepted

Get layer's unit of measurement as string in PyQGIS

You can get it with the following code : QgsUnitTypes.encodeUnit(self.cmbLyrSelect.currentLayer().crs().mapUnits())
user avatar
  • 11.8k
0 votes

What epoch does Proj4 for EPSG:4326 assume? NMEA stream provides Lat Lon at measurement epoch

The WGS 84 ensemble (EPSG:4326 for 2D CRS, EPSG:4979 for 3D CRS and EPSG:4978 for geocentric (ECEF) Cartesian 3D CRS and also some deprecated EPSG codes, but those do not matter) are all static, they ...
1 vote

Custom projection in Centimetre (based on EPSG:28350)

the one that works is this proj4 def (changed units to cm) +proj=utm +zone=50 +south +ellps=GRS80 +units=cm +no_defs based on epsg:28350 proj4 def +proj=utm +zone=50 +south +ellps=GRS80 +units=m +...
user avatar
  • 511
0 votes
Accepted

C# transform coordinates using specific transform code using GDAL

It has to go through a GEOGCS conversion first, in my example from ED50 TM0 (EPSG:23090) > ED50 (EPSG:4230) then run the URN pipeline transformation from ED50 (EPSG:4230) > WGS84 (4326). public ...
user avatar
0 votes

Adding Le Pouce Grid coordinate system to ArcMap

@StainlessSteal, thanks for posting your drawing. The drawing has 2 parts: linework showing the lot and linework showing the site. The lot looks close to where it should be in EPSG:3337, but the site ...
user avatar
  • 319
0 votes

Enlarging gBuffer width changes polygon shape

The gBuffer function needed a joinStyle = "MITRE" argument.
user avatar
  • 9
0 votes

Definition of "SR-ORG" CRS

https://spatialreference.org/about/ “There are currently three distinct catalogs of spatial reference systems on this website -- EPSG, ESRI, and User-defined (called sr-org)” SR-ORG is shorthand for ...
user avatar
3 votes

EPSG:4326 vs EPSG:3857 - Why are they aligned?

As mentioned in the documentation: By default, QGIS reprojects data “on the fly”. What this means is that even if the data itself is in another CRS, QGIS can project it as if it were in a CRS of your ...
user avatar
  • 5,642
1 vote

How to convert crs_wkt starting with GEOGCRS to a ogc_wkt string (to use rioxarray.reproject_match)

Looks like "Pole rotation (netCDF CF convention)" was added in PROJ 8.2: https://github.com/OSGeo/PROJ/pull/2835 https://github.com/pyproj4/pyproj/issues/948 Requires PROJ >= 8.2 to ...
user avatar
  • 5,642
2 votes

Geocoordinates to pixel value in GDAL?

Such a program "gdallocationinfo" exists already and it is one of the standard GDAL utilities. The usage is documented in https://gdal.org/programs/gdallocationinfo.html and the source code ...
user avatar
  • 51.5k
1 vote
Accepted

Extracting the pixel value of a geocoordinate from a GeoTIFF

First inspect the EXIF tags of geoTiff file using gdal info <geotiff_file> and get the coordinates of the corners of your GeoTiff. In case you want to do this programatically read this thread ...
user avatar
  • 536
1 vote

How do I change my view to nadir in QGIS?

I had the wrong CRS. After searching through, and trying several of the WMS services with NADIR imagery mentioned by Vince without significant changes, I tried making a new project altogether and ...
user avatar
1 vote

Using EPSG 4326 in ArcMap

Yes, 4326 is the standard epsg code for latitude/longitude geographic coordinate system using the WGS 84 datum. Looking at the link you included I think you may be confusing geographic vs projected ...
user avatar
  • 2,084
0 votes

Problems with uploading point coordinate British National Grid data into QGIS as .csv file

The "NT 9791" format would be a 1km x 1km square. You can download shapefiles here: https://github.com/charlesroper/OSGB_Grids/find/master As for the 6-digit number, I don't know.
user avatar
  • 319
1 vote
Accepted

Reprojecting from EPSG:3857 to EPSG:28992

I managed to solve the problem by converting from EPSG:4326 instead of from EPSG:3857. My understanding is that EPSG:4326 is a "flattened/map" version of EPSG:3857, so they are equivalent.
user avatar
  • 33
1 vote

GeoPandas input file to shapefile with coordinate system information

Set the CRS when you create the GeoDataframe: gdf = gpd.GeoDataFrame(df, geometry=etc..., crs="the CRS of your xy coordinates") Note, this requires you know the CRS your coordinates are in ...
user avatar
  • 51.8k
0 votes

ArcPy GridIndexFeature doesn't calculate square Grids - But using the tool in ArcGIS Pro desktop it does

You are not using the same parameters in the Desktop application as in the arcpy script. Re-run the script using the same parameters and see if you get the same results. Eg, your arcpy script does not ...
user avatar
3 votes
Accepted

Is there projected version of EPSG:4269?

All of the UTM's across North America have a NAD83 version. https://epsg.io/26914
user avatar
  • 319
2 votes

Is an SRID unique?

EPSG:4326 isn't unique even for EPSG codes. Within the EPSG Geodetic parameter dataset, EPSG:4326 can either refer to a CRS or an Extent. EPSG:4326 (extent) EPSG:4326 (CRS) For uniqueness you need ...
user avatar
  • 11.8k
0 votes

Converting points to latitude and longitude values in R

I've used the following code to make a test raster that prints like yours (but is about 1/100000000 of the size): testlat = raster(matrix(1:100,10,10)) extent(testlat) = c(xmin=13000,xmax=278000,ymin=...
user avatar
  • 51.2k
0 votes

GDAL raster shifting with python error: 'NoneType' object has no attribute 'GetGeoTransform'

To solve that, just put an r before the path: rast_src = gdal.Open(r'C:\temp\raster1.tif', 1)
user avatar
1 vote

What coordinate system is this with North and East?

I think they may be in EPSG::3067, ETRS89 / TM35FIN(E,N), with truncated values. The first coordinate pair (in E,N) may actually be 420000, 7190000. The GKFINxx zones don't cover enough zones to match ...
user avatar
  • 18.2k
0 votes

Imported CSV data points showing in wrong location

I was just about to add on to what @MortenSickel suggested. Exactly how it should go. Glad you have it sorted. I also use 2020/55 for my projections. Food for thought - when dealing with DMS etc.. ...
user avatar
0 votes

Adding Le Pouce Grid coordinate system to ArcMap

I plugged in all the primary and secondary control from the Mauritius Control Network They look to be in the right place. Note: I don't have ArcMap. I used QGIS, but I assume that ArcMap has EPSG:3337....
user avatar
  • 319
4 votes
Accepted

Imported CSV data points showing in wrong location

The steps for getting the projections right in this case: Set the right projection on your csv file in the import dialouge Do not touch the projection of google maps (or any other prepared geodataset ...
user avatar
  • 2,650
4 votes
Accepted

Where are the predefined CRS stored in my QGIS installation?

They aren't saved as a text file, but as a SQLite database file. You can search inside it with a SQL statement, or browsing its data, with DB Browser for SQLite. You can connect to it with the QGIS ...
user avatar
2 votes

CRS is empty but GeoDataFrame plotted correctly when reading shapefile with GeoPandas

The shapefile is a bundle of files. One of them is the .prj. When reading a shapefile with GeoPandas and either: .prj is an empty file no .prj file at all the output will be None. So, in this case, ...
user avatar
  • 23k
1 vote

Convert rworldMap to UTM

This does not change the numbers in the data, it only sets the coordinate system metadata: > projection(map.cp.utm) <- CRS(proj4string(CP.utm)) see how the extent is still in degrees: > map....
user avatar
  • 51.2k
0 votes

Spherical Mercator projection invented by Google and not conformal?

The spherical Mercator projection is a conformal projection from the sphere to the plane but it is not a conformal projection from an (non-spherical) ellipsoid to the plane. However, the errors are ...
user avatar
0 votes

QGIS rotate CRS around a new centre

I suggest to use a hotine oblique mercator WKT PROJCRS["WGS 84 / NZ Oblique Mercator", BASEGEOGCRS["WGS 84", DATUM["World Geodetic System 1984", ELLIPSOID[&...
user avatar
0 votes

projectRaster on a list of rasters

You are reprojecting to the last raster's CRS. It's a bit confusing, I recommend you to: Project to a known CRS (it could be stored in an object as well) Use a raster as template Also, since this ...
user avatar
  • 11.8k
0 votes

Convert coordinates in ArcGIS

I agree with comment above - just take out the commas - ArcGIS thinks that those columns are strings not numbers. I always double check in excel, formatting the columns specifically to be numbers, ...
user avatar
1 vote

If NAD83 is a datum and not a projection, what is happening when set CRS = EPSG:4269 (NAD83)

The proj4 string of EPSG:4269 reveals it all "+proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs ": which means that the "projection" is longlat, so it is unprojected on NAD83 datum, ...
user avatar
1 vote

Convert coordinates in ArcGIS

Mind the language settings of your ArcGIS and how the csv is saved. It is possible that your program does not recognise the Easting and Northing columns as numbers, but as text instead, and therefore ...
user avatar
  • 35
6 votes
Accepted

Can I work with two polygons that were created in a different CRS?

The 'on the fly' CRS transformation allows you to handle layers with different CRS without having to manually change their CRS !
user avatar
  • 1,439
7 votes
Accepted

Google Maps not being quite top-down in QGIS

I assume you have zoomed in quite a bit to see such details. It is a common misconception that you can see satellite imagery in Google Maps' Satellite mode. If you zoom in enough to see details you ...
user avatar
  • 850
2 votes

How to speed up rasterio.transform.xy

It appears like the authors of rasterio know about the slow conversions and have a WIP branch to address this issue. You could implement xy method yourself as the example in the first link ...
user avatar
1 vote

Orthographic projection for South America in ArcMap

You'll have to define a custom coordinate system. Open data frame properties (right click the map or Layers in the table of contents and select properties) and select the Coordinate System tab. Select ...
user avatar
  • 18.2k
1 vote
Accepted

GDAL coordinate transformation using C/C++ API (OGRCoordinateTransformation)

The gdal_retile python script https://gdal.org/programs/gdal_retile.html does not warp. It creates a new tiled file structure in the same coordinate system than the original. The -s_srs option is used ...
user avatar
  • 51.5k
0 votes

Winkel II fails while reprojecting on the fly in QGIS 3.12

According to https://github.com/qgis/QGIS/issues/35512 this will work without issues in any modern QGIS that uses PROJ >= 7.1.
user avatar
  • 8,163
2 votes

UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect

Reproject your data to a coordinate system with meters as units using to_crs. Right now you are buffering with a distance of 0.001 degrees which means the distance will not be the same depending on ...
user avatar
  • 51.4k
3 votes

UserWarning: Geometry is in a geographic CRS. Results from 'buffer' are likely incorrect

You need to understand Geographic coordinate systems and Projected coordinate systems. Geographic CRS(such as 'EPSG:4326') is not suitable for measuring distance. Converting degree to distance ...
user avatar
  • 578
0 votes

Should I use different projections to calculate distance and displaying on the Map?

From what I can gather of your question, you are mixing the display of the map with the output of the calculations, which will be a table. ...the output of the query... is tabular?? Use the ...
user avatar
  • 427

Top 50 recent answers are included