texts: Select Lines that have points on edge but not on end

texts:
I am looking for an easier way (if it exists) of selecting all lines that have points snapped to them but where the points are not on the end-points of the line.

My current method is to:

Generate temporary end-points for all lines (Feature Vertices to Points using BOTH_ENDS option
Select all points (from actual point layer) that intersect my lines
Remove from selection any points that intersect my temporary end-point features
Select by location any lines that intersect the remaining selected points

This seems quite long-winded to me, and as I'm having to do this process regularly I'm keen to know if there is a simpler way to do this.  I will be creating a model or arcpy tool to do the above for me, but would prefer to have a quicker process to follow.

I have 500+ points and 15,000+ lines, so my current process takes a while to run each time.

EDIT: My data is in a geometric network, and EVERY line end-point has a point attached as required by the network.  In this situation I am trying to find all lines with specific types of points that are snapped to the lines but are not on the ends of the lines.  These are not errors, so won't be flagged by network rules.


texts: Latitude - Longitude bounds from GeoTIFF info

texts:
here's the content of a GSFC-generated GOES TIFF image:

GOES 13 Imager frame  66 at UTC 12:45:53.520 day 312 of 2013
Vis pixels :   9061, 22912    lines :   2845, 10152
Vis Lon    : 999999,-23.1566 Lat   : 999999,-20.9474
Imager ch3:count(0,1023) => [(modea-0)/1]^1 =Uchar> (255,69)
(xscale,yscale) => (0.190476,0.333333)
Vis pixels:  3650,  28850 lines:  2250,  14250
(xstride,ystride) => (21,12)
IMC status:      1
nadir pixel,line:  15352,  8050
Ref lon,d,lat,yaw:  -74.99990305,   42164.17478,         0,         0
Attitude roll,pitch,yaw:          0,         0,         010491171

The TIFF file is here: http://goes.gsfc.nasa.gov/goeseast-lzw/overview/ir3/latest.tif

¿How can I get/calculate the geographic bounds from this information or from the GeoTIFF file itself?

¿Can I do this somehow opening the file in QGIS?  Do I have to develop some script?

I only need to do this one time.

Thanks to anybody dedicating his/her time on answering this!


texts: How to add automatically a color for a query made in Access

texts:
I made an Access database with my polygons and my Exceltable with data about floods. I linked the specific row to my shapefile table and made a query table of it.

With the manual way I gave the shapefile with "Convert symbology to Representation" automatically the color red (red is linked to a colorstyle based on the amount 18 floods). With this function the shapefile will be opened with the right (red) color and not with for example blue.

I would like to add a function like this to the query shapefile. The attribute table of this shapefile automatically changes when I add a flood to it in Access, but this shapefile still opens with the blue color instead of the red one (the Convert symbology to representation tool doesn't work for the query). When I already changed the color to red and add a new flood, the colorstyle also stay at 18, but the total amount (where the color is based on) is already 19.

I there a way to solve this problem?


texts: QgsPointDisplacementRenderer: setLegendSymbolItem creates dump

texts:
trying to set up a nice layout using the QgsPointDisplacementRenderer,
but I get a dump everytime I try to load the layer I create through the following renderer.. The dump only occurs when I use the setLegendSymbolItem lineidea anyone?

renderer=QgsPointDisplacementRenderer()
renderer.setCenterSymbol(svgsymbol)
renderer.setLegendSymbolItem("",markersymbol)    # RuleKey as String?!?
layer.setRendererV2(renderer)
layer.triggerRepaint()

Looks like the problem is the 'rule key', which Python wants as a string, but which is an integer, when I request it:

lsi = layer.rendererV2().legendSymbolItemsV2()
for l in lsi:
    print l.ruleKey()

so, until now,

QgsMapLayerRegistry.instance().addMapLayer(layer)

creates a dump.

The symbols were defined as follows:

svgStyle = {}
svgStyle['name'] = r"https://upload.wikimedia.org/wikipedia/commons/1/10/Drama_Masks.svg"
svgStyle['size'] = '8'
svgsymbollayer = QgsSvgMarkerSymbolLayerV2.create(svgStyle)
svgsymbol = QgsSymbolV2.defaultSymbol(layer.geometryType())
svgsymbol.changeSymbolLayer(0, svgsymbollayer)

markerprops = {'name': 'round', 'color': 'white'}
markersymbol=QgsMarkerSymbolV2.createSimple(markerprops)


texts: Google Earth Engine Greenest Pixel Mosaic Landsat 5 - Error: Image.select: Pattern 'NDVI' did not match any bands

texts:
Using the standard Earth Engine Quality Mosaic function I have successfully compiled Greenest Pixel mosaic for multiple years using Landsat 8 and Landsat 5. However recently I encountered an issue with Landsat 5 images prior to 2003 when the task is run (no errors in the console).
In the Run output I get the error:

"Error: Image.select: Pattern 'NDVI' did not match any bands"

The source script has worked perfectly for other years (including Landsat 8 TOA). Is there an issue with using the Greenest Pixel quality mosaic for pre-2003 Landsat 5 images?

Map.addLayer(fc);
Map.addLayer(bounding);
var start = ee.Date('2002-07-01');
var finish = ee.Date('2003-06-30');
var L8_Mackay_2002_2003_Greenest = ee.ImageCollection('LANDSAT/LT05/C01/T1_TOA')
.filterBounds(fc)
.filterDate(start,finish);
var count = L8_Mackay_2002_2003_Greenest.size();

var cloudlessNDVI = L8_Mackay_2002_2003_Greenest.map(function(image) {
  // Get a cloud score in [0, 100].
  var cloud = ee.Algorithms.Landsat.simpleCloudScore(image).select('cloud');

  // Create a mask of cloudy pixels from an arbitrary threshold.
  var mask = cloud.lte(20);

  // Compute NDVI.
  var ndvi = image.normalizedDifference(['B4', 'B3']).rename('NDVI');

  // Return the masked image with an NDVI band.
  return image.addBands(ndvi).updateMask(mask);
});
var greenest_Mackay_2002_2003 = cloudlessNDVI.qualityMosaic('NDVI');
var greenest_Mackay_2002_2003_image = ee.Image(greenest_Mackay_2002_2003)
.select (['NDVI']);

// Export the image, specifying scale and region.
Export.image.toDrive({
  image: greenest_Mackay_2002_2003_image,
  description: 'greenest_Mackay_2002_2003_image',
  scale: 30,
  region: bounding
});

Greenest Pixel Mosaic Earth Engine Code


texts: Select GeoServer vector tiles feature using OpenLayers and show features attribute in popup?

texts: I have installed vector tiles extension in GeoServer and published a PostGIS tables using this format. The requested vector tiles from GeoServer using Openlayers works fine. Now, I want to select and highlight a feature of a polygon layer. When I select a feature, it will be highlighted by color and then, open a popup in which show the feature attributes information. I have no idea, how can I do it?


texts: Stopping default deselect on click of map OpenLayers

texts:
I am trying to create a box selection where I don't want to unselect until I click on an action button.
Right now its getting unselected as soon as I click somewhere on points or anywhere on the map. Also, at the same time select and unselect gets propagated which I tried to stop using option (remove condition, but ideally I am not adding any condition, so it won't work.)
here is my experimental code

I tried to stop event propagation on click of map after selection
but didn't find it of any use as I exactly don't know which event to
stop.
I also though of setting deselected points back to selection
as an undo option. But no method to set into select.
And lastly in select I thought of option but no such option found.


texts: x, y to decimal degrees on mouse click in ArcObjects?

texts:
I am converting an old VBA code I used to use in ArcMap 9.2, into VB.NET for ArcMap 10. It should extract  x, y coordinates from OSGB and converts to decimal degrees while copying to clipboard. I converted it and does not work. Nothing happens when I click on button.

Protected Overloads Sub OnClick(ByVal sender As System.Object, ByVal e As System.EventArgs, ByVal shift As Long, ByVal X As Long, ByVal Y As Long)

    Dim pMxDoc As IMxDocument
    Dim pPoint As IPoint
    Dim pClone As IClone
    Dim pGeometry As IGeometry
    Dim pSpatialRefFactory As ISpatialReferenceFactory
    Dim pSpatialRef As ISpatialReference
    Dim pGeographicCoordSys As IGeographicCoordinateSystem

    'Get the point where the user clicked
    pMxDoc = My.ArcMap.Application.Document
    If pMxDoc.CurrentLocation.IsEmpty Then Exit Sub
    'Clone the point because we don't want to alter
    'the actual document's current location point
    pClone = pMxDoc.CurrentLocation
    pPoint = pClone.Clone
    pGeometry = pPoint 'QI

    'Create a new geographic coordinate system to use in the conversion
    pSpatialRefFactory = New SpatialReferenceEnvironment
    pGeographicCoordSys = pSpatialRefFactory.CreateGeographicCoordinateSystem(4277)
    'esriSRGeoCS_OSGB1936
    pSpatialRef = pGeographicCoordSys 'QI
    pSpatialRef.SetFalseOriginAndUnits(-180, -90, 1000000)

    pGeometry.Project(pSpatialRef)
    MsgBox(pPoint.X & ", " & pPoint.Y, , "Decimal Degrees")

    Dim MyData As New DataObject

    MyData.SetText(pPoint.Y & ", " & pPoint.X)

    MyData.PutInClipboard()
    MsgBox("Data copied to clipboard successfuly")

End Sub

I get only 2 warnings not errors with message:

  A reference was created to embedded interop assembly 'ESRI.ArcGIS.Carto' because of an indirect reference to that assembly from assembly 'ESRI.ArcGIS.ArcMap'. Consider changing the 'Embed Interop Types' property on either assembly.   ArcMapAddin2

UPDATE:

I got the tool working by changing from button to tool.

Protected Overrides Sub OnMouseDown(ByVal arg As ESRI.ArcGIS.Desktop.AddIns.Tool.MouseEventArgs)
    MyBase.OnMouseDown(arg)

The tool works and copies to clipboard. But when I press Ok in dialog box it closes ArcMap down. but why?


texts: Processing.runalg is working in qgis python console but not producing output when executed in cmd

texts:
This is my code :

import sys

from qgis.core import *

# Initialize QGIS Application
QgsApplication.setPrefixPath("C:\\ProgramFiles\\QGISBrighton\\OSGeo4W64\\apps\\qgis", True)
app = QgsApplication([], True)
QgsApplication.initQgis()

# Add the path to Processing framework
sys.path.append('C:\\Program Files\\QGIS Brighton\\apps\\qgis\\python\\plugins')

# Import and initialize Processing framework
from processing.core.Processing import Processing
Processing.initialize()
import processing

#Convert Raster to vector

processing.runalg("modeler:vectorize","D:\\Extracted zip files\\LC_hd_global_2001.tif","15",'C:\\Users\\Administrator\\Desktop\\Raster_to_vector2BAT.shp')
print 'Raster_to_VectorBAT run successful'


texts: Calculating distance between two Lat Long points with different altitudes?

texts:
If we assume that the altitude is held constant, the great circle distance is correct. However if the radius of the circle is not constant how can the distance be calculated? In my case the altitude changes linearly from point1 to point2.

What is the correct way then, to calculate the distance between two Lat Long points with different altitudes?


texts: Erasing small polygons from bigger overlapped polygons using ArcGIS for Desktop?

texts:
I have different bigger polygons (buffers) which are overlapped (A,B,C,...) and the inner smaller ones (a,b,c,...). I want to erase smaller ones from bigger ones in the way the results be A-a, B-b, C-c, ... . If I use erase it erases smaller ones from all of overlapped bigger ones. Do you know how I can fix it?

My software is ArcMap 10.


texts: Spatial filtering technique for Image enhancement

texts:
In Spatial filtering, the output image has fewer rows and columns than the input
image because it has an unfiltered margin corresponding to the
top and bottom rows and the left and right columns of the input
matrix that the filter window cannot reach.
For example: 2 rows (top row and bottom row) and 2 columns (left-most column
and right-most column of input data) are lost during the filtering process for a 3 x 3 window.
Is there any specific formula or technique through which we can easily find out the number of rows and columns lost, like what If someone wants to find the number of rows and columns lost in 8 x 8 window of an input image or m x n window of an input image?


texts: Using the ESRI javasript API, is there a way to set the distance the map will move when panning?  When I fire the map.panUp(), for example, it pans, but about twice as far as I would like it to.


texts:
I know this has been asked before (Where can I find block (sub-district) level boundaries for India?) and it seems it is still not possible to get better than district level data as a shapefile, but are there paid services (I am thinking mapmyindia and such) that provide such data?

Has anyone tried this and how accurate will be their data?


texts: Need to iterate field values from multiple feature classes

texts:
I have a dataset with two fields: 1. Day 2. Number
I need to filter the first by the number and create feature classes based on that number as they name etc. I have successfully done this with iterate field values. But, I am trying to add a submodel (or at least create another) that will do the same thing to multiple feature classes that were created to generate feature classes by the day. So, my goal is to filter by number, create FCs for those and then filter those FCs by Day. The image below works, but stops after filtering only one of the FCs created by Number. Any ideas on how to filter these in modelbuilder?


texts: Spatial Analyst's CostPath toolbox/function has PATHCOST of 0 when using in_memory workspace

texts:
I have a point feature and am interested in finding the least cost path between two points in the feature using a raster as the cost surface. I first calculate the cost distance and cost backlink via the CostDistance tool. I then use the CostPath tool to get the actual path and then look at the PATHCOST attribute for the raster cell's with value 3.

What I do not understand is if I run the CostPath tool with output workspace pointed at a GDB, the PATHCOST comes out as expected. When I run the CostPath tool using in_memory (this is actually part of a script) as the output workspace, the PATHCOST values always come out to zero. I've debugged this until I was able to pintpoint this as culprit, happening in ArcMap as well, not just the script environment.


texts: Difference between an interaction and eventListener in OpenLayers

texts:
I am wondering when should I use an interaction instead of an eventListener.

For example I can have a handler function like this:

function handler(evt) {
   if (evt.type === 'click') {
      console.log(evt);
   }
   return true;
}

And then add it as a click listener to my map:

myMap.on('click', handler);

Or have an interaction like this:

var interaction = new ol.interaction.Interaction({
   handleEvent: handler
});

And add the interaction to my map:

myMap.addInteraction(interaction);

That gives me the exact same result and exact same ol.MapBrowserPointerEvent is logged in the console.

So what is an interaction for? (except predefined interaction like zoom etc.)


texts: I want to share offline my data on my ArcGIS Online account using Python. I found this as a result of my research.(https://desktop.arcgis.com/en/arcmap/latest/map/publish-map-services/publishing-a-map-service-with-python.htm).I get this error while doing the first step (ie CreateMap Draft example 1). But will this fix the problem for offline publishing?


texts: OL3 How to select features, cluster inside a polygon

texts:
I'm trying to check features inside a polygon:

function howinpol(evt) {
    var polygon_extent = evt.feature.getGeometry().getExtent(); // get polygone extent
    var starr_uli = "", arr_fdgdg = [];

    mySource8.forEachFeatureIntersectingExtent(polygon_extent, function(feature) {
        xl = parseFloat(feature.get('lon')); // get lon
        yl = parseFloat(feature.get('lat')); // get lat
        feature.setProperties({ visible: 'visible' });  // add prop
        feature.setProperties({ inpol: 'in' }); // add prop
        arr_fdgdg.push(feature); // add features
        // list
        starr_uli += '
' + feature.get('name_kontr') + '&nbsp;&nbsp; id:' + feature.getId() + ''+ feature.get('billing_address_city') +', '+ feature.get('billing_address_street') +'
';
    });

    $('#objects ul').replaceWith('
' + starr_uli + '
'); //name list

    $('#cntobj').text(arr_fdgdg.length); //count in list

    mySource8.clear(); //clear
    mySource8.addFeatures(arr_fdgdg); //add select features
}

but it is beyond the scope if draw a circle or freehand polygon:

I tried to wrap if(evt.feature.getGeometry().intersectsCoordinate([xl, yl])):

function howinpol(evt) {
    var polygon_extent = evt.feature.getGeometry().getExtent(); // get polygone extent
    var starr_uli = "", arr_fdgdg = [];

    mySource8.forEachFeatureIntersectingExtent(polygon_extent, function(feature) {
        xl = parseFloat(feature.get('lon')); // get lon
        yl = parseFloat(feature.get('lat')); // get lat

        if(evt.feature.getGeometry().intersectsCoordinate([xl, yl])){
            feature.setProperties({ visible: 'visible' });  // add prop
            feature.setProperties({ inpol: 'in' }); // add prop
            arr_fdgdg.push(feature); // add features
            // list
            starr_uli += '
' + feature.get('name_kontr') + '&nbsp;&nbsp; id:' + feature.getId() + ''+ feature.get('billing_address_city') +', '+ feature.get('billing_address_street') +'
';
        }
    });

    $('#objects ul').replaceWith('
' + starr_uli + '
'); //name list

    $('#cntobj').text(arr_fdgdg.length); //count in list

    mySource8.clear(); //clear
    mySource8.addFeatures(arr_fdgdg); //add select features
}

It worked, but it loops through a very long time.

have any thoughts?


texts: Calculating population size from raster layer using QGIS?

texts:
Is it possible (using QGIS) to calculate the population size from a raster layer within a certain area?

For example:

I have a population density raster layer (100x100 grid cell) for South Africa, and I want to determine what the population of a certain district is.


texts: Create new shapefile for each unique attribute using for-loop with PyQGIS?

texts:
I have a roads shapefile. In its name field, I'll create a new shapefile per attribute. I already got each unique attribute under name. The for-loop that iterates over the unique attributes and selects the features isn't working.
from qgis.core import QgsProject

def split_shp():
    # Parameters
    layer_name = 'bicolRoads_primary_secondary'
    field_name = 'name'

    layer = QgsProject.instance().mapLayersByName(layer_name)[0] # road shp
    shp_index = layer.fields().indexOf(field_name) # road shp index
    unique_values = layer.uniqueValues(shp_index) # unique attributes under 'name' field

    for i in unique_values: # for-loop that iterates over the unique attributes
        dog = layer.selectByExpression('field_name=i') # selects features matching i/current attribute
        print(dog)

        #shp_writer = QgsVectorFileWriter.writeAsVectorFormat(layer, new_shp, 'utf-8', \
        #driverName='ESRI Shapefile', onlySelected=True)

split_shp()

The result is this


texts: How to create random point in according polygon for each dataset entry in a list (with R)

texts:
I have an area divided into polygons, and a dataset of observations which has a common field with the polygon layer.

I now would like to generate a random point for each dataset entry within the according polygon.

R was already a great help in digesting the dataset and even to produce shapefiles with summed values of the dataset, so I dare to post this here in the gis.stackexchange.
It would be nice to have a solution in R, cause i have to create multiple subsets which can so easily be done there by a loop, so i do no have to filter the data manually in QGIS all the time.

Other solution are welcome as well.

Polygonlayer: 39 Polygons

Polyid
A1
O19
A13

Datalist: 1360 Data entries

Polyid, value1, value2, value3
A17, ..., ..., ...
O6, ..., ...
A3, ...
A4, ...
O6, ...
A8, ...
A17, ...

Goal: -> Having a point-shape with 1360 random points, which reside in the according polygon and have all attributes of the datalist.


texts: Logarithmic raster with a non-logaritmic legend in R

texts:
Following the example under:

library(sp)
library(automap)

loadMeuse()
demo(meuse)
# Ordinary kriging
kr = autoKrige(log(zinc)~1, meuse, meuse.grid)
spplot(kr$krige_output,'var1.pred')

I get a plot with an interpolated map and a legend. The map and the legend have logarithmic vales. In many cases, a logarithmic scaled map is what I want. However, a logarithmic legend is not very useful.

I would like to replace the legend numbers to exp(numbers). More like this(GRASS GIS):

How can I can keep the logarithmic color scaling and get a readable legend in R?


texts: QGIS 3.4. Saving Layer Status (ON/OFF) for printing different Layouts

texts:
if i have a project with a lot layers (20-30) and i want to print different map layouts is there a possibility to save different status of which layer is on or off?
So for example:

Map Layout 1:
Layer 1:ON
Layer 2:Off
Layer 3:Off

Map Layout 2:
Layer 1:Off
Layer 2:On
layer 3:On

My Problem is, when i have to change something after 2 months not working on the project and plot the map layouts again, i dont remember which layer to switch on/off in which map layout...

In Autocad you can use the "layer status save"function to save different states of layer on/off for different map layouts...

Or is there a other possible workaround for that question? (i am still a beginner)


texts:
I need to get a detailed map for my own purposes. I need it detailed as much as I can see the buildings (like google map hybrid mode). I went to http://open.mapquest.com/ but I couldn't find the way to download a map for Dubai city as example!
How can I do that?
I need the map data to apply a vector layer on it.


texts:
I recently realized that when closing out of ArcMap, the process itself remains running in the background until my computer is either restarted or I manually kill each ArcMap process individually. I could not figure out what the issue was. I did some digging and read that it could be an issue with my ArcMap template (Normal.mxt) or with an add-in that was running and holding onto some reference files.
So, I tried both options and deleted my Normal.mxt as well as any ESRI Add-Ins that I had installed. In my case the only ESRI Add-In I had installed was the Pictometry toolbar.
My issue was resolved but I am trying to determine if this is something anyone else has experience with the Pictometry add-in or was my issue was simply resolved by deleting my Normal.mxt file?


texts:
I am trying to zoom to a feature by finding its name. I have the following code and executed then I am getting error : "Cannot read property 'geometry' of undefined". Please suggest how to solve it and where I have to modify the code.

    Simple Map

      html, body, #map {
        height: 100%;
        width: 100%;
        margin: 0;
        padding: 0;
      }

</style>

<input type="text" id="county" />
<input type="button" id="search" value="Get Details" />


texts: Force Bing logo to bottom of map OpenLayers

texts: Any way to force the logo position of a Bing layer in OpenLayers?  See this image for the strange location:


texts:
I have very long curvey line which I want to smplify. I have used the simplify tool on QGIS but it is still not simplified enough.
Tgose are the parameters I have used:

The thing is that I still have way too many vertices especially on the curve area:

I have also tried to put lower andhigher numbers in the Tolerance but the results, surprisingly for me, were almost the same.

How can I reduce the vertices even more?Is there any way to do it not manually?


texts: Qgis Spatial Join: numeric values change to negative values

texts:
I run the function "Join attributes by location" in Qgis 2.2 between a point vector layer and a polygon vector layer. When I open the resulting table of contents of point layer there is a column of numbers (fields from 10 to 13) that fail, because there is the same value in every row.

In "input" a list of unique numeric values (example: 100000000001, 100000000002, 100000000003...), in output results the same value in every row (example: -12345678, -12345678, -12345678...).


texts: Selecting the closest lines to a point in QGIS

texts:
I need to select the closest (and only the closest) lines to a point in QGIS. The lines may belong to different layers. The idea is to be able to combine the selected lines to form the smallest polygon that contains the point used to select the lines. This is an example of my situation:

The yellow point is the one i should use to make the selection; the lines can be blue, red, yellow, green or black, and they may belong to different layers. I tried the NNJoin plugin, but it doesn't work for me. Also the Nearest Neighbor tool doesn't work (it returns a List Index Out of Range error).

Does anybody have a suggestion?


texts: I'm new user of this software, please can somebody guide me, how to calculate stream order in QGIS or GRASS GIS, step by step.... mean commands and how to calculate watershed in this.


texts:
I need to create a TIN using available cross sections, which I already have.  I can use Create TIN, but the issue is that it would interpolate across to a nearby cross section which could be from another stream instead of using the  the next available cross section on the same stream.

Is there a way to interpolate TIN in sequence using cross section station value?


texts: Points in the 4326 coordinate system are not projecting correctly onto the 3857 coordinate system

texts:
I have imported two layers into QGIS: a series of points in the 4326 coordinate system and a google hybrid tile background in the 3857 coordinate system.

I have On the Fly projection enabled with the project coordinate system set to 3857, but my points in the 4326 coordinate system on not projecting into the correct location on the 3857 coordinate system.

Example point coordinates:

Latitude -> 33.2139961833333
Longitude -> -86.7717682

Layer info


texts: How to get WMS layer feature attribute name

texts:
Hereby i have included my WMS layer displaying code.I referred many sites i can't able to get WMS getfeature info.Please anyone help me how to get WMS layer feature info give source code please

    WMS - GetFeature

    var osm, map, state_boundary;
    var state_code;
    var vectorLayer, vectorSource, featureRequest;

        // set up the map
        osm = new ol.layer.Tile({     // basemap layer
            title: 'OSM',
            type: 'base',
            source: new ol.source.OSM()
        });
        state_boundary = new ol.layer.Image({
            title: 'india',
            source: new ol.source.ImageWMS({
            url: 'http://localhost/geoserver/workspace/wms?',
                params: { 'LAYERS': 'workspace:IndiaState' },
                serverType: 'geoserver'
            })
        });

        map = new ol.Map({
            layers: [osm, state_boundary],
            target: 'map',
            view: new ol.View({
                center: ol.proj.transform([78.188860, 23.176633], 'EPSG:4326', 'EPSG:3857'),
                zoom: 5
            })
        });
</script>


texts: Cut polygon into new polygons using lines in ArcGIS

texts:
I'm new to GIS and I work with ArcMap.

I have a layer with an orthophoto of a graveyard. On top of that, a layer with an old scanned map of the same area showing all graves.

I would like to make a layer of all the graves.

I do this by first drawing some help lines along the borders of each grave, then I draw a polygon feature, and then use the Cut Polygon Tool to cut it along the help lines.

Is there some way to speed up this process? Like a tool that cuts the polygon automatically along the help lines. That sort of recognizes all lines that cross the selected polygon and then cuts it into smaller polygons along the lines.


texts: How to make ocean drainage map using Python

texts:
I'd like to make ocean drainage map like this:

In particular I'm interested in Russia.
A relatively explicit data source is hydrosheds.org, however from there i can get data to display only endorheic basins. As a last resort I thought to try to highlight ocean basins by the touch of the basins to the oceans.


texts: Clipping a raster produces bigger raster in QGIS

texts:
I need to cut a raster layer (size 500 MB about) using a vector mask.

Going to Raster -> Extraction -> Clipper and using a vector mask in the same coordinate system of the raster, QGIS produces a .tiff raster bigger than the original raster (about 1.5 GB) and without the two and .prj file.

I cannot explain this: is there a way to have a clipped raster smaller than the original one and with the two and .prj associated files?


texts: Importing gpx track points as vector layer using PyQGIS

texts:
I have a gpx track and want to import it to QGIS using PyQGIS.
When I add the gpx manually (add vector layer),
I select track points,

and get all the points and attributes with the new layer:

When using PyQGIS(method described here), tracks import fine, but route points seem replaced by tracks.

gpx_names = ["track_points", "tracks"]
for entry in inFolder:
    for name in gpx_names:
        layerName=((os.path.splitext(entry)[0])+' '+name)
        gpx_layer=iface.addVectorLayer(os.path.join(dataPath, entry)+"?type="+name, layerName, "gpx")

(Using "|layername=" instead of "?type=" didn't help
How can I fix this/add my gpx as route points using PyQGIS?


texts: How to compute the ls factor equation for USLE model in qgis grass r.mapcalc

texts:
To compute the LS factor for the USLE model, the following equation easily works in ArcGIS.  How can I perform the same in QGIS under GRASS r.mapcalc?

LSfator = Power (("flowaccNW") * resolution/22.1, 0.6) *
Power(Sin("slope_degnw") * 0.01745) / 0.09, 1.3


texts: Aerial photo stitching and georeferencing?

texts:
I need to georeference 53 air photos from 1981 (1:25,000 scale) and was wondering what would be the most efficient, defensible, and accurate option. For context, the air photos are following a highway taken along a valley bottom and therefore the relief would not be significant as would be the case if mountains were involved. Essentially, I want to quantify the change in surface area around this highway from 1981 photos to current day; how do I process the former images?

Typically, one would georeference each air photo using current satellite imagery. However, in my case, I feel this would be difficult because the area of interest has changed through time, and therefore I fear I will not have enough reference points from current satellite imagery to georeference each air photo from 1981. Therefore, I thought I could first stitch the 53 air photos together using control points from the adjacent air photo taken the same day; essentially creating a mosaic from 1981. Then I could georeference the entire mosaic using current satellite imagery; I think on a larger scale I could find enough reference points to accurately geofrerence the mosaic. Would this process cause too much local distortion that I would no longer be able to accurately quantify surface area? Lastly, I'm curious if orthorectifying the image would significantly improve the surface estimation? I'm not looking for complete precision, but would like it to be accurate enough for a publishable manuscript.


texts:
I want to create a raster file with only the lake data with the r.lake tool. The data should only cover the lake that is red in the picture.
The inputs required for l.lake are:

Raster file with elevations
Water level
Seed point coordinates

However, I have two problems:

When I introduce the water level of the lake (highest elevation point of the lake, 180 m) then the file created covers all the surface that is below 180 m. (As in the picture, ). But how do I do it if I only want a specific area?

What is the meaning of Seed point coordinates? Is it the maximum height of the lake?

I guess that if I get the correct answer of those questions, I will be able to create a raster file with only the lake data.


texts: Given a layer of counties and a layer of roads, summarize the total mileage of roads by road name

texts:
Given a road layer and a county layer. I need to know the mileage of roads by road name in each county.

How can I do that?


texts:
How can I select by expression from a date field in a shapefile, dates that are NULL or >= 2020-01-01?
This doesn't work, it only selects those with dates >= 2020-01-01, not the NULLs
"field" is NULL OR "field" > '2020-01-01'


texts: How to install OpenLayers plugin on QGIS 1.8?

texts:
After having installed QGIS 1.8 I am not able to find the OpenLayers Plugin. In "the old days" I know I was supposed to check "add 3rd party repositories", but that button has disappeared. I tried adding the following repositories without luck:

http://pyqgis.org/repo/contributed
..openlayers.org/api/2.12
..github.com/sourcepole/qgis-openlayers-plugin.git
..github.com/nvkelso/qgis-openlayers-plugin.git
..plugins.qgis.org/plugins/plugins.xml
..plugins.QGIS.org/plugins/plugins.xml
..github.com/openlayers/openlayers

I hope someone can help me locate the OpenLayers plugin as I desperately need it.


texts:
Taking some OSM (*.osm) data downloaded from Geofabrik I have created a Spatialite database using spatialite_osm_map. Now I have problems, because some of the data that is supposed to be in the form of polygons appears in the database as polylines, which I can't style properly. As an example, objects tagged "natural=wood" appear correctly as polygons in the DB while objects tagged "landuse=forest" appear as lines.

I think this must be something to do with what I've seen called the "multipolygon problem". At present I am loading the affected objects into QGIS, exporting them as shapefiles and then polygonising these, but there must be a better way of working. Any help on this topic would be gratefully received.

I ought to mention that what I know about Spatialite and OSM data wouldn't fill a thimble.

Added later: here is a screenshot illustrating the problem. The "polygon" selected by the info tool is in fact a polyline. Nick.


texts: QGIS measure tool does not continue to measure within point symbol

texts:
We are using QGIS more and more, however, today a colleague noticed something strange when using the measure line tool.

We have point features and different symbols of different sizes and colors are stacked on each other to provide all the information for quality control. One of the main aspects is position and therefore the distance between points and of the points with respect to our framework data (streets, lots, etc.).

When measuring distance, the rubberband tool stops measuring when entering the area of the point symbol, which at display scales of ca. 1:350, can represent around 4 m.

Does someone know how to correct this?


texts:
I have been trying to convert this raster data using QGIS's Polygonize function: https://daac.ornl.gov/VEGETATION/guides/Decadal_LULC_India.html

The converted geopackage is returning many geometries with negative areas. So much so that the data is highly inaccurate, and some categories are having negative areas with x^10 sq. meters in magnitude.

The features basically resemble the ones I found in this stack overflow answer: https://gis.stackexchange.com/a/243707/161470 , where the holes are located outside the geometry and are larger than the geometry itself.

I wanted to know if the error had to do with the quality of data or the size of data we are trying to polygonize?


texts:
My ASTER GDEM2 has a average elevation error of -10m for my area. (I testet it using 90 defferential GPS reference points and compared it with old Topo25, SRTM, and GoogleMaps).

I would like to correct it by +10m for every pixel. Any idea? Whats the reason for that "error"? It is measured at a different height reference?


texts: Creating slippy map from large shapefiles?

texts:
I am trying to create a map of the UK based on a number of open data sources (Ordnance Survey Open Data, OSM, some CSV files etc).

My end goal is to create a map of the whole of the UK that I can view in OruxMaps or some other similar slippymap style offline client.

I am really struggling to work out the best way to do it though.

I have tried importing the data into QGIS, styling it and then using QGIS Server, but I can't figure out a good way of creating the slippy map from there.

I have tried TileMill, but it can't cope with the 2.5Gb contour .shp file.

I have looked around at various other things too, but just cannot work out the best way of doing this.

Any thoughts?


texts: How to add leaflet-measure-path addon to Leaflet.Editable?

texts:
I want to see a measures, generated by leaflet-measure-path on polylines, dynamically created via Leaflet.Editable 's control.
How do I do it most correctly?


texts: Arcpy get values in the 2nd field of every table in a geodatabase

texts:
I have a geodatabase in Oracle with around 50 tables. I want to get the values in the 2nd field of every table in my geodatabase.  I will be using the values in the 2nd field as an input for a geoprocessing tool. The 2nd field name is different in every table but the values are what I need for the tool. This is the code i'm working with so far:

import arcpy
aWS = r"C:\Users\dav\AppData\Roaming\ESRI\Desktop10.3\ArcCatalog\Connection to OracleGDB.sde"
arcpy.env.workspace = aWS

arcpy.env.overwriteOutput = True

luList = arcpy.ListTables("*LK")
for lu in luList:
    print lu

    Desc = "DESCRIPTION"

    cursor = arcpy.SearchCursor(lu)
    for row in cursor:
        print(row.getValue(Desc))

       #Get the values in the 2nd field of lu and put it in the "2ndField" input

       #Run Table to Domain function
       arcpy.TableToDomain_management(lu, 2ndField, Desc, "NewGeo.gdb", 2ndFieldName, "DESCRIPTION")


texts: Why force retile when zoom change in OpenLayers 2?

texts:
I'm new to GIS and OpenLayers, and I can't understand why is it that OpenLayers 2 force the base layer to retile when zoom is changed.

I'm talking about moveTo function in Grid.js, where following code exists:

// if grid is empty or zoom has changed, we *must* re-tile
var forceReTile = !this.grid.length || zoomChanged;

It is some piece of code which was added 9 years ago in version 2.5 and stayed there until now (2.13.1), so I don't think it would be a bug, but I just don't understand why it is needed. I mean, if zoom changed and I'm going to load new tiles, why to request for current tiles in the first place, if I'm going to change them anyway?

Under which conditions (if any) would it be safe to change the previous code to:

// only if grid is empty, we *must* re-tile
var forceReTile = !this.grid.length;

in order to reduce requests in my GIS?


texts: Obtain attribute value from a specific field

texts:
I have a point shapefile with a high number of attribute fields. For a single point, there is only 1 of these attribute fields completed, the others have a "NULL" value for that point. The field name for which a value is available is furthermore defined in the field "Blockname".

So, for point 1, the value in the field "Blockname" is OVSTNP_180. In the field "OVSTNP_180" (and only in this field), there is a value listed (e.g. 741).

In order to label the points, I would like to create a field "label" where this value (741 in the example above) is calculated. Which expression do I need to use (in ArcGIS)?

I assumed that, for example in ArcGIS, I could state something like this:
Label= ["[Blockname]"]
but that does not work...


texts: NOAA vdatum GDAL Variable Paths for Linux Ubuntu

texts:
I want to do vertical datum conversions on geotiffs using NOAA vdatum
The conversion status menu shows no progress and no conversion happens.
I see that Microsoft Windows needs GDAL environmental variables set.
The Windows paths are...

C:\Program Files\GDAL

C:\Program Files\GDAL\bin

C:\Program Files\GDAL\bin\gdal\apps

C:\Program Files\GDAL\bin\gdal-data

C:\Program Files\GDAL\bin\gdal\plugins

I am running Linux Ubuntu 16.04 with GDAL installed.  Do I need to set the environmental variables for Linux?  If so, I understand how to set paths (sudo -H gedit /etc/environment).

What would the paths be for Linux?


texts: Finding best path through set of many points using ArcGIS Pro

texts:
Assume I have many points, say 100 or more, that represent a station to visit.  I want to choose one of them as a starting point, and then navigate to each and every one thereafter to take a sample of something.  Essentially looking to create a table to list out and order the stations to  visit based on the shortest path from a particular starting point that I choose.
The points are not one a road or line or any but transit sort of problem, but are merely random places in say, a forest.  In sum, looking for way to choose one to start out and then navigate (walk) to each one thereafter by a shortest distance.  Forget the ups and downs, not important.


texts: Looking for Historical wind data (direction/speed) for E. Africa or greater

texts: Does anybody know of  a source for wind data ~04 - 14 for E. Africa or a larger area. I have looked at a number of sources but have not found anything for these consecutive years.


texts: How to create a function by passing polygon datatype as a parameter

texts:
I am getting stuck while executing the function,Function is getting created but while i am passing the values getting error.

CREATE OR REPLACE FUNCTION INS_CCX_ply(P_ID INTEGER,P_POLY_ID INTEGER,P_POLY POLYGON)
RETURNS voiD AS
 $BODY$
BEGIN
INSERT INTO CCX_poly(id,poly_id,poly) values(P_ID,P_POLY_ID,P_POLY);
END;
 $BODY$
LANGUAGE plpgsql;

SELECT * FROM INS_CCX_PLY(1,11,'point((1,1),(8,6),(6,8),(10,9),(1,1))');


texts: oracle 10 connection on arcgis 10, windows 7

texts:
I have problem odbc connection on arcgis 10, oracle 10 and windows 7. odbc connection on windows, running well. but when arcgis conection via arcatalog, appear error "Test connection error because of an error in initializing provider. specific driver could not be loaded due to system error 5: acccess dedied. (oracle in OraClient10g_home1, D:\Oracle\product\10.1.0\client_1\BIN\SQORA32.DLL)."

anyone can help me. many thanks before.

regards,
widi


texts: How to use the mouse select multiple features on a web based map

texts:
Let me know if this question isn't detailed or specific enough.

I have a map on a webpage.  I have a series of tools on the webpage where I select just a single feature on the map and then perform an action on that feature. i.e. reshaping the feature, deleting the feature, displaying a popup with detailed attributes about the feature.  All these work great with just one feature.

But now I have a need to select multiple features and perform the same operation on them as a group.  Since all the features will be located in close proximity, what I'd like to do is be able to select them by using the mouse to drag a rectangle over those features.  Ideally this would give me an array of features to iterate over or perhaps pass back to the server for server side operations.

Here is what I've tried so far:

        require(["esri/toolbars/draw"], function (draw) {
            tb = new draw(map);
            tb.on("draw-end", getFeatures);
            tb.activate(draw.EXTENT);
        });

    getFeatures = function (evt) {

    require(["esri/geometry/Polygon", "esri/tasks/query"], function (polygon, query) {
        var extentPoly = new polygon(evt.geometry);
        var extentQuery = new query();
        extentQuery.geometry = extentPoly.getExtent();

        dojo.forEach(map.graphics.graphics, function (graphic) {
            if (extentQuery.contains(graphic.geometry)) {
                selectedFeatures.push(graphic.getContent());
            }
        });
    });

However, when I get to extentQuery.contains, it says that contains is not a member of extentQuery.


texts:
I saw the latest update of QField 1.8 wherein users are allowed to draw/digitize in freehand lines/polygons geometries.
I was wondering how can it be enabled in QField?
We badly need it in our fieldwork/geotagging activities of natural resources.

EDIT: I have finally found out how to enable it. It could be seen in the QField Settings. The "Consider mouse as a touchscreen device" should be disabled for the Freehand Digitizing function to work.


texts: Using Overpass Turbo to return all cemeteries (nodes, ways and relations) with their municipality

texts:
I am trying to use the Overpass Turbo to retrieve from OSM all cemeteries in a given area (say, Carinthia), as well as the municipality (or other administrative unit) each is in.

Getting all cemeteries is easy enough − I made this query (adapted from taginfo) which works fine.

What I am struggling with is to add geocoding details − I want for each cemetery to know the city (or the county, or other adminstrative unit) it is in.
I have been trying to draw from examples on the wiki (particularly this one) for this, but could not figure it out.

(Also, if I understood correctly, the is_in operator can only be applied on nodes, and cemeteries can be mapped as nodes, ways or relations)

EDIT:

Building uponh @mmd answer in Listing churches with their administrative location using Overpass, I made this

area[wikidata="Q37985"]->.boundaryarea;
(
  node["landuse"="cemetery"](area.boundaryarea);
  way["landuse"="cemetery"](area.boundaryarea);
);
  foreach(
    node(w)->.n;
    .n is_in->.a;
    area.a[name][boundary=administrative][admin_level~"^[2-9]$"] -> .a;
    out center;
    convert way ::=::,
              ::id = id(),
              is_in=a.set("{" + t[admin_level] + ":" + t[name] + "}");
    out body;
  );

This does not process relations correctly yet though.


texts:
How to add JPEG images as vector data?

My requirement is if I should able to add few people data as a vector data on different layers such as Name. Ph.No, Village name, zip along with his photo. I use Erdas Imagine , ArcGis, Global Mapper and Gogle Earth. My requirement is add this information on Google Earth.


texts: Can radiometers measure radiations in the infra-red at night?

texts:
I started reading about passive sensors and radiometers in particular.
I know that microwave radiometers work day and night and in cloudy conditions. Thermal radiometers can measure energy during the night also, which makes sense since the heat is radiated from the Earth at night.

My question is whether near-infrared energy can be captured at night by radiometers from space or onboard an aircraft, on the contrary of visible passive sensors obviously?


texts: how to export multiple layers in single click to other folder

texts:
i have more then 100 layers which need to exported to another folder for further analysis is there any way i can do that in a single click.  I tried to use package map but it says georeferencing error. i want to export then as "use same coordinate system as data frame".

the next thing is, this the topographical map data of 100 sheet which obviously have same name as transportation , river, land use etc ., i do tried to make a personal data base through the command "FeatureClass To Geodatabase (multiple) " the personal data base was created as topo 1 topo 2 but they all have same feature information one is replacing all 99 ie river 1 to river 99 have same data of river 1
now suggest me what to do


texts: Cost Path tool fails saying "No Spatial Reference Exists"? All inputs are set to UTM 19N

texts: I've already run the least cost path tool several times successfully. But, now that I changed the destination input to a different shapefile (PCS: UTM 19N just like all the others) it keeps failing saying that No Spatial Reference Exists. Spatial reference is defined for all inputs though, so why is this happening in ArcMap?


texts: I've built a web application with Esri's Portal for ArcGIS for a client and it has a view count on the page before entering the application. I'm working to collect some metrics and I'm wondering if that's the total requests the application has had (zoom, resize, pan, search and zoom to a feature) or if it's just the number of times it's been opened by users?


texts: Modifying the the Find Task in ESRI-Leaflet?

texts:
I'm trying to figure out how to modify the Find Task in ESRI-leaflet so that when the user types in a new text for search, the old selection is cleared.

Currently, when you "find" one feature (For example, search for any state name here: http://esri.github.io/esri-leaflet/examples/finding-features.html), the only way to clear the highlighting is reload the entire page. If you perform another search, the second selection is added to the first, instead of replacing the first.

I have examined the documentation here: http://esri.github.io/esri-leaflet/api-reference/tasks/find.html, but can't find any indication of how to clear highlighting. I'm also pretty new to Javascript and Leaflet, and I'm wondering if there is some Leaflet methods I should be incorporating here.


texts: "EOFError: EOF when reading a line" Using raw_input() in ArcMaps Python Console

texts:
I am writing a python script to automate some ArcMaps data processing.  I am getting hung up at raw_input().  I am using raw_input to get input from a user.  A folder will be created with the folder name based on the user input.

I don't understand why this isn't working. I've viewed tutorials, read the manual, etc... and I still can't figure out why I am getting this error for raw_input().

I've managed to be able to get user to select files and directories suing the Tkinter toolkit should I look into using this toolkit to get text from a user prompt as well?


texts: How to convert multiple CSV Files to Shapefiles/Feature Classes using ModelBuilder

texts:
I am trying to convert multiple csv files into separate shapefiles using ModelBuilder in ArcMap 10.5. I have generated the following model:

In this model I am, specifying a folder containing multiple csv files. I then use an 'Iterate Tables' iterator to go through all the csv files in that folder. Then the 'Copy Rows' tool is used to write the rows from the input csv table to a new table in a Geodatabase. I am then using the 'Make XY Event Layer' tool to convert the table into an XY Layer. To make this an actual Feature Class and not just a layer file I use the 'Copy Features' tool.

When I run the model, I get an error message when the model tries to 'Copy Features':

ERROR 000210: Cannot create output C:\Users\jkhangura\Documents\Trial\Shapefile\Test.gdb\MiamiBeach
Failed to execute (Copy Features).

Can someone help me understand where I am failing with this model?


texts: How come the distance tool between QGIS and ArcMap gives a different value, even when using the same shapefile and CRS?

texts:
I have created a fishnet grid in Arcmap with size 512*512 meters over a certain other shapefile than I want to display on the web. With this in mind I used the Web Mercator CRS (EPSG:3395) to reference the data in.

I then loaded in the same fishnet grid in QGIS, using the same coordinate reference system. Now when I use the distance tool in QGIS, the width is about 312 metres, instead of the 512 meters defined (and measured) in ArcMap.

I was in the understanding that when using a projected CRS the values are in meters, so the coordinates would stay the same in both. Especially when using the exact same CRS.


texts: PyQGIS 3 Disable save project in print layout and report manager (layout menus and gui buttons)

texts:
I would like to know how I can disable the 'save project' option in the QGIS print layout and report manager layout menus. I would also like to disable the save button on both gui's in PyQGIS 3.

For example:
The following code snippets can disable elements of the project menu of the QGIS interface, I need to do the same for the print layout and report manager sub-gui's:

iface.actionCreatePrintLayout().setEnabled(False)
iface.actionShowLayoutManager().setEnabled(False)

Yields:


texts: Number of polygons in shapefile too high

texts:
I downloaded administrative shapefiles from www.gadm.com/country (China). The Python (Basemap) code below tells me that there are 2365 polygons in it. However the .csv and .dbf files included with the shapefiles tell me that there are only 344 administrative divisions at this level (which corresponds to real geography).

Forgive my lack of experience with shapefiles, but what could explain this discrepancy? Is it possible that this would be caused by small islands being considered a different polygon? If so, are there specific techniques to get rid of them?

My worry is that this increases 8-fold the computation time of point-to-polygon spatial joints, since there are 8 times more polygons.

from mpl_toolkits.basemap import Basemap
china_shapefile="C:/Users/.../CHN_adm_shp/CHN_adm2"

m = Basemap(projection='merc',llcrnrlat=15,urcrnrlat=55,llcrnrlon=70,urcrnrlon=140,lat_ts=20,resolution='c')
m.readshapefile(china_shapefile, 'provinces', drawbounds=True)

print len(m.provinces)


texts: Changing default Data Frame Order using ArcMap

texts:
I am using ArcMap 10.1.
The background: I have an .mxd that I intend to roll out across the company as our default template for all future maps/figure production. This .mxd contains two Data Frames; the first ("Main") is larger (occupying almost the entire page) and is designed to be used to hold the majority of mapped data whilst the second ("Inset") is much smaller and is to be used as an overview map.
The problem: When choosing to import this new .mxd to an existing .mxd (using the Change Layout button on the Layout toolbar) the wizard automatically selects the smaller ("Inset") Data Frame to occupy the larger part of the page. The screenshot below shows this interface with a red box indicating the two possible locations of each Data Frame:

Whilst it is possible to manually alter the order of the Data Frames using the Move Up or Move Down buttons, how do I change the default option for this setting (i.e. so that the "Main" Data Frame always appears in the larger space - "2")?
The reasoning is that I would like the process of using this new .mxd template to be as simple as possible, this is especially important as we have a wide range of users with varying abilities.
Does anyone know how to choose/change this default setting? It does not appear to be in the Data Frame Properties > Data Frame menu where I would expect it to be.


texts: Geomesa accumulo CURD data operations using WFS geoserver

texts: I have created geomesa accumulo datastore and can query features using command line. Now i want to perform data operations using Open Geospatial Consortium's (OGC) Web Feature Service (WFS) for creating, modifying and exchanging vector format geographic information. I don't want to create proxy client or deal with thrifts for programatically operating with accumulo storage. Instead what are other techniques to insert and read using filters for geomesa accumulo storage data.


texts: Table join not displaying in createLayer function in cartodb.js

texts:
I am trying to build on the cartodb tutorial to develop a map of fibre networks in Africa.  The example works fine when my SQL select is very simple i.e.

            var layerSource = {
                    user_name: 'afterfibre',
                    type: 'cartodb',
                    sublayers: [{
                        sql: "SELECT * FROM af_fibrephase",
                        cartocss: $("#category").text()
                    }]
                }

But when I attempt a more sophisticated query with a table join, such as below, the layer no longer renders.

            var layerSource = {
                    user_name: 'afterfibre',
                    type: 'cartodb',
                    sublayers: [{
                        sql: "SELECT af_fibrephase.the_geom, af_fibrephase.live, af_organisation.name FROM af_fibrephase, af_organisation WHERE af_fibrephase.operator_id = af_organisation.organisation_id",
                        cartocss: $("#category").text()
                    }]
                }

The SQL query tests ok via the API
https://afterfibre.cartodb.com/api/v2/sql?q=SELECT%20af_fibrephase.the_geom,%20af_fibrephase.live,%20af_organisation.name%20FROM%20af_fibrephase,%20af_organisation%20WHERE%20af_fibrephase.operator_id%20=%20af_organisation.organisation_id

I sense I am doing something obvious wrong but can't work out what.


texts:
I am trying to produce an earthquake map of the south-west pacific ocean area around New Zealand and its surrounding islands.

In order to produce a map which crosses the antimeridian, I took two ASTER dtm products and reprojected them from Pseudo Mercator EPSG:3857 to PDC Mercator EPSG:3832 which is centred on the pacific ocean

I added my earthquake data in csv format as points in EPSG:3857 and used 'save-as' to save as shapefile and then 'save-as' again to reproject them to EPSG:3832, however they still plot in exactly the same place, completely off from my reprojected dtm and at the wrong scale

The csv format looks like this (in qgis)

and they plot perfectly when i'm not trying to cross the antimeridian

Hopefully someone has an idea of what's going on and how to solve it..

I have also tried to use a layer extent to help when reprojecting the shapefiles but all this seems to do is cause me to lose all the data that is stored in the attribute tables!


texts: What projection system is used in China which is aligned with WGS84

texts: I have some GPS tracks in Beijing. I want to project it and convert the lat long degree values in meters for calculation in euclidian space. I was wondering what is the most suitable coordinate and projection system  China uses that closely resembles WGS 84 (should be one that is geocentric, I guess)?


texts: Placing Label in QGIS with EasyCustomLabeling?

texts:
I've converted the labels for a feature class I have using the EasyCustomLabeling tool in QGIS 3.2. Other labels within the layer are not having the same problem. I'm trying to place the label in an area that has no other overlapping labels or features. However, the label will get close but not stay quite where I want it, instead snapping far either to the right or the left. Does anyone have insight into why this is happening, or what I can do to remedy it?

Pic 1 and 2 shows the two places it loves to snap, and pic 3 is the location where I want it to be but it's phobic of.


texts:
I have a feature class that contains polygon features with a unique ID column. This shapefile is to be joined with a CSV table that I've exported to a file geodatabase table. The unique identifier in the table is the tuple (IDX, UID), there exists multiple columns with the same UID.

I would like to join the polygon feature and the table based on the relation ID=UID using ArcPy. That relation however is a one-to-one relation and I would like it to be a one-to-many relation so that one record in the polygon feature class table is joined with multiple records in the table.

Any ideas on how to do this using ArcPy?


texts: ArcMap Dynamic Text: Is it possible another layout text element in the same layout?

texts:
Is it possible to use ArcMap dynamic text in one layout text element that references another layout text element? Essentially I want to just type something in one text element, and then that text to be automatically copied into the other text element.

So far I have come up with this tag:

<dyn "type="document" name="VERSION"/>

where "VERSION" is the name of the text element I want to reference. I have also tried:

<dyn type="layout" name="VERSION"/>

Both outputs are "[empty]".

The closest ESRI help page on the topic I could find was this one:
http://pro.arcgis.com/en/pro-app/help/layouts/add-and-modify-dynamic-text.htm


texts: ImageCollection contains Images with different CRS

texts:
I'm using Google Earth Engine Python API in Jupyter. I'm working with an ImageCollection from Sentinel 2 Surface Reflectance dataset for a given area on a single day and filtering  by a region of interest (roi) polygon.

The ImageCollection only contains 2 images and one of these appears to be in a different CRS (EPSG: 32631) the other being in the expected EPSG: 4326. This is causing problems for further processing. The code for creating the image collection is:

roi = ee.Geometry.Polygon(
[[[3.94, 51.6],
    [5.3, 51.6],
    [5.3, 52.15],
    [3.94, 52.15],
    [3.94, 51.6]]], None, False)

dataset = ee.ImageCollection("COPERNICUS/S2_SR")\
    .filterDate('2020-04-20', '2020-04-21')\
    .filterBounds(roi)

Here is an screenshot of the collection in Jupyter:

I have two questions:
1. How can I ensure all images in a collection are in the same CRS (and what is the most robust method for achieving this)?
2. Why are images from one dataset in different CRSs?


texts: Implementing Access Token in Tile Server?

texts:
I am using Tileserver-GL for rendering vector tiles to web and mobile. For web, I am using Mapbox JS API GL and for Mobile, I am using Mapbox android SDK. Currently, I am using flat URL given by Tileserver-GL like

http://1.1.1.1:8080/data/v3/{z}/{x}/{y}.pbf

or for mobile use

mapView.setStyleUrl("http://1.1.1.1:8080/styles/bright/style.json");

but as you can see this URL is open and anyone can access it.

Since we are not using Mapbox Access Token, how can we protect our Url from public to abuse?


texts: How to change the PostGIS schema in raster2pgsql?

texts:
I was sucessfull when I used this mode to store Raster in Postgis (raster2pgsql.exe):

C: \ Program Files \ PostgreSQL \ 9.4 \ bin>
Raster2pgsql.exe -s 32373 -I -C -M
C: \ Users \ Geo \ Raster \ Stockpile.jp2
-F -t 100x100 Stockpile_MRM |
Psql -U postgres -d GEOSURVEY -h 192.168.1.14 -p 5432

But it goes directly to the public schema.
I used this mode to change the schema (schema_name.table), but it does not work. I saw in - http://postgis.net/docs/manual-2.2/using_raster_dataman.html, but nothing changes.

C: \ Program Files \ PostgreSQL \ 9.4 \ bin>
Raster2pgsql.exe -s 32373 -I -C -M
C: \ Users \ Geo \ Raster \ Stockpile.jp2
-F -t 100x100 BD_Maputo.Stockpile |
Psql -U postgres -d GEOSURVEY -h 192.168.1.14 -p 5432

what i´m doing wrong here?


texts: Replace a character in the GetFeatureInfo response

texts:
I'm having a lot of floating point raster datasets in GeoServer. GeoServer automatically uses a dot (.) to separate the decimals of a value.

Is it possible to change this into a comma (,)? Either within GeoServer, or by modifying the GetFeatureInfo response so it replaces . with ,


texts:
I am basic user of GIS. I have a project in Rwanda and have received CAD information in a coordinate system I am unaware of. In addition this isn't the first time I have come across similar references for similar areas. Can anyone work out the coordinate system from these coordinates:

decimal degrees -2.2710, 30.2612
unknown system 529041.54, 4748727.20


texts:
Error while building OSMBright

when i run the make.py

the following error is coming

Traceback (most recent call last):
  File "./make.py", line 123, in <module>
    build()
  File "./make.py", line 41, in build
    template = loads(templatefile.read())
  File "/usr/lib64/python2.7/json/__init__.py", line 338, in loads
    return _default_decoder.decode(s)
  File "/usr/lib64/python2.7/json/decoder.py", line 365, in decode
    obj, end = self.raw_decode(s, idx=_w(s, 0).end())
  File "/usr/lib64/python2.7/json/decoder.py", line 381, in raw_decode
    obj, end = self.scan_once(s, idx)
ValueError: Expecting , delimiter: line 9 column 9 (char 304)

What i have done wrong


texts: How to display partially a WFS layer? (just the content within a box)

texts:
I need to know how to display a WFS layers, but not the whole content, only the part within a rectangle.

I have this code:

var SMap = new OpenLayers.Layer.Vector("Status map", {
    strategies: [new OpenLayers.Strategy.BBOX()],
    protocol: new OpenLayers.Protocol.WFS({
        url: "http://demo.opengeo.org/geoserver/wfs",
        featureType: "states"
    }),
});
SMap.setVisibility(true);
map.addLayers([SMap]);

And I also have two coordinates from this code:

controlS = new OpenLayers.Control();
OpenLayers.Util.extend(controlS, {
    draw: function () {
        this.box = new OpenLayers.Handler.Box(controlS,
            {"done": this.notice},
            {keyMask: OpenLayers.Handler.MOD_SHIFT});
            this.box.activate();
    },
    notice: function (bounds) {
        var ll = map.getLonLatFromPixel(new OpenLayers.Pixel(bounds.left, bounds.bottom));
        var ur = map.getLonLatFromPixel(new OpenLayers.Pixel(bounds.right, bounds.top));
        alert(ll.lon.toFixed(4) + ", " +
              ll.lat.toFixed(4) + ", " +
              ur.lon.toFixed(4) + ", " +
              ur.lat.toFixed(4));
    }
});
map.addControl(controlS);

What can I do to display only the WFS layer's content within the imaginary box formed by 2 pairs of coordinates?

Thanks and regards,
Rafael


texts: How to flip, or invert attribute tables with respect to row ID arcgis

texts:
I need to invert the fields of an attribute table relative to the row ID field. As a very simple example of what I need; if the row IDs of 1, 2, and 3 corresponded to identical values in another field, I would need these values to be reordered so that the row IDs of 1, 2, and 3, corresponded to the inverted values of 3, 2, and 1, respectively. This is part of a model I am making so I cannot do this manually, I need some way of doing this using ArcGIS tools in modelbuilder.

I am using ArcGIS 10.2


texts: Nominatim search result difference between IE and Chrome

texts:
When searching for specific addresses in Nominatim I see different results in IE compared to Chrome. For example:

Im Singes 1 63225 Langen Deutschland

http://nominatim.openstreetmap.org/search.php?q=Im+Singes+1+63225+Langen+Deutschland&polygon_geojson=1&format=xml

the Url in both browsers is identical. In IE I get the following place as result

<place
    type="residential"
    class="highway"
    importance="0.63"
    display_name="Im Singes, Linden, Langen, Landkreis Offenbach, Regierungsbezirk Darmstadt, Hessen, 63225, Deutschland"
    lon="8.68072"
    lat="49.9842571"
    geojson="{"type":"LineString","coordinates":[[8.6814915,49.9844593],[8.6812535,49.9844064],[8.6810239,49.9843473],[8.68072,49.9842571],[8.6800859,49.9840417]]}"
    boundingbox="49.9840417,49.9844593,8.6800859,8.6814915"
    place_rank="26"
    osm_id="141557253"
    osm_type="way"
    place_id="97786309"
/>

and in Chrome it is

<place
    place_id="96505241"
    osm_type="way"
    osm_id="141557254"
    place_rank="26"
    boundingbox="49.9845947,49.9849051,8.6816261,8.6848355"
    geojson="{"type":"LineString","coordinates":[[8.6848355,49.984848],[8.684728,49.9848833],[8.6846498,49.9848986],[8.6845552,49.9849051],[8.6844579,49.9848991],[8.6831001,49.9847536],[8.6816261,49.9845947]]}"
    lat="49.9847536"
    lon="8.6831001"
    display_name="Im Singes, Linden, Langen, Landkreis Offenbach, Regierungsbezirk Darmstadt, Hesse, 63225, Germany"
    class="highway"
    type="residential"
    importance="0.52"
/>

Two different records with slightly different coordinates. Any idea?


texts: Calculating average stream slope at each location along stream using ArcGIS Desktop?

texts:
I'm trying to calculate the time of concentration for every point along a stream network. To calculate the time of concentration, I need the upstream channel length and channel slope.

I've been able to calculate the upstream channel length at each pixel using the Hydrology/Flow Length (upstream) tool (See image - red corresponds to longer cumulative channel length):

However, I have not been able to estimate the overall channel slope at each pixel (i.e. the slope between the pixel in question to the furthest-upstream headwater pixel). I know you can do this for individual basins (i.e. for one pixel), but my study area is very large (longest channel is 500+ km), so I dont want to generate a "watershed" for every single stream pixel, as this will be too time consuming.

I have also tried calculating the flow length with a slope raster as the "input weight raster," (and then dividing the result by the flow length raster (with no weight raster) to get the average slope along the stream, but this is far overestimating the overall stream slope.

Essentially, I need to find a way to make a map of the headwater elevation corresponding to any pixel in the stream network. Then I can subtract the elevation at each pixel and divide by the upstream channel length to get the overall stream slope.

(I'm using ArcGIS 10.1, Basic License)

Here is the error message I'm getting when I use the python script presented below (after running for 1 hour and getting through half of the points):


texts: Changing colors of polygon based on numbers in separate Excel sheet using ArcGIS Desktop?

texts:
I have a shapefile and a separate excel Sheet.

I want this shapefile to interact with the excel sheet. I have assigned numbers associated with colors on this sheet, what I would like to do is have the shape color change with the numbers in the excel sheet.

I believe this is something I can accomplish in ArcGIS 10 ModelBuilder but I don't know how to go about it.


texts: How to sample elevation data at 30m resolution at every 3m interval?

texts: I want to sample DEM data, 30m resolution, at every 3m under a drawn line and wondering if there is any built-in tool within ArcGIS 10.1. My take on this is to convert line (feature) to point and somehow define point interval to be 3m and then retrieve elevation values for each of those points. Please suggest on how to control the point interval or any other way of handling this task.


texts: Calculating the area of oyster drone imagery

texts:
I'm working on analyzing the area of an oyster reef and its development over the past years. The shown ortophoto is already loaded in QGIS. The unit of this file is degrees. That's where my main problem lies. For my work, I need to have the area of just the oysters in square meters. When I used the SCP classification tool, it showed me an area in the unit degrees^2. Is there a way or a factor to transform degrees^2 to m^2? Or do I have to use another tool to get my area in m^2? What's the best way to find out the area?
The targeted area can be just the pure oyster area or the whole oyster area not regarding the empty spaces in the inside of the reef. Both ways work!


texts: Extract Elevation of Points from DEM using Python QGIS or GRASS

texts:
I am utterly new to GIS, so please bear with me if I am providing incomplete information here.

I have a DEM file (.dt2). What I am looking to do is as follows:

take a series of points specified as lat, long
calculate their elevations using the DEM data
Output these points in as a list of ECEF coords

I need to do this in either GRASS or QGIS (using Python scripting, not the GUI).

I don't need a complete working code. I am only looking to be pointed in the right direction.


texts: Extract by mask, default naming and arcpy

texts:
The extract by mask tool in ArcGIS has only two required arguments. The first is inRaster and the second is inMaskData. The filename and format of the output raster are determined automatically. The following command will achieve this:

extracted_raster = arcpy.gp.ExtractByMask_sa(input_raster, extraction_buffer)

The output raster automatically gets a name like so: Extract_{in_raster_name}. It gets dumped into arcpy.env.workspace as expected.

In my case, I want my output raster to get dumped into in_memory, but I also want arcpy to automatically name it and pick the format for me. I tried doing this using two methods. The first consists of changing the workspace to in_memory before I run the above command:

   arcpy.env.workspace = 'in_memory'
   extracted_raster = arcpy.gp.ExtractByMask_sa(input_raster, extraction_buffer)

This throws an error, since it cannot find my input_raster, which is not located in_memory:

  ERROR 000865: Input raster: .\input\elsus.tif does not exist. Failed
  to execute (ExtractByMask).

The second method is to hardcode the output_raster location to in_memory. But using this method does not allow me to force arcpy to name the file:

arcpy.gp.ExtractByMask_sa(input_raster, extraction_buffer, r'in_memory\test')

I really don't care about the file name for this intermediate data. I prefer to store the features/rasters as objects and use the variable names instead.How can I force arcpy to handle the naming of the file, such that I only have to worry about the output object variable name?


texts: Moving data frame to particular X.Y location using real world coordinates?

texts: My mxd file has multiple dataframes. How do I move a given dataframe to a particular location within the main dataframe in the layout view to a particular x,y location given by user. The location x,y would dicatate the center of data frame to be moved.


texts: SQL WKT draws well with geometry, but not with geography data type

texts:
I have a following piece of code

DECLARE @g geometry;
DECLARE @borders geography;
SET @g = geometry::STGeomFromText('SOME WKT', 0);
SET @borders = GEOGRAPHY::STGeomFromText(@g.MakeValid().STAsText(),4326)
SELECT @g;
SELECT @borders;

Since it's too long to paste it here, WKT Can be fond at this link: https://justpaste.it/6qw0a

Can someone please explain to me, why it displays well when I draw it as geometry, but when I try to draw it as a geography, it displays entire world instead of a small group of islands.

Here's the screenshot:


texts: Publish imagepyramid on GeoServer SRS "unknown"

texts:
I want to publish an ImagePyramid created by GDAL.
gdal_retile.bat -v -r bilinear -levels 4 -ps 2048 2048 -s_srs "+proj=utm +zone=32 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs" -co "TILED=YES" -co "COMPRESS=JPEG" -targetDir "C:/Temp/pyramid/result" "C:/Temp/pyramid/image.ecw"

After adding the imagePyrimid source to GeoServer a shapefile is generated. A grid with image directory. If I add it to QGIS it is placed correctly but in GeoServer the native SRS in GeoServer is "unknown" so I fail to properly publish it.
I have followed the tutorial on geoserver.org as closely as possible.
https://docs.geoserver.org/stable/en/user/tutorials/imagepyramid/imagepyramid.html
I tried to explicitly specify the source srs_def without luck. I also tried without the -s_srs tag.
Is there something obvious I am missing here? I wonder why QGIS handles the coordinate system okay, when GeoServer struggles a bit.


texts:
I'm trying to print raster batch count number in PyScripter and ArcGIS 10.3
I'm using the next code:

import arcpy

carpeta = arcpy.env.workspace = r"C:\Base_de_datos\Cuencas_raster"

for raster in carpeta:
    with arcpy.da.SearchCursor(raster, ("COUNT")) as cursor:
        for row in cursor:
            count = row[0]
            print raster + "," + count

PyScripter give the next error:

Traceback (most recent call last):
  File "", line 16, in
RuntimeError: cannot open 'C'

Any idea to how resolve this?


texts: "Select by location" instead of SQL query for labeling?

texts:
Within ArcMap, is there a way to indicate what features to label via a "select by location" method instead of a SQL query?

I would normally approach the issue by adding an additional field in the feature class' attribute table, then populate it with something like a "Y" or "N" based on the features' location, indicating whether or not the feature should be labeled. However, in this instance, the feature class' attribute table is not editable and a new feature class cannot be exported.


texts: Find the intersecting area of more than 2 polygons in QGIS

texts: I have a map of 7 species distributions. I want to create polygons of the area where all 7 overlap, 6 overlap, 5 overlap etc. I created a polygon of where any 2 species overlap using the intersect tool under geoprocessing tools. However it would be incredibly tedious to do this for the area where more than 2 polygons intersect because I would have to merge all possible combination of shp. files. Is there a way to do this besides that?


texts: Tile server with dynamic mapnik rendering, caching, and pre-emptive rendering

texts:
I've been making some pretty (but computationally expensive) tiles with TileMill and I'd like to serve them up.

First, I served them straight from TileMill. Every tile gets rendered dynamically, which is simple, but slow.

Next, I exported a region as MBTiles then served it with MBTiles-server. Every tile is served statically, which is fast, but exporting all the tiles is too slow (multiple days).

I'm looking at TileStache, which combines dynamic rendering with caching. I imagine that this means the first time a region is looked at, it will be extremely slow, then quick thereafter.

I have plenty of spare CPU time (8 cores, not being used for anything else) and plenty of disk space (terabytes). I'd really like a server that can:

a) serve cached tiles if available
b) dynamically render tiles if necessary
c) continuously render tiles in the background to maximise the cache hit rate

Is mod_tile the only option? I'd prefer to avoid Apache, and mod_tile's documentation looks messy.

For context, I want to serve tiles for the whole of Australia, down to zoom 15.


texts: Deriving length values from GPX files in QGIS

texts:
I’m missing something real obvious and I apologize for revisiting a topic addressed in several Q&A’s here. I’m an experienced GIS user (not programmer) who’s new to QGIS (can’t afford Arc in retirement). I was given GPX files from a race organizer who wants to know whether they truly have a 50K course laid out. These were generated by Motion-X GPS on an iPhone and that program can also produce KMZ files.

The files had a number of small dropouts and had to be merged into a total race route, so I performed quite a few edits before moving on to determining the total distance of the current race route.

I have followed the manual, tutorials, and a few answers here several times now. In essence, the correct steps appear to be:

import the GPX tracks or KMZ file in CRS 4326
save to a shapefile in an equidistant CRS
add a new field to the attribute table for length, populated by
$length

All values in that field remain NULL. I tried bringing in the trackpoints and converting those to a line using point to one then repeated the steps above without success. I also tried dumping the coordinates from the GPX to a text file then reading that in as a CSV file before attempting the steps above. I tried importing a KMZ instead, without success. Every alternate attempt or workaround hasn’t achieved what I need.

I’m feeling pretty dim right now. The imported map layers (GPX or KMZ) display properly and register with other base geodata, of course, so it’s extra frustrating.

First few lines of a GPX:

<gpx xmlns="http://www.topografix.com/GPX/1/1" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" version="1.1" xsi:schemaLocation="http://www.topografix.com/GPX/1/1 http://www.topografix.com/GPX/1/1/gpx.xsd" creator="MotionXGPSFull 22.2 Build 4844R">
<trk>
<name>50K - 0 to 1</name>
<desc>Mar 28, 2015  7:38 am</desc>
<trkseg>
<trkpt lat="39.1847345" lon="-78.5107859">
<ele>308.437</ele>
<time>2015-03-28T11:38:28.034Z</time>
</trkpt>
<trkpt lat="39.1847345" lon="-78.5107859">
<ele>308.358</ele>
<time>2015-03-28T11:38:29.413Z</time>
</trkpt>


texts: Leaflet popup not pointed to marker when collapse

texts:
In my application I am using leaflet popup to display more details regarding the marker (Vehicle). User should be able to collapse and expand the popup.

Fist user click on the marker it will display the pop as

When click on collapse button, popup content will be collapse using jquary. But it will not pointed to marker correctly.

To do this I only change the css on the popup content using jqury. My collapse function be like.

function collapseDetail() {
    $("#serviceTable").css("display", "none");
    $("#serviceDataDownIcon").show();
    $("#serviceDataUpIcon").hide();
}

How can I avoid this and pointed the popup to the marker correctly ?


texts: one label for two-lane roads (osm, qgis, postgis)

texts:
I labeled my osm-lines-layer (PostGIS) with symbols for primary roads and motorways. While the symbols for the primary roads are well placed - with enough distance between each - the two-lane motorway is labeled with one label per lane (as you can see in the image).

Is it possible, to set something like a radius from each label, within no other label of the same type can be placed? Or can I just say: label one lane of a two-lane road?

I'm using qgis. The OSM-Data is imported to PostGIS with osm2pgsql.

EDIT:This new screenshot shows more detail:


texts: How to request features in a QGIS plugin?

texts: How does one go about requesting a feature to a certain plug-in?  The plug-in author is no longer active but he said I could request/implement the feature to someone who maintains plug-ins for QGIS.  The plug-in is 'Rectangles Ovals Digitising'.


texts: Split GeoTIFF into multiple cells with Rasterio

texts:
I'd like to create training data for a machine learning pipeline. To do this I want to take a large raster and split it into multiple equal sized cells.

I have a very basic loop for generating the pixel-wise coordinates, how do I generate a cell from these when splitting up a GeoTIFF?

def splitImageIntoCells(img, filename, squareDim):
 numberOfCellsWide = img.shape[1] // squareDim
 numberOfCellsHigh = img.shape[0] // squareDim
 x, y = 0, 0
 for hc in range(numberOfCellsHigh):
     y = hc * squareDim
     for wc in range(numberOfCellsWide):
         x = wc * squareDim
         # Need some method from Rasterio here
         # to crop at the given x and y...


texts:
Wondering if there is a code or tool that will allow me to format the ranges and labels for classes for Natural Breaks (Jenks) classification in ArcGIS?

I understand that if I alter the ranges it defaults to "manual" classification, however I'm just interested in getting the class ranges rounded to the nearest whole numbers so that its easier to read on the map.
In the images below, I like to round-off the ranges and labels 1,000-3,500, 3,500-7,000, 7,000-12,000, 12,000-25,000, etc. So this will be a "pseudo-Jenks" classification, but the goal is to make the numbers easy to read but still get the trends of a natural breaks classification.
Is this possible?


texts: Shifting overlapping lines with the same attributes in QGIS?

texts:
I have a layer of cables which overlap each other. I want to divide it so I can see them conserving their attributes.

I've just tried this method (Shifting display of overlapping lines in QGIS?), and it didn't help much because I can't select the lines, so I think maybe adding each line (with the new position) in a new shapefile, conserving the same attribute or by adding some instructions in the program code.

The only way I see for now is to really create different lines like in the picture above, with the same attribute.

it looks just like this but i can't select the lines or view its attributes


texts: Managed database won't validate in ArcGIS Server Manager

texts: I am running ArcGIS Server Manager and Portal for ArcGIS.  Unfortunately, some of my map services aren't working correctly, and I think it has to do with the Managed Database called ArcGIS_Data_Store.  It was just working fine the other day.  What should I do to fix this?


texts: Reading Sentinel-2 L2A product footprint (valid data) from metadata as shapely polygon?

texts:
Is the valid data footprint of Sentinel-2 products stored somewhere in the SAFE products metadata and if yes, how do I read it and convert it to a shapely.geometry.Polygon object in python?

EDIT: by valid data I mean the portion of the frame that is not 0


texts: Can't open specific layers with Qgis nor OGR

texts:
I'm working with mutliple sets of layers in order to create a new one.
As I need to load them directly from the internet, I use both /vsizip/ and /vsicurl/

But I'm stuck at loading one of them.
Here is what I do:

from osgeo import ogr
address_shp = ogr.Open("/vsizip//vsicurl/http://bano.openstreetmap.fr/BAN_odbl/BAN_odbl_01-shp.zip")
if address_shp is None:
    print ('Address : download successfull')
else:
    print ('Address : download failed')

parcels_shp = ogr.Open("/vsizip//vsicurl/https://cadastre.data.gouv.fr/data/etalab-cadastre/latest/shp/departements/01/cadastre-01-batiments-shp.zip")
if parcels_shp is None:
    print ('Parcels : download successfull')
else:
    print ('Parcels : download failed')

Output:

>Address : download successfull
>Parcels : download failed

I got the same result directly in QGis.
So, does /vsicurl/ handle https protocol? Or is the problem elsewhere?

OS : Windows 7 / Qgis 2.8.26 & 3.4.2


texts:
I am drawing a polygon with the drawing function within ArcGIS 10.1.  Is there a way to snap orthographically at 90 degrees while using the polygon draw command.

I added a picture of what i would like to do but as you can see the callout is not drawn at 90 degrees.


texts:
Could you please answer to these questions:
How to add TMS layer to Geoext?
Is it possible?

I couldn't find any information about it....


texts: Selecting features from a .shp file based on rasters pixel values

texts:
I am on windows 10, using QGIS 3.18.
I have a shapefile with multiple entries (polygons) and a raster with a few pixel values. I would like to use a "select by ... " to select the entries within the '.shp' that intersect with specific pixel values.
For example, try answering the question: Which entries (lines in my attribute table) intersect with "1"s in raster 'raster_name' (raster with 0 and 1s). In example 2, how to do the same but with rasters with a few different values (let's say 1 to 10) when I would like to 'select by ...' intersect with pixels with values equals '9'.
I got as far as:
'select by expression'
intersects($geometry, layer_property('raster_name', 'extent'))

With that, I could select every entry that intersects with a raster extent.
I also tried:
intersects($geometry, raster_name@band = value)

And it did not cam as a valid expression.


texts: Select by Attributes - If string has a letter

texts:
I'm running ArcMap 10.1.

I'm using a File Geodatabase.

I have a string coordinate field that I want to clean and some records have letters in it (mostly 'N' or 'E' or 'n' or 'e' e.g. '6213210N' or 'N6584312' but also some refer to map a grid reference e.g 'm29').

I want to select the records that have a letter regardless of it's position or what letter.

Why wouldn't this work?

UPPER("CoordinateField") LIKE '%[A-Z]%'

As far as I can see it works in for similar solutions but not in ArcGIS. As seen here.

Am I using incorrect notation for something stored in a File Geodatabase?

Hoping someone can shed some light on this. Much appreciated.


texts:
I would like to add the following package to my Python Enviroment: https://github.com/terraformer-js/terraformer/blob/master/packages/arcgis/README.md .
I have cloned the default eviroment. Via "Add Packages" I can only add packages from the list. Pip or conda do not work in the python console either. Does anyone have any ideas or experience with this?


texts: How to maintain the data at XML document after changing schema with Diagramer?

texts:
After exporting XML with schema and data I did some changes to the schema with ArcGIS Diagramer. now i need to import the new schema and all the data to a new empty file GDB. what are the best practices for this ?
 using "Publish" on diagramer save only the schema. How can i load the data most efficiently (i have more the 40 layers) from XML to geodatabase ?


texts:
The topojson file with my code shows only one propriety with one path,my topojson file looks like the sample . what I am doing wrong ?
I dont know if the problem with the topojson or the code as this my first project

the output show only one polygon where it must shows five

my code:

createMap("QasabatAmmanTempWithoutMadina.json");

function createMap(topojsonFile) {
    d3.json(topojsonFile, function(error, tj) {
    svg.selectAll("path").remove();

        //console.log(key in tj.objects);

        for(key in tj.objects) { features = tj.objects[key];console.log(tj.objects[key]); }
        console.log(features);
        var estados = topojson.feature(tj, features);
            console.log(estados)
        var b = path.bounds(estados);
        var s = 0.95/ Math.max((b[1][0] - b[0][0]) / width, (b[1][1] - b[0][1]) / height);
        var t = [(width - s * (b[1][0] + b[0][0])) / 2, (height - s * (b[1][1] + b[0][1])) / 2];

        projection
            .scale(s)
            .translate(t);

        g.selectAll("path")
            .data(topojson.feature(tj, features).features)
            .enter().append("path")
            .attr("d", path)
            .attr("class", "feature");
            //.on("click", click);

        g.append("path")
            .datum(topojson.mesh(tj, features, function(a,b) { return a !== b; }))
            .attr("class", "mesh")
            .attr("d", path);
    });
}

sample of topjson file:

  {"type":"Topology","arcs":[[[5629,4900],[26,-20],[21,-38],[-19,-32],[-11,-49],[-40,-82],[-44,-37],[-69,-61],[-53,-23],[-57,-10],[-74,3],[78,32],[78,30],[66,20],[84,6],[39,-14],[31,-58],[43,-68],[44,-36],[33,-30],[43,-37],[52,-47],[57,-40],[43,-36],[47,-16],[47,-18],[68,-6],[72,-6],[65,-4],[74,0]],[[4419,5688],[-31,4],[-113,106],[-12,29],[-20,25],[-27,16],[-38,5],[-171,-8],[-232,76],[-31,-61],[80,-4],[164,-72],[17,1],[180,-67],[27,-23],[34,-17],[36,-32],[24,-34],[9,-38],[4,-29],[6,-29],[6,-31],[-8,-18],[15,-44],[-8,5],[-32,16],[-14,7],[-41,-3],[-36,-13],[-24,-10],[-19,-11],[-24,-3],[-42,14],[-34,10],[-39,10],[-40,20],[-21,18],[-15,19],[-12,27],[-10,17],[-16,18],[-24,5],[-25,-1],[-22,-2],[-24,-18],[20,5],[28,-10],[19,-11],[15,-14],[11,-14],[-5,-20],[-53,32],[-98,-36],[-245,72],[-43,36],[-94,62],[-157,28],[-153,54],[-100,83],[-35,4],[-86,-21],[-63,7],[-55,1],[-98,4],[-101,24],[-72,28],[-36,11],[-28,11],[-15,3],[-98,-6],[-66,1],[-29,9],[-58,26],[-58,37],[-74,65],[-46,62],[-45,84],[-47,136],[-30,135],[-6,69],[6,61],[24,138],[64,122],[77,111],[159,158],[258,228],[68,59],[98,62],[61,25],[80,11],[129,-8],[47,10],[58,23],[99,47],[51,29],[85,58],[95,62],[42,21],[103,39],[193,62],[-26,-196],[3,-134],[0,-57],[5,-114],[-4,-137],[-12,-58],[0,-43],[34,-103],[53,-144],[57,-74],[59,-52],[163,-104],[39,-47],[69,-67],[35,-54],[7,-41],[-32,-103],[1,-20],[35,-121],[16,-29],[87,-83],[24,-15],[26,-2],[30,15],[13,-37],[39,-38],[4,-23],[49,-128],[-3,-30],[19,-46],[17,-70]]],"transform":{"scale":[0.000012395198638604345,0.000011753131339408293],"translate":[35.86398878517448,31.89607962076447]},"objects":{"Zahran_Trusted":{"type":"GeometryCollection","geometries":[{"arcs":[[0]],"type":"Polygon","properties":{"Level":3,"Name_1":"Amman","Name_2":"Qasabat

  Amman","Name_3":"Zahran"}}]},"YarmukAmman_Trusted":{"type":"GeometryCollection","geometries":[{"arcs":[[1]],"type":"Polygon","properties":{"Level":3,"Name_1":"Amman","Name_2":"Qasabat
  Amman","Name_3":"Yarmuk
  Amman"}}]},"RasAlEin_Trusted":{"type":"GeometryCollection","geometries":[{"arcs":[[2]],"type":"Polygon","properties":{"Name_1":"Amman","Name_2":"Qasabat
  Amman","Name_3":"Ras Al
  Ein"}}]},"BadrAmman_Trusted":{"type":"GeometryCollection","geometries":[{"arcs":[[3]],"type":"Polygon","properties":{"Level":3,"Name_1":"Amman","Name_2":"Qasabat
  Amman","Name_3":"Badr
  Amman"}}]},"Abdali_Trusted":{"type":"GeometryCollection","geometries":[{"arcs":[[4]],"type":"Polygon","properties":{"Level":3,"Name_1":"Amman","Name_2":"Qasabat
  Amman","Name_3":"Abdali"}}]}}}]

]


texts:
I'm trying to delineate floodplains in a catchment using "Stream WSE from point WSE Measurements" as in (https://blogs.esri.com/esri/arcgis/2014/03/17/floodplain-delineation-using-only-arc-hydro-models) but I got strange artifacts in the resulting flood depth raster (in the areas outlined with the rectangles, the orange line is the river center line):

In the biggest rectangle at the bottom of the image, the depth between each side of the "breaking line" rises sharply of 8 meters. I don't really understand because the underling DEM is smooth at this location with only a difference of 1 meter.

Any idea from where this kind of artifacts could come from?


texts: Is there still CoverageStack object in GeoTools?

texts:
This is the first time I'm dealing with Geotools 18. I need to work with a GridCoverage Stack object. I saw that there was one in a very old release (until v2.7.4 CoverageStack).

Is there a similar object in the nowadays release?


texts: Determining Austrian roads which are accessible for pedestrians from OpenStreetMap?

texts:
I've looked at the attributes/tags of the OSM highways but couldn't find any information about my issue.

Is there a way to receive the information which roads (of the OSM highway key) are accessible by pedestrians? I know there is the "pedestrian" value, but I would like to find out ALL roads where pedestrians are allowed to walk (they can not only walk on pedestrian roads, sometimes also on secondardy roads, etc, wich are accessible for cars and for pedestrians both).

Is it possible to get this information using OSM data? Or is there a good workaround?

My area of interest is the country of Austria.
My target is to receive all roads where pedestrians can walk...so not only footways...they can also walk on roads which are usually for cars (and pedestrians). Excluded for example should be all motorways (and other roads).


texts: Conflict standart leaflet polygon and polilyne vs leaflet-dvf plugin

texts:
I found some bug: if you are using leaflet-dvf plugin and want to draw standart leaflet polygon, you will get error in your browser's console like

Uncought RangeError: Maximum call stack size exceeded

it will reference to leaflet-dvf.markers.min.js:6
TextFunctions._updatePath

it is 1838th string in here
https://cdnjs.cloudflare.com/ajax/libs/leaflet-dvf/0.2/leaflet-dvf.markers.js

How i can handle that error?


texts: How do I calculate a xyz-position of one GPS position relative to another?

texts:
I have two GPS positions (lon, lat, alt) and want to get the position of one relative to the other one. (in meters, as x, y, z)

How do I do that?

Edit: The locations have a distance of up to three kilometers
And: I do NOT just want the distance.


texts: Adding query layer to MapView in Jupyter Notebook returns "RuntimeError: Cannot infer layer: will not be added to map"

texts:
I am working within a Jupyter notebook, trying to display a query layer (querying a SQL database) into a map view within the notebook. I receive a Runtime Error (described below) when I try to add this layer to the map view.
Using the following code (with specific names changed for privacy reasons) I can successfully create a map and add a layer from my organization's enterprise portal.
# Import Statements
from arcgis.gis import GIS
import arcpy

# Instantiate Connection to Portal
gis = GIS("portal_url", username="my_username")

# Get item from portal
portalItem = gis.content.get('itemID')

# Create map and add portal item to it
m = gis.map()
m.add_layer(portalItem.layers[0])

I can also successfully create my query layer using the following code, including testing that makeQueryLayer returned a feature class with geometry as expected. This same query was copied directly from a successful query layer being created in ArcGIS Desktop 10.7.1.
# Create query layer input variables
out_layer_name = "queryLayerOutput"
sql_Query="""SELECT OBJECTID
                    ,Field1
                    ,Field2
                    ,Shape.STEnvelope() as MBR
                FROM SQLDatabase.Table
                WHERE Shape.STIntersects(geometry::STGeomFromText('LINESTRING(""" + lat_long_string + """)', 4326).STEnvelope()) = 1""" # lat_long_string being the coordinates of my area of interest in proper 'E S, W N' format
oid_fields="OBJECTID"
shape_type="POLYGON"

# Creating the query layer
MBRQueryLayer = arcpy.MakeQueryLayer_management(geodatabase,
                                                out_Layername,
                                                sql_Query,
                                                "OBJECTID",
                                                "POLYGON",
                                                "4326",
                                                spatial_reference=spatialReference)

desc = arcpy.Describe("queryLayerOutput")
print("Feature class type: " + desc.featureClass.featureType) #returns 'Simple'
print("Shape field name: " + desc.featureClass.shapeFieldName) #returns 'MBR'
print("Shape type: " + desc.featureClass.shapeType) #returns 'Polygon'

Where I encounter problems is when I go to add the query layer to the map.
m.add_layer("queryLayerOutput")

Which returns the following error, the key line of which is "RuntimeError: Cannot infer layer: will not be added to map":

The documentation for MapView.add_layer() implies that I should be able to add a FeatureLayer, which to my understanding is what makeQueryLayer() produced. Does anyone understand why my layer is unable to be added? My best guess would be that somehow my data is not in an acceptable format for add_layer() to use. Somehow not the "right kind" of FeatureLayer???
I am working in the ArcGIS Jupyter Notebook environment (ArcGIS 10.7.1) on Windows 10.


texts:
I have created one script which set a defination query at certain location in script.
I use lyr.definitionquery = "some expression"

After running the script i want to remove the definition query . But not getting any hint to do it.


texts: Google Earth Engine - Sentinel NO2 mean, min, max chart

texts:
I'm using the code bottom to create NO2 Charts of a region.
Please, do you know if it possible to group the three graphs into one?

var collection = ee.ImageCollection('COPERNICUS/S5P/NRTI/L3_NO2')
.select('NO2_column_number_density')
.filterDate('2019-01-01', '2020-04-30');

var mean=ui.Chart.image.series(collection, geometry, ee.Reducer.mean(), 30).setOptions({title: 'NO2 mean'});
var max=ui.Chart.image.series(collection, geometry, ee.Reducer.max(), 30).setOptions({title: 'NO2 max'});;
var min=ui.Chart.image.series(collection, geometry, ee.Reducer.min(), 30).setOptions({title: 'NO2 min'});;

Link to code

Thanks in advance to those who can help me