title: Display Background Value as hollow, for more than one set of RGB values

text: I have two overlapping .ecw layers of the same area, though one is clipped in the middle, with a white edge. When I set the background value to 255,255,255, and define it as hollow, it only hollows that specific RGB value, when in reality, there are other "white" values in my layer, such as: 254,254,254; 253,253,253.. all the way to 249,249,249. Is there a way to define the display background value to include more that one unique RGB set of values? I know I can work around this (clipping, raster calculator, etc.), but I want to retain only one original .ecw file without producing duplicates, therefore I only want to use Arcmap symbology or image rendering tools. **Edit** : basically, what I need is exactly like the `Custom transparency options` in QGIS, but for Arcmap.


title: mapnik generate_tiles.py error using .xml from TileMill

text: I wanted to render png tiles (like form tile server of osm). I made postGIS db, succesfully run TileMill with OSMBright style. I read that I need now export mapnik xml from TileMill. I checked Mapnik as in python tutorial- everything works fine. I've changed generate_tiles.py:   * changed the .xml file to read    * changed the bbox (I have checked bbox of small area while exporting to PNG in TileMill)   * changed the render line to `render_tiles(bbox, mapfile, tile_dir, 15, 16 , "mymap")` but an error occures: `Failed to load fonts from: ./fonts in Map at line 3 of '/home/myacc/bin/mapnik/osm.xml/'` I think the .xml is not made properly (it has written path as reference to current path) but maybe I am wrong. My questions:   1. I am doing it wrong (generate xml from tilemill, put it into ~/bin/mapnik and changing generate_tiles.py)   2. Is there a better way to make this tiles   3. If not- what I should do to make this work (where are those fonts it cant find (I think there will be more errors with reading files) ) **edit** **1** I have checked : `python -c "import mapnik;print mapnik.fontscollectionpath"` and changed the fontpath in osm.xml, but there is no such fonts registered by Mapnik . I printed out known fonts by: `python -c "from mapnik import FontEngine as e;print '\n'.join(e.instance().face_names())"` but there isn't a `font face 'Free Sans Semibold'`


title: Can a DXF File contain raster images (TIFF/JPG)?

text: I am maintaining a large Program with a graphical editor. Our Program can export DXF and it can show Background Raster Images of various types (TIFF, JPG, BMP etc.). Is it possible to embed these Raster Images in the DXF export so that they will be displayed when such a file is opened with AutoCad?


title: Enhanced vegetation index using ENVI 4.7 bandmath

text: I computed Enhanced Vegetation Index using band math in ENVI 4.7 software. The expression used was: 2.5*((float(b1)- float(b2))/(float(b1)+ (6*float(b2)) +(7.5* float(b3)))) BYTSCL(float(b1), min = -1.0, max = 1.0) The output EVI file has values varying from 0 to 255 rather than -1 to +1. Can anyone suggest an edit to the expression or a method.


title: SQL Server 2008 Distance and Performance Problem

text: I had been working on PostGIS to do some geographic processes and due to some reasons, I had to start to use SQL Server 2008. I ported same code for SQL Server but performance is almost dying when I sort according to distances (ORDER BY dist). I dont know if i am missing something on SQL server because i am new on it.               DECLARE @myPoint geography;     SET @myPoint = geography::STGeomFromText('POINT(28.8801844444445 41.0299561111109)',4326);     SELECT TOP 8 id,     @myPoint.STDistance(geom.ToString()) as dist     FROM MY_DB     ORDER BY dist      This query takes almost 4 seconds on SQL Server while it is taking 120 ms on postgis. Thanks for answers.


title: Need a Python field calculation that works like the DatePart( ) VB Function

text: I have a field called Abandoned Year (short integer data type) that I am trying to populate with a four digit year from a column called Decommission Date (Date data type). I have done this in the past using the DatePart( ) VB function but how do I accomplish this same task in a Python field calculation? Here is the code written in VB: AbandonedYear = DatePart ("YYYY", [DecommissionDate]) Thanks


title: OpenLayers and custom Proj4 Projection

text: ok so I'm a newbee in using OpenLayers. I try to map some planetary data on my webMap server using OpenLayers. I defined a proj4 projection as follow:               Proj4js.defs["MYDATUM"] = "+proj=longlat +a=1740000 +b=1740000 +units=m +no_defs";      In the html file I also call the proj4js.js, and then I added the following lines:                map = new OpenLayers.Map('map', {                         controls: [                             new OpenLayers.Control.Navigation(),                             new OpenLayers.Control.PanZoomBar(),                             new OpenLayers.Control.ScaleLine({bottomOutUnits: ''}),                             new OpenLayers.Control.MousePosition(),                             new OpenLayers.Control.Graticule()                              ],                                   Projection:  new OpenLayers.Projection("MYDATUM")                          });      My issue is that the scaleline does not show the correct distance and seems not to take into account my own projection. Actually I tried a couple of example and it seems that the line                  projection:  new OpenLayers.Projection("MYDATUM")      does not have any effect on the ScaleLine size. How can I display correctly a scale bar using my own projection/datum ? Again, I'm not working with terrestrial data, so I do not want wgs84 or any of it's friend. thank you for any ideas/feedbacks/help A.


title: Issue clipping vector file with QGIS (geoprocessing tools)

text: I am trying to clip wind speed data (.shp format) to Canada's provincial boundaries (also .shp format). I am using the clip feature in vector > geoprocessing tools. For some reason the function always stops about three quarters of the way through and never completes the clip. See screenshot: ![enter image description here](http://i.stack.imgur.com/Bz6ep.png). Has anyone experienced this before and is anyone aware of any solutions? Thanks for your help! Best, Laura


title: model to connect points with line and set the digitize direction according to height

text: I'm trying to develop a tool that connects a point layer (manholes) with lines (waste water pipes) and sets the line direction according to an height attribute acording to gravity(from high to low). The first part will be using generate near table and then xy to line but with the second part I need help - I prefer a solution with model builder, is it possible using model builder? Or do I have to use python? And how to make the tool work only on selected features? Thanks. ![enter image description here](http://i.stack.imgur.com/ypytO.png)


title: Conversion to SDE changes names

text: I could not find anything on this topic already discussed. We converted our file-based system in ESRI to an SDE instance on Sql Server. We made changes in our code to open using the connection string and we change ArcMap to use the SDE as the default database. We cannot open our GP services to republish them which is an indicator that ArcMap cannot find the schemas mentioned in the GP service ParameterInfo. When I look at the database through ConfigurationManager the names have been changed from {schemaname} to {dbname}.DBO.{schemname} and ESRI does not let me change the names back. So do we have to alter every program to use the new schema name? And are we locked into a given database name? (If our client published to a different named server I can see us losing our ability to find the schemas once again) Personally I think it is shortsighted to put the server name into the name of all of the fields in the ESRI database. Surely there is a way around this???


title: How to return feature.shape from server by using WCF and ArcGIS Silverlight API?

text: I am using the ArcGIS Silverlight API. I all I need to do is to get the geometry from a feature by objectid to zoom to it:               getFeature(objectid);      So I used WCF service to get the geometry but I received the following error:               Type 'System.__ComObject' cannot be serialized      Can anyone help?


title: How to run SDE view commands in Oracle DB?

text: I would like to create SDE view to compare two tables in my oracle spatial database. I already prepared the view commands. But I dont know where and how to run this? I unable to run this commands in Toad SQL editor like other sql statements. Please help me if any one have idea.


title: Arcgis 10.x Create equidistant polyline that are perpendicular to an axis

text: I have a polyline that represents a cyclone track and I want to create equidistant polyline that are perpendicular to the cyclone track. Many thanks for your help. Elie


title: Not seeing the actual data value in getObservation response from MapServer 6

text: I'm experimenting with MapServer and SOS but I'm not getting back any useful data from the GetObservation call I'm doing. My .map file uses a Postgres query to get data and when I run that same query in PGAdmin I get the four columns I expect: id, geom, date and total_alkalinity_value. All the columns have the right data. Below is the getObservation response in which I see only the gml:boundedBy and ms:msGeometry values, no sign of the id or total_alkalinity_value.               <om:ObservationCollection xmlns:gml="http://www.opengis.net/gml" xmlns:ows="http://www.opengis.net/ows/1.1" xmlns:swe="http://www.opengis.net/swe/1.0.1" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:sos="http://www.opengis.net/sos/1.0" xmlns:ms="http://mapserver.gis.umn.edu/mapserver" xmlns:om="http://www.opengis.net/om/1.0" gml:id="WQIS_TOTAL_ALKALINITY" xsi:schemaLocation="http://www.opengis.net/om/1.0 http://schemas.opengis.net/om/1.0.0/om.xsd http://mapserver.gis.umn.edu/mapserver http://wqis-sos.test.niwa.co.nz/cgi-bin/wqis-sos?service=WFS&version=1.1.0&request=DescribeFeatureType&typename=total_alkalinity">     <om:member>     <om:Observation>             <om:samplingTime>         <gml:TimeInstant>           <gml:timePosition>1989-04-12</gml:timePosition>         </gml:TimeInstant>       </om:samplingTime>       <om:procedure xlink:href="urn:ogc:def:procedure:80"/>       <om:observedProperty>         <swe:CompositePhenomenon gml:id="total_alkalinity_value" dimension="3">           <swe:component xlink:href="urn:ogc:def:property:OGC-SWE:1:id"/>           <swe:component xlink:href="urn:ogc:def:property:OGC-SWE:1:date"/>           <swe:component xlink:href="urn:ogc:def:property:OGC-SWE:1:total_alkalinity_value"/>         </swe:CompositePhenomenon>       </om:observedProperty>       <om:featureOfInterest xlink:href="http://wqis-sos.test.niwa.co.nz/cgi-bin/wqis-sos?service=WFS&version=1.1.0&request=DescribeFeatureType&typename=total_alkalinity"/>       <om:result>         <gml:featureMember>           <ms:total_alkalinity>             <gml:boundedBy>               <gml:Envelope srsName="urn:ogc:crs:epsg:4326" srsDimension="2">                 <gml:lowerCorner>167.294382 -46.243591</gml:lowerCorner>                 <gml:upperCorner>167.294382 -46.243591</gml:upperCorner>               </gml:Envelope>             </gml:boundedBy>             <ms:msGeometry>               <gml:Point srsName="urn:ogc:crs:epsg:4326" srsDimension="2">                 <gml:pos>167.294382 -46.243591</gml:pos>               </gml:Point>             </ms:msGeometry>           </ms:total_alkalinity>         </gml:featureMember>       </om:result>     </om:Observation>     </om:member>     </om:ObservationCollection>      And here is the relevant stuff from the layer part of the .map file LAYER .... METADATA .... "sos_timeitem" "date" "sos_procedure" "ALKALINITY_SENSOR" "sos_procedure_item" "total_alkalinity_value" "sos_offering_id" "WQIS_TOTAL_ALKALINITY" "sos_observedproperty_id" "Total Alkalinity" "sos_describesensor_url" "http://wqis- sos.test.niwa.co.nz/ALKALINITY_SENSOR.xml"               TYPE POINT     STATUS ON     NAME "total_alkalinity"     DUMP TRUE     DEBUG 0          CONNECTIONTYPE POSTGIS     CONNECTION XXXXX     DATA "geom from (         SELECT se.id AS id,              l.geom AS geom,              CASE WHEN COALESCE ((se.start_date).month, 0) = 0                   THEN                     (se.start_date).year::CHAR(4)                   WHEN COALESCE ((se.start_date).day, 0) = 0                   THEN                     (se.start_date).year || '-' || (se.start_date).month || '-01'                   ELSE                     (se.start_date).year || '-' || (se.start_date).month || '-' || (se.start_date).day              END::DATE AS date,              av.value::NUMERIC AS total_alkalinity_value         FROM XXXXXXXX     ) AS se USING UNIQUE id USING SRID=4326"      So my question is: How do I get id and total_alkalinity_value in my getObservation response? *updated to show more of the .map file and add sos_dateitem* *Now that I have the sos_procedure_item set to total_alkalinity_value the om:procedure element contains the data value (in this response it is 80), but that seems like totally the wrong place for it*


