Problem in loading the legend of MapServer in ArcGIS

I have published a mapserver on my ArcGIS Server 9.3. I am developing a silverlight application using ArcGIS Silverlight API. Now I can load my map into my application and also the list of my layers that I can toggle their visibility. I just want to load my map legend in order to know the legend of each layer. The solution which I find in arcgis silverlight example is just as follow:               <esri:Legend Map="{Binding ElementName=Map}"          LayerIDs="BaseLayer"         LayerItemsMode="Tree"         />      In visual studio 2010, I have an error for LayerIDs part which says that "Unable to cast from....". The definition of my map is also:               <esri:Map x:Name="Map"  Background="White" MouseClick="QueryPoint_MouseClick" IsLogoVisible="False">         <esri:OpenStreetMapLayer ID="OpenStreetMap" Style="Mapnik">         </esri:OpenStreetMapLayer>              <esri:ArcGISDynamicMapServiceLayer ID="BaseLayer"                    Url="http://sina-pc/ArcGIS/rest/services/Test/MapServer"         Initialized="ArcGISDynamicMapServiceLayer_Initialized" />                 </esri:Map>      It also doesn't show any of my legends for BaseLayer map. Am I missing something or is there any problem with my map? Thanks for your helps.


Is it possible to export a layer with a related table to kml/kmz and view related data in Google Maps?

I have a polygon shapefile with a related table in ArcGIS 10. When I use the "Layer to KML" tool to export the shapefile and then load it into Google Maps, it doesn't see the related data. I would like the related data to appear in the balloon pop-ups in Google Maps. Is this possible? Thanks, Vanessa


Run FME Workbench with Python

I would like start/run FME workbench from one external python script. Is there some simple syntax/command to do this?


Why won't lines from GPX simplify/generalize/smooth?

I am working with GPX data. I imported them to ArcGIS by 'GPX to feature' and then I created lines from these points. I am trying to smooth (or simplify) the lines now, but the tools 'smooth line' and 'simplify line' are not working. (Maybe it is because of the big amount of points..) Do you know any tools or hints how to smooth these lines created from GPX? Thank you for your answers


HP Plotter 1050c plus

I have an (1998 manual) old plotter HP PLotter 1050c Plus (only works on Windows XP)... It prints fine, but the Source/Receivers are not showing up on the map nor the legend. Everything else is displayed. It is printing off the Adobe PDF created by ArcGIS 10.1 Not sure what the problem is, I have increased the size of the Src/Rec, but not luck. Troubleshooting manual is no help either. Has anyone experienced this problem and fixed it?


Base station coverage analysis in QGIS

