texts: Using ArcGIS 10 on a virtual machine with OSX

texts:
I'm currently using ArcGIS 10 on a Windows 7 64 bit Macbook Pro (2.53 Ghz and 8GB of RAM) using VMWare. However, when I am using ArcGIS I find it quite slow even with Windows 7 (64-bit) optimized for performance (no shiny graphics). My work involves making some python scripts and testing them, and whenever I run them, that is also slow, but even opening ArcMap and adding a small shapefile takes a long time (30s to 1 minute).

Is it possible to make my python program in OSX, import the ArcPy module and somehow run them from the Mac side? I realize that ArcGIS runs on Windows, but I was not sure whether any of the ArcPy part can be run separately. Or, could it be quicker to use the 'Unity' feature on Vmware to run these scripts?

I would also appreciate any tips on how much memory/processors to allocate for the virtual machine. At the moment it is 1 processor, and ~4gb of RAM and I would expect it to be quite a bit quicker.


texts: Extract All Pixel Values for each band in a Multi Band Raster

texts: I have a few multi band rasters I'm working with, and I need to perform statistical analysis (in another program) on the data. My mentor told me to determine the value for each pixel in each band, and from there move to the other program. But I'm not sure how to go about this. I've tried a lot of things, including pixel inspector and zonal statistics, but I don't think either of these give me exactly what I'm looking for. Does anyone have any other ideas?


texts: I have a polygon shapefile that corresponds to a point shapefile in the same projection.  I need to move the polygons so that the vertices line up with the points.  Right now they are about 3 feet off as someone shifted them by accident and saved. Is there a way to do this by snapping the vertex from one of the polygons to the point that it is associated with?


texts: Web Appbuilder consuming REST services without Web Map?

texts:
How do you load REST Services into Web AppBuilder programmatically with the Map class?

What file within Web AppBuilder allows developers access to the map class/map object?

Basically, I want to pull in (consume) REST services into ESRI Web AppBuilder without a pre-configured web map.


texts: Using marker symbol to point at the polygon instead of the leader line?

texts:
I am labelling polygons in ArcMap (10.5, standart licence, using Maplex). In case the label does not fit in the polygon I want to use marker symbol (triagle, like an arrow) to point! at the polygon the label refers to - to use this instead of the leader line

for what I want it to look like see green triangles on the right picture (triangle being next to the label and pointing the direction where the "label's" polygon is!!) (vs clasical arrow leader line - see black line with red arrow on the left picture):

I tried  to use leader with line decoration (like on the left pict.) to make it arrow leader and chosen the color of the leader line to white and the smallest line width ("no color" led in leader nor triangle showed up!). That actually helped but the tringle stayed (as expected) at the same position - but it needs to be next to the label!
I thought of finding a way how to set a leader line length and setting that to a really small number so that it was not (or barely) visible - but did not find the way to make this possible to set;
see also Way to set a label leader line length (not using annotation) in ArcGIS Desktop?

I want just one marker sign per label.
The solution should be useful also for situations where all the labels cannot be at the same direction from the polygon.
I cannot use annotations for this unless the whole thing stays dynamic.
I do not want the labels to be too far from the polygon border, say max. offset around 25 pts (to be able to place the arrow too).


texts: Get all overlapping feature information using WMS GetFeatureInfo, OL3, Geoserver

texts:
I've built a simple ol3 map with wms-polygon overlay functionality from a local geoserver 2.8.0. I added a working "WMS GetFeatureInfo" from the ol3 examples.
My problem is that I have some identical overlapping polygons with diffrent attributes. When I click on a polygon (with 10 identical copies) I just get the info from one feature (maybe the top one).

Instead of this I need a results from all of the identical overlapping features.

Does anybody know how to implement this, below is my attempt?

This is my code:

  var wmsSource = new ol.source.TileWMS({
    url: 'http:localhost:8082/geoserver/wms',
    params: {'LAYERS': 'gdi:M4_Porjektkarte'},
    serverType: 'geoserver',
   });

  var wmsLayer = new ol.layer.Tile({
    source: wmsSource
  });

  var view = new ol.View({
    center: [0,0],
    zoom: 3
  });

  var map = new ol.Map({
    layers: [wmsLayer],
    target: 'map',
    view: view
  });

  map.on('singleclick', function(evt) {
    document.getElementById('info').innerHTML = '';
    var viewResolution = /** @type {number} */ (view.getResolution());
    var url = wmsSource.getGetFeatureInfoUrl(
        evt.coordinate, viewResolution, 'EPSG:3857',
        {'INFO_FORMAT': 'text/html'});
    if (url) {
      document.getElementById('info').innerHTML =
          '<iframe seamless src="' + url + '"></iframe>';
    }
  });


texts: What does the 'EE' in 'Quantum GISEE' mean?

texts:
I'm doing some research on QGIS history and found in the German wikipedia article the term 'Quantum GISEE' (https://de.wikipedia.org/wiki/QGIS) but not in the English wikipedia article (https://en.wikipedia.org/wiki/QGIS). I fairly remember 'Quantum GIS' but not 'Quantum GISEE'.

Can somebody clarify this? Or is this just a typo?


texts:
I have two different single band images which I have processed in GEE for two different indices. I need to combine these two using PCA in GEE. And my output values should be scaled between 0-100 to assist in my interpretation. How can I do this?
The script below is what I've used:
// Compute AVI, BI AND SI over L8.
var img = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR');
// Function to cloud mask from the pixel_qa band of Landsat 8 SR data.
function maskL8sr(image) {
  // Bits 3 and 5 are cloud shadow and cloud, respectively.
  var cloudShadowBitMask = 1 << 3;
  var cloudsBitMask = 1 << 5;

  // Get the pixel QA band.
  var qa = image.select('pixel_qa');

  // Both flags should be set to zero, indicating clear conditions.
  var mask = qa.bitwiseAnd(cloudShadowBitMask).eq(0)
      .and(qa.bitwiseAnd(cloudsBitMask).eq(0));

  // Return the masked image, scaled to reflectance, without the QA bands.
  return image.updateMask(mask).divide(10000)
      .select("B[1-9]*")
      .copyProperties(image, ["system:time_start"]);
}
// Map the function over one year of data.
var img = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR')
   .filterDate('2019-01-01', '2019-12-31')
   .map(maskL8sr);
 var img = img.median();
 Map.addLayer(img.clip(table), {bands: ['B5', 'B4', 'B3'], min: 0, max: 0.3});

// Create an AVI image, define visualization parameters and display.
 var avi = img.expression(
        'cbrt(Band5*(1-Band4)*(Band5-Band4))**(1/3)',  {
      'Band5': img.select('B5'),       //450-510nm, NIR
         'Band4': img.select('B4'),    // 640-670nm, RED
         });
             var avi = avi.rename(['avi']);
  Map.centerObject(avi);
//*Map.addLayer(avi.clip(table), {min: -1, max: 1});
Map.addLayer(avi.clip(table), {min: -1, max: 1, palette: ['FF0000', '00FF00']});
Export.image.toDrive({
          image:avi.clip(table),
          region:table,
          description:'AVI',
          scale: 30,
         maxPixels: 1e9
       });

// Create an BI image, define visualization parameters and display.
var bsi = img.expression(
           '((Band6+Band4)-(Band5+Band2)/(Band6+Band4)+(Band5+Band2))*100+100', {

        'Band2': img.select('B2'),       //450-510nm, BLUE
        'Band4': img.select('B4'),    // 640-670nm, RED
        'Band5': img.select('B5'),    // 850-880nm, NIR
        'Band6': img.select('B6')     //15700-16500, SWIR1
           });
 var bsi = bsi.rename(['bsi']);
Map.centerObject(bsi);
Map.addLayer(bsi.clip(table), {min: 0, max: 200, palette: ['FF0000', '00FF00']});
Export.image.toDrive({
          image:bsi.clip(table),
          region:table,
          description:'BSI',
          scale: 30,
         maxPixels: 1e9
       });


texts: Creating Single Field address locator in ArcGIS for Desktop?

texts:
I'm trying to create a single field address locator for places with a feature layer that resolves with the name of the polygons which represents buildings in my map, I want to add alternate names for those buildings.

How can I do it?

I've read about the alternate name table but how does this work? Does the JOINID must be the same from the Primary Table in the locator?

Here's the screenshot of the current configuration:

UPDATE:

I have a map of my university and I want to create an application that tells you the shortest path to the building you want to go:

The issue here is that I created the single field locator with the features to use geocoding and reverse geocoding. The buildings can have more than one name so that's why I'm trying to figure it out how to add more names to my locator.

If you need more information let me know.


texts: Adding Pop-Ups to vector tiles layers (ArcGIS Online)

texts:
I've been trying to get Pop-Ups to work on a published vector tile layer, but can't seem to get them to work.

I followed this help document - https://pro.arcgis.com/en/pro-app/help/data/services/use-vector-tiled-layers.htm - specifically the portion under "Enable pop-ups using an associated web feature layer" but don't get the option to "Set Pop-up Source" when I right click the vector tile layer under contents.

I've published the a feature class layer online to my content, but I can't get to the point of setting the pop-up source to the vector tile layer as I can't even find this option anywhere.

Any Ideas?


texts: pgrouting to calculate distance using a chosen highway system

texts:
I am working on a research to calculate the distance of firms (scattering around the nation) to several metropolises via selected national highway system.

Then I overlay the highway map with firms' location. I have a sample map that is zoomed in to Beijing's map that encompasses both highway and firms:

The purple lines are national highway system captured by the city map (Beijing ) and black dots represent firms' location. My research wants to calculate the distance of firms to another major city (now as a test, I just need to calculate the distance of firms in Beijing to Shanghai.

I follow a post by underdark to calculate the shortest path from a point to another point via a selected path. I am a newbie to pgrouting and the programming. I have some trouble to generate nodes ID and calculate the shortest path.

This is the table of my highway data:

And this is my coding to generate nodes ID:

        CREATE VIEW highway_ext AS
    SELECT *, ST_StartPoint(geom), ST_EndPoint(geom)
    FROM highway

CREATE TABLE node1 AS
   SELECT row_number() OVER (ORDER BY foo.p)::integer AS id,
          foo.p AS geom
   FROM (
      SELECT DISTINCT highway_ext.ST_StartPoint AS p FROM highway_ext
      UNION
      SELECT DISTINCT highway_ext.ST_EndPoint AS p FROM highway_ext
   ) foo
   GROUP BY foo.p;

This is my codes to do a shortpath calc:

   CREATE TABLE network1 AS
   SELECT a.*, b.id as ST_start_id, c.id as ST_end_id
   FROM highway_view AS a
      JOIN node AS b ON a.ST_StartPoint = b.the_geom
      JOIN node AS c ON a.ST_EndPoint = c.the_geom;

      SELECT * FROM pgr_dijkstra('
   SELECT gid AS id,
          ST_start_id::int4 AS source,
          ST_end_id::int4 AS target,
          shape_leng::float8 AS cost
   FROM network',
1,
5110,
false,
false);

As a result, my nodes table is like nothing

Any idea how should I generate node ID and Edge ID? What else do I need to address them?


texts: WGS84 projection is invalid for aggregations

texts:
I am trying to extract the pixel value from the last image of an image collection. The problem is that when I do so, I get an error stating the following:

This is my code for extracting the pixel value of the image:

var precipCol = imageCollection5.select(['hourlyPrecipRateGC'])
    .filterBounds(geometry)
    .filterDate(getFirst(new Date()), currentDate);

var lastPrecip = precipCol.reduce(ee.Reducer.last());

