Using osgeo modules in ArcGIS10.1 and QGIS 1.8 python shells

This is related to an old question of mine - Issues with osgeo - qgis.core and arcpy modules Now in new versions of all the software QGIS Lisboa (1.8), ArcGIS 10.1 and pyScripter 2.5.3. They all use python 2.7 I can't get the following to work.               >>> import osgeo     >>> from osgeo import ogr      arcgis 10.1 and pyscripter fail here.               >>> import os     >>> os.chdir(r'L:\Data_Admin\QA\test_data\Boundaries\Local_Govt\TRC\TAB')     >>> driver = ogr.GetDriverByName('MapInfo File')     >>> fn='trc_boundary.tab'     >>> dataSource=driver.Open(fn,0)     >>> layer = dataSource.GetLayer(0)      > > > > feature = layer.GetFeature(0) extent=layer.GetExtent() extent > (-12895.844282306107, 376810.87105960253, 7944810.84780738, > 8222992.984638059) numFeatures=layer.GetFeatureCount() geometry = > feature.GetGeometryRef() Traceback (most recent call last): File "", line 1, > in AttributeError: 'NoneType' object has no attribute 'GetGeometryRef' QGIS python shell fails here. Not sure why - I checked http://www.gdal.org/ogr/classOGRGeometryCollection.html but it is not helpful. my > PYTHONPATH=C:\Python27\ArcGIS10.1;C:\OSGeo4W\apps\python27; C:\Program Files > (x86)\Quantum GIS Lisboa I have installed the gdal python bindings for arcgis but then where is ogr? also I can't seem to run               from osgeo import gdal     dataset = gdal.open("c:/abc.img")      > Traceback (most recent call last): File "", line 1, in AttributeError: > 'module' object has no attribute 'open'


OpenLayers GetFeatureInfo for GeoServer Layer Groups

I've created a GeoServer layer group. I'm using `OpenLayers.Layer.WMS` to retrieve the layer group, and this works fine. In addition, I'd like to pull feature info for the layer. For this, I'm using OpenLayer's `GetFeatureInfo`. But the `GetFeatureInfo` only returns feature info for only one of the layers in the layer group. Is there a way I can retrieve feature info for both layers in the group? And, is this an optimal way of dealing with layer groups? A section of my code is shown below:                   var layer_group = new OpenLayers.Layer.WMS(             "Sample Layer Group",             "http://localhost:8081/geoserver/wms",             {layers: "layer_group_name", transparent: true},             {isBaseLayer: false, visibility: false}         );         map.addLayers([layer_group]);      to retrieve the layer group. And:                   map.events.register('click', map, function (e) {             var url =  layer_group.getFullRequestString({                 REQUEST: "GetFeatureInfo",                 EXCEPTIONS: "application/vnd.ogc.se_xml",                 BBOX: layer_group.map.getExtent().toBBOX(),                 X: e.xy.x,                 Y: e.xy.y,                 INFO_FORMAT: 'text/html',                 QUERY_LAYERS: layer_group.params.LAYERS,                 WIDTH: layer_group.map.size.w,                 HEIGHT: layer_group.map.size.h});                  window.open(url,                     "getfeatureinfo",                     "location=0,status=0,scrollbars=1,width=600,height=150"                 );         });      To retrieve the feature info. My environment setup: GeoServer 2.1.3, PostGIS 1.5, OpenLayers 2.11. EDIT: The features shown are only for one layer ![enter image description here](http://i.stack.imgur.com/ABsBe.png) Increasing `FEATURE_COUNT` only increases the number of instances(rows) returned. The instances are over 1M. What I'd like to achieve is have one row having feature info for all the layers in the layer group.


How to select features in OpenStreetMap

I'm using OSM (online) with OpenLayers to show my map (I don't have a mapserver) but when I do this:               map = new OpenLayers.Map( 'map',  {numZoomLevels: 7, units:"km", minZoomLevel: "10" } );     map.addLayer(new OpenLayers.Layer.OSM("Main Map"));     map.zoomTo(3);      The map comes with too much information. I need to show only roads and countries. Any way to do this?


How to create paths within a boundary?