title: How to export 2D vector graphic from ArcScene to Illustrator?

text: I need to know how to export 2d vector graphic from ArcScene so I can edit it with Adobe Illustrator. I've tried changing the options when I export 2D scene as AI file or PDF file, but they always remain as a unique image that does not allow me to edit the lines and polygons. Anyone can help?


title: How to select a subwindow from a source image in EPSG 3857?

text: I need to select a subwindow from an image projected to epsg 3857. If I use, for example:                gdal_translate -projwin ${UL_x_3857} ${UL_y_3857} ${LR_x_3857} ${LR_y_3857} source_3857.tif source_cut3857.tif      where the corners are given in the 3857 projection. I got                Input file size is 1691, 1699      Computed -srcwin 54 40 1258 -2147483648 from projected window.      -srcwin 54 40 1258 -2147483648 falls outside raster size of 1691x1699      It works fine with UTM coordinates.. what am I doing wrong?


title: get segments of constant length from Linestring MS SQL SERVER

text: For Example I have simple linestring 'LINESTRING(25 50, 100 125, 150 190)' and I need to divide it to segments of length 10. In postgis i know there is a function ST_Segmentize but what about MS SQL Server?


title: Null output from field calculator code block

text: I have a simple calculation to compare two columns and populate a new column with the lower value but I am getting null output...any ideas? expression               re_score(!Signal_Strength!,!Signal_Strength_1!,out1)      code block               out1='99'     def re_score(f1,f2,out1):      if f1<f2:       out1 = f1      else:       out1=f2      return out1      > WARNING 000405: No records within table Succeeded at Tue Mar 06 15:13:45 > 2012 (Elapsed Time: 7.00 seconds)


title: Conflating Two Network Datasets in ArcGIS

text: I am looking to merge two network datasets in ArcGIS 10.0. I have an ArcInfo license and have trawled this forum but am looking for some clarification. I want to merge two road networks but only make changes to one of the datasets. The two tools I have been investigating are Integrate and Spatial Adjustment. From this post: Advise with 'Integrate' tool in ArcToolbox It seems that Integrate actually changes both datasets. Is it possible to weight one dataset so high that it does not change? Secondly, there doesn't seem to be a way to bulk create links to spatial adjust a very very large network. Does anyone have any recommendations? Thank you all in advance, Kel


title: AHP Workflow in GIS

text: We are to apply AHP method to determine waste disposal areas using GIS, the AHP methodology supposes, from what I gathered, alternatives known before hand (as explained in this very straightforward document). This will be the first time we work with AHP, thus, I am here to take away some confusion and see how this methods combines with GIS ( we work under ArcGIS 10.1 which sadly do not handle AHP, so we'll have to go with the usual geoprocessing tools and do the AHP calculations in some other software). Now, here's a flowchart I've prepared, again, this is only my first understanding of it : ![enter image description here](http://i.stack.imgur.com/Bai6n.jpg) What I am really most confused about is the extraction of the alternatives (which are here spatial areas), is the raster resulting from rasters weighted sum sufficient to indicate preliminary alternatives ? Could anyone forward me to a more "viable" workflow for using AHP in GIS ? Thanks.


title: Split linestring to smaler lines

text: Hi I need help splitteing a linestring into smaller linestrings, using `postgis`. I need to route to closest point on the linestring, so i need to have more points on the line               Sample Linestring = LINESTRING(11.959813031048 55.3269904906477,11.9613429171704 55.3257458375313,11.9656473425314 55.325719907258,11.9675402524794 55.3251753715195)      I have no luck with the following               select st_astext(st_segmentize(st_linefromtext('LINESTRING(11.959813031048 55.3269904906477,11.9613429171704 55.3257458375313,11.9656473425314 55.325719907258,11.9675402524794 55.3251753715195)'     ),2) ) as geom       But something is missing? I would like to have my line split into lines max 10 meters between points. I hope somebody could get me back on track. * * * I then dump the line into points               select st_astext((st_dumpPoints(st_segmentize(st_linefromtext('LINESTRING(11.959813031048 55.3269904906477,11.9613429171704 55.3257458375313,11.9656473425314 55.325719907258,11.9675402524794 55.3251753715195)'     ),10 ))).geom) as geom


title: Geoserver Tile Generating

text: I'm a complete newbie at this so please excuse my lack of knowledge. I need to be able to place hundreds of thousands of points on a map. I'm using OpenLayers with Google Maps base layer and I just installed OpenGeo Suite. I also need to be able to add/delete points and manipulate them real-time. Placing those points as a vector or marker layer doesn't work so I think I need to have the server generate tiles? Where do I start with this? I've scoured the documentation and this website and can't find where to begin with generating tiles and doing it dynamically with my data. Thanks!


title: How to check whether a fc is registered as versioned or not using python

text: How can I check whether or not a feature class is versioned using python? I'm trying to register feature classes using a python script. Before I do this I need to check whether a feature class is versioned. Are there keywords that we can use in coding while looping through the feature classes? Thanks in advance for the help.


title: Is there a Python Lib for requesting WMS/WFS and saving as image/PDF?

text: i am wondering whether there's Python open-source GIS lib which has APIs to support call WMS/WFS from another GIS server (e.g., GeoServer) and then save the response data(WMS Basemap and WFS layer) as pictures. any recommendations? thanks for any inputs! **UPDATE** : what I am trying to do is a Map Printing service, by using OpenLayers as the front-end and Django as the server; Client user set the extent and layers and then send the print request (which refers to the parameters, i.e., map extent, names of layers) to server, then server takes over this request and call WMS/WFS again by using request parameters, save the response as PDF, export this PDF link to client. The difficult part is that how the server call WMS/WFS and combine/overlay these responses together (i.e., put these map/layers together, since WMS is usually the base map, WFS points to the feature layers), finally save this combined object as Image. in current answers, urllib seems a good one, but i am not sure how to combine these responses (WMS, WFSs) together; OWSLib also seems another good option, but it indicates it's a client programming tool, I am a little confused that whether it's appropriate for my use... any other further inputs??? appreciate!


title: How is the default label field chosen by ArcMap?

text: I have 2 questions: I know how to manually go to Layer Properties>Labels tab>Label Field and choose my label field. I also know that I could set up my primary label field and save the setting as a layer (*.lyr) file. This question, though, is more about how the label field is chosen by the program prior to user interaction. First: **When a new layer (shapefile, ArcSDE layer, personal/file geodatabase featureclass) is added to ArcMap, how is the default label field chosen?** From my experience, it isn't chosen by field order, nor alphabetical by field. My best guess so far is that it _might_ chose the first field of text/varchar (alphanumeric) type. Second: **Is there a way to define what the "primary label field" will be for a layer within the shapefile, sde, or gdb files** , or are you stuck with saving layer files or letting the system pick for you.


title: Adding/updating field to polygon feature class that lists bordering (neighbor) polygons?

text: I have a dataset of zipcodes. I want to add a new field "NEIGHBORS" that includes bordering polygons in a comma delimited manner such as attached image: ![bordering polygons added to field](http://i.stack.imgur.com/BX1fB.jpg)


title: Setting the Color of a TextSymbol

text: public void DrawText( IActiveView view, string text, string fontName, int fontColor, float fontSize, IPoint point ) {                  ITextSymbol textSymbol = new TextSymbol();                  stdole.IFontDisp font = textSymbol.Font;             font.Name = fontName;             font.Size = (decimal)fontSize;             textSymbol.Font = font;                  IColor color = textSymbol.Color;             color.RGB = fontColor;             textSymbol.Color = color;                  ITextElement textElement = new TextElementClass();             textElement.Text = text;             textElement.Symbol = textSymbol;             textElement.ScaleText = true;                  IElement element = textElement as IElement;             element.Geometry = point;                  view.GraphicsContainer.AddElement( element, 0 );                  view.PartialRefresh( esriViewDrawPhase.esriViewGraphics, element, null );         }      Can anyone suggest to me why this code always sets the font and size correctly but sometimes fails to set the color correctly? I always see the correct color in the 'textSymbol.Color' property after the line that sets it to 'color.' On the map (within a layer) it sometimes appears to 'combine' with that layer's fill color and off the map it sometimes appears 'blue.'


title: How to attribute a polygon with the names of polylines it contains?

text: I have a set of polyline shapefiles and a polygon shapefile. The task is to split the polylines by the polygon border and write the names of polyline objects located within the polygon to the attribute table of the polygon shapefile. It's clear to me how to split lines by a polygon; **the problem is how to get the names and add them to the attributes of the polygon shapefile**. Is there any way to do it in ArcGIS (ArcView license)? Thank you.


title: How to convert 6 digit Latitude & 7 digit longitude to degrees

text: Sorry, but I'm a bit knew to the whole GIS world, and though I am making great headway drawing Polygons and such in Bing Maps, I am currently stuck. In our SQL Spacial Data, I came across MULTIPOLYGON data, where all the lat/long values are very large numbers (6 digits for Latitude and 7 digits for Longitude, before the decimal) I found a sql method to dismantle the MultiPolygon into several Polygons, but those large values do not work in Bing Maps or Google Maps. So I've been searching for an algorithm to convert those large numbers down to, what I have perceived as, regular Latitude/Longitude numbers, such as 45.19167074460551 -93.583288812804 (which does work on online maps). With that said, could somebody please either show me or point me in a great direction to learning how to convert 631456.83356859628, 3816200.7961099017 to something I can use online. I know the coordinates fall somewhere in Arizona. Your answers are MUCH appreciated. Jonathan


title: ArcGIS Geocode Single Address Python

text: How can I geocode a single address in python using arcgis 9.3? I found http://webhelp.esri.com/arcgiSDEsktop/9.3/index.cfm?TopicName=Geocode_Addresses_%28Geocoding%29 - but this does a whole table, and I was forced to use a MDB file as the feature class and then inport and export from it.


title: ArcObjects edit/update ODBC Oracle table

text: Can an ODBC Oracle table attached by OLE within ArcMap be edited/updated by using ArcObjects? If so, what are the procedural steps for doing this? Are there code examples out there showing how to do this? Thank You


title: Are there torrents that contain large free spatial datasets?

text: I've been having a lot of issues downloading entire datasets of SRTM or ASTER GDEM2 or similar free datasets, mainly because of poor services. Does anyone know of any torrents which contain such free datasets?


title: OpenLayers PanZoom control customization

text: I try to customize PanZoom control and have some questions. 1) What are OpenLayers.Control.PanZoom.X and OpenLayers.Control.PanZoom.Y constants? 2) Method "draw" takes "px" argument as input. I can not find where this method is called and what value of "px" it is passed? 3) Is it possible to change buttons div size through css? I've discovered that size and filenames of images are hardcoded.


