Given a .tiff topographic raster image with mountains and deep sea, I am looking for the x,y,z specifics of the highest or lowest point.  The x,y,z values being:   altitude (z-coordinate) latitude (y-coordinate) longitude (x-coordinate)   While gdalinfo raster.tif would give a summary of the bounding box and the highest / lowest elevation values , and note comments below), I would like to find the elevation, latitude, longitude of the highest and/or lowest pixels of my raster.  Would be better if these values are got as isolated output integer. I think raster processing iterating on all pixels of the raster and keeping values of the z-highest point may be a good way to go.
120739
Given a .tiff topographic raster image.   How to get the altitudes of the highest and lowest point/pixel ?
117666
I want to process a raster with values 0-16. I want to set all values between 1 and 15 as 1, a value of 16 as 0 and cells with a value of 0 as 0.   How do I maintain a cell value of 0?  I have tried the following:  ("raster_name@1" !=0 AND "raster_name@1" &lt;= 15*1)+("raster_name@1" = 16*0)   Results in: Values of 16 as 0, values 1-15 as 1, values 0 also = 1  I also tried:  ("raster_name@1" !=0 AND "raster_name@1" &lt;= 15*1)+("raster_name@1" 16*0)+("raster_name@1" = 0*0)   Results in: Values of 16 as 0, values 1-15 as 1, values 0 = 2  Please, can you explain how to retain cells with a 0 value as 0?
90058
How do I find the coordinates of the centre / centroid of a raster?   I have tried Zonal Geometry but this does not give me the coordinates.  I need to find the centroids of two rasters and create points and measure the distance and bearings between them. i.e. distance and direction between two raster centroids.  I am using ArcGIS Desktop 10.2.
105638
I have around 1,000 satellite images in tiff format, and I want to create a shapefile which will serve as an index to the rasters. This is something similar to a raster catalog, but I do not want to build a raster catalog.  Some obstacles I can forsee, is that the image are georefrenced, so they are not rectangular in shape (I am talking about the data area).  To clarify, I require the polygon to cover only the non-zero (or non nodata) pixels of the raster, and not the entire rectangular raster. Most of the answers so far, give a rectangular polygon, which covers the data, as well as the non-data pixels.  My Image     Result given by  tools I have examined (like raster catlog, various Arcscripts, custom Python script given in one of the answers):     Result that I want:
3759
i am very new to qgis. I want to create a hillshade file and also calculate the inclination.  The structure of the digital elevation model data looks like this:  5390.0125 5659.9875 251.8 5390.0125 5659.9625 250.6 5390.0125 5659.9375 247.5 5390.0125 5659.9125 240.6 5390.0125 5659.8875 236.1 5390.0125 5659.8625 234.8 5390.0125 5659.8375 235.6   The file format is *.dat. The coordinate system is Gauss Krüger, Bessel, 5. meridian.  I have no idea how to import the data in qgis.
39834
I have downloaded SRTM GDEM (~90m resolution). I am using ArcGIS 10. I have tried to use the Spatial Analyst tool to compute for slope. However, I cannot compute for the slope. The output values has only two ranges 0 and 0.1-90. I am not really sure what the problem is.
6180
I have a GeoTIFF published by Geoserver 2.2 and being used by Leaflet via WMS. How can I get Geoserver to return the raster value at a lat/long point?   I'm guessing it has to do with using WCS, but no examples seem to exist for doing this!  This will be similar to querying PostGIS with ST_Value(raster, ST_SetSRID(ST_Point(lon,lat),4326)). I dont really want to import the same raster into PostGIS just to find the raster value at a point. Or is this the recommended solution?
114635
I was wondering if anyone knows of a plugin, extension, process, etc. for mapping shadow locations based on date/time using raster data? I have searched and what I have been able to find all seems to require vectors or multipatch data. I have DSM rasters produced from LiDAR and would ultimately like to use those as the source for identifying locations of shadows, as I would like to the surrounding features in the raster to be used as sources of impact (i.e. trees). I am open to all options, but QGIS, Grass, and ArcGIS are preferred.
63804
I have a raster height model and would like to calculate the horizontal geographic direction (north, east, etc. - orientation) of each cell, then calculate then the mean value of all these cells in one region (polygon). This is to get information about to what direction a vineyard is exposed to. How do I do this?
44745
I have been attempting to use the Spatial Analyst Curvature tool with mixed results.  I originally tried to use it with a raster that was projected using a 1983 Teale Albers projection and the results were black, indicating the lowest values.  I then transformed the raster into UTM projection and ended up with grid-like lines on the raster which have values that are obviously not accurate.  Does anyone know why these lines have appeared or what I'm doing wrong.  Thanks.
23912
I have two layers in ArcGIS: One is a vector point and the other is a raster layer. I want to append the value of the raster data grid to the vector point. The problem is the vector point is away from any existing raster grids with value, so the command "Sample" return no value for those vector point.    So my question is: How can I find the nearest raster cell given a vector point and extract the value from the cell?   Actually there are two things I need to do:   Identify the outliner within a ranger (Usually 1-2 cells away from my existing raster) Assign a value from the nearest raster to them, if they are in a range     Edited: I have about 3000 vector points and my problem is this:   I need to append the raster value based on locations to the vector points. I used the tool "Sample" and it works well for most of the points.  However, there is a situation like this   The right point can get the value from the raster ("Sample" works) but the left point cannot due to the alignment issues.  Raster to polygon is not working because I get error message that I cannot turn them into vector because they are out of the domain. Also I have points like this  Those upper points should not get any values because they are too far away from the raster.  I thought of someways to work around:    Step 1. Do the "sample" first    Step 2. Select the null value after the "sample" operations  Step 3. construct buffer zones (within 1 decimal degree) based on the null points    Step 4. use zonal statistics??? or other tools to loop through all the raster cells in the buffer, find the nearest one, extract its value and put it in the vector point.    I am stucked at step 4. I don't know which tool in ArcSDK I can use to achieve this function.   Or else...Do you have better idea to deal with that?    I have 20 raster layers like this and I would like to make it automatically (using model builder and arcobject).
51826
I am trying to loop through multiple rasters and extract cells that are above 0.1 The code I am trying to use is:  import arcpy  arcpy.CheckOutExtension('Spatial') arcpy.env.workspace=r'F:\Sheyenne\Atmospherically Corrected Landsat\Indices\Main\NDWI\main_master' #pathway to all rasters in workspace directory rasters = arcpy.ListRasters('*.tif*') for raster in rasters:      #give new file a name      outraster = raster.replace('.tif','_water.tif')      arcpy.gp.RasterCalculator_sa("""raster &gt;= 0.1""",outraster) print('Done Processing')   For some reason I can't copy and paste the error (python actually shuts down when I run this code) but here is a screen shot of part of it.    EDIT: My parameters have changed and I now need everything >=-0.4 and but !=0.00.  I am trying to use:  import arcpy  from arcpy.sa import * arcpy.CheckOutExtension('Spatial') arcpy.env.workspace=r'D:\Sheyenne\Atmospherically Corrected Landsat\Indices\Main\NDWI\main_master' #pathway to all rasters in workspace directory rasters = arcpy.ListRasters('*.tif*') for raster in rasters:     #Save temp raster to disk with new name     ras = Raster(raster)     outraster = Con(ras &gt;= -0.4 &amp; ras !=0.00, ras)     outraster.save(raster.replace('.tif','_water.tif')) print('Done Processing')   but this returns:  ValueError: The truth value of a raster is ambiguous. Invalid use of raster with Boolean operator or function. Check the use of parentheses where applicable.   EDIT:  I think I just needed to change this line:   outraster = Con((ras &gt;= -0.4) &amp; (ras !=0.00), ras)
104815
If I have a data set that produces a graph such as the following, how would I algorithmically determine the x-values of the peaks shown (in this case three of them):
16009
I have downloaded Globcover data in GeoTIFF format. I want to create a new grid which only has values where the Globcover data (the only band in the TIF file) has one certain value (==190).  Can I do this with GDAL?  If not, which way should I go? The procedure should be callable from a Shell script (Linux). I have QGIS installed (but no clue how to use it).
99349
I have created a point layer with ~150 points. I have a DEM beneath it, which is where I have taken the X and Y values from.   Does anybody know how I can retrieve the Z values, without using the info button?  I am using ArcGIS Desktop 10. I have licences for the spatial analyst etc tools.
33940
I'm new to ArcGIS Desktop. I have created a layer representing zip codes and overlayed it with a layer containing some data points. I want to find the min and max value of the data points in each zip code. Currently, I am manually checking the values if all the data points in each zip code but this is becoming very time-consuming and cumbersome.   Is there a faster way to do this?
115395
I'm looking for a tool or a plug-in for QGIS that would allow me to shift rasters to new coordinates. An equivalent of ArcGIS'  tool would be the best option.  I know I can edit the world file that comes along with the raster to do this, but some raster formats (like geotiff) contain their metadata inside their binary format.