Feature layer not loaded

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); }


Creating DSM and DTM from a LAS using a freeware solution

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).


Edit Attributes in User Form

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


Reuse same style settings for chloropleth maps based on different datasets

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


Not empty geometry having blank 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


gdal2Tiles output differs from MapTiler

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!


calculate azimuth from two latitude/longtiude points

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?


FME - calculating distance to points along a line

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?


Raster projection fails without error

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.


Legend text wrapper (ArcGIS)

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


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

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?


Is Geoslavery complicity a violation of GISP code of Ethics?

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.


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

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).


Using java with mapinfo for integrated mapping?

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


Can i use conditional statements in cartodb?

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?


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

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.


Importing a DEM in ASCII format into Arcmap - Projection issues

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


QGIS r.cost / r.drain and errors

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


Copy ESRI Map by Value (not Reference) in JavaScript

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


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

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?


How to find the nearest point by using PostGIS function?

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.


How to import SRTM data for the whole planet into TileMill

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?


Can I use Google Maps as a layer in Arcmap?

Is it possible?


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

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.


Use arcpy to connect to database without connection file

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.


Unique referencing based on Co-ordinates

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?


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

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?


How to find loops in a network

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.


How to make turn feature class from line feature class?

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?


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

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)


topojson editor/viewer

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


.Net to Python - ArcGIS Add-in

**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...)


What is the difference between triangulation and trilateration?

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


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

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?


How to install SAGA 2.1 with QGIS dufour on Windows?

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


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

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?


Shapefile point coordinates to lonlat

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)


How do I execute a script in the arcpy window?

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.


How to turn off layers in ArcMap via arcpy?

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...


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

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 ?


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

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.


OpenLayers feature styling change by the baselayer

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?


Raster Reclassify fails upon Completion in ArcGIS

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?


Get feature from two overlapping KML layers doesn't work

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;      }


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

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.


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

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.


QGIS 2.2: Attributes table keeps crashing

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


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

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


Where does Open-Streetmap forest data come from?

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?


render geotiff tiles as a map layer

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?


How do I pin all labels in QGIS?

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


Base Layers switching of OpenLayers

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


Why is osm2pgrouting not importing?

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


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

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?


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

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);


search the entire globe for locations with very steep hills

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.


QgsMapLayerRegistry().removeMapLayers(wfs_layer) causes crash

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?


Zoom to location on map creation android

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.


Needed shape file for Australia and New zealand

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?


PostGIS layer attribute table empty in Quantum 2.0.1

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.


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

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.


River drainage basin boundary

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!


Buffer in ArcGIS v. ArcINFO 9.3

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?


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

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]


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

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


GIS Solutions And ROI For Local Government

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.


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

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!


How to get parks in American Samoa from GeoNames

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!


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

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.


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

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


Circle features missing in PostGIS after import from ArcGIS using ogr2ogr

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!


RuntimeError from r.viewshed in graphic modeler

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


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

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?


enhanced grid dojo

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.


Moving points onto lines (~neighborhood)

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


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

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.


How to calculate the percentage of a Line through Polygons?

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.


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

![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.


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

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)


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

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.


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

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


Using Union on multiple shapefiles to output single result

