texts: Using the PgRouting driving_distance function can we see the edges?

texts:
I have a directed graph and am using PgRouting's  driving_distance function.  My application needs it (or a related PostreSQL procedure) to return a full set of EDGES emanating out of the source, and not the just the VERTICES.

The Drivetime function (the essential "Dijkstra algorithm") returns the following columns:

Vertex_id
CostToSource (from that vertex)
Edge_id (where edge.source = vertex_id )

Of course a vertex can have several edges associated with it and only one edge per vertex is returned in this table (for some reason).

Qu: Can anybody think of a way to relate the above information with the network so that I can get what I need, which is:

Related Pair ( Vertex_ID, Edge_ID)
CostToSource (From vertex_id)

Hint: the obvious problem with joining vertex_id to the "source" in the network table is that some of the links do not belong in the result set. Any link heading back to the source should be excluded.

Thanks all.


texts: How to get layer variables using QGIS Python APIs?

texts: How can I access the layer variables defined by the user from QGIS Python APIs?


texts:
Is there a shortcut key to increase or decrease layer transparency, or to add the opacity slider to a layer in QGIS?
I couldn't find one in the existing shortcut key list and I tried loading the transparency file to that list, but it doesn't like the format.
A while back I noticed there Python script that can do something similar, add transparency to new layers imported into the workspace, but I am looking for a simpler tool.


texts: how to use "components" argument in geopy GoogleV3 geocoder?

texts:
I would like to use GeoPy to geocode a number of locations; I would like to use the GoogleV3 geocoder as it is supposed to allow me to restrict the results so that, for instance, I can specify that I only want to look for Springfield in Illinois and not in any other state. I think the "components" argument should allow me to do this.

However, I have not been able to figure out from the documentation what information I need to provide or how. All I have found is:
"components (dict) –
Restricts to an area. Can use any combination of: route, locality, administrative_area, postal_code, country."

Does anyone have an example of how this dictionary should be composed and of how it is called?


texts: Ending Python if-then using triple double quote gives EOL while scanning string literal

texts:
I keep getting the same error, I think it's because of how I end my if-then statement. It basically looks like this, except it's much longer (like, 300 lines long).
myIfThenVariable="""def myIfThen(field, fieldB):
   if field=="Thing One":
      return "T1"
   elif field=="Thing Two":
      return "T2"
   elif field=="Something else":
      return fieldB
   else:
      return field """"

The error I keep getting is:
SyntaxError: EOL while scanning string literal (MyRhubarb.py, line 374)

