Uniform Annotation Font Size in Multiple Data Frames

I'm making a map that has five inset maps that display a smaller map scale where point features are more dense. I have exported my labels to annotations in the geodatabase so that I could fit the labels that maplex wouldn't draw and move things around for the final touches. I am using the same annotation layer across all the insets so my GDB isn't cluttered with multiple annotation classes. The problem is that because these data frames use different map scales the annotations are drawing at different sizes because it is set to a single reference scale. I can't set it to 0 because it is in a geodatabase. Using Symbol Substitution changes the font size, but it doesn't scale the font to the same sizes (10pt in one data frame displays at a different size than 10pt in another data frame). My question is: **How do you display the same geodatabase annotation class in multiple data frames, with diffent scales, at the same page layout font size?**


Why Does Cold Start Take So Long?

I know very little about how GPS works. I noticed that after disabling GPS mode for a few days, even while standing on a street in the city, it can take several minutes for GPS to show my location, even after using GPS Status or GPS Test to reset and download A-GPS data over a 3G connection, and enabling Location in the Settings section. Is that to be expected, or should I try something else? FWIW, I use a Galaxy Nexus running CyanogeMod 11. Thank you.


The Hobbit Visualisation - Separate Infowindows

I'm being pushed to the max of my (albeit limited) Javascript skills so any help would be most appreciated! Forking from the the Hobbit Visualisation at http://cartodb.github.io/cartodb.js/examples/TheHobbitLocations/, I am attempting to create a set of buttons where users can select a type of filming location (e.g. villages, rivers, etc). Once one of these buttons have been clicked, I then want the hoverover infowindow on the bottom right to be unique to that specific filming location type (for example, rivers would show information such as "Length" and "Source" whereas Natural Parks would show "Year Designated" and "Terrain"). I have created the different infowindow instances, but can't seem to get the buttons to show only one infowindow at a time (at the moment they all come up one above the other!). The JS code for the buttons are as follows:               // defining button actions           $('.button').click(function() {         $('.button').removeClass('selected');         $(this).addClass('selected');         LayerActions[$(this).attr('id')]();         });          var LayerActions = {         nationalparks: function(){         sublayers[0].set({         sql: "SELECT * FROM hobbit WHERE type='National Parks' ORDER BY ranking DESC",      cartocss: "#hobbit{ marker-opacity: 0.9; marker-line-color: #FFF; marker-line-width: 1.5; marker-line-opacity: 1; marker-width: 20; marker-allow-overlap: true;}}"     });      Any help would be most appreciated! Thanks, Ben **Update** Sorry, just one more question, but the infobox on hover currently doesn't work on iPhone or iPad. Is there any workaround for this?


How to find which table has layer in mapinfo?

Is it possible to find if table has layer Or print list of layers? Edit:: I have a table when I selete some data from table like:               Select * from earth where population >100      Then it create new table but if want to hide layer closing table would not work so to hide layer I need to know table that are used to create layer?


Shapefile to table within python

