title: Draw a line and a buffer around it onto the spheroid

text: I need to draw a line between two points and create a buffer around this line. The only difficulty is that I need to draw not a "cartesian" line but a line that is laying on the spheroid. I tried to use ST_MakeLine but it has only geometry parameters so it draws only a straight line, not a curve (which is a straight line on spheroid). Unfortunately ST_LineFromText also does the same. Here is an example from http://workshops.opengeo.org/postgis- intro/geography.html. ![](http://workshops.opengeo.org/postgis- intro/_images/lax_cdg.jpg) I want to draw the red curve. (Btw the length of this curve can be calculated by the ST_Distance(geography, geography)) Thank you! Peter


title: How to do separability analysis in eCognition?

text: Can I use eCognition to check which bands are more helpful to separate individual classes and to get output result? If so, how can I do this?


title: Assertion Error when attempting arcpy.mapping.AddLayer

text: Whenever I attempt to use AddLayer it gives me: Runtime error <type 'exceptions.AssertionError'> This is what I've been using:               arcpy.mapping.AddLayer(MainDF, r"\\GEOGDATA\campusmapping\CampusWide\LayerFiles\CAD\Floor2\Casanova2", "BOTTOM")      Casanova2 is a .lyr file. I've also attempted to use this with a few other layer files. Every time it gives me an AssertionError. Do you know what I'm doing wrong?


title: Turn Table To Turn feature class

text: I will convert a Turn.dbf file with 'Turn Table To Turn feature class' Toolbox so I can make me one way online. But I get a warning: Edge features are not connected. Turn OID = 8 Edge features are not connected. Turn OID = 8 "C: \ DOCUME ~ 1 \ Tae \ LOCAL ~ 1 \ Temp \ turn4.fid" contains the full list of turns not comfortable to be imported. Does anyone know what this comes from?


title: Can Google Earth App for the iPad specify a KML using a URL without dealing with Google MyMaps?

text: Can Google Earth App for the iPad specify a KML using a URL without having to go through with one's Google My Maps Account?


title: Coordinates of upper left corner in the world file?

text: According to Wikipedia and ESRI last two lines in the world file represent x and y coordinates of **center** of upper left pixel. Now, if I have a world file like this: 32.0   0.0   0.0   -32.0   691200.0   4576000.0 ...can I say that the "real" upper left corner of the image (upper left corner of the upper left pixel) is: X = 691200.0 - (32.0 / 2) = 691184.0   Y = 4576000.0 - (-32 / 2) = 4576016.0 ...and should I use these coordinates if I want to use the image with Bing Maps or Google Maps (after converting to Lat/Long off course), or should I stick with original ones (pixel center coordinates). Thanks, Anil **UPDATE:** I just realized that the .tab file associated with the map has coordinates of the upper left pixel (0,0) calculated the same way as in my question above. Not sure about the exact formula but resulting coordinates are the same - shifted left for half of the pixel width and up for half of the pixel height. Now I'm pretty confident that this is the right way to do it.


title: Open Geospatial Consortium (OGC) standards

text: I was trying to understand about the OGC and OGC Standards but i got confused that are they provide standards like ieee or they provide standards that should be implemented with your application i am looking for some one to give me details about the ogc,It is hard to understand it on there site


title: How to add buttons to QT ui for qgis feature forms to loop through attributes?

text: I'm looking for a bit of guidance. I'm happy creating a standard feature form using QT designer (as excellently documented here) for QGIS but I was wondering how you add a forward and back button to a feature form, so that it would loop through all the identified features in QGIS? Below is a mockup of what I'd like to achieve but need some help. You'll see that I haven't added in any 'line edits' in my tabs and that was purely to produce a mockup for the purposes of this question. ![enter image description here](http://i.stack.imgur.com/xR1Dm.png)


title: How to use ArcGIS functionality in Python without running ArcMap (e.g. separate editor/environment)

text: I'm trying to take the next step in programming ArcGIS with Python and see if its actually possible to run my script/code without having to load ArcMap and paste them into the Python window. I imagine I would do this view some sort of IDE specific to Python that loads something from ArcGIS. How would I got about doing this?    **Update: Summary of the solution** I went the Eclipse + PyDev route. Here are my notes on installing everything and getting it working: **Download Eclipse** This is slightly confusing given all the different flavors that are available (Java, Java EE, C/C++, etc.). I went with C/C++ since that's what familiar to me. It doesn't really matter since PyDev gets installed after this step. Pick what you want. **Install PyDev** PyDev will be installed via Eclipse. The instructions linked worked perfectly. **Configure PyDev** The above question gives a screenshot of a Windows configuration, while this link from an ESRI blog has more info for Unix systems although it is from 2008. **Ran a test script** Fired up a script I had been working on and it worked. There seems to be some differences in executing code in this way in how to access shapefiles and the workspace, etc, more explicitness is required. More to learn here.


title: ArcMap v10 - Polygon Layer - Keeping one polygon above another

text: I have a shapefile with polygons of differing colors. I am using a white hash pattern to represent partial ownership. In short - Polygon #1 is Green and Polygon #2 is a White Hash pattern. **Problem:** The Hash Pattern is not staying on "top" of the Color Pattern. What do I need to do to keep the hash pattern on "top"? Thanks


title: How to track logged in user in geoserver?

text: Is there any way I can track the logged in users in geoserver(somewhere in log file). I have same users in geoserver and tomcat-realm.As I can see what URI are requested, IP address URI requested from, time and other info in access log of apache. But I need to find user name who logged in to request those URI. Thanks


title: Remove Holes Existing Within a Single Polygon in ArcGIS 10

