Show Fields in infowindow from my custom layer

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


"Add Comment" functionality in Flex Viewer

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?


R Raster package loop over pixels

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


Simplifying TopoJSON causes leak of Natural Earth in d3

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.


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

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.


Dynamically have symbology change as a due date nears?

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?


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

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.


Selecting MAX values within records with the same name

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


How to create a movable point in OpenLayers?

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?


Change zIndex of Vectors in openlayers

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


Pywps Timeout long processes

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


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

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.


Indian Latitude Longitude values

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.


query / filter on attribute tables with external tables joined

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


Is there a parellel line function for slds?

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


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

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!


Network Analyst not working

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)


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

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.


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

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.


Are the public domain vehicle tracks available?

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?


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

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


How Could I Calculate Shortest Path ?

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


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

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?


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

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?


draw shp files or kml files on ipad?

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.


Raster Combine - ESRI Spatial Analyst

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


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

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


Need information about the Envisat satellite

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?


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

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


GDB maintenance of muliple editors without ArcSDE

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


How to convert CAD Annotation to Shapefile with attributes?

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.


QGIS 2.0.1. missing Google Mercator CRS

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?


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

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?


GLONASS Satellite Position Prediction using Almanac Data

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.


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

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?


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

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?


Error Reporting Registry Location of ArcGIS for Desktop?

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?


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

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.


What is GIS for Indigenous Peoples?

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


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

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.


Problem compiling UI file generated by QGIS Plugin Buidler

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


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

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.


Visualizing Spatial Query

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


Assigning a nearest polygon to a point

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?


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

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?


Best way to train GIS professional working in Telecom

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.


Buffer aggregation method in OpenCarto library

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


Installing pgRouting in Windows 7

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.


Web Mapping Services from Europe besides OSM

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?


Clipping soil survey map service to feature class

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?


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

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?


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

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


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

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


Display information in infowindow using KMLLayer and kmz file

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>


Is Landsat 4-5-7 L1G geometric corrected?

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?


convert from HDF to GeoTiff

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)


Popups over markers not working in OpenLayers

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


Calling Web Geoprocessing Services From Other Programs

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.


Differences between Arc 10.0 and Arc 10.1 Layer files

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?


Polygon grid shape file for landmass

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?


Print A0 sized map

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.


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

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


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

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?


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

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!


Parsing a GeoJSON file with jQuery

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


Clickable points with openlayers

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.


Bathymetry of Quarry Lake in QGIS

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


How to increase mask buffer and then offset it using r.mapcalc in GRASS?

I'm in the process of learning on how to remove clouds and cloud-shadows from Landsat images using GRASS. So far I was able to remove a lot of clouds from those images using `i.landsat.acca`, but it seems like they still require a lot of refinements: ![Removing Landsat 7 clouds](http://farm9.staticflickr.com/8335/8113597796_9c1d97f120_b.jpg) I followed the GRASS manual for `i.landsat.acca` to obtain the image above: http://grass.osgeo.org/manuals/html70_user/i.landsat.acca.html. Basically, the process is to run different Landsat band images for the same scene with `i.landsat.toar`, and then processed with `i.landsat.acca` to get the cloud areas. Then use the cloud area as a mask using `r.mapcalc`:               r.mapcalc "MASK = if(isnull(clouds))"      **Problem#1** : the cloud edges is still visible. I need to increase the cloud-mask buffer to catch those renegade cloud edges. **Problem#2** : the cloud-shadows is still visible. If I include the `-s` flag to include a category for cloud-shadow when running `i.landsat.acca`, a large amounts of shadow from the terrain will consequently be removed as well... So to catch the cloud-shadows, I believe I need to make a new mask for them - by offsetting the current cloud mask until they cover the cloud-shadow areas. To achieve both objectives, I believe I need to modify and play with the `r.mapcalc` command, but I can't figure out how to write it (due to my bad Math and my-even-worse programming skills).. So, my ultimate question: _**Question#1_** : How do I edit `r.map.calc` to increase the mask area (let's say by a buffer of 10 units)? **Question#2** : How do we offset a mask in `r.map.calc` (in this case a little to the south-west direction)? If I know the answer to both of these: I can combine both cloud-mask and cloud-shadow-mask, to get a relatively cloud-free image. Then it's a matter of overlaying it to images of the same scene of different dates. I appreciate any help that I can get. Many thanks!


