Raster is a data format consisting of regular grids of values, usually stored in an image-like format.

learn more… | top users | synonyms

10
votes
4answers
1k views

Problems with NA values when reading .DEM file with R 'raster' package in Windows

I am trying to read a raster file in a .DEM format on windows using the 'raster' package in R. I get problems with NA values, when loading the data into R in Windows 7, but I do not have the problem ...
10
votes
6answers
1k views

How do I iterate through every cell in a continuous raster?

See this link for more details. The Problem: I want to loop through a continuous raster (one that has no attribute table), cell by cell, and get the value of the cell. I want to take those values ...
10
votes
5answers
1k views

Multicriteria Analysis : how to?

I have to do a multi-criteria analysis to answer the question : "which is the best lot to develop". A few of the criterias are : distance of the nearest bus stop (point layer with bus stops) ...
10
votes
2answers
2k views

Grid vs TIF vs IMG

As the title suggests both of the above three formats are raster files. In ArcGIS, apparently there's no difference between them. So what's the difference between these three formats?
10
votes
1answer
447 views

How do I randomly alter a raster map of habitat types?

I have a raster of habitat types for a specific area in Scotland. I need to create future habitat scenarios with changes in habitat to assess the population viability of a bird species. For example, ...
10
votes
3answers
1k views

How to speed up queries for raster databases?

I have a raster database in postgresql/postgis with these columns: (ID, rast, data_of_data). 'rast' is the column that has raster files in WKT format. An example query to find the DN value of a ...
10
votes
2answers
2k views

Aggregating points to grid using R

I have a question with regard to spatial aggregation in R. What I am trying to do is aggregate a point dataset to a grid. I am unsure however how to do this as I have little experience with this sort ...
10
votes
2answers
1k views

PostGIS raster value of a lat/lon point

I have a tiff raster in my PostGIS 2.0 database. I want to get the value of the raster at a specified lat/lon point. I have: postgresql-9.1 postgis-2.0-svn with raster support a raster imported ...
10
votes
3answers
324 views

From polygon to continuous measure of diversity?

First of all; I have tried to find a similar question, without success. Maybe it's because I'm quite new to GIS and I don't really know what exactly I'm looking for. If someone points me to a similar ...
10
votes
2answers
276 views

Count raster value in a bounding box (spatial extent)

I have simple question: given a bounding box, how can I count the number of specific values (say, I am interested in the number of value == 1 ) in a raster in 1) arcpy, or 2) raster package in R?
10
votes
1answer
802 views

What is some theory on the number of GCP when georeferencing?

Just wondering, is it always better to have more than 3 GCPs? Should the GCPs be spread out evenly? Where can I find some general resources on the theory of this? Other than trial an error, which ...
10
votes
2answers
423 views

What are the advantages/disadvantages of populating a PostGIS database with Raster information?

...as opposed to keeping the rasters as locally stored GEOTIFFs? I have a large project that will require making available approximetely 50,000, 10mb GEOTIFFs as WMS to an OpenLayers based web map. ...
10
votes
3answers
2k views

How can I mosaic images in ArcMap without it affecting the colour?

This is 4 images (NE, SE, NW, SW) sections of imagery. They have some overlap between them. Working with 1m 3-band orthoimagery. Different parts are sometime flown at different times/sources hence ...
10
votes
2answers
2k views

Intersecting a raster with a polygon using PostGIS - artefact error

I'm using PostGIS2.0 to do some raster/polygon intersections. I'm having difficulty understanding which operation I should use, and what the quickest way of performing this is. My problem is as ...
10
votes
2answers
1k views

How can I improve my Georeferencing results?

Background This is my second question related to georeferencing naked raster maps in order to re-visualize them on different coordinate systems and in conjunction with other data layers. The previous ...
9
votes
7answers
2k views

how to overlay shapefile and raster?

I have a shapefile with polygons. And I have a global raster file. I want to overlay the shapefile's polygons onto the raster grid and calculate the mean raster value for each polygon. How can I do ...
9
votes
7answers
18k views

Get boundary of raster image as polygon in ArcGIS Desktop?

I have a rectangular raster image in ArcGIS. My question is: how can I produce a polygon which is just the extent of the raster image? That is, I want to have a polygon layer with a single ...
9
votes
4answers
7k views

