combining map data of similar dimension to produce new maps using algebraic operations

learn more… | top users | synonyms

4
votes
2answers
143 views

Changing raster calculation values - shorten this script?

i'm doing something quite straightforward regarding raster algebra but am struggling either to find the correct function or make a function work correctly; Basically I have 2 rasters, representing ...
1
vote
2answers
487 views

QGIS from console: raster algebra

I'm newbie to programming in pyqgis. I need to perform a simple "GeoTiff * A" type of calculation(where A is a constant), but I am struggling to do so. I tried ...
0
votes
1answer
1k views

NameError: name 'Con' is not defined

I have 365 rasters with the same extent, and I am trying to sum them together by month. I have the script below, which was working until I added the line out1 = Con(IsNull(raster), 0, raster) which is ...
3
votes
1answer
28 views

Problem with algebra formulas for work with gdal_calc.py

I'm trying to update some values of a raster based on another raster. 1st raster (A): landuse data (values from 10 to 240) 2nd raster (B): mask data (only nodata values and a given value, for ...
0
votes
1answer
137 views

ArcGIS Raster Calculator wont work in main model

I've created a Submodel in ArcGIS, which stretches my rasterdata to new values. The strange thing is, that the model works, if I run it itself - as soon as i implement it in the main model, I get ...
0
votes
1answer
49 views

Can I multiply fractions in map algebra?

I have a couple of rasters that show % of success in different things. I want to calculate the overall probability of success. The numerical values are integers so they represent the percentage. This ...
1
vote
3answers
441 views

Nodata issues in map algebra

A dummy question. I need to change all values in my rasters into ones. I've been able to do this fine in the Rater Calculator just by multiplying the layer by zero and adding one (my_layer*0+1), but ...
1
vote
1answer
61 views

Map Algebra syntax

I'm trying to create a new raster from an existing slope raster where slopes are between 5 and 20 degrees. However, I'm receiving the error "ValueError: The truth value of a raster is ambiguous. ...
1
vote
0answers
41 views

Calculating percentage of cells under given site using ArcGIS for Desktop? [duplicate]

I am using ArcGIS 10.3 for Desktop. I have a file with three hundred of points (shapefile), corresponding with archaeological sites. For each site (point) I have the height, and want to calculate ...
5
votes
2answers
3k views

For Loop - Map Algebra (ArcGIS) using Python

ORIGINAL QUESTION Im relatively new to python and these fourms but Im am having trouble utilizing for-loops for spatial analyses within ArcGIS. It does not appear to recognize the variable, ...
2
votes
2answers
81 views

PostGIS raster MapAlgebra callback optimization

I'm trying to speed up an ST_MapAlgebra callback function that essentially takes a raster with four bands, checks each band for values less than a given number, and where the check is true, return the ...
12
votes
1answer
303 views

Frequency of overlapping isochrones with map algebra in PostGIS

I have a table of polygons representing travel time isochrones on particular days. For each origin point, there are five isochrone geometries (stored in separate rows). For each origin point, I want ...
6
votes
3answers
563 views

Performing raster map algebra in open source web GIS?

Map algebra for raster is easy with ArcGIS Raster Calculator (or QGIS). I had thought this would be easy with web GIS servers such as GeoServer and MapServer plus OpenLayers etc. as well. What I want ...
1
vote
1answer
43 views

Changing raster value based on point feature

I am working on some landslide modeling. Within ArcMap I have a raster layer, representing the area of the county I am working in, and a point layer that shows landslide locations. I reclassified ...
1
vote
1answer
626 views

Using Raster Calculator expression in Python

I've posted about this problem before but I've got a few new questions. So here's the problem. I have two rasters containing bathymetry data: one containing the bathymetry of the sea floor, we'll ...
1
vote
1answer
56 views

Con statement not working in Raster Calculator of ArcGIS for Desktop?

I'm trying to combine two rasters in 10.3 using a con statement with the raster calculator in map algebra. (ie - con (raster 1 ==1, raster 1, raster 2)) where raster 1 is a smaller extent and does ...
2
votes
1answer
210 views

Map algebra from files in multiple directories

