title: Just view topology errors in ArcGIS

text: I have a basemap (parcels), which is to be used as target for several conflation tasks. This basemap MUST NOT be modified in any way. The basemap ist part of numerous topology rules. Each time topology is validated ArcGIS inserts vertices in the basemap and/or the other layers. I want ArcGIS to simply show errors, rather then trying to dissolve them. Assigning the highest rank to the basemap doesn't prevent ArcGIS from changing it. Finding the appropriate solution for error correction (e.g. verification of measurements in the field) or accept the error as exception is part of MY job. How can I make ArcGIS to simply show errors without any attempt to fix them?


title: leaflet geojson contextmenu

text: I want to add a context menu on right click for various elements on my geojson layer (I'm doing a road map so on a right click on the road at any part I want to show the context menu). I've managed to get the left click working fine by using the onEachFeature and doing the following               function onEachFeature(feature, layer) {             layer.on({                 click: showAssetInfo,                 contextmenu: contextreg             });              }              function showAssetInfo(e) {             AssetMouseClick(e.target.feature.properties.objectid, e.latlng);         }      For the context menu I have followed the example here . The context menu library is found here I have the following that gets called on the document ready (jquery)               $.contextMenu({         selector: 'path.leaflet-clickable',         zIndex: 99999,         callback: function (key, options) {             var m = "clicked: " + key;             window.console && console.log(m) || alert(m);         },         items: {             "edit": { name: "Edit", icon: "edit" },             "cut": { name: "Cut", icon: "cut" },             "copy": { name: "Copy", icon: "copy" },             "paste": { name: "Paste", icon: "paste" },             "delete": { name: "Delete", icon: "delete" },             "sep1": "---------",             "quit": { name: "Quit", icon: "quit" }         }     });      I've tested it and the selector does return the GeoJson features, also if it attach the same menu to something else it works correctly. Is there something I am missing here? Also is there a good way to pass in the objectid to the menu when it starts up so I can use it when calling the different options of the menu I have created this fiddle to demonstrate http://jsfiddle.net/Q3L4c/22/


title: Reproject raster catalog - ArcGIS 9.3/FME 2012

text: I have a couple of raster catalogs in different file geodatabases and I need to re-project them to a new coordinate system. Now there is no option to project a whole raster catalog, nor can I use the 'Project Raster' tool, as it won't accept a raster catalog as file type (for selecting single rasters within the catalog). There are two options that I thought about so far:   * use tool 'Raster Catalog to Raster Dataset', reproject the raster dataset and bring it back into a catalog. Now the problem is on the one hand that this process takes ages (incl. building of pyramides + calculating statistics) and I loose the single items within the raster catalog (after loading the dataset back into a raster catalog), as it just mosaics the datasets together.   * use Python and loop through raster datasets in raster catalog and reproject them individually. I have not investigated this version, but I doubt whether that will work. In case it would, that could be a solution. Anyone did something like that before? My question therefore is, whether there is a proper/efficient way of reprojecting a raster catalog or the raster datasets within a raster catalog respectively from one coordinate system to another? I basically want to end up with raster catalogs that look exactly as the originals (incl. file names of raster datasets within them), but they do need to have a different projection. Any tips appreciated as my problem is quite urgent! **UPDATE:** I have access to FME 2012 now, never used it before though. Does anyone have an idea whether this exercise could be accomplished using FME? (it looks like FME can only access SDE raster catalogs, so that may not be the full solution either..)


title: viewshed analysis for observers wearing visors?

text: I'm starting with a DEM, and looking to identify every point that can be directly "seen" from any point in a particular polygon on an overlaid layer defining an area on that DEM - like any other viewshed analysis, but with the additional stipulation that the observer points can only "see" at certain downward angles. To help motivate what I'm looking for, my real goal is to identify run-out zones for avalanches. It's easy to identify "slideable slopes": just everything that's, say, 30-45 degrees. But if I want to identify the areas which snow would slide from the slideable areas into, the above quasi-viewshed starts to look like what I'm after. My exploring the surface and hydrology tools (and help files for the tools) in Arc hasn't turned up anything that I think will do the trick. Any ideas - either in Arc or otherwise? Thanks in advance for any insight.


title: add new SRID to sql server

text: Is it possible to add a new SRID to sql server 2012? I've tried               insert into sys.spatial_reference_systems values (4938, 'GEOCCS["GDA94",DATUM["Geocentric Datum of Australia 1994",SPHEROID["GRS 1980",6378137.0,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0.0,0.0,0.0,0.0,0.0,0.0,0.0],AUTHORITY["EPSG","6283"]],PRIMEM["Greenwich",0.0,AUTHORITY["EPSG","8901"]],UNIT["m",1.0],AXIS["Geocentric X",OTHER],AXIS["Geocentric Y",EAST],AXIS["Geocentric Z",NORTH],AUTHORITY["EPSG","4938"]]', 'metre', 1)      But I get > Ad hoc updates to system catalogs are not allowed. Haven't had much luck on google with any answers as well


title: Adjacent DEM files have narrow triangular gaps when joined

text: Newbie here. I've downloaded several DEM files from GeoCommunity which correspond to adjacent regions. When I join them together to form a larger relief map there are triangular 'gaps' on all four sides of each 'tile'. It turns out this is because the elevations are all the same in these triangular regions. Visually, it looks like trimming the triangular regions would allow the remainder to fit properly. Is this because the regions are not in fact rectangular? Something to do with UTM projection? If so how do you know which areas to trim?


title: How to understand google map overlay type from object (polygon or polyLine)

text: var Line = new google.maps.Polyline({         path:predectionArray,         editable: false,         strokeColor: '#FF0000',         strokeOpacity: 1.0,         strokeWeight: 2       });               var Rect=new google.maps.Polygon({         path:predectionArray2,         editable: false,         strokeColor: '#FF0000',         strokeOpacity: 1.0,         strokeWeight: 2       });      the object Line is poly Line and Rect is Polygon .The question is how can i retrive it's type from Object name ,Like Rect.type or Line.type ?


title: Geodesic measurements for short distances throughout US?

text: My ArcGIS experience has been limited to working on large-scale maps with planar projections. I'm moving to a project in which I want to calculate the distance between pro sports facilities throughout the U.S. and their city's respective downtown/central business district (both defined as geocoded points). My primary interest is in relative distances, and my question is whether using a geographic coordinate system for US map points and using "Near" to calculate geodesic distances between stadia and downtowns is a sound approach to avoiding issues that would arise from calculating distance based on a projection of the US?


title: ArcPy Empty Buffers

text: I have written an ArcPy script to process a DTM into a point feature class to be used in Landform classification and ranking. It successfully produces what I need most of the time. Randomly the buffer polygon step produces an empty feature class. If I run the same process "by hand" it works. If I run the script again, or from a value just before it failed, it also works. The issue is not consistent, it may do it on the 3rd loop, or the 300th loop. I've attempted to mitigate this issue by checking my delete variable statements (memory leak management), and by writing a get count followed by an "if" statement that runs a "compact", "repair geometry", "add spatial index", "Delete the last buffer feature class", and "re buffer". But to no avail. The first part of the script is below.. Any ideas would be welcome..               import arcpy     import winsound          ######## Check out any necessary licenses     arcpy.CheckOutExtension("spatial")     from arcpy.sa import *          ######## Script arguments     FolderLocation = "C:\\Data\\Basin\\"     InRaster = "C:\\Data\\DEM\\Basin1"     GDB = "C:\\Data\\Basin1.gdb"          LandForm = GDB+"\\LandForm"     TempP = GDB+"\\TempP"          TopVal = 1037     ProcessingVal = 772     BotVal = 150     Range = TopVal-BotVal     Range2 = ProcessingVal - BotVal     Range3 = Range2     print "There are: "+ str(Range) +" values in this raster. With: " +str(Range2) + " to process"     print " "          #### Naming Conventions     SingleValue = "s_"     BandValue = "b_"     CountValue = "r_"     LandFormIValue = "u_"     PNeigbours = "p_"     Frequency = "f_"     BufferName = "buf_"     MultiName ="m_"     IntName = "i_"     DisName = "d_"          MegaBuffCount = 0     MegaMultiCount = 0     MegaCount = ProcessingVal +1          while (MegaCount) > (BotVal):         MegaCount = (MegaCount-1)         ## Reclass Single Values         SingleValueName = SingleValue + str(MegaCount) +".tif"         SingleValueNamePath = FolderLocation + SingleValueName         SingleValuePolygonName = SingleValue + str(MegaCount)         SingleValuePolygonNamePath = GDB+"\\"+SingleValuePolygonName         SingleValueCountName = CountValue + str(MegaCount)+".tif"         SingleValueCountNamePath = FolderLocation + SingleValueCountName         LandFormIncrementName = LandFormIValue + str(MegaCount) +".tif"         LandFormIncrementNamePath = FolderLocation + LandFormIncrementName              RC1 = str(BotVal)         RC2 = str(MegaCount -1)         NODATA = "NODATA"         RC3 = str(MegaCount)                RC4 = str(TopVal)         Val = 1         Equ = "{0} {1} {2}; {3} {4} {5}".format (RC1,RC2,NODATA,RC3,RC4,Val)         print "Process Number: " + str(Range2) + "   Height Value: "+ str(MegaCount) + " Count Processed: " +str(int(Range3)- int(Range2))         print "A: Raster S Single Band:" + str(Range2) + " Equ: " +str(Equ)         Range = Range -1         Range2 = Range2-1              Raster = arcpy.sa.Reclassify(InRaster,"Value",Equ, "NODATA")         Raster.save(SingleValueNamePath)              ## To Poly         arcpy.RasterToPolygon_conversion(SingleValueNamePath, SingleValuePolygonNamePath, "NO_SIMPLIFY", "Value")         print "B: Conversion Polygon: "+ SingleValuePolygonNamePath              ####### TEST DIAGONALS         print "      B1: Processing Diagonals"         BufferNearName = BufferName + str(MegaCount)         BufferNearNamePath = GDB+"\\" + BufferNearName         MultiPartName = MultiName + str(MegaCount)         MultiPartNamePath = GDB+"\\" + MultiPartName         IntersectName = IntName + str(MegaCount)         IntersectNamePath = GDB+"\\" + IntersectName         DissolveName = DisName + str(MegaCount)         DissolveNamePath = GDB+"\\" + DissolveName              #Spatial Ref: GDA 94 Z50         sr = arcpy.SpatialReference(28350)              arcpy.Buffer_analysis(SingleValuePolygonNamePath,BufferNearNamePath,"100 decimeters","FULL","ROUND","ALL","")         bufcount = arcpy.GetCount_management(BufferNearNamePath)         result_value = bufcount[0]         buffercount = int(result_value)              MegaBuffCount = MegaBuffCount + buffercount              if buffercount == 0:         print "DANGER!!!!!!!!"             arcpy.Delete_management(BufferNearNamePath)             arcpy.Compact_management(GDB)             arcpy.RepairGeometry_management (SingleValuePolygonNamePath)             arcpy.AddSpatialIndex_management(SingleValuePolygonNamePath,0)             arcpy.Buffer_analysis(SingleValuePolygonNamePath,BufferNearNamePath,"100 decimeters","FULL","ROUND","ALL","")             winsound.PlaySound("SystemHand", winsound.SND_ALIAS)                  print "      B2: Buffer Error: " +str(bufcount) +" polygons. This Processing run: " + str(MegaBuffCount)                  bufcount2 = arcpy.GetCount_management(BufferNearNamePath)             result_value2 = bufcount2[0]             buffercount2 = int(result_value2)                  if buffercount2 == 0:                 arcpy.Delete_management(BufferNearNamePath)                 arcpy.Delete_management(SingleValueNamePath)                 arcpy.Delete_management(SingleValuePolygonNamePath)                 winsound.PlaySound("SystemHand", winsound.SND_ALIAS)                 winsound.PlaySound("SystemHand", winsound.SND_ALIAS)                 STOP # Deliberate error              else:             print "      B2: Buffer: " +str(bufcount) +" polygons. This Processing run: " + str(MegaBuffCount)              arcpy.MultipartToSinglepart_management(BufferNearNamePath,MultiPartNamePath)         Multicount = arcpy.GetCount_management(MultiPartNamePath)         Mulresult_value = Multicount[0]         MulCount = int(Mulresult_value)         MegaMultiCount = MegaMultiCount + MulCount              print "      B3: Multi Part: " +str(Multicount) +" polygons. This Processing run: "+ str(MegaMultiCount)


title: How can I do edit-and-continue in ModelBuilder?

text: I'm constructing a model in ModelBuilder, and I ran into an error with one of my tools. I fixed the error and want to continue running my model from that step onward without having to repeat the entire model. Is there any way to do this? Thanks!


title: Does anyone have aquifer depletion/recharge data for China, India, or any other country that isn't the United States?

text: Im doing a project for my capstone class in Geoenvironmental Studies and i am having a hard time finding data for aquifers in other countries.


title: Looking for web service to get State, County, and Place FIPS codes from an address

text: I have been searching for a few days to find a web service that will provide the State, County, and Place FIPS codes from an address input. Most of the services I've found only provide a State, County, and Census block code. For example, for the address 4139 S 143rd Cir, Omaha, NE 68137 would return: State: 31 County: 05 Place: 537000


title: Arcobjects initializing license error

text: I have ArcObjects .NET code that is running on a machine with ArcGIS for Server 10.1. The application worked on the dev server. But when I moved to Test server, I am not able to initialize the licenses. This is the code I am using               RuntimeManager.Bind(ProductCode.Server);     bool succes = false;     AoInitialize arcObjectInitEnv = new AoInitializeClass();     if (arcObjectInitEnv.IsProductCodeAvailable(esriLicenseProductCode.esriLicenseProductCodeStandard) == esriLicenseStatus.esriLicenseAvailable)     {         status = arcObjectInitEnv.Initialize(licenseCode);     }     else     {         throw exception(...);     }      I am sure I authorized the server right. and I am able to run ArcGIS server manager. And test the services and they are working fine. I compile as "Any Cpu" , which is the same I used for Dev. Is there tool that give me "availablity" for server, like desktop? I didn't install ArcObjects SDK on the test server, because I assume it is not required there. Is that true? How can I troubleshoot this?


title: Open Source Engine for Geocoding

text: I would create a service to geocode from existing shapefiles. I know that ArcEngine does not allow the creation of Windows Services. Is there an open source option someone could recommend?


title: Calculating length with Python and the PYTHON_9.3

text: (Please critic wholly if possible) If I have polyline data in NAD83 and I want to calculate geometry using the PYTHON_9.3 parameter from the arcgisscripting.create(9.3), I am getting mixed results. If I calculate geometry in ArcMap, I set the dataframe to the state plane that my data resides, then calculate geometry in the attribute table view but I don't know that I am doing it right. I don't project the data to the state plane that the polyline record resides that I want to calculate geometry on but instead just set the dataframe and let the project on the fly by arcmap work or not work. Now trying to automate the process, I do a gp.selectbylocation against state plane data then clip the polyline data for each selected state plane and then project the data to that state plane and do the "PYTHON_9.3" in gp.calculatefield_management piece. What environment settings should I consider or research more. Thanks guys!


title: ESRI Javascript v API 3.3: Format Info Window Content

text:
I'm attempting to query a parcel layer and display the information in a popup infoWindow when I click on a parcel, and I'm having some issues dealing with different datatypes.               var content = "Address: ${ADDR}" + "
Owner Name: ${OWNER}" + "
Parcel ID: ${APN}" + "
City: ${CITY}" + "
Acres: ${TOTAL_ACRES:NumberFormat}" + "
County Assessor Site"        var popUpTemplate = new esri.InfoTemplate("Parcel", content);      The ADDR, OWNER, APN and CITY datatypes are all string so they appear in the infoWindow without issue, but I'm having trouble getting TOTAL_ACRES and COUNTY_LIN to appear. TOTAL_ACRES is a double, and i thought that using the dojo.number.format method as outlined on the esri website would help, but it doesn't seem to do anything. COUNTY_LIN is a url, and has a data type of text, but it isn't appearing and I'm not exactly sure why. I've included screenschot of the relevant data info from inside ArcGIS, hopefully they are helpful. Any help is appreciated! I'm new at this and very stuck as to what I'm doing wrong. ![APN](http://i.stack.imgur.com/nHEZb.png) ![TOTAL_ACRES](http://i.stack.imgur.com/sEnj4.png) ![COUNTYLINK](http://i.stack.imgur.com/iowaS.png)


title: Using PHP to create a GIS Web application

text: I've developed many GIS Applications using ESRI ArcGIS 9.3 and 10.0 for ASP.NET . I'm wondering if i can do a GIS application using my SHAPE files and Pure PHP scripting . I've googled for that purpose , found many Map Servers solutions , but they were a little bit complicated , my question is : can i have a little Web Application to mange the shape files ,create tiles for them ,deals with Aerial Images, and then display map layers using OpenLayers for example ? i need a lite solution (No php extensions installing required - NO special setup ...etc) . Thank you .


title: How can I Enable Google Earth 3D features in QGIS?

text: I want to ask if there is any plugin or any other way where I can use Google Earth features in QGIS. When I'm saying Google Earth features I mean to be able to see roads, towns, etc and to be able to use the 3D option (Street view tool) I have already installed the OpenLayers Plugin but I did not find something to do it.


title: Copying netCDF coordinate system

text: **UPDATE: 7/26/2013** The issue I'm encountering is a known bug with the netCDF-Java API described in ticket TDS-391: _A user would like to convert radar data to netCDF format. Unfortunately, the data files have structure data, which fails badly for writing netcdf-3 (illegal data type, as netcdf-3 does not support structures) or netcdf-4 (unimplemented) files._ The issue is targeted for a fix in API ver. 4.4, which will hopefully will be out this summer, as per an email from Unidata netCDF Java Support received on July 8: _We are going to implement a solution for this in version 4.4, but it won't be out before our training workshop later this month. We may limit the implementation to NetCDF-4, and offer to skip trying to write structures in NetCDF-3._ * * * I have many level-III NEXRAD precipitation netCDF files which I created from the raw files using the netCDF-Java API (ver. 4.3.17). I am now trying to compute monthly average rainfall rates and returning a new netCDF file which I will then need to convert to ESRI GRID format. I am new to working with netCDF files, but have read up a bit on creating dimensions and variables using the `scipy.io.netcdf` module. I am using the first netCDF file in the directory as a sort of template file from which I'm grabbing the x and y locations and projection information. I have 3 dimensions:               In [4]: ncfile.dimensions     Out[4]: {'textStringSize1': 1, 'x': 131, 'y': 131}      I also have 4 variables:               In [3]: ncfile.variables     Out[3]:      {'FlatEarth': <scipy.io.netcdf.netcdf_variable at 0x67527d0>,      'PrecipArray_0': <scipy.io.netcdf.netcdf_variable at 0x6752690>,      'x': <scipy.io.netcdf.netcdf_variable at 0x67526d0>,      'y': <scipy.io.netcdf.netcdf_variable at 0x6752430>}      I have successfully created my `x` and `y` dimensions and variables and populated the `precip_array_0` variable with my monthly average values. The last piece of the puzzle, I think, is to copy the projection info which I believe is what the `FlatEarth` variable holds. I'm not sure how to look at what's in the variable:               In [5]: ncfile.variables['FlatEarth'][:]     ---------------------------------------------------------------------------     IndexError                                Traceback (most recent call last)     C:\Users\jbellino\Desktop\<ipython-input-5-aec6027d7a77> in <module>()     ----> 1 ncfile.variables['FlatEarth'][:]          C:\Python27_epd32\lib\site-packages\scipy\io\netcdf.pyc in __getitem__(self, index)         837          838     def __getitem__(self, index):     --> 839         return self.data[index]         840          841     def __setitem__(self, index, data):          IndexError: 0-d arrays can't be indexed.               In [6]: ncfile.variables['FlatEarth']     Out[6]: <scipy.io.netcdf.netcdf_variable at 0x67527d0>               In [10]: ncfile.variables['FlatEarth'][()]     Out[10]: ' '      Calling help on the object returns help for the entire `variables` object, not just the `FlatEarth` variable:               In [11]: ncfile.variables['FlatEarth']??     Type:       dict     Base Class: <type 'dict'>     String Form:{'y': <scipy.io.netcdf.netcdf_variable object at 0x06752430>, 'x': <scipy.io.netcdf.netcdf_variab <...> e object at 0x067527D0>, 'PrecipArray_0': <scipy.io.netcdf.netcdf_variable object at 0x06752690>}     Namespace:  Interactive     Length:     4     Docstring:     dict() -> new empty dictionary     dict(mapping) -> new dictionary initialized from a mapping object's         (key, value) pairs     dict(iterable) -> new dictionary initialized as if via:         d = {}         for k, v in iterable:             d[k] = v     dict(**kwargs) -> new dictionary initialized with the name=value pairs         in the keyword argument list.  For example:  dict(one=1, two=2)      I'm assuming that the `FlatEarth` variable uses the `textStringSize1` dimension, but I'm not sure how to create this variable in the new netCDF file or populate it once it's been created. Any help is much appreciated!


title: Using ArcPy for Copy Parallel, Move, etc (from Editor toolbar)?

text: Since VB can activate the Editor, then use the tools listed in the Editor, such as copy parrallel, I am wondering whether Python can do the same job? I have been searching the website all day. Still can't find the answer from early questions-and-answers. Can someone help me? Thanks in advance.


title: Which BC Albers should I use and is there any difference?

text: I was about to calculate the area of some polygons and before I did I did some research. I found that the area of the polygon will be represented in the units of the projection used. So I did some digging to find that the projection I need is BC Albers, but when I went into the list of projections in QGIS, I found that there were two that could the correct ones. The first is EPSG:3005 NAD83 BC ALBERS, and the second one is EPSG:3153 NAD83(CSRS)/BC ALBERS I have projected the polygon as both and saw no difference in the result. I was looking on http://prj2epsg.org/epsg/3005 and it says: "This CRS name may sometimes be used as an alias for NAD83(CSRS) / BC Albers. See CRS code 3153" Just wondering if anyone knows why there are two seemingly identical projections and if there is any difference. Also what does CSRS stand for? Thanks!


title: How to determine optimal windmill placement to maximize production?

text: How can I add automatically the maximum number of points (dots) with 6m diameter and 100m distance each point from other into an area. The dots will present the windmill base (6m) and I want the maximum number of points with 100m distance (Ratio) each in order to reach the highest level of production.


title: clipped data not going into the folder

text: I created a model in python that clips and converts different data sets into DWG based on the the grid number. I will like the DWG data (Output) to go into the folder I created as part of the geoprocessing.               # Import system modules     import sys, string, os, arcgisscripting          # Create the Geoprocessor object     gp = arcgisscripting.create()          # Load required toolboxes...     gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx")     gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")     gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Analysis Tools.tbx")          # Set the Geoprocessing environment...     gp.overwriteOutput = 1     gp.autoCommit = "1000"     gp.XYResolution = ""     gp.XYDomain = ""     gp.scratchWorkspace = ""     gp.MTolerance = ""     gp.randomGenerator = "0 ACM599"     gp.outputCoordinateSystem = ""     gp.ZDomain = ""     gp.snapRaster = ""     gp.configKeyword = ""     gp.outputZFlag = "Same As Input"     gp.qualifiedFieldNames = "true"     gp.extent = "MAXOF"     gp.XYTolerance = ""     gp.MDomain = ""     gp.spatialGrid1 = "0"     gp.outputZValue = ""     gp.outputMFlag = "Same As Input"     gp.geographicTransformations = ""     gp.spatialGrid2 = "0"     gp.ZResolution = ""     gp.spatialGrid3 = "0"     gp.MResolution = ""     gp.ZTolerance = ""     mapNO = "035A2"           # Script arguments...          Expression = "MAPNO = '" + mapNO + "'"      # Local variables...     grid = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\grid"     GRID_200SCALE = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\GRID_200SCALE"     BuildingAnno = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\BuildingAnno"     Building = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\Building"     Building_DWG = "C:\\Export\\Completed\\Building_" + mapNO + ".DWG"     BUILDING__2_ = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\BUILDING"     BUILDING_ANNO = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\BUILDING_ANNO"     Contours__2_ = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\Contours"     ContourAnno = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\ContourAnno"     Contours_DWG = "C:\\Export\\Completed\\Contours_" + mapNO + ".DWG"     CONTOURS = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\CONTOURS"     Rec = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\Rec"     RecAnno = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\RecAnno"     RecAnno__3_ = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\RecAnno"     Rec_DWG = "C:\\Export\\Completed\\Rec_" + mapNO + ".DWG"     CONTOUR_ANNO = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\CONTOUR_ANNO"     Rec__2_ = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\Rec"     CommTowers = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\CommTowers"     COMMUNICATION_TOWERS = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\COMMUNICATION_TOWERS"     Comm_DWG = "C:\\Export\\Completed\\Comm_" + mapNO + ".DWG"     CultAnno = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\CultAnno"     CultAnno__2_ = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\CultAnno"     Cult = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\Cult"     Cultural = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\Cultural"     Cultural_DWG = "C:\\Export\\Completed\\Cultural_" + mapNO + ".DWG"     HydroAnno = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\HydroAnno"     HYDROLOGY_ANNO = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\HYDROLOGY_ANNO"     HydroL = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\HydroL"     HydroP = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\HydroP"     HYDROLOGY_LINE = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\HYDROLOGY_LINE"     HYDROLOGY_POLY = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\HYDROLOGY_POLY"     Hydro_DWG = "C:\\Export\\Completed\\Hydro_" + mapNO + ".DWG"     Rail = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\Rail"     Rail__2_ = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\Rail"     RailAnno = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\RailAnno"     RailAnno__2_ = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\RailAnno"     Rail_DWG = "C:\\Export\\Completed\\Rail_" + mapNO + ".DWG"     Road = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\Road"     ROAD__2_ = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\ROAD"     RoadAnno = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\RoadAnno"     ROAD_ANNO = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\ROAD_ANNO"     Road_DWG = "C:\\Export\\Completed\\Road_" + mapNO + ".DWG"     Trans = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\Trans"     Transmission = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\Transmission"     Trans_DWG = "C:\\Export\\Completed\\Trans_" + mapNO + ".DWG"     Veg = "C:\\Export\\Geodatabases\\Planimetrics.gdb\\Veg"     Veg_DWG = "C:\\Export\\Completed\\Veg_" + mapNO + ".DWG"     Vegetation = "C:\\Export\\BaltimoreAppendedFINALDB.gdb\\PLANIMETRIC\\Vegetation"     MAPNO____001A3_ = "C:\\Export\\Completed\\MAPNO = '" + mapNO +"'"          #### what does this do ?????????????????     Completed = "C:\\Export\\Completed"##########This should contain the expression/concatination that inserts data into its folder                                                                                # Process: Select...     gp.Select_analysis(GRID_200SCALE, grid, Expression)          # Process: Clip (3)...     gp.Clip_analysis(BUILDING__2_, grid, Building, "")          # Process: Clip (2)...     gp.Clip_analysis(BUILDING_ANNO, grid, BuildingAnno, "")          # Process: Create Folder...     gp.CreateFolder_management(Completed, Expression) ######### This should go on top /// Replace this with a create folder function I already created ?????????????          # Process: Export to CAD...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\Building;C:\\Export\\Geodatabases\\Planimetrics.gdb\\BuildingAnno", "DWG_R2000", Building_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")          # Process: Clip (5)...     gp.Clip_analysis(CONTOURS, grid, Contours__2_, "")          # Process: Clip (4)...     gp.Clip_analysis(CONTOUR_ANNO, grid, ContourAnno, "")          # Process: Export to CAD (2)...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\Contours;C:\\Export\\Geodatabases\\Planimetrics.gdb\\ContourAnno", "DWG_R2000", Contours_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")          # Process: Clip (6)...     gp.Clip_analysis(Rec__2_, grid, Rec, "")          # Process: Clip (7)...     gp.Clip_analysis(RecAnno__3_, grid, RecAnno, "")          # Process: Export to CAD (3)...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\Rec;C:\\Export\\Geodatabases\\Planimetrics.gdb\\RecAnno", "DWG_R2000", Rec_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")          # Process: Clip (8)...     gp.Clip_analysis(COMMUNICATION_TOWERS, grid, CommTowers, "")          # Process: Export to CAD (4)...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\CommTowers", "DWG_R2000" , Comm_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")          # Process: Clip (9)...     gp.Clip_analysis(CultAnno__2_, grid, CultAnno, "")          # Process: Clip (10)...     gp.Clip_analysis(Cultural, grid, Cult, "")          # Process: Export to CAD (5)...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\CultAnno;C:\\Export\\Geodatabases\\Planimetrics.gdb\\Cult", "DWG_R2000", Cultural_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")          # Process: Clip (11)...     gp.Clip_analysis(HYDROLOGY_ANNO, grid, HydroAnno, "")          # Process: Clip (12)...     gp.Clip_analysis(HYDROLOGY_LINE, grid, HydroL, "")          # Process: Clip (13)...     gp.Clip_analysis(HYDROLOGY_POLY, grid, HydroP, "")          # Process: Export to CAD (6)...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\HydroAnno;C:\\Export\\Geodatabases\\Planimetrics.gdb\\HydroL;C:\\Export\\Geodatabases\\Planimetrics.gdb\\HydroP", "DWG_R2000", Hydro_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")          # Process: Clip (14)...     gp.Clip_analysis(Rail__2_, grid, Rail, "")          # Process: Clip (15)...     gp.Clip_analysis(RailAnno__2_, grid, RailAnno, "")          # Process: Export to CAD (7)...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\Rail;C:\\Export\\Geodatabases\\Planimetrics.gdb\\RailAnno", "DWG_R2000", Rail_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")          # Process: Clip (16)...     gp.Clip_analysis(ROAD__2_, grid, Road, "")          # Process: Clip (17)...     gp.Clip_analysis(ROAD_ANNO, grid, RoadAnno, "")                                                                                # Process: Export to CAD (8)...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\Road;C:\\Export\\Geodatabases\\Planimetrics.gdb\\RoadAnno", "DWG_R2000", Road_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")          # Process: Clip (18)...     gp.Clip_analysis(Transmission, grid, Trans, "")          # Process: Export to CAD (9)...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\Trans", "DWG_R2000", Trans_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")          # Process: Clip (19)...     gp.Clip_analysis(Vegetation, grid, Veg, "")          # Process: Export to CAD (10)...     gp.ExportCAD_conversion("C:\\Export\\Geodatabases\\Planimetrics.gdb\\Veg", "DWG_R2000", Veg_DWG, "Ignore_Filenames_in_Tables", "Overwrite_Existing_Files", "")


title: What are Differences between Feature and Image Services from ArcGIS Online?

text: Esri has recently released what they call their **Landscape Layers** which look to be mostly image services, rather than map or feature services. These image services look to have feature level detail where you can change the symbology to many different attributes, like a soil type name for example and even show them in a table. I believe they are planning to release feature services for each of or most of these layers. What are the main differences, advantages or disadvantages of these image services vs. feature services? There looks to be a lot of overlap between the two, or am I missing something here?


title: How to create contours from smoothed DEM instead of using Terrain to Raster?

text: I'm fairly new to working with lidar and have been assigned the task of creating 2ft contours for the entire county using USGS-provided lidar tiles. Is there a significant difference in the accuracy of contours when simply generating them from a smoothed DEM versus building a terrain dataset, creating a new raster from the terrain and then generating contours from the derived raster? The contours will become one of our data layers for purchase and may be used for anything from architectural design to wetland mitigation, so I would like them to be as accurate as possible in reducing noise and retaining real information. Any help would be greatly appreciated.


title: In MapBasic, how do I allow a user to choose a folder?

text: I have used the FileOpenDlg() function to allow the user to choose a file. I would like to do the same thing to allow them to choose a folder so I can save files to the folder without the user having to interact every time a file saves. Is there a function for this?


title: Where to find Data with site that have hazardous, toxic or radioactive substances?

text: Does anybody know where to find data that lists the location of sites that have hazardous, toxic or radioactive substances? I'm working on a HUD project and need to be able to find out if sites for new housing is within a certain distance of gas stations or other places that are potentially hazardous. Specifically interested in Florida, USA.


title: Do any of the SQL Server 2008 R2 SP 2 cumulative updates contain the MAXDOP spatial index fix?

text: I am troubleshooting ArcGIS 10.1 performance issues with SQL Server 2008 R2 SP 2 (10.50.4000) and the issue appears to be this MAXDOP issue: https://connect.microsoft.com/SQLServer/feedback/details/625896/spatial-index- not-used-on-multiproc-machine-unless-maxdop-set-downward The issue is that the query optimizer is not using the spatial index when it should. Setting a `OPTION (MAXDOP 1)` query hint on a SQL statement will cause the spatial index to be used as will an index hint, but neither are an option for the Arc-clients. I've read that SQL Server 2008 R2 SP 1 - Cumulative Update 3 and higher has the fix but I'm seeing it at SP2. And I don't see the fix listed in the documentation for any of the CUs, but I'm hoping it is undocumented and present in one of the cumulative updates. Has anyone had this issue at SP2 and fixed it with one of the SP2 cumulative updates?


title: Editing of vector layer programmatically with QGIS

text: I'm trying to edit vector layer using Python console and have couple of questions. Why this approach doesn't work and I have to use `canvas.refresh`? If I will use this code not in python console but in script file should I also use `canvas.refresh`?               # This approach doesn't work     layer = qgis.utils.iface.activeLayer()     layer.startEditing()     layer.moveVertex(7402946,6186167,0,0)          # This one works fine     canvas = qgis.utils.iface.mapCanvas()     layer = qgis.utils.iface.activeLayer()     layer.startEditing()     canvas.refresh()     layer.moveVertex(7402946,6186167,0,0)      What is the right way to perform such operation? I've found that I can make the same things using `QgsVectorLayerEditUtils`:               layer.startEditing()     eu = QgsVectorLayerEditUtils(layer)     canvas.refresh()     eu.moveVertex(7402946,6186167,0,0)


title: Overlapping Polygons in same feature class using Arcobjects

text: I have a polygon feature classs with millions of records, I need to find list of overlapping or intersecting polygons in the same feature class, and how to do it efficiently?   1. I need a ESRI Arcobjects code snippet for the same.    2. I want to know how much percentage overlap is occured.   3. Object ID's of overlapping polygons   4. Attached image and requirement in pictorial view. ![enter image description here](http://i.stack.imgur.com/6LiFb.jpg)


title: Making development patterns more efficient with ZSH

text: A common pattern I find myself encountering is converting a CSV to a projected GeoJSON as described here. That entails 3 steps, and I'd like to get it down to the following:               GEOIFY {file} {originalProjection} {newProjection} =         csv2geojon {file.csv} > {file.geojson} &&         ogr2ogr              -F geojson              -s_srs epsg:{originalProjection}              -t_srs epsg:{newProjection}             {reprojectedFile.geojson} {file.geojson}      If anyone has any good resources on how to create functions and use variables with ZSH (which is what I think I need to do this) that would be beautiful. I must be using the wrong vocabulary, because my googling hasn't yielded anything so far. _As always, thanks!_


title: Geographic Transformations using ArcPy vs. Coordinate Systems used by ArcGis online basemaps

text: I have 3 shape files with undefined projection, but I know which one it was used (Lisboa_Hayford_Gauss_IPCC). I need to 1st Define Projection. It works. I need to transform it from WGS_1984, so I've used:               import sys, arcpy     #globals set env.workspace and other used paths     import globals as g               def main():         # set local variables         inDataset = ["rede_Coimbra_node.shp", "rede_Coimbra_link.shp", "rede_Coimbra_zone.shp"]              prjFile = g.coveragePath + "Lisboa Hayford Gauss IPCC.prj"         prjFileDatum73 = g.coveragePath + "Datum 73 Hayford Gauss IPCC.prj"         prjFileWGS = g.coveragePath + "WGS 1984.prj"         # Geographic transformation -          transformation1 = "Datum_Lisboa_Hayford_To_Datum_73_1"         transformation2 = "Datum_73_To_WGS_1984_4"              sr1 = arcpy.SpatialReference(prjFile)         print sr1.name         #"spatialRef.name = Lisboa_Hayford_Gauss_IPCC"         #sr2= arcpy.SpatialReference(prjFileDatum73)         #print sr2.name         #print (arcpy.ListTransformations(sr1, sr2))         #spatialRef.name = "Datum_73_Hayford_Gauss_IPCC"         sr3 = arcpy.SpatialReference(prjFileWGS)         print (arcpy.ListTransformations(sr1, sr3))         #spatialRef.name = "GCS_WGS_1984"              try:             print "..."             for fc in inDataset:                 fcFullPath = g.coveragePath+fc                 #arcpy.DefineProjection_management(fcFullPath, prjFile)                 outfc= fc.strip(".shp")                 outfc = arcpy.env.workspace+"/"+outfc+"_wgs1984"                 print outfc                      arcpy.Project_management (fcFullPath, outfc,prjFileWGS, u'Datum_Lisboa_Hayford_To_Datum_73_1 + Datum_73_To_WGS_1984_4',prjFile)                           #arcpy.Project_management (fcFullPath, outfc,prjFileDatum73, transformation1,prjFile)                 #arcpy.Project_management (outfc, outfc+"_wgs84",prjFileWGS, transformation2,prjFileDatum73)              except Exception, e:             # If an error occurred, print line number and error message             tb = sys.exc_info()[2]             print "-->Line %i" % tb.tb_lineno             print e.message          if __name__ == '__main__':         main()      I also have a layer from OSM, which I also projected to WGS_1994 CS and it feets perfectly to the World Imagery provided by ArcGIS online or OSM basemap. However, my projected 3 shapefiles (now converted to feature classes) don't seem to fit so well: they have a deviaton from 1-2m to the raster basemap. Did I did something wrong or is it because the online basemaps are in other CS and I can't change that?


title: Dividing map into features (whole squares)

text: I'm using QGIS 1.8.0. I divided a map of a state with its municipalities into many features using a GRID (small squares), but instead of creating one ID per feature (one ID per square), it creates one ID per feature considering also the features divided by the municipalities lines, which I don´t want. I mean if a square is divided by the borderline between two municipalities, it considers two features. I would like to have the map of the state with its municipalities divided into many "whole" squares, and each and every one of those squares should have one unique ID. Could anyone help me? Thanks in advance!


title: QGIS on Windows Surface

text: I was wondering if anyone has had any success running QGIS on the Windows Surface tablet? Specifically with the Windows RT OS. Thanks...


title: Saving geotiff in R - why resolution on x and y is not the same?

text: I have some data and i want to save it as a geotiff in R. This is not a problem but even if in my command line i set the resolution on x as 1 and on y as 1 after the geotiff is generated the resolution on x is 0.999 .... something instead on 1. Resolution on y is 1 as i set it previously. I saw same behavior in ArcGIS as well when i generate geotiffs form data. I am wondering why this happens. Any explanation will be welcome. PS. this time my raster in R has no projections - i have user define "coodinates" which can be considered metric. Thanks, Monica OK here it is some generated data and the R code to generate the geotiff:                      set.seed(333)         v1 <- rnorm(100, mean=0.23, sd=0.47)         v2 <- rnorm(100, mean=0.17, sd=0.83)         v3 <- rnorm(100, mean=0.34, sd=0.68)         v4 <- rnorm(100, mean=0.59, sd=1.25)         v5 <- rnorm(100, mean=-0.27, sd=0.91)         v6 <- rnorm(100, mean=0.03, sd=1.37)         test.m <- cbind(v1, v2, v3, v4, v5, v6)         test.m <- t(test.m)         r2a <- raster(test.m, xmn=1, xmx=dim(test.m)[2], ymn = 1, ymx=7, crs=NA)         r2a              class       : RasterLayer          dimensions  : 6, 100, 600  (nrow, ncol, ncell)         resolution  : 0.99, 1  (x, y)         extent      : 1, 100, 1, 7  (xmin, xmax, ymin, ymax)         coord. ref. : NA          data source : in memory         names       : layer          values      : -3.397683, 4.136409  (min, max)      So, i expected to have resolution 1 x 1 and not 0.99 x 1. It is not really a big deal in the large scheme of things for what i am doing but i have to recognize - i am very curious why. Thanks again, Monica


title: Adjusting GPS collected data to a polyline shapefile

text: I am collecting data from a GPS that is attached to a car. When I display this data in a map, the set of collected points are supposed to match with the street shapefile. However what I see in my map is points overlying places where the car cant go, Ex. Buildings, Parcels etc. I understand this is a GPS accuracy issue. My question is: Is there any script or any way to spatially adjust this points to place them on top of the line so the collected data display on its right position in the map. Thanks


title: Is it possible to redesign OpenStreetMaps tiles?

text: We have an application working with open layers and are using open street maps for the map tiles However they don't look great Is there anyway to style these maps? I read something about Mapbox but it uses Leaflet API, however we use OpenLayers and our app is heavily using it what options do i have? thanks


title: Whats the name of vertices order property?

text: I'm trying to find the next vertices on polyline using python script in order to calculate the angle between them. Whats the name of vertices order property ? Vertices have an attribute for X,Y and for his number- at polyline the number is also the sequence of the vertices. how can i use this property ?   Thank Geog


title: Plotting X, Y coordinates generated from MySQL database in QGIS

text: Currently running QGIS 1.8.0 on Mac OS X 10.8.4. I can plot points on a map using the "Add delimited text layer" plugin to get x-y coordinates from a text file. However, I'd like to be able to extract the x-y coordinates from a MySQL database directly rather than having to create an intermediate text file. Can this be done? The database consists of numerous tables containing info about postcodes, eastings, northings, and other information. I can select "Add vector layer" and connect to a MySQL database but I'm not given the option to provide the SQL query necessary to extract the coordinate data or to define which variables should be defined as x and y coordinates.


title: Python Delete point in a distance of x meters

text: I took this code from this source Select maximum number of points more than x meters apart It seems to work for other users, but for me it's an endless loop. I do not understand why, though the code makes sense. The code itself gives me the result I want to get. Can you help me out of this endless loop: #               import arcpy, sys          feature = arcpy.GetParameterAsText(0)     distance = 500          arcpy.gp.overwriteOutput = True          #running NEAR analysis - every point gets attribute of a distance to the nearest point      #in same feature class     arcpy.AddMessage("running first near analysis")     arcpy.Near_analysis(feature, feature)          arcpy.AddMessage("inserting cursor")     cur = arcpy.UpdateCursor(feature)     row = cur.next()          arcpy.AddMessage("starting loop")          i=0     while row:         i+=1         #fids list will store list of deleted points so if any other point will have          #deleted one as the nearest and distance < 150 will not get deleted as this         #distance is no longer true         fids = []         while row:             if row.NEAR_DIST < distance:                 try:                     #it seems I didn't know if .. in .. at the time ;) such a fun to dig                     #this script up! index throws an exception if element is not in the                      #list                     fids.index(row.NEAR_FID)                     arcpy.AddMessage("OBJECTID = " + str(row.OBJECTID) + " is listed!")                 except:                     arcpy.AddMessage("deleting OBJECTID = " + str(row.OBJECTID))                     fids.append(row.FID)                     cur.deleteRow(row)                     d = 1             row = cur.next()         del cur, row, fids         try:             #this idiotic test is to break the loop when no points will have              #NEAR_DIST < 150, shameful - I know!             if d == 1:                 pass         except:             sys.exit()         d = 0         arcpy.AddMessage("loop iteration " + str(i))              #and again we go..         arcpy.Near_analysis(feature, feature)         cur = arcpy.UpdateCursor(feature)         row = cur.next()      The fonctionnal CODE:               import arcpy, sys          feature = arcpy.GetParameterAsText(0)                    def nearRoutine():         #calculate the distances using the current dataset         arcpy.Near_analysis(feature, feature)              #iterate through any features which are within the distance         cur = arcpy.UpdateCursor(feature, '"NEAR_DIST" < 500')         row1 = cur.next()         while row1:                  #this point is within the distance of its neighbor, so delete it             cur.deleteRow(row1)                  #now re-run this routine on the new dataset             del row1, cur             cur = arcpy.UpdateCursor(feature, '"NEAR_DIST" < 500')             row1 = cur.next()             nearRoutine          #call the recursive routine. It will get progressively faster to run as it will loop through fewer points each time     nearRoutine()      Thanks @Stephen Lead


