title: Create image (TIFF/PNG) from geodetic (lat/long) data

text: Using Java, I need to render a bunch of geodetic data to an image, ideally TIFF or PNG. Specifically this is weather data, with each sample having a latitude, longitude, and amplitude value. I'd like to avoid reinventing the wheeling by converting polar coordinates (lat, long) to Cartesian coordinates (x, y in `BufferedImage` or `Graphics2D`). Furthermore the amplitude value needs to rendered as a color gradient. I've looks at GeoTools and NASA's World Wind, but it's unclear and not obvious if it has this capability. Is there a Java library that can do this?


title: Removing multiple categories from GRASS vector feature

text: Pardon my lingo if this is not entirely correct - I am a software developer, not a GIS person. We're working a project that uses `v.overlay` to merge a map of slopes with a map of soil tests to produce a sloped soil test map. There's a number of steps involved around this, but this where it seems to go wrong. What I'm ending up with is a DBF file full of data, most of it seems to be duplicated. Columns copy over and I end up with a_*, b_*, etc... This is fine, but I don't know how to detach a category from a feature. Moreover I'm not really sure how these are associated in the first place. If I go in and delete all the data from the DBF I still end up with duplicate polygons (area type?) that I'm unable to remove. Basically I just need pointers on how to clean up a vector map after using `v.overlay` as this seems to be the source of the problems. Sorry if this is not clear. Ask for more information about something specific if it would be helpful and I'll see what I can provide. Thanks! ![Duplicate categories in GRASS 6.4](http://i.stack.imgur.com/GBU9r.png)


title: Geographically georeference image with limited information

text: It is possible to semi-accurately georeference an image given only the following information? I'd like to do this programmatically.   * Center point of terrain, Lat/Lng (WGS84)   * Pitch of camera   * FOV   * Altitude My intuition says there's not enough reference points, but I hope I'm wrong. If it's possible, any pointers to the math would be appreciated!


title: QGIS Spatial join many points to single vector

text: I've got two shape files, one contains many point coordinates and the other vector shapes. There are many point coordinates contained within each vector shape. I would like to perform a spatial join and create a new shapefile of the coordinates but with the attributes of the overlying vector shape. QGIS spatial join does not seem to work for this. Any ideas? Also setting up a spatial database would be the ideal solution but is not possible in the time frame required. Thanks Richard


title: Can i extract data from different vector layers with QGIS?

text: I'm working on a project where I have a point vector layer and a polygone vector layer I have demographic information associated with the polygone layer whereas the point layer is only localization of stores I have created buffers around each store as a "potential sales zone" but now I would like to see in an attribute table (or else) the list of stores with the data of each correponding polygone or here could be buffer Is it possible? and how do I do it? I'm using QGIS 2.0.1 Thank you


title: Plugins to guide inexperienced users: updating shapefiles and managing data

text: I am attempting to find a way to block adding features while certain attributes are not yet added. This would be useful in order to guide inexperienced clients and employees with the process of adding entities and updating shapefiles while managing data. For example, add feature in edition mode disables 'Okay' button until certain attributes are filled. Any suggestions on either a plug-in that exists or where to start for programming one?


title: ArcGIS legend columns

text: I am making a legend in ArcGIS 10, and have a few layers with a large number of classes. When I try to display the legend with two columns for that layer, the second layer always lines up with the heading/layer name. This doesn't look proper, and I want the first item in each column to line up with each other. I have used a few work-arounds, like splitting the legend, or converting to graphics and manually adjusting. Is there a place in the legend properties where this can be done?


title: Using ArcGIS 10 ModelBuilder to perform In-line variable substitution for input data path?

text: Does anyone know how/if you can use in-line variable substitution for an input path name in ModelBuilder for ArcGIS 10? I used "get field value" to retrieve a code from one feature class, and I want to substitute that text string in two places in another input variable pathname. For instance, if I used "get field value" to retrieve the code "ABC123", I would use that to call up a raster (make raster layer) with this path name: "D:\Data\Rasters\%ABC123%\el_%ABC123% Does anyone know if this will work? ModelBuilder won't recognize this as a valid path name right now.


title: Finding zones with nodata in ArcGIS 9.3

text: I have a DEM raster with some small and large areas of NoData. I also have a raster defining zones (the zones consists of a rectangular grid). I want to create a new raster marking the zones where there exists any nodata values in the DEM raster. My initial approach was to use the Zonal Statistics tool, using any statistics type (for example Maximum) and uncheck the "Ignore NoData" option. According to what I can read in the documentation, this should give me an output raster with NoData in the entire zones where the input raster contained any NoData pixels. (I could then perform some trivial operations to give me the desired final result). However, this doesn't work as expected. The output raster from the Zonal Statistics tool have some holes, but there are several places where there are definitely NoData in the input raster and data in the output. So, are there any other possibilities of achieving what I want? Or have I misunderstood how I should use the Zonal statistics tool?


title: How to find layer ID and name from Identify query using ArcGIS Web ADF?

text: I need to retrieve the layer ID of every layer returned by the _IQueryFunctionality.Identify_ function (I am using ADF with .NET) Identify function spits out a _DataTable_ array for layers which have answered to the query. Since there's no way to know what layers have not answered, you need to check the _TableName_ (which is the layer name) for every _DataTable_ in the array. But what if you have two layers with the same name? You can't know what layer answered and what layer didn't The closest answer I could find on the internet is a **link at esri forums** but the topic got disbanded before reaching a solution (which may not exist perhaps). This sounds like a problem to me, can someone please help me?


title: Identifying Coordinate System of Shapefile when Unknown?

text: I have a Shapefile but its coordinate system is Unknown, and there is no *.prj file. How I can identify it now? Are there any tools that could help?


title: Identifying projection for shapefile with unknown coordinate system, no z values and no measures?

text: I am attempting to create a GIS Project with multiple Brazilian cities. I have been collecting, and am interested in, shapefiles that delineate the neighborhood boundaries of said cities. One city's shape file that I have acquired is for Curitiba, Paraná (PR), Brazil. I obtained the city-neighborhood shape file from here: Institute of Research and Urban Planning: Curitiba. You will find the file under "Divisa de Bairros" (Neighborhood Divisions). I am using Quantum GIS to overlay these city shape files over Google Maps, because ultimately I may need to incorporate street or highway information in to my analysis. **_I am ambivalent ultimately about which source I use, and am open to solutions using Open Streets if possible._** I have read the GIS help files and followed the process to view the metadata in ArcMap to see what the coordinate system and/or projection system is for this shapefile. However, on reviewing the metadata it says that the coordinates have "No Z Values" and "No measures" and there is not defined coordinate system. **_My question is: How can I fix this coordinate/measurement issue so that I can correctly project the location of the Curitiba city shapefile on either the Google Maps or Open Streets projections so that I may move forward with my project?_** **_Note: I have already contacted the creator(s) of the files and asked if they may provide any additional information. I will post or answer my own question as need be._**


title: List of methods of determining the projection of a shapefile when .prj is missing

text: > **Possible Duplicate:**   >  Best practices for identifying the unknown coordinate system of a Shapefile It would be helpful to mention what information is needed to perform each technique. What are situations in which this may happen, and the best ways to deal with it?


title: Can anyone help me to determine the coordinate system of an image I need to georeference?

text: > **Possible Duplicate:**   >  .shp with unknown coordinate system I need to capture spatial information off of an old map and I can't determine the coordinate system of the image. I also need to examine/cross-reference some information via overlay. The subject area is Kirkland Lake Ontario, the scale of the map is 1 inch to 400 feet, and an example of the coordinates shown is **86,000E - 90,000N** (on the southwest corner of the map). I took a few guesses at the coordinate system (as did my co-worker) and we haven't had any success. We figured it might be something of the old NAD27 variety. I also tried using the road, rail, utility line and river layers from the province to georeference the old map, but I'm not having good results that way either. Accuracy is preferred too, so using the grid to georeference this map would be the preferred method. I fear this map may have been generated with some sort of custom coordinate sytem that the local mining companies were using many decades ago. Any ideas? Thanks in advance.


title: Guessing a map projection by its values

text: This is more a theoretical question and not an actual problem that I have.   **Is it possible to identify or classify projections by their numerical coordinates?** Or alternatively rule out in some sort of a decision tree which projection some coordinates have? UTM and Lat-Long projection are normally easy to identify. What about the others?   Of all the possible (and likely many more) projection forms: http://en.wikipedia.org/wiki/List_of_map_projections Assume the situation that you got a shapefile containing points without any .prj file or other information from the data provider. Is there some sort of pattern I could follow to rule out certain projections from the start? Maybe someone can give some insight into the diversity and properties of coordinate projections.


title: Using Affine Transformation to expand polygons in QGIS

text: I have a shape file, and I'd like to grow each feature by 13.7% or 14%. The program crashes every-time I try to use affine transformation. Any recommendations? I don't have a license for ARCgis. Link to the shapefile is here https://www.dropbox.com/sh/frys1pur7g1mts1/ZlfKXM-Go6


title: Equal (ish) areas Map Projection at the Tropics

