title: ArcObjects : Out of Server Memory when reconciling a large edit session

text: One of our users has been periodically experiencing 'Out of Server Memory' errors with ArcSDE 10.0 when trying to reconcile a large edit session. ![ArcObjects error -2147216115](http://i.stack.imgur.com/hZW8Z.jpg) ![ArcObjects error -2147216115](http://i.stack.imgur.com/AXy9I.jpg) VMware ESX Instance: Windows Server 2008 R2 Datacenter Service Pack 1 Intel Xeon E5-2660 @ 2.20GHz 8 GB RAM I've set up a performance monitor to track memory usage but I'm curious if this is a problem anybody else has experienced with versioned editing? Our RDBMS is Oracle and I have found this page by ESRI: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//002n00000026000000 However, this seems to assume that ArcSDE and Oracle are on the same machine (for us they are not). Should I remove the 'unlimited' table space associated with the editing user? Also, we make direct connections to your database using the Oracle client. Does this mean that I'm perhaps looking at the wrong server (ArcSDE) and should be speaking with the DBAs ? I'm thinking ArcObjects error means ArcSDE problem but perhaps somebody can correct me on this? Thanks in advance, ~jess~


title: Extract pixels from NDVI in ArcMap

text: I have a NDVI that I created from a Landsat 5 image. Is there a way to extract all the pixels above a certain threshold and save it as a shapefile? I am trying to create a vegetation mask of sorts to filter out areas in a city that have vegetation, leaving only the structures. Thank you.


title: Fix Landsat 7 SLC-off in ENVI?

text: I have a couple of Landsat 7 images and they all have striping as a result of the SLC failure. Does anyone have any suggestions as to how to destripe the image using ENVI? I plan to use these images for NDVI analysis. When I downloaded the data from GloVis, it also came with a Gap Mask folder...is that somehow related to the SLC-Off? Thank you.


title: Cause of MapServer Internal Server Error: Premature end of script headers: mapserv.exe?

text: I have changed MapServer version up to 6.2, changed folders Apache/cgi-bin and Apache/php, and replaced file cgi-bin/libmap.dll with ignore/oracle/libmap.dll. When I type 'mapserv -v' in cmd it shows me 'input=oraclespatial', but when I try to open in browser localhost/cgi-bin/mapserv.exe there is 'Internal Server Error' and 'Premature end of script headers: mapserv.exe' in logs


title: How to import large features to Google Fusion Table from KML?

text: I am trying to import rather large features to Google Fusion Table from KML files (I have also shapefiles). However, Google has a certain limitation: "a cell of data in Fusion Tables supports a maximum of 1 million characters". Guess I am hitting this barrier, since conversion fails on a polygon that looks big. "Looks" is not very quantitative, so here are my question:   1. How can I find which features are hitting this limit? I need something capable of batch processing (R, Python, QGIS or ArcGIS), since I have 6000 features to check.   2. Once I find features breaking the limit, what options do I have? Should I try to split the polygons or perhaps there is a way of limiting the precision, so that the effective size is smaller?  Here are the limits explained: Google site Some more info concerning the problem, together with my data, can be found on Google product forum. Any hints will be much appreciated!


title: Plugin Ui interation with user and Qgis

text: A downsized hydrogeology plugin ( GraphicGroundwater = GGW ), instilled the need of a reduced combination of Carson's _'fTools doVectorGrid'_ and Pryet's _'Qgridder'_ plugin(s). The GGW Ui is displayed via toolbar button, but the user forms and save file selections are empty and un-editable. It seems I am missing a major communication between the Ui and MainWindow, or possibly I have implemented the dialog improperly into the GGW plugin being built. Plugin File(s) Directory (Please bear in mind, this posts author retains novice python and GIS skills.) Highly appreciate any clues or advice on how to properly connect plugin Ui to the map canvas. Thanks. ![Empty inputs](http://i.stack.imgur.com/sZVwa.png)


title: Which feature class(es) is/are used by a service

text: I have several feature classes and a lot of services. When I try to make my feature class versioned I get the message that the feature class is locked which is obvious. But it's not telling me by which service it's locked. Is there a possibility to work through all feature classes with ARCPY and find out in which service it's referenced? Thnxs for all answers!


title: how to create a cloropleth world map in arcgis 10.2

text: I am new to arcgis and I am trying to create a cloropleth map to display ratio values for different countries. I believe that my basemap will need to have a country layer so that the shading of country polygons will be possible but I don't understand how to link this with my data. Will the world topographic map have an attribute table for the country layer in which I can add a new column containing my ratio values or do I need to create and import an excel file contains this information? If I need to create an excel sheet how do I link this so that the country polygons fill to represent the ratio values? I have found some useful tutorials to show how I can view layer properties and select quantities from the symbology tab however I don't understand have to get my data to this stage. Please can someone point me in the right direction Thanks!


title: Creating a loop for Lines data (spatial data)

text: I am trying to build a lines shapefile with the contour lines from the 20 meters isobath from along the East African coast. This contour line of 20 meters has different segments and I have this information on a list, and I am trying to convert this list into a shapefile and merging the different segments of the line. How could I create a loop to join the different lines of the contour line?               library(maptools)     library(sp)          # ctlns: contour lines     ctlns     [[1]]     [[1]]$level     [1] -20     [[1]]$x     [1] 39.11667 39.10909 39.11667 39.12754 39.11667     [[1]]$y     [1] -5.327083 -5.316667 -5.303509 -5.316667 -5.327083          [[2]]     [[2]]$level     [1] -20     [[2]]$x      [1] 39.21667 39.20889 39.21667 39.21852 39.22667 39.23068 39.23333 39.24369 39.25000     [10] 39.25311 39.25173 39.25202 39.25000 39.23763 39.23333 39.22281 39.21667     [[2]]$y      [1] -5.053763 -5.050000 -5.035417 -5.033333 -5.016667 -5.000000 -4.998276 -4.983333      [9] -4.977907 -4.983333 -5.000000 -5.016667 -5.022381 -5.033333 -5.040000 -5.050000     [17] -5.053763          for(i in 1:length(ctlns))     {       ctlns[[i]]       isolines<- data.frame(ctlns[[i]])         iso <- cbind(isolines$x,isolines$y)       L1 <-Line(iso)       Ls1 = Lines(list(L1), ID = "a")       iso_spatial<-mapply(cbind, list(L1))     }      I am new using spatial data on R, and I would appreciate any help.


title: Load a vector only for a specific zoom

text: I have a map with 4-5 vectors. At the beginning, the loading of the web application is very very slow, but was fine before I've added a new one vector, which is used only for little zooms. **Is it possible to launch the loading of a vector only when the specific zoom is reached?** I'm currently using openLayers **2.2** , and as beginner in openlayers, the only way I found to load a vector was like that :               var olmap = new ol.Map({         layers: [             new ol.layer.Tile({                 source: new ol.source.OSM({                     url: 'php/map/getOsm.php?z={z}&x={x}&y={y}'                 })             }), vector, vector_country, vector_countryName, vector_term, vector_port, vector_ramp         ],         renderer: ol.RendererHint.CANVAS,         overlays: [overlay],         controls: ol.control.defaults().extend([             new ol.control.FullScreen(), new app.CustomExport()         ]),         view: mapView     });      Thank you for your help. **EDIT** That was quite easy :               olmap.addLayer(nameofthelayer);


title: Is there an open source Geocoding tool which can be used commercially?

text: Is there any Geocoding tool which is open source and can be used in commercial application? We tried Google geocoding API's but cannot continue with it because of some licensing limitation. Now we are exploring GisGraphy but haven't finalized it. Please let us know if any other better alternative is available. **EDIT:**   Currently we want to GeoCode US data only but in future we may extend it to other countries.


title: Is there a straightforward way to get the source/target point geometry of an edge?

text: I was wondering if there is a straightforward way to get the source or target point geometry of an edge? It is still not clear to me how the source-target works. For example, if I use ST_StartPoint() to my edge, will it get the source or the target of the edge? Does it matter which direction I used when I digitised the line? Also, I have set False on both 'directed' and 'reverse_cost' fields; does this mean that pgrouting will check an edge both ways? source/target and vice versa? I hope I am making sense here! Thanks in advance. A


title: Convert Features to Graphics

text: I want to convert selected features to graphics and keep the featuere symbology. I was thinking I could just call the UID command below, however it seems to not be working?                 Dim pUID As New UID         Dim pCmdItem As ICommandItem         ' Use the GUID of the Save command         pUID.Value = "esriArcMapUI.FeaturesToGraphicsCommand"         pUID.SubType = 3         pCmdItem = My.ArcMap.Application.Document.CommandBars.Find(pUID)         pCmdItem.Execute()      Why does this command not work? This code does the conversion, however it does not return the feature symbology to the graphic. What is the best method to perform this conversion while keeping the feature symbology in the new graphic? Thanks


title: ArcObjects .NET 9.3 upgrade to 10 breaks IRasterDataset.Copy(...) code

text: I have to upgrade an application from ArcObjects 9.3 (C# .NET) to 10. Unfortunately, there seems to be some incompatibility regarding the IRasterDataset.Copy method. The following code worked fine in 9.3:               my_raster_dataset.Copy("my_copied_raster", my_sde_workspace);      `my_raster_dataset` is a raster dataset called "INPUT_RASTER" loaded from a file geodatabase. `my_sde_workspace` is an ArcSDE 9.3 geodatabase running on SQL Server 2005. Using ArcObjects 10, I get the following error message: System.Runtime.InteropServices.COMException: The table was not found. [VAT_INPUT_RASTER] at ESRI.ArcGIS.Geodatabase.IRasterDataset.Copy(String copyName, IWorkspace copyWorkspace) Any ideas or suggestions? Do I have to upgrade ArcSDE to 10 to make this work? Thanks in advance. **EDIT - Found the answer** Found the problem. I was checking out `esriLicenseProductCodeEngine` instead of `esriLicenseProductCodeEngineGeoDB`. Strange, that the code worked with ArcObjects 9.3 since 2010 and suddenly stopped working with ArcObjects 10.


title: Problems with upload waypoints to GPS unit

text: How can I upload waypoints from QGIS to my Garmin GPSmap 62 by using GPS Tools in plugins? I only manage to do it with an older Garmin GPS. Regards Kaja


title: Looking for political maps (administrative division)

text: I'm looking for an interactive map that shows ONLY administrative divisions of a country. I need to personalize it so when the user clicks a region, the map shows some information about it. If it comes in colors and animations, it would be better. Is there a way I can achieve this? I've seen around and I found very good maps but maybe they show additional geographic information that is really good but might not be relevant for my web application (PHP). One more thing: I'm avoiding the use of Flash and Java Applets. Thanks.


title: Mapnik Initialization does not work

text: I create geodatabase "training"               sudo su     su postgres     createdb -E UTF8 template_postgis2     createlang -d template_postgis2 plpgsql     psql -d postgres -c "UPDATE pg_database SET datistemplate='true' WHERE datname='template_postgis2'"     psql -d template_postgis2 -f /usr/share/postgresql/9.1/contrib/postgis-2.1/postgis.sql     psql -d template_postgis2 -f /usr/share/postgresql/9.1/contrib/postgis-2.1/rtpostgis.sql     psql -d template_postgis2 -c "GRANT ALL ON geometry_columns TO PUBLIC;"     psql -d template_postgis2 -c "GRANT ALL ON geography_columns TO PUBLIC;"     psql -d template_postgis2 -c "GRANT ALL ON spatial_ref_sys TO PUBLIC;"     createdb training -T template_postgis2       ![geodatabase training](http://savepic.su/3233935.jpg) Then I changed "default.style"               sudo cp "./default.style" "/usr/share/osm2pgsql/default.style"      Then I load shapes               sudo apt-get install subversion      mkdir "/home/oleg/mapnik"     svn co "http://svn.openstreetmap.org/applications/rendering/mapnik" "/home/oleg/mapnik"     bash /home/oleg/mapnik/get-coastlines.sh       Now I want create style               python /home/oleg/mapnik/generate_xml.py osm.xml my_osm.xml --dbname training --user postgres --password postgres --accept-none      ![link](http://savepic.su/3233934.jpg) Please, tell me how create style. _**Update:_** Now I have this problem.               oleg@oleg-pkvb:~/mapnik$ ./generate_xml.py osm.xml osm_new.xml --dbname training --user osm --world_boundaries /home/oleg/world_boundaries/ --accept-none     Traceback (most recent call last):       File "./generate_xml.py", line 204, in <module>     serialize(template_xml,options)       File "./generate_xml.py", line 80, in serialize     mapnik.load_map(m,xml,True)     RuntimeError: Postgis Plugin: FATAL:  Peer authentication failed for user "osm"          Connection string: ' dbname=training user=osm connect_timeout=4'       encountered during parsing of layer 'landcover' in Layer at line 341 of 'osm.xml'      _**Update 2:_**               oleg@oleg-pkvb:~/mapnik$ ./generate_xml.py osm.xml osm_new.xml --dbname training --host localhost --port 5432 --user osm --password postgres --world_boundaries /home/oleg/world_boundaries/ --accept-none     Traceback (most recent call last):       File "./generate_xml.py", line 204, in <module>         serialize(template_xml,options)       File "./generate_xml.py", line 80, in serialize         mapnik.load_map(m,xml,True)     RuntimeError: ERROR:  relation "planet_osm_polygon" does not exist      _**Update 3_** I try download data in my database "training"               oleg@oleg-pkvb:~$ sudo osm2pgsql -U osm -d training -W postgres -H localhost -P 5432 Moscow.osm.gz     osm2pgsql SVN version 0.70.5          Password:     Using projection SRS 900913 (Spherical Mercator)     Setting up table: planet_osm_point     NOTICE:  table "planet_osm_point" does not exist, skipping     NOTICE:  table "planet_osm_point_tmp" does not exist, skipping     PREPARE get_way (int4) AS SELECT AsText(way) FROM planet_osm_point WHERE osm_id = $1;      failed: ERROR:  function astext(geometry) does not exist     LINE 1: PREPARE get_way (int4) AS SELECT AsText(way) FROM planet_osm...                                          ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.          Error occurred, cleaning up      _**Update 4:_** Now I use osm2pgsql 0.83, postgis 2.1, postgres 9.1, mapnik 300000. I create user 'gisuser' and database 'gis'. But I have the same problem ((               postgres=# ALTER USER gisuser WITH PASSWORD 'postgres';     ALTER ROLE     postgres=# \q     oleg@oleg-pkvb:~$ sudo osm2pgsql -U gisuser -W postgres -d gis -H localhost -P 5432 Moscow.osm.gz     osm2pgsql SVN version 0.83.0 (64bit id space)          Password:     Warning: program compiled against libxml 209 using older 207     Using projection SRS 900913 (Spherical Mercator)     Setting up table: planet_osm_point     NOTICE:  table "planet_osm_point" does not exist, skipping     NOTICE:  table "planet_osm_point_tmp" does not exist, skipping     SELECT AddGeometryColumn('planet_osm_point', 'way', 900913, 'POINT', 2 );      failed: ERROR:  AddGeometryColumn() - invalid SRID     CONTEXT:  SQL statement "SELECT AddGeometryColumn('','',$1,$2,$3,$4,$5, $6)"     PL/pgSQL function "addgeometrycolumn" line 5 at SQL statement


title: Flip over min/max values in a stretched symbology for raster data

text: In ArcMap 10.1, we can render a raster in a colour ramp (stretched) symbology to draw a raster layer based on each pixel value. However, the legend for the layer by default starts from the max value (top, red) down to the min value (bottom, green). How can I draw it like form min (top, green) to max (bottom, red) instead? Tried Invert checkbox but it didn't help as it was the colour that was inverted. Thanks in advance! Cheers, Alex


title: Deleting points that fall outside a polygon

text: I have a polygon that is my region of interest. I also have a database of observation points that I have imported into QGIS as a points layer. I'd like to exclude all the observation points that fall outside my region of interest (the polygon). Also happy to delete those within the polygon if that is easier (I have polygons of both region of interest and region not of interest if that makes sense). Ultimately I'm just wanting to identify those within the region and export that information (at least for now). I can do this manually by deleting each point, but is there an automated way I can do this in QGIS? I'm afraid I haven't managed to find anything to really try, what I've found relates to clipping rasters or polygons, not points. Any suggestions much appreciated - thanks


title: polyline_GetAzimuth 999999 : Error executing function

text: > **Possible Duplicate:**   >  Is there a (better) way to debug Visual Basic Code Block in calculate field > using ArcGIS? I am graduate ESCI student, and my project involves a lot of ArcGIS use. My project is related to assessing coastal erosion rates and causes near the Arctic coastlines of Barrow, I am currently trying to figure out the azimuth of a plolyline. I think I know how, but when I run it in the field calculator, I get an error, called 999999 : Error executing function. I tried several other calculations from Easy calculate just to see if they will work, but it give me the same error. I am using ArcMap 10, and I downloaded the EC50 (easy calculate). Can anyone help me out with this, and direct me on how to solve this error ? thanks


title: Populating Arcgis.com Large Thumbnail Field

text: I'm consuming content from arcgis.com to build a custom map gallery. I see that the response objects have a field called "largeThumbnail", but in my case all the values are null. Anybody know how to configure content so the largeThumbnail property is populated? Thanks in advance for any suggestions.


title: PoI Range query

text: I have a dataset of 100,000 PoIs represented by (id, name, longitude, latitude).. I would like to perform a large number (around 10,000) search operations on this dataset. each search operation should return the list of PoIs that exist within a certain distance _range_ from a given point defined by its longitude, latitude. The points locations are defined as values, where the search range distance is in Kilometers radius from a given longitude, latitude point. I tried performing an exhaustive search to find the results of search, but it takes a lot of time to finish. My dataset is stored in CSV file. I think it might be more efficient if certain structures like KD-tree or a database with Spatial Index is used for the search operation. So, what methods do you think might be most efficient for performing this function ?


title: Cleaning shapefiles with islands in GRASS GIS

text: This may be very simple, but I'm new to GRASS GIS. I have a large shapefile containing global species ranges, which needs topological cleaning. I imported this into GRASS GIS 6.4.3 using v.in.ogr (a very slow process!). When using v.in.ogr errors were encountered:               > (Mon Feb 03 14:30:37 2014)                                                 > v.in.ogr dsn=C:\GIS_Phil\IUCN_all_ranges_raw.shp     > layer=IUCN_all_ranges_raw output=IUCN_all_ranges_raw -e Projection of     > input dataset and current location appear to match Layer:     > IUCN_all_ranges_raw Counting polygons for 726 features... Importing     > map 726 features...     > ----------------------------------------------------- Building topology for vector map <IUCN_all_ranges_raw_tmp>... Registering     > primitives... 35000 primitives registered 2128395 vertices registered     > Number of nodes: 29730 Number of primitives: 35000 Number of points: 0     > Number of lines: 0 Number of boundaries: 35000 Number of centroids: 0     > Number of areas: - Number of isles: -     > ----------------------------------------------------- WARNING: Cleaning polygons, result is not guaranteed!     > ----------------------------------------------------- Break polygons:     > ----------------------------------------------------- Remove duplicates:     > ----------------------------------------------------- Break boundaries:     > ----------------------------------------------------- Remove duplicates:     > ----------------------------------------------------- Clean boundaries at nodes:     > ----------------------------------------------------- Break boundaries:     > ----------------------------------------------------- Remove duplicates:     > ----------------------------------------------------- Clean boundaries at nodes:     > ----------------------------------------------------- Break boundaries:     > ----------------------------------------------------- Remove duplicates:     > ----------------------------------------------------- Clean boundaries at nodes:     > ----------------------------------------------------- Break boundaries:     > ----------------------------------------------------- Remove duplicates:     > ----------------------------------------------------- Clean boundaries at nodes:     > ----------------------------------------------------- Merge boundaries:     > ----------------------------------------------------- Change dangles to lines:     > ----------------------------------------------------- Remove bridges:     > ----------------------------------------------------- Building topology for vector map <IUCN_all_ranges_raw_tmp>... Building areas...     > 481585 areas built 12187 isles built Attaching islands... Number of     > nodes: 568734 Number of primitives: 2323622 Number of points: 0 Number     > of lines: 0 Number of boundaries: 2323622 Number of centroids: 0     > Number of areas: 481585 Number of isles: 12187 Number of areas without     > centroid: 481585     > ----------------------------------------------------- Find centroids for layer: IUCN_all_ranges_raw     > ----------------------------------------------------- Write centroids: WARNING: 432364 areas represent more (overlapping) features, because     > polygons overlap in input layer(s). Such areas are linked to more than     > 1 row in attribute table. The number of features for those areas is     > stored as category in layer 2     > ----------------------------------------------------- 14410 input polygons Total area: 2.55953E+014 (481585 areas) Overlapping area:     > 2.17134E+014 (432364 areas) Area without category: 2.10092E+012 (18832 areas)     > ----------------------------------------------------- Building topology for vector map <IUCN_all_ranges_raw>... Registering     > primitives... 1429509 primitives registered 3434995 vertices     > registered Building areas... 481585 areas built 12187 isles built     > Attaching islands... Attaching centroids... Number of nodes: 960111     > Number of primitives: 1429509 Number of points: 0 Number of lines: 0     > Number of boundaries: 966756 Number of centroids: 462753 Number of     > areas: 481585 Number of isles: 12187 Number of areas without centroid:     > 18832     > ----------------------------------------------------- WARNING: Errors were encountered during the import Try to import again, snapping with     > at least 1e-013: 'snap=1e-013' Default region for this location     > updated Region for the current mapset updated (Tue Feb 04 08:22:58     > 2014) Command finished (1072 min 21 sec)      The next step was to export the (hopefully) cleaned file using v.out.ogr. Many of the species range polygons overlap one another, or are islands within the ranges of other species. There are also likely to be holes in the ranges of some species that are topologically represented as islands. Should I be using the -c flag when exporting this data? With the -c flag I got this message:               > (Tue Feb 04 14:36:00 2014)                                                 > v.out.ogr -c input=IUCN_all_ranges_raw@PERMANENT type=area     > dsn=C:\GIS_Phil\grassoutput\grassoutisl.shp Exporting 481585 areas     > (may take some time)... WARNING: 18832 features found without category     > were skipped v.out.ogr complete. 3723175 features written to     > <IUCN_all_ranges_raw> (ESRI_Shapefile). (Tue Feb 04 16:02:58 2014)     > Command finished (86 min 57 sec)      Without the -c flag I got this message:               > (Tue Feb 04 16:28:33 2014)                                                 > v.out.ogr input=IUCN_all_ranges_raw@PERMANENT type=area     > dsn=C:\GIS_Phil\grassoutput\grassoutni.shp WARNING: The map contains     > islands. To preserve them in the output map, use the -c flag Exporting     > 481585 areas (may take some time)... v.out.ogr complete. 3809331     > features written to <IUCN_all_ranges_raw> (ESRI_Shapefile). (Tue Feb     > 04 17:56:31 2014) Command finished (87 min 58 sec)      Are either of these actually exporting the species ranges correctly? Are the features without category that are skipped when using the -c flag holes in species ranges, in which case they should be skipped, or are they useful data (eg portions of actual species ranges) that need to be retained? Cheers, Joey


title: Select Expression and Field Calculator in ONE STEP with QGIS

text: I have an Open Street Map Shapefile and i want to use Field Calculator to add information to a field called **a800900** according to the field **TYPE** of street, by different area (field called **ID_2** ). The field TYPE of street has the following values   * motorway    * trunk   * primary   * secondary   * tertiary   * residential And the field **ID_2** is a numeric code for 85 municipalities:   * 8527   * 8300   * 8302   * ... ![enter image description here](http://i.stack.imgur.com/jILK8.png) The information added to **a800900** it's based on the **TYPE** of street and will be different per each **ID_2**. So, I've been selecting each ID_2 and then performing calculations with Field calculator in two steps. STEP 1:"Select By Expression" `ID_2 IS '8283'` STEP 2: "Field Calculator" `CASE` `WHEN TYPE IS 'motorway' THEN 2405.396*0.01847` `WHEN TYPE IS 'trunk' THEN 2305.396*0.0154` `WHEN TYPE IS 'primary' THEN 1807.226*0.013` `END` The following figure show this process: ![enter image description here](http://i.stack.imgur.com/AimJn.png) Then, i do similar calculations for another **ID_2** , varying the values in Field Calculator: STEP 1:"Select By Expression" `ID_2 IS '8302'` STEP 2: "Field Calculator" `CASE` `WHEN TYPE IS 'motorway' THEN 2405.396*0.75` `WHEN TYPE IS 'trunk' THEN 2305.396*0.34` `WHEN TYPE IS 'primary' THEN 1807.226*0.47` `END` At the end I'll get values in **a800900** different for each **TYPE** of street and **ID_2**. As example, if **ID_2** is 8300 and **TYPE** of street is motorway, then **a800900** is 2450, if **ID_2** is 8302 and **TYPE** of street is motorway, then **a800900** is 2300, and so on. I have to make the estimation with Field Calculator for each **ID_2** , being 85 areas, and thats's why i'm looking a faster procedure. I would like to know if there is a ONE STEP procedure with Field Calculator, like: `CASE` `WHEN ID_2 IS '8283'AND WHEN TYPE IS 'motorway' THEN 2405.396*0.01847` `WHEN ID_2 IS '8283'AND WHEN TYPE IS 'trunk' THEN 2305.396*0.0154` `WHEN ID_2 IS '8283'AND WHEN TYPE IS 'primary' THEN 1807.226*0.013` `WHEN ID_2 IS '8302'AND WHEN TYPE IS 'motorway' THEN 2405.396*0.75` `WHEN ID_2 IS '8302'AND WHEN TYPE IS 'trunk' THEN 2305.396*0.34` `WHEN ID_2 IS '8302'AND WHEN TYPE IS 'primary' THEN 1807.226*0.47` `END` Is there a way in one step? Thank you so much


title: Problems with Buffer tool. Looking for help please

text: I am trying to create a buffer around different series of lines. I have a polylines file with 1 location shown below: ![enter image description here](http://i.stack.imgur.com/G9cnU.jpg) I want an equal spaced buffer that is dissolved around all the lines at this location spaced 812 meters from any given end point (end points of along outside lines). I imagine this should look almost like a pillbox. However, when I run the buffer tool, this is what I get: ![enter image description here](http://i.stack.imgur.com/ayiFt.jpg) Obviously that does not look correct. It is not symmetric. When I have a large shapefile that contains many sets of lines, all of the buffers have slightly different shapes. All groups of lines have the same spacing and distances. I am not sure if this has to do with orientation because if the lines run N/S or E/W, the buffer is almost correct: ![enter image description here](http://i.stack.imgur.com/PlqZN.jpg) This is almost what it should look like, but the spacing on the sides is a lot closer than the spacing on the top. I thought it may be an issue with dissolving the group of lines together so I tried making a rectangle around the lines so it would be 1 object: ![enter image description here](http://i.stack.imgur.com/B1zPV.jpg) However, this just created the same exact incorrect buffer I provided earlier. Does anyone have any suggestions as to how I can create this buffer I am looking for? As a side note: I need this buffer around all the transects together and cannot just take the centroid of these lines and make a circular buffer.


title: No values in overlay tab, only in diagram tab

text: I have basically the same issue like in this post. I joined a table to a .shp with the intention to display those values (integers) as a column overlay. Producing pie charts (Diagram tab) works fine, the feature from the joining is included and available in the drop down menu (so i assume the NULL values are no problem). In the overlay tab however, they are not within the drop down menu. Why is the feature not available in the overlay tab? Using QGIS 1.8.0 on Win7. Thanks.


title: importing shp file in matlab

text: I have digitize 3 rivers on google earth and converted into shp file to use in matlab. I have checked the shp file in Arcgis and they are working. But while opening in Matlab, matlab shows some error. Can anyone help me out...   1. Unsupported shape type PolyLineZ (type code = 13).   2. headerTypeCode = readHeaderTypeCode(shpFileId,callingFcn); The code I am giving is below : water = shaperead('rivers', 'UseGeoCoords', true); geoshow( water , 'Facecolor', 'blue')


title: Cost Distance Paths in Python

text: I have a feature dataset with 317 points (UpperPoints) and another with 394 (LowerPoints) points. Now I want to calculate the cost paths as in this example here. However there are multiple files used and i just want to iterate trough my shapefiles. In the end I want the cost paths from every single UpperPoint to every single LowerPoint. I know that I have to calculate the cosDistance rasters for every lower Point which will need a lot of computation time, but that's not an issue. My main problem is how to extract the points from the shapefiles and iterate trough these to create the CostDistance rasters. I assume that my solution, just using the for loop is too simplistic...               PointsHIGH = "C:\\Dogon\\lcp.gdb\\Pts_rangeH_high"     PointsLOW = "C:\\Dogon\\lcp.gdb\\Pts_rangeH_low"     costsurface = ("C:\\Dogon\\lcp.gdb\\Slope_asc")          outputCostDist = "C:\\Dogon\\lcp.gdb\\CstDstLOW"     outputCostBackLink = "C:\\Dogon\\lcp.gdb\\BckLnkLOW"     outputCostPath = "C:\\Dogon\\lcp.gdb\\CostPath"          for points in PointsLOW         outputCostDist = arcpy.gp.CostDistance_sa(points, costsurface, outputCostBackLink)         outputCostPath = arcpy.gp.CostPath(points, outputCostDist, outputCostBackLink,"BEST_SINGLE")      Further I want an output of all CostDistance rasters and CostPaths, how could I do this?


title: Get the geotransform in Mercator given a pixel size away from the origin

text: I have a radar image with the following information:   1. The projection is Mercator (Mercator_1SP in GDAL)   2. The projection origin is 0º lat 0º lon   3. At a given geographic coordinates (e.g. lon=10º lat=40º), the scale is 1kmx1km each pixel. I know at wich pixel in the image this coordinates are (they are usualy at the center of the image)   4. The datum can be WGS84 or others, so it's not always the spherical earth How can I determine the geotransform? I don't know how to transform this pixel scale distance into the projection's scale at the given point. \--EDIT-- So I need the upper left pixel coordinates and the pixel width and height in the Mercator coordinates from the data explained above. Thank you.


title: Implement ArcGIS Custom Features with C#?

text: Is it possible to implement custom features in .NET? If so, how? I suspect the issue is being able to do COM aggregation with .NET, but I can't find anything on that ... Does anyone know if it is possible to create a ~~runtime~~ COM callable wrapper for a .NET class that supports COM_INTERFACE_ENTRY_AGGREGATE_BLIND ?


title: How do I use python to access the attributes of features within the current extent? (ArcGIS 10)

text: The script I'm writing requires a specific attribute of the building polygon that the data frame extent is currently focused on. Is there a way to grab this data?


I am trying to create a user, say **dev** , for cartodb running in local machine. When i run the script create_dev_user               sh script/create_dev_user dev      It produces this `RuntimeError: ERR unknown command 'HMSET'` error. I am on ubuntu 10.04.


title: turn a Result object into a number

text: **import arcpy     from arcpy import env     from arcpy.sa import *     import math       # Check out any necessary licenses     arcpy.CheckOutExtension("spatial")        # OverwriteOutput     arcpy.env.overwriteOutput = True        # Set Geoprocessing environments     arcpy.env.autoCommit = "1000"     arcpy.env.scratchWorkspace = "D:\\NEO\\GISMODELLING\\BestLocationsModel.mdb"     # arcpy.env.extent = "214386.1775 3837166.9542 256099.9616 3907138.0305"     arcpy.env.extent = "D:\\NEO\\GISMODELLING\\BestLocationsModel.mdb\WasteWaterPlant\\Tart_UTM"     arcpy.env.cellSize = "10"     arcpy.env.mask = "D:\\NEO\\GISMODELLING\\BestLocationsModel.mdb\\WasteWaterPlant\\Tart_UTM"     arcpy.env.workspace = "D:\\NEO\\GISMODELLING\\BestLocationsModel.mdb\\WasteWaterPlant"     EucDist_Aggl = "D:/NEO/GISMODELLING/BestLocationsModel.mdb/EucDist_Aggl"     # Process: Get Raster Properties     Max_Aggl = arcpy.GetRasterProperties_management(EucDist_Aggl, "MAXIMUM", "")     print Max_Aggl     delta = int (Max_Aggl)     print delta     a = delta/1000     print a**      when I running this code I have this errore > 44250.27734375 > > Traceback (most recent call last): File "D:\NEO\test.py", line 28, in delta > = int (Max_Aggl) TypeError: int() argument must be a string or a number, not > 'Result' I need to use the value of Max_Aggl in loops I try to learn Python as I go


title: Who should be responsible for editing GIS data?

text: I work for a municipality with an enterprise GIS. What would be the argument for performing all data editing in the GIS department? Conversely, what would be the argument for data updates being handled at the point of entry (e.g. DPW: sewer & water; Assessing: parcels). I am interested in the ideal case, regardless of staffing and money.


title: OSM Database created: now what happens?

text: I imported an OSM file into a Postgres DB using the Osmosis tool. Now, I need to extract some data but I can't find any schema description (I'm really a newbie in GIS). Where can I find it? Do you know if there is a tutorial to learn how to query an Osmosis Postgres DB? I need to extract the names of all streets available within a given bounding box. How can I do that? Cheers


title: How to assign new values to a raster attribute table - by comparing two fields?

text: I have created a new raster from two raster datasets using combine operation in ArcGIS. The new attribute table has two fields from the input rasters. What I need is, to compare these two fields and assign meaningful values to a new field. For ex: see the attached image (I don't know how to past a table in here) ![enter image description here](http://i.stack.imgur.com/cBVEA.jpg) ... this list goes to few hundred. “New Class” is the field that I need to assign the new values. I have been using “code block” in field calculator using python. But it is too long that I have to write individual statements for each unique combination. I am looking for an efficient way other than my method. Something like reading from a reclassification table? Thanks for your helps!


title: vector cropping failure

text: I’m trying to crop a vector contour layer using the instructions given at http://lists.osgeo.org/pipermail/qgis-user/2012-February/015294.html . I created a polygon in a newly-created vector layer, outlining the area I want to keep. (My polygon layer has the same CLS as the contour layer I want to crop.) However, when I finish I get a shapefile with nothing left in it. All of the contour lines get cropped out, inside the polygon and outside of it. Any idea what the problem might be? Thnx TVZ


title: How do I specify a CRS for a fictional game map in QGIS?

text: I am about to start playing a steampunk RPG game called The Iron Realm, which uses a fictional world map. I have recently started using QGIS v1.7 and have no real experience with GIS sysstems. I have imported the map as a raster layer into a new QGIS project and have created some vector layers from it, including a polygon layer for the regions, a point layer for the towns, and a line layer for some added railway tracks and roads. The game master has just added a scale and some lat/long information to the game map. I assume from this that I can now add a CRS to my QGIS project. I have no idea how to do this and would really appreciate some help, so that I can use my QGIS project with the game effectively. I assume that I will need to create a custom CRS, but I am completely ignorant about how to do this. I think it has something to do with proj.4 but I have no idea how to use it.


title: Map Grid Index / Reference for Canadian Topo Map

text: Does anyone know where I can find some kml at different resolutions to show the grid reference for the Canadian topographic maps (CanTopo and toporama). Much like is displayed on the map here: http://wms.ess- ws.nrcan.gc.ca/wms/indexcouverture_cantopo As you zoom in the grids get more detailed so I imagine I'd require multiple kml files to provide the different levels of detail. Thanks in advance for any suggestions.


title: When should I use ST_Buffer?

text: I would like to know what is the sense of using st_buffer operation. In PostGIS reference guide there is only the syntax and a short description of this operation. Could you give some examples to understand the advantages of this operation?


title: Free border data and C++ code to display state/region borders with user-defined-colour areas?

text: **Situation:** What I really need cannot be shortly described in the title. We have a commercial application where geographical maps (borders only) are used or for selecting areas (input) or for colouring the areas by a colour from the scale based on the value calculated for the area. The implementation is very old and it uses very bitmaps with very low resolution. The borders where just painted in Paintbrush or the like simple tool, the colouring is done via flood-fill algorithm from the reference point kept separately from the image. There is not geo-position information captured for the maps; however, we uses intevals of ZIP codes. **The first goal:** My first goal is to find the cheap solution (ideally a LGPL or similar) software to reimplement that part of the application (native C++ app, Microsoft Visual Studio 2013, i.e. C++11). I have found the OpenStreetMap exists, but I did not find the concrete answer to whether I can use (the tiny part of) it just for displaying region borders and for colorizing these polygons. **My expectations for future:** I would like to choose the solution where the polygons would be callibrated so that having a geoposition I could show where it is, to detect what polygon it is inside, etc. Later, I may want to add some simple the street layout or whatever the future will require. Can you confirm whether I the OpenStreetMap data/tools is suitable for that purpose. Are there any other free alternatives, examples, border polygons and the related code available? Thanks


title: Postgis st_intersects fails on an obtained point on the line with the line itself

text: I am trying to split a line from its intersection with another line. I know I can use `st_split` directly here but I wish to split a line from several positions where it intersects other lines. What I discover is the following: I know that lines 1 and 31 intersect                   select st_intersects(a.geom,b.geom) from dumped_roads as a,                  dumped_roads as b where a.gid=1 and b.gid=31;      Output:                   st_intersects          ---------------         t         (1 row)      Also, `st_intersection(a.geom,b.geom)` returns an intersection point. But when I do,                   select st_intersects(a.geom, st_intersection(a.geom,b.geom)) from             dumped_roads as a,dumped_roads as b where a.gid=1 and b.gid=31;              st_intersects           ---------------         f         (1 row)      I cant understand whats the problem. Although `a.geom && st_intersection(a.geom,b.geom)` returns true but that wouldnt help me in splitting the line with the point obtained from intersection. If its indeed how it is, any way out on breaking a line from an obtained point which physically lies on the line but is not in the linestring's point set?


title: How to auto-fill a new point vector in a vector layer?

text: When I work _in-situ_ I manually add points in my shape-file. I would like that when I add a new point some attributes would already filled with values calculated programmatically. E.g in my first column a progressive number (first point automatically add `1`, second `2` etc. etc.) like an auto `$rownum` second column **id of the instrument measure** usually is something like `$rownum+n`where `n`is a number. I'm able to to that only in post processing with the command I've mentioned but for me would be very helpful have a quick look to labeled point when I'm taking them.So is there a way to pre-fill determinate values of a new point vector?


title: How can I list version name in my silverlight application?

text: I need to populate combobox with arcsde version name in my silverlight 5 application. Please give me some suggestions.


title: ET Geowizards: Table Already Exists Log Error

text: I have a license for ET Geowizards and I keep coming up with an error using various tools. Most recently I have been using the Global Snap Polylines tool. I ran it a few times and then it randomly threw the follwoing error in the log file: !Start Snap Polyline Layer Wizard Date: 2/27/2013 Time: 3:37:44 PM **The table already exists.** ETGW 10.2 for 10.1, 21 July 2012; ArcGIS 10.1.3035: ET_ExportToPGDBFCDropZM w =4 !Error: Snap Polyline Layer Wizard Date: 2/27/2013 Time: 3:37:44 PM The table does not already exist. I've tried using different output names and even different output workspace locations and the error still persists. I'm wondering if I need to clear the ET Geowizards cache somehow or maybe it's another issue altogether? If anyone knows what is causing this, I would like to hear. Thanks, Mike


title: Proxy trouble with GeoExt Query Panel?

text: I got the GeoExt QueryPanel working so far, it recognizes the layers available on my Geoserver but it is not able to get the attribute columns for the dataset. I configured an osgeo proxy on my Apache Server, I use the override- ext-ajax.js and have the OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url="; class configured. A DescribeFeatureType request is sent when selecting the layer I want to filter. The same happens when trying to choose the attribute column in the query panel, but there are none displayed. The sent request URLs look like this:               http://localhost/cgi-bin/proxy.cgi?url=     http%3A%2F%2Flocalhost%2Fcgibin%2Fproxy.cgi%2Fwfs%3F     VERSION%3D1.1.0%26REQUEST%3DDescribeFeatureType%26TYPENAME%3Dcases%26query%3D      The strange part here is, that the request is sent via the proxy.cgi twice. The request URL should look like this, going passed geoserver´s WFS:               http://localhost/cgi-bin/proxy.cgi?url=     http%3A%2F%2Flocalhost%3A8080%2Fgeoserver%2Fwfs%3F     VERSION%3D1.1.0%26REQUEST%3DDescribeFeatureType%26TYPENAME%3Dcases%26query%3D      Does anybody have a clue what´s going wrong here and why the request URL is wrong? The strange thing is that even if I remove the Open.Layers.ProxyHost class and the override-ext-ajax.js from my application, the request is sent that way. Can it be related to my server? Other requests sent in the application (GET & POST) work fine. The (relevant) code looks like this, maybe you guys can spot the error...               OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";              var owsUrl = "/geoserver/wfs";              var symbolizer = OpenLayers.Feature.Vector.style["default"];              symbolizer.graphicName = "square";              var store;              Ext.onReady(function() {                          store = new GeoExt.data.WFSCapabilitiesStore({                     url: Ext.urlAppend(owsUrl,                        "VERSION=1.0.0&REQUEST=GetCapabilities"),                     protocolOptions: {version: "1.1.0"},                     autoLoad: true,                     listeners: {                         load: initApp                     }                     });                 });                 function initApp() {                 var vecLayer = new OpenLayers.Layer.Vector("Query", {                     eventListeners: {                         featuresadded: function() {                             var extent = vecLayer.getDataExtent()                             mapPanel.map.zoomToExtent(extent)                         }                     }                 });                       var queryPanel = new gxp.QueryPanel({                     title: "Query",                     region: "center",                     width: 350,                     bodyStyle: "padding: 10px",                     autoScroll: true,                     layerStore: store,                     bbar: ["->", {                         text: "Query",                         handler: function() {                             queryPanel.query();                         }                     }],                     listeners: {                         query: function(panel, store) {                             var rule = new OpenLayers.Rule({                                 filter: panel.getFilter()                             });                             rule.symbolizer[queryPanel.symbolizerType] = symbolizer;                             var style = new OpenLayers.Style("", {rules: [rule]});                             wms.mergeNewParams({                                 sld_body: new OpenLayers.Format.SLD().write({                                     namedLayers: [{                                         name: panel.selectedLayer.get("name"),                                         userStyles: [style]                                     }]                                 })                             });                                  store.bind(vecLayer, {                                 initDir: GeoExt.data.FeatureStore.STORE_TO_LAYER                             });                         }                     }                 });      Any help is appreciated! Cheers, Philippe


title: Are Online Certs worth their weight?

text: No GIS Certification Programs in my area...Are Online Certs worth their weight? Background Information: Currently a Computer Science student, with 3 semesters remaining. My college doesn't offer any GIS courses, and there aren't any local certificate programs available. Would it be worth trying an Online Certification Program? Or should I try to teach myself as much information as possible and try to get into a position after graduation without any formal GIS info on my resume?


title: Is it possible to maintain primary keys on versioned featureclasses without arcobjects code?

text: Since it is not a good idea to use ObjectIDs as primary keys in relationships, some other column will need to be used. Is it possible to use DBMS sequences to populate primary key columns in versioned featureclasses, or will I need to write arcobjects code to do this? I've heard Global IDs can be used, but non-GIS DBA's maintaining related tables often frown at this idea.


title: Versioning commercial datasets in ArcSDE + Oracle

text: We are building an application with several spatial tables, provided by a commercial vendor, at its core. The vendor provides updates quarterly. We maintain several tables of our own, derived from these commercial tables. Our application would benefit from versioning the commercial tables along with our derived tables, such that we can see a consistent view of things at any point in time. Upon each update from the vendor, I would like to "diff" the changes and apply them to our database, then update our derived tables accordingly. Can this be accomplished with ESRI's implementation of versioning? Would this require registering our tables with the GDB? (I prefer to do the heavy lifting with PL/SQL, so I'm not sure how I feel about this.)


title: Using the GEOTIFF Model Transform Tag

text: I have managed to figure out how to map Lat/Long to pixels in FAA charts that use Tie Point tags. However, I am now working with charts that use Model Transform tags. So far I am at a loss on how to do Pixel=>Lat/Long and Lat/Long=>pixel transformations. The GEOTIF spec appears to say that Lat/Long Vector = Transform-Matrix * Pixel-Vector Here is an example transformation matrix:               798.320542 0.00959398695 0 -6084843.93      0.0203216712 -798.350061 0 3716484.45      0 0 0 0      0 0 0 1      If I have a transposed [X, Y, Z, 1] position vector, the -6084843.93 and 3716484.45 values in the matrix send me way out in left field. Here is an example for another chart:               92.4530932 5.28522981 0 455424.367      5.28258406 -92.4582959 0 554879.032     0 0 0 0      0 0 0 1      Obviously, I am missing something. I have all the My questions are 1) Given a vector of pixel coordinates [x, y, z, 1] and a Model Transform matrix (or its inverse), what is the process to determine Lat/Long. 2) Given Lat/Long coordinates and a Model Transform matrix (or its inverse), what its he process for determining pixel coordinates. Thanks


title: Is there any way to use open source tools to make sure that lines on a map dont overlap?

text: I am a generalist developer getting the hang of basic gis. I made this geojson map by clipping and converting different shapefiles. The problem is that the precincts of New Orleans extend into the Mississippi river. I want them to stop at the river. Is there any easy tool or method to get the shapes from not overlapping? The borders of the precincts should end at the river. The precincts should not extend into the river. I am using all open source tools.


title: How do I download specific OpenStreetMap data by tag?

text: I'd like to write an iPhone app which uses OSM data to map canals. The data is tagged in OSM as canals, but I'm not entirely sure how to download it. Libraries like route-me don't seem to provide a way to download the data. Ideally I'd like to give the user the option to download certain counties, then scrape certain OSM data tag types. How would I go about this?


title: Does QGIS support .sld style files for raster layers?

text: I'm currently looking at the documentation for geoserver, and I'm trying to load a raster and .sld file in QGIS to see if I can edit .sld files for raster layers. The .sld and raster are located here: http://docs.geoserver.org/2.4.1/user/styling/sld-cookbook/rasters.html Unfortunately I'm unable to load the .sld file in QGIS 2.2. I am only allowed to load .qml files when I try to load a style for that layer. It appears that QGIS only loads .sld files for vector layers. Is there anyway to use an .sld file on a raster layer?


title: (Professional) Label Positioning?

text: Just a general question here. I've implemented labeling in my code and solved a fair amount of problems that seemed insurmountable at first. Now that the labeling is working I have directed my attention to the positioning of the labels and am extremely frustrated with the capabilities offered by the base ArcEngine. I have no experience with Maplex (or any other tool) and suspect that that extension (or a 3rd party tool) is the only way to do decent labeling in the ArcGIS product line. Is this a true picture of the general experience? If so, what solution would you, dear reader, offer? Is Maplex the best and only way? Are there acceptable 3rd party solutions? Any suggestions in this area will be most appreciated!


title: problem georeferencing a TIFF QGIS

text: I'm trying to georeference a TIFF in QGIS2.0, 32Bit. The problem appears when QGis adds the map in my project, the map transforms into a black square and it goes to the wrong place. In the next image you can see the configurations of the transformation and the final result. in the other image you can see the where it is without georeferencing and where i want to place it. ![enter image description here](http://i.stack.imgur.com/iuo2r.jpg) Do you know why this is happening? Or other alternative way to do this?


title: How can I change all the entries in a field to be entirely lower case?

text: I have a field with a bunch of species names, but for a single species there a bunch of case combinations (e.g. Cat, CAT,cat). I want to make a new field which has all the entries in all lower case. I tried using the field calculator with the script = !SPECIES!.lower() but keep getting **ERROR 000539**. Anyone have any idea where I am going wrong?


title: Own layerswitcher value not defined

text:
I am trying to write my own layerswitcher because I want a dropdown box (and later a second one) My Webinspector says: `Uncaught TypeError: Cannot read property 'value' of undefined` but I cannot find what is wrong. In my opinion it is defined. Thanks to everyone trying to help me. Here is the code:

     <select onchange="auswahlLayer1()">             <option value="Einwohner">Einwohner</option>             <option value="Energieverbrauch">Energieverbrauch</option>     </select>


title: How to split polygons with another polygon and inheritage data?

text: Let's say that we have two tables with polygons representing areas. Table A and Table B. Geometries from table A represents both the valid and invalid areas. And the geometries in table B represents the invalid ones. For example, lets say that the geometries from table A represent, harvesting fields, and table B represents roads, shacks and generally areas that the harvesting is impossible for w/e reason (as polygons). Table B represents some kind of mask that's has being digitized differently. The problem arises that the mask, can pass through many fields. So I want to clip that mask (Table B), with the field polygon (Table A) to create a new set (Table C) consisting with the geometries that represents the masks 'per field' along with the name of the field it was cut from. This can be done by `create tableC as (select st_intersection(a.geom,b.geom) geom from tableA as a, tableB as b);` But I have troubles on how to select the tableA.name as well and transferring it to tableC. the selection part of an ideal querry would be like this: `create tableC as (select st_intersection(a.geom,g.geom) geom, tableA.name as name from .... );` A sample defintion of table C is the following:               Column     |   Type      |      gid        | int         |     geom       | geometry    |  <- st_intersetion(a.geom,b.geom)     name       |varchar(255) |  <- name of the feature of table A which took part in the st_intersection       Table A is defined as following:               Column     |   Type      |      gid        | int         |     name       | varchar(255)|     geom       | geometry    |      And Table B as follows:               Column     |   Type      |      gid        | int         |     geom       | geometry    |      And for a visual representation : ![enter image description here](http://i.stack.imgur.com/3gfWS.png)


title: Can I Open and edit a point shape file from another point shape file

text: My question is, If out in the field gathering points in one shape file with its own attribute table is it possible to create a form that will allow you to gather a point with a different attribute table ex.. If I am gathering points on a pipeline which has its own attributes and I come across a valve which has its own attributes can I open the valve layer for editing from the Pipeline layer? I have eleven different shape files but they are all associated with pipeline I would like to make a form that lets me fill in the attributes for a pipeline then click a button for the valve layer or any of the other layers then go back to pipeline


title: no Get More in Plugin Manager in qgis 2 + ubuntu 13.04

text: After upgrade of qgis 1.8 to qgis 2.0 I couldn't find the choice "Get more" in Plugin Manager and I can't install plugins. ![enter image description here](http://i.stack.imgur.com/D1ApV.png) ![enter image description here](http://i.stack.imgur.com/KZET9.png)


title: Shapefiles from Nationalgrid without a .prj file - what projection are they in?

text: Im working with these shapefiles: http://www2.nationalgrid.com/uk/services/land-and-development/planning- authority/shape-files/ They show pipelines & electricity grid of the UK. They do not have a .prj file included in the shapefile. They display fine in QGIS when layered on top of each other, but I need to know which projection they are in so I can display them on a world map. The x and y values are very large so I tried the Universal transverse Mercator projection "WGS84/UTM zone 30N ESPG:32630" because the UK is in zone 30.   This placed the shapefiles at the right scale, but just above the equator by Ghana (48 degrees too far south and 3 degrees too far west). ![enter image description here](http://i.stack.imgur.com/EJOlP.jpg) This seems to imply that the shapefiles' coordinates are zeroed half way between UTM zone 30U and 31U.   http://en.wikipedia.org/wiki/File:Utm-zones.jpg Does anyone know which projection I should be using? Thanks for your help.


title: Shortest Path Dijkstra

text: I would like to calculate the shortest path with dijkstra algorithm. sample query `SELECT seq, id1 AS node, id2 AS edge, cost FROM pgr_dijkstra(' SELECT gid AS id, source::integer, target::integer, length::double precision AS cost FROM ways', 30, 60, false, false);` I imported an osm map, it has proper network. The name of the table is `at_2po_4pgr`, colums are `id, osm_id, osm_name, osm_meta, osm_source_id, osm_target_id, clazz, flags, source, target, km, kmh, cost, reverse_cost, x1, y1, x2, y2, geom_way`. in result I would like to see the shortest path and streets name in a row where route was planning. Thanks


title: How to extract one band from a three band raster in ArcMap?

text: I'm fairly new to ArcGIS and I get into trouble extracting a single band layer from a raster image .img in ArcMap that I want to use after in zonal statistics. So far, all help I found is to do the operation in the other way around and merge bands and some explanation by using the python code. But I would like to use ArcMap interface. I tried to create a single layer from ArcCatalog and import it to ArcMap but I get a layer and not a raster. I also tried to import my image.img in ArcMap but then I don't find anything to extract only a specific band. If somebody could give me a clue, it would be very much appreciated.


title: Measuring area of Kernel Density plot?

text: Are there any tools or methods in ArcMap that can measure the area of each level of the kernel density output? Note this is just a kernel density plot (not a shape file or polygon). Kernel Density Analysis (dark green (10%) to red (90%) (1) I have tried a few things but i need it to be accurate. ![Kerenel Density Analysis \(dark green \(10%\) to red \(90%\)](http://i.stack.imgur.com/n9CsW.jpg)


title: Bug in ArcGIS WebADF's IMapFunctionality.DrawExtent?

text: I'm maintaining an old ArcGIS WebADF-based application that has recently been upgraded to use the 10.1 libraries. I'm not too optimistic about finding a solution as I think I'm the only person in the world still using this framework. However... There seems to be a bug in the function               IMapFunctionality.DrawExtent      whereby a graphics layers will not draw at the requested extent. It seems to draw at whatever the map's current extent is, not the desired extent that is passed in, and this may not be consistent. Non-graphics layers draw at the requested extent no problem. All online examples I've seen that use this function simply ask for the map's current extent to be rendered so I'm not too surprised that this hasn't come up before. The application is generating a PDF that combines the base layer and the graphics layer, which includes user-redlined features, in one image. The image has to be drawn at a specific extent which is not always the same as the map's extent. I thought there was a workaround available in asking the Map control to use web-tier image blending mode, so that the 'base layer' IMapFunctionality would include the graphics, but this isn't helping - the image only shows the base layer. Anyway, if anyone has any experience with this problem, or any insight on how to work around it, I would be very appreciative. Cross-posted here


title: Can .Net version of ArcGIS Server expose Oracle spatial data without ArcSDE, or is this limited to Java version?

text: I'm planning a new deployment of ArcGIS Server and I want to avoid using ArcSDE if possible. My vendor requires that we store our spatial data in Oracle and we prefer to use Microsoft server technologies whenever possible. Our ideal environment is to host ArcGIS Server on Windows Server and expose the REST services over IIS. Can the .Net version of ArcGIS Server expose Oracle spatial data without ArcSDE or is this limited to the Java version?


title: Calculate Latitude and Longitude from Range, Azimuth, and Elevation

text: I was wondering if there is a way to calculate the Latitude and Longitude of a given point provided you know the range, azimuth, and elevation of said point?


title: MapProxy installation

text: I'm tryign to install MapProxy on a windows manchine. My hope is that I can install it and test it on my own computer before moving it to a server. At the moment I can seem to find mapproxy-util and to be honest I don't understand how I am supposed to call it. I have read the tutorial at mapproxy.org but if anyone know of some other online tutorials that have been helpful I would really like to hear about it. Any advise is welcome.


title: Using Tilestache for serving maps to iPAD

text: I am planning to use tilestache as the tile server (map server) for serving maps to both iPAD and web clients. I plan to generate tiles by Tilemills and store tiles in a Sqlite database (mbtiles) and serve it using Tilestache. **My questions are:** **(A)** How scalable is the solution, for example the project may have 200 IPAD users in field? **(B)** Is client side local caching supported or do we have to pre-load tiles in the local storage? **(C)** Most important, can any such mobility solution using tilestache be sited to see the performance results?


title: How to create a data package for arcgis 9.3 in arcgis 10?

text: I need to create a data package in ArcGIS 10 for a company that is using ArcGIS 9.3. Is there any way to downgrade the versions like GDBs. I'm grateful if someone has a solution for this and help me. Thanks.


title: Heat map which will be auto-updated from Google Docs

text: I want to develop a heat map (standardized choropleth map) which will be auto- updated with data pulled from GoogleDocs. Which tools are the most helpful for developing that kind of map? Thank you.


title: "Dataset xxx does not exist or is not supported" using Table to Domain in Python

text: I have a directory of hundreds of *.csv files, containing code/description pairs for geodatabase domains. I'm trying to automate the creation of these domains using Python.   1. If I manually run the Table to Domain tool using a *.csv file in ArcToolbox, the domain is correctly created   2. If I create a Model using the Table to Domain tool, hard-coding in the input file (and all other parameters) then export this to a Python script, the domain is correctly created when I run the Python script   3. If I try to iterate through all *.csv files in the directory, running the Table to Domain tool for each file, I get the ""Dataset xxx does not exist or is not supported" error message: > Failed to execute. Parameters are not valid.ERROR 000732: Input Table: > Dataset C:\temp\01-01.csv does not exist or is not supported ERROR 000308: > Invalid field type Failed to execute (TableToDomain). I've placed a package here containing the template geodatabase, the working and broken python scripts, a few of the CSVs and the model which was used to export the working python script. I suspect the problem is due to the path not being correct for the *.csv file, but I can't work out how to get it working. `os.path.exists` shows that the file does exist, so it may be an internal ArcGIS thing. Note that the *.csv file does work correctly via ArcToolbox and the hard-coded Python script, so it's unlikely that the field type really is invalid. Any help gratefully received. Steve


title: Can QGIS wrap labels on second space?

text: The answer to this question tells us how to wrap labels on spaces. Is it possible to wrap labels on defined spaces? For example, wrap on just the first space or on just the second space. Secondly, is there a way of controlling the amount of vertical space inserted when wrapping labels? Wrapping on a character (or on a control-character) inserts too much vertical space. Nick.


title: Is it possible to have a customized zoom level?

text: I am setting up my own Map Tile Server according to this guide and it is working well now: http://switch2osm.org/serving-tiles/building-a-tile-server- from-packages/ The map currently has 19 zoom levels and it is initialized in the OpenLayer object "map" as below code shows:               map = new OpenLayers.Map ("map", {                 controls:[                     new OpenLayers.Control.Navigation(),                     new OpenLayers.Control.PanZoomBar(),                     new OpenLayers.Control.Permalink(),                     new OpenLayers.Control.ScaleLine({geodesic: true}),                     new OpenLayers.Control.Permalink('permalink'),                     new OpenLayers.Control.MousePosition(),                                         new OpenLayers.Control.Attribution()],                 maxExtent: new OpenLayers.Bounds(-20037508.34,-20037508.34,20037508.34,20037508.34),                 maxResolution: 156543.0339,                 numZoomLevels: 19,                 units: 'm',                 projection: new OpenLayers.Projection("EPSG:900913"),                 displayProjection: new OpenLayers.Projection("EPSG:4326")             } );      I am wondering if I want to zoom beyond the maximum zoom level 19, say zoom into level 23 or even better, have my own zoom level reconfigured, how to do that? This is because I want to display an indoor room map on my current map server. The indoor room map is small it needs zoom level higher than 19. Please point me to some resources to achieve this task as well. Thank you.


title: Stripes in Landsat-image after topographic correction (SAGA)

text: I got stripes in my topographic corrected Landsat images (see second image below). I used SAGA GIS 2.0.8 at xubuntu 12.04 Image and DEM are in the same projection, tried it with EPSG:32636 and EPSG:4326 ![original Landsat image](http://i.stack.imgur.com/obdPf.png) ![topo corrected Landsat image](http://i.stack.imgur.com/wmdt8.png)


title: lonlat from geometry point gives strange result

text: I added a vector layer to my main map with the controls like the example from openlayers: http://openlayers.org/dev/examples/draw-feature.html Now when I want to get the resulting lonlat from the point that was drawn I get strange results:               ({id:"OpenLayers.Geometry.Point_81", x:6418264.3901563, y:7161906.3780079,            bounds:{left:6418264.3901563, bottom:7161906.3780079, right:6418264.3901563,                   top:7161906.3780079}})      The main map is constructed like this:               map = new OpenLayers.Map( $('map'), {         sphericalMercator: true,         //maxResolution: '360/512',         maxExtent: new OpenLayers.Bounds(-20037508.34,-20037508.34,20037508.34,20037508.34),         projection: "EPSG:900913",         displayProjection: new OpenLayers.Projection("EPSG:4326"),         tileOrigin: new OpenLayers.LonLat(-180, -90),         controls: [new OpenLayers.Control.Navigation({dragPanOptions: {enableKinetic: true}}),                        new OpenLayers.Control.MousePosition(),new OpenLayers.Control.Zoom()],     eventListeners:{"moveend":mapEvent}                                              });      Can somebody give me a hint, it's probably an easy thing but I'm quite new to openlayers and javascript in general :-) Something with the conversion of projection ?


title: Dissolve ploygons in QGIS or Grass to new polygons of defined population sizes

text: I have a shapefile with 3800 polygons representing the smallest administrative areas in Ireland. Many of these, in rural areas, have tiny populations, under 100 people. In urban areas, all would have adequate populations. For each area I have a population figure, and various other attributes, all of which can be intelligently summed to give valid results for areas produced by combining the smaller areas. I need to merge the areas to get say, 500 or 600 areas of reasonable population. I need to leave the larger (in population) areas alone, where possible. I would greatly prefer to do this automatically. This is necessary for model fitting purposes, primarily. I've looked at a number of the dissolve questions here. 'Batch Dissolve Polygons based on Attributes in QGIS or FWTools' I think would require me to know which areas to merge, when this is what I need to calculate. 'Practical way of managing polygons and unions of them' and 'How do dissolve polygons from shapefile using open source tools?' have the same issue. Any suggestions?


title: How to read each pixel of each band of a multiband GeoTiff with GeoTools Java

text: I am new to working with Raster in GeoTools-very experienced in java. My goal is to read in a geotiff that has multiple bands, then create a HashMap for each pixel, where the key is simply the lat/lon of the pixel, and the value is a sortedMap in which teh key is the band name, and the value is the real value of the pixel for that band. I am then going to do some fascinating stuff with that map. So in basic terms, I need to read in a tif, get each band, and access then pixel's lat lon in each band. I have a rather pathetic attempt to read in a tiff below, but I can't find any straightforward way online to get each band, and get each pixel. Any help at all would be greatly appreciated.               GridCoverage2D g;// = new GridCoverage2D(null, null);                GeoTiffReader r = new GeoTiffReader(new File("c;\\temp\\thetiff.tif"));           GridEnvelope originalGridRange = r.getOriginalGridRange();      UPDATE: So I can now read in the file with the following code:                    public static void main(String[] args) {             try {               AbstractGridCoverage2DReader reader = new GeoTiffReader(new FileInputStream(new File("C:\\my.TIF")));//...; // Creating a reader               String[] coverageNames = reader.getGridCoverageNames();               // At this point, coverageNames may contain, as an instance, "pressure,temperature,humidity"                    String requestedCoverageName = coverageNames[0]; // e..g, "temperature"                    // Getting the coverage's properties               final GeneralEnvelope envelope = reader.getOriginalEnvelope();               final GridEnvelope gridRange = reader.getOriginalGridRange();                    // reading the coverage     //FAILS HERE               GridCoverage2D coverage = (GridCoverage2D) reader.read(requestedCoverageName, null);               System.out.println("test");             } catch (Exception ex) {               Logger.getLogger(GeoTiffvectorizer.class.getName()).log(Level.SEVERE, null, ex);             }           }      problem now is I know the image has 8 bands, but I only got one item in coverageNames[] and I get the following stacktrace error Error: One factory fails for the operation "ImageRead" Occurs in: javax.media.jai.ThreadSafeOperationRegistry java.lang.reflect.InvocationTargetException at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method) at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:62) at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43) at java.lang.reflect.Method.invoke(Method.java:483) at javax.media.jai.FactoryCache.invoke(FactoryCache.java:122) at javax.media.jai.OperationRegistry.invokeFactory(OperationRegistry.java:1674) at javax.media.jai.ThreadSafeOperationRegistry.invokeFactory(ThreadSafeOperationRegistry.java:473) at javax.media.jai.registry.RIFRegistry.create(RIFRegistry.java:332) at javax.media.jai.RenderedOp.createInstance(RenderedOp.java:819) at javax.media.jai.RenderedOp.createRendering(RenderedOp.java:867) at javax.media.jai.RenderedOp.getWidth(RenderedOp.java:2179) at org.geotools.gce.geotiff.GeoTiffReader.read(GeoTiffReader.java:560) at org.geotools.coverage.grid.io.AbstractGridCoverage2DReader.read(AbstractGridCoverage2DReader.java:239) at dgi.vectorspaceclustering.GeoTiffvectorizer.main(GeoTiffvectorizer.java:38) Caused by: java.lang.RuntimeException: javax.imageio.IIOException: I/O error reading header! at com.sun.media.jai.imageioimpl.ImageReadCRIF.create(ImageReadCRIF.java:317) ... 14 more Caused by: javax.imageio.IIOException: I/O error reading header! at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.readHeader(TIFFImageReader.java:458) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.locateImage(TIFFImageReader.java:465) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.seekToImage(TIFFImageReader.java:589) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.getImageTypes(TIFFImageReader.java:1104) at javax.imageio.ImageReader.getRawImageType(ImageReader.java:681) at com.sun.media.jai.imageioimpl.ImageReadOpImage.layoutHelper(ImageReadOpImage.java:228) at com.sun.media.jai.imageioimpl.ImageReadOpImage.(ImageReadOpImage.java:473) at com.sun.media.jai.imageioimpl.ImageReadCRIF.create(ImageReadCRIF.java:309) ... 14 more Caused by: java.io.EOFException at javax.imageio.stream.ImageInputStreamImpl.readShort(ImageInputStreamImpl.java:229) at javax.imageio.stream.ImageInputStreamImpl.readUnsignedShort(ImageInputStreamImpl.java:242) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.readHeader(TIFFImageReader.java:404) ... 21 more 2014-07-14T16:30:36.845-0400 SEVERE null javax.media.jai.util.ImagingException: All factories fail for the operation "ImageRead" at javax.media.jai.OperationRegistry.invokeFactory(OperationRegistry.java:1687) at javax.media.jai.ThreadSafeOperationRegistry.invokeFactory(ThreadSafeOperationRegistry.java:473) at javax.media.jai.registry.RIFRegistry.create(RIFRegistry.java:332) at javax.media.jai.RenderedOp.createInstance(RenderedOp.java:819) at javax.media.jai.RenderedOp.createRendering(RenderedOp.java:867) at javax.media.jai.RenderedOp.getWidth(RenderedOp.java:2179) at org.geotools.gce.geotiff.GeoTiffReader.read(GeoTiffReader.java:560) at org.geotools.coverage.grid.io.AbstractGridCoverage2DReader.read(AbstractGridCoverage2DReader.java:239) at dgi.vectorspaceclustering.GeoTiffvectorizer.main(GeoTiffvectorizer.java:38) Caused by: java.lang.reflect.InvocationTargetException at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method) at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:62) at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43) at java.lang.reflect.Method.invoke(Method.java:483) at javax.media.jai.FactoryCache.invoke(FactoryCache.java:122) at javax.media.jai.OperationRegistry.invokeFactory(OperationRegistry.java:1674) ... 8 more Caused by: java.lang.RuntimeException: javax.imageio.IIOException: I/O error reading header! at com.sun.media.jai.imageioimpl.ImageReadCRIF.create(ImageReadCRIF.java:317) ... 14 more Caused by: javax.imageio.IIOException: I/O error reading header! at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.readHeader(TIFFImageReader.java:458) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.locateImage(TIFFImageReader.java:465) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.seekToImage(TIFFImageReader.java:589) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.getImageTypes(TIFFImageReader.java:1104) at javax.imageio.ImageReader.getRawImageType(ImageReader.java:681) at com.sun.media.jai.imageioimpl.ImageReadOpImage.layoutHelper(ImageReadOpImage.java:228) at com.sun.media.jai.imageioimpl.ImageReadOpImage.(ImageReadOpImage.java:473) at com.sun.media.jai.imageioimpl.ImageReadCRIF.create(ImageReadCRIF.java:309) ... 14 more Caused by: java.io.EOFException at javax.imageio.stream.ImageInputStreamImpl.readShort(ImageInputStreamImpl.java:229) at javax.imageio.stream.ImageInputStreamImpl.readUnsignedShort(ImageInputStreamImpl.java:242) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.readHeader(TIFFImageReader.java:404) ... 21 more Caused by: java.lang.reflect.InvocationTargetException at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method) at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:62) at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43) at java.lang.reflect.Method.invoke(Method.java:483) at javax.media.jai.FactoryCache.invoke(FactoryCache.java:122) at javax.media.jai.OperationRegistry.invokeFactory(OperationRegistry.java:1674) at javax.media.jai.ThreadSafeOperationRegistry.invokeFactory(ThreadSafeOperationRegistry.java:473) at javax.media.jai.registry.RIFRegistry.create(RIFRegistry.java:332) at javax.media.jai.RenderedOp.createInstance(RenderedOp.java:819) at javax.media.jai.RenderedOp.createRendering(RenderedOp.java:867) at javax.media.jai.RenderedOp.getWidth(RenderedOp.java:2179) at org.geotools.gce.geotiff.GeoTiffReader.read(GeoTiffReader.java:560) at org.geotools.coverage.grid.io.AbstractGridCoverage2DReader.read(AbstractGridCoverage2DReader.java:239) at dgi.vectorspaceclustering.GeoTiffvectorizer.main(GeoTiffvectorizer.java:38) Caused by: java.lang.RuntimeException: javax.imageio.IIOException: I/O error reading header! at com.sun.media.jai.imageioimpl.ImageReadCRIF.create(ImageReadCRIF.java:317) ... 14 more Caused by: javax.imageio.IIOException: I/O error reading header! at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.readHeader(TIFFImageReader.java:458) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.locateImage(TIFFImageReader.java:465) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.seekToImage(TIFFImageReader.java:589) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.getImageTypes(TIFFImageReader.java:1104) at javax.imageio.ImageReader.getRawImageType(ImageReader.java:681) at com.sun.media.jai.imageioimpl.ImageReadOpImage.layoutHelper(ImageReadOpImage.java:228) at com.sun.media.jai.imageioimpl.ImageReadOpImage.(ImageReadOpImage.java:473) at com.sun.media.jai.imageioimpl.ImageReadCRIF.create(ImageReadCRIF.java:309) ... 14 more Caused by: java.io.EOFException at javax.imageio.stream.ImageInputStreamImpl.readShort(ImageInputStreamImpl.java:229) at javax.imageio.stream.ImageInputStreamImpl.readUnsignedShort(ImageInputStreamImpl.java:242) at it.geosolutions.imageioimpl.plugins.tiff.TIFFImageReader.readHeader(TIFFImageReader.java:404) ... 21 more


title: Web mapping with Openlayers

text: On a pedestrian wayfinding scenario, I need to drape some extruded buildings close to the road as the user navigates along the road on a 2D Map Background. I hope I can use OpenLayers and ReadyMap for this. I have vector data in shape file format. I am a newbie to web-mapping. Any ideas on this to start with?


title: Using AndNav tiles storage as input for OpenLayers.TMS

text: Is there a possibility to use map tiles from AndNav (created wit Mobile Atlas Creator) that are in format: Mapnik\12\2241\1369.png.andnav (I know that 12 is the zoom level but don't have idea what is 2241 and 1369...) as an input for OpenLayers.TMS layer? I found information that I could use local tiles with TMS: http://trac.osgeo.org/openlayers/wiki/UsingCustomTiles Unfortunatelly, I can find any place where AndNav map format is documented and reading sources didn't helped me too much. Would anyone could help me?


title: OpenLayers, zoom to vector layer extent

text: It seems quite simple thing to me, but I get an error. # map declaration               map = new OpenLayers.Map('map_element', {...     projection : _projObj.mercator,     displayProjection : _projObj.wgs84     }      # layer declaration               uses_layer = new OpenLayers.Layer.Vector("US Layer", {         strategies : [new OpenLayers.Strategy.Fixed()],         projection : _projObj.wgs84,         visibility : true,         protocol : new OpenLayers.Protocol.WFS({             version : '1.0.0',             url : '[...]wfs.map&service=WFS',             featureType : USType         })     });      # center to layer extent After the layer gets successfully loaded, I want to zoom to the vector layer extent:               map.zoomToExtent(uses_layer.getDataExtent());      and I get:               Uncaught TypeError: Cannot call method 'getCenterLonLat' of null     OpenLayers.Map.OpenLayers.Class.zoomToExtentOpenLayers.js:512     addSelectedLayermap.js:87     (anonymous function)46.105.19.68:36     f.event.dispatchjquery-1.7.2.min.js:3     f.event.add.h.handle.i      When I pass the same commands from JS debug command line (embedded in Chrome), the map get zoomed correctly; nonetheless,               uses_layer.getDataExtent();      from command line, gives me               bottom: 5134852.1720288     centerLonLat: null     left: 1744085.9879363     right: 1787837.9987162     top: 5150209.4139235      Now, `centerLonLat` is `null`, but in the docs it is said that the only reqested value is a bound with coordinates array, that IMHO fits with the `uses_layer.getDataExtent();` request. Any hint? Thanks :)


title: Alternative to IMapServer2.ExportScaleBar

text: Does anyone know a good alternative to the IMapServer2.ExportScaleBar function? I want to do essentialy the same task, but outside the ArcGIS Server environment. For a given map scale, predefined scale bar style, scale bar size and image resolution, I'd like to generate a scale bar image. I know I could add a scale bar to the page layout and then export the whole layout, but I do not want to clutter and/or modify the page layout if at all possible. This is the closest I was able to get, but works only if the resolution used is 96 DPI (screen resolution).               IMapSurround mapSurround = null;     IMap map = GetMap(); // returns the map          var originalScale = map.MapScale;          try     {         var scaleBar = DefaultSymbols.ScaleBar; // returns a scalebar with some predefined properties         var scaleBarUid = new UIDClass() {Value = scaleBar.GetType().GUID.ToString("B")};              mapSurround = map.CreateMapSurround(scaleBarUid, scaleBar);              var resolution = 300;              map.MapScale = 1000;              mapSurround.Refresh();         ((IScaleBar)mapSurround).UseMapSettings();              var querySize = (IQuerySize)mapSurround;         double widthPoints = 0;         double heightPoints = 0;              querySize.QuerySize(ref widthPoints, ref heightPoints);              var widthInches = 3.0; // scale bar is 3 inches wide         var heightInches = widthInches * (heightPoints / widthPoints); // determine height by respecting the aspect ratio              var pixelsWidth = (int)(widthInches * resolution);         var pixelsHeight = (int)(heightInches * resolution);              using (var bitmap = new Bitmap(pixelsWidth, pixelsHeight, PixelFormat.Format32bppRgb))         {             bitmap.SetResolution(resolution, resolution);                  using (var graphics = Graphics.FromImage(bitmap))             {                 graphics.Clear(Color.White);                      var simpleDisplay = new SimpleDisplayClass();                 simpleDisplay.DisplayTransformation.Resolution = resolution;                 simpleDisplay.DisplayTransformation.Units = esriUnits.esriInches;                 simpleDisplay.DisplayTransformation.Bounds = new EnvelopeClass                      {                          XMin = 0,                          YMin = 0,                          XMax = widthInches,                          YMax = heightInches                      };                 simpleDisplay.DisplayTransformation.VisibleBounds = simpleDisplay.DisplayTransformation.Bounds;                      var tagRect = new tagRECT {left = 0, top = 0, right = pixelsWidth, bottom = pixelsHeight};                 simpleDisplay.DisplayTransformation.set_DeviceFrame(ref tagRect);                      simpleDisplay.StartDrawing(graphics.GetHdc().ToInt32(), 0);                      var drawBounds = new EnvelopeClass                      {                          XMin = 0,                          YMin = 0,                          XMax = widthInches,                          YMax = heightInches                      };                      mapSurround.Draw(simpleDisplay, null, drawBounds);                      simpleDisplay.FinishDrawing();             }             bitmap.Save("c:\\image.bmp", ImageFormat.Bmp);         }     }     finally     {         map.MapScale = originalScale;              if (mapSurround != null)         {             map.DeleteMapSurround(mapSurround);         }     }      For resolutions other than 96 DPI this does not work correctly, the scale bar IS drawn to the bitmap, but the measurements it shows are not correct. Is there a way to do this? Any help and ideas are greatly appreciated.


title: How to Geocode a City

text: First, a little disclaimer. I know that I can just call one of any number of geocoding APIs to perform this for me. This is more of a curiosity than anything else. With that aside, how does someone (like Google) turn a city name into a latitude and longitude? From playing around with Google Maps it appears that most of the time a city is pinpointed to it's downtown area, but there seems to be no reason why one street corner is picked over any other street corner. I am just curious what algorithm a geocoding service provider would utilize for performing such an action. I hope that makes sense, and thanks for satisfying my curiosity.


title: Get the vertices on a LineString either side of a point (Sql Server 2008)

text: I have a point that sits on top of a LineString, both Sql Geometry types.               POINT('0.5 0.5')          LINESTRING('0 0 NULL 1, 1 1 NULL 1')      I want to find the vertices of the LineString that lie either side of the point so that I can interpolate the Z and M values. Based on the values above my desired result would be:               POINT('0.5 0.5 NULL 0.5')      Can this be achieved in Sql Server 2008 or is there an add-on available to support this? In reality I would be querying a LineString containing many points, not just two as shown in the example above.


title: How to update constraints if the table name was changed?

text: I modified the name of a Table in postGIS. But the constraints contains old table name. How can I change the constraints? Please help me.


title: How Can I Install "QGIS Server" in a Windows 7 PC?

text: I am trying to install the new "QGIS Server" feature that comes with QGIS 1.6 on my Windows PC. The directions for this at the QGIS website only show instructions for Mac OS, and those in other GIS blogs show instructions for installing the feature on Linux machines. I'm wondering if instructions for installation on Windows are available anywhere, or whether anybody can provide me with the steps for installing the feature on the Apache web server on my PC. Thank you in advance.


title: VBA field calculator: concatenating fields

text: I'm trying to get one column of text to move into another column that already has text. I DON'T want to replace it, I want the transferring text to be added to the already existing text in a shp file attribute table. ### Example               permit   acct id     123      567     324      456      I want to take the "acct id" numbers and bring them into the permit column and have it look like this...               permit     123567     324456      thank you, Z


title: What can I do in qgis to manage a WMS frequently unreachable?

text: In my project I have a layer linked to an external WMS server (provided by a local institution). The problem is that it is somewhat unreliable: it is down rather often. Is there something I can do to correctly manage this situation? At the moment it simply removes the layer, so if I were to save the project I should re-add the WMS later (when it is back on, unpredictable). A far better behaviour would be just ignoring it and show some kind of warning when trying to get data out of it


title: QGIS 1.8 Desktop and Ubuntu 11.10

text: I'm working with QuantumGis 1.8 Desktop on Ubuntu 11.10. Two days ago the program crashed while I was working with some WMS layers and then it hasn't been working anymore. When I open QuantumGIS Desktop the program starts, the logo appears but nothing more happens and then Ubuntu suddenly quit the application. I have tried several times both to remove QuantumGIS and to reinstall it, but the problem remains. Could you help me? Greetings, Jacopo


title: Understanding arcpy log files in the scratchWorkspace

text: When I run a tool through arcpy, a file simply called 'log' with no extension is created in the scratchWorkspace. This is quite different to the history log files which are XML and I cannot find any documentation on it. A line of this file typically looks something like:               201402261409   0     0     0Johnsmith C:\temp\tmpeoqn0r\SnapPou_fA551 = SNAPPOUR(C:\temp\tmpeoqn0r\t_t324, C:\Test Data\SampleFlowAcc\FlowAcc.tif,  10000.0)      So it seems the first column is a data/time stamp of the tool execution, then there are 2 numbers, followed by what must be the windows user who executed the tool and the tool command (with output & input locations) I would like to know what the number columns mean? More often than not, these are both 0's, but sometimes the first column is a 1 and the second column is anything between 1 and 9... Any ideas? Here is a bigger example of the file contents:               201402261410   1     7     0Jhonsmith C:\temp\tmpeoqn0r\Watersh_Flow1 = WATERSHED(C:\Test Data\SampleFlowDir\FlowDir_Fill11.tif, C:\temp\tmpeoqn0r\SnapPou_fA641)     201402261410   0     0     0Johnsmith in_memory\f55D3000E_BE6E_4FF6_A15E_3B3ACA024F7E = GRIDSHAPE(C:\temp\tmpeoqn0r\ws_6, WEED)     201402261410   0     0     0Johnsmith C:\temp\tmpeoqn0r\SnapPou_f90E1 = SNAPPOUR(C:\temp\tmpeoqn0r\t_t3211, C:\Test Data\SampleFlowAcc\FlowAcc.tif,  10000.0)     201402261410   0     9     0Johnsmith C:\temp\tmpeoqn0r\Watersh_Flow1 = WATERSHED(C:\Test Data\SampleFlowDir\FlowDir_Fill11.tif, C:\temp\tmpeoqn0r\SnapPou_f90E1)     201402261410   0     0     0Johnsmith in_memory\fAC52111C_9237_4398_897D_1E5D0A460157 = GRIDSHAPE(C:\temp\tmpeoqn0r\ws_7, WEED)     201402261410   0     0     0Johnsmith C:\temp\tmpeoqn0r\SnapPou_f8541 = SNAPPOUR(C:\temp\tmpeoqn0r\t_t3212, C:\Test Data\SampleFlowAcc\FlowAcc.tif,  10000.0)


title: gp.Describe() on sde table

text: Python newbie needing help with a simple task. I'm trying to create a list of stand-alone tables in an sde workspace. Code below works just fine to list feature class names and paths and table names.               fcList = gp.ListFeatureClasses()     for fc in fcList:         descFC = gp.Describe(fc)         print "FC Name: " + fc + ", Path: " + descFC.CatalogPath          tblList = gp.ListTables()     for table in tblList:         print table      However, I am not able to use the Describe method to grab the tables' catalog path.               tblList = gp.ListTables()     for table in tblList:         print table         descTBL = gp.Describe(table)         print descTBL.CatalogPath      I think my issue is in the gp.Describe line. What's wrong with gp.Describe(table)? From looking at code samples fc = feature class and fd = feature dataset. What abbreviation is used for table?


title: Developing Web Application For Map

text: I have huge data of maps usually most of the table are co-ordinates and the other are raster maps stored in SQL server 2008 and PostGres Server. I am using desktop application to analyse the data. Usually the desktop end is QGIS or mapinfo. Now I think It would be easier if end user could do the same with the same dataset. So I am building a web application so that those data could be queried and manipulated by end user as required and view as layer. How to get started with this? I had a look at this and this. > The requirement :   1. User should be able to query data from server and display it as map.   2. Since there are multiple tables they should come as multiple layers and user can disable and enable the layer.   3. There is no relation what so ever with google/bing/or third party maps ( may be this could be quite a feature for version 2) Actually I want to build simple web application that acts as basic GIS application. Well tell me if I am reinventing the wheel. I found this link helpful but could not provide step by step solution to the problem


title: How can I calculate driving distance between 2 lat/long points

text: I want to build a .NET application that can get an input of 2 lat/long points and return an output that specified the **driving** distance between the points. Appreciate suggestions for API's or methods to perform that task. Thanks


title: Convert PS/EPS/AI to shapefile?

text: Are there any programs which can convert from Postscript, Encapsulated Postscript, or Adobe Illustrator file formats to a shapefile? What about adding geo-referencing? I'm also open to recipes using a number programs if there isn't a single one to fit the bill.


How do I figure out the longitude and latitude coordinates necessary to create a circular 1 mile radius around one location? I don't mind how many coordinates that takes. For instance: Latitude = 28.4789 Longitude = -81.4682 What mathematical **theorem** or **formulas** would you use to accomplish a task like this? (I'm **not** using google maps, I just need to generate the necessary longitude and latitude to form the radius)


title: Difference between Geoserver stores: ImageMosaic, WorldImage, ImagePyramid

text: What are the differences between the Geoserver stores ImageMosaic, WorldImage, and ImagePyramid? I noticed that ImagePyramid is no longer included in the latest Geoserver. Is one of these more efficient than the others?


title: Permanently reorder fields in data (shp, gdb)

text: If I have a shapefile or geodatabase featureclass that has fields that I want to permanently reorder, is there a tool, script, or utility that will help me quickly accomplish this task? Generally I know you have to create a new field, then populate it with the old values, then delete the old field. I am hoping there is something already created out there that will help with this. Assuming you don't have FME or another ETL tool available, what is the next best suggestion? This is NOT a view level change (ie- drag/drop the fields in order), but rather write out to a new shapefile or featureclass. Thanks.


title: digital elevation model data help?

text: I know little about geographic data. I am trying to get a file of elevation data for the United Kingdom so that I can import it into a CAD program and convert it to an *.stl or *.dxf format. I intend to use a 3D CNC router to make a physical model of this data. Can anyone point me anywhere or explain how I would acquire data in this format of data I can convert?


title: How to add a TIF raster layer programmatically?

text: This is my code whenever i try to load .tif image it shows layer is not valid,`QgsRasterLayer` always return invalid layer (I am using linux 12.04).               QString myFileName = QFileDialog::getOpenFileName(this,tr("OpenFile"),"/home/ravi",tr("*.tiff"));     std::cout<<"File Name:"<<myFileName.toStdString();     QFileInfo myRasterFileInfo(myFileName);     QgsRasterLayer * mypLayer = new QgsRasterLayer(myRasterFileInfo.filePath(),     myRasterFileInfo.completeBaseName());     if (mypLayer->isValid())     {        qDebug("Layer is valid");     }     else     {        qDebug("Layer is NOT valid");        return;     }


title: How to group a multiple rasters in MapServer?

text: I had a map files for each state (raster file) of India. Now I want group all mapfiles to show a single raster or single map file? How I can do in using qgis or any platform?.