texts: Creating lines from points preserving gaps in QGIS

texts:
I have CSV sets of navigation points, with the following scheme:
Easting, Northing, Line_name, Fiducial.

I can connect the lines with Points to Path, using the Line_name as the Group field and Fiducial as the Order field.

However, I also need to create multiple segments of the same line, if the points are too far apart.
Example:
We have points every 0.5 m along a line, but there are two instances, where the following point is 10 m away. As an end result, I need three line segments with no connection across the 10 m gap. Ideally, those line segments will preserve the Line_name in their attributes.

I need a tool that can batch process very large files with thousands of lines.

I was thinking of calculating the distance to the previous point and every time this distance is above a given critical limit, I would increment an indicator field, which I can then use to create subgroups, but this is a bit involved.

Is there an easier way?


texts: Can a QGIS plugin install its own svg symbols?

texts:
I am building a QGIS plugin that can request point features from a WFS, based on criteria defined in the plugins dialog.

Once the results are in, I want to display them with different symbols, based on a field called 'appropriate_symbol'. If a result point-feature has the field value 'fox', it gets a custom made icon of a fox, if the result has a field value 'elephant', it gets a custom made elephant-icon. And so on.

I know that QGIS comes with many svg symbols pre-installed. But I want to use customized symbols and ideally they should be shipped and installed together with the plugin. Is that possible? Can I execute a python command within the plugins __init__ that checks the QGIS installation for the location of svg symbols and adds my own symbols if they are not present already (how would that be implemented)? Or is there are more elegant way?

It is fine if the symbols are available to the user when the plugin is not activated. And it is a rather small user base that would install the plugin from a private repository. But it is important that their custom symbols are implemented, which I why I am looking for an elegant way to deliver them together with the plugin.


texts:
I am using Refactor fields from the Processing toolbox, since Table manager gives me a warning not to use it anymore.
I only want to change the order of the columns/fields in a small table.

Using Refactor fields gives the error below.
Can someone help me solve this?

"An error has occurred while executing Python code: See message log (Python Error) for more details.
File"C:/PROGRA~1/QGIS2~1.16/apps/qgis/./python/plugins\processing\gui\AlgorithmDialog.py", line 290, in closeEvent QgsMapLayerRegistry.instance().layerWasAdded.disconnect(self.mainWidget.layerAdded)
TypeError: 'instancemethod' object is not connected


texts: How to identify and simplify point clusters with regard to time in PostGIS?

texts:
I just started to work with spatial databases and I want to write a SQL(PostGIS) query for automatic generalizing of raw GPS-tracks (with fixed tracking frequency). The first thing I am wokring on is a query which identifies points of standstill in form of query like "x points within a distance of y meters" to replace massive point clouds by representative points. I already realized to snap points within a certain distance and count the snapped ones. In the picture below one can see a raw example track (small black points) and the centers of snapped points as colored circles (size = number of snapped points).

CREATE table simplified AS
 SELECT count(raw.geom)::integer AS count, st_centroid(st_collect(raw.geom)) AS center
   FROM raw
  GROUP BY st_snaptogrid(raw.geom, 500, 0.5)
  ORDER BY count(raw.geom) DESC;

I would be as quite satisfied with this solution, but there is the time-problem: Imaging the track as a full-day-track in a city the person can return to places already visited before. In my example, the dark-blue circle represents the person's home which he visited twice but my query of course ignores that.

