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

learn more… | top users | synonyms

2
votes
2answers
631 views

Reclassify rasters using GDAL and Python

Using @SS_Rebelious very useful answer, I am trying to reclassify rasters using GDAL and Python. I have a grayscale raster (with values ranging from 0 to 255) and I want to reclassify it to a binary ...
2
votes
3answers
10k views

How to calculate Raster Area per Polygon?

In ArcMap, I have a raster file that contains 6 classes of landcover and a shapefile that is the census data (OA boundaries and their population). I want the area of each class in each OA? For ...
8
votes
1answer
834 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
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 ...
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?
6
votes
1answer
1k views

Counting NoData values in raster using ArcGIS Desktop?

I have created a raster in ArcGIS Desktop 10 with a large number of holes (NoData) in it. I want to fill the holes using something like the the following: rOut = ...
5
votes
3answers
2k views

Convert an arbitrary meta-data-free map image into QGIS project

My background: experienced programmer, amateur historian, complete beginner to GIS (so I'm using QGIS). I did learn to use the QGIS georeferencer tool though. The input: suppose I have an image like ...
4
votes
2answers
322 views

What is the source of horizontal and vertical striping in USGS DEMs?

When processing 30m and 10m DEM data downloaded through the National Map Viewer from the National Elevation Dataset, we noticed horizontal and vertical striping not only in the produced results, but ...
4
votes
1answer
1k views

Getting a constant background color for transparency

I am creating a series of tif files at 1:10k scale over a vast area (760 tiles) using arcpy.mapping.ExportToTIFF and then stitch them together to create a single ecw. The output looks ok but when I ...
4
votes
2answers
2k views

How to find the area of raster values within polygon zones?

I'm using ArcGIS 9.3 and I just can't seem to get this right - I've tried tabulate areas, zonal stats, converting to polygons, etc. I'm missing something somewhere. I have a USGS impervious surface ...
3
votes
2answers
1k views

In ArcMap 10.2, exporting raster as .tif with “use renderer” and “force RGB” still produces inaccurate colors

I have a DEM and my aim is to change its color ramp using ArcMap, then export this colored raster as a geotiff (8 bit). I checked "use renderer" and "force RGB" options but the results are not good ...
3
votes
2answers
1k views

How to calculate the area of rasters in QGIS?

I have been working with some rasters and would like to know how to calculate the area occupied by one element of a raster in Qgis.
2
votes
0answers
246 views

Problem loading output layer QGIS

Every time I run operations with temporary file output from the graphical modeler, I get the following error message: Oooops! The following output layers could not be open CLIP: ...
2
votes
1answer
628 views

Simple band arithmetic on GDAL virtual rasters (VRT)

I want to calculate the difference of two bands. I have the individual bands stored somewhere, but I was thinking that maybe there's a way of producing a VRT file that is the result of this ...
1
vote
1answer
53 views

how to produce a RASTER representing X coordinates [duplicate]

In a journal article I have read that two rasters, representing X and Y coordinates respectively (see attached figure), have been entered as predictors in a GIS-based Logistic Regression for a given ...
1
vote
0answers
143 views

Correlation between point feature (disease) and elevation

I have a problem! How can I calculate the correlation between point feature and Dem raster? Hypothesis: disease prevalence in high altitude is more than low altitude. I must do this in ArcGIS 10 and ...
1
vote
1answer
533 views

How to eliminate small features in a grid?

(Update) My final goal is to take some logged point data (1 point about every 20 ft or so there are over 15,000 points) and interpolate that point data. Then I would like to convert the interpolated ...
0
votes
3answers
71 views

Deleting polygon defined area from raster layer in ArcGIS for Desktop?

Which tool do I use in ArcGIS 10.1 for Desktop to exclude areas covered by vector polygon layers from a raster layer? I want to keep portions of the raster layers that exclude the portions covered ...
24
votes
3answers
6k views

What is Lanczos resampling useful for in a spatial context?

GDAL includes a resampling method beyond the normal mix of nearest neighbor, bilinear, cubic and splines: "Lanczos windowed sinc resampling". I understand that its a convolution filter, but unlike ...
27
votes
7answers
2k views

Work flow for georeferencing imagery using open source tools?

Does anyone have a recommendation for a work flow for georeferencing high resolution aerial photography using open source tools? Comparison of user friendliness and robustness between Quantum GIS ...
14
votes
6answers
1k views

Creating simple image from vector data

Can anyone suggest a quick (and preferably open source) way to create an image from a shapefile. For example, I have a shapefile containing parcel boundaries, I want to create a png file depicting ...
20
votes
5answers
4k views

What is Change Detection and how I can perform such analysis with open source tools?

From wikipedia page: Change detection for GIS (geographical information systems) is a process that measures how the attributes of a particular area have changed between two or more time periods. ...
13
votes
2answers
3k views

How to smooth/interpolate a raster in python using GDAL?

This is my first post so please feel free to critique where necessary. Background: I am developing in Python and using GDAL from OSGEO to manipulate and interact with rasters and shapefiles. ...
12
votes
1answer
2k views

Random sampling of raster using R

Is there a straightforward way of randomly sampling a raster so that the output of the process is a raster? I'm using an example that I found on the r-sig-geo list and I have also tried the ...
7
votes
2answers
984 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 ...
10
votes
3answers
6k 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 ...
17
votes
2answers
2k views

What resampling technique should be used when projecting aerial photos?

I'm doing some time-intensive projecting of aerial photos, and I'm curious - what resampling technique is best to use on aerial photos? In ArcMap, my options are NEAREST, BILINEAR, CUBIC, and ...
13
votes
2answers
6k views

Differences between Raster Dataset, Raster Catalog, and Raster Mosaic

I have about 40 scanned, colorful Geology maps that are GeoTiffs. I want to bring them into an Esri FileGeodatabase to lower their size, combine them as one and display them with a mostly transparent ...
11
votes
2answers
15k views

How can I join several raster files using QGIS?

I'm new to QGIS. How can I combine multiple raster layers into one layer? The raster images are of different areas with very little overlap. The goal is to bring numerous town plats into one project. ...
11
votes
1answer
884 views

Method to shade or overlay a raster map to reflect time of day and ambient light

I have a daytime raster of my study area. I have 20+ counters in the study area that count "stuff" with a timestamp. Using PyEphem I determined the type of ambient light for each collected data ...
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') # ...
3
votes
1answer
323 views

Random-Forest Classification of 10cm Imagery for species-distribution in R (no point-shapes)

I'm actually trying to classify 10cm-resolution images with the RF-package. I do this to separate some species in the forest, seen on my images. While I was searching for an answer, I found the ...
13
votes
5answers
724 views

Dikes extraction from DEM

I do have DEM of 1x1 meter and original LiDAR point cloud in .las, that DEM is made of. I need to extract river dikes (highest points of the dikes) to vector feature (point, polyline). Any idea for ...
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
6answers
7k 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): ...
6
votes
2answers
2k views