text: What map projection (non-interrupted) preserves areas well for the entire tropical region? I am thinking of using a Cylindrical Equal Area or a sphere based version of a Cylindrical Equal Area or Mollweide (http://spatialreference.org/ref/esri/54009/). My goal is to minimize area distortion.


title: How to crop coverage with WPS chain in sld?

text: I had 2 layers:               custom:Temperature: Point layer.     custom:Region: Polygon layer.      I would like to perform a gs:CropCoverage of the Raster Layer which is the result of gs:BarnesSurface of vector layer 'custom:Temperature' against the FeatureCollection collected from the vector layer 'custom:Region'. ![Process Chain](http://i.stack.imgur.com/Ts2DI.png) The SLD file is below:                   <?xml version="1.0" encoding="ISO-8859-1"?>         <StyledLayerDescriptor version="1.0.0"          xsi:schemaLocation="http://www.opengis.net/sld StyledLayerDescriptor.xsd"          xmlns="http://www.opengis.net/sld"          xmlns:ogc="http://www.opengis.net/ogc"          xmlns:xlink="http://www.w3.org/1999/xlink"          xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">           <NamedLayer>             <Name>CropCoverage</Name>             <UserStyle>               <Title>Crop Coverage</Title>               <FeatureTypeStyle>                 <Transformation>                   <ogc:Function name="gs:CropCoverage"> <!-- crop coverage -->                     <ogc:Function name="parameter">                       <ogc:Literal>coverage</ogc:Literal>                            <ogc:Function name="gs:BarnesSurface"> <!-- barnes surface -->                         <ogc:Function name="parameter">                           <ogc:Literal>data</ogc:Literal>   <!-- parameter data: should be the layer custom:Temperature -->                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>valueAttr</ogc:Literal>                           <ogc:Literal>TEMP</ogc:Literal>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>dataLimit</ogc:Literal>                           <ogc:Literal>2000</ogc:Literal>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>scale</ogc:Literal>                           <ogc:Literal>0.05</ogc:Literal>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>convergence</ogc:Literal>                           <ogc:Literal>0.2</ogc:Literal>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>passes</ogc:Literal>                           <ogc:Literal>3</ogc:Literal>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>minObservations</ogc:Literal>                           <ogc:Literal>0.1</ogc:Literal>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>maxObservationDistance</ogc:Literal>                           <ogc:Literal>1</ogc:Literal>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>pixelsPerCell</ogc:Literal>                           <ogc:Literal>4</ogc:Literal>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>outputBBOX</ogc:Literal>                           <ogc:Function name="env">                             <ogc:Literal>wms_bbox</ogc:Literal>                           </ogc:Function>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>outputWidth</ogc:Literal>                           <ogc:Function name="env">                             <ogc:Literal>wms_width</ogc:Literal>                           </ogc:Function>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>outputHeight</ogc:Literal>                           <ogc:Function name="env">                             <ogc:Literal>wms_height</ogc:Literal>                           </ogc:Function>                         </ogc:Function>                         <ogc:Function name="parameter">                           <ogc:Literal>queryBuffer</ogc:Literal>                           <ogc:Literal>0.5</ogc:Literal>                         </ogc:Function>                       </ogc:Function> <!-- end of barnes surface -->                          </ogc:Function>                          <ogc:Function name="parameter">                       <ogc:Literal>cropShape</ogc:Literal>                       <ogc:Function name="gs:CollectGeometries"> <!-- collect geometries -->                         <ogc:Function name="parameter">                           <ogc:Literal>features</ogc:Literal> <!-- parameter features: should be layer custom:Region -->                         </ogc:Function>                       </ogc:Function> <!-- end of collect geometries -->                     </ogc:Function>                   </ogc:Function>                 </Transformation>                 <Rule>                   <Name>rule1</Name>                   <RasterSymbolizer>                     <Geometry>                       <ogc:PropertyName>STATION_GEOMETRY</ogc:PropertyName>                     </Geometry>                     <ColorMap type="ramp">                       <ColorMapEntry color="#000084" quantity="-40" label="label" opacity="0"/>                       <ColorMapEntry color="#0000f1" quantity="-32" label="values" opacity="1"/>                       <ColorMapEntry color="#004dff" quantity="-24" label="label" opacity="1"/>                       <ColorMapEntry color="#00b1ff" quantity="-16" label="label" opacity="1"/>                       <ColorMapEntry color="#29ffce" quantity="-8" label="label" opacity="1"/>                       <ColorMapEntry color="#7aff7d" quantity="0" label="label" opacity="0"/>                       <ColorMapEntry color="#ceff29" quantity="8" label="values" opacity="1"/>                       <ColorMapEntry color="#ffc400" quantity="16" label="label" opacity="1"/>                       <ColorMapEntry color="#ff6800" quantity="24" label="label" opacity="1"/>                     </ColorMap>                   </RasterSymbolizer>                 </Rule>               </FeatureTypeStyle>             </UserStyle>           </NamedLayer>         </StyledLayerDescriptor>      The data of gs:BarnesSurface should be layer 'custom:Temperature' and feature of gs:CollectGeometries should be layer 'custom:Region'. The question is: How should i set 'data' and 'feature' when i request wms with the SLD above?               .../wms?layers=custom:Temperature,custom:Region&style=CropCoverage&... ?      Thanks


I have a polygon, imported from a shapefile in QGIS. I would like transform it to a linestring instead so I can save it as a GPX file later. How can I transform a polygon to a linestring in QGIS?


title: Only Draw Geometric Objects Above 8th Parallel North

text: I have many geometric objects (in both SRID:4326 - the_geom - and SRID:3857 - the_web_geom - in PostGIS extension tables of PostgreSQL), including country boundaries, for the country of Ghana and I want to make a hard cut-off of the 8th parallel north latitude - i.e. I only want to draw the_web_geom (SRID:3857) object if it falls above this line. Is there a way in either Mapnik, by adjusting the extent property of the layer, or by amending the PostgreSQL query to only draw the geometric features that fall above this line of latitude? I believe one solution, for the point objects, would be `...WHERE ST_Y(the_geom) > 8.0`, but not sure if this is the best approach or how this would work with polygon objects.


title: Which type in entities for PostGIS/Hibernate Spatial

text: I already manage to store, read and update PostGIS rows with Hibernate Spatial with a Point and a MLineString in a very simple case : a MLineString defining a road, display a point along this road function on the time (the M dimension). However, I'm still confused about which type to use. There is the model packages from PostgreSQL driver, PostGIS plugin (PGxx), Hibernate, HibernateSpatial, JTS ... What type set from which library should I use in my entities.


title: How to use arcpy to ungroup a group layer in an mxd

text: I am looking for a way to ungroup a group layer using arcpy... (instead of having to open the mxd and right-click on the layer to choose "Ungroup"...)               mxd = arcpy.mapping.MapDocument(pathToMyMxd)     df = mxd.activeDataFrame          newLayer = arcpy.mapping.Layer(pathToMyLayerFile)     arcpy.mapping.AddLayer(df, newLayer, "BOTTOM")          mxd.save()          layers = arcpy.mapping.ListLayers(df)     for layer in layers:         if layer.isGroupLayer:             # layer.ungroup()???


title: Calculating distances between points two feature datasets

text: I will do my best to explain my problem. I have two feature datasets. One set has about 100 points that are "city centers." The other dataset has points for each census tract in each of these 100 cities. THere is a row of data per census tract as well as an identifier number for each city (which is also in the city center feature dataset). I need to measure the distance between each city center and each point/census tract in that city and then have these distances in a table with the census tracts numbers....I was hoping there was some way to use python (I'm a beginner at python!) to loop through each point and measure the distance to the city center. ANY suggestions on how to best approach this would be great and the more detailed your response, the better! I've used point distance tool before when I only had 5 cities so I could just run each one manually to get the distances (each city was also a separate file)....


title: Will ActiveSync 4.1 work with ArcPad 10.0.2?

text: Because Microsoft's ActiveSync 4.5 will not work on machines with Windows XP (SP3), from what I've read I need to install ActiveSync 4.1. Has anyone tested this on ArcPad 10.0.2?


title: GIS "hack": ways arbitarily assigning coordinate in order to match the raster and vector layer in QGIS

text: Before I raise my question, I have to admit that I am a total newer on GIS, I don't have much background on GIS. But suddenly I need to take over the GIS job from a GIS specialist who just left the office. And since time is a luxury for me right now, I am looking for some "hack" to fix my immediate problem. The task is below: I need to spot out which bed in an elderly home belongs to a male or female resident, and present it visually on a map. The materials are below: - A raw map from the elderly home, marked the bed no for each bed - A spreadsheet contain the bed no and the gender of the resident on that bed What I have done till now: - "digitize" the raw map: I have used MapWindows and created a vector layer consist of beds, something like creating a new layer on photoshop, and then trace the bed on the new layer - created a new project in QGIS, import the raw map as raster layer, the bed layer (vector layer created by MapWindows), and the spreadsheet as csv - in QGIS, I have joint the bed layer with the csv file correctly - made adjustment in the "Style" tab of the bed layer (vector) so that the gender of each bed can be shown I think I am almost there, but problem arise: - in Mapwindows, both the raw map and the vector layer matched perfectly (i.e. I draw the trace the bed on the new layer and both the raster and vector layer align to each other perfectly) - this is not the case in QGIS I am thinking of something like "georeferencing" issue, which the GIS specialist keep saying before, but I don't know what the hack it is. I need to "realign" the raster and vector layer, how can I do it? Thanks.


title: Drawing a circle of a given area in MapInfo

text: I would like to draw a circle in MapInfo exactly 0.9Ha in area. I've done it by working out the corresponding radius should be, drawing a circle of arbitrary size, double clicking on it and then entering the radius. **Is there a tool to do this?** If there isn't a tool for this then I'm going to make one, so I have another question: **how do you programmatically get the geographic coordinates of a mouseclick?**


title: How to update a Mapinfo layer at certain time intervalls?

text: I have point which represents stress signal. The signal changes at certain time. How do I implement this in mapinfo ? I did this using time() function and increment it to certain seconds but whole application becomes busy which it is running the procedure. Is there other way? question on mapinfo-l


title: ArcGIS Conversion Error While Converting ED50 Datum to WGS84 Datum

text: With ArcGIS for Desktop, I tried a datum conversion of Transverse Mercator ED50 Datum to Transverse Mercator WGS84 Datum. Normally, the difference between two projections is approximately 180 meters but when the conversion is completed with ArcGIS' "Project" tool, the shift of the two features becomes 80 meters, not 180 meters. On a projected data frame, ArcGIS shows the exact same place for two datums, which is an expected issue, but in fact the WGS84 one is not at the correct location. After completing this process with a different software I opened two features on the same data frame with an ED50 datum. Features are not on the exact locations and there is approximately 100 meters visualization error (it means, the datum conversion error of the ArcGIS is 100 meters): ![True datums and coordinates \(Data Frame is ED50\)](http://i.stack.imgur.com/S07JF.png) A data frame without a projection and datum shows the real difference between the datums, which is approximately 180 m: ![enter image description here](http://i.stack.imgur.com/dNZ0b.png) Is this a known issue with ArcGIS or is there a solution for converting ED50 to WGS84 with "Project" tool of ArcGIS correctly?


title: How to write SLD rules with data-defined symbol size and color?

text: I have shapefile with points. They are classified - Graduated colors. I made SLD file with ArcMap2SLDConverter. I have a task to manual edit SLD file to show symbol diameters relating to field Area. I need to create one rule (for symbol COLOR) relating to one attributes ("Class" column), and second rule (for symbol SIZE) relating to other attributes ("Area" column) for the same geometry. So now I have points with rules for symbol colors, but I don't know how to edit my .sdl file to insert rules for symbol size considering they are related to other attribute than those rules for symbol colors. I don't know how to have two rules relating to two different attribute columns for the same shapefile? Can someone help me with that task? Thank you.


title: Is it possible to leave a field null when using the calendar widget in QGIS?

text: I'm creating and editing data with QGIS, the data is stored in a Spatialite database but I've run into this issue with shapefiles as well. I have a field I'm storing dates in and I want to use QGIS' calendar widget to populate that field when there is a legitimate value to be stored for that record, otherwise I need the field to stay NULL. The problem is the widget is populating the field with the current date by default. There's no way on the data entry form to clear it out and leave the field NULL. Am I missing something, perhaps a checkbox somewhere?


title: Using multiple wildcard parameters with arcpy.ListFeatureClasses

text: I know there has to be a way to do this but I'm having trouble figuring out how to do so. I would like to create a list of feature classes to perform operations on however, I want to have the list contain feature classes which have 2 separate and distinct file name identifiers. The code would look something like:               FCS = arcpy.ListFeatureClasses("X_*" OR "*_Y")      I know OR is incorrect however, is there some way to do this without having to re write code to work with the second wildcard. Thanks for helping gout, i know there is a selection method in the Select By Attributes which allows for multiple wildcard use using the OR statement.


title: Is NAD 83 the same as EPSG:4326?

text: I'm getting mixed messages about these projections. For example -Is this projection below, NAD 83, 4326?               GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",     SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],     UNIT["Degree",0.0174532925199433]]


title: What projection to use with Mapbox base/reference layers?

text: excuse me if this has been asked before, but I can not find answer after heavy googling. What coordinate system/projection does shapefiles need to be in to be alligned correctly with mapbox reference layers? This page states that the only projection supported is web mercator. However, when I project my shapefile to web mercator in ArcMap, the shapefile is displayed over Africa? Am I simply missing false easting and northing values? Anyway, HERE are examples of alignments with two different projections. As displayed above I've convert my shapefiles into different projections, without managing to align them with the reference layer. Help please?


title: How to assign attributes from hydrologic sinks to upstream features?

text: I have a large hydrology geometric network with more than 300 sinks and roughly 150 barriers. I'd like to have the SinkID added to the attributes of the barriers (some individual networks have a single barrier, some have multiple) for later analysis but haven't seen any obvious way of doing so. I could assign the SinkID to the first stream segment to which it's attached but some barriers may be 10 or more stream segments upstream. ![enter image description here](http://i.stack.imgur.com/uwpHe.jpg) Thoughts?


title: Use Round in ArcMap 10.x Field Calculator on field with both numeric and NULL values

text: With older versions of ArcMap (9.3, 9.2, etc) I was always able to round numbers down to 2 decimal places using field calculator even if that field had both numeric values and a bunch of NULL values. It used to just ignore the NULL values. With ArcGIS 10.x this is no longer the case. `round([UPSTREAMINVERT],2)` Rounding a field that has any NULL values leads to error 999999, error executing function - invalid use of Null: 'round' I realize it works just fine when you remove NULL values from the selection. However, I'd prefer not to have to do this for every field I'd like to have rounded. Has anyone else experienced this and found a work-around? I've tried to do an if-then statement in field calculator to only return the rounded value if the value is not NULL. This gives the same error message. Thanks.


title: Should I use OpenLayers or Leaflet?

text: I was debating with one of my collegues on OpenLayers v/s Leaflet. I made a point that OpenLayers is much better API if we wish to build a project, where you need direct connectivity to the Geoserver and PostGIS. Then I found Open Data Kit, which looks pretty new but has the features of connectivity with the Geoserver and PostGIS. So my project details are as follows,   1. Use the map interface to fetch Feature Info   2. Create a customized tool that takes the lat/lon from user as to where he/she clicks on the map and then fetches the Climate Data from the raster (which is handled by a py script on the server)   3. Allows user to upload excel, which is sent to the py script, which returns a GeoJSON, which creates Vector Features on the map   4. Allow user to create vector polygons, which will fetch the Features it intersects from the WFS Layer.   5. Fetches Layer from the PostGIS Datastore on GeoServer and displays the layers on the map. So now I am confused on which is better and why using OpenLayers over Leaflet makes more sense or not?


title: Seeking guidance for ArcGIS customization with VBA

text: I am looking for guidance, suggestions, tips, or insights in performing a sequence of operations using VB macros for ArcGIS. **A) Creating a new table** How to write a macro in VB that creates a new table with its associated fields, loads the table to the current data frame, and adds the required records? For instance, a parcel table with fields               P_ID (Long)     L_Use_Code (Long)     L_Use_Name (string)     Owner_ID (Long)      **B) Creating a new feature class** How to write a VB macro that creates a new shapefile feature class called “Parcel” and adds it to the map window? **C) Relating and joining data to layers** How to write a macro in VB that relates a table (as in (A)) to a layer (as in (B))? **D) Identification** How to write a macro in VB that shows a feature cursor operation to display the values of a given field for each row in a feature layer table?