var precipNum = lastPrecip.reduceRegion(ee.Reducer.last(), geom)
    .getNumber('hourlyPrecipRateGC_last').getInfo();

When I click on the inspector tab, the code editor does display the number as it should. The value I get when I click in the inspector is what I need to get from the code above.

The variable geometry is a polygon consisting of the following vertices:

The variable geom is the coordinate: (-106.224911, 31.548092).

The variable currentDate holds a string with the current date in YYYY-MM-DD.

function getFirst(date) returns current week's first date.

imageCollection5 is the following source:

What can I do?


texts:
In GeoServer I configured a layer containing positions that change in time (they are smart phone positions indeed). I'am able to insert fresh position in GeoServer datastore (a PostGIS database), but I'm not able to automatically reload the layer in client side using OpenLayers.

In addition I wish to load from the layer only the information contained in the visualized BBOX (i.e. what the user will see on its screen).

Could anyone provide me a suggestion?


texts: Batch Exporting Survey123 GDB Photos with custom filenames

texts:
I'm trying to export photos from Survey 123 with custom filenames by joining fields in another feature class.

I worked through ESRI's documentation "How To: Batch export attachments from a feature class in ArcMap" and got that working. Then I tried to follow the answers on an identical question asked here on StackExchange to use custom filenames.

After sitting with the top answer for a few hours I decided a simpler route. Instead, I decided to just do a join myself in ArcMap and exported it out as a table. Even still, I haven't been able to get it up and running. I just want to change ESRI's script to use a name based on a field that I joined to the table:

import arcpy
from arcpy import da
import os

inTable = arcpy.GetParameterAsText(0)
fileLocation = arcpy.GetParameterAsText(1)

with da.SearchCursor(inTable, ['OBJECTID_1', 'DATA', 'ATT_NAME', 'field_site_number']) as cursor:
    for item in cursor:
        attachment = item[1]
        filenum = str(item[0])
        filename = str(item[3]) + " " + str(item[2])
        open(fileLocation + os.sep + filename, 'wb').write(attachment.tobytes())
        del item
        del filenum
        del filename
        del attachment

But I keep getting an AttributeError that says:

Traceback (most recent call last):
  File "C:\Users\userDesktop\ExportAttachments.py", line 13, in <module>
    open(fileLocation + os.sep + filename, 'wb').write(attachment.tobytes())
AttributeError: 'NoneType' object has no attribute 'tobytes'

The 'DATA' field in the attachments variable is just populated as Blob. I have no idea why it's prompting this error or how to fix it.


texts: Symbolize Flow Direction on Line Feature in ArcGIS Online

texts:
I seems that ESRI has a serious limitation as it seems you can't symbolize the direction on water lines in ArcGIS Online.  With server you can.

From the reading on Esri forums this is a known problem - caused by complex symbology needed on the line in the publishing MXD.

Just wondering if anybody has found a good work around they would share.


texts: How can ST_VoronoiPolygons be used instead of the broken built-in algorithm in QGIS?

texts:
I tried to use QGIS to create voronoi polygons from a set of points with an attribute (group) and then merge the polygons with the same value for that attribute. (This last step is what QGIS calls Dissolve.) It worked in principle, but since the built-in voronoi algorithm produces erroneous output, I am trying alternatives, especially ST_VoronoiPolygons.

I managed to load a table of points (my_points) into PostgreSQL. It has a geom column (automatically created for PostGIS) and a group column.

I read that I should use the following commands:

drop table if exists voronoi_output;
create table voronoi_output as
    select ST_Dump(ST_VoronoiPolygons(ST_Collect(geom))) from my_points;

But the newly created table voronoi_output is missing the group column, which is needed for the next step (dissolve). How can this be fixed?

(If ST_VoronoiPolygons would accept polylines as input, I could just feed it a bunch of polylines instead of "groups" of points. Then the "dissolve" step would not be needed), but I do not expect that to work.)


texts: Editing ArcGIS Online Hosted Layers in Secured Web Application

texts:
Goal: Develop a web application using ArcGIS Javascript API that has a single login experience for the end user and that enables them to edit feature layers hosted on ArcGIS Online.

Setup: We are using a closed Wordpress instance (registration is off, credentials provided to a small user group) to handle user login/logout. Some layers consumed by the application are hosted on ArcGIS Online, others are hosted on our own instance of ArcGIS Server.

Background: According to this documentation, there are two authentication options for protected layers: (1) user login, and (2) application authentication. We must use the first option because of the limitations associated with application authentication, specifically the inability to edit features.

In order to maintain our single login experience, we are trying to figure out the best way to bypass the ArcGIS Online signin dialog that is associated with user login.

Based on some research, one way to go about doing this is to (a) generate a token using a "generic AGOL user" that has access to the necessary data layers and then (b) use that token with every request of a secured layer.

We can use a post request to successfully generate this token using hard-coded credentials in a PHP file on the server. But, as this post suggests, it is unsafe to pass the token to the client, so it is best to use a proxy to do everything.

Question: We want to use the proxy service that generates a token from hard-coded AGOL credentials and then appends that token to every request to a secure layer made by our Javascript app. Are there any examples of this available for us to reference?


texts: Availability of MXD file for OpenStreetMap data in Esri Geodatabase?

texts:
From what I understand, an ArcGIS mxd file contains the "styling" information on how to show data.

I downloaded and imported OpenStreetMap data into my geodatabase. Now I need to show them.

Is there any available mxd file match the OpenStreetMap fields/values that I can use?

I check this question(How to apply Mapnik style for OSM data in ArcMAP?), but the links are all dead.

ps: right now the all the highway/roads/trail display as a single line in my map. I want to use mxd to show them in a better way, like Mapnik or google.


texts: Need to download Census 2010 population,aging and housing information at block group level

texts: I would like to download 2010 census population,aging and housing information at blockgroup level.Such as Total population,age_5_9...,age_85_up,Household Count,Housing_units,SQML etc


texts: Access control on GeoNetwork catalog services (CSW)

texts: Is it possible to implement access control for the GeoNetwork CSW service? We have two user groups 1) all user and 2) authorized user. We would like to limit some content that can be seen by all user in CSW.


texts:
I am trying to convert a QGIS2 Plugin to QGIS3; specifically the Quickfinder plugin (I am aware of QGIS3 Locator and no I do not feel it provides the same or similar capability just thought I would get that out of the way). I used a Python 27 to Python3 converter to check the code and I changed some import module lines and qt4 to qt5 and resolved some initial errors.

The error I am getting now is

  no module named 'core' C:...\qgis-ltr./python\qgis\utils.py line 335.

Can anyone help me understand this one and where I would start looking for to fix it?


texts:
I'm trying to make a pretty shaded relief map using the technique described here.

I've done this process several times in the past without any problems but today am perplexed about what is going on.  I have a DEM (IMG format) that I've set the color ramp on.  The next step is to export the DEM to an RGB image so that it can be pan-sharpened with the hillshade.  However, when I run the export tool (right click on the layer, Data, then Export
Data...) the resulting layer is not even close to the color scheme I want.  I've made sure that I've checked the "Use Renderer" and "Force RGB" boxes as described in the article.

Here is my example.  The color scheme I want is on the left, what I'm getting from the export is on the right:

I've tried to adjust the symbology on the exported raster to get back to what I wanted but I don't know enough about what is going on to get it set right.

What am I missing here?  How do I get my exported image to match the colors I choose in the DEM?

Using ArcGIS 10.2.1, Windows 7.

-- EDIT --
Here is another screenshot of the whole DEM and the whole exported raster.  This one using a different color ramp (Elevation 1) on the left and the exported raster on the right:


texts: Identify task: Input polygon into IdentifyParameters

texts:
I wish to find all the features within a polygon/circle drawn by my user.

May I know if I should use IdentifyTask

https://developers.arcgis.com/en/javascript/jsapi/identifytask.html

Or is there a more recommended way?

Use Case: User drawn a circle/polygon and it will display all the address within the area.


texts: Anyone using JAMB (Just Another Map Book) on ArcGIS10?

texts: I wanted to experiment with JAMB for creating Plan and Profile figures that our company currently produces in AutoCAD.  I'm having trouble installing JAMB on ArcGIS10, probably because it has not been updated yet, and the jury is still out as to whether it will be.  One of the options I'm looking at is installing VirtualBox and running 9.3.1 side by side with 10.  I contacted ESRI support and they assured me that 9.3.1 can check out 10 licenses, so that looks like a viable option. However it would be easier if JAMB could run on ArcGIS 10.  Any thoughts?


texts: How to create a virtualgpkg in spatialite

texts:
I am currently using OGC GeoPackage features and tiles which is another spatial data format that is also a SQLite3 container.
I am using it in a mobile app but GeoPackage does not support spatial queries and spatial functions like SpatiaLite. I see there is a new feature of SpatiaLite that's a virtualgpkg. Like their virtual Shapefile.
Trying not to reinvent the wheel and take advantage of this cool feature and all the functionality.

In the reference of SpatiaLite 4.2.0 the syntax for creating virtualgpkg  is:

CREATE VIRTUAL TABLE tablename USING module(arg1, ...);

I need more examples of the syntax. How it references a GeoPackage Filename and the GeoPackage Tablename.

Can someone provide the code to create a virtual GeoPackage table in SpatiaLite?


texts: QGIS Temporal Controller - last/most recent object within a temporal frame

texts:
i am working with temporal controller for just a cluple of days - before i've been using time manager. right now i am searching for a way to display the last/latest/most recent feature within a time interval which is defined by the "Frame" settings of the "Temporal Controller" GUI.
finally it should look like img_01 (see below):
the yellow trace is implemented already and works fine.
but the red dot is the thing i am asking for ...
within my time frame of 15min the temporal manager returns a set of 3 (sometimes 4) points. and i would like to implement kind of a selection algorithm that returns just this one point, which is the last one within the frame i.e. the one with the MIN value of an expression that could look like @map_end_time - "points_time"
i have been trying differnet things but

i could not find out how to get a collection/set of geometries within an QGIS expression by using "geometry generator"
and even if i could find out how to get that collection i would not be able to select the one last point with the max-value of the time-field
i've tried to use the temporal controller variables @map_start_time and @map_end_time in QGIS Query Builder to select features on the layer level - but i could not find a way to get along with the syntax for handling date and time SELECTs on geopackage-data within the Query Builder-environment  e.g. how a select of to_datetime(@map_start_time + '2 Hour') could be translated into QGIS Query Builder syntax ? (the expamle works fine in QGIS expressions -> see labels/last line out of 4 in img_02)

img_01:

img_02:


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


texts: Selecting and deselecting features with PyQGIS

texts:
I'm trying to apply the algorithms "Random Points along Lines" at the feature level. I did some PyQGIS code to iterate through the features to select each individual instance and apply the algorithm in the layer with the selection, but it's not working properly.

In fact, the script works fine until I try include the layer.removeSelection() to clear the selection. It triggers an error below (in the run method though). If I don't include that line, the script works fine (but makes selection cumulative.)

I've tried few alternatives, but none seems to work. Does anybody have a suggestion?

Error:

Traceback (most recent call last):
  File "C:\OSGEO4~1\apps\Python37\lib\code.py", line 90, in runcode
    exec(code, self.locals)
  File "<input>", line 1, in <module>
  File "<string>", line 27, in <module>
  File "C:/OSGEO4~1/apps/qgis/./python/plugins\processing\tools\general.py", line 96, in run
    return Processing.runAlgorithm(algOrName, parameters, onFinish, feedback, context)
  File "C:/OSGEO4~1/apps/qgis/./python/plugins\processing\core\Processing.py", line 183, in runAlgorithm
    raise QgsProcessingException(msg)
