title: hardcoding port number in Arcgis License Manager

text: In the Arc10 Q&A for the 2011 user conference it says a performance gain may be realized by hardcoding the license manager to port 27000, and then instructing the firewall and Windows Defender to ignore that port. That same doc references the online help for the specifics of how to implement that on the license manager server, namely edit `services.txt` to look like the below (the trailing numbers on first two lines used to blank):               SERVER this_host ANY 27000     VENDOR ARCGIS PORT=1234     FEATURE ACT ARCGIS 1 permanent 1 7ED49106D630      This is all well and good to this point, but I'm confused about what is happening with the server port and the vendor ports. Which way does the communication flow go? If the ArcGIS Vendor port is hardcoded, then does every client need to be updated to explicitly match that port? (e.g. do I have to visit every client machine and set their ArcGIS Desktop Adminstrator to use port 1234? or is it 27000?) Do I set the firewall to ignore one port or both ports? in or out? And is #1234 actually a safe number, or range of numbers?


title: Extract many or clustered cells attributes around points

text: My question is similar to Extract raster values at points, but perhaps more complex. How to extract the attributes of clustered cells around a point/points? For example, I have many vegetation plots in 100m x 100m and have to examine the spectra at the points over SPOT Data (satellite image) which has 10 m resolution. I use ArcGIS 10 but QGIS is quite easy to follow. Thanks in advance for your help. Regards, Wim


title: Linear referencing: Assigning RouteIds to polyline feature class

text: I have a Linear Referencing Route System in ArcGIS Geodatabase format (polylines with M values) and a separate polyline feature class of short segments (called TMC) that do not have RouteIds currently. The Linear Referencing Routes are directional. So, there will be two long lines on top each other in most locations in the Route layer where there are two-way routes, with RouteIds like 1001 and 1002 and these are digitized in opposite directions. How can I assign the appropriate RouteIds to the TMC segment layer based on the digitized direction of the TMC segments? I only need to assign RouteIds when the TMC is approximately parallel to an existing Route and within 100 ft of the Route? We need to ignore the TMCs which do not meet these criteria. So, if a TMC segment intersects a Route, but is perpendicular to it, it should not be assigned a RouteId. Thanks.


title: Is it possible to add photo-realistic clouds in ArcGlobe animations?

text: This question asked about the best software to produce professional fly- through animations, and included a link to a nice video of the Grand Canyon. That video includes photo-realistic clouds in the sky. Is it possible to add clouds like this in ArcGlobe?


title: Where did my Create Features window go?

text: I closed the window and now I can't find it again. I've tried the Editor toolbar dropdown and Create Features is not even on there (at all, not even greyed out). Suggestions?


title: How to fix Graphserver warning "osmdb has no edges!"?

text: I'm really new to GIS. I have some shapefiles data that I'm trying to import in graphserver. I'm doing:               gs_osmdb_compile output3.osm output3.osmdb      And the result is: `WARNING: osmdb has no edges!`. I'm using a shapefile that may not be ready for import. I'm converting it to OSM, then I'm trying to import it using `gs_osmdb_compile`. I guess they are related to the v.net function in GRASS, but I'm not sure to understand what should be the end result. Thank you very much!


title: What is the difference when measuring in units rather than meters?

text: I have a function that takes a point and compares it with lines that are stored in PostGIS and then gets the line closest to the point. Initially I did not transform the data because I did not think that it would matter when I want to get the minimum distance. I thought that small distances (up to 100m) would not need transformation but I was wrong! Sometimes it gets a line that is not the closest one even if the compared distances are 1.2m and 2.5m. I discovered that and I re-designed the SQL command to transform the point and line to the local UTM zone and now its working fine. I would appreciate it if someone could explain to me how is this possible at such a small scale. Why did PostGIS think that a line 2.5m away from my point is closer than a line 1.2m away? Perhaps an emphasis on UNITS vs metric unit usage? Thank you for your time and effort! Regards, A


title: How to do kernel density estimation with a physical barrier/ boundary in QGIS?

text: Does anyone know of a programme that will allow kernelling that takes a boundary into account in its calculation rather than simply masking areas that are impossible? So far I have found: GME add on for Arc10 (Hawthorne Beyer) - I get an error everytime I specify a boundary shapefile. I've tried lots of different shapefile types and boundary complexities etc. It works fine when I don't specify a boundary. AdehabitatHR package in R (Calange 2011) - this works well but the boundary you specify has to be very simple - line segments 3x kernel bandwidth in length, and not too tortuous. For my data this is a big oversimplification. So I'm wondering if any other software can do it, GRASS or QGIS for example. Thanks


title: Where I can find world terrain data?

text: I need free terrain data for the whole world which I can import into TileMill. Is there some source where I can download them? Thank you.


title: How to represent an imaginary flat world in QGIS without using a round-world CRS

text: Like one or two other users, I'm interested in using QGIS to map an imaginary flat world. A previous answer suggested simply refraining from defining a CRS for such a world, but I cannot seem to make that work well with QGIS 1.7 on Mac OS X. Pretty much every time I interact with a layer, the software warns me that no CRS is defined, and if I ever inadvertently open the CRS dialog for any reason, QGIS appears to apply the default CRS to my map. That happens quite often when I'm fumbling around looking for a way to find a setting that will establish the appropriate scale for distances on the map. The fictional world is a flat plane of unknown (and possibly infinite) extent. The map is centered on a location that is locally taken to be the center of the world. Although the world is of unknown extent, the mapped area is about ten thousand kilometers north-to-south, and the same east-to-west. I can map the world adequately using tools like Photoshop and Illustrator, but obviously they aren't GIS systems. I'd like to use QGIS for the myriad additional data-handling features it would provide, but I'm pretty well stumped by my inability so far to find or create a CRS --or to get QGIS to accept working without one. I'm close to concluding that I'll have to write my own quasi-GIS system for this purpose that obeys the (admittedly idiosyncratic) rules of the fictional world.


title: Is it possible to script ISIS Mapper?

text: Does anyone have any experience of scripting or batch processing with ISIS Mapper? I am working on a large modelling and mapping project that will require the processing of significant numbers of 1d models through ISIS Mapper. Has anyone had any experience of scripting tools in ISIS mapper to batch process the extraction of model results or the interpolation of 1D results? A google search and hunt through the information provided by Halcrow doesn't provide any suggestions. Thanks Ed


title: How to install TimeManager for QGIS?

text: I'm trying to install TimeManager for QGIS. I can't find it under Plugins/Fetch Python Plugins and I have verified that the official repository is listed in the Plugins menu, so I'm wondering if there's some other way that I'm supposed to install it. Underdark, am I missing something obvious here?


title: Where to find PySAL toolbox for ArcGIS?

text: I installed PySAL but i can't find the toolbox for ArcGIS.


title: Changing Data Source of Layer in QGIS

text: Is there a way to source a shapefile to a layer file in QGIS much like how you do this in the properties of an ArcGIS layer? (i.e. Go to layer properties, source tab and then just depress the 'Set Dat Source' button) I've been all over the layer properties in QGIS and I don't see this anywhere... Thanks, Mike EDIT: I've received comments explaining that you cannot do this with shapefiles, but if I were to change the name of a shapefile and then open a pre-existing project file which contains this shapefile as a layer, I will get the 'Handle Bad Layers' dialogue, which allows me to navigate and re-source the layer to any shapefile I want to. I find it odd that you can resource in that scenario, but you cannot do it manually through the layer properties dialogue box.


title: Trying to display GeoServer layer controls but nothing is coming?

text:
I am a beginner with GeoServer trying to display layer controls but nothing is coming.                               <title>GeoServer Layer</title>         <link rel="stylesheet" href="http://192.168.0.7:8080/geoserver/openlayers/theme/default/style.css" type="text/css" />                               Geosever Layer

location

         Click on the map to get feature info


title: How to use paths with non ASCII characters for subprocess calls with ArcGIS

text: I know that it is possible to create calls using unicode character encoding like \xf8 for non ASCII characters but how do I do it writing a path in a script just like: 'e:\gis\python_ø\notepad.exe' using windows and ArcGIS 10.1? The code below shows three different methods which all works outside ArcGIS but only method 2 works with a script tool within ArcGIS. In a python shell 'e:\gis\python_ø\notepad.exe'.decode('latin-1') returns u'e:\gis\python_\xf8\notepad.exe' so I thought method 3 would work in ArcGIS but it returns the error: > Traceback (most recent call last): File "E:\GIS\pythonø\su.py", line 26, in > File "C:\Python27\ArcGIS10.1\Lib\subprocess.py", line 493, in call return > Popen(*popenargs, **kwargs).wait() File > "C:\Python27\ArcGIS10.1\Lib\subprocess.py", line 679, in **init** errread, > errwrite) File "C:\Python27\ArcGIS10.1\Lib\subprocess.py", line 893, in > _execute_child startupinfo) WindowsError: [Error 2] The system cannot find > the file specified Does anybody know how I can write non ASCII characters in paths for subprocess calls without the need to write them like \xf8?               # -*- coding: latin-1 -*-     import subprocess     ###    Method 1(str): Works with IDLE     fn1 = 'e:\\gis\\python_ø\\notepad.exe'     subprocess.call([fn1])     ###    Method 2(unicode): Works with IDLE and ArcGIS     fn2 = u'e:\\gis\\python_\xf8\\notepad.exe'     subprocess.call([fn2.encode('latin-1')])     ###    Method 3(unicode): Works with IDLE     fn3 = 'e:\\gis\\python_ø\\notepad.exe'.decode('latin-1')     subprocess.call([fn3.encode('latin-1')])


title: Invalid spatial reference Error

