texts: Creating multilayered classification map on Google Earth Engine?

texts:
I've been working on a classification scheme, and have 6 different classes i wish to map. Everything has been already done in terms of finding the classifications. the only thing I am struggling on now is how to import all 6 classes (6 classes but 7 vars) onto the same map, without them going on top of one another.

My end goal is for it to look like this:

Code used:

var image = ee.Image(landsat.filterDate ('2018-01-01', '2018-03-01').filterBounds(point).sort('CLOUD_COVER').first());

//print image
print('a landsat scene', image);

//define parameters

var green = image.select('B3');
var swi = image.select('B6');
var mndwi = green.subtract(swi).divide(green.add(swi)).rename('MNDWI');

var mndwipara = {min: 0, max: 0.7, palette: ['white', 'blue']};

Map.addLayer(mndwi, mndwipara, 'MNDWI');

//apply threshold to select only positive pixels
var lake_mask= mndwi.gt(0);
Map.addLayer(lake_mask);

//image reducer
var meanDictionary =image.reduceRegion({
reducer: ee.Reducer.mean(),
geometry: image.geometry(),
scale: 30,
maxPixels: 1e9
});
// The result is a Dictionary.  Print it.
print(meanDictionary);

//classification attempt - scum
var scum=image.expression(
'(nir)*0 + (nir < 0.018)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
});

 //classification attempt - sediment
var sediment=image.expression(
'(nir)*0 + (nir > 0.018 && red > green && green/nir < 1.3 && red/green> nir/red)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
});

//classification attempt - lowbiomass
var lowbiomass=image.expression(
'(nir)*0 + (nir > 0.018 && red < green && green > 0.065 && nir < red)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
});

//classification attempt-scum
var scum=image.expression(
'(nir)*0 + (nir > 0.018 && red < green && green > 0.065 && nir > red && nir/red > 4.07)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
  });

 //classifcation attempt bleached scum
  var bleached_scum=image.expression(
  '(nir)*0 + (nir > 0.018 && red > green && green/nir > 1.3 && green > 0.065)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
});

//classification attempt - microphytobenthos
var microphytobenthos=image.expression(
'(nir)*0 + (nir > 0.018 && red > green && green/nir < 1.3 && red/green< nir/red)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
});

//classification attempt - highbiomass
var highbiomass=image.expression(
'(nir)*0 + (nir > 0.018 && red < green && green > 0.065 && nir > red && nir/red < 4.07)*1', {
'nir': image.select('B5'),
'red': image.select ('B4'),
'green': image.select('B3'),
 });


texts: Can't remove the little islands when creating a svg map with kartograph.py

texts:
I'm new to Python and to Kartograph. I'm using it to create a map of France, Italy, Portugal and Spain. But I don't want the little islands (mainly the French ones), causing the generation of a big useless map, for me.

I used a Python script to create my map ( thanks to rmbianchi ):

# file: makeMap.py
from kartograph import Kartograph
K = Kartograph()

def myfilter(record):
    return record['iso_a3'] in ["FRA","ITA","ESP","PRT"]

config = {
    "layers": {
        "FTC": {
            "src": "ne_50m_admin_0_countries.shp",
            "filter": myfilter,
            "attributes": {"iso_a3":"iso_a3", "name":"name", "id":"iso_a3"},
        "bounds": {
    "mode": "polygons",
    "data": {
        "layer": "FTC",
        "filter": {
            "filter": ["ISO_A3", "in", ["FRA", "ITA", "ESP", "PRT"]]
        },
        "min-area": 0.50
    }
}
        }
    },
}

K.generate(config, outfile='ftc.svg')

But, little islands are still displayed, how can I avoid that?


texts: Converting .las files to LAS datasets using ModelBuilder iterator?

texts:
I have 1500 individual LAS files and I need to move them to separate LAS datasets (so 1500 LAS datasets).

I've used iterator (iterate files) in ModelBuilder and it creates 1500 individual LAS datasets with correct individual names but they all have points just from first .als file.

There is nothing wrong with .las files since I can open them in Global Mapper.

I am using ArcGIS 10.2.


texts: Geoserver read only error for updating wfs layer

