title: PostGIS Query: All rows within 1 hour of each other which intersect

text: I'm using postgres (postgis) and i want to perform a query that returns all rows which geometries intersect and have a timestamp within 1 hour from each other. **What would the SQL statement look like?** table looks like this:               id   { integer }         name { character varying }     bounds { geometry }       timestamp { timestamp with time zone }          // purposely left out bounds in the table example below (not necessary)          id   name       timestamp     ---+----------+--------------------------------     1    one        "2010-09-24 21:10:39.515+00"     2    two        "2010-09-16 09:21:09.362+00"     3    three      "2010-07-08 00:00:46.549+00"     ...      **EDIT #1** Not including the intersect matching here is an example: for every row give me any other rows that are within 1 hour:               select * from myTable order by t           id |               t                    ----+-------------------------------       9 | 2011-07-15 18:20:20.05+02      10 | 2011-07-15 19:05:00.05+02      11 | 2011-07-15 19:40:20.05+02      13 | 2011-07-15 20:31:01.05+02      14 | 2011-07-15 20:35:11.05+02     (5 rows)          result of needed query:           id |  matchid |  origTime                  | matchTime     ----+----------+----------------------------+------------------------------       9 |   10     | 2011-07-15 18:20:20.05+02  | 2011-07-15 19:05:00.05+02      10 |    9     | 2011-07-15 19:05:00.05+02  | 2011-07-15 18:20:20.05+02      10 |   11     | 2011-07-15 19:05:00.05+02  | 2011-07-15 19:40:20.05+02      11 |   10     | 2011-07-15 19:40:20.05+02  | 2011-07-15 19:05:00.05+02      11 |   13     | 2011-07-15 19:40:20.05+02  | 2011-07-15 20:31:01.05+02      11 |   14     | 2011-07-15 19:40:20.05+02  | 2011-07-15 20:35:11.05+02      13 |   11     | 2011-07-15 20:31:01.05+02  | 2011-07-15 19:40:20.05+02      13 |   14     | 2011-07-15 20:31:01.05+02  | 2011-07-15 20:35:11.05+02      14 |   11     | 2011-07-15 20:35:11.05+02  | 2011-07-15 19:40:20.05+02      14 |   13     | 2011-07-15 20:35:11.05+02  | 2011-07-15 20:31:01.05+02     (10 rows)


title: How to manipulate annotation text outline color in ArcEngine?

text: I can get the color of my annotation text like this:               IColor TextColor = lpLabelEngine.Symbol.Color;      How can I get / set the outline color? Something like this, I want a black letter with a white outline, and the white outline is 3 pixels wide.


title: Convert to DXF not displaying correctly in CAD

text: I am trying to export to DXF in QGIS 2.2.0, when I import in (AllyCAD 3.6) the shape of the polygons do not match my reference shape. I presume it has something to do with projections, but I have tried various projections and cannot seem to get this right. Any help would be appreciated.


title: Does ArcGIS 10 have a "Match Properties" tool like AutoCAD

text: My question is short-and-sweet.. I have text, pointers & lines to be added to a series of maps. ** **Creating spatial data via an edit session is not an option** ** This project is HUGE, and there will be hundreds of maps for each facet of the project. Consistency is important. My question .... is there a tool in ArcView with the same functionality as [Match Properties][1]? It would save me SO MUCH TIME!! This "one touch" functionality would be far superior to changing properties via the multi-stage process in ArcView (especially for rotation angles). I appreciate an feedback. Take care all


title: Calculating volume with Cut Fill tool

