title: Plotting a Route on Google Maps using a string passed to a function as a parameter

text: I'm trying to plot a route, using Google Maps. I need to be able to pass a string to the function, to build the Polyline. I can do the following and it works fine:               var flightPlanCoordinates = [            new google.maps.LatLng(18.9173333333, 117.2265),            new google.maps.LatLng(18.9178333333, 117.226833333),            new google.maps.LatLng(18.9203333333, 117.228333333),            new google.maps.LatLng(18.9328333333, 117.235833333)         ];      If I bring this into a function, I never see the route. Javascript shows no errors. My code looks like this:                       var Route;         var map;              function initialize() {             var lLat = 18.9173333333;             var lLong = 117.2265;             var iniZoom=13;             var myLatLng = new google.maps.LatLng(lLat, lLong);             var mapOptions = {                 zoom: iniZoom, center: myLatLng,                 mapTypeId: google.maps.MapTypeId.ROADMAP, panControl: false,                 zoomControl: true, scaleControl: true,                 streetViewControl: false, overviewMapControl: false             }             map = new google.maps.Map(document.getElementById("map_canvas"), mapOptions);             lMapRoute='18.9173333333,117.2265,18.9178333333,117.226833333,18.9203333333,117.228333333,18.9328333333,117.235833333';                  PlotRoute(lMapRoute);         }              function PlotRoute(pRoute) {             var RouteOptions = {                 geodesic: true,                 strokeColor: '#FF0000',                 strokeOpacity: 1.0,                 strokeWeight: 2             }             Route = new google.maps.Polyline(RouteOptions);             Route.setMap(map);             var RouteTaken=[];             Route.setPath(RouteTaken);                  var RoutePoints = pRoute.split(",");             for (var i = 0; i < RoutePoints.length-2; i=i+2) {                 RouteTaken.push( new google.maps.LatLng( Number(RoutePoints[i, 0]), Number(RoutePoints[i, 1]) ) );             }         }      Can anyone see why this would not work. Many thanks in advance. Colin


title: How to convert Shapefile geometries to WKB using OGR?

text: I got the geometry of a feature of the shapefile, and i want to save that geometry in the postgis (in WKB formats like it happens when we import shapefiles using shp2pgsql and psql commands). How do I convert then To get Geometry i have used OSGeo OGR library eg:               feat = layer.GetFeature(0)     geometry = feat.GetGeometryRef()      and I have got               <osgeo.ogr.Geometry; proxy of <Swig Object of type 'OGRGeometryShadow *' at 0x0096A2D8> >      so how will convert it into the WKB geometry? I am using Python for this.


title: Open MapGuide MrSID rasters with FDO GDAL. invalid extent. Not drawing

text: When I try and save the resource I get: Spatial context extent appears to be invalid (or Default). The XML:               <?xml version="1.0"?>      OSGeo.Gdal.3.4 Default This coordinate system has been overridden. LOCAL_CS[" _XY-MT_ ",LOCAL_DATUM[" _X-Y_ ",10000],UNIT["Meter", 1],AXIS["X",EAST],AXIS["Y",NORTH]] LOCAL_CS[" _XY-MT_ ",LOCAL_DATUM[" _X-Y_ ",10000],UNIT["Meter", 1],AXIS["X",EAST],AXIS["Y",NORTH]] Dynamic -10000000 -10000000 10000000 10000000 0.001 0.001 # These sids are 205mb a piece and theirs 53 of them (sat images). I decoded them to geotiff and they came out 2+gig ea. What do I need to do to get these to display in MG? Whats the smartest method to get them web ready?


title: Where can I find the latest WMX (JTX) Object Model Diagrams for ArcGIS 10?

text: I'm looking for a version that includes ESRI.ArcGIS.JTXUI--I'm hoping someone knows of an updated and online version of `C:\Program Files\WMX\Developer Kit\Diagrams\ObjectModel.pdf`. For reference, I'm on ArcGIS 10 SP2:   * ObjectModel.pdf is dated: 10/27/2010    * SDK.chm is dated 6/4/2010 _So these weren't updated by either service pack._


title: ArcObjects C#: Create Random Points fails with ERROR 999999: Error executing function

text: I am having trouble running some ArcObjects C# code that calls create random points. I would not be doing this in ArcObjects, but I am suffering with the same issue using Arcpy as seen here: http://forums.arcgis.com/threads/12782-Python-scripting-and-creating-random- points?p=85866&viewfull=1#post85866 Any help as to the cause of this error would be greatly appreciated (Indeed, an official confirmation that create random points is bugged would also be appreciated). code excerpt:               ESRI.ArcGIS.RuntimeManager.Bind(ESRI.ArcGIS.ProductCode.Desktop);     //ESRI License Initializer generated code.     m_AOLicenseInitializer.InitializeApplication(     new esriLicenseProductCode[] {      esriLicenseProductCode.esriLicenseProductCodeArcInfo },      new esriLicenseExtensionCode[] { esriLicenseExtensionCode.esriLicenseExtensionCodeSpatialAnalyst });     //ESRI License Initializer generated code.          //For some reason the license thing above doesn't always do the job so this is here     IAoInitialize ao = new AoInitialize();     ao.Initialize(esriLicenseProductCode.esriLicenseProductCodeArcInfo);          Geoprocessor gp = new Geoprocessor();     gp.OverwriteOutput = true;     gp.LogHistory = true;          //Input parameters     string constrainingFeatures = "lineRoute";     string outFeatures = "test";     string workspace = "\\\\1.2.3.4\\output\\Routes.gdb";          //Open workspace     IWorkspaceFactory workspaceFactory = new FileGDBWorkspaceFactoryClass();     IWorkspace gdb;     try     {         gdb = workspaceFactory.OpenFromFile(workspace, 0);     }     catch     {         Console.WriteLine("Unable to open geodatabase");         return;     }     IFeatureWorkspace fws = (IFeatureWorkspace)gdb;          //Get constraining feature class     IFeatureClass fc = fws.OpenFeatureClass(constrainingFeatures);          //Create random points - DOES NOT WORK WITH A PATH TO WORKSPACE OR THE IFeatureWorkspace     //CreateRandomPoints crp = new CreateRandomPoints(fws, outFeatures);     CreateRandomPoints crp = new CreateRandomPoints(workspace, outFeatures);     crp.number_of_points_or_field = "2000";     crp.minimum_allowed_distance = "1 Meter";          //DOES NOT WORK WITH A PATH TO cfc OR the IFeatureClass     //crp.constraining_feature_class = fc;     crp.constraining_feature_class = workspace + "\\" + constrainingFeatures;          //this out_feature_class does not seem to make a difference     crp.out_feature_class = workspace + "\\" + outFeatures;               try     {         gp.Execute(crp, null);     }      Here is the log:               <ResultViews>     <ResultView Tool="Create Random Points">     <CommandLine>CreateRandomPoints_management \\1.2.3.4\Output\Routes.gdb \\1.2.3.4\Output\Routes.gdb\test \\1.2.3.4\Output\Routes.gdb\lineRoute '0 0 250 250' 2000 '1 Meters' POINT 0</CommandLine>     <ToolSource>C:\Program Files (x86)\ArcGIS\Desktop10.0\ArcToolbox\Toolboxes\Data Management Tools.tbx\Feature Class\CreateRandomPoints</ToolSource>     <StartTime>Wed Mar 16 14:31:01 2011</StartTime>     <Parameters>     <Inputs>     <Parameter Label="Output Location" Type="Dataset">\\1.2.3.4\Output\Routes.gdb</Parameter>     <Parameter Label="Output Point Feature Class" Type="Scalar">\\1.2.3.4\Output\Routes.gdb\test</Parameter>     <Parameter Label="Constraining Feature Class" Type="Dataset">\\1.2.3.4\Output\Routes.gdb\lineRoute</Parameter>     <Parameter Label="Constraining Extent" Type="Scalar">0 0 250 250</Parameter>     <Parameter Label="Number of Points [value or field]" Type="Scalar">2000</Parameter>     <Parameter Label="Minimum Allowed Distance [value or field]" Type="Scalar">1 Meters</Parameter>     <Parameter Label="Create Multipoint Output" Type="Scalar">false</Parameter>     <Parameter Label="Maximum Number of Points per Multipoint" Type="Scalar">0</Parameter>     </Inputs>     <Outputs>     <Parameter Label="Output Feature Class" Type="Dataset">\\1.2.3.4\Output\Routes.gdb\\\1.2.3.4\Output\Routes.gdb\test</Parameter>     </Outputs>     </Parameters>     <Environments>     <Environment Label="Precision For New Coverages">SINGLE</Environment>     <Environment Label="Auto Commit">1000</Environment>     <Environment Label="Minimize memory use during analysis on terrains">false</Environment>     <Environment Label="Compression">LZ77</Environment>     <Environment Label="Coincident Points">MEAN</Environment>     <Environment Label="Random number generator">0 ACM599</Environment>     <Environment Label="Raster Statistics">STATISTICS 1 1</Environment>     <Environment Label="Level Of Comparison Between Projection Files">NONE</Environment>     <Environment Label="Output has Z Values">Same As Input</Environment>     <Environment Label="Maintain fully qualified field names">true</Environment>     <Environment Label="Tile Size">128 128</Environment>     <Environment Label="Pyramid">PYRAMIDS -1 NEAREST DEFAULT 75</Environment>     <Environment Label="Default TIN storage version">CURRENT</Environment>     <Environment Label="Output Spatial Grid 1">0</Environment>     <Environment Label="Cell Size">MAXOF</Environment>     <Environment Label="Output has M Values">Same As Input</Environment>     <Environment Label="Output Spatial Grid 2">0</Environment>     <Environment Label="Output Spatial Grid 3">0</Environment>     <Environment Label="Maintain Spatial Index">false</Environment>     <Environment Label="Precision For Derived Coverages">HIGHEST</Environment>     </Environments>     <Errors>     <Error>ERROR 999999: Error executing function.</Error>     <Error>Failed to execute (CreateRandomPoints).</Error>     </Errors>     <EndTime>36.00 seconds)</EndTime>     </ResultView>     </ResultViews>      The only other output message is:               Error HRESULT E_FAIL has been returned from a call to a COM component.      I also posted this here: http://forums.arcgis.com/threads/26029-Create-Random- Points-fails-with-ERROR-999999-Error-executing-function. Thanks!


title: Freelancers, how do you handle deliveries to remote clients?

text: Freelancers and small business owners, how do you handle product delivery and support to your remote clients? Example: I have a client who knows enough about ArcMap to add a shapefile, turn layers on and off, print a map, etc. - GIS isn't their thing, they uses it to make maps and manage digital data storage. When I have a deliverable, I spend alot time on the phone troubleshooting the deployment. Sure I can use GoToMeeting, get onto the remote box, and install everything, but its $50 a month, man! I'm also thinking of setting up a FTP site, which will be easy enough and will help with delivery of the product. But are there good options to GoToMeeting that anyone here can provide feedback on? EDIT: Lots of good stuff being shared here, but please keep in mind I need something _incredibly simple_ for my client to install and work with, before they can hand control over to me.


