title: Openlayers Base Layer is disappeared after adding more than 15 Overlays from Geoserver

text: I have developed a webmap application using Geoserver 2.4.0,Tomcat7 on Windows 2008 32bit Server. JVM Client. I have enabled GeoWebCache to speedup the map loading and configured separate directories for GeoWebCache and Geoserver data directory. In my application, I am loading GeoTiff images (Entire UK) as base layer 1st and then Postgis data (Vector data for plymouth) as overlays. Problem is When 20 Overlays are loaded, the base layer got disappeared. If I turnoff the 15 Overlay's visibility, then the base layer came to visible. I am unable to findout the problem as I am a newbie. Randomly I kept 4 to 5 Overlays On to see that base layer. i.e. I feel that there is no issue in Overlays because I can view all those layers (base and overlays) in Geoserver preview as well as GeoWebCache demo page. but is there any limitation to adding number of Overlays? 20 Overlays means is it exceed that limitation? Please help share your suggestions/solutions. It will help me lot to fix this issue. Many Thanks in Advance. Here is my code for adding Base and Overlays: Base Layers:                geoLays = new OpenLayers.Layer.WMS(                             ODC Overview,                             "http://" + geoIpAddress + ":8080/geoserver/gwc/service/wms",                             {                                 LAYERS: cite:ODC Overview,                                 STYLES: '',                                 format: 'image/png'                             },                             {                                 tileSize: new OpenLayers.Size(256, 256),                                 isBaseLayer: true                                                 }                         );                     osMap.addLayers([geoLays]);       Overlays:                 for (var i = 0; i < overlays.length; i++) {              var geoLayName = overlays[i].split(':');             var workSpace = geoLayName[0];             var layerName = geoLayName[1];             geoLays = new OpenLayers.Layer.WMS(                             layerName,                             "http://" + geoIpAddress + ":8080/geoserver/gwc/service/wms",                             {                                    LAYERS: overlays[i],                                 STYLES: '',                                 format: 'image/png'                             },                             {                                 tileSize: new OpenLayers.Size(256, 256),                                 isBaseLayer: false                             }                         );                                         osMap.addLayers([geoLays]);      }      I found some error on Geoserver log but I am not sure due to this only it's happend The error is: Request: getServiceInfo 2014-03-28 16:41:00,103 INFO [geoserver.wms] - Request: getLegendGraphic Strict = false Format = image/png Rules = [] Locale = null Height = 20 Width = 20 Layers = [null] Styles = [StyleImpl[ name=green]] FeatureType = null Env = {} Transparent = false Scale = -1.0 Exceptions = SE_XML LegendOptions = {} Version = 1.0.0 Request = GetLegendGraphic Get = false RawKvp = {STYLE=green, VERSION=1.0.0, STRICT=false, FORMAT=image/png, HEIGHT=20, REQUEST=GetLegendGraphic, WIDTH=20} BaseUrl = http://xxx.xx.xxx.xx:8080/geoserver/ RequestCharset = null 2014-03-28 16:41:00,118 INFO [geoserver.flow] - Running requests: 0, processing through flow controllers: 0 2014-03-28 16:41:10,290 WARN [wicket.Localizer] - Tried to retrieve a localized string for a component that has not yet been added to the page. This can sometimes lead to an invalid or no localized resource returned. Make sure you are not calling Component#getString() inside your Component's constructor. Offending component: [MarkupContainer [Component id = subform]] 2014-03-28 16:41:10,290 WARN [wicket.Localizer] - Tried to retrieve a localized string for a component that has not yet been added to the page. This can sometimes lead to an invalid or no localized resource returned. Make sure you are not calling Component#getString() inside your Component's constructor. Offending component: [MarkupContainer [Component id = subform]] 2014-03-28 16:41:20,868 ERROR [io.IObjectStreamFactory$DefaultObjectStreamFactory] - error writing object [Page class = org.geoserver.web.data.resource.ResourceConfigurationPage, id = 9, version = 1, ajax = 0]: null java.lang.NullPointerException at org.geoserver.gwc.layer.GeoServerTileLayerInfoImpl.getParameterFilters(GeoServerTileLayerInfoImpl.java:360) at org.geoserver.gwc.layer.GeoServerTileLayerInfoImpl.writeReplace(GeoServerTileLayerInfoImpl.java:153) at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method) at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source) at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source) at java.lang.reflect.Method.invoke(Unknown Source) at java.io.ObjectStreamClass.invokeWriteReplace(Unknown Source) at java.io.ObjectOutputStream.writeObject0(Unknown Source)   2014-03-28 16:41:20,868 ERROR [lang.Objects] - Error serializing object class org.geoserver.web.data.resource.ResourceConfigurationPage [object=[Page class = org.geoserver.web.data.resource.ResourceConfigurationPage, id = 9, version = 1, ajax = 0]] java.lang.NullPointerException


title: ESRI REST Api - need to download batches of tiles

text: So I wanted to use the ESRI REST Api to fetch map tiles for an application. Right now I am calling for each tile as I need it via its X/Y and zoom level(which is generating a bunch of http calls), but I'm wondering if there is a way to call for a batch of tiles in a single call. I've been looking through the API documentation but hadn't found a command to get more than a single tile, I was wondering if anyone with more experience might know how.


title: Arcgis9.3 Problem with search cursor and update cursor in two different table

text: I am making a python code for ArcGIS 9.3 : -I have one folder with a lot of shapefile polygon inside \- Ech polygon shapefile have a field "Average" and one ID \- I have a table with each field corresponding to the ID of shapefile I want to open each shapefile, take the "Average" and put it in the good field. So for shapefil 1, I take all the averages and i put them in line 1 of the table. To do it I ha one SearchCursor in my shapefile and one UpdateCursor in my table. But it doesn't work and it seems to come from the update function. The code is as following and the error come after "Marqueur 6.1" which is link to the SetValue method and to the update cursor. The mesage error is : Marqueur 6.1 : ERROR 999999: Error executing function. Failed to execute (OpeTable2).               for infile in glob.glob(os.path.join(path, '*.shp')):              MyShape = os.path.join(path, infile)         Name = os.path.basename(MyShape)         fileName, fileExtension = os.path.splitext(Name)              gp.AddMessage(Name)              if TestFieldCreate == False :             TestFieldCreate = True             rowsP = gp.UpdateCursor(MyShape)             rowP = rowsP.Next()             #row est une ligne et non une colonne             while rowP:                 ID_BVs = rowP.GetValue(Field_ID_BV)                  ID_BV_STR = str(ID_BVs)   #Converti en string                 NameBV = "BV"+ ID_BV_STR                 rowsP.updateRow(rowP)                 rowP = rowsP.Next()                 gp.addfield (MyTable, NameBV, "double")                  #On place le curseur dans le shp              rowsInShape = gp.SearchCursor(MyShape)         rowInShape = rowsInShape.Next()              while rowInShape:                  countStr =str(count)             gp.AddMessage( "Iteration N°" + countStr)             Moyenne = rowInShape.GetValue(Field_Moyenne)                  if Moyenne != 0 :                 gp.AddMessage("Marqueur 6")                      # On récupère les valeur pour AfoisP et le num de BV que l'on converti de la même manière que le nom des champs                  BV = rowInShape.GetValue(Field_ID_BV)                  BV_STR = str(BV)   #Converti en string                 FieldName = "BV"+ BV_STR                 gp.AddMessage("Marqueur 6.1")                 rowInTable.SetValue(FieldName, Moyenne)                 rowInTable.SetValue("NomFich", fileName)                 rowsInTable.updateRow(rowInTable)                 gp.AddMessage("Marqueur 7")                  gp.AddMessage("Marqueur 7.1")             rowInSHape = rowsInShape.Next()             count=count+1              gp.AddMessage("Marqueur 7.2")         rowsInTable.updateRow(rowInTable)         rowInTable = rowsInTable.Next()         gp.AddMessage("Marqueur 8")         NumStr = str(Num)         gp.AddMessage(Name + " Fichier N° " + NumStr)         Num = Num +1         gp.AddMessage("Marqueur 9")      Thanks a lot


What is best the opensource webclient alternative To Geocortex essentials? For someone without much coding experience


title: Link to an individual info window in CartoDB?

text: I'm a new user, so thanks for your patience and assistance. I just published a very large parcel dataset on CartoDB. Another website would like to link their parcel dataset to ours. Meaning, if you click on a parcel on their site and want to know more information about that parcel that's published on our site, they can link directly to the individual parcel. That info is shown in a CartoDB info window on our site. Is there a way to link directly to that individual pop-up -- perhaps in the URL?


title: Problem Initializing License for ArcEngine application

