title: Using Union on multiple shapefiles to output single result

text: Using QGIS 2.2 and Python, I would like to use the Union algorithm on multiple shapefiles (contained in a folder _TEST_ ) which would then produce single shapefile containing all the data. Currently the following code only outputs the union of the first and last file, how can I modify it so that the output of 2 layers is saved and then use that output with another layer etc? UPDATE: The shapefiles are being processed in terms of the `clip_function` however I can't load the temp files when running the `qgis:union` algorithms. Also, an error occurs when it tries to output the Result.               ##Model=name     ##Select_folder=folder     ##Cellsize=number 1000     ##Result=output vector          import os     import glob     import tempfile     from qgis import utils     from qgis.core import QgsVectorLayer, QgsMapLayerRegistry, QgsVectorFileWriter, QgsCoordinateReferenceSystem          #   Create grid with defined parameters     outputs_1=processing.runalg("qgis:creategrid", Cellsize, Cellsize, 24108, 18351.157175, 258293.802316, 665638.226408, 1, 'EPSG:7405', None)          def clip_function():          #   Set directory and search for all .shp files         os.chdir("C:\Users\TEST\DATA\\")         for fname in glob.glob("*.shp"):           #   Clip .shp files within directory to grid and save files as temp             outputs_2=processing.runalg("qgis:clip", outputs_1['SAVENAME'], fname, None)             outputs_3=processing.runalg("qgis:fieldcalculator", outputs_2['OUTPUT'], 'Rank_' + fname, 1, 10, 0, True, 1 , None)                      def binary_alg_wrapper(list_of_files,final_file_name):         run_count=0         temp_path=tempfile.gettempdir()         temp_input_name="temp_input_layer"         temp_output_name="temp_output_layer"         temp_input_layer_filename=str(os.path.join(temp_path,temp_input_name+".shp"))         temp_output_layer_filename=str(os.path.join(temp_path,temp_output_name+".shp"))         for file_name in list_of_files:             EPSG_Code_for_file=4326             run_count+=1             if run_count==1:                 #Copy the first file into a temp location                 temp_output_layer = QgsVectorLayer(file_name, "temp_output_layer", "ogr")                 Add_result = QgsMapLayerRegistry.instance().addMapLayer(temp_output_layer)                 write_error = QgsVectorFileWriter.writeAsVectorFormat(temp_output_layer,temp_output_layer_filename,"system",QgsCoordinateReferenceSystem(EPSG_Code_for_file), "ESRI Shapefile")             QgsMapLayerRegistry.instance().removeMapLayer(temp_output_layer.id())         if run_count>1:             #Copy the result of the last algorithm into a temp location             old_output_layer = QgsVectorLayer(temp_output_layer_filename, "input_layer",     "ogr")             Add_result=QgsMapLayerRegistry.instance().addMapLayer(old_output_layer)             write_error = QgsVectorFileWriter.writeAsVectorFormat(old_output_layer,temp_input_layer_filename, "system", QgsCoordinateReferenceSystem(EPSG_Code_for_file), "ESRI Shapefile")             QgsMapLayerRegistry.instance().removeMapLayer(old_output_layer.id())                  #Read in a new file and the file stored in the temp location             input1_layer = QgsVectorLayer(temp_input_layer_filename, "input1_layer", "ogr")             input2_layer = QgsVectorLayer(file_name, "input2_layer", "ogr")                  #Run the binary algorithm and store the result into a temp location or output location             if run_count==len(list_of_files):                 processing.runalg("qgis:union", input1_layer, input2_layer, temp_output_layer_filename)             else:                 processing.runalg("qgis:union", input1_layer, input2_layer, final_file_name)                     QgsMapLayerRegistry.instance().removeMapLayer(input1_layer.id())             QgsMapLayerRegistry.instance().removeMapLayer(input2_layer.id())               path_res = tempfile.gettempdir()     file_list = [shp for shp in glob.glob(path_res + "*.shp")]     Result=(os.path.join(path_res, "result.shp"))          if Select_folder:         clip_function()         binary_alg_wrapper(file_list, Result)


title: Feature count of current features in Extent based on attribute?

text: I have one points layer in a cached map service, symbolised with a unique renderer on one field.   e.g. petrol stations, symbolised by the company who owns them. I want to create a fairly basic application in the Esri Javascript API (or potentially with the Google Maps API, and this). One of the main components I need to have on this application is a dynamic legend that shows a feature count for the current extent.   e.g. Shell - 34, BP - 21, Texaco - 7, etc Has anyone seen any code that could get me started on this or can provide some references/pseudocode on how one would go about this?


title: Problems with WMS QGIS 1.8

text: Hi I have loaded WMS layers previously but i seem to have hit a brick wall. I am trying to load from either the Geological Survey Irelands website (near the bottom of the page): http://www.gsi.ie/mapping.htm or using the Tellus Project: http://www.tellusborder.eu/Data+Downloads/ I enter the location eg http://spatial.dcenr.gov.ie/wmsconnector/com.esri.wms.Esrimap/GSI_Bedrock_Geology? in server search, the search but neither are presenting layers to load. What am i doing wrong?


I am using Google maps for my application. For getting some custom private map data, I have set up a WMS using Mapserver. How can I add this Custom UMN Mapserver WMS layer over Google API (without using Openlayers)?


title: What extensions are included in ArcGIS 10.2 Basic Version

text: I wanted to purchase ArcGIS 10.2. i used its trial version for last 2 month. now i am thinking to have Basic version. My work is based on Spatial analysis. so i dont need 3d and business analyst. for that i wanted know what extension will come with Arcgis 10.2 basic version.?. what are the functions available? and Spatial analyst extension come with ArcGis basic? or i have to purchase spatial analyst separately if anyone have info about this plz tell me.


title: Creating a list of unique street names from a list of address

text: I would like to generate a unique list of street names from a list of UK addresses. But obviously there will be a number of streets with the same name. Which is fine, but I want to distinguish between streets that are genuinely different. On AtoZ, the street names on the index is written as:               Bell street, RH1     High Street, RH1       and so on. I have yet to come across a street name with more than one outcode. So based on this,               34, Bell street, RH1 6HH and      36, Bell street, RH1 6HJ      would refer to the same street, but say you had               34, Bell street, RH1 6HH and     36, Bell street, RH2 6HJ      can I be certain that these are two different streets? Thanks


title: How can we remove the ESRI logo from ArcGIS Viewer for Flex 2.2

text: Logo layer is set to be visible by default. Is there a way to turn it off? Thank you


title: How can I transform GCS_North_American_1983 to GCS_WGS_1984?

text: I posted a thread today describing a shift in some of the NRC data (not all). Apparently, this shift has to do with conversion between Geographic coordinates, latitude and longitude and Transverse Mercator Grid coordinates. Can someone please provide a guide on how to do this transformation in ArcGIS? I would like to transform that data to either `Clarke_1866_Lambert_Conformal_Conic` or `GCS_WGS_1984`. Source Data Coordinate System:               GCS_WGS_1984     Angular Unit: Degree (0.017453292519943299)     Prime Meridian: Greenwich (0.000000000000000000)     Datum: D_WGS_1984     Spheroid: WGS_1984     Semimajor Axis: 6378137.000000000000000000     Semiminor Axis: 6356752.314245179300000000     Inverse Flattening: 298.257223563000030000


title: Python ArcGIS 10.0 Expanding select layer by location from point (walking selection)

text: I have an old ArcView 3 tool which I've been tasked with converting to Python. I'm breaking it down into stages, but there's a "walking selection" I can't trap. To put the query in context, the user selects a site (point data) - selection A. All sites within X kilometres of selection A are then selected (selection B), then all sites within X of kilometres of selection B are selected, and so forth until no new selections are made. I've tried to use a While... loop with a Get Count to no avail. The problem I think is because the Get Count will always have a value, so I think I need to do something like Get Count <> Previous Get Count?               #Import ArcPy site package     import arcpy          # Overwrite files     arcpy.env.overwriteOutput = True           try:          # Make a feature layer with all the farms          arcpy.MakeFeatureLayer_management(r"shapefile.shp", "Farms")              # Make a feature layer containing only the farm of interest          arcpy.MakeFeatureLayer_management(r"shapefile_select.shp","Farms_select")              # Apply a selection to farms layer         fcFarms = arcpy.SelectLayerByLocation_management ("Farms","WITHIN_A_DISTANCE","Farms_select", "2000 Meters")              # Apply a further selection to farms layer         fcFarms2 = arcpy.SelectLayerByLocation_management("Farms","WITHIN_A_DISTANCE",fcFarms, "2000 Meters", "ADD_TO_SELECTION")              while getCount > 0:             fcFarms3 = arcpy.SelectLayerByLocation_management("Farms","WITHIN_A_DISTANCE", fcFarms2, "2000 Meters", "ADD_TO_SELECTION")             fcFarms4 = arcpy.SelectLayerByLocation_management("Farms","WITHIN_A_DISTANCE", fcFarms3, "2000 Meters", "NEW_SELECTION")             getCount = int(arcpy.GetCount_management(fcFarms4))          except:         print arcpy.GetMessages()


title: Creature individual Feature Classes from Features in a Feature Class

text: I have a feature class called Tracks2012. In the feature class I have 4 fields (MMSI, Ship_Type, Vessel_Name, IMO) as well as the usual Object_ID, Shape and Shape_Length field (as these features are polylines) Essentially what I am trying to do is create an individual Shapefile from each feature using arcpy. I'm not all that familiar with arcpy, and as oppose to cluttering up this question with what I attempted, I'll show the first little part where I set up where I am going to put my feature classes.               >>> import arcpy     >>> arcpy.env.workspace = "c:/GIS/mmsi.gdb"      Once I have this completed, I need to apply a buffer of 10km to each individual polyline shapefile followed by a reclassification of each shape file to have either a 0 or a 1 in relation to the presence of the buffered polyline (which will now be a polygon as a result of the output of the buffer). I am not sure whether or not all of these operations should be maintain within one arcpy script or several. I imagine this all seems quite confusing (the way I have written it up), and I am new to arcpy, so please before down voting my question perhaps suggest alternative to improving its structure. Thanks!


title: Tiles causing problems with dasharray

