title: Offline Phonegap/Openlayers App

text: I have a mobile app that uses Phonegap with OpenLayers and Geoserver. Now, I need to make the app available completely offline. The amount of data points can get very large depending on what client is loaded into the app, and each point has related properties. So, I need to be able to pull the data down to the phone when it is online and then use the local data when it is offline. I like the speed of GeoJSON and would like to use it if possible. I have built a GeoJSON file that loads on the map like I want it to. Now, I have to save edits/moves to the GeoJSON file and I'm stuck. Using a save strategy isn't working and from what I read it looks like I have to have a server script to make that work. How do I save the edits back to my GeoJSON file? Or am I going about this all wrong and need to use something else for offline data with OpenLayers?


title: Find Perpendicular Distance and Minimum Perpendicular Vector between a Point and a Line using GeoTools and JTS

text: I have a line formed by 2 lat/lon pairs, and the lat/lon of a point. I would like to find out the perpendicular distance between the line and the point on the Earth surface (can assume Earth as great sphere), and the minimum perpendicular vector (i.e. the projected "cross point" on the line). I am trying to use Geotools 8.0 and JTS for this. Below captured my testing code:                   //Coordinates in lon, lat         Coordinate linePt1 = new Coordinate(-5.71472, 50.06639);         Coordinate linePt2 = new Coordinate(-3.07000, 58.64389);              //Multiply all longitudes by the cosine of latitude.         //http://gis.stackexchange.com/a/29713/10772         linePt1.x = linePt1.x * Math.cos(linePt1.y);         linePt2.x = linePt2.x * Math.cos(linePt2.y);              LineString line = createLine(new Coordinate[]{linePt1, linePt2});              Coordinate pt1 = new Coordinate(-6, 54);         pt1.x = pt1.x * Math.cos(pt1.y);         Point point = createPoint(pt1.x, pt1.y);              double distanceOp = DistanceOp.distance(line, point);         System.out.println("Distance = " + distanceOp);              //Find the minimum perpendicular vector using "closestPoints()"         for (Coordinate c : DistanceOp.closestPoints(line, point)) {             System.out.println("=== " + c);             //verify if the point is on the line             System.out.println(CGAlgorithms.isOnLine(c, new Coordinate[]{linePt1, linePt2}));         }      the createPoint() and createLine() methods:               public static LineString createLine(Coordinate[] coordinates){         GeometryFactory factory = new GeometryFactory(new PrecisionModel(                 PrecisionModel.FLOATING), WGS84_SRID);         LineString line = (LineString) factory.createLineString(coordinates);         return line;     }          public static Point createPoint(double longitude, double latitude) {         if (longitude < -180 || longitude > 180) {             throw new IllegalArgumentException(                     "Longitude should be between -180 and 180");         }         if (latitude < -90 || latitude > 90) {             throw new IllegalArgumentException(                     "latitude should be between -90 and 90");         }         GeometryFactory factory = new GeometryFactory(new PrecisionModel(                 PrecisionModel.FLOATING), WGS84_SRID);         Point point = (Point) factory.createPoint(new Coordinate(longitude,                 latitude));         return point;     }      However, the result of "isOnLine()" return false. I was wondering if there is anything wrong. Is there something wrong with my verification method, or actually the way I used to find out the perpendicular distance and the "cross point" on the Earth surface is not correct? Please kindly advice


title: Making an access plan with openstreetmap

