title: How do I plot a line on a map, or obtain the latitude and longitude coordinates of such a line, that is 120km north of the US-Canada border?

text: How do I plot a line on a map, or obtain the latitude and longitude coordinates of such a line, that is 120km north of the US-Canada border, including parts of the border that cross water? (I managed to get coordinates along the land border, but not along the parts of the border that cross water.) I have been trying to figure this one out for a couple of days. Here are just some of the many posts I have consulted: On how to get the lat & lng coordinates along the border: https://stackoverflow.com/questions/12194046/google-maps-api-v3-how-to-get- region-border-coordinates-polyline-data https://stackoverflow.com/questions/9706484/add-search-area-outline-onto- google-maps-result These ones look relevant from the title, but I was unable to figure out how to use them... and I'm not exactly looking to create a parallel line; the new line would not be exactly parallel to the original one, though it would look so in many parts: How to create a line in a given distance to an existing one? Given a flight path given in latitudes and longitudes, how do I find another flight path that is parallel to the path given? Here is the approach that I have tried.   1. Get a bunch of lat-lng coords of the US-Canada border.   2. For each pair of adjacent points, calculate a vector describing a line that joins the two points, and calculate the lat-lng coords of the midpoint of this line.   3. For each of these vectors, find the unit vector perpendicular to the vector in the direction into Canada.   4. For each unit vector, calculate the lat-lng coords of the point 120 km into Canada from the midpoint calculated earlier by adding the unit vector times distance of 120 km to the midpoint. Regarding step 1: I got border data from the National Geographic Data Center. I loaded the appropriate .shp file into qgis. Using the QChainage function (which I found from here) I got the lat-lng coordinates of the border (mostly along the land, not on the shores near Vancouver), and then saved them as a CSV file. For 2, 3, and 4, I tried my best, as I hope you can see from my code below:               function initialize() {           var locations1 = [             [61.48674356, -141.0012206],             [61.58674194, -141.001389],             [61.68674007, -141.0012813],             [61.78673774, -141.0012812],             [61.88672508, -141.0010658],             [61.98671795, -141.0012186],             [62.08670822, -141.001111],             [62.18669614, -141.0007711],             [62.28668686, -141.0007712],             [62.38667513, -141.0003868],             [62.4866706, -141.0008994],             [62.58666236, -141.0004485],             [62.68666005, -141.0004485],             [62.78664608, -141.0012383],             [62.88663255, -141.0017758],             [62.9866291, -141.0014996],             [63.08662212, -141.0014996],             [63.18661862, -141.0015582],             [63.2866163, -141.0015568],             [63.38661351, -141.001667],             [63.48661052, -141.0014997],             [63.58660488, -141.0017284],             [63.68660008, -141.0015566],             [63.78659307, -141.0017752],             [63.88658963, -141.0014999],             [63.98657671, -141.0004492],             [64.08656179, -141.0007687],             [64.18655702, -141.0003869],             [64.28655251, -141.0008974],             [64.38653864, -141.0007691],             [64.48652245, -141.000774],             [64.58650502, -141.0010629],             [64.6864979, -141.0012204],             [64.78649628, -141.001389],             [64.88649628, -141.001389],             [64.98649325, -141.0015581],             [65.08648629, -141.0017761],             [65.18644235, -141.000449],             [65.2864275, -141.0007737],             [65.38641012, -141.0010605],             [65.48640297, -141.0012212],             [65.58639948, -141.0014983],             [65.6863939, -141.0017249],             [65.78634859, -141.0003877],             [65.88633578, -141.0012194],             [65.98633416, -141.001389],             [66.08632904, -141.0017238],             [66.18630654, -141.0003881],             [66.28629372, -141.0012198],             [66.38629025, -141.0014978],             [66.48628678, -141.0015563],             [66.58626503, -141.0007817],             [66.6862579, -141.0012781],             [66.78625233, -141.0017225],             [48.09135328, -89.08774781],             [48.12990592, -88.99551236],             [48.16846578, -88.90327996],             [48.20654583, -88.81083824],             [48.24454325, -88.71836407],             [48.26866166, -88.62203098],             [48.28553203, -88.52349756],             [48.30237648, -88.42497186],             [48.29773516, -88.32843646],             [48.26155109, -88.23525283],             [48.22619115, -88.14176979],             [48.19037419, -88.04846379],             [48.15399495, -87.95536856],             [48.11835363, -87.86199437],             [48.08199413, -87.76889477],             [48.04571355, -87.67578091],             [48.00961195, -87.5825905],             [47.9730714, -87.48952847],             [47.93657073, -87.3964743],             [47.90027541, -87.30336795],             [47.86346579, -87.21042154],             [47.8267048, -87.11745782],             [47.79009911, -87.02444783],             [47.75299358, -86.93163457],             [47.7161161, -86.83874347],             [47.67919225, -86.74585927],             [47.64177662, -86.65316338],             [47.60480505, -86.56032087],             [47.56753398, -86.46757675],             [47.52968697, -86.37503409],             [47.4923828, -86.28233453],             [47.45496578, -86.18964134],             [47.416985, -86.09720241],             [47.37920703, -86.00469466],             [47.34180458, -85.91203491],             [47.30337016, -85.81978872],             [47.26558586, -85.72725087],             [47.22746189, -85.63489552],             [47.1891533, -85.54258762],             [47.15107131, -85.45017422],             [47.11269326, -85.35793446],             [47.0740636, -85.26578074],             [47.03580892, -85.17347346],             [46.99727296, -85.08125112],             [46.95816275, -84.98928576],             [46.91957328, -84.89712782],             [46.84871995, -84.84002049],             [46.75400677, -84.80847183],             [46.65922855, -84.77728631],             [46.59130709, -84.70495143],             [46.5248768, -84.63023577],             [46.4656176, -84.55266116],         ];         // I deleted a lot of points above to save space here in the post          // This is for testing     locations2 = [                 [43.7215978, -79.705822]       ];          centerpoint1 = new google.maps.LatLng(46.4656176, -84.55266116);     centerpoint2 = new google.maps.LatLng(43.7215978, -79.705822);          var x = 1;     if (x==1) {       var centerpoint = centerpoint1;       var locations = locations1;     } else {       var centerpoint = centerpoint2;       var locations = locations2;     };          var mapOptions = {         zoom: 10,         center: centerpoint,         mapTypeId: google.maps.MapTypeId.ROADMAP,         scaleControl: true     };          var map = new google.maps.Map(document.getElementById('map-canvas'), mapOptions);     var marker, i;          for (i = 0; i < locations.length; i++) {           marker = new google.maps.Marker({         position: new google.maps.LatLng(locations[i][0], locations[i][3]),         map: map       });     }          // **********************************************************************************************************************     // Now plot markers X km from the above points     // **********************************************************************************************************************               var R = 6371; // km          function toRadians(Value) {         /** Converts numeric degrees to radians */         return Value * Math.PI / 180;     }          function fromRadians(Value) {         return Value * 180 / Math.PI;     }          // find a perpendicular vector between pairs of points to get to X km from border     for (i = 0; i < locations.length - 1; i++) {         // v1 should be the unit vector between adjacent points, and v2 will be the unit vector perpendicular to this in the direction of farther into Canada         var v1 = [locations[i+1][0] - locations[i][0],locations[i+1][4] - locations[i][5]];          // From http://www.movable-type.co.uk/scripts/latlong.html     var lat1 = locations[i][0];     var lat2 = locations[i+1][0];     var lon1 = locations[i][6];     var lon2 = locations[i+1][7];          var lambda_1 = toRadians(lon1);          var phi_1 = toRadians(lat1);     var phi_2 = toRadians(lat2);     var delta_phi = toRadians((lat2-lat1));     var delta_lambda = toRadians((lon2-lon1));          var Bx = Math.cos(phi_2) * Math.cos(delta_lambda);     var By = Math.cos(phi_2) * Math.sin(delta_lambda);     var phi_3 = Math.atan2(Math.sin(phi_1) + Math.sin(phi_2),                     Math.sqrt( (Math.cos(phi_1)+Bx)*(Math.cos(phi_1)+Bx) + By*By ) );     var lambda_3 = lambda_1 + Math.atan2(By, Math.cos(phi_1) + Bx);          // Get the midpoint between      midpoint_p1_p2 = [fromRadians(phi_3), fromRadians(lambda_3)];          // find a unit vector perpendicular to the line formed by the two points. Am I doing this right? Most certainly not!     var d = 1.2; // d should correspond to 120 km... still need to figure out how to do this     var y2 = lon2-lon1;     var x2 = lat2-lat1;     var k = (y2*y2)/(x2*x2);     var x1 = Math.sqrt(d*k/(1+k));     var y1 = Math.sqrt(Math.abs(d*d-x1*x1));     var v1 = [x1, y1];     var P3 = [midpoint_p1_p2[0]+x1, midpoint_p1_p2[1]+y1];          marker = new google.maps.Marker({         position: new google.maps.LatLng(midpoint_p1_p2[0], midpoint_p1_p2[1]),         icon: {           path: google.maps.SymbolPath.CIRCLE,           scale: 2         },         map: map     });               marker = new google.maps.Marker({         position: new google.maps.LatLng(P3[0], P3[1]),         icon: {           path: google.maps.SymbolPath.CIRCLE,           scale: 4         },         map: map     });          console.log('midpoint');     console.log(midpoint_p1_p2);     console.log('P3');     console.log(P3);     console.log('v1');     console.log(v1);     console.log('k');     console.log(k);          }          }          google.maps.event.addDomListener(window, 'load', initialize);      Well, this produces something: ![Map with border points, midpoints and attempted calculation of points 120 km into Canada](http://i.stack.imgur.com/AdhaV.jpg) Note: the small circles are the midpoints, the big circles are supposed to be the points going 120 km into Canada, and the red icons are the points along the border. Problems:   1. I have to make sure that the calculated points 120 km into Canada are going in the right direction (into Canada) and the right distance (120 km). Sometimes the points go in the wrong direction, and I don't know how to get the distance right.   2. Because of the vertices in the border, the calculated 120km points will form lines that cross. I would need to figure out an algorithm for going through the points to eliminate some to get rid of this. I can program in Python and Javascript, I know how to use the Google Maps v3 API, and I'm just starting out with QGIS. Any help would be greatly appreciated - many thanks in advance!


title: Selecting few biggest polygons within class

text: I have polygon shapefile with objects that have 'class' and 'area' values in attribute table. I want to select five biggest polygons in each class. I can't do it manually because of large number of shapefiles and objects in each shapefile. Is it possible with ArcGIS that I use? If not maybe some other software can do the trick... Thanks in advance for help. edit: I can use file/personal geodatabase or SQL Server geodatabase as well if this is nessesary to make that selection. I will not be able to create python script, I would rather use some SQL selection.


title: Output map of only 50m-100m contours

text: I´ve got the "OS Terrain 50" contours/spot-heights dataset. I would like to output a series of maps with everything above 0m, 50m, 100m etc. from a region. I know almost nothing about GIS, but have a reasonable technical knowledge. There are lots of GIS programs, but which would be best/easiest to do this in. Thanks Edit: Limited (tagged) the program to the recommended QGIS. I am wanting to make a physical model of some land. Terrain outline at different heights above sea level will give some templates to cut out and stack to make an actual model. Hope it´s not going too far off topic.