I am looking to create movement paths of fish based on observed locations. Since I am looking at movements of fish within rivers and lakes, simply connecting the points to form a polyline will not work as many of the paths would be over land. I need some way to restrict the movement paths to within the water boundaries. I am not a programmer and rely on the tools in the Arc toolbox to perform the operations. A new column in the FAT with distance between sequential points would also be useful. I am using ArcMap 10. Does anyone have an suggestions as to how to generate the movement paths? UPDATED 9/19/12 Thanks so much for the assistance, many of the suggestions look promising. A bit more information on the data; the majority of the locations are collected from remote receivers which tend to generate a lot of observations (one of my projects has over 3 million detections and growing, each detection record includea ID, datetime, lat and long). If a fish stays within the detection range of a receiver it is detected roughly once every two minutes which results in a lot of observations (points), so my first step is to average the lat and long over a time period in this case 1 day. Averaging the detections in the lake habitat typically works well, but doing so in the river portions can result in detections outside the river/lake boundary. So the first thing I need to do is "snap" the averaged locations into the River/Lake Boundary and then I would like to create a movement path that is restricted to within the river/lake boundary, ideally this path would include distance for each line segment. My rivers are polylines and the lakes are polygons, but I can convert the rivers to polygons if necessary. Alternatively, I could use the raw unaveraged data and all my detections would be in the boundary, although creating movement paths restricted to within the boundary is still problematic, but this would result in a lot of data, maybe too much for ARC? The attached image shows the daily average positions for two different fish. ![enter image description here](http://i.stack.imgur.com/j65pZ.jpg)


Upright/Horizontal lables when labeling polygons at perimeter?

I have several polygons (circles) I want to label with their ID number. The label should be outside the circles so I use the positioning options "Use perimeter" and the tickbox "Next to the line (or similar, my qgis talk spanish...). QGIS now automatically aligns the labels with the curvature of the line. Is there an option or a way to make them oriented horizontally? Martin Using QGIS 2.2


Problem using pgRouting 2.0

I follow this great post:pgRouting 2.0 for Windows quick guide, and I have loaded the test data from GeoJSON, like this: ![enter image description here](http://i.stack.imgur.com/UTZNZ.png) But I don't know how to load this layer into my database, and my DB manager doesn't look the same as in Anita's post: ![enter image description here](http://i.stack.imgur.com/wyHJi.png) How should I do to load and test this layer into my PostGIS database?


How to run a 2.0 ESRI DLL in a 4.0 Addin application

We have a complex set of operations running as Add-in's and services that span both ArcMap, ArcServer, and ArcGIS Explorer for Desktop. Currently we are compiling everything under .NET 3.5 but we may need to move to .NET 4.0 to resolve an issue. All of our code uses ProgressHelper from Explorer for Desktop 10.1 which is strictly a 2.0 framework applicaiton. Normally I'd add the command to the config file but the Add-In's are under ESRI control, not the .NET framework. Is there a way to modify the ESRI configuration files so that the .NET framework knows that 4.0 is the supportedRuntime environment and allows me to load the 2.0 framework DLL? As a testing note: ArcMap does have it's own configuration file and in that I do have the Startup section defined for both the 2.0 framework and the 4.0 framework.


Why does ArcMap edit session show different vertex count on feature than arcpy/ ArcToolbox

I am looking at features in an ArcSDE geodatabase via ArcMap. The feature shows to have four vertices by viewing the vertices in the map and by looking at the sketch properties. However if I run a feature vertex to point tool i get 5 vertices. If I get an arcpy.SearchCursor row object and do row.pointCount I get 5 vertices reported. I even exported the feature class out to a shapefile and got the same results. Whats going on here? Thanks.


How to create label with variable size in ArcMap?

I'd like to create a circle label in ArcMap in the middle of each Polygon feature in my map document to shows the value of a particular attribute associated with that feature (like number of available stations). I'd like to represent the number of stations by the size of the label, so that the larger number of station in a particular feature is represented by the larger circle in that feature. Is there any possibility in ArcMap to do so? any help would be much appreciated.


How to get ESRI Addin to read the System.ServiceModel configuration?

I have a complex set of add-in's that we need to run in ArcMap, ArcServer, and ArcGIS Explorer for Desktop. We are experimenting with changing from the old Web Reference service connection to a WCF Service Reference connection. When I try to create the client it cannot find the configuration settings to create the client. When I add just the URI to the embedded Settings.settings file, and create an EndpointAddress and Binding object to pass to the client, I get a client but my first attempt to access the server fails. So how do I access the WCF service configuration data from inside an ArcMap Addin?


Why are some values missing from the imported CSV?

Bit of a strange one: i'm importing a CSV where the rows are the WKT data for regularly progressing rectangles from a grid, then the rectangle name, then values (ranging between 2 & 254). When I import the CSV file it displays correctly as 'simple symbol' with field names, i.e. it can bulk-colour the rectangles and name them, but the data aren't available for (e.g.) graduated colouring. HOWEVER: if I use 'Categorized' as a style then classify, the specific value classes populate, including a blank final field. When I apply, ALL the rectangles are the blank. Checking metadata for the value row, the results for column names 'type / type name / length / precsion' are double/double/0/0. For a VERY similar WKT CSV import where this worked, those values were double/Real/24/15. Any thoughts why this might be? (v2.1 alpha on Win7. Booting into linux to see if it's a master build problem) EDIT: works fine on 2.1.0-master on ubuntu. My windows version is MORE up to date, however... Screenshot of the CSV ![here](http://i.stack.imgur.com/1ZWXQ.png) Thoughts: absent cells in columns 'i' onwards, commas form part of the WKT but it's saved as a CSV (however it seems to save as tab separated, since it opens fine and the cells display fine). Further info / visual explanation: even after saving the excel file as tab delimited text (to avoid comma problem risk) AND putting the WKT as the final column just in case, the problems persist. In 'fields', none have any value but all rows (corresponding to excel columns) are present![.](http://i.stack.imgur.com/rDgZr.png) However in 'styles' one can see that the 'Grand Total' column - like all columns - produces a mystery blank field (red) which all the cells get coloured as, AND there's a mystery blank data column after 'Grand Total' in the dropdown menu which seems to come & go![.](http://i.stack.imgur.com/Po6f0.png) All super weird, and seemingly only in windows. I'm not sure if importing and presenting correctly in ubuntu then saving the project & opening in Windows will cause a problem (I dual boot & have Qgis on both). Unrelated, but while I remember: I have to change the years to "Year: YYYY" rather than just YYYY as they were before, otherwise QGis doesn't like it, and changes them to "column i", etc.. I know this is a separate question, but just in case anyone knows offhand if there's a solution. EDIT: link to file which suffers similar problem but isn't the same file: For ref: this problem isn't there in 1.8.0. with the same files.


Zoom to selected features with ArcObjects

I am trying to figure out what is the best way to set the active view extent to be able to view all the selected features on the Arc Map. 1.The Map has only one layer. 2.The features are filtered with attributes. Below is the code where I tried,stuck with the part where I can loop through the features in the feature selection set but not able to set the active view extent to zoom to show all the selected ones. P.S: this is similar to doing the right click on the attributes table and selecting multiple rows and doing a Zoom to selected.               private void ZoomToMultipleDGVSelection(List<int> selectedRightIDs)         {                  IMxDocument pMxDoc = ArcMap.Document;             IMap pMap = (IMap)pMxDoc.ActiveView;                  ESRI.ArcGIS.Carto.ILayer layer = GetLayersClass.GetFieldBoundaryLayer;             if (layer is ESRI.ArcGIS.Carto.IGroupLayer)             {                      ESRI.ArcGIS.Carto.IGroupLayer groupLayer = layer as    ESRI.ArcGIS.Carto.IGroupLayer;                 ICompositeLayer pCompositeLayer = layer as ICompositeLayer;                 int layers = pCompositeLayer.Count;                 ILayer pLayer = pCompositeLayer.Layer[0];                 IFeatureLayer pFeatureLayer = (IFeatureLayer)pLayer;                      IFeatureSelection pFeatureSelection = (IFeatureSelection)pLayer;                 ISelectionSet pSelectionSet = pFeatureSelection.SelectionSet;                      IFeatureClass pFeatureClass = pFeatureLayer.FeatureClass;                 IQueryFilter pFilter = new QueryFilterClass();                      foreach(int ID  in selectedRightIDs)                 {                 pFilter.WhereClause = "RightID = " + ID.ToString();                      IFeatureCursor pFeatureCursor = pFeatureClass.Search(pFilter, false);                 IFeature pFeature = pFeatureCursor.NextFeature();                      pFeatureSelection.Add(pFeature);                      }                           //if (pFeature == null)                 //{                 //    System.Windows.Forms.MessageBox.Show("This section doesn't exist");                 //    return;                 //}                           IGeometry pgeom = (IGeometry)pFeature.Shape;                 pMap.SelectByShape(pgeom, null, false);                 IEnvelope pEnv = pgeom.Envelope;                 pMxDoc.ActiveView.Extent = pEnv;                 pMxDoc.ActiveView.Refresh();             }         }      Thank you in advance!


Can you make a travel region polygon with with Google Maps API?

I am working on a local design project and need to determine if someone is local to my area to allow them to participate. Local, as defined by the group, means that the person is within 2 hr driving time based on Google maps to 32 Vassar St., Cambridge, MA. While I can figure out how to get travel times from point to point, I cannot seem to figure out if there is a way to create a polygon that displays the region which we consider local. Any ideas on how to do this using Google tools?


What tools in Python are available for doing great circle distance + line creation?

I need to use Python to create a great circle distance -- both a number, and preferably some kind of 'curve' that I can use to draw in a client-side map. I don't care about the format of the curve -- be it WKT, or a set of pairs of coordinates -- but just want to get the data out. What tools are there out there? What should I use?


How to use a variable in Geoserver viewparams?

I just managed to create a parametric view that is working ok when i define the parameter value as, for instance, **viewparams:'value:31'** However, when i try to pass a value that comes from a variable, for instance, **viewparams:'value:variable'** , the wms layer just wont load. The variable is defined as **var variable= 31;** Any ideas on why this happens? I'm using OpenLayers + Geoserver + PostGIS Thanks EDIT 1 - Just to explain a little bit better what i'm trying to achieve. I have a web map application developed with Openlayers where i have a set of wms layers that are being created by Geoserver. Despite the fact that i already have some functionalities implemented, such as feature info in pop ups, images on those pop ups, etc, i thing that it would be really important to allow users to query the layers, to allow the map to show only what they wanted to see. The query would be done through a php form, where the values selected by users would act as variables on the Parametric views created on Geoserver, and by doing so, the features showing up on the map would be only the ones selected by the user. The parametric view i created for testing purposes works, if i "hardcode" the value for the parameter:               > var rov_dives = new OpenLayers.Layer.WMS(     >            "ROV Dives",     >                    "http://host:8080/geoserver/wms",     >                    {'layers': 'parametric_view', 'format':'image/png', 'transparent':'true', viewparams:'test_parameter:50'},     >                    {'opacity': 1, 'isBaseLayer': false, 'visibility': false}     >                 );      In this example the map only shows the features that comply with the condition stated on the sql script that defines the parametric view, in this case, all features that have a number of measurements higher than 50. If i could pass that 50 as a variable coming from the aforementioned php form, i could implement the solution i thought of... The definition of the variable and WMS layer i'm trying to use right now is:               var variable = 50;          var rov_dives = new OpenLayers.Layer.WMS(                "ROV Dives",                        "http://host:8080/geoserver/wms",                        {'layers': 'parametric_view', 'format':'image/png', 'transparent':'true', viewparams:'test_parameter:variable'},                        {'opacity': 1, 'isBaseLayer': false, 'visibility': false}                     );      Despite all my efforts i'm yet to find a solution for this. Hope this explanation helps to understand the problem! thanks in advance, Filipe


Where can I find Canadian postal code geometries?

Where can I find the shapefiles for Canadian postal codes?


Display Arabic text in ArcGIS

I have my GDB where names for a few fields in feature class is written in Arabic. When I label them, it is displaying the labels but unable to understand the text as it is symbolized/coded. Can anyone suggest how to add Arabic text and display the labels properly to make them readable?


What kind of coordinate system has values like 48569.800mE, 706602.032mN?

I am seeing the following set of coordinates, and I was wondering what kind of system is this? `48569.800mE, 706602.032mN`


Query to find all lines connected to other

I have problem with one spatial (Oracle) query. Simplified overview of the problem is: I have two tables:               CREATE TABLE thin_lines (       id NUMBER PRIMARY KEY,       geom SDO_GEOMETRY),       color VARCHAR2(32);          CREATE TABLE thick_lines (       id NUMBER PRIMARY KEY,       geom SDO_GEOMETRY),       color VARCHAR2(32);      Both tables contain 2D lines. The example picture: ![enter image description here](http://i.stack.imgur.com/7U9ZB.jpg) The problem: I want to create SQL query or PL/SQL procedure to find ids of all blue thin lines (whole blue branch comprising of blue lines) connected to the thick green lines.


Layer just for annotations

Im relatively new to QGIS and would like to know whether or not I can save my form annotations to a layer so I can turn them on or off.


Unable to access GeoServer from another computer in the same domain

I installed GeoServer 2.5.2, GeoNetwork 2.10.3 and Tomcat 7 on my Windows 7 computer {MyComputer} in our company's network. Now I would like to access GeoServer, GeoNetwork and Tomcat from other computers in the network. However, I could only access via the web browser GeoNetwork and Tomcat but not GeoServer. **NO** \- http://{MyComputer}.{domain}.com:8070/geoserver/web/ **YES** \- http://{MyComputer}.{domain}.com:8040/geonetwork/srv/eng/main.home **YES** \- http://{MyComputer}.{domain}.com:8080 Is there any configuration that I missed?


Working with Fulcrum data in QGIS

I have been sent some shapefiles, gathered using the Fulcrum app. However, when I load them into my project on QGIS 2.2, they wont display. They will display if I go to 'Zoom to Layer Extent" but under no other circumstances - even if I remove all other layers in the project.


How can you tell if a RDBMS is an ArcSDE Geodatabase?

I want to know if there is any way I can query the _RDBMS_ to determine whether it is an _ArcSDE Geodatabase_ , or if I have to do it through _ArcObjects_. I am wondering if it would be enough to check if the `SDE` schema exists on the target database, and then maybe check if the `SDE.VERSION` table also exists and then conclude that it's an _ArcSDE Geodatabase_. The _ArcSDE Geodatabases_ are `Oracle`, and I use `ArcObjects SDK 10 .NET` in `C#`.


converting vectors to raster (burning streams)

I converted my stream vector to a raster file (using polyline to raster conversion tool). I anticipated that I would get a No Data row in the raster file attribute table, so that I could use a conditional statement to burn the streams into the DEM: con(stream_raster>0, (focalmin_dem - 2), original_dem) However, since I did not get a No Data row, all this gives me is a a DEM img directly beneath the streams (not the remaining extent of the image). What am I missing?


OpenGeo GXP Examples

Hi looking at GXP examples here They do not list the dependencies within the JS file, so I was wondering how to make the files work. I would assume just copy the js file, e.g. viewer into my SDK app.js file and build that. However it just loads a white screen when I load it. Any hints on how to sue these examples? Thank you


overriding showMap

This code is not working as needed. I need just to make a window with addlayer button.               import java.awt.BorderLayout;     import java.awt.event.ActionEvent;     import java.awt.event.ActionListener;     import java.io.File;     import java.io.IOException;          import javax.swing.JButton;     import javax.swing.JFrame;     import javax.swing.JPanel;     import javax.swing.WindowConstants;          import org.geotools.data.FileDataStore;     import org.geotools.data.FileDataStoreFinder;     import org.geotools.data.simple.SimpleFeatureSource;     import org.geotools.map.FeatureLayer;     import org.geotools.map.Layer;     import org.geotools.map.MapContent;     import org.geotools.styling.SLD;     import org.geotools.styling.Style;     //import org.geotools.renderer.GTRenderer;     //import org.geotools.renderer.lite.StreamingRenderer;     import org.geotools.swing.JMapFrame;     import org.geotools.swing.JMapPane;     import org.geotools.swing.data.JFileDataStoreChooser;                    public class App2 {          public static void showMap(final MapContent map) throws IOException {                 final JMapPane mapPane = new JMapPane(map);           // mapPane.setMapArea(map.getLayerBounds());                 JFrame frame = new JFrame("cis");            frame.setLayout(new BorderLayout());            frame.add(mapPane, BorderLayout.CENTER);            JPanel buttons = new JPanel();                /* JButton zoomInButton = new JButton("Zoom In");            zoomInButton.addActionListener(new ActionListener() {public void actionPerformed(ActionEvent e) {  mapPane.repaint(); } });            buttons.add(zoomInButton);                 JButton zoomOutButton = new JButton("Zoom Out");            zoomOutButton.addActionListener(new ActionListener() { public void actionPerformed(ActionEvent e) {  mapPane.repaint();  } });            buttons.add(zoomOutButton);                 JButton pamButton = new JButton("Move");            pamButton.addActionListener(new ActionListener() { public void actionPerformed(ActionEvent e) {   mapPane.repaint(); } });            buttons.add(pamButton);            */            JButton addLayer = new JButton("Add Layer");            addLayer.addActionListener(new ActionListener()  { public void actionPerformed(ActionEvent e) {                try{                    File file = JFileDataStoreChooser.showOpenFile("shp", null);                 if (file == null) {                     return;                 }                 FileDataStore store = FileDataStoreFinder.getDataStore(file);                 SimpleFeatureSource featureSource = store.getFeatureSource();                 Style style = SLD.createSimpleStyle(featureSource.getSchema());                 Layer layer = new FeatureLayer(featureSource, style);                 map.addLayer(layer);                  }catch(IOException f){}             } });            buttons.add(addLayer);                 frame.add(buttons, BorderLayout.NORTH);            frame.setDefaultCloseOperation(WindowConstants.EXIT_ON_CLOSE);            frame.setSize(600, 400);            frame.setVisible(true);         }              public static void main(String s){                     MapContent map = new MapContent();                 map.setTitle("Climate information system");                 JMapFrame.showMap(map);                        }     }


Software needed to recreate VBA forms using python for Arc 10.0

I'm new to programming and need to convert VBA scripted tools for the 10.1 upgrade. I need to recreate the forms and use Python scripting. How do I do this? I need the form to function in 10 so that the upgrade will go smooth. A list of specific software that I will need would be so helpful.


Stop custom attribute form from resizing in Qgis

When I open a custom attribute form for a feature in Qgis it always gets resized to the default size of the original Qgis attribute form. I've tried to set the sizePolicy in _QT Designer_ to _Fixed_ instead of _Preferred_ to prohibit this, but it didn't help. If I resize the form when I open it by clicking on a feature in Qgis and then resize it, it still opens with Qgis default height/width when I click the next feature. Is there any kind of setting to tell Qgis to use the size that the form is designed to? (Running Qgis 1.8.0, QT Designer 4.8.4 on Windows 7)


How do I import AutoCAD [DWG] drawings to QGIS [SHP] using the right CRS settings?

We did some topographical surveys in Indonesia which are drawn in AutoCAD. It's our aim to publish them with QGIS and its webclient. I thought the best way is to georeference the data in AutoCAD 2014. Eventually, I opened the DWG-file in AutoCAD Map 3D 2014 (see AutoCAD- Map-3D-Screenshot) because I know of it's possibility to export as SHP-file. Finally, I imported the SHP-file in QGIS (nightly built) but the result is not the expected one, see QGIS-screenshot : The white dot(s) represents GPS-data which was correctly imported and rendered in the southern part of Nias. The blue dot(s) represents the imported SHP- files which should actually be located at the same location but are -- as you can see -- in the middle of nowhere. Do you have any clue about a mistake I made? While exporting/ importing the data I payed attention to the CRS and set it to "DGN95.UTM-47N" in AutoCAD Map 3D and to "EPSG:23867 - DGN95 / UTM zone 47N" in QGIS -- which, in my opinion, should be the same one and lead to a right result. Unfortunately, it doesn't. Or is there something I missed? I already thought about the "Data source encoding" which I can set in QGIS. Right now it's set to "System". I guess Debian Wheezy (which I use for QGIS) and Windows 7 (which I use for AutoCAD [Map 3D]) use a different one -- but I'd also suppose that I would not see anything at all if I import data with a wrong encoding. Some ideas? Any help is very appreciated. Thank you.


GEOS: how to access Noder functionality through the C API

All this stuff looks rather useful: http://geos.osgeo.org/doxygen/classgeos_1_1noding_1_1Noder.html However, from the GEOS README file: `The C++ interface should not be used directly; the geos project views it as a bug for another program to use the C++ interface or even to directly link against the C++ library.` Is there a C interface to the Noding functionality that I'm missing? I can't find it in `geos_c.h`.


How to add a custom SymbolLayerV2 to a Vectorlayer in QGIS

I'm trying to add a symbolLayer with a line-decoration (line-arrow) to a vectorlayer of linestrings. what i am trying at the moment is to create a new layer and set it's QgsLineSymbol with a list of symbollayers.               # Symbology             simple_props = { 'width' : '1', 'color' : '255,255,255' }             deco_props = { 'width' : '0.26', 'color' : '0,0,255' }                  simple = QgsSymbolLayerV2Registry.instance().symbolLayerMetadata("SimpleLine").createSymbolLayer(simple_props)             deco = QgsSymbolLayerV2Registry.instance().symbolLayerMetadata("LineDecoration").createSymbolLayer(deco_props)             s = QgsLineSymbolV2([simple, deco])                  self.ll.setRendererV2( QgsSingleSymbolRendererV2( s ) )       Every time QGIS trys to set the QgsLineSymbolV2 it crashes with a win32-exception. what am i doing wrong, or is there a possibilty to just add a new symbollayer to a existing vectorlayer?


How to calculate Gradient from DEM?

I have a digital elevation model for a catchment area and I would like to calculate the gradient by cells. Is that possible and does anybody know the steps?


Converting a Map to graphics and annotations

i want to export a map for later editing before print, so i thought about converting all the annotations and and graphics and save them in the .mxd map. when i export graphics after exporting labels , the graphics cover some labels, and there is no option to make the graphics in back level. but when i export the graphics before exporting labels, the labels makes some displacement due to some maplex config, dispite making a maplex lock on labels, it cannot be fixed. i would like to know if there's a way to manipulate graphics and annotation as layers despite saving them in the map. if there is an alternative solution to edit map in a CAD like style before print is welcome. i'm using ArcMap 10.1.


Closest Facility/Delete Identical Tool Help

I am having a bit of a time figuring out a work around for a problem I have. Essentially, I ran a Closest Facility operation on a set of about 200 points. Each point is an entrance to a building, and the underlying network is walking paths. I have named each of my facilities and incidents based on building number, and ran the operation so that each building was a facility as well as an incident, producing routes to and from each building. What I am looking for is for any given path, (1-2, 1-3, etc.), I would like the shortest route (in time). In the Route output, there are routes for each entrance. For example if building 1 has four entrances, there will be four routes from building 1 to building 2. I only want the shortest (time-wise) for each route. This is where I thought I would be able to use the Delete Identical tool. I advanced sorted the table, first based on the name of the route, then on the travel time. After this sort, the first unique record for each route name (1-2) had the shortest time. However, Delete Identical appears to delete all but the lowest Object_ID. This is a problem because as the routes were created and subsequently exported, they did not go in the order that I want now, so the lowest Object_ID is not always the shortest time for a specific building to building route. Any suggestions/comments would be much appreciated. Thanks, eric I'm using ArcGIS 10, with an ArcInfo license.


How to add a toolbar via python?

Through tutorials I learned how to add a toolbutton to the plugins-toolbar via python. Now I wonder how to add a complete toolbar with toolbarbuttons via python. Can anybody give some example-code? Thanks and regards


Simple script in ArcPy to analyze building proximity to restaurants

I am trying to make a script that will go through buildings in a shapefile row-by-row and do a count of points within .25, .5,...1 mile of it. The other points are restaurants.               import arcpy       arcpy.env.workspace = "C:\Users\zafard\Desktop\Building.gpd"            #need to do 4 iterations for the distance radii     buffDist = [0.25,0.5,0.75,1]                                   for Dist in buffDist:         print ("FOR A %s MILE RADIUS: " % (Dist))         inTable = "C:\Users\zafard\Desktop\Building\XYBuild.shp"              #set up cursor in table         rows = arcpy.SearchCursor(inTable)         row = rows.next()              #set up loop to analyze every row of the table         While row:             print row.OBJECTID                    #use select by distance tool                 arcpy.SelectLayerByLocation_management("FRest","WITHIN_A_DISTANCE", row, Dist, "NEW_SELECTION")                  #count the hopefully selected points for the layer.             count = int(arcpy.GetCount_management("FRest"))                  print count             row = rows.next()      I am now getting an error with running the tool in the while loop. "Object: Error in executing tool" Please help! This is mindwrakking!


How to generate legend in leaflet webmap?

Although I have checked the "Create Legend" option in the **qgis2leaf** plugin, it did not generate any legend in the output webmap. What am I missing here ? ![enter image description here](http://i.stack.imgur.com/JbBA8.png)


Set QGIS (on android and/or linux) to use a custom .qgs file by default when opened.

Need to have QGIS open with all my settings/layers/projections etc by default. Created a .qgs file that has all this info in it. Using Android and Ubuntu. Tried a shortcut to the .qgs file in android and that didn't work.


How to create a report of all feature classes in a geodatabase?

I am not getting any response from this script and no error displayed at all I am writing this script to help me display the feature classes in my geodatabase Here is my script               import arcpy     arcpy.env.workspace = r"Q:\basemaps\ALTA\FileGeodatabase\AB_IHS_GIS_DATA.gdb"     logfile = open(r"C:\tmp\log.txt ","w")     fcLst = arcpy.ListFeatureClasses()     for fc in fcLst:          logfile.write(fc + "\n")          logfile.close()


Combining Geological Grids

I have a situation where I have 3 Geological grids, each showing a different horizon + the Seabed Grid. I need to extrude between these grids in order to make the 3D polygon for the geological unit which is between these horizons. The problem occurs however when (due to Shallow gas or lower geological units protruding upwards) there are gaps in my grid; sometimes to the extent that one of the horizons only covers a small segment of the sea bed area at the surface. The solution therefore (I think) is to combine say, H3 + H2 for example, and H2 +1, and then: Extrude between [(H3+H2) (where H2 fills in the gaps of H3)] and [(H2 +H1)(Where H2 Fills in the Gaps of H1)]. The Situation becomes very complicated when some of H3 actually rises to the Seabed and the Seabed fills in some gaps in H3. A very complex series of combinations of grids is therefore needed involving probably H3, H2, H1 and the Seabed. The only way I have found to do this is in Global Mapper, but version I have does not have the combine grids tool, unlike the expensive newest version. The software I Have access to include Arc GIS 10.2 and freeware such as Map window, QGIS, GRASS etc. The solution may involve raster calculator conditional function, however I have no experience in how to accomplish this task using this. Any help is much appreciated.


Can someone explain me how to do the layout of routes in ArcMap? I have a shape of points with garbage containers and want to do the route of each truck but don't know how. I've tried with the network analyst, created a network dataset and I can't trace the route.


std::bad alloc when writing simple script in QGIS

I'm running QGIS 1.8 on Scientific Linux and I'm attempting to write a script (obviously this will be built upon to do something useful) which simply writes to a file.               ##raster=raster     value=5     ##value=output file      I select the raster, and the output file and I get the following; ![exception](http://i.stack.imgur.com/py8dI.png) and it crashes QGIS. I've seen countless bugs mentioning this, but nothing describing a fix. Is this something fixable? **EDIT** Just as an update with further information. It's not just the scripts I'm writing, it's any script that requires me to output a file. So I have just tried the Landscape Ecology Statistics package and attempted to output a file and that did the same.


How to use Field Value from arcpy.GetParameterAsText()?

I want to to is have a user select a feature class, which then populates the fields that can be selected (this works), create a new field in that feature (works), and then populate the new field with values based on the inputField values. Right now, it just populates the new field with the field name of the inputField. It seems like this question has been asked a thousand times, but I still can't get this to work.               import os     import arcpy     from arcpy import env     env.overwriteOutput = 1            inputFC = arcpy.GetParameterAsText(0)     inputField = arcpy.GetParameterAsText(1)     inputFieldLabel = "GeoLabel"     output = arcpy.AddFieldDelimiters(inputFC, inputField)     inputFieldName = arcpy.GetParameterAsText(2)          #Add the GeoLabel field to the Selected Feature Class     arcpy.AddField_management(inputFC, inputFieldLabel, "TEXT")          #Calculate the field values for the selected field, using UniueID field and InputFieldName     arcpy.CalculateField_management(inputFC, inputFieldLabel, output, "Python_9.3")


QGIS crashes when digitizing new fetaures

I can't currently digitize polygon features in QGIS 2.01 or 2.3. I can draw a set of points but as soon as I right click to finish the system crashes. It still works fine in 1.8. Windows 7 64


How to return a map catalog from a WPS process with GeoServer?

I guess everything is in the title, I would like to be able to select a map and have it returned.   This is because I'd like later on to return a GeoTIFF which would have been generated by the server, but I'm making small steps.   I already have a process that works fine on the GeoServer WPS demo, but when trying with other clients, instead of offering me a list of maps I have hosted to choose from, I am asked for a URL. (I have no idea as to what I should enter as URL, and I don't care since it is not very user-friendy...)   I think that the problem here is that GeoServer "cheats" when using his own demo client, because it knows the hosted maps, which other clients don't. I am therefore looking for the way, either for the process to return a catalog of the available maps (if it works), or another solution you could offer me. Thanks in advance !


Mapnik not generating tiles for tilecatche

imported data by using Osm2pgsql, to get rid off all problems related to data conversion, i stored this as latlong (-l parameter).I created stylesheet, and i tested with python script and it was working (so both data and stylesheet should be fine). Here comes problem..., when i tried to use it with open layers, i got my basic layer covered by pinkish "image not found" tile, but when i try to access tile by calling cgi script "http://localhost/tc/tilecache.cgi/1.0.0/osm/5/31/26.png" i have nicely generated tile is any one have idea what i'm doing wrong? files: osm.xml                   <Map background-color="#00000000" srs="+proj=latlong +datum=WGS84">       <Style name="roadStyle">         <Rule>                <LineSymbolizer stroke="#ff0000" stroke-width="100" />              </Rule>       </Style>            <Layer name="roads" srs="+proj=latlong +datum=WGS84">         <StyleName>roadStyle</StyleName>      <Datasource>          <Parameter name="table">            (select way,highway from planet_osm_line where highway in ('motorway','motorway_link','trunk','trunk_link','primary','primary_link','secondary','secondary_link','tertiary','tertiary_link','residential','unclassified') order by z_order) as roads              </Parameter>         <Parameter name="type">postgis</Parameter>              <Parameter name="user">www-data</Parameter>         <Parameter name="password">gis</Parameter>         <Parameter name="dbname">gis</Parameter>         <Parameter name="estimate_extent">true</Parameter>            </Datasource>       </Layer>          </Map>      tilecatche.cfg               [cache]     type=Disk     base=/tmp/tilecache          [osm]     type=Mapnik     mapfile=/var/www/tc/osm/osm.xml     extension=png               [basic]     type=WMS     url=http://labs.metacarta.com/wms/vmap0     extension=png      definition of my layers               var map = new OpenLayers.Map('map',           {             numZoomLevels: 12,             projection: new OpenLayers.Projection("EPSG:4326"),             displayProjection: new OpenLayers.Projection("EPSG:4326")           });               var layer = new OpenLayers.Layer.WMS( "VMap0",                      "tilecache.cgi?", {layers: 'basic', format: 'image/png' } );     var tms = new OpenLayers.Layer.WMS( "osm","tilecache/tilecache.cgi?",           { layers:'osm',  format: 'image/png'},{ isBaseLayer: false}  );


How should I connect from the ArcGIS Desktop to a virtual machine's ArcGIS Server?

Hy everybody! I have made an OpenVPN connection to an "inner network", my school's private network. I can use the remote desktop application after the connection, so it's surely work. The ArcGIS Server is in a virtual machine, in a Windows Server 2008 R2. The Server's version is 10.1 (I did not have refreshed that yet), and the ArcGIS Desktop version is 10.2. The Desktop application is installed on my computer. I have added myself as administrator on the ArcGIS Server Manager. How should I connect to the server? (add that GIS server to the ArcCatalog) I have wrote in the URL:               http://nameofthecomputerwithArcGISServer:6080/arcgis      But it says it cannot connect to the server. What should I do now?


Trouble with ArcGIS Server 10.2 web adaptor when Amazon Image start/stop

Every time we stop our ArcGIS Server 10.2 Amazon Image we need to reconfigure the web adaptor. Anybody have the same trouble or know of a solution? Thanks in advance.


Defining projection in arc toolbox (JSON to shapefile) in python

I have a arc toolbox which converts JSON txt file to a shapefile but we have to define the projection. How can I hard code the defined projection so that user don't have to input the projection ? I need to define projection parameter to web mercator and to add a function to reproject the shapefile to NAD 83 UTM. I am editing in Python.


QGIS adding or removing values in the style tab (layer properties)

I'm new to QGIS and have recently been testing the software as my org plan to migrate from ArcGIS to QGIS in the near future. I often use the 'Add values', 'Remove values' and the group values feature in ArcGIS (Symobology¬ Categories¬ Unique Values) to display Symbology for selected values. Is it possible to do this in QGIS? With categorized selected in the style menu the only options I can see is to add or delete. Any assistance would be greatly appreciated! Thanks


how can I assign an id column in sql query that merges geometries

I have a layer that consists of a number of features. Each features is a line segment and has an attribute named zone_name which takes one of 7 values. What I want to do is to use the atlas generation to print out a map for each zone (which consists of multiple features/lines) but (assuming I understand things correctly) the atlas generator iterates over all features in the layer -- first time I tried I got a large number of very small maps ;) After a bit of research I came up with this:               SELECT zone_name,          ST_Multi(ST_Collect(f.the_geom)) as singlegeom          FROM (SELECT layer_name, (ST_Dump(the_geom)).geom As the_geom                     FROM                     tracks_master, layers where feature_gid = gid ) As f     GROUP BY zone_name      in the database manager SQL window. When I run the query it produces a table with a single feature for each value of zone_name. The only problem is that there is no integer id field so I can not load this as a layer. I.e. when I select "load as layer" there isn't anything valid to put in "column with unique integer values". There must be a simple way of doing this but my SQL isn't up to it.


Reading non-spatial Access tables with ogr2ogr

I am trying to read a whole Access database into Postgis using the ogr2ogr (Version 1.9.1) command. For some reason however, it only reads spatial tables and no non-spatial tables. All my attempts to overcome this problem so far have failed^. I have tried configuring options such as "-nlt None" or "SET PG_LIST_ALL_TABLES=YES" How is it possible to read ALL tables and not only the spatial tables? This is the command I've used: ogr2ogr --config PG_LIST_ALL_TABLES YES -f "PostgreSQL" PG:"host=localhost user=postgres dbname=ele_schnottwil password=bsb port=5432" C:\ele.mdb -a_srs EPSG:21781 -overwrite -skipfailures I know PG_LIST_ALL_TABLES is a Postgis driver command but I read this could be used for the ODBC driver, too. Can anyone see the problem or help me out. Thanks so much in advance. Best regards, Arne


Enrich a dataset of POIs with OpenStreetMap

I have a list of POIs, some with a full description and some with only a few data entries, like the following:               6.9441000   50.9242000  [50677] (Ital) Casa di Biase [Köln]     6.9373600   50.9291800  [50674] (Ital) Al Setaccio [Köln]      However, I need the full dataset. Can I get this somewhere? If I have all the position data, is it possible to find the rest?   * a. name of the street   * b. name of the town So for example, the data should finally look like this:               10.5346100  52.1613600  [38300] (Chin) Wanbao Kommissstr.9 [Wolfenbüttel]     13.2832500  52.4422600  [14167] (Ital) LaPergola  Unter den Eichen 84d [Berlin]     13.3177700  52.5062900  [10625] (Chin) Good Friends  Kantstr.30 [Berlin]      Can I do this with OpenStreetMap? Should I parse OpenStreetMap data? Or OpenBabel?


How can I have my map open in the center of my visitors IP address?

I have a map of all restaurants in Florida. I would like to have the center of the map be the address of the visitors IP. How can I do this?


How to join attributes by location with command line

I would like to joint the attributes of 2 shapefiles by location but with a command line (linux). This feature is easy in QGIS but it is VERY slow. Is there a way of doing this with ogr or pyqgis ? Thank you


how to disable duplicate items in the legend for ArcGIS JavaScript API

I have a problem with duplicate items in the legend panel! ![Legend](http://i.stack.imgur.com/fvR2U.png) Can anyone help me? This is my feature Layer with popup info:               // Feature Layer             var msFeatureLayer = new esri.layers.FeatureLayer("url",{                 mode: esri.layers.FeatureLayer.MODE_SELECTION,                 outFields: ["*"],                 infoTemplate:new esri.InfoTemplate({                     title:"infos",                     content:"Code:${Code}
"                   }),                 displayOnPan:true,                 visible:true,                 opacity:1.0             });      And here i adding the layer to the map and to the legend:               legendLayers.push({layer:msFeatureLayer,title:'Ms'});             map.addLayers([msFeatureLayer]);                  // control when add a Layer to the map             dojo.connect(map,'onLayersAddResult',function(results){                 // Activate all Features from a Feature Layer                 showFromFeatureLayer(msFeatureLayer,"1=1");                           // Legend init                 new esri.dijit.Legend({                     map:map,                     layerInfos:legendLayers                 },"legendDiv").startup();                      // Add Checkboxes with label                 dojo.forEach(legendLayers, function(layer){                              var checkBox = new dijit.form.CheckBox({                         name: "checkBox" + layer.layer.id,                         value: layer.layer.id,                         checked: layer.layer.visible,                         onChange: function(evt) {                 var clayer = map.getLayer(this.value);clayer.setVisibility(!clayer.visible);                            this.checked = clayer.visible;                         }                     });                     dojo.place(checkBox.domNode,dojo.byId("toggle"),"after");                     dojo.place("
",dojo.create('label',{'for':checkBox.name, innerHTML:layer.title},checkBox.domNode,"after"),"after");                 });             });      thanks for your help !


Appending using PyShp

I am looking to do a simple append using Pyshp and keep running into walls. The editor object could not read the dbf header for the input shapefile so I tried using reader and writer as suggested but have run into some problems with writer. Initially I wanted to copy entire ShapeRecord objects but did not see a way to use writer to write them, so I tried a more granular approach like iterating and copying shape and record attributes which failed as well:               reader = shapefile.Reader('c:/temp/input.shp') #Open input file     w = shapefile.Writer(shapefile.POLYGON)        #Create writer     w.field('FLD1','C','20')                       #Create fields     w.field('FLD2','C','20')     w.field('FLD3','C','20')     w.field('FLD4','C','20')     w.field('FLD5','N','20')     w.field('FLD6','N','20')     w.field('FLD7','C','20')     for sr in reader.shapeRecords():         #iterate through shapeRecords in the input         w.poly(sr.shape)                     #write shape to writer         w.record(sr.record)                  #write record to writer     w.save('c:/temp/appendoutput.shp')       #save writer to file      I end up with two errors:   1. IndexError: tuple index out of range for `w.record(sr.record)` File "C:\Python26\ArcGIS10.0\lib\shapefile.py", line 959, in record [record.append(recordList[i]) for i in range(fieldCount)]   This perplexes me since print sr.record shows me that there are seven items in the record list that was read in by the reader and I have seven fields.   2. TypeError: iteration over non-sequence for the `w.poly(sr.shape)` Thanks again for any recommendations :)


Group and export to pdf based on an attribute field using DDP - Python Script

Using Python, I want to export a map series. Each series needs to be **grouped based on an attribute** field which is not being used in DDP settings (page num, page index etc...). For example: my shapefile RoadName.shp has a field called **RoadType**. Now when I ask my DDP to export maps, it should **group diff RoadTypes** into one single (multipage) pdf. Any quick solution please?


The best way to edit PostGIS/Geoserver data via web

I've got PostGIS installed and would like to create a web application to do some simple editions:   * Draw manually (over a Basemap) points, lines and polylines, using snapping points (points, end of lines, intersections, center of circles)    * Delete points, lines and polylines    * Rotate a line or group of lines around a point    * Move a line or group of lines with a base point    * Move polylines and points (manually and to a point given by the coordinates)   * Edit position of vertices of polylines and lines (manually and to a point given by the coordinates)    * Create buffers around points, lines and polylines  Which is the best way to do this via web? Is there any already-built JavaScript library? Thank you very much.


QGIS, SpatiaLite and GPS integration

I have QGIS an, access database and a Haicom gps all connected and running on a windows 8 machine. I want to integrate all 3 items so that while in the field i can enter information into the database ie bird names, numbers etc AND automatically update XY location from the GPS ie a site AND see where i am in QGIS relative to other layers such as vegetation aerials etc. status - at the moment QGIS easily logs a route and shows me where I am. my question is - should i migrate access data to spatial database such as SpatiaLite or similar and then work on automatic logging of GPS data with a field entry or write a macro/module in the access database to GET XY data from the GPS


Correct way to extract pixel data using GDAL. Current method gives repeated but changing value.

I am trying to open a spectral file and move bands of data from GDAL to OpenCV's Mat data structures. As of now I am just trying to make sure that I can correctly access each band using GDAL and pull the pixel data. So I want to explain what I think the following code is doing and if I am wrong(must be) then please correct me. It opens the file and loads the dataset. Then I get raster band 1 assigned to currentBand. The float pointer is then made to point at enough memory for one line of pixel data. RasterIO gets this data and places it in the newly allocated space. I then traverse the array printing the value that has been placed in each element. (Code in c++)               GDALDataset  *dataset;     GDALRasterBand *currentBand;     int blockXSize, blockYSize;     GDALAllRegister();     dataset = (GDALDataset *) GDALOpen( specImag->filename.c_str(), GA_ReadOnly );          currentBand = dataset->GetRasterBand(1);     currentBand->GetBlockSize(&blockXSize, &blockYSize);          float *fScanline;     fScanline = (float *) CPLMalloc(sizeof(float)*specImag->widthx);          currentBand->RasterIO( GF_Read, 0, 0, specImag->widthx, 1,                               fScanline, specImag->widthx, 1, GDT_Float32, 0, 0 );               for(int i=0; i<specImag->widthx; i++){             std::cout << "Data Point ";             printf("%d", i);             std::cout << " is: ";             printf("%d", fScanline[i]);             std::cout << std::endl;         }      (the code is not complete and will not run as is) When I run the code I get the same number (ex. 421076 or 4210432) for every element in the array. If I wait a little bit and run the program again I get a different number than the prevous run but still the same in each element. What am i missing when it comes to accessing the pixel/data in each band?


geoprocessing service keeps running with default value for choice list, instead of actual value selected

I have published a script tool to a geoprocessing service. The tool is supposed to take user input from a choice list (which contains values from a field), and then that is loaded into select by attribute tool to query the data. When I run the service, it keeps using the default value (which is the value I chose when I ran the script before publishing), even when I choose another value. Any idea why this might be? For my parameters, I have made the input feature class a constant value as well as the input field, and have a string parameter with the choice list. Then an 'any value' parameter that is derived from string parameter/choice list. Here's the code:               import arcpy, arcgisscripting, os     MapIndexLyr = 'MapIndexSelected'     ImageOutput = 'Image'     MI_FC_field = 'filename'     MapIndexSelect = r'in_memory\Selected_Feature'          InputFC = arcpy.GetParameterAsText(0)     QueryField = arcpy.GetParameterAsText(1)     MapName = arcpy.GetParameterAsText(3)          arcpy.MakeFeatureLayer_management(InputFC,MapIndexLyr)          arcpy.AddMessage('selecting features in Map Index layer...')     arcpy.SelectLayerByAttribute_management(MapIndexLyr,'NEW_SELECTION', QueryField+'=' +"'"+MapName+"'")          for row in arcpy.da.SearchCursor(MapIndexSelect,MI_FC_field):         selected_val = row[0]         img_name = selected_val.split('.')[0]         arcpy.AddMessage('iterating over selected Map Index features. ..')         arcpy.AddMessage('looking for image...')         for dirpath, dirnames, filenames in os.walk(r'D:\Records\GIS\MapSheetImages'):             for file in filenames:                 if file.endswith('.tif') or file.endswith('.sid'):                     if img_name in file:                         Out = arcpy.MakeRasterLayer_management(os.path.join(dirpath,file))                         arcpy.SetParameterAsText(4,Out)         del row, selected_val, img_name


ArcGIS Export to PDF with Thin White Lines

When I go to export to pdf from File -> Export Map, I see thin white lines that show up when I open up the map in Adobe. The lines appear only in the northern hemisphere in Russia, Canada and the United States. These lines do not appear when it is saved as a .ai file or .jpg, only when saved as a pdf. When this pdf is opened in Adobe Illustrator, those thin white lines split the polygons. Has anyone faced this problem before?


What gps app to use for mapping bus routes?

I'm looking for a GPS app for android to capture waypoints and lines. I'm going to map the route buses and the stop points. Which GPS app is recommended for marking “points” and marking the “lines”? I'll be using ArcGIS.


AGS Tiled Map Services and the .jpg file extension using js api

I have an old app that is still using esri js api 1.4 (i know, it's dated) but that app works and it can pull tiles from a 10.0 tiled map services. When it does that it calls for the tiles at the url with the \000\111\222 type format. It also appends a '.jpg' to the end of the request url. That works fine - i can even knock off the '.jpg' and still get the tile image. But now, my tiled service is being updated to 10.2 and I'm having some issues. When the request goes out with the .jpg I am no longer able to get the image - i get a 400 invalid url. If i lop off the .jpg then i get the image. The problem is that i'm not the one adding the .jpg to the request. That is the old api interacting with the ags service to produce that call. It's a black box to me. I load the layer and the js api does it's thing. Am i missing something or is my js api just too old to keep working?


Using Python to shorten polylines of varying lengths by 44% of their original length (or by a value contained in a field)

I have a feature class containing tens of thousands of polylines of varying lengths that were constructed from only a start and end point (of note, they are all at different angles, and many lines overlap one another). That said, I also have the original point file if it would be easier to work with. So....I am looking for a way to shorten these lines by 44% of their original length. Perhaps by creating a new point 44% "down the line" from the original starting point, then maybe removing the original end point and then redrawing the line? Seems like that would work. Or, suggestions welcomed, please. Of note,the data are projected, albers_NAD83. These are just straight lines, ie with a start point and an end point. The points also contain additional attributes that I would need to retain. Any guidance would be greatly appreciated. Thanks


How to use Leaflet to efficiently show parts of a massive dataset?

I keep seeing non-GIS web developers running into this problem, and I'm not sure what the solution is.   1. There is some dataset of thousands of items.   2. We want to show a map to the user, with the visible subset of them shown as interactive, clickable elements. What methods are there for doing this? I can think of these, but they're not very satisfactory, so I'm wondering what else there is:   1. Store all the data in a GeoJSON file, transfer it to the browser, and let Leaflet display it. Problem: doesn't really work with large datasets. TopoJSON raises the limit a bit. It also causes a big delay at page load.   2. Use Mapbox, store all the data in an interactive layer on Mapbox, and use Mapbox.js to display it. Works great, but costs money, and you can't host it yourself.   3. Use GeoServer to access a PostGIS database, use the WFS-geojson leaflet plugin to access the data from there. It probably works, but the WFS-geojson Leaflet plugin doesn't seem to be maintained anymore.   4. Use CartoDB, store all the data in a CartoDB table, and use CartoDB.js to display it. Works great, but can get very expensive. It's possible to host it yourself, but installing CartoDB is non-trivial. All of this makes me think there must be some much better, free way that I'm missing. What is it?


ArcMap redraws every layer when changing just one

I have 10 layers or so loaded in a data frame in ArcMap 10.2. Some are complex such as NHD hydrology and Tiger roads. Why does ArcMap redraw all layers when I change the symbology of just one? Is there a way to stop ArcMap from doing this behavior?


Increasing speed of Python scripts using arcpy

**Update 4/11/2014** It looks like the script was getting hung up on the Delete Features tool, so I switched to Truncate Table, as suggested in the answer below. I also removed the unused variables from the append tool. **Update 4/10/2014** I ran this script on my co-worker's computer (his machine has more memory AND contains ArcGIS 10.0/Python26) and it ran quickly. Hooray! Once my tech support finds the ArcGIS 10.0 CD, I will install and test to see if that improves the speed on my machine. To be clear, we are running the same script, our network drive and database connection are mapped identically, and the print statements are the same. I'll post an update here once that happens. **End updates** I need to increase the speed of some Python scripts that perform updates on an Oracle database. I had these Python scripts running well for a year+, via scheduled tasks and batch files to initiate the scripts. Last week I moved from an XP to a Windows 7 machine and ArcGIS 10.0 --> 10.1. Since then the scripts have become terribly slow. If I run this script using a small feature class (containing ~20 features) it runs in 30 seconds. If I use a medium feature class (~80,000 records) it runs in 15 minutes. The feature class I really need to be able to transfer quickly contains about 1,000,000 records - the script only goes as far as the print statement to check if the files exist (if statement in code below). This process would take only 35 minutes to complete on my XP/ArcGIS 10.0 machine. Below is the simplified code I have been testing with. Does anyone have suggestions on what I can do to increase the speed? Thanks, Patty               import arcpy, os, sys     from arcpy import env     arcpy.env.overwriteOutput = True     from datetime import datetime     import smtplib     import string     import urllib          #Define variables     inWorkspace = "O:/LANDING_PAD/BOE/example.gdb"      lpFeatures = inWorkspace + os.sep + "fc1"     outWorkspace =  "Database Connections\\THIS.sde"     arcpy.env.workspace = outWorkspace     workspace = ""     copyFC = outWorkspace + os.sep + "SDE.fc1_1" #The feature class the script will update via delete and append     schema_type = "NO_TEST"     fieldMappings = ""     subtype = ""     t = datetime.now()     print "This script began at: " + str(t)          if arcpy.Exists(lpFeatures) is True and arcpy.Exists(copyFC) is True:         print "Both files exist. Beginning delete..."         arcpy.DeleteFeatures_management(copyFC) #(copyFC)              print "ALL DONE DELETING!"              arcpy.Append_management(lpFeatures, copyFC, schema_type, fieldMappings, subtype) #Append data from landing pad to db         print "ALL DONE APPENDING!"         record_count = arcpy.GetCount_management(lpFeatures)         print record_count         r = datetime.now()         print "This script ended at: " + str(r)


Controlling page count increment of Data Driven Pages using ArcPy?

I was wondering if there is a way to add coding to Data Driven Pages with python. I would like to be able to take the "Data Driven Page with Count" and increment the count by one or in some cases two. For example, I have 10 data driven pages and need to add two more pages to the series i.e. making the final count 12 pages.


Conversion of point grid to a raster with open source tools (Gdal, QGIS ...)

I'm looking to do something which seem simple to me, but I haven't found a solution since 2 days I looked for. I've got a Postgis DB with lot of points (WGS84) spaced of 100 meters which cover a territory. For each of these points, I've got a value associated. I want to convert all theses point to a raster Geotiff. The ideal process will make a Geotiff file where each pixels corresponds to the value of the corresponding points of the db. I tried the heatmap plug-in of QGIS, but It doesn't do what I want, it interpolate lot of new points I don't need, I want in the tiff file only values I've got in my database. Any suggestion about how to easily generate this file ?


Prediction Population

I want to build a continuous data set with the population for 1970, 1981, 1991, 2001, 2011. I already have the population values for the same polygon, at sub-sections level (the lowest), for 1991, 2001 and 2011. The use of georreferenced databases only started (in Portugal) on the 1991 census so it's impossible to have that level of data for 1970 and 1981. However, I do have the number of buildings for each subsection for all the decades of the XX century. The questions are: How can I predict/estimate the population values at a sub- section level for 1970 and 1981? If possible, are the number of buildings helpful for the estimation?


ArcGIS Javascript API: Get Map Object from DIV

So I am currently writing an application using the ArcGIS JavaScript API, and I want to access the map object that is held inside of my div from a js module I have. I was thinking it would be a simple JQuery call but it doesn't seem to be the case, as when I try to use `map.centerAndZoom()`, it states that it is not a member of the object that returned by that call. Is there anyway to be able to reference the map object that was created in my main index.html elsewhere?


Cannot properly import MODIS MOD07_L2 .hdf files into R using rgdal

I've been trying to accurately import and process some .hdf files from the MODIS Atmospheric Profile product (MOD07_L2) in R for several days now. Still, there's something going wrong during data import. The below code can be reproduced using one exemplary .hdf file (MOD07_L2.A2013001.0835.006.2013001192145.hdf) that can be downloaded from Dropbox.               library(rgdal)          # Extraction of metadata via `GDALinfo`         filename <- "MOD07_L2.A2013001.0835.006.2013001192145.hdf"     gdalinfo <- GDALinfo(filename, returnScaleOffset = FALSE)     metadata <- attr(gdalinfo, "subdsmdata")          # Extraction of SDS string for parameter 'Skin_Temperature' (formerly 'Surface_Temperature')         sds <- metadata[grep("Skin_Temperature", metadata)[1]]     sds <- sapply(strsplit(sds, "="), "[[", 2)          # Raster import via `readGDAL`        sds.rg <- readGDAL(sds)      So far, so good, but here comes the confusing part:               > summary(sds.rg$band1)        Min. 1st Qu.  Median    Mean 3rd Qu.    Max.    NA's       -14870  -14850  -14850  -14840  -14840  -14820   53529       Considering the fact that Skin_Temperature has an official valid range from 150 to 350 K (see MOD07_L2:Format & Content), the mean would inherit a value of               > (-14840 - (-15000)) * 0.01     [1] 1.6      after considering the corresponding add_offset (-15000) and scale_factor (0.01). Note that we're still talking about Kelvin, not °C. Extracting SDS No. 8, i.e. Skin_Temperature, using               library(gdalUtils)     gdal_translate(filename, dst_dataset = "tmp.tif", sd_index = 8)      and opening the resulting file called "tmp.tif" in QGIS resulted in seemingly reliable values centered around 15000, i.e. roughly 27 °C. However, importing "tmp.tif" back into R using `raster` again resulted in values comparable to the ones shown above:               > summary(raster("tmp.tif"))                   tmp     Min.    -14867.31     1st Qu. -14848.13     Median  -14845.89     3rd Qu. -14840.53     Max.    -14819.93     NA's         0.00      I've been searching the internet and stumbled across similar problems related to rgdal. However, when I tried to cast `toUnSigned` on band 1 of my previously generated 'SpatialGridDataFrame', I received the following error message:               > toUnSigned(sds.rg$band1, 16)     Error in toUnSigned(sds.rg$band1, 16) : band not integer      Apparently, the data imported into R is not even of type integer (what it is supposed to be), but numeric:               > sds.rg$band1[1:5]     [1]        NA        NA -14839.40 -14840.25 -14839.26      Is there an apparent mistake in my code, or is there any point I miss when importing the .hdf and .tif files using `rgdal`? I would be extremely grateful for any kind of help.


Is it possible to create shapefiles in ArcGIS and convert into SpatialNET?

I'm looking at the prospect of working for a client who has implemented a SpatialNET solution. I have never used SpatialNET and primarily have experience with ESRI and Autodesk products. The client and Synchronoss (who owns SpatialNET) have both told me that the only solution is to purchase a SpatialNET license and VPN into their database. When I asked them about creating files in AutoCAD or ArcMap, they said it would not work. I find it hard to believe that the SpatialNET system is completely incompatible with them, especially considering that SpatialNET appears to be based off of AutoCAD. Anyone here have experience with this product?


Exporting well data from IHS Kingdom in a format that plays nice with ArcGIS

A coworker of mine has a bunch of well data (including georeferencing and attribute tables) in IHS Kingdom. We need to convert this data into a format that works with ArcGIS (e.g. .shp files). Unfortunately, it looks like the only options for exporting are as follows: .wtv2 - Tobin Well .asc - Landmark Well 2.x/3.x .asc - Landmark Well OWX .wel - Charisma .gds - GDSII .dat - Well Production Data As far as I can tell, none of these are compatible, and each of them are extremely poorly documented online. As is Kingdom IHS, for that matter. The only solution we've managed so far would be to export the map as a large PDF before georectifying and manually digitizing each of the wells in ArcGIS. However, there are too many wells to make this a viable option.


How to set up a QGIS project in MVS 2008?

I have built QGIS using Visual Studio 2008. I would now like to start using the QGIS API for Windows Development. Of course I have Qt installed on Studio I hope correctly. I am new to Visual Studio and can't find any "How To" guides on setting up a QGIS/Qt project. I am trying to go with the tutorial "Write a simple mapping application using the new QGIS Canvas API". I have created a new project called QGisSample and added the code to main.cpp. When I attempt to build the project I get all sorts of errors. This has to be a configuration problem setting up the project. But I just do not know where to start. There are, a few of, the errors:                   1>Compiling...         1>main.cpp         1>D:\QGIS.1.8.0-2\Source\src\core\qgis.h(31) : error C2470: 'QGis' : looks like a function definition, but there is no parameter list; skipping apparent body         1>D:\QGIS.1.8.0-2\Source\src\core\qgis.h(160) : error C2146: syntax error : missing ';' before identifier 'GEOPROJ4'         1>D:\QGIS.1.8.0-2\Source\src\core\qgis.h(160) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>D:\QGIS.1.8.0-2\Source\src\core\qgis.h(168) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>D:\QGIS.1.8.0-2\Source\src\core\qgis.h(168) : error C2146: syntax error : missing ';' before identifier 'GEO_EPSG_CRS_AUTHID'         1>D:\QGIS.1.8.0-2\Source\src\core\qgis.h(168) : error C2371: 'QString' : redefinition; different basic types         1>        D:\QGIS.1.8.0-2\Source\src\core\qgis.h(160) : see declaration of 'QString'         1>D:\QGIS.1.8.0-2\Source\src\core\qgis.h(168) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>D:\QGIS.1.8.0-2\Source\src\core\qgis.h(168) : error C2440: 'initializing' : cannot convert from 'const char [10]' to 'int'         1>        There is no context in which this conversion is possible         1>D:\QGIS.1.8.0-2\Source\src\core\qgsapplication.h(30) : error C2470: 'QgsApplication' : looks like a function definition, but there is no parameter list; skipping apparent body         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C2146: syntax error : missing ';' before identifier 'tr'         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C2433: 'QFile::QString' : 'inline' not permitted on data declarations         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : warning C4183: 'tr': missing return type; assumed to be a member function returning 'int'         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C2146: syntax error : missing ';' before identifier 'trUtf8'         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C2433: 'QFile::QString' : 'inline' not permitted on data declarations         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : warning C4183: 'trUtf8': missing return type; assumed to be a member function returning 'int'         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C2146: syntax error : missing ';' before identifier 'tr'         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C2433: 'QFile::QString' : 'inline' not permitted on data declarations         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : warning C4183: 'tr': missing return type; assumed to be a member function returning 'int'         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C2146: syntax error : missing ';' before identifier 'trUtf8'         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C2433: 'QFile::QString' : 'inline' not permitted on data declarations         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int         1>d:\qt.os.4.8.2\include\qtcore\../../src/corelib/io/qfile.h(68) : error C4430: missing type specifier - int assumed. Note: C++ does not support default-int      There is plenty more where that came from. I set this project up as a Qt Application. I also added the library where the core .h files are and the QGIS lib files that were generated when I built QGIS. I'm just fumbling around with this trying to Google the various problems. For now it appears I have left something out that is fundamental to the QGIS project. Any suggestions are gladly appreciated.


How to define image map with pins in Javascript?

I have image map in .png format. All I need to do is to add a sidebar where users click on for example "restaurants" and pins with restaurants will show on an image map. When user mouseover a pin, it will show some more informations. I would like to do it in Javascript, but not with Google Maps. I have really small amount of pins for four different categories. What is the best way to do it in Javascript?


MFC Dialog Application using ESRI ArcEngine 10.0

I'm new to ESRI ArcObjects and I'm trying to create a basic MFC dialog application. I add the License and MapControl to the dialog resource. I compile and run the app, but I get the following runtime error:- Could not bind to a valid ArcGIS installation. CoCreateInstance of OLE control {D7EAD71E-B663-4DCE-82C7-AEE2793DF5F3} failed. > > > Result code: 0x80040111 Is the control is properly registered? Warning: > Resource items and Win32 Z-order lists are out of sync. Tab order may be not > defined well. MFCDlg2.exe has triggered a breakpoint I have registered the OCX's and I can see them in the registry. I'm running Win7 (64 bit) and Visual Studio 2008. Please can anyone help Thanks.


How do you use GRASS's v.kernel?

I am flummoxed on how to use GRASS's v.kernel. I have a vector layer of around 2.5 million points. I want to make a heat map using v.kernel to show concentrations, since I have variable instances with overlapping points, sometimes huge overlaps. I've already gotten this vector layer in GRASS, and it displays just fine. I've tried using GRASS's v.kernel command based on what I've seen here and on other forums, and I can't get it to do anything besides output a raster that's just a pink square. Here's the command I'm using:               v.kernel --verbose input=master_grass7 output=master_grass7a_heatmap stddeviation=.0001      I've varied the **stddeviation** to all sorts of values from 1000000 to .000001, and it had no effect. I've read the v.kernel documentation repeatedly and don't really understand what it's getting at. At least, the instructions are on esoteric concepts, nothing practical. I've also checked the source code, and I'm not really understanding it, either. Yes, I can read C. The problem is it depends on a lot of stuff defined elsewhere in GRASS GIS. I've also done a lot of Google searching, and I can't find a comprehensive guide. All that I'm getting are scattered copies of the v.kernel doc/man page or people who apparently got it to work without a fuss. I've also checked up on the concept of kernel density estimation (KDE), and even then I don't see how to use the v.kernel command. That command appears to be a specific interpretation of KDE; its switches don't appear to correspond well to generic KDE concepts. So back to the main question here: how can someone who is not intimate with GRASS product development use the v.kernel command? Is there a plain language translation available?


In Qgis, Is it possible to edit the datasource of a Postgis layer?

After creating a layer with RT_SQL_Layer plugin, I want to change some details in the query. Is it possible? or I have to aways create a new layer every time I want to change something in the query? More specifically, I want to change the "Source for This Layer" as shown in example bellow: ![enter image description here](http://i.stack.imgur.com/AL0e9.jpg)


Rotating geometry to specific angle using pyqgis?

I do not know how to rotate a geometry using pyqgis. Is there any sample/example code available to show me how?


Problem: Arcpy Graph function empty graph output

**Background:** I have some python code that I need help with parameter passing into the MakeGraph_management arcpy function. I would like to create a two-series vertical line graph for every pair of input point feature classes using the 'zip' function (ie. input fgdb1/fc1 and fgdb2/fc1 >> MakeGraph >> graph1.png; input fgdb1/fc2 and fgdb2/fc2 >> MakeGraph >> graph2.png, ...). **My Problem:** Now I am trying to create graphs in batch using ArcPy/python 2.7.3, but have been struggling with the way these ArcPy functions work around the template file. I created a graph template in ArcGIS 10.2 with everything setup (axis labels, titles, footers, legend names, etc.) I just need to get the data input working to correspond with Series1 and Series2. The code below successfully iterates and creates graphs for every pair of input datasets as planned. However, the graphs contain no plotted data. I thought I could just use MakeGraph_management and SaveGraph_management since I have a template instead of also working with the ArcGIS Graph class. ![mock transect graph, no data](http://i.stack.imgur.com/Pmwpt.png) **My Code:**               # Import system modules     import arcpy     from arcpy import env          # enable overwrite outputs     arcpy.env.overwriteOutput = True          # input feature directory: yr0_transects     dir_yr0_transects = r"W:\Working\abcdefg_working\PointDistance\data\point_distance\year0\output\output_year0.gdb"          # input feature directory: yr1_transects     dir_yr1_transects = r"W:\Working\abcdefg_working\PointDistance\data\point_distance\year1\output\output_year1.gdb"          # output directory     outWorkspace = r"W:\Working\abcdefg_working\Profiles"          # Create list of year 0 Transect features     arcpy.env.workspace = dir_yr0_transects     yr0_transects_list = sorted(arcpy.ListFeatureClasses())          # Create list of year 1 Transect features     arcpy.env.workspace = dir_yr1_transects     yr1_transects_list = sorted(arcpy.ListFeatureClasses())          # Set workspace     arcpy.env.workspace = r"W:\Working\abcdefg_working\Profiles"          # Set local variables     graph_tee = r"W:\Working\abcdefg_working\Profiles\Profile_graph_template.tee"                    for yr0_transect, yr1_transect in zip(yr0_transects_list, yr1_transects_list):         try:             #Cast to string because ListFC() returns unicode.             string_yr0_transect = str(yr0_transect)             string_yr1_transect = str(yr1_transect)                  # create title for graph             graph_title = "Profiles, B.F.E." + string_yr0_transect[0:17] + string_yr0_transect[47:]                       # Create the graph object             graph = arcpy.Graph()                  # Add a vertical bar series to the graph             graph.addSeriesLineVertical(yr0_transect)             graph.addSeriesLineVertical(yr1_transect)                  # Specify the title of the Graph             graph.graphPropsGeneral.title = graph_title                  # Create graph output name             output_graph_name = string_yr0_transect[0:17] + string_yr0_transect[47:] + "_profile" + ".png"                  # Output a vertical line graph, which is created in-memory             arcpy.MakeGraph_management(graph_tee, graph, "temp_graph")                  #Execute SaveGraph             mygraph = arcpy.SaveGraph_management("temp_graph", output_graph_name, "MAINTAIN_ASPECT_RATIO", "1024", "768")                  del mygraph              except:             print arcpy.GetMessages()


Displaying a box with text on a leaflet map

Is there a way to display a label or textbox on a specific geo-place on a map produced by leaflet.js.? I know about the Marker- and Polygon Method But I can't find a way to put text in the polygon. I want the text and box to be the same size on every zoom level. And when dragging the map the box should be stuck to the map; like a sticky note. Is there a way to achieve this with Leafletjs or with SVG?


How to calculate Distance, Azimuth and Dip from two XYZ coordinates

**I have two XYZ coordinates and I would like to calculate the Distance, Azimuth and Dip between them.** I have come across a similar question on this site that calculates it the other way around (with a given starting coordinate and a given Distance, Azimuth and Dip and it calculates an end coordinate.) I have been trying to figure it out and I am sure its very simple. Here is that other question page. It says that this is a conversion between a spherical and Cartesian coordinates. How to convert Distance, Azimuth, Dip to XYZ?


What are some 3D visualisation tools that successfully work with GIS data?

There are a large number of 3d modelling tools available (both open and closed source). I would like to be able to:   * Import GIS data into a 3d environment and retain georeferenceing.    * Retain the quality and resolution of the source data.   * View the data in real time with game like speed. I have tried a number of tools but they never quite seem to be able to do these three points successfully. Do any tools exist that can reliably satisfy these requirements?


Difference spatial interpolation and trend surface analysis

What is the difference between spatial interpolation and trend surface analysis? For example, I read a scientific study that applied a neural network to estimate a surface solely based on the x and y coordinates, claiming that they do spatial inteprolation. However, to my understanding this is trend surface analysis, because Tobler's law is not considered. Am I right?! Is the difference between these two that spatial interpolation takes neighboring observations into account?


Can I connect  (& edit) to a postgresql database with arcgis without arcSDE?

After trying in vein to get spatialkit and zgis working with my postgresql database in arcgis 9.3 and 10 , I was wondering if anyone has had any luck with these products. My postgresql database has been serving geospatial data for 2 years through geoserver and can be edited in quantum gis fine....so why not in ArcGIS? Does anyone have any other recipes for using postgresql data with arcgis but WITHOUT ArcSDE? yours, Rob


How to export polygon coordinates?

How can I export each polygon-vertice's coordinate (every point coordinate of every polygon) of a polygon Geodatabase-feature class to an ascii text file (csv or txt)? I don't want to convert it to a point FC first. I cannot use ET Geowizards! I don't have FME (or is the little FME encluded in ARC able to do this? How?) I don't want to export to XML format - except if someone can tell me how to convert XML to a simple text file.


Error when looping Polygon to Raster with python through different fields ArcGIS 10

This is the first time I have posted a question on here so please have some patience! I've been trying to create a python script to loop through numerous shapefiles converting various fields to raster outputs. I've been working off a lot of what has been said in this previous thread on the same issue -> Python loop Polygon to raster not working But can still not get it to work. I keep getting the dreaded 999999 error. The code I've been trying is shown below:               #import arcpy module     import arcpy     from arcpy import env          #set working environment     env.workspace = "C:\\Users\\s175850\\Documents\\GIS\\GIS_Data\\Source_Methodology\\Pesticides\\Application_Term\\2008"     Dir = env.workspace          #List feature classes     fclist = arcpy.ListFeatureClasses()          #loop     for fc in fclist:            output = Dir + "\\r" + fc + "_ras"       cell = 5            #process polygon to raster       arcpy.PolygonToRaster_conversion(fc, "January", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "February", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "March", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "April", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "May", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "June", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "July", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "August", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "September", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "October", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "November", output, "CELL_CENTER", "", cell)       arcpy.PolygonToRaster_conversion(fc, "December", output, "CELL_CENTER", "", cell)          print "finished polygon to raster"      I'm not sure if there would be issues with outlining the 12 polygon to raster conversions for the 12 fields like I have done so I tried it with just one polygon to raster run for one field name. I still managed to get the same error. I am using ArcGIS 10. Thank you in advance for your help.


Some text-sizes of labels are not working in TileMill

I am using TileMill to create some offline tiles. Currently I have problem that some of the text sizes are not working and they are displayed as weird signs. Do you have any idea what could be wrong?                 #waterways {       line-smooth: 1;       line-join: round;       line-cap: round;              //label styling and placement       text-name: [name];       text-face-name: 'Calibri Italic';       text-fill: @water_labels;       text-placement: line;        text-halo-fill: fadeout(white, 0%);       text-halo-radius: 1.8;              [type = 'river'] {line-color: @water; line-width:12;          [zoom=15] {text-size: 15;}         [zoom=16] {text-size: 16;}         [zoom=17] {text-size: 17;}         [zoom=18] {text-size: 18;}         [zoom=19] {text-size: 19;}       }     }      ` Picture a) text-size: 18 ![text-size: 18;](http://i.stack.imgur.com/CE3pj.png) Picture b) text-size: 17 ![text-size: 17;](http://i.stack.imgur.com/GFoFb.png)


Looking for a simple hatching algorithm

If I have an arbitrary polygon, and wish to hatch it with lines at a given angle and spacing, is there a simple low-level algorithm which will return the set of hatch lines inside the polygon?


Quality of free geodata compared to commercial sources

What are the sources of excellent geo-location data as possessed by Google Maps ? I am wondering that, even with large amount of volunteer work involved, why is OSM data not comparable to theirs in terms of quality ?


How do I get the Coordinates of a shapefile into .csv

I have followed all the instructions on this site, but my output always shows Area, and Perimeter of each polygon. Am I missing something? Should I be changing a setting somewhere to select coordinates rather than area/perimeter?


Arcgis basemap with proxy.ashx not loading correctly

I have a .Net application using the JavaScript API to draw a map. I use the proxy.ashx file to proxy requests to secure map services (dynamic tokens), configured in proxy.config. I added an entry into proxy.config for ArcGIS Online:-               <serverUrl url="http://services.arcgisonline.com/" matchAll="true" />      My JavaScript is as follows:-                   esri.config.defaults.io.alwaysUseProxy = true;         esri.config.defaults.io.proxyUrl = "/proxy.ashx";              map = new esri.Map("map", {             extent: new Extent({                 xmin:174.5, ymin:-39, xmax:174.7, ymax:-38.8,                 spatialReference:{wkid:4326}             }),             basemap: "topo"         });      If I exclude the basemap, my map draws as expected but if I include the basemap (as above), the basemap is drawn partially and I get an infinite loop of null requests occurring:                   GET /proxy.ashx?undefined/tile/NaN/NaN/NaN?blankTile=false HTTP/1.1      I have stripped out all other code to eliminate that as a cause. If I comment out the two lines of JavaScript to use the proxy the basemap works perfectly so the problem has something to do with my proxy.config and/or the proxy.ashx file. Possibly I have an older version old proxy.ashx. Has anyone seen this error before, or point me in the right direction to get a basemap working with the proxy configuration? I am using JavaScript API 3.8


Django application "GeoNode" gives message "DLL load failed" on windows

I'm trying to install GeoNode which is a Django application on windows through these simple steps on this page. I've reached a point where is the GeoNode is installed and I've created my postgres database successfully. The problem is when I try to test the installation via the command `paver test` or sync my database using the command `python manage.py syncdb` I get this message `ImportError: DLL Load Failed: The specified module could not be found.`   My Environment is as follows: > Windows 7 32bit   >  Python 2.7.5   >  Django 1.5.5   >  GeoNode 2.0   >  Postgresql 9.1.3   >  Postgis 2.0


ArcGIS 10.2 Query Layer on SQL Server performance

I'm using a Query Layer on SQL Server in ArcMap. The Query Layer executes instantly in in SQL Server but but takes so long to draw in ArcMap that the system appears unresponsive for about 10 minutes or longer. During the ArcMap draw the one of the CPU's is maxed out on the SQL Server process. My Query is the STIntersects of a buffer on a line feature (Shannon) against a polygon feature class (Townlands) , as follows;               SELECT TOWNLANDS.TL_ID,TOWNLANDS.Shape FROM dbo.TOWNLANDS as townlands     with(index(FDO_Shape))      JOIN dbo.Shannon on townlands.Shape.STIntersects      (Shannon.Shape.STBuffer(2.0))=1      The query returns 186 rows instantly. These can be draw in the SQL Server Management Studio Spatial pane without a problem When I build a Query Layer in ArcMap with exactly the same syntax the system becomes unresponsive but does draw eventually. It looks as though, perhaps, ArcMap s not using the spatial index or is doing so differenlt from SQL Server causing an inefficient query on SQL Server that takes an age to return. Can anyone advise on a remedy? Thanks               ArcGIS Desktop: 10.2     ArcSDE: 10.2     RDBMS: Database and version: SQL Server 2008     OS: Windows Server


Measuring polyline lengths running through a raster - is this possible?

I have a reclassified raster layer with 4 values and a vector layer of polylines. I would like to discover the length of the lines running through each raster value, or the percentage of the lines within each raster value. Can anyone tell me if this is possible please? I was informed Zonal Statistics may be able to do this, but to no avail. Thanks.


How to get National Geographic basemap in QGIS?

I'm proficient in ESRI's products and recently started using QGIS. I've used the national geographic world map ( **NatGeo_World_Map** ) in ArcMap, but am struggling to replicate this in QGIS. The map is also available **here**.


Limit what fields can be edited GeoExplorer OpenGeoSuite

I have an issue with limiting what fields can be edited by the user. If I link my layer to my PostGIS table, the user is able to edit all the fields for a specific point. -If I create a view I can limit what they can see, however views are read only? Any obvious solutions I am not thinking about to allow me to limit which fields they can change or fill in when modifying a layer? Using OpenGeo Suite, GeoExplorer. Thanks


geoext popup + geoserver content.ftl

I am trying to include a jpg to a geoext popup. There are various shapefiles, but only one needs to load pictures. With only openlayers it was easy to do, I just used the content.ftl from geoserver 2.4. However, we switched to geoext and sadly I am not able to figure out how to do it now. Can somebody give me a hint how to force GeoExt.Popup to use the different content.ftl of geoserver - or is there another way?


Polygon to raster without losing any polygons

I'm trying to convert polygons to raster. But the output raster dataset doesn't have polygons which went into as input data. I need all the date in the attribute to follow with to the raster dataset as well. What am I doing wrong? I'm using ARCGIS 10.1.


Visualise data belonging to an area

I am looking for a way in QGIS (1.8.0) to visualise data, which is related to an area (polygon on a shape file). There are no x-y-coordinates with the data, just the info of the belonging area, means at an arbitrary position inside the area. Example: Samples x1,x2 and x3 are located in area a1, y1 to y5 in area a2. As the samples have certain properties on which I intend to do queries on, distinct locations seems somehow necessary (e.g. sample x1 and x3 have a certain property and therefore should be shown/highlighted ). I initially thought of assigning each sample some random coordinates from inside the areas, but what a hassle.. So, how to approach that??


What are the best ways to connect two coordinates on mapinfo?

I have with the help of this forum drawn a coordinate system with longitude and latitude at different location. I am trying to connect these sites what is the best way to connect. Is it just drawing lines or something else?