texts:
This question came from another topic: QGIS - creating VRT does nothing.
MacOS High Sierra (10.13), QGIS 3.10 Coruña, GRASS 7.8, GDAL/OGR 2.4.1, Python 2.7
The installation of GRASS was tricky (qgis config of several path to GRASS location)
I would like to create a VRT ("menu:raster:miscellaneous:Build virtual raster") or with toolbox:GDAL. When I open the window to select my files and hit "ok", both VRT window and select box close and it does nothing. When I try to merge several .asc file I get the same result.
When I try to create contour, it runs but log says that the result layer was not correctly generated.
Is the computer (8go Ram) just out of memory? (each asc file has 1 million points, it made a lot, but with just 4 tiles it does the same "nothing")
I guess a Python problem with a bad configuration of QGIS/Python. Got this "Python warning":

warning:/Applications/QGIS3.10.app/Contents/MacOS/../Resources/python/qgis/utils.py:737:
DeprecationWarning: the imp module is deprecated in favour of
importlib; see the module's documentation for alternative uses
mod = _builtin_import(name, globals, locals, fromlist, level)

I can't figure out which "utils.py" is failing and where it is located, inside qgis or python.
Do i have just to upgrade Python or Qgis? Or the plug-in which look to be the latest?


texts: Generating and populating point shapefiles from polyline shapefiles

texts:
Normally we have two shapefiles, Point and Polyline. We then populate my polyline layer with the point data ID where the lines start to where they end; eg. create two columns (from and to).

However in a GIS system representing electrical network, we have only one shapefile Polyline. No shapefile for Points. How can we generate and populate shapefile for points.

X and Y coordinates are available for start node and end node.

However in case of two line segments  of a polyline, the common point / joint should have should be represented by a single point only.
How can we generate and populate point shapefile and establish two columns in polyline with IDs of two points ?


texts: Getting well known text from varchar in redshift.

texts:
As you know, redshift does not support geometry types. There is a column in a dataset on my redshift db that has varchar column and the format looks like:

01060000000000000000___________________________

on and on for 765 characters. Since it looks like at some point it was a geometry type, I want to convert this into WKT with hopes of then converting to geom using 'ST_GeomFromText'. I feel like I may be complicating this more than it has to be, I've never worked with WKT so if you have suggestions on how to achieve the same result with a more simple approach I'd prefer that.

I've added a new column using:

SELECT AddGeometryColumn ('public','wkt_regions','geomm',4326,'Polygon',2);

and later:

ALTER TABLE wkt_regions ALTER COLUMN geom TYPE geometry(Polygon) USING st_astext(left(geometry, 764));

Without the left function to trim the geometry I get an error that says " Invalid hex string, length (765) has to be a multiple of two!"
But with it, I get an error that says "WKB structure does not match expected size!"


texts: Grass v.dist tool error when measuring distances from line segments to lines

texts:
I am using QGIS version 3.18.1-Zürich.
I am trying to measure the distance from various lines I created using the intersection tool, but I keep getting this error when using the v.distance tool:
WARNING: The combination of types is not supported by all formats.
WARNING: No attribute table found -> using only category numbers as attributes
Exporting 33990 features...
5..11..17..23..29..35..41..47..53..59..65..71..77..83..89..95..100
v.out.ogr complete. 33990 features (Unknown (any) type) written to  (GPKG format).
C:\Users\Matthew\Documents>v.out.ogr type="auto" input="vector_6070f4eda0a0f311" output="C:/Users/Matthew/AppData/Local/Temp/processing_WyqhMr/6573a748916e4ed884cb14f26484b549/from_output.gpkg" format="GPKG" --overwrite
ERROR 1: Wrong field type for fid
ERROR: Unable to create column

when I run v.clean on the pipe_incre line layer, as shown in picture below,  I get:
ERROR 1: Wrong field type for fid
ERROR: Unable to create column
C:\Users\Matthew\Documents>v.out.ogr type="auto" input="errora163fd5da5aa4b04b4590ec4cb90f3a6" output="C:\Users\Matthew\AppData\Local\Temp\processing_WyqhMr\ef8d0e593e064397973c1fb3065e7b93\error.gpkg" format="GPKG" --overwrite
WARNING: No attribute table found -> using only category numbers as attributes
Exporting 4918 features...
5..11..17..23..29..35..41..47..53..59..65..71..77..83..89..95..100
WARNING: 4918 features without category were skipped. Features without category are written only when -c flag is given.
WARNING: Output layer is empty, no features written


texts: How do you convert SpatialPoints object to PPP object?

texts:
I have a shapefile layer of polygons (houses in a city) that I want to run a kernel density on using the density() function in spatstat in R. To get my polygons as points I turned them into SpatialPoints using:
houses_points <- SpatialPoints(houses)
Now I have my houses polygons as points, but I now need to convert these SpatialPoints to a PPP object, so I can plug them into the density() function. How can I do this?
I have tried using:
houses_points <- as(houses, "ppp")
and get the following error message:
Error in as(houses, "ppp") : no method or default for coercing “SpatialPolygonsDataFrame” to “ppp”


texts: Locating DTEDs from one layer that intersect polygon in another layer using ArcPy?

texts:
I'm new to the GIS world, and using ArcGIS 10, so please forgive me if this is a simple question.  I have DTED files (level 2) that I have placed into a mosaic dataset, and I'm trying to determine which DTEDs intersect a polygon in another layer.  I'm using python, and used the SelectLayerByLocation method, but am having trouble getting the results from it.

I can't use the CopyFeatures method (which is what's in all the examples) because it's probably returning a mosaic layer instead of feature.  I also tried the Copy method, but that didn't work either.

Any suggestions?

Here's that part of my code:

   arcpy.MakeMosaicLayer_management('CreateMosaicDataset', 'mosaic_lyr')
   arcpy.SelectLayerByLocation_management("mosaic_lyr", "INTERSECT", buffershape)
   # If features matched criterias write them to a new feature class
   arcpy.Copy_management("mosaic_lyr", "DTEDs_for_Shape")

My solution was:

  # Create a Mosaic Layer from the mosaic dataset
    arcpy.MakeMosaicLayer_management('CreateMosaicDataset', 'mosaic_lyr')
    # Select DTEDs intersecting the polygon
    arcpy.SelectLayerByLocation_management("mosaic_lyr/Footprint", "INTERSECT", buffershape)
    # write selected features to a new feature class
    arcpy.CopyFeatures_management("mosaic_lyr/Footprint", "DTEDs_for_Shape")


texts:
I use the GeoTools API to read and render a compressed JPEG GeoTiff file:

 AbstractGridFormat format = new GeoTiffFormat();
    format.accepts(rasterFile);
    //this is a bit hacky but does make more geotiffs work
    Hints hints = new Hints(Hints.FORCE_LONGITUDE_FIRST_AXIS_ORDER, Boolean.TRUE);

   // InputStream stream = new BufferedInputStream(new FileInputStream(rasterFile), 40*1024);

    reader = format.getReader(rasterFile, hints);

    // Initially display the raster in greyscale using the
    // data from the first image band
    Style rasterStyle = createGreyscaleStyle(1);

    // Connect to the shapefile
    FileDataStore dataStore = FileDataStoreFinder.getDataStore(shpFile);
    SimpleFeatureSource shapefileSource = dataStore
            .getFeatureSource();

    // Create a basic style with yellow lines and no fill
    Style shpStyle = SLD.createPolygonStyle(Color.YELLOW, null, 0.0f);

    // Set up a MapContent with the two layers
    final MapContent map = new MapContent();
    map.setTitle("ImageLab");

    Layer rasterLayer = new GridReaderLayer(reader, rasterStyle);
    map.addLayer(rasterLayer);

The issue is that for bigger GeoTiff Files (e.g. file size 30 MB) the rendering takes a couple of seconds on an up to date machine, this is especially problematic if the user uses functions, like zooming.

Is there any way to use the API differently to speed up the process ? I did a rough profiling but only have seen that the renderer works a couple of seconds if the image gets refreshed.

The GeoTiff was created with gdal

gdal_translate -of GTiff -co COMPRESS=JPEG -co PHOTOMETRIC=YCBCR -co TILED=yes -a_srs EPSG:4326'


texts:
I've decided to give this map design thing a go - I downloaded the USA counties shapefile from National Atlas, exported them to SQLite from QGIS, and then linked the SQLite file in Tilemill, and this is the result I get:

I know that something with my SRS string is off, but when I selected WGS84, it still looks too wide. What should I be doing to make sure my map projection remains proper?

(The reason I exported to sqlite was that there are like 3000 nodes in that shapefile, and mapbox would only show 500 when I imported the shapefile directly)


texts: I am setting this font to be only 30% opaque but the display and printing of the font is not changing.  What am I missing?


texts:
I had a feature class made up of 1,700,000 polygons. I used Geopandas to create a geodataframe:

state = "MD"
state_gdb = r"C:\Projects\Pop_Alloc\{}_Data.gdb".format(state)
join_feat = "{}_Ftprnt_CB_Join".format(state)
bldg_feat_df = gpd.read_file(state_gdb, layer=join_feat)

No problem; it took maybe 5-10 minutes to run. I have another feature class; let's call it 'parcels.' It has around 2,600,000 polygon features. I tried to do the same thing; make a geodataframe.

parcel_gdb = r"C:\Projects\Pop_Alloc\Parcels_by_state.gdb"
state_parcel = "{}_Parcels_merge".format(state)
parcel_feat_df = gpd.read_file(parcel_gdb, layer=state_parcel)

It has now been running for several hours. Is there a reason for this? Do I simply not have enough in memory to create this geodataframe. Is there a way to resolve this issue (generator, chunking)?


texts: Flow data distribution within features of overlapping layers

texts:
I want to improve/adjust/comprehend the formula (see below) that was used in my previous question: Transferring flows (connections + values) between polygons for calculating flow values of mobile data respecting the spatial/geographic gist.

where
 - absolute value that is conveying
 - area of feature 'A' inside a cell 'I'
 - area of a cell 'I'
 - area of feature 'B' inside a cell 'II'
 - area of a cell 'II'

On the image below, there is a transfer of a value 100 from the cell 'I' to cell 'II'.

My question which factors can alternatively be used, besides the $area, in the formula for achieving better results in calculating the new flow value between feature 'A' inside a cell 'I' & feature 'B' inside a cell 'II'?

Moreover, I will be glad if someone can illuminate the non-meaningfulness/meaningfulness of this formula.

Aim: To compare new values from another source with values that are already known.

Besides the $area, I also own

addresses where inhabitants live
number of inhabitants per each address
thematic layers as polygons, e.g. Forest, Green areas, Streets, Buildings, Water etc.

From the above factors, I tried using the ratio of the area without forest, the number of inhabitants, and the number of addresses. I have access to QGIS 2/3, Python, R.

I previously tried asking it as Transferring values between polygons based on point criteria in QGIS?
, however, at this moment I suspect that it was too broad and I did not show enough of my own contribution.


texts: Cannot pass python string outputs as inline variables in modelbuilder

