3
votes
0answers
39 views

Subset two rasters with different extent and origin (and resolution) only where there are data in both in R

in R, I have two rasters with different extents and origins (and possibly even resolution). How to efficiently take the subset of each of them only where there are data in both of them? Here is a ...
3
votes
0answers
129 views

SAGA is not working within sextante (QGIS processing installed in Ubuntu)

Recently i decided to start working in Ubuntu (version 14.04). My first step was to add the ubuntugis-unestable repository, and then proceed to install SAGA, GRASS, OTB, QGIS and the bindings between ...
3
votes
0answers
51 views

point feature to raster with count of unique points

I am using ArcGIS 10.1 with Python 2.7 I have a point layer with thousands of points that each hold an id value of its owner. One owner can have multiple points, meaning that points can have the same ...
3
votes
0answers
43 views

GDALRasterizeGeometries doesn't burn correct Values

What I have is a vector shape file with a string field. What I've done is written a CSV file with the field ID, field_value and now I would like to rasterize the vector layer as a byte raster with the ...
3
votes
0answers
22 views

execute binary from grass shell?

I am new to grass and could really use some help with my problem. I am using hydrologic modeling software, and it include a grass program that I need to run from the grass shell. My question is, how ...
3
votes
0answers
22 views

Preparing data for semiovariogram analysis in ArcGIS with large sample sizes?

I am working on a statewide project to create surface water quality grids from point samples. In some instances, we have large sample sizes (n>50,000) which is problematic when creating ...
3
votes
0answers
67 views

Batch merging shapefiles from different folders, different subfolders

I tried some of the scripts I found but none worked so far. What I am trying to do is batch merging hundreds, maybe thousands of shapefiles from different folders/sub-folders into a single shp or ...
3
votes
0answers
69 views

DEM from contours

I have problem to create a DEM from contours. I know it is very simple but I can not make it work. I tried everything but nothing has happended! The shapefiles I've got are a con_lin and a con_pol ...
3
votes
0answers
69 views

Merge Raster Layer Specific Order in Graphic Modeler

I'm using the Graphic Modeler/Processing Modeler to automate a process. The last step in this process is Merging three different Raster Layers. I'm using the "Merge Raster Layers" from Saga 2.0.8. ...
3
votes
0answers
74 views

What is the openlayer.js alternative of unloadinvisibletiles of leaflet.js

I am using open layers with a tile server. I want to load maps while panning and destroy the loaded region when out of screen (Alternative to unloadinvisibletiles of leaflet.js) How am i suppose to ...
3
votes
0answers
55 views

Frontal Detection of sst & Chlorophyll image

I am working on sea surface temperature satellite images and looking for delineating gradients in open source on the basis of temperature. Is there any tool like Cayula and Cornillon front detection ...
3
votes
0answers
82 views

How to set up graduated colors for a layer using ArcObjects?

I have to render a feature layer with a graduated colors renderer, using the ArcGIS API for Java (or VB, .NET). As the feature layer data is created dynamically using webservices, I don't in advance ...
3
votes
0answers
57 views

Creating map gores in ArcGIS for Desktop?

Does anyone have clear instructions on creating gores for a globe in ArcGIS? I have tried to work with these: http://mappingcenter.esri.com/index.cfm?fa=ask.answers&q=1725 But have become ...
3
votes
0answers
46 views

Appropriate resolution of images for image classification?

