texts:
I am very new to GIS.

I have a dataset in Excel containing XY coordinates of points.  These coordinates are said to be in WSG84, but have different decimal places lengths. Longest of them have 9 decimal places. They have 3 digits before decimal points. So a 9 digit coordinate is XXX.YYY.ZZZ or a 6 digit one is XXX.YYY. I want them all to be in decimal coordinate system. Is there way that I can do this automatically in Excel?

When I try to convert the data to XY points, ArcMap disregards the decimals even though both the shape layer and the points layer are in WSG84 format. Thus, drawing the points to a totally unrelated space.

Example Coordinates from the dataset;

362.935.938 339.309.883
363.010.152 3.396.492
404.244.154 264.695.215
40.63   26.83
41.128.471  29.145.429
41.305.968  286.615.491
411.788.295 291.018.963
4.102.785   2.906.803
4.112.674   2.915.039
41.019  28.571
413.621.671 285.677.379


texts:
Using the GDAL Georeferencer QGIS 3.10.10 I have georeferenced an old map, so now it is rotated a bit. The NoData fields for the "frame" are now black (RBG = 0,0,0). I would like to get rid of the frame, but if I turn the RGB-Value of 0, 0, 0 to 100% transparency, some areas of the actual map, which are also black will become totally transparent too.
I vaguely remember solving this problem some years back by loading the georerenced TIFF into GIMP and filling the "frame" with some other colour that is not in my map image and then setting this colour to transparent in QGIS. I tried this this time, but the georeferencing is lost.
Does anyone know, how to remove the black frame without loosing the black pixel of my map?


texts: Snapping point to point based on attribute?

texts:
I'm using ArcMap.

I have two identical point feature classes, one snapped to a line and the other not. Each point represents a telecommunication pole and the line represents fiber optic cable. I sequence the poles to the fiber optic cable line using the 'Route Events GeoProcessing Wizard' which by design, snaps all points to the fiber cable. I then have to manually snap the points that don't belong on the fiber line manually, which is very time consuming. Is there a way to snap two identical point feature classes with unique ID's together, instead of snapping based off proximity?


texts:
I need to somehow combine my Spatialpolygons with my depth contour, but not sure how do to so. Maybe create a SpatialPolygonsDataframe?

This is what I've done so far:

> class(subarea0) #This is my area of interest (Eastern Canadian Arctic Sea)
[1] "SpatialPolygons"
attr(,"package")
[1] "sp"

> extent(subarea0)
class       : Extent
xmin        : -82.21997
xmax        : -57.21667
ymin        : 60.2
ymax        : 78.16666

library(marelac)
data("Bathymetry")#World bathymetric data in library (marelac)
names(Bathymetry);class(Bathymetry);str(Bathymetry)
[1] "x" "y" "z"
[1] "list"
List of 3
 $ x: num [1:359] -180 -179 -178 -177 -176 ...
 $ y: num [1:180] -89.5 -88.5 -87.5 -86.5 -85.5 ...
 $ z: num [1:359, 1:180] 2853 2873 2873 2873 2873 ...

  raster_bath<-raster(Bathymetry)#Transformed into a raster layer (not sure if this is the best way to achieve what I want)
    extent(raster_bath) <- extent(subarea0)#Transform the extend of my raster to the extend of my SpatialPolygons
v <- extract(raster_bath, subarea0)#Extract data from my Raster Layer for the locations of my SpatialPolygons

v is a list with depth and I am not quiet sure how/under what form to rebind this info with my spatial polygon... Ideas? Other suggestions?

EDITS

I apologize in advance, looks like I've been having a hard time to explain clearly what I was hoping to get. I don't want to summarize my depth info for each of my polygon but rather provide a detailed bathymetry within my area of interest.
I am definitely not sure if what I did up to now is the way to go, but ultimately what I would like to achieve is to end-up with a plot of my polygon with my depth information. Something like this  (thanks to Simbamangu), but without the depth info outside of my polygon. Basically, I just want the bathymetry info comprise in my overall area at the finest scale possible.
Should I first transform my SpatialPolygons into one and then into some grid format to then be able to clip with my raster layer?

Thank you!


texts:
I have a postgis table containing polygons. I want to transform my polygons  to 4 points polygons : retrieve the smallest MBR.

Here comes an illustration of what I want :

My original polygon is in pink color.
I guess, it will be quite easy to get the MBR (with ST_Envelope and BBox function)
But is there an easy way to get what I want?

Any ideas? Many thanks


texts: Running QGIS3 on Mac, multiple plugin problems, ModuleNotFoundError: No module named 'osgeo'

texts:
Couldn't load plugin 'Qgis2threejs' due to an error when calling its classFactory() method.

ModuleNotFoundError: No module named 'osgeo'
Traceback (most recent call last):
  File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 337, in startPlugin
    plugins[packageName] = package.classFactory(iface)
  File "/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/Qgis2threejs/__init__.py", line 27, in classFactory
    from .qgis2threejs import Qgis2threejs
  File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 674, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/Qgis2threejs/qgis2threejs.py", line 29, in
    from .q3dviewercontroller import Q3DViewerController
  File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 674, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/Qgis2threejs/q3dviewercontroller.py", line 23, in
    from .export import ThreeJSExporter
  File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 674, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/Qgis2threejs/export.py", line 28, in
    from .exportdem import DEMLayerExporter
  File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 674, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/Qgis2threejs/exportdem.py", line 30, in
    from .propertyreader import DEMPropertyReader
  File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 674, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/Qgis2threejs/propertyreader.py", line 27, in
    from .qgis2threejscore import calculateDEMSize
  File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 674, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
  File "/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/Qgis2threejs/qgis2threejscore.py", line 25, in
    from osgeo import gdal
  File "/Applications/QGIS3.app/Contents/MacOS/../Resources/python/qgis/utils.py", line 674, in _import
    mod = _builtin_import(name, globals, locals, fromlist, level)
ModuleNotFoundError: No module named 'osgeo'

Python version: 3.6.5 (v3.6.5:f59c0932b4, Mar 28 2018, 03:03:55) [GCC 4.2.1 (Apple Inc. build 5666) (dot 3)]
QGIS version: 3.2.2-Bonn Bonn, exported

Python Path:
/Applications/QGIS3.app/Contents/MacOS/../Resources/python
/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python
/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins
/Applications/QGIS3.app/Contents/MacOS/../Resources/python/plugins
/Library/Frameworks/Python.framework/Versions/3.6/lib/python36.zip
/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6
/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/lib-dynload
/Library/Frameworks/Python.framework/Versions/3.6/lib/python3.6/site-packages
/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python
/Users/indyatblake/Desktop/QGIS/Workshop
/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/vrtbuilderplugin
/Users/indyatblake/Library/Application Support/QGIS/QGIS3/profiles/default/python/plugins/vrtbuilderplugin/vrtbuilder


texts: Deleting records from table where ST_DIstance from another record in another table is greater than a certain amount

texts:
I'm currently learning SQL on the run, and have come across a situation where I'm trying to delete records from a table depending on a function result. I can select the records of interest with the following query -

select ST_Distance(a.reporter_location, b.report_location) as dist, b.hex, b.flight, b.report_timestamp, b."isMLAT" from reporter a, planereports b where a.name like 'Home1%' and ST_Distance(a.reporter_location, b.report_location) > 400000.0 order by b.report_timestamp;

but I'm unsure how to place this into a conditional for a delete statement. Any ideas? Brain is turning to mush.


texts: Rendering a map in QgsComposition on Win 7 / 8 in QGIS 2.2.0 with QgsComposerMap gives blank map

texts:
I'm writing a plugin for QGIS 2.2.0 which should render a pdf file with a map and some text.

The code is working under Linux (current Arch) and Windows XP, but it's not working on Windows 7 or 8. Output PDF on Win 7/8 has all labels, but no map, only blank space.

I don't know if it's relevant, but I want to render only raster layers.

The plugin is to be runned from QGIS, so it has access to self.iface variable, but I can also use standalone QgsMapRenderer.

What am I missing?
Please help me! ;)

Here is the code I use:

mapRenderer = self.iface.mapCanvas().mapRenderer()
mapRenderer.setLayerSet(input_map_layers) #list of map layers id's
mapRenderer.setOutputUnits(QgsMapRenderer.Millimeters)
mapRenderer.setOutputSize(QSize(map_width_on_print_mm,
map_height_on_print_mm), 300)
mapRenderer.setExtent(rect)

c = QgsComposition(mapRenderer)
c.setPlotStyle(QgsComposition.Print)
c.setPaperSize(210,297)

composerMap = QgsComposerMap(c, 14.7026,113.85,map_width_on_print_mm,map_height_on_print_mm)
composerMap.setFrameEnabled(True)
composerMap.setFrameOutlineWidth(0.5)
c.addItem(composerMap)

printer = QPrinter()
printer.setOutputFormat(QPrinter.PdfFormat)
printer.setOutputFileName(unicode(input_data_fields['output_pdf_filename']))
printer.setPaperSize(QSizeF(c.paperWidth(), c.paperHeight()),
QPrinter.Millimeter)
printer.setFullPage(True)
printer.setColorMode(QPrinter.Color)
printer.setResolution(c.printResolution())

pdfPainter = QPainter(printer)
paperRectMM = printer.pageRect(QPrinter.Millimeter)
paperRectPixel = printer.pageRect(QPrinter.DevicePixel)
c.render(pdfPainter, paperRectPixel, paperRectMM)
pdfPainter.end()


texts:
I have a MSSQL (2008 R2) SQL view that is visible and displayable within QGIS (LTR ESSEN 2.14.15) itself. However if I attempt to style the results by classification the map display disappears. If I add any form of querying to the table it disappears.

Loading the Feature table for review
shows the inherent issue- all the fields are returning an error within QGIS.

The view behaves perfectly within SQL Management Studio and within QGIS itself
as you can see above, it is up to a point. So I'm discounting //most// connection issues.

But I can't apply any styling or additional queries to the table whilst it is returning a table full of Errors.

Sadly this doesn't appear to be specific to this one particular view, but any view I'm pushing through to QGIS.

Has anyone ever encountered this type of behaviour before? And any concept of how to fix?


texts: Route with distance marker for each Km. (ArcGIS)

texts:
In need to create a route from polyline where each Km in the route is identified.

This is what I did using the create routes tool:

This is the output from the tool:

I was expecting a listing with a unique ID for every 1Km segment.

Any help here.


texts: What kind of high level semantic free GIS layers out there?

texts:
I'm looking for free GIS data that contains some high level information such as type of vegetation or type of soil or some other visual characteristics of the neighborhoods. I need more information than Google/Bing/Yahoo Maps. Do you know any GIS server that has this kind of information in it?

Thanks,


texts: Rendering features based on relationships using ArcGIS for Server?

texts:
We are trying to determine an efficient way of rendering features depending on values in a different table. Background is this: We have a few thousand features in a MapService on an ArcGIS Server and access that via ArcGIS Runtime apps. Let's say those features represent buildings. From time to time, there is work done in those buildings (repairs, cleaning, ...). This work is being logged in a second table with a reference to the building (via the objectid) and also structured into work packages. These work packages are either recurring or can be ad hoc (e. g. if there's a sudden equipment failure that needs to be dealt with immediately).