text: I am getting the error: > Invalid spatial reference. Spatial reference must match map's spatial > reference. Clear the map layers collection prior to changing the spatial > reference when I try to Zoom in on a feature I have selected in an infogrid. I will try to explain what I am doing and my code will be at the bottom of this question. I have an application that allows the user to map whatever account they are viewing based on the account number. The query is working and is showing up in the infogrid, but it is not being marked on the map with a renderer like it should and it will not allow me to zoom in and will give the message stated above. The odd thing is, this works perfectly as long as I allow the user to click a button to query. The code is almost identical. The only thing is that the results are not displayed as soon as the map is loaded. I am going to put the code that works and the code that doesn't. Any help is appreciated. C# code (Not working): This code is fired when the user clicks a button on a DIFFERENT control.               public void ShowMap(string account_num, string street, string city, string state, string zip)             {                 LocateAccount(account_num);             }                  void LocateAccount(string account_num)             {                      _accountGraphicsLayer.Graphics.Clear();                           QueryTask queryTask = new QueryTask("http://xx.x.xxx.xxx/arcgis/rest/services/Calvert_City_Test/MapServer/2");                 queryTask.ExecuteCompleted += ShowMapLocatorTask_AccountToLocatonCompleted;                 queryTask.Failed += ShowMapQueryTask_Failed;                      var dirty = DateTime.UtcNow;                 Query query = new Query();                 query.OutFields.Add("*");                 query.ReturnGeometry = true;                 query.Where = "CIS_METR LIKE '%" + account_num + "%'";                 query.OutSpatialReference = MyMap.SpatialReference;                 queryTask.ExecuteAsync(query);             }                  private void ShowMapLocatorTask_AccountToLocatonCompleted(object sender, ESRI.ArcGIS.Client.Tasks.QueryEventArgs args)             {                 // Clear previous results                 AccountDataGrid.Visibility = Visibility.Visible;                 GraphicsLayer graphicsLayer = MyMap.Layers["Account"] as GraphicsLayer;                 graphicsLayer.ClearGraphics();                      // Check for new results                  FeatureSet featureSet = args.FeatureSet;                 if (featureSet.Features.Count > 0)                 {                     // Add results to map                     foreach (Graphic resultFeature in featureSet.Features)                     {                         //resultFeature.Symbol = ResultsFillSymbol;                         graphicsLayer.Graphics.Add(resultFeature);                     }                 }                 else                 {                     MessageBox.Show("No features found");                 }             }                  public void ShowMapQueryTask_Failed(object sender, TaskFailedEventArgs e)             {                 MessageBox.Show("Query failed: " + e.Error);             }      C# Code (Working): This code is fired when the user clicks a button that is inside my mapping control.               private void AttributeSearchButton_Click(object sender, RoutedEventArgs e)             {                 RunQuery();             }                  private void RunQuery()             {                 _attributeQueryGraphicsLayer.Graphics.Clear();                      QueryTask queryTask = new QueryTask("http://xx.x.xxx.xxx:6080/arcgis/rest/services/Calvert_City_Test/MapServer/2");                 queryTask.ExecuteCompleted += GeneralQueryTask_ExecuteCompleted;                 queryTask.Failed += GeneralQueryTask_Failed;                      var dirty = DateTime.UtcNow;                 Query query = new Query();                 query.OutFields.Add("*");                 query.ReturnGeometry = true;                 query.Where = "CIS_METR LIKE '%" + QueryParametersTextbox.Text + "%' OR COMMENTS LIKE '%" + QueryParametersTextbox.Text + "%'";                 query.OutSpatialReference = MyMap.SpatialReference;                                queryTask.ExecuteAsync(query);                            }                  // Draw results when query is complete             private void GeneralQueryTask_ExecuteCompleted(object sender, ESRI.ArcGIS.Client.Tasks.QueryEventArgs args)             {                 // Clear previous results                 QueryDataGrid.Visibility = Visibility.Visible;                 GraphicsLayer graphicsLayer = MyMap.Layers["QueryResults"] as GraphicsLayer;                 graphicsLayer.ClearGraphics();                      // Check for new results                  FeatureSet featureSet = args.FeatureSet;                 if (featureSet.Features.Count > 0)                 {                     // Add results to map                     foreach (Graphic resultFeature in featureSet.Features)                     {                         //resultFeature.Symbol = ResultsFillSymbol;                         graphicsLayer.Graphics.Add(resultFeature);                     }                 }                 else                 {                     MessageBox.Show("No features found");                 }             }                  // Notify when query fails             private void GeneralQueryTask_Failed(object sender, TaskFailedEventArgs args)             {                 MessageBox.Show("Query failed: " + args.Error);             }      XAML: Everything in the XAML is set up exactly the same, it just uses different IDs. Any help would be appreciated. If you need more information, please do not hesitate to ask.


title: Career advice needed

text: I am a final year software engineering student, and I'm keen to make a career in GIS. I am a beginner-intermediate at C++. What should I work upon to possess the adequate skills? Also can you recommend some good colleges for a GIS-oriented Postgraduate course?


title: Multiuser Editing with ArcGIS through LAN

