title: Arctic EPSG:4326 map projection?

text: Does anyone know of a source for an Arctic EPSG:4326 map projection? Obviously I'm a GIS newbie, so gentle pointers in the right direction appreciated! It's a long story. I have some arctic ice data that I'd like to visualize with Tableau 8.0 which accepts external wms maps, so long as they're EPSG:4326 (Plate Carree) projections. Obviously such a projection tangent to the equator won't work for high latitudes. But there's nothing that says you couldn't do a mercator-like projection using the prime meridan as the 'equator' so that close to the pole x-y placements would be reasonably accurate. If such a thing existed I could upload it into Tableau 8.0 and plot my points. Or, if there was a recipe for reprojecting an existing map into this format using grass or gdwalp or some such I could do it myself. Anybody ever hear of such a thing? Or have any advice? Thanks.


title: Dynamically set zoom level based on a bounding box

text: I have a leaflet map whose size depends on the size of the browser window. I would like the zoom level to be dynamically chosen so that it is as zoomed in as possible while showing the entirety of the bounding box. Right now, I just have the zoom level hardcoded and the center point based on an average of points.               map = new L.Map('map', {       center: new L.LatLng(         latitudeSum/locations.length,         longitudeSum/locations.length       )       zoom: 9     })      Instead, I'd like to give it a bounding box (two islands) and have the zoom level chosen based on the size of the window.


title: Drop down menu for Coded Domain not available in edit session

text: Using ArcGIS Desktop 10.2.2, Advanced license. I created a new coded domain for a string field. I added the domain to an existing field. When I edit the features, i would expect to get a list of the coded values from the domain as a dropdown menu. This is not happening. What am I missing - is there a simple check mark or radio button setting?


title: GEOExplorer - Get Feature Info doesn't work in IE or Firefox

text: I published a map **here** using GEOExplorer and I can only get the "Get Feature Info" to work in Google Chrome. Nothing happens when I use Internet Explorer or Firefox. Can anyone help fix this?


title: Buffer and Intersect differences

text: “What is the difference between ‘Select by location using a buffer or intersection’ and the ‘Buffer’ or ‘Intersect’ Analysis tools?


title: How do I request an output from a WMS service in a non-standard projection?

text: I want to use a map service which is available in WMS format. The server is Mapserver, and the spatialreference of the map service is 4326. I wish to use it in a slippy map (web based Javascript Map), with a custom projection. This projection does not have a EPSG code. Is it possible to request the output, by providing the well known text for the projection? * * * _Update_ I can see that the GetMap request has a SRS paramter, and usually an EPSG code is supplied, like this: **SRS=EPSG:4326**. Can we pass something else to this SRS Paramter?


title: MapFish Print latitude, longitude innacuracy?