Using QGIS 2.2 and Python, I would like to use the Union algorithm on multiple shapefiles (contained in a folder _TEST_ ) which would then produce single shapefile containing all the data. Currently the following code only outputs the union of the first and last file, how can I modify it so that the output of 2 layers is saved and then use that output with another layer etc? UPDATE: The shapefiles are being processed in terms of the `clip_function` however I can't load the temp files when running the `qgis:union` algorithms. Also, an error occurs when it tries to output the Result.               ##Model=name     ##Select_folder=folder     ##Cellsize=number 1000     ##Result=output vector          import os     import glob     import tempfile     from qgis import utils     from qgis.core import QgsVectorLayer, QgsMapLayerRegistry, QgsVectorFileWriter, QgsCoordinateReferenceSystem          #   Create grid with defined parameters     outputs_1=processing.runalg("qgis:creategrid", Cellsize, Cellsize, 24108, 18351.157175, 258293.802316, 665638.226408, 1, 'EPSG:7405', None)          def clip_function():          #   Set directory and search for all .shp files         os.chdir("C:\Users\TEST\DATA\\")         for fname in glob.glob("*.shp"):           #   Clip .shp files within directory to grid and save files as temp             outputs_2=processing.runalg("qgis:clip", outputs_1['SAVENAME'], fname, None)             outputs_3=processing.runalg("qgis:fieldcalculator", outputs_2['OUTPUT'], 'Rank_' + fname, 1, 10, 0, True, 1 , None)                      def binary_alg_wrapper(list_of_files,final_file_name):         run_count=0         temp_path=tempfile.gettempdir()         temp_input_name="temp_input_layer"         temp_output_name="temp_output_layer"         temp_input_layer_filename=str(os.path.join(temp_path,temp_input_name+".shp"))         temp_output_layer_filename=str(os.path.join(temp_path,temp_output_name+".shp"))         for file_name in list_of_files:             EPSG_Code_for_file=4326             run_count+=1             if run_count==1:                 #Copy the first file into a temp location                 temp_output_layer = QgsVectorLayer(file_name, "temp_output_layer", "ogr")                 Add_result = QgsMapLayerRegistry.instance().addMapLayer(temp_output_layer)                 write_error = QgsVectorFileWriter.writeAsVectorFormat(temp_output_layer,temp_output_layer_filename,"system",QgsCoordinateReferenceSystem(EPSG_Code_for_file), "ESRI Shapefile")             QgsMapLayerRegistry.instance().removeMapLayer(temp_output_layer.id())         if run_count>1:             #Copy the result of the last algorithm into a temp location             old_output_layer = QgsVectorLayer(temp_output_layer_filename, "input_layer",     "ogr")             Add_result=QgsMapLayerRegistry.instance().addMapLayer(old_output_layer)             write_error = QgsVectorFileWriter.writeAsVectorFormat(old_output_layer,temp_input_layer_filename, "system", QgsCoordinateReferenceSystem(EPSG_Code_for_file), "ESRI Shapefile")             QgsMapLayerRegistry.instance().removeMapLayer(old_output_layer.id())                  #Read in a new file and the file stored in the temp location             input1_layer = QgsVectorLayer(temp_input_layer_filename, "input1_layer", "ogr")             input2_layer = QgsVectorLayer(file_name, "input2_layer", "ogr")                  #Run the binary algorithm and store the result into a temp location or output location             if run_count==len(list_of_files):                 processing.runalg("qgis:union", input1_layer, input2_layer, temp_output_layer_filename)             else:                 processing.runalg("qgis:union", input1_layer, input2_layer, final_file_name)                     QgsMapLayerRegistry.instance().removeMapLayer(input1_layer.id())             QgsMapLayerRegistry.instance().removeMapLayer(input2_layer.id())               path_res = tempfile.gettempdir()     file_list = [shp for shp in glob.glob(path_res + "*.shp")]     Result=(os.path.join(path_res, "result.shp"))          if Select_folder:         clip_function()         binary_alg_wrapper(file_list, Result)


Feature count of current features in Extent based on attribute?

I have one points layer in a cached map service, symbolised with a unique renderer on one field.   e.g. petrol stations, symbolised by the company who owns them. I want to create a fairly basic application in the Esri Javascript API (or potentially with the Google Maps API, and this). One of the main components I need to have on this application is a dynamic legend that shows a feature count for the current extent.   e.g. Shell - 34, BP - 21, Texaco - 7, etc Has anyone seen any code that could get me started on this or can provide some references/pseudocode on how one would go about this?


Problems with WMS QGIS 1.8

Hi I have loaded WMS layers previously but i seem to have hit a brick wall. I am trying to load from either the Geological Survey Irelands website (near the bottom of the page): http://www.gsi.ie/mapping.htm or using the Tellus Project: http://www.tellusborder.eu/Data+Downloads/ I enter the location eg http://spatial.dcenr.gov.ie/wmsconnector/com.esri.wms.Esrimap/GSI_Bedrock_Geology? in server search, the search but neither are presenting layers to load. What am i doing wrong?


I am using Google maps for my application. For getting some custom private map data, I have set up a WMS using Mapserver. How can I add this Custom UMN Mapserver WMS layer over Google API (without using Openlayers)?


What extensions are included in ArcGIS 10.2 Basic Version