I was just curious as to whats the best way within a python script to get the output a dbf file/ table while the input is a shapefile. (if possible which i'm assuming it is) Is there a way of converting a shape file into a table, or by selecting the attribute table of the shapefile and just outputting that. Or am i going to need to create a table and then import the values?


Adding a button from an add-in to a TOC context menu

I would like to add a button to a group layer context menu in the table of contents. The button is from an add-in that also contains an extension. Should be straight-forward enough. The extension monitors the IActiveViewEvents_Event events and checks for a group layer being added by a specific name. When a layer is added, the IActiveViewEvents_ItemAddedEventHandler() is fired and calls a method that adds the button to the group layer context menus if the name is right. That mechanism is working. The snippet of code that is failing is below:               // Set a uid for the group layer context menu.     UID gUID = new UID();         gUID.Value = "esriArcMapUI.GroupLayerContextMenu";   // "{863A0D98-73DC-4331-8658-ED0E22247E36}"          // Get the command bar denoting that context menu.     ICommandBar groupLayerContextMenu = ArcMap.Application.Document.CommandBars.Find(gUID, false, false) as ICommandBar;          // If there is a group layer context menu.     if (groupLayerContextMenu != null)     {         // Change the UID to the add-in command.         gUID.Value = ThisAddIn.IDs.TestContextMenuButton;              ICommandItem cmd = groupLayerContextMenu.Find(gUID);              if (cmd == null)         {             cmd = groupLayerContextMenu.Add(uid, ref missing);  'Command not available' error thrown here...             cmd.Refresh();         }      }      To me, it seems that the command is not found because ICommandBar.Add(uid) is trying to dig for it in the registry but because it resides in an add-in it has not been added to the registry. Is there another way to add a button, or any item for that matter, when it lives in an add-in? Or am I missing something simple?


Can individual feature properties be retrieved from the CartoDB.js API?

I'm building an application that uses CartoDB to show users how a proposed redistricting plan would affect them. The user enters an address and that fires a ST_Intersects query to CartoDB, the result of which is always two polygons: one is the old ward, one is the proposed ward. Type "40 Bradley Street, Burlington, VT" into the search box to see how this works. ![appscreen](http://i.stack.imgur.com/ioW6z.jpg) I'd like to push answers to the indicated questions (e.g. "Where will my polling place be?") as well. However, this requires the ability to retrieve individual attribute values from the query results, and I can't figure out how to do this. Here is the structure I'm trying to use after running cartodb.createLayer();               var 1990ward = cartodbproperty.ward[1];         var 2010ward = cartodbproperty.ward[0];              if (1990ward == 2010ward) {             $('#wardchange').text("No. You will remain in ward " + 2010ward + " under the proposed plan.");         }         else {             $('#wardchange').text("Yes. You are currently in ward " + 1990ward + ", but you would be in ward " + 2010ward + " under the proposed plan.");         }              $('#distchange').text("The " + cartodbproperty.district + " District");         $('#pollchange').text(cartodbproperty.pollingplace);      The "cartodbproperty" variables are dummies because I can't figure out how I would use the cartodb.js API to select the "ward" attribute value of the query-result feature where "year" = 2010. I believe I could pull this off with a parallel call to the SQL API - which could then be parsed as JSON - but that seems horribly inefficient. How can I get this data from the CartoDB.js API? Thanks as always for the great toolset!


Why does CalculateField set values to <Null> when output extent is in a different coordinate system?

I'm quite new to ArcGIS, and now I have to work on a script that goes through a list of feature classes to:   1. add a field to their attribute table and   2. fill that field with a number that I retrieve from somewhere else When I execute the following statements in the ArcMap Python window or in a minimal standalone script, it works just fine:               table = "JS4A"     arcpy.AddField_management(table, "BK_Code", "LONG")     arcpy.CalculateField_management(table, "BK_Code", "1", "PYTHON_9.3", "")      I can then open the attribute table of that feature class in ArcMap and see that this field has been added and populated with all 1s, just as expected. Now the weird thing is, when I run the exact same code in the complex production script, the field is added as well, but when I open the attribute table, all the values are `<Null>`. I managed to track the problem down to the following line of code in my production script, that comes before the part with AddField and CalculateField:               env.extent = desc.extent      where `desc` is the return value of an `arcpy.Describe()` call to some other, unrelated shape file. All `CalculateField_management` calls before that line will work fine, all calls after that line produce `<Null>`. When I replace the above line with a line that sets the output extent to a hardcoded value, it works. I don't understand this at all. Is this an ArcGIS bug? Unfortunately I can't omit the above line because it is required in the production script.


ArcGIS - make linear regression across multiple raster layers

I have a series of 5 rasters which represent a value for 5 consecutive years. Is there a way to create a linear regression for each cell over the 5 layers and output the slope? The result would be a raster with cell values equal to the slope of the linear regression line. Essentially it is exactly what is done in this post: How to represent trend over time? thanks


create a map from scratch

I am a newbie in GIS, so be indulgent if my question sounds stupid! I am working on a mapping project for scuba diving sites, and I need to create 3d maps viewable by users on android tablets, iOS, websites, etc... So far, I achieved the map result I expected using GoldenSoftware "Didger+surfer". The problem now is to eport the 3d map in a format manageable by a vrml java viewer. Is there a way to do this view QuantumGIS, or any other recommendation? Thanks for your help


How to run Arc Hydro in ModelBuilder?

I tried to run a ModelBuilder program but it gave the error below. ![enter image description here](http://i.stack.imgur.com/Zc5EM.png) I am working in ArcMap 10.2, windows 7 64 bit pro and the raster format is Grid. When installing ArcGIS 10.2, it required me to activate extension (customize/extension/select all extension). After the install the Arc Hydro toolbar does not add to extension to enable it to be active. I think the above error is because Arc Hydro has not been activated. What is your opinion?


GeoDa crashes when loading shapefiles

I'm trying to load around 6000 shapefiles (all the postcodes for a city) into GeoDa but it crashes everytime ("OpenGeoDa.exe has stopped working"). Their troubleshooting page details this problem but I've not been able to resolve the issue, as stated here: http://geodacenter.asu.edu/node/404#load I've tried to code a new id field in Arcmap but that did not work (was I right in recoding the new ID variable as a short integer?). I then tried in Arcmap to see if the field names were the problem. I cleaned up most of them up, making them shorter and removing underscores. However, I have been unable to change 5 of the fields (could they be read only?) and the names are over 30 characters long. Could this be the problem? Please see the screen grab below for an overview of my data and the length of one of the field names. Any help would be very much appreciated! ![enter image description here](http://i.stack.imgur.com/2mKvB.png) Edit: Screen grab number 2. Note: I reduced the number of fields to 5 and only selected a few polygons to see if size was the issue. The GeoDa error message on the right is displayed (still the same message). ![enter image description here](http://i.stack.imgur.com/FvOLh.png)


Tell Dojo to look for minified classes with ArcGIS JavaScript API

I've got an application build using the ArcGIS JavaScript API. I have several classes built with Dojo define. Each class is in a separate file, and everything is working fine. However, I've minified my classes (using Visual Studio Web Essentials), and I'd like to tell Dojo to use the minified files. I'm using a dojoConfig like the following:               var dojoConfig = {       isDebug: false,       packages: [{         name: 'js',         location: '/js'       }]     };      And dojo loads my files like '/js/MyClass.js'. What do I have to do the get Dojo to load MyClass.min.js ? This would be a successful resolution to my problem. But, even more, I have gzipped versions of those files, so how can I tell Dojo to download MyClass.min.js.gzip? I'd even like to take it a step further and bundle all of my classes into one download, say MyBundle.min.js.gzip. That not a problem. But, again, how to explain to the Dojo AMD loader that /js/MyClass1 and /js/MyClass2 are both in MyBundle? All this, given that I'm not using the Dojo build system.


ArcGIS JavaScript API: How to add a CSV as a Layer?

I've created a map on ArcGIS online, adding a layer by uploading a csv file. This works fine, see map at link http://www.arcgis.com/home/webmap/vi...1ba43a0256c06b. Now I'm trying to simulate this in JavaScript API. First I create a Web API service that returns the same data set I've got in csv file. This is published on my local webserer. In my web page, I create a webmap, adding a layer that references this web API:                function init() {      var options = {      basemap: "gray",      center: [153.115672, -30.296312],      zoom: 14      }      var mapServiceURL = "http://[myserver]/TRWDataManagerAPI/API/OzRegs";           // Create map and call Init function      map = new esri.Map("mapDiv", options);           map.addLayer(new esri.layers.ArcGISDynamicMapServiceLayer(mapServiceURL));      }      dojo.ready(init);      The webmap is displayed, but the layer doesn't have any effect. What's missing? Am I at the right track? Your help is greatly appreciated.


High Resolution Sinusoidal Interrupted Projection of the Earth

My grandfather built a rather large globe 10 years ago and never got to complete it (he started drawing about a fifth of it by hand and just gave up after that). I have been looking for a sinusoidal interrupted projection, with a traditional linework (but any would do) that I could get printed and apply to his globe to finish his work for him (as a surprise!). I have contacted a few companies but they are charging around 3000 EUR (linework and background topography image) which is way too much for me. I am a bit lost with the science of it and would really appreciate if someone who understands that stuff could help me out of bit : Basically where kind I find a high res image of the earth with country and city names, borders, etc... and how do I make it into a sinusoidal interrupted projection? Data I got from my grandpa a couple years ago, when I was shopping around : \- globe circumference = 2 meters \- diameter = 0,63662m \- scale = 1/20 M \- 10° L=111km111111/20M =5cm555 Thanks!


Building cartograms using Web GIS?

I would like to build a cartogram map using Web GIS. I have a basemap and a set of raw data in excel. I wondered what are the steps to build up a cartogram map using Web GIS? **Note** : Originally this Question solicited Answers for both Web GIS and ArcGIS for Desktop. A new Question has been asked focussed on ArcGIS for Desktop (Creating contiguous pixel based cartogram using ArcGIS for Desktop?) and so this one has now been focussed on Web GIS (which itself is still quite broad).


Raster Over Layer in Leaflet from GeoServer

I want to add a raster over layer onto the base layer and the raster layer is published by a Geoserver WMS. I set the raster layer's geographic coordinate system as the same as the cloudmate tiles. The raster source is a layer of the United States in the GeoServer sample data. I can get the layer in Leaflet and do spatial queries, like query a part of the raster by bounding box successfully, which means the WMS is good. But the raster image never overlays correctly. I've set the same GCS and projection. Why does that happens? Is there any successful instruction or sample about it?


Web based spatial data editor

Is possible to set up web access to GIS data such that a user could login and see the collected data and be allowed to suggest/draw changes but not be able to edit the data files? My idea is that community members could see data that has been collected through land use interviews, then be able to suggest changes to the project maps through a web based editor. A GIS tech would then review the suggestion to see if the change would be appropriate. One issue we have is that our data is considered highly proprietary and confidential. So security is very important. Thanks.


Is there a quick way to invert a node selection in QGIS?

is there any way to invert a node selection in QGIS? For example, I want to extract a small section of a long line - if I select the nodes I want to keep, is there any way to invert that selection so I can delete the unwanted nodes? The only way I can do it in v1.8 is to zoom all the way out, select all the nodes, zoom back in and deselect the ones I want to keep, then hit delete...which is a bit slow going... Just wondered if there's a quicker way with a tool or keyboard shortcut?


Openlayers getDataExtent not working

I am following this tutorial http://workshops.boundlessgeo.com/tutorial- wordmap/ on the web about using setting up sql view on geoserver. I have set up my sql view correctly and have the javasctipt part working. When I enter a location in the textbox and press enter, the map will only display features of that particular location. However, I would like to go further so that my map can also zoom into the features of that location. This is my code:               // Start with map of startWord     wms11.mergeNewParams({viewparams: "word:"+startWord});               //Add all layers to the map      map.addLayers([osm,wms11]);          // Text field component. On 'enter' update the WMS URL     var textField = new Ext.form.TextField({         value: startWord,         listeners: {         specialkey: function(fld, e) {             // Only update the map when the user hits 'enter'             if (e.getKey() == e.ENTER) {                 wms11.mergeNewParams({viewparams: "word:"+fld.getValue()});                 var bounds = wms11.getDataExtent();                 alert(bounds);                 map.zoomToExtent(bounds);});                 }             }         }     });       The map was zoomed into an incorrect area and the alert(bounds) shows the bounds is undefined. What is going wrong here? Thanks in advance!


Vector/Raster Intersection - Speedup and replacement for ST_Centroid

I have some raster data which is the altitude above sea level and also some vector data which represents building outlines. What I want to do is grab some buildings outlines and fetch all the raster points inbetween those buildings and their value - so all raster points that do not fall into buildings basically. What I have is this function, which should give me all points and all values for the raster points within some precalculated polygon, where the inner rings are the building outlines actually.               with vars as (select ST_GeomFromText('POLYGON((...),(...),(...))',...) as poly)     select ST_AsText(ST_Centroid(ST_Transform((inside.geomval).geom,...))) As geomwkt,     (inside.geomval).val from (select ST_Intersection(rast, (select poly from vars))      AS geomval from rasterlayer where ST_Intersects(rast, (select poly from vars))) as inside;      Note that the ST_Centroid(...) is required to reduce the polygons retrieved to single points - this seems a bit weird to me, but it works. Explain analyze gives me the following:               "Subquery Scan on inside  (cost=0.05..5630.94 rows=719 width=32) (actual time=279.400..6075.326 rows=2846 loops=1)"     "  CTE vars"     "    ->  Result  (cost=0.00..0.01 rows=1 width=0) (actual time=0.003..0.004 rows=1 loops=1)"     "  ->  Seq Scan on rasterlayer  (cost=0.04..5618.35 rows=719 width=32) (actual time=279.363..6020.237 rows=2846 loops=1)"     "        Filter: st_intersects(rast, $2, NULL::integer)"     "        InitPlan 2 (returns $1)"     "          ->  CTE Scan on vars  (cost=0.00..0.02 rows=1 width=32) (actual time=0.005..0.005 rows=1 loops=1)"     "        InitPlan 3 (returns $2)"     "          ->  CTE Scan on vars  (cost=0.00..0.02 rows=1 width=32) (actual time=0.021..0.023 rows=1 loops=1)"     "Total runtime: 6076.956 ms"      Although I imported the raster data in an indexed fashion (and rasterlayer consists of 50x50 meter tiles), I feel that the query is pretty slow. If I call my function for a large area, it takes hours or days. Can anyone provide me with a more efficient solution for the ST_Centroid(...) issue and tell me how to improve performance? Any help in speeding up the whole process is highly appreciated.


Openstreetmaps API

I want to get few results from Openstreetmaps data. Such as I will submit a city/postal code and in return I want to get all street name within the city/postal code. This can easily done by full OSM setup and PostGIS but for this few need I do not want to setup huge database. Is there any way to retrieve the data? There is an old API called openstreetmap API v.0.6. But I am not sure whether it can fulfill my purpose or not.


Data security for GIS Data (Esri shape files,Auto-cad files)

This is regarding security for GIS Data (in Esri & open source) (Esri shape files,Auto-cad files etc.) Does anyone know how to assign a password to or encrypt shapefiles OR CAD Data? I have found an ESRI script on this site but I think its too old..So if anyone knows such information/tool then please let me know... Thanks in advance : )


How to get Google Street View panoramas on the same street as the requested address?

The Google Maps API geocodes an address to a point. I can use that point to find the nearest panorama (street view) image. I can then determine the heading based on the two points and show the panorama pointed at the right house. The problem is that in some cases such as corner lots, the closest panorama might be on a different street. For example, a house with an address of 654 NW 4th Street is on the corner of 4th and Greenwood. The closest panorama is showing on Greenwood instead of 4th so I get a picture of the side of the house instead of the house. Am I doing something wrong here? Is there a way to get a list of panorama's so that I can find one that is on the correct street? It's odd because the correct panorama shows if using Google directly. E.g.: `https://maps.google.com/maps?q=654+NW+4TH+ST,+REDMOND,+OR+97756` Here is a JSFiddle demonstrating the issue.


json to xml transformation with x2js

I am trying to transform a JSON object to XML, but i have somehow get stuck in jscript (not an expert). Here is the code:               ...     var GEOJSON_PARSER = new OpenLayers.Format.GeoJSON();     var vectorLayerAsJson = GEOJSON_PARSER.write(vlayer.features);         ...     var x2js = new X2JS();         ...     var xmlDocStr2 = x2js.json2xml_str(vectorLayerAsJson);     alert(xmlDocStr2);         ...     // JSON arrays to string     var xmlDoc1 = x2js.json2xml_str({         "type": "FeatureCollection",         "features": [{             "type": "Feature",             "properties": {},             "geometry": {                 "type": "Polygon",                 "coordinates": [                     [                         [22580, -24812],                         [24772, -28076],                         [27732, -24204],                         [22580, -24812]                     ]                 ]             }         }]     });     alert(xmlDoc1);      In the first case the json2xml_str(vectorLayerAsJson) function returns something like a JSON object however in the second it returns the correct XML. Both objects are identical "vectorLayerAsJson" is the one declared as `"type":"FeatureCollection ..."` in the second case. I have also tried alternate ways of calling `x2js.json2xml_str` such as `x2js.json2xml_str($.parseJSON($("vectorLayerAsJson").val()));` with no success! Any ideas?


Equivalence of EPSG:3857 for MapInfo 9.5

I have copy a raster projected on EPSG:3857 with QGIS. I have edit image registration in Mapinfo but there is no projection for EPSG:3857. Do you know the projection who I need to use? Thank you. Jérôme


Can I convert raster to shape by specify class?

I use the method "IMultivariateOp.MLClassify",so I get many class such as the figure below, ![Convert](http://i.stack.imgur.com/Xxr5E.png) can you give some code in C# to convert each of them to shape?Thank you. one more question, why I just get the first three item?here is my test code:               IRasterLayer rl = (IRasterLayer)map_con.get_Layer(0);      IRaster r = rl.Raster;      IRasterBandCollection rbc = r as IRasterBandCollection;      for(int rc = 0; rc < rbc.Count; rc++)      {         IRasterBand rb = rbc.Item(rc);         ITable t = rb.AttributeTable;         int fileIndex = t.Fields.FindField("Count");         for(int i = 0; i < t.Fields.FieldCount; i++)          {           writer.WriteLine(t.GetRow(i + 1).get_Value(fileIndex));          }      }


Snap line layer to network in QGIS or PostGIS

I have GPS data I've taken of bus routes, and now I'd like to snap it to my road network. Both layers are line layers in a PostGIS DB. I'd like to use either QGIS or PostGIS, but if I have to go with GRASS or ArcMap, that's OK, too. Thanks! To clarify, I'm trying to snap lines to lines, not points to lines.


What function should I use to assess the degree of elevational change (rocky cliffs etc.) of a DEM?

I have a DEM of my research area and I want to assess how the topography/elevation changes. I am specifically looking to see what areas may be good suitable bat-roosting habitat so I am trying to identify which areas might be rocky hillsides or cliff-faces. I am wondering what kind of function should I use to best assess this? Slope? Ruggedness? I am also wondering HOW to go about utilizing the best function to do the job. I am a beginner GIS student/ecologist, any advice or help would be fantastic! Thanks in advance!


What is a radar projection?

Trying to overlay a radar image(radar projection) on a mercator projected map, but have no idea about what a radar projection is. So here are my quesions: What is a radar projection? What's its forward/backward formula? Does proj4 or gdal/ogr support this projection?


Zoom Commands as Mouse Gestures

Is it possible to use the zoom in, zoom out, and previous extent functions as mouse gestures in ArcMap? For example, say a user held down LShift then left click and dragged diagonally to the lower right. I would like that to function the same as zooming in, using the zoom in tool. Dragging in a different direction would do one of the other functions. A final function would be escaping from the command after a drag has been started, this would be the same as hitting the Esc button. I assume this would have to be done in ArcObjects. I've done some searching around for zooming functions in ArcObjects but most of the stuff I've found refers to older versions. I'm not at all familiar with ArcObjects, so some pointers would be very helpful. I have access to both Arc 9 and 10.


Why does my 113th Congress Shapefile match up with the wrong districts across the country?

I downloaded a Shapefile from the Census on the 113th congress broken down by district nationally: `http://www2.census.gov/geo/tiger/TIGERrd13_st/nation/tl_rd13_us_cd113.zip` Within the .dbf file, there were no good indicators of how the districts are normally laid out (Eg. PA-15); each state was assigned a number and then a congressional district, so I figured out what the state codes meant (based on other census shapefiles and how many districts were in each state and made a new column that I'm 99.9% sure is right based on the information given. The first four rows are what came with the shapefile and the last highlighted row is what I added: Modified DBF.jpg And here's what it looks like when I run a query for the districts that I've defined as Pennsylvania: "Pennsylvania" GIS.jpg Pretty sure that's not what Pennsylvania looks like! Can't figure it out -- because, again, I'm really sure that I defined each state's number as the proper actual state. Any ideas as to what might be wrong? And sorry about all the links to pictures -- since this is my first post, I'm not allowed to post any images or more than two links yet. Thanks!


Basic Question: What is amp datasets?

I have a job interview tomorrow for a summer student position working with GIS. I will be mainly responsible for assisting in the acceptance, storage and dissemination of geographic information, amp datasets,and associated metadata. What are amp datasets? Is it a fancy way of saying data? Also, what would i be expected to know about metadata? Thanks!


How to know the exact country or province from the longitude and latitude information automatically?

To be more specific, how to find all of point (from 1,000,000 items) located in England, or located in California? I need to process over 1,000,000 items, so a offline and automatic way is needed. The geographical information is like this: 39.2903848 -76.6121893 Is there any database can deal with this?


Arcpy: print to PDF with Vector shapefiles, text etc

I have a small python script which exports a set of data driven pages to PDF. This works fine, including appending each page to an 'all' pdf. My issue is that it creates rasters of each layer in the map. I can use something like CutePDF to print to a vector map, and that works great. Much smaller filesizes, but this is not easily done in a batch process. It requires you to choose the file name of each page manually. My question is this: Can you make Arcpy's ecport to pdf method actually use vector format for the content of the map. If not, is there a way to automate the printing to file, using cutePDF or something else? Cheers, Alex


Geoserver startup page

I have been looking for this for a while and can't find the answer: when you startup geoserver a webpage is opened at localhost:8080/geoserver/web but if you have changed your port config in jetty.xml then this is the wrong address... still, when you start up geoserver opens the page at 8080, which fails now of course... here is my question: where is this behaviour set in geoserver? while i am at it, is it possible for geoserver to autostart? thanks in advance update: this is for running geoserver on an OSX Machine


Which format is more convenient to publish in Erdas Apollo?

I am from Bolivia and I need to publish some Satellite images in Erdas Apollo. My question is Which format is more convenient for publish in this software?, because I have Landsat5 TM images in TIFF format and i'm watching the possibility to change or keep the format.


Using Query Layers with ArcGIS Server 10.0 for Java on Linux

I have an ArcGIS Server 10.0 for Java running under Red Hat Linux. I would like to be able to publish a Map Document that I have created which has a working Query Layer connecting to Oracle 11g. In ArcMap the layer works as expected. I have installed all the Oracle client software on the ArcGIS Server (the server is actually running Oracle also). I am able to connect via SQL Plus as the ArcGIS server user. I have set all the environment variables to the best of my knowledge. However when I publish the Map Document I get the following in the log WARNING The Layer:'XXXXXXXX' in Map:'Layers' is invalid. No SDEHOME environment variable set WARNING The Layer:'XXXXXXXX' in Map:'Layers' is invalid. No SDEHOME environment variable set INFO3 Server Object instance on machine XXXXXXXX is created with errors. The Layer:'XXXXXXX' in Map:'Layers' is invalid. No SDEHOME environment variable set INFO3 Server Object instance on machine XXXXXXXX is created with errors. The Layer:'XXXXXXX' in Map:'Layers' is invalid. No SDEHOME environment variable set I don't have ArcSDE installed and I thought the idea of Query Layers was that you didn't need it? Very confused and Google didn't turn up much :-( Any help would be most appreciated. Regards Colin


ArcMap Caching Character Problem

i have a basemap which was created with arcmap and i want to use it for my application with cached. in my viewport in arcmap there is no any Turkish character problem on it but when i cache my work it makes character problem. how can i accomplish this subject with settings ? thank u.


custom control on google maps API

I added 4 buttons to a map which zooms in to particular state in country on click using google maps api. I am a beginner and did it using example from google api custom control documentation. What i want to do is to create/load more controls/buttons for 10 big cities in each state on each button. i have data for each city in mysql. i will try to load the data on click then. Thank you in advance


ArcGIS Geoprocessing "None Type" Second time

I am trying to publish an ArcGIS geoprocessing service in 10.2 I have a script setup that will run multiple times through Pycharm, but when I run it in ArcGIS it will run once and return this error: > Traceback (most recent call last): File > "C:\Users\Desktop\git\GeoprocessingService.py", line 477, in main() File > "C:\Users\Desktop\git\GeoprocessingService.py", line 372, in main cur) File > "C:\Users\Desktop/git\GeoprocessingService.py", line 136, in value_query > results = cur.fetchall() File "C:\Python27\ArcGIS10.2\lib\site- > packages\psycopg2\extras.py", line 81, in fetchall res = > super(DictCursorBase, self).fetchall() File > "C:\Python27\ArcGIS10.2\Lib\decimal.py", line 543, in **new** m = > _parser(value.strip()) TypeError: 'NoneType' object is not callable > > Failed to execute (Script). The script itself is 800 lines, but here is the query I have been running that errors out:               def value_query(id_field, id, value_fields, start_date, end_date, date_field, table, cur):         query2 = "SELECT "         for field in value_fields:             query2 += "NULLIF({0}, '-999')::NUMERIC(6,2),".format(field)         query2 += "{0} " \                 "FROM {1} " \                 "WHERE {2} = '{3}' AND {0} BETWEEN '{4}' AND '{5}' " \                 "ORDER BY {0} ASC;".format(date_field, table, id_field, id, start_date, end_date)         #print "This is Query:"         #print query         arcpy.AddMessage("{0}".format(query2))         cur.execute(query2)              results = cur.fetchall()              #print "These are value query results:\n\n {0}".format(results)         return results      Like I said it runs once in ArcGIS then errors out. However, if I close Arc completely and run the exact same script again, it runs once then stops. I feel like there may be some sort of caching or something that Arc does, but am completely stumped.


QGIS/PostGIS: Select Counties within Region via Spatial Query

I'm using GQIS 2.2 and PostGIS on Postgresql 9.3. I have 5 regions in the contiguous United States. Each region has states, and each state has counties. I would like to think that the data I'm given has each region, state, and county line up cleanly, but one can never be sure. :) Here is an example: ![Regions have States, and States have Counties](http://i.imgur.com/xronc1n.png) I'm at the point in the project where I need to select all counties within a region. Success is getting all the counties within a given region (say, the Upper Midwest). None of my attempts at QGIS or PostGIS give me the results I want. It may either be because my queries are poor (I'm a beginner), or because my data is dirty, or both. ![An image with the midwest highlighted. I want to get all its counties](https://i.imgur.com/5C3MZKf.png) Here's what I've done so far: I've created PostGIS queries and also tried using QGIS's interface. ST_Within: ![ST_Within results](https://i.imgur.com/CXysYrW.png) Using QGIS's interface, we can see that ST_Within misses the border counties. ST_Touches: ![ST_Touches results](https://i.imgur.com/pnytmEw.png) Touches seems to fail pretty badly - none of the counties are actually the counties that I want. They seem to be outside the selection! This is what leads me to think that I may have to clean my data somehow. ST_Intersects: ![ST_Intersect](https://i.imgur.com/a4LFuoQ.png) Intersects is another step in the right direction, but it overselects. I also tried to get fancy using centroids, but I hit a few interesting edge cases involving centroids falling outside of county polygons and islands. ![Using Centroids](https://i.imgur.com/tugqhXl.png) I'm open to any advice. Of course, I could always, you know, just use manual selection, but I've had far too much coffee today to just walk away from an interesting problem. If the problem is just of selection, then fine. But I think the real problem is that my data is off. I'm open to any data cleaning techniques you have. I assume that I'm missing a general technique. I've read How to select by polygon in PostGIS with misaligned layers, and while the solution was indeed very cunning, it doesn't strike me as solving the underlying problem :) However, I am open to any suggestions. Thank you for reading.


Address to Zip+4 ArcGIS Address Locator or batch service? Free?

I have 200,000 addresses in MD for which I need zip+4. I am searching for an address locator for ArcGIS or a (free?) batch service that will output 9-digit zip (Zip+4) from address. I have already geocoded these addresses so I have lat/long from street or point address for a good portion. I do have Esri Street Map Premium service and that allows for 9 digit zip to Lat/long but not Address to Zip+4. Odd. There is a "PostalExt" field in the output but i does not populate. very annoying. Reference data must be available for the boundaries (or at least centroid) for this data but it is not available for use in the street map premium FGDB. Since I'm only interested in one state (MD), my plan was to reverse geocode all permutations of possible 9 digit zips (about 6 million) to get their lat/longs, but it seems quite dirty to do a spatial join on these points to the previously geocoded 200k address points. Any other suggestions very welcome. I have seen that USPS does not have or make available polygons for zip+4 but you can purchase tables of the address ranges per +4. Not sure if that would even be that helpful for me. Thanks :)


ArcGIS Geoprocessing Service gives Error 000816: Tool not valid?

I have developed few ArcGIS python toolboxes. These tools work fine on ArcGIS for Desktop (ArcMap) but gives **Error 000816** when I publish these as geoprocessing services. Can some one help me with the cause of the error?


How to dissolve naturalearth vectordata with gdal?

trying to use the workflow described by @aflaciano at How to render area that crosses 180°? I did the following: 1) downloaded and installed gdal-111-1600-x64-core.msi from http://www.gisinternals.com/sdk/ on my windows 7 64 bit pc. 2) downloaded data from http://www.naturalearthdata.com 3) used the following batch-file adapted from @aflaciano:               rem Clip world.shp between -169° and +180° of longitude in order to exclude the Chukotka Peninsula (1st part):     rem ==========================================================================================================     ogr2ogr world_clip.shp ne_110m_admin_0_countries.shp -clipsrc -169 -90 180 90     rem Shift and clip world.shp between +180° and +191° of longitude in order to include the Chukotka Peninsula (2nd part):     rem ====================================================================================================================     ogr2ogr Chukotka_Peninsula.shp ne_110m_admin_0_countries.shp -dialect sqlite -sql "SELECT sov_a3, ShiftCoords(geometry,360,0) FROM ne_110m_admin_0_countries" -clipsrc 180 -90 191 90     rem Merge the two parts:     rem ====================     ogr2ogr world_shifted.shp world_clip.shp     ogr2ogr -update -append world_shifted.shp Chukotka_Peninsula.shp -nln world_shifted     rem Finally, dissolve countries boundaries of world_shifted.shp obtaining world_output.shp:     rem =======================================================================================     ogr2ogr world_output.shp world_shifted.shp -dialect sqlite -sql "select ST_Union(Geometry),sov_a3 FROM world_shifted GROUP BY sov_a3"      at step 3 last line of the batch-file I must have done an error, because dissolve does not work for me, but I get no error message: ogr2ogr world_output.shp world_shifted.shp -dialect sqlite -sql "select ST_Union(Geometry),sov_a3 FROM world_shifted GROUP BY sov_a3" I tried also other variables for dissolve , e.g. sovereignt or admin, geounit, adm0_a3 ![enter image description here](http://i.stack.imgur.com/0m5Yo.jpg)


ArcObjects 10 - copying fields (table structure) to a new feature class

I'm creating a new feature class with the attribute table schema (field names, types etc.) copied from another feature class. Then, I'd like to modify this schema slightly (e.g. add an additional field). What is the best way to achieve this? Should I loop through all input fields and create respective output fields? Thanks in advance, Jan


Load Featureclass from ArcMap to a Database (sde) with Arcobjects

I have to develope an extension for ArcMap 10.2 that has to get spatial data from an Oracle Database, edit it and finally upload the changes. I create a scratch featureclass whit the data, I edit the layer in ArcMap whithout problems, but when I have to load the features in the database I can´t. I have tried IObjectLoader´s LoadObjects method, IFeatureDataConverter´s ConvertFeatureMethod and IFeatureClass´s CreateFeature method but I get a "Field not nullable" error. I think it would be a lock problem, but I cannot check it because when I do a QueryInterface to the ISchemaLock interface with the featureclass I have an error because is not a ISchemaLock (I don´t know why). This is the code:                           IFeatureClassName pSourceFCName = new FeatureClassNameClass();                 pSourceFCName.FeatureType = sourceFeatureClass.FeatureType;                 pSourceFCName.ShapeFieldName = sourceFeatureClass.ShapeFieldName;                 pSourceFCName.ShapeType = sourceFeatureClass.ShapeType;                 IDatasetName pSourceDSName = (IDatasetName)pSourceFCName;                 pSourceDSName.Name = sourceFeatureClass.AliasName;                      IFeatureDatasetName pTargetFDSName = new FeatureDatasetNameClass();                 IDatasetName pTargetDSName = (IDatasetName)pTargetFDSName;                 IFeatureClassName pTargetFCName = new FeatureClassNameClass();                 pTargetDSName = (IDatasetName)pTargetFCName;                 pTargetDSName.WorkspaceName = pTargetWorkspaceName;                      IFields pTargetFields;                 IFields pSourceFields;                 IFieldChecker pFieldCheck = new FieldCheckerClass();                 IEnumFieldError o;                      pSourceFields = sourceFeatureClass.Fields;                 pFieldCheck.Validate(pSourceFields, out o, out pTargetFields);                           IFeatureDataConverter pFeatureDataConverter = new FeatureDataConverterClass();                 pFeatureDataConverter.ConvertFeatureClass(pSourceFCName, null, null, pTargetFCName, null, pTargetFields, "", 1000, 0); --->> ERROR      Any help? Thanks in advance Víctor


How to open a shapefile using Python in QGIS?

I have tried to open my shape files with the help of a tutorial (http://spatialgalaxy.net/2012/01/27/qgis-running-scripts-in-the-python- console/) but unfortunately I couldn't open my shapefiles for my first work in QGIS. I put the directory of the shapes and the script in my python path how the tutorial recommended. So above are two pictures showing the problem by script and in command prompt: ![enter image description here](http://i.stack.imgur.com/gJbZ4.jpg) ![enter image description here](http://i.stack.imgur.com/ObwQp.jpg) Are there any other solutions (in a form of python language) to open shapefiles? I am having headeches because when I'm trying to do like the website exposes I couldn't open the shaoes and doing the commands of importing Loader I got messages of errors. Hope someone here help me! /////////////////////////////////////////////////////////////////////////////////////////// I have tried this afertnoon reproduce all the process dicussed here. I put the the script in the folder of my python path I think is correct. As I said any times later when I use the command "sys.path" many directories appears showing where are my "python paths" and after that I resolve put the script in "C:\OSGeo4W\bin" but I am not sure if here is the correct local. The photo below show what happens when I try to load the shapefiles by python: ![enter image description here](http://i.stack.imgur.com/RcYME.jpg) @gene, I follow your instruction, but how I update the question a mistake persists in the running process. First I have many "python paths" and I don't know where do I have to put the script, second when I use the command ldr.load_shapefiles('F:/qgis_1/exercicio_buffer_L1') I recive an error message like this in the photo. What do I have to do? ///////////////////////////////////////////////////////////////////////////////////////////


Static map visualization of multiple categories/attributes in a point?

I have a point data with multiple categorical values for each point location. For example: > point1 = a,b > > point2 = a > > point3 = b,c > > point4 = a,b,c,d,e > > and so on Note that the values are categorical/nominal, therefore, I cannot create some form of ranking among data points. What map visualization options can I do with this type of data?


Android ArcGIS Javascript 2.2

The website claims that touch gestures are supported for Android browsers. I'm running 2.2 and have not been able to get touch gestures to work. What can be done to enable it? This works for iPhone but not for Android browsers. How would one go about making it possible to enable touch support on Android browser? Here's a link you can try on your mobile Android browser: http://help.arcgis.com/en/webapi/javascript/arcgis/demos/mobile/mobile_simplemap.html


How can I programmatically project the Z value of a z-aware feature point, using ArcObjects in .NET/C#?

I discovered the following code leaves the Z coordinate of featurePt untouched:               // if not WGS84 project it     if (!SpatialRefsEqual(featureSpatialRef, srWgs84))     {         // featurePt here is a Z-Aware instance of IPoint         // use the Geometry Interface to do the actual project         featurePt.SpatialReference = featureSpatialRef;         featurePt.Project(srWgs84);     }      I see the IVerticalCoordinateSystem interface but it unclear to me how to use it to do a projection. For example, simply plugging it in via:               featurePt.Project(myInstanceOfIVerticalCoordinateSystem)      still leaves Z untouched. It also seems like I could just write a routine to convert the Z-coordinate in the linear units of the layer to my desired linear units (meters in this case). Is there a best practice for projecting the Z-coordinate? Thanks in advance, Josh **EDIT** According to mkennedy at ESRI, ArcGIS does not provide a true vertical transformation (thank you for answering!). If anyone knows of any outside library to vertically transform among vertical datums please let me know. Thanks!


Use OpenLayers WMSGetFeatureInfo on Mapserver Layer?

I'm able to use OpenLayers' OpenLayers.Layer.Mapserver declaration to generate a raster from GRIB data. I've read in the MapServer docs and on forums that it should be possible to get point information back from the raster using OL's WMSGetFeatureInfo. I've tried a few combinations of others' reported solutions. However, as the control is for a regular WMS and also nearly all of the solutions involve GeoServer, I've pretty much reached a wall. So this is my code:               <script type="text/javascript">          var map, info, wmsmap;               function init() {              map = new OpenLayers.Map("map", {                  projection: new OpenLayers.Projection("EPSG:3857"),             units: "m",             maxResolution: 156543.0339,             maxExtent: new OpenLayers.Bounds(                 -20037508, -20037508, 20037508, 20037508.34),             displayProjection: new OpenLayers.Projection("EPSG:4326"),             controls: [                 new OpenLayers.Control.Navigation(),                 new OpenLayers.Control.KeyboardDefaults(),                 new OpenLayers.Control.PanZoomBar(),                 new OpenLayers.Control.Scale(),                 new OpenLayers.Control.Attribution()             ]              });              var osmLayer = new OpenLayers.Layer.OSM("OpenStreetMap");          // Test MapServer WMS layer              wmsmap = new OpenLayers.Layer.MapServer( "grib plot", "/cgi-bin/mapserv",                  {mode: "map",                  map: "wxmap/glmpvis.map", layers: "mosaic", format:'image/png',              transparent:true, reaspect: "true"},                      {singleTile: true, ratio: 1, isBaseLayer: false, opacity: 1, visibility: true,                  projection: new OpenLayers.Projection("EPSG:3857")});          // add layers to map        map.addLayers([osmLayer, wmsmap]);          // add controls             map.addControl(new OpenLayers.Control.LayerSwitcher());        map.addControl(new OpenLayers.Control.Permalink());             map.addControl(new OpenLayers.Control.MousePosition(            {div: document.getElementById('ex_llmousectl'),             id: "ll_mouse", formatOutput: formatLonlats} ));        map.addControl(new OpenLayers.Control.MousePosition(             {div: document.getElementById('utm_mousectl'),              id: "utm_mouse", prefix: "UTM ", displayProjection: map.baseLayer.projection, numDigits: 0} ));        map.setCenter(new OpenLayers.LonLat(-117, 33).transform(             new OpenLayers.Projection("EPSG:4326"),             map.getProjectionObject()),6);                  info =  new OpenLayers.Control.WMSGetFeatureInfo({                               url: "/cgi-bin/mapserv",                          maxFeatures: 1,                          title: 'Identify features by clicking',                          queryVisible: true,                          layers: [wmsmap],                          vendorParams: {                                   map: "wxmap/glmpvis.map",                               },                          eventListeners: {                              getfeatureinfo: function(event) {                                  map.addPopup(new OpenLayers.Popup.FramedCloud(                                          "chicken",                                          map.getLonLatFromPixel(event.xy),                                          null,                                          event.text,                                          null,                                          true                                   )); //map.addPopup                                        } // getfeatureinfo                          } // eventListeners             });             map.addControl(info);        info.activate();             function formatLonlats(lonLat) {             var lat = lonLat.lat;             var long = lonLat.lon;             var ns = OpenLayers.Util.getFormattedLonLat(lat);             var ew = OpenLayers.Util.getFormattedLonLat(long,'lon');             return ns + ', ' + ew + ' (' + (Math.round(lat * 100) / 100) + ', ' + (Math.round(long * 100) / 100) + ')';         }               } // end of init          </script>      And my mapfile (some url and file locations changed for privacy, but not necessary for this question)               MAP       NAME "testgrib"       IMAGETYPE      PNG       EXTENT         -14000000 3000000 -7000000 7000000       STATUS         ON       SIZE           2145 1377            SHAPEPATH      "../shapefiles"       SYMBOLSET      "../symbols/symbols35.sym"       FONTSET        "../fonts/fonts.list"       DATAPATTERN "^.*$"       IMAGECOLOR     255 255 255            PROJECTION        "init=epsg:3857"       END       WEB         IMAGEPATH "test/img/tmp/ms_tmp/"         IMAGEURL "http://localhost:8080/wxmap/test/img/tmp/ms_tmp/"         METADATA           "wms_title"                  "WMS Test " ## REQUIRED           "wms_onlineresource"         "http://localhost:8080/cgi-bin/mapserv?" ## Recommended           "wms_srs"                    "ESPG:3857 EPSG:4326 EPSG:4269 EPSG:3978 EPSG:900913" ## Recommended           "wms_abstract"               "This text describes my WMS service." ## Recommended           "wms_enable_request" "*"           "ows_sld_enable" "true"         # testing            "wms_feature_info_mime_type" "text/html"         END       END            OUTPUTFORMAT         NAME "png"         DRIVER GD/PNG         MIMETYPE "image/png"         IMAGEMODE RGBA         EXTENSION "png"         TRANSPARENT ON       END                 LAYER # GRIB attempt         NAME         mosaic         STATUS       ON         TYPE         RASTER              DATA         ../grib/ds.vis.bin              CLASSITEM "[pixel]"         METADATA          "gml_include_items" "all"          "wms_include_items" "all"         END          # testing          DUMP TRUE         HEADER  "../templates/test_header.html"         TEMPLATE "../templates/test_body.html"         FOOTER "../templates/test_footer.html"              PROJECTION          "proj=lcc"          "lat_1=25"           "lat_2=25"          "lat_0=0"           "lon_0=-95"           "x_0=0"          "y_0=0"          "a=6371200"          "es=0.0"          "+no_defs"         END              PROCESSING "NODATA=32129"         PROCESSING "BANDS=01"         PROCESSING "SCALE=0,11200"         LABELITEM "[pixel]"              CLASS         NAME "< 1/4"         EXPRESSION ([pixel] < 400 )         STYLE             COLOR 100 0 100         END         END         CLASS         NAME "1/4 - 1/2"         EXPRESSION ([pixel] >= 400 AND [pixel] < 800 )         STYLE             COLOR 200 200 0         END         END         CLASS         NAME "1/2 - 3/4"         EXPRESSION ([pixel] >= 800 AND [pixel] < 1200 )         STYLE             COLOR 150 0 0         END         END         CLASS         NAME "3/4 - 1"         EXPRESSION ([pixel] >= 1200 AND [pixel] < 1600 )         STYLE             COLOR 200 0 200         END         END         CLASS         NAME "1 - 2"         EXPRESSION ([pixel] >= 1600 AND [pixel] < 3200 )         STYLE             COLOR 100 100 0         END         END         CLASS         NAME "2 - 3"         EXPRESSION ([pixel] >= 3200 AND [pixel] < 4800 )         STYLE             COLOR 255 153 0         END         END         CLASS         NAME "3 - 5"         EXPRESSION ([pixel] >= 4800 AND [pixel] < 8000 )         STYLE             COLOR 0 0 200         END         END         CLASS         NAME "5 - 6"         EXPRESSION ([pixel] >= 8000 AND [pixel] < 9600 )         STYLE             COLOR 0 100 0         END         END         CLASS         NAME "6 - 7"         EXPRESSION ([pixel]  >= 9600 AND [pixel] < 11200 )         STYLE             COLOR 0 200 0         END         END         CLASS         NAME "> 7"         EXPRESSION ([pixel] >= 11200 )         STYLE             COLOR 250 250 250         END         END            END # End of mosaic layer          END       And I need template files here for the response: test_header.html               <!-- MapServer Template -->     <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"         "http://www.w3.org/TR/html4/transitional.dtd">     <html>       <head>         <title>MapServer Template Sample</title>       </head>      test_body.html               <body>       band1: [value_0]     </body>      test_footer.html               </html>      I'm not positive that the test_body.html is correct. But, the main issue is that when I click on the image, nothing happens. Firebug doesn't indicate any query is sent to MapServer for the layer. Inserted a breakpoint in the WMSGetFeatureInfo code at the getfeatureinfo function but the execution didn't go into the code. As I realize that perhaps the issue centered around not declaring the layer as OpenLayers.Layer.WMS instead of MapServer, I also attempted to create the layer using OpenLayers.Layer.WMS. But, that created issues of its own. Kind of a toss up then which question to ask. Might ask the other one later. So, if anyone has any thoughts on how to use this with the Mapserver layer or perhaps a different control, I'd really appreciate the help. Thanks!


How to join external tables with a shapefile's attribute table in QGIS?

I need to join a .csv table into an attribute table of a shapefile. I wanted to do this with QGIS and i don't seem to find anyway of doing it. I know that the mmqgis plugin does it, but my shapefile has 6000 records and this plugin crashes at 2000 and something... Any ideas?


Running database operations against a stored GIS geometry

I'm using Quantum GIS and Postgres with the GIS extensions to store shapefiles that I'd like to run point in geometry and point in union geometries against. While it's easy to save and load the vector layers to the database, I'm having problems figuring out how I can run operations against them. How do I use the database to SELECT WHERE lat/long values in another table are contained with a saved polygon geometry? How do I union a couple of geometries and use that in the WHERE clause? I'm having a lot of trouble figuring out what the schema of the saved polygons actually is. Edit: The shapes have a number of holes in them, if it matters for performing the selects and the unions.


how to create a new postgis database for usage with QGIS under windows

After installing postgresql with postgis under windows xp I have 3 databases:               postgis     postgres     template_postgis      now I try to create a new spatially enabled database in pgAdminIII: I select template_postgis in the pgAdminIII object browser and via rightclick in the context-menu I select `new object ->new database`. when I try to use this database in QGIS with the spit-plugin I get the following error message: > the chosen database does not have PostGIS installed, but this is required > for storage of spatial data. How do I create a working spatial postgis database? please help a rookie thanks Kurt


Why is ArcGIS crashing when I try to load .adf raster data?

I am pretty new to GIS, so forgive me if this is a stupid question. I downloaded a bunch of land cover files from USGS seamless warehouse of Florida. Every time I try to 1) view the .adf file (which I am assuming is the raster data image) in arcCatalog or 2) add the .adf file in arcMap the program crashes. Any reason for this? I am running the trial version of ArcMap 10 on my home PC in order to avoid using the slow ones at my universities lab, so I figured maybe I am missing some kind of extension or something. I am not really sure, I do not have much experience working with raster data so far.