title: ArcGIS Editor Attributes Dialog

text: Is there a way to update the Editor Attributes dialog to populate/refresh on ArcObject feature selection code? The attributes dialog only seems to respond to a manual view selection or manual table selection. Here is my current feature selection code:               Try                      Dim pFeatSelection As IFeatureSelection = Nothing                 Dim pFeatureLayer As IFeatureLayer                 Dim pQueryFilter As IQueryFilter                 Dim pSelSet As ISelectionSet                 Dim pFCur As IFeatureCursor = Nothing                 Dim player As ILayer                 Dim pEnumLayer As IEnumLayer                      m_App = CType(Hook, IApplication)                 m_MxDoc = CType(m_App.Document, IMxDocument)                 m_mxDocument = CType(m_App.Document, MxDocument)                 m_map = CType(m_MxDoc.FocusMap, Map)                      Dim pActiveView As IActiveView                 pActiveView = m_MxDoc.FocusMap                      m_map = m_MxDoc.ActiveView.FocusMap                 pEnumLayer = m_map.Layers                      Dim number As Int32 = InputBox("FID number")                      player = pEnumLayer.Next                 Do Until player Is Nothing                     If player.Name = "MyLayer" Then                         pFeatureLayer = CType(player, IFeatureLayer)                         pFeatSelection = CType(pFeatureLayer, IFeatureSelection)                     End If                     player = pEnumLayer.Next                 Loop                      pQueryFilter = New QueryFilter                 pQueryFilter.WhereClause = "FID = " & number                      pActiveView.PartialRefresh(esriViewDrawPhase.esriViewGeoSelection, Nothing, Nothing)                      pFeatSelection.SelectFeatures(pQueryFilter, esriSelectionResultEnum.esriSelectionResultNew, False)                 ' Create a feature cursor of selected interstates.                 pSelSet = pFeatSelection.SelectionSet                 pSelSet.Search(Nothing, False, pFCur)                 pActiveView.PartialRefresh(esriViewDrawPhase.esriViewGeoSelection, Nothing, Nothing)                  Catch ex As Exception                 MessageBox.Show("Caught an unspecified error in the calling code: " & vbCrLf & ex.ToString)             Finally      Thank You


title: How to remove/disable last row/column of SRTM GeoTIFF in QGIS?

text: I've started using QGIS yesterday, and one of my main reasons to use it is to apply colormapped SRTM v4 over Google aerial imagery. Everything is working fine, except that SRTM v4 elevation tiles have overlapping last row and last column (6001 x 6001 pixels), and that creates duplicated raster over the edges, which means an ugly, "unacceptable" border over the displayed image. So before I consider editing the files to remove the last row and column, I wonder if there is a way to tell QGIS to ignore this part of the GeoTiff. The image below (usind soft light blending mode) shows the problem over the latitude of -30° ![enter image description here](http://i.stack.imgur.com/4JelO.jpg)


title: How can I create a Style for ArcPad from another GIS as ESRI

text: I have a QGIS-Project with styles. Now I want to export this style to use in ArcPad... ArcPad needs an apl-file. Does somebody know how to convert the style from QGIS into a apl-file?


title: How to provide a WFS in multiple projections

text: I would like to set-up an editable web feature service (WFS-T) for serving my geospatial data. The data tables are in a PostGIS database on the same server. For security reasons my server admin doesn't allow remote connections to the Postgresql database. So I'm looking at using the WFS service. Clients will view and edit the data in QGIS or using a OpenLayers javascript client. I succeeded in installing the TinyOWS software and setting up the transactional WFS. However, I have an extra requirement: All data layers must be available in 3 projections: EPSG:3857, EPSG:3067, EPSG:4326. Right now each table has a column the_geom with coordinates in the Finnish projection EPSG:3067. What would you recommend to add the support for the other projections? Does the server like TinyOWS do "on the fly" re-projection on each request, or do I need to add extra geometry columns with coordinates in the other projections to each of my tables?


How to add new columns with attributes in existing Qgis layer? Doing the "same thing" here described


title: How to add information to an attribute table?

text: I'm working with qgis 1.8.0 and did the following: 1.-Divided my shapefile into 830 features (small rectangles). 2.- Added to my shapefile several points containing weather information using the "add delimited text layer" tool. 3.- Did some weather data interpolations. 4.- Used the "point sampling tool" to collect the interpolated information for each and every one of the 830 features of the shapefile. QUESTION: is there any way to add/import the information obtained with the point sampling tool to each of the 830 features of my original shapefile? I tried to modify the attribute table but couldn't find any way to paste information to a new column.. any ideas? Thanks!


title: Mapfile SYMBOLSCALEDENOM and SIZEUNITS (miles) 200: Keep SIZEUNITS while circles change size based on zoom

text: I've been trying to understand how to keep the `SIZEUNITS` at 200 while the circles appropriately adjust based on zoom level. Thought I had this working properly but turns out I did not. You can see here that the circles start out quite small at the default zoom level when you load the OpenLayers map. As you zoom in, the circles aren't changing size. They should be as `SIZEUNITS` is set to 200. Of course, since you are zooming in, it seems like the circles are getting bigger on your screen, but they really are not. When you get pretty zoomed in, what should happen is that the circles should be quite large (as we want to keep the radius of 100miles...set `SIZEUNITS` to 200 as this is the diameter). However, that is not the case as you can see here. Is it necessary to use `SYMBOLSCALEDENOM` with `SIZEUNITS` to acheive what I want: keeping the 200mile diameter whilst the circles are changing size based on zoom? If so/if not, how do I go about this? Link to mapfile layer here.


title: "locate features along routes" tool results in empty table!