title: How to prepare or normalize raw shape (.shp) files for Openlayers use

text: Hi all it seems like I am having difficulties in normalizing raw map files' resolutions to openlayer's use most probably in parameters mapext and imgext. I hope you can help me with this. Using raw (.shp, .shx, .xml, .prj, .dbf) and associated raster files for landuse-related data included in an ArcGIS 10 project (.mxd) that I've obtained from a credible private source, I have used Mxd2Qgs to read .mxd project to .qgs. Then I used quantum GIS' Mapserver Export plugin to export selected layers to MAP FILE. I can view the generated map on browser but cannot do so when I use it in openlayers. The generated map file has an extent: EXTENT 272813.696194 1546212.773835 340292.781152 1587773.606956 (its quite huge). I have tried adding to openlayer's map initialization properties maxResolution, numZoomLevels, units, projection and displayProjection but there is no map rendered, although the map has been queried. Btw I have just guessed the maxResolution here, copying from other huge datasets online :) The map has never been queried when I added the map file's EXTENT to maxExtent: OpenLayers.Bounds(...). The same thing goes when I used the files included in the MS4W tutorial #1 data set. My questions would be:   * How to know a given map's maxResolution parameter for openlayers? Can this be guessed?   * Are there any special processing needed for map/shapefiles before they can be used in openlayers?   * How should a map file be written so that it could be used in openlayers? Thank you so much in advance for your tips. EDIT: Attached MAP FILE and javascript files.


