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

learn more… | top users | synonyms

8
votes
2answers
288 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 ...
8
votes
3answers
2k views

Determining statistical relationships between rasters using R vs ArcGIS Desktop?

I am trying to analyze how sets of rasters relate to each other using some statistical techniques. As, I don't have much experience using the spatial statistics tools in ArcGIS I was exporting my ...
8
votes
2answers
1k views

Resample binary raster to give *proportion* within new cell window

I would like to go from a binary raster of forest/non-forest at 30m resolution to a raster at 240m resolution with the value of each cell being the proportion of that cell forested - i.e. the data ...
8
votes
3answers
3k views

Using polygons to apply unique value to all raster cells within those polygons?

In ArcGIS 9.3, I have a raster file that I would like to "mask" with polygons, giving a single value to all the raster cells that overlap the polygons. How can I do this? Here's a screen shot, ...
8
votes
1answer
6k views

Converting a polygon into a raster using R

I'm trying to convert a shapefile into a raster within R. My approach is to read in the raster as follows: library(rgdal) # Loads SP package by default demo <- readOGR('F:/data/', 'shapefile') # ...
8
votes
2answers
766 views

Creating randomly shaped clumps of cells in a raster from seeds of 1 cell/pixel?

As my title says, I wish to "grow" clumps of cells from seeds within a raster. My base raster is full of 1's and 0's, 1's indicate land and 0's sea/NA areas. From the 1's I wish to select 60 random ...
8
votes
2answers
598 views

How to create a solar intensity map for a garden plot?

I'm interested in creating a raster surface that depicts the average solar intensity per unit area for a small garden plot. I imagine the inputs into the model would be the following: Plot ...
8
votes
1answer
4k views

How to reproject a raster file in QGIS with datum transformation?

I want to reproject a tiff raster file from projection A (ETRS89/UTM Zone 32) to projection B (DHDN/Gauss-Kruger zone 2). The two projections have a different datum so I have to use a datum ...
8
votes
6answers
1k views

seeking “Available Sky” implementations

I'm looking for implementations of an algorithm called "Available Sky". ArcGIS (Spatial Analyst or GRID) is preferred, but solutions in GDAL, SAGA GIS or others are perfectly acceptable. The ...
8
votes
2answers
1k views

How to keep spatial reference using arcpy.RasterToNumPyArray?

I am using ArcGIS 10.1, and want to create a new raster based on two preexisting rasters. The RasterToNumPyArray has a good example which I want to adapt. import arcpy import numpy myArray = ...
8
votes
1answer
516 views

Reduce extent of raster processing in ArcPy or ArcGIS Desktop?

I need to reduce the processing extent of my raster by some sort of automatic procedure in python. Currently the borders of the dataset I've been using has anomalies (due to the classification of the ...
8
votes
1answer
2k views

Is there a fast way to copy large raster datasets between File Geodatabases?

Why does it take so long to copy a large raster dataset from one file geodatabase to another? I can copy the entire database in windows explorer in minutes but copying a single raster in ArcCatalog ...
8
votes
2answers
1k views

Extract the area of a DEM above where 60% of the region lies?

I have a DEM of a mountainous region. I'd like to extract the areas of the DEM above which 60% of the total DEM area lies. So, essentially I'd like to extract the areas representing upper portions of ...
8
votes
3answers
875 views

Best way to manage large number of geoTiff rasters

I have around 2-3000 1-bit geoTiffs of size 11000x9400. Color 0 should be transparent, and QGis gets it wrong in the first place, so I used a .vrt file to set all the parameters in one single place. ...
8
votes
1answer
939 views

How to apply a Majority filter in GRASS GIS?

I was wondering if someone could help me write a majority filter for GRASS GIS. I have this map (see below) which I need to run a filter over to make the red areas a bit more dominant. I know I can ...
8
votes
2answers
903 views

How to convert a vector boundary to a raster map?

I need to output a PNG file of a raster map overlayed with boundaries present in a vector map. The problem is, the r.out.png script only exports a raster map. My boundaries are in a vector format. ...
8
votes
4answers
1k views

Looking for a fast, open-source raster cost-distance function to use in code