What we want do now is render the buildings differently based on their entries in this second table (e. g. render them red, when the last repair failed; render them green when they were cleaned this week, render only those in a specific work package, ...).

Obviously, a simple renderer is unable to do this, as we are only able to define it on the attributes of the feature (the buildings in this case).

Right now we're wondering how to effectively render these features. We are unable to change the data structure, so we cannot just add some value to the feature itself. The options we see are the following:

Prerender the layers. This would be done via a batch program or being intiated by the user. This would mean that we'd have to publish the MapService with those layers inside of it, which is a problem because we'd be unable to quickly add new work packages as they'd need their own layers.
Render them on the device. This might have an impact on performance, though I don't think that's going to be a real problem. One problem though, is that by doing this, we'd have differing presentation depending on which device you use (e. g. the people in the back office use a desktop app, while the teams in the field are mostly equipped with Android devices). Also, there would be the problem of how to efficiently render them on the device as one would have to query the database on how to render which feature.
Use a Server Object Interceptor. This way, we could inject a virtual field into the features and use that for rendering on the device. We don't have any experience with SOIs, though and are not totally convinced in regards to performance as we'd have to do additional database queries on each request.

Is there any obvious option that we're missing? If not, does anyone have experience with our choices so far or is otherwise able to shed some light on which one's the way to go?


texts: Clipping where entire polygon is outside of bounds using QGIS?

texts:
I'm trying to use QGIS atlas to generate a whole bunch of maps, based on a grid of polygons. What I planned to do was create a new grid square polygons within QGIS, clip the polygons that don't fall within the base layer, and then generate the atlas files from there.

However, when I try to Clip, where a polygon that I want to use as a reference for the atlas settings partially falls outside the base area, that polygon is clipped to the exact bounds of the base layer.

Is there a way that I can clip where any polygon that entirely or partially is covered by the base layer, the entire square polygon is retained, and only those that fall completely outside the base layer are clipped?


texts: Terrain correction (flattening) of Sentinel 1 images

texts:
Using the script available in this forum (belowlink), I tried for terrain correction (flattening) of Sentinel 1 images

Getting local incidence angle from Sentinel-1 GRD image collection in Google Earth Engine

Here is my working script.

https://code.earthengine.google.com/0a0ae5e730afd766d0c2eb4c4a07c937

Unfortunately, why the resulted gamma0dB is all white?

Could someone look at above link and solve it?

I really want to return terrain corrected (flattened) gamma0dB.


texts: GeoServer map layers/tiles Loaded in the screen or Not

texts:
I'm using OpenLayers - 2.11,GeoServer 2.13, ExtJS 3.4, and Java for developing GIS based web applications.
Please guide me to check if all my GeoServer map layers/tiles Loaded in the screen or not.
While dragging map is taking time to load next tiles. So that time i want to mask all the map body and unmask after loading all the map tiles.


texts: Creating tiles in the mapon from geojson

texts:
I generate a tile from postgis. I'm trying to superimpose data from geojson. But when data is superimposed, an error occurs:

Mapnik LOG> 2018-02-12 14:33:14: feature_style_processor: Style=way_style required for layer=way_layer does not exist.

Tile is created, but without data from geojson. Here is the code I get the data from the GeoJson-file:

def _get_graph_way(self):
    s = mapnik.Style()
    r = mapnik.Rule()
    way_line = mapnik.LineSymbolizer()
    way_line.stroke = mapnik.Color('#ffffff')
    way_line.stroke_width = 4
    r.symbols.append(way_line)
    s.rules.append(r)
    self.my_map.append_style('way', s)
    way_layer = mapnik.Layer('way_layer')
    way_layer.styles.append('way_style')
    ds = mapnik.GeoJSON(file="geo.json")
    way_layer.datasource = ds
    self.my_map.layers.append(way_layer)

This is an example of data in a GeoJson-file:

{
    "type": "FeatureCollection",
    "features": [
        {"id": 1,
            "type": "Feature",
            "geometry": {"type": "LineString", "coordinates": [[50.192642, 53.224538], [50.193331, 53.224360]]}, "properties": {"is_exit": true, "is_entry": true, "relations": [2, 3]}}
    ]
}


texts: How to use map service legend on ArcGIS Online?

texts:
How can I use the legend of a mapservice as part of a AGOL generated webmap?

This is the map I have created which does not have a legend

And this is the map service in use
[code]https://203.177.0.139/ArcGIS/rest/services/Leng/2013_Road_Data_elf/MapServer/[/code]

Pretty sure it has a legend because this has the correct content
[code]https://203.177.0.139/ArcGIS/rest/services/Leng/2013_Road_Data_elf/MapServer/legend[/code]

Could this still be a certificate issue? Our certificate is not signed or anything. So we are placing the map service address as an exception to an internet browser for it to allow access to "unsigned" content. By doing this we are able to use the AGOL webmaps with mapservice from an unsigned source.

What doesnt seem to make sense is why, after doing the exception list, would AGOL allow features to be drawn but not the legends (or the info configured for the pop-up dialog, as we also recently found out). What specific about these that makes them excluded?


texts: how to set the geotransform using GDAL in python

texts:
I am trying to create an image from a numpy array and georefenrece it so that I can import it into qgis. I have the latitude and longitude in degrees for the top left corner of the array as well as the spacing for the grids in metres.
My understanding  is that I can use the  SetGeoTransform function to set these values and then set the projection which, in my case is lcc. The problem is that the image lands on a completely different part of the canvas.
For the

  SetGeoTransform((originX, pixelWidth, 0, originY, 0, pixelHeight))

is it correct to set the origin in degrees and the width in metres?

Can you see what elese I am doing incorrectly?

THis is an example of the code:

import osr
import numpy
import gdal
format = "GTiff"
driver = gdal.GetDriverByName( format )
dst_ds = driver.Create( 'dst_filename.tif', 134, 164, 1, gdal.GDT_Float64)
dst_ds.SetGeoTransform( [ -10, 8000, 0, 45, 0, -8000 ] )
srs = osr.SpatialReference()
srs.ImportFromProj4('+proj=lcc +lat_1=45 +lat_2=45 +lat_0=45 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs')
dst_ds.SetProjection( srs.ExportToWkt() )
raster = numpy.zeros( (164, 134), dtype=numpy.uint64)
dst_ds.GetRasterBand(1).WriteArray(array)
# Once we're done, close properly the dataset
dst_ds = None


texts:
I have three 'layers' I would like to render on a map.

1) County-level estimates for various ACS 5-year variables (acquired through the tidycensus package) (in pseudocode below: data)

2) State boundaries, generated by aggregating 1) to the State-level (in pseudocode below: states)

3) The Great Lakes (acquired from here) (in pseudocode below: lakes)

When I map 2) and 3) using the following:

  tm_shape(states) +
    tm_borders(col = "black",
               alpha = 0.5) +
  tm_shape(lakes) +
    tm_fill("skyblue")

...it properly renders the state outlines and lakes. However, when I render 1), 2), and 3) using the following...

tm_shape(data) +
    tm_polygons("variable") +
    tm_legend(legend.position = c("right", "bottom")) +
    tm_layout(title = "var name",
              title.size = 1.1,
              title.position = c("center", "top"),
              inner.margins = c(0.06, 0.10, 0.10, 0.08)) +
 tm_shape(states) +
    tm_borders(col = "black",
               alpha = 0.5) +
 tm_shape(lakes) +
    tm_fill("skyblue")

...the lakes are always omitted, no matter the ordering. Has anyone encountered this before? Am I missing something about the ordering of 'layers' in tmap?


texts:
A friend and I are currently creating a pretty straight forward LINQ to ESRI library that sends HTTP queries to an ArcGIS MapServer and then parses the JSON response appropriately.

Everything is currently working beautifully except one thing; multiple Where statements. This seems like a simple syntax issue, but I can't seem to get it working, or find out if it's supported. I'm trying to do something like this.

Doesn't Work:

MapServer/22/query?where=MY_FIELD='Whatever'&where=MY_OTHER_FIELD='Whatever'&f=json

Does Work:

MapServer/22/query?where=MY_FIELD='Whatever'&f=json

Also Works:

MapServer/22/query?where=MY_OTHER_FIELD='Whatever'&f=json

I have tried several different syntax versions, but can't seem to get a multiple Where statement to work. This is my primary source for breaking down these queries. If anyone knows the syntax to do this, or knows the source that says it's not possible.


texts: How I can insert point to LINESTRING in PostGIS? For example, I have line like LINESTRING(0 0, 10 0) and point which lies on line segment: POINT(5, 0). I need to add this point to line. Expected result is LINESTRING(0 0, 5 0, 10 0). It looks like ST_AddPoint(geometry linestring, geometry point, integer position); can helps me but I don't know how I can determine position value in case if line will not so simple.


texts: What is difference between OGC and ISO standards - why do some open specs carry both standards?

texts: A good example is the Simple Feature Access - "OpenGIS® Simple Features Access (SFA), also called ISO 19125". Why use one over the other?  And why have both?


texts: Google Earth Engine: unexpected behavior of .mosaic()

texts:
I'm trying to get the Snow Cover extent from Sentinel-2 imagery using GEE.
Here is my script: https://code.earthengine.google.com/1fbe39044e73c262fe935db00dc2adb3

The script performs the following:

Create a Cloud Mask (expression based on S2 bands)
Calculate the NDSI (expression based on S2 bands)
Add Sentinel-2 imagery and sort them chronologically.
Use mosaic.() to create a mosaic with the most recent images
Create a Snow Cover Extent layer (expression based on the above
created bands)

The problem is when I visualize the results I get strange values for pixels in the image overlap areas. The visualized values are different from the band values (used Inspector to check).
For example, the white (cloud) strip at the right side of the image (see script) have a value of 1 (cloud) which is correct but the visualization is wrong (should be in red).

Even stranger, when I zoom in and out the visualization colors change!!

This makes me doubt what the mosaic () method really does. I thought only the most recent image is used.

can someone please look at the code and help me figure out what is going on ?

Thanks


texts: Breaking closed line when intersecting with polygon using SAGA GIS?

texts:
I have a number of polygons (in lines shape) which represent the new condition of an area and a numder of polygons (in polygon shape) which represent the old condition of the same area.

How can I break the new polygons (line shape) at the point that intersect with the old polygons (polygon shape) using SAGA GIS?


texts: Exporting a web map with Google Maps Satellite using Qgis2Web?

texts: I am really new to QGis, is it possible on QGis export a map using qgis2web with the Google Map Satellite layers? I tried several times but is not working.


texts: Coordinate system changes after convertion from feature to raster?

texts:
For an hydrologic model I'm supposed to add a DEM as a .tiff data and another raster as .tiff data for roughness. So I exported my DEM as a .tiff without any problems. I have the roughness Information as a .shp file. So I used the "Feature to Raster" tool, to convert my .shp into a raster and then I exported the raster as a .tiff. When I try to add my roughness .tiff I get the error message  "Coordinate System of Raster does not match DEM".