Joining hdf files with multi bands and no attribute tables to a shape file

I have a hdf file showing raster data and I want to join this to a shape file (in Acrgis) so can aggregate the values to areas rather than cells. Unfornuately the hdf has no attribute table so I cannot join. When I try to create one it says it can't be done as isn't a single band file. Is there a way around it so I can join the data? Thanks!


Metadata for paper maps

In GeoNetwork/ISO 19115, how do you describe physical access parameters for a physical paper map (i.e. of type CI_PresentationFormCode==mapHardcopy)? For example I want to be able to find the map through a room number, archive box and signature number associated to the map. Also interesting would be to store the size (A0, A1, ...) or weight. Storing online resources works quite well with GeoNetwork and is good documented, for example you can enter a online location (URL) via CI_OnlineResource, But what about physical parameters of "hardware" maps? Are there other standards than ISO19115 that are better suited for physical maps? **Clarification:** I'm aware that in ISO-19115, one can define that a data **IS** a paper map by setting CI_Citation.presentationForm to mapHardcopy (an element from the CodeList CI_PresentationFormCode). But: If my data is a paper map, how do I describe properties of it that are very specific to paper maps:   * weight   * size   * book shelf   * material health   * library number The only thing I can imagine is the library number which is storable in CI_Citation.identifier. What about the others?