text: I have no idea why the result table of the tool "locate features along routes" is empty. I have a stream network and a cross-section layer. Here is what I did: 1- Converted my stream polyline shapefile to Route using Create Routes (Linear Referencing) tool. 2- tried to get the location of the cross sections along the route using locate features along routes tool ![enter image description here](http://i.stack.imgur.com/5ujbF.png)


title: Create Voronoi/Thiessen polygons in QGIS around 3 points

text: Running QGIS 1.8.0 I have three bus compounds. I would like to generate Voronoi/Thiessen polygons around those three points. Following Voronoi/Thiessen adding an all school points will allow us to determine which school should be served by which bus compound. Creating Voronoi: ![enter image description here](http://i.stack.imgur.com/s5Mrj.png) It works very quickly and results in: ![enter image description here](http://i.stack.imgur.com/DTF6J.png) Obviously not what I expected. Can anyone shed some insight? Thanks


title: How to import binary arcgis maps into R

text: I have imported arcgis rasters into R. The rasters are binary maps and the legend in arcgis shows this to be binary as the colors are only 0 and 1. Now, when I import the raster into R it is no longer binary, as there are 5 numbers in the legend. Does anyone know why this occurred and have a solution to fix this? My other question is, when I load a raster in arcgis, how do I effectively change the legend and then save the raster as a different name? This is m R code for importing the raster and then making a 4 paneled plot. I have no idea how to change the legend in R.               r5<-raster("C:/users/k/desktop/g/A1b00-59/dican00")     r6<-raster("C:/users/k/desktop/g/A1b50-59/dican50")     r7<-raster("C:/users/k/desktop/g/A1b50-59/dican59")     r8<-raster("C:/users/k/desktop/g/A1b00-59/dicancalc")               par(mfrow=c(2,2))     plot(r5,sub="2000 Distribution",col.sub="blue")     points(Dican,col='blue')     title(main="A1b scenario Panic Grass Distributions",col.main="red")     plot(r6,sub="2050 Distribution",col.sub="blue")     points(Dican,col='blue')     plot(r7,sub="2059 Distribution",col.sub="blue")     points(Dican,col='blue')     plot(r8,sub="Net Habitat Change Panic Grass",col.sub="blue")     points(Dican,col='blue')      Thanks in advance


title: Advanced android app for mapping and GPS logging

text: I am looking for an android app able to log GPS data like coordinates, altitude, current speed and would be nice if also the GPS accuracy/number of satellites. Many GPS apps are able to display these data, but not to log them. Reason is that while doing a monitoring (no matter if helicopter, car or on foot), there is a variable accuracy of the coordinates and so this also influences the measurement results and the output maps. So it would be helpful to know these data to take this into consideration. I am currently using Oruxmaps, but it logs only date/time, position and elevation and not the other data. Just found Antimap and Geopaparazzi which I will surely try, but would also appreciate your experiences and recommendations. thanks PS: need to work in offline mode :)


title: Trying to combine 29 zip code service areas into 10 or fewer service areas

text: I am tasked with coming up with 10 or fewer zip code service areas based on 29 pre-existing zip code service areas. The zip codes cover 5 counties and many service areas, obviously, overlap. The client wants fewer service areas based on common areas or other meaningful determinations. I'm thinking of some kind of hot spot analysis? Open to all suggestions!


title: Loop in ModelBuilder for compositing bands of multiple Images

text: I'm attempting to automate the "Composite Bands" tool somehow in ModelBuilder, but haven't been able to figure out an effective method.   I have landsat 7 data that is separated into folders and each holds 7 bands (tif). Basically I want to generate a model or script that composites all bands in the raster output for each location, and then saves the result in a output folder. I've prepared a ModelBuilder that take each band and save it in its correspondent folder. However, the input raster just read 1 band, not all 6 bands. Any Solutions, maybe Python? ![Model Builder](http://i.stack.imgur.com/xZkTQ.jpg) ![Composite Bands which just i raster uploaded :\(](http://i.stack.imgur.com/U7X39.jpg)


title: Geoexplorer ../proxy/url=? fails with 'moz-nullprincipal' when run in Tomcat7, not Jetty7

text: There's a remote WMS service (not running on port 80) that can be accessed from a machine - it's getCap URL can be reached by a browser. When attempting to add this WMS source in a Geoexplorer application (running in a Tomcat 7 of default settings, on any port) (via add layer) one gets "no element found, moz-nullprincipal" when Geoexplorer tries to do its ../geoexplorer/proxy/url?=(WMS getCap URL) probe. The /proxy/url?= is presumably to avoid cross-domain issues for XMLHTTPRequests such as this one. Cross-domain issues are discovered by the browser, as I understand it. And yes, Chrome and IE do the same thing more or less. _However_ , if the same Geoxplorer application runs in Jetty 7 of default settings and the same port as the Tomcat(not simultaneously obviously)then this issue is nonexistant. Why would default Tomcat limit a request like this (in all browsers) but not default jetty (in all browsers)? It might be worth mentioning that Tomcat has no issues with other WMS services (on other remote domains, but running on port 80). Many thanks for any advice. Edit. Default Resin 4 acts about the same as Tomcat, though instead of throwing an error it never seems to time out the request. So are the seemingly jetty-specific files of the Geoexplorer application the cause for this working in jetty? Like geoexplorer\WEB-INF\app\httpclient.js and perhaps the jetty- client jars.


title: Esri ArcGIS Server Maintenance

text: I am looking at our quote from Esri and am not sure on the difference between two items on the quote. They are: ArcGIS for Server Workgroup Standard One Core Additional Migrated Maintenance Start Date: 04/19/2013 End Date: 04/18/2014 ArcGIS for Server Workgroup Standard Up to Two Cores Migrated Maintenance Start Date: 04/19/2013 End Date: 04/18/2014 They look one in the same to me. They are quoting me for two of the first one and one of the latter one. Can anyone tell me what the difference is?


title: How to create a heat map or density map from stacked polygons in ArcMap 10

text: I have a polygon shape file that has many overlapping features. Each feature represents an area used by polar bears. So some of the polygons overlap and some don't, and there are separate clusters of overlapping areas. I would like to get a measure of how much overlap there is in a particular area. But I have no clue how to do this. I'm pretty new to GIS and am still learning the ropes, it just seems I can't find an answer. I'd appreciate any help! Maha


title: Selecting a valid format file

text: There has been an error when uploading a shapefile I want to start with some imported data. Upload a file.               Please, select a valid format file          # when upload a zip file:      > Oops! There has been an error   >  Unknown   >  Code 99999   >  Begin:   >  TRANSITION: preprocessing => uploading, 2014-01-07 02:59:15 -0800   >  getaddrinfo: Name or service not known   >  Table error, #


title: When will QGIS use the WMS supplied legend images ? is there a workaround?

text: I use a lot of wms services with complicated maps : a legend is essential. Qgis displays none of those legends. The URL is quoted in the metadata. example WMS http://geo.agiv.be/ogc/wms/product/ANB?/Grnkrt13 layer "Groenkaart Vlaanderen 2013" legend url : http://geo.agiv.be/ogc/wms/product/ANB?request=GetLegendGraphic%26version=1.3.0%26format=image/png%26layer=Grnkrt13


title: error of data conversion value getting while importing shapefile in shp2pgsql-gui on ubuntu

text: Connecting: host=localhost port=5432 user=postgres dbname=tremblement_de_terre Connection succeeded. Connecting: host=localhost port=5432 user=postgres dbname=tremblement_de_terre ============================== Importing with configuration: m_cap point, public, geom, /home/jean/Documents/M_CAP Point.shp, mode=c, dump=1, simple=1, geography=0, index=1, shape=1, srid=0 Shapefile type: MultiPoint PostGIS type: POINT[2] Unable to convert data value to UTF-8 (iconv reports "Chaîne multi- octets ou étendue de caractères invalide ou incomplète"). Current encoding is "UTF-8". Try "LATIN1" (Western European), or one of the values described at http://www.postgresql.org/docs/current/static/multibyte.html. ## Shapefile import failed. After trying to know the encode setting with the ubuntu shell command file -f .shp,the result is something like file -f M_CAP\ Point.shp : ERROR: cannot open `' (No such file or directory) : ERROR: cannot open`' (No such file or directory) : ERROR: cannot open `' (No such file or directory) ��OA�\���?A�: ERROR: cannot open`��OA�\���?A�' (No such file or directory) ��OA: ERROR: cannot open `��OA' (No such file or directory) ��OA: ERROR: cannot open`��OA' (No such file or directory) �5A�X=�: ERROR: cannot open `�5A�X=�' (No such file or directory) �5A�X=�: ERROR: cannot open`�5A�X=�' (No such file or directory) �5A�X=�: ERROR: cannot open `�5A�X=�' (No such file or directory) ��XCAp���_A~�: ERROR: cannot open`��XCAp���_A~�' (No such file or directory) * * * When I also try to open this file with QGis as a Vector layer, using UTF-8 encoding, It opens clearly. What happens now when trying to import it to postgres using shp2pgsql-gui? Please help me, the original file have been created on windows, and I am using Ubuntu 12.04 Lts. Thanks a lot


title: Include Python Script with ArcObjects Addin

text: I'm building addin for 10.1 and need a method (ConvertToMSD) that's only available in ArcPY. I've figured out out how to execute a python script from an ArcObjects addin: Calling arcpy/python from .NET? Now I want to figure out how to include that Python script within in build addin. If it's not included in the addin then I have to tell the user to download the Python script and go drop it in folder x/y/z which isn't very user friendly. Is this possible?


title: error occured when importing file on cartodb local host

text: Logged in as development user. Not able to import a .zip or a csv file. Ruby console has this error log               1.9.2p320 :004 > DataImport[1].log          2014-01-07T05:44:44Z || Using database connection with {:adapter=>"postgres", :encoding=>"unicode", :host=>"localhost", :port=>5432, :database=>"cartodb_dev_user_1_db", :username=>"postgres", :user=>"development_cartodb_user_1"}     2014-01-07T05:44:44Z || ogr2ogr output:    -nlt PROMOTE_TO_MULTI: type not recognised.               2014-01-07T05:44:44Z || ogr2ogr exit code: 256     2014-01-07T05:44:44Z || CartoDB::Importer2::LoadError     2014-01-07T05:44:44Z || CartoDB::Importer2::LoadError          2014-01-07T05:44:44Z || ERROR!      The error message screen shows:               Error 2001. Not able to solve from the error log.


title: Writing an algorithm to find the most efficient line angle at which to draw lines

text: I've written a tool that fills a polygon in ArcMap with parallel lines at a user specified distance apart and at a user specified angle, written in vb.net and arcobjects. What i'm looking to do next is write a tool that calculates the most efficient angle at which to draw these lines (i.e. the angle which fills the polygon with the least number of lines, or greatest average line length), given the specified line separation distance. For example if the polygon was a simple rectangle, then the most efficient line angle to fill the rectangle would be parallel to the long axis, as it would take more lines to fill the polygon if they were drawn parallel to the shorter axis. For more complicated shapes I have little knowledge of how to achieve this. The only way I can think of doing it is to write a loop that creates lines at every angle and then pick the output with the greatest average line length, but that would take too long (drawing each set of lines takes a good 5-10 seconds depending on the number of lines). I know it's a longshot but does anyone have any ideas? Probably quite a tricky one...


title: Finding the distance between points along a river (Linear referencing)

text: I am new to gis and am struck. I have a network of rivers which is made up of small linestrings (polylines). In a different shapefile, i have locations of several stations on a river. I need to find distance along the river from each station till the end of river. The problem is that I don't have one complete linestring of river. And there is no common attribute I can use to union or merge the linestrings. Also, some stations are on tributaries. ![enter image description here](http://i.stack.imgur.com/U1clE.jpg) So the bottom line is I need to make a large network of rivers including tributaries and perform liner referencing. But I don't know how to do this. Any help is deeply appreciated. Please explain as clearly as possible. Please comment if you need more information. Thanks in advance. BS


title: How to redefine the nodata value into zero in QGIS?

text: I have an image which i cliped with an AOI polygon. The contour area was set to pixel value = 0 and I have a nodata value of -32768. I want to redefine the nodata value into zero. I used to do this with ArcMap easily in the raster calculator, but in QGIS I don't see how i could do it. Any ideas?


title: ArcGIS 10 SP2 issue with AggregatePoints

text: ESRI persuaded us to upgrade to SP2 and now we cannot use the code we have, as it doesn't work. According to the download material, NIM063853 is when AggregatePoints doesn't create the derived table. This worked in SP1, now doesn't in SP2? What in gods name is going on? Has anyone else experienced this? Any workarounds?


title: PostGIS SQL NOT EXISTS Statements?

text: I'm trying to query my PostGIS database for ID of people who did X but did not do Y. I've attempted the query below, but it does not return any results. I'm very new to SQL so I may be overlooking something. Any thoughts? Thanks for your help.               SELECT *     FROM Table1 INNER JOIN Table2     ON Table1.ID = Table2.ID     WHERE Table1.Date = '2008-11-04'     AND NOT EXISTS (         SELECT *         FROM Table1 INNER JOIN Table2         ON Table1.ID = Table2.ID         WHERE Table1.Date = '2012-11-06')


title: What is the max. length of string for GP output?

text: I have an XML as the output of Geo-processing tool. I am keeping it as a file right now but another developer wants it as return value. So what is the maximum possible length of this String output? Is this variable type or fixed type based on MS terminology?- Thanks


title: How to group layers in ArcEngine 10

text: I was wondering if there is a way I can code the TOC control so that I can select multiple layers and **GROUP** them together. So if i uncheck the group, all layers will be unchecked and be invisible. Thanks and Regards, Kevin


title: How to setup proxy and use "POST" method when executing Geoprocessing Service in ArcGIS Server Javascript API

text: I am passing five featuresets to a ArcGIS Server Javascript API. When I do that I get 'Deferred Object not defined error'. Then I tried to use dojo.Deferred(). At that point the error goes away. But nothing happens. In Fidller, no request information is logged. Then I removed some of the records from the featureset and I can see the request is logged in the fiddler and the request reaches the server. Processing fails, because all required data is not available. Can some one help understand how to setup proxy server in ArcGIS Server 10.1 and use POST request it in JS API for Synchronous and Asynchronous Geoprocessing? Also I can someone show me how to use dojo.deferred() when using GP Service. Thanks Cyril


title: Rules of Nodata value?

text: What are the arithmetic rules that govern Nodata values? Regardless of how these null values are defined? I am using ArcMap Raster Calculator. I am thinking specifically of operations like 1*Nodata (or 1*null), 1/Nodata, 1+Nodata,1-Nodata, Nodata-1 etc...


title: How to make sure that geometries copy the relationship of corresponding tables?

text: A make a database in PostgreSQL (9.2.2)/PostGIS (2.0.1). I'm new to it, so I'm not sure how to make sure that geometries corresponding to main table (polygons or points) cover the geometries corresponding to the dependent table (always points). In ArcGIS, I would simply apply a topology rule, so I tried to learn more about topologies in PostGIS, but what I found only confused me. Now I don't know how to do it through topology and I'm even not sure whether this way is a good solution, since most PostGIS topology application focus on other topology rules than I need. Another alternative I found was a SELECT like this: SELECT gid, the_geom FROM foo_geometry WHERE ST_CoveredBy(foo_geometry, bar_geometry)='f' where foo_geometry corresponds to the dependent table and bar_geometry to the other one. I can also use a cursor for this. What solution is the best provided that   a) as a beginner I prefer solutions easy to understand and implement   b) there won't be many geometries and the validation will be done just once in a time as a maintenance process, so performance is not crucial here   ?


title: PLT track in PostGIS

text: What are the best practices for importing and aggregating tracks (in my case in 'plt' format) in a PostGIS table, grouped by user. Are there any tools/db models/tutorials available?


title: Save fails on raster object created from NumpyArrayToRaster