text: What tool would be best to create an access plan (with directions written on roads, local zooms, etc.) from OpenStreetMap (preferably using the MapQuest Open map)? ![example of access plan](http://i.stack.imgur.com/bpbeQ.gif)


title: How to create CGA rules in esri cityengine

text: I'm new in using esri cityengine 2011 . i red it's tutorials . it was very interesting and helpful for but there is no references or tutorials for how to create cga rules . how can i find tutorial or help about this subject?


title: Using String Functions (LEFT, RIGHT, MID) in Select By Attributes against a File Geodatabase using ArcMap 10

text: I finally (albeit reluctantly) made the switch from a personal geodatabase (pGDB) to a file geodatabse (fGDB) not too long ago. The main reason it took so long to switch has to do with how Select by Attributes deals with string functions. For example, I used to be able to easily select all records where the first character of a particular field starts with "A".               left([Some_Field],1 = "A")      This works without a hitch when querying against a pGDB. However, try to do the same thing against a fGDB.               left("Some_Field",1 = "A")      This returns the following error: There was a problem selecting An invalid SQL statement was used. I realize I can use the field calculator to calculate a new field equal to the first character of some other field. But who wants to do that every time? Is there a way to use string functions inside Select By Attributes?


title: Rhombus shaped pixels when viewing an image referenced in WGS84

text: Following the Question What is the projection of WGS84 in GIS software?, in the comments section of the answer http://gis.stackexchange.com/a/73578/5256, Jake very well explains that > Because the coordinates -180°,-90° in WGS84 are in angular units, and they > need to be mapped to cartesian coordinates. In this case, that's done using > the trivial equirectangular projection, which simply exchanges the angular > units with linear units, but that process is still a projection. I think that a related visual example of an image referenced in Geographic Coordinates is illustrated below. Perhaps my guess is completely wrong. I would like to understand what causes the "skewing" of pixels, which happens in both OpenEV as well as in ArcGIS (using data from the same satellite acquisition, packed in an NTF file). So, the question is **What is this skewing effect?** ## A visual example A small fragment from a WorldView-2 image, purely in Geographic Coordinates, for which gdalinfo reports the following               Coordinate System is `'                                                                                GCP Projection =                                                                                GEOGCS["WGS 84",                                                                                DATUM["WGS_1984",                                                                                SPHEROID["WGS 84",6378137,298.257223563,                                                                                AUTHORITY["EPSG","7030"]],                                                                                TOWGS84[0,0,0,0,0,0,0],                                                                                AUTHORITY["EPSG","6326"]],                                                                                PRIMEM["Greenwich",0,                                                                                AUTHORITY["EPSG","8901"]],                                                                                UNIT["degree",0.0174532925199433,                                                                                AUTHORITY["EPSG","9108"]],         AUTHORITY["EPSG","4326"]]     GCP[  0]: Id=UpperLeft, Info=               (0.5,0.5) -> (-16.8511111111111,13.0025,0)     GCP[  1]: Id=UpperRight, Info=               (9215.5,0.5) -> (-16.6588888888889,13.03,0)     GCP[  2]: Id=LowerRight, Info=               (9215.5,7167.5) -> (-16.6588888888889,12.8911111111111,0)     GCP[  3]: Id=LowerLeft, Info=               (0.5,7167.5) -> (-16.8511111111111,12.8641666666667,0)      when displayed in OpenEV as _Non-Georeferenced_ (squared pixels, term derived from OpenEV's related Viewer icon/button) ![small fragment from a WorldView-2 image, Non-Georeferenced viewing](http://i.stack.imgur.com/gxmUW.jpg) and Georeferenced (rhombus-shaped pixels!) ![small fragment from a WorldView-2 image, Georeferenced viewing](http://i.stack.imgur.com/ub7az.jpg). As explained in OpenEV's manual, for the _Georeferenced Display_ icon/button: **_Non-Georeferenced_** means that _the view is being displayed in raw pixel/line coordinates for the source raster_ and **_Georeferenced_** means that _the view is trying to display in georeferenced units (if any are available)_. The same skewing effect appears also in ArcGIS when viewing another fragment from the same acqusition (credits for the screenshot to a friend): ![skewing effect in ArcGIS](http://i.stack.imgur.com/VBQ2y.jpg)


title: Error .ecw is not a supported raster data source

text: I am a user from Spain and a I have recently installed QGIS 2.0 in my computer which has Ubuntu 13.04 and when I try to use an .ecw file it says: Error, .ecw is not a supported raster data source. I am new in Linux and in QGIS, so please could you help me and explain what could I do step by step. Thank you very much.


title: How to drag label using DragFeature or ModifyFeature?

text: I am new to openlayer. I have seen the examples in which DragFeature is able to move my points, polygon and lines. But when it comes to test the Label moving in IE8, then the dragging functionality doesn't work. though I am able to drag the label in Firefox. Please provide me help here.


title: How to read BMP file location of image used for picture symbol in LYR file?

text: I have a number of *.lyr files which are pointing to various .bmp files which are used as marker symbols in ArcGIS for Desktop 10.1 I was wondering if there is a way to extract the file paths of the *.bmp files each *.lyr file uses for its picture symbol. Is it possible to do this using Python?


title: Categorized Style : showing all possible values, not just the ones in use

text: I am new to QGIS and experimenting before rolling this out to my colleagues. I have a shapefile with some codes, and i have a csv with those codes and other information that i need. I have joined them, and that works fine - BUT is there a way to stylize the shapefile so that only records present in the shapefile show up in the layer manager and legend? At present, my shapefile is nicely coloured by "description" but the layer props shows every single code/description in the CSV, not just the ones in use in the shapefile. In ESRI (for example) it would only show records with non- null values - I only have 4 or 5 codes in the shapefile, and, i mean, some of the values listed i know are not even polygons! But they are displayed, because the csv file has them.


title: calculate average commute time within 15 miles of city center

text: I am trying to figure out how to calculate the average commute time to work within 15 miles of the city center (and then map it in a form of continuous symbols (red to green gradient) with a legend, scale and everything. That would be called “average commute times to work). Could you help me for the first part of the problem? Thank you very much


title: keep the same field headers after joining a table to a shape file in QGIS

text: I noticed that after joining a table to a shape file the headers of each fields in the table used for the join are changed by QGIS and they are attributed to a different name. For instance, i had a shape file with vegetation_codes and a table with vegetation_codes and veg_definition (these were the headers). After joining the table to the shape file and also saving the shape file and creating a new one, the vegetation_classes header was changed into vegetati_1. Is there a way to keep the original headers? Notice: i am using QGIS dufour on windows 7 and teh table that i join to the shape file is csv (tried with dbf as well and the problem is the same)


title: QGIS 1.8 not loading after install

text: I have just installed QGIS 1.8 on Windows 7. Once the install completed and I tried to open QGIS (via the desktop icon) nothing happened. I had this issue when installing QGIS 1.7. I was able to get 1.7 to open by using the qgis.bat file, but it is not in the the bin folder in the QGIS 1.8 version. I began to have this issue after installing OpenGeoSuite 4W. I thought I uninstalled everything I needed. I am worried I may need to re-install Windows. Does anyone have any ideas as to what I can do to get QGIS 1.8 to load and open. Thanks...


title: Does GPS Navigator send data to Satellite?

text: Does GPS Navigator send data to satellite? Is this true about regular devises available in the market?


title: ArcGIS 10.2 Release Date

text: Anyone have an idea on an ArcGIS 10.2 release date? I wans't even aware of it until stumbling upon this article: UC Agenda Sneak Peek: Part I


title: Unable to explode multipart polygon

text: I have a taxlot that is giving me troubles in a Union. When I try to run the union I get the tried and true ''Geometry is not M-Aware'' message. (there are no Polygon-M features). Based on readings on this site and the ESRI forums, I'm now trying to explode my taxlots prior to the union. However I seem to have come across a mystical indestructible shape. I've tried check and repair geometry with no results. I've tried the Multipart-to-Singlepart tool and receive this message:               Cannot process feature with OID ...     Failed to execute      So I tried the explode function in the Advanced Editing Toolbar and go this lovely message: ![explode fail](http://i.stack.imgur.com/hEl0v.png) I calculated (!SHAPE!.partCount) the number of parts into a field and came up with 22. Here is what the shape looks like:![enter image description here](http://i.stack.imgur.com/Hg2Lu.png) Help?


title: How to convert from dgn to shp or dwg files?

text: I am a frequent recipient of microstation dgn files and don't like it because I don't have software to convert it to autocad. I do have ARCVIEW, but I can't edit the files. I have tried other open source programs to convert, but they all are lacking. I'm wondering if anyone has used QGIS to convert from dgn to shp or dwg file types? I'm skeptical because I noticed that QGIS supports dgn file types, but **_doesn't_** support dwg file formats. This is odd... Anyways, please post up, I'm probably going to install QGIS this week on a production machine! Thanks


title: What is the best way to handle geometries between Python functions?

text: What is the best way to handle geometries between different Python functions? For example: I have a function that exports a line from grass.script. The line will be a parameter of another function where it will be passed in. Is it best to export them to WKT, or safe them in a GPX, or GeoJSON? The formats OGR offers are endless. But which one is most recommended?


title: GeoTools Oblique Mercator Example

text: I would like to use the Oblique Mercator projection in GeoTools. I'm currently using a WGS84 Mercator projection like this:                       final CoordinateReferenceSystem wgs84CoordSystem = CRS                     .decode(LATLON_EPSG_CODE);             final CoordinateReferenceSystem worldMercatorCoordSystem = CRS                     .decode(MERCATOR_EPSG_CODE);             final MathTransform latLonToMercator = CRS.findMathTransform(wgs84CoordSystem,                     worldMercatorCoordSystem);             final MathTransform mercatorToLatLon = CRS.findMathTransform(worldMercatorCoordSystem,                     wgs84CoordSystem);      The Oblique Mercator EPSG code doesn't appear to be supported (and I want to be able to specify the latitude of centre and azimuth anyway). Could anyone give me an example of how to create a WGS84 Oblique Mercator projection and MathTransform to allow me to convert to/from Lat/Lon please?


title: Spherical Mercator projection invented by Google and not conformal?

text: I had some certitudes on the Google Spherical Mercator also known under EPSG code 900913 and 3857 before looking at this presentation, I was sure that this projection was in the Mercator family (so conformal). With this input, I'm now not sure about conformal status of Spherical Mercator. I didn't dive enough in mathematical "stuff" here to be able to qualify validity of the doc. I also want to confirm / infirm this exact projection was never there before Google Maps landing. Any opinions on those 2 questions?


title: Kernel Density output labeling

text: I have been looking for a way to try and make the output of the **Kernel Density** function in ArcGIS for Desktop make a little more sense to our administration. I ran the function on a sample dataset with 500 points, and the 'density' where only 1 point was within the search radius was given a value of 6.2342 at the peak. I then divided the entire raster by this value. The scale is now at least in relation to the peak value found at one point. Is this methodology sound?   Has anyone found better ways to label the scale that go beyond 'more dense, less dense'? EDIT: Just wanted to reference How to interpret GRASS v.kernel results? as the method for checking the values of the output is beneficial. This helps to clarify how these values are computed and provides a way to validate results which helped me to understand the assigned values.


title: Can't save result of "Extract Values to Points (Spatial Analyst)" inside the working feature dataset!

text: I have a file geodatabase in which I have a feature dataset. I have some feature classes inside the feature dataset. I wonder why when I attempt to save the result of Extract Values to Points (Spatial Analyst) inside my working feature dataset, the result is genereated inside the file Geodatabase but not under my feature dataset. No matter what I do I can't make it save the result inside my feature dataset. I can easily save the results of for instance Extract Values to Points (Spatial Analyst) inside my working feature dataset. But "Extract Values to Points (Spatial Analyst)" seems to be different in that regard. Any idea why this is happening? Edits: I am using Arcmap 10.0 calling this tool from Arc toolbox. All the layers are in the same coordinate system. Here is a snapshot before I run the tool: ![enter image description here](http://i.stack.imgur.com/ssRoY.jpg) But as can be seen here the result is generated outside the feature dataset: ![enter image description here](http://i.stack.imgur.com/Frgba.jpg)


title: I just want to make a simple loop in PostgreSQL

text: In the PostGIS SQL Editor, I just want to make a loop to make this function work This query is to give the UTM Zone number to one column of Global DTM raster table               For Loopid = 0 to 1000, then execute the following statments:          UPDATE public.globaldtm      SET "UTM" = loopid      WHERE rid IN (SELECT rid FROM globaldtm WHERE ST_Intersects(rast,(select geom from utm where gid =loopid)));      ![enter image description here](http://i.stack.imgur.com/mqyEJ.png) How to write a simple code inside this SQL editor to make the loop work. Thanks.


title: PostGIS - Interpolate point between contours

text: Using PostGIS, what is the best method to interpolate the z value for points (the orange dot in the image below) located between each contour layer (blue lines and red lines in the image below)? So for the blue contours it's going to be ~ 917 and for the red contours it'll be ~ 819. ![example-map](http://i.stack.imgur.com/UUmeS.jpg) The blue contours are a 2D Multipolygon and the red contours are a 2D Multilinestring (so for both, the contour value is an attribute field. If necessary I could update the geometries to have a z value). Thanks!


title: Geometry Operands - BBOX vs Within - What's the difference?

text: What is the difference between the "BBOX" and "Within" spatial operators. Can anyone please shed some light?


title: How to serialize IStyleGallery Object?

text: > **Possible Duplicate:**   >  Can I store an ArcObject inside a BLOB or XML field? I am saving shape file in database and now i want the symbology selected for the feature layer also get stored in database. I am facing an issue I cannot serialize IStyleGallery, and if its not serialized than how can i store that in database. Anyone got any good solution?


title: Pan sharpen Landsat 8 using Grass 7

text: I have followed this great tutorial on pan sharpening using Grass GIS 7 (I'm using beta 3). Everything in the tutorial seemed to have worked, but my final output is not as expected (pictured below). ![PanSharp RGB output](http://i.stack.imgur.com/v5hx0.png) I tried with various sharpening methods, colour balancing stages and with and without atmospheric correction. I also tried with different images. Any advice on where I'm going wrong? Here's the output from r.univar: ## Of the non-null cells: n: 41623319 minimum: 0.089549 maximum: 1.17079 range: 1.08124 mean: 0.161591 mean of absolute values: 0.161591 standard deviation: 0.110987 variance: 0.0123181 variation coefficient: 68.6836 % sum: 6725971.35730981


title: Recommendations for creating contour map from depth points

text: I am trying to put together a safe routes map for a local lake. I have an old paper map that has a grid of various depths as points. I would like to draw contour lines through these points and create a topographic map of the lake depth. From here I would like to apply a graduated shading pattern for some specific depths to show safe routes etc. I've mostly used MapInfo, GlobalMapper and such tools, and mainly for generating heat maps on targeted demographic data. I've never worked with contours or other topographic data before. Any recommendations on tools/plugins or online guides that might help me on my project? Thanks!


title: How can I use GDAL to batch define a projection?

text: I have a fodler with 1500 ECW rasters, no spatial reference information defined. Can I use GDAL to batch define a projection for all 1500 in the folder?


title: Clip raster layer with *raster* mask layer in QGIS

text: Given two input layers   1\. The clip layer - an arbitrary layer (though in my case it's a float 1 band geotiff)   2\. The mask layer - a paletted _raster_ in which every pixel is either NO DATA or a single "color" value, same for every pixel. I want to create an output raster whose overall extent matches the extent of the mask layer and each pixel is either NO DATA if the corresponding pixel of the mask layer is NO DATA or the value of the corresponding pixel from the clip layer if the corresponding pixel of the mask layer is not NO DATA. I know I can mask the clip layer with a polygon, but the mask layer is a raster. Polygonizing the raster layer is possible in principle, but would result in a polygon with a lot of very short edges. Is there a direct way to do this in QGIS? I suppose there's probably a path that involves creating a new layer and using the raster calculator, but I'm hoping there's something more direct. Further complicating the matter is that the two rasters might not have the same resolution and origin. I believe this is equivalent to Extract by Mask in ArcMap. A gdal command solution is equally acceptable to QGIS. A solution in R is OK, but not preferred.


title: Data source for determining postal codes given location?

text: What is the best data source for determining international postal codes given a latitude & longitude?


title: d.out.file is not working

text: I'm trying to use 'd.out.file' to save the displayed map to a PNG file. It failed and produced this message: "Command 'd.out.file' not yet implemented in the GUI. Try adding it as a command layer instead." How can I resolve this? It also seems that I cannot run other d.* commands


title: qgis - slope from DEM - errors due to projection and oversimplified 'z factor'

text: I'm trying to generate a slope raster from a DEM raster. It seems there is a pretty significant error due to the 'z-factor' option. I'm fairly new to this, but some other posts and responses indicate a general workflow: 1) obtain a DEM. The ones I got (NED) are unprojected: horizontal coords are lat-lon, and vertical coords are meters. Two possible branches here, A or B: A2) project the DEM to UTM so that horizontal and vertical coords are in the same units A3) then do slope analysis on the UTM-projected DEM, using z-factor = 1.0 OR B2) do slope analysis on the unprojected DEM, using z-factor = 111120 (which has been suggested in some posts as a number of meters per degree - see error discussion below) B3) then project the raster to UTM Problem with option A: the projection process introduces gridded artifacts, and these get amplified in the slope raster: Original DEM (NAD83, unprojected; zoomed in to one peak; min/max set to show high contrast here): ![Original DEM \(unprojected\)](http://i.stack.imgur.com/Z0s4A.png) DEM projected to NAD83 UTM zone 10, note gridded artifacts: ![enter image description here](http://i.stack.imgur.com/EyGrn.png) slope analysis output from NAD83-UTMz10-projected DEM, note the same artifacts: ![enter image description here](http://i.stack.imgur.com/r7sGk.png) Problem with option B: 111120 seems vastly oversimplified and produces significant errors, see error table below - though, gridded artifacts are greatly reduced since they don't exist in the DEM that you run slope analysis on. First, number of meters per degree of latitude is not the same at all latitudes, due to the earth not being a perfect sphere. Second, and more obvious, number of meters per degree of longitude is a function of latitude, varying from ~111km at the equator to 0 at the poles. For our area (Nevada County, California) here are some numbers from an online calculator, referenced in other places as a trusted source (but feel free to verify these): Nevada County S extent ~ 39.0deg (WGS84) Nevada County N extent ~ 39.5deg (WGS84) Nevada City ~ 39.25deg (WGS84) 39.0 deg: 1 deg lat = 111015.45 meters 1 deg lon = 86626.37 meters 39.25 deg: 1 deg lat = 111020.23 meters 1 deg lon = 86320.71 meters 39.5 deg: 1 deg lat = 111025.01 meters (about 9.5 meters longer than at S extent of Nevada County) 1 deg lon = 86013.39 meters (about 613 meters shorter than at S extent of Nevada County) (at equator, for reference:) 0.0 deg: 1 deg lat = 110574.27 meters 1 deg lon = 111319.46 meters So, making some straightforward assumptions about how qgis uses the 'Z Factor' argument, here are some resulting errors (note this is in Nevada County but not the same location as the above images), along with Range-Bearing Line from Terrain Navigator Pro which is treated as gospel in this case; also note that this is just for one sample set of from/to; error would vary by compass direction of the from-to line also: ![enter image description here](http://i.stack.imgur.com/2WrZu.png) (z-factor is entered by hand for each row; all the other columns are calculated from it) So, any ideas? Am I calculating something wrong, or is there a better workflow, or is there just no really good way to get around the error?


title: What is the maximum attachment file size in Esri file and ArcSDE geodatabases?

text: I am trying to figure out the maximum size of attachment in a file and ArcSDE geodatabases while working on a project specification details. I am not able to find these figures documented by Esri anywhere. Running own tests, it seems as a file geodatabase has a limit of ca 600MB per file that is being attached. When trying to attach everything that is larger than 600MB, I am getting an error (where `Points3D` is the name of the feature class I am trying to attach a file to). ![enter image description here](http://i.stack.imgur.com/QGFIz.png) Just to see whether it is a file geodatabase only limit, I have run some tests with a SQL Server enterprise geodatabase and here I am getting stuck at somewhere between 700 and 800MB and either get the same error message (except that instead of the feature class name I have my SQL Server database filegroup name) or have ArcMap unexpectadly crashed and closed. **What is the official/documented maximum file size limit that one can upload to an Esri file geodatabase?** I also welcome comments if someone has succeeded uploading anything larger than 1GB to a geodatabase as an attachment.


title: What does shading indicate in home range plot?

text: I've simulated a home range in R using the clusthr function from the adehabitatHR package. The home range is made from 50 x,y coordinates. The output is attached. ![enter image description here](http://i.stack.imgur.com/AQYIE.png) I'm trying to get clear on what the shading indicates. The shading is clearly related to the size of the polygon. Smaller polygons = darker shading, larger polygons = paler shading. But how exactly is a polygon assigned a shade? Does the shading represent a percentage of the area covered?


title: PostGIS stored procedures

text: I am wondering if there are open source applications that make extensive uses of the special functions (the ST_* stored procedures) provided by PostGIS? I found many apps that support storing map data in PostGIS, but they seem to just fetch and store data from PostGIS and then operate on them in application code. I would like to learn about how to use those special functions as part of the queries by looking at sample applications. Thanks!


title: OpenLayers code didn't response popup featureinfo

text: Hey I'm trying to popup the featureinfo of WMS layer using OpenLayers code. But it does'nt work. Instead the cursor moves in the state of continous moving when I click on map. Here is my code:               <title>Untitled 1</title>             <script src="http://localhost:8080/geoserver/openlayers/OpenLayers.js" type="text/javascript">         </script>         <script type="text/javascript">             var map;              var wms;              var wms1;                  function init(){                      alert("inside function");                 var bounds = new OpenLayers.Bounds( 143.853, -43.648,                     148.479, -39.574);                       var options = {             controls: [],             maxExtent: bounds,             maxResolution: 0.0180703125,             projection: "EPSG:4326",             units: 'm'             };                           map = new OpenLayers.Map('map', options);                 wms = new OpenLayers.Layer.WMS("Peoples Colony",                                     "http://localhost:8080/geoserver/wms",                                    {                                         LAYERS: "topp:tasmania_state_boundaries",                                        STYLES: '',                                        format: "image/jpeg",                                        transparent: "true",                                        tiled: 'true'                                         },                                    {                              isBaseLayer: true,                         yx : {'EPSG:4326' : true}                                         }                                    );                                    wms1 = new OpenLayers.Layer.WMS("tasmania_roads",                                     "http://localhost:8080/geoserver/wms",                                    {                                         LAYERS: "topp:tasmania_roads",                                        STYLES: '',                                        format: "image/jpeg",                                        transparent: "true",                                        tiled: 'true'                                         },                                    {                              isBaseLayer: false,                         yx : {'EPSG:4326' : true}                                         }                                    );             map.addLayers([wms,wms1]);             map.zoomToMaxExtent();                       map.addControl(new OpenLayers.Control.Navigation());                   map.addControl(new OpenLayers.Control.Scale($('location')));                       var info = new OpenLayers.Control.WMSGetFeatureInfo({             url: "http://localhost:8080/geoserver/wms",              title: 'Identify features by clicking',             queryVisible: true,             eventListeners: {                 "getfeatureinfo": function(event) {                     map.addPopup(new OpenLayers.Popup.FramedCloud(                         "chicken",                          map.getLonLatFromPixel(event.xy),                         null,                         event.text,                         null,                         true                     ));                 }             }         });         map.addPopup(info);         info.activate();                              }               </script>         Anyone plz help me out? I think that I didn't mention the layers here. Is this the problem?


title: Is it possible to have measurements in centimeter and decimeter in 2D flat (geometry)?

text: Is it possible to have measurements in centimeter and decimeter in 2D flat (geometry) for ArcGIS?


title: Download MODIS data for Mexico

text: I would like to download MODIS npp images for Mexico but I cannot figure out which files correspond to the country. How can I associate each tile to its location without looking at the million of browsing images. I am an experienced python and R user if that helps. I am manually downloading the images from this link: http://e4ftl01.cr.usgs.gov/MOLT/MOD17A3.055/2000.01.01/


title: XMLHttpRequests clogging/pending in Leaflet.js

text: If I have a lot of active layers and eg. zooom in and out, dozens of tile xmlhttprequests are sent. This causes a problem in most browsers, clogging it so to speak, and leaving the map blank and dozens of Pending resource requests in Network Console (in eg. Chrome). Any known work-around for this issue? I'm creating a map with dozens of layers for a client, and with some activity on the map (zoom, adding lots of layers simultaneously, etc), this issue will appear. What is needed, I believe, is some internal queue in Leaflet.js, as the queue in eg. Chrome browser is buggy. Any advice much appreciated.


title: Are there any open source network analysis tools?

text: A researcher in my field* recently created an "Urban Network Toolbox", an ArcGIS extension built upon Network Analyst. It calculates a number of measures on urban networks, such as betweenness and reach. He was criticized for relying on ArcGIS, a commercial tool not available to many researchers, especially outside the US. But ArcGIS was a convenient platform, since it contains the Network Analyst extension, simplifying the work required for the analysis. **Are there open source network analysis tools that do the work of the Network Analyst extension?** Are there other methods for creating the _network datasets_ required as input for the toolbox? My understanding of the toolbox is that it uses the network dataset to calculate a distance matrix for a given set of points on the network. The toolbox does not use any of the higher level problem solving tools in Network Analyst. * _Quantitative Urban Morphology, perhaps? No one seems settled on a name for the broader field. The particular brand of research I do is with Space Syntax._


title: Find point close to points from postgis table

text: I have a postgis point table with a large list of cities. I need to find duplicates in it, but only if they are 10km close each other, which means that is the same city. For example, I have 5 points called 'Eneby', and supose that 4 of this are 10Km close each others, and the 5th is 300Km away, what will mean that is another Eneby town. The query in this case should only have to show 4 points. So, I've made this query, but is not showing the expected result:               Select p1.id, p2.id, * from points as p1, points as p2 where p1.name='Eneby'     and st_dwithin(p1.geom,p2.geom,0.1, true)     and p1.id=p2.id      The result is the same 5 points, no matter what distance (in degrees) I set. I'm trying with different distance values, but is always the same result. Does anybody know if is there any mistake in this query? I've been looking through this forum, and found some simmilar questions, but can't find exactly this problem. Thanks.


title: How to properly use PostGIS with Maps

text: First time using PostGIS in an app. Target:   \- Get Heat Map of customers from address table using leaflet.js and heatmap.js Have DB called store_dev. Installed OSM Bright from Mapbox in it. Have a table called "addresses" with street, city, state, zip, lat and lng. When user enters Address I get lat / lng from openstreetmap nomination service. Do I need to store anything else?   Nomination returns things like boundingbox and polygonpoints, any advantage on storing this?


title: Add a time dimension in ESRI mosaic raster dataset

text: I am in the process of creating a script to pull daily .asc files from a ftp server, convert them to a raster format, and push them into a mosaic raster dataset. Ultimately, the raster dataset will be updated with a new image every day, and I would like a time slider to appear when the dataset is viewed in ArcMap. I think I will want to use:               AddRastersToMosaicDataset_management(in_mosaic_dataset, raster_type, input_path,                                                {UPDATE_CELL_SIZES | NO_CELL_SIZES},                                                {UPDATE_BOUNDARY | NO_BOUNDARY},                                                {NO_OVERVIEWS | UPDATE_OVERVIEWS},                                                {maximum_pyramid_levels},                                                {maximum_cell_size},{minimum_dimension},                                                {spatial_reference},                                                {filter}, {SUBFOLDERS | NO_SUBFOLDERS},                                                {INCLUDE_DUPLICATES | EXCLUDE_DUPLICATES},                                                {NO_PYRAMIDS | BUILD_PYRAMIDS},                                                {NO_STATISTICS | CALCULATE_STATISTICS},                                                {NO_THUMBNAILS | BUILD_THUMBNAILS},                                                {operation_description},                                                {NO_FORCE_SPATIAL_REFERENCE | FORCE_SPATIAL_REFERENCE})      This seems to be used for a 'batch' of raster files, however (e.g. fiter='*.tiff'). I will only be adding one raster at a time (1/day after processing the received .asc files). It is not clear to me exactly how / where I'll add the time information. The rasters will all be overlapping, so having a time slider is rather critical... or it won't be so interesting! Does someone have suggestions on how to accomplish what I am trying to achieve?


title: Model not editable on a different PC

text: I am using ArcGIS 10.1 model builder on windows 8 to create create models for my processes. But when I open the model using Arc GIS 10.1 open model for editing with Arc GIS 10.1 running on Windows 7 it appears not to be editable. I can run the model but I cannot edit it. When I open the model in ArcGIS 10.1 running on windows 7 the model appears as on the screenshot attached below...!Does anyone have an idea what the problem might be Could? I am yet to check if I will have the same problem when I try and edit the model on another windows 8 platform... The model does not use additional tool from extensions like "spatial analyst and so on..." it only uses data management tools. Thank you in advance. ![enter image description here](http://i.stack.imgur.com/9kknf.jpg)


title: Correct coordinate systems, but trustworthy distance between points?

text: I have limited GIS experience with projections and I have some questions about trying to accurately display point data from two different sources in ArcGIS 10.1 in order to do a proximity analysis. The first data set, which already contains Latitude and Longitude coordinates, is from a state government agency and its metadata lists Horizontal Coordinate Scheme (projection): Not Applicable, Ellipsoid: GRS80, Horizontal Datum: NAD83, Horizontal Units: Feet, Distance Resolution: Not Applicable, and Geographic Coordinate Scheme: Decimal Degrees. The second data set contains addresses that must be geocoded. I used the 10.0 US Streets Geocode Service (ArcGIS Online) from the Geocoding toolbar to create a point layer based on this second dataset. Everything appears to line up perfectly well, but since the goal is an accurate proximity analysis, I hoped to confirm that the distances among points seemed plausible before doing more work. I have tried different approaches for projections, ranging from assuming the geocoded points were WGS84 and then projecting as NAD83 to trying to re- project both of the layers to a Projected Coordinate System such as a relevant State Plane NAD83 system. Again, the points look fine and the "Source" coordinate system information matches across both layers and the data frame. The problem is that, no matter what I do, the distance tool is giving me what seem to be extremely low distance values compared to an online tool that measures linear distance between points. For example, the online tool reports a linear distance of 6,214 feet between 2 reference points, while different projection approaches with ArcGIS are giving me distances of 4,418.753176 feet and 4,384.41136 feet between the same 2 reference points. Here are my questions: 1) Does the 10.0 US Streets Geocode Service automatically display points with WGS84? If so, is the correct strategy to then define the Geographic Coordinate System for the points layer as WGS84 and, after that, use the Projections and Transformations->Feature->Project tool to display in NAD83? 2) Is it a mistake to think of a website that uses Google Maps to calculate linear distance between points as some sort of verification? Are the reported distances on the website also flawed due to the "Web Mercator" projection used for Google Maps? If it is a mistake, can anyone suggest a better resource? Does an 'authoritative' resource even exist? Any feedback or advice would be greatly appreciated. Apologies if I'm incorrectly using some of the terminology or missing something very basic. Thank you!


title: Get info from WMS layer with Openlayers mobile

text: I am trying to build a webapp for iPad using OpenLayers Mobile. One of the features it should have is that you can tap on a tree displayed on a basemap and get info about this tree. These trees are stored as a WMS service in geoserver and the info should be displayed in a popup. In the standard version of OpenLayers I would use WMSGetFeatureInfo to achieve this, but this seems not to be working with the mobile version. When using OpenLayers mobile, I get an error message saying: TypeError: OpenLayers.Control.WMSGetFeatureInfo is not a constructor Can you help me fix this? Thnx in advance.


title: How to set the length of a line using openlayers

text: I'm trying edit a polygon by changing the length of some lines by code. Found this method to calculate based on the bearing. But I didn't find anything on how to calculate the bearing, before I start research my math notes and code this bearing calculation, any recomendations ? thanks


title: plot path of points

text: I have a sample of the address histories of households.               ID  ZIP1   ZIP2   ZIP3     01  XXXXX  XXXXX  XXXXX     02  XXXXX  XXXXX  XXXXX     03  XXXXX  XXXXX  XXXXX     04  XXXXX  XXXXX  XXXXX     ...      It's easy enough for me to find the centroids of the ZIP codes. What I would like to do is link those centroids together to create an image showing household `i`'s path from `ZIP1` to `ZIP2` to `ZIP3`. Does anyone know any tools to make this easy in QGIS (or even better: R)? I have the feeling there must be some tool, but I don't know enough to think of what the name would be.


title: How to manage user privilege on layers to Arcgis 10 Desktop?

text: **Problem:** I am a Computer Analyst on a town hall. Here we have the sector A and B. Both work in ArcGIS 10 Desktop. The sector A has granted to access/edit everything. The sector A has all the privileges! The sector B can see everything. But, the sector B can only edit some Layers/Dataset. I'm having problems with the Oracle user B sector in SDE. **What I tried:** I tried to grant privileges in Oracle. To the user of sector B I grant to access/edit all. I only blocked the sector B some tables (layers). That he can't change. The user of the sector B can access the ArcGIS Desktop and can see all Layers. But, the user of sector B get an error when, by the ArcGIS Desktop, he try to edit a Layer that he has privileges. I can't understand, because I grant all the tables, including the temp tables of ArcGIS, just don't grant some tables (of some layers). The error occurs when giving a "Start editing" a layer that user B has full permission. **That is the description of the error:** " Error: No editable layers ArcGIS 10 None of the layers in your map can be edited. You cannot start an edit session on this data. Solution: Make sure that the data you want to edit is in a format that can be edited in ArcGIS and that you have the proper permissions to edit it. More information: In ArcGIS, you can edit data from shape files; geodatabase feature classes; stand-alone table formats, such as INFO or .dbf; and the attribute tables of certain rasters. If your map has data in any other format, you are unable to edit it. Make sure that the data you are attempting to edit is in one of these supported formats. If the data is not from any of these workspaces, it may be possible to export it to one of these formats to be able to edit it. Many times, you encounter this error along with layer warnings that may provide more specific reasons why the layer cannot be edited. For example, if you have data from a geodatabase in the map but the feature classes are read-only, you also receive a warning message telling you that you need to check the permissions. " **Question:** Anyone have any suggestions? How do you control your users in ArcGIS 10 Desktop? What I did wrong? Can I manage in ArcGIS Server Manager users of ArcGIS Desktop? The idea is simple, the user sector B in ArcGIS Desktop can not access some layers. But I do not know do it. **Thank you!**


title: How to join nodes to network lines ( and not the other points)

text: I'm trying to do a shortest path analysis on my network. In order to prepare the network and attach polygon centroids to the roads, I used GRASS' V.NET operation=connect. the problem with this approach is that instead of connecting points to nearest line, it connect nodes to nearest line OR node which in turn results in unrealistic paths. For example, as you can see in the screen shot, nodes in the central region instead of being connected to the nearest road, are connected to the line that is three time further away. Is there a way to fix this? Thanks P.S. I'm using Grass 6.4 on Mac. ![enter image description here](http://i.stack.imgur.com/3OIRP.png)


title: Arcpy error message: 'Layer' object is not subscriptable?

text: What the code does successfully so far, is gets a list of the layers in the first data frame of the current map. What it doesn't succeed in doing is, joining the layer names to the gdb location using os.path.join(), so I can then do a clip against an existing feature class. I think I need to convert the lyr names to a string before I can use it in a Clip_analysis, how do I do that? Here is what I've got so far:               import arcpy     import arcpy.mapping     import os          #Set the current map     mxd = arcpy.mapping.MapDocument("CURRENT")          # Set the dataframe     df = arcpy.mapping.ListDataFrames(mxd)[0]          # Get list of Layers from toc     layers = arcpy.mapping.ListLayers(mxd,"",df)          #Set the output workspace     outWorkspace = r'C:\Users\jsommerville\Documents\ArcGIS\California\ProjectData.gdb'          for lyr in layers:          outFeatureClass = os.path.join(outWorkspace, lyr)          arcpy.Clip_analysis(lyr, "extent", outFeatureClass)      The error message I'm getting: 'Layer' object is not subscriptable


title: How to configure 3rd parties for the sextante plugin in qgis?

text: I'm currently exploring the enourmous potencial that the sextante plugin for qgis has to offer. The model builder will surelly be a very handy tool. The problem is I'm not beeing able to run 3rd party (like SAGA or Orfeo). I know I have to configure the paths, but i'm not sure how. I did it like this 1 but still don't get tools running. ![enter image description here](http://i.stack.imgur.com/MDcAJ.png) I'm getting this error: ![enter image description here](http://i.stack.imgur.com/G5jEB.png) ERROR with SAGA:               global name 'false' is not defined Traceback (most recent call last):        File "C:\Users\vbnunes/.qgis/python/plugins\sextante\core\GeoAlgorithm.py", line 115, in execute     self.processAlgorithm(progress)        File "C:\Users\vbnunes/.qgis/python/plugins\sextante\saga\SagaAlgorithm.py", line 153, in processAlgorithm     self.calculateResamplingExtent()        File "C:\Users\vbnunes/.qgis/python/plugins\sextante\saga\SagaAlgorithm.py", line 109, i      Any help? Thanks


title: ERDAS Imagine: negative z-coordinates for automatically set tie points

text: I'm a new ERDAS Imagine user, so it's probably a beginner problem I'm facing: I get negative elevations in my DTM after automatic tie points setting.. I am geo-referencing a stereo-pair of ASTER images to a reference image (with x,y,z, coordinates) in order to generate a DTM from them. When I set the GCPs in the point measurement tool of LPS _without_ the use of tie points, then triangulate, I can generate a quite useful DEM. However, if I use the option 'automatic tie generation' and check the coordinates after triangulation, I get, what seems, random elevations for the tie points (although x and y seem ok, both in the ASTER-images and the reference image): basically any value between -30'000m and +50'000m. Consequentially, the resulting DEM is useless (but this also confuses me somewhat, as I don't expect the 'DTM extraction' process to use the z-coordinates of the tie points, but only the ones of the GCPs). Anybody has experience with this? Any clue what I am missing? I'm glad for every hint!


title: How do I convert a text file with coordinates to database table in Python?

text: How do I convert a text file with point coordinates to database tables in the Dbase.dbf format in Python?


title: Cloud and Snow Differentiation

text: I am doing unsupervised classification of a mountaneous terrain and I have used 500 classes to classify the land use. In my result, **cloud** and the **snow** are categorized in the same class. How do I differentiate them. Can someone please help me? I want to eliminate the cloud but not the snow. Possible solution to it in either ERDAS IMAGINE or ArcGIS would be appreciated. Thanks!


title: Attribute error:'NoneType' object has no attribute

text: I tried to run a python script for getting the coordinate system of a feature/shapefile as follow               >>> import arcpy       >>> from arcpy import env       >>> env.workspace="D:\data"       >>> fc="river.shp"       >>> dsc=arcpy.Describe(fc)       >>> dsc.spatialReference.Name       >>> print dsc        i'm getting following error at line 5 > Attributeerror:'NoneType' object has no attribute Name.


title: What is primary difference between ArcGIS for Desktop and ArcGIS for Server?

text: What is the primary difference between ArcGIS for Desktop and ArcGIS for Server?


title: Implementing non-default keyboard shortcuts in Openlayers

text: Following on from this thread, I am trying to implement @Aragon's suggestion and add some keyboard shortcuts to OpenLayers via a handleKeyEvent method in `OpenLayers.Handler.Keyboard`. This ought to be a simple task that many people implement but for some reason I'm just not getting it.               var keyboardControl = new OpenLayers.Control();          var keyHandler = new OpenLayers.Handler.Keyboard(keyboardControl, {         handleKeyEvent: function(evt) {                 alert("You pressed a key!");                 switch(evt.keyCode) {                     case 87:                         doSomething([0, 1]);                     case 83:                         doSomething([0,-1]);                     case 65:                         doSomething([1,0]);                     case 68:                         doSomething([-1,0])                 }             },     });          function doSomething(vals){         alert("Something will happen here eventually...");     };          map.addControl(keyboardControl);      I get no errors but neither of my alerts get fired. What am I doing wrong? Many thanks!


title: Convert and export lat longs from text layer - QGIS

text: I have imported a `Text Delimited Layer` that was a `csv` file with three fields `id, lat, long`. The layer shows up correctly and the values are in `Feet`. I need to convert the `Feet` to `Decimal Degrees` and then export the file. None of the menu options seem to point towards this goal. Can this be done? I have been searching for solutions on the Internet without luck so far. Edit: I tried right clicking on the layer and `Save As` to csv, but it has the same values even after changing the CRS to a Decimal Degree projection.


title: Mapnik: mapnik.has_pycairo() returns False. How to fix this on Windows for Python27?

text: I'm trying to modify rundemo.py (comes canned with Mapnik Windows binaries) to render output files at a higher resolution, but it seems I still have issues with my Mapnik-to-Cairo setup. I've got mapnik and its python bindings working, and I've also got cairo and these python bindings working (..I think). If I start Python, I can do the following without issue:               E:\x_FOSS\mapnik\v2.0.0\mapnik-2.0.1rc0\demo\python>python.exe     Python 2.7.1 (r271:86832, Nov 27 2010, 18:30:46) [MSC v.1500 32 bit (Intel)] on win32     Type "help", "copyright", "credits" or "license" for more information.     >>> import mapnik     >>> import cairo     >>>      Meanwhile, back in rundemo.py line 351, it fails on this condition:               if HAS_PYCAIRO_MODULE and mapnik.has_pycairo():      And this.. `print mapnik.has_pycairo()` **# returns False**. But this conflicts with what I'm reading in the Mapnik blog which says these Windows binaries are built with Cairo backend support. Does anyone know what conditions have to be in place for `mapnik.has_pycairo()` to return True? * * * **[Update, are these build, or runtime instructions?]** Per Zero_Qualms' advice, I downloaded the mapnik source and found reference to `HAVE_PYCAIRO` in %\mapnik-v2.0.1\bindings\python\mapnik_python.cpp. I'm not familiar with C++, so I'm not sure if these are build instructions or runtime instructions. First, line 86 has this block following several `#include` statements:               #if defined(HAVE_CAIRO) && defined(HAVE_PYCAIRO)     #include <pycairo.h>     static Pycairo_CAPI_t *Pycairo_CAPI;     #endif      Then there are a couple insignificant calls to `defined(HAVE_PYCAIRO)`, followed by this `has_pycairo()` property at line 347:               // indicator for pycairo support in the python bindings     bool has_pycairo()     {     #if defined(HAVE_CAIRO) && defined(HAVE_PYCAIRO)         Pycairo_CAPI = (Pycairo_CAPI_t*) PyCObject_Import(const_cast<char *>("cairo"), const_cast<char *>("CAPI"));         if (Pycairo_CAPI == NULL){             /*               Case where pycairo support has been compiled into               mapnik but at runtime the cairo python module                is unable to be imported and therefore Pycairo surfaces                and contexts cannot be passed to mapnik.render()              */              return false;         }         return true;     #else         return false;     #endif     }      That block comment seems to imply my mapnik-pycairo linkage just isn't correct, and that it can be fixed. Next I searched for `PyCObject_Import` and found the following block in this file: %\mapnik-v2.0.1\bindings\python\python_cairo.cpp:               void register_cairo()     {         Pycairo_CAPI = (Pycairo_CAPI_t*) PyCObject_Import(const_cast<char *>("cairo"), const_cast<char *>("CAPI"));         if (Pycairo_CAPI == NULL) return;              boost::python::converter::registry::insert(&extract_surface, boost::python::type_id<PycairoSurface>());         boost::python::converter::registry::insert(&extract_context, boost::python::type_id<PycairoContext>());     }      Blast! It seems reasonable that I could create the registry entries manually, but ..what entries?


title: what is the visualization technology behind cartoDB?

text: does cartoDB uses mapserver/geoserver/other? what is the technology behind its visualization? thanks allot!


title: How to do geo coding with openlayers?

text: How can I achieve Geocoding with the JavaScript OpenLayers API, in the same fashion as the Google Maps API?


title: Using group by to set a PostGIS datasource using QGIS Python API

text: I am writing a QGIS Plugin and wish to obtain data from a PostGIS table using the Group By clause using something like the following code:               uri = QgsDataSourceURI()     uri.setConnection('localhost', '5432', 'postgis', 'postgres', 'xxx')     Query="select param_id, st_union(the_geom) as the_geom from myTable group by param_id"     uri.setDataSource('', '(%s)' % postGisQuery, "the_geom",'','gid')     vlayer1 = QgsVectorLayer(uri.uri(), 'testlayer1', "postgres")       This gives the error "No key field for query given." because it is looking for the gid key field. But if gid was put in the select statement of the Query, we won't get aggregated rows, just distinct rows. If I try an aggregation function on gid, I get the same error. If I try aliasing the aggregate, I get an "Erroneous query" error (even `count(gid) as gid`), same as if I put 'param_id' as the key field in setDataSource, because 'param_id' is not a primary key, I guess. In this case, forcing the key field in the uri seems to make the group by clause useless. Is there another workaround I haven't thought of or is this the expected behaviour? EDIT: in QGIS version 1.7.3 and my plugin code like below works, but not in my dev version (1.9.90):               Query="select param_id, st_union(the_geom) as the_geom from myTable group by param_id"     uri.setDataSource('', '(%s)' % postGisQuery, "the_geom",'','param_id')     vlayer1 = QgsVectorLayer(uri.uri(), 'testlayer1', "postgres")       So it allowed the use of param_id as a key, even though it isn't a primary key. Has something changed in recent versions that doesn't allow such code? In fact there doesn't seem to be an error message in the Log Messages, QgsVectorLayer just doesn't load the layer - using `QgsMapLayerRegistry.instance().addMapLayer(vlayer1)` Does anyone know if such code should work in the latest dev version (maybe there is a more recent version that allows this or the feature has been deliberately discontinued?) Ian


title: OpenGeo Dashboard.app won't start on Mavericks

text: I've installed the new OS X Mavericks and OpenGeo Dashboard 3.1 won't launch. My console has this message: "2013/10/25 12:25:36,909 PM OpenGeo Dashboard[11811]: The function `CGFontSetShouldUseMulticache' is obsolete and will be removed in an upcoming update. Unfortunately, this application, or a library it uses, is using this obsolete function, and is thereby contributing to an overall degradation of system performance." followed by "2013/10/25 12:25:38,254 PM com.apple.launchd.peruser.501[208]: (org.opengeo.suite.50368[11811]) Exited with code: 1" A problem with the dashboard and the new Mavericks? I have reinstalled OpenGeo Dashboard twice...


title: QGIS select feature, "add to current selection" feature

text: I am trying to perform a feature selection on qgis similar to "add to current selection" interactive selection method in Arcgis Desktop, is that possible? An example: I perform a spatial query and select a set of features, afterwards I need to manually add a couple of features to this set, if I use the feature selection button I lose my first selection, and start a new one. Any way of adding features to a selection? Many thanks everyone in advance,


title: Problem when creating database based on template_postgis

text: I have created a PostGIS template but when I try to create a database based on this template this is what happens: Connecting to the template:               \c template_postgis;          You are now connected to database "template_postgis" as user "postgres".          \dt                    List of relations       Schema  |      Name       | Type  |  Owner        ----------+-----------------+-------+----------      public   | spatial_ref_sys | table | postgres      topology | layer           | table | postgres      topology | topology        | table | postgres     (3 rows)          create database prueba with TEMPLATE template_postgis;     CREATE DATABASE      Connecting to database prueba:               \c prueba;     You are now connected to database "prueba" as user "postgres".          prueba=# \dt                   List of relations      Schema |      Name       | Type  |  Owner        --------+-----------------+-------+----------      public | spatial_ref_sys | table | postgres     (1 row)      You can see only the public schema was copied, no topology schema layer or topology were copied. I had no errors when I created the template:               #createdb template_postgis           #createlang plpgsql template_postgis           #psql -d template_postgis -f /usr/pgsql-9.2/share/contrib/postgis-2.1/postgis.sql            #psql -d template_postgis -f /usr/pgsql-9.2/share/contrib/postgis-2.1/postgis_comments.sql            #psql -d template_postgis -f /usr/pgsql-9.2/share/contrib/postgis-2.1/spatial_ref_sys.sql            #psql -d template_postgis -f /usr/pgsql-9.2/share/contrib/postgis-2.1/rtpostgis.sql           #psql -d template_postgis -f /usr/pgsql-9.2/share/contrib/postgis-2.1/raster_comments.sql            #psql -d template_postgis -f /usr/pgsql-9.2/share/contrib/postgis-2.1/topology.sql            #psql -d template_postgis -f /usr/pgsql-9.2/share/contrib/postgis-2.1/topology_comments.sql      Could you please help me? Thanks in advance


title: Postgres casting between varchar and double precision suddenly stopped working

text: I have been importing data into PostgreSQL tables all day using ogr2ogr and then spatially enabling the data with no issues up to now. I have been cutting and pasting the code and this is what I have been using:               ALTER TABLE mytable ADD COLUMN geom geometry(POINT,4326);     UPDATE mytable SET geom = ST_SetSRID(ST_MakePoint(cast (lon AS double precision),cast (lat AS double precision ) ),4326);      On the last table, I suddenly get this error:               ERROR:  invalid input syntax for type double precision: "#VALUE!"      Nothing has changed, at least as far as I know, I was importing a csv file, same as all the others. I checked the csv file for spaces or unusual characters but it looks in order > > Column lon, lon are of type **character varying** here is a sample of the data               id,name,name_alt,lat,lon,description     1,Abbott Seamount,,31.8,174.3,A seamount located north of Midway Islands (US) in the North Pacific Ocean; name origin unknown.     2,Abra Canyon,,17.51667,120.36667,An undersea canyon located off the northwest coast of the Philippine island of Luzon in the South China Sea; named after the Abra River     3,Abraham Canyon,Abraham Sea Valley,52.61667,172.33333,An undersea canyon in the area of the Aleutian Islands in the North Pacific Ocean; named in association with Abraham Bay; Attu Island; Alaska (US).      4,Abrolhos Bank,,-18.5,38.75,A bank located off the coast of Brazil in the Atlantic Ocean; upon which sits the Arquipélago dos Abrolhos (archipelago).      5,Abrolhos Seamounts,Abrolhos Ridge; Cadeia dos Abrolhos,-17.33333,36.5,Seamounts located off the coast of Brazil in the Atlantic Ocean; named in association with the Abrolhos Bank and Arquipélago dos Abrolhos (archipelago).      6,Abubacer Ridge,,36.8,1.65,An undersea ridge located off the coast of Spain in the western Mediterranean Sea; named for Ibn Tufayl (known as Abubacer to Europeans); 12th-century Spanish-Arab philosopher and physician.      7,Aby Canyon,,4.25,-3.66667,An undersea canyon located off the coast of Côte d'Ivoire; named in association with Lagune Aby; Côte d'Ivoire (lagoon)      What am I missing here?


title: How much permgen space should i add for using 5 geoserver

text: My computer is windows 7 64 bit with 4 gb ram. I am using 5 geoserver in a single apache. How much amount of space i should added..


title: How to Lock Grid and Scalebar in QGIS Atlas

text: Have just started using Valmeira and I am looking at generating an Atlas. I am using irregularly sized features as the index. The challenge I have, is when I set up the template page and scalebar and grid, if I go to another page (with a different sized feature) the resulting scale change makes the scalebar push out, and the grid become too small too read. There are two workarounds -   1. Use regularly sized index features - then Grid and Scalebar size are not an issue, as the scale doesn't change between pages.   2. Group the Features by their approximate scale to reduce the amount of changes in scale, and break the project into a number of atlases. My preferred option would be if there were a 'Disconnect Object from Atlas' style of option. This would mean that I could set the Scalebar and Grid for that page, then disconnect them from the Atlas automaitc refresh once I have completed them.. Cheers.


title: TileStache : Unknown geometry type in geoJSON file

text: I'm currently trying to set up a **TilSstache** server on Ubuntu 12.04 to serve some GeoJSON tiles from a shapefile/JSON file (administratives boundaries) for a web mapping project using the Leaflet framework, but I keep having **weird geoJSON files** with **"Unkown" geometry type and "null" coordinates** :               "geometry": {"type": "Unknown", "coordinates": null}.      Here's the complete story : I struggled a bit to install everything, but the Tilestache instance seems to be up and running. When using the following `proxy` provider in `tilestache.cfg`, **the OpenStreetMap tiles are correctly displayed** :               "osm":     {         "provider":         {         "name": "proxy",         "url": "http://tile.openstreetmap.org/{Z}/{X}/{Y}.png"         }     }      The next step was therefore to create the geoJSON tiles from either the shapefile or the JSON file. I'm using one of the "official" shapefile for french administratives boundaries, provided by the IGN (the National Geographic Institute of France), available HERE. This shapefile is provided in `EPSG:2154` (Lambert-93), a projection not directly supported by TileStache. QGIS 2.2.0 was used to project it to `EPSG:4326` (WGS84), save it as a new shapefile and also export it as a JSON data file. The resulting shapefile (`departement.shp` and the associated files) and JSON data (`departement.json`) are available for download HERE. I then set up the following config file to test both approaches (file paths were modified) : For the JSON data :               "test":     {       "provider":       {         "name": "vector",         "driver": "GeoJSON",         "parameters": {"file": "/XXX/XXX/departement.json"},         "properties": ["CODE_DEPT"]       }     }      For the shapefile :               "test2":     {       "provider":       {         "name": "vector",         "driver": "shapefile",         "parameters": {"file": "/XXX/XXX/departement.shp"},         "properties": {"CODE_DEPT": "code_dept"}       }     }      **The TileStache server run smoothly with this parameters (no errors thrown in the log file)**. However, when trying to access the relevant URLs (`http://127.0.0.1/test/0/0/0.geojson` or `http://127.0.0.1/test2/0/0/0.geojson`), i get a geoJSON like this :               {"type": "FeatureCollection",           "features": [     {"geometry": {"type": "Unknown", "coordinates": null}, "type": "Feature", "properties": {"NOM_DEPT": "AIN", "CODE_DEPT": "01"}},     {"geometry": {"type": "Unknown", "coordinates": null}, "type": "Feature", "properties": {"NOM_DEPT": "AISNE", "CODE_DEPT": "02"}},      {"geometry": {"type": "Unknown", "coordinates": null}, "type": "Feature", "properties": {"NOM_DEPT": "ALLIER", "CODE_DEPT": "03"}},      ...     ...     ]}      I thought it may be related to projection problem, so **I tried to add the following parameter for the two providers, to specify the projection used** :               "projection":"WGS84"      However :   * **The unknown geometry type / the null coordinates is still there.**   * **Many features are now absent in my geoJSON** (many departements are missing when I access `http://127.0.0.1/test/0/0/0.geojson` or `http://127.0.0.1/test2/0/0/0.geojson`). I'm totally new to GIS (I didn't know about QGIS, TileStache or even tiles 2 weeks ago), so I may have missed something really obvious... Thanks !


title: Problem validating parent model with nested models in modelbuilder

text: I have a parent model with several nested models. I'm having trouble getting the parent model to validate. The issue is with sub model 4. The two screenshots are of the parent model and sub model 4. The output from sub model 3 feeds into parameter “Merged_BBK_Assets” in sub model 4. There are no other input model parameters in sub model 4. If I run the first three sub models then open, validate and run sub model 4 from within its model window, it runs fine. The Merged_BBK_Assets is created in sub model 3 and does not exist till then.   ![Sub Model 4](http://i.stack.imgur.com/Thvth.png) ![Parent](http://i.stack.imgur.com/XavLt.png)


title: problem transforming pgr route to a line - unexpected extra lines from St_MakeLine

text: I have created a line geometry from a pgr_dijkstra querry and I am getting an unexpected looking line in qgis. The following image illustrates the problem: ![enter image description here](http://i.stack.imgur.com/yFc70.png) The magenta line is the line created from mapping the route from my table where I have created a line geometry using the following code:               update analysis.routing     set     routeline = (select st_makeline(pgr.geom) ,     from (select * from pgr_dijkstra('ways',61789,64493)) as pgr)      The dotted red line is the line created in qgis from calling the same select statement in the DB manager and adding the results to canvas:               select * from pgr_dijkstra('ways',61798,64493)      Does anyone have any recommendations on how to fix this problem with ST_MakeLine. That is what appears to be the problem. It looks like ST_MakeLine is connecting additional points which is adding the additional lines above that are not part of the pgr network. Thanks for your help.


title: Finding normal vector to DEM

text: I have a DEM grid (Surfer, or whatever format) and an x,y position on that grid. I also have a series of depths from the surface. I need to calculate, for each depth, the x,y,z of the projected point, normal to the DEM. The application is the display of a new type of geophysical dataset where the energy returns are normal to the surface. The dataset has 256 depths for each x,y position, and there are thousands of x,y positions. Hoping that there's a simply way to do this? Thanks :)


title: How do I convert coordinates to raster (bsq, bil) lookup?

text: I have a .bil file, which has `interleave = bsq` in the hdr file. I'm having trouble finding any insight into how I would essentially do a lookup on the file. I have no problem reading a specific chunk of bytes randomly from the file, but am unsure of how I would go about converting a lat/lon pair to a specific offset within the file. I've only seen solutions in R, but ideally I'd like to be able to do this in python or javascript, using proj4 or something to that effect. Any advice?


title: How to query DEM elevation values in a web application based on cursor/click point location?

text: I work for a city government that has LiDAR derived DEM rasters which we would like to make accessible to other departments in the city. The problem is that many departments don't have anyone with GIS skills and/or licenses. We are looking to develop a simple, easy to use application that will just return elevation values from the DEM based on cursor click/point location. We would like to do this without having to serve out the high resolution DEM dataset for speed and ease of use of the application. We currently have GeoPDFs that can return elevation data as long as users have the TerraGo toolbar installed. This is expensive and extremely cost- prohibitive so we are looking for a cheap option that can be used by employees and the public alike. We are just looking to query based on XY location of a point which will return the elevation value based off of the DEM. Has anyone done this or have advice on the best way to proceed? We do have access to Google Maps Engine (the Elevation API is SRTM so resolution isn't good enough for our purposes), as well as ArcGIS Online. Also worth mention: We are currently running Oracle Spatial but are transferring to SQL Server 2012 with ESRI's LGIM for our spatial data(geodatabases and feature classes).


title: Why are legends absent from ArcGIS Online map service?

text: With reference to an older thread I posted, (How to use map service legend on ArcGIS Online?) We have installed a security certificate for our host machine, but the absence of a legend using the AGOL templates persists. As an example, this is a **public map** we have created. Aside from the gray canvas basemap, it uses a map service hosted in our machine for " **road classification** ": What would be the rationale behind this legend issue and what can be done?


title: Which are the most useful NMEA sentences for location based services

text: I am pretty new to GIS, I have been learning about all the different NMEA sentences from this website http://aprs.gids.nl/nmea/#allgp What is $GPGSV and $GPGSA and why are they useful to location based services system?


title: No clue of the Georeferencing tool in gvSIG

text: I have kept trying to locate the georeferencer in gvSIG but I could not find it. Some tutorials say it is under Tools from the main menu! I installed several versions of gvSIG, both on Windows and Linux but the georeferencer does not appear anywhere. I also tried to install a georeferencing plugin from `http://downloads.gvsig.org/download/gvsig- desktop/pool/org.gvsig.georeferencing/` but after I installed it, the gvSIG won't start. Any idea?


title: Conversion From Transverse Mercator To Longitude & Lattitude

text: I have to transform (7-parameters datum transformation) coordinates from northings & eastings to lattitude & longitude. Following are details of source coordinate system(i.e. northings & eastings one) :               1) Projection Details:       1.1)Projection Name      1.2)Origin Lat/Lon      1.3)Origin North/East      1.4)Scaling Factor     2)Datum Name: AGD84     3)7-parameter Datum Details:      3.1)Rotn X:0.292      3.2)Rotn Y:-0.443      3.3)Rotn Z:-0.277      3.4) Trans X:117.763      3.5) Trans Y:-51.51      3.6)Trans Z:139.061      3.7)Scaling Factor:-0.03939657799319541      Please note that I am using `"Proj4Js"` for the above mentioned transformation for which following is the "defs" settings:               ***Proj4js.defs["<EPSG-CODE>"] = "+proj=<someProjectionImplementation>       +towgs84=-117.763,-51.51,139.061,0.292,-0.443,-0.277,-0.03939657799319541 +no_defs";***      Could anyone confirm that is there any issue in this setting? Thanks, Gajendra


