Tagged Questions
2
votes
1answer
42 views
Tabulate Area crashes Python and ArcGIS
I have a script that runs Tabulate Area on a single feature class multiple times with different rasters. In total Tabulate Area runs 5 times. 2 of the input rasters are provided as arguments to the ...
3
votes
1answer
38 views
Is it possible to batch clip rasters in Quantum GIS?
I have been able to batch clip rasters in one directory using python in ArcGIS. However, I am not sure how to do the same in Quantum using the python console. Is this possible? If there are different ...
6
votes
3answers
108 views
How to efficiently read and reclassify many rasters in R
I've been tasked to create a suitability analysis of wave conditions in the Gulf of Mexico. I have 2 thousand or so raster files that are about 8 MB each (2438 columns, 1749 rows, 1km cell size). ...
2
votes
1answer
87 views
How to calculate volume from a raster file?
I have a raster file containing the depth of water. I need to find the volume of water using that raster file. I need to find the area of each cell and multiply with the depth. How can I do it in ...
2
votes
2answers
46 views
Raster format to read with python
I currently have a DEM raster stored in an ESRI geodatabase. I've been asked to create a purely python (e.g., no arcpy functions) program that works with the values. Therefore, I need to convert the ...
0
votes
1answer
34 views
Maintain 8-bit pixel depth in Raster Calculator
Would anyone know how to maintain 8-bit pixels when computing raster calculations in ArcMap?
0
votes
1answer
49 views
How to find locations for hydro power using raster analysis?
I have a stream network and DEM file both in raster format. I need to extract values from DEM to stream network. For the analysis in stream network after that, I need to find out the places for hydro ...
0
votes
2answers
75 views
Raster Calculator question for Multiple Rasters
I am having trouble for writing script in Python (ArcGIS 10.1). The task is generally divide each raster by 100. I have many of these raster files....
So what I've done is the following:
from ...
5
votes
1answer
128 views
ModelBuilder: How to extract data by date, export to separate tables and interpolate tables to rasters
I've been using ArcGIS for a couple years, but when it comes to model builder I've only worked with very simple models. I'm hoping someone that is much smarter than myself can help me in building a ...
2
votes
1answer
168 views
GDAL_translate: converting ESRI GRID to Geotiff in batch
The Challenge:
To batch convert 1000 rasters in ESRI Grid binary to Geotiff using either using a batch file or Python 2.7 without ArcGIS/Arcpy module.
The Tools:
Win7Pro, GDAL 1.9.2 (via OSGeo4W), ...
4
votes
4answers
110 views
How to display a “Save As” dialog for raster files?
I'm writing a QGIS plugin in python. I need to show a "Save as" dialog window in order to allow the user to select the path of the output raster file. Do I need to create a QFileDialog and manually ...
1
vote
0answers
53 views
Creating Raster with values based on function output
I have a function that calculates a price based on the input slope and distance. I want to write the price to a raster as the rastervalue. How do I do that? OpenSource and ArcMap solutions would work.
...
0
votes
0answers
138 views
Figuring out Latitude and Longitude from XY Image [closed]
I'm having difficulty calculating the latitude and longitude values of an image that I am hosting inside of a div tag on my webpage. I can easily figure out the x,y (pixel) coordinates of the image ...
3
votes
0answers
73 views
How do I write XML metadata with Python?
Some Python software I have written generates raster files via GDAL. I would like the software to create an XML metadata file for each raster so that I can include details about how the raster was ...
2
votes
2answers
116 views
Calculate Percent coverage using another raster, with dif resolution and orientation
I have a 0.5 meter raster that defines presence or absence (cells are either 1 or 0). I have another 5 meter raster, which I created. I would like to calculate a percent coverage of the underlying 0.5 ...
3
votes
2answers
138 views
Majority filter (minimum mapping unit) in ArcGIS with larger window size
I want to run the majority filter in ArcGIS with a window size of 4 or 5 i.e. greater than the default 3. Is it possible? Or some other command (or tool) which might help me achieve that?
thanks!
1
vote
0answers
47 views
How to clip Raster using multiple polygons and Python? [duplicate]
Possible Duplicate:
clipping a raster by multiple datasets or polygons
I have polygon layer and DEM and I want to write a program using python to cut the raster into multi-raster based on ...
1
vote
1answer
46 views
cannot print listRasters
I'm trying to run this very simple script to check a problem in another more comlpex one.
import arcpy
arcpy.env.workspace= 'C:\\Users\\Umberto\\Desktop\\Prove_Script\\Karakoram\\temporary_01'
rst= ...
1
vote
1answer
82 views
How to add multiple rasters to a layer file?
I'm trying to make a .lyr file to bring into ArcMap with arcpy.mapping. I can bring a single raster into a .lyr with the script below. I'm unsure how to proceed to get multiple rasters from the same ...
7
votes
2answers
172 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 ...
1
vote
3answers
60 views
Apply rotation before FeatureToRaster_conversion
I am converting polygon shapefiles to rasters. With arcpy I am using the function FeatureToRaster_conversion.
gp.FeatureToRaster_conversion(InFeatures, InField, OutRaster, InCellSize)
This function ...
2
votes
1answer
125 views
Saving rasters in a Python for loop fails only on last iteration
I'm a newcomer to Python scripting and I'm having what's seemingly a very simple problem but one that's stumped me enough to send me posting to these forums for the first time.
I'm using a for loop ...
3
votes
1answer
484 views
ArcGIS 10 - Raster cell values and raster attribute tables - reassign raster cell values from another attribute
It seems like this should be relatively easy to accomplish, but I've been messing with this issue off and on for a couple days and have not successfully figured it out.
Problem: I have a large (size ...
5
votes
1answer
91 views
GDAL/OGR Determine if dataset is raster or vector
I am writing a script using the GDAL/OGR Python bindings and it is possible that the user may provide a raster dataset or a vector dataset as input. If the user provides a vector dataset, the script ...
5
votes
1answer
186 views
How to load raster data with hyperlink path in ArcGIS MXD?
I have a problem/question and hopefully someone of you is able to help me:
I have thousands of georeferenced raster files which have corresponding polygons (outer borders of rasters) with a hyperlink ...
3
votes
1answer
142 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 ...
7
votes
3answers
475 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 ...
4
votes
5answers
277 views
Python loop Polygon to raster not working
I want to convert a large number of different shp polygons to raster in ArcGIS 10.1, using an existing field column "class" for the new raster classes. I am not very fluent in python but I thought ...
3
votes
2answers
567 views
Extract raster by mask output does not completely overlap the mask
I have a flow accumulation grid and a shapefile containing catchment areas. Now I want to do some calculations on them so I need to extract data from cells in the raster adjacent to certain catchment ...
7
votes
2answers
246 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 = ...
3
votes
1answer
149 views
Problem editing raster and saving with projection and transformation intact
I am trying to make a custom tool to Arcmap using Python. It is supposed to read in a raster DEM, a point shape file and a distance, and output another raster.
The output raster should be ...
8
votes
1answer
321 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 ...
3
votes
1answer
109 views
raster2pgsql: Strange error when creating a table
Hi GIS Stackexchange people,
When I try to load a asc file into my psql setup, I get a strange error message.
I have created the sql input file via
raster2pgsql -s 4236 -I -C -M gis_data.asc -F ...
10
votes
3answers
843 views
Writing numpy array to raster file
I'm new to GIS, and I'm lost with something I feel ought to be relatively simple. I have some code that converts infrared images of Mars into thermal inertia maps, which are then stored as 2D numpy ...
3
votes
1answer
80 views
Probability distribution from list of rasters
I have 100s (exactly 224) of wind distribution rasters. I want to create a probability distribution raster of wind speed using all these rasters, so we can predict the possible wind speed in any ...
10
votes
2answers
161 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?
3
votes
3answers
196 views
How to compare a set of images with a shapefile?
I have around 150 images (rasters), each one is just a map of certain "classes", e.g. pixel value 1=forest.
I also have a shapefile that consists of roughly the same area as the 150 images constitute ...
7
votes
1answer
240 views
Reduce extent of raster processing
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 ...
0
votes
1answer
176 views
zonal histogram fails to create output table when looping
Using python I'm having problems when trying to attempt a FOR LOOP with the ZonalHistogram function of ArcGIS 10.1. It works fine if I use it only for one image but it gives me this error while for ...
4
votes
3answers
380 views
How do I extract a raster's extent in python?
I am looking for an easy way to extract a raster's extent in python.
Currently, I am using gdalinfo and then parsing out the corner coordinates. The problem with this approach is that gdalinfo is ...
6
votes
2answers
738 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 ...
1
vote
2answers
148 views
Cannot get arcpy.MosiacToNewRaster to work with list
Below I have a snippet of a much larger code that I cannot get to work and I am not sure why. It seems that I always have trouble using arcpy functions with rasters. Can anyone help me get this code ...
1
vote
0answers
115 views
Reclassify ESRI raster using variables from a SQL Server table (connecting with pyodbc)
I am trying to reclass a raster file (representing vegetation and land cover types) using values stored in a SQL Server Express table. I have a set of ~400 species records in the table, with each row ...
0
votes
0answers
18 views
Python syntax problem - multiple conditional in raster processing [duplicate]
Possible Duplicate:
Python problem in conditional combining of multiple rasters and produce one output
could anyone suggest the correct syntax for this statement:
outCon1 = Con (((inRaster1 ...
2
votes
2answers
309 views
Python problem in conditional combining of multiple rasters and produce one output
I have two rasters for same spatial extent and need to compare both rasters cell by cell and produce one output raster. Each has 10 land cover classes. Based on those 10 classes the comparison should ...
5
votes
1answer
482 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 ...
5
votes
3answers
521 views
How to smooth / “reinterpolate” a raster with GDAL?
Short:
I want to change the resolution of a raster and smooth the gray ramp like shown in the images bellow. The preference is to use GDAL, PIL or Numpy.
Description:
I'm kriging points into rasters ...
6
votes
3answers
745 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
182 views
What Open Source (command line/API) method is there to orthorectify and georeference an aerial image given telemetry data?
How can I take an image and orthorectify it automatically with open source software?
The information associated with each image is:
camera specs
latitude, longitude, elevation of the camera
...
5
votes
2answers
499 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 ...
