Load information about a feature

I must mention first of all that I am very new to this area... I am developing a GIS application. I created a simple map using geoserver with openlayers, it is all working through apache tomcat - and is connected to postgis/postgresql database. When I click on a point (the points represent weather stations - and stored in database) it is set to display the tablewith all the data like: fid, station_name etc. What I want to do is - to store station_name(I clicked on) in a variable and use for my needs. In this case I was trying to break response.responseText in parts, and get it from there but it never worked... response.responseText - returns all the information and creates the tables. I need to only get ONE value from database. Any ideas!? Any help greatly appriciated!!! Here is the code: In this part it is set to read a click and make a response:               map.events.register('click', map, function (e) {                          document.getElementById('nodelist').innerHTML = "Loading... please wait...";                     var params = {                         REQUEST: "GetFeatureInfo",                         EXCEPTIONS: "application/vnd.ogc.se_xml",                         BBOX: map.getExtent().toBBOX(),                         SERVICE: "WMS",                         INFO_FORMAT: 'text/html',                         QUERY_LAYERS: map.layers[0].params.LAYERS,                         FEATURE_COUNT: 50,                         Layers: 'map_with_stations',                         WIDTH: map.size.w,                         HEIGHT: map.size.h,                         format: format,                         styles: map.layers[0].params.STYLES,                         srs: map.layers[0].params.SRS};                          // handle the wms 1.3 vs wms 1.1 madness                     if(map.layers[0].params.VERSION == "1.3.0") {                         params.version = "1.3.0";                         params.j = parseInt(e.xy.x);                         params.i = parseInt(e.xy.y);                     } else {                         params.version = "1.1.1";                         params.x = parseInt(e.xy.x);                         params.y = parseInt(e.xy.y);                     }                          // merge filters                     if(map.layers[0].params.CQL_FILTER != null) {                         params.cql_filter = map.layers[0].params.CQL_FILTER;                     }                      if(map.layers[0].params.FILTER != null) {                         params.filter = map.layers[0].params.FILTER;                     }                     if(map.layers[0].params.FEATUREID) {                         params.featureid = map.layers[0].params.FEATUREID;                     }                     OpenLayers.loadURL("http://localhost:82/geoserver/wms", params, this, setHTML, setHTML);                     OpenLayers.Event.stop(e);                 });             }                  // sets the HTML provided into the nodelist element             function setHTML(response){                      document.getElementById('nodelist').innerHTML = response.responseText;                            };


Why do my map labels disappear when publishing as a service on ArcGIS server?

I am developing a map service, rather basic, with some polygons and labels. Surprisingly, when I deploy this service on the production server, the labels disappear (see here). I have a second local server - for this one, no problem: the labels are well and truly there. Do you have an idea why my labels vanish that way? Is there some special tuning I should apply on the production server to have them back? Many thanks for your help. **Precision:** It happens with labels with field information. If I use a static text as a labels, the labels are there. The server seems unable to include object fields info within labels. **Additional info:** Apparently, this is a well-known problem in ArcGIS server 10.1. See this thread on ArcGIS forum. No solution is known yet...


How to obtain map service url from within Server Object Extension (.NET)

I'm having trouble figuring out how to obtain the full url to the map service a Server Object Extension is running on top of. I bascially want:               IMapServer4 mapServer = (IMapServer4)serverObjectHelper.ServerObject;     string fullUrl = mapServer.url       Which would yield: http://sampleserver6.arcgisonline.com/arcgis/rest/services/Census/MapServer Any Ideas? Thanks for any suggestions.


Help with Python Update Cursor syntax

I'm having trouble running a loop to update the values in a particular field in attribute table by reading from a csv file. Any input is greatly appreciated. The following is the sample code:               import arcpy          featureclass = "C:\\xyz"     table1 = "C:\\pqr.csv"     alpha= list(csv.reader(open(table1,'r')))          crows=arcpy.UpdateCursor(featureclass) # guess here it should be like cursor( , )          for crow in crows:         for row in alpha:             if crow.STNAME==row[3] and crow.ZIP==row[6]:                 if int(crow.STNUM)>=int(row[13]) and int(crow.STNUM)<=int(row[14]):                     if row[4]==crow.STTYPE and row[15]=="B":                        # fill a field in attribute table with a particular row[] value from csv                     else:                         if int(crow.STNUM)%2==0 and row[4]==frow.STTYPE and row[15]=="E":                       # fill the field in attribute table with a particular row[]value from csv


How to create a 4096x4096 png from a 1440x720 data set using GDAL?

I'm having a bear of a time trying to do something that seems fairly straightforward. I'm trying to take a netcdf file from ftp://eclipse.ncdc.noaa.gov/pub/OI-daily-v2/NetCDF/2011/AVHRR-AMSR/ and turn it into a 4096x4096 png so that I can tile it up for use w/google maps. I've been using the following commands but it always seems to create blank areas on the east and west borders:               /usr/local/bin/gdal_translate -a_srs EPSG:4326 -outsize 4096 4096 -a_ullr 0 90 360 -90 NETCDF:"input.nc":sst -of netCDF output.nc     /usr/local/bin/gdaldem color-relief output.nc -b 1 -alpha colorscale.txt -of PNG output.png     /usr/bin/convert output.png -transparent #FFFFFF output.png      Am I just approaching this the wrong way? Am I using the wrong tools? FWIW - gdal_translate is v1.8.0 Thanks in advance. **UPDATE:** "Answer" below from Dija - colleague of mine - is actually a response to the comments on the question itself. **UPDATE** Accidentally added an answer instead of a comment. Moving it here. I thought I'd weigh in here. I've been working on this issue for a little while now. I've tried many iterations of the process. I scaled these images to the 4096x4096 as it actually creates a better image than if not. Meaning, if I leave the image in it's native size, 1440x720 I get a larger gap at the prime meridian (for this data set) and at the 180th meridian for other data sets. I scaled the output image to be a multiple of 256 (to fit better with google) and get various results. If I scale it to 1536x768 the gap shows up to the right of the 0 meridian. If I scale to 1280x640 it's to the left. If I leave it alone, it's pretty much splits it. I have taken the base image that gdaldem spits out, cut it in half and matched the edges up...they line up perfectly. This leads me to think that there's something weird going on with gdal2tiles.py I wrote my own tiler that I'm getting ready to turn loose on our images to see if the gap remains. I'll put together some screen shots to show what I am talking about with how the gap as well as report on what my tiler does. Thanks for any input you might have. It's appreciated. **UPDATE 2** I ran the maptiler I wrote against the original raster and it works fine with out a gap. gdal2tiles.py seems to be the problem. It introduces a transparent edge to the last tile that touches the 180th meridian.


How do use the Shapefile Converter in Oracle 12c to load a Shapefile?

**Update:**   I resolved the service name issue. At some point I put in the right combo of parameters (or maybe I omitted some) and this message was displayed:   `The following key/value is missing: -s db_sid OR -sn db_service_name`   The `-sn` option isn't in the docs as far as I can tell. In my case, I replace the `-s GDSRV` with `-sn pdborcl`.               Now I'm getting all sorts of errors:       Connecting to Oracle Database using...       localhost, 1521, null pdborcl, GISD, ********, shapes, tl_2013_us_state, null, 0       Dropping old table...          Table not previously created.       Creating new table...       java.sql.SQLSyntaxErrorException: ORA-00902: invalid datatype            java.sql.SQLSyntaxErrorException: ORA-00942: table or view does not exist          Record #1 not converted.       java.sql.SQLSyntaxErrorException: ORA-00942: table or view does not exist          Record #2 not converted.       java.sql.SQLSyntaxErrorException: ORA-00942: table or view does not exist          ...            Record #56 not converted.     56 record(s) not converted.     0 record(s) converted.     Done.      **Original Question:** I am attempting to use the shapefile converter in Oracle12c and I'm not having much success and I think it is due to the container database (CDB) vs pluggable database (PDB). I need to somehow specify my PDB: `pdborcl` If I provide the SID, I get an invalid username/password....as expected. `C:\temp>java oracle.spatial.util.SampleShapefileToJGeomFeature -h localhost -p 1521 -s GDSRV -u GISD -d gisd -t shapes -f tl_2013_us_state -r 0 -g geom`               host: localhost     port: 1521     sid: GDSRV     db_username: GISD     db_password: ********     db_tablename: shapes     shapefile_name: tl_2013_us_state     SRID: 0     db_geometry_column: geom     Connecting to Oracle Database using...     localhost, 1521, GDSRV null, GISD, ********, shapes, tl_2013_us_state, null, 0     java.sql.SQLException: ORA-01017: invalid username/password; logon denied                  at oracle.jdbc.driver.T4CTTIoer.processError(T4CTTIoer.java:450)             at oracle.jdbc.driver.T4CTTIoer.processError(T4CTTIoer.java:392)             at oracle.jdbc.driver.T4CTTIoer.processError(T4CTTIoer.java:385)             at oracle.jdbc.driver.T4CTTIfun.processError(T4CTTIfun.java:938)             at oracle.jdbc.driver.T4CTTIoauthenticate.processError(T4CTTIoauthenticate.java:480)             at oracle.jdbc.driver.T4CTTIfun.receive(T4CTTIfun.java:655)             at oracle.jdbc.driver.T4CTTIfun.doRPC(T4CTTIfun.java:249)             at oracle.jdbc.driver.T4CTTIoauthenticate.doOAUTH(T4CTTIoauthenticate.java:416)             at oracle.jdbc.driver.T4CTTIoauthenticate.doOAUTH(T4CTTIoauthenticate.java:825)             at oracle.jdbc.driver.T4CConnection.logon(T4CConnection.java:596)             at oracle.jdbc.driver.PhysicalConnection.<init>(PhysicalConnection.java:715)             at oracle.jdbc.driver.T4CConnection.<init>(T4CConnection.java:385)             at oracle.jdbc.driver.T4CDriverExtension.getConnection(T4CDriverExtension.java:30)             at oracle.jdbc.driver.OracleDriver.connect(OracleDriver.java:564)             at java.sql.DriverManager.getConnection(Unknown Source)             at java.sql.DriverManager.getConnection(Unknown Source)             at oracle.spatial.util.SampleShapefileToJGeomFeature.main(SampleShapefileToJGeomFeature.java:341)     Exception in thread "main" java.lang.NullPointerException             at oracle.spatial.util.SampleShapefileToJGeomFeature.main(SampleShapefileToJGeomFeature.java:345)      From SQLPLUS, I would use the PDB name instead to connect and I'm in. How do you connect to the PDB when calling oracle.spatial.util.SampleShapefileToJGeomFeature? `java oracle.spatial.util.SampleShapefileToJGeomFeature -h localhost -p 1521 -s pdborcl -u GISD -d gisd -t shapes -f tl_2013_us_state -r 0 -g geom`               host: localhost     port: 1521     sid: pdborcl     db_username: GISD     db_password: ********     db_tablename: shapes     shapefile_name: tl_2013_us_state     SRID: 0     db_geometry_column: geom     Connecting to Oracle Database using...     localhost, 1521, pdborcl null, GISD, ********, shapes, tl_2013_us_state, null, 0     java.sql.SQLException: Listener refused the connection with the following error:     ORA-12505, TNS:listener does not currently know of SID given in connect descriptor                  at oracle.jdbc.driver.T4CConnection.logon(T4CConnection.java:673)             at oracle.jdbc.driver.PhysicalConnection.<init>(PhysicalConnection.java:715)             at oracle.jdbc.driver.T4CConnection.<init>(T4CConnection.java:385)             at oracle.jdbc.driver.T4CDriverExtension.getConnection(T4CDriverExtension.java:30)             at oracle.jdbc.driver.OracleDriver.connect(OracleDriver.java:564)             at java.sql.DriverManager.getConnection(Unknown Source)             at java.sql.DriverManager.getConnection(Unknown Source)             at oracle.spatial.util.SampleShapefileToJGeomFeature.main(SampleShapefileToJGeomFeature.java:341)     Caused by: oracle.net.ns.NetException: Listener refused the connection with the following error:     ORA-12505, TNS:listener does not currently know of SID given in connect descriptor                  at oracle.net.ns.NSProtocolStream.negotiateConnection(NSProtocolStream.java:272)             at oracle.net.ns.NSProtocol.connect(NSProtocol.java:263)             at oracle.jdbc.driver.T4CConnection.connect(T4CConnection.java:1360)             at oracle.jdbc.driver.T4CConnection.logon(T4CConnection.java:486)      Here's the docs on the 12c converter: http://docs.oracle.com/cd/E16655_01/appdev.121/e17896/sdo_shapefile_converter.htm#SPATL1427


pgShapeLoader: how to load a shapefile in a different schema than 'public'

I'm trying to load a shapefile in POSTGIS using pgShapeLoader. It works but I can load table only in the "public" schema and i don't know how to configure, the tool to upload my shapefile in another schema. Does it possibile? Any suggestion for workaround? Anyway I can use shp2pgsql and then psql from command line in POSTGIS .... Thank you very much in advance! Cesare


model only runs open "modelbuider"

I am using ArcGIS 10.2.2 for Desktop. I have a model made in ModelBuilder which currently does:   1. input XY coordinates excel table    2. data display XY event    3. converting a format event shape    4. add field "name1"    5. add field "name2"  So far so good but my problems are:   1. the model only works so edition    2. I cannot add tool to make a selection of my new shapefile ![enter image description here](http://i.stack.imgur.com/42wh5.png)


Change in behaviour of layer extents

I have a class that inherits a vector memory layer, it all works fine a previous version of QGIS, but on 2.2 it no longer works. It either does not know that the layer is WGS84 or getting the extents is somehow failing.               class DotTrail(QgsVectorLayer):         def __init__(self,name, colour):             QgsVectorLayer.__init__(self, "multipoint?crs=WGS84", name , "memory")             self.pr =self.dataProvider()         self.featId = -1          self.points = []                       # set the colour             rendererV2 = self.rendererV2()             layerSymbol = QgsSymbolV2.defaultSymbol(self.geometryType())             layerSymbol.setColor(QtGui.QColor(colour))             layerSymbol.setAlpha(1)             rendererV2.setSymbol(layerSymbol)             QgsMapLayerRegistry.instance().addMapLayers([self])              def addPointUpdate(self, p):             self.points.append(p)             self.update()              def addPoint(self, p):             self.points.append(p)              def clear(self):             self.points = []             #self.update()              def update(self):             print 'Update'         if self.featId == -1:             feat = QgsFeature()             feat.setGeometry(QgsGeometry.fromMultiPoint(self.points))             res, fid = self.pr.addFeatures( [feat] )             self.featId = fid[0].id()             for feat in self.pr.getFeatures():                 print feat.geometry().asMultiPoint()                 else:                 self.pr.changeGeometryValues({self.featId: QgsGeometry.fromMultiPoint(self.points)})                 feat = QgsFeature()                  self.updateExtents()             self.updateFields()         print 'Dot layer extents'         print self.extent().toRectF().getCoords()      It adds the following points Update [(151.193,-33.8895), (151.193,-33.8895), (151.195,-33.8877), (151.195,-33.8913), (151.191,-33.8913), (151.191,-33.8877), (151.195,-33.8877), (151.195,-33.8877)] But outputs the following extent Dot layer extents (0.0, -33.89132551272559, 151.1951840971511, 0.0) I can't see anything wrong with it.


Use of postgis with hibernate

We have an existing postgis database. Our application is developed in Java, but queries database through a jni dll, which appears to be not very effective for our purpose. I've managed to migrate to jdbc, and it's working well, but the next step would be to rely on a framework such as hibernate-spatial, which i don't know at all. How mature is it? I've seen some people here using it, what do you think of it? It looks that it supports only cartesian coordinates and it may be a concern for us for distance computation.


Modify raster layer display and symbology properties programmatically with arcpy or ArcObjects

I have over 1500 NITF files that I need to be able to modify the layer display and symbology of programmatically from within ArcMap. I need to set the properties outlined in red below: ![enter image description here](http://i.stack.imgur.com/4WZFa.png) ![enter image description here](http://i.stack.imgur.com/AMmhj.png) I've looked through the arcpy docs and it doesn't look like these properties are exposed - honestly, I didn't expect them to be. I have a feeling I can get to them through ArcObjects. Anyone had any experience/luck modifying these properties programmatically?


A global, grid-like projection for creating heatmaps

I'm putting together an application where I need to create a vector grid which will be used to store and display a heatmap. It has the following requirements:   * Can cover the entire planet.    * The vast majority of grid squares won't have values.   * I don't want to have to store the grid itself; I'd like to calculate it on the fly.   * The scale of the data used with the grid could vary massively.   * I anticipate wanting grid squares anything from 1km across to 100km across. (I'm aware how many this will be (~510million for 1km, ~51,000 for 100km)).   * Values will be accumulated/aggregated for each grid square.   * Ideally I'd be able to easily use smaller grid cells to calculate values for larger ones rather than store the larger grid cell values.   * I'll be using OpenLayers to drape it over OpenStreetMap.   * I'll be storing it in SpatiaLite or SQLite, so preferably supported by those natively (i.e. for SpatiaLite = a supported CRS; or for SQLite = a pure numbers based system). So my question is: What projection should I use for this grid? Also - is there a good way to designing this? Does anyone know of a good potential solution to this problem or have solved similar before? Or can point me in a useful direction. **Edit** Use-case - basically I'm aggregating up bounding boxes of various different shapes and sizes. They can be anything from a few hectares to thousands of square km in size. They may also be in different projections. Below is a bespoke version of the sort of thing I'm aiming to accomplish automatically on a larger scale. ![enter image description here](http://i.stack.imgur.com/YEvt9.jpg) Many thanks.


QGIS projects - collaboration on different platform

on my computer, I have both - Ubuntu and Windows - operating systems. Since I work on both of them and QGIS save path to particular layers using absolute path, it's very inconvenient to store projects on local disk. Moreover, I wish to access QGIS projects and layers from other computers. Is there any possibility to store them on the Internet? What will be the best, free of charge, solution? I thought about Dropbox, but I'm not sure if there's more convenient way.


How to change the name of "Base layer" and "Overlays" in LayerSwitcher?

I would like to change the name of "Base layer" and "Overlays" in my LayerSwitcher div to something else. I've noticed that the div including the texts only has a class and no ID which means it's little bit tricky to access. I managed to change "Base layer" name with               document.getElementById("layerswitcher").firstChild.firstChild.innerHTML = "new name";      but I don't like this solution. Best would be if I could define the names directly with OpenLayers functions. Is there one? I read something about using language translations but it's not what I'm after.


How can I use a neural network in GIS?

I'd like to learn how to use a neural network in GIS. I surfed in the internet to find good resources or a tutorial about this but I couldn't find any helpful resources. Are there any good references or tutorials about using neural networks in GIS? I'd like to use an artificial neural network (ANN) in mapping hazard zones, for example flood zones or landslides. I use ArcGIS 10 to process the layers but don't know how can use an ANN.


How to calculate the track length of GPX files with PostGIS?

I am trying to calculate the length of a GPX track. The original GPS data looks like this:               <?xml version='1.0' encoding='UTF-8' standalone='yes' ?>     <gpx creator="nl.sogeti.android.gpstracker"        version="1.1" xmlns="http://www.topografix.com/GPX/1/1"        xmlns:gpx10="http://www.topografix.com/GPX/1/0"        xmlns:ogt10="http://gpstracker.android.sogeti.nl/GPX/1/0"        xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"        xsi:schemaLocation="http://www.topografix.com/GPX/1/1        http://www.topografix.com/gpx/1/1/gpx.xsd">         <metadata>             <time>2014-04-25T19:20:39Z</time>         </metadata>         <trk>             <name>Track 2014-04-25 19:20</name>             <trkseg/>             <trkseg>                 <trkpt lat="52.51025036" lon="13.3907098">                     <ele>93.30000305175781</ele>                     <time>2014-04-25T19:59:24Z</time>                     <extensions>                         <gpx10:speed>2.75</gpx10:speed>                         <ogt10:accuracy>7.0</ogt10:accuracy>                         <gpx10:course>236.5</gpx10:course>                     </extensions>                 </trkpt>                 <trkpt lat="52.51021522" lon="13.390611">                     <ele>94.0</ele>                     <time>2014-04-25T19:59:26Z</time>                     <extensions>                         <gpx10:speed>2.5</gpx10:speed>                         <ogt10:accuracy>7.0</ogt10:accuracy>                         <gpx10:course>235.1999969482422</gpx10:course>                     </extensions>                 </trkpt>                 <!-- And many more -->             </trkseg>         </trk>     </gpx>      I imported the data into PostgreSQL/PostGIS. The track points have been imported into the table `track_points`:               gpxdata=# \d track_points               Table "public.track_points"            Column       |           Type           |     --------------------+--------------------------+      ogc_fid            | integer                  |       wkb_geometry       | geometry(Point,4326)     |       track_fid          | integer                  |       track_seg_id       | integer                  |       track_seg_point_id | integer                  |       ele                | double precision         |       time               | timestamp with time zone |       magvar             | double precision         |       geoidheight        | double precision         |       name               | character varying        |       cmt                | character varying        |       desc               | character varying        |       src                | character varying        |       link1_href         | character varying        |       link1_text         | character varying        |       link1_type         | character varying        |       link2_href         | character varying        |       link2_text         | character varying        |       link2_type         | character varying        |       sym                | character varying        |       type               | character varying        |       fix                | character varying        |       sat                | integer                  |       hdop               | double precision         |       vdop               | double precision         |       pdop               | double precision         |       ageofdgpsdata      | double precision         |       dgpsid             | integer                  |       gpx10_speed        | double precision         |       ogt10_accuracy     | double precision         |       gpx10_course       | double precision         |       Now I use the following query which I found here.               SELECT gps.ogc_fid, tracks.name,      ST_Length_Spheroid(ST_MakeLine(gps.wkb_geometry),      'SPHEROID["WGS 84",6378137,298.257223563]') AS track_len         FROM (             SELECT ogc_fid, time, wkb_geometry              FROM track_points              ORDER BY ogc_fid, time         )          AS gps, tracks         WHERE gps.ogc_fid = tracks.ogc_fid         GROUP BY gps.ogc_fid, tracks.name         ORDER BY gps.ogc_fid;      The result of the query is:                  | ogc_fid |          name          |    track_len     |        | integer |    character varying   | double precision |     ----------------------------------------------------------      1 |    1    | Track 2014-04-25 19:20 |        0         |      * * * The following query posted by Brad ..               SELECT ogc_fid, time, ST_AsEWKT(wkb_geometry) FROM track_points LIMIT 5;      ... produces this result:               ogc_fid |           time           | st_asewkt     integer | timestamp with time zone | text     ------------------------------------------------------------------------------        1    |  2014-04-25 19:31:13+02  | SRID=4326;POINT(13.38866662 52.54784071)         2    |  2014-04-25 19:31:14+02  | SRID=4326;POINT(13.3886988 52.54778812)        3    |  2014-04-25 19:31:15+02  | SRID=4326;POINT(13.38873776 52.54773285)        4    |  2014-04-25 19:31:16+02  | SRID=4326;POINT(13.38877592 52.5476608)        5    |  2014-04-25 19:31:17+02  | SRID=4326;POINT(13.38881506 52.54760139)


Regression analysis - creating dataset

![enter image description here](http://i.stack.imgur.com/895yf.png) I have a vector-based dataset of soils/parent material. The statistics have already been done, with coefficients for each category of parent material. I am not sure how to use this information to create a soils raster layer for analysis within Raster Calculator. For example, to create the model result raster, I'm thinking: = 0.158544748 + (-0.422151664 * LAT) + (01.184392010 * LONG) + (0.001245377 * ELEV) ...and so on for all the variables. BUT, how do I create the raster for the soils/parent material information, and how would that look in the above equation for the raster calculator? (they are labeled PARMAT AL, PARMAT GR, PARTMAT GRHW... in the table) Thanks! EDIT- I'm using ArcGIS 10.1, have Spatial Analyst and many other extensions (not sure what's relevant). Also, I was previously advised to 'create a "ParamEst" column for each of the PARMAT (soils) values, to be populated with the parameter estimates, then convert to raster. Each PARMAT type would then be a parameter estimate and can be added as a standalone layer.' Therefore, would each PARMAT type be exported out as its own raster, with its parameter estimate, and simply added to the model? So instead of + (-1.121263096 * PARMATAL), would it simply be + (PARMATAL) in the equation, because the derived raster is already a parameter estimate? I hope I'm making sense - I only have a basic grasp of statistics, so I'm not sure if I'm asking things properly. ***I'm unable to reach that person again, which is why I haven't had further clarification.


postgis -- return records sorted by size of st_intersection

I have two sets of spatial records in two different tables. I'm in rails so I'll show my code in that way seeing that its source: one set is CensusPlace records, which are US Census Tigerline shapefiles. Note that Census Places are typically towns and cities (with some sub- municipalities) that can exist in multiple counties. the other set of records are County records. It would be very handy if I could run a postgis command that would return my st_intersects results for a given city like so:               CensusPlace.where(:state => "SC").each do |cp|         puts "\n---------------census place: #{cp.name}"         County.where{st_intersects(:county_shape, cp.census_place_shape)}.each_with_index do |c, i|           puts "#{c.name} county"           puts "\toverlap ranking: #{i+1}"              end      end      such that the actual returning of results would be ordered by the the overlap size. It would be natural for the results to actual return a new "virtual column" or what have you with a 0-1 value of that overlap percentage as well. Thanks


How to set symbology based on an attribute of a 1-N table

Now, it is possible to edit geographic features and add elements in child tables using 1-N relationships. I want to construct a stacked - symbology which allows me to represent the multiple attributes. For example if I have a line feature "linear park" and my linear park has a linked table (1-N) called "users", i want a stacked or multi-level symbology (with some offsets) that allows me to depict all the "users" of that section of park. Is it possible in QGIS? I know it was possible in ArcGIS using relates. I appreciate your time!!


How to indicate a range of zoom on the map?

Basically what I need and I want to know is how you can do to tell the map that respects me a zoom range to display the data, ie, I want to see the map from afar without information, and as the zooming you can go see the information. Can anyone help me with this?


How to convert .shp into .gpx using QGIS?

I'm newbie using QGIS. I would like to convert .shp into .gpx. So, using Qgis, I open my shapefile and I tried to "save as" gpx but it gives an OGR error: `creation of field ObjectId failed (OGR error: Field of name 'ObjectId' is not supported in GPX schema. Use GPX_USE_EXTENSIONS creation option to allow use of the <extensions> element.)` How can I fix this error? Thanks.


QGIS - Save as GPX

I would like to convert .shp into .gpx. I tried to "save as" gpx but it gives an OGR error: creation of field ObjectId failed (OGR error: Field of name 'Name' is not supported in GPX schema. Use GPX_USE_EXTENSIONS creation option to allow use of the element.) GPX_USE_EXTENSIONS=Y does not work. Checking 'Skip Attribute Creation' does not work.


How to open QGIS using VBA?

Using VBA, I want to open QGIS (the actual GUI) and execute a simple script through the Python console or the editor.


crust index in Erdas Imagine

hey people, i have the trouble with CI index in Erdas, when i am useing the image from landsat 5 tm with atmospheric correction it gives me the whole white picture in the results, the formula CI=1-((red-blue)/(red+blue)) also i tryed to change R band with NIR, but......does not work at all when i tryed to insert in the model the bands seperately the reslusts are the same, just white pic, can anyone help with this type of issue?


What's a good library for displaying many small maps on a web page

I would like to create a visualization on the web similar to the following: http://namemapper.babynamewizard.com/namemapper/ Are there any good libraries for doing this? Is OpenLayers light-weight enough to handle many maps at once?


TIGER/Line processing

Does anyone know of scripts for merging TIGER/Line shp files? It is for extracting the road features only, but for the whole of USA. They are provided as shp files for one county at a time here: http://www.census.gov/geo/www/tiger/ If anyone sells the data, a link is also appreciated.


Comparing two featureclasses using C#

I'm sorry, I know this is a pretty simple question but am having trouble getting this to work. Basically, I have 2 linear featureclasses and they both have the same fields and spatial reference. FC1 has 675 records and FC2 only has 505 records. So how can I check what Prop_ID(s) are missing from FC2? I have some code that will copy the missing record from FC1 to FC2 but must first determine the missing Prop_ID. many thanks


resize image resolution from 1000 m to 250 m

I'm beginner of ARCGIS and RS. I have many problems with the procedure. One of them is   I have NDVI images and Soil map. Both are at the different resolution: NDVI is 250 m. soil map is 1000 m. So I wanna ask how can I change the resolution from 1000 m to 250 m. Thanks a lot


Qgis 1.8 fails to install

QGIS 1.8 will not start on Windows 7 Basically a similar issue as the above post. Version 1.8 fails to install with the standalone installer; with OSGeo4w, while it doesn't run into trouble during the installation, ultimately fails as well (same error message as below when opening the app). Standalone installer has the following problems: At the end of the installation, it reaches "qgis-common.bat" execution and stalls; meanwhile, Task Manager shows "find string" GREP utility running in the background and repeatedly crashing every minute or so. Ending the process results it immediately starting back up. If the installation is forced to quit or completed, trying to open Qgis results in the error message that **qwt5.dll** is missing. QGREP crash signature:               Problem Event Name:APPCRASH     Application Name:   findstr.exe     Application Version:    6.1.7601.17514     Application Timestamp:  4ce78c0a     Fault Module Name:  findstr.exe     Fault Module Version:   6.1.7601.17514     Fault Module Timestamp: 4ce78c0a     Exception Code: c0000005     Exception Offset:   00003c76     OS Version: 6.1.7601.2.1.0.768.3     Locale ID:  1033     Additional Information 1:   0a9e     Additional Information 2:   0a9e372d3b4ad19135b953a78882e789     Additional Information 3:   0a9e     Additional Information 4:   0a9e372d3b4ad19135b953a78882e789


How to capture data for mapping

Well it might be an newbie question but I can't find it In my project I need the paths in a campus but these paths are inexistents so I'm going to walk around and mapping (I don't know if it's the correct term) the points to build the correspondent paths, so my question is if exists an easy way to do it or how do you think I can accomplished the work, please write all yours thoughts.. Thanks in advance..


QGIS Newbie Rendering Time

I am new to QGIS and have just installed 1.8.0 on a Gateway laptop that is about 2 years old. It is running Windows 7, has an Intel Core i3 2.13 GHz processor, 4 GB of RAM. I've just brought in my first layers (235 MB of data in 3 shapefiles) to my first project, and it takes forever to render. By forever I mean more than one minute for each task: zoom, move, rectangle select, changing the order of layers. The data is PLSS data from the state of Oregon as provided through the federal government through BLM. Although I downloaded only the Oregon data, when I right-click on the twnshp layer onto attribute table, I do see records relating to CA and ID as well as Oregon. Will I not be able to run QGIS on this computer? Thanks so much, Yeye


Create a dbf table by copying the structure of another dbf table

I have a dbf table and I wish to automatically create another dbf table with a different name but with exactly the same structure i.e.fields as the existing table. Is there any ready code in VBA and ArcObjects for that? Thanks


Categorized style won't display in QGIS2.2

I'm a qgis novice trying to put together a building characteristics map. While 44 out of 46 of my layers have their data categorized and displayed correctly, two layers will not display different colours. Instead I am getting the same colour all over my map, which is the colour assigned to an extra unvalued category that appears at the bottom of my list, as shown in the screenshot. ![Styles](http://i.stack.imgur.com/KcpTH.jpg) Any ideas? Thanks in advance


Can you programmatically change the button image for a ESRI.ArcGIS.Desktop.AddIns.Button?

I'm new to the AddIn customization model in ArcGIS 10 and I'm trying to figure out how to change a button image for a button. I'm looking for the equivalent to the BaseCommand UpdateBitmap method and haven't had any luck. Is this a limitation of the AddIn model? Or did I miss something obvious?


How to list raster bands?

I am interested in looping through a directory and listing only the red bands (Band 1) from 4-band NAIP imagery. From there I would like to mosaic all of the band 1 layers into a new single band raster image using Mosaic to New Raster. Is there a way to specify band 1 using `arcpy.ListRasters()` or should I find a new approach?


QGIS Spatial Queries

I was wondering if I could do the following within QGIS: I want to update a column in my point layer with the attribute information from a specific column in a polygon layer - so I have a layer of postcodes and I want to give it a column of say School Catchment and I want to populate this column with the value of the school name column in the catchment area polygon layer. Any help would be greatly appreciated. Thanks! Steve


converting coordinates from unusual DMS format to decimal

So I've been shared a dataset with what to me is an unusual format for degree min second data. In their own words:                   The coordinates provided are longitude and latitude in degrees minutes and seconds.         So for example, coordinates recorded in the spreadsheet as          X_Long: -10.0135 Y_Lat:  53.5726, can be taken as Longitude: 10º 01'35" W,          Latitude: 53º 57'26" N using standard notation.      I need to get this into decimal. Any quick way to do this ? Edit: in trying to solve this I made a following code:                   # Fake data         X_Long <- as.numeric(c(-2.4564, -1.0134, -8.01))              long.temp <- abs(X_Long - trunc(X_Long))         long.min <- trunc(long.temp*100)         long.sec <- trunc(as.integer((long.temp*100 - long.min)*100))      Its giving me weird rounding errors:                   > X_Long         [1] -2.4564 -1.0134 -8.0100         > long.temp         [1] 0.4564 0.0134 0.0100         > long.min         [1] 45  1  0         > long.sec         [1] 63 34 99      This seems to be coming from the trunc function - compare results here                   > long.temp*100         [1] 45.64  1.34  1.00         > trunc(long.temp*100)         [1] 45  1  0      Anyone know why this is happening ?


GPS tracking and Web mapping?

We are developing a tracking software and we are the stage where we test the reliability of our location stigma on the following web-maps   * Google Map   * Nokia HERE    * Open Street Map (OSM) We observe that our location appears accurate on Google maps, however when switched to Nokia or Openstreet there is a discrepancy of around 40-50 metres from our actual location. Most importantly our GPS position accuracy is around 2m. Could it be inaccuracy of map calibration or any similar problems known by anyone?


Compare features in PostGIS with updating attributes

I have two polygonal layers: "existing" and "modified". ![enter image description here](http://i.stack.imgur.com/GRyFn.jpg) Attribute table of "existing" and "modified" layers have same column "feature_id" - number of feature. Please help me create query which write attributes 't' (if features matches) and 'f' (if features not matches) in 'compare_result' column of "modified" layer using ST_Equals function. Now I receive attributes 't' (if features matches) and 'f' (if features not matches) in pgAdmin result window using query:               SELECT  ST_Equals(modified.the_geom, existing.the_geom) AS compare_result FROM  modified, existing      WHERE modified.feature_id=existing.feature_id      but I don't can create query? which write result in attribute table. Thanks! **EDIT:** Solution without creating additional layer:               UPDATE         modified     SET         compare_result = compare_result_temp          FROM         (             SELECT     r."feature_id" AS idmodified,     ST_Equals(r.the_geom, m.the_geom) AS compare_result_temp             FROM          "modified" AS r,          "existing" AS m     WHERE r."feature_id" = m."feature_id"         ) AS foo      WHERE       foo.idmodified = modified."feature_id";


IGeometryBridge.GetPoints throws NotImplementedException. Am I doing something wrong?

I am working on some .NET code that processes lots of polygon geometries at the vertex level. In order to reduce COM interop overhead, I want to call `IGeometryBridge2.GetPoints` once instead of calling `IPointCollection.Point` repeatedly:               // using ESRI.ArcGIS.Geometry;     static IPoint[] GetVertices(this IPolygon singlepartPolygon)     {         var pointCollection = (IPointCollection4)singlepartPolygon;         var points = new IPoint[pointCollection.PointCount];         var geometryBridge = (IGeometryBridge2)new GeometryEnvironment();         geometryBridge.GetPoints(pointCollection, 0, ref points);         return points;     }      `geometryBridge.GetPoints(…)` throws a `NotImplementedException`. (I have tried initialising the `points` array with empty `Point` instances before calling `GetPoints`, but this didn't make any difference. I tried `geometryBridge.QueryPoints(…)`, too, and this failed with a different exception.) Am I doing something wrong, or is this more likely to be a bug?


Bootstrap3 Resize map height for phone

I have a leaflet map in a bootstrap 3 html application. I have content both above and below the the map which has a height of 600px. When I access the application from a phone, the width of the map adjusts to the the width of the screen. However, the height of the map does not update and you end up with a seemingly endless map sandwiched between two chunks of content. The user must then endure endless scrolling. You can see what i mean here: http://punkmap.com/examples/fauxnmap/index.html Anyone know what's the best way to handle this?


Comparing field values to update new field using ArcGIS for Desktop?

I am currently doing an analysis in ArcMap 10 using python. I have two shape layers, one is a point layer and the other one is a line layer. The point layer have attributes for point's unique ID, as well as X and Y coordinates (expressed as coord_X, coord_Y). The line layer include line's ID, and its has coordinates for starting point and ending point. Here they are expressed as FromID_X, FromID_Y, ToID_X, ToID_Y. There are also two empty column in line layer's attribute table as FromID and ToID for points. All coordinates in both point and line layer are already know by field calculation, but the FromID and ToID of starting point and ending point in line layer is empty. I want to figure out the ID of point from point layer and fill them in the table of line layer, i.e. FromID and ToID. I have an idea to compare the coordinates of FromID_X (line's layer) with Coord_X(point's layer), if they are the same, output the Point's ID to line layer as FromID and repeat this for ToID. But now, I am not familiar with python in ArcMap 10, I am wondering if you can tell me some idea or grammar which I can look into so I could use python in ArcMap to do this work? Or you know some tools for such kind of work?


How do I automatically export a release of the EPSG registry?

I am trying to use the EPSG Registry API to export a release of the registry for offline use but am receiving a registry error. After following the Developer Guide I am HTTP POSTing the following XML request to the API URL:               <?xml version="1.0" encoding="UTF-8"?>     <GetRecords          xmlns="http://www.opengis.net/cat/csw"          xmlns:ogc="http://www.opengis.net/ogc"          startPosition="1" maxRecords="100"          outputFormat="application/xml; charset=UTF-8"          resultType="results">       <!--           Retrieve the ReleaseObject releated to the current Version History record.         Then use the ReleaseObject identifier to obtain the actual compressed          GML Dictionary file containing the EPSG Dataset via a getRepositoryItem          request.       -->        <Query typeNames="Association=a RegistryObject=release">         <ElementSetName typeNames="release">full</ElementSetName>         <Constraint version="1.1.0">           <ogc:Filter>             <ogc:And>               <!-- find the correct association -->                <ogc:PropertyIsEqualTo>                 <ogc:PropertyName>$a/targetObject</ogc:PropertyName>                 <ogc:Literal>urn:ogc:def:version-history:EPSG::6.15</ogc:Literal>               </ogc:PropertyIsEqualTo>               <ogc:PropertyIsEqualTo>                 <ogc:PropertyName>$a/associationType</ogc:PropertyName>                 <ogc:Literal>                   urn:x-ogp:def:AssociationType:EPSG:ReleaseFor</ogc:Literal>               </ogc:PropertyIsEqualTo>               <!-- now relate it to the ReleaseObject -->                <ogc:PropertyIsEqualTo>                 <ogc:PropertyName>$a/sourceObject</ogc:PropertyName>                 <ogc:PropertyName>$release/id</ogc:PropertyName>               </ogc:PropertyIsEqualTo>            </ogc:And>           </ogc:Filter>         </Constraint>       </Query>     </GetRecords>      The XML is copied from _section 5.3 'How To... Export a Release'_ of the guide, and the URL to send `GetRecords` requests to is suggested by _section 4.1.1 'GetRecordById'_. The XML error response I get is as follows:               <?xml version="1.0"?>     <ows:ExceptionReport xmlns:ows="http://www.opengis.net/ows" version="2.0">       <ows:Exception exceptionCode="InvalidRequest">         <ows:ExceptionText>XML Schema validation error</ows:ExceptionText>         <ows:ExceptionText>     Error:Error      Line:    2      Column:  193      Message: cvc-elt.1: Cannot find the declaration of element 'GetRecords'.         </ows:ExceptionText>                                                      </ows:Exception>     </ows:ExceptionReport>      I am using the `curl` command line program to run the query as follows:               curl -H 'Content-Type: appliation/xml' --data @request.xml 'http://www.epsg-registry.org/indicio/query'      Can you see where I'm going wrong? If not is there any other way to automatically download the registry? I am aware that I can login to the Registry HTML interface and download the registry from there but I would like to automate the process.


Australian Topo Map GeoTIFFs

Does anyone know if Australian topo map GeoTIFFs are made available anywhere (preferably for free)? LINZ makes them available in New Zealand - http://www.linz.govt.nz/topography/topo-maps/index.aspx ,I was hoping a similar service was available for Australia. Cheers, Gavin


Technical comparison of Oracle and PostgreSQL

I have searched on google and the archives for a technical comparison of Oracle and PostgreSQL (also MySQL), unsuccessfully. Is there any comparison that covers more than just the basic differences? Maybe you can give me some comparison?


geoserver and tiger roads layer overlay problem

Been playing with open layers, geoserver and the tiger roads shapefile but cannot get it to line up correctly over new york. I have looked online and followed the instructions regarding setting the projections and setting google maps to sphericalmercator = true however it still doesnt line up correctly. Any ideas?                       var map;             var untiled;             var tiled;             var pureCoverage = false;             // pink tile avoidance             OpenLayers.IMAGE_RELOAD_ATTEMPTS = 5;             // make OL compute scale according to WMS spec             OpenLayers.DOTS_PER_INCH = 25.4 / 0.28;                  function init(){                 format = 'image/png';                      var bounds = new OpenLayers.Bounds(                     -20037508, -20037508,                     20037508, 20037508.34                 );                      var options = {                     projection: new OpenLayers.Projection("EPSG:900913"),                     displayProjection: new OpenLayers.Projection("EPSG:4326"),                     maxExtent: bounds,                     //maxResolution: 0.00075764453125,                     units: 'm'                 };                     map = new OpenLayers.Map('map', options);                     var GoogleMap = new OpenLayers.Layer.Google(                     "Google Satellite", {                     type: G_SATELLITE_MAP,                     maxZoomLevel: 30,                     "sphericalMercator": true                     }                );                        // setup tiled layer                 wmsLayer = new OpenLayers.Layer.WMS(                     "tiger_roads - Tiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'tiger:tiger_roads',                         styles: '',                         format: format,                         tiled: 'true',                         transparent:'true',                         isBaseLayer: false                     }                 );                      map.addControl(new OpenLayers.Control.MouseDefaults());                  map.addControl(new OpenLayers.Control.Navigation());                 map.addControl(new OpenLayers.Control.ZoomBox({alwaysZoom:true}));                 // Coordinate display at bottom of map                 map.addControl(new OpenLayers.Control.MousePosition());                  map.addLayers([GoogleMap, wmsLayer]);                     var point = new OpenLayers.LonLat(-73.96543, 40.78885);                 // Need to convert zoom point to mercator too                point.transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject());                 map.setCenter(point, 12);                 }


Loading graph (.grf fle) into MXD using ArcPy?

I have hundreds of pairs of FGDB and associated MXD files. I need to insert a graph in each mxd linked to its corresponding FGDB. No problem creating graph, but stuck on how to get graph into mxd file using arcpy. If the graph creation code is run inside the mxd manually, the graph appears, but I need to do it in the python script loop.               import arcpy          workspace = r'C:\gdb_dir'     workspaces = arcpy.ListFiles('*.gdb')        # list of files to process          for workspace in workspaces:         name_prefix = workspace[0:-4]         current_MXD = name_prefix + '.mxd'       # this file will already exist              dataSrc = r'C:\name_prefix\Sample.gdb\DetailFeatures'         fieldX = 'NEAR_X'         fieldY = 'NEAR_Y'         graph_template = r'C:\temp\ScatterGraphTemplate.tee'         out_graph_name = "ScatterTest"         out_graph = r'c:\graphs\' + name_prefix + '.grf'   # How to insert this into mxd?              graph = arcpy.Graph()              graph.addSeriesScatterPlot(dataSrc, fieldY, fieldX)         graph.graphAxis[0].title = "Near Y"                  graph.graphAxis[2].title = "Near X"         graph.graphPropsGeneral.title = "Test of Flake Property Graph"         arcpy.MakeGraph_management(graph_template, graph, out_graph_name)              arcpy.SaveGraph_management(out_graph_name, out_graph, "MAINTAIN_ASPECT_RATIO", 600, 375)      How do you insert out_graph into current_MXD?


Meaning of red cross icon next to Geoprocessing tool/model?

I am using ArcGIS 10.1 for Desktop. When I double-click on the toolbox command line, it is not working and appears as a red cross symbol in command tool icon. So please, what can I do?


Does EO-01 Hyperion L1GST data show DN values, radiance or reflectance?

I am working on Hyperion hyperspectral imagery and I cannot figure out whether L1GST level does do show DN values, radiance or reflectance. * * *


Conversion from Namibian Coordinates to Google Maps (WGS84)

I have coordinates in Namibia, but don't know how to convert them to use them in Google Maps. The system is Lo22/19°. Then Y: +55 435,89 X: -498 157,28 Can anyone help with a conversion?


Apparently MOBAC is not generating MBTiles correctly

I am trying to generate MBTiles with MOBAC, to be used in an offline map Android application (trying to use OSMDroid). Now I generate the atlas, but something goes wrong because I get the error > android.database.sqlite.SQLiteException: no such table: images (code 1): , > while compiling: SELECT tile_data FROM images LIMIT 0,1 This is my MOBAC setup: ![enter image description here](http://i.stack.imgur.com/K0i7i.jpg)


Erdas python scripting help

I want to write scripts for Erdas in Python. I downloaded some Python examples here but I could not make much. Where can I get the full Help material for Erdas python scripting?


Creating a map cache without ArcServer?

Is there a way to create a tiled map cache without using ArcServer? I have aerial photography in .sid format that takes a while to load, whereas cached jpeg tiles will load much much faster. Any way I can work around this issue?


Polygon with the same points over and over

I am in the process of reading shape files (provided to us by the vendor dealing with GIS products) using a C# program and loading them into a SQL database. I first read the text for the shape and then use that text to update my geometry and geography types. And I am doing this as part of a SSIS package. My SSIS package failed on one shape file and threw me an exception: > A .NET Framework error occurred during execution of user-defined routine or > aggregate "geometry": System.FormatException: 24305: The Polygon input is > not valid because the ring does not have enough distinct points. Each ring > of a polygon must contain at least three distinct points. After running some manual updates, I found the record that was giving me grief. Then I took out parts of the multi-polygon and ran individual selects on them, creating individual polygons, eg.               SELECT geometry::STGeomFromText('POLYGON(( 118.501323586697 -20.3203577291617, 118.504216161911 -20.3220101539757, 118.502671059623 -20.3212136027048, 118.501323586697 -20.3203577291617 ))', 4326).MakeValid()      Finally, I found the polygon that was the problem shape: `'POLYGON(( 118.860739531873 -20.2274797478397, 118.860739531873 -20.2274797478397, 118.860739531873 -20.2274797478397, 118.860739531873 -20.2274797478397 ))'` I am trying to understand how to fix this so my SSIS pckage runs and translates and loads my SQL tables without me having to fix up anything manually. Any help will be greatly appreciated. Thanks


Problem with Coordinate Tolerances when I export my .shp files to .mdb

![enter image description here](http://i.stack.imgur.com/GLRXh.jpg) When I export my shapefile to a personal geodatabase using ArcMap, the result is not the same as the original database. I do not know how to solve this problem.


vb.net arcobjects validate topology function only works after I use arcmap topology tool validate topology first

I find a strange thing occasionally, I have a validate topology function developed by vb.net arcobjects. It work not properly. The case is:   1. open arcmap   2. start editing   3. use my validate topology function, nothing happen, no error find (there really several errors)   4. set the extent to global and use arcmap topology toolbar 'validate topology in current extent' tool to validate topology   5. open error inspector and get errors list (4 errors)   6. change the current extent and use arcmap topology toolbar 'validate topology in current extent' tool to validate topology   7. open error inspector and get errors list (2 errors)   8. remain the same extent as step 6, use validate topology function developed by vb.net arcobjects, all errors displayed in error inspector. I want to know is it my code bug or operator issue? My code:                   ' Step 1 - access a topology     Dim topoUiD As UID = New UIDClass     topoUiD.Value = "esriEditorExt.TopologyExtension"     Dim m_application As IApplication = My.ArcMap.Application     Dim topologyExt As ITopologyExtension = CType(m_application.FindExtensionByCLSID(topoUiD), ITopologyExtension)  'original code from esri     Dim topology As ITopology          topology = CType(topologyExt.CurrentTopology, ITopology)     MsgBox("successfully access current topology")          ' Step 2 - Finding topology errors     Dim geoDS As IGeoDataset = CType(topology, IGeoDataset)     Dim errorContainer As IErrorFeatureContainer = CType(topology, IErrorFeatureContainer)     Dim eErrorFeat As IEnumTopologyErrorFeature     eErrorFeat = errorContainer.ErrorFeaturesByRuleType(geoDS.SpatialReference, esriTopologyRuleType.esriTRTLineCoveredByLineClass, Nothing, True, False)     Dim topoError As ITopologyErrorFeature          ' check the eErrorFeat length to determine whether has error     topoError = eErrorFeat.Next          topologyExt.ClearActiveErrors(esriTEEventHint.esriTENone)          topologyExt.DelayEvents(True)          Do While Not (topoError Is Nothing)         MsgBox("Topology error exists!")         topologyExt.AddActiveError(topoError, esriTEEventHint.esriTENone)         topoError = eErrorFeat.Next     Loop          topologyExt.DelayEvents(False)      Another thing is if I am not in editing session and use validate topology tool developed by vb.net arcobjects, the map will crash.


Select by Attribute --Top 500 values selection in attribute table - ArcGis

Is there any way we can select top 500 values in attribute table? I tried by sorting area field ascending/descending... but I can't count big numbers to select. In microsoft excell we can do know, similarly in arcgis, is there any code to keep in attribute selection ??


How to call ESRISignAddIn Silently (no UI)?

I am creating a build workflow and I ran into some problems. When I call **ESRISignAddIn.exe** from the command prompt it is throwing up a dialog saying "Add-In successfully signed". I was expecting it to report that back within the cmd dialog and not open up a new messagebox. Have anyone else seen this or have worked out how to sign an add-in "silently"? The /s switch does not work. This issue occurs with the 10.2 version of the exe (latest at time) and earlier versions.               "c:\Program Files (x86)\Common Files\ArcGIS\bin\ESRISignAddIn.exe" "C:\temp\my.esriAddIn" /c:"C:\temp\newCert.pfx" /p:myPass


Wildlife surveys using tablet and GPS, import/export to/from ArcMap

I've found some info about this on a previous post (Advanced android app for mapping and GPS logging), but it didn't provide enough info, I'm hoping some other tech experts can help! I have a tablet (Android) and it has integrated GPS, I would like to run a (preferably free) program into which I can input polygon layers (boundary units) and then track our flight path as we move through those polygon units. We are talking fairly large scale (units are at least 30 km²). There's no wifi, and we require something that keeps up to 100kph. Our current GPSs do this fine, but we can't import polygons into the GPSs. We have used ArcPad in the past, but it is too slow. To summarize - we want a tablet-based program to which we can import GIS polygon layers and then watch our track on the tablet as we fly transects in the polygon units. We would like to export these track to ArcMap when we're done, so pretty much any format will work. I've looked into a couple options (Ozi, Google Earth), but I'd like to know if anyone has any experience with tracking GPS on a tablet at higher speeds and across large areas. Thanks!


Trilateration algorithm for n amount of points

I need to find algorithm that can calculate centroid **A** (aka gravity center, geometric center, center of mass) from the figure where circles T1,T2,T3,T4,T5,..,Tn intersect AND length of line **R** from centroid to farthest corner of mentioned figure Following information is given:   * T1 Latitude = 56.999883 Longitude = 24.144473 Radius = 943   * T2 Latitude = 57.005352 Longitude = 24.151168 Radius = 857   * T3 Latitude = 57.005352 Longitude = 24.163356 Radius = 714   * T4 Latitude = 56.999042 Longitude = 24.168506 Radius = 714   * T5 Latitude = 56.994226 Longitude = 24.15709 Radius = 771 Result should look like this: A Latitude = XX.XXXXXXX Longitude = XX.XXXXXXX Radius = XX ![enter image description here](http://i.stack.imgur.com/kng2D.gif) As you probably already figured out, I am working on software that can find device location by closest Wifi Access Points or Mobile Base stations, as number of access points or base stations might change, I need an algorithm that can adapt to uncertain amount of points. There are some similar questions here and here, but none of them exactly answers to my question. Thanks!


Synching Two ArcMap JavaScript webmaps to respond to panning

I am looking for some advice on how to approach a problem I am having. I am looking to build an ArcGIS webmap using JavaScript and HTML 5 that reacts to a user panning on two separate, synchronous maps. Modest Maps provides an excellent example of Loading Two Maps and Synching Their Panning that helps visualize what I need to do. However, I need to accomplish the _Loading Two Maps and Synching Their Panning_ with ArcMap products. The Overview dijit in Arcmap seems to be capturing the essence of what I am trying accomplishing. However, the overview dijit in the widget overviewmap seems to react rather then pan synchronously. The dijit css feeds into this map and can be used as a resource if needed. I couldn't find anything in it that would accomplish synchronous panning. I am open to any ArcGIS suite methodology to accomplish the same as the Modest Maps look and feel. If you have any suggestions of comments I would greatly appreciate it.


KML file seems to load in Openlayers without error but doesn't show vectors

I'm trying to load a KML file into an Openlayers map. It seems to successfully run but I am not seeing anything show up on the map. Vectors should show up in an area of New Mexico if it's working correctly. Here is my code (here is a link to the KML file )               var map;          var baseMaps = [];     var sensors = [];          //create map and center over location     var init = function () {         // create map              map = new OpenLayers.Map({             div: "map",             theme: null,             projection:new OpenLayers.Projection("EPSG: 4326"),             controls: [                 new OpenLayers.Control.Attribution(),                 new OpenLayers.Control.TouchNavigation({                     dragPanOptions: {                         enableKinetic: true                     }                 }),                 new OpenLayers.Control.Zoom()             ],             layers: [                 new OpenLayers.Layer.OSM()             ]         });              map.setCenter(new OpenLayers.LonLat(-106.7653, 32.3197).transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject()),9);     };          //show kmlFile     var showSensor = function(sensorVO)     {         if(sensorVO.selected == true)         {             var kmllayer = new OpenLayers.Layer.Vector("KML", {                 strategies: [new OpenLayers.Strategy.Fixed()],                 protocol: new OpenLayers.Protocol.HTTP({                     url: sensorVO.url,                     format: new OpenLayers.Format.KML({                         extractStyles: true,                          extractAttributes: true,                         maxDepth: 2                     })                 })             });                  map.addLayer(kmllayer);             var select = new OpenLayers.Control.SelectFeature(kmllayer);             map.addControl(select);             select.activate();                    map.setCenter(new OpenLayers.LonLat(-106.639261,32.649908).transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject()),10);                  var obj = {"vo":sensorVO, "kmlLayer":kmllayer};             sensors.push(obj);         }              else         {                  for (var i = 0; i < sensors.length; i++)              {                 var vo = sensors[i].vo;                 if(vo.id == sensorVO.id)                 {                     map.removeLayer(vo.kmlLayer);                     sensors.splice(i, 1);                 }                  }              }          }          init();          setTimeout(function(){         var obj = {"id":1, "selected":true, "url":"kml-files/1858-Vegitation-Map.kml"};              showSensor(obj);     }, 1000);


How do I center Openlayers map over New Mexico

I'm trying to center open layers over NM and having a hard time to get it working. Here is my code. I'm trying to use decimal values and not the pixel values openlayers uses by default.               var init = function () {         // create map              map = new OpenLayers.Map({             div: "map",             theme: null,             projection:"EPSG:4326",             controls: [                 new OpenLayers.Control.Attribution(),                 new OpenLayers.Control.TouchNavigation({                     dragPanOptions: {                         enableKinetic: true                     }                 }),                 new OpenLayers.Control.Zoom()             ],             layers: [                 new OpenLayers.Layer.OSM()             ],             center: new OpenLayers.LonLat(105.9372, 35.6869),             zoom: 5         });          };


Legend labels for Vector layer symbols - how to name them (GeoEXT)

I'm coding with GeoExt and I have some problems trying to display the names of diferent features in a vector layer. This is my code:               mapPanel = new GeoExt.MapPanel({         map: map     });          legendPanel = new GeoExt.LegendPanel({         defaults: {                 labelCls: 'mylabel',                 style: 'padding:5px'                 },         bodyStyle: 'padding:5px; border:1px solid #8db2e3; border-top-width:0px;border-left-width:0px;border-right-width:0px',         width: 348,         autoScroll: true     });      In my legend I have 6 sublayers and all of them are titled with "Untitled" + 1,2,3,4,5, ... I would like to change the name of these different features, removing Untitled 1, 2, 3, ... Does anybody can help me to do that? I think it is easy but I don't know how. Thanks in advance! Rafael


How to modify two polygons which have common point and line at the one time

I would like to modify two polygon at the one time because these two data has same point and line. Is it possible?? In the present circumstance, I have to modify polygon one by one. First, I modify one polygon. After that,another polygon will be modified with snap function for overlapping to another polygon. Thank you so much in advance.


Postal areas for Germany? - Shape?

does anyone of you know how and where I can download postal area districts of Germany in shape files (download free of charge if possible). best regards desputin


Changing CRS crashes QGIS every time

When I open QGIS, and use the tools, setting -> project properties, QGIS will crash. When I also use the setting -> option -> CRS , Select new CRS also causes a crash. It looks like couldn't use the Coordinate Reference System Selector. Can I fix it by myself?


GDAL: Convert specific RGB (GeoTIFF) color to transparent

I have an RGB GeoTIFF file and I want to convert a specific RGB colour to transparent. I would like to use GDAL for doing this, but external Python scripts are also an option (though I haven't learned the art of Python scripting yet). My thoughts are:   1. Add a 4th band to my GeoTIFF file and fill with 255.    2. Insert 0 into the 4th band where RGB colour == _target colour_ The first part I think I can solve using               gdalwarp -dstalpha source.tif target.tif      The second though I don't know how to solve. **matt wilkie** suggest using a Python script he wrote, but I think that is done separately for _each band_ and not for a _specific colour_. **How can I assign a value to the 4th band of my TIFF when the RGB colour is equal to a specific RGB colour?**


ArcGIS problem accessing fields attribute table

I encounter a strange problem. I have a simple shape file to which I added some fields and then populated them. Now I wanna access those fields from a script to run some buffer operations. I wanna do this per line segment so I execute a `Select_analysis`. The tool runs without throwing an exception but not all fields are present in the newly generated shape file, namely only those from **before** I added the fields. For testing purposes I tried to access the fields also with a `SearchCursor`-method, which returns `none`. Strange enough when the shape is used within ArcGIS directly everything works fine (`Select_analysis`, calculating on the newly added fields,...). Any idea what is going on? Every hint is very much appreciated! EDIT: I'm using ArcGIS 9.3.1. on WindowsXP. I believe posting the script isn't of much help as the problem seems to be with the ESRI software. But I do it anyway, maybe someone does spot something strange.                # ---------------------------------------------------------------------------     # zonering_leggerV2.py     # docu still to come...     # ---------------------------------------------------------------------------          # Import system modules     import sys, string, os     from GeoSql import ArcSql     import GIS_tools as gt          # Create the Geoprocessor object     # Geoprocessor-object (ArcMap 9.1 and older)     try:         import win32com.client         gp = win32com.client.Dispatch('esriGeoprocessing.GpDispatch.1')          # Geoprocessor-object (ArcMap 9.2 and newer)     except:         import arcgisscripting         gp = arcgisscripting.create(9.3)               try:         # global gp stettings         gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx")         gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Analysis Tools.tbx")         gp.OverWriteOutput = 1         gp.CheckOutExtension("Spatial")         msg = gp.getMessage(2)         gp.addmessage(msg)          except:         msg = gp.getmessages()         print msg         gp.addmessage(msg)                    def buff(source, dest, field_name_zone, side, type_end):              # Process: Buffer (2)...         try:             gp.Buffer_analysis(source, dest, field_name_zone, side, type_end, "NONE", "")             msg = gp.getmessage(2)             print msg, '\n'         except:             msg = gp.getmessages()             print msg, '\n'                    def merge(zone_l, zone_r, out_data):         # Process: Merge...         merge_str = "%s;%s" % (zone_l, zone_r)         try:             gp.Merge_management(merge_str, out_data)             msg = gp.getmessage(s)             print msg, '\n'         except:             msg = gp.getmessages()             print msg, '\n'                    #------------------------------------------------------------------------------------------------------------     # --> programme parameters          type_end = 'FLAT'     field_name_zone = {'KR': 'KERN_R',                       'KL':'KERN_L',                       'BR':'BINNEN_R',                       'BL': 'BINNEN_L',                       'f_name_zone_new': 'code_zone',                       'KERN': "_KZ",                       'B_BINNENKANT': "_BBI",                       'B_BUITENKANT': "_BBU"}          #------------------------------------------------------------------------------------------------------------     #--> file declarations and data gathering          dir_base = "D:\\opdrachten\\legger\\dijkvakken\\"         dir_gp_tmp = dir_base + "gp_tmp\\"          vakken_shp = "%sbijlmerring_042012.shp" % dir_base          column = "PKM_NAME"     #vakken = gt.get_data_from_table(vakken_shp, fields = [column])     vakken = ["A2035-001"]     invalid =('-')          for vak_name in vakken:         if invalid in vak_name:             vak_name_mutat = vak_name.replace(invalid, '_')                   selected_vak = "%s%s.shp" % (dir_gp_tmp, vak_name_mutat)              # Local variables...         buffer_KL = "%s%s_%s.shp" % (dir_gp_tmp, vak_name_mutat, field_name_zone['KL'])         buffer_KR = "%s%s_%s.shp" % (dir_gp_tmp, vak_name_mutat, field_name_zone['KR'])         buffer_BR = "%s%s_%s.shp" % (dir_gp_tmp, vak_name_mutat, field_name_zone['BR'])         buffer_BL = "%s%s_%s.shp" % (dir_gp_tmp, vak_name_mutat, field_name_zone['BL'])         #buffer_shapes = [buffer_KL, buffer_KR, buffer_BR, buffer_KL]                   dat_zones ={'kern_l': (buffer_KL, field_name_zone['KL'], "LEFT"),                    'kern_r': (buffer_KR, field_name_zone['KR'], "RIGHT"),                    'binnen_l': (buffer_BL, field_name_zone['BL'], "LEFT"),                    'binnen_r': (buffer_BR, field_name_zone['BR'], "RIGHT")}              merge_kern = "%s%s_mergeKern.shp" % (dir_gp_tmp, vak_name_mutat)         merge_binnen = "%s%s_mergeBinnen.shp" % (dir_gp_tmp, vak_name_mutat)              dissolve_kern = "%s%s_dissolveKern.shp" % (dir_gp_tmp, vak_name_mutat)              erase_out = "%s%s_erase.shp" % (dir_gp_tmp, vak_name_mutat)              result_tmp = "%s%s_zones.shp" % (dir_gp_tmp, vak_name_mutat)              sql = ArcSql(column, [vak_name]).statement()         print sql              # Process: Select...         try:             gp.Select_analysis(vakken_shp, selected_vak, sql)             msg =gp.getmessage(2)             print msg, '\n'         except:             msg =gp.getmessages()             print msg                        for i in dat_zones.itervalues():             buff(selected_vak, i[0], i[1], i[2], type_end)              merge(dat_zones['kern_l'][0], dat_zones['kern_r'][0], merge_kern)              try:             gp.Dissolve_management (merge_kern, dissolve_kern, column)             msg =gp.getmessage(2)             print msg, '\n'         except:             msg = gp.getmessages()             print msg, '\n'              gt.add_field(dissolve_kern, field_name_zone['f_name_zone_new'], 'TEXT', length_txt='10')         gt.update_fields(dissolve_kern, sql, [field_name_zone['f_name_zone_new']], [field_name_zone['KERN']])              for f in [dat_zones['binnen_l'][0], dat_zones['binnen_r'][0]]:             gt.add_field(f, field_name_zone['f_name_zone_new'], 'TEXT', length_txt='10')             if '_L' in f:                  gt.update_fields(f, sql, [field_name_zone['f_name_zone_new']], [field_name_zone['B_BINNENKANT']])             elif '_R' in f:                 gt.update_fields(f, sql, [field_name_zone['f_name_zone_new']], [field_name_zone['B_BUITENKANT']])             else:                 'error at update field binnenzone...!\n'                 sys.exit()                   merge(dat_zones['binnen_l'][0], dat_zones['binnen_r'][0], merge_binnen)               try:             gp.erase(merge_binnen, dissolve_kern, erase_out)             msg = gp.getmessage(2)             print msg, '\n'         except:             msg = gp.getmessages()             print msg, '\n'              merge(erase_out, dissolve_kern, result_tmp)


How to find the two furthest points between a group of points?

I have a table with 20000 records and with those columns Shop, longitude, latitude, geom I would like to find a way to find the pair of ths farthest shops between a selection of shops For example we have outletnames : A,B,C,D,E,F,G,H,I,J,K,L,M,N,O... I would like to know what is the bigest distance between all those shops A B C D and E A-B = 4 km A-C- = 2 km A-D = 1 km A-E = 1,3 km B-C = 20 km C-D = 50 km .... I would like to have an optimised query that would give me as a result 50 km Thank you


Is it possible to make a routable graph from Polygons?

I would like to know if it is possible to make a graph with pgrouting 1.05 on polygons. What I mean is that I have to make a graph on a database without roads but with polygon geometry. I would like to make a graph on this geometry. Thanks for your help and I'm sorry for my english


QGIS raster clip error: Cutline not of polygon type

I have a raster file and want to clip it, but I'm getting this error: "ERROR: Cutline not of polygon type" what does this error mean and how can I resolve it?


How do I edit joined parcels to eliminate gaps in a parcel fabric?

I have a parcel fabric layer that has a gap in it: ![enter image description here](http://i.stack.imgur.com/KluAC.jpg) I want to edit the selected parcel on the bottom to fill the gap. I don't want to construct a new parcel to do so. How do I modify the parcel in such a way that it maintains its attributes and fills in the gap with the parcel above it? * * * **EDIT:** Based on answers and comments, I created a topology and placed the parcel fabric inside it. It appears that you cannot use topological tools to edit a parcel fabric. As such, any answers using topology will not work.


Is it possible to create Thiessen polygons around nodes in QGIS? e.g. http://en.wikipedia.org/wiki/Voronoi_diagram


change default font of labels in QGIS

I've just opened a feature request, but maybe someone knows a temporary solution. In QGIS 2.2 one can change the default font of the Print Composer (through the Menu Option -> Setttings). But I think it is not possible to do the same with labels. This is a little bit boring when one has to export a project from a Linux to a Windows machine. In Linux (debian in my case) the default font in Cantarell, font that is not recognized by Windows. This cause a boring error message. No big deal, just annoying. Thanks


Get a list of property elements using GeoServer

I'm using GeoServer with openlayers and I want to build a filter check-box dynamically, depending on the elements of this column: ![Column](http://i.stack.imgur.com/Isl8G.png) How can I query GeoServer so that I get a list of the different elements in this column? Or at least all the information of all the columns in XML so that I can parse it?


QGis export part of map (box with coordinates) as raster

I would like to export a map of Switzerland (as part of googlemaps) to use as a background for a graph, but I need to know exactly where each corner of the map corresponds to (wrt spatial coordinates). Is there a plugin to do this in QGIS? Ideally I would like to be able to manually select a bounding box (while using the OpenLayers underneath) and get a picture exported together with the bounding box coordinates. Thanks


GWR Model with Zero Intercept

In GWR of ArcGIS, if we get negative intercepts that give negative predicted values of dependent variable that are undesirable. Then , is there any possibility to make intercept zero?


style creation in QGIS

![style model dash double dot dash](http://i.stack.imgur.com/jAZtS.jpg) I want to create a line style in QGIS (model is attached here with), Any body can help me?, please help me.


Are there more events in OpenLayers 3, besides "moveend"?

I am building an API on top of OpenLayers 3. I can't display here the whole code for professional reasons. This is a bit how my classes look like.               // external file xgis.js          xgis = {};               // xgis.map     xgis.map = function(options) {         // map definitions ...     };     // Inherits from ol.Map     goog.inherits(xgis.map, ol.Map);      The goal is to build an SDK, which I named here "xgis". We want to implement some callbacks. For example, when I perform a zoomIn, I have a callback named onZoomIn that executes some code. I know that it is done with the event "moveend", like in this fiddle. Now, I want to implement some callback functions:   * onZoomIn --> fires after a zoom in   * onZoomOut --> fires after a zoom out   * onMapReady --> fires when the map is ready (has loaded all the layers)   * onPanEnd --> fires when the drag/pan has ended Are there other events in OpenLayers 3, like "zoomend" or something? Thanks


ArcGIS Server Error: "Feature service requires a registered database"

I'm attempting to publish a **feature service** to ArcGIS server. I'm running SQL Server Express and ArcSDE. When I try to publish the service I get the error "00090: Feature service requires a registered database". ![enter image description here](http://i.stack.imgur.com/pJaOt.png) When I validate the database registration it appears to succeed, but I still can't publish the service. The features I am publishing are located inside of the geodatabase that I am registering. This geodatabase is located on the same server that ArcGIS Server resides, but I am attempting to publish from another machine. Am I missing some crucial step? Here is a bigger version of the screenshot.


Problem with QGIS 2.0 installation from OSGeo4W set up

I'm trying to install QGIS 2.0 on Windows XP but when the download begins, suddenly OSGeo4W show me this message:                Unable to extract/2       #define FT_SUBLYPH_FLAG_ROUND_XY_TO_GRID       #define FT_SUBGLYPH_FLAG_SCALE      -- the file is already in use.       Please stop all OSGeo4W processes and select "Retry", or select "Continue" to go on anyway       (You will need to reboot)      Nothing happens neither if I choose retry nor if I choose continue, and I must to shut down the installation Hint? :-)


Join attributes by date in Qgis

I learned (thanks!) to use the join feature built into the layer properties, but in my case it doesn't work because data do not match perfectly. Infact shape refers to geotagged photos, csv is a log file which contains orientation data (azimuth); I would like to join them by date field (day, hours, seconds): I would need to join the "nearest" data, non spatially but temporally!


Shape length distance wrong in Robinson projection

I have a world map in the WGS 1984 Geographical Coordinate System. I want to measure the lengths of borders between countries in kilometers. To get the units in meters, I have first used the Define Projection tool and then the Project tool to change to the Robinson Projected Coordinate System. After that, I am using the Intersection tool and setting the output type to "line" to get a Shape_length column where I will have all the borders and their lengths. But, the value in the Shape_length column appears to be wrong as it measures the borders to be only a few meters. The Measure tool does the same thing. I can't get my head around what I am doing wrong. I've attached a print screen of the data frame. What could the problem be? ![enter image description here](http://i.stack.imgur.com/8G6c8.jpg)


Merge/mosaic categorical rasters in QGIS/GDAL

I am attempting to make yearly bins of some monthly MODIS data for all of the US. The cells of the raster are defined with Julian dates indicating the date at which a burn event took place. After merging the data using the gdal script listed below, I have full coverage of the US but I suspect that some of the data was lost. Being new to QGIS and GDAL, I haven't been able to figure this one out yet, but in ArcGIS there was an option when mosaiking to keep the (max, min, mean, blend, etc) value of overlapping cells. With this data, a max (date) should be the appropriate answer as long as it falls within the julian range (1-366). After figuring out how to get raster cell counts in QGIS, it looks like just the last couple of months are represented in the data and the cells with low Julian dates have been washed out or lost. (scroll over on code)               for YR in $YEAR; do     gdal_merge.py -o MCD45yr_${YR}_burndate_2.img -of 'HFA' -a_nodata '-32768' -n '900' -n '9998' -n '9999' -n '10000'  MCD45monthly.A${YR}???.Win03.051.burndate.tif;      -32,768 is the NoData value, and the other values listed with `-n` are data anomalies. Any suggestions for this issue? I would expect the majority of the Julian dates to have at least a couple of cells with data in that range.


Why are Data Projections relevant?

Sorry for my ignorance - especially after being in this area for a long time. Why is it useful to store data in unusual projections and datums? I understand the value of projections as **output** constructs, because of distortion, etc: that much is justifiable. However, I don't understand why, for instance, states use state plane projections in data. We have decimal accuracy: why not just store precise values in EPSG:4326? Is this entirely a remnant of pre-auto-reprojection days, or is there a value proposition I'm missing? edit: I'll restrict the scope of this question to **vector data only** to make it more concrete. Thanks whuber for the tip.


Shapefile attributes to Autocad Annotation

How can I export the attribute information for the shapefile/feature class to Autocad Annotation/Text. I'm able to export the shapefile to autocad but unable to get the attributes in annotation. Need you help its really urgent. Also Let me know if there any any way we can append attributes in autocad without using blocks. Thanks and have a great day. regards, Ramakrishna Billakanti


Tilecache + Mapnik + Apache generates incomplete tiles (missing layers)

I have TileCache + Mapnik working with a custom osm.xml stylesheet. When I run tilecache_http_server.py the tiles are served OK. But when I run TileCache under Apache (with mod_wsgi) I see incomplete tiles (with missing layers). What can go wrong? Furthermore, generating tiles is very slow. Can this have something to do with the missing layers? Thanks.


ArcMap crashes on TabbedFeatureInspector sample

The TabbedFeatureInspector, (VB.NET), sample crashes ArcMap when I do the following:   1. Open ArcMap.   2. Open map document with feature class with class extension loaded.   3. Select a feature.   4. Start an edit session.   5. Click the Attributes button to open the Attribute editor.   6. Click the Create Features tab. At this point the Create Features tab displays as a blank tab and ArcMap locks up. The wierd thing is that this only happens with the initial selection of a feature and opening of the Attribute editor. If you open the attribute editor, and then select a different feature before selecting the Create Feature tab ArcMap doesn't crash. I initially found this by building a tabbed inspector into an application using the sample as the template. After discovering the problem I tried it in the TabbedFeatureInspector sample and found the same result. The offending code appears to the be SetParent line in Inspect event where the inspector parent window is set to the groupbox on the user control.               Public Sub Inspect(ByVal objects As ESRI.ArcGIS.Editor.IEnumRow, ByVal Editor As ESRI.ArcGIS.Editor.IEditor) Implements ESRI.ArcGIS.Editor.IObjectInspector.Inspect              SetParent(Me.m_inspector.HWND, Me.GroupBox1.Handle.ToInt32)         ShowWindow(Me.m_inspector.HWND, SW_SHOW)      Any suggestions on how I can correct this?


Is it possible to reveal the relief parameters “slope inclination” and “slope exposure” by dint of ASTER GDEM?

I’m trying to find information in technical and scientific articles – is it possible to reveal the relief parameters “slope inclination” and “slope exposure” by dint of ASTER GDEM? Or we can define from it the only parameter of relief – the height above sea level? I can’t find any articles that confirm this possibility, nor reject it. If it’s possible, could you tell me the way to do it or the direction from which I can start searching for this information?


How to express not within/inside predicate in spatial database with good performance

In our application we want to express some spatial predicates as objects not within some areas, but unfortunately it has very bad performance. The following SQL is generated by our application:               SELECT *     FROM geom_table     WHERE NOT (SDO_INSIDE(geom_column, SDO_GEOMETRY(...), srid) = 'True'     AND geom_column IS NOT NULL); I've been searching on internet for a while but surprisingly fail to find any information talking about it. I think it's not a non-reality use case, or is it? I wonder if there is a better way to express the SQL statement?


QGIS labels are different on canvas, composer and from PDF

I am running QGIS v2.2 off a 64 bit Windows server. I have created a map with a number of local government area polygons, road polygons (to which I have applied a Layer Source to limit the labelling to major roads), and some town point objects. The polygons and points labels are OK, however, the road labels appear differently on the Composer and different again when I export to PDF or as an image. Some labels are missing or repeated. Is there something wrong with my labelling? I am using Arial 4.25 with a 0.5mm buffer, curved and centred on the line, inside/outside angle for curved characters is 29, High priority, showing all labels, merging connected lines to avoid duplicates. Even if there is a problem with my labels, I don't understand why they would appear differently between the three views of the same map.


Using GeoBrain WCS4DEM in QGIS?

I added the GeoBrain WCS server to the server list:   http://geobrain.laits.gmu.edu/cgi-bin/gbwcs-dem?version=1.0.0& I connect to the server and get a list of layers. I select the first one listed - SRTM 30MUSA and select ADD. My map window is displaying a region in NM, the project CRS is EPSG:4326. The response is ![as here](http://i.stack.imgur.com/xXMGi.png) Note that the bbox is nonsense. Any idea what I'm doing wrong? Edit: This question is closely related, asking how to use a clip region to set the bbox in a WCS request, and the example is also GeoBrain.


Using ST_MakeEnvelope over the dateline

I'm using ST_MakeEnvelope in PostGIS to try and make geometries that cross the international dateline, but they are making a shape that wraps the other way around the globe. The geometry I want is one which covers the eastern half of Australia and the dateline. It looks something like:               ST_MakeEnvelope(133, -49, -170, 0, 4283)      But:               ST_Intersects(ST_SetSRID(ST_MakePoint(150, -25), 4283), ST_MakeEnvelope(133, -49, -170, 0, 4283));      returns false. The shape goes the wrong way around the world! If I just take:               ST_Intersects(ST_SetSRID(ST_MakePoint(150, -25), 4283), ST_MakeEnvelope(133, -49, 180,v0, 4283));      I get true, as expected. How can I set a direction to my envelope? Or should I create a polygon some other way? (If you're wondering: I'm trying to select geometries in a table for a map view, and when the view includes the dateline, nothing is shown on the map) [ETA] If it helps, the box actually arrives in Spherical Mercator, looking something like:               14855319.859171, -6294627.0589708, 21149946.918142, 0      But calling ST_Transform:               st_transform(st_makeenvelope(14855319.859171, -6294627.0589708, 21149946.918142, 0, 900913), 4283)      Results in geometries with the negative-valued longitudes:               SRID=4283;POLYGON((133.447608799776 -49.1163523261112,133.447608799776 0,-170.006794251165 0,-170.006794251165 -49.1163523261112,133.447608799776 -49.1163523261112))


How to convert x y coordinates to longitude latitude using QuantumGIS?

I'm working with this downloadable shapefile that has X Y coordinates. They are akin to `'3672187.92698000, 534175.72095400'`. I would like to convert them to longitude latitude so they are more like `'-90.097017, 29.963176'`. I've seen this question tackled using ArcMap however I don't have that software. I was able to download and install Quantum GIS but I am unfortunately perplexed by its complicated interface. Would like to do the conversion with it, if possible.


Mapserver dynamically line

I'm a beginner for GIS especially for Web GIS that's implementation to mapsever. I created a project that I've uploaded at dedyfuadi.web.id . For the map features, I made by QGIS 2.0 and for framework I choose pmapper for displaying map dynamically. But, I have some problem for create line between 2 points directly by syntax (I wont manually adding the lines from QGIS). What should I do?


Is there any instructions to download OSM in postgres database? I use geosuite (geoserver, postgres, postgis). I'd like to download USA and Russia (1-12 levels). I also want to download sld styles. (Is there any standard way). So, there is intranet. We use geosuite. And I'd like to download OSM with sld styles.


How do I add feature from WMS to a different layer in OpenLayers?

I have a WMS layer I serve from SQL View. When I GetFeatureInfo I receive JSON presentation of the feature.               {"type":"FeatureCollection","features":[{"type":"Feature","id":"Trips.fid--3a106d6c_13f7fcd4d00_-7dab","geometry":{"type":"MultiLineString","coordinates":[[[34.78731593803529,31.31823531417712],[34.78843141500525,31.336622682320776],[34.79414493381831,31.355201806839506],[34.794207563998924,31.371163684259038],[34.79296734448457,31.375581921229493],[34.78922614608776,31.390623411898762],[34.78827495215694,31.408575008023988],[34.78636890354401,31.42299930464733],[34.78632872060815,31.42193697540446],[34.78629782738711,31.422105913663092],[34.78628031267043,31.421992838029546],[34.786190597348,31.42158272671166],[34.78616186179376,31.421905927160644],[34.784898002487694,31.424062281290652],[34.773401914392345,31.434990008872095],[34.76559077292494,31.446642271859293],[34.77342394487551,31.466799593779264],[34.77975547041099,31.490976936335215],[34.78263023215561,31.506823501958817],[34.792679583505226,31.531054387260554],[34.80048435469499,31.55136414506193],[34.80654617304027,31.574946250797954],[34.813199698392914,31.59725795152748],[34.81470088275501,31.619577056013927],[34.809541545607935,31.657018778594132],[34.80638283046715,31.674769044083895],[34.81569739824674,31.69878758883678],[34.82117079659263,31.721351512377144],[34.822676820645704,31.74065914102294],[34.82640625397537,31.760366608188313],[34.83277797643619,31.795874550345438],[34.84013377052116,31.80349376296454],[34.85764369385628,31.809682398342275],[34.87655834935133,31.809001817509152],[34.89977578898199,31.81087415167135],[34.91555026001387,31.814641705371283],[34.92078885084598,31.819327945165384],[34.94872743329605,31.826006370096646],[34.96865364149424,31.831077555261356],[34.977436258261335,31.8348796713872],[34.988078880595516,31.834460104995458],[34.99298330536876,31.826840136006965],[35.00026754690508,31.826135241324245],[35.005568854286416,31.82592323790467],[35.01774054730234,31.8215162498958],[35.02155196628523,31.818288829834177],[35.02691352351924,31.813341112120234],[35.03143746024918,31.80979583988096],[35.03815953441021,31.80619724694477],[35.045420139243944,31.80474599411263],[35.05317138181679,31.804966181835955],[35.06486714082564,31.804141661709476],[35.07584180677903,31.800281337346696],[35.08838316507283,31.800341312456926],[35.10640165800501,31.801455597207628],[35.12610124373968,31.800681742893538],[35.14697634526269,31.79722164231699],[35.16219627795028,31.79329681701213],[35.171594439546396,31.797658570113885],[35.187671634314036,31.797928602082614],[35.188362774999014,31.79653341639412],[35.1887256549307,31.7961346714363],[35.19047234637513,31.794962708212374],[35.19256271051755,31.79461237466848],[35.19270785200568,31.79363734891888],[35.193583458825614,31.79312975679172],[35.194747974527196,31.793406487587923],[35.196095096006964,31.79164621288934],[35.196955729849186,31.790773244292048],[35.1978874456438,31.790160915552967],[35.198711620213004,31.78987941300083],[35.19929508848269,31.787229378909604],[35.19850752277577,31.785375424789216],[35.19941639636708,31.784687679081976],[35.199241688426945,31.78231397727036],[35.198391385845284,31.781402127099458],[35.201245278823734,31.77605635342819]]]},"geometry_name":"geom","properties":{"Id":1,"StartDateTime":"2013-06-23T01:49:57Z","EndDateTime":"2013-06-23T03:23:32Z","Length":94.57}},{"type":"Feature","id":"Trips.fid--3a106d6c_13f7fcd4d00_-7daa","geometry":{"type":"MultiLineString","coordinates":[[[35.20076873304309,31.77473722793188],[35.204546738625524,31.774719539653546],[35.208965075961345,31.77618510615382],[35.208640148891355,31.78202931583896],[35.20975512430616,31.78246578222451],[35.21137792550667,31.781691930524566],[35.213549710915274,31.780813498217512],[35.21476900302076,31.78168291327713],[35.21469230369844,31.780946053301573],[35.21257227234834,31.781205360240914],[35.20864390470302,31.783054725633296],[35.20829895121278,31.78302020473956],[35.203248270523765,31.78222448311647],[35.19854766107952,31.786334626536405],[35.19984531073088,31.78751601702663],[35.20014413481394,31.788041155898213],[35.19725069211512,31.790485693913038],[35.19110776321938,31.795067815242238],[35.18161040176461,31.801656830805612],[35.1687011405429,31.79567693197168],[35.159244819963,31.795794581845737],[35.11942475024883,31.798775201079778],[35.08039358590231,31.800856179659302],[35.063209196166355,31.805088407881986],[35.04805678839844,31.804031639631894],[35.032927570632985,31.80847877245744],[35.02092427914747,31.81877639146],[34.99012978758847,31.829259174717787],[34.98804795187599,31.839040689444193],[34.977978798901404,31.835101996624207],[34.9498888706375,31.82596783070529],[34.92736197558343,31.823111392882268],[34.91695636557349,31.815625240675068],[34.903688415937665,31.811906632219383],[34.88296941318597,31.810547012219338],[34.86002652341891,31.81005225691859],[34.8400240232646,31.803567891209664],[34.83015973660191,31.79428193887949],[34.832461959577664,31.781116935605244],[34.82711462326763,31.762582638502867],[34.82278272110514,31.743736272787874],[34.82054030348634,31.720919484741522],[34.81420911419586,31.695453086060326],[34.80627534157327,31.675224452093136],[34.81001464412418,31.653220249705704],[34.814807584458386,31.623882775804965],[34.81329122844819,31.60270186644005],[34.8080051217612,31.582810947810554],[34.80303812730152,31.560726326313713],[34.797198877391246,31.54151166321618],[34.786754234512,31.52322945275979],[34.78295267962934,31.504485052320216],[34.77614342612191,31.48399237576942],[34.77342018413583,31.467251749582854],[34.76645583238203,31.450437751509558],[34.76590528960303,31.442231443933725],[34.78525964273001,31.422175881020067],[34.78907352850801,31.401593669076053],[34.79129725926219,31.38162192285955],[34.795717357038875,31.365271563573277],[34.792079987685526,31.350441990935515],[34.786453724639145,31.329829087974357],[34.78839683025119,31.309436171497104],[34.79081076415792,31.289914211434404],[34.7976451665849,31.28794303366183],[34.80366544178834,31.286754464613146]]]},"geometry_name":"geom","properties":{"Id":2,"StartDateTime":"2013-06-23T04:46:14Z","EndDateTime":"2013-06-23T06:03:47Z","Length":102.14}}],"crs":{"type":"EPSG","properties":{"code":"4326"}}}      How can I add this item in OpenLayers to new vector layer (to simulate SELECTED style)?


Full extent of a map

Is there a way to get the full extent of a map? I am looking to use this extent to zoom out to whole world view on a button click.