I'm a rookie in python programming. I want to do map algebra for multiple rasters but these rasters are in different folders. Here's the code ive created so far: #Importing modules import arcpy from ...
2
votes
1answer
42 views

How to weigh the size of a spatial event with respect to the number of its hotspots?

This is probably a map algebra problem. I have a 2D map with values ranging 0.0 to 1.0 spread all over it. The algorithm that produces them is an application of the fractional Brownian motion. The map ...
0
votes
0answers
23 views

Map algebra with raster [duplicate]

I am doing some map algebra with a raster however the result is not correct. I think it is because the cellsize (x,y) is in decimal degree (0,00027777778, 0,00027777778) and the z values are in meter. ...
-1
votes
1answer
70 views

Adding Rasters with multiple Field Values using Map Algebra in ArcGIS for Desktop? [duplicate]

Embarrassed that I haven't looked at this for some time (can't quite get the syntax down). (ArcGIS 10.2, Win7x64). Within Map Algebra (Raster Calculator) Trying to add the cell values of 3 rasters, ...
2
votes
1answer
766 views

R giving error when attempting map algebra

I'm trying to manipulate some downloaded rasters using R, and keep getting the same error when I try to do any type of map algebra. The example below uses a raster file downloaded from ...
8
votes
2answers
4k views

Using individual bands from multi-band image in ArcGIS raster calculator?

I'm trying to use individual bands from a multi-band TM image (IMG file) within the raster calculator in ArcGIS 10 (SP2/3). I can add individual bands to my MXD, but raster calculator will only use ...
1
vote
2answers
91 views

Alternative to nested Con statements in python

I have a nested con statement that I'd like to put into a script, as it's hard to follow in its current form. The raster it produces is incorrect, so it'd be easier to debug in a script. The ...
0
votes
1answer
105 views

Calculating the difference between two grids in percent with QGIS

I'm using the QGIS GDAL Raster Calculator to get the difference in percent between two grids. Grid A contains an integer, and grid B contains an integer which may be the same as or larger than the ...
2
votes
0answers
87 views

Calculating Aerosol optical thickness from Landsat 7

EDIT I have added question down I am reading one paper and I found method to retrive AOT from landsat 7 image This should be simplifield equation, because they assumed observation zenit as 0. Than ...
5
votes
2answers
186 views

Map algebra in QGIS

How do you calculate standard deviation for each raster if you have set of rasters (100)? I was searching in raster calculator, but I didnt find it. Also google redirect me to using R or python. But I ...
1
vote
2answers
59 views

How can I stop my map algebra being limited to the size of my smallest raster?

I have five input rasters and I'm using map algebra to find areas where all of them have a value over 75, areas where all of them are over 25, and areas where one or more are under 25. I can do this ...
0
votes
1answer
147 views

How to speed up raster/map algebra in PostGIS?

I am testing map algebra with PostGIS (2.1 with PostgreSQL 9.3) on two images, and found it to be rather slow. For example, dividing two images of 7200 * 3600 pixels took 110 ~ 120s on a desktop with ...
1
vote
1answer
271 views

Arcpy.Delete_management() not deleting folder?

I've written a script in Python that takes in two ascii files that represent surfaces, then convert them to Rasters, and finally Negate their values (i.e. multiply all cells by -1). In the script, I ...
4
votes
1answer
2k views

Map Algebra expression in a model fails on open workspace

I have model which runs 2 processes that create rasters and a third process that (attempts to) add the two of them using Map Algebra. The expression is very simple, just "%Raster_A%" + "%Raster_B%". ...
1
vote
0answers
69 views

How to correctly input complex multiband functions in Raster Calculator?

I'm trying to compute a MTVI2 (Modified Triangular Vegetation Index - Improved) and EVI (Enhanced Vegetation Index) from Landsat 8 imagery. My formulas are as follow (I'm using ModelBuilder): MTVI2= ...
5
votes
1answer
535 views

Error using con in Map Algebra in Python script

I am experiencing a painful problem in the middle of a Python script I am writing. When I get to a Map Algebra function, I have tried two different options and both return a different error. At first ...
1
vote
2answers
332 views

Access individual bands and use them in map algebra

I have a list of multiband rasters, consisting of: Red=Band 1, NIR=Band 2, SWIR1=Band 3 and SWIR2= Band 4. What I want is to access only the red and NIR bands and save it to new file. Then, I'll be ...
1
vote
0answers
39 views

