title: LineSymbols not rendering in ArcGIS Runtime for WPF SDK 10.1.1 without AcceleratedDisplay

text: I've posted this over on ESRI's ArcGIS Runmtime for WPF forums but I thought I'd repost here to get some more developer exposure. Any input would be greatly appreciated. Original post here: http://forums.arcgis.com/threads/82721-LineSymbol-and-SimpleLineSymbol-not- rendering-without-AcceleratedDisplay With the AcceleratedDisplay enabled, SimpleLineSymbols and LineSymbols both render correctly. If I turn it off, both fail to render completely. Nothing else has changed such as their geometries or symbologies. Is this a known issue and/or am I doing something wrong? Additionally, SimpleFillSymbols, PictureMarkerSymbols, and TextSymbols seem to work just fine regardless of AcceleratedDisplay settings. That is to say, if I render the symbols as SimpleFillSymbols they will show up (BorderBrush with no FillBrush). Code Behind:                   private void MakeSimpleSymbol()     {                 //Create a simple line symbol                 SimpleLineSymbol sls = new SimpleLineSymbol();                 sls.Color = new SolidColorBrush(Colors.Blue);                 sls.Style = SimpleLineSymbol.LineStyle.Solid;                 sls.Width = 2;                      //Create a polygon object                 Polygon p = new Polygon();                 p.SpatialReference = new SpatialReference(102100);                      //Populate it with some points                 ESRI.ArcGIS.Client.Geometry.PointCollection pointCollection = new ESRI.ArcGIS.Client.Geometry.PointCollection();                 ObservableCollection<ESRI.ArcGIS.Client.Geometry.PointCollection> rings = new ObservableCollection<ESRI.ArcGIS.Client.Geometry.PointCollection>();                 MapPoint wgsPoint = new MapPoint(0,0);                 MapPoint wgsPoint2 = new MapPoint(50,0);                 MapPoint wgsPoint3 = new MapPoint(0,50);                 pointCollection.Add(wgsPoint);                 pointCollection.Add(wgsPoint2);                 pointCollection.Add(wgsPoint3);                 rings.Add(pointCollection);                 p.Rings = rings;                      //Create a graphic and set its properties                 Graphic graphic = new Graphic();                 graphic.Geometry = p;                 graphic.Symbol = sls;                      //Create a graphics layer and add the graphic to it                 GraphicsLayer gLayer = new GraphicsLayer();                 gLayer.Graphics.Add(graphic);                      //Add the graphics layer to the map                 _map.Layers.Add(gLayer);     }      XAML:               <Window x:Class="ESRIWPFTestApplication.MainWindow"          xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation"             xmlns:x="http://schemas.microsoft.com/winfx/2006/xaml"             xmlns:esri="http://schemas.esri.com/arcgis/client/2009"             Title="MainWindow" Height="350" Width="525">         <Grid>             <!-- AcceleratedDisplay is OFF -->             <esri:Map x:Name="_map" UseAcceleratedDisplay="False">             </esri:Map>         </Grid>     </Window>


title: Best Program to use for DEM data?

text: I need to estimate the elevation at various access points of some rivers in South Carolina. I use the elevation information to calculate the average gradient (feet/mile) between two access points. The gradient between points A and B on a river is equal to the elevation difference between A and B divided by the distance along the river between A and B. Usually I get the elevation information from contours on the USGS topo maps. However many South Carolina topo maps have only 20-foot contours when they reach the lower Coastal Plain and within about 20 miles of the coastline. Most states have topo maps with 5-foot contours near the coast. I learned that South Carolina has been working on acquiring accurate elevation data for the entire state. Digital Elevation Models (DEMs) are available for many counties. The data is available here: http://www.dnr.sc.gov/GIS/lidar.html They mentioned I would need some kind on software to read the DEM data and draw the contours with. They suggested that I try Global Mapper. I tried the trial version but got the following error: “ Unable to load overlay gdb -- This release of the GeoDatabase is either invalid or out of date. The ESRI File Geodatabase library does not support GDB files from ArcGIS v9.x or earlier. You will need to convert your GDB data set to a newer version using ArcGIS v10.x or later to make it readable.” Next I looked at ESRI’s ArcGIS Explorer Desktop. I tried to download the trial version, but the larger file was corrupt twice, so I requested a disk be mailed to me. My question to this group -- is there another program available to allow me to see the DEM contours? I do not want to pay for an expensive program considering I need relatively little from it. I would appreciate any advice from users with experience in this area.


title: How to begin with tool development for GRASS gis or QGIS GRASS

text: I wanted to begin with customized tool development in Grass GIS or QGIS grass. How should i begin in to develop the tool. I search a lot in google search engine for begin with tool development in Grass GIS or QGIS grass. Can any one please help me in the tool development process in Grass.


title: Computational most efficient way to convert Cartesian to Geodetic Coordinates

text: As far as I understand the literature there are several ways to convert a set of geoecentric cartesian coordinates to geodetic coordinates   * Bowring, 1976   * Borkowski, 1989   * Lin & Wang, 1995   * Toms, 1996   * ... What would be the most efficient way to convert a cartesian coordinates to geodetic coordinates - most efficient meaning fastest (and most direct) when being implemented on a computer?


title: What is CRS:84 projection?

text: Does anybody know is `CRS:84` in WMS v.1.3.0 `getCapabalities` response is same projection as WGS84 ~~(EPSG:2180)~~ (EPSG:4326)? And is `[westBoundLongitude],[eastBoundLongitude],[southBoundLatitude],[northBoundLatitude]` in `EX_GeographicBoundingBox` respectively equivalent with `[minx],[miny],[maxx],[maxy]` in earlier version of WMS response (1.1.1)?


title: How to convert DTED Level 1 into SRTM DTED

text: We have a GIS tool that handles only SRTM DTED level 1 terrain files. SRTM level 1 terrain files are limited to 60 deg lat. So I need to fill up the gap from 60 lat up to 90 deg lat. I have DTED level 1 data covering the area 60 deg lat and upwards that could be used for that. Is there a way to convert that DTED level 1 data into SRTM data? Very much appreciated, Yvan


title: Indoor positioning system hardware

text: I'm a software developer, dont know much about the hardware specifications and dont know where to start. (Even not sure is that the right stackexchange I should ask this question.) So, I want to know about (1) the cost and (2) document, link, etc for the following specifications: I want to make a project which I will put a chip (?) to a user and track that users 3d location with at least 10 cm accuracy, 10 times a sec, in a particular approximately 25x40 meters place and use that 3d location in my C, C++, python, java application. I read something about triangulation, and it makes sense. However thats all i know. I dont know any device that I can buy to do such.. Can you send me something to read?


title: WMS getMap request with POST and XML

text: It's possible to make WFS queries by POSTing query as XML. Does GeoServer support this also for WMS getMap request? If yes, I'm looking for information on XML syntax for such query. Sounds like a "open documentation and start reading" question, but I haven't found any information. It's possible to POST getMap query with same syntax as regular GET queries have, but I'm looking for XML solution.


title: How to populate SQL options for parameter in ArcGIS 10 tool?

text: I'm trying to create a tool where individuals can select a variety of records using a SQL expression that is then passed to a calculation. I can add a parameter under the tool, and select SQL expression, but when I click on the SQL expression box its not "populated" by potential fields, and therefore unable to get unique values. I thought this might be because the tool needed to know which .shp file it was calling, so I added a higher level parameter to specify the .shp file, but still no luck, the SQL parameter remains entry. Any suggestions would be much appreciated,


title: How to change the WMS Background Color in OpenLayers?

text: I am using geoserver to serve up an image pyramid basemap. I want to have all areas outside of the data extent be black (instead of white). I have the map displayed inside of a div tag and changing the css on that tag or the container (another div) it resides in has no effect. I have checked both the geoserver and openlayers documentation, but have not seen any way to change the background color. Any suggestions on where to continue searching? Do I need to generate a black tiles and serve them?


title: in python, srs.ExportToProj4() returns an empty string, why?

