title: What might cause "ERROR 000728: Field <fieldName> does not exist within table" in this code?

text: I have seen ERROR 000728 several times and I could never figure out why it happens. Here is a snippet from my code:               str_Dissolve = arcpy.Dissolve_management(Intersect_str, "in_memory/str_Dissolve", "FID_Con_st", "", "MULTI_PART", "DISSOLVE_LINES")       I do have a field named "FID_Con_st". There is no space in the field name and it doesn't seem to be any refresh bug, because it persists after closing the script. These were possible causes suggested by an earlier Question and an ArcGIS Discussion Form Posting. Are there other possible causes that anyone can suggest?


title: Populating Dropdown based on Geoserver service

text: I am new to web programming here. I have a table which has multiple columns including: Restaurant Name,Restaurant Id, Lat, Lon, City, District, Country. My requirement is to build a web client using OpenLayers and JavaScript and allow users to search a Restaurant. I am using GeoServer to create a layer from this table. My problem is : I have to add dropdowns i.e. one for country, another for District /City based on the country selected. This city/district/country mapping has to come from the table which is in Teradata. Is there any way I can get the data to populate the dropdown from GeoServer services Or do I need to create a new service of my own. If I create my own service, can it be hosted from GeoServer. Any inputs?


title: Convert a FeatureLayer to graphics in Arcobjects

text: I want to convert a Map or FeatureLayer to graphics on Arcobjects 10.1, Any help or code examples **Problem Solved** this is the working code:                   private void Convert_to_graphics( IMapDocument doc)         {                  IMap pmap;             IFeature pFeature;             IElement pElement = null;             IMarkerElement pp = new MarkerElementClass();             IGraphicsContainer pGraphicsContainer;             IActiveView pActiveView;                  pmap = doc.get_Map(0);             pActiveView = doc.ActiveView;              //pGraphicsContainer = doc.ActiveView.GraphicsContainer;                       int[] array = {8 , 7 , 6 , 5 , 1}; //layers indexes ,                   pGraphicsContainer = (IGraphicsContainer)pmap;              foreach (var i in array)             {                           IGeoFeatureLayer fl = (IGeoFeatureLayer)pmap.get_Layer(i);                 IFeatureClass fc = fl.FeatureClass;                      IFeatureCursor fcur =   fc.Search(null, false);                      IFeatureRenderer fr = fl.Renderer;                      pFeature = fcur.NextFeature()  ;                       while (pFeature  != null)                 {                          if (pFeature.Shape.GeometryType == ESRI.ArcGIS.Geometry.esriGeometryType.esriGeometryPoint)                     {                         IMarkerElement pmkElement = new MarkerElementClass();                                  ISymbol sym = fr.get_SymbolByFeature(pFeature);                                  pmkElement.Symbol = (IMarkerSymbol)sym;                                  pElement = (IElement)pmkElement;                               }                     if (pFeature.Shape.GeometryType == ESRI.ArcGIS.Geometry.esriGeometryType.esriGeometryPolyline)                     {                         ILineElement pmkElement = new LineElementClass();                         pmkElement.Symbol = (ILineSymbol)fr.get_SymbolByFeature(pFeature);                         pElement = (IElement)pmkElement;                     }                     if (pFeature.Shape.GeometryType == ESRI.ArcGIS.Geometry.esriGeometryType.esriGeometryPolygon)                     {                         IFillShapeElement pmkElement = new PolygonElementClass();                         pmkElement.Symbol = (IFillSymbol)fr.get_SymbolByFeature(pFeature);                         pElement = (IElement)pmkElement;                     }                               if (pElement != null)                     {                         pElement.Geometry = pFeature.Shape;                              pGraphicsContainer.AddElement(pElement, 0);                     }                          pFeature = fcur.NextFeature();                 }                               }                      pActiveView.Refresh();                   }


title: Do I need an Add-in or a Stand-Alone Application for a custom form?

text: The type of application I want to build for ArcMap is to be run from inside of ArcMap and on a click of a button a custom form appears in front of User in which User enters certain parameters and then by some processing, the application returns resulting image to user in the same ArcMap Document... So will this application be an Add-in or Stand-Alone Application..??


title: How to add a large inset svg or pdf in QGIS 2.3 print composer?

text: I have a 48x36 wall map that will include 5 large inset maps, the largest of which is almost 11x17. Due to a bug with html annotations, I need to export each inset individually then include them as images on the final large wall map. I'd like to include them as vector art (svg or pdf). svg has known clipping problems and does not look good in this case - and I want to see the grid coordinate labels so using a layer level masking scheme won't work. So, I export a pdf which looks perfect in an external pdf viewer. But, trying to include that pdf as an image in the composer just shows a large red X. The implication is that pdf is not a valid format for including images in the composer - is that true? As a third option, I try to export the inset map as a raster image file (png or such) but it won't do it because the resolution and the image size are too big and it gives a memory allocation error. Reasonable but too bad. Any other ideas? The next try will be just writing an image of the html annotations and nothing else, and displaying them on top of all other layers in the composer. Surely there's a better way... Granted, if the html annotation bug (does not respect the 'Lock layers' setting in composer) were fixed, this would all be a moot point, but there's no response on that one currently, I think they have bigger fish to fry at the moment. Thanks in advance for any ideas! **UPDATE:** workaround is working so far but is lengthy: 1) export an svg of just the html annotations (all other map layers off) in just the inset area 2) edit that svg to set fill-opacity to 0 for the relevant objects (easier said than done - largely trial and error) so that the svg can be used as an overlay 3) add the edited svg as an image to the final wall-map composer So far, the above is working for one inset. Next steps to attempt: repeat this for each of the 5 insets, and also for the main map since it has its own annotations that you don't want to show in any of the insets; then make sure to turn off all annotation layers for the final composer rendering. It would be nice if the svg would export with the appropriate transparency so it can be used as an overlay in this manner, but, I'm not even sure how to word that as a feature request. Can anyone help clarify? **UPDATE 2:** (more notes on the workaround) when saving the svg from each inset-only composer:   * check 'background' for the map object, and use some color that will make it easier to search-and-replace lines in the svg file after export (i.e. look for lines with that fill= color and fill-opacity="1", and change to fill-opacity="0"   * when adding the image in the main composer, make sure to uncheck 'Background' so you can see the effective transparency of the svg itself - otherwise you will go on a wild goose chase trying to find the right fill-opacity setting to change!   * when editing the svg, search and replace by color as above, and you will also need to change fill-opacity from "1" to "0" for at least one more object, usually the fifth entry in the svg; if that's not it, use trial-and-error; can just hit enter on the filename of the image in composer to instantly reload and redisplay it   * frame: do not display frame or grid (or coordinate labels) in the overlay svg; display those in the map object in the wall map composer


title: ArcMap runs custom script once and then it crashes on second run

text: So I made a script that takes data from an MSSQL database and inserts it into a featureclass. I added the script to a toolbox and am able to run the tool from within ArcMap 10.1. The problem is that I cannot re-run the tool. When I try and run it a second time, ArcMap crashes. I'll start ArcMap again, run the tool and when I try and run it the second time ArcMap crashes again. By calling `arcpy.AddMessage()` coupled with `time.sleep(5)` I have narrowed the problem down to one line of code: `pyodbc_cursor.execute("{call dbo.GetData(?)}",(DataID))` It seems like even though I `close()` and `del` all my pyodbc variables (connection, cursor, row) after I'm done with them the cursor is causing a problem when the script is executed the second time around. EDIT: This seems like some kind of bug between ArcMap, arcpy and pyodbc


title: Using multiple address locators in ArcGIS

text: I am geocoding millions of addresses for my research (every registered voter in six states). I am using two geocoders: one that is really fast and one that has fewer unmatched addresses when I geocode. Is there a way to use ArcMap or ArcCatalog to use the second address locator on only the unmatched addresses from the first? Right now I am importing the attribute table using the R programming language, exporting only the unmatched addresses as a CSV, re- importing them to ArcGIS, and geocoding, but I figure there has to be a faster way.


title: How to display selected features in ArcGIS Identify Dialog using Visual Studio 2010/ArcObjects?

text: I'm brand new to ArcObjects SDKs and am struggling. I have a Python Add-in performing a query to select records (works great)and now trying to call the identify dialog via an .NET addin button that displays the identify dialog box to show attributes of the selected records. Below is the code I have at this point. I currently have the identify dialog displaying, but no records appearing. I know I need to input the selected records somewhere....but not sure where. Any thoughts would be appreciated. (I'm using Visual Studio/Microsoft Visual Basic 2010 and ArcGIS 10.2.1)               Imports ESRI.ArcGIS.ArcMapUI     Imports ESRI.ArcGIS.Carto          Public Class Identify_Button     Inherits ESRI.ArcGIS.Desktop.AddIns.Button     Dim pMxDoc As IMxDocument     Dim activeView As IMap          Public Sub DoIdentify(ByVal activeView As ESRI.ArcGIS.Carto.IActiveView, ByVal x AsSystem.Int32, ByVal y As System.Int32)     pMxDoc = My.ArcMap.Application.Document     activeView = pMxDoc.FocusMap     If activeView Is Nothing Then         Return     End If          Dim map As ESRI.ArcGIS.Carto.IMap = activeView.FocusMap     Dim identifyDialog As ESRI.ArcGIS.CartoUI.IIdentifyDialog = New ESRI.ArcGIS.CartoUI.IdentifyDialogClass     identifyDialog.Map = map          'Clear the dialog on each mouse click     identifyDialog.ClearLayers()     Dim screenDisplay As ESRI.ArcGIS.Display.IScreenDisplay = activeView.ScreenDisplay          Dim display As ESRI.ArcGIS.Display.IDisplay = screenDisplay ' Implicit Cast     identifyDialog.Display = display          Dim identifyDialogProps As ESRI.ArcGIS.CartoUI.IIdentifyDialogProps = CType(identifyDialog, ESRI.ArcGIS.CartoUI.IIdentifyDialogProps) ' Explicit Cast     Dim enumLayer As ESRI.ArcGIS.Carto.IEnumLayer = identifyDialogProps.Layers     enumLayer.Reset()          Dim layer As ESRI.ArcGIS.Carto.ILayer = enumLayer.Next          Do While Not (layer Is Nothing)         identifyDialog.AddLayerIdentifyPoint(layer, x, y)         layer = enumLayer.Next()     Loop          identifyDialog.Show()          End Sub     Public Sub New()          End Sub          Protected Overrides Sub OnClick()         DoIdentify(activeView, 300, 100)     End Sub          Protected Overrides Sub OnUpdate()         Enabled = My.ArcMap.Application IsNot Nothing     End Sub


title: Python coding question, looping through feature classes, AddField, CalculateField, then errors (might be schema lock)