title: Working with temporary files in ArcObjects using Geodatabases?

text: My current add-ons rely on shapefiles for working with temporary files. For example, if I run a spatial analyst tool, I create a new shapefile to store the output geometry. ArcGIS stores all temporary files in a default GDB. Can someone point me in the direction of some .NET code or resources that would help me to structure my temporary file handling in the same way? It would be nice to have everything in one GDB that could be cleared, rather than in a folder containing many shapefiles and rasters.


title: Openlayers: Check the map content and tell if the feature exists in the map bound

text: In my map there are numberof layers. I want my custom legends to show when the related layer feature is visible in the current map view. I am able to turn on/off the display of the legend according to display scale but I need to hide the legend if the feature is not in the map view even if it is in the range of display scale. I go through the layer collection of the map. I tried to debug when the vector layer feature is visible in the map view but I get layer.features.length as always zero. I need it total in terms of openlayers script. I do not want to utilize GeoExt. Can anyone help on this? Thanks


title: intersection between a geometry and a layer in Arcobjects VBA

text: i have to make an intersection between a polygon and a layer of parcels, the polygon used in the clip is not belonging to any layer and the result of this intersection should be all the parcels inside and cuted by the polygon (clip) so i thought to use ITopologicalOperator but i haven't understood how to work with this tool , also it seems to give one feature as result but in my case i need a set of features. can you give me the right tool to use in Arcobject for my case ? any advice is welcome :) thanks in advance


title: Would it be possible to use geohash for proximity searches?

text: I'm looking to optimize point proximity geo searches time. My input is lat,lng point and I'm searching on a precomputed set of locations to n nearest points. I don't care how much time/space the building of the precomputed index of locations will take but I do care the queries will be super fast. I'm thinking about using geohash as the search key, where I would first check if I get results for X chars of the key and then continue to trim down chars from the end of the key until I start to see results. To my (very sparse for now) understanding of geo index techniques this approach should be able to produce the fastest results (in terms of query time) compared to all other known implementations (such as R Tree and co.) // this post should probablly be tagged: geohash, hilbert-curve, research


title: creating one tree to many observations in qgis

text: Looking to create a database for QGIS where each point represents a tree (with name attribute) and another table of observations about each tree over time. Ideally I can update tree location in QGIS without breaking the relationship (sometimes they are inaccurate).


title: Find all of the ways that intersect a linestring

text: I'm new to Open Street Map and PostGIS, and am trying to figure out how to find all of the ways that intersect a linestring (i.e. all of the streets that intersect a route). Is there a way to do this that doesn't involve calling ST_Intersects on every row of the table?


title: OpenLayers map initialised but not displaying

text: I am using OpenLayers with Google maps. I have created my map object successfully - there are no errors raised and console.log(map) shows that it is indeed an OpenLayers class - but nothing is showing on screen. All I am seeing is a white rectangle. Does anyone know what might cause this? Any tips for finding the bug that is causing this? With no errors or obvious problems I don't know where to begin. It should be working!


title: How is ST_PointOnSurface calculated?

text: The PostGIS documentation states that `ST_PointOnSurface` returns "a POINT guaranteed to lie on the surface". It seems like this function could be trivially implemented to give results that satisfy the documentation but provide little real-world utility, though I'm certain that PostGIS provides a non-trivial implementation. This introduction to PostGIS provides a nice comparison and contrast of `ST_Centroid` with `ST_PointOnSurface` and says that "[ST_PointOnSurface] is substantially more computationally expensive than the centroid operation". Is there a more thorough explanation of how `ST_PointOnSurface` is calculated? I've been using `ST_Centroid`, but have encountered some edge cases in my data where the centroid is outside the geometry. I believe that `ST_PointOnSurface` is the correct substitute, but the function name and documentation leave room for uncertainty. Further, is the computational expense of `ST_PointOnSurface` incurred even if the centroid does lie within the geometry already?


title: KML to Shapefile using GDAL (LIBKML) - MultiGeometry data not transforming

text: I am trying to convert a KML file to Shapefile using GDAL (ogr2ogr). Regular KML files convert fine but KML files containing MultiGeometry features (eg: Geometry containing both Point and Polygon) do not get transformed. The output is simply a blank shapefile containing no shapes or attributes. I have tried this command:               ogr2ogr -overwrite -f "ESRI Shapefile" -where OGR_GEOMETRY='MultiGeometry' <ouputfolder> <input kml>      For an example, KML file which has MultiGeometry, you can see the us_states.kml hosted by Google. US States KML - Google.com Other info: I am using 'GDAL 1.10.1, released 2013/08/26' for windows.


title: Calculating minimum/maximum width for overlay of two polygon strips in ArcGIS or some Opensource GIS

text: Is there some solution for calculating minimum and maximum width for overlay of two polygon strips? When I clip these strips the result is another strip with variable width about 1 km and 20 km long. How can I get minimum and maximum width for the clipped strip in shorter axis? I try script from Calculating Polygon Width in ArcGIS but it doesn´t suite with this issue. Is there some solution for ArcGIS 10.1? example: clipped strip with red line for minimum width, blue for maximum Or is it some solution solving the problem straightway without clipping inputs?


title: Is it possible to create standalone installers for QGIS plugins?

text: Plugin in QGIS are managed by the internal plugin manager, and this requires active (and fast) internet connection. Another way to do it offline is answered here How to install a QGIS plugin when offline? which requires technical knowledge. Is there not any way we can get standalone installers for individual plugins just like the qgis package?


title: how to work with netcdf data in Qgis 2?

text: Once it is displayed as a raster, how to see data of different time steps? eg. meteorology datamonthly mean and the likes... where to change the time step? how to calculate seasonal average, annual variations etc..


title: PostGIS ST_Within returns zero data