title: How to convert coordinates from meters to decimal degrees?

text: I have a shape file. The coordinates (Lat/Long) are projected as meter (555000, 155000). But, I want to re-project it to decimal degrees(51.523666, -0.07467). How can I do it with ArcMap or Quantum GIS? Thanks, In advance


title: How to reduce trajectories with Bellman algorithm?

text: I need an algorithm to reduce trajectories with Bellman algorithm. I haven't found anything searching in Google. Can someone help me? It would be cool have it in MATLAB. Thank you


title: Query to find lat/lon of a house by street.name & house.number

text: How can I get house coordinates (latitude, longitude) by street.name and house.number? If it is possible, how I can get it? Many thanks!


title: How to debug ogr2ogr ERROR 1: Did not get at least 2 values or invalid number of set of coordinates?

text: I'm trying to convert a government data contour GML (http://geo.data.linz.gv.at/katalog/geodata/hoehenschichten/2011/1m-Hoehenschichten.gml) to Shapefile. It seems like it was created by FME:               <?xml version="1.0" encoding="UTF-8"?>     <gml:FeatureCollection xmlns:gml="http://www.opengis.net/gml" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:fme="http://www.safe.com/gml/fme" xsi:schemaLocation="http://www.safe.com/gml/fme 1m-Contour.xsd">     <gml:boundedBy>     <gml:Envelope srsName="AUT-GK31-5" srsDimension="3">     <gml:lowerCorner>67741.3127211605 341935.72972973 245</gml:lowerCorner>     <gml:upperCorner>79959.5666279546 360502.420825918 602</gml:upperCorner>     </gml:Envelope>     </gml:boundedBy>     <gml:featureMember>     <fme:Contour gml:id="ida7f8be19-1a86-42be-886c-2238b778e00d">     <fme:_elevation>434</fme:_elevation>     <fme:_length>2987.20354766336</fme:_length>     <gml:curveProperty>     <gml:LineString srsName="AUT-GK31-5" srsDimension="3">     <gml:posList>68485.5237522595 355869.11751307 434 68487 355867.757575758 434 68488.6860465116 355866 434 68492 355862.369426752 434 68493.6929133858 355861 434 68497 355858.324840764 434 68498.921052631$     </gml:LineString>     </gml:curveProperty>     </fme:Contour>     </gml:featureMember>      Ogrinfo (OGR 1.9) recognizes multi line strings:               > ogrinfo 1m-Hoehenschichten.gml      Had to open data source read-only.     INFO: Open of `1m-Hoehenschichten.gml'       using driver `GML' successful.     1: Contour (Multi Line String)      But                ogr2ogr -f "ESRI Shapefile" -nlt MULTILINESTRING25D  -skipfailures  out 1m-Hoehenschichten.gml Contour      fails with               ERROR 1: Did not get at least 2 values or invalid number of     set of coordinates <gml:posList>68485.5237522595 355869.11751307 434 68487 355867.757575758 434 68488.6860465116 355866 434 68492 355862.369426752 434 68493.6929133858 355861 434 68497 355858.324840764 434 68498.9210526316 355856 434 68502 355851.9375 434 68502.5818965517 355851 434 68503.5340909091 355849.465909091 434 68505.6853448276 355846 434 68507 355843.896551724 434 68509.2459893048 355841 434 68512 355836.846774194 434 ....      My interpretation so far is that OGR does not recognize that the file contains 3D geometries. How can I find out if the file is broken or whether I'm dealing with an OGR bug/missing feature?


title: Recommended Eclipse Plugins for development with OpenLayers and GeoExt?

text: we plan to create an map oriented online application and muse about the 'best' / helpful plugins for Eclipse. The used technology is OpenLayers and GeoExt. For example is the Javascript development tools (JSDT) helpful? Or - are there any other tools, not necessarily Eclipse plugins, which help during development? Can you share your experiences? D3


title: Newbie to GIS Development

text: Greetings everyone, I am new to GIS development and was wondering if there are tutorials available for development in ArcGIS 10. I have some experience writing simple scripts in python and fooling around in model builder, and I want to take the dive and make an effort to become a real programmer. I need some more info to get a better grasp on doing some real GIS programming to make some custom tools in ArcGIS 10. Where/how did you guys learn all your GIS development skills? Any book/web site recommendations? Just an FYI I will being doing my coding in visual studio using the .net framework. Thank you all for your recommendations.


title: Alternatives to Windows Explorer and ArcCatalog when copying a file gdb

text: I need to copy a file gdb from my local machine to a server using a remote desktop connection. I have been doing this using a simple copy/paste operation in Windows Explorer. I know that it is best practice to do this using ArcCatalog, but in this situation I can't. What I would like to know is how risky this procedure is. I know there is a chance that there may still be locks on the datasets, or other processes accessing it, which could cause the data to be copied partially/incorrectly. I have indeed noticed this where just having ArcCatalog open on my local machine (and not inside the folder the containing the gdb) could cause the paste to fail/complete incorrectly. What I have found is that it is quite inconsistent. Yesterday I had an mxd open which was referencing the file gdb. I then copied it through Windows Explorer and pasted it on the server. I was able to work with the contents on the server without a problem. Today, the first thing I did after logging in was to copy and paste the same file gdb onto the server again. When I tried to open it in ArcMap, it was identified as a normal folder `test.gdb` and not an actual geodatabase with the name `test`. Would using something like TeraCopy help in this case, or do I just have to hope that all the data will get carried over each time? I see a similar question noted that ArcCatalog is the safest for moving the gdb (I agree, I would not dare to cut and paste the gdb through Windows).


title: 1'' DHM: Hillshade is too rough - how to "generalize"?

text: I got this hillshade from my DHM (I used the hillshade-module from qgis). For my purpose there is too much relief detail (roughness). How can I get rid of this in QGIS or GRASS? ![hillshade from 1" DHM](http://i.stack.imgur.com/oINrc.png) Thanks! Martin


title: ESRI Flash Command in custom code

text: I would like to implement the same flash feature command that ESRI uses when accessing the context menu on a table row. I have used the following successfully to flash a selected feature:               IActiveView activeView = ArcMap.Document.ActiveView;     IFeatureIdentifyObj featIdentify = new FeatureIdentifyObject();     featIdentify.Feature = feature;     IIdentifyObj identify = featIdentify as IIdentifyObj;     identify.Flash(activeView.ScreenDisplay);      This code works, but does not have the 'crosshairs' additional effect that the ESRI flash command does. I was looking possibly executing a command with ICommand.Execute, but could not find an ID for the Flash Command from the table context menu. Anyone have success with implementing the same Flash Command that ESRI uses within custom code? Does anyone know where to find the ID used for the Flash Command?


title: Local variance image in python using gdal and a running window approach

text: I want a local variance image with a 3x3 of a geospatial raster image using python. My approach so far was to read in the raster band as an array, then using matrix notation to run a moving window and write the array into a new raster image. This approach worked well for a high pass filter as described in this tutorial: http://www.gis.usu.edu/~chrisg/python/2009/lectures/ospy_slides6.pdf Then I tried to calculate the variance with several approaches, the last one using numpy (as np), but I just get a gray image with the same value everywhere. I am open to any kind of solution. If it gives me the average local variance in the end, that would be even better Thanks, Mace                 rows = srcDS.RasterYSize       #read in as array       data = srcBand.ReadAsArray(0,0, cols, rows).astype(np.int)            #calculate the variance for a 3x3 window       outVariance = np.zeros((rows, cols), np.float)       outVariance[1:rows-1,1:cols-1] = np.var([(data[0:rows-2,0:cols-2]),         (data[0:rows-2,1:cols-1]),         (data[0:rows-2,2:cols]  ),         (data[1:rows-1,0:cols-2]),         (data[1:rows-1,1:cols-1]),         (data[1:rows-1,2:cols]  ),         (data[2:rows,0:cols-2]  ),         (data[2:rows,1:cols-1]  ),         (data[2:rows,2:cols]    )])       #output        outDS = driver.Create(outFN, cols, rows, 1, GDT_Float32)       outDS.SetGeoTransform(srcDS.GetGeoTransform())       outDS.SetProjection(srcDS.GetProjection())       outBand = outDS.GetRasterBand(1)       outBand.WriteArray(outVariance,0,0)       ...


title: How to change the color of the layer after searching

text: My task is searching a single house from a layer and then zoom to the center. I did this part but i can't clearly identify this house. I am try to do this by adding a fillcolor to the searched area. but i can't.               var point1 = new OpenLayers.Geometry.Point(mins.x, mins.y);                 point1.transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));                  point1.style={fillColor: '#FF0000'};// coloring the point     var point2 = new OpenLayers.Geometry.Point(maxs.x, maxs.y);                 point2.transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));                  point2.style={fillColor: '#FF0000'};// coloring the point     var bounds = new OpenLayers.Bounds();                 bounds.extend(point1);                 bounds.extend(point2);                 bounds.toBBOX();                      map.zoomToExtent(bounds);       but it can't change the previous color eg: see this page https://gis.fatih.bel.tr/webgis/default.aspx


title: How to customize logic for defining default speed limits for roads in osm2po

text: I need to customize osm2po logic in defining possible speed limits of ways. osm2po takes into account only class of road while I need to consider one more parameter: is road inside or outside some city, i.e. for example: speed limit for higway=trunk inside city = 60km/h while outside it should be equal to 90km/h. What is a proper way to implement such logic? Is it feasible at all to do withing osm2po? Or may be it better/easier to load graph data into Postgres/PostGIS using osm2po generated sql-script and do speed/cost corrections there (assuming that we already have all city borders loaded into Postgres/PostGIS db) and use pgRouting as routing engine? Any advices highly appreciated.


title: How to determine the longest side of a polygon in QGIS?

text: Is there a possibility to determine the longest side of a polygon with QGIS? As I need to do that for hundreds of polygons I would prefer an automized method (like a plugin). So far my idea was:   * transforming a polygon into a line   * cut it on the vertices   * measure the length   * here I became stuck, how to become the longest line without checking the whole bunch of lines one by one? looking forward to your answer


title: Add metadata to dataset (feature class, coverage or shapefile) from a .csv using arcpy?

text: I am looking for a script that will take metadata information provided in an excel table, convert it to a .csv and import the values to populate some of the dataset's ArcGIS metadata elements. I came across this post today and it is essentially the reverse of what I need (it grabs metadata elements and exports them to a .csv) and thought it might help: Creating a table containing all filenames (and possibly metadata) in a File Geodatabase I am a novice programmer (with some python under my belt) and am new to this site.


title: How to find the most east or west longitude around the date line

text: I've got a python script that is comparing extents of rasters in various projections. Generally it works, but sometimes, especially when the rasters are global in extent or cross the dateline it fails. In effect I need an algorithm that will return if a longitude is 'more east' of another one. A simple greater than breaks down is some cases. For example:               lon1 = 175     lon2 = 177     lon3 = -181          lon1 > lon2 #returns True lon2 is 'more east' than lon1       lon1 > lon3 #falsely returns False even though lon3 is more east than lon 1        And I don't think I can do a simple transform of my longitude to 0-360 as I don't have control over the incoming coordinate system. For example, in one custom Mollweide projection I'm working on 'x' is in meters and range from -18040094 to 18043982


title: Do Label Expressions Slow Down Drawing Speed?

text: Further to my earlier question about using different fonts on the same PC, I am wondering if label expressions would slow down drawing speed for a MXD? Here is my label expression:               def FindLabel ( [GEONAME] ):           import os           if str(os.getenv("COMPUTERNAME")) == 'GIS-HP':             return ("<FNT name='Adobe Fan Heiti Std B' size='8'>" + [GEONAME] + "</FNT>")           else:             return ("<FNT name='Arial' size='8'>" + [GEONAME] + "</FNT>")      This works like a charm and I haven't noticed anything yet, but I feel like adding this to all of my layers will dramatically slow down the drawing speed on my MXD. Drawing speed is one of the more important things for me as I share my MXD's with other users for them to do basic analysis. Usually I save a separate MXD, open it and change all the labels to something easy before sharing. If this doesn't slow down the drawing speed, that would be fantastic and save me some time down the line. Anyone have any ideas?


title: Chorpleth Data visualization in Mapbox Leaflet API

text: I am using Mapbox Leaflet API. I have followed the following tutorial of Mapbox. https://www.mapbox.com/mapbox.js/ (Please change to Data visualization Tab) Onclick its zooming on state and hover its showing popup. I am trying to get the state name on click to show more data from database in another div. If I can get the state name on click with zooming then my problem can be solved. How can I do this? I used the following code. In which the hover data shows in top right corner. Can I show it in a fixed div?                           var info = L.control();                           info.onAdd = function (map) {                 this._div = L.DomUtil.create('div', 'info');                 this.update();                 return this._div;             };                  info.update = function (props) {                 this._div.innerHTML = (props ?                     '<b>' + props.NAME + '</b>'                     : 'Hover over');             };                  info.addTo(map);


title: Geodatabase add current date and time to attribute

text: Is it possible in Arcmap 10.0 to create a geodatabase feature class that has a date field which auto populates with the current date and time, every time a new feature is created? Is this possible without programming? Have checked ESRI help but nothing. The feature class is not in ArcSDE so it is not versioned.


title: How to style virtual vertices when using the ModifyFeature control with OpenLayers?

text: I am building an Android application that uses OpenLayers embedded in a WebView. Being mindful of the imprecision of touch interfaces I am trying to style all controls to be large and easy to operate. I have followed the documentation on feature styling but have been unable to control the size of the virtual vertices when using the createVertices mode of the ModifyFeature control. By way of example consider the following style with a pointRadius of 5:               var customStyle = new OpenLayers.Style({                         fillColor: "black",                         strokeColor: "red",                         pointRadius: 5,                         strokeWidth: 30                     });          var customStyles = new OpenLayers.StyleMap({                         "default": customStyle,                         "select": customStyle,                         "temporary": customStyle                     });          var polyLayer = new OpenLayers.Layer.Vector("Polygon_Layer", {styleMap: customStyles });      ![small grips](http://i.stack.imgur.com/eymJl.png) Meanwhile the following has a pointRadius of 50:               var customStyle = new OpenLayers.Style({                         fillColor: "black",                         strokeColor: "red",                         pointRadius: 50,                         strokeWidth: 30                     });      ![large grips](http://i.stack.imgur.com/oEH1I.png) As you can see this has no impact on the size of the virtual vertices midway along each edge. Does anyone know a technique to control their size?


title: QGIS Map server on MacOS -- xml errors

text: I have just installed the map server on my mac OS laptop running 10.9 (maverik). It all seemed to install OK but I did note the warnings that the Map server build was not tested on 10.9. when I point the browser to: `http://localhost/qgis- mapserv/qgis_mapserv.fcgi?` I get:               This XML file does not appear to have any style information associated with it. The document tree is shown below.      `<ServiceExceptionReport xmlns="http://www.opengis.net/ogc" version="1.3.0"> <ServiceException code="OperationNotSupported">Please check the value of the REQUEST parameter</ServiceException> </ServiceExceptionReport> ` The site www.opengis.net is home to PURL so I am very confused. I was following the tutorial http://hub.qgis.org/projects/quantum- gis/wiki/QGIS_Server_Tutorial. It was very good a covering installation of the components but I could not figure out the configuration section. It seems to assume you understand how the whole thing works, it does not even give the names of the config files (I did figure that out, I think). The more I look with Google the more confused I become as there appear to be several ways of doing this all with different instructions. Long term I will not be using MacOS to host this but it is my development platform.