_core.QgsProcessingException: There were errors executing the algorithm.

Code:

layer = qgis.utils.iface.activeLayer()

temp_layers = []

for f in layer.getFeatures():

    if f['remaining']:
        layer.removeSelection()
        layer.select(f.id())

        #layer.selectByIds([f.id()]) # alternative methods, same error

        print("selected: "+ str(layer.selectedFeatureCount()))
        print("selected: "+ str(layer.selectedFeatureIds()))
        print("new points: "+ str(f['remaining']))

        params = {
            'INPUT' : QgsProcessingFeatureSourceDefinition(layer.id(), True),
            'MIN_DISTANCE' : 5,
            'OUTPUT' : 'memory:',
            'POINTS_NUMBER' : f['remaining']
        }

        results = processing.run('qgis:randompointsalongline', params)
        #results = processing.runAndLoadResults('qgis:randompointsalongline', params)

# merging temp layers
params = {
    'LAYERS': temp_layers,
    'CRS': 'EPSG:3857',
    'OUTPUT': 'memory:'
}

processing.runAndLoadResults('qgis:mergevectorlayers', params)


texts:
When I open a raster file in QGIS (for example a GeoTiff file), I want to quickly check the value of a clicked raster cell. I found there's an 'identify' tool in the toolbar, but when I click with identify I only see a value (Derived) (Clicked coordinate) Band 1

Is there any tool in QGIS to quickly inspect the value of any clicked raster cell?


texts: Can I specify an external WMS service with THREDDS?

texts:
I know that when you use THREDDS, you can specify the internal ncWMS-based WMS service thusly at the top of your catalog:

 <service name="wms" serviceType="WMS" base="/thredds/wms/"/>

The only problem is that the internal WMS service doesn't work for our unstructured (e.g. triangle based grids).

Is there a way to specify an external WMS service here?


texts:
I'm working with a geodataframe of points and I need to create a grid of rectangular polygons.

I know how to have bounds of the geodataframe with gdf.geometry.bounds. Now using this information and with 2 variables length and width  (L, W) of rectangular polygon, how can I create a grid ?


texts: How to convert 2d polylines to 3d polylines in ArcGIS desktop?

texts:
I have some normal 2d polylines for river cross sections. These polylines are 2d and do not contain elevations. I have profile elevation for each cross sections. something like this :
cross section 1:

       station   elevation
         0       963.04
        27.2    963.04
        32.64   963.02
        38.08   962.85
        43.52   962.71
        54.4    962.71
        59.84   962.61
        65.28   962.47
        70.72   962.23

        cross section 2 :
     Station   Elevation
        0       959.43
        5.39    959.43
        10.77   959.41
        16.16   959.22
        21.55   959.1
        48.48   959.1
        53.86   959.15
        59.25   959.28
        64.64   959.43
        134.66  959.43

    cross section 3
    and so forth....

How can I impose these points on my 2d polylines and transfer them to 3d ones ?


texts: How do use gdal.Warp() and gdal.warpOptions() of python API?

texts:
According to http://gdal.org/python/,the gdal.Warp() and gdal.WarpOptions() can be used to warp image from HDF to geotiff fromate by geolocation arrarys, but there is no any example.it would not run, even I try to use .vrt file including lat/lon and database layers sites.

the vrt file :

<?xml version='1.0' encoding='UTF-8'?>
  <VRTDataset rasterXSize="" rasterYSize="">
  <Metadata domain="GEOLOCATION">
    <MDI key="LINE_OFFSET">0</MDI>
    <MDI key="LINE_STEP">1</MDI>
    <MDI key="PIXEL_OFFSET">0</MDI>
    <MDI key="PIXEL_STEP">1</MDI>
    <MDI key="SRS">GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"]]</MDI>
    <MDI key="X_BAND">1</MDI>
    <MDI key="X_DATASET">HDF4_SDS:UNKNOWN:"C:/Users/TLQ/Desktop/11/MYD09.A2011093.0410.006.2015217030905.hdf":38', '[2030x1354] Longitude (32-bit floating-point)</MDI>
    <MDI key="Y_BAND">1</MDI>
    <MDI key="Y_DATASET">HDF4_SDS:UNKNOWN:"C:/Users/TLQ/Desktop/11/MYD09.A2011093.0410.006.2015217030905.hdf":37', '[2030x1354] Latitude (32-bit floating-point)</MDI>
  </Metadata>
  <VRTRasterBand dataType="Float32" band="1">
    <SimpleSource>
      <SourceFilename relativeToVRT="0">HDF4_SDS:UNKNOWN:"C:/Users/TLQ/Desktop/11/MYD09.A2011093.0410.006.2015217030905.hdf":45', '[2030x1354] 1km Band 3 Path Radiance (16-bit integer)</SourceFilename>
      <SourceBand>1</SourceBand>
    </SimpleSource>
  </VRTRasterBand>
</VRTDataset>

the code:

file = "C:/test/MYD09.A2011093.0410.006.2015217030905.hdf"
ds = gdal.Open(file)
subDatasets = ds.GetSubDatasets()
gdal.WarpOptions(geoloc=True, format='GTiff')
gdal.Warp("C:/test/MYD09.A2011093.0410.006.2015217030905.tif",'C:/test/tel.vrt',\           options=gdal.WarpOptions(geoloc=True, format='GTiff'))


texts: Setting icon for the submenu items in QGIS

texts:
I had added a new menu item as well as submenu item using the below code. I had succeeded in setting icon for all the submenus(HeatMapAnalysis,Buffer and Intersect.) But I can't set icon for the main submenu "MySubMenu". I tried setting the icon for the same using the code which is in bold below in the code, but no luck.

def initGui(self):
    """Create the menu entries and toolbar icons inside the QGIS GUI."""

    icon_path = ':/plugins/DSGeneric/icon.png'
    self.add_action(
        icon_path,
        text=self.tr(u'Ds Generic'),
        callback=self.run,
        parent=self.iface.mainWindow())
    self.custom_menu = QMenu(QCoreApplication.translate("mytool", "Mytool"))
    self.iface.mainWindow().menuBar().insertMenu(self.iface.firstRightStandardMenu().menuAction(), self.custom_menu)
    self.custom_submenu = QMenu(QCoreApplication.translate("mysubmenu", "&MySubMenu"))
    **self.custom_submenu.setIcon(QIcon("/icons/logo_icon.png"))**
    self.custom_add_submenu(self.custom_submenu)

    # Heat Map submenu
    icon = QIcon(os.path.dirname(__file__) + "/icons/logo_icon.png")
    self.heat_map_action = QAction(icon, "HeatMapAnalysis", self.iface.mainWindow())
    QObject.connect(self.heat_map_action, SIGNAL("triggered()"), self.heatmapAnalysis)
    self.custom_submenu.addAction(self.heat_map_action)
    # Buffer submenu
    icon = QIcon(os.path.dirname(__file__) + "/icons/logo_icon.png")
    self.buffer_action = QAction(icon, "Buffer", self.iface.mainWindow())
    QObject.connect(self.buffer_action, SIGNAL("triggered()"), self.buffer)
    self.custom_submenu.addAction(self.buffer_action)
    # Intersect submenu
    icon = QIcon(os.path.dirname(__file__) + "/icons/logo_icon.png")
    self.intersect_action = QAction(icon, "Intersect", self.iface.mainWindow())
    QObject.connect(self.intersect_action, SIGNAL("triggered()"), self.intersect)
    self.custom_submenu.addAction(self.intersect_action)


texts:
I'm trying to merge 24 DEM files,

after applying the tool (Mosaic to New Raster) only 11 out of the 24 merged.

All the original files are in the same projection, pixel type, band, etc. I also tried DEM to raster tool first before using Mosaic to New Raster and the results were the same.


texts:
is there a way to enter a numger in the field and the data would fill up according to my specifications?? example : i have to name things that has fixed names only different numbers that are related ..such as the format is "FATx-FDTxx-OLTxxx" and id write only 36 and the data would show "FAT36-FDT9-OLT"

more explaning:

im trying to have a fixed format input in qgis with simple changes .. for example in one of the layers has an attribute field that is practically the same "FAT1-FDT1-OLT1" "FAT6-FDT2-OLT1" and so on

so what im trying to accomplish is having a default input of the "FAT - FDT -OLT" and then i enter the numbers .. i used the text default in the fields properties
but those numbers are related as in the FAT numbers go through 64,the FDT number cycle from 1-4... so now im trying to figure out an easy way to write in the input only the number and the rest would just show up .. example : id write only 36 and the data would show "FAT36-FDT9-OLT"

is there a way to do so?


texts: What to do in CartoDB so that null/zero values are not shown on the map?

texts: I've created a map and want to present some data with the bubbles (supply-demand), but in some locations the supply is zero. I don't want the bubble to be on the map in those cases. How can I do it? Is there any way to set it up in the options?


texts: Converting a list of lists to table and export in CSV in Google Earth Engine

texts:
I want to convert data presented in a list of lists (JSON) in a table so that I can export it in CSV format. The following is the link to my code.

https://code.earthengine.google.com/d44153ff2b315d314f68fcfe1cc90fff

There's a list called 'DataList' which I want to convert to table and export in CSV.

    var collection = ee.ImageCollection(sentinel.filterDate('2016-11-01','2016-11-30')
                         .filterBounds(gulbarga)
                         .sort('CLOUD_COVER_ASSESSMENT')
                         .first());
print(collection.size());
print(collection);

var data = collection.getRegion(gulbarga,30);
print(data,'DataList');


texts:
I am working on Oracle Spatial 11g and I am referring to Oracle docs as a learning resource.

Does any one know any more Oracle Spatial tutorials?

We want to integrate Oracle Spatial data with ArcGIS and QGIS.


texts: Calculating the mean between several raster layers, with NoData values

texts:
I have several Landsat 8 scenes from which i have derived the surface temperature. To obtain an average, I would like to generate the mean between all my scenes. However, prior to this I am planing to delete the cloud cover from each image.

My next step would then be to perform the average of each cell value, ignoring NoData values. That is, if I have 6 scenes in total, and one of them has the value NoData, that cell should calculate the mean using only the other 5 scenes.

How could I get this done? If I am right, the Raster Calculator tool is not an option since it will only perform calculations in those areas where all the layers overlap.


texts: Comparing two Digital Elevation Models (DEMs) from LAS files?

texts:
I have two LiDAR files (.las), one is original let's say with X points. And the other is copy of the first .las file but with Y points, where Y is less than X.

Now, I want to compare the Digital Elevation Models (DEMs) of these two .las files and visualize how different they are.

I want to get information like RMSE, standard deviation, among other types of comparison.

I would appreciate, if anyone could tell me what softwares, and ways to get the comparison info.


texts: I'm plotting various structural geologic SVG symbols and am wondering how to have only a few of the symbols per line - ie., when I use an anticline symbol, I get a repetition of the symbol and when I change the size, I see still many of the same symbol, only larger. What sizing constraint am I missing?


texts: Schema Locks Starting As Soon As User is Disconnected

texts:
I am having continual issues with scheme locks on our ArcSDE database, as soon as users are disconnected schema locks start appearing.  We have 20 or so feature classes and it will show anywhere from 30-90 locks at a time.  The timestamps show that they start up as soon as the user is disconnected, either by a system restart or the DBA disconnecting the user.  I am the only user; however, am not the DBA.

We are running ArcGIS for Server 10.4.1 on SQL Server R2 2012 and I am connecting using ArcGIS Desktop Advanced.  The feature classes in our enterprise geodatabase are published to our Portal as feature services and we do not use the data store.

