texts:
I have a large QGIS projekt with about 50 layer, several of which are QuickMapService Base maps or other WMS-Layers. I work with this project for over a year now and had no problems.
I run QGIS on a Manjaro-Linux platform. Recently my QGIS was upgraded to 3.12 against my will and was very unstable afterwards. I reinstalled the LTR 3.10.6 from the AUR. Since then the following problem occurs very frequently.
I can load my project without issues and work for a while, but when zooming QGIS appear to have issues loading the maps. The map window turn blank and the status bar indicates constant activity, for as long as you want to wait. Neither raster nor vector data are rendered. When trying to close QGIS it indicates, that a process is still running. When restarting QGIS, the project works fine until the issue starts again.


texts: OpenLayers 4 adding a georeferenced image

texts:
I'm trying to add a georeferenced image layer to my map object. I'm using Node, with the v4.2.0 Openlayers module.

My code to add the georefenced image is as follows. Within the layers array of the map object (there is another OSM WMS background layer):

new ol.layer.Image({
  source: new ol.source.ImageStatic({
    attributions: '© some note',
    url: __dirname + '/image.tiff',
    projection: 'EPSG:3857',
    imageExtend: [486303.9778100852272473,
                  800886.4222691040486097,
                  486822.4427326705190353,
                  6801434.6447209259495139]
  })
})

Using the ol-debug.js library I get the following error:

Uncaught TypeError: Cannot read property '2' of undefined at Object.ol.extent.intersects > at ol.source.ImageStatic.getImageInternal > at ol.source.ImageStatic.ol.source.Image.getImage > at ol.renderer.canvas.ImageLayer.prepareFrame > at ol.renderer.canvas.Map.renderFrame > at ol.Map.renderFrame_ > at ol.Map.<anonymous>.

What could be wrong here?
Meanwhile I have used a workaround by creating a tiled XYZ image which works, but I would still like to know what is going wrong here.


texts: TimeManager on PyQGIS: InvalidTimeLayerError: ValueError('min() arg is an empty sequence',)

texts:
I am running a map showing multi-temporal data using PyQGIS and the Time Manager plugin.

I believe that these lines in the code may be a problem.

settings = layer_settings.LayerSettings()
settings.layer = earthquakesLayer
settings.startTimeAttribute = "DATE_AND_TIME"
settings.endTimeAttribute = "DATE_AND_TIME"
timeLayer = timevectorlayer.TimeVectorLayer(settings,iface=iface)

The last line of the code is giving the following error.

Traceback (most recent call last):
  File "<input>", line 1, in <module>
  File "C:/Users/Aaron/.qgis2/python/plugins\timemanager\timevectorlayer.py", line 100, in __init__
InvalidTimeLayerError: ValueError('min() arg is an empty sequence',)

Is there anything wrong with the arguments I have put in the Time Manager?

Here is the whole code I am trying to run.

from PyQt4.QtGui import*
from timemanager import*
#set background of map to a certain colour
iface.mapCanvas().setCanvasColor(QColor(196,229,248))
iface.mapCanvas().refresh

countries = "C:\Users\Aaron\Documents\pygis\world_boundaries.json"
earthquakes = "C:\Users\Aaron\Documents\pygis\Earthquakes.geojson"
basemap = QgsVectorLayer(countries, "World", "ogr")

#colour the basemap to a light grey colour
gray = {'color':'220,221,222'}
mapSym = QgsFillSymbolV2.createSimple(gray)
renderer = QgsSingleSymbolRendererV2(mapSym)
basemap.setRendererV2(renderer)

earthquakesLayer = QgsVectorLayer(earthquakes, "Earthquakes", "ogr")
#QgsMapLayerRegistry.instance().addMapLayer(earthquakesLayer)
earthquakesLayer.setLayerTransparency(30)

magnitude = (("Light",4.5,4.7,"253,255,22","2"),("Moderate",4.7,5.0,"253,190,22","5"),("Strong",5.0,5.4,"253,122,22","7.5"),("Major",5.4,6.2,"253,99,22","9.75"),("Great",6.2,7.2,"253,0,22","12"))

ranges = []
for label, lower, upper, color, size in magnitude:
    props={}
    props["color"] = color
    props["size"] = size
    props["outline_width"] = '0'
    props["outline_color"] = "0,0,0,0"
    sym = QgsMarkerSymbolV2.createSimple(props)
    rng = QgsRendererRangeV2(lower, upper, sym, label)
    ranges.append(rng)

field = "mag"
renderer = QgsGraduatedSymbolRendererV2(field, ranges)
earthquakesLayer.setRendererV2(renderer)

QgsMapLayerRegistry.instance().addMapLayers([earthquakesLayer,basemap])

settings = layer_settings.LayerSettings()
settings.layer = earthquakesLayer
settings.startTimeAttribute = "DATE_AND_TIME"
settings.endTimeAttribute = "DATE_AND_TIME"
timeLayer = timevectorlayer.TimeVectorLayer(settings,iface=iface)
try:
    ctrl = Qgis.utils.plugins['timemanager'].getController()
except:
    ctrl = qgis.utils.plugins['TimeManager'].getController()
tlm = ctrl.getTimeLayerManager()
tlm.registerTimeLayer(timeLayer)
tlm.setTimeFrameType('months')
ctrl.setLoopAnimation(True)
start = time_util.str_to_datetime(timeLayer.getMinMaxValues()[0])
end = time_util.str_to_datetime(timeLayer.getMinMaxValues()[1])
tlm.setCurrentTimePosition(start)
ctrl.setAnimationFrameLength(1000)
ctrl.toggleAnimation()
ctrl.toggleAnimation()


texts: Shapefile vector point to ppp spatstat with polygon windows

texts:
I always get error to this when converting shapefile point to ppp spatstat that windowed by a polygon shapefile
# polygon that to be window
neighborhoods <- st_read("neighborhoods/neighborhoods.shp")

# convert CRS to planar projection as recommended by (https://stackoverflow.com/questions/59597078/use-sf-polygon-object-as-window-in-spatstat)
neighborhoods_3857 <- st_transform(neighborhoods, crs = 3857)

# point that to be PPP spatstat
trees <- st_read("trees/trees.shp")

# convert to planar projection
trees_3857 <- st_transform(trees, crs = 3857)

The problems, the "trees_3857" doesn't have dataframe columns that represent in EPSG3857 coordinates, so Feature column of "trees_3857" doesn't have x and y columns that respect to EPSG 3857
q <- ppp(x=?, y=?, win=neighborhoods_3857)

what I have done but error
first =using as.ppp instead of ppp to convert a vector shp to ppp spatstat
test1 <- as.ppp(trees_3857, win=neighborhoods_3857)
Error in as.ppp.sf(trees_3857, win = neighborhoods_3857): unused argument (win = neighborhoods_3857)

second = add geometries of "trees_3857" to feature columns with QGIS named "xcoord" & "ycoord"
test2 <- ppp(trees_3857$xcoord, trees_3857$ycoord, win = neighborhoods_3857)
Warning message:
"65217 points were rejected as lying outside the specified window"


texts: Print composer layout without opening project in GUI (no iface solution)

texts:
I would like to print some composers without open the files that contains them in QGIS.

without using templates
not using iface
using,

app = QgsApplication(sys.argv, gui_flag)
app.initQgis()
project.read(QFileInfo(project_path))

for QGIS 2.18.x (then I can't use QgsLayoutManager)

How can I call a composer without iface and layout manager?


texts: Getting Elevation Data for Elevation Attribute using ArcGIS Pro

texts:
In ArcGIS Pro, how can I insert the elevation data into an attribute table?

There is already an elevation field created, and the feature class is set to where all features are on the ground, so you would think that when features are created that they would automatically take on the ground elevation. I feel like this is the case, anyway, but I just can't figure out how to get this data into the attribute table.

There isn't a Z coordinate system set on the map, but the XY is in North Carolina StatePlane.

I have elevation data for other points in the feature class, but that's because I got them straight from surveyors.


texts: ALOS DEM shows continious no data how to over come?

texts:
How to overcome No Data in ALOS DEM which acquired at location of Ernakulam, Kerala. Attached .JPG screen shot of my screen.

Why it is showing NODATA??


texts: Recording scale in attributes when drawing a polygon in QGIS

texts:
In QGis I would like to know how to get the scale of the vue when drawing a new polygon?

The idea is to be able to estimate the quality of the polygons drawn by my colleagues and me and to know which ones can be used for official purpose.


texts: How to batch classify a multi-temporal raster image stack?

texts:
I am using dense time stacks of remotely sensed imagery to analyse vegetation change in the Arctic from 1984-2010 for one scene. I will be needing to classify each image, likely using the Maximum Likelihood Classifier tool in ArcMap 10.2. In the Batch mode of the tool, it just seems to allow importing multiple raster bands, then exporting to a raster output. Is there a way to feed multiple raster .tiffs in to the same classification schema? I read a paper (Castrence et al 2014) where they performed "a single classification over the entire time span, thereby [avoiding] problems associated with multiple independent classifications of images from different dates." Could you build one signature file and then use it for all the rasters? I am not experienced with Pyton scripting, but I'm getting the feeling this will be needed...I know that the images must be relatively calibrated against a master. Basically I'm lazy and I want to avoid doing things manually.

I also have PCI Geomatica, if anyone knows of a batch classification tool/model with EASI.


texts: Dynamic Content and Infowindow-Size & -Position

texts:
I have to display Dynamic Content (aka embeds … tweets, youtube-videos, facebook-posts etc.) inside a Infowindow … now how can I make sure to resize and position the infowindow properly so the following will not happen:

Pic 1: This is perfect … (but has no embeds…)

Pic 2: … and this not so much. Because of the embeded Content the Infowindow gets resized to the bottom … the «arrow» or the Infowindow-«Bottom» should be @ the height where you can spot «Rapp» on the map

How can I achieve that the infowindow gets resized (more height to top, not to bottom), so it stays positioned correctly.


texts: Building georeferenced GeoTIFF from mbtiles automatically?

texts:
I want a georeferenced raster.In that way I don't have to georeference this raster manually .
I know that geotiff has the georeference inside it , I got the mbtiles of the map how do I build the georeferenced geotiff.

my mbtile mbtile file tronsform it to Geotiff


texts:
I have downloaded a US federal lands data layer (fedlanp010g.shp.tar.gz) from http://nationalmap.gov/small_scale/atlasftp.html?openChapters=chpbound#chpbound

In R, I am trying to create a buffer around each polygon, However gBuffer crashes when the code below is run. I added a pause (sys.sleep(1)) to add a delay, and I received the error message "Value of SET_STRING_ELT() must be "CHARSXP" not a 'NULL'" at row 14.

Here is the code I have been using (I am processing with a loop):

library(rgeos)
library(raster)

fed = readOGR("downloaded file from above here")
i = 1
temp = fed[i,]
buff = gBuffer(temp,byid=TRUE,width=25)
result = erase(buff,temp)

for(i in 2:nrow(fed)){
    temp = fed[i,]
    buff = gBuffer(temp,byid=TRUE,width=25)
    temp_buff = erase(buff,temp)
    result = rbind(result,temp_buff)
  }


texts:
I have a feature class and table. The table contains about 18000 data. the polygon feature class contain 269 data. Both common column are [MAHALLE_AD] and [Mahalle]. Also Both are  same personal geodatabase. I want to transfer from table data (18000) to feature class using common column.

I get comment in verifying expression ( The expression was verified successful, but no records were returned). I couldn't find a solution.


texts: GRASS GIS 7.4 in Linux for command line usage - setting up environment (.bashrc file)

texts:
I need to use GGRAS GIS 7.4.0 on Linux server (Debian 8.10 Jessie) via command line (with flag --exec) and I'm still having troubles with it.

Description: I compile and install GRASS GIS 7.4.0 with no problem on that machine. When I type

grass74

in terminal, the text GUI launch. But when I type some command ie.

grass74 v.import input=some_shapefile.shp output=some_file

or

grass74 /root/grass74/test_loc/PERMANENT v.import input=some_shapefile.shp output=some_file

it gives an

root@maptest:~# grass74 /root/grass74/test_loc/PERMANENT v.import input=/root/test_shp/test_point.shp output=bod
Cleaning up temporary files...
Starting GRASS GIS...
ERROR: </root> is not a valid GRASS Location because PERMANENT Mapset is missing
Exiting...

So I tried to set up .bashrc file for setting up GRASS GIS environments like at this page Working with GRASS without starting it explicitly.

In bashrc file, I added this lines:

export GISBASE=$GISBASE:/usr/local/grass-7.4.0
export PATH="$GISBASE/bin:$GISBASE/scripts:$PATH"
export LD_LIBRARY_PATH="$GISBASE/lib"
export GRASS_LD_LIBRARY_PATH="$LD_LIBRARY_PATH"
export PYTHONPATH="$GISBASE/etc/python:$PYTHONPATH"
export MANPATH=$MANPATH:$GISBASE/man

generate GISRCRC
MYGISDBASE=$HOME/grass74
MYLOC=test_loc
MYMAPSET=PERMANENT

# Set the global grassrc file to individual file name
MYGISRC="$HOME/.grassrc74"

echo "GISDBASE: $MYGISDBASE" > "$MYGISRC"
echo "LOCATION_NAME: $MYLOC" >> "$MYGISRC"
echo "MAPSET: $MYMAPSET" >> "$MYGISRC"
echo "GRASS_GUI: text" >> "$MYGISRC"

Then restart server and type "grass74" again, but got this error:

root@maptest:~# grass74
Cleaning up temporary files...
Traceback (most recent call last):
  File "/usr/local/bin/grass74", line 2003, in <module>
    main()
  File "/usr/local/bin/grass74", line 1884, in main
    clean_temp()
  File "/usr/local/bin/grass74", line 1621, in clean_temp
    call([gpath("etc", "clean_temp")], stdout=nul)
  File "/usr/local/bin/grass74", line 245, in call
    return subprocess.call(cmd, **kwargs)
  File "/usr/lib/python2.7/subprocess.py", line 522, in call
    return Popen(*popenargs, **kwargs).wait()
  File "/usr/lib/python2.7/subprocess.py", line 710, in __init__
    errread, errwrite)
  File "/usr/lib/python2.7/subprocess.py", line 1335, in _execute_child
    raise child_exception