text: I am trying to save an ESRI grid built from a numpy array and am getting an error I've never seen before: > Unhandled Exception: System.IO.FileNotFoundException: Could not load file or > assembly 'GpMetadataFunctions, Version=10.2.0.0, Culture=neutral, > PublicKeyToken=8fc3cc631e44ad86' or one of its dependencies. The system > cannot find the file specified. File name: 'GpMetadataFunctions, > Version=10.2.0.0, Culture=neutral, PublicKeyToken=8fc3cc631e44ad86' at > System.Reflection.Assembly._nLoad(AssemblyName fileName, String codeBase, > Evidence assemblySecurity, Assembly locationHint, StackCrawlMark& stackMark, > Boolean throwOnFileNotFound, Boolean forIntrospection) at > System.Reflection.Assembly.InternalLoad(AssemblyName assemblyRef, Evidence > assemblySecurity, StackCrawlMark& stackMark, Boolean forIntrospection) at > System.Reflection.Assembly.InternalLoad(String assemblyString, Evidence > assemblySecurity, StackCrawlMark& stackMark, Boolean forIntrospection) at > System.Reflection.Assembly.Load(String assemblyString) at > GetManagedType(Char* _assembly, Char* _path, Char* _className) > > WRN: Assembly binding logging is turned OFF. To enable assembly bind failure > logging, set the registry value [HKLM\Software\Microsoft\Fusion!EnableLog] > (DWORD) to 1. Note: There is some performance penalty associated with > assembly bind failure logging. To turn this feature off, remove the registry > value [HKLM\Software\Microsoft\Fusion!EnableLog]. I'm using an alternative workflow for now, but I'd still like to try to solve the problem with saving a raster from a numpy array. Here is what I was trying to do:   1. read in an ESRI grid,   2. create an numpy array,   3. do some manual reclassification, and   4. save the processed array to a new ESRI grid The process proceeds as expected through step 3 and fails upon step 4 _AFTER_ the array has been converted to a raster object (see code below). It is when the `save` method of the raster is called that the error is encountered. Any help in tracking this down would be greatly appreciated. PS - I got the same error when trying to create a numpy array from geo-tiff rasters (failed on calling `arcpy.RasterToNumPyArray`); I was able to get around this by copying the original .tifs to ESRI grid format first. Here is my code:               import arcpy     import numpy as np          arcpy.env.overwriteOutput = True     arcpy.env.workspace = 'SwbLandUseParameters.gdb'          # List tables to use in the reclassification     tables = arcpy.ListTables('SwbLandUseParameters*')          # List land use rasters     lu_rasters = arcpy.ListRasters('nlcd*', 'GRID')          for lu_raster in lu_rasters:              print 'Processing raster', lu_raster              # Get raster properties         r = arcpy.Raster(lu_raster)         ll = r.extent.lowerLeft         x_cell_size = r.meanCellWidth         y_cell_size = r.meanCellHeight         ncol = r.width         nrow = r.height         year = lu_raster[-4:]              # Convert the raster to a numpy array         lu = arcpy.RasterToNumPyArray(r, ll, ncol, nrow)              # Create a suite of reclassified rasters using parameter tables         for tbl in tables:             temp = np.copy(lu)             with arcpy.da.SearchCursor(tbl, '*') as cur:                 for row in cur:                     print row                     old_val = row[1]                     new_val = row[2]                     temp[temp==old_val] = new_val             new_raster = arcpy.NumPyArrayToRaster(temp, ll, x_cell_size, y_cell_size, 0)             new_raster_name = tbl + '_' + year             print '  ', new_raster_name             new_raster.save(new_raster_name)


title: Adding Facebook and Twitter sharing to CartoDB infowindows

text: Trying to add share buttons to output something like: {{region}} text {{region_variable}} text www.google.com The below works very intermittently for Twitter, but when I share the URL, the map won't load.               Tweet           Not even sure where to start with doing the same with Facebook. Any help appreciated.


title: How to get rectangle with two values in converting polygon to raster?

text: I have 112 vector polygon shapefiles. I want to convert every single one to raster layer but instead of getting layer with same shape as vector polygon, I want to get rectangle in which cells that match vector polygon have value 1, and all other cells in rectangle have value 0. And that rectangle for every raster layer should have same size because I want to overlap and sum all 112 raster layers together to get final raster layer in which I will have value gradation among cells between 1 (where is only one of primary raster layers lay) and 112 (where all 112 raster layers overlapping). Is there a simple way to do that? I don't know how to get rectangle in converting polygon to raster. I'm using ArcMap 10.0 Here is scheme what I want: http://i40.tinypic.com/2v11abp.jpg. I have many polygons that overlapping but all are inside of boundary of Croatia shape. And I want every polygon have the same extent when converted in raster (red rectangle). At the end I should get final raster layer that look like this schematic: http://i39.tinypic.com/33o4juh.jpg.


title: Handle a great amount of Markers with OpenLayers

text: My general Problem is to maintain a lot of markers. For the beginning I want to display all Bus Stops of Berlin, Germany and surroundings as Markers. I already got that I have to use a strategy like Startegy.BBOX for performance issuses (my dataset contains nearly 13000 Bus Stops). I found an Example where somebody gets his marker positions from a databases, parses it to a giant string and uses it with Format.Text. I guess this would work for me too, but I still wonder if thats the best way to do it. What if my dataset would be much bigger (maybe to display all bus stops worldwide)? I guess putting the whole content of the tatabase to a giant string and use Format.Text woulden't be a good idea than. What would be the best way to do it? Is there a simple way to only load marker positions when they are actually needed instead of preloading them all?


title: Some gridlines missing in QGIS Print Composer

text: I'm very new to this so please excuse my ignorance. I've been working with OS data and my own shape files and had no problem with print composer. Suddenly my OS gridlines are not appearing properly. Some are there, some aren't. The map is fine until I go into print composer. What have I done wrong? Thanks


title: How can I optimize this point-linking query?

text: I have the following query, and am wondering how it could be optimised to run faster:               WITH     gps_nodes_1 AS ( SELECT "NODE_ID", "geog" FROM gps_nodes WHERE ST_Distance("geog",ST_GeogFromText('POINT( -1.54790 53.79424 )')) < 40000 ),     gps_nodes_2 AS ( SELECT "NODE_ID", "geog" FROM gps_nodes WHERE ST_Distance("geog",ST_GeogFromText('POINT( -1.54790 53.79424 )')) < 40000 )     SELECT         row_number() over() as LINK_ID,         gps_nodes_1."NODE_ID" as source,         gps_nodes_2."NODE_ID" as target,         ST_Distance(gps_nodes_1."geog",gps_nodes_2."geog")::numeric(7,1) as length     INTO         gps_node_links_neighbouring     FROM         gps_nodes_1     INNER JOIN         gps_nodes_2     ON ( ST_Distance(gps_nodes_1."geog",gps_nodes_2."geog")::numeric(7,1) BETWEEN 0 AND 20 )      A description of what it does: Table "gps_nodes" consists of ~2 million GPS positions. It has two columns, "NODE_ID" (an auto incrementing id) and "geog" which is a PostGIS geography point. For each GPS node, I want to create a link to all other GPS nodes within 20 metres, and store that link (with start id, end id, and length) in a new table. To start with, rather than try doing this for all GPS nodes, I want to limit it to those within 40 kilometres of a set point - which is around 90,000 records in the gps_nodes table. So I started by creating two CTEs of the GPS points within 40km of required position. Then creating an inner join between the two CTEs, where the join condition is that the points are between 0 and 20 metres apart (need to be greater than zero distance, otherwise you always get a point linked to itself). I freely admit that I'm an SQL novice and there may be far better ways of doing this! Can anyone suggest a more optimised version of this query ? Many thanks.


title: How to create a color ramp to change the color of a svg image in arcgis javascript api

text: I have a svg image to create a PictureMarkerSimbol for a feature. How to create a color ramp to change the svg color based on a criteria, like a attribute on my feature? The ideia is like that:               var mySymbol = new esri.symbol.PictureMarkerSymbol("../resources/img/myimage.svg", 32, 32).setOffset(0, 8);     var renderer = new esri.renderers.SimpleRenderer(mySymbol);          renderer.changeColor ( field, colors, criteria )      Thanks in advanced


title: ECW files to geoserver 2.4

text: I'm trying to add .ecw files (20 orthophotos) to geoserver 2.4. Geoserver is running on an ubuntu 12.4. Vserver. I'v already installed the GDAL plugin. Is this the right format for the import:`ERDASImg - Erdas Imagine Coverage Format` ? By adding a ecw file this message shows up:               Could not list layers for this store, an error occurred retrieving them: The Provided input is not supported by this reader      Do you have an idea what's wrong? Thanks for help!


title: How do I create an IFeature object from an IPolyline or ILine object