I want to find base station coverage for a instant location. I have building data. And i determine optimum coverage as 500m far away from a base station (buffer). So how can i find optimum base station location for this example. Looks like that ![enter image description here](http://i.stack.imgur.com/i1PHi.png)


Setting More Than One Workspace

I have a number of tables that contain known values in a field. For example, I have a soil type table with one field in it. That field contains a list of known soil types. I have another table that has a list of known vegetation types. I have various other tables that contain lists of other known types of other data. What I want to do is:   * Read in each table, do a cursor search and add the values to a list via `tableList.append`   * Now that I have built each list, I want to compare the values from a separate feature class table to the known values in the list I have just created using `if row.SOIL_TYPE in SoiltableList:`   * In order to do this, I need to reset the workspace I guess, because you cannot set a path parameter for listing data in any of the arcpy Listing methods.              env.workspace = r"C:\BASE_TABLES.gdb"          SoiltableList = []     VegtableList = []          for table in arcpy.ListTables():           if table == "SOIL_TYPE":           for row in arcpy.SearchCursor(table)               SoiltableList.append(row.SOIL_TYPE)        if table == "VEG_TYPE":           for row in arcpy.SearchCursor(table)               VegtableList.append(row.VEG_TYPE)          # more iterations of table          # Reset workspace     env.workspace = r"C:\TEST_GDB.gdb"          for fc in arcpy.ListFeatureClasses():        if fc == "SOIL_FEATURECLASS":           for row in arcpy.SearchCursor(fc):              if row.SOIL_TYPE in SoiltableList:                 print "YES"              else:                 print "No"       # More iterations of if/else logic to compare values from fields from various feature classes to the known values in the list I built.      The 2 questions I have are:   1. Do I need to reset the workspace the way I did above? is there a better way of doing this or is what I am doing completely acceptable?   2. Does my code make sense? Am I over complicating what I am trying to do? I though that maybe I would define a function or 2 to simplify the code, but I'm still in the process of figuring out what else I want the code to do. Thanks for any suggestions. Mike


"The tool not licensed" ERROR on a geoprocessing service on ArcGIS Server 10.1

I have a python script which runs well as geoprocesing tool. When published on ArcGIS Server 10.1 as geoprocessing service, however, it throws the following error:               Unable to complete operation.     Error executing tool.: Failed to execute. Parameters are not valid. ERROR 000824: The tool is not licensed. Failed to execute (Kriging). Failed to execute (NewCuringModelExpPy). Failed to execute (New_Curing_Model_ExpPy).      The script runs without any error on ArcCatalog with Spatial Analyst extension authorized. The code also has got proper license checked out as below:               arcpy.CheckOutExtension("spatial")      I don't think ArcGIS Server requires Spatial Analyst extension authorized, does it? Otherwise, is it to do with any ArcGIS Server configuration? Any idea is appreciated. By the way, I am running on ArcGIS Server 10.1 Advanced so Spatial Analyst extension should have been part of it. Cheers,


Openlayers – Empty spaces when resizing map in height (opera and safari)

The following test site is supposed to let the map extend to the maximum available space and to ensure a minimal height of 300px. This works perfectly in firefox. But in opera and safari on all apple devices (Mac, IPhone, IPad), after resizing the window a few times will create empty areas (screenshot).   Js fiddle
Jemand musste Josef K. verleumdet haben, denn ohne dass er etwas Böses getan hätte, wurde er eines Morgens verhaftet. »Wie ein Hund!« sagte er, es war, als sollte die Scham ihn überleben. Als Gregor Samsa eines Morgens aus unruhigen Träumen erwachte, fand er sich in seinem Bett zu einem ungeheueren Ungeziefer verwandelt. Und es war ihnen wie eine Bestätigung ihrer neuen Träume und guten Absichten, als am Ziele ihrer Fahrt die Tochter als erste sich erhob und ihren jungen Körper dehnte. »Es ist ein eigentümlicher Apparat«, sagte der Offizier zu dem Forschungsreisenden und überblickte mit einem gewissermaßen bewundernden Blick den

                  Ps: At the time of writing the OSM map server seems to be down. (Working again)


How to make Story Map with two photos that can swipe?

I am trying to make a story map with some heritage data. I want to have a picture of each heritage site two times(Swipable: that i can swipe between the two photos), one photo for the current and one for earlier years. It can be a swipe or any other tool that is available as long as it is not the same heritage sites for two points. I am open to suggestions on the platform but i am familiar with ArcGIS Online.


arcpy: how to access element in list?

I would like to know how to get the elements from the list, maybe using for loop?... here are the code that can save all the angles into myList[]               myList = []          for num in range(1,totalStop+1):         in_features = "C:/Users/Desktop/fyp/fyp.gdb/sampleline"         out_feature_class = "C:/Users/Desktop/fyp/output.gdb/samplelineTemp"             where = '"OBJECTID" = ' + "%d" %num              # Execute Select         arcpy.Select_analysis(in_features, out_feature_class, where)                   for row in arcpy.da.SearchCursor(out_feature_class, ["OID@", "SHAPE@","streetnam", "angle"]):                 # Print the current line ID                 print "Walk along",row[2]                      #Set start point                 startpt = row[1].firstPoint                      #Set Start coordinates                 startx = startpt.X                 starty = startpt.Y                      #Set end point                 endpt = row[1].lastPoint                      #Set End coordinates                 endx = endpt.X                 endy = endpt.Y                      #print startx, starty, endx, endy                              value = row[3]                                      #append all the angles of line into myList                 if value not in myList:                     myList.append(value)      after append(value), now myList contain all angles(azimuth) of the lines. mylist[]= -91.3478982503 -3.51652064715 -124.295923979 31.1849858238 43.8914254459 -59.4335738611 111.914929118 69.5778228682 The thing I would like to calculate is the first angle(myList[0]) minus the second angle (mylist[1]), return the value. Then the second angle(myList[1]) minus the third angle (mylist[2]), return the value. May anyone kindly tell me how to perform the above steps? Loads of thanks!


Google Map With Search Box and Dragable Pin

I'm trying to build a map with a search and with drag-able pin and the results will give an address and Latitude Longitude. I'm trying this code but something doesn't work and I cant find it... Please help                          <meta name="viewport" content="initial-scale=1.0, user-scalable=no" />     <link href="http://code.google.com/apis/maps/documentation/javascript/examples/default.css" rel="stylesheet" type="text/css" />
     <input id="address123" type="textbox" value="Sydney, NSW">          <input type="button" value="Geocode" onclick="codeAddress()">

         Marker status:
Click and drag the marker.
         Current position:

         Closest matching address:


Obtain normalized raster average per polygon?

I have a raster that contains percentages of total forest killed plus a polygon layer that defines different zones. My goal is to find a good statistical representative that describes the forest killed situation per zone. So what I did first is I used zonal statistics as table with the raster as the input raster and the zones as polygons as the feature zone. This calculated me an average per zone based on the raster cells in the zone. However, this destorts the result. One zone might only have 5 raster cells in there but all of these cells have a value of 50%, which gives me an overall average of 50% for this zone. Another zone might have 1000 raster cells and obtains also an overall average of 50%. So that doesn't show that the latter zone has more raster cells with fores killed. So I need to normalize it somehow. What I thought I could do is try and calculate a weighted average. So I would need to know how many raster cells there are in total and how many cells there are within each zone. With that I could divide the amount of raster cells per zone by the total amount of raster cells in the whole region and then times that by the average of the raster cells per zone? Is that the right approach to not misrepresent the data? Also...how would I do that in ArcMap10??


How to add rasters of different sizes with Arcmap or Qgis raster calculators?

I've got 8 rasters all of different sizes. What is the proper syntax, for the raster calculators of either Arcmap or Qgis, that will output a raster that covers the extent of all my data not just the area of overlap? Thanks.


How to get $length in meters in QGIS?

I read a lot of posts concerning that subject, nevertheless I cannot solve to problem of getting $length in the field calculator in meters rather than in decimal degress in QGIS. The CRS of original shapefile is WGS84 (EPSG:4326). I reprojected the shape to World Euidistance Cylindrical (Sphere) (EPSG:3786). I did the same in project settings (with on the fly transformation enabled). However I am still getting the length of the polylines in decimal degrees. Any ideas?


How to render only those Leaflet layers which are in current map view?

Is it possible to render Leaflet layers which are in current map view so that other layers which are not in that map view doesn't render? I have many polgons (buildings) and it is really slow rendering if they all renders so I thought it would be ok to divide buildings into many shapefiles which will render only in that map view. Thanks.


Looking for ArcPy tool for calculating statistics on field with summary field AND boolean fields

What I'm searching for is a tool that will sum a field by a secondary unique field AND booleans... Let's say I have 12 table each with 4mil records and the following fields:   * Donation_Amount ($ amount of a donation)   * County_FIPS (county from which the donation was made)   * Republican (y/n)   * Democrat (y/n)   * Independent (y/n)   * In_district (y/n)   * Individual (y/n)   * PAC (y/n) I want to summarize donations per boolean for each unqiue county fips. Ideally the result I'd like is a table of all unique FIPS and an amount sum and count for Repub, Dem, Ind, In_District, Indiv, PAC. Tried to get it to work with iterate in model builder and couldnt... Also read about summarize and frequency tools but Summarize seems to be limited with how many secondary fields you can use, and Frequency tool is only for ArcInfo (I need for ArcView 10.1) Essentially I know this is simply a database process to summarize data based on secondary fields... I have heard of tools in R that do this (in Plyr and Reshape2: cast/melt) but would ideally like to do it in a program with GUI and even more ideally, in ArcMap so it wouldnt require too much shuffling of data... Thanks all! Much appreciated.


Method of extracting the intersecting areas of multiple overlapping polygons (5) and classifying intersects based on spatially joined points

I have concave hulls surrounding clusters of data points which has led me to a lot of overlap and fuzzy boundaries. I am attempting to intersect the overlapping areas but have a minimum of 5 different layers overlapping. Finally I would like the count to determine the classification of the intersecting polygons based on the dominating count. Any help?


How to find X and Y of clustered features?

How can I find longitude latitude values of clustered features? Please look the image. I hope you understand my question. I need to find out the X and Y of clustered features on every zoom levels. ![enter image description here](http://i.stack.imgur.com/hdrgn.jpg) Sample Code:               var defaultStyle = new OpenLayers.Style({             strokeColor     : "#FFFFFF",             strokeOpacity   : 1,             strokeWidth     : 2,             fillColor       : settings.color,             fillOpacity     : 0.9,             cursor          : "pointer",             pointRadius     : "${point_radius}",             pointerEvents   : "visiblePainted",             label           : "${label}",             fontColor       : "#000000",             fontSize        : "12px",             fontFamily      : "verdana",             fontWeight      : "bold",             labelAlign      : "center"          }, {             context: {                 label: function(feature) {                     var sumValue = 0;                     if(feature.attributes.count>1) {                         for (var i=0; i<feature.cluster.length; i++) {                             sumValue++;                         }                         return sumValue;                     }else{                         return 1;                     }                 },                 point_radius: function(feature) {                     var pix = 10;                     if(feature.cluster) {                         pix =Math.min(feature.attributes.count, 30) * 2.5 + 8;                     }                     return pix;                 }             }         });              var onSelectStyle = new OpenLayers.Style({             strokeColor     : "#000000",             fontColor       : "#FFFFFF",             pointerEvents   : "visiblePainted",             fontSize        : "15px"         });              var strategy = new OpenLayers.Strategy.Cluster({distance: 30, threshold: 2});              var ly=new OpenLayers.Layer.Vector(GSMapLayerName,{                 strategies: [strategy],                 styleMap:new OpenLayers.StyleMap({                             "default":defaultStyle,                             "select": onSelectStyle                         }),                 displayInLayerSwitcher:false             });         map.addLayer(ly);      Note: I am using openlayers 2.11 Thanks for any reply


How to fetch openlayers layers from pyqgis?

My objective is to fetch maps provided by openLayers, as in openLayers_plugin, from a stand alone python script. I would like to get it superposed to my shp file. As can be read on http://www.qgis.org/fr/docs/pyqgis_developer_cookbook/loadlayer.html               layer = QgsVectorLayer(data_source, layer_name, provider_name)     if not layer.isValid():print "Layer failed to load!"      provide a good way for loading layer in a stand alone python script. And it's possible to do it as well for wfs data source, or other data source types.   I do it with shp files, and it works well, but what is the uri related to openLayers datas and maps ? I also tried to get openLayers_plugin ...               sys.path.append(r"C:\usr\python\plugins\openlayers_plugin")     qgis.utils.loadPlugin('openlayers')          qgis.utils.startPlugin('openlayers')     print qgis.utils.isPluginLoaded('openlayers') # return False      Is it possible to use openlayers plugin in stand alone pyqgis script ? cf [ Can I use QGIS plugins in PyQGIS as a stand alone script without iface? ] which is deprecated now. No it's not ? thank you for any help.


Cannot Display all KML attribute data in Google Maps

I've linked a hosted KML file in to my custom Google Maps test website, the KML file is something I created from a clipped shapefile. I am struggling to view the attribute info once I have clicked on a KML feature, in my case a polygon. I opened the KML file in word pad and all of the attribute info is present. In Google Earth I can see the attribute data, in Google Maps I have to click the actual polyline of the KML data, only then does a small popup appear but there is no attribute data. I have copied the code from my JavaScript file, if you need my html and css I can copy that in too. Please forgive my ignorance, I am a GIS tech and I have limited knowledge when it comes to scripting languages. The KML file is hosted here: http://dodgerza.net23.net/GM_WARDS_Clipped_Final.kml Thank you.               // JavaScript Document              // initialise variables         var startLat = 53.479243;         // initial centre of the map         var startLong = -2.244239;         var startZoom = 12;                   // inital scale of the map         var my_map;         // stores the text details for each point of interest         var POI = new Array(15);         // establish the main icon         image = "POI_icon.png";         // stores Place 'objects' in an array         var marker_info = [];         // marker for Town Hall         var Town_Hall;         // infowindow varaible         var infowindow;              // directions service variables         var directionsDisplay;         var directionsService;                             // add weather service, toggle on and off              function initialize() {             var latlng = new google.maps.LatLng(startLat, startLong);             var weatherOn = true; //starts off false, weather layer is not on by default              // setting up map options         var mapOptions = {             center: latlng,             zoom: startZoom,              };              // assigning map to its variable         my_map = new google.maps.Map(document.getElementById("map_box"),             mapOptions);         // add weather layer in as service          var weatherLayer = new google.maps.weather.WeatherLayer({             temperatureUnits: google.maps.weather.TemperatureUnit.CELSIUS,         });         // setting a listener that will toggle the weather layer         google.maps.event.addDomListener(document.getElementById("weatherToggle"), 'click', function () {             if (weatherLayer.getMap() != null) {                 weatherLayer.setMap(null);             } else {                 weatherLayer.setMap(my_map);             }         });                   // adding KML layer         // census Ward data         var wardLayer = new google.maps.KmlLayer({             preserveViewport: true, // preserve location of viewport               suppressInfoWindows: false,             url: 'http://dodgerza.net23.net/GM_WARDS_Clipped_Final.kml'             //             });         wardLayer.setMap(my_map);         // setting a listener that will toggle the ward layer         google.maps.event.addDomListener(document.getElementById("wardToggle"), 'click', function () {             if (wardLayer.getMap() != null) {                 wardLayer.setMap(null);             } else {                 wardLayer.setMap(my_map);             }         });                             // add traffic         var trafficLayer = new google.maps.TrafficLayer();         trafficLayer.setMap(my_map);         // setting a listener that will toggle the traffic layer         google.maps.event.addDomListener(document.getElementById("TrafficToggle"), 'click', function () {             if (trafficLayer.getMap() != null) {                 trafficLayer.setMap(null);             } else {                 trafficLayer.setMap(my_map);             }         });         //                }     google.maps.event.addDomListener(window, 'load', initialize);                    // function that creates a user-defined 'Place' object     function Place(name, url, latitude, longitude) {         this.name = name;         // 'this.property' is the same as 'Place.property         this.url = url;         this.latitude = latitude;         this.longitude = longitude;     }          // Set up the data for the Places          POI[1] = new Place('The Royal Exchange Theatre', 'www.royalexchange.co.uk', 53.4829138, -2.244502);     POI[2] = new Place('Cathedral Gardens', 'www.manchester.gov.uk', 53.485237, -2.242042);     POI[3] = new Place('Piccadily Gardens', 'www.manchester.gov.uk', 53.481049, -2.236826);     POI[4] = new Place('St. Johns Gardens', 'www.manchester.gov.uk', 53.477666, -2.252317);     POI[5] = new Place('YMCA Training', 'www.ymcatraining.org.uk', 53.48039, -2.23418);     POI[6] = new Place('Convention Complex', 'www.manchestercentral.co.uk', 53.477194, -2.24641);     POI[7] = new Place('John Rylands Library', 'www.library.manchester.ac.uk', 53.480274, -2.24876);     POI[8] = new Place('The Hidden Gem', 'www.hiddengem.catholicfaith.co.uk', 53.480185, -2.246432);     POI[9] = new Place('Manchester Art Gallery', 'www.manchestergalleries.org/', 53.47882, -2.241822);     POI[10] = new Place('Deansgate Station', 'www.nationalrail.co.uk/stations/dgt/details.html', 53.474402, -2.251033);     POI[11] = new Place('Manchester Oxford Rd. Station', 'www.nationalrail.co.uk/stations/MCO/details.aspx', 53.473961, -2.241806);     POI[12] = new Place('Manchester Piccadilly Station', 'www.nationalrail.co.uk/stations_destinations/MAN.aspx', 53.47732, -2.230718);     POI[13] = new Place('Salford Central Station', 'www.nationalrail.co.uk/stations/SFD/details.aspx', 53.48336, -2.255255);     POI[14] = new Place('Manchester Victoria Station', 'www.nationalrail.co.uk/stations_destinations/MCV.aspx', 53.487382, -2.242649);     POI[15] = new Place('NCP Parking lot 1', 'www.ncp.co.uk/', 53.479443, -2.246413);     POI[16] = new Place('NCP Parking Lot 2', 'www.ncp.co.uk', 53.479544, -2.248215);                    //  ---------------------------------------          function draw_the_map() {              var latlng = new google.maps.LatLng(startLat, startLong);         var myOptions = {             zoom: startZoom,             center: latlng,             mapTypeId: google.maps.MapTypeId.ROAD,             panControl: true, //make pan control visible             panControlOptions: {                 position: google.maps.ControlPosition.LEFT             }, //place the pan control on the left             zoomControl: true, //make zoom control bar visible             zoomControlOptions: {                 position: google.maps.ControlPosition.LEFT             }, //place the zoom control on the left             scaleControl: true, //make scale bar visible... position fixed at bottom right             streetViewControl: true //enable street view.. position automatically appearing above zoom controls         };              // my_map = new google.maps.Map(document.getElementById("map_box"), myOptions);                   // the following code adds marker for Manchester Town Hall         latlng = new google.maps.LatLng(53.479243, -2.244239);         Town_Hall = new google.maps.Marker({             position: latlng,             map: my_map, //map object             title: "Manchester Town Hall"         });              directionsDisplay = new google.maps.DirectionsRenderer();         directionsService = new google.maps.DirectionsService();                   for (var i in POI) {             marker_info[i] = new Object;             desc = '
 ' + POI[i].name + '
 Website:' + '' + POI[i].url + '
';             addMarker(i, POI[i].latitude, POI[i].longitude, desc);         }     }          // function to create markers     function addMarker(i, markLat, marklong, description) {         var latlng = new google.maps.LatLng(markLat, marklong);         var marker = new google.maps.Marker({             position: latlng,             map: my_map, //map object             title: POI[i].name, //tooltip title             icon: image, //set icon         });         marker_info[i].marker = marker;         marker_info[i].des = description;         // create the infoWindow...         google.maps.event.addListener(marker_info[i].marker, 'click', function () {             if (!infowindow) {                 infowindow = new google.maps.InfoWindow();             }             // setting the content             infowindow.setContent(marker_info[i].des);             // open info window             infowindow.open(my_map, marker);                  calcRoute(marker_info[i].marker);         });          }          // called when a Place is chosen from the radio button list     function openInfoWindow_radio(Place_number) {         if (infowindow != undefined)             infowindow.close();              infowindow = new google.maps.InfoWindow();         infowindow.setContent(marker_info[Place_number].des);         // open info window         infowindow.open(my_map, marker_info[Place_number].marker);         calcRoute(marker_info[Place_number].marker);          }          // called when the 'reset the map' button is pressed     // closes an open window and re-centres the map     function resetmap(startLat, startLong, startZoom) {         my_map.setCenter(new google.maps.LatLng(startLat, startLong), startZoom);         infowindow.close();         directionsDisplay.setMap(null);         directionsDisplay.setPanel(null);          }          function calcRoute(ori_marker) {         // intialise         directionsDisplay.setMap(my_map);         directionsDisplay.setPanel(document.getElementById("directions_box"));         var request = {             origin: ori_marker.position,             destination: Town_Hall.position,             travelMode: google.maps.TravelMode.DRIVING,                   };         directionsService.route(request, function (result, status) {             if (status == google.maps.DirectionsStatus.OK) {                 directionsDisplay.setDirections(result);             }         });     }


Is there anywhere I can download the new USGS Topo TIFS?

Is there anywhere I can download USGS Topo TIFS? The new ones? The USGS site just gives me GeoPDF and their web service restricts the scale you can see them at. Edit : I am looking for the the new USGS Topo maps (c 2012).


rooftop modeling

I´m trying to model building rooftops based on a LiDAR dataset. I've already interpolated using Delaunay's triangulation a DEM with 1m resolution. Due to LiDAR data is not accurate in building edges, I managed to get cadastre data so I have building boundaries. Does anyone know how to build a 3D model of rooftops based on the previously mentioned inputs (DEM and cadastre data). Any help would be highly appreciated Thanks Adrian


Can Facebook help show (statistics) how people move from place of living to place of work?

Maybe some of you guys know if you or FB can help, or who? Perhaps together with some University? I got this idea that FB probably holds quite a lot of location information on people (that are using FB with location enabled) during different days of week. Using this data it should be possible to calculate a (partially representative) statistical view of travel in cities, such as Stockholm, that could help to improve traffic planning, and possibly (which is rarely discussed) how to avoid the need for travel. The latter however, I guess, would require that people's occupation is identified in some form and mapped to corresponding company locations (incl size?) in various locations. Regards, Fredde


Find Nearby GPS Navigation Route

I am developing a car pooling website. For that, I have created the database which has the different route options selected by different user. These routes are saved in the Database as Geography datatype. For data access, I used Entity Framework. In the C# code, I used DBGeography datatype for process the routes. I have a requirement that, the user can select one route, and identify the users who selected the same/nearby route which is saved in the database. How do I achieve this?


Is it possible to rotate the entire layout view, not just the data frame?

I have a project that requires the use of a CAD-like map template that is Landscape orientation, but with the title block along the right edge, with all the text in it rotated 90 degrees. Something like this: ![enter image description here](http://i.stack.imgur.com/XQ2xO.jpg) Except that the drawing, sheet #, etc. in mine are rotated the same way the Project Name, Drawing Name, etc. are rotated in that example. I already have this template created, but in Portrait orientation. I could rotate the data frame 90 degrees, render it to PDF, rotate the PDF page back to Lanscape, and it will look like the example above. But what I'd **like** to do is rotate the entire layout view so that my Portrait template is presented to me in Landscape on screen, so I can look at and manipulate the contents of the data frame in the same orientation they'll appear when printed. Rotating all the text elements and text boxes in an ArcGIS layout is a pain and usually requires lots of small adjustments to line everything up properly after it's been rotated, so I'd rather not have to create a new template and rotate everything to match the orientation above. I don't believe ArcGIS can actually do this, but let me know if I'm wrong..


How to select among overlapping linear features by create date?

I have a table which has overlapping linear data; for example:               A----------C and B---------C.        Of these two records I would only like to return the record with the latest create date; in other words, the record created last. How do I set up up my Search Cursor to produce the desired result? Really what I'm attempting to perform is dynamic segmentation depending on the create date. My end result can be `A---B` then `B---C` or if the create date of the first record (which is from `A-----C`) is newer than the create date for `B-----C`, I should just return record `A----C`.


Determining closest representation of polygon as collection of points+radii

Have a bit of a tough question, which I realize there may not be an easy solution to. Figured what better place to ask. I'm attempting to transform a polygon into a series of arbitrary points with arbitrary radii (circles!) that best represent the area covered by the polygon. Minor over- and underflow is acceptable, as is circle-overlap, in hopes of achieving an efficient solution (i.e. fewest possible points, no gaps). Ideally, a given poly would be represented by a few large circles, and several smaller circle on the perimeters. Essentially, the problem is that I have a dynamic number of polys that get hit on geo-spatial queries given a specific gps coordinate, however, we are required to move to a system wherein I will not be able to utilize a point- within-poly query, but will have to rely on point-within-distance queries. Hopefully someone has at least attempted something similar, and, if not, hopefully someone is willing to throw some ideas around! Open to most languages, but this needs to be done programmatically! * * * Update per clarification in comments: My point is that I feed indexes to the system and I get alerted when they are hit. I don't perform any queries on this system myself (black box), so I wouldn't have control enough to negate a query. That's the reason I need to transform the poly to a representation of points.


Center of a bounding box

I have bounding box coordinates of my shapefile   (fetched from `layer.GetExtent()` Python with GDAL/OGR):               (-38.6102467739459, -38.017601026186576, 33.01563382506067, 33.624945228091406)      So I assume it's SW and NE points of the rectangle. I would like to find out the center of this rectangle (in lat/lon).   How should I proceed?


find latlngs that are some x distance away from a point to south to north to east to west

I have a latlng i want to find four points with equal distance from this latlangs in all directions east, west, south an north is there any formula . I am trying to build polygons around higways so draw polygon around poly line the distance will be less than 1 km.


GeoServer Web Admin Interface: 404 error

I'm working through the OpenGeo workshop and I'm stuck with setting up a GeoServer. When I try to access the Web Admin Interface `localhost:8080/geoserver/`, Tomcat reports a "HTTP 404 Status" error page. The same holds for the GeoWebCache configuration page. Some additional info:   * My system: Ubuntu 13.04 with OpenGeo Suite 3.0.2 installed as described in the official reference. Things work fine on windows by the way.   * I started Tomcat 6.0.35 via `$/etc/init.d/tomcat6 start` as the installation details point out.   * The log-file of tomcat at `/var/log/tomcat6` reports a "severe" FileNotFoundException for the file masterpw.xml:              SEVERE: Exception sending context initialized event to listener instance of class org.geoserver.platform.GeoServerContextLoaderListener     org.springframework.beans.factory.BeanCreationException:     Error occured reading security configuration;     nested exception is java.io.FileNotFoundException:     /usr/share/opengeo-suite-data/geoserver_data/security/masterpw.xml (No such file or directory)      At the end of geoserver.log in `/var/log/geoserver` the following warning is reported:              WARN [support.DisposableBeanAdapter] - Invocation of destroy method failed on bean with name 'geoServerLoader':     org.springframework.beans.factory.BeanCreationNotAllowedException:     Error creating bean with name 'GWCLifeCycleHandler':     Singleton bean creation not allowed while the singletons of this factory are in destruction (Do not request a bean from a BeanFactory in a destroy method implementation!)      I'm new to web mapping and the above steps together with an internet recherche didn't resolve my problem. Do you have any suggestions? * * * **Update:** Setting up a GeoServer outside of the OpenGeo framework works just fine. The key is that there is an OS-independent binary provided on GeoServer.org using a webapp bundled in Jetty instead of Apache Tomcat. The steps it took were:   1. Download and unzip `geoserver-2.3.2-bin.zip`.   2. Move the folder to your preferred installation path and cd into it.   3. Add the environmental variable JAVA_HOME to your PATH:              export JAVA_HOME=/usr/lib/jvm/java-1.6.0-openjdk-amd64/jre      **Note** : This path is distribution dependent (here: Ubuntu 13.04)   4. Shut down the Apache Tomcat servlet engine:              sudo service tomcat6 stop        5. Start up the Server by executing the provided java archive in the folder:              java -jar start.jar        6. The Web Admin Interface is now available via web browser under              localhost:8080/geoserver      However, I suspect that if one knows about how to assign port numbers etc it might still work with Apache as installed by the OpenGeo Suite. Suggestions are welcome.


Field Value Mapping QGIS

I am wondering if there is a way to do field value maps that are filtered based on another attribute field? I know how to associate value maps, but not that are connected to another field. I have two fields in a table Type and Subtype. Type will have three values: Trees, Bushes, and Flowers. The Subtypes will have different values. ie:                Subtype for trees would be conifer and deciduous      Subtype for flowers would be roses, daisies, and sunflowers      Subtype for bushes would be blueberry, raspberry, and blackberry      What I want is when Trees is selected from Type, only conifer and deciduous are available from subtype; etc... Thanks...


To select one stream segment at a time for processing

I have attached my stream order network image. I have five different order streams. There are several segments in each order. Now my problem is to select one segment of stream at a time for the further processing. I think this is not a big problem. But I am in trouble, I m using ArcGIS 10 and Python for processing. It would be great if you give some guidance to select each segment of stream at a time using python and GIS10. ![enter image description here](http://i.stack.imgur.com/ebL2I.jpg)


Vector Mapping for iOS?

Are there any vector map data providers with API's, that i can use to build my app on iOS natively (not on a webview). Specifically im looking for:   1. For a library/framework that uses existing map servers as datasources to get vector maps. Even a HTTP based API that returns vector map data would be helpful!   2. Any tools that would help me render (on openGL) such a map on iOS. Mapnik (http://mapnik.org/) looks good, but im still unsure on how i would go about using it!  Things I've looked and why those wont work for me (tell me otherwise!):   1. Google Maps & Apple Maps : Cannot use as a datasource for vector map format.   2. Mapbox : Coming out with map vector support. Not yet out there.   3. TileStache : Looks promising. Would like more information on this if anyone has used it?


How to dynamically change SLD style of WMS layer being served by GeoServer from PostGIS?

How can I dynamically change a SLD style of a WMS layer being served by GeoServer from PostGIS?


Geocode data base of addresses using google maps api

I have a table with over a thousand addresses in postgresql, I want to make a java class that returns me the longitude and latitude of the addresses to write to the database. Does anyone have an example of how to do this? I will use the google maps api.


How to surround a Polygon object with a corridor of specified width?

I'd like to build a corridor surrounding a polygon, so that the corridor's outer bounds all have the same distance to the inner bounds. Like this: ![enter image description here](http://i.stack.imgur.com/gjXjU.png) I've made this lousy graphic in Paint, since I don't even know how this operation is called. I hope you get the idea. To avoid confusion and misplaced suggestions: the outer polygon is NOT a scaled version of the inner. I know I need some transformation for this, I just don't know what kind of. I'm implementing it in Python, I thought of shapely library, but I'm not sure if that's the one to do it, less so which method might be useful. Thanks in advance


Tasseled Cap in ArcGIS for Desktop 10.1 on Landsat TM 4-5 .tif files

I am attempting to conduct a tasscap Transformation via "Raster Property Functions" on a composit file that I have created from unprocessed, raw data files (meaning tiff-files Bands 1-5 and Band7) via the "composit bands"-tool from Landsat TM 4-5. The error message occuring when I try to apply the function is as follows: "The function could not be inserted. The Input or Output arguments are invalid." Thanking for any helpfull comment!


Loading QGIS and files from Python's use of os.system?

I am trying to open QGIS and load two files using the following code in Python:               os.system('qgis D:\data\vector.shp D:\data\raster.tif')      However it returns this error message:               'qgis' is not recognized as an internal or external command, operable program or batch file.      This line works when I enter it from the command line (withouth the 'os.system' bit) - it opens QGIS and loads the two passed files. I can load other commands such as 'gdal_merge' using the 'os.system' approach with Python (Python 2.7 in Spyder in case it's relevant). Am I missing certain bindings? I couldn't find anything online to point me in the right direction. All I can find is using the Python console in QGIS which I amn't trying to do. Thanks, Deirdre


Simultaneously adding multiple search results from Catalog into ArcMap

System spec: ArcInfo 10.1 I have performed a search on an indexed file geodatabase, looking for feature classes with the word "building" in their name. I either need to add all of these to a blank map, or copy/export them all to a new location. I was hoping that after performing the search I would be able to select all of the returned results and add them to the map, but it appears that I have to add them one at a time. There's thousands of feature classes returned, so adding them one by one isn't really an option. Am I missing something really simple? Or is there a better way to do this?


Google Maps Polygon to WKT

I need to convert a Google Maps polygon object (can has multi paths and holes) to wkt string (multipolygon) at client. Then i'll able to save polygon as geography in Sql Server 2012. I tried Wicket (https://github.com/arthur-e/Wicket) but it has a bug converting multipolygons. Is there any ready to use library? Any ideas, please? Thanks.


Exploring ArcObjects through C# or VB from Python background

While I understand that this may be off topic to an extent, I am in a bit of a bind and could use some help directing my thought process and research. I have little programming background aside from:   1) Scripting in a language similar to C   2) Practical use of Python for use in geoprocessing with ArcGIS This has worked just fine for me so far, but there appears to be a need for me to move into a .NET-based language so that I can start making use of the tools exposed in ArcObjects that ESRI has not yet made available through Python. A decision is being made by my management (with some direction from me) on whether the company should move forward and create a standard for coding in ArcObjects that is aimed at using VB or C#. The main considerations (since I am the only person who will be taking part in this) are future extendability, "safety" of code, and readability. Most of these concerns appear to be things that can be dealt with in the writing of the code as opposed to the language itself, but I don't know enough to be certain. I've read an article or two comparing C# to VB in terms of .NET, but I would like to ask for further reasoning to consider when deciding on one language over another for a person of my background. What are the most important things to keep in mind when choosing which of the two languages to begin study in for use with ArcObjects? Sorry if this is too verbose! I simply want my need to be understood. Thanks. **EDIT** : Responses have been great. I went with Dan's answer since it covered all the ground I asked about and went the extra step to point me at IPython.


Display a raster stack with several bands in ArcGIS?

I have already made a raster stack (5 bands in one TIFF) of my multitemporal Landsat images in R and have also clipped a smaller region of every raster stack for each image. Now I want to do some classification not in R but in ArcMap. But I am having problems just displaying the data. As far as I understand I have to create a raster catalog with different raster datasets... each raster dataset shall contain a raster stack for one year at one point in time. How can I make ArcGIS display the bands correctly, e.g. for the band combination RGB=4,3,2?


European Mountains Shapefile Data Sources

I am trying to locate a GIS freely downloadable shapefile (preferably polygon) of European main mountain chains (not an elevation data), eg Alps, Carpathian Mountains. I tried variuos sources , including Natural Earth, EEA, and many others, but couldn't find the right data. Can anybody help?? Many thanks in advance, Magda


Unexpected output when using gdal_rasterize

Ok, I used the following code for converting my layer into raster format but I get a tif file that is completely black or white...               from osgeo import gdal, ogr          # Define pixel_size and NoData value of new raster     pixel_size = 1     NoData_value = 100          # Filename of input OGR file     # vector_fn = '/Users/PeterVanvoorden/Desktop/GRBgis_322/Shapefile/Adp322.shp'     vector_fn = '/Users/PeterVanvoorden/Desktop/GRBgis_322/New_Shapefiles/Sbne40179.shp'          # Filename of the raster Tiff that will be created     raster_fn = '/Users/PeterVanvoorden/Desktop/test.tif'          # Open the data source and read in the extent     source_ds = ogr.Open(vector_fn)     source_layer = source_ds.GetLayer()     x_min, x_max, y_min, y_max = source_layer.GetExtent()     print x_min, x_max, y_min, y_max          # Create the destination data source     x_res = int((x_max - x_min) / pixel_size)     y_res = int((y_max - y_min) / pixel_size)     target_ds = gdal.GetDriverByName('GTiff').Create(raster_fn, x_res, y_res, 1, gdal.GDT_Byte)     target_ds.SetGeoTransform((x_min, pixel_size, 0, y_max, 0, -pixel_size))     band = target_ds.GetRasterBand(1)     band.SetNoDataValue(NoData_value)          # Rasterize     gdal.RasterizeLayer(target_ds, [1], source_layer, burn_values=[0])      STARTING QUESTION: I am trying to rasterize vector data using the GDAL_rasterize function. (Started here: syntax error when using gdal in python console QGIS) But with all methods they told me were possible, I always get something I don't need. Below the 3 methods I tried... * * * Directly in python IDLE:               >>> import os     >>> os.system('gdal_rasterize -tr 7999.8 9999.53 -te 170000.2 168000.27 177999.8 177999.8 -a Pixelkost /Users/PeterVanvoorden/Desktop/GRBgis_322/Shapefile/Adp322.shp /Users/PeterVanvoorden/Desktop/Test.tif')      Than I get the value: 32512 as result but it should output a tif file on my desktop... * * * When I do the same in a python file and run it through IDLE:               import os          shapefile = '/Users/PeterVanvoorden/Desktop/GRBgis_322/Shapefile/Adp322.shp'             (shapefilefilepath, shapefilename) = os.path.split(shapefile)     (shapefileshortname, extension) = os.path.splitext(shapefilename)                      outrastername = shapefileshortname + '.tif'     outraster = '/Users/PeterVanvoorden/Desktop/' + outrastername          os.system('gdal_rasterize -tr 7999.8 9999.53 -te 170000.2 168000.27 177999.8 177999.8 -a Pixelkost' +  shapefile + ' ' + outraster)       In this case nothing happens at all. * * * Running it in terminal (OSX):               mbp-van-peter:~ PeterVanvoorden$ gdal_rasterize -tr 7999.8 9999.53 -te 170000.2 168000.27 177999.8 177999.8 -a Pixelkost /Users/PeterVanvoorden/Desktop/GRBgis_322/Shapefile/Adp322.shp /Users/PeterVanvoorden/Desktop/Test.tif      Than I get: `-bash: gdal_rasterize: command not found` * * * I assume, looking at the last method, that there is a problem with my gdal package? I reinstalled it after having a few problems and installed GDAL 1.10 using my terminal just with a easy install like this: sudo easy_install /Users/PeterVanvoorden/Downloads/GDAL-1.10.0 Any help is appreciated!


Create basemap tiles for a custom shape (for reducing the tiles count)

in order to create a world wide tiles (for a basemap) in acceaptable zoom level one must have a huge hd (see in this post). but, most of the world is sea... so i want to create x levels of the whole world (x is small say 5-10) and then y levels of only the mainland (y is big say 20) and by that save allot of cpu time for the tile creatiion and hd space for saving unneeded tiles. i use tilemill to create tiles and you can only define a bounding box rectangle. can anyone recommend a solution (or maybe a different approach?) thanks! edit: the tiles are created from osm data


Long Time for conversion raster file to vector in QGIS

I'm a new user of QGIS, and I'm tring to convert a .tif file in a vector. To do that I have used the "Raster conversion" option, but it's taking long time. I'm waiting for the conversion since 20 minutes. This .tif file is of 1 142 007 Kb. So my question is, is it normal? Is there something that i could do?


Labeling features with font dependent on polygon area?

i want to label my island data with variety of size of text based on the island area for example island that have more than 100km, i want to label it using Arial 12 while island less than 100km i want to label using Arial 10..im label it using maplex.


downloading waypoints from Garmin GPSmap 60cx to QGIS2.2. Am I missing a step?

QGIS2.2 64 on W7 64 is failing to download waypoints from Garmin GPSMap 60cx USB port. GPSBabel installed and the GPS Tools dialog is OK: GPS device: “Garmin serial”, Port: ”USB”, Feature type “Waypoints”, Layer name is specified as is GPX output file. But when press OK the layer remains blank and its attribute table unpopulated, however the table does have the GPS field headings so QGIS is clearly recognising the device. I know that the Garmin GPS, its driver and port connection are all OK because all the Waypoints download just fine into ManifoldGIS. Am I missing a step somewhere? I don't really want to have both 32bit and 64bit instances of QGIS but thanks for the observation. Further testing: I've tried changing the Garmin "Serial data format" from "GARMIN" (default) to "NMEA In/NMEA out" but same problem. I have noticed that under "GARMIN" there is a little message saying " _The Garmin serial data format is disabled when USB is connected_ ". Is this significant for QGIS? Clearly no problem for the ManifoldGIS GPS console though. The other thing is that QGIS becomes unresponsive and crashes after a while and indeed immediately if you try to refresh the connection. Is this a bug? I know I can download the waypoints into a separate program, create a gpx or shp and then import into QGIS but this is a nuisance and awkward in the field as well as negating any possibility of live tracking. One would hope that in 2014 with the ubiquity of GPS devices one wouldn't have to do that. I'm prepared to upgrade to Garmin 62 or 64 (different protocols?) if anyone has had success connected either of these models to QGIS2.2 64bit.


Why do simplify polygon results differ between ArcGIS and PostGIS?

I am mainly interested in why and how the Simplify Polygons tool in ArcGIS handles the simplification process, including preventing and resolving topological errors better than the equivalent ST_SimplifyPreserveTopology in PostGIS destroys topology when trying to simplify to similar levels. Use case scenario: PostGIS database queried for vectorr GeoJSON polygon overlay in Leaflet. One can pass the ST_SimplifyPreserveTopology to the geometry field but my tests have shown that preprocessing data in ArcGIS gives me smaller datasets and topological correctness. Here is an example: ## Case 1 - Original data ![enter image description here](http://i.stack.imgur.com/Lw5CP.png) ## Case 2 - `ST_SimplifyPreserveTopology(the_geom,0.001)the_geom` ![enter image description here](http://i.stack.imgur.com/vQicT.png) note: querying from CartoDB Further, when simplifying Case 1's data in ArcGIS, there are no topology errors produced (ArcGIS _can_ resolves all errors).


Do pop-ups really work in Flexviewer 2.4?

I am trying to configure pop-ups using this code:               <layer label="Lansdslide Layers" type="dynamic" visible="true" alpha=".75"                     url="http://yoururl/ArcGIS/rest/services/Hazards/SR2/MapServer/"/>                     <sublayer id="1" popupconfig="popups/PopUp_Historic_Landslide_Points.xml"/>                     <sublayer id="2" popupconfig="popups/PopUp_Scarps.xml"/>                     <sublayer id="3" popupconfig="popups/PopUp_Scarps_Flanks.xml"/>      that call simple xmls like PopUp_Demographics_BlockGroups.xml found in the sample flexviewer, but with the fields named appropriately for my layers I've also tried following this esri forum with closing the sublayer element and get the error:               http response status: 200      Fault code: Client.CouldNotDecode Fault info: Error #1085 Fault details: null Does anyone have a url that shows pop-ups really working?


Attributes that characterize a geodetic track

I need a way to index and retrieve a large amount of track data. To do this I need to derive from every track as many attributes as I can. For now, I derive   * Name   * Tortuosity   * Speed (average, max)    * Covered Distance    * Length (in time)   * Bound (min lat/long and max lat/long)   * Trip mode (car, bus, plane)   * Presence of tunnel (estimated length) But I need extra ideas, spread your imagination ^_^


ArcGIS Flex API: Loop thru map graphics and make them visible

Using the flex viewer 2.4 and Flex 4.5 SDK. I have a bunch of map graphics stored in a map layer. Initially, each graphic is set to visible = false. At the click of a button, I want to loop over each graphic, set it to visible = true, wait three seconds, then onto the next graphic and do the same. Problem: the graphics don't become visible in the loop - they become visible all at once after the button click event handler has returned. I've tried doing a refresh() on the graphics and the graphic layer in the loop, but that does not help. Any ideas? Thanks.


Loop through records

I have written a python script that is supposed to basically loop through selected records of a feature class. Here's a sample of what a table looks like:               OBID   STID FL  TL           1   1   0   9            2   1   8   19           3   2   3   23           4   3   0   0            5   1   21  29           6   3   0   11           7   1   27  31           8   3   13  19      Here's what I'm trying to do:   1. Create a SearchCursor (for main file) and InsertCursor (for new file)   2. Sort all records by STID, FL and TL in order with the SearchCursor    3. Record TL for the first row and move onto next row   4. Record STID and FL in a variable   5. Compare current FL to previous TL   6. If FL<=TL, export both rows to new table with InsertCursor    7. Continue until current STID is equal to previous STID. If not, go back to step 3 (ie, repeat the process for new STID) And here's the code I wrote which is not working:               import arcpy     arcpy.env.workspace = "c:\\temp\\test.mdb"     abc = arcpy.SearchCursor("road","","","STID;FL;TL","STID A;FL A;TL A")     first = 1     for row in abc:          if first != 1:              spid = row.getValue("STID")              frl = row.getValue("FL")                  if spid == prev:                      if frl <= trl:                         print spid, frl                         prev = row.getValue("STID")                         trl = row.getValue("TL")                  else:                       prev = row.getValue("STID")                       trl = row.getValue("TL")           else:              prev = row.getValue("STID")              trl = row.getValue("TL")              first = 0     del row     del abc


Where can I deploy a JavaScript app for collaborative testing?

I'm developing a new ArcGIS Server JavaScript application, and I'd like to show the progress to the client, so he can provide feedback as I develop it. The application consists of an Index.html file, a *.JS file containing my code, a *.CSS file for the styles and a *.JSON file containing the map configuration. I don't have access to a public-facing web server in my organisation. **Are there any sites (preferably free) where I can upload my files, so he can run and test the application?** I looked at github but from what I can gather it's a _repository_ , rather than a place where my application will actually _run_ (eg, test the Identify tool and verify that it works) - please let me know if that's not correct. Thanks


How to clean vector segment intersect errors in GRASS

I'm using GRASS to process a vector map. When I export back to .SHP and do a 'Check Geometry Validity' in Quantum GIS it lists several errors of this type...               3,segments 199 and 203 of line 0 intersect at 129334.076193, 669460.784787      I've tried using v.clean to get rid of them with no effect. Can anyone give me some hints for the correct v.clean tool to use please? Also is there a vector geometry checker in GRASS so I dont have to keep exporting back to shapefile to check for errors? Here are the commands I run. The errors start appearing in 'pass2' The map is free of errors prior to import into GRASS (as far as the Quantum GIS checker is concerned at least).               v.in.ogr -o dsn=E:/geography/uk2/gb_pass0.shp output=pass0     v.clean input=pass0 output=pass1 --o tool=rmarea,rmdac,rmdupl,rmline thresh=100000     v.generalize input=pass1 output=pass2 --o method=chaiken threshold=100     v.generalize input=pass2 output=pass3 --o type=boundary method=reduction threshold=100 -r     v.generalize input=pass3 output=pass4 --o type=boundary method=douglas threshold=40 -r     v.clean input=pass4 output=pass5 --o tool=prune thresh=40     v.clean input=pass5 output=pass6 --o type=area tool=rmarea thresh=7000000     v.clean input=pass6 output=pass7 --o tool=rmdangle,break,rmdupl thresh=1000     v.clean input=pass7 output=pass8 --o tool=rmbridge,rmsa,rmdac,bpol,rmdupl     v.clean input=pass8 output=pass9 --o tool=break,rmdupl,rmline     v.centroids input=pass9 output=pass10 --o     v.out.ogr input=pass10@mapset1 type=area layer=1 format=ESRI_Shapefile dsn=E:/geography/uk2/gb_boundary olayer=gb_boundary -e


working with multiple wfs-t layers

I am using multiple wfs layers and using editing tools in that like edit,draw,save etc. The code is shown here......               var district = new OpenLayers.Layer.Vector("district(wfs)", {                         strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],                         //projection: new OpenLayers.Projection("EPSG:26910"),                     protocol: new OpenLayers.Protocol.WFS({                     version: "1.1.0",                     // loading data through localhost url path                     url: .../geoserver/wfs",                     featureNS :  "http://www.openplans.org/topp",                     //maxExtent: mapextent,                     // layer name                     featureType: "district",                     // geometry column name                     geometryName: "the_geom",                     schema: ..../geoserver/wfs/DescribeFeatureType?version=1.1.0&;typename=topp:district"                     })                     });                                          district.style =           OpenLayers.Util.applyDefaults({strokeColor: "#0000ff"},                                          OpenLayers.Feature.Vector.style["default"]);          var major_location = new OpenLayers.Layer.Vector("major_location", {                         strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],                         //projection: new OpenLayers.Projection("EPSG:26910"),                     protocol: new OpenLayers.Protocol.WFS({                     version: "1.1.0",                     // loading data through localhost url path                     url: .../geoserver/wfs",                     featureNS :  "http://www.openplans.org/topp",                     //maxExtent: mapextent,                     // layer name                     featureType: "major_location",                     // geometry column name                 geometryName: "the_geom",                 schema: .../geoserver/wfs/DescribeFeatureType?version=1.1.0&;typename=topp:major_location"                 })                 });                          major_location.style = OpenLayers.Util.applyDefaults({strokeColor: "#0000ff"},                      OpenLayers.Feature.Vector.style["default"]);      Now, when i am doing editing in these layers then that editing is shown on only one layer.But,I want to edit both layers simultaneously.For editing the code that i am using is:                           var panel = new OpenLayers.Control.Panel(             {'displayClass': 'customEditingToolbar'}             );              var navigate = new OpenLayers.Control.Navigation({             title: "Pan Map"             });              var draw = new OpenLayers.Control.DrawFeature(             district, OpenLayers.Handler.Polygon,             {             title: "Draw Feature",             displayClass: "olControlDrawFeaturePolygon",             multi: false             }             );                   var edit = new OpenLayers.Control.ModifyFeature(district,{             title: "Modify Feature",             displayClass: "olControlModifyFeature"             });                   var del = new DeleteFeature(district, {title: "Delete Feature"});              var save = new OpenLayers.Control.Button({             title: "Save Changes",             trigger: function() {           if(edit.feature) {             edit.selectControl.unselectAll();             }             saveStrategy.save();             alert('saved');             },             displayClass: "olControlSaveFeatures"             });                       panel.addControls([navigate, save, del, edit, draw]);             panel.defaultControl = navigate;             map.addControl(panel);      So,please tell me is it possible to edit multiple layers?


Filter by attribute to show only records with whole numbers

Is there a way to show points that only contain whole numbers in a column (eg. 009, 010, 011 and NOT 9,1; 9,2; 9,3...) in ArcGIS 10?


Why would a map lost the coordinate system when I use it in ArcGis?

I have a file in AutoCad with ETRS89 coordinates, and when I use the same map in ArcGIS it loses coordinates. With other files do not happen.


How to find lat/lon values for every pixel in a GeoTIFF file?

I want to a store GeoTIFF file pixel by pixel in the database. I need to store pixel values along with the lat/lon of that pixel on the GeoTIFF image. Is it possible to do that?


Is there any way to view maps in "perspective" with ArcGIS (ArcMap do not)?

I need to work with ArcGis by the amount of data that I use and the answers I need, but would like to see my maps in perspective and I can't find a way. Does anyone knows how to do it?


Rasterize vector with gdal

_"I would like to calculate the Euclidian distance from a residence to a public service of interest e.g. the distance from a residence to a green space (as the crow flies) in QGIS. After browsing the Internet it seems that I need to rasterize my vector before I can use the Proximity tool provided by GDAL. When I try to make the conversion from vector to raster my thought was to take the shp.file which contains the polygons. The function "rasterize" tells me to chose an "output file for rasterized vectors (raster)" but I have no idea what to chose here. Should I create a raster for this? And how do I this? I'm using QGIS 1.8 Mac and GDAL 1.9. Thanks in advance."_ Ok, it seems like I have not done my homework good enough. Let me start over: I have a vector layer with points from which I want to calculate the distance to vector layers with lines or polygons. Still I want to calculate the Euclidian distance (as the crow flies) from the points in the first layer to the lines or polygons in the other layers: The distance from a residence (point) to a stream (line) or forest (polygon). Additionally, I would like to put the distances in the attribute table of the vector layer with points in a new column. Sorry for the inconvenience.


Batch Clip (Editing Session)

I'm trying to correct some topology errors in my geodatabase, now focusing on two feature classes: grass and buildings. They must not overlap. However, there are plenty of spaces where the grass was generalized, thus, covering entire buildings. There are no automatic fixes for this errors - unfortunately the subtract function deletes the overlapping part of both features, when I'd only like to delete the grass feature class and preserve the building feature class. In this way, I'm using the "Clip..." function on the Editor Menu to use the buildings to cut out this portion of the grass. However, there's no way to select more than one building at a time to apply this function. I cannot use the Erase Function from ArcToolBox, because the same inconsistence can happen also in the other way: generalized buildings covering grass areas. Also, I wouldn't like to have an output feature class, once I just need to change the shape of my features, nothing else. **Long story short, is there any way I could use the "Clip..." function from the Editor Menu with multiple polygons, by selecting them all at a time?**


Well-known text (wkt) validator

Is there anything like http://jsonlint.com or http://geojsonlint.com for WKT? How can I validate a complex WKT string, say e.g. a Multipolygon?


Visualization of poin.shp in QGIS 2.0 with manually selected min max values

I would like to visualize a poin.shp file in QGIS 2.0 with manually selected min max boundaries. The data consists of many different points lined up as a line. The aim is to visualize the different data values from each point with a different color. The problem is that I can't find an option to cut the data range only in the visualitation and not in the data. I have something in mind like "generate new color map" for grid data, where you can define mode, classes and min/max. Any advice would be greatly appreciated, thanks!


Are unicode symbols supported in QGIS labels?

I want to use unicode symbols in a QGIS label - is this supported? I.e. I have this expression for labelling a line:               case when xat(-1) - xat(0) > 0 then       'GF ->'      else       '<- GF'     end      and want to replace the '->' and '<-' with the unicode characters representing left and right arrows. Is there a way I can achieve this?


How to draw Geodesic lines in Openlayers?

Hi i want to display multiple lines between two coordinates. i will pass no of lines and color of the line, it has to display like below image. please see the below image.this is my previous post ![enter image description here](http://i.stack.imgur.com/RNqJH.png)


What's the best way to run a visibility analysis in PostGIS 2?

I am looking for the best way to run a visibility analysis in PostGis 2, i.e. given 2 points and obstacle polygons, are these 2 points visible from each other ? I did not find anything relevant for now. Are you aware of :   * such a function   * a library I could use to perform this   * or do I have to implement it by myself ? Thanks in advance


Openlayers infobox?

is there any infobox available in Openlayers same like google maps infobox google maps infobox


Calculating area gets a different result

I've calculated the area of polygons with QGIS by using `ctr-i - geometry - $area`. I've set it on real. The results seem fine, but I've noticed they differ from the area given when I use the i-function. For example: polygon 4 gives an area of 7,662 km2 when using the i-function but when I `Calculate` it, it gives an area of 7682816.27364349. Does anybody know what I'm doing wrong? Layer CRS is Amersfoort / RD new EPSG:28992


Where is the File Geodatabase Raster Dataset file stored?

I have a esri file geodatabase with a raster dataset. How can I find out which physical file the dataset is referring to? (I am retrieving the file geodatabase through the internet. It all seems to go well, but the raster dataset is missing in the file database after retrieval. Zipping the file geodatabase before "transport" solves this problem. To debug the issue I want to see which physical files are missing/different between the zipped and not zipped databases.) ![enter image description here](http://i.stack.imgur.com/MJUbp.jpg)


Problem when compiling QGIS (with customizations) for Windows

I was able to customize a few features of QGIS in Ubuntu. Now I am trying to compile the same code on Windows, by following the installation guide given for Windows. I have downloaded all of the dependencies using OsGeo4W and struck at SIP. I tried setting environment variables to point to sip site- pacakges too ( C:\OSGeo4W\apps\Python27\lib\site-packages ) Any guidance would be of lot of help. I am using visual studio compiler libraries. Here is the cmake trace               Quantum GIS version: 1.8.0 Lisboa (10800)     Could not find GRASS     Touch support disabled     Found Proj: C:/OSGeo4W/bin/proj.dll     Found GEOS: C:/OSGeo4W/bin/geos_c.dll     Found GDAL: C:/OSGeo4W/bin/gdal18.dll (1.9.2)     Found Spatialindex: C:/OSGeo4W/bin/spatialindex1.dll     Found Qwt: C:/OSGeo4W/bin/qwt5.dll (5.2.1)     Found Sqlite3: C:/OSGeo4W/bin/sqlite3.dll     Found PostgreSQL: C:/OSGeo4W/lib/libpq.lib     Found SpatiaLite: C:/OSGeo4W/bin/spatialite.dll     Pedantic compiler settings enabled     Found Python executable: C:/Python27/python.exe     Found Python version: 2.7     Found Python library: C:/Python27/libs/python27.lib     CMake Error at cmake/FindSIP.cmake:54 (MESSAGE):       Could not find SIP     Call Stack (most recent call first):       CMakeLists.txt:476 (FIND_PACKAGE)          Configuring incomplete, errors occurred!      Thank you. Srikanth


How to install GDAL 1.10 with OGR SQLite driver support on Ubuntu?

How can I install GDAL 1.10 with OGR SQLite driver support, and JPEG and PNG drivers per http://www.gdal.org/frmt_mbtiles.html


postgis values within XYZ viewport in openlayers

I got this XYZ server where I get my maptiles. Now, when the user is at zoom level 1 ( closest). I want it to execute a command that gets the current viewport from the map and Queries my PostGIS database to get the geometry points that are contained within that viewport. How could this get done? and what is the best way to display this data as fast as possible on the map? Thank you in advance :)


How can I identify gaps in polylines and polylines with more than one end point using ArcGIS 10?

I have a multipart polyline shapefile containing a large number of lines, some of which are broken by gaps and some of which have more than one start and end (by mistake - see attached). I need to have a multipart shapefile output, which contains only unbroken lines with one start and one end. Is there a way to identify these errors (can't seem to do it with the the topology rules in Arc)? Once identified they are easy to fix manually. Thanks.


Creating space at intersection of two links

I have a shapefile consisting of two links which are crossing over each other. I want to create an empty space at the point of intersection of the two links. As of now, I have been manually doing this like this:   1. Split the line at two points near the intersection. And then remove that feature.   2. Similarly carry out this for the other line. So by doing this, I create a empty space at the point of intersection. I want to know whether there is any better way to achieve this using a plugin or any other way? Can a plugin can be created specifically for this task ? (I'm ready to write one, if this is achievable.) Below is an image where one link touches other link. Now, I want to create an empty space there. ![enter image description here](http://i.stack.imgur.com/cEEDp.png) This is the output I'm expecting: ![enter image description here](http://i.stack.imgur.com/PU8Gl.png)


Adding time dimension to layer in GeoServer

I have a GeoServer installation that is used as a proxy in front of a THREDDS Data Server to be able to add additional metadata to WMS layers (which I haven't found any support for in THREDDS). I have configured a GeoServer WMS store and added the the WMS GetCapabilities URL to the THREDDS server to the new WMS store. I have then added a layer from the new WMS store to the list of layers from GeoServer. So far so good. GeoServer can now act as a proxy for the layer in THREDDS. The problem comes when I try to enable Dimensions for the new layer. I go into edit mode for the layer and then click the "Dimensions" tab. I then enable the "Time" dimension and try to select the correct "Attribute". However no attributes are found in the list and I am not allowed save the new configuration since "Attribute" is a required value. Any clue to what I can do to solve this problem is greatly appreciated. When not accessing the layers via GeoServer the time dimension works fine.


Can I use a Esri geocoding service with my custom basemap?

I am looking at the JavaScript samples in here. I am searching some street by search bar. But not zooming to street, it only viewing city. I have a custom tiled map service layer in my server. Can I use this geocoding service and my custom basemap?


How to write raster field values into attribute table using ArcPy?

Is it possible to use ArcPy (10 or 10.1 )to write raster field types into a featureclass? either SDE or filegeodatabase?. I'm hoping to use arcpy as means of automating the process of adding rasters as an attribute to a featureclass. Cheers.


adding geometries to a source, target table for multiple rows in pgrouting

I am having an issue trying to perform the pgr_dijkstra algorithm and getting back a line in my table for multiple rows. I have a table that has a columns of source, target node id's but I cannot figure out how to feed them into the pgr_dijkstra function properly. My table basically looks like this, but I need to fill in the routeline column with the shortest path line connecting the source and target. ` ID | source | target | routeline 1 | 30 | 217 | geometry(linestring) 2 | 217 | 2143 | 2 | 15 | 43 | ` I would eventually like to make a linestring out of several of the rows from the roueline column, as I am basically trying to follow survey respondents around a network to all the stop locations they visited, but I can't seem to find a pgr function for multiple ordered nodes. For now, i just figured I would connect the two nodes and then connect the lines together later. Please let me know if you have a better approach! I have been following the tutorial here : http://workshop.pgrouting.org/chapters/geoserver.html My problem is when I run my code:               select a.uniqueid, a.source, a.target , st_makeline(route.geom) as routeline      from           analysis.routing as a,          (select geom from pgr_dijkstra('ways',a.source,b.target)     order by seq) as route;      I get the following error, and can't figure out how to pass the source and target from each row into the function.               ERROR:  invalid reference to FROM-clause entry for table "a"     LINE 1: ...uting as a, (select geom from pgr_dijkstra('ways',a.source,b...                                                                  ^     HINT:  There is an entry for table "a", but it cannot be referenced from this part of the query.      Does anyone know what is going on, I can't figure out how to reference my source from the reference table. Eventually, I need to insert this into an update statement such as               update analysis.routing set     routeline = (select ***the querry above with results as linestring)


How can I add a grid (from projection MTM8) on my layers

I need to add a grid on a map I created with my own layers and GPS data. How can i do this? I know how to do it in Arcgis, but not in QuantumGIS. Can you help me?


Question on Exporting data from 'collect values' output in ArcGIS model builder

I have a question regarding an answer to "Exporting data from 'collect values' output in ArcGIS model builder". I'm trying to do the same thing as the author to that question and I've followed the steps and the code in that answer. The code I am using is below:               multivaluesToCsv(r"%CSV File%", "%Feature Class Name%", "%NNZScore Value%", "%NNExpected%", "%NNObserved%", "%PValue%", "%NNRatio%")               import csv          def multivaluesToCsv(csvfile, featureclassnames, nnzscorevalue, nnexpectedvalue, nnobservedvalue, pvaluescore, nnratiovalue):     with open(csvfile, 'wb') as f:      w = csv.writer(f)      columnNames = ["FC", "NNZ", "NNE", "NNO", "PV", "NNR"]      w.writerow(columnNames)      rows = zip(*map(lambda x: x.split(';'), [featureclassnames, nnzscorevalue, nnexpectedvalue, nnobservedvalue, pvaluescore, nnratiovalue]))      w.writerows(rows)      My problem is that the CSV file is only being populated with the final row of values. The titles all appear but only the last iteration's data is included. Any suggestions?


How can I set a layer visible dynamically using Flex Viewer?

I have a set of layers in operationallayers tag in my Flex Code connected to ArcGIS Server 10, and I want to make them visible dynamically in a ComboBox Change Handler. Is it possible? Thank you


vector tool missing when install QGIS 1.7.4 on new computer

I installed QGIS version 1.7.4 a few months ago. Under the menu Vector > Data Management Tools I have two options to make a join: 'Join Attributes' and 'Join attributes by location' (see screenshot). ![enter image description here](http://i.stack.imgur.com/SJsqv.png) My colleague installed the same version a couple of weeks ago. In here version, the option 'Join attributes' is missing. Today I installed QGIS 1.7.4 on some computers in a classroom and the same option was missing. It is a handy tool to join a shapefile with a DBF, instead of doing it through Layer properties > tab Joins. Can I find it somewhere else?


mulitlinestring to linestring

I am using postgres with postgis and I have a set of points which I would like to interpolate with a set of lines. Now my problem is that the set of lines are of type multilinestring. Is there a way to convert them to linestring?


Why won't osm2po setup table correctly?

OK, I have no idea what's going on. I'm using a Mac and successfully installed pgRouting. I used osm2po on some data, which created a sql file. I imported the sql file into a pgRouting database and all is well, except the fact I only took into account routes for cars. After Carsten helped me with configuring `osm2po.config`, I ran osm2po again that ran as expected. From the terminal, I run the command:               psql -U [username] -d [dbname] -q -f "/Users/John/Desktop/osm2po-4.6.9/wmCar/wmCar/wmCar_2po_4pgr.sql      I was expecting it to be successful like last time, as to my knowledge the conditions are the same. However I get the error:               psql:/Users/John/Desktop/osm2po-4.6.9/wmCar/wmCar_2po_4pgr.sql:101739: ERROR:  INSERT has more expressions than target columns     LINE 2: ..., -2.0046492, 52.5307149, -2.0026647, 52.5310823, >'010200002...      With pgAdmin I look at the empty table that was created, and compared to the first sql file I imported, the new table was missing a column. That being `geom_way`. I find this odd as I haven't changed anything. I threw osm2po folder to the trash, and unzipped a new osm2po folder. To test whether it was going to work, I ran osm2po on the original pbf file with the original `osm2po.config` settings and tried to import that into the same pgRouting database. The error appeared again. For some reason, the last column, `geom_way`, isn't being added to the table. Would anyone know why it's started doing this?


Calculate Field based on field value in a second table

I have a table in a MSSQL Server Enterprise SDE database. I need to join a table from a MS SQLExpress Workgroup SDE database and then calculate 11 fields in the Enterprise table based on values in the Workgroup Table then remove the join. I have tried to do this in model builder add it takes forever to run. I thought of selecting only those records that don't match before performing the calculate, but then thought maybe there were other ways to tackle this problems. So now I am thinking that there are three alternatives that will do what I need in a quicker manor. My question is can anyone give insight on which would be faster or easier to do. Option 1 Create a python code in ArcGIS 10.1 to perform the join select the mismatched records and perform the calculation then move on to the next field, when all done remove the join. Option 2 Use FME Desktop to perform the same steps outlined above. Option 3 Use SQL and perform the same steps above in Database An input or suggestions would be appreciated.


openlayers and labeling

Okay, I've been banging my head against the wall on this and I'm pretty sure I'm missing something silly. So I'm hoping you guys can help. The idea here is to allow the user to type in the label they want. It appears on the map in the center and then they can move the label to the location they want it to be. (Then I'll add functionality for styles, but basic functionality first, frills later.) The problem is that I can't move the labels. I can't select them either. In fact, I can't select anything on the layer, but it is on top (or at least it should be). And yes indeed, 'labelSelect' is set to True. Anyway, here is the code to initialize the layer and the controls:               editableLayer = new OpenLayers.Layer.Vector("User Layer");     editableLayer.setZIndex(4000);     layerArray.push(editableLayer);          map = new OpenLayers.Map({     div: "map",     layers: layerArray,     controls: [         new OpenLayers.Control.LayerSwitcher(),         new OpenLayers.Control.PanZoomBar(),         new OpenLayers.Control.Navigation(),         new OpenLayers.Control.ScaleLine({geodesic: true}),         new OpenLayers.Control.KeyboardDefaults(),         new OpenLayers.Control.MousePosition({             displayProjection: new OpenLayers.Projection("EPSG:4326"),             numDigits: 5,             emptyString: '',             prefix: "x: ",             separator: " y: "         }),         new OpenLayers.Control.EditingToolbar(editableLayer),         new OpenLayers.Control.SelectFeature(editableLayer, {toggle: true, clickout: true}),         new OpenLayers.Control.DragFeature(editableLayer),         overviewMapControl     ]     });      Here is the code for creating the label:               var labelElement = document.getElementById("LabelText");     var labeltext = labelElement.innerHTML;          var labelFeature = new OpenLayers.Feature.Vector(                 new OpenLayers.Geometry.Point(                     map.getCenter().lon,                     map.getCenter().lat                 ),                 {                     text: labeltext                 },                 {                     label: labeltext,                     labelSelect: true                 }             );     editableLayer.addFeatures([labelFeature]);      So, what did I mess up? I tried ModifyFeature with the mode set to DRAG as well, that didn't work either. I suspect the problem is somehow in my layer ordering, but even when I add code for map.setLayerIndex it doesn't help and when I inspect the elements in the browser it seems to be on top. Thanks.


How to convert a symbol style set for ArcGIS for use in QGIS?

Is there a way to convert a symbol style set for ArcGIS for use in QGIS? There are symbols used in Incident Command I would like to use in QGIS, but can only find a style set for use in ArcGIS. The style set I have is saved as a .style file.


Not able to convert GPS coordinates received from a software into Latitude and Longitude

I am working on a software which returns GPS coordinates in this format: Lat : '48684499' Lon : '5203474' I can not understand what type of coordinates are these and not able to convert them into proper Latitude and Longitude like : Lat : '-41.295929' Long : '174.77573'


Big IRowBuffer problem

I have a small software that connects to a database (SQLServer) reads a few values and transfers them to a SDE Oracle instance (11g). The problem is: I have over half million records that I need to transfer. My code reaches for the records without an issue, but when I try to insert...(my guess is that you guys already got throught this) I get a lovely COM error, screaming at me about Single Threaded and Multi Threaded apartments. Here is a snippet:                       IFeatureWorkspace workspace = GetWorkspace() as IFeatureWorkspace;             ITable newTable = workspace.CreateTable("blevers", fields, uid, null, "");                  IRowBuffer buffer = newTable.CreateRowBuffer();             ICursor cursor = newTable.Insert(false);                  IMultiuserWorkspaceEdit workspaceEdit = workspace as IMultiuserWorkspaceEdit;             IWorkspaceEdit2 workspaceEdit2 = workspaceEdit as IWorkspaceEdit2;             workspaceEdit.StartMultiuserEditing(esriMultiuserEditSessionMode.esriMESMNonVersioned);             workspaceEdit2.StartEditOperation();                  for (int i = 0; i <= stringList.Count - 1; i++)             {                 buffer.set_Value(0, stringList[i]);                 cursor.InsertRow(buffer);                 Console.WriteLine("Escrevendo linha " + i.ToString() + " de " + stringList.Count);             }             cursor.Flush();                       workspaceEdit2.StopEditOperation();             workspaceEdit2.StopEditing(true);      I usually get the error inside the for loop. How can I insert big big cursors without problems? EDIT: i've also tried using a SchemaLock, but a COM error explodes in my face when converting the ITable to ISchemaLock. I'm out of ideas!


Cache service issue

i have a map service that is cached and whe I try to access this service, I get the message displaying that downloading picture and I can see that it is trying to download and then finally get no display. When I looked in fiddler, I am getting 500 error but there are some tiles that got generated. I am not sure what is the issue. Thanks Jay ![enter image description here](http://i.stack.imgur.com/Utgmr.jpg) ![enter image description here](http://i.stack.imgur.com/x6bv8.jpg)


Migration from .gdb to .mdb - how to change automaticly all querys at MXD from " " to [ ]?

after migration from File Geodatabase (.gdb) to Personal Geodatabase (.mdb) - how to change automatically all field names at query definitions at the MXD from " " to [ ] What's else is important when conducting such transition ? I'm using version 10.1 sp1 and I have a MXD with dozens layers and dozens of queries, also at the symbol classes of the labels there are even more SQL syntax to change. Thanks in advance!


Are there any good resources on WFS support in ArcGIS? The only documents I found so far are these: How_to_work_with_WFS_using_ArcGIS_Desktop and Consommer des services/données WFS dans ArcGIS 9.2 But I want to know how it works! I'm interested in how ArcGIS can be a client to a WFS Server. Seems like the Interop conecction is downloading all data on the WFS server into a local Cache, is there any way to get it to fetch the features on- demand? **Update:** I voted for this idea Make ArcMap an effective WFS Client to get better support for WFS.


What is a seamless table?

A client has asked me to provide a dataset in seamless table format. Could anybody describe to me what a seamless table is?   I understand that it is a format used by mapinfo, but can I create a seamless table in some other way, by FME for example?   Thanks.


Self hosted, highly accurate, address geocoding?

I'm working on an app that requires Rooftop level accuracy coupled with very high volumes. The client would prefer to self-host the solution which might not be a bad idea considering our volume and response time requirements. I'm not seeing how the usual suspects (Google, MapQuest, ...) might help and the geocodes we get from address correction software are inconsistent quality- wise. Any pointers or suggestions to vendors and/or solutions are appreciated. FWIW we're currently using PostGIS.


How to perform enhancement operations on an satellite imagery in QGIS?

I am a new user of QGIS please suggest how to perform image enhancement operations on an satellite imagery like ratioing, contrast stretching, edge engancement, filtering etc. in QGIS?


How to programatically check for a mouse click in QGIS

I want to know how to check for a mouse click in QGIS. I am trying to write a python plugin and want to provide functionality similar to the "Select Single Feature" tool that already exists in QGIS. I checked the QGIS api docs and found               QgsMapCanvas::CanvasProperties::mouseButtonDown      This sounds promising. I have a QgsMapCanvas object but I can't see how to access the mouseButtonDown attribute. I am completely new to the QGIS api so any and all help is greatly appreciated. Thank you in advance.


it seems like OpenLayers.Filter.Comparison.LIKE doesn't filter special characters! for example,               filter = new OpenLayers.Filter.Comparison({          type: OpenLayers.Filter.Comparison.LIKE,          property: 'name',              value: 'blue'          })       this will work fine. but:               filter = new OpenLayers.Filter.Comparison({          type: OpenLayers.Filter.Comparison.LIKE,          property: 'name',              value: 'light-blue'          })       this will cause an error because of the - it also can't deal with dots what can I do to fix that ? **EDIT** I tried:               var filter = new OpenLayers.Filter.Comparison(       { type: OpenLayers.Filter.Comparison.LIKE,         property: 'name',         value: input }     );     filter.value2regex( "-","." );       I wanted to escape those two characters from any value. but nothing changes.


Viewing relationship classes ties for multiple features

I've been working with a complex ArcGIS 9.3.1 geodatabase with multiple feature classes (FC's) interconnected by many relationship classes. My question: supposing I have a FC with features a,b,c and another FC with features 1,2,3 - **I want to see which featuers of the first FC are to connected to which features of the second FC** (of course in reality, the geodatabase is much more complex with more FC's, more features and more relations). Now, I know that the **IDENTIFY** tool does it, but it does it to one feature at a time, and I can't export the ties to an outer file (it is only visual). Is there any tool I can use to do that for whole layers at a time? (I've seen Geodatabase Diagrammer very briefly, but as far as I understood, it is useful for the planning of the structure of the geodatabase, but not it's content. correct me if I'm wrong).


Merging Issues doubling up

I have two shapefiles with some of the parcels (surveying) overlapping each other. How do you get those overlapping parcels from doubling up? Is there a way you can have it dissolved or union etc... so that second parcel with the same "tax id number" from being included in the merge? Thank You.


Why does GeoServer copy my PostGIS tables?

I have some geospatial data in a PostGIS database. I want to serve that data to a web map app with GeoServer. But every time I add a table as a layer in GeoServer, it creates a duplicate table in PostGIS. For example, when I try to add the table "district_boundary", I end up with a new table, named "district_boundary0" that contains all the data from the original table, duplicated for no apparent reason. I don't want this! There's no reason to have multiple copies of the same data - and there are lots of reasons not to! _So how can I make GeoServer use the table I tell it to and not make its own?_ Thanks in advance! _Note:_ This is probably not a matter of my tables not having the correct structure to be served up. I added the GeoServer-created "district_boundary0" to check this, and just ended up with "district_boundary00"...


animation of polylines in ArcGIS 10

I would like to think with all the animations out there thsi would be relatively simple but it hasn't been. I would like to animate the movement of a stream (polyline) from one end to another, with all its forks and tributaries (smaller polylines). I have been reduced to trying the Time Stamp property in ArcGIS 10. I have a number of questions related to this, such as why is it that the Time Stamp property always defaults to hours, or that if I try to alter the time property to a different time step, it causes the one segment I have been working with to disappear? If anyone has worked with the time stamp property in ArcGIS 10 successfully for a whole layer I would love to hear from you. Dan Brenner Eureka, CA


Where can the GDAL executables be found in a QGIS installation?

I run QGIS 1.8.0 Lisboa on a Windows 7 machine. Yesterday I attempted to install the TauDEM plugin using these instructions: http://lists.osgeo.org/pipermail/qgis-developer/2012-June/020807.html Since doing this the GDAL tools for QGIS do not work with error: > The process failed to start. Either the invoked program is missing, or you > may have insufficient permissions to invoke the program. I understand this means they cannot find the GDAL executables. However, neither can I. Where actually are the GDAL executables in the QGIS folder? I assume I then have to add this manually to the Environmental Variables or into the GDAL tool settings in QGIS - is this correct? PS. I can run GDAL tools through the OSGEO4W Shell so GDAL itself is alive and well.


spatialite sql and possibly using case

i have a database that has got tables with images and some of the tiles for the images is repeated and i want to eliminate the duplicates based on which imagery is oldest so as to remain with the newest images, i have run the following to get the duplicates from the database but do not know how to delete the row based on the age of the imagery               select "index",path from ngi_catalog_aug2012 where "index" in      (SELECT "index"      FROM ngi_catalog_aug2012     GROUP BY "index"     HAVING ( COUNT("index") > 1 )     ) order by "index"      this gives results like               index       path          2228BD_23    /data/ngi_old/images_30k_sid/2228/2228BD_23_498_469_12_0007_PAN_RECT.sid          2228BD_23   /data/ngi_old/images_30k_sid/2228BD/2228BD_23_498_469_11_0100_PAN_RECT.sid      and i would like to delete for instance the second row


OpenLayers or OGCServer request optimization

I have an OpenLayers-OGCServer-Mapnik setup through a WMS layer. I noticed when i use the openlayers UI in a quick manner (fast zoom with mouse scroll, a couple of consecutive pans with mouse) the response from WMS slows down, because its obviously processing all the tiles between the previous "static" state and the last one. Is there a switch in openlayers or perhaps mapnik ogcserver to optimize this behavior by either postponing requests until the user pauses for a moment on the same location (clientside) or canceling previous requests from the same IP on the OGCserver (serverside)


Problem transforming projection on outlying features from shapefile

I am fairly new to GIS data and as such please feel free to ask for extra data I might not have provided. I have some regional data from Natural Earth, more specifically this: http://www.naturalearthdata.com/downloads/10m-cultural- vectors/10m-admin-1-states-provinces/ (States and Provinces, first download). This uses a longlat-projection, which I wish to turn into a miller projection. However, several of the regions get "mangled" by the transformation. It seems to only affect those sitting at the border of the map (far west / far east), more specifically: (In the format "Country - Region name") Antarctica - ATA-00 (Antarctica aggregation) Brazil - Goiás Fiji - Northern Russia - Chukchi Autonomous Okrug Zambia - Eastern If I manually go in and edit the nodes touching the edge, it seems to fix the problem. However, I am trying to automate the process, so I can automatically generate new data from the shapefile when a new release is made. Is there any smart way of overcoming this issue? The error presents itself as if the feature is trying to combine all the polygons into one, drawing ugly lines across the entire map to connect them. An example can be seen here: http://i.imgur.com/bzJ5V.png (taken from QGIS). The result is the same whether through a script (using GDAL) or through the QGIS program.


Fetching field values from uncommitted features on Qgis

I'm trying to fetch a feature's field value on Qgis, I use      `QgsFeature.attributeMap().value(cont int & key).toDouble()`      This works fine for committed features but it always returns " 0 " if the feature is not yet committed to the vector layer !!


Connnect to PostGIS db using QGIS - when not on localhost

I"m just getting started with PostgreSQL/PostGIS db (v9.2) and I'm wanting to connect to db while not on the localhost using QGIS. I've been able to connect successfully on localhost, however I'm not sure what I need to do to make this non localhost connection. I'm able to ping the db server from my desktop. In pgAdmin within the db properties I defined the following privileges: ![enter image description here](http://i.stack.imgur.com/O8V7p.jpg) Then in QGIS here are my connection settings (using postgres username/password, but logged in as my domain user): ![enter image description here](http://i.stack.imgur.com/e9mvj.jpg) I've tried searching the web for this, however all the hits are just coming up with examples of how to connect on localhost. I'm thinking maybe I have to add my domain info as user login? Thanks