text: This is not a SDE question. File GDB. Hi, I have a general question on looping through feature classes in a file geodatabase. I have (attached) a python script that adds a ton of fields and calculates them. It seems to run okay for 2-3 feature classes but by the time it hits the 4th or 5th one it often breaks down to Error 99999 or a dreaded schema lock error. Then if I comment out the first 2-3 fc's and run it just for 4 and 5 one at a time it runs fine. Is there some way to 'clear' the loop at the end or should I MakeTableView do everything, export the table, then Delete the TableView? Should I be using environments too instead of using the long paths? I don't have many people look at my code and have no real training so if there are any very very obvious things I'm doing wrong please let me know.               import arcpy, time, datetime, csv     from arcpy import env     env.overwriteOutput = True          calcfiles = [     "E:/SpiderOak/projects/streetview/pedestrian_safety/processing/tables/census_joins.gdb/nyc_0250m_census",     "E:/SpiderOak/projects/streetview/pedestrian_safety/processing/tables/census_joins.gdb/nyc_1000m_census",     "E:/SpiderOak/projects/streetview/pedestrian_safety/processing/tables/census_joins.gdb/phi_0250m_census",     "E:/SpiderOak/projects/streetview/pedestrian_safety/processing/tables/census_joins.gdb/phi_1000m_census",     "E:/SpiderOak/projects/streetview/pedestrian_safety/processing/tables/census_joins.gdb/saj_0250m_census",     "E:/SpiderOak/projects/streetview/pedestrian_safety/processing/tables/census_joins.gdb/saj_1000m_census"     ]          prefld = ["hkq",     "hk1",     "hkq",     "hk1",     "hkq",     "hk1"     ]          print 'Calc Apportioned Census Loop started at this time: ' + time.strftime('%c')           for cfname, prefieldname in zip(calcfiles, prefld):              cenvar = [                 "totpop",         #'!SET001001!*!pctorgarea!', #Total Population                      "sbatot",         #'!SET005001!*!pctorgarea!', #Total Population:                      "sbamal",         #'!SET005002!*!pctorgarea!',   #Male:                 "sbam05",         #'!SET005003!*!pctorgarea!',         # Under 5 Years                 "sbam09",         #'!SET005004!*!pctorgarea!',         #5 to 9 Years                 "sbam14",         #'!SET005005!*!pctorgarea!',         #10 to 14 Years                 "sbam17",         #'!SET005006!*!pctorgarea!',         #15 to 17 Years                 "sbam24",         #'!SET005007!*!pctorgarea!',         #18 to 24 Years                 "sbam34",         #'!SET005008!*!pctorgarea!',         #25 to 34 Years                 "sbam44",         #'!SET005009!*!pctorgarea!',         #35 to 44 Years                 "sbam54",         #'!SET005010!*!pctorgarea!',         #45 to 54 Years                 "sbam64",         #'!SET005011!*!pctorgarea!',         #55 to 64 Years                 "sbam74",         #'!SET005012!*!pctorgarea!',         #65 to 74 Years                 "sbam84",         #'!SET005013!*!pctorgarea!',         #75 to 84 Years                 "sbamel",         #'!SET005014!*!pctorgarea!',         #85 Years and over                      "sbafem",         #'!SET005015!*!pctorgarea!',   #Female:                 "sbaf05",         #'!SET005016!*!pctorgarea!',         #Under 5 Years                 "sbaf09",         #'!SET005017!*!pctorgarea!',         #5 to 9 Years                 "sbaf14",         #'!SET005018!*!pctorgarea!',         #10 to 14 Years                 "sbaf17",         #'!SET005019!*!pctorgarea!',         #15 to 17 Years                 "sbaf24",         #'!SET005020!*!pctorgarea!',         #18 to 24 Years                 "sbaf34",         #'!SET005021!*!pctorgarea!',         #25 to 34 Years                 "sbaf44",         #'!SET005022!*!pctorgarea!',         #35 to 44 Years                 "sbaf54",         #'!SET005023!*!pctorgarea!',         #45 to 54 Years                 "sbaf64",         #'!SET005024!*!pctorgarea!',         #55 to 64 Years                 "sbaf74",         #'!SET005025!*!pctorgarea!',         #65 to 74 Years                 "sbaf84",         #'!SET005026!*!pctorgarea!',         #75 to 84 Years                 "sbafel",         #'!SET005027!*!pctorgarea!',         #85 Years and over                      "ractot",         #'!SET013001!*!pctorgarea!',   #Total Population:                 "racwht",         #'!SET013002!*!pctorgarea!',      #White Alone                 "racblk",         #'!SET013003!*!pctorgarea!',      #Black or African American Alone                 "racnat",         #'!SET013004!*!pctorgarea!',      #American Indian and Alaska Native Alone                 "racasn",         #'!SET013005!*!pctorgarea!',      #Asian Alone                 "racpac",         #'!SET013006!*!pctorgarea!',      #Native Hawaiian and Other Pacific Islander Alone                 "racotr",         #'!SET013007!*!pctorgarea!',      #Some Other Race Alone                 "ractwo",         #'!SET013008!*!pctorgarea!',      #Two or More races                      "histot",      #'!SET014001*!pctorgarea!',    #Total Population                 "hisnoh",      #'!SET014002*!pctorgarea!',      #Not Hispanic or Latino:                 "hisnwh",      #'!SET014003*!pctorgarea!',          #White Alone                 "hisnbl",      #'!SET014004*!pctorgarea!',          #Black or African American Alone                 "hisnna",      #'!SET014005*!pctorgarea!',          #American Indian and Alaska Native Alone                 "hisnas",      #'!SET014006*!pctorgarea!',          #Asian Alone                 "hisnpa",      #'!SET014007*!pctorgarea!',          #Native Hawaiian and Other Pacific Islander Alone                 "hisntr",      #'!SET014008*!pctorgarea!',          #Some Other Race Alone                 "hisntw",      #'!SET014009*!pctorgarea!',          #Two or More races                 "hisyes",      #'!SET014010*!pctorgarea!',       #Hispanic or Latino:                 "hisywh",      #'!SET014011*!pctorgarea!',          #White Alone                 "hisybl",      #'!SET014012*!pctorgarea!',          #Black or African American Alone                 "hisyna",      #'!SET014013*!pctorgarea!',          #American Indian and Alaska Native Alone                 "hisyas",      #'!SET014014*!pctorgarea!',          #Asian Alone                 "hisypa",      #'!SET014015*!pctorgarea!',          #Native Hawaiian and Other Pacific Islander Alone                 "hisytr",      #'!SET014016*!pctorgarea!',          #Some Other Race Alone                 "hisytw",      #!SET014017*!pctorgarea!',          #Two or More races                      "edutot",         #'!SET025001*!pctorgarea!',   #Population 25 Years and over:                 "edulhs",         #'!SET025002*!pctorgarea!',      #Less Than High School                 "eduhsg",         #'!SET025003*!pctorgarea!',      #High School Graduate (includes equivalency)                 "edusoc",         #'!SET025004*!pctorgarea!',      #Some college                 "edubac",         #'!SET025005*!pctorgarea!',      #Bachelor's degree                 "edumas",         #'!SET025006*!pctorgarea!',      #Master's degree                 "eduprf",         #'!SET025007*!pctorgarea!',      #Professional school degree                 "edudoc",         #'!SET025008*!pctorgarea!',      #Doctorate degree                      "emptot",         #'!SET037001*!pctorgarea!',   #Civilian Population In Labor Force 16 Years And Over:                 "empemp",         #'!SET037002*!pctorgarea!',      #Employed                 "empune",         #'!SET037003*!pctorgarea!',      #Unemployed                      "medhhi",         #'!SET057001!',                     #Median household income (In 2010 Inflation Adjusted Dollars)                      "pastot",         #'!SET080001*!pctorgarea!',   #Households:                 "paswit",         #'!SET080002*!pctorgarea!',      #With public assistance income                 "pasnpu",         #'!SET080003*!pctorgarea!',      #No public assistance income                      "housun",         #'!SET093001*!pctorgarea!',    #Housing units                      "tentot",         #'!SET094001*!pctorgarea!',   #Occupied Housing Units:                 "tenown",         #'!SET094002*!pctorgarea!',      #Owner Occupied                 "tenren",         #'!SET094003*!pctorgarea!',      #Renter Occupied                      "occtot",         #'!SET095001*!pctorgarea!',   #Housing units:                 "occocc",         #'!SET095002*!pctorgarea!',     #Occupied                 "occvac",         #'!SET095003*!pctorgarea!',     #Vacant                      "grttot",         #'!SET102001*!pctorgarea!',   #Renter-occupied housing units with cash rent:                 "grtl03",         #'!SET102002*!pctorgarea!',      #Less than $300                 "grtl06",         #'!SET102003*!pctorgarea!',      #$300 to $599                 "grtl08",         #'!SET102004*!pctorgarea!',      #$600 to $799                 "grtl10",         #'!SET102005*!pctorgarea!',      #$800 to $999                 "grtl12",         #'!SET102006*!pctorgarea!',      #$1,000 to $1,249                 "grtl15",         #'!SET102007*!pctorgarea!',      #$1,250 to $1,499                 "grtl20",         #'!SET102008*!pctorgarea!',      #$1,500 to $1,999                 "grta20",         #'!SET102009*!pctorgarea!',      #$2,000 or More                      "rntmed",         #'!SET104001!',   #Median Gross Rent                      "povtot",         #'!SET113001!*!pctorgarea!',      #   Families:                 "povpov",         #'!SET113002!*!pctorgarea!',      #      Income in 2010 below poverty level:                 "povcwc",         #'!SET113003!*!pctorgarea!',      #         Married Couple Family: With Related Child Living  Bellow Poverty Level                 "povcnc",         #'!SET113004!*!pctorgarea!',      #         Married Couple Family: No related children under 18 Years                 "povmal",         #'!SET113005!*!pctorgarea!',      #         Male Householder, no wife present:                 "povmwc",         #'!SET113006!*!pctorgarea!',      #            With related children under 18 Years                 "povmnc",         #'!SET113007!*!pctorgarea!',      #            No related children under 18 Years                 "povfem",         #'!SET113008!*!pctorgarea!',      #         Female Householder, no husband present:                 "povfwc",         #'!SET113009!*!pctorgarea!',      #            With related children under 18 Years                 "povnoc",         #'!SET113010!*!pctorgarea!',      #            No related children under 18 Years                 "povabv",         #'!SET113011!*!pctorgarea!',      #      Income In 2010 at or above poverty level                      "trnw16",         #'!SET128001!*!pctorgarea!',      #   Workers 16 Years and over:                 "trncar",         #'!SET128002!*!pctorgarea!',      #      Car, truck, or van                 "trntax",         #'!SET128003!*!pctorgarea!',      #      Public transportation (Includes Taxicab)                 "trnmot",         #'!SET128004!*!pctorgarea!',      #      Motorcycle                 "trnbik",         #'!SET128005!*!pctorgarea!',      #      Bicycle                 "trnwlk",         #'!SET128006!*!pctorgarea!',      #      Walked                 "trnotr",         #'!SET128007!*!pctorgarea!',      #      Other means                 "trnhom",         #'!SET128008!*!pctorgarea!',      #      Worked at home                      "timw16",         #'!SET129001!*!pctorgarea!',      #   Workers 16 Years and over:                 "timwka",         #'!SET129002!*!pctorgarea!',      #      Did not work at home:                 "timu10",         #'!SET129003!*!pctorgarea!',      #         Less than 10 minutes                 "timu20",         #'!SET129004!*!pctorgarea!',      #         10 to 19 minutes                 "timu30",         #'!SET129005!*!pctorgarea!',      #         20 to 29 minutes                 "timu40",         #'!SET129006!*!pctorgarea!',      #         30 to 39 minutes                 "timu60",         #'!SET129007!*!pctorgarea!',      #         40 to 59 minutes                 "timu90",         #'!SET129008!*!pctorgarea!',      #         60 to 89 minutes                 "tim90p",         #'!SET129009!*!pctorgarea!',      #         90 or More minutes                 "timhom",         #'!SET129010!*!pctorgarea!',      #      Worked at home                      "nattot",         #'!SET133001!*!pctorgarea!',      #   Total Population:                 "natnat",         #'!SET133002!*!pctorgarea!',      #      Native Born                 "natfor",         #'!SET133003!*!pctorgarea!',      #      Foreign Born:                 "natntz",         #'!SET133004!*!pctorgarea!',      #         Naturalized Citizen                 "natnoc",         #'!SET133005!*!pctorgarea!',      #         Not a Citizen                      "ownocc",         #'!SET142001!*!pctorgarea!',      #   Owner-occupied housing units                 ]                varexp = [                 '!SET001001!*!pctorgarea!', #Total Population                      '!SET005001!*!pctorgarea!', #Total Population:                      '!SET005002!*!pctorgarea!',   #Male:                 '!SET005003!*!pctorgarea!',         # Under 5 Years                 '!SET005004!*!pctorgarea!',         #5 to 9 Years                 '!SET005005!*!pctorgarea!',         #10 to 14 Years                 '!SET005006!*!pctorgarea!',         #15 to 17 Years                 '!SET005007!*!pctorgarea!',         #18 to 24 Years                 '!SET005008!*!pctorgarea!',         #25 to 34 Years                 '!SET005009!*!pctorgarea!',         #35 to 44 Years                 '!SET005010!*!pctorgarea!',         #45 to 54 Years                 '!SET005011!*!pctorgarea!',         #55 to 64 Years                 '!SET005012!*!pctorgarea!',         #65 to 74 Years                 '!SET005013!*!pctorgarea!',         #75 to 84 Years                 '!SET005014!*!pctorgarea!',         #85 Years and over                      '!SET005015!*!pctorgarea!',   #Female:                 '!SET005016!*!pctorgarea!',         #Under 5 Years                 '!SET005017!*!pctorgarea!',         #5 to 9 Years                 '!SET005018!*!pctorgarea!',         #10 to 14 Years                 '!SET005019!*!pctorgarea!',         #15 to 17 Years                 '!SET005020!*!pctorgarea!',         #18 to 24 Years                 '!SET005021!*!pctorgarea!',         #25 to 34 Years                 '!SET005022!*!pctorgarea!',         #35 to 44 Years                 '!SET005023!*!pctorgarea!',         #45 to 54 Years                 '!SET005024!*!pctorgarea!',         #55 to 64 Years                 '!SET005025!*!pctorgarea!',         #65 to 74 Years                 '!SET005026!*!pctorgarea!',         #75 to 84 Years                 '!SET005027!*!pctorgarea!',         #85 Years and over                      '!SET013001!*!pctorgarea!',   #Total Population:                 '!SET013002!*!pctorgarea!',      #White Alone                 '!SET013003!*!pctorgarea!',      #Black or African American Alone                 '!SET013004!*!pctorgarea!',      #American Indian and Alaska Native Alone                 '!SET013005!*!pctorgarea!',      #Asian Alone                 '!SET013006!*!pctorgarea!',      #Native Hawaiian and Other Pacific Islander Alone                 '!SET013007!*!pctorgarea!',      #Some Other Race Alone                 '!SET013008!*!pctorgarea!',      #Two or More races                      '!SET014001!*!pctorgarea!',    #Total Population                 '!SET014002!*!pctorgarea!',      #Not Hispanic or Latino:                 '!SET014003!*!pctorgarea!',          #White Alone                 '!SET014004!*!pctorgarea!',          #Black or African American Alone                 '!SET014005!*!pctorgarea!',          #American Indian and Alaska Native Alone                 '!SET014006!*!pctorgarea!',          #Asian Alone                 '!SET014007!*!pctorgarea!',          #Native Hawaiian and Other Pacific Islander Alone                 '!SET014008!*!pctorgarea!',          #Some Other Race Alone                 '!SET014009!*!pctorgarea!',          #Two or More races                 '!SET014010!*!pctorgarea!',       #Hispanic or Latino:                 '!SET014011!*!pctorgarea!',          #White Alone                 '!SET014012!*!pctorgarea!',          #Black or African American Alone                 '!SET014013!*!pctorgarea!',          #American Indian and Alaska Native Alone                 '!SET014014!*!pctorgarea!',          #Asian Alone                 '!SET014015!*!pctorgarea!',          #Native Hawaiian and Other Pacific Islander Alone                 '!SET014016!*!pctorgarea!',          #Some Other Race Alone                 '!SET014017!*!pctorgarea!',          #Two or More races                      '!SET025001!*!pctorgarea!',   #Population 25 Years and over:                 '!SET025002!*!pctorgarea!',      #Less Than High School                 '!SET025003!*!pctorgarea!',      #High School Graduate (includes equivalency)                 '!SET025004!*!pctorgarea!',      #Some college                 '!SET025005!*!pctorgarea!',      #Bachelor's degree                 '!SET025006!*!pctorgarea!',      #Master's degree                 '!SET025007!*!pctorgarea!',      #Professional school degree                 '!SET025008!*!pctorgarea!',      #Doctorate degree                      '!SET037001!*!pctorgarea!',   #Civilian Population In Labor Force 16 Years And Over:                 '!SET037002!*!pctorgarea!',      #Employed                 '!SET037003!*!pctorgarea!',      #Unemployed                      '!SET057001!',                      #Median household income (In 2010 Inflation Adjusted Dollars)                      '!SET080001!*!pctorgarea!',   #Households:                 '!SET080002!*!pctorgarea!',      #With public assistance income                 '!SET080003!*!pctorgarea!',      #No public assistance income                      '!SET093001!*!pctorgarea!',    #Housing units                      '!SET094001!*!pctorgarea!',   #Occupied Housing Units:                 '!SET094002!*!pctorgarea!',      #Owner Occupied                 '!SET094003!*!pctorgarea!',      #Renter Occupied                      '!SET095001!*!pctorgarea!',   #Housing units:                 '!SET095002!*!pctorgarea!',     #Occupied                 '!SET095003!*!pctorgarea!',     #Vacant                      '!SET102001!*!pctorgarea!',   #Renter-occupied housing units with cash rent:                 '!SET102002!*!pctorgarea!',      #Less than $300                 '!SET102003!*!pctorgarea!',      #$300 to $599                 '!SET102004!*!pctorgarea!',      #$600 to $799                 '!SET102005!*!pctorgarea!',      #$800 to $999                 '!SET102006!*!pctorgarea!',      #$1,000 to $1,249                 '!SET102007!*!pctorgarea!',      #$1,250 to $1,499                 '!SET102008!*!pctorgarea!',      #$1,500 to $1,999                 '!SET102009!*!pctorgarea!',      #$2,000 or More                      '!SET104001!',                      #Median Gross Rent                      '!SET113001!*!pctorgarea!',      #   Families:                 '!SET113002!*!pctorgarea!',      #      Income in 2010 below poverty level:                 '!SET113003!*!pctorgarea!',      #         Married Couple Family: With Related Child Living  Bellow Poverty Level                 '!SET113004!*!pctorgarea!',      #         Married Couple Family: No related children under 18 Years                 '!SET113005!*!pctorgarea!',      #         Male Householder, no wife present:                 '!SET113006!*!pctorgarea!',      #            With related children under 18 Years                 '!SET113007!*!pctorgarea!',      #            No related children under 18 Years                 '!SET113008!*!pctorgarea!',      #         Female Householder, no husband present:                 '!SET113009!*!pctorgarea!',      #            With related children under 18 Years                 '!SET113010!*!pctorgarea!',      #            No related children under 18 Years                 '!SET113011!*!pctorgarea!',      #      Income In 2010 at or above poverty level                      '!SET128001!*!pctorgarea!',      #   Workers 16 Years and over:                 '!SET128002!*!pctorgarea!',      #      Car, truck, or van                 '!SET128003!*!pctorgarea!',      #      Public transportation (Includes Taxicab)                 '!SET128004!*!pctorgarea!',      #      Motorcycle                 '!SET128005!*!pctorgarea!',      #      Bicycle                 '!SET128006!*!pctorgarea!',      #      Walked                 '!SET128007!*!pctorgarea!',      #      Other means                 '!SET128008!*!pctorgarea!',      #      Worked at home                      '!SET129001!*!pctorgarea!',      #   Workers 16 Years and over:                 '!SET129002!*!pctorgarea!',      #      Did not work at home:                 '!SET129003!*!pctorgarea!',      #         Less than 10 minutes                 '!SET129004!*!pctorgarea!',      #         10 to 19 minutes                 '!SET129005!*!pctorgarea!',      #         20 to 29 minutes                 '!SET129006!*!pctorgarea!',      #         30 to 39 minutes                 '!SET129007!*!pctorgarea!',      #         40 to 59 minutes                 '!SET129008!*!pctorgarea!',      #         60 to 89 minutes                 '!SET129009!*!pctorgarea!',      #         90 or More minutes                 '!SET129010!*!pctorgarea!',      #      Worked at home                      '!SET133001!*!pctorgarea!',      #   Total Population:                 '!SET133002!*!pctorgarea!',      #      Native Born                 '!SET133003!*!pctorgarea!',      #      Foreign Born:                 '!SET133004!*!pctorgarea!',      #         Naturalized Citizen                 '!SET133005!*!pctorgarea!',      #         Not a Citizen                      '!SET142001!*!pctorgarea!',      #   Owner-occupied housing units                 ]                for fname, expre in zip(cenvar, varexp):                 try:                     #arcpy.AddField_management(cfname, fname, "FLOAT")                     #arcpy.CalculateField_management(cfname, fname, expre, "PYTHON")                     arcpy.AddField_management(cfname, prefieldname+fname, "FLOAT")                     arcpy.CalculateField_management(cfname, prefieldname+fname, expre, "PYTHON")                       # this is with the prefieldname gonna add that later, don't neeed now                     print 'Loop ' + prefieldname + " - field: " + fname + ' done at this time: ' + time.strftime('%c')                  except:                     print 'Script for ' + cfname + ' screwed up at this time: ' + time.strftime('%c')           print 'Script ended at this time: ' + time.strftime('%c')