title: Tilemill label and Shield are overlapping, how can this be solved?

text: I am trying to use tilemill labels and Shield together for my road designing. But there is a problem with that. Please see the following image for better understanding. The shield is ovarlapping the Road name. ![enter image description here](http://i.stack.imgur.com/OX3Z1.jpg) I used the following code for my .mss               #road::Labels {           text-name: [ROAD_NAME];           text-face-name: "Candara Bold";           text-size: 16;           text-fill: #000000;           text-placement: line;            text-halo-fill: #f6f3f3;           text-halo-radius: 1;           text-spacing: 400;           }         #road::N_Shield {             [ROAD_TYPE =1]{             shield-file:  url("C:\wamp\www\shield.png");             shield-name: [ROAD_NO];             shield-face-name:"Arial Bold";             shield-size:16;             shield-fill: black;             shield-line-spacing:5;             shield-min-distance:1;             shield-allow-overlap:true;             shield-placement:line;             shield-dx:150;             shield-spacing:450;             }          }      Anyone can help me what to do.


title: How to use PostGIS to handle complex geoprocessing workflows?

text: Our organization is considering moving our geoprocessing workflow to PostGIS. We are currently using ArcGIS, with a plethora of custom Python tools used in ModelBuilder. We're moving most of our data into PostGIS to be consumed by a variety of apps, and now we're asking if it also makes sense to perform the data processing there as well. We process data to be compatible with our software. A customer purchases our software, gives us their data, and we process it to be optimized for use in our software. This requires us to build a variety of tools to handle varying qualities of input data. We can't expect to receive data in a particular format or schema, so we build tools to map input fields to output fields, parse single fields into multiple fields, merge multiple datasets, etc. We also perform spatial joins, intersections, trim whitespace and concatenate fields, and many other common operations. PostGIS appears to be perfectly capable of performing all of our processing needs. For those of you who use PostGIS to do your data processing, do you have any advice for organization, tools to use, etc.?   * do you use it in conjunction with QGIS python processing?   * are people using a Python ORM for non-web processing? I've been leaning towards using GeoDjango since it has a Python ORM for PostGIS. Our initial test of using PostGIS to process data has many large SQL text blocks in Python code and we are thinking that the GeoDjango ORM may help with creating more manageable and readable code. There's also the GeoAlchemy ORM that interacts similarly with PostGIS, and doesn't appear to be as web-specific as Django is. I haven't heard of people using PostGIS to do geoprocessing as much as I see people using QGIS or ArcGIS, so I want to know if it is a comparable alternative.


title: Raster Data to Polygon

text: 1. I need to "translate" raster information to a set of polygons.   2. The raster to polygon tool will not work.   3. I created a fishnet with the same resolution of the raster and converted it to polygons.   4. I now have a "Grid" of polygons with the same resolution as the raster.   5. I would like to populate the attribute table with the value of the cells in the raster data In the images below you will see my "Grid" of polygons and my raster with the "grid" overlay-ed. Now how do i get the raster values into the attribute table of the polygons? Any Suggestions?!["Grid" Of Polygons](http://i.stack.imgur.com/YHl8M.png)![Raster with "Grid" of polygons overlayed](http://i.stack.imgur.com/od6MU.png)


title: Delete photos with no location in EXIF information

text: I can extract the locations of my photos in QGIS using Photo2Shape or Geotag and import photos, but I want to be able to delete the photos that do not have location information with them (they are backed up elsewhere). At the moment I have to look in the attribute table which photos get imported with a latitude and longitude value of 0. Then I have to find these files manually and delete them. This is very time consuming and it would be great if it was possible to do this automatically somehow. I guess a another possible workaround would be for the photos that have locations to be copied to a sub folder. I just wondered if anyone knew of a tool (GIS or otherwise) that could achieve this before I set about writing a script to do it for me.


title: QGIS - load data by SQL window

text: I have problem with QGIS (2.0.1), PostGIS (2.0.4) combination. \- A can load tables or views by Add PostGIS layer. This is OK. \- But a can't load data by DB Manager and SQL Window. Executing query run but load result is problem. Error message is: > _PostgreSQL layer has no primary key_ But Primary key exist in table. DB Manager find out them itself. Do you have any idea where is problem?   Could my problem be with different schema then public? Edit: Do this when I click to _Load Now_ only. And realy do this outside public schema. In public schema it's OK.


title: KML to SHP with elevation not equal to Z value

text: Can I convert a Google Earth KML file to a shapefile with the elevation coordinates not as Z values? I can add a vector layer using a KML file and then save it as a shapefile but apparently the elevations as saved as Z values which I cannot have in the application into which I need to import the resulting shapefile. Thanks


title: Linux Client for OpenGTS server

text: I've installed an OpenGTS server and successfully linked it with the GPS2OpenGTS_Pro application for android. The buses I'm tracking offer 3G with a wireless AP, however whenever the bus stops, the AP is shut down to save battery. The android device has trouble reconnecting all the time, so I've opted for a Raspberry Pi that can connect to the AP through ethernet. I've got gpsd to read the NMEA output from a BU-353, but I'm stuck trying to figure out how to stream the GPRMC information to my OpenGTS server. I haven't been able to find an application that uploads GPRMC data, and couldn't find a tutorial for how to write a script to do the same (possible using "cat /dev/ttyUSB0" and filtering for just GRRMC lines?). Is there an easy way to do this in python or PHP? Example HTTP POST: http:// 192.3.203.99:8080/gprmc/Data?acct=ACCOUNTNAME&dev=VEHICLENAME&gprmc=$GPRMC,020114,A,3128.7540,N,14257.6714,W,000.0,000.0,131113,,*7


title: ArcMap, Joining information from lyr with attribute table

text: I'm working on a landcover project with ArcMap 10.1 right now and use the ModelBuilder a lot. So far I extracted a table from a certain layer using the Summary Statistics tool. However, one field of the mentioned table does only contain a code that isn't of much use to me. I do have a lyr file containing the right legend, but I cannot simply join the lyr file with my extracted table. Is there some way to connect these sources of information? I do need to write some documentation on what tools I used during the project (using the ModelBuilder) and therefore don't want to resort to just write the correct code down manually within my extracted table. Any help is highly appreciated! Best Regards, T.


title: Installing GRASS on MINT 14?

text: I'm very new to mint/linux/GRASS and have no clue how to deal with these errors. Can I just change the URL somewhere? Here is the simplest reproducible command to produce a failure: sudo apt-get install mysql-client-core-5.5 Here is the output from my attempt to run/install grass:               mint@mint ~ $ grass     The program 'grass' is currently not installed. You can install it by typing:     sudo apt-get install grass-core     mint@mint ~ $ sudo apt-get install grass-core     Reading package lists... Done     Building dependency tree            Reading state information... Done     The following extra packages will be installed:       gdal-bin libarmadillo3 libdap11 libdapclient3 libepsilon0 libfreexl1       libgdal1 libgeos-3.3.3 libgeos-c1 libgfortran3 libhdf4-0-alt libhdf5-7       liblapack3 libmysqlclient18 libnetcdfc7 libodbc1 libogdi3.2 libpq5 libproj0       libspatialite3 libxerces-c28 mysql-common odbcinst odbcinst1debian2 proj-bin       proj-data python-numpy     Suggested packages:       python-gdal grass-doc grass-gui grass-dev e00compr avce00 gnuplot gpsbabel       gpstrans netpbm python-rpy2 xml2 libhdf4-doc libhdf4-alt-dev hdf4-tools       libnetcdf4 libmyodbc odbc-postgresql tdsodbc unixodbc-bin ogdi-bin       proj-ps-doc python-numpy-doc python-numpy-dbg python-nose python-dev       gfortran     The following NEW packages will be installed:       gdal-bin grass-core libarmadillo3 libdap11 libdapclient3 libepsilon0       libfreexl1 libgdal1 libgeos-3.3.3 libgeos-c1 libgfortran3 libhdf4-0-alt       libhdf5-7 liblapack3 libmysqlclient18 libnetcdfc7 libodbc1 libogdi3.2 libpq5       libproj0 libspatialite3 libxerces-c28 mysql-common odbcinst odbcinst1debian2       proj-bin proj-data python-numpy     0 upgraded, 28 newly installed, 0 to remove and 130 not upgraded.     Need to get 1,043 kB/33.6 MB of archives.     After this operation, 94.8 MB of additional disk space will be used.     Do you want to continue [Y/n]? Y     Err http://archive.ubuntu.com/ubuntu/ quantal-updates/main mysql-common all 5.5.28-0ubuntu0.12.10.1       404  Not Found [IP: 91.189.92.156 80]     Err http://archive.ubuntu.com/ubuntu/ quantal-updates/main libpq5 amd64 9.1.6-1ubuntu1       404  Not Found [IP: 91.189.92.156 80]     Err http://security.ubuntu.com/ubuntu/ quantal-security/main mysql-common all 5.5.28-0ubuntu0.12.10.1       404  Not Found [IP: 91.189.92.200 80]     Err http://security.ubuntu.com/ubuntu/ quantal-security/main libmysqlclient18 amd64 5.5.28-0ubuntu0.12.10.1       404  Not Found [IP: 91.189.92.200 80]     Failed to fetch http://security.ubuntu.com/ubuntu/pool/main/m/mysql-5.5/mysql-common_5.5.28-0ubuntu0.12.10.1_all.deb  404  Not Found [IP: 91.189.92.200 80]     Failed to fetch http://security.ubuntu.com/ubuntu/pool/main/m/mysql-5.5/libmysqlclient18_5.5.28-0ubuntu0.12.10.1_amd64.deb  404  Not Found [IP: 91.189.92.200 80]     Failed to fetch http://archive.ubuntu.com/ubuntu/pool/main/p/postgresql-9.1/libpq5_9.1.6-1ubuntu1_amd64.deb  404  Not Found [IP: 91.189.92.156 80]     E: Unable to fetch some archives, maybe run apt-get update or try with --fix-missing?     mint@mint ~ $


title: how to rename the band-names of a layer stack?

text: I created a layer stack of 20 images (all GTiff-files) using the gdal_merge tool in QGIS, available at:               raster > miscellaneous > merge      The result is a GTiff layer stack of the 20 images, but the band names are band1, band2, band3...band20, looking at the properties. Is there a way to rename the band-names, either afterwards, or during creation of the stack? They are all from different sensors/dates and will be more meaningful to work with other names than band1,band2,band3. Thanks!


title: 3D Interpolation using point cloud in GRASS GIS

text: I've imported a set of 3D points (z,y,z) into GRASS GIS. I'm having trouble: a) visualising my dataset in 3D using nviz or the standard data viewer (I can only get a 2D plan view of my data) and, b) interpolating them into a solid 3D voxel object. My dataset is pretty straightforward and there's only 1520 points. The GRASS user interface is not that user friendly and I cant find a suitable tutorial to follow. Any ideas or pointers would be very much appreciated? I'm sure I should be using the v.vol.rst tool, but when I input my parameters, I get the following error message: ERROR: Column type of wcolumn is not supported (must be integer or double). Any ideas? Cheers Gene Hammond