text: I am using PostgreSQL 9.3.3 & PostGIS 2.1.1 and I have a table Entries which has a column called "latlng" of geographic type (this column stores entry's latitude and longitude coordinates). I am trying to find all entries that belong to a given polygon - like this:               entries = Entry.find_by_sql("         SELECT *         FROM Entries         WHERE ST_Within(             Entries.latlng::geometry,             ST_GeomFromText(                 'POLYGON((                     #{swLat} #{swLng},                     #{swLat} #{neLng},                     #{neLat} #{neLng},                     #{neLat} #{swLng},                     #{swLat} #{swLng}))',                 4326)         );")      Which translates into the following SQL statement:               SELECT *     FROM Entries     WHERE ST_Within(         Entries.latlng::geometry,         ST_GeomFromText(             'POLYGON((                 45.61490708797995 -121.32521530382634,                 45.61490708797995 -121.32338469617367,                 45.615492910571655 -121.32338469617367,                 45.615492910571655 -121.32521530382634,                 45.61490708797995 -121.32521530382634))',             4326)      Even though I am sure that my table has entries that have coordinates that belong to such polygon, I get zero results back. Any ideas on what could be causing this?


title: How to add an independent panel in the ArcMap layout view?

text: In ArcMap 9.3, I want to make a map in which I start from a wider location (for example, a country) and then, I give some zooms until one specific area. I have all the data I need: country administrative boundaries, state/province boundaries and city boundaries (in Brazil and Paraguay). I already got to make my first panel, but how do I plot an independent "frame/panel" into my ArcMap layout view? I added a screenshot below, where I drew in red (paintbrush) the region I would like to zoom and plot in a separate panel. Also, I would appreciate pointers to any tutorial about how to do this kind of location map (start in broader area and zoom to specific place). Tks. ![enter image description here](http://i.stack.imgur.com/RztAA.png)


title: How to use ArcMap for HydroPower Potential

text: I am trying to determine the sites for small hydropower. I have a elevation dem, rainfall data, Curve number data (from soil and landuse shapefiles). I have already used the DEM to fill sinks, find flow direction and accumulation...up to the creation of catchment polygons and slope. I am stuck here. Is there a way to use this to generate points for potential sites?


title: How to sum up polygon area by category?

text: I have a vector layer of habitats. Each feature belongs to a category and has its area calculated. I'm trying to find out the total area of each type of habitat. I've looked at QGIS "Statistics by categories" which gives everything but the Sum of each category. The GroupStats plugin looks like it would do it, but I can't get the "calculate" button to activate. I'm new to QGIS tools/scripts etc but happy to try to modify an existing tool but haven't yet found where they are or how to do this. If anyone can help either directly or point me in the right direction, I'd appreciate it.


title: Batch create fishnet using FME

text: I want to create a fishnet based on the extent of individual polygons from a list of polygons. I plan on processing each resulting fishnets further. I am maintaining a worldwide data set on Quarter Degree Grid Cells split into country level grids. I have previously used a python script to generate the grid cell sets. My workflow has been as described in an article on my blog. The point which breaks my workflow is first generating a world level fishnet, certainly when I am operating on a 1/8th degree cell level and below. I want to create the necessary fishnets on a country level instead. To support that workflow I need to know the following: Does FME provide functionality to create fishnet polygons with a defined lon/lat starting point? And if so does that method accept scripted inputs so that this can be done based on a vector data source with n countries/regions defined?


title: Cannot install ArcGIS for desktop 10.2.1

text: I've checked the system requirements and meet them all. When I press SETUP (image below) a small black box will come up for a few seconds (on back left hand side of image) saying that the programme is too big to fit in memory. Ive checked how much memory I have and it is more than enough. ![enter image description here](http://i.stack.imgur.com/vyOkX.jpg)


title: SQL-query should return GeoJSON FeatureCollection

text: I'm trying to return a GeoJSON FeatureCollection object from a SQL-query with PostgreSQL 9.3 and PostGIS 2.1. This is my query so far:               SELECT row_to_json(fc) AS p      FROM (SELECT 'FeatureCollection' AS type, array_to_json(array_agg(f)) AS features       FROM (SELECT 'Feature' AS type, ST_AsGeoJSON("wkb_geometry")::json AS geometry      FROM "4547cd00-b570-4504-9796-cf2f61dbebe4") AS f) AS fc;      I'm getting the error:               ERROR:  parse error - invalid geometry      However, if I just run this query:               SELECT ST_AsGeoJSON("wkb_geometry")::json AS geometry      FROM "4547cd00-b570-4504-9796-cf2f61dbebe4";      It works and print out a nice geometry. If I ignore the ST_AsGeoJSON it also prints out a result, but not something I can use as a GeoJSON FeatureCollection, as I need the type and coordinates seperated into JSON.


title: How can I get the circle of a sphere (on Earth ellipsoid) around a point using the Latitude and Longitude?

text: Specifically finding points along the circle. If possible without converting into Northing/Easting beforehand. I'm looking for the fastest possible way to calculate points in Latitude/Longitude in a circle around a point on the Earth ellipsoid. I've seen ways to figure out a circle of a sphere, however not while using Latitude and Longitude angles directly and not on an ellipsoid.


title: Parse GML with Openlayers

text: In our project we do a GET request to a WFS service. This service returns some GML. To parse this GML, we are playing with the GML parser of Openlayers. This is working as long the GML isn't 'complex'. It isn't returning everything. This is how a GML object looks like:               <gml:featureMember>       <mag:keywords gml:id="TopografischElement_Trefwoord.fid--255d1722_146cd32a493_b9c">         <mag:keywordid>1</mag:keywordid>         <mag:keyword>Vlees</mag:keyword>         <mag:keywordtopografischelement>           <mag:keywordtopolink gml:id="1">              <mag:keywordid>1</mag:keywordid>              <mag:topografischelementid>4</mag:topografischelementid>              <mag:typologie1>Wegennet</mag:typologie1>              <mag:typologie2>Straat</mag:typologie2>              <mag:typology1>Road network</mag:typology1>              <mag:typology2>Street</mag:typology2>           </mag:keywordtopolink>        </mag:keywordtopografischelement>      </mag:keywords>     </gml:featureMember>      This the code we use to parse it:               var content = _this.format.read(data);      and _this.format is:               this.format = new OpenLayers.Format.GML();      So when looping through "content", each object has a data property. When looking that property, we can find the keywordid and keyword (in this example 1 and 'Vlees'). The problem is, i need the data below that, i need the mag:topografischelementid. How can i make it parse that data too, and not only the first 2 fields? Thank you


title: Add value in Attributestable without reducing the attributetable

text: I am working with the modelbuilder to create a model which can be used by people who do not have a lot of experience with ArcMap. I have a Shapefile with a huge attributetable. I created a new column for a calculation. What I want to do is set the value in the newly created column, depending on a value in another column and field to a special value. This is how it works (simplyfied). The column a is the column I have already. Depending on the value of a field in a, I want to populate a field in column b with a specific value.               a    b     1         2    OHBN     3    OFX     4     5    HGF     5    HGF     2    OHBN     2    OHBN     1      I tried already the select tool and the the field calculator, but that recudes my attributes table to the selected columns. I dont want that to happen. Any suggestions?


title: QGIS pixel/square cms using scale, dpi

text: The data (vectors) I'm using are 5cm in resolution. I think that resolution refers to accuracy of the data and not the image resolution. I've produced maps using the QGIS print composer and exported them as images. . The resolution I've used for these images are at 100 dpi and the scale is at 410 (1:410). Now, if I want to know how many pixels are per square cms in real world length, how can I do that? Thanks, Kaushik


title: ArcGIS Engine Map Control application

text: ![Map Control Application Snapshot](http://i.stack.imgur.com/Sa3p3.png) I am new to ArcGIS. I am trying to create a map control application but when I am adding my digitized map in map control it is only showing raster data (i.e. picture map) not the digitized features and layers. whereas in TOC control it is showing the name of all the layers as you can see in the picture.


title: Extract Administrative Boundaries into OSM poly Format

text: I need to extract some administrative boundaries from big OSM extract, and put them in .poly file. Are there any tools for this task? I've tried Perl script by Gary68 (http://svn.openstreetmap.org/applications/utils/gary68/boundaries.pl), but it returned something like ERROR: at least one needed node missing in osm file. ERROR: at least one needed way missing in osm file. ERROR: at least one needed node or way missing. Also I've installed osmium, but I don't see tool for my task. What would you advise me? (Basically I need to get a list of all houses on all streets of a certain city, and do it automatically for number of cities)


title: ArcMap 10.2 SplitLineAtPoint setSplitRule not calculating RATIO

text: I am trying to set a length field to use the RATION split policy. I set the split policy for the feature layer but when I use the SplitLineAtPoint_management it does not calculate the RATIO length. In the MakeFeatureLayer_management help there is a warning saying "Geoprocessing tools do not honor geodatabase feature class or table field split policies." Does this mean that the SplitLineAtPoint_management will not use the RATIO to calculate the length? This is what I am trying to do (and thanks for any help):               inLayerName = r'D:\Data\Data.gdb\Roads'     inNodes = r'D:\Data\Data.gdb\intersect'     lengthField = 'LENGTH'     outNameLyr = arcpy.MakeFeatureLayer_management (inLayerName, inLayerName + '_lyr', '', 'in_memory')     outNodePntLyr = arcpy.MakeFeatureLayer_management (inNodes, inNodes + '_lyr', '', 'in_memory')     desc = arcpy.Describe(outNameLyr)     field_info = desc.fieldInfo     for index in xrange(0, field_info.count):             if field_info.getFieldName(index) == lengthField:                     field_info.setSplitRule(index, 'RATIO')                     print(' {0} field Split Rule set to {1}\n'.format(field_info.getFieldName(index),field_info.getSplitRule(index)))     arcpy.SplitLineAtPoint_management(outNameLyr, outNodePntLyr, os.path.join(outGDB,outName) + '_splt')


title: Using the Rule-based style to calculate each grid cell from overlaps

text: I am using the Rule-based styling properties (QGIS 2.2) to determine the colour of grid cells for a merged shapefile. Each grid cell defined by its _ID_ may be overlapped by several layers as defined by _Layer_ and each _Layer_ has its own _Score_ (a single column). Focusing on ("ID" = 80) as shown in the image, I want to create a case statement where this _ID_ has an overall _Score_ of 6 or above which will make that cell become red. I have tried something similar to the following:               case when "LAYER" IN ( 'Site_11' , 'Site_12' , 'Site_18' )     then "SCORE" +"SCORE" + "SCORE" >=6     END      This is wrong but it's the kind of logic I'm following at the moment. Any advice welcome! ![enter image description here](http://i.stack.imgur.com/qCT0c.png) ![enter image description here](http://i.stack.imgur.com/cTvbt.png)


title: Building PostGIS with SFCGAL from source fails

text: I'm trying to build latest PostGIS (2.2dev) with SFCGAL. SFCGAL is installed in /usr/local/lib, and includes are in /usr/local/include/SFCGAL I've run configure as:               ./configure --with-sfcgal=/usr/local/bin      and it finds SFCGAL (even if it doesn't output the version, which should be 1.0.4):               -------------- Compiler Info -------------      C compiler:           gcc -g -O2     C++ compiler:         g++ -g -O2     SQL preprocessor:     /usr/bin/cpp -traditional-cpp -P     -------------- Compiler Info -------------      C compiler:           gcc -g -O2     C++ compiler:         g++ -g -O2     SQL preprocessor:     /usr/bin/cpp -traditional-cpp -P     -------------- Dependencies --------------      GEOS config:          /usr/bin/geos-config     GEOS version:         3.3.8     GDAL config:          /usr/bin/gdal-config     GDAL version:         1.10.0     SFCGAL config:        /usr/local/bin     SFCGAL version:            PostgreSQL config:    /usr/bin/pg_config     PostgreSQL version:   PostgreSQL 9.1.13     PROJ4 version:        48     Libxml2 config:       /usr/bin/xml2-config     Libxml2 version:      2.8.0     JSON-C support:       yes     PostGIS debug level:  0     Perl:                 /usr/bin/perl     --------------- Extensions ---------------      PostGIS Raster:       enabled     PostGIS Topology:     enabled     SFCGAL support:       enabled      Everything seems ok, but it fails during the linking of utilities (e.g. raster loader or, in case raster is disabled, the next one is shp2pgsql):               /bin/bash ../libtool --mode=link \     gcc -I ../liblwgeom -g -O2  -fPIC -DPIC  -Wall -Wmissing-prototypes    shpopen.o   dbfopen.o getopt.o shpcommon.o safileio.o shp2pgsql-core.o shp2pgsql-cli.o   ../liblwgeom/liblwgeom.la -o shp2pgsql  -lc      libtool: link: gcc -I ../liblwgeom -g -O2 -fPIC -DPIC -Wall -Wmissing-prototypes  shpopen.o dbfopen.o getopt.o shpcommon.o safileio.o shp2pgsql-core.o shp2pgsql-cli.o -o  .libs/shp2pgsql  ../liblwgeom/.libs/liblwgeom.so -lc -Wl,-rpath - Wl,/home/giova/bin/postgis/lib     ../liblwgeom/.libs/liblwgeom.so: undefined reference to `sfcgal_solid_num_shells'     ../liblwgeom/.libs/liblwgeom.so: undefined reference to `sfcgal_point_y'     ../liblwgeom/.libs/liblwgeom.so: undefined reference to `sfcgal_point_x'     ../liblwgeom/.libs/liblwgeom.so: undefined reference to `sfcgal_geometry_collection_geometry_n'     ../liblwgeom/.libs/liblwgeom.so: undefined reference to `sfcgal_linestring_point_n'     ../liblwgeom/.libs/liblwgeom.so: undefined reference to `sfcgal_polyhedral_surface_polygon_n'     ../liblwgeom/.libs/liblwgeom.so: undefined reference to `sfcgal_triangle_set_vertex_from_xyz'     ../liblwgeom/.libs/liblwgeom.so: undefined reference to `sfcgal_solid_shell_n'     ../liblwgeom/.libs/liblwgeom.so: undefined reference to `sfcgal_solid_create_from_exterior_shell'     (and other similar undefined references)      What's going wrong with liblwgeom linking? Why liblwgeom has all the SFCGAL references undefined? The SFCGAL building seems ok, and if I run readelf on libSFCGAL.so everything seems ok. I past the readelf output run on liblwgeom.so, which shows the undefined symbols Thanks for your precious help. Giovanni                   00066048  00000107 R_386_JUMP_SLOT        00000000   sfcgal_solid_num_shells         000660d0  00000a07 R_386_JUMP_SLOT        00000000   sfcgal_linestring_add_point         00066174  00001407 R_386_JUMP_SLOT        00000000   sfcgal_linestring_num_points         000661a4  00001707 R_386_JUMP_SLOT        00000000   sfcgal_point_y         000661c8  00001907 R_386_JUMP_SLOT        00000000   sfcgal_geometry_collection_num_geometries         00066260  00002207 R_386_JUMP_SLOT        00000000   sfcgal_triangle_vertex         00066380  00002f07 R_386_JUMP_SLOT        00000000   sfcgal_triangle_set_vertex_from_xy         00066390  00003007 R_386_JUMP_SLOT        00000000   sfcgal_polygon_interior_ring_n         000663d0  00003407 R_386_JUMP_SLOT        00000000   sfcgal_solid_create_from_exterior_shell         00066484  00003e07 R_386_JUMP_SLOT        00000000   sfcgal_triangle_set_vertex_from_xyz         00066490  00004007 R_386_JUMP_SLOT        00000000   sfcgal_version         000664b8  00004407 R_386_JUMP_SLOT        00000000   sfcgal_polygon_exterior_ring         000664ec  00004707 R_386_JUMP_SLOT        00000000   sfcgal_point_create_from_xy         000664f4  00004807 R_386_JUMP_SLOT        00000000   sfcgal_geometry_collection_add_geometry         00066540  00004e07 R_386_JUMP_SLOT        00000000   sfcgal_linestring_create         000665d8  00005807 R_386_JUMP_SLOT        00000000   sfcgal_linestring_point_n         0006663c  00005d07 R_386_JUMP_SLOT        00000000   sfcgal_geometry_is_3d         00066660  00006007 R_386_JUMP_SLOT        00000000   sfcgal_geometry_is_empty         000666a0  00006407 R_386_JUMP_SLOT        00000000   sfcgal_multi_linestring_create         000666b8  00006707 R_386_JUMP_SLOT        00000000   sfcgal_point_x         000666d0  00006907 R_386_JUMP_SLOT        00000000   sfcgal_point_z         000666ec  00006a07 R_386_JUMP_SLOT        00000000   sfcgal_geometry_collection_geometry_n         00066738  00007007 R_386_JUMP_SLOT        00000000   sfcgal_polyhedral_surface_polygon_n         0006673c  00007107 R_386_JUMP_SLOT        00000000   sfcgal_multi_point_create         00066744  00007207 R_386_JUMP_SLOT        00000000   sfcgal_polygon_create         0006678c  00007807 R_386_JUMP_SLOT        00000000   sfcgal_multi_polygon_create         00066790  00007907 R_386_JUMP_SLOT        00000000   sfcgal_geometry_type_id         000667cc  00007d07 R_386_JUMP_SLOT        00000000   sfcgal_polyhedral_surface_add_polygon         000667d4  00007e07 R_386_JUMP_SLOT        00000000   sfcgal_polygon_num_interior_rings         000667fc  00007f07 R_386_JUMP_SLOT        00000000   sfcgal_geometry_collection_create         00066800  00008007 R_386_JUMP_SLOT        00000000   sfcgal_polygon_add_interior_ring         0006680c  00008307 R_386_JUMP_SLOT        00000000   sfcgal_polygon_create_from_exterior_ring         000668b8  00008e07 R_386_JUMP_SLOT        00000000   sfcgal_point_create_from_xyz         000668d0  00008f07 R_386_JUMP_SLOT        00000000   sfcgal_polyhedral_surface_create         0006694c  00009607 R_386_JUMP_SLOT        00000000   sfcgal_triangle_create         0006697c  00009907 R_386_JUMP_SLOT        00000000   sfcgal_triangulated_surface_create         0006699c  00009b07 R_386_JUMP_SLOT        00000000   sfcgal_triangulated_surface_num_triangles         000669d4  00009c07 R_386_JUMP_SLOT        00000000   sfcgal_triangulated_surface_triangle_n         000669e0  00009d07 R_386_JUMP_SLOT        00000000   sfcgal_geometry_delete         000669ec  00009e07 R_386_JUMP_SLOT        00000000   sfcgal_polyhedral_surface_num_polygons         00066a4c  0000a507 R_386_JUMP_SLOT        00000000   sfcgal_triangulated_surface_add_triangle         00066a64  0000a707 R_386_JUMP_SLOT        00000000   sfcgal_point_create         00066a6c  0000a807 R_386_JUMP_SLOT        00000000   sfcgal_solid_shell_n              1: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_solid_num_shells             10: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_linestring_add_point             20: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_linestring_num_points             23: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_y             25: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_collection_num_geometries             34: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangle_vertex             47: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangle_set_vertex_from_xy             48: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_interior_ring_n             52: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_solid_create_from_exterior_shell             62: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangle_set_vertex_from_xyz             64: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_version             68: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_exterior_ring             71: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_create_from_xy             72: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_collection_add_geometry             78: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_linestring_create             88: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_linestring_point_n             93: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_is_3d             96: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_is_empty            100: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_multi_linestring_create            103: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_x            105: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_z            106: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_collection_geometry_n            112: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polyhedral_surface_polygon_n            113: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_multi_point_create            114: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_create            120: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_multi_polygon_create            121: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_type_id            125: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polyhedral_surface_add_polygon            126: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_num_interior_rings            127: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_collection_create            128: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_add_interior_ring            131: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_create_from_exterior_ring            142: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_create_from_xyz            143: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polyhedral_surface_create            150: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangle_create            153: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangulated_surface_create            155: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangulated_surface_num_triangles            156: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangulated_surface_triangle_n            157: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_delete            158: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polyhedral_surface_num_polygons            165: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangulated_surface_add_triangle            167: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_create            168: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_solid_shell_n            543: 0004be90    25 FUNC    GLOBAL DEFAULT   11 lwgeom_sfcgal_version            758: 0004c840   198 FUNC    GLOBAL DEFAULT   11 lwgeom_sfcgal_noop            923: 00000000     0 FILE    LOCAL  DEFAULT  ABS lwgeom_sfcgal.c            989: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_solid_num_shells           1030: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_linestring_add_point           1079: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_linestring_num_points           1095: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_y           1105: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_collection_num_geometries           1162: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangle_vertex           1255: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangle_set_vertex_from_xy           1263: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_interior_ring_n           1286: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_solid_create_from_exterior_shell           1312: 0004c840   198 FUNC    GLOBAL DEFAULT   11 lwgeom_sfcgal_noop           1343: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangle_set_vertex_from_xyz           1348: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_version           1362: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_exterior_ring           1383: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_create_from_xy           1386: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_collection_add_geometry           1412: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_linestring_create           1464: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_linestring_point_n           1492: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_is_3d           1497: 0004be90    25 FUNC    GLOBAL DEFAULT   11 lwgeom_sfcgal_version           1503: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_is_empty           1522: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_multi_linestring_create           1533: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_x           1542: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_z           1549: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_collection_geometry_n           1575: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polyhedral_surface_polygon_n           1577: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_multi_point_create           1580: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_create           1600: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_multi_polygon_create           1601: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_type_id           1618: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polyhedral_surface_add_polygon           1621: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_num_interior_rings           1635: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_collection_create           1636: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_add_interior_ring           1639: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polygon_create_from_exterior_ring           1696: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_create_from_xyz           1703: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polyhedral_surface_create           1743: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangle_create           1759: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangulated_surface_create           1770: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangulated_surface_num_triangles           1788: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangulated_surface_triangle_n           1791: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_geometry_delete           1794: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_polyhedral_surface_num_polygons           1825: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_triangulated_surface_add_triangle           1831: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_point_create           1833: 00000000     0 NOTYPE  GLOBAL DEFAULT  UND sfcgal_solid_shell_n


title: How to increase FUSION-LTK lidar processing productivity?

text: I am processing Lidar data with FUSION-LTK. My study area contains 5 blocks of 1000 las tiles per block. Each block is broken into approximately 4 subblocks of 250 las tiles. I am calculating grid metrics and strata metrics to assess vegetation structure. I wrote simple batch files to process entire blocks (e.g. block 1 subblock1; block 1 subblock2;...) using the command line tools in FUSION-LTK. The problem I am encountering is that I am only able to process one subblock at a time. When I try to simultaneously start processing on another block, the current session of FUSION-LTK stops the processing and moves to the new block. How can I improve the processing workflow so that more than one subblock are processed simultaneously, or at least enhance the multiprocessing capabilities of LTK?


title: What is the closest CRS to "Meridiano de Greenwich"

text: I have a 1906 scanned map that I am trying to georeference in QGIS for Mac. The map has French and English lat-long, indicating them as either "Meridiano de Greenwich" or "Meridiano de Paris." I would like to know what the closest approximation to "Meridiano de Greenwich" for choosing the CRS might be.


title: Retrieve path of running plugin in QGIS

text: In my plugin's directory I've got some XML settings files that I'd like to parse at a particular point during execution of the plugin. Is there a way to get the absolute path where the plugin's .pyc file is located?


title: Buffer Zip Region Polygon with GeoTools

text: What I want to achieve is to buffer a zip region geometry with a distance in meters and than get the coordinates of the new (expanded) polygon. The original geometry polygon is parsed from a kml file. This is the code that I wrote so far: Parse the kml file:               public Geometry getPolygon(String plz) {       String poly = storage.get(plz);            if (poly != null) {         try {           InputStream stream = new ByteArrayInputStream(poly.getBytes("UTF-8"));           Parser parser = new Parser(new KMLConfiguration());           SimpleFeature f = (SimpleFeature) parser.parse(stream);                return (Geometry) f.getAttribute("GEOMETRY");         } catch (UnsupportedEncodingException e) {           e.printStackTrace();         } catch (SAXException e) {           e.printStackTrace();         } catch (ParserConfigurationException e) {           e.printStackTrace();         } catch (IOException e) {           e.printStackTrace();         }       }            return null;     }      Transform the polygon:               public Geometry getBufferedZipPolygon(Geometry zipPolygon, double distanceMeters) {       CoordinateReferenceSystem wgsCRS = null;       CoordinateReferenceSystem distanceMetersCRS = null;       try {         wgsCRS = CRS.decode("EPSG:4326");         distanceMetersCRS = CRS.decode("EPSG:32630");              Geometry targetGeometry = zipPolygon.transform(distanceMetersCRS);              Geometry buffer = targetGeometry.getBuffer(distanceMeters);         buffer.transform(wgsCRS);              return buffer;       } catch (FactoryException e) {         e.printStackTrace();       } catch (TransformException e) {         e.printStackTrace();        }            return null;     }      I have two Questions now. Is this transformation valid and how can I get the new Coordinates from my transformed Geometry Object? Answer: With the help from iant I finally figured out this solution. Parser:               public Polygon getPolygon(String plz) {       String polyString = storage.get(plz);            if (polyString != null) {         try {           InputStream stream = new ByteArrayInputStream(polyString.getBytes("UTF-8"));           Parser parser = new Parser(new KMLConfiguration());           Polygon polygon = (Polygon) parser.parse(stream);           return polygon;         } catch (UnsupportedEncodingException e) {           e.printStackTrace();         } catch (SAXException e) {           e.printStackTrace();         } catch (ParserConfigurationException e) {           e.printStackTrace();         } catch (IOException e) {           e.printStackTrace();         }       }            return null;     }      Transformation:               public Geometry getBufferedZipPolygon(Polygon zipPolygon, double distanceMeters) {       CoordinateReferenceSystem wgsCRS;       CoordinateReferenceSystem distanceMetersCRS;       try {         wgsCRS = CRS.decode("EPSG:4326", true);         distanceMetersCRS = CRS.decode("EPSG:32630", true);              MathTransform transform = CRS.findMathTransform(wgsCRS, distanceMetersCRS);         MathTransform reverseTransform = CRS.findMathTransform(distanceMetersCRS, wgsCRS);              Geometry transformedZipPolygon = JTS.transform(zipPolygon, transform);         Geometry bufferZipPolygon = transformedZipPolygon.buffer(distanceMeters);              return JTS.transform(bufferZipPolygon, reverseTransform);       } catch (FactoryException e) {         e.printStackTrace();       } catch (TransformException e) {         e.printStackTrace();       }            return null;     }      Get the coordinates:               bufferedPlzPolygon = (Polygon) getBufferedZipPolygon(plzPolyStorage.getPolygon(zip), range);     Coordinate[] coords = bufferedPlzPolygon.getCoordinates();


title: *.lyr to *.style conversion in ArcGIS

text: Do you know how it's possible to convert a layer symbology or an *.Lyr file to an *.style file? Thanks, Jérôme


title: How to create a storm swath map?

text: I can't really find the words to explain what I'd like to achieve so I attached a picture. ![enter image description here](http://i.stack.imgur.com/QKIxt.jpg) The top picture was made of two layers, the layer with circles was made/styled using the instructions from an answer to a question that I posted a few days ago (How to scale with radius attribute). The other layer was created by converting that same layer using the Points to Paths plugin. The bottom image is what I want to achieve. For background info, the data is about wind information from typhoon Haiyan that struck the Philippines last November 2013. I want to get/generate the storm swath ( **the lower picture from the image above** ) that kind of look like the screenshot that I took from USAID (image below). ![typhoon Yolanda storm swath](http://i.stack.imgur.com/ltQw9.jpg) How can I achieve this or is this possible? Thanks in advance.


title: How can I remove default zoom in/out [+/-] buttons from openlayers to keep the navigation bar on?

text: I want to remove the default zoom +/- buttons from openlayers. Here is my code                   <script type="text/javascript">         var map, ia_wms;         function init(){     map = new OpenLayers.Map('map', {zoomDuration: 1,projection: 'EPSG:3857'     ***/*controls: [         new OpenLayers.Control.Navigation(),         new OpenLayers.Control.ArgParser(),         new OpenLayers.Control.Attribution()         ]*/***         });         var sr_all = new OpenLayers.Layer.WMS("EOL","http://192.168.0.233/cgi-bin/mapserv?map=/var/www/mstest/plant.map",{'layers':"map_outline",'transparent':true,'format':'image/png'}                ,{                 isBaseLayer: false,opacity: 100                 },                 {singleTile: true, ratio: 1}                 );         var gmap = new OpenLayers.Layer.Google(         "Google Streets", // the default         { numZoomLevels: 25,visibility: false}         ,{             isBaseLayer: true         }         );               map.addLayers([sr_all,gmap]);         map.addControl(new OpenLayers.Control.PanZoomBar());              map.zoomToExtent(new OpenLayers.Bounds(69.728426,22.311238,69.774227,22.356185));         map.setCenter(new OpenLayers.LonLat(69.754227,22.332185).transform(         new OpenLayers.Projection('EPSG:4326', 'EPSG:3857'),         map.getProjectionObject()         ), 14);         }     </script>      when I added, controls: [ new OpenLayers.Control.Navigation(), new OpenLayers.Control.ArgParser(), new OpenLayers.Control.Attribution() ] to the map = new OpenLayers.Map, to remove default zoom +/- buttons, the map is not loading. Please see the image also. ![ Both the navigation buttons exist. :\( ](http://i.stack.imgur.com/s1LVE.png) Both the navigation buttons exist. :( **I want to remove the small- 2 button- one.** I changed the code like this:                       map = new OpenLayers.Map('map', {zoomDuration: 1,projection: 'EPSG:3857'         controls: []         });      and added                       map.addControl(new OpenLayers.Control.PanZoomBar());         map.addControl(new OpenLayers.Control.KeyboardDefaults());      But, the map is not loading. Do not know any syntax error..


title: Fitting a word/tag cloud into geographic boundaries

text: Are there any available tools or a recommend methods for displaying a tag or word cloud (also know as wordle) onto a world map, much like in this example: Map: United States of surnames? (There may be other examples on the NY Times.) Note that that in the example above, word are located by fixed coordinates, whereas text clouding might be implemented with the wordle algorithm (as described on SO) in an automatic manner. So, is there a software that let you put words at random on a given country, or do you think it's best to tweak the wordle original algorithm (e.g., by adding constraints based on predefined frontiers)?


title: MapInfo cannot join after deleted some rows with "del"

text: I tried to join my table with another one, with the update column. But it gives me an error, that it cannot be updated. Actually it does update the table to a certain point. I deleted some rows before I tried to update the table, manually by clickling left on the row and pressed the key "del", so the object got removed, the row still kinda exist, but only in grey. But that is the point where the join stops, from there on there are no further entries. I tried to close and reopen the table, but the grey rows are still in there, any idea?


title: Districts of Israel - free source

text: Do you know of any free GIS sources of the Israeli administrative divisions: districts (mechozot, מחוזות) and subdistricts (nafot, נפות) which rely on official government information? I heard there is a some free US government source, but failed to Google it. **Update (2014-02-11):** The Israeli government has opened a national GIS website. The UI will make your eyes bleed, the interface is only in Hebrew and there is no way to download the source files, only "save map as image". On the other hand, it's the only official source. The map can be zoomed and downloaded chunk-by-chunk and georeferenced using QGis or alike. Depending on the required accuracy, this process can take anywhere between 30 minutes to a few hours and produce good-enough results. Permalink to the Administrative boundaries layer.


title: Which license esri server did I get the license from?

text: I have a esri license server setup with a failover (as described in the ArcGIS documentation). All the ArcGIS clients have the environment variable `ARCGIS_LICENSE_FILE` set to `27001@mainlicenseserver;27001@failoverlicenseserver`. My question is: Is there a way a can see from which license server a ArcGIS client (for instance ArcMap) got the license?


title: Point of Interest coordinate data base

text: Can any one help me to find data a database for Points of Interest (with coordinates) for countries in Europe? I got some data from POI PLAZA, but I am not getting sufficient data for Scandinavian countries. Can any one help me to extract the data?


title: Add multiple fields to attribute table by python, with a prompt window for parameters

text: I know how to add field to arcmap shapefile using Python:               import arcpy     arcpy.env.workspace = "H:\pythonexercise\NZL_rds\NZL_roads.shp"     arcpy.AddField_management("NZL_roads", "LENGTH", "TEXT","","","50")      But how to make the python code prompt a window for "field", "data type" and "length"? Also, what if I want add many fields at once like below:   * field "SPEEDLIMIT" data type as "num" and length "3".    * field "BREADTH" datatype as "num" and length "4"   * field "LANE" datatype as "num" and length "4" It should be dynamic, not fixed like my code.


title: Add hot point in android application based on Arcgis SDK

text: I am working on an application based on Arcgis android SDK. Now I wonder if I can add the hot point like google map and etc. For example, not matter you use google map in pc or your phone, when you touch the text or symbol of the poi, the detail information of the poi will be displayed by inforwindow(pc) or something else. I wonder if we can make it in android?


title: QGIS GISBASE specify location

text: I am just a beginner at this QGIS and GRASS and I am using this QGIS 2.0.1 Dufour which throws a error at the startup that > " **QGIS couldn't find your GRASS installation. Would you like to specify > path (GISBASE) to your GRASS installation?** " I cannot find the root cause of the problem or any way to solve it. > From other forum posts, I can see that many others have seen this problem > but no one has given a clear cut solution to it. If someone can help me with this issue, it will be really helpful. My project work is at halt, so please help me solve it. * _The installer is *_ > QGIS-OSGeo4W-2.0.1-3-Setup-x86 and > WinGRASS-6.4.3-1-Setup Also one more update is that it runs fine on my laptop which intel i7 and 4GB RAM with Windows 8 but not on my desktop which is an INTEL DUAL CORE and 2GB RAM which has Windows 7 Enterprise **UPDATE** _my "postinstall.log"_ **FILE** :               `C:\Program Files\QGIS Dufour\postinstall.log     OSGEO4W_ROOT=C:\PROGRA~1\QGIS Dufour     OSGEO4W_STARTMENU=C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour     OSGEO4W_ROOT_MSYS=/C/PROGRA~1/QGIS Dufour          C:\PROGRA~1\QGIS Dufour>"C:\PROGRA~1\QGIS Dufour"\bin\make-bat-for-py.bat              Generating .bat files for all .py files in C:\PROGRA~1\QGIS Dufour\bin               C:\PROGRA~1\QGIS Dufour>set ICON=C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_osgeo.ico           C:\PROGRA~1\QGIS Dufour>set ICON_CMD=C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_cmd.ico           C:\PROGRA~1\QGIS Dufour>set ICON_TCLTK=C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_tcltk.ico           C:\PROGRA~1\QGIS Dufour>set BATCH=C:\PROGRA~1\QGIS Dufour\bin\grass64.bat           C:\PROGRA~1\QGIS Dufour>textreplace -std -t "C:\PROGRA~1\QGIS Dufour"\bin\grass64.bat           C:\PROGRA~1\QGIS Dufour>textreplace -std -t "C:\PROGRA~1\QGIS Dufour"\bin\grass64           C:\PROGRA~1\QGIS Dufour>textreplace -std -t "C:\PROGRA~1\QGIS Dufour"\apps\grass\grass-6.4.3\etc\fontcap           C:\PROGRA~1\QGIS Dufour>mkdir "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\GRASS GIS 6.4.3"            C:\PROGRA~1\QGIS Dufour>xxmklink "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\GRASS GIS 6.4.3\GRASS 6.4.3 GUI.lnk" "C:\PROGRA~1\QGIS Dufour\bin\grass64.bat" "-wx" \ "Launch GRASS GIS 6.4.3 with wxGUI" 1 "C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_osgeo.ico"            XXMkLink     ver 1.00    (c)2005 Copyright  Pixelab, Inc.          The shortut created as follows          Shortcut path:     C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\GRASS GIS 6.4.3\GRASS 6.4.3 GUI.lnk     Target object:     C:\PROGRA~1\QGIS Dufour\bin\grass64.bat     Arguments;         -wx     Working Directory: \     Description:       Launch GRASS GIS 6.4.3 with wxGUI     Display Mode:      Normal Window (1)     Icon file:         C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_osgeo.ico               C:\PROGRA~1\QGIS Dufour>xxmklink "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\GRASS GIS 6.4.3\GRASS 6.4.3 Old TclTk GUI.lnk" "C:\PROGRA~1\QGIS Dufour\bin\grass64.bat" "-tcltk" \ "Launch GRASS GIS 6.4.3 with the old TclTk GUI" 1 "C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_tcltk.ico"            XXMkLink     ver 1.00    (c)2005 Copyright  Pixelab, Inc.          The shortut created as follows          Shortcut path:     C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\GRASS GIS 6.4.3\GRASS 6.4.3 Old TclTk GUI.lnk     Target object:     C:\PROGRA~1\QGIS Dufour\bin\grass64.bat     Arguments;         -tcltk     Working Directory: \     Description:       Launch GRASS GIS 6.4.3 with the old TclTk GUI     Display Mode:      Normal Window (1)     Icon file:         C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_tcltk.ico               C:\PROGRA~1\QGIS Dufour>xxmklink "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\GRASS GIS 6.4.3\GRASS 6.4.3 Command Line.lnk" "C:\PROGRA~1\QGIS Dufour\bin\grass64.bat" "-text" \ "Launch GRASS GIS 6.4.3 in text mode" 1 "C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_cmd.ico"            XXMkLink     ver 1.00    (c)2005 Copyright  Pixelab, Inc.          The shortut created as follows          Shortcut path:     C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\GRASS GIS 6.4.3\GRASS 6.4.3 Command Line.lnk     Target object:     C:\PROGRA~1\QGIS Dufour\bin\grass64.bat     Arguments;         -text     Working Directory: \     Description:       Launch GRASS GIS 6.4.3 in text mode     Display Mode:      Normal Window (1)     Icon file:         C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_cmd.ico               C:\PROGRA~1\QGIS Dufour>xxmklink "C:\ProgramData\Desktop\GRASS GIS 6.4.3.lnk" "C:\PROGRA~1\QGIS Dufour\bin\grass64.bat" "-wx" \ "Launch GRASS GIS 6.4.3 with wxGUI" 1 "C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_osgeo.ico"            XXMkLink     ver 1.00    (c)2005 Copyright  Pixelab, Inc.          The shortut created as follows          Shortcut path:     C:\ProgramData\Desktop\GRASS GIS 6.4.3.lnk     Target object:     C:\PROGRA~1\QGIS Dufour\bin\grass64.bat     Arguments;         -wx     Working Directory: \     Description:       Launch GRASS GIS 6.4.3 with wxGUI     Display Mode:      Normal Window (1)     Icon file:         C:\PROGRA~1\QGIS Dufour\apps\grass\grass-6.4.3\etc\gui\icons\grass_osgeo.ico               C:\PROGRA~1\QGIS Dufour>del "C:\PROGRA~1\QGIS Dufour"\bin\grass64.bat.tmpl           C:\PROGRA~1\QGIS Dufour>del "C:\PROGRA~1\QGIS Dufour"\bin\grass64.tmpl           C:\PROGRA~1\QGIS Dufour>for %f in ("C:\Users\ab281182\AppData\Local\Temp") do call set TEMPDRIVE=%~df           C:\PROGRA~1\QGIS Dufour>call set TEMPDRIVE=C:           C:\PROGRA~1\QGIS Dufour>cd C:\Users\ab281182\AppData\Local\Temp           C:\Users\ab281182\AppData\Local\Temp>"C:\PROGRA~1\QGIS Dufour\bin\vcredist_2005_x86.exe" /q /t:C:           C:\Users\ab281182\AppData\Local\Temp>if errorlevel 3010 echo Dufour\etc\reboot 1>C:\PROGRA~1\QGIS           C:\Users\ab281182\AppData\Local\Temp>del "C:\PROGRA~1\QGIS Dufour\bin\vcredist_2005_x86.exe"           C:\Users\ab281182\AppData\Local\Temp>"C:\PROGRA~1\QGIS Dufour\bin\vcredist_2008_x86.exe" /q           C:\Users\ab281182\AppData\Local\Temp>if errorlevel 3010 echo Dufour\etc\reboot 1>C:\PROGRA~1\QGIS           C:\Users\ab281182\AppData\Local\Temp>del "C:\PROGRA~1\QGIS Dufour\bin\vcredist_2008_x86.exe"           C:\Users\ab281182\AppData\Local\Temp>"C:\PROGRA~1\QGIS Dufour\bin\vcredist_2010_x86.exe" /q           C:\Users\ab281182\AppData\Local\Temp>if errorlevel 3010 echo Dufour\etc\reboot 1>C:\PROGRA~1\QGIS           C:\Users\ab281182\AppData\Local\Temp>del "C:\PROGRA~1\QGIS Dufour\bin\vcredist_2010_x86.exe"           C:\PROGRA~1\QGIS Dufour>mkdir "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour"      A subdirectory or file C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour already exists.          C:\PROGRA~1\QGIS Dufour>xxmklink "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\MSYS Shell.lnk"       "C:\PROGRA~1\QGIS Dufour\apps\msys\msys.bat" " " \ "Minimal SYStem" 7 "C:\PROGRA~1\QGIS Dufour\apps\msys\msys.ico"           XXMkLink     ver 1.00    (c)2005 Copyright  Pixelab, Inc.          The shortut created as follows          Shortcut path:     C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\MSYS Shell.lnk     Target object:     C:\PROGRA~1\QGIS Dufour\apps\msys\msys.bat     Arguments;               Working Directory: \     Description:       Minimal SYStem     Display Mode:      Minimized  (7)     Icon file:         C:\PROGRA~1\QGIS Dufour\apps\msys\msys.ico               C:\PROGRA~1\QGIS Dufour>xxmklink "C:\ProgramData\Desktop\MSYS Shell.lnk" "C:\PROGRA~1\QGIS Dufour\apps\msys\msys.bat" " " \ "Minimal SYStem" 7 "C:\PROGRA~1\QGIS Dufour\apps\msys\msys.ico"           XXMkLink     ver 1.00    (c)2005 Copyright  Pixelab, Inc.          The shortut created as follows          Shortcut path:     C:\ProgramData\Desktop\MSYS Shell.lnk     Target object:     C:\PROGRA~1\QGIS Dufour\apps\msys\msys.bat     Arguments;               Working Directory: \     Description:       Minimal SYStem     Display Mode:      Minimized  (7)     Icon file:         C:\PROGRA~1\QGIS Dufour\apps\msys\msys.ico               C:\PROGRA~1\QGIS Dufour>textreplace -std -t apps/msys/etc/fstab           C:\PROGRA~1\QGIS Dufour>dllupdate -oite -copy -reboot "C:\PROGRA~1\QGIS Dufour\bin\libeay32.dll"      C:\PROGRA~1\QGIS Dufour\bin\libeay32.dll: Version=1.0.1.3     C:\Windows\system32\libeay32.dll: Version=0.0.0.0     Target does not exist, no action.          C:\PROGRA~1\QGIS Dufour>dllupdate -oite -copy -reboot "C:\PROGRA~1\QGIS Dufour\bin\ssleay32.dll"      C:\PROGRA~1\QGIS Dufour\bin\ssleay32.dll: Version=1.0.1.3     C:\Windows\system32\ssleay32.dll: Version=0.0.0.0     Target does not exist, no action.          C:\PROGRA~1\QGIS Dufour>textreplace -std -t bin\pyuic4.bat           C:\PROGRA~1\QGIS Dufour>set OSGEO4W_ROOT=C:\\PROGRA~1\\QGIS Dufour           C:\PROGRA~1\QGIS Dufour>textreplace -std -t apps\Python27\Lib\site-packages\PyQt4\pyqtconfig.py           C:\PROGRA~1\QGIS Dufour>call C:\PROGRA~1\QGIS Dufour\bin\make-bat-for-py.bat      'C:\PROGRA~1\QGIS' is not recognized as an internal or external command,     operable program or batch file.          C:\PROGRA~1\QGIS Dufour>call "C:\PROGRA~1\QGIS Dufour"\bin\o4w_env.bat           '"C:\PROGRA~1\QGIS"' is not recognized as an internal or external command,     operable program or batch file.     OSGEO4W home is C:\PROGRA~1\QGIS Dufour     Synchronizing CRS database with GDAL/PROJ definitions.     QFSFileEngine::open: No file name specified     QFSFileEngine::open: No file name specified     QFSFileEngine::open: No file name specified     QFSFileEngine::open: No file name specified     QFSFileEngine::open: No file name specified     Loading epsg.wkt     0 WKTs loaded     CRS update (inserted:0 updated:0 deleted:3631 errors:0)     No CRS updates were necessary.          C:\PROGRA~1\QGIS Dufour>textreplace -std -t bin\qgis.bat           C:\PROGRA~1\QGIS Dufour>textreplace -std -t bin\qgis-browser.bat           C:\PROGRA~1\QGIS Dufour>REM get short path without blanks           C:\PROGRA~1\QGIS Dufour>for %i in ("C:\PROGRA~1\QGIS Dufour") do set O4W_ROOT=%~fsi           C:\PROGRA~1\QGIS Dufour>set O4W_ROOT=C:\PROGRA~1\QGIS Dufour           C:\PROGRA~1\QGIS Dufour>if not 1 == 0 mkdir "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour"      A subdirectory or file C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour already exists.          C:\PROGRA~1\QGIS Dufour>if not 1 == 0 nircmd shortcut "C:\PROGRA~1\QGIS Dufour\bin\nircmd.exe" "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour" "QGIS Desktop 2.0.1" "exec hide C:\PROGRA~1\QGIS Dufour\bin\qgis.bat" "C:\PROGRA~1\QGIS Dufour\apps\qgis\icons\QGIS.ico"           C:\PROGRA~1\QGIS Dufour>if not 1 == 0 nircmd shortcut "C:\PROGRA~1\QGIS Dufour\bin\nircmd.exe" "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour" "QGIS Browser 2.0.1" "exec hide C:\PROGRA~1\QGIS Dufour\bin\qgis-browser.bat" "C:\PROGRA~1\QGIS Dufour\apps\qgis\icons\QGIS.ico"           C:\PROGRA~1\QGIS Dufour>if not 1 == 0 nircmd shortcut "C:\PROGRA~1\QGIS Dufour\bin\nircmd.exe" "~$folder.desktop$" "QGIS Desktop 2.0.1" "exec hide C:\PROGRA~1\QGIS Dufour\bin\qgis.bat" "C:\PROGRA~1\QGIS Dufour\apps\qgis\icons\QGIS.ico"           C:\PROGRA~1\QGIS Dufour>if not 1 == 0 nircmd shortcut "C:\PROGRA~1\QGIS Dufour\bin\nircmd.exe" "~$folder.desktop$" "QGIS Browser 2.0.1" "exec hide C:\PROGRA~1\QGIS Dufour\bin\qgis-browser.bat" "C:\PROGRA~1\QGIS Dufour\apps\qgis\icons\QGIS.ico"           C:\PROGRA~1\QGIS Dufour>set OSGEO4W_ROOT=C:\\PROGRA~1\\QGIS Dufour           C:\PROGRA~1\QGIS Dufour>textreplace -std -t "C:\PROGRA~1\QGIS Dufour\apps\qgis\bin\qgis.reg"           C:\PROGRA~1\QGIS Dufour>"C:\Windows\regedit" /s "C:\PROGRA~1\QGIS Dufour\apps\qgis\bin\qgis.reg"           C:\PROGRA~1\QGIS Dufour>set OSGEO4W_ROOT=C:\\PROGRA~1\\QGIS Dufour           C:\PROGRA~1\QGIS Dufour>textreplace -std -t "C:\\PROGRA~1\\QGIS Dufour"\bin\qt.conf           C:\PROGRA~1\QGIS Dufour>xxmklink "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\Setup.lnk" "C:\PROGRA~1\QGIS Dufour\bin\setup.bat -R C:\PROGRA~1\QGIS Dufour" "" "" "Installer" 1 "C:\PROGRA~1\QGIS Dufour\OSGeo4W.ico"           XXMkLink     ver 1.00    (c)2005 Copyright  Pixelab, Inc.          The shortut created as follows          Shortcut path:     C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\Setup.lnk     Target object:     C:\PROGRA~1\QGIS Dufour\bin\setup.bat -R C:\PROGRA~1\QGIS Dufour     Arguments;              Working Directory:      Description:       Installer     Display Mode:      Normal Window (1)     Icon file:         C:\PROGRA~1\QGIS Dufour\OSGeo4W.ico               C:\PROGRA~1\QGIS Dufour>textreplace -std -t bin/setup.bat           C:\PROGRA~1\QGIS Dufour>textreplace -std -t bin/setup-test.bat           C:\PROGRA~1\QGIS Dufour>mkdir "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour"      A subdirectory or file C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour already exists.          C:\PROGRA~1\QGIS Dufour>xxmklink "C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\OSGeo4W.lnk" "C:\PROGRA~1\QGIS Dufour\OSGeo4W.bat" " " \ "OSGeo for Windows command shell" 1 "C:\PROGRA~1\QGIS Dufour\OSGeo4W.ico"           XXMkLink     ver 1.00    (c)2005 Copyright  Pixelab, Inc.          The shortut created as follows          Shortcut path:     C:\ProgramData\Microsoft\Windows\Start Menu\Programs\QGIS Dufour\OSGeo4W.lnk     Target object:     C:\PROGRA~1\QGIS Dufour\OSGeo4W.bat     Arguments;               Working Directory: \     Description:       OSGeo for Windows command shell     Display Mode:      Normal Window (1)     Icon file:         C:\PROGRA~1\QGIS Dufour\OSGeo4W.ico               C:\PROGRA~1\QGIS Dufour>xxmklink "C:\ProgramData\Desktop\OSGeo4W.lnk" "C:\PROGRA~1\QGIS Dufour\OSGeo4W.bat" " " \ "OSGeo for Windows command shell" 1 "C:\PROGRA~1\QGIS Dufour\OSGeo4W.ico"           XXMkLink     ver 1.00    (c)2005 Copyright  Pixelab, Inc.          The shortut created as follows          Shortcut path:     C:\ProgramData\Desktop\OSGeo4W.lnk     Target object:     C:\PROGRA~1\QGIS Dufour\OSGeo4W.bat     Arguments;               Working Directory: \     Description:       OSGeo for Windows command shell     Display Mode:      Normal Window (1)     Icon file:         C:\PROGRA~1\QGIS Dufour\OSGeo4W.ico               C:\PROGRA~1\QGIS Dufour>set OSGEO4W_ROOT=C:\\PROGRA~1\\QGIS Dufour           C:\PROGRA~1\QGIS Dufour>textreplace -std -t apps\Python27\Lib\site-packages\sipconfig.py      `      So what do you think now? * * *


title: how manage mysql spatial polygons in parent child style

text: I want to store location-polygon data in a parent-child relation. I am using spatial data in MySQL and Google map for polygon. Eg:   1. World      * 1.1 US     * 1.2 China     * 1.3 Japan     * 1.4 South Africa        * 1.4.1 Pretoria       * 1.4.2 Cape Town     * 1.5 India        * 1.5.1 New Delhi       * 1.5.2 Chennai Table in mysql:               CREATE TABLE IF NOT EXISTS `geom` (       `id` int(11) NOT NULL AUTO_INCREMENT,       `g` geometry DEFAULT NULL,       `location` varchar(100) NOT NULL,       PRIMARY KEY (`id`)      My questions are:   1. How can I identify India is the parent of New Delhi?   2. How to get children of India? Any spatial query?


title: Wrong sorting of elevation classes QGIS

text: I'm a totally beginner with QGIS and maybe I have a very easy problem, but I can't help myself. I got an elevation shapefile with classes of elevations like that: 1-100m, 101-200m, ... When I try to classify this dates to display them with a continuous shift of colours, QGIS always sorts the classes like that: 1-100m, 1001-1101m, 101-200m, 1101-1200m, ... So the colours of the elevation classes are mixed and don't display the continuous shift. Is there an easy way to sort the classes in the attribut table in the right way? Thanks a lot!


title: Can I use Point as the geom type in pgrouting?

text: in my table i have the field geom(POINT, srid) and its populated, this table also have the source an target field witch is blank(just like pgrouting guide). I tried the function pgr_createTopology() and it returns ok, created the vertex_pgr table but it is empty. My question is if the geom field must be a geom(LINESTRING)? My problem is: i have some points and need to trace routes beetwen them. Do i have to make line string between my points? i appreciate any response. Thanks


title: Automated merge of polygons up to a limit

text: I have a polygon feature class with a field, let's say it's population. I need to find an automated way to merge these polygons up to a given population number. E.g. I have polygons with population from 0 to 211. And I need to create districts with about 400 habitants. (Between 380 and 420.) ![enter image description here](http://i.stack.imgur.com/i0uJc.png) I work with ArcGIS 10.


title: Programming Geoserver 2.0.2 to add a raster data store and layer without the UI

text: I have a directory of raster imagery that will be updated continually. From this imagery, I am making Image Pyramids using the Geotool's PyramidBuilder utility. I need to setup a cron job to automatically add raster datastores and layers to Geoserver without using the UI. After looking at the REST section of the Geoserver manual I was able to add my workspace, "testWS", but trying the create an ImagePyramid datastore did not work. Since I have access to the datastore, I expanded on the shapefile example and tried to               curl -u admin:geoserver -XPUT -H 'Content-type: text/plain' \         -d '/opt/geoserver_data_dir/2.0.2/data/test_pyramid.pyr' \         "http://localhost:8080/geoserver/rest/workspaces/testWS/datastores/test_pyramid.pyr/external.imagepyramid?configure=all"      Where test_pyramid.pyr is the base of my ImagePyramid at this location. This gave me an error of "No such datastore: test_pyramid" Is there a better way to add a new datastore and layer to Geoserver without manually adding each one via the UI? I need help in crafting the proper REST statement that will add an extisting ImagePyramid as a datastore and layer. Is there some Java code that could do this? I looked at the Geoserver python extensions but they did not have this either.