Setting transparency in an SLD file for a 4 channel GEOTIFF

I'm struggling to set the a background colour to transparent in an SLD file for a 4-channel GEOTIFF using GeoServer 2.1.4. I think I've set my background colour to transparent alpha layer in a graphics package correctly. I've set up an SLD file to read Red/Green/Blue channels of my GEOTIFF = this resolves colours correctly for the RGB image = but still shows the background black colour (see below). When I do a GetFeatureInfo request on the WMS, I can see an ALPHA_BAND=0 where I want the transparency and ALPHA_BAND =255 where I do not want transparency (see attached figure). ![enter image description here](http://i.stack.imgur.com/DDH4g.png) I was hopeing to do something like setting the ALPHA_BAND=0 to Opacity=0 in the SLD or something. However, there's no option in SLD schema that I can see ? I've read the SLD cookbook and the SLD reference documentation. Alpha transparency is referred to but only with ColorMap' property. Can I mix and match 'ColorMap' and ChannelSelection' properties in the SLD ? Am I going about this the right way ? Thanks in advance for assistance! This is the SLD I'm using:                   <?xml version="1.0" encoding="ISO-8859-1"?>         <StyledLayerDescriptor version="1.0.0" xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc"           xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"           xsi:schemaLocation="http://www.opengis.net/sld http://schemas.opengis.net/sld/1.0.0/StyledLayerDescriptor.xsd">           <NamedLayer>                 <Name>bathy example</Name>                 <UserStyle>                   <Name>RGB image 3 channels</Name>                   <Title>RGB image</Title>                   <FeatureTypeStyle>                         <Rule>                           <RasterSymbolizer>                                 <ChannelSelection>                                   <RedChannel>                                   <SourceChannelName>1</SourceChannelName>                                   </RedChannel>                                   <GreenChannel>                                   <SourceChannelName>2</SourceChannelName>                                   </GreenChannel>                                   <BlueChannel>                                   <SourceChannelName>3</SourceChannelName>                                   </BlueChannel>                                             </ChannelSelection>                           </RasterSymbolizer>                         </Rule>                   </FeatureTypeStyle>                 </UserStyle>           </NamedLayer>         </StyledLayerDescriptor>


Remove duplicate names in FeaturePopups cluster tooltips and popups

I am using FeaturePopups control to show tooltips and popups for clusters. I have features connected (ends of road segments) with the same names. Is it possible to remove them from tooltips and popups, so to have a list with unique names? Thank you, Bogdan


Line intersection with circle on a sphere (globe or earth)