title: Road Symbology Issues

text: I have preparing road symbology using Arcmap 10.1 version. When I was prepared attribute based symbology as a layer shown on symbols were overlapped. I have attached an example image to show the problem.![enter image description here](http://i.stack.imgur.com/JdFiB.jpg)


title: Writing multiple file geodatabase tables into csv file with arcpy

text: I was asking a question about Merging geodatabase previously, (link: Merge from arcpy.ListTables() producing duplicate rows?) and although the script worked, it took more than 9 hours to run when merging >3000 tables with 1 rows in each table. I was wondering whether instead of merging the tables it would be faster to write the contents of each geodatabase table into a csv file and keep appending tables to it. If so, how to do it?? I now need to merge tables of different geodatabase (it has 39 fields and varying number of rows) and I have a feeling that if I use the script from my previous question it will take even longer than 9 hours. Basically, how to iterate through 3000 tables and write the attribute tables of each to the same csv and keep appending to it? Here's the screenshot of the geodatabase tables I talked about above: ![enter image description here](http://i.stack.imgur.com/Cri5T.jpg)


title: How are people using Python data structures and classes in ArcPy?

text: This question may expose my ignorance on programming but I'm curious about how people are using different python data structures within ArcPy. This page lists the data structures in Python. I understand how lists can be implemented in GIS (list of feature classes, list of feature types, list of data frames, etc). I understand how sets can be used as well (to remove duplicates). How are people implementing tuples, dictionaries, and other data structures within ArcPy? Also, are there other examples of lists and sets which I haven't listed? Furthermore, no doubt, people are creating custom classes in ArcPy. Under what circumstances and situations do you require these? Can you provide examples? Is anyone creating custom classes which inherit from the built-in arcpy classes? I don't require answers to all of these questions, I'm just curious how people are using Python in GIS and what workflows require these customizations.


title: What is the best FME Transformer for Text search and replace in attribute fields

text: I have a number of feature classes which have fields containing text which has somehow become improperly encoded. The fields contain german text which should be either ä, ö, ü, ß but instead contains for instance KindergÃ¤rten instead of Kindergärten ...etc What transformer best allows me to check the field for "Ã¤" and replaces it with ä?


title: Is it possible to display filtered results on OpenStreetMap?

text: is it possible to goto www.openstreetmap.org and do a search for a place like Auckland (New Zealand) or Melbourne (Australia) -but- only show Cities or Neighbourhoods or Suburbs? I understand that it's possible to use filter tags to filter results .. but does this work on the main website? or is it only used for downloading data. if it's for downloading data .. can someone please link to some information about how to do the following   * Neighbourhoods in NZ   * Cities in NZ thanks :) NOTE: This is related to these (here and here) previously asked GIS questions.