ArcGIS 10 Python script to search directories and subdirectories for intersects and clip them

I'm trying to develop a tool that queries our company's data stores for any data that intersects a feature. So when we acquire a new property we can run this tool and it will find all of the data we have that intersects the new property, clip it, and save it in a new location. Any help would be greatly appreciated.


Can anyone point me to a dataset (shapefile or other format) for the North American Continental Divide?

I'm looking for a more detailed dataset than the 1:2,000,000 shapefile available from the National Atlas, preferably one suitable for use at 1:24,000.


How to query latitudinal bounds on 95% of a polygon in PostGIS?

I'm going to say thanks in advance. I'm very new to GIS and this forum is a tremendous help. My simple questions have been solved; here's one that might be a little trickier. Given a polygon in a PostGIS geometry table, I want to find the northern and southern latitudes that would bound an arbitrary portion of the polygon's area (say 95%), with the remainder trimmed equally from the north and the south. I think what I want may be called a bounding box, though for a large polygon it might seem more like a conic section. I've noticed a number of functions relating to boxes in the documentation, but haven't figured out how they could be defined as I've described.


Looking for fastest solution for Point in Polygon analysis of 200 million points

Apologies for asking a question that has been asked before -- I've read through answers, and I'm left wondering how to implement the recommendations. I've never used SQL, Python, C, and have only used ArcGIS once before -- I'm a complete beginner. I have a CSV containing 200 million observations with the following format:               id,x1,y1,x2,y2,day,color     1,"-105.4652334","39.2586939","-105.4321296","39.2236632","Monday","Black"     2,"-105.3224523","39.1323299","-105.4439944","39.3352235","Tuesday","Green"     3,"-104.4233452","39.0234355","-105.4643990","39.1223435","Wednesday","Blue"      For each set of coordinates (x1/y1 and x2/y2), I want to assign the US Census Tract or Census Block that it falls within (I downloaded the Census tract TIGER shapefile here: ftp://ftp2.census.gov/geo/tiger/TIGER2011/TRACT/tl_2011_08_tract.zip). So, I need to do a point-in-polygon operation twice for each observation. It's important that the matches be very accurate. What is the fastest way to do this, including time to learn the software? I have access to a computer with 48GB of Memory -- in case that might be a relevant constraint. Several threads recommend using PostGIS or Spatialite (Spatialite looks easier to use -- but is it as efficient as PostGIS?). If those are the best options, is it imperative to populate an Spatial Index (RTree?)? If so, how does one do so (e.g. using the Census Tract Shape File)? I would be extremely grateful for any recommendations that include example code (or a pointer to example code). My first attempt (before finding this site) consisted of using ArcGIS to do a spatial join (x1/y1 only) of subsample of the data (100,000 points) on US Census Block. That took over 5 hours before I killed the process. I'm hoping for a solution that can be implemented on the entire dataset in under 40 hours of computing time. Thank you!


How to color markers based on polygon area?

I have a map with around 317 markers and each marker is the center point of a polygon. The polygonal area varies throughout, so I want to represent the markers in such a way that the polygon with the least area has the lightest shade of the marker color and the marker representing the largest area has the strongest shade of the same marker color. Can I dynamically produce such colored icons on google maps api??


Join attribute table

I am using QGIS 1.7.3 on Mac and when I try to use the tool "join" in the properties toolbox, I have no possibility to add a name for either the file or the field used to join tables. Has anyone encountered this problem and solved it ???


arcobjects update geometry Z values of line vertices with a given number

I am trying to develop arcgis 10.0 addin to update Z values of vertices of a selected line. Scenario: The tool scans all features in a given shape file and populates the list of features which has at least one vertex with a Z vlaue = 0. User then selects a feature from the list and clicks a button to update Z values to a given height - this last part does not work. I have tried the following for updating Z values of a selected polyline feature with no luck:                   Dim pFeatureSelection As IFeatureSelection                 pFeatureSelection = pFLayer                 Dim pEnumID As IEnumIDs                 pEnumID = pFeatureSelection.SelectionSet.IDs                 Dim id As Long                 id = pEnumID.Next                 Dim pPointCollection As IPointCollection                 Dim pFeature As IFeature                 Do While Not id = -1                     pFeature = pFLayer.FeatureClass.GetFeature(id)                     pPointCollection = pFeature.Shape                          For lPnt = 0 To (pPointCollection.PointCount - 1)                          SetZValueOnPoint(pPointCollection.Point(lPnt), Convert.ToDouble(txtZValue.Text))                          Next                          id = pEnumID.Next                 Loop      Please your help in providing any hint or suggesting error in the code above is much appreciated! many thanks meryloo


ajax not working with postgresql

I know is a silly question but I thought I should ask since you guys know more about postgresql. Sorry for being analytical, but I think its the only way to make my question clear. I'm using postgresql 9.1 /postGIS 2.0 /openlayers/geoserver and apache I've managed to get a feature's fid when clicked from a vector layer. I edited it using javascript and kept only the number, for example 64. Then I set this as a value of a hidden html form. Using ajax I m trying to pass this value to a PHP value and perform a query. Openlayers, javascript and ajax work fine but php never gets the value. Here is my code                           function selected_feature(event){     var sf = event.feature.fid;         var sfs = sf.split(".");         //that's the number out of fid         var sff = sfs[1];         //feed it to the form         document.getElementById("blahhh").value=sff;         //alert it-works         alert(document.getElementById("blahhh").value);         //call ajax         loadXMLDoc();         }                   <!--the form-->         <form id="hifo" name="hifo" method="post">         <input type="hidden" id="blahhh" name="blahhh" />         </form>                   //outside openlayers code                 function loadXMLDoc(){             var xmlhttp2;             alert("inside");//works             var t2lFname=document.getElementById("blahhh").value;             var parb="blahhh="+t2lFname;                  if (window.XMLHttpRequest)               {// code for IE7+, Firefox, Chrome, Opera, Safari               xmlhttp2=new XMLHttpRequest();               }             else               {// code for IE6, IE5               xmlhttp2=new ActiveXObject("Microsoft.XMLHTTP");               }             xmlhttp2.open("POST","map.php",true);             xmlhttp2.setRequestHeader("Content-type", "application/x-www-form-urlencoded");             xmlhttp2.send(parb);                  xmlhttp2.onreadystatechange=function()               {               if (xmlhttp2.readyState==4 && xmlhttp2.status==200)                 {             alert("suces");//works                 }               }                    //and then php     $blah = $_POST['parb'];      echo $blah;      $queryd='       SELECT        pins.p_name        FROM         pins         WHERE         pins.p_id='.$blah;      //fetching to fill array      $resultd=pg_query($conn, $queryd);      $nord=pg_num_rows($resultd);      $d=0;       while($arrayd=pg_fetch_array($resultd)){           $name[$d]=$arrayd['p_name'];         $d++;             echo $arrayd['p_name'];                     echo $arrayd[$d];             echo $d;         }      So basically use the fid as an id, via ajax to query the db. All this code is in the same file (named map.php).Problem is php's echos not working and query not working at all. I ve searched the web, I red tutorials and still no luck. Everything works except php. What should I do? Please advise Thanks PS I also tried OpenLayers.Request, found here and still no luck


XY event layer doesn't work with ArcGIS Server 10.1

I have a geoprocessing tool written in Python that contains operations on an XY event layer. It worked perfectly when running from ArcCatalog. When published as geoprocessing service on ArcGIS Server 10.1, however, it never worked. The JSON response when executed as REST service was               "error": {       "code": 400,       "message": "Unable to complete operation.",       "details": [        "Error executing tool.: ERROR 999999: Error executing function.\nThe XY event source name has not been setup correctly.\nFailed to execute (CopyFeatures).\n\nFailed to execute (NewCuringModelExpPy).\nFailed to execute (New_Curing_Model_ExpPy)."       ]      }      I have been struggling with this for days. Has anyone else ever come across this issue before or how was that solved? Thanks! Alex


How to cut a shape file?

I want to use a shape file of the Corine Land Cover. This covers all of Europe, but I am only interested in the area of Scandinavia. How can I cut this file using free software (preferably a commandline tool)? I tried _ogr2ogr_               C:\Program Files\GDAL>ogr2ogr.exe -f "ESRI Shapefile" c:\Users\Christoph\geo\corine\clc06_c312_coniferous\clc06_c312_dalarna.shp c:\Users\Christoph\geo\corine\clc06_c312_coniferous\clc06_c312.shp -clipdst 0 60 20 70      but without any usable result (just almost empty files).


finding nearest point to a point

I work in GIS with VBA. I have 2 Points layer.One of them is City points and another is hydro station points.I want write code in vba for selecting one point from station layer and finding closest city to this selected point.On the other hand I should find nearest point to a selected point. Please help me.


ArcGIS REST export map features as geoJSON

Is it possible to export the features from a map as geoJSON using an ArcGIS REST service? I'm trying to make a vector layer in OpenLayers based on the output from the map's present view based on a layer.


Parameters for ED50 -> ETRS89/WGS84 seven-parameter transformation which are good for Denmark?

I am looking for parameters for a seven-parameter transformation (Rx, Ry, Rz, Tx, Ty, Tz, s) from ED50 UTM (zone 32U) to ETRS89 and/or WGS84. The best match I could find is here: http://geodesy.spacecenter.dk/~geod/REF/3drefdk.htm Unfortunately, it seems like these parameters produce coordinates which are off by some 50 meters or so. Ideally, the parameters should lead to a transformation within a couple of meters of what the online calculator at http://valdemar.kms.dk/trf/ produces.


Publish a Feature class on ArcDesktop with C#?

All, I have a feature class in ArcCatalog. I need to save the feature class to a Map and then publish the map to the Arcserver. All programatically and hopefully in C#. I can't seem to find any examples at ESRI or by googling. Any help and examples of this would be absolutely amazing. But just leading me down the right path is cool too. Thanks guys!


Get the lines between all points of a polygon in PostGis : avoid nested loop?

I am trying to get all the lines between the points of a polygon, by writing PostGreSQL functions. I have written the two following functions that do the job fine :   * the first function to get all the points from the polygon               CREATE OR REPLACE FUNCTION PointsFromPolygon(polygon geometry)      RETURNS SETOF geometry AS     $$     DECLARE         point geometry;     BEGIN     FOR point IN SELECT DISTINCT points.geom FROM ( SELECT (ST_DumpPoints(polygon)).* ) AS points LOOP         RETURN NEXT point;     END LOOP;     END;     $$     LANGUAGE plpgsql ;        * the second function to create all segments between these points.               CREATE OR REPLACE FUNCTION AllSegmentsFromPoints(polygon geometry)     RETURNS SETOF geometry AS     $$     DECLARE         point1 geometry;         point2 geometry;         i integer;     BEGIN     i:=1;         FOR point1 IN SELECT * FROM PointsFromPolygon(polygon) LOOP             FOR point2 IN SELECT * FROM PointsFromPolygon(polygon) OFFSET i LOOP                 RETURN NEXT ST_MakeLine(point1,point2);             END LOOP;         i:= i+1;         END LOOP;     END;     $$     LANGUAGE plpgsql;      Is there a better way to perform this ? Perhaps by doing a cartesian product between the points ?


What Professional / Industry Associations are available that are related to GIS?

This was a well liked (but not hugely rated) question in the beta, so here goes. Please provide a brief synopsis and links if possible. I'd also like to hear about things that are in the planning stage. What group or association do you belong to? What are you working towards? Why?


Prevent layer from WMSGetFeatureInfo request on openlayers web map

I use this example code get feature info from my openlayers web map.I add 4 layers to map including raster image as base map. But i don't need feature info from that raster layer.But when clicking on the map to get information It gives information of all layers. How I modify above mentioned example code to prevent information from raster layer? Pls help me,Any idea welcome.Thanks


looking to integrate R-trees into AlchemyDB

I am the author of AlchemyDB: http://code.google.com/p/alchemydatabase/ I would very much like to add geospatial indexing in the form of R-trees (most likely) and need some advice on whether or not there any C/C++ packages that have simple stripped down R-trees. I am a big fan of small codebases that are easy to read. If anyone has any suggestions or knows of any experts that I could bug (I can pay for their time), I would be grateful


Publish PostGIS raster to geoserver

I have a database of raster and vector and I want to use them in my webmapping. There is not a problem to publish the vector using geoserver but I can't publish my raster. Is there possibility to publish postgis raster to geoserver?


How to print gwc layer?

I have gwc layer on my Geoserver.   Now i want print this layer using Geoserver print module but dont know how to do this.   When i try wms layer all works fine. But how to print gwc layer? **MORE INFO** In cliet side i have a gxp and try to print using gxp_print plugin. My problem in that i cant set url to gwc layer.   I read this and this: so i try this url:               layer.setUrl("/geoserver/gwc/service/wms?BBOX=7637927.0361051%2C6845332.4516735%2C7750900.9638949%2C6971453.5483265&HEIGHT=256&     WIDTH=256&TRANSPARENT=true&map_resolution=305.74811309814453&VERSION=1.1.1&     FORMAT=image%2Fpng&TILED=true&SERVICE=WMS&format_options=dpi%3A100&LAYERS=gisurfo:atd&     REQUEST=GetMap&STYLES=&SRS=EPSG%3A900913");      But geoserver return me error:               400: Requested horizontal resolution: 441.30440542890574 , best match:      305.74811309814453 exceeds 10% threshold. Perhaps the client is configured with an      incorrect set of scales (resolutions), or the DPI setting is off compared to the one in GWC ?      This same like in links above. So how to print gwc? **UPDATE** Good news. I change scales in my `config.yaml` file to scales what i get from gwc layer (but i get only left part of each scale). Now i try to print gwc layer but it not printing. In Geosrever log i found Bad Request:               java.io.IOException: Error (status=400) while reading the image from      http://geourfo.ru:9091/geoserver/gwc/service/wms?BBOX=6866019.0%2C7636631.0%2C6950767.0%2C7752196.5     &HEIGHT=625&WIDTH=458&VERSION=1.1.1&FORMAT=image%2Fpng&SERVICE=WMS&     format_options=dpi%3A75&LAYERS=gisurfo%3Aatd&REQUEST=GetMap&STYLES=&SRS=EPSG%3A900913     : Bad Request      After i change width and height to 256/256 and get my tile!   So question: why iprint plugin take tiles with HEIGHT=625 and WIDTH=458 but not with 256x256?


Can i do Custom Styling on Geoserver

i want to use Geoserver WMS server for custom styling. lets say i have an area as polygon in a shape file. But i dont have this area's temperature as an attribute. From a web service i gather temperature points inside that area. Now i want Geoserver to color this area according to these dynamic temperature points here is an example ![Temperature Map](http://i.stack.imgur.com/xJ1U8.png) if it is possible, i even want to place these borders on image and some explanations about colors on right side. How can i send these temperature points to geoserver on wms. And can i do some custom coloring according to these dynamic parameters? I dont know how to convert these point data to area colors. My next question is similar but without points, can i colorize a polygon with dynamic attributes which isn't on shape file. here is an example i want to colorize these polygons according to their product type, but as earlier question, these product types are gathered dynamically ![enter image description here](http://i.stack.imgur.com/P9CCy.png) I hope these are possible


IMapControlEvents2.OnKeyDown not hit on Alt + key

I want to catch the Alt key plus an addionally hit key in the MapControl's `OnKeyDown`. I did already set MapControl's `KeyIntercept` to `esriKeyIntercept.esriKeyInterceptAlt` but the event handler doesn't get raised on Alt or Alt + key. It does get hit on other keys. Any suggestions?               _mapControl.KeyIntercept = (int)(esriKeyIntercept.esriKeyInterceptAlt);      Handler:               void _mapControl_OnKeyDown(object sender, IMapControlEvents2_OnKeyDownEvent e)     {       Trace.WriteLine("hit");     }


focal statistics with variable radius

is there an easy way to do focal statistics on a raster, with a variable radius for the neighbourhood to be searched? As in: that search radius would be stored in another raster? i tried something like               f_max = arcpy.sa.FocalStatistics(Tau,arcpy.sa.NbrCircle(rad_ras,"CELL"), "MAXIMUM","DATA")      it returns               Traceback (most recent call last):       File "K:\Informatik\tools\arcpy\GKPROZ\test_3.py", line 383, in <module>         ufereros_gq(r"F:\25sg\25107\TG7\access\gis_logger_tg7.mdb",gq_nummer, False)       File "K:\Informatik\tools\arcpy\GKPROZ\test_3.py", line 300, in ufereros_gq         f_max = arcpy.sa.FocalStatistics(Tau,arcpy.sa.NbrCircle(rad_ras,"CELL"), "MAXIMUM","DATA")       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Functions.py", line 4796, in FocalStatistics         ignore_nodata)       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Utils.py", line 47, in swapper         result = wrapper(*args, **kwargs)       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Functions.py", line 4783, in wrapper         neighborhood = Utils.compoundParameterToString(neighborhood, ParameterClasses._Neighborhood)       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Utils.py", line 75, in compoundParameterToString         return str(parameter)       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\sa\ParameterClasses.py", line 202, in __str__         self.units])       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\sa\CompoundParameter.py", line 39, in _toString         userProvidedPounds = [value == "#" for value in values]       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Functions.py", line 3598, in EqualTo         in_raster_or_constant2)       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Utils.py", line 47, in swapper         result = wrapper(*args, **kwargs)       File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\sa\Functions.py", line 3595, in wrapper         return _wrapLocalFunctionRaster(u"EqualTo_sa", ["EqualTo", in_raster_or_constant1, in_raster_or_constant2])     RuntimeError: ERROR 000732: Input Raster: Dataset # does not exist or is not supported      although the raster(s) exists. I assume this means that the Nbr function can't work with a raster, as the same line works fine for a fixed value for the radius of the circle. I have found a workaround in that i do the focal statistics 10 times with a different radius each time and then use arcpy.sa.Pick to choose from those 10. This is however painfully slow. Is there a better way to do it? any help would be appreciated. i use ArcGis 10.0 SP4 on an info license on Windows 7 64


Converting a shp file to a list of latitudes and longitudes

I have a shp file with the latitudes and longitudes for the boundaries for various geographic regions. I'd like to convert the shp file into lists of just latitudes and longitudes for the regions. I'll then use these lists in my program to determine which region a specific point resides in. Is there a quick way to do this conversion or will I have to read up on the APIs? Even converting the shp file into GML would be a help as I could easily access that programmatically.


add marker to map from list of coordinates

Forgive me if this question has an incredibly obvious answer. I'm very new to ESRI and I'm trying to add markers to my map via gps coordinates. I have not found an example that does something this simple. When I try to adapt the more complicated ones, my code breaks. Would this be an image layer? a SimpleMarkerSymbol? a PictureMarkerSymbol? What would really help is an extremely streamlined example of a marker on a map so I can see how it works and then adapt my code from there. It seems all the samples I'm finding have all kinds of fancy extras that I'm not familiar enough to play with. EDIT I've actually got some markers on my map! yay! But now, I'm having a tough time with the coordinates showing up where I expect. As you can see from the screenshot, the highlighted marker uses gps coordinates (42.10560 -88.380801) that should put it somewhere in northern Illinois but for some reason, its all the way off the east coast. The second marker shows up in Illinois but I guestimated the coordinates for that one. ![sandbox screen shot](http://i.stack.imgur.com/QDn7q.jpg) I think it has something to do with the SpatialReference class but I've tried changing the WKID attribute to various things that some of the mapping experts in my office have suggested but nothing seems to change anything. Here is my code for playing around with.               <!DOCTYPE html>                     <meta http-equiv="Content-Type" content="text/html; charset=utf-8">                  <meta name="viewport" content="initial-scale=1, maximum-scale=1,user-scalable=no">         <title> Topographic Map with graphics</title>     <link rel="stylesheet" href="http://js.arcgis.com/3.7/js/dojo/dijit/themes/claro/claro.css">     <link rel="stylesheet" href="http://js.arcgis.com/3.7/js/esri/css/esri.css">

      I must be missing something. Please help!


How to get around the 10-character limit when converting from kml to shp with ogr2ogr

I'm trying to create custom .shp files by drawing polygons on a new layer in Google Maps Engine, exporting that .kml file, and then converting that file to a .shp file with ogr2ogr. If I run a command like this: `ogr2ogr -f 'ESRI Shapefile' foo.shp foo.kml` I get the following printed to stdout: `Warning 6: Normalized/laundered field name: 'Description' to 'Descriptio'`. This causes a problem downstream because when I load the .shp file into a r-tree and then query it to perform reverse geocoding, I get back a field "Descriptio" instead of "Description", which just looks odd. Moreover, this happens despite the fact that when I open the .kml file in vim, I don't see **any** "description" tag anywhere. So my question is: is there any way around this? Can I alter the KML file in any way to prevent `ogr2ogr` from complaining that it needs to normalize/launder the field `Description`? Please let me know if this question is written poorly. It's my first time on GIS Stack Exchange and I'm not at all trying to be difficult/annoying/insufferable.


Postgis add Geography to function()?

i was wondering im getting different results when using geography() on different positions within a query. So, whats the most accurate way of forcing geography() over a functions with geometry columns? example: * * * SELECT polys.id, sum(st_length(st_intersection( **geography** (lines.the_geom), **geography** polys.the_geom))))/1000 total_lengths from schema.table1polys polys schema.table2lines lines on ST_Intersects(geography(lines.the_geom),geography(polys.the_geom)) * * * ..or.. * * * SELECT polys.id, sum(st_length( **geography** (st_intersection(lines.the_geom,polys.the_geom))))/1000 total_lengths from schema.table1polys polys schema.table2lines lines on ST_Intersects(geography(lines.the_geom),geography(polys.the_geom)) whats the most accurate? i have sort of 2 / 3 Km of difference between both querys.


Is it possible to link a QGis project with an excel file?

I have an Excel file with coordinates and I have used this file to plot the coordinates in QGis. I was wondering if it was possible to permanently link the excel file to the QGis project, so that when I add new coordinates to my excel file, these will show up in QGis when I open the project. Any ideas how to solve my problem?


How to improve Triangulated Irregular Network (TIN) display speed?

I am preparing Triangulated Irregular Network for large area (5000 Sq Km) in ArcGIS, The TIN models takes too much time to open in ArcMap, and if I move or change anything again its starts to rendering, It takes too much time to do small changes, how to stop the rendering, and how to keep it in Builded TIN in ArcMap. Number of Nodes in TIN is 3633846 Number of Triangles is 7267645 it takes three minutes to open fully.


Definite way to confirm a geocode as correct

So you tested few geocoding services for the same street address. You got results varying by 0.7 to 3 kilometers. How can you tell which geocoding service was more right, other than letting your user mark their real location, compared to where you show they are, on a map? feels a bit like a chicken and egg problem, unless you let the user mark their location on the map.


How do I convert a QGIS 1.8 plugin to QGIS2.0?

I have developed a plugin under QGIS1.8 and now I want to make it work on QGIS2.0 I have chage qgisMinimumVersion into 2.0(originally 1.8) but QGIS2.0 still give me the error say               This plugin is incompatible with this version of QGIS     Plugin designed for QGIS 1.x      how to fix it??