How to calculate raster statistics for polygons?

how can I calculate sums, averages etc of raster-points (multi-band) per polygon of a vector-layer. I was told that this is called "zonal statistics". I tried that with QGIS first. There is a way to ...
9
votes
3answers
5k views

How to fully load a raster into a numpy array?

I have been trying to check my filters on DEM raster for pattern recognition and it is always resulting in missing last rows and columns(like..20). I have tried with PIL library, image load. Then with ...
9
votes
6answers
6k views

Rasterfile clipping error: Cannot compute bounding box of cutline

I have a raster file and want to clip it, a vector file is supposed to be the mask layer. I use the most obvious method (Raster->Extraction->Clipper). No matter what I do, I get this error (see pic): ...
9
votes
3answers
3k views

Best method to aggregate 1m DEM mosaic to 10m DEM

I have a mosaic of a 1m LiDar derived DEM. I need to output a subset of the data as a 10m DEM. I am currently using the aggregate tool in ARCGIS 10 to produce a mean value for each new 10m pixel. ...
9
votes
2answers
2k views

Difference between Geoserver stores: ImageMosaic, WorldImage, ImagePyramid

What are the differences between the Geoserver stores ImageMosaic, WorldImage, and ImagePyramid? I noticed that ImagePyramid is no longer included in the latest Geoserver. Is one of these more ...
9
votes
2answers
234 views

Are “gridded data” and “raster data” the same thing?

I'm working on a project that will involve both "raster data" (aerial imagery and DEMs) and "gridded data" which I understand to be other kinds of numerical measurements taken at grid points, and ...
9
votes
3answers
3k views

How to get hypsometric curve from elevation data?

I have my elevation data but cannot make an attribute table from it and so I'm not sure how to get any of that data into a table in order to make a hypsometric curve.
9
votes
5answers
665 views

Raster diff: how to check if images have identical values?

Is there a means to check to see if any 2 given raster layers have identical content? We have a problem on our corporate shared storage volume: it's now so big that it takes over 3 days to conduct a ...
9
votes
3answers
884 views

Creating raster by randomly choosing cell value from multiple overlapping rasters?

I'm using ArcGIS Desktop 10 with its Spatial Analyst extension. How do I combine multiple rasters into one, always choosing randomly from the values of overlapping cells? I have an image that may ...
9
votes
4answers
7k views

How can I export raster data from Arcmap to an R *.grd file format?

I want to work with raster data in R, but I have problems exporting my raster dataset from ArcMap. By "export data"/GRID I get a file with an *.aux extension. What I need is *.grd. Can anybody help ...
9
votes
3answers
2k views

Make the nodata area of a resampled orthophoto overview white?

I have a set of orthophoto images which are each fully filled with data, but when placed together the coverage does not form a rectangle so there are nodata areas. In MapGuide, I need to be able to ...
9
votes
2answers
816 views

Python, GDAL and building raster attribute tables

I have an integer raster for which I would like to build a raster attribute table using Python and GDAL. I can create a GDAL raster attribute table in Python as follows: >>> rat = ...
9
votes
2answers
272 views

Methodology for Hansen classification of Global Forest Watch?

I have just discovered amazing Hansen´s classification about forest loss and re-growth, available at http://www.globalforestwatch.org/, published in Science, 2013 as: Hansen, M. C., Potapov, P. V, ...
9
votes
1answer
2k views

Converting point data into gridded dataframe for histogram analysis using R?

I am very new at using GIS data and only modestly experienced with R. I've been reading about how to analyze spatial data using the spatial-analyst.net PDF book, so I'm not completely lost, but I ...
9
votes
3answers
1k views

How to combine corner-touching polygons into a single feature

I have created a polygon layer from rasters, and only have a single classification ID (and calculated acres) for each polygon. I would like to "dissolve" all touching features into larger, individual ...
9
votes
2answers
902 views

GeoTIFFs made in GDAL don't have min/max range when loaded in Arc/QGIS

I'm creating some floating point rasters in GDAL with the GeoTIFF driver. When I load the resulting images into QGIS or Arc the default symbology puts the min at -3.40282e38 and max to 3.40282e38, so ...
9
votes
1answer
452 views

