title: KML from PostGIS Geometry into Openlayers

text: im getting the KML value from my postgis database.. `st_askml()` And the output is like this: `"<Polygon><outerBoundaryIs><LinearRing><coordinates>-99.999999990686774,48.499999959021807 -99.999999990686774,47.999999988824129 -99.500000020489097,47.999999988824129 -99.500000020489097,48.499999959021807 -99.999999990686774,48.499999959021807</coordinates></LinearRing></outerBoundaryIs></Polygon>"` Now i need to display them on a Openlayers Map i got. Layer: `var KML_Layer = new OpenLayers.Layer.Vector('KML');` Formatting code:               function GetFeaturesFromKMLString (strKML) {     var format = new OpenLayers.Format.KML({             'internalProjection': options.projections.source, //EPSG:3785/900913             'externalProjection': options.projections.display //EPSG:4326         });          return format.read(strKML);     };      Adding the feature to the layer: `KML_Layer.addFeatures(value.st_askml);` Debuggin info: value.st_askml = The output format i pasted earlier in topic. After running it trough the `GetFeaturesFromKMLString()` its set to `[]` (No value here) How can i get this to display correctly on the map?


title: Link to google maps ignores url parameters for streetview

text: I have a web site where I need to generate link to Google Maps in Street View mode. I have available the coordinates to generate the URL. When I try this URL in a desktop browser, it does what I want, and takes me to a street in Norway, and activating street view:   https://maps.google.com/maps?spn=0.008443,0.022767&ctz=-60&t=h&z=16&layer=c&cbll=58.819625,5.75462&cbp=12,44.85,,0,18.41 When I try exaclty the same URL in a **mobile browser** , it just ignores the parameters, and redirects to http://maps.google.com . When I try another link that is not streetview, it works. I noticed it first in an iPad, and then tried in a Samsung Galaxy using Firefox, and same there. I even tried using Firefox on desktop using the mobile device emulator plugin emulating an iPad, and same there. So it is Google that does this if the user agent reports itself as a mobile device. Then I tried to navigate to a specific location in Google Maps on the mobile device, and clicked on the street view icon. Then a new window was opened with this URL:   http://maps.gstatic.com/m/streetview/?panoid=vZ_OD_yZLEJYC5wQQ9q5zg&cbp=0,247.65935203825757,,0,0 This works for this location, but I need to generate the URL based on coordinates, so I don't have a valid panoid to use. Is panoid the only way to use url parameters for displaying street view on mobile device? If so, is there any API available for generating a panoid based on coordinates?


title: Mosaic large dataset with R

text: I need to merge together 20 32 bit tiff. For that I tried ArcGis but it keeps crashing (the individual tiff are about 1GB). I need a finale Tiff in order to run some analysis and I figured that having 20 different file would be messy for the user who is going to do the analysis. I found a mosaic tool in the raster package in R but I can't really figure out how to use it. My first idea was to create a list of the tiff and then ask the mosaic tool to process it but it gives me error. > Error in function (classes, fdef, mtable) : unable to find an inherited > method for function "mosaic", for signature "list", "missing" I found this script online but it seems that many function doesn't work anymore               require(raster)     input.rasters <- lapply(list.files(pattern="^TmB50.*[.]tif$"), raster)     full.extent <- unionExtent(input.rasters)     bounding.raster <- raster(full.extent,crs=projection(input.rasters[[1]]))     res(bounding.raster) <- res(input.rasters[[5]])     resampled.rasters <- lapply(input.rasters, function(input.raster) {     target.raster <- crop(bounding.raster, input.raster)     resample(input.raster, target.raster, method="bilinear")     })     raster.mosaic <- mosaic(resampled.rasters, fun=max)      source http://www.nceas.ucsb.edu/scicomp/usecases/createrasterimagemosaic Any idea how I could do that?


title: Where can I find a map with the location of towns and cities in Southern Sudan?

text: For a project I have some name of towns and cities in the Southern part of Sudan, the region South Kordofan. Many of these towns I can't find on Google Maps, but I need to know their location. Example of towns I can't find are: Kuma, Mendi, Kalkutta, Um Serbeba and Al Nugra. Hope someone knows a good source?


title: Using OpenLayers to render custom tiles

text: I have a very large image which I want to display using OpenLayers so visitors can zoom in and pan. I have tiled it (1 image for the first layer, 4 for the second, 9 for the third, etc.) and I am trying to display it using a TMS layer. Following this information I have it working more or less. However, on some zoom levels, I can not pan all the way to the right and the image seems cut off (even though the tiles exist). I have set my image height and width in px as the layer's maxExtent. Is there a better way to do this? Thanks.


title: Declarations required for Mapbasic api calls

text: I am trying to read a serial port by using api calls createfile, readfile, getcommstate, setcommstate etc... this all requires the proper declaration of a structure "dcb"               typedef struct _DCB {      DWORD DCBlength;      DWORD BaudRate;      DWORD fBinary  :1;      DWORD fParity  :1;      DWORD fOutxCtsFlow  :1;      DWORD fOutxDsrFlow  :1;      DWORD fDtrControl  :2;      DWORD fDsrSensitivity  :1;      DWORD fTXContinueOnXoff  :1;      DWORD fOutX  :1;      DWORD fInX  :1;      DWORD fErrorChar  :1;      DWORD fNull  :1;      DWORD fRtsControl  :2;      DWORD fAbortOnError  :1;      DWORD fDummy2  :17;      WORD  wReserved;      WORD  XonLim;      WORD  XoffLim;      BYTE  ByteSize;      BYTE  Parity;      BYTE  StopBits;      char  XonChar;      char  XoffChar;      char  ErrorChar;      char  EofChar;      char  EvtChar;      WORD  wReserved1;     } DCB, *LPDCB;      My first question is how would you declare this in mapbasic? Then how to determine the structures length?


title: Measuring distance between points along a river system

text: I have a forked river with multiple x,y data points along it and I need to measure the distance between these points when travelling along the river. However the forked part is throwing a wrench into things. I need the total distance as you move from 1-2, 2-3, 3-4 etc. The problem is the points are not sequential along the river and some points are in one fork while the next points are in the other. This is in ArcMap10. Thanks


title: Image display/scaling in OpenLayers vs Google Maps

text: When I plot an image on OpenLayers, it looks like that image is stretched a little bit as compared to Google Maps. I was wondering if it is normal or there is something wrong with my code. OpenLayers: ![enter image description here](http://i.stack.imgur.com/uRTDd.png) Google Maps (maps.google.com): ![enter image description here](http://i.stack.imgur.com/Ru6Xr.png) Following is the XML :               <?xml version="1.0" encoding="UTF-8"?>     <kml xmlns="http://www.opengis.net/kml/2.2">       <Folder>         <name>Ground Overlays</name>         <description>Examples of ground overlays</description>         <GroundOverlay>           <name>Large-scale overlay on terrain</name>           <description>Overlay shows Mount Etna erupting                on July 13th, 2001.</description>           <Icon>             <href>http://code.google.com/apis/kml/documentation/etna.jpg</href>           </Icon>           <LatLonBox>             <north>37.91904192681665</north>             <south>37.46543388598137</south>             <east>15.35832653742206</east>             <west>14.60128369746704</west>             <rotation>59.1556640799496235</rotation>           </LatLonBox>         </GroundOverlay>                 </Folder>     </kml>      And following is my JavaScript code:               xml=jQuery.parseXML(oRequest.responseText)  //oRequest.responseText contains above xml              groundOverlays=jQuery(xml).find('GroundOverlay');              jQuery.each(groundOverlays,function(index, value){             var groundOverlay=new Object();             parentFolder=jQuery(groundOverlays).parent('Folder').find('name').text();                  groundOverlay.href=jQuery(value).find('href').text();             groundOverlay.name=jQuery(value).find('name').text();             groundOverlay.description=jQuery(value).find('description').text();             groundOverlay.north=jQuery(value).find('north').text();             groundOverlay.south=jQuery(value).find('south').text();             groundOverlay.east=jQuery(value).find('east').text();             groundOverlay.west=jQuery(value).find('west').text();             groundOverlay.rotation=jQuery(value).find('rotation').text();                            var options = {                    'opacity': 1.0,                  'isBaseLayer': false,                 numZoomLevels : 20 ,                 //aspectRatio:50,                  };             var image = new OpenLayers.Layer.Image(                 groundOverlay.name,                 groundOverlay.href,                 new OpenLayers.Bounds( groundOverlay.west, groundOverlay.south,  groundOverlay.east, groundOverlay.north),                 new OpenLayers.Size(1, 1),                 options                 );                  map.addLayers([image]);                   });


title: Can image layer be rotated in OpenLayers?

text: While working on Ground Overlays, I noticed that it has a tag called `<rotation>`                <?xml version="1.0" encoding="UTF-8"?>         <kml xmlns="http://www.opengis.net/kml/2.2">           <GroundOverlay>             <name>Large-scale overlay on terrain</name>             <visibility>1</visibility>             <description>Overlay shows Mount Etna erupting on July 13th, 2001.</description>             <Icon>               <href>http://kml-samples.googlecode.com/svn/trunk/resources/etna.jpg</href>             </Icon>             <LatLonBox>               <north>37.91904192681665</north>               <south>37.46543388598137</south>               <east>15.35832653742206</east>               <west>14.60128369746704</west>              <rotation>11.9556640799496235</rotation>             </LatLonBox>           </GroundOverlay>         </kml>      My question is - can an image layer be rotated in 2D view? Following is the code I am using to add an Image and set its bounding box.               var options = {                    'opacity': 1.0,                  'isBaseLayer': false,                 numZoomLevels : 20                   };             var image = new OpenLayers.Layer.Image(                 groundOverlay.name,                 groundOverlay.href,                 new OpenLayers.Bounds( groundOverlay.west, groundOverlay.south,  groundOverlay.east, groundOverlay.north),                 new OpenLayers.Size(1, 1),                 options                 );                  map.addLayers([image]);      Following the screen shot when above KML is displayed on my local system without rotation. ![enter image description here](http://i.stack.imgur.com/wHKN3.png) Google Maps does not seem to rotate the image, but Google Earth plugin does. ![enter image description here](http://i.stack.imgur.com/A0xJ5.png) ![enter image description here](http://i.stack.imgur.com/lAyqz.png)


title: Automatically apply topology rule to only selected portions of a .shp file?

text: I am looking for some examples/tutorials on how to apply a topology rule to only a selected portion of a shp file using a script in ArcGIS/GRASS. Using ArcGIS seems to always require an interactive approach - with the topology editor able to highlight exeption and the user can mark them as exeptions. GRASS seems to have much less detailed topology editing capabilities compared to ArcGIS. I would like to have a script that automatically fixes exceptions or not depending on something like an sql statement or the proximity of exceptions to another vector file. For example, I want to fix dangles in roads everywhere but downtown, because downtown they are purposeful but out of downtown they are mistakes.


title: Graduated Symbol classification not working with a decimal range

text: I am classifying an .adf wind resource map which I vectorised using the Polygonise tool (which worked very well). I have classified the data using the Graduated symbol renderer using Natural Breaks (Jenks) mode. Unfortunately although the Range column will take decimals, only columns with whole numbers will render. Does anyone know if there is a way to do this. Am using QGIS 1.7.3.


title: How to disable Layer Preview in geoserver, which is available by default for public?

text: I am using geoserver for publishing my spatial data. The problem I am facing is, Any public user can download my spatial data without my permission using default URL of geoserver. This is obviously undesirable in any context. Please provide guidelines.


title: How to write addressToLocations onError?

text: How can I write a code for checking the error condition in Javascript (ArcGIS) for addressToLocations(params, callback?, errback?). if somebody entered an invalid address- how can I call error function? Please share your thoughts! my code below:- -->               <meta http-equiv="X-UA-Compatible" content="IE=EmulateIE7" />     <title>Find Address</title>     <link rel="stylesheet" type="text/css" href="http://serverapi.arcgisonline.com/jsapi/arcgis/3.2/js/dojo/dijit/themes/claro/claro.css">     <link rel="stylesheet" type="text/css" href="http://serverapi.arcgisonline.com/jsapi/arcgis/3.2/js/esri/css/esri.css" />


title: How do you search for a particular attribute from the attributes table?

text: I am working with a postcode shapefile which has over 8000 attributes and wish to search for particular attributes without having to scroll through thousands of others, is there a way of searching and selecting particular attributes?


title: Most efficient way to tile/pyramid 4 million images using gdal for geoserver WMS

text: My scenario consist of 4 million geolocated images each measuring 1536x1536 pixels in size. What is the most efficient way to pyramid these images? My end goal is to serve these images using geoserver via WMS requests. I have tried creating a virtual dataset of all images and putting them into gdal_retile.py. I have also tried merging each tile with a transparent "null" image that covers the globe then pyramid the result. Both methods seem like they will take forever and make me think there is a better way. I also have the requirement of adding additional images in the future. Is there an easy way to do this without having to regenerate everything again? I am a geoserver and gdal noob so any any help is greatly appreciated!!


title: Openlayers wfs vector layer

text: I'm trying to add a vector layer with a wfs layer and it's not working.Can anyone point a solution? Openlayers code:               var wfs = new OpenLayers.Layer.Vector(                     "Stavros Features",                     {                         strategies: [new OpenLayers.Strategy.BBOX()]                         ,projection: new OpenLayers.Projection("EPSG:4326")                         ,protocol: new OpenLayers.Protocol.WFS({                             version: "1.1.0"                             ,srsName: "EPSG:4326"                             ,url: "http://localhost:8080/geoserver/wfs"                             ,featureNS :  "http://postgis.org"                             ,featureType: "polygon"                             ,geometryName: "the_geom"                         })                     });      In the geoserver i have a workspace _'myWorkspace'_ , a store _'myStore'_ and a layer _'myLayer'_ and wfs is enabled. Also in the _'publishing'_ tab of my layer in the _'WMS Settings Default Style'_ i have selected _'polygon'_ as myLayer is multipolygon. What am i doing wrong?   Am i suppose to specify the workspace name in the 'featureNS' property? Is there a good tutorial for open layers?   I have read the docs which are very poor and the examples doesn't help me always. I have also read the user manual of geoserver and i found it very good. Thanks.


title: How to create jpg image of map in Geomedia using vb.net?

text: I am using Geomedia. I want create jpg image of the map. How can it be done using vb.net?


title: How to export a raster into ascii without nodata values?

text: I need to export a grd raster file into ascii .xyz without nodata values. This is not a problem in itself...I´m using Saga-GIS export grid to xyz module.(ArcGIS 10.01 produces a very strange table of xyz values with lots of fields??)...but the person who wants to use the data would like it without the no-data values. The project area is a state border - so there are always going to be no-data values within the bounding box due to the unregular nature of the border. The data will be used in WINDPRO (which I do not have so I can´t test it).   1. is there a way to export a raster without getting no-data values?   2. Does anyone know windpro, and what could the problem be with no-data values? thanks for any suggestions, rob


title: are there any shapefiles for a simple road network of the UK that are similar to what you would find on google maps?

text: I wish to add a shapefile to ArcMap that displays a road network which is similar to what you may find on google maps. I dont require a detailed ordnance survey map, just the roads


title: Using an InfoWindow with a Feature Layer and a DialogWindowBase

text: I need to have an infoWindow that displays the Outfields of a FeatureLayer. The infoWindow is popping up, but it is popping up no matter where you click on the map and it does not contain any of the content. I have looked through the documentation and the only difference is that instead of using a UserControl, I am using a DialogWindowBase. Below is my XAML and C#. XAML:               <DataTemplate x:Key="MyFeatureLayerInfoWindowTemplate">                         <StackPanel Margin="7">                             <TextBlock Text="{Binding [WVAL_ID], StringFormat='ID: \{0\}'}" FontWeight="Bold" Foreground="White"  />                             <TextBlock Text="{Binding [SIZE], StringFormat='Size: \{0\}'}" FontWeight="Bold" Foreground="White" />                             <TextBlock Text="{Binding [WVAL_USE], StringFormat='Use: \{0\}'}" FontWeight="Bold" Foreground="White" />                             <TextBlock Text="{Binding [WVAL_TYPE], StringFormat='Type: \{0\}'}" FontWeight="Bold" Foreground="White"  />                             <TextBlock Text="{Binding [OPENS], StringFormat='Opens: \{0\}'}" FontWeight="Bold" Foreground="White"  />                             <TextBlock Text="{Binding [TURNS], StringFormat='Turns: \{0\}'}" FontWeight="Bold" Foreground="White"  />                             <StackPanel Orientation="Horizontal" >                                  </StackPanel>                         </StackPanel>                     </DataTemplate>          <esri:Map x:Name="MyMap" Background="#FFE3E3E3" WrapAround="True" MouseClick="MyMap_MouseClick" MouseRightButtonDown="MyMap_MouseRightButtonDown" Extent="-9834972.92753924,4441899.425293319,-9833977.88119163, 4442762.485358352">                    <esri:Map.Layers>                              <esri:ArcGISTiledMapServiceLayer ID="MyLayer" Url="http://services.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer" />                         <esri:ArcGISDynamicMapServiceLayer ID="Calvert_City" Url="http://localhost:6080/arcgis/rest/services/CalvertCity_Test_2/MapServer"/>           <esri:FeatureLayer ID="Water_Valves" Url="http://localhost:6080/arcgis/rest/services/CalvertCity_Test_2/MapServer/0" Renderer="{StaticResource MySimpleRenderer}" OutFields="WVAL_ID,SIZE,WVAL_USE,WVAL_TYPE,OPENS,TURNS" MouseLeftButtonDown="FeatureLayer_MouseLeftButtonDown">         </esri:FeatureLayer>                         </esri:Map.Layers>           </esri:Map>                      <esri:InfoWindow x:Name="MyInfoWindow"                                  Padding="2"                                  CornerRadius="20"                                  Background="{StaticResource PanelGradient}"                                  Map="{Binding ElementName=My_Map}"                                  ContentTemplate="{StaticResource MyFeatureLayerInfoWindowTemplate}"                              MouseLeftButtonUp="MyInfoWindow_MouseLeftButtonUp">                      </esri:InfoWindow>      C#: Here is the C# code that I have that will at least build and allow the infoWindow to appear:               private void MyMap_MouseClick(object sender, Map.MouseEventArgs e)             {                      FeatureLayer featurelayer = MyMap.Layers["Water_Valves"] as FeatureLayer;                 System.Windows.Point screenPoint = MyMap.MapToScreen(e.MapPoint);                      MyInfoWindow.Anchor = e.MapPoint;                 MyInfoWindow.IsOpen = true;             }      Here is the way the documentation does it:                private void MyMap_MouseClick(object sender, ESRI.ArcGIS.Client.Map.MouseEventArgs e)             {                 FeatureLayer featureLayer = MyMap.Layers["MyFeatureLayer"] as FeatureLayer;                 System.Windows.Point screenPnt = MyMap.MapToScreen(e.MapPoint);                      // Account for difference between Map and application origin                 GeneralTransform generalTransform = MyMap.TransformToVisual(Application.Current.MainWindow);                 System.Windows.Point transformScreenPnt = generalTransform.Transform(screenPnt);                      IEnumerable<Graphic> selected =                     featureLayer.FindGraphicsInHostCoordinates(transformScreenPnt);                      foreach (Graphic g in selected)                 {                          MyInfoWindow.Anchor = e.MapPoint;                     MyInfoWindow.IsOpen = true;                     //Since a ContentTemplate is defined, Content will define the DataContext for the ContentTemplate                     MyInfoWindow.Content = g.Attributes;                     return;                 }             }      When I try to run it while using the generalTransform is says that is cannot do it because it is not an ancestor. Any help will be appreciated. Thanks!


title: How to use python to export ArcGIS Server Map Service Properties?

text: What arcpy module can I use to write a python script to mine for fields of my map service properties and export them to a .tx file? I assume we need to tap into GIS Servers within the ArcCatalog interface?


title: OpenGeo SDK: Add New Server trouble with proxy/?url=

text: I'm currently developing an application with the OpenGeo SDK. I wish to have the "Add New Server" functionality and I've followed the example from the documentation: http://suite.opengeo.org/opengeo- docs/webapps/sdkintro/viewer/addnewserver.html In my app.'s file, I've added proxy: "proxy/?url=" and also * @require overrides/override-ext-ajax.js When I test the application on my server (e.g., suite-sdk debug path/to/myapp) everything works fine. However, when I deploy the application (I use a Tomcat container) then the "Add New Server" functionalities doesn't work any more. I've always the following error message when I try to add a WMS: Error getting WMS capabilities (Layer source not available.). Please check the url and try again. I don't know what I'm doing wrong. Thanks in advance for your help :-)


title: why my tile rendered and Httpd(mod_tile) server not communicating?

text: it's been weeks since i have been trying to build my own openstreet maps tile server after following dozens tutorials i have been able to install renderd and mod_tile and i have successfully uploaded data in postgresql-server. Now i am stuck at this problem my mod_tile is running and my renderd is also running but both are not communicating? if i run generate_tiles.py it will generate tiles successfully but everytime i request tiles by hhtp://myserver/osm_tile/0/0/0.png it reposnd wiht 404. My renderd.conf file               [renderd]     socketname=/var/run/renderd/renderd.sock     num_threads=4     tile_dir=/var/lib/mod_tile     stats_file=/var/run/renderd/renderd.stats               [mapnik]     plugins_dir=/usr/lib64/mapnik/input     font_dir=/usr/lib64/mapnik/fonts     font_dir_recurse=1          [default]     URI=/osm_tiles/     TILEDIR=/var/lib/mod_tile     XML=/home/gis/mapnik/osm.xml     HOST=localhost     TILESIZE=256      and here is mod_tile.conf file                   # This is the Apache server configuration file for providing OSM tile support     # through mod_tile          LoadModule tile_module /usr/lib64/httpd/modules/mod_tile.so          <VirtualHost *:80>         ServerName localhost         ServerAlias a.tile.openstreetmap.org b.tile.openstreetmap.org c.tile.openstreetmap.org d.tile.openstreetmap.org         DocumentRoot /var/www/html          ModTileTileDir /var/lib/mod_tile          LoadTileConfigFile /usr/local/etc/renderd.conf          ModTileEnableStats On          ModTileBulkMode Off          ModTileRequestTimeout 3          ModTileMissingRequestTimeout 10          ModTileMaxLoadOld 16          ModTileMaxLoadMissing 50          ModTileVeryOldThreshold 31536000000000          ModTileRenderdSocketName /var/run/renderd/renderd.sock          #   ModTileRenderdSocketAddr renderd.mydomain.com 7653               #ModTileCacheExtendedHostname cache.tile.openstreetmap.org     #ModTileCacheExtendedDuration 2592000          ModTileCacheDurationMax 604800          ModTileCacheDurationDirty 900          ModTileCacheDurationMinimum 10800          ModTileCacheDurationMediumZoom 13 86400          ModTileCacheDurationLowZoom 9 518400          ModTileCacheLastModifiedFactor 0.20               ModTileEnableTileThrottling Off     ModTileEnableTileThrottlingXForward 0     ModTileThrottlingTiles 10000 1           ModTileThrottlingRenders 128 0.2              LogLevel debug      This is output when i run my renderd                /usr/local/bin/renderd -f -c /usr/local/etc/renderd.conf     renderd[29317]: Rendering daemon started     renderd[29317]: Initiating reqyest_queue     iniparser: syntax error in /usr/local/etc/renderd.conf (7):     -> ;[renderd01]     iniparser: syntax error in /usr/local/etc/renderd.conf (14):     -> ;[renderd02]     iniparser: syntax error in /usr/local/etc/renderd.conf (33):     -> ;** config options used by mod_tile, but not renderd **     iniparser: syntax error in /usr/local/etc/renderd.conf (44):     -> ;[style2]     iniparser: syntax error in /usr/local/etc/renderd.conf (51):     -> ;** config options used by mod_tile, but not renderd **     renderd[29317]: Parsing section renderd     renderd[29317]: Parsing render section 0     renderd[29317]: Parsing section mapnik     renderd[29317]: Parsing section default     renderd[29317]: config renderd: unix socketname=/var/run/renderd/renderd.sock     renderd[29317]: config renderd: num_threads=4     renderd[29317]: config renderd: num_slaves=0     renderd[29317]: config renderd: tile_dir=/var/lib/mod_tile     renderd[29317]: config renderd: stats_file=/var/run/renderd/renderd.stats     renderd[29317]: config mapnik:  plugins_dir=/usr/lib64/mapnik/input     renderd[29317]: config mapnik:  font_dir=/usr/lib64/mapnik/fonts     renderd[29317]: config mapnik:  font_dir_recurse=1     renderd[29317]: config renderd(0): Active     renderd[29317]: config renderd(0): unix socketname=/var/run/renderd/renderd.sock     renderd[29317]: config renderd(0): num_threads=4     renderd[29317]: config renderd(0): tile_dir=/var/lib/mod_tile     renderd[29317]: config renderd(0): stats_file=/var/run/renderd/renderd.stats     renderd[29317]: config map 0:   name(default) file(/home/gis/mapnik/osm.xml) uri(/osm_tiles/) htcp() host(localhost)     renderd[29317]: Initialising unix server socket on /var/run/renderd/renderd.sock     renderd[29317]: Created server socket 4     renderd[29317]: Renderd is using mapnik version 2.2.0     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSans.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSerif-Italic.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSans-Oblique.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSans-ExtraLight.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSans-BoldOblique.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSerifCondensed-Bold.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSansCondensed-Bold.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSerif-BoldItalic.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSerifCondensed.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSansMono-BoldOblique.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/unifont-5.1.20080907.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSansMono.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSansMono-Bold.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSerif-Bold.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSerif.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSansCondensed-BoldOblique.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSerifCondensed-BoldItalic.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSans-Bold.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSansCondensed-Oblique.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSerifCondensed-Italic.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSansMono-Oblique.ttf     renderd[29317]: DEBUG: Loading font: /usr/lib64/mapnik/fonts/DejaVuSansCondensed.ttf     Running in foreground mode...     debug: init_storage_backend: initialising file storage backend at: /var/lib/mod_tile     debug: init_storage_backend: initialising file storage backend at: /var/lib/mod_tile     renderd[29317]: Starting stats thread     debug: init_storage_backend: initialising file storage backend at: /var/lib/mod_tile     renderd[29317]: Loading parameterization function for      renderd[29317]: Loading parameterization function for      debug: init_storage_backend: initialising file storage backend at: /var/lib/mod_tile     renderd[29317]: Loading parameterization function for      renderd[29317]: Loading parameterization function for      renderd[29317]: Using web mercator projection settings     renderd[29317]: Using web mercator projection settings     renderd[29317]: Using web mercator projection settings     renderd[29317]: Using web mercator projection settings      Help me out here and point any mistake in configuration or anything?


title: How to parse serial date in ArcGIS?

text: I have imported a shapefile that contains a serial date number that looks like this: 41256.260799. I would like to create a new field that shows only the hour and second like this: 06:33. I've tried the available VBS and Python date functions, but nothing seems to work. Would I make my new field type text? Any help appreciated.


title: Are there any GIS software that can generate the input file for Air Quality Model(e.g. Caline4)?

text: I want to find a GIS software that can be used to graphically input some data, like geometry of emission points, receiver points and roads, by clicking in the base map (.dxf). And then I can reformat the data, saving them into a text file with a format which can be used by other air quality models (e.g. Caline4, ISCST3). Please give me some clue to achieve it. I would only need to know the name of the software and the documentation of reformatting data.


title: Can't install CartoDB on Ubuntu Lucid - mapnik dependency issue

text: Previously asked this question at askubuntu - but maybe here is the right place. I'm trying to install CartoDB, and have run into the dreaded "you may have requested an impossible situation" message:               $ sudo apt-get install -y libmapnik-dev python-mapnik mapnik-utils     Reading package lists... Done     Building dependency tree            Reading state information... Done     Some packages could not be installed. This may mean that you have     requested an impossible situation or if you are using the unstable     distribution that some required packages have not yet been created     or been moved out of Incoming.     The following information may help to resolve the situation:          The following packages have unmet dependencies:       libmapnik-dev: Depends: libmapnik (= 2.2.0+dev20131008.git.fe06c95-1~lucid1) but it is not going to be installed                      Depends: libboost-thread-dev (>= 1.42) but it is not going to be installed                      Depends: libboost-filesystem-dev (>= 1.42) but it is not going to be installed                      Depends: libboost-regex-dev (>= 1.42) but it is not going to be installed                      Depends: libboost-python-dev (>= 1.42) but it is not going to be installed                      Depends: libboost-system-dev (>= 1.42) but it is not going to be installed                      Depends: libboost-program-options-dev (>= 1.42) but it is not going to be installed       mapnik-utils: Depends: libboost-program-options1.49.0 (>= 1.49.0-1) but it is not installable                     Depends: libboost-system1.49.0 (>= 1.49.0-1) but it is not installable                     Depends: libmapnik (= 2.2.0+dev20131008.git.fe06c95-1~lucid1) but it is not going to be installed       python-mapnik: Depends: libboost-python1.49.0 (>= 1.49.0-1) but it is not installable                      Depends: libboost-thread1.49.0 (>= 1.49.0-1) but it is not installable                      Depends: libmapnik (= 2.2.0+dev20131008.git.fe06c95-1~lucid1) but it is not going to be installed     E: Broken packages      I've tried[1] most of the suggestions in the How do I resolve unmet dependencies question, but no dice. There are a few things that make this situation particularly hard to diagnose and fix (for me):   * It's Lucid (and it has to be, for CartoDB to work)   * I'm just following instructions (a pretty long list of them)   * The instructions didn't quite work, so I've had to vary a couple of things, such as installing Python 2.7.3 (using altinstall, so it doesn't overwrite 2.6.x)   * There are lots of custom ppa's. The commands I ran are documented in this Gist: https://gist.github.com/stevage/9593847 And the full (huge) log is here: https://gist.github.com/stevage/9593835 As requested:               $ cat /etc/apt/sources.list     ## Note, this file is written by cloud-init on first boot of an instance     ## modifications made here will not survive a re-bundle.     ## if you wish to make changes you can:     ## a.) add 'apt_preserve_sources_list: true' to /etc/cloud/cloud.cfg     ##     or do the same in user-data     ## b.) add sources in /etc/apt/sources.list.d     ## c.) make changes to template file /etc/cloud/templates/sources.list.tmpl     #          # See http://help.ubuntu.com/community/UpgradeNotes for how to upgrade to     # newer versions of the distribution.     deb http://archive.ubuntu.com/ubuntu/ lucid main     deb-src http://archive.ubuntu.com/ubuntu/ lucid main          ## Major bug fix updates produced after the final release of the     ## distribution.     deb http://archive.ubuntu.com/ubuntu/ lucid-updates main     deb-src http://archive.ubuntu.com/ubuntu/ lucid-updates main          ## N.B. software from this repository is ENTIRELY UNSUPPORTED by the Ubuntu     ## team. Also, please note that software in universe WILL NOT receive any     ## review or updates from the Ubuntu security team.     deb http://archive.ubuntu.com/ubuntu/ lucid universe     deb-src http://archive.ubuntu.com/ubuntu/ lucid universe     deb http://archive.ubuntu.com/ubuntu/ lucid-updates universe     deb-src http://archive.ubuntu.com/ubuntu/ lucid-updates universe          ## N.B. software from this repository is ENTIRELY UNSUPPORTED by the Ubuntu      ## team, and may not be under a free licence. Please satisfy yourself as to     ## your rights to use the software. Also, please note that software in      ## multiverse WILL NOT receive any review or updates from the Ubuntu     ## security team.     # deb http://archive.ubuntu.com/ubuntu/ lucid multiverse     # deb-src http://archive.ubuntu.com/ubuntu/ lucid multiverse     # deb http://archive.ubuntu.com/ubuntu/ lucid-updates multiverse     # deb-src http://archive.ubuntu.com/ubuntu/ lucid-updates multiverse          ## Uncomment the following two lines to add software from the 'backports'     ## repository.     ## N.B. software from this repository may not have been tested as     ## extensively as that contained in the main release, although it includes     ## newer versions of some applications which may provide useful features.     ## Also, please note that software in backports WILL NOT receive any review     ## or updates from the Ubuntu security team.     # deb http://archive.ubuntu.com/ubuntu/ lucid-backports main restricted universe multiverse     # deb-src http://archive.ubuntu.com/ubuntu/ lucid-backports main restricted universe multiverse          ## Uncomment the following two lines to add software from Canonical's     ## 'partner' repository.     ## This software is not part of Ubuntu, but is offered by Canonical and the     ## respective vendors as a service to Ubuntu users.     # deb http://archive.canonical.com/ubuntu lucid partner     # deb-src http://archive.canonical.com/ubuntu lucid partner          deb http://security.ubuntu.com/ubuntu lucid-security main     deb-src http://security.ubuntu.com/ubuntu lucid-security main     deb http://security.ubuntu.com/ubuntu lucid-security universe     deb-src http://security.ubuntu.com/ubuntu lucid-security universe     # deb http://security.ubuntu.com/ubuntu lucid-security multiverse     # deb-src http://security.ubuntu.com/ubuntu lucid-security multiverse     deb http://ppa.launchpad.net/fkrull/deadsnakes/ubuntu lucid main     deb-src http://ppa.launchpad.net/fkrull/deadsnakes/ubuntu lucid main      and               $ cat /etc/apt/sources.list     ## Note, this file is written by cloud-init on first boot of an instance     ## modifications made here will not survive a re-bundle.     ## if you wish to make changes you can:     ## a.) add 'apt_preserve_sources_list: true' to /etc/cloud/cloud.cfg     ##     or do the same in user-data     ## b.) add sources in /etc/apt/sources.list.d     ## c.) make changes to template file /etc/cloud/templates/sources.list.tmpl     #          # See http://help.ubuntu.com/community/UpgradeNotes for how to upgrade to     # newer versions of the distribution.     deb http://archive.ubuntu.com/ubuntu/ lucid main     deb-src http://archive.ubuntu.com/ubuntu/ lucid main          ## Major bug fix updates produced after the final release of the     ## distribution.     deb http://archive.ubuntu.com/ubuntu/ lucid-updates main     deb-src http://archive.ubuntu.com/ubuntu/ lucid-updates main          ## N.B. software from this repository is ENTIRELY UNSUPPORTED by the Ubuntu     ## team. Also, please note that software in universe WILL NOT receive any     ## review or updates from the Ubuntu security team.     deb http://archive.ubuntu.com/ubuntu/ lucid universe     deb-src http://archive.ubuntu.com/ubuntu/ lucid universe     deb http://archive.ubuntu.com/ubuntu/ lucid-updates universe     deb-src http://archive.ubuntu.com/ubuntu/ lucid-updates universe          ## N.B. software from this repository is ENTIRELY UNSUPPORTED by the Ubuntu      ## team, and may not be under a free licence. Please satisfy yourself as to     ## your rights to use the software. Also, please note that software in      ## multiverse WILL NOT receive any review or updates from the Ubuntu     ## security team.     # deb http://archive.ubuntu.com/ubuntu/ lucid multiverse     # deb-src http://archive.ubuntu.com/ubuntu/ lucid multiverse     # deb http://archive.ubuntu.com/ubuntu/ lucid-updates multiverse     # deb-src http://archive.ubuntu.com/ubuntu/ lucid-updates multiverse          ## Uncomment the following two lines to add software from the 'backports'     ## repository.     ## N.B. software from this repository may not have been tested as     ## extensively as that contained in the main release, although it includes     ## newer versions of some applications which may provide useful features.     ## Also, please note that software in backports WILL NOT receive any review     ## or updates from the Ubuntu security team.     # deb http://archive.ubuntu.com/ubuntu/ lucid-backports main restricted universe multiverse     # deb-src http://archive.ubuntu.com/ubuntu/ lucid-backports main restricted universe multiverse          ## Uncomment the following two lines to add software from Canonical's     ## 'partner' repository.     ## This software is not part of Ubuntu, but is offered by Canonical and the     ## respective vendors as a service to Ubuntu users.     # deb http://archive.canonical.com/ubuntu lucid partner     # deb-src http://archive.canonical.com/ubuntu lucid partner          deb http://security.ubuntu.com/ubuntu lucid-security main     deb-src http://security.ubuntu.com/ubuntu lucid-security main     deb http://security.ubuntu.com/ubuntu lucid-security universe     deb-src http://security.ubuntu.com/ubuntu lucid-security universe     # deb http://security.ubuntu.com/ubuntu lucid-security multiverse     # deb-src http://security.ubuntu.com/ubuntu lucid-security multiverse     deb http://ppa.launchpad.net/fkrull/deadsnakes/ubuntu lucid main     deb-src http://ppa.launchpad.net/fkrull/deadsnakes/ubuntu lucid main     ubuntu@tm-21:/mnt/cartodb$ cat /etc/apt/sources.list.d/*     deb http://ppa.launchpad.net/cartodb/gis/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/gis/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/mapnik/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/mapnik/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/nodejs/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/nodejs/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/postgresql/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/postgresql/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/redis/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/redis/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/varnish/ubuntu lucid main     deb http://ppa.launchpad.net/cartodb/varnish/ubuntu lucid main      Any suggestions, advice how to proceed? I really don't grok the APT system well, unfortunately. It's a server, so I can't benefit from many of the GUI based tips. Thanks for any help! (I've also posted to the CartoDB mailing list) [1]                  46  sudo apt-get install -yV libmapnik-dev python-mapnik mapnik-utils        47  sudo apt-get install -f libmapnik-dev python-mapnik mapnik-utils        48  sudo apt-get clean        49  sudo apt-get install -f libmapnik-dev python-mapnik mapnik-utils        50  sudo apt-get autoclean        51  sudo apt-get install -f libmapnik-dev python-mapnik mapnik-utils        52  sudo apt-get install -f        53  sudo apt-get install -f libmapnik-dev python-mapnik mapnik-utils        54  sudo dpkg --configure -a        55  sudo apt-get -f install        56  sudo apt-get install -f libmapnik-dev python-mapnik mapnik-utils        57  sudo apt-get -u dist-upgrade        58  sudo apt-get -o Debug::pkgProblemResolver=yes dist-upgrade        59  sudo apt-get install -y libmapnik-dev python-mapnik mapnik-utils


title: Postgis performance: st_dwithin vs static buffer and st_contains

text: I have a table with a polyline column. I need to find which of those polylines are within 20 meters from a bunch of moving points. I can do this no problem with st_dwithin(). What I would like to know is if there would be a performance benefit if I were to create a static buffer around the polylines, like a materialized view, and use st_contains to test the points.


title: City Engine - 3D zoning envelope & street orientation

text: Is there a way to change the front & back angles that generate the 3D zoning envelope according to street orientation? Ideally there should be 6 cases which describe angular ranges and for each case a different angle is assigned to front,back and sides of the lot. Would this require scripting? I am a complete novice in this and any help would be appreciated.


title: How to get point x and y geometry from a shapefile?

text: I have a shapefile of US Cities and I want to get the X and Y coordinate of EACH city in that shapefile. I've tried this:               for city in city_cursor:         geom = city.Shape         point = geom.getPart()         citylist_City_Name.append(city.CITY_NAME)         citylist_Country.append(city.CNTRY_NAME)         citylist_Admin.append(city.ADMIN_NAME)         citylist_Population.append(city.Population)         citylist_X_Coor.append(geom.point.X)         citylist_Y_Coor.append(geom.point.Y)      but I get the error: Traceback (most recent call last): File "C:/Users/workd.py", line 43, in citylist_X_Coor.append(geom.point.X) AttributeError: 'PointGeometry' object has no attribute 'point' I don't really understand the error message? How can I fix it? Thanks!


title: Using arcpy, is it possible to check if a thumbnail exists for a specified file?

text: I have written a very simple tool that builds thumbnails for every map document in a folder : `root_fld` stands for root folder and `rec` tells if search must be recursive.               import arcpy     import os     import fnmatch     import glob          # getting list of mxds     if rec:         # going recursive         matches = []         for root, dirnames, filenames in os.walk(root_fld):             for filename in fnmatch.filter(filenames, '*.mxd'):                 matches.append(os.path.join(root, filename))     else:         # going flat         matches = glob.glob(os.path.join(root_fld, '*.mxd'))          # create thmb     for m in matches:         mxd = arcpy.mapping.MapDocument(m)         mxd.makeThumbnail()         mxd.save()      Works fine but some thumbnails are a pain to create and i'd like my script to filter and ignore mxds that already have a thumbnail :               ...     if not mxd.hasThumbnail():         mxd.makeThumbnail()         mxd.save()      According to MapDocument class documentation (http://resources.arcgis.com/fr/help/main/10.1/index.html#//00s30000000n000000), this property doesn't seem to exist. But maybe there's a workaround ?


title: Using OpenLayers how can I Extract some information relative to the whole KML file?

text: I am creating a KML layer like this:               new OpenLayers.Layer.Vector("KML", {                     projection: myDisplayProjection,                     strategies: [new OpenLayers.Strategy.Fixed()],                     protocol: new OpenLayers.Protocol.HTTP({                         url: myKmlUrl,                         format: new OpenLayers.Format.KML({                             extractStyles: true,                             extractAttributes: true                         })                     })                 });      I am extracting styles and attributes. I have some attributes attached to the document (as well as the placemarks/features). Will this information be on the resultant layer (I think not) or how do I get this information other wise? The information relevant to the whole layer. Thanks


title: Defining the appropriate Projected Coordinate System

text: I am a beginner at ArcGIS and have a technical question. I need to assign a projected coordinate system to my boundary/polygon data in order to calculate square miles. There are many to choose from. Which should I be using? Here are my geographic coordinate system details: Angular Unit: Degree (0.017453292519943295) Prime Meridian: Greenwich (0.000000000000000000) Datum: D_North_American_1983 Spheroid: GRS_1980 Semimajor Axis: 6378137.000000000000000000 Semiminor Axis: 6356752.314140356100000000 Inverse Flattening: 298.257222101000020000 Thanks in advance for any consideration.


title: Converting LinestringM to Polygon PostGIS

text: I'm working with PostGIS and plotting polygons on a Google Earth plugin. I've found a few examples via Google on how to convert a LinestringM to a Polygon for rendering, but none of them, including the PostGIS docs, define very well the params for a LinestringM. Can someone identify what's happening here? Specifically, I need to know what each of the four params of LinestringM are supposed to represent. Snippet:               ....     parameters.put("wkb_geometry", "ST_MakePolygon(ST_GeomFromText('LINESTRINGM(75.15 29.53,77 29,77.6 29.5, 75.15 29.53)'))");     ....


title: Calculate bearing between two decimal GPS coordinates

text: I'm building a virtual tour application that uses custom panoramas in conjunction with Google's StreetView service. With StreetView, you can have 1 to n links at each particular location to move to, and these links are based on compass bearings. Each panorama has an associated GPS location to six decimal places acquired via iTouchMaps. My issue is when I calculate the compass bearings from each link location, it returns a bearing that is a tiny (as in 0.00001 degree difference) from true north. Given a map like                                      4                            |                            |                            |                            2                            |                            |                            |                            1      where each number is a location and has the following values:               1- 43.682213, -70.450696     2- 43.682194, -70.450769     4- 43.682179, -70.450837      What I'd like is a method to determine programatically that the bearing from 1 to 2 and 2 to 4 is 70 degrees, and 2 to 1 is 250 degrees. Using the code from here as a starting point, I wrote the following function               //starting lat/long along with converting degrees to radius         var endLat = location.lat();  //for rhumb calc         var endLong = location.lng();              //loop over response, calculate new headings for links and add link to array         for(var i=0; i<data.length; i++){           //this link's lat/long coordinates           var startLat = data[i].lat;           var startLong = data[i].lon;                //get the delta values between start and end coordinates           var dLong = endLong - startLong;                //calculate rhumb bearing to           var dPhi = Math.log(Math.tan(endLat/2+Math.PI/4)/Math.tan(startLat/2+Math.PI/4));           if (Math.abs(dLong) > Math.PI ) dLong = dLong > 0 ? -(2*Math.PI-dLong) : (2*Math.PI+dLong);           var bearing = (toDeg(Math.atan(dLong, dPhi)) + 360) % 360;      Using the above function and substituting               location for position 2     data[0] for position 1     data[1] for position 4      I get the bearings               359.9958174081038 degrees to 4 (should be 70 degrees)      and               0.0038961130024404156 degrees to 1 (should be 250)      Why is this function returning such different values than expected?


title: Insert Geography data into SQL Server 2008

text: Can someone guide me to or give me an example on how would one go about insert/query Geography data in Sql Server 2008 ( Polygon with 4 points and Point in specific) using C# EF or Linq2Sql. I've been searching for few days now and i can't figure it out :( Thanks in advance!


title: How to increase developer productivity with ESRI technology

text: We are a small team of .NET developers. We have ample GIS experience, and none of us are new to software/database development or system administration. We have technical degrees and many years of industry experience. We have attended ESRI developer summits. ESRI technology - primarily ArcGIS Server, ArcSDE, and ArcObjects - plays a small but necessary role in all the software we develop. Despite ESRI's minority status in our technology stack, we spend an inordinate amount of time troubleshooting elusive bugs, crafting workarounds, deciphering its vague error messages, tracking down performance problems, and recycling processes. Typically our problems stem from genuine bugs, poor exception handling, limiting design/architectural decisions, lack of documentation, instability, or some combination thereof. (I'm speaking of the ESRI stack here.) From a project manager's perspective, I'm very concerned about team productivity. This costs us a lot of time. We do not have time to learn every idiosyncrasy of the ESRI stack, but we still need to get things done. (Can't live with it, can't live without it.) _What pragmatic suggestions do you have for increasing developer productivity with ESRI in the mix?_ Edit: I'm not looking for suggestions on alternative technology stacks.


title: QGIS Zoom from python console

text: What are the commands to apply all types of zoom (full, on selected layer, etc.) in python console? I'm searching both on pyqgis cookbook and API documentantion, but I'm not finding anything interesting.


title: Gathering GIS data from the masses

text: I am looking at collecting spatial information from a range of stakeholders to identify opportunities within a landscape. However, the data that we are collecting are not presently mapped and so we wish to use an interactive method to get this information onto a map. Therefore, my question is: Which online websites/methods would people recommend for circulating to a non- technical audience to try and get them to map what they know?? Where the emphasis is on ease of use by the stakeholder. Many thanks in advance for any advice anybody can offer!


title: How to change the base map in Ushahidi?

text: Fascinated by the idea of a 'crowd map' and 'citizen journalism', I've been perusing Ushahidi and it should serve me well. I would like to emulate it except for - the base map. I would prefer a personal base map(s) served from Geoserver for instance and not the Google Map. I see OSM can also be used but it lacks the detail I require. Any ideas on How-To or alternative projects? [I've literally scrounged the web via Google search for answers without success. The Ushahidi forum has not been helpful as well]


title: Why does GeoNames say there are only 73 parks in Mexico?

text: I did this GeoNames query http://api.geonames.org/search?featureCode=PRK&username=genadinik&country=MX&style=full&maxRows=10 and it says there are only 73 parks in all of Mexico. Am I doing something wrong, or GeoNames mostly covers the United States? Thanks!


title: How to change the default color used by openlayers?

text: I am very, very new to web mapping.. I have an OpenLayers map made from shape file from PostGIS database, have a vector layer and it has default style, but I want to change the style of vector layer. I was searching online, and I did not find simple and clear examples for this... Is there anyone who can help me, or maybe someone who wanted to do the same and found simple tutorial for this? Any help greatly appreciated.. Thanks in advance!!               var map;                  function init()             {                      var bounds = new OpenLayers.Bounds                 (                         68.089442, 6.752729,                         97.407576, 37.072537                     );                     var options = {                     controls: [ new OpenLayers.Control.Navigation(),                                 new OpenLayers.Control.PanZoom()                         ],                         maxExtent: bounds,                         maxResolution: 0.11843675,                         projection: "EPSG:4326",                         units: 'degrees'                     };                            map = new OpenLayers.Map ("map", options );                      var india = new OpenLayers.Layer.WMS(                     "cite:india_state - Tiled", "http://localhost:8080/geoserver/cite/wms",                     {                         LAYERS: 'cite:india_state',                         STYLES: 'style',                         //format: format,                         tiled: true,                         tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                     },                     {                         buffer: 0,                         displayOutsideMaxExtent: true,                         isBaseLayer: true,                         yx : {'EPSG:4326' : true}                     }                 );               var style = new OpenLayers.Style();     //rule used for all polygons     var rule_fsa = new OpenLayers.Rule({     symbolizer: {          fillColor: "#ff9a9a",     fillOpacity: 0.5,     strokeColor: "#000000",     strokeWidth: 1,     strokeDashstyle: "dash",     label: "${name}",     labelAlign: "cc",     fontColor: "#333333",     fontOpacity: 0.9,     fontFamily: "Arial",     fontSize: 14}     });     var rule_highlight = new OpenLayers.Rule({     filter: new OpenLayers.Filter.Comparison({     type: OpenLayers.Filter.Comparison.EQUAL_TO,     property: "classification",     value: "1",     }),     symbolizer: {     fillColor: "#FF7144",     fillOpacity: 0.6,     strokeColor: "#FF0000",     strokeWidth: 2,     strokeDashstyle: "solid",     label: " ${name}",     labelAlign: "cc",     fontColor: "#000000",     fontOpacity: 1,     fontFamily: "Arial",     fontSize: 16,     fontWeight: "600"}     });     style.addRules([rule_fsa, rule_highlight]);          var polygon = new OpenLayers.Layer.Vector("Polygon", {                 Style: 'style',                 rendererOptions: {zIndexing: true}             });                      map.addLayers([india,polygon]);                       map.zoomToMaxExtent();                    }


title: How to display the map.InfoWindow outside of the map container?

text: I need to show an infoWindow outside the map (i.e. have it owned by the body of the HTML document, not the map div). If i subclass InfoWindowBase I can do this myself but I really don't want to do that. Is there any way of doing this? I have tried the following without success:-                    var infoWindow = new esri.dijit.InfoWindow({           domNode: dojo.create("div", null, dojo.body())         });     map.InfoWindow = infoWindow;      and map.infoWindow.place("", dojo.body()); thanks, Phil


title: OpenLayers WFS protocol: work with filters if the protocol is defined out of layer constructor

text: follow question:wfs-load-feature-check the scenario is: i am trying to load a WFS layer, since it takes a while to load all features, I want to give user a notice while waiting for response. the process employs spatial filter (BBOX strategy) and attribute filters. I am wondering whether there is a good way to make WFS protocol work with these filters, in other words, only request the filtered features. thanks a lot! **code sample:**                var wfsLayer = new OpenLayers.Layer.Vector("WFS", {             strategies: [new OpenLayers.Strategy.BBOX()],             styleMap: new OpenLayers.StyleMap({                 strokeWidth: 3,                 strokeColor: "#333333"             })          });          var filt = new OpenLayers.Filter.Logical({         type: OpenLayers.Filter.Logical.OR,         filters: [             new OpenLayers.Filter.Comparison({                 type: OpenLayers.Filter.Comparison.EQUAL_TO,                 property: "TYPE",                 value: "highway"             }),             new OpenLayers.Filter.Comparison({                 type: OpenLayers.Filter.Comparison.EQUAL_TO,                 property: "TYPE",                 value: "road"             })         ]     });          var prot =  new OpenLayers.Protocol.WFS({         url:  "http://demo.opengeo.org/geoserver/wfs",         featureType: "tasmania_roads",         featureNS: "http://www.openplans.org/topp"     });          var _CallBack = function(resp) {         if(resp.error) {             console.log('error');             return -1;         }         wfsLayer.addFeatures(resp.features);         console.log('success');         return 1;     };          var response = prot.read({         callback: _CallBack     });


title: TIN Classes to a layer ARCMAP 10

text: Just wondering if there is a method of creating a layer for each classification range in a TIN. Another way to state my problem: I would like to make each colour code a layer that I can "turn on and off" i.e. it will have a check box at the front. Thanks David


title: How do I change a raster's geotransform with GDAL?

text: I'm trying to use GDAL to set the geotransform of a new raster that I've created by following the steps outlined in the GDAL API Tutorial.               # create the new dataset     driver = gdal.GetDriverByName('GTiff')     dataset = driver.Create('test_gt.tif', 60, 60, 1, gdal.GDT_Float32)          # check the default geotransform     print dataset.GetGeoTransform()  # prints (0, 1, 0, 0, 0, 1)          # try to alter the geotransform and ensure that it has been set     dataset.SetGeoTransform([0,1,0,0,0,-1])     print dataset.GetGeoTransform()  # prints (0, 1, 0, 0, 0, -1)          dataset = None  # closes the dataset          # Try reopening the dataset now and see if the geotransform has been set.     ds = gdal.Open('test_gt.tif')     print ds.GetGeoTransform()  #prints (0, 1, 0, 0, 0, 1)      The dataset.SetGeoTransform() documentation says that this should set the affine transformation coefficients (which, according to the `dataset.GetGeoTransform()` documentation is set to `[0, 1, 0, 0, 0, 1]` by default), but as you can see from my above code, the changes don't seem to actually take effect when I try to change them. I've even tried flushing the new raster's cache to disk by `dataset.FlushCache()`, but this doesn't seem to save the changes either. How can I have GDAL actually save the altered geotransform to disk? I'm using GDAL version 1.6.3 installed from PyPI (via the command-line tool `pip`).


title: How to combine corner-touching polygons into a single feature

text: I have created a polygon layer from rasters, and only have a single classification ID (and calculated acres) for each polygon. I would like to "dissolve" all touching features into larger, individual polygons. The 3 areas circled in red (below) contain green polygons/pixels that only touch on the points of the squares. Intuitively (from an end user perspective) all of the green polygons within a circled area should be connected and treated as a single polygon because there are no "breaks" where the features/cells do not touch. What would be the best method to accomplish "joining" or "dissolving" all of these features using ArcGIS. I have ArcInfo and Spatial Analyst extension. If there is a better methodology for doing this in Spatial Analyst, I can pre- process the polygons that way, too. ![Raster to Polygon](http://i.stack.imgur.com/cfKoe.jpg)


title: Openlayers using Google Maps as base layer: possible to zoom out farther than max extent?

text: See attached picture. This is how it is when the page loads. I'm currently using `map.zoomToMaxExtent();` ![enter image description here](http://i.stack.imgur.com/dLExg.png) Everything is working pretty well **but I'd like to zoom out more/see the whole world using OSM/Google Map Layers: is this even possible?** Note that the zoom bar on the left is set to the lowest bar when the page loads (if you can't see the map well, just right click the image and paste into a new browser tab). This post was somewhat helpful but even when I do `MIN_ZOOM_LEVEL` that doesn't do much. Also note I am not setting any bounds, extents, or resolutions. Toying with those did not seem to do the trick either. Any input is appreciated. Edit: here's a codeshare if you want to take a look: http://codeshare.io/GLC4i **UPDATE:** http://i.imgur.com/mSpTchY.png This is not perfected yet (and it makes me think something may be wrong with my OL code itself) but currently I am using this `var extent = new OpenLayers.Bounds(-1.32,81.71,-1.18,-61.80).transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:3857"));` I tried `var extent = new OpenLayers.Bounds(-20037508.34279,-10018754.17139,20037508.34279,10018754.17139);` (both with and without transform) but this did not do anything. **UPDATE 2** : http://i.imgur.com/XzuGxrn.png `wrapDateLine: false` works fine for now: can zoom out farther, pan around, etc etc. Now I would like to cut off everything to the left of the date line on the left side of the map (the 'first' portion), and cut off everything to the right of the date line on the right side of the map (the '3rd' portion) WHILE still keeping the map full screen. I've tried bounds, resolutions, but no go. Is this even feasible? Something like this but it should take up the full screen: possible to do? See below: http://i.imgur.com/feBEnx3.png


title: How to find the polygon containing a given coordinate using ArcGIS Geometry Service?

text: Say I have an ArcGIS Server GeometryService running, and a map service containing ZIP code polygons (all ArcGIS 10.1). I want to pass a coordinate, and have it return which ZIP code contains that coordinate. There seem to be many ways to approach "does THIS point fall within THIS polygon?", but how would I find out _which_ polygon among them all contains a point? In other words, there's this example where, given two points and a polygon, it will return which of the points is contained by the polygon. But what if you only know the point, and the objective is to determine which polygon?


title: How can I get timestamps from Garmin's .gdb files?

text: I have got the some Gps data in the form of Garmin's .gdb files. I have converted it to CSV files using GPSBabel. The CSV output contains Latitude, Longitude, and some kind of free form text remarks. I do not see any timestamps. Do Garmin's .gdb files contain timestamps? This is the first time, I'm working with this file format, so my understanding of the format/data is spotty at best. Where can I find more information about this format? How can I get the TimeStamps of the Gps points?


title: Using 3d Polyline in QGIS to calculate Volume?

text: This question may be completely off the mark. I have a survey drawing in DWG format, I use AutoCAD 2d but the survey was created in LSS. The Countour lines in the DXF have XYZ coordinates. Firstly how do i use this layer (i can save it out seperately) in QGIS and how do I display it? Secondly, is it possible to run a plugin or something (i dont know python) to calculate the volume of this area.... thirdly, i know that the maximum depath of the area is 7m below the highest point. How to create a new layer with the a z coordinate 7 m below maximum height...eg. max height is 110 therefor the z would be 103. Lastly is any of this possible and am i completely wrong in assuming this can be done? Any pointers or advice is appriciated.


title: Converting OSM_ID to coordinates

text: I have downloaded some open street map data from Webgeo, but instead of coordinates it came with an "OSM_ID". I cannot seem to find any data dictionary or information about this variable on the web. Does anyone know what it is and how do I convert it to coordinates? An example of the data is as follows:               osm_id,access,addr:housename,addr:housenumber,addr:interpolation,admin_level,aerialway,aeroway,amenity,area,barrier,bicycle,brand,bridge,boundary,building,capital,construction,covered,culvert,cutting,denomination,disused,ele,embankment,foot,generator:source,harbour,highway,historic,horse,intermittent,junction,landuse,layer,leisure,lock,man_made,military,motorcar,name,natural,office,oneway,operator,place,poi,population,power,power_source,public_transport,railway,ref,religion,route,service,shop,sport,surface,toll,tourism,tower:type,tunnel,water,waterway,wetland,width,wood,z_order     1432784186,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,St Aubin,,,,,village,,,,,,,,,,,,,,,,,,,,,,,     1787222125,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,London Berth,,,,,locality,,,,,,,,,,,,,,,,,,,,,,,     1787082727,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,,South Pier,,,,,locality,,,,,,,,,,,,,,,,,,,,,,,      Thanks for any help


title: how to tell google geocode() that you are looking for an administrative_area address and not a city?

text: If you were to do a reverse lookup for a region (States or Province) and that region happen to have a city with the same name, then Google geocoder will return the information on the city instead of the region. How would you pass parameters to Google telling it that you are indeed interested in the region and not the city? Calling Google WebAPI would be something like this: http://maps.google.com/maps/api/geocode/json?components=administrative_area:Durango|country:MX&sensor=false But instead, we want to call geocoder.               _geocoder = new google.maps.Geocoder();     _geocoder.geocode(request_dict, function(results, status) { ... }      How do you form the "request_dict"?


title: Trying to find a .NET GeoJSON library

text: I'm trying to find a .NET GeoJSON library which I can use to export my SqlServer2008 spatial boundaries (in binary format -> eg. `STAsBinary()`) from my .NET code to my HTML frontend to display on a google map. At first I felt that it might be a StackOverflow question, but this site I feel is more specialised to helping me out. Does anyone know of any? NOTE: This dll does NOT have to be imported into SqlServer2008. In fact, it will not. I'll use it in my C# MVC WebSite.


title: Hiding Rest Service Url in Javascript

text: While users using my application, they can view resource from browser and see my arcgis rest service URLs. So If I use token based security, security token can view. Can I hide this from users?


title: multipatch features are not visible in 3D ArcGlobe

text: I have created some 3D features in ROCKWORKS software, which represents lithology data. Then I converted it to 3D shapefile, and then I added it to ArcGlobe as well as ArcScene. It is clearly visible in 3D in ArcScene but in ArcGLobe it is not visible as 3D. Instead it is just visible as a 2D Feature class. Then later I have tried some extrusion options, but that graphic rendering is not good like ArcScene, and I have tried some other options like some sketch up files I have converted to multipatches but the same problem . Is there any way to add 3D symbology from ArcScene to ArcGlobe? Actually, I want to show lithology data like in the following link, and I want to create one globe service, how to do it..... http://blogs.esri.com/Info/blogs/arcgisexplorerblog/archive/2011/04/15/3d-subsurface- profiles.aspx Thanks and regards Srinivas


title: Help with joining a large dbf file to shape file

text: My overall objective is to create a timeline animation depicting when a person is working in a particular region. The .dbf table consists of over 6000 record entries of shift data, including date start and date stop times along with an assignment ID. I have approximatly 100 different shapes/regions also with assignment ID's. I would like for that region to turn on or off based of the table shift data and show when someone is working that region. When I join the table to the shape file, only the first shift entry gets attached and displayed. Do I have to have 6000 shapes for each of the 6000 record entries? There has to be a better way, please help. Thank you!


title: saving vectors in PostGis

text: I need to save my vectors layer back to postgis, which if i read correctly i can do using geoserver however im confused. All the examples ask you to import a shapefile. I dont have a shapefile. This will be a empty vector layer ontop of a couple of base layers which the user will add polys, points etc etc. How do i create the tables to hold the vector data. Also how do i know the structure of the tables. e.g columns, column names, primary / foreign key, relationships etc etc Thanks in advance


title: How to export multiple featureclasses to shapefiles?

text: I have plenty features in a geodatabase. I need to save/export them into shapefiles, because my mate works with Manifold, which is not able to read geodatabases. Any Ideas? Click on every single layer -> export is not a good option, due to the amount of features in that database.


title: Importing multiple shapefiles with tablespace and inheritance defined

text: In the UK we are very fortunate to have the Ordnance Survey who provide the whole of Britain in vector format. The data is provided as over 60 tiles, each of which contains something like 16 shapefiles for each feature type. Each shapefile is prefixed with the tile ID. There are very good threads here already that deal with the general case for importing all this lovely data into Postgis. For example:               find . -name *.shp -exec ogr2ogr -f PostgreSQL PG:'dbname=my-gis-database        active_schema=public host=localhost user=postgres password=my-password' '{}' -s_srs        EPSG:27700 -a_srs EPSG:27700 -overwrite -skipfailures -progress \;      My question is this. Is it possible to import all of the data in one simple script like this so that:   * Each tile inherits from a parent table e.g. HP_rivers INHERITS (rivers)   * The tablespace is provided to ensure it all loads onto a disk big enough to take it all in one go rather than my small but speedy SSD My goal is to have a single table structure which allows me to interrogate the whole of Britain or from a given tile, and to store the whole structure on a tablespace on slower but bigger disk. Any ideas very gratefully received!


title: Find arc endpoint with bearing and radius counting with curvature of the earth?

text: So my problem is that I have start point [lat,long] for arch, bearing, radius and I need compute end point [lat,long] ? End point should be in 180 degrees angle from start point and I should be aware of curvature of the earth ... So let say I have start point [54.01,104.05], bearing is 170 and radius is 2 km ... Any ideas? Thanks in advance.


title: How to select polygons associated with (not necessarily containing) corresponding points?

text: I have a list of polygons, each associated with the x/y coordinates of a point. These points may not intersect the polygons. How can I select the polygons that do contain the corresponding point? I only have a ArcView license.


What kind of unit of measurement for scale does SQL Server 2012 Spatial use?


title: arcpy.da.SearchCursor - unique value get second field value

text: I have a feature class and using arcpy.da.SearchCursor I'm getting unique values of a field (e.g. Name). How can I use this method to get the associated row value of a different field (diff field = Name2)? So far I use this to get unique values on only field Name:               values = [row[0] for row in arcpy.da.SearchCursor(fc, ("Name"))]         uniqueValues = set(values)         for value in uniqueValues:        When searching through the feature class how do I also get the value of another field, field Name2? Thanks


title: How to add several layers to the table of contents in ArcMap using Python in one shot?

text: I am already able to add individual layers to ArcMap table of contents using the following codes:               # Adding the layer to the table of contents     mxd = arcpy.mapping.MapDocument("current")     df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]     path = os.path.abspath(LayerName)     addLayer = arcpy.mapping.Layer(path)     arcpy.mapping.AddLayer(df, addLayer)     arcpy.RefreshActiveView()     arcpy.RefreshTOC()     del mxd, addLayer      However I can't figure out how I can load multiple layers at the same time. For example the following codes don't work and only only of the layers gets added:               # Adding mulitple Layers to the table of contents     mxd = arcpy.mapping.MapDocument("current")     df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]     path = os.path.abspath(FirstLayer)     path2 = os.path.abspath(SecondLayer)     path3 = os.path.abspath(ThirdLayer)     .     .     .     addLayer = arcpy.mapping.Layer(path)     addLayer = arcpy.mapping.Layer(path2)     addLayer = arcpy.mapping.Layer(path3)     .     .     .     arcpy.mapping.AddLayer(df, addLayer)     arcpy.RefreshActiveView()     arcpy.RefreshTOC()     del mxd, addLayer      Any thought/ idea is much appreciated.


title: Question related with NDVI value on water

text: There is a surreal thing on my NDVI images. NDVI values on lake is ranging from -0.04 to 0.7. (The lake surrounded by forest and has freshwater). Theoretically, negative value of NDVI is correspond to water. I want to mask NDVI image to remove pixels with negative values which include waterbody. But there are still positive values on lake. So I cann't set zero mask to remove whole waterbody. So what I'm gonna ask are: (i) why NDVI value on water has so high positive value? is there any explanation for that? (ii) how can I mask waterbody in order not to remove the other pixels? By the way, the NDVI images extracted from MOD13Q1 vegetation index product. Thank you so much I have already converted reflectance value into NDVI value as dividing by scale factor (10000). I uploaded 2 images. first image is NDVI image in april where value of lake is negative (correct). In April, lake is fully covered by snow. This image doesn't have any problem. ![enter image description here](http://i.stack.imgur.com/SUEWi.jpg) This one is NDVI image in June when no snow cover on lake. ![enter image description here](http://i.stack.imgur.com/XKXKs.jpg). when I check values of other lakes on image 2, they are in correct (negative) values. Only this lake has positive value. The name of lake is Khuvsgul with freshwater. So I'd like to ask whether green plant at the bottom of the lake can effect on the NDVI value during the period no snow cover on lake or not? I don't have any idea other than this guess.


title: PostGIS Geometry Query returns "Error: Operation on mixed SRID geometries" only for certain values

text: I have PostGIS table with two geometry columns, both defined with SRID 4326. I can insert into the table without problem, using the following `INSERT` statement (where `lng` and `lat` are values passed in programmatically):               INSERT INTO pad_meta (         uuid, created, updated, name, origin, radius, area, expiry, creator     ) VALUES (         $1, now(), now(), $2, ST_GeomFromText('POINT(lng, lat)', 4326), $3,          ST_Buffer(ST_GeomFromText('POINT(lng, lat)', 4326), $4), $5, $6     )      But when i query for an intersection using ST_Intersects, **depending on the value of the point** I get `ERROR: Operation on mixed SRID geometries`. For example, this query works:               SELECT * FROM pad_meta where ST_Intersects(         'POINT(-122.334172173172 46.602634395263560)'::geometry, area::geometry      ) ORDER BY created DESC;      And this errors out:               SELECT * FROM pad_meta where ST_Intersects(         'POINT(-122.334172173172 47.602634395263560)'::geometry, area::geometry      ) ORDER BY created DESC;      Note, **they are identical queries except the value of the longitude.** I have experimented with different values, but not identified a clear transition point between the queries that work and don't. I think I am fundamentally misunderstanding something. For the moment, I have solved/corrected/worked around the problem by re-formatting the query to use `ST_GeomFromText` and explicitly specifying the SRID:               SELECT * FROM pad_meta where ST_Intersects(         ST_GeomFromText('POINT(-122.334172173172 47.602634395263560)', 4326), area     ) ORDER BY created DESC;      But I honestly don't really understand what the difference is, or if this truly "the" solution. **My question is: Why am I getting an error only for specific values, and what is the proper way to format this query?** Here is my table definition for reference:               CREATE TABLE IF NOT EXISTS pad_meta (       uuid CHAR(32),       created TIMESTAMP,       updated TIMESTAMP,       name VARCHAR(128),       origin GEOMETRY(Point, 4326),       radius INTEGER,       area GEOMETRY(Polygon, 4326),       expiry TIMESTAMP,       creator CHAR(32),       PRIMARY KEY (uuid)     );      I have also verified that there is only one type of SRID in the geometry_columns:               SELECT f_table_name, f_geometry_column, srid FROM geometry_columns;     f_table_name | f_geometry_column | srid     --------------+-------------------+------      pad_meta     | origin            | 4326      pad_meta     | area              | 4326      Help/advice appreciated. Thank you! (Note: I have also seen this question, but since I am already explicitly defining my geometry SRIDs when inserting into the table, it seems like that is not what is happening.)


title: same or similar search word with field calculator

text: I am using arcgis 10.1 I have two tables A and B need to find attributes of table B data are on the table A I'm doing a join between table A and B now I want to use the field calculator via a script to compare field "name" field of table A and table name B I need help on how to make this script


title: Find all census tracts N miles from a lat, lon point

text: Using R, I'd like to find all census tracts within a given distance of a point (lon,lat). i.e. All census tracts within 20 miles of 34.0522, 118.2428 I have used the sp library and over command to locate the tract of a given point, but would like to use it to find all the tracts covered by a "disc". Ideas?


title: Create a point x meters from origin using st_translate in postgis

text: How do I create a point 1000 meters away from a reference point? The ST_Translate only supports geometry type. I also tried using ST_GeomFromText but that returns a unit and not in meters. Here is my sample query:               SELECT ST_AsText ( ST_Translate ( ST_GeogFromText('SRID=4326;POINT(-70.01 43.27)'), 1000, 0 ) );      Is there a way that I can convert the output of ST_Translate into meters?


title: How can I convert Polygon to LinearRing Geometry?

text: How can one convert a Polygon to LinearRing geometry using OpenLayers? I'm thinking that I can do this by getting bounds of the polygon, then use this to construct a Points array. Finally, use OpenLayers.Geometry.LinearRing(Points_Array) constructor to create the LinearRing. My attempts so far have flopped. Any code samples on how to achieve this will be highly appreciated. UPDATE: I've accepted @drnexgis' answer for what it taught me on converting geometry to GML format. It only partially answered the question, but it was helpful.


title: Why are there no more error messages after importing QGIS modules into python?

text: I have a problem when writing python code whilst including the qgis modules. As can be seen below, when I usually call a variable that does not exist (for example), I receive the usual error message. Just what you'd expect to happen. But after having imported qgis.core and qgis.utils into python, all error messages disappear (see below). This obviously makes code debugging extremely tedious, as I never get any indication, why and where a code failed.               >>> print this_variable_does_not_exist     Traceback (most recent call last):       File "<stdin>", line 1, in <module>     NameError: name 'this_variable_does_not_exist' is not defined     >>> from qgis.core import *     >>> from qgis.utils import *     >>> print this_variable_does_not_exist     >>>                          (<< why is there no error here?)      Does anyone have an idea, how I could stop this "error-suppressing" behavior? Thanks in advance, I appreciate any help!


title: Can I reproduce this classic cartographic hillshade relief effect (hachures) electronically using modern DEM data?

text: I own a number of old cetral-European maps and I always liked the classic hand-drawn cartographic hillshade relief presentation that these maps feature. I am not wanting to reproduce the shading (which is similar to swiss shading) but the accompanying lines. These are not contours but rather **down-slope "gradient" lines** that (often with a degree of imagination) create a hillshade effect. I don't have my maps on me right now so I can't scan a good detailed example of this but I found some examples online. As far as I know, and at least in Czechoslovakia, these types of maps were being produced up until the second world war. I would like to produce some modern maps with an old feel and featuring this type of effect. **Is there any way that this could be done in a modern GIS (Preferably ArcGIS) using an existing DEM?** Illustrator or Photoshop solution would also be acceptable. Not a part of the question but I would also be interested to find out: What is this particular hillshading called? Where could I get some more information about how it was done and what rules it adheres to? ![enter image description here](http://i.stack.imgur.com/jHCnc.jpg) ![enter image description here](http://i.stack.imgur.com/FwZ4o.jpg)


title: Generating vector GeoPDF from ArcGIS Server services

text: I am exploring possibilities to generate vector GeoPDF from data which is accessible from ArcGIS Server services. I have ArcGIS Server MapService with complex styles. They are very simplified in FeatureService - so it is not an option. I have already tried to use **Export Map** operation of MapService and was lucky to get nicely styled vector PDF from it. Of course, I have a coordinates of BoundBox (this info actually is used for querying). So, now I need to include also additional data (some labels, borders and other pre-print graphic) and generate result PDF. It seems like a task for FME Desktop. But! I have not found how I can import input PDF from http query. In fact it seems that I can't even open PDF from disk in FME, it is capable only to write it. The question is: what possible solutions do you see for this task? The key thing is to keep complex styling from ArcGIS Server. I can also access to mxd which was used for creation of service. PS I need solution which is possible to run from scripts, as it will be scheduled task. **EDIT1** May be there is any solution to combine several PDF files into one in Python? Some advanced Python lib for editing PDFs? Thanks!


title: How to close multipart polyline gaps in ArcMap

text: I have a multiple multipart contour lines that originally came from CAD drawing. I converted these into a polyline shapefile in ArcGIS. The issue I have is with gap in each contour lines that are direct result of CAD text annotations for contour value. These gaps differ in size. Is there a way to snap multiple lines to their closest nodes at once in ArcMap (or using any other tools)? I attach a screenshot of the contours from ArcMap: ![Screenshot](http://i.stack.imgur.com/rWqJm.png) Many thanks in advance, Magda


title: Is there an OpenLayers event when reaching bounds?

text: I have an app using OpenLayers and OpenStreetMap tiles. There are more zoom levels for the capital and less zoom levels for the entire country. Let's say the user is in the capital, in his car, in traffic, with localisation active. The app is showing the map at the maximum zoom level. Them map is following him like a GPS, keeping current position in the center and marking it with a graphic symbol. When the car reaches the bounds of the city the map have no tiles and show pink color. I need to change the zoom level before the pink missing tiles appears. I can calculate an area inside the bounds of the capital and when the current position is inside this area the app will zoom out to the maximum zoom level available for the country. Is this the best solution? The same problem exists when panning map with the finger. I thought to restrictExtend to the capital but then the user will be captive inside and the single chance to see the country again is to zoom out. Do you know an algorithm or an OL function/event to help me find the best solution? Thank you, Bogdan


title: Looking for a FREE Cloud based 'Web Mapping Stack' Host

text: Any pointers to free cloud-based geoweb hosts? Use case: To display a basic interactive web map in Leaflet/ OpenLayers with data hosted in PostGIS served by GeoServer. I am looking at something very close to the OpenGeoSuite since I am playing with that at the moment already. [This is basically a teach- yourself webmapping exercise and there are zero personally hosted server resources on my sleeves]. I've found capnix to host the website but now need the 'Geo'-portion. CartoDB's 5MB data limit is too small for my project needs and OpenGeo's Cloud Edition is too short (3 days!). Any ideas?


title: How to create floating-only window in ArcMap?

text: I'm looking for a way to create a floating window in ArcMap. To give an example just look at the window of the Identify tool. Floating means that it stays in front of the map document at all times and that the user can continue to work with ArcMap. I know that the interface IDockableWindowDef can be used to create dockable windows which can also float, **but I don't want them to dock**. To my knowledge it's not possible to prevent forms created by IDockableWindowManager from docking if e.g. the user pushes them to the border of the ArcMap window. Any ideas? [edit] The solution was to search for keywords like child window and MDI. HTH [edit2] The problem's solution seems to be as easy as this (see llcf's answer):               MyForm form = new MyForm();     form.Show(NativeWindow.FromHandle(new IntPtr(m_application.hWnd)));


title: What raster-based GIS systems really work?

text: Raster GISes like GRASS, ArcGIS/Spatial Analyst, and Idrisi can perform a rich set of data processing and analytical procedures loosely known as "map algebra." In today's computing environment it is becoming commonplace to maintain rasters of 100,000,000 cells or more in many different formats and to demand relatively complex calculations like viewsheds, watersheds, and terrain identification, as well as image processing capabilities. It seems that many open source, free, and inexpensive solutions are out there. But which ones really hold up in practice? That is, which ones can handle large grids efficiently, can easily get data in and out, are reasonably bug free, and offer a full complement of analytical procedures? What are the pitfalls or hidden limitations that you don't find out until you have invested lots of time in learning these systems? (This last question is one that is not easily answered with Web searching and where I am hopeful respondents can offer valuable advice.) I am especially interested in solutions that can both integrate well with and compete with the popular (but expensive) commercial systems (which means Windows compatibility is important).


title: How do I create a negative to positive value raster with neutral 0?

text: I have a raster layer with a band with values from -137 to 9229. I want to create a color scheme showing how far from 0 the positive numbers are in one color, and how far from zero the negative numbers are in another color, with a neutral 0 (that, unmistakenly, doesn't lean towards the positive or negative colors, nor diverges from the color scale). Moreover I'd like a "-55"-area to be just as far from the neutral color (0) as a "+55"-area, but in different colors. So if neutral would be #333, then -5 might be #338, and +5 would be #833, even though the array of values, as in this case, is heavily biased towards positive numbers. Hope that explains my problem. I'd greatly appreciate any help I get! This is for my Bachelor thesis project on the creation of the metro and highway system in DC during the 50s to 80s. EDIT: Oh, gosh. I forgot to tell you I'm using QGIS 1.7.4 on an OSX Lion. And have very limited access to Cartographica and ArcGIS full suite. Best regards, Henrik


title: OpenGL samples for use with ArcEngine? Inspired by video I linked in post

text: Inspired by this video: http://www.youtube.com/user/mresri#p/a/u/0/KjcDWG8nV1g What are some good OpenGL resources for use with ArcEngine? C# ArcEngine 10


title: Suggestions for simple mapping platform for alumni web site

text: I'm working on a proof of concept alumni site with a WordPress backend. I want to have a very simple map that shows where alumni are living. I'm thinking something like the Bing API or a simple open source solution - points on a map right now. Not sure yet how I will be getting locations, perhaps a user submission system somehow. Any ideas on good, solid, simple ways to accomplish something like this? I'm on a LAMP stack hosted with WebFaction, and I know just enough PHP to be dangerous. EDIT: Here's what I ended up with. Jlivni got me looking the Google route and some hacking on this method gave me this: ![alt text](http://i.stack.imgur.com/nZM6B.png) It's powered on the back end by a publicly viewable Google Spreadsheet. Pretty easy to implement and hack up. And I can pass of the editing of the spreadsheet to just about anyone and they can carry out that task.


title: Is there a plugin that collects and allows you to just "turn on" open layers?

text: I'm thinking on a similar approach as gapminder.org, only that this plugin would contain GIS-data, not only statistics. I'm specifically interested in if there is a plugin like this for QGIS. If there isn't one, why is that? Isn't it genius? =} Thanks


title: What should I pass to ILayerFactory.Create for a WMSLayerFactory?

text: I was expecting ILayerFactory.CanCreate to return true for an IWMSLayerDescription. I recursively traverse all the descriptions in a WMS connection but CanCreate always returns false.


title: Problem viewing SQL 2008 "Date" field

text: I am using ArcMap Basic 10.1 SP1 and I'm connected to a permitting database stored in a SQL Server 2008 R2 database. I am trying to pull the tabular data over into my ArcGIS Desktop Project so that I can do a join with the permit boundary information based on permit numbers.   One of my fields in the SQL Server database is called "DateIssued" and it is of SQL datatype "date" rather than "datetime". When I drag over my table from my "Database Connections" area of Catalog, all the fields come over except those marked as "Date" and work fine. The "Date" items are simply not listed as a field. I tried to get around this by making some SQL views to try to CAST or CONVERT the SQL data to be a datetime or even a varchar, but then when I drag the view over, it gives me the error message that "Warning - Could not add the specified data object to the map. Invalid SQL syntax" To be honest, I am not sure that ArcGIS 10.1 even supports using the SQL views as I seem to get that message on any views. Has anyone run into this problem with the SQL "date" data type? I thought I might generate some temporary tables to do the date to datetime conversion to see if that might work. I'll post my findings here as well.


title: GeoDjango WebMapping Application Update

text: According to the below GIS.SE questions, GeoDjango was the preferred framework for creating simple and quick webmap applications. But these questions are now all 2-3 years old. Can we verify whether or not consensus says that GeoDjango still remains top among webmapping application frameworks? What are the advantages of using Rails for a GIS/mapping application? What are the Python tools/modules/add-ins crucial in GIS? Looking for a developer-friendly web gis


title: WMS Pop Up not working....?

text: I have created a web based GIS using Geoserver-Postgis-Openlayers... Every thing was working fine till I changed the server path from 'Testing to LIVE.' WMS POP up Function has stopped working in LIVE Server Environment & showing the following error message.. > "Message: Access is denied. > > Line: 177 ; Char: 9 ; Code: 0 ; URI: > ....../geoserver/OpenLayers-2.8/lib/OpenLayers/Request/XMLHttpRequest.js " WMS Pop Up code is as follows (Working perfectly at my Test Server)               info = new OpenLayers.Control.WMSGetFeatureInfo({                 url: '..../geoserver/LT_WEB/wms',                  title: 'Identify features by clicking',                 queryVisible: true,                 eventListeners: {                     getfeatureinfo: function(event) {                         map.addPopup(new OpenLayers.Popup.FramedCloud(                             "",                              map.getLonLatFromPixel(event.xy),                             null,                             event.text,                             null,                             true                         ));                     }                 }             });             map.addControl(info);             info.activate();                  map.addControl(new OpenLayers.Control.LayerSwitcher());             map.zoomToMaxExtent();      SERVER - Microsoft Windows 2003 R2. Web Browser - Internet Explorer 8.0. In both Servers (Testing & LIVE) Geoserver,Postgis & Openlayers 2.8 installed & all paths Verified. Awaiting for your valued response. Thanks in advance. Sandipan


title: Split polygons programatically

text: I have looked through several other questions and am still having trouble. I have a set of 103 polygons that span my entire map. Each of those polygons is divided into 36 equal sized polygons. Now I need to iterate through each of those 36 squares, and split them into 100 equal squares. I've tried getting in contact with the person that split the first 103, but all efforts have failed. Can anyone point me in the right direction? Fishnet seems to be out of the question as it spans the entire feature class, not the polygons. I'm guessing I will be writing a python script. Any pointers would be appreciated! ![What I have to work with.  103 Townships each split into 36 sections.  I need to split the 36 sections each into 100 equal squares](http://i.stack.imgur.com/S3Haa.png)


title: Is there a way to "disperse markers" across multiple point layers?

text: I have a bunch of point layers which often have overlapping features. The disperse markers geoprocessing tool works within a layer but is there a way to do it across more than one layer, so none of the points are overlapping? Using ArcMap 9.3.


title: How to geocode addresses that are outside Address Ranges of source layer?

text: Im using TIGER lines to geocode a table of addresses. Say for instance the address range for a street is 1100-1200. What would be the best way to match an address such as 1094?