title: Event for Selected Feature when Editing?

text: Are there Events in the ArcGIS JavaScript API that I can use to get the selected features being used by the attribute inspector and when features are added to the map? I am using the editor widget in combination with the Attribute Inspector. I have also created a Template picker and used that with the editor widget. In my application, I am editing points, lines, and polygons from a Feature Server using the On Demand Mode. I have seen examples that show selected features using the On Selection mode, but was not pleased with having to add both the map server and the featureServer layers so I am trying to find another way. I would like to be able to start a function that draws a graphic for the point features when they are selected and displayed in the attribute inspector, as well as when the features are first completed but don't know which event is firing. I appreciate any assitance. Thanks, Mele


title: How can I change the coordinate system of gpx-points?

text: I'm working with qgis, v. 2.0.1. I have some points from a gps device (Garmin) and they already are in the .gpx format. Now I have two problems: First, I can't open the gps-Points with the "gpx-tools" plugin by simply choosing "load gpx file". I always have to use the gpx-converter, otherwise the points don't appear in the map window. But actually it should work with "load gpx file", shouldn't it? The second problem is that I want to change the coordinate reference system of the points from WGS84-EPSG:4326 to CH1903/LV03-EPSG:21781. But as soon as I do this, the points disappear...What can I do? Thanks.