title: Mapnik style editor

text: Is there any tool or editor which can shorten the cycle for experimenting on different styles for mapnik. It would be wonderful with something like cloudmade's editor, but with the possibility to download the result as xml or cascadenik. What i want is:   * Short iteration time (ie. few steps and fast generation time for newly changed style) preferably live edit   * Tool which gives help to generate the styles (both graphical tools and word completion would be great alternatives. Any ideas?


title: show message if map service unavailable

text: If any of my ESRI REST map services are down I want to be able to tell the user, at minimum with an alert or message. I thought I could use onLayerAddResult(layer, error) but it doesn't seem to be getting called if the layer based on a stopped/down service is added to the map, only when a layer is successfully added. I see some examples using esri.request but it doesn't seem like I should have to do an extra call to the service since in the Firebug console I see these calls like this to each service: > http: > //myserver/ArcGIS/rest/services/Streets/Mask/MapServer?f=json&dpi=96&transparent=true&format=png8&layers=show%3A-1&callback=dojo.io.script.jsonp_dojoIoScript2._jsonpCallback With a response of either: > *dojo.io.script.jsonp_dojoIoScript2._jsonpCallback({"serviceDescription" : > null,"mapName":"Layers","description":"","copyrightText":"","layers":[{"id":0,"name":"LMBCO_Mask","parentLayerId":-1,"defaultVisibility":true,"subLayerIds":null}],"tables":[],"spatialReference":{"wkid":2263},"singleFusedMapCache":false,"initialExtent":{"xmin":806353.163236593,"ymin":139639.306949227,"xmax":1165380.94101437,"ymax":327833.751393672,"spatialReference":{"wkid":2263}},"fullExtent":{"xmin":824755.941014371,"ymin":38597.6402825604,"xmax":1146978.16323659,"ymax":323319.862504783,"spatialReference":{"wkid":2263}},"units":"esriFeet","supportedImageFormatTypes":"PNG32,PNG24,PNG,JPG,DIB,TIFF,EMF,PS,PDF,GIF,SVG,SVGZ,BMP","documentInfo":{"Title":"","Author":"","Comments":"","Subject":"","Category":"","Keywords":"","AntialiasingMode":"Normal","TextAntialiasingMode":"Force"},"capabilities":"Map,Query,Data"});* if the service is up and running, or: > > *dojo.io.script.jsonp_dojoIoScript2._jsonpCallback({"error":{"code":400,"message":"Unable > to complete operation.","details":["Service 'NYCStreets/LMBCO_Mask' of type > 'MapServer' does not exist or is inaccessible."]}});* if the service is not available. Is there a way to handle this reponse that is already coming back? Thanks for your help! Chad Wilcomb