texts:
I'm trying to pass a couple of strings from python as parameters in the extract by mask tool in modelbuilder. Unfortunately it doesn't recognize them. When I type them in the parameters it thinks I'm trying to add files called %shp% and %tif% from and to my default gdb. This is my first ever attempt at passing values from python to modelbuilder so I have no idea how to go about trying to fix it. This is also my first question on StackExchange so hopefully everything is clear. Here's the code:

import arcpy
import os
invalue = arcpy.GetParameterAsText(0)

#Getting the basename of the input file
newbasename = os.path.basename(invalue).strip("'")
newname = os.path.splitext(newbasename)[0]

#Creating the strings for the output parameters
outshp = newname + ".shp"
outtif = newname + ".tif"

#Setting the strings as parameters
arcpy.SetParameterAsText(1, "C:/Base_Data/Topo/boundary/" + outshp)
arcpy.SetParameterAsText(2, "C:/Base_Data/Topo/geotif/" + outtif)

and the model:

As you may have guessed I want to batch process a bunch of these rasters where I have a shapefile of the same name as the mask and want the output to be the same name as the input but in a different location.


texts:
I have a bunch of points in a text file that are represented using lat long (116.3158516 39.9748016). I import them into my GIS enabled database and create trajectories using st_makeline() as geometries projected in EPSG:4326. I'm using geometry here because all the GPS data falls into an area that is 50km^2: correct me if this is my first mistake.

Now I'm doing things like calculating trajectory length (st_length) and simplifying the geometry using st_simplifypreservetopology() but according to the documentation the results for the length and the tolerance for the simplification are in the coordinate system of the data; which in my case is degrees. This obviously doesn't make sense. I need it to be in meters.

Do I need to change the way I'm storing my data or do I need to convert my trajectories and points into something else before feeding them into the st_length() and st_simplifypreservetopology() functions?

Screenshot for Mike T.


texts:
I downloaded the sample mbtiles data from  https://www.mapbox.com/help/data/trails.mbtiles

In the metadata of this file I found that the format is pbf.

But if I directly query the tile_data column in images table of it, and use some tools such as vector-tile-base to decode it, it will return errors with this message: google.protobuf.message.DecodeError: Wrong wire type in tag. (I have tested the tool on some mvt files and it works with no problem)

So how should I get the decoded pbf data from .mbtiles files?

Appendix:
Decoding code here:

import sqlite3
import vector_tile_base as vtb
tile_db = sqlite3.connect('trails.mbtiles')
cursor = tile_db.cursor()
one_example = cursor.execute('select tile_data from images limit 1').fetchall()[0][0]
vtb.VectorTile(one_example)


texts: mbtiles generated by tilestache does not open in qgis

texts: I have tilestache-seed.py working fine for a tms layer and also outputting mbtiles but the mbtiles generated does not open in qgis and according to mapbox the tiles can be separated into three tables but  the tile table does not contain tile_id  which field should I use as a tile_id and will it open in qgis.  I also want to export the tms into an mbtiles using mb-util so I need to know where do I put the metadata.json file in the tms layer.  any  guidelines are appreciated.  I have read all the documentation for tilestache, MB-Util


texts: Does an openlayers vector layer work across browsers?

texts:
I was under the impression that it was.

Here is my example:

<html lang='en'>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
<title>GeoRSS</title>
<script src="http://maps.google.com/maps/api/js?v=3.5&amp;sensor=false"></script>
<script src="http://openlayers.org/api/OpenLayers.js"></script>
<script src="google-v3.js"></script>

<script type="text/javascript">

var map;

function init() {
    var restrictedExtent = new OpenLayers.Bounds(-180, -90, 180, 90);

    map = new OpenLayers.Map("map_element"
    );
    map.addControl(new OpenLayers.Control.MousePosition());
    map.addControl(new OpenLayers.Control.LayerSwitcher());

    var layer = new OpenLayers.Layer.Google("Google", {"sphericalMercator": true});

    var usgs = new OpenLayers.Layer.Vector("Earthquakes", {
        projection: new OpenLayers.Projection("EPSG:4326"),
        strategies: [new OpenLayers.Strategy.Fixed()],
        protocol: new OpenLayers.Protocol.HTTP({
            url: "http://earthquake.usgs.gov/earthquakes/catalogs/7day-M2.5.xml",
            format: new OpenLayers.Format.GeoRSS()
        })
    });
    map.addLayers([layer, usgs]);
    map.setCenter(new OpenLayers.LonLat(-11682133.90525, 4676768.13795), 4);
}

</script>

</head>
<body onload="init()">
    <div id="map_element">
    </div>
</body>
</html>

For me this is only working on IE.  I would like to keep this xml file on the USGS server and call it dynamically rather than save it locally. Is there anyway to get this to work, on other browsers (besides IE) without using a proxy script?

Thanks!


texts: Visualize given temperatures on image/map/plane

texts:
I am looking for a way to visualize some data of the form :

latitude longitude temperature

as image(s). Have many of these 3-tuple lines uniformly distributed on a region.

Ultimately, I am looking for a dynamic way, but even a statically creation of an image would be fine.

I need the image in order to use it as an overlay for google maps (or something similar) and looking for an end result like the one here.


texts: Converting UTM coordinates to shapefile in QGIS

texts:
I have an excel sheet with UTM coordinates (Easting and Northing) for four points. I want to create a shapefile (rectangle) in QGIS, but am very new to the program and following all the guides I can't seem to get it to work.

I have converted the excel spreadsheet to a Tab Delimited Text file, and gone on to QGIS and gone on to "Add a delimited text layer". I have imported the .txt file, changed the file format to 'custom delimiters', where it is checked on 'tab'. But the issue comes when I try to set the x and y coordinates. I can see in the sample data box at the bottom that it has not imported my .txt file as a table at all (which is how it is appearing on all the tutorials). So as there no longer two columns with my x and y data, it will not let me set the correct titles to allocate the x and y data. I can see that certain tabs have also merged when converted to a .txt i.e. it is combining one of my coordinates with the tab title to the left of it.

I was then going to use the delimited text layer created, and then convert it to a shapefile from there. Would this be correct?

Sorry if this is a very simple question or is not very clear, I am very new to QGIS and mapping.


texts: How to prepare or normalize raw shape (.shp) files for OpenLayers2 use

texts:
It seems like I am having difficulties in normalizing raw map files' resolutions to OpenLayer's use most probably in parameters mapext and imgext. I hope you can help me with this.

Using raw (.shp, .shx, .xml, .prj, .dbf) and associated raster files for landuse-related data included in an ArcGIS 10 project (.mxd) that I've obtained from a credible private source, I have used Mxd2Qgs to read .mxd project to .qgs. Then I used Quantum GIS' MapServer Export plugin to export selected layers to a MAP file.

I can view the generated map in browser but cannot do so when I use it through OpenLayers. The generated map file has an extent:

EXTENT 272813.696194 1546212.773835 340292.781152 1587773.606956
# (it's quite huge).

I have tried adding to OpenLayer's map initialization properties maxResolution, numZoomLevels, units, projection and displayProjection but there is no map rendered, although the map has been queried.  I have just guessed the maxResolution here, copying from other huge datasets online :)

The map has never been queried when I added the map file's EXTENT to maxExtent: OpenLayers.Bounds(...). The same thing goes when I used the files included in the MS4W tutorial #1 data set.

My questions would be:

How to know a given map's maxResolution parameter for OpenLayers? Can this be guessed?
Are there any special processing needed for map/shapefiles before they can be used in OpenLayers?
How should a map file be written so that it could be used in OpenLayers?

EDIT: Attached MAP file and JavaScript files.


texts: Raster clipping without producing no data values

texts:
I'm trying to clip a raster in ArcMap using a polygon, and I'm wondering if there's any way to have the resulting raster not include any "no data" values (I apologize for the poor grammar). My process in pictures is as follows:

Input Raster:

Polygon as clipping extent:

Clip tool options:

Resulting raster with No Data shown in red:

I just want the resulting raster to consist of the cells that have values. Any hints or tips on how to do this?

P.S. I've tried the Extract by Mask tool and it doesn't help.


texts: Transforming geographical points to geometry in PostGIS?

texts:
I have a PointField(geography=True) field in Django, which is translated into the following Postgres type: geography(Point,4326).

I set the geography field to True because I thought it made sense for storing coordinates of physical locations in the world, but inserting points received from Google Maps causes them to be transformed inside the DB.

So for instance, I put the following point into Postgres: SRID=4326;POINT (30.2484381 -97.7503308):

UPDATE "table" SET "location" = ST_GeogFromText('SRID=4326;POINT (30.2484381 -97.7503308)') WHERE "uuid" = '1d604100-b838-4cd3-9ca3-f600648ed3d8';

Now, when I try to retrieve the point I just stored, it comes out transformed:

SELECT ST_AsEWKT("location") FROM "table" WHERE "uuid" = '1d604100-b838-4cd3-9ca3-f600648ed3d8'

This is the result: SRID=4326;POINT(30.2484381 -82.2496692).

I have lot of points stored that have been transformed into a geographical system, and now I want to transform them back into the geometrical points that were inserted originally.

Is there a PostGIS function to do this?


texts: QGIS vector layer labels are not showing from python script

texts:
I have a python module, where a QGIS memory vector layer is created, within the same method that creates this layer I have constructed labels using:

text_format = QgsTextFormat()
label = QgsPalLayerSettings()
label.fieldName = 'Description'
label.enabled = True
label.setFormat(text_format)
labeler = QgsVectorLayerSimpleLabeling(label)
layer.setLabelsEnabled(True)
layer.setLabeling(labeler)
layer.triggerRepaint()

And then subsequently from another Python module I add this layer to the QGIS workspace using:

QgsProject.instance().addMapLayer(layer)

The labels are not showing in QGIS, however, if I open the "Labels" dialogue from within QGIS and press "Apply" then the labels show.  So I assume that the label field and text is set properly.

I have repeated the problem described above in the Python console within QGIS v3.2.2 using the following script:

from qgis.core import *

layer=QgsVectorLayer('MultiLineString?crs=epsg:4326&field=ID:string&field=Description:string', "MultiLines", "memory")
pr=layer.dataProvider()
record = QgsFeature()
record.setAttributes(['1','test'])
record.setGeometry(QgsGeometry.fromWkt('MULTILINESTRING ((1 1, 6 1))'))
pr.addFeature(record)

text_format = QgsTextFormat()
label = QgsPalLayerSettings()
label.fieldName = 'Description'
label.enabled = True
label.setFormat(text_format)
labeler = QgsVectorLayerSimpleLabeling(label)
layer.setLabelsEnabled(True)
layer.setLabeling(labeler)
layer.triggerRepaint()

QgsProject.instance().addMapLayer(layer)

Any suggestions on what I need to add to my python script to have the labels show?

QGIS v3.2.2-Bonn was used.


texts:
This question is about staying organized and current on all your to dos and project tasks.  For me I am an avid user of GTD and a variety of online tools such as evernote, toodledo, and gmail.

I would love to hear about organizational tools that are used by my fellow GIS practitioners.  What tools do you use to stay on task and get your work done?