What could be causing all these schema locks to occur, and is this something on the publisher and software user side, or on the database side?


texts: Debugging ArcObjects SDK 10.3 addins with ArcMap 10.4

texts:
I recently upgraded my version of ArcMap from 10.3 to 10.4.1, including all extensions, patches, etc.  I did not remove the 10.3 ArcObjects SDK, as I wanted to continue to develop 10.3 compatible addins for our organization, as we have not all upgraded yet.  I am using Visual Studio 2013 to build the addins.

I first modified MyAddinProject.csproj.user file to point to my new Desktop10.4\ArcMap.exe.

At my first attempt to debug the addin, after changing no code, ArcMap crashes shortly after splash screen and I get the following error in VS:

A first chance exception of type 'System.Runtime.InteropServices.COMException' occurred in System.Windows.Forms.dll
Additional information: Error HRESULT E_FAIL has been returned from a call to a COM component.

I've tried modifying the Arcmap.exe.config file include the following:

<?xml version="1.0" encoding="utf-8" ?>
<configuration>
  <startup useLegacyV2RuntimeActivationPolicy="true">
    <supportedRuntime version="v4.0" sku=".NETFramework,Version=v4.0"/>
    <!--<supportedRuntime version="v4.0.30319"/>-->
    <!--<supportedRuntime version="v2.0.50727"/>-->
  </startup>

The first and second 'supportedRuntime' lines produce the above error.  However, if I use the v2.0 supportedRuntime, ArcMap runs and the addin works, but I am unable to debug.

Is there something that should be modified to enable me to debug 10.3 addins in 10.4? Or is this not possible?

Related: ArcGIS and .Net compatibility (Extending ArcObjects)


texts: Better approach for selecting related records in multiple featureclasses based on query of related table? (sde/SQL server)

texts:
I have an enterprise geodatabase (SQL) that contains a table ("linktable") which has a key field (Main_ID) that is also in up to a dozen related feature classes (via a relationship class).

I want to run a script in Arctoolbox that allows the user to specify a query statement on the linktable.  This would then gather the Main_IDs of the resulting records and go through the feature classes open in an mxd one by one looking for the related Main_ID field in each of them.  I want the related records in each feature class selected.  The user can pick the selection method (new, add_to, remove, etc).

There are over 180,000 records in linktable.  Potentially tens of thousands in each fc.

I came up with a method 1-2 years ago that involved creating a long SQL statement with IN on each fc.  It works up until that list gets way too long.
It was suggested I repost this issue here and try to solve it using more db-centric ideas.  I'm just not a full-time developer and my python only goes so far.  I would love to find a more efficient solution.  I'm posting the current, inefficient code here as well as a link to my most recent post on this problem

import arcpy, os, string
import arcpy.mapping as MAP

#all relevant fcs as well as Link table should be loaded into current mxd
mxd = MAP.MapDocument("CURRENT")
df = MAP.ListDataFrames(mxd)[0]

#Get Link tables from input
linkTable = arcpy.GetParameterAsText(0)
#Make lists of all feature layers in dataframe
layers = MAP.ListLayers(mxd, "", df)

#get sql expression (sql calculator, query of linktable)
sqlExp = arcpy.GetParameterAsText(1)
arcpy.AddMessage(sqlExp)

#set selection type (new, add_to, etc)
selectType = arcpy.GetParameterAsText(2)

#key field contains Main_IDs that potentially relate Link table to fcs
linkKeyField = "Main_ID"

#Main code to select all records in fc with matching Main_ID in Linktable
if len(layers) > 0:
    # Set the SearchCursor to look through the selection of the linkTable
    sourceIDs = set([row[0] for row in arcpy.da.SearchCursor(linkTable, linkKeyField, sqlExp)])
    # Add DBMS-specific field delimiters
    fieldDelimited = arcpy.AddFieldDelimiters(arcpy.Describe(linkTable).path, linkKeyField)
    # Add single-quotes for string field values
    valueList = ["'%s'" % value for value in sourceIDs]
    # Format WHERE clause in the form of an IN statement
    whereClause = "%s IN(%s)" % (fieldDelimited, ', '.join(map(str, valueList)))
    arcpy.AddMessage("SQL Clause: {}".format(whereClause))
    for lyr in layers:
        if len(arcpy.ListFields(lyr, "Main_ID")) > 0:
            # Process: Select Layer By Attribute
            arcpy.AddMessage("Querying related records in {0}".format(lyr))
            arcpy.SelectLayerByAttribute_management(lyr, selectType, whereClause)

else:
    arcpy.AddMessage("No availble layers for selection.")
    sys.exit(0)

del sqlExp, mxd, df, linkTable, linkKeyField, layers


texts: Vector TileGrid (fishnet) Shapefile to X Y min / X Y max txt file

texts:
I'm trying to achieve a python script to take a TileGrid (fishnet style) of polygons. First I want to add 5 new fields to the shapefile. TileID (Text), YMin (double), YMax (double), XMin (double), and XMax (double). Then the final goal would be to spit a txt file with those fields in the end. Ideally this text file would not include "" and : between the fields. Although, I'm not sure that will be possible.

Now I was able to achieve this (with some success) in ModelBuilder. However, I'd like to put it into an ArcPy code. I had ModelBuilder send out its version of the Python code (attached). However, it is very messy. Would anyone please help me clean this up some? I am very new to Py So I am really insterested in how this should be done.

# Import arcpy module
import arcpy

# Local variables:
Shapefile = ""
Output_Feature_Class__7_ = ""
Output_Feature_Class = ""
Output_Feature_Class__3_ = ""
Output_Feature_Class__4_ = ""
Output_Feature_Class__5_ = ""
Output_Feature_Class__6_ = ""
Output_Feature_Class__2_ = ""
Output_Feature_Class__8_ = ""
Output_Feature_Class__9_ = ""
v2016_lake_cook_county_qquads_shp__2_ = ""

# Process: Add Field
arcpy.AddField_management(Shapefile, "TileID", "TEXT", "", "", "", "",            "NON_NULLABLE", "NON_REQUIRED", "")

# Process: Add Field (2)
arcpy.AddField_management(Output_Feature_Class__3_, "MinX", "DOUBLE", "",      "", "", "", "NON_NULLABLE", "NON_REQUIRED", "")

# Process: Add Field (3)
arcpy.AddField_management(Output_Feature_Class__4_, "MaxY", "DOUBLE", "",     "", "", "", "NON_NULLABLE", "NON_REQUIRED", "")

# Process: Add Field (4)
arcpy.AddField_management(Output_Feature_Class__5_, "MaxX", "DOUBLE", "", "", "", "", "NON_NULLABLE", "NON_REQUIRED", "")

# Process: Add Field (5)
arcpy.AddField_management(Output_Feature_Class__6_, "MinY", "DOUBLE", "", "", "", "", "NON_NULLABLE", "NON_REQUIRED", "")

# Process: Calculate Field
arcpy.CalculateField_management(Output_Feature_Class__7_, "MinX", "!shape.extent.XMin!", "PYTHON_9.3", "")

# Process: Calculate Field (2)
arcpy.CalculateField_management(Output_Feature_Class, "MaxY", "!shape.extent.YMax!", "PYTHON_9.3", "")

# Process: Calculate Field (3)
arcpy.CalculateField_management(Output_Feature_Class__2_, "MaxX", "!shape.extent.XMax!", "PYTHON_9.3", "")

# Process: Calculate Field (4)
arcpy.CalculateField_management(Output_Feature_Class__8_, "MinY",    "!shape.extent.YMin!", "PYTHON_9.3", "")

# Process: Calculate Field (5)
arcpy.CalculateField_management(Output_Feature_Class__9_, "TileID", "  [TILE_NAME]", "VB", "")

# Process: Add Attribute Index
arcpy.AddIndex_management(v2016_lake_cook_county_qquads_shp__2_,   "TileID;MinX;MaxY;MaxX;MinY", "OutputTest.txt", "UNIQUE", "NON_ASCENDING")


texts: ArcMap get frequency of different values in text field

texts:
I am using ArcMap 10.5. I have a large dataset of species records and ultimately want to display a breakdown of the data by different taxon group, e.g. 'bird' 'mammal' 'plant' etc., using something like a pie-chart.

I have a field in my table named  Taxon_Group which lists the taxon group for each record. I have tried using graph wizard to make a graph/chart from the data directly from this field, but it will only accept numerical values and graph the actual values of each record, rather than the count of records that fall into a particular class. I also tried to obtain the 'COUNT' statistic using the Summary Statistics tool, but again, it seems this will only accept numerical values.

The only way I have figured out how to do what I want is to use ModelBuilder, and use the Select tool to select each taxon group with a definition query, and then the Get Count tool to count the rows of the selection for each taxon group. But this is time consuming as I have many taxon groups and will probably involve some manual copy and pasting of frequency values.

Does anyone know of an efficient way to achieve this?


texts: How to create a custom Attribute Inspector for mobile app using the ESRI JS API

texts:
I am making a mobile version of a JS web application I developed.  I'm noticing though that the AttributeInspector widget by default uses standard dojo/form controls.  This means that the form within the widget does not work correctly on a mobile browser.

Is there a way to change the AttributeInspector so that it uses the dojox/mobile controls?  Right now, I am simply providing the LayerInfos object to constructor without specifying the FieldInfos object.  I'd like to avoid that if I can because I'd like the app to be able to use different layers with different fields.

Here is how I create the AttributeInspector:

var layerInfos = this.layers.map(function(layer) {
    return {
        featureLayer : layer.layerObject,
        showAttachments : true,
        isEditable : layer.layerObject.isEditable()
    };
});

this.attInspector = new AttributeInspector({
    layerInfos : layerInfos
}, this.inspectorDiv);
this.attInspector.startup();

Here is what it looks like on my Android phone:

As you can see, the fields get cut off.  Also, the dropdown lists do not appear as they should on a mobile device.

Is there anyway to do this without writing a lot of code to account for different field types and domains?


texts: How can I remove a layer from layerswitcher in OpenLayers?

texts: I have several vectors and WMS layers in layer switcher. But I don't need these layers to be displaying in layer switcher. How can I remove these layers from layer switcher?


texts: How to send WPS request using OpenLayers

texts:
I have generated the WPS request below using GeoServer's WPS Request Builder. I wish to append more parameters to this request (geometry selection), and then send it either using GET or POST from within my OpenLayers JavaScript code. I can't seem to get information on how this can be done. How can I append geometry selection parameter to the request, and how can I send the request using OpenLayers? Any assistance will highly appreciated.

<?xml version="1.0" encoding="UTF-8"?>
<wps:Execute version="1.0.0" service="WPS" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.opengis.net/wps/1.0.0" xmlns:wfs="http://www.opengis.net/wfs" xmlns:wps="http://www.opengis.net/wps/1.0.0" xmlns:ows="http://www.opengis.net/ows/1.1" xmlns:gml="http://www.opengis.net/gml" xmlns:ogc="http://www.opengis.net/ogc" xmlns:wcs="http://www.opengis.net/wcs/1.1.1" xmlns:xlink="http://www.w3.org/1999/xlink" xsi:schemaLocation="http://www.opengis.net/wps/1.0.0 http://schemas.opengis.net/wps/1.0.0/wpsAll.xsd">
  <ows:Identifier>gs:Aggregate</ows:Identifier>
  <wps:DataInputs>
    <wps:Input>
      <ows:Identifier>features</ows:Identifier>
      <wps:Reference mimeType="text/xml; subtype=wfs-collection/1.0" xlink:href="http://geoserver/wfs" method="POST">
        <wps:Body>
          <wfs:GetFeature service="WFS" version="1.0.0" outputFormat="GML2">
            <wfs:Query typeName="learning_workspace:pop_layer"/>
          </wfs:GetFeature>
        </wps:Body>
      </wps:Reference>
    </wps:Input>
    <wps:Input>
      <ows:Identifier>aggregationAttribute</ows:Identifier>
      <wps:Data>
        <wps:LiteralData>pop</wps:LiteralData>
      </wps:Data>
    </wps:Input>
    <wps:Input>
      <ows:Identifier>function</ows:Identifier>
      <wps:Data>
        <wps:LiteralData>Sum</wps:LiteralData>
      </wps:Data>
    </wps:Input>
    <wps:Input>
      <ows:Identifier>singlePass</ows:Identifier>
      <wps:Data>
        <wps:LiteralData>Yes</wps:LiteralData>
      </wps:Data>
    </wps:Input>
  </wps:DataInputs>
  <wps:ResponseForm>
    <wps:RawDataOutput mimeType="text/xml">
      <ows:Identifier>result</ows:Identifier>
    </wps:RawDataOutput>
  </wps:ResponseForm>