I wanted to purchase ArcGIS 10.2. i used its trial version for last 2 month. now i am thinking to have Basic version. My work is based on Spatial analysis. so i dont need 3d and business analyst. for that i wanted know what extension will come with Arcgis 10.2 basic version.?. what are the functions available? and Spatial analyst extension come with ArcGis basic? or i have to purchase spatial analyst separately if anyone have info about this plz tell me.


Creating a list of unique street names from a list of address

I would like to generate a unique list of street names from a list of UK addresses. But obviously there will be a number of streets with the same name. Which is fine, but I want to distinguish between streets that are genuinely different. On AtoZ, the street names on the index is written as:               Bell street, RH1     High Street, RH1       and so on. I have yet to come across a street name with more than one outcode. So based on this,               34, Bell street, RH1 6HH and      36, Bell street, RH1 6HJ      would refer to the same street, but say you had               34, Bell street, RH1 6HH and     36, Bell street, RH2 6HJ      can I be certain that these are two different streets? Thanks


How can we remove the ESRI logo from ArcGIS Viewer for Flex 2.2

Logo layer is set to be visible by default. Is there a way to turn it off? Thank you


How can I transform GCS_North_American_1983 to GCS_WGS_1984?

I posted a thread today describing a shift in some of the NRC data (not all). Apparently, this shift has to do with conversion between Geographic coordinates, latitude and longitude and Transverse Mercator Grid coordinates. Can someone please provide a guide on how to do this transformation in ArcGIS? I would like to transform that data to either `Clarke_1866_Lambert_Conformal_Conic` or `GCS_WGS_1984`. Source Data Coordinate System:               GCS_WGS_1984     Angular Unit: Degree (0.017453292519943299)     Prime Meridian: Greenwich (0.000000000000000000)     Datum: D_WGS_1984     Spheroid: WGS_1984     Semimajor Axis: 6378137.000000000000000000     Semiminor Axis: 6356752.314245179300000000     Inverse Flattening: 298.257223563000030000


Python ArcGIS 10.0 Expanding select layer by location from point (walking selection)

I have an old ArcView 3 tool which I've been tasked with converting to Python. I'm breaking it down into stages, but there's a "walking selection" I can't trap. To put the query in context, the user selects a site (point data) - selection A. All sites within X kilometres of selection A are then selected (selection B), then all sites within X of kilometres of selection B are selected, and so forth until no new selections are made. I've tried to use a While... loop with a Get Count to no avail. The problem I think is because the Get Count will always have a value, so I think I need to do something like Get Count <> Previous Get Count?               #Import ArcPy site package     import arcpy          # Overwrite files     arcpy.env.overwriteOutput = True           try:          # Make a feature layer with all the farms          arcpy.MakeFeatureLayer_management(r"shapefile.shp", "Farms")              # Make a feature layer containing only the farm of interest          arcpy.MakeFeatureLayer_management(r"shapefile_select.shp","Farms_select")              # Apply a selection to farms layer         fcFarms = arcpy.SelectLayerByLocation_management ("Farms","WITHIN_A_DISTANCE","Farms_select", "2000 Meters")              # Apply a further selection to farms layer         fcFarms2 = arcpy.SelectLayerByLocation_management("Farms","WITHIN_A_DISTANCE",fcFarms, "2000 Meters", "ADD_TO_SELECTION")              while getCount > 0:             fcFarms3 = arcpy.SelectLayerByLocation_management("Farms","WITHIN_A_DISTANCE", fcFarms2, "2000 Meters", "ADD_TO_SELECTION")             fcFarms4 = arcpy.SelectLayerByLocation_management("Farms","WITHIN_A_DISTANCE", fcFarms3, "2000 Meters", "NEW_SELECTION")             getCount = int(arcpy.GetCount_management(fcFarms4))          except:         print arcpy.GetMessages()


Creature individual Feature Classes from Features in a Feature Class

I have a feature class called Tracks2012. In the feature class I have 4 fields (MMSI, Ship_Type, Vessel_Name, IMO) as well as the usual Object_ID, Shape and Shape_Length field (as these features are polylines) Essentially what I am trying to do is create an individual Shapefile from each feature using arcpy. I'm not all that familiar with arcpy, and as oppose to cluttering up this question with what I attempted, I'll show the first little part where I set up where I am going to put my feature classes.               >>> import arcpy     >>> arcpy.env.workspace = "c:/GIS/mmsi.gdb"      Once I have this completed, I need to apply a buffer of 10km to each individual polyline shapefile followed by a reclassification of each shape file to have either a 0 or a 1 in relation to the presence of the buffered polyline (which will now be a polygon as a result of the output of the buffer). I am not sure whether or not all of these operations should be maintain within one arcpy script or several. I imagine this all seems quite confusing (the way I have written it up), and I am new to arcpy, so please before down voting my question perhaps suggest alternative to improving its structure. Thanks!