According to the sources they have the same coordinate system. Before I started converting and exporting my layers they also have the same coordinate system.
My theory is, that with the convertion process, the coordinate system has been changed without indicating it. Is that possible?


texts: Modifying GeoNode metadata - can't overwrite template

texts:
I am following this Developer-Workshop-Sheets to add an additional metadata field to my GeoNode.
https://geonode.org/dev-workshop/#/5 (this might be lost in the future, but isn't that important to the question, so I won't copy it).
I am stuck at step 5 where I have to copy a builtin geonode template, modify it and put it in my apps template directory. Usually this is how it works to modify the stock geonode. But this time I have no luck. Putting site-packages/geonode/layers/templates/layouts/panels.html into mygeonode/mygeonode/templates/layouts/panels.html won't be enough to replace it. My changes are just ignored.
I dig into the geonode code and learned that layers/templates/layouts/panels.html is called from layers/templates/layers/layer_metadata.html where it says:
  <div id="mdeditor_form" class="col-md-12 form-controls">
    {% form layer_form using "layouts/panels.html" %}
    {# layer_form|as_bootstrap #}
  </div>

My guess, as a total newby to django is, that this way of pulling one template into another isn't following the "regular" template replacement that django is usually following and is just using the path literally (but it can't be a relative path either so I am a little bit lost...). But I can't find anything about the term "using" in the documentation, mainly because that general term is used on every single page.
If my guess is somewhat right, could somebody point me into a direction on how to make shure I can use the usual template replacement logic here?
And if I am totally on the wrong track, and this might be an issue of first-found-first-served (and a wrong order of custom provided templates before stock templates), where in GeoNode would I find the definition of template serving order for namespaces?
I am trying to accomplish this with GeoNode 3.0 (based on Django 2.2).


texts:
I am hoping to open the following polygon layer in QGIS:
https://www.arcgis.com/home/item.html?id=764b72f3f6754e8fb1360c95bdf18918
The download is a 'pitem' file.
Are these files exclusively for ArcGIS, or is it possible to open in QGIS?


texts: Problem in ubuntu loading geoserver wms layer in local server with leaflet

texts:
I am just starting to work with geoserver. I am trying to do the same openlayers-geoserver example but using leaflet to view the map instead of openlayers:

Publishing a shapefile

This code is not working:

<!DOCTYPE html>

<link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.7.3/leaflet.css" />

<title>Leaflet Web Map</title>

I just get this clean-empty screen in linux ubuntu browser:

What could be the problem?


texts: Why is QGIS crashing when I export this image?

texts:
I'm trying to export a PNG of the OS 250k raster with OS Open Rivers over the top. Trying to export the entirety of the national grid caused it to crash, so I tried again with the bottom half of Great Britain and it's still crashing.
Here's the crash report:

Crash ID: 1d5ed672dc1b64fc24118b59e560969ad8cbd8bb
Stack Trace
QVectorPath::convertToPainterPath :
QVectorPath::convertToPainterPath :
QVectorPath::convertToPainterPath :
QVectorPath::convertToPainterPath :
QRasterPaintEngine::fillRect :
QRasterPaintEngine::drawImage :
QgsRasterDrawer::drawImage qgsrasterdrawer.cpp:174
QgsRasterDrawer::draw qgsrasterdrawer.cpp:105
QgsRasterLayerRenderer::render qgsrasterlayerrenderer.cpp:274
QgsMapRendererCustomPainterJob::doRender qgsmaprenderercustompainterjob.cpp:312
QgsMapRendererCustomPainterJob::renderPrepared qgsmaprenderercustompainterjob.cpp:233
QgsMapRendererTask::run qgsmaprenderertask.cpp:240
QgsTask::start qgstaskmanager.cpp:78
QThreadPoolPrivate::reset :
QThread::start :
BaseThreadInitThunk :
RtlUserThreadStart :
QGIS Info
QGIS Version: 3.10.14-A Coru�a
QGIS code revision: 8374282d2a
Compiled against Qt: 5.11.2
Running against Qt: 5.11.2
Compiled against GDAL: 3.1.4
Running against GDAL: 3.1.4
System Info
CPU Type: x86_64
Kernel Type: winnt
Kernel Version: 10.0.19042

The size of the image is 24000 * 20000 pixels at 1:250k resolution and 254 dpi (I found that 254 worked for a 16000 * 12000 file in terms of getting it to fit the canvas properly, but I don't know if it works for this size).
I have 32 GB of RAM, so it should be fine, right? Please bear in mind I'm new to GIS.


texts: Some objects missing after importing a .shp into PostGIS

texts:
I have a .shp file which I am trying to import into a PostGIS database using the "PostGis SHP and DBF Loader". This works smoothly but some (maybe 5%) of the records have an empty geometry field. (I have double checked the original files and there should be geometry on every row). I can't see anything special about these rows - it appears random but presumably isn't.

Can anyone tell me why this might happen?


texts: Editing Collector Field So That It Displays All Options Without Drop-Down List

texts:
I'm attempting to make a Collector form so that all options for a new feature are visible without having to click on it and open a drop down list.
The below image is what I want it to look like:

However, the below images show what I have now:

Once you click on an option, for example, "EM new feature", it opens to this:

And then, once the "Equity Elements" option is chosen, opens to:

My goal is to have everything visible as in the first image, without having to go through the extra steps to get to the desired option. Though, nothing I've found online or tried has led me to being able to obtain this in the collector app.
Does anyone know what I'm missing or doing incorrectly?


texts: Whitebox GAT not starting under Windows 10 - General JAVA or Whitebox problem?

texts:
I am currently running the DEM preprocessing on a relatively large (50000 km2 at 30 m resolution) and steep DEM. The objective is to extract the stream network and some morphologic attributes (e.g., slope) of the extracted streams. the problem is that the DEM is quite noisy on the valley bottoms, hence, if I use the standard fill tool in ArcGIS, a majority of valley bottoms is filled, and the information that was available there originally, is lost.

So much for the context. I was using WhiteBox GAT for similar problems in the past, because its breaching Algorithm was great for that objective.
Now, I was forced to update my computer to Windows 10, recently. Now I am not able to start the WhiteboxGIS.jar, any more. I have the newest Java version (8.51) from Oracle installed. The problem persist both by clicking on the jar file, or by opening it from the console. Using either method, just nothing happens: no error, no messages in the console, nothing.

In the past I had a similar problem with Whitebox (under Win 7) that was solved by repairing the association between java and the jar files using this tool: http://johann.loefflmann.net/de/software/jarfix/index.html. Now, unfortunately, that method is not successful. I also tried the standards, e.g., reinstalling java.

Unfortunately I also do not have access to another PC with a different operating system to run it there. I am also not an expert on Java, so unziping the jar did not help me understand potential causes for the problem.

Did anybody already trying to run Whitebox under Win 10, or has anybody who is more familiar with JAVA any idea for how I could get this work?


texts: How to select geography features within a certain distance given in miles?

texts:
I am just taking the first steps in using POSTGIS, but having trouble finding a simple example.

I have a database with a Key field, lat and lon fields and a GEOGRAPHY field.

I just want to have a query that gives me all of the Key fields with 100 miles of the following Lat/Lon 52.474702, -132.373041.

A pseudo query would be something like

select Key from GIS where GEOM is within 100 miles of 52.474702, -132.373041


texts: Qgis2web - how to make 'hybrid' satellite layer be able to switch on/off after exporting?

texts:
As the title indicate, How to make 'hybrid' satellite layer be able to switch on/off after exporting?
All other layers have blue ticked box to switch on/off.
Only hybrid layer does not have it.
Is there any setting that I missed out?


texts:
How do you compute vertical transformations in QGIS?

Ideally, the procedure / plugin will allow the loading of sundry geoid models.


texts: NAD 1983 Illinois State Plane- not projecting correctly

texts:
I'm having trouble with CAD drawings not projecting correctly within ArcMap 10.2.  I am working in NAD 1983 State Plane Illinois East FIPS 1201 (US Feet) and want my CAD drawing to be located in the city of Chicago near my other layers.

Below are the three scenarios I've tried.

1.) I've opened ArcMap, added an existing shapefile whose coordinate system is NAD_1983_StatePlane_Illinois_East_FIPS_1201_Feet.  I checked my Data Frame propeties, and it is also this coordinate system.  I saved my map, closed it and opened ArcCatalog.  In ArcCatalog, I tried to "define projection" for my CAD drawing.  I set this to NAD_1983_StatePlane_Illinois_East_FIPS_1201_Feet.  However, when I import this layer into ArcMap, it is floating in Missouri. I right clicked, opened the properties and it reads "Projected Coordinate System: NAD_1983_StatePlane_Illinois_East_FIPS_1201_Feet"

2.) In AutoCAD Map 3d, I open the CAD drawing. I chose "Map Setup">"Coordinate System">"Assign".  I then navigated to USA, Illinois and chose "IL83-EF" with the decription "NAD 83 Illinois State Plane, East Zone, US Foot".  I save this and add the CAD file to my ArcMap document.  It is also floating in Missouri but reading "Projected Coordinate System: NAD_1983_StatePlane_Illinois_East_FIPS_1201_Feet".

3.) I had a technical call with ESRI, and the support tech advised me to do the following. In ArcMap, I add the CAD drawing (with no coordinate system applied) to the map, execute the "CAD to Geodatabase" Conversion tool. Within this dialogue box, the tech asked me to clikc on the "Spatial Reference" tab at the bottom.  I chose this, and was then instructed to search for "Ohio" and choose the north state plane zone. This function adds a point, polyline, polygon and annotation layer to my GDB.  However, if I right click and choose "Zoom to Layer", I am again brought to Missouri. However, there is no layer visible.

I am wondering if anyone has experience with CAD files not projecting correctly?   Is this normal, and will I just have to use the "shift" in the Georeference toolbar?  I am also very curious to know if th


texts: QGIS and SpatiaLite database: Setting user access permissions

texts:
SpatiaLite databases seems to match all my desires to store GIS layers and table.

But an important question is: is it possible to limit writing to only one user ? If yes, what is the way to do it? I read the cookbook and searched the internet without finding an answer.

I know that PostGIS is better for multi-access purposes but there are only three cartographers in my company and they don't have enough needs in GIS to switch to PostGIS, which is much more complicated than SpatiaLite...

Thank you.


texts: Seeking accurate transformation library to convert between EPSG 4326 and 3857 for C#?

texts:
I've used Proj.Net library to convert X, Y in EPSG 3857 to latitude and longitude coordinates in WGS84 and vice-versa, but the conversion of latitude is not accurate compared to http://spatialreference.org/ref/sr-org/6864/.

For example, if the input is:

5676207.469887, 4291510.258087 then outputs are:

50.990239260333, 36.116815741834 using spatialreference.org and
50.990239260328906 36.116603266233213 using Proj.Net

Below is my code:

private const string Wkt3857 = "PROJCS[\"Popular Visualisation CRS / Mercator\", " +
                                       "GEOGCS[\"Popular Visualisation CRS\", " +
                                       "DATUM[\"Popular Visualisation Datum\", " +
                                       "SPHEROID[\"Popular Visualisation Sphere\", 6378137, 0, AUTHORITY[\"EPSG\",\"7059\"]]," +
                                       "TOWGS84[0, 0, 0, 0, 0, 0, 0], AUTHORITY[\"EPSG\",\"6055\"]]," +
                                       "PRIMEM[\"Greenwich\", 0, AUTHORITY[\"EPSG\", \"8901\"]], " +
                                       "UNIT[\"degree\", 0.0174532925199433, AUTHORITY[\"EPSG\", \"9102\"]], " +
                                       "AXIS[\"E\", EAST], AXIS[\"N\", NORTH], " +
                                       "AUTHORITY[\"EPSG\",\"4055\"]], " +
                                       "PROJECTION[\"Mercator\"], " +
                                       "PARAMETER[\"False_Easting\", 0], " +
                                       "PARAMETER[\"False_Northing\", 0], " +
                                       "PARAMETER[\"Central_Meridian\", 0], " +
                                       "PARAMETER[\"Latitude_of_origin\", 0], " +
                                       "UNIT[\"metre\", 1, AUTHORITY[\"EPSG\", \"9001\"]], " +
                                       "AXIS[\"East\", EAST], AXIS[\"North\", NORTH], " +
                                       "AUTHORITY[\"EPSG\",\"3785\"]]";

        public static void XyFromLatLng(double lat, double lng, out double x, out double y)
        {
            var gcsWgs84 = GeographicCoordinateSystem.WGS84;
            var cf = new CoordinateSystemFactory();
            var f = new CoordinateTransformationFactory();
            var sys3857 = cf.CreateFromWkt(Wkt3857);

            var transformTo3857 = f.CreateFromCoordinateSystems(gcsWgs84, sys3857);
            double[] fromPoint = { lng, lat };
            double[] output = transformTo3857.MathTransform.Transform(fromPoint);
            x = output[0];
            y = output[1];
        }

        public static void LatLngFromXY(double x, double y, out double lat, out double lng)
        {
            var gcsWgs84 = GeographicCoordinateSystem.WGS84;
            var cf = new CoordinateSystemFactory();
            var f = new CoordinateTransformationFactory();
            var sys3857 = cf.CreateFromWkt(Wkt3857);

            var transformToWgs84 = f.CreateFromCoordinateSystems(sys3857, gcsWgs84);
            double[] fromPoint = { x, y };
            double[] output = transformToWgs84.MathTransform.Transform(fromPoint);
            lat = output[1];
            lng = output[0];
        }

So I need a library in C# to transform these two CRSs to each other as accurate as possible.

I used DotSpatial for projection and the values are accurate compared to PostGIS. But I don't find out the reason of latitude difference required using Proj.Net.


texts:
I am trying to import a TIFF file that contains 2 bands. When using the following code in R, it seems only the first band is being recognised.

S1<-"my/path/"
S1<-list.files(S1, full.names = TRUE, pattern="tif$")
S1<-lapply(1:length(S1), function (x) {raster(S1[x])})

S1 being at this point:

> S1
[[1]]
class       : RasterLayer
band        : 1  (of  2  bands)
dimensions  : 3865, 6899, 26664635  (nrow, ncol, ncell)
resolution  : 14.83, 14.83  (x, y)
extent      : 361363.5, 463675.7, 5760647, 5817965  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=32 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
data source : /shared/Training/EARSEL0918_UrbanClassification_Germany/Processing/Coherence__20180412_20180424_TC.tif
names       : Coherence__20180412_20180424_TC

My objective is to create a raster stack containing both layers (for further processing in R). If I run

S1<-stack(S1)

> S1
class       : RasterStack
dimensions  : 3865, 6899, 26664635, 1  (nrow, ncol, ncell, nlayers)
resolution  : 14.83, 14.83  (x, y)
extent      : 361363.5, 463675.7, 5760647, 5817965  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=utm +zone=32 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0
names       : Coherence__20180412_20180424_TC


texts: How is range calculated by a GPS receiver?

texts:
I wanted some help on understanding how range is calculated by a GNSS/GPS receiver(rx), I tried looking on google, but can't find the required results.