I've been scouring the internet for a fast, open-source cost-distance function that I can embed inside my code. The functionality should basically work like the ArcGIS implementation where I can ...
8
votes
1answer
806 views

Calculating distance from points to different class in raster

I am looking for some advice on how to calculate the distance from a point to a different class of a raster. Specifically, I have observations (points) of animal use overlaid on land cover data. For ...
8
votes
3answers
2k views

Poor image quality of raster pyramid files in ArcGIS Desktop (no anti-aliasing)?

There are several GeoTIFF files (size ~7000x7000 px) which I try to open in ArcGIS Desktop 10.0. I have already created raster pyramids for them (.ovr files) using geoprocessing tool but their quality ...
8
votes
1answer
881 views

How to add a band to an existing GeoTiff using GDAL?

This seems like it should be straightforward but I can't for the life of me figure it out. Given an existing GeoTiff containing a single band, how can you add another band to that GeoTiff? Is the ...
8
votes
1answer
2k views

How to make a new raster using con() script in Python?

I'm new to Python and have been trying to code a raster calculator style script. Basically, I want to convert a raster that has values 0 to -10,000 into a raster that is populated by 1s between ...
8
votes
1answer
3k views

GDAL RasterizeLayer Doesn't Burn All Polygons to Raster

I am trying to burn a shapefile to a raster using GDAL's RasterizeLayer. I pre-create an area of interest raster from a different shapefile, given a specific pixel size. This AOI then serves as a base ...
8
votes
4answers
306 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 ...
8
votes
1answer
2k views

Smoothing raster features into polygon

I have raster data representing urbanized areas, which I'd like to trace and convert to a polygon. Using ArcMap's Raster to Polygon tool, I'm able to get something close to that. Unsimplified: ...
8
votes
3answers
2k views

calculate area of vegetations in a raster file

How do I calculate area of several vegetation types in a raster (.img file type) image. They are color coded. Am using QGIS (1.8.0) and ERDAS 9.2.
8
votes
4answers
405 views

Identify corrupt shapefiles or rasters via Python or ArcObjects

I'm interested in a way to identify rasters (e.g., ESRI GRID format) or shapefiles that are corrupt that would prevent them from being used in an analysis tool. I have a large number in structured ...
8
votes
3answers
383 views

R, rasterVis and RGL: how to prevent reflective surfaces?

I have been experimenting with digital terrain model files in ESRI ASCII grid raster format. Despite zero experience with this type of data I found it very easy to load into R, as per code below. The ...
8
votes
1answer
553 views

How to grid unevenly sampled categorical data?

I am looking for a method for gridding categorical data. I have extracted from nautical charts and field sheets a collection of points specifying the nature of the surface of the bottom of the ocean. ...
8
votes
1answer
359 views

Where is this raster offset being introduced in the QGIS processing toolbox?

I'm using the (awesome!) processing toolbox in QGIS 2.0 to reclassify some NoData values in a raster from -9999 to 0, but at some point in the process there's an offset of about 0.2 degrees being ...
7
votes
2answers
2k views

How to build an artificial 'dam' across a river?

I have a DEM. Created all files related to watershed including accumulation, direction, streams etc... I'm stuck creating a dam. I request you to please let me know how to build an artificial 'dam' ...
7
votes
3answers
3k views

How to change null values raster values to 0 in QGIS?

I need help... I'm having problems with QGIS. I have some rasters that appear to have null values, even though I reclassified them to obtain a boolean raster :( How can I change the NULL values to 0?
7
votes
2answers
6k views

Find the intersection of two rasters

I have two raster images, Boolean classified (1 or 2). They overlap in given areas. I need a raster that is only the overlapping sections. Using Arcinfo 10. Answer can only use raster methods, no ...
7
votes
3answers
3k views

Simple Question: Raster from ArcGIS to Geoserver

we have a couple raster files which could be open in ArcGIS. They come with the extension of .aux and .rrd and they are the GRID file format from ESRI. We want to put them on Geoserver to overlay them ...
7
votes
3answers
6k views

How to add Raster legends in QGIS map composer?