How can I find the intersection point of a line and a circle on a sphere (on earth) Please, if any one knows a good solution help. Thanks in response :) ![Imagine this on a globe or sphere](http://i.stack.imgur.com/yyxPS.png)


Does MySQL only support bounding boxes calculations instead of the real polygon geometry?

Within our company we are deciding which DBMS we will use to store spatial data. IT wants MySQL with the Spatial extension because there is already knowledge about handling MySQL DB. We are doing quite simple calculations with the data. The only thing I am worried about is the polygon handling of MySQL. Is it still true that MySQL only supports bounding box calculations instead of the real polygon geometry? Or is this problem gone with newer versions of MySQL?


How can I find start and end points from a road layer?

I have a shapefile with geometries MultiLineString and Point. I have loaded the shapefile into PostGIS. From that table I cannot find any (lon,lat) except geometrical values with some SRID,.... columns. How about knowing of nodes, edges, start point and end point for a particular layer that has above geometries; and to create a start point and end point from a shape file I have tried few shapefiles, finally stopped with               Error:Start vertex was not found.


Using log4net in ArcMap addin solution

I would like to use log4net as the logging framework for my addin. However, now I'm having trouble in enabling logging in debug mode **and** in the released esriAddin file. It works in either one, not both at the same time. I have a dedicated configuration file for log4net (called log4net.xml). When the **Build Action** of the file is set to `Content`, the logging works in the debug mode but not in the released addin. When Build Action is `AddinContent`, it's the opposite. The reason for this is that when in debug mode, the solution uses the bin folder as the base directory, but when the addin is compiled the folder is the "esriaddin" file. So when in debug mode and `Build Action`=Content" the log4net.xml is at `C:\Users\%username%\AppData\Local\ESRI\Desktop10.0\AssemblyCache\{9FFCCA4E-6DD6-81B5-17C3-7683C725F70D}`. When I change Build Action to "Addin Content" (while still in debug mode), the log4net.xml is not copied to the directory! My problem is how can I configure my solution so that the log4net config file would be in a **static location** regardless of the Build Action? I need to do this because I want to automate the building process and changing the **Build Action** manually is not an option. I've checked this thread but its solution is to edit the `ArcMap.exe.config` file which is not an option for me (many users for the addin, the file may not be empty). Any suggestions?


How To find the Catchment Area for a Dam?

I have a DEM of my study area, and polyline shapefile of the Dams in that area. I have run the standard Hydrology tools on it( Fill, Flow Direction & Flow Accumulation). I know that you can find the Watershed contributing to a pour point, but how do you use the Damlines, to get the catchment Area contributing water to the dam's reservoir?


3D Building Generator Tool

I am using ArcGIS desktop 10 and I would like to change a 2D polygon into 3D building. Does anyone have any idea from where I can download the **3D building generator tool** for ArcGIS Desktop? I read about this tool here. This tool gives various features and can build complex roof tops for building sites. With this tool I can develop a rounded roof in 3D ..etc Moreover, any suggestions how I can import this 3D buildings in a web service using asp.net and C#? Any help please?


ESRI ArcMap Topology rules on multi story buildings

I would like to run some topology rules to find snapping errors in our polygons but I am running in to the problem where we are working on multi story buildings with every polygon regardless of floor is stored in one feature class... We separate the floors out by adding layers with definition queries... From what I can see I can only create topology rules (like "boundaries must be covered by") for a feature class and not a queried layer... In the screen shot provided the grey polygon is not snapped correctly to the North West most corner of the only polygon completely in frame of the screenshot... Running a boundary covered by will easily spot this error but only if I disassemble the feature class so that each floor is on its own feature class... Needless to say the more floors involved the more tedious this becomes... Not to mention that I must then reassemble the separated feature classes back into one as this is how the deliverable is required to be... There is one feature class for the room polygons and one for the CAD line work... There must be an easier way to find gaps like this without separating the feature classes to run the topology then putting the separated classes back together to produce the deliverable... Any help would be greatly appreciated... Please ask if you need further clarification... ![Screenshot of a problem area](http://i.stack.imgur.com/qGgZP.jpg)


Geoprocessor Execute HRESULT COM Error

Recently, I started having problems with my C# projects running ArcGIS10.1 geoprocessing tools. My tools point to custom scripts in project toolboxes. The code in question has run without issue in the past. However, something in my environment has changed and even the simplest code throws errors. Please see sample below. The code fails on `gp.Execute()`. I've triple checked, labels, names, paths, parameters etc. I've tried every thing I can think of to track down this problem including re- installing Arc. I'm dead in the water if I can't solve this nasty problem. Thanks in Advance                   private void runTool()         {             try             {                      Geoprocessor gp = new Geoprocessor();                 gp.AddToolbox(@"C:\0_Test\Test.tbx");                      IVariantArray parameters = new VarArrayClass();                 parameters.Add(@"MyStringArg");                      // Execute the tool.                 gp.Execute("test", parameters, null);                  }             catch (Exception ex)             {                 System.Windows.Forms.MessageBox.Show(ex.Message, "TestGP Error", System.Windows.Forms.MessageBoxButtons.OK);                 return;             }         }


gdal_translate produces tif with different colors as original

I am using SDE as my raster data source with the ultimate goal to tile this dataset: I am running the following gdal_translate command to generate the initial VRT:               gdal_translate  -projwin 390605 5588320 410128 5568805 -a_srs EPSG:25832 -a_nodata none -of vrt SDE:server,5153,db,user,pw,table,image clippedvrt.vrt      Tiling the resultant VRT produces tiff with weird colors. The colors are either intensified or totally different from the original as you can see from the test results Here is the result of gdalinfo on the resultant VRT:               Band 1 Block=128x128 Type=Byte, ColorInterp=Palette     Color Table (RGB with 256 entries)     0: -1,-1,-1,255     1: 0,0,-1,255     2: -1,22873,22873,255     3: 26214,13107,0,255     4: -1,22873,22873,255     5: 26214,-17734,-1,255     6: -1,-3342,22873,255     7: 0,-6683,0,255     8: 0,0,0,255     9: -19790,26214,0,255     10: -13108,-1,-13108,255     11: 0,0,0,255     12: -11052,-11052,-11052,255     13: -13108,-1,-1,255     14: -3342,-10538,-21846,255     15: 26214,13107,0,255     16: -19276,-32640,26214,255     17: -1,-1,-1,255     18: -1,-1,-1,255     19: -1,-1,-1,255     20: -1,-1,-1,255     21: -1,-1,-1,255     22: -1,-1,-1,255     23: -1,-1,-1,255     24: -1,-1,-1,255      could someone shade some light as to what might be happening. Thank you Moses


Python script for identifying duplicate records (follow up)

I am trying to calculate a field that marks duplicates and sequentially numbers them. I am using the code below that I found here: Generating sequential numbering for duplicate values in a field sorted by ObjectID. It works, but I get a duplicate starting value before it begins sequentially numbering UNLESS I only run the field calculator on the selected values. When I calculate after manually selecting the duplicate rows and running it on only those values, it works fine. Here's the code:               prevFieldValue = ''     counter = 0     def GetDuplicateCounter(myFieldValue):       global prevFieldValue       global counter       if myFieldValue == prevFieldValue:         counter += 1       else:     counter = 1       prevFieldValue = myFieldValue       return counter      Results running it on all roads (nothing selected): ![calculated without selection](http://i.stack.imgur.com/WvUY0.jpg) Results running it after manually selecting the duplicated fields (not too efficient!). ![calculated on selected records](http://i.stack.imgur.com/EeLdp.jpg) Any help would be appreciated! Thanks! b


Add image from my computer for background in Print Composer?

Add image either by dropdown or icon does not produce any kind of dialog box where I could upload an image. Item properties is checked and showing in the right hand panel, but it's blank (no options showing at all). Do I need to set properties for this somewhere else (don't see how to do it inside Print Composer)? QGIS is new to me. It looks to be a powerful suite of tools. Just can't get this working.


Can't connect to WMS service? An Amazon-EC2 problem?

I'm running GeoServer 2.3.5 on a Windows instance in AWS for the first time. The data is on the server in PostGIS. The map layers have been published and I can see them locally without any problems. However, when I try to access the WMS service for these layers from another machine, I can't connect. I'm wondering if there is some configuration I may be missing because this is set up through AWS. Any ideas?


Trying to use Python code and user defined inline parameters to build a SQL expression

I am trying to build a tool where the end user can run a model and make a selection from a feature class based on 3 different fields, calculate the total shape length, & export the results to a csv file. I have everything else figured out except for one part. I want the user be able to make a selection based on all three fields **or** based on one , two, or all three fields without them having to edit the selection query. As you probably guess when I don't populate one of the user defined parameters it pops up with an error. I think the way to do this would be to use the "calculate value"in model builder and use Python to build a sql expression and use that expression for the selection. I was thinking of having 3 user defined parameters with drop down list of field values and using inline variable substitution in the calculate value tool to create a sql expression based on what the user defined in the parameters. This is the expression I have right now for them to select using all three field values:               "DISTRICT" = '%District%' AND "SURFACE" = '%Surface%' AND "JURISD" = '%Jurisdiction%'      I have a basic understanding of Python. Any help or suggestions would be greatly appreciated. Thanks!


Request WFS from GeoServer with Jquery AJAX using INTERSECT cql_filter error

I am having difficulty troubleshooting this issue and would appreciate some ideas. Thanks in advance! I am trying to make an AJAX request from a Leaflet map with Jquery to GeoServer to return a WFS selection using the INTERSECT sql_filter param. I have it working fine using the BBOX sql_filter. **Successful request using BBOX filter**                  var geojsonLayer = new L.GeoJSON();        var currSelectionLayer = 'DRCOGPUB:rea_drcog_jurisdiction_county';        var sw_lat = 39.93290692296977;        var sw_long = -105.39459228515624;        var ne_lat = 40.094882122321174;        var ne_long = -105.150146484375;                  function loadGeoJson(data) {                 console.log(data);                 geojsonLayer.addData(data);             }        var geoJsonUrl = "http://gis.drcog.org/geoserver/DRCOGPUB/ows?service=WFS&version=1.0.0&request=GetFeature&typeName="+ currSelectionLayer +"&srsName=EPSG:4326&maxFeatures=10000&outputFormat=json&format_options=callback:getJson&cql_filter=BBOX(the_geom," + ne_long + "," + ne_lat + "," + sw_long + "," + sw_lat + ")";              $.ajax({                 url: geoJsonUrl,                 dataType: 'jsonp',                 jsonpCallback: 'getJson',                 contentType: 'application/json',                 success: loadGeoJson             });        map.addLayer(geojsonLayer);      However, when I apply an INTERSECT filter I get this error: > Resource interpreted as Script but transferred with MIME type text/xml Here is my code for trying to get the INTERSECT filter to work: **Un-successful request using INTERSECT filter**                       var geojsonLayer = new L.GeoJSON();             var currSelectionLayer = 'DRCOGPUB:rea_drcog_jurisdiction_county';                  function loadGeoJson(data) {                 console.log(data);                 geojsonLayer.addData(data);             }                  var geoJsonUrl = "http://gis.drcog.org/geoserver/DRCOGPUB/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=" + currSelectionLayer + "&srsName=EPSG:4326&maxFeatures=10000&outputFormat=json&format_options=callback:loadGeoJson&cql_filter=INTERSECT(the_geom,%20POINT%20(-105.36369323730467%2040.10013461308659)";             $.ajax({                 url: geoJsonUrl,                 dataType: 'jsonp',                 jsonpCallback: 'getJson',                 contentType: 'application/json',                 success: loadGeoJson             });             map.addLayer(geojsonLayer);      I have tried many variations on the AJAX request and tweaking my Apache config but have had no luck. Any ideas would be appreciated. Please let me know if you need more info. **UPDATE:** I enabled the mime_module (mod_mime.so) in Apache and added the content type with this config code:               AddType text/javascript .js      **UPDATE - Final Fixed Code:**                       var geojsonLayer = new L.GeoJSON();             var currSelectionLayer = 'DRCOGPUB:rea_drcog_jurisdiction_county';                  function loadGeoJson(data) {                 console.log(data);                 geojsonLayer.addData(data);             }                  // using INTERSECT strategy             var geoJsonUrl = "http://gis.drcog.org/geoserver/DRCOGPUB/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=" + currSelectionLayer + "&srsName=EPSG:4326&maxFeatures=10000&outputFormat=json&format_options=callback:getJson&cql_filter=INTERSECTS(the_geom,%20POINT(-104.986145%2039.733066))";             $.ajax({                 url: geoJsonUrl,                 dataType: 'jsonp',                 jsonpCallback: 'getJson',                 contentType: 'application/json',                 success: loadGeoJson             });             map.addLayer(geojsonLayer);


how does a personal-gdb work from a MS Access perspective?

I am trying to get a better general understanding of how ArcGIS works with MS Access when using a personal geodatabase. If I create a new pgdb and then import a very simple shapefile and finally look at the result through Windows explorer (rather than ArcCat), the only visible item from the pgdb creation process is one mdb file. Opening the mdb file in Access shows there are 11 tables (in the case of a shp file that had one point). Generally, I can find very little documentation aboutlooking at a pgdb through Access, and if there is any sense in doing so (even if it is just for the benefit of trying to find out how it all works). Most resources about pgdb design don't go say too much compared the the other two recommended gdb alternatives. However, Access is still extremely popular, so I think there is some benefit to looking at the topic in more detail. So is it possible to know how the mdb tables that ArcGIS creates are related, or are the relations defined by files that neither Access nor Windows can see? Is it possible to know how ArcGIS uses these tables to provide the "spatial" component to Access? Is it ever benefitial to use MS Access to work with a shp file stored as a pgdb? Or is this all futile because pgdb technology is proprietary and no detailed publicly available knowledge exists? Any insight/links that answer the above and related types of questions would be appreciated.


how to learn arcgis server programming?

I'm new in programming c# and and arcgis server . i just use arcgis server wizard to create services and web applications but now i'd like to learn more about arcgis server and programming with asp.net . I just found programming arcgis server with asp.net book. i didn't buy this book .i coudln't find any website about learning arcgis server . could you please recommend another resources to learn arcgis server programming for beginners?


How to automatically place parcels on suitability map?

I need to make a map that randomly places parcels (groups of pixels or polygons) on a suitability map by first filling the most suitable areas, then the next most suitable and so on, until a pre-defined number of parcels have been placed on the map. This is probably best explained with an example. I have a shapefile with polygons on it (counties) and each polygon has a value (number of acres for a crop). I also have a suitability map (0-9 suitability for that crop). I need to place the number of acres defined in each polygon onto the suitably map. I want to first fill the most suitable lands (9) and if there is still land left fill the next most suitable land (8) and so on until the number of acres for that county are distributed. Each polygon has a different number of acres and a different range of suitability, e.g. some will have 0-9, others might only have 3-7. Within each suitability layer the placing of the parcels can be random. I have just started looking into Python scripting for ArcGIS. It seems like I might be able to do this with Python but I'm not really sure where to start. I would be interested in hearing any suggestion or advice on how to accomplish this.


Are there any common field names for the FGDC street address standard fields?

I have been looking into implementing the FGDC street address data standards so that my address data is potentially consistent with other US data I obtain in the future. However, the FGDC documentation doesn't suggest actual names to name the fields. I currently am using the following fields:               555 1 SW Main St, Ste B           - ST_NUMBR:  555      - ST_NSFX:   1      - ST_PREFIX: SW      - ST_PRETYPE:      - ST_NAME:   Main      - ST_TYPE:   ST      - ST_POSTDIR:      - ST_UNIT:   Ste B      To be consistent, I should probably change the `ST_PREFIX` field to `ST_PREDIR`, but if I'm making changes I might as well know if there are other conventions I should be following. Are there any standard or common names for such address fields?


Automating Point Data Collection

I'm supposed to go through Google Maps and manually record names and addresses of dentists in Texas for geocoding, by zip codes. This task is going to be immensely tedious and I know that there has got to be a better way to go about it. Spending two weeks recording, geocoding and verifying point locations is something I'm trying to avoid having to do at all costs but I'm lost as to what else I can do. I've thought about writing a perl or python script to either scrape the map or hook up with Google's API but can't really seem to find any advice. On top of that, scraping is against Google's TOS. Any ideas? I'm bangin' my head against my desk here. I know that this isn't a traditional GIS question but I'm not sure where else to ask.


DEM local adjustments - terraforming

I'm looking for software that can perform local adjustments in DEMs with tools similar to paintbrushes used in graphics editing programs but streamlined to respect geomorphology. To give you an idea of what I have in mind check out this video of SimCity 4 terraforming mode. I would give a lot to have such a tool in ArcGIS environment. Any local DEM modification in ArcMAP is a hassle. Let's say I have a DEM with river bathymetry imprinted in it but there is a river port along that wasn't covered and has the height of the water surface (that's the usual case with LIDAR derived DEMs - water bodies are filled). The only information about the harbour is the maximal depth. So to transfer this information to my DEM, I erase the entire harbour area (null values), imprint a small raster with maximal depth in the center of the harbour basin and then I fill the remaing null cells with a moving window filter. The results of this operation are usually sufficient but the whole process is complicated. If only I had a "lower terrain paintbrush" :)


Adding CSV upload browse button to Ext.Action