C/A codes repeat every 1 ms and the rx delays internally generated sequence untill it correlates with the sequence it receives from the satellite. Thus the `delay' can only be between 0 to 1 milli second, which means the distance between satellites and rx can be [0, 300) km. Where as the actual distance is around 20,000 km.

What am i missing?


texts: Labeling with two decimal places (rounded) even when second is zero in ArcMap?

texts:
I have a benchmark data that need to label.

The data where in the field has attribute height like 12.340, 12.400, 13.690.

I have to label the height for examples 12.34, 12.40, 13.69.I already used vb script Round([ZOR],2).

However, for 12.34 and 13.69 it was okay but for 12.400 it just appear 12.4, i need to it appear 12.40.

So what should i do?


texts:
I have a polyline that I would like to query features that are on it's left, right, top, bottom. I am not looking for the features surrounding it but for features that are on a certain location from the polyline. Example: I am trying to query features that are on the left side of the polyline. The pictures below demonstrate what I am looking for.

I tried ItopologicalOporator, but the buffering takes all the features surrounding it which is not what I want.
any Ideas  ?


texts: Problems with labels while using Value relation

texts:
I'm using the value relation widget to populate fields from a .csv-table with values corresponding to a chosen value.

I don't know why, but the label caption displayed on the map is showing the key column and not the value which is related. Any ideas what I'm doing wrong here?

I'm using the current_value() filter.


texts: ESRI REST Api - need to download batches of tiles

texts:
So I wanted to use the ESRI REST Api to fetch map tiles for an application. Right now I am calling for each tile as I need it via its X/Y and zoom level(which is generating a bunch of http calls), but I'm wondering if there is a way to call for a batch of tiles in a single call.

I've been looking through the API documentation but hadn't found a command to get more than a single tile, I was wondering if anyone with more experience might know how.


texts: Counting number of roads between two points using QGIS

texts:
I use QGIS 2.18 and I have a point shapefile of many animal locations and a shapefile of roads.

I want to use QGIS to tell me how many roads are between each pair of animal locations.

So, basically how many roads would they have to cross to get to the other point.


texts: Assigning circle colours to map using QGIS?

texts:
I am trying to make a map showing which zip codes are served by which hospitals. There are five hospitals and each can serve one or more zip codes. I'd like to show each zip code with a dots of different colors to represent the different hospitals serving that zip.

I'm trying to assign colour circles to a map with zip codes. Each zip code has hospitals assigned to it and I would like the hospitals to be represented as colour circles on the assigned zip codes. So, for instance if a zip code has three hospitals I would like the hospitals to be represented by three different colour circles each.

The hospitals do not have addresses since they are assigned to zip codes in which they do not reside. So, I used the zip codes to get the Latitude and Longitudes. In most cases the hospitals have similar latitude and Longitude because they were assigned primarily by zip code.

I've tried layer properties and categorized symbol but it just assigns one color circle per zip code. Do you know what I'm doing wrong?


texts: How to define pour points from a stream network at a barrier for watershed calculations?

texts:
I'm looking for a way to check if cells from one raster file have contact or get close to cells from a second raster file. For example the blue raster is a barrier and the coloured is the flow accumulation and I want to know where it is possible for the flow accumulation to pass over in to the barrier raster which could be a river. Afterwards I want to define the watershed from this over passing points. I don't know how to find this points to compute with them, so I'm just looking for some suggestions how to start not a complete solution.

Using the con tool like Hornbydd suggested

Flowaccumulation raster

Extracted Flowaccumulation raster

Stream Network


texts: How to Display the World Map in Projected Coordinate Systems

texts:
I am trying to create a simple presentation map and I download the global level country boundaries shapefile from http://www.diva-gis.org/Data. The layer coordinate system is in Geographic Coordinate System: GCS_WGS_1984. Can you please let me know how I can display the map on UTM or any projected coordinate System using ArcMap?

Thanks


texts: Oracle spatial SDO_RELATE: why better performance result from UNION ALL/INTERSECT combined individual specified mask

texts:
Initially I was trying to find out why it's so slow to do a spatial query with multiple SDO_REALTE  in a single SELECT statement like this one:

SELECT * FROM geom_table a
WHERE SDO_RELATE(a.geom_column, SDO_GEOMETRY(...), 'mask=inside')='TRUE' AND
SDO_RELATE(a.geom_column, SDO_GEOMETRY(...), 'mask=anyinteract')='TRUE';

Note the two SDO_GEOMETRY may not be necessary the same. So it's a bit different from SDO_GEOMETRY(a.geom_column, the_same_geometry, 'mask=inside+anyinteract')='TRUE'

Then I found this paragraph from oracle documentation for SDO_RELATE:

  Although multiple masks can be combined using the logical Boolean
  operator OR, for example, 'mask=touch+coveredby', better performance
  may result if the spatial query specifies each mask individually and
  uses the UNION ALL syntax to combine the results. This is due to
  internal optimizations that Spatial can apply under certain conditions
  when masks are specified singly rather than grouped within the same
  SDO_RELATE operator call. (There are two exceptions, inside+coveredby
  and contains+covers, where the combination performs better than the
  UNION ALL alternative.) For example, consider the following query using the logical
  Boolean operator OR to group multiple masks:

SELECT a.gid   FROM polygons a, query_polys B   WHERE B.gid = 1   AND
SDO_RELATE(A.Geometry, B.Geometry,
                   'mask=touch+coveredby') = 'TRUE';

  The preceding query may result in better performance if it is
  expressed as follows, using UNION ALL to combine results of multiple
  SDO_RELATE operator calls, each with a single mask:

SELECT a.gid
      FROM polygons a, query_polys B
      WHERE B.gid = 1
      AND SDO_RELATE(A.Geometry, B.Geometry,
                   'mask=touch') = 'TRUE' UNION ALL SELECT a.gid
      FROM polygons a, query_polys B
      WHERE B.gid = 1
      AND SDO_RELATE(A.Geometry, B.Geometry,
                   'mask=coveredby') = 'TRUE';

It somehow gives the answer for my question, but still it only says: "due to internal optimizations that Spatial can apply under certain conditions". So I have two questions:

What does it mean with "internal optimization", is it something to do with spatial index? (I'm not sure if I'm too demanding on this question, maybe only developers in oracle know about it.)
The oracle documentation doesn't say anything about my original problem, i.e. SDO_RELATE(..., 'mask=inside') AND SDO_RELATE(..., 'maks=anyinteract') in a single SELECT. Why does it also have very bad performance? Does it work similarly to SDO_RELATE(..., 'mask=inside+anyinteract')?

EDIT:

After some research, I've got some additional information. When explaining the execute plan for ANDed SDO_RELATE case, it shows that only one of predicate has been applied with spatial index, the image of the plan is enclosed.

One way to improve the performance I just found is to put one of the SDO_RELATE into sub-query, and that will make both predicate utilize the spatial index. That is:

SELECT * FROM geom_table a
WHERE SDO_RELATE(a.geom_column, SDO_GEOMETRY(...), 'mask=inside')='TRUE'
AND a.id IN (SELECT b.id FROM geom_table b WHERE SDO_RELATE(b.geom_column, SDO_GEOMETRY(...), 'mask=anyinteract')='TRUE');

And the execute plan is:

But I wonder if there is better options, because here the HASH JOIN seems to be an extra cost.


texts: QGIS 3.10 create DEM from height contours - the maximum height values of the grid to do not match the maximum height values of the contours

texts: I have created a DEM using contours spaced 10 metres via QGIS grid ITN interpolation method based on height (metres) cell size 10 pixels/metres in QGIS, and another grid of 5 pixels cell size. When I load the DEM into QGIS, the maximum height is below that of the highest contour.  Would I have to adjust the grid cell size or other parameters to ensure a closer match with the maximum height of the contours? Or use another interpolation method eg nearest neighbour or IDW?


texts: ExecuteError: ERROR 000358 - Using a variable within arcpy.Select_analysis()

texts:
I have looked at Including variable in where clause of arcpy.Select_analysis()? and while it is very similar to the question I'm asking I ran into the same issue after implementing its syntax. I feel like I'm still missing something related to how to implement the correct syntax for the below.

where_query = '"field_name" <= ' "{0}.format(year_range)"

I'm trying to make what into a geoprocessing tool that can be reused in ArcMap, but I'm receiving the following error. It looks like an issue with the way that I'm formatting the where_query variable, but I can't seem to get it right.

Traceback (most recent call last):
  File "C:\ArcPy\mean_center_drift.py", line 19, in <module>
    arcpy.Select_analysis(in_feature, year_out_name, where_query)
  File "c:\program files (x86)\arcgis\desktop10.3\arcpy\arcpy\analysis.py", line 84, in Select
    raise e
ExecuteError: ERROR 000358: Invalid expression 1
Failed to execute (Select).

and here is the code in its current form.

import os
import arcpy

in_feature  = arcpy.GetParameterAsText(0)
out_features  = arcpy.GetParameterAsText(1)
origin_year = arcpy.GetParameterAsText(2)
field_name = arcpy.GetParameterAsText(3)

for x in range(10, 140, 10):
    year_range = int(origin_year) + x
    where_query = '"field_name" <= ' "{0}.format(year_range)"

    year_out_name = os.path.join(out_features, "Years_{0}".format(x))
    mean_out_name = os.path.join(out_features, "Mean_{0}".format(x))

    arcpy.Select_analysis(in_feature, year_out_name, where_query)

    arcpy.MeanCenter_stats(year_out_name, mean_out_name, "#", "#", "#")


texts:
my problem is regarding routing.
I have a network based on bicycle graph links (SRID: 4326), with bidirected and directed ways.
I want to give the user the functionality to draw a route (shortest path) from a given point to a certain location (the input parameters are sent from my app, as two geometric POINTS).

So to approach this problem i did the following steps:

Added reverse cost column and assigned an arbitrary high value for the ONE WAY links (Geometry LINESTRING)
Used the ST_ClosestPoint function to get the nearest source node from the start position, and the nearest target node from the final destination.
Once i have my source node and target_node i launched the Dijkstra Shortest Path function of pgrouting (taking account the reverse_cost for ONE WAY links):

SELECT bz_bike_graph_line.* FROM bz_bike_graph_line JOIN
    (SELECT * FROM shortest_path('SELECT gid as id, start_node::int4 as
                     source, end_node::int4 as target, distance::double
                     precision AS cost, reverse_cost::double precision
                     AS reverse_cost FROM bz_bike_graph_line' ,source_node,
                     target_node, true, true)) AS route
     ON bz_bike_graph_line.gid = route.edge_id

QUESTION:

All seems to be working in the right way BUT when returning the route the starting point and the starting route node never coincides (idem target and final dest.), obviously i should say, due the fact we return the closest points.
Is there a way to resolve it or solve the task in other way?
Some GIS expert here :) ?

P.A.
Thanks!


texts:
I have two GPX-Vector-Layers with points. One "Todo-File" (Locations that I want to visit) and one "Done-File" (locations that I had visited extract from the GPS-Data from phote I had taken).

I have these two layers in QGIS and what I want is a seperate layer, that includes all points, that dont have a (fuzzy) match (=all Points that I really dont visited).

Is there a chance to get this? In best case with a tolerance (Done-Point is/is not within a radius of 100m from a Todo-Point)


texts: Erase tool for polylines with ArcGIS Standard License

texts:
I have a very long polyline representing a pipeline.  I also have another feature class with smaller polylines which are coincident with the pipeline.  I am trying to cut the pipeline with the smaller segments to essentially get the inverse of the feature class of smaller segments.

So, I have this:

and I'm trying to get this:

I cannot use the Erase tool because I am using the Standard License of ArcGIS desktop.  I cannot install third-party software (ET-Geowizards) or QGIS due to IT policies.

Normally, I would use the Union tool but it only works for polygons.

One thought I had was to use linear referencing and create an event table of the gaps and generate line segments from that.  Another thought was to use search and insert cursors to build the lines based on the coordinates of the existing lines.  But I was hoping to use some workflow of the analysis/overlay tools. Surely, this has been done before using built-in tools.

I don't want to do this manually because I have hundreds of segments.

Any thoughts on a workaround?


texts:
How can we make 2 lines that meet each other (sitting on top of another), intersect each other in QGIS? How can we make a point between them?

The reason why I want to make an intersection is that I want to make routing work. (I want the route to go through that point instead of taking longer rotues.)


texts: Leafet - Loading and Instantiating WFS data on Demand

texts:
I am designing a WebGIS that will give the user the option to load a large range of WFS layers from geoserver. I am keen to avoid the overhead of having to load and instantiate all of these WFS datasets on pageload and would rather only load them when they are toggled on in the layer control.
Will this potentially avoid downloading MB's of data that may not even get added to the map by the user?
I figured that If I could attach a function to the layers toggle to either instantiate (if it doesn't already exist) the layer using the WFS feed or remove it from the map depending on whether it is being toggled on or off.


texts: Converting sign for a 32-bit GeoTIFF with an open source tool?

texts:
I have a 32-bit GeoTIFF that represents floating point elevation values, and I'd like to change the sign so that the values are negative instead of positive.

What is the easiest way to do this using open source tools?


texts:
I'm trying to use a map of Austrias districts from the GADM project.
The Shapefile can be found here: AUT_adm.zip (All other formats provided by the project appear exectly equal.)

  The coordinate reference system is latitude/longitude and the WGS84 datum.

Comparing this map with the official (but not as detailed) NUTS3 shapefiles from "Statistics Austria" and with Google Maps (all done within QGis) shows differences between the borders of the administrative areas.

Google Maps (and OpenStreetMap, ...) are nearly equivalent with the files from Statistics Austria.
Unfortunately, the maps from GADM are not congruent with the other sources. Additionally the distortion grows from West to East.

West of Austria (GADM=red, StatA=green, Google maps as reference):

East of Austria (GADM=red, StatA=green, Google maps as reference):

Of course I tried to import the original data as well as projections to the same CRS (epsg:3857). As expected this did not alter the results.

My questions:

Am I using a wrong projection for the GADM data? (Is the documentation on the homepage and the files not corret.)
How can I get the right projection and get rid of the differences?


texts: How to get extent of multiple geometries?

texts:
I can query a geometry extent like this:

select ST_Extent(geom) from tableName;

But sometimes I am querying geometries more than one. I am getting 5 points that different locations. So I wanto to zoom my map that containing this 5 points extent. Is this possible?

select geom from tableName where category='airport';

This gives me 5 airport. and I want to zoom this 5 point extent.


texts: Register as versioned with move edits to base - when not to?

texts: What limitations do I get, when I chose to register as versioned with the option to move edits to base?


texts:
I am trying to rasterize a spatialite database using gdal_rasterize. I am using gdalUtils in R to do so and have run the following

gdal_rasterize -burn 1 -sql SELECT * FROM plants WHERE binomial="Albies alba" -ts -dialect SQLITE 129600 64800 -ot Byte ./plants.sqlite ./test_r.tif

This is giving me a status 1 error and I am wondering if I my SQL syntax is wrong.


texts: FME Error when writing to relationship class in ArcSDE Geodatabase writer type

texts:
I'm working in the FME workbench through ESRI Data Interoperability extension (ArcGIS 10.4.1). I'm writing to an SDE geodatabase with relationship classes established. When attempting to setup the relationship between a table data stream and feature classes in FME, I get the following error:

"Unexpected execution path. Please contact Safe Software
GeodatabaseRelationshipFeaturesPipeline::GeodatabaseRelationshipFeatures(SortFactory): A fatal error has occurred. Check the logfile above for details.
Some error has occurred closing the pipeline factory, see previous warning
None of the edits made during the edit session were saved."

I'm fairly certain that I have all feature type parameters set correctly to write feature relationships from the source fgdb fc's and sql server non-spatial table data to the destination relationship class. The relationship specific attributes set include:

-geodb_feature_has_relationships = Yes
-geodb_oid - set on each dataset from counter transformers
-geodb_type - set to geodb_relationship
-geodb_rel_origin_oid - renamed from fc geodb_oid
-geodb_rel_destination_oid - renamed from table data geodb_oid

A feature merger transformer was used to build relationship (O:M), then the merged output port was connected to the relationship class writer.

The translation was working prior to setting up the relationship (when connecting the source fc data streams to the writer polygon fc, and the source non-spatial table data stream to the writer table.  Now the translation fails even after disabling the relationship workflow portion of the workspace.  The same error is thrown every time.


texts:
I'm trying to calculate the amount of pixels with the value of 0.3 to 0.7 between two different rasters (each raster of a different year) to calculate the loss of green spaces trough time.
Is it possible to do this in QGIS or SAGA GIS?


texts: Add custom maps in Mapbox as tile layer in Folium

texts:
I was trying to create a map in Folium using a custom map built tile layer created with Mapbox Studio (using 3 different datasets and a satellite basemap). I was wondering if there is a simple way to add it with Folium? I read in the documentation that style layers and datasets had to be called separately. I tried with the following (and a few variations), but each time it loaded the Python data with no base map.

folium_map = folium.Map(location=[-8.907970, 33.433200],
                        zoom_start=7,
                        tiles='http://{s}.tiles.mapbox.com/v4/(username).(mapid)/{z}/{x}/{y}@2x.png?access_token=(access token)',
                        attr='Mapbox')

I was using Mapbox (instead of loading the datasets independently) because my datasets are pretty big (up to 100MB, upwards of 100,000 points) and that causes the Leaflet map to lag tremendously.


texts: How to doing spatial join for several rasters

texts:
I have four rasters, each of 32 bits (single band), with the same cell size, the same number of rows and columns and identical extension. I need to make something like a "Spatial join", I mean, I need one of these four raster contains in its attribute table cell values of the other three raster.

I am programming this with Python and I've tried several tools, including "Combine" that generates me a raster with the combination of the four raster and its attribute table contains columns with the values of each raster, however, I want the cell value is not altered and hold the raster MLWS, only adding to its attribute table three columns with the values of the other three raster (HAT , LAT , MHWS ). For this I need to generate some kind of spatial Join, looking for the corresponding cell in the other raster and bring value to my raster attribute table.
This is an image for the obtained with the "Combine" tool:

It works perfectly, I just need the final raster represents only MLWS values and not the combination of the four rasters, containing in the attribute table columns with the values of the other three raster.


texts: Why does my mouse cursor jump to another screen location when I try to edit a Layer?

texts:
I'm using QGIS 3.18 Zurich installed with OSGeo4W64 installer on an Dell Latitude 5290 2-in-1 tablet running Windows 10 and I found a consistent bug. This bug does not occur on my Windows 8.1 desktop running the same version of QGIS.
Problem:
The mouse jumps to another screen location when I attempt to click the Toggle Editing "pencil" icon. That will cause a bug forcing the mouse cursor to jump to the bottom left hand corner of the status bar.
I'm also unable to select the Check QGIS version drop down menu item under Help when I attempt to hover the mouse over it. That will also cause the mouse cursor to jump away.
This bug happens even when there are no layers loaded yet, but moving the cursor to the Toggle Editing bug still occurs. I've tried different version of QGIS stand-alone, Hanover, Zurich, and OSGeo4W64 installer versions. Same problem.
Computer Details--
Tablet: Dell Latitude 5290 2-in-1,
Processor:  Intel(R) Core(TM) i5-8350U CPU @ 1.70GHz   1.90 GHz,
Installed RAM:  8.00 GB (7.86 GB usable),
System type:    64-bit operating system, x64-based processor,
Pen and touch:  Pen and touch support with 10 touch points
Windows Specification--
Edition:    Windows 10 Pro,
Version:    20H2,
Installed on:   3/28/2021,
OS build:   19042.985,
Experience: Windows Feature Experience Pack 120.2212.2020.0


texts: What could be making QGIS 2.10 crash when loading a project file?

texts:
On QGIS 2.10, running on Ubuntu 15.04 Linux, I have a project file that has a fairly large number of mainly shapefile layers, which gets to 99% on the slider to the right of where it has a message saying 'Loading layer name of last layer' and then just hangs there, for a long period of time (sometimes an hour) before the qgis process crashes.

I haven't had this fault before, even when working with other projects with large numbers of complex layers, and this fault affects one particular project file for no obvious reason.

I have already tried a couple of different methods of resolving this, firstly by removing this last layer by manually editing the .qgs project file in a text editor, and secondly by disabling all plugins and restarting qgis in case it could be a plugin causing the problem. Neither of these affected the issue. If I edit the .qgs file to remove the last layer, the program hangs instead as it displays 'Loading layer...' for the preceding layer.

I have no reason to think that the data of any of the actual layers are missing or corrupted, and in any case if this had happened I would have expected the "Handle bad layers" message instead.

Does anyone know what kind of situation is likely to have caused this, or how to remedy it?


texts: What is the recommended method to install or upgrade to QGIS 3 on centos 7/8?

texts:
On https://www.qgis.org/en/site/forusers/alldownloads.html#rhel-centos-scientific-linux
There is reference to use http://elgis.argeo.org/ for RHEL/CENTOS... but this is only about RHEL/CentOS/SL 6 and 5
In epel7 the current qgis version is 2.14.9, and it does not appear available in epel8
I ask this because a user complains that a recent GeoPackage file can apparently not be opened with qgis 2.14.9
Thanks in advance.


texts: Lat/long values in a WMS GetFeatureInfo request

texts:
I have tried to send manual WMS GetFeatureInfo request to GeoServer with the parameter values X and Y (pixel values) taken from the map object of OpenLayers.

I have got the response with the wrong values which is not matched with the values I have seen in the GeoServer preview with the same long/lat.

The requested URL was:

  geoserverUrl/wms?REQUEST=GetFeatureInfo&EXCEPTIONS=application/vnd.ogc.se_xml&BBOX=getextentValue&SERVICE=WMS&INFO_FORMAT=text/html&QUERY_LAYERS=xxxxx:atm&Layers=xxxxx:atm&WIDTH=mapwidth&HEIGHT=mapheight&format=image/jpeg&styles=&srs=EPSG:4326&version=1.1.1&x=eventxvalue&y=eventyvalue";

I have referred with these posts also: Can Geoserver return the raster value of a lat/lon point.

Is there any possible way to give the clicked long/lat values given instead of X and Y pixel values in a GeoServer WMS GetFeatureInfo request?

or

What are the factors which are to be clear when making this request manually. More particular width, height and the BBOX values?

In client side I have made an AJAX call:

map.events.register("click",this, function(e){

    var clickedLatLong =map.getLonLatFromViewPortPx(e.xy);
    var clickedPixel =  map.getPixelFromLonLat(clickedLatLong);
    //alert(clickedLatLong);
    //alert(map.size);
    if(!map.getMaxExtent().containsLonLat(clickedLatLong)) {
        return;
    }else{

        alert(clickedPixel);
        var bbox = map.getExtent();
        alert(bbox);
        $.ajax({
            type: 'GET',
            url: 'getMultiLayersValue.html',
            dataType: 'html',
            data:{long:clickedPixel.e.x,lat:clickedPixel.e.y},

            success : function(result){
                alert(result);
                $("#DisplayResponse").html(result);
                $("#DisplayResponse").dialog( "open" );

            },
            complete:function(xhr, txtStatus){

            },
            error: function(xhr, error){
                showLoadingMaskDialog(false,loadingText);
                alert("Taking More Time to Respond\n Please Try Later or Contact System Administrator");
            }
        });//ajax
}

In server side I have made this URL request:

  "localhost:8080/geoserver/wms?REQUEST=GetFeatureInfo&EXCEPTIONS=application/vnd.ogc.se_xml&BBOX=getextentValue&SERVICE=WMS&INFO_FORMAT=text/html&QUERY_LAYERS=xxxxx:atm&Layers=xxxxx:atm&WIDTH=mapwidth&HEIGHT=mapheight&format=image/jpeg&styles=&srs=EPSG:4326&version=1.1.1&x=eventxvalue&y=eventyvalue";

Eventually no error messages were found only I got the wrong values for my clicked long/lat. May be the values of pixel getting wrong but I unable to find the solutions.


texts: Create a table that is not registered with the geodatabase

texts:
I am doing some testing that requires having a Oracle table that is not registered with the geodatabase. For instance, Adding an Incrementing ID Field can only be performed on a table that is not registered with the geodatabase.

How can I create a table this not registered with the geodatabase?


texts: Positioning esri.symbol.PictureMarkerSymbol at the center when clicked?

texts:
I use ArcGis Javascript api 3.2. I have a symbol:

candidateSymbol=new esri.symbol.PictureMarkerSymbol('image/PushPin.png', 30, 30);

I want to center it when it is being clicked.

PS: There will be many such symbols on the map. I want to center the screen on the symbol that is being clicked.


texts: Inserting a new field to a memory layer python GDAL

texts:
i'm having an issue with creating a new field in a layer create on a memory data source and adding values to it : here is my snippet code .

newfield = ogr.FieldDefn("aa", ogr.OFTReal)
erasedlayer.ResetReading()
erasedlayer.CreateField(newfield)
for i in range(erasedlayer.GetLayerDefn().GetFieldCount()):  # iterate on parcelle features to get the names
    print(erasedlayer.GetLayerDefn().GetFieldDefn(i).GetName())
erasedlayer.ResetReading()
for feature in erasedlayer:
    feature.SetField("aa", feature.GetGeometryRef().GetArea())

and this is the error i got :

Traceback (most recent call last):
  File "/home/kail/PycharmProjects/project/file.py", line 190, in <module>
    feature.SetField("aa", feature.GetGeometryRef().GetArea())
  File "/usr/lib/python3/dist-packages/osgeo/ogr.py", line 4808, in SetField
    return _ogr.Feature_SetField(self, *args)
RuntimeError: Invalid field name: 'aa'

i verified with the loop for the field names on my layer and the "aa" field exist


texts: The select-by-location tool of ArcMap won't work properly!

texts:
I have faced a weird problem caused by the select-by-location tool of ArcGIS. My goal is to identify the polygons intersecting a reference (colored polygon). Accordingly, I applied the select-by-location tool, but it wrongly included some polygons that there was not any overlap between them and the reference polygon (as the following image).

The interesting thing is that when I use the same function in QGIS, I don't encounter such a problem.

What is going wrong there? Why does the ArcMap fail in correctly determining the intersecting polygons?


texts: Adding GeoJSON API live feed as a layer in my AGOL web map

texts:
I've only ever added static CSV and GeoJSON content. How can I add one of the following links to my AGOL map as a dynamically updated live feed point feature layer?

All Incident Data (.csv format)
https://www.fire.ca.gov/imapdata/mapdataall.csv
JSON API (Optional parameters: year, inactive=true/false)
https://www.fire.ca.gov/umbraco/api/IncidentApi/List?inactive=true
GeoJSON API (Optional parameters: year, inactive=true/false)
https://www.fire.ca.gov/umbraco/api/IncidentApi/GeoJsonList?inactive=true


texts: ArcGIS Server Geoprocessing service for saving data

texts:
I am using ArcGIS 10.2.2 Pack with ArcGIS Server.
I have two polygon feature layers and I created a model to intersect them and save output to new feature file in database. The model works then I press start.
Is it possible to create a geoprocessing service in ArcGIS server, which in given time interval takes two feature layers, intersect them and creates output?
I tried to make it, but it seems that service does not do any work, just stands on.


texts: Translate (convert format): Set the size of the output file

texts:
I have a raster .asc file with Pixel Size = 0.000863923,-0.000863923 (as given in Rightclick -> Properties -> Metadata -> Properties -> Pixel Size) and Dimensions = X: 5788 Y: 2315 (where X is width and Y is Height).

I want to export this raster (with a new name) after changing the properties into Pixel Size = 0.000833333,-0.000833333 and Dimensions = X: 6000 Y: 2400.

For this, I want to use Translate (convert format) in the Processing Toolbox -> GDAL/OGR -> [GDAL] Conversion.

In Translate (convert format), I am able to see an option :

Set the size of the output file (In pixels or %)

What value should I give in order to get a raster with Pixel Size = 0.000833333,-0.000833333 and Dimensions = X: 6000 Y: 2400 (from a raster with the same extents; but with Pixel Size = 0.000863923,-0.000863923 and Dimensions = X: 5788 Y: 2315) ?

Also, I am able to see another option :

Output size is a percentage of input size

Do I have to select (select by ticking or by giving 'Yes') this option or not ?

If I give Set the size of the output file (In pixels) = 100, then the python code will be :

processing.runalg("gdalogr:translate","C:/Users/Sreeraj/Desktop/1.asc",100,False,"",0,"","96.993358471,98.9933584709,14.0067404817,17.0067404815",False,5,4,75,6,1,False,0,False,"","C:/Users/Sreeraj/Desktop/2.asc")

You can see the value 100 in the above line of code (after the input file 1.asc).

If I want to give -outsize 6000 2400 ; then how can I change the given below python code (change from Set the size of the output file (In pixels) = 100 to -outsize 6000 2400) ?

processing.runalg("gdalogr:translate","C:/Users/Sreeraj/Desktop/1.asc",100,False,"",0,"","96.993358471,98.9933584709,14.0067404817,17.0067404815",False,5,4,75,6,1,False,0,False,"","C:/Users/Sreeraj/Desktop/2.asc")


texts: Creating holes from overlapping polygons in ArcObjects and C#?

texts:
I have a polygon layer which contains some overlapping polygons.

I want Holes instead of overlapping polygons.

i tried some code

                  ESRI.ArcGIS.ArcMapUI.IMxDocument mxd = ArcMap.Application.Document
         as ESRI.ArcGIS.ArcMapUI.IMxDocument;
            if (mxd.SelectedItem is ESRI.ArcGIS.Carto.IFeatureLayer)
            {
                ESRI.ArcGIS.Carto.IFeatureLayer flay = mxd.SelectedItem as ESRI.ArcGIS.Carto.IFeatureLayer;
                ESRI.ArcGIS.Geodatabase.IFeatureClass fcls = flay.FeatureClass;
                ESRI.ArcGIS.Geodatabase.ISpatialFilter spatialFilter = new ESRI.ArcGIS.Geodatabase.SpatialFilterClass();

                spatialFilter.SpatialRel = ESRI.ArcGIS.Geodatabase.esriSpatialRelEnum.esriSpatialRelOverlaps;
                ESRI.ArcGIS.Geodatabase.IFeatureCursor fcur = flay.Search(spatialFilter, true);
                ESRI.ArcGIS.Geodatabase.IFeature fea = fcur.NextFeature();
                while (fea!=null)
                {
                    ESRI.ArcGIS.Geometry.IGeometry geo= fea.Shape;
                    ESRI.ArcGIS.Geometry.IGeometryCollection gcoll = geo as ESRI.ArcGIS.Geometry.IGeometryCollection;

                    fea = fcur.NextFeature();
                }

How can I do this?


texts: How to fix error: Computed value too large, when performing supervised classification using Google Earth Engine

texts:
I am trying to run a supervised object-based image analysis on a Sentinel 2 image using a 13 class shapefile to train the classifier, but I am having trouble running the classification part of it and get the error: computed value to large.  I am using the following code (based on https://developers.google.com/earth-engine/classification)

    Map.addLayer(training, {color: 'red'}, 'training');

  // Get the values for all pixels in each polygon in the training.
var trainedpolygons = image.sampleRegions({
  // Get the sample from the polygons FeatureCollection.
  collection: training,
  // Keep this list of properties from the polygons.
  properties: ['classvalue'],
  // Set the scale to get Sentinel pixels in the polygons.
  scale: 10
});

// Create an SVM classifier with custom parameters.
var classifier = ee.Classifier.svm({
  kernelType: 'RBF',
  gamma: 0.5,
  cost: 10
});

// Train the classifier.
var trained = classifier.train(trainedpolygons, 'classvalue');

 // Classify the image.
 var classified = image.classify(trained);

// Display the classification result and the input image.
Map.addLayer(image, {bands: ['B4', 'B3', 'B2'], max: 0.5, gamma: 2});
Map.addLayer(classified,
            {min: 0, max: 1, palette: ['red', 'green']});`