title: How to resolve pgrouting dijkstra_sp() function error?

text: insert into dijsktra_results(geom) select ST_UNION(the_geom) from dijkstra_sp('roads_layer',                 (select source_id from roads_layer where gid=9975)::integer,                 (select target_id from roads_layer where gid=2967)::integer);      This query to perform the routing operations. "dijsktra_results" is a newely created table for routing operations with gid and geom columns of 2D(MultiLineString). This is the error I am getting               ERROR:  column "the_geom" does not exist     LINE 1: SELECT gid,the_geom FROM shortest_path('SELECT gid as id, so...                        ^     QUERY:  SELECT gid,the_geom FROM shortest_path('SELECT gid as id, source::integer, target::integer, length::double precision as cost FROM roads_layer', '1' , '20142' , false, false), roads_layer where edge_id = gid      CONTEXT:  PL/pgSQL function "dijkstra_sp" line 13 at FOR over EXECUTE statement          ********** Error **********          ERROR: column "the_geom" does not exist     SQL state: 42703     Context: PL/pgSQL function "dijkstra_sp" line 13 at FOR over EXECUTE statement      Can anyone resolve this problem or guide for an extent am I gone wrong anywhere.Till now everything worked fine.I calculated vertex_id,edge_id,length.... And how this, dijkstra_sp() function works and calculates the shortest path.


title: Show Fields in infowindow from my custom layer

text:
I want to display 2 fields from the POIData layer. What i get now is ${ToegangsNr} ${Baanvak} in text. Does anyone know how to do this in an DynamicMapServiceLayer. Thanks in advance               POIData = new esri.layers.ArcGISDynamicMapServiceLayer("https://gis.strukton.com/ArcGIS/rest/services/POI_Data/MapServer", { id: "POI", "visible":false });             map.addLayer(POIData);                                                        dojo.connect(map, "onClick", function(evt) {                                     map.infoWindow.resize(175, 125);                                         var symbol = new esri.symbol.SimpleMarkerSymbol().setStyle(esri.symbol.SimpleMarkerSymbol.STYLE_X).setSize(12);                                         symbol.outline.setWidth(4).setColor("blue");                                              dojo.connect(map, "onClick", function(evt) {                                                 map.graphics.clear();                                                 var graphic = new esri.Graphic(evt.mapPoint, symbol);                                                 map.graphics.add(graphic);                                                      map.infoWindow.setTitle("${ToegangsNr}");                                                 map.infoWindow.setContent("
Baanvak: ${Baanvak}
	");                                                           map.infoWindow.show(evt.screenPoint, map.getInfoWindowAnchor(evt.screenPoint));                                         });                                              dojo.connect(map.infoWindow, "onHide", function() {                                                 map.graphics.clear();                                         });                                                                         });


title: "Add Comment" functionality in Flex Viewer

text: I want to add some kind of commenting functionality for users to leave notes on a map. I know that Esri deprecated the "add pop up" element in 2.3.1 because it didn't save. My current thinking is to direct users to create a circle feature and add their comment as an attribute using the "Edit" widget, however that seems needlessly complicated. Does anyone have a simpler solution for saving location-aware comments on the map?


title: R Raster package loop over pixels

text: I want to perform a trend analysis for each pixel over it's time dimension. Therefore I wrote a loop over each pixel doing the calculation but this is extremely slow. I am sure there must be a faster way to perform the computation. Also the dimensions of the file are 1440,720,480. Is this a doable size for R? The problem is that i cannot find the statistical function implemented in IDL which I otherwise use.               gridfile = 'test.nc'     data = brick(gridfile)          require(fume)          mk = matrix(numeric(0), 720,1440)           for(i in 1:72) {       for(j in 1:144) {       temp = mkTrend(as.vector(data[i,j,]))       mk[i,j] = temp$`Sen's Slope`       }     }


title: Simplifying TopoJSON causes leak of Natural Earth in d3

text: I'm working with topoJSON and D3 to display Natural Earth data on a globe. Obviously, we need to simplify our geometries while still keeping borders consistent. We've created a script as based on Mbostock's Map Tutorial. Here's the link to our gist. Once we run the script, it outputs all the countries as .json files. ![json files](http://i.imgur.com/6QcT2yT.png) However, when we load the USA, we encounter a problem. When we mouse over the USA, the entire globe is highlighted. ![United States without mouse hover](http://i.imgur.com/dxbgUZm.png) United States without mouse hover. ![United States on mouse hover](http://i.imgur.com/O2y3CZf.png) United states on mouse hover Is there something in our script that would cause this? We need to simplify our geometries for performance purposes, and we need to preserve boarders. TopoJSON was the clear choice, but we could be calling it incorrectly.


title: How can I copy columns from layer 1 to layer 2?

text: I have two layers, the first layer is pipe: that has some fields (column) title like : length, id, material. Also I have another layer that is empty (without any fields). How can I copy fields (column) title from layer 1 to layer 2? I am new to QGIS.


title: Dynamically have symbology change as a due date nears?

text: I need to make a countdown based on symbolization. Let's say I have several groups of sections with a Date field in it. I have when the section is due in the date field. Can I get it to dynamically hatch close lines as the date nears?


title: Where can I get the location of all Universities and Colleges in North America?

text: This seems fairly straightforward but I'm having trouble finding the right data. All I need is the name, address, and lat long of all universities and colleges in North America. I've come across a list of many universities in America here: http://answers.google.com/answers/threadview/id/578065.html but just by name.


title: Selecting MAX values within records with the same name

text: I'm working in ArcMap 9.3 and I've never dealt with that kind of queries... I need to Select the MAX value from column [A] for each set of records in the table which have the same value in column [B] For example               [A]  [B]          100  aaa             101  aaa          102  aaa          999  aaa          100  bbb          101  bbb          102  bbb          999  bbb so the selection would be record#4 (where [A]=9999 which is the biggest value of all 'aaa') and record#8 (where [A] = 999 which is the biggest of all 'bbb') Is it possible to do that using queries and subqueries when working with databases? Any help appreciated! Dan


title: How to create a movable point in OpenLayers?

text: I have the following code:                   function mapCreate(lon,lat){         map = new OpenLayers.Map("map1");         var osm = new OpenLayers.Layer.OSM();         vectors = new OpenLayers.Layer.Vector("Vector Layer");         map.addLayer(osm);         var center = new OpenLayers.LonLat(lon,lat).transform(             new OpenLayers.Projection("EPSG:4326"),             map.getProjectionObject()         );              point = new OpenLayers.Geometry.Point(center.lon,center.lat);         vectors.addFeatures([new OpenLayers.Feature.Vector(point)]);         drag = new OpenLayers.Control.DragFeature(vectors);         //map.addLayer(vectors);         map.addControl(drag);         drag.activate();         map.setCenter(center, 15);         map.addLayer(vectors);         point.events.register('moveend',point, function(evt){             alert('hello');         });          }      The part that does not work is:                       point.events.register('moveend',point, function(evt){             alert('hello');         });      Can anyone help me with creating an event listener for a point in open layers?


title: Change zIndex of Vectors in openlayers

text: I have a bunch of linestring features on a vector layer. Some of them overlap. If I trigger the hover event on those features, I want them be in front. I did find some functions to set the z-index of layers, but nothing for manipulation individual features. I tried using the individual style maps, but that doesn't work. What am I doing wrong?                   // Add vectory layer         this.vector = new OpenLayers.Layer.Vector("Track", {             renderers: ['Canvas', 'SVG', 'VML'],             rendererOptions: { zIndexing: true },             styleMap: myStyle         });         tbMap.map.addLayer(this.vector);      In a loop various lines are added like so:                       //make the new tracksegment             myTrack = new OpenLayers.Feature.Vector(new OpenLayers.Geometry.LineString(path), {hash: hash});                  // Draw vector on the layer             this.vector.addFeatures([myTrack]);      And the stylemap:               var myStyle = new OpenLayers.StyleMap({         "default": new OpenLayers.Style({             pointRadius: config.pointRadius,             strokeColor: config.strokeColor,             fillColor: config.fillColor,             fillOpacity: config.fillOpacity,             strokeWeigth: config.strokeWeigth,             strokeWidth: config.trackWeight,             strokeOpacity: config.strokeOpacity,             fontColor: config.fontColor,             fontSize: config.fontSize,             fontFamily: config.fontFamily,             graphicZIndex: 0         }),         "select": new OpenLayers.Style({             pointRadius: config.pointRadius,             strokeColor: config.strokeColor_selected,             fillColor: config.fillColor,             fillOpacity: config.fillOpacity,             strokeWeigth: config.strokeWeigth,             strokeOpacity: config.strokeOpacity,             fontColor: config.fontColor,             fontSize: config.fontSize,             fontFamily: config.fontFamily,             graphicZIndex: 100         }),         "hover": new OpenLayers.Style({             pointRadius: config.pointRadius,             strokeColor: config.strokeColor_hover,             fillColor: config.fillColor,             fillOpacity: config.fillOpacity,             strokeWeigth: config.strokeWeigth,             strokeWidth: config.trackWeight,             strokeOpacity: config.strokeOpacity,             fontColor: config.fontColor,             fontSize: config.fontSize,             fontFamily: config.fontFamily,             graphicZIndex: 100         })      })


title: Pywps Timeout long processes

text: I'm executing a long process WPS of around 180sec, and when it start after 60sec the WPS (response) is abandoned, and the process continue to run until the end (checked on the log file). I tried the solution proposed by Pywps wiki of reducing the timeout of the Apache server (http://wiki.rsg.pml.ac.uk/pywps/Async_issue) without success ! I found the response which seems to be the only issue : http://lists.wald.intevation.org/pipermail/pywps-devel/2013-April/001598.html But I didn't find how can I insert my "call process" by self.cmd("process.sh...).


title: Where can I get a comprehensive list of Rules expression operators?

text: I'm selecting and styling my points layer on a contour map of Philadelphia using Rules Based Styling. (see: Can I use rules based rendering in QGIS to select features?). The filtering process using [string_field = 'string'] turns out not to be accurate for one of my three rules. I have 'Homocide - Justifiable' and it returns points with (at least) 'Thefts'. I'd like to try another expression. In any case it would be nice to have a list of all the possible expressions.


title: Indian Latitude Longitude values

text: I am in need of around 10000 and odd Indian latitude and longitude values for my openlayers project. I have to load them in a postgis db. so a sql file or atleast text file would be better than a link.


title: query / filter on attribute tables with external tables joined

text: i tried to do some queries on joiend table fields ... i've got/found some answers concerning version 1.8 - which result in "use select + export to create a combined dataset" so i just did some tests with Layer / query in qgis 2.2 valmiera and 2.3-master(build from 2014-05-02)) and ended up with the same problems: while   * "layer properties / style / column expression" tool is working fine based on "external" fields - even if i build complex expressions   * and "table - select faetures by expressionn" doesn't seem to have any problems with where-clauses "layer / query" returnes strange results with my demo-data; i.e.:   * no problems with expressions unsing fields from the layer's attribute table    but - as soon as i address "external" fields   * but - as soon as i address "external" fields     * the test feature returns an error: a simple query like "zsp_experimente_xls_WFL01_AREA" >= 1.52" (where 1.52 is valid value out of the sample-list) returns    "OGR 3 error 1: type mismatch or improper type of arguments to >= operator"      but "zsp_experimente_xls_WOHNFL01_AREA" LIKE "%" returns all records when testing the query ("the where clause teturned 1364 rows")     * the confusing thing about that:    the the sample tool privides correct results ! ! !   ((which means to me that the select part of the statement works fine)) so finally i am still argueing that one of the essential benefits of the join- concept is, that you do not have to mix up geometry and attribute data -> which was a really great idea !!!!!!!      to select data by expression and export geometry and "added" additional data really **undermines** this **concept** by **duplicating data** !!! so: does anybody have any idea how to deal with that problem ? best regards robert


title: Is there a parellel line function for slds?

text: I am attempting to make an sld style which symbolizes a line with an offset. I have tried with filters, functions online graphics but can achieve only ugly results. The line consists out of 2 line symbolizers. one thick underlying line which is displayed at an offset westerly of the achse. The other thinner line lies above the thicker line overlapping the easterly edge of the thicker line. an example is here: http://maps.zgb.de/rrop/sld_img/style_hatchedline_offset.png If anyone knows what I can do to produce this sld I´d be most grateful, yours, Rob


title: Where can I find polygon data of wine (viticulture) areas and their main grape or varietal?

text: Where is coordinate data on the geographic boundaries of wine (vitaculture) areas and their main type of wine (grape or varietal)? a GIS shapefile or geoJSON file format would be perfect. Ideally I would like the **polygon areas** for wines of USA or Europe along with a few attributes corresponding to each area. But even just data for at least one major country or varietal would be nice. I can handle the geographic coordinates in any format. Thanks!


title: Network Analyst not working

text: I have network analyst installed as an extension in my arcGIS 9.3, but it just does not work!!! when I go to ArcCatalog, the option of new network data set also does not appear... :( I have to do this job of finding routes & I need this to work. I am zero at programming & I don't have arc server as well. Thanks ![enter image description here](http://i.stack.imgur.com/oyVWJ.jpg) ![enter image description here](http://i.stack.imgur.com/ykLfr.jpg)


title: Why do I get different output from gdal2tiles vs. maptiler?

text: I am trying to generate map tiles using gdal2tiles for an area much larger than MapTiler can handle. However, my tile caches are not showing up on the maps. I checked gdal2tiles by running the same file through it and MapTiler - and the outputs come out differently. Here's my gdal2tiles command: > /usr/local/Cellar/gdal/1.9.0/bin/gdal2tiles.py -p mercator --s_srs 3857 -z > 1-9 -v myFile.tiff For instance, in the maptiler output, I get the file > 9/85/197.png While gdal gives me different files, like > 9/85/313.png The output should be identical, right? Is there something I'm missing? My projection on the original file is definitely correct. I have also tried gdal2tiles with 900913. I am running gdal2tiles on MacOS X Mountain Lion, gdal 1.9.0.


title: Determine the furthest driving distance from a road layer to the nearest point of a point layer using OSM/QGIS

text: In a specific geographic area, I have a road layer and roughly 30 points dispersed around the geographic area as part of a point layer. Both layers are in SHP and PostGIS formats. The records from the point layer correspond to "stops" on the road layer. **How can I calculate the "furthest driving distance at any point in the network to the nearest 'stop' distance is ..."?** and as a follow up, can I save that route, or show the furthest points as a different layer? I'll be using QGIS and OSM or my own road layer data for this.


title: Are the public domain vehicle tracks available?

text: I'm looking for recorded vehicle tracks that can be replayed to (load) test my application. Is there such data available? Best would be a european data basis but everything else would do the job as well. Or are there (free) tools available to generate tracks that can be used to simulate a vehicle driving on a road and sending a position every 10 seconds?


title: How to replace values in rasters in directory 1 based on condition on other rasters in another directory?

text: Having 10 raster files in one directory and another 10 rasters in another directory and both have the same dimensions. directory 1:                   files1 <- list.files("C:\\2010", "*.envi", full.names = TRUE)         s <- stack(files1)      directory 2:                files2 <- list.files("C:\\2012", "*.envi", full.names = TRUE)      s1 <- stack(files2)      I would like to filter s based on s1. I want to replace values in all layers in s by NA whenever corresponding layers in s1 have values greater than 0.6(corresponding pixels).


title: How Could I Calculate Shortest Path ?

text: I am working on QGIS - 1.8 Version and frequently getting an error while using Road Graph and Shortest Path Opetions. Kindly Suggest or provide solution.. With Regards Prashant Persai


title: When editing feature-linked annotation with font color white, the annotation disappears

text: When moving or editing the annotation, the annotation disappears. We are editing on a black background, so the text isn't just blending in to the background. The attribute table does not show any abnormality, but the "Zoom to Selected" option is grayed out. It can be selected in the table, but does not display in the display window. Anyone have this problem? Found any solutions?


title: Why are the tile coordinates different in Google maps and Google earth?

text: I have a .NET application that typically uses the standard Google Maps JavaScript API V2. Within this application I am using a GTileLayerOverlay to add my own custom dynamic tile layer.                 function addJourneyLineOverlay(assetId, startDate, finishDate) {                  var path = "http://<%= Request.Url.Host %>";             if ("<%= Request.Url.Port %>" != 80) {                 path += ":<%= Request.Url.Port %>";             }             path += "/Private/polyLineTileLayer/Tile.aspx";                  var myCopyright = new GCopyrightCollection("(c) ");             myCopyright.addCopyright(new GCopyright('Demo', new GLatLngBounds(new GLatLng(-20, 130), new GLatLng(-30, 155)), 0, '@2012 IMT Pty Ltd'));                  var journeyPathLayer = new GTileLayer(myCopyright);             journeyPathLayer.getTileUrl = function (tile, zoom) {                 return path + '?zoom=' + zoom +                   '&x=' + tile.x +                       '&y=' + tile.y +                           '&assetId=' + assetId +                               '&startDate=' + startDate +                                   '&finishDate=' + finishDate;             };                  journeyPathLayer.isPng = function () {                 return true;             };             journeyPathLayer.getOpacity = function () {                 return 1.0;             };                  journeyPathOverlay = new GTileLayerOverlay(journeyPathLayer);             map.addOverlay(journeyPathOverlay);         }      This overlay requests it tiles from an aspx page that gets some points from a database, determines if the points are in the tile and then returns an image with the points on it. I developed it based on the same pattern as Gheat.net. And I am also using the Gmap.Net library for its Mercator projection, for when I need to work out which points should be in which tiles. This works really well when I am using the standard Google maps API. However, recently I have to integrate my application with our clients Google Earth Enterprise server. This means that I am now getting my map imagery and API from our clients server and not from Google proper. This has caused a problem with my overlay. My overlay is not longer position where it should be relative to the map location. They are appearing a bit lower that where they should be. Using firebug I examined the request for my tile layer and found out that the tile X and Y where different than what I was expecting. At zoom level 12, using Google maps a particular LatLng location is in tile X-3392, Y-2312. But in Google Earth the same location at the same zoom level is in tile X-3392, Y-2305. that's a 7 tile difference in the Y axis. What can account for this difference? Does Google Earth use a difference projection, if so is there a convenient .net library for it?


title: draw shp files or kml files on ipad?

text: I need to be able to draw a polygon in the field while the landowner is there to approve it. Is there an app that is easy to use so all of the field hands can operate it? I also need to be able to export the file to be used in the equipment used to spray.


title: Raster Combine - ESRI Spatial Analyst

text: I am trying to combine a whole bunch of rasters with Esri Spatial Analyst. Seems to be doing what I want EXCEPT the calculated output is only the area where all 4 rasters overlap. I am a complete noob at this, don't know what I'm doing wrong!! Using 9.3.1


title: Large point to Polygon by buffer-join-buffer-dissolve issues

text: I have a few large datasets (about 10 million points) which I am buffering, then joining, dowgrading and finally buffering and dissolving. This is all in a model which works fine with sample smaller datasets. Details are in Buffer point dataset based on compared value Is there a way to split the point dataset and then process it through to the dissolved buffer and then re-join the files? The key question here is to split a "square" dataset into blocks so that the re-join does not create jagged edges. OR do I just have to swith to another software as in ArcGIS Buffer/Dissolve performance? The model builder python code is below. BUT FIRST THE ERROR IS > Executing (Buffer): Buffer > P:\2012\Job_044_DM_Radio_Propogation\Working\FinalPropogation\TRC_Longlands\TRC_Longlands.gdb\TRC_Longlands_pts > P:\2012\Job_044_DM_Radio_Propogation\Working\FinalPropogation\TRC_Longlands\TRC_Longlands.gdb\TRC_Longlands_buf > "30 Meters" FULL ROUND NONE # Start Time: Thu Mar 08 12:54:21 2012 Succeeded > at Thu Mar 08 14:42:55 2012 (Elapsed Time: 1 hours 48 minutes 34 seconds) > Executing (Spatial Join): SpatialJoin > P:\2012\Job_044_DM_Radio_Propogation\Working\FinalPropogation\TRC_Longlands\TRC_Longlands.gdb\TRC_Longlands_pts > P:\2012\Job_044_DM_Radio_Propogation\Working\FinalPropogation\TRC_Longlands\TRC_Longlands.gdb\TRC_Longlands_buf > P:\2012\Job_044_DM_Radio_Propogation\Working\FinalPropogation\TRC_Longlands\TRC_Longlands.gdb\TRC_Longlands_join > JOIN_ONE_TO_ONE KEEP_ALL # WITHIN "30 Meters" # Start Time: Thu Mar 08 > 14:42:56 2012 ERROR 000426: Out Of Memory Failed to execute (Spatial Join). > Failed at Thu Mar 08 15:01:31 2012 (Elapsed Time: 18 minutes 35 seconds)               # ---------------------------------------------------------------------------     # xy_downgradedSS.py     # Created on: 2012-03-08 11:05:07.00000     #   (generated by ArcGIS/ModelBuilder)     # Usage: xy_downgradedSS <Imported_XY> <Buffered_tower_range> <XY_Feature_Class> <INPUT_txt_file>      # Description:      # ---------------------------------------------------------------------------          # Import arcpy module     import arcpy          # Load required toolboxes     arcpy.ImportToolbox("Model Functions")     arcpy.ImportToolbox("P:/2011/Job_237_BushfireHazard_Datasets/Working/post_erratum/ProcessingModels.tbx")          # Script arguments     Imported_XY = arcpy.GetParameterAsText(0)     if Imported_XY == '#' or not Imported_XY:         Imported_XY = "%Path%\\%Name%.gdb\\%Name%_layer" # provide a default value if unspecified          Buffered_tower_range = arcpy.GetParameterAsText(1)     if Buffered_tower_range == '#' or not Buffered_tower_range:         Buffered_tower_range = "%Path%\\%Name%.gdb\\%Name%_buf" # provide a default value if unspecified          XY_Feature_Class = arcpy.GetParameterAsText(2)     if XY_Feature_Class == '#' or not XY_Feature_Class:         XY_Feature_Class = "%Path%\\%Name%.gdb\\%Name%_pts" # provide a default value if unspecified          INPUT_txt_file = arcpy.GetParameterAsText(3)     if INPUT_txt_file == '#' or not INPUT_txt_file:         INPUT_txt_file = "P:\\2012\\Job_044_DM_Radio_Propogation\\Working\\FinalPropogation\\TRC_Longlands\\trc_longlands.txt" # provide a default value if unspecified          # Local variables:     v_Name__join = Buffered_tower_range     v_Signal_Strength__2_ = v_Name__join     Signal_Strength_Added__2_ = v_Signal_Strength__2_     Buffered_tower_range__3_ = Signal_Strength_Added__2_     Source = XY_Feature_Class     Source_Added = Source     v_Signal_Strength = Source_Added     Signal_Strength_Added = v_Signal_Strength     Path = INPUT_txt_file     Name = Path     v_Name__gdb = Name     v_Path_ = "%Path%"          # Process: Parse Path     arcpy.ParsePath_mb(INPUT_txt_file, "PATH")          # Process: ParseName     arcpy.gp.toolbox = "P:/2011/Job_237_BushfireHazard_Datasets/Working/post_erratum/ProcessingModels.tbx";     # Warning: the toolbox P:/2011/Job_237_BushfireHazard_Datasets/Working/post_erratum/ProcessingModels.tbx DOES NOT have an alias.      # Please assign this toolbox an alias to avoid tool name collisions     # And replace arcpy.gp.ParseName(...) with arcpy.ParseName_ALIAS(...)     arcpy.gp.ParseName(Path)          # Process: Create File GDB     arcpy.CreateFileGDB_management(v_Path_, "%Name%", "CURRENT")          # Process: Make XY Event Layer     arcpy.MakeXYEventLayer_management(INPUT_txt_file, "longitude", "latitude", Imported_XY, "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]];-400 -400 1000000000;-100000 10000;-100000 10000;8.98315284119521E-09;0.001;0.001;IsHighPrecision", "")          # Process: Copy Features     arcpy.CopyFeatures_management(Imported_XY, XY_Feature_Class, "", "0", "0", "0")          # Process: Add Field (4)     arcpy.AddField_management(XY_Feature_Class, "Source_Ch", "TEXT", "3", "2", "50", "", "NULLABLE", "NON_REQUIRED", "")          # Process: Calculate Field (4)     arcpy.CalculateField_management(Source, "Source_Ch", Name, "PYTHON", "")          # Process: Add Field (2)     arcpy.AddField_management(Source_Added, "Signal_Strength", "FLOAT", "3", "2", "10", "", "NULLABLE", "NON_REQUIRED", "")          # Process: Calculate Field (3)     arcpy.CalculateField_management(v_Signal_Strength, "Signal_Strength", "re_score(!Rx!,f2)", "PYTHON_9.3", "f2='99'\\ndef re_score(my_field,f2):\\n\\n if (my_field*1) < 14:\\n  f2 = '0'\\n elif (my_field*1) < 24:\\n  f2 = '1'\\n elif (my_field*1) >= 24:\\n  f2 = \"2\"\\n \\n return f2")          # Process: Buffer     arcpy.Buffer_analysis(Signal_Strength_Added, Buffered_tower_range, "30 Meters", "FULL", "ROUND", "NONE", "")          # Process: Spatial Join     arcpy.SpatialJoin_analysis(v_Signal_Strength, Buffered_tower_range, v_Name__join, "JOIN_ONE_TO_ONE", "KEEP_ALL", "", "WITHIN", "30 Meters", "")          # Process: Add Field (3)     arcpy.AddField_management(v_Name__join, "Downgraded_SS", "FLOAT", "3", "2", "10", "", "NULLABLE", "NON_REQUIRED", "")          # Process: Calculate Field (2)     arcpy.CalculateField_management(v_Name__join, "Downgraded_SS", "re_score(!Signal_Strength!, !Signal_Strength_1!)", "PYTHON", "def re_score(f1,f2):\\n    if f1 < f2:\\n        out1 = f1\\n    else:\\n        out1 = f2\\n    return out1")          # Process: Buffer (3)     arcpy.Buffer_analysis(Signal_Strength_Added__2_, Buffered_tower_range__3_, "30 Meters", "FULL", "ROUND", "LIST", "Downgraded_SS")


title: Need information about the Envisat satellite

text: I am new to remote sensing and I want to know about sensor design of Envisat satellite. Can anybody help me or preferably direct me to some good references?


title: Calculate field of view / angle of a point to a shapefile

text: I am trying to automate a calculation of the blocked fraction of a field of view from a point to a shapefile. I.e. find α in ![](http://i.stack.imgur.com/Q7ccV.png) , such that α / 360 gives me the fraction of the field of view covered by the shapefile in reference to the total field of view. I know ArcGIS features a 'field of view' based on a DEM/raster (right?), but I am looking for a simpler application, between a simple point and polygon shape. Is anyone aware of software / a function that is able to do this? It is very easy to do by hand, but I need to process roughly 250 cases, so I thought I'd better ask for other suggestions :)


title: GDB maintenance of muliple editors without ArcSDE

text: I'm working with ArcGIS 10.1 SP1. I have finished 18 months of collecting data and digitization and built an impressive GDB with a lot of tables, relationship classes and 100 feature classes. Now i need to deliver it to the costumer. The GDB need to be updated twice a year by me and the customer can make some changes at tables but not at feature classes. What are the best practices for GDB maintenance ? I thought about editor tracking and using complex relationship classes. More ideas are most welcome. Thanks a lot, Geog


title: How to convert CAD Annotation to Shapefile with attributes?

text: How to add AutoCAD Mtext to GIS attribute table using FME? What transformer I should use? The reader file in CAD and I want to transform it to GIS Shapefile. In this case I want to expose the attribute from text (fme_text_string) to DC_LAND_PARCEL_BOUNDARY(polygon geometry). So how to expose the attribute along with the value, so that in ArcGIS i only turn on the labeling options to show the annotation.


title: QGIS 2.0.1. missing Google Mercator CRS

text: I've downloaded the latest QGIS version and I cannot add basemap. I installed OpenLayers plugin and while adding any map from list QGIS crashes. I tried to enable Google Mercator CRS first, yet it's not in the list of available CRSs. Does anybody know how to fix this and how can start adding basemaps? Or is it possible to add National Geographic map as a basemap?


title: Should the import of local spatial vector files in different formats be handled Server side?

text: I'm working on a Esri Silverlight Web API app and wondering the best way to handle the import of local spatial vector data files. Sometimes these will be your standard shape file but often these will be in vendor specific proprietary file formats that will have to use special API's to parse out the data within them. My question is whether it would likely be a better practice to handle the data conversion on the server side vs. attempting to do it on the client side?


title: GLONASS Satellite Position Prediction using Almanac Data

text: GLONASS ICD has a section named "Algorithm of calculation of satellite motion parameters using almanac". It gives the formulas in order to calculate the position and velocity of the satellite. Unfortunately I cannot get the same results as in the given example, and I don't know why. Can you suggest any other documents, sample code or calculations about this problem? There is a paper by Wang et al. called "Predicting GLONASS Satellite Orbit Based on an Almanac Conversion Algorithm for Controlled Ionosphere Scintillation Experiment Planning". In this paper the proposed algorithm converts GLONASS almanac to GPS almanac, and the procedure that is used for GPS position calculation can be applied. However, there are undefined parameters in the given algorithm, namely `Dt^2`, `dt2`, `Ohm` and `Ohm_e` (Figure 4). Does anyone have an idea about this? Thanks in advance.


title: Creating individual points from count field in point feature class using ArcGIS for Desktop?

text: I have a point feature class with a count field associated with each of the points. I would like to create a new feature class where each of the points are duplicated by the value of the count field, so that rather than being an aggregate count, they represent individual events (I need to run a Monte Carlo on them so individual events are required). Does anyone have any idea how to do this in ArcGIS for Desktop?


title: Looking for US National Grid to Lon/Lat Converter in Python, C, or C++

text: I see here that there is Fortran source code for going from USNG to lon/lat: http://www.ngs.noaa.gov/PC_PROD/USNG/ Is there some python or C/C++ code that does the same? My understanding that this is _not_ a function provided by Proj4/pyproj/GDAL, correct?


title: Error Reporting Registry Location of ArcGIS for Desktop?

text: I need to stop ArcMap from error reporting during an Arcpy loop. The error report dialog stops my loop and I need for it to go away, permanently I found this link which references the registry location as 'HKCU\Software\ESRI\Settings\ErrorReports'. The actual folder hierarchy in my regedit is much different and I cannot locate the ErrorReports registry ![enter image description here](http://i.stack.imgur.com/M7VnJ.jpg) What is the actual location of the registry responsible for error reporting in ArcMap 10.2.2?


title: How to access data from a table that has been related?

text: I have a data driven page table that I am trying to associate with route line data. The object of this exercise is to print the associated data on the data driven page AS or LIKE a data driven page attribute. I have a district table which is a line route that has start and end stationing of the line through the districts. It is a very simple table with only 4 districts. So each event has a from and to station. The map center table is point based. There are a number of maps within each district, so it is a many to one relationship. Essentially, each district is a range of stations and the map centers station locations fall within those, not necessarily associated with any particular start or end point. I have tried spatially joining the route to the map centers, but did not get anything useful. It basically just added the map centers fields to the route table and only listed the route records, but there was no map center data. When I tried to spatially join the map centers to the Route table, it just gave me several locations from the last district, but there was nothing in the other districts. Note: I am using a UTM coordinate system, so low coordinate values on the Eastings are in the west, however the stationing begins in the east and gets progressively higher towards the west. The Last district of course is in the west. I get nothing when I try to relate the tables either way, other than a comment or something in the table drop down menu that identifies the name that I gave the relate. It doesn't seem to do anything useful at all. I am using Arcmap 10.2 desktop, with a standard license.


title: What is GIS for Indigenous Peoples?

text: I am working on a short course for a group of indigenous 10-15 year old children in a couple of weeks and want to start the program with a short intro (90s) of what GIS is and how it can help them. I can't find any videos like this online. I have drafted a script and want to send it out to the GIS community for your feedback before I go in a create the video. It will be animated like http://www.youtube.com/watch?v=6EAMqKUimr8 using videoscribe. Please see https://docs.google.com/document/d/10HsuZgDx04LlVPazqUfI8SMyzKpajJ2S6X5Xv79TSqA/edit?usp=sharing and comment in this site. **EDIT/UPDATE** I am stuck on 2 points, namely writing 2 more examples of how GIS can help in understanding the different elements of a community (in Landuse planning) and a few more relevant global examples of how spatial thinking helps to solve issues. See 7 and 8 in the doc. Thanks in advance for your time...the final product will be creative commons and free to use


title: Suggestions for how to secure a custom ArcGIS extension using a license file or install key

text: I'm looking for a simple methodology for licensing custom ArcGIS extensions with an installation key or license file. When the user attempts to enable the extension I want to check that they have a valid license before enabling the tools. I’m not sure if this is the best forum for this question but I bet some of my fellow GIS developers have tackled this before. Using an installation key would be the preferred approach because it seems less hassle, both for me and for the user. I’m sure there are many of ways to accomplish this, but I’m trying to work out the basic process and I have a number of questions about it.   1. I’m assuming that once the user enters the key during the installation a value gets written to the registry. The registry value is then checked to enable the extension. Is there a better way?    2. What are the options for generating the installation key? What rules can you create to determine a valid installation key vs. an invalid installation key?   3. What value should you write to the registry? Perhaps encrypt the installation key and write the encrypted value to the registry. Then when checking that key to enable the extension, decrypt the registry key, and then make sure it meets the same requirement of a valid installation key? I’m not looking for Fort Knox here, but it seems you want the registry key value to be something not completely apparent to someone snooping around. I've written extensions for ESRI software for years, (back to ArcView 3), without concern for validating licensing. Usually the extension goes to a handful of clients with an associated license agreement text file, and that was it. And, honestly, I prefer it that way. If someone wants to break a licensing agreement they will find a way, so why waste time and energy trying to stop it. However, for a variety of reasons I now need to apply a basic licensing scheme on a few extensions. So, any suggestions or guidance would be greatly appreciated.


title: Problem compiling UI file generated by QGIS Plugin Buidler

text: I have used the plugin builder QGIS 1.8 to create an empty plugin. The next step is to compile the ui file using pyuic. pyuic.py runs but halts at the line with `from PYQt4 import QTCore` and an error message of 'No module named PyQt4' is displayed. I am using an OSGeo Install and I also have an install of last weeks master. Using the plugin builder in the master produces the same issue. It could be a path issue and if anyone replies, setting the path is part of the answer, knowing where to set the path is the other part. Could someone also briefly explain what python is looking for when is is looking for a module. I am guessing that it is looking for a dll but it could be looking for an API or even a py / pyc file. Unless I can first check that the module exists with Windows explorer and can find its path, I will not later be able to set a path. An alternative would be use gnu make or other compiler to compile the plugin using the make file produced by the plugin builder. Is there a compiler that makes this simple (a pre-compiled windows installer that comes with instructions on how the compiler with a make file)? Can QT Designer compile the ui files that it edits? Andrew M


title: Map extent changes - stack overflow exception - ArcGIS Silverlight 3.1 API

text: In our Silverlight application in many work flows we are changing map extents. The below method is called multiple times, each time the application throws a stack over flow exception. We can not find the reason for this issue. Can anyone provide me with the best way to handle this issue? Please help this is an urgent production issue.               public void SetMapExtent(Envelope envelope)     {     if (envelope.SpatialReference == null)     envelope.SpatialReference = mapSection.SpatialReference;     mapSection.Extent = envelope;     }      Please help us on this issue.


title: Visualizing Spatial Query

text: I'm new to GIS applications and trying to figure out what options I have for visualizing spatial query results that combine raster and vector data using OpenLayers/GeoExt/GeoServer/PostGIS. A hypothetical query that could be created on a web client by a user might be: _"give me all raster pixels with an index value of 1 that are within a distance of 120 meters of any river (line vector feature) and touching or inside a selected watershed (polygon vector feature)"_. I think I have a good idea how I could write a where clause for a spatial SQL query such this, but it's not at all clear where I'm submitting this type of query from my Javascript web client. Some PHP examples I've found seem to suggest that a server side process will have to submit the SQL query, presumably directly to the PostGIS database? Is this correct? It's also not super clear how I can visualize mixed raster/vector query results. I've seen some purely vector examples that create a Vector map layer in open layers and when the vector features come back to the client, those features are simply pushed onto the Vector map layer. Assuming a mixed vector+raster query is possible in PostGIS, should I be expecting to follow the same process (get vector features back, then add them to a Vector OpenLayers layer) for visualization? Lastly, does it even sound like I'm thinking along the right lines here or am I attempting to do something that is difficult or even impossible given the feature sets available in Open Source GIS tools these days? Thanks in advance for any advice or direction


title: Assigning a nearest polygon to a point

text: I have a point layer of households and a polygon layer of tax lots. I want to assign a nearest tax lot polygon to each household point. Note that as the household points are coded on road, not within a tax lot, ST_Within does not work properly. Do you have any suggestions?


title: ArcSDE/ST_Geometry: Can I store polygon and point geometries in same row?

text: In the system I'm building, we've got a facility that can be represented as either a point, or a point within a polygon (the point in the polygon may or may not be the polygon's centroid). I'm using Oracle/ArcSDE with the ST_Geometry storage type. Is it possible to store both point and polygon data for a single feature? Or does each require its own record?


title: Best way to train GIS professional working in Telecom

text: This question is not about any code or bugs. This is about GIS profession. I am GIS professional working with **Mapinfo, Openlayer, SQL SERVER , Oracle and System integration .** I learned everything from stackoverflow.com and gis.stackexchange and google.com and thanks a lot to you guys. Now I want to have some training in this field related to telcom and GIS. What sort of training should I take that can be applicable to development of telecom and help use vital data to be presented as useful information? As for any technical jobs the requirement changes and need for update is essential. I want to be expert in the field of GIS related to Telecom with current simple knowledge and without training I don't think it's possible.


title: Buffer aggregation method in OpenCarto library

text: What is the use case of this function? Can it be used to aggregate building polygons as well?


title: Installing pgRouting in Windows 7

text: Can you please help me out in Installing pgRouting in Windows 7 64-bit version. I have been trying it for days. Since I am completely new to this please provide me the initial requirements and steps in installing pgRouting.


title: Web Mapping Services from Europe besides OSM

text: Wondering if there are any other Web Mapping Services based in Europe or even Asia? Fed up with the Bings and Googles and even OSMs of the WMS world. This is 2012 not 2002, shouldn't we have a much better FREE view of the world by now?


title: Clipping soil survey map service to feature class

text: I'm guessing this is not possible but is there any way to clip a Esri map service such as the Soil Survey Map for a another feature polygon? http://server.arcgisonline.com/ArcGIS/rest/services/Specialty/Soil_Survey_Map/MapServer I know it's a not a feature service but what is confusing is how if you zoom in far enough it will show the soil series information. What I would like to do is clip it to a smaller area like less than 1000 acres based on a polygon feature class. Maybe someone knows of another feature service available that would show the soil survey down to the soil series detail like this Esri map service does?


title: How can i topologically sort 2D triangles and put in an array?

text: I have a set of 2D triangles in a Geometry array which together forms a connected strip when placed them according to edge adjacency. However, in the array they are not in topological order based on the edge adjacency. How can i topologically sort them and put in the array?


title: Openlayers - how to display base layer before overlays and display overlays as they are being loaded?

text: This is my first question here. I already found a lot of good information here about openlayers so I guess it is the right place to ask. I am not a programer, I just usually adapt what I can find to suit my needs for my personnal website. I am using openlayers to display about 20 gpx files over a basemap. Each gpx file is about 600ko,each is displayed in its own overlay. As it was, everytime I added one more gpx overlay, it took longer before anything is displayed on my webpage. So I guess Openlayers is loading everything before displaying anything. So I am looking for a way to display the base layer first and then load and display one gpx overlay at a time so that the user does not have to wait for too long before he sees anything happening on his screen. Here is my webpage with the map underwork at the bottom of the page : http://www.untourenvelo.ch/carnet-de-voyage/suisse/apres- And here is what I tried to do but does not work. I put a "loadend" event on the last created overlay so that the next overlay is not created before the previous is not completely added to the map. But this does not work : The first one shows correctly (of course) but the second one does not show because the (i-1) is "undefined" according to console log. The third one is displayed and ... that's it :( If anyone can help. Here is my initOSM function :                   function initOSM() {          map= new OpenLayers.Map("CarteOSM");     var mapnik = new OpenLayers.Layer.OSM();     var fromProjection = new OpenLayers.Projection("EPSG:4326"); // Transform from WGS 1984     var toProjection = new OpenLayers.Projection("EPSG:900913"); // to Spherical Mercator Projection     var position = new OpenLayers.LonLat(13.41,52.52).transform( fromProjection, toProjection);     var zoom = 2;      map.addLayer(mapnik);     //map.zoomToMaxExtent();          map.setCenter([25,45],zoom);               var cyclemap = new OpenLayers.Layer.OSM("cyclemap","http://tile.opencyclemap.org/cycle/${z}/${x}/${y}.png",     {'sphericalMercator': true, isBaseLayer:true}     );     map.addLayer(cyclemap);           StyleGPX = new OpenLayers.StyleMap({     "default": new OpenLayers.Style({     //Icone en fonction du type de waypoint, le nom du fichier png doit être identique au texte dans la balise sym du fichier gpx     externalGraphic: "../../images/GPS/${sym}.png",     pointRadius: 2,     //graphicHeight: 60,     graphicYOffset: -32,     graphicXOffset: -5,     label: "${name}",      labelAlign: 'rb',     fontSize: 11,     fontFamily: "Arial",     fontColor: "black",     labelYOffset: 0,     labelXOffset: -10,     fillColor: "ffcc66",     strokeColor: "#ffe400",     strokeWidth: 1.5,      strokeOpacity: 1     }),     "select": new OpenLayers.Style({     externalGraphic: "../../images/GPS/${sym}.png",     graphicYOffset: -32,     graphicXOffset: -5,     pointRadius: 20,     fillColor: "#66ccff",     strokeColor: "#3399ff",     graphicZIndex: 2,     label: "${name}",     labelAlign: 'rb',     fontSize: 11,     fontFamily: "Arial",     fontColor: "black",     labelYOffset: 0,     labelXOffset: -10     })     });               StyleGPX.styles["default"].defaultStyle.externalGraphic = "";     StyleGPX.styles["default"].defaultStyle.pointRadius = 2;     StyleGPX.styles["default"].defaultStyle.label = "";               //Tableau des fichiers gpx a afficher     FichiersGPX = ["Europe","Turquie","Iran","Asie Centrale","Bolivie","Chine", "Chili_Argentine", "Indonesie"];          //Create layers from array of files' names     CalquesGPX[0] = new OpenLayers.Layer.Vector(FichiersGPX[0], {     protocol: new OpenLayers.Protocol.HTTP({                 url: "../../images/GPS/"+FichiersGPX[0]+".gpx",                 format: new OpenLayers.Format.GPX({extractWaypoints: true, extractRoutes: false, extractAttributes: true})     }),     strategies: [new OpenLayers.Strategy.Fixed()],     visibility: true,     styleMap: StyleGPX,     projection: new OpenLayers.Projection("EPSG:4326"),     eventListeners: {"featureselected": onFeatureSelect, "featureunselected": onFeatureUnselect},     'displayInLayerSwitcher':false      });     map.addLayers([CalquesGPX[0]]);               for (var i=1; i<=FichiersGPX.length-1; ++i) {     CalquesGPX[i-1].events.register('loadend',CalquesGPX[i-1],function() {     CalquesGPX[i] = new OpenLayers.Layer.Vector(FichiersGPX[i], {     protocol: new OpenLayers.Protocol.HTTP({                 url: "../../images/GPS/"+FichiersGPX[i]+".gpx",                 format: new OpenLayers.Format.GPX({extractWaypoints: true, extractRoutes: false, extractAttributes: true})     }),     strategies: [new OpenLayers.Strategy.Fixed()],     visibility: true,     styleMap: StyleGPX,     projection: new OpenLayers.Projection("EPSG:4326"),     eventListeners: {"featureselected": onFeatureSelect, "featureunselected": onFeatureUnselect},     'displayInLayerSwitcher':false      });     map.addLayers([CalquesGPX[i]]);          });     }               //Selection des elements pour affichage de pop up dans un overlay On declare le tableau complet des calques     selectionElement = new OpenLayers.Control.SelectFeature(CalquesGPX,{     clickout: true, toggle: false,     multiple: false, hover: false,     toggleKey: "ctrlKey", // ctrl key removes from selection     multipleKey: "shiftKey" // shift key adds to selection     });     map.addControl(selectionElement);     selectionElement.activate();          //Affichage de l'echelle de la carte     var echelle = new OpenLayers.Control.ScaleLine({ topInUnits: 'km', bottomInUnits: 'm'});     map.addControl(echelle);           //Affichage du selecteur de calque     var layer_switcher = new OpenLayers.Control.LayerSwitcher();     map.addControl(layer_switcher);     layer_switcher.maximizeControl();           }


title: Is there any method to analyze data from fire protection events?

text: I want to analyze data from events in fire protection a create some illustration of maps. In short: my task is to find methods for data visualization. Data are represented by point layer (each event is localized and except the coordinates, each event has some information in attribute table). Can you tell me some methods I can use (except classic chloropleth map and diagram map)? I analyzed data with kernel density tool, or I have used "count points in polygon" tool from Hawths Analysis tools for creating cartogram maps. I used this tool for Quadrat Analysis too - creating vector grid (polygons) over the extent and counting the events in each polygon (quadrat). Last method was line diagram map (creating buffer from line data and counting points in each polygon segment of buffer, this information was necessary to visualize events that happened on the roads. Is there some other method or tool I can use? Thanks for help


title: Display information in infowindow using KMLLayer and kmz file

text: I added the kmz layer to my map www.spc.noaa.gov/products/md/ActiveMD.kmz However when I click on it, all I get in the infowindow is "MD 1558". (This may change for you since it is updated in near realtime) If I download the kmz file and rename to a kml and open in Notepad I see the following code below. Is it possible to include the http://www.spc.noaa.gov/products/md/MD1558.kmz link in the infowindow? If not, is it possible to configure an infoWindow with the KMLLayer so that I could add my own link? Thank you.               <?xml version="1.0" encoding="UTF-8"?>     <kml xmlns="http://www.opengis.net/kml/2.2">     <Document>     <Folder>     <name>SPC Active MDs as of Fri Aug 15 18:11:02 UTC 2014</name>     <open>1</open>     <description>     <![CDATA[     Visit http://www.spc.noaa.gov/products/md/ for more info.     ]]>     </description>     <NetworkLink>     <open>1</open>     <name>MD 1558 Active Till 1915 UTC</name>     <visibility>1</visibility>     <refreshVisibility>1</refreshVisibility>     <description>     <![CDATA[     MD 1558 Active Till 1915 UTC.<br />     Visit http://www.spc.noaa.gov/products/md/md1558.html for more info.     ]]>     </description>     <Link>     <href>http://www.spc.noaa.gov/products/md/MD1558.kmz</href>     <refreshMode>onInterval</refreshMode>     <refreshInterval>60</refreshInterval>     </Link>     </NetworkLink>     </Folder>     </Document>     </kml>


title: Is Landsat 4-5-7 L1G geometric corrected?

text: Im using Landsat L1G for our project and need L1G to be geometric corrected for later use. One thing I know for sure that L1G need to be employed a Digital Elevation Model (DEM) for topographic accuracy. I follow a guide to use ArcMap (from ArcGis) to correct the uncorrected Landsat image. But the process required DEM and DOQ (Digital orthophoto quadrangle) image to run. DEM is available but DOQ is not avaiable outside US. As I read, ArcMap using DOQ to establish the GCPs for uncorrected image. Since the DOQ is not avaible, I wondering if only using DEM to correct the L1G image is good enough? In other hands, can anyone know that the Landsat L1G is corrected with GCPs? **Updated** : thanks all for your help. we decided to use data from L1T source instead. But I face another problem cant find any information about sattelite observation zenith angle/azimuth of Landsat image. We need it to make sure two different image has the same condition for further calculation. Any one have idea where I can find that information?


title: convert from HDF to GeoTiff

text: I am having some problems with the conversion from HDF to GeoTiff. I am following the steps described in this webpage. The problem is that I can't find what is the NDVI subdataset name (XXXX) in order to run the gdaltranslate command.               gdal_translate -of GTiff HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf"      :MOD_Grid_monthly_CMG_VI:XXXX modis_ndvi01.tif      This is the gdalinfo output for one of the files(MOD13C2.A2001001.005.2007078152825.hdf):               E:\GDAL>gdalinfo MOD13C2.A2001001.005.2007078152825.hdf     Driver: HDF4/Hierarchical Data Format Release 4     Files: MOD13C2.A2001001.005.2007078152825.hdf     Size is 512, 512     Coordinate System is `'     Metadata:       HDFEOSVersion=HDFEOS_V2.9       LOCALGRANULEID=MOD13C2.A2001001.005.2007078152825.hdf       PRODUCTIONDATETIME=2007-03-19T19:28:25.000Z       DAYNIGHTFLAG=Both       REPROCESSINGACTUAL=reprocessed       LOCALVERSIONID=5.2.1       REPROCESSINGPLANNED=further update is anticipated       SCIENCEQUALITYFLAG=Not Investigated       AUTOMATICQUALITYFLAGEXPLANATION=No automatic quality assessment is performed in the PGE       AUTOMATICQUALITYFLAG=Passed       SCIENCEQUALITYFLAGEXPLANATION=See http://landweb.nascom.nasa.gov/cgi-bin/QA_WWW/qaFlagPage.cgi?sat=terra for the product Science Quality status.       QAPERCENTMISSINGDATA=0       QAPERCENTOUTOFBOUNDSDATA=0       QAPERCENTCLOUDCOVER=0       QAPERCENTINTERPOLATEDDATA=100       PARAMETERNAME=CMG 0.05 Deg Monthly NDVI          (...)          Subdatasets:       SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly NDVI       SUBDATASET_1_DESC=[3600x7200] CMG 0.05 Deg Monthly NDVI MOD_Grid_monthly_CMG_VI (16-bit integer)       SUBDATASET_2_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly EVI       SUBDATASET_2_DESC=[3600x7200] CMG 0.05 Deg Monthly EVI MOD_Grid_monthly_CMG_VI (16-bit integer)       SUBDATASET_3_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly VI Quality       SUBDATASET_3_DESC=[3600x7200] CMG 0.05 Deg Monthly VI Quality MOD_Grid_monthly_CMG_VI (16-bit unsigned integer)       SUBDATASET_4_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly red reflectance       SUBDATASET_4_DESC=[3600x7200] CMG 0.05 Deg Monthly red reflectance MOD_Grid_monthly_CMG_VI (16-bit integer)       SUBDATASET_5_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly NIR reflectance       SUBDATASET_5_DESC=[3600x7200] CMG 0.05 Deg Monthly NIR reflectance MOD_Grid_monthly_CMG_VI (16-bit integer)       SUBDATASET_6_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly blue reflectance       SUBDATASET_6_DESC=[3600x7200] CMG 0.05 Deg Monthly blue reflectance MOD_Grid_monthly_CMG_VI (16-bit integer)       SUBDATASET_7_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly MIR reflectance       SUBDATASET_7_DESC=[3600x7200] CMG 0.05 Deg Monthly MIR reflectance MOD_Grid_monthly_CMG_VI (16-bit integer)       SUBDATASET_8_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly Avg sun zen angle       SUBDATASET_8_DESC=[3600x7200] CMG 0.05 Deg Monthly Avg sun zen angle MOD_Grid_monthly_CMG_VI (16-bit integer)       SUBDATASET_9_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly NDVI std dev       SUBDATASET_9_DESC=[3600x7200] CMG 0.05 Deg Monthly NDVI std dev MOD_Grid_monthly_CMG_VI (16-bit integer)       SUBDATASET_10_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly EVI std dev       SUBDATASET_10_DESC=[3600x7200] CMG 0.05 Deg Monthly EVI std dev MOD_Grid_monthly_CMG_VI (16-bit integer)       SUBDATASET_11_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly #1km pix used       SUBDATASET_11_DESC=[3600x7200] CMG 0.05 Deg Monthly #1km pix used MOD_Grid_monthly_CMG_VI (8-bit unsigned integer)       SUBDATASET_12_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly #1km pix +-30deg VZ       SUBDATASET_12_DESC=[3600x7200] CMG 0.05 Deg Monthly #1km pix +-30deg VZ MOD_Grid_monthly_CMG_VI (8-bit unsigned integer)       SUBDATASET_13_NAME=HDF4_EOS:EOS_GRID:"MOD13C2.A2001001.005.2007078152825.hdf":MOD_Grid_monthly_CMG_VI:CMG 0.05 Deg Monthly pixel reliability       SUBDATASET_13_DESC=[3600x7200] CMG 0.05 Deg Monthly pixel reliability MOD_Grid_monthly_CMG_VI (8-bit integer)


title: Popups over markers not working in OpenLayers

text: I have a map developed based on google map API and OpenLayers. It has some markers on it. I'm trying to create a popup for each marker. The markers are in a markers layers. Besides the markers layer, there are multiple other layers in the map as well. I'm trying to use OpenLayers.Control.SelectFeature API to do it. But the popups are not working. Following is my code. Any hints are appreciated!               function createPMap(){     if(pMap != null) pMap.destroy();          // Create pMap     pMap = new OpenLayers.Map({       projection: new OpenLayers.Projection("EPSG:900913"),       displayProjection: new OpenLayers.Projection("EPSG:4326"),       units: "m",       maxExtent: new OpenLayers.Bounds(-20037508, -20037508, 20037508, 20037508.34),       allOverlays: false,       theme: null,       controls: []     });          // Add control     pMap.addControl(new OpenLayers.Control.LayerSwitcher());     geocoder = new OpenLayers.Control.Geocoder();     pMap.addControl(geocoder);          // Create gplayer     gplayer = new OpenLayers.Layer.Google(        "Google Physical",        { type: G_PHYSICAL_MAP,  isBaseLayer: true,sphericalMercator: true, minZoomLevel: 0 }     );     pMap.addLayer(gplayer);          // Create markers layer     markers = new OpenLayers.Layer.Markers( "Sites" );         var epsg4326 = new OpenLayers.Projection('EPSG:4326');     var epsg900913 = new OpenLayers.Projection('EPSG:900913');               var size = new OpenLayers.Size(21,25);     var offset = new OpenLayers.Pixel(-(size.w/2), -size.h);     var icon = new OpenLayers.Icon('http://www.openlayers.org/dev/img/marker.png',size,offset);                          markers.addMarker(new OpenLayers.Marker(new OpenLayers.LonLat(-88.63023333,36.00091667).transform(epsg4326, epsg900913),icon));                 var icon1 = icon.clone();         markers.addMarker(new OpenLayers.Marker(new OpenLayers.LonLat(-73.66586,42.297826).transform(epsg4326, epsg900913),icon1));         var icon2 = icon.clone();         markers.addMarker(new OpenLayers.Marker(new OpenLayers.LonLat(-110.802188,38.17243402).transform(epsg4326, epsg900913),icon2));         var icon3 = icon.clone();         markers.addMarker(new OpenLayers.Marker(new OpenLayers.LonLat(-110.9660481,38.43077243).transform(epsg4326, epsg900913),icon3));                  pMap.addLayer(markers);              pMap.zoomToMaxExtent();              selectControl = new OpenLayers.Control.SelectFeature(                 markers,                 {onSelect: onFeatureSelect, onUnselect: onFeatureUnselect});              pMap.addControl(selectControl);         selectControl.activate();        }          function onPopupClose(evt) {         selectControl.unselect(selectedFeature);     }          function onFeatureSelect(feature) {         console.log("We are at select.");         var selectedFeature = feature;              var popup = new OpenLayers.Popup.FramedCloud("SitesPopup",              //new OpenLayers.LonLat(-98.63023333,36.00091667).transform(epsg4326, epsg900913),             feature.geometry.getBounds().getCenterLonLat(),             new OpenLayers.Size(100,100),             //"<h2>"+feature.attributes.name + "</h2>" + feature.attributes.description,             "Test a popup",             null, true, onPopupClose         );         feature.popup = popup;         pMap.addPopup(popup);     }          function onFeatureUnselect(feature) {         console.log("We are at unselect.")         if(feature.popup) {             pMap.removePopup(feature.popup);             feature.popup.destroy();             feature.popup = null;         }     }


title: Calling Web Geoprocessing Services From Other Programs

text: Where can I go to find information about calling ArcGIS Server geoprocessing services from other programs like a .NET Framework program or possibly even a SQL Server job? So far, I've found Using ESRI geoprocessing services from a non-ESRI client (Python)? which seems to indicate that I'll have to put together a program to call the REST interface from scratch.


title: Differences between Arc 10.0 and Arc 10.1 Layer files

text: We have a program that cannot import Layerfiles from 10.1, but instead changes them to 10.0 layer files. I'm wondering if there is any Layer functionality that is lost in this transition. Right now only a couple of people are using 10.1 in our organization and they have been saving their layer files as 10.0, no issues. But going forward if we all change to 10.1 will this automatic saving to a 10.0 version of a layer file lose anything that we would notice?


title: Polygon grid shape file for landmass

text: I am interested in generating a polygon grid shape file for the landmass of the world. What I have presently is two tables in postgreSQL. Table 1: Coastal area shapefile (http://openstreetmapdata.com/data/coastlines) : Consisits od Linestrings column names : gid, fid, geom Table 2: Polygon grid shape file created for the above Table 1 with the help of QGIS Grid is (1 lat x 1 lon) in size.   column names : gid_grid, id_grid, coord_grid, geom_grid Now I want to find the polygon grids in Table 2 which are in the land area (even if a small part of coastal line falls inside the grid, it has to be marked as land area). My initial thought is to find the polygon grids which intercept with any of the coastal lines, and mark it as a land area (using postgis queries). The problem is, it will not find any polygon grids which totally falls inside the land area. Any suggestions please?


title: Print A0 sized map

text: Currently I am exploring an option/software to print A0 sized maps on web and also on desktop, can you guide me if any software support this printing option. Note: ESRI will be expensive therefore I need an alternative, I have Autodesk Map guide studio 2014 license.


title: How do I make vectorlayer with sql on PostGis via python code

text: I need to make a vectorLayer and fill with the SQL's results on PostGIS via python code. I'm working with python 2.7, QGIS 1.8 and PyQt 4.5 How do I it ? thanks


title: Using ArcGIS to filter cells where 7 or 8 of the neighbors are different

text: I have land cover data that I would like to clean using some generalization tools. Majority filter does almost exactly what I want, but is missing one key option. Majority filter allows you to change cells where 5 out of 8 of the neighbors are a different class. The problem is that in my data there are many instances when 5 out of 8 is too aggressive. I would like to only target those cells that are surrounded by 7 or 8 different neighbors. Is there a way to do this in Arc?


title: finding a location on a line using a point and a relative distance with arcgis js api

text: The problem I'm trying to solve is to find a point on a line given another point on the same line and a distance from the same point. I've been looking on the various RouteTask examples but they mostly tackle the problem of finding a route given two points. Can anyone advise if RouteTask is the right approach for this kind of problem, and how do I do it? Thanks!


title: Parsing a GeoJSON file with jQuery

text:
I'm trying to iterate through a GeoJSON file (below) and eventually populate an array with a few attributes within "properties". I keep coming up short working through various examples I can find. Do I need to embed another $.each() to get to properties. If it isn't obvious I'm new to this and I've hit the wall. What I have so far:               $(document).ready(function () {                     $.getJSON('testpoint.geojson', function (data) {                     var items = [];                     $.each(data.features, function (key, val) {                         items.push('
' + val + '
');                     });                     $('
', {                         'class':'my-new-list',                         html:items.join('')                     }).appendTo('body');                 });             });      And my json is as follows:               {         "type": "FeatureCollection",         "features": [             {                 "type": "Feature",                 "properties": {                     "gid": 1,                     "PID": "9-18-3",                     "BCT": "BCT126",                     "OWNER_TYPE": "A",                     "LOCNO": 0,                     "LOCEXT": "",                     "STREET": "CROSBY LANE",                     "ACQUIRED": "5/7/2010",                     "GRANTOR": "John A. SPARGO",                     "UPLAND": 0,                     "WETLAND": 3.96,                     "TOTAL": 3.96,                     "HABITAT": "salt marsh"                 },                 "geometry": {                     "type": "Point",                     "coordinates": [                         -70.03209,                         41.78278                     ]                 }             }         ]     }


title: Clickable points with openlayers

text: I'm trying to have some very basic functionality with vector points with Openlayers 2.12 and I can't figure out how to track a click event on a point. I create a map using the following code and add one point:               var map;     function init(){         map = new OpenLayers.Map('map');              var gsat = new OpenLayers.Layer.Google(             "Google Satellite",             {type: google.maps.MapTypeId.SATELLITE, numZoomLevels: 22}         );                      map.addLayer(gsat);              var vectorLayer = new OpenLayers.Layer.Vector("vector");              vectorLayer.addFeatures(             new OpenLayers.Feature.Vector(                 new OpenLayers.Geometry.Point(-57.635021, -25.276987).transform(                 new OpenLayers.Projection("EPSG:4326"),                 map.getProjectionObject()                 )             )         )               map.setCenter(new OpenLayers.LonLat(10.2, 48.9).transform(             new OpenLayers.Projection("EPSG:4326"),             map.getProjectionObject()             ), 2         );                      map.addLayer(vectorLayer);      Then, I've tried to identify when the point is clicked:                   // detect when points are clicked         pointClick = new OpenLayers.Control.GetFeature(vectorLayer, OpenLayers.Handler.Point,          {eventListeners:{"click": clickNotice}};                                                     map.addControl(pointClick);              // detect point events          function clickNotice(evt) {             alert('Point clicked');             point.deactivate();          }                    }      This is not working, (although I don't get any errors in Firebug), but I'm not sure of the correct syntax for handling points.


title: Bathymetry of Quarry Lake in QGIS

text: I've been asked to do a project that maps the bathymetry of a quarry lake. It seems like I've seen a feature in QGIS that would help with this if I took a bunch of depth points in the lake. Any thoughts on where to find this? And how to do it? And there's not a chance this data could already exist somewhere from remote sensing, soil surveys, etc.? Thanks