OSError: [Errno 2] No such file or directory

Does anybody knows what I'm doing wrong? Bad setting up environmental variables (bashrc file)? Or something else?

EDIT:
I tried to modify command according to mankoff answer. I have created location at:

/root/grass74/pokus1/PERMANENT

and the command is:

grass74 EPSG: /root/grass74/pokus1/PERMANENT --exec v.import input=/root/test_shp/test_point.shp output=test_point

but I still get:

Cleaning up temporary files...
Starting GRASS GIS...
Executing <v.import input=root/test_shp/test_point.shp output=test_point> ...
Traceback (most recent call last):
  File "/usr/local/grass-7.4.0/scripts/v.import", line 291, in <module>
    sys.exit(main())
  File "/usr/local/grass-7.4.0/scripts/v.import", line 173, in main
    **vopts) == 0:
  File "/usr/local/grass-7.4.0/etc/python/grass/script/core.py", line 413, in run_command
    ps = start_command(*args, **kwargs)
  File "/usr/local/grass-7.4.0/etc/python/grass/script/core.py", line 380, in start_command
    return Popen(args, **popts)
  File "/usr/local/grass-7.4.0/etc/python/grass/script/core.py", line 74, in __init__
    subprocess.Popen.__init__(self, args, **kwargs)
  File "/usr/lib/python2.7/subprocess.py", line 710, in __init__
    errread, errwrite)
  File "/usr/lib/python2.7/subprocess.py", line 1335, in _execute_child
    raise child_exception
OSError: [Errno 2] No such file or directory
Execution of <v.import input=root/test_shp/test_point.shp output=test_point> finished.
Cleaning up temporary files...

EDIT - OGRinfo output

root@maptest:~# ogrinfo /root/test_shp/test_point.shp
INFO: Open of `/root/test_shp/test_point.shp'
      using driver `ESRI Shapefile' successful.
1: test_point (Point)


texts: Setting label at top right of MultiLineString in OpenLayers 4

texts:
I  would like to set the label a the top right of the geometry because i want to see it entirely.

I dont want the label overlapping the geometry.

I tried with the extent without success.