</wps:Execute>

My environment: GeoServer 2.1.3, OpenLayers 2.11, PostGIS 1.5.


texts: Networkx Graph from two separate shapefiles

texts:
I have two shapefiles that represent a railway network.

The first shapefile (tracks.shp) contains LineStrings representing tracks
The second shapefile (crossings.shp) contains Points representing rail crossings.

My intention is to create a Graph network where the nodes are these crossing points and edge the tracks between these crossings.

Here is a visual of both shapefiles (using QGIS): The green points are the crossings (vertex)

I saw this answer which is similar but in this case, the vertex for the Networkx graph has to come from another shapefile.

This was the code I wrote:

def generate_graph():
    G_nodes = nx.read_shp("crossings.shp")
    positions_nodes = {k: v for k,v in enumerate(G_nodes.nodes())}

    G_edges = nx.read_shp("tracks.shp")
    positions_edges = G_edges.edges()

    X = nx.Graph()
    X.add_nodes_from(positions_nodes.keys())
    X.add_edges_from(positions_edges)

The issue is I do not know how to sync both graphs. For example, If I try to display the network, I get a KeyError due to discrepancies between both data sources.

X = nx.Graph()
X.add_nodes_from(positions_nodes.keys())
X.add_edges_from(positions_edges)
nx.draw_networkx_edges(X,positions_nodes, positions_edges)

Any pointers? Both on this error and my attempt to sync both shapefiles into one NetworkX graph.


texts:
I'm trying to create a geometry column in a table which I created.

I found the following code:

cur.execute('SELECT addgeometrycolumn(%s, %s, %s, %s, %s)' %(table_name, geometry_column_name, EPSGcode, 'POLYGON', 'XYZ'))

but I always get the same error : no such column table_name

It seems that the addgeometrycolumn is not the same function that the one I found during my research. So please can someone help me to add a geometric column to an existing table?


texts:
I have some really basic newbie questions. I have following data with me:

ESRI Shape files (.shp, .shx, .dbf, .prj, .sbn and .sbx)
KML/KMZ files
GPS coordinates (~400k) of some cars (~50) tracked for a week (csv file with lat, longs and timestamps)

End Goal :
1) First, I want to overlay this data on some web-based system (like google maps), so that I can embed it in my website (I have the flexibility of running a server if a solution requires that).
2) On top of this layer, I want to overlay the gps coordinates based route for each of the cars. The end user should be able to select some of the cars and the routes for only those cars should be visible.
3) Finally, I should be able to simulate the routes, something similar to this http://hobbycoders.com/demos/gpssimulation/

Problem : If someone can direct me to the right set of tools/technologies that I should use to implement this. As of now I am just hopping from one technology to another without knowing what I should choose. The shapefiles are for some hypothetical place in the middle of pacific ocean.


texts: Sharing specific location on map using CartoDB?

texts:
I have created a map of B&Bs/inns in the US.  I am trying to figure out how to send someone a link to the exact location on the map when I zoom in on their blue marker identifying their property.

If I click the share button provided by Carto, I am only able to link to the default zoomed out version of the map.

I would like to generate the custom URLs automatically by clicking a button and not have to do it manually.

Here is the map: http://www.innfan.com/view-map/


texts: Transform stereographic projection using GDAL

texts:
How can I read/wrap a JPG image without any embedded projection/coordinate system to an autoreferenced GEOTIFF with Mercator or similar coordinate system using GDAL?

The image is 500x500 px and I know it uses a polar stereographic projection with center at 40ºN,3ºW and that it's spation resolution is 3.2Km x 3.2Km.

I suppose I've to go with gdalwrap with +proj=stere but it's not clear to me how to specify the other parameters.


texts: Command line-tool to convert geojson file to line-delimited geojson

texts: I'm loving the ndjson-cli tools for working with newline-delimited GeoJSON files. But I can't seem to find a tool that converts a standard GeoJSON (whether minified or prettified) ond converts it to newline-delimited. Is there one?


texts: Record link analysis / validation of records (URL checks) don't use web proxy configuration

texts:
If we try to analyze the record links, geonetwork server ignore the webproxy configuration and try to connect on direct way.

The proxy is set in the admin settings and also with JVM parameters.

For example Harvesting is working over proxy, it only fail at URL checks.

We also tried to set a global system proxy on our Red Hat system but nothing worked.

geonetwork-version: 3.10.1


texts: Issue with loading zoom slider and info template styling in ArcGIS Javascript API

texts:
I have been working on a web app that displays points based on a date query, and an infotemplate pops up when you click on a point. When the user clicks on a point, it retrieves the HTML contents of the infotemplate, but it appears as plain text on the map (no dijit window with an 'x' button in the corner). I assumed that I was loading the tundra CSS theme incorrectly, but that is not the case as I have removed the CSS reference from the HTML of other web maps yet they still retain the infotemplate box. I'm assuming it's an option in the map object but I have had no luck.

As you can see, the zoom bar and infotemplate are not being styled correctly. Does anyone know how to resolve this, or know where these styles are derived?


texts: How to update a postgres/postgis table periodically with multiple data?

texts:
I have created a postgis table from point.shp with 5 fields containing 20 records. Now I want to populate the table with another 150 records keeping the previous 20 records so that the total records count reach 170.

I have tried 'QGIS to Postgis shp loader (SPIT - Shape to PG Import Tool)' and also 'Postgis shapefile & DBF loader' but in both cases the tools are asking to overwrite the previous records in the table.


texts: Web map viewer for QGIS Server and PostGIS

texts:
I have succesfully installed and configured QGIS Server. I tried WMS and WFS via QGIS Desktop and it works perfectly.

Now i am looking for a web map viewer to display (only read (WMS viewer i guess)) the layers being served in QGIS Server. I tried qgis2leaf which generates a leaflet html with a qgis project layers. Quite powerful. I tried also a similar plugin that exports the layers to a openlayers html file.

The problem is that what this does is to EXPORT the data in the layers in the QGIS project to JSON files that are read in the generated html file.

In my particular case, i am using a postgis database and layers' information may change as other users edit it. I want to avoid to regenerate the html files every night or whensoever (I can have 20 or more maps) so that the data is exported again and people can see the changes on the web map viewer. I just want to export the project once and the data of the layers stay linked to the postgis tables.

Is there a good free web map viewer that respects this condition?


texts: Assign multiple domains to multiple feature classes within a geodatabase