text: I'm new to GIS so this may be trivial: I'm using MapFish Print to create some `PNG`. I have issues to place the viewport. When I use a random GPS coordinate: `-20.904522,55.475106` which seem to point roughly (<20m accuracy) to the same place in google and in openstreetmap.org, my viewport in MapFish Print is way up north (more than 100km). What did I miss ? Here's my config (I've commented some lines that I did try):               {       layout: 'A5 landscape',       title: '',       units: 'dd',       #srs: 'EPSG:900913',       srs: 'EPSG:4326',       #srs: 'EPSG:3857',       #geodetic: true,       #geodetic: false,       outputFilename: 'mapfish-print',       outputFormat: 'pdf',       layers: [         {             type: 'OSM',             format: 'image/png',             baseURL: 'http://tile.openstreetmap.org/',             ## This one seems incorrect and sets the 0,0 several             ## thousand kilometer in the north             #maxExtent: [-180,-90, 180, 90],             maxExtent: [-180,-180, 180, 180],             #srs: 'EPSG:900913',             #geodetic: false,             tileSize: [256, 256],             resolutions : [               1.40625, 0.703125, 0.3515625, 0.17578125, 0.087890625,                0.0439453125, 0.02197265625, 0.010986328125, 0.0054931640625,                0.00274658203125, 0.001373291015625, 0.0006866455078125,                0.00034332275390625, 0.000171661376953125, 0.0000858306884765625,                0.00004291534423828125],             extension: 'png',         },       ],       pages: [         {             center: [55.475106, -20.904522],             #center: [0, 0],             #srs: 'EPSG:4326',             srs: 'EPSG:900913',             #srs: 'EPSG:3857',             geodetic: false,             #geodetic: true,             scale: 500000,             dpi: 190,             mapTitle: "",             comment: "",         }       ]     }


title: Universial addressing systems

text: There are some grid-based global addressing systems. NAC (Natural Area Code) is licensed (expensive), but I've seen one from Japan which used more interesting honeycomb-like divisions. I cannot find reference to the latter one, can anyone help with this?


title: What is the more accepted practice for defining tolerance values?

text: **The Context**   Due to restrictions, I have an application that requires line features to be selected if they are within a user defined tolerance from a points data set. I can not currently use any third party API to perform so have written an efficient algorithm that completes this task effectively. The application will be used by people that may or may not have a sufficient grasp on spatial concepts such as the defined unit of measure for their data set, or parameters requiring 'acceptable' values in map units.   **My Befuddlement**     1. The use of the Tolerance parameter and whether I should keep the tolerance in a commonly known unit such as Centimetres, Meters or Kilometres and perform the required calculations to massage the input into an appropriate value; or    2. Assume the user knows the data, or has the ability to ascertain the data sets map units and enter an intelligent, useful parameter value.    **The Question:**   What is the more accepted practice?   **The Implications**   If _Befuddlement 1_ is the accepted practice or considered more user friendly, can you proved the some direction towards ascertaining the Datum and other relevant coordinate information to perform the required function. If _Befuddlement 2_ is more generally accepted, then case solved. Though it does make the solution easier and more straight forward, I feel forcing the user to go through a discovery process may lessen the user experience of the software. Any thoughts on the matter would be greatly appreciated.


title: Mosaic to New Raster tool fail

text: I'm trying to merge 24 DEM files, ![24 DEM](http://i.stack.imgur.com/QvzEr.jpg)   after applying the tool (Mosaic to New Raster) only 11 out of the 24 merged. ![New Raster](http://i.stack.imgur.com/as0gC.jpg)   All the original files are in the same projection, pixel type, band, etc. I also tried DEM to raster tool first before using Mosaic to New Raster and the results were the same.


title: Georeferencing CAD Causing Close Out/Shut Down of ArcMap?

text: I am currently running ArcGIS 10.2.2 for Desktop. I have no other georeferencing issues regarding TIFFs shapefiles, or otherwise. However, when I introduce a CAD (or DWG) file into the map and attempt to georeference the file with links, ArcMap immediately freezes for a moment and closes out. If I move the CAD file using shift or scale, there are no issues, only when I attempt to use links.


title: How to make feature display after converting Features to Graphics?

text: I used the "Convert Features to Graphics" tool on a polygon layer. It worked correctly with my feature no longer shown but a graphic being displayed instead. However now, I'd like to display my feature again. The feature layer is turned on and the symbology is set, just as it was before. I tried deleting the Graphic group from the Annotation Groups in the Data Frame Properties, and I also tried deleting the graphic, with no success. How do I enable my layer to show on the map again? I'm using ArcMap 10.2. If I add the shapefile again using the Add Data button, it appears in the map. I could use this workaround but it seems like there should be a way to salvage my existing layer, rather than having to abandon a layer each time you convert its features to graphics.


title: How to customize locator file to geocode unusual address?

text: I am very new to geocoding. I am trying to create address locator and geocode a bunch of address using ArcCatalog 10.1, and I get problem for the following addresses:   1. Road Type is in Street Name: eg: Glenview Heights Lane    2. There is direction in street name: eg: East Street   3. I have an address of 75 Romfield Cir (cir represents circuit), somehow it doesn't recongnize it even though I have added cir, circuit as a street type. I am thinking it is probably conflicting with the street type circle (cause they all start with cir).


title: How to load style from postgresql?

text: Is there any way to load a postgresql style for a vector layer using python? There is a function: `loadNamedStyleFromDb (const QString &db, const QString &theURI, QString &qml)`.Does it support loading of postgresql styles?


title: What is the best way to sum the area from a number of polygons in .NET/SQL Server?

text: I have an Mapbasic/.NET application that read data from SQL Server 2008 R2. I have a table with polygons and need to sum up the area of a subset of the polygons. This should work dynamic so the user should be able to sort on different parameters and get the sum of the areas from their selection. My options are:   1. MapBasic   2. C# .NET   3. SQL SERVER Mapbasic seems to be to slow for this, personally I would prefer to do this direct in SQL Server but have not found any way to sum up the area with T-SQL. Have anyone done this before or is there any reason to do this at any of the levels?


title: How to set up a QGIS project with SQL Server which can be opened on another computer?

text: I am using QGIS with SQL server 2008. I need to save project on one computer and then copy project to another 20 clients! After I save project and try to open it on another computer I am getting error. ![HANDLE BAD LAYERS](http://i.stack.imgur.com/RQrzx.png) How I can setup project that can be used on another computer?


title: Problem with Pgrouting 1.5

text: Hi all and thanks in advance. I have Postgis 2.0.1 and PgRouting latest version. I'm using osm2pgrouting to import OpenStreetMap data everything goes ok the trouble is that I always get source and target columns without data and when i run `SELECT assign_vertex_id('ways', 0.00001, 'the_geom', 'gid');` by hand on PGSQL I get the following error: > ERROR: syntax error at or near "vertices_tmp" LINE 1: SELECT > ST_addGeometryColumn(''vertices_tmp'', ''the_geom'', ... ^ > > _**_*** _**Error **_** _*_** > > ERROR: syntax error at or near "vertices_tmp" SQL state: 42601 Character: 31 After hours of reading on the net I've realized that the function assign_vertex_id() had to change somehow because it references other functions like StartPoint that it's called ST_StartPoint and i've changed that but i still get errors. Where can I get the .sql file with the proper SQL to make the upgrade and postgis/pgrouting can work together? Best regards from Portugal, Tiago


title: Difference between Positioning and Localization

text: I would like to know what are the differences between **Positioning** and **Localization**. In most papers they are used interchangeably. **Are they the same? Does one come after another?** I always thought in this order:   1. Localization: GPS coordinates found.   2. Positioning: Found you in a map.   3. Navigation: Now you can navigate thru a map. **Is this correct?**


title: Update table using contains function

text: I have a table that containing placename regions and another table of roads. I have another table that containing centroid of both placename region and roads.The table structure of third table is placename,Roadname,x,y. I want to update the placename field of the roads which are passing through the corresponding placename region. I used the following query:               Update haryana set placename = placename_region.name from placename_region where haryana.the_geom in (select st_centroid(ch_roads.the_geom) from ch_roads,placename_region where st_contains(placename_region.the_geom,ch_roads.the_geom))      But it went wrong. wrong placenames updated in placename field. Please help me to solve this


title: How to stop gdalwarp creating world-spanning outputs near the dateline?

text: I'm using gdalwarp to manipulate SRTM tiles near the the dateline (i.e. 180°, aka the antimeridian). SRTM tiles have a very slight (1/2 pixel) overlap with the meridian. You can see this using gdalinfo:               gdalinfo S16W180.hgt     Driver: SRTMHGT/SRTMHGT File Format     Files: S16W180.hgt     Size is 1201, 1201     [...]     Lower Left  (-180.0004167, -16.0004167) (180d 0' 1.50"W, 16d 0' 1.50"S)     Upper Right (-178.9995833, -14.9995833) (178d59'58.50"W, 14d59'58.50"S)     [...]      So the source spans the dateline by a tiny amount. This causes problems with gdalwarp, which ends up creating huge globe-spanning outputs.               gdalwarp -t_srs "epsg:900913" S16W180.hgt test.tif     gdalinfo test.tif     Driver: GTiff/GeoTIFF     Files: test.tif     Size is 1703, 5     [...]     Lower Left  (-20037508.330,-1806798.473) (180d 0' 0.00"W, 16d 7'13.00"S)     Upper Right (20032839.451,-1689152.120) (179d57'29.01"E, 15d 5'45.84"S)      Note the longitudes span (almost) the entire globe, and also the number of lines is unexpectedly small (5) Is this a bug in gdalwarp? If not, what are the correct options to pass to gdalwarp to get a sensible output?


title: Feature layer not loaded

text: I am developing an Android app using a Webmap with features created on ArcGIS Online. This webmap has a feature layer. I tried to load the basemap and then adding the feature layer, or tried to load the only webmap by its URL. In both ways, I can't display my map. Please, could you help me to find what's the matter? The Layout is:               <TabHost         android:id="@+id/tabHost"         android:layout_width="match_parent"         android:layout_height="match_parent" >              <LinearLayout             android:layout_width="match_parent"             android:layout_height="match_parent"             android:orientation="vertical" >                  <TabWidget                 android:id="@android:id/tabs"                 android:layout_width="match_parent"                 android:layout_height="wrap_content" >             </TabWidget>                  <FrameLayout                 android:id="@android:id/tabcontent"                 android:layout_width="match_parent"                 android:layout_height="match_parent" >                      <LinearLayout                     android:id="@+id/tab2"                     android:layout_width="match_parent"                     android:layout_height="match_parent"                     android:orientation="vertical" >                          <LinearLayout                         android:id="@+id/mapLayout1"                         android:layout_width="match_parent"                         android:layout_height="match_parent"                         android:layout_alignParentLeft="true"                         android:layout_alignParentTop="true" >                              <LinearLayout                             android:id="@+id/buttonLayout"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:layout_alignParentTop="true" >                                  <Button                                 android:id="@+id/centerButton"                                 style="?android:attr/buttonStyleSmall"                                 android:layout_width="wrap_content"                                 android:layout_height="wrap_content"                                 android:text="CENTER"                                 android:textColor="#000000" />                         </LinearLayout>                              <com.esri.android.map.MapView                             android:id="@+id/map"                             android:layout_width="fill_parent"                             android:layout_height="fill_parent"                             android:layout_alignParentLeft="true"                             android:layout_alignParentTop="true">                             <!-- mapoptions.MapType="Topo"  -->                         </com.esri.android.map.MapView>                     </LinearLayout>                          <LinearLayout                         android:layout_width="wrap_content"                         android:layout_height="wrap_content"                         android:layout_alignParentBottom="true"                         android:layout_alignParentLeft="true"                         android:layout_marginBottom="51dp"                         android:layout_marginLeft="126dp"                         android:orientation="vertical" >                     </LinearLayout>                          <ZoomControls                         android:id="@+id/zoomControls1"                         android:layout_width="wrap_content"                         android:layout_height="wrap_content"                         android:layout_alignParentBottom="true"                         android:layout_centerHorizontal="true" />                 </LinearLayout>                      <LinearLayout                     android:id="@+id/tab3"                     android:layout_width="match_parent"                     android:layout_height="match_parent"                     android:background="#158"                     android:orientation="vertical" >                          <LinearLayout                         android:id="@+id/LayoutRow5"                         android:layout_width="match_parent"                         android:layout_height="wrap_content"                         android:layout_marginTop="35dp"                         android:orientation="horizontal" >                              <TextView                             android:id="@+id/bullet1"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:gravity="left"                             android:text="Scarica il regolamento del parco!"                             android:textColor="#00fb4c"                             android:textSize="26dp"                             android:typeface="serif" />                     </LinearLayout>                          <LinearLayout                         android:id="@+id/LayoutRow6"                         android:layout_width="match_parent"                         android:layout_height="wrap_content"                         android:layout_marginTop="20dp"                         android:orientation="horizontal" >                              <TextView                             android:id="@+id/bullet2"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="Informazioni turistiche"                             android:textSize="26dp"                             android:typeface="serif" />                     </LinearLayout>                          <LinearLayout                         android:id="@+id/LayoutRow7"                         android:layout_width="match_parent"                         android:layout_height="wrap_content"                         android:layout_marginTop="50dp"                         android:orientation="horizontal" >                              <TextView                             android:id="@+id/bullet3"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="Eventi della settimana"                             android:textSize="26dp"                             android:typeface="serif" />                     </LinearLayout>                          <LinearLayout                         android:id="@+id/LayoutRow8"                         android:layout_width="match_parent"                         android:layout_height="wrap_content"                         android:layout_marginBottom="68dp"                         android:layout_marginTop="20dp"                         android:orientation="horizontal" >                              <TextView                             android:id="@+id/bullet4"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="Avvisi ai naviganti"                             android:textColor="#f60b00"                             android:textSize="30dp"                             android:typeface="serif" />                     </LinearLayout>                 </LinearLayout>                      <LinearLayout                     android:id="@+id/tab1"                     android:layout_width="match_parent"                     android:layout_height="match_parent"                     android:background="#158"                     android:orientation="vertical" >                          <LinearLayout                         android:id="@+id/LayoutRow1"                         android:layout_width="match_parent"                         android:layout_height="wrap_content"                         android:layout_marginTop="35dp"                         android:orientation="horizontal" >                              <TextView                             android:id="@+id/longTV"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="Longitude : "                             android:textColor="#f66"                             android:textSize="30dp"                             android:typeface="serif" />                              <TextView                             android:id="@+id/longTVvalue"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="999.99° E"                             android:textSize="30dp" />                     </LinearLayout>                          <LinearLayout                         android:id="@+id/LayoutRow2"                         android:layout_width="match_parent"                         android:layout_height="wrap_content"                         android:layout_marginTop="20dp"                         android:orientation="horizontal" >                              <TextView                             android:id="@+id/latTV"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="Latitudine : "                             android:textColor="#f66"                             android:textSize="30dp"                             android:typeface="serif" />                              <TextView                             android:id="@+id/latTVvalue"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="999.99 °N"                             android:textSize="30dp" />                     </LinearLayout>                          <LinearLayout                         android:id="@+id/LayoutRow3"                         android:layout_width="match_parent"                         android:layout_height="wrap_content"                         android:layout_marginTop="50dp"                         android:orientation="horizontal" >                              <TextView                             android:id="@+id/vesselTV"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="Nome barca:"                             android:textColor="#f66"                             android:textSize="30dp"                             android:typeface="serif" />                              <TextView                             android:id="@+id/vesselTVvalue"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="TEST BOAT"                             android:textSize="30dp" />                     </LinearLayout>                          <LinearLayout                         android:id="@+id/LayoutRow4"                         android:layout_width="match_parent"                         android:layout_height="wrap_content"                         android:layout_marginBottom="68dp"                         android:layout_marginTop="20dp"                         android:orientation="horizontal" >                              <TextView                             android:id="@+id/userTV"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="Nome utente:"                             android:textColor="#f66"                             android:textSize="30dp"                             android:typeface="serif" />                              <TextView                             android:id="@+id/userTVvalue"                             android:layout_width="wrap_content"                             android:layout_height="wrap_content"                             android:text="blue-thread"                             android:textSize="30dp" />                     </LinearLayout>                 </LinearLayout>             </FrameLayout>         </LinearLayout>     </TabHost>      The onCreate method of the main activity is written in two different ways. 1st: @Override public void onCreate(Bundle savedInstanceState) {               super.onCreate(savedInstanceState);     setContentView(R.layout.tabs_layout2d);          tabHost=(TabHost)findViewById(R.id.tabHost);     tabHost.setup();     TabSpec spec1=tabHost.newTabSpec("TAB 1");     spec1.setContent(R.id.tab1);     spec1.setIndicator("INFO IMBARCAZIONE");           TabSpec spec2=tabHost.newTabSpec("TAB 2");     spec2.setIndicator("VEDI MAPPA");     spec2.setContent(R.id.tab2);         TabSpec spec3=tabHost.newTabSpec("TAB 3");     spec3.setContent(R.id.tab3);     spec3.setIndicator("SITI");     tabHost.addTab(spec1);     tabHost.addTab(spec2);     tabHost.addTab(spec3);     centerButton = (Button) findViewById(R.id.centerButton);          mercPoint = ToWebMercator(Posizione_dummy_lat, Posizione_dummy_long);          // Set the Esri logo to be visible, and enable map to wrap around date line.     mMapView = new MapView(             this,             "http://www.arcgis.com/home/item.html?id=81d2dcf906cf4df4889ec36c8dc0c1f9",             "", "");         mMapView.inflate(getApplicationContext(), R.id.map, null);          mMapView.setEsriLogoVisible(true);     mMapView.enableWrapAround(true);     mMapView.centerAt(mercPoint, true);      } 2nd: @Override public void onCreate(Bundle savedInstanceState) {               super.onCreate(savedInstanceState);     setContentView(R.layout.tabs_layout2d);          tabHost=(TabHost)findViewById(R.id.tabHost);     tabHost.setup();     TabSpec spec1=tabHost.newTabSpec("TAB 1");     spec1.setContent(R.id.tab1);     spec1.setIndicator("INFO IMBARCAZIONE");           TabSpec spec2=tabHost.newTabSpec("TAB 2");     spec2.setIndicator("VEDI MAPPA");     spec2.setContent(R.id.tab2);         TabSpec spec3=tabHost.newTabSpec("TAB 3");     spec3.setContent(R.id.tab3);     spec3.setIndicator("SITI");     tabHost.addTab(spec1);     tabHost.addTab(spec2);     tabHost.addTab(spec3);     centerButton = (Button) findViewById(R.id.centerButton);          mercPoint = ToWebMercator(Posizione_dummy_lat, Posizione_dummy_long);          mMapView = (MapView) findViewById(R.id.map);     mMapView.setMapOptions(new MapOptions(MapType.TOPO));      ArcGISFeatureLayer webLayerParco = new ArcGISFeatureLayer("http://www.arcgis.com/home/item.html?id=2fd577ea8e534354a13afe236a4d0441", ArcGISFeatureLayer.MODE.ONDEMAND);      mMapView.addLayer(webLayerParco); mMapView.setEsriLogoVisible(true); mMapView.enableWrapAround(true); mMapView.centerAt(mercPoint, true); }


title: Creating DSM and DTM from a LAS using a freeware solution

text: Is it possible to export a surface and terrain model files from an aerial, LiDAR point cloud (LAS) set? I am using SAGA GIS for very basic operations, but I couldn't find this option (and I have a hunch SAGA GIS should support both operations).


title: Edit Attributes in User Form

text: I am trying to convert my old ArcGIS 9.3 (VBA) code to ArcGIS 10 using Visual Basic 2008 Express Edition. So far I have sucessfully created an pop-up window that displays the attributes of a selected feature in text boxes on the form. * * * I added a reference to the editor and made some adjustments to my code. My form still crashes when I try to press the button. Am I missing something or am I just referencing the attribute table in correctly? THANKS!!!               Private Sub cmdUpdate_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles cmdUpdate.Click         pMxDoc = My.ArcMap.Document         Dim pMap As IMap         pMap = pMxDoc.FocusMap         Dim pFLayer As IFeatureSelection         Dim Count As Integer              'Get a reference to the editor.         Dim editorUID As ESRI.ArcGIS.esriSystem.UID         editorUID = New ESRI.ArcGIS.esriSystem.UID         editorUID.Value = "esriEditor.Editor"         Dim editor As IEditor         editor = My.ArcMap.Application.FindExtensionByCLSID(editorUID)              'Get FeatureLayer         pFLayer = pMap.Layer(Count)         Dim pFCursor As IFeatureCursor = Nothing         pFLayer.SelectionSet.Search(Nothing, False, pFCursor)         Dim pF As IFeature         pF = pFCursor.NextFeature              'Get Dataset from FeatureLayer         Dim pEditDataset As IDataset         pEditDataset = pFLayer              'If not editing, start editing dataset's Workspace         If editor.EditState <> esriEditState.esriStateEditing Then             editor.StartEditing(pEditDataset.Workspace)         End If              'Initialize first (only) feature         Dim pFeat As IFeature         pFeat = pFLayer.NextFeature              'Update attributes with user input         pFeat.Value(pFeat.Fields.FindField("YEAR")) = txtPrjYear.Text         pFeat.Value(pFeat.Fields.FindField("PROJECTYPE")) = txtPrjType.Text         pFeat.Value(pFeat.Fields.FindField("LOCATION")) = txtLocation.Text         pFeat.Value(pFeat.Fields.FindField("AGENCY")) = txtAgency.Text         pFeat.Value(pFeat.Fields.FindField("DESCRIPTION")) = txtDescrip.Text         pFeat.Value(pFeat.Fields.FindField("COST")) = txtCost.Text         pFeat.Value(pFeat.Fields.FindField("DATEEDITED")) = txtDate.Text         pFeat.Value(pFeat.Fields.FindField("USEREDITED")) = txtUser.Text              'Stop editing and save edits         editor.StopEditing(True)              'Refresh map         pMxDoc.ActiveView.Refresh()              'Close form         Me.Close()     End Sub


title: Reuse same style settings for chloropleth maps based on different datasets

text: I am new to QGIS. I have been using MapInfo for a couple of years, but I am considering to switch to QGIS. So far I like it, but there are some things that I haven't found out yet. I want to make multiple graduated chloropleth maps based on a varying set of data but with the same graduated/rule based style. For example, I have a polygon vector maps with postal code areas and would like to give the polygons a color based on the drivetime to a certain location A. I can make this style as a graduated style or as a rule based style and save it. Next, I need to make the exact same map for a different dataset, say drivetime to location B. I want to reuse the same settings (colors en cut-off points for the categories). In mapinfo this was very easy: Import an excel file with different datasets in each column-> make 1st map and save the thematic map settings -> Make a new map based on a new column but with the saved thematic map settings. How do I do this in QGIS without -having to remake the graduated style each time (saving a graduated style remembers which column was used and is gone when you change it) or -redefining all rules in the rule-based style? Or is there a way to make a batch of maps that does work in QGIS 1.8? The tutorial on batching maps does not work anymore in 1.8 (http://www.qgis.nl/2012/07/28/series-kaarten-genereren/?lang=en ) Hope someone can help me. THanks in advance, niels


title: Not empty geometry having blank text

text: i'm having problem with geometry column, for some rows i have geometry with blank text, this query will describe more the problem :               select geom, ST_AsGeoJSON(geom), ST_Isvalid(geom), ST_IsEmpty(geom) from polygons      ![enter image description here](http://i.stack.imgur.com/3JK3H.png) so the rows do not have empty geom but i can see their geometries using qgis. any help to solve this problem


title: gdal2Tiles output differs from MapTiler

text: I'm in need to process a large batch of Landsat images with MapTiler and since the version I'm using (1.0 Beta 2) does not support multiple input files, I created a simple python script that calls gdal2tiles for every file. The problem is that using, apparently, the same settings, I get different coordinates. This is the correct output from MapTiler:                <?xml version="1.0" encoding="utf-8"?>         <TileMap version="1.0.0" tilemapservice="http://tms.osgeo.org/1.0.0">           <Title>mar_apr.tif</Title>           <Abstract></Abstract>           <SRS>EPSG:900913</SRS>           <BoundingBox minx="-26.36147186954573" miny="-65.83667775218426" maxx="-21.96598523780907" maxy="-62.31334995428857"/>           <Origin x="-26.36147186954573" y="-65.83667775218426"/>           <TileFormat width="256" height="256" mime-type="image/png" extension="png"/>           <TileSets profile="mercator">             <TileSet href="6" units-per-pixel="2445.98490468750010" order="6"/>             <TileSet href="7" units-per-pixel="1222.99245234375010" order="7"/>             <TileSet href="8" units-per-pixel="611.49622617187504" order="8"/>             <TileSet href="9" units-per-pixel="305.74811308593752" order="9"/>             <TileSet href="10" units-per-pixel="152.87405654296876" order="10"/>             <TileSet href="11" units-per-pixel="76.43702827148438" order="11"/>             <TileSet href="12" units-per-pixel="38.21851413574219" order="12"/>             <TileSet href="13" units-per-pixel="19.10925706787110" order="13"/>           </TileSets>         </TileMap>      And this is the incorrect output from gdal2Tiles:               <?xml version="1.0" encoding="utf-8"?>         <TileMap version="1.0.0" tilemapservice="http://tms.osgeo.org/1.0.0">           <Title>mar_apr_CLIP1.tif</Title>           <Abstract></Abstract>           <SRS>EPSG:900913</SRS>           <BoundingBox minx="-25.32504747468283" miny="1.94797423573108" maxx="-21.33272919903126" maxy="5.10698976386579"/>           <Origin x="-25.32504747468283" y="1.94797423573108"/>           <TileFormat width="256" height="256" mime-type="image/png" extension="png"/>           <TileSets profile="mercator">             <TileSet href="6" units-per-pixel="2445.98490468750010" order="6"/>             <TileSet href="7" units-per-pixel="1222.99245234375010" order="7"/>             <TileSet href="8" units-per-pixel="611.49622617187504" order="8"/>             <TileSet href="9" units-per-pixel="305.74811308593752" order="9"/>             <TileSet href="10" units-per-pixel="152.87405654296876" order="10"/>             <TileSet href="11" units-per-pixel="76.43702827148438" order="11"/>             <TileSet href="12" units-per-pixel="38.21851413574219" order="12"/>             <TileSet href="13" units-per-pixel="19.10925706787110" order="13"/>           </TileSets>         </TileMap>      Does MapTiler perform any other operations besides running gdal2Tiles? Thanks in advance!


title: calculate azimuth from two latitude/longtiude points

text: I have two points with longitude and latitude. If I connected those two points with an arc (on the earth's surface), is there any way to calculate:   1. angle between geographical North and this line    2. angle between magnetic North and this line mathematically from just those points with lat/long?


title: FME - calculating distance to points along a line

text: I am new to FME, and trying to get my head around what should be pretty straightforward issues, but I find myself going around in circles Basically I want to determine the distance (chainage) to a series of points on a water main, which consists of 8 connecting line segments. I have managed to create the chain of points every 10m by using the Labeller and then the GeometryCoercer (fme_point). I also joined up the 8 line segments into 1 by using LineJoiner. The question is now how do I get the chainage to each of the points along the line?   * I've tried MeasureGenerator, but that gives me the chainage to the vertices not to the points   * I've tried ShortestPathFinder but they all come out of the Unused port. I also tried it with the TopologyBuilder first, but this had no effect   * I've tried PointOnLineOverlayer but this splits my line.   * I checked FME - Calculate distance of point along a line but couldn't get the answer either. I think I'm close, but anybody got any suggestions?


title: Raster projection fails without error

text: I have a raster in one coordinate system (SC State Plane NAD83 NSRS2007) and need to project it to another (SC State Plane NAD83 HARN). When I run Project Raster, it fails silently: the dialog box output says "Succeeded," but no output is produced. Also, clicking "Close" on the dialog box brings the tool dialog back up, as it would if a tool failed to complete. Without an error message, I am struggling to understand what the cause of the problem is, and how to circumvent or resolve it.


title: Legend text wrapper (ArcGIS)

text: I'm interested in wrapping long legend text (in ArcGIS 9.3.1) so that the text does not exceed the layout extent. I would like to call the "GetFormattedString" Function within the "UIButtonControl1_Click" procedure. Any pointers on how to get the legend text to wrap at 25 characters would be great. Thank you! Here is the code:               Public Sub UIButtonControl1_Click()          'Refresh legend              Dim pMxDoc2 As IMxDocument         Dim pPageLayout2 As IPageLayout         Dim pGC2 As IGraphicsContainer         Dim pElem2 As IElement         Dim pMSF2 As IMapSurroundFrame         Dim pMS2 As IMapSurround              Set pMxDoc2 = ThisDocument         Set pPageLayout2 = pMxDoc2.PageLayout         Set pGC2 = pPageLayout2         pGC2.Reset              Set pElem2 = pGC2.Next             Do Until pElem2 Is Nothing             If TypeOf pElem2 Is IMapSurroundFrame Then                 Set pMSF2 = pElem2                 Set pMS2 = pMSF2.MapSurround                          Dim strFormattedString As String                     strFormattedString = GetFormattedString("TextContentblahblah", "Arial", 5, 12)                  pMS2.Refresh             pMxDoc2.ActiveView.Refresh                  End If             Set pElem2 = pGC2.Next                  Loop          End Sub          Private Function GetFormattedString(strTextContent As String, strFontName As _     String, sngFontSizeInPoints As Single, sngMaxWidthInPoints) As String              Dim varWordArray As Variant         Dim pMxApp As IMxApplication         Dim pAppDisplay As IAppDisplay         Dim pTransformation As ITransformation         Dim pTextFont As IFontDisp         Dim pTextSymbol As ITextSymbol         Dim dblXSize As Double         Dim dblYSize As Double         Dim strGoodWidth As String         Dim strFinalString As String         Dim strTestString As String         Dim i As Integer              '------ Split the string into an array of words         varWordArray = Split(strTextContent, " ", , vbTextCompare)              Set pMxApp = Application 'QI         '------ Get a pointer to the application display         Set pAppDisplay = pMxApp.Display 'QI         '------ Set up the Transformation         Set pTransformation = pAppDisplay.DisplayTransformation         '------ Set up the Font         Set pTextFont = New StdFont         pTextFont.Name = strFontName         pTextFont.size = sngFontSizeInPoints         '------ Set up the TextSymbol         Set pTextSymbol = New TextSymbol         pTextSymbol.font = pTextFont              '------ Add each word into the test string and test for width         pAppDisplay.StartDrawing pAppDisplay.hdc, 0         For i = LBound(varWordArray) To UBound(varWordArray)             If strTestString = "" Then                 strTestString = varWordArray(i)             Else                 strTestString = strTestString & " " & varWordArray(i)             End If             '------ Get the TextSize             pTextSymbol.GetTextSize pAppDisplay.hdc, pTransformation, strTestString, dblXSize, dblYSize             '------ If the word added is < max width keep adding to the line, else make new  line by inserting vbcrlf             If dblXSize < sngMaxWidthInPoints Then                 strGoodWidth = strTestString             Else                 strTestString = varWordArray(i)                 strFinalString = strFinalString & vbCrLf & strGoodWidth             End If              Next i         strFinalString = strFinalString & vbCrLf & strGoodWidth         pAppDisplay.FinishDrawing         'remove first vbcrlf and pass back         GetFormattedString = Right$(strFinalString, Len(strFinalString) - 2)          End Function


title: Any low-cost Virtual Private Servers able to run GeoServer?

text: I'm looking at virtual private server hosting options for GeoServer. At a recent conference, the most often cited cost of such servers was around $50/month. I've found a few places that advertize VPSs for much less, some as low as $17/month. The specs appear to be adequate for installing and running GeoServer, at least for my initial experimental needs. The reason I'm looking at the lower cost options is that I'm not going to be generating any income here. I'm curious to hear from others that may have tried this ultra-low-cost GeoServer route. How did it work out for you, and if so, which low-cost VPS providers seemed to work out the best, if any did at all?


title: Is Geoslavery complicity a violation of GISP code of Ethics?

text: Would developing applications to support Geoslavery be considered a violation of the GIS Professional (GISP) code of ethics ? > Human tracking devices, however, introduce a new potential for real-time > control that extends far beyond privacy and surveillance, _per se_. As a > result, society must contemplate a new form of slavery characterized by > location control. **Geoslavery** now looms as a real, immediate, and global > threat. \- Jerome E. Dobson and Peter F. Fisher, IEEE Technology and Society > Magazine, Spring 2003 Dobson defines geoslavery as: > a practice in which one entity, the master, coercively or surreptitiously > monitors and exerts control over the physical location of another > individual, the slave. What about a applications that support _Volunteered_ Geoslavery? > It is therefore understandable how inciting fear has made the trade-off of > spatial data privacy for security appear on balance to be a good bargain for > many people. - Nancy Obermeyer Thoughts on “Volunteered (Geo)Slavery”, > Indiana State University.


title: How do you deal with (repair/avoid) broken MXD\LYR paths?

text: In relation to this question ... In the past, whenever I cleaned up a bit and moved layers around (or changed their names) my MXD's paths were broken and it was time-consuming to correct them. I'd like to hear from your experiences and maybe learn a few tips from you: How do you deal with broken MXD paths? (I'm looking for solutions for one computer or a simple network drive, so if your solution deals with larger networks or the web please say so).


title: Using java with mapinfo for integrated mapping?

text: Can we use integrated mapping with java? Are there any examples?


title: Can i use conditional statements in cartodb?

text: I want to code my Cartodb.js map so that for each marker, if a column is not empty, it embeds a video into the infowindow, and if the column is empty, then it does noting.I saw this example online and tried to copy it, but it didn't work. http://www.w3schools.com/js/js_if_else.asp for my code, inside the {{content}} area, I said: if ({{video_embed_link}} IS NOT NULL) {< iframe width="430" height="242" src="{{video_embed_link}}" frameborder="0"> < /iframe>} any suggestions on what to do instead?


title: Avoiding upside-down data-driven labels from layers when rotated in print composer

text: I'm making a printed map using QGIS. I have a bunch of roads in one layer. I've labeled them with one field "name" from the layer properties. I've set the placement to "curved", and the results look very good in the canvas. Most roads are labeled, labels follow the streets, and no labels are upside-down. In my print composer, I've rotated the whole map 45 degrees. (I need to print it with this 45 degree rotation for various reasons.) The problem now is that some of my labels are upside-down. Not all of them, since I only rotated 45 degrees. But some. Not a huge deal, but I'd like to figure out a way to fix it.


title: Importing a DEM in ASCII format into Arcmap - Projection issues

text: I've been trying for the last couple of days to import an ascii grid file into Arcmap. However, I'm really battling to project the data. Can anyone please assist? It would be very much appreciated. I've tried defining the projection but to no avail. The data should be located at roughly the following coordinates... 29. S 26. E. Please see .asc header information below... ncols 990 nrows 1140 xllcenter -97500 yllcenter -3237675 cellsize 25 NODATA_value -9999


title: QGIS r.cost / r.drain and errors

text: I'm trying to find the optimal route of several archaeological sites with Qgis 2.2. I've a DEM and r.walk or r.viewshed executed successfully, but r.cost always fails creating the output, however r.cost executed from python script works perfectly... Once I've the r.cost map, it's time to use r.drain! I selected r.cost file as input, a vector layer with initial points, the output file... clicked "run" and QGIS gives me a python error. I've tried r.walk layer, differents point layers, enter coordinates manually, use qgis 2.0 etc. and always I get the same error. Any idea? Thanks,               Traceback (most recent call last):       File "C:/PROGRA~1/QGISVA~1/apps/qgis/./python/plugins\processing\gui\AlgorithmExecutionDialog.py", line 208, in accept         self.setParamValues()       File "C:/PROGRA~1/QGISVA~1/apps/qgis/./python/plugins\processing\gui\AlgorithmExecutionDialog.py", line 143, in setParamValues         self.paramTable.valueItems[param.name]):       File "C:/PROGRA~1/QGISVA~1/apps/qgis/./python/plugins\processing\gui\AlgorithmExecutionDialog.py", line 190, in setParamValue         return param.setValue(value)       File "C:/PROGRA~1/QGISVA~1/apps/qgis/./python/plugins\processing\parameters\ParameterMultipleInput.py", line 73, in setValue         s += self.getAsString(layer)       File "C:/PROGRA~1/QGISVA~1/apps/qgis/./python/plugins\processing\parameters\ParameterMultipleInput.py", line 147, in getAsString         layers = dataobjects.getVectorLayers(self.datatype)       File "C:/PROGRA~1/QGISVA~1/apps/qgis/./python/plugins\processing\tools\dataobjects.py", line 83, in getVectorLayers         if shapetype == ALL_TYPES or layer.geometryType() in shapetype:     TypeError: argument of type 'float' is not iterable


title: Copy ESRI Map by Value (not Reference) in JavaScript

text: Any way that you know of in JavaScript to copy an ESRI map object by value, not by reference? I know that I can do this:               var foo = new esri.Map("foo");     var bar = foo;     foo.id; //"foo"     bar.id; //"foo"     foo.id = "oof";     bar.id; //"oof"      what I want is to make a copy of foo so that               bar.id = "bar";     foo.id; //"foo"     bar.id; //"bar"      Any ideas? So far the best I've come up with to to construct two different maps, and then copy foo's graphics, features, etc over to bar:               var foo = new esri.Map("foo");     //do some stuff to foo, adding layers etc     var bar = new esri.Map("bar");     bar.addLayers(foo.getLayersVisibleAtScale());     //etc


title: Does Arcpy work for Engine or Server the same as Desktop?

text: I am contemplating a project where there is not going to be access to ArcGIS Desktop and wondering if Engine and/or Server install ArcPy and Python? Is this possible? Does anyone have an example?


title: How to find the nearest point by using PostGIS function?

text: table A:               lat | long | the_geom | code | sign      13.8433095 | 100.6360357 | 0101000020E61.... | ABC | start_point     13.7544738 | 100.5459646 | 0101000020E6..... | ABC | end_point     13.4124215 | 100.6232332 | 0101000020E61.... | DEF | start_point     13.2423438 | 100.2324426 | 0101000020E6..... | DEF | end_point      table B:               lat | long | the_geom | code         13.7546285 | 100.5458729 | 0101000020E.... | ABC         13.7546698 | 100.5458513 | 0101000020E.... | ABC     13.7547107 | 100.5458233 | 0101000020E.... | DEF         ...      I would like to find the nearest point(s) of each points (start and end point of table A) compare with every points with the same 'code' in table B ? What's the best PostGIS function/PostgreSQL query to solve this ? What I have tried are .. ST_Distance_Spheroid and ST_DWithin and ST_Distance But I have no idea how to combine them within a single query.


title: How to import SRTM data for the whole planet into TileMill

text: I need to get SRTM elevation data for the whole planet into TileMill. I already walked through the instructions on Steve Bennett's blog and it worked fine for small areas. The problem is now that for many SRTM files the merged shapefile of contours exceedes the maximum size for shapefiles. On the other hand if I extract one shapefile for each SRTM file there would be too many files to be imported into TileMill. Is there any best practice of importing SRTM data for the whole planet into TileMill?


title: Can I use Google Maps as a layer in Arcmap?

text: Is it possible?


title: ArcEngine in WPF -> publish -> creates .msi install file -> installer fails

text: I am trying to install my custom ArcEngine application on another machine. I just uninstalled and reinstalled ArcGis Engine Runtime 10 and authorized it. "Unable to install or run the application. The application requires that assembly ESRI.ArcGIS.Framework Version 10.0.0.0 be installed in the GAC first" The weird thing is that a previous version of my application still runs. Even the most long-shot type advice is welcome. I'm hitting a wall. * * * I removed the reference to ESRI.ArcGIS.Framework, I guess I wasn't using it. So I thought that would fix the problem, but actually the problem persists. So I assume that one of the assemblies I use is dependant on Framework. On the target machine, I went to the assembly directory and there is no Framework assembly on that machine and there is on my machine. So short term quick fix is to get the assembly on that machine, I am installing devkit on the target to see if that gives me the assembly I need.


title: Use arcpy to connect to database without connection file

text: It seems like arcpy can only work with connection files. This creates one of two problems:   1. Deployment requires manual creation of connection files.   2. We have to write a bunch of code to recreate the connection file on each run. (Creating it only the first time is not feasible since the connection information may change.) As a software developer who does not have access to the production environment, I want to minimize manual steps in the deployment, and I want to keep my code simple. I also want this information to be very easy to change and not require opening up ArcMap (which would be manual), so I want to store it in a simple text config file and establish connections on the fly in my script. Is there a way to create a connection to a database without creating a connection file? An "in memory only" connection (like is possible with pretty much any other database technology, _even Oracle_ , or a number of other libraries, like GDAL)? I'm using ArcGIS 10.2, and I will only be establishing query layers.


title: Unique referencing based on Co-ordinates

text: I have a list of 400 nodes, all split across 'maps' which are gridlines 1000x1000m. Each grid area is referenced by the first 3 digits of the X value of the centre and first 3 digits of the Y value of the bottom left corner of the box e.g. if the corner is 205000,195000 then the reference to that map would be 205195, since the trailing zeroes are irrelevant. I can confidently get this to calculate automatically for each map in Field Calculator. However, I want to give a reference to each node based on which map it is in and a unique identifier within that starting from 001. For example; Node 1 [205150,195150] would be 205195-001 Node 2 [205284,195330] would be 205195-002 Node 3 [206202,195330] would be 206195-001 The benefit of this is that each map will work individually and the user can simply refer to each as Node 001, 002...but the central record can identify each node. I can work out the map number bit fine using field calculator, but the unique number afterwards is proving a bit difficult. I would need to check the existing nodes added to that map, find the current highest number and add 1 to that to produce the next one. Any ideas?


title: QGIS - How can I convert a raster file to a vector file without making same value cells into one polygon?

text: How can I convert a raster file to a vector file without making same value cells into one polygon? ![enter image description here](http://i.stack.imgur.com/nCS0e.jpg) In the attached image, the values "2124", "2269", "2321", "2634", and "2920" are the polygons automatically combined. Is there any way to convert each pixel to one polygon?


title: How to find loops in a network

text: I am a new in this field and need help, how to code using SQL spatial a function that retrieve loops in a network.The input, flag: Geography and barrier: Geography and two spatial tables, first one will list the edges(edgekey, from_junction, to_junction, polyline) and the second table containing juctions(junctionkey, point). The output will be a spatial table containing (edgekey, polykey, from_junction, to_junction, polyline), a single polygon will have the same polykey for each from_junction to_junction value.


title: How to make turn feature class from line feature class?

text: I have a road network of one county in Norway . To this road network, I will make myself a turn feature class of all the road barrier in the road network. I have a feature class with all the road barriers in a point feature class that shows where the road barriers are. Is there anyone who can help me to create a model in model builder or a python script that creates a turn feature class from my road barrier point feature class? There are too many road barriers inn the county to do this manually, so I have to automate the method to do it . I have started to make me a buffer zone of about 1 meter around all my points in my point feature class . Then I clip all the buffer zones against the road network . I do this to get a line feature class with the road barriers that I will make the turn feature class from. But my problem starts here : I do not know how to make a turn feature class from my line feature class that is connected to my road network . Can someone please try to help me with this problem?


title: How to define disturbed/clipped contours elevation for Volume Polygon?

text: I want to calculate surface volume for different layer on TIN surface. I am working on an area in Norway which is 1 km * 1 km. The aim is to calculate surface volume of different imaginary elevation. Because it is a clipped area of a big Island,some contours have been clipped as well ! Now the problem is those areas that have been clipped considering as zero or plane surface, although they are almost high hills and have around 140 m height ! Is there any way to prevent this error or fix the disturbed contours?? I attached the picture for become more clear! ![enter image description here](http://i.stack.imgur.com/zu7ws.jpg)


title: topojson editor/viewer

text: I'm converting a map from geojson to topojson, and I'd like to validate that the conversion is ok. Unfortunately qgis doesn't seem to support topojson yet. I've found this excellent application http://geojson.io/, that worked really well, but now is giving me a strange error. So I'd like to know what editors/viewers are there available with topo


title: .Net to Python - ArcGIS Add-in

text: **Background:** I have been tasked with porting a 9.3 .NET toolbar to ArcGIS 10 / 10.1. I am a dyed in the wool Python programmer and, while digging through someone else's .Net is interesting, it isn't the most fun in the world! Once ported I have a number of improvements that I want to implement. The existing tool is a toolbar with 4 buttons. Each button launches a menu that accepts parameters, performs some processing* and generates output files (shapefile in GDB and text logs). **Question** Is it possible to create this same style workflow, `toolbar > button > menu > run`, using python add-ins? This needs to be shippable to a wide range of users (10 & 10.1) so I do not want to generate windows using WxPython (for example). *By in large the processing is not digging into ArcObjects geoprocessing functionality that is not accessible via Python. Anything geoprocessing that is, is pretty easily recoded in pure Python. Edit: Both posters so far have noted that the functionality I am really looking for is ArcGIS10.1. I should have noted that I realize this, but still wanted to hear some options. This should be targeted at 10.1+ (for when 10.2 rolls...)


title: What is the difference between triangulation and trilateration?

text: Looking around I noticed that many people exchange the terms (tringlulation and trilateration) for the same sense. What's, actually, the correct sense of Triangulation and what are the differences from Trilateration? thx


title: How to access current job id from geoprocessing script running on esri server?

text: I want to programmatically access the current job id of a geoprocessing script running on ESRI 10.2. The script is what needs to know the job id, not the caller of the script. While searching, I have seen plenty of examples of how the submit job request returns with the job id. This isn't what I want. In my script that is running on the server, I want to include the job id in the logs, however I haven't found how to obtain the job in in the script. I initially checked `arcpy.env`, but I didn't see anything. Where else should I look?


title: How to install SAGA 2.1 with QGIS dufour on Windows?

text: I need that someone tell me how i can install saga in qgis 2.0 correctly in windows, because if I install QGIS 2.0 and when I try to use SAGA I get this Missing dependency.This algorithm cannot be run :-( This algorithm requires SAGA to be run.Unfortunately, it seems that SAGA is not installed in your system, or it is not correctly configured to be used from QGIS


title: Layer map plots on top of one another in R/ggplot2

text: I'm making a choropleth map in R/ggplot2. I have the colored part working just fine. This code               p <- ggplot() +        geom_map(data=c_1860,                 aes(map_id = GISJOIN, fill = AHL002),                map = shp_1860) +       expand_limits(x = shp_1860$long, y=shp_1860$lat) +       scale_fill_gradient2(low="white", high="red") +       theme_tufte()      print(p)      produces this map: ![choropleth](http://i.stack.imgur.com/Nuwoj.png) I can also make an outline of the United States without a problem. This code               us     <- map_data("usa")     p <- ggplot() +        geom_path(data=us,                 aes(x=long, y=lat, group=group),                 color="black", size=0.5) +       theme_tufte()      print(p)      produces this map: ![outline](http://i.stack.imgur.com/J4VUR.png) I want to draw the outline around the choropleth to give the map a shape, but that doesn't work. This code just produces the same output as the choropleth map.               p <- ggplot() +        geom_map(data=c_1860,                 aes(map_id = GISJOIN, fill = AHL002),                map = shp_1860) +       expand_limits(x = shp_1860$long, y=shp_1860$lat) +       scale_fill_gradient2(low="white", high="red") +       geom_path(data=us,                 aes(x=long, y=lat, group=group),                 color="black", size=0.5) +       theme_tufte()      print(p)      How can I put the outline of the US around the choropleth map?


title: Shapefile point coordinates to lonlat

text: I'm currently trying to parse a shapefile (http://www.esri.com/library/whitepapers/pdfs/shapefile.pdf) with c++. While I have no problems parsing all the data I need, there is a problem with the coordinates. I got a description file for the data, which states that the used projection is "ETRS89 (Ellipsoid: GRS80), UTM (Zone 33)". I need all points to be in lon/lat format. I spent the whole day googling around to find a solution for this problem, but there doesn't seem to be one. data example:               33354891.477852, 6060028.67074774             ^               ^            X (E)           Y (N)


title: How do I execute a script in the arcpy window?

text: I've been writing python in my text editor so I have a python script, script.py. I'm still developing it so don't want to put it in my toolbox yet, and at the moment this means copying and pasting it from my text editor to the python window in arcgis. Is there a command that lets me call it from the arcgis python window? Something along the lines of:               run script.py      It seems logical that it should exist, but perhaps not.


title: How to turn off layers in ArcMap via arcpy?

text: I'm working on a script in ArcMap 10.1 that takes a screenshot and I need to know how to extinguish (turn off) a layer in my table of contents.


Is there a way to run st_split from within QGIS? I am looking to split a polygon with a line feature. My data is stored in a PostGIS database, and I would like to perform this operation from within QGIS. I am aware of PG_Query for QGIS, but st_split is not a tool within the plugin. Thanks...


title: Can I add GeoExt or ArcGIS JS API to Openlayers + Dojo App?

text: I've written a basic map app with Dojo and Openlayers (using titlebar , toggle buttons , checkbox tree and some more). I've seen ArcGIS JS Api has nice tooltip / balloon and GeoExt has nice Legend Panel which I want to add both to my app. ArcGIS looked a little bit complicated to add (it interrupted Dojo) so I passed it. Can I add only GeoExt just for the legend Panel or will I need to choose either GeoExt or Dojo in my app ?


title: Why do I get huge memory and CPU spikes setting up the geocode tool

text: I am running the geocode tool against a database table which contains values for address, city, state, zip. I am using an address locator that I downloaded locally http://gis.ats.ucla.edu/ Whenever I bring up the window for the tool and then select the address locator and then subsequently select the fields I get a huge CPU and Memory spike in ArcMap. Almost 1GB of RAM and 100% CPU. This last about 10 seconds for each field and the UI is completely unresponsive during this time. The process eventually works but it looks broken. Using Geocode Address from the right click menu doesn't give these problems, only the tool itself. It's not a huge problem but it makes demonstrations look very shaky. Thanks for any help.


title: OpenLayers feature styling change by the baselayer

text: How can I have the style of the features on a vector layer depending on the OpenLayers' "base layer"? For example, some base layers are dark images and I need (for the contrast) the features to be displayed bright in the vector layer upon the base layer. And as soon as base layer changes to a bright one, features should be redrawn dark. How can I achieve this with OpenLayers code?


title: Raster Reclassify fails upon Completion in ArcGIS

text: When attempting to reclassify a raster in a student edition of ArcInfo 10, the program crashes without an error upon reaching 100%. The reclassified raster is written, but with no spatial info. I've got Spatial Analyst and the latest service pack. I suspected it might be a RAM issue, so I split the raster, but I'm having the same problem even with the much smaller dataset. Edit: The plot thickens - I copied the file and mxd to a different location on the desktop, and now the reclassify works! Could having the original file located inside too many folders possibly have caused the problem?


title: Get feature from two overlapping KML layers doesn't work

text: I did an implementation similar to http://openlayers.org/dev/examples/sundials.html The difference is that I have two or more kml layers at the same location and want to show their information combined. The problem is that the _featureselected_ layer's event fires only once, for the top most kml layer. Does anyone know how to get this working, or a possible workarround? Here is my code:               function createFetureInfoService() {     var wmsServerArray = [];     var dictionary = {};     var items = [];     var location;          for (var i = 0; i < wmsControls.length; i++) {         wmsControls[i].deactivate();     }          for (var i = 0; i < kmlControls.length; i++) {         mapa.removeControl(kmlControls[i]);         kmlControls[i].deactivate();     }          kmlControls.length = 0;     wmsControls.length = 0;     totalWMSServerCount = 0;     currentWMSServerCount = 0;     totalKMLCount = 0;     currentKMLCount = 0;          for (var i = 0; i < mapa.getNumLayers(); i++) {         layer = mapa.layers[i];         if ((layer instanceof OpenLayers.Layer.WMS) && !layer.isBaseLayer) {             if (!dictionary.hasOwnProperty(layer.url)) {                 totalWMSServerCount++;                 wmsServerArray.push(layer.url);                 dictionary[layer.url] = "";             }         }         else if ((layer instanceof OpenLayers.Layer.Vector) && (layer.hasOwnProperty("options")) && (layer.options.hasOwnProperty("protocol")) &&         (layer.options.protocol.hasOwnProperty("format")) && (layer.options.protocol.format instanceof OpenLayers.Format.KML) && !layer.isBaseLayer) {             layer.displayOutsideMaxExtent = true;             select = new OpenLayers.Control.SelectFeature(layer);             kmlControls.push(select);             layer.events.on({                 featureselected: function (e) {                     var feature = e.feature;                     location = e.xy;                     items.push({                         xtype: "propertygrid",                         autoHeight: true,                         title: feature.layer.name,                         source: feature.attributes                     });                     currentKMLCount++;                     showGetFeaturePopup(items, location);                     feature.layer.selectedFeatures.length = 0;                 }             });             mapa.addControl(select);             select.activate();         }     }          for (var i = 0; i < wmsServerArray.length; i++) {         var wmsFeatureInfo = new OpenLayers.Control.WMSGetFeatureInfo({             url: wmsServerArray[i],             map: mapa,             infoFormat: "application/vnd.ogc.gml"         });         wmsControls.push(wmsFeatureInfo);         wmsFeatureInfo.events.on({             getfeatureinfo: function (e) {                 location = e.xy;                 Ext.each(e.features, function (feature) {                     items.push({                         xtype: "propertygrid",                         autoHeight: true,                         title: feature.fid,                         source: feature.attributes                     });                 });                 currentWMSServerCount++;                 showGetFeaturePopup(items, location);             },             beforegetfeatureinfo: function (e) {                 this.layers = filterQueryableLayers(this.url);                 },             nogetfeatureinfo: function (e) {                 currentWMSServerCount++;                 showGetFeaturePopup(items, location);             }         });         wmsFeatureInfo.activate();     }       }               function showGetFeaturePopup(items, location) {     totalKMLCount = 0;     var layers = mapa.getLayersBy("visibility", true);         for (var i = 0, len = layers.length; i < len; i++) {         if ((layers[i] instanceof OpenLayers.Layer.Vector) && (layers[i].hasOwnProperty("options")) && (layers[i].options.hasOwnProperty("protocol")) &&         (layers[i].options.protocol.hasOwnProperty("format")) && (layers[i].options.protocol.format instanceof OpenLayers.Format.KML) && !layers[i].isBaseLayer) {             totalKMLCount++;         }     }          if ((currentWMSServerCount + currentKMLCount) < (totalWMSServerCount + totalKMLCount))         return;     currentWMSServerCount = 0;     currentKMLCount = 0;          if (items.length == 0)         showInfoPopup(mapa.getLonLatFromPixel(location), 'Nenhuma informa\u00E7\u00E3o dispon\u00EDvel');     else {         var popUp = new GeoExt.Popup({             title: "Informa\u00E7\u00E3o sobre as Camadas WMS",             layout: "accordion",             anchored: false,             resizable: true,             map: mapa,             location: location,             items: items,             autoScroll: true,             width: 600,             height: 600,             unpinnable: false,             draggable: true,             closeAction: 'hide'         });         popUp.show();     }     items.length = 0;      }


title: PyQGIS - Start/Stop Editing with data Provider - Commiting changes

text: I want to edit a specific field for a `QgsVectorLayer`, within a QGIS application. I am facing inconsistent results from my code, ie: changes are not alway commited. Please could someone point me to where it is wrong within this loop:               #vlayer is a QgsVectorLayer, imported from SQLITE DB     vlayer = QgsMapLayerRegistry.instance().mapLayer(layer)     caps = vlayer.dataProvider().capabilities()                     polygonGeom = QgsGeometry( rubberBand.asGeometry() )      #if vlayer is capable of edition     if caps & QgsVectorDataProvider.ChangeAttributeValues:              vlayer.selectAll()              for feature in vlayer.selectedFeatures():                  if feature.geometry() == None:                 continue                  if feature.geometry().intersects(polygonGeom):                 fid = feature.id()                 idx = vlayer.fieldNameIndex('some field')                 attr ={idx:'some value'}                 print fid,idx,attr                 vlayer.dataProvider().changeAttributeValues({fid:attr})              vlayer.invertSelection()      Here is a post of my console showing that some features are selected and should be edited, but again, some are not properly updated to the underlying database.               2 16 {16: 4}      12 16 {16: 4}      13 16 {16: 4}


Postgres: How do I get the closest distance for each element within a point coordinate table from a polyline table? And if this is possible what version of Postgres is this doable.


title: Def Query to not display feature if it isn't attributed

text: I have a bunch of features that I don't want to keep checking if something is attributed. I I want to query it out to not show the geometry if a column has the default string 'NOT ATTRIBUTED' to it. I can get it to run "Name = 'NOT ATTRIBUTED'" to work, but Name = NOT 'NOT ATTRIBUTED' doesn't work. It won't run at all using the 'IS' It has to use the equals sign. It's on a SQL server.


title: QGIS 2.2: Attributes table keeps crashing

text: I am currently working on an intensive project that requires me to edit the fields of literally thousands of line segments. While updating the shapefile, it crashed on me. Today, qgis has probably crashed on me upwards of five times now. Originally, I thought it was just my computer but now it's happened to one of my group members. I tried opening the project file and re-adding a slightly older version of the file, but qgis crashed once more when I tried to look at the attribute table to see any saved changes. Both of us are running qgis 2.2 on mac laptops. The crashing happened in the version of qgis before this one, which is why I upgraded. I thought that maybe there was a bug that had been addressed in the latest version. Has anyone else ran into this issue of editing large files and being unable to open their shapefiles' attribute tables without crashing? Does anyone know any way around this issue? Thanks


title: Boundless Opengeo Suite: Geoserver fails to launch due to Java heap size

text: When I go to start Boundless OpenGeo Suite 3.1 from the OpenGeo Suite Dashboard, Geoserver will launch if I point it to the default geoserver data_dir, but when I point it to mine, which I have used successfully for > 8 months, javaw.exe fills all remaining memory on the machine and Geoserver fails to launch. I am having some trouble interpreting the logs to find the cause of this issue; does anybody know how I might go about finding and fixing the problem? I am new to this community, so if there is any way that I can improve this submission, please let me know and I'd be happy to fix it. Steps taken:   * Restarted Machine   * Uninstalled/Reinstalled OpenGeo Suite   * Pointed Geoserver to the default data_dir Machine:   * Windows Server 2003 Service Pack 2   * AMD Opteron 280 2.41GHz   * 3.83 GB RAM Startup Log (first block): * * * ## \- GEOSERVER_DATA_DIR: e:\path_to_install\data_dir 2014-06-11 16:13:58,423 WARN [referencing.factory] - Axis elements found in a wkt definition, the force longitude first axis order hint might not be respected: PROJCS["WGS84 / Simple Mercator", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS_1984", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree", 0.017453292519943295]], PROJECTION["Mercator_1SP_Google"], PARAMETER["latitude_of_origin", 0.0], PARAMETER["central_meridian", 0.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["m", 1.0], AXIS["x", EAST], AXIS["y", NORTH], AUTHORITY["EPSG","54004"]] 2014-06-11 16:14:02,033 WARN [org.geoserver] - Ignoring store directory 'layergroups' 2014-06-11 16:14:02,595 WARN [org.geoserver] - Ignoring store directory 'styles' 2014-06-11 16:14:04,282 WARN [org.geoserver] - Error connecting to 'Taxlots'. Disabling. 2014-06-11 16:14:04,876 WARN [org.geoserver] - Ignoring store directory 'styles' 2014-06-11 16:14:07,439 WARN [org.geoserver] - Ignoring store directory 'layergroups' 2014-06-11 16:14:07,439 WARN [org.geoserver] - Ignoring store directory 'styles' 2014-06-11 16:14:07,501 ERROR [geotools.jdbc] - There's code using JDBC based datastore and not disposing them. This may lead to temporary loss of database connections. Please make sure all data access code calls DataStore.dispose() before freeing all references to it 2014-06-11 16:14:08,829 WARN [org.geoserver] - Ignoring store directory 'styles' 2014-06-11 16:14:09,454 WARN [org.geoserver] - Ignoring store directory 'styles' 2014-06-11 16:15:04,873 WARN [support.DisposableBeanAdapter] - Invocation of destroy method failed on bean with name 'geoServerLoader': java.lang.NullPointerException 2014-06-11 16:15:04,873 ERROR [context.ContextLoader] - Context initialization failed org.springframework.beans.factory.BeanCreationException: Error creating bean with name 'gwcFacade' defined in URL [jar:file:/E:/path_to_install/OpenGeo%20Suite/webapps/geoserver/WEB- INF/lib/gwc-2.2-SNAPSHOT.jar!/applicationContext.xml]: Cannot resolve reference to bean 'gwcStorageBroker' while setting constructor argument; nested exception is org.springframework.beans.factory.BeanCreationException: Error creating bean with name 'gwcStorageBroker' defined in URL [jar:file:/E:/path_to_install/OpenGeo%20Suite/webapps/geoserver/WEB- INF/lib/gwc-2.2-SNAPSHOT.jar!/geowebcache-core-context.xml]: Cannot resolve reference to bean 'gwcMetaStore' while setting constructor argument; nested exception is org.springframework.beans.factory.BeanCreationException: Error creating bean with name 'gwcMetaStore' defined in URL [jar:file:/E:/path_to_install/OpenGeo%20Suite/webapps/geoserver/WEB- INF/lib/gwc-2.2-SNAPSHOT.jar!/geowebcache-core-context.xml]: Instantiation of bean failed; nested exception is org.springframework.beans.BeanInstantiationException: Could not instantiate bean class [org.geowebcache.storage.metastore.jdbc.JDBCMetaBackend]: Constructor threw exception; nested exception is java.lang.OutOfMemoryError: Java heap space


title: Where does Open-Streetmap forest data come from?

text: If you open this url http://www.openstreetmap.org/?mlat=50.6524658203125&mlon=9.4482421875&zoom=13 you are in a forest in germany. I downloaded the complete Open-Streetmap Data for Germany. I can find the waynodes for the tracks in this area... But not for the forest itself Cna anyone help? Which way (id) defines the forest under the marker? Or is it just the Background image for areas where's nothing defined?


title: render geotiff tiles as a map layer

text: i have the following problem: i have downloaded tiff files with gMapMaker and transformed them to GeoTiff using the associated .map files in QGIS. now i want to render these GeoTiff files to use them as a map (like using Google Maps API) and then apply some geographic info (as an SVG element on a particular geographic point) for example using d3.js or OpenLayers.. The important thing is that i have to do all of this offline (yes, i'm trying to re-create a sort of Google Maps offline).. until now, i tried to use Javascript but i didn't find any library to load GeoTiffs; moreover i can use also other languages (i have no restrictions from this point of view)... can you help me?


title: How do I pin all labels in QGIS?

text: How can I pin all my labels in QGIS 2.2/2.4 (Data Defined Labels) without clicking every single label?


title: Base Layers switching of OpenLayers

text: I want to know what happen when we switch between the baseLayers of an openLayers Map, I mean when we switch the baseLayer, does it use the layer.setVisibility() to false for all the remaining layers or it uses map.removeLayer() for all others


title: Why is osm2pgrouting not importing?

text: I installed pgrouting successfully and did a little import with a 20Mb file. Works good! Now I'm trying to import a file with the size of a country and it doesn't even start to import. The last 2 lines in the console are: "Trying to load data - Trying to parse data" and then it's done. Why is it not working like it should? It's not even a out-of-memory bug, it's just not working! Thanks


title: How to upload a series of OS map sheets using QGIS and merge together?

text: I need to upload historic OS sheets in tiff format, and 'stitch' them together, as some districts are spread over a number of sheets. I have no idea how to do this. Can anyone help?


title: converted SHP to geoJSON - how to put this on OSM map?

text: I'm a beginner at this. I converted a shapefile to a .js file but I am having trouble adding this data to a map using leaflet. The osmscarb.js: var osmscarb = { "type": "FeatureCollection", "features": [ { "type": "Feature", "properties": { "Id": 0, "POINT_X": -8826927.77461, "POINT_Y": 5417539.05834, "INTXN": "Luttrell/Danforth", "FILENAME": "Luttrell Loop.jpg" }, "geometry": { "type": "Point", "coordinates": [ -8826927.774612015113235, 5417539.058339669369161 ] } } I'm putting this in my map:               var myLayer = L.geoJson().addTo(map);     myLayer.addData(osmscarb);


title: search the entire globe for locations with very steep hills

text: What are some databases that provide high-resolution data about altitude of the ground? It seems that with such a database (with resolution to about 10 meters or less), I could run some kind of simple algorithm to calculate "steepness" at any given point then group these points together to find areas, say, bigger than 1 sq km with high "steepness" and then just give me the top of the list.


title: QgsMapLayerRegistry().removeMapLayers(wfs_layer) causes crash

text: I have a WFS layer, served by Geoserver, loaded in QGIS using the standard 'Add WFS layer' dialog. When I try to remove it using:               QgsMapLayerRegistry().instance().removeMapLayer(wfs_layer.id())      I get a full application crash, with no Python errors. This also happens when I use `wfs_layer.reload()` and `wfs_layer.dataProvider().reloadData()`. These methods previously worked on the same WFS layer. The service and content of the layer has not changed. The methods work when manipulating a (PostgreSQL) database layer. I have made no changes to the code that I can see would affect the ability to reload data. Any ideas or experiences?


title: Zoom to location on map creation android

text: So currently I have an application that zoomed to a certain extent. The client would like it to zoom to location. The problem I'm running into is the mapview is declared in the xml and I'm having a hard time figuring out how to make the map zoom to location after it's been created and before a user clicks a button that will zoom to location. I've centerAndZoom(), zoomToResolution(), centerAt() with setScale(). Nothing seems to work. The function I created is getting called and everything seems to be running as it should, but the map is not centering and zooming. Any ideas would be greatly appreciated. I also tried to set the mapOtions programmatically with no success.


title: Needed shape file for Australia and New zealand

text: I was trying to find a shape file which has Australia and New zealand together(states) so that i could convert to jvectormap for ANZ. where i can find it?


title: PostGIS layer attribute table empty in Quantum 2.0.1

text: I recently installed QGIS 2.0.1 (Dufour) on a Windows 7 PC (64-bit). Today I tried, for the first time since Dufour was installed, to load PostGIS layers into QGIS. I successfully set up a connection and I could load layers into QGIS (using both DB Manager and the AddPostGIS layers tool). However, although in DBManager the tables show that they contain a number of features, I cannot see them in the Preview tab (of DBManager). Moreover, the attribute tables of the PostGIS layers I loaded into QGIS are empty (0 features, and the "show selected only" box is not ticked). To make sure it wasn't a problem with my database (PostgreSQL 9.1, PostGIS 2.1.0SVN), I uninstalled QGIS 2.0.1 and re-installed 1.8 (Lisboa). In Lisboa, I can see my layers rendered in QGIS and their attribute tables are not empty. Is this a known problem with version QGIS 2.0.1 or have default settings changed since version 1.8? Thanks for your help.


title: Can I import symbology from an .axl file into ArcMap?

text: I'm setting up new Arcserver services, and some users would like the symbology in the new services to match the old IMS services. With that in mind, is it possible to either: a) import the symbology from an axl file into ArcMap b) convert the axl file to an avl file so the symbology can be imported Thanks in advance for any help, any ideas appreciated.


title: River drainage basin boundary

text: I'm looking for a set of outlines of the river basins in the UK - I don't need a particularly detailed or highly accurate map - just a map saying point X will drain into the Thames or point Y will just go into the sea! Just the main river names would be fine. I have the full Environment Agency rivers set and the OS DTMs if I want... but how do I convert this into a catchment boundary? ps. Sorry about the poor tagging I can't find suitable existing and I don't have the rep to create new ones. Thanks!


title: Buffer in ArcGIS v. ArcINFO 9.3

text: I have ArcGIS 9.3 and my question is regarding Buffering. Is that only allowed in ArcINFO or is there a way to use a buffer in ArcGIS?


title: LatLon or LonLat, what's the "right" way to display coordinates and inputs?

text: I'm trying to get a sense if this is an issue for others or every input/output should be labeled so user is not confused and just go with it? I think almost everyone pronounces it as "LatLon". Who started it? Is it because it's in alphabetical order compared to "LonLat"? Mapping Lat and Lon to Cartesian plane Lon is "x" and Lat is "y" so since we say "(x,y)" it should be said as "LonLat". And now for display of information. Should the status bar on a mapping application display La,Lo or Lo,Lat? Should it just be labeled as one way and let user deal with it? And same with input, what's the right way to order the fields? KML's format is Lon,Lat,Altitude. While other apps is Lat,Lon and so have to be very vigilant about when converting formats. Is there a standard? [this might be a wiki question since I don't really expect to have _a_ right answer, but maybe get a consensus of best practices or a discussion going]


title: How configure ArcGIS 10 Server for .NET at http://server/rootContext/ArcGIS rather than http://server/ArcGIS

text: We are developing a Flex application which will communicate with ArcGIS server via a reverse proxy. The reverse proxy configuration is such that the host and path portions are separate configurations. Does anyone know if you can configure ArcGIS to respond to a URL like: http://serverName/rootContext/ArcGIS rather than the default: http://serverName/ArcGIS Has anyone done this? Any suggestions on how that might be done? Thanks, Carl


title: GIS Solutions And ROI For Local Government

text: In your experience, what GIS solutions provide the greatest ROI for local government's responsibilities? Examples: Parcel management, managing CRM cases from the public, emergency response, etc.


title: Getting all populated places per state when number of places over 6000

text: I am able to query for the number of places. The problem becomes when there are over 6000 places, the system gives me an error message that the requested start row is too big. Here is an example of a request which would generate that: http://api.geonames.org/search?featureCode=PPL&username=genadinik&country=US&style=full&adminCode1=VA&maxRows=1000&startRow=6000 How do people usually get around this? Thanks!


title: How to get parks in American Samoa from GeoNames

text: I have this query which does not return anything: http://api.geonames.org/search?featureCode=PRK&username=genadinik&country=US&style=full&adminCode1=AS&maxRows=1000 Should I change it? Is American Samoa considered its own country? Un-american Samoa? lol. Thanks!


title: How to modify the colors of a raster map to use only one color in varying shades?

text: I have a raster map in GRASS containing elevation values. I want to display the elevation values in only one color (specificaly blue) but in varying shades of blue. I tried the options in r.colors, but none of those displays the map in monochrome.


title: VBA to VB.NET (Visual Studio 10) Migration - ArcMap Focus

text: I will be migrating several VBA applications currently in ArcMap 9.3.1 to Visual Studio 10 in the near future, and wanted some feedback on the general migration steps listed below. I wanted to make sure that I did not miss any critical steps and that the order of the steps is correct. 1.) Run the ArcGIS Code Converter tool to update old Object names. 2.) Create a Visual Basic project and add all the mxd specific modules, and classes to the new VB project. Re-create any forms or command buttons that existed in VBA. 3.) Run Code Advisor for Visual Basic 6.0. 4.) Open the VB project in Visual Studio 10 and run the Visual Basic Upgrade Wizard. 5.) Manually update any common dialogs, MsgBox, and Win32 API code sections. 6.) Add ESRI references and libraries. 7.) Register new Visual Studio dll. Here are some website references that I have used so far. Migrating VBA/VB6 ArcObjects Applications to .NET http://resources.arcgis.com/gallery/video/arcgis- desktop/details?entryID=A721AC1F-1422-2418-341B-DFD761945F47 How to migrate from VB6 to VB .NET (wiring events) http://resources.esri.com/help/9.3/arcgisengine/dotnet/ca0956fb-b3d8-48a0-8c85-3559ae02dc0f.htm ArcGIS Application Migration: Using the Visual Studio 2008 Upgrade Wizard http://resources.arcgis.com/gallery/video/arcgis- desktop/details?entryID=93EC2C14-1422-2418-880E-0CFF2043C7FD Thank You


title: Circle features missing in PostGIS after import from ArcGIS using ogr2ogr

text: I have a ArcGIS feature class of geometry type polygon. It was created as a convex hull with a 200 meter buffer on a point dataset, aggregating rows on a particular field. The number of points used to generate each polygon is variable, starting at 1. I use ogr2ogr to upload from a personal geodatabase (.MDB) to a remote PostGIS database, and MapServer 5.4 to display. All features are uploaded successfully, but those polygon features that were generated from a single point do not plot using MapServer. Based on how the polygons are generated, all those based on a single point ought to be circles with radii of 200 meters. In the PostGIS table the geometry field is empty for all those features, and following a tip I found here I tried "select st_npoints(wkb_geometry) from merged_dist_areas;" and it appears all the "missing" polygon features have no points at all, which may explain why they do not plot. What can I do to plot these circles? Thanks!


title: RuntimeError from r.viewshed in graphic modeler

text: after using r.los in QGIS want to switch to r.viewshed in GRASS GIS 7. i need to do a lot of viewshed analysis form given points of a vectorlayer. if i try to add r.viewshed to the graphic modeler i get the following output:                   Traceback (most recent call last):           File "C:\OSGeo4W64\apps\grass\GRASS-7\etc\gui\wxpython\gmodeler\frame.py", line 698, in OnAddAction         completed = (self.GetOptData, action, action.GetParams()))           File "C:\OSGeo4W64\apps\grass\GRASS-7\etc\gui\wxpython\gui_core\forms.py", line 2363, in ParseCommand         get_dcmd = get_dcmd, layer = layer)           File "C:\OSGeo4W64\apps\grass\GRASS-7\etc\gui\wxpython\gui_core\forms.py", line 458, in __init__         frame = self)           File "C:\OSGeo4W64\apps\grass\GRASS-7\etc\gui\wxpython\gui_core\forms.py", line 1562, in __init__         param = p)           File "C:\OSGeo4W64\apps\grass\GRASS-7\etc\gui\wxpython\gui_core\gselect.py", line 2148, in __init__         switcher = self._giface.GetMapDisplay().GetToolSwitcher()           File "C:\OSGeo4W64\apps\grass\GRASS-7\etc\gui\wxpython\gmodeler\giface.py", line 23, in __getattr__         return getattr(self._giface, name)           File "C:\OSGeo4W64\apps\grass\GRASS-7\etc\gui\wxpython\gmodeler\giface.py", line 23, in __getattr__         return getattr(self._giface, name)           File "C:\OSGeo4W64\apps\grass\GRASS-7\etc\gui\wxpython\gmodeler\giface.py", line 23, in __getattr__         return getattr(self._giface, name)           File "C:\OSGeo4W64\apps\grass\GRASS-7\etc\gui\wxpython\gmodeler\giface.py", line 23, in __getattr__     ***...and so on, till***         return getattr(self._giface, name)         RuntimeError         :         maximum recursion depth exceeded while calling a Python object      is there a way to solve this propblem? i'm on windows 7. thx - niels... edit: System Info   GRASS Version: 7.0.svn   GRASS SVN Revision: 59235M   Erstellungsdatum: 2014-03-12   GDAL/OGR: 1.10.1   PROJ.4: 4.8.0   GEOS: 3.4.2   SQLite: 3.7.17   Python: 2.7.4   wxPython: 2.8.12.1   Platform: Windows-7-6.1.7601-SP1


