texts:
I am working with both EPSG:4326 and EPSG:4230 in OpenLayers 3 and I am using my WMS service which only supports EPSG:4326 maps. What I need to do is show a map in EPSG:4326 and transform it to EPSG:4230 whenever needed and vice-versa.

Now, I get maps from WMS in EPSG:4326 and define the EPSG:4230 projection with proj4js, then create a view with this newly defined projection. Then I set the map's view. So I expect to see the map in EPSG:4230. Am I doing the right thing by only setting the views projection to EPSG:4230?

My other question is: when I do it this way and perform a zoom in, the map view request maps from my WMS server in EPSG:4230. However my server does not support this projection and, so, I get errors server side and it does not work. What should I do in that situation? Is there any request projection type parameter for zooming?

Here is my code:

EPSG:4230 define;

proj4.defs("EPSG:4230", "+proj=longlat +ellps=intl +towgs84=-87,-98,-121,0,0,0,0 +no_defs");

EPSG:4326 to EPSG:4230

$( "#ed50" ).click(function() {

    var extentTransform = ol.proj.transformExtent(extentOriginal, 'EPSG:4326', 'EPSG:4230');
    var projection2 = ol.proj.get('EPSG:4230');
    projection2.setExtent(extentTransform);

    var view2 = new ol.View({
       projection: projection2,
       center: ol.extent.getCenter(extentTransform)
       zoom: 1
    });

    map.setView(view2);

});

EPSG:4230 to EPSG:4326

$( "#wgs84" ).click(function() {

    var projection2 = ol.proj.get('EPSG:4326');
    projection2.setExtent(extentOriginal);

    var view2 = new ol.View({
        projection: projection2,
        center: ol.extent.getCenter(extentOriginal),
        zoom: 1
    });

    map.setView(view2);
});


texts: understanding bounding box example in openlayers

texts:
I was just checking out the angular-openlayers-directive and came across the following example of bounding box , now see the setting for the bounding box , the code is as follows:

var app = angular.module("demoapp", ["openlayers-directive"]);
        app.controller("DemoController", [ "$scope", "olHelpers", function($scope, olHelpers) {
            angular.extend($scope, {
                offset: 0,
                cairo: {
                    lat: 30.0047,
                    lon: 31.2586,
                    zoom: 10,
                    bounds: []
                }
            });

            $scope.$watch("offset", function(offset) {
                $scope.cairo.bounds[0] += parseFloat(offset, 10);
                $scope.cairo.bounds[1] += parseFloat(offset, 10);
                $scope.cairo.bounds[2] -= parseFloat(offset, 10);
                $scope.cairo.bounds[3] -= parseFloat(offset, 10);
            });
        }]);

I am not quite understanding what exactly is the use of adding bounding box in this example as , even if i take of the bounding off the bounding box settings , nothing really changes. Can anybody explain ?


texts:
In the Python GEE API I'd like to create a .tif file of the selected country boundaries. This file will then be the input of a R script.
I can easily select the country but I don't understand how to use the reduceToImage function to transform it into a raster.
country_code = "CG"
country = ee.FeatureCollection('USDOS/LSIB_SIMPLE/2017').filter(ee.Filter.eq('country_co', country_code))

country.getInfo()

#country = country.reduceToImage({
#    'properties': [],
#    'reducer': ee.Reducer.first()
#})

#task_config = {
#    'image': country,
#    'description':'congo',
#    'assetId': ee.data.getAssetRoots()[0]['id'] + '/congo'
#}
#
#task = ee.batch.Export.image.toAsset(**task_config)
#task.start()


texts:
I need help formatting an outline on a layer. The problem is that i have a layer with this attribute table:

And the aspect of this layer is like that:

What I need is to show an outline not around every "Dicofre", but around the entire layer, i.e. around all "SGRU".
Basically I want this but with an outline around all layer:

I hope I made myself clear.
Thanks in advance.


texts: Encoding categorical variable for random forest with sklearn

texts:
I have a random forest model that works pretty well, taking a bunch of vanilla remote sensing raster data as input. I think it could be improved with addition of some information that I currently have stored as categorical variables (for example: geological substrate, landform, etc.). Because my RF model is built using sklearn, I need to encode these categorical variables numerically.

I've been convinced after reading this and this that one-hot encoding is a bad idea with tree-based methods because it can reduce a single important feature into lots and lots of less important features--bloating the model and reducing performance at the same time. Plus it sounds like a nightmare to add 20+ new raster bands to my image stack. But I'm struggling with which of the alternative encoding techniques to use.

Label Encoding: I understand the critique that this introduces the "dog > cat > horse" problem, but it seems like decision trees are designed such that this isn't actually a problem.

Hashing Trick: As far as I can tell, this is functionally equivalent to assigning largish random numbers to each class. It's not really clear to me how this doesn't suffer from the same ordinality problem as label encoding, but the internet seems to think it's great.

My question is this: If I have a moderate number of categorical variables (less than 40), what is the best way to encode them specifically for use in a random forest?

---EDIT---
I made a vocabulary error in the original post. I am not trying to include a large number of categorical variables/features. Instead, I am trying to include a very small number of categorical variables/features (for example: lithology, landforms), each of which has a small-to-moderate number of classes/levels within my ROI (e.g. "sandstone, granite, limestone"). So the "less that 40" statement above applies to the number of classes/levels within each of a very small number of variables/features, rather than to variables per se. So that is what I am trying to figure out how to encode.


texts: Visualize 3D polygon and TIN stored in postgis

texts: How can we visualize 3D polygon and TIN data stored in POSTGIS? My input is 3D mesh data generated by AR application. I have created polygon by adding all the vertices of an object. I can visualize my objects in the unity, but I want to visualize the data in the form I have stored in postgis. (mostly to verify that the I have stored it is correct) Also, while am running query like ST_Intersects I am getting wrong output based on the visualization in unity. I do not have shape files. I have also, tried viewing it using X3D viewer, by generating X3d file using ST_AsX3D from postgis, but didn't work.Please suggest any tool/ option to visualize this data.


texts: How to elevate cross-sections in HEC-RAS 5.0.3?

texts:
Then creating the cross sections in HEC-RAS it appears that all of the cross sections are in the same elevation. I have tried to change each elevation then creating each cross section, but I get elevations errors.

How do you create cross sections that show that there is a height difference between the first point and the second?


texts: Making an access plan with openstreetmap

texts: What tool would be best to create an access plan (with directions written on roads, local zooms, etc.) from OpenStreetMap (preferably using the MapQuest Open map)?


texts: How to find intersected line is from left or right on another line?

texts:
There is a main LineString and its direction is known. I want to find the side of intersected lines to this main line.

Example image:

Black line is main LineString and red one is intersected LineString.

How to find it is on left side?

Any math solutions, recommendations or code examples can be helpful.


texts: Blue Marble rasters and NaturalEarth shapefiles

texts:
I am working on a personal 3D project which combines metadata from shapefiles, and rasters for topography, true color, bathymetry, etc...

I am using the Blue marble images (which I understand are using the WSG Datum WGS84 and the Plate Carrée projection) with shapefiles from natural-earth (which also use the WSG Datum WGS84).

It was my understanding that both are using the same projection (which could be referenced with epsg:4326) and indeed they seem to match to a point.

However when looking at lands near the the poles, some inconsistencies are very apparent (see screenshot of Greenland below)

Looking at the Blue Marble Navigator (blue-marble.de) which offers an interface to navigate through the Blue Marble images, I noticed the exact same discrepancy (see screenshot below, with the mismatched black zones)

Should I re-project the raster/shapefiles with a different projection, or are the Blue Marble images altogether unfit to be correctly matched to epsg:4326 shapefiles?


texts:
I know that ol3 view accepts center and zoom parameters to properly position the view. Is there a way to initialy render map that fits to provided extent [xmin,ymin,xmax,ymax]?

If no how to calculate proper zoom level that would fit to provided extent?

I know that there is fit method defined in view but that would require the map to be rendered and then zoomed.


texts: Python-Scripting in ArcGIS: Point as Input?

texts:
Is there any possibility to add a Point-Feature as an input parameter in a ArcPy-/ArcGIS-Python-script? If I choose "Point" as datatype in the parameter list, I get 2 input-fields for a X- and a Y-coordinate.
What I want is something like a button to get a crossfade cursor for selecting "on-click" a coordinate in the map.
Is there anything like this?


texts: Need to select and delete features in ModelBuilder

texts:
I'm new to ModelBuilder and I'm planning to take some class of it, but in the meantime I have a project with a relative easy task.

The high-level plan of this project is:

Load shapefile into modelbuilder [USER PARAMETER]
Select by attribute and delete selected features
Load another shapefile into modelbuilder [USER PARAMETER]
Make a buffer of 15 meters of this recently loaded shapefile
ERASE first shapefile with the second shapefile (buffered)
Add Field to the output of the ERASE process.
Use Field Calculator to calculate meters on this new field(!shape.length@meters!)
Select all features that are <200 meters and delete them
Make a Centroid of the result that is >200 meters.

My problem is that I can't select and delete features after the Field Calculator step.

Does someone have an idea of what I did wrong?

A snapshot of what I have so far:


texts:
I have an OpenLayers3 map on my local filesystem, and I want a point layer to use SVG icons. This seems to fail silently:

var style = new ol.style.Style({
    image: new ol.style.Icon({
        size: [95, 95],
        src: "styles/transport_aerodrome.svg"
    })
});

Is this because the SVG, since it is a local file, does not have content-type, but OL3 needs SVGs to have content-type: image/svg+xml? Or is there another reason this fails, and a way to fix it?


texts: Adding multiple fields in join attributes by location in graphical modeler

texts:
I want to add multiple fields in "fields to add" using the join attributes by location algorithm in the graphical modeler.

It works if I add only one field, like DISTANCE, but when I want to add others it gives me errors. I tried separating with ,, ; and '', tried with brackets too. Is it possible to add other fields but not all existing fields?

I'm using QGIS 3.6.


texts: Retrieving user information in ArcGIS Server token protected geoprocessing service using ArcPy?

texts:
I have a geoprocessing service on ArcGIS server 10.3 that is protected with tokens.

This GP Service is a python script. Is it possible to retrieve the user information, such username and password, from inside the script using the token passed as parameter?

Is there a way with arcpy to retrieve user information from a valid token?


texts: how can 2 geotiff rasters be aligned and forced to have same resolution in python using gdal?

texts: I have 2 raster files of different resolution (having different pixel size and origin, but both files are georeferenced). I want to align the 2 files so that each pixel of first raster coincides with the corresponding pixel of second raster. In QGIS, this can be done using the align tool, but I wish to write a script in python for this purpose.


texts: Remove all properties from Image in Google Earth Engine

texts:
How does one go about removing properties from an Image object in Google Earth Engine?