title: Contents of Downloadable Data from ESRI website?

text: I'm looking for the simple grey basemap to use in a map I'm working on, but want to download it. It's one of the options with their dynamic data. ESRI has no real order on their download page. Anyone with a list of contents inside these data downloads? I figure the newer ones just have more updated data. ![enter image description here](http://i.stack.imgur.com/mrkUy.png)


title: How to add new data frame(s) to map using ArcPy?

text: I want to add a 2nd data frame to an existing "*.mxd" file. Right now, i see option to add layers and edit properties of "CURRENT" data frame but could not figure out how to add a new one and export outputs from 1st frame to 2nd one. Any suggestions?


title: Is this correct behavior for ST_ConcaveHull (PostGIS)?

text: I'm getting this error on the `ST_ConcaveHull` function call and after reading the docs, I'm unclear as to whether this Error is correct behavior or not. Here is the business end of the query (vx are points):   `select id , ST_ConcaveHull(ST_Collect(vx), 0.80)` And here's the error:               ********** Error **********       ERROR: ExteriorRing: geom is not a polygon       SQL state: XX000       Context: PL/pgSQL function st_concavehull(geometry,double precision,boolean) line 67 at assignment      Since I'm not passing polygons of any sort, I assume the `ST_ConcaveHull` function is creating them internally. I assume that there must be 1 or more grouping of points that don't produce an expected polygon. And this wouldn't surprise me, but I'd rather just get back a null geometry rather than an error that borks the query. Is the behaviour as-designed, and if so, what is the best way to isolate the data that is causing the error (I assume I would need to pre-filter the data somehow)? I can also tell you that   `select id , ST_ConcaveHull(ST_Collect(vx), 0.99)` completes successfully. PostGIS Version: > POSTGIS="2.0.4" GEOS="3.4.2-CAPI-1.8.2 r0" PROJ="Rel. 4.8.0, 6 March 2012" > GDAL="GDAL 1.10.0, released 2013/04/24" LIBXML="2.7.8" LIBJSON="UNKNOWN" > (core procs from "2.0.4" need upgrade) RASTER (raster procs from "2.0.4" > need upgrade) Postgresql Version: > PostgreSQL 9.2.2, compiled by Visual C++ build 1600, 64-bit