texts:
I have a geodatabase that has multiple domains (50 +) and I want to assign them to specific feature class fields (there's over 50 feature classes) in an automated way because I'd like to run this script on a nightly basis.

For example, I have a feature class named "Hydrant" and two domains (Hyd_MAKE, Hyd_SIZE) that I want to assign to the "Hydrant" feature class fields "Make" and "Size" respectively. In addition to the "Hydrant" feature class, I have several other feature class fields (around 50 feature classes) that I want to assign other domains to. Is there a Python script that I can write to accomplish this or should I go with another route?

I know I can use AssignDomainToField to assign domains to a feature class field but I want to automate this with a For Loop or some other iterative process; That hasn't been working out too well.

Here is what I have so far (I am open to a completely different direction; this is just what I originally thought of):

import arcpy
aWS = r"GDBPath"
arcpy.env.workspace = aWS
arcpy.env.overwriteOutput = True

# Assign domain to field
fcList = arcpy.ListFeatureClasses()

for fc in fcList:
    print fc

    fields = arcpy.ListFields(fc)
    for field in fields:
        print "Field: " + field.name

        # Hydrant feature
        if fc == "Hydrant":
            fh = fc

            # name of fields
            f_fh_list = ["Make", "Size"]
            f_fh_list.sort()
            for f_fh in f_fh_list:
                print f_fh

                domains = arcpy.da.ListDomains()
                ##Loop through all domains
                for domain in domains:
                    fh_dn = domain.name
                    fh_dom_list = ["Hyd_MAKE","Hyd_SIZE"]
                    fh_dom_list.sort()
                    if fh_dn in fh_dom_list:
                        arcpy.AssignDomainToField_management(fh, f_fh, fh_dn)

                else:
                     print "Next domain"

        else:
            print "Next feature


texts: Creating WMS proxy without capabilities document?

texts:
I've recently been tasked with making a custom written WMS server be usable with newer clients and apps such as QGis.

The original server was written years ago by a consultancy firm in C++ and who are now no longer in business.

My client, the owners of said custom server at the time it was written never got the source code for the product, and it's quite a specialist application (It connects to custom knowledge base systems throughout their systems) so a simple re-write is out the question.

The biggest issue is that at the time it was built to use a very, very early version of OpenLayers (or similar) and deliver the map imagery into a fairly simple browser based user interface.

Beacuse of the custom nature of everything, all the creators ever implemented was the "GetMap" service call, and not the "GetCapabilities" call.

Beacuse of this, apps like QGis and Esri refuse to talk to the server, and many of the tool kits I'm familiar with such as SharpMap and MapWindow won't allow me to use it as a WMS source either.

I've been looking at using MapProxy to potentially act as a go between that will provide a caps document, but get it's tiles via normal WMS calls to the custom server, but I've no experience with MapProxy so I don't know if that would solve my problem.

Right now, I'm going to have a play with "MapServer" and see if I might be able to construct a custom map file to use for pass through, failing that I'm thinking I might have to build a stand alone server that looks like a regular WMS service to the client apps, but behind the scenes provides essentially a made up on the spot caps doc and uses direct HTTP requests to get tiles to pass through.

Can MapProxy be set up to use WMS as a source, but get by without a capabilities document?

I've been reading the MapProxy docs, and there's nothing in the WMS section that addresses this question.

Update

All great answers so far, I tried the MapProxy approach, and while it worked, it was a bit heavy going on the VM it was running out of.

The proxy is being implemented as a virtual instance sat between the real server and the newer clients.

Installing Java apps on the real server as I've found out is a no/no, so anything Java has to be in a VM, which leaves me with a standard Apache2 and already installed MapServer binary on the real server, so my current approach is to actually leverage MapServer-cgi to talk to the custom service on the same box.

I did actually wonder if I could just get mapserver to deliver the imagery that the custom server was delivering, but no such look as the original developers have encoded the graphics in some weird way.

The old server was definitely a case of, if you want anything doing to this in future, your gonna have to call us to do it type of implementation :-(

I'm still testing different approaches however, so there may still be a few more updates to come.

Final Update

So I actually ended up having to go the custom MapFIle with MapServer route, it's a bit slow, but it works and it works well, at least for the time being, until we can build a completley new server to replace it.

Given that I initially asked this question

"Can MapProxy be set up to use WMS as a source, but get by without a capabilities document?"

I'm going to mark @jgrocha's answer as the answer to the question, as he answered the Q directly, and as I've tested it since, and it would have worked for me had I been allowed to use Java, then I would have most likley gone with that reply.


texts: Understanding ArcMap automatic rematching?

texts:
I am geocoding fairly large data using ArcGIS. I had a considerate amount (16%) of locations had tied results after geocoding using a composite address locator using multiple fields. I am curious what criteria ArcMap use when you hit the "Rematch automatically" button.

Is it choosing the one with the highest score? Will the results of matched records, especially the matched with tied candidates, change after I hit the "Rematch automatically"?

Just want to make sure I know how it is processing the automatic rematch before I use it since this process cannot be undone.


texts: Error while performing kriging interpolation in QGIS 3.8.3

texts:
I am using QGIS 3.8.3 -Zanzibar in order to perform kriging interpolation.
I get this error massage:

"algorithm 'Ordinary kriging' starting…
Input parameters:
{ 'BLOCK' : True, 'DBLOCK' : 100, 'FIELD' : 'lb_acre', 'LOG' : True, 'POINTS' : 'file:///C:/Users/%D7%A9%D7%A8%D7%90%D7%9C/Google%20Drive/%D7%9E%D7%97%D7%9C%D7%A7%D7%94%20%D7%90%D7%92%D7%A8%D7%95%D7%A0%D7%95%D7%9E%D7%99%D7%AA/%D7%A4%D7%A8%D7%95%D7%99%D7%A7%D7%98%D7%99%D7%9D/%D7%A9%D7%A8%D7%90%D7%9C/R&D/wo%20pistachio/yield%20map%20401016-3/4010_16_3_1.csv?type=csv&detectTypes=yes&xField=Lon_adj&yField=latitude&crs=EPSG:4326&spatialIndex=no&subsetIndex=no&watchFile=no', 'PREDICTION' : 'C:/AMD/rg.sdat', 'SEARCH_DIRECTION' : 0, 'SEARCH_POINTS_ALL' : 0, 'SEARCH_POINTS_MAX' : 20, 'SEARCH_POINTS_MIN' : 16, 'SEARCH_RADIUS' : 1000, 'SEARCH_RANGE' : 0, 'TARGET_USER_FITS' : 0, 'TARGET_USER_SIZE' : 1, 'TARGET_USER_XMIN TARGET_USER_XMAX TARGET_USER_YMIN TARGET_USER_YMAX' : None, 'TQUALITY' : 0, 'VARIANCE' : 'C:/AMD/retg.sdat', 'VAR_MAXDIST' : -1, 'VAR_MODEL' : 'a + b * x', 'VAR_NCLASSES' : 100, 'VAR_NSKIP' : 1 }

statistics_kriging "Ordinary Kriging" -TARGET_DEFINITION 0 -POINTS "C:/Users/שראל/AppData/Local/Temp/processing_28777ed756654e7089095e07186c04f2/75a66ffbc5de4eada6fad02bf13dd4d6/POINTS.shp" -FIELD "lb_acre" -TQUALITY 0 -LOG true -BLOCK true -DBLOCK 100.0 -VAR_MAXDIST -1.0 -VAR_NCLASSES 100 -VAR_NSKIP 1 -VAR_MODEL "a + b * x" -TARGET_USER_SIZE 1.0 -TARGET_USER_FITS 0 -SEARCH_RANGE 0 -SEARCH_RADIUS 1000.0 -SEARCH_POINTS_ALL 0 -SEARCH_POINTS_MIN 16 -SEARCH_POINTS_MAX 20 -SEARCH_DIRECTION 0 -PREDICTION "C:/AMD/rg.sdat"
statistics_kriging "Ordinary Kriging" -TARGET_DEFINITION 0 -POINTS "C:/Users/שראל/AppData/Local/Temp/processing_28777ed756654e7089095e07186c04f2/75a66ffbc5de4eada6fad02bf13dd4d6/POINTS.shp" -FIELD "lb_acre" -TQUALITY 0 -LOG true -BLOCK true -DBLOCK 100.0 -VAR_MAXDIST -1.0 -VAR_NCLASSES 100 -VAR_NSKIP 1 -VAR_MODEL "a + b * x" -TARGET_USER_SIZE 1.0 -TARGET_USER_FITS 0 -SEARCH_RANGE 0 -SEARCH_RADIUS 1000.0 -SEARCH_POINTS_ALL 0 -SEARCH_POINTS_MIN 16 -SEARCH_POINTS_MAX 20 -SEARCH_DIRECTION 0 -PREDICTION "C:/AMD/rg.sdat" -VARIANCE "C:/AMD/retg.sdat"

Execution completed in 6.57 seconds
Results:
{'PREDICTION': 'C:/AMD/rg.sdat', 'VARIANCE': 'C:/AMD/retg.sdat'}

Loading resulting layers
**The following layers were not correctly generated.

C:/AMD/retg.sdat

C:/AMD/rg.sdat

You can check the 'Log Messages Panel' in QGIS main window to find more information about the execution of the algorithm."

I would be gratefull for solution.


texts:
I am trying to do a spatial join on two tables, the join runs successfully but there's no data in the new table. what am i doing wrong?

CREATE TABLE us_countyparts_2008_with_count_v5 AS(
  SELECT cp.statefp, cp.countyfp, cp.cd110fp, sum(bd.pop00)
  FROM countyparts2008_with_cd AS cp
  JOIN us_block_2000_with_data AS bd
    ON ST_Contains(cp.geom, bd.geom)
  GROUP BY cp.statefp, cp.countyfp, cp.cd110fp
)

I did IsValid on both and they both returned as valid geometries

This is some of the data from countyparts

and this is some data from us_blocks

I'm pretty sure the error is something to do with the extent
this is the extent for the us_block file

"BOX(-2356113.74289801 -783964.164622752,1044933.24146018 845925.228214008)"

and this is the extent for the county file
"BOX(-124.041983 24.3963080000001,-69.3723069999999 48.540437)"

If this is the problem then how do I fix the extent for the us_block file


texts: How can I fix badly specified geographic coordinate system with units in arc-seconds?

texts:
Check out this depth-to-rock soil data.  It is linked from here.  The spatial reference system is a latlon, but the units are arcseconds.  I could not reason out how to override the SRS even after staring at this page in the Proj4 wiki for a while and trying to convert the origin coordinates using cs2cs with a +towgs84 argument.

Here's the gdalinfo output for that file:

Driver: E00GRID/Arc/Info Export E00 GRID
Files: rockdepm.e00
Size is 6936, 2984
Coordinate System is:
GEOGCS["NAD27",
    DATUM["North_American_Datum_1927",
        SPHEROID["Clarke 1866",6378206.4,294.978698213898,
            AUTHORITY["EPSG","7008"]],
        AUTHORITY["EPSG","6267"]],
    PRIMEM["Greenwich",0,
        AUTHORITY["EPSG","8901"]],
    UNIT["degree",0.0174532925199433,
        AUTHORITY["EPSG","9108"]],
    AUTHORITY["EPSG","4267"]]
Origin = (-449100.000000000000000,177840.000000000000000)
Pixel Size = (30.000000000000000,-30.000000000000000)
Corner Coordinates:
Upper Left  ( -449100.000,  177840.000) (Invalid angle,Invalid angle)
Lower Left  ( -449100.000,   88320.000) (Invalid angle,Invalid angle)
Upper Right ( -241020.000,  177840.000) (Invalid angle,Invalid angle)
Lower Right ( -241020.000,   88320.000) (Invalid angle,Invalid angle)
Center      ( -345060.000,  133080.000) (Invalid angle,Invalid angle)
Band 1 Block=6936x1 Type=Int32, ColorInterp=Undefined
  Min=0.000 Max=191.000
  Minimum=0.000, Maximum=191.000, Mean=123.286, StdDev=38.233
  NoData Value=-2147483647

Notice that the pixel size is 30 arcseconds.  If you divide a given set of coordinates by 3600 you get something that makes sense in degrees.

In the course of composing this question I have discovered that the data is available with other SRSs, but in case those are poorly defined as well I thought I would post anyhow.


texts:
I'm trying to extract areas that have the tag name = "Maroc Telecom using Overpass Turbo, to do so I executed this query:

area
 [name="Maroc Telecom"];
out geom;

When I look at the data tab I find an xml file which contains the features of my query yet on the map I can't see these areas to export them as GeoJSON.

Here is the resulting XML from Overpass Turbo query:

<?xml version="1.0" encoding="UTF-8"?>
<osm version="0.6" generator="Overpass API 0.7.55.4 3079d8ea">
<note>The data included in this document is from www.openstreetmap.org. The data is made available under ODbL.</note>
<meta osm_base="2018-08-27T14:38:02Z" areas="2018-08-27T13:42:02Z"/>

  <area id="2502675189">
    <tag k="amenity" v="telephone"/>
    <tag k="building" v="public_building"/>
    <tag k="name" v="Maroc Telecom"/>
    <tag k="shop" v="mobile_phone"/>
  </area>
  <area id="2643523107">
    <tag k="building" v="yes"/>
    <tag k="name" v="Maroc Telecom"/>
  </area>
  <area id="2651330415">
    <tag k="addr:city" v="taza تازة"/>
    <tag k="name" v="Maroc Telecom"/>
    <tag k="office" v="telecommunication"/>
  </area>
  <area id="2681225147">
    <tag k="addr:city" v="Tanger"/>
    <tag k="building" v="yes"/>
    <tag k="name" v="Maroc Telecom"/>
    <tag k="office" v="telecommunication"/>
  </area>
  <area id="2778063776">
    <tag k="building" v="yes"/>
    <tag k="name" v="Maroc Telecom"/>
    <tag k="office" v="telecommunication"/>
  </area>
  <area id="2799103379">
    <tag k="addr:city" v="Ouarzazate"/>
    <tag k="addr:postcode" v="45000"/>
    <tag k="building" v="commercial"/>
    <tag k="name" v="Maroc Telecom"/>
    <tag k="office" v="telecommunication"/>
  </area>
  <area id="2854549079">
    <tag k="building" v="yes"/>
    <tag k="name" v="Maroc Telecom"/>
    <tag k="office" v="telecommunication"/>
  </area>

</osm>


texts:
I'm looking for an open source GIS server or desktop package that supports as many file formats as possible, is fast, runs on Linux or Mac and can be controlled via an API.

Planning to do very basic editing such as intersects so I guess most GIS platforms will be able to do this.

Does anyone have any recommendations as to a suitable package?


texts: Process for resolving error messages in NEXTGIS

texts:
I'm uplifting a map to NEXTGIS and I get the errors:

Response status code is 502
{"exception": "ValidationError", "message": "GDAL library was unable to open the file.", "attr": "source", "serializer": "raster_layer"}

What do these mean?
How do I resolve the problem?
What is the general procedure for resolving any error messages from NEXTGIS uploads?


texts:
I have a shapefile of 20 cities in ArcMap 10.

I need to position a sheep farm so that it is not within 2 miles of any of these cities.

So, I created a buffer, with 2 miles as distance.

I used the polygon to raster tool.

Now, I need to make this raster so that all space within any of these circles is [0] and all space outside the circles is [1].

I’m not sure how to do this part! I’m new to ArcMap and have been using raster calculator, reclassify tools, no success however.


texts: I am drawing hydro features with the freehand tool. Does anyone know how to delete vertices (i.e. back up) in ArcMap 10 while drawing with freehand in the construction tools? I was just curious if there was a way to do this without having to go ahead and complete the polygon and then go back and edit the vertices post completion. Thanks!


texts:
Any help would be appreciated. I am trying to write a python script to determine flow direction on a water pipe network. I have calculated the mid point of each section of pipe and created a point feature class. I now need to determine the flow direction.

Originally used the ESRI cartographic line feature to do this but need the flow direction to be seperated from the polyline.

When I drew the pipe network it was digitized in flow direction, so all i really need to do is access the from/to field of the polyline and populate the an angle field within the point feature class.


texts: C/C++/C# Library for raster to vector conversion

texts:
I am looking for an "open source" C/C++/C# Library able to convert raster to vector data (Polygonize).
I know gdal_polygonize.py is doing the job, but I can't use Python.
Portrace is an other option, but first it is a tool (and not a Library) and second it is not focussed on GeoData.


texts: Optimizing workflow for FUSION's PolyClipData?

texts:
I have multiple large LiDAR datasets (.laz) from different areas, where I need to extract buildings. I have prepared a shapefile containing a buffer around the buildings I care about, and are using the FUSION tool PolyClipData to clip the point clouds.

I am relatively new to working with ALS data and have a couple of questions to hopefully save me some time as I start processing other datasets. My first data set (after having thinned it) is 780 million points, of which about 15 million points are within the buffer zones. It takes my computer more than five hours to use the PolyClipData tool. Is that normal?

I afterwards will create DEM and DTMs using the clipped data set. In the end, all I care about is the height of the buildings and extracting roof shapes.

I have considered a couple of things that might make the tool faster, but unfortunately don't have the time to test them, and I am hoping to get help knowing if and why any of these (or further) solutions might be helpful.

a) Is it more efficient to have square polygons than rounded polygons to cut the point clouds with?
b) Is it better to dissolve the polygons into a multipolygon or is it better to have multiple polygons?
c) Is it faster to creta DEMs and DTM for the entire data set and then cut those using the polygons? I figure that the interpolation will take much longer with more points. But on the other hand, cutting the output raster from the DEM and the DTM should be fast.
d) I am using .laz files and using the laszip.dll tool to let FUSION work directly on those files. Would I in total save time if I thin my dataset (ThinData) directly into .lda (FUSION native format) or .las?