title: Is it possible to create road-type dependent 'shields' in QGIS?

text: In QGIS, I need to place road numbers on a map (they are point labels) in solid rectangular boxes to make them stand out. Ideally they need to be the same colour as the road they refer to, so a motorway number should be in a blue box because the motorways are blue etc. The boxes also need to scale with the label size, as some roads are for example 'A9', whereas others are 'A4567'. The roads and their labels have codes defining their type so this could be used to set the colours. Is this possible?


title: enhanced grid dojo

text: I have a widget which consists of two `panels`.I need to place a dojo `enhancedGrid` in both the panels . However,the `enhancedGrid` appears only in one of the two panels despite the fact that I have created two `enhancedGrids` , one for each panel. I am not able to figure out why this is happening and i cannot post my code because my nature of work doesnot permit it. Any help will be appreciated.


title: Moving points onto lines (~neighborhood)

text: I do have two vector layer, of which one is a point layer based on "events" by remote sensing and the second one is a line layer from local research. In my case these are earthquakes and tectonic faults, but I guess one could simply choose "car-accidents and roads" as a general example. So what I'd like to do is move/copy the points onto the closest point of the lines, as long as its within a tolerance distance (say 1-2km or 0.0xx°), with the new point layer (+attr moved y/n). Any ideas ? Linux, QGIS 1.8