Is there a plugin or method to add legends for raster layers (with a custom colourmap) as an item in the print composer?
7
votes
2answers
1k views

Creating a multispectral image from scratch

I want to make a multispectral image from cero to do some tests on it. Something really simple like 5 completely uniform bands with salt and pepper noise on them or a square of different values at the ...
7
votes
3answers
1k views

Merging hundreds of Aerial Photos into one big file?

I have tried to merge my 405 Aerial Photos into one. I have FME and ArcInfo 10.0 at my disposal on Windows 7 64bit. I have tried the following: create a mosaic; export all images from a ...
7
votes
3answers
2k views

GDAL and Python: How to get coordinates for all cells having a specific value?

I've got an Arc/Info Binary Grid---specifically, an ArcGIS flow accumulation raster---and I'd like to identify all cells having a specific value (or in a range of values). Ultimately, I'd like a ...
7
votes
1answer
4k views

Best Way to Handle Many Raster Files in QGIS?

I have a project with a large number of raster files, but it is becoming very slow. I wonder which is the best way to handle all this. I'm simply georeferencing these files creating .wld file (geoTiff ...
7
votes
2answers
1k views

Extracting roads out of Landsat raster images

For a study into the expansion of road networks in the rainforest I am trying to extract roads out of Landsat images. We already have sharp and cloud-free composites on which the roads are clearly ...
7
votes
3answers
4k views

Raster mosaic not displaying raster values ArcGIS 10

I have eight DEM tiles that I downloaded from the USGS National Elevation Dataset. The individual rasters are displaying correctly, but when I try to combine them in a Raster Mosaic I get arbitrary ...
7
votes
3answers
1k views

Creating raster or vector raster from DEM point file?

i have a point shape of 75.000 elevation points in a regular 5 m grid. We now want to create a polygon grid with the elevation points as "centroids". Unfortunately, all attempts to do so fail: I ...
7
votes
4answers
4k views

How to install “1-band color raster color table”?

I would like to use the "1-band raster colour table" extension. I'm using QGIS Wroclaw 1.7.4. In the plugin installer I can't find this extension, even I quote for third party extensions thing. Does ...
7
votes
2answers
252 views

Finding specific pixel values ​​in ArcGIS using Python

I have a raster file containing some integer values. I need to find all pixels whose values are not defined in list like this one: l = [1,2,4,8,16,32,64,128]. Is there any easy way to solve that ...
7
votes
2answers
5k views

How to extract values from rasters at location of points in R?

My question is about extracting values from rasters at the location of points. With the function extract this is very easy, and the function gives me a dataframe with the values of all the variables ...
7
votes
2answers
6k views

How to clip a vector feature to the extent of an existing raster dataset using arcpy

I have a clipped raster dataset with polylines (from the contour tool) overlayed (see attached). I am trying to clip the polylines to the extent of the raster dataset, which seems simple enough, but ...
7
votes
3answers
3k views

Creating a raster surface of overlaps from vector polygon shapefiles?

Basically what I'm trying to do (in ArcGIS 10) is to compare three vector files in order to create a 'comparison' raster layer. I tried using Data Comparison (under Data Management tools) but it ...
7
votes
4answers
3k views

Garmin etrex raster maps support

I currently use an excellent set of printed topographic route maps for mountain biking. I would like to scan the maps, georeference them, and upload the maps to a handheld GPS device. Does the etrex ...
7
votes
2answers
2k views

QGIS and Raster Mapping in 'Seamless' layer

I would like to use Ordnance Survey base mapping in QGIS. I would like to set up something like a seamless table in MapInfo. I think that the equivalent in QGIS is a Virtual Raster Table. The files ...
7
votes
2answers
3k views

Modify raster layer display and symbology properties programmatically with arcpy or ArcObjects

I have over 1500 NITF files that I need to be able to modify the layer display and symbology of programmatically from within ArcMap. I need to set the properties outlined in red below: I've looked ...
7
votes
3answers
803 views

R - How to sampleRandom() from large raster without NA values

I have a number of vary large rasters which need to be randomly sampled with the return value being a matrix of x, y, and value. The raster package sampleRandom(raster,n, na.rm=TRUE, xy=TRUE) will do ...