I thought maybe I could remove properties one at a time by setting them to null (which apparently you can do for Features, based on Property with null-value disappearing from feature when creating feature-collection in Google Earth Engine?

How do you remove multiple properties from an image?

Example:

var im = ee.Image("LANDSAT/LT4_L1T_TOA/LT40280331982317XXX01");

var imcc = im.get('CLOUD_COVER');
print(imcc);

var im2 = im.set('CLOUD_COVER', null);
var im2cc = im2.get('CLOUD_COVER');
print(im2cc);
print(im2);

Which sets the property CLOUD_COVER to null, but doesn't actually remove the property itself:

I'm only trying to do this as an exercise to learn the structures of objects in Google Earth Engine and how to manipulate them.


texts: Placing a table in Layout View in ArcMap 10.5

texts:
I'm trying to place a table in the Layout View for a map that I am making.

However, when I click on Add Table to Layout from the table options, it only adds the column names to the layout... How can I resolve this issue?

Image showing the issue I'm talking about.


texts: QGIS - Missing GDAL Tools plug-in (merging Raster TIFFs of satellite image)

texts: I am trying to install the GDAL Tools plug-in but it is missing from the Plug-ins list. I am looking to merge multiple Raster TIFFs which make up my satellite image basemap, exported using Atlas, which I can use offline.


texts: Is it possible to get coordinates for point from satellite image?

texts: Is it possible to get coordination for point from satellite image using any programming language?  For example, if I clicked on a specific point on a satellite image hosted on a web application, I would like to retrieve Lon/Lat for that point.


texts: How can I enable Shapefile GUI Loader in pgadmin3?

texts:
I have installed PostgreSQL 9.1 with PostGIS 1.5.3 on Ubuntu 11.10.

I also installed pgadmin3 1.14.0-beta1, following this instructions.

If I remember correctly, when I installed PostGIS on Windows pgadmin automagically came with Shapefile GUI Loader plugin already installed. On Ubuntu however, the plugin drop down menu is empty.

Instructions I found seems to be Windows specific.

Is there a way to install / enable it on Ubuntu?


texts: How to fit bounds after adding multiple markers

texts:
I am using this piece of code to display map markers pulled from a MySQL database using leaflet js on an open street map, but the map bounds are obviously hard coded at the top (with setView), and I have no idea how to set the bounds to include all iterated markers from the database. Can anyone help? I'm guessing I need to keep an array of the markers and update the bounds based on that?

var map = L.map('map').setView([51.505, -0.09], 13);

        L.tileLayer( 'https://cartodb-basemaps-{s}.global.ssl.fastly.net/light_all/{z}/{x}/{y}.png', {
          maxZoom: 18,
          attribution: 'CartoDB | &copy; OpenStreetMap contributors'
        }).addTo(map);

        $( document ).ready(function() {
          addCompanies();
        });

        function addCompanies() {
          for(var i=0; i<companies.length; i++) {
            var marker = L.marker( [companies[i]['latitude'], companies[i]['longitude']]).addTo(map);
            marker.bindPopup( "" + companies[i]['company']+"
Details:" + companies[i]['details'] + "
Telephone: " + companies[i]['telephone']);
          }
        }

        var companies = JSON.parse( '<?php echo json_encode($companies) ?>' );


texts: Why doesn't requested tile get cached in browser?

texts:
I'm trying to optimize the load time of a mapping project using QGIS server and Openlayers 4 and have a question about browser cache.

PNG tiles from tile data with URLs in XYZ format get cached correctly and show up in the Networking tab with Status Code 304 Not Modified.

But tiles from a WMS server get requested every time (Status Code 200 OK) and don't get cached by the browser. The request URL doesn't change and looks like that:

http://host.com/cgi-bin/ffr/qgis_mapserv.fcgi?map=/home/qgis/map.qgs&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&FORMAT=image/png&TRANSPARENT=true&LAYERS=ocean&WIDTH=256&HEIGHT=256&CRS=EPSG:3857&STYLES=&BBOX=-7514065.628545966,5009377.085697312,-6261721.357121638,6261721.35712164

Any chance to cache this tiles in the browser so they don't get requested each time from WMS server?


texts: Enable 3D viewer in geo explorer-geo server

texts: I am trying to enable 3D viewer(google earth) in Geo explorer-geo-server by installing google earth plugin but doesn't work . any ideas ?


texts: Extracting single area with multiple layers from two different data sources

texts:
I'm very new to GIS.

I am using ArcMap.  I have a state map of congressional districts from one data source.  From another data source, I have blank census tracts.  From yet a third data source (this one non-GIS data from a CSV file) I joined it to the census tracts, matching data to the census tracts.  So far it works great, I can create layers based on census data that matches these census tracts.  Here's the final issue:

I need a "clean" map of certain congressional districts, that shows the census tracts and all their layered data within just a single district.  If I try to clip a district, it is just the outline of the district.  The only solution I have come up with is to clip the district and carefully manually select ~20 census tracts, then I have to completely recreate all the layers to replicate my statewide map.  This would be insanely time consuming for all the districts I need to pull.  I also tried select by location but it just copied the district and ignored all the other data.

Any ideas or suggestions?


texts: How to query annotations from a GDB file in .Net

texts:
I have annotations in a GDB file that I need to extract to load into a non-ESRI spatial store.  It's my understanding that there isn't any tool to export annotations to a shape file.  I need to copy points, text, and style info.

I have not worked with annotations before.  My prior experience with ESRI products involves writing a shape file parser in C#.Net and writing .Net code that interacts with the C API for ArcSDE to query ArcSDE.  It's my understanding that ArcSDE SE_annotation methods don't really return very much detail about annotations, so I would be open to working with ArcObjects as well.  I've been reviewing that API, however, and I'm not clear on how I would query annotations.

If anyone has experience with this sort of thing and can suggest ideas or point me to relevant documentation I'd be very grateful.


texts: How to calculate area using QgsDistanceArea?

texts:
I am trying the following code as given in the QGIS documentation:

d = QgsDistanceArea()
d.setProjectionsEnabled(True)

print "distance in meters: ", d.measureLine(QgsPoint(10,10),QgsPoint(11,11))

But getting the following error:

  " AttributeError: 'QgsDistanceArea' object has no attribute 'setProjectionsEnabled' "

What is wrong here?


texts: Is it possible to read a subset of points with the PDAL C++ api?

texts:
I am using the PDAL C++ API to read LAS/LAZ data. We're finding some of the tile files are very large (160M+ points), so I wish to take a sample of them. Although the point order in the tiles is arbitrary, it definitely isn't random, so we probably want random sampling. Alternative methods would appear to be too slow and/or require an index. Index building is the specific limit I'm hitting with 160M points!
I could easily do this myself, but it looks like PDAL reads all points into memory first. This seems wasteful for such large datasets.
Is it possible to tell PDAL to "only read a random selection of N valid points from disk"?


texts: How can I smooth a discrete raster in ArcGIS?

texts:
I started with a continuous spline raster of hotspots, classified by standard deviations from the mean. This displays as 7 classes, but since it's a continuous raster in ArcGIS 10.2, it appears nice and smooth. However, my goal is to take the splined continuous raster that looks nice and smooth, and make it appear 3D using cartographic representations. That means I needed to convert this spline to vector.

Now, this next part is where I'm not sure about my approach, is there are better way to do this? What I did was I ran the raster calculator on my splined continuous raster for each class. I did calculations like where spline_raster >=2.58, for example. I want to do this for hotspots and coldspots, so I have 6 classes I need to extract through raster calculator. I did this, and raster calculator converted the output into discrete rasters. So where my spline was >= 2.58, it was now 1, and everything else = 0. The issue wasn't so much this, but that the output looked much more "stairstepped" or pixelated compared to my continuous spline. I understand the difference between discrete and continuous, but I'm wondering if there's a best practice for at least making the discrete raster's pixels so much smaller that at the same scale as the continuous spline, they appear almost as smooth? I tried going into the environment settings and decreasing the cell size from the original 1414 to 5, and it hardly changed how the raster calculator output looked. I don't understand this part very well. I researched this a little, and figured I would try resampling my rasters. I tried resampling using nearest neighbor, bilinear, and cubic...they all still look the same. My conclusion is that there is some part of this process I don't quite understand.

Since I would then be vectorizing the raster calculator output (using Raster to Polygon tool), the vectorized polygons look all jagged because of the stairstepping, and I have been using the Smooth Polygon tool to make them look smoother. I feel this approach is very roundabout and not the best way to do it, but I don't know a better way. Because of my field of research and the types of data I've used, I don't typically do a lot of raster analysis.


texts: How can I make a reference plane to measure surface volume?

texts:
I am looking to make a reference plane like the links below. I want to measure surface volume of an area in different elevation layer.

I just need to know how to make the reference plane in ArcGIS 10.1 !

http://resources.arcgis.com/en/help/main/10.1/index.html#//00q900000037000000

Best,
Mahnoush


texts:
I want to use eVis but cannot find it in the plugin manager. I know it was going to be merged into the core of QGIS.
Has it changed name and location?


texts: Error when using r.series in a python script

texts:
I’m writing a tool script in python and at a certain point I want to use the GRASS tool r.series. I want the algorithem to read all raster files from a folder and calculate the maximum with r.series. The thing is I’m getting the error: “Wrong parameter value: None” when using the code below:

##RasterMaxFinder=name
##Select_directory=Folder
##Save_results=Folder

import glob, os, processing

list = ( Select_directory + "/" + "*.tif" )

processing.runalg('grass7:r.series', list ,False,6,'-10000000000,10000000000',None,30.0, Save_results + "/" + "Output.tif")

It must have something to do with the extend, but I can't figure out what. When I'm using the tool directly in QGIS, erverything works fine. Any Ideas?


texts:
I've been using the WmsLayer for my WMS Service needs, rather than my old way of doing things by sussing out a GetMap URL and posting that image to a GraphicsLayer. It's been working well but I seem to have found a service that breaks it.

Here is the getCapabilities URL I'm using:

http://ags.pdc.org/arcgis/services/global/pdc_active_hazards/MapServer/WMSServer?SERVICE=WMS&REQUEST=GetCapabilities

I'm able to to get the Name of the layers available (they're using simple numbers to ID them). I provide the WmsLayer the  getCapabilitiesURL, the layer ID number & a Proxy Url based upon the proxy handler from the documentation here at ESRI. Unfortunately that doesn't seem to be enough. When debugging, I crash to an unhandled exception as soon as I add it to my map control.

Category: ManagedRuntimeError
Message: System.InvalidOperationException: Sequence contains no elements

Using Firebug to look at my hosted application, it reports this error:

Unhandled Error in Silverlight Application Data at the root level is invalid. Line 1, position 1.   at System.Xml.XmlTextReaderImpl.Throw(Exception e)
   at System.Xml.XmlTextReaderImpl.Throw(String res, String arg)
   at System.Xml.XmlTextReaderImpl.ParseRootLevelWhitespace()
   at System.Xml.XmlTextReaderImpl.ParseDocumentContent()
   at System.Xml.XmlTextReaderImpl.Read()
   at System.Xml.Linq.XDocument.Load(XmlReader reader, LoadOptions options)
   at System.Xml.Linq.XDocument.Parse(String text, LoadOptions options)
   at ESRI.ArcGIS.Client.Toolkit.DataSources.WmsLayer.client_DownloadStringCompleted(Object sender, DownloadStringCompletedEventArgs e)
   at System.Net.WebClient.OnDownloadStringCompleted(DownloadStringCompletedEventArgs e)
   at System.Net.WebClient.DownloadStringOperationCompleted(Object arg)

The XML returned by the GetCapabilities call seems sound to me - at least it appears like all of the other WMS services I've been communicating with - so I'm not sure what parsing errors it is running into. Any ideas? I can provide further info if there is something I haven't thought to check.


texts: How to determine which EditEvent is fired BEFORE the feature is created?

texts:
I have written an AddIn-Tool and want to handle the OnSketchFinished-Event.
The Way

((IEditEvents_event)_myEditSketch).OnSketchFinished += OnFinishedSketch;

doesn't really work for me, cause I need the get the event before the feature is created.
I rather need something like the OnFinishSketch that was fired by the IEditTask.

Does anyone have an idea how to solve this problem? I'm very thankful for any suggestions.


texts:
Since we know about convolution filter in image processing, it is basically a matrix applied to an image and a mathematical operation comprised of integers and it works by determining the value of a central pixel by adding the weighted values of all its neighbors together.
Question: In Envi software, we apply a convolution filter under:-
the Filter  >  Convolutions and Morphologytab
Now there we have a default 3×3 LPF (low-pass filter) and HPF(High pass filter), On the basis of their usage/kernels, how exactly these filters work?


texts:
Does anyone one know where I can get the Geoportal Server 1.1.1 search widget for Flex 2.5?

I have been tasked with installing this to an existing geoportal server and Flex viewer does anyone know how to do this. I have found some instructions but they seem to require reloading the viewer. Any Ideas or advice would be appreciated.


texts: GeoServer's "infinite loading" on new layer screen

texts:
I deployed a GeoServer 2.10.1 version on tomcat 8.x. then I created a new connection from a Postgis database and published a layer. But, when I tried to publish a new layer, the GeoServer's stayed on the "infinite loading" screen, and   the publish option not work properly.

Does someone has the same problem?


texts: Setting maximum slope/flow length in flow accumulation using SAGA?

texts:
I like to limit the slope length in SAGA GIS. In other words, the flow accumulation should only account for a specific length until a subsurface flow initiates  and surface flow stops. Actually, I have flow length up to 2000 m.

In SAGA GIS, I only find tools to calculate the flow length but no tool to set a maximum slope/flow length threshold.

Do you have any idea for a workaround?


texts: Why can't I enter decimals in QGIS r.water.outlet?

texts:
I am using the latest version of QGIS 2.14.  I have a raster DEM file. I ran the r.watershed and it created all the raster layers for flow direction, stream, etc.

My next step is to run r.water.outlet to subdivide a basin.  In the r.water.outlet input window, more specifically the northing and easting input text box, the built in filter does not allow me to enter numbers past the decimal point.

Example: for a point at -73 35 45.5 E when I tried to enter -73.35455 the most that I can enter is -73. and no numbers past the decimal point.  Am I doing something wrong?


texts: Is there a setting to force ignore invalid geometries when performing spatial analysis in QGIS 3.0?

texts: When performing spatial joins I get no results due to geometry errors in one of the layers, but there is a note saying to force ignore geometry errors in the log. Where can I activate this setting?


texts: How to use gdal_translate to transfer all 400 gcp's at once

texts:
With gdal_translate we can add GCP from one source to another source

gdal_translate -gcp pixel1 line1 X1 Y1 -gcp pixel2 line2 X2 Y2 ... source.tif sourceGCP.tif
gdalwarp -r bilinear -t_srs EPSG:xyz sourceGCP.tif sourceCRS.tif
gdalwarp -t_srs EPSG:4326 sourceCRS.tif destCRS.tif

Writing 400 or even more GCP's the format is time-consuming, is there any wild card in which we can mention to translate all the GCP's from source to destination at once?


texts: Managing outages to make ArcSDE updates?

texts:
I would like options to manage ArcSDE outages for making updates to schema, which I am not able to do when users are just reading the feature classes.  I would like to not have to stay late to get these changes and I would like to not have downtime for an hour.  I imagine well thought design would mitigate this.

I need to delete a few objects and create a few objects.


texts: Add table from SDE to ArcMap TOC via ArcMap python window

texts: I would like to add a specific table or shape file to my TOC via the ArcMap python window from the SDE connection I'm currently connected to.  Since I'm already connected to the SDE connection in ArcMap I'm assuming I don't have to restate the user name and password in the python script right? How would I going about doing this?


texts:
I'm wondering how to edit multiple layers, maybe standalone.

I have 2 shapefiles.
1 line and
1 point

What i want is the following:
When i select a line, then the line and the point shape has to go in de edit mode. The selected line changes from color (say red to blue). Now i want to insert serveral new points on the line, wich gets a id number from the line shape.

I think this is possible with some py code, but i can figure out how.

Any ideas?


texts: Is there an EPSG code to specify WGS84 in (lon,lat) order instead of (lat,lon)?

texts:
Is there an EPSG code like EPSG:4326 or EPSG:4329, but that specifies (lon,lat) order instead of (lat,lon)?

I've tried exploring the websites
http://spatialreference.org and
http://epsg-registry.org/

and done a bunch of googling, but no luck so far.

I'm interested in this because the Geoportal Server CSW I'm using assumes (lon,lat) but seems to ignore the CRS, while the PyCSW CSW I'm using assumes (lat,lon), but honors the CRS request.  Thus if I had an EPSG CRS for (lon,lat) order, I could specify this to both services and get uniform response.


texts: Feature linked annotations not updating when editing attributes in features?

texts:
I am using ArcGIS 9.3.

I have a copy of an original SDE database with feature classes, annotations, and relationship classes for those annotations. The copied geodatabase is an exact replica of an ArcSDE database called PRDGIS.SDE. Versions of this database are copied automatically overnight by python script and preserved for about a month, then deleted. The most recent version from our SDEarchive folder was pulled for this test. Relationship classes and topologies were removed. A geometric network was then created.  Using the blueprints of the original relationship classes, I re-created the relationship classes for all annotations after the network was built. I am testing to see if the features are linked as well as the annotations. Everything is linked up correctly (the right annotations belong to the right features), however, when I attempt to update attribute information for a feature and have the annotation reflect that change, the annotation simply disappears.  The annotations were related by a composite, forward direction, one to many relationship.


texts: How to estimate the average width of a sidewalk?

texts:
I have a table with two polygons representing urban block with and without the sidewalk:  geomBlock and geomBlockWithSidewalk.

On my city we can approximate the polygons of geomBlockWithSidewalk to a geomBlock buffer of 2 meters. The standard OGC function that generates this buffer is

geomBlockWithSidewalk_false = ST_Buffer(geomBlock,2.0)

It is "false" because reality is very very different than this ideal scenery. My problem now is estimate the sidewalk width w, that is How to estimate w on the best approximation for st_buffer(geomBlock,w)?

Here an illustration of typical urban blocks shapes and sidewalk variations.

.

NOTICE:

If the sidewalk representation are stored as "rectangular strips",  the problem of this question is similar to
"How can I calculate the average width of a polygon?", because both solutions need to solve the problem of "estimation of the average width of a rectangular strip".

Another usual representation for sidewalks, used here, is aggregating the shape of the sidewalks around a block with the block...

The rectangular strip here can be also a ring (a polygon with a hole): the sidewalk is the polygon resulted from the spatial difference between geometries of the blockWithSideWalk (g2) and the block (g1)  "g2 - g1", that is ST_Difference(g2,g1).

Imagining cut the ring: the cut make it a "rectangular strip"... But not exactly, it is more precisely a "trapezoidal strip".

It is a new question. After @whuber closed this question as "duplicated", I posted the general question: "Is there a st_buffer inverse function?", where we can see that all the questions are very closely related,

the "width of the beach of a little island" problem,
the calculation of "the average width of a street" problem, and
this (here) "average width of a sidewalk represented with the block" problem.

The solution can use the "hypothesis of retangular blocks with sides L and H >>w". It is a function,

CREATE FUNCTION sidewalk_width(geometry, geometry) RETURNS float AS $$
    -- calculations using geomBlockWithSidewalk and geomBlock
    -- ST_Within(geomBlock,geomBlockWithSidewalk) is true.
$$ LANGUAGE SQL immutable;


texts: River flow accumulation: % upstream deforestation area as input

texts:
I am trying to create a 'flow accumulation' raster whereby each river pixel tells me the % upstream area that has been deforested.
This screenshot shows the HydroSHEDS Flow Accumulation dataset. I have drawn four circles (ignore big red circle) around four pixels at the start of one river. So each pixel essentially has its own catchment area, with a certain # of 30x30m deforested pixels (blue in this case) inside it. Is this possible in GEE?


texts: Update column/Join in QGIS similar Mapinfo

texts: Looking somewhat similar to UPDATE COLUM option ( for join) MapInfo in QGIS. This exists ?


texts: Merge records with identical ID in ArcGIS Desktop?

texts:
I need to merge records with the same PIN, and populate the null fields with values if they exist. The highlighted ones in the image are what I'm talking about.

This is just a test run; for the final product I will have 2,000+ records and 14 additional fields, instead of 4.

Does anyone know a good workflow to accomplish this?

If my description is not clear enough, this is what I need.

No more duplicate records, each parcel should have one record, and PIN 35.00071.01 (the first highlighted record) should read:

P12_acres: 79 | P12_amount: 10666 | P17_acres: 160 | P17_amount: 24000


texts: Split/Cut polygons within feature classes based on a polygon mask

texts:
I am looking for a way to split features based on a mask. The end goal is the deleting of the features within the mask, without creating an output.

The polygons are bigger than my intented mask, but I just need to delete the area within the mask. Therefor I cannot delete based on location or attribute. By splitting these polyongs based on the mask, I can then select them and delete them. Afterwards I want to "paste" new data in the deleted area.

I am looking for a way to do this within modelbuilder or Arcpy, though my knowledge of Arcpy is beginner at best.

I am using ArcGis 10.3.1.

Clarification of the question:


texts: Correcting Blender Error in ArcGIS Desktop?

texts:
I have a map featuring 4 data frames with the same/similar files in each data frame: esri world imagery basemap, and various shapefiles and feature classes pulling from the network.  I reopened the mxd recently and I'm getting the following drawing error:

One or more layers failed to draw:

Blender Error. Not enough memory to draw the map. Release memory by ungrouping base map layers or decelerating raster basemap layers.
Blender Error. Not enough memory to draw the map. Release memory by ungrouping base map layers or decelerating raster basemap layers.
Blender Error. Not enough memory to draw the map. Release memory by ungrouping base map layers or decelerating raster basemap layers.
Blender Error. Not enough memory to draw the map. Release memory by ungrouping base map layers or decelerating raster basemap layers.

The map looks correct in the mxd, however when I export to a PDF, the basemap is incorrect.

How do I remedy this? I'm not sure how to release memory or decelerate basemap layers and I have nothing grouped.

Using ArcMap 10.2.2


texts: Converting SQL geometry data queries into R spatial dataframes

texts:
I'm interested in querying data with a geometry field from SQL server and using it in R as a spatialPolygonsDataFrame.

I've seen plenty of examples of using rgdal to import a shapefile, but shapefiles are far from the ideal storage medium for a map layer.

Can anyone point me to the correct syntax?

Does the geometry data have to be converted on it's way out of sql server in order for R to recognize it? to WKT or WKB for example?

In the end, I want a live connection to this data so that my maps are pulled directly from SQL and not from an exported shapefile.


texts: Converting Shapefile to GeoJSON. Both OGR and Mapshaper converters give errors

texts:
I'm trying to convert a large Shapefile (120 MB) to GeoJSON.

OGR converter gives an error that says process time took over 15 seconds, so process was cancelled.

Mapshaper gives the following error: File size of .shp doesn't match size in header

What can I do now?  I'd prefer not installing any desktop software to do a simple conversion.


texts: British National Grids projection does not add Ordnance Survey grids to shapefile?

texts:
I have created a shapefile in ArcCatalog by right-clicking the folder I want it in, then selecting new, then shapefile, then changing the type to "polygon", and using the British National Grids under the projected co-ordinate system.

I then tried to add the shapefile in ArcMap by clicking the Add Data button and selecting the shapefile I created.

My problem is that when I open the shapefile in ArcMap, the map is still blank and I cannot see the projected shape of the British National Grids.


texts: Adding points at certain distances along lines

texts:
I have encountered the following problem: I have a set of lines in ArcMap. I need to place points at specified (not regular) distances along these lines. The distances are not the same between lines.

I feel its a waste of time to manually measure each distance and Place a Point there.

Any tips for a quicker way?


texts: Can't cut raster with a polygon shapefile

texts:
I'm trying to cut a CLC raster (several, actually) with a shapefile (polygon) boundary from a country (Portugal, continent only), bu i keep getting this error:

QGIS version: 3.8.2-Zanzibar
QGIS code revision: 4470baa1a3
Qt version: 5.11.2
GDAL version: 2.4.1
GEOS version: 3.7.2-CAPI-1.11.0 b55d2125
PROJ version: Rel. 5.2.0, September 15th, 2018
Processing algorithm…
Algorithm 'Clip raster by mask layer' starting…
Input parameters:
{ 'ALPHA_BAND' : False, 'CROP_TO_CUTLINE' : True, 'DATA_TYPE' : 0, 'INPUT' : 'D:/d_araujo/Documents/SIGs/Dados/CLC/CLC_CHA_12_18/clc2018_cha1218_v2018_20_raster100m/CLC2018_CHA1218_V2018_20.tif', 'KEEP_RESOLUTION' : True, 'MASK' : 'D:/d_araujo/Documents/SIGs/Projetos/UTAD/PROJETO_FINAL_QGIS/DADOS/SHAPE/boundary_continent.shp', 'MULTITHREADING' : False, 'NODATA' : None, 'OPTIONS' : '', 'OUTPUT' : 'TEMPORARY_OUTPUT', 'SET_RESOLUTION' : False, 'SOURCE_CRS' : None, 'TARGET_CRS' : QgsCoordinateReferenceSystem('EPSG:3035'), 'X_RESOLUTION' : None, 'Y_RESOLUTION' : None }

GDAL command:
gdalwarp -t_srs EPSG:3035 -of GTiff -tr 100.0 -100.0 -tap -cutline D:/d_araujo/Documents/SIGs/Projetos/UTAD/PROJETO_FINAL_QGIS/DADOS/SHAPE/boundary_continent.shp -cl boundary_continent -crop_to_cutline D:/d_araujo/Documents/SIGs/Dados/CLC/CLC_CHA_12_18/clc2018_cha1218_v2018_20_raster100m/CLC2018_CHA1218_V2018_20.tif C:/Users/d_araujo/AppData/Local/Temp/processing_fdc14054ef044d2ebbef8d7e8da929b0/50223fb5352f4f44a621aba9eea2a396/OUTPUT.tif
GDAL command output:
Copying raster attribute table from D:/d_araujo/Documents/SIGs/Dados/CLC/CLC_CHA_12_18/clc2018_cha1218_v2018_20_raster100m/CLC2018_CHA1218_V2018_20.tif to new file.

Creating output file that is 3415P x 5686L.

Processing D:/d_araujo/Documents/SIGs/Dados/CLC/CLC_CHA_12_18/clc2018_cha1218_v2018_20_raster100m/CLC2018_CHA1218_V2018_20.tif [1/1] : 0ERROR 1: IllegalArgumentException: point array must contain 0 or >1 elements

ERROR 1: Cutline polygon is invalid.

Using internal nodata values (e.g. 0) for image D:/d_araujo/Documents/SIGs/Dados/CLC/CLC_CHA_12_18/clc2018_cha1218_v2018_20_raster100m/CLC2018_CHA1218_V2018_20.tif.

Copying nodata values from source D:/d_araujo/Documents/SIGs/Dados/CLC/CLC_CHA_12_18/clc2018_cha1218_v2018_20_raster100m/CLC2018_CHA1218_V2018_20.tif to destination C:/Users/d_araujo/AppData/Local/Temp/processing_fdc14054ef044d2ebbef8d7e8da929b0/50223fb5352f4f44a621aba9eea2a396/OUTPUT.tif.

Execution completed in 1.52 seconds
Results:
{'OUTPUT': 'C:/Users/d_araujo/AppData/Local/Temp/processing_fdc14054ef044d2ebbef8d7e8da929b0/50223fb5352f4f44a621aba9eea2a396/OUTPUT.tif'}

Loading resulting layers
Algorithm 'Clip raster by mask layer' finished

I feel the problem relies only on the boundary shapefile, but in ArcGIS i've used it sucessfully for various tasks.

I tried Check validity on the shapefile, and with GEOS method it doesn't return any error, and with QGIS method, it just freezes, and doesn't produce any result. I tried Fix geometries, but it returns nothing useful. I tried to buffer with a 0 m distance, and this shape also produces the error above. When i create a new shapefile with a simple rectangle on it, i can cut the CLC with it, so the raster doesn't seem to be the issue.


texts: Selecting maximum Value based on other field using ArcMap?

texts:
I have problems with a SQL query. The table looks like this, it is a point feature class both in postgres and file geodatabase, I am working in ArcMap.

name: elev_p
HEIGHT AREAid
100    1
200    1
300    2
500    2
100    1

So, what I want is the highest HEIGHT within each AREA_id. If there are 200 points within each polygon, I only want one, and the one with maximum HEIGHT.

Is this possible in ArcMap or do I have to use SDE layer and pgAdmin?

This AREA_id comes from polygon layer which I joined with elev_p, I tried other way around and I got the maximum HEIGHT but I need to know which exactly point it is, so its selection I need.

WHAT I HAVE TRIED: because one copy (sde version) can handle subquery I used HEIGHT in (SELECT MAX(HEIGHT) FROM elev_p group by AREA_id). But that query selected all Height values in entire feature class which had the highest value in AREA_id, not one for each AREA_ID.

SOLUTION: SELECT * FROM elev_p a INNER JOIN (SELECT MAX(HEIGHT) as bheight, AREAid FROM elev_p GROUP BY AREAid) b ON a.AREAid = b.AREAid AND a.HEIGHT = b.HEIGHT;

note: this is for spatial data in database like oracle or postgres, arcmap solution I do not have

also, if there is 2 points with same value wich is highest within AREAid then it will choose both of them, in case someone needs only exactly one


texts: st_crop() follow latitudinal circle on polar projection

texts:
I am generating a polar projection with the rnaturalearth::ne_countries() dataset. Because I'm only interested in high-latitude regions, I want to crop the shape above a particular latitude. However, when I crop, transform, and plot, the southernmost extent of each country that remains in the dataset does not follow latitudinal curves, because only two points represent the southern bounds of the data. For example, the US west coast and east coasts are connected by a straight line, as highlighted in the image below (red rectangle added manually). Ideally, I want the southern "borders" of countries that were cut in half to follow the thin blue dotted line (added manually).

How can I get a latitudinal crop + transform to follow curves? Spatial* or sf* solutions are both welcome.
Reproducible example:
library(sf)
library(tidyverse)

countries = rnaturalearth::ne_countries(returnclass = "sf") %>%
  st_crop(y = st_bbox(c(xmin = -180, ymin = 35, xmax = 180, ymax = 90)))
ggplot() +
  geom_sf(data = countries)

countries2 = st_transform(countries,
             crs("+proj=laea +lat_0=90 +lon_0=0 +x_0=0 +y_0=0 +ellps=sphere +units=m +no_defs"))

ggplot() +
  geom_sf(data = countries2)


texts:
Is there anyway to add together GetCount Results? I am using GetCount to find the number of selected features for multiple layers but when I try to add them all together I get an error message saying it is an unsupported operand type for +: 'Result' and 'Result'.

SampledFaild = arcpy.GetCount_management('SamplingFaild')
SampledValid = arcpy.GetCount_management('SamplingValid')
Total = SampledFaild + SampledValid


texts: Does proj4js is provided with an EPSG dataset?

texts:
Is each version of the proj4js project delivered with/against a different EPSG dataset version? (like in the PROJ project)
In other words, if I'm using a specific proj4js version, does it mean that I'm also using a certain EPSG dataset version?
If so, where can I find the reference to it?


texts: How to write multiple .las files from polygon clipping with PDAL?

texts:
I am using PDAL (through Python) to crop multiple polygons from a point cloud and store them as separate files. For efficiency reasons I would like to define this in a single pipeline. However, using the approach shown below, PDAL only seems to store the data related to the last defined polygon. Is it possible to write multiple files with a single pipeline?

pipeline_definition = {
    'pipeline': [
        {
            'type': 'readers.las',
            'filename': input_file,
            'tag': 'read'
        },
        {
            'type': 'filters.crop',
            'polygon': poly1.wkt,
            'inputs': ['read'],
            'tag': 'poly1'
        },
        {
            'type': 'filters.crop',
            'polygon': poly2.wkt,
            'inputs': ['read'],
            'tag': 'poly2'
        },
        {
            'type': 'writers.las',
            'filename': 'branching_1.las',
            'inputs': ['poly1']
        },
        {
            'type': 'writers.las',
            'filename': 'branching_2.las',
            'inputs': ['poly2']
        }
    ]
}

pipeline = pdal.Pipeline(json.dumps(pipeline_definition))
pipeline.validate()
pipeline.execute()


texts: Measuring area using correct CRS and units in QGIS

texts:
I have uploaded a shapefile and I am simply trying to calculate some simple areas of provinces (and I know their areas from Wikipedia). My map is from Zimbabwe, and I have saved the project, reprojected all layers, and saved everything as WGS 84 / UTM Zone 36S. The units are all in metres.

When I use $AREA function in field calculator, the areas come around very very small. I expect areas in order the thousands of square kilometers, but I get numbers like 3m2.

Why are my results so far out?

Edit: I've added images of the property settings


texts: Spatial Join - Definition of 'Polygon is within X distance to point"?

texts:
What is the exact definition used by ArcGIS?
Does the entire polygon need to be within X to the target point? Or only at least point within the polygon?


texts:
I have recently downloaded and installed the OpenGeo Suite from Boundless, and I am running through a few tutorials. I'm using GeoExplorer to manipulate layers defined in GeoServer (running from my localhost) which is in turn connected to a PostGIS db (again stored locally).

I am successfully logged on as the admin and I have successfully published my map in GeoExplorer to a html page.

I want to now export it to a pdf, but my print button is greyed out.

Does anyone know why the print button is greyed out?


texts: I am using the extract by mask tool to fit my raster to a mask (province of Nova Scotia). I've noticed that the raster doesn't quite fit the mask; the north and south boundaries of the province is cut-off. Is there a way I can extend the boundaries of the raster to fit the mask? Below is a picture of the problem:


texts: In ESRI Collector, is there a way download/sync a layer that has attachments without attachments, but still allow for the adding of them?

texts:
We have data collectors using ESRI collector all across a state.  Every morning they have to download a new work area - World Topo basemap, plus slice of one layer (attribute data focused), and a slice of another layer (picture/attachment focused).

We wanted to try and eliminate the download portion each morning, giving our teams and extra 15-30 minutes every day.  We found we could create a basemap of the entire state at a high enough resolution (Problems using REST API to generate Tile Packages from World Topo (For Export) ), and we can download the attribute specific layer for the entire state (takes 5-10 minutes) but we can't download the picture focused layer (pics are tie to the attribute layer through IDs we create and maintain / geographic location), because of the sheer size (GBs) of pics that would need to be downloaded. Is there a way to download an attachment enabled layer within collector, such that:
We don't download any of the attachments when defining the work area and doing the initial pull/sync
But can still upload new attachments to the points created in the layer that get synced at the end of the day?

Anyone know?


texts: Converting raster tif to point shapefile using Python?

texts:
I've been working on a data that is in .tif raster format. I need to convert that into a points shapefile (.shp) to proceed. I went through GDAL in Python but there's only gdal.Polygonize(). So how do I convert it into a points .shp?
I've used this code so far:

from osgeo import gdal,ogr
import sys
import os
gdal.UseExceptions()
os.chdir(path)

src_ds = gdal.Open("IND_ppp_2015_v2.tif")

if src_ds is None:
print ("Unable to open Worldpop data")
sys.exit(1)

try:
    srcband = src_ds.GetRasterBand(1)
drv = ogr.GetDriverByName("ESRI Shapefile")
dst_ds = drv.CreateDataSource(dst_layername + ".shp")
dst_layer = dst_ds.CreateLayer(dst_layername,srs = None)

gdal.Polygonize(srcband, None,dst_layer, -1,[],callback = None)

Can you suggest a way to convert into points .shp?


texts:
I am trying to separate a multipolygon into equal parts. I am aware of the divide polygon plugin and that worked well for me in previous examples where a single polygon is involved. Unfortunately I have a multipolygon with multiple parts that are not connected. When I dissolve the multipolygon it shows me only one attribute in the attribute table, so far so good. But when I try and split that up into equal parts it splits it up irregularly based on the different parts of my polygon.
Is there a way I can split that multipolygon with its different parts into equal sized sections? All I found online was focused on connected single polygons.


texts: Is it possible to retrieve the data directory path via the REST API?

texts:
I'm writing a client for the GeoServer REST API and I want to create a data store for a directory of shapefiles. Since the REST API does not offer this functionality (only stand-alone shapefiles can be uploaded), I want to provide the functionality in my own client.

To do that, I need to know the path to the data directory, which is determined by GeoServerResourceLoader.lookupGeoServerDataDirectory(...). I've looked at the documentation and source code, and it doesn't seem to be exposed by the REST API. Have I missed something or is it indeed not possible to retrieve the path?


texts: Separate management of data and metadata policy in GeoNetwork

texts: I am working on my GeoNetwork and I am filling an ISO19139 for a spatial dataset. I associated the zip file of the dataset as “associated resource”. The metadata is public and automatically the zip file too so I can download  metadata and zip file too. I need to manage the download of the zip file with a separately data policy from the the metadata (restricted for data, open for metadata). Is it possible?


texts: Can QGIS 3.0 be installed and run side-by-side with QGIS 2.x?

texts: I'm interested to know if QGIS 3 can be installed and run side-by-side with a QGIS 2.x installation on windows?


texts: How to get a total amount of cells per classification type in a raster in ArcGIS 10.1?

texts:
I did an supervised classification for for this watershed.  The attribute table doesn't have the count for total number of cell for each classification.  I've created a field titled 'Count' with type 'Long', but I haven't figured out how to populate it.  The photo is just to better illustrate what I'm trying to do.

-Using ArcGIS 10.1

-Supervised classification done on Landsat7 images that had multiple bands.  I don't know why the table has a red, green, and blue column


texts: Using QgsExpression in a custom validation?

texts:
We use our own validation in QGIS,
which can be called from different areas of the plugin. In addition, reports can be generated from the results.

Now we have switched to our own attribute forms in QGIS and wanted to integrate our validation there as well.
That works fine, but sometimes the QGIS validation gets in the way. Therefore, we wanted to integrate the QGIS validation into our process.

We've written a wrapper that converts the QgsExpressions and constraints of a field into our rules and calls them to process the result.

Here is a small example of how we want to use the QgsExpression:

def validateImpl(self) -> bool:
    utils = QgsExpressionContextUtils()
    scope = utils.formScope(self.feature)
    context = QgsExpressionContext()
    context.appendScope(scope)
    # scope.setFeature(self.feature)
    self.expression = QgsExpression(self.expressionStr)
    value = self.expression.evaluate(context)
    return bool(value)

Unfortunately, it does not work self.expression.evaluate always returns 0 as a result, even though the validation would definitely fail.

I think that's because when we change a TextField in the attribute form and call validation on the signal textChanged (), then the actual attribute value in the feature has not been set, in other words, a change of the attribute does not seem to be noticed by QgsExpression.

If I set the attribute in the feature for every textChanged (), then the validation works:

self.feature.setAttribute (attribute, value)

So if the rule says that the attribute 'name' must contain the value 'QGIS', and setAttribute was called with (name = 'QGIS'), the rule would evaluate positively.

If you now close the attribute dialog with 'Ok' and save the layer, the value will not be saved.
For example, if the value 'Test' was entered as 'name' and has now been changed to 'QGIS', then after saving and reopening the attribute dialog, 'Test' will still be displayed.

If I remove the code (setAttribute), the validation will not work anymore, but the value will be saved back to the feature.

At the moment I do not understand what the problem is, how can I use a QgsExpression which checks the currently entered value in the attribute form?

I've seen that there's a function "current_value" in the scope I get from utils.formScope (self.feature), but I have no idea how to use it.


texts: Get distance in meters between lat/lng coordinates, divided in North and East direction

texts:
I've got 2 couples of lat/lng coordinates and I would like to get distance between them in meters, but divided in X (east) and Y (north) directions.

I've found this function in a  Stack Overflow answer to get distance between 2 points, but I don't know how to separate the North and the East distance:

var rad = function(x) {
  return x * Math.PI / 180;
};

var getDistance = function(p1, p2) {
  var R = 6378137; // Earth’s mean radius in meter
  var dLat = rad(p2.lat() - p1.lat());
  var dLong = rad(p2.lng() - p1.lng());
  var a = Math.sin(dLat / 2) * Math.sin(dLat / 2) +
    Math.cos(rad(p1.lat())) * Math.cos(rad(p2.lat())) *
    Math.sin(dLong / 2) * Math.sin(dLong / 2);
  var c = 2 * Math.atan2(Math.sqrt(a), Math.sqrt(1 - a));
  var d = R * c;
  return d; // returns the distance in meter
};


texts: Cross-validation error threshold in ArcGIS

texts:
I'm learning the Geostatistical Wizard toolset in ArcMap 10.4.1 and I am performing the cross-validation step in an IDW interpolation.

I have done some research regarding the process of cross-validation, however I can't find any information regarding the threshold of the prediction error
in the auto-correlation model.

What is the threshold of the prediction-error, which excludes sampling points in the auto-correlation model?


texts:
I've been using the Geometric Network tools for quite sometime now to trace electric distribution networks for my job and all of a sudden when try to trace the circuits and connected features from each particular transformer the trace won't work and the features aren't selecting. I have it set to Find Connected and the box in the Analysis dropdown and options is set to Selected Features. I can't figure out what has gone wrong. Also there hasn't been any ArcMap update since I'm still using 10.3.1.

Can anyone suggest a fix or know what's wrong?


texts: Filter and show images from an imageCollection in Google Earth Engine

texts:
once filtered an imageCollection in Google Earth Engine, how can I add as layers the images resulting by the filter process?

I need to show all the images (S2_L2A) available for a specific point (I've already filetered the collection using a geometry) in a specific time range (aready done) and with a CLOUDY_PIXEL_PERCENTAGE less than 5%.

this is my script

var PNA: Point(15.90,38.12)
var S2L2A: imageCollection"Sentinel-2 MSI: MultiSpectral Instrument, Level-2A"
var start = ee.Date('2018-01-01');
var finish = ee.Date('2019-10-01');
var imageSet = S2L2A.filterBounds(PNA).filterDate(start, finish).filterMetadata('CLOUDY_PIXEL_PERCENTAGE','less_than', 5);
var count = imageSet.size();
print ('size of collection imageSet', count);

Now...How can I add as layers the resulting images?


texts:
I have a layer loaded in QGIS and I would like to change the color. I've seen some code but for the previous versions and they don't work in QGIS 3.4.

How can I do this?


texts:
I am having some issues with adding xy data from a .csv file. As you can see in the picture below the points are clustering. I am not sure if this is a excel format issue or a arcmap issue.

I am simply adding the .csv to arcmap> display xy data> no z field and coordinate system of: Geographic Coordinate System: Name: GCS_WGS_1984


texts: Building a DEM with both a raster and vector points

texts:
So I'm struggling here, trying to get a decent DEM for the Azores Islands because I want to include bathymetric data on it.

I have a raster with elevation for the islands (land) with a 30 meter resolution, which is actually more detail than I would ask for but, unfortunately, the bathymetry I have managed to acquire are vector points 3000 meters apart.

I converted the raster into contour lines and tried to make an interpolation using both elements (points and lines), needless to say I run into problems around sea level.

My best results were obtained using a spline interpolation with different smoothing for lines and points on GRASS, but I get artifacts whenever the distance between a vector and a line is close to 3000 meters.

Any tips on how to get a decent DEM when using data of different types and with such different resolutions?


texts: Python script: Access data from databases

texts:
I need to access data that is inside an Oracle database using a Python script. It's a Unix server running oracle 10g. I don't have much access to the database server.

It works fine on model builder when the connection is open and I add my data, but exporting it to Python means all the connection properties are gone (they are just "TABLE NAME" etc.)

Basically the script grabs this information based on a query, let's say "SELECT id WHERE Status = 0", and it writes it to a new table in a filegeodatabase. I need to do it through python as I am processing it further, but that's not my issue.

I just can't get it to connect. I've tried to use pyodbc plugin using the syntax below:

cnxn = pyodbc.connect('DRIVER={SQL Server};SERVER=rtsprod;DATABASE=DBNAME;UID=USERID;PWD=PASSWORD')

Now I know it's an Oracle database but can't find any documents on how exactly to connect using the pyodbc plugin. Is there a better way?


texts: Massachusetts Woodland Shapefile / GDB Not Rendering Correctly QGIS

texts:
I am working on a personal project that utilizes USGS Land Cover Woodland data sets and I seem to have an issue in QGIS where part of the state of Massachusetts is missing / not rendering correctly. The first image is what QGIS renders and the second image is the sciencebase.gov thumbnail of the layer. Both the Shapefile and the FileGDB have this issue.
This issue occurs in a new blank QGIS 3.16.3 project with EPSG:4269 projection.
https://www.sciencebase.gov/catalog/item/5b6fa9f4e4b0f5d57878e73d
https://www.sciencebase.gov/catalog/item/5b6fa9f3e4b0f5d57878e73b


texts: pythonaddins.OpenDialog doesn't work properly

texts:
I am developing a toolbar by the pythonaddins, in regards to define all the input raster files. I wrote the follwoing lines so the user enters how many data he wants to read and accordingly, an input box will be displayed

# Import
import arcpy, os,sys, pythonaddins

Folder_Location = arcpy.GetParameterAsText(0) # please select a folder to save the project files
ScN=arcpy.GetParameterAsText(1) # please write the number of scenarious to be compared
#-----------------------#
ScN=int(ScN)
for i in range (1,int(ScN)+1):
    MessageD='Select the raster layer of flood depth for scenario No   '#+ str(i)
    Depth_rstr[i]= pythonaddins.OpenDialog(MessageD,False,r'C:\','Select')

unfortunatly, the OpenDialog window just display folders without allowing me to double click it to select any file

hoowever when I write the same lines in python window, it works well

I notices when I operate this puthonaddins- button, I got this error in the python window

      TypeError: GPToolDialog() takes at most 1 argument (2 given)

Could any one advice me, how to solve this problem? I want this could to select files.


texts: Difference between postgis.net and postgis.org for PostGIS documentation?

texts:
What is the difference between two sites hosting different PostGIS docs?

For example:

http://postgis.net/docs/ST_Intersects.html
http://postgis.org/docs/ST_Intersects.html


texts: Creating concave polygon with constraint angles in QGIS

texts:
I'm having point data. By using that I want to create concave polygon with some constraint angles (i.e., >45 and <225 degree).

Is it possible to do in QGIS?


texts: Add multiple Layer files in subdirectories to mxd using python

texts:
I have more than 300 layer files located in many sub-Folders and sub-sub-Folders. All Sub Folders are located in one large directory.I read Use python to add layers to TOC
and I try, with arcpy, to detect all the layers and add them to mxd.
Is it possible?


texts: How can I further query the result of an Esri REST endpoint query in JavaScript?

texts:
I am querying Esri REST endpoint to get features from a dataset, but how can I further query the resultset that comes back without making a call back out to the server?

I make an initial js call very broad so that I get all of the features I could possibly need. Then on that data set, in JavaScript, I want to make a sub query

Do I use Esri syntax against the local dataset somehow?
Or is there a JavaScript query language?
Or do I blindly loop over the dataset?


texts: QGIS exports slightly inaccurate GeoJSON coordinates

texts:
I noticed in my map that when I export to GeoJSON, some smaller feature have their coordinates shifted. I have set the export COORDINATE_PRECISION to 99 (max), tried various coordinate systems but no luck. In QGIS my layer coordinate systems are 3857 and I export to 4326.

Here are the map in QGIS:

Here is the imported GeoJSON in Mapbox (the result is the same in geojson.io and others)

To eliminate project setting errors, I've started a fresh new project. Then added OpenStreetMap layer with coordinate system 3857. Then created a new shapefile layer with the same coordinate system and drew parallel rectangles. With just an export to 4326 GeoJSON, the rectangles get distorted.

I also exported via online Shapefile to GeoJSON converter using ogre.adc4gis.com to eliminate QGIS export function having a problem during projection, but the result is the same.

Edit 1:
Following Jakob Miksch's comment, I created a new shapefile layer with EPSG:4326. I exported it and in most cases this fixed the issue. However, minor inconsistencies are still there as below. The inconsistencies mostly happen in certain areas of the building.

Edit 2:
Here is the shapefile I've been working on: Shapefile

Edit 3:
Switching the coordinate system to custom system centered on the location I'm working on seems to help a little bit with the output, but I still see interest challenges. For example, here, all these boxes share an edge (from a guideline), but the box above doesn't have the middle vertex, and therefore is straight, while, the top vertex of the smaller boxes below utilize a slightly off vertex. Interestingly, if I recreate the box above, but add a vertex, it removes those ugly overlaps (still not everything is perfectly parallel, however).

Edit 4:
Here is my latest GeoJSON file. It utilizes the custom coordinate system as recommended here, which minimizes the errors significantly. However in other renderers outside of QGIS (geojson.io, MapBox), the edges won't remain parallel/perpendicular.


texts:
I was given 8 files that supposedly contain GIS data. They appear to be a file split into parts and are named: "stuff.zip [1_8].dat" , "stuff.zip [2_8].dat" , etc..

I tried stripping the ".dat" but that didn't help to unzip the files. I know this is probably a zip problem rather than GIS problem but I was hoping someone might recognize how the heck someone would send GIS data in this form. Opening with a text editor shows the following type of data, any idea if this is an old ESRI file type that was given the wrong extension or something? Thanks for any ideas.

  UEsDBBQAAgAIAJlIQTmhyPOiMJdSALDsUgAXAAAAU2V3ZXJfSW52ZW50b3J5X01hcC5wZGYAM4DM
  fyVQREYtMS41DQoxIDAgb2JqDQo8PA0KL0xlbmd0aCAyIDAgUg0KL0ZpbHRlciAgL0ZsYXRlRGVj
  b2RlDQo+Pg0Kc3RyZWFtDQp4nOy8WbOqSNcuej8j/A8LUBBQ6QUVe7ABGxCxp7FB7BGbc2LXxfnt

Edit: it looks like this is a problem with how the data was encoded for the email message? After searching around a bit, I'm not sure how to deal with this. Obviously this is in the completely wrong forum at this point but thanks for any help. Here is a header before the alphanumeric string begins:

  Content-Type: application/octet-stream;
      name="stuff.zip"
  Content-Transfer-Encoding: base64
  Content-Disposition: attachment;
      filename="stuff.zip"


texts: 1km circles around lat-long-points in many places in world

texts:
I have hundreds of lat-long-points spread around the world, and must create circle-polygons around each of them, with radius exactly 1000 meters. I understand the points must first be projected from degrees (lat long) to something with meter-units, but how can this be done without manually searching and defining the UTM-zones for each point?

Here is a mwe for first point in Finland.

library(sp)
library(rgdal)
library(rgeos)
the.points.latlong <- data.frame(
  Country=c("Finland", "Canada", "Tanzania", "Bolivia", "France"),
  lat=c(63.293001, 54.239631, -2.855123, -13.795272, 48.603949),
  long=c(27.472918, -90.476303, 34.679950, -65.691146, 4.533465))
the.points.sp <- SpatialPointsDataFrame(the.points.latlong[, c("long", "lat")], data.frame(ID=seq(1:nrow(the.points.latlong))), proj4string=CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))

the.points.projected <- spTransform(the.points.sp[1, ], CRS( "+init=epsg:32635" ))  # Only first point (Finland)
the.circles.projected <- gBuffer(the.points.projected, width=1000, byid=TRUE)
plot(the.circles.projected)
points(the.points.projected)

the.circles.sp <- spTransform(the.circles.projected, CRS("+proj=longlat +ellps=WGS84 +datum=WGS84"))

But with second point (Canada) it doesn't work (because wrong UTM-zone).

the.points.projected <- spTransform(the.points.sp[2, ], CRS( "+init=epsg:32635" ))

How can this be done without manually getting and specifying UTM-zone point per point? I don't have any more info per point than lat long.

Update:

Using and combining the great answers from both AndreJ and Mike T, here is the code for both versions and plots. They are different on 4th decimal or so, but both very good answers!

gnomic.buffer <- function(p, r) {
  stopifnot(length(p) == 1)
  gnom <- sprintf("+proj=gnom +lat_0=%s +lon_0=%s +x_0=0 +y_0=0",
                  p@coords[[2]], p@coords[[1]])
  projected <- spTransform(p, CRS(gnom))
  buffered <- gBuffer(projected, width=r, byid=TRUE)
  spTransform(buffered, p@proj4string)
}

custom.buffer <- function(p, r) {
  stopifnot(length(p) == 1)
  cust <- sprintf("+proj=tmerc +lat_0=%s +lon_0=%s +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs",
                  p@coords[[2]], p@coords[[1]])
  projected <- spTransform(p, CRS(cust))
  buffered <- gBuffer(projected, width=r, byid=TRUE)
  spTransform(buffered, p@proj4string)
}

test.1 <- gnomic.buffer(the.points.sp[2,], 1000)
test.2 <- custom.buffer(the.points.sp[2,], 1000)

library(ggplot2)
test.1.f <- fortify(test.1)
test.2.f <- fortify(test.2)
test.1.f$transf <- "gnomic"
test.2.f$transf <- "custom"
test.3.f <- rbind(test.1.f, test.2.f)

p <- ggplot(test.3.f, aes(x=long, y=lat, group=transf))
p <- p + geom_path()
p <- p + facet_wrap(~transf)
p

(Not sure how to get the plot into the update).


texts:
I am using the Geotools library to develop a Java Swing application. I am having a problem creating a layer with lines between points. Can someone help me?

Here my code:

import com.vividsolutions.jts.geom.Coordinate;
import com.vividsolutions.jts.geom.GeometryFactory;
import com.vividsolutions.jts.geom.LineString;
import java.awt.Color;
import static org.geotools.coverage.processing.operation.Interpolate.TYPE;
import org.geotools.feature.DefaultFeatureCollection;
import org.geotools.feature.simple.SimpleFeatureBuilder;
import org.geotools.geometry.jts.JTSFactoryFinder;
import org.geotools.map.FeatureLayer;
import org.geotools.map.Layer;
import org.geotools.styling.SLD;
import org.geotools.styling.Style;
import org.opengis.feature.simple.SimpleFeature;
import org.opengis.feature.simple.SimpleFeatureType;

public class DrawLine {
    public Layer getLayerLineByCoord(Coordinate[] coords) {
        GeometryFactory geometryFactory = JTSFactoryFinder.getGeometryFactory();
        LineString line = geometryFactory.createLineString(coords);

        SimpleFeatureBuilder featureBuilder = new SimpleFeatureBuilder((SimpleFeatureType) TYPE);
        featureBuilder.add(line);
        SimpleFeature feature = featureBuilder.buildFeature("LineString_Sample");

        DefaultFeatureCollection lineCollection = new DefaultFeatureCollection();
        lineCollection.add(feature);

        Style style = SLD.createLineStyle(Color.BLUE, 1);
        return new FeatureLayer(lineCollection, style);
    }
}


texts:
I want to push/write some generated raster layers to my postgres database directly in R. In theory this could be supported by the rgdal package, but sadly it isn't. In the latest version only copy, but not write is supported (Can be queried by gdalDrivers() ).

Earlier this year it was announced that there might be a coding project as part of the google summer of code competition
Link

Does anyone know what has happened there and if write via writeGDAL is now possible?


texts: How to reset the length of the control of scaleLine in OpenLayers?

texts:
I want to change the style of the control of the ScaleLine,but I have no idea where to change it,anybody can give me some advice?

map.addControl(new OpenLayers.Control.ScaleLine());


texts: Modifying ArcGIS Network Analyst route when route created is flawed?

texts:
I'm analyzing an optimized routes for garbage collection. There are roughly 250 orders per route. I've ran the Network Analyst tool both with time and length impedances. In both cases, the output solution will indicate the route will 'drive' past an order without stopping. I'm looking to figure out how to correct this in my analysis but without having to closely examine each of the 1,000 pickup points.

I'm using the Route tool, partially because when I do a Fleet Analysis, the problem still exists, but GIS takes longer to solve. Impedance is set to Minutes, I've checked Reorder Stops to Find Optimal Route, as well as preserve First and Last Stop. U-Turns are allowed only at Dead Ends.

The attached graphic is the first example I could find. The route travels from 28 to 29, passing 40 in the process.


texts: Automate filling of attribute table using QGIS?

texts:
I am new to QGIS and wish to automate my current process.

In my project I have to open the attribute table of the batch_16 layer having 5 fields (date, checked_by, td_created, td_catagor and op_comment).

There are 4 additional layers also (woody, sand, builtup, baresoil)

The logic is, if I make any new features in any of these 4 layers after jumping to a grid tile of batch_16 layer as shown in the pic, then I have to assign some attributes to the selected grid tile attribute field as below:

td_catagor: what layers I have worked on after pressing "zoom to selected" to the grid tile and before moving to another selection. woody - W, baresoil - B, sand - S, builtup - BU. If there are more than one layers I have created new features then each letter should be separated by comma as depicted in the image.

How can I automate this process?


texts: WFS GetFeature request with different SRS

texts:
I'm having problems with a WFS GetFeature request. I'm using OpenLayers with an OSM map as base layers which use the Google projection. Overlaying the layer with WMS works perfect. When I do a request to return all features within the bounds of a polygon, no features are found.

My layer uses EPSG:28992.
the base layer uses: EPSG:900913

My layer configuration in GeoServer:

When I set the declared SRS to EPSG:900913 and set "SRS handling" to "Reproject native to declared" it works fine, but I need OpenLayers to reproject it, not GeoServer.

Below is an example request I would do. I think the problem is that GeoServer does not convert the spatial filter to the native SRS.

http://mydomain.com/wfs?request=GetFeature&version=1.1.0&
srsName=epsg:28992&
typeName=Namespace:Layer1&
propertyName=(*)&
CQL_FILTER=year = 2011 and INTERSECTS(geometry,POLYGON((565585.61613069 6816272.786200799,565585.61613069 6817032.3791692,564362.62367835 6817032.3791692,564362.62367835 6816272.786200799,565585.61613069 6816272.786200799)))&
outputFormat=GMl2&

Update

I just tried re-projecting some points from OpenLayers (EPSG:900913) via proj4js to EPSG:28992 (PostgreSQL and GeoServer). I took 3 points and  compared them. They all showed an offset.

Converted: POINT(233492,24349424947 582191,0681508556)

In Database: POINT(233530,852147064 582315,680920578)

Difference x: 38,60865281453
Difference y: 124,6127697224
Converted: POINT(133459,1980658749 455770,7481867363)

In Database: POINT(133486,356985081 455879,67223342)

Difference x: 27,1589192061
Difference y: 108,9240466837
Converted: POINT(176920,43262818802 317930,36182729155)

In Database: POINT(176954,154083767 318023,515324756)

Difference x: 33,72145557898
Difference y: 93,15349746445

If the offset was constant I could correct it. Now I'm not sure if the offset is caused by GeSserver or by proj4js. I'm guessing proj4js since the WMS layer does not show any offset.


texts: postgresql 8.4 and postgis 1.5 error: type modifier is not allowed for type 'geometry'

texts:
I would like to generate a buffer (within a limit of 3km  )  where the user clicked point.
In pgAdmin, select my database and open a query editor and create 2 PostGIS functions (UTM Zone Function and ST_Buffer_Meters Function)

Used following SQL command to create another column in to my table. my table is “schools” with lat and lon fields.
ALTER TABLE schools ADD COLUMN geom_buffer geometry(POLYGON,4326);

but appear error: type modifier is not allowed for type 'geometry'

Why is that..??/ Pls. help me anyone. ..


texts: Is there a way to filter a list of layers under the drop down whenever selecting a layer to run a geoprocessing tool? Because I have 40 some layers displayed in the Table of Contents of ArcMap which can be pretty messy.


texts: How to create postcode areas for country with none?

texts:
How are postcodes created?

I know postcodes were first created by the Royal Mail in the United Kingdom.

If I want to develop and create postcode areas for a new country, how do I go about it?

What is the first step?
What type of professional people do I have to engage?
What software can I use?
Where  can I get more information?

I'm planning to  work on this project.


texts: 2D/3D Area parameter differences in Surface Volume calculation dialog box?

texts:
In dialog box of volume surface calculation and also in ASCII file there are 2 out come with names 2D Area and 3D Area. What is difference between these two numbers? What is included in 3D that is not in 2D? For example :

 Plane_Height, Reference, Z_Factor, Area_2D, Area_3D, Volume
 45,00, ABOVE, 1,000000, 644929,78213426, 716708,47208592, 17433798,859008

Note: Decimal point is  a comma.


texts: Effects of heterogeneous vegetation in mapping canopy cover using randomForest in R

texts:
I have a 5.8 m resolution satellite image (LISS IV aboard IRS-P6, in the wavelength (micrometers): 0.52 - 0.59 (band 2),0.62 - 0.68 (band 3) and 0.77 - 0.86 (band4). The swath covers a densely forested region in the Eastern Himalayas. My objective is to predict regions of canopy cover. Regions of canopy vs. non-canopy.

I intend on using several predictor variables in light of: Incorporating terrain data to predict canopy cover using randomForest in R
and The stability of randomForest models after increasing predictor variables.

The area covered in my image ranges from 700 m to 3500 m and comprises at-least 3 forest types excluding high altitude grasslands and other landcover (water, rock etc.). As I want to map canopy cover over this entire region, what would be the ideal set of training data to use?

Should it cover all spectral classes of vegetation? As my objective is just canopy I suspect the training points should constitute a wide range of tree crown types.

Edit: Wish to map canopy vs. non-canopy within classes of canopy coverage (canopy cover percentage per pixel).

Definitions: by canopy I specifically mean the crown area of a primary stand. Secondary stands with a crown are not the target.


texts: Cesium Polygon entity get all positions and holes

texts:
I could not find a nice way to get the all positions and holes from Cesium polygon (see the coordinates structure in attached image), and as far as possible not to use recursive function to get all the coordinates. Also checked API and there is nothing that I can use, ex. "getCoordiantes" or something like that.

Usage to check all the coordinates height for comparing.

Here is also the GeoJSON example:

{
    "type": "FeatureCollection",
    "features": [{
        "type": "Feature",
        "properties": {
            "dataType": "AirspaceType",
            "geometryExtension": {
                "extruded": {
                    "maxHeight": 10000,
                    "minHeight": 0
                }
            },
            "id": "d50b2df0-9447-4a27-8710-eb62dff66ade",
            "name": "Unnamed Layer",
            "category": "default"
        },
        "geometry": {
            "type": "Polygon",
            "coordinates": [
                [
                    [0, 0, 0],
                    [0.2, 0, 0],
                    [0.2, 0.2, 0],
                    [0, 0.2, 0],
                    [0, 0, 0]
                ],
                [
                    [0.03504, 0.096817],
                    [0.035384, 0.18917],
                    [0.18132, 0.18917],
                    [0.174109, 0.020599],
                    [0.03504, 0.096817]
                ],
                [
                    [0.059077, 0.100937],
                    [0.059077, 0.173034],
                    [0.157283, 0.173034],
                    [0.157283, 0.100937],
                    [0.059077, 0.100937]
                ],
                [
                    [0.07968, 0.154838],
                    [0.13153, 0.156555],
                    [0.136681, 0.124626],
                    [0.094102, 0.118446],
                    [0.07968, 0.154838]
                ]
            ]
        }
    }]
}


texts: Adding a Pegman of street view with a Leaflet map

texts: I want to add a Pegman of street view to my Leaflet map. How do I add it for free and open source way to have the possibility to see street view of every place in the map?


texts: All feature classes are not shown when creating new topology in ArcCatalog?

texts:
I am generally a self taught GIS guy and I typically figure things out with extensive research. I have been reading through the ArcGIS Help online, specifically on Creating A Topology. The link provides a step by step process on creating the a new topology.
I have a variety of layers related to a sewer system and I would like to create a topology so that I can ensure everything is connected appropriately prior to setting up a geometric network. The process I have followed in ArcCatalog is as follows:

Right click the relationship class in the geodatabase, select New, then select Topology.
Click next, in the dialogue that pops up I set the name of the topology and the cluster tolerance. Click next.
On the next screen I have to select the feature classes that will participate in the topology. Unfortunately all the feature classes listed and selectable are not all those that are in the relationship class of the geodatabase. See the following screenshot:

Turns out that only one (1) of the features classes that would be part of the new topology is selectable. I need many of the others that are not listed and selectable in the dialogue box.


texts: Calculate XY Coordinate with "CalculateField" tool

texts:
I used model builder to build a model which can calculate the XY coordinates
and it's perfect!

However, when I export it to VBscript and run it in VBA it crashes....

The error happens at

 gp.CalculateField_management RasterT_Kriging1_shp__2_, "x", "pPoint.X", "VB", "Dim Output As Double"

and says object required.

Following is my code to calculate the X:

Dim Output As Double

Dim pPoint As IPoint

Set pPoint = [Shape]

Output = pPoint.X

Can anyone help me?

Ok, following is my code:

Dim gp As Object
Dim shape As Object
Set gp = CreateObject("esriGeoprocessing.GPDispatch.1")

' Check out any necessary licenses
gp.CheckOutExtension "spatial"

' Load required toolboxes...
gp.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Spatial Analyst Tools.tbx"
gp.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Conversion Tools.tbx"
gp.AddToolbox "C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management Tools.tbx"

' Local variables...
Temp = "C:\temp"
Kriging_poin4 = "C:\temp\New File Geodatabase.gdb\Kriging_poin4"
Output_variance_of_prediction_raster = ""
Dim Temp As String
Dim Kriging_poin1 As String
Dim Temp__2_ As String
Dim RasterT_Kriging1_shp__4_ As String
Dim RasterT_Kriging1_shp As String
Dim RasterT_Kriging1_shp__2_ As String
Dim RasterT_Kriging1_shp__6_ As String
Dim RasterT_Kriging1_shp__5_ As String
Dim point2_Events_shp As String
Temp = "C:\Temp"
Kriging_poin1 = "C:\Temp\Kriging_poin"
Output_variance_of_prediction_raster = ""
Temp__2_ = "C:\Temp"
RasterT_Kriging1_shp__4_ = "C:\Temp\RasterT_Kriging.shp"
RasterT_Kriging1_shp = "C:\Temp\RasterT_Kriging.shp"
RasterT_Kriging1_shp__2_ = "C:\Temp\RasterT_Kriging.shp"
RasterT_Kriging1_shp__6_ = "C:\Temp\RasterT_Kriging.shp"
RasterT_Kriging1_shp__5_ = "C:\Temp\RasterT_Kriging.shp"
point2_Events_shp = "C:\Temp\point2 Events.shp"

' Process: Feature Class To Shapefile (multiple)...
gp.FeatureClassToShapefile_conversion "'C:\temp\New File Geodatabase.gdb\point2'", Temp__2_

' Process: Kriging...
 gp.Kriging_sa point2_Events_shp, "x", Kriging_poin1, "Spherical 249.101224", "250",   "VARIABLE 12", Output_variance_of_prediction_raster

' Process: Raster to Point...
gp.RasterToPoint_conversion Kriging_poin1, RasterT_Kriging1_shp__4_, "Value"

' Process: Add Field...
 gp.AddField_management RasterT_Kriging1_shp__4_, "x", "DOUBLE", "12", "4", "", "", "NON_NULLABLE", "NON_REQUIRED", ""

' Process: Add Field (2)...
gp.AddField_management RasterT_Kriging1_shp, "y", "DOUBLE", "12", "4", "", "",   "NON_NULLABLE", "NON_REQUIRED", ""

' Process: Calculate Field...
gp.CalculateField_management RasterT_Kriging1_shp__2_, "x", "pPoint.X", "VB", "Dim Output As Double"

Dim pPoint As IPoint

Set pPoint = [shape]

Output = pPoint.X

' Process: Calculate Field (2)...
gp.CalculateField_management RasterT_Kriging1_shp__6_, "y", "pPoint.Y", "VB", "Dim    Output As Double"

Output = pPoint.Y

End Sub

I am following the steps you told me last time in

Using kriging with data in GDB

but it's still doesn't work!


texts:
I tried to import a building-shapefile (ax_gebaeude) into the postgis database (Postgresql 11). I normally use Qgis 3.12 with drag & drop for that. In this case I get an error message:
Failed to import some layers!

Creation of fields failed

The layer can be openend by QGis or ArcGisPro without any problems. I also tried to use the tool fixe geometries before. I renamend all columns in small letters. I tried to export the file via ArcGis, to get rid of the problems but the problem remains. Any ideas?
The field names are:


texts: How to import from E00 Interchange File with ArcGIS 10.8?

texts:
I have a dataset in old ArcInfo E00 Interchange format.  In the past I could extract the data from this by using the "Import From Interchange File" tool in the Coverage toolbox, however this toolbox doesn't appear to be present in ArcGIS 10.8.
I searched around for more info and found the ArcGIS Coverage toolbox licensing page that said I needed ArcInfo Workstation installed to use these old tools.  Unfortunately I don't have access to ArcInfo Workstation anymore.

Import From Interchange File - Requires ArcInfo Workstation installed

How does one go about extracting from E00 in more recent versions of ArcMap that don't have ArcInfo Workstation?


texts:
What is the difference between a "Feature Class" and an "SDE Feature Class"?
Or a "Table" and an "SDE Table"?
When browsing data in ArcCatalog, hosted on the same enterprise gdb, I see these data types listed and am not sure what the difference(s) between them are.

Also, in the attached ArcCatalog screenshot below are two line feature classes.  Their icons are different.  Why are different icons signifying?


texts: How to multiply two field values in pyqgis?

texts:
In a vector layer table there are two fields F1 and F2 with float type. I want to multiply these fields and insert results to new field F3.
In below code new field created but without any value. What is problem?

layer = iface.activelayer()
F1 = self.ui.cBox_Width.currentText()
F2 = self.ui.cBox_Length.currentText()

new_Field = "F3"
provider = layer.dataProvider()
caps = provider.capabilities()

idx = provider.fieldNameIndex(new_Field)
try:
    if idx == -1:
        if caps & QgsVectorDataProvider.AddAttributes:
            res = provider.addAttributes([QgsField(new_Field, QVariant.Double)])
except:
    return False

layer.startEditing()
bufferfield = layer.fieldNameIndex('buffer_km')
formula = QgsExpression("\"F1\" * \"F2\"")
formula.prepare(layer.pendingFields())
for f in layer.getFeatures():
    f[idx] = formula.evaluate(f)
    layer.updateFeature(f)

layer.updateFields()
layer.commitChanges()


texts:
I'm creating a plugin to guide a user though a GIS analysis process.  One of the things the user wants is the ability to display "charts" of the resulting data sets.  I'm looking to build simple statistical charts such as, histograms, bar charts, pie charts, etc of the fields in the resulting data set.  The charts to be displayed are not geo-specific so the built in Diagrams are not quite what I'm looking for.

I'm looking for an approach to display plots of the result data set as a layer or a using either a python (matplotlib, ggpot for python) or javascript (d3, vega (vega-lite), plotly, etc) plotting package.


texts: QGIS select feature, "add to current selection" feature

texts:
I am trying to perform a feature selection on qgis similar to "add to current selection"  interactive selection method in Arcgis Desktop, is that possible?

An example: I perform a spatial query and select a set of features, afterwards I need to manually add a couple of features to this set,  if I use the feature selection button I lose my first selection, and start a new one.

Any way of adding features to a selection? Many thanks everyone in advance,


texts: Maximum zoom levels for tilestache rendering mbtiles

texts:
I have been given an mbtiles file which contains tiles for zoom levels between 16 and 20. When I preview the mbtiles on tilestache, I am not able to zoom in beyond the 18th level. If I request the map at zoom level 19, tiles of zoom level 18 are being rendered. I have viewed the mbtiles files and sure enough there is text "minzoom16" and "maxzoom20".

There is another mbtiles file with tiles from zoom levels 11 to 16. It is rendering tiles requested for zoom level 18 too; beyond zoom level that no tiles.

This makes me wonder if there is any internal setting in the tilestache code that limits the server to rendering tiles up to a certain zoom level. I have searched the source, but couldn't find anything in particular; not that I am a python programmer, but still.

So, does tilestache limit the zoom levels to a certain number? If yes, how can it be changed?

Thank you.

EDIT:

"bholakpur":
{
   "provider": {"name": "mbtiles", "tileset": "/home/***/****/tiles/bholakpur.mbtiles"},
   "preview": {"zoom": 20, "ext": "png"}
}

EDIT 2:
Bing layer is the base layer for the map and as far as my tests go Bing has only 18 zoom levels. May be that is the reason the other layers are not requesting tiles beyond the 18th zoom level.


texts: Creating shapefiles using Java and Accumulo

texts: I have tracking data with latitude and longitude stored in an Accumulo database table.  We want to move toward using Geoserver for some layer display features on the client side.  Right now I'm looking through the Python script project that will generate shape files.  I'm wondering if anyone has other suggestions on ways I can get down to the basic functions needed to take Lat/Long trails and convert them to shapefiles for geoserver to pick up.


texts: ArcGIS: Maintaining the same layer styles over multiple data frames

texts: Let's say you have a map with three or four data frames. In each data frame there is the same road layer. Is there a way to ensure that all the symbology/labeling/annotation is the same for that layer across all data frames. So let's say I modify the symbology of the road layer in data frame 1, is it possible to automate that change across all data frames?


texts: Rename long column names for dbf conversion

texts:
Is there an ArcObjects API call that will rename long column name? If a column name is too long for conversion to a .dbf file, I want rename to a valid length.

I have searched ArcObjects help. Either such a function does not exist, or my search terms are not getting it done.

I did find a field validation function.

If no such function exists, is there a algorithm or renaming standard I could reference to roll my own?


texts: Required Windows installer for "Geoserver 2.15.2"

texts:
Checking latest version of GeoServer 2.15.2 & found that Windows installer is not available

However "war" file option is available ,also downloaded "Platform Independent Binary" but it looks like pure Windows installer is not available as of now OR changed their deployment strategy for Windows not sure.

Does anyone know?


texts: Image collection filtering doesn't work in JavaScript client library Google Earth Engine

texts:
I am trying to create a beautiful and intuitive interface for working with Google Earth Engine. In the process, I needed to filter the collection of satellite images. I get an error when calling any filter function.
What could be the mistake?
TypeError: Cannot use 'in' operator to search for 'Object' in undefined
at Function.ee.promote_ [as promoter_] (browser.js?9dfd:23586)
at ee.ApiFunction.ee.Function.promoteArgs (browser.js?9dfd:21333)
at ee.ApiFunction.ee.Function.apply (browser.js?9dfd:21310)
at ee.ApiFunction.ee.Function.call (browser.js?9dfd:21307)
at Function.ee.ApiFunction._call (browser.js?9dfd:21405)
at Function.ee.Filter.lt (browser.js?9dfd:21979)

Code example:
const defaultShot = new ee.ImageCollection('COPERNICUS/S2_SR')
const filterShots = defaultShot.filter(ee.Filter.lt('CLOUDY_PIXEL_PERCENTAGE', 20))

I am writing a SPA application in Vue CLI 3.


texts:
I am editing road network to create scenarios but Shortest Path is not recognizing the edits. I am assuming the edited street network is not connecting with other streets at intersections. I have extended yellow line to meet the red line but Shortest path is not recognizing the connection.
Could anyone please tell me how to solve this issue please?


texts:
I just try to visualize some points using GeoServer to publish them as WMS and bind this WMS to my Openlayers website using the GeoServer internal Heatmap process to visualize it. I don't want to use the Openlayers Heatmap function.

But there seem to be some problems, look here:

I'm using this Openlayers implementation to bind my WMS URL:

https://openlayers.org/en/latest/examples/wms-tiled.html

Changing the params to:

params: {'LAYERS': 'test:points','TILED': true,'STYLES': 'heatmap'},

And using this as SLD to define the Heatmap:

http://docs.geoserver.org/stable/en/user/styling/sld/extensions/rendering-transform.html

Any ideas why this happens?


texts: Is it possible to get the style of a tab file I open in QGIS?

texts:
I'm trying to duplicate a map made in mapinfo pro with QGIS in order to create shp files for geoserver.
When opening the tabs in QGIS they get some random style instead the one defined in mapinfo pro.
Is it possible to auto create the same styles or do I need to manually redefine them


texts:
I am trying to add Python script as a tool in ArcMap 10.6. I have the result or outcop3 saved in the folder Outputs. I added a output folder in the parameter and made the direction as output but i am unable to find this Outputs folder when i click on the parent folder. All i can see is the main folder called "Test". Why is it so? I am adding the script i used here.
import arcpy
import os
import sys
from arcpy import env
from arcpy.sa import *
arcpy.CheckOutExtension("Spatial")
arcpy.env.overwriteOutput = True

Raster1 = arcpy.GetParameterAsText(0)
......
Outputs = arcpy.GetParameterAsText(4)
env.workspace = r"C:\Test"
arcpy.env.cellSize = "MINOF"
Raster1 = Raster("Raster_one")
.....
outcop3 = Con(Raster4 == 1,2)

outcop3.save(r"C:/Test/Outputs/outcop3")
´´´


texts: Can privacy be achieved when using Bing Maps API for local authority use case?

texts:
We have recently won a contract to provide a website for project run for a local authority. As part of the website, there should be a page with a simple map showing the project's locations, when you click one you get info text and an illustration photo.
We have implemented it using Bing maps, but had to take it off because the local authority is "concerned that the use of Bing Maps is breaching Privacy".

I couldn't find any similar website that uses Bing Maps, but still can't' understand their concern.

Many websites that belong to the authorities are using Google Maps which I'm sure has the same Privacy policy as Bing. Example is the journey planner in BVG website, and Germany was the one to force Google to create the anonymizeIP function in Analytics.
I seek any input about using Bing Maps in a local authority project website.


texts:
How to convert .geotiff files to .tiff files (ordinary plain images)?

Basically, I have cropped the image to a rectangle of interest and
knowing its left-lower location and having chosen the pixel size,
I could combine several data sources as ordinary .tiffs (and could
finally forget the tiresome ArcGis world..)

Any working method goes. We will employ usual image processing tools,
and the location information is useless at that phase.

Thank you for any help! =)

UPDATE:

This is what I got when I tried to read GeoTiff in GIMP 2.8.4:


texts: Fetching Tiles for Zoom Level from Data Base

texts: I have 6.7 million tiles of 200 x 200 meter each stored in a remote database, Can i plot them on the web dynamically as my zoom level changes or i move towards a new area my DB is queried and files matching this zoom level are fetched, Using GeoServer.


texts: Why are a polygon's vertices jumbled after importing KML into MySQL?


texts: Migrating web app from ArcGIS Online to ArcGIS Server?

texts:
I have a webapp built in ArcGIS Online using the available templates, I want to move it to ArcGIS Server, I cannot find Esri's 'official' way to do it.

Where can I find it?