I am attempting to perform image classification on some maps that I have in PNG format. This classification is failing due to there being too many many pixels in the training samples (I appear to be ...
3
votes
0answers
43 views

Is there a way to expose Spherical Mercator Lat Lon in Google Map API?

The map tiles in the Google Maps are in Spherical Mercator projection (srid 900913), but in the API the methods report coordinates in Lat/Lon WGS84, e.g. map.getCenter(). Is there a way to return the ...
3
votes
0answers
172 views

Geometric Network for Electric Utilities

I'm not expert in the use of geometric networks nor in Python scripting and I need help with the following problem. I'm using ArcGIS 10.0. I have defined a geometric network for an electric ...
3
votes
0answers
134 views

Problem importing KML file to OpenLayers - Projection System

I know that many have this problem,but none of the answers helped me so i suppose it's a different problem.This is my code: <!DOCTYPE html> <html> <head> ...
3
votes
0answers
70 views

kriging of compositional data (grain size)

I have posted this already in the GEONet Forum but without any suggestions, now I am hoping for some feedback here: I have a problem concerning interpolation of compositional data of grain sizes. I ...
3
votes
0answers
93 views

Preserve datatypes when using SQLite and QGIS

Following Nathan's SpatiaLite/DBManager solution in this answer, I find that once the layer is loaded from DB Manager into QGIS (2.6.1), any calculated fields are appearing as text values instead of ...
3
votes
0answers
215 views

Changing workspace path of query layer using ArcPy?

I'm trying to change the workspace path of a query layer using the code below and it always breaks. The layer is connected to a view in an Oracle Spatial 10g database. The database is SDE enabled ...
3
votes
0answers
89 views

How to deal with data which are not normally distributed when kriging?

I'm a novice with advanced spatial interpolation, and I'm trying to use the Geostatistical Analyst to interpolate a surface based on borehole data, using Ordinary Kriging. The help file explains how ...
3
votes
0answers
43 views

Subtracted 2 DTM layers to find the height of trees but all we have is the pixel values, can I convert this to height?

For a piece of coursework we have an aster layer that captures the tops of trees and a panorama layer that sees ground height (both are raster) so we have subtracted one from the other to get tree ...
3
votes
0answers
103 views

How to display the slopes that only facing west using ArcGIS for Desktop?

Here is what I am struggling with ... In ArcGIS for Desktop I have to Derive and symbolise areas considered susceptible to avalanche risk based on the following simple rules: areas with a slope ...
3
votes
0answers
73 views

Why do each of my buffer-polygons have different areas

I`m using Qgis 2.2 and 2.4 and loaded a wms-map using the coordinate system ERTS89 UTM 32N. I set 8 points (all in on shp-fileand in UTM 32N) on various locations on the map. After that i created a ...
3
votes
0answers
37 views

Prevent PostGIS Topology from splitting edges on line insert

I'm dealing with a large amount of road data represented with PostGIS Topology types. Anytime I add an edge over another edge, the edge is automatically split and a new node is created at the ...
3
votes
0answers
217 views

QGIS 2.6 Brighton - Graduated Symbols

I can't get graduated symbols to show in QGIS 2.6. This is what I have done. I have two fields of type integer, saved as a shape file which I have tried to import as a shape file and as a secondary ...
3
votes
0answers
48 views

How to geocode a building with different entrances

I am working on a routing map for a university campus. the routes are modeled and network analyst works fine. The problem: Many of these buildings have several entrances. Is there a way to geocode ...
3
votes
0answers
41 views

QGIS - assign different colors for SVG markers also in categorial / graduated symbols?

I understand it is possible to customize the color of an SVG marker via param(fill) However, I have to set the color manually, and cannot automatically assign it through categorical or graduated ...
3
votes
0answers
23 views

QGIS isn't displaying dates in the attribute table

I have a QGIS project that links to data in a Postgres database. For some reason, QGIS is not displaying any data in date fields. The date fields appear to be fine in Postgres and I have performed ...
3
votes
0answers
184 views

Is there any reason not to upgrade to Windows 8.1?

I have the option to upgrade to 8.1, but before I do I want to see if anyone has had any compatibility issues. I mainly run ESRI, but occasionally use QGIS, and other open source products. Has anyone ...
3
votes
0answers
80 views

ArcGIS Viewshed keeps crashing

I've been trying to run a viewshed analysis in ArcGIS 10.1 and everytime I've tried running it, it crashes the whole program, with no error. See the image below. I've made sure it has OFFSETA field ...
3
votes
0answers
56 views

Definition queries deleting themselves in ArcMap10.2

After switching to ArcMap 10.2 I've noticed that my definiton queries will seemingly randomly remove themselves. I have no idea what triggers the deletion since I have noticed it after reordering ...
3
votes
0answers
24 views

Labelling Max coverage by Area

I'm struggling to find a quick way to identify and label an area by the maximum coverage of an overlapping polygon. Has anyone done a process like this lately? I'm currently playing with the Pivot ...
3
votes
0answers
58 views

How to export two related layers to KML with multigeometry using QGIS?

From QGIS, I would like to create a KML file, for GoogleEarth, which contains multigeometry in that KML file. By example: consider two layers in QGIS -- the first layer contains polygons, and the ...
3
votes
0answers
31 views

How to set m parameter of clusterR function in R?

I want to manage a clustered computation on raster and thus use the "clusterR" function in R (package 'raster'). However, there's a tuning m parameter to set, but the only explanation about this ...
3
votes
0answers
81 views

How can I programmatically remove a layer from the QGIS composer legend without removing it from the map canvas?

I am working on a QGIS Python plugin which includes automated map production. I use a series of map templates, each of which has a main map, an overview map (with a layer set different to the main ...
3
votes
0answers
64 views

What is the area that OGR extracts from the (geometry) shp files for polygons?

The ESRI July 1998 technical documentation for polygon shapefiles at page 10 makes no mention of either area measures or length measures (apart from byte/record-layout matters). Yet OGR enables me to ...
3
votes
0answers
36 views

QGIS Cross Section Profile Distance using Ellipsoid?

I am working with QGIS 1.8 and 2.6 and can't get any of the Cross Section Profile tools to provide the ellipsoid ("proper") distances. All the distances use geographic coordinates which are too ...
3
votes
0answers
18 views

Action: Move pointfeature to other layer - Point becomes Multipoint

I'm trying to write an action that moves a pointfeature from one layer to another. It works fine until i try to save the layer with the new feature on. He says he can't save a MULTIPOINT to a ...
3
votes
0answers
72 views

Dilating a polygon by a circle sector

I want to dilate a polygon by a circle sector in ArcMap. Because my python skills are not that good to figure out how I can use the R prototype of How to create an oriented buffer using arcpy? . Is ...
3
votes
0answers
50 views

Diminishing performance of gdal_polygonize in a python subprocess for GRASS to PostgreSQL

I have written a python script to polygonize GRASS rasters into a PostgreSQL database using the gdal utility gdal_polygonize: sub.call('gdal_polygonize.py -8 gdb/quercus/{0}/cellhd/{1} -f ...
3
votes
0answers
182 views

SetNull in Python script for Arcgis assign 0 value

I am writing a script in Python where I assign to few Tiff defined in a function as inputRaster a NoData value as follows: whereClauseSetNull = "VALUE > 100" outSetNull = SetNull(inputRaster, ...
3
votes
0answers
48 views

MapServer - Os Open Data - Tutorial

I have recently received my OS OpenData discs from Ordnance Survey, and for the last two evening tried to get this data into MapServer for ultimate use as a WMS for OpenLayers. Safe to say I have ...
3
votes
0answers
46 views

Editing Cursor for GPS live tracking in QGIS 2.4

I am interested in using QGIS 2.4 for vessel navigation due to its ability to display NOAA vector charts or BSB's as well as user created shapefiles. I currently oversee a project that requires ...
3
votes
0answers
80 views

arcpy.mapping.ExportToTIFF fails to render layer?

Has anyone else had this problem? I have an mxd file with two layers in it, one of which is fairly 'heavy' (i.e. takes a minute or so to draw). All layers are visible by default. I have the ...
3
votes
0answers
33 views

Generate a series of ellipses on a line

I am trying to create a series of ellipses (mocked up in green below) of defined dimensions along a line. I have got as far as generating regularly spaced points along the line that are attributed ...
3
votes
0answers
238 views

Select and copy attributes features in ArcMap using Python add-in tool?

I'm trying to add a tool to an add-in toolbar to select features from an existing feature class and copy attributes to another feature class,from point to point. TopoHeight_Field_Name='Comment' def ...
3
votes
0answers
109 views

Does QGIS work with Anaconda?

I am trying to get QGIS to work with the Anaconda distribution of Python but I do not know how to get this to work. In one of the answers it was mentioned that QGIS does not work with Anaconda and I ...
3
votes
0answers
99 views

How to set Vertical Coordinate system when creating a 3D shapefile?

I'm working with GDAL v1.11 and am seeing the vertical CS stripped when exporting a 3D shapefile (from oracle spatial). This is the command I'm running: ogr2ogr -f "ESRI Shapefile" B_POINT.shp ...
3
votes
0answers
162 views

Hexagon binning using Google Maps

I have a bunch of geolocated points that I am looking to cluster on a google map using hexagons. These clusters will be stored in a db for quick retrieval and will be reduced to a single point and ...

15 30 50 per page