I am trying to create a CSV File upload on GeoExt Map App. I need to place the upload function within the Ext.Action, so that I can add it to the toolbar of the GeoExt Panel. I am trying to implement this example. Here is my code,               action = new Ext.Action({             text: "Upload Excel",             control: new Ext.create('Ext.form.Panel', {                     title: 'Upload a CSV File',                     width: 400,                     bodyPadding: 10,                     frame: true,                     renderTo: Ext.getBody(),                     items: [{                         xtype: 'filefield',                         name: 'csv',                         fieldLabel: 'CSV Upload',                         labelWidth: 50,                         msgTarget: 'side',                         allowBlank: false,                         buttonText: 'Select CSV File'                     }],                          buttons: [{                         text: 'Upload',                         handler : function() {                             var form = this.up('form').getForm();                             if(form.isValid()){                                 form.submit({                                     url: 'file-upload.py',                                     waitMsg: 'Uploading the CSV File...',                                     success: function(fp, o) {                                         Ext.Msg.alert('Success', 'Your csv file "' + o.result.file + '" has been uploaded.');                                     }                                    });                              }                            }                     }]             }),              map: map,             // button options             tooltip: "Upload CSV File",             // check item options             group: "newTool"         });         actions["upCSV"] = action;         toolbarItems.push(action);      Firebug keeps giving me this error, `TypeError: b[d.xtype || e] is not a constructor` Am I declaring the function incorrectly within the Ext.Action?


Where can i find free mapnik map styles with focus on waterways?

I am using Google maps for my map tiles. These tiles are very road oriented. As i am making an application for waterways i would like to change my map so that its main focus is on waterways. I have a little bit experience with making custom map tiles with Mapnik, but it still requires me to make a lot of styling decisions and styling is not my strongest point. So it would be great if there are some free stylings available that let me render map tiles with a focus on waterways, like the ANWB (Dutch) water map. I have found a site where i could change and test styles on the fly, but they had to serve the tiles and i think it also cost money to use their service. I am able to host such a tile service myself once i have the correct tiles.


Harvest data from an ArcGIS Server

I would like to include some state agency data in an application I am creating. The agency will not share the data so I have to find another approach. What I am interested in are the points on the webmap below along with the attributes of each point (onClick). I know I can easily digitize the points but would prefer not to take this approach http://fwp.mt.gov/gis/maps/huntPlanner/?bma=true&region=7 Is there an easier way instead of digitizing? Is there a way to access and export the dataset from this webmap instead of digitizing? thanks, -j


Add a OpenStreetMaps road data layer on Google Maps API v3

I've successfully embedded a Google Maps API v3 map on my site, and added OpenStreetMaps (OSM) as the base layer. Is it possible to use Google's tiles as the base layer, and then overlay _street data_ on top of that, rather than the entire OSM layer?               <script type="text/javascript" charset="utf-8">         $(function(){                  google.maps.MapTypeId.push("OSM");                  var osmMapType = new google.maps.ImageMapType({                 getTileUrl: function(coord, zoom) {                     return "http://tile.openstreetmap.org/" + zoom + "/" + coord.x + "/" + coord.y + ".png";                 },                 tileSize: new google.maps.Size(256, 256),                 isPng: true,                 alt: "OpenStreetMap layer",                 name: "OpenStreetMap"             });                  $("#map-canvas").gmap3({                 action: 'init',                     options:{                       center: [22.49156846196823, 89.75802349999992],                       zoom:2,                       mapTypeId: google.maps.MapTypeId.ROADMAP,//google.maps.MapTypeId.OSM,                       mapTypeControl: true,                       mapTypeControlOptions: {                         style: google.maps.MapTypeControlStyle.DEFAULT                       },                       navigationControl: true,                       scrollwheel: true,                       streetViewControl: true                     }             });         });     </script>


Reduce extent of raster processing

I need to reduce the processing extent of my raster by some sort of automatic procedure in python. Currently the borders of the dataset I've been using has anomalies (due to the classification of the landsat scanning lines). ![ScanLines Classified](http://i.stack.imgur.com/97zi3.jpg) I would like to reduce the extent of the raster processing so that when I retrieve the information above by a SetNull function the borders are not processed. Something like ![Area](http://i.stack.imgur.com/bTjzW.jpg) Now this can easily be done with an extract by mask function (i.e. figure2) although that would require manual selection of the masking region. Is there a menthod I can utilize to calculate the extent region (ignoring nodata values) minus the border width in the x and y direction of each row/column and then utilize the `arcpy.env.extent` in Python?


Programs to log NMEA on Nomad 900?

I am interested in logging NMEA data on a Trimble Nomad 900. Unfortunately the SatViewer records a maximum of 64,000 characters, too little for my application. I don't want to have to purchase Terrasync/Pathfinder because I already have tools to parse NMEA and those pieces of software are overkill for my intended use. Thanks in advance.


Has anyone seen this error in modelbuilder?

I have a model that runs thru a rather lengthy set of processes. The screenshot below shows an error(bottom left of screenshot) that pops up several times while the model is running. It seems to popup between tools in the model but I can’t say for sure. If I click OK on the error each time it pops up, the model continues running successfully. As you can see, it is pretty non-descript so I’m just wondering what’s causing it. ![Screenshot](http://i.stack.imgur.com/t2ejz.jpg) Edit: PolyGeo, now that you mention it, I was initially using file gdb tables and feature classes as inputs to the model. IIRC, I started seeing this issue once I switched the inputs to an Oracle SDE gdb and a SQL Server View. Its an off and on thing where sometimes it pops up and sometimes it doesn't pop up at all. I had an issue once when running this model where the DB admin had to kill the connection...he didn't say why. I'll have to investigate further but I think its definitely related to the DBMS. I'm using Desktop 10.2 and the parent model has seven submodels. All models contains standard geoprocessing tools and no stand alone python scripts.


Buffer Selected Line Segments Using ArcGIS 10 or Python

I haven't quite found a solution for buffering selected line segments of a road network without extporting the selected line segments to a new file and then buffering. I want to avoid creating a bunch of new datasets for this task. Either environment will do (ArcGIS application or Python). Any suggestions?


How to create a field with a specified slope projecting from a 3D polygon

Im trying to figure out how to create a elevation raster from a known location (3D polygon) that projects at a specified slope (1%) for a specified distance(1-2 km). Does anyone have any idea how to do this?


Preparations for a transition to ArcSDE

So I'm a GIS Intern for a city government, we are preparing to purchase ArcSDE and implement it to our current GIS setup. I'm curious as to what should be done before it is purchased. Our GIS consists of several file geodatabases for parcels, road centerlines, hydrology, etc. as well as several other purposes. As well as shapefiles for the previously listed. So the question is, for experienced ArcSDE users out there, is there any sort of preprocessing/preparation work I should do before trying to implement ArcSDE? Any tips would be greatly appreciated. Yeah I've tried looking around for tutorials and things and can't find any useful sources :(. I'm asking about making changes to feature classes or to the file geodatabases themselves. I wanna say we are getting Enterprise and plan on doing Check-In/Check-Out replication We are currently using ArcGIS Desktop 10.2.1


Security alert dialogbox when launching arctoolbox

Every time I launch arctoolbox from arccatalog, the following dialog box appears: ![enter image description here](http://i.stack.imgur.com/veLTW.png) Then, whatever I answer, it takes ages (~36s) before the box disappear and I can finally use arctoolbox. What could be done to remove that?


can we update imagemosaic layer with geoserver rest api?

I am new to GIS and geoserver. I am trying to add and modify imagemosaic layers to geoserver with geoserver rest api. Actually I could add the layers. But I have modify them with rest api while adding or after adding. The parameters I have to set is: Layer name, srs handling, declared rss, background values, use jai image read Can I set these values with rest api? If there is documentation about it i will read but I couldn't find any. And if we can set these values I have to solve out how? Thanks,,


Write a cpp-Plugin for QGis

I want to write a plugin for QGIS 2.0 Dufour in C++. Do I have to take the sources/API from Version 1.8 (Lisboa)? It is right, that the actual API only can be used for the actual Releases(Weekly Master)? In this case I have to download the actual sources and recompile it everytime I want to deploy it? Thanks for any suggestions.


How to load memory output from QGIS processing?

Many processing algorithms have an option to save output as memory. If ran from toolbox, this works fine, because you can tick the "Open output file after running algorithm" box. Looking at the source code for QGIS, the code for loading the memory layer seems to be defined in the function `Postprocessing.handleAlgorithmResults`. The function uses access to an `alg.outputs` list, and does the following:               for out in alg.outputs:         progress.setPercentage(100 * i / float(len(alg.outputs)))         if out.hidden or not out.open:             continue         if isinstance(out, (OutputRaster, OutputVector, OutputTable)):             try:                 if out.value.startswith('memory:'):                     layer = out.memoryLayer                                # NOTE!!                     QgsMapLayerRegistry.instance().addMapLayers([layer])   # NOTE!!                 else:                     # ...       When you run processing algorithms from the console, is there a way to load the layer without access to this object? I can run               processing.runalg("qgis:intersection", layer1, layer2, "memory:")      or even               processing.runalg("qgis:intersection", layer1, layer2, "memory:myLayerName")      I can however not find a way to grab the resulting output.


How to connect contour lines from two different sources?

How to connect contour lines from two different sources (LiDAR & Photogarametry) in order to obtain continuous and uniform contour lines layer ? how can i interpolate them together ? How can i Snap lines only if they have the same Height value ? ![contours from 2 sources](http://i.stack.imgur.com/lXrQl.jpg)


Labeling only the latest feature

I need to label only the latest feature in a particular feature class. I actually would like to automate it, so that when I will add any new feature to this feature class, the label will automatically be assigned to the latest one. Could anyone please help me with line of code?


What is the best way to apply the correct ZIP Code when Geocoding

I need to batch geocode a set of addresses, some of them are missing the zip code causing the lat/lng to be placed incorrectly, is there a good way easily correct this? I am using a few geocoders, gps visualizer, esri's us streeets geocoder, and OSM geocoder in qgis.


How to add a legend and toggle the legend with a layer in Openlayers?

I'm currently trying to add a legend to my map and to show/hide it, when the corresponding layer (an image layer) is shown/hidden.   My approach was to add a custom button with a .png of my legend. That worked without problems. The image layer ("light") is deactivated when loading the map:               light.setVisibility(false);      Now, I want to show the legend when activating the layer and hide it when deactivating the layer and so on. And that works as well, but only once using these lines:               (function legend() {         light.events.register("visibilitychanged", light, function() {         if (light.getVisibility() === true) {             map.addControl(button);             } else {             map.removeControl(button);             }         });       }());      I guess i have to reload the function somehow after the button is removed? All things I tried didn't work and that was the best I was able to come up with. Any help is much appreciated.


Why is it not advisable to perform spatial analysis on layers with source datasets having same datum but different projections?

I have read that it is not advisable to perform spatial analysis on two layers with source datasets having the same datum but different projections loaded in ArcMap. Since the second layer should supposedly line up with the first one due to lining up data on the fly in ArcMap, I was wondering what could possibly go wrong?


How to set drawing style for QgsRasterLayer in QGIS 2?

I am using PyQGIS.   In QGIS1, I can use:               rlayer.setDrawingStyle ( QgsRasterLayer.PalettedColor )        but In QGIS2, I got some errors. I follow pyqgis developer cookbook, put some demo scripts in qgis python console without success.               rlayer.drawingStyle() #got AttributeError     rlayer.setDrawingStyle(QgsRasterLayer.SingleBandPseudoColor)      #got AttributeError from QgsRasterLayer, change it to QgsRaster     rlayer.setDrawingStyle(QgsRaster.SingleBandPseudoColor) #got TypeError


How to add a search box (OpenGeo Suite Client SDK)?

I'm using the OpenGeo Suite Client SDK to create a web mapping application. To perform searches I'm currently using the QueryForm plugin. Instead of using the QueryForm I would like to add a search box to "map.tbar" - the results should ideally still be listed in the east panel's FeatureGrid. The layer being queried is a polygon shape file, containing buildings and real properties. The layer is for querying only, its only fields are FID, shape and ID (adresses and real properties). What would be the way to do this? I've been looking at the GeocoderComboBox and AutoCompleteComboBox, but I'm not quite sure how to implement them the way I wish to do (still new to programming). Any help is appreciated. Thanks, Markus.


How to implement signals on add or modify feature?

I've created a test plug-in in QGIS successfully, but have had a hard time implementing a signal on a vector layer event (add or modify feature). I've add this to my initGUI method: QObject.connect(self.iface.activeLayer(),SIGNAL("featureAdded(QgsFeatureId)"),self.addedGeometry) and this definition:               def addedGeometry(self, intValue):             QMessageBox.information( self.iface.mainWindow(),"Info", "Hi!!!" )      I get the following error when I load my plug-in:               QObject.connect(self.iface.activeLayer(),SIGNAL("featureAdded(QgsFeatureId)"),self.addedGeometry)     TypeError: arguments did not match any overloaded call:       QObject.connect(QObject, SIGNAL(), QObject, SLOT(), Qt.ConnectionType=Qt.AutoConnection): argument 1 has unexpected type 'NoneType'       QObject.connect(QObject, SIGNAL(), callable, Qt.ConnectionType=Qt.AutoConnection): argument 1 has unexpected type 'NoneType'       QObject.connect(QObject, SIGNAL(), SLOT(), Qt.ConnectionType=Qt.AutoConnection): first argument of unbound method must have type 'QObject'      If anyone could help I would be very grateful! I can't find a lot of examples of implementing signals on QGIS objects.


Using arcpy from the command line on ArcGIS Linux Server 10.1

We installed ArcGIS SERVER (not Desktop!) 10.1 on a RHEL (Red Hat) machine. ArcGIS Sever is supported for Linux. However, it appears they don't make it easy to script in python, from the command line, making use of the arcpy module. There are many things awkward about this. The arcpy library that they use is actually for Windows, with some external patching. In fact, it comes with its own implementation of python. So pathnames look a little weird. You have to be the arcgis user to even get access to python or arcpy. So here is how my code looks:               import arcpy     """                                                                                From http://resources.arcgis.com/en/help/main/10.1/index.html#//0155000005mm000000          The Python libraries installed by an ArcGIS Linux Server installation is Windows 64-bit Python. Therefor, stand-alone Python scripts should always use the Win\     dows path separator (\) when accessing data and other Python modules. Also, be aware that there is a Z:\ which is mapped to the Linux systems' /. Within a Pyt\     hon script you can access data and other Python modules using Z:\ mapping. For example, data in /myuser/myproject/data can be accessed as z:\myuser\myproject\\     data.                                                                                This does not apply to Python scripts you publish from ArcGIS for Desktop because the publishing process converts all paths to the correct format.                       """     """                                                                                I'm using:                                                                                sudo su arcgis                                                                                /opt/arcgis/server/tools/python ./arcpytest.py                                                                                """               dataDir=r"z:home\arcpy-tmp"+"\\"     print dataDir          shp='tgr39105lkA.shp'     arcpy.FeatureToLine_management(dataDir+shp, 'temp_links')      What I get back is:               $ /opt/arcgis/server/tools/python ./arcpytest-1.py     z:home\arcpy-tmp\     Traceback (most recent call last):       File "./arcpytest-1.py", line 28, in <module>         arcpy.FeatureToLine_management(dataDir+shp, 'temp_links')       File "C:\Program Files\ArcGIS\Server\arcpy\arcpy\management.py", line 2343, in FeatureToLine         raise e     arcgisscripting.ExecuteError: Failed to execute. Parameters are not valid.     ERROR 000732: Input Features: Dataset z:home\arcpy-tmp\tgr39105lkA.shp does not exist or is not supported     Failed to execute (FeatureToLine).      Has anyone tried to use arcpy in this direct way? It's a REALLY powerful server, so we put ArcGIS server on it, but want to use arcpy for programming... Any help or hints appreciated!


LASToMultiPoint_3d

Trying to use the arcgis python command **LAStoMultipoint_3d** as seen here. I'm running the code snippet against one of my LAS files. It creates a `.shp, sbx, shp.xml, shx and dbf` file. However, the files are empty. The code is as follows:               #import script modules     import arcpy     from arcpy import env          arcpy.CheckOutExtension("3D")     env.workspace = "D:\myFolder\"          arcpy.LASToMultipoint_3d("D:/myFolder/myFile.las",                              "D:/myFolder/Output" , 1.5)          print 'done'      The code executes and creates empty files. I have tried this with a half-dozen different files varying in size and complexity.


how to buffer inside a polygon in Javascript

I'd like to know if it's possible to buffer a polygon on the INSIDE only using OpenLayers and other Javascript tools (e.g. jsts). Here's what I found about the topic in general: Buffering INSIDE a polygon with ArcGIS (possible): ArcGIS: inside buffer with attributes? Buffering a polygon on BOTH sides with JSTS (possible): Creating a Polyline Buffer in OpenLayers But I need to buffer a polygon on just the inside. Either a "real" buffer with the output geometry, or a "fake" buffer with just a line thickness and style. However line thickness in OL goes on both sides of the line. End goal is to get something like this: ![enter image description here](http://i.stack.imgur.com/Zt3k2.png) Side-question - Does any know if it's possible to do the country outline (like above) via Google Maps API?! Have not been able to find any reference. So I'm trying to duplicate with OpenLayers and a geojson outline of the world.


Is it possible for me to perform a map of New Mexico where deaf person lives in a county from the US Census database?

I would like to get some feedback from those who have been doing or perform on US Census. Map. My question here is there a way for me to get information on Deaf person living in New Mexico. Would US Census be helpful? On the other hand, do I need to find more information from the state of New Mexico? I wanted to show a map of where each county where a deaf person lives? Is that possible or is it the word "deaf" to broad to be included in the US Census database? Thanks for your input, Robert


Open my GeoServer WMS in ArcMap

I set up GeoServer on one of my old laptops a few weeks ago with a connection to the internet over my home cable internet line (see: http://bit.ly/e9OUNP and http://bit.ly/hj1a6J) I'm slowly learning the ins and outs of GeoServer, GeoExplorer, and other aspects of serving maps on a website. I assumed with this setup that I (or anyone else) would be able to access my GeoServer layers in ArcMap through WMS. I installed the OpenGeoSuite Community Edition in the default directory. I have published layers with the proper settings and check boxes enabled. I can load the layers in a GeoExplorer map embedded in a web page (see: http://bit.ly/ecpGE2). I have tried to load these layers into ArcMap using the following request (and many variations of it): http://donmeltz.dyndns.org/geoserver/wms?SERVICE=wms&VERSION=1.1.1&REQUEST=GetCapabilities but I cannot get this to work. A couple of hours of searching has not turned up a solution to my problem yet. I'm hoping someone on here can point me in the right direction, and provide me with some ideas of what to look at next. Thanks.


Where is qgis getting the NO DATA value from?

I have a project that contains an 8-bit paletted layer (call it the mask layer) and a float 32 single band layer (call this the clip layer). I create a new raster by using the raster calculator, defining the extent to match the clip layer, and the value is defined as (mask@1 >= 0) * clip@1 The resulting raster, when viewed with gdalinfo, has a NoData value of -3.4028234... e+038. Yet the input clip layer has no NoData value. Where is QGIS getting this NoData value from?


How to add WMS Layer from openweathermap?

I've tried to add some WMS Layers to my webmap using ESRIs Java Script API. I'm very interested in using weather WMS layers and I've already found some usefull. If there is anything better free available WMS service (for Africa or World), please give me a note ;)) One working example: http://developers.arcgis.com/en/javascript/samples/layers_wms/ Unfortunately it won't work on jsfiddle (I think because of proxy configuration) To add a custom one, I've tried to use this WMS Layer:               var wmsLayer = new esri.layers.WMSLayer("http://wms.openweathermap.org/service?");          //set visible layers - this is a world temperature layer         wmsLayer.setVisibleLayers([6]);           wmsLayer.setImageFormat("png");         map.addLayer(wmsLayer);       }      This also works. But how can I add other layers like these: http://openweathermap.org/tutorial/openlayers and combine them in webmap using ESRIs JavaScript API? Would be awesome to get some tipps ;) Thanks!


How to split multiple lines into equal parts or overcome split tool limitation using Basic or Standard license?

I would like to split multiple lines into equal **aproximately** 10m lengths. (Line length / 10m = number of segments) The split tool can only split one selected feature at a time. To overcome this I can merge multiple lines into a single multigeometry feature then run the split tool on the merged line. Works fine until I run into what seems to be a limitation - the **split tool can only split into a maximum of 999 equal parts**. When I enter a higher number the OK button goes dim. Is there a workaround to overcome this limitation or another way to split mupliple lines into equal parts? Basic or Standard licence only. Figure 1 - Split Tool dialog showing 999 equal parts: ![enter image description here](http://i.stack.imgur.com/gtlJE.png) Figure 2 - Split Tool dialog showing 1000 equal parts: ![enter image description here](http://i.stack.imgur.com/U8tCz.png)


Problem with reading data from PostGIS as KML file using PHP

I am using the following PHP code in order to read data from PostGIS and save it in KML format as a file. When I load the page a document with ".PHP" extension is returned (not .kml) which inside has the data in KML format. The page returns no error, but the console of the developer tools gives the following note: > Resource interpreted as Document but transferred with MIME type > application/vnd.google-earth.kml+xml: http://mydomain.com/.... Does anyone have any idea what the problem could be? The PHP code:               <?php          $dbconn = pg_connect("dbname='nyc' user='postgres' password='****' host='localhost'");     if (!$dbconn) {         echo "Not connected : " . pg_error();         exit;     }          $sql = 'SELECT *, ST_AsKML(ST_Transform(geom,26918)) AS kmlgm FROM nyc_census_blocks LIMIT 5';           # Try query or error     $query_result = pg_query($dbconn, $sql);     if (!$query_result) {         echo "An SQL error occured.\n";         exit;     }          $kml ="<kml><Document>";          for ($i = 0; $i < pg_numrows($query_result); $i++) {              $townname = pg_result($query_result, $i, 0);         $townkml =  pg_result($query_result, $i, 1);              $kml .= "<Placemark><name>Town: ".$townname."</name><description>".$townname."</description>".$townkml."</Placemark>\n";     }          $kml .= "</Document></kml>";          header('Content-Type: application/vnd.google-earth.kml+xml');          echo $kml;          pg_close($dbconn);     ?>      and the kml data returned inside a file with ".php" suffix!!:               <kml><Document><Placemark><name>Town: 1</name><description>1</description>360850009001000</Placemark>     <Placemark><name>Town: 2</name><description>2</description>360850020011000</Placemark>     <Placemark><name>Town: 3</name><description>3</description>360850040001000</Placemark>     <Placemark><name>Town: 4</name><description>4</description>360850074001000</Placemark>     <Placemark><name>Town: 5</name><description>5</description>360850096011000</Placemark>     </Document></kml>


problems using v.net.allpairs

I am using GRASS 6.4.3, and I am trying `v.net.allpairs` to do network analysis. I have a line shapefile and a point shapefile, and after I use v.net.allpairs, I type the command `d.vect network`, and I got this： ![enter image description here](http://i.stack.imgur.com/nnzRY.png) The cross in the photo is one of the point in the shapefile, but I do not want the **straight line** that links to the left curve which is my line shapefile because the **connect** operation in the algorithm may add the straight line first and then calculate the cost through that additional line, thus may have inaccuracy(error). All I want to ask is：Can I avoid this kind of result？ And how？ Because I have seen the manual and did not see other operations that I want.


How to display label in symbol in ArcGIS?

As you see on the map below, there are labels in circle coloured symbols. How can I show them in this format in ArcGIS? ![example map](http://i.stack.imgur.com/c0na2.png) If you look at the legend of this map it shows the different categories with different size and color, so my request is a little bit complicated. And there are 3 different labels displaying on the map. How can I do this kind of map? ![legend of map](http://i.stack.imgur.com/QW0Gx.png)


problems using GRASS v.net and v.net.allpairs

I am using GRASS 6.4.3, and I am trying to use it to do network analysis. I have a line shapefile and a point shapefile, and here is my code(I take v.net.allpairs as my reference)：               v.net tc_bf25 points=tc_sta_point op=connect thresh=1000 out=network nlayer=1     v.net.allpairs network out=allpairs     v.db.select allpairs      Here's some of my output： ![enter image description here](http://i.stack.imgur.com/2H535.png) and I compare this to the results I got previously using **pgRouting** and **roadgraph** , except the influence of `connect` operation, the result of GRASS is very different than other two, and I think GRASS is giving me wrong answers. But I have checked the code above and run the script more than two times, the GRASS results are all the same. Did I make mistakes that I didn't recognize？ Or in my case, pgRouting is a better tool to do the network analysis？


Is it possible to scale an ArcGIS feature in a multipatch feature class using Arcpy?

I've imported a number of features in Python using the following method: - > Import3DFiles_3d Unfortunately they are all extremely large so I would like to scale them using arcpy. I have seen that it is possible to scale from the ArcMap interface using the Editor Toolbar however I'd really like to be able to do this using Python. I've currently started looking at the following link but I'm not getting anywhere fast: - http://resources.arcgis.com/en/help/main/10.1/index.html#//018w00000005000000 If anyone could advise on scaling imported features, I'd be most grateful


Any .NET library to calculate polygon area?

Is there a .NET library to calculate area enclosed by polygon made up of geo coordinates (lat longs)?


How to create Vector Polygons at the same amazing speeds GISCloud is able to render them?

**EDIT (outlining new bounty added):** I have added the largest bounty possible to this question in hope of getting a very detailed answer. Your solution should outline a step-by-step guide utilizing open source tools (it believe giscloud.com is using some custom code combined with tools like tilestach, leaflet possibly postgis v2.0,++??...) and any code/customization to get the same rendering performance that giscloud.com is getting without using flash and while being visible on the same browser versions as giscloud achieved. Please note: the existing answers below provide some good info but point to an answer which one can implement. Bonus points for anyone who can outline how we can get QGIS to render layers using the same/similar method :) * * * I have been looking for a solid solution which would allow me to **create a web map and overlay vector polygons without taking forever to load such data** with the goal of allowing me to make each polygon display a different color on a hover event. As far as I am aware there are 3 specific options to achieve this through either canvas, SVG, Flash. Flash seems like it would be the best solution if it would work on apple iphones/ipads as its seems to provide the fastest rendering and cleanest display. Canvas seems to be the second best choice but takes VERY long if you have hundreds of polygons being displayed on a map whereas SVG takes even longer to render. I _almost lost hope in finding a solution to this problem_ but today I came across a company called GISCloud http://www.giscloud.com (currently in beta with free signup). This company has SOMEHOW managed to figure out an amazing way to render hundreds of vectors on a map in near real-time. I was amazed with their approach and my question to the community relates to how we can replicate their approach for use with existing technologies such a leaflet, openlayers, wax... Take a look for yourself by viewing this amazing demo: http://www.giscloud.com/map/284/africa Make sure you hover over any of the polygons on the page and test the zoom controls to see that these polygons are indeed vectors. What I have noticed by looking at requests with firebug is that the map is requesting specific json files. It seems that depending on the zoom level/area there are multiple json files being requested. EDIT: I should also mention here that once giscloud loads the data on the page hovering over a vector immediately changes the color without creating a new request. _EXAMPLES:_   * http://cft1.giscloud.com/t/1316509973/map284/layer1156/3/3/3.json   * http://cft1.giscloud.com/t/1316509973/map284/layer1156/3/5/3.json   * http://cft1.giscloud.com/t/1316509973/map284/layer1156/3/4/4.json   * http://cft1.giscloud.com/t/1316509973/map284/layer1156/3/3/4.json   * http://cft1.giscloud.com/t/1316509973/map284/layer1156/3/5/4.json I am assuming the url structure follows the standard tiling service logic (for example the 3rd to last folder being the zoom level...). In any case I have analysed the actual data of these json files and it seems the logic they are using follows some type of logic by which they create their vectors just based off these data values:   * **width/heith:** they define the width and height of the data being served in each json request   * **pixels:** here they define pixels values which I am assuming somehow relates to some general x/y pixel coordinates for generalized point levels? I am guessing they somehow have a way of automatically simplifying the region depending on the zoom level. I am assuming by them using pixels coordinates I am guessing they are dramatically reducing the size of the data that needs to be loaded compared to lat/long data.   * **styles:** here they define two RGB css values. "F" representing the polygon file color and "S" representing the polygon border color.   * **geom:** here is where I am guessing they are somehow defining specifically defining each polygon within the tile being loaded where such data is being defined based off the map container window. Whats also interesting is that each entry has a "S" value which I am assuming is used as an optional attribute or feature link value and at the end of each entry here there is an area which seems to define a specific per vector ID along with the layer ID which I am guessing is utilized to somehow join the data from each json tile request being called. EDIT: I am also assuming they somehow have figured out a way to automatically determine and split up the data which needs to be loaded for each tile depending upon the size of the data which would need to be loaded for the requested tile. Here is an extracted breakdown of one of these requests:               {"width":256,"height":256,"tile":     {"pixels":     [0,6461,-1,0,5,148,0,509,-1,10715,-1,1,-1,251,-1,1,-1,1,-1,251,-2,3,-1,255,-1,249,-2,5,-2,247,-1,509,-3,251,-1,2,-2,253,-2,252,-2,254,-1,255,-1,254,-1,255,-1,1276,-2,13,-1,233,-1,2,-1,253,-1,1,-1,255,-1,247,-1,1306,-1,1533,-1,1269,-1,1276,-1,2303,-1]},          "styles":     [{"f":"rgb(99,230,101)","s":"rgb(5,148,0)","lw":"0"}],          "geom":     [     {"s":0,"p":[4,143,5,144,3,146,1,146,2,143,4,143],"c":"layer1156_5098"},     {"s":0,"p":[-2,143,0,140,2,141,2,144,1,146,-2,144,-2,143],"c":"layer1156_5067"},     {"s":0,"p":[7,143,5,144,4,143,2,143,2,141,5,138,6,139,5,141,7,143],"c":"layer1156_5051"},     {"s":0,"p":[10,141,11,137,12,137,14,137,12,142,9,143,9,142,10,141],"c":"layer1156_5041"},     {"s":0,"p":[1,136,0,140,-2,143,-2,136,1,136],"c":"layer1156_5038"},     {"s":0,"p":[8,143,5,141,5,137,8,136,10,137,10,141,8,143],"c":"layer1156_5033"},     {"s":0,"p":[5,137,2,141,0,140,1,136,1,136,2,135,3,136,5,137],"c":"layer1156_5028"},     {"s":0,"p":[10,134,12,136,11,138,8,135,10,134],"c":"layer1156_5020"},     {"s":0,"p":[-2,133,0,136,-2,136,-2,133],"c":"layer1156_5005"},     {...}     ...     ]     }      I would really appreciate any input on how we can replicate the same (or similar) type of speed using postgis (which I what they seem to be using as well).


ArcObjects .NET: Get more info out of COM Exceptions - how to get description info from ErrorCode?

Ok, I can trap the COM exception and convert the ErrorCode to the proper enum, (e.g. esriGeometryError.E_GEOMETRY_WRONGTYPE). But is there a way to get the error description from the enum? For example, when I mouse over esriGeometryError.E_GEOMETRY_WRONGTYPE in Visual Studio I get a popup that says "The geometry parameter was of the wrong type for the method." It would be great to retrieve that information programatticaly once the exception is trapped. Any help greatly appreciated!


GeoServer: Best way to publish 2500 TIFF or 71 ECW files?

I have ortophotos of an area that I need to serve as a background map for my vectors. I got it as 2500 files (71,5MB each) in raw TIFF format with corresponding TFW word file - 180GB of data. Coordinate system is local and it matches with my vectors (there is no EPSG code for it but I named it as “32805” and put a proper definition). For desktop use in MapInfo I converted them to ECW (with some tool that comes with MapInfo) and retile as much bigger to have only 71 file because opening 2500 tiff files is overkill. I just merged 49 TIFF’s (7x7) into one ECW - 35000x35000pixels – largest is around 200MB) It works great and its really fast in MapInfo. Now I am quite confused - HOW to serve them in GeoServer ? I have publish one TIFF and one ECW for comparison. ECW is much faster in browser preview (I am aware of ECW server licensing issue but that shouldn’t be a problem). I have found a presentation “GeoServer on steroids” and read about ImageMosaic, ImagePyramid, retiling, adding overview etc it was quite informative but still don’t know what to do. My question is : how should I do it ? Mosaic or Pyramid, and if is one of the answers positive, I need your advice or some hints. I would really like to be ECW because of disk space so there will be no need to keep 180GB of tiff’s on a server. Data will be served through LAN with maximum 20 users connected at peak hour.SQL server have not so huge amount of data. Sorry if I miss some other info, but I will send it if needed. * * * > Geoserver 2.1.4, Windows 7 32bit, 2GB System Memory, (1.7.0_09 (Java > HotSpot(TM) Server VM), Native JAI+Native JAI ImageIO = true * * *               Original TIFF     gdalinfo D:\75720-47970.tif     Driver: GTiff/GeoTIFF     Files: D:\75720-47970.tif            D:\75720-47970.tfw     Size is 5000, 5000     Coordinate System is `'     Origin = (7572000.000000000000000,4797500.000000000000000)     Pixel Size = (0.100000000000000,-0.100000000000000)     Metadata:       TIFFTAG_SOFTWARE=Adobe Photoshop 7.0       TIFFTAG_DATETIME=2006:10:09 13:02:57       TIFFTAG_XRESOLUTION=72       TIFFTAG_YRESOLUTION=72       TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  ( 7572000.000, 4797500.000)     Lower Left  ( 7572000.000, 4797000.000)     Upper Right ( 7572500.000, 4797500.000)     Lower Right ( 7572500.000, 4797000.000)     Center      ( 7572250.000, 4797250.000)     Band 1 Block=5000x1 Type=Byte, ColorInterp=Red     Band 2 Block=5000x1 Type=Byte, ColorInterp=Green     Band 3 Block=5000x1 Type=Byte, ColorInterp=Blue     -------------          ECW file which is retiled from 7x7  original tiffs          gdalinfo D:\OF-45.ecw     Driver: ECW/ERDAS Compressed Wavelets (SDK 3.x)     Files: D:\OF-45.ecw     Size is 35000, 35000     Coordinate System is:     LOCAL_CS["LOCAL - (unsupported)",         UNIT["Meter",1]]     Origin = (7571500.000000000000000,4798500.000000000000000)     Pixel Size = (0.100000000000000,-0.100000000000000)     Corner Coordinates:     Upper Left  ( 7571500.000, 4798500.000)     Lower Left  ( 7571500.000, 4795000.000)     Upper Right ( 7575000.000, 4798500.000)     Lower Right ( 7575000.000, 4795000.000)     Center      ( 7573250.000, 4796750.000)     Band 1 Block=35000x1 Type=Byte, ColorInterp=Red            Overviews: 17500x17500, 8750x8750, 4375x4375, 2187x2187, 1093x1093, 546x546, 273x273, 136x136     Band 2 Block=35000x1 Type=Byte, ColorInterp=Green       Overviews: 17500x17500, 8750x8750, 4375x4375, 2187x2187, 1093x1093, 546x546, 273x273, 136x136     Band 3 Block=35000x1 Type=Byte, ColorInterp=Blue       Overviews: 17500x17500, 8750x8750, 4375x4375, 2187x2187, 1093x1093, 546x546, 273x273, 136x136


Process every 100 records in list

An arcpy merge process is failing due to the size of the datasets and I need to just merge every 100 tiles before it goes on to the next 100 tiles. I wrote the following snippet but it is just putting the same tif into the list. It's probably something simple that I am doing wrong but would appreciate your advice....               arcpy.env.workspace=root+"\\Stage1\\"+name+"\\GeoTif"     tifList=arcpy.ListRasters('*.tif')     tifListCount=len(tifList)          tempTifList=[]     x=0     oldx=0     while x <= tifListCount:         for tif in tifList:             while x<(oldx+100):                 tempTifList.append(tif)                 x=x+1                  print "List "+str(oldx)+" is"+ str(tempTifList)             # Add merge process on tempTifList             tempTifList=[]             oldx=x      Gives > , u'PFTI_Roads2_Z31.tif', u'PFTI_Roads2_Z31.tif', u'PFTI_Roads2_Z31.tif', > u'PFTI_Roads2_Z31.tif', u'PFTI_Roads2_Z31.tif', u'PFTI_Roads2_Z31.tif', > u'PFTI_Roads2_Z31.tif', u'PFTI_Roads2_Z31.tif', u'PFTI_Roads2_Z31.tif', > u'PFTI_Roads2_Z31.tif'] List 75900 is[u'PFTI_Roads2_Z32.tif', > u'PFTI_Roads2_Z32.tif', u'PFTI_Roads2_Z32.tif', u'PFTI_Roads2_Z32.tif', > u'PFTI_Roads2_Z32.tif', u'PFTI_Roads2_Z32.tif', u'PFTI_Roads2_Z32.tif', > u'PFTI_Roads2_Z32.tif', u'PFTI_Roads2_Z32.tif', u'PFTI_Roads2_Z32.tif


Wanting to build Web Mapping Application

I am wanting to build a web mapping application and am trying to gather all the pieces that I'll need. Is it possible to build REST services with opensource Geoserver from an ArcMap document and tools built in ESRI like I would using Arc Server Manager? Basically what I would like to do is build a website with a Flex Mapping application with geoprocessing widgets from my Arc Map document, build the rest services for it with Geoserver, and host it on something like Amazon Ec2. Is this possible? -Thanks Adam


Getting an" 001270 error : Consolidating the data failed." while publishing as map service in arcgis server 10.1

I am getting an error "001270 error" . I am inserting table for sql-server database and joining few records to polygon layer also I have lat,Lon data in few tables so I am displaying data using "Display x,y values" and publishing all data into arcgis server 10.1 (using ArcGIs Desktop 10.1) While publishing I was getting above error.. also my database & folders are registered in ArcGIS Server (through ArcCatlog I did this) I found following article which explain that its bug Link : http://support.esri.com/en/bugs/nimbus/TklNMDgzMjc5 ..but I am not sure.. So any help regarding this will be great..thanks in advance : )


Arcobjects in Python: assign interface to objects with feature (shape) geometry

I am trying to use ArcObject in Python. I do have problem with assigning an interface (i.e. IZ.InterpolateZsBetween) to my object with geometry data, code is:               from Comtypes.client import CreateObject, GetModule     import arcgisscripting          gp = arcgisscripting.create(9.3)     gp.setproduct("ArcInfo")          esriSystem = GetModule("C:/Program Files (x86)/ArcGIS/com/esriSystem.olb")     esriGeometry = GetModule("C:/Program Files (x86)/ArcGIS/com/esriGeometry.olb")     esriDataSourcesGDB = GetModule("C:/Program Files (x86)/ArcGIS/com/esriDataSourcesGDB.olb")     esriGeoDatabase = GetModule("C:/Program Files (x86)/ArcGIS/com/esriGeoDatabase.olb")          sPath ="c:/temp/test.gdb"     pWSF = CreateObject(esriDataSourcesGDB.FileGDBWorkspaceFactory,\                         interface=esriGeoDatabase.IWorkspaceFactory)     pWS = pWSF.OpenFromFile(sPath, 0)     pFWS = pWS.QueryInterface(esriGeoDatabase.IFeatureWorkspace)     pFClass = pFWS.OpenFeatureClass("test_z")     pFCursor = pFClass.Search(None, True)     pFeat = pFCursor.NextFeature()     pIZ = pFeat.QueryInferface(esriGeometry.IZ2)      Why this is not working? The output is:                Traceback (most recent call last):       File "<pyshell#175>", line 1, in <module>         pIZ = pFeat.QueryInterface(esriGeometry.IZ2)        File "C:\Python25\lib\site-packages\comtypes\__init__.py", line 1078, in QueryInterface      self.__com_QueryInterface(byref(iid), byref(p))      COMError: (-2147467262, 'Taki interfejs nie jest obs\xb3ugiwany.', (None, None, None, 0,       None))      So I tried:               pShape = pFeat.ShapeCopy     pIZ = pShape.QueryInterface(esriGeometry.IZ2)     pShape.InterpolateZsBetween(0, 0, 0, 4)      And this is working. But how to stick it back and save it in my feature? I am learning Python and ArcObjects, so please be understandful. Any help will be appreciated.


How could I turn my data into a heatmap / intensity map?

I have database that tracks a large number of entities. Each entity has a 'coverage area' or 'catchment area' - a radius in miles from a specified UK postcode, eg               Leyton office     4.3m    E11 4LL     Stratford office  5.2m    E15 6ZZ     East Ham office   2.1m    E4 8QQ      I want to turn this into a geographic intensity map - a map of a region which is shaded or coloured to indicate the number of entities that cover each point on the map. Segmentation into small regions is also acceptable. This is probably a one-off project, so it's unlikely my organization will provide money to buy software. That means zero or at least no-cost solutions are preferable. Thanks in advance,


Fastest way to spatially join a point CSV with a polygon Shapefile

I have a 1 billion point CSV file and a shapefile with around 5,000 polygons. What would be the fastest way to spatially join points and polygons? For each point, I need to get the containing polygon id. (Polygons don't overlap.) Usually, I'd load both data sets into PostGIS. Is there a faster way to get the work done? I'm looking for an open-source solution.


Selecting the MAX value from a feature using ArcObjects?

I have a feature (point) with a field called 'Sequence'. The values in this column range from 1 to _n_. I can use IQueryFilter to get the first value, 1.                                   Dim pQry As IQueryFilter = New QueryFilter                         pQry.WhereClause = """Seq"" = 1"                         pCursor = pFC.Search(pQry, True)                         If TypeOf pCursor Is IFeatureCursor Then                             pFeat = pCursor.NextFeature                             If TypeOf pFeat Is IFeature Then                                 fStartElevation = pFeat.Value(nElevationFld)                             End If                         End If      I cannot figure out how to use this method to select the MAX value, since it varies across features. Is it possible? I'm thinking I just need the correct syntax. I have tried a few variations of `pQry.WhereClause = """MAX""Seq"""` after doing some forum reading regarding aggregate functions in ArcGIS, but I haven't found a solution that works. If this isn't possible (I figured it should operate similarly to SQL), is there another preferable method for performing this task? I am aware of Select By Attribute, but the syntax for this is tricky (examples online do not work in Arc10) and still requires the attribute table to be read. EDIT: Here is a working solution using @George's suggestion based on IDataStatistics. I may just turn this into a small function that takes in a Featureclass and field name for the sake of re-usability.                                   pCursor = pFC.Search(Nothing, False)                         If TypeOf pCursor Is IFeatureCursor Then                             Dim pData As IDataStatistics = New DataStatisticsClass                             pData.Field = "Sequence"                             pData.Cursor = pCursor                                  Dim pStatResults As IStatisticsResults = pData.Statistics                             Dim intMaxSequence As Integer = pStatResults.Maximum                         End If


JS/browser API library to show/interact with many markers with the best performance?

What's the best Javascript library that can show about 50000 (and more) clickable markers without resorting to things like clustering markers,etc? The markers are clickable and I want to associate them with some action. The main issue is coping with browser performance, so I'm already only considering performance in Chrome and Firefox. Yes, I also understand browser performance is dependent on the PC performance but still, there are better and worse libraries. I have experience with GMaps v3 but with so many markers it is very slow. Changing to GEarth plugin on the browser solves the problem. Currently I'm using OpenLayers serving from my GeoServer WMS. I guess I need a JS lib like Openlayers tailored for very high performance, for example using WebGL to make use of the PC's GPU. I know openlayers v3 would have webgl but I really want to hear the feedback from experienced users for these kind of big data on a map. Also it should be opensource and free to use and be able to connect to my WMS to serve the tiles, like openlayers does. I also don't mind if it's not JS API but some plugin I have to install on my PC / browser, like GEarth, if that improves the performance radically.


ArcGIS 10 and IGeometryServer2

Has anyone already use IGeometryServer2 or IGeometryServer within C# code ? I'm looking for a sample since I have not succeeded in using it. I have already selected:   One Polyline (IPolyline) : for cutting   One Polyline (IGeometry) : the shape to be cut Thanks in advance. Best regards,


Tilemill OSM Bright Installation problems

I followed the OSM Bright installation guide and patiently imported germany as osm data into tilemill (→ see here: https://www.mapbox.com/tilemill/docs/guides/osm-bright-mac- quickstart/#step_0_download__install_required_software) Now I run in some issues upon first launch: – I want to restrict the zoom level and center point to a smaller area. – When I open the project settings, I can change (decrease) the zoom level and center point, but when I click save, nothing happens. I f I do this repeatedly, I get a 'Forbidden' error message. Going back to map mode with cancel works, and maybe the settings are saved anyways… – Overall, the performance of the map is very slow and not much is rendered after 20min. Definitely not the performance needed to edit the map… Am I missing something here? Is this a matter of the .osm size or what? Anybody had the same problem (I have done the install numerous times now, without improvement) Cheers Lars


Cartodb.js tooltip on a viz

I am having a hard time implementing a UI tooltip in a CartoDB.js viz. Could someone demonstrate how to implement using **this fiddle** on the variable {{NAME}}?


Problem with result of query in PostGIS

I run the following query and irrelevant of the point coordinates and radius I always recieve all the records of the table I'm querying. Can anyone point out what the problem could be with this query?               SELECT * FROM Buildings WHERE ST_DWithin(geom, ST_Buffer(ST_GeomFromText('POINT(0 0)', 26918), 2000, 'quad_segs=8'),0);      Just to be clear, my aim is to select all buildings that fall into a selected circular buffer. Here the (0,0) are the coordinates of the center of circle and 2000 is the radius. No matter what value I put for these, the query returns all the records of table. Please help! The following image shows the results of 12 records: ![enter image description here](http://i.stack.imgur.com/5gfPl.png)


How can i find shrines nearby location

How can i find shrines/complex/institute nearby location in qgis and highlight them. The below image is similar to my shape files which contains all information of all attributes,i don't know python, i was using spatial query to find location within area or beside roads, now i want to find complex, malls near to my location and highlight them. ![Location image](http://i.stack.imgur.com/GdNil.png)


Arcgis custom Legend rest service

There is a description for Legend services on arcgis web site. > The legend supports the following layer types: ArcGISDynamicMapServiceLayer, > ArcGISTiledMapServiceLayer, FeatureLayer and KMLLayer. If the layer is an > ArcGISDynamicMapServiceLayer or ArcGISTiledMapServiceLayer created using > ArcGIS Server 10.0 SP1 or higher the legend is generated using the REST Map > Service legend resource. > > httpp://mapservice-url/legend > > If the layers are version 10 or lower the legend is created using the > ArcGIS.com legend service. In order to use the ArcGIS.com legend service > your map service needs to be publicly accessible and your application must > be able to access ArcGIS.com. My version is 10. So I want to a rest service that returns my legends. Is there a sample app about this. How can I return legend images. Should I use arcgis dll references?


QGis 2.0.1 mapserver extension missing

I tried to export my Qgis projects to mapserver, using RT mapserver extension, because there is not any other extension in my version of QGis.(QGis 2.0.1 Dufour on ubuntu 12.04 LTS). After installing the extension, when i try to export, there is the error:                 Une erreur est apparue lors de l'exécution du code python :      Traceback (most recent call last): File "/home/jannes/.qgis2/python/plugins/rt_mapserver_exporter/plugin.py", line 65, in run from .mapfileexportdlg import MapfileExportDlg File "/usr/lib/python2.7/dist-packages/qgis/utils.py", line 453, in _import mod = _builtin_import(name, globals, locals, fromlist, level) File "/home/jannes/.qgis2/python/plugins/rt_mapserver_exporter/mapfileexportdlg.py", line 30, in import mapscript File "/usr/lib/python2.7/dist- packages/qgis/utils.py", line 453, in _import mod = _builtin_import(name, globals, locals, fromlist, level) ImportError: No module named mapscript Version de Python : 2.7.3 (default, Sep 26 2013, 20:26:19) [GCC 4.6.3] Version de QGIS : 2.0.1-Dufour Dufour, exported Chemin vers Python : ['/usr/share/qgis/python/plugins/processing', '/usr/share/qgis/python', '/home/jannes/.qgis2/python', '/home/jannes/.qgis2/python/plugins', '/usr/share/qgis/python/plugins', '/usr/lib/python2.7', '/usr/lib/python2.7/plat-linux2', '/usr/lib/python2.7/lib-tk', '/usr/lib/python2.7/lib-old', '/usr/lib/python2.7/lib-dynload', '/usr/local/lib/python2.7/dist-packages', '/usr/lib/python2.7/dist-packages', '/usr/lib/python2.7/dist-packages/PIL', '/usr/lib/python2.7/dist-packages/gst-0.10', '/usr/lib/python2.7/dist- packages/gtk-2.0', '/usr/lib/pymodules/python2.7', '/usr/lib/python2.7/dist- packages/ubuntu-sso-client', '/usr/lib/python2.7/dist-packages/ubuntuone- client', '/usr/lib/python2.7/dist-packages/ubuntuone-control-panel', '/usr/lib/python2.7/dist-packages/ubuntuone-couch', '/usr/lib/python2.7/dist- packages/ubuntuone-installer', '/usr/lib/python2.7/dist-packages/ubuntuone- storage-protocol', '/usr/lib/python2.7/dist-packages/wx-2.8-gtk2-unicode', '/usr/share/qgis/python/plugins/fTools/tools', '/home/jannes/QGis workspace'] How can I solve this problem please?


Is it necessary to master openlayers to use geoserver?

Is it necessary to learn openlayers to use geoserver? Can i completely ignore OpenLayers so long as i master using geoserver to publish data on my website because i understand openlayers is contained in geoserver? I would like to set up a web GIS but i find learning openlayers or leaflet difficult. I would like to create a web application that receives data from, say a mobile phone and combines it with data served through GeoServer from my PostGIS database


Converting Excel Table to Geodatabase with Repeating Numbers Results in all Zeroes

I'm trying to convert an Excel file to a DBF file. I understand how to do it correctly using the Table to Table tool. However I am having an issue where a column of repeating numbers (labeled ZoneNum) are converted into nothing but zeroes. There doesn't seem to be any other issues with the other columns such as street names, ID number, work done, etc. I believe it has to do with the repeating numbers. For example 16 streets might be under ZoneNum 1 so their attribute would be 1. The next 7 streets might be under ZoneNum 3. I've tried to store the number as text instead but then nothing showed up. I've formatted each column as numbers or text. I've tried converting to a csv file instead to then convert to DBF file with no results. Please help! Here's 3 sample lines from my file (minus the linebreaks):               104,14,Yosemite Drive,Everglade Circle,South End,Cape Seal          105,14,Zion Avenue,Jefferson Street,Carlsbad Street,Slurry Seal          106,17,Aberdeen Drive,Venice Avenue,West End,"1.5"" Edge Grind & Overlay"


Issue with `-sds` naming convention in gdal_translate

I'm converting ~100 .hdf files to .tiff using gdal and getting some unfortunate naming conventions when using `-sds` `for %A in ("...\*.hdf") do gdal_translate -of GTiff -sds "%A" "%A.tiff"` However, it appears that it is first converting the files and then separating the subdatasets, therefore giving me file names like .tiff1, .tiff2, etc (completely useless). Any suggestions on how to fix this within the same piece of code or with a simple second step?


How do I redefine the projection of a DEM?

I have a DEM that was exported from AutoCAD civil 3D 2013 but it's in the wrong projection. I'm using an area specific projection that was not provided as an option in CAD DEM export dialog box. I choose an other projection to finish the export and now I want to associate the proper projection to the DEM. Is there a way to redefine the projection of the DEM in ArcGIS 10.1 to the proper one?


Van der Grinten IV projection. Is it possible and how?

I like the look of the Van der Grinten IV projection with it's two hemispheres and would like to use if for a map of the world. I could only find information on how to use the Van der Grinten I projection. Can I transform my data to the VDG-IV projection and use that for my project? If the Van der Grinten IV projection is impossible are there then other types of similar globular projections available?


How to get power lines and substations data for Mexico

Hello guys. I am looking for an electric lines shapefile for Mexico and other Central and South American countries. Do you have any idea where I can get this data?? Thanks, Gracias, Salamat, Asante Sana


Analysis of GPS error

I'm wondering if there have been any studies which attempt to analyze or characterize GPS error distribution over different situations/conditions (e.g. surrounded by tall buildings/mountains, precipitation/clouds/fog, rural areas, etc.) to understand how these errors may change as compared to 'perfect' conditions.


Cached map performance when viewed in ArcMap

ArcGIS server cached map service that displays at the following scales: 1:20,000 1:10,000 1:5,000 1:1,000 Refresh takes 1 second when using the Javascript/silverlight viewer, but 4-7 seconds when consuming the service in ArcGIS Desktop. Any thoughts to as why this might be? Thanks John


Extrude Inclined 3D Boreholes

I need to extrude the 3D Borewell data using the inclined angle in ArcScene. I've seen only the options for extruding them in one direction(i.e., perpendicular downwards or upwards). I got point data with Direction, distance and I need to use both the attributes for extruding below the surface.


Adding .bil raster file to GeoServer and viewing it on GeoServer or GeoNode

Ok so my Boss gave me this job. Uptil now i was working on ESRI products and this is my first venture into OpenGEO. He has given me a link, http://www.worldclim.org/current, from where i need to use the 10 arc minute Bioclim raster which are in .bil format and add them to GeoServer. I am absolutely clueless how to do that? Pls help me


Open Source Flexviewer application

I would like to know if there are any open source flex viewers some thing on the same line as ArcGIS Flex Viewer ??


Is ArcGIS Desktop-ArcMap 10 affected by the stacked toolbars bug (NIM003161) when opening an MXD on multi-monitor configurations?

In ArcMap 9.3, we had issues with our toolbars always stacking up instead of remaining where they were placed. This was mainly an issue on dual monitor workstations and we would occasionally see the issue on a single monitor workstation. When we reported the issue to ESRI, they said it was a known bug: NIM003161. I won't have an opportunity to install ArcGIS 10 on a dual monitor workstation for a couple of months and I was wondering what the user experience is when opening saved MXDs on dual (or triple) monitor workstations? Are toolbars staying put? I looked up bug #NIM003161 and its status is on-hold and last modified on March 25, 2010 and the version fixed is still blank. That part worries me a little, but I'm hoping ESRI just hasn't gotten around to updating that bug. Or maybe ESRI will read the answers to this question on http://gis.stackexchange.com and update their bug info accordingly!


hDC of DynamicDisplay always 0

When I try to access the hDC of my dynamic layer the result is always 0. Is there a way to get a non zero hDC from a layer that inherits BaseDynamicLayer?               public class BalloonCalloutLayer : BaseDynamicLayer     {          public override void DrawDynamicLayer(esriDynamicDrawPhase DynamicDrawPhase, IDisplay display, IDynamicDisplay dynamicDisplay)     {          int hdc = display.hDC; //This is always 0     }          }


ArcSde Oracle Maximum number of grids per feature (8000)

I have a problem with the import of data (lines) to ArcSde Oracle. 2011-07-14 12:10:04| 22.9| 0.0|ERROR |Error while executing stream insert operation for table 'CHEMN'.:SDE_ERROR_CODE(-51) Underlying DBMS error. SDE_EXT_ERROR(29875) SDE_ERROR_MSG1(ORA-29875: échec d'exécution de la routine ODCIINDEXINSERT ORA-20092: Maximum number of grids per feature (8000) exceeded. ORA-06512: à "SDE.ST_DOMAIN_METHODS", ligne 1982 Have you a solution ? Thanks


Changing DEM in Arc's on-the-fly orthoprocessing

ArcGIS has the ability to do on-the-fly orthorectification once you set it up in a mosaic dataset. You point to the DEM when you set up the dataset. That's cool, but let's say you want to change the underlying DEM that you're rectifying to once the dataset has been created. Is there a way in Arc to do this "on the fly" or do you have to start from the beginning?


distance to polygon boundary

i have a layer of points, and a vector layer of county boundaries for the UK. I would like to be able to find the nearest county boundary to any given point and the distance using QGIS. any ideas? I've installed processing, and trying to use those tools -- but no luck so far. I tried converting the polygon to nodes, but the nodes are spaced too widely for this to be accurate.


Draw order (zIndex) of OpenLayers polygons and points in the same layer

My question is similar but slightly different from this one: http://stackoverflow.com/questions/5003349/openlayers-vector-features-z- indexing I have an OpenLayers Vector layer which contains polygons and points from PostGIS. I do not want to create a separate table for each entity. I have a StyleMap set for the layer and have zIndexing set to true. My points are using an external graphic and my polygons just a normal style. Problem: I cannot get the external graphic to render above the the polygons. If I used a normal point then it does render above the polygons. I believe I've tried all combinations of graphicZIndex and zIndex for the point and polygon layers. Help? Layer def:                       var lyrJamoats = new OpenLayers.Layer.Vector("Jamoats", {                     strategies : [new OpenLayers.Strategy.Fixed()],                     styleMap : Styles.styleJamoat,                     rendererOptions : {                         zIndexing : true                     },                     protocol: {protocol stuff}         });      Style def:               styleJamoat : new OpenLayers.StyleMap({         "default" : new OpenLayers.Style(null, {             rules : [new OpenLayers.Rule({                 symbolizer : {                     "Point" : {                         fillColor : "#ffcc66",                         strokeColor : "#ff9933",                         strokeWidth : 2,                         graphicZIndex : 1000,                         //zIndex: 99999999,                         externalGraphic : "./img/jamoat_pin.png",                         pointRadius : 30                     },                     "Polygon" : {                         strokeWidth : 2,                         strokeOpacity : 1,                         strokeColor : "#FFFFFF",                         fillColor : "#002776",                         fillOpacity : 0.75,                         zIndex: 0                         //graphicZIndex: 0,                     }                 }             })]         }),      Cheers, James @geo_james


QGIS image export changes the map center and scale

I've been trying to export a few map layers using QGIS print composer. First, I loaded the layers and opened the composer using `Project -> New Print Composer`: ![enter image description here](http://i.stack.imgur.com/tDsZa.jpg) Gave it a proper name: ![enter image description here](http://i.stack.imgur.com/0p3UI.jpg) Added a new map: ![enter image description here](http://i.stack.imgur.com/tCJno.jpg) Exported as image: ![enter image description here](http://i.stack.imgur.com/nUYA5.jpg) However, the outcome is shifted a few dozens of kilometers westwards: ![enter image description here](http://i.stack.imgur.com/XMMLw.jpg) This specific west shift is inconsistent; sometimes I experience a smaller shift or a shift in another direction. Moreover, the scale often changes too - which makes map exporting impossble. I'm using the newest QGIS version: ![enter image description here](http://i.stack.imgur.com/WHN2u.jpg) **Any idea how to export the image properly? I'm currently using screenshots, which are a poor solution. Thanks a lot in advance!**


How to find which subtype using which domain?

I'm trying to delete domain, but i can't because it is in use at an attribute domain. How can I find out at which layer and which subtype ? I have a big geodatabase with more than 30 layers an each one have few subtypes ?


Setting border-style to 'no pen' using QgsGraduatedSymbolRendererV2

I am developing a QGIS plugin that, after a series of inputs, will produce a thematic map using the QgsGraduatedSymbolRendererV2. Because the polygon are very small, the boarders make it such that when zoomed out the entire map almost looks completely black. I want to remove the boarders to do this but no documentation or form I have found have demonstrated how to do this. Below is what I have to create the thematic map, but somewhere in there I need to set the border style to 'no pen'. But how!?                 uri = QgsDataSourceURI()       # set host name, port, database name, username and password       uri.setConnection("localhost", "5432", "public", "postgres", "password")            uri.setDataSource("susip", str(filename)+"_mce_results", "the_geom")       layer = QgsVectorLayer(uri.uri(), str(filename)+"_mce_results", "postgres")            fieldName = "score"       fieldIndex = layer.fieldNameIndex( columnName )       provider = layer.dataProvider()       numberOfClasses = 10       color1 = QColor("red")         color2 = QColor("green")       color3 = QColor("white")       ramp = QgsVectorGradientColorRampV2(color1, color2)            renderer = QgsGraduatedSymbolRendererV2.createRenderer(layer, columnName, numberOfClasses, QgsGraduatedSymbolRendererV2.Quantile, QgsSymbolV2.defaultSymbol(layer.geometryType()), ramp)            layer.setRendererV2( renderer )       QgsMapLayerRegistry.instance().addMapLayer( layer )


How to sum pixel values within a polygon zone in QGIS?

I have two layers: a raster layer and a vector layer (polygon). I can see the values of the raster pixels but I don't know how to combine them with information from the polygon layer. More specifically, I would like to :   1. sum values of all pixels within a the borders of a polygon   2. count the number of pixels with a specific value within the borders of a polygon Can anybody help and advise on how to do this in QGIS? Many thanks!


Embedding slippy map in Access form

I have an access database which holds several sites. This is linked to a MapInfo TAB file which holds geometry for each site. I would like for our non-GIS users to be able to view the sites in an interactive way (i.e. a slippy map with layers, not just a static image - I've already got static images in there) rather than having to boot up MapInfo or, worse, ask someone else to... I have embedded a MapInfo mapper inside a form but distributing this has proved to be more expensive than I originally expected (though to be fair my expectations may have been unreasonable!). As far as I know it is not possible to embed HTML in a form - otherwise I would use OpenLayers or similar. I could include a link to a web-based map but that isn't ideal... Any suggestions? What software is available for this sort of thing?


Is the term "Quantum GIS" trademarked?

We need to link to Quantum GIS on a major official website. My legal department asks me if "Quantum GIS" is trademarked. Is it? Thanks -- Alex


How to display Radius after Resizing a Feature in Openlayers

@captDragon has helped me with calling the coordinates of a polygon after I have resized or dragged the feature. Ineed the radius now, i have done a radius alert on my other map,Which works correctly, I,m trying to replace the following code:                   vectors.events.on({           featuremodified: featureModified           });               function featureModified(event) {              var bounds = event.feature.geometry.getBounds();         var answer = "bottom: " + bounds.bottom + "\n";         answer += "left: " + bounds.left + "\n";         answer += "right: " + bounds.right + "\n";         answer += "top: " + bounds.top + "\n";         alert(answer);         }      With this code                   function onFeaturesAdded(circle) {            var area = circle.geometry.getArea();            var radius = 0.565352 * Math.sqrt(area);            alert(radius);            }      I tried                      function featureModified(event) {            var area = event.geometry.getArea();            var radius = 0.565352 * Math.sqrt(area);            alert(radius);            }      But this does not work. Below is my working file:                                                                                <link href="css/style.css" rel="stylesheet" type="text/css" />

         <input type="radio" name="type" value="none" id="noneToggle"                onclick="toggleControl(this);" checked="checked" />                navigate

           <input type="radio" name="type" value="regular" id="regularToggle" onClick="toggleControl(this);" />         draw regular polygon

            <input type="radio" name="type" value="modify" id="modifyToggle"                onclick="update();" />                    modify feature

                   <input id="rotate" type="hidden"                         name="rotate" onChange="update()"  checked="checked"/>

                 <input id="resize" type="hidden"                        name="resize" onChange="update()"  checked="checked"/>                 <input id="keepAspectRatio" type="hidden"                        name="keepAspectRatio" onChange="update()" checked="checked" />

                 <input id="drag" type="hidden"                        name="drag" onChange="update()"  checked="checked"/>


Removing duplicate features in PyQGIS

I'm trying to remove duplicate features from PyQGIS. Manually I can do it by using "Edit -> Merge Selected Features" and then "Multipart to Singlepart" from fTools in QGIS, but can this functionality be reached from PyQGIS? Can't find anything about Merge in the API docs. I tried using QgsGeometryAnalyzer.dissolve but it seems to keep the duplicate parts in the multipart features created. I hope someone can describe an easy way to accomplish what I'm after. Regards, Andreas


Is it possible to georeference an existing, un-georeferenced pdf?

I was wondering if there was anyway to georeference a pdf directly without first converting it to an image. I have access to ArcGIS 10.1 but haven't been able to find any information that would indicate it is possible. I am willing to try other open source software if they have a solution. Thanks! Edit: I receive site plans in PDF form, often generated from AutoCAD. Currently I save the pdf as a jpg to import into ArcMap, geo-reference it and then digitize things like building footprints. I'm just wondering if there is a way to skip the conversion to jpg step.


climatic data from MODIS or other sources

I'm looking for climatic information in the form of raster like MODIS daily LST (Land Surface Temperature) data. I'm looking for other parameters like precipitation, humidity, air pressure and so on. Does anybody have information where could I find MODIS (or other sources) data for these parameters and what is the best choice for each?


UK National Grid references and GeoDjango

Newbie question. I'm trying to insert some UK National Grid references into a Django PointField defined as follows:               oscode = models.PointField(srid=27700, null=True, blank=True)      My question is this: can I store them in their original form ('TR3241', for example) or do I have to convert them to lat/lng first? I don't know how to format them as WKT. This is what I get if I try inserting them as `POINT(TR3241)`, which is obviously wrong:               INSERT INTO places (placeid, structidx, subidx, county, name, oscode) VALUES ('10', '1', '1', 'Kent', 'Dover', 'TR3241');     LINE 1: ...'1', 'Kent', 'D1', 'Eastry', 'Bewsbury', 'Dover', 'POINT(TR3...                                                                  ^     HINT:  "POINT(" <-- parse error at position 6 within geometry      So: is there a way that I can store points like 'TR3241' in GeoDjango, or do I have to convert them to standard lat/lng first? I do have code to convert National Grid > lat/lng, but I was hoping GeoDjango would be able to store them in the original format, and do any conversion for me.