text: I have a shapefile with the following .prj file: PROJCS["NAD_1927_Albers",GEOGCS["GCS_North_American_1927",DATUM["D_North_American_1927",SPHEROID["Clarke_1866",6378206.4,294.9786982]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Albers"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-96.0],PARAMETER["Standard_Parallel_1",29.5],PARAMETER["Standard_Parallel_2",45.5],PARAMETER["Latitude_Of_Origin",23.0],UNIT["Meter",1.0]] When I run the following python code:                   srs = osr.SpatialReference()         if ( srs.ImportFromWkt( prjText ) != 0 ):             return "error"              srsProj4 = srs.ExportToProj4()         sourceProjection = pyproj.Proj( srsProj4, preserve_units=True )      the variable `srsProj4`, as returned from `srs.ExportToProj4()`, is an empty string. When I look here: http://www.gdal.org/ogr/classOGRSpatialReference.html it says that `ExportToProj4()` returns an empty string in the case the srs is in a local (non-geo-referenced) coordinate system. But as far as I understand, the text in the .prj file (above) is a valid, regular geo projection. This code works with other .prj files I've used. Can you tell what's special about this one that's making it fail? Thanks.


title: Can I move selected features in ArcMap with arrow keys?

text: I would like to be able to move features by small (preferably predefined) increments using the keyboard arrow keys? Can this be done?


title: Incorrect date displayed with Esri JavaScript API v.2.6

text: I am working on a site that is using version 2.6 of the Esri Javascript API. I have an info window that displays some date information, and the date is one day behind the date in the database. For example: If I identify a feature in ArcMap, the date returned is 8/15/2007. When I query for the same feature through the Rest endpoint I get 2007/08/15 00:00:00 UTC. When I identify the feature in the web application, the date returned is 8/14/2007. In the javascript code I am formatting the date as follows:               ${SALEDT:DateFormat(selector: 'date', fullYear: true)}      What is causing this, and what can I do to have the correct date returned?


title: How to parse polygons from GML in OpenLayers?

text: I'm reading a GML file from an URL. This file contains polygons. I would like to draw the polygons on a map, where each should have a different color. But all examples I found on the web so far seemed to be outdated and didn't work. Any Hints? Here the parsing (that's the part that already works):                   var theParser = new OpenLayers.Format.GML();         theParser.internalProjection = null;         theParser.externalProjection = null;         theParser.extractStyles      = false;         theParser.extractAttributes  = true;         var features = theParser.read(response.responseText);         if(features)         {              for(var i=0;i<features.length;i++)             {                 var feature = features[i]            //Each feature is an instance of OpenLayers.Feature.Vector                //Each has an attribte geometry, which is of type                 //OpenLayers.Geometry.MultiPolygon and which contains an                //OpenLayers.Geometry.Polygon             }     }


title: How to add an existing feature class to a new feature dataset

text: I have a feature class.               IFeatureClass AOIfc = AOIfws.OpenFeatureClass("BrigadeAOI");      I need to be able to put this feature class into a feature dataset that is newly created. The following line of code doesn't work for me.               IFeatureDataset AOIfd = AOIfc.FeatureDataset;      How do I accomplish this? Thanks!


title: Install Shapely with GEOS on heroku

text: I have a web mapping app running on Heroku. My app is written in python. Using pip I was able to install various pure python libraries that were needed. I would now like to install Shapely too. This is simple enough (with pip) but I am unsure how to install the geos pre-requisite within the virtualenv (venv) directory. I would much appreciate any pointers. Thanks, John


title: What is wrong with how I am using ST_Within?

text: I am learning postgis 2.1. I created a postgis table from a shapefile and imported it into the DB. It made a table like this:               db=# \d districts           ...          geog       | geography(MultiPolygon,4326) |       I also created a table "callsforservice" importing tabular data from a CSV, which I then converted into points with a SRID = 4326.               db=# update callsforservice set latlongpointdetails = ST_SetSRID(latlongpoint, 4326);      I want to know which `callsforservice.latlongpointdetails` are included in which `districts.geogs` I tried with this query, which takes very, very, very long time. The limit 1 is only to shorten the time the query is running. After 2 min or so I kill it:               select * from callsforservice, districts where ST_Within(geometry(districts.geog), callsforservice.latlongpointdetails) limit 1;       If I add an additional condition to the where clause (`where districts.district='1-1'`) -- this query returns immediately with 0 rows, even though districts does have a row with this ID. Do I have the syntax right for finding out which calls came in what districts? The docs for ST_within say "the source geometries must both be of the same coordinate projection, having the same SRID." If I examine the two tables with the `\dt` command then the geog column in districts shows type `geography(MultiPolygon,4326)` but the `latlongpointdetails` column in callsforservices only shows type geometry -- without parenthesis like in the districts table. Maybe this is the reason for the failures of ST_Within? Where am I going wrong?


title: About proxy in javascript api

text: I am new about Arcgis javascript API. While I am doing a sample application, I confused. I am analyzing this sample application. This application is using the http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Louisville/LOJIC_LandRecords_Louisville/MapServer/0 as FeatureLayer. And using a proxy                   esri.config.defaults.io.proxyUrl = "/proxy";         esri.config.defaults.io.alwaysUseProxy = false;      I created a new visual studio project and proxy pages. And I changed FeatureLayer and Base Map. Program is working. But I opened Mozilla Firebug and traced the traffic. My Application using multiple Http POST data, while zoom in. But Arcgis application Only 2 GET **This is my custom application screencast.** (This is screenCast of my application i.stack.imgur.com/F1C17.png) This is Arcgis Sample Application Here (help.arcgis.com/en/webapi/javascript/arcgis/samples/widget_measurement/) (This is Arcsgis Applcation screencast i.stack.imgur.com/muZV9.png)


title: Geographic Transformations in Alaska Between NAD83 and WGS84

text: I have done a ton of research on the most accurate transformation to use in mainland Alaska to convert between the NAD83 and WGS84 datums, but I just can't seem to find a definitive answer. ArcMap gives 3 options for my area of interest:   1. WGS_1984_(ITRF00)_To_NAD_1983      * This one I am really unsure of and I have not found much information on when it should be used. I know it has something to do with epochs or realizations of datums, but I am in over my head on this one. Can someone please enlighten me?   2. NAD_1983_To_WGS_1984_5      * This is the one I think might be most applicable, but I am unable to find any definitive information to say whether or not it should be used for Alaska.     * Esri says it is applicable for Alaska, but I can't find confirmation of this anywhere else     * EPSG doesn't list this transformation when I search with coordinates in Alaska (or at all for that matter)   3. NAD_1983_To_WGS_1984_1      * I am fairly certain that this is a "null" transformation that doesn't actually do anything, and would therefore be least accurate. Especially in Alaska where everything shifts with time. I should also note that we are using NAD83 (NSRS 2007) in our geodatabase and converting GPS survey data from WGS84 to NAD83 (NSRS 2007) with a custom transformation developed by a geodesist for our area of interest which shifts everything by a meter or so in x and y (the Esri transformation for this is a null transformation that doesn't actually do anything). What I am trying to do is load legacy/reference data that is currently in NAD83 by first transforming to WGS84 (my question here) and then transforming to NAD83 (NSRS 2007) via our custom transformation. I need to define a definitive procedure for the project team to follow and it will be used on potentially hundreds of layers so I need to get it right.


title: geojson from st_dwithin query

text: Im trying to query all the objects within 100m and I want the result to be in geojson format not list of geojson. Im using PostgreSQL 9.1, PostGIS 2.0. Tried this query:               SELECT ST_AsGeoJSON(geocolumn) FROM geotable WHERE ST_DWithin(geocolumn, 'POINT(1000 1000)', 100.0);      but the result is a list of geojsons. Is there a way to convert the result of my query to a single geojson?


title: Why new feature disappears when editing polygon?

text: I try to create internal polygon within existing polygon. After I created a ring, I add a new feature inside the ring. Suddenly, this new feature disappears. I confirmed that this feature went to the origin( 0,0). The polygon and project file are used UTM as projection. If somebody knows, why it is happening, let me know the reason and solution. QGIS:1.8.0 OS:Windows 7


title: arcgis 10.0 extensions not getting authorized after reinstall

text: I have reinstalled ArcGIS 10.0 and trying to authorized ArcEditor and other extensions. I am able to authorized ArcEditor whereas extensions are not getting authorized. Any idea? Thanks in advance.


title: ArcGIS Custom Toolbox crashes ArcMap 10.0

text: I have prepared Custom ArcToolBox using ArcObjects in java. It was working fine before though was giving NullPointerException but this something java related error and that is fine. Later I have reinstalled ArcGIS10.0 and trying to run the same tool. The moment I start the process my ArcMap10.0 is getting closed. Interesting part is if I remove my toolbox related jar from /Desktop10.0/java/lib/ext folder it is working fine again. Does anyone has any idea what could be the issue? Thanks in advance.


title: How to use multiple wildcards within an iterator in ArcGIS ModelBuilder?

text: I would like to know if it is possible to use multiple wildcards to select specific shapefiles when using the Iterate Feature Classes tool in ArcGIS modelbuilder. E.g. I use this to select all shapes beginning with T but wonder how I would ask the iterator to select shapes beginning with R,S & T for example. ![Iterate Feature Classes](http://i.stack.imgur.com/k0bXl.png) Any ideas much appreciated.


Is Python Toolbox (.pyt) newly introduced in ArcGIS 10.1? Was it there in ArcGIS 10.0? Thanks


title: Update WMS layer with time dimension

text: I have a time series of MODIS images that I serve as WMS with the time dimension option. I successfully created the store following this tutorial and I can load the layers and set the time variable. The question is: How should I update this store / layer when I receive a new MODIS image? Should I just add a new record to the Postgis table that controls the store? Or do you use the REST API? I tried using the example given here for updating the imagemosaic but I get an error saying POST is not an allowed method (shown bellow). I've tried allowing for post to work but could not succeed.               curl -v -u user:my_password -XPOST -H "Content-type: text/plain" -d "file:///home/davi/novo_teste/rs.A20140101_IVP_NDVI.tif" "http://localhost:8181/geoserver_lab/rest/workspaces/maranhao/coveragestores/banana/external.imagemosaic"     * About to connect() to localhost port 8181 (#0)     *   Trying ::1... connected     * Connected to localhost (::1) port 8181 (#0)     * Server auth using Basic with user 'davi'     > POST /geoserver_lab/rest/workspaces/maranhao/coveragestores/banana/external.imagemosaic HTTP/1.1     > Authorization: Basic ZGF2arTolY2hpbmcrxMw==     > User-Agent: curl/7.19.7 (x86_64-redhat-linux-gnu) libcurl/7.19.7 NSS/3.14.0.0 zlib/1.2.3 libidn/1.18 libssh2/1.4.2     > Host: localhost:8181     > Accept: */*     > Content-type: text/plain     > Content-Length: 54     >     < HTTP/1.1 405 Method Not Allowed     < Allow: GET, PUT     < Date: Wed, 12 Mar 2014 12:18:51 GMT     < Server: Noelios-Restlet-Engine/1.0..8     < Transfer-Encoding: chunked     <     * Connection #0 to host localhost left intact     * Closing connection #0


title: detect signal where delete polygon in qgis?

text: I am working with QGIS and i want when i delete a polygon and then i commit this action I detect the signal of this operation and use it in python script, I tryed this but dont work:               def delete():       sql= ("delete from zoneshumidesdelimitation where idzoneshumideslocal  = " + str(featureid))      executeQuery(sql)      layer.committedFeaturesRemoved.connect(delete)       Could someone help me please!


title: ArcGIS Javascript API - Looking For Styled Navigation Plugin/Widget

text: I am looking for a fancy navigation plugin to use with the ArcGIS Javascript API. Something along the lines of the Navigation Widget (see image below) available in the Flex Viewer would be great. Anyone have any ideas/suggestions? ![Navigation Widget](http://i.stack.imgur.com/z6lrD.png)


title: Programatically adding popups into OpenLayers map

text: I'm developing a web service and want to inset multiple popups dynamically to a web page according to the request. What I did was this               <script> //injected programmatically     var popups=[];     popups.push('11.1117177 46.077171   Title One   Description one,Second line.    24,24   0,-24');     ....     </script>     ....               function showPopups() {         var vectorLayer = new OpenLayers.Layer.Vector(         "Popup layer", {         style : reconstructedTrajectoryStyle         });              var parser = new OpenLayers.Format.Text();         var src = new OpenLayers.Projection('EPSG:4326');         var dest = new OpenLayers.Projection('EPSG:900913');              var s, wktLine, features, reprojectedFeatures;         for (s = 0; s < popups.length; s++) {         wktLine = popups[s];         var wktLine1="lat   lon title   description icon    iconSize    iconOffset \n"+wktLine;         features = parser.read(wktLine1);         reprojectedFeatures = features.geometry.transform(src, dest);         vectorLayer.addFeatures(new OpenLayers.Feature.Vector(         reprojectedFeatures));         }              osmmap.addLayer(vectorLayer);         osmmap.zoomToExtent(vectorLayer.getDataExtent());     }          $(document).ready(function() {     initialize();     ....     showPopups();          });      What I got at the end was not successful. I got the error on the browser console saying Uncaught TypeError: Cannot call method 'transform' of undefined ز What am I doing wrong? Should I use another easier technique?


title: Retrieve common paths of polylines ArcObjects

text: I have a polygon featureClass and a polyline featureClass. I now want to handle all geometries within the latter that are contained in the boundaries of polygon features, that means they shell simply touch the polygons boundary. This accords to the simple "select features from that share a line segment from " -operator within the "select by location"-tool. But how to achieve this with ArcObjects? I already tried the following:               IPolyline reshaper = ...;       ITopologicalOperator topo = (ITopologicalOperator) singlePolygon;       IGeometry boundary = topo.Boundary;       IRelationalOperator relOp1 = (IRelationalOperator) boundary;       if (relOp1.Contains(reshaper)) {...}      I also worked with the horrible Relation-Method which is intended to use on more complex queries. I tried it with `relOp1.Relation(reshaper, "g1 contains g2")` but simply received a meaningless HRESULT-Error of kind `ARCWEB_E_INVALID_RING_CALCULATION_TYPE` if I looked up the errorlist from ArcoBjects correctly (which I cannot ensure...). So what now?


title: How do I solve Gdalwarp error 'too many points failed to transform' for remapping Geostationary to Lambert conformal?

text: I am trying to remap from Geostationary to Lambert conformal using gdalwarp. My input data is in netcdf, and is in geographic coordinates (degrees) and I would like to output the remapped data to netcdf. I have created a corresponding vrt file for the input netcdf data. Gdalwarp will output the netcdf file, but the output data is all zeros and I receive the following error:               Creating output file that is 5120P x 5120L.     Processing input file netcdf.vrt.     ERROR 1: Too many points (441 out of 441) failed to transform,     unable to compute output bounds.     Warning 1: Unable to compute source region for output window 0,0,5120,5120, skipping.     0...10...20...30...40...50...60...70...80...90...100 - done.      I attempted the following command:               /usr/bin/gdalwarp -s_srs "+proj=geos +h=35785831 +lon_0=-75 +x_0=-0.151844 +y_0=0.151844 +a=6378140 +b=6356754.99999591 +units=degrees +no_defs" -t_srs "+proj=lcc +ellps=clrk66 +a=6378137 +b=6378137 +e=0.0818191910435 +lat_0=24.9999 +lon_0=-95 +lat_1=24.9999 +lat_ts=25.0001 +units=meters +no_defs" -te -1952976.3246 -828316.5944 3248431.6754 4373091.4056 -of netCDF -geoloc -overwrite -r bilinear -ts 5120 5120 netcdf.vrt out.nc      Can gdalwarp remap from geographic coordinates to projected? Or do I need to translate geographic to projected first? Also, can gdalwarp read projection information straight from netcdf or do you NEED to write to .vrt first? Here is what gdalinfo outputs from the input file: (it is a GOES 13 file from CLASS)               Size is 512, 512     Coordinate System is `'     Metadata:       NC_GLOBAL#Conventions=CF-1.4       NC_GLOBAL#Satellite Sensor=G-13 IMG           NC_GLOBAL#Source=McIDAS Area File     Subdatasets:       SUBDATASET_1_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":auditTrail       SUBDATASET_1_DESC=[3x80] auditTrail (8-bit character)       SUBDATASET_2_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":data       SUBDATASET_2_DESC=[1x665x2036] data (32-bit floating-point)       SUBDATASET_3_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":lat       SUBDATASET_3_DESC=[665x2036] lat (32-bit floating-point)       SUBDATASET_4_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":lon       SUBDATASET_4_DESC=[665x2036] lon (32-bit floating-point)     Corner Coordinates:     Upper Left  (    0.0,    0.0)     Lower Left  (    0.0,  512.0)     Upper Right (  512.0,    0.0)     Lower Right (  512.0,  512.0)     Center      (  256.0,  256.0)      And additional gdal info on the data variable:               Driver: netCDF/Network Common Data Format     Files: goes13.2013.100.174518.BAND_04.nc     Size is 2036, 665     Coordinate System is `'     Metadata:       data#coordinates=lon lat       data#long_name=0-255 Brightness Temperature       data#type=VISR       NC_GLOBAL#Conventions=CF-1.4       NC_GLOBAL#Satellite Sensor=G-13 IMG           NC_GLOBAL#Source=McIDAS Area File       NETCDF_DIM_EXTRA={time}       NETCDF_DIM_time_DEF={1,4}       NETCDF_DIM_time_VALUES=1365615900       time#long_name=seconds since 1970-1-1 0:0:0       time#units=seconds since 1970-1-1 0:0:0     Geolocation:       LINE_OFFSET=0       LINE_STEP=1       PIXEL_OFFSET=0       PIXEL_STEP=1       SRS=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"]]       X_BAND=1       X_DATASET=NETCDF:"goes13.2013.100.174518.BAND_04.nc":lon       Y_BAND=1       Y_DATASET=NETCDF:"goes13.2013.100.174518.BAND_04.nc":lat     Corner Coordinates:     Upper Left  (    0.0,    0.0)     Lower Left  (    0.0,  665.0)     Upper Right ( 2036.0,    0.0)     Lower Right ( 2036.0,  665.0)     Center      ( 1018.0,  332.5)     Band 1 Block=2036x1 Type=Float32, ColorInterp=Undefined       NoData Value=9.96920996838686905e+36       Metadata:         coordinates=lon lat         long_name=0-255 Brightness Temperature         NETCDF_DIM_time=1365615900         NETCDF_VARNAME=data         type=VISR      Any help will be greatly appreciated!


title: Where to download ASTER scene boundaries?

text: I am looking for a shapefile of ASTER scene boundaries for the Kansas and Oklahoma region. The attached screenshot from GLOVIS is what I'm after, except I need polygon footprints with relevant attributes such as scene ID and year. Is such a product available for download? ![enter image description here](http://i.stack.imgur.com/kyaG1.jpg)


title: GeoExt: Turning off layer title in legend

text: I am using GeoExt for my web mapping project. My web map includes a legend generated by GeoExt. Is it possible to turn off the layer title for a specific layer in legend? I am aware that there is a showTitle property for GeoExt.WMSLegend. However setting this property to false will turn off the layer titles of all the layers in the legend. Is there any way to turn off the layer title of just one single layer?               // Add a legend panel     legendPanel = new GeoExt.LegendPanel({         map: map,         title: 'Legend Panel',         renderTo:'description',         border: false,         autoScroll: true,         defaults: {             showTitle: false,             style: 'padding:5px',                        baseParams: {                 FORMAT: 'image/png',                 LEGEND_OPTIONS: 'forceLabels:on;fontName=Verdana;fontSize:11'             }         }     });


title: What is a good metric for mapping accuracy

text: I have recently received a dataset that supposedly mapped the altitude a given region (a DEM). I would like to determine the accuracy of the new map by comparing it with a previously generated DEM. I am not sure what kind of metric I can use to ensure that the performance is good enough. (Good enough means that the horizontal deviation of a feature can not be off too far and that the vertical deviation can not be off too far) One way to attack this would be to descritize the space into small chunks and then look at the distribution of the vertical residuals. If that block has vertical residuals that are less than the vertical deflection then I can assign that block a TRUE value and if it does not I could assign the block a FALSE value. I can then simply calculate the mean of the TRUE/FALSE values across all the blocks to get an overall value such as 67% or something which would tell me how good my mapping performance was. There are huge problems with this method:   1. It does not seem to factor in the horizontal deflection requirement   2. It seems as though it would violate the iid (identical independently distributed) assumption normally required of the mean calculation. This is the same as stating that the individual boxes would be correlated.   3. There does not seem to be a defensible method to choose how big the descritization box size should be. An alternative could be that I simply place special markers with known 3D locations and then simply measure their vertical and horizontal displacement. This seems like it could share some of the same problems as the earlier method. (how many markers should I place? are the measurements correlated if they are too close?) Before running off and doing something silly, I wanted to see if the community has a golden metric or method that has good statistical properties. Something that is continuously measured across the DEM, takes into account both the vertical deflection and the horizontal deflection, and I can look at the distribution and make a defensible argument? Please let me know if this question is unclear and thanks for any suggestions. P.S. Could somebody could comment on the need of the two individual mapping accuracy levels (one for vertical deflection and one for horizontal deflection)? It seems to me that a horizontal mis-mapping would be reflected in a vertical deflection. Is it standard to have two separate requirements for mapping?


title: Adding feature to remote geodatabase takes a long time

text: I have a program that I must add some features to a `feature class` that is in a sql server remote database. But when I want to open connection and add them, it takes a long time, about one or two minutes and it is not good. I tried to add them manually using SQL statements and I cannot do that if I enable geodatabase on my database. So if I want to add features manually, I should not enable geodatabase. On the other hand, I want to edit shapes using arcobject so I must enable geodatabase. Here is the question:   * If I dont enable geodatabase, I can load and add features very fast and it is great but I cannot edit them using arcobject.   * If I enable geodatabase I can edit features but it is too slow. What are your suggestions? Is there an alternative way? Are there any libraries that I can use to edit shapes directly from none geodatabase?


title: A library to edit geometry in C#

text: For my windows application, I am looking for a method to edit vertexes and geometry of features that are stored in a SQL Server database. I just want a simple function for dividing polygons into 3 parts and so on. Are there any open source library for these functions?


title: How to convert to web mercator projection using an open source tool?

text: I've got a KML file which I need converted to a web mercator projection. This file will go into Tilemill. What are some good open source tools to convert files to the web mercator projection?


title: How To Turn Off Map Drawing When Using Python in ArcMap?

text: I have a Python script that creates layers in ArcMap. The script works fine, but when a new layer is created a map is always automatically drawn, which makes the script take longer to process. I really don't need the layers to be drawn on the map. Is there any way to set my script to not draw the layers after they have been created? During a search I found a similar question: How to disable and enable drawing in ArcMap 10 using Python? Attempting that solution did not work; "No module named SendKeys" I'm using ArcGIS 10.2 for Desktop.


title: ArcGIS JS API Relationship Queries

text: Has anyone ever put the results of a RelationshipQuery() into a map.infoWindow that's been replaced with a dijit.popup? Here's the workflow: User clicks a point, and all identify results are put into popup. <--this part is fine I'm struggling with calling a relationshipquery, then putting the results into the infoTemplate. The function findRelatedRecords performs 3 relationship queries and assembles an array of objects that are related to the single point. I cannot get this data back into the executeIdentifyTask function. The function getWindowContent takes the relationshipquery array and converts it to html and constructs a dijit.tabcontainer & 2 contentpanes (tabs). It return the tabcontainer.domNode. Any help? I'm really struggling!!               function executeIdentifyTask(evt) {         facSel.clearSelection();         identifyParams.geometry = evt.mapPoint;         identifyParams.mapExtent = map.extent;              var deferred = identifyTask.execute(identifyParams);              deferred.addCallback(function(response) {                // response is an array of identify result objects               // Let's return an array of features.           return dojo.map(response, function(result) {             var feature = result.feature;             console.log(result);             //feature.attributes.layerName = result.layerName;             if(result.layerName === 'parcel_area2'){               //console.log(feature.attributes.site_code);               var template = new esri.InfoTemplate("", "${site_code} <br/> Owner of record: ${facil_loca}");               feature.setInfoTemplate(template);             }             else if (result.layerName === 'logistics2'){               var template = new esri.InfoTemplate("", "Parcel ID: ${Type}");               feature.setInfoTemplate(template);             }             else if (result.layerName === 'query2'){                 findRelatedRecords(feature,evt) //<--get relationship query results here                     var template = new esri.InfoTemplate("", getWindowContent(data,evt.graphic));                 feature.setInfoTemplate(template);                            }             return feature;           });         });         // InfoWindow expects an array of features from each deferred         // object that you pass. If the response from the task execution         // above is not an array of features, then you need to add a callback         // like the one above to post-process the response and return an         // array of features.         map.infoWindow.setFeatures([ deferred ]);         map.infoWindow.show(evt.mapPoint);            }


title: Symbology rotation

text: I am using arcgis 9.3, windows xp and want to know if it's possible to have symbology in my polygons not rotate with data frame when I switch between mapsheets in my map series? Any help would be greatly appreciated Oh yeah -I do NOT have MapAtlas or PLTS. I DO have ds mapbook installed, and that is where the trouble starts. Thanks!


title: qgis rgb raster layer: min/max-values definied by user not working?

text: I tried several times to define the min/max-values for the rgb-display under properties>style: Benutzer min/max-Werte (user min/max-values) but it is not working. By chance it worked once but now it is not working anymore. The other functions as to invert the color representation or the selecton of the channel are working, hm. The raster file are created by the following process: Scanned as an pdf>imported in photoshop cs5 (modus rgb)>saved as (in photoshop) tif (color profile not embedded, respectivel sRGB IEC61966-2.1 not activated). Does anybody have an idea where the problem is? Thank's in advance! phil


title: How does one publish a map to the web?

text: Just starting out with ArcGIS server at work, I've had no prior experience but I been given the responsibility to display our map on the web. I have a ArcGIS server running on a server and I am using ArcMap 10 on my desktop. I created a Folder in the server where I saved the map document. I can add a new service using this document and I'm given a url when I do so. But when I go to that URL I get this error: `Server Error in '/ArcGIS/Services' Application. Runtime Error` And, when going to the same url from the server I get: `Server Error in '/ArcGIS/Services' Application. No Content` I'm using ArcMap10 and ArcGIS server 10


title: Determine if a polygon intersects itself in OpenLayers

text: I let the users draw polygons and then I try to save them out to the SQL Server 2008 SQLGeography field, but it is giving me an error: > Error 24200: The specified input does not represent a valid geography > instance. I think the issue is because the polygon intersects itself and SQL doesn't accept this behavior. Is it possible to do a check in OpenLayers to see if the WKT from the polygon would be accepted by the SQLGeography.STGeomFromText function?


title: select feature info on user defined polygon

text: Ok so I have a Openlayers App which is displaying few of the features and I now want to give the user the freedom to create a Box or polygon. All the features which are overlapped by this user defined box, I wish to export the data and then display it onto a table in the html.   1. Is there any OL example to such application requirement?   2. Which would be better to extract the feature info, WFS or WMS?  I know it can be done in GMaps but I have OSM and GMaps both as base maps so I require a OL functionality.


title: Is Data Interoperability Extension for ArcGIS Server 10 required to extract data in .dgn format?

text: I am trying to create a Geoprocessing Task in my .mxd to in order to extract data to .SHP and .DGN formats which will be utilized in a data extraction tool for a web application. Right now, there is a version conflict between ArcMap 10 / ArcGIS for Server 9.3 which will be solved next week the server is upgraded to 10. Is the Data Interoperability (DI) extension required for ArcGIS Desktop to create the GP task? More importantly, is the DI extension/license required for ArcGIS Server 10 for the .DGN extraction tool to work on the application?


title: How to recreate Mapinfo labels in QGIS?

text: I try to convert a mapinfo tab-file containing textobjects wich are positioned and rotated in mapinfo into a shapefile or something similar without loosing all information. All information is located in the mapfile and the attribute table is empty. Just opening the tab-file in Qgis gives a point with an id containing NULL - the text-value lies inside the geometry and the same happens when trying to open the file as an exported mid/mif(thought the point geometry is not really the same). When opening the mif-file in a texteditor it looks like               "label"     x1 y1 x2 y2     Font (font and color info)     Angle 0.8     Justify Center      where the x1y1 and x2y2 is the bbox (bottom left and top right points) of the text, angle is the bbox rotation and the text itself is contained inside the "label". Can I somehow create a polygon from the bbox values and the angle and attach the label and the angle as attributes somehow?


title: Can we direct connect to PostgresSQL without ArcSDE

text: I have postgreSQL installed on my machine.My machine already has ArcGIS 10. No Sde installed. Can we create geodatabase in PostgreSQL with out installing ArcSDE? * * * Jay Cummins is right. I am intended to use create esri geodatabase. Can we do this with direct connect by simply installing postgreSQL on same/remote machine? Can we connect to Postgres with ArcCatalog and create featureclasses. Can we load raster in it? I have seen in many articles mentioning post- installer to setup the ArcSDE objects in the database even for direct connect. Should we need ArcSDE?


title: How to compute a new attribute based on changes in another attribute?

text: I am attempting to classify a set of gps time-encoded point data into behaviors based on different attributes. I have created an attribute that is 0 for home and 1 for away based on location, and now want to number the trips away from home (a set of points `01111111111110` would be one trip because it started and ended at home). I have added the attribute field that will have the trip numbers, but don't know how to calculate the field so it is based on the home/away field. **Here is an example** of the GPS data (using "*" to indicate irrelevant information and simply indexing times as 1, 2, etc.), the "Home/Away" indicator described above, and the desired trip indicator, "Trip", which I need to compute:               Time Lat Lon Home/Away Trip        1   *   *         0    0        2   *   *         1    1        3   *   *         1    1     ....       12   *   *         1    1       13   *   *         0    0       14   *   *         0    0       15   *   *         1    2       16   *   *         1    2     ....        34   *   *         1    2       35   *   *         0    0       36   *   *         0    0       37   *   *         1    3     ....      My data set is too large to manually go through and number each trip in the attribute table, so is there any way to calculate the field based on how the home/away attribute is ordered and each "clump" of away points is designated as a trip? UPDATE: these are the bare bones of what Python code might look like (I'm not experienced with code). ### Expression:               trip = Reclass(!home!)      ### Codeblock:               def Reclass(home):         if (home = 0):           return 0          elif (home = 1 and lastValue = 0):           return _(incremental numbering?)_         elif (home = 1 and lastValue = 1):           return lastValue        _**UPDATE:_**   After using matt wilkie's recommended script I've made some changes so that my first trip is number 1, my second is 2, etc etc. Here is the code modified from matt's:               import arcpy     rows = arcpy.UpdateCursor("test2")          trip = 0     for row in rows:         if row.home == 0:             prev = row.home             row.TRIP = trip             rows.updateRow(row)              elif row.home == 1 and prev == 0:             trip += 1             prev = row.home             row.TRIP = trip             rows.updateRow(row)             rows.next()              elif row.home == 1 and prev == 1:             prev = row.home             row.TRIP = trip             rows.updateRow(row)             rows.next()              row.TRIP = trip         rows.updateRow(row)               del row, rows      Then I just select for home = 0 and calculate my trip field back to 0. Neatly ordered trips. Thanks a ton everyone!


title: How to create a vertical scale bar in ArcMap?

text: Wondering if anyone out there knows how to create a vertical scale bar in ArcMap? With text I can adjust the angle. Is this possible with a scale bar as well?


title: What is the maximum number of vertices for polygons & polylines in PostGIS?

text: What is the maximum number of vertices for polygons and polylines (including multipart) in PostGIS? Thanks


title: Does Metadata for ArcGIS Python script tools in custom toolbox work?

text: * I've created several python scripts and created tools from them.   * I have a custom toolbox that contains all of my tools.   * I add this toolbox to the ArcToolbox window in ArcCatalog.   * I want to create metadata and set parameter explanations for each tool. I've noticed that, when I navigate to my toolbox using the ArcToolbox window, I cannot edit metadata by highlighting the tool in the toolbox window and choosing the 'Description' tab in the catalog window. I can, however, right click in the tool and choose 'Item Description', which opens another window and I can edit metadata in there. Next, I go to the Catalog Tree Window and navigate to the location of the toolbox. In this situation, I can highlight the tool and choose the 'Description' tab and edit the metedata there. This, however, does not translate the edited metadata back to the same tool in the ArcToolbox window. If I open the tool in ArcToolbox window and highlight a parameter, I still get a description message of 'No description available'. So, it's like I have to edit the metadata twice. Once in the Arctoolbox window and then again in the Catalog Tree Window. I've noticed this about the Description tab in ArcCatalog. It only honors selection in the Catalog Tree Window. If I a have tool selected in the ArcToolbox window, the description is still set to whatever is highlighted in the tree window, not the toolbox window. Is this how it's supposed to work? I'm OK with only having to be able to edit metadata in the Catalog tree window, but the fact that metadata does not translate back to the same tool in ArcToolbox window is odd to me. Am I doing something wrong?


title: Publishing Maps to Mobile Devices and the Web

text: I have a bunch of agricultural information stored in a local postgis database. We use qgis to manipulate/analyze this information and to print output. I would like to publish this information (both **raster and vector** layers) so that it can be loaded as a layer on **both mobile devices and web browsers**. The information to be displayed needs to be a **function of the user that is logged in**. It also must be able to **interact with the gps on end users mobile devices**. I would like to host all of the data to be published on the web offsite because my upload bandwith is terrible. In the past I have saved vectors as .kml files and loaded them into google maps for users, but this strategy seems to have problems when **displaying large datasets**. I'm comfortable doing a large amount of development work but am also willing to pay for solutions that save me time and headaches. After researching this project for a couple of days I've found that there are about a million ways to do this. Any help narrowing down my options would be much appreciated.


title: python script error- with GetParameterAsText

text: I have a seemingly simple code which iterates through a list of feature classes and calculates an ID field in the attribute table that is equal to the first 5 characters of feature class name. The script works fine if I hard code the workspace. As soon as I change it to "Get Parameter as Text" I receive this error: ERROR 000539: Error running expression: expression : name 'expression' is not defined Failed to execute (CalculateField). The script is:               import arcpy, sys, array, time, datetime          #Get feature dataset     inws = arcpy.GetParameterAsText(0)          #Set geoprocessing environments     arcpy.env.workspace = inws     arcpy.env.overwriteOutput = True          #Iterate through the list of feature classes          fcList = arcpy.ListFeatureClasses()          for fc in fcList:         print str("processing" + fc)          #Define field and expression     field = "SID"     expression = str(fc[:5])          #Create a new field name and expression     arcpy.AddField_management(fc, field, "TEXT")     arcpy.CalculateField_management(fc, field, "expression", "PYTHON")      does anyone know what is going on here? Thanks!


title: pour point user click input over map for watershed delineation

text: I have a drainage density map and i want to interactively retrieve watershed by clicking a point on map. Once I zoomed to an area of interest and click the highest accumulation point of my choice then i get in return watershed only encompassing that outlet. I am wondering if anyone had already done that manually or automatically. Any suggestion to do it will be helpful. Thanks.


title: How to label using just numbers from start of text field?

text: In ArcMAP 10.2, I am trying to format the labels on my map. If the label is 1402 Bowie Road But I want only to display numbers: 1402 What is the expression to create this? * * *


title: Posting UserCredentials with AsyncTask of ArcGIS for Android?

text: My problem is that I want to get JSON data from ArcGIS for Server For example: http://xxxx.com/arcgis/rest/services/Map/MapServer/legend?f=pjson I can get and read JSON data from this URL but if ArcGIS for Server has a password I couldn't so , must use UserCredentials but every sample example use UserCredentials for identity task or querytask but I read data from URL. how can I post UserCredentials with my AsyncTask ?


title: Automated testing for GIS apps

text: There are several commercial and opensource tools for regression automation. While many of them could be employed to test GIS apps, they fall short, at times, dealing with the unique characteristics and use-cases of GIS apps and their imagery. What are some tools for validating GIS apps? I came across one such tool called Eggplant. Check out their video link in the bottom. They claim to not only compare data sets but also compare images. Seems interesting! Are there other such commercial or opensource tools we can evaluate? Thanks,   Ramesh


title: How to export a shapefile through the Wikimapia API?

text: I am trying to download a kml file from Wikimapia containing information on quaters (districts), the very same information that is nicely available on Wikimapia. What I want to get at the end of the day is a shapefile containing as many polygons as there are districts in Milan to use in my GIS. This is my URL for the query I pasted in my web browser.               http://api.wikimapia.org/?function=box&key=[mykey]&lon_min=8.936&lat_min=45.282&lon_max=9.421&lat_max=45.652&format=kml&category=quarters      From this I get a tiny kml.part file, that I then rename .kml. The file doesn't open in Google Earth or QGIS. If I open it with a text editor I can just few lines of xml, but no information concerning my geometries... Am I doing something logically wrong?


title: How to import a map consisting of pictures

text: I want to built some sort of Fantasy Map with Quantum GIS. The whole map is built out of 128x128 single pictures, 16384 in total. It should be a flat map. How do I import these pictures and align them in a grid into QGIS?


title: Create vector layer from feature collection in OpenLayers

text: I have used a WPS process `gs:UnionFeatureCollection` to compute the geographical union of two collections of WFS vector features. The process returns a collection of features which I intend to display on the map, first, by constructing a layer and adding it to the layer switcher.                wpsClient.execute({         server:'local',         process:'gs:UnionFeatureCollection',         inputs:{             first:features1,             second:features2         },         success:function(output){             var features = output;             localStorage.setItem("Feature", JSON.stringify(JSON.decycle(features)));             window.location.reload();         }     })      This is the code to execute the process. I am saving the resulting feature collection to browser LocalStorage and refreshing the browser window. During the loading a new vector layer is being created and added to an array of layers that will be added to the map.               if(localStorage.getItem("Feature") != null){         var features = JSON.retrocycle(JSON.parse(localStorage.getItem("Feature")));         var vectorLayer = new OpenLayers.Layer.Vector("Union Layer");         vectorLayer.addFeatures(features);         layers.push(vectorLayer);     }      Here is my problem, the layer does appear on the layer switcher but not on the map even when enabled. Is defining features for the layer sufficient to display it or do I need to add it to the geoserver? (I'm pretty new to web mapping but not so much to javascript)


title: How can I find out Persional GeoDatabase version?

text: I am working ArcMap software. How can I find out Personal Geodatabase Version? I have client file. But I want to find out which version they have to used. I am looking for your guidance


title: How do you remove isolated roads from a bigger network?

text: I'm trying to perform an OD Cost Matrix calculation on a set of points using a road network from OpenStreetMap data. The network contains lots of roads that don't connect to the rest of the network. In the case where a point snaps to such a place on the network, it can't reach any other sample points. How do I remove these sections of road? Or if that's not possible, how could I systematically identify those areas and connect them to the rest of the network?


title: Adding multiple WMS layers in OpenLayers

text: I have an array of layers that I am passing to `map.addLayers()`. When I pass the array to the function in `$(document).ready(function (){})`, it doesn't add them to the map (i.e. `map.layers` does not have any of the new layers added. But if I open up the JavaScript console and run the function, it works and my added layers appear. Here's my code:               function crear( n ){             $.getJSON('json/capas.json', function(data) {                 console.log(data);               for (var i=0; i < n.length; i++) {                         capaObjetos[i] = new OpenLayers.Layer.WMS(                                       data[0][i]["id"],                                       data[0][i]["url"],                                       {                                       "layers": data[0][i]["layers"],                                         transparent: true,                                         format:"image/png",                                         CQL_FILTER: queryCQL(i),                                         projection: map.projection,                                         displayOutsideMaxExtent: false                                       },                                       {                                         singleTile: true,                                         ratio: 1,                                         buffer: 0,                                         visibility: false                                       }                                   );                       // console.log(capaObjetos[i]);                         }                      for (var key in capaObjetos) {                       var value = capaObjetos[key];                       agregadas.push(value);                 }                  });          }          agregadas = [ ];               var crasicaule, microfilo, rosetofilo, submontano, perennifolia, caducifolia, encino, pino, mezquite, tascate, mesofilo, oyamel, inducidos, halofilo,      natural, cultivados, gipsofilo, riego, temporal, mezquital, tular, sinvegetacion, asentamientos, carretera;          var capaObjetos = {               "0": crasicaule,               "1": microfilo,               "2": rosetofilo,               "3": submontano,               "4": perennifolia,               "5": caducifolia,               "6": encino,               "7": pino,               "8": mezquite,               "9": tascate,               "10": mesofilo,               "11": oyamel,               "12": inducidos,                "13": halofilo,               "14": natural,               "15": cultivados,               "16": gipsofilo,               "17": riego,                "18": temporal,               "19": mezquital,               "20": tular,               "21": sinvegetacion,               "22": asentamientos,               "23": carretera     };          $(document).ready(function () {                 capaBase();                 TreeMenu();          // ======= ADD LAYERS TO THE MAP ===========                     map.addLayers(agregadas);     });      OpenLayers is initialized in `capaBase()`. I've tried adding them individually in a `for` loop but that also doesn't work either. Any ideas? **UPDATE** If I change my create function and move my JSON locally, it works just fine.                 function crear( n ){           console.log("crear capa: " + n);                for (var i=0; i < n.length; i++) {                  capaObjetos[i] = new OpenLayers.Layer.WMS(                               capas[0][i]["id"],                               capas[0][i]["url"],                               {                               "layers": capas[0][i]["layers"],                                 transparent: true,                                 format:"image/png",                                 CQL_FILTER: queryCQL(i),                                 projection: map.projection,                                 displayOutsideMaxExtent: false                               },                               {                                 singleTile: true,                                 ratio: 1,                                 buffer: 0,                                 visibility: false                               }                           );               console.log(capaObjetos[i]);         }         // push capas de 'USO DEL SUELO Y VEGETACION' a una arregla, agregadas                   for(var key in capaObjetos) {                       var value = capaObjetos[key];                       agregadas.push(value);                   }       }      Here's the JSON I'm using.               [{               "0": {"id":"crasicaule", "class":"MATORRAL", "field":"CVE_UNION", "values": [{"0":"MC"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "1": {"id":"microfilo", "class":"MATORRAL", "field":"CVE_UNION", "values": [{"0":"MDM"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "2": {"id":"rosetofilo", "class":"MATORRAL", "field":"CVE_UNION", "values": [{"0":"MDR"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "3": {"id":"submontano", "class":"MATORRAL", "field":"CVE_UNION", "values": [{"0":"MSM"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "4": {"id":"perennifolia", "class":"SELVAS", "field":"CVE_UNION", "values": [{"0":"SAP"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "5": {"id":"caducifolia", "class":"SELVAS", "field":"CVE_UNION", "values": [{"0":"SBC"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "6": {"id":"encino", "class":"BOSQUES", "field":"CVE_UNION", "values": [{"0":"BQ"},{"1":"BQP"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "7": {"id":"pino", "class":"BOSQUES", "field":"CVE_UNION", "values": [{"0":"BP"},{"1":"BPQ"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "8": {"id":"mezquite", "class":"BOSQUES", "field":"CVE_UNION", "values": [{"0":"BP"},{"1":"BPQ"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "9": {"id":"tascate", "class":"BOSQUES", "field":"CVE_UNION", "values": [{"0":"BJ"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "10": {"id":"mesofilo", "class":"BOSQUES", "field":"CVE_UNION", "values": [{"0":"BM"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "11": {"id":"oyamel", "class":"BOSQUES", "field":"CVE_UNION", "values": [{"0":"BA"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "12": {"id":"inducidos", "class":"PASTIZALES", "field":"CVE_UNION", "values": [{"0":"PI"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},                "13": {"id":"halofilo", "class":"PASTIZALES", "field":"CVE_UNION", "values": [{"0":"PH"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "14": {"id":"natural", "class":"PASTIZALES", "field":"CVE_UNION", "values": [{"0":"PN"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "15": {"id":"cultivados", "class":"PASTIZALES", "field":"CVE_UNION", "values": [{"0":"PC"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "16": {"id":"gipsofilo", "class":"PASTIZALES", "field":"CVE_UNION", "values": [{"0":"PY"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "17": {"id":"riego", "class":"AGRICULTURA", "field":"CVE_UNION", "values": [{"0":"RA"},{"1":"RAP"},{"2":"RAS"},{"3":"RP"},{"4":"RS"},{"5":"RSP"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},                "18": {"id":"temporal", "class":"AGRICULTURA", "field":"CVE_UNION", "values": [{"0":"TA"},{"1":"TAP"},{"2":"TAS"},{"3":"TP"},{"4":"TS"},{"5":"TSP"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "19": {"id":"mezquital", "class":"VEGETACION", "field":"CVE_UNION", "values": [{"0":"MKX"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "20": {"id":"tular", "class":"VEGETACION", "field":"CVE_UNION", "values": [{"0":"VT"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "21": {"id":"sinvegetacion", "class":"VEGETACION", "field":"CVE_UNION", "values": [{"0":"DV"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "22": {"id":"asentamientos", "class":"INFRASTRUCTURA", "field":"CVE_UNION", "values": [{"0":"AH"}], "layers":"mx:usv250s5_union", "url":"http://localhost:8082/geoserver/wms"},               "23": {"id":"carretera", "class":"INFRASTRUCTURA", "field":"EDO_CLAVE", "values": [{"0":"29"}], "layers":"tlx:carreteraTLX", "url":"http://localhost:8082/geoserver/wms"}}]      For maintenanance reasons, I would like to keep the JSON as a separate file and to have it loaded into the function. So my updated question is: Why would loading the JSON from an external file be treated differently in OpenLayers, than loading the JSON locally?


title: How to convert .shp file to zip file with mapnik?

text: I am creating android application for campus indoor map. I research so much think.I decide use osmdroid. I write here what I have done 1) I draw indoor map with JOSM. 2) I convert file.osm file to file.shp with Qgis I found this code               import mapnik     m = mapnik.Map(600,300)     m.background = mapnik.Color('steelblue')     s = mapnik.Style()     r = mapnik.Rule()     polygon_symbolizer = mapnik.PolygonSymbolizer(mapnik.Color('#f2eff9'))     r.symbols.append(polygon_symbolizer)     line_symbolizer = mapnik.LineSymbolizer(mapnik.Color('rgb(50%,50%,50%)'),0.1)     r.symbols.append(line_symbolizer)     s.rules.append(r)     m.append_style('My Style',s)     ds = mapnik.Shapefile(file='ne_110m_admin_0_countries.shp')     layer = mapnik.Layer('world')     layer.datasource = ds     layer.styles.append('My Style')     m.layers.append(layer)     m.zoom_all()     mapnik.render_to_file(m,'world.png', 'png')     print "rendered image to 'world.png'"      **but it convert .png. How can I convert it zip file for showing it with osmdroid on my android app. this ?**


title: Proper workflow for LAS to DEM conversion

text: I regularly come across LiDAR (.LAS) files with elevation data. When I do, I always struggle to find a way of converting them to a DEM for hydrological analysis. I have developed several workflows, and even if they usually (but not always!) output what I need, they don't seem very optimized. Also, I would like to have one workflow that I can rely on, and not have to try three different every time. These are my current workflows: _FME_ : LAS reader --> PointCloudCombiner --> RasterDEMGenerator --> Geotiff writer. (reference time 5 min). _ArcMap_ :   1. LAS to multipoint --> Multipoint to Terrain --> Terrain to raster (10 min).    2. LAS to multipoint --> Create TIN (since Topo to raster can't read Point.Z values) --> TIN to raster (15 min).    3. Create new Mosaic dataset and add LAS as raster. Export to ESRI GRID. _LAStools_ : LAStoTXT --> Make XY event layer --> To Shapefile --> Topo to raster (3-4 hrs) I also know about the LAS dataset to raster and LAS dataset to TIN etc, but they usually won't recognize my LAS files. From what I understand #1 in Arcmap is the ESRI preferred method (?). All I want is a **hydrologically correct raster** that I can use for further analysis. What would you use? I have access to ArcGIS Standard 10.1 (soon 10.2), 3D analyst, Spatial analyst. Opensource solutions are welcome as well as scripting if necessary (Python).


title: How to Project a point on feature in PyQgis?

text: In Arcgis IProximityOperator.QueryNearestPoint Method was available to get the nearest point. Similarly Is there any function available in pyQgis? Thanks in Advance.


title: Version problem When getting GML data with Mapxtreme WFS Client

text: When producing GML data with MapXtreme WFS Client in .Net, I have introduced the problem which is parameter of "gmlversion" is working on "GML3", but not working on "GML2". Has Anyone an idea about this problem ? My code as at the below; MapInfo.Data.MultiFeatureCollection mfc = MapInfo.Wfs.Client.WfsClient.GetFeature(tbWfsServerUrl.Text, new string[] { featureTypeName }, null, **gmlVersion** , 100000, null);


title: Is there a limit for how many layers can be put into a group?

text: I am Running QGIS Lisboa 1.8.0 on a MacBook Pro. I have several vector layers with similar information which clutter the Layers panel. I created a group to put them in but when I drag and drop the layers into the new group I get just so many into the new group and then the layers panel becomes unstable. I can start to drag another layer and the new group disappears. QGIS then exits generating a large error report. I am using the method given in the manual for QGIS 1.7.0 and it appears to work up to the point where it all goes wrong. Is this a known problem with this version of QGIS? I havent been able to find a reference to this particular problem. Nigel


title: Date Parameter in model builder - ArcMap

text: Basically I'm building a model where the user has to select the data from a drop down list and the result is executed. I have a problem when trying to pass a date as a parameter. Below you can find the query builder I'm using. ![enter image description here](http://i.stack.imgur.com/WestO.png) The part which is not highlighted is working successfully. The problem is with the DATE_ parameter. When verifying the statement "an invalid SQL statement was used" error is being given. Is the problem regarding date conversion? Any idea for a possible solution for this problem is much appreciated. Thanks


title: Color blending/compositing in web mapping frameworks?

text: The color/layer blending capabilities of QGIS 2 produce marvelous results. Blended rasters are especially nice to look at. Blending is, to me at least, a better way of combining layers than mucking about with the transparency. I'm wondering how to achieve the same result in a web mapping framework. Mapbox/TileMill supports a wide range of compositing operations but these don't work with base maps. I'm looking for a mechanism that does this in the browser. OpenLayers 3 seems capable of performing such operations but I can only find examples of how to dynamically control brightness/contrast. Maybe the rest is in the making? I have a single-color raster that I want to blend/composite with the base map. I realize that blending changes the color but that's not a problem as I'm interest in perserving and communicating the saturation. PS One of the options is to blend the layer and base map in QGIS, save as image, georeference image and upload to server as raster. This, of course, is a very undesireable way of proceeding.


title: Cache generation failed. Progress gets stuck at 0%

text: I'm trying to publish a map service containing a large .tif raster (48GB). I managed to publish a small cached raster successfully (for testing), but when I tried to publish the large one, the cache generation got stuck at 0% even though it kept saying it's in progress, and again at 24%. Now I can't even publish the small one because the application is not responding right after I select "Publish". **Details:** It's an ArcGIS Server Enterprise 10.1 on a Windows Server 2008 OS machine.   * the machine is on a Microsoft Azzure cloud.    * the machine has 4 CPU cores (so I set the number of instances to 5).   * it's a raster TIF file with LZW compression.   * I chose 5 scales, and tried different formats for the tiles.  ![enter image description here](http://i.stack.imgur.com/6nJtx.jpg) Thanks for your help


title: Creating a GeoExt Action from OpenLayers Function

text: I am trying to create a GeoExt function out of this OpenLayers function, which displays the getFeatureInfo. I was able to create this GeoExt.Action() but it does not give me any output. Any idea where am I going wrong? This is the function addToPopup(). This is my application, `http://128.196.142.12/geo/test/test_new.html`.


title: google maps aerial photography

text: Ive been playing around with geoserver, google maps and openlayers and have been reasonably successful however when using the aerial photography from google i dont seem to be able to zoom in as far as i should be able to. Any ideas? Thanks in advance


title: Viewing Esri ASCII Raster file in ArcGIS for Desktop?

text: i have some problem viewing an ascii file in a structured way in arcGIS. i have used IDL ENVI too to solve this issue but all data is not displayed. I'll really appreciate your help. please check this link for files. https://drive.google.com/?authuser=0#folders/0B8iAKVklcaMUVU5qY2NNVFBRbDQ


title: Mapfish/Geoserver printing: layer problems

text: I am in the process of migrating my map application to another server with a newer Geoserver version (from 2.4 to the latest version available (master) using the 2.7-SNAPSHOT version of the printing extension (for testing reasons)). My print button in my application was working so far on the old server, but on the new server I keep getting this kind of error:               Error while generating PDF:     com.itextpdf.text.DocumentException: org.mapfish.print.InvalidJsonValueException: spec.layers[0].baseURL has an invalid value: http://a.tile.openstreetmap.org/ (spec.layers[0].baseURL has an invalid value: http://a.tile.openstreetmap.org/)     at com.itextpdf.text.pdf.PdfDocument.add(PdfDocument.java:786)     at com.itextpdf.text.Document.add(Document.java:278)     at org.mapfish.print.config.layout.Page$1.add(Page.java:76)     at org.mapfish.print.config.layout.MapBlock.render(MapBlock.java:81)     at org.mapfish.print.config.layout.Page.render(Page.java:73)     at org.mapfish.print.config.layout.MainPage.render(MainPage.java:63)     at org.mapfish.print.config.layout.Layout.render(Layout.java:59)     at org.mapfish.print.output.AbstractOutputFormat.doPrint(AbstractOutputFormat.java:30)     at org.mapfish.print.output.PdfOutputFactory.print(PdfOutputFactory.java:58)     at org.mapfish.print.MapPrinter.print(MapPrinter.java:173)     at org.mapfish.print.servlet.MapPrinterServlet.doCreatePDFFile(MapPrinterServlet.java:366)     at org.mapfish.print.servlet.MapPrinterServlet.createPDF(MapPrinterServlet.java:184)     at org.mapfish.print.servlet.MapPrinterServlet.doPost(MapPrinterServlet.java:112)     at javax.servlet.http.HttpServlet.service(HttpServlet.java:647)     at javax.servlet.http.HttpServlet.service(HttpServlet.java:728)     at org.springframework.web.servlet.mvc.ServletWrappingController.handleRequestInternal(ServletWrappingController.java:159)     at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153)     at org.springframework.web.servlet.mvc.SimpleControllerHandlerAdapter.handle(SimpleControllerHandlerAdapter.java:48)     at org.springframework.web.servlet.DispatcherServlet.doDispatch(DispatcherServlet.java:923)     at org.springframework.web.servlet.DispatcherServlet.doService(DispatcherServlet.java:852)     at org.springframework.web.servlet.FrameworkServlet.processRequest(FrameworkServlet.java:882)     at org.springframework.web.servlet.FrameworkServlet.doPost(FrameworkServlet.java:789)     at javax.servlet.http.HttpServlet.service(HttpServlet.java:647)     at javax.servlet.http.HttpServlet.service(HttpServlet.java:728)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:305)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.geoserver.filters.ThreadLocalsCleanupFilter.doFilter(ThreadLocalsCleanupFilter.java:27)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:243)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:74)     at org.geoserver.wms.animate.AnimatorFilter.doFilter(AnimatorFilter.java:70)     at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:70)     at org.geoserver.filters.SpringDelegatingFilter.doFilter(SpringDelegatingFilter.java:45)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:243)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.geoserver.platform.AdvancedDispatchFilter.doFilter(AdvancedDispatchFilter.java:49)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:243)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:311)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)     at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:116)     at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:83)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)     at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)     at org.springframework.security.web.access.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:113)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)     at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.geoserver.security.filter.GeoServerAnonymousAuthenticationFilter.doFilter(GeoServerAnonymousAuthenticationFilter.java:53)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)     at org.springframework.security.web.authentication.www.BasicAuthenticationFilter.doFilter(BasicAuthenticationFilter.java:150)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)     at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)     at org.geoserver.security.filter.GeoServerBasicAuthenticationFilter.doFilter(GeoServerBasicAuthenticationFilter.java:82)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)     at org.springframework.security.web.context.SecurityContextPersistenceFilter.doFilter(SecurityContextPersistenceFilter.java:87)     at org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1.doFilter(GeoServerSecurityContextPersistenceFilter.java:52)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)     at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.springframework.security.web.FilterChainProxy.doFilter(FilterChainProxy.java:173)     at org.geoserver.security.GeoServerSecurityFilterChainProxy.doFilter(GeoServerSecurityFilterChainProxy.java:134)     at org.springframework.web.filter.DelegatingFilterProxy.invokeDelegate(DelegatingFilterProxy.java:346)     at org.springframework.web.filter.DelegatingFilterProxy.doFilter(DelegatingFilterProxy.java:259)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:243)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:75)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:243)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:42)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:243)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.geoserver.filters.SessionDebugFilter.doFilter(SessionDebugFilter.java:47)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:243)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.geoserver.filters.FlushSafeFilter.doFilter(FlushSafeFilter.java:43)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:243)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:109)     at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:243)     at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:210)     at org.apache.catalina.core.StandardWrapperValve.invoke(StandardWrapperValve.java:222)     at org.apache.catalina.core.StandardContextValve.invoke(StandardContextValve.java:123)     at org.apache.catalina.authenticator.AuthenticatorBase.invoke(AuthenticatorBase.java:472)     at org.apache.catalina.core.StandardHostValve.invoke(StandardHostValve.java:171)     at org.apache.catalina.valves.ErrorReportValve.invoke(ErrorReportValve.java:99)     at org.apache.catalina.valves.AccessLogValve.invoke(AccessLogValve.java:953)     at org.apache.catalina.core.StandardEngineValve.invoke(StandardEngineValve.java:118)     at org.apache.catalina.connector.CoyoteAdapter.service(CoyoteAdapter.java:408)     at org.apache.coyote.ajp.AjpProcessor.process(AjpProcessor.java:200)     at org.apache.coyote.AbstractProtocol$AbstractConnectionHandler.process(AbstractProtocol.java:589)     at org.apache.tomcat.util.net.JIoEndpoint$SocketProcessor.run(JIoEndpoint.java:310)     at java.util.concurrent.ThreadPoolExecutor.runWorker(ThreadPoolExecutor.java:1145)     at java.util.concurrent.ThreadPoolExecutor$Worker.run(ThreadPoolExecutor.java:615)     at java.lang.Thread.run(Thread.java:724)     Caused by: org.mapfish.print.InvalidJsonValueException: spec.layers[0].baseURL has an invalid value: http://a.tile.openstreetmap.org/ (spec.layers[0].baseURL has an invalid value: http://a.tile.openstreetmap.org/)     at org.mapfish.print.map.readers.HTTPMapReader.checkSecurity(HTTPMapReader.java:136)     at org.mapfish.print.map.readers.HTTPMapReader.<init>(HTTPMapReader.java:76)     at org.mapfish.print.map.readers.TileableMapReader.<init>(TileableMapReader.java:42)     at org.mapfish.print.map.readers.OsmMapReader.<init>(OsmMapReader.java:50)     at org.mapfish.print.map.readers.OsmMapReader$Factory.create(OsmMapReader.java:43)     at org.mapfish.print.map.readers.MapReaderFactoryFinder.create(MapReaderFactoryFinder.java:36)     at org.mapfish.print.map.MapChunkDrawer.renderImpl(MapChunkDrawer.java:120)     at org.mapfish.print.ChunkDrawer.render(ChunkDrawer.java:55)     at org.mapfish.print.ChunkDrawer.tableLayout(ChunkDrawer.java:49)     at com.itextpdf.text.pdf.PdfPTable.writeSelectedRows(PdfPTable.java:786)     at com.itextpdf.text.pdf.ColumnText.goComposite(ColumnText.java:1765)     at com.itextpdf.text.pdf.ColumnText.go(ColumnText.java:882)     at com.itextpdf.text.pdf.ColumnText.go(ColumnText.java:877)     at com.itextpdf.text.pdf.ColumnText.go(ColumnText.java:866)     at com.itextpdf.text.pdf.PdfPRow.writeCells(PdfPRow.java:549)     at com.itextpdf.text.pdf.PdfPTable.writeSelectedRows(PdfPTable.java:764)     at com.itextpdf.text.pdf.PdfPTable.writeSelectedRows(PdfPTable.java:889)     at com.itextpdf.text.pdf.ColumnText.goComposite(ColumnText.java:1774)     at com.itextpdf.text.pdf.ColumnText.go(ColumnText.java:882)     at com.itextpdf.text.pdf.ColumnText.go(ColumnText.java:877)     at com.itextpdf.text.pdf.ColumnText.go(ColumnText.java:866)     at com.itextpdf.text.pdf.PdfDocument.addPTable(PdfDocument.java:2602)     at com.itextpdf.text.pdf.PdfDocument.add(PdfDocument.java:727)     ... 98 more     Caused by: org.mapfish.print.InvalidJsonValueException: spec.layers[0].baseURL has an invalid value: http://a.tile.openstreetmap.org/     at org.mapfish.print.map.readers.HTTPMapReader.checkSecurity(HTTPMapReader.java:133)     ... 120 more      Can somebody point me in the right direction what to look for in order to get rid of this error? Thanks Jan


title: SDO_Inside equivalent in PostGIS for 3D objects

text: I am looking for a solution for querying a point whether it is located inside a 3D object or not. If the point is located on the surface of 3D, I could use ST_3DIntersects or ST_3DIntersection. But, if the point does not touch the surface, but is located inside the exterior, how could I query it? In short, I am looking for SDO_Inside equivalent in PostGIS for 3D object. ST_Contains does not support 3D. Here is an example on using ST_3DIntersects. I chose a point that are inside the region. Result is false. SELECT ST_3DIntersects(checkpoint,region) As checkpointinregion FROM (SELECT ST_GeomFromText('POINT(2604000 1205000 1500)', 2056) As checkpoint, (SELECT ST_GeomFromText('POLYHEDRALSURFACE(((2600000 1200000 1000, 2600000 1200000 2000, 2600000 1300000 2000, 2600000 1300000 1000, 2600000 1200000 1000)), ((2600000 1200000 1000, 2600000 1300000 1000, 2700000 1300000 1000, 2700000 1200000 1000, 2600000 1200000 1000)), ((2600000 1200000 1000, 2700000 1200000 1000, 2700000 1200000 2000, 2600000 1200000 2000, 2600000 1200000 1000)), ((2700000 1300000 1000, 2700000 1300000 2000, 2700000 1200000 2000, 2700000 1200000 1000, 2700000 1300000 1000)), ((2600000 1300000 1000, 2600000 1300000 2000, 2700000 1300000 2000, 2700000 1300000 1000, 2600000 1300000 1000)), ((2600000 1200000 2000, 2700000 1200000 2000, 2700000 1300000 2000, 2600000 1300000 2000, 2600000 1200000 2000)))', 2056)) As region) as myobjects; ##  checkpointinregion f (1 row)


title: ArcObjects conform to .NET standards Error

text: So, I'm following the "Walkthrough: Migrating VB6 to VB.NET for ArcGIS 10" to convert my VB code to .NET (using Visual Studio 2010), and I'm running into this error, "NullReferenceException was unhandled by user code", object reference not set to an instance of an object. It suggests that I use a "new" keyword to create an object instance. I imported all my namespaces under Option Explicit (for this example the GeoDatabase namespace was imported), and turned Option Strict On, and applied all the CType conversions that VS required. Since this walkthrough example is for ArcGIS 10, and I'm using 9.3.1, could this be the issue? I tried to pass the "new" keyword below, but I'm still getting the error. Here is the code. Any input will be greatly appreciated.               If SelectionSet.Count < 1 Then             MsgBox("Please select one or more wells", MsgBoxStyle.Critical, "No Selected Wells")             Exit Sub             'Retrieve all features             'Set pFeatCursor = pFClass.Search(QFilt, False)         Else             'Retrieve just the selected features             'This statement moves the selected features into a feature cursor.  All             'subsequent work is done using the feature cursor.             'UPGRADE_WARNING: Couldn't resolve default property of object pFeatCursor. Click for more: 'ms-help://MS.VSCC.v90/dv_commoner/local/redirect.htm?keyword="6A50421D-15FE-4896-8A1B-2EC21E9037B2"'             SelectionSet.Search(QFilt, False, CType(FeatCursor, ICursor))         End If              'Get the first feature         Dim Feature As IFeature = Nothing          Feature = FeatCursor.NextFeature '<<<ERROR      Thank You! **Update** So, I have my ICommand Sub that calls my CreateStick Sub with the new ByRef Cursor As ICursor exception.                Private Sub ICommand_OnClick() Implements ICommand.OnClick             Me.CreateStick()         End Sub           Public Sub CreateStick(ByVal QueryFilter As IQueryFilter, ByVal Recycling As Boolean, _         ByRef Cursor As ICursor)      I stil have two errors on this line:               SelectionSet.Search(QFilt, False, ByRef CType(FeatCursor, ICursor))      ByRef and ICursor are underlined. ByRef = Expression expected and ICursor error = ICursor is a type and cannot be used as an expression.


title: Elevation Data not showing on Globe

text: Using ArcGLobe I'm attempting to display elevation data. I'm using Digital Terrain Elevation Data (DTED) data and trying to load it in the following fashion.               IRasterLayer rasterLayer = new RasterLayerClass();     rasterLayer.CreateFromFilePath(myFilePath);     Globe.AddLayerType(rasterLayer , esriGlobeLayerType.esriGlobeLayerTypeElevation , true);      Nothing shows up on the globe. Thoughts? Ideas?


title: QGIS 2.4 Processing (Zonal Statistics) error

text: In QGis 2.2 all my models (processing model builder) are working fine. When I start the same models in 2.4 I get errors like this: > Prepare algorithm 5: Zonal Statistics Running Zonal Statistics [6/8] > Parameters: RASTERLAYER_RASTER > =C:\Temp\processing\779310258f5d4cc5a0d3f0c747b2627f\RESULT.tif, > VECTORLAYER_GRENZEN > =C:/Tempordner_GIS/ASL_EK/Bereinigt/018-01/DGM5_UTM32_6stell.shp Prepare > algorithm 0: Zonal Statistics Running Zonal Statistics [1/1] Parameters: > INPUT_RASTER > =C:\Temp\processing\779310258f5d4cc5a0d3f0c747b2627f\RESULT.tif, RASTER_BAND > =1, INPUT_VECTOR > =C:/Tempordner_GIS/ASL_EK/Bereinigt/018-01/DGM5_UTM32_6stell.shp, > COLUMN_PREFIX =_, GLOBAL_EXTENT =True Failed Failed Error executing > algorithm 5 Error executing algorithm 0 'NoneType' object has no attribute > 'GetGeoTransform' See log for more details This is the error log: > Uncaught error while executing algorithm Traceback (most recent call last): > Traceback (most recent call last): File > "C:/PROGRA~1/QGISCH~1/apps/qgis/./python/plugins\processing\core\GeoAlgorithm.py", > line 217, in execute self.processAlgorithm(progress) File > "C:/PROGRA~1/QGISCH~1/apps/qgis/./python/plugins\processing\algs\qgis\ZonalStatistics.py", > line 78, in processAlgorithm geoTransform = rasterDS.GetGeoTransform() > AttributeError: 'NoneType' object has no attribute 'GetGeoTransform' Has anyone a tip for me? Thank you.


title: How to choose which band is used to calculate DN when converting a raster to vector?

text: I currently have a raster layer of flood water levels (elevation) which I need to convert to a shapefile for future use. My raster is 4 bands: RGB + one band with values from 0 to 10, representing 10 classes of water height. What I did : Raster > Conversion > Polygonise (Raster to Vector). It created a new shapefile with a column 'DN' but when I look at it I see the following values : 0, 63, 124, 185 and 255. It matches with the values contained in the first band of my raster. My guess is that when the shapefile was created, it automatically took the values of the first band. How can I force it to take the values of the 4th band instead ? Thanks ! **Edit:** If it can help someone else, here is what I finally did : Raster > Conversion > Polygonise (Raster to Vector)   Fill all the fields properly and simply edit the code block and add "-b" and the number of the band which you want to keep the values. ![Raster to vector band choosing](http://i.stack.imgur.com/gWdqg.png)


title: arcpy symbology colors changing on RefreshActiveView()

text: I have a problem that I'm not sure is functionality not built into arcpy or I'm just doing something wrong. I'm using arcpy.mapping.UpdateLayer() to apply symbology to a layer, but afterwards some elements of the symbology needs to be removed depending on the attributes of the value field. I've achived this programmatically, determining which symbology values need to be removed. I'm using something similar to the following code to remove undeeded values from the symbology.               lyr = arcpy.mapping.Layer(<layer name>)          a = lyr.symbology.classValues     b = lyr.symbology.classLabels          a --> [u'1', u'2', u'3', u'4', u'5']     b --> [u'Low', u'Moderate', u'High', u'Very High', u'Extreme']          a.remove('5')     b.remove('Extreme')          a --> [u'1', u'2', u'3', u'4']     b --> [u'Low', u'Moderate', u'High', u'Very High']          lyr.symbology.classValues = a     lyr.symbology.classLabels = b          arcpy.RefreshActiveView()      The method seems to work great. The problem is the colors set by the symbology change after the arcpy.RefreshActiveView(), which defeats the whole purpose of applying symbology in the first place. Is there a way to keep the colors the same after removing symbology values? Thanks in advance!


title: What's next if our KML is too big/complex for the Google Maps API?

text: Our web app includes simple mapping capabilities (currently just markers and KML overlays on an embedded Google Map). This works pretty well; the only real limitation we face is KML overlays that go beyond Google's size & complexity restrictions for KML. We're considering standing up our own server (e.g. GeoServer or ArcGIS Server); but that seems like a huge step just to serve up (for example) 15 MB of KML when Google's limit is 10 MB. I need a sanity check: Is there some middle ground between Google's free and easy API for KML overlays, and setting up my own tile server?


title: Determine tile identity of tile set displayed on Google Maps/OpenStreetMap

text: I'm not sure if this is best way of approaching this, so I will explain where I'm coming from first. I want to load data onto my webpage containing a map based off a lat/long point. The only way I know how to do this is to constantly be getting the a lat/lng point in the displayed map. Obviously this isn't very efficient since I either need to keep track of this point, or constantly be polling it. Each time I had the updated this point, I'd have to do some sort of processing to determine it's map location, which would add even more inefficiency. When I had that, I could then proceed to load my data. My idea was that if I could pull down an id of map tiles of a particular region at a set zoom level, then I wouldn't need to do the crazy/complicated lat/lng processing. But, my stumbling block is that I don't know of any way this is possible.


title: Saving selection as new layer in QGIS 2.4?

text: I've selected a number of polygons from one layer and want to create a new layer. Because I can't right click using a mac how do I find "Save Selection As". Its not in the list in the "Layer" drop-down from the top toolbar. I'm using QGIS 2.4 with Mac OS X 10.9.4.


title: OpenLayers 3, Markers appear in different places when viewing with iPad

text: I have an OpenLayers3 map showing static map image an ontop of that a vector layer containing markers (features). The features display nicely when viewing the map on a laptop or a pc, but when viewing the map on an iPad the markers markers appear to have a some kind of offset (so that they dont't align correctly with the map). I am using pixel projection. Does anybody have any idea why this is happening? I've tested on multiple different browser (Safari and Chrome on the iPad), and I have came to a conclusion that it is not a browser issue.


title: domain for transverse mercator projection, using proj.4

text: I have to apologize for the naive question -- it's my first encounter with projections and this software package, and I intend to use it for a non- geographical application (hence I work on a unit sphere) - when using the transverse mercator transformation, I find that it is not defined for latitude values between pi/2 and 3*pi/2, as well as -pi/2 and -3*pi/2 and thus discontinuous - why should that be, and how do I make this compatible with "my definition" of spherical coordinates, as (r,theta,phi) with -pi<=theta<=pi and -pi/2<=phi<=pi/2? It's easy to see that for phi I can just substract pi/2 (thus set the equator to 0), but what about theta? Thanks a lot, Ben               import matplotlib.pyplot as plt     import numpy             as np          from pyproj import Proj          pi=np.pi          tmerc = Proj("+proj=tmerc +a=1")               def f(a,b): return tmerc(a,b,radians=True,errcheck=True)               nice   =[]     naughty=[]          for i in     np.arange(-2*pi,2*pi,.2):         for j in np.arange(-2*pi,2*pi,.2):             try:    a,b=f(i,j) ; nice.append((i,j))             except: naughty.append((i,j))          plt.plot(np.array(nice   )[:,0],np.array(nice   )[:,1],color='green',marker='o',markeredgewidth=0.,linestyle='None')     plt.plot(np.array(naughty)[:,0],np.array(naughty)[:,1],color='red'  ,marker='o',markeredgewidth=0.,linestyle='None')


title: Filegeodatabase fanout and Python problem

text: I have fme wb which does fanout with gdb. It means that in the end as output sets of data I get few new created gdb. But If I start this fme wb within python script It does not do fanout. Just write data in one gdb. i do not understand where could be a problem?


title: Import 3D Analyst Functions With ArcPy

text: This really seems like I should be able to do this, but I am running into some problems. I want to import 3D Analyst with python, in the same way I import Spatial Analyst, but it returns an error. For Spatial Analyst I use the following:               from arcpy.sa import *      This allows me to use the spatial analyst functions in python scripts outside of ArcMap. I want to do the same with 3D analyst, but I am running into an error. I really feel like it should just be:               from arcpy.3D import *      I also tried arcpy.3d, but either way I get the error:               Parsing error <type 'exceptions.SyntaxError'>: invalid syntax (line 1)      Any help would be most appreciated, thanks in advance.