New answers tagged coordinate-system
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())
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 +...
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 ...
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 ...
0
votes
Enlarging gBuffer width changes polygon shape
The gBuffer function needed a joinStyle = "MITRE" argument.
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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 ...
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.
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.
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 ...
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 ...
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
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 ...
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=...
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)
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 ...
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..
...
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....
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 ...
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 ...
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, ...
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....
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 ...
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[&...
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 ...
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, ...
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, ...
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 ...
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 !
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 ...
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 ...
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 ...
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 ...
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.
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 ...
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 ...
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 ...
Top 50 recent answers are included
Related Tags
coordinate-system × 7751qgis × 1863
arcgis-desktop × 966
python × 511
raster × 439
gdal × 438
r × 403
coordinates × 371
shapefile × 330
arcmap × 325
proj × 319
wgs84 × 317
postgis × 278
epsg × 267
latitude-longitude × 229
utm × 219
arcpy × 191
convert × 180
datum × 167
georeferencing × 166
openlayers-2 × 156
openlayers × 147
gdalwarp × 138
geoserver × 133
geotiff-tiff × 129