texts:
I have a data set which is an extract from OSM. I'm only using the roads shapefile and adding labels. The extract is the state of Pennsylvania. I then tried to import this layer to PostGIS from the DB manager within QGIS. I noticed the program was unresponsive so I figured it was working on the import and I left it alone all night to do its thing. The next morning QGIS and DB manager was still unresponsive so I forced the program to quit.

I then deleted parts of the dataset I did not need in attempt to make it smaller and easier to upload and import. Tried to import VIA DB manager again, added a primary key of my choosing, specified the SRID's and clicked ok. Once again DB manager and QGIS go unresponsive.

Is this normal for QGIS when working with OSM extracts? I understand the extracts can be rather large. Should I just let DB manager do its thing and hope it finishes soon?


texts: Calculating vegetation indexes with GDAL raster calculator?

texts: I want to calculate a few of vegetation indexes [ example: NDPI (red band-green band)/(red band + green band)] in a model with the GDAL raster calculator but I don't know how put the equation here:


texts: Is code completion available in QGIS Python Console?

texts: I wanted to know if the autofill functionality is available in the python console for QGIS as the one available in ArcGIS Desktop's Python Console or am I supposed to hardcode it.


texts: Statistically comparing surface slope of 2 different areas?

texts: I could use the statistical tests like t-test and ANOVA, however I am afraid those test will easily show significant differences if I use the slope values of all cells from each area.


texts: Adding legend to KML file uploaded to ArcGIS Online

texts:
I am using ArcGIS desktop 10.7 and I have a free ArcGIS Online account.
I converted a raster to KML file in order to be able to add the raster image to the ArcGIS online account. (I used the layer to KML tool). When I add the KML to the online map there appears to be no way to create a legend or somehow add the raster legend from my original map document. The link below describes a way to do this using the ArcGIS server but I am not sure exactly how to do this or if it is even possible with the basic account.
https://enterprise.arcgis.com/en/server/latest/publish-services/windows/authoring-maps-for-publishing-as-kml.htm


texts: How to perform volume calculation in QGIS

texts:
EDIT 07/12/12

Ok, meanwhile I came to another problem I want to discuss.
I just tried around with the volume calculation with some layers and came to a point where I receive the following error.
Screenshot:

I did some research and discovered this explanation:
http://blog.webforefront.com/archives/2011/02/python_ascii_co.html
The problem is, that I have absolutely no Idea about python programming.
So my question is: Does anyone know how to solve this described error problem?

Note: I left the older questions below in order to see the process I made.

I m currently trying to make a volume calculation in QGIS. I realised, that
using GRASS or Sextante is the right choice to do this. But since there are barely useful tutorials, that go into depht nor any other guides or hints that may help in this context, I want to seek help here.
So I want to ask these questions:

1)I want to calculate the volume of a specific map area. So I discovered the algorithm "r.volume."

1.1.)How does it work exactly?The Help text in that menue is not very useful for me as a Grass beginner.

1.2) Why do I need raster layers for this and no vector layers?

1.3) What certain criteria do these rasters have to fulfill?

1.4) I added a screenshot. What exactly are these "clumps" and what rasters have to be chosen in the drop down list?

1.5) Can I have an accurate example of how to perform such an r.volume algorithm?I d be thankful for any detailed help/ raster files to test this.

2)I learned, that GRASS (in QGIS) was integrated into the Sextante Plugin for QGIS.

2.1) Does this mean I can only use Sextante and still have all the functions of Grassgis?I d welcome that because I always had trouble with importing data into Grassgis. Sextante doesn t (obiously) require this.(?)

2.2) A strengh of gvSIG is Sextante, so I read; Is the Sextante Plugin for QGIS equal to gvSIG Sextante concerning performance/ possibilities?

thanks a lot,
Shepherdjo

EDIT:07/11/12 8.08am:

I know this is not easy to answer but it would be great too, if you can help me out with the other questions (2.1/2.2).


texts:
I can't figure out how to use the processing tool "merge divided roads" to summarize all feature information. I have a two-line road, sharing a road number. Both lines are separate and have their unique OBJECTID and also an attribute of interest to me - Flow volume. Flow volume tells me how much goods has travelled on the line of the road.

I would like to make these two lines into one, and I can do that geometrically with merge divided roads sucessfully:

But the new line only carries the flow attribute from one of the original lines. In the tool, an "out_table" can be generated, supposedly enabling me to generate the sum that i need. I just cant figure out how to use it.

It looks like this:

I can't figure out what to do with that. I guess I should use a join, but I've tried and I can't get it right.


texts:
I am trying to convert a shapefile to GML 3.1.1 . I am using gdal 1.9.1 however the GML conversion it generates for a geometry feature is GML 2 and it also swaps the long/latitude coordinates. Are there simple workarounds to keep the coordinate positions the same, while changing the output to GML 3.1.1?

Thanks.


texts: Landsat util error when downloading data

texts:
I am having trouble when trying to download data using Landsat util. Using the following command:
landsat download LC81710792016253LGN01

I get the following issue:
landsat download LC81710792016253LGN01
LC81710792016253LGN01 is not available on AWS S3 or Google Storage
Time spent : 1.69 seconds

I found the scene ID by searching therefore using Landsat-util so I know it is available. Some of the stuff I tried was to update USGS using:
pip install -U usgs

And to uninstall landsat-util, the update USGS and reinstalling without any luck. Some other forums suggest to use my credentials which might be the issue because when I use:
landsat -h

I get the following message:
--username          USGS Eros account Username (only works if the account has special
                                    inventory access). Username and password as a fallback if the image
                                    is not found on AWS S3 or Google Storage

--password          USGS Eros account Password

So I tried the command again with my credentials:
jaun@jaun-pc:~$ landsat --username jaun****@gmail.com --password Ja*****9_ download LC81710792016253LGN01
usage: landsat [-h] [--version] {search,download,process} ...
landsat: error: argument subs: invalid choice: 'jaun****@gmail.com' (choose from 'search', 'download', 'process')

It might be that I use the --username and --password command wrong but there is nowhere on the internet where I can get examples thereof.


texts: Using Python to backup ArcSDE geodatabase on SQL Server Express?

texts:
Is it possible to use Python to do a backup of ArcSDE databases stored in SQL Server Express?  I have been using a custom executable that uses Arcobjects

IDataServerManagerAdmin.BackupGeodatabase()

but, I would prefer to use Python. The executable makes our IT staff a bit nervous.


texts: Is there a list of possible value for SDO_GEOM.VALIDATE_GEOMETRY

texts: When the geometry is valid return 'TRUE' but when is invalid return an error code, You can google each one, but I want to have a system user-friendly where when I found a wrong geometry already said what is the problem so user proceed to fix it.


texts:
I'm a (reasonably) proficient GIS user (mostly QGIS, but some Arc GIS) who has been thrown the conundrum of converting/importing CAD urban planning drawings (2D drawings produce in AutoCAD 2014) to GIS shapefiles (or similar).

I have had some success in importing the drawings, but have realised that all of the drawings have been produced as lines (not polylines or polygons) which means I am losing a layer of detail that would be really very useful in GIS.

For example, I have a site plan of a new housing layout which shows hundreds of individual plots of land, which 'look' like polygons but are actually represented in AutoCAD as lines. As I understand it from reading through forum posts here it's probably best that these plots are converted to polygons in AutoCAD before we attempt to export/import them to GIS?

I'm really just looking for a solid enough workflow to produce 2D architectural plans in AutoCAD and import them to GIS so that I can now attach additional attributes to individual polygons (e.g. buildings, plots) and polylines (roads, rivers etc).

As anyone come across any useful resources that would point me in the right direction?!


texts: How to change PostGIS feature attributes via openlayers 3+?

texts:
I am trying to find a working example of editing the attributes of a feature in a PostGIS table published in GeoServer (2.8), using Javascript libraries like openlayers (3 & 4).

The closest thing I can find is the Modify Feature example  from the openlayers documentation, but it only touches geometry not attributes. In addition it is based on GeoJSON, not a database backend.

The next closest is an OpenGeo blog, but it's just a blog with a screenshot without and code snippet. Also relevant is a similar GIS.SE question: Map editor R & leaflet?. The answer there suggests that one has to program it from javascript.

Before reinventing some wheels, can anyone give me some pointers how to implement such a function for setting the attribute of a feature?

FYI, I have copied the said openlayers 4 example below as a starting point.

<!DOCTYPE html>

    <title>Modify Features</title>
    <link rel="stylesheet" href="https://openlayers.org/en/v4.0.1/css/ol.css" type="text/css">

-- EDIT --

This is not a duplicate of the said question How can I use WFS-T in openlayers?, as that question is about OpenLayers 2. As said in the beginning, my question is about OpenLayers 3&4, which is not backwards compatible with Openlayers 2, and missing a large number of features in OpenLayers2. So the said question/answers don't apply here. I've also used the openlayers tag per tag description to differentiate this from openlayers2 questions.


texts: Converting raster file from float to integer type in ArcGIS for Desktop?

texts: I have a float type raster file (an NDVI made from ERDAS IMAGINE) and would like to convert it to an interger type raster file in order to use Focal Statistics function (Majority). Is there any way to do this?


texts: Need help in ArcGIS Java Script API : Installation & configuration with different browsers

texts:
We are using ArcGIS Java Script API (version 3.2, 3.0, 2.8) according to our requirement we are using various API’s.

Currently ESRI launched new version 3.3 so if I want to view old examples references i.e version 2.8, 3.0 how do I view the old examples (Because they removed old examples)?
Similarly, Can I install multiple SDK( 2.8,3.0,3.2) on single machine? In addition, I need recommended browser settings as if one of my application is not displaying the legend in IE 8.0

Therefore, what are the recommended settings for all versions (2.8 to 3.2)

Any suggestions /discussions will be great !!!

Thanks in advance : )


texts: How to have a window popup as a toolbar on the right side (with fixed size) by default in QGIS?

texts:
I have created a plugin which first shows login window and then user needs to enter details and press submit. After details are verified at back-end, a new window appears with buttons for different tasks.

I wish to have this window appear as a toolbar on the right side in QGIS by default such that the user do not need to move this window again and again.

For example, below is AutoField Plugin Window on the right side of QGIS:

I have tried to do the same with below code:

dock_widget = QDockWidget(self.iface.mainWindow())
dock_widget.setWidget(MainWindow) dock_widget.setFixedWidth(356)
dock_widget.setFixedHeight(600)
self.iface.addDockWidget(QtCore.Qt.RightDockWidgetArea, dock_widget)

However, half of the bottom of this toolbar appears blank. How to have this toolbar fixed with particular size on the right hand side?


texts: Calculating distortion on Equirectangular Projection

texts:
I am trying to calculate distortion so I can distort overlaying text and forms to precisely match an image of an equirectangular projection.

How does one calculate the distortion at a given latitude on an equirectangular projection 1:45,000,000 (say, 2000 pixels wide x 1000 pixels high)?

I've been trying to figure out this post and its links to no avail:
How to create an accurate Tissot Indicatrix?

I am not a professional, just a very interested amateur.