In this case, the sophisticated query should only collect points with contiguous timestamps (or id's), so that it would produce two representive points here. My first idea was a modification of my query to a 3d-version (time as third dimension), but it does not seem to work.

Does anybody have any advice for me? I hope that my question is clear.

Thank you for the line-idea. I realized to make and simplify a linestring as you can see in the screenshot below (dots are original points).

What I still need is to determine the places of rest (> x points in < x meters radius), ideally as one point with an arrival time and a leaving time... any other ideas?


texts: Set parameter of field from number to string with Python

texts:
Since I realised Python is a very useful tool in ArcMap I don't want to do any workarounds any more. I have a field with is numeric, but I want to perform this calculation in the field calculator on it

def change(GRIDCODE):
    if GRIDCODE == "1":
      return "XZ"
    elif GRIDCODE == "3":
      return "XY"
    elif GRIDCODE == "4":
      return "XX"
    else:
      return "None"

But, the field is a numeric one and it won't let me change it to a string. I suppose I have to change the parameter of the field first. Of course I could add a field, set the parameter of that field to string and then run the code above, a little bit modified, on the shapefile.

But is there a solution where I change the parameter of the field within the code?


texts: Creating fishnet where each point from point file is located in middle of individual grid cells?

texts:
I have a point dataset in ArcGIS 10, and wish to create a fishnet, where each point is in the centre of an individual grid cell. Thus each point must be in the centre of a square.

I have used the Create Fishnet Tool, creating a fishnet with a cell size equivalent to the distance between my points. Then I have defined the projection of my fishnet to be the same as my point dataset.

My problem is that my points are located in the corners of the grid cells instead of in the centre. Is there any way of fixing this?


texts: How to have a drop down list for feature editing (attributes) in Quantum GIS?

texts:
I'm editing a lot of polygons of a shapefile in QGIS and i need to classify them, so i have to write the classes a repetedly a lot of times. In ArcMap there is an option to have domains set for the fields. How do i do this in QGIS?
Thanks


texts:
I'm trying to initialize the pyproj transformer as laid out here: Repeated transformations. This seems like it should be pretty easy. I'm trying to go from NC stateplane to lon/lat.

https://epsg.io/3358

https://epsg.io/6318

import pyproj
trans = pyproj.Transformer.from_crs(3358, 6318)`
Traceback (most recent call last):
  File "/home/MyUser/anaconda3/lib/python3.7/site-packages/IPython/core/interactiveshell.py", line 3331, in run_code
    exec(code_obj, self.user_global_ns, self.user_ns)
  File "<ipython-input-10-0c6306cf1b50>", line 1, in <module>
    pyproj.Transformer.from_crs(3358, 6318)
  File "/home/MyUser/anaconda3/lib/python3.7/site-packages/pyproj/transformer.py", line 311, in from_crs
    CRS.from_user_input(crs_from),
  File "/home/MyUser/anaconda3/lib/python3.7/site-packages/pyproj/crs/crs.py", line 440, in from_user_input
    return CRS(value, **kwargs)
  File "/home/MyUser/anaconda3/lib/python3.7/site-packages/pyproj/crs/crs.py", line 296, in __init__
    super().__init__(projstring)
  File "pyproj/_crs.pyx", line 2338, in pyproj._crs._CRS.__init__
pyproj.exceptions.CRSError: Invalid projection: epsg:3358: (Internal Proj Error: proj_create: cannot build projectedCRS 3358: cannot build geodeticCRS 4152: SQLite error on SELECT name, ellipsoid_auth_name, ellipsoid_code, prime_meridian_auth_name, prime_meridian_code, area_of_use_auth_name, area_of_use_code, publication_date, deprecated FROM geodetic_datum WHERE auth_name = ? AND code = ?: no such column: publication_date)

Something seems off, but I can't figure out what I'm doing wrong.


texts: Importing a shapefile and keep getting error with PostGIS manager

texts:
I am new to PostGIS and using pgAdmin. I have tried to import a shapefile into PostGIS manager and it gives me this message:

Importing with configuration: ex_hzfi99mczhvqbnkrya6ryz35mveyt_osm_point, public, geom, C:\Users\rhiem\Desktop\Database file\Slovenia\ex_hzfi99McZHVQbNkryA6Ryz35mveYt_osm_point.shp, mode=c, dump=1, simple=0, geography=0, index=1, shape=1, srid=0
Shapefile type: Point
PostGIS type: POINT[2]
Unable to convert data value to UTF-8 (iconv reports "Invalid argument"). Current encoding is "UTF-8". Try "LATIN1" (Western European), or one of the values described at http://www.postgresql.org/docs/current/static/multibyte.html.

Shapefile import failed.


texts: Connecting nodes to existing network lines using QGIS

texts:
I have road network (polyline layer) and point layer:

I am trying to connect nodes to the existing network lines. I able to connect nodes with lines using connect nodes to lines in Networks in QGIS as shown in first image.

The problem here is: It creates seperate line to connect node to line as shown with green circle in the image. I am trying to connect node in the existing line instead of creating new line as shown in Image 2.
I have searched existing answers to:

Joining nodes to network lines ( and not other points) using GRASS
Connecting points with lines using QGIS
Connecting points with line using virtual layer in QGIS

but the listed solutions are only for conventional point to line connections.

Image 2:


texts:
I am working on  a long riperian buffer that crosses the boundaries of 20 cities. My buffer is subdivided in the same layer, in 20 sections that correspond to the limits of the cities. Each section is divided into multiples sections that represents the recreational potential of the site. That makes roughly 300 polygons. In my attribute table, i have a field « ID » that is empty. I need to row that field according to the cities corresponding to my polygons.

This is what tried to do. At first, i am trying to row my « city » field. Then i try to row the « ID » field according to the order newly created by the row « City » command. It doesn't work. When i ask QGIS the request " City = @row_number " it mix everything up. My polygons that were, up to that point, accurately identified, are modified. For instance, a polygon identified to the city X is now identified to the city Y.

Do you know i could row all my polygons according to their cities ?

Ex :  ID 1 / City X / Section 1
      ID 2 / City X / Section 2
      ID 3 / City Y / Section 1

P.S Excuse me for possible errors in the translation. The QGIS I am using is in French.


texts: Create a polygon with a list of coordinates gave by QInputDialog in pyqgis console

texts:
I have a problem when i run this script. It gives me:

  TypeError: QgsGeometry.fromPolygon(unknown-type): argument 1 has unexpected type 'list'

Here is my code:

# -*- coding: cp1252 -*-
from qgis.core import *
from PyQt4.QtGui import *
from PyQt4.QtCore import *
from qgis.PyQt.QtCore import QVariant
from qgis.utils import *
from PyQt4.QtGui import QInputDialog #this is for your dialogs
#
qid = QInputDialog()
canvas = iface.mapCanvas()
#
N=0
X=0
Y=0
input, nba=QInputDialog.getText( qid, "Enter data", "N",
QLineEdit.Normal, "N" )
if nba:
   N=int(input.split()[0])
   print N
list=[]
q=1
while q<=N:
    input, ok= QInputDialog.getText( qid, "Enter data", " X,Y ",
    QLineEdit.Normal, "X" + "," + "Y")
#
    if ok:
        X=float(input.split( "," )[ 0 ])
        print X
        Y=float(input.split( "," )[ 1 ])
        print Y
#
    print "q",q
    print "N",N
    list.append('QgsPoint('+str(X)+","+str(Y)+')')
    q+=1
## cretae polygone
point=[]
print "Liste:",list
# space of entities
scale=10
rect = QgsRectangle(float(X)-scale,float(Y)-scale,float(X)+scale,float(Y)+scale)
canvas.setExtent(rect)
pt = QgsPoint(float(Y),float(Y))
canvas.refresh()
# Specify the geometry type
layer = QgsVectorLayer('Polygon?crs=epsg:4326', 'polygon' , 'memory')
prov=layer.dataProvider()
point= ', '.join(list)
point1="["+point+"]"
print 'point:',point1
#
#Point3 = [QgsPoint(100,100), QgsPoint(100,200), QgsPoint(200,200), QgsPoint(200,100)]
feat = QgsFeature()
feat.setGeometry(QgsGeometry.fromPolygon(Point1))
prov.addFeatures([feat])
# Update extent of the layer
layer.updateExtents()
# Add the layer to the Layers panel
QgsMapLayerRegistry.instance().addMapLayers([layer])


texts: Loading raster from PostGIS DB into QGIS: succeeded yesterday, fails now

texts:
I have a PostgreSQL(v.9.2)/Postgis(v.2.0.3) database which includes raster layers (added through raster2pgsql).

To open them in QGIS (v.1.8.0-Lisboa), the "Add Postgis layer" menu will not work; apparently it does not recognize the "raster" type as a geometry.
But I was able (following this little gem of advice) to add the raster layers to the canvas by opening the QGIS DB Manager, right-clicking the layer, and selecting "Add to Canvas".
Yesterday, all this worked fine.

Today, when I try this same procedure for the same layers in the same DB, nothing happens, they are not added to the canvas.
If I try to drag them onto the canvas, I get an error message "The layer is not a vaild layer and can not be added to the map".

The layers in the DB have not been changed in any way since yesterday (except perhaps by occasional Vacuuming of the whole DB).
They do have a "rid" and a "rast" attribute, and the latter is filled with long strings (looking similar to geometry values), thus probably still correct, at least not empty.

Dropping (deleting) one of the layers and recreating it from scratch (using exactly the same batch files as before) did not help; they were again flagged as "invalid layers".

I then created raster layers from the same data (and still using the same batch files) in a freshly created new DB, and here it all works again: both right-click > "Add to canvas", or dragging from the DB Manager will add the layer to the QGIS canvas just fine. Let's see for how long.

Is there any chance at all to fix the issue in the older DB?
If so: In which direction would I even have to look to find out what the problem is? DB corruption? QGIS settings? GDAL stuff?

Thanks for any pointers.


texts: Why are Postgis sequences not automatically updating in QGIS?

texts:
I'm adding new features to a polygon datasets in QGIS. The data is held in PostGIS but when i have a look at the attribute table for the polygon the gid field is not automatically updating, the entry reads NULL

I'm editing under the postgres login which looking at the sequeneces in PGadmin says that

ALTER TABLE parking_enforcement.restrictions_gid_seq OWNER TO postgres;

Is there anything special that i need to be doing in Postgres?


texts:
I am trying to create a symbol in MapServer. In this occasion it's a cross over a church. So I found a .png image of a black cross and include it in my symbolset.sym file:

SYMBOLSET
  SYMBOL
    NAME 'cross'
    TYPE PIXMAP
    IMAGE "cross.png"
  END
END

and in my mapfile:

 CLASS
  NAME "Church"
  EXPRESSION (("[building]" IN "church,chapel"))
  STYLE
      COLOR 208 202 208
      MAXSCALEDENOM 8000
  END
  STYLE
      SYMBOL "cross"
      SIZE 100
  END
END

The effects produced are weird though. My layer includes SIZEUNITS meters so the more I zoom in, the smaller the symbol get. And apart from that the symbol is repeated inside my polygon.

How can I maintain the size of the symbol and stop the repeating effect?


texts: MySQL normalizing points with latitude > 90

texts:
Facing very strange problem, I've points table for storing location of users in location(POINT) column. Problem with st_distance_sphere function where I've stored some specific location.

update points set location = point(30.492039, 76.587216) WHERE id = 1;

its Working fine but when I've changed lat,lng to

update points set location = point(52.051503, 113.471190) WHERE id = 1;

its called:

Incorrect arguments to st_distance_sphere

Table:

CREATE TABLE `points` (
  `id` int(10) NOT NULL,
  `name` varchar(45) ,
  `location` point,
  PRIMARY KEY (`id`),
)

Query:

    SELECT id,ST_Distance_Sphere(point(52.408089, 113.489844),location)/1000
AS distance FROM points HAVING distance < 100;


texts: osgeo4w not installing processing tools in QGIS toobox.

texts:
I have uninstalled and reinstalled QGIS several times from the OSGeo4W application.  I have done this with no problem on several other times, but now I cannot get the tools to load.  Most of the help simply say to just use the OSGeo4W application and not to worry about anything else.  What do you do when the OSGeo4W application is not working.

Do I need to change the Windows environment variables manually?  I have never had to do this before.


texts: Export Alberta Oil Sands Information Portal data from ESRI MapServer

texts:
I'm trying to get the Alberta Oil Sands outlines from https://maps.alberta.ca/genesis/rest/services/OSIP-Oilsands_Business_Data/Latest/MapServer

Until now I've tried the KML and KMZ export – both are empty. The JSON is ESRI JSON, which only contains metadata.

The online view and the WMS work – but these are just SVG shapes (js) or JPG tiles (WMS).

I thought that this would open Alberta goverment data. If not: How can I trace the outline to GeoJSON?


texts: Extracting relative risk values from a CARBayes CAR model in R

texts:
I'm working with the CARBayes package in R, modeling disease maps. I have a working convolution prior (BYM) CAR model, spatial_mod, which contains fitted values based on my model. My model takes on the form:

 O[i] ~ Poisson(mu[i])
log(mu[i]) = log(E[i]) + a[i] + bX1[i]+  b[i] +  h[i]

Where O[i] are my observed disease counts per region, and E[i] are my expected disease counts per region. The b[i] and h[i] terms are my structured and unstructured random effect terms.

This is the model I run in CARBayes:

form<-O~1+offset(log(E))

model_spatial<-S.CARbym(formula=form, data=data,family="poisson", W=W.mat, burnin=10000, n.sample=16000,thin=10)

I can extract the fitted values from my model as:

model_spatial$fitted.values

To find the RR for each region, i, would I simply exponetiate these fitted values?

exp(model.spatial$fitted.values)

In winBUGS, I know that the RR[i] for a given region will be exp(a[i] + bX1[i]+  b[i] +  h[i])

But, I don't know if this is true in CARBayes.

I apologize for not providing reproducible data.


texts: Error in GUID fields using "Simple Data Loader"

texts:
I have a replicated table in my SQL Server (PE_AMB_ANP), that table has a field named "rowguid" (rowguid is auto generated)

When I try yo load data, the field appear matched with 'none', but when the process continues it throws an error "duplicate key {0000-..}"

I don't want to update that field using ArcMap.

How can exclude "rowguid"?


texts:
I have digitized three rivers on Google Earth and converted into shapefile to use in MATLAB. I have checked the shapefiles in ArcGIS and they are working. But while opening in MATLAB, it shows the following errors:

Unsupported shape type PolyLineZ (type code = 13).
headerTypeCode = readHeaderTypeCode(shpFileId,callingFcn);

The code I am giving is below:
water = shaperead('rivers', 'UseGeoCoords', true);
geoshow( water , 'Facecolor', 'blue')


texts: How to restore a back up file into pgadmin III

texts:
My backup file containing about 50-60 tables , but I am able to just load only 6 tables for a new db and 31 tables for postgis dbase . And this was the error when I restore a backup file.

pg_restore: [archiver (db)] Error while PROCESSING TOC:
pg_restore: [archiver (db)] Error from TOC entry 1393; 0 0 SHELL TYPE box2d postgres
pg_restore: [archiver (db)] could not execute query: ERROR:  type "box2d" already exists
    Command was: CREATE TYPE box2d;

I should I resolve this problem. I am using PostgreSQL 9.1 and please help me how can I load all tables without any errors.my backup file contains 90 MB.


texts:
Can someone point me in the direction of some high spatial resolution and reasonably recent satellite imagery for central Canada (Manitoba)? I am trying to determine whether I can generate accurate data of tree lines for residential properties for land surveys. The current procedure is for the guys in the field to take a few ambiguous shots of trees and for the CAD guys interpret them - often causing issues. Does anyone know if it is possible to determine tree lines using remotely sensed data? I am hoping I can manipulate the data in ArcGIS to create polygons of cover, and then overlay the tree shots on to the surveyed data.

Key requirements:

Relatively recent (no good if the trees have been pulled down since image was taken!
High enough resolution for residential lot level
Free (well at least until I work out whether it is feasible)


texts: Programmatically Populate Value List in ArcGIS Tool Interface?

texts:
I'm just trying to set up the parameters in a tool I've built in ArcGIS. What I want to do is read in a feature class or table and get the values in a specific field called "SITE". The field contains a bunch of site names. What I want to do is read in the values of that field and populate them into the 'Value List' tool parameter property

I've been looking at the tool properties and I though I can use the 'Obtained From' property, but it looks like I can only obtain field names, but not that values in the field.

Is there anyway else to do this?

Is there anyways to read in the values from the SITE field and populate it to a value list programmatically?


texts: Python Library for converting GeoJSON Multi-polygon to polygon

texts:
I would like to know if there is a Python library that can be used to convert GeoJSON Multi-polygon to a Polygon.

Note: i want to convert .geojson file Multi-Polygon feature type to Polygon feature type.


texts: Can we create topographical map of private land using GPS elevation data?

texts:
Is there a way to create a relatively precise topographical map of a private land using GPS elevation data? As the land is a little bit over 100,000 square meters, the available regional topographical maps in Canada lack the precision I'd like to have for this project.

I have abundantly walked the land with a GPS and hence have thousands of points with elevation. Is there any way I could turn this into a somewhat interesting DEM or contour lines with QGIS?


texts: How to query relations and visualise it on Overpass turbo

texts:
I would like to know the reason the below posted code shows ways and nodes?
As shown below, I want to visualise only the relations but when I run the code on overpass turbo, it shows ways and nodes?
Code:
area["ISO3166-1"="DE"][admin_level=2]->.germany;
area["addr:city" = "Potsdam"](area.germany)->.potsdam;
(
   //node["sport" = "climbing"](area.germany);
   //way["sport" = "climbing"](area.germany);
   relation["amenity" = "university"]["addr:postcode"="14469"](area.potsdam);
)->.s;

out geom qt;


texts: Using FeatureLayer as pure client side query in ArcGIS API for JavaScript?

texts:
I've created this wrapper around feature layer where I create the empty feature layer as is in the API help page :

  new FeatureLayer(featureCollectionObject, options?)

providing the first parameter an object (not url) and some options as the second paramater. then I fill this feature layer with my data that I get from my service (not the esri map service btw)...

So the features are on the map (what a beauty) and then I want to query them, but it seems like this api doesn't do the client side queries with where clause... using queryFeatures or selectFeatures doesn't matter,  it always ends up with

  "FeatureLayer::_query - query contains one or more unsupported
  parameters"

I can create my own wrapper and loop through the feature atrributes, and then set them selected based on my search criteria, but I want feature layer to do it for me, without having to send a request to a url (which doesn't exist in our case).

Simple test, now with a proper feature layer (instantiated using the url) with the SNAPSHOT mode, one would expect that query with a where clause would be performed on the client side, but no, it actually sends a request to a server (what for when I have all the data on client?)

I've reused someone else's jsfiddle here -> jsfiddle .Just click in the map to run simple query ("resolution = 3") and you can see it fails.

Can anyone shed more light into this for me please?


texts:
I am currently evaluating the use of PostGIS (version 2.0) for:

a project that will store 3D coordinates (x, y, z) data
the data will have it's own local coordinate system which does not correspond to any datum (SRID)

Question: Can this data be stored without any SRID? If so, can I still take advantage of the 3D spatial query capabilities in PostGIS?

Thanks for you time,
Craig


texts: Create point geometry in a feature class with all the attributes prepopulated

texts:
So here is the problem. I have a text file with space delimited values. A standard text file with the headings separated by spaces and then on each next line values for these headings (columns) also separated with spaces. Each row represents a certain geographic point feature; however no coordinates (XY-pair) are stored in this text file.

My ultimate goal is to view each row in ArcMap individually (like one can do with the attribute table in ArcMap), edit the attribute values (if needed) and specify the point geometry of the feature by clicking on a map. And here is my issue – it seems as I cannot create a feature class without SHAPE field, but because I don’t have any coordinates in the source text file I cannot populate SHAPE field upon creating a new feature class.

What is the most efficient way of achieving this? I am a bit puzzled due to the workflow when all the attributes are populated yet the geometry is to be created. Ideally the solution would be to let user see the attribute table (obtained from the text file), then user selects a certain row user wants to work with, she digitizes a point on the map and the SHAPE field is populated with this point geometry. I understand that this is not how ArcMap was designed, but how one can come closest to this without investing too much in writing an advanced add-in or using ArcObjects?

I have considered so far:

creating a feature class from this text file with some random
coordinates, but it is not an applicable alternative because to edit
a point feature a user would need to zoom to this point and then move
it to another geographic area (using the ArcMap editing session).
Since the geographic extent of the points in this collection might be
quite large, it would be nearly impossible to use this approach.
creating a geodatabase table from the text file and then let user
click on the map to digitize a map point which would represent a
certain row in the table (however creating a point feature in a
separate feature class). Thereafter, user can choose to transfer all
the attributes of a geodatabase table row to the newly created
feature (was thinking of using arcpy.da.UpdateCursor here or running
a Python script with the join table & feature class > export logic).
However, this would imply that when user digitizes a point feature,
he would need to enter a custom ID that would match the ID of the row
in the table. Right now it seems as there is no maintained unique ID
heading in the text file and since multiple text files can be
received independently, this approach looks also pretty cumbersome.

Any other suggestions folks?


texts: Geocoding in 64 bits with background geoprocessing?

texts:
I have 10.1 SP1 installed plus the background geoprocessing patch for 64 bit geoprocessing.
Im trying to geocode a large file, so wanted to take advantage of 64 bits processing. My Surprise when i try to change properties of  locators  to make use of runtime memory limits to more the 2024 it gives me an error saying out of range limits. How can i process set up locators as described in this blog post

I also have enabled the background geoprocessing in ArcCatalog > Geoprocessing options and created a model to run the job, but when looking at taskmanager at run time its always a process of arccatalog (32bits) running and nothing related to 64 bits.


texts: arcmap mapping information of 2 layers, how to get rid of non-related results

texts:
I have 2 layers on the map: 1 drain line layer and 1 segment layer. The segment layer is covered by the drain line layer.

Now I want to get the corresponding attribute values on drainline layer which at the same place of the segment.

In order to do that, I do the following things:
(1) create a buffer at the segment layer
(2) create an intersect at the buffer and the drainline layer

After that I can read the mapping information of the darinline and the segment from the intersect. But there are extra feature that not from the segment also be involved in the intersect. The result of the intersect is as follows:

I don't want the The features within the blue circle because it not relevant to the segment layer.

Is there a way to avoid the useless feature from the drain line layer? Or is there other ways to mapping the data of the 2 layers?

I am using arcmap 10.2.


texts: Merging multiple rasters with the SAME spatial extent

texts: So I began with lidar elevation data 1m resolution, and I performed a terrain analysis to extract specific ranges of elevation from the raster (1st picture is 0-1m elevation, 2nd picture is 1-2m elevation). These extractions resulted in 2 different rasters, but for further analyses I would like to combine them back into one raster but now with just these two elevation size classes (I want the third picture but as a new raster I can process further). Mosaic is not working because the two new rasters are the same spatial extent. I have both arcGIS 10 and LIDAR Analyst.


texts: Resizing polygon features in ArcGIS Desktop based on attribute value?

texts:
I have a polygon dataset representing all the municipalities within a certain county.

Is there is a way to resize each feature in the dataset based on an attribute value?

Obviously, the features are naturally sized by their area, but I'd like to be able to resize them by - for example - population or median income for visual effect.

I know how to add and use the Scale tool on the Editor toolbar; this function is exactly what I want to do (i.e., resize / scale while retaining the shape), but I wonder if there isn't a more precise way to do this (i.e., based on the value of a certain attribute).

I know it's possible to do what I ask by selecting each feature individually then setting the scaling factor to whatever attribute value, but I'm curious if there isn't a more efficient way to do this.


texts: Creating different size service areas for single facility using ArcGIS Network Analyst?

texts:
When running NA Service Area analysis for a single point location (facility), I see a different service area polygon returned depending on whether Overlapping and Non-Overlapping polygons. For a single facility I would expect the same results back as there is no overlap likely to occur (for single extent of service area).

Is there a subtle difference in the processing which would mean even single facilities are processed  differently in one or other approach?


texts: List all overlapping polygons by objectID

texts:
How can I derive a table of overlapping polygons?

ArcPy
One feature class
Minimum Output: ObjectID, and Overlapping ObjectIDs

These don't quite do it

Count Overlapping Features - It doesn't identify which ones are actually overlapping


texts:
What data type do I need to create a LiDAR 3D surface model?
I have downloaded LiDAR/DSM data from Environment Agency, both as a WMS and as individual ASCII files, and it looks good but I want to view this data in 3D rather than 2D. Similar to what a Scene has by default in ArcGIS Pro for the Ground Terrain 3D surface (image attached showing what mine looks like). Basically I want a 3D DSM.

I have done this sort of thing in the past with other software such as QGIS, and even perhaps ArcMap, but I am struggling in ArcGIS Pro but I don't think it is a GIS problem, I am not convinced I have the right data or I am going about it the right way! I am by no means an expert in 3D and LiDAR data, but I do want to use it and understand it more.
I firstly tried to download data from the Environment Agency, and this gives me ASCII files, but I think these are basically raster files (?) and I probably need some sort of vector dataset, with a XYZ so I can plot each height value by it's attribute?
I am sure in the past I got LiDAR data as a shapefile or some other vector format, but when I look on EA website, or the Defra website (I have an account), then it only gives me options for WMS and so on, which I don't think are any good if I want to make this data 3D.
I have been looking into this for a couple of days now and I thought I would ask how other people are doing it. I have watched some online guides and people have created a LAS dataset from LAS files, but I can only seem to get hold of ASCII files or Web Map Services.
Not sure if I just need to find somewhere to download the appropriate file types, or if I am over-complicating the whole thing!


texts: Adding/viewing attachments to features using only ArcGIS DEsktop Basic license level?

texts:
Does anyone have an elegant solution to adding/viewing attachments, like a table/photos for example from within ArcMap using only the ArcGIS basic license level?  Or any suggestions?

I'm looking for a way to deliver an .mxd or file geodatabase with features that are clickable and link to tables of info and photos.  The above feature of 10.0 is great; however I am limited by the license we have.


texts: Is there a way to schedule a GDB on an SQL server to automatically compress once a day? As of right now I need to remember to open catalog and compress a Geodatabase manually everyday.


texts: How to process raster maps to get the same number of rows and columns?

texts: For modeling purposes, I need to use two raster maps that contain exactly the same number of rows and columns. How should I process my maps to obtain this result? For now, my maps don't have the same number of rows and columns...


texts: arcpy CopyFeatures_management fails when trying to copy annotation layers

texts:
I have a python code (using arcpy) trying to copy data into memory.
the data is an annotation layer (which probably contains special characters and/or chars)

I have the following code:

arcpy.MakeFeatureLayer_management(u'GIS.annotation', "TestFL")
arcpy.CopyFeatures_management("TestFL", u'in_memory\\___tst') # fails

the second line fails and throws an exception comnplaining about characters. throws something like:

  ERROR 001156: Failed on input OID 1735, could not write value
  '\xd7\x90' to output field TEXTSTRING.

Any ideas on why this happens or how to fix this?


texts: ESRI addProxyRule causing issues with back button functionality in Firefox

texts:
Currently, i'm using the addProxyRule to set up the proxy to be used when calling arcgis. I'm then taking that result and generating a url to take the user to a different part of our site pertaining to where they click on the map. What is happening though is that in Firefox, the back button does not take the user back to the previous page, it takes the user back 2 pages instead and doesn't seem to cache the previous page for the user to return to.

I've looked at the document.referrer in the console and it has the correct referring link so i'm not sure what is happening here.

the jquery code is below and i'm providing a link to the site (http://geospatial.com, if you click on the map on the page, it is calling the code below to direct the user to the country they select) where this issue is happening in case it helps solve this mystery.

function init() {

// setup proxy
var proxyUrl = "/desktopmodules/evg/evgstore/proxy.ashx";
esri.addProxyRule({ urlPrefix: "gis.geospatial.com", proxyUrl: proxyUrl });

$(document).ready(function () {

    if ($("#evg-home-map").length > 0) {
        // draw map

        map = new esri.Map("evg-home-map", {
            basemap: "gray",
            wrapAround180: true,
            nav: false,
            logo: false,
            slider: true,
            minZoom: 2,
            center: [20, 50]
        });

        dojo.connect(map, "onClick", mapClicked);
        dojo.connect(map, "onMouseOver", mapHover);
    }

});

function mapClicked(e) {
    var pt = e.mapPoint;

    // define query
    var query = new esri.tasks.Query();
    query.returnGeometry = true;
    query.geometry = pt;
    query.outFields = ["EVC_3_LETTER, C_NAME_CAPS"];

    // execute task
    var queryTask = new esri.tasks.QueryTask(ArcGisCountry);
    queryTask.execute(query, function (fs) {

        //esri.hide(loading);

        //// assign featureSet array to a single variable.
        if (fs.features.length > 0) {
            var url = "/products/countries/";
            var countryAbbr = fs.features[0].attributes["EVC_3_LETTER"];
            var countryName = fs.features[0].attributes["C_NAME_CAPS"];

                url = url + countryName.slugify("-") + "-" + countryAbbr.slugify("-") + "/";
                window.location.href = url;
        }

    }, function (errorObject) {

        alert("MapClicked Error: " + errorObject);
    });
}


texts: Silently installing v-late ArcMap extension

texts:
Is there a way to install the v-late ArcMap extension silently?
I'm also looking to deploy linkage mapper.
https://www.arcgis.com/home/item.html?id=69963ab422d04e649b64ac11cbadafed
https://circuitscape.org/linkagemapper/
I'm not sure how much more detail I can give.  Both extensions only come with zip files, no exe files, and instructions on how to do it in the graphic interface.  The v_late author suggested looking at ESRIRegAsm:  http://resources.arcgis.com/en/help/arcobjects-net/conceptualhelp/index.html#//0001000004n6000000
I actually caught these two commands running multiple times with sysinternals process monitor.  But it hasn't worked when I tried it myself.  The second command requires temp files anyway.
C:\Program Files (x86)\Common Files\ArcGIS\bin\ESRIRegAsm.exe" /p desktop /v 10.8 /s "C:\Program Files (x86)\ArcGIS\V_LATE_arcgis104\V_LATE_2_arcgis104.dll"

"C:\Program Files (x86)\Common Files\ArcGIS\bin\ESRIRegAsm.exe" /p desktop /v 10.8 /f:C:\Users\admin\AppData\Local\Temp\arc7CB7\10bD929.tmp /s "C:\Program Files (x86)\ArcGIS\V_LATE_arcgis104\V_LATE_2_arcgis104.dll"

It requires an elevated prompt.  If I could allow it without prompt elevation, that would be progress.  I can easy install ArcPro because it's an msi.  But these extensions seem impossible to deploy in a setting like a school classroom.
Is there a registry entry like this that would allow non-admins to install an extension like this?
HKLM\SOFTWARE\Wow6432Node\ESRI\Desktop10.8\Settings\AddInFolders

Other possibilities?
Put in C:\Users\<username>\Documents\ArcGIS\AddIns\Desktop10.8
Or run this:  C:\Program Files (x86)\Common Files\ArcGIS\bin\ESRIRegAddIn.exe c:\temp\myaddin.esriAddIn /s
And I don't know how to uninstall it for testing purposes.
How to install v_late by hand:
How to add the command:

Unpack the files:
V_LATE_2_argis104.tlb
V_LATE_2_argis104.dll
V_LATE_2_argis104.xml
V_LATE_Help.html
V_LATE_Help.jpg
V_LATE_categories.gif
Readme_V_LATE.txt

in the same folder.

Start ArcMap and choose "Customize Mode..." from the Customize Menu. You have to start ArcMap with administrator rights (right mouse click on the ArcMap icon and choose "Run as administrator"). This has to be done only once (allowing the registration of the extension).

Go to the "Commands" tab and choose "Add from file". Navigate to the folder where the "V_LATE_2_argis104.tlb" is located. Select the file (V_LATE_2_argis104.tlb) and click Open.

The Added Objects dialog box appears, reporting which new objects have been registered with ArcMap. Click OK. If you get the Message box "no new objects added", you may not have started ArcMap with administration rights, which is necessary for the first (!) installation.

In the "Categories" list on the left side click the new entry "V-LATE 2.0v_ArcGIS_104". Click and drag the "V-LATE 2.0" command from the Commands list and drop it on any (existing) toolbar.

Click Close

If the new Command is greyed out, you have to load a shapefile into ArcMap.

Select a polygon shapefile (feature classes from a geodatabase are not fully supported) and click on the "V-LATE 2.0" command


texts: Is there a method who break a shapefile into multiple files based on an attribute?


texts:
I'm new to GIS. I'm using the CMV template for my applications, but I can't get the Editor widget to work. The problem is that all the functionalities form the Editor Toolbar are not working. The only thing that works is adding a new point to the map, but then you can't delete nor edit it...
Here is my code:

return declare([_WidgetBase, _TemplatedMixin, _WidgetsInTemplateMixin], {

    templateString: template,
    i18n: i18n,
    widgetsInTemplate: true,
    editor: null,
    isEdit: false,
    mapClickMode: null,
    layerInfos: [],
    postCreate: function() {
        this.inherited(arguments);
        this.own(topic.subscribe('mapClickMode/currentSet', lang.hitch(this, 'setMapClickMode')));
        if (this.parentWidget && this.parentWidget.toggleable) {
            this.own(aspect.after(this.parentWidget, 'toggle', lang.hitch(this, function() {
                this.onLayoutChange(this.parentWidget.open);
            })));
        }
    },
    toggleEditing: function() {
        if (!this.isEdit) {
            var ops = lang.clone(this.settings);
            ops.map = this.map;
            ops.layerInfos = this.layerInfos;

            array.forEach(ops.layerInfos, function(layerInfo) {
                var controlOptions = layerInfo.controlOptions;
                console.log("ControlOptions: " + controlOptions);
            });

            var con = domConstruct.create('div', {
                innerHTML: '...',
                'style': 'text-align:center;'
            }, this.containerNode, 'only');

            require(['esri/dijit/editing/Editor'], lang.hitch(this, function(Editor) {
                this.editor = new Editor({
                    settings: ops
                }, con);
                console.log("click on new Editor => " + ops.layerInfos);
                this.editor.startup();
                this.editor.editToolbar.on('graphic-first-move', function(e) {
                    console.log('Move has begun');
                });
                ops.map.enableSnapping({
                    /*snapPointSymbol: symbol,*/
                    tolerance: 20,
                    snapKey: keys.CTRL
                });

            }));
            this.toggleBTN.set('label', this.i18n.labels.stopEditing);
            this.toggleBTN.set('class', 'danger');
            this.isEdit = true;
            topic.publish('mapClickMode/setCurrent', 'editor');
        } else {
            this.endEditing();
            topic.publish('mapClickMode/setDefault');
        }
    },
    endEditing: function() {
        if (this.editor && this.editor.destroyRecursive) {
            this.editor.destroyRecursive();
        }
        this.toggleBTN.set('label', this.i18n.labels.startEditing);
        this.toggleBTN.set('class', 'success');
        this.isEdit = false;
        this.editor = null;
    },

    onLayoutChange: function(open) {
        // end edit on close of title pane
        if (!open && this.mapClickMode === 'editor') {
            this.endEditing();
            topic.publish('mapClickMode/setDefault');
        }
    },
    setMapClickMode: function(mode) {
        this.mapClickMode = mode;
        if (mode !== 'editor') {
            this.endEditing();
        }
    }
});

This is the configuration for the widget:

editor: {
    include: true,
    id: 'editor',
    type: 'titlePane',
    path: 'gis/dijit/Editor',
    title: '<i class="fa fa-pencil-square-o"></i>&nbsp;&nbsp;' + i18nViewer.Editor_menuLabel,
    open: false,
    position: 7,
    options: {
        map: true,
        mapClickMode: true,
        editorLayerInfos: true,
        settings: {
            toolbarVisible: true,
            showAttributesOnClick: true,
            enableUndoRedo: true,
            createOptions: {
                polygonDrawTools: ['freehandpolygon', 'autocomplete']
            },
            toolbarOptions: {
                reshapeVisible: true,
                cutVisible: true,
                mergeVisible: true
            }
        }
    }
},


texts: Using ArcPy for Copy Parallel, Move, etc (from Editor toolbar)?

texts:
Since VB can activate the Editor, then use the tools listed in the Editor, such as copy parrallel, I am wondering whether Python can do the same job? I have been searching the website all day. Still can't find the answer from early questions-and-answers.

Can someone help me? Thanks in advance.


texts: CartoDB - can we overlay 2 visualizations

texts: Suppose we make two map visualizations from a loaded polygon shapefile in the CartoDB free online system.  Could we overlay these 2 maps on top of each other in a web page using the CartoDB.js API? The reason is, we need to make about 28 maps but in some cases the user will want to combine several maps.  Instead of making even more combinations we'd like to have transparent maps that could be overlaid.


texts: Calculate length of paths with defined vertexes

texts:
I'd like to calculate length of some paths with defined vertexes and each route is characterized by some records.

For example, I need the length of all combinations of lines (AB, BC, CD, AD, AC, BD).

Is there a Geoprocessing tool or QGIS Plug-in to do this, avoiding merging the records of the single path manually?


texts: Generate XYZ Tiles with Python not working

texts:
I am currently trying to generate XYZ Tiles with Python from a .tif file.
I tried using the QGis GUI and after some time the tiles were generating correctly.
After that I saw on another link, that I could see the processing command in the Processing History in the GUI. My folder output in the GUI is the following:

The corresponding history command in the GUI is:
processing.run("qgis:tilesxyzdirectory", {'EXTENT':'-443.962200000,456.037800000,-4759.144700000,-3659.144700000 []','ZOOM_MIN':11,'ZOOM_MAX':11,'DPI':96,'BACKGROUND_COLOR':QColor(0, 0, 0, 0),'TILE_FORMAT':0,'QUALITY':75,'METATILESIZE':4,'TILE_WIDTH':256,'TILE_HEIGHT':256,'TMS_CONVENTION':False,'OUTPUT_DIRECTORY':'C:\\Users\\xxx\\Documents\\FP\\QGIS EPSG','OUTPUT_HTML':'TEMPORARY_OUTPUT'})

After that I tried to reproduce the same output in python. The code itself runs without errors but the folder structure is complete different and it contains many negative Y folder and the code itself doesnt really stop.
The folder output generated from the python script is the following:

My python script for the tiles is:
import requests
import numpy as np
import wradlib as wrl
import sys
import matplotlib.pyplot as pl
import os

from qgis.core import *

QgsApplication.setPrefixPath('C:\\OSGeo4W\\apps\\qgis\\', True)
qgs = QgsApplication([], False)
qgs.initQgis()

import processing
from processing.core.Processing import Processing
Processing.initialize()

path_to_tif = "geotiff30101055.tiff"
rlayer = QgsRasterLayer(path_to_tif, "Standard raster layer")
QgsProject.instance().addMapLayer(rlayer, False)

processing.run("qgis:tilesxyzdirectory", {
    'EXTENT': '-443.962200000,456.037800000,-4759.144700000,-3659.144700000 []',
    'ZOOM_MIN': 11, 'ZOOM_MAX': 11, 'DPI': 100,
    'TILE_FORMAT': 0,
    'QUALITY': 75,
    'METATILESIZE': 4,
    'TILE_WIDTH': 256,
    'TILE_HEIGHT': 256,
    'TMS_CONVENTION': False,
    'OUTPUT_DIRECTORY': 'C:\\Users\\xxx\\Documents\\FP\\First', 'OUTPUT_HTML': 'TEMPORARY_OUTPUT'
})

I hope you can help me now because I didnt find a solution on my own so far.
Thanks in advance


texts: How do I cut a single line with polygons (in several segments) so that each cut line take data from the polygons? Something like a centroid takes.


texts: How do I automate the production of a map series?

texts:
I would like to produce a series of maps in which each map contains the same layers, with each layer of a given type having the same styling from map to map, but the source of the half of the layers differs from map to map. To be concrete, the map layers would be

common_source/layer1
common_source/layer2
common_source/layer3
map_specific_source/layer4
map_specific_source/layer5
map_specific_source/layer6

All the layers are vector layers.
Layer 6 is a 3x3 grid of cells. I need to select the middle cell and extract just that cell to use as the extent for exporting the map. Unfortunately that cell does NOT have the same object ID (or any other field value) that is constant from map to map.
Finally, I want to export the map, restricted to that middle cell, as a pdf.

The common source layers and styling can be stored in a template, so that part is easy. But how do I load the map specific layers from the correct source and style them?

I looked at the graphical modeler, but it doesn't seem to be designed for this sort of task. Is there another tool to look at, a way to do this in templates, or do I need to do this in Python?

Edit:
I should have mentioned that the map specific layers are grouped together in zip files. This means one approach would be make one map pointing to a zip called something like the_map.zip, then create each new map by copying each zip file in turn over the_map.zip, opening the project, and picking the correct center cell.


texts: Compare image in imageCollection with another image GEE

texts:
I want to calculate the number of monthly in a year that have less than 10% of the annual precipitation of that year, using the CHIRPS daily precip dataset on GEE.

My logic with code below (initiated for YY-1998 and annual precip for 1998 code not shown)-

1) Calculate the monthly precip for year YYYY and save it as an imagecollection by_monthYY.

2) Write a conditional statement to compare each image in by_MonthYY with annual precipitation of same year, such that if the by_MonthYY <=10% annual precipitation, a value of 1 is assigned, else 0, all saved in an other image collection by_MonthdryYY.

3) Reduce by_MonthdryYY by sum, to get the total number of dry months for that year.

I need to do above for multiple years and finally reduce by mean to get the average number of dry months over a time period.

var chirps_timeperiod = chirps_complete.filterDate('1998-01-01', '1998-12-31'); //YY=1998

var month_sum98 = ee.List.sequence(0, 1*12).map(function(n) { // .sequence: number of years from starting year to present
var start = ee.Date('1998-01-01').advance(n, 'month'); // Starting date
var end = start.advance(1, 'month'); // Step by each iteration

return ee.ImageCollection(chirps_timeperiod)
    .filterDate(start, end)
    .sum()
    .set('system:time_start', start.millis());
});
print(month_sum98);
var collection98 = ee.ImageCollection(month_sum98); //monthly precip for 1998

var replacement = ee.Image(1);
var conditional = function(image) { //image is map image collection of monthly precip by year separetly
return image.where(image.lte(0.1*(annualprecip_1998)), replacement);
};
var dryseason98=month_sum98.map(conditional);

However, I get an error saying image.lte is not a function. I also used image.ee.Filter.lessThanOrEqual with error image.ee is undefined. What function should I use to compare the two images?


texts:
How can I implement the Wrapper of PostgreSQL?

I found this useful looking wrapper from a FOSS4G routing with pgRouting Workshop but I do not know how to implement it!


texts: ArcGIS automatically convert string variables to numeric?

texts: If arcgis is automatically coverting string values to numeric variables in a text document generated by stata, how do you change those variables back into string? I tried opening up the fields under properties in arccatalog and it says I cannot change them because the tables are not empty


texts: Small unreadable toolbar icons of QGIS 2.16.3

texts: I have encountered with problem of small icons of toolbar of QGIS 2.16.3. I have ASUSS with high resolution (3200x1800). Changing resolution did not help. Also to apply solution by Dan Antonielli (http://www.danantonielli.com/adobe-app-scaling-on-high-dpi-displays-fix/) did not help as well. May be anybody solved the problem?


texts: Problem of "null" values added for fields with no data when appending data or adding new features to file geodatabase

texts:
I am working with a bicycle and pedestrian facility file geodatabase where I am continually appending new records to feature class attribute tables.  When I append new records, if there is no data for a certain field, ArcMap will automatically populate the records as "null" for that field.

Is there any way to prevent the records from being populated as null?  I am afraid that it will send the signal that there are errors in the database for database users.  I am considering using "Find and Replace" or possibly a script to replace the "null" values but wanted to see if there was a way to prevent this since I don't want to constantly have to delete the null records every time I append data or add features.


texts: Make Leaflet MarkerClusters stay open all the time

texts:
I use the Leaflet MarkerCluster plugin and when clustering is clicked, previous clustering closes, how can I solve this?
How can I keep it open, no matter which cluster is clicked?


texts: Disable editing of geometry in Geoexplorer, allowing only the editing of attributes for specific layer?

texts:
I wonder if, for a specific layer, it is possible in Geoexplorer to disable the editing of the geometry but still been able to edit the attributes related to the geometry (attributes edit popup).

I know that it's possible to exclude some of the attributes using excludeFields

(excludeFields
Array Optional list of field names (case sensitive) that are to be excluded from the property grid of the FeatureEditPopup.
http://gxp.opengeo.org/master/doc/lib/plugins/FeatureEditor.html)


texts:
I have a few .adf's that I aquired from a fellow using ArcGIS (as I understand, this format is a specific ArcGIS format).
I can open it up in qgis but their relative geographic position is off, seems the information is lost and all layers "upper left corner" are stuck at 0,0.
The cell sizes should be 2,-2 and 50-50, but all cell sizes for all rasters are set to 1,-1.
Is there a way I can manage these files in QGIS without losing this information?


texts: Identifying points on the intersection of three roads in PostGIS

texts:
I have two tables point and roads in my PostgreSQL 9.5 database. Some of the points lie on the intersections of three roads (line segments). The sample scenario is as follows:

I want to identify all such points that lie on the intersection of three roads. I came across this question on GIS-SE, but the code returned me all points on the intersections of individual line segments which is not desired. Can someone suggest me how do I identify the points that lie on the intersection of three segments (street crossing) only?


texts: How I can export (import) data from TAB format to MIF and vice versa?


texts: QGIS Saving grid to SpatiaLite database efficiency

texts:
I wanted to make vector grid in QGIS basing on canvas extent and save it to SpatialLite db.

Therefore, I am creating new database:

Layer > Create Layer > New SpatiaLite Layer > ... > Save.

Then, I am starting grid generation:
Processing > Toolbox > Vector Grid.
I am choosing shapefile extent as Grid extent, setting X and Y spacing and Grid type.

At the bottom I am choosing recently created database as output.
Process starts, and grid creates but it take long time in comprison with saving to shapefile.

Is there any way to speed up this process?

My QGIS version is 2.14 Essen


texts: Accessing the config.json in tileserver-gl-light from Ubuntu 20.04 server via command line

texts:
I am trying to access the config.json file for Tileserver-gl-light from a command line on an Ubuntu 20.04 Server. I am new to this operating system and just cannot seem to find the file. From my understanding I can change the info via sudo nano but I cannot find the file in my system. I have a working server using tileserver-gl-light, but intend to change the default map of Zurich to one that better suits my project.
This is the tutorial I used when setting up the server, only deviating to set python3 to work as python.


texts: Convert Areal Interpolation raster into polygons weighted by population (ArcGIS Pro)

texts:
Using the Areal Interpolation function in ArcGIS Pro, I have an interpolated raster layer of disease counts. This was derived from a district-level dataset.
I am now wanting to portion this raster into my other layer that has smaller areal census units - these have population attributes.
The automatic Areal Interpolation Layer to Polygons function does not give an option to weight the raster by an attribute when converting to polygons (I want to weight by population).

What would you recommend?


texts:
I have created a form in C# which collects data such as name, date, workspace and several other fields from the user and stores this data in a temporary .txt file to be used later by other functions.

I am trying figure out how to make an ArcMap MXD that when opened will automatically open and run this form and solicit this data from the user.


texts: Use Microstation wireframe model for OpenStreepMap

texts:
I've got an entire city in Microstation DGN format.
I'm able to convert it into Google Earth (KML), DXF, SVG and other formats, and then in .osm with ogr2osm.

However, the 3D model is build entirely with wireframes (triangulated polygon strings).

This does not seem suitable for .OSM.

Does anyone have experience with such problems and is it possible, and will it solve it, to change the model from wireframes into solid surfaces or solid blocks.


texts: How to get Esri *.gdb file into QGIS on Mac?

texts: I have tried the above options, but not being too computer savy, some I didnt get. I tried the homebrew option, but terminal says an error occurs when trying to install. I need this 1 file for a work project.


texts: Calculating the distance from each point in one layer to the nearest point in another layer

texts: I am trying to find the average distance from several points in a layer to points in another layer. I am using ArcMap Basic - I know there is a way to do this using the Advanced version of ArcMap but was wondering if anyone knew a workaround for Basic (if it exists).


texts: CSV file not joining to Shapefile using QGIS?

texts:
I am trying to join a .csv onto shapefile data. This csv file is a non-geometric file. When I go to join the two, the shapefile reads that the columns are null. Both files have matching columns.

I saved another file with the same name as a .csvt ("String","Integer","Real"), but nothing seemed to work. In my .csv I have a total of three columns: text, number, %.

I am operating on QGIS 2.8.

Shapefile:

CSV:


texts: Why can't an "AI export" from ArcGIS 10 be opened in CorelDRAW X5?

texts:
It doesn't matter whether the .ai file is 1.5GB or 5MB, whether it's complex or simple; I always get the "The file cannot be opened." critical error.

I use CorelDRAW X5 not Adobe Illustrator but there should not be any issues since CorelDRAW opens and saves to and from AI without any issues whatsoever. Except for .ai files exported from ArcMap. There seem to be no issues with the PDF export. It opens in CorelDRAW fine and I can manipulate and edit vector graphics. However I would like to be able to use the AI output for my publishing work. (Layers, named objects, etc.)

Any workarounds?


texts: Calculate field based on user defined parameters in ModelBuilder?

texts:
I want to use ModelBuilder to dissolve polygons by a field whose values are concatenated strings based on 3 other fields. The general process is this:

Add field, type string, called "concat". Field calculate that field to be [field1]&""&[field2]&""&[field3]. Dissolve polygons by field "concat".

However, I want to make field1, field2, and field3 user defined. How can I do this? In the field calculate box I can pick the three fields I want but these three fields could be different on another dataset. So when running the tool I want the user to be able to pick the polygon dataset, pick the three fields for the dissolve, then hit go.


texts: GeoServer ImageMosaic coming out blank with overlapping images

texts:
I'm trying to combine multiple "strips" of imagery into a single mosaic. Each of these strips is composed of many .tif files, and I have a directory with multiple subdirectories, each which contains one of these strips (big group of .tif files). The .tif files within any one of these strips do not overlap, but the strips themselves do overlap to varying degrees.
According to GeoServer's documentation, I should be able to point the ImageMosaic data store at the root directory to combine several strips. The exact overlap (I think) would be determined by the merge behavior parameter in the layer configuration, as well as the granule sorting parameter. So far, however, every attempt to mosaic all these strips' worth of .tif files together has failed, resulting in a blank image. I'm ultimately wanting to serve this mosaic with a WMTS - and the tiles do get "created" - but they're completely blank. I can request them, and I get an all-white image.
This happens whether I let GeoServer auto-create the shapefile index and properties file, or I make them myself.
I CAN get any of the individual strips' constituent .tif files to mosaic together, and I can serve those images via WMTS. The problem only occurs when trying to merge overlapping bundles of images. If I run the WMTS into ArcMap, zooming way out reveals nothing, while zooming in to certain areas seems to reveal a few actual tiles. I'm guessing there are the areas in which there were no overlapping pixels at all.
My only guess is that without a specific parameter in the shapefile index to reference for granule sorting, GeoServer just decides to serve up nothing in a given area. However, I did make my own ShapeFile index with an arbitrary ranking parameter to test whether I could get the sorting to work and imagery to show up. Still a no go.
Has anyone here had success with mosaicking overlapping imagery? Something major I'm missing?


texts: Proportional circles as function of value NOT range

texts:
I have two point layers, each containing one column of values.

When using the graduated style in layer properties, the size of the points is relative to the total range of values  within that layer . Classification (natural jenks) breaks the value range into certain intervals, they are assigned sizes, and that is what you get.

However, when you have different ranges, the circle sizes become meaningless when represented on the same map.  For example:

Here, we see the same circle sizes representing different value ranges.  How can I make circle size (or radius) in layer properties a function of the  actual value  assigned to the point, as opposed to classifying the values into intervals which are then assigned sizes?


texts:
i use mapinfo .tab for show the map, and i want to use SQl Server To select the other data, can mapserver have 2 connection in layer? or there are other ways this my .MAP layer configuration.

LAYER
    NAME layer_name
    TYPE POLYGON
    STATUS ON
    CONNECTIONTYPE OGR
    CONNECTION "./map/1172.TAB"
 >>>CONNECTION "connection sql"<<<
    DATA "select from mapinfo join to sqlserver"
    METADATA
      "wms_title" "Sabang"
    END

    PROJECTION
      "init=epsg:32646"
    END

    CLASS
        NAME "Elevation Poly 5"
        STYLE
            COLOR        232 232 232
            OUTLINECOLOR 32 32 32
          END
    END
END # Layer


texts: Acquiring pixel data from TIFF-file into CSV-file in QGIS

texts:
I have a TIFF-file which gives me a value of forest age per pixel in the whole Amazon basin.
I am trying to do a spatial-statistical analysis between forest age and climatic variables.
For this my goal is to have an attribute table that shows the forest age, climate variable, lat/long per pixel.
This TIFF file did not come with a text file or CSV, and I am not able to extract an adequate attribute table.
How can I acquire a CSV file with all the pixels values in QGIS for my study area? (age of forest, lat/lon, etc.)


texts: Display numeric scale in two different units, eg 1 inch = 1000 feet

texts:
How do I display a numeric scale for a map in two units, for example
1 in = 1000 ft ?

I know it's possible to display a scale bar as "Numeric," which provides only a ratio, in this example 1:12,000.

It seems like it should be possible to insert an expression into a label, only I can't figure out how to reference a map item in the expression builder. The map item has a variable "map_scale" but when I put this variable into a text item, it just displays as a string.


texts:
I am using GRASS r.skyview command (v.7, Windows) to determine "urban canyon effect" on pollution due to the presence of buildings.
The command nicely works, unfortunately I need to process a very large raster (representing buildings in an entire municipality) with a high resolution (10260x6762 pixels) and the command, which process every single pixel of the raster, run very slowly (about 24h to run on the entire area).

Since in the end I need the value of the skyview factor only for some specific points in space, I am wondering if it is possible to automatically do de following:

1) import in GRASS a vector layer representing point locations

2) select a single point

3) set GRASS region to a 100m x 100m square around selected point

4) launch r.skyview command

and loop through all points in vector layer.


texts:
We are using GeoServer to show a river/stream network as a tile layer.  Our data store comes from a shapefile, and the layer is styled to display the stream names as a label that follows the stream path. Tile caching is enabled, which is providing excellent performance.

One requirement is that users can click a stream and have it highlight, as well as return information about the stream. I'm currently implementing this as a GetFeatureInfo request that returns the feature's data as well as its geometry as GeoJson.  I'm then using the GeoJson to add a polyline in our google map, which renders on top of the tile layer.

The problem with that is that the polyline covers up the labels.  I would like instead to bring back the highlighted stream as part of the tile layer.

What I'm thinking of doing is having there be two layers, say streams_unselected and streams_selected.  Selected would be styled differently, say with a thick yellow stroke.  Then I would pass filters for each tile layer.  By default, i.e. without any stream selected, the streams_unselected layer would have a filter that always passes, and streams_unselected would have a filter that never passes.

Then, when the user selects a feature, I would modify the filters of each layer based on the feature_id (we'll call it X) of the feature selected.  The filter for streams_unselcted would specify that the feature_id != to X, and for streams_selected it would specify that the feature_id == X.

First question: Does this seem like a sane approach?

Second question: I'm currently using GeoWebCache, and specifying a URL to /geoserver/gwc/service/gmaps.  I don't think this will work with filtering.  I do love the performance, however.

I've seen posts online about GeoWebCache direct integration with WMS (http://geoserver.geo-solutions.it/edu/en/caching_data/caching_maps.html).  It looks like this may solve my problem.  Am I right in assuming that I can use a WMS GetMap url instead, passing my filters in the url, and GeoServer will use GeoWebCache in the background to serve tiles from that if they are cached, and generate new tiles if not?

I've also seen something on GeoWebCache Parameter Filters, maybe that's the better route.


texts: how to configuring geomajas wih oracle spatial db?

texts:
I have just done intial setup of the geomajas project and trying to config its data source with oracle database.
I am following the given example for sqlserver and POSTGIS.
My oracle layerroad xml file is like this:

<bean name="layerRoads" class="org.geomajas.layer.geotools.GeoToolsLayer">
              <property name="layerInfo" ref="layerRoadsInfo"/>
              <property name="parameters">
                     <list>
                           <bean class="org.geomajas.configuration.Parameter">
                                  <property name="name" value="namespace" />
                                  <property name="value" value="Oracle" />
                           </bean>
                           <bean class="org.geomajas.configuration.Parameter">
                                  <property name="name" value="user" />
                                  <property name="value" value="mdsys" />
                           </bean>
                           <bean class="org.geomajas.configuration.Parameter">
                                  <property name="name" value="passwd" />
                                  <property name="value" value="mdsys" />
                           </bean>
                           <bean class="org.geomajas.configuration.Parameter">
                                  <property name="name" value="database" />
                                  <property name="value" value="mdsys" />
                           </bean>
                           <bean class="org.geomajas.configuration.Parameter">
                                  <property name="name" value="host" />
                                  <property name="value" value="localhost" />
                           </bean>
                           <bean class="org.geomajas.configuration.Parameter">
                                  <property name="name" value="port" />
                                  <property name="value" value="1521" />
                           </bean>
                           <bean class="org.geomajas.configuration.Parameter">
                                  <property name="name" value="dbtype" />
                                  <property name="value" value="Oracle" />
                           </bean>

In web.xml file I have written this

WEB-INF/example/geotools/oracle/layerRoads.xml

In pom.xml

<dependency>
                     <groupId>org.geotools.jdbc</groupId>
            <artifactId>gt-oracle-spatial</artifactId>
            <version>2.5.5</version>
              </dependency>
              <dependency>
                      <groupId>com.oracle</groupId>
               <artifactId>ojdbc14</artifactId>
              <version>10.2.0.2.0</version>
              </dependency>

When I am running it it gives me this error

[Console scanner] WARN  o.g.layer.geotools.DataStoreFactory - Datastore factory PostGIS(PostGIS Database) is not available
[Console scanner] WARN  o.g.layer.geotools.GeoToolsLayer - The layer could not be correctly initialized: layerRoads
ception: No datastore found. Possible causes are missing factory or missing library for your datastore (e.g. database driver).
Available() method of your DataStoreFactory class to find out which libraries are needed.
factories : PostGIS


texts: Where is layer stacking from merge raster tool of qgis?

texts:
In the newest version of Qgis (2.18) there isn't layers stack option in order to merge raster bands. So tried to merge Landsat bands in only one file and I couldn't.

What can I do to stack bands of landsat in only one file in Qgis 2.18 now?


texts: Loading an EM Survey Data into QGIS (des, dfn,dat and gdb file sets)

texts:
I have downloaded an EM survey data pack recently and within the zip file it contains:

Final_EM.des - header file describing contents of the located data
Final_EM.dfn - Channel definition file containing info on located data in the *.dat file
Final_EM.dat - Flat ascii file containing located EM, magnetic and digital terrain data
Final_EM.gdb - Geosoft database file containing located EM, magnetic and digital terrain data

I am trying to upload the EM data files onto QGIS so that I can visualise it against some other geological data files.
I have trying dragging and dropping, tried uploading as a vector layer, raster, it does not seem to be able to load up.
Just wondering if there is anyone here that could help me out?


texts: Cut road polygon layer into small polygons?

texts:
I need a model to split that road stream at the intersections to have a road polygon to every street. I have tried to buffer the streets and use the buffer to cut the road polygon but it failed and I don't want to do it manually by the cut polygons have more than 1200 streets here it will take time and effort.

I want to split the purple polygon to smaller polygons that surround every street with its road polygon.

I am using ArcGIS 10.3.


texts: Supervised Classification within the bounds of a FeatureCollection

texts:
I am a bit new to classification within GEE and was hoping to get some help regarding classifying within a specific region.

var bands = ['NDVI'];

var landSat = ee.Image(ee.ImageCollection('LANDSAT/LC08/C01/T1_8DAY_NDVI')
  .filterBounds(roi)
  .filterDate('2014-5-01', '2014-5-10')
  .select(bands)
  .sort('CLOUD_COVER')
  .first());

var landSatTest = ee.Image(ee.ImageCollection('LANDSAT/LC08/C01/T1_8DAY_NDVI')
  .filterBounds(MY FEATURE COLLECTION POLYGONS)
  .filterDate('2014-5-01', '2014-5-10')
  .select(bands)
  .sort('CLOUD_COVER')
  .first());

//Merging imports of sample regions into one feature collection
var trainingFC = not_fallowed.merge(fallowed);

//Training data
var training = landSat.select(bands).sampleRegions({
  collection: trainingFC,
  properties: ['landcover'],
  scale: 30
});

//Random Forest Classifier Training
var trained = ee.Classifier.randomForest().train({
  features: training,
  classProperty: 'landcover',
  inputProperties: bands
});

//Run classification
var classified = landSatTest.select(bands).classify(trained);

//Displaying results
Map.centerObject(iqROI, 7.5);

//70FF00 (green) not fallowed, FF2D000 (red) fallowed
Map.addLayer(classified, {min: 0, max: 1, palette: ['70FF00', 'FF2D00']},
'classification');

So rather than defining a new ROI within landSatTest I would like to apply it to my FeatureCollection in order to only classify that area as I am only interested in that area's outputs. Otherwise it seems I get an area that adds a lot of extra space.


texts: View HTML field data with QT4 browers as a QGIS Run Feature Action Tool

texts:
I've successfully populated a sqlite table with an HTML-encoded-table as a field attribute, and so there is literally a Table within a single field value in sqlite, for each record.  This is going to join with spatial features in a spatialite table to replicate a One polygon to Many records relationship.  Map Tips can already successfully pop up the table when the mouse moves over the feature.

Now I'd like to increase functionality so that I can copy from the table, and scroll down through larger tables. I want it to open up in a browser box ( The largest tables don't show properly because the Map Tips pop-up window has its top margin aligned outside the monitor space, and the html table itself is aligned to the top of that map tips pop-up widow -- which I'm not worrying about trying to resolve )

I'm looking at a QGIS-Python-QT4-Brower example from this site: http://linfiniti.com/2012/03/a-python-layer-action-to-open-a-wikipedia-page-in-qgis/
titled "A python layer action to open a wikipedia page in QGIS"
I saw that I might be able to do what I want with Python and QT4-Brower Window that's already a part of the QT4 packages that come with QGIS.  I just need to run an Action on the field that contains my HTML table and have it pop up in that QT4 Browser, which would be more versatile than the Map Tips pop-up.
The only difference with what I want to do is that I want to view the html-data in my field with the Qt4-browser rather than have the Qt4 browser open a designated wikipedia page.  But I can't figure out how to pass the data to the browser.. I need it to already be in a file.

My knowledge of Python is limited to basic scripting tasks, iterating through Lists with cursors etc.  I figure the "A python layer action to open a wikipedia page in QGIS" might provide part of the code that someone else could more easily come up it. Does anyone know how to pass data into a Qt4 web-browser, rather than opening a file with the data already in it?

UPDATE:
I'm now very close to solving this by scripting Python to write output to a myTable.html file.  But, when passing the string of data from my attribute table field [% MyFieldName %] I get an error

PYTHON

mytable = [% MyFieldName %]

(that causes: "  TypeError: expected a character buffer object")

I can't pass the string of data from the [% MyFieldName %] to the mytable variable without a Python Error box popping up instead of Qt4 Browser Renderer, but using this test script below (as a single line) I can get the Qt4 Browser to popup with the text "my html table will go right here":

mytable = "my html table will go right here";
tabfile = open('C:\MyFolder\myTable.html', 'w');
tabfile.write(mytable);
tabfile.close();
from PyQt4.QtCore import QUrl;
from PyQt4.QtWebKit import QWebView;
myWV = QWebView(None);
myWV.load(QUrl('file:///C:/OSGeo4W/myTable.html'));
myWV.show()

The funny thing is that by changing the value of the mytable variable above to
= [% MyFieldName %] causes the QGIS Action to throw up a popup PYTHON ERROR Dialog Box THAT ACTUALLY RENDERS MY HTML STRING INTO A TABLE PROPERLY, JUST LIKE I WANT IT, but followed by syntax error message and other details of the program and entire script.  And so I know that [% MyFieldName %] does pass the data from QGIS (to atleast) the Python Error Dialog box  which can render html tables  :-/

Now when I only do a one line test script setting a Python Variable to  someXYZ = [% MyFieldName %] then that generates the TypeError: expected a character buffer object


texts: Is it possible to focus on different (but related) features on the same Atlas page?

texts:
I have a map showing a city and its divisons. I have then several features in other layer with the style set to show only the feature designed by the atlas by $id = $atlasid

So far so good. Then I want to create a print composer with two maps: the first one will be set to be "controled by atlas" with 10% margin around it (this I know how to do). Then the second one need to show the same feature but need to be centered by the city division that it belongs to. In other wors I need to input a formula that says whitch layer that map need to use to center the image. (like in the image)

Any ideas, Plugins or solutions?


texts: How to refer to another table in a CASE WHEN Statement in PostGIS?

texts:
I'd like to set column3 in  table1 to a specific value by checking if a certain threshold (Queried by intersection) is exceeded. My first table holds geometry information with additional columns,the same is true for my second table. I finally ended up with this code:

UPDATE table1
    SET column3 = CASE

                    WHEN t.x >15 THEN 1
                    ELSE 0

                  END

    FROM LATERAL (SELECT table2.column7 AS x FROM table2 WHERE ST_INTERSECTS(table1.geom, table2.geom)) AS t;

But I get a message: ERROR:  invalid reference to FROM-clause entry for table "table1",

with a hint: HINT:  There is an entry for table "table1", but it cannot be referenced from this part of the query.

What is wrong about my syntax?


texts: Calculate area of adjoining raster class cells inside polygon

texts:
I am calculating the size of deforestation in private properties. My two data sets are a raster layer of deforested patches, each containing a number as the year it has been deforested, and a polygon layer of the region for the properties.
Using Tabulate Area I have the total area deforested each year on each property.
I would like now to analyse the size of the polygons deforested, ie. , calculate area of raster cells of same value that are adjoining on each property.

Is there a way to accomplish this on ArcGIS?


texts: Sharing local *.nmf's that are stored on separate drive for ArcGIS Explorer Desktop?

texts:
I am trying to use ArcGIS Explorer to give administrators basic access to maps that are stored in a drive that is specifically for GIS use. We don't want the data stored on the drive that they use because of the possibility of them accidentally altering the maps. We have no server access and cloud services are out because of sensitive data.

I have application configurations set up and linked to the map on my network drive and set the ArcExplorer map to store relative paths.

The problem is that when the map is opened through Explorer, it always opens as "default map" and all the layer file links are broken. I have tried changing the order in which I set it all up but so far, the results remain the same.


texts: Merging Attributes of points into polygon using ArcMap?

texts:
I am using ArcMap 10.

I have points layer which the attribute looks like this:

FID Shape ID   Lat       Long
0   Point 2065 -5.51666  106.48333
1   Point 2319 21.93333  120.73366

I also have polygons layer :

FID Shape ID
0   Polygon 0
1   Polygon 0

How could I join both attribute file so that it become

FID Shape ID   Lat       Long          FID2    Shape

0   Point 2065 -5.51666  106.48333     1       Polygon
1   Point 2319 21.93333  120.73366     1       Polygon
2   Point 2304 22.93333  129.73366     3       Polygon

The actual case is, there are points that fall within one polygon.


texts: arcpy symbology colors changing on RefreshActiveView()

texts:
I have a problem that I'm not sure is functionality not built into arcpy or I'm just doing something wrong. I'm using arcpy.mapping.UpdateLayer() to apply symbology to a layer, but afterwards some elements of the symbology needs to be removed depending on the attributes of the value field. I've achived this programmatically, determining which symbology values need to be removed. I'm using something similar to the following code to remove undeeded values from the symbology.

lyr = arcpy.mapping.Layer(<layer name>)

a = lyr.symbology.classValues
b = lyr.symbology.classLabels

a --> [u'1', u'2', u'3', u'4', u'5']
b --> [u'Low', u'Moderate', u'High', u'Very High', u'Extreme']

a.remove('5')
b.remove('Extreme')

a --> [u'1', u'2', u'3', u'4']
b --> [u'Low', u'Moderate', u'High', u'Very High']

lyr.symbology.classValues = a
lyr.symbology.classLabels = b

arcpy.RefreshActiveView()

The method seems to work great. The problem is the colors set by the symbology change after the arcpy.RefreshActiveView(), which defeats the whole purpose of applying symbology in the first place. Is there a way to keep the colors the same after removing symbology values?

Thanks in advance!


texts: Is it be possible to re-render parts of a tile pyramid when postgis table updates?

texts:
I have this idea of using postgis together with a tile renderer for serving my organization with tiled maps. As I believe my data will be partially altered almost on a daily basis it does feel like a performance issue if I have to render the whole cache every night as the changes are quite small but still significant.

Would it be possible to re-render parts of a tile pyramid or a grid in a cache instead of a complete re-rendering?

And have anyone tried to execute this rendering in for example mapnik from inside postgresql/postgis using a database trigger? I believe the first issue would be to determine what tiles are affected and then re-render them specifically.

Any ideas?


texts: Unable to reference GPS waypoints correctly in ArcGIS

texts:
I'm working on a project. Waypoints were taken in the field with a Garmin Dakota 20.
I've tried two different ways of putting them into ArcGIS. One way is through converting the GPX file to a shapefile via DNRGPS. This results in the coordinate system being wrong (even if changing it in the program) and the waypoints end up at the wrong place. The other is with the data interoperability tool which results in the waypoinst missing spatial reference. Even after adding such in through the catalogue this doesnt move the Waypoints anywhere. What am I doing wrong? I want them matched to the current data frame with its maps (Sweref99)

EDIT: After editing once again through the ArcCatalog i tried changing it to Geographic Coordinate System instead of the projected (still Sweref99) and now, instead of them being nowhere, I get them at the same spot I did through dNRGPS (quite a bit off, but at least on the map).


texts: How to add multiple python expressions in field calculator

texts:
How would I add multiple python expressions in the ArcMap Field Calculator if possible?

I tried separating each expression on a new line and also separating with a semi-colon but neither worked.

Here are the expressions I am trying to execute at once:

!field!.replace("+","")

'_'.join(num.zfill(4) for num in !field!.split('_'))

!field![:9]


texts: cutting waterlines at the edge of polygons

texts: Using ArcGIS 10.3, I would like to cut the river lines at the edge of water body polygons bordering the region. Using the intersect tool and split line to points tool resulted with the disappearance of the lines located outside the water body polygons. Any idea how to split these lines at the border without erasing them?


texts:
How to pass slope attribute to the command

   "gdaldem slope manj.tif outpiu.tif -p use percent 45 -s 111120/370400"

  It's giving error like FAILURE: Too many command options.

    give any working example for my understanding perpose.


texts: What is the range of values of Sentinel-2 Level 2A images?

texts:
I found in Level 2A Product Format Specifications Technical Note page29 that:

"The L2A Quantification Value is aligned with the L1C Quantification Value of the L1C product from which the L2A product is generated"

When I refer to Sentinel-2 Products Specification Document page 403 I found the following:

Why would the range of values of Level 1C images be 1 to 65535 if the images are supposed to be 12 bits. And what is the range than of Level 2A, is it 0 to 10000?


texts: Installing Opengeo Suite on Ubuntu 13.04 - " trying to overwrite '/usr/bin/shp2pgsql-gui' "?

texts:
I am trying to upgrade Opengeo Suite from 2.0 to 4.0, specifically for the upgrade of PostGis 2.0 to 2.1, I followed the instructions outlined on this page: http://suite.opengeo.org/opengeo-docs/installation/ubuntu/upgrade.html
But had no success.

I realized I hadn't fully installed opengeo-client or postgis 2.1, so I used apt-get install and received this:

The following packages have unmet dependencies:
opengeo-client : Depends: postgis-2.1 (>= 2.1.0) but it is not installed
postgresql-9.3-postgis-2.1 : Depends: postgis-2.1 but it is not installed
E: Unmet dependencies. Try using -f.

So I used -f and received..

The following extra packages will be installed:
  postgis-2.1
The following NEW packages will be installed:
  postgis-2.1
0 upgraded, 1 newly installed, 0 to remove and 3 not upgraded.
28 not fully installed or removed.
Need to get 0 B/701 kB of archives.
After this operation, 1,839 kB of additional disk space will be used.
Do you want to continue [Y/n]? Y
dpkg: warning: files list file for package 'postgresql-9.2-postgis2' missing; assuming                     package has no files currently installed
(Reading database ... 384786 files and directories currently installed.)
Unpacking postgis-2.1 (from .../postgis-2.1_2.1.0-1+opengeo_amd64.deb) ...
dpkg: error processing /var/cache/apt/archives/postgis-2.1_2.1.0-1+opengeo_amd64.deb (--unpack):
 trying to overwrite '/usr/bin/shp2pgsql-gui', which is also in package postgis 1.5.3-2ubuntu1
dpkg-deb: error: subprocess paste was killed by signal (Broken pipe)
Errors were encountered while processing:
 /var/cache/apt/archives/postgis-2.1_2.1.0-1+opengeo_amd64.deb
E: Sub-process /usr/bin/dpkg returned an error code (1)

I was recommended to go back and forth with apt-get upgrade and apt-get -f install to isolate the issue. But I am assuming that Postgis 1.5 is still in use, but I cannot find where to remove it. (If it is possible to do a complete wipe of the Opengeo suite then I can also do that if it would be simpler, all of my data is backed up)

I am not experienced in Linux. Most of my work is on the GIS/Web side.

My Ubuntu release is: Ubuntu 13.04 Raring


texts:
This might not be strictly GIS related, but here goes.

I am looking for a script (DOS based) or application (Windows) that does the following:

In a filesystem there are a collection of raster maps, for instance 2528aa_2008_ed1.tif, which gives the geographic location of the map, the year issued, and the edition.

We then receive new maps, so that for instance we now receive a map that is 2528aa_2010_ed3.tif.

The script or application must look at the file names in the folders and replace 2528aa_2008_ed1.tif with 2528aa_2010_ed3.tif, so that the folder now only contains the new updated map.

Has anyone had to do this, or can anyone give me guidance on scripts or applications out there that can do this.


texts: Unable to connect to database within Conda Data Cube environment

texts:
I followed the steps listed at https://cubes.readthedocs.io/en/v1.0.1/install.html
and did the following:

Installed Python

Installed PostgreSQL
ssudo sh -c 'echo "deb http://apt.postgresql.org/pub/repos/apt $(lsb_release -cs)-pgdg main" > /etc/apt/sources.list.d/pgdg.list'
wget --quiet -O - https://www.postgresql.org/media/keys/ACCC4CF8.asc | sudo apt-key add -
sudo apt-get -y install postgresql

Installed Miniconda and created a conda datacube environment
conda config --add channels conda-forge
conda create --name cubeenv python=3.6 datacube
source activate cubeenv

datacube system check

Gives me an error stating

psql: error: could not connect to server: could not connect to server: No such file or directory
Is the server running locally and accepting
connections on Unix domain socket "/tmp/.s.PGSQL.5432"?

Problem is that PostgreSQL is accessible outside the environment but not getting connected once the environment is activated.
A screenshot of the same is attached.


texts: QGIS missing documentation about GetFeatureInfo result template

texts:
How/where to change HTML of GetFeatureInfo result?

GeoServer's documentation has GetFeatureInfo Templates about changing GetFeatureInfo result template.

ArcGIS's documentation has Customizing a WMS GetFeatureInfo response about changing GetFeatureInfo result template.

Could not find appropiate information in QGIS documentation. Found someting about "writing a filter" in this page. It is mentioned here that Currently, the HTML response template is hardcoded and there is no way to customize it and here I did not see any mentions about GetFeatureInfo response template

For example, I want to remove first - layer name - row from the table.


texts:
I am having a hard time finding any detailed documentation defining ESRI's JSON specs.  I am hoping someone can shed some light on the following two questions.

Similar to the GeoJSON spec, the first array in the 'paths' collection is always an exterior ring, and all subsequent arrays in this collection are interior rings (holes)?
Is above is true, how does the ESRI JSON spec deal with multi-polygons?


texts: Publishing GeoTIFF using WCS service with QGIS Server?

texts:
I need to publish some GeoTIFFs using QGIS Server (2.2.0-1+precise1).  I already have installed and setup QGIS Server to serve those geotiffs via WMS.  But I need to serve rather the raster data to be able to use it in raster calculations.

I don't know how to set up or call the service from an URL.  I tried with the parameter SERVICE=WCS instead of SERVICE=WMS, to no avail.

Can QGIS server handle WCS?


texts: Looking for an OpenEV replacement, i.e a fast geotiff viewer

texts:
I have in the past used OpenEV to have a fast preview of geotiff files. I tried to install it in Ubuntu 16.04, but it isn't in the ubuntu-gis ppa nor could I find another simple way to install it. The software seems to be abandoned. I know I can also pipe raster files into qgis, but qgis is often too heavy for a quick look and other Image viewers (e.g. Eye of Gnome) often struggle with geotiff files.
So i was wondering if there is a simple FLOSS alternative for openEV that can be used in newer versions of Ubuntu and can be used in the command line.


texts: Converting NOAA's World Ocean Database NetCDF files to tabular data using Python

texts:
I'm trying to convert variables from NetCDF files to tabular data, but can't understand how to group observations by some index or key variable.

For example, when loading data from WOD queries using the xarray library

import xarray as xr

ds = xr.open_dataset(r'ocldb1573479266.24760_OSD.nc')
ds

I get something like

<xarray.Dataset>
Dimensions:                         (Alkalinity_obs: 32749, Ammonia_obs: 15252, CFC11_obs: 16256, CFC12_obs: 16317, Chlorophyll_obs: 58045, DeltaC14_obs: 1456, DeltaHe3_obs: 1929, Helium_obs: 2143, Neon_obs: 377, Nitrate_obs: 108954, Oxy18_obs: 409, Oxygen_obs: 300929, Phosphate_obs: 183091, Pressure_obs: 581435, Salinity_obs: 581435, Silicate_obs: 167977, Temperature_obs: 581435, Tritium_obs: 2015, biosets: 1149, casts: 51401, numberofpis: 6079, pCO2_obs: 8513, pH_obs: 42365, tCO2_obs: 24672, z_obs: 581435)
Coordinates:
    lat                             (casts) float32 ...
    lon                             (casts) float32 ...
    time                            (casts) datetime64[ns] 1890-06-27T02:00:00.000214528 ... 1988-03-05
    z                               (z_obs) float32 ...
Dimensions without coordinates: Alkalinity_obs, Ammonia_obs, CFC11_obs, CFC12_obs, Chlorophyll_obs, DeltaC14_obs, DeltaHe3_obs, Helium_obs, Neon_obs, Nitrate_obs, Oxy18_obs, Oxygen_obs, Phosphate_obs, Pressure_obs, Salinity_obs, Silicate_obs, Temperature_obs, Tritium_obs, biosets, casts, numberofpis, pCO2_obs, pH_obs, tCO2_obs, z_obs
Data variables:
    country                         (casts) |S40 ...
    WOD_cruise_identifier           (casts) |S40 ...

(. . .)

    Temperature                     (Temperature_obs) float32 ...
    Temperature_sigfigs             (Temperature_obs) int8 ...
    Temperature_row_size            (casts) float64 7.0 6.0 7.0 ... 58.0 45.0
    Temperature_WODflag             (Temperature_obs) int8 ...
    Temperature_origflag            (Temperature_obs) float32 ...
    Temperature_WODprofileflag      (casts) int8 ...
    Temperature_Access_no           (casts) float64 ...
    Temperature_Scale               (casts) |S170 ...
    Temperature_Instrument          (casts) |S170 b'' b'' b'' ... b'' b''

(. . .)

When inspecting shapes

print('CASTS          : ', ds.casts.to_dataframe().shape )
print('Temperature    : ', ds.Temperature.to_dataframe().shape )
print('Temperature_obs: ', ds.Temperature_obs.to_dataframe().shape )
print('time           : ', ds.time.to_dataframe().shape )

I get the following

CASTS          :  (51401, 4)
Temperature    :  (581435, 1)
Temperature_obs:  (581435, 1)
time           :  (51401, 4)

So, as expected, it seems for each cast there is a set of Temperature observations. But I'm not familiar with xarray (or NetCDF4-python) methods for mask, group, merge, etc. So my question is: how can I create something like this?

|     | cast | lat | lon | temperature | ... |
|:---:|:----:|:---:|:---:|:-----------:|:---:|
|  0  |   0  |     |     |             |     |
|  1  |   0  |     |     |             |     |
|  2  |   0  |     |     |             |     |
| ... |      |     |     |             |     |
|  21 |   0  |     |     |             |     |
|  22 |   1  |     |     |             |     |
|  23 |   1  |     |     |             |     |

I already tried to use to_dataframe() method, but it did not work.


texts: Trying to change the coordinate system of feature classes

texts:
I'm trying to change the coordinate system of a batch of feature classes to 'Geodetic datum GDA94(code=4283) using Define Projection. I receive no python error message and no change to the coordinate system, just lists returned from ListFeatureClasses and ListWorkspaces. My code is below. Can someone suggest why Define Projection isn't converting the data it receives?

import arcpy
from arcpy import env
import os

#set path where the shapefiles are located
arcpy.env.workspace = "Z:/temp"

#returns @string list of GDB folders
workspaces = arcpy.ListWorkspaces('*', "FileGDB")
for item in workspaces:
    print item
    env.workspace = item

#returns @string list of feature classes
    fcs = arcpy.ListFeatureClasses()
    for file in fcs:
        print '\t', file

 # set local variables
        for file in fcs:
          inData = file
          coordinateSystem = arcpy.SpatialReference(4283)
          arcpy.DefineProjection_management(inData, coordinateSystem)


texts: Getting SRID Well Known Text (WKT) metadata from PostgreSQL?

texts:
The SRID metadata is showed by, e.g., this web page:

 Definition: Well Known Text (WKT)
 PROJCS["UTM Zone 24, Northern Hemisphere",
GEOGCS["GRS 1980(IUGG, 1980)",
DATUM["unknown",
SPHEROID["GRS80",6378137,298.257222101],
TOWGS84[0,0,0,0,0,0,0]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",-39],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["Meter",1],
AUTHORITY["epsg","6211"]]

This information is on my local database, but how can I see it?
A function like SRID_AsText(srid). Or some complex PostgreSQL query over all PostGIS tables to see the SRID metadata?


texts: Unable to Drag to rearrange layer drawing order in Table of Contents of ArcMap?

texts:
I'm new to working in ArcMap, but I seem to have a very small issue.

I'm using ArcMap 10.2.2 on Windows 7 Ultimate and for some reason, I can't rearrange the layers I'm working with in the table of contents by clicking and dragging them in the order I want.

Is there some sort of option I have to enable?


texts: Connect QGIS to the custom application server

texts:
I'm new to QGIS and investigating how to migrate one home made GIS application on QGIS.

The biggest difficulty for me now is to find the way how to plug QGIS to the application server that is being used now to fetch and save data.

What is the right direction to begin with?