Full code can be found here: https://code.earthengine.google.com/8f8a4ffeafd1e8c94869099583c53d4f


texts: Where can I find examples of fuzzy logic applications in ArcGIS?

texts:
The last ArcGIS version (10) includes two arctoolbox tools for fuzzy logic (Fuzzy Membership, and Fuzzy Overlay). Does anyone know where can I find some examples showing how they can be applied?

(From Desktop 10 help.)


texts: psycopg2 not working because column "geom" does not exist

texts:
I was planning to post this question on Stack Overflow before realizing that it probably was more related to my poor knowledge of GIS (I'm an absolute beginner in GIS and my knowledge of Python is moderate).

(config: pgAdmin 4, PostgreSQL 10)

I'm trying to use the solution provided in this repository but I have the following error:

  File
  "C:\Users\user\anaconda\Anaconda3\lib\site-packages\pyris\api\extract.py",
  line 47, in _query
      cu.execute(q, params) psycopg2.ProgrammingError: column "geom" does not exist

Where the 47th line is located within this piece of code:

def _query(q, params=None):
    """Carry out a SQL query

    Only fetch one result
    """
    Logger.debug("processing query '%s'", q)
    with psycopg2.connect(database="pyris",
                          user=DATABASE['USER']) as cnx:
        with cnx.cursor() as cu:
            if params is not None:
                cu.execute(q, params)
            else:
                cu.execute(q)
            return cu.fetchall()

I assume it comes from the fact that my table (created from this file) only has the following columns:

  gid, depcom, nom_com, iris, dcomiris, nom_iris, typ_iris, origine

I tried to manually add a "geom" column whose type would be "geometry" while the value would be "4326" but I'm not even sure I should use that number (as I said, I'm a noob).