Add Zoom Level to Map within TileLayer

I used ArcGIS javascript API 3.0 version and I use Bing map. I need to be able to zoom to 22., 23., 24., .. levels. So I was added an addNumLevels function to Map;                esri.Map.prototype.addNumLevels =function(num){              for(var i = 0; i < num; i++ ){             var ln = this.__tileInfo.lods.length;             this.__tileInfo.lods[ln] = {};             this.__tileInfo.lods[ln].level = ln;             this.__tileInfo.lods[ln].resolution = this.__tileInfo.lods[ln-1].resolution / 2;             this.__tileInfo.lods[ln].startTileRow = this.__tileInfo.lods[ln-1].startTileRow;             this.__tileInfo.lods[ln].startTileCol = this.__tileInfo.lods[ln-1].startTileCol;             this.__tileInfo.lods[ln].endTileRow = this.__tileInfo.lods[ln-1].endTileRow * 2;             this.__tileInfo.lods[ln].endTileCol = this.__tileInfo.lods[ln-1].endTileCol * 2;         }     }      But I when I upgraded the ArcGIS js api to version 3.2, my addNumLevels function stopped working. Can anybody help me? Thanks.


arcgis10: create perpendicular transects to stream at specified intervals

I am trying to do what my title says. The goal is to input a line, segment it in an evenly distributed way by a specified distance, convert the vertices to points, then draw lines perpendicular to the segment at a specified distance. These lines will then be converted to point in a similar fashion and return the maximum elevation value both east and west of each point. This whole process is meant to extract samples of an approximation of bank height for a river. I have found tools and processes that do the line segmentation (even though annoyingly it seems one cannot access the nifty Editor toolbar processes in python), but I am having trouble finding a script that makes the perpendicular transect lines. Has anyone done anything like this before? Note: I have found the script here, but it has too many malfunctions and I gave up trying to fix it (http://forums.arcgis.com/threads/49206-Perpendicular- transects-at-regular-intervals). Thank you!


Is there an se_toolkit version for ArcSDE 10?

I have the one for 9.3 downloaded--is there a version for ArcSDE 10? So far my 9.3 version of se_toolkit is working against 10 geodatabases, but I've only tried a couple of the tools(sdeping and sdestats).


Edit End Points in Polyline (Python/ArcMap10)

I am trying to make a polyline endpoint snapping tool in python. using a Search Cursor I read in the values to an array of line objects (which contain the feature id of the line, and two "EndPoint" objects - first and last - which contain the X and Y coord of each), then process some data interpolation of the points to give all points their correct values. So far, that part works perfectly. However, when I try to go back and edit the shapefile, I can see the data that I am changing, but I can't seem to edit it. Here is the code for updating the fields:               # --------------------------     #     # Update the file with the new EPs     #     # --------------------------          # Create update cursor     #     rows = arcpy.UpdateCursor(outputDirectory + "\\" + trails_fc + ".shp")          # Enter for loop for each feature/row     #     i = 0     for row in rows:         # Create the geometry object         #         feat = row.getValue(shapefieldname)              partnum = 0              if feat.getPart(partnum)[0].X != allLines[i].startEP.x:             arcpy.AddMessage("Change: " + str(feat.getPart(partnum)[0].X) + " to " + str(allLines[i].startEP.x) )             feat.getPart(partnum)[0].X = allLines[i].startEP.x              rows.updateRow(row)              arcpy.AddMessage("Now is: " + str(feat.getPart(partnum)[0].X))              i+=1      My readout consists of statements like this: Change: -105.512166832 to -105.699533165 Now is: -105.512166832 For some reason, the rows are not updating. And for the life of me, I can't find a tutorial or instructions on how to edit a particular point in a polyline. I can only find how to edit a point as a field in a point shapefile. Anyone have any ideas? Thanks a lot!


Arcpy cursors, WHERE clauses, and date/time fields

I'm having trouble making a selection of rows that meet certain date criteria using my arcpy.UpdateCursor(). I need to update records in a field in a geodatabase table with an ascending value based on the date field. Esentially, what I want is for each date on or after 6/1/2004 and before or on 9/30/2004 to have a value starting at 1 ascending to 122 (the number of days between those dates, inclusively). The trouble seems to be with the syntax in my WHERE clause, I get a syntax error that says: > Syntax error (missing operator) in query expression '( "VAL_DATE" >= > date'2004-06-01 00:00:00' AND "VAL_DATE" <= date'2004-09-30 00:00:00' )'. Can anyone point me to the "missing operator" the error is throwing in my face?? FYI, here is the bit of code of I'm using:               import arcpy     from datetime import datetime          table = r'removed...\noaa_precip_master.mdb\noaa_precip_master'          start_sim_date = datetime.strptime('06/01/2004', '%m/%d/%Y')     end_sim_date = datetime.strptime('09/30/2004', '%m/%d/%Y')          sim_date_counter = 1          where = '"VAL_DATE" >= date\'%s\' AND "VAL_DATE" <= date\'%s\'' % (start_sim_date,end_sim_date)     ## where = '"VAL_DATE" > date\'%s\'' % (start_sim_date)     print where     rows = arcpy.UpdateCursor(table,where)     for row in rows:         date = row.VAL_DATE         print date, sim_date_counter         row.SIM_DATE = sim_date_counter         rows.updateRow(row)         sim_date_counter += 1     del rows      UPDATE: Solved, see change to "WHERE" statement below:               where = '[VAL_DATE] >= #%s# AND [VAL_DATE] <= #%s#' % (start_sim_date,end_sim_date)


apache vs lighttpd vs nginX as OpenLayers / GIS Server

How is using one of the above technologies to host an Openlayers-based GIS application? My investigation with apache and lighttpd showed that apache is not the best for heavy load environments. In some cases there have been suggestions to try out the ngin-X route. Has anyone tried the three against each other and made a benchmark or a test? Thanks!


ArcGIS 10 IGPFunction2 and IGPFunctionFactory Implementation

I'm creating a custom ArcToolbox tool in C# .NET in ArcGIS 10 SP2. Though I seem to be having two issues which maybe related: Firstly, my tool is run as a background task i.e. the GP dialog does not show, instead the progress bar displays on ArcMap's status bar. Can I get the GP progress pop up dialog to display instead when my GP tool is executed? Secondly, I want to use IGPUtilities3:AddToMap or IGPUtilities3:GetMap() in the execute method of the GPFunction class but neither work. Either my dataset doesn't get added or GetMap() return null. If I check the InApplication or CanAddToMap properties they return false. I'm just doing a IGPUtilities3 utils = new GPUtilitiesClass. I guess it hasn't tied into the ArcMap application but how do I fix this in a GPFunction class? Thanks Rob


Will mousedown events fire with Touch events on OpenLayers mobile?

I can't get the mousedown event to fire on a mobile device when I click on a marker. I'm using openlayers 2.11 with TouchNavigation. The mousedown event fires just fine in Chrome on my Desktop. Are there touch specific events for openlayers? If yes, where can Ifind the list. The openlayers documentation is driving me crazy. The code I use:                 var start_point = new OpenLayers.LonLat(latlngs[0].x,latlngs[0].y);       var start_marker = new OpenLayers.Marker(start_point,green_pin);       start_marker.events.register('mousedown', start_marker, function(evt) { alert('Ok'); OpenLayers.Event.stop(evt); });       StartFinishLayer.addMarker(start_marker);


Accessing "many to many" relationship class using ArcGIS Identify?

I am trying to create a many to many relationship class in ArcGIS 10.0. When Identifying the features though, I don´t see the associated table values. Here is the Class schema: Destination Table name = "TypeDef" ( a table containing a list of Type definitions)               TypeID* | Typ_Name      Origin Table name = "Points" ( A point geometry)               ElementID* | Element_Name      Association Table name ="ElementTypes" (join table with ElementIDs and TypeIDs)               ElementID* | TypeID*      My ManyToMany Relationship class is defined as:               SourceTable     PK: ElementID     FK: ElementID          DestinationTable     PK: TypID     FK: TypeID      However, when I use the Identify Tool on a point, only the Relationship name is shown and not the attributes. What am I doing wrong? ![Many2Many](http://i.stack.imgur.com/NM2WY.png)


Why is local hosting of ArcGIS Online webmap giving error?

I have been trying to locally host this webmap I created using AGOL (to modify it): http://dpwh.maps.arcgis.com/apps/MapTour/index.html?appid=7bd06d8330254b2688577d82f81bb92e It uses the Map Tour template. I did step 1 of the readme file, and since I already had built the application inAGOL, skipped to step 5 and configured configOptions to use the appid, such as it now looks like this               configOptions = {                     // Enter either a web map ID or a web application ID                     // - Using a Map Tour application ID created through the interactive builder in ArcGIS.com allow to reuse                      //   the settings like the header configuration without having to edit any configuration here                     // - Web map and application can be dynamically specified through url parameters (webmap and appid)                      //   if the authorizedOwners property is specified (see below)                     // - If an application id is specified, the webmap configuration isn't used                     webmap: "",                     appid: "7bd06d8330254b2688577d82f81bb92e",                     // Optional list of authorized application or web map owners                     // Has to be specified when using appid or webmap url parameters (e.g. ["user1"], ["user1", "user2"])                     authorizedOwners: [""],                     // Select the application layout between "three-panel" and "integrated"                     layout: "three-panel",                     // The application header title, if no title is specified, the webmap's title is used                     title: "",                     // The application header subtitle, if not specified the ArcGIS.com web map's summary is used                     subtitle: "",                     // Optionally force a web map layer to be used as the tour points data source (can be overridden through URL)                     // Otherwise the upper visible point layer is used (excluding map notes)                     sourceLayerTitle: "",                     // Use the first data record as an introduction instead of a point (can be overridden through URL)                     firstRecordAsIntro: true,                     // Optional zoom level to be applied for the story points following introduction (can be overridden through URL) (e.g: 0 to 16 for Gray Canvas, 0 to 19 for Street Map, -1 to disable)                     zoomLevel: -1                 };      Loading the local version however only produces a "Fatal error: Failed to get portal configuration" error. This is the case even if I only use the webmap id, or just load the page unmodified. Help please to make this work?


How can i keep browser from caching map services?

I am using esri's js api with arcserver 10.1 map services. My issue is that sometimes another application deletes a feature from the geodatabase that is the basis for my map service. But my browser does not react quickly to the change. It seems to be caching the map service response. This is evident in ie and firefox. Even a browser refresh will still show the feature that should have been deleted already. Is there a way to force the cache to flush and get new map service data?


How to Calculate Geolocation Intersection of Two Vectors?

I'm unfamiliar with much GIS terminology but am tasked with solving a specific geolocation problem. For what it is worth, this is for an iOS project, so the data I can collect for analysis is that which is provided by iOS 7.1. I'm a little out of my depth here, but I feel the answer must be trivial. I have two GPS devices, A and B. Let us say that both have broadcast a start location and a current location. All I have are latitude and longitude readings, nothing fancy, and I cannot rely on any heading data I have available. If I know A-start and A-current, as well as B-start and B-current, I should be able to extrapolate vectors from those coordinates and determine if those two devices' paths will cross and at what geolocation. My understanding on this, given that my base use-case is in-city, is that I can consider the world to be flat for this purpose. As well, I do not need a super-accurate answer, just a good, fast, best-guess. Presumably I would need to convert these coordinates to a standard XY cartesian system, calculate, and convert the answer back to real-world geolocation? For that matter, any resources that collect these types of calculations would be appreciated. As I am unfamiliar with the problem space, I'm having trouble getting good Google results in my searches. As well, I am open to any libraries that may be recommended for such data analysis.


Difference between python IDLE and eclipse on accessing GRASS

I have seen these posts:GRASS and Python, Using Eclipse to develop GRASS Python programs. I used to use python IDLE to control GRASS, but I can't make python and GRASS connected, as I described in this post:Problem with python script to control GRASS GIS from outside - How to import grass.script under Win 8.1? So now I am thinking of using eclipse to control GRASS, trying to avoid problems I faced as I used python IDLE. But I don't know whether I will get the same problem after I use eclipse, and I can't tell the difference between using python IDLE alone and using eclipse. I am working with GRASS 6.4.3 and python 2.7.6 under WIN 8.1. Please help me to make decision and furthermore, to solve the connecting problem, thank you!


How to get higher resolution images from USGS into openlayers

Is there a way to get better resolution maps from USGS? Thanks,


How do GIS web apps work?

I am completely new to GIS and therefore have questions regarding how the GIS Web development should be made, because the entire concept still confuses me. So, imagine I have a map. And by this, i mean the shape files for the layers and so on. How i should i use this in my web application? I have read about MapServer, OpenLayers, PostGIS, but can't quite understand how and why this all enters the question. If i have my map, my files, shouldn't the problem be just on to represent it in the web page? As you probably have already realized, i'm quite clueless on this, and if you could explain this to me like i'm really dumb, i would be appreciated :D


Transit Routing using pgrouting as core routing engine?

i'm a student and i'm working on a project. the aim is given two sets of coordinates(origin & destination) route through the map to determine the optimum path using public transportation and of course walking edges. i read about GTFS, but its too complex,i dont care about transit schedule, time, etc. my only aim is to provide the user a graphical representation of the shortest path. and the cost will only be the length of edges multiplied by the fare of mode of public transportation (e.g. length in kilometer * bus fare per kilometer), edge length multiplied by the average speed of the mode of transportation(e.g. legnth in kilometer * bus average speed by kilometer). besides, my country hasn't yet provide any gtfs data yet. i've decide to use postgresql, postgis and pgrouting to do this. Here's my idea and what i've done so far. first: download an openstreet map extract taken from http://metro.teczno.com/#manila. then: convert the XML OSM data to a routable database using osm2pgrouting tool. I've done this using cgywin since i'm using windows. encountered some error though, i think?. see my results below: > reference nd=[some digits] has no coressponding node entry..   > ..   >  Nodes table created   >  2create ways failed:   >  Types table created   >  Way_tag table created   >  Relations table created   >  Relation_ways table created   >  Classes table created   >  Adding tag types and classes to database...   >  Adding relations to database...   >  Adding nodes to database...   >  Adding ways to database...   >  Creating topology... NOTICE: CREATE TABLE will create implicit sequence > "vertices_tmp_id_seq" for serial column "vertices_tmp.id"   >  CONTEXT: SQL statement "CREATE TABLE vertices_tmp (id serial)"   >  PL/pgSQL function "assign_vertex_id" line 14 at EXECUTE statement   >  '#########################   >  size of streets: 1072190   >  size of splitted ways : 2281124   >  finished dont know exactly if this kind of result should concern me but i tried to do some routing. seems working so far. if you have any comments on my osm2pgrouting result, i would appreciate the help. =) anyway after converting the XML data, load QGIS application and load the ways table using DB Manager PostGIS plugin to convert the table to a shapefile, i will then enable labeling to see the edge id of each edge. then i will start tracking which edges are used for public transportation, (given that i already have a physical map which contain routes) then: i will create a new table or multiple tables that will establish relationship between the "way" table and my public transportation datas. i haven't decided their structure yet, but thats the idea. now using pgrouting. im going to execute a shortest path query which will only route on edges that are also transit edges and walkable by foot (footway). now my question is, Am I on the right track here? Am i missing something or there is something wrong with my idea? please i would really appreciate suggestions, corrections etc. and if there are things that would help me in making light to this project of mine. Thanks.


Python error when running SpatiaLite with OSGB(27700)

I have a project open in QGIS and am trying to import an Ordnance Survey Address Point layer into SpatiaLite. I have recently installed the Qspatialite plugin. I open SpatiaLite in the DB Manager then use the Import QGIS layer to choose the Address Point layer. I then get the following Python error code data and have no idea of what is happening. An error has occured while executing Python code: Traceback (most recent call last): File "/Users/Nlaw/.qgis//python/plugins/QspatiaLite/importQgis/importQgis.py", line 64, in run if not self.parent.db.uploadQgisVectorLayer(layer, tableName=self.tablenames[i], srid=self.srid, selected=self.selected): File "/Users/Nlaw/.qgis//python/plugins/QspatiaLite/Classes/Database.py", line 249, in uploadQgisVectorLayer header,data=self.executeQuery("""INSERT INTO "%s" VALUES (%s,%s)"""%(tableName,','.join([unicode(value).encode('utf-8') for value in values_auto]),','.join('?'*len(values_perso))),tuple([unicode(value) for value in values_perso]),commit=True) File "/Users/Nlaw/.qgis//python/plugins/QspatiaLite/Classes/Database.py", line 45, in executeQuery cursor.execute(query,params) IntegrityError: AddressPoint.Geometry violates Geometry constraint [geom-type or SRID not allowed] Python version: 2.7.2 (default, Jun 20 2012, 16:23:33) [GCC 4.2.1 Compatible Apple Clang 4.0 (tags/Apple/clang-418.0.60)] QGIS version: 1.8.0-Lisboa Lisboa, Python path: ['/Applications/QGIS.app/Contents/MacOS/../Resources/python', '/Users/Nlaw/.qgis//python', '/Users/Nlaw/.qgis//python/plugins', '/Applications/QGIS.app/Contents/MacOS/../Resources/python/plugins', '/Library/Frameworks/SQLite3.framework/Versions/3/Python/2.7', '/Library/Frameworks/GDAL.framework/Versions/1.9/Python/2.7/site-packages', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python27.zip', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat- darwin', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat- mac', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat- mac/lib-scriptpackages', '/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib- tk', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib- old', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib- dynload', '/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/PyObjC', '/Library/Python/2.7/site-packages', '/Applications/QGIS.app/Contents/Resources/python/plugins/fTools/tools', '/Users/Nlaw']


How do you explain what GIS is?

Whenever I do a presentation about anything related to GIS, I consider the level of knowledge of GIS of the audience. The worst times are definitely when I need to explain GIS to a group that has almost no knowledge about it. I said virtually, because by know everyone has interacted with some type of electronic map in one way or another (elephant in the room is Google Maps), so I usually start from there and explain everything that is necessary to create Google Maps. I have found that works quite well. I know there are several formal definitions of GIS. The questions are:   * Which formal definition of GIS do you use?   * How do you explain it to somebody that doesnt know much about the field in elevator pitch style? (30 seconds)   * How do you explain it when you do have more time? (10 min)


Viewing panoramic imagery in ArcGIS

Is there an ArcGIS extension or utility that will view 360 degree panoramic images? I've found utilities that are part of larger software packages for managing mobile scanning data, but all I'm after is the spherical viewing capability. If there isn't a tool out there, does anyone know of a .NET control that I could incorporate into an ArcMap extension that will handle spherical display?


Minimum field in attribute table

I have two fields in an attribute table and I want to add a third one which will contain the minimum of these two fields. For instance, Field 1 Field 2 Field 3 100 50 50 20 50 20


Union with Concatenated Values

I have a polygon feature class in which the polygons frequently overlap each other. I want to find all of the unique combinations of overlapping polygons within that feature class. I believe I want to perform a Union on the feature class, but I also want to be able to concatenate the names so that each polygon has a unique name made up of the concatenated name of all the polygons that make it up. Any ideas on how to do this?


What is a good Javascript web API (preferably open) for building web maps that require little/no maintenance?

I have recently volunteered to make a web map for a local non-profit. The map will be pretty basic. It needs to display Lane county, Oregon split into multiple districts, with a "central office" location for each district. I plan on building this in my spare time over the summer, and I want to be able to hand this off to them with very little risk of it needing any sort of future maintenance. I have built a number of web maps using the ArcGIS Javascript API with ArcGIS Server, but I am planning on learning a different API since this organization doesn't have any ESRI software (or any GIS for that matter). I have been considering OpenLayers/OpenGeo. One thing to keep in mind is that this organization doesn't/won't have a server set up to host the layers of this map so I'm not sure how I would store and access the data. JSON maybe? I'm currently unfamiliar with most non-ArcGIS solutions, so any advice would be appreciated.


Feature Layer with ONDEMAND mode

I can't understand how to display all features from the feature layer. In the code bellow I specified 2 features in the layer and set in feature layer's mode to esri.layers.FeatureLayer.MODE_SNAPSHOT. According esri documentation this options means that all features will be selected and displayed. But they are not. The code is:                  var jsonFS = {                     "displayFieldName": "Name",                     "fieldAliases": {                         "Name": "Name"                     },                     "geometryType": "esriGeometryPoint",                     "spatialReference": {                         "wkid": 102100 //WGS_1984_Web_Mercator_Auxiliary_Sphere                     },                     "fields": [{                         "name": "Name",                         "type": "esriFieldTypeString",                         "alias": "Name",                         "length": 255                     }],                     "features": [{                         "attributes": {                             "Name": "bla"                         },                         "geometry": {                             "x": -8919439.31450887,                             "y": 4928270.761925456                         }                     },                     {                         "attributes": {                             "Name": "bla bla"                         },                         "geometry": {                             "x": -8155495.379532158,                             "y": 5075380.311392084                         }                     }]                 };                 var fs = new esri.tasks.FeatureSet(jsonFS);                      var featureCollection = {                     layerDefinition: {                       "geometryType": "esriGeometryPoint",                       "fields": [                       {                         "name": "Name",                         "type": "esriFieldTypeString",                         "alias": "Name"                       }                       ]                     },                     featureSet: fs                 };                  var featureLayer = new esri.layers.FeatureLayer(featureCollection, {                     mode: esri.layers.FeatureLayer.MODE_SNAPSHOT                 });


ArcObjects - Get units from ISpatialReference

How do I get the units (meter/feet) from an ISpatialReference object in ArcObjects?


Change output file type in QGIS georeferencer?

I have successfully georeferenced a raster file via the Georeferencer plugin (Georeferencer GDAL Version 3.1.9) in QGIS 1.7.0. However, I noticed that the output georeferenced raster file is in TIFF format. Is there a way for me to change it to something like JPEG (or others) to save disk space and memory usage? Thanks. BTW, I am running this in Mac OS X 10.6.7.


catchment area using SAGA and Arcgis 10.2

I would like to calculate catchment area. I used D8 algorithm and running in both arcgis 10.2 and SAGA GIS 2.1.2. However, the results from two software are different. The lowest value of catchment area in arcgis = 0; high value = 109325, while SAGA GIS lowest value= 625, highest value=47820. I think the result should be similar and maybe I do something wrong. Could you please explain what I did wrong? These are what i am doing, the input: DEM25 In arcgis, to calculate catchment area: Dem25, then run tool flow direction, then run tool flow accumulation, then 25 (cell size) * 25 * flow accumulation In SAGA GIS, terrain analysis -- hydrology --catchment area (parallel), elevation = dem25, sink routes=no set, weight = not set, catchment area = create, Step =1, linear flow = no tick, linear flow threshold =500. convergence = 1.100000001 Thanks so much for your help.


OpenLayers Vector Layer not aligned with Google Oblique Imagery

Check http://jsfiddle.net/vasile/BKvcX/6/ In the first map, the yellow line should follow the street curve, however if you pan the map up-down, you clearly see the misalignment due to what ? Different projection used by the oblique imagery ? The same line, on a normal Google Maps integration (second map) is rendered "normal" How can the misalignment in the first map be fixed ?


Displaying "numberoffeatures" in grid

How to display numberoffeatures in map using _geoext_ on grid?                Ext.onReady(function() {                             app = new Ext.Viewport({                     layout: "border",                     items: items                 });               items.push({             xtype: "grid",             ref: "capsGrid", // makes the grid available as app.capsGrid             title: "Available Layers",             region: "north",             height: 150,             viewConfig: {forceFit: true},             store: new GeoExt.data.WFSCapabilitiesStore({                     url: "http://localhost:8080/geoserver/wfs?request=GetFeature&typename=GisData:sites&version=1.1.1&resultType=hits&filter=(%3CFilter%3E%3CPropertyIsEqualTo%3E%3CPropertyName%3EVendor%3C/PropertyName%3E%3CLiteral%3EvendorTE%3C/Literal%3E%3C/PropertyIsEqualTo%3E%3C/Filter%3E)",                     protocolOptions: {version: "1.1.1"},                     autoLoad: true                 }),             columns: [                 {header: "numberOfFeatures", dataIndex: "numberOfFeatures", sortable: true}             ]         });      I get xml response in firebug: ![enter image description here](http://i.stack.imgur.com/03baY.png)


ArcMap field calculator VBScript

I want to do some advanced level scripting with VBA in the field calculator. I have done a simple coding in the code block. Here is the codes.               dim X       if [Name]="University" then         X=101101       elseif [Name]= "College" then         X=101102       elseif [Name]= "Medical College" then         X=101103       elseif [Name]= "School" then         X=101104       else         X=""       end if      This is working correctly. But I want to do something more. Those are-   1. If the Name field contain the word University or College or Medical College or School then X return the expecting value.   2. Is it possible to create new fields through VBScript? If its possible then the X will split into two from its middle (3-3) and put the value in the newly generated column/fields. If this possible it will help me a lot. Thanks in advance


Getting lane data from OSM

I want to extract road network data from OSM for a small city. From where can I get the lane information of any road in OSM. I tried exporting the data in the OSM XML format and there was no lane fields in that file.


Getting number of tables loaded with integrated mapping

I am integrating a MapInfo map into a form in an Access database. I am trying to find out the number of tables loaded but NumTables() and MapperInfo() are both getting "unrecognized command" errors. Is there another way to get the number of tables?


Point to Polygon

> **Possible Duplicate:**   >  Create new polygon layers from postcodes latitude & longitude I'm running QGIS 1.8.0 I have a list of points that I need to convert to one polygon. So far I have been using the **Point2One** plugin which produces a zigzagged polygon, instead of one full polygon. To change the output result, I modified the numerical order of the **ID** column in the dbf to connect points in order yet I still get the same result. The "Create output features based on input field" or "Sort points by this field" yields an error I've looked at Create a polygon layer from point data but my problem seems a bit different. I thought converting these points to a csv or txt format might work. This seems like fairly easy problem to solve but I'm hung up. Any suggestions? Thanks all!


Get Point at distance on a LineString using NetTopologySuite

I am quite sure there is a function in NTS that actually gets us a point at a distance along a line starting at first vertex. I am using NetTopologySuite version 1.12.1. In Shapely, the python version of JTS, it is quite simple to use               >>> LineString([(0, 0), (0, 1), (1, 1)]).interpolate(0.75, normalized=True).wkt     'POINT (0.5000000000000000 1.0000000000000000)'      I am trying to do the same thing using NTS but i hightly doubt they could be using a jargon word for this. I tried LinearReferencing namespace, but could not find a suitable method. Does anyone have idea on this?


Relative Hyperlinks in QGIS without end user's software dependency

I am trying to work out how can I add a relative path hyperlink to a shapefile in QGIS. I have found several great tips on various forums but only one has so far potential to be implemented. The overall goal is for the end user to be able to access a hyperlinked pdf/png/jpg etc without me having to set up any software dependencies under layer properties Actions tab in QGIS. I assume that the end user will have the software required to open a document but I would like the end user's computer to choose a default program to open the files. One of the possibilities would be to use a 'Python' action which I located on one of the forums which after several modifications did not work for me: http://lists.osgeo.org/pipermail/qgis-user/2011-May/012117.html. Here is the screenshot of my layer's action properties: ![enter image description here](http://i.stack.imgur.com/1xF4F.png) In the attribute table, the column with hyperlink is called 'PDF and the relative hyperlink is: /doc/document.pdf My folder structure is:                - /HyperlinksTest      -- /data      --- Renewable_Energy_Zones.shp      -- /doc      --- document.pdf      HyperlinksTesting.qgis      When I identify the feature, I get the following Python error: ![enter image description here](http://i.stack.imgur.com/odkDG.png) Am I using wrong variables in the Python code? Is there any other way of using relative paths in hyperlinking data to external documents? I noticed there is eVis extension but as far as I see it relies on specyfying software to open a document with, which I'd like to avoid. Many thanks, Magda


PostgreSQL query selection in MapServer is not working?

When I make a selection in the query in Mapserver I got a blank page. The query I use is: DATA "geom from (select * from bagmutlijst_201210.vo_gepland_gemeente where gemcode = '0141') as selection using unique gemcode using srid=28992" The data is correct because if you use the following query (without a selection) I got a correct map. DATA 'geom from bagmutlijst_201210.vo_gepland_gemeente using unique gemcode using srid=28992' Why is the selection not working?


OpenLayers.Layer.Vector not working?

I'm struggling to get my vector layer working. The wfs works as I have tested it by creating a wfs layer and adding it to the map, now I'm trying to display a vector layer with the results of a search but nothing is showing up. The results of the search appear in the gridpanel but not in the vector layer. Here is the code for the search:               var features = [];     var wfs = new OpenLayers.Layer.Vector("Search");     var select = new GeoExt.grid.FeatureSelectionModel();          var protocol = new OpenLayers.Protocol.WFS({         url:  "http://localhost:8080/geoserver/wfs",     featureType: "location_points",     featureNS: "http://www.openplans.org/topp",     geometryName: "the_geom"     });               formPanel = new GeoExt.form.FormPanel({     title: "Species Search",     height: 150,     region: "north",     protocol: protocol,     items: [{         xtype: "textfield",         width: 200,         name: "scientific__like",         fieldLabel: "name",         allowBlank: false,         minLength: 4     }],     listeners: {         actioncomplete: function(form, action) {             features = action.response.features;             store.loadData(features);             vm=map.getLayersByName("Search");             if(vm.length===0){                 wfs = new OpenLayers.Layer.Vector("Search");                 map.addLayers(wfs);                 store.bind(wfs);                 select.bind(wfs);             }         }     },     buttons: [{text: 'search',         handler: function(){             formPanel.search();         }     }],     keys: [{ key: [Ext.EventObject.ENTER],          handler: function() {             formPanel.search();         }     }]     });      Firebug shows the following error: layer.map is null layer.map.addControl(selectControl); (FeatureSelectionModel.js (line 174)) I followed an example on line (http://ian01.geog.psu.edu/geoserver_docs/apps/gaz/index.html). I appologise if someone has already asked this, I've looked everywhere for the answer and found nothing, I'm quite new to this, so any help would be very much appreciated. Thanks Will


buffer using topological operator

IMap map = axMapControl1.ActiveView.FocusMap;             IFeatureSelection fSelection = map.get_Layer(0) as IFeatureSelection;                  ISelectionSet selSet = fSelection.SelectionSet;                       ICursor cursor;             selSet.Search(null, true, out cursor);                  IFeature feat = ((IFeatureCursor)cursor).NextFeature();      Now i want to create buffer to my cursor. How can i do it?


How to map a range of values to a single value in FME

A dataset with a set of points has an attribute called HEIGHTS which is filled with numbers like, 10.8, 0.4 and 150.6 etc... I want to transform this dataset to a shape file. In the transformation I want to add a new column to the shape file which will need to contain a color code.   How can I map the height values to this color code?   For example all heights between 5-10 should have the color dark blue, all heights between 10 and 20 should have the color light blue etc...... I tried to use the valuemapper transformer, but it only seems to work for 1:1 mappings. I am using FME 2009.   This means that I cannot use the AttributeRangeMapper and AttributerRangeFilter transformers.   Can I conclude that what I want is not possible in FME 2009 and I should look at a newer version on FME. Thanks for your efforts.