text: I have written an ArcEngine application using vb.net and Visual Studio 2010. When I initially started the code, I indicated an ArcInfo license was available during the project setup. This worked fine for a while, and then one day for no apparent reason it stopped...but at the time I didn't need ArcInfo, I just needed ArcView, so I changed the LicenseInitializer to use ArcView and it kept on working, so I didn't give it any thought. (How I wish I'd written down what might have changed!) Now I do need ArcInfo for what I'm trying to do, and I've tried changing the requested license back using the Project->Add ArcGIS License Checking tool just in case I accidentally messed something up, but it won't work. What is particularly confusing is that the error message I'm getting says ArcInfo isn't initialized, but then if you scroll down, it says it is (see image). Note that I'm not trying to use any of the extensions, just ArcInfo itself. ![License Failure](http://i.stack.imgur.com/HA3tl.jpg) I've stepped through the code, and the InitializeApplication function on the Startup handler (code below) returns true. The error occurs when my form with an embedded AxMapControl actually attempts to use the license (BeginInit method on AxMapControl).               Private Sub MyApplication_Startup(sender As Object, e As Microsoft.VisualBasic.ApplicationServices.StartupEventArgs) Handles Me.Startup         'ESRI License Initializer generated code.         m_AOLicenseInitializer.InitializeApplication(New esriLicenseProductCode() {esriLicenseProductCode.esriLicenseProductCodeArcInfo}, _         New esriLicenseExtensionCode() {})     End Sub      I know that I have an ArcInfo license on this computer, and now that I've regenerated the license for my code with the License Checking Tool I'm confident I haven't messed with any normally hidden files I shouldn't be touching. Any ideas as to what might be causing this error??


title: How do I extract beginning/ending points for polyline shapes?

text: I need to extract the beginning/ending coordinates of every shape within a polyline shape file. It would be best to actually add them as values within the shape files .dbf (I can live without doing this). What is the easiest way to do this in ArcMap? I know that XTools has a plugin that does this really easy, but I don't want to purchase software for something that seems like it should be fairly trivial. Thank you! **Edit** I am currently using ArcMap 10.


title: How to change Geoserver Mapfish print option from pdf to gif

text: I am an absolute novice when it comes to Geoserver and Mapfish but have set myself the task of obtaining a GIF print output rather than PDF. From reading documentation and forums it seems this is possible, and I have a config.yaml with output formats defined (snippet below) but they are not being displayed when I go to my local Geoserver installation. Config.yaml #=========================================================================== # allowed DPIs #===========================================================================               dpis: [75, 150, 300, 600]          #===========================================================================     # the allowed scales     #===========================================================================     scales:       - 25000       - 50000       - 100000       - 200000       - 500000          #===========================================================================     # the list of allowed formats     #===========================================================================     outputFormats: ['pdf', 'png', 'gif']      JSON result from my localserver/geoserver/pdf/info.json {"scales":[{"name":"1:25,000","value":"25000"},{"name":"1:50,000","value":"50000"},{"name":"1:100,000","value":"100000"},{"name":"1:200,000","value":"200000"},{"name":"1:500,000","value":"500000"},{"name":"1:1,000,000","value":"1000000"},{"name":"1:2,000,000","value":"2000000"},{"name":"1:4,000,000","value":"4000000"}],"dpis":[{"name":"75","value":"75"},{"name":"150","value":"150"},{"name":"300","value":"300"},{"name":"600","value":"600"}],"layouts":[{"name":"A4 portrait","map":{"width":440,"height":483},"rotation":true},{"name":"Legal","map":{"width":440,"height":483},"rotation":false}],"printURL":"localserver/geoserver/pdf/print.pdf","createURL":localserver/geoserver/pdf/create.json"} From looking at the documentation at http://www.mapfish.org/doc/print/protocol.html I would expect these formats to be defined in this JSON string? I have not set about creating a request yet to see if I can print a GIF so any thoughts/pointers would be appreciated. The version of Geoserver is 2.2.5 patched to 2.2.xx with a 2.2-SNAPSHOT version of Mapfish. I have the demo from Geoserver working correctly and pdfs being generated, but really need GIFS otherwise MapFish is a none starter. Apologies in advance if I have incorrectly formatted parts or missed some vital piece of information out.


title: Can I use a PNG point marker?

text: I want to style a point layer from png file. Can i do this in quantum gis?


title: Custom Buttons for Google Drawing Manager

text: can we design my own custom button for the Drawing manager? According to API https://developers.google.com/maps/documentation/javascript/reference#DrawingManagerOptions we can customize only the symbols which is drawn out of the drawing manager. for example i need to make this drawing manager look like Google Map maker http://www.google.co.in/mapmaker


title: OSX 10.9 + QGIS 2.0.1 - Missing osgeo module

text: Running the following:   * `QGIS 2.0.1` (via Kyngchaos)   * `OSX 10.9`   * `GDAL/OGR 1.10.1` (via Homebrew). Recently, I've become unable to use the OpenLayers plugin, which is complaining that I'm missing the `osgeo module`. Has anyone else dealt with this? As always, thanks.


title: How to correctly set up Oracle connection parameters?

text: Just testing out the new QGIS2 install here and trying to connect it to our huge SDE/Oracle spatial database. Im struggling with the connection parameters. Here is the connection properties to the databse from within ArcGIS10.1:               Instance = sde:oracle10g:g01l10g      I can't figure out how this maps across to the connection fields in the QGIS2 popup. Can anyone help me please?


title: Developing QGIS Plugin on windows

text: I have developed a small plugin which acts as a plugin for QGIS. On linux, QGIS (customized) as well as the plugin were built using QT creator and worked well. Same code I am trying to port on to windows - I have built QGIS using Visual Studio and then built the plugin in Release using QT creator (mingw). When loading the plugin I get the following error. 'qgis.exe': Loaded 'F:\21stMarch2013\qgis_Build\output\plugins\RelWithDebInfo\RTools.dll', Binary was not built with debug information. 'qgis.exe': Unloaded 'F:\21stMarch2013\qgis_Build\output\plugins\RelWithDebInfo\RTools.dll' "Failed to load F:/21stMarch2013/qgis_Build/output/bin/RelWithDebInfo/../../plugins/RelWithDebInfo/RTools.dll (Reason: Cannot load library F:/21stMarch2013/qgis_Build/output/bin/RelWithDebInfo/../../plugins/RelWithDebInfo/RTools.dll: The specified module could not be found.)" The file 'RTools.dll' is very much present in the plugins directory. I have also copied all the dependencies into the plugins folder so as to make sure there is no issue with the dependencies. Any advice would be of great help. \--Srikanth


title: How to make Mapnik render Erdas images in RGB instead of grayscale?

text: I have some orthophoto files (Erdas Imagine format) enclosed in a vrt file. When I view it in any desktop GIS it has proper RGB colors, but when I use the vrt in Mapnik to render some tiles, I get grayscale images as a result. Why? Update: If I convert the Erdas Imagine files into png-s (I need transparency where there is no data) Mapnik render works ok, so I guess it's something to do with how Mapnik interprets Erdas files and its bands.


title: Problems with matching layers with different projections

text: I have two layers with the following projections: Monte Mario IT1 UTM* Anyhow, I followed the instructions (open the project, set the projection, select on the fly projection, load the layers but no way, the two layers does not match each other. I could only obtain the result of picture 1 where in red are the transport flows, in blue the transport network. The crs of the two layers are: red--> +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=1500000 +y_0=0 +ellps=intl +towgs84=-104.1,-49.1,-9.9,0.971,-2.917,0.714,-11.68 +units=m +no_defs blue--> +proj=tmerc +lat_0=0 +lon_0=9 +k=0.9996 +x_0=500053 +y_0=-3999820 +ellps=intl +towgs84=-87,-98,-121,0,0,0,0 +units=m +no_defs If I set one of the layer with the same crs as the other, the two networks becomes very far from each others (see picture 2) Which do you think might be the problem? Thank you, Selini ![in red transport flows in blue the transport network](http://i.stack.imgur.com/dzuA5.jpg) ![enter image description here](http://i.stack.imgur.com/wFklQ.jpg)


title: Find Direction/Location of point in polygon using Python ArcGIS10.1

text: I have a point located inside a polygon, I want to determine its location or direction i.e. if it is in the northwest corner or southwest. I have the x,y of the point. Any idea on how to do this? I have the shape file of the point and polygon. Any help on this is most welcome.


title: Is there a method to create offline geocoder using OSM xml data?

text: I am trying to develop an offline android map that makes use of OSM map tiles of 15 zoom level, For this purpose i used an open source library : Osmdroid Now, i m looking into the possibility of creating an offline geocoding/ reverse geocoding that can be integrated with my application, can i use Osm xml data for that purpose? if so , then can anyone suggest/explain how to use it to create SQlite db.. to be used wit my app While doing this, i hope i m not violating any data usage terms or conditions? Thanks


title: How to delete a map from the GRASS interface?

text: I want to delete a map that I generated from GRASS. I tried removing a map from the "Map Layers for each display" tab in the GRASS GIS Layer Manager, but it only removed the map from the list and not deleted. I want to delete the map totally that it won't appear as an option in when I load the maps in the workspace (Workspace>Load Map Layers). How is it done?


title: Copying Network Analysis Layer - Service Area Polygon in Python Script

text: I have built a geoprocessing tool using arcpy modules in python. Screenshot shown below: ![enter image description here](http://i.stack.imgur.com/SaAlx.jpg) In this tool I am creating Service Area Polygons and selecting census boundaries within those Service Area Polygons to extract the demographic data. The part where I am creating the service area polygon and copying them to output is shown below:               #PART 3: CREATING SERVICE AREA POLYGON USING DRIVING DISTANCE AS COST FACTOR AND INTERSECTING WITH CENSUS LAYER AND CREATING OUTPUT          Service_Area = "Service Area"          # Process: Make Service Area Layer     arcpy.na.MakeServiceAreaLayer(road_network, Service_Area, "Length", "", distance)          # Process: Add Facility to the Network     arcpy.na.AddLocations(Service_Area, "Facilities", "facility", "", "")          # Process: Solve the Network     arcpy.na.Solve(Service_Area)          # Copy Polygons to Service Area output     arcpy.CopyFeatures_management(Service_Area, service_area_output)          # Process: Create a new Feature layer - Census layer     arcpy.MakeFeatureLayer_management(cen_blk,"blocks_lyr")          # Process: Create a new Feature layer - Census layer     arcpy.MakeFeatureLayer_management(service_area_output,"service_area_output")          # Process: Select Layer By Location...     arcpy.SelectLayerByLocation_management("blocks_lyr", "HAVE_THEIR_CENTER_IN", "service_area_output")     if census_output=="#" or census_output=="":         census_output==None          else:         # Process: Copy the intersected census layer to output FC and add it into current workspace         arcpy.CopyFeatures_management("blocks_lyr", census_output)          gp.AddMessage (" ")      I cannot seem copy the service area polygon to a feature class. I am getting a following error message. Is CopyFeature not the right method to copy the service area polygons? ![enter image description here](http://i.stack.imgur.com/2bNvd.jpg) Thank you for your help in advance. Any feedback is appreciated.


title: Qgis - taudem plugin - python script not found

text: I installed the sextante toolbox on qgis. my os is fedora. grass works fine but i've got problems with taudem. it seems like some python files are missing. for instance when i try to run the pitremove algorithm the sextante info log yields: "HYDU_create_process (./utils/launch/launch.c:69): execvp error on file /home/me/.qgis/python/plugins/sextante/taudem/pitremove (No such file or directory)" indeed, there's no pitremove.py in the plugins/sextante/taudem directory. so did I missed one step in the plugin installation ? thanks for your help


title: Geotool SwtMapPane draws in view but not in TitleAreDialog

text: I've written a very basic rcp app to show points and lines in a view. When I tried to move this code to a wizard page, it gave the error below. When I use code that simply loads and displays a shapefile, it works for both cases. How do I resolve this issue, or what should I do differently? Jun 12, 2014 2:58:36 AM org.geotools.renderer.lite.StreamingRenderer fireErrorEvent SEVERE: org/geotools/filter/text/cql2/CQLException java.lang.NoClassDefFoundError: org/geotools/filter/text/cql2/CQLException My code:               public void createPartControl(Composite parent) { //returning main for TitleAreaD         final Composite main = new Composite(parent, SWT.BORDER | SWT.RESIZE);         main.setLayout(new GridLayout(1,false));         main.setLayoutData(new GridData(SWT.FILL, SWT.FILL, true, true));         Button b = new Button(main,SWT.NONE);         b.setText("load");         b.addSelectionListener(new SelectionListener() {                 @Override             public void widgetSelected(SelectionEvent e) {                 try {                     displayFeatures(mapPane,content);                 } catch (Exception e1) {                     e1.printStackTrace();                 }             }             @Override             public void widgetDefaultSelected(SelectionEvent e) {}         });              content = new MapContent();         content.getViewport().setBounds(new ReferencedEnvelope(-180,-90,180,90,DefaultGeographicCRS.WGS84));             mapPane = new SwtMapPane(main, SWT.NO_BACKGROUND | SWT.BORDER | SWT.DOUBLE_BUFFERED);         mapPane.setMapContent(content);         mapPane.setBackground(Display.getCurrent().getSystemColor(SWT.COLOR_WHITE));              StreamingRenderer renderer = new StreamingRenderer();         mapPane.setRenderer(renderer);     }          private void displayFeatures(SwtMapPane mapPane, MapContent map) throws Exception {         SimpleFeatureCollection features = createFeatureCollection();         Style style = createStyle("shape");         FeatureLayer layer = new FeatureLayer(features, style);         map.addLayer(layer);         mapPane.reset();         mapPane.redraw();     }          private SimpleFeatureCollection createFeatureCollection() throws Exception {         SimpleFeatureType TYPE = DataUtilities.createType("testtype","shape:Geometry:srid=4326, name:String");              List<SimpleFeature> features = new ArrayList<>();         SimpleFeatureCollection fc = new ListFeatureCollection(TYPE, features);         SimpleFeatureBuilder fb = new SimpleFeatureBuilder(TYPE);              Geometry g;         WKTReader reader = new WKTReader();              g = reader.read("POINT(10 10)");         features.add(fb.buildFeature(null, new Object[]{g, "p1"}));              g = reader.read("POINT(20 20)");         features.add(fb.buildFeature(null, new Object[]{g, "p2"}));              g = reader.read("POINT(30 -30)");         features.add(fb.buildFeature(null, new Object[]{g, "p3"}));              g = reader.read("POINT(40 -70)");         features.add(fb.buildFeature(null, new Object[]{g, "p4"}));              g = reader.read("LINESTRING(10 120, 40 -50)");         features.add(fb.buildFeature(null, new Object[]{g, "line1"}));              g = reader.read("LINESTRING(1 0, 40 30)");         features.add(fb.buildFeature(null, new Object[]{g, "line1"}));          return fc;     }          private Style createStyle(String geomPropertyName) {         StyleFactory sf = CommonFactoryFinder.getStyleFactory(null);         FilterFactory2 ff = CommonFactoryFinder.getFilterFactory2(null);              Rule pointRule = sf.createRule();         Function getGeomFn = ff.function("geometryType",ff.property(geomPropertyName));         Filter pointFilter = ff.equals(getGeomFn, ff.literal("Point"));         pointRule.setFilter(pointFilter);         pointRule.setElseFilter(false);              Graphic graphic = sf.getDefaultGraphic();         graphic.graphicalSymbols().clear();         graphic.graphicalSymbols().add(sf.getSquareMark());              Symbolizer sym = sf.createPointSymbolizer(graphic, geomPropertyName);         pointRule.symbolizers().add(sym);              Rule lineRule = sf.createRule();         lineRule.setIsElseFilter(true);         Stroke stroke = sf.createStroke(ff.literal(Color.RED),  ff.literal(3.0f));         Symbolizer lineSym = sf.createLineSymbolizer(stroke, geomPropertyName);         lineRule.symbolizers().add(lineSym);              FeatureTypeStyle fts = sf.createFeatureTypeStyle(new Rule[]{pointRule, lineRule});         Style style = sf.createStyle();         style.featureTypeStyles().add(fts);         return style;         }


title: Direction and Distance Point Placement

text: I'm working with 19th century land tax records that provide a direction (expressed as NW, SE, etc.) and a distance in miles from the court house. The info in the Land Tax Records (LTR) is in effect a version of polar radiation. For each reading you start at the court house as your 0,0. The records have: LastName; Direction, Distance or LastName, Distance, Direction. Such as Smith, 4, SW that translates to Smith, 4 miles southwest of the court house. Is there a way to place a symbol/marker/whatever at that point using QGIS 2.0? The end result should roughly result in the shape of the county in 1860. It is possible and very cumbersome to do that with CAD by using the direction and distance alignment on a line and then placing a loci at the end of it. I can use the direction and distance measures by placing them into Excel, plotting them in Vectorworks, and then exporting them to QGIS. I'm trying to cut out the intermediate steps by having QGIS do the work. I can convert the directions to azimuths in Excel and the distances to whatever I need. But having 400 separate lines for a county is cumbersome if there's a more elegant solution.


title: How to place multiple symbols on one polygon centroid?

text: I am using QGIS for visualizing the health conditions of elderly living in an old age home, some of the elderly have multiple medical device with them and I want to show them on the map. Suppose I have one bed for each elderly, I have created a centroid for each bed using `ftool`, then I join this label with my csv table, and then use `Style` tab in the layer properties. In the csv table, I have 4 columns named `A`, `B`, `C`, `D` and each represent the presence/absence of a medical device. Currently I am using `rule-based` in the `Style` tab and use `Offset X,Y` to fine tune the position of the dots so that they don't overlap with each other completely. I don't like the `Offset X,Y` method, because if I need to change the position of those symbols, I need to do it manually one by one, and I have 30 elderly homes that I need to work on. Is there a better way?


title: Leaflet GeoJSON features not showing up in map

text: New to GeoServer and Leaflet...loving it so far, but I am weak in JavaScript. I have just built a simple web page that should load points from a GeoServer WFS service using GeoJson. However, even though the page loads without errors and I can see the GeoJSON returned from GeoServer using Firebug, the points are not showing up on the map. Can anyone look at my code and see what I am missing? Thanks in advance for any assistance!                               <link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.5/leaflet.css" />       <title>title name here</title>


title: How can I GROUP BY linear spatial clusters with PostGIS?

text: I have a series of points in a horizontal line, arranged into discrete spatial clusters by the underlying data pattern. Here's an example of what they look like: ![enter image description here](http://i.stack.imgur.com/wY4ml.jpg) How can I use PostGIS to group them spatially into the highlighted sets A, B and C? This is a harder problem than it initially seemed to me. So far I've tried:   * Using ST_X() to order them, then lag() to identify points between which there are gaps. Unfortunately this just groups the points into "gap" groups and "interior" groups.   * Using convex hulls around arbitrary spatial thresholds. This both failed to separate all clusters and had the side effect of creating mid-cluster separations. (For future reference, this was a dumb idea.)   * Using ST_GeoHash() to cluster them. This doesn't accurately catch the linear nature of the geometry. Any assistance greatly appreciated.


title: QGIS Custom Form: "Other" category in combo box replaced by line edit

text: I am trying to create a custom form for QGIS 2.0 with comboboxes containing an "other" category. I wonder how to provide a line edit to enter info replacing the "other" choice if chosen. Not sure if I can sync a combo box and a line edit at the same time under field properties. In psuedocode:               myfield = dialog.findChild(QComboBox, "Myfieldname")     otherfield = dialog.findChild(QLineEdit, "Myfieldname")          if myfield.currentIndex() == -1:         set myfield to otherfield.text()      This thread's idea seems very similar, but I'm having trouble adapting it: How to automatically populate fields instantly? my custom form code so far:               Naturefield = None     Nat_autrefield = None          def formOpen(dialog, layerid, featureid):         global myDialog         myDialog = dialog         global Naturefield, Nat_autrefield         Naturefield = dialog.findChild(QComboBox, "Nature")         Nat_autrefield = dialog.findChild(QLineEdit, "Nat_autre")         Naturefield.currentIndexChanged.connect(autreChanged)               def autreChanged(value):         #intending to take on form of line edit only if 'Autre' chosen         if Naturefield.currentText() == 'Autre':             Naturefield.setLineEdit(Nat_autrefield)      setLineEdit() currently transforms my combo box to make the list editable for all choices if 'autre' is chosen and the combo box stays a line edit once 'autre' is chosen--not exactly what i'm aiming for. Instead, I'd like the combo box to receive the input from a line edit only if 'autre' is chosen, and revert to normal combo box if other choices are chosen. This is all in the interest of data quality control. I'm new to QtDesigner and Python, so any guidance would be appreciated! Thank you in advance.


title: Cannot Display Data From POSTGIS

text: I try to display the data that I have converted shp into PostGIS database with MapServer and OpenLayers. When I try to show directly shp file is successfully displayed in OpenLayers, but when I try to display the data shp that has been stored in PostGIS can not be displayed in OpenLayers. I had several times looking for this solution in the forums but the result is still an error. Below is the map file that I use.   1. **postgis.config** for connection                CONNECTIONTYPE POSTGIS     CONNECTION "host=localhost dbname=postgisprime user=postgres port=5432 password=admin"     PROCESSING "CLOSE_CONNECTION=DEFER"   1. config.inc.map for projection               CONFIG PROJ_LIB "/ms4w/proj/nad/"  3\. post.map for layer               LAYER           NAME    Base         TYPE    polygon         STATUS  ON         DUMP    TRUE         INCLUDE "postgis.config"         DATA    "geom FROM \"provinces\" USING UNIQUE gid using srid=1"          METADATA         ows_title "Base"         gml_include_items "all"         wfs_include_items "all"         ows_featureid "gid"         ows_enable_request   "*"     END          CLASS                COLOR 242 239 232       END # class      END # layer         END 4.postgisbase.map               MAP       INCLUDE "config.inc.map"     NAME "postgisbase"  #name to give your map service     EXTENT 221238 881125 246486 910582     UNITS meters     CONFIG      "MS_ERRORFILE" "D:/ms4w/tmp/ms_error.txt"          PROJECTION         "init=epsg:4326"     END        WEB         MINSCALEDENOM 100         MAXSCALEDENOM 100000     METADATA             ows_title "PostGIS "             #ows_onlineresource "http://www.postgis.us/demo/chapter_11/GetPAMap.ashx?"             wms_version "1.1.1"             wms_srs   "EPSG:2249 EPSG:4326 EPSG:26986 EPSG:3785 EPSG:900913"             wfs_version "1.0.0"             wfs_srs   "EPSG:4326"         END       END #End Web     INCLUDE "post.map"     END # Map File      and this the script for call the postgisbase.map in openlayers               var Ina = new OpenLayers.Layer.WMS( "Peta Dunia", "../../../cgi-bin/mapserv.exe?MAP=D:/ms4w/apps/AIS/postgisbase.map", {layers:'base'},{wrapDateLine: true},{sphericalMercator: true}); and the error when i try to load the page is               [Wed Oct 26 13:52:32 2011].171000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:32 2011].312000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:32 2011].375000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:32 2011].718000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:32 2011].734000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:32 2011].750000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:32 2011].750000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].0 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].78000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].250000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].328000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].359000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].437000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].515000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].687000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].812000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].906000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].953000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].984000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:33 2011].984000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:34 2011].140000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)     [Wed Oct 26 13:52:34 2011].156000 loadLayer(): Unknown identifier. Parsing error near (provinces):(line 46)      Any suggest for this error ?


title: Extracting data from a raster to a polyline based on start and end points?

text: ![enter image description here](http://i.stack.imgur.com/ZOxBD.png) I have a polyline shapefile with the attributes seen above. startx and starty represent the coordinates of the beginning of each segment of a polyline. endx and endy represent the coordinates of the end of each segment of a polyline. I have a raster with values of Flow Accumulation, known as FA. I would like to extract the values on the raster at the coordinates of the start points and place them in this attribute table at "StartFA" and at the end points and place them in this attribute table at "EndFA". I have so far performed "Feature Vertices to Points" using point types "START" and "END". I then performed "Extract Values to Points". How do I put the values I am looking for into my polyline shapefile attribute table? I am using ArcGIS 10.1 with an ArcInfo license.


title: ArcEngine BaseTool derived classes stops receiving OnKeyDown messages

text: I have two BaseTool derived classes in a custom ArcEngine application. Both tools do the same thing: start a INewBezierCurveFeedback on a first mouse click and then stop it and return the resulting polyline on a second mouse click. Both also have "Overrides Sub OnKeyDown" handlers to allow users to "cancel" the tool by pressing the Escape key. **Test Scenario:** Start with Tool #1, click on the map to start the line feedback and press the escape key and OnKeyDown gets fired, great! Do it again and get same result. **Problem:** Now select Tool #2, and repeat. I don't get the OnKeyDown message here. Go back to Tool #1 and not longer get it there either! I can also reverse this test, starting with Tool #2 and produce the same results. **What's going on here?!? Any help?**


title: how can I add all feature classes of the same name to an mxd in python

text: I have a bunch of folders which all have (almost) identical folder structures and file geodatabases of identical structures which in this folder structure. I need to add all of the same type of feature class to an mxd and merge them into a new feature class. There are a lot of them so I thought this seemed like the sort of thing I could do in python. The folders are all in the same folder and there is a common path to the feature class in question, for example:               project_name/Subfolder1/subfolder2/project_name.gdb/featuredataset/featureclass.      so I'm thinking of trying to loop through the folders and extraxt the feature class from each of them using a method like this, but I'm not quite sure how arcpy interacts with, and loops through windows folders. Since I'm looking for any feature class of this name within the root folder could I use a more brute force method that searches in the root folder and subfolders for the name and adds it to the mxd?


title: How to create an address locator style?

text: Does anyone have any idea about creating or modifying an existing address locator style in order to use it for esri geocoding? I need this to create my own geocoder that respects the specifics of my country, Morocco.


title: How to generate Report with Unicode Character. (Gujarati)

text: I am using ArcMap 10.1 and want to generate report using Report Designer. In Report header i use **Label Control** for display unicode character(Gujarati). but when i copied gujarati text from word doc and assign to label text property, it display as question marks, 'garbled' or 'nonsense'. Is there a simple way to display Gujarati characters in an ArcMap label?


title: Field calculator problem

text: I have an issue withe the field calculator. It is impossible for me to multiplicate 2 columns (area * KZ). The result is always "0". I used the formula area * KZ....is this the right formula? LINK to shape file Here the steps I did:   1. Using the "advanced python field calculator" to replace the strings "A"into 1, "B" into 2, "C" into 3 and so on   2. Using the field calculator (processing plugin) with the formula KZ * area without success ![screenshot field calculator](http://i.stack.imgur.com/FIL5D.jpg)


title: How to get page number of Feature in GeoServer 2.4?

text: I have GeoExt-based client and FeatureGrid with PagingToolbar in it. Also I have a big vector layer (6018 features), so I have to use grid with multiple pages. I want to synchronise grid with map. If I click on feature which is in current grid's page, there's no problem. If the feature is in my database, but not on the current page, I need to go to the page with this feature and choose it in the grid. The problem is I can't get this page. I can get all 6018 features and calculate the page number in JS, but it's bad. So, I should send request to GeoServer to get the page number (with page size specified) for feature with ID specified, but I can't find how to do it. Maybe there is an ability to use something like COUNT the number of features SORTED_BY SMTH, WHERE ID less than MY_ID?


title: PGRouting analyzegraph

text: It shows the message when I try to execute the command `SELECT pgr_analyzeGraph('Highway','the_geom',0.0001);`. I confirmed that the pgrouting, postgis, postgis_topology does including into the database.               function pgr_analyzegraph(unknown, unknown, numeric) does not exist     LINE 1: SELECT pgr_analyzeGraph('Highway','the_geom',0.0001);                    ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.      What is the possibility of this message?


title: how to capture the clicked point before moving a vertex

text: I need the location of clicked point before moving a vertex. IVertexMoved event will provide me the coordinates of the point after vertex is moved but I need before moving. Please help. Regards, Sheethal


title: How to set the cursor to busy mode?

text: How can I set the cursor to busy-mode/in progress from my qgis plugin (cpp)? I couldn't find anything in the API documentation.


title: How do you fix Overlapping Lines?

text: I have a line dataset that contains some lines where the "end points" are actually found in the middle of the line. In the image below, the nodes are shown in green and one of the end points in red. The line doubles back on itself. How would I ensure the end points are actually at the end of the line, and remove the overlap? ![complex line](http://i.stack.imgur.com/hJ2bx.png) These lines are both valid, and complex - as they have overlaps, as confirmed using Shapely.               >>> from shapely.wkt import loads     >>> wkt = "LINESTRING (0 0, 2 0, 1 0)"     >>> l = loads(wkt)     >>> l.is_valid     True     >>> l.is_simple     False      Ordering the coordinates would work in some cases:               >>> sorted(l.coords)     [(0.0, 0.0), (1.0, 0.0), (2.0, 0.0)]      But not in others as it would change the geometry of the line:               >>> wkt = "LINESTRING (0 0, 2 0, 1 -1)"     >>> l = loads(wkt)     >>> sorted(l.coords)     [(0.0, 0.0), (1.0, -1.0), (2.0, 0.0)]      Attempts at buffering the line with a 0 buffer, and merging lines have not changed the coordinate order.


title: Contours: Interrupt contour lines beneath a elevation number

text: Is there a way to interrupt a contour line beneath a elevation number? ![enter image description here](http://i.stack.imgur.com/p1RCy.png) Thanks!


title: Is it permissible to redefine raster resolution (not resample!) in order that it be included in a SDM?

text: I'm fitting a SDM using R. A requirement is that all environmental predictor layers have the same extent, origin and resolution. I have high resolution bathymetry data (25m) that I will resample to 100m (precision of occurrence records). However, some environmental predictors pulled from a hydrodynamic model are at a much lower resolution (~7km) e.g. temperature / salinity. I am aware of the introduction of error in resampling from low -> high resolutions. However, resampling the from 25m -> 7km is not suitable for my requirements. Is it permissible, using QGIS, to define the resolution of the low resolution layers by saving the resolution to that of the higher res layers? Please note I am not resampling - just changing the number of columns & rows in the saved tiff so that R sees the layer as having an equal resolution to the other layers.


title: How can I convert SRTEXT to PROJ4TEXT with QGIS

text: I have a SRTEXT like this: `GEOGCS["My Data",DATUM["Not specified (.....` I want to get PROJ4TEXT of this STREXT. Can I get it from QGIS. Because I want to insert it in Postgresql SPATIAL_REF_SYS table. (or I need PROJ4TEXT data to use this projection)


title: Run the Sub procedure by selecting an object from the map

text: In Mapbasic made custom menu selected item to run sub procedure, what i want to make the Sub procedure to be ready to run by selecting an object from the map and run again by selecting another object, and so on


title: Move and zoom map-layer while editing

text: I want to create lines with the help of "CAD-like improved line capture", make a point, pretend an angle or distance want to put the second point. To this end, I have to zoom out on the map, select a different part of the map and there zoom in again. Here, however, I lost the tool and the firt point/vertex ist lost. Is there a way to prevent this? thank you very much * * * original German: Ich möchte mit Hilfe von "CAD-like improved line capture" Linien erstellen, erstelle einen Punkt, geben einen Winkel oder Distanz vor, möchte den zweiten Punkt setzen. Dazu muss ich aber aus der Karte zoomen, einen anderen Auschnitt wählen und dort wieder hineinzommen. Dabei allerdings habe ich das Werkzeug wohl "abgewählt" und der Punkt ist verloren. Gibt es eine Möglichkeit dies zu verhindern? herzlichen Dank


title: ArcGIS Android - Get LatLon position of the graphic

text: I'm adding marker to my map using the following code:               mvHelper.addMarkerGraphic(Double.parseDouble(friendData.getLatitude()), Double.parseDouble(friendData.getLongitude()), friendData.getFriendName(), friendData.getStatus(), null, icon, false,                     0);      This is my listener:               mvHelper.setOnGraphicClickListener(new OnGraphicClickListener() {             @Override             public void onGraphicClick(Graphic graphic) {                 Log.d("MVHELPER GRAPHIC LISTENER TAG",graphic.toString());                 calculateDirection((float)((Point)graphic.getGeometry()).getX(),(float)((Point)graphic.getGeometry()).getY());             }         });        I'm trying to get the lat long position of the clicked marker. But the listener doesn't provide it. So I tried the above way and got this error: `com.esri.core.io.EsriServiceException:Unable to complete operation. Location "Location 2" in "Stops" is unlocated. Need atleast 2 valid stops. "Stops" does not contain valid input for any route.` and this is the method where I get the route:               calculateDirection(float x, float y){     Log.d("Marker X and Y",String.format("%f // %f",x,y));     final Point loc = mapView.toMapPoint(x, y);     Log.d("Marker X and Y in MapPoint",String.format("%f // %f",loc.getX(),loc.getY()));     Thread t = new Thread() {             @Override             public void run() {                 try {                          RouteParameters rp = mRouteTask                             .retrieveDefaultRouteTaskParameters();                     NAFeaturesAsFeature rfaf = new NAFeaturesAsFeature();                          Point p = (Point) GeometryEngine.project(loc, wm,                             egs);                     Log.d("Marker X and Y after projection",String.format("%f // %f",p.getX(),p.getY()));                     StopGraphic point1 = new StopGraphic(mLocation);                     StopGraphic point2 = new StopGraphic(p);                     rfaf.setFeatures(new Graphic[] { point1, point2 });                     rfaf.setCompressedRequest(true);                     rp.setStops(rfaf);                          rp.setOutSpatialReference(wm);                          mResults = mRouteTask.solve(rp);                     mHandler.post(mUpdateResults);                 } catch (Exception e) {                     mException = e;                     mHandler.post(mUpdateResults);                 }             }         };              t.start();      } **Update 1:** After trying @jdONeill answer. I get the Lat and Long values.This is what I did.               Point temp = (Point)graphic.getGeometry();          Point p = (Point) GeometryEngine.project(temp, wm,egs);     calculateDirection((float) p.getX(), (float) p.getY());      This is the result:               04-30 16:17:54.408  17908-17908/com.test.esrisample D/Marker X and Y﹕ -93.203285 // 44.881104     04-30 16:17:54.408  17908-17908/com.test.esrisample D/Marker X and Y in MapPoint﹕ -10379311.931913 // 5607101.545301     04-30 16:17:54.458  17908-17908/com.test.esrisample D/dalvikvm﹕ GC_FOR_ALLOC freed 484K, 22% free 10018K/12760K, paused 16ms, total 16ms     04-30 16:17:54.719  17908-18235/com.test.esrisample D/Marker X and Y after projection﹕ -93.238945 // 44.908331      Why is there a slight change in the latitude and longitude values?


title: Can I create a polyline with one long segment and three shorter underlying segments?

text: We currently have a street centerline polyline feature class where the street segments correspond to our pavement management database. Each segment has a unique ID that matches a record tracking pavement type/year. Segment breaks occur at intersections, but can also occur mid-block if a subdivision was built in multiple additions. We would like to start adding Left From/To and Right From/To fields to the street centerline feature class for geocoding purposes. I would like to be able to have a street segment start at an intersection and end at an intersection for address coding, but still keep the underlying multiple segments that correspond to the breaks in the pavement management database. Is there an easy way to do this? Do I need to have two seperate line files or is there way to have this in one feature class?


title: Calculate distance between X/Y coordinates by individual in R

text: In R, I'd like to calculate the straight-line distance between each capture of an animal in my dataset - no need to take into account the curvature of the earth. I have a dataset that looks like this:               record.nbr<- c(1:8)     individual<-  c(1,2,1,2,1,3,1,2)     x<-c(167685,167945,167685,153985,167685,158675,167645,167667)     y<-c(9876548,9879248,9876838,9596548,9926548,9878578,9876548,9166548)     julian.date<-c(125,126,127,127,128,129,130,130)     captures<-data.frame(cbind(record.nbr, individual, x, y, julian.date))      Here each "individual" represents a unique animal, and x/y points represent locations where the animals were captured. x/y points are UTM coordinates. Some individuals are caught many times, while other individuals are only captured once. I'd like to calculate the distance (and the mean date) between recaptures by individual. In the output table, each row would summarize the distance between consecutive recaptures of an individual, plus the mean date between recaptures, like so:               individual first.rec  second.rec  distance  mean.date     1          1             3        V         126     1          3             5        W         127.5     1          5             7        X         129     2          2             4        Y         126.5     3          4             6        Z         128


title: How to use Scale tool when Georeferencing Raster in ArcMap?

text: I'm georeferencing a bunch of maps in ArcMap 10.0. The Shift and Rotate commands are pretty easy to use, no problems there. But whenever I use the Scale command, the image completely changes position. Is there a way to put an anchor down (similar to the rotate tool)? I feel like I must be missing something but haven't figured out what it is yet. I don't want to rubbersheet them into place, so I'm trying not to use control points until I have it in position, but first I need to make sure the scale bar on the map matches the distance it's supposed to. Edit: Asking a question does wonders in finding an answer, even when nobody has answered yet. I found the Rescale_management tool and will use that to set the scale properly first, then position. If the length of the scale bar is 1256ft and it should be 1200ft, I'll use a factor of .9554 in the X and the Y directions. Is there a better way?


title: How to assign node IDs to links in a network?

text: I have two shapefile layers in ArcGIS for Desktop 10. I have one line layer and node layer. My node layer(shp) has a "code_1" field and line layer has a "code_node" field. The intersection point of the line layer has only one node. I want get the code from the "code_1" field of the node layer for the line's "code_node" field For example, for one line, the first node has code_1=1 and the node at the end of line has code_1=2. In this case the `code_node` attribute in that line must be "12",. How can I automatically get the code from the nodes for each line? ![enter image description here](http://i.stack.imgur.com/vVegb.gif)


title: Can i access UTF grids without a server?

text: I'm creating a slippery map that will run off of a usb stick or cd (client request). I'm able to load up the map, display pre-cooked tiles and some some simple vectors. The client would like some more interactivity and the level of detail they want (cell values from rasters) is far to large for json files. I've made UTF grids using tilestache and have created a directory for them in the same folder as the map html. However, the UTF grids wont display and i suspect it is because it violates the same-origin policy, so my questions are as follows:   1. Can i access and dispaly UTF grids in local storage without a server?   2. Is there another route i should be taking other than UTF grids? Thank you


title: What is the difference between Sub and a Function in Mapbasic?

text: In MapBasic it seems like Functions can have side effects so I don't really see the difference between using the following: `Declare Sub MyProcedure (myParameter As String)` and then use `myParameter` after `Call MyProcedure(myParameter)` or `Declare Function MyProcedure () As string` and do something like `myParameter = MyProcedure()` Am I missing something?


I am trying to set up Geoportal to harvest services hosted by Geoserver. I tried follow this instructions:   1. Open your Geoportal Server.    2. Go to the Administration tab Choose the 'Add' option    3. Choose 'Register resource on the network' and click Proceed    4. In the URL resource form, paste the OGC GetCapabilities URL (including all parameters)    5. Keep the options for allowing the resource to be found and for it to be synchronized at an interval    6. Choose the interval.    7. Select if you want to be notified of synchronization results.    8. Click Save.    9. In the Manage page, approve this newly registered resource. Now I did the following in Geoserver:   1. Registered a ArcSDE shapefile in Geoserver that I can see in "Layer Preview" and I can get the common formats like OpenStreet map working for the layer. But what is an OGC GetCapabilities URL??? Where do i declared and link it to the layer I registered? Many thanks for any help Carlos.


title: Moving from MapInfo to web mapping

text: I use MapInfo and have a lot of data in .tab format that I would like to put into a database in order to access it through a website. I would like to continue using MapInfo as my GIS programme and export to the web database with the minimum of fuss. What would be a good type of database for this? I was considering MySQL as I have used it a lot before and it has native spatial data types. However I have read that MapInfo can't connect directly to MySQL and that its spatial data support is a bit patchy. Alot of discussion on the subject is several years old though, so can anyone confirm or deny this? If not MySQL, what is a good choice? PostgreSQL/PostGIS?


title: Is ANZLIC Metadata supported using the new approach of ArcGIS 10.x?

text: I like the new approach to metadata at ArcGIS 10 because I fall into the camp that prefers minimal metadata for many purposes, but like to have the option to work with more comprehensive standards when they are appropriate. However, I'm keen to hear whether people think ArcGIS 10 works well for those trying to adhere to ANZLIC Metadata Profile 1.1? My understanding is that the ANZLIC Metadata Profile 1.1 (based on ISO 19115) is very similar to the ISO 19139 Metadata Implementation Specification so can the ArcGIS 10 metadata style called ISO 19139 Metadata Implementation Specification simply be used? Or is the ArcGIS/ANZLIC community awaiting the release of an ANZLIC specific metadata style for ArcGIS 10 by Esri or a third party? A reply to an ArcGIS Discussion Forum thread suggests that it might be.


title: Maps, black and white/ gray scale printing

text: I'm working on Mapserver to do beautiful maps with AGG rendering Until now, I've only thought about the best way to get color maps. Now, I have to think how to make color ones compatible with a grey scale printing http://en.wikipedia.org/wiki/Grayscale How do you deal with this? Any advices, tips or resources will be welcome. * * * Edit : My question is really more relative to visual perception than Mapserver. It seems I will have to read Semiology of Graphics by Jacques Bertin (reedited this year http://www.esri.com/news/releases/10_4qtr/bertin.html ) or look on Edward Tufte various book/articles http://www.edwardtufte.com/tufte/


title: Cutting up a shape file into smaller parts

text: I've downloaded the 350MB coastline shape file from Openstreetmap, but I really would like to work with a small part of them, is there a way to generate a shapefile of a small area. I'm using MacOSX or Unix for this. Edit: Apparently the coastline file from Openstreetmap wasn't one big polygon so everything worked fine. Without cutting polygons.


title: Updating a layer's bounding box in GeoServer after WFS-T insert

text: I'm creating a data store in GeoServer from an OpenLayers DrawFeature control, using WFS-T to save the geometry into GeoServer. This works great, except the data store's bounding box doesn't get set in GeoServer. However, if you create a data store from a shapefile, the bounding box gets automatically set. I know you can go into the GeoServer admin page and tell it to calculate the bounding box, but I need a way to do it programmatically.


title: Enable FIND in ArcGIS Server 10.1

text: Trying to get a FindTask working via this sample and the only difference between my code and the samples that I can see is that the sample data has a find operation enabled. I can't find in my server manager where to enable this. And is it even necessary? Right now the search isn't returning any results. ![ESRI Service](http://i.stack.imgur.com/9amAl.png)               <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN" "http://www.w3.org/TR/html4/strict.dtd">      Find Task               <style type="text/css">       html, body { height: 100%; width: 100%; margin: 0; padding: 0; overflow:hidden; }       body { font-size: 0.9em; font-family: Geneva, Arial, Helvetica, sans-serif;       } .details { font-weight: bold; } #grid { border:       1px solid #333;}     </style>     <script type="text/javascript">       dojoConfig = {         parseOnLoad: true       };     </script>     <script type="text/javascript" src="http://serverapi.arcgisonline.com/jsapi/arcgis/?v=3.2"></script>     <script type="text/javascript" language="Javascript">       dojo.require("dojox.grid.DataGrid");       dojo.require("dojo.data.ItemFileReadStore");       dojo.require("esri.map");       dojo.require("esri.tasks.find");       dojo.require("dijit.layout.BorderContainer");       dojo.require("dijit.layout.ContentPane");            var findTask, findParams, map;            function init() {         var initExtent = new esri.geometry.Extent({"xmin":-9199621.530456403,                 "ymin":3365235.5214724857,                 "xmax":-9077475.159256855,                 "ymax":3443965.6606061114,                 "spatialReference":{"wkid":102100}});         map = new esri.Map("map",{extent:initExtent});         var censusMapLayer = new esri.layers.ArcGISDynamicMapServiceLayer("http://gisprod2/arcgis/rest/services/Dynamic/Zoning/MapServer", {         });         map.addLayer(censusMapLayer);              //create find task with url to map service         findTask = new esri.tasks.FindTask("http://gisprod2/arcgis/rest/services/Dynamic/Zoning/MapServer");              //create find parameters and define known values         findParams = new esri.tasks.FindParameters();         findParams.returnGeometry = true;         findParams.layerIds = [4];         findParams.searchFields = ["ZONECLASS"];              dojo.connect(map, 'onLoad', function(theMap) {           //resize the map when the browser resizes           dojo.connect(dijit.byId('map'), 'resize', map,map.resize);         });       }                 function execute(searchText) {         //set the search text to find parameters         findParams.searchText = searchText;         findTask.execute(findParams, showResults);       }            function showResults(results) {         //symbology for graphics         var markerSymbol = new esri.symbol.SimpleMarkerSymbol(esri.symbol.SimpleMarkerSymbol.STYLE_SQUARE, 10, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new dojo.Color([255, 0, 0]), 1), new dojo.Color([0, 255, 0, 0.25]));         var lineSymbol = new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_DASH, new dojo.Color([255, 0, 0]), 1);         var polygonSymbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_NONE, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_DASHDOT, new dojo.Color([255, 0, 0]), 2), new dojo.Color([255, 255, 0, 0.25]));              //find results return an array of findResult.         map.graphics.clear();         var dataForGrid = [];         //Build an array of attribute information and add each found graphic to the map         dojo.forEach(results, function(result) {           var graphic = result.feature;           dataForGrid.push([result.layerName, result.foundFieldName, result.value]);           switch (graphic.geometry.type) {           case "point":             graphic.setSymbol(markerSymbol);             break;           case "polyline":             graphic.setSymbol(lineSymbol);             break;           case "polygon":             graphic.setSymbol(polygonSymbol);             break;           }           map.graphics.add(graphic);         });         var data = {           items: dataForGrid         };         var store = new dojo.data.ItemFileReadStore({           data: data         });         grid.setStore(store);       }       dojo.addOnLoad(init);     </script>      Find State/City/River: Layer Name Field Name Value


title: How to publish my map from geoserver

text: > **Possible Duplicate:**   >  How to develop a map with layers in Geoserver? I'm trying to create a simple map for my city. I have a shape file my city, but I can't create on the web. I tried to take from geoserver to use geojson format but I didn't. Can you help me please.


title: How to Display Live Incoming Smartphone Images with Ushahidi on a Map

text: I have not yet used Ushahidi and I could not find the answer on their website to my question. Can I easily do the following?   1. setup an Ushahidi server/map,    2. go outside to make photos with my Android/iPhone smartphone   3. upload these photos to Ushahidi which should be saved in a database   4. display the photos on the map according to the coordinates. Also few lines of comments would be nice. When I now visit my site and open the map, I would like to see some dots (actual position of an entry/uploaded image). By double clicking on the dot, the original uploaded image+comments should be displayed(pop up). Is this possible with Ushahidi? edit: nobody working with Ushahidi?


title: British National Grids projection does not add Ordinance Survey grids to shapefile

text: I have created a shapefile in ArcCatalog by right-clicking the folder I want it in, then selecting new, then shapefile, then changing the type to "polygon", and using the British National Grids under the projected co- ordinate system. I then tried to add the shapefile in ArcMap by clicking the Add Data button and selecting the shapefile I created. My problem is that when I open the shapefile in ArcMap, the map is still blank and I cannot see the projected shape of the British National Grids I need. Please help.


title: How to add more events like Control.Click and Handler.Point for a single variable in OpenLayers

text: I have two points which are marked dynamically on the map and when I mark the start point it is returning lon,lat, but I want to assign this coordinates to the variable that marked/handled the start symbol and a point. Ex:               measureControls={                   start: new OpenLayers.Control.DrawFeature(start, OpenLayers.Handler.Point),                   stop: new OpenLayers.Control.DrawFeature(stop, OpenLayers.Handler.Point),      for the same start and stop, I would like to add event 'click' and assign the point when clicked. Thank you


title: Is it possible to use symbols instead of colours to visualise raster cell values?

text: I have a raster image (GeoTIFF) with _integer_ values. I know how to visualise raster values in QGIS using colormaps or similar, but I would like to assign symbols to each colour value instead - just like a colormap, but with symbols. **Example:** ![enter image description here](http://i.stack.imgur.com/IRWBO.png) **Question:**   Is this possible and in that case how?


title: Only return the NumberOfFeatures in a WFS query

text: I'm running a GetFeature WFS query in Geoserver (looking for features within a polygon) but all I actually want is the total number of features. Is there a way of just returning this info (in the hope that the query will run a little faster) rather than all the records? The query is returning over 10,000 features so I hope it would make a difference if possible. Many thanks, Steve


title: QGIS project on QGIS server -- cannot load wms layer

text: I've created a qgis project on QGIS 1.8, copied it onto the QGIS Server, and am trying to load the project as a WMS Layer back into QGIS. There is one vector layer in the project. Paths have been set to absolute and the CRS for both the layer and project is EPSG 27700. In the "Add Layer(s) from a Server" dialog box, the project is showing, but not the vector layer. Clicking on the project, at the bottom of the dialog box it says "No Common CRS for selected layers". Any help is very welcome -- I have searched the forums and tried the solutions suggested QGIS Server Capabilities Response doesn't contain any Layers and QGIS Server Capabilities Response doesn't contain any Layers (Fedora 15 64-bit) but still cannot get it to work. Many thanks!


title: Greyscale Raster (tif) in ArcMap

text: I've searched all over google and have found nothing. I want some tiff files I have in ArcMap to be in greyscale with the rest of the layers in colour, similar to the style override > greyscale option in MapInfo. I will be printing A1 posters and do not want the background imagery in colour. I've tried changing the symbology to stretched (as suggested on a forum), but it looks weird, like this: ![enter image description here](http://i.stack.imgur.com/PuzKv.png)


title: Cropping elements from one layer to an element from another

text: I'm trying to make a hex map of a country (Russia, in this case). To do this, I've opened a map of Russia in QGIS (from Natural Earth), set up my projection, and created a hex grid with MMQGIS. So far, so good. Now I want to crop the hexes with the shape of Russia. I've tried using the intersect option, but I get this error:               "Input CRS error: One or more input layers missing coordinate      reference information, results may not be as expected."      and an empty shapefile for my troubles. What am I doing wrong? UPDATE: Ok, so now I have them both in the same coordinate system (based on an Albers projection using meters), but I'm still not able to get an intersection to work: > Warnings: Input CRS error: One or more input layers missing coordinate > reference information, results may not be as expected. GEOS geoprocessing > error: One or more input features have invalid geometry. Some output > geometries may be missing or invalid. The output shapefile is empty. I _was_ able to generate a version that almost works using WGS84 coordinates, but that had some problems with the 180° meridian... ![](http://f.cl.ly/items/3E0p2V3Y3d1F2X0C3T24/QGIS.png)


title: automatic templates KML

text: I have imported several points from QGIS into Google Earth. When I click over an icon point, I get a ballon with all the attributes, OK. One of the attributes is called Text. Ok, so what if i want to take the Text value to be used in ot modify the balloon? I don't want to have to manually write, point by point, the value. i have seen in this URL that you can use extendedData to create a template. The KML generated by QGIS, instead of this example snippet (used in Google tutorials)               <Placemark>     <styleUrl>#golf-balloon-style</styleUrl>     <ExtendedData>       <Data name="holeNumber">         <value>1</value>       </Data>     </ExtendedData>      QGIS creates this...               <Placemark>     <styleUrl>#golf-balloon-style</styleUrl>     <ExtendedData>     <SchemaData schemaUrl="#punts_R">       <SimpleData name="Text">my text</SimpleData>         </SchemaData>     </ExtendedData>      I wonder if this QGIS created schema there is some way the parameters replacement may work (I didn't get it work, but yes with the Google 'standard'). It would save me a lot of work...


title: Why isn't StackBuilder installation of PostGIS working?

text: I've installed PostGIS (ver. 2.0.1) on my Mac using EnterpriseDB's StackBuilder, which I also used for my PostgreSQL (9.2.2) installation. Although StackBuilder says that PostGIS is now installed, and although I have a PostGIS directory in /usr/local/PostgreSQL/9.2/, it doesn't appear to be actually installed (see queries below).               gregory=# select postgis_full_version();     ERROR:  function postgis_full_version() does not exist     LINE 1: select postgis_full_version();                    ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.          gregory=# create extension postgis;     ERROR:  could not open extension control file "/usr/local/PostgreSQL/9.2/share/postgresql/extension/postgis.control": No such file or directory      and `select * from pg_available_extensions;` does not return a record for PostGIS. Does anyone know what's up with this and how to get it to work?


title: Unable to display FeatureLayer using ArcGIS server

text: I can't figure out why this FeatureLayer doesn't show up on top of my base map. I'm using ArcGIS Server 10.0 and the ArcGIS JS API 2.6. Here's the information from the service's REST parameters page:               Feature Layer: Inspection Locations (ID: 5)     Display Field: NAME     Type: Feature Layer     Geometry Type: esriGeometryPoint          Min. Scale: 250000     Max. Scale: 0          Extent:      XMin: -8879044.1973     YMin: 4938432.3627     XMax: -7999181.0204     YMax: 5622923.0498     Spatial Reference: 102100 (3857)     .     .     .     Supported Interfaces:   REST          Supported Operations:   Query   Add Features   Update Features   Delete Features   Apply      Edits ...and here's the code:                 <script type="text/javascript"> var djConfig = { parseOnLoad: true }; </script>       <script type="text/javascript" src="http://serverapi.arcgisonline.com/jsapi/arcgis/?v=2.6"></script>       <script type="text/javascript">djConfig = { parseOnLoad:true };</script>            <script type="text/javascript">            dojo.require("esri.map");       dojo.require("dijit.layout.ContentPane");       dojo.require("esri.layers.FeatureLayer");       var map;            function init() {              var initExtent = new esri.geometry.Extent({           "xmin": -8879044.1973,           "ymin": 4938432.3627,           "xmax": -7999181.0204,           "ymax": 5622923.0498,           "spatialReference": {             "wkid": 102100           }              });         map = new esri.Map("map", {extent: initExtent });              var basemap = new esri.layers.ArcGISTiledMapServiceLayer("http://server.arcgisonline.com/ArcGIS/rest/services/World_Topo_Map/MapServer");              var featureLayer = new esri.layers.FeatureLayer("http://gis4/ArcGIS/rest/services/NY/NYeFind/FeatureServer/5",{          mode: esri.layers.FeatureLayer.MODE_SNAPSHOT,          outFields: ["*"]        });              map.addLayers([basemap, featureLayer]);              var resizeTimer;         dojo.connect(map, 'onLoad', function(theMap) {           dojo.connect(dijit.byId('map'), 'resize', function() {             clearTimeout(resizeTimer);             resizeTimer = setTimeout(function() {               map.resize();               map.reposition();             }, 500);           });         });              } // end init              dojo.addOnLoad(init);      Oddly enough, if I hit the "View in ArcGIS.com Map" link, it works perfectly.


title: How to increase RAM allocation to QGIS?

text: I am using QGIS 1.8.0 and it runs slowly (with huge layers the size of NSW, Australia) and it appears that QGIS is limited in its RAM allocation to 128.9MB. How can I increase the speed with which it deals with such big files by allocating a larger amount of my computers power to QGIS? I am running Windows 8 Pro (64-bit) on: \- Intel Pentium P6200 @ 2.13GHz (dual core); \- 2.00GB RAM \- AMD Radeon HD 6300M Series   * Engine clock speed: 500-750 MHz   * Processing power (single precision): 80-120 GigaFLOPS   * Polygon throughput: 125-187.5M polygons/sec   * Data fetch rate (32-bit): 16-24 billion fetches/sec Texel fill rate (bilinear filtered): 4-6 Gigatexels/sec Pixel fill rate: 2-3   * Gigapixels/sec Anti-aliased pixel fill rate: 8-12 Gigasamples/sec   * Memory clock speed: 800-900 MHz Memory data rate: 1.6-1.8 Gbps Memory bandwidth: 12.8-14.4 GB/sec


title: count lines in polygons

text: I found the following link to a post saying that a count lines in polygons had been added to ftools for QGIS: http://hub.qgis.org/issues/4712 However I can't find how to access this feature nor any online discussion of how to do this - can anyone help? Alternatively, are there any other tools that will allow this? I have a hexagonal grid and want to count the number of animal track paths that go through each cell for a migratory distribution map. Thanks


title: Switching career from GIS Technician/Analyst to GIS Web Map Developer/Programmer?

text: I am a GIS Analyst at work that dabbles in python scripting in ArcGIS and the Google Maps JavaScript API V3. My professional goal is to either be in web mapping development or desktop/server GIS application development. I have seem to hit a plateau and I am not sure how to progress with my learnings or professionally to reach my goal. I ask those of you out there that made the jump from technician/analyst how did you do it? How can I land a job in web mapping/application development with little to no experience.


title: Binding IActiveViewEvents::ItemReordered to method in my ATL class

text: May you help me please to know how I can wire this event to my method? I have an ATL class and I need to know when a layer added to map object so that I can update my list of layers in this class.My class implements IExtension and I've put this code in Startup method of this interface: //Some code to obtain a reference to Map object __hook(&IActiveViewEvents::ItemReordered,pMap/ _pointer to map_ /,&MyClass::MyEventHandler); and in Shutdown method of IExtension I unhook it,but I receive compile error says the source and handler method not compatible while both signatures,arguments and names are exactly same. Is it true a all this method?any simpler method for this wiring?or any reference on net? Regards


title: How to interpolate daily climate data using ArcGIS for Desktop?

text: I want to interpolate daily climate data available from 3 locations. I am not sure how to interpolate the data for daily values for 30 years from the available 3 locations.


title: How to display coordinates of a point in a projected cooridante system?

text: I want to check the coordinates of a point. On selecting this point in ArcGIS it displays a latitude and Longitude . In ArcGIS I found an option in 'Identity Tool' that the Latitude and Longitude can be changed into desired units say for eg: Decimal degrees, decimeters etc. But when the same file is opened in QGIS and the same point is selected it shows me a Latitude and Longitude point but its unit can not be adjusted. In QGIS, is there a method to change the display units of the point, similar to how ArcGIS does this?


title: No difference between distances measured in projected and non-projected data

text: I'm performing measurements of distances in Km between a lot of localities by using the "Point Distance" in ArcGis 9.3. I realized that the results obtained after using projected coordinate system are the same of those obtained without projecting data. I change projected coordinate system in "Data frame properties". Where am I wrong? I am a newbe in Gis and in ArcGis too. Any advices will be invaluable. Thanks in advance.


title: Route optimisation for multiple vehicles

text: I have 100 destinations and 5 vehicles and I need to code a solution that efficiently routes the vehicles to each destination, so each destination is visited by 1 of the vehicles. Some of these destinations could have time windows to be visited during as well. I'm using PostgreSQL and PostGIS, and was looking to use pgrouting but I'm not sure it's suitable for multiple vehicles - from my limited knowledge of Dijkstra etc I think they're designed for a single vehicle. Any thoughts on whether pgrouting can solve this problem, and if so, any code examples? If not, are there any open-source alternatives that can do it?


title: How to sensibly digitise this raster, scriptably, in ArcGIS

text: I have the raster below which I want to convert to a set of vector lines: ![Raster Lines](http://i.stack.imgur.com/m1b8k.png) As you can see, there are small breaks in the lines and some bits are double- width. I want to be able to reliably import these sort of images to a vector in ArcGIS. I've spent a lot of time trying to get tools like Integrate, Extend Lines, Trim Lines and Unsplit Lines to do what I want, but I just can't seem to achieve it. I guess this problem must have been experienced before, does anyone have any ideas for a workflow that would do what I want? I thought of maybe trying to clean up the raster before-hand, but I'm not sure how to go about that. Again, any ideas would be much appreciated.


title: Maximum number of Tiled Layers in Open Layers

text: I'm developing an OpenLayers site where the client wants a dozen layers that can be individually turned on and off, which is a bit of a bummer in terms of performance. The data is PostGIS and static enough to cache tiles. The area is quite small geographically but requires 9 zoom levels to allow detailed views. I'm using a cached tiled basemap which is fast and efficient. I've been using Single Tile WMS layers for the other dozen layers, but am wondering if I should used Tiled, cached versions of these. There are other efficiency gains to be made by limiting the buffer, and adding sub-domains such that images can come from tiles.example.com, tiles2.example.com, etc, however all other efficiency strategies aside, I'm wondering at what point the OpenLayers performance difference is for 12 untiled layers, verses 12 tiled layers? Said another way, at what number of tiled layers in OpenLayers will the performance degrade such that Single Tile layers are preferable? If I have 12 tiles per view, and 12 layers, can OpenLayers handle 144 tiles without bogging down?


title: Creating fixed width polygons in ArcMap

text: I've been adding easements, which typically have a fixed width. Are there any existing tools or scripts that would allow me to define one edge and automatically create the other? The buffer tool is poor, the line offset tool is for lines. I want to be able to use the same tools (sketch, intersect, endpoint arc) and have actual curves,not point approximations (what buffer does).


title: Where does "Update measurements" put its updated measurements? (MeemsTools)

text: I have several .shp files loaded into MapWindows GIS. I go to plug-ins > Meems Tools and click on "UpdateMeasurement Plugin", and it loads. I click on the Update measurements, and in the the Update measurements window that appears, I click the "Update measurements" button, without changing any of the options in the window. The "select attribute" section grays out, and the bottom of the window lists "Finished updating measurement values", next to the "MapWindow Wiki help" link. > Where does it put these updated measurement values?    possibly relevant information: I'm using Windows 7 version 6.1 .   The files are loaded as multiple layers in the same window of MapWindow GIS.   The files are simplifications using MapShaper.org of the "district" files from "Boundary" on http://hawaii.wr.usgs.gov/maui/data.html .


title: Multi-Colored Polygon in Thematic Map

text: I have to make a thematic map with an attribute table with duplicated features but different attributes. See below:               FEATURE   ATTRIBUTE     -------------------     1         bl     1         gr     2         gr     2         re     3         re     4         bl     5         ye     5         ye      Does anyone knows how to paint a feature (polygon) with the corresponding colors? I mean like for feature 1 paint half blue and half green, for feature 2 half green half red... Thanks.


title: Need gps trackers for social science research

text: I'm interested in studying people's movement over time (24 hours), and would like to do so with devices that also measure elevation very accurately. I understand there's stuff out the where I can buy these for about $20 a piece, as I need about 30 of them. The data will be used in arcgis. Can someone recommend a specific brand to buy?


title: Bing Map Aerial Question

text: I am aware of Microsoft Bing Maps Aerial licensing changes with ArcGIS for Desktop at the end of this year. My question here is why I have two different ArcMap I saved had them in it. One was able to show it and the other would not allowed me to show me? What is the reason for that? ArcGIS 10.1 SP 1 Here is screenshot: ![Here is the Sample of the screenshot](http://i.stack.imgur.com/oZGlQ.jpg)


title: How to perform Shapefile to Raster with ArcMap 10.1 while preserving resolution and values?

text: I am new to ArcMap and I would really appreciate some help. I want to create a new raster using a value in a column of a point shapefile. I have a point shapefile with points equally spaced and each placed at the center of each pixel of a raster layer. (I created this shapefile using 'Raster to point' tool) Now I want to create a raster with the same resolution of the original raster but using as value the values of column in the shapefile table. How can I do that? I tried running 'Point to raster'(black squares)but it's not doing what I am looking for. ![enter image description here](http://i.stack.imgur.com/qaPra.png) This is what happens running 'point to raster' setting the cell size as the original raster: ![cell size set equal to original raster](http://i.stack.imgur.com/Agr2q.png)


title: ArcSDE installation; "sdesetup: command not found"

text: I am trying to install ArcSDE component for postgresql database from A DVD file. I have copied files of following folders on the server :   1. `command_64`   2. `pg_64` in both of theses files there are `install` files. for both I used :               ./install -load      yes     mount point: DVD location or wherever you are installing from     install directory: enter to accept /home/sde     Package numbers to load: enter to accept all     yes     no      now I have sdeexe101 file which contains bin, etc, liband locale folders in it. I also As the root user, copied the st_geometry.so, libsg.so, and libpe.so files from $SDEHOME to the lib directory in the postgres location.               [root@machine pgsql]# cp /home/sde/sdeexe101/lib/st_geometry.so /usr/lib/pgsql     [root@machine pgsql]# cp /home/sde/sdeexe101/lib/libsg.so /usr/lib/pgsql     [root@machine pgsql]# cp /home/sde/sdeexe101/lib/libpe.so /usr/lib/pgsql      when I tried to execute sdesetup to create geodatabase and authorize it using:               sdesetup -o install -d POSTGRESQL -u sde -p sde...      I got this: **`sdesetup: command not found`** what am I missing here? Are pg_64 and commands_64 files that I needed? I appreciate any help!


title: GeoServer and GeoExplorer issues setting up per-layer security

text: I can't get per-layer security working... I have the following set up:   1. A user USER, with no group membership and one role EXAMPLE_ROLE and a password   2. the role has no parent role and no properties   3. under data security I have:      * _._.r set to *      * workspace.mylayer.w set to EXAMPLE_ROLE   4. All the services security is set to _._ * But, when USER logs in to geoexplorer's saved map with mylayer in there, USER can edit, but saving returns the error message "{server}mylayer is read-only. If I set the second secirity to workspace.mylayer.w * then that user can edit. What's going on? What am I missing? EDIT: using version 3.1 of OpenGeoSuite, GeoServer 2.2-SNAPSHOT EDIT2: on looking into it further, I noticed that I don't get the login button on geoexplorer. I also get this error on Chrome's console:               Failed to load resource: the server responded with a status of 404 (Not Found) http://tobler-v:8080/geoserver/rest/imports      and I wonder if this is something related to this: GeoExplorer provokes "401 Unauthorized" request and lacks login


title: How do I get an accurate time attribute when importing a GPX layer as a vector

text: When I import a GPX file in QGIS (as a vector as track points) and I check the time attribute it appears to be 4 hours later than it should be. Interestingly when I take this same GPX file and use GPSBabel to convert it to a CSV file I find that the time listed in the CSV is the correct time. So am I doing something wrong or is this a bug in QGIS. Thanks for the help. David


title: How to generate an image from OSM on command line?

text: Is there a GNU/Linux command line tool that allows to generate an image file from OSM, preferably as an SVG?


title: Which framework is best to implement an application robustly, using ArcGIS Server JS API?

text: Over on the Esri forum, user Matt Moyles suggested that the approach used in the Esri JS samples isn't suitable for robust development of a web mapping application using HTML, JavaScript and CSS: > ESRI's suggested approach to javascript application development is ancient > and outdated. I would not recommend following the samples. Dojo 1.7 supports > AMD with Asynchronous dependency loading. I would start with the dojo > boilerplate template and "work" the arcgis api into that. The samples are > not suitable for serious application developers. They are mostly just proof > of concept snippets. > > A serious application should be developed with the aid of some sort of > framework to help structure things. I have been using dojox.mvc with great > success! But other options include things like backbone.js, spine.js, or > even javascript MVC. > > Dojo Boilerplate - https://github.com/csnover/dojo-boilerplate   * does anyone agree/disagree with this statement?   * are there any online examples of better approaches to web application development using the ArcGIS Server JS API?