title: OSM_id missing from ways table created by Osm2pgrouting

text: GIS newbie here so excuse my ignorance if this has a simple answer. I have been playing about with pgRouting and using Osm2pgrouting to load OSM files from Geofabrik. However I just noticed that the OSM_id isnt being created by Osm2pgrouting. I'm trying to exclude some roads in pgRouting that are tagged "motors=no" in OSM but cant because I cant match them without the osm_id field?


title: What is the best colour-ramp to use for elevation?

text: This is maybe rather subjective as it concerns style, but I would like to have a few decent colour ramps for my DEM when it is layed over my hill-shading. I have created the Swiss hillshade model as shown here, but I'm not too convinced that the colour scheme is good. Any examples or screenshots would be great.


title: Adding a delimited text layer in QGIS 1.8 only imports a fraction of the total rows

text: I am trying to add a delimited text layer, and I think I've formatted it correctly as a CSV file. There are no errors when I import the data. But even though there are 1300 rows total, only 189 appear in the table when I right click on the layer and open the attribute table. Is there any reason why it would import so few rows? Thanks!


title: how to create a network from gps tracks?

text: I'm trying to build a network from several GPS-trakcs, to digitalize hiking routes on my land. I have several GPX-files from routes I've tracked with the gps. I managed to display them in Arcmap 10.2 and also create lines between the points. But my knowledge/skills kind of fall a bit short from here on. so can somebody give me an easy step-by-step how to create a network from gps- tracks? one of my problems so far is this: several of the routes have some line segments overlapping, how do i make the overlapping segments in to one(preferably the averaged/middle one) and create a junction where they split? see picture "Q1"![black red and green symbolizes three different track logs](http://i.stack.imgur.com/61KRI.png) **My second challenge is this:** I want to add an attribute with information about average speed between my gps point, this is a simple operation but since there is a lot of gps points i need to be able to do it fairly automatic. So far what I've tried is "points to line" to get a line, then i used "split line at vertices" to get separate line segments between each gps point. Then I tried to join tables, between the line segments and gps-points, so that the line segments would have time information on its start and end, then it would be easy to calculate the speed for every row. So i tried using "spatial join", the problem here is that it creates separate rows for the join of "line segement and its start point" and the "line segement and its end point" see table: ![enter image description here](http://i.stack.imgur.com/zjDuB.png) marked in yellow is one line segment and its join to two gps points(see Shape Length is the same and DateTime is different) My goal is as mentioned to find the average speed the used on this line segment. If my approach to the problem is totally off I'm open to try other methods!


title: ModelBuilder: How to extract data by date, export to separate tables and interpolate tables to rasters

text: I've been using ArcGIS for a couple years, but when it comes to model builder I've only worked with very simple models. I'm hoping someone that is much smarter than myself can help me in building a complicated model. I've got a very large table of data (about 300k rows) that contains points with lat/lon data, wind speed and date/time. There are about 60 different points with wind values that have been recorded for each hour during 2012. What I'd like to do is build a model that will:   1. Iterate through each row and select/export the data associated with first hour (00Z) Jan 1st 2012 (about 60 rows) to a separate table, then the 2nd hour (01Z) of Jan 1st 2012 to its own separate table and so on through all of 2012.    2. Perform Inverse Distance Weighting on all of the hourly tables with lat/lon and wind values that were created from step 1 above   3. Reclassify the rasters so that points in the raster that have wind values of 5 or greater are reclassified to "1" and all other cell's wind values are reclassified to "0"   4. Then iterate through the rasters created in step 3 and use something like "Collect Values" & "Cell Statistics" to sum up all of the rasters into a final summed raster for all hours in 2012. Was going to include screenshot of my data, but I guess since this is my first post I'm not trusted yet. I've broken the date down several different ways (just date, just time, date-time, and date-time in text format) because I've had issues trying to iterate through and select the data in my original, very large table only by date. Hopefully I provided enough information, but if not I will try to go into more detail. Thanks in advance for any assistance.


title: Conversion from Shapefile to TopoJSON format results in incorrect map when rendered using d3js

text: I am trying to build a map of Canada showing FSA boundaries using d3js. I got shape data from Statistics Canada [http://www12.statcan.gc.ca/census- recensement/2011/geo/bound-limit/bound-limit-2011-eng.cfm]. To convert the SHP file to topojson, I am following the process from Mike Bostock's Tutorial 'Lets Build a Map' [http://bost.ocks.org/mike/map/]. First convert the map to GeoJSON using ogr2ogr, then using topojson in Node.js, convert the GeoJSON file to a TopoJSON file. I am executing the commands:               > ogr2ogr -f "GeoJSON" -s_srs EPSG:21781 -t_srs EPSG:4326 output.geojson input.shp      Followed by the command:                > topojson --id-property CFSAUID -p name=PRNAME -p name -o output.topojson output.geojson      The resulting map that is generated renders some provinces properly, however, the following provinces show up as jumbled lines.   * Ontario   * Newfoundland    * Quebec    * Manitoba    * Saskatchewan    * Alberta    * British Columbia Here is a basic version of my d3 call.               var g = svg.append("g").attr("id", "fsa");     d3.json("can/output.json", function (error, json) {         var counter = 0;         var subunits = topojson.object(json, data);         g.selectAll("path")         .data(subunits.output)         .enter()         .append("path")         .attr("d", path)             });      This is what the map looks like: https://www.dropbox.com/s/k776hg2doynurxi/canada.png


title: Getting: org.opengis.referencing.NoSuchIdentifierException: No transform for classification "Transverse Mercator Zoned Grid System" with geotools

text: Using geotools I've been trying to create a referenced envelope after converting/reprojecting from WGS 84 to UTM with specific zones and it seems that either I need to do something with authority factories or figure out some kind of a workaround for this. Here are the errors:               Exception in thread "main" org.opengis.referencing.NoSuchIdentifierException: No transform for classification "Transverse Mercator Zoned Grid System".     at org.geotools.referencing.operation.DefaultMathTransformFactory.getProvider(DefaultMathTransformFactory.java:290)     at org.geotools.referencing.operation.DefaultMathTransformFactory.getDefaultParameters(DefaultMathTransformFactory.java:316)     at org.geotools.referencing.factory.epsg.DirectEpsgFactory.createCoordinateOperation(DirectEpsgFactory.java:2610)     at org.geotools.referencing.factory.BufferedAuthorityFactory.createCoordinateOperation(BufferedAuthorityFactory.java:1041)     at org.geotools.referencing.factory.epsg.DirectEpsgFactory.createCoordinateReferenceSystem(DirectEpsgFactory.java:1955)     at org.geotools.referencing.factory.BufferedAuthorityFactory.createCoordinateReferenceSystem(BufferedAuthorityFactory.java:783)     at org.geotools.referencing.factory.AuthorityFactoryAdapter.createCoordinateReferenceSystem(AuthorityFactoryAdapter.java:801)     at org.geotools.referencing.factory.FallbackAuthorityFactory.createCoordinateReferenceSystem(FallbackAuthorityFactory.java:640)     at org.geotools.referencing.factory.FallbackAuthorityFactory.createCoordinateReferenceSystem(FallbackAuthorityFactory.java:644)     at org.geotools.referencing.factory.AuthorityFactoryAdapter.createCoordinateReferenceSystem(AuthorityFactoryAdapter.java:801)     at org.geotools.referencing.factory.ThreadedAuthorityFactory.createCoordinateReferenceSystem(ThreadedAuthorityFactory.java:731)     at org.geotools.referencing.DefaultAuthorityFactory.createCoordinateReferenceSystem(DefaultAuthorityFactory.java:179)     at org.geotools.referencing.CRS.decode(CRS.java:489)     at org.geotools.referencing.CRS.decode(CRS.java:417)      Here is the block of code that caused the error, it is related to trying to use CRS.decode():                   MathTransform mathTransform = CRS.findMathTransform(aoi.getCoordinateReferenceSystem(), CRS.decode("EPSG:" + epsg));         DirectPosition2D srcDirectPosition2D = new DirectPosition2D(aoi.getCoordinateReferenceSystem(), aoi.getMinX(), aoi.getMinY());         DirectPosition2D destDirectPosition2D = new DirectPosition2D();         mathTransform.transform(srcDirectPosition2D, destDirectPosition2D);         DirectPosition2D srcMaxDirectPosition2D = new DirectPosition2D(aoi.getCoordinateReferenceSystem(), aoi.getMaxX(), aoi.getMaxY());         DirectPosition2D destMaxDirectPosition2D = new DirectPosition2D();         mathTransform.transform(srcMaxDirectPosition2D, destMaxDirectPosition2D);      This is what I got when I printed out the coordinate reference system:               PROJCS["WGS 84 / UTM zone 18N",      GEOGCS["WGS 84",      DATUM["World Geodetic System 1984",        SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],        AUTHORITY["EPSG","6326"]],      PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],      UNIT["degree", 0.017453292519943295],      AXIS["Geodetic latitude", NORTH],      AXIS["Geodetic longitude", EAST],      AUTHORITY["EPSG","4326"]],      PROJECTION["Transverse_Mercator", AUTHORITY["EPSG","9807"]],      PARAMETER["central_meridian", -75.0],      PARAMETER["latitude_of_origin", 0.0],      PARAMETER["scale_factor", 0.9996],      PARAMETER["false_easting", 500000.0],      PARAMETER["false_northing", 0.0],      UNIT["m", 1.0],      AXIS["Easting", EAST],      AXIS["Northing", NORTH],      AUTHORITY["EPSG","32618"]]      Apparently I may need to set a comparison tolerance when I use CRS related functions, but I'm not sure how to go about doing that. http://docs.geotools.org/latest/javadocs/org/geotools/referencing/CRS.html


title: Unable to cast COM object (IMapServer --> IMapServerObjects)

text: I have just deployed an AGS-based application to our web server. This app worked great in development, but after deployment to the server I get the following errror: `Unable to cast COM object of type 'System.__ComObject' to interface type 'ESRI.ArcGIS.Carto.IMapServerObjects'. This operation failed because the QueryInterface call on the COM component for the interface with IID '{651521F3-AAFF-42B9-8C2D-1980625A13BF}' failed due to the following error: No such interface supported (Exception from HRESULT: 0x80004002 (E_NOINTERFACE)).` The C# code generating this error is:               IServerContext context = som.CreateServerContext(ConfigurationManager.AppSettings["mapService"], "MapServer");     IMapServer mapServer = (IMapServer)context.ServerObject;     IMapServerObjects serverObjects = (IMapServerObjects)mapServer;      For some reason, the `IMapServer` instance cannot be casted to `IMapServerObjects`. What could be wrong?


title: Which interpolation method for noise mapping?

text: I want to create a noise map with measured point data. But i want to have a planar representation, so I need an interpolation method. Which one is the 'best' one for measured noise data? I prefer to 'IDW' interpolation, but I´m not sure about that. Any ideas? thanks in advance!


title: Main OpenStreetMap-based sites?

text: I don't like the way the standard OSM maps look (too much details), and would like to check if some OSM-based services are better than Google Maps. "List of OSM based Services" doesn't list them by popularity. What are the main sites that use OSM maps for general purposes (ie. alternative to GoogleMaps)? Thank you.


title: arcpy update cursor: can't connect to a database table

text: I am using updateCursor to retrieve some rows from a non spatial table using a connection of type "Spatial Database Connection" (That creates a .sde file) from ArcCatalog. However, it doesn't work, I have the following error: > (RuntimeError('ERROR 000972: Table must have Object IDs in order to be > sorted.',))) Even if my table is not spatial (there is no objectid field). Also, when I try to connect to another non spatial table from another database, It works and I can retrieve the rows. Does anyone have a solution please? D.M


title: Third party base layer other than Google map

text: Can we use third party base map (images tiles) for displaying vector data?, like we use google map as base layer for displaying vector data in OpenLayers application.


title: What's the Semantic Web got to do with GIS?

text: A few years ago there was a lot of discussion about the Semantic Web. Are any of the concepts underlying the Semantic Web coming into fruition with respect to GIS? Or is it just another hippie dream? > I have a dream for the Web [in which computers] become capable of analyzing > all the data on the Web – the content, links, and transactions between > people and computers. A ‘Semantic Web’, which should make this possible, has > yet to emerge, but when it does, the day-to-day mechanisms of trade, > bureaucracy and our daily lives will be handled by machines talking to > machines. The ‘intelligent agents’ people have touted for ages will finally > materialize. \- Tim Berners-Lee Sounds like Tim may have been fishing the same streams as Richard Brautigan, whose poem published in 1963, All Watched Over by Machines of Loving Grace has everything but the cloud:               I'd like to think (and     the soner the better!)     of a cybernetic meadow     where mammals and computers     live together in mutually     programming harmony     like pure water     touching clear sky.


title: Where can i find oceanic borders defined as a set of lat/long points?

text: I'm trying to write some code that, given a latitude and longitude, will report which ocean (if any) this point is in. Is there a public resource available that defines the borders of the oceans as a set of lat/long points that I can download?


title: How to get easting & northing of cell in raster using ArcPy?

text: It might be very simple task. But I would appreciate if someone can comment on How I can get easting and northing of a single cell raster using arcpy as this value is used for in further steps. Thanks!


title: How to make sure OnSketchFinished always works?

text: In my VBA application i have a listener setup that listens for OnSketchFinished, when the polygon is created a messagebox appears and then a userform. My problem is that the listener will work 90% of the time, not 100%. Here are the snippets of code in ThisDocument:               Private WithEvents EditEvents As Editor     Private m_pEditor As IEditor     Public Sub SetEvents()               ' This is the listener for the Editor     Dim pID As New UID     pID = "esriEditor.Editor"     Set m_pEditor = Application.FindExtensionByCLSID(pID)     Set EditEvents = m_pEditor          End Sub               Private Function MxDocument_OpenDocument() As Boolean     Call SetEvents     End Function               Private Sub EditEvents_OnSketchFinished()         Call frmCityDev.EditEvents          End Sub      In my form I have code for EditEvents which brings up a msgbox to notify that a new polygon has been created. Am I missing something to make sure that the listener will always work? Cheers


title: How to correctly use cursor.execute() with PostGIS?

text: I do not understand how the operator works some times... Here is an example that works,               "SELECT ST_Distance(ST_GeomFromText('POINT(-3.16496271127842 55.9262620593642)'),ST_GeomFromWKB(a.geometry)),ST_AsText(ST_Centroid(a.geometry)) FROM <table> As a;"      and here ia an example that fails,I use a simple execute methos like this,               cc = conn.cursor()     a = ('POINT(-3.16496271127842 55.9262620593642)')     b = "SELECT ST_Distance(ST_GeomFromText(%s),ST_GeomFromWKB(a.geometry)),ST_AsText(ST_Centroid(a.geometry)) FROM <table> As a;"     cc.execute(b,a)      sometimes, the operator works fine and sometimes it does not. I have read documentation about the operators and how to use python with SQL but I have not found WHY it acts like this in a random kind of way. Does anyone have any hints and tips? Cheers, A


title: ArcGIS 10 MapControls in MS Access, can it be done?

text: In the past (9.3 ArcGIS) I was able to develop forms in MS Access where I could drop a MapControl on and have GIS functionality all within Access as long as ArcMap was installed on the same machine. I now have ArcGIS 10 and Office 2003 and when I start Access go to VB editor and create a form I am unable to add a MapContol, it says I do not have the appropriate license. I've had a search on the ESRI website and Help and I have not found a clear answer. Am I correct in understanding that they have changed something behind the scenes with regards to the MapControl, one cannot simply use it anymore? I've tried dropping an ESRI LicenseControl on the form and I get the same message rejecting me as I have no license to use the LicenseControl! It would be good if someone clarifies this and explain how one can embed a MapControl into a MS Access form?


title: QGIS FIle Read Error when opening an existing project

text: My computer will not open a project i have been working on for weeks! it is coming up with an error of > Project file read error: unexpected end of file at line 1 column 1 for file > C:/Users/xxxx/Documents/Nxxxx/Exxxx.qgs Is there anyway to get it back?


title: How to drop ArcSDE connection using Python scripting?

text: I want to be able to control connections to ArcSDE via python. Here is my script so far using ESRI's example.               import arcpy     from arcpy import env     env.workspace = r"SDEConnection"      fcList = arcpy.ListFeatureClasses()     print str(fcList) + "\n"     env.workspace = ""      arcpy.ClearWorkspaceCache(r"SDEConnection")     print arcpy.GetMessages() + "\n"      But this isnt working .. it crashes my Python window :(


title: Display Table of contents of the map in ArcMap in the online published Map/Web map application

text: The following is the sample layer information of the map in ArcMap 10.1: ![enter image description here](http://i.stack.imgur.com/eOYZM.jpg) How can I have the same Table of contents in the published map online or in web map application, so that user can select/unselect layer's content? Is there any feature I need to enable to have all contents appear for online user? Do I need to implement JavaScript API?Is there any sample available? Thanks


title: Dissolve polygons in QGIS (GEOS) with tolerance

text: I have a polygon vector layer representing the area inundated in a flood. The layer has been produced in a such a way that some of the areas are composed of multiple smaller touching polygons (see screenshot below, labelled "original data"). I would like to dissolve the polygons where they touch. If I use the dissolve geoprocessing tool in QGIS's fTools, which uses the GEOS library, I get some unexpected results (see "dissolve with QGIS fTools"). It seems that GEOS doesn't use a threshold in determining if polygons are touching, as mentioned in these issues: http://trac.osgeo.org/qgis/ticket/3126 http://trac.osgeo.org/qgis/ticket/3184 http://hub.qgis.org/issues/2806 The same approach in ArcGIS gives the expected result (see "dissolve with ArcGIS 10.1"). Apparently this is a "feature" of GEOS, rather than a bug. Is there any way to fix the data I have to produce the expected result? A minor loss of precision (a few cm?) could be acceptable for me. Ideally the solution would be something I could write into a QGIS plugin, as a pre-processing step. Related question (with no answers): How to set tolerance value of Geoprocessing tool in QGIS? Link to sample data as a shapefile: https://www.dropbox.com/s/c8ns7wwfv289ad0/sample.zip ![dissolve polygons](http://i.stack.imgur.com/Ayteo.png) ![v.dissolve](http://i.stack.imgur.com/AA59l.png)