texts: Calculate the distance between two coordinates (WGS84) in ETRS89

texts:
I have an Excel-file with coordinates (WGS84). From these coordinates I create a line (shp) and transform it (ETRS89). After the transformation I calculate the distance from point to point. I do this using a Python-script and ArcGIS.

However, I would like to calculate the distance between the points before I run my script and without an ArcGIS license. I found several answers to that question (i.e. John Cook,Heversine I,Heversine II and Heversine III). But they are all not accurate enough. I would like to get exact the same value as I get using "Calculate Geometry" in ArcGIS in ETRS89) without using ArcGIS. There is only a difference of a few centimeters, but in sum it is a few kilometers. Is there any option to do that? Any special libaries, which allow to transform coordinates and calculate the distance?


texts: How do I create a line with fixed size from a single point in QGIS?

texts:
My question pretty much sums up my problem. Basically, I have a layer containing points:

And I need to convert them into lines with fixed size (25 meters). I don't care if the point will be taken as starting point or middle point of this line, and I don't care about the direction of each line. By the end I would like to have something like this:

Do you guys know how should I proceed?


texts: Identify Route Locations tool to display Z coordinate

texts:
Identify Route Locations tool working with routing which is connected to M coordinate.

Is there a ArcGIS tool which display route location with Z coordinate? It could be similar aproach - identify nearest vertexes and interpolate Z coordinate.


texts:
I have implemented the style windbarbs as described in the Geoserver documentation here:
https://docs.geoserver.geo-solutions.it/edu/en/multidim/accessing_multidim/rtx/wind_barbs.html
However the barbs are pointing in the wrong direction, away from where the wind is blowing from. They should point TOWARDS where the wind is coming from. To find out what might be wrong I changed the parameter direction to a literal of 180 and then the barbs are pointing due south, indicating a southerly wind which means the barbs are indeed rotated to the correct wind direction.
The next step was to add 180 to the calculated value and then the barbs are actually pointing in the correct direction, but this feels like the wrong way to go about it. Is the atan2 function assuming to add 180 degrees?
<ogc:Add>
 <ogc:Literal>180</ogc:Literal>
 <ogc:Function name="toDegrees">
   <ogc:Function name="atan2">
    <ogc:PropertyName>u</ogc:PropertyName>
    <ogc:PropertyName>v</ogc:PropertyName>
   </ogc:Function>
 </ogc:Function>
</ogc:Add>

The layer store source is a GeoTIFF with the u & v gradient bands extracted/converted from NOAA grib2 file using GDAL:
gdal_translate -b 149 -b 150 -tr 0.25 0.25 -r bilinear gfs.t06z.pgrb2.1p00.f000 wind_250mb.tif


texts: Batch Cutting a Raster File using a Moving Extent in QGIS

texts:
I have a LIDAR file with 20 GB size, and I would like to cut it into smaller areas for gravity terrain correction. It is more like a moving window processing, with LIDAR as the windowed data, then the window is moved each time to output cropped LIDAR files. I thought that I could do that using batch processing with this kind of algorithm:

Prepare the LIDAR file, close "no data" gaps using SAGA > Close Gaps.
Go to raster calculator, put the LIDAR file there as inputs, define output resolution.
Define the initial extent and output file naming.
Keeping all similar, we use something like expression builder to incrementally change the extent for the subsequent 2nd cropped LIDAR, 3rd cropped LIDAR, and so on.

The only thing I could not understand is how I can build such expression in the batch processing window?
Adding a row there does not work for me.


texts: Comparing unclassified super-objects with classified sub-objects

texts:
I have two shapefiles exported from ecognition. The first shapefile contains the land cover classification and the second one contains unclassified superobjects of the first shapefile. I want to compute the percentage of each landcover classes within my superobjects. I also want to get the NDVI percentage of the superobjects based on the vegetation class.

I found this link http://resources.arcgis.com/en/help/main/10.1/index.html#//000800000044000000
but I can't seem to find this in ArcMap. I don't know if I can use this to achieve what I want.


texts: Getting shapefile of river from OpenStreetMap?

texts: I do not want to download the whole map of an area. I just want rivers and other waterbodies present in an area in shapefiles so that I can use them in QGIS. I tried geofabrik but it is giving whole map and that also of full country while i need data of a city. Right now I am trying JOSM. will update you when it will works.


texts: Bing map v7 adding infobox to georss parser module

texts:
Using this plugin (source) to load pins from an xml file, pins get added correctly, but now trying to add infoboxes to each respective pin.

In the function that receives the array of pins I'm trying to loop through each item and add infoboxes accordingly.  Any ideas as to what I'm missing here?

      function GeoRSSImportedCallback(items, bounds) {
          map.entities.push(items);

          for (var i = 0; i<items.length; i++) {
            items[i].Title = items[i].title;
            items[i].Description = items[i].description;
            Microsoft.Maps.Events.addHandler(itmes[i], 'click', displayInfoBox);
          }

          map.entities.push(infoboxLayer);
          map.setView({ bounds: bounds });
      }

    function displayInfoBox(e) {
        pinInfoBox.setOptions({title: e.target.Title, description: e.target.Description, visible:true, offset: new Microsoft.Maps.Point(0,25)});
        pinInfoBox.setLocation(e.target.getLocation());
}

I do initialize the following in the GetMap() function

        pinInfoBox = new Microsoft.Maps.Infobox(new Microsoft.Maps.Location(0, 0), { visible: false });
        infoboxLayer.push(pinInfoBox);


texts: What are MapFish Printing Capabilities for WMS and WFS?

texts:
we are evaluating whether MapFish Printing module fit our needs or not.. From a 2010 presentation http://2010.foss4g.org/presentations/3234.pdf, it says MapFish printing has following readers:

Web Map Server (WMS)
Google Maps
Tile Map Server (TMS)
Image
KaMap
KaMapCache
Open Street Map (OSM)
Tile Cache
GeoJSON

I am wondering does it support

ArcGIS Server WMS?
OpenLayers.Layer.ArcGISCache layer?

besides, what about WFS?
It lists GeoJSON, does this mean it supports GeoJSON-Based WFS?
again, ArcGIS Server WFS?


texts: Maintaining ArcSDE for PostgreSQL database?

texts:
I'm getting familiarized with using PostgreSQL for ArcGIS. I've only ever used SQL Server.

I was successful in setting up an ArcSDE database within PostgreSQL 9.2, but what maintenance should I do for PostgreSQL database?


texts: ArcGIS Desktop Supervised Classification problem

texts:
I have a complex classification problem to solve:

I have multiple shapefiles containing the polygons of all the
parcels of different kinds e.g. a shapefile with all polygons with
maize, another with rice etc
I have a time series multispectral imagery of those parcels and indexes e.g.NDVI calculated
I have made a shapefile with all the training polygons of the classes needed to describe the parcels according to their description: 10 polygons maize, 10 rice etc.
I have constructed a signature file with the above training polygons merged accordingly to the classes, containing all the different layers' mean values: Layers-> All the NDVI rasters, R,G,B layers from all times

Now with these components -or with others if suggested so-, I need to classify all the parcels to a class according to the signature.
For the sake of example lets assume I have 10000 parcels, 20 classes/cultivation types and 30 layers. As a result the Signature file has "size" 20x30.

I need a supervised classification approach, like Maximum Likelihood.

What can I do to assign each parcel to a class?


texts: Is there an opensource way of making 'streetview' panoramas?

texts:
I'm interested in making streetview panoramas and am curious if there are any opensource approaches that I should consider for both making and publishing these panoramas.

I've been playing with the Android functionality on my phone to make Android Photospheres (using a camera app SDK - here for the curious). These can be published in a streetview-like panorama via Google. This works, and is nifty but it requires Google to approval the panorama for it to be published online.

Another approach is to generate my own tiled grid of photos and use the google stretview api:
https://developers.google.com/maps/documentation/javascript/streetview#StreetViewOverlays

I'd much rather produce these using some opensource tools, so I'm not completely locked into Google. In this interview with Steve Coast from OpenStreetMap, he refers to a device you can put on your phone which will result in streetview like imagery. Is there a recommended way for this imagery be collected, stored, stitched together and shared?


texts: Creating points from variable without long/lat using MapBasic?

texts:
Can anyone help me to create a points without long/lat in mapbasic if i have a table with variable integer.

I will have popupmenu for this table and everytime i select one it will create a points equally separated (10m or 15m) to each other (vertically separated) and name each point with suffix of how many points was created.. ex. Loc_1_1, Loc_1_2.... up to Loc_1_10 (if i choose the first row).

The View map would look like below,


texts: How to efficiently get the length of linestrings intersecting a huge raster file?