text: I am working on an ArcObjects application (C#), which is supposed to carry some customized editing operations, the database is ArcSDE based on Postgresql via (accessible via a LAN network). One feature we're keen on is **multi-user-edition (MUE)** : multiple users should be able to edit the same layer (or objects in a layer) simultaneously. The problem is that the application is frequently (once every week) overwhelmed (congested) thus disabling the MUE and leaving only one user to carry all the edition workload. I am not very sure this is specifically related to the application (though there is a fair likelihood it is), so I had to wonder whether such an issue is experienced by users accessing a DB (via ArcSDE) and using ArcGIS desktop as an editor ; *** in other words, my hope is to get an idea about how solide is ArcGIS when it comes to MUE (how well does it deal with conflicts, scalability and data size when MUEing)***. (Heuristically, if it turns out that ArcGIS is infallible or near-infallible in this regard, I'll then know that I have to optimize somehow the application) Thanks in advance for your feedback.


title: How connect C# to Postgresql in host j.layershift.co.uk

text: I have installed postgresql database in http://postgres- project-1241043.j.layershift.co.uk/ host. I want to connect to the database using C#. I use Npgsql with following connection string.               connectionString = @ "Server = postgres-project-1241043.j.layershift.co.uk, Port = 5432, User Id = postgre; Password = abcdef; Database = dbluanvantn;";      But I am not able to connect to the server and get error:               Npgsql.NpgsqlException: Failed to a connection to 'postgres-project1241043.j.layershift.co.uk'.      Am I using correct connection string?. Help me fix it.


title: Are latitude and longitude independent of each other or correlated?

text: I want to develop a bivariate gaussian kernel in order to handle the random noise in my GPS tracks where each point has three parameters (lat, long and time stamp). I was wondering if latitude and longitude are correlated or independent of each other? My common sense says they are independent because no causal relationship exists between them.


title: ArcGIS Geoprocessing and 32/64 bit architecture Issue

text: I am coding a VB.net console application using ArcObjects to call ArcGIS 10.1 geoprocessing tools. I have several geoprocessing calls that actually work, but that pop up an irritating message box when I call them from code. This is a problem because I need to run my application in batch at the command line and this message box pauses the application while it waits for input. ![enter image description here](http://i.stack.imgur.com/F6C1u.png) The message box appears consistently on certain geoprocessing calls, but not for all calls. The two main culprits are "Calculate Field" and "Raster Calculator". Below is an example of the former:               Dim GP As New Geoprocessor     Dim CalcFieldTool As New CalculateField     CalcFieldTool.in_table = pFClass     CalcFieldTool.field = "Elevation"     CalcFieldTool.expression = "!RASTERVALU!"     CalcFieldTool.expression_type = "PYTHON"     CalcFieldTool.code_block = " "     GP.Execute(CalcFieldTool, Nothing)      The messages does NOT appear if I use VBA syntax for the calls instead of Python. However, my application has to run against ArcServer which does not (and cannot) include VBA, and so I am tied to Python Syntax. Also, the message box does not appear when I run the same geoprocessing tools in ArcGIS Desktop. This code works and the Elevation field is indeed updated with the values from the RASTERVALU field. I am just need to avoid the message box popping up. Also note that this is well tested code and used to work without the popup until recently... perhaps ArcGIS 10.1 service pack 1. (I cannot upgrade to 10.2.) I am now almost certain that this is related to an architecture mis-match somewhere. My computer is 64 bit, ArcGIS is 32 bit, my console application is compiled as 32 bit (not Any CPU) and I have both the original 32 bit version of Python installed as well as the ArcGIS 10.1 Service Pack 1 64 bit background geoprocessing download installed. But I am unsure how to direct my code to use the latter. Suggested work arounds would be appreciated.


title: Best Book or other Resources on the Business Benefits from GIS

text: I know this posting is an excellent source of info for the Top unsung business benefits for GIS. My question expands on that for any industry books related to this subject. I know there have already been several postings for recommended books but I am trying to find books that explicitly explain the value of GIS to a company. In short, I am looking for books that will clearly explain the business benefits (financial and otherwise) that GIS can provide for a company. If there are any resources on the web for this information, please let me know. Thanks, Seth


title: Is editing feature service through the REST API an okay approach?

text: I have a database that receives real-time data from alarms throughout the day, or on request. I'd like to create a web map that displays the "bad" alarms. My current idea is to create a Python script that queries the database every 15 minutes and gets a unique ID and a lat/long of the bad alarms. I will have a feature service setup and use the REST API (through Python) to add or delete alarms based on the previous query. I'll then have a basic web map set up to display this feature service with other related layers. Is this a good approach? This data will only be viewed internally. I've never used the REST API before, but it seemed like a straightforward solution. Thanks for any guidance.


title: How can I speed up slow-loading, template-generated MXDs?

text: I have an ArcMap template (.mxt) file which causes .mxd documents created with it to take 5-10 minutes to load (if it ever loads). Even after stripping out all the data sources and tables, the issue still persists. Once it loads, it handles fine. Any general speculation of what could be causing this? Thanks


title: Reclassify many rasters in ModelBuilder

text: How do I reclassify many rasters in ModelBuilder when each raster has its own unique values? For example, when you create a model in ModelBuilder and you add the Reclassify tool you may want to reclass several rasters with the Iterate Raster tool as an input. However, when you open the Reclass toolbox you choose a value for reclassification and then name it as an output raster. The problem is how can you, when iterating the raster, reclass with several reclassification values for each raster? This is what I am trying to do. * * * @Michael Miles-Stimson Here is the python I copied and paste here.               # -*- coding: utf-8 -*-     # ---------------------------------------------------------------------------     # Reclassify2006.py     # Created on: 2014-06-23 16:22:40.00000     #   (generated by ArcGIS/ModelBuilder)     # Description:      # ---------------------------------------------------------------------------          # Import arcpy module     import arcpy          # Check out any necessary licenses     arcpy.CheckOutExtension("spatial")          # Load required toolboxes     arcpy.ImportToolbox("Model Functions")               # Local variables:     Final_NDVI = "E:\\FY2013_projects\\PMC_project\\USGS EarthExplorer\\2006\\Final_NDVI"     v_outFGDB___Value__new = "E:\\FY2013_projects\\PMC_project\\USGS EarthExplorer\\2006\\Vectors\\NDVI_reclassify.gdb"     v_Name_ = "E:\\FY2013_projects\\PMC_project\\USGS EarthExplorer\\2006\\reclassify\\%Name%"     v_Value__new = "E:\\FY2013_projects\\PMC_project\\USGS EarthExplorer\\2006\\Vectors\\NDVI_reclassify.gdb\\%Value%.shp"     Name = "ast_l1b_00302032006175500_20101030104009_12125_ndvi"     v_Value__new__2_ = "E:\\FY2013_projects\\PMC_project\\USGS EarthExplorer\\2006\\Vectors\\NDVI_reclassify.gdb\\%Value%.shp"     ast_l1b_00302032006175500_20101030104009_12125_ndvi_img = "E:\\FY2013_projects\\PMC_project\\USGS EarthExplorer\\2006\\Final_NDVI\\ast_l1b_00302032006175500_20101030104009_12125_ndvi.img"          # Process: Iterate Rasters     arcpy.IterateRasters_mb(Final_NDVI, "a*", "", "NOT_RECURSIVE")          # Process: Reclassify     arcpy.gp.Reclassify_sa(ast_l1b_00302032006175500_20101030104009_12125_ndvi_img, "Value", "-0.17647059261798859 0.022727273404598236 1;0.022727273404598236 0.10313901305198669 2;0.10313901305198669 0.1428571492433548 3;0.1428571492433548 0.16883116960525513 4;0.16883116960525513 0.19117647409439087 5;0.19117647409439087 0.21538461744785309 6;0.21538461744785309 0.24324324727058411 7;0.24324324727058411 0.29870128631591797 8;0.29870128631591797 0.56451612710952759 9", v_Name_, "DATA")          # Process: Parse Path     arcpy.ParsePath_mb(Name, "NAME")          # Process: Raster to Polygon     tempEnvironment0 = arcpy.env.parallelProcessingFactor     arcpy.env.parallelProcessingFactor = "ast_l1b_00302032006175500_20101030104009_12125_ndvi"     arcpy.RasterToPolygon_conversion(v_Name_, v_Value__new, "SIMPLIFY", "VALUE")     arcpy.env.parallelProcessingFactor = tempEnvironment0          # Process: Add Field     arcpy.AddField_management(v_Value__new, "Acres", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")          # Process: Calculate Field     arcpy.CalculateField_management(v_Value__new__2_, "Acres", "[Shape_Area]", "VB", "")


title: How do I get a json list of REST SOE's supported by a Mapservice?

text: The html returned by a rest mapservice url includes "supported extensions" down at the bottom. However the json returned does not list the extensions. Is there a way to get this list of SOE's in a json format? I'd include a sample url, but I can't find any rest SOE's when I search - does anyone know of a sample that has been published? **Update** Thanks to Bill, here's an ESRI sample that illustrates the issue: The html page returns Supported Extensions, but when formatted as json, it doesn't. **Update 2** ESRI support filed it as an enhance request (NIM062551).


title: How do I find software for my GPS?

text: I recently unearthed a GPS unit I purchased a few years. It's a Garmin GPSmap 60 (picture) and has a usb port on it. Is there a way that I can link it to my mac? I do have windows installed also, and could use that if necessary. I tried plugging it in, and downloading data via google earth pro without success. I also hunted around online to see what I should install but did know what I was searching for exactly (a driver? a software program? ). * * * All suggestions appreciated (and it's for a worthy cause as I hope to start gathering data for OpenStreetMap) - thanks!


title: How to use both Parse Path and Calculate Value in model?

text: Am I doing this wrong? I need to create hillshades for a set of raster, but I want the hillshades to retain the raster's original name. I used parse path for that, however there are names that are more than 13 characters so I used calculate value to edit the names. I used the expression r"%Value%".replace("-",""). It worked, I could see from the output that the value changed, however the hillshade does not adopt the naming even if I used %Value%. ![enter image description here](http://i.stack.imgur.com/K4vOY.png) Here is my model. I preconditioned everything so that the hillshade will be the last process. ![enter image description here](http://i.stack.imgur.com/8jFi4.png)


title: Retrieving polygon shapes from postgis database and display it to the browser

text: I have a codes that will convert polygon shapes that is stored in the database into kml format and count all the rows in it and display it to the browser. models.py #my model in creating fields in my table               from django.db import models     from django.contrib.gis.db import models as gismodels          class Butuan_Parcel(gismodels.Model):      NewPin = models.CharField(max_length=254)     geometry = gismodels.MultiPolygonField(srid=32651)     objects = gismodels.GeoManager()          def __unicode__(self):         return self.NewPin;      views.py               from django.shortcuts import render_to_response     from django.contrib.gis.shortcuts import render_to_kml     from butuan_parcel.models import Butuan_Parcel          def all_kml(request) : #functions that convert to kml format     locations = Butuan_Parcel.objects.kml()     return render_to_kml("gis/kml/placemarks.kml", {'places' : locations})          def map_page(request) : #functions that count rows stored in shapefile table     lcount = Butuan_Parcel.objects.all().count()     return render_to_response('butuan_parcel/map.html', { 'location_count' : lcount})      urls.py               from django.conf.urls import *      from butuan_parcel.views import *      from django.views.static import *     from django.conf import settings         from django.contrib import admin           admin.autodiscover()      urlpatterns = patterns('',     (r'^admin/', include(admin.site.urls)),     (r'^kml/', all_kml), #calling the all_kml functions from views.py      (r'^$', map_page)) #calling map_page functions from views.py      map. hmtl               <html>     <body>     var map, base_layer, kml;      function init() {      map = new OpenLayers.Map ('map');       base_layer = new OpenLayers.Layer.WMS ("Openlayers WMS", "http://labs.metacarta.com/wms/vmapo", { layers: 'basic'});      kml = new Openlayers.Layer.GML("KML", "/kml/", { format: Openlayers.Format.KML }); map.addLayers([base_layer, kml]); map.zoomToMaxExtent(); }      </script>          There should be a total of {{ location_count }} Locations. To see this data in Google Earth, OPen as KML          </body>     </html>      and this is the output in localhost:8000               There should be a total of 105 Locations. To see this data in Google Earth, OPen as KML      now, i want to query the polygon shapes and display it to the browser, not to count anymore. Is there anyone who knows how to retrieve the polygon and display it to the browser.. :) ?? Your help is very much appreciated.. ,..,


title: Any good open source GIS library or API for Mobile

text: I am looking for GIS library for tracking on mobile application similar to ESRI mobile library but should be open source and freely customizable preferable java base. I am making Android app which will be display a Map and I have a webservice which will return a JSON response containing Location values ( Lat, Lon),So I want to show this information on the map and perform some task such as routing, for example:   1. What is the distance between point A to point B   2. When will point A reach the destination thanks


title: Looking for a pythonic way to calculate the length of a WKT linestring

text: I was quite unsatisfied with Calculating Length of Linestrings in WGS84 in Miles. It kept me wondering if there is a more convenient, Pythonic way to calculate the length of a WKT linestring according to a given SRID. I have in mind something like:               srid="WGS84"     line="LINESTRING(3.0 4.0, 3.1 4.1)"     print length(line, srid)      I'm looking for an accurate answer, not `sin\cos` approximations. Any ideas?


title: How to change the datasource with python?

text: How do I change the datasource of all layers in a map document using arcpy and python


title: How to join two LinestringM into one without dropping M value?

text: I have, for example, two Linestrings with M value:   `LINESTRINGM(1 0 2, 2 0 4, 4 0 7)`   and   `LINESTRINGM(4 0 7, 6 0 9, 10 0 10)` How can I merge them into one Linestring and preserve the M value? Expected result: `LINESTRINGM(1 0 2, 2 0 4, 4 0 7, 6 0 9, 10 0 10)`


title: If road lengths are raster layer, how to calculate road lengths for different buffers?

text: A question related calculation of road layer lengths using raster layer for land use regression. If raster layers are used instead of vector layer for calculation of land use area and road lengths, how can we represent the road lengths layer as a raster layer and how can we calculate road lengths in raster layer for different buffers? How can we build a model for this using ModelBuilder?


title: Trouble getting metacarta to work with qgis?

text: I have a basic issue with getting qgis to accept the metacarta wms service and wonder if anyone can guide me. I add the "basic layers" to the map, but still can't see the tiles there when I go to the map. I even do the right click -> zoom to layer extent, but there is no visual data in qgis. Only complete white. I must be missing something basic so I appreciate any help or advice. Thank you! ![alt text](http://i.stack.imgur.com/1f73s.jpg)


title: ArcGIS 10.0 Python SearchCursor file locking

text: I have a script that gets a value from a shapefile's field to return to the user. It seems that only when the arcpy.SearchCursor is called ArcMap locks the file and it doesn't get removed after the script is finished running. To get the lock off I have to close ArcMap. In the script I delete the SearchCursor object after using it as well as the row object. The way I have the script working is that it tries to delete the workspace folder on subsequent runs but can't because of the lock...until of course I close ArcMap. Is there any advice on getting this lock to go away? This seems to be a new occurence for me. Thanks.


title: CartoDB simple process

text: I'm attempting to create a simple process with CartoDB, I have a layer with number of polygons, I want to to run a SQL with a multipoint line to report which polygons were intersected.


title: Remove Collar from ECW File (Orthophoto Mosaic)

text: Maybe i'm missing a basic concept when it comes to manipulating/removing values from an RGB orthophoto mosaic. Is it possible to completely eliminate the collar (black background) from an orthophoto mosaic that is in an ECW format. I was able to do in a TIFF format but the file needs to be in an ECW format as it is very large. I've had limited success, using FME, to accomplish this but the resultant file always has areas that are not removed. Hope all this make sense. I'm open to other formats that will work, but as long as it performes to the same level as an ECW. I have access to FME and ArcInfo.


title: GPS post processing - Differential Correction

text: I have a large set of collected GPS points that I'm pulling into Arcgis from SQL and creating an xy event layer to then map them. The fields for the points are as follows: recorded, lat, long, altitude, speed, hAccuracy, vAccuracy. Overall the points are fairly accurate, but of course there are a few cases where the accuracy is off by quite a bit. I have worked with differential correction before, but it was with the same software/hardware (Trimble/Pathfinder) used to collect the data. Is it possible to correct the points after I have created a shapefile? If not, is it possible to correct with the Pathfinder software them before I bring them into GIS? I also realize that if in fact I can correct them I will have to pull use multiple base files as the GPS points collected are over multiple dates. Any ideas or help is greatly appreciated!


title: Projecting shapefile with transformation using OGR in python

text: I want to change coordinate system of a shapefile from British National Grid to WGS84 (in order to be able to convert it to KML later on). To get the best results of projecting I would like to use OSGB_1936_To_WGS_1984_NGA_7PAR transformation (or later on define different transformation depending on location in UK). The only description of this transformation that I know comes from ArcMap and looks like this: Coordinate Frame - dx=446,000000 dy=-99,000000 dz=54,000000 rx=-0,9450000 ry=-0,261000 rz=0,435000 s=-20,892700 .. and I have no clue how to use it with osr.CoordinateTransformation(BNG, WGS84) Any help highly appreciated!! Jarek


title: Add-Ins Category missing from Customize / Commands / Categories - Can't Use Add-Ins

text: I've created a Add-In button that says "Hello World" when clicked. I want to add button to a toolbar. I've tried creating my own "Category" in the Add-In wizard, no success. I've also accepted the default "Category" called "Add-Ins", no success. The Add-In appears to load but is not assigned to a visible "Category". So, I can't use it. I've run the Add-In via Visual Studio, Installed via the ESRI ArcGIS Add-In Installation Utility, and loaded the file via the Customize / Customize Mode / Commands. What is the issue? Do I need to enable Add-Ins or turn on an extension? Any help is appreciated. Thank you


title: Spatial join - many to one

text: I'm trying to do a spatial join that takes the attributes of multiple polygons and applies them to a single polygon in another layer. The data is currently in a geodatabase, but can be exported to shp if needed. In the example below I want to end up with polygon 1 also having the attributes of the features that intersect it so A, B and C. The output doesn't have to remain as a single feature so the output layer would have features 1A, 1B and 1C. ![Many to one join](http://i.stack.imgur.com/ccUr9.png) I have access to ArcGIS 10.1 or QGIS 2.0 **EDIT: To clarify, the output would need to maintain the geometry of 1 for all 3 output features so that they have the same area and 3 different sets of attributes rather than splitting the feature into 3 parts.**


title: How to do supervised image classification in saga with the help of R script or Python script

text: I am using Ubuntu 12.04 and want to do image classification in SAGA. SAGA has numerous classification algorithms which is not present in GRASS and QGIS. I am new to SAGA, but more accustomed to work on R or python. I installed RSAGA package in R, but could not find any command to do image classification. Please let me know how to do image classification in SAGA through an R/python/shell scripts. Thanks, Karthik K.


title: How to properly establish a PostgreSQL connection using QgsCredentials

text: I'm trying to establish a Postgres connection using **QgsCredentials** as described here: How to access user credentials in database connection? The only problem is that I cannot establish a working connection unless I "physically" add the layer to the map canvas (which I don't want because I only want a connection to the DB and a reference to the layer which will be edited later on in the code) When clicking a button the _run_ function reads:               uri = QgsDataSourceURI()     connInfo = uri.connectionInfo()          (success, user, passwd) = QgsCredentials.instance().get(connInfo, None, None)          if success:         # Isn't this enough to establish the connection and get a valid reference to the layer?         uri.setConnection("HOST", "PORT", "DB_NAME", user, passwd)         uri.setDataSource("SCHEMA", "TABLE_NAME", "the_geom")         LYR = QgsVectorLayer(uri.uri(), "LYR", "postgres")              # Unless I add the following line, I am not able to use properly the object         # (which asks the user the credentials until he/she inserts the right ones and gets         # a working connection with a valid reference to the layer)         QgsMapLayerRegistry.instance().addMapLayer(LYR)      What am I doing wrong? Thanks in advance, m.


title: how do i create clusters of nearby points?

text: I m trying to create clusters of nearby points (which would be the adresses of my customers around a store). Buffer is not an option cauz i need the zone to be limited by customers, so a polygone shape not round I've tried the convex hull but then my attibute table is empty and I need all the associated data inside the shape I've also tried contour but i have 500 stores and i need to do it for each store separately The point to polygone extension considers all the points together Does someone have an idea of how to do this? I've tried looking into the questions but i can't find a specific answer to my problem I'm working on QGIS 2.0.1 Thank you


title: Should I consider any other geocoding services?

text: I need to do geocode lookups so that I can show members of my website other members who are in their vicinity. My approach will be to obtain their latitude and longitude from their mobile device and then do a proximity calculation using the Haversine formula. If they won't allow me to use their location or if they're on a computer, I'll do a geocode lookup using either the postal code or city/region/country they provided when they signed up to obtain their latitude and longitude. I want to obtain global data, i.e. not just North America and Europe, but Asia and South and Central America if possible. My platform is Python and Django. I can't spend a lot of money. I'd also like to be able to obtain some type of support (phone, email, chat). I've done a brief assessment of the following geocode services and am wondering if there's anyone else I should consider?   * Cloudmade: Yes. They appear to meet my requirements.   * MapQuest: Yes. They appear to meet my requirements.   * Google: No. Their terms state that I have to display a map which I don't want to do.   * Yahoo: No. Their business doesn't depend on this service. Who knows how long it will be around?   * Microsoft: No. Bing may get sold. Again, Microsoft doesn't rely on this service so I consider them too risky.   * ArcGIS: No. They meet my requirements but they cost at least $1,500 for their Desktop product.   * deCarta: No. Their service looks good but their site says they don't have geocode data for Africa, Asia, and most of South and Central America.   * Cartodb: Yes. They appear to meet my requirements.   * Texas A&M Geocoding Services: No. Their data appears to be U.S.-centric.   * Data Science Toolkit: No. Support seems questionable.   * NetAcuity: No. This is overkill for my purposes and would likely be far too expensive.   * MaxMind GeoIP: No. This requires IP addresses which aren't always accurate.   * TwoFishes: No. This isn't a commercial product. No support.   * Boundless: No. Another expensive enterprise-grade product.   * Pitney Bowes: No. Enterprise-grade Would Cloudmade and MapQuest be my best bets? Is there anyone else I should be looking at?


title: History of the web map tile

text: I'm looking for a brief history of a 256x256 web map tile, in particular who used it first (came up with the idea) and who popularized it. Thanks!


title: Looking for opensource projects for mobile offline trip planning for public transit and walking routes

text: I'm looking for initiatives pursuing mobile offline trip planning using a combination of public transit data and walking-routes. Thus far, without much success. Osmand, seems to be closest, but the part for offline tripplanning isn't completely funded yet. I may contribute, but haven't got a clue about its viability. I'm looking for any projects pursuing this, or papers (comfortable with reading the academic math-stuff) outlining this. Also roadblocks preventing this (or rather: making quick implementation unlikely), etc. welcome of course. I'm 'only' looking to do offline trip planning within one city at a time, which obviously brings down complexity lot. P.S: adapted question from question posted on opentripplanner-forum: https://groups.google.com/forum/?fromgroups=#!topic/opentripplanner- users/m6PGVRBA0eU


title: Google Maps and Search for KML URL

text: Can anyone tell me how I could add the functionality as seen on google maps that allow users to specify a URL to a KML file in the search bar? I would like to include this in my own embedded map but I'm really lost on how to execute it. I've seen examples for KML URL search on other sites but they don't seem to work. Maybe google has dropped this feature from embedded maps/api3 maps?


title: How to create feature dataset by pgAdmin(postgres) that can by use in sde conection?

text: How to create feature dateset in postgres(pgAdmin) that can by use in sde conection (arcgis)?


title: Generate straight lines between river confluences

text: I have the shapefile containing polylines that represent the river. I need to draw generalized straight line approximations of stream lines. River lines have the following attributes:   * river_id is a unique identifier;   * down_id is the river_id of the downstream line;   * down_node is the downstream node identifier;   * up_node is the upstream node identifier;   * type, indicates whether the line is a source, sink or intermediate line. Any ideas how to do it?


title: How to display custom context menu during create feature task?

text: I have a created a custom context menu using `ICommandBar` and am displaying it on create feature mouse-down right click. I am constructing a point collection on each mouse down and then creating a polyline. For each point I am refreshing using `EditSketch.RefreshSketch()`. On right click, I created a custom context menu to finish the sketch. Sometimes, I get the default "Sketch Tool/Select Elements" context menu, which I do not want to get on right click. The code that I'm using to display the context menu:               public override void OnMouseDown(int Button, int Shift, int X, int Y)     {       if (Button == 2)       {         _objContextMenu.Popup(0, 0);       }     }


title: MapServer and IIS Express configuration

text: I hope you can help me with the following: We need to create a stand alone asp.net based site that will display some maps. We are using IIS Express to run a site and now we need to add the map server part. We have map files from the existing Apache/php & OpenLayers based site and need to add them to our site. My MapServer understanding is basic so any help is appreciated. Also where I can find any good references/tutorials? Thanks in advance.


title: Problem with interpolating a vertical profile in 3D Analyst (ArcMap 10.1)

text: I have a problem when I want to create a vertical profile in 3D Analyst from a DEM file. I want to interpolate the data and show points from 20 to 20 meters, as I specified in the option panel (see figure 1). After I press the button to show the Profile Graph, the X-axis do not show values between 20 meters as I specified in the option panel (see figure 2). What can I do to solve this issue? Note: the Y-axis is displayed correctly in terms of altitude. Datum: WGS_1984 ![enter image description here](http://i.stack.imgur.com/i8Y49.png) ![enter image description here](http://i.stack.imgur.com/ICM4d.png)


title: Will killing ArcMap process while it's copying features corrupt the file geodatabase?

text: So..., I made the unfortunate decision to try and copy/paste 22,000+ address points into the same file geodatabase. It's taking soooooooo long.... So my question is, how much of a chance is there to corrupt a file geodatabase during a copy/paste operation if you kill the ArcMap process midway thru?


title: Upgrade QGIS from 1.8 Lisboa to 2.0?

text: I have QGIS 1.8 and I want to upgrade to 2.0 Dufour... Should I just download the file and launch the installation? Or first uninstall 1.8? Thank you! Lily


title: "Official" definition to hyperspectral sensor

text: I was checking in the literature about a formal definition for hyperspectral sensors, and the most complete definition I found was from The Center of Space Research (@ University of Texas in Austin): > Hyperspectral data sets are generally composed of about 100 to 200 spectral > bands of relatively narrow bandwidths (5-10 nm), whereas, multispectral data > sets are usually composed of about 5 to 10 bands of relatively large > bandwidths (70-400 nm). Is there any consensus about it? Does exist an international body or forum to deal with remote sensing definitions?


title: gdalwarp won't export AAIGrid

text: I'm trying to use gdalwarp to reproject an ascii grid file. According to http://www.gdal.org/formats_list.html this shouldn't be a problem and indeed I can happily use gdal_translate to create grd files (AAIGrid) but not gdalwarp:               gdaltest 14:40:08 $ gdal_translate -of AAIGrid rj_DEM_1978.grd junk.grd     Input file size is 401, 301     0...10...20...30...40...50...60...70...80...90...100 - done.     gdaltest 14:40:25 $ rm junk.*     gdaltest 14:40:41 $ gdalwarp -of AAIGrid rj_DEM_1978.grd junk.grd     Output driver `AAIGrid' not recognised or does not support     direct output file creation.  The following format drivers are configured     and support direct output:       VRT: Virtual Raster       GTiff: GeoTIFF       NITF: National Imagery Transmission Format       HFA: Erdas Imagine Images (.img)       ELAS: ELAS       MEM: In Memory Raster       FITS: Flexible Image Transport System       BMP: MS Windows Device Independent Bitmap       PCIDSK: PCIDSK Database File       ILWIS: ILWIS Raster Map       SGI: SGI Image File Format 1.0       Leveller: Leveller heightfield       Terragen: Terragen heightfield       netCDF: Network Common Data Format       HDF4Image: HDF4 Dataset       ISIS2: USGS Astrogeology ISIS cube (Version 2)       ERS: ERMapper .ers Labelled       RMF: Raster Matrix Format       RST: Idrisi Raster A.1       INGR: Intergraph Raster       GSBG: Golden Software Binary Grid (.grd)       GS7BG: Golden Software 7 Binary Grid (.grd)       PNM: Portable Pixmap Format (netpbm)       ENVI: ENVI .hdr Labelled       EHdr: ESRI .hdr Labelled       PAux: PCI .aux Labelled       MFF: Vexcel MFF Raster       MFF2: Vexcel MFF2 (HKV) Raster       BT: VTP .bt (Binary Terrain) 1.3 Format       LAN: Erdas .LAN/.GIS       IDA: Image Data and Analysis       GTX: NOAA Vertical Datum .GTX       NTv2: NTv2 Datum Grid Shift       CTable2: CTable2 Datum Grid Shift       ADRG: ARC Digitized Raster Graphics       SAGA: SAGA GIS Binary Grid (.sdat)          gdaltest 14:40:52 $ which gdal_translate     /Library/Frameworks/GDAL.framework/Programs/gdal_translate     gdaltest 14:41:05 $ which gdalwarp     /Library/Frameworks/GDAL.framework/Programs/gdalwarp     gdaltest 14:41:10 $       As you can see, the format list is somewhat shorter than normal. From the two `which` commands at the bottom you can also see that both `gdalwarp`and `gdal_translate`are located in the same library. Does anyone have an explanation for this? I can work around it easily enough but this shouldn't be happening, should it?


title: Create Cost Array ignores NoData!

text: After a great deal of research here is what is happening: We select two points on a hydrography map with an island inbetween. It is expected that the CostDistanceWithPolyline would return a line that traverses around the island. We are now finding that the Cost Array that is computed ignores the existance of SOME of our NoData areas, but not all of them! What is wrong here??? We cannot find the issue!!! The hydrography map with island in question: note the value is NoData for a point I chose inside the island geometry. ![Hydrography map with Info display of NoData](http://i.stack.imgur.com/nGPy7.png) The resulting shortest path. Kinda hard for a boat to do that! ![Path thru the island](http://i.stack.imgur.com/dSX7e.png) The cost array and it's value at the island location. ![cost array with Nodata ignored!](http://i.stack.imgur.com/GoiBE.png) The call is so stupid for getting a cost array that I cannot determine how the call could be messed up. As per code that ESRI did for us: we pass the hydrography map, the target point, and call GetCostDistanceFull Here is our call to get the cost distance and backlink arrays.                   private static IGeoDataset[] ComputeFullRaster(IFeatureClass origin, RasterDataset costArray, IGPMessages messages)         {             List<IGeoDataset> results = new List<IGeoDataset>();             IDistanceOp distanceOp = new RasterDistanceOpClass();                  // Create the cost distance             messages.AddMessage("Creating the distance cost raster...");             object objectMissing = Type.Missing;                  results.Add( distanceOp.CostDistanceFull(                 (IGeoDataset)origin, (IGeoDataset)costArray,                  true, false,                 false, ref objectMissing, ref objectMissing) );             messages.AddMessage("Creating the backlink cost raster...");             results.Add(distanceOp.CostDistanceFull(                 (IGeoDataset)origin, (IGeoDataset)costArray,                 false, true,                 false, ref objectMissing, ref objectMissing));             return results.ToArray();         }


title: Calculate difference between values in linear order using ArcGIS

text: The best way to describe this without knowing Python or VB jargon is to add these screenshots. I would like to emulate this Excel calculation in ArcMap10.1. MEAS is the distance along my route and ET_Z is the elevation at the route vertex. I just need the difference from the next point down the line from the previous one. ![Excel calculation](http://i.stack.imgur.com/5XZv1.jpg) I'm not sure the best approach to this, but figured that there may be simple Python solution for this. Thanks in advance!


title: how can i export a shapefile

text: hi I want to select the cities that exist in Morocco and export it to a shape file and I used this code but its not working, it gives me an error: code:               import arcpy     from arcpy import env     arcpy.env.workspace= r"E:\cities"      Mo= r"E:\cities\cities.shp"      rows = arcpy.SearchCursor(Mo,'"CITY_NAME"=\'Rabat\'')     row = rows.next()     row.getValue("CNTRY_NAME")     where='"CNTRY_NAME" ="Morocco"'     arcpy.SelectLayerByAttribute_management("cities","New_SELECTION",where)     arcpy.RefreshActiveView()          error:     Traceback (most recent call last):       File "E:\cities\scripts\selection_villes.py", line 9, in <module>         arcpy.SelectLayerByAttribute_management("cities","New_SELECTION",where)       File "c:\program files\arcgis\desktop10.1\arcpy\arcpy\management.py", line 6461, in SelectLayerByAttribute         raise e     ExecuteError: ERROR 000358: Expression non valide     Échec de l’exécution de (SelectLayerByAttribute).               Échec de l’exécution de (selectionnerVilles).      please help! * * * thank u so much!that was very helpfull.it works now thank u again :)


title: How to find lowest value that is not null using Python in ArcGIS for Desktop?

text: I have four fields in my table, (A,B,C,D). The field type is double with some fields being NULL. I used min(A,B,C,D) to return the smallest value and it works well if there are no NULL value in any of the fields . What would be the easiest way to return the minimum value that is not NULL? I think I found a very simple way: min(filter(None,[!A!,!B!,!C!,!D!]))


title: EPSG 4267 (NAD27) boundaries

text: I need the bounding box coordinates for NAD27 (EPSG/SRID 4267) to check the user's coordinate input. http://spatialreference.org/ref/epsg/4267/ shows the area that doesn't cover North America at all. The EPSG Geodetic Parameter Registry gives the positive longitude as well. What's the trick? What coordinates should I use or where I can find them? Thanks, yup


title: How do you extract values from a map using a transect (just a line) in Idrisi

text: So I have a base map containing Data, I want to make a transect of that map, showing the data of the base map along the transect. I tried making a raster line (each point has a random value) I can now use extract to get the values of the base map along this line. Problem is extract displays them alphabetically. And to give each point of my transect an identifier i used random, so when it get's displayed alphabetically, it loses all meaning. Can anyone help me.


title: openLayers and geoserver: WMSGetFeatureInfo does not return feature beyond wrapdateline

text: I built a map made of several layers. I need to display the layers features fields at the place where the user clicks (see the controler code hereafter). The layers are correctly displayed and go across the dateline. If the user click on a feature, its fields values are written in a div. However, if I click on a feature beyond the dateline, my controler returns nothing. How can I force it (or force geoserver?) to consider clicks beyond the dateline?               [...]     var options = {minResolution:minResolution, maxResolution:maxResolution,             controls:[new OpenLayers.Control.PanZoom(), new OpenLayers.Control.NavToolbar(), layersSwitcher, graticule]};     var map = new OpenLayers.Map("map-id", options);     [...]          var map = new OpenLayers.Map("map-id", options);          var infoGnrl = new OpenLayers.Control.WMSGetFeatureInfo({                    url:'http://mysite.org/geoserver/ocean/wms',                    title:'identify feature by clicking',                    output:'features', infoFormat:'application/vnd.ogc.gml',                    format: new OpenLayers.Format.GML,                    eventListeners: {                     getfeatureinfo: function(event) {                         eezName='';                         lmeName='';                         if (typeof(event.features[0])=='undefined'){document.getElementById('infoGeneral').innerHTML=' '; return};                         for (ii=0; ii< event.features.length; ii++){                             thisEEZNAME=event.features[ii].attributes['eez'];                             thisLMENAME=event.features[ii].attributes['LME_NAME'];                             if (typeof(thisEEZNAME)!='undefined') {if (eezName=='') {eezName=thisEEZNAME} else {eezName+=', '+thisEEZNAME}};                             if (typeof(thisLMENAME)!='undefined') {if (eezName=='') {eezName=thisEEZNAME} else {lmeName+=', '+thisLMENAME}};                         }                       document.getElementById('infoGeneral').innerHTML=eezName;                     }                    }                 });


title: Import CSV file into Postgis DB with Location-field formatting?

text: I am working to import CSV into table in my Postgis Database. First, I created the table               Create Table test (     name varchar,     created_date Date,     location double precision     );      When I excuted the command to               COPY test FROM 'data/file.csv' DELIMITERS ',' CSV HEADER;      I got this error:               ERROR:  invalid input syntax for type double precision: "(40.74053344654042, -74.00778384953068)"     CONTEXT:  COPY test, line 2, column location: "(40.74053344654042, -74.00778384953068)"      What should I make as a type when I create the model ? Float, double precision or what exactly ? Data are similar to (40.74053344654042, -74.00778384953068) UPDATE: CSV File look like this:               name,created_date,location     "Coffe",01/02/2007,(40.74053344654042, -74.00778384953068)     "School",01/02/2007,(40.72324713800021, -74.00495699358042)     "Building 234",01/02/2007,(40.842741313574706, -73.83840584215893)     "Building 4",01/02/2007,(40.842741313574706, -73.83840584215893)     "Building 2435",01/02/2007,(40.75433132495244, -73.99262239963087)     "Building 255",01/02/2007,(40.74482004786735, -73.98511337722212)


title: Add a layer as input with code of OpenLayers / GeoExt using GeoServer

text: I try to develop a GIS app with raster data. I have integrated the code of the following openlayers example in my app. http://dev.openlayers.org/releases/OpenLayers-2.13.1/examples/wps.html I have built my app's interface, using the code from the following GeoExt's example: http://api.geoext.org/1.1/examples/wms-tree.html I have uploaded my raster data on GeoServer. The piece of code from the openlayers' example that add the raster image in the input form is the following:               function addRasterInput(input) {         var name = input.identifier;         var field = document.createElement("input");         field.title = input["abstract"];         var url = window.location.href.split("?")[0];         field.value = url.substr(0, url.lastIndexOf("/")+1) + "data/tazdem.tiff";         document.getElementById("input").appendChild(field);         (field.onblur = function() {             input.reference = {                 mimeType: "image/tiff",                 href: field.value,                 method: "GET"             };         })();     } ` I'm trying to modify the above code so that, when the user check a raster image from the GeoServer Demo WMS, (the folder that retrieve data from an xml GetCapabilities file and it is displayed to the left side of the screen) to be added to the input form. How can I do this ? Thank you in advance, George


title: Why does the time taken for a process or tool to operate differ depending on changes made in a Model, downstream ? (ArcGIS, ModelBuilder)

text: Has anyone experienced this issue, do you have any advice to share? An example situation in a model, Model = Input is a feature class, a field is added, the field is calculated using calculate field. This process takes for example 4 mins. The feature class has 2.9 million records / rows. My issues is - using the exact same input data, same computer, same extent etc and all same other settings in the model, that when additional tools and processes are added to the model downstream , e.g. AFTER this section of the model that now when the model is run these first few steps of the model - the Add Field and Calculate Field now take 11 mins! My question is, does ArcGIS somehow scan the whole model before it is run and does ArcGIS then set memory limits? e.g are different amounts of memory available to the model to use dependent on the tools and process being lined up to run in the model? or is there something else that can impact on tool or process run times, that is caused by changes to the model after the tool or process in question? any ideas or thoughts appreciated, as i am hoping to speed up the models as much as possible.


title: Progress-bar not appearing on Leaflet-MarkerCluster

text:
I am having problems with the progress-bar of leaflet markercluster. I am using it in a geodjango application which uses the django-leaflet app to handle the leaflet requirements. Everything works but the progress bar (it doesn't show up). I am following this example. Here is the template:

     {% leaflet_map "main" callback="main_map_init" %}


title: Where can I get jsde_sdk.jar for ArcSDE 10?

text: I had an environment working with no problems, with `Geonetwork` `harvesting` metadata from `ArcSDE` _9.3_. It was done trhu the following libraries, copied from the ArcSDE installation to the Geonerwork `lib` dir:   * `jsde_sdk.jar`   * `jpe_sdk.jar`   * `concurrent.jar`   * `icu4_3_2.jar` After upgrading to `ArcSDE` _10_ , the harvesting stopped working. It is certainly because I'm using `*.jar`s for version _9.3_ , while having a _v10_ `ArcSDE`, but the new installation doesn't have those `jar`s. Does anyone knows where to get them?


title: What are the Pros and Cons of hstore data type in postgis?

text: PostgreSQL has a module named hstore that implement the use of key/value pairs (it's something like the OpenStreetMap model, but in a single table). What are the Pros and Cons of **hstore** regarding GIS compatibilities, data adding and retrieving?


title: Oracle tracing, bind variables, and ArcSDE

text: Oracle tracing with TKPROF is the recommended way to troubleshoot data access performance problems with SDE/Oracle. Unfortunately, queries that differ in any way are considered "different" and summarized independently. For our long- running geoprocessing task, millions of queries like this are generated:               SELECT A, B, C FROM TABLE WHERE (OBJECTID = 123);      The TKPROF output for this query shows stats for _every single call_ of this type, because `OBJECTID` differs in each case. For example:               SELECT A, B, C FROM TABLE WHERE (OBJECTID = 123);          call     count       cpu    elapsed       disk      query    current        rows     ------- ------  -------- ---------- ---------- ---------- ----------  ----------     Parse        1      0.00       0.00          0          0          0           0     Execute      1      0.00       0.00          0          0          0           0     Fetch        1      0.00       0.00          0          3          0           1     ------- ------  -------- ---------- ---------- ---------- ----------  ----------     total        3      0.00       0.00          0          3          0           1          ...          SELECT A, B, C FROM TABLE WHERE (OBJECTID = 124);          call     count       cpu    elapsed       disk      query    current        rows     ------- ------  -------- ---------- ---------- ---------- ----------  ----------     Parse        1      0.00       0.00          0          0          0           0     Execute      1      0.00       0.00          0          0          0           0     Fetch        1      0.00       0.00          0          3          0           1     ------- ------  -------- ---------- ---------- ---------- ----------  ----------     total        3      0.00       0.00          0          3          0           1          etc...      Millions of stats like this are pretty useless. They need to be summarized for every query that follows the pattern `SELECT A, B, C FROM TABLE WHERE (OBJECTID = x)`. The simple `OBJECTID = 123` filter is generated by `GetRow()`. Here is an example of the C# code driving other problematic queries:               IQueryFilter filter = (IQueryFilter)gdbContext.CreateObject("esriGeodatabase.QueryFilter");     filter.WhereClause = "REF_IN_ID  = " + nodeID + " OR NREF_IN_ID = " + nodeID;     streets.Select(filter, esriSelectionType.esriSelectionTypeSnapshot, esriSelectionOption.esriSelectionOptionNormal, null);      If Oracle bind variables are used, however, the TKPROF output is more helpful. For example:               INSERT INTO TABLE (A, B, C)     VALUES      ( :a1, :a2, :a3)          call     count       cpu    elapsed       disk      query    current        rows     ------- ------  -------- ---------- ---------- ---------- ----------  ----------     Parse      687      0.01       0.01          0          0          0           0     Execute    687      0.08       0.10          4        689       4851         687     Fetch        0      0.00       0.00          0          0          0           0     ------- ------  -------- ---------- ---------- ---------- ----------  ----------     total     1374      0.09       0.12          4        689       4851         687      So the question is: is there any way to force ArcSDE/ArcObjects to use bind variables when constructing Oracle queries? If not, how do you recommend performance tuning for a long-running geoprocessing routine like this?


title: Grass GUI integration

text: Is it possible to embed the GRASS map display to a C# application? I know that a PNG image can suffice but to make the application more interactive it could perhaps be integrated. In essence I would just be re-skinning the GUI of grass and adding some simulation buttons etc.


title: How to Create Cost Surface Raster for Least Cost Path Analysis?

text: I have a DEM raster and would like to carry out a cost path analysis. Among other factors such as slope, proximity rivers etc. I would also like to take elevation into consideration. According to the model I am trying to build (using ArcGIS 10.1), the cost of building at high altitude locations would be much more costly compared to low altitude places. I was wondering how I could create such a cost surface? ![enter image description here](http://i.stack.imgur.com/psd8M.jpg)


title: Converting Lambert Conformal Conic to WGS84 using Proj.4

text: I am attempting to transform data from Lambert Conformal Conic lat/long and I am having some issues with the results. The shapefile I am attemtping to transform can be found here: http://tfr.faa.gov/save_pages/0_8326.shp.zip Since I am using proj.4 I used GDAL/srs to convert the WKT to a proj.4 init string. In order to do that I had to change the WKT projection from "Lambert_Conformal_Conic" to "Lambert_Conformal_Conic_SP2". The WKT I used:               PROJCS["North_America_Lambert_Conformal_Conic",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic_SP2"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-77.03638889],PARAMETER["Standard_Parallel_1",38.69259533705237],PARAMETER["Standard_Parallel_2",39.02625053604863],PARAMETER["Latitude_Of_Origin",40.0],UNIT["Meter",1.0]]      The resulting proj.4 string from the WKT:               "+proj=lcc +lat_1=38.69259533705237 +lat_2=39.02625053604863 +lat_0=40 +lon_0=-77.03638889 +x_0=0 +y_0=0 +ellps=GRS80 +units=m +no_defs"      The proj.4 string I use for WGS84 lat/long:               "+proj=longlat +ellps=WGS84 +no_defs"      When I get a Lambert Conformal Conic point from the shapefile, I transform it to Lat/Long using this line of code:               pj_transform(pj_lcc, pj_latlong, 1, 1, &fX, &fY, NULL );      Where pj_lcc is my Lambert Conformal Conic projection info, pj_latlong is my WGS84 lat/long projection info and fX, fY are the Lambert Conformal Conic point. An example result from the above transform is:               lcc point: {4842.47, -71287.5} --- long, lat: {-1.34356, 0.686926}      The lat, long should be around 38º51'34"N, 77º02'11"W (center of the shape according to http://tfr.faa.gov/save_pages/detail_0_8326.html ) Any clues as to what I am doing wrong?


title: WeedMap pro users - Can weed points and polygons be extracted for use in other GIS clients?

text: Are there any users out there of the WeedMap Pro product? I know a lot of local councils in Australia use this for weed management. I am involved with a regional weeds committee involving a number of councils. Each council is putting forward their weeds infestation layers to produce a region wide map, however the councils are using WeedMap pro and apparently they cannot extract their data for use outside of WeedMap. This doesn't sound right? You would think there has to be some way of extracting the spatial representations of these infestations! I have been in contact with the WeedMap people (because we are actually looking at the product ourselves - its quite cool, but this kind of issue annoys me!) and they said the product uses a SQL server 2005 and stores the geometry as text, not a db specific spatial data type. Also that they don't allow direct access to the database due to security. Is there anyone out there using WeedMap pro that knows a way around this? Even if all I end up with is a list of coordinates for points, and polygon centroids that would be acceptable. Thanks Ando


title: Can I specify a point and get all features within a certain radius of this point?

text: I have recently shifted to QGIS from ArcGIS and am looking for some help/pointers with the following: Does QGIS have a method by which I can specify a lat-lon coordinate (or a utm x,y coordinates, coordinate system not imp), and get all features that lie within a certain radius of this point? I downloaded the 3rd party plugins for QGIS and after looking through all of them, didn't find anything with this functionality (though I could have very well not understood something completely). I am considering writing a plugin in case something doesn't already exist Edit: I looked at the option under: View -> Select -> Select Features By Radius. This seems to be similar in concept to what I am aiming to do (except that the center needs to be specified).


title: How to import/export geographic information (gml or csv file) to geodatabase?

text: I am searching for C#/Java libraries that can import and export geographic information (gml/csv file formats) from and to geodatabases. If not yet available, is there other way to do this task? Any guidance on the right direction is appreciated.


title: How to overlay map layers (shp and csv) in R?

text: This is a crosspost from stackoverflow. I know very little about GIS software and am doing all my mapping in R. Apologies in advance if this is too basic question. Let's say I have two shapefiles from different sources but with different attributes. Say, one is for adminstrative boundaries of Texas (`boundaries.shp`) and the other is for Texas rivers (`rivers.shp`). I also have a third file `towns.csv` showing the locations of towns within the state. After reading in the files, I can overplot town locations on the adimistrative boundaries in the `maptools` package:               plot(boundaries); points(towns$lon, towns$lat)      But how can I overlay all three? Surely there is an easy way of doing this?


title: How to assign a specific IP to Apache Tomcat server

text: I have installed Apache tomcat 6 on windows server 2008, i want to assign a specific IP to apache tomcat server, where i do that so that Apache should listen on that specific IP ?


title: How to measure length of line accurately?

text: In postgis 1.5 for the below given query the answer is 1540651.62674101 meters ,               SELECT ST_Length(ST_GeomFromText('LINESTRING(77.83042478941326      9.98463896136159,78.03469506809812 23.90398909552571 )',4326),true);      but in google earth the same line gives 1538923 meters in its properties, around 0.1% error between the two readings. Which one is accurate, if Google earth reading is accurate, how to get it in postgis1.5 ... any idea... please help.


title: Sort a collection of points based on proximity using pythagorean (2D) distance

text: I would like to present a user with a list of points, sorted by proximity to the user's location. Is there any problem with using the pythagorean formula to calculate the two dimensional distance between the point and the user's location for this? Specifically I am wondering if I could potentially get an incorrect ordering of points, due to the fact that I am not projecting them, or using a formula for three dimensional distance. In my case, the actual distance is not important, just the relative distance of one point vs another. My code in F# is:               query {                         for feature in set do             let x2 = double feature.Longitude             let y2 = double feature.Longitude             let distance = sqrt ((x2 - x1) ** 2.0 + (y2 - y1) ** 2.0)             sortBy distance             select feature }


title: Height of the geoid in ECEF coordinates

text: I know that the relationship of geoid height `N` and ellipsoid height `h` to orthometric height `H` is defined by `H = h - N`, but what if I have some position in Earth-Centered-Earth-Fixed (ECEF) coordinates? Specifically, I am looking to relate the height of the geoid (for some particular position on the geoid) to the earth-centered height. (Using WGS84 here.)


title: Insert points into SQL Server?

text: I am collecting X and Y values from a web service (Twitter) via a python script. In the long run, this will run over a period of months and I intend on stopping at around the 6 million point mark. The original coords im getting are geographic WGS84, but I will need to convert these to projected WGS Web Mercator. Ill be later publishing this table to an ArcGIS Server map service and caching it. This is a personal project to learn python with no deadline and was wondering if it would be a good idea to solely make use of the native spatial types from SQL Server? My current untested plan:   * CREATE a table with SSMS, with a GEOMETRY field setup (and some other attributes)   * In my python script, make use of arcpy or pyproj to convert the lat/lons in WGS84 to WGS84 Web Mercator (or can I avoid this somehow and its all achievable with SQL?)   * Make use of pymssql to INSERT records, and insert the points into the GEOMETRY field in the table. My question is, what would be a good, simple and efficient approach to take a pair of lat/lons in WGS84, and then insert them into a SQL Server table making use of SQL Server spatial types and have a resulting points layer that is in WGS84 Web Mercator, so that I can render/query them in ArcGIS Desktop 10.1? I do have access to arcpy/ArcSDE 10.1 if need be but was hoping to use this as an example of not requiring ArcSDE. Make use of comments if you need more information and I will edit it into my Q.


title: Why does ftools' difference tool fail?

text: I have an area polygon with rocks on it i want to use difference to map the soft (ground, not rocks) cover. to do that i try using the (ftools) vector>geoprocessing>difference getting an empty difference layer. I tried performing this analysis both on windows and linux using shapefiles and sqlite. a zipfile containing the data can be downloaded at:https://www.dropbox.com/s/s1lrepa95lwdm9e/stackexchange.zip I have a few questions: 1\. Can anyone tell me why? 2\. Where can i find the qgis log so i might be able to figure this out 3\. Is there an alternative to ftools to perform such analysis? 4\. it seems to me that some of the ftools function dont work for me, does anyone know of a solution for that? 5\. might this be somehow encoding related (some of the layers were produced in windows some linux) OS: windows xp, linux ubuntu 12.04 qgis: 1.8


title: How to check what application uses my services - arcgis server 9.3

text: I serve a lot of map services via my arcgis server 9.3. Can I determine somehow which webapplication or desktop clients are using or were using my services? I know that I can view apache access log and by this way execute this task, but I think that it is not most elegant solution.


title: How to edit and use text produced by Feature InfoTool?

text: How can I access and use text from Feature InfoTool dialog. I'm working with JMapFrame. I found this snippet:                        InfoTool tool = new InfoTool() {              @Override              public void onReporterUpdated() {                  String text = getTextReporterConnection().getText();                  System.out.println(text);              }          };      but I don't see how to link this to my JMapFrame. So basically I don't know how to get InfoTool from JMapFrame.


title: Question about sizes of input rasters for raster calculation

text: To calculate any formula (conditional statements/relational operation etc) in the raster calculator in ArcGIS using two different sources of raster images, I'd like to know input rasters must be the same size (row, col num) or not. Generally, cell sizes of input rasters are the same but col and row numbers are slightly different. Can the raster calculator operation manage their sizes automatically?


title: MapServer / OpenLayers: Is there a limit of 16 zoom levels?

text: I set up a MapServer (WMS) with a raster basemap like this (Mapfile):               MAP       PROJECTION        "init=epsg:900913"       END       WEB         METADATA           "wms_title" "mymap"           "wms_srs" "epsg:900913"         END       END       LAYER         METADATA           "wms_title" "myraster"           "wms_srs" "epsg:900913"           "wms_enable_request" "*"         END         PROJECTION          "init=epsg:900913"         END         NAME         myraster         DATA         "/path/to/myraster.tif"         STATUS       ON         TYPE         RASTER         PROCESSING   "BANDS=1,2,3"         OFFSITE      0 0 0       END      END      Now I display this map with OpenLayers:               var mylayer = new OpenLayers.Layer.WMS( "My Map", "http://127.0.0.1/cgi-bin/mapserv?map=/path/to/mymapfile.map", {layers:'myraster'} );      That works fine so far. However, there seems to be a limit of 16 zoom levels and I would like to be able to zoom in deeper than that. I tried to add the option `numZoomLevels:22` but that did not have any effect. My raster image should be large enough - opened in an image viewer the resolution is much higher than displayed with OpenLayers (I also tried using a larger raster image). Is there a way to enhance the maximum number of zoom levels? EDIT - added OpenLayers code:               var wgs84 = new OpenLayers.Projection("EPSG:4326");     var sphericalMercator = new OpenLayers.Projection("EPSG:900913");          map = new OpenLayers.Map( 'map', { projection: sphericalMercator });          layer = new OpenLayers.Layer.WMS( "OpenLayers WMS",             "http://127.0.0.1/cgi-bin/mapserv?map=/path/to/mapfile.map",             { layers: 'myraster' } );          map.addLayer(layer);     map.setCenter(new OpenLayers.LonLat(7.7, 48).transform(wgs84,sphericalMercator), 13);


title: Where can I find High Quality Road Network Information

text: I am carrying out a research on how complex city road networks get and I need to get accurate map information. I started with openstreetmaps however I ran into a couple of issues: some of the dumps contain duplicate information - for example two ways which have the same name tag have distinct way ids and this skews my representation. Moreover, when I converted their data into a graph, the graph was not connected, infact about 30% of the roads are not connected to the giant component. Could this be a mistake in my representation or something to do with openstreetmaps? I have also tried dumps from the US geological survey however roads are split into segments without enough information to recombine them again. This makes it impossible to get full road information. Can you help with road network data suggestions?


title: How to activate and deactivate popup feature in OpenLayers?

text: I am a beginner in OpenLayers and I want to activate and deactivate a popup feature, which has to show lat,long of the click, which is initiating the popup [no markers]. I am currently using OpenLayers.Control.WMSGetFeatureInfo and a variable to switch between getinfo and popup, But it has some issues, whenever I am clicking always it is making calls to WMS server, so I want separate control, which I can enable and disable as I need... Any ideas welcome... Thanks in advance...


title: Is there a way to auto-update symbology between .mxds for the same shapefile? (Python?)

text: I have a variety of .mxds for any one project I'm working on, which are all different but all use the same set of base shapefiles for my project. I know I can save a .lyr version of my initial shapefile and import this into subsequent .mxds to retain the original symbology I have set, but if I then change the symbology of the original shapefile, even if I save over the .lyr file with the new symbology, nothing updates in any of my existing files? My question is, is there an easy way to auto-update this symbology, or is it a manual task? If it's a case of manual re-linking then that's a lot of updating I'm going to have to do on a routine basis - so I'm keen to find a simple way of doing this? I have no python experience but if this is the way to do it then I am eager to learn!! In terms of what I have found so far for a Python solution: I found this link here on StackExchange, which I think might be a good starting point in terms of Python, as it seems to be a script to auto-update data sources [ How to update shapefile data source in multiple dataframes and mxds using lyr.replaceDataSource and Python dictionary? ] And also I found this thread on ArcGIS Idea [ http://ideas.arcgis.com/ideaView?id=087300000008FNBAA2 ] albeit dated 2010, which seems to ask the same question as me, and someone has provided a link to some arcpy.mapping sample scripts which again look very similar to what I need, though they seem to do one .mxd at a time which is about the same effort as relinking them manually. I did download the sample python scripts to have a look at, and I've installed Python (command line), but I have no experience with this and am struggling a bit with knowing where to start to get these to run! Any help or advice in terms of a Python solution would be very much appreciated, but I'm not sure if there's a simpler way of going about this [ **EDIT** : which AUTO links rather than re-running a script with every edit] - it seems odd to me that this is not a function built in to ArcGIS? [ **EDIT** : I am using ArcMap 9.2]


title: Openlayer GML Format fail at special character

text:
I am working on openlayer and I want to extract features from gml file but it fails at special character and return 0 features. Please check the sample code -               <!DOCTYPE html>

                                         In above example testsuccess.gml print arlert as feature length 25 and for testfail.gml it fails at special character(``) and return nothing how to solve this. I try to load in QGIS it render successfully but in error log as `OGR[3] error 1: XML Parsing Error: Invalid character reference at line 1, column 13467` I want to extract all features using openlayer , can you help me to solve problem... url for test http://wikisend.com/download/552992/testfail.gml http://wikisend.com/download/969314/testsuccess.gml Thanks Amit


title: Bug in Qgis 2.0 Field Calculator with Values from Joined table

text: I have a data of hospitals. I have two excel sheets, one has no of patients examined in each hospital and other has no. of patients treated in each hospital. I want to find percentage of treated patients to examined. I have joined the two sheets to hospital layer in QGIS. Using field calculator i used this expression to calculate the percentage.               toreal( "EIHSsitewise data on STI treated_Total") / "EIHSsitewise data on STI      examined_Total" *100      But i get null result. I have tried `multiplying by 1.0 instead of using toreal()`. I have found that, when save the layer and then perform the same operation i get the accurate results. So i am assuming that there is a problem on using values from joined tables. Has anyone else faced this problem? Is this a bug? and is there a way to get the results without saving everything in new layer?


title: Bug when viewing data from SQL Server in QGIS

text: I'm having an issue with viewing data in QGIS 2 based on a spatial query in SQL Server 2008 R2. The data appears to load, polygons appear and I can add thematic shading based on data. However when I view the attribute table for the layer in QGIS, it seems like the same value is being replicated over and over again. ![QGIS Data Screengrab](http://i.stack.imgur.com/aKBDE.png) If I add labels, only one label appears despite there being multiple features with different attributes in the original SQL Spatial query. ![QGIS Labels Screengrab](http://i.stack.imgur.com/7eaMl.png) When I export the layer to a shapefile and open it in QGIS, all the data appears correctly and labelling works without any problem. ![enter image description here](http://i.stack.imgur.com/d5PIr.png) Why is the layer imported from SQL Server behaving so strangely?


title: No way to De-authorise ArcGIS Server? How to get license back?

text: At v10 of ArcGIS Server, from what I understand of the new license model, it is not possible to de-authorise ArcGIS Server (like you can with Desktop) Clients have visibility over their licences with the handy Customer Care Portal tool.   \- From here you can see how many licenses you have authorised and how many you still have available. Scenario: \- Purchase 1 license for ArcGIS Server (license level irrelevant) CCP: 0 Auth, 1 Avail) \- Install & Authorise in a dev environment CCP: 1 Auth, 0 Avail) \- Roll over to production, install and attempt to authorise = Fails as no license available = No ability to release license, so have to get in touch with ESRI Inc Customer Care to reset the license This might not be too tricky if your over in States, but what about the rest of the world? Get in touch with distributor --> to get in touch with Redlands/Customer Care --> and then back again = Takes lot of time. Just wondering if I am either missing something or if there is a better workaround to the above? At 9.3.1 and prior, it was possible to install using the same ECP several times, and through the online registration process, if Redlands saw that an ECP with one server license was getting authorised lots of times, then the client would get a phone call.   = I imagine the new model tries to take the work away from ESRIs end, but can be inefficient when trying to roll over.   * I also assume that if my machine died, and I had to reinstall again, still would have to get in touch with Customer Care to reset license?


title: sqlite to spatialite

text: I have a sqlite database with a table that has a column for longitude and one for latitude. Is there a way to make this data base 'spatial'?


title: Converting google map api coordinates into meters

text: How would I go about converting the results of my map bounding box into meters within javascript? For example, if I have the coordinates 51.52366625446907, -0.07467269897460938 decimal degrees, can I project that into something like 555000 155000 meters? (If it is easier, I could try to export my raster from ArcGIS in decimal degree format; I don't know which option would be better).


title: CartoDB - update existing table from external data (i.e. Mysql database)

text: I would like my CartoDB data to be updated, using external data (mysql database, I run and manage locally using php scripts), so that the CartDB table gets updated with the most recent data, whenever the map gets loaded. Is that even possible?


title: Is there an IRC channel for ArcGIS or for GIS itself?

text: Is there an IRC channel or chat where people can talk about GIS systems, or about ArcGIS?


title: Computing radius of gyration in PostGIS

text: The radius of gyration is a measure of the standard deviation of distances between points (and their center of mass), i.e., how frequently and how far points are moving. A formal definition is: ![radius of gyration equation](http://i.stack.imgur.com/9sLoz.png) I have a large table of localization data (x-y coordinates, not a geometry field - might be transformed if faster?). What is the best way to compute the radius of gyration with PostGreSQL / PostGIS? For the center of mass, `ST_Centroid(ST_Union(ST_MakePoint(x,y)))`? Then `ST_Distance`? How to aggregate all this? `ST_Distance(MakePoint(x,y),ST_Centroid(ST_Union(ST_MakePoint(x,y))))` does not work (`ERROR: column "x" must appear in the GROUP BY clause or be used in an aggregate function`). Sorry, maybe a stupid question, I'm not so familiar with PostGreSQL.


title: DEM: How to convert a mosaic dataset to TIN?

text: I merged several DEMs together in on single DEM with mosaic dataset. Now I wanted to convert this mosaic dataset in Vector with TIN. But when I use the tool (3D Analyst Tool) Raster to TIN comes an Error... Have I to convert the mosaic Dataset first to another format? (How?) Or do you have any other ideas to convert mosaic dataset to vector??


title: Rasterization tools for vector datasets with styles

text: I have hundreds of Personal Geodatabases and I want to create one georeferenced TIFF for each one. Each PGDB has the same feature classes ('CotourLines', 'Hidrography', etc) for what I have some layers with the styles. The style (or symbol) for each feature could depend on many fields. And the final TIFF files must have a fixed pixel size. I was searching for tools I can use. From ArcToolBox there are some tools like 'Feature to Raster' but doesn't fit to my requirements. Maybe should I make an automation? Having an MXD containing the layers with the styles, loading the feature classes of each PGDB, and using 'Export Map' to create the TIFF file? This is the only solution I have in mind. I think there should be an straightforward solution I could use. Maybe with an ArcGis extension or an open source tool. I could export the PGDB to ShapeFiles (or a PostGIS database) and have the styles in SLD files.


title: Python in ArcGIS: How to use a variable as part of the output path

text: I´m trying to use a variable as part of the output path of a shutil command. Unfortunately I couldn´t work it out how to do it. Can anybody help?               import arcpy     import shutil     Name = arcpy.GetParameterAsText(0)     shutil.copyfile("C:\\test.txt", "C:\\%Name%.txt")      This question here is similar, but was never answered: ArcGIS 10 Model Builder In-line variable substitution for input data path


title: Can I use in a postgis db within the same table a combination of CRS?

text: I'm experimenting around with my postgis2.0 database. I want to create a basic table to store all my point data and in the process,I realized that the coordinates are in different CRSs. So naturally I'm wondering now if the `the_geom` column can contain a variety of CRSs - or Im confined to a single CRS per table. To visualize it a bit, Is the following possible?               point_id | the_geom      |     1        | (geom_string1) <- CRS EPSG:4326     2        | (geom_string2) <- CRS EPSG:2100


title: Alternative vegetation indices

text: I have been using NDVIs with limited success to identify trees in the central Great Plains region of the USA. The problem I have been encountering is that reflectance from farm fields/pastures have essentially the same spectral signature as the trees I am identifying. Is there a vegetation index which can be generated from 4-band NAIP imagery that can do a better job at isolating tree cover mixed throughout agricultural areas? Perhaps a pre/post processing step may be most effective? Either way, thanks for the advice. ![NDVI Example](http://i.stack.imgur.com/RjUSf.jpg)


title: OpenLayers custom request for WMS GetFeatureInfo

text: In OpenLayers, there is a `WMSGetFeatureInfo` control that can be attached to a layer and listen to `click` or `hover` events on the map, requesting info on the feature being currently clicked (or pointed) at. However, how to make a custom `GetFeatureInfo` request, with the following input?   * layer   * [x,y] mouse position OR [lon,lat] position in the map Is it possible to do this elegantly with the `WMSGetFeatureInfo` control, or do I have to ajax the request myself?