title: Testing for selection in ModelBuilder

text: I'm working on a model that will create different buffers for different selections. So I'm using Make Feature Layer and adding an expression to select a subset of features and then using the buffer tool on that subset. But there may be times when a selection doesn't return any features, in which case I wouldn't want to run the buffer. Is there a way in ModelBuilder to test whether or not any features are selected and only run the subsequent step if there are? Alternatively I could export to a Python script and handle this testing there but I'm not quite sure how to accomplish this in Python either.


title: Set Minimum Zoom for View in OpenLayers 3?

text: How can I set a minimum zoom level for a view in OpenLayers 3?               map = new ol.Map({         layers: [imageLayer, vectorLayer],         renderer: 'canvas',         target: 'map',         view: new ol.View2D({             projection: pixelProjection,             center: ol.extent.getCenter(pixelProjection.getExtent()),             extent: pixelProjection.getExtent(),             zoom: 3.5,             maxZoom: 4, // there is option for maxZoom, but no minZoom?         })     });


title: How can I export MBTiles from TileMill 2?

text: In TileMill 1 it was easy to export MBTiles for a selected region with my own style. In TileMill 2 I can only package the style or export vector tiles from a previously imported datasource. How can I export MBTiles with my custom style? Thanks for any hints. I haven’t been able to Google anything relevant.


title: How to create contour from NetCDF?

text: I have downloaded Bathymetry data from **GEBCO**. Its in **NetCDF** format. I want to create contour from it. I tried so many tools to convert the file into raster but couldn't do so. How can I create raster from netCDF?


title: Projection - Some of my data are in odd position

text: Some of my data has this position lon: 181 / lat: 91 So when I try to draw them on the Google map:               var geographic = new OpenLayers.Projection("EPSG:4326");     var mercator = new OpenLayers.Projection("EPSG:3857");          var wfsProtocol = new OpenLayers.Protocol.WFS.v1_1_0({                          url: GeoServerURL,                                      featureType: "Vector_Layer",             featureNS: "http://mapmap.org",             srsName: "EPSG:3857",                     version: "1.1.0",                     extractAttributes: true,                     isBaseLayer: false,                     visibility: true                 });          wfsProtocol.read({                     filter: filter,                               callback: processTheQuery,          scope: [new OpenLayers.Strategy.BBOX({                         ratio: 1                     })],       })                    Layer = new OpenLayers.Layer.Vector("test", {      maxExtent: new OpenLayers.Bounds(-180, -85, 180, 85).transform(geographic,mercator),     displayOutsideMaxExtent: true,      displayInLayerSwitcher: false               })                    function processTheQuery(request) {       if (request.data && request.data.bbox) {     var b = request.data.bbox;        sExt = new OpenLayers.Bounds(b[0],b[1],b[2],b[3]);      Layer.destroyFeatures();       Layer.addFeatures(request.features);     map.zoomToExtent(sExt)      I got this error:               <?xml version="1.0" encoding="UTF-8"?>     <ows:ExceptionReport version="1.0.0"       xsi:schemaLocation="http://www.opengis.net/ows http://myDomain/geoserver/schemas/ows/1.0.0/owsExceptionReport.xsd"       xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ows="http://www.opengis.net/ows">       <ows:Exception exceptionCode="NoApplicableCode">         <ows:ExceptionText>java.lang.IllegalStateException: A transformation exception occurred while reprojecting data on the fly     A transformation exception occurred while reprojecting data on the fly     Latitude 91°00.0'N is too close to a pole.</ows:ExceptionText>       </ows:Exception>     </ows:ExceptionReport>


title: Geoserver transform post request to get?

text: I have the following geoprocessing post request:               <wps:Execute xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.opengis.net/wps/1.0.0" xmlns:wfs="http://www.opengis.net/wfs" xmlns:wps="http://www.opengis.net/wps/1.0.0" xmlns:ows="http://www.opengis.net/ows/1.1" xmlns:gml="http://www.opengis.net/gml" xmlns:ogc="http://www.opengis.net/ogc" xmlns:wcs="http://www.opengis.net/wcs/1.1.1" xmlns:xlink="http://www.w3.org/1999/xlink" version="1.0.0" service="WPS" xsi:schemaLocation="http://www.opengis.net/wps/1.0.0    http://schemas.opengis.net/wps/1.0.0/wpsAll.xsd">     <ows:Identifier>gs:BufferFeatureCollection</ows:Identifier>      <wps:DataInputs>        <wps:Input>         <ows:Identifier>features</ows:Identifier>         <wps:Reference mimeType="text/xml; subtype=wfs-collection/1.0" xlink:href="http://geoserver/wfs" method="POST">         <wps:Body>           <wfs:GetFeature xmlns:medford="http://medford.opengeo.org" service="WFS" version="1.0.0" outputFormat="GML2">           <wfs:Query typeName="medford:firestations"/>           </wfs:GetFeature>         </wps:Body>        </wps:Reference>       </wps:Input>       <wps:Input>         <ows:Identifier>distance</ows:Identifier>         <wps:Data>          <wps:LiteralData>0.024</wps:LiteralData>         </wps:Data>       </wps:Input>     </wps:DataInputs>     <wps:ResponseForm>       <wps:RawDataOutput mimeType="text/xml; subtype=wfs-collection/1.0">       <ows:Identifier>result</ows:Identifier>      </wps:RawDataOutput>     </wps:ResponseForm>     </wps:Execute>      How can I transform it to get request? Thanks!


title: How to create shortest route and multiple routes service using openlayers?

text: I created a digital database of many important features and locations of my city. and I want to create a web application in which shortest distance between two locations can be calculated and displayed. Is it possible with openlayers as client scripting. My shapefiles are in postgis. What I am expecting is :   * The database contains all roads, point features etc.   * When the user types a location, the required location name can be listed out from the database and then user can similarly select another location   * Then my technique is to find the nearest road points for the two locations (may be using buffer method) and then    * finding the shortest and possible routes between the two road points along the road vector features. But I don't know how to write the code and algorithms.


title: When we type an address in search bar, in some web mapping applications like google maps, a list of words automatically displayed from the database.

text: What is actually happening behind the code and how should be the code? when we search for a specifc word, will it look in the database (here it is postGIS) and lists the similar words.


title: How do I rotate the polygon about an anchor point using python script?

text: I know there would be no tool in arcpy module as such but mignt have in third party modules. I however tried taking the vertices of the polygon and moving using the below formula. but did not get the correct output coordinate values and polygon thus generated is wierd. I have tried this function               def RotateAxis(AnchorX,AnchorY,inputx,inputy, WindDirection):         x = inputx - AnchorX         y = inputy - AnchorY         resultx = (x * math.cos(WindDirection)) - (y * math.sin(WindDirection)) + AnchorX         resulty = (x * math.sin(WindDirection)) + (y * math.cos(WindDirection)) + AnchorY         return (resultx,resulty)      where Anchor is the origin and WindDirection is the angle which I tried taking in radians as well as degrees but could ot get correct results. I tried with projected coordinates that didn't work either. Is there any suggestion?


title: Searching for GIS data sources relating to the Soviet GULAG system

text: I'm desperately searching for GIS data related to the Soviet prison camp system, the GULAG. There are quite a few very good maps out there, but I cannot seem to find data in order to make my own. Thinking of locations of the camps, administrative divisions of the penal system etc. Any ideas, anyone?


title: How to use SRTM30 Plus data with QGIS or other GIS tool?

text: I am trying to access the new Bathymetry data published here: http://topex.ucsd.edu/WWW_html/srtm30_plus.html Is there a way QGIS or similar tool can read these files. Nothing seems to be working, but I am very new to using GIS tools and their rainbow of formats. Any examples or ideas would be appreciated. Thanks Eric


title: OpenStreetMap display as simple as Google Maps?

text: I'm glad OpenStreetMap exists at all, but I find their display has way too much details. Here's an exemple of the same location: ![enter image description here](http://i.stack.imgur.com/5nvoO.png) Is there a way to view OSM contents with the same level of details as Google Maps, especially on smartphones? Thank you.


title: Which landmark to detect a CBD?

text: Anyone knows what shapefile (ArcGIS) available to detect a CBD (Central Business District), nationwide? Some use city hall as the landmark for a CBD. Is it applicable to use Courthouse? Where can I get these data?