texts:
I created a linestring wfs layer in postgis table and published to geoserver. I am able to display the layer on the openlayers map but while editing the layer getting error:

  {http://localhost:8090/geoserver/postgis_24_sample}sampleLine is
  read-only

below is my code for posting the edit data. Please suggest me how to clear this issue.

var transactWFS = function(mode, f) {
          var gml = new ol.format.GML({
            featureNS: 'http://localhost:8090/geoserver',
            featureType: 'sampleLine',
            srsName: 'EPSG:3857'
         });
         var node;
         switch (mode){
            case 'insert':
                node = formatWFS.writeTransaction([f], null, null, gml);
                break;
            case 'update':
                node = formatWFS.writeTransaction(null, [f], null, gml);
                break;
            case 'delete':
            node = formatWFS.writeTransaction(null, null, [f], gml);
            break;
          }
        var payload = xs.serializeToString(node);
        $.ajax('http://localhost:8090/geoserver/postgis_24_sample/wfs', {
            type: 'POST',
            dataType: 'xml',
            processData: false,
            contentType: 'text/xml',
            data: payload
            }).done(function() {
                mySampleLineVectorLayer.getSource().clear();
            });
        }


texts: Chart inside pop up window in Javascript

texts:
How to show a chart inside pop up window? I'm using Leaflet JS, which has support for raphael plugin http://dynmeth.github.com/RaphaelLayer/. Is it possible to create a div inside popup window? I was thinking about http://softwarebyjosh.com/raphy-charts/

If it's not possible with Leaflet, I am opened to only Raphael map solution.

Thanks


texts:
I have my raster layers of Digital Elevation Model, aspect, slope and curvature for a certain region. In this region I want to apply a grid with a certain spacing in its nodes, and then get the corresponding values of DEM, aspect, slope, curvature for the grid nodes sites. I have tried in ArcMap 10.2 to apply the fishnet, but it doesn't seem to give me this option (get raster values exactly on grid nodes).

Any idea on how can I perform this in ArcMap or R?


texts: Spatial data? Geodata? Geographic Data? Geospatial data?

texts:
What word indicates spatial data within the geographical and geological domains? In particular, what word is precise enough for a science journal article (i.e., formal prose)?
Candidates might include "spatial data," "geodata," "geospatial data," and "geographic data." "Spatial data" seems too generic (not necessarily within the geographical and geological domains).
What is the consensus, if any exists?


texts:
I have tried to apply styles (stroke,fill..) in a vector layers(kml,geojson) but I never see changes in the map. This is my code:

<!DOCTYPE html>
<html>
<head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">

<script src="https://openlayers.org/en/v3.20.1/build/ol.js" type="text/javascript"></script>

<link rel="stylesheet" href="https://openlayers.org/en/v3.20.1/css/ol.css" type="text/css">

<style>

.map {
    height: 510px;
    width: 100%;
}

</style>

</head>
<body>

<div id="map" class="map">
</div>
<script type="text/javascript">

_myStroke = new ol.style.Stroke({
     color : 'rgba(255,0,0,1.0)',
     width : 1
});

_myFill = new ol.style.Fill({
     color: 'rgba(255,0,0,1.0)'
});

myStyle = new ol.style.Style({
     stroke : _myStroke,
     fill : _myFill
 });

var raster1 = new ol.layer.Tile({
      source: new ol.source.Stamen({
      layer: 'toner-lite'
      })
});

var isocrona = new ol.layer.Vector({
      source: new ol.source.Vector({
          url: './datos/isocronas_json.geojson',
          format: new ol.format.GeoJSON(),
          extractStyles: false,
          style: myStyle
          })
      });

var map = new ol.Map({
      target: 'map',
      layers: [raster1, isocrona],
      view: new ol.View({
          projection: 'EPSG:4326',
          center: [-4.4328, 36.7277],
          zoom: 12

          }),

});

</script>
</body>
</html>


texts: WMTS - get tiles only for a specific area of interest

texts:
Multiple basemaps of Austria can be achieved via a WMTS Service provided by basemap.at. The GetCapabilities interface is https://www.basemap.at/wmts/1.0.0/WMTSCapabilities.xml

Working in a small area of interest and having all data loaded takes a while when zooming in and out frequently. Thus, I would like to get the basemap only for a specific region.

I know that it should be possible to set a bounding box or retrieve only specific tiles, but no matter what specifications I try, it does not work.  I do not get an error message, but it always loads the basemap of whole Austria. I work within ArcMap and connect to the service via the GIS Servers - Add WMTS server.

I found this question (Getting a WMS-like bounding box from WMTS services) and tried to make the request similar.

For the layer I'm interested in, I have the following information:
Tile matrix sets is google3857, Identifier is bmaporthofoto30cm, Bounding Box (WGS84) is 8,782379, 46,358770, 17,500000, 49,037872.
But I am unsure which value to set for the parameter TileMatrix.

That is my modified request: https://www.basemap.at/wmts/1.0.0/WMTSCapabilities.xml?service=WMTS&request=GetTile&version=1.0.0&layer=bmaporthofoto30cm&style=default&format=image/png&TileMatrixSet=google3857&TileMatrix=0&TileRow=1&TileCol=1

Any idea why it does not work and what I can try?


texts: CartoDB - "sql applied has errors" message - where can I see the error text?

texts:
I'm still getting my bearings with the interface it seems...

When running a SQL query in CartoDB against a table... it seems like a get the correct result, but I also get a see a "sql applied has errors" message when I hover over the SQL icon which has both a red and a green dot displayed on it.  Is there any way to see the particular error message?  I don't see a place to get that info...


texts: How much do leaves on/off affect the number of LiDAR ground returns?

texts:
We would like to collect a generic purpose topographic LiDAR dataset over a large region in Belgium. Our minimum requirement is to achieve at least 1 point/m² on the ground. Ideally, we would like to have images from spring or summer because it would then be possible to get some extra information about shrubs AND more chances to have good flight conditions, but this is not our priority (the density of soil return must be sufficient). Alternatively, we could work during winter in order to have less vegetation cover in crop fields and deciduous forest, but then we would have less information on understory shrubs. This decision is difficult to take without knowing how much "ground return" we would lose when the vegetation cover is at its maximum.

Therefore, I am wondering if there are some shared experience about the proportion of ground returns under deciduous forest with or without leaves. I guess that this proportion is also depends on the LiDAR footprint, but I could not find any comparative studies.


texts: Obtaining ArcView 3.2 Spatial Analysis 1.1?

texts: Is it still possible to obtain and use a copy of the Spatial Analyst 1.1 extension to ArcView 3.2 in order to finish a final project at a University?


texts: Drag marker in "repeated" map, wrong latlng

texts:
When you render a map if you zoom out far enough it will repeat the map infinite times.

The problem is that the coordinates will increase if you drag something outside the "original" map.

Example:

and if you drag to the left

The 'lat' seems to be fine, since it only duplicates horizontally, but 'lng' is way off.

Is there a way to correct the position?

L.vesion = "1.0.0-rc.3+27f05a2"


texts:
Encountering a GeoServer / OpenLayers WFS layer not visible issue, I found this interesting answer. Solution about HTML page on the server surprised me very much as I tried previously to load an HTML page from the Jetty GeoServer without any success, encountering always an error:

HTTP ERROR: 404 NOT_FOUND RequestURI=/web/ws3.html

Of course, I tried several other places on my http://localhost:8080/geoserver.

Can someone confirm that its possible to load any HTML page from GeoServer, and if a special configuration is required to reach this status?


texts: How do you line up a pre-existing public domain map with csv coordinate points on QGIS?

texts: I have been struggling with this simple problem all day. I found a Florida map with the counties on Google for GIS. I uploaded the map to QGIS with no problem. However, whenever I try to upload points (both latitude/longitude and x,y coordinates) of areas in Florida, the points do not line up with the map. I have tried different coordinates and I have experimented with the CRS of the maps and I have lined everything up with the same references, but I still have no luck. Has anyone ever tried to do this before?


texts: Sieve in QGIS returns unexpected results

texts:
I have a rasterimage with values between 1 and 5.

I want to sieve it, using QGIS - Raster - Analysis - Sieve.

But then I get this:

The values changed. Even worse, when I color it, the whole raster image is one value.
The weird part, I created the image with an R script. When I sieve another output of the same R Script, it works...

What is going wrong?


texts: Can QGIS be used to calculate the up-slopes (total area draining to a specified point) of multiple separate points simultaneously?

texts: I am trying to calculate the area draining to a point of 240 separate points using QGIS. At the moment I'm using the Saga - Upslope Area algorithm to individually work this out for each of the 240 points using their coordinates and a DEM input. This is very time-consuming so is there is a quicker more efficient algorithm which would allow me to input the whole vector file for the point data to simultaneously output the upslopes for all the points at once?


texts: Shortest path to draw all polyline routes in geometric network without repeating

texts:
My problem is related to the "shortest path" problem, however it's slightly different in that I don't want to draw one path but to draw all paths without repeating any routes (as much as is possible). I'm trying to work out what this geometric problem is called (if it has been covered at all) in order to find an algorithm to work with.

Say there are 6 points in a network that are attached via a web of 10 edges (non-weighted). I must make a route starting form any one single point and move down every edge in any order but without repeating any of these edges. When I find myself in a position where there are no possible new edges forward I can either take one of the previously travelled edges back or I can go "off-road" directly to a new point, however these are considered "lost time". What I want is to find an algorithm (or even the name of this geometric problem) that will minimise "lost time" As I attempt to move down all routes.

An easier way to explain might be a printer that draws out all of these paths: while the printers pen is down (i.e. drawing) we are not wasting energy, however as the printers pen will not repeat any drawn lines it may have to move between points with the pen up, which we consider as wasted energy. I want to minimise energy waste.

Can anyone help point me in the right direction?


texts: ArcMap, creating lines from a set of points that connect to the other line

texts:
I have a bunch of points lining up on one side of the highway.

I would like to create line from each point that reaches to the other side of the highway. The line should be (almost) perpendicular to both sides of the roads. Or, at least it connects the shortest distance between highways. I would also like to let the newly created lines to have the same attribute fields of the points.

Can anyone give me some inputs?


texts: SQL queries in PostGIS database layer using Python and QGIS

texts:
I want to create some layers queries in my QGIS plugin using Python.
I have find here a question and I try to follow with success (works fine for a shapefile).

from osgeo import ogr

filepath = 'C:/Users/username/Desktop/path/clip.shp' # set the filepath
layer_name = filepath[:-4].split('/')[-1] # get the layer name
driver = ogr.GetDriverByName ("ESRI Shapefile")
ogr_ds = driver.Open(filepath)
TEST=2
sql = "SELECT myfield FROM %s WHERE OBJECTID=%s" %(layer_name,TEST)
layer = ogr_ds.ExecuteSQL(sql)
feat = layer.GetNextFeature()
val = feat.GetField(0)

print (val.decode('utf-8'))

but that question is for shapefile what I need more if that is layer in PostGIS database.

I know only the connection like this:

import psycopg2

try:
    conn = psycopg2.connect("dbname='template1' user='dbuser' host='localhost' password='dbpass'")
except:
    print "I am unable to connect to the database"


texts: Splitting up a single polygon into many polygons of equal point count

texts: I have point data split into 15 different values and want to split up the (green) polygon into a number of smaller polygons each containing the same number of points so that I can make a chloropleth map showing which sort of value is dominant in each area. I'd also like to maintain the structure of the islands whilst I'm doing that. Was thinking something along the lines of finding the nearest x neighbors of each point and then clustering from that and then editing with the Cut Polygons Tool but was wondering if there were any tools / script that did it more efficiently?


texts: How can I get a single part dissolve using ST_Union or ST_Aggr_Union?

texts:
I am trying to get dissolve functionality using ST_Union using this query :

SELECT c.fid, ST_Union(c.boundaryshape) FROM c GROUP BY c.fid;

Than I found instead of using ST_Union I should use ST_Aggr_Union :

SELECT c.fid, ST_Aggr_Union(c.boundaryshape) FROM c GROUP BY c.fid;

I am getting results like dissolve but it looks like this is multi-part.

Can somebody let me know how can I get the single-part dissolve?


texts:
I have this sample of a tree with the table 't_cheminement' :

My schema is a table with : letter, startPointNumber(cm_ndcode1), endPointNumber(cm_ndcode2) and the geom.

I have some issues when I try to extract all ends of my tree : C,G,E,I

Some end are well oriented but sometimes not. The blue arrow is the same way as start and point numbers, they aren't inverted.

I try to use ST_Startpoint and ST_Endpoint :

select t1.cm_code, t1.cm_ndcode1, t1.cm_ndcode2 from t_cheminement as t1
inner join t_cheminement as t2 on ST_Equals(ST_StartPoint(t1.geom), ST_StartPoint(t2.geom)) is false
inner join t_cheminement as t3 on ST_Equals(ST_StartPoint(t1.geom), ST_EndPoint(t3.geom)) is false
inner join t_cheminement as t4 on ST_Equals(ST_EndPoint(t1.geom), ST_EndPoint(t4.geom)) is true or ST_Equals(ST_EndPoint(t1.geom),
ST_StartPoint(t4.geom)) is true
UNION
select t1.cm_code, t1.cm_ndcode1, t1.cm_ndcode2 from t_cheminement as t1
inner join t_cheminement as t2 on ST_Equals(ST_EndPoint(t1.geom), ST_StartPoint(t2.geom)) is false
inner join t_cheminement as t3 on ST_Equals(ST_EndPoint(t1.geom),
ST_EndPoint(t3.geom)) is false
inner join t_cheminement as t4 on ST_Equals(ST_StartPoint(t1.geom),
ST_EndPoint(t4.geom)) is true or ST_Equals(ST_StartPoint(t1.geom),
ST_StartPoint(t4.geom)) is true

I also tried to use start and endpoint numbers :

select t1.cm_code, t1.cm_ndcode1, t1.cm_ndcode2 from t_cheminement as t1
inner join t_cheminement as t2 on t1.cm_ndcode1 <> t2.cm_ndcode1
inner join t_cheminement as t3 on t1.cm_ndcode1 <> t3.cm_ndcode2
inner join t_cheminement as t4 on t1.cm_ndcode2 = t4.cm_ndcode2 or t1.cm_ndcode2 = t4.cm_ndcode1
UNION
select t1.cm_code, t1.cm_ndcode1, t1.cm_ndcode2 from t_cheminement as t1
inner join t_cheminement as t2 on t1.cm_ndcode2 <> t2.cm_ndcode1
inner join t_cheminement as t3 on t1.cm_ndcode2 <> t3.cm_ndcode2
inner join t_cheminement as t4 on t1.cm_ndcode1 = t4.cm_ndcode2 or t1.cm_ndcode1 = t4.cm_ndcode1

Is it possible to help me to fix my query because it's taking a lot of time to execute and maybe I'm wrong about all 'inner join'.


texts: Is it possible to use Collector and Survey123 to generate new fields and append to existing records simultaneously?

texts:
What I am trying to do is make a webmap for Collector that field staff will use to collect a single piece of information at the same locations monthly.  In this case they are gathering groundwater elevation.  I was able to set up the webmap for Collector and link a survey from Survey123.  However, this creates a whole new record each time.  I was hoping to make the survey add {month_year} and they could add their measurement.
I think I am overreaching the capabilities of Survey123/Collector.  My backup is just using Collector and adding a new field prior to the fieldwork each month, unless someone has a better suggestion.


texts: WPS process importing and returning output

texts:
I'd like to import to the GeoServer catalog the result of one of my processes (in addition to returning it to the client), which is a shapefile.
I have managed to find some kind of import code, but it doesn't publish the layer. I also have to return a URL that links to that newly imported shapefile.

My code so far : (I have simplified a little bit, I understand that the way it is presented here (with a direct path) is useless.)

@DescribeProcess(title = "AddSHP", description = "Returns a SHP")
public class AddSHP implements GSProcess {

private Catalog catalog;

public AddSHP(Catalog catalog) {
    this.catalog = catalog;
}

@DescribeResult(name="ShapeFile", description="ShapeFile", type = String.class)
public RETURN_TYPE execute(
        throws IOException {
    String shp;
    shp = "Path/to/my/shape.shp";
    String shpname = shp.substring(shp.lastIndexOf('/') + 1);
    addShapeFileDataStore(shpname, shpname, shp);
    return (SOMETHING);
}

private boolean addShapeFileDataStore(String title, String description, String pathToFile)  {
    DataStoreInfo dsInfo = catalog.getFactory().createDataStore();
    dsInfo.setName(title);
    dsInfo.setDescription(description);
    dsInfo.setEnabled(true);
    dsInfo.setType("Shapefile");
    dsInfo.getConnectionParameters().put("create spatial index", true);
    dsInfo.getConnectionParameters().put("charset", "ISO-8859-1");
    dsInfo.getConnectionParameters().put("filetype", "shapefile");
    dsInfo.getConnectionParameters().put("cache and reuse memory maps", true);
    dsInfo.getConnectionParameters().put("url", "file:" + pathToFile);
    dsInfo.getConnectionParameters().put("namespace", "http://");
    try {
        catalog.add(dsInfo);
    } catch (Exception e) {
        System.out.print("Error adding data store to catalog");
    }
    return true;
}
}

Could anybody help ?
Thanks


texts: Trying to understand updateCursor() for python

texts:
I am using ArcMap 10.2.2. I am attempting to write a script that will cycle through the classes in a database and change the attribution of 1 row.

Is my variable "lyr" being recognized as the classes in the dataset?

When referencing an attribute to change, do you reference the alias or the actual name of the attribute?

Code sample:

import arcpy

dataset = arcpy.GetParameterAsText(0)
arcpy.env.worspace = dataset
classList = arcpy.ListFeatureClasses()

for lyr in classList:
     rows = arcpy.UpdateCursor(lyr)
     for row in rows:
         othField = row.getValue("OTH")
         if othField == "No Information":
             row.OTH = "noInformation"
             rows.updateRow(row)
         else:
             pass
del row
del rows


texts: Calculating mean value of polygon from raster in PostGIS?

texts:
I started with a NetCDF .gri and .grd raster file of the UK provided by a colleague. I clipped it in R to be only London, exported and converted to a ASC file, and then imported it to PostGIS using the following commands in R:

library(raster)
uk_raster <- raster("AnnMean2011.grd")
london_area <- extent(-720000.0,-630000.0,-50000.0,25000)
london_raster  <- crop(uk_raster, london_area)
writeRaster(london_raster, filename="AnnMean2011.asc", format="ascii")

And then on the Ubuntu command line:

raster2pgsql -I -C -s 10001 -t 20x20 AnnMean2011.asc annualmean | psql -d james_traffic

I now have a raster table in PostGIS. The SRID of 10001 is the following by the way, again, provided by a colleague:

INSERT INTO spatial_ref_sys(srid, auth_name, auth_srid, proj4text)
VALUES (10001,'CMAQ_Urban',10001,'+proj=lcc +a=6370000 +b=6370000 +lat_1=35 +lat_2=65 +lat_0=52 +lon_0=10 +x_0=000000 +y_0=00000');

In the same database I have a polygon file, SRID 27700, which covers London. I would like to calculate the average value within each polygon, from the raster.

I'm trying something like this, but it's not right:

select polygons.postcode, avg(st_value(joined_data.rast))
    from (
       select (ST_Intersection(raster.rast, 1, polygons.geom)).*
    from raster, polygons
       where ST_Intersects(raster.rast, 1, polygons.geom)
  ) joined_data
group by polygons.postcode

How would I go about this please?

It seems to get the polygon and the raster aligned properly, I have to convert them both to WGS84 I think.


texts:
I am developing a custom script to create something like the Split tool but in Basic License level.

So, at first, I Would like to create an Index Grid Feature Tool over my rectangle (this rectangle entire hover my raster) ,but.... The problem is that this Tool automatically DONT calculate the Polygon Width and Polygon Height.... But in arcMap 10.2 it DOES!? Why?

Have i problem with some configuration of my ArcGIS?

Is something changed in this tool from version 10.0 ?

Or can someone help me how can i get the right polygon width and height other way?


texts: Gradient fill based on scores (similar to elevation)

texts:
I'm rather new to using GIS (QGIS, specifically). I've mapped out a cohort of around 550 patients across Long Island, New York based on their full address by using MMQGIS plugin. This data is represented as "Single Symbol" "simple marker"s across the map. In the same layer, their 'PCL Score' (integer ranging from 17-85) is also assigned to each patient.

How can I either edit the current layer or add a new layer so that I have one continuous gradient of color showing how high or low the PCL trend is based on location? Basically, this is much like an elevation map where higher elevations are darker color and lower elevations are lighter color.


texts:
This is my first time using python in ArcGIS. I have a gdb of about 100 feature classes, and need to get Zonal Statistics for each. I'm using a textbook, and got this so far:

import arcpy
from arcpy import env
env.workspace = 'D:\Documents\Archaeology\Dissertation\Data\GIS\Places_Bounds\Places\Exploded_Bilinear\Argolid.gdb'
fclist = arcpy.ListFeatureClasses ("ID_*")
for fclist:
    arcpy.sa.ZonalStatistics ("fclist", "place_id", "pelop1arcm34b", "MEAN", "NODATA")

The error is coming from line 6. I am assuming my trouble is defining the input for the zonal statistics -- I want it to refer to the geodatabase set in the workspace, and loop through the whole dataset (about 100 feature classes). For reference: place_id is the unique identifier, and pelop1arcm34b is the DEM.


texts: Fill in the gaps from a ST_Union of multiple polygons

texts:
I have the following polygons:

01060000000300000001, 01060000000400000001, 01060000000400000001, 01060000000200000001, 01060000000100000001, 01060000000200000001, 01060000000100000001, 01060000000300000001, 01060000000100000001, 01060000000100000001, 01060000000100000001, 01060000000100000001

And I'd like to merge them using:

SELECT st_union(ARRAY(SELECT geom FROM areas))

And I get: https://pastebin.com/V18JtCz8 (I should paste here but it's above the char count - if you're reading this and the link doesn't work, it's simply the output of the above)

When I plot this output, there are holes (deliberate) but I'd like to just take the outline of the polygon as the new polygon.

How would I go about this?


texts: Iterating over a classification using different sample sizes in google earth engine zx

texts:
I am trying to loop through an RF classification and a confusion matrix using different sample data but I keep getting the following error before generating the matrix even though B1_median does exist: List (Error) Property 'B1_median' is missing.
var bands =['B1_median', 'B2_median', 'B3_median', 'B4_median','B5_median','B6_median','B7_median','B8_median','B8A_median','B9_median','B10_median','B11_median','B12_median','NDVI','EVI','Entropy','asm','contrast','corr','var','idm','savg','svar','sent','ent','dvar','dent','imcorr1','imcorr2','maxcorr','diss','inertia','shade','prom'];
var features = ee.List([Sample05, Sample10, Sample15, Sample20, Sample25, Sample30, Sample35, Sample40]);

var classifier = features.map(function(t)
  {
  return ee.Classifier.randomForest({numberOfTrees: 10})
  .train({
    features: t,
    classProperty: 'Class',
    inputProperties: bands })
  classified = kruis_m.select(bands).classify(t)
  print ('calssifi', classified)
  })


texts: Use of Postgis raster datasets for visualization and generation of GeoJSON

texts:
I am currently using Boundless/OpenGeo 4.02 with PostGIS 2.1. I would like a little advice on the development of specific functionalities. I am currently attempting to begin with a number of GeoTiffs, and would ultimately like to use them to for display in an OpenLayers map, and also provide the ability for users to draw a polygon (or enter coordinates) on this map and extract a GeoJSON file of pixel or point values from the underlying raster. My thinking is to import the geotiffs as PostGIS raster files, open them as a data source in OpenGeo server, and then display them in the OpenLayers map. Users can then query the map by drawing a polygon, which would query the PostGIS table and return the selected pixel values formatted as GeoJSON.

First, does this sound like a reasonable avenue to take given the required functionality?

Second, where can I find the most up to date documentation on correctly reading PostGIS rasters into OpenGeo server? I have seen references to Image Mosaic JDBC, but can't seem to find how I build in this capability into the version of OpenGeo (4.02) that I am currently using....


texts: Calculating distance between coordinates

texts: I would like to calculate the distances between a list of coordinates (classic GPS coord.). From coordinate A to B from B to C from C to D and so on, and so forth. Is there any  program, or website, or any QGIS plugin i could use? Sadly im not familiar with R. Any suggestions? Thanks!


texts:
I am trying to make a plugin for QGIS.
This plugin has lots of forms and based on the condition the form calls a child Form. The following is an example

I am using PyQt for creating forms. I want the Parrent and QGIS not be accessible(totally locked) while the CHILD is active.
and the user either can select or close the child.

Can someone guide me on how can it be done?


texts:
I'm using a WFS that is openly available as a demo from the MapServer web site.

http://demo.mapserver.org/cgi-bin/wfs?SERVICE=WFS&VERSION=1.0.0&REQUEST=GetCapabilities

And now I want to make a client application that can consume this service in C++ (QT) using the GDAL libraries.

I have compiled GDAL with Xerces and curl as the docs say but apart from that they are not very clear about how to use it for a WFS, I'm sure the functionality is there because it says it supports it but I can't find any WFS specific documentation, so far I have just tried to open a connection and get this error:

GDALAllRegister();
const char *pszFilename = NULL;
pszFilename = "http://demo.mapserver.org/cgi-bin/wfs?";
m_dataset = (GDALDataset*) GDALOpen(pszFilename, GA_ReadOnly);

if( m_dataset == NULL )
{
    std::cout << "Dataset is empty" << std::endl;
    exit( 1 );
}

Error:

ERROR 4: `/vsimem/http_1/file.dat' not recognised as a supported file format.

Dataset is empty
ERROR 4: Failed to create temporary file:/tmp/file.dat

I have also tried instead of the URL, an XML document that contains the URL like the docs say this should work.

<OGRWFSDataSource>
    <URL>http://path/to/WFS/service[?OPTIONAL_PARAMETER1=VALUE[&amp;OPTIONNAL_PARAMETER2=VALUE]]</URL>
</OGRWFSDataSource>

But this just says the file isn't a compatible type, any help would be appreciated here.


texts: Feature selection of one vector layer based on another vector layer: Any advice for using PostGIS or Python with PostGIS together?

texts:
I have a data set (more than 50,000+) of points (addresses) and Polyline (street) shapefiles of the area of interest. The street shapefile contains Annual Daily Traffic (ADT) data of each street for a long time period (1990 - 2010) for example: 1990, 1995, 2000, 2005 and so on.

I need to select streets within 30m buffer based on address locations with highest Annual Daily Traffic (ADT). Considering the complex scenario of streets (depicted by two examples below), the closest road to the address location can be a minor road or a major road but it must contain highest ADT value. If condition is met withing 30 m with highest ADT, then I need to shift my point feature to that street location otherwise the buffer is increased to 50m. This whole process needs to be iterated for all address locations (point features)enter image description here

I am looking for some advice so that I may come to know about a suitable tool i.e., Python, PostgreSQL and PostGIS etc. because ArcGIS model builder took a lot of time for this task. Any advice for using PostGIS spatial queries for this task or integrating Python with PostGIS would be more suitable?


texts:
I have a geodetic WGS84 position, lat, lon, alt, and direction cosines from the same point, uN, uE, uD. I'm brand new to geodetic/ecef conventions. How do I take a vector (say, a unit vector with magnitude 1), with direction uN, uE, uD, and represent it in ECEF?

I can find a lot of information online about going the reverse way: ECEF to NED, but nothing about NED to ECEF.

I was already able to convert my geodetic position to ECEF, but I'm pretty sure I'll need a direction cosine matrix to convert the NED vector to ECEF, right?

But isn't latitude geodetic?

Would I need to convert to geocentric to perform the basis transformation and create the direction cosine matrix for NED -> ECEF?

Can't for the life of me figure out how to represent vectors like these in X,Y,Z


texts: Creating queries to select by attribute with arcpy

texts:
I'm quite new with arcpy and I'm finding selecting by attribute a bit of a nightmare. There's a problem with my query on the script below [query = "\'Location\' = observer"], but I can't tell what it is. 'Location' is a field of the feature class I'm dealing with (fc_copy).

with arcpy.da.UpdateCursor(fc_copy,field_names) as cursor:
    for row in cursor:
        observer = row[0]
        arcpy.MakeFeatureLayer_management(fc_copy, 'fl_copy_OBS')
        arcpy.MakeFeatureLayer_management(fc_copy, 'fl_copy_TAR')
        query = "\'Location\' = observer"
        arcpy.SelectLayerByAttribute_management('fl_copy_OBS', 'NEW_SELECTION', query)
        arcpy.SelectLayerByAttribute_management('fl_copy_TAR', 'NEW_SELECTION', query, invert_where_clause = 'INVERT')

This is the error I get => Invalid expression
                           Failed to execute (SelectLayerByAttribute).

Is there any documentation to learn the rules to follow to build queries for arcpy?


texts: Apply Symbology from not identical layer

texts:
I've been using an ArcPy script for a little while now that applies the symbology ('unique values, many fields') from a template layer to a newly uploaded layer to avoid having to set the symbol, font, unicode, and size, every single time I upload a new version of the file.

This code has worked perfectly fine in the past. It creates 9 different values from 2 value fields. However, today I ran into a problem when the data I uploaded needed only 8 of the 9 values. The symbology didn't display at all and I had to manually change the symbology for each of the 8 values.

I was wondering if there was a way to modify the code so that it would check to see if a value needed to be created or not, and didn't create it if it wasn't needed, avoiding this error.

For example: If only 6 values were needed, the symbology for 6 values would be applied rather than all 9 (thus breaking the symbology for the layer).

I'm not sure if you need my code but I'll provide it anyway. It's  pretty simple:

addLayer = arcpy.mapping.Layer("symbology template layer")
arcpy.mapping.AddLayer(df, addLayer, "BOTTOM")

symbLayer = r"smybology template layer"

arcpy.ApplySymbologyFromLayer_management("layer that needs symbology", symbLayer)

for lyr in arcpy.mapping.ListLayers(mxd, "", df):
  if lyr.name.lower() == "symbology template layer":
    arcpy.mapping.RemoveLayer(df, lyr)


texts: Creating 3D geometry column by combining existing 2D geometry column with "elevation" column?

texts:
I want to create a table that contains 3D data.
I currently have 2D data stored as EWKT, which is multipoint geometry in one column, and another column that has elevation from sea level in another column.
Is there a way that I can combine these two columns and output a single column with 3D data?

Example:

Geom column ("010400000001000000010100000000001038011323410000E9631AFC5641")

Elevation column ("1346.9495")

My goal is to combine the Geom column with the Elevation column, to create a 3D geometry column.

My thought is to possibly use subquerries, but I don't know how.

I am using pgAdmin 4.


texts: Raster loaded in ArcMap 10.2.1 not showing decimal places for value

texts:
I have an IMG format raster that I am loading into ArcMap 10.2.1. The pixel type and depth, as reported by the Source tab in Layer Properties is 32-bit floating point. The values, as reported in the Symbology tab and Identify window are integers, where I was expecting non-zero values to the right of the decimal point. For example, ArcMap gives me a value of 21.000000, whereas if I open the same image in a different GIS package like QGIS, the value for the same pixel is listed as 21.345664.

Has anyone encountered this issue?


texts:
I'm new to QGIS and am more familiar with MapInfo.

I've been playing around with QGIS to do some habitat mapping based on aerial photography and was wondering if there is a feature in QGIS which automatically draws a poly line or polygon whilst you move the cursor over the screen?

MapInfo has a handy feature called 'auto node' which does this by dropping nodes constantly as the cursor moves. It's useful as you can get very smooth lined polygons which are difficult to achieve if you have to click the mouse everytime to drop a node.

Can QGIS do this???

Thanks!

Steve


texts: How to set the workspace to a featureclass from an unknown file gdb?

texts:
When a user loads an mxd, I am getting the IWorkspace object from the featurelayer's dataset:

IFeatureClass featClass = (IFeatureClass)featLayer.FeatureClass;
IDataset dataSet = (IDataset)featClass.FeatureDataset;
IWorkspace workSpace = (IWorkspace)dataSet.Workspace;

But this obviously does not work if the user has a layer from a file geodatabase (dataset is null).

I can check to see if it's a gdb layer but nothing in the featurelayer or its featureclass tells me the source of the gdb.

How do I go about setting the workspace object to a featurelayer that comes from a gdb?


texts: Percentage of area within a buffer in QGIS

texts: I have a layer containing the land use of a specific area,and a buffer. I would like to calculate the percentage of each category contained in the buffer. For example in the picture i want to know the percentage of the gray parts and orange parts within the buffer. I have tried with intersect but it doesn't work


texts:
I posted a  similar question about this a month or so ago and I was still in early stages of researching the problem thus not to many results came from it .After some time and effort I have come close to having my desired results. The goal is to have a search box that queries the table on my cartodb account and then results get sent back to the API. The code below works as described however that is because I have the value of the form hard coded. I did this to show that the function is correct and that the value does get passed on the the query.My problem thought is once i delete the hard coded value and try to type the text in, the  query results  on my Console  are blank.  So some where the value is getting dropped and I am not quite sure where. (I posted the whole code so that you can copy and paste to see the results...If desired)

<!DOCTYPE html>

<title>Bounds example | CartoDB.js</title>
<meta name="viewport" content="initial-scale=1.0, user-scalable=no" />
<meta http-equiv="content-type" content="text/html; charset=UTF-8"/>
<link rel="shortcut icon" href="http://cartodb.com/assets/favicon.ico" />

   <link rel="stylesheet"
      href="http://libs.cartocdn.com/cartodb.js/v3/themes/css/cartodb.css" />

<form action="#" method="post"  id="form">
    <input type="Text" name="userInput" id="test" value="ARNOLD RD">
    <input type="submit" value="Submit" onclick="processFormData();" />
</form>


texts: OpenLayers 3.5.0, clear highlighted feature after removing feature from source

texts:
I want the user to be able to delete the highlighted feature. This works well, but the feature is not cleared until I move the pointer outside the feature geometry.

I am using the OL 3.5.0 library.

This is the code:

var selectPointerMove_Highlight = new ol.interaction.Select({
  condition: ol.events.condition.pointerMove
  });
map.addInteraction(selectPointerMove_Highlight);

var selectClick_Delete = new ol.interaction.Select({
  condition: ol.events.condition.click
  });
map.addInteraction(selectClick_Delete);

selectClick_Delete.getFeatures().on('change:length', function(e) {
  var feature = e.target.item(0);
  vectorSource.removeFeature(feature);
  selectClick_Delete.getFeatures().clear();
  selectPointerMove_Highlight.getFeatures().clear();
  });

Here is a JSFiddle to show the behaviour. I would expect the highlight to vanish after the clear function is called.

http://jsfiddle.net/goldrydigital/5azhb1o3/

Edit:

I tried a map.renderSync() to force a re-render of the map, but that didn't do it.


texts:
Styling layers in GeoServer can be tedious to say the least. While Styler can relieve some of that tedium, it doesn't yet work in chrome. It can wreck existing styles if not carefully used and still involves a fair amount of manual labor.

I was wondering if there were any flex/AIR SLD generators floating around?


texts:
I have tried several times to run “extract by expression” algorithm in Model Builder with the following expression, but system always crashes, maybe there is something wrong:
CASE
WHEN @maximumnumberofbuildings <= minimum("id_count", group_by:= "id") THEN "id_count" = minimum("id_count", group_by:= "id")
WHEN @maximumnumberofbuildings >= maximum("id_count", group_by:= "id") THEN "id_count" = maximum("id_count", group_by:= "id")
WHEN (@maximumnumberofbuildings > minimum("id_count", group_by:= "id") AND @maximumnumberofbuildings <=median("id_count", group_by:= "id")) THEN "id_count"<= median("id_count", group_by:= "id")
WHEN (@maximumnumberofbuildings < maximum("id_count", group_by:= "id") AND @maximumnumberofbuildings >=median("id_count", group_by:= "id")) THEN "id_count">= median("id_count", group_by:= "id")
END

maximumnumberofbuildings, id_count, id are integer64 type. I have tried also to run it from processing toolbox (with a number instead of maximumnumberofbuildings) but it works very slowly only with one WHEN, so an easier expression, but the output is a empty layer.
Operating system: Windows 10 version 1903
QGIS version: 3.10.6-A Coruña
GDAL version: 3.0.4
GEOS version: 3.8.1-CAPI-1.13.3
PROJ version: Rel. 6.3.2, May 1st, 2020


texts: PostGIS function ST_SymDifference between two segment groups produces unexpected results


texts: Can I statistically smooth attributes spatially (using distance/surrounding polygons) in ArcGIS?

texts:
I have Census Areas with some of the Census attributes. Is it possible to take an attribute, for example, percentage of people in full time employment, and derive a spatially smoothed version of the attribute? If so, how would I go about doing this?

I want to do this to stop getting big shocks from one area to another e.g. where one area has 80% unemployment when all the surrounding areas have 10% unemployment.

By spatially smoothed I'm open to different methods e.g. all polygons surrounding the polygon to be smoothed, or looking at all polygons within a set distance and using some distance-weighting. I will also want to iterate through all of the polygons so that they are all smoothed based on the original attributes.

I'm using ArcMap 10.2 but could use ArcMap 10.4, QGIS, or anything else within reason. Census areas are polygons and are all within one feature class.


texts: Read contents of KML file into CSV format using Python (BeautifulSoup)

texts:
I have a KML file called Layers.KML which has multiple sub folders - Tree_Heights_Elevations, 20190516 - Chain Link Fence, 20190516 - Clearing Limits, 2016-0516 - Solar Array .

Tree_Heights_Elevations subfolders contain multiple points, each of which contains tree heights. Here is the screen shot of the KML file

As you can see in the image above -

Layers.KML has another folder which is Tree_Heights_Elevations which has multiple points. Eg: ASSPOTX which has some information. I want to extract value of Elevation.

I am using BeautifulSoup Python package but unable to extract the Elevation. Could anyone guide me the right way to extract data. I do not know if this is the simplest way to extract data into CSV

import csv, os
from bs4 import BeautifulSoup

inputfile = "Layers.KML"
with open(inputfile, 'r') as f:
  soup = BeautifulSoup(f)

  for node in soup.select('Tree_Heights_Elevations'):
       print(node)


texts: Is there a database/service that allows me to get the terrain type to a coordinate?

texts:
I hope i'm not totally wrong with this kind of question at this stackexchange page. Is there any sort of service avaiable that gives me (or allows me somehow to generate) the areatype to a coordinate? Something like (lat,long) -> position is at a beach. I need this kind of data on a more or less global scale, but "only" at those places, where you have at least a gsm (not satallite) connection - so "only" those areas that are more or less close to populated places.

The closest thing i've found would be to analyse map pictures and store the results, which doesn't seem very elegant to me and probably would take alot of time. Is their any better option avaiable?


texts:
I used to embed images and even websites in feature's form in QGIS 2 using Control Type WebView in field properties.

I don't seem to find the WebView Control Type in QGIS 3 anymore.

So, how can I display an Image in QGIS forms?

(The URL is stored in a field imageLink for each feature).


texts: Can you use the US Census ACS API with TIGER Shapefile GEOID?

texts:
The US Census makes available the TIGER Shapefiles for geographic area like States, Counties, Tracts, etc. (https://www.census.gov/geo/maps-data/data/tiger-line.html). They also supply an API for accessing the survey data. (http://www.census.gov/data/developers/data-sets/acs-survey-5-year-data.html).

I assumed you would be able to simply pass over the GEOID value to this API and get the detail for the the provided geography. Yet is seems the only way is to construct a URL is with the geographic hierarchy like so:

http://api.census.gov/data/2013/acs5?get=NAME,B01001_001E&for=tract:000100&in=state:01+county:073.

It would appear the shapefiles do provide the necessary data to construct this URL but passing GEOID would be so much simpler.

1) Is this possible to use this API providing only GEOID and am I missing it?

2) Am I overlooking something that would not make this possible? For instance, GeoIds not being unique across all geographies.


texts: GeoServer pre-generated tiles not working

texts:
I have a GeoServer running with a layer that I can load in leaflet using WMS.  Currently "caching on demand" seems to work: If I load the map for the second time it loads much faster than loading it for the first time. However I would like to pre seed the map. But it doesn't seem to work: once the map is seeded, it still loads as slow the first time. Then, the second time, as "caching on demand" seems to kick in, it loads faster again. When I inspect the GeoServer data folder, files and folders are created after seeding, so at least tiles are being generated. But after the first request, additional files are created. So to me it looks like preseeding is working, but it does not create all the tiles necessary and so it still needs to create new tiles when I request tiles. Here is a screenshot of the settings that I use in GWC to generate the tiles:

The layer is originally defined in SRS EPSG:28992, but it is reprojected to EPSG:900913 as defined in the settings:

I am using the following code using Leaflet to load the map:

  var mymap = L.map('mapid').setView([52.296904, 4.852940], 15);

var wmsLayer = L.tileLayer.wms('http://localhost:8080/geoserver/wms?',
{
layers: 'acme:pand',
tiled: true

}).addTo(mymap);


texts: Resetting Layer Extent in MXD file using ArcObjects and C#?

texts:
When I change data source properties of a layer using ArcObjects in C#, I lose the extent set for that layer.
How can I reset the layer extent?

Below is the code in which I am resetting database connection properties for a layer causing to loose layer extent values like top, left, bottom, right.

            ESRI.ArcGIS.Carto.ILayer layer = pLayer;

            var dataLayer = (ESRI.ArcGIS.Carto.IDataLayer)layer;
            var datasetName = (IDatasetName)dataLayer.DataSourceName;
            var newWorkspaceName = (IWorkspaceName)((IDataset)targetWorkspace).FullName;

            datasetName.WorkspaceName = newWorkspaceName;
            dataLayer.DataSourceName = (IName)datasetName; ;
            layer = (ESRI.ArcGIS.Carto.ILayer)dataLayer;


texts: Inherited EditorServerCommandAction double posts ReturnMessage

texts:
I have a bunch of inherited EditorServerCommandActions in my WebADF solution and whenever I use the ReturnMessage inside EditorServerAction() it generates two javascript popups at the client side.

Is this a common know error or am I simply missing something?

public class CustomCommandAction : EditorServerCommandAction
{
    protected override void EditorServerAction()
    {
        ReturnMessage("Foo");
    }
}


texts: Duplicate Global IDs within same feature class

texts:
I received a geodatabase from a contractor and I observed there are multiple feature classes with column of Global ID type.

These Global ID type columns, when I checked data i realized there are multiple records with in same feature class having duplicate global ids.

This seems to be against global id column integrity so does this mean this record was generated outside Esri environment?

Can I replace duplicate global ids with new global ids ?


texts: What needs to be done to ArcGIS Server if the server machine/host name is changed?

texts:
We have a Windows Server 2008 where ArcGIS Server 10 was installed.  Not much has been done with the machine yet, but the name of the machine had to be changed after the ArcGIS installation.  The old server name was DAYTONA and the new name of the server is PIGIRON (names changed to protect the innocent).

In the server activity log there is a message:
"Machine daytona is not a valid server container. Machine daytona is not found."

I also see an error in the system event log that reads:
"DCOM was unable to communicate with the computer daytona using any of the configured protocols."

I found a reference to daytona in the file C:\Program Files (x86)\ArcGIS\Server10.0\server\system\Server.dat, but when I remove the machine reference in this file, save it, and restart the SOM, the machine reference gets put back into this file.  This leads me to believe that there are other things that need to be updated because of the server name change.

I found something on ESRI's website about doing this on a Linux or AIX box and the three step solution was uninstall ArcGIS, rename the server, install ArcGIS.  I might end up doing that, but thought there might be a way to do it with out reinstalling everything.

The SOC is not listed in ArcCatalog under this GIS Server.  We tried adding it anyway, and then removing it, but that didn't work.


texts: Finding four states each of six pairs of which have common boundary using PostGIS?

texts:
A Mathematics SE question on common boundary points of connected sets inspires this one.

Construing the word "state" broadly (so that, for example, each province or territory of Canada is a "state") is there any place on earth where, among four states, each a contiguous region, every one of the six pairs among them share a common boundary that is not just a point?

(As a matter of geometry, this can happen with four regions but not with five or more.)

A map of the United States shows Colorado, New Mexico, Arizona, and Utah meeting at an isolated point, and in Canada the provinces of Manitoba and Saskatchewan and the Northwest Territories and the territory of Nunavut similarly share an isolated point. I don't know any others of that sort, but I suspect some exist.

(If the four states are A, B, C, D, then the six pairs are AB, AC, AD, BC, BD, CD.)

The following diagram shows a simplified topology.

The shared boundaries AB, AC, AD, BC, BD and CD should all be lines.


texts: How to make web mapping project run (serve map) faster?

texts:
I'm implementing a project with openlayers 2.12, Geoserver 2.1.3 and PostgreSQL 9.1/PostGIS 2.0

I also use HTML5 Geolocation API and MapQuest's Geocoding API.

In the database I have data (osm_line, osm_point, osm_polygon, osm_road) of a county I downloaded from OpenStreet Map.

The user can also search for points. That's implemented with JavaScript and vector layers.

The problem is , after a certain zoom level (about 16/20) the tiles load very slow. Especially in FireFox and Internet Explorer.

Also, when panning, tiles are blended in Internet Explorer.

Please advice, how can I fix this?

Is there a way to make the whole system load faster?

EDIT

As for the IE problem, with tiles blending, a friend of mine Skyped me this

As for the loading speed issue, I checked the CPU usage and Geoserver goes up to 60-70, when postGis and browsers are about 08-10.

And I also remember that I edited the settings of Geoserver so I increased the SLD file limit.

Also the base layer is WMS, containing the four tables together
{layers: 'cultumap:planet_osm_polygon, cultumap:planet_osm_line, cultumap:planet_osm_roads, cultumap:planet_osm_point'}

Any thoughts?


texts: How do I find content in ArcGIS Online which is not in a folder

texts:
I have a lot of content on AGOL, most of it is divided into folders. However there are some datasets which are not in folders and only exist under All My Content. How do I find the datasets which are not in folders so I can clean up My Content?

For individual items I can see which folder they are stored in by looking at the item's properties, however this will take too long to do manually.

I can't use the filter to exclude folders.

I would rather not pay $1,000 for CleanMyOrg


texts: QGIS and Chrome Browser showing different resolution for same WMTS?

texts:
I have created WMTS using GeoServer and this is tiled service. When i am visualizing this in browser, this is showing me labels and geometry perfectly fine. But when i open this in QGIS using Add WMS/WMTS Tiled Layer, it is blurring/distorting the geometry and labels in QGIS.

How can i fix this problem?

I have checked layer in QGIS2 and QGIS3, there wasn't any difference.

This is WMTS Loaded in Browser

This is QGIS WMTS loaded


texts: Seeking MODIS Land Surface Temperature reference?

texts:
Can anyone point me toward an article that converts the LST values on the MODIS LST Image to degrees Celsius?

My purpose is find a reference (article, book or theses) for this formula that converts DN to degree celsius (C=DN*.02-273.15)


texts: GeoTools: grid transformation using NTv2Transform - error locating grid file

texts:
I'm trying to use GeoTools library to calculate the grid coverage based on NTv2 .gsb file.

However, I'm unable to create the instance of NTv2Transform.

When I use this piece of code:

File gridFile = new File("C:\\Dev\\VTS\\Grids\\WA\\wohpgn.gsb");
NTv2Transform transform = new NTv2Transform(gridFile.toURI());

I get this exception:

Exception in thread "main" org.opengis.referencing.NoSuchIdentifierException: Could not locate NTv2 Grid File file:/C:/Dev/VTS/Grids/WA/wohpgn.gsb
    at org.geotools.referencing.operation.transform.NTv2Transform.<init>(NTv2Transform.java:105)
    at com.frequentis.example.Main.transformFromNad27UsingGrid(Main.java:75)
    at com.frequentis.example.Main.main(Main.java:90)

Am I using the syntax properly?
Is there some other way of how to achieve that?
Do you know about some workaround?

I was trying it with versions 18.4 and 19.1 with Maven.


texts: How to create 3D boxes within a KML file without ground contact?

texts:
I'm actually trying to visualize the output of an atmospheric model within Google Earth using a KML file. The model is divided into several 2.5 x 2.5 x 2.5 km boxes and each box should get a color like green, red, yellow and so on... This already works fine within Matlab like this plot shows:

But now, my question is, how could I visualize this within Google Earth by a KML file? For the polygons I only found the altitude modes of "relative" and "absolut", but I didn't find a solution to create boxes with a height of e.g. 500 m without contact to the ground. So the box should start at e.g. 250 m above ground...

Has anybody an idea how I could realize this? :)


texts: How to divide a polyline by intersecting isochrones

texts:
I have two shapefiles, the first representing a path (a polyline with a single category), the second representing isochrons from the same origin (point 0) which intersect the above mentioned path.
I wonder if it is possible to divide the polyline (path) by the intersecting isochrones (at point 1, point 2, etc.), thus obtaining a segmented polyline with different categories.

QGIS has an useful tool called "line intersections" (vector - analysis tools - line intersections), but it only creates a further shapefile with points at the intersections of the two layers.

Can someone suggest a solution preferably in QGIS or GRASS GIS.


texts: Watershed / drainage network analysis; floating point or integer raster for input?

texts:
I am a little desperate to find a solution to my problem. As I have detailed in this post I am unable to create a detailed drainage network (and catchments) from high resolution LiDAR data.

I think I narrowed down the issue to the use of integer raster vs floating point raster as the input for my analysis.

Both methods yiled dramatically different results, both methods are unasable for a detailed model. Using an integer raster as input yield a very generalized stream network and using a floating point raster yield a detailed, very disconnected and clustered stream network.

Which type of raster should I be using for input and why?

Current method:

convert to integer raster
fill sinks
flow direction
flow accumulation
stream network /stream links
stream order
watersheds (stream link input)


texts: QGIS: WFS LAYER with BBOX query loads all data (URL works fine in Browser/FME workbench)

texts:
I want to pull data from Berlins Geo portal (FIS BROKER).
When I do the query in the browser (firefox) or through FME Workbench, it returns fine:

Eg.:
https://fbinter.stadt-berlin.de/fb/wfs/data/senstadt/s_wfs_alkis_ortsteile?service=WFS&version=2.0.0&request=GetFeature&TYPENAMES=fis:s_wfs_alkis_ortsteile&srsName=urn:ogc:def:crs:EPSG:6.9:25833&BBOX=390812,%205821130.01,%20392562,%205819077.24

When I add a layer in QGIS, it loads everything…

Could somebody help me with this?


texts: Describing weather data in ArcMap area solar radiation analysis

texts:
I am writing up the methodology of the analysis I did on solar radiation availability for a city. I used the tool 'area solar radiation' with ArcMap 10, using the time configuration 'TimeWholeYear(2008)' and latitude 51.

Now I want to write a short paragraph about the weather data used in this analysis (i.e. type of data, data source, etc.) but I cannot find anything about it on the GIS help. Does anyone have any ideas of where to find a description of the weather data used in ArcMap?


texts: I'm a regular user of the 'Eliminate Sliver Polygons' tool in QGIS. I've just switched from 2.18 to 3.4.4 and the 'Eliminate Sliver Polygons' tool isn't available. Has it been removed from this version? Are there any alternatives available?


texts: Adding WMS layer in QGIS standalone Python script

texts:
How do I add a WMS layer on a map, in a standalone Python script with QGIS 2.4?

The methodology explained at the bottom of the official QGIS documentation page seems to be outdated and does not work: http://docs.qgis.org/testing/en/docs/pyqgis_developer_cookbook/loadlayer.html


texts: QGIS Virtual Layers from DB query load as 'missing' when opening project file

texts:
I'm using QGIS 3.16 and I have a layer from a database query that loads as a table, 'Data_vFoodLocations'. I then created a bunch of Virtual Layers with queries like the one below:
select FieldX, FieldY, MakePoint(map_lon, map_lat, 4326)
from Data_vFoodLocations

Whenever I load the project, those Virtual Layers all show up as 'Missing' and broken.
If I open them with 'Edit virtual layer' and re-add them with no changes, they work fine. It's just annoying to do every time I open the project.
I assume(?) it's something along the lines of the DB query loading after the virtual layers? Is there anything I can do to my project file to fix this?


texts: How to make colour tables in QGIS to use in multiple projects?

texts:
I would like to know how to make a colour table that can be saved and used across many projects. I am using this to colour rock types, for example:
Rock one=blue
Rock two= yellow
Rock three= red
etc. I have many rock types and want to ensure that the colouring is the same for all projects. Is there a simple way to apply a universal colour table across many projects?


texts: Zoom into country polygon with Leaflet with mouse click

texts:
I have a GeoJSON file and I'm trying to create a click event to zoom into the specific country. I would also like to zoom out. So far the click events is not working, what am I doing wrong?
    // Creating map object
    var myMap = L.map("map", {
      center: [10.0522, -10.2437],
      zoom: 3
    });
    // Adding tile layer
    L.tileLayer("https://api.mapbox.com/styles/v1/mapbox/light-v10/tiles/256/{z}/{x}/{y}?access_token={accessToken}", {
      attribution: "© Mapbox © OpenStreetMap Improve this map",
      tileSize: 512,
      maxZoom: 18,
      zoomOffset: -1,
      id: "mapbox/light-v10",
      accessToken: API_KEY
    }).addTo(myMap);
    // Load in geojson data
    var geoData = "static/data/map.geojson";
    var geojson;
    // Grab data with d3
    d3.json(geoData, function(data) {
      // Create a new choropleth layer
      geojson = L.choropleth(data, {
        // Define what  property in the features to use
        valueProperty: "coordinates",
        // Set color scale
    scale: ["#ffffb2", "#b10026"],
    // Number of breaks in step range
    steps: 10,
    // q for quartile, e for equidistant, k for k-means
    mode: "q",
    style: {
      // Border color
      color: "#fff",
      weight: 1,
      fillOpacity: 0.8
    },
    // Binding a pop-up to each layer
    onEachFeature: function(feature, layer) {
      layer.bindPopup("Name: " + feature.properties.name + "
Nr Loans:
" + feature.properties.name);
    }
  }).addTo(myMap);
  markers.on("click", function (event) {
    // Assuming the clicked feature is a shape, not a point marker.
    map.fitBounds(event.layer.getBounds());
  });
  markers.addTo(map);


texts:
I have the centerline of a river and multiple lines perpendicular to it on Qgis. I would like to extract the cross-sections in correspondence to those lines, how could I do?
I tried qprof but the .csv in output is not correct (if I try to plot it with excel it doesn't give back the section)


texts: Resampling raster to even coordinates (no decimals) using ArcGIS Desktop?

texts:
I am resampling a raster from one CRS to another, and from 0.5 m resolution to 10 m.

But I want the new cells to have even coordinate extents, instead the old coordinate decimals remain. How do I fix this?

I want the new raster to start at xMin 6875000.000,

but the new raster has 687499.521251

How do I change so new raster has even extents?


texts: calculating share of layer refer to another layer

texts:
I have two vector layers (1) with 107 subdistricts (2) soil map with classification of the soil value. I like to calculate the share of the area of the soil map with values classified as "very_high" for each subdistrict. I would like to calculate it with postgresql code. I tried it with

`select g.gid, g.name as subdistrict, round(ST_area(ST_union(ST_INTERSECTION(g.geom,b.geom)))::numeric/100/g.area_ha, 2) as "share_value_very_high",
ST_union(ST_INTERSECTION(g.geom,b.geom)) as geom_union
FROM soil.soilmap AS b, structure.subdistrict as g
WHERE ST_Intersects(g.geom, b.geom) and b.value = 'very_high'
group by g.gid, g.name, b.value;`

It works well and the calcualtion was fine but the resulting table contains only those subdistricts , all in all 78, having soils of very high value. So I tried a left join to get all subdistricts and join the soil map as following

select g.gid, g.name as subdistrict, round(ST_area(ST_union(ST_INTERSECTION(g.geom,b.geom)))::numeric/100/g.area_ha, 2) as "share_value_very_high",
ST_union(ST_INTERSECTION(g.geom,b.geom)) as geom_union
FROM structure.subdistrict as g
left join soil.soilmap AS b
on ST_within ((b.geom, g.geom)
WHERE b.value = 'very_high'
group by g.gid, g.name, b.value;

But then the calculation was completely wrong e.g. while the correct share would be 90 % the resulting calculation was only 1 %. Also, I have now only 66 subdistricts. Obviously, I have done a big mistake but where is it?.


texts: gdal_translate fails to add raster to existing vector geopackage

texts:
I have an existing GeoPackage soils.gpkg with a single vector feature set.

I am attempting to use Gdal to add a tiff raster to the GeoPackage; however, it fails in the following way:

C:\Users\Me\Documents\Stuff\GIS> gdal_translate -of GPKG tiff_test_wgs84.tiff soils.gpkg -co APPEND_SUBDATASET=YES
0ERROR 1: sqlite3_exec(INSERT INTO gpkg_tile_matrix_set (table_name,srs_id,min_x,min_y,max_x,max_y) VALUES ('tiff_test_wgs84',-1,-77.42268164174656,38.75415988758648,-77.41131759799439,38.76552393133866)) failed: no such table: gpkg_tile_matrix_set
...10...20...30...40...50...60...70...80...90...100 - done.
ERROR 1: failed to prepare SQL INSERT OR REPLACE INTO 'tiff_test_wgs84' (zoom_level, tile_row, tile_column, tile_data) VALUES (5, 0, 0, ?): no such table: main.gpkg_tile_matrix
ERROR 1: failed to prepare SQL INSERT OR REPLACE INTO 'tiff_test_wgs84' (zoom_level, tile_row, tile_column, tile_data) VALUES (5, 0, 1, ?): no such table: main.gpkg_tile_matrix
ERROR 1: failed to prepare SQL INSERT OR REPLACE INTO 'tiff_test_wgs84' (zoom_level, tile_row, tile_column, tile_data) VALUES (5, 0, 2, ?): no such table: main.gpkg_tile_matrix
ERROR 1: failed to prepare SQL INSERT OR REPLACE INTO 'tiff_test_wgs84' (zoom_level, tile_row, tile_column, tile_data) VALUES (5, 0, 3, ?): no such table: main.gpkg_tile_matrix

Does Gdal not try to create the gpkg_tile_matrix and gpkg_tile_matrix_set tables? I know that it does when creating a fresh GeoPackage from the same tiff file.

As far as I know, I've supplied the correct arguments to gdal_translate for appending a raster file to an existing GeoPackage.

I've verified that soils.gpkg is a valid GeoPackage, looked at it in QGis etc. No idea what I'm doing wrong.


texts:
I am working on a RTK application where I use proj7.1 to make necessary transformations. The application is written in python3.7 using subprocess.Popen to call cct and cs2cs which may have performance impact as the software is running on Raspberry PI 3B hardware. I have tried to use the pyproj instead but so for no luck. Maybe someone can give me a pointer how this could be done using my scripts shown below.
I have two scripts first I convert from ITRF2014 epoch 2020.4 to ETRF97(1999.5)
pipeline2 = (
"echo %f %f %f |cct -t 2020.4 -c2,1,3 +ellps=GRS80 +proj=pipeline "
"+step +proj=cart "
"+step +init=NKG:ITRF2014 +inv "
"+step +init=NKG:SE "
"+step +proj=cart +ellps=GRS80 "
)
After that I transform to SWEREF99TM+RH2000 with this command
"echo %f %f %f | cs2cs  +init=epsg:7928 +to +init=epsg:5845 +geoidgrids=SWEN17_RH2000.gtx "
Best regards Assar W.


texts: Sharing ArcGlobe bookmarks with ArcScene and vice versa?

texts:
I am working on a flyover animation and would like to be able to seamlessly go between ArcGlobe and ArcScene and the other way around. I tried saving and using the bookmarks in but this does not seem to work.

Is there another way to synchronize a specific location in ArcGlobe and ArcScene?


texts: Table of statistics of rasters input in QGIS

texts:
I work with QGIS GDAL,I have many raster image 'geotiff' with one band for many dates, I want to create a table that contains statistics for all rasters that I have. I tried the algorithm 'raster layer statistics' in QGIS but it gives just an 'html' file for one raster in input, but what I need is to have as result a table that contains statistics for many raster of many dates. I want result like table (csv,txt format..) that gives:
dates mean min max..

Any help to do this with QGIS?


texts: I clipped a small scene out of a larger orthophoto and the output file is 3 times as big (7.5 GB) as the original tiff-file (2.5 GB). How to get it small while keeping the original resolutuion?


texts: Arcpy fails to refresh page layout after changing the definitionQuery of a layer

texts:
I am trying to automate an export of map images in JPEG format. My script creates the map and exports it, but the export does not honor the changes made to the definitionQuery of the layer in the map (only one layer currently). The code is:

for i in range(int(minuid),int(maxuid)+1):
    layer.definitionQuery='UID > {}'.format(int(i))
    dtFrame.scale=28000
    newExtent = dtFrame.extent
    newExtent.XMin=12345.678
    newExtent.XMax=12345.678
    newExtent.YMin=12345.678
    newExtent.YMax=12345.678
    dtFrame.extent = newExtent
    arcpy.RefreshActiveView()
    mapdoc.save()
    with arcpy.da.SearchCursor(layer,('UID','attr1'),where_clause=None,sql_clause=None) as rq:
        for point1 in rq:
            nameval=point1 [1]
            uidval= point1 [0]
            defq='UID = '+str(int(uidval))
            defq = unicode(defq,"utf-8")
            layer.definitionQuery = defq
            mapdoc.activeView = dtFrame
            mapdoc.activeView=dtFrame.name
            mapdoc.save()
            arcpy.RefreshActiveView()
            pictype="mapexport"
            arcpy.mapping.ExportToJPEG(mapdoc, r"location/"+nameval+str(int(uidval))+" "+pictype)


texts: How does "Show only visible features" work in Map Composer?

texts:
If I tick show only visible, I get no features.

I have created an attribute table in the Map Composer in QGIS trunk
QGIS version 1.9.90-Alpha QGIS code revision 4d6cc6d
I note recent feature update http://hub.qgis.org/issues/4496
but this doesn't explain my experience. I suspect it is misuse on my part.

I would like to select a feature in the main attribute table in QGIS, then have it show only the label for that feature and show only that row in the attribute table in Print Composer.

I explored the options considerably I think, when I uncheck the box I get as many rows as the setting will allow, no relationship to visible. I have checked the Map 1 in the dropdown.

I have written this much, maybe it will help someone else...
I am going to explore reloading the layer, it is in fact a delimited text layer since I am using the read only status of that layer to provide some versioning control. Will see what happens when I reload it as a shapefile.

This is cutting edge cartography as far as my skillset is concerned....


texts: Thematic maps on Non Spatial Data in MapGuide

texts:
I need to create the thematic map on a polygon layer. But the data is not available in the layer; I am using Postgres database for that. It means I need to create the themes on layer using non spatial data.

Please give me any hint to this, it is very urgent.


texts: How to let user upload a raster in google earth engine

texts: I'm working on a Google Earth Engine app where we would like to allow the end user to upload a raster file for use in the calculations. Is there a good way to do this? We already know that a user can manipulate the geometry imports with the drawing tools, but is there a way to just upload a file and then pass the file to an algorithm for processing?


texts: How to use Leaflet to efficiently show parts of a massive dataset?

texts:
I keep seeing non-GIS web developers running into this problem, and I'm not sure what the solution is.

There is some dataset of thousands of items.
We want to show a map to the user, with the visible subset of them shown as interactive, clickable elements.

What methods are there for doing this?

I can think of these, but they're not very satisfactory, so I'm wondering what else there is:

Store all the data in a GeoJSON file, transfer it to the browser, and let Leaflet display it. Problem: doesn't really work with large datasets. TopoJSON raises the limit a bit. It also causes a big delay at page load.
Use Mapbox, store all the data in an interactive layer on Mapbox, and use Mapbox.js to display it. Works great, but costs money, and you can't host it yourself.
Use GeoServer to access a PostGIS database, use the WFS-geojson leaflet plugin to access the data from there. It probably works, but the WFS-geojson Leaflet plugin doesn't seem to be maintained anymore.
Use CartoDB, store all the data in a CartoDB table, and use CartoDB.js to display it. Works great, but can get very expensive. It's possible to host it yourself, but installing CartoDB is non-trivial.

All of this makes me think there must be some much better, free way that I'm missing. What is it?

EDIT

Maybe I wrote off the WFS-geojson plugin too easily. There's a fork that still sees some activity (4 months ago): https://github.com/johanlahti/azgs-leaflet


texts:
I'm looking for good quality, fine resolution 4-band satellite photos of Sweden for academic use (and I'm willing to pay for it). What I'm after is as detailed as possible photos of certain part of Sweden where I would like to classify vegetation with commonly used indices. In order to do so I need NIR band, which is missing in the data provided by Swedish authorities. They instead offer readymade vegetation maps, which are a bit outdated for my purposes. Photos should be taken after 2005, preferably for a couple of years in different seasons.

In addition I'm looking for good quality weather data for Europe and Africa, as accurate as possible in order to evaluate the conditions birds face during migration from wintering grounds.

Could anyone help?


texts: Copy and transfer only the maximum values of data from one column to another in the field calculator

texts:
In the Attribute Table there is a field with the diameter and width of the pipe. I want to transfer all the data from this field to a new field. While in one line there can be several values of the diameter and width as shown in the screenshot, I need only transfer the maximum diameter, in the field calculator.

For example - in the screenshot, move to the new field only the maximum values of 1020 * 10.5 in the new field,  And the values of 720 * 11 should not be transferred to a new field.


texts: Using CARTO VL on local Ubuntu installation

texts:
I'm trying to use CARTO VL on local Ubuntu 16.04 installation.
To connect to local CARTO installation, I used the following config in my index.html:

carto.setDefaultConfig({
  serverURL: 'http://me.localhost.lan:3000'
});

and then simply connect to dataset like this:

const source = new carto.source.Dataset('compact_regions');

I used "me" as the subdomain here and I can use the same dataset from my trial CARTO account but I can't do the same with my local CARTO.
Could anyone please give me some hints or what not? Do I need to do anything with config / https to make this work?


texts: arcpy Problem with updateRow() when trying to update fields with UpdateCursor

texts:
I have a list of values which I am trying to add to an existing shapefile. The field where the values need to be inserted has been added via arcpy.AddField_management and thus the values of the rows are empty or 0. The variable expression1 is a list of floats with the length of the rows in the shapefile. To update the rows I use UpdateCursor:

import arcpy
import math
from arcpy import env
env.workspace = "S:/Black grouse modelling/shapefiles"
inFishnet = "fishnet_label.shp"
rows = arcpy.UpdateCursor(inFishnet)
fields = arcpy.ListFields(inFishnet)
for row, i in zip(rows, expression1):
    for field in fields:
        if field.name == "dist":
            print i
            row.setValue(field.name, i)
            rows.updateRow(row)
del rows
del row

The output I get is as follows:

1265.60109093

Traceback (most recent call last):
  File "<pyshell#218>", line 6, in <module>
    rows.updateRow(row)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\arcobjects\arcobjects.py", line 102, in updateRow
    return convertArcObjectToPythonObject(self._arc_object.UpdateRow(*gp_fixargs(args)))
RuntimeError: ERROR 999999: Error executing function.

The value I want it to update with gets printed as requested (1265.60109093). And I think there isn't a problem with the row.setValue() line. The error comes when I am trying to update the row. I have thought it might be an issue with the float numbers. As these might be longer than the field allows for. However, I have tried rounding the numbers and then inputting those values into the rows but that doesn't work either. I have also tried the following code without setValue() but this doesn't work either:

for row, i in zip(rows, expression1):
    print i
    row.dist = i
    rows.updateRow(row)

Any help is appreciated.

EDIT:

as requested the values of expression1:

expression1[:5]
[1265.6010909259928, 1247.9398102685327, 1232.0557693355686, 1218.0185004559901, 1205.8924977596705]

which is calculated via the following function:

def distFrom(self, other):
        ox = other.x
        oy = other.y
        xDist = self.x - ox
        yDist = self.y - oy
        return (xDist**2 + yDist**2)**0.5

Obviously I do not need these values to be so accurate, but rounding them to 2 decimals for example does not solve the problem.

I have also tried it with a different variable where the values are either True or False and converting these to strings. This does also not solve the problem.

I am using python 2.6.5 and arcGIS 10.0


texts: How do I create a shapefile using the latitude and longitude information in my dBASE table? I am using ArcMap 10.1


texts: What kind of 3D models can QGIS's 3D View handle?

texts:
I'm trying to get QGIS to load a 3D model made in sketchup, as seen here:

However, the model is not loading in the QGIS 3D Viewer.  Any suggestions?  QGIS does not give any input on what type of file they're looking for.


texts: I still can't seem to get this (Add an API key to QGIS for QOSM layer) to work.  I need to add an API key to QGIS for OpenCycleMap. I'm trying to follow the path Web > OpenLayer plugin > OSM/Thunderforest > Set API Key as described in the linked question, but my menu does not have the "OSM/Thunderforest" option. How do I update the OpenLayers Plugin?


texts: GTFS network dataset not working with ArcPy?

texts:
I prepared a GTFS network dataset for Sydney following guidelines of Add GTFS to a Network Dataset tool.

I setup a multiprocessing python script which is working perfectly fine with 64-bit Python in my local computer.

Recently I got access to a powerful VM (Windows 10 Enterprise).

I installed GTFS tool on the VM and registered TransitEvalutor.dll for 64-bit geo processing.

But the GTFS data is giving me an IOError on the VM: "C:\Users\sahab\kitchen\30mincity\GTFS.gdb\SEP18\SEP18_ND" does not exist  ---both for Python 32- and 64-bit. But it is working with Network Analyst tool.

Network Analyst Solve geoprocessing tool only is working when background geoprocessing is disabled.

Can you give a workaround to resolve this issue?


texts: Nearest Neighbor search in PostGIS using cases

texts:
I am using the following PostGIS nearest neighbor query to determine how close a point is to an area (polygon) that has been flooded and it's working well:

SELECT
     ST_AsText(ST_SetSRID(ST_MakePoint(531673.0, 180848.2),27700)),
     "flooding".gid,
      ST_Distance(ST_SetSRID(ST_MakePoint(531673.0, 180848.2),27700),
         "flooding".geom) AS distance
FROM
  "flooding"
WHERE "flooding".geom && ST_Expand(ST_SetSRID(ST_MakePoint(531673.0, 180848.2),27700), 400)
ORDER BY ST_Distance(ST_SetSRID(ST_MakePoint(531673.0, 180848.2),27700), "flooding".geom) ASC
LIMIT 1;

However I would like to categorize the distance as follows and then output the result:

CASE
  WHEN distance > 300 THEN 'OUTSIDE'
  WHEN distance > 50 AND distance  <= 300 THEN 'CLOSE'
  WHEN distance < 50 THEN 'INSIDE'
END as flood_result

I know I need to use the "case when else" format to accomplish the above conditional statement, however I do not know how to do this with regard to ST_Expand.


texts:
I'm trying to read in a shapefile in R, but keep getting this error:

  nul character not allowed (line 1).

I have tried various shapefiles, all which work fine in QGIS and I have even tried open-source ones, all giving the same error.

I have tried various packages, they all give the same error, but here is an example with readOGR:

myShape <- readOGR("C:\", "ne_110m_coastline")

What am I doing wrong please?


texts: Resolving coordinate systems in AutoCAD and ArcMap

texts:
I have been correcting geological maps using Bing imagery in AutoCAD. The imagery is good enough to see the features. I want to transpose to ArcMap to access good topographic information. I'm having trouble resolving at hundred mile scale. Unfortunately, a lot of work was drawn in AutoCAD with an unusual coordinate system (AZ83-WF) that has no ArcMap counterpart. Still, it seems the Bing coordinates should be absolute and splines drawn over it should also be absolute.

I agree that AZ state plane western foot should correspond to AZ83-WF. I have tried that with no luck.The problem seems to be that auto-cad projects the base map, in my case bing maps, and leaves the drawing as is. ArcGis also refuses to reproject autocad objects, and further insists that base maps remain in WGS84. I played with arcgis for autocad today. It has a promising feature called assign and list coordinate system. Unfortunately, it simply ignores my AZ state plane requests.

Ultimately, this seems to be a case of bad luck. I arrived at the AZ state plane trying to resolve a Geological Society of America map (without success). It was just left there when I discovered bing. All you can draw on is WGS84.


texts: Avoid Python conflicts with the system environment (in Windows)

texts:
When are you new to QGIS and PyQGIS what works in the background you might be confused by the different Python Interpreters and customization by settings so called environment variables.
How do avoid such conflicts?


texts: Recent Documents list in ArcMap is not being updated

texts:
On a colleague's machine, the list of Recent Documents isn't being populated when closing an MXD.

she opens map1.mxd
saves the document and closes ArcMap
re-opens ArcMap

map1.mxd isn't shown in the > Existing Maps > Recent section on the welcome screen, or under the > File menu.

Where is the list of recent documents stored? (Interestingly, her welcome screen does list 2 documents, just not the most recently-accessed ones.)

Of course this is all within a corporate system with network drives, locked-down permissions, etc.... ArcMap 10.2.1 on Windows 7.


texts: Buffering with rectangular end caps and not dissolving

texts:
I am using QGIS to attempt to achieve what should be a simple task.

I need to buffer a polygon layer with flat end caps, and without dissolving the resulting layer. In QGIS there are 3 tools for buffering:

v.buffer GRASS command
Fixed Distance buffer in Vector > Geoprocessing
GDAL Buffer vectors tool

v.buffer allows controlling the end cap style and whether the resulting layer is dissolved or not.

Fixed distance buffer in Vector > Geoprocessing outputs an undissolved buffer layer, however the edges are round and overlap. It doesn't seem possible to control the end cap type for the edges.

GDAL buffer vectors tool presents a similar problematic to Fixed distance buffer in Vector>Geoprocessing tools.

How do I achieve a discontinuous (undissolved) buffer with flat edges in the QGIS environment?


texts:
I'm doing a map to describe trekking paths and I need a button for the geolocation function. I saw ol examples about custom control, button and geolocation, and using the drupal geolocate button I managed to get a map. It works, but I need a marker to locate the position, and I added the following code at row 86, but it doesn't work, I'm doing something wrong?

positionFeature.setStyle(new ol.style.Style({
    image: new ol.style.Circle({
        radius: 6,
        fill: new ol.style.Fill({
        color: '#3399CC'
        }),
        stroke: new ol.style.Stroke({
        color: '#fff',
        width: 2
        })
    })
}));


texts: Converting MapInfo to interactive layered PDF to access the full map?

texts:
I want to convert the full map I have in Pathloss 5.0/MapInfo 9.5.1 format into PDF so that I can access the whole map in a PDF file.

If I zoom in, the PDF file will be zoomed into the smallest section of the map (similar to the Pathloss/MapInfo). I can't upgrade my MapInfo version.


texts: Standalone python app that uses qgis.core without installing qgis

texts: I am building a standalone python application that has its own GUI, i would like to use qgis.core/pyqgis as a map display only. But from my searches on this topic, it seems that the full QGis software has to be installed on both developer and prospective users' machine to have this work. Is there not a way to embed Qgis map display without installing the full version of Qgis software?


texts: (ArcMap) ECW to LZW-compressed TIF - why does the filesize pentuple?

texts: I have a few ECW tiles that I want to combine and save as 8 bit unsigned TIF. Together, they're about 180 MB. When I use mosaic to new raster, the resulting TIF ends up being a whopping 1.8GB big. Is ECW simply a far more efficient file format or is there a way to get a comparably sized TIF?


texts: Geographic coordinates; Different latitude & longitude format?

texts:
I am having trouble identifying what format these coordinates are in. Is it DMS?

Latitude    Longitude
N49-46.7    W086-56.4
N51-28.7    W000-27.7
N40-51.0    W074-03.6
N01-21.6    E103-59.4

How do I change it to decimal?


texts: Where do I find the covariance matrix in Erdas Imagine?

texts:
I've just run a principal components analysis on a Landsat image and am trying to calculate my factor loadings. I'm using the textbook Digital Image Processing by John Jenson as a guide and it says that I need to plug my eigenvalues, eigenvectors, and the covariance matrix into an equation to get the factor loading for each component and band. However, while I had the eigenvalues and vectors output to a file, I can't find the covariance matrix anywhere.

Can anyone tell me where I can find the covariance matrix? I also have access to ENVI and ArcMap in case either of those programs can give me it too.


texts: Can I correct position data from a device suffering from the 4/6/19 roll-over problem?

texts:
I work for a company that sells apparatus containing GPS receivers. Over the years a lot of this gear has been sold and deployed. A number of machines in the field contain GPS receivers that feature the week number roll-over issue (WNRO).

Replacing all GPS receivers is costly. Software updates are easier. The receivers however are black boxes to us so if we were to do anything about it with a software update it would have to be on the position data level, the lat-lon data that we read from the device.

I am thinking about using translation tables that map position data that would have been accurate some 19 years ago to position data that is accurate now. I imagine this is possible because the difference in location of satellites between now and 19 years ago would be known. On the other hand, I would think this may be hard because the correction would be different for different regions and  because we do not know what satellites the device is fixed to at any particular time.

Is this possible at all? Do these tables exist? Is there a reasonably accurate way of using data from a rolled-over device to determine ones position?


texts: Graph/Network building and analysis of linked polygons in ArcGIS for Desktop?

texts:
This is more of a follow-up question to Dividing polygons into *n* number of groups of equal counts with ArcPy? and Grouping village points based on distance and population size using ArcGIS Desktop?, particularly the answers by FelixIP.

I am trying to group a set of polygons spatially with a max population constraint. (The Grouping Analysis tool in ArcGIS is similar to desired results, but can only do temporal or spatial constraints, no population -- so some network analysis with ArcPy is needed.)

I am mostly following the steps outlined in FelixIP's answer in the first question, but am stuck at the part where he suggests "[fi] and [ti] are sequential number of connected nodes. To populate this table search this forum on how to assign from and to nodes to link." Unfortunately, there seems to be no such answer (or I'm phrasing poorly), hence my question:

How do I create a network/graph that can be used to calculate the number of sequentially connected nodes?
And in fact, do I correctly understand that these [fi] and [ti] fields represent the total number of nodes that each node can be reached from (fi) and the number of neighbor nodes each node can reach (ti)?

So far I have everything else needed until that point: a Point node layer, and a Polyline link layer (with from and to node). Visual aid below for context:

(source: maks.is)


texts: Get dimension of rectangular polygon PostGIS

texts:
Is there a way to calculate the width and height of a rectangular polygon such as the following?
POLYGON((530934.26 161504.25,
         530941.99 161515.92,
         530945.36 161513.76,
         530937.66 161502.06,
         530934.26 161504.25))

I was trying to go on the root of calculating the diagonal ST_MaxDistance(geom, geom) and area  st_area(geom) and their by get to the dimensions in a mathematical way but doesn't look like it is straight forward.


texts: Rotate Mapbox GL JS map after flyTo transition

texts:
I can't seem to have much luck chaining a rotateTo transition after a flyTo transition, here are 2 examples the first one chains another flyTo transition after the 'moveend' event is fired, this works but if I try to rotate instead, no luck. Am I trying to do something silly?

Fly after transition

Rotate after transition


texts: How to identify gaps in linear point data in FME?

texts:
I have been given a spreadsheet which contains coordinates of welds along a pipeline.  I'm using FME to create points from the coordinates.  On average, the points are roughly 11 metres from each other, with some exceptions.
See below image:

What I would like to do is identify when the gaps are larger than about 15 or 20 metres apart.  For example, the image below shows a gap of about 90m.

As this pipeline is not yet fully built, there will be many such gaps along its length.  Here's an example of a couple of gaps which are several kilometres long (this image is zoomed out so the points look like lines):

One idea I had was to connect the points that are within 15 or 20 metres of each other.  Or perhaps group them together somehow by giving each group a unique ID.  Any change in that ID, will identify a gap in the points.  However, I'm not really sure how I could do either of those in FME.
Any ideas?


texts: WMS scale denominators and zoom levels: why are these tables different?

texts:
I'm trying to understand how WMS scale denominators translate to tile zoom levels.  (Specifically, Google Maps zoom levels.)  I found this answer to What ratio scales do Google Maps zoom levels correspond to? and this OpenStreetMap reference, which are in the same ballpark but differ by about 10%.

Why are they different?

Does it depend on the projection used?

Does it depend on latitude?

The width and height of the WMS request?

Anything else?

Note that I am not interested in true scales here.  Given the min and max scale denominators from the WMS capabilities XML, and using using this approach to calculate the size and position of tiles, I want to know what zoom levels I can expect the WMS endpoint to provide images for.


texts: On the flow generated Symbols in QGIS (2.14), Function

texts:
I'd like to use on the flow generated / assembled symbols in QGIS 2.14, but am not able to code with Python at all  (will later learn it).

Could someone give me a code-snippet with short comments on how to start an external program/perl/script incl catching the the return value (in this case the path to a new created icon).

Guess something alike should be possibel via
 -> Symbols -> SVG -> Data-Defined -> Function

Just to explain my idea..... I started leaning Perl (to handle CSV and JSON) to use it for WeatherUnderground-API calls for my current GIS project, to document and evaluate my recent Experimental Archeological project Living Mesolithic 2.0 (see also e2e.stoneage-network.com).
Due to svg is also only XML, it shouldnt be that difficult to 'assemble' new symbols (icons) which provide more information.


texts: ArcGIS fails to import all records from huge CSV file to file geodatabase table despite being within FGDB table size limits

texts:
I'm using ArcGIS 10.0 on Windows 7 64-bit with 4GB of RAM.

I have some very large tables in CSV format to import to ArcGIS, they all have about 30 fields, upwards of 5 million records per table (a few have double that or more), and file sizes up to about 5 GB. I am trying to import each of them into a file geodatabase as separate tables so I can, ultimately, link them to a feature class and analyze the results in the tables according to their location.

The problem is that ArcGIS seems to just quit importing records at a certain point. I'm using the "Table to Table" tool under Conversion > To Geodatabase, but the "Copy Rows" tool has the same problem. Even if I just add the CSV file directly to ArcGIS without trying to convert it to an FGDB table first, the problem is the same. One of my tables has about 11 million records, and ArcGIS only imports about 10 million of them. ArcGIS doesn't tell me that any error has occurred, the tool just finishes as if nothing is wrong.

I've tried it a few times now and the number of records that make it into the FGDB table is always the same, and doesn't appear to be a file size limit I've ever heard of (not a square of 2 or 16). ArcGIS was able to import another CSV with about 6 million records and all the records came through (though with the problems I'm having with the larger table, the smaller one is kind of suspect now too). ESRI's web site lists the following size limits in a file geodatabase, and I'm far from hitting any of them:

  File geodatabase size: No limit
  Table or feature class size: 1 TB (default), 4 GB or 256 TB with keyword
  Number of feature classes and tables: 2,147,483,647
  Number of fields in a feature class or table: 65,534
  Number of rows in a feature class or table: 2,147,483,647
  Geodatabase name length: Number of characters the operating system allows in a folder
  Feature class or table name length: 160 characters
  Field name length: 64 characters
  Text field width: 2,147,483,647

All I really need to do to these tables are add a couple fields, delete a couple others, and generate values for the new fields (sums of a few of the existing fields). I'm using ArcGIS for it because I'm familiar with the field calculator and I know (or knew, until now) that it could handle tables consisting of millions of records, whereas most other desktop software I have handy (MS Access/Excel) chokes on that many records. So I'm open to using some other piece of software to manipulate the original table and then exporting the (much smaller) resulting table to ArcGIS. Really, the fact that I'm having this problem and that ArcGIS is not giving me any errors or warnings that the problem is even occurring makes me want to handle this data outside ArcGIS as much as possible.


texts: What alternatives are there for hosting shapefiles too large for GeoCommons?

texts: I've got a ~500M shapefile I'd like to place where it is publicly accessible and can be found via search. Specifically this dataset is the Bernalillo County, NM assessor parcel and tax database, which would be of some interest to others, but it way too large for GeoCommons. I'd prefer a geo-data repository site rather than a generic file repository (eg, RapidShare and the like).