Here's the long story.

I am visualizing/mapping data using the Processing programming language and would like to have the 2D mapped data (different sized fonts and circles) appear undistorted when wrapped to a 3D globe. The data is mapped using equirectangular x, y's and the maps I want to use as backdrops are all this projection, so I'm assuming I want to "match" this distortion (e.g. by calculating distortion via latitude using Tissot equations?). Using the programming language I can precisely distort both the text and the circles. I think all I need are the equations to do it correctly.

Here is the original 2D data map:

When wrapped it looks distorted, like this:

How can I make my 2D image look undistorted when wrapped to the 3D sphere?

For reference, here's the same question asked differently on the Processing forum.

I'm not sure I want to reproject to an orthographic projection.  I want my 2D data map to wrap to a 3D sphere model that can be interacted with (i.e. spun).

I am using a 3D modeling program (Cinema 4D) to wrap a sphere with a 2MB "Blue Marble" image (equirectangular projection) from NASA.

When wrapped it appears undistorted from all hemispheres (not just one hemisphere, as an orthographic projection would be?), see: still from 3D model above.  (The modeling program is doing the orthographic projection for me as I rotate the object, I suppose.)  Therefore, I think that if I distort my 2D data map in a similar way it too will appear undistorted on the 3D sphere.  Here's a shot I took with an equation that approximates equirectangular distortion.  You'll notice the egg shaped ellipses from the 2D image look like a circle when wrapped to the 3D sphere.  Similarly, the Tissot ellipses also appear as circles on the 3D sphere.

This is why I was looking at the Tissot equations...to more precisely figure out the distortion of the equirectangular projection at different latitudes so I could distort my overlay accordingly.

Perhaps I should use a GIS program.  I just downloaded Cartographica and will see if I can figure it out.

Any Mac software suggestions for someone new undertaking this task?


texts: Getting Latitude and Longitude from a Linestring

texts:
I have a linestring geometry in my Postgres database. I wanted to extract the latitude-longitude coordinates of each point inside the linestring. Example

 "0102000020E610000002000000CCDA5CEB41D323406E2A9CEE86C54A40CCDA5CEB41D323406E2A9CEE86C54A40"
"0102000020E610000002000000AF5FB01BB6CD22407C7E18213C8E4A40AF5FB01BB6CD22407C7E18213C8E4A40"

Till now, I have executed the query

SELECT st_astext(geometry) from warning_vmz;

result

"LINESTRING(7.17999 52.66463,7.17899 52.66563,7.17857 52.66604)"
"LINESTRING(7.2860666 53.1664885,7.2860666 53.1664885)"

so I have tried the another function

Select ST_DumpPoints(ST_GeomFromText(st_astext(geometry))) from warning_vmz;

output

"({1},01010000004C89247A19652240E8A4F78DAF8D4A40)"
"({2},010100000012BD8C62B965224076543541D48D4A40)"

If I use ST_DumpPoints(ST_GeomFromText(geometry)) then I get an error

ERROR:  parse error - invalid geometry
HINT:  "0102000020E6100000711300004C" <-- parse error at position 28 within geometry
********** Error **********


texts:
Through this question and answer I've learned about the nature of Multipolygon shapefile, which was defined as multiple exterior rings.
In short, (1) it contains multiple exterior rings (polygons), and (2)  its attribute table shows one record for these multiple features.

So my expectation was, by choosing Multipolygon in Create New Temporary Scratch Layer dialog I could produce such file, i.e. multi-geometry plus single record.

Choosing "Multipolygon" and create polygons with "Add feature",

then saving the scratch layer as shapefile, however, produced a polygon shapefile with multiple records. The picture below shows a purple shapefile named "test", which was created by "Save As" from the Scratch Layer. (RHS: Attribute table).

What did I do wrong? As multipolygon I anticipated a single record attribute table, but I see 4 records. Or, is this a bug?

(Note) Before saving, the Scratch layer (cyan "New scratch layer") also had 4 features in its attribute table, without attribute column.


texts: What routing system support off-road routing for Leaflet integration?

texts:
I have a collection of regionally specific linestrings and multi-line strings which are off-road paths. I would like to be able to integrate these into a Leaflet routing systems which would exclusively use these paths and no other existing roads/trails. From the research so far, it looks like I could use perliedman leaflet-routing-engine to hook into ORSM or Valhalla, but then I would need to feed the custom off-road paths into these systems. I would like to keep these paths private for at least a period of time.
Is building these out in OSM as private the best option for feeding the routing system?
Is there a better way to just feed these paths into ORSM or Valhalla directly?


texts: ArcGIS Javascript API Feature Layer retrieval

texts:
I would want to develop an application where we could query the layer based on a result set of certain attributes of the features and highlight features.

Approach
Currently I tried using the feature services of ArcGIS server to do so. The feature service has around 50k features.

Issue:
While trying retrieve and highlight the feature using a query task it takes quite a long time for the process to complete and the browser literally hangs.

Question:
1) Is there a faster way to retrieve and highlight the features?
2) Could I use a better approach rather than using a feature service to achieve my end goal? Please do suggest.

Have provided my code below.

  esri.config.defaults.io.proxyUrl = "proxy.ashx";
  dojo.require("dijit.layout.BorderContainer");
  dojo.require("dijit.layout.ContentPane");
  dojo.require("esri.map");
  dojo.require("esri.layers.FeatureLayer");
  dojo.require("esri.tasks.query");
  var map;
  var featureset=[];
  var infoTemplate = new esri.InfoTemplate();
  infoTemplate.setTitle("${ROADNAME}");
  infoTemplate.setContent( "ROAD NAME: ${ROADNAME}
"
                         + "CAT${CAT}"+ "LINK_ID${LINK_ID}");

  function init() {

 try {
         var initExtent = new esri.geometry.Extent({"xmin":103.55,"ymin":1.13,"xmax":104.16,"ymax":1.56,"spatialReference":{"wkid":4326}});

map = new esri.Map("map",extent:esri.geometry.geographicToWebMercator(initExtent)});

var basemap = new esri.layers.ArcGISTiledMapServiceLayer("http://server.arcgisonline.com/ArcGIS/rest/services/World_Topo_Map/MapServer");
var basemap_sing = new esri.layers.ArcGISTiledMapServiceLayer("http://karthikpc:8399/arcgis/rest/services/Carriage_Mercantor/MapServer",{ displayLevels:[0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19]});

var featureLayer = new esri.layers.FeatureLayer("http://karthikpc:8399/arcgis/rest/services/Carriage_Mercantor/FeatureServer/0", {
            mode: esri.layers.FeatureLayer.MODE_SNAPSHOT,
            infoTemplate: infoTemplate,
            outFields: ["*"]
        });

  var highlightSymbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_CROSS, new    dojo.Color([255,0,0]), 3), new dojo.Color([125,125,125,0.35]));
   map.addLayer(basemap);
   map.addLayer(basemap_sing);
   dojo.connect(map, "onLoad", new_function);
   }
catch(e){
    alert('An error has occurred: '+e.message);
    }

    }

 function new_function(){
    var queryTask = new   esri.tasks.QueryTask("http://karthikpc:8399/arcgis/rest/services/Carriage_Mercantor/FeatureServer/0");
    var query = new esri.tasks.Query();
    query.returnGeometry = true;
    query.outFields = ["LINK_ID","ROADNAME","CAT"];
    query.where = "CAT='CATA'";
    dojo.connect(queryTask, "onComplete", function(featureSet) {
          map.graphics.clear();
          var symbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_CROSS, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new dojo.Color([255,255,255,0.35]), 1),new dojo.Color([125,125,125,0.35]));
   dojo.forEach(featureSet.features,function(feature){
   var graphic = feature;
   graphic.setSymbol(symbol);
   graphic.setInfoTemplate(infoTemplate);
   map.graphics.add(graphic);

   });

   });
   queryTask.execute(query);
   }

   dojo.addOnLoad(init);


texts:
I'm having difficult using basic date functions in cartodb.

This could be my formatting or syntax, but I'm stuck.

I simply want to show customers on my map that have not been visited in the last 90 days.  I have a table of geocoded customer points and a column 'last_visit'.

The following sql statements fail:

SELECT  datediff('d', last_visit, currdate()) FROM table_2016_10_26_geocoded_pids_y

'fucntion currdate() does not exist'

SELECT  datediff('d', last_visit, curdate()) FROM table_2016_10_26_geocoded_pids_y

'fucntion curdate() does not exist'

SELECT  datediff('d', last_visit, today()) FROM table_2016_10_26_geocoded_pids_y

'function today() does not exist'

SELECT  datediff(d, last_visit, curdate()) FROM table_2016_10_26_geocoded_pids_y

'unknown column d'

SELECT  datediff('dd', last_visit, now()) FROM table_2016_10_26_geocoded_pids_y

function datediff(unknown, timestamp with time zone, timestamp with time zone) does not exist

'column d does not exist'


texts: R colored raster from QGIS to ggplot or tmap

texts:
I am fairly new in doing GIS work in R, but I thought this would be a nice situation to get my feet wet.

Here is the situation: I need to create multiple plots with the same background, the same polygon layer from a shapefile, but a different selection of points (another shapefile), each time an other animal species. So I would think to create a loop where I run through a list of species, create the plot (either ggplot of tmap, or something else) and save the plot as a PNG.

I loaded the shapefiles with st_read() (from the sf package), and making the selection and put them in the plot is fairly straight forward. But I have a problem when I want to add the background map.

A few things I have tried:
Download OSM tiles using the read_osm() function from tmaptools. But my shapes are in a local Projection (Amersfoort/RD New) and the map is slightly misaligned with my shapefile for some reason.

Download tiles with get_map() from the ggmap package, but this fails with a '400 Bad Request'

Export a TIFF of the background map I want in QGIS is and try to import it into R using the raster function. This solves the alignment, but the layer is not in colour.

I think the last route might be the simplest, but than my question is: how do I export a map in QGIS that I can import in R to display as a coloured background map in ggplot, ggmap or tmap.

But also other suggestions, I am very new to this, and there is a lot of different information out there so I might be doing this all wrong...


texts: Creating catchments/basins around stream feature class using ArcGIS Desktop?

texts:
I have a reasonably accurate steam network for a watershed.  This network includes the main channel, major forks, and tributaries.  Each channel, fork, and tributary are individual records with single part geometries.  The lines are topologically connected such that each tributary is snapped to its downstream geometry.  Each line's geometry equals the flow direction of the watercourse it represents (lines point downhill).  I also have a surface model with the same extent as the overall watershed.  While these data currently exist in a vector datatype I can transform them into a raster datatype if needed.

I would like to create a new  dataset such that each division delineates the (approximate) drainage for each tributary, fork, or main channel.

@WhiteboxDev states that  Several [QGIS processing tools ] take stream network inputs and map subbasins in Creating catchment from point using DTM raster in QGIS? but I am unaware of what tool those would be.

Using ArcGIS Desktop, how do I create areas around my existing stream features where each area approximates the drainage around each stream feature?