title: Best practices for identifying training sites

text: Prior to image classification, collection of training sites and obtaining a signature file is the most critical step. Field visits with a reliable GPS Unit is usually the first step followed by overlaying the points on the to-be- classified image. Other than this when identifying training sites on satellite images, especially the coarse resolution ones what are the rules normally followed by professionals - like minimum number of pixels in a training site, distance between two training sites for same/different classes, how improvements are made after preliminary classification etc. Standard references are also welcome.


title: Create WFS with Mapguide Studio Enterprise 2011

text: i need create a WFS with MGS 2011, i used the help users mapguide for create the WFS but that didnt work, can anyone help me please? i am new in this tecnology.


title: Geoserver installed successfully but WFS doesn't work

text: I installed GeoServer 2.1 RC1. It installed successfully, but when i click on WFS in services category, i got the following message. Oops, something went wrong... Sorry, something unexpected happened on the server. Here's an error report you can include in a JIRA bug report about this issue: org.apache.wicket.WicketRuntimeException: Can't instantiate page using constructor public org.geoserver.wfs.web.WFSAdminPage() at org.apache.wicket.session.DefaultPageFactory.createPage(DefaultPageFactory.java:212) at org.apache.wicket.session.DefaultPageFactory.newPage(DefaultPageFactory.java:57) at org.apache.wicket.request.target.component.BookmarkablePageRequestTarget.newPage(BookmarkablePageRequestTarget.java:298) at org.apache.wicket.request.target.component.BookmarkablePageRequestTarget.getPage(BookmarkablePageRequestTarget.java:320) at org.apache.wicket.request.target.component.BookmarkablePageRequestTarget.processEvents(BookmarkablePageRequestTarget.java:234) at org.apache.wicket.request.AbstractRequestCycleProcessor.processEvents(AbstractRequestCycleProcessor.java:92) at org.apache.wicket.RequestCycle.processEventsAndRespond(RequestCycle.java:1250) What is wrong ? How can i fix this problem ?


title: How to make MapInfo able to edit geographic data from PostGIS?

text: I'm using MapInfo 10.5 to read geographic data from PostGIS 1.5.2. After a pg_dump/pg_restore to another PostGIS-database on another server I'm trying the same via a VPN-tunnel. When I open a table I get the message "Map objects cannot be edited" and then I can't edit the geometries. The other attribute data is editable. I've tried to create my database and tables both by hand and with easyloader with the same result. At the same time when I restore the database on my local computer I only get the message for one of the tables. Any idea what could causing this problem? Regards Ulrik


title: How can i stop rendering Qgis

text: In my python plugin ,i select features from a vector layer dynamically using a loop and executing : VectorLayer.select(feature.id()) when i start the loop, the selection become very slow because the rendering . can i suspend and reactivate rendering in qgis programatically in python.


title: Why am I getting "sequence size must match size of the row" error?

text: I'm trying to use an Insert Cursor to insert a point geometry into an existing point feature class. I'm finding the coordinate and appending them to a list called rowValues, ex:               rowValues.append("((({0}, {1}), 34))".format(X, Y))      This is giving me a list of value in this format:               ((3656320.51246, 1613747.99484), 34)      Here is the cursor I am using:               cursor = arcpy.da.InsertCursor("Points",("SHAPE@XY","QuarterSections"))      When I run this code:               for row in cursor:         cursor.insertRow(row)      I get this follow error:               Traceback (most recent call last):       File "Quarter.py", line 110, in <module>         cursor.insertRow(row)     TypeError: sequence size must match size of the row      Can I not use the string formatting? I've been trying to figure this out for awhile now, but I'm stumped.


title: Does Geoserver support WMTS protocol?

text: Does Geoserver support WMTS protocol? I find only WMS WFS and WFS-T support.


title: I have a problem to display the legend in QGIS Web Client

text: I have installed QGis Map Server and QGis Web Client on Centos 6.4 64-bit operating system, everything works fine except the legend box which is filled with garbage data instead of text (it seems that some fonts are missing), please have a look of the screenshot attached. OS: CentOS release 6.4 (Final) Kernel: 2.6.32-358.2.1.el6.x86_64 #1 SMP Software packages: qgis-mapserver-1.8.0-9_0.el6.x86_64; qgis-1.8.0-9_0.el6.x86_64; mod_fcgid-2.3.7-1.el6.x86_64; fcgi-2.4.0-10.el6.x86_64; httpd-2.2.15-28.el6.centos.x86_64; postgresql92-libs-9.2.4-1PGDG.rhel6.x86_64; postgresql92-9.2.4-1PGDG.rhel6.x86_64; postgresql92-server-9.2.4-1PGDG.rhel6.x86_64; ![enter image description here](http://i.stack.imgur.com/KewzI.jpg)


title: Open --> Setting --> Project Properties will shut down Qgis

text: When I want change my project and I select setting and the project properties . It's will give me some message, like "Exception - Bad allocation" How do fix the problem ? http://imgur.com/XEHcs


title: Polygon inside polygon invalid in QGIS

