I have made a raster elevation map of Greenland, and have information of the region, cities and more as vectors. The information is not in the same place as my elevation map.  I have selected EPSG:5938 on the bottom right in QGIS and on my informations and elevation map I have set the CRS to the same like the picture    How do I get it all on the same position?
120503
To start off with, I'm very new to QGIS. I'm attempting to place a shapefile with locations for 3 substations on a map of Peniche, Portugal. The locations for the substations are given in latitude and longitude coordinates.   When I start off by importing the substation shapefile, I make sure that both the layer CRS and project CRS are set to EPSG:32629, which gives each substation point the correct coordinates, however the scale is at 1:1, which is not correct. As soon as I then add a raster map of the region or a shapefile of Portugal using the same CRS the position of the substations also becomes incorrect.   I'm sure there's a simple solution I'm not getting, but I've been trying for a while and have had no success. Here's a dropbox link to the files
86909
I am in QGIS. I have a raster layer of mountains and a lake. I have a vector layer (shapefile) of the lake. It sits nicely where it supposed to be over the raster level. I am going to convert this height map to an STL so that I can cnc this out of wood, so I'd like the actual lake to be something like 200 feet lower than it shows on the raster layer so that I can fill it with blue epoxy. I have read a lot of clipping and everything else, and I am lost.
100956
I've generated a raster DEM in ACADCIVIL in TIFF format: name.tif &amp; name.tif.aux.xml. ArcGIS reads the data correctly but QGIS sets the whole data to the same value = minimum height Looks that QGIS 3.4 doesn't support this format. Same happens on 2.18  gdalinfo output:  Driver: GTiff/GeoTIFF Files: Surface1.tif        Surface1.tif.aux.xml Size is 405, 528 Coordinate System is: LOCAL_CS["POSGAR_07_F5",     UNIT["metre",1,         AUTHORITY["EPSG","9001"]]] Origin = (5512262.832999999634922,6109047.309000000357628) Pixel Size = (1.000000000000000,-1.000000000000000) Metadata:   AREA_OR_POINT=Area   TIFFTAG_SOFTWARE=Autodesk Civil3D 2008 Image Structure Metadata:   INTERLEAVE=BAND Corner Coordinates: Upper Left  ( 5512262.833, 6109047.309)  Lower Left  ( 5512262.833, 6108519.309)  Upper Right ( 5512667.833, 6109047.309)  Lower Right ( 5512667.833, 6108519.309)  Center      ( 5512465.333, 6108783.309)  Band 1 Block=405x100 Type=Float32, ColorInterp=Gray   Min=38.088 Max=39.473    Minimum=38.088, Maximum=39.473, Mean=38.565, StdDev=0.301   NoData Value=3.4028234663852886e+38   Offset: 38.082,   Scale:0   Metadata:     STATISTICS_COVARIANCES=9.08270199017587E-02     STATISTICS_MAXIMUM=39.473079681396     STATISTICS_MEAN=38.564526756118     STATISTICS_MINIMUM=38.088008880615     STATISTICS_STDDEV=0.30137388124225     STATISTICS_VALID_PERCENT=52.86   I've changed scale to 1 --> same result File Properties --> Info Band  Number Band    No-Data    Min   Max 1    Band 1  3.40282e+38  n/a   n/a
86533
I'm new to QGIS and working with GIS in general (student of architecture), so please bear with me.  I have obtained shapefiles with topographic curves of an ocean floor. My intention is to export these as CAD files and render as 3D models of the ocean floor in Revit. I have done this procedure earlier with terrain shapefiles. My problem now is that these files mapping depth does not contain z values, so when exporting them into CAD they are completely flat rather than 3D models.  The depth information in meter exists as attributes in the file (see attached image). What I want is to transform this depth value into real z coordinates for each line so that the curves have different heights when importing as CAD file in Revit.  Thankful for any pointers in the right direction.
31731
I have a project and a topomap in WGS84 (EPSG4326). On the fly CRS transformation is activated. I open Google satellite layer from the OpenLayers plugin. The Google layer is WGS84 / Pseudo Mercator (EPSG3857). But then the whole project has also become WGS84 / Pseudo Mercator (EPSG3857).  I thought on the fly transformation would reproject to EPSG4326. Why is that? Is it better to work directly in WGS84 / Pseudo Mercator (EPSG3857) systematically when we use Google Earth?
106025
The problem   I have tens of thousands of small polygons for which I want to create buffer/doughnuts for, but I don’t get the units right. I want meters and not degrees, and changing to (what I think is) the right CRS my layer doesn’t longer show. The polygons are small farming plots. I work in QGIS 3.4.15.   Let me take you through it  I start by adding a vector layer. I opt for ‘Directory’ as a source type (Encoding: UTF-8 by default) and under Source (Type: UK. NTF2 by default) I point to the folder in which I store the shapefile.   The default CRS is set to EPSG:4326 – WGS 84, for both layer and the project (checked in bottom right), and unit is thus in degrees. I add the layer world map (EPSG:4326 – WGS 84 by default) by taping “world” in the coordinate box in the middle bottom, I also add Open street map (EPSG:3857 – WGS 84 / Pseudo-Mercator by default) from XYZ Tiles. The polygons appear nicely on the map considering some are just adjacent to the nearby lake, many nicely adjacent to roads, etc.   When attempting to do buffers the distance is obviously set to degrees. I have tried to change the CRS, both for the layer(s) and the projects, but without luck. After researching (e.g. here: ) I believe that the correct CSR would be Arc 1960 / UTM zone 36S EPSG:21036 (it is my understanding that either 36N or 36S would work as the polygons appear both south and north of, but close enough to the equator). When setting the mentioned CRS for my polygon layer, the polygons disappear from the map, and remain gone even after setting the project CSR (by bottom right) to the same. It seems that whenever I try with different CRS in UTM, the layer doesn’t show. Does this mean that the polygons are somehow not compatible with units in meter? Is there an issue with importing the .prj file?  Other things I've tried:   I have made a good attempt searching for solutions. I have ensured that the polygon layer is on top of the layers list. I have ensured the file names are identical and tried all different suitable CRS I can think of. I have seen similar questions been solved by saving the .shp 'as', but I don't think that's my issue as it's already a .shp file. Can it be the case that I need to perform a datum transformation? I gave it a shot using a transformation code (1122) provided by above link, but then the polygons were projected on the Antarctica on the world map, so obviously something went quite wrong...   Note that when I have changed to what I believe is the right CRS I can execute buffers in units of meters, however, they do not appear on the window…
81880
I received several shapefiles from an colleague that he told me is in WGS:84. If I look at the metadata in Qgis it tells me the following:  In layer spatial reference system units xMin,yMin 3572722.82,5453423.31 : xMax,yMax 3582448.25,5462926.22 Layer Spatial Reference System +proj=longlat +datum=WGS84 +no_defs   Now when I try to save it in another projection say WGS:84 Pseudo Mercator (EPSG: 3857) it gives me the following:   Export to vector file failed.     Error: Failed to transform a point while drawing a feature with ID '1'. Writing stopped. (Exception: forward transform of     (62425.896446, 95242.109906)     PROJ.4: +proj=longlat +datum=WGS84 +no_defs +to +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs     Error: latitude or longitude exceeded limits)   A similar message is returned in R with the rgdal package. I suspect that the initial projection is somehow wrong.   Is it truly WGS84?   The problem occurs with all shapefiles.  ogrinfo -so    tells me:  INFO: Open of `XXX.shp'       using driver `ESRI Shapefile' successful. 1: XXXABC (Polygon)
76218
I have a script which contains several r objects which are CRSs that I use for transforming spatial objects. After updating R and many packages, there is a warning which wasn't there before. LCC &lt;- CRS(&quot;+proj=lcc +lat_1=49 +lat_2=77 +lat_0=0 +lon_0=-95 +x_0=0 +y_0=0 +ellps=GRS80 +units=m +no_defs&quot;)  Warning message: In showSRID(uprojargs, format = &quot;PROJ&quot;, multiline = &quot;NO&quot;, prefer_proj = prefer_proj) :   Discarded datum Unknown based on GRS80 ellipsoid in CRS definition  I understand that this is due to a change in how spatial information is stored... PROJ4 to PROJ6 upgrade and &quot;Discarded datum&quot; warnings ..but beyond that, I don't understand what this error means or if I need to do anything in response. Any thoughts?
62433
I have a point layer and a polygon layer that are all in UTM51N projection but do not align. Just want to know when should I use Assign, Reproject, or Define.  Experimented with a polygon in EPSG:4326 - WGS 84 - Geographic, which is the project CRS since this polygon is the first layer I loaded. Target projection is UTM51N.  Reproject - aligned  Assign - not aligned  Define - not aligned  Save layer as - aligned  Using 3.4 so there is no "on the fly projection" feature.
96879
How do I change the coordinate value of a vector layer (13293061,-572595) to 119.413,-5.13677 (vector data &amp; Google Earth do not match)?   If I output my data using Google Earth as background, the vector data does not match with the underlying google layer.   Printing the result to PDF also has this same problem.
29176
I am working on QGIS. I have a raster of EVI, in which values are comprised between 0 and 1. I need to vectorize this raster in order to make an intersection with another vector. My issue is that when I vectorize my raster, the values are either 0 or 1, but I want the values detailed between 0 and 1. The best would be to have the values every 0.1 (0.1, 0.2, 0.3, ...). Do you have any idea of how I can proceed? I was thinking that I could multiply the values of the raster by a hundred, but is that even possible?
39146
I learned on ArcGIS, and I'm familiarizing myself with the use of QGIS.  So far I love it, but the process of re-projecting and transforming datums seems less intuitive than in Arc.    I recently geocoded 63 locations, pulled in the csv and transformed into a points shapefile. I then set the projection for that points layer as the same for my block groups multi-polygon (EPSG 102121- NAD_1983_Michigan_GeoRef_Feet_US).    The 2 layers are still misaligned, even if I save each layer with the specified projection and try to reopen in a blank project.  Can anyone help?
72125
Excuse me if this is a stupid question, but I've honestly tried to find the answer in a few threads here on SE, and I cannot seem to follow the suggested steps correctly.    Project is in  (WGS 84 / Pseudo-Mercator) BaseMap (OpenStreetMap) layer is set to EPSG:3857 Vector data has been imported from CSV, and has Lat/Long coordinates, set to EPSG:3857.    The layer appears right in the middle of the map, around 0,0 obviously.   Why is this, what am I doing wrong?
44385
I am using Windows 7 enterprise and QGIS version 2.12.2 Lyon.  Have found an unexpected result when I use QGIS to generate a topographic position index (TPI) raster from a 30m by 30m DEM.  I find a pattern of what I believe to be aberrantly valued pixels forming a square grid pattern that are approximately 900 m on each side.  I expected the resulting raster to not have this grided appearance.  UPDATE: I have looked at my hillshade, slope, aspect, terrain ruggedness index and roughness rasters generated using the same tool and they all have this same "grid" (although it is much less apparent than in the TPI raster).  Here's a screenshot of the original DEM I used:   Here's a screenshot of the resulting TPI.  Notice the grid pattern formed over what I suspect to be an accurate TPI raster layer.   Here's how I created the TPI Raster: Raster – analysis – DEM:
91633
I have a shapefile and I loaded it into QGIS using Add vector layer. After that, I can see points in different positions. I open the Attribute table and I can see columns &quot;X&quot; and &quot;Y&quot; with values such as (151449.965057, 204642.790761). Now I want to convert them to latitude and longitude. How can I do that? Therefore, I use the option Save as to export the data to csv format using the CRS WGS 84. However, nothing changes in the CSV file (&quot;X&quot; and &quot;Y&quot; stay the same as the original data).  I read and tried several suggestions from this forum without help, e.g., .
40823
I'm fairly new to QGIS, but I'm getting the hang of it. I have two layers - one DEM and one raster made from shape file (Buildings), but need one with higher Elevation where the buildings are for terraflow. But I can't seem to get the raster calculator to work with if Statement. I found a Workaround by assigning a fixed elevation value to the building raster, that is higher than any Point on DEM.  But on DEMs with high Elevation difference this means that some Buildings are over 50m tall!!!  would like to make an if statement like this:  A  B A = DEM  B = Buildings  if(B&gt;A,A+10,A)
113684
In qgis in the raster menu there is an item "convert vector to raster".  This behaves strangely - I do not understand how it is supposed to work.  I can select a vector file, an attribute and an output raster. But, for some reason the output file is only accepted if it exists and even if I select an existing raster nothing seems to happen.  Also I do not understand how the program knows what kind of raster to create, what cell-size, projection etc.  How can I convert vector to raster?