How should I manage PostGIS Raster data with different projections?

I have a requirement to store and manage archaeological geophysics data that is collected as a rectangular array of samples -- a raster image. Each raster will usually 20x20 or 30x30 floating-point ...
9
votes
3answers
705 views

How to extract breaklines from a given DEM?

Breaklines are very useful for describing a terrain - for example to build an exact TIN. My question is, how to extract these breaklines from a given DEM? I calculated the curvature (vertical) - ...
9
votes
4answers
284 views

How to determine if point is surrounded using raster processing?

I am attempting to improve a currently extremely cumbersome vector/python process for a natural hazard model. At the moment we have a lengthy script which generates distance/bearing lines from a given ...
9
votes
4answers
2k views

How to modify single Pixel values in QGIS?

I would like to modify a raster (pixel values) in QGIS. The plugins "value tools" can easily provided me with the value of the pixel of my raster. However I know know how to modify it. I have found ...
9
votes
1answer
615 views

How to view a 155 MB raster stored in PostGIS?

I am using the following command to load a raster into a postgres database: raster2pgsql -I -C -e -Y -F -s 3086 -I -C -M myraster.tif myraster -F -t 30x30 | psql -U postgres -d database -h localhost ...
9
votes
4answers
1k views

How to “clean” and vectorize a raster map?

I'd like to make use of a geological map as shown in the example below, which is only a pretty small and downsized part of the map, so editing/redrawing doesn't seem to make sense in the end. My idea ...
9
votes
2answers
764 views

Scriptable alternative to ArcScan (particularly vectorising lines within tolerances)

Sadly the answer to this question suggests that it is not possible to script ArcScan from within a python script. Therefore, does anyone know of any alternatives to ArcScan (either within ArcGIS, or ...
9
votes
1answer
380 views

Efficient raster sampling of billions of polygons (bounding boxes)

How can a raster be computed efficiently (in Python), given a set consisting of billions of bounding boxes (read sequentially from a file), and given that the raster values for each cell should give ...
9
votes
2answers
489 views

What is the correct way to convert a USGS Orthoimage from WorldImage format to tiled GeoTiff for GeoServer?

USGS provides different resolution orthoimage quadrangles in the WorldImage format (*.tif with an accompanying *.tfw world file). Some of these quadrangles can be used by GeoServer (2.1.0 - 2.1.2) by ...
8
votes
7answers
5k views

Extract raster values at points

How can I extract values from a raster by points. I prefer not in Arcgis. I prefer in Qgis or Mapwindow or other open source gis.
8
votes
6answers
16k views

How to get X Y coordinates and cell value of each pixel in a raster using Python?

I’m really new to Python and I would like to know whether there is a quick method to get cell values of a raster pixel by pixel and the coordinates (map X Y coordinate of centre of each pixel) using ...
8
votes
5answers
330 views

How to ensure uniform colour for DEMs?

This is a silly question, but one I've been trying to figure out for a while. I have two Aster DEM tiles, one next to the other. How do I ensure uniform colour for both?
8
votes
5answers
1k views

Finding minimum bounding extent of given pixel value within raster?

I am wondering if there is a way to find the minimum bounding extent for a raster with a particular value. I clipped a raster from a global image and the extent is set as the global extent with a lot ...
8
votes
3answers
4k views

How to find the locations of the highest values in a raster?

Using ArcGIS 10, I have a raster where I would like to find the pixel with the maximum value in the raster and return its location (center of pixel) in decimal degrees. I would like to iterate through ...
8
votes
2answers
2k views

How do I layer stack rasters in QGIS?

I have some multi-spectral data, where the individual bands are separate Raster files. How do I layer stack these rasters so that I get one raster with 3 bands instead of 3 different rasters? This ...
8
votes
2answers
530 views

How to add different sized rasters in GDAL so the result is only in the intersected region

I'm writing a Python method that adds two rasters and generates a single raster output. For reasons beyond my control, the extents of the input rasters are different, but they do overlap. Is it ...
8
votes
2answers
253 views

More accurate way to calculate area of rasters

In my daily work, I am constantly asked to calculate areas of global raster datasets in geographic projection at 30 arc second resolution. These datasets are normally the result of a Combine operation ...