Automating Raster Calculator [duplicate]

I am seeing mixed answers on other posts so I thought I would ask here. Is there a way to automate the raster calculator? I am trying to create a series of inundation maps at fine increments. It ...
0
votes
0answers
90 views

multiply a rang of values in a floating raster to a constant

I have a floating raster layer. I want to multiply a rang of its values to a constant. So, I used the "Con" tool in map algebra with the following script. Con(("Extract_rast1" >= 0) & ...
1
vote
1answer
2k views

Piping output of Iterate Raster tool into Raster Calculator tool using ModelBuilder?

ArcGIS v10.1 ModelBuilder Using ModelBuilder and map algebra, I want to compute a ten year average grid from a stack of 10. The problem is that only 1 grid (%raster%), the last one through the ...
2
votes
0answers
232 views

Plotting Vector Arrows from U and V Component with ArcGIS for Desktop?

I have two raster layers representing the U and V components of a vector wind. I could generate the magnitude raster by applying Square root in Map Algebra in ArcGIS. I also generated a raster layer ...
2
votes
0answers
212 views

Arcpy script for summing pixel values in loop can't handle large number of rasters

I want to sum the pixel values of large numbers of rasters of different extents (ArcMap v10.2). 16GB RAM. I have the following code (with help from this site) which solves the different extent ...
12
votes
4answers
353 views

What are some Domain Specific Languages (DSLs) for GIS?

I'm curious what sort of domain specific languages (DSLs) are used in GIS. The only one I can think of is the DOCELL language that was provided in ArcInfo Workstation GRID, but is no longer ...
2
votes
0answers
182 views

How to compute the ls factor equation for USLE model in qgis grass r.mapcalc

To compute the LS factor for the USLE model, the following equation easily works in ArcGIS. How can I perform the same in QGIS under GRASS r.mapcalc? LSfator = Power (("flowaccNW") * ...
4
votes
1answer
515 views

How to do loops on raster cells with python console in QGIS?

I'm new here and beginner with python. I'm looking for a solution to my problem for some hours but now you are my last chance... I'm using the python (2.7.5) console in QGIS (2.4.0). Is there a ...
9
votes
3answers
1k 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 ...
8
votes
1answer
1k views

Partially updating raster (DEM)?

How to update an .img DEM in ArcView 9.3? I do have DEM which need to be updated only if few parts. The area to be updated has been produced in separate DEM with Spline Interpolation with Barriers. ...
0
votes
1answer
268 views

How to calculate cumulative sum conditionally across raster?

I don't know if this is possible without coding, but I need to: Create a raster file that, for each value in the input raster, sums all pixel values across the entire extent that are less than or ...
0
votes
1answer
549 views

How to compare two rasters and reclass the greater value(python script)?

I am working on an assignment that requires us to reclass a raster(via python) from two rasters. We need to compare the power received at one cell(on raster1) and compare it to the same cell on raster ...
1
vote
1answer
1k views

Raster Multiplication - Can't Multiply Sequence by Non-Int of Type “Float”

I am trying to multiply a list of GRID raster datasets using a different value for each raster. The value is in the raster name. Some examples are [wse_1979_5,wse_1980_0,wse_1980_5], where 1979_5, ...
1
vote
1answer
1k views

Map Algebra - Grass GIS (r.mapcalc) merge two raster with NoData value

I have a slope raster and a road raster. The black pixels in the road raster are NoData values and the pink pixels have a value of 1. They don't cover exactly the same area. I want to create a new ...
0
votes
1answer
146 views

How to iterate RGB raster layers in ModelBuilder?

I have an RGB raster that I'm calculating color-based indices from. When I do the calculations manually, I can select a certain band of the RGB raster to calculate. I need to do the same operations on ...
1
vote
0answers
36 views

Map algebra from QGIS python console? [duplicate]

What command do I use to do map algebra on rasters from the QGIS python console?
1
vote
0answers
194 views

How can I ignore NoData values in map algebra? [duplicate]

I need a way to exclude NoData from calculations. I don't want any NoData cells from any of the rasters to be used in pick or con statements.