title: What Python package to use to calculate midpoint?

text: I asked a question earlier about how to calculate a midpoint. I was wondering if there is a Python library that does this alreay?


title: Data management in GRASS from command line?

text: I'm trying to import a large set of raster files into GRASS. The GUI import tool works fine for a handful of files, but it appears to require that I click the "Layer id" checkbox for each individual file that I want to import, and that doesn't work well at all for large sets of file. So my question is, how do you manage data in GRASS from the command line? I realize that this is a broad topic, so please feel free simply to direct me toward the documentation for the relevant tools. I'd also be happy with an approach using the GUI tool in this case, but as I'd prefer a command-line approach because I'm just learning GRASS and need to script almost everything I'm doing, and have read that the Python GUI for GRASS has issues with Mac OS 10.8 (I haven't run into those issues, yet). UPDATE: It does not appear that a GUI solution will work in this case. Having clicked the checkboxes for about 700 files that I wanted to import, the process does not appear to have succeeded due to a "to many files open" error. The tail end of the log is given below.               startupinfo, creationflags)       File "/System/Library/Frameworks/Python.framework/Versions     /2.6/lib/python2.6/subprocess.py", line 616, in __init__          errread, errwrite) = self._get_handles(stdin, stdout,     stderr)       File "/System/Library/Frameworks/Python.framework/Versions     /2.6/lib/python2.6/subprocess.py", line 978, in _get_handles          c2pread, c2pwrite = os.pipe()     OSError     :     [Errno 24] Too many open files


title: Which census geography boundaries do congressional districts preserve?

text: Tried searching for this information elsewhere and couldn't find a clear answer in the TIGER documentation, so I thought I'd turn to you good folks. Obviously, congressional districts can cut across counties and 'cities', but which census geography boundaries do they preserve (if any)? Tracts? Block Groups? Blocks? I suppose I'm most interested in learning what is the largest type of census geography whose boundaries are preserved by congressional districts in every case. Thanks in advance!


title: Mapserver and Tilecache configuration issue

text: Currently I am using **mapnik2** with **tilecache** for map generation. As there is shift in some of the tiles generated,I am trying to configure **mapserver** with tilecache for caching tiles. I have installed **mapserver 6.2.1** and **tilecache 2.11**. In the **tilecache.cfg** file , I configure as below to access map data:               [Layername]     type=MapServerLayer     bbox=255923,904698,1050731,1437039     data_extent=255923,904698,1050731,1437039     layers=layer1,     mapfile=/path to mapfile/     srs=EPSG:32643      When I check if tiles are generated by using the url: http://tilecache.org/docs/README.html an xml file is generated with error being,               <?xml version='1.0' encoding="ISO-8859-1" standalone="no" ?>     <!DOCTYPE ServiceExceptionReport SYSTEM "http://schemas.opengis.net/wms/1.1.1/exception_1_1_1.dtd">     <ServiceExceptionReport version="1.1.1">     <ServiceException code="LayerNotDefined">     msWMSLoadGetMapParams(): WMS server error. Invalid layer(s) given in the LAYERS parameter. A layer might be disabled for this request. Check wms/ows_enable_request settings.     </ServiceException>     </ServiceExceptionReport>     Cache-Control: max-age=10, must-revalidate     Content-type: text/plain          An error occurred: msWMSLoadGetMapParams(): WMS server error. Invalid layer(s) given in the LAYERS parameter. A layer might be disabled for this request. Check wms/ows_enable_request settings.       File "/var/www/tilecache/TileCache/Service.py", line 343, in cgiHandler         format, image = service.dispatchRequest( params, path_info, req_method, host )       File "/var/www/tilecache/TileCache/Service.py", line 208, in dispatchRequest         return self.renderTile(tile, params.has_key('FORCE'))       File "/var/www/tilecache/TileCache/Service.py", line 138, in renderTile         data = layer.render(tile, force=force)       File "/var/www/tilecache/TileCache/Layer.py", line 444, in render         return self.renderTile(tile)       File "/var/www/tilecache/TileCache/Layers/MapServer.py", line 50, in renderTile         wms.loadOWSParameters(req)       File "/usr/lib/python2.7/dist-packages/mapscript.py", line 1770, in loadOWSParameters         def loadOWSParameters(self, *args): return _mapscript.mapObj_loadOWSParameters(self, *args)      Please help me in solving this issue.Thanks in advance.


title: How to use DateDiff() Function across day transition in ArcGIS Field Calculator?

text: I am trying to calculate the time difference between two fields using the DateDiff() function in minutes. I have successfully calculated the times that fall within the same day but I am having trouble with the times that are occurring in day transition (after 11:59 PM). For example: My date fields are as such, MM\dd\yyyy hh:mm:ss tt So between these two times 01\17\2007 11:57:08 PM and 01\18\2007 12:15:01 AM The field calculator would apply a huge negative number. I am thinking a conditional statement is needed for basic arithmetic but am at a loss. Please be as detailed as possible. I am not a developer.


title: Is there general Earth Projection in ArcGIS for Desktop?

text: I'm starting a new personal project and I need a little bit of helping on setting it up. First off I'm trying to determine the area of a rectangle that is 1 mile in width and the circumference of the earth minus a mile in length(like a great circle minus a mile). I understand I can solve this with an equation but I want to map it. Anyways I don't know what projection I need for this and how to get started. Is there a general earth projection in ArcMap? Do I even need a projection if it is the whole earth? Do I need to worry about equal area distortion? and etc.


title: How can I utilize NumPy arrays to optimize big data geoprocessing?

text: I'm interested in learning how to utilize NumPy arrays to optimize geoprocessing. Much of my work involves "big data", where geoprocessing often takes days to accomplish certain tasks. Needless to say, I am very interested in optimizing these routines. ArcGIS 10.1 has a number of NumPy functions that can be accessed via arcpy, including:   1. NumPyArrayToFeatureClass (arcpy.da)   2. RasterToNumPyArray (arcpy)   3. TableToNumPyArray (arcpy.da) For example purposes, let's say I want to optimize the following processing intensive workflow utilizing NumPy arrays: ![enter image description here](http://i.stack.imgur.com/bMVvx.jpg) The general idea here is that there are a huge number of vector-based points that move through both vector and raster-based operations resulting in a binary integer raster dataset. **How could I incorporate NumPy arrays to optimize this type of workflow?**


title: ArcObjects: Equivalent/Similar object to old Avenue ObjectTag from ArcView 3.x?

text: I have annotation where a single annotation feature consists of one or more TextElements and RectangleElements grouped together in a GroupElement and then stored in the annotation feature. To edit the contents of the annotation feature later, it would be nice to be able to know which elements are related to each other. The elements have natural groupings, but those groups cannot be grouped together and then grouped into the final GroupElement for storage. Due to a bug, the ESRI optimized map service cannot handle nested groups in annotation. (See this posting.) So all elements just have to be thrown in together. An identifier object, holding, say, a group name string, attached to each element would make identification of elements belonging to each subgroup much easier to identify. The Avenue ObjectTag object in old ArcView GIS 3.x (I guess it was a property) used to be great for such things. You could stick a tag on just about any object in the whole model, including dialog controls. Really handy. I haven't thought about this -- until recently -- since the old Avenue days, but it has now come to my consciousness again. I haven't been able to find anything like the old ObjectTag in ArcObjects. Any ideas?


title: Calculate line distances inside of polygon

text: I'm working with two data sets:   * NHD polyline stream layer for my entire state   * HUC Watershed boundary Polygons for the entire state I am looking for a way to calculate the total length of streams within each watershed boundary polygon, and would ideally like to have each polygon then attributed with the sum of stream miles within its boundary. Can anyone help me?


title: After SDE Conversion, cannot open our GP Services to republish them

text: We have several services that were originally published against a local file- based database (with ArcMap 10.1). We spent the past week converting the data into an SDE database and we are running under Arcmap 10.2.2. Now it is my turn to convert the services to use the database and get them published to the Arc Server. All of the internal code just calls one of two methods and those all got changed. In Arcmap I removed all of the folder shares, renamed the gdb local database to hide it, and pointed Arcmap at the SDE server. Then I made that database the default. Next I rebuilt our DLL, reregistered it to the Desktop, deleted all of the old service toolbox entries and recreated them. I should specify here that almost every one of our services specify schema names in the database for our parameters in the ParameterInfo method. When I go to publish the services the execute window does not open except for our one service that does not use schema's to define the parameters. The implication is that ArcMap cannot find the schemas in the database. Yes I've tried with our schema names simple (eg. Simple_Polygon) as well as fully qualified (eg. OurServer.DBO.Simple_Polygon) but neither attempt has been successful. I cannot see any other thing I need to do in order to get ArcMap to see our schemas so that I can publish the service. What am I missing??? I cannot find any ESRI documentation that talks about this aspect of publishing a geoprocessing service.


title: QGIS says "osgeo [python-gdal]" isn't installed on Mac

text: I know similar questions have been asked (like this one) -- but I'm encountering a unique problem. I installed GDAL complete 1.11 from KyngChaos, but when I open QGIS, I get the following error message: ![screen shot](http://i.stack.imgur.com/wpYYt.png) Interestingly enough, when I look at the Libraries directory in Finder, there isn't a subfolder called Frameworks. Which tells me something went wrong. Any help would be greatly appreciated. I'm running Mac OS X 10.9.2.


title: Is there a way to embed the filter UI into a map (for the user to use)?

text: Is there a way to embed the dynamic filter UI into a map for the user to use? I want them to be able to filter data in the same way that I can do it from my UI. Is the first map on the following page just a "photoshopped" image? Or is there a way to let the user filter the data on their end? http://cartodb.com/analyze


title: Is is possible to run the QGIS python module server side?

text: I'm pretty new to programming but I want to use the QGIS python module on a server. Basically I am trying to create a JPEG map programatically on a server. I already have a python script to create routes server side I am just trying to find the best way to create a jpeg map with an OSM (or Google) background map, labels, and a legend. QGIS seems like it might be able to do this I am just not sure how difficult or if it's even possible to do on the server using python. I know I could probably do it via Leaflet but since I don't actually need a web map I was hoping QGIS might be a solution. Also it's a Unix server if that matters.


title: How to parse KML data using geotools?

text: Geotools contains some packages dealing with KML: org.geotools.kml and org.geotools.kml.bindings. I assume it was designed to parse KML data. Do you know how to use it?


title: Differentiating a building from a road on a map?

text: This question was originally asked here by Steven and I felt the question was better asked here. Steven would also like to know the sources Google buys it maps from. I remember seeing a post about it on this forum, but couldn't find it when I searched. Please do help...


title: Future for geospatial on App Engine?

text: 2 years ago, this question was posted here at Gillies's blog, on date of 2009-05-15. Future for geospatial on App Engine? At that time there was no news relating to any apps of GIS web services using Datastore (Google's App Engine). Is there any new information on GIS and App Engine over last 2 years? In particular, I am looking for a WFS that uses Datastore on the backend and a browser and a javascript API on the client that will allow people around the world to collaboratively make maps, which can be represented in structured way such as with geoJSON or a shapefiles. In general, what cool map making and spatial querying web services use Google's App Engine Datastore instead of typical relational database ?


title: Arcpy Spatial Join Running Extremely Slow

text: I have a simple arcpy script that runs a series of spatial joins. However, the script is running extremely slow: one set of joins takes over 12 hours. There are a few for loops, but I have timed these and they are not particularly slow. I'm pretty sure the problem is the spatial join function itself. The data sets are fairly large (> 100k observations) but the joins are much faster if run in ArcGIS itself. However, this is not an option as I need the process to be callable from within other code. Has anyone experienced this problem before? I've seen a few posts related to extremely slow processing time with arcpy but mostly related to problems with the structure of loops or compile time (both of which are okay here). Are there any tricks to speeding up arcpy? Code is very simple: i.e." set up join => set up fieldmappings => run join => repeat." I attached the portion of the code that is running extremely slow below.               ### import necessary modules     import arcpy     import os     import csv     import rename          ### set paths     workGDB = r".\spatialJoinDatabase.gdb"     exportGDB = r".\joinExport.gdb"     temporaryGDB = r".\temp.gdb"     exportDirectory = r".\exportFiles"          ### allow arcpy to overwrite existing files     arcpy.env.overwriteOutput = True          ... Other Spatial Joins Here => produces "temp3"          ### Join with waterways and lakes          print("Locals Set")     targetFeatures = os.path.join(temporaryGDB, "temp3")     joinFeatures = os.path.join(workGDB, "lakes")     outputFeatures = os.path.join(temporaryGDB, "temp4")     print("Begin Waterway Operations")          ### only want to keep the "distnaceNAME" and "Name" fields          keep = arcpy.FieldMappings()     keep.addTable(targetFeatures)     newMap = arcpy.FieldMap()     newMap.addInputField(joinFeatures, "NAME")     newColumn = newMap.outputField     newColumn.name = "lakeName"     newColumn.aliasName = "lakeName"     newMap.outputField = newColumn     keep.addFieldMap(newMap)          arcpy.SpatialJoin_analysis(targetFeatures, joinFeatures, outputFeatures,                                "JOIN_ONE_TO_ONE", "KEEP_ALL", fieldNames, "CLOSEST", "", "distanceLake")          targetFeatures = os.path.join(temporaryGDB, "temp4")     joinFeatures = os.path.join(workGDB, "water")          keep = arcpy.FieldMappings()     keep.addTable(targetFeatures)     newMap = arcpy.FieldMap()     newMap.addInputField(joinFeatures, "NAME")     newColumn = newMap.outputField     newColumn.name = "majorWaterName"     newColumn.aliasName = "majorWaterName"     newMap.outputField = newColumn     keep.addFieldMap(newMap)          arcpy.SpatialJoin_analysis(targetFeatures, joinFeatures, outputPath,                                "JOIN_ONE_TO_ONE", "KEEP_ALL", keep, "CLOSEST", "", "distanceMajorWater")          # now fix features on water          inputData = os.path.join(temporaryGDB, "temp5")          # set condition to identify points on the water     condition1 = '("distanceLake" = 0 OR "distanceMajorWater" = 0) AND'     # set condition to identify points that satisfy condition 1 and are close to an Urban Area     condition2 = ' ("distanceUrban" BETWEEN 1 AND 500 OR "distanceCity" BETWEEN 1 AND 500)'     # update values so that distance to nearest Urban Area is zero => name of Urban Area should be correct          query = condition1 + condition2     print(query)          tableRow = arcpy.UpdateCursor(inputData, query)     print(len(tableRow))          for row in tableRow:         row.distanceUrban = 0         row.distanceCity = 0         tableRow.updateRow(row)          del tableRow          print("Operation Complete")


title: Projected maps don't overlap

text: I am very new to arcGIS and so apologise if this is a really simple question, but it has had me (and what I can find on the internet beat for 2 days now. I am trying to build a series of layers of data in arcGIS 9.3.1 for a river catchment in Leicestershire, using data from several sources. All data covers at least the area of interest. I have maps using British_National_Grid and NAD_1983_Albers as the projected co-ordinate system. When projected together one British National Grid layer doesn't all appear with the others. The Alders layers match the majority of the National Grid layers. As far as I can tell the National Grid layers data (see below) is the same for all. > Data Type: Shapefile Feature Class Shapefile: C:\luke-work\Analysis 2 > data\aquatic landscapes\AquaticLandscapes.shp Geometry Type: Polygon > > Projected Coordinate System: British_National_Grid Projection: > Transverse_Mercator False_Easting: 400000.00000000 False_Northing: > -100000.00000000 Central_Meridian: -2.00000000 Scale_Factor: 0.99960127 > Latitude_Of_Origin: 49.00000000 Linear Unit: Meter > > Geographic Coordinate System: GCS_OSGB_1936 Datum: D_OSGB_1936 Prime > Meridian: Greenwich Angular Unit: Degree The other two projections data are > Data Type: Personal Geodatabase Feature Class Location: C:\luke- > work\Analysis 2 data\catchment\catchment.mdb Feature Dataset: Layers Feature > Class: DrainageLine Feature Type: Simple Geometry Type: Line > > Projected Coordinate System: NAD_1983_Albers Projection: Albers > False_Easting: 1500000.00000000 False_Northing: 6000000.00000000 > Central_Meridian: -100.00000000 Standard_Parallel_1: 27.50000000 > Standard_Parallel_2: 35.00000000 Latitude_Of_Origin: 18.00000000 Linear > Unit: Meter > > Geographic Coordinate System: GCS_North_American_1983 Datum: > D_North_American_1983 Prime Meridian: Greenwich Angular Unit: Degree Any suggestions on what I am doing wrong would be greatly appreciated. Luke


title: shape area is not updated

text: I made some holes (difference) in a shape file with Q-GIS 2.2.0, than I imported the shape file in ArcGIS 10.0 and I exported the shape file to a feature class in a File Geodatabase. In the feature class the "shape area" doesn't show the right value. I added a new field and I used "calculate geometry". ArcGis calculated the right surface and it is different from the old one in "shape area". Why ArcGIS 10 is not updating the value? Is there a way to recalculate or edit the standard field "shape area"?


title: How to zoom to selection in arcpy.mapping and only see selected features?

text: I selected data by using select layer by attribute, then I need to only select these selection and zoom in to the selection. my code is this, it can zoom in base on the scale I set, but I see other adjacent polygon that has other Code_No, why?               Layer = arcpy.mapping.ListLayers(mxd, "lyr", df)            arcpy.SelectLayerByAttribute_management(lyr, "NEW_SELECTION", "Code_NO = '%s'"%mycodeNo)                df.extent = lyrLayer1.getSelectedExtent()     df.scale *= 2


title: QGIS view bounding box on Add vector layer GeoGSON protocol?

text: Our server can serve GeoJSON with an optional bounding box HTTP parameter. e.g. `http://our_gis_server.org.au/layername?bbox=xmin,ymin,xmax,ymax` In QGIS I can add this vector layer as a Protocol and point directly at the server. Is it possible to pass the current view bounding box dynamically as part of this query to reduce load on the server? I.E. can QGIS pass in the map bounding box as parameters to the protocol call?


title: Calculate size of area where data values exceed certain limit

text: I am using Python and **Basemap** to plot sea ice concentration within the Arctic region. I need to calculate the total area (square meters) where the sea ice concentration exceeds a certain limit (e.g. 15%), but I am unsure as how to do that. ![Sea ice concentration in the Arctic \(%\)](http://i.stack.imgur.com/R1qbp.jpg) My code for plotting the map is:               levels=np.arange(np.min(mydata),np.max(mydata),5)     mymap = Basemap(projection='npstere',lon_0=0,boundinglat=65)     x, y = mymap(lon,lat)          CS1 = mymap.contourf(x,y,mydata,levels,                        cmap=mpl_util.LevelColormap(levels,cmap=cm.RdBu_r),                        extend='max',                        antialiased=False)          mymap.drawparallels(np.arange(-90.,120.,15.),labels=[1,0,0,0]) # draw parallels     mymap.drawmeridians(np.arange(0.,420.,30.),labels=[0,1,0,1]) # draw meridians     mymap.drawcoastlines()     mymap.drawcountries()     mymap.fillcontinents(color='grey')      My guess is that I need to convert the contours to polygons and then calculate the areas, but I would really appreciate some guidance as to how to do this.


title: Upload / Import to GeoServer

text: Can I upload a shapefile to the GeoServer Import Data page, or does it only import data locally on the server? When I click Browse to select the data to import, it seems to only browse the server filesystem.


title: Openlayers documentation for WPSClient and layer.events.on

text: I am having trouble finding the following documentation:   * The openlayers wpsclient documentation: eg: what are the requirements for the wpsclient javascript to run JTS:overlaps or JTS: Difference              wpsClient.execute({     server: 'local',     process: 'JTS:splitPolygon',   //if this changes to another function     inputs: {     polygon: mypolygon,            //what is required here??     line: myline     }        * The documentation for Openlayers layer.events.on({... _sketchcomplete and more but i dont know what they are and where to find the info_ ). What options are their? I cant find this in the openlayers api documentation. Im sure im missing something obvious here!! so i appoligise if thats the case.


title: Data interoperability json refresh

text: I am using ArcGIS 10.2.2 with the data interoperability extension. I am also building a web app using python, django and django-rest-framework. I'm building a json api that i can add to arcmap using data interop. but i can't figure out how to get arcmap/catalog to dump the cached json file it downloads and replace it with an updated version (as i change the api output). Any ideas?


title: Data Entry Validation and Google Sync

text: I'm relatively new to CartoDB, but am excited about its features. I am trying to do something relatively basic but am having trouble. I want to build a database of points representing locations. Each location has attributes I want to store. I would do it directly through CartoDB's user interface (as it is nice to place the point on the map directly where I want it as opposed to relying on geocoding), but when editing the fields (in map view or table view) there is no input validation. **Q1: Is there a way to create validation in CartoDB?** To get around this, I tried using a Google Form to structure, simplify, and validate the data entry (using multiple choice, etc.). I even followed Andrew's Color World example to append data from each form submission to my CartoDB table. I was able to get it working, but I don't know js and was having issues with unanswered questions and text with special characters (eg , /). Submissions with such entries would not append to my CartoDB table. So... my next hope was to use CartoDB's sync feature with my Google Form Responses table. I was able to set it up, but am disappointed with the infrequent syncing (only every hour?), the read-only/one-way sync, and the inability to (permanently) specify data types (other than strings). In addition, this method did not allow me to automate the geocoding before or after migrating/syncing it to CartoDB. CartoDB did not automatically geocode it. And, even if it did, it only updates/syncs the data every hour. I want my data entry person to be able to quickly verify that the address got geocoded correctly. Once I have the database created, I will also need to periodically update it with new project sites/location entries. Q2: Can you provide a solution for what I want to do or point me to an example or tutorial that can help me? **Need data validation Easy, basic data entry system Need geocoding and ability to verify correctness of placement after each data entry (direct placement on a map would be even better)** (I don't know much html, js, or css, but I could possibly adapt some code if I had an example to work with. However, it needs to operate seamlessly and without my data entry person needing to know any code. This post was helpful, but didn't provide enough information or tools.)


title: Performance issue with SQLite + Spatialite + C#

text: I am having a performance issue with SQLite database (.db) which extends to SPATIALITE I am trying to update Geometry of 1,00,000 records in database (.db) which taking around 50 minutes. Too much slow. my code is like below ::                       for (int q = 0; q < list.Count; q++)                  {              ArrayList castarraylist = new ArrayList();              castarraylist = (ArrayList)(list[q]);                   using (var cmd = new SQLiteCommand(con))                       using (var transaction = con.BeginTransaction())                  {                      cmd.Transaction = transaction;                           for (int y = 0; y < castarraylist.Count; y++)                      {                              cmd.CommandText = Convert.ToString(castarraylist[y]);                                 cmd.ExecuteNonQuery();                      }                      transaction.Commit();                      GC.Collect();                  }              }       Here each castarraylist contains 5000 records. which updating into database with transaction. so loop go through 20 times and complete the update all. While I manually check the time it's increasing the time at each iteration for 5000 records. like 1st 5000 records processing time > 1:11 minute 2nd 5000 records processing time > 1:25 minute 3rd 5000 records processing time > 1:32 minute 4th 5000 records processing time > 1:40 minute 5th 5000 records processing time > 1:47 minute 6th 5000 records processing time > 1:52 minute ... ... ... 17th 5000 records processing time > 3:32 minute 18th 5000 records processing time > 3:44 minute 19th 5000 records processing time > 4:02 minute 20th 5000 records processing time> 4:56 minute Why this happening I don't able to understand. My sourcecode written in C# and my laptop configuration is i5 2.6 GHz, 4 GB RAM, 500 GB HD. I made connection like below ::               SQLiteConnection con = new SQLiteConnection("Data Source=" + fullPath + ";Version=3;Count Changes=off;Journal Mode=off;Pooling=true;Cache Size=10000;Page Size=4096;Synchronous=off");       (*fullpath - is my database path)


title: MapWindow - Load basemap image given extent?

text: I'm new to MapWindow and am attempting to display an image given the geographic extents of an area. Tiling would be ideal, but I imagine that retrieving a simple image from a WMS would be easiest. Are there any simple examples available for this? I'm using MapWinGIS and .NET.


title: Labels aren't showing up when I convert from .mxd to .kml

text: I've been trying to convert a map from .mxd to .kml for google earth. The map has municipalities labeled as well as train stations. When I save it and open it in Google Earth, only the train station labels appear. I have tried moving the municipal layer up in the table of contents before saving. I have also tried taking out the rail stations labels and only having the municipal labels. Neither has done the trick. Any help is greatly appreciated.


title: how to fill voids in ASTER stereo-images for DEM generation?

text: I'm trying to use a stereo-pair of ASTER images to generate a DEM (in ERDAS Imagine). However, the mountains are rather cloudy and the images can't be used as they are. But the relevant parts (valley bottoms) are well visible. I had this idea of cutting the cloudy parts from the stereo-images and inserting the SRTM (in ArcGIS). With this I probably could generate a DEM. But I'm not sure how much sense this makes because I would loose the stereo-view where I insert the SRTM (as I would have to insert the exactly same SRTM-DTM in both ASTER-images). So my questions: do you think it makes sense to insert SRTM-data where there are clouds in ASTER? Alternatively I could use a different ASTER-scene to insert where I have those voids in my original ASTER-set. But then probably there is the question of the viewing angle of the sensor, if it was not exactly the same, this idea is also useless..?!


title: Can I customize a Google Maps layer in terms of content?

text: I don't want all features that come with the Google street layer. I just want basic infrastructure like roads, national parks, cities, streets with their names, etc. Everything that is public is OK. I don't want private companies and shops to appear on the map. I am using OpenLayers to display the Google Maps street layer. So maybe it's possible to choose what to display by using the Google Maps API instead of OpenLayers..


title: Looking for python/VB script to insert a new uniqe field

text: Is there any phyton/VB script which will add a field to shapefile which will be corresponding to first three characters from the shapefile name ? I am working in modelbuilder and arcgis 10. Thanks in advance.


title: Raster calculation confusion

text: I need to do some raster calculation, I've been trying to make it work but I am not sure everything is working properly. I need to apply a formula, but I am unsure if I am doing it properly. The formula is needed to find out the luminance: > L(Bx) = (CN(Bx) – 1) * G(Bx) - Q0(BX) Bx is the band, G is 0.64 and Qo is -5, CN is the number of the band ( in this case 3), the name of the file is "B30@1" Should it look like this ? > ( 3B30@1 - 1 ) * 0.64B30@1 - -5B30@1 After I apply this formula, the image become a gray square, is this normal? Have I applied the formula correctly


title: Where can I find a data model for land use?

text: Does anyone know where I can find a data model for land use? I know about data models available on Esri Support site. Are there any other resources? Thank you in advance.


title: Exporting .jpg files from data driven pages in Arcmap

text: I have generated about 90 maps through data driven pages(ArcGIS 10.0). I now want to export them as .jpg files in either landscape or portrait format. I tried python but it is not working, it worked to export as.png though. Any suggestions?


title: ArcGIS for Server 10.1, Registering a database data store fails with Underlying DBMS error

text: When I registered a database data store to ArcGIS Server in ArcCatalog, the following error occurs: > The connection property set was missing a required property or the property > value was unrecognized. Underlying DBMS error [Unknown errorNo extended > error.] I set up the ArcGIS Server and Oracle on two computers as follow 192.168.1.51               ArcGIS for Server 10.1 64bit     ArcCatalog 32bit     Oracle Instant Client 11.1.0.7 32bit     Oracle Instant Client 11.1.0.7 64bit      192.168.1.52               Oracle Database 11g Enterprise Edition 11.1.0.7.0 64bit      I have successfully made the Oracle 11g "Enabled Enterprise Geodatabase" through the ArcCatalog, and it seems nothing wrong. I added the publisher database connection with the following variables               Database Platfrom: Oracle     Instance: 192.168.1.52:1521/orcl     Authentication Type: Database authentication     User name: sde     password: sde     [v] Save user name and password      And I have added the **Oracle Instant Client** directory to the PATH and environment variable, I am sure that the 64-bit DLL is found before the 32-bit one. But I still get the error message. I am wondering if there is anything else I have to set up?


title: How to format raster legend in QGIS?

text: Using QGIS 2.0, I have created a color legend for raster values in the print composer. By default each color swatch appears to the left of its label in the legend but I want the color swatchs to appear to the right of their labels. I'd also like to be able to specify how the items are aligned with the legend's frame (e.g. right/left). I've looked in the item properties for the raster legend in the print composer, but haven't found any controls to effect this behavior.


title: QGis - QgsComposerScaleBar - nearest km

text: I added a QgsComposerScaleBar to a QgsComposition. I wish to have segment rounded to the nearest km but I don't know how to achieve this. Thank you in advance.


title: Good projection for covering most of the world

text: I'm writing an application for building navigation routes thru the city streets, and not quite sure what projection would be best. I tried EPSG:3587, however, the tools I'm using - GeoTools and PostGIS - seem to disagree on how to interpret this one ( GeoTools and PostGIS with ESPG:3857 ) Could someone recommend a similar projection that would cover most of the world, and be appropriate for looking at the data at the street level only? The data I'm working with is OSM for maps and a lot of WGS84 lat,lon coordinates for way points.


title: pgrouting php problem "pg_fetch_result(): Unable to jump to row 0 on PostgreSQL result index 3"

text: I'm trying to build a simple routing web-app using pgrouting, followed some tutorial and the latest workshop. Because I'm using PostgreSQL 9.3(x64), Postgis 2.1.1, and pgrouting 2.0 so I did some modification to the PHP script so it can be used for my version. Here are the php code:               <?php          // Database connection settings     define("PG_DB"  , "db_mapserver");     define("PG_HOST", "localhost");      define("PG_USER", "postgres");     define("PG_PSSWD", "postgres");     define("PG_PORT", "5433");      define("TABLE",   "jalan1");           $counter = $pathlength = 0;          // Retrieve start point     $start = split(' ',$_REQUEST['startpoint']);     $startPoint = array($start[0], $start[1]);          // Retrieve end point     $end = split(' ',$_REQUEST['finalpoint']);     $endPoint = array($end[0], $end[1]);          // Find the nearest edge     $startEdge = findNearestEdge($startPoint);     $endEdge   = findNearestEdge($endPoint);          // FUNCTION findNearestEdge     function findNearestEdge($lonlat) {          // Connect to database     $con = pg_connect("dbname=".PG_DB." host=".PG_HOST." port=".PG_PORT." user=".PG_USER."    password=".PG_PSSWD);          $sql = "SELECT gid as id, source::int, target::int, geom,               st_distance(geom, ST_GeomFromText(                   'POINT(".$lonlat[0]." ".$lonlat[1].")', 4326)) AS st_dist              FROM ".TABLE."               WHERE geom && st_setsrid(                   'BOX3D(".($lonlat[0]-200)."                           ".($lonlat[1]-200).",                           ".($lonlat[0]+200)."                           ".($lonlat[1]+200).")'::box3d, 4326)              ORDER BY st_dist LIMIT 1";          $query = pg_query($con,$sql);            $edge['gid']      = pg_fetch_result($query, 0, 0);       $edge['source']   = pg_fetch_result($query, 0, 1);       $edge['target']   = pg_fetch_result($query, 0, 2);       $edge['geom']     = pg_fetch_result($query, 0, 3);            // Close database connection     pg_close($con);          return $edge;     }          // Select the routing algorithm     switch($_REQUEST['method']) {          case 'SPD' : // Shortest Path Dijkstra             $sql = "SELECT rt.gid, st_AsText(rt.geom) AS st_wkt,                     st_length(rt.geom) AS length, ".TABLE.".id                  FROM ".TABLE.",                      (SELECT gid, geom                          FROM pgr_dijkstra(                             '".TABLE."',                             ".$startEdge['source'].",                             ".$endEdge['target'].",                             3000)                      ) as rt                  WHERE ".TABLE.".gid=rt.gid;";       break;          case 'SPA' : // Shortest Path A*             $sql = "SELECT rt.gid, st_AsText(rt.geom) AS st_wkt,                       st_length(rt.geom) AS length, ".TABLE.".id                    FROM ".TABLE.",                        (SELECT gid as id, geom                            FROM pgr_astar(                               '".TABLE."',                               ".$startEdge['source'].",                               ".$endEdge['target'].",                               3000)                        ) as rt                    WHERE ".TABLE.".gid=rt.gid;";         break;       } // close switch            // Database connection and query       $dbcon = pg_connect("dbname=".PG_DB." host=".PG_HOST." port=".PG_PORT." user=".PG_USER." password=".PG_PSSWD);            $query = pg_query($dbcon,$sql);             // Return route as XML       $xml  = '<?xml version="1.0" encoding="UTF-8" standalone="yes" ?>'."\n";       $xml .= "<route>\n";            // Add edges to XML file       while($edge=pg_fetch_assoc($query)) {              $pathlength += $edge['length'];            $xml .= "\t<edge id='".++$counter."'>\n";       $xml .= "\t\t<id>".$edge['id']."</id>\n";       $xml .= "\t\t<wkt>".$edge['wkt']."</wkt>\n";       $xml .= "\t\t<length>".round(($pathlength/1000),3)."</length>\n";       $xml .= "\t</edge>\n";       }            $xml .= "</route>\n";            // Close database connection       pg_close($dbcon);            // Return routing result       header('Content-type: text/xml',true);       echo $xml;            ?>      I only modified the php code so it matches the postgis 2.x function and my own road network database. When I get into firefox and tried to compute the route, I got this error from Firebug:               pg_fetch_result(): Unable to jump to row 0 on PostgreSQL result index 3      Am I missing something? The route ran normally in QGIS DB Manager and pgAdmin. glad to have a suggestion from you all regards UPDATE(3:22 pm GMT+7): I rebuilt my whole system, installing PostgreSQL 8.4, PostGIS 1.5, and pgrouting 1.x, using the original tutorial it still gave me the same error. I'm using openlayers as my map renderer.


title: Determining spatial cluster of park types in ArcGIS

text: Question 1: I am trying to determine which (if any) of my different park types are spatially clustered across my study area (city boundary). There are 162 total parks, 9 nominal categories (which represent statistically determined clusters or park 'types' based on 33 social, built, and ecological variables). Also, would I run this analysis on the polygons (park boundaries) or points (based on centroids), or some other way? Question 2 (related but different): I am also wondering if I can run the entire analysis (including spatial) in Arc? The cluster groups I have now reflect the results of an SPSS process including 1. a factor analysis on the 33 variables to reduce redundancy then 2. a two-step cluster analysis (Hierarchical to determine clusters then k-means to form groups). Thank you very much for your assistance.


title: Probably base64 encoded WKB - Reading in QGis/ArcGIS?

text: I have a probably base64 encoded WKB string. How would I load that feature represented by the string into QGis or ArcGIS? I found the QGis QuickWKT plugin, but that does not handle my base64 encoded string. I have some Python programming skills. Get your own WKB from http://nibis.lbeg.de/net3/public/ogc.ashx?Service=WMS&PKGID=22&cardo3SessionGuid=C3_2d80a40c-529e-4574-912b-53a3858b173c&layers=L160&query_layers=L160&width=1&request=GetFeatureInfo&height=1&srs=EPSG%3A31467&version=1.1.1&bbox=3585999%2C5743699%2C3586001%2C5743701&info_format=text%2Fplain&y=1&x=1 if you like.


title: How to accurately place geographic object into photo?

text: My task is to create before/after photo visualization of a wind turbine which is planned to build.   A photo is taken from a point of known coordinates, and coordinates of a turbine are also known.   How to accurately place this turbine into the photo? What I already do is to identify some object on the photo and on an aerial photo. Then I measure distance on the aerial photo map to this object and estimate its heiht. Next I measure a height of this object in the photo in pixels and by proportion calculate height of the turbine.   It's very rough calculation because real object heiht estimation is biased. Next problem is with a direction where to place the turbine. I tried some angular calculations, but they were rather misleading and I ended up with drawing line from point to turbine and further, and looking on what background I should place the turbine. If you have any ideas how do this smarter and more accurately, please share.


title: accessing the output of a python script tool dialog ArcGIS 10.1

text: Hi I have a a fairly simple script tool dialog which has been setup to get 6 parameters from the user and return them as the tool output as a list. the aims is that the out put would then be used elsewhere within a python add-in tool. I know i could revert to storing these values as a text file but that's not my preferred option. I can see the result in the results window but just cant get to it. Below is the script to dialog calls on closure and the seventh parameter is set to an output to store the first 6 parameters. I have been using the pythonaddins module to call the tool through GPToolDialog(Toolbox,Toolname) What i really want is to parse the output parameter from the script tool to the python addin. SCRIPT CALLED BY PYTHON TOOL DIALOG               import arcpy     para1 = arcpy.GetParameter(0)     para2 = arcpy.GetParameter(1)     para3 = arcpy.GetParameter(2)     para4 = arcpy.GetParameter(3)     para5 = arcpy.GetParameter(4)     para6 = arcpy.GetParameter(5)     paraList = [para1,para2,para3,para4,para5,para6]     arcpy.SetParameter(6,paraList)


title: Geoserver: gs:Reproject in SLD fails

text: I want to reproject data via sld:               <?xml version="1.0" encoding="ISO-8859-1"?>     <StyledLayerDescriptor version="1.0.0"      xsi:schemaLocation="http://www.opengis.net/sld StyledLayerDescriptor.xsd"      xmlns="http://www.opengis.net/sld"      xmlns:ogc="http://www.opengis.net/ogc"      xmlns:xlink="http://www.w3.org/1999/xlink"      xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">       <NamedLayer>         <Name>RAIN_SLD</Name>         <UserStyle>           <Title>RAIN</Title>           <FeatureTypeStyle>             <Transformation>               <ogc:Function name="gs:BufferFeatureCollection">                 <ogc:Function name="parameter">                   <ogc:Literal>features</ogc:Literal>                   <ogc:Function name="gs:Reproject">                     <ogc:Function name="parameter">                       <ogc:Literal>features</ogc:Literal>                     </ogc:Function> <!-- end parameter -->                     <ogc:Function name="parameter">                       <ogc:Literal>forcedCRS</ogc:Literal>                       <ogc:Literal>EPSG:4326</ogc:Literal>                     </ogc:Function> <!-- end parameter -->                     <ogc:Function name="parameter">                       <ogc:Literal>targetCRS</ogc:Literal>                       <ogc:Literal>EPSG:900913</ogc:Literal> <!-- end parameter -->                     </ogc:Function> <!-- end targetCRS -->                   </ogc:Function> <!-- end reporject -->                 </ogc:Function> <!-- end parameter -->                 <ogc:Function name="parameter">                   <ogc:Literal>distance</ogc:Literal>                   <ogc:Function name="env">                     <ogc:Literal>BUFFER_DIS</ogc:Literal>                     <ogc:Literal>10000</ogc:Literal>                   </ogc:Function>                 </ogc:Function> <!-- end parameter -->               </ogc:Function>             </Transformation>                  <Rule>               <Name>BUFFER</Name>               <Title>BUFFER</Title>               <ogc:Filter>                 <ogc:PropertyIsEqualTo>                   <ogc:PropertyName>SN</ogc:PropertyName>                   <ogc:Literal>A0001</ogc:Literal>                 </ogc:PropertyIsEqualTo>               </ogc:Filter>                    <PolygonSymbolizer>                 <Stroke>                   <CssParameter name="stroke">#d3d3d3</CssParameter>                   <CssParameter name="stroke-width">2</CssParameter>                 </Stroke>               </PolygonSymbolizer>                    <PointSymbolizer>                 <Graphic>                   <Mark>                     <WellKnownName>circle</WellKnownName>                     <Fill>                       <CssParameter name="fill">#ff0000</CssParameter>                     </Fill>                   </Mark>                   <Size>6</Size>                 </Graphic>               </PointSymbolizer>             </Rule>           </FeatureTypeStyle>         </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>      When i requested the layer with the SLD above, i got an error: > Could not convert the value EPSG:4326 into the expected type interface > org.opengis.referencing.crs.CoordinateReferenceSystem for parameter > forcedCRS. I found an example: Geoserver: Georectification via SLD fails It is useful. But in my situation, 2 CRSs(forcedCRS and targetCRS) must be set. How could I do with it?


title: GRASS raster map types: UInt16 and Float64

text: i am in the processing of using r.out.gdal to export DEMs from GRASS into QGIS this is the **output** of the **first** **try** : using **UInt16** as the data type. r.out.gdal input=DEM0maslTopog@UTM_35S format=GTiff type=UInt16 "output=C:/Users/xxx/Documents/xxx/xxx/GIS/Working/UTM35S/DEM G2Q UInt16/0maslDEM G2Q UInt16 R2" Exporting to GDAL data type: UInt16 Precision loss: Raster map of type DCELL to be exported as UInt16. This can be avoided by using Float64. Checking GDAL data type and nodata value Exporting to GDAL raster r.out.gdal complete. Successfully finished Since the program suggests using **Float64** , this i try all to end up with error 6, which is not detailed in the GRASS manual (both in-Grass and at http://grass.osgeo.org/grass65/manuals/r.out.gdal.html): **Output from Try 2** r.out.gdal input=DEM0maslTopog@UTM_35S format=GTiff type=Float64 "output=C:/Users/xxx/Documents/xxx/xxx/GIS/Working/UTM35S/DEM G2Q UInt16/0maslDEM G2Q Float64 R2" Exporting to GDAL data type: Float64 Checking GDAL data type and nodata value ERROR 6: SetColorTable() only supported for Byte or UInt16 bands in TIFF format. Exporting to GDAL raster r.out.gdal complete. Successfully finished my query: would the error impact significantly on the export (and results on subsequent DEM usage) Bearing in mind that i would like to clip the DEMS to a vector area and use two DEMS (produced and exported in exactly the same way) to calculate volume.


title: GDAL_CALC works but I get a python error at the end of each process that prevents automation

text: Python 3.4.1, GDAL 1.11.0, OSGEO4W, Windows 7, OSGEO4W Windows Shell All 64-bit (ArcGIS is installed as well 10.2 with Python 2.7) So gdal_translate works fine. gdal_calc provides perfect outputs but at the end of each process I get the error below. This is a problem as I cannot automate anything due to this pop up error. I hit close and the error goes away, the output file is 100% perfect but then the next gdal_calc line gives me the same error. This is a sample line of code that gives the error (almost all my gdal_calc commands give the error but not gdal_translate).               gdal_calc.py -A 2000025-2000032.s0481pfv50-sst-16b.hdf.tiff_C.tiff -B C:\temp\QUAL\2000025-2000032.m0481pfv50-qual.hdf.tiffBINARY.tiff  --outfile=A25_32_SSTQUAL4.tiff --calc="A*B"      and this is the error. As I say the process actually works and will continue when I close the error box. After the word DONE appears the error pops up, I close the window and then the next line of code runs, then the error pops up again and so on. ![enter image description here](http://i.stack.imgur.com/aOlah.png) ![enter image description here](http://i.stack.imgur.com/Xdrq0.png)               Additional Error message generated          > Problem signature:       Problem Event Name:   APPCRASH       Application Name: python.EXE       Application Version:  0.0.0.0       Application Timestamp:    5193f3af       Fault Module Name:    ntdll.dll       Fault Module Version: 6.1.7601.18247       Fault Module Timestamp:   521eaf24       Exception Code:   c0000005       Exception Offset: 0000000000053290       OS Version:   6.1.7601.2.1.0.256.48       Locale ID:    1033       Additional Information 1: 8c64       Additional Information 2: 8c64dfac0942d27b36722f7434c64847       Additional Information 3: e156       Additional Information 4: e156d5603c95b65c33088aa70929b4be


title: CartoDB Point Clustering?

text:
I need to create clusters in Cartodb. I reffered to this code https://gist.github.com/andrewxhill/8695353. I am still getting an error.               <!DOCTYPE html>                     <meta charset="utf-8">                  <title>CartoDB Point Clustering</title>         <meta name="description" content="">         <meta name="author" content="">         <meta name="viewport" content="width=device-width, initial-scale=1, maximum-scale=1">         <link rel="stylesheet" href="http://libs.cartocdn.com/cartodb.js/v3/themes/css/cartodb.css" />                           <link rel="stylesheet" href="main.css">

Point clustering

This is a demonstration of point clustering using CartoDB. The method uses an advanced mix of SQL, CartoCSS, and CartoDB.js, continue at your own risk...

Cluster size [48]: <input type="range" name="size" id="cluster_size" value="48" min="20" max="60" step="2">

What's here?

Click a point to find out!


title: How can one get the value of the point where the AxMapControl is centered at?

text: I am using the AxMapControl to display a map and I want to record a point where I am and want to be able to return to the same point later on using the CenterAt() method.


title: ArcGIS/arcpy: Input of in-memory geometry object to Dissolve tool in arcpy script throws Error 000732, input doesn't exist or is not supported

text: In an arcpy script, I run a Union on a list of geometry objects. The output is set to an empty geometry object. The Union runs fine. What I am trying to do is then feed that output into the Dissolve tool. No matter how I format the input, I get the following error: Parameters are not valid. ERROR 000732: Input Features: Dataset in_memory\f8537CACB_AB51_4BDC_B8E7_0895F11E6AF0 does not exist or is not supported ESRI help files indicate that tools set to output to an empty geometry object return a list of geometry objects. The output from the Union tool gives `<Geometry object at 0x26d25910[0x27397068]>` when queried with repr(). I have tried putting the Union result directly into the Dissolve tool and likewise the 0 index item from the result, seeing as it supposed to be a list. Both result in an error, the latter saying that a geometry object is not indexable. Here is a code snippet:               if len(lstHSGeoms) > 1:                 #Union geometries into one geometry                 lstHSUnionGeom = arcpy.Geometry()                 arcpy.Union_analysis(lstHSGeoms, lstHSUnionGeom, "ONLY_FID")                 arcpy.AddMessage('          lstHSUnionGeom:  %s' % repr(lstHSUnionGeom))                 lstHSOneGeom = arcpy.Geometry()                 arcpy.Dissolve_management(lstHSUnionGeom,lstHSOneGeom, None, None, "SINGLE_PART")                 hsOneGeom = lstHSOneGeom[0]      The list lstHSGeoms is a Python list of polygon geometries, which have checked out as valid. The Union runs, but the Dissovle doesn't. I'm stymied. BTW, I tried using an empty string for the fourth parameter of Dissolve (seeing that it is supposed to be a string type), but it didn't make a difference. Any ideas would be greatly appreciated.


title: Unclear on proxy used for AJAX in Javascript

text: I do understand how the proxy helps in transferring the `XMLHttpRequest` between the ArcServer and the Server (IIS/Apache) but where I am confused is that the first time the Html Doc loads, doesn't the Server need to request the ArcServer then as well. How does that request complete without the `proxy`? I mean the first time the page loads it does not need that `proxy` and still the map is loaded in the Map container. But later if we send a `XMLHttpRequest` to the Geoprocessing Tool, the response from the ArcServer can only be received by the Server if we have setup the proxy, right? So how does it happen the first time and not the second time?


title: Clustering Point for StatePlane WKID 2236

text:
Does anyone have an example to point me in the right direction to utilize Clustering with the ESRI JS API working with a State Plane wkid: 2236? I want to provide a solution for an internal application just like the demo on ESRI JS website. Though, all our data and the application I'm supposed to update is defined wkid:2236. I've tried modding the demo code though either I'm overlooking something or it's not possible. If this isn't possible any suggestions? Thanks. * * * Thank you. I had overlooked the "spatialReference" option within the constructor. This gets me to the point of displaying the cluster symbols on the map. Though, when I click the symbol I received the following error and I included the call stack. Any ideas? CHANGES TO ESRI POINT CLUSTERING EXAMPLE
     var initExtent = new esri.geometry.Extent({ "xmin": **432163.89125**, "ymin": **1451074.4432**, "xmax": **710435.747749999**, "ymax": **1626999.3528**, "spatialReference": { "wkid": **2236** } });
     ...
     var basemap = new esri.layers.ArcGISTiledMapServiceLayer("**http://gisserver/ArcGIS/rest/services/BaseMap/MapServer/**");
     ...
     var photos = esri.request({
                         "url": "**js/testdata.json**",
                         "handleAs": "json"
                     });
     ...
      clusterLayer = new extras.ClusterLayer({
                     "data": photoInfo.data,
                     "distance": 100,
                     "id": "clusters",
                     "labelColor": "#fff",
                     "labelOffset": 10,
                     **"spatialReference":"2236"**,
                     "resolution": map.extent.getWidth() / map.width,
                     "singleColor": "#888",
                     "singleTemplate": popupTemplate
                 });

      **ERROR**   Unhandled exception at line 34, column 167548 in http://serverapi.arcgisonline.com/jsapi/arcgis/?v=3.1 0x800a01b6 - Microsoft JScript runtime error: Object doesn't support property or method 'toJson'      **CALL STACK**   > toJson [arcgis] Line 34 Script   >  updateHighlight [PopupBase.js] Line 19 Script   >  onSelectionChange [Popup.js] Line 19 Script   >  select [PopupBase.js] Line 19 Script   >  _updateFeatures [PopupBase.js] Line 19 Script   >  setFeatures [PopupBase.js] Line 19 Script   >  _2b5 [arcgis] Line 15 Script   >  setFeatures [Popup.js] Line 19 Script   >  _addSingles [eval code] Line 275 Script   >  onClick [eval code] Line 126 Script   >  advice [arcgis] Line 15 Script   >  _504 [arcgis] Line 15 Script   >  _onClickHandler [arcgis] Line 34 Script   >  Anonymous Function [arcgis] Line 15 Script   >


title: Add your own data to ArcGIS.com javascript template and host locally?

text: ArcGIS.com has a basic viewer javascript template where you can create a map using ArcGIS.com and then download the files/template to be hosted locally. I won't be able to upload my data to ArcGIS.com so I want to know how hard it is to add my data to the template to be hosted locally. I currently have the template downloaded and hosted locally with an ArcGIS.com basemap. I would like to add my own ArcGIS Server map services to it. Thank you


title: Python Script within ArcMap running GME Command passes a double forward slash when writing the log file

text: I thought I made some serious headway there today. My script works perfectly fine when run from a Python Shell, but in ArcMap, whether the script is run from the Python Command line, or within a toolbox, GME fails/produces an error when loading. Another person is experiencing the identical behavior and has a thread started on the ESRI Forums. Here is a screenshot of his error, my OS seems to cut off more of the error: ![enter image description here](http://i.stack.imgur.com/9B8aI.png) Has anyone else experienced this error? Are you able to execute GME commands from within an ArcMap session? There is only one config file which, while related to logging, I don't believe has any bearing on this. I'm not aware of any other switch when you execute GME from a command line which disables logging.


title: How to measure without coordinates?

text: My current Project is to "map" a turntable in qgis. I know, very unlikely but that's what my teacher wants. My teacher also wants me to measure something in this project (Very specific, isn't it?). And that's where the problem starts. Qgis wants me to define geographical coordinates or something like that. The problem is, that i don't have any coordinates because it's a picture of a turntable and not an actual map. Any ideas? Greeting from Germany.


title: Looping Dataframe in Python doesn't Return Expected Results

text: I'm creating a python addin script to be used on some mxd's I have created as templates. Each map contains 2 dataframes. I need to loop through each dataframe and do some processing on each. I set my loop as follows:               for df in arcpy.mapping.ListDataFrames(mxd, "DFMAIN"):                         for lyr in arcpy.mapping.ListLayers(mxd, "", df):                                # do some processing here     for df in arcpy.mapping.ListDataFrames(mxd, "DFINSET"):                         for lyr in arcpy.mapping.ListLayers(mxd, "", df):                                # do some processing here      In this scenario, if the second dataframe is the 'active' dataframe, it gets chosen first and then some weird stuff happens. I get layers from my first dataframe copied into the second which baffles me and I am not even using any functionality that would copy layers (i.e. arcpy.mapping.AddLayer or anything similar) I just update some definition queries. If I change the template so that the first dataframe is the active one, then everything works fine. I tried indexing the dataframe loops:               for df in arcpy.mapping.ListDataFrames(mxd, "DFMAIN")[0]:      but I get the same weirdness if the second dataframe is active. I guess the easy answer is, make sure the template has the first dataframe as the active one and save it that way, but I have users who may change this unknowingly. I want the geoprocessing to take care of the active dataframe issue. I can't seen anything in the help menu to set the dataframe as the active one.


title: How to determine latitude longitude point from a pixel within a MODIS tile?

text: I need to know latitude / longitude of each pixel in a MODIS tile image in HDF in order to reproject from sinusoidal (MODIS) to mercator (Google Maps). I Can read values from HDF in my C++ program, I can reproject with proj4 having latitude / longitude of a point. I found a MODLAND Tile calculator  with source code in c++ but it have too many dependences in Linux libraries and I'm programming a Windows aplication. Any suggestion will be appreciated


title: how can i install Semi-Automatic Classification plugin manually?

text: I try to install with Fetch Python Plugins. But fetch python plugins show only those which are installed? can i install manually when direct installation failed repeatedly?