var labelStyle = new ol.style.Style({
  geometry: function(feature) {
    var geometry = feature.getGeometry();
    if (geometry.getType() == 'MultiLineString') {
      // Only render label for the widest polygon of a multipolygon
      var polygons = geometry.getLineStrings();
      var widest = 0;
      for (var i = 0, ii = polygons.length; i < ii; ++i) {
        var polygon = polygons[i];
        var width = ol.extent.getWidth(polygon.getExtent());
        if (width > widest) {
          widest = width;
          geometry = polygon;
        }
      }
    }
    return geometry;
  },
  text : new ol.style.Text({
    font: '12px Calibri,sans-serif',
    fill: new ol.style.Fill({color: 'black'}),
    //stroke: new ol.style.Stroke({color: 'white', width: 1}),
    overflow: true,
    backgroundFill : new ol.style.Fill({color: 'orange'}),
    placement : "point",
    textBaseline : "top"
  })
});

  var styles = [geomStyle, labelStyle];

  var vector = new ol.layer.Vector({
    source: vectorSource,
    style: function(feature) {
      labelStyle.getText().setText(feature.get('QUALIF') +' '+feature.get('DENOM'));

var styles = [geomStyle, labelStyle];

var vector = new ol.layer.Vector({
  source: vectorSource,
  style: function(feature) {
    labelStyle.getText().setText(feature.get('QUALIF') +' '+feature.get('DENOM'));

    var geometry = feature.getGeometry();
    var ext = geometry.getExtent();
    var offsetX = (ext[2]-ext[0]);
    var offsetY = (ext[3]-ext[1]);
    console.log(offsetX);
    /*
    labelStyle.getText().setOffsetX(offsetX);
    labelStyle.getText().setOffsetX(offsetY);
    */
    return styles;
  },
  declutter: true
});

The full code on GITBUB


texts: Copying fields from one attribute table to another in QGIS

texts: I have two layers. One is point, the other one a line. They both share same x and y coordinate. How do I transfer one field from the point layer to another field in the line layer based on the same coordinates.


texts:
I have a point feature class and a polygon feature class, and each point is located within one polygon. Is there a way to programmatically calculate the distance from each point to the farthest location on its containing polygon's edge (see example below)?

I have an ArcGIS Desktop basic licence.


texts: gdal_translate SQLite driver has no raster capabilities

texts:
I am trying to convert a raster to SQLite using rasterlite2 and gdal a la:
gdal_translate -f SQLite raster.img raster.sqlite
but after running the command I get the following error:
ERROR 1: SQLite driver has no raster capabilities.

I am using Ubuntu 20.04 and have installed the following with apt install

gdal-bin
libgdal-dev
sqlite3
spatialite-bin
libspatialite-dev
rasterlite2-bin
librasterlite2-dev

I thought these are the only dependencies, but maybe I am missing something?
Does anyone spot anything missing?


texts: Clipping raster with vector boundaries using QGIS?

texts:
I am using QGIS.  I would like to clip a raster precipitation layer using an admin boundary layer that is vector data.  However the geoprocessing tools seem to be usable only for vector data.

How can I clip this precipitation layer?


texts: Clipping a polygon from a raster without georeference

texts:
Is it possible to clip a polygon from a raster layer without being georeferenced?
When I try I get this error message
 "ERROR 1: Unable to compute a transformation between pixel/line
and georeferenced coordinates"
and am not really sure what it means. I'm speculating that it wants me to georeference the raster layer but I can't as it is a photo, not a map (not a photo of a map) and if I just georeference arbitrarily to try and placate it then the pictures warps.


texts:
I have a layer with several shapes and from here I want to delete the lines. I tryed Toggle editing and Select features by polygon but when closing the polygon, it dissapears. Also I tried to delete all the nodes that I don't want by selecting them with the Node Tool buton, but nothing happened. What do you suggest me to do?

This is my example:


texts:
I have the following code in Python for a tool I want to build. For some reason I viewshed don't exist in the folder I have chose (out). Also when I am trying to define the output workspace every time the file I have created inside the Input folder is deleted so I created one more outside of the folder where my input data are. Why is that? I am using ArcGis 10.2.2

import arcpy, os
from arcpy import env
from arcpy.sa import *

env.workspace = arcpy.GetParameterAsText(0)
out = arcpy.GetParameterAsText(1)
fc = arcpy.ListFeatureClasses("Pnt*", "Point")
ras = arcpy.ListRasters("clip*", "GRID")
point = "Pntclip_pol1"
dem = "clip_pol1"
i = 1
for shp in fc:
    for raster in ras:
        if (shp == 'point' and raster == 'dem'):
            inRaster = raster
            inObserverFeatures = shp
            outViewshed = Viewshed(inRaster, inObserverFeatures, "")
            outViewshed.save(out + "view" + str(i))
            i = int(i) + 1
            point = "Pntclip_pol" + str(i)
            dem = "clippol" + str(i)


texts:
Is there a standardized web feature (WFS) or web processing query to get a single pixel value from a map?

I don't need a whole tile, which seems to be what the WFS and WMS standards are designed for. I just want the pixel value at a specific latitude / longitude.

Is there a standard for this type of query?


texts: QGIS extract part of description to another attribute table column

texts:
I have an attribute table with the fields "name" and "desc".
In the "desc" field i have the following string: use:xxx population:xxx.
The xxx are examples given that it is a very large dataset.

I want to extract the xxx after the population to another field called "population" using field calculator.

In the image bellow you can see the real format of the field "desc". My intention is to extract the number in float after the population (just the number).


texts: Combining 2 different symbologies in QGIS

texts: I've got a layer (road axis) with 2 different attributes : one with the amount of trafic, another with the number of people living next to this road. As symbology for the trafic amount, I like use different thicknesses. For the number of people, I want to use a different colors. How can I make a combination of these 2 symbologies.


texts:
I am having an issue with OGR's ExecuteSQL, but not with the actual SQL query.

What I am intending to do is to take a large number of polygons ~500k and find out if they intersect with a multipolygon.

The SQL query takes approximately 4 minutes (Which some may say is slow, please point out if I'm doing it wrong).

The critical issue here is, when I try to work with the results of the query, it is incredibly slow. It takes approximately 2hrs30min for it to loop through all of the results and add it to a new file.

This issue persists even when I replace the contents of the loop with a simple "pass" statement. To further the issue, trying to "print len(Intersection_Results)" takes so long I stopped the code after 5 minutes.

Please see my code below:

Multipolygon_file= ogr.Open(r"multipolygon.sqlite",1)#This file contains both the polygons and the multipolygon in two different layers

SQL_COMMAND = "SELECT polygons.* FROM polygons, multipolygon WHERE
ST_Intersects(polygons.geometry, multipolygon .geometry)"
Intersection_Results = Multipolygon_file.ExecuteSQL(SQL_COMMAND, dialect = "sqlite")

Results_File= ogr.Open( r"Results_Table.sqlite",1)
Results_File_layer = Results_File.CreateLayer("Intersect1000yr_Results",projection,geom_type=ogr.wkbPolygon)

laydefn = Results_File_layer.GetLayerDefn()
Results_File_layer.StartTransaction()
for feature in Intersection_Results :
    outFeature = ogr.Feature(laydefn)
    outFeature.SetGeometry(feature.GetGeometryRef())
    Results_File_layer.CreateFeature(outFeature)

Results_File_layer.CommitTransaction()
Multipolygon_file.ReleaseResultSet(Intersection_Results)
Results_File.FlushCache()

Just to further add to the list of other things I have tried to fix this issue:
I have limited the results in the SQL query to only give me the first 1000 results so that I can check the outputs to make sure they are what I expect: a layer containing only the polygons that intersect with the multipolygon. It was exactly what I expected

I also run this EXACT query before running it with the multipolygon, I use a bounding box with the same extent of the multipolygon to cut down the data from  approx 5.8m polygons to 800k. This takes approimately 1min to execute the SQL query and a further ~1min to write the results to the multipolygon.sqlite file.


texts: Intersection of raster and vector - road file and maps (tif)

texts:
ArcGIS Desktop 10.6.1

I have a railroad vector layer (with 10 rows, that is 10 segments) and 30 USGS historic maps (tif files) that cover that route. I'd like to find the intersection, such that one of the attributes in the vector layer is where it crosses these tif files. That is, which maps do I need to use to digitize each section of the road.

I'd like the attribute table in the vector file to list the name of the tif file. I have attached an image.


texts: Error while removing gaps with v.clean in QGIS

texts:
I have a polygon layer (Geopackage with 378 features in EPSG:3857) with many very small gaps between them (topology checker lists more than 5000). The gaps are in a range of several hundred meters to fractions of millimeters in map units (meters). Using QGIS 3.16, I tried Grass v.clean. However, it is not able to repair the geometries. No output is generated.
Error messages in the protocol mention:

ERROR 1: osgeo::proj::io::WKTParser::Private::buildEllipsoid: non
double value
WARNING: Unable to determine input map's vector feature type(s).
ERROR: Unable to create OGR spatial reference

I tried also QGIS native repair geometries as well as fill gaps algorithm - the resulting output has the same gaps.
What else could I try to remove these gaps?


texts: How can we improve the OpenLayers documentation?

texts: I have difficulty understanding the documentation. Is it just me?  Are there alternative sources of information?


texts: Spatially join point and polygon attributes in FME

texts:
I'm new to FME, and most of the time I work with ArcGIS Pro. I have a few building polygons, and a few points containing additional information. Now, in ArcGIS Pro I can spatially join these two, creating one layer with the polygon "look" but added with the point information. This then needs to be exported to a new shapefile.

How can I recreate this in FME? I now used the PointOnAreaOverlayer, but when I export the area data, it still only contains the attributes from the input area, and the attributes from the points are not added. When I look at the area data with the inspector, it does show an attribute table wit the combined data, but when writing this to a new file, it does not.

Can someone help me out?


texts:
I was trying to install QGIS on centos 6.3 and follow the instruction here.

However after running yum install qgis qgis-python qgis-grass qgis-mapserver I got below errors:

    ---> Package libvpx.x86_64 0:0.9.0-8.el6_0 will be installed
    ---> Package python-setuptools.noarch 0:0.6.10-3.el6 will be installed
    --> Finished Dependency Resolution
    Error: Package: grass-6.4.2-8.el6.x86_64 (elgis)
               Requires: libgeos-3.3.8.so()(64bit)
               Available: geos-3.3.8-2.el6.x86_64 (elgis)
                   libgeos-3.3.8.so()(64bit)
               Installed: geos-3.4.1-1.el6.x86_64 (@opengeo)
                   Not found
               Available: geos-3.2.2-1.el6.rf.x86_64 (rpmforge)
                   Not found
               Available: geos-3.2.2-2.el6.x86_64 (elgis)
                   Not found
               Available: geos-3.3.0-1.el6.x86_64 (elgis)
                   Not found
               Available: geos-3.3.1-2.el6.x86_64 (elgis)
                   Not found
               Available: geos-3.3.2-1.el6.i686 (epel)
                   Not found
               Available: geos-3.3.6-1_0.el6.x86_64 (elgis)
                   Not found
    Error: Package: grass-libs-6.4.2-8.el6.x86_64 (elgis)
               Requires: libgeos-3.3.8.so()(64bit)
               Available: geos-3.3.8-2.el6.x86_64 (elgis)
                   libgeos-3.3.8.so()(64bit)
               Installed: geos-3.4.1-1.el6.x86_64 (@opengeo)
                   Not found
               Available: geos-3.2.2-1.el6.rf.x86_64 (rpmforge)
                   Not found
               Available: geos-3.2.2-2.el6.x86_64 (elgis)
                   Not found
               Available: geos-3.3.0-1.el6.x86_64 (elgis)
                   Not found
               Available: geos-3.3.1-2.el6.x86_64 (elgis)
                   Not found
               Available: geos-3.3.2-1.el6.i686 (epel)
                   Not found
               Available: geos-3.3.6-1_0.el6.x86_64 (elgis)
                   Not found
     You could try using --skip-broken to work around the problem
     You could try running: rpm -Va --nofiles --nodigest

my yum repolist shows:

    yum repolist
    Loaded plugins: fastestmirror, refresh-packagekit, security
    Loading mirror speeds from cached hostfile
     * epel: ftp.jaist.ac.jp
     * rpmforge: ftp.heanet.ie
    repo id                                                      repo name                                                                                status
    base                                                         CentOS-6 - Base                                                                                6,381
    elgis                                                        EL GIS 6 - x86_64                                                                                198
    epel                                                         Extra Packages for Enterprise Linux 6 - x86_64                                                     9,757
    extras                                                       CentOS-6 - Extras                                                                                13
    opengeo                                                      opengeo                                                                                52
    rpmforge                                                     RHEL 6 - RPMforge.net - dag                                                                        4,638
    updates                                                      CentOS-6 - Updates                                                                                1,326
    repolist: 22,365

How can I proceed to install the QGIS?

I tried the --skip-broken method but QGIS cannot be loaded.


texts:
I've been looking around but haven't found an answer to match.

I have a number of endpoint rates (EPR, I'm working on shorelines), and I want to know if it's possible to create a buffer based on the EPR numbers. So say a point has moved 2 meters, but another point has moved -0.5, is there a way to calculate a buffer based on each unique EPR number?
Hope it makes sense, and that I have clarified my question.

Edit (17.04.2019):
I'm using ArcMap, working with the DSAS software-extension. My plan is to try and create a disaster-assessment map, and I want to show how the shoreline could potentially move, based on past movement.
So I have a set of EPR numbers, representing whether the coast has eroded or accreted, and I want to, based on the EPR numbers, visualise where those points could be located within a number of years if that makes sense.


texts: What might cause "ERROR 000728: Field  does not exist within table" in this code?

texts:
I have seen ERROR 000728 several times and I could never figure out why it happens.

Here is a snippet from my code:

str_Dissolve = arcpy.Dissolve_management(Intersect_str, "in_memory/str_Dissolve", "FID_Con_st", "", "MULTI_PART", "DISSOLVE_LINES")

I do have a field named "FID_Con_st".

There is no space in the field name and it doesn't seem to be any refresh bug, because it persists after closing the script.  These were possible causes suggested by an earlier Question and an ArcGIS Discussion Form Posting.

Are there other possible causes that anyone can suggest?


texts: How to use scale-rule based renderer in QGIS?

texts:
I wish to render a shapefile of buildings such that at scales below 1:5000, they have a border, and above 1:5000, they do not have a border.

For this, I went into layer properties, selected Rule-Based, and then selected the current Rule and then went to Refine Current Rules>>Add Scales to rule. After this I set the symbol for scales >5000, like I wanted.

Here is my layer properties.

You can see that there is an initial symbol. How Do I remove this? This causes problem because it is present behind the symbol; for example You can see the boundary here even though the scale is greater than 1:5000:

I have found that that if I remove the boundary of the first symbol(the one in blue), the layer appears as expected at scales >1:5000. This seems like a workaround rather than the straightforward procedure. What is the correct procedure to use scale -rule based renderer in QGIS?


texts:
I have orthophoto(GeoTiff) and 4 GPS coordinates that represent rectangle.

How can I crop orthophoto by using GDAL? (I need to remove areas outside of rectangle)

UPDATE:

I tried gdal_translate

gdal_translate -projwin top_left_corner.longitute top_left_corner.latitude bottom_right_corner.longitute bottom_right_corner.latitude orthophoto.tif crop.tif

But seems it expects top left corner and bottom right corner in other coordinate system?


texts: GeoServer CSS label along polygon boundary - overflow

texts:
I labeled a polygon along it's edge by using the boundary of an inside buffer:

label-geometry: [boundary(buffer(GEOM, -0.01))];

This works, but my problem is, that at certain zoomlevels/positions the label runs out of the polygon instead of being labeled around the corner or on the (long enough) straight line.
Moving the map extend only a few pixels changes this again.

My CSS looks like this (inner black stroke is to depict the labeled boundary):

*{
  fill: #E30613 ;
  fill-opacity: 0.2;

  font-family: "SansSerif.plain";
  font-fill: #194383;
  font-size: 20;

  label-geometry: [boundary(buffer(GEOM, -0.01))];
  stroke-geometry: [GEOM],[boundary(buffer(GEOM, -0.01))];

  label: 'Really very very very long label';
  label-follow-line: true;
  label-max-angle-delta: 180;
  label-remove-overlaps: true;
  label-conflict-resolution: true;
  label-allow-overruns: false;

  stroke: #194383, black;
  stroke-width: 2, 0.1;
}

Does anyone know a way to force the label to alway be on the boundary?


texts: Calculate QUARTER from Date field using ArcGIS Field Calculator?

texts:
I currently have a Date field in the format of MM/DD/YYYY.

My goal here is to convert the Dates into 1 of 4 quarters, along with the Year.  The field to be created will be called QTR_YR.

My first question is what data type do I need to make the QTR_YR field?  Should it be Integer (Short or Long), or Date, or Double?

I am old school so I'd prefer VBA but I can use Python.


texts:
I'm importing datasets to CartoDB manually but the import queue seems to be stuck in 'Pending...'. I cannot delete these never ending imports as I can't see any button and when I try to upload more datasets it says:

Unfortunately you can only import up to 3 files at the same time

How can I resolve this?


texts: Displaying a PNG file over a OpenLayers v6 map

texts:
I am having difficulty finding a clear solution to this, I may just be over complicating it but I am not sure how to just add a static PNG file over top of a OpenLayers map with a opacity of 0.5. I've been searching the documentation and other related question but none have seemed to clear things up.

Could someone point me in the right direction, the data bounds are N: 85.044, S: -85.044, E: 180, W: -180. The image is of SMAP (Soil Moisture Active Passive) data I have processed with Python and exported as a PNG to display over this map in OpenLayers. Would leaving it as girdded data (simply stored in a matrix) make it easier/simpler to project?

This is  my image to add to the mapand it's size is 3856x1624. I know the black of the image will make it hard to see teh rest of the map, however right now I am just trying to get an image to overlay on the map.

Below is my most recent attempt to get it too work:

// Create the map object
const map = new Map({
target: 'map',
layers: [
new TileLayer({
  source: new OSM
  })
],
view: new View({
  projection: 'EPSG:4326',
  center: [0, 0],
  zoom: 0
  })
});

// Extent of the SMAP data bounds
imageextent = [-180, -85.044, 180, 85.044];

// Creating a new image layer using my PNG image as a static source
var image_layer = new ImageLayer({
  source: new Static({
    url: 'dataImage.png',
    imageExtent: imageextent,
    projection: 'EPSG:4326',
    crossOrigin: ''
  }),
  className: 'Image_Layer',
  opacity: 0.5,
  extent: imageextent,
  map: map  // This should add the image as an overlay
})

So when I run this using Node.js I have a map in EPSG:4326, but no image on top.


texts: ArcGIS JSAPI Query Related Records/Relationship Table

texts:
I am using this example https://developers.arcgis.com/javascript/3/jssamples/fl_query_related.html

to develop an application which queries related features.  I have two feature classes and have created a N-M relationship class.  The objective is to click on an electric feeder and query the related customers.  There are 26 electrical feeders  and over 400 customers.  My relationship class was based on the OBJECTIDs being the primary keys and the foreign keys being the "feederID" that the structures have and the that the customers are assigned to.

Ideally, I'd expect to see several customers per feeder, but based on my logic in this code sample, I am only returning the matching ObjectIDs per feeder.

What should be changed to return each customer which shares the same feeder?

 function findRelatedRecords(evt) {
        var features = evt.features;
        var relatedTopsQuery = new RelationshipQuery();
        relatedTopsQuery.outFields = ["*"];
        relatedTopsQuery.relationshipId = 0;
        relatedTopsQuery.objectIds = [features[0].attributes.OBJECTID];
        wellFeatureLayer.queryRelatedFeatures(relatedTopsQuery, function(relatedRecords) {
          console.log("related recs: ", relatedRecords);
          if ( ! relatedRecords.hasOwnProperty(features[0].attributes.OBJECTID) ) {
            console.log("No related records for ObjectID: ", features[0].attributes.OBJECTID);
            return;
          }
          var fset = relatedRecords[features[0].attributes.OBJECTID ];
          var items = array.map(fset.features, function(feature) {
            return feature.attributes;
          });
          //Create data object to be used in store
          var data = {
            identifier: "OBJECTID",  //This field needs to have unique values
            label: "OBJECTID", //Name field for display. Not pertinent to a grid but may be used elsewhere.
            items: items
          };

          //Create data store and bind to grid.
          store = new ItemFileReadStore({ data:data });
          grid.setStore(store);
          grid.setQuery({ OBJECTID: "*" });
        });
    }


texts: Conducting common-attribute dissolution ONLY on adjacent features using QGIS

texts:
I am working in QGIS. I have a shapefile of city parcels, which I am trying to dissolve where parcels have matching owners AND are touching. "Owner" is an attribute field.
How can I accomplish this?


texts: Interchange qgis._core.QgsVectorLayer  osgeo.ogr.Layer

texts: Is it possible to convert qgis._core.QgsVectorLayer instance to osgeo.ogr.Layer instance and vice versa in python script ?


texts: QGIS Atlas generation for multiple map extents?

texts:
I'm trying to create 5 maps in composer showing the same layer features but different extents of a map (attached an image to give everyone an example). I've opted to use atlas to save me manually creating 5 maps... but I have hit a barrier.

I think where I am going wrong is that I am trying to create an atlas based off a point layer (which is my coverage layer) so it is creating individual maps (about 60 +) for the points displayed.

However... I would like atlas to generate 5 maps for the different geographical extents that cover those points?

I right in thinking that I'll have to go and create a polygon shapefile that creates say a bounding box around 5 different areas to contain the points of interest? & then select this as my coverage layer when generating an atlas?

Does anyone know a way to achieve what I am trying to do? Any views or general help would be great.


texts: Adding categories classification to rule

texts:
How do I add sub categories to a rule with  automatic classification as I modify a field?
Basically I want to give the chance to display each single new feature.


texts: Mapping data: error in classIntervals: var is not numeric

texts:
When executing the 'classIntervals' command I got the error while mapping data from Eurostat.

library(rgdal)
library(RColorBrewer)
library(sp)
library(GISTools)
library(classInt)
library(maptools)
library(SmarterPoland)
library(reshape)
library(reshape2)

# create a new empty object called 'temp' in which to store a zip file
# containing boundary data
temp <- tempfile(fileext = ".zip")
# now download the zip file from its location on the Eurostat website and
# put it into the temp object
  download.file("http://epp.eurostat.ec.europa.eu/cache/GISCO/geodatafiles/NUTS_2010_60M_SH.zip", temp)
# now unzip the boundary data
unzip(temp))
plot(EU_NUTS)

proj4string(EU_NUTS)

EU_NUTS <- spTransform(EU_NUTS, CRS("+proj=merc +a=6378137 +b=6378137  +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs"))

plot(EU_NUTS)

EurostatTOC <- getEurostatTOC()

head(EurostatTOC)

data <- getEurostatRCV(kod = "lfst_r_lfu3rt")

unique(data$age)
unique(data$sex)
unique(data$time)

sub_data <- subset(data, (age == "Y20-64") & (sex == "T") & (time == "2012 "))

mapdata <- cast(data, geo ~ time + age + sex)

head(mapdata)

head(EU_NUTS@data)

EU_NUTS@data = data.frame(EU_NUTS@data, mapdata[match(EU_NUTS@data[, "NUTS_ID"],
                                              mapdata[, "geo"]), ])

my_colours <- brewer.pal(5, "BuGn")

breaks <- classIntervals(EU_NUTS@data$X2012._Y20.64_T, n = 5, style =   "fisher", unique = TRUE)$brks

## Error in classIntervals(EU_NUTS@data$X2012._Y20.64_T, n = 5, style = "fisher",  :  var is not numeric


texts: Mosaic Raster Function Error ArcGIS Pro 2

texts:
I have a class of GIS students who are attempting to use the Mosaic Raster Function in ArcGIS Pro to stitch together three small raster clips.  For most, the function is working fine but some are experiencing a general error where the function just does not work.  The error just states that it could not be completed.  They all seem to be using the same imagery as raster inputs.  They are all selecting the same operation and they are working on the same type of computer with the same version of ArcGIS Pro.

Is there some other variable that we should be considering which could be causing this very general error?


texts: Removing las2dem streaks from raster DEM?

texts: I've created a series of DEM files using the las2dem tool from LAStools; however since I don't have a licensed version they have black streaks running across them. Does anyone know a good method for removing these streaks?


texts: Reducing coordinate decimal places in WKT using C#?

texts:
I have imported several shapefiles into a Microsoft SQL Server 2017 database
When I use the STAsText() function on geography objects I get the following WKT:

SELECT [Adm0GeoBoundaryID]
      ,[Adm0Code]
      ,[Adm0Name]
      ,[Coords].STAsText() as WKT
      ,[CodeISO3166_1A3]
FROM [dbo].[Adm0GeoBoundaries]

This is the output:

MULTIPOLYGON (((-62.979628399999974 17.509175600000155,
-62.979600899999951 17.509182200000058,
-62.979600899999888 17.509228199999981,
-62.979575599999933 17.5093026000002,
-62.979671999999788 17.509398900000182,
-62.979789100000083 17.509464600000079,
....... [cut]......
-62.967030600000029 17.464950999999978,
-62.966103599999975 17.464639899999998)))

This output is then converted to a GeoJSON to show it on the web.

As you can see the precision of the coordinates is exaggerated as we have 15 decimal places that correspond to precision to the millimeter and we are talking about showing boundaries over a map!

So it's enough to cut to 4 decimal places corresponding to ~1m of precision.

Reference: http://blog.brainsen.com/post/geographical-coordinate-precision

With the Reduce(tolerance) method I can simplify the geometry with the Douglas-Peucker algorithm and that's fine, but it's completely useless to millimetre precision, that in the end results to double my GeoJSON file.

So I'm wondering if there is an easy way to decrease the precision of the decimal places for each coordinate of a certain amount.

I was thinking also to a dirty solution like a RegEx on the WKT...


texts: Changing lat long to UTM in csv using Save As in QGIS

texts:
I have tried to change the lat long to UTM using "Save As" in QGIS 2.18.13 and tried again in QGIS 3.0.1, I cannot get it to work.

I have a .csv file in the form

  id Base_latitude   Base_longitude
  L001   -33.913912  150.862774
  L002   -33.942172  150.857137
  L003   -34.055889  150.747552

I add it to QGIS with the CRS set to ESPGS:4326 - WGS 84

I then right click on the layer, click "Save As" and set the CRS to ESPG:32756 - WGS 84 / UTM zone 56S and save it as a new .csv.

According to the tutorials and previous stackexchange answers this should produce a .csv file in the form

  X      Y      id   Base_latitude   Base_longitude
  number number L001 -33.913912  150.862774
  number number L002 -33.942172  150.857137
  number number L003 -34.055889  150.747552

Instead the output .csv is exactly the same as the input. Am I missing a step? I've tried in two different QGIS versions (2.18.13 and 3.0.1) so I don't think it's a bug.


texts: Are there other values than "wms" and "WFS" for providerType() for a WMS/WFS layer?

texts:
According to the docs, providerType() returns the provider type (provider key) for a layer.

According to https://gis.stackexchange.com/a/262050/105333, the values for a WMS/WFS layer are wms and WFS.

I found no documentation why wms is in lowercase and WFS is in uppercase. Also there seems to be no list of all existing provider types available.

This raises the question, if there are additionally other values like e.g. WMS or wfs around.


texts: After merging, the file becomes dataframe. Then R does not allow me to use geom_sf

texts:
I am merging attributes of provincical government spending to the China map.
I have a shp file and a csv file.
After inner_join, the class becomes dataframe.
Then R shows this error message: "stat_sf requires the following missing aesthetics: geometry"
How do I change the merged file into geometry?
china_shp <- st_read("C:/Users/wenra/Box Sync/china_map/chn_admbnda_adm1_ocha.shp")
class(china_shp)
str(china_shp$ADM1_EN)

gdp_china<-readr::read_csv("./province_gdp.csv")
map<-inner_join(gdp_china,china_shp)
class(map)
ggplot(map)+
  geom_sf(aes(fill=gdp))


texts: ArcGIS 10.2 Windows 7 only uses 4Gb of RAM

texts: I have 16Gb and have installed 64-bit background processing, but ArcMap seems to only be able to access 4Gb. I can use more in total with other programs. Is there a setting somewhere in Windows or ArcGIS that is only allowing it to access 4GB? Thanks.


texts: Find nearest polygons of a multi-line string

texts:
I am trying to find the nearerst (100-200 meters) parcels of the green line. Given The photo below, they database should return 12 records. I am using ST_DWithin postgis function. Is this the correct way?
FROM
gis_streets as gs --green line
left join "parcel" as p --black polygons
on ST_DWithin(gs.geom,p."parcel_geometry",3000)


texts: How to achieve an "inner glow" effect on polygons drawn using the Google Maps v3 javascript API?

texts:
Is it possible to produce an "inner glow" effect on polygons drawn using the Google Map v3 javascript API?

If not possible using the vanilla Google Maps API, is there a 3rd party tool that can do it? Or perhaps some other technique I could use to accomplish my goal?

The inner glow effect I desire would resemble this:

Thanks!


texts: Select by location using for-loop in Python Terminal

texts:
I have two shapefiles ('drenagem' and 'rio_pri').

I need to loop through every feature from 'rio_pri' and save all those features from 'drenagem' that intersect 'rio_pri'.

So I wrote the following code to run in Python Terminal.

import processing
drenagem = QgsVectorLayer('/home/tiago/Documentos/Rios_mata_disol.shp','Rios_mata_disol','ogr')
rio_pri = QgsVectorLayer('/home/tiago/Documentos/rios_linha_costa.shp','rios_linha_costa','ogr')
# Loop
for r in rio_pri.getFeatures():
    processing.runalg("qgis:selectbylocation", drenagem, r,['intersects'],0,0)
    r_cod = r['COCURSODAG']
    path = r_cod+'.shp'
    processing.runalg('qgis:saveselectedfeatures', drenagem, path)

Then for every feature from 'rio_pri' I get this message at the console.

Unable to execute algorithm
Wrong parameter value: <qgis._core.QgsFeature object at 0x7ffb6620b770>

I guess I should not use the 'r' as an argument of qgis:selectbylocation, but I could not figure out what should I use instead. Any ideas?

I'm using QGIS 2.18.17 on Linux Mint 19.


texts:
I am wanting to join a csv with an shp file any ideas on the methods used to complete this.

I've seen a few different ways none of which are working for me?


texts: Extract the lowest value from an attribute in QGIS

texts:
I am not to sure if this is possible?  What I am looking to do is.

1 - I have polygon and point data. I want to use the polygon to select the point data that falls within in the radius of the polygon.

2- The point data has attributes linked to a common id. I want to then extract the lowest value from that common ID then place it into a field in the polygon.

E.g.  point data looks like

field1 200C063D094B09   field2 BLSK field3 2345
field1 200C063D094B09   field2 CDMN field3 2456
field1 200C063D094B09   field2 GTNG field3 2432

Then I want to add Field 1 and 3 to the polygon.

Is this even possible? I am not good at scripting so this would be manual


texts:
Is there a freeware alternative to Tableau 8? One that allows me to use my own data (shapefiles)?
Tableau is a comercial software that permits you to create a dashboard were you can show at same time maps e diferent kind of graphs. Besides it is very limited for a GIS, it could be very useful to make simpler visual analysis and communicate results in a interactive environment. As my specialty is environmental analysis, like vulnerability analysis, my main interest is to find a tool to communicate the results and that allows anyone to do simple analysis, like correlation graphics (see the link for example http://esl.jrc.ec.europa.eu/). There is some tools in the market but or they are expensive, or they do not permit that you use your own data.

Tableau 8 http://www.tableausoftware.com/new-features/8.0


texts: Performing Ordinary Kriging in SAGA via QGIS SEXTANTE Plug-in?

texts:
I'm trying to learn how to apply the kriging analysis in QGIS using the SEXTANTE plug-in and SAGA tools.  It appears that I have my SEXTANTE and SAGA folders configured properly according to the user guide.  I've never used a kriging function before and I can't find a good resource from SAGA explaining the different input parameters.  Does anyone have experience using the ordinary kriging tool in SAGA that can provide advice or a reference?  I'm currently getting the following error message in QGIS:

An error has occured while executing Python code:
Traceback (most recent call last):
File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\gui\ParametersDialog.py", line 137, in accept if self.setParamValues():
File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\gui\ParametersDialog.py", line 91, in setParamValue if not self.setParamValue(param, self.paramTable.valueItems[param.name]):
File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\gui\ParametersDialog.py", line 129, in setParamValue return param.setValue(widget.getValue())
File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\parameters\ParameterExtent.py", line 11, in setValue if value is None: NameError: global name 'value' is not defined


texts: PostGIS Geocode function taking a lot of time some address, but not others

texts:
I setup a PostGIS/Tiger database and have been going through trying to optimize my queries. At this point I have am having issues with geocoding certain addresses, and I just can't seem figure it out. I wanted to ask for some options to explore.

Some addresses seem to geocode quickly (20%), while others (80%) seem to be incredibly slow (10-20 seconds).

These queries are lightning fast ( < 1 second in many cases). Even if I change the address around a lot, to avoid query caching, it still comes back. Both queries below for example are ones I never ran before in Postres. There are more, but you get the idea.

SELECT rating FROM geocode('1300 Pennsylvania Avenue Northwest,
                            Washington, DC 20500', 1)

SELECT rating FROM geocode('Southold, New York 11971', 1)

If I run these to geocode, however, they all take 20-60 seconds to run.

SELECT rating FROM geocode('296 Ridgewood Ave, Glen Ridge, NJ 07028')

SELECT rating FROM geocode('2520 S Decatur Blvd, Las Vegas, NV 89102')

SELECT rating FROM geocode('Haverford St, Boston, MA 89102')

SELECT rating FROM geocode('1419 Westwood Blvd, Los Angeles, CA', 1)

SELECT rating FROM geocode('Forest Hills, NY 11375', 1)

A few notes / things I have tried:

I have every state loaded into the database (whole thing is 130GB or so) so I understand that adds overhead, but shouldn't it be uniform in its latency?
There isn't a ton of space left on the server after the setup (maybe 1GB or less). Not sure if that would affect performance, especially considering some queries run fast.
I do recall running the command select install_missing_indexes();, and I ran it a few other times but recall getting errors on creating unique indexes (presumably b/c they are already created). When I go to certain tables, they certainly do have indexes.
In many similar queries that ran fast, they were new geocode queries (I never ran them before for that specific address or location). It seems some run faster than others.
It could just be me, but it seems like Washington DC addresses, malformed or otherwise, always come back in < 1 second. Anywhere in California seems to be an issue. NY is mixed results.

I've hit a wall - Could someone offer me some thoughts on what I could investigate at this point?

this is the error I get when I run(rerun)

select install_missing_indexes();

It takes about 20-30 minutes to run, and finally spits this out:

geodb=# select install_missing_indexes();

ERROR:  could not create unique index "uidx_tiger_data_dc_faces_tfid"

DETAIL:  Key (tfid)=(210417114) is duplicated.

CREATE UNIQUE INDEX uidx_tiger_data_wy_faces_tfid ON tiger_data.wy_faces
USING btree(tfid);";tree(soundex(place));) varchar_pattern_ops););ty) varchar_pattern_ops);

PL/pgSQL function install_missing_indexes() line 4 at EXECUTE

@lr1234567 identified the issue was that my indexes were not installed and I had to reload specific tables to get it working.

Once the indexes loaded, however, the tracts were not coming in. I suspect something was broken, or there was a bug, or maybe I didn't do something right, but regardless the tracts would not load.

I ultimately created an AWS public image from the Ansible Playbook, and then did a postgres data dump of all the tract tables for all 50 states, and then imported them into my solution.

After that things were better, and my queries are all under 5 seconds (there are a few random ones that take 18 seconds, but for the most part all is well)


texts: Building single sided buffer around categorized polygon layer

texts:
I need to build a 100m buffer around a categorized polygon layer using QGIS 3.2.

The layer has four classes and no holes. The difficulty here is, that the buffer needs to be split into the 4 categories at the edges of the original polygons. Hard to explain in words, so I drew a quick sketch:

It doesnt have to follow a specific rule and it can be a multistep solution as well.

How can this be done using QGIS 3.2? And is there a special term for single sided buffers like this?


texts:
Able to connect to  Oracle Spatial 11g and print tables with SDO_GEORASTER attribute both in QGIS python console and standalone application.I get below error when loading rasters from oracle Table.

  Traceback (most recent call last):   File "", line 1, in
   Exception: unknown

    strProjetName = "C:/OSGeo4W/apache/htdocs/QGIS-Web-Client-master/projects/myworld.qgs"
    def add_LayersFromDB():
        QGISAPP = QgsApplication(sys.argv, True)
        QgsApplication.setPrefixPath(r"C:\OSGeo4W\apps\qgis", True)
        QgsApplication.initQgis()

        #Connect to Oracle and Fetch Table Names
        con = cx_Oracle.connect('myoradb/myoradb@189.53.61.140:1521/newdev')
        print con.version
        cur = con.cursor()
        cur.execute(u"select TABLE_NAME from user_tab_columns where data_type='SDO_GEORASTER'")

        tables = cur.fetchall()
        print tables

        QgsProject.instance().setFileName(strProjetName)
        print QgsProject.instance().fileName()
        render = QgsMapRenderer()

        uri = QgsDataSourceURI()
        uri.setConnection("189.53.61.140", "1521", "newdev", "myoradb", "myoradb")

        render = QgsMapRenderer()

        for t in tables :
           print str(t[0])
           uri.setDataSource('',  str(t[0]), "RAST")
           uri.uri()
           vlayer = QgsRasterLayer(uri.uri(),str(t[0]), 'oracle') #FAILS HERE
           QgsMapLayerRegistry.instance().addMapLayer(vlayer)
           lst = [ vlayer.id() ]
           render.setLayerSet(lst)
           rect = QgsRectangle(render.fullExtent())
           rect.scale(1.1)
           render.setExtent(rect)

        QgsProject.instance().write()
        cur.close()
        con.close()
        QgsApplication.exitQgis()

    add_LayersFromDB()


texts:
I want to "put together" single objects of a layer (shapefile) in order to create a new layer out of the already selected objects.

In my case, I work with a shapefile of Mexico that shows its States (estados) and I want to select a few states, make a new layer out of those states and visualize them in another colour.


texts: Find a path between an origin and destination given scattered points

texts:
We have a postgis database containing a lot of ordered (by date/time) nautical point (lat/long). Each of these points are affected to a leg between an origin and a destination (both origin and destination are on a coast).

Based on this we are trying to find every possible paths (hopefully only one in most case) for each origin/destination combination (several legs for a given combination), do anyone know if there already exists some sort of algorithm/function that does this ?

EDIT: "affected to a leg" means that we can map every point we have to a combination origin-destination (so for instance we know that Px (-122°, 85°) has been passed by a vessel doing the voyage between origin A and destination B). What we are tying to do is find the possible geographical path between A and B based on the positions we know have been visited by vessels doing this voyage.


texts: Transforming GRIB file into GeoTIFF that is correctly georeferenced using Python GDAL?

texts:
I am trying to transform a GRIB file into a GeoTIFF to be used in a ArcGIS, but am having trouble getting the image to display in the correct location. I have been able to create a GeoTIFF, using GDAL in Python, that shows the data but is not showing up in the correct location when brought into ArcGIS. The resulting image is below.

The data I am working with can be downloaded from: https://gimms.gsfc.nasa.gov/SMOS/SMAP/L05/

I am trying to project the data into WGS84 Web Mercator (Auxiliary Sphere), EPSG: 3857

Note: I have tried bringing in the data via ArcMap by creating a Raster Mosaic which should be able to work with GRIB data, but I didn't have any luck. I have also tried using the Project Raster tool, but ArcGIS does not like the default projection that comes from the GRIB file and gives an error.

The code I'm using:

import gdal

src_filename = r"C:\att\project\arcshare\public\disaster_response\nrt_products\smap\20150402_20150404_anom1.grib"
dst_filename = r"C:\att\project\arcshare\public\disaster_response\nrt_products\smap\smap_py_test1.tif"

#Open existing dataset
src_ds = gdal.Open(src_filename)

#Open output format driver, see gdal_translate --formats for list
format = "GTiff"
driver = gdal.GetDriverByName( format )

#Output to new format
dst_file = driver.CreateCopy( dst_filename, src_ds, 0 )

#Properly close the datasets to flush to disk
dst_ds = None
src_ds = None

I am not very well versed in using GDAL or GDAL in Python.


texts:
How to match two raster layers, I have one raster and a shapefile and I want to convert the shapefile into raster BUT I want it to have the same cellsize and its cells overlay (same positions) to the first raster.
I already set as a 'mask' the first raster in the Environments options, but the resulting raster cells doesn't overlay with the first raster.


texts: Collecting multipolygons into a single multipolygon with ST_Collect() in PostGIS

texts:
I would like to collect 50,000 individual multipolygons into one big multipolygon. Unfortunately when I perform ST_Collect(geom) it bundles them into a geometry collection and cannot be cast as ::geometry(MultiPolygonZ, 2193) which is of no use to me.

How do I get prevent this from happening? I have tried wrapping the ST_Collect() in a ST_CollectionExtract() but that merely returns the original 50,000 multipolygons.

Here is the code:

CREATE TABLE new_table AS
SELECT
  1::INT as id,
  ST_Collect(geom)::geometry(MultiPolygonZ,2193)
FROM
  old_table


texts: Disconnected Satellite map tiles usable with OpenLayers or other JavaScript API

texts:
I'm looking for satellite map, usable with OpenLayers (or with another JavaScript API), with a good scale (at least 100m, MapQuest stops at 1000m) and available for a disconnected mode (available for download).

If I'm right Google forbids to download the tiles and only allows you to keep it in the browser cache, which in my case isn't efficient enough or too risky if the user clears his browser cache.

I'm not really good with terms of service, if possible something free or not too expensive.


texts: GeoKettle Spatial Reference Transformations

texts:
I'm attempting to use GeoKettle to transform some Ordance Survey shapefiles from British National Grid projection into WGS84 EPSG:4326; however, there is a discrepancy in the location of approx 80m between the output from GeoKettle compared to the same data that has been reprojected in ArcGIS or QGIS.

In the SRS Transformation dialog the 'source spatial reference system' is correctly detected as British National Grid (EPSG 27700) and then I select EPSG 4326 from the right-hand pane so I don't see anything that I could be doing wrong - is this a known error?


texts: Is there a DTM available for free somewhere?

texts:
To go a little bit further than the semantic question regarding DEM/DTM/DSM, I'd like to know if there is a free (of charges and use) DTM datasets? And what about resolution?

As far as I know,only DSM are available to public for free.


texts:
Rough English (Translated by a user):

  When trying to access plugins I am told Unable to open the Python
  console.  When I start qgis, the message also shows up but I don't know what it means.
  The problem is that I have no access to the Add Delimited Text Layer plugin.

Original Spanish:

  al intentar acceder a complementos me dice que No se pudo abrir la
  consola de Python .
  al iniciar qgis tambien aparece el mensaje pero lo ignoro .
  el problema esta en que no puedo acceder a complementos de Add Delimited
  Text Layer.
  ayuda por favor.


texts: GeoDjango geometry transform method is swapping x with y

texts:
I need to convert a geometry from srid 4326 to srid 3857:
myPoint4326 = Point(x=-73,y=40, srid=4326)
myPoint3857 = myPoint4326.transform(3857,clone=True)
print(myPoint4326.wkt)
print(myPoint3857.wkt)

The transform method seem to swap my point x with y:
POINT (-73 40)
POINT (4452779.631730943 -12123477.91689172)

When I draw the "POINT (4452779.631730943 -12123477.91689172)" on https://maththinking.com/gisvisualizer/# , it is not the same location as POINT (-73 40).
Am I missing something?


texts: What are offline/client geometry capabilities of ArcGIS API for JavaScript?

texts:
Lately I've been going through ArcGIS JS API samples, but from what I've seen, it looks like all geometric functions always go through server, and use server side REST functions. There is no option to even calculate area of polygon using only JS API without server side funcions.

I have a case, where I have feature layer added from JSON file and I want to use some geometric functions un features inside. Server connection is not an option, everything needs to be done on client side.

Is there a way to do this with ArcGIS JavaScript API only, without using server side tasks or maybe some other alternative frameworks, that could use features from JSON file?


texts:
I have query by distance working in a framework that generates a list of features within a fixed distance of a clicked feature. I would like to make the queried distance variable based on either:

Pixel Distance from clicked feature (something like):

select cartodb_id, label, title from full_list where st_distance( the_geom, st_GeomFromText('POINT("+lon+" "+lat+")', 4326), true ) < ("+pixels+")

Or current zoom level-based resolution (something like):

select cartodb_id, label, title from full_list where st_distance( the_geom, st_GeomFromText('POINT("+lon+" "+lat+")', 4326), true ) < (SELECT CDB_XYZ_Resolution("+zoom+"))

I wonder if this is possible. Here is the jsbin wherein CDB_XYZ_Resolution uses a fixed zoom variable: http://jsbin.com/vifaleri/1/edit

To update the distance variable when a feature is clicked, I'm not quite sure where to include something like this:

map.on('zoomend', function(){
    var zoom = map.getZoom();
    console.log(zoom);
  });

to be reachable by the CDB_XYZ_Resolution("+zoom+") part of the feature click query.

Thanks and apologies for my ignorance. If grids are the solution, I have been unable to implement grid queries of points in the same cell as the selected point. I've tried this, among other similar grid queries that return no records:

WITH hgrid AS (SELECT CDB_HexagonGrid(ST_Expand(!bbox!, greatest(!pixel_width!,!pixel_height!) * 15), greatest(!pixel_width!,!pixel_height!) * 15) as cell) SELECT hgrid.cell as the_geom_webmercator, 1 as cartodb_id, count(i.cartodb_id) as point_count FROM hgrid, full_list i WHERE i.cartodb = "+data.cartodb_id+" AND ST_Intersects(hgrid.cell, i.the_geom_webmercator)

Thanks.


texts: Possible error when adding objects date atributte in ArcGIS rest

texts:
I have published in ArcGIS server 10.04 (on win 2008 r2 64 bit + sde + oracle) table ACCESS_LOG on FeatureServer with esriFieldTypeDate field. Then i try to add record with /addFeatures:

[
  {
    "attributes" : {
      "WHO" : "127.0.0.1",
      "WHEN" : 1320859183000,
      "WHAT" : "topoMap"
    }
  }
]

and it work as expected - the record is in rest. But when I change date attribute from "WHEN" to "when" (set it lower case) - i get an error:

'features' param is invalid

It gets better: when I set all parameters lower case and "WHEN" upper case - the record can be added.

I've tested and encountered this behavior in ArcGIS server 10.2, 10.4 on oracle, postgres, mssql. For the record i check in arcgis sample servers and couldn't reproduce the error.

Did you encountered similar problem?
Could you in your infrastructure create simple table, try to reproduce the problem and give the feedback?


texts:
Possible Duplicate:
  Bulk load multiple shapefiles into PostGIS

Is it possible to bulk upload shp files into postgis. Im currently playing around with postgis, qgis and geoserver (as a potential arcgis server, arcmap and arcsde replacements) and want i want to be able to do is run a batch file (I guess) that will upload a number of shp files (exported from arcsde) that overwrite existing shp/tables files in postgis.

In effect postgis would be a slave database (for the time being).


texts: How can I zoomToExtent and then zoom out?

texts:
I use openlayers 2.12

I use HTML5's geolocation API to get the coords and then I insert them in a layer as feature (point). After that, I zoom in , since the layer has no other feature.

map.zoomToExtent(geoloc.getDataExtent());

geoloc is the name of the layer...

My problem is, how do I zoom-out automatically after the zoomToExtent?

Because the map zooms in too much and I can see just the point and very little details. It would be more user-friendly to zoom out a bit

Thanks


texts: Loading QGIS external processing tools into Python environment

texts:
I'm looking to run some QGIS algorithms in a separate python environment (not the one built in to QGIS) and would like to use some externally provided processing tools (ie. FUSION, LAStools, UMEP etc). I've been able to get an environment working (I'm using Python 3 in Spyder through Anaconda, and QGIS 3.16) and can access the built in processing tools, but can't figure out how to add these additional tools.
Thoughts?


texts: Count number of features that matches a (minimum) value for each group (id)

texts:
Context
In FME I want to count for each group (defined by a "group by" id in any kind of aggregation-like transformer, such as the StatisticsCalculator), how many occurrences equal to a certain value can be found in a field (aka attribute).
In my situation, I'd like to figure out how many values in a field named hour matches the minimum and 06 (in most of the records, (i.e. rows), but not necessarily for all, 06 is the minimum) for each id group.
The StatisticsCalculator can compute the minimum for each group, but I'm not able to figure out how to actually count how many values match this, neither how many values match a given one, e.g. 06`.
Sample
group_id hour
       1    6
       1    6
       1    6
       1    7
       2    7
       2    7
       2    7
       2    9
       3    6
       3    6
       3    7
       3    8

Desired outputs
This table:
group_id hour_min count
       1        6     3
       2        7     3
       3        6     2

and this one:
group_id hour_eq6_count
       1              3
       2              0
       3              2

or this single one:
group_id hour_eq6_count hour_min hour_min_count
       1              3        6              3
       2              0        7              3
       3              2        6              2


texts:
According to the Leaflet docs the marker's autoPan option works like this:

  Set it to true if you want the map to do panning animation when marker
  hits the edges.

What does 'the marker hits the edges' mean?

When I set the option to true and add a new marker right on the edge or outside of the current view, the map does not automatically move.


texts: Can a cached map service automatically Track feature class changes and update only the changed locations

texts:
Is there any option I can enable in an arcgis 10.2 cached map service, so as to automatically update the tiles for the locations that have been changed ?

If it is not that simple, I suppose I will have to:

  Target your cache updates to cover only the areas where data was edited (as described here)

and as referred here to :

  create a custom tool to track all the changes since the last update, then export these changed areas to a feature class and update tiles based on the feature class boundary.

However how could I track the changes in the feature class ?


texts: How do I visualize a many-to-many relationship?

texts:
Preferred software use: ArcGIS Pro and ArcGIS Online

For a client I have to visualize a building plan that has a 'many-to-many'-link with attribute data. As with many-to-many relationships; sometimes an attribute applies to multiple rooms, and rooms themselves are in most cases linked to multiple attributes.

For this I want to create an attractive visualization. Preferably a bit like this:

The screen is split in half, the left is the list of attributes. Clicking on one attribute will make a change in the map on the right, where only the rooms will lighten up that apply to this attribute. Clicking in the map will make a selection on the list on the left side. Showing only the affected attributes for that specific room.

I'm not sure if this is possible in ArcGIS Online WebAppDeveloper at the moment. There are many widgets, but my knowledge about these tools is currently very limited.


texts: GDAL Raster: How to obtain the cell ( and also the position) where maximum occurs?

texts:
I can obtain the maximum/minimum value of a raster via the ComputeRasterMinMax code:

  using (var ds = Gdal.Open(tiffFile, Access.GA_ReadOnly))
  {
    var nCol = ds.RasterXSize;
    var nRow = ds.RasterYSize;
    var allValues = new int[nCol * nRow];
    Band band = ds.GetRasterBand(1);
    band.ReadRaster(0, 0, nCol, nRow, allValues, nCol, nRow, 0, 0);
    double[] minMax = new double[2];
    band.ComputeRasterMinMax(minMax, 0);

  }

As per documentation, minMax[0] is the minimum value, and minMax[1] is the maximum value of the band.

But how to get the cell that these maximum and minimum occurs, besides actually looping through the array and compare the item one by one?


texts:
I want to see the coordinates of the standard WMTS grid at various zoom levels so that I can see their "coordinates" as well as their extents. I am pretty sure that I had used a plugin for this at one point but it was a long time ago and I cannot find it again.
Whenever I try to search for help on this I get a deluge of questions about WMTS layers in general, but what I really want is just a grid with the coordinates and zoom levels that I can overlay on the canvas
Is there a plugin for this, or at least a vector dataset containing polygons with this information embedded in them?


texts: Collect results of iterative processes in QGIS graphical modeler

texts:
I'm using QGIS graphical modeler to build an iterative process that is supposed to:

extract features from vector layer 1 (VL1) overlapping with features from vector layer 2 (VL2) (iterative per features in VL2)
then process these extractions individually
then collect/merge the separate outputs for further processing (e.g. merge all outputs to a single layer, then proceed with other algorithms)

How do I realize step 3 in the graphical model? Is that even possible?
Step 1 and 2 work fine as an "iterative process" when using Vector Features as input (instead of Vector Layer), as was suggested here: When running the model I can tap the "iterative" button:

But subsequently using the "Algorithm output" will always be just for the output of one "feature-iteration" (generating an output for each feature in VL2), but not the collection of all of them (generating one output).
Below is a picture of the model I need it for: Extract by polygon (yellow circle) has to be applied per feature (not per layer). The "Merge vector layers" (pink circle) merges the output of the various algorithm outputs, but not per feature iteration.


texts: Adding image watermark to map via ArcGIS API for JavaScript?

texts:
i need to use a png as a watermark for a map that is generated via javascript api. Don't know how to do that.

I'm using ArcGIS 10.1


texts: What is None argument in sql clause in UpdateCursor?

texts:
I have a bit of code that is a groupby iterator.  It works well, and I understand most of it, but I am wondering what the None argument does within the sql_clause in the code below.  I know this is probably a really simple explanation and I have seen this in other SQL clauses as well, but have not been able to find a good explanation of why it is needed.  Below is my snippet of code:

case_fields = ["unique_id", "ELCODE", "SeasonCode"]
max_field = "PERCENTAGE"
sql_orderby = "ORDER BY {}, {} DESC".format(",".join(case_fields), max_field)

with arcpy.da.UpdateCursor(intersect, "*", sql_clause=(None, sql_orderby)) as cursor:
    case_func = itemgetter(*(cursor.fields.index(fld) for fld in case_fields))
    for key, group in groupby(cursor, case_func):
        next(group)
        for extra in group:
            cursor.deleteRow()

Just looking for an explanation or a source that describes what the None argument does within the UpdateCursor where clause.


texts: Loop checking if attribute field exists in PyQGIS

texts:
I'm building my first script and I need help to create a loop instead of my repetitiv functions. The aim is to check if every field already exists, and to build it only if it doesn't.

This script does what I want:

layer = iface.activeLayer()
 dpr = layer.dataProvider()
if layer.fields().indexFromName("apr") == -1:
    dpr.addAttributes([QgsField("apr", QVariant.String)])
if layer.fields().indexFromName("iso") == -1:
    dpr.addAttributes([QgsField("iso", QVariant.String)])
if layer.fields().indexFromName("exp") == -1:
    dpr.addAttributes([QgsField("exp", QVariant.String)])
if layer.fields().indexFromName("foca") == -1:
    dpr.addAttributes([QgsField("foca", QVariant.String)])
if layer.fields().indexFromName("foca35") == -1:
    dpr.addAttributes([QgsField("foca35", QVariant.Double)])
if layer.fields().indexFromName("shsp") == -1:
    dpr.addAttributes([QgsField("shsp", QVariant.Double)])
if layer.fields().indexFromName("fyd") == -1:
    dpr.addAttributes([QgsField("fyd", QVariant.Double)])
if layer.fields().indexFromName("gyd") == -1:
    dpr.addAttributes([QgsField("gyd", QVariant.Double)])
if layer.fields().indexFromName("ra") == -1:
    dpr.addAttributes([QgsField("ra", QVariant.Double)])

if layer.fields().indexFromName("pi") == -1:
    dpr.addAttributes([QgsField("pi", QVariant.Double)])

Here's my first attempt to convert it into a loop. This code doesn't work because of synthax, probably this partQVariant.fields_arr[i][1]

Also, I would like to know if there is another way to optimize it? The loop's action performs addAttributes every time.
Could not we optimize that by using addAttributes function only one time after the loop?

layer = iface.activeLayer()
 dpr = layer.dataProvider()
fields_arr =  [["apr", "Double"],
                ["iso", "Double"],
                ["exp", "String"],
                ["foca", "Double"],
                ["foca35", "Double"],
                ["shsp", "Double"],
                ["fyd", "Double"],
                ["gyd", "Double"],
                ["ra", "Double"],
                ["pi", "Double"]]
for i in range(0,9):
    if layer.fields().indexFromName(fields_arr[i][0]) == -1:
        dpr.addAttributes([QgsField(fields_arr[i][0], QVariant.fields_arr[i][1])])


texts: Automating plan report generation using ArcGIS Desktop?

texts:
i have shapefile that contains  a lot of parcels "parcel has like fields (n° parcelle,proprietaire, cin) see picture :

and i want for all parcels to generate this plan :

take a parcel and generate this. (red text is clarification)

software or language programming not important, i'm open to any idea or answer that can solve this problem(product presentation format for all parcels with the format in picture 2 ).


texts: Handling warning about DX and DY in gdal_translate from .Gtiff to .ASC?

texts:
When I try to convert Gtiff to .asc file, QGIS shows the warning:

Producing Golden Surfer style file with DX and DY instead
of CELL SIZE since the input pixels are non-square. Use the
FORCE CELL SIZE = TRUE creation option to force use of DX for
even though this will be distorted. Most ASCII Grid reader
(ArcGIS included) do not support the DX and DY parameters

Can I solve this problem?


texts: Can polyline shapefile hold geometries with 0 vertices?

texts:
I have a polyline shapefile and a few geometries (5 out of 2000+) have no vertices at all.

I am a bit surprised and was wondering whether my shapefile is broken (so I have a bug when creating the shapefile) or my shapefile reader is broken (and does not consider for shapes with no vertices)?

I am actually reading the shapefile through a Python library (shapelib) trying to extract some information.


texts: Extracting area in orthophoto high resolution image tab file?

texts:
How can I extract an area in an orthophoto high resolution image tab file?

I tried to translate it in MapInfo to other data files but cannot.


texts: Export table structure information to Excel QGIS or MapInfo

texts:
I want to create a data dictionary of the tables of our primary layers. So I want a way to export the table structure to Excel. I don't need the actual data, I just want to know what columns are in each table and the data types for each column.

Our layers are in MapInfo, and I also have access to QGIS.

I can use the MapInfo Table Manager to export to a .met file, but when I copen that with Excel, the information is imported into one column - then requiring manual processing.


texts: How to add GDAL output formats support to WCS requests in GeoServer?

texts:
I am running GeoServer 2.18.1 and I installed GDAL extension as well as GDAL native libraries (I tried several versions with the same result).
I followed this guide : https://geoserver.geo-solutions.it/downloads/releases/2.8.x-ld/doc/community/gdal/index.html.
I manage to get the GDAL formats in the Raster data sources so I assume GDAL extension is correctly installed.
I've set all the necessary environment variables.
But no new format appears in the ServiceMetadata section of the WCS 2.0 GetCapabilities document.
Is there something I am missing?


texts: PostGis ST_ConcaveHull Zero Longitude Error


texts: How to create a line by passing two points and save in feature class in arcobjects using vb.net

texts:
I want to create a line from two points in acrobjects. but i don't know how to save that line to a feature class

public static ILine CreateSegment(IPoint pFrom, IPoint pTo)
 {
    ILine line = new LineClass()
    line.PutCoords(pFrom,pTo)
    return line as ILine;
 }


texts: What is the Difference between Geometric and Geographic columns?

texts: I am a beginner in GIS and Learning PostGIS. I came across this two terms Geom and Geog, Could anyone explain in brief. What is the best to use with map data and Do all functions support both Geom and geog.


texts: Where can I can find a UTM 10 x 10 km grid?

texts: I am looking for a UTM 10x10km grid map for Portugal and I can’t find it. Does somebody know where I can find this shapefile?


texts:
So reading this GRASS/Shell documentation I decided to test my own for repeat using this lines

CATS="(v.db.select -c map=training columns=cat)"
echo $CATS
for i in $CATS
do
    where="cat = $CATS"
    outputname=region_uspo_$CATS
    regionname=subset_uspo_$CATS

v.extract -o input=training type=area where=$condition output=$outputname
done

but I keep getting an error message with v.extract

  ERROR: v.extract: Sorry,  is not a valid flag
  ERROR: v.extract: Sorry,  is not a valid flag
  ERROR: v.extract: Sorry,  is not a valid parameter
  ERROR: v.extract: Sorry,  is not a valid parameter

now I've also tried this
CATS="$(v.db.select -c map=training columns=cat)"

Can anyone help me figure out where I did wrong?


texts: ArcGIS Server Manager keeps getting COM errors

texts:
I am using ArcGIS Server 10 on IIS6.  While publishing out mapping services/restarting services/deleting services, I keep getting various COM errors like the one below..after viewing the Event log it looks like the ServiceWizard.aspx file is generally the root cause.

However, usually if I wait thirty minutes or go to the restroom or other seemingly non-related activity and then try again it that I was trying to do will work.

The errors are a little to scattered to describe fully.  The question is..Is there something I may be doing wrong or should I restart the SOM every few days or another option?  FYI, there are a lot of mapping services on this instance.  Thanks.

"An error has occured on the server. For details please check the Event (Application) log on the web server.
Error HRESULT E_FAIL has been returned from a call to a COM component."


texts:
I want to label my polygons with information, and I want the labels to be in two rows - name on row 1 och volume on row 2.
My problem is, that I get an extra "quare" at the end of row 1?
That has never happened before.

EDIT: Here is the expression I wroght:

And, it works, I just get that extra quare, like this:

I use QGIS 2.10 and Windows 10, could it be because of the Windows 10? I never had this problem on Windows 8.

EDIT 2: I realize that I need to clarify my question, I have two different columns. I want the information in column 1 to be in the first row and the information from column 2 to be in the second row of the labels.
But, my expression makes the labels with an extra square at the end of row 1.

I use encoding windows-1252.


texts: How Can I Add a Summary Row to Graphic Table Element Using ESRI Production Mapping Extension?

texts: I currently have a map book (data driven pages) with 160 pages, and counting, consisting of domestic water pipes for a city. I would like to create a Graphic Table Element with a summary row totaling the pipe lengths on each page.  I have the ESRI Production Mapping extension.  How can I add a summary row to the Graphic Table Element?


texts: Calculating coefficient of variation from gstat output to create uncertainty map in R

texts:
I was asked by my supervisor to create uncertainty maps of kriging interpolation based on the coefficient of variation ((sd/mean)*100%). I used krige function from gstat package to perform the interpolation.

#interpolation using kriging with external drift
krig1 <- krige(xSO4.2009. ~ easting+lograin.2009., dat1, grid.uk, model=fitvar1)

The output produce the prediction values (var1.pred) and the prediction variances (var1.var).

If I want to create the uncertainty map based on the coefficient of variation (COV), does it mean that I just need to use the prediction values and the variance to calculate COV then map it?

pred<-krig1@data$var1.pred
var<-krig1@data$var1.var
krig1$cov<-(sqrt(var)/mean(pred))*100


texts: DIY easy and ~accurate small-scale landscape survey

texts:
Summary
Is there any equipment that I can rent 'inexpensively' (less than $400) that would allow 1-2 people to accurately (+/- 6") capture the topography of hilly terrain of perhaps 250' on a side in an afternoon?
Details
My wife and I are about to embark upon a redesign of the landscape around our home in the mountains. Our lot started out on a generally 1 in 3 slope, and this has (of course) been changed in many interesting ways during building. There are small hills, patios, semi-flattened areas followed by large falloffs.
As part of our designing, we want to capture exactly what the landscape looks like now, digitally. (Then we'll sculpt in the computer what we want, taking into account existing boundaries and surrounding contours, and see roughly how much dirt we need or will need to remove.)
We performed a ghetto survey this weekend with strings and poles, roughly sampling a grid of elevations every 10'. The results we got produced a reasonable-looking representation of the terrain. However, the 10' sampling was a bit too coarse for my tastes and—more importantly—we were only able to sample about 1/4 of the area we need before our measuring poles weren't tall enough or we ran out of reference objects to attach to.
So: what hardware exists that I can rent to do this task for me? Ideally I'd like to take a pole, walk all over the property and every few feet stick it on the ground, press a button, hear it go "beep!" and have it record an XYZ sample that is accurate within +/- 6". After doing this for a couple hours I'd plug it into the computer and download some DEM, CSV, or other accessible data of all the points in a non-proprietary format.
It looks like Trimble offers some neat gadgets in this arena, and I see that some companies rent this sort of gear out, but not being a landscape survey engineer professional I cannot figure out exactly which, if any, products meet my needs. Can someone who knows this product space please help suggest what I should be looking for? Thanks!
Note: I originally asked this question on Home Improvement, but as suggested there and approved here on meta I'm cross-posting it to this site.


texts: How to cycle thru items in print composer?

texts: In print composer I have many items placed on the "paper".  These items include borders, one map, bar scale, logo and an arrow.  When the items overlap each other I am only able to select the item on top.  If I want the item below I have to select the first item and send it to the back.  Is there an easier way to select the item below another, without having to send the first item to the back?


texts: How to Embed a table into QGIS Graphic Modeler?

texts:
I am automating a process in QGIS. So far, I have built the following model.

The HSG Lookup.csv is a user input right now. But I want it to be a predefined table(or csv file) embedded into the model and not the user input.

Is there a way to do that within the modeler?


texts:
I am new to using QGIS and I am currently trying to put roads onto my canvas.

I have tried adding only 'highway' and 'oneway' polyline tags from the OSM database, but there is still a few roads missing (see picture, circled in orange)

I kind of assumed high way will cover all the roads (apparently not!)

Is there an efficient way to find out which tag I have missed? (instead of trial and error with all the polyline tags)


texts: Using IPython if set ArcGIS's python as default kernel in spyder?

texts:
i install anaconda first then set the arcgis python 2.7.8 win32 as interpreter just like
Importing arcpy in Spyder said,but why i cannot open my IPython then?


texts: Why is GDALPolygonize so much slower than ArcGIS Raster to Polygon?

texts:
I am attempting to polygonize a raster using GDALPolygonize() in a Python script.  The script began polygonizing yesterday at 5pm and is still polygonizing now at 9:30 am.  I have no clue how far along it is, but I know it is still going because when I refresh my Windows explorer I can see the file size change for the output shapefile.

My raster is rather large, but I still don't expect it to take this long.  My raster is 35,486 Columns and 23,682 Rows with a 1 meter cell size.  It is a binary raster where a value of 1 represents data and 0 is NoData.

When I polygonized in ArcGIS using Raster to Polygon in the Conversion Toolbox it took 56 seconds.  The resulting shapefile is 200mb while the shapefile still being created by GDALPolygonize is still only 100mb.  That makes me think GDAL is about half way done after running all night.

Specs:
Windows 7 64bit,
8gb RAM,
GDAL 1.10 64bit,
ArcGIS Desktop 10.2,
64bit Background Geoprocessing for ArcGIS Desktop,
Python 2.7.3 64bit

UPDATE
Day 2 - GDALPolygonize is still running.  It has gone overnight 2 nights in a row and through a whole day without completing.  ArcGIS took 56 seconds.


texts: Can you both export and use a variable in the same GEE script?

texts:
Let's say I have too files, library, and main

content of library:

exports.someVariable = 5;

content of main:

var library = require('path/to/library')
print(library.someVariable) // correctly prints 5

But let's say I want to also use someVariable in library, maybe something like this:

content of library:

exports.someVariable = 5;
exports.someOtherVariable = someVariable+3; //will say it doesn't exist in this scope

What I tried
content of library:

var someVariable = 5;
var someOtherVariable = someVariable+3; //it works so far
exports.someVariable;
exports.someOtherVariable;

content of main:

var library = require('path/to/library')
print(library.someVariable) // undefined
print(library.someOtherVariable) // undefined

Any suggestions?


texts: Geometry from wkt/polygon using geotools

texts: I have polygon string ,how can I find geometry by giving polygon as input,Is there any function like we have ST_GeometryFromText('POLYGON ....) in postgis ,which give geometry.


texts: Creating custom line symbol with text in it using ArcGIS Desktop?

texts:
I have some layers for which I have to use special symbology. Most of these symbols are lines (continous, dashed, dotted, etc.) but they have different text in them for example:

----GAS----GAS----GAS----

Can I create a custom line symbol with text in it like in AutoCAD?

So far I've created a line symbol with multiple symbol layers. One Cartographic Line Symbol layer for the line and several Marker Line Symbol layers for each character.
However, it's very complicated to setup the templates, the offsets and so on for each layer. And I have dozens of layers with such line symbol, so it's really time-killing and hard to document properly.

Is there any simplier way to create such linetype?

I'm using ArcGIS for Desktop 10.3 with Advanced license.

Examples:

How I want the lines look like:

Some layers I have to style:


texts: .carto file extension is not allowed in CARTO

texts:
I want to import a .carto file but I get the following error message: unfortunately this file extension is not allowed.

I didn't make this file. Do you know if there is any tool that would help me understand what's the problem with it? Unless it's a problem from CARTO as they say they don't recognize there own file extension.

I have a trial subscription.

Update

I tried BeaR's answer: to zip it. It was uploading, it lead me to the layer dashboard before giving me a new error ...

It doesn't seem to be a new issue as they already talked about it in Github in 2017. I don't know if the answer was solved but even when renaming my file with only lowercase I wasn't able to load it.


texts: Plotting Bare Soil Line on landsat data in RStudio?

texts:
I am using RStudio with landsat package for computing and plotting the Bare Soil Line (BSL) on landsat tif image running on my HP Pavilion i7 laptop. The codes for running the BSL() function is as under:

b3 <- as.data.frame(raster("LT51500371993170ISP00_B3.TIF", na.rm=TRUE))
b4 <- as.data.frame(raster("LT51500371993170ISP00_B4.TIF", na.rm=TRUE))

bsl.1993 <- BSL(b3, b4, "quantile", ulimit=0.99, llimit=0.005)

I have received the error message when I write the code without na.rm=TRUE as by default argument is FALSE for na.rm. After using na.rm=TRUE I received the same message as under. But one thing I want to mention here although my system is 64 bit machine with 8GB RAM and running R-studio 32 bit version when I execute the code for whole tile of landsat RAM reaches up to approximatively 7GB, it gives the message to running the memory low message after 5 minutes of the execution then the same message appeared the second time as under:

Error in quantile.default(ratio43, llimit) :
 missing values and NaN's not allowed if 'na.rm' is FALSE

Please suggest me the solution to solve this matter and how do I improve the low memory. I have also running the linux machine on VMware Workstation.


texts: Variance interval charting in Google Earth Engine

texts:
How do I create an annual MODIS NDVI time series chart by day-of-year (DOY) that shows the interquartile range as an interval envelope with the current year's observations overlayed?

I want to put this year's observations in the context of the expected range - something similar to ggplot2's geom_ribbon() function:

I've not seen an example for how to make an interval envelope chart with Earth Engine. The closest example I have seen is to simply plot the time series of all years as separate lines using the doySeriesByYear() function.

// Define area of interest and display on map.
var aoi = ee.Geometry.Polygon(
        [[[-122.10333508394413, 36.980921854180906],
          [-122.10333508394413, 36.962269988961935],
          [-122.06179303072147, 36.962269988961935],
          [-122.06179303072147, 36.980921854180906]]], null, false);

Map.centerObject(aoi, 13);
Map.addLayer(aoi, {color: 'white'}, 'AOI');
Map.setOptions('SATELLITE');

// Load MODIS vegetation index data and subset a decade of NDVI images.
var band = 'NDVI';
var ndviCol = ee.ImageCollection('MODIS/006/MOD13A1')
                  .filter(ee.Filter.date('2010-01-01', '2020-07-01'))
                  .select(band);

// Plot NDVI as an annual time series, year over year; emphasize 2020 with a
// thicker, darker line.
print(ui.Chart.image
          .doySeriesByYear({
            imageCollection: ndviCol,
            bandName: band,
            region: aoi,
            regionReducer: ee.Reducer.mean(),
            scale: 500,
            sameDayReducer: ee.Reducer.mean()
          })
          .setOptions({
            hAxis: {
              title: 'Day of year',
              titleTextStyle: {italic: false, bold: true}
            },
            vAxis: {
              title: 'NDVI (x1e4)',
              titleTextStyle: {italic: false, bold: true},
              viewWindow: {min: 4000, max: 8000},
              baseline: 0
            },
            lineWidth: 0.5,
            colors: ['#4caf50'],
            series: {10: {lineWidth: 5, color: '#087f23'}},
            legend: {position: 'none'},
          }));


texts:
I am trying to build QGIS from source on x64 Windows machine.

QGIS: qgis2.99.0
VS: 2015

I downloaded latest cygwin, OSGeo4W etc. Now, CMake gives an error:

Looking for gaiaStatisticsInvalidate in C:/OSGeo4W64/lib/spatialite_i.lib - not found
CMake Error at cmake/FindSpatiaLite.cmake:75 (MESSAGE):
  Found SpatiaLite, but version is too old.  Requires at least version 4.2.0
Call Stack (most recent call first):
  CMakeLists.txt:251 (FIND_PACKAGE)

spatialite_i.lib ( I checked, rechecked and reinstalled couple of times ) is 4.3.0. If I remove spatialite check from cmake, I can create VS project, but then it gives errors on Link stage.