creating a raster of the residuals of a regression between two rasters

I am trying to create a raster of the residuals of a regression between two rasters. i.e. I would like to carry out a regression of one raster against another of the same extent and plot the residual ...
6
votes
1answer
4k views

How to import a raster into PostGIS?

I'm trying to follow the instructions here to load a raster into a PostGIS database: python raster2pgsql.py -s 4269 -I -r *.tif -F myschema.demelevation -o elev.sql I understand that I replace * ...
6
votes
2answers
1k views

Creating shapefile of satellite image footprints?

I'm looking for a tool or a workflow to compute satellite image footprints, for a large number of images across directories, and produce a shapefile. Ideal work flow would be: fire up tool point to ...
4
votes
3answers
2k views

How can I specify the CRS of a raster layer in pyqgis?

I manipulating netcdf raster data from the QGIS python console. The metadata of my netcdf file does not provide the Reference system. However I know it is 'EPSG:4326' I am writing a python script to ...
3
votes
2answers
4k views

How to convert .geotiff --> .tiff?

How to convert .geotiff files to .tiff files (ordinary plain images)? Basically, I have cropped the image to a rectangle of interest and knowing its left-lower location and having chosen the pixel ...
12
votes
3answers
23k views

Extract Raster Value into Polygon Attribute

I have Raster Data for a State which has Elevation Data in its pixels. I also have Vector layer which has the 300 divisions of State. what I need is, I want to get the average elevation value of ...
9
votes
2answers
3k views

How to reclassify a very large land cover dataset?

Consider the NLCD2001 Land Cover dataset for Alaska (download link). I need to reclassify this dataset so that only pixels of value 41, 42, and 43 are preserved; all other pixel values should become ...
7
votes
1answer
322 views

How to create points around edge of viewshed results?

Is there a tool or module, or multiple tools, in ArcGIS that will allow you to automate the creation of points around the outer edges of a viewshed? What I am trying to do is gather data about the ...
6
votes
1answer
943 views

Python equivalent of gdalbuildvrt

Is there a way to perform the same task as the gdalbuildvrt utility using GDAL Python bindings? So far I have not found any way to do this other than creating a vrt of a single dataset and manually ...
6
votes
3answers
3k views

Raster: How to get elevation at lat/long using python?

I was wondering if anyone has some experience in getting elevation data from a raster without using ArcGIS, but rather get the information as a python list or dict? I get my XY data as a list of ...
5
votes
2answers
2k views

How to Merge Raster Datasets of Different Resolutions?

I have two DEMs in .tif format and I am looking to merge the two. The main problem with merging these are that they are different resolutions (one is 10m and one is ~20m). Also, another problem with ...
5
votes
1answer
2k views

Loading a raster using QGIS into a PostGIS2.0 enabled database

Following on from my attempt to install a postgresql 9.1 db with postgis 2.0 on windows 7, I am trying to figure out how to load a raster. I have successfully loaded a shapefile, and am trying to do ...
5
votes
1answer
351 views

Unknown Coordinate System on old drawing

I work with many old civil engineering documents, and occasionally I'll need to georeference one which has coordinates pairs listed at many given points but does not list which coordinate system these ...
4
votes
1answer
1k views

How to store raster data in spatialite?

I want to cache my raster data(large geotiff files) on local database. I thought to use spatialite. I found librasterlite which provide functionality to do that. Can anybody point to some tutorial on ...
4
votes
1answer
1k views

Why does QGIS's Clipper tell me “Computed -srcwin falls outside raster size” and how do I address it?

I'm trying to clip a heightmap to a small area that I'm interested in. When I attempt to use Raster > Extraction > Clipper I get an error dialog with the message ""Computed -srcwin falls outside ...