I have my if-then book-ended with triple quotes.
I tried adding a backslash after the word field ( field\""").
I copied and pasted the same code into the field calculator in ArcMap and was able to run it without issue.


texts: What is the W_{variablename} variable in GeoDa's spatial lag model?

texts:
What is the W_{variablename} variable in GeoDa's spatial lag model output?

The geoda workbook does not give any explanation.
As far as my knowledge goes a spatial lag regression model does not need an additional independent variable.


texts: What's the default histogram manipulation of QGIS for displaying RGB images?

texts: When you open an RGB geotiff (e.g. bands 4-3-2 of Sentinel-2) as raster dataset in QGIS - what histogram manipulation is applied by default? Of course I can look at the min-/max-values, but I'd like to know the underlying formulas/algorithm.


texts: If-then expression and font color in Mapinfo label expressions

texts:
How might I obtain the following, other than having two identical tables open to achieve: I would like to label an object using data from two columns from one table?

The first half of the label (from layer EN_FALL15:) will always be yellow bold font with black halo.

The second half of the label (from layer Cat_Priority:) would depend on the data in the layer:

When BLUE is entered, I would like blue bold font with white halo;

When RED is entered, I would like red bold font with white halo;

When YELLOW is entered, I would like yellow bold font with white halo.

Can mapinfo distinguish this difference and label accordingly?


texts:
In rLiDAR, you can segment trees by using CHMs ASCIIs produced by Fusion. Can the same be done in lidR? I know in lidR you have to use the .laz file for tree segmentation:
LASfile <- system.file("extdata", "MixedConifer.laz", package="lidR")
las <- readLAS(LASfile, select = "xyz", filter = "-drop_z_below 0")
# Using Li et al. (2012)
las <- segment_trees(las, li2012(R = 3, speed_up = 5))

So basically, instead of using the .laz file, can I just use its representative CHM.asc produced by Fusion to segment/find trees and then make convex hulls. The reason for using lidR is that I have done crown delineation using rLiDAR, but since rLiDAR does not use more than one CPU cores, it takes a lot of time. However, lidR does support parallelization, hence it's way faster than rLiDAR.


texts: create image out of geojson (gdal_rasterize?)

texts:
Short version: I need to know how I can convert a .geojson-file to an according .png-file (or any other image format)

long version:
I have a .geojson which was created by gdal_polygonize from a .png . I want to make that geojson to a .png again.
I tried gdal_rasterize but had an error:

C:\>gdal_rasterize source.geojson destination.png
At least one of -3d, -burn or -a required.

Usage: gdal_rasterize [-b band]* [-i] [-at]
       [-burn value]* | [-a attribute_name] [-3d]
       [-l layername]* [-where expression] [-sql select_statement]
       [-of format] [-a_srs srs_def] [-co "NAME=VALUE"]*
       [-a_nodata value] [-init value]*
       [-te xmin ymin xmax ymax] [-tr xres yres] [-tap] [-ts width height]
       [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/
             CInt16/CInt32/CFloat32/CFloat64}] [-q]
       <src_datasource> <dst_filename>

I looked here for more information but I don't understand these options. I just want a PNG file. It does not have to contain any information about the SRS and so on. The source.geojson contains many polygons and each polygon has an additional field called 'DN' in which is an HEX-Code of the color of that polygon.

Or do you know another simple method to create an image from a geojson-file?

Edit:
As user30184 described in the comments I tried this:

C:\>gdal_rasterize -burn 255 -burn 0 -burn 0 -ts 4950 4090 source.geojson target.tif

And it crated a picture (which is 463 MB..). I could not view it with the Windows Photo Viewer nor with IrfanView but I could inport it in QGIS as a raster file. But the colours are not correct. All polygons have the same color(red) and the "background" is black. I need the colors of each polygon in the image so that I can compare the original image and this created image.

This is the structure of my source.geojson:

{
"type": "FeatureCollection",
"crs": { "type": "name", "properties": { "name": "urn:ogc:def:crs:OGC:1.3:CRS84" } },

       "features": [
            { "type": "Feature", "properties": { "DN": "#FFFFA5" }, "geometry": { "type": "Polygon", "coordinates": [ [ [ 13.28, 52.66 ], [ 13.29, 52.65 ], [ 13.29, 52.66 ], [ 13.28, 52.66 ], [ 13.28, 52.66 ] ] ] } },
            { "type": "Feature", "properties": { "DN": "#FFFFAA" }, "geometry": {    "type": "Polygon", "coordinates": [ [ [ 13.39, 52.64 ], [ 13.39, 52.64 ], [ 13.397, 52.648 ], [ 13.39, 52.65 ], [ 13.39, 52.64 ] ] ] } }
       ]
}

You can see the colour in the 'DN'-field


texts: How can I upload an annotation featureclass to ArcSDE using C#?

texts: I can upload the shape files to ArcSDE database successfully,but I fail to upload annotation featureclass.


texts: Change from LineString to MultiLineString when saving layer edits in QGIS

texts:
In QGIS, when a user creates a new single part line in an open editing session for a “Multi-Line” SpatiaLite dataset, the QGIS data type is initially a “LineString”.  But, when the user commits the new record to the Multi-Line QGIS data set by pressing the “Save Layer Edits” button, QGIS changes the parsed data type from “LineString” to “MultiLineString”.  A similar change from “LineString” to “MultiLineString” occurs when committing a new single part line record to an ESRI Shapefile of type “PolyLine”.

Is the change from “LineString” to “MultiLineString” controlled by QGIS (e.g. try saving as “LineString”, except save as “MultiLineString”), or, is it controlled by relevant data provider?  If it is controlled by QGIS, does that imply it will be default logic for any line or multiline dataset?


texts:
I created some vector polygons, added there field height, fill it with values. After it using polygon to raster converted it to polygon (specified that value field to height).

Than I got rasters.

But when I called Viewshed. I got error: All cells in raster have the NODATA value.


texts: Masking out places with elevation below 5 m from a DEM

texts: I have a DEM of a study region; I want to modify the DEM in order to reflect the coastline outlook when the sea level was 5 meters higher than today. What I have done so far is to use 'Raster Calculator' with the quite simple following formula: DEM >= 5. Then, in the returned raster, I have set to 'hollow' the 1 values (representing cells with elevation >= 5m), and to 'blue' the 0 values. How and by which tools can I use the obtained raster as a mask to modify the original DEM. In other words, I want to mask out (i.e., remove) places with elevation below 5 m from the original DEM. I am using arcGIS 10.1 with advanced license.


texts: Creating geom object from list of lat, long coordinates


texts: Left to right labels disappearing in ArcMap

texts:
I try to write this word 'אזוה"ת' (Left to right label) with the value from field "nam" in the Label Expression with no success (the label expression dissapear):
u'אזוה"ת ' + [nam]

the wanted result should be: אזוה"ת כפר קאסם
this is my field details:

when I wrote the expression this way:
u'אזוהת ' + [nam]

it works but not as I need (need to add " Apostrophes sign):

I know I can create another field with the value 'אזוה"ת' , but i need the solution with the way I wrote above
I also read Integrate right to left text label expression with fields in ArcGIS 10.3  and didn't found any answer.


texts: What does this error mean- 'Underlying Connection was closed'

texts:
I am receiving this error message on one of our server while publishing a mobile project in mobile project center to ArcGIS server (10.1).

What connection is this one referring to and how to resolve this?


texts:
I have a KML file with Waypoints and several Tracks in it. When I open the file with QGIS, I see the LineString for each Track as well as a Points array for each Track. My problem is: QGIS shows the same Points array for each Track, even though there are discrete Points arrays for each Track in the KML file.

My process for opening the file was to start QGIS, click on "Add Vector Layer" and choosing "Select All" in the dialog that appeared listing each Folder in the KML file. I have also tried opening the file and importing only a single Points array instead of everything; I still get only the first Points array specified for the first Track in the KML file instead of the Points array for a specific Track.


texts: Estimating slope by fitting a quadratic surface (Evans 1980) in ArcGIS Pro

texts:
I am needing to calculate slope using the quadratic surface method of Evans (An integrated system of terrain analysis and slope mapping, 1980), but not sure how to implement it in ArcGIS Pro. A GRASS version of the code is located in this thesis, and the source code is included in GRASS as r.param.scale, but I don't really speak GRASS or C.

In ArcGIS, slope is estimated as explained here. Essentially the equation is the sqrt(rise^2 + run^2). This is different from how Evans (1980) estimated slope.

Evans (1980) estimated slope from a fitted surface, with the formula

z = ax^2 + by^2 + cyx + dx + ey + f, with a,b,c,d,e,f as constants.

Slope is sqrt(d^2 + e^2), while curvature is 2a + 2b.

Does anyone have any knowledge as to how to implement this equation in ArcGIS Pro?

I apologize for a lack of a reproducible example, but if someone points me in the in the right direction, I will definitely add a reproducible example.


texts: Perform calculations for each Polygon in a MultiPolygon object (Earth Engine)

texts:
I have layer that contains many polygons. I want to run calcuation for each polygon in the layer so for each polygon I have the mean and the std. deviation and from this I can run calculation for each polygon in the layer individually.

Right now the situation is that when I calculate for 1 polygon only I get the mean and std. deviation for it so my calculation is good because the data is based on the statistics of this one polygon:

But if I have layer with more than one geometry, as you can see, It takes the mean and std deviation from all the polygons and nor for each one so my calculation is not as I want it to be:

As you can see, the right side polygon got diffferent values because the mean and std deviation are based on all the  polygons.

Is there any way to make GEE to take it separately for each polygon and not only for  the whole collection?
I read about reduce region but as you see it seems like it doesn't work...

My end goal is to get statistic for each polygon separately so I can run my calculation and get info about each polygon .

This is the calculation part:

//calculation
var std2 = ee.Number(tableWithStats.get('NDVI_stdDev')).divide(2);
var mean1 = ee.Number(tableWithStats.get("NDVI_mean"));

print('std2',std2,'mean1',mean1);

// the classes borders
var negBorder=mean1.subtract(std2);
var posBorder=mean1.add(std2);

print('negBord:',negBorder,'posBord:',posBorder);

//create the layers
var imageNDVI=MyImage.select('NDVI');
var gtPOS=MyImage.gt(posBorder).selfMask().rename('range');
var ltNEG=MyImage.lt(negBorder).selfMask().rename('range');
var betMEAN=MyImage.gt(negBorder).and(imageNDVI.lt(posBorder)).selfMask().rename('range');

var pos=gtPOS.select('range').multiply(3);
var neg=ltNEG.select('range').multiply(2);

var ndviClassCol = ee.ImageCollection.fromImages([neg,betMEAN, pos]);

// Mosaic the ImageCollection.
var ndviClassImg = ndviClassCol.mosaic();

// Display the classified mosaic to the Map.
Map.centerObject(geometry, 16);

Map.addLayer(pos,{palette:['9370DB']});
Map.addLayer(neg,{palette:['0fff50']});
Map.addLayer(betMEAN,{palette:['fff200']});

Map.addLayer(ndviClassImg, {palette: ['ff0000', '00ff00', '006622'], min: 1, max: 3},'mosaic');

Edit: When I try the code that Daniel suggested that is not client side   on layers with many plots, I get different results.

This is when I check one individual plot out of 2000:

and when I run the code on all the 2000 plots, I wouls excpect to get the sam resulr but I get something different:


texts: QGIS 2.18.4 - Select by location returns no selection

texts:
I have a point layer (green points) and a polygon layer (orange line) which a number of the points fall within.

I have double checked and both layers are in the same CRS.

I have tried multiple functions in the normal 'select by location tool', within, intersect etc but the tool runs and returns no selection.

I have also downloaded the 'Select Within' plugin. This also does not work giving me the error below.

I am at a loss as to why the select by location tool returns no results and also do not understand the error generated by the 'Select Within' plugin.


texts: pgRouting cannot find extension control file

texts:
I am trying to use pgrouting on OS X. I installed the following:

brew install postgres
brew install postgis
brew install pgrouting

When I try to create the extension pgrouting I get the following error:

ERROR:  could not open extension control file "/usr/local/Cellar/postgresql/9.4.0/share/postgresql/extension/pgrouting.control": No such file or directory

Also when I list all installed extensions I only get the following:

\dx
                                     List of installed extensions
  Name   | Version |   Schema   |                             Description
---------+---------+------------+---------------------------------------------------------------------
 plpgsql | 1.0     | pg_catalog | PL/pgSQL procedural language
 postgis | 2.1.5   | public     | PostGIS geometry, geography, and raster spatial types and functions
(2 rows)

How to I have to install pgrouting in order to use it?
I am on OS X 10.10.1


texts: Google Earth Engine: Tasseled Cap Index calculation and extraction of the mean value for each polygon of a shapefile

texts:
I have been trying to use a shapefile consisting of bounding boxes of cities in order to extract the mean value of the tasseled cap vegetation index. I inserted the shapefile as a FeatureCollection and filtered the Image Collection of Landsat 8 based on it.

What I want to achieve is to find a way to calculate the mean value of the index for each polygon and then export them in a dictionary or a table or any other format that I could further analyze it in GIS software. This means that the values should be connected to an ID attribute that each polygon of the shapefile has.

  var boxes = ee.FeatureCollection("users/dimmegkio/CH_GEE_preprocessing");
  var calculateTasseledCap = function (image){
  var b = image.select("B2", "B3", "B4", "B5", "B6", "B7");
  //Coefficients are only for Landsat 8 TOA
    var brightness_coefficents= ee.Image([0.3029, 0.2786, 0.4733, 0.5599, 0.508, 0.1872])
  var greenness_coefficents= ee.Image([-0.2941, -0.243, -0.5424, 0.7276, 0.0713, -0.1608]);
  var wetness_coefficents= ee.Image([0.1511, 0.1973, 0.3283, 0.3407, -0.7117, -0.4559]);
  var fourth_coefficents= ee.Image([-0.8239, 0.0849, 0.4396, -0.058, 0.2013, -0.2773]);
  var fifth_coefficents= ee.Image([-0.3294, 0.0557, 0.1056, 0.1855, -0.4349, 0.8085]);
  var sixth_coefficents= ee.Image([0.1079, -0.9023, 0.4119, 0.0575, -0.0259, 0.0252]);

    var brightness = image.expression(
            '(B * BRIGHTNESS)',
            {
                'B':b,
                'BRIGHTNESS': brightness_coefficents
                }
            );
  var greenness = image.expression(
    '(B * GREENNESS)',
            {
                'B':b,
                'GREENNESS': greenness_coefficents
                }
            );
  var wetness = image.expression(
    '(B * WETNESS)',
            {
                'B':b,
                'WETNESS': wetness_coefficents
                }
            );
  var fourth = image.expression(
      '(B * FOURTH)',
        {
          'B':b,
          'FOURTH': fourth_coefficents
          }
        );
  var fifth = image.expression(
      '(B * FIFTH)',
        {
          'B':b,
          'FIFTH': fifth_coefficents
          }
        );
  var sixth = image.expression(
    '(B * SIXTH)',
    {
      'B':b,
      'SIXTH': sixth_coefficents
      }
    );
  brightness = brightness.reduce(ee.call("Reducer.sum"));
    greenness = greenness.reduce(ee.call("Reducer.sum"));
    wetness = wetness.reduce(ee.call("Reducer.sum"));
    fourth = fourth.reduce(ee.call("Reducer.sum"));
    fifth = fifth.reduce(ee.call("Reducer.sum"));
  sixth = sixth.reduce(ee.call("Reducer.sum"));
  var tasseled_cap = ee.Image(brightness).addBands(greenness).addBands(wetness)
                             .addBands(fourth)
                             .addBands(fifth)
                             .addBands(sixth).rename('brightness','greenness','wetness','fourth','fifth','sixth')
  return tasseled_cap;
};

var start_date = "2018-12-30"
var end_date = "2018-12-31"
var cloud_cover = 10

var select_2018 = ee.Image("LANDSAT/LC08/C01/T1_RT_TOA/LC08_217073_20180613");
var BRD = boxes.geometry();
Map.centerObject(BRD);

var landsat8_collection = ee.ImageCollection('LANDSAT/LC08/C01/T1_TOA')
        .filterDate('2013-09-01', '2013-12-31')
        .filterMetadata('CLOUD_COVER', 'less_than', cloud_cover)
        .filterBounds(BRD)

var landsat8_tasseled_cap = landsat8_collection.map(calculateTasseledCap);
console.log(landsat8_tasseled_cap.getInfo())
Map.addLayer(landsat8_tasseled_cap,{},'Landsat 8 Tasseled Cap');
Map.addLayer(landsat8_tasseled_cap,{min: 0, max:1, bands:['brightness']},'brightness');
Map.addLayer(landsat8_tasseled_cap,{min: 0, max:1, bands:['greenness']},'greenness');
Map.addLayer(landsat8_tasseled_cap,{min: 0, max:1, bands:['wetness']},'wetness');
Map.addLayer(boxes,{},'Chinese cities bounding boxes');


texts: Script to batch export pdf causing incorrect inset map

texts:
I have written a python script to batch export several pdfs from a single mxd, switching layers on/off and naming the output pdf based on values supplied by the user in a csv file. We import the required datasets into a template mxd I've set up containing pre-named groups and the company sidebar. The template has a large data frame in the centre of the page for the main data and a smaller "inset map" data frame in the sidebar below the legend to show the datasets' position relative to known geographic locations. The geographic location of the main data frame is shown in this inset map.

The issue we're having is that the first pdf export behaves as expected, but the inset map in the rest of the exports shows a larger "main data frame" extent than it should. I've attached two images to illustrate what we're seeing.

The python script is quite long (as it's part of a larger toolset script), but if you need to see it I can add it (though all it does is read a csv, switch layers on/off and export as pdf - no scaling or moving dataframes).

Do I need to refresh / regenerate some values before each export?


texts: Building progress dialog with pythonaddins in ArcMap?

texts:
As a component of my add-in, I compute a bunch of summary statistics when the user clicks a button.  This can take a while to process and it can appear as though the program has stopped responding.  I would like to add a dialog box which can indicate to the user that the program is still processing.  I was interested in using wxPython, but I don't necessarily want to rely on users installing dependencies that aren't included with the standard bundle.

Any good ideas?

I was considering something like this (except for my concern above):


texts:
I have to create some viewsheds. I want to use the grass gis r.los.
DEM as elevation, coordinates of a point collected via "coordinate capture" and the area is restricted by a second raster-layer.

My problem is this error:

Does somebody know what to do?


texts:
I have a problem with the infowindow of one of my sublayers in a Carto map. First of all I'll show my code. I create the layer options and then I get my sublayer. I set its options and a first version of the infowindow. After that, I create the event onclick, what changes the infowindow based on a variable I've declared before.

var subLayerOptions2 = {
          sql: "SELECT * FROM zonas WHERE cartodb_id IN (SELECT zonaid FROM devices_zonas WHERE deviceid IN (SELECT deviceid FROM users_devices WHERE userid=" + this.auth.getUserId() + "))",
          cartocss: "#layer {polygon-fill: " + ConstantesColores.colores[3] + "; polygon-opacity: 0.25; line-color: " + ConstantesColores.colores[3] + "; line-width: 1; line-opacity: 1;}",
          interactivity: ['cartodb_id']
}

// Creacion e introduccion del layer en el mapa
var sublayer2 = layer.getSubLayer(2);
sublayer2.set(subLayerOptions2);
sublayer2.infowindow.set('template', infow);
sublayer2.on('featureClick', (e, latlng, pos, data, layer) => {
              sublayer2.infowindow.set('template', infow.replace("#DispAsociados#", arrayZonasConDisp[data.cartodb_id].info));
              alert("Hey! You clicked " + data.cartodb_id);
});
sublayers.push(sublayer2);

The problem is that, the first time I click an element of the map, the infowindow is not updating well, the variable is empty, but when I click a second time it appears well. If I change the element I click, the first time appears the info corresponding to the last click, not this one. It's like a delay. It always show the info of the previous one.


texts: Can't import multiple features or rasters into file geodatabase

texts:
I am trying to use the import multiple features into geodatabase option in ArcCatalog (ArcGIS 10.1). Or, similarly, I try to use the import raster option, which screen layout is similar to the multiple features screen, allowing you to import multiple rasters at once.

The process goes through with a green check mark, but when I expand the geodatabase, nothing is there, not if I refresh, not if I close the program and reopen.

If I use the option to import a single feature class into a geodatabase, it goes fine. That is what I have been doing as a work around. However, there is no such "import single raster" screen, so now I'm stuck as I need to import rasters.

Any ideas of what could be the problem and how to fix it?

I'm teaching a GIS course at a university, so am scratching my head about this problem, as my students are experiencing the same problem on their computers as well.

Here is the dialogue after it is finished--nothing unusual that I can note (used ~~~ to replace personal info here):
Executing: FeatureClassToGeodatabase C:\Users\~~~\Documents\ArcGIS\Data\Canada\ontario_coastline.shp;C:\Users\~~~\Documents\ArcGIS\US_coastline.shp C:\Users\~~~\Documents\ArcGIS\~~~.gdb
Start Time: Tue Oct 13 22:20:03 2015
Running script FeatureClassToGeodatabase...
Completed script FeatureClassToGeodatabase...
Succeeded at Tue Oct 13 22:20:03 2015 (Elapsed Time: 0.00 seconds)

And this is the output from the results menu from another attempt with just one file, but using the import multiple command:


texts: Importing multiple raster images in PostGIS

texts: I have several raster images for one area and I am trying to import these multiple raster images in PostGIS. Is there a efficient way to import these raster images at once?


texts: Referencing Qt UIs external to plugin with python

texts:
I'm developing a plugin that sets up a working environment for a user: loading layers (from db), applying styles and several custom designed (Qt) UI files etc.

Primary user interaction is in the Attribute editor (Attribute Table) where there's a mix of Autogenerated and custom designed UIs.

Some of the custom UIs contain Qt objects/widgets that I need to interact with e.g., a QPushButton and a QTabWidget with QWidgets.

Since the plugin is only setting up the working environment, letting the user work freely afterwards, how can I react to user actions/events calling python functionality?

From the following posts I get the notion that my plugin needs to act as the widgets launch platform to be able to send signals and listen on slots.

How to reference to my in QT designer created buttons in python?
Connect not working for custom QGIS form

So, in short, I'm trying to interact with custom (Qt) forms external to my plugin.

How can I accomplish this?


texts: Select Features by Polygon - Not Displaying the Shape Drawn

texts:
I am a relative newbie to QGIS 1.8, have looked through the previous answers, but have found nothing that resolves my issue.

I have a number of datazones as polygons on a layer.  One of the datazones contains two islands that I need to split.  With the layer selected and visible, editing turned on, I choose Select Features by Polygon.  I click on a point in the sea around the island, move my mouse to another point, click again, a line appears, repat this until I get to the point where I want to put the last dot and right click.  What I understand should happen at this point is that a filled shape should appear, but it doesn't, although the islands become selected.

Is this a bug or am I missing something obvious?  The same is true of all of the vector layers from other shapefiles.

Thanks a lot

Andy


texts: Esri Maps for SharePoint installation hangs?

texts:
We're trying to install Esri Maps for SharePoint 4.0 using the installation guide's PowerShell commands.

The guide has the following commands:

    Add-SPSolution -LiteralPath C:\Software\ESRI.ArcGIS.Mapping.SharePoint.wsp
Install-SPSolution ESRI.ArcGIS.Mapping.SharePoint.wsp -WebApplication http://<web application name> -GACDeployment
Enable-SPFeature -identity "ESRI.ArcGIS.Mapping.SharePoint.wsp13_ArcGIS Map WebPart" -URL http://<web application name>/sites/MySiteCollectionName
WebContentInstaller.ps1 -install
UpdateCustomBCSConnectors.ps1 –ServiceContext http://<web application name>/sites/MySiteCollectionName

Note well - this is a poor installation approach, as it should wait for the install to complete before enabling the feature.

The Install-SPSolution command hangs forever, or at least for an hour.

I'm guessing that the solution has an event receiver within it that tries to contact ESRI (why it would do that during a solution install I don't know), however I can't be sure.  Because of our firewall setup this request hangs forever - essentially freezing the install.

Am I on the right track here, or is there another solution to this?  At this stage it's looking like we have to:

Identify what site it's trying to access
Create a new hosts file entry denying access rather than freezing with the Firewall


texts: MakeXYEvent setting values in CSV with numbers as text in result?

texts:
I thought I had a working ETL script going that was writing a CSV from a non-spatial Oracle DB to ArcGIS gdb when one day it broke. I whittled it down to the MakeXYEventLayer process.  Basically I have a field from the CSV that has a bunch of null values and some numerical values.  Depending on how the tool reads the CSV, if it sees a null value first it sets the field as text.  If it sees a number in that row it'll set it as integer.

-Would a sort on the CSV work to help solve this issue?
-Or is there anyway to force the MakeXYEventLayer results into an already existing shapefile with field set to integer?


texts: QGIS rendering TIFF color with wrong color

texts:
When I try to display a GeoTIFF image with single band pseudocolor, the style seems to be wrong for values out of the range 255.

The legend is shown as

The green pixels have a value larger than 255. I think the problem might be related with integer range but my GeoTIFF has a data type of int32.


texts: I want to change pixel values of elevation data. I want to change values meter to feet. So i need to multiply pixel values with 3.2808. How can i do that in arcmap?


texts: Using Dimension Feature in ArcGIS for Desktop?

texts:
I would like to use the Dimension Feature tool in ArcGIS - this button  as I want to have a label on every side of some polygons (they are all reasonably simple polygons).

When I click on my polygon to add the dimensions (and don't move or click the mouse again), they come up like this which is how I want them:

However as soon as I click the left mouse button to place them, they go like this:

Am I missing an option somewhere? it does not allow me to select 'Dimension Feature' in the Construction tools section of the create features dialog. Its only available as a button (that I found by adding the button manually). The closest construction tool I can select is Dimension Edge, but that means I have to add dimensions for each edge of the polygon. Its a slower process.


texts:
I have a problem concerning the selection by attribute with python on arcmap.

I use this code:

import arcpy
from arcpy import env
arcpy.env.workspace= "E:\cities"
Mi=(r"E:\cities\rivers.shp")
rows = arcpy.SearchCursor(Mi,'"NAME" = \'Mississippi\'')
row = rows.next()
largeur = row.getValue("MILES")
arcpy.selectLayerByAttribute_management("rivers","New_SELECTION",'"MILES" >' + str(largeur))

And i got this error:

  Traceback (most recent call last):   File
  "E:\cities\scripts\Selection.py", line 8, in
      arcpy.selectLayerByAttribute_management("rivers","New_SELECTION",'"MILES"

    ' + str(largeur)) AttributeError: 'module' object has no attribute 'selectLayerByAttribute_management'

  Échec de l’exécution de (SelectionByAttributes)

Can u help me please??


texts: How can I reduce the size of my file geodatadase?

texts: We're using file geodatabase for tax mapping. we save building photos and there are 6000 buildings. The size of this data 2GB. So loading the data takes a lot of time we are not unable to do work easily. If there is any options to solve this problem, please share them with me .


texts: GeoServer PostGIS layer with mixed SRS possible?

texts:
I'd like to setup a layer in GeoServer based on a table in PostGIS DB. The geometry column in the table has no SRID set (Find_SRID for this table returns 0).
Geometries inserted into the table don't have the same SRID. Is it still possible to create a GeoServer layer from this table?

When I try to setup the layer, GeoServer sets "Native SRS" to one of the SRS present in the table data and I cannot remove it.

If its not possible, is the some possibility in GeoServer to serve global GIS data with mixed SRS?


texts: QGIS: Distance units in CRS layer: WGS 84

texts:
First post so please bear with me.

When I run a vector analysis on two shapefile layers with multiple points, I get a distance as my answer.  What are the units on this distance?


texts:
I have a dataframe with around 40 points. They are grouped by an ID variable (their name), which four points in each group. I'm trying to turn these points into polygons, but I've had a lot of trouble with different methods in both sf and sp — currently I have them as a large sf dataframe, but I can't seem to find any function that would use a grouping variable to create polygons.

I've tried using the st_cast method, the general idea being

polypoints <- polypoints %>% group_by("ID") %>% st_cast("MULTIPOINT") %>% st_cast("MULTILINESTRING") %>% st_cast("MULTIPOLYGON")

but the step from multipoint to multiline string leaves the geography empty — what are the steps for a process like this?


texts: Most efficient way to batch clip in FME?

texts:
I have one large shapefile and 100 smaller shapefiles that are smaller areas within the large shapefile.

Using the Clipper transformer, I want to clip the large shapefile using each of the smaller shapefiles (I want to do the clip 100 times), and then write out 100 resulting shapefiles.

What is the most efficient way to do this in FME?


texts: Compute for shortest/perpendicular distance between track point and centerline

texts:
I am working on a python school project that uses gpx tracks from different devices. Latitude, Longitude, Elevation and Timestamp are the only information that are constantly available. HDOP, Speed, etc. information may or may not exist.
Projection used by the devices: WGS84. The ultimate purpose of the program is the computation of travel time given a specific route.

CURRENT NEED: Snap gps track points to the nearest centerline.

CURRENT TO DO: Compute for shortest/perpendicular distance, value of which will be used in deciding which line the point snaps to.

QUESTION 1: What would be more accurate to do/use?

QUESTION 2: What would be faster computationally?

OPTION 1: Convert to UTM first.

OPTION 2: Use Haversine's Formula.

OPTION 3: Use Vincenty's Formula.

CURRENT READINGS:

Calculate distance, bearing and more between Latitude/Longitude points: http://www.movable-type.co.uk/scripts/latlong.html
Distance Calculation Algorithms:
http://www.ga.gov.au/scientific-topics/positioning-navigation/geodesy/geodetic-techniques/distance-calculation-algorithms
Vincenty vs Haversine Distance Calculations?
http://www.researchgate.net/post/Vincenty_vs_Haversine_Distance_Calculations
Geographic Information Technologies for analysing the digital footprint
of tourists: [http://www.agile-online.org/Conference_Paper/cds/agile_2014/agile2014_91.pdf4
Adventures in GIS Programming and Development: a Geospatial Tutorial: http://www.toptal.com/gis/adventures-in-gps-track-analytics-a-geospatial-primer
Is there a good GPS track analysis library? Is there a good GPS track analysis library?

WHAT I HOPE TO SOLICIT BY POSTING THIS QUESTION: Experience-based answers.


texts: What Windows OS does ArcGIS Desktop perform best on?

texts:
I'm planning to install ArcGIS Desktop 10.6 on a completely clean system that will only be used for GIS. The system meets the requirements for 10.6, but as stated by ESRI, 10.6 can be installed on Windows 7, 8, 10 or Windows Server 2008, 2012 or 2016. Since I'm almost totally flexible, I'm interested what OS Desktop performs best on.

Is there any evidence to support one version of Windows over another when choosing to install ArcGIS for Desktop on a Windows machine?


texts: Self hosted, highly accurate, address geocoding?

texts:
I'm working on an app that requires Rooftop level accuracy coupled with very high volumes. The client would prefer to self-host the solution which might not be a bad idea considering our volume and response time requirements.
I'm not seeing how the usual suspects (Google, MapQuest, ...) might help and the geocodes we get from address correction software are inconsistent quality-wise.
Any pointers or suggestions to vendors and/or solutions are appreciated.
FWIW we're currently using PostGIS.


texts:
I'm trying to execute the raster calculator function in VB.NET but I'm having trouble with the gp.execute line. Sorry if the code is confusing, I'm fairly new to VB.net. I keep getting this error:
An unhandled exception of type 'System.Runtime.InteropServices.COMException' occurred in ESRI.ArcGIS.Geoprocessor.dll
Additional information: Error HRESULT E_FAIL has been returned from a call to a COM component.

Here's my code. I thought it might have been the expression, but the error occurs at the execute command.
Private Sub ExecuteRasterCalc()
    Dim raster1 As String = "T:\TeamProject\Group2\Ca Slope and Aspect\CA_ASP_RISK1.tif"
    Dim raster2 As String = "T:\TeamProject\Group2\Ca Slope and Aspect\CA_SLP_RISK1.tif"
    Dim raster3 As String = "T:\TeamProject\Group2\Ca Landcover Fire Intensity\CaFIRE_INTEN11.tif"
    Dim raster4 As String = "T:\TeamProject\Group2\Ca Biomass Load\FUEL_LOAD1.tif"
    Dim expression As String = "'" & raster1 & "'+'" & raster2 & "'+'" & raster3 & "'+'" & raster4 & "'"
    Dim outputpath As String = "T:\TeamProject\Group2\Test.gdb\rastercalc"
    Dim gp = New ESRI.ArcGIS.Geoprocessor.Geoprocessor()

    gp.OverwriteOutput = True
    gp.AddOutputsToMap = True

    Dim parameters As IVariantArray = New VarArray
    parameters.Add(expression)
    parameters.Add(outputpath)

    Dim result As IGeoProcessorResult = gp.Execute("RasterCalculator", parameters, Nothing)
End Sub


texts: How are satellite images mapped to terrain coordinates so precisely?

texts:
Many geoinformation systems have both satellite images and terrain features (houses, roads, etc) with coordinates. So they can show terrain features drawn over the satellite image or they can display an "it's here" mark on a satellite image given certain coordinates.

How does mapping happen (in layman terms)? A satellite flies around Earth and makes a gazillion satellite images, these shots are then sent to map producer. How does the map producer know which point on a satellite image correspond to which coordinate?


texts:
I am trying to load a series of csv files as vectors in PyQGIS. I run this code to try to get the first one:

uri='file://D:/Saved Stuff/gulags/1.40.csv?delimiter=;&yField=y&xField=x'
layer = QgsVectorLayer(uri, '1.40', 'delimitedtext')
layer.isValid()
False

Obviously (?) it won't load when I run QgsMapLayerRegistry.instance().addMapLayer(layer). When I add the csv manually to a map it loads fine. Here is what that looks like

  name;y;x;total

  AMURSKII ZhELEZNODOROZhNYI ITL;51.375889;128.134147;84053

  BELOMORO-BALTIISKII ITL;62.912774;34.456850;78918

  BIRSKII ITL;48.999752;132.461997;11634

Can anyone help me with this?


texts: Find NoData values inside a polygon area from a raster

texts:
I have a raster (representing a polygon in black) which contains NA values (represented by white pixels) that are located inside the polygon and outside the polygon.

How can I get the index value for the NA values inside the polygon (delimited by the red contour) in the raster ? Next, I would like to use this index to replace the NA values as in the question Change the values of NA cells in a raster by using a geographic subset of the raster as a condition. In this question, the function "cellsFromExtent" is used to get the index value for the NA values in the extent. In my case, I would like to extract NA values within a polygon.


texts: Selecting features that do NOT intersect in ArcGIS Desktop?

texts:
I used the "create spatially balanced points" in Geostatistical Analyst in ArcGIS 10 to generate a large set of points with equal inclusion probability but a "spatially balanced" dispersal. I then buffered these using different buffer settings, with the idea of explicitly testing the influence of sampling unit size on population estimates (in my case, tree cover) using a fixed area sampled. This latter detail is key - I sum of the area for the sum of the selection of buffered points to be the same between comparisons. A toy example: what's the difference sampling 10 polygons 1 ha each in area (10 ha total) versus 100 polygons, each measuring 0.1 ha (again, 10 ha total). The point layer created by the tool qualitatively has the desired dispersion, but when you start buffering points, some start to overlap, reducing the total area (plus violating any possibility of independence between samples).

What I'm looking to for is an efficient way to find features that DO NOT intersect. I saw a post asking a similar question for PostGIS, but was looking for a way to approach the question using ArcGIS.


texts: Finding nearest points to an interpolated point along a Shapely linestring


texts: Getting layer style "name" L.control.layers on map click?

texts:
I am adding map styles in L.control.layers like this:

L.control.layers({
        'Mapbox Light': L.mapbox.tileLayer('mapbox.light',{name:"base"}),
        'S_S': L.mapbox.styleLayer('mapbox://styles/geotrek/s***',{name:"F1"}).addTo(map),
        'S_F': L.mapbox.styleLayer('mapbox://styles/f****',{name:"F2"})
    }, null).addTo(map);

I can get the layer name on change

 map.on('baselayerchange', function(l){
        console.log("baselayerchange",l);
        console.log("baselayerchange",l.layer.options);
        console.log("CONTROLS:", controls)
    })

How can I get the current layer style name like S_F or S_S or Mapbox Light or options name on map CLICK event?


texts: QField showing position is outside raster image?

texts:
I am a QGIS beginner and I have been asked to map an area that is being reforested. The project itself is quite easy, I have to generate a shapefile containing the polygons that delimit areas with different tree species, using aerial images as backgound. The map should be displayed in QField, to facilitate the operation in the field. They did not give me any georeferenced aerial images, so I made a collage with screenshots from satellite immages using Photoshop, and I georeferenced the raster (WGS 84) with the georeferencer plugin of QGis, using coordinates that I got from Google Earth. Before making the shapefile I wanted to try the project in QField to make sure that the localization was working well but my position is outside the raster image, where it should be.

What could be the problem?


texts:
I want to convert these points to lines, however they are not sorted sequentially along a line (by ID or FID), i.e. the converted lines zig-zag across the obvious trail of points.

How can I join these point trails based on proximity to points?

I don't have access to the Near (Analysis) tool, which I thought might help me move froward.


texts: What specifically is an 'attribute' and how does it work with the related chart?

texts: I'm a new QGIS user (and GIS technology in general) and I'm trying to learn its fundamental features. I have a big dilemma about vector shapefiles, and points in particular: basically, I don't understand what specifically is an 'attribute', and how it works with the related chart (I guess that the 'ID' is required for a search within that chart) Please help me! Sorry for my bad english.


texts: Viewing in-line variable value in model parameter dialog box of ModelBuilder?

texts:
I've developed a model that summarizes road attributes (e.g., volumes) across a digitized polyline. The model has input/parameter dialog boxes for a road feature class, workspace, and attribute summary table.

These parameters are entered every time the model runs. I've provided default values but users want to save their latest parameters. I'm using ArcGIS 10.3 with Windows 7.

Instead of providing the entire model, I've provided a part of it below. The image shows my dilemma. I can store and retrieve model parameters in a table called C_Parameters.

Though I can provide the NAME of a parameter variable (e.g., %RoadFeatureClass%)  in the input/parameter dialog box in a ArcToolbox tool, I can't display the parameter variable VALUE in the dialog box.

Basically, I wish the Input Features box of the Copy Features tool would be C:\SEWRPC\CORDEN.gdb\Base_out_Link3 instead of RoadFeatureClass.

I've tried using Python's object methods (e.g., %RoadFeatureClass%.title() and print(%RoadFeatureClass%)),which didn’t help.

Do you have any suggestions how to display a variable value in a ModelBuilder input dialog box?


texts: QGIS Composer causing QGIS to crash when opening a certain template

texts:
Every time I try to open a certain template in the Print Composer, I get a blank window and no response from the program.
This then causes the whole program to crash.

This particular template is rather important and I need to access it, one of the possible causes of the crash could be the big amount of maps in the template and the rendering required. However, this issue has only just arisen as last I tried to open the template it was working fine.
Any advice or potential fixes?
Thanks


texts: Missing classes from qgis.core module only in PyCharm's editor

texts:
I would like to use PyCharm Community Edition 2020.1.2 x64 on Windows 10 as IDE for writing QGIS plugins. I need to base on old version - QGIS 2.18, so I need python 2.7, qt4, and also core modules from QGIS 2.18. I downloaded QGIS 2.18 from official site. It provides a ready to use python environment as .bat script, that sets all needed python paths. So I used path D:\Program Files\QGIS 2.18\bin\python-qgis-ltr.bat to set main python interpreter for my project.
After that I was able to run python console in PyCharm, import all QGIS modules, and any classes from QGIS standard collection like QgsFeature, QgsVectorLayer etc. Everything was fine untill I tried to import this classes from PyCharm editor in .py file as from qgis.core import QgsFeature. In editor PyCharm sees the package qgis.core and I'm able to import it, but PyCharm's code completion inside the module sees only classes that are saved in D:\Program Files\QGIS 2.18\apps\qgis-ltr\python\qgis\core\__init__.py e.g. QgsEditError class, but all standard QGIS class underlines on red.
Python console ran on exactly the same interpreter see the same sources as are visible in "External Libraries" tab in PyCharm, so command:

qgis.core.__path__ returns ['D:\\Program Files\\QGIS 2.18\\apps\\qgis-ltr\\python\\qgis\\core'],
qgis.core.__file__ returns 'D:\\Program Files\\QGIS 2.18\\apps\\qgis-ltr\\python\\qgis\\core\\__init__.pyc,
dir(qgis.core) returns full list of QGIS standard core classes.

If I run the .py script, in which PyCharm's code validator is showing incorrect import object by red underline, the script is executed without any error and I can easly print the imported QgsFeature object. Why my PyCharm's editor see different content of qgis.core module then console ran also from IDE on exactly the same interpreter?


texts: How to add the geocode filtering to gxp.plugins.GoogleGeocoder

texts:
The Google geocode provides the capability of applying filters/restraints to the search, e.g. restricting the search to an area with a specific postal code, etc... I want to accomplish the same with the gxp_googlegeocoder plugin. Looking through the API I cannot seem to find how this can be implemented.

Has anyone else been able to accomplish filtering of sort on this plugin, if so what needs to be added/done?

Here is the link to the Google API.
Here is the link to the GXP API.


texts: Using QGIS' "network protocol" to import data from OPeNDAP?

texts:
After looking at the OPeNDAP documentation
I assumed I could import data into QGIS by downloading the .ascii file OPeNDAP generates. However the downloaded OPeNDAP .ascii doesn't seem to be compatible with the ERSI .asc import methods on QGIS.  This website says that you can import OPeNDAP and GeoJSON vector data via QGIS "network protocol" but future searches on google and the QGIS documentation (module 11 online resources) hasn't helped.

Does accessing data via network protocol mean using the Add WFS function of QGIS?

Can you only open OPeNDAP vector data and not raster data on QGIS?

This is the ascii file I need from OPeNDAP. It should be a "Basin Characterization Model" that is a raster with multiple layers.

Does anyone have experience getting data from OPeNDAP into QGIS?


texts: How to create an editable personal copy of a read-only PostGIS layer in QGIS 3.4

texts:
I have a few vectory layers in a PostGIS database. I want to give the other users of the database read access to these layers. I don't want them to be able to alter them, though. Instead, if anybody wants to make their personal adjustments, they should create their personal copy, of which they are the owner. In PostgreSQL,

CREATE TABLE personal_table as public_table;

does what I want. Of course, this can also be done from the DBManager.

Is there any other way to accomplish that within the GUI of QGIS 3.4 for users, who are not comfortable with learning code commands?

Usually, the users drag and drop the layers they need from the PostGIS menu in the browser. Rightclick==>duplicate only creates a copy of the layer in the project but doesn't create a new table with ownership, so the user still cannot edit.


texts: Python Question - how do I extract a part of a string

texts:
I'm building a tool in python and for that I have a question:

I have a string such as "Denver.dwg Group Layer\Denver.dwg Annotation".

How do I assign to a variable (for later print) only the first part of the string, "Denver.dwg" ?? I need to provide a general solution, as the city name differs from case to case (it might be shorter or longer) and so counting the places until the g of dwg (10 places) won't help.

(I'm writing the python tool for arcgis 9.3.1)


texts: How to update a text column of selected features in QGIS?

texts:
Let's say I have a bunch of postcodes. After making a selection and having the records highlighted in the attribute table, how to update all of them and give them the same category?

In ArcGIS, after having the features selected, you simply had to use field editor and say for example: name = asd.

In QGIS that doesn't work. How to do it in QGIS?


texts:
I have

linestring that is calculated on a spheroid.
Shpere distance of where the point should be on the linestring
Distance spheroid / sphare.

On one lib, its all calculated in spheroid (Normal map).
Other lib uses sphare.
Lets say the total distances are:
Spheroid: 76363.06884857068
Sphere: 82238.6290291082
The sphareoidic distance from the start would be 409, I need to get distance from spheres perspective.
Any tips where should I look?


texts:
I created form-value related dropdowns in a QGIS project that work perfectly. The Idea is to collect data of points in the field. All types ("Art") have subtypes ("Subart") that are dependant on their type.

Here some examples:

As I opened the project in QField, the dropdowns were blocked. Is QField not capable to handle related dropdowns?


texts: Is it possible to use spatial function to retrieve non-spatial data?

texts:
I'm using PostGIS  in postgresql.I wish to calculate points within the radius from particular geocode.For that i have created a table with spatial column geocode having datatype as geometry and inserted some values using following query.

INSERT INTO aspatial_address( addressline,country,geocode)
  VALUES ('','', ST_GeomFromText('POINT(14.083333 70.583333)', 4326));

As geocode having geometry datatype inserted values of geocode are stored in hex format.Ex. for "POINT(10,19.1)" it stores '010100000000000000000024409A99999999193340' as geocode.

I can retrieve these value in normal geocode format using ST_AsText(geocode)geocode function using :

select addressid,ST_AsText(geocode)geocode from aspatial_address as A where ST_Point_Inside_Circle(a.geocode,14.083333,19.583333,7)

Till everything is fine.

Now i wish to select geocode values from non-spatial postgres tables.
for that i have tried same query but their in incompatibility between datatype of geocode.As spatial have geometry datatype and non-spatial have character varying datatype it gives error.Can i convert datatype through query using any function or i have to change query?


texts: How to tile an Image created in Adobe Illustrator

texts:
I'm a beginner in this and don't understand something.

I have a simple map in Adobe Illustrator, and I want to tiling it to using in webglearth (https://www.webglearth.com/). How I can do it?


texts: How to keep only the coastlines on TIFF file (No ELEV)

texts:
I have a PNG library (524 288 Tiles of 256 * 256px) which I georeferenced via gdal (gdal_translate & gdalwarp => FILE.tiff)
The sea is black (RGB 0, 0, 0) and the land is white (RGB 255, 255, 255).
I would like to keep only the "limit" between land and sea (= coastlines).
I saw that there is gdal_contour but it is always a question of an elevation (ELEV) that I do not have and therefore I do not know how to get these coastlines ...
gdalinfo returns me this structure for each TIFF
Driver: GTiff/GeoTIFF
Files: map_11_480_100.tiff
Size is 181, 313
Coordinate System is:
PROJCRS["WGS 84 / Pseudo-Mercator",
    BASEGEOGCRS["WGS 84",
        DATUM["World Geodetic System 1984",
            ELLIPSOID["WGS 84",6378137,298.257223563,
                LENGTHUNIT["metre",1]]],
        PRIMEM["Greenwich",0,
            ANGLEUNIT["degree",0.0174532925199433]],
        ID["EPSG",4326]],
    CONVERSION["Popular Visualisation Pseudo-Mercator",
        METHOD["Popular Visualisation Pseudo Mercator",
            ID["EPSG",1024]],
        PARAMETER["Latitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8801]],
        PARAMETER["Longitude of natural origin",0,
            ANGLEUNIT["degree",0.0174532925199433],
            ID["EPSG",8802]],
        PARAMETER["False easting",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8806]],
        PARAMETER["False northing",0,
            LENGTHUNIT["metre",1],
            ID["EPSG",8807]]],
    CS[Cartesian,2],
        AXIS["easting (X)",east,
            ORDER[1],
            LENGTHUNIT["metre",1]],
        AXIS["northing (Y)",north,
            ORDER[2],
            LENGTHUNIT["metre",1]],
    USAGE[
        SCOPE["Web mapping and visualisation."],
        AREA["World between 85.06°S and 85.06°N."],
        BBOX[-85.06,-180,85.06,180]],
    ID["EPSG",3857]]
Data axis to CRS axis mapping: 1,2
Origin = (-1252344.271424327744171,7331600.036586093716323)
Pixel Size = (215.927731168155702,-215.927731168155702)
Metadata:
  AREA_OR_POINT=Area
Image Structure Metadata:
  COMPRESSION=LZW
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (-1252344.271, 7331600.037) ( 11d15' 0.00"W, 54d50'37.50"N)
Lower Left  (-1252344.271, 7264014.657) ( 11d15' 0.00"W, 54d29'33.52"N)
Upper Right (-1213261.352, 7331600.037) ( 10d53'56.08"W, 54d50'37.50"N)
Lower Right (-1213261.352, 7264014.657) ( 10d53'56.08"W, 54d29'33.52"N)
Center      (-1232802.812, 7297807.347) ( 11d 4'28.04"W, 54d40' 6.87"N)
Band 1 Block=181x45 Type=Byte, ColorInterp=Palette
  NoData Value=255
  Color Table (RGB with 256 entries)
    0: 0,0,0,255
    1: 255,255,255,255
    2: 0,0,0,255
       [...]
  255: 0,0,0,255

I have this :

and I would like to get this: (by photoshop)

Can you tell me the method?


texts: QGIS: how to label a continent map with distinct values from shapefile?

texts:
I have a continent shapefile containing the names of its different countries. Unfortunately, many of these countries have more than one entry (i.e. more than one row in the attribute table).

Therefore, when I am trying to add labels of the country names to the map, some of country names are printed out more than one time.

I read Getting list of distinct values from shapefile field using QGIS?, but was not able to apply it to my specific problem.

What can I do about this?


texts: Error in creating output ... table already registered?

texts:
I am afraid I am a complete newbie to GIS and ArcGIS and it has fallen on my shoulders to try and resolve a problem with our ArcGIS installation.  Assume I know nothing about GIS or the software please...until 2 weeks ago I hadn't even touched a GIS system.  Anyway we have ArcCatalog 9.1 and ArcMap 9.1 installed on a very old server.

I can see the tables/views created in the database but the user's layers are not being populated (red exclamation marks in ArcMap).

I have opened up ArcCatalog and I can see the database connections on the server under Admin User's app data folders.  If I open up these (there are a number) some or all have at least one of the tables missing or a necessary view missing.

I think the error is related to registering the tables with the sde.

Sorry for any mistakes I have made, I have spent a bit of time on this and will edit/correct the above as I receive further information.

EDIT (content removed above as I try and focus in on the true problem)

I have recreated the table the scripts are failing to create from a database backup.  I can then see this in one of the database connections on ArcCatalog.  However when I try and import the feature class into the other connections it comes up with an error stating:

Table Already registered.....violation of UNIQUE KEY constraint....

Obviously the table exists now but I cannot see anything within the DB Connections relating to the same name.  I believe I need to unregister or re-register the table with the sde...or something?!

Edit - Full error when trying to "copy features":

Error in creating output c:\....... table already registered
[Microsoft OLE DB Provider for SQL Server: Violatio of UNIQUE KEY constraint
'registry_uk2'.  Cannot insert duplicate key inobject 'sde_table_registry'.]
Failed to execute (Copy Features)."


texts: Reverse Geocoding: Add municipality to resultant shapefile?

texts:
I have a list of points with Lat & Long coordinates that I've plotted on the map. I used reverse geocoding to attach the corresponding address. The resultant shapefile returns the lat & long fields as well as street address. However, I will need to add the Municipality to make the data more usable. For example, I need to differentiate the same addresses in different cities (i.e. 123 Test St exists in both City A and City B). Basically, I need my output to contain four things: lat, long, address, municipality.

I've looked into the reverse geocde tool, but I can't seem to select the fields I want in the output. Is there a way to add this to my output or do I need to use a different tool to achieve this?

My address points that I'm using for my address locator has the municipality information that I'm looking for.

I'm running ArcGIS 10.1 (I have the basic license).


texts: I have GeoNetwork (v 3.0.4) and I want to add the Supplemental Information attributes to the default metadata view. How can it be done?


texts:
I'm currently trying to import a 3D contour file into a QGIS project.

How can I export the DXF file from Microstation with height information?

I already have imported it to the QGIS project but can not create a TIN because the height info is missing.

I'm using QGIS Lisboa 1.8.0


texts: Converting a shapefile from decimal degrees to UTM (km) in R

texts:
The code that I have below seems like it works in converting my vector roads which is in decimal degrees into a UTM CRS so that I can more easily work on the data in kilometers rather than decimal degrees.
However, upon inspecting roads.UTM I can see that the extent of the SpatialLinesDataFrame is still in decimal degrees, extent: 24.41195, 32.79441, -17.92888, -9.358345  (xmin, xmax, ymin, ymax), and yet the CRS says it's in +proj=UTM.
How can I convert roads into UTM so that the extent is in UTM and not decimal degrees?
If you want to download the roads shape file you may do so here: It is a shape file from the World Bank, which you can download by click on the link: "Download Zipped Shapefile Here: Zambia Roads".
Here is the code I have:
rm(list=ls())
library(sf)
roads<-readShapeSpatial("Zambia_Roads.shp")
crs(roads)
proj4string(roads) <- "+init=epsg:32735 +proj=utm +zone=35 +south +datum=WGS84 +units=m +no_defs
+ellps=WGS84 +towgs84=0,0,0"
roads.UTM<-spTransform(roads, CRS("+init=epsg:32735 +proj=utm +zone=35 +south +datum=WGS84 +units=m +no_defs
+ellps=WGS84 +towgs84=0,0,0"))
roads.UTM


texts:
Is there a way to find the centroids of highest density areas in the density raster (red areas in the figure) created with the Optimized Hot Spot Analysis in ArcMap10.3.
Theoretically I think I should do something like:

Create poligons of different areas (green, yellow, orange, red  in
the figure)
Select red areas
X and Y of the red areas

I literally have no idea how to do it  the stage 1 in ArcMap while I I know how to work out part 2 and 3.
Or is there an easiest way to find the density of population (for example) and find out the coordinate of the highest density centers without using the Advanced Hot Spots tool?


texts: How to prevent not visible objects (scale rules) from being identified?

texts: In Symbology tab I am using rule-based symbols and I have defined visibility rules dependent on the minimum scale. Above some scale threshold some features are not visible as I expected. But when I am clicking on not visible feature with identify tool it is still listed in the identify results. How to prevent it?


texts:
I'm working on Qgis 3.8 , I have a mold filled with squares, I would like to burn the average value within each square of other vector or raster file.

This is the vectoriced file with all the data that I want to burn in those squares

This is the original raster file

In the past I have done it (this isn't the final file, it's similar to how it should look), but now I can't remember how, neihter find the post that helped me with that.

Lowering the resolution of the raster file won't be usefull, because I am comparing different rasterfiles that have to have the same resolution and squares with data and without data in order to be compared in R Studio.


texts: How to use select by attributes and edit selected features in a Feature-Class in ArcMap (10.4) ModelBuilder

texts:
I have a Feature-Class with alot of fields (it is the product of a union of 8 different layers) in which I would like to select by attributes and then edit a field (actually at least two fields in the later process).

At the moment I am stuck because "select layer by attribute" does not work with a Feature-Class nor with a Shapefile and "extracting data" does not work because I need the Feature-Class containing all features for further selection and editing of fields.

Is there a tool that works with ModelBuilder I overlooked?

I am specifically not looking for a solution involving VBA/Python.


texts: Removing white background from GeoTIFF in Geoserver

texts:
I'm having trouble with removing white background/bounding box from raster in Geoserver. When I add WMS service in ArcMap or Autocad I'm seeing this white bounding box and to be fair I don't know how to remove it.

Here is example:

I even tried changing SLD styles in Geoserver but without any luck. Here is my .xml style:

<?xml version="1.0" ?>
<sld:StyledLayerDescriptor version="1.0.0" xmlns="http://www.opengis.net/sld" xmlns:gml="http://www.opengis.net/gml" xmlns:ogc="http://www.opengis.net/ogc" xmlns:sld="http://www.opengis.net/sld">
    <sld:UserLayer>
        <sld:LayerFeatureConstraints>
            <sld:FeatureTypeConstraint/>
        </sld:LayerFeatureConstraints>
        <sld:UserStyle>
            <sld:Name>te</sld:Name>
            <sld:Title/>
            <sld:FeatureTypeStyle>
                <sld:Name/>
                <sld:Rule>
                    <sld:RasterSymbolizer>
                        <sld:Geometry>
                            <ogc:PropertyName>grid</ogc:PropertyName>
                        </sld:Geometry>
                        <sld:Opacity>1</sld:Opacity>
                        <sld:ChannelSelection>
                            <sld:RedChannel>
                                <sld:SourceChannelName>1</sld:SourceChannelName>
                            </sld:RedChannel>
                            <sld:GreenChannel>
                                <sld:SourceChannelName>2</sld:SourceChannelName>
                            </sld:GreenChannel>
                            <sld:BlueChannel>
                                <sld:SourceChannelName>3</sld:SourceChannelName>
                            </sld:BlueChannel>
                        </sld:ChannelSelection>
                    </sld:RasterSymbolizer>
                </sld:Rule>
            </sld:FeatureTypeStyle>
        </sld:UserStyle>
    </sld:UserLayer>
</sld:StyledLayerDescriptor>

Image is in TIFF format with 3 bands. Is it problem in style format, image itself or Geoserver?
Or is it even possible to host raster images without white bounding box as seen in picture above?
Any advice? :)


texts: Needing to find parent zip code of a bunch of missing zip codes in the US

texts:
We're using MapPoint to export zip code / territory combinations. Unfortunately, the export doesn't include PO Box zip codes which mean that we're missing a lot of zip codes

The PO Box zip codes are contained within a parent zip code and so I just need to work out what the parent zip code is for each of the missing zip codes.

I have 2 SQL tables:

ALL zip codes with Lat/Lng and Geography column. Also has a GeometryShape column
standard zip codes with Lat/Lng and Geography column.

Couple of ways I was looking to do it...

either find the closest Geography point to each missing zip code and make the closest one the parent zip (using STDistance)
OR
Convert the geometry shape to a geography shape of the standard zips, then use STIntersects to find which shape the missing zip lat/lng belongs to.

Having a bit of bother with this, first time I've worked with geospatial querying in SQL


texts: Using ArcGIS locators in R to bulk geolocate dataset?

texts: I have ArcGIS on my computer, but my laptop is entirely too slow to do even basic ArcGIS functionality.  I need to mass geolocate a dataset, and have access to the .loc locators provided by ArcGIS.  Is there a library in R, or a programmatic method to use the locators for geoprocessing?


texts: Why osmId is not the same across the same road?

texts:
I download .shp and .dbf files for roads, extract it and save to database. I display points with the same name of street on my own Google Map. This is how it looks like:

I was interested in how osmID works, so on each marker I marked the osmID along with its coordinates. And I found out that there are several different osmIDs along the same street with the same name. Here you can see it (the first line is osmID):

But I would like to get all points across the same road, not partial content. Is there any way to get all points across the same road ? There are several roads with the same name across my country. How osmID really works? Distance between this two points show on a screen is 180 meters, in my project it is too great a distance to be able to implement the idea.


texts:
I have a network data set for footpaths which I have generated through the collation of path networks from a variety of sources.

I want to compare this with another dataset which also contains a path network, to see which is more complete / where there might be features in one dataset and not the other.

It is likely that the features representing the same path will be non-identical between the datasets, but they will be within a couple of metres of each other.

How could I go about this?


texts: How to create network by GTFS data in ArcGIS Network Analyst?

texts:
I want to create my network by using a GTFS file. I already opened a network dataset but when I load the GFTS to network dataset tool it gives me some errors.

Can someone please help? Below is a screenshot of the error:


texts: Setting Layer rendering Blending Mode of Feature using Python

texts:
I'm trying to set the blending mode of a layer like on the picture, Darken on Feature, not on layer, with Python.

Here is the code that i used to modify the "Layer" Blending mode.

layer=iface.activeLayer()
target_field = 'delta_t'
myRangeList = []
symbol = QgsSymbol.defaultSymbol(layer.geometryType())
myRange = QgsRendererRange(0, 10.0, symbol, '0 a 10')
myRangeList.append(myRange)
myRenderer = QgsGraduatedSymbolRenderer(target_field, myRangeList)
myRenderer.setMode(QgsGraduatedSymbolRenderer.Custom)
layer.setRenderer(myRenderer)
layer.setBlendMode(QPainter.CompositionMode_Darken)

How can i modify the "Feature" Blending Mode on "Darken"?


texts: Dynamic segmentation with respect to a control factor

texts:
Is there a way either with python scripts or other solution to dynamically segment a table with linear measures (from-to fields)
Each segmented section (from-to) should have all the item numbers written to a field that is within a certain range.
e.g:

Input table

row1.  A-----c ===> itemID =1400
row2.  B-----F ===> itemID =1401
row3.  D-----F ===> itemID =1402

Output Table

row1.  A-----B ===> itemID =1400
row2.  B-----C ===> itemID =1400; 1401
row3.  C-----D ===> itemID =1401
row4.  D-----F ===> itemID =1401; 1402

Note: the input table is just a dbf file...


texts: How to find nearest neighbor in Oracle's ST_GEOMETRY?

texts:
I am trying to find nearest road segment to a point in Oracle. I can do this with SDO_GEOMETRY type i.e

select
/*+ LEADING(pois) INDEX(route ALL_ROUTES_SPTAIL_INDEX)  */
route.*
from DEVICE_STATUS pois,ALL_ROUTES route where
MDSYS.sdo_nn(route.shape, pois.shape,'sdo_num_res=1') = 'TRUE'
and pois.objectid=729415

How can I find nearest neighbor when both tables have ST_GEOMETRY type? or is there any way to convert
ST_GEOMETRY to SDO_GEOMETRY


texts:
I have a grid polygon layer and a line layer.

Once imported into PostGIS (let's call the tables poly and line), how can one build a query (with ST_Intersect for example) that counts the number of line intersections with polygons (the same as the QGIS tool:  MMQGIS > Spatial Join > Intersects > First)

and returns the poly table with an extra column containing the sum of line intersections per polygon?

In some cases, the same line will intersect the same polygon more than once, and in that case, it should only be counted once.

So far, I have this:

 CREATE TABLE sp_join AS
 SELECT polygons.* COUNT (l.*)
 FROM  polygons as p
 JOIN lines l ON  ST_Intersects(p.wkb_geometry, lines.wkb_geometry)


texts:
This seems basic but I am new to QGIS and I can't find how to do it. I have three shape files each with the same attribute. Each shape file is imported into its own layer. Each shape file has attribute values in the range 1-100, but different sub-sets of those values are present in each file.

I want a color ramp where attribute value 62 (for example) has the same color on all three layers, but the colour ramp I am using assigns a different color to value 62 on each layer. I think it is because the range of values in each layer is different. I think my color ramp is a continuous range and not discrete colours tied to a discrete value.

How do I build a ramp where the full range of values is covered and a given value always has the same color on different layers?


texts: ArcPy Error : Input raster: is not the type of Composite Geodataset, or does not exist

texts:
Below is a script I am working on to extract time series NetCDF to raster. The script is working fine without mask layer. When I tried to assign a mask file to extract ROI before I save it to tiff, it fails with a 000860 error.

for dimension in nc_Dim:
        if dimension == "time":
            top = nc_FP.getDimensionSize(dimension)
            for i in range(0, top):

                dimension_values = nc_FP.getDimensionValue(dimension, i)
                nowFile = str(dimension_values)
                nowFile = nowFile.translate(None, '/')

                if int(nowFile[-4:]) == 2005:

                    dv1 = ["time", dimension_values]
                    dimension_values = [dv1]

                    NetCDF = arcpy.MakeNetCDFRasterLayer_md(inNetCDF, variable, x_dimension, y_dimension, nowFile, band_dimension, dimension_values, valueSelectionMethod)
                    print "success"
                    outname = outLoc + nowFile
                    mask = "D:\data\NEW.shp" #(1)
                    maskdata = arcpy.sa.ExtractByMask (mask, nowFile) #(2)
                    ##after executing line 1 and 2 its shows below error.

                    arcpy.CopyRaster_management(maskdata, outname + ".tif", "", "", "", "NONE", "NONE", "")

                    print outLoc, outname

Error Messgae :

Traceback (most recent call last):
  File "<module1>", line 40, in <module>
  File "C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 7051, in ExtractByMask
    in_mask_data)
  File "C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Utils.py", line 53, in swapper
    result = wrapper(*args, **kwargs)
  File "C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\sa\Functions.py", line 7047, in Wrapper
    out_raster)
  File "C:\Program Files\ArcGIS\Desktop10.3\ArcPy\arcpy\geoprocessing\_base.py", line 504, in <lambda>
    return lambda *args: val(*gp_fixargs(args, True))
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000860: Input raster: is not the type of Composite Geodataset, or does not exist.
Failed to execute (ExtractByMask).


texts: Convert feature to raster, not just shape using ArcGIS Desktop?

texts:
I am trying to convert a shapefile into a raster, but ArcGIS gives values only to outside cells, not inside´s, so I get a raster with a hole in the inside with the shape of the feature.

How can I inverse the results?


texts: Reprojecting Natural Earth dataset to World Eckert I

texts:
My question is connected with QGIS 2. Could anybody explain me how I ought to reproject the Natural Earth dataset there in order my world map to look correctly? When I strive to reproject these data I get the result that you can see on the image below.

I was using the next parameters for my projection:

+proj=eck1 +lon_0=120 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

That is a modified "World Eckert I" projection.


texts: How to Select polygon feature in ArcMap and pass Selection into Python Addin?

texts:
I have a script that uses a arcpy's 'SelectByAttribute' syntax to choose a record in a polygon feature class. Following that, I do a 'SelectbyLocation' to select points from a feature class that lie within the selected polygon. I do some processing on the points afterwards.

What I want to know is, can I somehow use the 'Select Features' tool from the 'Tools' toolbar to select the polygon and then use that selection to pass into a python addin button instead of using the 'SelectByAttribute' noted above? I hope this makes sense. I'm new to python addins.


texts:
Can you tell me the procedure in ArcGIS 9.3 to create 3D lines which connect the 3D points using the elevation of the points.
For example, if pointA has 20 elevation and pointB has 30, then the line would connect pointA to pointB vertically such that one end point of the line would have 20 elevation and the other would have 30.
Or, i would like to know how to create a 3D line directly without using the points, where i can define the x,y,z of the end points.


texts:
I have a line vector file and want to create raster file from it.
Here, the raster file's each pixel corresponds with shortest distance from the line.

How can I do it with QGIS ?


texts: How do I make my Cartodb info windows wider?

texts:
According to my read of andrewxhill (more readable version) I should be able to change the info window size  by adding CSS, like so:

  x

Non-valid picture URL

{{description}}

description

{{description}}

{{name}}

But this has no impact on the width of my info window.


texts: External Web Map Server (Tile WMS) for geoserver

texts:
Can I add to geoserver Tiled WMS layer?

For example Tile WMS:

https://example-server.maps.eu/layer1/11-1115-700


texts: How do I set a target layer in the Editor Toolbar?

texts: I am using ArcMap 10.4 and have not been able to figure out how to set a target layer.


texts: Produce a filtered version of a GRIB file

texts:
How can you produce a filtered GRIB file containing a subset of information from an original "bigger" GRIB file?

Let's say I have a GRIB file covering pressure, temperature and wind direction for a large area and for 5 days. I'd like to be able to produce a file with only the data matching: {pressure, days 1 and 2, a geo-rectangle}. How could I approach this?

I found some information about GDAL and specifically gdal_translate. Looks like I could use commands to achieve this but I'm struggling to digest gdal's documentation.

One more option I'm considering is to build my own library (it would be in C#) to parse and write GRIB files, but it seems like a lot of work.

How can this problem be tackled?


texts:
I'm using OpenLayers - 2.11,GeoServer 2.13, ExtJS 3.4, and Java for developing GIS based web applications.

I have added some point features in one vector layer.
It has id, name, address, etc.

Can I filter some points with 'name' from that vector layer and only show matched features in vector layer?
I have to filter some points and later I have to show all hidden data.

eg:


texts: How to Batch Import Collada (.dae) files to Esri MultiPatch format?

texts:
I was wondering if anyone has come across a tool, or can point to code snippets, for batch importing Collada (.dae) files into a MultiPatch format.  I have found a good "How to", but it is all manual.  Also with this manual process the Collada file is not imported into is exact geospatial location.

Here is the youtube link to the "How to"
http://www.youtube.com/watch?v=udhm3yxVinc


texts: OGR Resetting SetAttributeFilter and SetSpatialFilter

texts:
This is more of a behavior question about OGR. When I set the spatial and attribute filters like so:
    self.ogr_layer.SetSpatialFilter(clip_geom)
    self.ogr_layer.SetAttributeFilter("MyField = 'some value")

How do I UNset these filters if I want to loop over the features again? I'm looking for clues in the docs but can't quite figure it out.
Do I need to close and re-open the ogr.Layer?


texts: osm2pgsql stops after importing several extracts

texts:
I want to import several city extracts from the bbike website into my postgis database. I use windows on a desktop pc with 16GB RAM and a quadcore processor. The first couple of extracts imported fairly well. I used this command:

osm2pgsql --append -s --database <database>--username <username> <path to pbf>

After around 10 extracts the import looks like it stopped (Speed 0.0k) without any error message or it is extremely slow. When I add the cache flag and set it to 12000 it has now effect.

Looking at pgadmin there still seems to be activity on the db. This is what my pgAdmin panel looks like:

Did anyone experience the same and got a solution?
What could be the reason for this?


texts: Spanish municipality level shapefile and population data

texts:
I'm struggling with the Spanish municipality level.

I need to incorporate population data with the lat/long for each Spanish municipality and below are the details of my current files:

I cannot find one master CSV having population for each municipality,
instead I have 52 files having population for municipalities within
each Spanish province;
I have the lat/long for each municipality in a separate CSV, and
The problem is, I need a master CSV with all data in one place and my CSV files don't have an id I can use to map the data.


texts: Closest point to road intersections with GeoPandas and Shapely

texts:
I have a geopandas dataframe, where each record's geometry field is made up of a LineString.  Each LineString represents the points which make up a road.  I'd like to find the locations of all road intersections from this data.  Currently, I have nested for-loops which take each record and pairwise compute the intersection of the two LineString objects.  My questions are:

Is this the best way to go about doing this?  Even 10 pair-wise computations is yielding ~78000 intersection points.  Perhaps I need to store in a spatial database like rtree if there are indeed that many intersections, but this seems inefficient...
My ultimate goal is, I have spatial points, for which I'd like to find the closest intersection point.  If there is a better or more efficient method than first finding all intersections and then finding the closest point, I'd love to hear suggestions.


texts: GDAL to project from pixels to geo coordinates

texts:
I am using GDAL library in C++. On my dataset, I have obtained geo transform with GetGeoTransform to GT array (GCP is NULL).

Than I use this values for this equation (Source):

Xgeo = GT(0) + Xpixel*GT(1) + Yline*GT(2)
Ygeo = GT(3) + Xpixel*GT(4) + Yline*GT(5)

I also add 0.5 offset to Xpixel and Yline, to calculate coordinate in the "center" of pixel.

If I put [0,0], I got correct result in lat/lon. However, if I change Yline to something other than 0, my longitude is messed up - instead of 44.0, I got -349.2... (latitude is still correct). What am I doing wrong ?

My data informations  -
origin, pixel size and top left are taken from GT array used in equation above; projection is GetProjectionRef

Projection is `GEOGCS["Coordinate System imported from GRIB file",DATUM["unknown",SPHEROID["Sphere",6371200,0]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433]]'
Origin = (0,0) - GT(2), GT(4)
Pixel Size = (0.5,-33.268) - GT(1), GT(5)
Top Left = (5.75,66.634) - GT(0), GT(3)


texts: How to zoom to a selected feature with pyQGIS

texts:
I want to create a function that selects a feature and zooms to it (similar in QGIS). Therefore there is following function:

QgsMapLayerRegistry.instance().addMapLayer(self.vlayer)

def zoomTo(self):
    layer = self.vlayer
    atable = self.ui.table

    selectList=[]
    for i in atable.selectionModel().selectedRows():
        ID = atable.item(i.row(),0).text()
        selectList.append(int(ID))

    layer.setSelectedFeatures(selectList)

The selected features are highlighted on the map. But I have no idea how to make a "zoom" to the selected features or some kind of focus them in the middle of the map.


texts: How to use Prolong Line feature with an input layer in WGS84?

texts:
I have a source point table in SRID:4326, and to make it look geographically accurate, I project my project file to SRID:2877 (State Plane CO) (layer still reports source CRS of SRID:4326).

In order to use the Prolong Line feature in Easy Custom Labeling, the project and layer CRS's must match.

But when I set the point layer CRS to SRID:2877 (overriding SRID:4326) the layer moves way over to the middle of Arizona (should be in Denver).

This happens with any of the State Plane CRS's for Colorado Central in the list.

In addition to drawing in the wrong place, it also greatly affects the 'map unit' size I'm using for the symbols and points in my map.

I've tried various datasets from my database and from the City of Denver, with all the same results.

Is there a setting I'm missing to do this accurately?


texts: Is it possible to run the arcpy module in the Python for SPSS plugin to run the search and update cursors?


texts: n:m relations of lines including directions

texts:
For trail planning and route signalization projects, I am looking for a solutions which allows me to connect "physical" trail segments (available as .gdb from authorities) with overlaying routes, whereas the routes may follow the segments in both directions.

The routes should be stored in a Postgres DB.

I assume this is a n:m relation, but how should the directions be manageed?
Routes normally consist of >1 segment and there are likely to be more than one route on one segment.

I use QGIS 2.14 on Mac OSX, postgres 11.


texts: Is there a quick way to calculate the change in an attribute over distance?

texts:
I am working in QGIS 3.14.0 and have a set of point data input via latitude and longitude in a CSV file; this file also includes data taken at these points. I would like to visualize the change in this data from point to point as a sort of slope (i.e., since it's nutrient data, micromolar per meter)
Currently I am using distance matrix to get the distance from point to point in meters, field calculator to get the slope, multipoint to point, then point to path all in order to display this in the way that I want, but I am wondering if there is a quicker way to do this. I have uploaded some images just as an example of my data here. The lines in my images are what I want them to look like, however, when I go through the steps outlined above, the lines only go from point 1 to 1 and 2 to 2 and so on, rather than 1-2 2-3, etc.


texts: Can I create dynamic style functions in OpenLayers with class breaks obtained at run time?

texts:
I have a GeoJSON object with ~1000 polygons and a number of numeric attributes on each polygon. I want to be able to style the polygons based on these attributes using dynamically calculated class breaks. I already have a method for calculating the class breaks.

I can create style functions for the polygon layer if I know the class breaks ahead of time, eg:

var bg_style = function(feature, resolution){
    var medinc = feature.get('medinc');
    if (medinc<36000) {
        fillColor='#ffffcc';
    }
    else if (medinc>=36000 && medinc<62000){
        fillColor='#c2e699';
    }
    else if (medinc>=62000 && medinc<92000){
        fillColor='#78c679';
    }
    else if (medinc>=92000 && medinc<140000){
        fillColor='#31a354';
    }
    else if (medinc>=140000){
        fillColor='#006837';
    }
    style = new ol.style.Style({
    fill: new ol.style.Fill({
      color: fillColor
    }),
    stroke: new ol.style.Stroke({
      color: [0,0,0,1],
      width: 0.5
    })
  });
    return style;
};

but if I don't know the class breaks ahead of time (they depend on the data), I'm at a bit of a loss as to how to use the dynamically calculated breaks in the style function.

It seems like if I want to use variables for the class max/mins & color hex codes, they're out of scope since the style function is not taking any arguments beyond (feature, resolution). I could possibly attach them to the data as additional attributes, but this seems sleazy and creates a few thousand pieces of redundant information.

Basically what I want to do is replace the min/max and hex color for each range in the code snippet above with the class breaks generated at run time. Also no guarantee the number of classes is fixed. Any suggestions???


texts: Postgis 2.0 - How to get a Point from given LinestringM and a measure?

texts:
If I have a cable route linear referenced and stored in a LinestringM with M value is the length (in meters) from the start of the route - for example:

LINESTRING M (0 0 0, 10 0 20, 12 0 40, 20 0 50, 21 0 70)

Now I know the cable route is broken at 45 meters from start point, that would be in the middle of two points (12 0 40) and (20 0 50) (see the M value). With Postgis, how can I get the broken point coordinate - which should be POINT (16 0) in this case?


texts: Is there any free source of Digital elevation model of China with the resolution of 30 meters or better and vertical accuracy of +-5m. Something along those lines


texts: Tearing of SpatialLinesDataFrame and problem with gClip function

texts:
I am trying to use ggmap with geom_map and geom_path commands to add a SpatialPolygonsDataFrame of African countries and a SpatialLinesDataFrame over a raster.  When I use the defined function gClip() which I gathered from the referenced question, I get the error from ggplot

  ggplot2 doesn't know how to deal with data of class SpatialLines

If I don't try to clip the SpatialLinesDataFrame, the precipitation lines plot on the map, but there is the expected "tearing" in places where contours continue outside of the bounding box.

I tried a few possibilities in modifying gClip() but can't figure it out.

Tearing of polygons using ggmap and readOGR

http://www.filedropper.com/archive_129

rm(list=ls(all=TRUE))

library(raster)
library(rgdal)
library(rgeos)
library(sp)
library(maptools)
library(ggmap)

setwd("~/folder")

africa_countries <- readOGR(dsn="shapefiles", layer="AfricanCountries",
                            stringsAsFactors=FALSE, verbose=FALSE)
precip_contours <- readOGR(dsn="shapefiles", layer="precip_contours",
                           stringsAsFactors=FALSE, verbose=FALSE)

precip_contours
class       : SpatialLinesDataFrame

gClip <- function(shp, bb){
  if(class(bb) == "matrix") b_poly <- as(extent(as.vector(t(bb))), "SpatialPolygons")
  else b_poly <- as(extent(bb), "SpatialPolygons")
  proj4string(b_poly) <- proj4string(shp)
  gIntersection(shp, b_poly, byid = T)
}

precip_contours <- gClip(precip_contours, matrix(c(-19,-36,52,38), ncol=2))

africa_countries <- fortify(africa_countries)
precip_contours <- fortify(precip_contours)

precip_contours <- fortify(precip_contours)
Error: ggplot2 doesn't know how to deal with data of class SpatialLines

# plot

africa <- get_map(location = c(-19,-36,52,38), zoom = 4, maptype = "watercolor")
africa <- ggmap(africa)
africa <- africa + geom_map(data=africa_countries, map=africa_countries, aes(x=long, y=lat, map_id=id), color ="white", fill ="orangered4", alpha = .4, size = .2)
africa <- africa + geom_path(data=precip_contours, aes(x=long, y=lat), colour="black", lwd=0.15)
africa


texts: GeoWebCache is able to serve layers in the following three standards: TMS, WMTS and WMS-C. Is there any remarkable advantage in using any of them in particular within an OpenLayers framework?


texts: Preserving formatting of Scales & Legends in Templates

texts:
My QGIS templates don't seem to preserve the formatting of my scales or legends. Is there a way to save these so they don't have to be redone each time? (I'm doing a series of maps that will be the same scale and have the same legends.)
Thanks
TVZ


texts: minimum distance from the shore (how to retrieve the "shore vertices")

texts:
I am student of Statistics, this is my first time using GIS., I am sorry for the stupidity of the question but I browsed online and could not find a satisfactory answer.

I have some spatial data in the format (lat, long, y).
I have noticed that the sea proximity influences the statistics y, Is it possible (and if it is how can you do it) to interpolate the coast and then calculate a minimum distance between the polyline interpolating it and the measurement point?

I can do the interpolation on R if it is complicated on GIS, but how can I retrieve the data of the "vertices" of the shore?
I understand I have to use polygons, but how does that work do I have to do it manually or is there an automatic way?

Please do not say to read the guides, I have read quite enough and I am clearly missing something.

I have downloaded GRASS GIS and have an ArcGIS online account.
My data about the Finland coast come from here. country: Finland, subject: Administrative Areas.

Thanks

Irene


texts:
I am working on a asp.net c# application by using google map control. I am displaying 80+ markers on the lat/ long receceiving from database. Now I want to display a route covering all the markers. I have written a code for route generation but its showing route on very limited markers.
I have search on google and some where its mentioned that there is a restriction of 10 location route generation, if it is then how could I accomplish this task?
GoogleMap1.Directions.Add(new GoogleDirections()
                {
                    Origin = new Location(preLat, preLong),
                    Destination = new Location(Convert.ToDouble(geo.Lat.ToString()), Convert.ToDouble(geo.Long.ToString())),
                    Draggable = false,
                    SuppressMarkers = true,
                    PolylineOptions = { StrokeColor = "red", StrokeWeight = 3 },

                });


texts: Python shortcut to find the distance between two coordinates points from different zones

texts: I am trying to find the distance between two points on different coordinate zones. I am referring to this Measuring distances when crossing UTM zones? and don't understand anything from it.


texts: I just updated to ArcGIS 10.3, and lost the ability to use all of the tools that are dependent on python.  I do have "C:Python27>ArcGIS10.3>python.exe".  The Python27 folder also had the ArcGIS10.2 folder - I deleted it and restarted my comp to be sure it was gone in case Arc was looking in the wrong place.  Could someone provide instructions on how to tell Arc to recognize that Python is present?


texts: Is there an ArcGIS 10 Amazon Machine Image (AMI) patched with SP4?

texts:
I know you can patch the SP3 AMI with SP4, but is there a SP4 AMI for ArcGIS server 10 and if it exists, what is the AMI Name (or AMI ID).

I know Esri has to specifically share an AMI with my Amazon account before it is discoverable in an AMI search.  But you wouldn't know to ask for it.


texts:
I am doing this tutorial
http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#/How_to_geocode_a_single_address/00010000038r000000/

but it requires the following file: UnitedStates.gdb

I downloaded DataMaps10 from ESRI but I do not see any .gdb files.


texts: Display label with both direction symbols

texts:
I have a line layer describing a road network, with traffic attribute values for both driving directions.
I would like to display label for both traffic values with direction symbols (above label is traffic in digit direction) :

   1500 >
----------->
   < 2200

Direction symbols option works perfectly but displays only one direction symbols :

   1500
----------->
   2200 >

I managed do display two direction symbols using an expression :

'<' || "dir_1" || '\n' || "dir_2" || '>'

but it's correct or not depending on the line orientation :

             1500 >                 1500 >
correct : <-----------   wrong : ----------->
             < 2200                 < 2200

Is there a workaround ?


texts:
I am having around million features stored in PostgreSQl table (with PostGIS extension).

I can generate vector tiles by performing below steps

Generate GeoJSON  from PG table using GDAL
Generate vector files in folder format using tippecanoe
Sync the features to S3

Then using Mapbox gl js layer is rendered in map

Now we are trying to add editing functionality. so that users can add, edit or delete features.

Is there a way to efficiently generate vector tiles with new created features using tippecanoe and keep the vector tiles up to date?

Is tippecanoe is right tool for this job?


texts: Installing GeoTrace plugin in QGIS on Mac OSX?

texts:
I am trying to install the experimental plugin called GeoTrace in QGIS 3.2 on Mac OSX.

Error message is as follows:

File "/Users/petergriffin/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/GeoTrace/tools/install_dependencies.py", line 116, in install
  import gttracetool
File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 665, in _import
  mod = _builtin_import(name, globals, locals, fromlist, level)
File "/Users/petergriffin/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/GeoTrace/tools/gttracetool.py", line 40, in
  import gttrace as trace
File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 665, in _import
  mod = _builtin_import(name, globals, locals, fromlist, level)
File "/Users/petergriffin/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/GeoTrace/tools/gttrace.py", line 33, in
  from skimage.graph import route_through_array
File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 665, in _import
  mod = _builtin_import(name, globals, locals, fromlist, level)
 ModuleNotFoundError: No module named 'skimage'

During handling of the above exception, another exception occurred:

Traceback (most recent call last):
File "/Users/petergriffin/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/GeoTrace/geo_trace.py", line 188, in open_trace
  self.dlg = GeoTraceDialog(self.iface)
File "/Users/petergriffin/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/GeoTrace/geo_trace_dialog.py", line 62, in __init__
  self.setup_gui()
File "/Users/petergriffin/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/GeoTrace/geo_trace_dialog.py", line 85, in setup_gui
  installer.Installer().install() #lol
File "/Users/petergriffin/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/GeoTrace/tools/install_dependencies.py", line 119, in install
  assert False, "An unexpected import error occurred"
  AssertionError: An unexpected import error occurred


texts: When I setup my layer properties I selected some of the attributes to be hidden on the Fields tab. I would like to change the EDIT WIDGET from hidden to a different function however some of them are no longer available in the layer properties dialog box. How do I make them available in the layer properties dialog box under the Fields tab again. Using version 1.7