texts:
I have a huge shapefile (roadnetwork of Europe) and a raster file (Europe, 100x100m resolution).
I want to calculate the length of the roadnetwork for each raster cell, resulting in a new raster layer.
The problem is that the number of raster cells is huge and the roadnetwork is very large.
My current approach in R is to build 60x60 sub-raster (using the SpaDES-package, convert each sub-raster to polygons (using stars and sf) and then intersect the polygons with the road network.
The problem is that this takes very long. Parallelizing the prozess does not work because of memory limitations (st_intersects/st_intersections cause doParallel to make copies of the whole road network for each thread).
Is there a more efficient way to perform the task? Is there some software with which this taks can be done efficiently?
Thanks!


texts: Osm2pgsql failed Cannot detect file format. Try using -r

texts:
I use Osm2pgsql 1.0 and this my action :

sudo osm2pgsql --slim -c --cache 29000 --proj 3857 russia-latest.osm.pbf  -H  xxxx.local -d osm -U osm_user -W -P 5432

was failed after some time:

Osm2pgsql failed due to ERROR: Cannot detect file format. Try using -r.

Log in console betwen starting osm2pgsql and error displaing:

Allocating memory for dense node cache
Allocating dense node cache in one big chunk
Allocating memory for sparse node cache
Sharing dense sparse
Node-cache: cache=29000MB, maxblocks=464000*65536, allocation method=11
Mid: pgsql, cache=29000
Setting up table: planet_osm_nodes
Setting up table: planet_osm_ways
Setting up table: planet_osm_rels
Using built-in tag processing pipeline
Using projection SRS 3857 (Spherical Mercator)
Setting up table: planet_osm_point
Setting up table: planet_osm_line
Setting up table: planet_osm_polygon
Setting up table: planet_osm_roads

Reading in file: russia-latest.osm.pbf
Using PBF parser.
Processing: Node(347886k 399.9k/s) Way(31888k 51.02k/s) Relation(871680 759.97/s)  parse time: 2643s

Reading in file: password=xxxxxxxx
node cache: stored: 347886451(100.00%), storage efficiency: 51.32% (dense blocks: 7163, sparse nodes: 309598190), hit rate: 100.00%

I was confused, cause big node(22Gi),ways and polygon (both 7+Gi) have been created...


texts: VendorOption parameters (graphic-margin,...) based on attribute

texts:
I want to add margins around polygon graphic fills based on an attribute but it doesn't seem to work. Maybe I'm doing something wrong or maybe it's not possible at all.
The same goes for other VendorOption parameters.

This works:

<VendorOption name="graphic-margin">10</VendorOption>

This does not work:

<VendorOption name="graphic-margin">
  <ogc:PropertyName>usuaris</ogc:PropertyName>
</VendorOption>

Nor this:

<VendorOption name="graphic-margin">
  <ogc:Div>
    <ogc:PropertyName>usuaris</ogc:PropertyName>
    <ogc:Literal>50</ogc:Literal>
  </ogc:Div>
</VendorOption>


texts: OpenLayers shows empty tiles if I toggle my layer

texts:
I display my custom tiles with OpenLayers and I have set .olImageLoadError{display: none !important;} in my CSS to avoid displaying pink empty tiles (as described here).

In this way empty tiles are not displayed at first, but if I switch off and on my layer, empty tiles are visualized, as depicted in the upper right corner of the following image:

I'm not the only one having this problem (see here), but I did not find a solution so far.

Have you got any idea on how to solve this?

Thank you.


texts: Determining zip code given lat/long coordinates?

texts:
I am trying to reverse geocode crime data from lat and long coordinates to zip code. I have a robust dataset with lat/long coordinates for each data point, as well as the zip code shape files.

How would I export a file that associates a zip code to each data point based on the initial lat/long coordinate?


texts: LAS sample file for electric utility poles

texts:
Does anybody know a souce where I can find a sample LAS file that has electric utility poles?
Like seen here - http://resources.arcgis.com/en/help/main/10.1/index.html#//015w00000040000000


texts: Where to find neighborhood boundaries for Israel?

texts:
Is there a place that has neighborhood boundary data for Israeli cities as polygons?

I have searched the web without results.


texts: How to publish a single (standardized) QGIS symbology which relates to different datasets of same schema?

texts:
Given a specific (locally extended/stored) symbology and styling in QGIS which refers to several different shapefiles or tables (GeoPackage/Spatialite), how can I easily move or share this with another QGIS desktop user?

Currently, the QGIS symbology is managed in the project file (.qgs) plus there are configuration settings like symbology and font paths. So the symbology is "hard wired" to the dataset and to local symbologies and fonts.

What I want is to publish a single "generic/template" symbology (i.e. signature/symbol settings) with a mapping to exchangeable datasets names - all in paths relative to the symbology (but without the dataset). ArcGIS knows a similar thing they call .MDX and "Map Package".

An example would be a topographic map style from an authority ("mytopostyle" with .qgs, symbole/text fonts, etc.) referring to mydata, which is then published and applied to different, separate regions containing all the same data schema (files mydata_region_A.sqlite, mydata_region_B.sqlite , ...).

Any hints?


texts:
I would like to make use of Signal and Slot while using PyQgis along with PyQt4.

Despite spending some time trying out tutorials and some somewhat related post on gis.stackexchange I can't figure out a proper way of using the canvasClicked signal returned from the QgsMapToolEmitPoint Class (API doc).

Here after is where I'm at. I am building my application from the few tutorials man can find over the internet.

    self.clickPoint_action = QAction(
        QIcon(":/ourapp/zoomin_icon"),
        "Click Point",
        self)

    # create toolbar
    self.toolbar.addAction(self.clickPoint_action)

    # connect the tool
    self.clickPoint_action.triggered.connect(self.clickPoint)

    # create the map tool(s)
    self.tool_clickPoint = QgsMapToolEmitPoint(self.map_canvas)

def clickPoint(self):
    self.map_canvas.setMapTool(self.tool_clickPoint)

So, I understand how to set my map tool (and successfuly use other built-in and custom MapTool), but can't figure out where to go from there to get the corresponding QgsPoint emitted after a mouse click.

I have few experience on Python and just dive into Qt and PyGIS but already have some background mostly wth java.

thanks for sharing!


texts: Function to calculate time series wise max and min

texts:
I have image collection from that I want to make a loop function where I want to call all images that have the same temporal resolution but date range varies (like 2000-06-09 to 2019-06-17) through the time series and calculates one max and one min for this date range. how to do this?

//Image listing
var modis = ee.ImageCollection("MODIS/006/MOD13Q1").filterBounds(table).select('NDVI');

var modis_list = [];
for(var i = 6; i <=10; i++) {
  for(var j = 2000; j <=2019; j++){
    modis_list.push(modis.filter(ee.Filter.calendarRange(i, i, 'month'))
                          .filter(ee.Filter.calendarRange(j, j, 'year')));
  }
}
print(modis_list);
var modis_filtered_flattened =
    ee.ImageCollection(ee.FeatureCollection(modis_list).flatten());

 print(modis_filtered_flattened);


texts:
How can I "extract by expression" the N first (10 highest values of one column, for example)
I need it for the graphical modeler.


texts:
I want to merge the cells[green] in ArcMap as shown in picture below.

Is there any way to do it automatically instead of manually?

There are 200 classes and every time I need to merge the cells for each class.


texts: Measuring distance to roads using raster data in ArcMap?

texts:
I am using ArcMap 10.1 on a Windows x64 system.

I am trying to analyse the driving force of urban land use land cover changes. I have 7 raster files represents distance to roads, distance to river, slope, elevation, population density, precipitation, and land use map within my area. I have been using the “combine” in Arcmap 10.1 to merge multiple raster layers into a new single one, and I have been exporting the attribute table to SPSS to do the binary logistic regression.

In binary logistic regression model, the built-up land is a binary dependent variable, while the slope, the elevation etc. are independent variables. See logistic regression result below: the coefficient (B) of the distance to river is 0, which means there is no relationship between the built-up land changes and the distance to river. Also, the distance to roads got the same result. I don't feel right the regression result, because generally the coefficient of distance from the roads should be negative, which means that the probability of built-up land increases near the roads than away from the roads.

I crated these two distance maps from the roads and the river polyline datasets, respectively, using Eculidean distance in arcmap 10.01, and then doing the regression. See the processing and the results below:

Is there something wrong with my method to calculate the distance to river and roads?

All maps are the same projection and cell size.


texts: Information about the setDataDefinedProperty

texts:
In the QGis C++ API v1.6 there is a function called setlowerValue to set a value for a specific Symbol

In my application i use QGis C++ API 2.8.2, i try to understand how the function setDataDefinedProperty works ?
Is it the same ?
In my application I have a list of symbol corresponding to a depth of the sea and a specific color, I want to categorized that

How can i porceed in 2.8.2 ?
Here is my code :

lrenderer = QgsFeatureRendererV2::defaultRenderer( QGis::Polygon );
QList<QgsSymbolV2*>         lSymbols    = lrenderer->symbols();
QList<QgsSymbolV2*>         lSymbolsToRender    = lrenderer->symbols();
QStringList lListValue;
lListValue  << "0" << "3" << "5"<< "10";

for(int i = 0 ; i < lListValue.size() ; i++) {
    QgsSymbolV2* newsymbol = lSymbols.at(0);
    QgsSymbolLayerV2List  lSymbolLayerList = newsymbol->symbolLayers();

    QgsSymbolLayerV2 * lSymbolLayer;
    lSymbolLayer = lSymbolLayerList.at(0);
    lSymbolLayer->setFillColor(QColor(115, 182, 239)); // couleur de remplissage bleu
    lSymbolLayer->setColor( QColor(125, 137, 140) ); // couleur courbe de niveau
    lSymbolLayer->setDataDefinedProperty("Value",lListValue.at(i));
    //newsymbol->setLowerValue( QString::number(i) ); // Version 1.6 old code
    //((QgsSingleSymbolRendererV2*)lrenderer)->insertValue( lListValue.at(i), newsymbol ); // Version 1.6 old code
    newsymbol->insertSymbolLayer(lListValue.at(i).toInt(), lSymbolLayer);
    lSymbolsToRender.append(newsymbol);
}

There is no information in the API DOC onLine


texts: How to import gdb file data into an Excel sheet

texts:
How can I extract the data from a gdb file into an Excel worksheet?
Is this possible at all?


texts: Why does the raster calculator only give out 1 and 0 values

texts:
I have a problem with the raster calculator of ArcMap 10.5.1 because it only gives out a raster with 0 an 1 values. I'd like to calculate the Leaf Area Index out of a NDVI-Raster by using this equation:

LAI = (((1+NDVI)/(1-NDVI))*NDVI)^0.5

I already have a raster with the NDVI values ranging from approx. -0.699 to approx. 0.823.
In the raster calculator I used the following expression:

(((1+"NDVI.img")/(1-"NDVI.img"))*"NDVI.img")^0.5

where "NDVI.img" is the raster file I named above. I looked at the graph of the equation and it should give out decimal values higher than 0, but I only get a raster file full of 1 and no stretched values. Can you tell me what I did wrong?

PS.: I already tried putting float(...) around the equation but that didn't help. Neither did using the float(...) command on the NDVI.img file beforehand.


texts: Scaling tiles for missing zoom levels in Leaflet?

texts:
I have two sets of tiles. One for zoom level 10 and one for zoom level 14. However, I want to be able to see tiles on my map at zoom levels 9 through 15. For the nonexistent sets of tiles I want to use the tiles from either level 10 or 14 and scale them accordingly. E. g. at the zoom level 9, 11 and 12 I want to have the tiles from the level 10 displayed, but scaled by factor 0.5, 2 and 4 respectively and at the zoom level 13 and 15 I want to have the tiles from the level 14 displayed and scaled by factor 0.5 and 2. I can achieve this for the levels 9 and 15 using the "maxNativeZoom" and "minNativeZoom" options, but how to do this with the other zoom levels?

In Open Layers, which I am using now I can solve this by specifying the "server resolutions" and "client resolutions" options. Now I want to try Leaflet - mostly because of better handling of overlays. That means I want to find a solution that works with overlays, like markers.

I've been looking at the setTransform function, which probably scales the tiles in the case of max/minNativeZoom, but my knowledge of Javascript and HTML is limited and I couldn't figure out how to use it.

Here's my scratch:

        <title>Leaflet map</title>
        <link rel="stylesheet" href="leaflet.css" />