text: I'm having a problem with a stroke-dasharray in Geoserver SLD. When previewing in single-tile mode, everything looks like I want it to, but when I preview in tile-mode I get gaps on the border of the tiles. Left image is tiled, right image is the (preferred) single-tile: ![Tiles](http://i.stack.imgur.com/zuWAZ.png) ![Single- tiled](http://i.stack.imgur.com/5ur03.png) You can see the middle tile here, for some reason it doesn't render the dasharray on either end, even though I have dashoffset set to 0. ![Middle tile](http://i.stack.imgur.com/0nRRH.png) Any idea how to fix this problem? I want to use tiles.


title: Is it possible to style a part of a label?

text: My label consists of two different things:               <ogc:PropertyName>user</ogc:PropertyName>     <ogc:PropertyName>room</ogc:PropertyName>      I would like to style the two parts differently, for example have the room number in bold. Is it possible?


title: Updating a column with line direction in MapInfo

text: I was wondering if there is a way to update a column in MapInfo with the digitised line direction. I have two TAB files, one is complete lines and the other is sections of the same lines (both have a unique ID). I was wondering if it would be possible to add a column in both files with the direction of the line so I can check if the directions are identical, e.g: ![enter image description here](http://i.stack.imgur.com/aPcKQ.png) Thanks, Alan


title: statistics by rectangle (raster image)

text: I have a shp file with some rectangles and I want to extract statistics within these areas from a raster image. I prefer ArcMap 9.3.


title: What kind of improvement can be done for raster tiles for mobile usage

text: We generate raster map tiles and deploy them in the server, then we use javascript in the client to request the tiles and re-organize them to map.(like google map and Openstreetmap). However this is the PC solution. And the tile size is 256x256(pixel) with 96dpi. Now we want to build apps in android and iOS. However we found that the same 256x256(pixel) tile will have a different visual effect in devices with different resolution and the tile is almost unreadable in a high-resolution device(The text and labels in the tile is too small). Then I wonder how do you make the tiles readable in the mobile device, is there any way to improve them?


title: ArcSDE - where to start?

text: As I long ago figured out that I need more complex query-asking abilities (other than the standard select provided by arcgis), and as I will someday have to admin a small scale SDE database that will import data from my corporate's big and country-wide SDE, I wanted to set up a small workgroup SDE on a standalone computer so I could deal with it a bit, to get the hang of dealing with a sde. But where do I start? I can't find a good manual to start with - The Arcgis 9.3.1 Help(on which I'm working) isn't much help. I know I need to Install SQL server, but I can't understand any of the other steps (creation, import of data and such). Do you know of any good guides and resources to guide me with the arcgis 9.3.1 SDE from scratch? (Oh, and does SDE 10 differ a lot from SDE 9.3.1?)


title: Converting MapInfo TAB to GeoJSON

text: Is there an easy way to convert MapInfo TAB files into GeoJSON? I am working on a web-based mapping system which will use GeoJSON to transfer data. I also use MapInfo and would like to be able to convert files from MapInfo into GeoJSON. I have tried the Universal Translator that comes with MapInfo but it doesn't support GeoJSON. If there isn't a way to do it I will make a MapInfo plug-in myself (and probably post it online). I just don't want to waste my time making this if there is already one out there. **Update 2012-03-19:** You can now export to GeoJSON through Quantum GIS.


title: Why does the buffer tool keep producing ellipsoid when I want a circle?

text: I want to create a 50 mile radius around a point. I've got the point as a layer. I've run the buffer tool found in Vector -> Geoprocessing Tools and set the distance to 0.5. Every time I run the tool on either a point or polygon I keep getting a ellipsoid 50 miles wide by 70 miles high. I've also tried this on both QGIS and ArcGIS and both produce the same problem. I'm sure I'm over looking something. If someone could point it out what I am missing that would be great. Thank you for your time.


title: Why does ogr2ogr silently skip some rows?

text: I have a CSV file called `office_locations.csv` that contains point data. I want to load it into a SQL Server geospatial database using ogr2ogr. The CSV file is encoded in UTF-8. It contains a header and 14902 records. The first lines look like this:               Id,Name,Lat,Lng,CountryCode,Address,CityName,Airport,AirportCode,RailwayStation,PostalCode,StateCode,StateProv,Country_Id     "7188","Tirana Airport","41.42108838","19.71271276","AL","Tirana Airport Muhamet Gjollesha Str., Muhamet Gjollesha Str., Tirana","Tirana","1","TIA","0","","","",0     "30768","Tirana Downtown","41.332","19.832","AL","Rruga E Durresit. Nr 61, Tirana","Tirana","0","","0","","","",0     "52400","Sheraton Italia Square","0.0","0.0","AL","Square Italia Hotel Sheraton , Tirana,","Tirana","0","","0","","","",0     "4650","Escaldes","42.511","1.548","AD","Avda. Miquel Mateu, 25, Escaldes-engordany , Andorra, Ad700","Andorra","0","","0","AD700","","",1     "152576","Andorra","42.508","1.522","AD","Avda D Enclar, 142, Andorra La Vella, Ad500","Andorra La Vella","0","","0","AD500","","",1      I'm working on a Windows 7 workstation. My GDAL version is 1.9.2. **ogr2ogr import skips records** To import CSV data using ogr2ogr you need to create a VRT file. I created one that looks like this:               <OGRVRTDataSource>         <OGRVRTLayer name="office_locations">             <SrcDataSource>office_locations.csv</SrcDataSource>             <GeometryType>wkbPoint</GeometryType>             <LayerSRS>WGS84</LayerSRS>             <GeometryField encoding="PointFromColumns" x="Lng" y="Lat"/>         </OGRVRTLayer>     </OGRVRTDataSource>      When I use the VRT file, ogr2ogr reads and writes only 4977 features:               $ ogr2ogr --debug ON -overwrite -f MSSQLSpatial  "MSSQL:server=.;database=OfficeImport;trusted_connection=yes" office_locations.vrt     OGR: OGROpen(office_locations.vrt/02F28328) succeeded as VRT.     OGR_MSSQLSpatial: EstablishSession(Connection:"server=.;database=OfficeImport;trusted_connection=yes")     ODBC: SQLDriverConnect(DRIVER=SQL Server;server=.;database=OfficeImport;trusted_connection=yes)     OGR: OGROpen(MSSQL:server=.;database=OfficeImport;trusted_connection=yes/02F3CFB8) succeeded as MSSQLSpatial.     OGR: OGROpen(office_locations.csv/02F3C808) succeeded as CSV.     OGR_MSSQLSpatial: Using column ogr_fid as FID for table office_locations.     OGR2OGR: 4977 features written in layer 'office_locations'     ODBC: SQLDisconnect()     VRT: 4977 features read on layer 'office_locations'.     CSV: 4977 features read on layer 'office_locations'.      Two thirds of the of the records have been silently skipped! **Finding the missing records** The id column in the CSV contains a unique id for each office. I can sort by that to find the missing records. To order the rows in the database I used a command like this:               $ sqlcmd -d OfficeImport -Q "SELECT id, name FROM office_locations ORDER BY CAST(id AS INT);" | more     id   name               ---- -------------------------------------------------------     11   Dublin - Airport     21   Dublin - Finglas     36   Brussels - Charleroi Airport     42   Crete - Chania Marna     79   Lesvos - Airport     83   Corfu - New Port     [...]      To order the rows in the spreadsheet I used Excel 2010's Sort feature on the Data tab. I discovered that office 47 is the first missing record:                "47","Fuerteventura - Airport","28.452","-13.87","ES","Aeropuerto Fue L.7, El Matorral.pto Del Rosario, Fuerteventura, 35610","Fuerteventura","1","FUE","0","35610","","",136      **Isolating the skipped row** I created a new CSV file called `test_import.csv` that contains just office 47 and a header:               Id,Name,Lat,Lng,CountryCode,Address,CityName,Airport,AirportCode,RailwayStation,PostalCode,StateCode,StateProv,Country_Id     "47","Fuerteventura - Airport","28.452","-13.87","ES","Aeropuerto Fue L.7, El Matorral.pto Del Rosario, Fuerteventura, 35610","Fuerteventura","1","FUE","0","35610","","",136      And a new VRT called `test_import.vrt` that looks like this:               <OGRVRTDataSource>         <OGRVRTLayer name="test_import">             <SrcDataSource>test_import.csv</SrcDataSource>             <GeometryType>wkbPoint</GeometryType>             <LayerSRS>WGS84</LayerSRS>             <GeometryField encoding="PointFromColumns" x="Lng" y="Lat"/>         </OGRVRTLayer>     </OGRVRTDataSource>      ogr2ogr debug output shows that 1 feature was read and written:               $ ogr2ogr --debug ON -overwrite -f MSSQLSpatial  "MSSQL:server=.;database=Argus;trusted_connection=yes" test_import.vrt     OGR: OGROpen(test_import.vrt/03488318) succeeded as VRT.     OGR_MSSQLSpatial: EstablishSession(Connection:"server=.;database=Argus;trusted_connection=yes")     ODBC: SQLDriverConnect(DRIVER=SQL Server;server=.;database=Argus;trusted_connection=yes)     OGR_MSSQLSpatial: Using column ogr_fid as FID for table argus_locations.     OGR_MSSQLSpatial: Using column ogr_fid as FID for table location.     OGR_MSSQLSpatial: Using column ogr_fid as FID for table location_utf8.     OGR_MSSQLSpatial: Using column ogr_fid as FID for table office_locations.     OGR_MSSQLSpatial: Using column ogr_fid as FID for table test_import.     OGR: OGROpen(MSSQL:server=.;database=Argus;trusted_connection=yes/0349CFA8) succeeded as MSSQLSpatial.     OGR: OGROpen(test_import.csv/0349D868) succeeded as CSV.     MSSQLSpatial: DeleteLayer(test_import)     OGR_MSSQLSpatial: Using column ogr_fid as FID for table test_import.     OGR2OGR: 1 features written in layer 'test_import'     ODBC: SQLDisconnect()     VRT: 1 features read on layer 'test_import'.     CSV: 1 features read on layer 'test_import'.      This time office 47 is imported, but last time it was skipped. Why is this row skipped in the full file but imported in the test file? Why doesn't ogr2ogr acknowledge that any rows were skipped? **How many columns does each table have?** MappaGnosis asked: "How many columns does the single office (47) import's attribute table have and how many does the original import have?" The input CSV file has 14 columns:               $ powershell -Command "(Import-Csv -Path 'office_locations.csv' | Get-Member -Type NoteProperty).Count"     14      The original import table is called `office_locations`, and the single office import table is called `test_import`. Each has 16 columns:               SELECT       OBJECT_NAME(object_id) AS table_name,       COUNT(*) AS column_count     FROM sys.columns     WHERE object_id IN (OBJECT_ID(N'test_import'), OBJECT_ID(N'office_locations'))     GROUP BY object_id;          table_name         column_count     -----------------  ------------     office_locations   16     test_import        16      ogr2ogr created one column for each of the fields in the CSV files, and generated two extra columns `ogr_fid` and `ogr_geometry`. **Importing just the ids and the coordinates** Darren Cope suggested: "Can you create a second .csv with all records, but strip out all of the columns except the ID and lat/lon columns? Try to import that, and see if all records come in. Sometimes there is an improperly delimited row somewhere that will confuse things, and is hard to pick out visually..." I used this PowerShell command to extract just the `id`, `lng`, and `lat` columns from the `office_locations.csv` file:               $ powershell -Command "(Import-Csv -Path 'office_locations.csv') | Select Id, Lng, Lat | Export-Csv -NoTypeInformation -Encoding UTF8 -Path 'office_minimal.csv'"      It created a new CSV file called `office_minimal.csv`. It contains the minimal amount of data to locate and identify each office. The first lines look like this:               "Id","Lng","Lat"     "7188","19.71271276","41.42108838"     "30768","19.832","41.332"     "52400","0.0","0.0"     "4650","1.548","42.511"     "152576","1.522","42.508"      I created a corresponding VR file called `office_minimal.vrt`. It looks like this:               <OGRVRTDataSource>         <OGRVRTLayer name="office_minimal">             <SrcDataSource>office_minimal.csv</SrcDataSource>             <GeometryType>wkbPoint</GeometryType>             <LayerSRS>WGS84</LayerSRS>             <GeometryField encoding="PointFromColumns" x="Lng" y="Lat"/>         </OGRVRTLayer>     </OGRVRTDataSource>      Using the new 'minimal' files, ogr2ogr imports all 14902 rows:               $ ogr2ogr --debug ON -overwrite -f MSSQLSpatial  "MSSQL:server=.;database=OfficeImport;trusted_connection=yes" office_minimal.vrt     OGR: OGROpen(office_minimal.vrt/03388328) succeeded as VRT.     OGR_MSSQLSpatial: EstablishSession(Connection:"server=.;database=OfficeImport;trusted_connection=yes")     ODBC: SQLDriverConnect(DRIVER=SQL Server;server=.;database=OfficeImport;trusted_connection=yes)     OGR_MSSQLSpatial: Using column ogr_fid as FID for table office_locations.     OGR_MSSQLSpatial: Using column ogr_fid as FID for table test_import.     OGR: OGROpen(MSSQL:server=.;database=OfficeImport;trusted_connection=yes/0339CFB8) succeeded as MSSQLSpatial.     OGR: OGROpen(office_minimal.csv/033A0F48) succeeded as CSV.     OGR_MSSQLSpatial: Using column ogr_fid as FID for table office_minimal.     OGR2OGR: 14902 features written in layer 'office_minimal'     ODBC: SQLDisconnect()     VRT: 14902 features read on layer 'office_minimal'.     CSV: 14902 features read on layer 'office_minimal'.


title: Postgis objects inside Leaflet circle

text: How can I retrieve database objects from postgis database which are contained inside Leaflet circle object? I assume circle is only a point with some radius and not a polygon. I wonder can it be somehow converted to the polygon so I can read objects inside it from postgis database. Thanks!


title: How do I copy the coded value for a field and not the raw value?

text: I am working with ArcGIS 10 and I have added a field in a feature class and am attempting to copy the text from another field's coded values, but the field calculator is simply copying the "code" text from the source field and not the "value" text from the domain rule. Thanks in advance!


title: How to leverage packages like Geocortex for small business GIS professionals?

text: Can anyone talk about how packages like Geocortex could be utilized by small business or working from home GIS professionals? I would like to specifically know if Geocortex replaces any other packages or sets of tools that would make running your small business platform that much easier.


title: Storing Photos in layer in ArcPad 10.2

text: I am using ArcPad 10.2 to do some field data collection. On of the functions I would like to be able to do is associate a photo (taken in the field with the built in camera) to an existing point layer. My assumptions is that I will need an empty field to store the photo, but I am unsure of what that field needs to be as a data type. Any help and info would be appreciated. Thanks...


title: mmqgis plugin not showing up in plugin dropdown menu

text: I have installed the plugin from the plugin dialog box, but I am unable to figure out exactly how to run the plugin.


title: How to sample a raster using the Point Sampling Tool in QGIS 2.0?

text: I have been trying to sample raster using the Point Sampling Tool in QGIS 2.0. Though the new shape file is getting generated it does not have the values unfortunately. How do I take care of this? This tool is really necessary for my work but somehow its not working.


title: How do I know if the Model Builder has updated.in ArcCatalog?

text: How do I know if the Modelbuilder has updated ? So when I worked on Modelbuilder and then make some modified and added or subtracted tools. Sometimes you would forgot when you come back to it in a few days later and you have two same Model Builder and you don't know which one is the most recently ? ![enter image description here](http://i.stack.imgur.com/lodDE.png)


title: How to compare files in two folders?

text: I'm pretty new user of ArcGIS. My intention is compare files (csv) from one folder with files in another folders. In those folders exist files with the same name but in second folder files have different endings. First Folder:   * almukantarat.csv    * almukantarat1.csv Second Folder:   * almukantarat_A.csv   * almukantarat1_A.csv I don't know, maybe it can be done in Python or only in ModelBuilder (using Compare files?). Please, could You help me how script or model should look like? I have ArcGIS Desktop 10.0. Thanks


title: CRS for submission to Minnesota DOT

text: Does anyone have experience using QGIS to submit shape files to the Minnesota DOT (MNDOT)? I am using EPSG:26915 but got rejected due to an error in the CRS. MNDOT Specs call for "UTM coordinates, Zone 15E extended, NAD83, GRS80, horizontal units = meters, vertical units (if used) = feet." Any help would be appreciated. Thanks.


title: POI database of tourist attractions, or POI database with popularity index

text: I'd like to build a list of popular points of interest and tourist attractions for many/every city in the world. Is there any sort of dataset that I can use for this which either contains only popular attractions, or something like the OSM data with a popularity index that I can filter through myself? Open or Free is preferable. I've tried playing with FourSquare's API, but it is against their TOS to save any of their data for longer than just caching. It was also common for things like a Starbucks to appear in the list of most popular locations, but this wouldn't be too troublesome to work around. (I guess because more people want to check-in at Starbucks than a Museum). I'm looking through OSM at the moment, but there doesn't seem to be any indication of the popularity or importance of a POI. The tags are nice, but not very helpful to me if there is no way to differentiate Disney World from a touristy taco-shack. Thanks in advance for any help or direction.


title: Openlayers + WMS layer. Use scale or resolution + DPI

text: I have a problem with displaying map with WMS layer in Openlayers. I'm able to display it but the image of the map is not sharp especially at higher zooms. I have a documentation from the map provider and below is part of description: The data scale set is: > 1:20 mil., 1:6 200 000, 1:3 000 000, 1:1mil., 1:500 000, 1:200 000, 1:75 > 000, 1:25 000, 1:12 500 , 1:6 250 For the best quality of displayed data while zooming is recommended to count with recalculation 96dpi/200dpi (Map DPI = 200), it means 0,48 x default scale set. I'm not sure if this DPI setting that is wrong in my code. I found out that Openlayers standard DPI is 72 DPI. Can you suggest anything how can I fix it? ## WMS getcapabilities:               <SRS>EPSG:32633</SRS>     <LatLonBoundingBox minx="-8.060897" miny="32.143274" maxx="56.352488" maxy="66.324824"/>     <BoundingBox SRS="EPSG:32633" minx="-1700000" miny="3800000" maxx="2300000" maxy="8000000"/>      ## My code:               <!DOCTYPE html>     <html>       <head>         <title>WMS map</title>         <script src="http://openlayers.org/api/OpenLayers.js"></script>          <script type='text/javascript'>                var map; function init() {          {         var mapoptions = {           maxExtent: new OpenLayers.Bounds(-1700000, 3800000, 2300000, 8000000),           numZoomLevels:14,           scale: [20000000, 6200000, 3000000, 1000000, 500000, 200000, 75000,            25000, 12500, 6250],           maxResolution: "auto",           projection: "EPSG:32633",           displayProjection: new OpenLayers.Projection("EPSG:32633")              }};         map = new OpenLayers.Map( 'map',mapoptions );                       var wms = new OpenLayers.Layer.WMS( "layer_name" ,'https://wms.dot.com',         {layers: 'layer_name',        format: "image/png"},         {}      );      map.addLayer(wms);      if(!map.getCenter()){ map.zoomToMaxExtent(); } }      </script>     </head>          <body onload='init();'>         <div id='map' style='width: 1000px; height: 1000px;'>     </div>     </body>     </html>


title: Linking Geoportal Map Preview to a Geoserver service?

text: I am very new to Geoportal and Geoserver I have an installation of Geoportal. In the metadata of a resource I can indicate a URL link. I have seen that depending of the linked data the "Preview link" is shown. I also have a Geoserver serving WCS,WMS and WFS services. What do I need to put in the link URL of a Geoportal resource so when I click preview the request is sent to Geoserver and a shape or grid is rendered in the preview Window? Or if it has to be WMS which image is supported by the preview window? Please let me know if I am totally confused about linking Geoportal and Geoserver. Many thanks, Carlos.


title: Splitting a line into two in OpenLayers

text: I have some WFS layers in OpenLayers, containing lines. I try to make a tool to split these lines and add a symbol in the breakpoint. Then save them to db(postgis). How can I implement it? Thanks.


title: what is the syntax for a Maximum operator in ArcGIS

text: The requested maximum operator is a function that produces a new raster with the highest values on a cell-by-cell basis from an array of coincident input rasters. It is the kind of Maximum Support classification performed by the 'HARDEN' tool of the Idrisi software. The raster calculator spatial analyst tool should be used. I tried max("raster_A","Raster_B","Raster_C","Raster_D") But got the following error: 'exceptions.ValueError': The truth value of the raster is ambiguous Thanks in anticipation


title: ArcGIS Package information?

text: I'd like to know how deal with the ArcGIS Package information?   I have ArcGIS 9.3.1 and would like to know how to work with map packages, and layer packages.   Also how do I access it?


title: How do I remove a join (between a standalone table and a feature layer) with ArcObjects?

text: I already know how to do a table join with ArcObjects, via the `IDisplayTable`, `IMemoryRelationshipClassFactory` etc. interfaces. Here are some links to resources that have example code for this:   * **"Advanced" table join options in ArcObjects** _(ESRI ArcGIS Desktop discussion forums archive)_   * **Export a join table to arcmap** _(question here on GIS Stack Exchange)_ What I need to know is, how do I undo, i.e. remove such a join? * * * Once I've joined a standalone table to a feature layer, I end up with references to:   * an `IRelationshipClass` (resulting from the join operation);   * an `ITable`/`IStandaloneTable`/`IDisplayTable` (the table that was joined to the feature layer); and   * an `IFeatureLayer` (the feature layer to which the table was joined). Do these interfaces, or rather the objects behind them, allow unjoining at all? I've seen that `IRelationshipClass` has various `DeleteRelationship…` methods, but I can't see how they could work towards that end.


title: How to use 'remap' in reclassify tool code in ArcGIS Python Toolbox

text: I am making a reclassify tool in ArcGIS python toolbox (.pyt) ... I am able to use the remap datatype but am not able to parse the information present in Remap in the Reclassify function's third argument ... Code is below:               import arcpy, string     import numpy     from arcpy import env     from arcpy.sa import *          class Toolbox(object):     def __init__(self):         """Define the toolbox (the name of the toolbox is the name of the         .pyt file)."""         self.label = "Toolbox"         self.alias = ""              # List of tool classes associated with this toolbox         self.tools = [Remap]          class Remap(object):     def __init__(self):         self.label = "Remap"         self.description = ""         self.canRunInBackground = False          def getParameterInfo(self):         feature_class = arcpy.Parameter( displayName="Input Raster ", name="input1", datatype="Raster Layer", parameterType="Required", direction="Input")          field = arcpy.Parameter( displayName="Input Field ", name="input2", datatype="Field", parameterType="Required", direction="Input")          remap = arcpy.Parameter( displayName="ReMap ", name="remap", datatype="remap", parameterType="Required", direction="Input")              out_class = arcpy.Parameter( displayName="Output Raster ", name="ouput1", datatype="Raster Layer", parameterType="Required", direction="Output")               field.value = "Value"           field.enabled = True              remap.parameterDependencies = [feature_class.name,field.name]              params = [feature_class,field,remap,out_class]         return params          def isLicensed(self):         return True     def updateParameters(self, parameters):              #if parameters[0].value:             #parameters[1].parameterDependencies = [parameters[0].name]             #parameters[2].parameterDependencies = [parameters[0].name,parameters[1].name]         return           def updateMessages(self, parameters):         return          def execute(self, parameters, messages):         inRaster = parameters[0]         outRaster = parameters[3]         remap2 = parameters[2]         arcpy.gp.Reclassify_sa(inRaster,"Value",remap2,outRaster,"DATA")         return


title: What's your view on the direction of GIS?

text: I am a long standing GIS developer and consultant. I am often asked, is there still a place for Desktop GIS, to which I, _currently_ , answer, of course; if there are jobs specific to individuals, you shouldn't be developing enterprise solutions for them. However, with the advent of the cloud (for want of a better descriptor), and the advances in delivering Web Services and functionality via the cloud, how long do you think it will be before enough functionality is deliverable via the cloud? Or do you think, as I do, there will always be a place for desktop GIS? I am interested, as it is a developing situation and most of my clients are interested in knowing where we (GIS) are heading. For now, I still believe that the application or resolution of complex GI problems will necessitate the need for Desktop apps and that it will be some time before all the necessary functionality is deliverable via the web AND, I believe it is still not worth the extra development for companies in delivering or developing complex solutions for power GIS users. What say you?


title: Convert .dwg file into .dwf file using Teigha tool

text: I am using teigha tool to convert my dwg file into dwf file. Can anybody help to explain the procudure for the same. Your prompt reply will be greatly appreciated!


title: Sort by geometry using INHERITS?

text: hope everyone's doing well. I have a statewide point data set (parent table: 600K rows, 500 columns) I need to analyze at multiple scales based on several different spatial divisions. I'll also need breakdowns by a few key attributes. Here's a list of planned child table types with counts:   * mainRegions (4 child tables, a few extra columns for demographics)   * subRegions (11 child tables, a few extra columns for demographics)   * counties (67 child tables, a few extra columns for demographics)   * KeyAttributeOne (15 child tables, no additional columns)   * KeyAttributeTwo (48 child tables, no additional columns) Concerning spatial divisions, we aggregate and report exclusively on the region, sub-region and county level of the state (see image) but analyze the data at the census tract or block level. There are only two attributes in the parent table that can be used to group the records in any meaningful way. ![State of Florida Divided Three Ways](http://i.stack.imgur.com/K5d2y.jpg) It's my understanding that implementing an inheritance strategy will allow for easy aggregation by geometry, attribute, or a combination of both. I wonder, however, if the easy ad-hoc queries will be worth the work required up-front. **PostGIS In Action** (awesome resource) gives a good example of a linestring table sorted by state_fips but I haven't heard any firsthand accounts of using inheritance with geometries. -Rob


title: Saving the GUID in arcgis?

text: In my database of points with lat and long. After clipping the data to get only the ones in the specific state, my ID column that I had before disappeared. How can I clip my data and still keep the GUID that I had before? I am clipping directly from Microsoft SQL server. On the server the GUID column is type uniqueIdentifier.


title: How to create a color ramp in ArcGIS10 where the zero value is represented by white?

text: I have a simple problem which (to me) seems to be difficult to solve in ArcGIS. The map I plot is a GeoTIFF with emission data. I use a default algorithmic color ramp (rainbow) to do so. Now I would like to replace the lowest value (=0.0) with white instead of blue. The only solution I have come up with yet is to create a new Mulltipart color ramp with the first part being preset and the second algorithmic. This way, the color steps in the first part are much coarser than in the second one. Is there a way to reduce the number of colors in the preset color ramp to 1? Any help is really appreciated - thanks!


title: Fastest Way to Intersect Envelope with Multiple Feature Classes ArcObjects

text: I have the situation where I am doing a large number of intersections and I need to do them fast. The situation is I have a iMap with 1 or more layers and a series of envelopes, lots of envelopes. I have tried the following:                       long nFeatureCount;             iExtent->PutCoords(ulX, ulY, lrX, lrY);             // I created an array of all the layers feature classes at the outset             for(int i=0; i<nFeatureClassCount; i++)             {                 iSpatialFilter->putref_Geometry(iExtent);                      iFeatureClasses[i]->FeatureCount(iSpatialFilter, &nFeatureCount);                      if(nFeatureCount > 0)                 {                     return true; //we have an intersection                 }                  }      I also tried doing a select at the map level like so:                       long nFeatureCount;             iExtent->PutCoords(ulX, ulY, lrX, lrY);             ipMap->ClearSelection();             ipMap->SelectByShape(iExtent, NULL, VARIANT_TRUE);             ipMap->get_SelectionCount(&nFeatureCount);                  if(nFeatureCount == 0)             {                 return false;             }             else             {                 return true;             }      Both of these work and the `iMap` seems to be a bit faster when I have 5 or so layers. I was just thinking someone had encountered this and might have a faster way. I saw someone had the idea of throwing all the geometries in a `IGeometryBag` but wouldn't that negate the spatial indexing?


title: How to reproject WGS84 to a metric coordinate system with own reference point in QGIS?

text: I have GIS data in EPSG:4326 consisting of 2D points, linestrings and polygons that are all close (i.e. less then a couple of hundred of meters) to a certain WGS84 coordinate point. I would like to convert/reproject the geometry from a WGS84 to a 2D system in which the coordinates are north/south and east/west distance from the reference point in meters. In other words, the reference point should become the origin of a new euclidic SRS and the axis unit should be meters. How I would I do that in QGIS? I assumed I would be able to generate a custom SRS with my reference point as a center but I have no idea how to express this in QGIS? Was I right to assume that I would need to dig into PROJ.4 expressions for that? Or should I some other tool? I have the whole osgeo4w toolbox at my disposal.


title: How to make a custom properties dialog for a custom text background

text: I'm developing a custom text background in ArcObjects like in this example. I would like to know if it is possible to make a custom properties dialog to set up parameters of this text background. It would be another "Type" instead of "Line Callout" in the properties window: ![line callout parameters dialog](http://i.stack.imgur.com/CsmxU.png) The custom properties dialog would be packaged in the same .esriaddin than my custom text background. Does the ArcObjects java API allows this ? Does anyone know globally how to do this ? Thanks. I saw this example: [Sample: Simple dispersal renderer][3]. It's possible to make a custom property page for renderers and layers extending BaseCustomRendererPropertyPage and BaseCustomLayerPropertyPage. But which class could be extended to make a property page for a TextSymbol and a TextBackground ?


title: How to import shape file into SQL Server 2008 R2 using only SQL scripts?

text: I am trying to import some shape files into SQL Server 2008 R2. I found this useful tool called Shape2SQL however I am not sure how I can import shape file directly using only SQL scripts as this is important for build automation and deployment. So far I haven't found anything useful just yet. Would anyone be able to point some directions? Thanks.


title: Import Arcpy using PythonWin

text: I tried to install pythonwin26-212 from ArcGIS10.0 Desktop CD. pythonwin26-212 stops working as soon as I press "Enter" in the interactive window. (ps. I downloaded pythonwin26-212 from sourceforge too. It crashes the same way.) So, I downloaded pythonwin26-218 from sourceforge, and installed it. The pythonwin26-218 wouldn't crash but it couldn't recognize arcpy. I doubt the reason would be the unmatched directories. pythonwin is supposed to be in the same folder as numpy (?) if it was installed from the CD, but it is now in a different folder: C:\Python26\ArcGIS10.0\Lib\site-packages\numpy C:\Python26\Lib\site-packages\pythonwin After I move the pythonwin folder into C:\Python26\ArcGIS10.0\Lib\site- packages\ I click the pythonwin .exe file, it wouldn't open but gives the following message: Traceback (most recent call last): File "", line 1, in File "pywin\framework\startup.py", line 49, in **import** (moduleName) File "C:\Python26\ArcGIS10.0\Lib\site- packages\pythonwin\pywin\framework\intpyapp.py", line 3, in import win32con ImportError: No module named win32con Anyone knows how to fix this? Thanks! Thanks, A.


title: Mapinfo 11.5 - Legend Designer - problem with symbols

text: First of all I must say I love new legend designer. It is great to have finally better tool for positioning and even some changing the size of lines and region symbols, hooray. Big step for Mapinfo to be one day user friendly professional tool for map layouts. But.. When I use Legend Designer all content of legend window is automatically converted to bitmap for layout print/export to PDF. Is there some option to turn this off and keep text/vectors? I found out this way I cant use it at all with symbols. ![Problem of Legend Designer](http://i.stack.imgur.com/A8GCL.png)


title: How to calculate the scale of current Zoom Level with parameters like extent, resolution?

text: I got the following parameters:               Width and height of Map: 450px, 560px          units : 'meters'          maxExtent : new OpenLayers.Bounds(806677.9759,802420.9858,8876817.5647,853497.8186405064)          resolutions: [199.51887828322785, 99.75943914161392, 49.87971957080696,      24.93985978540348, 12.46992989270174, 6.23496494635087, 3.117482473175435,      1.5587412365877176, 0.7793706182938588, 0.3896853091469294, 0.1948426545734647,      0.09742132728673235, 0.048710663643366174, 0.024355331821683087, 0.012177665910841544,     0.006088832955420772, 0.003044416477710386, 0.001522208238855193]          numZoomLevels: 18          OpenLayers.DOTS_PER_INCH = 90.71428571428572;      Who inspire me how to calculate the scale of current Zoom Level?


title: Difference between DynamicMapServiceLayer and TiledMapServiceLayer

text: I'm a bit of a noob to Esri still so I have a few questions regarding these two types: What is the difference between DynamicMapServiceLayers and TiledMapServiceLayers? Do separate services need to be published to use each type for a single map? Are both valid types one can use in Esri's FLEX API?


title: Can someone give the absolute lowdown on the following list of spatial operators

text: Ok, I have to admit, I don't really delve too deep into the semantics of the spatial operators list below, I'm more of a user, I write software to drive the SQL most of the time and don't really think about it too much. However, I have a situation where a spatial operation I'm doing is going dog slow, on a table that has great indexes on and has never caused a problem before now. As a result, I'm trying to figure out which of the following is my best bet to use, that will give me the best performance, however, sometimes I find my searches don't return anything when I know they should, and sometimes return stuff when they shouldn't. SO I'm reaching out, to ask if some one can give me the actual definition, such as criteria to match, fully within enclosing geometry, or within outer ring etc etc of each of the following operations:               Equals(Geom1, Geom2):int     Disjoint(Geom1, Geom2):int     Touches(Geom1, Geom2):int     Within(Geom1, Geom2):int     Overlaps(Geom1, Geom2):int     Crosses(Geom1, Geom2):int     Intersects(Geom1, Geom2):int     Contains(Geom1, Geom2):int     Relate(Geom1, Geom2):int      If you can give a concrete example of when they will and will not match too, that would be great. For reference, I'm searching on a table of Linestrings, using a rectangular polygon, rotated to be pointing in the direction of travel (It's part of a traffic managment app) For further reference, here is a sample of the SQL the app generates to perform a search:               SELECT recordID,AsBinary(geometry) AS geometry,Distance(GeomFromText('POINT(-1.84101 54.85078)',4326), geometry) AS distanceFromGps FROM linegrid WHERE Intersects(GeomFromText('POLYGON ((-1.8413149820810311 54.850782468607292, -1.8409507853094111 54.850952257034713, -1.8408279009723911 54.850864894077496, -1.8411920982612455 54.850695105650068, -1.8413149820810311 54.850782468607292))'), geometry)      You can see in this example I'm using **intersects** trying the others as mentioned doesn't quite give the results I expect, hence why I'm asking the question. Please note, I'm not particularly bothered about actual faster/other ways of doing this at the moment, BUT if you do happen to spot a better way do feel free to shout up :-) For now all I'm mostly trying to understand is the differences between each of the operations. As for doing it faster/better I intend to open a different question for that at a later date. For more reference, the Spatial engine this is working against is Spatilite 2.3.0, and I can't upgrade it due to the devices it's running on not having a more up to date build for it (Windows CE) # Update 27-1-14 I'm going to choose mintix'es answer here as that provided a super easy to understand grid and description of all the various operations, which answered my question. I'm still faced with getting performance on the routines up, but I do now have a few strategies to try on that one, I may still seek advice here from the group with regards to that. @Vince deserves some credit too, as that white paper is definately the definitive answer, but it's also quite challenging to read (I had to read through it several times) , for anyone wanting the real I am on these operations The Clementini paper is the way to go, if you need an easy to understand overview while developing solutions however, the FME link should be open in another browser tab. Cheers Shawty


title: How to georeference downloaded map using ArcGIS for Desktop?

text: I want to georeference a downloaded map. ![enter image description here](http://i.stack.imgur.com/7PtsS.jpg) In the downloaded map there is no graticule but I know the projection which was used to make this map. I also know the the location (x, y coordinates) of four points. However, I cannot identify this location on the map. Can anyone help with advice on how to register such a map?


title: Mapserver WMS, features clipped and positioned

text: How can I prevent features being clipped at tile boundaries in an WMS layer as per the picture ![enter image description here](http://i.stack.imgur.com/wqTmr.png) The arrow and line segments are being clipped when viewed in openlayers as a WMS layer. The map file is as follows:-                     MAP                   NAME "WMSROUTETEST"                   CONFIG "PROJ_LIB" "/opt/company/servers/halo/company-global/site/web/conf/mapfiles"                   FONTSET "/opt/company/servers/halo/company-global/site/web/conf/resources/fonts/fonts.list"                   SIZE 400 300                   EXTENT 517531.000000071 139263.350999558 566460.000000404 185808.999999165                   STATUS ON                   UNITS METERS                   #RESOLUTION 72                   SIZE 580 400                   MAXSIZE 8192                   IMAGETYPE png                   OUTPUTFORMAT                           NAME png                           DRIVER AGG/PNG                           IMAGEMODE RGB                           MIMETYPE "image/png"                           EXTENSION "png"                           FORMATOPTION "INTERLACE=OFF"                           FORMATOPTION "QUANTIZE_FORCE=on"                           FORMATOPTION "QUANTIZE_COLORS=256"                           FORMATOPTION "COMPRESSION=9"                           TRANSPARENT ON                   END                   PROJECTION                           "init=epsg:27700"                   END                   WEB                       METADATA                           WMS_TITLE            "department Route"                           "wms_srs"            "EPSG:27700"                            "wms_enable_request" "*"                        END                       IMAGEPATH "/opt/company/servers/halo/company-global/site/web/public/tmp/"                       IMAGEURL "/tmp/"                   END                        SYMBOL                           NAME "arrow"                           TYPE vector                           FILLED true                           POINTS                               0 0.4                               1 0.4                               1 0                               2 0.8                               1 1.6                               1 1.2                               0 1.2                               0 0.4                           END # POINTS                           ANCHORPOINT 1 0.5                   END # SYMBOL                        SYMBOL                     NAME "square"                     TYPE vector                     FILLED true                     POINTS                       0 4                       4 4                       4 0                       0 0                     END # POINTS                     ANCHORPOINT 0 0.5                   END # SYMBOL                                  LAYER                           NAME "departmentrouteworkerline"                           GROUP "departmentrouteworker"                           METADATA                               WMS_TITLE "department Route"                           END                                TYPE LINE                           CONNECTIONTYPE postgis                           #EXTENT -6 50 6 60                           CONNECTION "host=localhost port=5432 dbname=company user=postgres password=GOAWAYNASTYHACKERS"                           STATUS OFF                           #Must validate an item if to be used in query                           VALIDATION                             'worker' '^[0-9]+$'                                DATA "wkb_geometry from (SELECT work_item_sequence as seq, count(work_item_sequence), (st_collect(route_line)) as wkb_geometry FROM ng_department.worker_%worker%_route_table GROUP by work_item_sequence ORDER by work_item_sequence ASC ) AS FOO using unique seq using srid=27700"                           EXTENT 517531.000000071 139263.350999558 566460.000000404 185808.999999165                           LABELITEM 'seq'                           #This prevents labesl being repeated across tiles                           PROCESSING 'LABEL_NO_CLIP=True'                           PROCESSING 'POLYLINE_NO_CLIP=True'                           OPACITY 80                                CLASS                                   #Route Line Style                                    STYLE                                       COLOR 255 255 255                                       ANTIALIAS true                                       WIDTH 8                                       LINECAP round                                    END                                     STYLE                                       COLOR 0 0 0                                       ANTIALIAS true                                       WIDTH 3                                       LINECAP round                                       PATTERN 40 10 END                                     END                                    #Route arrows                                   STYLE                                     #GEOMTRANSFORM "start"                                     GAP -300.0                                     COLOR 128 128 128                                     OUTLINECOLOR 0 0 0                                     SYMBOL "arrow"                                     SIZE 18.0                                     ANGLE AUTO                                   END # STYLE                                END                             END                END      Additionally I have some symbols like the one below that do not sit perfectly on the line, can these be nudged in some way? ![enter image description here](http://i.stack.imgur.com/7fOHk.png)


title: Unknown coordinate system warning when programmatically adding layer

text: I have another interesting question for you. I have data with unknown coordinate system. When I add this data using a standard tool(File->Add Data) in ArcMap, a warning dialog appears. It looks like: ![warning dialog](http://i.stack.imgur.com/SENWR.png) But this dialog didn't appears if I add data programmatically using arcobjects. Is there a way to show this dialog, if I add layers programmatically? P.S. I found DatumCheckerClass, but that's not what I need.


title: How combine different geometry types into the same feature class

text: I've searched the question in this blog. Maybe I don't use the right words but anyway my problem is not present in any post so I'll try to explain my question by opening a new one. I have a polygon shape (eg states). I calculated the centroids and the function saved them in an another file (Feature To Point function). Now I have 2 separated layers, point and polygon shape (centroids and polygons), but I need a unique shape. Functions Merge and Union (if I understand what I've read) can only be made between shape with the same feature class (point-point, line-line, polygon-polygon). I need a layer that contains, for example, the polygons of the states and their centroids Is there a way to get the results I need? Or a procedure to get the centroids of the polygons on the same shape? I hope I explained Sorry for the bad english!


title: How to arcpy.Describe() a layer when script is run from outside of ArcMap?

text: The following code will work if I run it from the ArcMap Python dialog, or from a GP tool, but **not from PyScripter IDE** :               import arcpy          mxd = arcpy.mapping.MapDocument(r"c:\temp\<my_mxd_name>.mxd")     layerList = arcpy.mapping.ListLayers(mxd)          for lyr in layerList:       print "Layer = " + lyr.name            #Get all fields in the layer       fList = arcpy.Describe(lyr).fields       for f in fList:         print "Field = "  + f.baseName      When I run the exact same code from PyScripter, it fails on `arcpy.Describe(lyr).fields` with the error `IOError: "<TOC_layer_name>" does not exist`, even though it does exist in the referenced MXD table of contents. When PyScripter fails it opens up the ESRI file `_base.py` to the `describe` function. I am assuming there must be something different about describing a layer in ArcMap vs calling it directly from Python, but cannot figure out what it could be. How do I properly use the arcpy.Describe() function on a layer when I am **not running the script directly from ArcMap**? I have figured out that if I use `arcpy.Describe(lyr.dataSource)` it functions properly, but only returns the fields from the source data layer. However, **since I'm trying to read fields from a joined table** (which does work in the original code above when running from AcrMap), I really need to see the fields from all joined tables, not just the data source.


title: Net Primary production (NPP) calculation

text: When I calculate Net primary production (NPP) for a local scale, what kind of software should I use? I always use ENVI, ERDAS, ARCGIS and IDRISI. Is it possible to use these programs for NPP calculation? Thanks


title: How to change default image in tree node with GeoExt

text: I'm using GeoExt to create a tree. In each node I want the corresponding image at left instead of the default image ![default image](http://i.stack.imgur.com/3XI5p.png) like this site. How can I do it? **Edit** I found this example but I don't know how can I add childrens **Edit 2** This is part of my code:               var bpn = new OpenLayers.Layer.WMS("bpn",          url,         {             LAYERS: [ 'Año 2004', 'Año 2005', 'Año 2006', 'Año 2007', 'Año 2009'],             format: "image/png",             transparent: "true",             projection: 'EPSG:4326'          },          {              buffer: 0,               displayOutsideMaxExtent: true,              isBaseLayer: false,              displayInLayerSwitcher: false,              yx: {'EPSG:4326' : true}           });          var treeConfig = [         {             nodeType: "gx_layer",             layer: bpn,             isLeaf: false,                  loader: {                 param: "LAYERS",                  },             expanded: false         }];          tree = new Ext.tree.TreePanel({         border: true,         region: "west",         title: "Entre Ríos",         width: 250,         split: true,         collapsible: true,         collapseMode: "mini",         autoScroll: true,              loader: new Ext.tree.TreeLoader({             applyLoader: true,             uiProviders: {                 "layernodeui": LayerNodeUI             }         }),         root: {             nodeType: "async",             children: treeConfig         },              rootVisible: false,         lines: false,     });      Each layer (node)`'Año 2004', 'Año 2005', 'Año 2006', 'Año 2007', 'Año 2009'` has a different image. What should I change in `treeConfig`?


title: How to tell if a device supports Assisted GPS, and if how to get the Assisted GPS cooridinates?

text: 1) Is Assisted GPS implemented at the hardware and OS level, or is it a hardware and secondary software thing. I've never used this so wondering what is required to 'access' this. I noticed a wiki article so is seems it is very device specific, but wondering if this is full and complete list. 2) How are the aGPS value accessed? Are they streamed through a diff port (v. true GPS)? Is that configurable? I tried to ask appropriate questions to understand more about this, but I'm kind of lacking context so my questions may need some edits from some who knows more about this. thanks df


title: Changing from one coordinate system to another without re doing

text: I'm new to this. Can I change from the Nigeria Mid Belt system to the Nigeria West Belt system easily, or do I need to recreate my project? Thank you!


title: IElement, extract Element from GroupElement nested in GroupElement

text: As per title. I am trying to extract Element from nested GroupElement. I can not find any code sample for that. Please refer to one or point me in to right direction, how to do that. It sounds simple but turned out to be nightmare. I cannot figure out how whole mechonism of Group in the Group works. Any advice?


title: ODP.NET invoked member is not supported in a dynamic assembly

text: I am using ODP.NET to get access to a geometry column in an Oracle database. My code (derived from Oracle's sample) often (but not always) crashes on the last line where the geometry column value is converted to the SdoGeometry type. Any ideas?               OracleConnection oraConnection = new OracleConnection("User Id=scott;Password=tiger;Data Source=oracle");     oraConnection.Open();     OracleCommand oraCmd = new OracleCommand();     oraCmd.Connection = oraConnection;     oraCmd.CommandText = "select geometrie from gw_pro where pro_id=3357";     OracleDataReader oraReader = oraCmd.ExecuteReader();     oraReader.Read();     Geo.SdoGeometry geom = (Geo.SdoGeometry)oraReader.GetValue(0);


title: What is ArcSDE API used for?

text: To create applications that manipulate ArcObjects there is the VBA and Java runtimes. To create application that process data with ArcToolbox tools, there is ArcPy for Python. Today, I discovered that there is a C and Java APIs for ArcSDE! What are ArcSDE APIs for? Do they provide exactly the same functionalities as ArcSDE command lines? Are there any other ArcGIS API?


title: What is happening when I register a "Feature Class" with Geodatabase?

text: I'm registering a table w/ native geometry with the geodatabase and it is taking a long, long time. Actually I'm not sure if ArcCatalog is frozen or crashing. The table does have several million features (points). The database is a SQL Server 2012 database and there is nothing happening on the database when I run SQL Profiler so whatever is going on, I think it's on the client. What is ArcCatalog (10.2) doing that is making it slow. Is there a way to make it go faster?


title: How can I programmatically open an OnPoint map to specific coordinates?

text: We have an OnPoint server at work and I'm wondering how to open the mapping service and having it zoom in to my specified latitude and longitude. The reason I ask is because we have an application at work that requires our users to manually type in the latitude and longitude into our system for searching. Our application does have this information so we'd like to offer the users a link. Here is an example URL that you could plug into Google Maps. Here is an example of a GIS system running OnPoint. I do see if their promotional PDF clearly states that you can perform searches via URL queries but I see _zero_ documentation.                Launch OnPoint from external applications by sending URL to OnPoint Website


title: What gestures make sense in a GIS application?

text: Have been testing Microsoft Beta APIs for kinect with the kinect Toolkit. To create an appliation that catch a users gesture is very easy. As a GIS developer it feels like it would integrate perfectly with ESRI ArcGIS for WPF, to build for example kiosk applications where a user can interact with a map. What gestures would make sense for a user to navigate a map?


title: How do I correctly perform a spatial join between lines that travel along each other in ArcGIS

text: I have 2 layers, one representing streets, and one representing pipes. The pipe layer has a street field that I want it to be automatically filled (I want the pipe to show the street that it's on) I tried performing a spatial join based on 'closest' but I encountered the following problem: many pipes never intersect the center of the streets they're on (they are usually paralel), and, at the spatial join, those respective pipes don't get the needed data from the street they go along but from the ones that are perpendicular (cause they intersect the pipes, so they have one point that is closest). I would like some advice. Thanks a lot.


title: d.legend script is not working

text: I'm using the d.legend script in grass but it is not working. I can only show/hide the legend using the 'add legend' option in the Map Display window in GRASS. I need to make d.legend script to work in order for me to be able to call GRASS in another program. The script should be working in order for me to be able to use it and not using the add legend option in the window.


title: Is there a name for deformation-based maps?

text: I have already seen maps using deformation to visualize information. For example, this map shows temporal distance between some French cities when traveling by train: ![](http://www.nf-ameublement- pro.com/sc/images/new/newsletter/septembre/carte_france.jpg) Is there a name for this kind of map? Also, is this type of map available in any GIS software?


title: How can I access the subdatasets of a grib file with gdalwarp?

text: How can I access the subdatasets of a grib file? I am working with NLDAS data and I want to process only the first band of the grib file. I know how to do it with HDF data but I couldn't handle it with grib files. Any hints? For example this one has 11 bands in it. NLDAS_FORA0125_H.A20070702.1700.002.grb I want to gdalwarp the first band of it which is temperature. This is the gdalinfo output.               Driver: GRIB/GRIdded Binary (.grb)     Files: NLDAS_FORA0125_H.A20070702.1700.002.grb     Size is 464, 224     Coordinate System is:     GEOGCS["Coordinate System imported from GRIB file",         DATUM["unknown",             SPHEROID["Sphere",6371200,0]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433]]     Origin = (-125.000500000000002,53.000500000000002)     Pixel Size = (0.125000000000000,-0.125000000000000)          Corner Coordinates:          Upper Left  (-125.0005000,  53.0005000) (125d 0' 1.80"W, 53d 0' 1.80"N)     Lower Left  (-125.0005000,  25.0005000) (125d 0' 1.80"W, 25d 0' 1.80"N)     Upper Right ( -67.0005000,  53.0005000) ( 67d 0' 1.80"W, 53d 0' 1.80"N)     Lower Right ( -67.0005000,  25.0005000) ( 67d 0' 1.80"W, 25d 0' 1.80"N)     Center      ( -96.0005000,  39.0005000) ( 96d 0' 1.80"W, 39d 0' 1.80"N)          Band 1 Block=464x1 Type=Float64, ColorInterp=Undefined            Description = 2[m] HTGL (Specified height level above ground)       NoData Value=9999       Metadata:         GRIB_COMMENT=Temperature [K]         GRIB_ELEMENT=TMP         GRIB_FORECAST_SECONDS=0 sec         GRIB_REF_TIME=  1183395600 sec UTC         GRIB_SHORT_NAME=2-HTGL         GRIB_UNIT=[K]         GRIB_VALID_TIME=  1183395600 sec UTC     ...


title: How to select a feature and get its fields by mouse click using Python Addin in ArcGIS?

text: I'm building a python addin tool in arcgis and I would like to click a feature, say a polygon on the map, and get the polygon's attributes,i.e. ID=3,NAME = "CITY3"...It's similar as the arcmap tool Identify without a pop- up window. I can only think of a way by getting the mouse click XY first, and do a select by location query to get the feature.I'm wondering is there a better way of doing this? Or is it possible to grab the information from Identify tool using python addins? Thanks a lot.


title: How to transform a DEM to a cartesian coordinate system?

text: I am using the script r.inund.fluv. I run the script with all the required parameters and it produced a lot or errors. I emailed the author of the script and he said that _"I think you are not using a cartesian coordinate system... ("The DTM resolution is 0:00:03 x 0:00:03 meters") and the code work with a cartesian grid. "_ I'm using a raster map from SRTM for the elevation. Its DTM resolution is 0:00:03 x 0:00:03 meters. As what is required for the script to work, the author said that I should convert it to a cartesian coordinate system. How is it done for a raster SRTM map?


title: gdal2tiles format of argument for the '--srcnodata' flag

text: When you run gdal2tiles, by default it sets pixels outside your map area to black-transparent. I would like to set them to white opaque. Looking at the documentation... http://www.gdal.org/gdal2tiles.html it looks like the setting `--srcnodata=NODATA` might do it, but it doesnt explain what the allowed values for the 'NODATA' argument are. I have made various attempts but they had no effect e.g... `gdal2tiles -r near --srcnodata=255 -z12-12 in.tiff outDirectory`


title: In a versioned geodatabase, what impact do delta tables and the state tree have on query performance?

text: We have a versioned arcsde geodatabase (arcgis 9.3.1 on oracle 10g) with a fairly complex data model that includes about 100 featureclasses and non spatial tables, a geometric network and many relationship classes. The data is edited daily by 5 or 6 arcmap users utilizing sde versioning. In addition versions are created by automatic services that interface with other business systems to perform edits in the geodatabase. Query performance degenerates noticeably during the course of the day, so we have implemented a nightly script to achieve a full compress. On occasions when a relatively large number of edits are performed, the system can become unusable until after a full compress. Its been suggested that oracle as configured can not come up with decent execution plans when confronted with these volatile delta tables. Is this a reasonable explanation? What approach should be taken to resolve it? **Update in response to comments**   * By end of day, the state tree is very linear, with only a little branching.    * We compress nightly (get a full compress by deleting all versions).   * Business tables are analyzed regularly.   * Delta tables are not analyzed. They are locked (Attempt to analyze returns error "ORA-20005 object statistics are locked"). Neither are the volatile tables in the sde schema - STATES, STATE_LINEAGES.


title: Trouble viewing MrSID in ArcView 3.3

text: I have been using 3.3 for some time and now I have Arcmap10 on the same computer. My problem is that I cannot pull up the Mr Sid files. 3.3 can't find them ,the id has the right info but no image of the mr sid file.What has happened.They worked together for some time and then one day nothing.error loading cannot make a theme. Thanks for any info


title: How to read Autodesk Spatial Data File (.SDF) in PHP

text: I am working on an application using **Autodesk MapGuide**. The application is being developed on PHP and hosted in IIS 7.5. I had to create point, line and polygon geometries on-the-fly. I found a working code sample and customized it for my logic. The sample writes the geometry in an **SDF** (Autodesk Spatial Data File) file. The code writes only three fields - `ID, Text and Geometry`. The ID field is **"AutoIncrement"** Now, I have to store the remaining data about the feature drawn on map (and stored in SDF) into a table in **MySQL**. For that, I have to fetch the `ID` of the last entry into the SDF file, which is not exposed at any point. Also, I have to select a row in MySQL table based on the feature I select on map/drawing (whose geometry is stored in SDF file) for updates. For that also, I have to again read/query the SDF. How can I achieve that using PHP? Thanks


title: Where are ComReleaser & WebObject ManageLiftime in ArcGIS 10.2.1

text: I used to use ComReleaser in desktop, and WebObject in ArcGIS Server apps to release com objects - e.g. cursors and the like. I am undertaking an upgrade from 9.3.1 to 10.2.1 of an app - and for the life of me I cannot find the ComReleaser and WebObject to help me manage com object lifetimes! Are they just not needed anymore in 10.2.1?


title: Using Qt Designer to Create a Form in QGIS that can search for a specific point in a layer

text: Good Day! I would like to ask if you can create a form in Qt Designer to create a form that can search for a point on the map using data that the user will provide such as the Name. I am currently creating a map in QGIS for my school project and i was thinking if i can create a form that can search the layer for a specific point and show the result in a dialog box using Qt Designer. I would appreciate all your help. Thank you


title: Create ESRI shape file with attributes from mySQL

text: I have some geometry which can be perceived as regions within a country and wish to create a shape file which contains all the regions along with some attributes (e.g. name, population, area code, etc ) I have the geometry from the regions in x,y WGS84 lat lon (not polygons so cannot be sure they are closed loops which is priority) currently stored in mySQL. I wish to be able to create a single shp (with associated dbf and shx files) to contain this data. Any pointers on command line tools and methods to do this would be gratefully received.


title: Dojo Slider to control tiled map services

text: am creating a web map to show Sea Level Rise in my County. I want to use a Slider for a user select the increment of Sea Level Change that wish to see. I have been able to create Programatically create the Slider using:               var vertSlider = new VerticalSlider({            minimum: 0,            maximum: 10,            discreteValues: 1,            //value: 0,            intermediateChanges: true,            style: "height: 200px; width: 15px;",                 onChange: function (value) {                dom.byId("sliderValue").value = value;            }            }, "vertSlider");      The maps layers are created with:               var layer1 = new ArcGISTiledMapServiceLayer("http://url");     var layer2 = new ArcGISTiledMapServiceLayer("http://url");     ... and so on till 10 layers are added.      I need now to be able to pass the value of the slider to draw one of my map layers. So when the slider is moved from position 1 to 2, layer 1 would be turned off and 2 would display. How would i do this using JS or the ESRI API? There are many examples using the TimeSlider, but those do not translate to the Dojo Slider.


title: Select feature or graphics for editing geometry

text: i research about selecting feature or graphics for editing geometry and i found geometry editor sample but in this sample only you can add new feature after adding, you can not select and edit. Now i can draw point, polygon, polyline on graphiclayer but when i select one graphic, i could select and add existing graphiclayer but when i change editmode, i could not edit. Emergency i need to select feature or graphics for editing geometry. Are there any sample or code or ideas ? Also in ios sample it exists.


title: Stop QGIS opening in a new workspace in Ubuntu

text: Very basic issue I imagine other Ubuntu users will have encountered: every time you open QGIS, it does so in a new workspace - note the arrow in the sidebar is empty not full, just after opening QGIS. I can see how this could be useful, if you want full screen editing and stuff, but it can also be really annoying! Anyone else had this issues? Sure there's a very simple solution, I've just not found it yet! ![qgis-load](http://i.stack.imgur.com/NIJKw.png)


title: Problem about converting raster to polyline and calculating total length

text: I have raster file of line segment. I need to convert it to line and get the distance of that line. How can i do using Python and ArcGIS 10? I have tried with GIS tools Convert Raster to Polyline but the line segment is not straight and it is divided into no of segments, isn't it possible to get single line segment. ![enter image description here](http://i.stack.imgur.com/SgoG4.jpg) I have attached herewith pictures. My problem is to find out the width of dam. For that I have drawn cross section at a regular interval along the stream. I have also shown how I will find out the width.. If you have any suggestion please provide me. And talking about my problem, after raster analysis of taking the difference of two raster. I have a raster segment shown in right side of picture. I tried to Convert Raster to Polyline but line is not straight and break into several segment. I need to have one line segment which is straight.


title: How do I adjust transparency on heatmap so that I can see background raster map? QGIS 1.8.0

text: I have created a raster heatmap for the first time, but I cant seem to adjust the transparency settings to let me see another background OS raster map beneath it. I have attempted to swap the positions of the maps around in the layers bar but this does not make a difference.


title: Creating pop-ups on point layers

text: i have a postgresql database which is refreshed all 10 minutes and i need to create graphs in pop-up with the data. When i click, a pop-up (or annotation ?) appears. How can i make it possible ? Thanking you in anticipation


title: Heatmap processed on the whole layer / dataset

text: How can I process the geoserver heatmap on my whole database rather than on the displayed part of it? What I want is to keep the same color for each place, everytime. Here is an illustration of what I get currently. ![enter image description here](http://i.stack.imgur.com/3DE48.png) ![enter image description here](http://i.stack.imgur.com/JAXdg.png)


title: Which application generates *.MAG, *.MOB and *.CAR files?

text: I have following type of files (I have tried FME too.I think they are DGPS generated file?)   1. ".MAG" , ".MOB" and ".CAR" extension files many more   2. ".SES" ,"ERR", "LOG" several  I have attached files in the link               Or see the Snapshots      ![Snap1](http://i.stack.imgur.com/VnNvR.png) ![Snap2](http://i.stack.imgur.com/1ZdGE.png) ![Snap3](http://i.stack.imgur.com/OTyuj.png) Now i want to know:   1. Which application generates these files?    2. How to read these files (esp ".MAG", ".MOB" and ".CAR") and    3. How to generate spatial data from these files (esp ".MAG", ".MOB" and ".CAR")


title: How to identify parameters sent from ArcMap to geoprocessing service on ArcGIS Server?

text: I have a question regarding the use of ESRI geoprocessing services and maybe some of you can give me a hint for the following "problem"... I have set up a geoprocessing service which should solve a Vehicle Routing Problem on a given network. If I call the service via Rest URL with the expected parameters the service does not calculate a route, but does also not print any error messages. The service itself works correct. I have tested the functionality by calling the service via ArcMap by entering the expected parameters via the GUI. The parameters are send to the service on the server, the results are calculated there and finally visualized in ArcMap. Thus I think the service works correct. Now the question: Is there any software that can show me how the request is build up / how the parameters are send from ArcMap to the server? My idea is to copy the syntax and to use it in the Rest URL. I have already tried Charles (http://www.charlesproxy.com/) but it doesn´t show any paramters. Another idea was to look at the log files on the server, but also here I cannot identify any parameters. Thanks!


title: Finding points a certain distance from each other

text: I would like to find the maximum amount of points within a specified domain that are a certain distance from each other. For example, how many points within an area are 1000ft from each other? I envision each point with a radius of the desired distance, but I don't know if I must start with a point or what.


title: Where can I find Road shapefile/e00/etc for the entire world?

text: I am hoping that there is a single source for this data, ie. 1 geodatabase. If not, can someone point me towards the sources that would make up the entire world's road transportation network? I have come across references to the ESRI Digital Chart of the World, but I can't find any place to download it. Ideas?


title: How can I transform geometry data to bounds?

text: My SRID is 3137 and I want to transform geometry data(e.q : "0101000000904160258CC91841C8073D6723545141") to bounds. I need to bounds coordinates for use in maxExtend() function. Help me please!


title: QGis overlap labels

text: I have a csv file imported into QGis 1.8 and this layer is receiving real time data.Into this attribute table layer there is a field which values i use to label it but overlapping all values when showing. Is there any way to label only the last record when the attribute table is updating?. Thank you very much


title: Import TomTom OV2 into Google Maps as My Place Stars

text: Is there anyway to import my TomTom POVs (saved in an OV2 file) into my google maps account as favorites? i.e. My Places > Starred. I appreciate this may involve converting the OV2 into something else and possibly using google earth or other software. Thanks


title: Legends for WMS layer on Geoserver

text: Is there any way to get a representation of a WMS legend in a form that isn't an image? eg XML Using an image for the website in our system will mean the text is not sharp, it doesn't resize with the web site, and the format of it is limited to what WMS provides.


title: PyQGIS standalone gdal_ECW_JP2ECW.dll error

text: This is my code in a batch file My error messages are gdal_ECW_JP2ECW.dll error unable to open EPSG support file gcs.csv Thanks in advance               REM Change OSGEO4W_ROOT to point to the base install folder     SET OSGEO4W_ROOT=C:\OSGeo4W     SET QGISNAME=qgis     SET QGIS=%OSGEO4W_ROOT%\apps\%QGISNAME%     set QGIS_PREFIX_PATH=%QGIS%          REM Python Setup     REM set PATH=%OSGEO4W_ROOT%\bin;%QGIS%\bin;%PATH%     set PATH=%OSGEO4W_ROOT%\bin;%QGIS%\bin\gdalplugins;%QGIS%\bin;%PATH%     SET PYTHONHOME=%OSGEO4W_ROOT%\apps\Python27     set PYTHONPATH=%QGIS%\python;%PYTHONPATH%          python StandAlone_Field.py


title: Launching Open Layer

text: ![Road](http://i.stack.imgur.com/sOcQI.jpg)Bound Image....I currently have a few shapefiles of a place with a background map of that place. I have imported both shapefiles into postgresql and geoserver. After loading the layers onto geoserver I am being able to open the WMS layer in the OpenLayers preview of GeoServer. Is there any way I can add extra data to these shapefile and view those data in the WMS layer? Also how do I launch the open layer in localhost?


title: Restrict pan and zoom in ArcGIS Javascript API?

text: I'm creating a web map of a university campus. I have the map functioning, but I don't want to be able to pan off of the basemap, into whitespace. Do you know of a function that will auto-pan the extent back to the basemap?


title: How do you Turn on the Grass plugin in Quantum GIS?

text: I downloaded and installed Quantum GIS 1.5.0 with Grass using the OSGeo4W installer on my XP machine. Unfortunately, when I go to Plugins>Manage Plugins I cannot find the Grass plugin. When I browse to OSGeo4W\apps\qgis\plugins I can see grassplugin.dll, grassrasterprovider.dll and grassprovider.dll. So, I would appreciate it if someone could tell me how to make QGIS recognize the Grass install. Jeremy


title: How to find GIS consultancy/project work (rather than jobs)

text: The post GIS jobs online listed resources for finding a GIS **job** (as in a full-time or contract position). What resources are available for finding GIS **projects** or **consultancy work**? eg:   * a small project to integrate a map into an existing site;   * incorporate an address-search into a workflow;   * migrate legacy GIS code to a new format; These are one-off tasks which require a GIS professional, but which don't entail a full-time position. How do you find the right people for the task, or find the task if you're the right person? All the online resources I've seen are aimed at "jobs".


title: Is there a Web app for building GeoExt tree configurations?

text: I have a GeoExt application and looking for approaches for building tree configuration (in JSON format). Now tree configuration is made manually and looks like:                   var treeConfig = new OpenLayers.Format.JSON().write([         {             text: "Layers",             expanded: true,             icon: "./images/Layers.png",             children: [                 {                     text: "Nature",                     icon: "./images/globe.png",                     checked: false,                     group: true,                     children: [{                         nodeType: "gx_layer",                         layer: "Lakes",                         icon: "./images/lake.png"                     },                     {                         nodeType: "gx_layer",                         layer: "Bogs",                         icon: "./images/bog.png"                     }]                 }                 ...      But I would like to have web-application which I can use for building such configuration. All layers are published as WFS layers and it is possible to get the list of layers using GetCapabilities request, for example. Please give me any advice or may be known approaches how I can solve my problem.


title: Layer intersection

text: I'm trying to implement an intersection with ArcGIS 9.3.1, flex and rest api. I have two layers, I query the first one that returns a Graphic object, with this Geometry I want to query a second layer with an intersection. The problem is that the second query returns all the features even outside the intersection geometry. There's an example outside over the net? :)


title: Does anyone know when will ESRI re-introduce "answered thread, MVP, etc" features in the "new" forums?

text: There are a few simple features that made the "old" ESRI forums very useful. 1\. Marking a post as answered 2\. MVP point structure which Does anyone know when will ESRI re-introduce these features in the "new" forums?


title: Can I save MapInfo thematic*.thm files to a database?

text: I have an MS-SQl database. Is it possible to save mapinfo *.thm files to a database as I would save mapinfo *.tab files?


title: Piechart of a text column from attribute table in qgis

text: I have this problem: I would like to create a pie chart from a column of a attribute table, and I would like to see this pie chart above the map. The column contains names not numbers (I think fields aren't in a numeric format... on the other hand they're names!) I work with marine species distribution data and I built a database of records of many many species... In the specific, I have a column, called 'species', where there are many records (names) of some marine species. Several species may have many records, other species may have only a few numbers of records, so my objective is to graphically see the distribuition of records among the species. Please see the attachments and thank you for your help img1 http://postimg.org/image/rn56c8l4z/ img2 http://postimg.org/image/e6918ynmj/


title: Is there a way for me to load a KML files locally in Openlayers

text: Recently I posted this question KML file seems to load in Openlayers without error but doesn't show vectors because I was having trouble getting KML vector data to show up. After searching around a bit it seems there is an issue with loading KML files locally. I am developing a stand alone Application that will run on the desktop or on IOS. I am running openlayers via a WebKit view. I am loading my openlayers page from the local file system witch is part of the applications files when it's installed. My guess is there is some sort of issue loading local files. Is there a way around this? I would even consider converting all my KML files if needed. Ideally I would like to be able to point to either a local file or a remote kml file and load it. Here is the code I am using. It's from the Openlayers kml example that's online. The only difference is I am trying to load locally without running the html page on any sort of server.               var map = new OpenLayers.Map({         div: "map",         layers: [             new OpenLayers.Layer.WMS(                 "WMS", "http://vmap0.tiles.osgeo.org/wms/vmap0",                 {layers: "basic"}             ),             new OpenLayers.Layer.Vector("KML", {                 strategies: [new OpenLayers.Strategy.Fixed()],                 protocol: new OpenLayers.Protocol.HTTP({                     url: "kml-files/lines.kml",                     format: new OpenLayers.Format.KML({                         extractStyles: true,                          extractAttributes: true,                         maxDepth: 2                     })                 })             })         ],         center: new OpenLayers.LonLat(-112.169, 36.099),         zoom: 11     });


title: How to remove legends in cartoDB?

text: I am able to add and remove cartoDB layers over my google maps but after removing layer the legend stays . how to remove that legends ? here is code that i used   document.getElementById('school').onclick = function() { var layerUrl = 'http://astrosnehakul123.cartodb.com/api/v2/viz/4c13821a-e192-11e3-8f8f-0e73339ffa50/viz.json'; if (school.checked) { // add cartodb layer cartodb.createLayer(map, layerUrl) .addTo(map,0) .on('done', function(layer) { layer1 = layer; var sublayer = layer.getSubLayer(0); sublayer.on( function(e, pos, latlng, data) { cartodb.log.log(e, pos, latlng, data);                                   });                              sublayer.on('error', function(err) {                           cartodb.log.log('error: ' + err);                         });                            })                       .on('error', function() {                         cartodb.log.log("some error occurred");                       });                }             else{                      layer1.remove();                 layer1.clear();             }         }


title: Create a raster map from more than one shape

text: I've some shapefile loaded in Arcmap (steets, buildings...) with an appropriate symbology ad I need to convert them to unique raster, to obtain a map (like topographic map) with a specified resolution (cellsize) ad scale. There's a method, or plugin or program to do that? Thanks!


title: Reversing a "combine" function on a raster

text: I have a raster that I think was "combined" using an arcmap function. I am wondering if there is a way to extract just one raster layer from the combined dataset in order to use it separately, either in R or Arcmap. If you want to see the data it is available here It takes this form in R:               class       : RasterLayer           dimensions  : 10995, 10349, 113787255  (nrow, ncol, ncell)          resolution  : 100, 100  (x, y)          extent      : -446880, 588020, -612570, 486930  (xmin, xmax, ymin, ymax)          coord. ref. : +proj=aea +lat_1=34 +lat_2=40.5 +lat_0=0 +lon_0=-120 +x_0=0 +y_0=-4000000 +datum=NAD27 +units=m +no_defs +ellps=clrk66 +nadgrids=@conus,@alaska,@ntv2_0.gsb,@ntv1_can.dat           data source : G:\Faculty\Mann\Share\Fire\FRAP_FireThreat\FireThreat.tif           names       : FireThreat           values      : 0, 255  (min, max)          attributes          fields : ID Rowid   COUNT FUEL_RANK FROTCLASS URBAN THREAT THREAT2PEOPLE             min :  1     0       1        -1         0     0     -1            -1             max : 97    96 5360409         3         3     1      4           424


title: How can I permanently reorder fields in a feature class programmatically with ArcObjects, VB.Net?

text: Would someone be able to share the logistics on how to permanently programmatically reorder fields in a feature class with ArcObjects, VB.Net. I want to be able to do this in the geodatabase not in the IOrderedLayerFields Interface of the Carto assembly.


title: Open Street Map upload session error

text: While running Open street map file retriever I get the error - 174434708 / 7698676919 (7524242211 left(s)) Error : error retrieving or decompres file : java.net.SocketException: Connection reset Is there a way to re run the last two step without doing it all over again? The import has already been done and has taken 344 hours 7 minuts 48 seconds


title: Clustering not working on features in OpenLayers

text: I've built a mapping solution using OpenLayers that works well except I can't get the points to cluster. I've looked at the examples on the site, and I've even implemented those examples in my environment to make sure they work. When I try to do it with my map, however, all of my points show up, but they're not clustered.   Here is my code where I create the Vector Layer:               var layerVector = new OpenLayers.Layer.Vector("Regions",{                      strategies: [                                  new OpenLayers.Strategy.Cluster()                         ]                 });      I add all of the points to an array of features:               var feature_point = new OpenLayers.Feature.Vector(                                 this, {                                     "rName" : rObj.tName                                 }, {                                     'fillColor' : rObj.rColor,                                     'strokeWidth' : borderWidth,                                     'strokeColor' : rObj.rColor,                                     'fillOpacity' : regOp,                                     'pointRadius' : regPercent                                 });                                 fArray.push(feature_point);      Add the features:               layerVector.addFeatures(fArray);      And add the layer to the map:               window[mapName].addLayer(layerVector);      Seems pretty simple, but it seems like I'm missing something. In the example on the OpenLayers site, they add a styleMap to the layer, but I tried this and it doesn't do anything for me. I'm honestly not sure what it even means:                var style = new OpenLayers.Style({                         pointRadius: "${radius}",                         fillColor: "#ffcc66",                         fillOpacity: 0.8,                         strokeColor: "#cc6633",                         strokeWidth: "${width}",                         strokeOpacity: 0.8                     }, {                         context: {                             width: function(feature) {                                 return (feature.cluster) ? 2 : 1;                             },                             radius: function(feature) {                                 var pix = 2;                                 if(feature.cluster) {                                     pix = Math.min(feature.attributes.count, 7) + 2;                                 }                                 return pix;                             }                         }                     });                          strategy = new OpenLayers.Strategy.Cluster();                          clusters = new OpenLayers.Layer.Vector("Clusters", {                         strategies: [strategy],                         styleMap: new OpenLayers.StyleMap({                             "default": style,                             "select": {                                 fillColor: "#8aeeef",                                 strokeColor: "#32a8a9"                             }                         })                     });      Any help would be immensely appreciated.