texts:
I'm revisiting GIS after several years of not using it and have opted for QGIS (currently v2.14). All the example files I have are of the US (Alaska) but I am learning QGIS for work purposes. I live in the UK so want to find a shapefile of the UK mainland. It doesn't have to be anything fancy just an outline would do, and also ideally in the British National Grid CRS (although I understand you can reproject but I haven't figured out how yet!)

Does anyone know where I can get shapefiles of the UK for use in QGIS?


texts: How do I use UK ordnance survey maps with OpenLayers API without committing myself to using the OpenSpace API. I need to have generic OpenLayers code only.


texts: Leaflet maxBounds doesn't prevent zooming out

texts: In Leaflet 1.0.2, if you set maxBounds on a map, it prevents panning beyond those bounds, but does not prevent zoom out beyond them. Is this by design? It seems incorrect behaviour to me. I need to prevent display of the map beyond given bounds, and thought maxBounds would do it.


texts: Selecting numbers within the slash (/) using ArcMap Field Calculator?

texts: I have date data that I'm separating into three separate fields (all fields are string type, including the input field) by day, month, and year. Unfortunately, the date field is raw data and isn't consistent; there are eight variations for how the data in this field is organized (d/m/yy; d/m/yyyy; d/mm/yy; d/mm/yyyy; dd/m/yy; dd/m/yyyy; dd/mm/yy; dd/mm/yyyy). I figured out how to separate the day and the year into separate columns, but I'm having difficulty separating the month. I've tried the VB Script Mid() function by initially separating the different formats and performing the function, but that was extremely time-consuming, and I have close to 100,000 data points that I need to go through. Is there a python code that I can use, or a combination of different VB scripts, that can select the numbers between the slashes (/) for all of these? Or am I doomed?


texts: How to make graphs of latitude and longitude and generating summary table

texts: I am not from geography. I am from quantitative finance. I want a graph with longitude and latitude of various locations of universities of my academic literature. And want to generate a summary table. I am attaching both map and table. I have tried Leaflet library in R and Folium in Python. But not able to generate this kind of map.


texts: Choosing coordinate system to base coordinates off of Meridian of Mexico?

texts:
I am using a source that is very old and all of the latitude coordinates are based off the Meridian of Mexico.

How would I search for these on my map?

I am working in QGIS 3.0. Also here is a map using the same coordinate system:


texts: Listing grid cell numbers for outlet according to flow direction using ArcGIS Desktop?

texts:
I am trying to list grid cell numbers whose flows contribute to the red cell I draw in the figure according to the known flow directions. Directions and colors are shown in the figure.

What tool can I use to list these grid cell numbers?


texts: How do I resolve error DataGridView/Label1 is not declared?

texts:
' set the active display
        Dim disp As ESRI.ArcGISExplorer.Mapping.MapDisplay = ESRI.ArcGISExplorer.Application.Application.ActiveMapDisplay

    ' define the item that selected by user
    Dim SelItem As ESRI.ArcGISExplorer.Application.SelectedItemsCollection = ESRI.ArcGISExplorer.Application.Application.SelectedItems

    ' define the name of selected item
    Dim SIName As String

    ' define a featurelayer for selected item
    Dim FLitem As FeatureLayer

    ' get the name of user selected item
    SIName = SelItem.Item(0).ToString()

    'display the name on the label
    Label1.Text = SIName

    ' find the featurelayer that has the name which equals to the one user selected
    FLitem = disp.Map.FindByName(SIName)

    ' create a new bindtable and fill it with featurelayer's table
    Dim FLbind As TableBindingAdapter
    FLbind = New TableBindingAdapter(FLitem.Table)
    FLbind.Fill()

    ' using the datagridview to show the bindtable which is the selected rows
    DataGridView1.DataSource = FLbind


texts: How can I stream attitude (roll,pitch,yaw) into Google Earth?

texts:
I am streaming position data (lat, lon) of an airplane into Google Earth via NMEA 0183 ($GPRMC, $GPGGA sentences) and it works OK (I generate and stream these from my own software). Now I would also like to stream attitude (roll, pitch, yaw) but I can't find any NMEA sentence for this, except proprietary vendor sentences which Google Earth seem to ignore.

In Google Earth, I selected Tools -> GPS -> Real Time -> NMEA, but I haven't tried "Garmin PVT" yet as I can't find any reference to the protocol in connection with roll, pitch and yaw (found this document though).

Is there any way? I know .kml files can include this information, but I can't stream a .kml file into GE (only one-time import).

This is what I see now. Notice no pitch, roll or yaw data:

Roll, pitch and yaw:


texts: Exporting geographic coordinates after Project in ArcGIS Desktop?

texts:
I used the 'project' tool to project all layers in my map (ArcGIS v10.2) to the same coordinate system (UTM GDA 94 zone 56) to do a spatial join. However, the coordinates remain the same in the attribute table (I understand that this is normal). However, I need to export all the coordinates in UTM GDA 94 zone 56 to an excel file. I had thought I could use the tool 'Convert coordinate notation' to convert the points and then do an excel export to accomplish this, but I have not had success with the 'Convert coordinate notation' tool.

Is anyone able to assist?


texts:
I have generated tiles using gdal_retile command. The tiles do not have wld file.
Is there any way to generate a wld for each tile (bash processing is preferred since I have around 9K tiles stored).


texts:
I can import arcpy and have been using ArcGIS Desktop, but in the tutorial I am following

mxd = mapping.MapDocument("CURRENT")
NameError: name 'mapping' is not defined

This is frustrating because I can access arcpy.mp and other sub modules.


texts: Creating a 3D buffer polygons (plane) from a 3D polyline that follow line catenary curvatures

texts: I am looking to create a 3D buffer polygon that represents a plane from a 3D polyline feature that represents an electric power (transmission) line. This polyline has Z-values that give it catenary curvatures when viewed in ArcScene. I am currently using ArcGIS and would prefer an ArcGIS method to add my current workflow. This plane is intended to differentiate vegetation points that are taller than the electric power line that could potentially fall in a strike the line wires.


texts: Calculating sum of parts of column based on another column using QGIS Field Calculator?

texts:
I would like to calculate the sum of values in one column, based on values in another column. To illustrate my question I uploaded the picture below. The second column consists of zipcodes, the third total sales. I would like to know what the total sales per zipcode are? So in this case:

What is the sum of column three for column two's value 1023? Output then should be 4 (0+1+1+2+0). The sum for 1024 then should be 11,5 following the same logic.

As this is quite a large dataset I would like to have the QGIS field calculator (or any other tool) calculate a new field which states the total sales in the zipcode. So far I can't seem to get it working and was hoping someone here could help me.

As you can see I have some missing values, is this a problem in calculation? Perhaps good to mention is that it is a table.


texts: How do I make two layers of point data align perfectly?

texts:
I have two point data sets that represent the same business locations, however, the points are slightly off for one of the two layers. As in one layers points are in the correct position over the building, while the other layers points are slightly off (in the parking lot, neighboring road, etc.).

I need to match the dataset that is slightly off to align perfectly with the dataset holding correct point locations.


texts: How to eliminate polygon overlap areas in QGIS?

texts:
Is there a way to remove polygon overlaps within the same layer using QGIS, preferably removing the overlap area of the biggest polygon (having the greatest surface)?

All of my polygons are in the same layer, and I would like to perform this operation on the entire layer at once rather than manual editing.


texts: Interpolate a Smooth Contour From Points

texts:
I have points scattered around some origin point, each with a score. The score should be zero at the origin (similar to distance). I want to generate smooth contour lines from these points.

The red numbers are the score for a particular point, the blue icon is the origin.

I tried the Contour tool (from vector) and the result is far too jagged. Next I tried IDW Interpolation with a relaxed distance coefficient (p=2, max is 100). The result is closer to what I want, but still a bit overfit.

I'm considering removing the far-away point to the south and expanding the extent with a template polygon. Any other recommendations for making this look better?

EDIT:
Wanted to mention I tried a lot of other methods too. Heatmap KDE offers a way to scale the radius by my score in addition to the heatmap weights. I'm curious how this works. Is the radius used for the gaussian kernel selectively larger for those points with a larger score? How is this different than for the weight?


texts: Adding new field and updating field value from another shapefile?

texts:
I have a shapefile of route (ObjectID,Name,FirstStopID,LastStopID,length) and another shapefile of stops(ObjecID,Name,RouteName..).

I want to have the name of stops in the route shapefile, so could anyone tell me how I can add two new fields and copy the information from the stop shapefile?

All the route and stops are shapefiles exported from the Network Analyst extension.


texts: Creating maps programmatically using ArcObjects. Adding legends, scale bars, north arrows?

texts:
I'm looking for some examples related to the publishing of maps programmatically to image format using ArcObjects. More specifically I would like to add standard map elements such as a scale bar, legend, and north arrow.

Currently I am able to create an MXD, add several layers from paths stored in a dictionary, set the symbology based on layer files, set the extent, and export the map as a JPEG. Up to this point I haven't found any examples of ArcObjects code that also makes the map look more professional by adding additional elements (and sizing/positioning them).

Any help would be much appreciated.


texts: How can I keep the original name of the .kml file in the layer panel?

texts:
I have a folder with many .kml files.
Example:

John.kml
Jessica.kml
Donald.kml

When I add them to QGIS, the names shown in the layer panel changes
Example:

John.123.kml
Jessica jess.kml
Donald_ert.kml

Is it possible to keep the original names when I add the files?


texts:
I am trying to find all cells that are a certain distance and within a certain window from a cell of interest.

For example: I want all cells in a distance of 10m within a search window of 10 - 40 degrees (from North).

Any hint how to do that using either a GRASS function or shell?


texts: Measuring distance from a point layer to a raster layer

texts:
I would like to calculate the distance (in meters) away camera traps (point layer) are to a water source.

Now the water source data I have been given via an external source from the national park I am collecting the data from however I am unsure on how to determine the distance a camera trap is to the nearest water source.

Is this possible to do between a raster layer and a vector point layer?


texts: Network Analyst service area Tobreak value changed

texts:
MY problem is that when I was calculating the service area, the Tobreak value had changed.
I selected 500 meter as my default Breaks.

But the result shows me that some "Tobreak" value is smaller than 500

How can I solve this problem?


texts: How to get hold of the Activated Event for the ArcMap window?

texts:
I am using ArcObjects 9.3 .NET. In my ArcMap dockable window, I have a single element host hosting a WPF user control (That's the only means of using WPF for ArcMap extensions as far as I know).

While most things work fine, I am experiencing a known issue with WPF/Winforms interoperability that I can't quite solve. In particular, when I have focus on my dockable window (a textbox in my WPF control in the element host to be exact) and then alt-tab from ArcMap and go back into ArcMap again, the cursor continues to blink in the textbox within my dockable window, but no text appears when I type. In fact, the focus seemed as if it went to the main ArcMap window, and anything that I type appears to trigger the ArcMap menu shortcuts (F -> File, E-> Edit) etc.

Further research led me to this thread, with a similar problem:
http://social.msdn.microsoft.com/Forums/en-US/wpf/thread/1b50fec6-6596-4c0a-9191-32cd059f18f7/

The workaround proposed was to listen to the Activated Event of the parent form. Upon the Activated event, do a .Focus() on the WPF user control.

However as the parent form in this case is probably the ArcMap window itself, I am having trouble accessing it from the DockableWindow. I can get the dockable window's parent control via .Parent or .ContainerControl, but these controls don't expose the .Activated event as they are technically not Form objects.

Any ideas how I can listen to the activated event (aka when ArcMap gets refocused from an alt-tab event)?

Thanks in advance!


texts: Masking an image in Erdas using Thematic Image, the mask comes out looking wrong and completely different than input image

texts: I am trying to mask out shadows and clouds in a WV2 8-band image. I classified an image to identify the shadows and clouds and then using the mask took in erdas used the classified image to mask out the specified pixels in the input image using "setup recode". The first time I did this the image came out looking great. The second time it came out looking hazy, the color map was off, and just all around did not look good. I did the exact same steps as I had the first time and now I can't get the mask to come out right. Any ideas what I am doing wrong? I just don't understand because I did the exact same thing as I did when it came out right and yet I am getting a different output.


texts: Plugin development for QGIS 3.16 (Windows): how to handle dependencies on external packages?

texts:
I've been working on a plugin for use with QGIS 3.16. I have developed it locally and I want to now deploy it so others can use it. The problem or rather question I'm having is: How do I deal with dependencies on external packages?
I found this answer:
Development of a plugin which depends on an external Python library
and this discussion:
http://osgeo-org.1560.x6.nabble.com/A-pipinstall-plugin-is-possible-First-What-s-the-difference-between-the-Osgeo4w-Shell-td5107633.html
But I'm wondering how people handle this right now.
I'm pretty new to Python and QGIS both so maybe my concerns just stem from not being informed. Some of the packages I had to install are a little tricky to get working on Windows (geopandas for example), which is the first part I'm not sure how to handle. I've read that installing with conda seems to work reliably, but that doesn't seem to be a viable option since I have no idea whether a user will have conda installed (or am I missing something here?).
I'm basically looking for some documentation/recommendations to handle dependencies semi-gracefully. I would like to avoid changing the users setup if possible, since my plugin might need a different version of some package and another plugin might require a different version. Any tips/good reads on how to go about this?
Is it possible to have the plugin in its own virtual environment for example? Or can I somehow download the necessary packages and include them from a file? For geopandas, that would probably mean also downloading its dependencies, can I then force geopandas to use those versions of packages instead of what might be on the users path? Would you recommend such an approach?
I'm working on Windows by the way, and the plugin does not have to be cross-platform at least for now.
I realize this is a pretty broad question, but I'm hoping to get some tips on where to start looking for a solution.


texts: How to fix "Connection failed - Check settings and try again" when connecting to PostGIS?

texts:
I have Postgres running via Postgres.app on rails already and am trying to connect QGis to the database.

PostGIS is working already. (questions already found online deal with this issue only)

I do not have a password or username set up for this database, as indicated in my database.yml file for rails:

development:
  adapter: postgis
  encoding: unicode
  postgis_extension: true
  schema_search_path: public,postgis
ffsdfsdfsdf  pool: 5
  database: gisdatabase
  host: localhost

So I've gone into QGis and am trying to set up a new PostGIS connection:

Service:
Host: localhost
Port: 5432
Database: gisdatabase
SSLMode: allow
username:
password:

I've changed the port around and managed to get this error:

"could not connect to server: Connection refused..."

Changing the port back to 5432, I get a different error, thus indicating that the connection is not being refused:

"Connection failed - Check settings and try again."

I'm thinking that perhaps its the blank "Service" form or perhaps an error with the empty username and password fields that are causing the error, but I don't really know.

I've tried to find the correct Service name but have been unsuccessful


texts: Separating two lists within single object using ui.chart.array.values

texts:
I have the code below, however, when I try to plot my chart, I am getting one error message. I believe this is related to the fact that within my object VIIR2014List I have two lists, one for average radiance, another for coordinates.
How can I get rid of this error message?
////////////////------------------------------------------------ VIIRS DATA
var VIIRS = ee.ImageCollection("NOAA/VIIRS/DNB/MONTHLY_V1/VCMSLCFG");
VIIRS = VIIRS.select('avg_rad');

var latLon = ee.Image.pixelLonLat().select('longitude');

var VIIRS2014 = VIIRS.filterDate('2014-1-1','2015-1-1').mean().addBands(latLon);
var VIIRS2019 = VIIRS.filterDate('2019-1-1','2020-1-1').mean().addBands(latLon);

var VIIR2014List = VIIRS2014.reduceRegion({
  reducer: ee.Reducer.toList(),
  geometry: mydiagonal,
  scale: 1000,
  maxPixels: 1e9
});

print (VIIR2014List, "Lists 2014"); /// 2014

var VIIR2019List = VIIRS2019.reduceRegion({
  reducer: ee.Reducer.toList(),
  geometry: mydiagonal,
  scale: 1000,
  maxPixels: 1e9
});

print (VIIR2019List, "Lists 2019"); /// 2019


texts: substitute for landsat 7 for land classification?

texts:
I am doing a project where I need to do a land classification for a forest reserve in Honduras. The dates I need are 2005, 2010, and 2015. I'm good for 2015 with LS8, but I am having difficulty with imagery for the other dates, especially 2005 because of the SLC-Off issue with LS7. I have done some mosaicing with some of the LS7 images to destripe it but it has still left some bad pixels. Because this area is outside of North America I cannot find LS5 data for 2005.I also face issues with cloud cover so finding enough quality images for a mosaic is a challenge.

Any suggestions or ideas on other satellites I could use or how to deal with the LS7 data? I have access to Arcmap, ENVI and geomatica.


texts: Can I create, union and intersect 3D multipatch features using 3D Analyst extension (ArcGIS 10)?

texts:
I'm on ArcGIS 10. According to all the documentation I've read, 3D Analyst should be able to do what I am trying to get it to do, but I am running into problem after problem and I'm almost ready to give up and find a simpler way to answer the question I'm being asked.

In a nutshell, I have several layers I need to intersect or union in 3D. Some are points (thousands of them, which each represent a 3D volume, either a cube or a sphere, the dimensions of which are attribute values), some are polygons. I'm pretty sure the 3D Union and 3D Intersect tools only work on Multipatch features so I need to convert all of them to Multipatch features first.

The main problem I'm having so far is that the 3D Analyst Union and Intersect tools refuse to operate on Multipatch features that are not "closed", and I can't seem to force ArcGIS to create Multipatches that are closed. I've tried converting my one polygon layer (8 features) to a Multipatch both by extruding it in Arcscene using its height and converting that to multipatch and also using the "Extrude Between" tool, and they both produce some multipatches that are closed and some that aren't. I tried converting my points to spheres, using an example I found on ESRI's site here:

http://www.esri.com/news/arcuser/1010/3danalysis.html

Which tells you to symbolize the points as spheres, using an attribute as the radius, then use the "Layer 3D To Feature Class" tool to make them into multipatches. That also produces multipatches that are not closed and therefore useless. I haven't tackled converting points to 3D cubes yet.

So if anyone knows a way to reliably create closed multipatches from 2D polygon layers and/or center points (with attributes describing the dimensions), let me know. Any alternative software solutions would be great too, I'm not tied to ESRI. Indeed, after mucking about with this for the last few days and realizing that ESRI's own tools create multipatches that can't be used with many 3D Analyst tools, an alternative would be preferable.

ETA: What I'm ultimately after is to calculate the volume of space that is occupied by (for example) cubes with a an attribute value of X where they intersect spheres with a value of Y (of course, I have lots and lots of values of X and Y to run).


texts: Finding blank attributes and update using ArcPy?

texts:
I'm trying to find and update all empty attributes for all fields in shapefile with "-" by iterate through all fields in shapefile to find which fields are empty.

I tried tp update all empty attributes using the UpdateCursor and it is not done

import arcpy
from arcpy import env
import os

# Set the workspace for the ListFeatureClass function
env.workspace = "c:/base"
fcList = arcpy.ListFeatureClasses()
for fc in fcList:
     fields = arcpy.ListFields(fc)
         for field in fields:
            with arcpy.da.UpdateCursor(fc2, field) as cursor:
                for row in cursor:
                     **#here I donot know how to check blank values in all feileds#**

I am struck here and from here I don't know how to iterate all rows in all fields and update all empty rows with "-".

I am attaching example


texts: Creating layer that intersects two attributes using ArcGIS Pro

texts:
I have one polygon feature class.
I am trying within this feature class to show clashes of polygons by specific dates.
Does anyone know how the best method to do this within ArcGIS Pro and then ensure it is picked up in real time?
I was thinking to do this with FME and it is straight forward with FME but I am wondering if Pro has another way of doing this.
Is there a way to do this with ArcGIS online also using arcade expressions?


texts: Projected raster data has inconsistent extent but not shapefile

texts:
I want to calculate the average value of a raster file within each polygon boundary of my shapefile. Because I want to preserve area I wish to use the Goode Homolosine (Land) projection. The original shapefile has a geographic coordinate system WGS1984, and the original raster file is also in WGS1984.

I used the Project tool to project the shapefile to the Goode Homolosine (Land) projection, and the Project Raster tool for the raster file, also choosing Goode Homolosine (Land) as the projected coordinate system.

It is my understanding that because both of the original files have a defined coordinate system, i.e., WGS1984, I should not be using the Define Projection tool.

Problem: When I open the newly projected raster file in ArcMap, I receive the following error message:

  Warning, inconsistent extent!

  One or more of the added layers has an extent that is not consistent with the associated spatial reference information.

  Re-projecting that data in such a layer may lead to unexpected behaviour

Two things to note:

I do not receive this error message when I add the newly projected shapefile;
This error message occurs even when I add the newly projected raster file to an empty data frame (which, if I understand correctly, should automatically determine the coordinate system for me).

I have just spent a considerable amount of time reading on how to re-project data and cannot see any problem in what I am doing.

~~~~~~

Raster data: http://www.earth.columbia.edu/people/gmccord/sitefiles/file/malaria_ecology.zip
Shapefile is not publicly available, but metadata states GCS_WGS_1984


texts: Which bands or channels these NOAA AVHRR satellite images involve?

texts:
These two satellite images are from NOAA AVHRR, but I do not know which bands or channels exactly both images represent.

Can you tell me the band or channels that are involved?


texts: Visualizing raster coverage of PostGIS in OpenJUMP?

texts:
I am new to OpenJUMP and I am trying to use this software for visualizing PostGIS rasters.

I was able to display vector layers of PostGIS in OpenJUMP but I cannot display raster tables. Do you know how I can do it?


texts: Extra vertices being added to polyline during ArcMap editing of shapefile?

texts:
I have a layer of multiple polylines each of which had some vertices. While editing vertices in ArcMap 10.5.1, all lines have become features with too many vertices. Does anyone know what is the reason and how can I change them back

These are the steps and details on what I made:

Create feature; Draw polylines
Start "Editor": Edit Vertices, move some of vertices and change their curves
Double click to finish the edits
After some edits, the I encountered too many nodes on each polyline instead of the nodes I already draw.

The data is a shapefile (polyline)


texts: Why are edits in shapefile table not saved?

texts:
I have created a shapefile (point data) and wish to change some parameters in the files table data.  I can edit the data and see the coresponding change in my my symbology.  However, when I save the shapefile, the data and symbology revert back to the way it was before the edits.  What am I doing wrong?

I have a shape file of a series of points.  One column in my table has values between 1 and 6 that show up as color-coded points on my map.  The data for the column came from a spreadsheet that I had joined with the shapefile.  I was trying to update certain values in the column (I set the shapefile so I could amend it, opened the table and changed the value.  The color of the point on my map changed, however, when I tried to save the changes, the color (and value) reverted to the previous settings.  I did find a way to make it work.  I amended my original spread sheet in Excel and exported it as a dbf.  I then unjoined the dbf and shapefile and rejoined with the new data. Don't know if that's the correct way, but it worked.


texts: joining a table from zonal statistics as table with a layer, problem with reserved names

texts:
I have a raster file (percentage of total forest killed) and a shapefile with polygons (municipalities). I thought as a first step I would like to find out which municipality shows the greatest percentage of forest kills (along with other stats).

So I ran a Zonal Statistics as Table, with the idea that I can join the table after with the polygons, so that I can color them accordingly.

I have done that several times before but I keep having to do the same annoying steps and I wonder if there is a smarter way around it.

All the stats names form the Zonal Statistics as Table output are reserved words and cause the join to fail, so I always have to copy and rename all the fields (min, max, count, sum ..) and delete the old ones. Which is not too bad but I was wondering if someone has figured out a way around this?

I am using ArcMap10


texts: How to zoom to layer which is invisible at load time?

texts:
I have a layer whose maxResolution is set when it is added to the map, so that it will only be visible at a certain zoom level.

When the user clicks on a button called zoomto layer, I want to zoom to that particular layer. But when I try to take the extent of that layer, it's always getting infinity as it's not visible in the map.

My code:

let lay = this.mapService.getLayerByName(layer.layerName);
    let resolution = lay.getMaxResolution();
    let map = this.mapService.getMap();
    let zoom = Math.floor(map.getView().getZoomForResolution(resolution)+1);
    map.getView().setZoom(zoom);

    let layerExtent = lay.getSource().getExtent();
    console.log("layerExtent",layerExtent)

    this.mapService.fitExtent(layerExtent);
    map.getView().fit(lay.getSource().getExtent(), map.getSize());
    map.getView().setZoom(zoom);

Is there any way to zoom to a layer which is not visible at that moment?


texts:
Imagine a table of points of interest along with lat/long. The user identifies two of these points. I want to do the following:

Calculate a great-circle path between these points
Find all POIs that lie within "n" miles of the great-circle path from step 1.

I'm totally new to this stuff, what might be a good tool for accomplishing this? And how would you go about solving it using that tool?


texts: QGIS Merged Layers/ Cannot Find Windows "minidump" to Restore Past version off Dropbox

texts:
I've copied a layer style and was playing around with the display of a point layer when QGIS Valmiera 2.2 crashed. A recently-created polygon layers attributes are now appearing in the place of another polygon layer with 100x more features, although the displayed feature count is still 506 instead of the 4 in the table. As the program crashed, it flashed a message as to where the minidump was located- not sure how to find it.

I attempted to restore previous versions of the layer and all associated files from Dropbox, but received an error of "permission denied". I verified that I have full access to AppData\Temp, which makes me believe that Temp is not the folder I should be searching in. Please help me find this folder or reload the 506 features and hours of work!


texts: How to create an interactive based filter (checkboxes) with multiple layers using Carto VL

texts: I'm relatively new to coding, and have been developing a map using Carto VL and Mapbox. I've found a way to create a filter within a map where users can turn off and on different categories within that filter. However, I'm trying to create a layer selector that allows for users to turn off and on different layers using checkboxes, or something similar. I've been unable to figure out how to add just one checkbox that responds to an entire layer, as opposed to specific types of data within the layer. How would I be able to create this multi-layered layer selector?


texts: Get addresses from TIGER geocoder in polygon

texts:
I have the TIGER geocoder installed, setup and working in POSTGIS.  I also see the reverse geocoder lets me get an address from a point on the map.

But now I want to provide a polygon and get all the addresses within the polygon.  Is this possible without trying to iterate over all possible points in the polygon?  (How would that even be accomplished?)

Best,
Matt


texts: Delete Features from a feature class that participates in a topology

texts: I have an enterprise geodatabase and I would like to delete features using the tool from a feature class that resides in one of the feature datasets and participates in a topology.  I tried using the tool, but it just churns and churns, ArcCatalog locks up and stops responding.  Is this possible to delete all the features, not the feature class if it participates in a topology?  It is part of a replication.  We get data from somewhere else and we just load in the new data to keep our schema.  The schema is the same for the new data.  I tried doing it from a version and also in ArcMap with locks and not responding.


texts:
**import arcpy
from arcpy import env
from arcpy.sa import *
import math
  # Check out any necessary licenses
arcpy.CheckOutExtension("spatial")
   # OverwriteOutput
arcpy.env.overwriteOutput = True
   # Set Geoprocessing environments
arcpy.env.autoCommit = "1000"
arcpy.env.scratchWorkspace = "D:\\NEO\\GISMODELLING\\BestLocationsModel.mdb"
# arcpy.env.extent = "214386.1775 3837166.9542 256099.9616 3907138.0305"
arcpy.env.extent = "D:\\NEO\\GISMODELLING\\BestLocationsModel.mdb\WasteWaterPlant\\Tart_UTM"
arcpy.env.cellSize = "10"
arcpy.env.mask = "D:\\NEO\\GISMODELLING\\BestLocationsModel.mdb\\WasteWaterPlant\\Tart_UTM"
arcpy.env.workspace = "D:\\NEO\\GISMODELLING\\BestLocationsModel.mdb\\WasteWaterPlant"
EucDist_Aggl = "D:/NEO/GISMODELLING/BestLocationsModel.mdb/EucDist_Aggl"
# Process: Get Raster Properties
Max_Aggl = arcpy.GetRasterProperties_management(EucDist_Aggl, "MAXIMUM", "")
print Max_Aggl
delta = int (Max_Aggl)
print delta
a = delta/1000
print a**

when I running this code I have this errore

  44250.27734375

  Traceback (most recent call last):   File "D:\NEO\test.py", line 28,
  in
      delta = int (Max_Aggl) TypeError: int() argument must be a string or a number, not 'Result'

I need to use the value of Max_Aggl in loops

I try to learn Python as I go


texts: Legal description polygon mapping/creation with ArcPy?

texts:
What I have:

an excel file with the names of properties within an area that is part of the PLSS in which Township/Range, Section and Lot numbers are used.  The file contains a legal description for every property.  Here are some examples of the individual legal descriptions:

NW/2 SW/2 - meaning the Northwest half of the Southwest half of a given section.
Lots 1,2,3,4 - meaning lots numbered 1-4 of a given section.
E/2 - meaning the Eastern half of a given section.

The excel data contains these fields.
Township_Range, Section, Legal_Description, Property_Name

I also have the Townships/Range, Sections and Lots as separate feature classes.

What I want to do:

I would like to write a script that goes through the excel table and adds features to a feature class based on the legal descriptions of every property.  I feel like this could be done by selecting the appropriate lots(divisions of each section) use arcpy.copy and append to the new feature class. However, the problem lies in creating a good way to select lots based on legal descriptions alone.  There are no fields within the lot feature class that provide a legal description.

after using a search cursor to narrow down the specific twp_rng and section,
It would somehow have to work like so...

if LEGAL_DESC == "E/2":
    select the lots that lie in the eastern half of the section    #How?
    copy selection
    append to new feature class

Obviously, this is a pretty dirty and could be optimized by having the legal description as an object instead of writing a separate "if" statement for each possible scenario, but that's not what i'm worried about.  The logic of how to select the correct portion of lots is the problem.

Anyone have any neat ideas??


texts: How can i project a coordinate system for many rasters all at once?

texts:
help
many rasters i have and should be projected in UTM coordinate system.
what should i do?


texts: Raster from Ascii text file doesn't project properly

texts:
I created an ascii grid of annual precipitation for CONUS in MATLAB and imported it into Arcmap. The text file can be downloaded here:

https://www.dropbox.com/s/0k4x257x597ogij/meanTotalPmm.txt?dl=0

After using the "ascii to raster" tool, I attempted to assign a coordinate system but without luck. Even though the correct information are recorded under the source tab in the layer properties, the file definitely doesn't project properly. As seen in the picture, it should align with the Oregon layer which is in the NAD_1983_Albers coordinate system. I used the "define projection" and the "project"  tool, and I attempted to edit the information from within ArcCatalog without success.

Does anybody know what I am doing wrong?


texts: Automated georeferencing based on topological point pattern matching in ArcGIS

texts:
I have a lot of building plans (PDF/PNG format) which I would like to automatically georeference. I have the base map "OS Mastermap Mono", which is use for a lot of the building plans.

If I use ArcScan to vectorise the PDF/PNG, is there a way for ArcMap to automatically georeference this layer as it shares the same background map?

To provide an example:

I load in the basemap OS Mastermap Mono.

I wish to georeference the following map, which is existing development:

And then overlay the proposed development:

The purpose is to have a GIS that visualises plans in a easier way than our current system which is just PDF documents, and allows us to calculate area of existing/proposed building space.

I hope this is clear as to what I'd like to eventually achieve, perhaps there is an easier way.


texts: R: mask a raster by a polygon also keep cells partially covered

texts:
I want to mask a raster by a polygon, the R script was as followed. However, some cells with only partially covered weren't included in the masked raster as listed in the picture link (https://drive.google.com/open?id=1EVehD8lYsR4oTr0SJPlhZhtti9b4WXxG), is there any method to accomplish?

library(raster)
library(sp)

r <- raster(xmn=1, xmx=5, ymn=1, ymx=5, nrows=4, ncols=4)
r[] <- 1:length(r)

Sr1 = Polygon(cbind(c(2,4,4,1,2),c(2,3,5,4,2)))
Sr2 = Polygon(cbind(c(4.1,4,2),c(2,3,2)))

SpP = SpatialPolygons(list(Polygons(list(Sr1), "s1"), Polygons(list(Sr2), "s2")), 1:2)

plot(mask(r, SpP))
plot(SpP,  add=TRUE)


texts: Plotting trajectory points on maps in python

texts:
When I read the paper Characterizing Urban Landscapes using Geolocated Tweets, I'm attracted by the figures illustrated in it, like this:

How do I plot these figures for my own data?
Which software or package could be used to achieve this?
It seems like the basemap package in python doesn't contain the POI layer.
My working language is python and OS is Linux, while the data(containing ID, Lng, Lat, timestamp) is represented using Numpy arrays.


texts:
I'm trying to use gdalwarp to reproject an ascii grid file. According to http://www.gdal.org/formats_list.html this shouldn't be a problem and indeed I can happily use gdal_translate to create grd files (AAIGrid) but not gdalwarp:

gdaltest 14:40:08 $ gdal_translate -of AAIGrid rj_DEM_1978.grd junk.grd
Input file size is 401, 301
0...10...20...30...40...50...60...70...80...90...100 - done.
gdaltest 14:40:25 $ rm junk.*
gdaltest 14:40:41 $ gdalwarp -of AAIGrid rj_DEM_1978.grd junk.grd
Output driver `AAIGrid' not recognised or does not support
direct output file creation.  The following format drivers are configured
and support direct output:
  VRT: Virtual Raster
  GTiff: GeoTIFF
  NITF: National Imagery Transmission Format
  HFA: Erdas Imagine Images (.img)
  ELAS: ELAS
  MEM: In Memory Raster
  FITS: Flexible Image Transport System
  BMP: MS Windows Device Independent Bitmap
  PCIDSK: PCIDSK Database File
  ILWIS: ILWIS Raster Map
  SGI: SGI Image File Format 1.0
  Leveller: Leveller heightfield
  Terragen: Terragen heightfield
  netCDF: Network Common Data Format
  HDF4Image: HDF4 Dataset
  ISIS2: USGS Astrogeology ISIS cube (Version 2)
  ERS: ERMapper .ers Labelled
  RMF: Raster Matrix Format
  RST: Idrisi Raster A.1
  INGR: Intergraph Raster
  GSBG: Golden Software Binary Grid (.grd)
  GS7BG: Golden Software 7 Binary Grid (.grd)
  PNM: Portable Pixmap Format (netpbm)
  ENVI: ENVI .hdr Labelled
  EHdr: ESRI .hdr Labelled
  PAux: PCI .aux Labelled
  MFF: Vexcel MFF Raster
  MFF2: Vexcel MFF2 (HKV) Raster
  BT: VTP .bt (Binary Terrain) 1.3 Format
  LAN: Erdas .LAN/.GIS
  IDA: Image Data and Analysis
  GTX: NOAA Vertical Datum .GTX
  NTv2: NTv2 Datum Grid Shift
  CTable2: CTable2 Datum Grid Shift
  ADRG: ARC Digitized Raster Graphics
  SAGA: SAGA GIS Binary Grid (.sdat)

gdaltest 14:40:52 $ which gdal_translate
/Library/Frameworks/GDAL.framework/Programs/gdal_translate
gdaltest 14:41:05 $ which gdalwarp
/Library/Frameworks/GDAL.framework/Programs/gdalwarp
gdaltest 14:41:10 $

As you can see, the format list is somewhat shorter than normal. From the two which commands at the bottom you can also see that both gdalwarpand gdal_translateare located in the same library.
Does anyone have an explanation for this? I can work around it easily enough but this shouldn't be happening, should it?


texts: Part of layer cut/disappears after exporting to pdf in QGIS

texts:
This is driving me crazy. The picture says it all:

You can see whole layer in QGIS browser and Composer, but when i export to pdf some layers are cut. I have tried moving the cut layers to the top of a layer tree, tried disabling labeling, tried to export them one by one, tried exporting to image and every time the result is same.
This somehow happens to only a few random layers.
Why is this happening?


texts: Assign multiple colors to features within a single vector layer

texts:
I want to give multiple colors to different features within a shapefile based on attribute values , is this possible in QGIS ?

If yes what is the procedure ?

Is there a way to provide color by means of hexadecimal codes ?


texts: Connect GeoServer WMS data to Leaflet.js

texts:
I have WMS Layer in GeoServer, I've made it following the walk through tutorial in https://www.e-education.psu.edu/geog585/node/701

I'm confused on how to connect my WMS to leaflet.js.


texts: Merging multiple QgsRasterLayer objects to create virtual raster (.vrt file) using gdal

texts:
I am creating multiple QgsRasterLayer objects in pyqgis and I need to merge them into a virtual raster .vrt file which will eventually be loaded as a layer in QGIS.
The QgsRasterLayer objects are being created by passing imagery urls.

From what I could gather, I need to use gdalbuildvrt as mentioned at https://www.gdal.org/gdalbuildvrt.html but there the .vrt file is being created through list of .tif files.

I also checked the solution at pyQGIS Merge raster layers but there too they are loading the raster layers through .tif files.

Do I need to create .tif files from all the raster layer objects (as mentioned in Exporting raster to rendered image using standalone script in QGIS?) and then use them to create .vrt or is there some other way to achieve this?


texts: Servers that implement OGC Geo/Time OpenSearch

texts:
The OGC is has published the "OpenSearch Geo and Time Extensions". That specification is based on, but extends, the OpenSearch Geo extension and OpenSearch Time extension.

As it turns out, I'm working on an open source (C#) implementation of the client side for OpenSearch. I'd like to include support for the OGC extensions to query by Geo and Time.

There are publicly available servers that implement the Geo and Time extensions, especially in the earth observation area (e.g. NASA Echo). However they don't implement the "extra" bits defined by OGC.

Are there any publicly available servers that implement any of the OGC extensions to the original OpenSearch Geo/Time extensions? In particular, I'm looking for implementations that expose the spatial and temporal containment operations ("relation" parameter).

Publicly available is preferred because:

Saves work in setting up a test server
Unit tests can be written that use it directly
Avoids confirmation bias in testing

However in the absence of publicly available, anything implementing the draft OGC spec (even "git branch" level maturity) would be better than what I have now.


texts: I am using the QGIS as a WMS client. I'm loading the WMS layer but the legend characters (which are in Persian language) are not displayed correctly in the Layer tab. Is there anyway I can handle the character encoding in QGIS?


texts: Vizualise different projections of a vector using QGIS

texts:
How to vizualise different map projection for a vector layer at once with QGIS,  then export the two forms as Geographic WGS84

e.g. I have a shapefile for Europe, how to vizualise it in CRS Lambert Azimutha Equal Area, Lambert Conic Conformal, and GCS WGS84 at once, to see the shifting and etc... then export them all as WGS84

Please see

PS: I know how to reproject layer and save them in any CRS or GCS


texts: Error with filterMetadata of FeatureCollection in Google Earth Engine Python API

texts:
I tried to filter a FeatureCollection based on a property i.e. where "Condition" == "Suitable", and count the number of features returned. This is the code I tried.
suitable_Features = irrigation_Cond.filterMetadata('Condition',"equals",'Suitable')
    suitable_Features.size().getInfo()

I got this error
HttpError                                 Traceback (most recent call last)
~\Anaconda3\lib\site-packages\ee\data.py in _execute_cloud_call(call, num_retries)
    333   try:
--> 334     return call.execute(num_retries=num_retries)
    335   except googleapiclient.errors.HttpError as e:

~\Anaconda3\lib\site-packages\googleapiclient\_helpers.py in positional_wrapper(*args, **kwargs)
    133                     logger.warning(message)
--> 134             return wrapped(*args, **kwargs)
    135

~\Anaconda3\lib\site-packages\googleapiclient\http.py in execute(self, http, num_retries)
    914         if resp.status >= 300:
--> 915             raise HttpError(resp, content, uri=self.uri)
    916         return self.postproc(resp, content)

HttpError: <HttpError 400 when requesting https://earthengine.googleapis.com/v1alpha/projects/earthengine-legacy/value:compute?prettyPrint=false&alt=json returned "Collection.reduceColumns: Error in map(ID=000000000000000006b6):
Number.gte: Parameter 'left' is required.">

During handling of the above exception, another exception occurred:

EEException                               Traceback (most recent call last)
<ipython-input-80-2f4a623ea954> in <module>()
----> 1 suitable_Features.size().getInfo()

~\Anaconda3\lib\site-packages\ee\computedobject.py in getInfo(self)
     93       The object can evaluate to anything.
     94     """
---> 95     return data.computeValue(self)
     96
     97   def encode(self, encoder):

~\Anaconda3\lib\site-packages\ee\data.py in computeValue(obj)
    676           body={'expression': serializer.encode(obj, for_cloud_api=True)},
    677           project=_get_projects_path(),
--> 678           prettyPrint=False))['result']
    679
    680

~\Anaconda3\lib\site-packages\ee\data.py in _execute_cloud_call(call, num_retries)
    334     return call.execute(num_retries=num_retries)
    335   except googleapiclient.errors.HttpError as e:
--> 336     raise _translate_cloud_exception(e)
    337
    338

EEException: Collection.reduceColumns: Error in map(ID=000000000000000006b6):
Number.gte: Parameter 'left' is required.


texts: QGIS 2.16 - Installation problem under Windows 10

texts:
I installed twice Qgis 2.16 using OSGeo 4w install. First I uninstalled all the previuos version and installed 2.16. Everything worked fine but when I open Qgis I have a warning message:

Couldn't load QGIS utils.
Python support will be disabled.

Traceback (most recent call last):
  File "", line 1, in
  File "C:/OSGEO4~1/apps/qgis/./python\qgis\utils.py", line 20, in
    from future import standard_library
ImportError: No module named future

Python version:
2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)]

QGIS version:
2.16.0-Nødebo 'Nødebo', d0b3e39

Python path:
['C:/OSGEO4~1/apps/qgis/./python', u'C:/Users/Gerardo/.qgis2/python', u'C:/Users/Gerardo/.qgis2/python/plugins', 'C:/OSGEO4~1/apps/qgis/./python/plugins', 'C:\\OSGEO4~1\\bin\\python27.zip', 'C:\\OSGEO4~1\\apps\\Python27\\DLLs', 'C:\\OSGEO4~1\\apps\\Python27\\lib', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\plat-win', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\lib-tk', 'C:\\OSGEO4~1\\bin', 'C:\\OSGEO4~1\\apps\\Python27', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\PIL', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\jinja2-2.7.2-py2.7.egg', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\markupsafe-0.23-py2.7-win-amd64.egg', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\pytz-2012j-py2.7.egg', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\win32', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\win32\\lib', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\Pythonwin', 'C:\\OSGEO4~1\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode']

I uninstalled and installed again but I get the same error. Qgis starts but I have nothing python related as expected. I know about this

Couldn't load QGIS utils; Python support disabled; No module named future

and also this

https://hub.qgis.org/issues/14577

The problem is the same even though the last link states the issue is closed.

Any ideas on how to solve the problem? some solutions on the first link do not apply to a win system


texts: QGIS 3x setting scale-based label opacity

texts:
I would like to decrease the opacity of my layer when zooming out.
I prepared something like this:
 Case
 When @map_scale>1750 Then (1 - transparency) * 255
 End

It doesn't work.
The similar thread is here:
How to set transparency per polygon in a shapefile in QGIS?
but applies to the shapes
How to make the scale-dependant label transparency?


texts: How to programmatically close a QMessagebox without clicking ok or X?

texts:
So part of my script in PyQGIS is this:

QMessageBox.warning(None,"Warning","bla bla bla")

The box apperas just fine but i want it to close after a few seconds so i am looking for something like this:

QMessageBox.warning(None,"Warning","bla bla bla")
time.sleep(2)
(command that closes the box without clicking ok or hitting X)

Any ideas?