Thus I first added a column "geom" with "geometry" type but the program is now throwing me this error:

  File
  "C:\Users\user\anaconda\Anaconda3\lib\site-packages\pyris\api\app.py",
  line 105, in get
      res.update(coord) AttributeError: 'list' object has no attribute 'update'

Located in here:

@api.route("/search/")
class IrisFromAddress(Resource):
    @api.doc(parser=address_parser,
             description="Look for an IRIS for a specific address.")
    @api.marshal_with(address_fields, envelope='address')
    def get(self):
        args = address_parser.parse_args()
        query = args['q']
        Logger.info("Look for IRIS for address '%s'", address)
        coord = address.coordinate(query)
        Logger.info("Get coordinate (%s, %s)", coord["lon"], coord["lat"])
        Logger.info("For address '%s'", coord["address"])
        if coord['address'] is None:
            return []
        res = extract.iris_from_coordinate(coord['lon'], coord['lat'])
        res.update(coord)
        return res

And when I try on pgAdmin to add a value to that column, I always get the same error no matter what I try:

ALTER TABLE geoiris
set geom = 4326

ALTER TABLE geoiris
set geom = '4326'

ALTER TABLE geoiris
set geom = cast(4326 as geometry)

Which all result in:

  ERROR:  syntax error at or near "geom"

Long story short, the more I try to find a solution, the more I get confused. Which is why I'm requesting your help.

EDIT:

I loaded my file through PostGIS Shapefile Import/Export Manager. And, in "Options...", I didn't select the "Load only attribute (dbf) data" option. Which -- I think -- is the right thing to do because it "strips the geometry column from the loading process, leaving just the attribute columns" according to this website.


texts: Changing ArcGIS Server 10.1 port (6080 to 8080)

texts:
Due some network issue (Maps are not displaying in secure environment), I want to change the default port 6080 to 8080. I am going through Web help and few links in which they mentioned that install ArcGIS web adopter but I do not think that it is recommend for Amazon instance.

We have our ArcGIS Server 10.1 on Amazon instance (windows environment) and I am using the REST link in my Java script applications & sometimes the ArcGIS Port 6080 is blocked due to network permissions and hence my maps are not visible that’s why I need 8080 port

Here are the links I am going through…

Ports used by ArcGIS Server
http://resources.arcgis.com/en/help/main/10.1/index.html#/Ports_used_by_ArcGIS_Server/015400000537000000/

Deployment scenarios
http://resources.arcgis.com/en/help/main/10.1/index.html#//015400000488000000

Amazon configuration for ArcGIS Server 10.1
http://hackgeo.com/arcgis/arcgis-server/arcgis-server-10-1/arcgis-server-10-1-on-ec2-single-server-initial-setup

So any workarounds/suggestions will be great : )

Thank you in advance!!


texts: Using .BMP image as customizable polygon fill in QGIS

texts:
Normally in my office, I work with ArcMap. For making vegetation maps, I use one specific fill all the time (surface1.bmp). This fill needs to be a different colour in each polygon/feature. In this crisis-times I work from home, using QGIS instead. The desired fill is not available in QGIS, so I imported it in the symbology as an SGV. However, I want to be able to customize the colouring of this fill.
How is this possible?
I imported it using Raster Image Fill and import the .bmp. This did not work, so i converted to SGV using inkscape. It looks fine, but I still cannot change colours for each feature as one could do with any other pre-cooked QGIS symbol style.
See below for the .bmp picture.

I want to be able to make the black and white different colours in QGIS, without having to do this manual before uploading it in QGIS.


texts: Is it possible to connect to FEMA's WMS and WFS using QGIS?

texts:
I have just found out about the wonderful world of WMS and WFS. I have added the FEMA versions of both to QGIS. WMS and WFS. The WMS seems to connect and adds something to QGIS, yet it is not visible. The WFS connects, but I get an error that it is not a valid layer.

So... I guess that my question is - Is anyone willing to test these and maybe give some guidance on how to better use these?

I have also tried NYS WMS and the Nationwide Inventory WMS http://107.20.228.18/ArcGIS/services/FWS_Wetlands_WMS/mapserver/wmsserver? with no luck. Similar results to FEMA's WMS.


texts: Raster calculator ERROR 000539 when converting values to be from 0-1 using ArcGIS Pro

texts:
I am getting this error when trying to make values of GHI from 0-1
I only add a layer then I run the raster calculator as attached.
This is my code :
 ("Global horizontal irradiation"  -  "Global horizontal irradiation".minimum) / ("Global horizontal irradiation".maximum - "Global horizontal irradiation".minimum)

ERROR 000539: Traceback (most recent call last): File "", line 1, in  File "", line 8, in rcexec_DC684390_0C54_4F14_9EBF_B34DF3F2E1B1 RuntimeError: Unspecified error Failed to execute (RasterCalculator).


texts: Adding a pop-up confirmation window before executing a script tool in a ArcGIS Toolbox

texts:
Similar question to this one, but it's 8 years old so seeing if theres any solutions since:
pop up message box before a tool execution with arcpy/python
I have a custom ArcGIS Toolbox, with a script tool. After I click "Ok" when running my tool, I want a pop-up window to appear, something along the lines of "Are you sure you want to execute this tool" -> Yes will run it normally, no will exit/return to the inputs of the script tool.
Is this even possible? I am on ArcGIS Desktop 10.6.


texts: R convert sf object back to SpatialPolygonsDataFrame

texts:
How do I convert an sf object back to a SpatialPolygonsDataFrame?

I am running R 3.4.0 and sf version 0.4.1. The sf documentation mentions an as_Spatial function, but this isn't available in my install (for some reason).


texts: Seeking a GIS for indoor mapping and direction finding

texts:
I've seen a few options for indoor GIS platforms such as Google's indoor mapping, implementations of OpenLayers, Open Street Map, and many other plausible solutions. There's even a post on here asking a very similar question to mine but all the answers involved how to provide locations to users.

For my use I don't care to map where the user is but rather simply provide a usable map of our building and in my ideal world I would be able to provide static directions (think Google Maps without location for our building). Is there anything out there that can help us achieve this? I'm willing to learn as long as there's documentation.


texts: Connecting different points (more than two) with different lines in QGIS

texts:
I have different points (stations) which I made as point layers in my map, and now I want to create a line network of a public transportation system, as such I want to connect different points (more then two, it might be even more than ten) to different lines.
I want the information of the line numbers which are crossing certain station in the proberties section of the certain stations.
On the other site, I want the information of the stations which are crossing by certain lines in the properties section of the lines.
What is the best tool for this?
EDIT: Yes, I know, I could easily make a line layer and just write the information manually in the chart , but that's not the way I want to do it respectively we are allowed too. Also, I have to use QGIS for that operation.


