New answers tagged

2 votes

Is EGM2008 different to WGS84 ellipsoid?

In short I have a DEM dataset which describes the vertical reference used as "height above the WGS84 ellipsoid". As I understand it this assumes vertical values to have MSL as a reference ...
mins's user avatar
  • 192
0 votes

How does rioxarray reproject match handle nan values?

This turned out to a be a problem with the initial projection of the source raster (EPSG: 9822 is supposed to be Albers Equal Area Conic but it's not a CRS that is valid in QGIS and looking up the ...
MKF's user avatar
  • 561
1 vote

How many CRS / EPSG codes are there?

This helps explain why at least three of the CRS registers exist: https://committee.iso.org/files/live/sites/tc211/files/Resources/GuideToCRSRegistries3.pdf (Peter Parslow, chair of ISO/TC 211)
Peter Parslow's user avatar
3 votes
Accepted

Assign a CRS to a raster in R

The raster tells you its extent (the coordinate values at the corners) are: extent : 3388038, 3611062, 5266712, 5517788 (xmin, xmax, ymin, ymax) Then you do this: crs(r) <- st_crs(4326)$wkt # ...
Spacedman's user avatar
  • 62.5k
4 votes
Accepted

Is there a comprehensive list of mapnprojections implemented in R?

The projections and projection operations are not implemented in R directly - but in PROJ, which is a standalone library at system level. This library is integrated with R somewhat loosely - you can ...
Jindra Lacko's user avatar
0 votes

Converting EPSG:3346 to WGS84

This looks like a library that you need : https://github.com/dvdoug/PHPCoord
Vincnetas's user avatar
1 vote

Transforming from NAD83(2011) or from WGS84 gives the same result (pyproj)

You get that result because for Proj NAD83(2011) and WGS 84 mean the same. You can get a full report by testing the transformation with Projinfo. projinfo -s epsg:6318 -t epsg:4326 Candidate ...
user30184's user avatar
  • 63k
0 votes

Convert WGS84 to Local Point

You could convert WGS84 to ECEF coordinates, then choose a ref point (should get a matrix 4x4) to convert to ENU coordinates. But WGS84->ECEF is costly, you may directly use mapbox/cheap-ruler to ...
dvorak4tzx's user avatar
3 votes

Compute the area of differents buildings for differents regions in GeoPandas

You can use groupby each row, reproject it using .estimate_utm_crs, and measure the area: import geopandas as gpd file = r"/home/bera/Desktop/GIStest/worldwide_squares.gpkg" df = gpd....
BERA's user avatar
  • 68.7k
0 votes

Converting NZTM (New Zealand Transverse Mercator) to Lat/Long

Created a javascript module "nztm2000-latlng" and released it to NPM repository. The precision is comparable to that of the C implementation. Javascript is not precise enough in floating ...
Zhiwei liu's user avatar
3 votes

Calculating volume of fill needed using QGIS

Solution: reproject You use a CRS that is not suitable for measurements. Reproject to a local CRS, e.g. the UTM zone valid for your area of interest. See also here for backround. Details: why you get ...
Babel's user avatar
  • 67.2k
0 votes

How to set simple square grid as coordinate system in QGIS

All projected coordinate systems are just square xy grids. The cartographic deformation is between the ellipsoid and the grid, but if you don't know their real world location or don't need to locate ...
Gabriel De Luca's user avatar
1 vote

How to set simple square grid as coordinate system in QGIS

I suppose you have a layer grid with an arbitrary cell size and a layer points, both in the same CRS (can be any CRS, does not play a role for this solution to work). On the points layer, use this ...
Babel's user avatar
  • 67.2k
0 votes

ArcGIS Pro 3.0: Changing coordinate system in dynamic text to show XY coordinates not dms, dd etc

Get rid of the "units" tag, to get rid of the "m" in the displayed text Add the "showDirections" tag, do disable the "E" and "N" in the displayed text ...
rmaron's user avatar
  • 1
2 votes

How to produce a Cartopy map in UTM with north-east up?

I can't say how to do it with a UTM projection, but I was able to do something similar using a custom-modified Oblique Mercator projection where I've enabled setting the "+gamma" parameter ...
Ajean's user avatar
  • 131
0 votes
Accepted

Export image to it's nominal scale

I managed to solve it by specifying which image from the ImageCollection I wanted, like so: Map.centerObject(table); // Here I'm trying to set the Bits according to table 'Bitmask for QF_Cloud_Mask' ...
Nikos's user avatar
  • 934
1 vote

Export image to it's nominal scale

In your code, dataset is ee.ImageCollection. Then, you select one band for each image in dataset: var brdf = dataset.select('DNB_BRDF_Corrected_NTL'); so brdf is still ee.ImageCollection. The error ...
Oliver Lopez's user avatar
  • 1,174
2 votes

EPSG for ETRS89 / Lambert Conic in Spain?

Note/Disclosure: I'm a member of the Geodesy subcommittee that maintains the EPSG registry. I work for Esri as my day job. The short answer is that no one has asked for it. The geodesy subcommittee--...
mkennedy's user avatar
  • 18.8k
2 votes
Accepted

Differences in extracted values after projecting raster

When you use the "pipe" operator, it puts the thing on the left as the first argument of the function call on the right. You are doing: NDVI_2005 <- NDVI_2005 %>% projectRaster(...
Spacedman's user avatar
  • 62.5k
1 vote
Accepted

Overlaying map shaper geometries

To convert from a shapefile already projected to EPSG:27700 into TopoJSON using WGS84 projection I found the solution was to use init= as described here my example became the following. mapshaper map....
Daniel Harper's user avatar
0 votes

How to rotating a DXF world map along the longitude using QGIS?

The issue here was that I'm missing georeferencing information, the Reproject Layer tool works fine as expected.
Luiz Borges's user avatar
0 votes

How to rotating a DXF world map along the longitude using QGIS?

Please try in QGIS, on the General tab of the Layer Properties, there is an Update Extents button
Gurminder Bharani's user avatar
6 votes

Orientation differs between projections. What is the "Real" North?

As others have pointed out, orientations measured in GIS will always be relative to "grid north" for a particular projection. The relationship between grid north and true north varies ...
lambertj's user avatar
  • 2,962
4 votes

Orientation differs between projections. What is the "Real" North?

The North of a projection refers to the projection's own definition. Projections can be rotated or otherwise distorted in relation to the geometric body they project to a plane. If you need ...
bugmenot123's user avatar
  • 10.7k
-1 votes
Accepted

Incorrect transformation being used in GeoPandas?

Using EPSG:8999 or EPSG:7911 instead of EPSG:4326 makes the "shifted" data align with the "correct" data. That is really frustrating and annoying, both in ArcGIS and GeoPandas/...
Benny Jobigan's user avatar
0 votes

Importing DWG issue

You can do this by using a custom coordinate system (CRS) for these layers. Select option Custom projections from the Settings-menu. There: Add a new CRS (+ button) Enter a relevant name, say: BNGmm ...
peter's user avatar
  • 765
0 votes

Define latitude longitude grid from center point

Reproject to a projected coordinate system. Then create a point grid, square buffer and go back to lat/long: import geopandas as gpd import shapely import numpy as np lat, lon = 50, -10 # my point ...
BERA's user avatar
  • 68.7k
2 votes
Accepted

Issue with GPKG imported in QGIS

It seems the layer's CRS is not correctly recognized. When loading the layer to an empty project, I get the same error you describe in the comments. To solve it, proceed as follows: Create a new, ...
Babel's user avatar
  • 67.2k
2 votes
Accepted

autonavi / amap xyz tile layer appears shifted in OpenLayers (Chinese Web Mercator with GCJ-02 datum)

The Baidu layer will align with OSM (but not using https) by changing var y = tileCoord[2]; to var y = -tileCoord[2] - 1; (see https://github.com/openlayers/openlayers/issues/10374) For the amap layer ...
Mike's user avatar
  • 12.3k

Top 50 recent answers are included