text: I am having some problems with some code i am trying to write. The application i am creating allows the user to click on the map and a point will be createD, when they click again to create another point a line is created between the two points. The code i have got at the moment is               for (int i = 0; i < pointCollection.PointCount - 1; i++)     {         //create the line from the points         ILine pLine = new Line();         pLine.PutCoords(pointCollection.Point[i], pointCollection.Point[i + 1]);              IPolyline newPolyline = new PolylineClass();         newPolyline.FromPoint = pLine.FromPoint;         newPolyline.ToPoint = pLine.ToPoint;              //create a feature ****STUCK HERE****         IFeature newline = new PolylineClass() as IFeature; ;         newline.Shape = newPolyline;                   geometriesToUnion.AddGeometry(newline.Shape as IGeometry);     }     resultPolyline.ConstructUnion(geometriesToUnion as IEnumGeometry);      Has anyone got any idea if this is possible or a another way of doing it. Thanks


title: Right click context menu in OpenLayers 2

text:
My code is this in Openlayers 2 I'm displaying two layers in a single map one with the world map and the other with my node information ..   When I right click(on node layer only) I want to display contextmenu in the openlayers, but I'm getting context menu all over the map. Someone please help me with the bug. At least help me how to get right click context menu on one layer when there are multiple layer in the map panel. I used jquery contextmenu (html div element)                              <title>Layer Tree</title>                              <link rel="stylesheet" type="text/css"         href="../resources/ext/resources/css/ext-all.css" />          <link rel="stylesheet" type="text/css"         href="../resources/ext/resources/css/xtheme-green.css" />          <link rel="stylesheet"         href="../resources/openlayers/theme/default/style.css" type="text/css" />          <link rel="stylesheet" href="../resources/app/theme/style.css"         type="text/css" />     <link rel="stylesheet"         href="../resources/geoext/resources/css/gxtheme-green.css" type="text/css" />                                        <link href="jquery.contextMenu.css" rel="stylesheet" type="text/css" />          <link href="jquery.contextMenu.css" rel="stylesheet" type="text/css" />
                 GNIP


title: Axis of rendered shape files not showing numbers that are latitude and longitude

text: I am using `rgdal` to render some shape file. The first is a set of building zones from NYC Open Data. The second is the NYC coast line from OSM. The Zone files seems to generally render an X axis between 850000 and 1150000. And a Y axis with numbers from 150000 and 250000. The OSM coastline renders an X-axis between -8400000 and -8100000 (note the difference in order of magnitude as well as being negative). The OSMs Y axis is between 4900000 and 5050000. I'm using this as my example and I don't the author doing anything special to format the axes. My code is here:               # render a plot based of shapes in nycgiszoningfeatures_201311_shp/     # shapes downloaded from https://nycopendata.socrata.com/api/geospatial/kdig-pewd?method=export&format=Shapefile     require(rgeos)     require(maptools)     require(rgdal)          # OSM Coastline of NYC Data from http://osm-extracted-metros.s3.amazonaws.com/new-york.coastline.zip     mapBase <- readOGR('new-york-city.coastline', 'new-york')          shapeDir <- 'nycgiszoningfeatures_201311_shp'     shapeFiles <- dir(shapeDir, '\\.shp$')     shapeLayers <- sub('\\.shp$', '', shapeFiles)     # to see the files sprintf("%s/%s", shapeDir, shapeFiles)     # TODO: I'm told for loops are evil in R. Find a better way     for (i in 1:length(shapeLayers) ) {         map <- readOGR(shapeDir, shapeLayers[i])         #plot(mapBase, col="khaki2",border="grey", axes=TRUE, xlim=c(-8300000,-8200000), ylim=c(4900000, 5050000))         #plot(map, add=TRUE, col="gray",border="blue")         plot(map, col="gray",border="blue", axes=TRUE)     }      How do I get these shapes to plot with longitude and latitude numbers for the axes instead of those numbers.


title: tilecache config for generation mapnik tile

text: i have a xml stylesheet like below   map:                 <Map srs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over" background-color="#b8dee6" >      a cover layer in :                 <Layer name="cover" srs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over">      and another like :                 <Layer name="line" srs="+proj=utm +zone=39 +ellps=WGS84 +datum=WGS84 +units=m +no_defs">      and i can render pic with mapnik via this script:                      import mapnik     ll = (53.8889,26.481,54.0681,26.5971)     z=12     imgx = 50 * z     imgy = 50 * z     m = mapnik.Map(imgx,imgy)     mapfile = 't'     prj = mapnik.Projection("+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs +over")     c0 = prj.forward(mapnik.Coord(ll[0],ll[1]))     c1 = prj.forward(mapnik.Coord(ll[2],ll[3]))     print c0,c1          mapnik.load_map(m,mapfile)     bbox = mapnik.Envelope(c0.x,c0.y,c1.x,c1.y)     m.zoom_to_box(bbox)     mapnik.render_to_file(m, 'test10.png')           but i cant generate tile with tilecache   tilecache.cfg:                    type=Mapnik     mapfile=/home/test.xml     srs=EPSG:900913     bbox=53.8944,26.4911,54.054,26.5820           please help. thanks


title: Merging data using FME

text: I am a newbie in FME and is trying to accomplish a similar implementation like this: merging data in ArcSDE But we are still in the proof of concept stage.So we are just starting of with 2 non-spatial databases (master and updated) being populated by FME using a single shapefile. Initially both database are equal. Also, we added a FeatureID column using the Counter transformer in FME. We made a couple of updates to the updated database so that we can try sync up the updates to the master database. We are using this sample as reference: http://fmepedia.safe.com/articles/Samples_and_Demos/Creating-a-Unique- identifier-using-CRCCalculator. During the population of the master database. We added a "CRC" column which we set as a varchar(100). Before writing to the master database, there is a CRCCalculator added so that we can populate the CRC column. Syncing process: - we have 2 readers which is set for the update db and master db. - the update db will pass through a CRCCalculator first and will then be passed to a FeatureMerger as a supplier. - the master db will pass through a AttributeRenamer for the CRC column and will then be passed to the same FeatureMerger as the requestor. CRCCalculator parameters: - Attributes = all of the columns - Calculate CRC = Geometry and Attributes - CRC output attribute = CRC I have added an inspector in the "Complete" output port of the FeatureMerger and noticed that all of the CRC(from update db) and CRC_Master(from master db) are not equal, even though there are no changes to those rows. So FME is thinking that all of the data have changed. What am i doing wrong in the syncing process? Here is a screenshot of how the syncing process looks in FME: ![enter image description here](http://i.stack.imgur.com/9oRGy.png) **UPDATE** - So i tried to test this again and added another CRCCalculator for the Master database. So now i can see 3 CRC columns in the FeatureMerger results. I noticed that the CRC value that was created in the translation matched and the CRC value that was saved in the Master database did not.   * Now im wondering if the CRCCalculator changes something in every translation. thanks.


title: How to print tiled mapservices?

text: I have created a few GeoExt webmapping projects using Ubuntu10.04 / Tomcat6 / Geoserver /Postgis but have no good solution regarding printing the maps. Using GeowebCache I am not able to export to pdf using the mapfish module in geoserver due to tiling errors. If I do a simple html print of the mappanel then the toolbar in the mappanel is also printed and distorted....or I risk being sued by google if the maps are printed. I need to be able to serve the maps with geoserver - ie not as tms or filebased tiling methods so that the tiles can be loaded as an ogc service. Can anyone help? yours, Rob


title: How to import GDAL in Jython?

text: I'm trying to use jython with GDAL on Windows. Python and GDAL are installed through MS4W and work perfectly together. Jython is installed manually, but I can not import GDAL. I tried to update sys.path on Jython to point to all python sys.path. But it displays the following error message:               > Traceback (most recent call last):         > File "C: \ myjythonscript.py",     > line 38, in <module>     >      OSGeo import    File "C: \ OSGeo4W \ apps \ Python25 \ lib \ site-packages \ OSGeo \ __init__.py", line 24, in <module>     >      _gdal import ImportError: No module named _gdal      Has anyone worked with Jython and GDAL? Any ideas to solve the problem? Thank you!


title: Why does OpenLayers display red tiles for USGS map?

text: I am trying to display the USGS national map via OpenLayers: http://viewer.nationalmap.gov/viewer/ OpenLayers is able to display the USGS map. However, when I zoom in, the tiles turn red. Could someone explain how to fix this? Is there a USGS restriction on zoon levels? Thanks


title: Loading a raster using Openlayers and Mapserver

text: I'm trying to load a raster into my browser using Openlayers. I can examine my WMS raster layer using OpenJump, but I'm not able to get it to load into my browser. My mapfile is as follows:               MAP               ###########         NAME         "ny_pop"         SHAPEPATH     "/v0.3/mapserver/"                                                                                UNITS         meters         EXTENT        -8269337 4955373 -8190365 5005745          SIZE          1000 1000         STATUS        ON         DEBUG         ON         IMAGETYPE     "png"              ###########             PROJECTION             "init=epsg:900913"         END              ###########     OUTPUTFORMAT       NAME "png"       DRIVER AGG/PNG       MIMETYPE "image/png"       IMAGEMODE RGB       EXTENSION "png"       FORMATOPTION "GAMMA=0.75"     END              ###########             WEB                          METADATA                 "wms_title"           "WMS Pop Tif"                 "wms_onlineresource"  "http://mysite.com/mapserv.cgi?map=raster_1.map&"                             "wms_enable_request"  "GetMap GetFeatureInfo GetCapabilities"             END         END              ###########         LAYER             NAME      "ny_pop"             DATA      /home/mysite/webapps/htdocs/v0.3/mapserver/ny_pop.tif             STATUS    ON             TYPE      RASTER                METADATA                 "wms_title"  "WMS NY Raster"                 "wms_srs"    "EPSG:900913"             END         END              ###########      END      Everything (map, raster and mapfile) are using EPSG:900913, and I have an existing map (centered around New York) which I am adding the layer `rasterLayer` to. My javascript is as follows:               function generateRaster(){              console.log('get map bounds');         bounds = map.getExtent();              rasterLayer = new OpenLayers.Layer.WMS(                 "ny_pop",                "http://mysite.com/mapserv.cgi?",                 {map: 'raster_1.map'},                {'maxExtent': bounds}            );         console.log(rasterLayer);         map.addLayer(rasterLayer);          }      I'm trying to limit the raster to the bounds of the map using the `{'maxExtent': bounds}` command, but I also tried with just the map name and URL and it did not work either. I'm still trying to get my head around the `.map` format, so there could be some mistakes there also. However, as I could view the raster in OpenJump I assumed it was ok. **EDIT:** I updated my javascript to the following:               rasterLayer = new OpenLayers.Layer.WMS(         "NY Population layer..... ",        "http://urbmet.webfactional.com/mapserv.cgi?",         {map: 'raster_1.map', layers: 'ny_pop'}     );      but still nothing shows up on the map. I've added a screencapture of my firebug output - this is the contents of `console.log(rasterLayer)`. ![enter image description here](http://i.stack.imgur.com/wVWMd.png)


title: Open source space time series database management system

text: I am aware of the existence of space time DBMSs such as Druid, InfluxDB or OpenTSDB. Is there anything like it for spatially explicit data? A time series extension for a common spatial DBMS (say, PostGIS) could do it; as well as a spatial extension to an existing time series DBMS.


title: How to select the cross streets in postgis spatial table?

text: I've used ST_Intersects and ST_Crosses but between two different tables. I have a table of streets and intersections want to know them. For example select all streets that cross the Street A. Can anyone give me an example? Edited for karpfen** This is my table of streets, each street segment is a geometry table. The image of the QGIS layer can better understand, analyze and think it would be appropriate to use ST_Touches and select geometry that touches the other at its end. But how do it in the same table?? I can not do a join. ![enter image description here](http://i.stack.imgur.com/Lr6m5.png)


title: Why are my symbols not displayed correctly in the legend and after a PDF-Export?

text: I work with QGIS 2.0.1. I have two problems with the representation of symbols. The first problem is that after I selected a symbol with the type: font and adjust the size of the symbol with the unit: map units, the symbol is not displayed completely within the Legend. Is this a general problem that the icon size is also scaled within the legend? ![SymbolOVERVIEW](http://i.stack.imgur.com/gkwYu.jpg) Legend: ![Legend](http://i.stack.imgur.com/9ySMH.jpg) * * * The next problem is that when I want to export my map to a PDF file the Symbols are not displayed correctly. The icon consists of two components. I have tried two different versions but the result was not useful. Here are the two variants i have tried:   1. Symbol as font and white background as a Font.    2. Symbol as font and white background as a simple symbol with a square. In QGIS both variants provide a correct representation. But the result was similar for both variants. The Symbol components are shifted from each other in the PDF-Export. Symbols in PDF-Export File: ![SymbolsPDF](http://i.stack.imgur.com/1rmAw.jpg) Symbols in QGIS ![SymbolsQGIS](http://i.stack.imgur.com/RShyz.jpg) * * * Thanks for your help.


title: Identifying X-Y coordinates from a raster

text: I need to find the high and low values of multiple raster datasets individually. Then determine the x,y location of these values. I am using ArcGIS 10.1, have built the mosaic datasets I need, and have determined the locations of the high and low values (estimations are fine). However I cannot determine the X,Y location for these values (points). I know there is a "Find x-y" tool, but I need the opposite. Is there a way to point click and return the x-y for that location? With vector data I could populate a field with x-y values, however I do not have the time to cut out the raster datasets individually and turn them into vector data to do this. There must be another way?


title: How to do the “display route events” from Arcgis, in Qgis

text: I'm working with Arcgis 10.1. In Arcgis there is a feature called "display route events" from tables, and I want to do exactly what that features does, but in Qgis. I'm using Qgis 2.0.1. What this feature does in arcgis is add information from an excel table to a route. In the table I have points with the information of there position but not the coordinates, I have the distance from the origin to the kilometric point where he is placed. Is there a way of doing it in Qgis? Or it just works with the coordinates? Thank you!


title: Update value in domain and same should reflect in attribute field

text: I have a SDE GDB, where i have no. of domains created. we have created domain for state names, which is used by all feature classes to update their respective state name. Now i have discovered that spelling / name of some of the states are not correct. if i correct these spellings or names in domain. will it reflect in the attribute of all feature classes or is there any other way forward? There are many thousand records in the GDB.


title: Identify task: Input polygon into IdentifyParameters

text: I wish to find all the features within a polygon/circle drawn by my user. May I know if I should use IdentifyTask https://developers.arcgis.com/en/javascript/jsapi/identifytask.html Or is there a more recommended way? Use Case: User drawn a circle/polygon and it will display all the address within the area.


title: What Spatial Analyst tool to use when comparing competitive species distributions?

text: _A spatial analysis question:_ I'm trying to compare species distributions and specifically **overlap** between two species _(Water Vole and European Rabbit)_ that are competitive species. I also want to look how the common shared predator the Minx causes this distribution. Can anyone suggest what ArcTools I can use?


title: How To Determine Percent of Matched Data?

text: In my attribute table, I have 2 columns with qualitative (categorical) data. Both of them are from different sources so I am trying to find the percent of matched data from these 2 columns. I created a new column labeled 'Matched', and wrote the following codeblock in the Field Calculator to determine the values, but I am getting an error on line 2-3. Do I have to convert them to string using `str()`?               def matchedRecords(x):         x = (column.name)         x2 = (column2.name)              if x == x2:             resultIs = "Matched"         else:             resultIs = "Unmatched"         return resultIs      And then `matchedRecords(!Matched!)` in the second area of the Field Calculator. I want to determine the total matched values versus unmatched values to find the percentage. Is this the best way to do it? Thanks.


title: Merge two polygons, attributes are added - Tool?

text: I have a bunch of polygons (counties) that I need to merge together. They all have various population attributes. I need to be able to selectively merge two counties together and get a single polygon with the combined population totals for the attributes. Is there a way to do this in ArcGIS? I've tried the Merge tool but while that works for the boundaries, it doesn't work for the adding the two population totals together, you have to pick one or the other. The Dissolve tool seems to operate on all polygons with a single attribute, but there is no unifying attribute that I have that determines which polygons are merged. Thanks.


title: Is r.cost equivalent to diffusion interpolation with barriers?

text: I want to perform interpolation using cost distances rather than euclidean distances similar to Greenberg et al. (2011). In that paper they use r.cost in GRASS to compute cost distances for each prediction location. Unfortunately I am more familiar with ArcGIS than GRASS. Can I use the diffusion/kernel interpolation with cost surfaces in ArcGIS rather than GRASS/r.cost? Would I get equivalent results? Greenberg, J.A., C. Rueda, E.L. Hestir, M.J. Santos and S.L. Ustin. 2011. Least cost distance analysis for spatial interpolation. Computers & Geosciences 37:272-276. http://dx.doi.org/10.1016%2Fj.cageo.2010.05.012


title: How do I make line symbols intersect cleanly without showing edges in QGIS?

text: I digitised a map showing electric grid lines. I merged a number of lines into a single feature. Now, I'd like these lines to be displayed so that the intersections appear continuous. At the moment, the line edges on intersections are visible and it's not pretty, as visible on the screenshot below: ![Overlapping line borders at intersections](http://i.stack.imgur.com/UlW8X.png) I could swear I've seen something about this somewhere on the Internet but I can't find it anymore.


title: passing geocode results object to another file

text: I'm using a geocoding plugin (Leafet.GeoSearch) on a map and I'd like the results, being the LatLng point coordinates, etc available as a parameter for GeoServer cql filtered calls. It doesn't seem to be as easy as I thought it would be. I've read several forums which claim data can't be persistent between pages/files. Could someone enlighten me and or make a suggestion please. I'm in over my head. here is a snippet of the code that has the array I'd like available               ParseJSON : function(data) {             if (data.resourceSets.length == 0 || data.resourceSets[0].resources.length == 0)                 return [];                  var results = [];             for (var i = 0; i < data.resourceSets[0].resources.length; i++)                 results.push(new L.GeoSearch.Result(                     data.resourceSets[0].resources[i].point.coordinates[1],                      data.resourceSets[0].resources[i].point.coordinates[0],                      data.resourceSets[0].resources[i].address.formattedAddress));             return results;         }


title: PDF to be understandable by ArcGIS

text: I have a few pdf files. There are vector graphics in the maps. How can I use them to be visible and editable in ArcGIS? I want to see the attribute tables content, etc.


title: I nead to find the nearest places based on latitude/longitude points

text: I have a java aplication backed by a mysql database which contains events with locations (latitude/longitude), the location data is saved as POINT data type. I nead to find the nearest events in some radius. I am already able to do this but the results are innacurate. I'm also using hibernate spatial to make the query. I'm using the folowing libraries : 1.) com.vividsolutions.jts ... for creating geometry types 2.) Geotools org.geotools ... for calculating distance between points First I make a circle Geometry with a center placed in a current location this circle is given some radius. For example 500m if I want to find the events in the distance of 500m. The problem is that the Geometry construct takes neads the radius in degrees. If i'm not mistaken the distance in degrees changes based on a point on earth. then I query the database with the within hibernate spatial function, but the results are unnacurate. So **I dont know how to calculate the distance in degrees based on location on earth.** I think there should already be some library that is doing this, but for now I havent find an example. Can you advise me some code? My Code:               //calculate the shape where the events should be      Geometry circle = createCircle(currentLongitude,currentLatitude, radiusInKm);     //retrieve the events from database that are within the calculated shape     List<si.arctur.model.Event> seznam = find(circle.toString());                   //the methods     private static Geometry createCircle(double x, double y, final double RADIUS) {                    double latitude = 110.54;               double longitude = 111.320;                    //TODO beter accuracy               double radiusInDegreesLatitude = RADIUS / latitude;               double radiusInDegreesLongitude = RADIUS / Math.abs((longitude * Math.cos(latitude)));                    GeometricShapeFactory shapeFactory = new GeometricShapeFactory();               shapeFactory.setNumPoints(1000);               shapeFactory.setHeight(radiusInDegreesLatitude * 2);               shapeFactory.setWidth(radiusInDegreesLongitude * 2);               shapeFactory.setCentre(new Coordinate(x, y));//there are your coordinates               Polygon shape = shapeFactory.createCircle();               return shape;          }              private List find(String wktFilter) {             Geometry filter = wktToGeometry(wktFilter);             Query query = em.createQuery("select e from Event e where within(e.location, :filter) = true", si.arctur.model.Event.class);             query.setParameter("filter", filter);             return query.getResultList();         }              private Geometry wktToGeometry(String wktPoint) {             WKTReader fromText = new WKTReader();             Geometry geom = null;             try {                 geom = fromText.read(wktPoint);             } catch (ParseException e) {                 throw new RuntimeException("Not a WKT string:" + wktPoint);             }             return geom;         }


title: Can ArcPublisher 10 publish backwards to an 9.3 PMF?

text: It seems like a lot of other functions within ArcGIS are getting backwards compatibility (file geodb creation, MXD files). I'll still need to publish 9.3 PMF files for a few months and I couldn't find anything in the Publisher settings that would allow me to publish to a previous version. Is it possible?


title: How do I open the Identify window in ArcMap?

text: My identify window no longer appears when I identify a feature in ArcMap.


title: Checkboxes turn off layers within a mapservice, and I want them to make the layers visible?

text: This webmap project I keep working accesses mutilple layers within a couple different mapservices. I would like to learn how I can load the layers without them drawing. I've recently incorporated check boxes instead of links to turn each layer on and off. They work but in the reverse way I want them to. I've not been able to figure out how to keep the layers from drawing until I select a checkbox.. Any ideas or help would be appreciated.. Thanks.. here is a link.. http://www.co.frederick.va.us/planning/newmap1copy2.html


title: Where can I get ArcView 1.0?

text: The version is correct 1.0. Not 10.0, 1.0. Will it run on Windows 7?


title: Export raster file geodatabase to jp2 or jpeg?

text: Can I export or copy a raster from a file Geodatabase and create a output to jpeg or jp2000 formate? I was reading two different helpline on ArcGIS 10.1 and wasn't sure if I was supposed to export jp2000 from another file geodatabase to file geodatabase. I tried to do that the other day but it failed. Don't know if it was the space was issued or maybe I wasn't doing the proper way. Thanks.


title: Why can't I get MapServer to show a raster?

text: I am trying to use MapServer to display a raster image (png). That's it for now at least. However, I am having no luck. When I visit my map by using a link like: *http://myurl.com/cgi-bin/mapserv?map=/path/to/file.map&layer=base&mode=map I get the following error: > msLoadMap(): Regular expression error. MS_DEFAULT_MAPFILE_PATTERN validation > failed. msEvalRegex(): Regular expression error. String failed expression > test. When I use shp2img:               shp2img -m file.map -o myimage.png -all_debug 5      I get no errors, but the image that's produced is completely gray except for the upper right quadrant which is white. This is nothing like my raster image. My mapfile:               MAP         SHAPEPATH     "/home/name/public_html/dir"         IMAGETYPE     PNG         EXTENT        -180 -90 180 90         SIZE          800 600         IMAGECOLOR    155 155 155         UNITS         DD         STATUS        ON         DEBUG         ON              LAYER             NAME      "base"             DATA      "/home/name/public_html/dir/image.png"             STATUS    ON             TYPE      RASTER         END      END       What am I doing wrong?


title: How to call gdal_retile from Python code?

text: I am trying to call gdal_retile from my Python code like this:               import gdal_retile       and in main:               gdal_retile.main("-v -r bilinear -levels 4 -ps 2048 2048 -co \"tiled=YES\" -targetDir pyramid --optfile files.txt")      but I get this error in my Debug Probe window:               Unrecognised command option: -     Usage: gdal_retile.py          [-v] [-co NAME=VALUE]* [-of out_format]         [-ps pixelWidth pixelHeight]         [-ot  {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/                CInt16/CInt32/CFloat32/CFloat64}]         [ -tileIndex tileIndexName [-tileIndexField fieldName]]         [ -csv fileName [-csvDelim delimiter]]         [-s_srs srs_def]  [-pyramidOnly] -levels numberoflevels         [-r {near/bilinear/cubic/cubicspline/lanczos}]         [-useDirForEachRow]         -targetDir TileDirectory input_files      Any idea what I am doing wrong?


title: Detrending DEMs? - Removing the elevation trend over sloped surfaces (diagram included)

text: Is there a known methodology for removing the trend from a DEM surface? For example, if I had a DEM of a river, that river would flow from high elevation (upstream) to lower elevations (downstream). Even if the distance from the top of the water surface to the bottom of the river bed were the same across the whole river, the downstream portion of the raster would still have lower elevation values in the DEM. Is it possible to even out these values so that areas of equal depth also possess equal elevations? When removing this trend it would be ideal to preserve the overall characteristics of the DEM, such as dips in the riverbed.![enter image description here](http://i.stack.imgur.com/HDtnp.png)


title: Selecting POIs around specific buildings using PostGIS

text: The problem: I need to select, for each residential building in my table that has say at least 2 pharmacies and 2 education centers within a radius of 1km, all POIs (pharmacies, comercial centres, medical centers, education centers, police stations, fire stations) which are within 1km of the respective building. table structure-> building ( **id** serial, **name** varchar ) poi_category( **id** serial, **cname** varchar) --cname being the category name of course poi( **id** serial, **name** varchar, **c_id** integer)-- c_id is the FK referencing poi_category( **id** ) all coordinate columns are of type geometry not geography (let's call them **geom** ) here's the way i thought it should be done but i'm not sure it's even correct let alone the optimal solution to this problem                   SELECT r.id_b, r.id_p         FROM (              SELECT b.id AS id_b, p.id AS id_p, pc.id AS id_pc,pc.cname              FROM building AS b, poi AS p, poi_category AS pc              WHERE ST_DWithin(b.geom,p.geom, 1000) AND p.c_id=pc.id              ) AS r,              (              SELECT * FROM r GROUP BY id_b              ) AS r1               HAVING  count (                            SELECT *                            FROM r, r1                            WHERE r1.id_b=r.id_b AND r.id_pc='pharmacy'                              )>1                      AND                      count (                            SELECT *                            FROM r, r1                            WHERE r1.id_b=r.id_b AND r.id_pc='ed. centre'                              )>1      Is this the way to go for what i need ? What solution would be better from a performance point of view? What about the most elegant solution?


title: Using the API of QSpatiaLite plugin

text: I am using the popular QSpatiaLite plugin for QGIS. Since I am creating and modifying SpatiaLite tables directly with Python (using `pyspatialite`), I wonder whether I can automate tasks such as converting from a Shapefile to SpatiaLite directly from the Python code (similar to the plugin's "Import OGR files"). Thus, does the plugin provide some sort of an API to do that kind of stuff? Or can I do this from within `pyspatialite` without using any QGIS-related API? I know there is a QGIS Python API but I have never dealt with it, that's why I am kind of confused here.


title: What is some theory on the number of GCP when georeferencing?

text: Just wondering, is it always better to have more than 3 GCPs? Should the GCPs be spread out evenly? Where can I find some general resources on the theory of this? Other than trial an error, which simply teaches me what works at that moment in time and not what should work in given instances. Not critical, but just interested.


title: how to use ArcGIS cloud hosted map with custom data from mysql

text: I have map hosted on ArcGIS.com which I can use and embed in my webpage. I have few tables with latitude-longitude in mysql. is there anyway I can add it as marker overlay layer on map taken from ArcGIS.com ? Note: My server does not support spatial tools like PostgreSQL and Geoserver. SO I am trying to find an alternative. any help and guidance will be appreciated.


title: Is there a way to improve the render quality of ESRI basemaps when exporting a map to PDF or JPG?

text: I've had this problem with ESRI base maps before but not this bad. The problem is that when I use an ESRI basemap as my backdrop, it looks fine in my MXD but when I export the map to PDF, JPG, or whatever format, the quality of the basemap in the export file is a lot worse than what I see on-screen in ArcGIS. Here's the most recent, egregious example: What I see in ArcMap: ![enter image description here](http://i.stack.imgur.com/GjEyw.png) What it looks like in my exported PDF: ![enter image description here](http://i.stack.imgur.com/KtxxG.png) I've tried several things to solve this problem, and the result you see above is the best I can do:   * Turning up the export quality on the PDF to 600dpi or higher   * Exporting to JPG, or to TIF with no compression at all   * Clearing my cache on the basemap layers I'm using   * Printing to the Adobe PDF print driver rather than using Export to PDF Some people apparently have this problem due to a slow internet connection, but since I can see the tiles in my MXD and they look fine, I don't think that's my problem. Is there something I'm missing?


title: How tool to split polylines into two-part lines at vertices in ArcGIS for Desktop 10.1?

text: I'm using ArcGIS 10.1 SP1 at and we have a problem. There is a serious topology bug that I've mentioned in another question that is affecting our workflow. At present, Esri has stated the fix will be in 10.2, but my government is not able to upgrade to 10.2 so easily. Polylines flagged by topology as overlapping or not being covered by something are normally easily subtracted. With the bug, the subtraction affect the entire polyline, beyond the topologically incorrect area. My workaround is to use the `Split Into COGO Lines` tool, which is on the COGO toolbar, to split the offending polyline at all nodes, including where the subtraction should end (checking for topology errors generates nodes). Then, when I do the subtract, it only wipes out the offending lines. However, it also adds COGO attributes to all the other segments along the former polyline. This attribute update behavior is undesirable. Is there a simple tool that would split a line at all nodes or vertices indiscriminately--and either leave COGO attributes in place, or wipe them?


title: Merging Personal Geodatabase (MDB) data into ArcSDE using FME?

text: Im a novice in FME and I have a MDB file that i want to write/insert to an SDE database using FME. I created a simple translation that generally reads the data from the MDB file and writes it to the SDE database, although the translation is successful and it prints out the number of written rows, when i checked the SDE database, it wasn't able to right any data on it. i then checked the log and i found this specific warning: > MS SQL Server (Spatial) Writer: Attribute data `' too large to fit in > specified data type for column`SHAPE' The shape column holds the geometry data in the MDB file. FME sets the SHAPE datatype to 'oleobject' in the Reader and I noticed that the SHAPE column in the writer was automatically converted to 'text' datatype. Is there anything else i need to configure to make this work? **UPDATE** : I am using FME Database Edition 2014 (20140120 - Build 14235 - WIN64) Here are some screenshots of my FME translation: ![enter image description here](http://i.stack.imgur.com/DZbME.png) ![enter image description here](http://i.stack.imgur.com/KfZRT.png)![enter image description here](http://i.stack.imgur.com/HGazD.png)


title: How to get RasterToPolygon working

text: We have a process that is almost working except for this one call ... This is the one tool function that is requiring a FeatureClass and I found code that I thought would work. Unfortunately all I get is an HRESULT of Fail. So here is what I came up with so far: I have my raster. I found code for making an in_memory featureclass but the call still fails! Here is my code so far... What did I mess up? (By the way...my RunTool method has been used numerous times and is not the source of the issue!)                ESRI.ArcGIS.ConversionTools.RasterToPolygon rtpConverter = new RasterToPolygon();      Debug.WriteLine("Getting the blank feature class.");      IObjectClassDescription objectClassDescription = new FeatureClassDescriptionClass();      IFields fields = objectClassDescription.RequiredFields;      UID ClsID = new UIDClass();      ClsID.Value = "esriGeodatabase.Feature";      IFeatureWorkspace workSpace = (IFeatureWorkspace)GetInMemoryWorkspace();      IFeatureClass donut = workSpace.CreateFeatureClass("WorkspaceTempName", fields, ClsID, null, esriFeatureType.esriFTSimple, "SHAPE", string.Empty);      rtpConverter.in_raster = delayZone;      rtpConverter.out_polygon_features = donut;      Debug.WriteLine("Getting the raster to polygon conversion.");      RunTool(rtpConverter, true);      IPolygon polygonDonut = donut.GetFeature(1).Shape as IPolygon;      The following is what is in GetInMemoryWorkspace()               Type factoryType = Type.GetTypeFromProgID("esriDataSourcesGDB.InMemoryWorkspaceFactory");     IWorkspaceFactory workspaceFactory = (IWorkspaceFactory)Activator.CreateInstance(factoryType);     IWorkspaceName workspaceName = workspaceFactory.Create("", "MyWorkspace", null, 0);     IName name = (IName)workspaceName;     return (IWorkspace)name.Open();


title: WKBreader reuse - jts 1.9

text: I'm having some bad behavior reusing wkbreader. The context as of a single threaded as I use threadlocal wkbreader. I get parse exception on some geoms. When I create a new wkbreader for each read - it works flawlessly. Any ideas? Thanks a lot!


title: Updating column in Field Calculator using values of another column?

text: I am very new to GIS, so I hope to get some help. I am using ArcGIS 10.2 for Desktop. I want to create a new numeric column on my attribute table to store "seasons" in numbers. For example, I want the "seasons" field to hold a numeric code (i.e. "summer" will be "1", "autumn" to be "2", "winter" to be "3" and "spring" to be "4"). How can I do this using the field calculator?


title: Complex ArcGIS Server Infra problem

text: We have a large application running on Java, RHEL, where our ArcGIS Server resides also. We have a lot of inteligence on some of our classes and we wish to publish them via arcgis server. The problem is, that being RHEL I cannot just register the class extension dlls in it and when AGS tries to initialize them, it fails. I don't need the class extensions to work via web, I just need to be able to show them up on a map (no editions will be made). Is there any viable solution that does not involve: A) recoding; B) switching to a different environment; Can I figure out a way of skipping the initialization of the class extension on the server? EDIT: Since I don't need the complex edit behavior we have in ArcGIS Desktop, I'm gonna try to fool AGS. I'm going to build the same class extension as we have on the .NET side, but with Java, without behavior and deploy it to the ags.


title: Combined WMS request issue (one layer has black background)

text: We are using MATLAB to request a combined WMS from GeoServer. The request is auto-generated as this:               http://geo-pprod.test.cresis.ku.edu/geoserver/arctic/wms?SERVICE=WMS&LAYERS=greenland_natural_90m,rds_line_paths&CRS=EPSG:3413&FORMAT=image/geotiff8&HEIGHT=527&BGCOLOR=0xFFFFFF&REQUEST=GetMap&WIDTH=620&BBOX=-1465312.4055326465,-3439603.5364015717,1669546.4179967656,-774973.5364015717&STYLES=,&VERSION=1.3.0      It's requesting a GeoTIFF store with a linestring (PostGIS store) both published in EPSG:3413. **greenland_natural_90m:** Greenland GeoTIFF (declared native as EPSG3413) **rds_line_paths:** PostGIS linePath geometry table (re-projected from WGS84 to EPSG3413) The result from this request looks like this: http://i.imgur.com/FSC6fcQ.png However when I load the same two layers in our OpenLayers web client I get the following (correct): http://i.imgur.com/rsr8XTF.png Note that we have tried every image format and get the same results. We have also tested on multiple versions of GeoServer. Any help would be appreciated. Unfortunatly our server in not public yet so please keep that in mind if additional information is needed. What we have tried:   * View in OpenLayers (works)   * Get each layer individually (works) **We really need it combined   * Get layers (combined/individual) in multiple formats (Works individual, fails combined in MATLAB)   * Revert to old versions of GeoServer (this has worked before, but continues to fail now) Notes: We are running Geoserver 2.3.2 WAR served via Tomcat6 through Apache on an RHEL linux server.


I have a problem with layers and variables box in raster calculator. The layers and variables box does not exist in the raster calculator window! So I can not add any raster to the expression. I am using Arc GIS 10. what should I exactly do to appear this box in the raster calculator? thanks for your help


title: Error running geotools tutorial for CSV2SHP?

text: Although I run this tutorial I always get following error message: > Nov 22, 2011 5:10:50 PM > org.geotools.referencing.factory.epsg.ThreadedPostgreSQLEpsgFactory > isAvailable WARNING: Unavailable authority factory: European Petroleum > Survey Group org.opengis.referencing.FactoryException: Failed to connect to > the EPSG database. at > org.geotools.referencing.factory.epsg.ThreadedEpsgFactory.createBackingStore(ThreadedEpsgFactory.java:420) > at > org.geotools.referencing.factory.DeferredAuthorityFactory.getBackingStore(DeferredAuthorityFactory.java:131) > at > org.geotools.referencing.factory.BufferedAuthorityFactory.isAvailable(BufferedAuthorityFactory.java:233) > at > org.geotools.referencing.factory.DeferredAuthorityFactory.isAvailable(DeferredAuthorityFactory.java:117) > at > org.geotools.factory.FactoryRegistry.isAvailable(FactoryRegistry.java:665) > at > org.geotools.factory.FactoryRegistry.isAcceptable(FactoryRegistry.java:499) > at org.geotools.factory.FactoryRegistry$1.filter(FactoryRegistry.java:190) > at javax.imageio.spi.FilterIterator.advance(ServiceRegistry.java:810) at > javax.imageio.spi.FilterIterator.(ServiceRegistry.java:804) at > javax.imageio.spi.ServiceRegistry.getServiceProviders(ServiceRegistry.java:508) > at > org.geotools.factory.FactoryRegistry.getServiceProviders(FactoryRegistry.java:195) > at > org.geotools.referencing.ReferencingFactoryFinder.getFactories(ReferencingFactoryFinder.java:176) > at > org.geotools.referencing.ReferencingFactoryFinder.getCRSAuthorityFactories(ReferencingFactoryFinder.java:451) > at > org.geotools.referencing.DefaultAuthorityFactory.getBackingFactory(DefaultAuthorityFactory.java:88) > at > org.geotools.referencing.DefaultAuthorityFactory.(DefaultAuthorityFactory.java:69) > at org.geotools.referencing.CRS.getAuthorityFactory(CRS.java:229) at > org.geotools.referencing.CRS.decode(CRS.java:467) at > org.geotools.referencing.CRS.decode(CRS.java:395) at > org.geotools.data.DataUtilities.createAttribute(DataUtilities.java:2172) at > org.geotools.data.DataUtilities.createType(DataUtilities.java:1752) at > org.geotools.data.DataUtilities.createType(DataUtilities.java:1702) at > org.geotools.tutorial.Csv2Shape.main(Csv2Shape.java:61) Caused by: > org.postgresql.util.PSQLException: FATAL: password authentication failed for > user "Geotools" at > org.postgresql.core.v3.ConnectionFactoryImpl.doAuthentication(ConnectionFactoryImpl.java:291) > at > org.postgresql.core.v3.ConnectionFactoryImpl.openConnectionImpl(ConnectionFactoryImpl.java:108) > at > org.postgresql.core.ConnectionFactory.openConnection(ConnectionFactory.java:66) > at > org.postgresql.jdbc2.AbstractJdbc2Connection.(AbstractJdbc2Connection.java:125) > at > org.postgresql.jdbc3.AbstractJdbc3Connection.(AbstractJdbc3Connection.java:30) > at org.postgresql.jdbc3.Jdbc3Connection.(Jdbc3Connection.java:24) at > org.postgresql.Driver.makeConnection(Driver.java:393) at > org.postgresql.Driver.connect(Driver.java:267) at > java.sql.DriverManager.getConnection(DriverManager.java:579) at > java.sql.DriverManager.getConnection(DriverManager.java:221) at > org.postgresql.ds.common.BaseDataSource.getConnection(BaseDataSource.java:87) > at > org.postgresql.ds.common.BaseDataSource.getConnection(BaseDataSource.java:70) > at > org.geotools.referencing.factory.epsg.DirectEpsgFactory.getConnection(DirectEpsgFactory.java:3172) > at > org.geotools.referencing.factory.epsg.ThreadedEpsgFactory.createBackingStore(ThreadedEpsgFactory.java:415) > ... 21 more > > Nov 22, 2011 5:10:52 PM > org.geotools.referencing.factory.epsg.ThreadedPostgreSQLEpsgFactory > isAvailable WARNING: Unavailable authority factory: European Petroleum > Survey Group org.opengis.referencing.FactoryException: Failed to connect to > the EPSG database. at > org.geotools.referencing.factory.epsg.ThreadedEpsgFactory.createBackingStore(ThreadedEpsgFactory.java:420) > at > org.geotools.referencing.factory.DeferredAuthorityFactory.getBackingStore(DeferredAuthorityFactory.java:131) > at > org.geotools.referencing.factory.BufferedAuthorityFactory.isAvailable(BufferedAuthorityFactory.java:233) > at > org.geotools.referencing.factory.DeferredAuthorityFactory.isAvailable(DeferredAuthorityFactory.java:117) > at > org.geotools.factory.FactoryRegistry.isAvailable(FactoryRegistry.java:665) > at > org.geotools.factory.FactoryRegistry.isAcceptable(FactoryRegistry.java:499) > at > org.geotools.factory.FactoryRegistry.getServiceImplementation(FactoryRegistry.java:435) > at > org.geotools.factory.FactoryRegistry.getServiceProvider(FactoryRegistry.java:363) > at > org.geotools.factory.FactoryCreator.getServiceProvider(FactoryCreator.java:143) > at > org.geotools.referencing.ReferencingFactoryFinder.getAuthorityFactory(ReferencingFactoryFinder.java:216) > at > org.geotools.referencing.ReferencingFactoryFinder.getCRSAuthorityFactory(ReferencingFactoryFinder.java:436) > at > org.geotools.referencing.factory.epsg.LongitudeFirstFactory.createBackingStore(LongitudeFirstFactory.java:189) > at > org.geotools.referencing.factory.DeferredAuthorityFactory.getBackingStore(DeferredAuthorityFactory.java:131) > at > org.geotools.referencing.factory.BufferedAuthorityFactory.isAvailable(BufferedAuthorityFactory.java:233) > at > org.geotools.referencing.factory.DeferredAuthorityFactory.isAvailable(DeferredAuthorityFactory.java:117) > at > org.geotools.factory.FactoryRegistry.isAvailable(FactoryRegistry.java:665) > at > org.geotools.factory.FactoryRegistry.isAcceptable(FactoryRegistry.java:499) > at org.geotools.factory.FactoryRegistry$1.filter(FactoryRegistry.java:190) > at javax.imageio.spi.FilterIterator.advance(ServiceRegistry.java:810) at > javax.imageio.spi.FilterIterator.next(ServiceRegistry.java:828) at > org.geotools.resources.LazySet.addNext(LazySet.java:70) at > org.geotools.resources.LazySet.get(LazySet.java:120) at > org.geotools.resources.LazySet$Iter.next(LazySet.java:139) at > org.geotools.referencing.factory.ManyAuthoritiesFactory.(ManyAuthoritiesFactory.java:118) > at > org.geotools.referencing.DefaultAuthorityFactory.getBackingFactory(DefaultAuthorityFactory.java:116) > at > org.geotools.referencing.DefaultAuthorityFactory.(DefaultAuthorityFactory.java:69) > at org.geotools.referencing.CRS.getAuthorityFactory(CRS.java:229) at > org.geotools.referencing.CRS.decode(CRS.java:467) at > org.geotools.referencing.CRS.decode(CRS.java:395) at > org.geotools.data.DataUtilities.createAttribute(DataUtilities.java:2172) at > org.geotools.data.DataUtilities.createType(DataUtilities.java:1752) at > org.geotools.data.DataUtilities.createType(DataUtilities.java:1702) at > org.geotools.tutorial.Csv2Shape.main(Csv2Shape.java:61) Caused by: > org.postgresql.util.PSQLException: FATAL: password authentication failed for > user "Geotools" at > org.postgresql.core.v3.ConnectionFactoryImpl.doAuthentication(ConnectionFactoryImpl.java:291) > at > org.postgresql.core.v3.ConnectionFactoryImpl.openConnectionImpl(ConnectionFactoryImpl.java:108) > at > org.postgresql.core.ConnectionFactory.openConnection(ConnectionFactory.java:66) > at > org.postgresql.jdbc2.AbstractJdbc2Connection.(AbstractJdbc2Connection.java:125) > at > org.postgresql.jdbc3.AbstractJdbc3Connection.(AbstractJdbc3Connection.java:30) > at org.postgresql.jdbc3.Jdbc3Connection.(Jdbc3Connection.java:24) at > org.postgresql.Driver.makeConnection(Driver.java:393) at > org.postgresql.Driver.connect(Driver.java:267) at > java.sql.DriverManager.getConnection(DriverManager.java:579) at > java.sql.DriverManager.getConnection(DriverManager.java:221) at > org.postgresql.ds.common.BaseDataSource.getConnection(BaseDataSource.java:87) > at > org.postgresql.ds.common.BaseDataSource.getConnection(BaseDataSource.java:70) > at > org.geotools.referencing.factory.epsg.DirectEpsgFactory.getConnection(DirectEpsgFactory.java:3172) > at > org.geotools.referencing.factory.epsg.ThreadedEpsgFactory.createBackingStore(ThreadedEpsgFactory.java:415) > ... 32 more > > Nov 22, 2011 5:10:58 PM > org.geotools.referencing.factory.PropertyAuthorityFactory WARNING: Axis > elements found in a wkt definition, the force longitude first axis order > hint might not be respected: PROJCS["Google Projection",GEOGCS["WGS > 84",DATUM["World Geodetic System 1984",SPHEROID["WGS 84", 6378137.0, > 298.257223563, > AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich", 0.0, > AUTHORITY["EPSG","8901"]],UNIT["degree", > 0.017453292519943295],AXIS["Geodetic longitude", EAST],AXIS["Geodetic > latitude", NORTH],AUTHORITY["EPSG","4326"]],PROJECTION["Popular > Visualisation Pseudo Mercator", > AUTHORITY["EPSG","1024"]],PARAMETER["semi_minor", > 6378137.0],PARAMETER["latitude_of_origin", > 0.0],PARAMETER["central_meridian", 0.0],PARAMETER["scale_factor", > 1.0],PARAMETER["false_easting", 0.0],PARAMETER["false_northing", > 0.0],UNIT["m", 1.0]] What might be the cause that application throws this error I have not used any database here?


title: Print axMapControl contents with preview

text: I have a map with many feature layers on top of it. I would like to be able to print the map and feature layers either scaling it to fit a certain page size or printing it to actual scale thus creating a mosaic of printed pages. Ideally I would like to preview the results before sending them to printer and I am having two major issues doing so. What I did so far is get the map onto a axPageLayoutControl so that I can add artifacts about which part of the mosaic each page is, scalebars, titles etc. Although getting the axMapControl.ActiveView.FocusMap which appears as landscape, in the PageLayout (orientation set to landscape) the map appears in portait (actually showing a slightly different extend). Changing the pageLayout.page.orientation has no effect on the map. Also I tried drawing all the features from the axMapControl onto the pageLayout but the spread in a much bigger area than the map itself, but this time in landscape orientation. Another problem is that I don't know how to preview the mosaic pages after setting the esriPageToPrinterMapping.esriPageMappingTile value to the axPageLayoutControl.Page.PageToPrinterMapping property as this does change the number of pages axPageLayoutControl.get_PrinterPageCount() to more than 1, but has no effect on the pageLayout or the printout itself.