text: Looking at this map data from naturalearth (part of Oman & UAE to be specific), each country is defined as a multipolygon. In this case UAE (in red) contains two polygons and Oman (beige) three. ![Polygon inside polygon](http://i.stack.imgur.com/jZGAY.jpg) Problem is that QGIS claims UAE to be invalid, since the smaller polygon is "inside" the other one. I've tried making the small polygon a ring instead, but that's also not okay, since it does not touch the other ring. How should this be organized in order to be valid? Must these be separate features? Joined physically?


title: Scaling CartoDB to multiple machines

text: I'm interested to scale CartoDB installation to multiple machines to distribute computational resources. We start using torque with some big tables (more than 5 to 10M points) and the processing for each datacube tile is quite long. Is it possible to have some hints on what you are currently using to scale CartoDB and in particular PostgreSQL/PostGIS?


title: how to georeference my png using gdal

text: Hello I have created my own fantasy map, I have the map in a full size png and the map broken down into tiles ( with no geo data ) I installed MS4W and tried to use it to get the georeference data put is seems the program got stuck while trying to do the 4th step here this is my info I got from the 3rd step: -Size is 4961, 3231  -Coordinate System is `'   * Image Structure Metadata: -INTERLEAVE=PIXEL    * Corner Coordinates:    * 1st corner Upper Left( 0.0, 0.0 )   * 2nd corner Lower Left ( 0.0, 3231.0 )    * 3rd corner Upper Right ( 4961.0, 0.0 )    * 4th corner Lower Right ( 4961.0, 3231.0 )   * Center ( 2480.5, 1615.5 )    * Band 1 Block=4961x1 Type=UInt16, ColorInterp=Red    * Band 2 Block=4961x1 Type=UInt16, ColorInterp=Green   * Band 3 Block=4961x1 Type=UInt16, ColorInterp=Blue    * Band 4 Block=4961x1 Type=UInt16, ColorInterp=Alpha 4th step said something about -gcp not available or something I was just wondering if someone could get my png georeferenced an into a kml or csv file for me so I can continue my project or help me with using the right code to continue on here is my map png It would be most appreciated and thanked


title: Where can I find QGIS Cartogram plug-in?

text: Can anyone point me in the right direction of the cartogram plugin please? I've had a look on the net but can't find what I'm looking for. Thanks


title: Change symbology of layer of unique symbol to unique with python

text: I have a problem concerning how can I change the symbology of a layer with unique symbole to unique values. I don't know how to update it. I used this code:               import arcpy     mxd = arcpy.mapping.MapDocument("CURRENT")     lyr = arcpy.mapping.ListLayers(mxd,"region")[0]     if lyr.symbologyType == "UNIQUE_SYMBOLE":       lyr.symbology.valueField = "REGION"       lyr.symbology.addAllValues()     arcpy.RefreshActiveView()     arcpy.RefreshTOC()     del mxd      * * * it does not give me an error message and the version I'm working on arcgis 10.1 the problem is that I want to change the symbology that is a kind of a unique symble to category "uniques values"?


title: How to setup a SQL trigger in CartoDB

text: I need to setup 2 SQL triggers in CartoDB that will update my **inland** and **coastal** tables each time my **obs** table is updated. I have never setup a SQL trigger but I have looked at the documentation and it appears as though I need to have a function stored somewhere that the trigger activates. So far my first trigger looks like this:               CREATE TRIGGER update_inland     AFTER INSERT ON obs     FOR EACH ROW     EXECUTE PROCEDURE update_inland()      But I get an error message that " _function update_inland() does not exist_ ". I want my function to be this:               UPDATE inland SET lt_dispatch_level = obs.named_lt_dispatch_level      FROM obs      WHERE obs.created_at = (SELECT MAX(created_at) FROM obs)      AND inland.cartodb_id = 1      How can I create this function and subsequent trigger in CartoDB? I have been reading the documentation for both Postgres triggers and functions but I need some help.


title: How to import text only out of a string field to a new field

text: I want to import data out of a string field in my shapefile to a new field for labeling. The string field contains river names and codes and I just want to display the text i.e river names


title: Pass a variable from Modelbuilder into a Python script

text: I'm new to modelbuilder and python. I have used modelbuilder to iterate through all rows in a shapefile (Points), process them, and produce a seperate output shapefile for each row. The filename is generated using an inline variable (I think that is the correct terminology). i.e. all output files have identical filenames except that the value of the PointID from the original input file is appended to the filename to enable them to be distinguished. This works fine. I've also written a python script to take the above output files from the model one at a time and do additional processing. At present the input filename is hard coded into the script (i.e. I have to edit the script to change the input filename each time I run it). Again this works fine while just processing a single file. However, I want to integrate the model and the python script so that I can run the script for each row in the original feature class automatically. My question is therefore, how do I replace the hard coded name of the input file with a variable in the python script and how do I pass the filename to the script. Any assistance will be appreciated. Thanks Steve


title: Unattractive image artefact when exporting from QGIS's Print Composer

text: When I export my map from Print Composer in QGIS 1.7.4, an odd red band appears across the middle of the image. This happens when I export to both PNG and PDF. Its not visible on the screen before export, or present in the original raster image.![Images showing red band across the middle](http://i.stack.imgur.com/UuGnB.jpg) The band is to the left of the island halfway down the coast, below Pwani. Any ideas what's causing this and how to fix it?


title: Is it impossible to measure Right of Way in ArcGIS?

text: I'm helping bring CAD people into ArcGIS so that everyone throughout the company is on the same page and there is one database for everything. Right now, one of the biggest hurdles is making sure all of the functionality that they have with GIS is available and as easy to use as it is in what they're using now. I want this transition to be easy and with as little hesitation/bitterness as possible. I cannot find a measurement tool that snaps at a 90 degree angle. The minimum distance tool doesn't work, because it will find the minimum for the whole line segment, instead of taking in a piece of a line. I want to calculate the minimum at a specific point or in a very localized area. Also, I would prefer a button, similar to the measure tool. I don't want them to have to load data into the toolbox dialogs in order to do a simple measurement. The measure tool doesn't snap perpendicular. We're wanting to minimize the amount of possible error. I want to check this forum before I do one of two things--create a custom button or create a feature class that I can use as a temporary measurement-- they snap perpendicular to one side of the right of way line to the other side and then after they finish, they delete the line. Thank you for your help. *We are currently using 9.3, but are upgrading to 10 soon.


title: Point (of a linestring) within Polygon using ogr and Python

text: I'm currently working on a project in which I need to build a topological network out of the geometry features I find in shapefiles. So far using Ben Reilly's open source project I've managed to transform linestrings into networkx edges, as well as detecting close features (other linestrings say) and adding them to the closest point so that I can run shortest path algorithms. But that's fine for one shapefile. However, I now need to connect features from different shapefiles into a big networkx graph. So for instance, if a point is within a polygon I would connect it (by connect it I mean add a networkx edge - add_edge(g.GetPoint(1), g.GetPoint(2)) with the point in the next shapefile that also is within a polygon that shares a similar attribute (say, ID). Note that the polygons in the different shps only share the same IDs and not the coordinates. The points that fall within the polygons also do not share the same coordinates. My solution to this problem was to identify the point that resides in a polygon, store it, find the point in the next shapefile that resides in the polygon with the same id and then add the networkx edge between them. How to find if a point resides within a polygon? Well, there's a well known algorithm: RayCasting algorith that does that. This is where I actually got stuck though, because in order to implement the algorithm I need the coordinates of the polygon, and don't know how to access them right now, even after skimming through a documentation of the OGR's Geometry. So, the question that I'm asking is how to access the polygon points, or coordinates OR is there an easier way of detecting whether a point falls within a polygon? Using python with osgeo.ogr library I coded the following:                if g.GetGeometryType() == 3: #polygon                     c = g.GetDimension()                     x = g.GetPointCount()                     y = g.GetY()                     z = g.GetZ()      see the image for a better understanding of my issue. ![alt text](http://i.stack.imgur.com/hCKRq.png) **[EDIT]** So far I've tried storing all the polygon objects in a list with which I would then compare the **linestring** first and last point. But Paolo's example is related to using the Point Object reference and the Polygon Object reference, which would not work with a line object reference since not the **entire** line is within the polygon but rather the first or the last point of its linestring. **[EDIT3]** Creating a new Geometry point object from the coordinates of the first and last point of the linestring and then using that one to compare against the polygon geometry objects saved in a list seems to work just fine:               for findex in xrange(lyr.GetFeatureCount()):         f = lyr.GetFeature(findex)         flddata = getfieldinfo(lyr,f,fields)         g = f.geometry()         if g.GetGeometryType() == 2:             for j in xrange(g.GetPointCount()):                 if j == 0 or j == g.GetPointCount():                     point = ogr.Geometry(ogr.wkbPoint)                     point.AddPoint(g.Getx(j),g.GetY(j))                     if point.Within(Network.polygons[x][0].GetGeometryRef()):         print g.GetPoint(j)      Thanks Paolo and Chris for the hints.


title: Full list of ISO ALPHA-2 and ISO ALPHA-3 country codes

text: I'm searching for a **full list of countries with their country codes**. Like on this page (need full and valid): http://www.nationsonline.org/oneworld/country_code_list.htm Thanks.


title: Getting line FID from polygon ITopologyGraph

text: I want to export polygon Arc list to txt file for that I created topologygraph for personal geodatabase feature class (polygon) using ITopologyGraph::build. I got polygon FIDs but I am not getting IDs for individual edges of polygons. Please help me.


title: Using Standard and Maplex Labelling Engines at same time?

text: Is there any way to have one layer using the Maplex Labelling Engine for labeling and another layer using the Standard Labelling Engine at the same time in ArcMap? https://www.dropbox.com/sh/iwfuysjnx4n732x/AAC6LFDz07N-PfL2-1hJltXZa I have uploaded multiple screen shots. Mainly I am trying to be able to have all of the labels line up in the gorman shopping center. And I want to be able to see all the labels.


title: Best GeoDjango MySQL workaround for distance queries

text: I'm using mysql 5.1.50, so the distance QuerySet method is not available to me. Here is my model:               class Location(models.Model):         name = models.CharField(max_length='100')         point = models.PointField(spatial_index=False)              objects = models.GeoManager()      I would like to find all the locations where the point is within 5 miles of another point (eg. test_point). I've found some potential ways to do this but am unsure of all of them.               Location.objects.filter(point__within=test_point.buffer(1))      But it's not clear to me what the unit is for '1'. It looks like degrees. I can convert miles to degrees but it seems like there's a cleaner way to do this.


title: Python console not loading with QGIS1.7-Ubuntu 10.04.2 LTS lucid,

text: Any tips on how to get the Python console working on Q-GIS (Quantum GIS)? (Its not working from the software's menu.) System:   * Python 2.6.5    * QGIS1.7 OS   * Linux: Ubuntu 10.04.2 LTS lucid Error message from python running in the Terminal window:               Python 2.6.5 (r265:79063, Apr 16 2010, 13:09:56)          [GCC 4.4.3] on linux2         Type "help", "copyright", "credits" or "license" for more information.         >>> from qgis.core import *         Traceback (most recent call last):           File "<stdin>", line 1, in <module>         RuntimeError: the sip module implements API v7.0 to v7.1 but the qgis.core module requires API v6.0         >>>


title: Standard Colors for Landcover Maps

text: I have classified Landsat images into five classes - A variety of Palm, Water bodies, built-up/open area, other vegetation and miscellaneous. I used GRASS GIS for the process and the final classified maps use inappropriate colors. Now as I have to assign appropriate colors to these classes I was thinking about existence of internationally accepted standards for such purposes. As palm which covers only 10% of total area is the class of my interest(shown in yellow in the following screenshot) I want it to be clearly distinguishable in my A4 size paper output. In the most unlikely event of absence of standards proper color suggestions are greatly welcome. ![enter image description here](http://i.stack.imgur.com/F9bl3.png)


title: Calculate XY Coordinate with "CalculateField" tool

text: I used model builder to build a model which can calculate the XY coordinates and it's perfect! However, when I export it to VBscript and run it in VBA it crashes.... The error happens at                gp.CalculateField_management RasterT_Kriging1_shp__2_, "x", "pPoint.X", "VB", "Dim Output As Double"      and says object required. Following is my code to calculate the X:               Dim Output As Double          Dim pPoint As IPoint          Set pPoint = [Shape]          Output = pPoint.X      Can anyone help me? Ok, following is my code:               Dim gp As Object     Dim shape As Object     Set gp = CreateObject("esriGeoprocessing.GPDispatch.1")          ' Check out any necessary licenses     gp.CheckOutExtension "spatial"          ' Load required toolboxes...     gp.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx"     gp.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx"     gp.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx"               ' Local variables...     Temp = "C:\temp"     Kriging_poin4 = "C:\temp\New File Geodatabase.gdb\Kriging_poin4"     Output_variance_of_prediction_raster = ""     Dim Temp As String     Dim Kriging_poin1 As String     Dim Temp__2_ As String     Dim RasterT_Kriging1_shp__4_ As String     Dim RasterT_Kriging1_shp As String     Dim RasterT_Kriging1_shp__2_ As String     Dim RasterT_Kriging1_shp__6_ As String     Dim RasterT_Kriging1_shp__5_ As String     Dim point2_Events_shp As String     Temp = "C:\Temp"     Kriging_poin1 = "C:\Temp\Kriging_poin"     Output_variance_of_prediction_raster = ""     Temp__2_ = "C:\Temp"     RasterT_Kriging1_shp__4_ = "C:\Temp\RasterT_Kriging.shp"     RasterT_Kriging1_shp = "C:\Temp\RasterT_Kriging.shp"     RasterT_Kriging1_shp__2_ = "C:\Temp\RasterT_Kriging.shp"     RasterT_Kriging1_shp__6_ = "C:\Temp\RasterT_Kriging.shp"     RasterT_Kriging1_shp__5_ = "C:\Temp\RasterT_Kriging.shp"     point2_Events_shp = "C:\Temp\point2 Events.shp"          ' Process: Feature Class To Shapefile (multiple)...     gp.FeatureClassToShapefile_conversion "'C:\temp\New File Geodatabase.gdb\point2'", Temp__2_          ' Process: Kriging...      gp.Kriging_sa point2_Events_shp, "x", Kriging_poin1, "Spherical 249.101224", "250",   "VARIABLE 12", Output_variance_of_prediction_raster          ' Process: Raster to Point...     gp.RasterToPoint_conversion Kriging_poin1, RasterT_Kriging1_shp__4_, "Value"          ' Process: Add Field...      gp.AddField_management RasterT_Kriging1_shp__4_, "x", "DOUBLE", "12", "4", "", "", "NON_NULLABLE", "NON_REQUIRED", ""          ' Process: Add Field (2)...     gp.AddField_management RasterT_Kriging1_shp, "y", "DOUBLE", "12", "4", "", "",   "NON_NULLABLE", "NON_REQUIRED", ""          ' Process: Calculate Field...     gp.CalculateField_management RasterT_Kriging1_shp__2_, "x", "pPoint.X", "VB", "Dim Output As Double"                    Dim pPoint As IPoint          Set pPoint = [shape]          Output = pPoint.X               ' Process: Calculate Field (2)...     gp.CalculateField_management RasterT_Kriging1_shp__6_, "y", "pPoint.Y", "VB", "Dim    Output As Double"                    Output = pPoint.Y               End Sub      I am following the steps you told me last time in Using kriging with data in GDB but it's still doesn't work!


title: qgis print composer: more then one "overviews" to show overlap of more then two maps

text: I have a map composition with 3 pages (maps). I would like to draw the overlapping on every map. So the first map is overlapping the second, and the second is overlapped by the first an the third. To show the overlap, you can use the overview-feature but there you can just choose ONE map - the first OR the third. I need do choos both. Is that possible? Or is there a other way to draw the overlap?


title: How to use dbscan in python for gis coordinates?

text: I am curious how to input gis coordinates (lat long) in the http://scikit- learn.org version of dbscan for python. Thanks


title: Change "Default" label in label manager Arc 10

text: Does anyone know if there is a way to change the Default Label in the Label Manager for ArcGIS 10? Otherwise is there a way to mass label features. I have about 500 feature layers in this MXD and in each of the layer tables is a column named "Label". Because the locations were geocoded, the default uses the text string from the column "Loc_name". I'd like it to use my label column for all of them. In the label manager, I have the option to label everyone with the default, thats why i originally asked if that default can be changed. I've tryed to copy the text string from my column into "Loc_name" but the feild length is set at 14 and it is un-editable and 14 is not enough characters.


title: What annoys you about your GIS Environment?

text: Hello Spatially Enabled Folk, I have a question regarding your favourite GIS environment. Over the past few decades, there have been the usual major players, MapInfo, ESRI, Intergraph (Grass?), and we all love and hate aspects of the major players of GIS. The purpose of my question is to provide some forum (community wiki?) where we can place suggestions that better our working lives. So please be specific. Please try and be constructive, give descriptions of process, data, expected outcomes and scenarios. > **Constructively, What are your primary issues with your GIS package of > choice?**    _The Seed for the Question_   Looking at the Hadoop project by Apache and some of the impressive work completed there, which started the thinking cogs. Since GIS is such a computationally intensive profession, how has this technology not being incorporated in our favourite GIS software! There are some fantastic solutions that could enhance a GIS professionals life. The other purpose is to look at the possibility of taking some of these ideas and potentially incorporate them in FOSS4G projects (QGIS? GeoServer? New Project???), increase usability and functionality and give the big guys a run for their money ;)


title: How to enable MapTips in ArcGIS Engine? (ArcGIS 10)

text: My organization is working on a new version of our software using ArcGIS Engine. Previous versions of the software (running on ArcGIS Engine 9.3) displayed MapTips fine as long as they were enabled in the mxd. Our new version running on ArcGIS Engine 10 doesn't display MapTips (and they are fine in the mxd). **Is there any code or setting to enable MapTips, or any known reasons that would cause them to not display in an ArcGIS Engine 10 application?** I have already discovered that layers within a Basemap layer won't display MapTips, so that is not the cause. Thank you for any advice. **EDIT** We tried what iterationx suggested, but the MapTips wouldn't display regardless of the state of that flag. The suggestion brenth made was what we needed:               mapControl.ShowMapTips = true;      Thanks


title: Can I export a layer's field value to metadata ?

text: Is there a way around for me to export a layer that has the field name "Quadname" and export it to metadata to capture all the quads names? I am looking for that to work on my metadata using Arcgis 9.3.1. Thanks


title: Atributing fields in arcmap and excel

text: How do I populate a field in ArcMap (and Excel) with only a specific part of another field's attributes? For example: I want to populate the LATDEG field with only the _degree part_ of a latitudinal co-ordinate and not the entire value.


title: Converting WGS84 points to a 2D x,y point (plotting)

text: I have been trying to convert the WGS84 style geo-positions to x,y coordinates on a 2D plane. I have looked into many different links, websites, studies and I honestly gave up on finding a plane copy and paste solution. Things that I have observed:   * I select a point on the Openstreetmap let's say **48.8532, 2.3528**.   * I let that point to be presented as the origin on my 2D plane, (0,0) in x,y.   * I already know the ranges of latitude and longitude values.   * I define the range for the x,y values of the 2D plane. (ex. x=[0, 50000] y=[0,5000])   * I interpolate the geo-positions using linear interpolation to the numbers on the 2D plane. Now my questions:   1. Is this correct? Is the WGS84 representation a linear coordinate system?   2. If the answer to the previous question is no, then could you provide a code that does this correctly? (no plane formulas and no links to further studies please, just code in C, C++, C#, Java or Javascript. Note that I'm not asking for a distance, I'm asking for a point-to-point convertor code).


title: Using IRelQueryTable

text: I create a RelQueryTable like this to join a GeoFeatureLayer to an external table. Everything seems to work (no errors, I can look in the table and see the correct fields.) Silly question: what do I return from this that I can use as a GeoFeatureLayer? When I try to cast the table to GeoFeatureLayer it will not work. It was hard enough to riddle out how this even worked (and I still have many unanswered questions that the bathtub documentation doesn't help with while I'm adrift at sea.) Does anyone have intelligent samples for these APIs? (Instead of changing 9.x to 10.x they should have released some code.)               Type memoryRelationshipClassFactory = Type.GetTypeFromProgID             ( "esriGeodatabase.MemoryRelationshipClassFactory" );     IMemoryRelationshipClassFactory memoryRelationshipFactory =         Activator.CreateInstance( memoryRelationshipClassFactory ) as              IMemoryRelationshipClassFactory;          IRelationshipClass relationshipClass = memoryRelationshipFactory.Open(          "Join",     geoFeatureLayer.DisplayFeatureClass as IObjectClass,     layerField,     table as IObjectClass,     tableField,     "forward",     "backward",     esriRelCardinality.esriRelCardinalityOneToOne );          Type rqtFactoryType = Type.GetTypeFromProgID( "esriGeodatabase.RelQueryTableFactory" );     IRelQueryTableFactory rqtFactory = Activator.CreateInstance( rqtFactoryType ) as          IRelQueryTableFactory;          IRelQueryTable relQueryTable = rqtFactory.Open(         relationshipClass,         true,         null,         null,         EXTERNAL_FIELDS,         true,         false ) as IRelQueryTable;


title: can widget developed by ArcGIS API for flex as an add-in?

text: I am currently involved in a project needs me to develop a widget using ArcGIS API for flex. It will be integrated to an existing web application. I am confused now because I haven't touch the existing project and code, but only start a new project from scratch. I want to know whether the flex widget can be developed to an add-in which like the desktop add-in develop in .net, python or java? What should I start for developing an extra widget for a web application using ArcGIS API for flex? I am using arcgis server 10.2.2 and adobe flash builder 4.6. very very new in flex and arcgis api for it, totally no idea.


title: Dynamic defining of transformation type for ISpatialReferenceFactory.CreateGeoTransformation call

text: Suppose given a point (IPoint) in some coordinate system. Required to transform it to another coordinate system. If both of coordinate systems are predefined, we can use `IGeometry2.ProjectEx` with transformation specified by corresponded constant from `esriSRGeoTransformationType`, `esriSRGeoTransformation2Type`, or `esriSRGeoTransformation3Type` enumerations. The code looks like (source)               package arcgissamples.geometry;          import com.esri.arcgis.geometry.*;     import com.esri.arcgis.system.*;          public class ProjectGeometry {            public static void main(String[] args) {         try {           EngineInitializer.initializeEngine();           AoInitialize aoInit = new AoInitialize();                // Create a point with Geographic coordinates...           Point point = new Point();           point.putCoords(-100.0, 40.0);                System.out.println("");           System.out.println("Original coordinates: " + point.getX() + ","                     + point.getY());                // Create the SpatialReferenceEnvironment...           SpatialReferenceEnvironment spatialReferenceEnvironment = new SpatialReferenceEnvironment();                // Apply the initial spatial reference...           ISpatialReference geographicCoordinateSystem = spatialReferenceEnvironment                     .createGeographicCoordinateSystem(esriSRGeoCSType.esriSRGeoCS_NAD1927);           point.setSpatialReferenceByRef(geographicCoordinateSystem);                // Create the output projected coordinate system...           ISpatialReference projectedCoordinateSystem = spatialReferenceEnvironment                     .createProjectedCoordinateSystem(esriSRProjCSType.esriSRProjCS_NAD1983UTM_13N);                // Create the GeoTransformation...           //***************************************************           // !!!! How to define arg for createGeoTransformation??           IGeoTransformation iGeoTransformation = (IGeoTransformation) spatialReferenceEnvironment                     .createGeoTransformation(esriSRGeoTransformationType.esriSRGeoTransformation_NAD1927_To_WGS1984_5);                // Project the point...           point.projectEx(projectedCoordinateSystem, esriTransformDirection.esriTransformForward,                     iGeoTransformation, false, 0.0, 0.0);           System.out.println("Projected coordinates: " + point.getX() + " , "                     + point.getY());                System.out.println("Done!");                aoInit.shutdown();              } catch (Exception ex) {           ex.printStackTrace();         }       }     }      But problem is that both original coordinate system and output coordinate system are defining at the time of program execution. So I need a function to define constant for `ISpatialReferenceFactory.CreateGeoTransformation` method like               int defineTransformationType(ISpatialReference srcCoordSystem, ISpatialReference destCoordSys)      Is it possible? And how to implement this function?


title: Shortcut to copy current cursor coordinates in ArcMap

text: Is there a shortcut or an extension to ArcMap 10 that would allow me to quickly copy (to clipboard) the coordinates at the current map position/click? I could implement such a functionality easily enough, however I'm hoping someone has already done so. edit: Apparently the code is much simpler than I figured:               IPoint pt = ArcMap.Document.CurrentLocation;     Clipboard.SetData(DataFormats.StringFormat, pt.X + " " + pt.Y);


Where can I find a shapefile of zoning for the city of Phoenix and other municipalities? I would rather not have to geo-reference a picture and digitize over it, there must be some data already out there somewhere!


title: Importing OpenStreetMap (OSM) file into MySQL using Osmosis?

text: I've been trying to importing osm file into mysql using osmosis, but still failed.


title: How to pass distances to Multiple Ring Buffer in ModelBuilder?

text: I want to use Multiple Ring Buffer in ModelBuilder however I'm not able to pass it 'distances' values. It only accepts one block with a single double value. It's not accepting a list of double values or a string with separated values. How to pass multiple 'distances' values to Multiple Ring buffer script? EDIT for PolyGeo:      The model is really trivial: ![enter image description here](http://i.stack.imgur.com/MVjMe.png) As long as the distances block is a single double value I'm able to connect it to the script, if it's something else, it doesn't allow me to connect it as Distances. The script itself uses convertValueTableToList function on the entered value. However I can't find an appropriate block with ValueTable type.


title: QGIS: Turning Filter On and Off

text: When a filter is applied to a layer the layer cannot be edited. Is there any way of quickly turning the filter on or off without entering the Query Builder dialog box and removing the filter?


title: Trouble implementing Key authentication in geoserver

text: I'm trying to implement key authentication in geoserver 2.3.3.      The steps I followed are:     * Added the authkey plugin   * Created the authkey filter   * Secured the layers by assigning workspaces to different users   * Added authkey filter to the filter chain I'm passing the **authkey** (obtained from authkey.properties file) as a **parameter** along with the request while accessing the layers from the client side of my application. But still the authentication does'nt work and the layer's are inaccessible.   Is there anything that I'm missing out or anything wrong with my procedure.   Kindly help me with this. Thanks.


title: Login / Logout Geoserver Geoext

text: I am quite new to Sencha and I am trying to create a WEB GIS including GeoExt and Geoserver! In Geoserver I have created users and roles that;s why I need to create Login Logout but I can't find a solution, only this enter link description here but this works only when a page loads! So I need a button in a toolbar, or somewhere else in a panel input fields and button to work! Can you please help me!


title: How to programmatically associate a layer with a custom UI?

text: I know it is possible to associate a layer with a custom UI for editing feature attributes using the QGIS interface. But is this possible with Python? This link shows how to achieve the result manually. But I would like to be able to create layers automatically and have them use custom UI without having to redefine the ui file and method each time.


title: Setting correct units for Near tool?

text: I need to use the near tool to calculate distances between data points and a stretch of coastline. My problem is that when I run the function, the values that come up in the Field of the points' attribute table are in unknown units. When performing the function, the Units option is set to meters, but through comparison using the Measure tool I know the returns in the attribute table are not in meters. The distances are apparently measured in the units of the input features' (the datapoints) co-ordinate system. The co-ord system of the points is GCS_WGS_1984. Can anyone tell me what the units of this geographic co-ord system are, or how I can find out? Thanks a lot


title: How do I add coordinates to a polygon layer to be able to sort the polygons from north to south in QGIS?

text: I have a polygon file, where each polygon represents one property. I want to be able to sort the polygons from north to south or vice versa. I need them to be sorted in this way when I export them to a different program. In the other program I want the properties to appear in a north to south order when I visit them one by one. My first problem is, the polygons don't have information about there coordinates in the attribute table. How can I add this? I have tried using the field calculator, but I don't seem to get it right. And, next step, if I sort them in QGIS, can I somehow "save" that order so that they appear in the same order in the other program? Thank you for any help! EDIT:   I have tried the method in the answer below, but I don't get it to work on my data. Then after almost giving up, I read in the help of the field calculator that $y gives the y coordinate of a point. But I have a polygon layer? What do I do?


title: automatic save outputs variables in the folder workspace using pyhon for arc map 10.1

text: I want to create python script with aim to import him in environment arcmap 10.1 and to make a hydrologic analysis because the scipt meybe can have many outputs raster file or feature class would like to find a way for automatic save outputs variables in the folder workspace for example I want to import a raster file (Dem) and to take as export a 3 raster file from the command fill sinks, Fill Sinks,flow directions and flow accumulation from arc hydro toolbox cl = my import dem fill = output for Fill Sinks fdr = output for flow directions Fac = output for flow accumulation and i want export them(3 raster) automatically saved to the folder workspace without definitive output parameters from the user               import arcpy          arcpy.ImportToolbox("C:\Program Files (x86)\ArcGIS\Desktop10.1\ArcToolbox\Toolboxes\Arc Hydro Tools.tbx", "archydrotools")          # Local variables:     arcpy.CheckOutExtension("spatial")     arcpy.env.overwriteOutput = True     arcpy.env.workspace = arcpy.GetParameterAsText(0)     cl = arcpy.GetParameterAsText(1)     # Process: Fill Sinks     arcpy.FillSinks_archydrotools(cl, fill, "", "", "ISSINK_YES")     # Process: Flow Direction     arcpy.FlowDirection_archydrotools(fill, Fdr, "")     # Process: Flow Accumulation     arcpy.FlowAccumulation_archydrotools(Fdr, Fac)      any idea ?


title: How ensure data put into PostGIS gets proper CRS attached?

text: I recently had an intractable problem that we finally traced back to the fact that some points in a PostGIS table had a CRS associated with them and others did not. This resulted in some points being displayed in QGIS while other were not, even when they were they were in different layers and all the coordinates were actually using the same CRS. I ended up solving this problem by using postgis functions ST_SetSRID and ST_TRANSFORM to set all the point to the same CSR but I am still not sure if these functions correctly set the CRS and whether or not I will not have the same problem after I next add data into the table. My question now is how do I make sure that all data I put into PostGIS gets the proper CRS attached to the_geom? Doing this has two advantages as I see it:   1. you don't get prompted to set the CRS each time you build a layer from the table (I was puzzled why this happened with my 'old' data).   2. you don't get wierd problems with QGIS when you try and build layers from the data in the tables ;) The points that had no CRS associated with them were all entered into PostGIS using the SPIT plugin from shape layers that I had saved to disk. Recently some kind soul on stackGIS pointed out that I could short circuit the process by using the database manager to save a layer to PostGIS. This appears to save layer CRS with the data into PostGIS. From now on I intend to use the database manager to save data into PostGIS an the assumption that this will protect me from headaches like this in future. Is this a reasonable approach? Apart from getting answers to the question I was also motivated to write this up as an easier to read summary of the question and answers referenced above. In getting to the answers there we found several red herrings that have confused the basic issue and I have left those out here.


title: use raster properties output in attribut table via inline variable in model builder

text: I have a small model that should find the maximum value of a raster and use this number to fill a field in a shapefile. For some reason the variable of the rasters max value isn't accepted in the calculation of the new field in the shapefile. I just get a general 999999 error. I've attached an image for better understanding of the problem. The field I create in the model is of type DOUBLE because I read that the tool 'get raster properties' uses this datatype as output. I also tried using the modelbuilder only tool 'collect values' to save the maximum value in and feed it then into the calculation. Obviously that didn't work, otherwise I wouldn't ask here :( Can anyone please explain to me why this isn't working and how to solve the problem? Cheers, Thea ![enter image description here](http://i.stack.imgur.com/U02jw.png)


title: Is it possible to rearrange the items in qgis layer context menu?

text: I'm referring to the menu I get when I right-click on a layer in the TOC. I'd like to move Properties to the top. I can't find anything in the customization dialog that seems to apply.


title: How can I add Raster Map as permanent layer òr Level to Microstation .DGN file

text: Hei Am having Raster Map and would like to visualize permanently whenever the DGN Opens , I dont wish to attach them as reference in Raster Manager but rather save it as one of the Level in Microstation.DGN file thank you in advance for all your suggestion or solution.


title: Geoserver Print Module Help

text: I am using Heron Mapping Client 1.0.0 to build a web application which renders the shapefiles from Geoserver 2.4.4 . I am on Apache Tomcat 7.0 with geoserver 2.4.4 deployed as WAR file. I am trying to incorporate PrintDialog function in my html code, and I have installed Geoserver Print Module successfully. I can see the list of parameters when opening following links in the browser :                 http://localhost:8080/geoserver/pdf/info.json      This is my html code for print dailog according to heron's specifications               {type: "printdialog", options: {url: 'http://localhost:8080/geoserver/pdf', windowWidth: 360        , showTitle: true        , mapTitle: 'Test'         , mapTitleYAML: "mapTitle"                , showComment: true         , mapComment: 'Test'          , mapCommentYAML: "mapComment"         , showFooter: true          , mapFooter: 'Test'          , mapFooterYAML: "mapFooter"                 , printAttribution: true                  , mapAttribution: null                     , mapAttributionYAML: "mapAttribution"          , showOutputFormats: true          , showRotation: true          , showLegend: true         , showLegendChecked: true         , mapLimitScales: false         , mapPreviewAutoHeight: true           , mapPreviewHeight: 400     }}      This code is added as a tool and then clicking on the respective button will result in opening up of a dialog where a user can modify the printing data. However, when i click on the "print" button , I get " internal server error :500" as a message in the browser. while inspecting element in the browser, i get the following error message in console :               POST http://localhost:8080/geoserver/pdf/create.json 500 (Internal Server Error)      Can anyone help me with this ?? I dont think this problem is web-client related but in geoserver configuration instead.


title: save a WMS image on disk and apply a classitem (gdal commands?)

text: I need to save an image on disk trough WMS coloring the pixels applying a CLASSITEM (e.g. Mapserver CLASSITEM). Is it possible using gdal, if so, how? Thanks in advance, D


title: How can I Normalise Drunk Driving and Postcodes data?

text: I have been provided a dataset of drink driving incidents aggregated into postcodes. This is where the incident happened, not where the driver lives. Postcode 1 - 12 incidents   Postcode 2 - 163 incidents   etc. This is for an entire state in Australia, so just showing total counts will pretty much be a population map. Therefore I will need to do a bit of normalisation. With postcodes not being equal-area, my obvious choices are either area or population (which I have the information for. Population might not be accurate, as there might be factors like particular rural areas being notorious for drink driving due to no public transport, or other factors. Area seems incorrect, as lets say we have two equal area postcodes, one rural and one urban, is it fair to normalise by area alone? I then considered other options such as calculating the total lengths of roads within each suburb and using that. I don't have traffic data, but in an ideal world, would using some sort of traffic count data be also an option? The essential question I want the map to answer is, "Where are the hotspot areas for drink driving within this state of Australia?" How would you go about this? \--EDIT-- Final result http://au.news.yahoo.com/qld/video/watch/21428474/drug- and-drink-driving-hotspots-revealed/ The news guy didnt quite follow the script/advice I provided, and I really wanted to avoid 3D, but you know media. Anyhow, not perfect, but appreciate your help!


title: Select Duplicate Points PostGIS

text: I have a table made of merged point datasets in PostGIS. I want to select all of the features where two or more points have identical geom values. At present I just want to select so I can review the associated field values.


title: ERROR Add Basemap using Arcmap10

text: ![enter image description here](http://i.stack.imgur.com/zIE2N.png) When i have to use Bing Map Aerial, it will come NVIDIA error, then closed on main page. how to solve this problem? thanks to advance


title: QGIS – Batch maps (atlas plugin style) in same location

text: I've got more than 200 shapefiles of points representing different information for the same location. I would like to make as many maps as shapefiles with identical background and extent. I'd rather use QGIS, but ArcGIS 9.3/10 would be ok too. I found EasyPrint and Atlas, which are great, but both are zooming on the shapefile.. which I don't want to. Atlas is way better as it lets you use your own print composer.. but it only wants one shapefile (coverage) while I've got many. Any idea on what I could use? I don't mind a bit of programming, especially because that would save me a looooot of time! Thanks a lot!


title: How do i get ArcObjects not to ignore a layer definition that has been set?

text: I am creating a simple text search application, where the user can enter an address search and the result appear in a listbox, which they can then select it and zoom to the address. Now the problem I am having is that I have set layer definition, within the layer I am working with to only show Residential addresses ( and this works, as i can only see residential addresses in the attribute table), however when I enter a street (in the app), that contains not only residential address, the records appear in my list box when though I have set the layer definition. Is this behavior normal ? and what is the best way to get round this? Thanks


title: How to change the geometry type from Point to Multipoint within an existing table in PostGIS?

text: Is there a PostGIS function that can change geometry type for an existing table? We need to change from POINT to MULTIPOINT. The table will be empty when we change the geometry type and we cannot just drop/create the table.


title: Improving Map service performance

text: I have a map service that has more than 60 data layers. For faster performance, I have separated the data with daily updates as a dynamic service. I still have a map service that base data (50+ layers). I have implemented scale rendering to make sure that only relevant content is shown at a given scale. I am still not getting good speed and trying to figure out the best way to make this service even faster. I am looking for some MAJOR strategies that would allow to make it faster. These services are intended to be used within a Silverlight application. Appreciate any feedback. Thanks Jay


title: How do I change an existing srs in QGis?

text: I am showing postgis data with srid 28992 (Amersfoort / RD New). But when I plot the geospatial data on a map (with openlayers plugin), the data is shown a hunderd meters off because the datum shift is not taken into account. When I make a custom CRS with the correct datum shift and set my layer to that CRS is it shown correctly. I want to update QGis' standard interpretation of 28992 so it takes into account this datum shift (+towgs84 parameter) correctly and when I open new data it is shown correctly immediately. I found that QGis has a sqlite db called srs.db in /usr/share/qgis/resources/. But even when I update the record for 28992 it still doesn't show the data correct on the map. When I inspect the layer in QGis I see it still has no +towgs84 parameter. Does this mean that QGis does not actually use this sqlite db? Because the data I show is from a postgresql database I also edited spatial_ref_sys.sql to make sure that it also knows about the datum shift. I also edited the epsg en esri files in /usr/share/proj/ just to be sure, but no effect. At this moment I cannot think of another place where QGis gets its projection information from.


title: Shapefile with overlapping polygons: calculate average values

text: I have a very big polygon shapefile with hundreds of features, often overlapping each other. Each of these features has a value stored in the attribute table. I simply need to calculate the average values in the areas where they overlap. I can imagine that this task requires several intricate steps: I was wondering if there is a straightforward methodology. I’m open to every kind of suggestion, I can use ArcMap, QGis, arcpy scripts, PostGis, GDAL… I just need ideas. Thanks!


title: CQL filter causing slowness in Java \ Geoserver app

text: I have a Java application that includes an OpenLayers map with WMS tiles pulled from Geoserver. Geoserver and the Java app are looking at the same MySQL database. In the Java application I have a complex filtering system to select data and I would like only the selected data to be shown on the map. I have tried to do this by passing a CQL filter of the form "`ID in (123,3000,6001,5464, etc...)`" unfortunately Geoserver 2.3 takes around 12 seconds to render each tile as the list of ID's is so long and I can expect this to get worse as more data is added. I cannot pass through the SQL query that my Java app is using as it is a complex join that pulls in data from tables not available to geoserver. Also there are security issues as I do not want to expose anything more than the geometry and the feature ID to the public internet via geoserver. Is there a better way of passing the ID's across that is faster?


title: Error 'Invalid geography. Antipodal edge...' with Postgis

text: I have a table in postgis like this               CREATE TABLE locations (id int, name VARCHAR(255), pos geography(POINT, 4326));      with a couple of points. E.G.               INSERT INTO locations (name, pos) VALUES ('freiburg', ST_GeographyFromText('POINT(7.86 48.0)'));      But importing the table as a layer in Qgis results in multiple errors               Query: DECLARE qgisf0_0 BINARY CURSOR FOR SELECT      st_asbinary(st_force_2d("pos"::geometry),'NDR'),"id" FROM      "public"."locations" WHERE "pos" && st_makeenvelope(-180,-90,180,90,4326)      returned 7 [ERROR:  Invalid geography. Antipodal (180 degrees long) edge      (180 -90,180 90) detected, add a point between to make two edges that span      less than 180 degrees.]      and nothing gets drawn. My QGis is 1.9.0, postgis 2.0.3


title: Convert data to another Coordinate Reference System

text: Assuming I've a point like this : > POINT (751318 9760879) This point is in EPSG:27700 coordinate reference system.   How do I convert the point to display in a map with EPSG3857 coordinate reference system ?   ## Update I don't use any GIS application. I am developing a web application (using leafletjs ). So I don't want to convert data with a tool,I've to convert in code.


title: gdal python set projection of a raster - does not work

text: Gdal 1.9.2, ECW driver: ECW/ERDAS Compressed Wavelets (SDK 3.x) I am having problems with the script that supposed to set projection of the ECW to EPSG:27700, i.e. the British_National_Grid. Please could you have a look at the script and suggest what could be wrong? I created two versions of teh script and neither of them work :( I have checked and it does not seem to be the ECW file format/ECW sdk issue as the script does not updte the projection for Geotiffs either. Any help much appreciated! many thanks script 1:               import sys, os, os.path, string, array, math, shutil, glob, re, locale, math, traceback, time, datetime, csv     from osgeo import gdal, gdalconst, osr     from osgeo.gdalconst import *     import subprocess          driver = gdal.GetDriverByName("ECW")     driver = gdal.GetDriverByName("GTiff")          # set OSGB96     sr = osr.SpatialReference()     sr.SetProjection ("EPSG:27700")     sr_wkt = sr.ExportToWkt()     print sr_wkt          inDir = str(sys.argv[1])     for file in os.listdir(inDir):         ds = gdal.Open(inDir + "\\" + file, gdal.GA_Update)         if ds:             print('Updating projection for ' + file)             res = ds.SetProjection(sr_wkt)             if res != 0:                 print('Setting projection failed ' + str(res))             ds = None # save, close         else:             print('Could not open with GDAL: ' + file)      and script 2:               import sys, os, os.path, string, array, math, shutil, glob, re, locale, math, traceback, time, datetime, csv     from osgeo import gdal, gdalconst, osr     from osgeo.gdalconst import *     import subprocess          driver = gdal.GetDriverByName("ECW")     driver = gdal.GetDriverByName("GTiff")          ds_reference = gdal.Open('D:\\reference_geotif\\reference.tif', gdal.GA_ReadOnly)     prj_wkt = ds_reference.GetProjectionRef()     print prj_wkt          inDir = str(sys.argv[1])     for file in os.listdir(inDir):         ds = gdal.Open(inDir + "\\" + file, gdal.GA_Update)         if ds:             print('Updating projection for ' + file)             res = ds.SetProjection(prj_wkt)             if res != 0:                 print('Setting projection failed ' + str(res))             ds = None # save, close         else:             print('Could not open with GDAL: ' + file)


title: How to add an attribute field to an existing Shapefile via Python

text: I have a Python script that adds an attribute field to a Shapefile if doesn't exist. This is easy to do with ArcGIS (graphically or via Python), but I'm looking for something that doesn't depend on ArcGIS. I tried this unsuccessfully with OGR, since my Shapefile contains features. I've looked at pyshp, but similarly there is no way to modify the schema after it has been created. I haven't had a try with shapefile (for Python), but I don't see this feature advertised. I also can't see how this can be done by tinkering with the DBF file via dbfpy. Does anyone have any ideas? * * * **Solved** I found a solution using OGR and thanks to help from a previous question. Here is a complete example:               from osgeo import ogr          # Open a Shapefile, and get field names     source = ogr.Open('my.shp', 1)     layer = source.GetLayer()     layer_defn = layer.GetLayerDefn()     field_names = [layer_defn.GetFieldDefn(i).GetName() for i in range(layer_defn.GetFieldCount())]     print len(field_names), 'MYFLD' in field_names          # Add a new field     new_field = ogr.FieldDefn('MYFLD', ogr.OFTInteger)     layer.CreateField(new_field)          # Close the Shapefile     source = None      My problem was that I used `layer_defn.AddFieldDefn(new_field)` rather than `layer.CreateField(new_field)`. Many thanks to the help, and sorry for not throughly checking for the similar other question.


title: What ratio scales do Google Maps zoom levels correspond to?

text: Can anyone provide me with a link (or some details) on the actual ratio to "zoom level" figures for Google Maps? e.g. **Google Maps Level 13** = 1:20000


title: How to define a line's width graduated by its attributes?

text: I have a line shapefile with lots of lines that have different attribute values. I want to present these lines in different widths according to their attribute values without classes but a graduated symbolization based on each attribute value. Is there any way of doing this without defining different classes by myself? It is intended to create a radial flow map with different line widths as a combination of the following examples: Creating a Radial Flow Map and: Generating Distributive Flow Maps I would prefer to do it with ArcMap, but I am also open to trying it with QGIS.


title: Is there a GroupRenderer in ArcGIS for javaScript?

text: Is there a GroupRenderer or similar functionality in ArcGIS for JavaScript? What I would like to do is to show PictureMarkers with a label showing an attribute on the Graphics object.


title: Can concurrent processes be run in a single model?

text: In model builder (for Arc GIS 10) is it possible to run two or more geoprocessing operations concurrently? For example I have to run repair geometry on 2 files and then merge them. Currently it repairs topo of file A and then file B and then merges them. Each file takes about 30h to repair. Can I get it to concurrently repair A & B, wait for both to be complete and then merge them. At the moment I have to open two instances of arc catalog and run the repair tool on a file in each and once complete do a merge which means that I have to come in to office to check when both are complete...on a quad core, each repair process only uses 1 core at a time. The above example is simple but if multi-threading is possible where multiple operations are run concurrently (if they don't depend on each other) and wait for the depended files to be complete before the next process is run. Thanks,


title: Is there a way to get layers from a .msd using arcpy?

text: I need the name of each layer in an .msd, but I could not find support using arcpy (ArcGIS 10.0). Is there a method I missed? In the current workflow, since I am publishing the document, I use the REST service:               descriptionUrl = url + "/ArcGIS/rest/services/" + mapServer + "/MapServer?f=json"       restResponseFp = urllib.urlopen(descriptionUrl)     mapDescription =json.load(restResponseFp)     layerNames = [str(layer["name"]) for layer in mapDescription["layers"]]


title: Named ITable removal

text: I am dynamically creating a join between a GeoFeatureLayer and an external ITable derived from a SQL query and opened with OpenQueryClass(). It works. Now if I could only delete the Join!? But even worse than that: when I re- enter the code that creates a new join (after I run some code that seems to remove the existing joined fields from the GeoFeatureLayer) it fails to create the ITable the 2nd time because "The FeatureClass already exists." So it appears that not only must I reset the GeoFeatureLayer I must also remove the named ITable. Is this correct? If so, how?


title: Why is my Definition Query not returning any records?

text: I am currently working on a mapping project and only need 8 tenements out of many to appear. I have gone into the query builder to complete this and unfortunately it will not work. Into the builder a put Tenement = 'E59/1002' AND Tenement = 'E59/1421'. It verifies but no records are found? I am not sure why this is happening as this method has always worked for me before when i need to select multiple tenements? Any help? Cheers, Kat


title: How can i increase network analyst performance in mobile?

text: I'd like to use network analyst in mobile. I use arcgis sdk for mobile.the problem is that i use url service for routing .I have a network analyst service in ArcGIS server and use it in mobile.routing in mobile is time consuming.I think because of internet speed i have the problem . Is there any way to increase performance of routing? Thanks for advance


title: Remove item from legend using ArcPy?

text: I can't remove _item_ from the legend (lyr called "atikot") using arcpy with this code:               import arcpy     from arcpy import env          removeLyr = arcpy.mapping.ListLayers(mxd, "atikot")[0]     mxd = arcpy.mapping.MapDocument(r"C:\Project\Project.mxd")     legend = arcpy.mapping.ListLayoutElements(mxd,"LEGEND_ELEMENT")[0]     for lyr in legend.listLegendItemLayers():         if lyr.name == "atikot":             legend.removeItem(lyr)     del mxd      The code does not work. Does anyone have an idea as to what's wrong with the code?


title: Is MapProxy and Google's Static Maps API a possible alternative to the QGIS Openlayers plugin?

text: I've been looking for an alternative to the openlayers plugin as a means of adding imagery basemaps to QGIS. The plugin is great but there's a bug when exporting basemaps via QGIS map composer that causes vector layers to become offset from the underlying basemap, or vice versa. I saw mapproxy mentioned as a possible alternative here, leading me to create this config file which nearly works I _think_ , except for one thing: google's static map api wants tile (yx) location in a lat,long format, whereas mapproxy sends requests in a simple tile numbering system: 0,0 0,1 1,1. I wanted to avoid accessing imagery tiles directly as this is against google's tos. Is using mapproxy and google static maps / bing static maps feasible? Is there another way to get aerial imagery basemaps into QGIS that will work with map composer? Being new to web mapping, I suspect I'm barking up the wrong tree so hoped to get some help here =]


title: No get more option in QGIS

text: I have recently installed QGIS and trying to work with it but I don't see this option "Get More" at all. As suggested by some previous source, I also checked whether I am behind a proxy but that isn't the case. What could be the cause? thanks in advance.


title: Error issues when using Cut Fill operation in ArcGis 10.2

text: I wonder if someone is able to assist with an issue I have about a series of errors I get when using the _cut/fill_ tool in spatial analyst, ArcGis 10.2? Someone posted a question sometime last year entitled: What is causing - Arcmap errors (010302/010067) while using 3D analyst/CutFill? but as you can see it was closed for being too localized. This, however, is my exact problem, and I have seen that the page has been viewed many times, so I don't think I am the only one! The author of the question said that a miraculous solution to the problem manifested itself the following day; I have not been so lucky. Does someone know of a solution to this?


title: Read feature table with FME / Python

text: I try to convert a feature table Geodatabase to MapInfo with FME / Python But each test, my application crash. There is no geometry in my feature table in file geodatabase. Why my application crashes with attribute table ?                   writerTab=FMEWriter("MITAB")     writerTab.open(RepMapInfo)          schemaTab=FMEFeature()     while reader.readSchema(schemaTab):         #print writerTab.getSchemaFeatures()         log.logFeature(schemaTab)         writerTab.addSchema(schemaTab)          while reader.read(schemaTab):         log.logFeature(schemaTab)         writerTab.write(schemaTab)      The code is very simple. Help me :)


title: Online Routing Service (non-Google)

text: im currently blueprinting something i wanna develop. Thinking about the stuff that would be needed, i figured out that i would need to have a pgRouting/postGis Database somewhere. The program i wanna make, should calculate data and return the shortest path. Is there any online service who offer this kind of service (an accessible pgRouting/postGIS Database)?? even for a fee ?? Since i couldn't find one im asking here. Thanks a lot. * * * Update from comments: I want to develop an app that USES a routing service, to display the (elsewhere) routed data on the phone. Like oh i am here, thisService->please tell me how to get from x to y.


title: Measure each angle of a selected polygon

text: This post is also a continuation of a previous one. I have the code (in VBA and ARcObjects) shown below to measure each angle of a selected polygon in degrees. I have the following problems: 1\. Some angles are negative and I do not wish this since I measure the interior angles of a polygon. 2\. Some angles are measured wrong i.e. I get the reflex angle instead of the normal angle. e.g. 270 instead of 90 degrees. Could somebody see the code and help me for the neccessary modifications please? There is a procedure that calls a function shown also below. Thanks.                   'This procedure gives the angles between the vertices of a polygon      ' ----------------------------------------               Dim pMxDoc As IMxDocument     Dim pGeoFtrLyr As IFeatureLayer     Dim pFtrCls As IFeatureClass     Dim pFtrCsr As IFeatureCursor     Dim pFtr As IFeature     Dim pGeomColl As IGeometryCollection     Dim pSegColl As ISegmentCollection     Dim pCurve As ICurve     Dim pLine1 As ILine     Dim pLine2 As ILine     Dim dVtxAng As Double     Dim i As Long, j As Long     Dim PI As Double          '    ' Calc PI     PI = Atn(1) * 4          ' Get a cursor on the features in the first layer     Set pMxDoc = ThisDocument     Set pGeoFtrLyr = pMxDoc.FocusMap.Layer(12)     Set pFtrCsr = pGeoFtrLyr.Search(Nothing, False)               'Select a polygon          Dim pEnumFeature As IEnumFeature     Set pEnumFeature = pMxDoc.FocusMap.FeatureSelection     pEnumFeature.Reset     Set pFtr = pEnumFeature.Next     If pFtr Is Nothing Then        MsgBox "Please select a valid feature", vbInformation      Exit Sub     End If     If Not TypeOf pFtr.ShapeCopy Is IPointCollection Then     MsgBox "Not a valid feature", vbInformation     Exit Sub     End If                   ' Get all disjoint paths         Set pGeomColl = pFtr.ShapeCopy              ' Loop thru the paths         For i = 0 To pGeomColl.GeometryCount - 1             ' Get all segments that make up this path             Set pSegColl = pGeomColl.Geometry(i)             Set pCurve = pSegColl                  ' Loop thru the segments             For j = 0 To pSegColl.SegmentCount - 2                      ' Get the next two lines                 Set pLine1 = pSegColl.Segment(j)                 Set pLine2 = pSegColl.Segment(j + 1)                      ' Calculate the left side angle (in degrees)                 dVtxAng = CalcAngleBetweenLines(pLine1, pLine2, True)                  MsgBox Abs(dVtxAng * 57.2957795)                  Next j                  ' Check for a closed polylin                  If pCurve.IsClosed Then                 Debug.Print pFtr.OID; " is closed!"                 Set pLine1 = pSegColl.Segment(j)                 Set pLine2 = pSegColl.Segment(0)                 dVtxAng = CalcAngleBetweenLines(pLine1, pLine2, True)                      'get each angle of the polygon in degrees                 MsgBox dVtxAng * 57.2957795                  End If              Next i                       End Sub               Private Function CalcAngleBetweenLines(pLine1 As ILine, pLine2 As ILine, AsDegrees as Boolean) As Double               Dim dAng1 As Double, dAng2 As Double, dVtxAng As Double, PI As Double          dVtxAng = pLine1.Angle - pLine2.Angle               ' Get the reversed angle of the first line     pLine1.ReverseOrientation     dAng1 = pLine1.Angle     If dAng1 < 0 Then dAng1 = 2 * PI + dAng1     pLine1.ReverseOrientation          ' Get the angle fo the second line     dAng2 = pLine2.Angle     If dAng2 < 0 Then dAng2 = 2 * PI + dAng2          ' subtract 2 from 1 - that's the smallest angle between the two lines     dVtxAng = dAng1 - dAng2          ' Ensure we have the angle to the left of the lines          If dVtxAng >= 2 * PI Then     dVtxAng = dVtxAng - 2 * PI     ElseIf dVtxAng < 0 Then     dVtxAng = dVtxAng + 2 * PI     End If               CalcAngleBetweenLines = dVtxAng     End  Function


title: Is it possible to display the length of each segment with the Openlayers measure control?

text: I am trying to make a OpenLayers measure tool that displays the length of each segment on top of each segment and the total length of the measurement as shown below: ![Measure control with length of each segment](http://i.stack.imgur.com/p1mZJ.jpg) I have used this example http://openlayers.org/dev/examples/measure.html. I can see in the example that the geometry created is one big LineString. Therefore my general thought was to get the measure control to create a new line segment (LineString) for each of the user click except the first, then add an attribute with the segment length and display it with the label parameter in the style. But I don't know if it will conflict with the current design of the measure control or if this is the correct approach to my problem? All comments or guidance would be appreciated.


title: Right to left language layout for automatic KML popup tables

text: Using the `Layer to KML` tool in ArcMAP generates a pop-up table that has left to right orientation, i.e the title is on the left and the value is on the right, like the following exampl ![Left to right orientation table](http://i.stack.imgur.com/yxEFj.png) I need to create a right to left kml popup table, like this: ![Right to left orientation](http://i.stack.imgur.com/vriMQ.png) I know I can edit them manually to appear that way, but I have hundreds of elements. The ArcGIS tool doesn't seem to be able to do it, so is there any other KML tool that can (not limited to ArcGIS, any other platform will do)? My original data is a point feature layer, with several fields.


title: Converting .multiple shapefiles into Image Pyramid?

text: I have multiple shapefiles, and I want to import them as image tiles. How do I convert a shapefile into image format and then convert that single image file into an image pyramid such as that of OpenStreetMap (OSM)? I have worked with GeoServer for backend and OpenLayers to render the map but now I am asking about converting my own custom shapefiles into OSM tile pyramid.


title: How to read a feature table empty with FME / Python?

text: I would like to convert a feature table Geodatabase (empty) to MapInfo, but my feature table is empty. I would like to write the schema of feature table in MapInfo file. Is it possible to **force** the writer MapInfo ? Thanks


title: Seeking commercial web service for LIDAR/DSM data in USA

text: I apologize for not knowing all the proper terminology to ask this question. I am seeking a commercial service for providing pay-per-use high resolution elevation data in North America. It must come with defined ToS (terms of service -- a guarantee of availability). A paid service is best. The ideal service would return a chunk of elevation data, perhaps in ARC/INFO ASCII GRID format, on demand. Ideally this grid would be seamlessly centered on a given coordinate. For this application I need a regular grid of elevation points perhaps 2000 or 3000 grid elements in the x and y direction (if the data returned were richer LIDAR data, I would need to interpolate into a grid using las2dem or some such). I am aware of scattered sources of free LIDAR elevation data:   1. http://lidar.cr.usgs.gov/   2. http://www.opentopography.org/   3. http://calm.geo.berkeley.edu/ncalm/ddc.html I'm seeking something that is (a) offering wider coverage with less hassle (b) a web service returning data in close to real time (c) supplying seamless data centered on a point on earth, rather than tiles (d) commercially usable and supported with terms of service (e) command line scriptable (e.g. using wget or other scriptable tools outside of a GIS). For example:               wget http://ppv.example.org?lat=37.1234&lon=-122.1234     &rows=1024&columns=1024&format=ascii_esri_grid&data_type=dsm:best_available      Is this a crazy desire, or is there something that fits the bill? _Possible question keywords: WFS, web service,hosted LIDAR data server_


title: Multiple ring buffer with attributes

text: Building on ArcGIS: inside buffer with attributes?, I'd like to be able to create multiple ring buffers, while retaining (and possibly dissolving on) attribute values. How can this be done?


title: How to determine flow amount along a stream using ArcGIS 9.2?

text: I have DEM with resolution 20m and also stream network map in shape format, isohyetal map in shape format. I have also historic numerical data for rainfall, runoff and temperature from 1970-2000 for just one gauging station in the area 10x14 km = 140 sqm. can anybody give me idea how can I calculate the stream flow amount(Q) for entire area and especially along the main stream. Thanks in advance. Subija.


title: Visual Basic editor command icon is always disabled

text: I recently installed ArcGIS 9.3 in a PC with Windows 7.0. I opened an existing mxd file within which I developed some toolbars for a project. Each icon of a toolbar has behind it a VBA and ArcObjects code. Unfortunately, I noticed that icons are all missing and the symbol of each icon has been replaced by a red circle with a slahed line and the commnad of Visual basic editor under the Tools menu is always disabled-inactive, thus I cannot get in the VBA editor. Can anybody help me please?


title: QGIS 2.2 crashes when saving digitized data in ArcSDE 10.1

text: I am using QGIS 2.2 with ArcSDE 10.1 for Oracle 11gR2 I have created a polygon feature class with SDO_geometry, when trying to save the digitized data QGIS crashes. Please help, I want to edit the ArcSDE data in QGIS.


title: OpenLayers Select a feature from attribute or id

text: Using Openlayers.....Is it possible to select a feature on a map by an attribute? eg: Select a feature from a dropdown menu populated with the feature "id" or a table field attribute values.


title: 3D KNN index in PostGIS

text: I am looking for a way to do fast KNN queries on a dataset of 10+ million 3D points (and growing). I was hoping PostGIS would do the trick with its new <-> and <#> operators, but they don't seem to work in three dimensions. Does anyone know any workarounds that would extend the KNN GiST index to 3D? Right now I'm using a "normal" spatial index with the &&& operator, but that creates all the issues that led to the addition of the KNN index in the first place.


title: Using OpenLayers and OpenStreetMaps to display statistical information with colors : possible?

text: I'm planning to build a statistical tool using OpenLayers and OpenStreetMaps. This tool will have to display statistical data for french departments with colors. What I would like to do is the same as this : http://sd1878.sivit.org/geoclip_stats_o3/index.php?profil=FR#v=map3;i=ocs.p_c1_2006;l=fr;z=178717,7038714,864706,751147 Does OpenLayers/OpenStreetMaps allow this ? Or should I look for other alternatives ? Thanks


title: Anyone recognize this (reference-system related?) pattern of distortion?

text: I'm new to GIS and was making some progress until getting rather stumped on this one. I've imported what I believe to be lat/lon data based on US FIPS county codes in weather alerts (NOAA SPC). When I use the geodjango "simplify" method on my geometries, they render on a map according to my expectations. But when I attempt to perform relationship queries on the raw geometries, I get distance and overlap values that do not correspond to the real world geography. For example, the two upper blobs in this top image should be adjacent and much further south, and the lower blob should fit Florida better (the distances and overlaps in this image match my ST_Distance and ST_Overlaps values): ![Distorted geometry](http://cl.ly/image/323I1V1U2D37/Screen%20Shot%202013-06-19%20at%204.55.20%20PM.png) Here is a rendering of the "simplify" geometry, which shows how these spatial entities are truly related in the real world (pink regions): ![Correct geometry](http://cl.ly/image/1c2D1W2M3c2e/Screen%20Shot%202013-06-19%20at%204.54.56%20PM.png) Does this distortion look familiar to anyone? Thanks!!


title: Improve display of fine resolution raster

text: ![Map of yields in Iowa](http://i.stack.imgur.com/v45gF.png) In the above map, I am showing yields on some parts in Iowa at a pixel resolution of 100m. However, in this case due to the relatively small number of pixels (and high resolution), the map is not very clear. Is there any way I can smoothen out the raster so that it displays nicely using ArcGIS? Thanks to @ Michael Miles-Stimson for suggestions, the yields that I am displaying are the 'VALUE' of the raster. Therefore it is not possible to do something like a focal mean (which destroys the raster). I would rather avoid changing the attribute table to use the yields as 'VALUE' because there are several other variables I want to plot as well, and this would loose information and be very time-consuming for my raster.


title: ArcMap error: displays black screen and locks up

text: I have been getting an unusual error in ArcMap lately. The entire ArcMap window turns black (with the exception of the Table of Contents, Attribute Tables and any other tools). It has only happened in Layout View so far, and the black area covers the page, as well as the area behind it. This occurs when the map is refreshing and I've got an ESRI basemap in my map. ArcMap is unresponsive when this happens. When I try to "End Task" through the Task Manager, I get a message that the program is waiting for a response from me (usually to save). When I click "Save" in ArcMap, the program shuts down, but when I reopen the map, it has not been saved. This has been happening quite often, with a variety of map documents and basemaps. I've searched online, but haven't found anyone else having this problem or any possible solutions. Any ideas how to avoid this problem?


title: how to add data from a layer to postgis table with different field names

text: I am importing waypoints from GPS (in gpx format) using Layer>Add Vector layer and then I want to clean up the data and add some more attributes and finally add them to a master table in PostGIS. The import works fine but I get a table with a lot of fields I don't need and the fields I do need in some cases have the wrong names nor can I modify the table without doing a save as and reloading it. So I end up   1. importing the layer,    2. saving it as a in a format that I can modify    3. importing it again   4. clean up the data and add new attributes (using some of the fields i don't need) Now I want to update the master table. The only way I can see of doing this is to save the layer as PostGIS (or do this in step 2) and then use SQL to do the update so I can select which fields I want inserted in which columns. Is there a more stream lined way of doing this? PS. I am guessing the answer is write a script to do it. To this end I have just purchased The PostGIS cookbook and the Pyqgis book but it will be a few weeks while they travel around the world.


title: How to highlight on mouseover without wiping out select control

text:
I am trying to highlight features of a vector layer as the mouse moves over them. I can make this work. I'm also trying to make the features selectable and have a popup box when they're clicked. this works too. what I cant make work is both of them at the same time. The hoverselect seems to overide the popup select. I've based my code on http://dev.openlayers.org/releases/OpenLayers-2.11/examples/highlight- feature.html but with a few differences. Does anyone know why the hover select and select can't coexist?               select = new OpenLayers.Control.SelectFeature(vector_Layer, {clickout: true});              vector_Layer.events.on({                 "featureselected": onFeatureSelect,                 "featureunselected": onFeatureUnselect             });             map.addControl(select);             select.activate();                  select.handlers['feature'].stopDown = false;              select.handlers['feature'].stopUp = false;          //REMOVED AS OVERRODE SELECT selectHover = new OpenLayers.Control.SelectFeature(vector_Layer, {hover: true});          function onPopupClose(event) {         select.unselectAll();     }          function onFeatureSelect(event) {             var feature = event.feature;             var content = "
"+feature.attributes.SHEET.value + "
" + feature.attributes.description;             if (content.search("


title: Building SQL statement in javascript and casting variables as strings

text: I'm working with Postgres 8.3/PostGIS which is storing some polygons. I'm writing a javascript function which will take in a parameter of a watershed id from another function. The javascript will build up a sql statement and then send that along to a jsp script which will send that to the database and return some geojson which will then be output into Openlayers. The problem is that the geojson is not being returned because I think I have a problem with my select statement and some type (integer/text) problem. The table holding the polygon geometry contains a column set as varchar(50) but it contains some id's which have only numbers, and some id's which contain numbers and string characters. The function attempts to return a geojson from the database of a polygon of a given id. I've left off the last end of the function because I think my problem is somewhere in the building of the sql statement.               function displaySelectedShed(shedId) {         var json_url = "04-sql-to-json.jsp?sql=";         var sql_geom = "st_transform(ST_Force_2D(the_geom),900913)";         var sql = "select st_asgeojson(" + sql_geom + "), id ";             sql += "from medford.subshed2 ";             sql += "where subshed2.id = " + shedId +"";             json_url += encodeURIComponent(sql);      }      The select statement above looks like this:               select st_asgeojson(st_transform(ST_Force_2D(the_geom),900913)), id from medford.subshed2 where subshed2.id = 01616500      and I can test this by putting this:               http://localhost:8080/spatialdbtips/04-sql-to-json.jsp?sql=select%20st_asgeojson%28st_transform%28ST_Force_2D%28the_geom%29%2C900913%29%29%2C%20id%20from%20medford.subshed2%20where%20subshed2.id%20%3D%2001616500      into a browser window. But the error I get is this: Error: operator does not exist: character varying = integer I know this has probably something to do with the fact that my shed id is a number when inserted as a variable into the sql statement, but I've tried all sorts of cast combinations to get it to be a string, since the values in the database are stored as such, but I can't get them to work. I've tried both the CAST as well as the :: such as sql += "where subshed2.id = " + shedId +"::text"; but I can't get it to work. Does anyone have any suggestions to get my watershed id numbers (the shedId variable) to act like string in the sql where clause, and how to build it up in a way I've shown in javascript?


title: Why does ArcGIS 10 Calibrate Routes give Error 000443?

text: I am trying to calibrate routes using points along them. Those points are represent counts of hikers per hour. I have done a spatial join between my route layer and my points to have the RouteID in common. Some points do not match any route. Then, I tried to use the **Calibrate Routes** tool but the attribute table of the new shape is exactly the same as the one of the Route shape. And if I try again I get an error 000443. Would that be due to the fact that some of my points do not match any route? Should I calculate an average value of the counts for each point before I use the calibration? Thank your for your help!!


title: openlayers:transform from ITM to 900913 is not accurate

text: i did a custom projection ITM (Israel TM) and i then am doing transform to all the feature before im adding them to the map, from ITM to 900913 proj because i want to show the feature on the map. this is the string for the ITM proj:(i am using it with proj4js)               +proj=tmerc +lat_0=31.73439361111111 +lon_0=35.20451694444445 +k=1.0000067 +x_0=219529.584 +y_0=626907.39 +ellps=GRS80 +towgs84=-48,55,52,0,0,0,0 +units=m +no_defs       but, when i add the feature like this:                           features.geometry.transform(ITM, proj900913);                 vectors.addFeatures(features);      the feature is place 4 meter to the north and 4 meter to the east, and not in the place its suppose to be. its also happening when i am load KML layer. how can i place the point in the right place? thx


title: Are CartoDB Map views by month or by year?

text: We have 10000 map views with the free version of CartoDB. Are they reset every month or is it the maximum views for the all utilisation ?


title: How To calculate raster land cover inside overlapping buffer zones?

text: I am working with a raster land cover file that consists of several land cover types (e.g. forest, buildings, grass ect.) in total there is 40 different land cover types. I have 1200 bufferzonen (vector file) each on 500 meters. I would like, for each bufferzone (they have all a unik_ID) to calculate how many square meters that consist land cover type 1 (forst), 2 (buildings), 3 (grass) ect. Some buffer zones are overlapping.... In the end I would like one crosstable that contains all the 1200 bufferzones and the square meters of each landcover type inside all the buffers. I´m using ModelBuilder.


title: Add New Project -> Class Library (Engine) targets the .NET 2.0 framework

text: Any one have comments on why the ArgGis Extensions in the 10 framework still target the .NET 2.0 framework? For clarification do the following steps in Visual Studio: Add New Project -> Class Library (Engine) -> Click properties -> see that it targets the .NET 2.0 framework.