title: convert a dbf into a raster image using python? work around for arcpy

text: ArcGIS' arcpy does not have a function to add a table to a layout. I was thinking that maybe taking that table (e.g. a dbf), there may be some python script command to turn it into a raster (e.g. do something like a screen capture on the table). You could then add the screen capture raster image to the layout using arcpy. Anyone know of a function in python, or a few lines of python code to do this? Thanks for your advice.


title: How to calculate the percentage of a Line through Polygons?

text: I have a straight line which runs through a bunch of polygons, which don't have spacing between them (border to border). Is there an easy way with QGIS to calculate which polygons it passes and by what percentage of the total distance its within these polygons? I'm using Ubuntu, QGIS 1.8.


title: How can I solve "The Network Analyst solver could not determine a solution"?

text: ![enter image description here](http://i.stack.imgur.com/Ezdi2.png) I am using ArcGIS 10.0. I have network data for undulated terrain so there are lots of inaccessible parts of the network. When I was using this model, I get "The Network Analyst solver could not determine a solution" error. I'm not aware how to solve this error. I have attached an example image to show the problem. I am looking for your guidance. Thanks in advance.


title: Have editable label locations (dynamic) on a layer on which there is query

text: Is it possible to have dynamic labels, set by editable x y parameters in the attribute table, on a layer which also has a query applied, in QGIS? i.e. I want to set the query and get all the peripheral data off the map canvas, then align the remaining labels with the "move label" tools in QGIS 1.8. If I create a second shapefile just for labelling, later on when the master dataset has edits, I end up loosing the x y coords established for the labels.


How to rank polylines according to their similarity to a ref polyline? ![enter image description here](http://i.stack.imgur.com/0cOP4.png) **Update 1:** simulations show that simple $\sqrt{\sum{(case-ref)^2}}$ can find better similarity between polylines (as set of points) than Hausdorff distance! but not satisfied yet! (right-click >> view image >> for full resolution) ![enter image description here](http://i.stack.imgur.com/mUV0q.png)


title: CORS lets one ArcGIS JavaScript (REST) request through but not the other

text: I have a problem with Cross-Origin-Requests invoked by the ArcGIS Javascript API. I'm able to access my web app without any issues if I use `localhost:8080`. But a REST call (by the ArcGIS JS API) fails if I use `mymachine.local`. The `exportImage` REST call works properly but the `query` REST call invoked by `featureLayer.selectFeatures` fails with a _Not allowed by Access-Control- Allow-Origin_ error. I've gone through http://blogs.esri.com/esri/arcgis/2011/09/29/cross-origin- resource-sharing-cors-with-the-arcgis-api-for-javascript/ but I don't think there's a problem with the CORS setup as one request goes through and the other one fails.


title: Why does Raster Calculator return GRID ERROR-Unable to open the input grid: 1+FLOATC:\document...?

text: I'm calculating Simple Ratio (SR) from NDVI using Raster calculator in Arcgis. Here is expression I wrote in Raster calculator. SR=(1+float[NDVI_apr.tif])/(1-float[NDVI_apr.tif]) But when I evaluate this expression, the following error always appear on the screen. **GRID ERROR- Unable to open the input grid: 1+FLOATC:\document... Error in executing grid expression.** What is wrong with this expression? Thank you all