text: I have a shapefile that consists of a single polygon. I have gaps within the polygon that need to be eliminated. The resultant dataset shoulde be a single, solid polygon. I have tried implementing topology rules (rule: must not have gaps) on the dataset (it's in a geodatabase and within a feature dataset), but I don't think this is the appropriate solution/technique to accomplish this task. Any ideas? I'm on ArGIS 10 with an ArcInfo License.


title: Printing maps from OpenLayers with a Google Maps base layer

text: I need to implement an automatic pdf report generator on a web site. In the generated pdf file I need to insert images from maps extracted from OpenLayers with a Google Map base layer. On top of Google Map I have an overlay with some custom OWS (WMS/WFS etc). Is this possible? And is it allowed? I was reading the ToS of Google, and it is not totally clear to me.


title: How to run shapefile.py

text: I'm script challenged and trying to repair a shapefile with a missing .shx. This tool sounds ideal: https://code.google.com/p/pyshp/, but the four step instructions are not enough. Can someone give me a quick walk through how you would run this tool.


title: OSGeo4W and QGIS 2.0

text: I just ran the OSGeo4W installer to update my previous installation and I noticed that it was updating my QGIS installation from 1.8 to 2.0. I haven't been able to find any information on the web about QGIS 2.0 being released yet. The installer ran without any issues, but now that it is complete I cannot find executables for QGIS. The folder for QGIS exists within the `apps` folder under `C:/OSGeo4W/`, but there are no executables in there. I've already completely uninstalled all OSGeo4W applications and reran the installer using the Express Install option, but I had the same result. How can I use QGIS and has version 2.0 been released yet?


title: Proximity in space and time

text: I have some point data which represent daily lat-lon locations of an animal, with an associated timestamp. I would like to identify all points where STATIONARY = TRUE. A point qualifies as being stationary if a 100km buffer around it overlaps an additional (say) 5 **temporally adjacent** points. So if _day 10_ is my point of interest, I want to ask whether 5 temporally adjacent days are within a 100km buffer of this point. If days 5,6,7,8 & 9; OR days 11,12,13,14 & 15; OR days 8,9,11,12,13 (etc.) are within the buffer, then STATIONARY = TRUE. If, however, days 5,7,9,11 & 13 are within the buffer, but not the alternate (even) days in between, then STATIONARY = FALSE I think some kind of moving window buffer will provide the solution, but I don't know how to implement this. I have been trying to get my head around this problem in both ArcGIS and R, but have had no brain waves so far. This is the closest I've got to a solution, but it doesn't quite fit, I don't think: Identification of consecutive points within a specified buffer Here is some dummy data, which approximates my data structure (though in reality I have twice daily locations (midday and midnight) with some locations missing - but I'll worry about that later)               x<-seq(0,15,length.out=20)     y<-seq(10,-10,length.out=20)     t<-seq(as.POSIXct('2013-07-01'), length.out = 20, by = "days")     data<-data.frame(cbind(x,y,t=as.data.frame.POSIXct(t)))                           x           y          t     1   0.0000000  10.0000000 2013-07-01     2   0.7894737   8.9473684 2013-07-02     3   1.5789474   7.8947368 2013-07-03     4   2.3684211   6.8421053 2013-07-04     5   3.1578947   5.7894737 2013-07-05     6   3.9473684   4.7368421 2013-07-06     7   4.7368421   3.6842105 2013-07-07     ... ...         ...       ...


title: one-to-many join using arcpy

text: I have a feature class that contains polygon features with a unique ID column. This shape-file is to be joined with a csv table that I've exported to a file geodatabase table. The unique identifier in the table is the tuple (IDX, UID), there exists multiple columns with the same UID. I would like to join the polygon feature and the table based on the relation ID=UID using arcpy. That relation however is a one-to-one relation and I would like it to be a one-to-many relation so that one record in the polygon feature class table is joined with multiple records in the table. Any ideas on how to do this using arcpy?


title: Field Calculate Degree Minute Second in Different Format

text: I have a field that is in degree minutes seconds, populated from ET Geowizards COGO Inverse tool. It populates the values as 90-12-28.15 where 90 is the degree, 12 is the minutes and 28.15 is the seconds. I need to use this field to label on a map. The format of the label is calling for the degree minutes seconds characters (i.e. 90° 12' 28.15") so I want to add a new field and use VB/Python to field calculate these characters in and remove the dashes. I've been working on this for a while and I can't seem to figure it out, but I'm sure this is doable. Thanks, Mike


title: Is it possible to intall the opengeostack and QGIS nightly?

text: Did someone managed to get the opengeostack running alongside qgis so the plugin suite-qgis-plugin can be used?   I am trying it on ubuntu1304, installing opengeostack works ok, but then when installing qgis-nightly, i get an error: > dpkg: error processing > /var/cache/apt/archives/libgdal1_1.9.0-3.1ubuntu1_amd64.deb (--unpack): > trying to overwrite '/usr/lib/libgdal.so.1', which is also in package > libgdal 1.9.1-1 > > dpkg: error processing /var/cache/apt/archives/gdal- > bin_1.9.0-3.1ubuntu1_amd64.deb (--unpack): trying to overwrite > '/usr/bin/ogrtindex', which is also in package libgdal 1.9.1-1


title: Debugging FMW memory usage near Group-based Transformers

text: I have a large FME FMW file which is using a lot of memory at runtime. I have been tasked with trying to find ways to reduce the memory footprint. Early on in the workbench there is (what I assume is) a group-based transformer - LineOnAreaOverlayer, which looks for lines that intersect polygons. I suspect that the group-based nature of this transformer is pushing up memory usage, as the thousands of features which are fed into it need to be held in memory and processed all in one go (at least that's the theory). My aim is to halt processing when features are first output from this transformer to see how much memory is being used at this point. I'm using inspection points on the transformer's outputs. However, it seems I either don't understand group-based transformers or inspection points. I break processing on the first feature that is output from the transformer, and the connector feature counts in the rest of the FMW all show features having already gone right through the transformation and being written. I thought that all input features were backed-up in front of the LineOnAreaOverlayer before anything else happened? Can anyone help me clarify this please? Also, any other useful information on tracking memory usage in an FMW would be much appreciated - right now my only options are the log's peak memory usage numbers or trying to match a procmon log with the fme log, which will be a tedious task.


title: Adding rows of data to a non-spatial table in ArcMap with a SDE SQL Server for Metadata purposes

text: I'm attempting to use several non-spatial tables to better organize the large set of feature classes I continually add to a SDE Geodatabase (SQL Server Express). So I have several tables setup and the relationships between them setup. Basically the plan is to use these tables to simply add metadata so to speak so the feature classes can be easily reported on. The tables will show who the owner of the geographic area the data comes from and what type of data it is. None of this needs to spatial. Basically I want organization outside of needing to add this metadata type data in every record of the feature class. However, after creating these non-spatial tables there doesn't seem to be way to add any data to these inside of ArcMap. The Editor doesn't work since it since there are no spatial layers present. I've tried arcpy and I only get errors. Is the only way to use these non-spatial tables is to create some interface outside of ArcMap? I was think I could use some 3rd Python module to add data directly to SQL Server but have no idea where to find that. Also, how would I create a UI for the Python Scripts? Any ideas or direction would be appreciated.


title: Create a graphic

text: I'm working with ArcGis 10 and I have different polygon features in different layers with their attributes like shape_area. I would like to create a unique graphic where is possible to see the respective percentage of each one. Can I do this if I have different layers? if is not possible is there a good way to do this graph? Thank you.


title: Why does my gpx layer go to the wrong spot even if they have the same CRS?

text: Help! I am a QGIS beginner with a big problem! I am trying to load a gpx file from Garmin into my map. The project CRS is WGS 84, and the incoming data is in CRS WGS 84. When I look at the bottom of my map in QGIS or on the coordinate capture, all of my locations appear to be in the UTM format. However when I load in the gpx files, the points go to their long/lat locations, which is of course very far away from everything else because everything else as UTM coordinates in the thousands. I have on the fly enabled... but I am not sure how to fix this. Even if I save the GPX data in coordinate form, it still loads to its lat/long numbers as UTM coordinates. I am hoping there is an easy fix that I am too inexperienced to know about?


title: Trying to use Near to relocate points to the waterways

text: We are trying to relocate elements that are near the water but not quite on the water and relocation them to the nearest point. I'm trying to build the parameters to use the Near tool programmatically with a point on a small island. It should find the point right away but returns a FAIL error. Can anyone see why I still get a FAIL??? This is the latest update in code that we tried. The tool seems to need FeatureSets only so I added logic to create Featuresets from the FeatureClass.                   messages.AddMessage("Creating the SearchNear tool.");         IGPUtilities3 util = new GPUtilitiesClass();         UID ClsID = new UIDClass();         ClsID.Value = "{52353152-891A-11D0-BEC6-00805F7C4268}";         IFeatureWorkspace fws = (IFeatureWorkspace)GPUtilities.GetInMemoryWorkspace();         IObjectClassDescription objectClassDescription = new FeatureClassDescriptionClass();         IFields fields = objectClassDescription.RequiredFields;         int shapeFieldIndex = fields.FindField("SHAPE");         IField field = fields.get_Field(shapeFieldIndex);         IGeometryDef geometryDef = field.GeometryDef;         IGeometryDefEdit geometryDefEdit = (IGeometryDefEdit)geometryDef;         geometryDefEdit.SpatialReference_2 = GetSpatialReference();                    geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPoint;         IFeatureClass pointFeature = fws.CreateFeatureClass("tempPointFeatureClass", fields, ClsID, null, esriFeatureType.esriFTSimple, "SHAPE", string.Empty);         IFeatureBuffer buffer = pointFeature.CreateFeatureBuffer();         IFeatureCursor featureCursor = pointFeature.Insert(true);         IGeometry pointShape = ((IFeature)inputPoint).Shape;         buffer.Shape = pointShape;         featureCursor.InsertFeature(buffer);         featureCursor.Flush();         Marshal.FinalReleaseComObject(featureCursor);         IFeatureLayer pointFeatureLayer = new FeatureLayerClass();         pointFeatureLayer.FeatureClass = pointFeature;         ((ILayer)pointFeatureLayer).Name = pointFeature.AliasName;         pointFeatureLayer.DataSourceType = "Point Feature Class";         geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPolygon;         IFeatureClass polygonFeature = fws.CreateFeatureClass("tempPolyFeatureClass", fields, ClsID, null, esriFeatureType.esriFTSimple, "SHAPE", string.Empty);         buffer = polygonFeature.CreateFeatureBuffer();         featureCursor = polygonFeature.Insert(true);         ((IGeometry)water).SpatialReference = GetSpatialReference();         buffer.Shape = (IGeometry)water;         featureCursor.InsertFeature(buffer);         featureCursor.Flush();         Marshal.FinalReleaseComObject(featureCursor);         IFeatureLayer polyFeaturelayer = new FeatureLayerClass();         polyFeaturelayer.FeatureClass = polygonFeature;         ((ILayer)polyFeaturelayer).Name = polygonFeature.AliasName;         polyFeaturelayer.DataSourceType = "Polygon Feature Class";              Near relocateCompute = new Near()         {             in_features =   (FeatureLayer)pointFeatureLayer,             near_features = (FeatureLayer)polyFeaturelayer,             location = "LOCATION",             angle =    "NO_ANGLE",             method =   "GEODESIC"         };                while (incrementor < 4)           {         try         {             relocateCompute.search_radius = (SearchRadius * incrementor * 1852).ToString() + " Meters";             IFeatureClass newPoint = (IFeatureClass)RunTool((IGPProcess)relocateCompute, false);                 ....      The RunTool method basically encapsulates the whole process of running a tool and the boolean is to determine if a string is returned or a CreateFeaterClassFromString is returned.


title: Can you create labels based on an attribute field in a shapefile with Arcpy?

text: I was curious if anyone knew how to access the ArcGIS label engine with arcpy? I want to be able to label point features on my map in a python loop.


title: Can I ignore very low levels of spatial autocorrelation in my data?

text: I have >1 million postcodes clustered in ~7000 middle super output areas for England and Wales which I'm analysing using MLwiN. I have used a multilevel model with random effects and have clustered the postcode observations at the middle super output area. My concern is that when analysing the middle super output areas using a Queens 1st order contiguity weights matrix (in GeoDa), Moran's I came to 0.00324916 (p=0.013) (see graph below). While there is some spatial autocorrelation, this seems to be a very, very small amount - is it small enough for me to ignore? If I couldn't ignore this then I would have to cluster at the regional level leaving me with only 32 clusters. As I'm using MCMC estimation this would dramatically increase the computational intensity of my model runs (leading to days per run instead of hours), hence I would quite like to avoid this if possible. .![enter image description here](http://i.stack.imgur.com/J4mNX.png)


title: ArcSDE performance and table changes on upgrade from 9.3 to 10.1?

text: I'm using ArcSDE 9.3 and SQL Server 2008 R2 and I want to migrate from ArcSDE 9.3 to ArcSDE 10.1. Is there any performance differences between ArcSDE 10.1 and ArcSDE 9.3 ? Is there any change in tables of ArcSDE 10.1 and ArcSDE 9.3? I'm familiar with Online Help about how to upgrade ArcSDE for SQL Server from 9.3 to 10.1 . I couldn't find resources about performance changes in 9.3 to 10.1.


title: Custom OSM, styling points

text: I have my own OSM file, which is created by some application which I written. You can see some of nodes in my OSM file bottom of the post. I want to fill +/-0.0001 latitude, +/-0.0001 longitude around my point. I've tried these to style points, however, sadly, height and width attiributes seem to not working :/ What may be the other solution(s) to fill +/-0.0001 lat-long around my point? `layer-amenity-points.xml.inc`               <Style name="points">         <Rule>           &maxscale_zoom18;           <Filter>[amenity]='mycustomamenity'</Filter>           <PointSymbolizer height="52" width="52" file="&symbols;/overlay_img.png" placement="interior"/>         </Rule>         <Rule>           &maxscale_zoom17;           <Filter>[amenity]='mycustomamenity'</Filter>           <PointSymbolizer height="26" width="26" file="&symbols;/overlay_img.png" placement="interior"/>         </Rule>     </Style>          <Layer name="amenity-points" status="on" srs="&osm2pgsql_projection;">         <StyleName>points</StyleName>         <Datasource>           <Parameter name="table">           (select amenity           from &prefix;_point           where amenity is not null           ) as points</Parameter>           &datasource-settings;         </Datasource>     </Layer>      `custom.osm`               <node id="10000001" lat="0.476527778159" lon="0.999861111111" visible="true">         <tag k="amenity" v="mycustomamenity"/>         <tag k="myvalue" v="573"/>     </node>     <node id="10000002" lat="0.476805555937" lon="0.999861111111" visible="true">         <tag k="amenity" v="mycustomamenity"/>         <tag k="myvalue" v="576"/>     </node>     <node id="10000003" lat="0.477083333715" lon="0.999861111111" visible="true">         <tag k="amenity" v="mycustomamenity"/>         <tag k="myvalue" v="538"/>     </node>      **Update** With MarkerSymbolizer in the answer below, I couldn't find how to draw squares. There is just circles can be created by MarkersSymbolizer. **I want to draw squares.**


title: How to display image on the Graphics View widget in QGIS

text: The idea is to create a plugin in QGIS in which on mouse clicking of the layer show up a popup with Graphics View widget on it. I found this pure code:               self.scene = QtGui.QGraphicsScene()     self.scene.addPixmap(QtGui.QPixmap('example.jpg'))     self.dlg.ui.grfPicture.setScene(self.scene)      Maybe somebody have some experience with this and can help?


title: take record from ms-access using python and show it in Qt

text: I have a ms-access database and i try to make QGIS plugin that take record from table and show it in Qt, database have format like this :               NO  Skala 1     Skala2      Rekomendasi     ================================================     1.  1000        1000        Skala Sama     2.  1000        5000        Skala Bisa digunakan     3.  5000        10000       Skala Tidak Cocok     4.  5000        50000       Skala Terlalu Besar      what i try to make is, when i input skala1 (ex:1000) and skala2 (ex:5000) it automatic take "Skala Bisa Digunakan" from database and put it in Qt, i try to use argument like :                def Compskala(self):             skal1 = self.ui.lineSkl_1.text() #(input)             skal2 = self.ui.lineSkl_2.text() #(input)             bandskal = ''                    #(from table)             if skal1 and skal2:                 bandskal = ' How to direct this to table?'      but i really didn't know what i must do next to direct my code to take record from table Can someone help me? Sorry for my bad english


title: Update a field attribute in one column based on a change to another column in one table

text: When editing I want to be able to change the attributes of one field and upon doing this allow for the automatic update of another field in the same attribute table. Essentially I have one layer with a code, description and legend information. I would like to be able to edit the code and have the description and the legend information update with this edit with its correct description and legend information. Is there any way to do this within ArcMap's tools or an associated script? Thanks,


title: ogr2ogr reprojection works on command line but shifts output when same command run from inside python?

text: Reprojection with ogr2ogr works fine on command line, but shifts the polygon to northeast when same command is run from inside python using os.system(...) Any ideas why this could be? I hard-coded the full ogr2ogr path to avoid hitting different installations from each environment. PRJ files for both versions are identical so the difference must be in the "shp" file. e.g. Here is the command:               /opt/local/bin/ogr2ogr -s_srs epsg:28992 -t_srs epsg:4326 -a_srs epsg:4326 /mnt/output_4326.shp /mnt/input_28992.shp


title: Create a stacked bar chart with matplotlib

text: I'm trying to create a stacked bar chart in python with matplotlib and I can draw my bar one up the other               # -*- coding: utf-8 -*-     import psycopg2     import matplotlib.pyplot as plt     import numpy as np          # Connect to an existing database     conn = psycopg2.connect("dbname=RABATTEMENT user=postgres password=#######")          # Open a cursor to perform database operations     cur = conn.cursor()          # this is the query we'll be making     query = """             SELECT             trains, ideal, correct, deficient, without             FROM __matplotlib.rabat             """     # execute the query     cur.execute(query)          # retrieve the whole result set     data = cur.fetchall()     cur.close()     conn.close()          # unpack data     trains, ideal, correct, deficient, without = zip(*data)          N = len(trains)         #-- Number of trains     ind = np.arange(N)      # the x locations for the groups     width = 0.35            # the width of the bars: can also be len(x) sequence          p1 = plt.bar(ind, ideal, width, color='green')     p2 = plt.bar(ind, correct, width, color='yellow', bottom=ideal)     p3 = plt.bar(ind, deficient, width, color='red', bottom=correct)     p4 = plt.bar(ind, without, width, color="grey", bottom=deficient)          plt.ylabel('%')     plt.title('Qualite du rabattement par depart de train')     plt.xticks(ind+width/2., trains)     plt.yticks(np.arange(0,100,25))          plt.show()      My problem is that the documentation only show how to stack 2 bars. So, you have p1 (ideal) and p2 (correct) where bottom=ideal. If I want to add a 3rd bar, I thought that putting for p3 bottom=correct (p2), that you draw it correctly, but the 3rd bar is drawn behind the first 2, and it's bottom value is p2. ![Figure](http://i.stack.imgur.com/GrFnw.png) My original data (fields are in French, but I translated them in the right order: ![Database](http://i.stack.imgur.com/Plo72.png) I tried using int() or..., but since my variables are lists. What I'm trying to achieve is a 100% total each time as in : ![Real results](http://i.stack.imgur.com/2gEhS.png) Those are ok, but where made in Excel after a VERY tedious process than I'm trying to automate with a python script...


title: How to delete selected features in Qgis using Python

text: I have a shape file with only one field _DN_. I would like to delete all features which does not have specific value of the field. Let say the _DN_ can be 1,2 or 3. So far I was able to find out how to select Features and how to delete feature, but I did not succeed to combine them. Let say I want only features with DN = 3 to remain: Selection of the Features found here, however, I do not know if it does what suppose to since I did not find the way how to check what is inside.               select = layer.getFeatures(QgsFeatureRequest().setFilterExpression(u'"DN"!=3'))      delete features found here               res = Ilayer.dataProvider().deleteFeatures()      As far as I know, the `deleteFeatures()` requires IDs of the features, but I have no idea how to get them from from `select` which is `QgsFeatureIterator`.


title: infoWindow does not open on initial page load (ArcGIS Javascript API)

text: I am working on a webapp using the ArcGIS Javascript API. I noticed that when I first load the map (currently living on a development server), the infoWindow does not open (when clicking on a feature). When I refresh/reload the map, the infoWindow then works as it should. This behavior happens in both IE and Firefox. I am guessing there is some problem with my initial load, perhaps the 'init' function is loading before the libraries are finished loading somehow..? I appreciate any suggestions. Below is the the HTML setup. Currently I am using $(function () {                   dojo.ready(init);     });      to add the init setup in the main .js file. Thanks, Jason HTML:               <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN" "http://www.w3.org/TR/html4/strict.dtd">     <html>     <head>     <meta charset="utf-8" />     <meta http-equiv="X-UA-Compatible" content="IE=7", "IE=9" />          <link href="mainCopy.css" rel="stylesheet" type="text/css" />     <script src="jquery.js" type="text/javascript"></script>                    <link rel="stylesheet" type="text/css" href="http://localserver/arcgis_js_api/library/2.1/arcgis/js/dojo/dijit/themes/claro/claro.css">     <link rel="stylesheet" type="text/css" href="http://localserver/arcgis_js_api/library/2.1/arcgis/js/dojo/dojox/grid/resources/Grid.css">     <link rel="stylesheet" type="text/css" href="http://localserver/arcgis_js_api/library/2.1/arcgis/js/dojo/dojox/grid/resources/claroGrid.css">     <link rel="stylesheet" type="text/css" href="http://localserver/arcgis_js_api/library/2.1/arcgis/js/dojo/dojox/image/resources/image.css">                    <script type="text/javascript">     var djConfig = { parseOnLoad: true };     </script>     <!-- Tell the browse to include the dojo library -->     <script type="text/javascript" src="http://localserver/arcgis_js_api/library/2.1/arcgis/"></script>     <script src="application.js" type="text/javascript"></script>     </head>          <style type="text/css"> </style>               <body class="claro">      ......


title: How to use default identify tool in ArcMap from C#?

text: I made a custom "Find" tool(winform) which is called after specific button presed in ArcMap toolbar. In this new "Find" form I need to use identify for one or more selected items in DataGridRow. So far I failed to find a way to call default "Identify" tool from ArcMap. I could make my own tool, but that seams like a waste of time, as I dont need any custom functionality in it, besides the default. Is there a way to call ArcMap "Identify" form(with feature properties etc.) from my custom winform tool passing one or more features? I'm relatevly new to ArcMap as well as C#. Any help would be appreciated.


title: esri.dijit.popup custom layout

text: I'm looking for a way to customize the layout of popup. What is asket to me is to change the defined labels of navigation inserting a custom string and move and change the icons to navigate the first change is something like this: change the header title "(1 of 2])" to "element 1 of 2" the second change is to move the navigation tools from the header to the footer and to use labels instead of images (something like previous/next) Is that possible or do i need to implement a own popuptemplate? thanks


title: OpenLayers Vector KML draw black when change the style

text: I am trying to change the style of an OpenLayers.Vector using radio buttons and redraw it, but the result is a black vector. I tried own build OpenLayers11 and 12, beside api openlayers hosted.               <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/html4/loose.dtd">                <title>See kml data</title>     <meta http-equiv="Content-type" content="text/html;charset=UTF-8">     <meta http-equiv="Content-Style-Type" content="text/css">                      I tried beside an Openlayers.Vector, Openlayers.GML. I tried beside simple OpenLayers.Style, OpenLayers.StyleMap and server hosted kml style. The results are shown bellow: ![Initial default style](http://i.stack.imgur.com/wIIAx.jpg) ![User choosen style](http://i.stack.imgur.com/RtwBl.png) Any clue would help! Mihai Edit of the choosestyle(userstyles) function which on the alert returns a style object (checked with Firebug DOM):               function choosestyle(userstyles) {     //alert(userstyles.elements[0].value);     //var userstyle = document.forms[0].elements[iStyle];     //alert(userstyle.length);     for(var i = 0; i < userstyles.length; i++)        {        if (userstyles[i].checked)           {                 kml.style = userstyles[i].value;             kml.redraw();                 alert(kml.style);             userstyles[i].checked = false;           }        }     }      But the map has the same black color.


title: Can I open ArcView 3 project files (*.APR) with Qgis? Is there a converter?

text: I would like to open files from an old ArcView 3.x project with Quantum GIS to work with the data. But this project data is saved as an APR-file, which I can not open with Quantum GIS. Is there any possibility to convert the data (for example to QGS files) to work with it? Thank You.


title: Exporting vector to kml with ogr2ogr without extended data

text: If I convert a shapefile as kml in QGIS, GRASS or ogr2ogr directly, defining the name and/or description field with the -dsco tag, all columns in the attribute table of the input layer are exported to the kml as ExtendedData. Is there a way to create the kml with the Name and Description only, i.e., without the Extended data?


title: Opening Old MapInfo Workspaces

text: I have many "old" (10+ year) MapInfo workspaces that I want to reopen but cannot locate some of the tables. MapInfo will not open the workspaces without these tables and the data/objects I want only exist in the editable layer of these old workspaces. Is there a way of overiding MapInfo's inability to open the workspace without these tables??


title: FME: how to extract the coordinate of the North West corner

text: In an FME translation from a polygon geodatabase feature class to a DGN I am trying to place a label in the top left corner of a box. I have tried using CoordinateExtractor to find the coordinates of a corner (then 2DPointReplacer and TexAdder to add text to this coordinate) but this works on the number of the vertex, which unfortunately isn't constant as not all of my boxes are quadrilateral. What I need to find is a way to find the upper left vertex. Does anyone know how to do this?


title: Installing Spatialite-GIS

text: I was wondering if there is a way to install Spatialite-gis on a Windows & machine? I have found may posts and links about this, but I have yet to find an installer for spatialite-gis. This Google Group post links to different sites, but there are not installers available, am I missing something? https://groups.google.com/forum/#!topic/spatialite-users/5qHN1oUq5oQ This is the direct link from GAIA to spatialite-gis, but there are no installers. https://www.gaia-gis.it/fossil/spatialite_gis/index This is a link to other downloads, including spatialite-gui, but nothing for spatialite-gis http://www.gaia-gis.it/gaia-sins/windows-bin-amd64/ Any help or clarification would be useful. Thanks...


title: Using attributes of a polygon in one layer to identify polygons in another

text: So I'm in the process of developing a web application, I am using open source tools so qgis, postgis, play framework for the web app. And I have come across an interesting problem that I'm trying to work on if someone might be able to point me in the right direction. What I need to do is this I have one polygon layer A with attributes about the polygons, and polygon layer B with attributes also. So if a polygon in layer B falls within a polygon in layer A look up the attributes for the polygon in A, now compare those with the attributes of the polygon in layer B to determine if the attribute for the polygon in B is a multiple ( 2x, 3x, 4x, etc ) of the attribute of the polygon in A. I am going to be using the results of this in the web application for other things. So I have multiple options, but looking for some direction as to where to start. Do I need to use something like sextante? There are several of these operations that I need to do. It Should be repeatable, but only need to be repeated when updates to the layers come out 1 or twice a year! But it will be acting on several states data layers. Thank you in advance for any help that you all may bring. Thank you PolyGeo for suggesting to edit the post


title: Creating contours from partial data set

text: Is it possible to create contour lines based on part of the points and not the entire dataset?


title: Print City Map in High Resolution with TileMill (OSM Data)

text: I'm rather beginner in Print Map and GIS. I'd like to print the map of my little city using the data in OpenStreetMap. I've got an environment(local) with OSM DataBase and I'm using TileMill to create StyleSheet and export data. My problem is the definition of export. As you know, in TileMill you can export data in SVG, PDF and XML, and you can define just the BBox. The Resolution seems hard coded and set to 72dpi. To increase the quality and the details of export image I have to reduce the BBox (since to have a higher zoom level). Since I need to print a classic city map (the paper size is, more or less, 80cm X 50cm and the city area 2km X 1km), does anybody have an idea for print it with a good resolution? If i set 300 dpi as resolution, I'll solve my issue, but maybe I can't. As I said before, I have a local environment, but if someone knows some other ways ... I appreciate your help Any ideas? Thanks in advance and Regards, Paolo


title: Display DEM in WPF application

text: In my application which is based on .net4/C#/wpf, I need to display a decent sized terrain file (eg. GeoTiff). It is an elevation file which should be viewed in 3D. The user should be able to draw small polygon it and interact with the map (zoom, pan). I looked for these options:   * Global Mapper SDK: No 3D for free license. Win Forms SDK so some work to bring things in WPF.    * ArcGIS runtime SDK: I don't think 3D is free (or at all) here. Plus, runtime SDK needed.    * Helix3D toolBox: Gets slow with large data. Supports only (*.bt) files. I want to know some other option primarily free or reasonably expensive. Kindly feel free to ask further details.


title: converting UTM cordinates into map area

text: I have coordinates in UTM coordinates for 3 points. I do not have licenced ArcGIS software. Is there any available app or place on internet where , if the above 3 coordinates can be put and the result is a map of the area enclosed by the above points along with the area displayed (preferably in a format which can be sent directly to printer).


title: How can one read the response of a GetFeature request with Openlayers?

text: I am using Openlayers to send a GetFeature request to Geoserver. With the help of Firebug I can see that the response contains some of the attribute data I need. I am interested in showing some attributes once the user clicked on a feature. How can I make sure that geoserver sends all attributes and how can I read them in Openlayers? Here is my code, for now I only succeeded at displaying the feature selected on my map using an Openlayers example:               Feature = new OpenLayers.Layer.WMS(             "Ramsar sites",             "http://ramsardev.grid.unep.ch:8080/geoserver/wms",             {                 layers:'ramsar_sdi:features',                 transparent:'true',                  visibility: true,                 isBaseLayer: false             });                  map.addLayer(Feature);          select = new OpenLayers.Layer.Vector("Selection", {styleMap: new OpenLayers.Style(OpenLayers.Feature.Vector.style["select"])});          map.addLayer(select);          control = new OpenLayers.Control.GetFeature({             protocol: OpenLayers.Protocol.WFS.fromWMSLayer(Feature),             box: false,             hover: false,             multipleKey: "shiftKey",             toggleKey: "ctrlKey"         });          control.events.register("featureselected", this, function(e) {             select.addFeatures([e.feature]);              });          control.events.register("featureunselected", this, function(e) {             select.removeFeatures([e.feature]);         });           map.addControl(control);             control.activate();           The problem is that I can't see any of my attributes anywhere in the event object `e` when the event `featureselected` is triggered. Any help will be much appreciated!


title: qgis2 html annotation problem - follow-up

text: Referring to this question, I can't find out how to 'attach' an html annotation to an existing object, to take data from that object and use it in the placeholders in that annotation object. A test case video is here, also referenced in the comments on the first answer to that question. For me, that comment doesn't show up on the first page - you have to hit 'show 1 more comment' so maybe it never got seen? Any ideas? Thanks


title: How do I active "Open GRASS tools"-button?

text: I have opened to GRASS layers and need to use the "Open GRASS tools"-button. This button remains unactivated somehow. Whilst the "Edit GRASS vector layer" right next to it works perfectly fine. What can I do?


title: QGIS problems getting layers to line up correctly.

text: I am trying to add a raster layer and a vector layer to qgis. I have set both of their CSR to NAD83 Oklahoma North State Plane. It doesn't seem to work even with all the files set to those CSR. What is likely the issue is that the raster has no extent, and QGIS doesn't seem to have a way to set the extent on rasters. Coming from ArcGIS world where all these are set on the fly... Where else might I be going wrong.


title: how to extract a variable at specific time and depth from netcdf in r

text: I have the following netcdf file containing 7 variables and 4 dimensions.               [1] "file data.nc has 4 dimensions:"     [1] "depth   Size: 12"     [1] "time   Size: 37"     [1] "x   Size: 305"     [1] "y   Size: 432"     [1] "------------------------"     [1] "file data.nc has 7 variables:"     [1] "int stereographic[]  Longname:stereographic Missval:NA"     [1] "float longitude[x,y]  Longname:longitude Missval:1e+30"     [1] "float latitude[x,y]  Longname:latitude Missval:1e+30"     [1] "short chla[x,y,depth,time]  Longname:chla Missval:-32767"     [1] "short nitrat[x,y,depth,time]  Longname:nitrat Missval:-32767"     [1] "short phosphat[x,y,depth,time]  Longname:phosphat Missval:-32767"     [1] "short pbiomass[x,y,depth,time]  Longname:pbiomass Missval:-32767"      I would like to extract several months at a specific depth and converted to a ASCII raster. I am able to use this command myraster<-raster("data.nc", varname="nitrat", lvar=1, band=2, level=20) but I am confused about which depth and time I am actually extracting, it is unclear to me how to use the parameters of the raster function to point to the specific time and depth. Also I am not able to read a variable using the function get_att and I get this message               test<-open.ncdf("data.nc")     Error in ncatt_get(test, "time") :        Error, first passed argument must be an object of class ncdf4      I would like to be able to get a list of the different time and depth so I can use the correct parameters in the raster function.I have both "ncdf","ncdf4" packages installed.I am working in Windows with R 3.0


title: clustering points in postgresql to create contour map

text: UPDATE: To try to make this a little more simple, is it possible to create a query which selects a group of points such that each point in the group is at most X units away from at least one other point in the group? If that is possible, then can it be extended to return more than one row such that each row is a group of points as specified above? I'm starting to think this isn't possible, so I'm working on some java code to do it, but I still have some hope. Original post and info below:      I have a postgresql / postgis table with a grid of weighted points, and I'd like to essentially create a contour map of those points. My goal is to get one or more polygon that surrounds points of a certain weight, and then repeat for other weights. I'm going to convert those polygons to KML for display, which I can figure out later and is only relevant because I don't necessarily need this to be something dynamic or efficient - I can do it manually and save off the polygons at each step. I've read a lot of other posts about this but they all seem to require additional libraries or software, and my points are already in a grid so the ST_SnapToGrid solution doesn't make sense to me. Is there a strictly postgresql solution? My thought is that there would be a query that creates a polygon (using ST_ConvexHull or ST_ConcaveHull) around points that are within some distance from each other (which should be a set value because my points are already equally spaced on a grid). Example lat/lons/weights (for distance reference):               41.96892, -87.9351, 1904     41.97392, -87.9351, 1940     41.97892, -87.9351, 1978      Table:               CREATE TABLE temp_weighted_points     (       id serial NOT NULL,       latitude double precision,       longitude double precision,       weight double precision,       type character varying(10),       geom geometry(Point,4269),       CONSTRAINT temp_weighted_points_pkey PRIMARY KEY (id)     )      I've tried this link (among others) but can't get any of those solutions to work: http://stackoverflow.com/questions/13087288/separate-gis-points-into- groups-of-points-that-are-close-to-each-other Any suggestions? EDIT:   Regarding what the data looks like, it's a grid of drive times over the city of Chicago. To create the data, I found a boundary shape file of the city of Chicago, then used a postgresql query to create a grid of equally spaced points within that boundary / polygon (for pictures and the query, see How to create regular point grid inside a polygon in Postgis?). Each point is separate by .005 lat / lon units (degrees? not sure what it's called) as can be seen in the example lat/lons above. Once I had these points, I used to google's drivetime api to generate drive times from each of the points to four different locations, and then took the smallest drive time between the four locations and stored that data as the weight for the grid point (the weight is in seconds, so the first example point above has a drivetime of about 31 minutes). That's what I have at this time - points that are equally spaced over an area / polygon and that each have weights / heights associated with them. What I'm trying to get is a contour map of this data, which will basically look like the city of Chicago with 4 depressions (the locations) or an upside down mountain with four peaks, because points on the grid close to the locations will have very small / short drive times, while points farther away will be higher.


title: Do ArcGIS Online feature services support a DISTINCT or returnOnlyUniqueValues option?

text: ArcGIS Server REST API has:               returnDistinctValues         Description: If true, then returns distinct values based on the fields specified in outFields. This parameter applies only if the supportsAdvancedQueries property of the layer is true.     Values: <true | false>     Example: returnDistinctValues=true      http://resources.arcgis.com/en/help/arcgis-rest- api/index.html#//02r3000000r1000000 Is there an equivalent when querying agol feature service?


title: csv to shp python and formating a row

text: I´m doing a script to convert csv to shp, but I want to format one of the columns. It´s a string then I´d like to format like this a[8:10]+' / '+a[11:16]+' UTC' being a the value in the column. this is the code working               #Set up blank lists for data     nombre,codigo,y,x,time=[],[],[],[],[]          #read data from csv file and store in lists     with open('C:/archivecsvToRead.csv', 'r') as csvfile:         r = csv.reader(csvfile, delimiter=',')         for i,row in enumerate(r):             if i > 0: #skip header                nombre.append(row[0])                codigo.append(row[1])                y.append(float(row[2]))                x.append(float(row[3]))                time.append(row[4])               #Set up shapefile writer and create empty fields     w = shp.Writer(shp.POINT)     w.autoBalance = 1 #ensures gemoetry and attributes match     w.field('Longitud','F',10,8)     w.field('Latitud','F',10,8)     w.field('Date','D')     w.field('Date / Time','C',50)     w.field('ID','N')          #loop through the data and write the shapefile     for j,k in enumerate(x):         w.point(k,y[j]) #write the geometry         w.record(k,y[j],codigo[j], time[j], nombre[j]) #write the attributes          #Save shapefile     w.save(out_file)      How can I format the column time in order to show a[8:10]+' / '+a[11:16]+' UTC' being 'a' the value in the column ?


title: Script tool problems

text: I'm having trouble making a script tool in ArcGIS. I have developed a script that will edit certain values in the attribute table. Without incorporating a script tool this script works fine and edits the relevant attributes adequately. Here it is:               import arcpy          fc = "C:\Users\maureen\Documents\ArcGIS\EDRN\EDRN_LINK.shp"     fields = ["FID", "RIVERNAME", "SOURCE", "DRAINS", "RIVERTYPE", "LCC_USER"]          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == 4402:                 row[1] = "Calder"                 cursor.updateRow(row)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == 4402:                 row[2] = "2009 Aerial Photography"                 cursor.updateRow(row)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == 4402:                 row[3] = "NO"                 cursor.updateRow(row)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == 4402:                 row[4] = "1"                 cursor.updateRow(row)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == 4402:                 row[5] = "T0015685"                 cursor.updateRow(row)      When I edit the above and try to create a script tool my script is as follows:               import arcpy          fc = "C:\Users\maureen\Documents\ArcGIS\EDRN\EDRN_LINK.shp"     fields = ["FID", "SOURCE", "RIVERNAME", "RIVERTYPE", "DRAINS", "FRMMAINRIV",   "LCC_USER"]          f_id = arcpy.GetParameterAsText(0)     srce = arcpy.GetParameterAsText(1)     rvrnme = arcpy.GetParameterAsText(2)     rvrtyp = arcpy.GetParameterAsText(3)     drns = arcpy.GetParameterAsText(4)     fr_main_riv = arcpy.GetParameterAsText(5)     lcc_us = arcpy.GetParameterAsText(6)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == "f_id":                 row[1] = "srce"                 cursor.updateRow(row)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == "f_id":                 row[2] = "rvrnme"                 cursor.updateRow(row)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == "f_id":                 row[3] = "rvrtyp"                 cursor.updateRow(row)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == "f_id":                 row[4] = "drns"                 cursor.updateRow(row)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == "f_id":                 row[5] = "fr_main_riv"                 cursor.updateRow(row)          with arcpy.da.UpdateCursor(fc, fields) as cursor:         for row in cursor:             if  row[0] == "f_id":                 row[6] = "lcc_us"                 cursor.updateRow(row)      Once I have made the necessary edits and ran the script tool wizard the script doesn't work. The detail box states that the scrip has completed and that it was successful, but non of the specified edits have been occurred and only a blank space remains. The above makes me think that the problem lies in the script tool wizard, specifically the data types and parameter properties I am using. Can anyone suggest any changes I can make to the the script tool wizard or the amended script?


title: Best method to find routes from a set of points

text: I have lines table (postgis enabled database) representing the streets. I have collected few gps points along a route and I would like to match it with the streets in order. **Scenario** \- I am jogging along street A followed by street B, street C. Let's say I have 3 gps points collected when I am in running along each of those streets - total 9 points.      My desired output is a list - **[street A,street B,street C]**      **What I have done so far** > Get street name given a gps point. This is based on closest distance from > the point to any street in the vicinity. As expected when there are > overlapping streets, the results gets skewed.    **What I think can work** > Take pair of points and find the street(line) making smallest angle with > line connecting the pair of points. If this is a good approach, what should > I expect when I make transition from one street to another (like points > (3,4) or (6,7) in the scenario)? They might make a smaller angle with > another street close by. Is there a better way to achieve this or can you make my approach more robust?


title: "ERROR: invalid input syntax for integer" importing data from Spatialite into PostGIS through DB Manager

text: I am trying to import data from Spatialite into PostGIS through drag and drop in DB Manager but all I received was an error: `Feature write errors:   Creation error for features from #1 to #202. Provider errors was:   PostGIS error while adding features: ERROR: invalid input syntax for integer: ""` ![enter image description here](http://i.stack.imgur.com/m2U6t.png) Is there any way to fix this?   Thanks.


title: Correctly drawing wrapping circles using Basemap

text: I've been using the functions from here to draw circles on maps using Matplotlib Basemap. This works well, but breaks badly at the map edges: ![Wrong](http://i.stack.imgur.com/67Dgb.png) I'm wondering whether there's a general approach for detecting this problem and correcting the errors in the points? I've tried separating any `lon, lat` points with a `lon < 0`, constructing two separate polygons using Shapely, then plotting them separately using Descartes, but that's not a satisfactory solution, really.


title: Correlation between a raster and a point dataset

text: I have an interpolated surface (in raster form) with temperature values and I also have the measurements of temperature in certain locations that belong to the same region (with the raster). I want to plot in a scatterplot the relationship between the raster (estimated values) and the measurements (observed values) and compute the correlation between the two datasets and also the line that best describes the relationship among the two datasets. Any idea how can I do this in ArcMap or QGIS?


title: Any way to perform zonal statistics on column other than value or do reclass faster?

text: I have to perform ArcGIS zonal statistics on a raster, but the column I am interested in is not the 'value' column. Now, I can reclass the raster so that the column of interest becomes the 'value' column. However, the reclassification is super slow (reclassifying a 60 m raster for the US Midwest, length of reclass ASCII file ~500,000 rows.) Is there any way I can skip the reclassification and just perform zonal statistics on my column of interest? Alternatively, any tips on making the reclassification process much faster will help too. Currently the reclassification is on course to take 10 hours to finish.


title: Extracting all points within a certain distance and height difference from a Raster

text: I have a raster containing tens of millions of points all with X, Y and Z data defined. I need to figure out a simple way to extract and create a new raster or vector with ALL points at a 720 m distance from each other (X,Y) and at 120 m elevation difference (Z) from each other. I have ZERO knowledge of SQL and Python. I have been trying to do this on VBA and came up with a couple of algorithms but the processing time is unreasonable and unrealistic. I am sure there must be a simple GIS approach to accomplish this but cannot seem to find it. I am using ArcMap.


title: How to run GeoAdmin3

text: I coloned Geoadmin repo from github. I tried tu run this localy but no sucess. I have my own rc_ file but still it's not working. Anyone can help me?


title: Highest DEM pixel in polygon

text: I'm working with QGIS and have a big DEM raster layer with height information and a vector layer with many polygons, representing 1 x 1 degree tiles. Now I want to find for every polygon the highest value of the DEM pixels in it.


title: ArcMap add-in (vb.net) datagridview will crash map when right click on the title part

text: I find the problem accidentally. In my ArcMap add-in, there is a search function which allows the user to search the geodb based on the value of a field. It works fine but today when I right click the title in the datagridview casually, the ArcMap crashed. After that I simulate the scenario several times, every time ArcMap crashes. My code is listed bellow:               Sql = "SELECT ITEM_NO AS 'ITEM NO', ITEM_CATEGORY_NO AS 'ITEM CATEGORY NO', ARRIVAL_DATE AS 'ARRIVAL DATE' FROM testDB WHERE ITEM_NO = '" + item_No + "'"          Try             adapter = New SqlDataAdapter(Sql, connObj)             adapter.Fill(ds)             connObj.Close()             DataGridView1.DataSource = ds.Tables(0)             Catch ex As Exception                MsgBox(ex.ToString)     End Try      The result is: ![enter image description here](http://i.stack.imgur.com/PWLQE.jpg) Once I right click the part in the red frame, the map will crash. Can anyone please tell me what's wrong with my code?


title: How to download older version of arcgis?

text: I need to re Install ArcGis after formatting one computer, but on version 10.1 (not 10.2.2 that is available on line). were can i found this version for download ?


title: Popup from multiple services

text:
I am trying to have one infowindow/popup show from multiple services. I have successfully added the second service to my code, however, only 1 result will display. Ex. if I am at full scale, only the zoning layer is returned in the popup. If I zoom in, only the parcel layer is returned in the popup. How can I get both services to show within the popup? Additionally, how can I make it so only the visible layers show within the popup? Here is the code I have edited:               dojo.require("esri.dijit.Popup"); //Infowindow          var identifyTask,identifyParams;                function mapReady(map){             dojo.connect(map,"onClick",executeIdentifyTask);             //create identify tasks and setup parameters         identifyTaskZoning = new esri.tasks.IdentifyTask("https://www.marioncountyfl.org/arcgis/rest/services/Dynamic/Zoning/MapServer");        identifyTaskParcels = new esri.tasks.IdentifyTask("https://www.marioncountyfl.org/arcgis/rest/services/Dynamic/Parcels/MapServer");                  //Set Zoning Parameters        identifyParamsZoning = new esri.tasks.IdentifyParameters();        identifyParamsZoning.tolerance = 7;        identifyParamsZoning.returnGeometry = true;        identifyParamsZoning.layerIds = [0,1,2,3,4];        identifyParamsZoning.layerOption = esri.tasks.IdentifyParameters.LAYER_OPTION_VISIBLE;        identifyParamsZoning.width  = map.width;        identifyParamsZoning.height = map.height;             //Set Parcel Parameters        identifyParamsParcels = new esri.tasks.IdentifyParameters();        identifyParamsParcels.tolerance = 15;        identifyParamsParcels.returnGeometry = true;        identifyParamsParcels.layerIds = [0,1,2,3,4];        identifyParamsParcels.layerOption = esri.tasks.IdentifyParameters.LAYER_OPTION_VISIBLE;        identifyParamsParcels.width  = map.width;        identifyParamsParcels.height = map.height;             //resize the map when the browser resizes        dojo.connect(dijit.byId('map'), 'resize', map,map.resize);       };            function executeIdentifyTask(evt){      identifyParamsParcels.geometry = evt.mapPoint;      identifyParamsParcels.mapExtent = map.extent;      var deferred = identifyTaskParcels.execute(identifyParamsParcels);      deferred.addCallback(function(response){        if (response.length > 0) {     //         console.log(response.length)     //         response is an array of identify result objects         //         Let's return an array of features.         return dojo.map(response, function(result){           var feature = result.feature;           feature.attributes.layerName = result.layerName;           if(result.layerName === 'TaxParcel'){               console.log(feature.attributes.NAME);               var template = new esri.dijit.PopupTemplate({                             title:"Parcels",                             description:"Parcel ID: {Parcel Identification Number}
 Address: {Site Address}"                             });               feature.setInfoTemplate(template);             }                      return feature;         });         }       else {           identifyParamsZoning.geometry = evt.mapPoint;           identifyParamsZoning.mapExtent = map.extent;           var deferred = identifyTaskZoning.execute(identifyParamsZoning);           deferred.addCallback(function(response){             // response is an array of identify result objects                 // Let's return an array of features.             return dojo.map(response, function(result){               var feature = result.feature;               feature.attributes.layerName = result.layerName;             if(result.layerName === 'Zoning Classifications'){               console.log(feature.attributes.NAME);               var template = new esri.dijit.PopupTemplate({                             title:"Zoning",                             description:"Zoning: {Zoning Classification}
 Description: {Zoning Description}"                             });               feature.setInfoTemplate(template);             }             else if (result.layerName === 'Assisted Living Facilities'){               var template = new esri.dijit.PopupTemplate({                             title:"Assisted Living Facilities",                             description:"Parcel ID: {Parcel}
 Applicant: {BusinessApplicantName}
 No. of Residents {NbrResidents}"                             });               feature.setInfoTemplate(template);             }                         else if (result.layerName === 'Family Divisions'){               var template = new esri.dijit.PopupTemplate({                             title:"Family Divisions",                             description:"Parcel ID: {Parcel}
 Project Number: {project_nbr}
 Status: {Status}"                             });               feature.setInfoTemplate(template);             }                         else if (result.layerName === 'Policy 120'){               var template = new esri.dijit.PopupTemplate({                             title:"Policy 120",                             description:"Parcel ID: {Parcel}
 Business Name: {BusinessName}
 Date of Letter: {DateOfLetter}"                             });               feature.setInfoTemplate(template);             }                         else if (result.layerName === 'Zoning Changes'){               var template = new esri.dijit.PopupTemplate({                             title:"Zoning Changes",                             description:"Parcel ID: {Parcel}
 Request Number: {REQUEST_NBR}
 Property Owner: {PROPERTY_OWNER}"                             });               feature.setInfoTemplate(template);             }             return feature;             });           });       };         map.infoWindow.setFeatures([deferred]);         map.infoWindow.show(evt.mapPoint);                 });         map.infoWindow.setFeatures([deferred]);         map.infoWindow.show(evt.mapPoint);             };


title: PostGresql and PostGIS courses in Canada

text: I have been looking for good professional courses for learning more than just the basics for PostGresql nd PostGIS. I have found a number of good courses, but they seem to all be in the US. I was wondering if anyone knows of similar type courses being offered in Canada.               http://www.otg-nc.com/implementing-databases-postgresql/     http://www.otg-nc.com/postgresql-administration/      Thanks...


title: Can I use arcpy.Exists for domains?

text: A script I made creates domains within a selected geodatabase. I need it to not create the domains if they already exist because if they do, it gives an error saying invalid domain name. Is there anyway around this? I have overwrite on, but I assume that doesn't have anything to do with domains? I'm using 10.0 Thanks


title: How to apply schema changes on Esri geodatabase two-way replica set?

text: I have an Esri geodatabase two-way replica set. I need to update the schema on both databases in the replica set. Esri describes doing this using ArcGIS schema tools (compare/export/import replica schema - see here). I have control of both databases in the replica set and I would prefer to make a python script, which updates the schema and run that on both databases in the replica set. Are there potential problems with applying the schema changes to both databases instead using Esri's schema tools?


title: how to build and install pgrouting-2.0 on postgresql-9.3 on ubutntu 12.04 Lts precise?

text: I try to build following instructions from this link https://github.com/pgRouting/pgrouting after downloading from the github. I have done the cmake command and there is my error output:               ~/Téléchargements/pgrouting-master/build$ cmake --DWITH_DD=ON ..     -- PGROUTING_GIT_TAG: pgrouting-2.0.0     -- PGROUTING_GIT_BUILD: 0     -- PGROUTING_GIT_HASH: d6ed2cb     -- PGROUTING_GIT_BRANCH: master     -- PGROUTING_VERSION_REVISION_NAME: d6ed2cb master     -- PGROUTING_VERSION_REVISION_HASH: d6ed2cb     -- UNIX=1     -- WIN32=     -- POSTGRESQL_EXECUTABLE is /usr/lib/postgresql/9.3/bin/postgres     -- POSTGRESQL_PG_CONFIG is /usr/bin/pg_config     -- PostgreSQL not found.     CMake Error at CMakeLists.txt:122 (message):        Please check your PostgreSQL installation.     -- Configuring incomplete, errors occurred!      hOW MUST i DO PLEASE?


title: Setting parameters in Arcmap for dynamic python script

text: I have a simple python script that works when run from PythonWin, but when I try to make it dynamic to run inside ArcMap I keep receiving: : ERROR 999999: Error executing function. Failed to execute (CreateFeatureclass). It must be something I am doing when setting the parameters from the script properties from within ArcMap. I've created scripts with user input before, so for the life of me I can't figure out why it won't work.               import arcpy     import os          arcpy.env.overwriteOutput = True          path = arcpy.GetParameterAsText(0)     name = arcpy.GetParameterAsText(1)     type = "POINT"          m = "DISABLED"     z = "DISABLED"     prjFile = "C:/Program Files/ArcGIS/Desktop10.0/Coordinate Systems/Projected Coordinate Systems" + \               "/State Plane/NAD 1983 (US Feet)/NAD 1983 StatePlane North Carolina FIPS 3200 (US Feet).prj"     reference = arcpy.SpatialReference(prjFile)          arcpy.CreateFeatureclass_management(path, name, type, "", m, z, reference)      When creating the parameters I have the first input as follows: Data Type = workspace or feature dataset type = required direction = Input Environment = workspace The second input is: Data Type = Featureclass....I've also tried Shapefile as well Type = Required Direction = Output


title: Geocodes fails when published to ArcGIS Server

text: I have a composite address locator that was built using our streets line file and our address points that works great in ArcMap 10.1 sp1. But when published to ArcGIS Server (10.1 sp1) no results are ever returned, whether in a REST environment or within a map service. I've gone through the esri help files and it doesn't appear that I am doing anything wrong. (This is the first time I've ever published one out.) Has anyone else had this problem, and if so, how did you fix it? Is it due to it being a composite? As a side note, I did publish out the individual address locators and those too didn't work. Thank you in advance


title: shapefiles in QGIS and ArcGis

text: Hi I'm working with a set of shapefiles and have some problems with the "align" of them in QGis when I load Openlayeres to reference some points. OLs has an straight horizontal baseline but the shapefiles have a slighted inclination. I tried changing CRS-datum etc with no luck. However, when I open the same shapefiles in ArcGis they have a straight base line similar to that of OL in QGIS. Any ideas of the problem? If there is any problem. https://picasaweb.google.com/117688487090958721662/QGIS


title: Offset points vertical to a linear feature class

text: In the last few months I am experimenting with ArcGIS in order to create a synthesis of tools in model builder for my postgraduate thesis. The problem I am facing for a long time now is the following: I have one (or more it doesn't matter) linear objects (that represent faults) and with the help of "create routes" and "make route event layer" tools I place points on it (representing epicenters of an earthquake). I want to achieve a relocation (or offsetting) of all these points vertically to the trend-line of the linear feature I create with the "smooth line" tool. The only thing I have managed till now is a relocation by changing a small script I found in this site. > `def Shift_XY_Coordinates(shape, x, y): value = distance point = > shape.getPart(0) point.X = x + value point.Y = y + value return point` with the expression               Shift_XY_Coordinates(!shape!, !POINT_X!, !POINT_Y!)      I want to offset the points by 10 km. How am I supposed to add the correct distance every time when the angle of the linear feature class may be different every time? At the example I am using the trend line has an angle of 45 degrees so I tried with the Pythagorean Theorem and found out that I need to add approximately 7000 meters. I don't know If I was clear enough in what my problem is but I'd greatly appreciated if anyone could suggest me a solution or at least guide me at this. I also attach a picture of what I want to achieve. The starting location of the epicenters is the red line.![enter image description here](http://i.stack.imgur.com/FEi96.jpg) Thank you.


title: Use a variable to calculate a field with arcpy as a script in a toolbox?

text: stringvariable = "banana"     arcpy.CalculateField_management("c:\point.shp", "SUBDIRECT", "stringvariable", "PYTHON")      I get an error saying stringvariable is not defined what's the deal? I check stringvariable by printing it to screen and it contains what it should. I am trying to use a string variable to calculate a field with arcpy. I am stumped. Also it works like this when run as a script through idle or outside of arcmap, but now I am trying to add this script to a toolbox and I am getting this error. thanks


title: What Lidar processing tools are available in Python?

text: I have been using FUSION and the command line FUSION Lidar Toolkit (LTK) to process Lidar data. A broad Google search ("Lidar Python") yielded libLAS and pyLAS as Python lidar libraries, however, these appear to provide only read and write access to LAS data. I am particularly interested in creating intensity and density images in addition to canopy surface models from point clouds. Is there a generally accepted set of tools in Python that can accomplish the same sort of tasks FUSION LTK is capable of?


title: How long does geoserver store the cookie it generates from authentication via j_spring_security_check?

text: I have a website which connects to a geoserver, it fetches sessions from geoserver via j_spring_security_check. How long is the cookie given by geoserver valid and can I edit its expiration?


title: Connecting gvSIG with SQLSERVER

text: I am working for the first time with gvSIG and SQL SERVER, I want to connect gvsig with spatial database SQL SERVER. I don't know how to, Is it by using a protocol and which one ? plz help me


title: Error: basic_string::_S_create, (trying to use .ecw files)

text: I have a problem, after resolving my problem with: .ecw is not a supported data file. I am in a new one. When I try to load an .ecw file in QGIS it says: basic_string::_S_create


title: OpenLayers 3 Layer Visibility

text: I'm trying to upgrade my Openlayers 2.12 map to Openlayers 3 to take advantage of the fantastic transition effects on transparent layers (Something OL2 and Leaflet can't do attractively). This and I don't want to lag behind upgrading my sites when OL3 is officially released. In my current site (OL2.12) I use check boxes in a simple HTML menu to toggle layer visibility. I push each layer to an array (I hope I'm right in thinking OL3 now automatically generates an array for the layers called 'layers') and each checkbox calls this function (checkboxes are given a value that is representative of their layers array number):               function layerswitch(evt){         layers[evt.value].setVisibility(evt.checked);     }      In OL3 this no longer works, and I can't find any examples or documentation that details how to set layer visibility, any help would be greatly appreciated. Thank you for your help.


title: How to remove Raster from Raster Dataset programmatically?

text: I have Raster (it is 6 adf files). This raster was added to Raster Dataset in sde database.                           Geoprocessor geoProcessor = new Geoprocessor();                 Mosaic mosaic = new Mosaic();                      mosaic.inputs = str; // it is path to the directory with 6 adf files                 mosaic.target = Settings.Instance.RasterDatasetFullName;                      mosaic.mosaic_type = "LAST";                 mosaic.colormap = "FIRST";                      object outRaster = geoProcessor.Execute(mosaic, null);       How to remove this Raster from RasterDataset? Is there any classes for this purpose?


title: Methods to calculate Viewshed analysis

text: Does anyone know an algorithm to calculate Viewshed analysis? I calculated them in ArcGIS, MapInfo and SAGA GIS. I wish to calculate Viewshed analysis using algorithm or methods which I can use in MS SQL server or any offline service?


title: How to check if feature class or table contains field with specific name before adding it?

text: Im trying to write a script that finds out if the shapefile has a specific field and if not to add the field. I get it to show me the field, but the loop goes through each field, so it would create a new field each time the specified field is not there. But if it's theres I just want it to say it is there and if it's not to create new field. My code so far just says:               "not filled in yet"     "not filled in yet"     "not filled in yet"     "cities.shp"     "not filled in yet"      (prints out like this cause there are five fields) My code so far:               fieldList = arcpy.ListFields("cities.shp")          for field in fieldList:          if field.name == "CAPITAL":         print "cities.shp"          else:          "then add new field"


title: Use ArcGIS-server restapi to return geometries as Json objects

text: I'm very new to ArcGIS-server REST api, but have been using ArcGIS desktop and other web maps for a while. In one of the Buffering example, I found one GP service which takes geometres in as parameters and return buffered gemetries as json objects. This seems to be very interesting to me. I'm thinking about writing another REST GP service myself, which can take a string parameter and create a bunch of gemetries as Json object. (The reason for that is I'm creating some spatial elements on the fly with users' request) But it's very difficult even to find a place to start... Should I write a python script which can return a Json object and then publish that script as a service? If anyone can give me a hint on where to start, it would be much appreciated!!! Thanks


title: In ArcObjects, how do I change the spatial reference of a Map (Data Frame) in an existing MXD?

text: I need to change the spatial reference of an existing MXD programmatically. In ArcMap UI, I would take the following steps: 1) Open the data frame properties and go to the Coordinate System tab; 2) click on little globe button (Add Coordinate System) and choose New Projected Coordinate System; 3) choose Mercator; 4) enter custom name for the coordinate system; 5) change the Standard_Parallel_1 parameter to the appropriate value; and 6) Save/PressOK. Ultimately, I need to automate this process over hundreds of MXDs. I am really very confused about where even to start. My research has turned up a lot of references to IProjectedCoordinateSystemEdit and related classes and interfaces and to reading/writing parameters to arrays with IParamter and related, but I am finding a lot of the information confusing. So, _which Interfaces/Classes should I start with, and what is the most concise way of accomplishing automation of the process I describe above?_ My code currently opens the MXD to be processed in an IMapDocument and changes a suite of PageLayout and Map settings. What is left is the changing of the spatial reference. Thanks for any guidance you can offer. UPDATE (@Hornbydd) Here is the Python main body code. (BTW, when accessing ArcObjects from Python, one thing that Python can't do is implement interfaces, so responding to events and extending ArcGIS isn't possible. I have yet to find any com objects that it can't access, though.)               f = open(chartprops, 'r')          i = 0     for chrt in f:         if i == 0:  #skip the header             i = i + 1             continue         props = chrt.split(',')         chartnm = props[0]         scale = float(props[1])         latorig = props[2]         pixwidth = float(props[3])         pixheight = float(props[4])         if len(props) > 5:             pprwidth = float(props[5])             pprheight = float(props[6])             pprorienttxt = props[7]         else:             pprwidth = pixwidth / 254.0             pprheight = pixheight / 254.0             if pprwidth > pprheight:                 pprorienttxt = 'LANDSCAPE'             else:                 pprorienttxt = 'PORTRAIT'         if pprorienttxt == 'PORTRAIT':             pprorient = 1         else:             pprorient = 2              #copy template chart MXD to chart output location         #==== WARNING! ====         # If the mxd already exists in the detination folder,         # it will be overwritten         #==================         destmxd = outloc + os.path.sep + chartnm + '.mxd'         shutil.copy(templt, destmxd)              pMapDoc = NewObj(esriCarto.MapDocument, esriCarto.IMapDocument)         pMapDoc.Open(destmxd)         pPgLayout = pMapDoc.PageLayout         pPgLayout3 = CType(pPgLayout, esriCarto.IPageLayout3)         pPage = pPageLayout3.Page         pPage.Units(esriInches)         pPage.PutCustomSize(pprwidth, pprheight)         pPage.Orientation(pprorient)         pMap = pMapDoc.Map         pMap.MapUnits(esriInches)         pMap.SetPageSize(pprwidth, pprheight)         pMap.ReferenceScale(1.0/scale)         pMap.MapScale(1.0/scale)              pMapDoc.Close()         i = i + 1     f.close()


title: Georeferencing declassified CORONA images of coastlines?

text: I'm trying to georeference declassified CORONA images of coastlines in order to asses coastline changes and human impact. In most case, georeferencing is easy, but since images of coastlines can only have GCP's inland of the object of interest (in rare occasions, islands an artificial features seaward), namely, shore line and strip of coastal sands, there will be inherent georeferencing errors. How can I minimize errors? I've thought of two options, but have yet to try them:   1. Georeferencing long swathes of coastline at a time (roughly 2Km long by 200m wide) with a linear/1st order equation (Corona images aren't orthorectified), so any localized errors will have minute effect on the outcome. possible problems: not enough GCP in 200 meters wide coastal area.   2. Short swathes of coastline (roughly 5oom * 1Km), thus probably increasing the potential GCP's inland and reducing local errors. possible problems: cumulative error if I try to use one image as GCP for another. The reason I am inclined to georeference one complete CORONA strip is to save space, and ease workload for my computer. I am using ArcMap 10.2, and have a recent orthorectified high resolution image of the terrain


title: Plotting species distribution on map in R

text: I'd like to plot the distribution of a species on a map of Switzerland, following this tutorial.               library(maps)     library(mapdata)          map('worldHires', 'Switzerland', xlim=c(5, 12), ylim=c(45,48))      From this website I download the shapefile for e.g. _Pinus sylvestris_               psylvestris <- readShapePoly("Pinus sylvestris.shp")      I want to add the range to the map:               plot(psylvestris, add=TRUE, xlim=c(5, 12), ylim=c(45,48))      But I do not see anything (other than the map).. It appears to me the coordinates might be given in a different system (?)               > psylvestris     class       : SpatialPolygonsDataFrame      features    : 52      extent      : -1202142, 4710481, -925131.2, 6334053  (xmin, xmax, ymin, ymax)     coord. ref. : NA      variables   : 1     names       : Id      min values  :  0      max values  :  0       Unfortunately I do not know anything about spatial analysis and thus I have no idea how to fix this. Can you help me please?


title: 502 Bad Gateway

text: I am trying to access the featureInfo in my application and I keep getting 502 Bad Gateway Network error. In teh console log under the response tab I am getting this error,               This proxy does not allow you to access that location (castle.iplantcollaborative.org:8080).      The GET param is,               http://castle.iplantcollaborative.org:8080/geoserver/wms/?LAYERS=opengeo%3APimaPine,opengeo%3ABIEN_Species_Richness_Margalef,opengeo%3ABIEN_Species_Richness,opengeo%3ABIEN_Species_Richness_Margalef_Interpolated&QUERY_LAYERS=opengeo%3APimaPine,opengeo%3ABIEN_Species_Richness_Margalef,opengeo%3ABIEN_Species_Richness,opengeo%3ABIEN_Species_Richness_Margalef_Interpolated&STYLES=,,,&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetFeatureInfo&BBOX=-18295967.087793%2C-10684061.894102%2C18295967.087793%2C10684062.234102&FEATURE_COUNT=10&HEIGHT=546&WIDTH=935&FORMAT=image%2Fpng&INFO_FORMAT=application%2Fvnd.ogc.gml&SRS=EPSG%3A900913&X=174&Y=140      My `/usr/lib/cgi-bin/proxy.cgi` has the domain's url added to it's allowed hosts. The apache error.log file was giving this error,               Permission denied: exec of '/usr/lib/cgi-bin/proxy.cgi' failed      But when I changed the `chmod` for the proxy.cgi file to `root:root` it is not giving that error but it continues to give 502 Bad Gateway. I know I need to make some changes somewhere in the user permissions but not sure where? Any idea what I should correct?


title: Revise from repository on standalone Win7?

text: I'd like to fix the bug in QGIS whereby starting QGIS offline or on a slow connection can kill it, using this revision: http://hub.qgis.org/projects/quantum- gis/repository/revisions/bb623a8c94e49ed478f736a58df49732edd338b4 Problem is, I installed from the standalone installer, not OSGEO4W. What's the best way to go about this, and possible future repository-based installs? Save my files, uninstall QGIS & start again from OSGEO4W?


title: Why does the QGIS 'project properties' option disappear in Ubuntu?

text: For some reason my 'project properties' button has disappeared for good. I'm using Ubuntu 12.04 and QGIS 1.8. Sometimes it would straight up disappear, to return later, but now its just 'gone'. I'll try to re-install, but has anyone else noticed this? I have uninstalled the QT framework (to combat some bugs in QGIS and Scribus).


title: OpenLayers AnimatedCluster delete all features

text: I am using `OpenLayers.Strategy.AnimatedCluster` in my application in order to cluster and show points to the end-user. The example is here. All points are located in concrete vector layer. I need to remove all features from that vector layer, but when I use `removeAllFeatures` or `destroyFeatures` then it removes clusters only from the concrete zoom level and if I zoom in and then back I see all my clusters again. **How can completely remove them from the map?**


title: What is default compress number for jpg2000 when exporting?

text: I have a file geodatabase that has 635 GB Imagery I copied from Mosaic Dataset. I am trying to export it to jpg2000 to a new file geodatabase. I am thinking what is the default compress number for jpg2000 when I export it ? I thought it was 75 percent so I calculate it would compressed to about 158 GB ? Am I correct ? I am using ArcGIS 10.1 SP3


title: OpenLayers make feature not clickable

text: I have defined an `OpenLayers.Select.Control` \- I need only `hoveron` and `hoverout` events. I have implemented it using callbacks that and it works. But I need to disable just on feature clicking functionality. How can I do that?


title: Open Source Tool capable of reading Esri File Geodatabase (.gdb)?

text: I have an ESRI File Geodatabase (`.gdb`) file that I would like to convert to shapefile (or any other non-proprietary format). Are there any Open Source tools to do so?


title: Using ArcGIS 10 on a virtual machine with OSX

text: I'm currently using ArcGIS 10 on a Windows 7 64 bit Macbook Pro (2.53 Ghz and 8GB of RAM) using VMWare. However, when I am using ArcGIS I find it quite slow even with Windows 7 (64-bit) optimized for performance (no shiny graphics). My work involves making some python scripts and testing them, and whenever I run them, that is also slow, but even opening ArcMap and adding a small shapefile takes a long time (30s to 1 minute). Is it possible to make my python program in OSX, import the ArcPy module and somehow run them from the Mac side? I realize that ArcGIS runs on Windows, but I was not sure whether any of the ArcPy part can be run separately. Or, could it be quicker to use the 'Unity' feature on Vmware to run these scripts? I would also appreciate any tips on how much memory/processors to allocate for the virtual machine. At the moment it is 1 processor, and ~4gb of RAM and I would expect it to be quite a bit quicker.


title: GetLegendInfo SOAP call in ArcGIS Server 10

text: I am generating a legend for a Silverlight map using the GetLegendInfo method in the ArcGIS 10 SOAP Web Services. The array of MapServerLegendInfo objects that is returned does not contain annotation layers or raster layers that are in my map service. Does anyone know how I can get those layer types to return, or if there is a good workaround? I know that the ESRI Silverlight API v2.1 contains it's own legend control, however I am not able to use that. Thanks, Seth


title: Most common unmet GIS need for businesses

text: What would you say is the most common unmet GIS need for businesses, as far as technology goes? I'm a programmer and I'm trying to think of a business idea. I like GIS, so I think I might want to do something related to GIS. In order to identify opportunities, I'd like to know what the most common unmet GIS needs are for businesses.


title: Merge several DEM into one in ArcGIS 10

text: Is there any ArcGIS 10 tool to merge several DEM into a single one?


title: How to join a WFS layer to a stand-alone table in OpenLayers

text: I'm a complete novice at OpenLayers, and am checking whether the following is possible. I envisage a dataset of polygons in GeoServer, being served as a WFS layer into OpenLayers. The polygons won't have any attributes other than a unique identifier (eg State name). The attributes for the polygons will be provided by a stand-alone table, which should match to the polygons via the unique identifier. The join needs to be configured within the user's browser (rather than on the server) since the table will be unique to that user and that session - the table is generated dynamically as the result of the user's actions. So my questions are:   * is it possible to join a WFS layer to a stand-alone table?   * is there any sample code showing how to set this up? Thanks


title: SuperMap GIS in depth assessment

text: Does anybody have any idea about SuperMap GIS? \--What are its limitations and strengths? \--How does it compare to other GIS (especially ArcGIS)? \--Can a local Government Department rely on it?


title: 60GB .gdbtable file in my file geodatabase – what is it, and can I delete or make smaller?

text: The title says it all really – I have a single file named 'a0000000.gdbtable’ in my default file geodatabase that is 60 GB in size. I have a number of other files with .gdbtable extension that aren't nearly as big. Just wondering what it is and whether I can delete it or make it smaller? It's taking up too much space!


title: openlayers - wms - map not loading - what Bounds/Extent/Scale to use?

text: I simply want a map on a webpage where I can zoomin and zoomout and navigate around. I'm working on this for 4 days now, and can't get it to work: I want to display the map using a WMS-service here in Austria with Openlayers. Here's the getCap XML:               http://gis.lebensministerium.at/wmsgw/?key=6ae2a41178c77ccee1db7248c9f92b04&VERSION=1.1.1&REQUEST=GetCapabilities&SERVICE=WMS      so this is my openlayers code:               function init() {         map = new OpenLayers.Map("map", {             controls : map_controls                 });         wms = new OpenLayers.Layer.WMS("GIS", "http://gis.lebensministerium.at/wmsgw/", {             layers : 'Orthophoto,Luftbild_MR,Luftbild_1m,Luftbild_8m,Satellitenbild_30m',             key : '6ae2a41178c77ccee1db7248c9f92b04'                 }, {                     projection : new OpenLayers.Projection('EPSG:31287'),                     maxExtent : new OpenLayers.Bounds(110000, 250000, 690000, 590000)                 });         map.addLayer(wms);         if(!map.getCenter()) {             map.zoomToMaxExtent();         }}      and here's the link to take a look at the actual page: myMapSite but it should look something like this (made it with an desktop app): Link_to_Picture (I think the small rectangles are due to the free service and to not use it commercially) So basically I want to see the whole map and then zoom in. Since the only way to see a map (although the zoom-factor is way off) was to set the maxExtend with the boundingBox values from the getCapabilities-xml file, I'm assuming that my Bounds/Extent/Scale is wrong? But what values should I use? or is there another problem? I'm really desperate...


title: How to get a AddIn to appear in customize mode categories?

text: I was hoping that someone here might have tried this, it seems like it could be a simple error somewhere. I was running the sample AddIn that was provided by ESRI and it worked fine until I tried to add it to ArcMap. Once I added the Simple_Add_In esriaddin file to ArcMap it would compile and I would be able to add the toolbar component of the sample but all the buttons on the toolbar would appear as "missing" instead of buttons as well as the commands would not appear in the customize menu, and I cannot add them as buttons. I have installed ArcObjects into Eclipse IDE.


title: How can I apply a raster style to multiple layers using MapServer?

text: I'm using Mapserver to display a raster in browser using a WMS request from OpenLayers. At the moment I have just one raster, and used the following method of styling it:               MAP               ###########         NAME         "ny_pop"                                                                                UNITS         meters             SIZE          256 256         STATUS        OFF         DEBUG         ON         IMAGECOLOR    128 128 255         CONFIG        "MS_ERRORFILE" "ms.log"              ###########             PROJECTION             "init=epsg:900913"         END              ###########         OUTPUTFORMAT           NAME "png"           DRIVER AGG/PNG           MIMETYPE "image/png"           IMAGEMODE RGB           EXTENSION "png"           FORMATOPTION "GAMMA=0.75"         END              ###########             WEB                          METADATA                 "wms_title"           "WMS population layer"                 "wms_onlineresource"  "http://mysite.com/mapserv.cgi?map=population.map&"                             "wms_enable_request"  "GetMap GetFeatureInfo GetCapabilities"             END         END              ###########             LAYER             NAME      "new_york_pop"             DATA      /home/mysite/webapps/htdocs/data/ny_pop_den.tif             STATUS    OFF             TYPE      RASTER                METADATA                 "wms_title"  "WMS NY Raster"                 "wms_srs"    "EPSG:900913"             END                  CLASSITEM "[pixel]"                  # class using simple string comparison, equivelent to ([pixel] = 0)             CLASS                 EXPRESSION "0"                 STYLE                   COLOR 0 0 0                 END             END                    # class using an EXPRESSION using only [pixel].             CLASS                 EXPRESSION ([pixel] > 0 AND [pixel] < 5000)                 STYLE                   COLOR 215 25 28                 END             END                  CLASS                 EXPRESSION ([pixel] >= 5000 AND [pixel] < 10000)                 STYLE                   COLOR 253 174 97                 END          END          END      I want to include several more rasters layers in this mapfile. Is there a way that I can apply the style to multiple layers by creating a style object and applying it to each layer? I'm not sure if this is the correct terminology as I could not spot how to do this on the mapserver site. I thought about generating this using php, or copy/pasting this layer structure. If I go the copy/paste route the mapfile would become extremely long and hard to read.