text: This might be a very simple question - but it is quite vital and I want to make sure I am thinking this through right. I have used the cut/fill tool to calculate the difference in volume of a volcano cone. I have an attribute table, is it just a case of adding up all the -ve values to find the overall volume? (negative because the volcano has inflated i.e. I should be gaining volume). * * * So, This is an example of my attribute table. However, I think that I actually need both the negative and the positive values. So, are you saying that I would have to do (2*-110676.30...)+(1*17935.13...)+(1*-55080.16...)+etc? Thanks for your help. ![Volume attribute table](http://i.stack.imgur.com/q5DjA.png)


title: How to convert an XML to Shapefile using FME?

text: I am looking to convert XML to shapfiles using FME. The XML files are user generated while they draw points and lines (mainly) of interest on a web map. We generate the XML in the backgroung containing the spatial information of user clicks. I would like to ask whether we need to follow a specific format during XML generation so as to proceed to converting them to shapefiles using FME. Thank you **Here is an XML snippet**               <?xml version="1.0" encoding="UTF-8"?>     <Configuration>      <Note>       <NoteTypes>        <Type Desc="Point:a" Id="0"/>        <Type Desc="39.66195294811828" Id="1"/>        <Type Desc="20.84663271903992" Id="2"/>       </NoteTypes>      </Note>      <Note>       <NoteTypes>        <Type Desc="Point:b" Id="0"/>        <Type Desc="39.6634561263267" Id="1"/>        <Type Desc="20.84773778915405" Id="2"/>       </NoteTypes>      </Note>      <Note>       <NoteTypes>        <Type Desc="Point:c" Id="0"/>        <Type Desc="39.66210161555193" Id="1"/>        <Type Desc="20.8490788936615" Id="2"/>       </NoteTypes>      </Note>      <Note>       <NoteTypes>        <Type Desc="Line:l1" Id="0"/>        <Type Desc="39.66206857837209" Id="1"/>        <Type Desc="20.84651470184326" Id="2"/>        <Type Desc="39.66291928072009" Id="3"/>        <Type Desc="20.84694385528564" Id="4"/>        <Type Desc="39.66319183362618" Id="5"/>        <Type Desc="20.84708333015442" Id="6"/>        <Type Desc="39.66352219934388" Id="7"/>        <Type Desc="20.84763050079346" Id="8"/>       </NoteTypes>      </Note>      <Note>       <NoteTypes>        <Type Desc="Line:l2" Id="0"/>        <Type Desc="39.66354697670904" Id="1"/>        <Type Desc="20.84759831428528" Id="2"/>        <Type Desc="39.66295231749322" Id="3"/>        <Type Desc="20.8476734161377" Id="4"/>        <Type Desc="39.66248980123196" Id="5"/>        <Type Desc="20.84869265556335" Id="6"/>        <Type Desc="39.66220072699667" Id="7"/>        <Type Desc="20.84896087646484" Id="8"/>       </NoteTypes>      </Note>      <Note>       <NoteTypes>        <Type Desc="Line:l3" Id="0"/>        <Type Desc="39.66203554117646" Id="1"/>        <Type Desc="20.84653615951538" Id="2"/>        <Type Desc="39.66141609083402" Id="3"/>        <Type Desc="20.84785580635071" Id="4"/>        <Type Desc="39.66219246771504" Id="5"/>        <Type Desc="20.84893941879272" Id="6"/>       </NoteTypes>      </Note>     </Configuration>


title: "Upgrade Metadata" prompt in ArcCatalog

text: A few months we began seeing a pop up in ArcCatalog when hitting the "Description" tab asking if you want to upgrade metadata. Clicking yes or no results in the existing metadata being deleted. Clicking the x to close the pop up is the only way to keep the existing metadata intact. Is this a known bug or am I missing something? Searched around with Google but didn't find much. ![Upgrade Metadata prompt in ArcCatalog](http://i.stack.imgur.com/DpEeI.png)


title: Why does collapse dualline to centreline only create an outline of the dataset?

text: I'm trying to collapse all divided roads in one centreline and to integrate any other spurs or errors in the dataset to create one set of centrelines. I have used Merge Divided Roads but it leaves a lot to be desired (most of the roads I wanted collapsed did not merge). Here is the dataset: ![enter image description here](http://i.stack.imgur.com/853RY.jpg) And here is what happens when it is collapsed: ![enter image description here](http://i.stack.imgur.com/h2ACN.jpg) I have also repaired geometry, reprojected and made it into a singlepart feature but nothing seems to work. I have tried it on multiple computers.


title: how do you convert graphics to features and then associate those features with an existing layer?

text: Say for example, that you have water ties. These ties connect from a home to the main water line along the street. You want to draw these ties in ArcMap and give them attribute values regarding how long the ties are. However there is already an existing layer and many of the ties have already been drawn and connected to this layer. You are told to finish drawing the remaining lines and input the remaining measurements for those lines into the already established layer's attribute table. After drawing a line from the Polygon (home) to the waterline, how do you relate it to the existing layer so that you can enter its length information?


title: z-Coordinate and interpolation plugin in QGis

text: I've a csv file that includes weather stations locations (x,y,z) and precipitation values. I loaded the file by using add delimited text and then converted it into a point shapefile WS84. By using the interpolation plugin I managed to create the different rainfall zones using pseudocolor and later I extracted the contour. However, as I have z values for elevation I'd like to take this into consideration. The interpolation plugin includes the option to select z-Coordinate to load elevation points but I cannot select it. When I open the plugin without any layer loaded the check box works fine, but when the layer is present is not possible to select the check box. The values of elevation column (z) are in mm. Am I missing anything ? Qgis 1.8.0, interpolation plugin loaded from the repository, running Windows Vista. Addition: I believe the answer is below. z data must be within the shapefile not as a feature. Thanks Bernd V! "Also note that you can use the z-coordinate for the interpolation. This is the elevation of your data that is stored within the geometry itself and not in an attribute column. This is only possible if your GIS data is in 3 dimensions." in gistutor.com/quantum-gis/20-intermediate-quantum-gis- tutorials/… . So it seems you need to generate "real" 3D-data first, but this is a blank spot for me, sorry. – Bernd V


title: Count the number of corners of a polygon

text: How can I automatically count the number of corners that has a polygon uisng VBA and ArcObjects? Thanks


title: Measure the angle of each corner of a polygon

text: How can I automatically measure the angle for each corner of a polygon (in a vector GIS) using VBA and ArcObjets? Thanks


title: How to get multiple shortest paths with PgRouting in one query?

text: I want to run shortest path algorithm on multiple source&target pairs at once and get a result as a table and process it then. How do I do this? This query does not work:               SELECT a.source, a.target, paths.*      FROM all_to_all a, shortest_path_astar('...', a.source, a.target, false, false) paths;          ERROR:  function expression in FROM cannot refer to other relations of same query level      (btw, all_to_all does not literally mean all to all, :) it's just a number of random pairs) This does not work either:               SELECT *      FROM all_to_all a, (        SELECT * FROM shortest_path_astar('...', a.source, a.target, false, false) yyy     ) AS t2;


title: Adding a latitude/longitude graticule in CartoDB (e.g. WMS)

text: I'd like to add a basemap that has a graticule in latitude/ longitude to the Map View of a CartoDB project. This could be in the WMS format for example, but I haven't been able to find one that works. For my ArcMAP work, I first tried this NOAA WMS link, but got an "unspecified error" message. When I tried it in CartoBD, I get options for 1, 5 and 10 degree, but these end up with a blank background. http://maps.ngdc.noaa.gov/arcgis/services/graticule/MapServer/WMSServer?request=GetCapabilities&service=WMS I can email NOAA about the error, but ideally, I would prefer a satellite or terrain basemap that has lat/longs included in the map. So, I thought I'd ask the forum first. For example, I have an ArcGIS server link from NOAA that works great in ArcMAP: GIS Servers\arcgis on maps.ngdc.noaa.gov\arctic with many layers of interest.


title: torque cartodb heatmap at higer zoom levels

text: I am trying to use cartodb torque to display overtime location data using Cartodb Torque. Vizualization works great at zoom level 9, but when i zoom in further to level 13 instead of the intensity i see scattered points (rightly so because i dont have points to group at certain locations for that zoom level). same behavior can be seen in the demo as well: http://blog.cartodb.com/post/66687861735/torque-is-live-try-it-on-your- cartodb-maps-today My question is it possible to have density map displayed even at higher zoom levels when using Torque? Thanks


title: Creating curves using the Parcel Editor tool in ArcGIS 10.1

text: We are migrating parcel editing environments from AutoCAD to ArcGIS and I'm having trouble creating curves within a Parcel Fabric that are not tangent to lines. For example, what would be some methods for creating the curve (a=9.618, r=190.000) for this image. ![enter image description here](http://i.stack.imgur.com/CaHYb.png) I have yet to find good documentation for performing this task. Anyone know of good sources?


title: Cannot import gdal_array

text: I attempted to use GDAL without installing due to not having rights to install on this computer. I managed to get several gdal modules working, but when I attempted to use ReadAsArray() I received the import error               ImportError: No module named _gdal_array      I looked through the python folder and found that gdal and gdalconst both had pyd files and compiled files, but gdal_array does not. When I ran VersionInfo() I got the value 1910 for gdal if that helps. I also have python 2.6 and numpy Does anyone know a way I could fix this problem without using a gdal installer?


title: Export my map in various formats

text: I want to implement a sort of content negotiation so as to export my map in formats like pdf jpg etc. Is there any way to export my current map view (with all the enabled options so far) in various formats? With the term current map view i mean "what the end user sees at a specific moment". My map is generated by combining map Services from different web servers in the Web. As a result an option such as combining my services in arcMap, publish my map as a MapService and use the export capability provided cannot help: Ex. http://sampleserver6.arcgisonline.com/arcgis/rest/services/USA/MapServer/export?bbox=-241.20512466608835,-27.80850864596666,51.688012394341115,67.48206676752923 I want to get the data live directly from the servers. In fact would like to use a function like myMap.export("format") where myMap is the current map object shown by my browser.


title: Can't add SOC machine

text: I am trying to add a SOC machine to SOM server (distributed setup). From the server, I can see the service running. I have all the accounts properly setup. I may be still missing a super critical settings somewhere. This is the error that I am getting: > Machine is not a valid server container. Machine is not found. My environment is Windows Server 2008 R2 (both SOM and SOC machines), ArcGIS Server 10. I am curious if some had any idea what else to try. Is there a way to check if SOC is properly configured without using ArcCatalog or ArcGIS Server Manager. Thanks Jay


title: How to realign misplaced GeoTiffs?

text: I have bought a few (30) GeoTiff maps from my country's Geo-department - South Africa. They all have a date in the filename and the newest one was created in 2004, I see one is dated 1977 - that old! Anyway, they came with a `.tfw` file each. When I simply copy these into Google Earth, or use MapTiler to create a default-settings WSG84 tile set, the map is unaligned. Meaning that they are really close to their correct positions but a bit off center and the deviation from supposed-position differs everywhere on the image/map. Either the maps are wrongly geo-tagged or I use the wrong formats. Please note: I have no clue or experience with maps and their types (WSG84/UTM etc) I need to convert these maps into tiles for my iPhone app. I have also tried using maps of USA and they worked and were outputted perfectly. All I can tell about these maps which may help pro's is:   * It is a big white block with the map inside and the key/legend on the sides/bottom.   * The map inside has an arced top and bottom and sides/longitude-lines are skew. I would like to think it is the same as this image: (Only a small part thereof of coarse!) ![Random earth image](http://i.stack.imgur.com/D6l53.png) I have tried to realign/resize the map as an image-overlay on Google Earth, but then had no idea of how to output it or get the needed (and corrected) data from it. If anybody know how to use these, convert them to tiles or tips/tutorials that could help me fix this and get the maps straight, please tell! _If possible I would prefer to do it in MapTiler or at least on a Windows OS machine_ Thank you!


title: How can you find ArcGIS version programatically?

text: Is there a way using ArcObjects.net to find out what version of ArcGIS is installed on a machine (i.e. 9.3., 10.0, 10.1)?


title: What means the gdal_translate message: "Too many command options"?

text: I have a little question, regarding the gdal_translate-utility: I want to extract raster-files from database-stored gridfiles using the gdal_translate: > gdal_translate -ot Float32 -of SAGA PG: "host='localhost' port='5432' > dbname='$dbname' user='$dbowner' schema='grids' table='testrun' mode='2'" > /output/... But instead of delivering output it shows me "Too many command options." Then it lists usage and supported output formats. What does this message mean? I was searching the net already but found nothing which could be helpful. Thanks for an answer in advance, knutella


title: Python-arcpy : How make feature list from gdb and value list from field and use them in gaLayer function?

text: I'm not an expert in python, however I have to write a script to do geostatistical analyses with the GACreateLayer function from Arcgis. I have several data (>2000) and I want to use list of features and list of field values to perform cokriging. I have written a script, the first part return list of features in my gdb and the list of values in the field that I have chosen (it works and return field values and feature names) and the second part do cokring (it works for one feature and its field) However when I want to use both lists (inputDSet, end of script) in order to create a GA Layer it doesn't work.               File "C:\Program Files\ArcGIS\Desktop10.0\ArcPy\arcpy\ga.py", line 921, in GACreateGeostatisticalLayer     raise e     ExecuteError: Failed to execute. Parameters are not valid.     ERROR 045001: Input dataset(s) error. Table of inputs is not complete.     Failed to execute (GACreateGeostatisticalLayer)      I don't know how correctly use lists? May I transform them or something else? Here, is my code:               import sys, string, os, arcgisscripting     # Create the Geoprocessor object     gp = arcgisscripting.create()          def listFcsInGDB():         # list all Feature Classes in a geodatabase          gp.workspace = "C:\\INTERP\\PJ_class.gdb"              fcs = []         for fds in gp.ListDatasets('','feature') + ['']:             for fc in gp.ListFeatureClasses('','',fds):                 #yield os.path.join(fds, fc)                 fcs.append(os.path.join(fds, fc))         return fcs          fcs = listFcsInGDB()     for fc in fcs:         print fc          # List of value in fields     input_dataset =  fc        Atts = 'PJ_RacPJ'    #field with rainfall values                 rows = gp.searchcursor(fc)     row = rows.next()          NewList = []     for fc in fcs:         for row in gp.SearchCursor(fc):             fcValue = row.getvalue(Atts)             NewList.append(fcValue)     #print NewList          arcpy.CheckOutExtension("GeoStats")          # Load required toolboxes     arcpy.ImportToolbox("C:/Documents and Settings/ArcGIS/Toolbox.tbx")          # Local variables:     Cokriging_xml = "C:\\INTERP\\COK_.xml"     CK = "CK2"     ids1 = "C:\\INTERP\\data.gdb\\dem1000"               for fc in fcs:         for fcValue in NewList:             inputDset1 = fcValue          #Variables             inputDset2 = ids1            #Covariable             InputDset = "inputDset1;inputDset2"          # Create cokriging layer (code from model builder)              tempEnvironment0 = gp.autoCommit         arcpy.env.autoCommit = "1000"         tempEnvironment1 = gp.spatialGrid1         arcpy.env.spatialGrid1 = "0"              arcpy.GACreateGeostatisticalLayer_ga(Cokriging_xml,InputDset,CK) ##Function              arcpy.env.autoCommit = tempEnvironment0         arcpy.env.spatialGrid1 = tempEnvironment1      I need help. If someone has an idea? It would be great. Thanks


title: Build Raster Attribute Table with GDAL Utilities

text: Is there a way to build a raster attribute table using one of the GDAL utilities? I know you can view it using `gdalinfo` and I have successfully created one using Python, but it would be so much easier to just run a command line utility to build a raster attribute table on a raster.


title: Edit data using arcobjects by SOE Web Services on Arcgis Server 10.0

text: I trying to develop WebService which main purpose is enable editing data. My Web Service get as param information about records, tables etc. Code looks like below:               IFeatureWorkspace ifw = null;             IWorkspaceEdit iwep = null;             try {                 ifw = new IFeatureWorkspaceProxy(new IVersionProxy(currentUserVersion));                 iwep = new IWorkspaceEditProxy(ifw);                 iwep.startEditing(false);                 iwep.startEditOperation();                 performOperation(mergeData, ifw, iwep);     //this method make store on records             } catch (AutomationException e) {                 e.printStackTrace();             } catch (IOException e) {                 e.printStackTrace();             } catch (Exception e) {                 e.printStackTrace();             } finally {                 try {                     iwep.stopEditOperation();                     iwep.stopEditing(true);      } catch (AutomationException e) {                      } catch (IOException e) {                      }             }      Everything works fine for one single request, but above code is horrible slow for multirequest for different version. Because different people has own version which is used to editing. I think the most time take startEditOperation and startEditing. Everything starts as capability of MXD file. Anyone can give me a solution how to using soe web services give capablities for editing data which will have a quite good performance?


title: Where can I get aerial images for commercial use?

text: Related question: Free aerial map images for commercial use - Alternatives to Google maps I'm looking for a service that will allow me to download images at up to 1 meter resolution. Currently our system is based around the Mercator projection (pardon me if my terminology is incorrect) documented here. On this scale, we need imagery up to zoom level 17, which works out to be 1.2 meters per pixel. The only quote online that I could find was from Google Maps API for Business, which is $10,000 at the time of this writing. From what I can tell, the license doesn't allow images to be downloaded. Our system doesn't guarantee connection to the internet, so images will be pre-loaded. Currently I'm using MapQuest Open Aerial Tiles, but I need a little more coverage that that service provides. If it matters, we don't charge for the images, but we do charge for the software that uses the images. Images are only a small part of the overall software, so we consider it to be a free service provided with the product. What services exist that meet these requirements? I am willing to pay a license fee, provided it isn't ridiculous (<20,000 or so per year).


title: Where can I find the WKT standard or reference?

text: I had a little confusion about the representation of a geometry as well-known- text (WKT). I was searching for the WKT Specs to search for the exact representation, but I did not find the Standard. Is there a standards document or reference for WKT? where can I find it?


title: How to export attribute table to Excel using ArcGIS Desktop?

text: In previous versions I could just export as a .DBF file, open it with Excel and then convert it to an Excel file. However, now in ArcGIS 9.3.1 the .DBF file can no longer be open by Excel. I have also tried saving as a text file but to no avail.


title: Export table from ArcGis to Excel file

text: Is there a way to export a table with selected records to an Excel spreadsheet? I was trying to export just to .txt file but ArcGIS gives me an error message. The other problem I have is with coordinates. I want them only with one locus after the comma, I mean like the format: xx xx xx,x not like I have now xx xx xx,xxx. ![alt text](http://i.stack.imgur.com/AnQLd.jpg) Do you know any solutions? The version of my software is 9.3


title: GDAL minidriver and Tilecache

text: I'm trying to use GDAL minidriver to read tiles folders generated from Tilecache seeding. I have set this in my GDAL source `<ServerUrl>http://myserver.com/test/${z}/${x}/${y}</ServerUrl>`, and when I try to use gdal_translate with this minidriver file I have 404 error:               0ERROR 1: GDALWMS: Unable to download block 0, 0.     URL: http://myserver.com/test/8/0/0     HTTP status code: 404, error: (null).      The error is right because I have another kind of structure. In example I have `http://myserver.com/test/08/000/000/.` It seems to me that the output of tilecache seeding process does not have a standard TMS structure, and for this reason I have these kind of errors. Am I wrong? My source `SRS is EPSG:3004`. Could it be this the reason of a non standard TMS folders structure?


title: Existing ArcGIS Rest Services - What platform to use for building mapping application?

text: I have clients that have an existing ArcGIS Server. Map layers are published as rest services and I am going to build a mapping application that displays their datasets. Functionality should include users' ability to click on layers and retrieve attribute info. I've been using Google's JS API v3 for most of my projects. I am also using leaflet and CartoDB for other projects. I am trying to stay away from ESRI's JS API for server as it seems a little cumbersome. I'd prefer to use something like Google Maps or Leaflet for this project. Any suggestions as to what platform I should use to construct this application?   • Will it make my life easier to just use ESRI's JS API for server? • Any major downfalls to loading the REST services in Leaflet or Google ?


title: Export SRTM30 Dem file as 3d

text: i have downloaded some SRTM Dem files here: http://dds.cr.usgs.gov/srtm/version2_1/SRTM30//srtm30 They are from the USGS site. I would like to export them as 3d. I have ran through a few tutorials online, eg 1) http://www.brocku.ca/maplibrary/proc...iew/3dview.htm 2) http://prod.library.utoronto.ca/maplib/3D.pdf But still don't get it. Question: 1) so are the SRTM DEM files exportable to be viewed in 3d via Arcscene/ Arcmap? 2) how do you do it? Thank you!


title: Need advice on hardware/os upgrade

text: Forester, not a geopspatial or IT professional, currently using QGIS and SAGA for generally county level geoprocessing of multipspectral imagery and DEM's (cluster analysis, zonal statistics, stream channel modeling). Currently exceeding practical processing limits for laptop/windows configuration using HP EliteBook w/4GB Ram, Intel Core i5, 32 bit Win7 OS. Need advice on logical next step in hardware/OS/software. Future ambitions include deployment of PostGIS and derivation of canopy height models from stereo matched imagery. Thinking low end workstation and Ubuntu OS. Guidance/advice/considerations and recommendations much appreciated. I'm a one man shop trying to operate a bit beyond my capabilities. I do not speak fluent IT or GIS, need a balance of horsepower, idiot proof, moderate learning curve and intermediate support. Cost is important but functionality, flexibility and capability are primary. Thanks.


title: how to display polygons as points until certain zoom level?

text: can somebody tell me how to display polygons as points until you zoom to some level. I'm speaking about OpenLayers. So, for example, if the zoom level is less then 10 than it should display feature as a point and when you zoom to 10 or further it should display it as polygon with the shape it really has. I've tried something like this:               var my_style = new OpenLayers.Style({     'fillColor': '#9900cc',     'fillOpacity': .8,     'strokeColor': '#330066',     'strokeWidth': 1,     'pointRadius': 4     });     var vector_style_map = new OpenLayers.StyleMap({     'default': my_style     });          map.events.register("zoomend", map, function () {         if (map.getZoom() < 10) {             console.log("zoom level is less then 10");             layer.styleMap = vector_style_map;         }     });      But it just paint my polygon with defined colors, it doesn't change the shape to point. I'm probably missing something obvious but I would appreciate if you could help me. Thanks


title: How to update polygon longitude if negative?

text: I have a polygon, with longitude in [-180, 180], with 2 parts around the dateline. I want to make a map centered on the date line. My plane is to: 1./ update the longitude, so that if lon < 0 , lon=lon+360 2./ dissolve the two part to get a unique polygon. My problem is now 1./ I'm trying:               ogr2ogr nodateline\new_polygons.shp org_polygon.shp -dialect sqlite -sql "SELECT ShiftCoords(geometry, 360, 0) FROM org_polygon WHERE X(geometry)<0"      and the command returns the following error:               Warning 6: Normalized/laundered field name: 'ShiftCoords(geometry, 360, 0)' to 'ShiftCoord'      EDIT: I also tried with `ST_X(geometry)` instead of `X(geometry)`: same error. The error only occur with the "WHERE" clause. a./ what does the error mean? b./ what would be the correct syntax?


title: Is IGraphicTracker leaking memory?

text: It seems that when geometries are removed from an IGraphicTracker instance, the geometries remain in memory until IAoInitialize.Shutdown() is called. I have tried removing the geometries by calling the IGraphicTracker.RemoveAll() method and by storing the IDs returned by IGraphicTracker.Add() method and calling the IGraphicTracker.Remove() method for each ID. Is there something that I'm missing in the example provided below? I'm currently using ArcGIS 10.1 and Visual Studio 10 targeting .NET Framework 3.5. An example can be found here: http://forums.arcgis.com/attachment.php?attachmentid=20140&d=1356115382 This problem was originally noted on the ESRI Dev Forum: http://forums.arcgis.com/threads/72516-Memory-Leak-with-Graphic-Tracker


title: OGR: How to save layer from AttributeFilter to a shape filter?

text: I select a subset of features from my shape file like so:               node_lyr = nodes.GetLayer()     q_str = "%s = %i" % (field, value)     node_lyr.SetAttributeFilter(q_str)      What is the best way to make a shape file from this subset?


title: What is a non-geographical feature?

text: When using the term "feature" in GIS, what is often (if not always) meant is a _geographical_ feature which "is a representation of real world phenomenon associated with a location relative to the Earth" (ISO 19125-2:2004). Are there non-geographical features and if so, are there any examples of such? For example, could a compass or a scale bar be considered non-geographical features?


title: Conversion from NZMG to latitude and longitude

text: I want to convert `NEW ZEALAND MAP GRID` coordinates i.e. `Northing` and `Easting` into `WGS84` coordinates i.e. `Latitude` and `Longitude`. I have searched on internet but there is no proper explanation of how to do this, or an online calculator to do this. My final goal is to write a program in `C#` or `JAVA` , which will convert `NZMG` coordinates into `WGS84` coordinates.


title: What is the most "logical" route generating site (for Europe by car)?

text: What is the most logical route generating site do you think is out there? By logical I mean optimizing the driving time between two points. (Not always is the fastest way to proceed found along the shortest routes, making use of main-ways and highways when needed.) I tried Google Maps. It seems to make use of highways a lot. Bing Maps, Mapquest, and Map24 all have different calculation methods but they seem to not use highways as much as Google Maps. I mentioned "for Europe by car" because I would like to take a trip from _københavn(denmark)_ to _tekirdağ(turkey)_.


title: How to copy attributes from one polygon layer to another?

text: I have a problem that I can't seem to get my head around. I have two polygon layers:   * Polygon A - is a subset of polygon B with the same fields and has identical polygons to Polygon B   * Polygon B - has the attribute data I want to be in Polygon A How can this be done? I tried the QGIS tool "Join Attributes by Location" but as some of the polygons are within others, it tends to link to the first intersect it finds (the outer polygon).


title: How would I convert this point to WKID 4326?

text: This is a point in Brewster, MA, USA. How would it be converted to WKID 4326? 324299.3348436486, 832511.5939650714


title: Looking for advice on getting GIS data into existing AutoCAD site survey

text: I've got a variety of GIS data from Albemarle County in ShapeFile and DWG format. I'm looking to get the topographic layers and property boundaries together with my existing site survey created by a local surveyor. What's the best practice for doing this? Is it "just XREF it all into one master file" and scale/move around or is there a better way to make this data normalized? I'm in one of the corners of the grid squares that the county has created, so I'd probably need to have 4 of them combined. Edit: I have AutoCAD 2013. I have very little true GIS experience, but I have friends who have offered to help with QGIS.


title: How can you insert a table record without geometry into a feature?

text: I am working with ArcGIS and Python and trying to append table rows into a featureclass. You cannot use append and you cannot just make a tableview of the featureclass and then try to append it in. Any ideas? Thanks in advance.


title: How to draw rotated symbol RASTER type layer in Mapserver?

text: [PARTIALLY SOLVED]: I was on the wrong track I need to use a UVRASTER, see this part of the docs But doing leaves with with another question about how to get the angle correctly, please see this question I have a raster data source (GRIB2 file) which contains wind direction data for the whole world at an 1 degree interval. I use this mapfile:               MAP         EXTENT      -180 -90 180 90         CONFIG      "CPL_DEBUG" "ON"         CONFIG      "MS_ERRORFILE" "/var/www/mapserver/routeguard/poc.log"         SYMBOLSET   "symbols.sym"         FONTSET     "../fonts/fonts.txt"         IMAGETYPE   PNG8         IMAGECOLOR  255 255 255         SIZE        1600 900          # ------------------------------------------------------------------------------     # the Map Projection     # ------------------------------------------------------------------------------         PROJECTION             #which to use????             "init=epsg:4326" #54004 #4326             "units=m"         END          #----------------------------------------------------------------------------     #grib wave direction layer example     #----------------------------------------------------------------------------         LAYER             NAME "grib_wave_direction"             TYPE RASTER                  DATA "/var/www/mapserver/data/grib/routeguard/4nwcurpol99_2013110412.grb2" #Current direction @ Mean sea level             PROCESSING "BANDS=1.1"                  CLASS                 STYLE                     SYMBOL 'current_S'                     ANGLE [pixel]                     MINSIZE 20                     MAXSIZE 30                     COLOR 0 0 180                 END             END  # CLASS                  PROJECTION                 #"init=epsg:4326"                 "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +pm=-360"             END         END     END #end map file      This generates the following output: ![output from mapserver](http://i.stack.imgur.com/Y0oZj.png) Obviously, the symbol is not drawn, and also not rotated. **How can I draw the symbol and rotate it based on the pixel value?** The symbol is defined in my symbols.sym like so:               SYMBOL         NAME "triangle"         TYPE VECTOR         FILLED TRUE         POINTS             0 1             .5 1             1 1             0 1         END     END      Thank you! I'm on MS 6.4


title: Ogr2ogr does not support file name with non-ASCII character?

text: I tried to use the following script to convert shapefile to mif:               ogr2ogr -f "MapInfo File" D:\docs\arcgis_data\shp_mif\BOU.mif D:\docs\arcgis_data\shp\BOU.shp      it works. But I found that if the name of the shapefile has non-ASCII characters like this:               D:\docs\arcgis_data\shp\臺北.shp      Then the covert will fail with message:               Unable to open datasource `D:\docs\arcgis_data\shp\臺北.shp' with the following drivers     ......      Any idea to fix it?


title: OpenLayers: make polygons not zoom

text: I have a simple OpenStreetMap for which I want to display a series of circles. Each circle represents readership for that area, so its radius will change based on the numbers. However, when I create the circle and zoom in the map, the circle zooms in as well - it seems like OpenLayers assumes the circle is "glued" to the map and represents a piece of it, when really it's more like a programmatically created marker. How can I get the circle to stay the same pixel radius no matter what the zoom level of the map? Or is there some way to get a marker to be created as a Vector instead of an image?               var map = new OpenLayers.Map({       div: 'geomap',       projection: "EPSG:900913"     });     map.addLayer(new OpenLayers.Layer.OSM());          var epsg4326 = new OpenLayers.Projection("EPSG:4326"); // WGS 1984 projection     var projectTo = map.getProjectionObject(); // The map projection (Spherical Mercator)     // test point     var lonLat = new OpenLayers.LonLat(-74.265982, 40.665443).         transform(epsg4326, projectTo);          var vectorLayer = new OpenLayers.Layer.Vector("Overlay");     var point = new OpenLayers.Geometry.Point(lonLat.lon, lonLat.lat);          var mycircle = OpenLayers.Geometry.Polygon.createRegularPolygon         (             point,             50000,             40,             0         );          var featurecircle = new OpenLayers.Feature.Vector(mycircle);     vectorLayer.addFeatures([featurecircle]);          map.addLayer(vectorLayer);     map.zoomToExtent(vectorLayer.getDataExtent());      UPDATE: Here is the final answer (thanks Rene)!               var map = new OpenLayers.Map({         div: 'geomap',         projection: "EPSG:900913"       });       map.addLayer(new OpenLayers.Layer.OSM());            var epsg4326 = new OpenLayers.Projection("EPSG:4326"); // WGS 1984 projection       var projectTo = map.getProjectionObject(); // The map projection (Spherical Mercator)       var lonLat = new OpenLayers.LonLat(-74.265982, 40.665443).           transform(epsg4326, projectTo);            var vectorLayer = new OpenLayers.Layer.Vector('FSA', {         styleMap: new OpenLayers.StyleMap({           "default": new OpenLayers.Style({             pointRadius: "${radius}",             fillColor: '#ff0000',             fillOpacity: 1.0,             strokeColor: '#000000',             strokeWidth: 2,             strokeOpacity: 1.0           })         })       });       map.addLayer(vectorLayer);            var feature = new OpenLayers.Feature.Vector(           new OpenLayers.Geometry.Point(lonLat.lon, lonLat.lat),           {             radius: 10           }       );       vectorLayer.addFeatures([feature]);            map.zoomToExtent(vectorLayer.getDataExtent());


title: How to generate geodesic linestring from two points in C#

text: I'm looking for a **_FREE_** C# class, library or set of functions that will allow me to generate a **_geodesic_** linestring (or array of vertices, whatever) from two decimal degree points. I've found javascript functions that do what I want, but no luck with C# ones. I suppose I could translate it from the above javascript. But I'd prefer to use something tried and tested. Anyone know of such functions, classes or libraries?


title: Finding municipal-level administrative boundary data for India?

text: Does anyone know how to find a shapefile with boundaries of Indian municipalities? I have tehsil boundaries as for year 2001, but I would like boundaries of municipalities and, ideally, also historic ones. I believe Census of India has the information but I couldn't find any maps or shapefiles...


title: Create Friendly Labels for ArcGIS Print Service Options with JS API

text: How can I change the print options received from an ArcGIS print service to be more friendly and descriptive using the Esri sample print service?


title: How to create Online IDW(Inverse Distance Weighted) Interpolation?

text: I have some point data which store as JSON.As follow: `[["Lat":"12.12","Long":"23.23","ZValue":"123"],...]`. I tried to interpolate them as contour line and displayed them on my website ? I do know how to interpolate data in ArcInfo,so I tried to publish a geoprocessing task and using ArcGIS JavaScript API to generate them but it failed to published to ArcServer.![enter image description here](http://i.stack.imgur.com/44kp3.jpg) However,I search for a long time if there exits some other solution but found nothing about how to implement this. Thanks for any suggestion!


title: How to create a new service area analysis?

text: I need to create and run a new service area analysis. I'll need to do an analysis of how far you are coming by car from one hospital in 30 min, 60 min and 90 min. Is there anyone who can help me with a model of this in ModelBuilder?


title: Use Python to parse string

text: OK, I know that this is Python 101, but I can't seem to figure it out. I have a function setup that will return a string, specifically an azimuth in DMS. I need to take the return values and pull out only the value that appears prior to an ":" So I may get a value x = 123:34:12.4 I need x to then equal 123 and use that to updateRow in an attribute table. As usual... THANKS!


title: What is the most appropriate Projection for the Gulf coast region?

text: What are the most appropriate datum and projection I could use for the Gulf coast region - spanning Louisiana, Mississippi, Alabama, and Florida? I have to analyze block group data along the coast. Thanks


title: How can I edit and save feature attributes of a WMS layer in OpenLayers?

text: I am creating a data editing application in OpenLayers where the user will need to be able to edit the attributes of features (which are stored on a GeoServer installation). I have seen documentation for editing features in a WFS layer, but I haven't found anything detailing how to edit the attributes of a feature in a WMS layer. Is this feasible to do using a WMS layer? I am currently retrieving data from GeoServer using the GetFeatureInfo request. I would like to load this data into a div, modify it, then save it back to the server. I haven't seen any examples of OpenLayers maps doing this, but I feel like it should involve a simple POST request to the server.


title: Update GEOS to 3.4

text: I want to use the rgeos function gDelaunayTriangulation (which requires GEOS 3.4) on my Mac, but when I load rgeos it says:               > library(rgeos)     rgeos version: 0.3-2, (SVN revision 413M)     GEOS runtime version: 3.3.3-CAPI-1.7.4       I have installed and uninstalled both rgeos and GEOS (GEOS framework v3.4.2-2 from www.kyngchaos.com) a number of times and in different orders, but no matter how I try do I get rgeos to say the runtime is anything else than 3.3.3. What should I do to get GEOS 3.4?


title: How to vary the infoWindow size for different feature layers

text: I have multiple feature layers on the map, and I wish to display a different sized infoWindow for each layer. As far as I can tell, the map has one infoWindow, which can be resized using: `map.infoWindow.resize(width,height)` Each feature layer then provides the infoTemplate which contains the Title and Content, but no sizing information (right?). How can I show a large infoWindow when clicking on FeatureLayer1, and a small infoWindow when clicking on FeatureLayer2?


title: Title data from KML/SHP to google maps

text: Im trying to take land title data from a Shape file and convert it so it will display on Google maps. Works fine going to Google earth I can get the plots and be able to select each one and pull up the relevant data. But I am unable to find anything about doing this on Maps, so far all I have managed is to show the plot data, but cannot call up any data on it. I have been hunting for solutions for a week now and Im at my wits end. Any one have any ideas? Thank you


title: ST_Geometry vs SDE Binary: seeking recommendations, experiences

text: I've got some questions I've collected as we try to come to a decision about whether to use ST_Geometry or SDE native geometry for an Oracle/SDE datastore on a new project. The organization has quite a few legacy systems that use the SDE Binary type, and this would be the first to use ST_Geometry if we went that route. Apart from the ability to access spatial data in ArcSDE with SQL, what are the benefits and drawbacks of ST_Geometry? What sort of pitfalls have you experienced in implementing with one type or the other? Can more than one data type live in the same schema? Note: for this project the data in SDE will be all points and polygons.


title: Can I color different layers with the same pseudocolor color map?

text: I have multiple DEM layers in my map (downloaded from atlas.louisiana.gov/lidar), but I cannot seem to get the colors to line up. You can see in the attached image that the colors are not continuous between the layers, but they should be. (I've checked the values, and they are). I've tried several ways around this, but none of them seem to work. This data is rendered as Single Band Gray. 1) right clicking on one layer, selecting copy style, and then pasting it into the other layer. 2) saving the style from one layer and loading it into the other. 3) changing the color map to grayscle, inputting min and max, and then reverting back to pseudocolor. (that one seemed silly) It seems like the problem is that I can only do a pseudocolor with standard deviations, and not min/max values. Is that correct? There seems like there should be some way around that... Any ideas?


title: Distort image to display on leaflet.js based on known 16x16 point grid

text: I have an image that is a satellite picture, and with php gd I want to distort it so I can overlay it on a leaflet.js map. To make it a bit more simple, I have a grid of points, 16x16 which bind the image to the coordinates, like I chopped the image to 256 pieces and found appropriate coordinates for corners of each of those pieces. It's tedious and inefficient to draw 256 pieces simultaneously, hence I'd like to process the image on the server-side. Is it even possible to do it with just PHP's GD library or is it necessary to reinforce the arsenal? The projection I'm using (and would like to stick with) is standard leaflet.js EPSG:3857. How do I even start about such a task ?


title: ArcEngine draw a polyline or multipoint

text: Does anyone have a bit of code to plot a multiline? I just want to add some IPoints to a collection and display in a particular color and width, but I can't find a good way to do this. Right now I'm using this after a dynamic draw event: `dynamicDisplay.DrawMultipleLines(pointcollection);`   But there's no clear place to set the graphic attributes, and it plots incorrectly.


title: Can I Manually Add Wizards to a Map on CartoDB?

text: Some of my locations are not getting geocoded in spite of my checking spellings etc.. For instance, I'm trying to load some data on Chhattisgarh, an eastern Indian state and Lakshadweep, an island in the Arabian Sea which falls under the Indian territory. Need urgent help with this.


title: OpenLayers and the legend

text: Is it a possibility in OL to add a legend (symbology) directly to (under) particular layer (not only switch on/off the layer)?


title: QGis: Natural breaks (Jenks) - algorithm is buggy?

text: I don't know how the Natural Breaks (Jenks) algorithm works in QGis graduated symbolism, and I've deliberately not researched it in order to ask this question, as my point is based on what seems to be intuitively wrong, rather than what may be technically correct (but not useful or expected). So! If you have a range of values and constrain them to a few classes, Jenks does a nice job of splitting up the data:![A](http://i.stack.imgur.com/Mi2mz.png) But if you add further classes, Jenks just repeats the lowest value, often: ![B](http://i.stack.imgur.com/2ab2b.png) Since I don't know (yet - I will now research it) how Jenks works, it'd be remiss of me to say the algorithm is coded "wrongly", however the implementation is peculiar and unexpected. Possibly if it auto-constrained itself to the maximum number of independent classes? Or pooled identical- valued colour bands?


title: Inherited EditorServerCommandAction double posts returnMessage

text: I have a bunch of inherited EditorServerCommandActions in my WebADF solution and when ever I use the ReturnMessage inside EditorServerAction() it generates two javascript popups at the clientside. Is this a common know error or am I simply missing something?               public class CustomCommandAction : EditorServerCommandAction     {         protected override void EditorServerAction()         {             ReturnMessage("Foo");         }     }


title: Converting from DMS to DD using Python in Field Calculator?

text: I need to convert lat/long that is expressed as degrees, minutes, and seconds in the data _into_ decimal degrees. For example, in the data, they are listed as N335042.06 in the Latitude column, and W86031.04 in the Longitude column. I have done this problem before where I created a script which converted DMS to DD, and vice-versa, so I guess I could use bits from that. But the problem I'm having is how to ignore (for lack of a better word) the 'N' and 'W' in the data? Can I skip them? And DMS are listed all together without any symbols or spaces. Can I use `len()`, `range()`, `split()` to specify what part to read from the value? For example, can do the following? N335042.06 where, 33 = degrees 50 = minutes 42.06 = seconds ...? I came across this ESRI article, but it's in VB. Will probably use it as a reference but some of the terminology/syntax is different from Python. Final code that works!               # Pre-logic     def latDD(x):       D = int(x[1:3])       M = int(x[3:5])       S = float(x[5:])       DD = D + float(M)/60 + float(S)/3600       return DD          # Expression     latDD(!Latitude!)


title: How to convert a National Bridge Inventory(NBI) ASCII file to a shapefile

text: I have a problem. I have downloaded the NBI ascii file data from http://www.fhwa.dot.gov/bridge/nbi/ascii.cfm. Now I am not able to convert it into a shapefile. The ASCII file has a latitude longitude column. I understand that it can be imported in a GIS software using X,Y values. But I am not able to make out whether the latitude longitude is in DMS or decimal degrees. I am posting a sample here. Any help would be appreciated Sample Latitude: 40454200   Sample Longitude: 122190700 If there is any other way of converting it, I would like to know.


title: How to export data from DBF to Excel template using ModelBuilder?

text: I have created a model in ModelBuilder that uses the iterator tool to run through a polygon layer and clip a raster based on the polygon's unique name in the attribute table. The model then copies the rows of the clipped raster's attribute table and creates a dbf table saved with the polygon name. I now need the dbf table data to be exported or copied and pasted into an excel workbook that contains formulas in a separate worksheet to analyze the data. Is there any way to do this in ModelBuilder? Python? I have used the do.TransferSpreadsheet command in VBA to transfer data from MS Access to an Excel Template in past.. The command works great for automated analysis of complex data and would like something similar for ArcGIS.Thanks in advance!


title: Implementing WFS within .NET

text: There are several discussion about implementing (or not implementing) WFS with .net from scratch. However, if I do need to, how do I get started? I know it's a broad question. I apologize for that. I need help getting started. Any kind of walk through will be helpful. edited (8/2/13): thank you for all the comments and responses. I am looking for a step by step walkthrough with details assuming I'm just a GIS person with limited knowledge to a lot of programming concepts. Also, if this is not the right place to ask this question, please let me know.


title: Umlaute not displayed correctly when loading .mdb MS Access into QGIS (1.9.0-Master) with eVis

text: Good evening, I am loading an .mdb (MS Access 97 or 2003) to QGIS (1.9) via the eVis database conection. Connecting to the mdb, selecting a table and displaying it as a point-layer through the appropriate x/y-fields works fine. However in the attribute table the German Umlaute ä, ö, ü and also ß are displayed as inverted "?". After days of searching I don't really know how to tackle this one. I am grateful for any recommendation on how to fix this. The umlaute in the mdb / attribute table are displayed correctly in ArcMap. I would love to get this running in QGIS. Thank you.


title: Entity Framework & Sql Server 2012 spatial

text: I'm struggling to find any concrete performance of Entity Framework paired with SQL server 2012. We're looking at storing nearly a million points in a database an right now the most common scenario would be to find all points within radius of N. I've mocked up a small prototype using EF5 and it seems fine. However, a senior developer has mentioned that in his opinion EF is bloated and I should consider using a micro-orm or even pure ADO.NET as EF is a bloated monster (lol). Performance is crucial on this project and I would like to know if anyone has experience to prove this guy otherwise (or right for that matter)?


title: How to adjust raster size for usage in raster calculator?

text: I want to make a visibility analysis of planned wind turbines in my region. I have a DEM raster and a polygon describing the possible area where the wind turbines might be set up. People are complaining, that the view to a prominent mountain in the region will be disturbed. My "simple" plan is, to raise the DEM-values covered by the polygon up to 200 meters and then perform a visibility analysis from the top of that mountain to be able to identify the regions, where the wind turbines might "block" the view. Problem: I'm not able to successfully combine the vector and the DEM. Steps taken so far:   * Polygon was clipped from the DEM   * Raster values from the clipped part were set to 200 What doesn't work: When using the raster calculator, trying to add DEM + clip, the resulting image only covers the clip. As i understand, the extend of the two rasters should be the same, and that NoData values will prevent operations (googling for 3 hours without success). I fail in resizing my clip to the DEM and also do not get rid of the NoData values.(Really tried every tool under Raster to no avail) Unfortunately, the rastercalc tutorials i found do not cover these "basic" things, so i'm completely lost at the moment. Thanx Bernd p.s. using latest qgis master from osgeo4w


title: Can I make an ESRI checkout replica from a checkout replica?

text: If I create a checkout from a clients geodatabase, can I then turn around and put that replica on our ArcGIS Server and create a couple replicas from it? So, then, our replicas would be synced to the "replica" on our server which would then be exported back out to a personal geodatabase. Would that personal geodatabase be able to be synced back to the clients database or would doing all this ruin what makes it a replica? I hope this made sense. Basically, what I'm trying to accomplish is getting a replica from a client that multiple people will be working on, but before syncing it back to the client, we would like the replicas we worked on synced back to 1 database so we can check the connectivity before giving it back to the client.


title: How to compile QGIS 2.4.0 on Ubuntu Precise 12.04

text: Been trying to compile this for a couple of hours, using this download: http://qgis.org/downloads/qgis-latest.tar.bz2 and these instructions: http://htmlpreview.github.io/?https://raw.github.com/qgis/QGIS/master/doc/INSTALL.html#toc3. I have problems with `libproj-dev` - I can't install it, get this error: > The following packages have unmet dependencies: libproj-dev : Depends: > libproj0 (= 4.8.0-3~precise2) but 4.8.0-3~precise5 is to be installed E: > Unable to correct problems, you have held broken packages. So I installed the headers by hand - compiler accepted them - and have tried to link in /usr/lib/libproj.so as `PROJ_LIBRARY` in ccmake. The compiler doesn't like this, I get: > make[2]: *** No rule to make target`/usr/lib/libproj.so', needed > by`output/lib/libqgis_core.so.2.4.0'. Stop. make[1]: *** > [src/core/CMakeFiles/qgis_core.dir/all] Error 2 I also get 'Could not find GRASS' a little earlier in the build - don't know if this is significant. I'm a bit confused, `libproj.so` shouldn't be a target? Anyway, can anyone help? My ppa is the ubuntugis stable one - maybe I need unstable?


title: PostGIS 2.0 "No function matches the given name and argument types"

text: After a lot of research, I figured out how to build PostGIS 2.0 from source on a Turnkey Linux LAPP (Postgres 8.4) server. The main draw of PostGIS 2.0 for me is the ST_ConcaveHull function, aka "shrink wrap". Now that I have a functioning PostGIS 2.0, I'm trying to implement these functions. I'm working with the NWS Shapefiles for States and Forecast Zones. Most if not all of the geometries I'm dealing with are multipolygons. I've figured out how to draw these polygons, and this is the query I use to request data in the GeoJSON format:               SELECT ST_AsGeoJSON(geom) as geom FROM states WHERE id = 5;      Unfortounately, because I'm working with State sized data, some of the requests are HUGE! The GeoJSON for New Jersey is 2mb! I know PostGIS has some great functions to reduce the complexity of the geometries, but I'm having trouble implementing them. The functions I want to use are:   * ST_Simplify    * ST_Union    * ST_ConcaveHull Unfortunately, I'm getting some very unhelpful error messages when I run queries using these functions.               SELECT ST_AsGeoJson(ST_Simplify(geom)) as geom FROM states WHERE id = 5;          Error: SQLSTATE[42883]: Undefined function: 7 ERROR: function st_simplify(geometry)          does not exist LINE 1: SELECT ST_AsGeoJson(ST_Simplify(geom)) as geom FROM states W...      ^ HINT: No function matches the given name and argument types. You might need to add      explicit type casts.      I went through a LOT of time and effort to make sure that I had PostGIS 2.0 instead of 1.5 for the explicit purpose of using things like ST_ConcaveHull and ST_Simplify, and now Postgre is telling me they don't exist! I found another post on here that gave me a query listing all global functions available to Postgre:               SELECT  proname, proargnames, prosrc      FROM    pg_catalog.pg_namespace n     JOIN    pg_catalog.pg_proc p     ON      pronamespace = n.oid     WHERE   nspname = 'public'      Listed in the results are:   * st_union {geom1,geom2} geomunion   * st_union pgis_union_geometry_array   * st_union aggregate_dummy   * st_concavehull {param_geom,param_pctconvex,param_allow_holes}   * st_simplify LWGEOM_simplify2d   * st_simplifypreservetopology topologypreservesimplify So I'm kind of stuck and confused at this point. The query results say the function doesn't exist, but the output of the above query says they DO exist. Can anybody shed some light on this for me? What am I doing wrong?


title: Add a popup to a Vector LineString in OpenLayers?

text: Is it possible to add a popup (without a marker icon) to a Vector (LineString)? This is my code, if possible, how to change it to add the popup? I need the popup to display some additional information to the road segement.               <script>     var ytrajectory=[];     ytrajectory.push('LINESTRING (11.1119339 46.0772706, 11.1117177 46.077171, 11.111613 46.0770837, 11.1115376 46.07702)');     ytrajectory.push('LINESTRING (11.1115376 46.07702, 11.1114585 46.0768611, 11.1114219 46.0767673)');          ...     </script>     function showReconstructedTrajectoryy() {     var vectorLayer = new OpenLayers.Layer.Vector(     "Reconstructed trajectoryy", {     style : reconstructedTrajectoryStyleYellow     });          var parser = new OpenLayers.Format.WKT();     var src = new OpenLayers.Projection('EPSG:4326');     var dest = new OpenLayers.Projection('EPSG:900913');          var s, wktLine, features, reprojectedFeatures;     for (s = 0; s < ytrajectory.length; s++) {     wktLine = ytrajectory[s];     features = parser.read(wktLine);     reprojectedFeatures = features.geometry.transform(src, dest);     vectorLayer.addFeatures(new OpenLayers.Feature.Vector(     reprojectedFeatures));     }          osmmap.addLayer(vectorLayer);     osmmap.zoomToExtent(vectorLayer.getDataExtent());     }


title: Legend in FlexViewer - Only show classes that are visible in the current map extent

text: I am using ArcGIS Viewer for Flex to create web map. The legend widget can be set so that the legend in web map to only show **layers** that are available at the current map scale. This is great but I want to set the legend to only show **_classes that are visible in the current map extent_**. Does anyone know how to do it? Thank you


title: How to add logical operators in wfs openlayers

text: How i can add logical operators AND, OR both togehter in WFS openlayers.                   <script type="text/javascript">         var map;              function init() {                  var bounds = new OpenLayers.Bounds(                     2309872.5, 2399960.25,                     2388258.25, 2450114.75                 );                     var options = {                     controls: [                                 new OpenLayers.Control.Navigation(),                                 new OpenLayers.Control.PanZoom()                          ],                     maxExtent: bounds,                     maxResolution: 306.1943359375,                     projection: "EPSG:32644",                     units: 'm'                 };                 map = new OpenLayers.Map('map', options);                 var  map_wms = new OpenLayers.Layer.WMS(                     "cite:union - Tiled", "http://localhost:8080/geoserver/cite/wms",                     {                         LAYERS: 'cite:union',                         tiled: true,                         EPSG:32644,                         tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                     },                     {                         buffer: 0,                         displayOutsideMaxExtent: true,                         isBaseLayer: true                     }                  );                            var layerPolygon2 = new OpenLayers.Layer.Vector("Layer Information", {                                             strategies: [new OpenLayers.Strategy.BBOX()],                                             visibility: true,                                             protocol: new OpenLayers.Protocol.WFS({                                             version: "1.0.0",                                             srsName: "EPSG:32644", // this is the default                                             url:  "http://localhost:8080/geoserver/cite/wfs",                                             featureType:"union",                                             featureNS: "http://www.opengeospatial.net/cite",                                             extractAttributes: true,                                             transparent:true,                                             geometryName: "the_geom"                                                       })                          filter: new OpenLayers.Filter.Logical({                     type: OpenLayers.Filter.Logical.AND,                     filters: [                         new OpenLayers.Filter.Logical({                         type: OpenLayers.Filter.Logical.OR,                         filters: [                             new OpenLayers.Filter.Comparison({                             type: OpenLayers.Filter.Comparison.EQUAL_TO,                             property: "classes",                             value: "A?Pt"                             }),                             new OpenLayers.Filter.Comparison({                             type: OpenLayers.Filter.Comparison.EQUAL_TO,                             property: "classes",                             value: "P"                             })                             ]                         )}                         new OpenLayers.Filter.Logical({                         type: OpenLayers.Filter.Logical.OR,                         filters: [                             new OpenLayers.Filter.Comparison({                             type: OpenLayers.Filter.Comparison.EQUAL_TO,                             property: "luse_2005",                             value: "BU"                             }),                             new OpenLayers.Filter.Comparison({                             type: OpenLayers.Filter.Comparison.EQUAL_TO,                             property: "luse_2005",                             value: "CL"                             })                             ]                         )}                          ]                      })                           });                                map.addLayers([map_wms,layerPolygon2]);                      map.zoomToMaxExtent();      } when i try to run above code in browser nothing is being displayed and firebug showing ReferenceError: init is not defined init() error. thanks!!


title: Print XYZ layer with MapFish Print plugin

text: I have a XYZ layer(i hope i call it right) in GeoExt/gxp application:               ,{source: "ol",group: "background",fixed: true,type: "OpenLayers.Layer.XYZ",args:        ["ЕЭКО","http://maps.rosreestr.ru/ArcGIS/rest/services/BaseMaps/BaseMapAnno/MapServer/tile/$\{z\}/$\{y\}/$\{x\}",{visibility: true}]}      In MapFish documentation says that print plugin can print TMS and Tile Cache. But its print only WMS layers. I gonna set url to XYZ layer or something? **UPDATE** I found working XYZ print in here http://map.geo.admin.ch/. But There is something strange with layer which set to print. **UPDATE2** More interesting that GeoExt.PrintProvider not send XYZ layer to print plugin. But if i give to it gwc layer i can see gwc layer to post parameters: ![enter image description here](http://i.stack.imgur.com/60ras.jpg) But this layer not print. Why PrintProvider not send XYZ layer in post parameters?


title: MapInfo 10.0 - saving a MapInfo browser to as a table

text: I geocoded a table in MapInfo 10.0. In the browser window, the lat and long are viewable, however the columns are not saved to the original table. The "Save Table" option is not available. Are there any options of saving a browser as a table?


title: How to use the MMQGIS plugin through the Python console?

text: I know how to make Hub Lines using the GUI, I want to know how to do the same using the Python Console, a tutorial would be great. Thank you


title: Calculate "CO_FIPS" field based on "COUNTY" field

text: I have several shapefiles with "COUNTY" already populated, I need to populate the "CO_FIPS" field based on the county. BREVARD = 009, BROWARD = 011 etc..


title: STIntersects - Update field to give me count of points in polygon

text: I have a points layer (dbo.ptLayer)   * Around 1M points   * Spatial Geometry Type (dbo.ptLayer.geom)   * No spatial index just yet, but will create one once data gathering complete. I have a polygon layer (dbo.polygonLayer)   * Around 500 polygons.   * Spatial Geometry Type (dbo.polygonLayer.geom) Both have fields called ID. **I want to populate an empty integer field in the polygon layer, with a count of the total number of points within each polygon.** Although I have access to other software products, I am interested to learn what can be done purely within SQL and SQL Server. I believe I should be making use of STIntersects but would like to know what is the best way of doing an update to populate this field.


title: Generating Point from Centroid of Polygon Error

text: I have a feature class of polygons within a geodatabase and I'm trying to generate points based on the centroid of the polygon using the Feature to Point tool in ArcGIS 9.3. However, I end up with several points that fall just outside of the polygon. So ... is this the wrong tool to be using for this application? Also, am I making the wrong assumption about this tool generating points from the centroids of polygons?


title: How to move/edit labels?

text: I need to move/edit few labels in a job. How can I work on it? I am using QGIS 1.7.4 and the label option in the toolbar (not the layer options)


title: Select Move Label Tool in QGIS

text: > **Possible Duplicate:**   >  How to move labels in QGIS? I would like to manually move 6 labels but the Move Label Tool is unable to be selected. I have placed labels by selecting the layer and clicking the "ABC" Labeling button. After going into Toggle Editing, I'm able to Change Labels but not Move or Rotate Labels. I'm running QGIS 1.7.4. I've looked at: How to move labels in QGIS? Is there a way to manually place labels in QGIS? It's a very simple problem but I can't seem to crack it. Thanks all, Zach


title: Table Manager - Shapefile (Labels)

text: I have a point shapefile which is an output of the plugin table manager. Now I’m labeling the features, but i can´t seem to get the move label bottom to work, and is active. Is there any connection? Thanks in advance


title: How to add points to a point vector layer - QGIS

text: I have a point layer with multiple data points in QGIS, but I want to add some extra data points to it and adjust some info that we missed because the ground was wet and we needed to go around it. I can't figure it out. I would think you would go in edit mode and have an option to add points, but can't find it.


title: How to add a title in a GetLegendGraphic Request

text: The GetLegendGraphic request works well. I display the legend at a fixed position in my viewer using OpenLayers. But I would like to know if there is any way to add some custom text to the legend. The legend that geoserver produces is like this: http://docs.geoserver.org/latest/en/user/_images/samplelegend.png But I need to produce a legend with a title: http://minus.com/m4AP102Jj Does anyone know if the GetLegendGraphic request in GeoServer, or the SLD style have a way to add this? **UPDATE (23/08/12)** For now I´ll use a litle trick to add a title by adding in my SLD another custom rule above all others that not match any other rules: Example:               <sld:Rule>       <sld:Name>[NAME]</sld:Name>         <sld:Title>[TITLE]</sld:Title>         <ogc:Filter>           <ogc:PropertyIsEqualTo>             <ogc:PropertyName>[PROPERTY_NAME]</ogc:PropertyName>             <ogc:Literal>[SOME_VALUE_THAT_NOT_MATCH_OTHER_RULES]</ogc:Literal>           </ogc:PropertyIsEqualTo>         </ogc:Filter>         <sld:PolygonSymbolizer>           <sld:Fill>             <sld:CssParameter name="fill">[MATCH_BACKGROUND_COLOR]</sld:CssParameter>           </sld:Fill>         </sld:PolygonSymbolizer>           </sld:Rule>


title: Does the USDA SSURGO database have a web-based API that I can query?

text: Based on a previous question, Can I download the entire NRCS SSURGO database?, it appears that it is difficult to import the entire USDA SSURGO database. The Web Soil Survey provides the ability to draw a polygon and then download data. But the polygon has to be fairly small, and there is a lot of pointing, clicking, and waiting to get the desired data. There is also a Geospatial data gateway where it is possible download county-level data. However, these do not meet my needs. This is what I want to do:   1. I start with the lat/lon coordinates of the corners of a rectangle   2. find the area of each soil type in the rectangle   3. find % clay, sand, silt, organic matter as a function of depth for each soil type  Is there a way that I can do this programatically? (Solutions in R preferred, but any open-source solution will be useful)


title: How can I consume map service using static IP address for VPN

text: I'm current using the url (http://services.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer) for my iPad map service. However, we need to connect to VPN to access the internal network at the same time. And the VPN server blocks the public Internet and only allow me to setup IP Address based public Internet access. So my question is there static IP addresses associated to the above url I could specify? Thanks


title: Calculate Distance between a Polygon and a Point Along a Path

text: I am very new to ArcGis and am trying to figure out how to do this. Basically, this is what I am trying to do: I have a map of Canada and in it I have polygons representing areas (each about the size of a municipality). I also have points describing the location of every airport in Canada and the network of roads all across Canada. **What I need to do is find the distance from each area to the nearest airport using a road.** What I tried to do already was find the centroid of each polygon because I thought that I could use that to calculate the distance between two points along a road. But then I realized that the centroid is not necessarily connected to a road (especially in the northern parts of Canada) therefore I cannot use that. I would really appreciate some sort of guidance as to how I would accomplish this as I am very new to ArcGIS and have not found much to help me (without confusing me) online.


title: Using VBScript to label shapefile in ArcMap

text: If you're reading this for the first time, it might be best to focus on section beneath **ANOTHER UPDATE**. I'm a beginner with VBScript and I have a question about using an VBScript expression for a label. I want to label street address numbers and the fractional address number, if it exists, in my map. This data is stored in an address points shapefile. I go Properties > Labels > Expression... in Esri software to enter this expression:               [HSE_NO]&" "& [FRAC_NO]       With this expression the labels will be created next to the address points as follows:                412 1/2 412      See how both labels are on the same line? I would like to have the labels look like this:               412     412 1/2      I want the first label to be the address value without a fraction, and the following lines would be address values with fractions. I welcome your suggestions. **UPDATE** I think what I need is more complex than I realized. As I mentioned, I want first label to be address street value (HSE_NO) without a fraction (where FRAC_NO is NULL). I'm working on getting this to occur with what's below. It's not correct! Any ideas?               Function FindLabel ( [HSE_NO], [FRAC_NO] )     if (IsNull[FRAC_NO]) then         FindLabel = [HSE_NO]          End If     End Function      **ANOTHER UPDATE** I have a point shapefile of addresses. Some address points contain the address number and street name. Some points contain address number, address fractional number, and street name. Points containing fractional values will overlap points containing not containing fractional values. I want to label these overlapping points like this:               1159     1159 1/2      My current VBScript labels addresses in this manner.               1159 1/2 1159      Does anyone know how to get the desired result? Here's my current expression:               Function FindLabel ( [HSE_NO], [FRAC_NO] )       if (IsNull([FRAC_NO])) then         FindLabel = [HSE_NO]       else         FindLabel = [HSE_NO] & " " & [FRAC_NO]        End If     End Function      Let me know if there are questions. Thank you.


title: Trouble exporting from QGIS to Illustrator

text: I am having trouble getting my files to transfer or export to a vector file. When I open the file in Adobe illustrator there is a lot of information missing. I have tried to "save selection as a vector..." I have all my layer selected. I tried selecting individual layers, all the layers at once, I selected the layers individually. Next I tried to play with the CRS pull down menu. I used all the options and still get no results. I finally got what may be a file but when i open it in illustrator I get an error message that says "file cannot be opened because there is not enough RAM,, I have 12 gigs of RAM, should not be a problem. Am I missing something? A plug-in?


title: How to change the raster NoData value with ArcObjects?

text: I am trying to change the a raster's NoData value with ArcObjects and keep coming up short. Here is the code that I thought would work:               IRasterDataset _dataset = ...get the dataset...     object[] noDataValue = ...set the no data values by raster band...     var rasterbandCollection = (IRasterBandCollection) _dataset;     for(int i = 0; i < rasterbandCollection.Count; i++)     {         ((IRasterProps) rasterbandCollection.Item(i)).NoDataValue = noDataValue[i];     }      This value persists during the ArcMap session, but does not save the changes to disk it appears. Any ideas on how to make this permanent (without having to save a second raster dataset)?


title: ArcSDE: Specifying fields to not allow null values

text: I have a feature class in SDE that has various fields where the setting AllowNulls is set to FALSE. When creating a feature in ArcMap in this feature class it still allows nulls to be specified in the attributes, albeit an empty text string. I can then save my edits without any error messages. What's going on? Thanks Rob


title: Adding Google map layer in Openlayers - what about this function?

text: How to add Google map as a layer (something like " _var gmap = new OpenLayers.Layer.Google("Google Streets_ ", {visibility: false})"?) to this _mymap.js_ function i already have:               var map = new OpenLayers.Map({         div: 'map',     //  allOverlays: true,         theme: null,         projection: sphericalMercatorProj,         displayProjection: geographicProj,         units: 'm',         numZoomLevels: 18,         maxResolution: 156543.0339,         maxExtent: new OpenLayers.Bounds(             -20037508.34, -20037508.34, 20037508.34, 20037508.34         ),         controls: [             new OpenLayers.Control.Attribution(),             new OpenLayers.Control.Navigation(),             new OpenLayers.Control.PanZoom(),             new OpenLayers.Control.LayerSwitcher()              ],         layers: [             new OpenLayers.Layer.OSM('OpenStreetMap', null),             sundialsLayer,             sprintersLayer,             poisLayer         ],         center: new OpenLayers.LonLat(0, 0),         zoom: 2     });      Thanks. Note: suggested answer is fine -                   layers: [             new OpenLayers.Layer.OSM('OpenStreetMap', null),             new OpenLayers.Layer.Google(             "Google Streets", // the default             {numZoomLevels: 20},             )                  sundialsLayer,             sprintersLayer,             poisLayer         ],


title: Make Table View with field_info Output

text: I am using Make Table View to copy data from one table to a table view object and in the process hide and/or rename some fields. When I output the data to a dbf it hides and/or renames the fields accordingly. When I output the data to a geodatabase, it hides the specified fields but does not rename the ones that should be renamed. Any ideas why this is happening? EDIT: After some further testing, it looks like it is changing the field name, and not the field alias in the geodatabase. I tried changing field.name to field.aliasName, but that did not work.               import arcpy     from arcpy import env          # Set the enviroment     env.workspace = '//GIS-SDE/ParcelData.mdb'     env.qualifiedFieldNames = False          # Get the fields from the G_AccountInfo table in ParcelData.mdb     fields = arcpy.ListFields("G_AccountInfo")          # Create a fieldinfo object     fieldinfo = arcpy.FieldInfo()          # Iterate through the fields and set them to fieldinfo     for field in fields:         if field.name == "ACCOUNT_NUM":             fieldinfo.addField(field.name, "ACCOUNT", "VISIBLE", "")         elif field.name == "OWNER_NAME1":           fieldinfo.addField(field.name, "OWNERNAME", "VISIBLE", "")         elif field.name == "OWNER_ADDRESS_LINE1":           fieldinfo.addField(field.name, "OWNERADDRESS1", "VISIBLE", "")         else:             fieldinfo.addField(field.name, field.name, "HIDDEN", "")          arcpy.MakeTableView_management ("G_AccountInfo", "AccountInfo", '//GIS-SDE/ParcelData.mdb', fieldinfo)          # Give the output table a unique name     output = '//GIS-SDE/ParcelData.mdb/cleaned_data'          # Add the table back to the gdb     arcpy.CopyRows_management("AccountInfo", output)


title: Map Info Workspace problem

text: I have opened a workspace with a file on a CD and now it will only open when the CD is in the CD drive - I can't find a way around it - any help gratefully received , thanks


title: map of europe in list of lat/long points

text: I'm looking for a list of lat/long coordinates which will describe all countries of europe. Ideally the list will also describe on which country each closed path belongs to. I'm imagining something along the lines:                lat    long    cc  type     ------ ------- ---- ----     55.05   -7.25   IE  EU27     55.18   -6.97   IE  EU27     55.25   -6.48   IE  EU27     55.22   -6.1    IE  EU27     55.05   -6      IE  EU27     54.83   -5.7    IE  EU27     54.63   -5.92   IE  EU27     ...              55.23   -6.93   IE  EU27     55.05   -7.25   IE  EU27          58.63   -5      UK  EU27     58.55   -4.67   UK  EU27     58.45   -4.77   UK  EU27     58.58   -4.6    UK  EU27     58.45   -4.48   UK  EU27     58.6    -4.02   UK  EU27     58.63   -3.02   UK  EU27     58.37   -3.12   UK  EU27     57.85   -4.2    UK  EU27     57.85   -3.77   UK  EU27     57.6    -4.42   UK  EU27     57.6    -4.1    UK  EU27     ...               Regarding accuracy, I'd be happy with anything over 5.000 points in total. Where should I look at? I will convert those coordinates to x/y using `fromLatLngToContainerPixel` of google maps and then draw the map using javascript.


title: How can I convert geomagnetic coordinates to geographic coordinates without doing the math myself?

text: I have some data which are in geomagnetic coordinates; that is, they are latitude and longitude, but in reference to the magnetic north pole in Canada rather than the geographic north pole. As you might guess, I'd like to match them to other data which are in other coordinate systems. (The data originated as a computation of the auroral oval; i.e., they originated as geomagnetic and there is no geographic source I could consult, even in principle.) The data will be stored in PostGIS. The easy way to do this is to find an EPSG code for the geomagnetic coordinate system, tag the data with that code, and then everything is magically reprojected to whatever I need at the time. However, I can't find any EPSG codes on spatialreference.org containing either _geomagnetic_ or _magnetic_. Various Google searches also turn up nothing. For the level of precision I need, any latitude/longitude coordinate system (e.g., WGS84) modified to have a different north pole is good enough. So if there's an easy way to define such a coordinate system and load it into PostGIS, that would be OK too. **Clarification** For example, suppose that instead of converting from geomagnetic coordinates to geographic, I simply wanted to convert from UTM 15N to WGS84. I could save the UTM object into PostGIS unmodified, and tag it with the EPSG code 26915. Then, when I later made a query, I could say to PostGIS, "please return results in EPSG code 4326", and PostGIS would convert to WGS84 automatically. Similarly, I would like to tag the geomagnetic data with some coordinate system and let PostGIS do the reprojection behind the scenes, without me having to call a `reproject()` function manually.


title: Map Projection - where to begin

text: First of all, I realize that this question can be considered like "open" or maybe "too wide" (and definitely very basic for all of you). However, I am completely new in this and I am not looking for more then several hits where to begin, which topics I should search for. Also, I search for some similar questions and I found few open questions with similar meaning (Java, GIS knowledge…What next?, The Best books from ESRI, Where can I find QGIS tutorials and resources?) so I hope someone can help. **The Issues:** * * * I have data in geographical coordinates which I would like to project them into geomagnetic coordinates. * * * What do I need for that? Is information on the internet enough to solve this or I should search for a book? What should I search for?


title: Also Cant see WFS-T layer with OpenLayers and GeoServer (Windows)

text: I am in the process of learning for doing WFST - OpenLayers in web. I followed this tutorial in the GisTutor. But i stuck with 2 problems   * the WFS-T layer is not appearing   * the feature created in client cannot be saved (no response when save button is clicked, and accordingly there is no additional row in the PostGIS table of WFST_TEST) I followed this GIS.SE question and I also do not configured the Proxypass as instructed here (also because I do not have apache running the GeoServer). In the latter part of the discussion between Janne Mattila with Iant, Iant suggests that: > don't run from file://somefolder/mytestpage.html always use > localhost:28080/geoserver/www/mytestpage.html this should fix your problem. > – iant and then Janne seems to be happy: > Nice and simple, that fixed it! I don't fully understand why it works like > that, though. Can you explain what happens in OpenLayers when I access the > page using file:// and when I use localhost:28080 instead? – Janne Mattila I finally able to call my html by `localhost:8080/geoserver/www/rinotestpage.html` in my browser, by putting my `rinotestpage.html` and the `THEME` folder (for bringing button PNG files to the app) into the `C:\Program Files\GeoServer 2.3.2\data_dir\www`. But my 2 problems in the above are persisting (WFS-T layer not appeared and feature created in client is not saved) My question:   1. Janne and Iant, how did you crack this? How to solve the problem (WFS-T layer not appearing)?   2. Janne, did you put your `testpage.html` in the same folder as i did ? (`C:\Program Files\GeoServer 2.3.2\data_dir\www`). For additional infos:   1. I am using : PostgreSQL 9.1.9.1 Windows 32-bit + PostGIS 2.0.3.1 + GeoServer 2.3.2 + Windows Server 2003 in VMWare Player 5.0.2   2. I can see in my firebug that there is "POST wfs" request to my `http://localhost:8080/geoserver/wfs`, but the request took forever and no response. I am guessing this is why my WFS-T layer is not shown. ![enter image description here](http://i.stack.imgur.com/Bb3De.jpg) many thanks /Rino


title: Trying to display some data in OpenLayers - it's not erroring, just not loading :(

text: I'm trying to display some KML data on a map, using OpenLayers. Nothing complex - but it's just not rendering. Here's my javascript code               var map;          function init(element, kml) {         map = new OpenLayers.Map(element, {             projection: 'EPSG:3857',             layers: [                 new OpenLayers.Layer.Google(                     'Google Streets', // the default                     { numZoomLevels: 20 }                 ),                 new OpenLayers.Layer.Vector('Kml data')             ],             center: new OpenLayers.LonLat(153.132, -26.987)                 // Google.v3 uses web mercator as projection, so we have to                 // transform our coordinates                 .transform('EPSG:4326', 'EPSG:3857'),             zoom: 10         });              // Grab the features.         var features = convertKmlToFeatures(kml);              // Display it.         map.layers[1].addFeatures(features);     }          function convertKmlToFeatures(kml) {         var kmlFormat = new OpenLayers.Format.KML();         return kmlFormat.read(kml);     }      Basically,   * Create a map with two layers.   * Convert the Kml to features.   * Display the features on the map.  I just don't see what I've done wrong :( Oh - on a side note, the KML data is fine, because when I tried it to visualise it on this website, it worked. ![enter image description here](http://i.stack.imgur.com/6cvfa.png)


title: Problem generating Geodatabase Diagrammer

text: I am using ArcGIS 10.1 and Visio 2010. When I am trying to generate the diagram using the geodatabase diagrammer it has the following to error messages: Source at Microsoft.Office.Interop.Visio.DocumentsClass.Add(StringFileName) at GeodatabaseDiagrammer.VisioDiagram.OpenVisioDrawing(String fileName, String pathName) in C:\agis\Tools\ArcGIS\10x\GeodatabaseDiagrammer\AddInSourceCode\Version1.2\GeodatabaseDiagrammer\Helpers\VisioDiagram.vb:line 102 Msg: File not found. Then it would pop up another dialog box and say Source: at GeodatabaseDiagrammer.VisioDiagram.SectionHeader() in C:\agis\Tools\ArcGIS\10x\GeodatabaseDiagrammer\AddInSourceCode\Version1.2\GeodatabaseDiagrammer\Helpers\VisoDiagram.vb:line 160 Msg: Object reference not set to an instance of an object.


title: QGIS regular points around each polygon in a layer (vs ALL polygons)

text: I have a number of buffered points (some where polygons intersect & are dissolved. I'd like to create regular spaced points grid around EACH polygon -- not ALL polygons. When I do Vector -> Research Tools -> Regular Points --- I have to choose a layer as an Input Boundary, but it doesn't seem to loop through each shape to create a GRID just for that shape, and instead creates one large grid around ALL shapes. The end result: ![enter image description here](http://i.stack.imgur.com/2X4iY.jpg) How can I get a point grid for EACH individual shape vs ALL?


title: How can I add attributes contained in JSON data to an OpenLayers popup?

text:
I am working on a website using jQuery, Javascript, and OpenLayers. I am adding markers to an OpenLayers Vector layer using the following code:               var pointLayer = new OpenLayers.Layer.Vector("POI Markers", {projection: "EPSG:4326"});     $.getJSON('json/DATA.json', function(data) {       $.each(data.transmitters, function() {         var pointFeatures = [];         var px = this.longitude;         var py = this.latitude;         // Create a lonlat instance and transform it to the map projection.         var lonlat = new OpenLayers.LonLat(px, py);         lonlat.transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));              var pointGeometry = new OpenLayers.Geometry.Point(lonlat.lon, lonlat.lat);         var pointFeature = new OpenLayers.Feature.Vector(pointGeometry, null, {             pointRadius: 16,             fillOpacity: 0.7,             externalGraphic: 'marker.png',         });              pointFeatures.push(pointFeature);         pointLayer.addFeatures(pointFeatures);       });     });      I am using the following code to select markers and create popups:                var selectControl = new OpenLayers.Control.SelectFeature(pointLayer, {       hover: true,       onSelect: function(feature) {         var layer = feature.layer;         feature.style.fillOpacity = 1;         feature.style.pointRadius = 20;         layer.drawFeature(feature);                  var content = "
Feature:
" + feature.id +                   "

Location:
" + feature.geometry +"
";                  var popup = new OpenLayers.Popup.FramedCloud(             feature.id + "_popup",              feature.geometry.getBounds().getCenterLonLat(),             new OpenLayers.Size(250, 100),             content,             null,              false,              null);             feature.popup = popup;             map.addPopup(popup);         },         onUnselect: function(feature) {             var layer = feature.layer;             feature.style.fillOpacity = 0.7;             feature.style.pointRadius = 16;             feature.renderIntent = null;             layer.drawFeature(feature);                  map.removePopup(feature.popup);         }     });     map.addControl(selectControl);     selectControl.activate();      As you can see, the popup contains information regarding the feature id and geometry. What I would like to do is add information that is contained inside the DATA.json file referenced in the first code block to the corresponding popup. I am using a large JSON file. The following is a sample:               "LOCATION_A": {           "latitude": 10.2070,           "longitude": 25.3215,           "altitude": 1.2,           "status": "CLOSED",           "serverAddress": "192.168.0.1"      }      I would like to add the latitude, longitude and status information to the corresponding popup so when the user selects a popup, they view this information. In addition, it would be useful to change the "externalGraphic" file to reflect the status of each LOCATION. Continuing to use the above JSON file, if a POI were closed, the "externalGraphic" file should be one PNG file, while if the POI were open, it should be another. I am using jQuery in this project and find it fascinating! Any assistance with this problem is greatly appreciated. Thank you all so much.


title: i have an unhandled exception in asp.net

text: i am going to create a website and load a gis map to it. i am using arcgis 10 and i created a web application in arcgis server 10 i created a new WebSit in visual studio 2010 and add a map resource manager. add a MapResourceItem for this. after that add a map control and set MapResourceManager property. i config ArcGISIdentity too. but unfortunately when i run application occurs this error > Error processing server request > > [SoapException: Error processing server request] > System.Web.Services.Protocols.SoapHttpClientProtocol.ReadResponse(SoapClientMessage > message, WebResponse response, Stream responseStream, Boolean asyncCall) > +4159 System.Web.Services.Protocols.SoapHttpClientProtocol.Invoke(String > methodName, Object[] parameters) +328 > ESRI.ArcGIS.ADF.ArcGISServer.MapServerProxy.ComputeScale(MapDescription > MapDescription, ImageDisplay MapImageDisplay) +117 > ESRI.ArcGIS.ADF.Web.DataSources.ArcGISServer.MapFunctionality.GetScale(Envelope > extent, Int32 mapWidth, Int32 mapHeight) +564 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.get_Scale() +794 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.setExtent(Envelope value, Boolean > addSetExtentCallbackResult) +2556 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.set_Extent(Envelope value) +39 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.validatePrimaryMapResource() +2586 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map._OnLoad(EventArgs e) +115 > ESRI.ArcGIS.ADF.Web.UI.WebControls.Map.OnLoad(EventArgs e) +37 > System.Web.UI.Control.LoadRecursive() +134 > System.Web.UI.Control.LoadRecursive() +232 > System.Web.UI.Control.LoadRecursive() +232 > System.Web.UI.Page.ProcessRequestMain(Boolean includeStagesBeforeAsyncPoint, > Boolean includeStagesAfterAsyncPoint) +3412


title: Pointstack with different count-1-symbols / SLD / GeoServer

text: I try to use the pointstack sld to visualize my data at different scales. I used the example from http://suite.opengeo.org/opengeo- docs/cartography/rt/pointstacker.html, and everything is fine. However, I have a dataset of several thousand points, further defined by specific IDs. E.g. police stations (1402) or fire brigades (1401). Using the example from opengeo I can use only one symbol (or in this case an external vector graphic), but **I would like to have the specific symbol of e.g. fire brigade or police if the count is 1**. If I add a filter for the specific ID, the SLD-validation in GeoServer is fine but the filter is ignored. Therefore all points of the dataset (also non- police) are shown at the layer preview. So my idea is to set a filter at first (only use points with the ID 1401) and use the pointstack filter next. But this does not work. If I do not use the ID filter, all points of the dataset with the count = 1 are getting e.g. the police symbol. ![enter image description here](http://i.stack.imgur.com/qD5y5.jpg) The picture shows some stacked points and a lot of police stations (count = 1). However there aren´t so much police stations. Normally I would add a filter to check only points with the requested ID, so other points should be ignored for visualization (or) I can define point stacks for different IDs with different symbols. Maybe such a filter is not working in combination with the pointstack rule? I hope someone can help me :)


title: How to compare two rasters and reclass the greater value(python script)?

text: I am working on an assignment that requires us to reclass a raster(via python) from two rasters. We need to compare the power received at one cell(on raster1) and compare it to the same cell on raster 2. If the value of that cell in raster1 is greater than the value in raster2, the cell would be assigned a value of "1" while a value in raster2 greater than raster1 would get a "2". We also need to reclassify certain values as "0" since they are considered out of range. So the final raster would have three values: 0-out of range/no service 1-better value in raster1 2-better value in raster2 I can try to clarify further if needed.