Tiles causing problems with dasharray

I'm having a problem with a stroke-dasharray in Geoserver SLD. When previewing in single-tile mode, everything looks like I want it to, but when I preview in tile-mode I get gaps on the border of the tiles. Left image is tiled, right image is the (preferred) single-tile: ![Tiles](http://i.stack.imgur.com/zuWAZ.png) ![Single- tiled](http://i.stack.imgur.com/5ur03.png) You can see the middle tile here, for some reason it doesn't render the dasharray on either end, even though I have dashoffset set to 0. ![Middle tile](http://i.stack.imgur.com/0nRRH.png) Any idea how to fix this problem? I want to use tiles.


Is it possible to style a part of a label?

My label consists of two different things:               <ogc:PropertyName>user</ogc:PropertyName>     <ogc:PropertyName>room</ogc:PropertyName>      I would like to style the two parts differently, for example have the room number in bold. Is it possible?


Updating a column with line direction in MapInfo

I was wondering if there is a way to update a column in MapInfo with the digitised line direction. I have two TAB files, one is complete lines and the other is sections of the same lines (both have a unique ID). I was wondering if it would be possible to add a column in both files with the direction of the line so I can check if the directions are identical, e.g: ![enter image description here](http://i.stack.imgur.com/aPcKQ.png) Thanks, Alan


statistics by rectangle (raster image)

I have a shp file with some rectangles and I want to extract statistics within these areas from a raster image. I prefer ArcMap 9.3.


What kind of improvement can be done for raster tiles for mobile usage

We generate raster map tiles and deploy them in the server, then we use javascript in the client to request the tiles and re-organize them to map.(like google map and Openstreetmap). However this is the PC solution. And the tile size is 256x256(pixel) with 96dpi. Now we want to build apps in android and iOS. However we found that the same 256x256(pixel) tile will have a different visual effect in devices with different resolution and the tile is almost unreadable in a high-resolution device(The text and labels in the tile is too small). Then I wonder how do you make the tiles readable in the mobile device, is there any way to improve them?


ArcSDE - where to start?

As I long ago figured out that I need more complex query-asking abilities (other than the standard select provided by arcgis), and as I will someday have to admin a small scale SDE database that will import data from my corporate's big and country-wide SDE, I wanted to set up a small workgroup SDE on a standalone computer so I could deal with it a bit, to get the hang of dealing with a sde. But where do I start? I can't find a good manual to start with - The Arcgis 9.3.1 Help(on which I'm working) isn't much help. I know I need to Install SQL server, but I can't understand any of the other steps (creation, import of data and such). Do you know of any good guides and resources to guide me with the arcgis 9.3.1 SDE from scratch? (Oh, and does SDE 10 differ a lot from SDE 9.3.1?)


Converting MapInfo TAB to GeoJSON

Is there an easy way to convert MapInfo TAB files into GeoJSON? I am working on a web-based mapping system which will use GeoJSON to transfer data. I also use MapInfo and would like to be able to convert files from MapInfo into GeoJSON. I have tried the Universal Translator that comes with MapInfo but it doesn't support GeoJSON. If there isn't a way to do it I will make a MapInfo plug-in myself (and probably post it online). I just don't want to waste my time making this if there is already one out there. **Update 2012-03-19:** You can now export to GeoJSON through Quantum GIS.


Why does the buffer tool keep producing ellipsoid when I want a circle?

I want to create a 50 mile radius around a point. I've got the point as a layer. I've run the buffer tool found in Vector -> Geoprocessing Tools and set the distance to 0.5. Every time I run the tool on either a point or polygon I keep getting a ellipsoid 50 miles wide by 70 miles high. I've also tried this on both QGIS and ArcGIS and both produce the same problem. I'm sure I'm over looking something. If someone could point it out what I am missing that would be great. Thank you for your time.