texts: Missing areas when using mosaic after projectRaster or vice versa

texts:
I have 290 MODIS NPP tif tiles that I want to merge and project into EPSG 4326. mosaic_raster() works just fine (Fig 1) but when I project the global raster (with or without a template) I lose almost all the southern areas (Fig 2).

Fig 1  Fig 2 .

Fig1 metadata

Fig2 metadata:crs correct, also this ext would be sufficient

The data has been acquired from here and the original crs is

+proj=sinu +lon_0=0 +x_0=0 +y_0=0 +a=6371007.181 +b=6371007.181 +units=m +no_defs

For merging I used mosaic_raster as discussed here

e <- extent(-180,180,-90,90)
epsg4326 <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
template <- raster(e)
projection(template) <- epsg4326
writeRaster(template,
            file="MyBigNastyRasty.tif",
            format="GTiff", overwrite=TRUE )
mosaic_rasters(gdalfile=my_rastlist,
               dst_dataset="MyBigNastyRasty.tif",      of="GTiff")

and  for reprojecting

npp_ <- raster("~/MyBigNastyRasty.tif")
npp_reproj <- projectRaster(npp, crs=epsg4326)

Same results if I use templateraster with correct res and crs inside the projectRaster.

I also tried to project every tiff tile first and then merging, but then I had a different origin. To tackle this, I tried to use the one global templateraster when projecting the tiles. However, then it took 25hours for 290 tiff files although I used 6 cores in beginCluster()!

For this, I used code

rlist <- lapply(my_rastlist, raster)
beginCluster()
npp_individual_reproj <- lapply(rlist, function(r) {projectRaster(r, templateraster)})
endCluster()

Any ideas what is wrong?


texts:
I created a CSV with GEOID's and imported it into ArcGIS Pro. Arc is treating the GEOID's as numeric and the files as read-only. When I attempt to export the file to a geodatabase, converting the GEOID to text as part of the export, I get Error 000210. As a result I cannot join this file with a tract map where Arc has treated the GEOID as text.
Any ideas?


texts: Image math (reduction) within image collection

texts:
I want to generate monthly average Sentinel-1 images in GEE and adapted this answer here to get 12 mean radar images over a year. However, as radar imagery in GEE is in decibel, this returns the geometric mean and not the arithmetic mean. I have been able to convert a single decibel image to power using the equation:

var PWR = ee.Image(10).pow(image.select('VH').divide(10))

but I do not know how to implement this over each image in a collection prior to calculating the monthly average. Is it possible to calculate this when generating an image collection?

My code is:

var sentinel1 = ee.ImageCollection("COPERNICUS/S1_GRD");

var months = ee.List.sequence(1, 12);

var byMonth = ee.ImageCollection.fromImages(
    months.map(function (m) {
      return sentinel1
      .filter(ee.Filter.eq('instrumentMode', 'IW'))
      .filterBounds(geometry)
      .filterDate('2017-01-01','2017-12-31')
      .filter(ee.Filter.calendarRange(m, m, 'month'))
                .select('VH')
                .mean()
                .set('month', m);
}));

and a working example can be found here


texts: Creating polygon from leased unit metes and bounds calls using ArcMap?

texts:
I have a dataset containing with metes and bounds listing an area.
I want to build a layer from these calls and produce the lines (polygons) that show the area.

I know there may be a bit of inaccuracy due to ground based measurements, etc., but I have to start somewhere so I am looking for a process to create the layer and begin to build lines based on bearing and range, mostly.

I tried in ArcMap 10.6 to create a line layer in a Geodatabase, but I am guessing I need to create certain fields to make this work and I do not know what fields to start with.

Overall, what I am ASKING for is a simple tutorial that will start with nothing but an empty map, and allow a layer to be created and populated with lines, (or polygons,) based on inputs like this:

S 09-32 W 387.9 feet
S 61-40 E 275.1

One call appears to describe a curve:

  "S 20-14 E 9002.0 feet, with a curve to the left having a Radius of
  832.6 feet, an ARC length of 888.3 feet, a Chord length of 846.8 feet, and a Chord bearing of S 17-35 E to the Point of Beginning"

Be aware that there are many tutorials about this sort of thing, but I have not found one that will take me from nothing to the end, without making assumptions that other parcels or polygons are already created in a layer, or that a Geodatabase already properly exists.


texts: Adding text to label expression in QGIS?

texts:
Using QGIS 2.18

I have line feature and want to add a foot tick mark and a deg symbol to label
the label as of now  shows

"dis"   ||  "Azimuth"

339045

I would like it to show

3391' 45 deg

The expression is simple and i know the answer is simple as well.


texts: Can you insert a marker symbol at a specific scale?

texts:
I have a map that is consistently used at a 1:500 scale. I want to be able to add structures to the map to view and move them around to get the best fit. The end user is new to GIS and I am trying to limit the amount of editing they need to do and would like to use marker symbols. They can easily be added to the map, rotated and moved. I need the structures to come in at the proper scale. Currently this is being done using a paper copy of the map and scaled paper cut outs of the structures. Is it possible to do this using ArcMap 10 of 10.1?

Thanks


texts: Seeking tutorials on raster data formats and raster manipulation?

texts:
I am looking to improve my understanding of raster formats and how to manipulate them. I do not really want to go into the mathematics of it, but really want to know the differences between the different formats available, best practices, performance optimization etc.

Does anyone have any good sources which they would like to share?


texts: How to rotate a SVG marker via a data-defined override to the value of one field in a layer

texts:
I have a pyqgis stand alone app that display linear features. I put an arrow svg symbol (QgsSvgMarkerSymbolLayerV2) on the line, now I would like to rotate it based on the 'value' field of each line feature in the layer. I tried some variations to the suggested approach below, but they are not working:

sl = QgsSvgMarkerSymbolLayerV2('arrow.svg')
lDD = QgsDataDefined(True, True, 'CASE WHEN "value" < 0 THEN 180.0 ELSE 0.0 END CASE', 'value')
sl.setDataDefinedProperty('angle', lDD)


texts: Android GPS: average location resolution time

texts:
What is the average location resolution time on a SmartPhone(Android) GPS?

Assumptions :
 The user is under the open sky.


texts: Drawing roads with different widths from GeoJSON using C#?

texts:
I'm drawing roads with different widths. The road data is constructed using GeoJson. (Each road is represented with a LineString.)

When I draw the roads with almost the same width, a trivial, direct drawing solution, i.e., each road is drawn as a polyline, it will look good:

However, if some roads have huge width, things become complicated. Direct drawing will result into some offset at the junction:

I think I can get the contour path of the road with huge width, and compute the intersect point of the wide road and the roads that it crosses with.

Then an irregular polygon is drawn to make it the cross smoother. But this is not elegant.

I'd like to explore a more elegant way to fix this despite the performance.


texts:
I selected buildings from OSM (which represents my city). I insterted a new column 'height' and I filled with a random number exactly number 2, so every building have 2metres altitude. Buildings are in 4326 and these are stored in the geom field. I would like to extrude all building with the ST_Extrude function.

I tried : SELECT gid, ST_Extrude(geom, 0,0, height) as geom
FROM csong_buildings;

not working:


texts:
I converted a raster file (one band grey scale) to an numpy array and set a range to it. I Did like that:

import gdal
import numpy
ds = gdal.Open('H:/GDAL_docs/a.tif')
myarray = numpy.array(ds.GetRasterBand(1).ReadAsArray())
selection = numpy.logical_or(myarray >= 0.3557, myarray <= 0.0753)

That works fine. Now I have the problem to convert the numpy array varriable back to a raster file (my favorite format would be GeoTIFF) is that possible? I also struggle to install rasterio to my osgeo4w64 instalation.


texts: How do I calculate the Terrain Ruggedness Index for each county in the United States in R?

texts:
I have a shapefile of US counties and high-resolution elevation data that spans the entire contiguous United States. My goal is to calculate a terrain ruggedness index for each county. The functions (that I've been able to find, e.g. spatialEco::tri) all take raster layers as arguments.

Based on mdsummer's excellent answer and given a boundary shapefile and a raster layer of elevation data, it's easy to calculate zonal statistics:

require(sf)
require(tidyverse)

# Shapefile of US counties in California
calif <- USAboundaries::us_counties("1960-01-01", resolution = "high", states = c("CA")) %>%
  mutate(county_fips = as.numeric(fips)) %>%
  select(county_fips, geometry)

# Load elevation data (at a low resolution for now)
elev <- elevatr::get_elev_raster(as(calif, "Spatial"), z = 2, src = "aws")

# Group the elevation raster according to county_fips
polymap <- fasterize::fasterize(calif, elev, field = "county_fips")
elev[is.na(values(polymap))] <- NA

# Zonal statistics
# v <- raster::values
zonal_stats <- tibble(value = raster::values(elev),
                      county_fips = raster::values(polymap)) %>%
  group_by(county_fips) %>%
  summarize(mean_elev = mean(value))
map <- left_join(x = calif, y = zonal_stats, by = "county_fips")
plot(map["mean_elev"])

I'm having difficulty seeing how to apply a function that takes a raster layer to each county individually. If I run the following code:

# Terrain Ruggedness Index (entire state)
tri.calif <- spatialEco::tri(polymap)
plot(tri.calif)

tri.calif.crop <- crop(tri.calif, extent(calif))
plot(tri.calif.crop)
plot(st_geometry(calif), add = TRUE)

this calculates the TRI across the state using the default cell size of the tri function:

but obviously these calculations aren't happening strictly within each county. How do I apply a function (like tri) that takes a raster layer to the raster that's contained within each county individually?

Once I have that, it's easy enough to calculate the mean TRI across all cells within the county, for example, using the same zonal statistics approach described above?

Once I have that


texts: Converting Feature Class to Coverage in ArcPy?

texts:
Not much info out there these days on coverages in PYTHON.

I need to convert a feature class to a coverage but cannot find out how to do it in python.

All I have found is "http://resources.esri.com/help/9.3/ArcGISDesktop/com/Gp_ToolRef/conversion_tools/feature_class_to_coverage_conversion_.htm"
but this doesn't get me far.

I have all required licenses.

I have not scripted out the example yet as its formatting doesn't make sense to me.

With help from below, I was able create a coverage. That coverage still needs to have Labels created and the arcpy.CreateLabels_arc is not recognized by arcpy. Not much help out there regarding this either.

I do have ArcINFO installed. I am wondering if I have to create an arcinfo workspace as referenced in http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//0017000000ps000000 .

If this is the case how do I then implement this for my labels?

arcpy.FeatureclassToCoverage_conversion(env.workspace + '\\' + r + ' POLYGON', TEMP + str(FMUNAME) + '_Cov' , "", "DOUBLE")


texts: Overview map size related to the main map rule

texts:
My question is as simple as it's title implies.

I have a map scaled in 1:20.000 and print it in a A0 paper.

Also, I want to add an overview map in it.

I wanted to ask if there is any rule regarding the printed size of the frame of the overview map or should I put its size just from experience?