title: Implement 'Zoom To Selected Features' in VB.Net

text: Seems like this should be dead simple. I'm writing an add-in for ArcMap 10 in VB.Net. I need code that will reproduce the 'Selection --> Zoom To Selected Feature' menu option.


title: How do I serialize/deserialize .NET dictionary without key value elements?

text: I'm trying to get my C# class to serialize to a json that complies with GP Job, for example:                {           "jobId" : "jb10b153b77e74a839e31af53e5d4d0d1",            "jobStatus" : "esriJobSucceeded",            "results" : {             "output_Buffer" : {"paramUrl" : "results/output_Buffer"},              "Final_wellpoints" : {"paramUrl" : "results/Final_wellpoints"}           },            "inputs" : {             "Input_Features" : {"paramUrl" : "inputs/Input_Features"},              "Distance" : {"paramUrl" : "inputs/Distance"}           },            "messages" : [             {"type" : "esriJobMessageTypeInformative", "description" : "Submitted."},              {"type" : "esriJobMessageTypeInformative", "description" : "Executing..."},              {"type" : "esriJobMessageTypeInformative", "description" : "Succeeded."}           ]         }      But my dictionaries (results and inputs) serialize with key/value cruft:               //     {             "jobId":"jb10b153b77e74a839e31af53e5d4d0d1",              "jobStatus":"esriJobSucceeded",              "results":[                     {"Key":"output_Buffer",     "Value":      {"paramUrl":"results/output_Buffer"}},                     {"Key":"Final_wellpoints",     "Value":        {"paramUrl":"results/Final_wellpoints"}}],              "inputs":[                     {"Key":"Input_Features",     "Value":        {"paramUrl":"inputs/Input_Features"}},                     {"Key":"Distance",     "Value":        {"paramUrl":"inputs/Distance"}}],              "messages":[                     {"type":"esriJobMessageTypeInformative",     "description":"Submitted."},                     {"type":"esriJobMessageTypeInformative",     "description":"Executing..."},                     {"type":"esriJobMessageTypeInformative",     "description":"Succeeded."}             ]     }      **Update:** I've commented out the WCF code and replaced with Json.NET call.               [DataContract]     public class GPJob     {         public static string Test()         {                  Dictionary<string,taskParam> resDict = new Dictionary<string,taskParam>();             resDict.Add("output_Buffer",new taskParam(){ paramUrl = "results/output_Buffer"});             resDict.Add("Final_wellpoints",new taskParam(){ paramUrl = "results/Final_wellpoints"});             Dictionary<string,taskParam> inputDict = new Dictionary<string,taskParam>();             inputDict.Add("Input_Features", new taskParam(){paramUrl="inputs/Input_Features"});             inputDict.Add("Distance", new taskParam(){paramUrl="inputs/Distance"});                  GPJob job = new GPJob             {                 jobId = "jb10b153b77e74a839e31af53e5d4d0d1",                 jobStatus = "esriJobSucceeded",                 results = resDict,                 inputs = inputDict,                 messages = new List<message>()                 {                     new message(){type="esriJobMessageTypeInformative",description = "Submitted."},                     new message(){type="esriJobMessageTypeInformative",description = "Executing..."},                     new message(){type="esriJobMessageTypeInformative",description = "Succeeded."}                 }                  };             string json = JsonConvert.SerializeObject(job,Formatting.Indented);             return json;             /*             DataContractJsonSerializer ser = new DataContractJsonSerializer(job.GetType());             System.IO.MemoryStream ms = new System.IO.MemoryStream();             ser.WriteObject(ms, job);             string json = Encoding.Default.GetString(ms.ToArray());             return json;                 */         }              [DataMember(Order=0)]         public string jobId { get; set; }         [DataMember(Order = 1)]         public string jobStatus { get; set; }              [DataMember(Order = 2)]         public Dictionary<string, taskParam> results { get; set; }              [DataMember(Order = 3)]         public Dictionary<string, taskParam> inputs { get; set; }         [DataMember(Order = 4)]         public List<message> messages { get; set; }     }          [DataContract]     public class taskParam     {         [DataMember]         public string paramUrl { get; set; }     }     [DataContract]     public class message     {         [DataMember(Order=0)]         public string type { get; set; }         [DataMember(Order = 1)]         public string description { get; set; }     }      Here's the output:               {       "jobId": "jb10b153b77e74a839e31af53e5d4d0d1",       "jobStatus": "esriJobSucceeded",       "results": {         "output_Buffer": {           "paramUrl": "results/output_Buffer"         },         "Final_wellpoints": {           "paramUrl": "results/Final_wellpoints"         }       },       "inputs": {         "Input_Features": {           "paramUrl": "inputs/Input_Features"         },         "Distance": {           "paramUrl": "inputs/Distance"         }       },       "messages": [         {           "type": "esriJobMessageTypeInformative",           "description": "Submitted."         },         {           "type": "esriJobMessageTypeInformative",           "description": "Executing..."         },         {           "type": "esriJobMessageTypeInformative",           "description": "Succeeded."         }       ]     }      **Update2** Looks like I spoke too soon, deserialization of this string fails:               string json = "{\"jobId\":\"jeabfba358b69412abeafd63e415957bf\",\"jobStatus\":\"esriJobWaiting\"," +                               "\"results\":[],\"inputs\":[],\"messages\":[]}";     GPJob j= JsonConvert.DeserializeObject(json, typeof(GPJob),JsonSerializerSettings) as GPJob;      I get this exception:               Cannot deserialize JSON array into type 'System.Collections.Generic.Dictionary`2[System.String,Anadarko.iMaps.JobLib.taskParam]'      It works ok when the results and inputs are not empty. **Update3:** This string deserializes:               {"jobId":"jeabfba358b69412abeafd63e415957bf","jobStatus":"esriJobWaiting","results":null,"inputs":null,"messages":null}      So here's the hack:               // hack: replace [] with null     string hackedJson = json.Replace("[]", "null");     GPJob j= JsonConvert.DeserializeObject(hackedJson, typeof(GPJob)) as GPJob;      Of course this results in a null instead of an empty dictionary. Any idea where that might cause an issue (if so, I guess I could handle that in the setters).


title: Code-Point Polygon Data and Google Maps

text: I'm currently working with OS Code-point polygons and google maps coordinates. I have google maps point coordinates for a list of houses which are located at the center of each full postcode (according to google maps) they correspond to. The code-point polygon OS data has polygons for each full postcode marking their boundaries, but also has a fair few (5%-ish) small polygons, called "vertical streets" - small areas (usually tall buildings) within other post code polygons which have more than one postcode assigned to them. These each have a reference number and separate look-up table which gives the (usually) 2-3 postcodes that are assigned to that building. My issue is - How can I try and assign a postcode in buildings that fall in the vertical streets, given that I don't have a full address (only street name) for these properties?. I have checked online and can't find any indication of how this can be dealt with. Google maps postcode data is not the most accurate either - I'm not sure how to deal with this. Thanks -


title: Is there any way to perform spatial filter on multiple layers together?

text: I have a scenario where a featureclass[Child Version] is queried for some overlapping features on [SDE.DEFAULT] Version using ISpatialFilter, After getting few result, these features are looped and again queried the same way on [Child version]. (exactly opposite). This takes lot of time. I still hope that the compexity can be reduced somehow. Is there any ways that we can run spatial filter in some kind of combination, so as to increase performance.   I have tried Spatial Cache, but because extent is not known it takes lot of time to fill the cache. Any suggestions?


title: How to scrape data from a google map (in flash)?

text: I've contacted the author of this NYC homicide map from the NY Times several times to see if he could send me a spreadsheet of information the data is based off of, but I have not had any response (for a very long time). Is there anyway I could scrap the data from the google map directly? I downloaded the swf file and converted to xml, but even within there I can not tell where the information is I want. Here is the zipped up file of the converted xml. I'm hopeing someone with more knowledge of online mapping can quickly tell me if this possible or not (or where I should be looking in the xml file).


title: MapServer WMS not configured properly?

text: This is related to another question I asked just recently: Use OpenLayers WMSGetFeatureInfo on Mapserver Layer? unicoletti indicated that he/she was able to use OpenLayer's WMS Layer declaration to generate and query MapServer WMS. However, when I attempt the same, I get a error response: **mapserv(): Web application error. Query mode needs a point, imgxy and mapxy are not set.** This error makes sense in that when I use OL's MapServer Layer declaration, the url string includes these options with values filled in. But using the WMS Layer declaration with standard WMS GetMap options, i.e. BBOX, WIDTH, HEIGHT, generates this error. The image is not correct either which I describe next. Another symptom is that the entire extent of the data defined in the mapfile is displayed rather than what's indicated in the BBOX when just doing a standard WMS GetMap request via a URL rather than using OpenLayers. So for example, here's the first part of my mapfile:               MAP        NAME "testgrib"        IMAGETYPE      PNG        EXTENT         -14000000 3000000 -7000000 7000000      And here's a GetMap request using standard WMS options: http://localhost:8080/cgi- bin/mapserv?REQUEST=GetMap&SERVICE=WMS&VERSION=1.1.1&FORMAT=image%2Fpng&MODE=map&MAP=wxmap%2Fglmpvis.map&LAYERS=mosaic&SRS=EPSG%3A3857&TRANSPARENT=TRUE&STYLES=&BBOX=-13807095.590531,3308267.5857033,-12241665.251469,4482340.3400001&WIDTH=640&HEIGHT=480 As you can see, the BBOX extent values are different from those in the mapfile. But, MapServer returns an image with the extent values in the mapfile no matter what BBOX values are passed. I've changed the mapfile extent to match the BBOX values and sure enough, an image of only that area is returned. Hence this leads me to believe we've missed something during our MapServer installation, some setting or configuration. I've checked to see that WMS was supported and received SUPPORTS=WMS_SERVER when making a GetCapabilities request as well as a -v on the MapServer executable. I considered attaching the options to the query function that I can see are used to generate images using OL's MapServer Layer declaration to get around this issue. But, as we eventually want our customers to access our data through the WMS, solving that would seem to be the main goal here. I've searched the web and MS docs for a possible setting/configuration, but nothing really stood out as being critical. I do have "wms_enable_request" "*" in the map file and, as indicated, I am receiving images from MapServer, albeit for the max extent of the data vs the viewing area. And I'm receiving a response back to the query, so this doesn't look like a cross domain/proxy issue. Also, this occurs with displaying an image for raster data as well as an image of point data. The only thing that works is to use OL's MapServer declaration to generate the images properly for the viewing area. But, that's not what we want for a long term solution with our customers accessing the data via WMS. So, open to suggestions as to what we should either configure differently or what we might've missed to enable standard WMS in our MapServer. Thanks!


title: Is anyone using JSONP in Web mapping applications to get around the browser same origin policy?

text: Javascript based Web mapping apps that use ajax requests to get data from a server are limited by the Same Origin Policy to only request the data from a server in the same domain (subdomain and port). One method commonly used to get around this is to have a script that acts as a proxy located within the same domain. The proxy sends the request to the server outside of the domain and then forwards the results back to the requesting browser. I have heard of some non-spatial Web apps using JSONP "Javascript Object Notation with Padding" to get around the same origin policy. I haven't seen it in use in any Web mapping apps. Can anyone identify an example? What are the pros and cons?


title: Is it possible to integrate CA SiteMinder with ArcGIS Server 10?

text: I am researching whether or not it is possible and/or desirable to integrate CA SiteMinder with ArcGIS Server 10 in order to restrict access to the various web service endpoints based for authenticated users. I have googled for whitepapers, solutions, and documentation but every search is coming up blank.   * Has anyone had experience with a CA SiteMinder / ArcGIS Server integration?   * If so, is there any documentation on the integration?


title: Unwanted Italicization of Symbol Text (ArcMap 9.3.1)?

text: So I am running in to a problem that has got me stumped. I'm working on a road systems map, and there are a few symbols with text in the middle. The problem is that the text within the badge, triangle, etc. is randomly becoming italicized. I have gone into the symbol property editor, made sure that the text layer has no italics, rotation, anything. However the text is still being represented in italics when on the map? The symbol itself is not rotated or shifted in any way, just the text inside. I might be missing something simple, but I cannot figure it out for the life of me. I am using ArcMap 9.3.1, and if it matters, the MPS Atlas Extension.


title: How to delete a non-empty store in geoserver

text: I have a workspace in with many stores, while i'm trying to delete the store an error shows up saying that it is unable to delete non-empty store. I tried to delete it after adding shape file the store, but only the shapefile is deleted.


title: Geoserver raster (in geotiff format) becomes black when grouped by other layers in vector format?

text: Good day, I am a newbie in field of GIS systems. Currently I have an issue when grouping layers, raster and vector. I would like to use satellite image in geotiff format as a base layer. It works well separate but when I group it by other layers it becomes black. Could anyone explain, why? I dont know if it helps, here is the details of geotiff in geoserver: > Coverage Band Details Band Data type Null Values Unit RED_BAND Unsigned 16 > bits 256.0 GREEN_BAND Unsigned 16 bits 256.0   >  BLUE_BAND Unsigned 16 bits 256.0


title: What expression returns data defined color from a color field eg: (30 50 40 10) using color_cmyk

text: Hi I have a geology plan with a base color field. I need to write an expression in the data defined color option using "color_cmyk" command. I can not seem to get the expression written correctly. The data is in the format "30 80 20 10" Thanks Nick


title: How to split line at a specified distance and getting coordinates of start and end points of each line segment？

text: I tried **v.split** in GRASS： ![enter image description here](http://i.stack.imgur.com/vbUQd.png) and my notebook ran a lot of time but showed this error： **bad allocation** , then I click "OK", and got nothing. Are there other ways to split lines？and I hope to get 2 new columns in the same shapefile about coordinates of start and end points of each line segment.


title: How to make mapable geo data from a solid jpg/png

text: I have a bunch of hardiness zones, as you can see here: http://mangoly.herokuapp.com/workroom/googlemaps/zones/ They are solid colored PNG images. I know the geo coordinates of the bottom left hand corner of all the images (they are all in the same corner). I am able to add these to the Google Maps API as groundoverlays with no problem. However I would like to make these layers clickable. So my guess is I need to convert the images into some kind of vector/geo data format. Because since each image has white space, they are all stacked on top of each other. So when I click on a layer, it will think I am choosing top most image. I need to some how convert the image into a format that make google maps understand that I only care about the actual solid color and not the white space (polygons maybe?). Any ideas?


title: Trying to match delimited text layer to my shapefile

text: I'm pretty new to QGIS and having a bit of trouble getting a text layer to match my shapefile. I'm trying to plot Canadian cities on a shapefile I have of Canada provinces. The provinces shapefile is from: http://www5.statcan.gc.ca/access_acces/alternative_alternatif.action?l=eng&dispext=zip&teng=gpr_000b11a_e.zip&k=%20%20%20%2040968&loc=http://www12.statcan.gc.ca/census- recensement/2011/geo/bound-limit/files-fichiers/gpr_000b11a_e.zip and the city data I got from: geocoder.ca under there free data downloads. So I tried reading a few articles on this and enabled on the fly CRS transformation. Looking at the provinces shapefile for the metadata it says: In layer spatial reference system units : xMin,yMin -141.018,41.6814 : xMax,yMax -52.6194,83.1355 Layer Spatial Reference System: +proj=longlat +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +no_defs So when I plotted the city file I choose to project it using WGS 84. In layer spatial reference system units : xMin,yMin 0,0 : xMax,yMax 140.87,76.4935 Layer Spatial Reference System: +proj=longlat +datum=WGS84 +no_defs and unfortunately my city data doesn't look plotted right. It looked mirrored and there is a point that is way off (guessing the data is wrong). Not sure why it would be mirrored though. Suggestions? Since I don't have the rep yet to post images I uploaded an image of what my output looks to skydrive: http://sdrv.ms/YV2Whl Cheers


title: Using multiple LayerIDs in a Legend

text: This should be simple, but I cannot find any documentation on how to do it. I have a legend that is bound to my map that contains several layers. I do not want all of the layers to appear on my Legend so I am using LayerIDs. The problem is that I do not know how to use multiple ones in the LayerIDs. Below is my code: XAML(Legend):               <esri:Legend Map="{Binding ElementName=MyMap}"                              LayerIDs="CalvertCity, WaterValves"                              LayerItemsMode="Tree"                              ShowOnlyVisibleLayers="False">                         <esri:Legend.MapLayerTemplate>                             <DataTemplate>                                 <StackPanel Orientation="Horizontal">                                     <CheckBox Content="{Binding Label}"                                         IsChecked="{Binding IsEnabled, Mode=TwoWay}"                                         IsEnabled="{Binding IsInScaleRange}" >                                     </CheckBox>                                     <Slider Maximum="1" Value="{Binding Layer.Opacity, Mode=TwoWay}" Width="50" />                                 </StackPanel>                             </DataTemplate>                         </esri:Legend.MapLayerTemplate>                         <esri:Legend.LayerTemplate>                             <DataTemplate>                                 <CheckBox Content="{Binding Label}"                                 IsChecked="{Binding IsEnabled, Mode=TwoWay}"                                 IsEnabled="{Binding IsInScaleRange}" >                                 </CheckBox>                             </DataTemplate>                         </esri:Legend.LayerTemplate>                     </esri:Legend>      XAML(ESRI Maps):               <esri:Map x:Name="MyMap" Background="#FFE3E3E3" WrapAround="True"  MouseRightButtonDown="MyMap_MouseRightButtonDown" Extent="-9834972.92753924,4441899.425293319,-9833977.88119163, 4442762.485358352">                     <esri:Map.Layers>                         <esri:ArcGISTiledMapServiceLayer ID="MyLayer" Url="http://services.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer" />                         <esri:GraphicsLayer ID="CandidateGraphicsLayer">                             <esri:GraphicsLayer.MapTip>                                 <Grid>                                     <Rectangle Stroke="Gray"  RadiusX="10" RadiusY="10" Fill="#77FF0000" Margin="0,0,0,5" >                                         <Rectangle.Effect>                                             <DropShadowEffect/>                                         </Rectangle.Effect>                                     </Rectangle>                                     <Rectangle Fill="#DDFFFFFF" Stroke="DarkGray" RadiusX="5" RadiusY="5" Margin="10,10,10,15" />                                     <StackPanel Orientation="Vertical" HorizontalAlignment="Center" Margin="30,20,30,30">                                         <TextBlock Text="{Binding [Address]}" HorizontalAlignment="Left" Foreground="Black" />                                     </StackPanel>                                 </Grid>                             </esri:GraphicsLayer.MapTip>                         </esri:GraphicsLayer>                         <esri:ArcGISDynamicMapServiceLayer ID="CalvertCity" Url="http://localhost:6080/arcgis/rest/services/CalvertCity_Test_2/MapServer"/>                         <esri:FeatureLayer ID="WaterValves" Url="http://localhost:6080/arcgis/rest/services/CalvertCity_Test_2/MapServer/0" Renderer="{StaticResource MySimpleRenderer}"                                                       OutFields="WVAL_ID,SIZE,WVAL_USE,WVAL_TYPE,OPENS,TURNS" MouseLeftButtonUp="FeatureLayer_MouseLeftButtonUp"/>                     </esri:Map.Layers>                                 </esri:Map>      If you need any more information, just let me know. I will appreciate any help. Thank you in advance.


title: How to delete rings automatically using QGIS/Grass

text: I'm cleaning up a (very) messy shapefile in QGIS/Grass. I now have several large polygons with numerous (100s) of small rings. Is there a quick and easy way to remove/ fill all these small rings. Selecting and deleting each one manually is not feasible.


title: Is it possible to permanently save the Edit Widget for a Postgis layer?

text: Is there a way of configure the editing type for a POSTGIS layer through Edit Plugin and do not lose the changes after remove the layer or exit the program?


title: Easy Search of Road or Address in QGIS

text: Is there an easy way to search for a road name in a polyline layer in QGIS? I'm looking for something like what ArcMap 10.1 has, just a basic search box where I can type in a name and have it give me a list of matches, then have the option to zoom to it, or highlight it, or something


title: How do I reach a layer inside a GroupLayer?

text: I am having problems in ArcObjects. I am struggling to get to the individual (feature-)layers inside a GroupLayer. Does anyone know how to do that ?


title: Create a new vector layer that fills space between other polygons in QGIS

text: I inherited a clipped (500m radius) land cover map containing vector layers of different land cover categories such as buildings, grass, roads, etc from a colleague. ![Land Use Map](http://i.stack.imgur.com/ZaLUD.jpg) The empty (white) space in the map is actually paved surfaces in reality but is it possible to fill that "empty space" with another vector layer that creates a new landcover polygon class automatically? I'm under the impression that it is possible to create another layer with the symmetrical difference tool but I think I need two input vector files - the one with all the current layers and another layer. Am I supposed to create another vector layer that covers the exact same area (i.e. 500 m radius about a certain point?) and how do I go about doing that? \--------- EDIT ---------- Ok just to summarize, this is what I did to get the difference layer:   1. Add a vector layer with center-point to existing maps.   2. Create a buffer layer using this centerpoint. Vector > Geoprocessing tools > Buffer. Use the new vector layer with center-point as input layer.    3. With the new buffer layer, you can now use the Difference Tool (Vector > Geoprocessing > Difference). Use the buffer layer as your input layer, and the other layer being your main map.    4. You will have to edit the attributes of the new layer so it corresponds to the same type of attributes in your base map.   5. Merge the layers using the mmqgis tool (mmqgis > transfer > merge layers), and create a new shapefile with all the layers merged together.


title: How to use GDAL to burn a selected value from one raster into another?

text: I think this should be a relatively simple process, but I'm not quite sure how to approach it. I have a shapefile containing 1000's of polygons. I also have 19 rasters. All of the polygons fall withing the extent of the rasters. The rasters are equal in extent and resolution. Previously I have used GDAL to extract data for each of these polygons from each of the 19 rasters. This worked very well. However, I now wish to extract the pixels into xyz format or similar for all the areas **not** covered by these polygons. Currently I am using `gdal_rasterize` to `burn` a chosen nodata value into each of the rasters indicating where the polygons occur. I would then need to translate the raster removing all of the nodata values. This is a very heavy process and is taking an awefully long time for 1 raster let alone 19. My question is - is there a way to `gdal_rasterize` all of my polygons onto just 1 of the rasters, then use the nodata values from this raster as a mask, and apply these nodata values to the same locations on the other 18 rasters. I'm sure this would be quicker than my current method, but I don't know how to do this in GDAL. Any help will be very much appreciated!


title: How to visualize the water equivalent between two raster layers in arcmap 10.1?

text: I do have two raster layers (interpolated with Empirical Bayesian Kriging), that show the elevation of a glacier in April, and September. I want to show the change in elevation between the layers in water equivalent. Any suggestion on how I can do that?


title: Identifying clusters of points within a given distance

text: So I have a database of lats/longs representing incidents occurring on or near roadways. My goal is to identify the street block areas with the highest number of incidents in the database. I don't have any other map layers I'm using for this. What I'd like to do is, for each point, see how many other points are a given defined distance away (a quarter of a mile, for instance). The result would then show where the clusters are. Any recommendations on the best way to do this? I have copies of both QGIS (with the spatial query plugin) and ArcMap, although I'm much more familiar with QGIS.


title: How to find closest point on multi line String from a point in PostGIS

text: I want to find out the nearest point on my road network which is a multiline string from a point with lat lon value. I have tried the following codes but it giving a error message that my road network is not a line.               SELECT ST_AsText(ST_Line_Interpolate_Point(foo.the_line, ST_Line_Locate_Point(foo.the_line, ST_GeomFromText('POINT(90.402 23.752)'))))     FROM (SELECT ST_GeomFromText(geom) As the_line from roads) As foo;      Where `geom` is my geometry column with multiLineString   `roads` is the table name   Long `90.402` and Lat `23.752`   Projected Coordinate system UTM_Zone_46N


title: Can this WPS be used in QGIS?

text: is it possible to integrate the following WPS-LINK in QGIS? http://gdz.bkg.bund.de/wps_koordtrans?REQUEST=GetCapabilities&SERVICE=WPS


title: How to name functions in OpenLayers?

text: I don't know well enough OpenLayers syntax, so I need to ask here. I need to name every function in OpenLayers (for creating map - tiles, for drawing markers on the map etc.). So, I have this simple code:               map = new OpenLayers.Map('map');     var wms = new OpenLayers.Layer.WMS( "OpenLayers WMS",     "http://vmap0.tiles.osgeo.org/wms/vmap0?", {layers: 'basic'});      I would need to put that inside the function and name that function, something like this:               function createMap() {         map = new OpenLayers.Map('map');         var wms = new OpenLayers.Layer.WMS( "OpenLayers WMS",         "http://vmap0.tiles.osgeo.org/wms/vmap0?", {layers: 'basic'});                    }          createMap();      But that code doesn't work. So, can you tell me if this is possible and why my code doesn't work? Thanks.


title: Is it possible to rotate labels in Print Composer?

text: Is it possible to rotate a label in Print Composer (QGIS 1.8)? Thx.


title: Is there is a index that calculates for Sand in ENVI?

text: On a remote sensing project, I'm trying to identify areas that are sand using band math. I know there is an index for vegetation (NDVI). I was wondering if there is an index that calculates for sand. I tried googling it and I couldn't really find anything substantial. Thank you.


title: File Validation using ArcPy

text: I have a script which scans a directory and outputs basic raster data information such as the file name, format, number of bands, and etc. I need a way to make it so if the directory does not contain raster data (i.e., anything other than raster data), a message is displayed stating that the directory doesn't have the correct data type. I know ArcPy has a `Describe()` function that I could use to determine the type of data in a folder, but am not sure how to implement it. This is what I have so far:               rasterList = arcpy.ListRasters("*", "ALL")     filesType = arcpy.DataType('RasterDataset') # Can use `DatasetType` as well.                                                  # I've tested this function to describe                                                 # raster data and ArcPy prints out                                                 # 'RasterDataset', that is why I have it                                                  # there in the brackets.     for name in rasterList:         if rasterList == filesType:             print ("\nFilename:"), name         else:             print ("This directory does not contain any raster data.")      Any suggestions?


title: 'File is too big' error when uploading QGIS Plugin

text: I have written a simple QGIS plugin which zooms the map canvas to a user specified UK postcode. I would like to share the plugin but the csv file, which contains stripped down OS Code-Point open data is just over 40mb. This is too large to upload to the qgis repository and impractical to download and share. I have also tried storing the postcodes in a JSON file and pickle file, both of which are also to large a file size. Is there any alternative way of storing the postcodes which may result in a smaller file size? Thanks


title: Migrating code to 10.1 and recycled search

text: I am migrating 9.3 Desktop extension to 10.1 I am getting this bug when I am trying to update a row in the geodatabase "Cannot call Store on a recycled row while editing" And seems the reason is, I am doing recyclable search like this eatCursor = featLayer.Search(spatialFilter, true); So, I changed it to false recyclable. And it worked. But my question is why 10.1 is not able to update recycled row, where 9.3 was able to? what the reason of this change? Is there any reason 10.1 disabled this feature?


title: How to do a looping to Interpolate a huge Rainfall Time-series in GRASS GIS or R

text: In GRASS GIS I have access to a table in PostgreSQL Database with this structure: Table Name: rainfall_time_series Columns: id | year | month | day | value_mm | x | y | gauge_code Setting g.region to the Watershed limits, I need create for each day one interpolated raster, representing the rainfall. This times-series is very huge, has 17 years, my question is: is it possible to do this operation with a looping? it is probably easier in R, but I'm newbie in this software. Thanks in advance


title: How to show just some labels based on feature intersection of layers?

text: I'm trying to use "expression based" labeling in QGIS, but I'm stuck with a particular case. I managed to use it for a project with the help of other topics and FAQ but I'm not skilled enough to take care of this problem alone. I have 5 layers, one "project" layer with around 10 objects (10 different marinas around rivers in France) and 4 layers for particular water area (different scales) with also polygons. I want to only show labels of those 4 last layers of the polygons that intersect any of my 10 objects of layer 1. I hope to be clear enough. Any guess of the query I need to use? _Note: I already used the vector selection option and created a new layer with the selection. I'm looking for another option, label based because this is too long. I know there is also a "show/hide" option, but it is not very simple to use in my opinion._ ![enter image description here](http://i.stack.imgur.com/3EFlO.jpg) Thank you for your time!


title: QGIS crashes when using dxf to shp converter

text: this is my first question on Stackexchange. I use a GPS to capture data in .dxf, manipulate it in a CAD program that I know how to use, then the dxf2shp converter (I think Python is mentioned) to let me generate .shp files for the end user. Every time I've done this recently I've got a meesage from QGIS that it has stopped working and will close. I've tried re-installing QGIS (1.8.0), but still have the same message. Am I doing something wrong or is this a known problem. Please note I'm not really computer literate, so if you could avoid assuming I know what a technical term means, i'd be grateful.


title: Field Calculator Does not Work in Arc10

text: I'm having a problem, described below, with Field Calculator in ArcGIS 10. I wonder if there is a known fix at this point. I'm updating to 10.1 soon, but this work cannot wait until then. **Details:** I'm trying to divide two fields and then multiply the answer by 100 into a calculated field. I have created a field to receive the answers, "Grad_Per". I highlight the [Grad_Per] field, right click, select Field Calculator, and enter [elev_range]/[SLength]*100. Answers then populate the Grad_Per field down to record number 564 but no further. The answers up to that point are correct; I've checked them in MS Excel and by hand. It's just that I've got about 9500 more records to calculate and they have values of 0 after the calculation is complete. Thanks so much for your help, Vin


title: Issues for 3D length using GPS road tracking for undulated topography

text: We are tracking road network for undulated topography using Trimple Juno GPS. Is it any problems for 3D length using UTM projection? Thanks to Advance


title: How do I delete multiple data points from a .csv file through QGIS?

text: I have a .csv file with points in North America, but I only want to focus on a section through North America that excludes places such as Alaska and Greenland. I've selected the points below, but the Delete Selected button is greyed out. Is there any way to delete these points through QGIS or will I have to edit my .csv file and import it again? http://imgur.com/cjqOdFG


title: How to generate routes that follow underlying road network?

text: I have a set of points (~2300 bus stops) and I have a network of roads that these bus stops lie on. When I have a bus route that goes from A to B, how do I generate the route such that it follows the underlying road and not just the direct crowfly distance. ![enter image description here](http://i.stack.imgur.com/x4FwJ.png) The road data is from OSM and the stop data is user generated. At present the data is hosted on :routemaster.hydlab.in


title: How to add multiple OpenLayers Vector Layer markers with data parsed from JSON file

text: I would like to parse a local JSON file and add markers to a map based on the latitude and longitude parsed from the JSON file. Here is a sample from the local JSON file:               "LOCATION_A": {                   "latitude": 10.2070,                   "longitude": 25.3215,                   "altitude": 1.2,                   "status": "CLOSED",                   "serverAddress": "192.168.0.1"      }      I am unsure how to proceed beyond the following block of code:               var markers = new OpenLayers.Layer.Vector("Sites", {       projection: map.displayProjection,       strategies: [new OpenLayers.Strategy.Fixed()],       protocol: new OpenLayers.Protocol.HTTP({         url: "filename.json",         format: new OpenLayers.Format.JSON()         }),       });      I am unsure how to use the parser to do this. I would like to use jQuery for this exercise now. Any assistance would be appreciated. Thank you in advance.


title: How to delete points from CSV point layer?

text: I have a point layer produced from importing a Delimited Text Layer and I would like to delete some points. It seems that it is not possible to toggle editing in a point layer. How can I delete the selected points? (QGIS version 1.8.0-Lisboa)


title: Global Mapper Opening Error

text: I am having a problem when I am opening Global Mapper workspace that I had already created and saved it. Evey time I try to open it an error appear and nothing opens on the screen. Is there a way to open this .gmw file. Thank you.


title: ogr SimplifyPreserveTopology does not keep the Topology

text: I use `ogr2ogr` to simplify a layer of polyline, from the document: > -simplify tolerance: (starting with GDAL 1.9.0) distance tolerance for > simplification. This method will preserve topology, in particular for > polygon geometries. However, it seems that the topology are not kept. http://pbrd.co/SPt39N What's the problem? And the gdal version I use is `1.10.0`.


title: Georeferencing a column doesn't work in my local instance of CartoDB

text: I'm installing a instance of cartodb in AWS/Ubuntu 12.04, and I'm having a hard time trying to georeference a column with textual addresses. I know that the data I'm trying to georeferencing is ok since I uploaded it in an online cartodb account and it's working there, so it must be a local problem. I've tried with different cartodb versions, and I get different errors: v2.9.3 to v2.10.1 -> I upload the file, table is created, then I select 'georeference' in the proper column and seems that the job is processed by resque (progress bar appears in UI), two columns are added (the geom (all values as null) and cartodb_georef_status (all values as false)), but they are not updated later with proper values. Checking the log, I can see this:                   18:54:42 resque.1  | ** [18:54:42 2014-03-24] 7121: got: (Job{geocodings} | Resque::GeocoderJobs | [{"job_id"=>1}])         18:54:42 resque.1  | ** [18:54:42 2014-03-24] 7199: resque-1.23.0: Processing geocodings since 1395687282         18:54:42 resque.1  | ** [18:54:42 2014-03-24] 7121: resque-1.23.0: Forked 7199 at 1395687282         18:54:44 web.1     | (0.009665s) SELECT * FROM "geocodings" WHERE "id" = 1         18:54:44 web.1     | (0.003518s) SELECT * FROM "user_tables" WHERE ("user_tables"."id" = 1) LIMIT 1         18:54:44 web.1     | (0.004384s) SELECT * FROM "users" WHERE ("users"."id" = 1) LIMIT 1         18:54:44 web.1     | [pool] Creating a new connection for cartodb_dev_user_1_db:postgres (1)         18:54:44 web.1     | (0.005848s) BEGIN         18:54:44 web.1     | (0.000490s) UPDATE "geocodings" SET "updated_at" = '2014-03-24 18:54:43.110957+0000', "state" = 'started' WHERE ("id" = 1)         18:54:44 web.1     | (0.002168s) COMMIT         18:54:44 web.1     | (0.000164s) SET standard_conforming_strings = ON         18:54:44 web.1     | (0.000165s) SET client_min_messages = 'WARNING'         18:54:44 web.1     | (0.000161s) SET DateStyle = 'ISO'         18:54:44 web.1     | (0.050511s)          18:54:44 web.1     |         ALTER TABLE map_03_2          18:54:44 web.1     |         ADD COLUMN cartodb_georef_status BOOLEAN DEFAULT FALSE         18:54:44 web.1     |                18:54:44 web.1     | (0.001902s)          18:54:44 web.1     |             SELECT md5(trim(both from regexp_replace(regexp_replace(concat(address), E'[\n\r]+', ' ', 'g'), E'"', '', 'g'))) AS searchtext         18:54:44 web.1     |             FROM map_03_2         18:54:44 web.1     |             WHERE cartodb_georef_status IS false OR cartodb_georef_status IS NULL         18:54:44 web.1     |             GROUP BY searchtext         18:54:44 web.1     |             LIMIT 10 OFFSET 0         18:54:44 web.1     |                  18:54:44 web.1     | (0.000194s) DROP TABLE IF EXISTS geocoding_cache_1395687283         18:54:44 web.1     | [Rollbar] Scheduling payload         18:54:44 web.1     | [Rollbar] Sending payload         18:54:44 web.1     | [Rollbar] Got unexpected status code from Rollbar api: 401         18:54:44 web.1     | [Rollbar] Response: {         18:54:44 web.1     |   "err": 1,         18:54:44 web.1     |   "message": "unauthorized"         18:54:44 web.1     | }         18:54:44 web.1     | [Rollbar] Details: https://rollbar.com/instance/uuid?uuid=3a66799c-34fa-4733-bf78-0a251fd5bf24 (only available if report was successful)         18:54:44 web.1     | (0.000767s) COPY (         18:54:44 web.1     |           SELECT trim(both from regexp_replace(regexp_replace(concat(address), E'[\n\r]+', ' ', 'g'), E'"', '', 'g')) as recId, trim(both from regexp_replace(regexp_replace(concat(address), E'[\n\r]+', ' ', 'g'), E'"', '', 'g')) as searchText          18:54:44 web.1     |           FROM map_03_2         18:54:44 web.1     |           WHERE cartodb_georef_status IS FALSE OR cartodb_georef_status IS NULL         18:54:44 web.1     |           GROUP BY recId         18:54:44 web.1     |           LIMIT 1000000         18:54:44 web.1     |       ) TO STDOUT (FORMAT csv, HEADER)         18:54:44 web.1     | ETHON: Libcurl initialized         18:54:44 web.1     | ETHON: performed EASY url= response_code=0 return_code=url_malformat total_time=0.0         18:54:44 web.1     | ETHON: performed EASY url= response_code=0 return_code=url_malformat total_time=0.0         18:54:44 web.1     |          18:54:44 web.1     |          18:54:44 web.1     | Started GET "/api/v1/geocodings/1" for (IP) at 2014-03-24 18:54:44 +0000         18:54:44 web.1     | Processing by Api::Json::GeocodingsController#show as JSON         18:54:44 web.1     |   Parameters: {"id"=>"1"}         18:54:44 web.1     | (0.017074s) SELECT * FROM "users" WHERE ("username" = 'ec2-54-185-162-75') LIMIT 1         18:54:44 web.1     | (0.014046s) SELECT * FROM "geocodings" WHERE (("geocodings"."user_id" = 1) AND ("id" = '1')) ORDER BY "created_at" DESC LIMIT 1         18:54:44 web.1     | (0.004323s) SELECT * FROM "user_tables" WHERE ("user_tables"."id" = 1) LIMIT 1         18:54:44 web.1     | (0.000345s) SELECT * FROM "users" WHERE ("users"."id" = 1) LIMIT 1         18:54:44 web.1     | [pool] Found a connection for cartodb_dev_user_1_db:postgres (2)         18:54:44 web.1     | Completed 200 OK in 286ms (Views: 0.8ms | Models: 0.0ms)         18:54:44 resque.1  | 2/3         18:54:45 web.1     | (0.008280s) BEGIN         18:54:45 web.1     | (0.003983s) UPDATE "geocodings" SET "updated_at" = '2014-03-24 18:54:44.365433+0000', "remote_id" = 'c3489b58b38511e38e7722000aed033e' WHERE ("id" = 1)         18:54:45 web.1     | (0.002151s) COMMIT         18:54:45 web.1     | (0.000305s) BEGIN         18:54:45 web.1     | (0.000401s) UPDATE "geocodings" SET "total_rows" = 3, "processed_rows" = 2, "updated_at" = '2014-03-24 18:54:44.530939+0000', "state" = 'completed' WHERE ("id" = 1)         18:54:45 web.1     | (0.001879s) COMMIT         18:54:45 web.1     |          18:54:45 web.1     |          18:54:45 web.1     | Started GET "/api/v1/geocodings/1" for (IP) at 2014-03-24 18:54:45 +0000         18:54:45 web.1     | Processing by Api::Json::GeocodingsController#show as JSON         18:54:45 web.1     |   Parameters: {"id"=>"1"}         18:54:45 web.1     | (0.002891s) SELECT * FROM "users" WHERE ("username" = 'ec2-54-185-162-75') LIMIT 1         18:54:45 web.1     | (0.000345s) SELECT * FROM "geocodings" WHERE (("geocodings"."user_id" = 1) AND ("id" = '1')) ORDER BY "created_at" DESC LIMIT 1         18:54:45 web.1     | (0.003857s) SELECT * FROM "user_tables" WHERE ("user_tables"."id" = 1) LIMIT 1         18:54:45 web.1     | (0.000350s) SELECT * FROM "users" WHERE ("users"."id" = 1) LIMIT 1         18:54:45 web.1     | [pool] Found a connection for cartodb_dev_user_1_db:postgres (2)         18:54:45 web.1     | Completed 200 OK in 148ms (Views: 0.9ms | Models: 0.0ms)         18:54:47 web.1     | (0.000256s) BEGIN         18:54:47 web.1     | (0.000428s) UPDATE "geocodings" SET "updated_at" = '2014-03-24 18:54:46.543169+0000', "cache_hits" = NULL WHERE ("id" = 1)         18:54:47 web.1     | (0.001688s) COMMIT         18:54:47 web.1     |          18:54:47 web.1     |          18:54:47 web.1     | Started GET "/api/v1/geocodings/1" for (IP) at 2014-03-24 18:54:47 +0000         18:54:47 web.1     | Processing by Api::Json::GeocodingsController#show as JSON         18:54:47 web.1     |   Parameters: {"id"=>"1"}         18:54:47 web.1     | (0.016115s) SELECT * FROM "users" WHERE ("username" = 'ec2-54-185-162-75') LIMIT 1         18:54:47 web.1     | (0.008763s) SELECT * FROM "geocodings" WHERE (("geocodings"."user_id" = 1) AND ("id" = '1')) ORDER BY "created_at" DESC LIMIT 1         18:54:47 web.1     | (0.000302s) SELECT * FROM "user_tables" WHERE ("user_tables"."id" = 1) LIMIT 1         18:54:47 web.1     | (0.000353s) SELECT * FROM "users" WHERE ("users"."id" = 1) LIMIT 1         18:54:47 web.1     | [pool] Found a connection for cartodb_dev_user_1_db:postgres (2)         18:54:47 web.1     | PG::Error: ERROR:  relation "cdb.geo_c3489b58b38511e38e7722000aed033e" does not exist         18:54:47 web.1     | LINE 1: SELECT NULL FROM "cdb.geo_c3489b58b38511e38e7722000aed033e" ...         18:54:47 web.1     |                          ^: SELECT NULL FROM "cdb.geo_c3489b58b38511e38e7722000aed033e" LIMIT 1         18:54:47 web.1     | (0.013806s)          18:54:47 web.1     |         CREATE TABLE cdb.geo_c3489b58b38511e38e7722000aed033e (         18:54:47 web.1     |           recId text,          18:54:47 web.1     |           SeqNumber int,          18:54:47 web.1     |           seqLength int,          18:54:47 web.1     |           displayLatitude float,          18:54:47 web.1     |           displayLongitude float         18:54:47 web.1     |         );         18:54:47 web.1     | PG::Error: ERROR:  relation "cdb.geo_c3489b58b38511e38e7722000aed033e" does not exist         18:54:47 web.1     | LINE 1: SELECT NULL FROM "cdb.geo_c3489b58b38511e38e7722000aed033e" ...         18:54:47 web.1     |                          ^: SELECT NULL FROM "cdb.geo_c3489b58b38511e38e7722000aed033e" LIMIT 1         18:54:47 web.1     | (0.000213s) COPY cdb.geo_c3489b58b38511e38e7722000aed033e FROM STDIN (FORMAT csv)         18:54:47 web.1     | PG::Error: ERROR:  relation "cdb.geo_c3489b58b38511e38e7722000aed033e" does not exist         18:54:47 web.1     | LINE 1: SELECT NULL FROM "cdb.geo_c3489b58b38511e38e7722000aed033e" ...         18:54:47 web.1     |                          ^: SELECT NULL FROM "cdb.geo_c3489b58b38511e38e7722000aed033e" LIMIT 1         18:54:47 web.1     | (0.102498s)          18:54:47 web.1     |         UPDATE map_03_2 AS dest         18:54:47 web.1     |         SET the_geom = ST_GeomFromText(         18:54:47 web.1     |             'POINT(' || orig.displayLongitude || ' ' ||         18:54:47 web.1     |               orig.displayLatitude || ')', 4326         18:54:47 web.1     |             ),         18:54:47 web.1     |             cartodb_georef_status = true         18:54:47 web.1     |         FROM cdb.geo_c3489b58b38511e38e7722000aed033e AS orig         18:54:47 web.1     |         WHERE trim(both from regexp_replace(regexp_replace(concat(address), E'[\n\r]+', ' ', 'g'), E'"', '', 'g')) = orig.recId         18:54:47 web.1     |                18:54:47 web.1     | (0.000225s) DROP TABLE IF EXISTS geocoding_cache_1395687283         18:54:47 web.1     | [Rollbar] Scheduling payload         18:54:47 web.1     | [Rollbar] Sending payload         18:54:47 web.1     | Completed 200 OK in 185ms (Views: 0.9ms | Models: 0.0ms)         18:54:47 resque.1  | ** [18:54:47 2014-03-24] 7199: done: (Job{geocodings} | Resque::GeocoderJobs | [{"job_id"=>1}])        * v2.11.0, master -> I upload the file, table is created, then I select 'georeference' in the proper column and seems that an error happens even before the job is processed by resque (progress bar doesn't appear in UI). Checking the log, I can see this:               Started POST "/api/v1/geocodings" for (IP) at 2014-03-24 18:11:22 +0000     18:11:22 web.1     | Processing by Api::Json::GeocodingsController#create as JSON     18:11:22 web.1     |   Parameters: {"formatter"=>"{address}", "table_name"=>"map_03_2"}     18:11:22 web.1     | (0.000638s) SELECT * FROM "users" WHERE ("username" = 'ec2-54-185-162-75') LIMIT 1     18:11:22 web.1     | (0.007037s) SELECT * FROM "user_tables" WHERE (("user_id" = 'b249cedd-dc33-49eb-86f1-a53f1ca17731') AND ("name" = 'map_03_2')) ORDER BY "id" DESC LIMIT 1     18:11:22 web.1     | (0.012165s) SELECT "pg_attribute"."attname" AS "name", CAST("pg_attribute"."atttypid" AS integer) AS "oid", format_type("pg_type"."oid", "pg_attribute"."atttypmod") AS "db_type", pg_get_expr("pg_attrdef"."adbin", "pg_class"."oid") AS "default", NOT "pg_attribute"."attnotnull" AS "allow_null", COALESCE(("pg_attribute"."attnum" = ANY("pg_index"."indkey")), false) AS "primary_key", "pg_namespace"."nspname" FROM "pg_class" INNER JOIN "pg_attribute" ON ("pg_attribute"."attrelid" = "pg_class"."oid") INNER JOIN "pg_type" ON ("pg_type"."oid" = "pg_attribute"."atttypid") INNER JOIN "pg_namespace" ON ("pg_namespace"."oid" = "pg_class"."relnamespace") LEFT OUTER JOIN "pg_attrdef" ON (("pg_attrdef"."adrelid" = "pg_class"."oid") AND ("pg_attrdef"."adnum" = "pg_attribute"."attnum")) LEFT OUTER JOIN "pg_index" ON (("pg_index"."indrelid" = "pg_class"."oid") AND ("pg_index"."indisprimary" IS TRUE)) WHERE (("pg_attribute"."attisdropped" IS FALSE) AND ("pg_attribute"."attnum" > 0) AND ("pg_class"."relname" = 'assets') AND ("pg_namespace"."nspname" !~* 'pg_*|information_schema')) ORDER BY "pg_attribute"."attnum"     18:11:22 web.1     | Completed 422 Unprocessable Entity in 269ms (Views: 0.3ms | Models: 0.0ms)        * v2.11.1 -> I upload the file, but I cannot view the rows in the UI (I got an error). Checking the log, I can see this:               17:30:32 web.1     | Started GET "/api/v1/tables/map_03_2/records?rows_per_page=40&page=0&mode=asc&order_by=cartodb_id&filter_column=&filter_value=&api_key=66fc98d3bddf79a6c0938f4282a2140fd9b00360" for (IP) at 2014-03-24 17:30:32 +0000     17:30:32 web.1     | [Rollbar] Reporting exception: The action 'index' could not be found for Api::Json::RecordsController     17:30:32 web.1     | (0.000598s) SELECT * FROM "users" WHERE ("username" = 'ec2-54-185-162-75') LIMIT 1     17:30:32 web.1     | [Rollbar] Exception not reported because Rollbar is disabled     17:30:32 web.1     |      17:30:32 web.1     | AbstractController::ActionNotFound (The action 'index' could not be found for Api::Json::RecordsController):


title: Writing iso19139 compliant xml metadata from a spreadsheet

text: I am looking at a way to migrate the metadata from a spreadsheet into iso19139 compliant xml metadata documents for a collection of shapefiles. I can parse the data out of the spreadsheet into strings without issue however I am getting lost on how to take those strings and use them to populate the elements. This example: http://pymotw.com/2/xml/etree/ElementTree/create.html covers most of what I want to do, however I would need to generate the tree beforehand based on the iso19139 schema. PyXB seems to have schema bundles for OpenGIS but I would like some insight on how to connect everything.


title: Line to points at bends

text: I have a pipeline (shapefile) that I am trying to convert to points where the pipeline has a bend in it. I tried to convert the vertices to points in ArcMap 10.1, however this will not work in this case as my line has a vertex every 1 metre (even when the line is straight).


title: Using Tracking Server 10.0 with ArcGIS 10.1?

text: I'm using ArcGIS 10.1 for Desktop and ArcGIS 10.1 for Server. I want to use Tracking Server 10.0 in my Web GIS. My web services are created by ArcGIS 10.1 for Server and i want to show and manage real-time data in my webgis. Is there any way to use Tracking Server 10.0 with ArcGIS 10.1 for Server?


title: PostGIS select by lat/long bounding box

text: I have a minLat, minLong, maxLat and maxLong from a box drawn on a map. Right now, my query to get all the points in the box looks like this:               SELECT *     FROM geomTable     WHERE (longitude BETWEEN minLon AND maxLon)     AND (latitude BETWEEN minLat AND maxLat)      I want to use the geometry column instead of lat/long columns to get the results. I tried this query:               SELECT *     FROM mytable     WHERE mytable.geom && ST_MakeEnvelope(minLon, minLat, maxLon, maxLat, 4326);      from this post: Select bounding box using postgis but it is returning no results. Does anyone have an example of how to select all the points within a box created by min and max lat/longs using geometry in postgis?


title: ArcGIS javascript api - layer reordering - how to add new layers?

text: I am learning ArcGIS javascript API and wanted to make a reordering panel for my webapp. I looked through the ArcGIS sample. but I'm stuck at `addLake()` function. If I'm understand the sample correctly, to add the new lake layer               1. get the existing DynamicLayerInfos from the current map     2. create a new DynamicLayerInfo for the new layer     3. push to the DynamicLayerInfos     4. set the new DynamicLayerInfos to the map      The example uses `LayerDataSource` that reference a table as the new data source for the layer. And how do you add a new layer from a public layer on a MapServer/FeatureServer. Can someone give me an example?


title: ArcGIS 10.1: "ERROR 000223: Cannot extract event table properties" when iterating through multiple files/folder

text: I am trying to convert ~270k csv-files with lat/long-data to layerfiles. The files are split into folders by year, i.e. ~30 folders with ~9k csv-files in each. I have written a script to iterate through the files and it works well for test-samples of 100+ split into a couple of different folders. However, when I try to run it on the full set of files, it will either:   1. work well for the first folder (~9k files) and then when moving on to the second folder it will give the following error for all subsequent files:   arcgisscripting.ExecuteError: ERROR 000223: Cannot extract event table properties   Failed to execute (MakeXYEventLayer).   2. on my other computer it will usually stop after an arbitrary amount of file-conversions before it finishes the first folder (typically at +-8.5k of the conversions) and then give the same error for all subsequent files The error seems to be similar to this: http://forums.arcgis.com/threads/80161-Make-XY-Event-Layer-Error ,but my error pops up after an arbitrary amount of files has been processed. I've included my code for reference below. Note that I have also tried to split the list of files per folder into chunks of 100 and defined the layercreation as a function that is being called for each chunk, but the error still persisted in the same way. I suspect that it might be some setting in ArcGIS or on my computers that causes it, but I don't know what. Does anybody have any idea where I might start to solve this?               #setting paths for the project      path="C:/project"     savepath = path + "/outputlayers"     csvpath = path+ "/csv"          #lists the directories in the csv directory     yearlist = os.listdir(csvpath)          #checks the directory of each year:     for year in yearlist:         indir = csvpath + "/" + year         hourlist = os.listdir(indir)                  #runs through each file and converts it to a .lyr file in the outputdirectory             for infile in hourlist:                 try:                     # Set the local variables                     in_Table = indir + "/" + infile                     x_coords = "longitudex"                     y_coords = "latitudey"                     z_coords = ""                     out_Layer = infile[:-4]                     saved_Layer = savepath + "/" + year + "/" + infile[:-4] +'.lyr'                          # Set the spatial reference                     spRef = arcpy.SpatialReference(4301)                          # Make the XY event layer...                     arcpy.MakeXYEventLayer_management(in_Table, x_coords, y_coords, out_Layer, spRef, z_coords)                          # Print the total rows                     print arcpy.GetCount_management(out_Layer)                          # Save to a layer file                     arcpy.SaveToLayerFile_management(out_Layer, saved_Layer)                      except:                     # If an error occurred print the message to the screen                     print arcpy.GetMessages()


title: Merge different Landsat scenes in Grass GIS

text: I am new to Grass GIS, and I would really need some help in solving this issue. I have 17 scenes from Landsat8 covering an area. I have created color infrared images of bands 5,4,3 and now I want to merge all scenes into a mosaic, to further analysis. I can't seem to manage to do this, neither with r.patch or with i.imagine.mosaic. Does anyone know another way for merging? I tried also in QGIS, I exported from Grass as tiffs , loaded them into Qgis, and used the merge command, but the images become grayscale and with some white corners overlaying. I would really appreciate some suggestions, as I am stuck :) Thanks, Clara


title: ArcGIS: Merging and joining two line feature classes while maintaining attributes from both

text: I created an event layer using the `MAKE ROUTE EVENT LAYER` tool based on anticipated scheduled events. The road layer which is the base layer is intact. What I would like to achieve is combine layer 1 (base roads) with layer 2 (scheduled events) into one Result layer as shown in the figure below while maintaining attributes from both layers. ![line up layers](http://i.stack.imgur.com/WsaG5.png) I have tried using Join, data join, spatial join, and other aspects and I still can't get it to do it right. If ArcGIS managed to merge the layers, it usually has missing data but for the most part my efforts have been unsuccessful. what is the "best of practice way" to solve my problem?


title: Problem with Performance when accessing selected Features from PostGIS layers in QGIS using Python

text: I have 2 PostGIS layer in my QGIS project. They have about the same amount of data and a similar structure except that layerX has its PK on a single column and layerY has its PK on 4 columns. When I try to access the selected features from the layers using python, the performance differs significantly. Accessing the selected features from layerX performs approximately 10 times faster than accessing the features from layerY. The code I use to obtain the selected features is the following:                   iface.activeLayer().selectedFeatures()      I have the impression that this difference in performance is based on the use of Primary Keys in the tables. Am I right? Is there maybe a faster/ better way to access featureattributes from a layer?


title: How do you programmatically access a feature's attachments?

text: I am attempting to access a feature's attachments in ArcObjects and I'm following the docs for **Working with feature attachments**. Everything is fine until attempt to access IAttachment.Data because I am expecting there to be data, but that property is null: ![enter image description here](http://i.stack.imgur.com/L2ujT.png) which means I can't proceed past: `IMemoryBlobStream memoryBlobStream = attachment.Data;` I am thinking there is a missing step (or steps) in the docs that would fill the attachment object's `Data` property with actual data? How is this done? I am able to access the attachment via the Identify window within ArcMap, so I know the data is there: ![enter image description here](http://i.stack.imgur.com/wf3NU.png) ![enter image description here](http://i.stack.imgur.com/6xpZA.png) ![enter image description here](http://i.stack.imgur.com/MC8Ix.png)


title: GeoTools and image storage / retrieval

text: Please can someone let me know how I can create a schema (using SimpleFeatureTypeBuilder()), read and write image data? Simple data types (Integer/Float etc) work fine but how do I handle image data? Thank you! Martin...


title: how to print to stdout specific columns with values corelated for other columns from a txt file

text: I'm new to python and I have a txt file with 20 columns.   I need to select 3 columns which are correlated like this:               c1  c2  c3     aaa exy hfg     aaa mmm fgh     aaa abc asd      I have to find all `c2` and `c3` that have values corresponding to `c1` value `aaa`.


title: Repair topology before pgrouting / No results of pgrouting

text: My goal is to calculate lots of distances along the streets. I followed http://workshop.pgrouting.org/ tutorial except that I used osm2po instead of osm2pgrouting to import osm file. The only problem I'm facing now are gaps in results - Dijkstra algorithm seems to fail in finding the route between two nodes of the network in 10% of cases. It simply returns nothing, empty table. I suppose the problem is with the map, that it's not "routable enough" (I can't imagine any other reason). Did anyone have the same situation? Do you know how to fix it? ![enter image description here](http://i.stack.imgur.com/e1PzR.jpg)


title: How to bulk import ESRI projections into Postgis

text: Where can I get a file that would allow me to bulk upload all of the ESRI projections into the spatial_ref_sys table of postgis. I know that you can get insert statements on spatialreference.org; however one must look these up one by one --and I would rather just load them all at once.


title: Using OSM Buildings in Leaflet: how to deactivate layer on loading map

text:
I'm trying to use OSMBuildings in my Leaflet map.               <!doctype html>                   <meta charset="utf-8">       <title>OSM Buildings Test</title>       <link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.7.2/leaflet.css"/>

                         All is working but in my `L.control.layers`, I note that the layer "OSM Buildings" starts activated and I'd like to have it deactivated but I dont know how to do it. You can test it Any suggestions are appreciated.


title: Most efficient way to find if a polygon falls with in the boundary of another polygon using postgis

text: I have a polygon layer that contains some polygons inside of polygons. I need to be able to identify the polygons inside of polygons and discard the polygons inside of the polygons. What is the most efficient way of doing that?


title: How can I center Openlayers div on screen?

text: I have a div centered on screen (using width:75%;height:50%;margin: 0 auto;) When I init the map div with OL , it auto aligned to the left. How can I keep having the div in the center of the window?


title: How to calculate percentage area cover using Focal Statistics

text: I am trying to produce raster outputs in ArcGIS for later use in MaxEnt to model species distributions, and want to produce variables with statistics on the percentage cover of different habitat polygons (i.e. woodland) and linear features (i.e treelines or hedgerows) at different spatial scales (e.g within a moving analysis window). I have a fairly rudimentary understanding of ArcGIS. I feel that Focal Statistics tool should be able to achieve what i want if i convert my layer inputs to raster format, however percent is not listed as an option in the Statistics field.


title: Need help in ArcGIS Java Script API : Installation & configuration with different browsers

text: We are using ArcGIS Java Script API (version 3.2, 3.0, 2.8) according to our requirement we are using various API’s.   1. Currently ESRI launched new version 3.3 so if I want to view old examples references i.e version 2.8, 3.0 how do I view the old examples (Because they removed old examples)?   2. Similarly, Can I install multiple SDK( 2.8,3.0,3.2) on single machine? In addition, I need recommended browser settings as if one of my application is not displaying the legend in IE 8.0 Therefore, what are the recommended settings for all versions (2.8 to 3.2) Any suggestions /discussions will be great !!! Thanks in advance : )


title: How to insert reversed map with Google Maps API?

text: Is there a way to put map with Google Maps API to a website, but upside down? North must be down, and south up. Is there any possible way to do it? I really need to put map that way and I need really fast solution.


title: Is it possible to create a PostGIS sequence to populate Area column of a polygon

text: I've created Postgres/PostGIS sequences before to assign unique gid values when a feature is created in QGIS but is it possible to create a sequence that will automatically populate an area column? Equally would it be possible for a sequence to be created to calculate the length of a line and even better the X Y values of a point?


title: How do I fix out of bounds geometric extent error?

text: I'm new to GIS, and I'm roughly following the instructions here: http://mapbox.com/tilemill/docs/guides/postgis-work/ to try and import data I obtained from here: http://factfinder2.census.gov/faces/nav/jsf/pages/searchresults.xhtml?refresh=t into a Postgres database and then query from Tilemill. I'm getting the following error when Tilemill tries to run the query:               Detected out of bounds geographic extent (-180,4424208.8281,-13492834.8813,85.051) for layer 'acs'. Please ensure that the SRS for this layer is correct. Its native extent is '-13693889.5411,4424208.8281,-13492834.8813,4566130.4316'       Does this error have something to do with the parameter 4269 I used in the following command I used to populate Postgres?:               shp2pgsql -s 4269 140_00 census_tracts | /usr/local/bin/psql -d default


title: Is PhD in GIS useful in non-educational career path?

text: A bit of an open debate - I am curious for your opinions. I am considering commencing a PhD in GIS which may concern the application of surveys within the natural environments. Eventually, I wish to continue my career in professional consultancy and not to pursue academia beyond these studies. I wonder whether anybody has any experience or can just share their thoughts about how useful is to have a GIS PhD in a professional industry. My concern is that people have told me that for some professions (e.g. Law), having a PhD qualification can in fact reduce your employability. How is it seen by industry employers to have a PhD in GIS? Also, are there any specific sites or tips anyone can suggest? Thanks, Mandy


title: Leaflet Marker Cluster: Popups won't open after first

text:
I built a map using Leaflet and markercluster.js, and I'm having a strange issue where after I click on the first popup, I'm unable to open any other popup. This happens whether you close the first popup manually or not. When I don't use markercluster.js, the popups appear as normal, so I know it has something to do with that library. This is the second time I have used this library; the first time I had no issues. I'm getting no errors in the console. Also, I had this issue regardless of if the .openPopup() method was attached to the marker variable. Code below: > var puetrasList = new L.MarkerClusterGroup(); >   >   >               for (i=0 ; i < dataStore.length ; i++){ >                   var a, point; >                   a = dataStore[i]; >                   lat = parseFloat(dataStore[i][0]); >                   lng = parseFloat(dataStore[i][1]); >                   title = "
" + dataStore[i][2] + > "
" + >                                       "
" + >dataStore[i][3] + "
" + >                                       "
" + >dataStore[i][4] + ", " + > dataStore[i][5] + "
" + >                                       "
" + >dataStore[i][6] + "
"; >                   marker = new L.Marker(new L.LatLng(lat, lng)); >                   marker.bindPopup(title).openPopup(); >                   console.log(marker) >                   puetrasList.addLayer(marker); >               } >   >   >   >               //----- Summon Map ----> >               map.addLayer(puetrasList); >               map.addLayer(layer); >


title: OSM, Displaying a road with a line on a map from selected nodes

text: Along with a friend, we have started an educational side-project,for fun not as a commercial or trial software, to test what we have learned these last few months. We decided we would do a simple map navigation program that would take an input from a user (starting place + desired location) and draw the route based on the Open Street Maps data we have extracted.   For our first attempt we assume the user is in a car and implemented in Java a k-shortest path to find the route between them. I also used a short algorithm to find the nearest node to a user's present and desired location, which while not perfect, works neatly for the level we desire to achieve. And it's working really well, the routes given over to all test maps are always valid and generally fast, a few times even optimal.   The results are returned to an ArrayList of Nodes(id,ref,long,lat). What we now want is to be able to input the starting-node,finishing-node and all the routing nodes as a single road in a map view in a web page.   Our general idea is that the user will be able to click two points at the map , and a red line will be drawn along the nodes-ways our Java class has returned. How do you propose we do this? I have tried leaflet and js, have only managed to center the user's desired map and zoom level, as well as a central point, but the line drawing has proven problematic for more than 2 nodes, either failing to display entirely or being severly out of road or in wild shapes. This has lead me to think I am taking an entirely wrong approach. Could anyone please give me a suggestion as to what resources I should be looking at?   Thank you. As for the code, I have to work off my head now. I seem to be missing something right now, but the general idea was as follows for the prototype:               var map;     var ajaxRequest;     var plotlist;     var plotlayers=[];          function initmap() {         map = new L.Map('map');         var osmUrl='http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png';         var osmAttrib='Map data © OpenStreetMap contributors';         var osm = new L.TileLayer(osmUrl, {minZoom: 12, maxZoom: 15, attribution: osmAttrib});                    map.setView(new L.LatLng(37.929, 23.6975),14);         map.addLayer(osm);              var roadLayer = new OpenLayers.Layer.Vector("Roads");               map.addControl(new OpenLayers.Control.DrawFeature(roadLayer, OpenLayers.Handler.Path));                                                  var nodes = new Array(           new OpenLayers.Geometry.Point(37.92, 23.69),           new OpenLayers.Geometry.Point(37.93, 23.70),           new OpenLayers.Geometry.Point(37.92, 23.70)      );          var roadLine= new OpenLayers.Geometry.LineString(nodes);          var roadStyle = {        strokeColor: '#0000ff',        strokeOpacity: 0.5,       strokeWidth: 5     };     var drawnLine = new OpenLayers.Feature.Vector(roadLine, null, roadStyle);     lineLayer.addFeatures([drawnLine]);     map.addLayer(roadLayer);     }      The array of Geometry points is the Array that is returned by the Algorithm, but for the example I used a few static points. In the actual code, the line is returned severly out of bounds for a few nodes, while it is not returned at all after the nodes increase. However for 2-3 nodes it is working correctly. That's why I assume my approach is bad.


title: Unable to Add Features to PostGIS Layer in Quantum GIS 2

text: I am testing out Quantum GIS 1.9-master and I've noticed that I am only able to add features to a PostGIS 2.0 layer if that layer is in a schema with the same name as my login. For example, if my login names is "brett", I am able to add features to the table brett.test_table, but not fred.test_table (regardless of the table permissions (I do have INSERT, UPDATE, DELETE permissions on tables in both schemas). Strangely enough, I am able to perform updates to features in tables in other schemas (i.e. I can move or reshape a feature). I am able to add features fine to tables in other layers using Quantum GIS 1.8 (so I don't think it's a PostgreSQL permissions issue). Has anyone else experienced a similar issue? I'm wondering if this is a bug, or a new way of handling editing in Quantum GIS 1.9.


title: using "abs()" function in QGIS field calculator?

text: I'm pretty new to GIS in general and I'm trying to take an absolute value of the difference between two fields in field calculator when creating a new field, but using ABS () returns an invalid expression in Preview. ABS () is not listed on the supported expressions--is there a way around this? Thanks! John


title: How to create an ant task to build ArcGIS for Flex Viewer?

text: Has anybody created ant task to build ArcGIS for Flex Viewer? I want to put ArcGIS Flex Viewer into our Continuous Integration system (Hudson). So far it's a very difficult task to do.


title: Leaflet draggable Marker

text: I have setup a draggable marker on my map with Leaflet API. Its working fine. But the problem is when I am trying to drag it to the border area the map is not moving. I want a panning map with Marker dragging. Please let me also know how can able to control the map panning when dragging the marker. Thanks


title: Resources to integrate ESRI gis to SAP ?

text: I heard to this term to integrate GIS with SAP.What does this actually mean? Where can i find resources to integrate GIS to SAP. I have no idea on where to start and to proceed further?


title: Adding point shapefile to map shuts down program

text: I am trying to add a point file that I have used before to a map. Every time I add it the program (ArcGIS) crashes. I can add other files, including the same data for a different year, but this specific file crashes the program, so I assume there is something wrong with the file. It would be very difficult to recreate the file, so I need to diagnose what might be wrong and then repair the file...if possible.


title: Kriging with SAGA and QGIS

text: I have the "processing" toolbox installed on QGIS2 and it would appear to give me access to the SAGA variogram tools and and kriging tools. What I am missing and what I can't find via Google either, is some sort of help on how to use these packages. The variogram tool dumps out a csv file containing what appears to be an experimental variogram for multiple directions plus covariance etc. The kriging tool has stumped me, what is a saga grid and why if I create one (which seems for all the world to be a tiff) won't the kriging tool find it? It would appear that I can specify the nugget and sill but as the variogram tool isn't overly helpful (at least as far as I can tell without any help files) then I can't know what values to use. The SAGA site has a users guide but for most tools it seems to do nothing more than list their existence. I have also found a presentation of slides showing SAGA in use but as the windows in native SAGA don't match the QGIS implementations the slides aren't very helpful. Does anyone know of a good (or even not very good) documentation? Can anyone throw some light onto this?


title: Are there any publicly available ArcSDE Servers I can use for developing client applications

text: I have been asked by our client if I can add the ability to connect to instances of ArcSDE servers from our desktop application. I am sure that we can do this a number of different ways, what I was looking for was a server I could connect to whilst doing a bit of research / initial development. We currently don't have any ArcGIS licences or products at present. Cheers Andy


title: GeoServer GetCapabilities Not Working

text: I am trying to learn how to use the OpenGeo suite by following the tutorials on the site. The one I'm currently trying to follow is http://workshops.opengeo.org/geoext/stores/capabilities.html. Most of it works, but I can't get the grid to populate with the layers from GeoServer. This is my code:               <html>         <head>             <title>GeoExt Workshop Application</title>             <link rel="stylesheet" type="text/css" href="ext/resources/css/ext-all.css">             <script type="text/javascript" src="ext/adapter/ext/ext-base.js"></script>             <script type="text/javascript" src="ext/ext-all.js"></script>             <script src="openlayers/lib/OpenLayers.js"></script>             <script type="text/javascript" src="geoext/lib/GeoExt.js"></script>                  <script type="text/javascript">                  Ext.BLANK_IMAGE_URL = "ext/resources/images/default/s.gif";             var app, items = [], controls = [];                  Ext.onReady(function() {                 app = new Ext.Viewport({                     layout: "border",                     items: items                 });             });                  items.push({                 xtype: "gx_mappanel",                 ref: "mapPanel",                 region: "center",                 map: {                     numZoomLevels: 19,                     controls: controls                 },                 extent: OpenLayers.Bounds.fromArray([                     -122.911, 42.291,                     -122.787,42.398                 ]),                 layers: [new OpenLayers.Layer.WMS(                     "Medford",                     "http://localhost:8080/geoserver/wms",                     {layers: "medford"},                     {isBaseLayer: false}                 )]             });             controls.push(                 new OpenLayers.Control.Navigation(),                 new OpenLayers.Control.Attribution(),                 new OpenLayers.Control.PanPanel(),                 new OpenLayers.Control.ZoomPanel()             );                  items.push({                 xtype: "grid",                 ref: "capsGrid", // makes the grid available as app.capsGrid                 title: "Available Layers",                 region: "north",                 height: 150,                 viewConfig: {forceFit: true},                 store: new GeoExt.data.WMSCapabilitiesStore({                     url: "http://localhost:8080/geoserver/wms?SERVICE=WMS&REQUEST=GetCapabilities&VERSION=1.1.1&",                     autoLoad: true                 }),                 columns: [                     {header: "Name", dataIndex: "name", sortable: true},                     {header: "Title", dataIndex: "title", sortable: true},                     {header: "Abstract", dataIndex: "abstract"}                 ],                      bbar: [{                     text: "Add to Map",                     handler: function() {                         app.capsGrid.getSelectionModel().each(function(record) {                             var clone = record.clone();                             clone.getLayer().mergeNewParams({                                 format: "image/png",                                 transparent: true                             });                             app.mapPanel.layers.add(clone);                             app.mapPanel.map.zoomToExtent(                                 OpenLayers.Bounds.fromArray(clone.get("llbbox"))                             );                         });                     }                 }]             });                  </script>         </head>         <body>         </body>     </html>      I've been playing around on firebug to try to figure out what's going on, particularly on the request               http://localhost:8080/geoserver/wms?SERVICE=WMS&REQUEST=GetCapabilities&VERSION=1.1.1&      since this is the request that is supposed to get the layer info. The request header it sends is below. The response for some reason is 0 bytes. I'm not very familiar with web development, but I thought this should be a GET request instead of an OPTIONS request.               OPTIONS /geoserver/wms?SERVICE=WMS&REQUEST=GetCapabilities&VERSION=1.1.1& HTTP/1.1     Host: localhost:8080     User-Agent: Mozilla/5.0 (Windows NT 6.1; WOW64; rv:18.0) Gecko/20100101 Firefox/18.0     Accept: text/html,application/xhtml+xml,application/xml;q=0.9,*/*;q=0.8     Accept-Language: en-US,en;q=0.5     Accept-Encoding: gzip, deflate     Origin: null     Access-Control-Request-Method: GET     Access-Control-Request-Headers: x-requested-with     Connection: keep-alive      Here is the response header:               HTTP/1.1 200 OK     Allow: GET, HEAD, POST, PUT, DELETE, TRACE, OPTIONS     Content-Length: 0     Server: Jetty(6.1.8)      What is going wrong? And how can I fix it? EDIT: Here is part of the log file, which includes starting up OpenGeo and then loading the webpage I built in the tutorial. The rest is just the requests for the map tiles that look like               127.0.0.1 -  -  [23/Jan/2013:18:10:52 +0000] "GET /geoserver/wms?LAYERS=medford&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&FORMAT=image%2Fjpeg&SRS=EPSG%3A4326&BBOX=-122.87109375,42.36328125,-122.8271484375,42.4072265625&WIDTH=256&HEIGHT=256 HTTP/1.1" 200 6993      I went over the character limit by including them so I took them out.               127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "HEAD /dashboard/version.ini HTTP/1.1" 200 0      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces.json HTTP/1.1" 200 685      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/AENVTest/datastores.json HTTP/1.1" 200 155      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/AENVTest/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/medford/datastores.json HTTP/1.1" 200 262      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/medford/coveragestores.json HTTP/1.1" 200 178      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/earth/datastores.json HTTP/1.1" 200 268      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/earth/coveragestores.json HTTP/1.1" 200 166      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/opengeo/datastores.json HTTP/1.1" 200 146      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/opengeo/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/usa/datastores.json HTTP/1.1" 200 148      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/usa/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:43 +0000] "GET /geoserver/rest/workspaces/analytics/datastores.json HTTP/1.1" 200 138      127.0.0.1 -  -  [23/Jan/2013:18:10:44 +0000] "GET /geoserver/rest/workspaces/analytics/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:44 +0000] "GET /geoserver/rest/workspaces/world/datastores.json HTTP/1.1" 200 252      127.0.0.1 -  -  [23/Jan/2013:18:10:44 +0000] "GET /geoserver/rest/workspaces/world/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:44 +0000] "GET /geoserver/rest/layers.json HTTP/1.1" 200 3036      127.0.0.1 -  -  [23/Jan/2013:18:10:44 +0000] "GET /geoexplorer/maps/ HTTP/1.1" 200 556      127.0.0.1 -  -  [23/Jan/2013:18:10:45 +0000] "HEAD /dashboard/version.ini HTTP/1.1" 200 0      127.0.0.1 -  -  [23/Jan/2013:18:10:48 +0000] "HEAD /dashboard/version.ini HTTP/1.1" 200 0      127.0.0.1 -  -  [23/Jan/2013:18:10:50 +0000] "HEAD /dashboard/version.ini HTTP/1.1" 200 0      127.0.0.1 -  -  [23/Jan/2013:18:10:52 +0000] "OPTIONS /geoserver/wms?SERVICE=WMS&REQUEST=GetCapabilities&VERSION=1.1.1 HTTP/1.1" 200 0      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "HEAD /dashboard/version.ini HTTP/1.1" 200 0      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces.json HTTP/1.1" 200 685      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/AENVTest/datastores.json HTTP/1.1" 200 155      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/AENVTest/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/medford/datastores.json HTTP/1.1" 200 262      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/medford/coveragestores.json HTTP/1.1" 200 178      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/earth/datastores.json HTTP/1.1" 200 268      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/earth/coveragestores.json HTTP/1.1" 200 166      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/opengeo/datastores.json HTTP/1.1" 200 146      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/opengeo/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/usa/datastores.json HTTP/1.1" 200 148      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/usa/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/analytics/datastores.json HTTP/1.1" 200 138      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/analytics/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/world/datastores.json HTTP/1.1" 200 252      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/workspaces/world/coveragestores.json HTTP/1.1" 200 21      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoserver/rest/layers.json HTTP/1.1" 200 3036      127.0.0.1 -  -  [23/Jan/2013:18:10:53 +0000] "GET /geoexplorer/maps/ HTTP/1.1" 200 556      127.0.0.1 -  -  [23/Jan/2013:18:10:55 +0000] "HEAD /dashboard/version.ini HTTP/1.1" 200 0


title: QGIS Raster Calculater Syntax

text: I am trying to use the raster calculator in QGIS to combine a number of rasters. I want to calculate a mean, this is easy:               (raster_1@1 + raster_2@1 + raster_3@1) / 3      I then want to create a max. raster - i.e. for each pixel I want to compare raster_1, 2 & 3 and use whichever is the greatest of the three. According to the QGIS documentation I should be able to use the following code (modified from QGIS documentation example):               gt( raster_1@1, raster_2@1, raster_1@1 )      (All cells in raster_1 with value greater than (gt) the appropriate pixels in raster_2 will be replaced with the appropriate pixels from raster_1). I'd then have to repeat this for raster_3. This doesn't work however as the raster calculator says the expression is invalid. What am I doing wrong? None of the examples of conditional statements provided in the QGIS documentation seem to work. I have tested this with .tif and .sgrd files and the same thing happens. I'm using QGIS 1.8 on a Linux machine. Thanks in advance.


title: Learning Spatial Statistics

text: I have a significant background in statistics, and I'd like to learn more about spatial statistics (and statistics associated with geospatial data). Are there any especially good resources for someone looking for an overview of the field? I am not looking to solve a specific problem at this point. Just want to understand the field and how it differs with statistical applications on non- spatial data.


title: SAGA, GRASS Raster values to null

text: I want to make the cells which values "1" to null in QGIS with Grass or Saga tools (exactly Saga).


title: QGIS "transform error" after downloading from OpenStreetMap

text: Beginner's question. I am working on this tutorial: http://www.qgis.org/en/docs/training_manual/vector_analysis/basic_analysis.html After selecting WGS84 UTM zone 34S in Project Properties, I receive several of these transform error messages: Simplify transform error caught: forward transform of (0.366535, 1.570812) PROJ.4: +proj=longlat +datum=WGS84 +no_defs +to +proj=utm +zone=34 +south +datum=WGS84 +units=m +no_defs Error: latitude or longitude exceeded limits Can anyone guide me as to next steps? From other questions I've read on this site I assume it's a long/lat issue. Many thanks.


title: How do I implement Google Maps API V3 Image Map Type?

text: I'm new to the Google Maps API V3 and I am trying to add a set of map tiles via the image map type. I followed the steps in this article, and have had no dramas creating the Tileset. I used GDAL2TILES and the the tilesets looks as expected in the preview apps that are generated with the program. However I cannot get the tilset to show up in my google maps application. My implementation of the Image map type looks like this below                   //The index.html and 'sma_map_1989' folder are in the same directory          var smaMapOverlay = new google.maps.ImageMapType({           getTileUrl: function(coord, zoom) {             return 'sma_map_1989' + '/' +zoom+ '/' +coord.x+ '/' + coord.y +'.png';           },           tileSize: new google.maps.Size(256, 256)         });         map.overlayMapTypes.push(smaMapOverlay);      Am I missing something here? It's such a simple thing to do in leaflet and Openlayers.


title: changing layout of a gxp.plugins.WMSGetFeatureInfo response

text: Have a little problem with an app I’m building. I´m using the openGeo suite SDK to create an app. I added a WMSGetFeatureInfo tool following the Boundless tutorial and to that point it was all ok. Then I decided to change the Templates to get just the meaningful part of the data. I made a Template following instructions here. If I preview the layers using OpenLayers and click any feature I get the right answer but then if I go to my app and use the WMSGetFeatureInfo tool it ignores the config and shows me all the fields in the layer. What I’m I doing wrong?


title: Tasseled Cap Transformation Rapid Eye

text: I want to do tasseled cap transformation on rapid eye data. the main question would be about the parameters which are used. Somebody has a suggestion?


title: Zooming to specific polygon boundary Openlayers

text: I want to zoom to specific polygon boundary of some polygon layer. How can I do this.The map layer is WMS layer.


title: More than one WMS layer in QGIS

text: I have a problem when using WMS layers. I'm usin QGIS. I'm loading a WMS layer containing some aerial images. And then on top of it another WMS layer that contains information of the city, boundaries etc with labels. When I load the WMS layer with the labels, the background goes blank and I can't visualise the aerial images. I've tried the same thing in Arc and works fine but I have to use QGIS. Cheers.


title: How to convert Direction Of Travel OSM vectors from a "system" to another?

text: I have a painful doubt according Direction of Travel (DOT) vectors in OSM data. I'm using ArcGIS10, and a OSM data shapefile. I have a OSM road layer, that has many "Both ways" vectors and many "One way". I'll consider only One way for my study. Its DOT was performed as it was created the vector in its origin: The DOT is set from the start point of the vector to its finishing point. But, **I want to change this system to mine** , which is:   * Each vector has a Reference Node and a Non-Reference Node.    * The Reference Node is the **southernmost node** and if the latitudes of both end nodes are identical and their longitudes differ, the Reference Node is **the westernmost node**. An example: ![Example](http://i.stack.imgur.com/xjOJm.png) So, if I have a One way vector, I'm going to add its DOT with a letter T if it goes to the Reference Node or with a letter "F" to the Non-Ref node. An example: ![Example DOT](http://i.stack.imgur.com/z9T9Z.png) Considering that OSM one way vectors are set in another "system", **how can I have it into my system?** Thanks for your appreciate time.


title: Why is it so hard to download only a part of OSM data by attribute tag?

text: Is there any online service that allows to download only one thematic layer of the OSM data for a specific geographic region? For example all objects related to vegetation (trees, grassland,...), but not streets, railways etc. There must be an easy and straightforward way to do so (I hope?). I'd prefer to use the data directly in QGIS. I am aware of Geofabrik which provides some data packages. But I am interested if there is any way to specify the data I want before actually downloading it?


title: Communication between widgets (flex)

text: Default Communication between widgets hi, i'm new in flex, and im working on Webmapping application using flex 4.6, java, and sql server as database, i'm just trying to make a communication between wigdets, i would like to copy each selected item of my combobox in widget1, and show it in the textinput in widget2, i just followed this steps explaind by Robert here http://forums.arcgis.com/threads/138...+communication but no luck, the problem that i'm facing is that i can't call components of the first widget (textinput, button....) from the second widget, they are often underlined !! at jump you found the change i made, Widget 1               var layoutArr:ArrayCollection = new ArrayCollection();                 private function widgetCommunication1():void                 {                     addSharedData("Liste des chapitres", layoutArr);                  }      Widget 2               icationCompleteHandler(event)"                         initialize="basewidget1_initializeHandler(event)">          <fx:Script>             <![CDATA[                                  import com.esri.viewer.AppEvent;                 import mx.events.FlexEvent;         protected function basewidget1_initializeHandler(event:FlexEvent):void                 {                     ViewerContainer.addEventListener(AppEvent.DATA_PUBLISH, sharedDataUpdated);                 }           private function sharedDataUpdated(event:AppEvent):void                 {                     var data:Object = event.data;                                    if (data.key == "Liste des chapitres")                          {                         dplist.selectedIndex = ***IntChap***.text;                     }                  }             ]]>         </fx:Script>          <s:DropDownList id="dplist" x="10" y="10" width="282" height="26" labelField="IntituleChap" dataProvider="{RemoteRapportDao.Selection.lastResult}" change="dplist_changeHandler(event)" ></s:DropDownList>      IntChap which is bold & italic is the textinput in widget1 , i got it underlined in widget 2, even if i made this communication any one can tell me what i'm messing ?? any help or code will be appreciated !


title: GRASS - v.clean Conflation - How to snap a layer to a base layer

text: Someone told me in this post that GRASS GIS can handle multilayer topology, what I have not found any information about. I am not a specialist of GRASS GIS so I would need help with that if its really feasible. What I am trying to do is to snap a layer to a base layer, so when both share the sames boundaries they would be topologically snapped. The base layer is the official cities boundaries of my country. The layer to be snapped is made of cities subdivisions but may not always fit. ![enter image description here](http://i.stack.imgur.com/TGA1e.png) I have to be sure that the base layer would no be modified at all.


title: GeoTiff TMS Overlay

text: I'm using OpenLayers to display a set of TMS tiles I have on my webserver. Plain and simple, I was wondering how I would go about overlaying a set of GeoTiff files I have over the tiles so that when I click somewhere on the map, it would give me the information from the tiff files (altitude and such).


title: How to create the fuel models & live fuel moisture rasters for Wildfire simulation?

text: Has anybody done wildfire modelling in GRASS? Where could I obtain / how do I create the fuel models and live fuel moisture rasters to determine rate of spread (r.ros)?


title: Using arcpy to enable layer masking (Advanced Drawing Options)

text: I am looking for a way to enable and change layer masking (under Advanced Drawing Options) with arcpy... it seems you can do this if you have the Production Mapping extension, but I do not have it... anyone know of a way to do this (just using ArcGIS for Desktop Advanced 10.2.2)?


title: Renaming a field using arcpy

text: What I need to do:   1. rename a field name of a table/feature class   2. copy all values to the new field So far I have done following code as artwork21 suggests: **EDIT:**               import sys     import traceback     import arcpy     from arcpy import env          ## ARGUMENTS     # argv[1] = input table/feature class path     # argv[2] = input old field name     # argv[3] = input new field name            path = sys.argv[1]     oldFieldName = sys.argv[2]      newFieldName = sys.argv[3]           env.overwriteOutput = True          fields = arcpy.ListFields(path)     for field in fields:         if field.aliasName == oldFieldName:             if not oldFieldName == newFieldName:                 fieldType = field.type                 # Add new field                 arcpy.AddField_management(path, newFieldName, fieldType)                 #Calculates the new field based on old field values                 arcpy.CalculateField_management(path, newFieldName, "!"+oldFieldName+"!", "PYTHON", "")                 # Delete the old field (if necessary)                 arcpy.DeleteField_management(path, oldFieldName)      **Question:** How can I map `field.type` to `AddField_management` method's field type? And while the field is in a middle place, then field deleted from middle and added to the last. That does not looks like as if field name is renamed. Is there any better solution that help me to do these things? Thanks