texts: Is there an industry-standard official mapping of GLONASS satellites to "PRN" values?

texts:
In the U.S. GPS system, satellites are uniquely identified by the PRN number broadcast by each satellite.  This ID is usually visible on most GPS receiver platforms that provide access to individual GPS satellite signal information to identify the broadcasting satellite (e.g., accessible on Android via the GPSStatus.Listener and the GpsSatellite.getPrn() method (http://goo.gl/QLse3)).

However, the Russian GLONASS system does not currently have the concept of a PRN number broadcast by every satellite, as many of the legacy satellites are an FDMA system, not a CDMA system such as the U.S.  So, in GLONASS these satellites broadcast the same code on a different frequency.

However, an ID is still available for GLONASS satellites for many GPS receiver positioning platforms that support dual GPS/GLONASS positioning, often via the same access point as the PRN value for U.S. GPS satellites (e.g., on Android via the GpsSatellite.getPrn() method, with seemlingly standardized values across devices I've seen from 65 to 88.  See the first screenshot listed in the GPSTest app on Google Play (https://play.google.com/store/apps/details?id=com.android.gpstest) for an example list of U.S. and GLONASS satellites obtained from the GPSStatus.Listener on a Samsung Galaxy S3).

I haven't yet been able to find an official industry-standard source for how these ID values for GLONASS satellites are determined.  However, proprietary documentation for particular products have shown GLONASS satellites numbered from 65 to 88 [1] and 38 to 61 [2].

Is there any industry-standard official mapping between GLONASS satellites and "PRN" values?


texts: How to get value for specific coordinates in R

texts:
I have a bunch of raster files in different formats:
- NETCDF 4
- HDF 5
- GeoTIFF

Using rgal package I'm able to open any of them and use them as objects in R!
Creating stacks of them allows me to calculate correlation coefficient.

But what if I would like to mąkę some calculations on one pixel only (for example to create scatterplot)?

Am I able to extract value for particular coordinates?
Obviously I can use gdal2xyz tool, but I would like to do that entirely in R!


texts: Error when running BetterBusBuffers "Count Trips on Lines" tool in ArcGIS Pro

texts:
After running the "Count Trips on Lines" tool of the BetterBusBuffers toolbox, I encounter the following error after running step 2, "count trips on lines:"
"There is no transit service during this time window. No service_ids cover the weekday or specific date you have selected."
On further inspection of the GTFS package's calendar.txt file, I find that every one of the [day of week] fields is entered as zero, which seems unusual. Original GTFS package is here.
Screenshot is below:

Any ideas for a workaround?


texts: Does GDAL has any substitute function for the ExtractByMask in arcpy

texts:
I am now using arcpy: arcpy.sa.ExtractByMask(weight_raster,elev)
both weight_raster and elev are raster data in different size.

I want to use functions in GDAL to replace this arcpy methods. Does anyone know how to do it?


texts: Water features in ArcScene viewer not displaying correctly?

texts:
I am using ArcMap 10.4.1 and ArcScene 10.4.1.

My scene consists of dtm, multipatches, polygon features and point features. Together they work quite well, except for the water body (polygon) showing where it should not. You can see the outline of the island in the first picture (but with water where it shouldnt be), and when I zoom in, the island appears. In the third image, you can see how the polygon looks in ArcScene.

It seems like there is some problems with the polygon having a hole inside it, and that something isnt rendering correctly or something. I have tried extruding it on top of the terrain and setting a fixed height, but there seem to be something not quite right.

I've also tried setting rendering hierachy, and interpolating the polygon to match the terrain height, but with no success.

Could there be something with the terrain model that is not working 100%? Is there a best practice when combining these types of geodata into ArcScene Viewer?


texts: Selecting a Polygon based on the location of a point

texts:
Currently I have a polygon layer made up of three shapes. All the shapes are listed in the attribute table. One of these three shapes intersects a point. I would like to select just the shape that has an intersection of that point.

I am guessing I use the select tool with a piece of code to select the correct shape?

Thanks for reading

Daniel


texts:
attention: This might be more a question for "askubuntu.stackexchange" but it is GIS-related
Geonode 2.0 is available for Ubuntu 12.04 but not for Ubuntu 14.04. According to the docs it should be sufficient to make an

sudo apt-get install python-software-properties
sudo add-apt-repository ppa:geonode/release
sudo apt-get update
sudo apt-get install geonode

I have a fresh install of Ubuntu Server 12.04 with a default LAMP configuration and it already fails on the first line:

Depends: python-pycurl but is not installable

The first line works on Ubuntu 14.04 but geonode will not be installable on 14.04. I am somewhere stucked. Is there a Geonode release available for 14.04?


texts: Making combo chart in Google Earth Engine

texts:
How can I make combo chart in Earth Engine? I have a code that makes a chart based on image collection for given polygon. I would like to draw the mean value line in the chart as well. Is it possible to do it in Earth Engine?
var geometry = ee.Geometry.Polygon(
        [[[11.396256200031866, 51.10045908804461],
          [11.396256200031866, 51.06940343482601],
          [11.451187840656866, 51.06940343482601],
          [11.451187840656866, 51.10045908804461]]]);

var collection = ee.ImageCollection("MODIS/006/MOD13Q1")
      .filterDate(ee.Date.fromYMD(2019, 1, 1), ee.Date.fromYMD(2019, 12, 31))
        .select("NDVI");

var chart = ui.Chart.image.series(collection, geometry, ee.Reducer.mean())
      .setChartType('ColumnChart');

print(chart);

What I have now:

What I want to have:


texts: Deegree Configuration of existing tile cache (.png files) as datasource for wms and wmts

texts:
I am unsuccessful in configuring an existing png-file tile cache that I want to use as data source for the provision of wms and wmts services.

Despite having worked myself through the deegree documentation and the utah example data set (that uses a geotiff-tile cache) I did apparently not configure deegree correctly.

I have a tile cache having the structure zz/xx/yy.png
with a total amount of 9 zoom levels.

Within the datasource/tile/tilematrixset directory I have defined this tile cache as follows (example for two zoom levels):

<TileMatrixSet xmlns="http://www.deegree.org/datasource/tile/tilematrixset" configVersion="3.2.0">

  <CRS>urn:ogc:def:crs:EPSG::31467</CRS>

  <TileMatrix>
    <Identifier>0</Identifier>
    <ScaleDenominator>1111111</ScaleDenominator>
    <TopLeftCorner>427461.787361 103142.057325</TopLeftCorner>
    <TileWidth>256</TileWidth>
    <TileHeight>256</TileHeight>
    <MatrixWidth>1024</MatrixWidth>
    <MatrixHeight>1024</MatrixHeight>
  </TileMatrix>
  <TileMatrix>
    <Identifier>01</Identifier>
    <ScaleDenominator>555555</ScaleDenominator>
    <TopLeftCorner>427461.787361 103142.057325</TopLeftCorner>
    <TileWidth>256</TileWidth>
    <TileHeight>256</TileHeight>
    <MatrixWidth>512</MatrixWidth>
    <MatrixHeight>512</MatrixHeight>
  </TileMatrix>
</TileMatrixSet>

I do not know how to set up the required xml document in the xml document within the datasource/tile/ directory since I do not know if my cache qualifies as a filesystem cache.


texts: Is there an iOS app to display several offline raster maps with transparency?

texts:
There are several iPhone apps displaying offline maps.
Some of them can show custom maps, for instance loaded as MBTiles or sqlite files.

Galileo-app does this pretty well.

But it is only possible to display only one tiled map at a time.
I need to display several tiled raster layers (say 3 or 4) with a transparency (for instance with an opacity control slider for each layer).

My rasters are hand drawn tracks (black 'wire' on white/transparent background), old maps tiles and openstreetmap tiles. Goal: find my way by crossing information from those georeferenced maps.

TiledMap app looks to be doing this, but is missing from iTune (dev account not renewed).

Is there any iOS app doing this? (Otherwise I would have to code one, but this is too much work)

Thanks


texts: ArcMap 10.7, Is it possible to mosaic rasters using a method that prioritizes rasters by the proximity to it OWN center?

texts:
I am working with older imagery that has been manually scanned, and as expected, the fidelity of the images is more reliable towards the centers of the images. There is a good deal of overlap in the images that I need to mosaic, and it would be really helpful if I could mosaic them based on which raster's portion of the overlapping area is closer to its own center. Clearly this is not what the "closest to center" method does.
Another thing that could work is if I clipped out a portion of each image, leaving behind only a central area where the imagery is more reliable. However the images have already been georeferenced and this would involve creating a huge number of separate masks or rectangles. Since all of these images are about the same size, it would be ideal if there were a tool like "extract by rectangle," but where you could define the rectangle based on its position within the raster, rather than that of the overall data frame.


texts: How to filter for distinct values and apply multiple filters in QGIS attribute table?

texts: I am very new to QGIS. I am doing an assignment for class where we have to answer several questions about data obtained from Global Fishing Watch. One of the questions is how many vessels have a tonnage greater than 1200? I know how to filter for this. However, the data includes multiple voyages for each ship. As you can see in the attached screenshot, when I filter for tonnage, many of the ships appear more than once (the MMSI field is the ID number of the ships). My question is how do I filter for tonnage AND make it so that each ship will only appear once? I figure I need to filter by tonnage and filter the MMSI for distinct values, but I am unsure how to apply two filters to a search and I need to know how I can filter a field so that it will remove duplicates. I have a feeling that the "count_distinct" function is what I need, but so far it hasn't worked, I think my syntax is incorrect. As for applying two filters, just putting "and" between my expressions has not worked thus far. As I said I just started and do not yet have a hang of how the expression syntax works.


texts: How to obtain statistics from an intersection?

texts:
I am new to the earth engine. My question is on how to obtain statistics from the overlapping region of two areas.

I have uploaded a table that contains the locations of river basins (https://code.earthengine.google.com/?asset=users/basins/AqueductBasins). These commonly lie in several countries. I want to obtain the information on the mean population density for these basin areas in a specific country (so: mean for Basin1 in CountryA, mean for Basin1 in CountryB, mean in Basin2 in CountryA and so on).

I have taken this feature collection on the country borders:

var countries = ee.FeatureCollection('ft:1tdSwUL7MVpOauSgRzqVTOwdfy17KDbw-1d9omPw');

This image contains information on the people per cell:

var Pop2015 = ee.Image('JRC/GHSL/P2016/POP_GPW_GLOBE_V1/2015');

I want to print a table with the average people per cell in the basin area in the different riparian states.


texts: Creating a spatial PostGIS table as a result of a PostGIS query

texts:
When I want to store the results of a spatial query as another PostGIS table, what should I do? I am interested in identifying all residential buildings that fall within a set of parcels:

CREATE TABLE residentialBuildings AS
SELECT buildings.the_geom
FROM parcels as p
JOIN buildings as b
ON ST_Contains(p.the_geom, b.the_geom)
WHERE p.residentialarea > 0;

This seems to work, but it does not return a spatial result - instead the table has the binary representation of the geometry. I'm not sure what I should be doing to create a table as a result. Do I need to first make a table, and then convert the result into a geometry?

(I have read the documentation, but I feel like I'm missing the right search terms as I feel this must be a fairly obvious thing to do).

I realize that there are many commands to make a geometry

ST_GeogFromText
ST_GeographyFromText
ST_GeogFromWKB
ST_GeomCollFromText
ST_GeomFromEWKB
ST_GeomFromEWKT
ST_GeometryFromText
ST_GeomFromGML
ST_GeomFromKML
ST_GMLToSQL
ST_GeomFromText
ST_GeomFromWKB

but I don't know which option to choose from, or if this is the correct approach.


texts: QGIS DB Manager Import layer/file Oracle not working - Error 8 Loading of the layer failed

texts:
QGIS 2.12.0 - DB Manager 0.1.20 - Windows 7 - Oracle 11g 11.2.0.3.0

Import vector layer

Shapefile as well as a layer (vector) and with every possible combinations of settings, the results are the same. It also happens with a previously exported table (as shp) from same db/schema:

Error 8
Loading of the layer "myschema"."mytable" failed

An inspection of the "loaded" table indicates that only the 'geom' column was created (no other columns where created), it is however of UNKNOWN geometry type and the table is empty.

user_sdo_geom_metadata record was created correctly.

All other operations within DB Manager on existing tables work correctly.

Is there perhaps another way of loading existing vector layers onto Oracle (pyQGIS?)?


texts: I have 2 tables. Each table has an ID and i need to insert only one attribute from first table to a specific field in the second table using python script. This should be done using UpdateCursor, but i can't figure out how, Any tips?


texts:
Is there any way to have the legend text above the symbol and beside it, not just above?

What i do at present is to cut the end of the text off and insert a label then group them, however i want to resize the box and possibly have the text go into two lines, or preferably wrap around the symbol.

Excuse the sloppy work below, a rush job for illustrative purposes.


texts: Geocortex workflow 5 calculate service area

texts:
How can I use calculate service area activity in Geocortex workflow 5?
I have tried all the possible ways using ArcGIS NA Server service layer. When I am running the workflow error is displayed: invalid workflow inputs. I am unable to understand what is the issue.


texts:
I'm trying to query the database (DataStore) and want to let the datastore calculate a value for me.

In SQL I would do it like this:

SELECT name, ST_Area(ST_Intersection(the_geom, ....) AS myValue FROM ...

Is there any way, how I can do this with CQL. My current approach is this, but I do not know how to add this calculated field. I need this, as I want to sortBy this field.

DataStore ds = ...
Filter filter = CQL.toFilter("INCLUDE");
SimpleFeatureCollection features = source.getFeatures(filter);

Any idea, how to add my area function to the resultset?


texts: PostGIS - For each polygon in table A, find nearest point in table B and update table A with value

texts:
I have a table called sites of polygons, and a table called postcodes of points. For each row in sites I need to find the nearest point in postcodes, then update the sites column postcode from column pc in postcodes. Can anyone please help me write this in PostGIS?

This gets all the correct values, I just don't know how to do an update that can read these straight into the postcodes column:

SELECT DISTINCT ON(g1.key_column)  g1.key_column As key_column, g2.pc As
postcode
FROM sites As g1, postcodes As g2
WHERE g1.key_column <> g2.key_column AND ST_DWithin(g1.geom, g2.geom,
1000)
ORDER BY g1.key_column, ST_Distance(g1.geom, g2.geom)


texts: Selecting Specific Entities to Style Separately from the Same Layer in QGIS 3.2.4

texts:
I downloaded the roads of Rome from OSM and save them as a single Shapefile.

To style specific entities from the "highway" column I categorized this column in the layer styling panel.

Now I got all the 28 highway types, each in a different color.

However, I want to change the style for only one entities and to keep all the other in a different style.

To do it manually will take a lot of time, so I am looking for a faster way.


texts: How can I query for a list of specific road types in the same Overpass query?

texts:
I'm trying to query for a number of specific road types in the same request, and I can't get the syntax right at all.

Can anyone alter the following query to select "motorway", "trunk" and "primary" roads?

[timeout:25];
way["highway"="motorway"](53.5084,-2.9916,53.6,-2.8);
out body;
>;
out skel qt;


texts: Flip over min/max values in a stretched symbology for raster data

texts:
In ArcMap 10.1, we can render a raster in a colour ramp (stretched) symbology to draw a raster layer based on each pixel value. However, the legend for the layer by default starts from the max value (top, red) down to the min value (bottom, green).

How can I draw it like form min (top, green) to max (bottom, red) instead?

Tried Invert checkbox but it didn't help as it was the colour that was inverted.


texts: How to go about creating a custom tool/adding special functionality for a web app

texts: In my soon to be developed web app, I will have an image service and also a feature service that is essentially an index to the image service, with all kinds of attributes. I would like for the user to be able to select a feature in the index layer, and then get a prompt to view the corresponding image in the image service. I only want the user to see the images that they have 'selected' by selecting features in the index layer. I know python pretty well but am still new with javascript. I was wondering about creating a custom tool in python and then publishing that as a GP service, but I'm unsure of how to manipulate feature/map services in python. I would assume this could be accomplished in javascript but I dont know enough to know how to do that. In terms of the web app, I would try to build it either in arcgis online or leaflet (which I am just learning). Wondering if anyone has some advice on a good approach to this task? I am going to start by investigating how to use the ArcGIS REST API in python.


texts: Dissolving redundant fields in ArcGIS for Desktop?

texts:
I have a map of a small neighbourhood. I'm trying to label these but because of the grid being squares the objects are split and therefore showing two two or more labels of each location.

I'm trying to use the dissolve tool to fix this but I must be doing it wrong.

First I open dissolve from geo-processing.

I input "Administrative text" as input features (the layer containing the names).

In the dissolve fields I've tried numerous ways but nothing has worked.

The attribute containing locations names is "TEXT".


texts: QGIS - Buffer distance as portion of polygon variable width

texts:
Is it possible to create a buffer of for example 3/4 of the polygons local width?
I have a polygon of river body, somewhere the river has 20 meters in width, somewhere 70 m. I need to extract the cunette, which varies based on the rivers width. So i need to create a negative buffer that is 3/4 of the current river body polygon width, so it copies its current geometry.
I think this must be feasible, however i can't find a way...
Any ideas?


texts:
I am getting an error trying to create a bounding box of a geometry using ST_Box2D(). Is the function deprecated? If so, what is the equivalent function?

SELECT ST_Box2D(ST_GeomFromText('LINESTRING(1 2, 3 4, 5 6)'));


texts: Modifying available attributes when reading CAD feature classes using ArcObjects

texts:
I would like to read CAD data including the default attributes provided by direct-read CAD data models, as well as other 'optional' attributes, using ArcObjects.  The direct-read CAD data model can be found here:
Direct-read CAD Data Model
Basically, when reading in CAD feature classes, I have the following classes, with predefined attributes:

Annotation
Multipatch
Point
Polygon
Polyline

I've read up on class extensions and plug-in data sources, but I wasn't sure which, if either, would provide what I am seeking.
Is there a way to modify/extend these pre-defined classes, or define a new class, to allow me to read additional attributes not included in the direct-read model?


texts: QGIS 2.2 is very slow while saving the edited layers

texts:
I am using QGIS 2.2 64 bit in windows 7 system. I have an issue while saving the layers after editing. Its taking more time to save my edits. Some times it's working fine but very often i am getting this issue and also the QGIS will be in  Not Responding stage at that time. I tried in different postgres database and different types of layers(polyline and polygons) but still having the same issue.

Is anyone having any idea to resolve this issue?


texts:
My question is very similar to Cost Distance Analyst In QGIS?

However, I'm looking for a means by which I can create a raster displaying the cumulative cost of each pixel to a specified destination rather than simply the least cost path between sets of specified start & end points. In essence, creating a "friction" map in which each raster pixel's value is it's lowest "cost" to get from itself to a specified destination (polygon or single pixel).

I imagine there's a way to do this in QGIS using GRASS operations, but can't figure out how to produce such a map -- does anyone know how? I would imagine it to be the r.cost() function, but without the need to specify start points and stop points.


texts: Add an extended set of plugins to a fresh QGIS installation

texts:
There are a number of non-default QGIS plugins I use every day that are not included in a fresh installation of QGIS. QuickMapServices, QConsolidate are just to name a few - I think everyone has their own list. Unfortunatelly, every time I install QGIS to a new machine I have to manually download every item on this list. So I wonder if there is a way to add reference to such plugins to a fresh QGIS installation.

The best think I came up with so far is to restore my list with a backup of "plugins" directory.


texts: How to get the polygon vs. Polygon Intersection Coordinates using python

texts:
I want to find the intersection coordinates between two polygons.

aLayer = canvas.currentLayer()
selection = aLayer.selectedFeatures()
for geom1,geom2 in itertools.permutations(selection,r=2):
    intersect=geom1.geometry().intersection(geom2.geometry())
    intersect_list=intersect.asPolyline()
    count=len(intersect_list)

Its giving the intersection geometry but when i am converting that as a polyline to get the intersection vertices it is giving empty list. Here what i am doing wrong.


texts: how to show all features visible on map?

texts: i have restored on mssql a backup file and added it on qgis but it shows only some features NOT all features. i want all features to be shown on map.


texts: Using Python-Flask-MongoDB with Geolocation field running on App Engine?

texts:
I'm trying to write an app which would act as a address caching API, so basically it needs to store addresses (geolocations) and later on I could query on a random address and radius and get all the cached locations within the radius of that point.

I prefer to use Python on Google Appengine, I tried DataStore GeoPT but to do such a query it would end up being really expensive and not really efficient.
So I started reading about MongoDB Geospatial Query, $near option. But apparently there is no support for any Python-MongoDB-GAE that also supports GeoLocation fields e.g Flask-MongoEngine.

Does anyone know a solution to this or even another stack to work on which makes this implementation feasible?


texts: Performance of 'extract by mask' tool in PostGIS?

texts:
Context: A large raster is clipped to create a bunch of smaller rasters based on a bunch of shapefiles.

I've used st_clip() which seems to take forever looping through the different shapefiles, so I've switched to the much faster (but fast enough?) combination of  'where st_intersects()' and smaller tile sizes.

The Extract by Mask tool on arcgis seems to be faster.
Is there a difference in how they are both done?

As I understand it, it rasterizes the vector shapefile and creates a matrix of 0s and 1s and extracts the values from the base raster which intersect with the 1s.


texts: Returning centroid coordinates from footprint coordinates using QGIS

texts:
I have a set of footprint coordinates in KML. However I do not have the building coordinates.

How do I return building's centroid (single point) by using the footprint coordinates in KML using QGIS?


texts: Why does Leaflet with provider tiles not draw markers as expected?

texts:
I'm working with Leaflet maps in Shiny for R. I'm producing a map with some markers on it. When I use addTiles () for my base map, everything works as expected. However with addProviderTiles (), the markers don't appear on the map. Only after I trigger an action that causes them to be redrawn, they show up as I would expect them.

This doesn't seem to be mentioned in the vignettes, so I'm not sure if it's a bug or I'm missing some parameter options.

I made a reproducible example. If you run it as is, the markers will appear normally. If you comment line 13 and uncomment line 14, the markers won't work initially. I added an observer () that redraws the markers when you click on the map.

library(shiny)
library(sp)
library(leaflet)

ui <- bootstrapPage (
                leafletOutput("map", width="80%", height="800px"))

server <- function(input, output, session) {
  output$map <- renderLeaflet({
      print ("render tiles")
      b <- bbox (breweries91)
      leaflet (breweries91) %>%
      addTiles () %>%
#      addProviderTiles (providers$CartoDB.Positron) %>%
      fitBounds (b[1], b[2], b[3], b[4])
      })

  observe ({
      print ("make markers")
      leafletProxy ("map", data = breweries91) %>%
      addAwesomeMarkers ()
      cl <- input$map_click
      if (!is.null (cl))
          print ("click")
      })
}
shinyApp (ui, server)

Edit software versions:

Ubuntu: 16.04.2
Mozilla Firefox 52.0.1
R: 3.3.3
leaflet: 1.1.0
shiny: 1.0.0


texts: CartoDB Torque viz: make bubbles grow to show intensity/frequency?

texts:
I'm working on a project that maps the places of origin of the "letters to the editor" submitted to a magazine.

I'm using CartoDB's Torque viz to show the evolution of this data through time.

I would like to show that multiple letters originated from the same address by increasing the size of the bubble at that address. For example, the Marker Fill is 4 for one letter from a particular address, 6 for two letters, 8 for three letters, etc.

Does anyone have any suggestions about how I might accomplish this?

Many thanks in advance!


texts: Why does Q-Chainage sometimes offer meters and sometimes degrees?

texts:
Q chainage is a plugin which puts regular points on a line vector at regular intervals which you can specify.

For a shapefile I've been given (and converted to line vector), it offers me to specify interval in meters.

For a shapefile I've created (and converted to line vector), it offers me degrees.

Degrees are no use to me. I want to use MY created line vector, and be offered meters by the plugin.


texts: Creating relationship class with ArcPy using GlobalID?

texts:
I'm trying to create a relationship class with Python using a GlobalID.

When I use Python it gives me the error:

  Failed to execute. Parameters are not valid. ERROR 000800: The value
  is not a member of GLOBALID. Failed to execute
  (CreateRelationshipClass).

When I try to use the Create Relationship Class tool through the Toolbox, I get the same error. If I select GlobalID as the primary key, it will ONLY let me select a GlobalID as the foreign key. It will not let me select a GUID of my choice.

But when I'm in ArcCatalog and right click the Dataset and go to New > Relationship Class...

It works!

What's wrong here?


texts: OpenLayers: failed to create a Ext.list that wraps a Openlayers.request.Get function

texts:
I am using Sencha Touch + OpenLayers to develop a demo app. One feature is to get the catalog of feature data in GeoServer. I want to use Ext.List to display the name of feature layers. I tried following codes:

App.AddFeatureLayerList = Ext.extend(Ext.List, {
createStore: function(){
    Ext.regModel('WFSLayers', {
        fields: ['name', 'title', 'srs', 'featureNS']
    });
    var data = [];
    var request = OpenLayers.Request.GET({
        url: root+'/geoserver/ows?service=wfs&version=1.0.0&request=GetCapabilities',
        handler: function(response){
            var XMLformat = new OpenLayers.Format.XML();
            var xml = XMLformat.read(response.responseText);
            var CAPformat = new OpenLayers.Format.WFSCapabilities();
            cap = CAPformat.read(xml);

            for (var i = 0; i < cap.featureTypeList.featureTypes.length; i++) {
                var featureT = cap.featureTypeList.featureTypes[i];
                data.push({
                    name: featureT.name,
                    title: featureT.title,
                    srs: featureT.srs,
                    featureNS: featureT.featureNS
                });
            }
            return new Ext.data.Store({
                model: 'WFSLayers',
                sorters: 'title',
                data: data,
                getGroupString: function(record){
                    return record.get('title')[0];
                },
            });
        }
    });
},
initComponent: function(){
    this.store = this.createStore();
    this.itemTpl = new Ext.XTemplate('{name}');
    this.grouped = true;
    this.listeners = {
        itemtap: function(dataview, index, item, e){

        }
    };
    App.AddFeatureLayerList.superclass.initComponent.call(this);
}
});
Ext.reg('app_addFeatureLayerList', App.AddFeatureLayerList);

However, i got an error:

Uncaught DataView requires tpl, store and itemSelector configurations to be defined.

then if I remove the codes:

new Ext.data.Store({
            model: 'WFSLayers',
            sorters: 'title',
            data: data,
            getGroupString: function(record){
                return record.get('title')[0];
            },
        });

and put them to pass to this.store after call "this.createStore()", I will get a empty list. however, I traced array data and it turns out it has data, but the store is empty.

another funny thing is if I put 'createStore()' as a public function, and call it in list configs, e.g, 'store: createStore()', I will get an error message saying my IP is not allowed to access to target URL, which this is not correct since I use a proxy and I added other WFS layers from this server already...
I can't find what's wrong with it. Anybody could help me give me some hints?

Appreciate all help!!!


texts: How to ignore a zero in denominator when computing for mean of rasters?

texts: Supposed I have images A, B and C. They have pixels a, b and c respectively in index (0,0). I want to compute for the mean in index (0,0). Simply (a+b+c)/3. However I want to omit pixel that have nodata or less than or equal to 0. For example if I had c=0 then, (a+b+c)/2 will be the average for that index. I want to know how to do this in raster calculator in QGIS.


texts: Set background or text color depending on expression in an attribute table in QGIS composer

texts:
I am using QGIS 3.0.2.
I am working on a print composer including an attribute table.
I would like to set different colors to each rows, depending on an expression, refering to an attribute.
I didn't find a way.

Is that possible?


texts: How do I edit multiple features at one time?

texts: In my case, there are about 30 overlaid features. How do I select all and edit them at one time? I want to reshape the features, mainly using the node tool to reshape a straight line into a (fake) curve one.


texts: Comparing attributes of two shapefiles in ArcMap?

texts:
I have two identical shapefiles of Estonia. The only difference is that the attribute table of shapefile 1 has mean temperatures from the period of 1991-2000 and the attribute table of shapefile 2 has mean temperatures from the period of 2001-2010.

I want to compare these two shapefiles now and the result should be a shapefile which shows the difference in mean temperature between the two periods.

How can I do it?

I tried Feature Compare but it did not work. I got an error.


texts: Leaflet.js map - Pan too and Zoom in on point

texts:
This is the first time I am using Leaflet.js. I am trying to create a story map but have come across an issue.

Currently I have a map with a scrolling text bar on the left hand side. I would like users to be able to scroll through the text on the left and when they click on a part of the text it will automatically pan too and zoom in on the correct area on the map, at the same time displaying the marker.


texts: How to extract polygon data from Google maps?

texts:
Is there any way to extract polygon data from Google maps into kmz or shp?

Here the website: http://www.vodafone.co.nz/network/coverage/

note: I only want 3G coverage of the country


texts: Achieving more than basic labelling of Reference Grids (Graticules, Measured Grids) in ArcMap

texts:
I'm happy to be corrected, but at the moment I do not think it is possible to meet any of the following mapping requirements using either ArcMap or ArcPy:

Label measured grids using a format like 6828000mN where 68 is superscript, 28 is normal, 000m is superscript and N is normal
Label measured grids using a format that at least enables mN and mE to be included even if all the superscripting above cannot be supported
Label graticules using a Degrees, Minutes and Seconds (DMS) format where the minutes and seconds can be written as 00 rather than 0

Does anyone have any thoughts on how to achieve one or more of the above?

Note that I have excluded Grid and Graticule layers from consideration here, because even if they can support these grid and graticule labelling scenarios they are not dynamic and I need the labelling above in a Maps On Demand (MOD) / Data Driven Pages (DDP) like application.

As an update it is possible to do some of 1. and 2. using Corner Labels in a static map as long as it is displayed in a Projected Coordinate System.  However, this hangs ArcMap (and thus prevents it being saved) when the Coordinate System of the Data Frame is switched to Geographic which is necessary for an application like mine that uses a Geographic Coordinate System to determine the appropriate MGA Zone to use for each map at runtime.  Basic labels work for Projected Coordinate Systems just not the fancy ones required in 1.


texts:
I am an economics student working on a paper for which I need to find and analyze regional data for SO2 emissions.
I came across some .1x1 data files but I don't have any idea how to open/display these files.


texts: Batch geocoding 1.5 million data points in SQL Server

texts:
I am currently attempting to batch geocode 1.5 million points that are stored in a sql database.  The points are not very specific (ie only city, state points).  I have used ArcGIS in the past but due to the new ArcGIS online credit system, it is restrictively expensive to geocode with their database.  For that reason, I was able to access a composite geocoder that I attempted to connect to ArcGIS via the Create a geocoder feature, however, this system is not working as well (either it freezes with my database, or if I use a test dataset, it cannot match any addresses).

I am seeking suggestions on:

How to best setup the geocoder I have to run in ArcGIS or
Other methods that may work better (I am really open to anything!)

Let me know if there is more helpful information that I can provide.


texts:
I have a field that is a Long type which represents a date.  It is in YYYYMMDD format.  (Example 04/10/2018 is 20180410.)

How do I convert that Long to a Date type using ArcGIS Pro 2.1?


texts: Adding widget resizing feature to Configurable Map Viewer?

texts:
Is there any way of adding a widget resizing feature or implementing it?

I noticed that the layer widget for example expands according to the number of layers inside it, this can be used to implement widget resizing.


texts: ArcMap creates unnecessary file types when saving shapefile (.cpg .prj, lock, etc.)

texts:
When saving shapefiles, geocoding results or even exporting data in ArcMap, many different files are created. Most of the time when I save a shapefile, I want only the shapefile and not a .cpg, .prj, lock file, etc.

Is there a way to tell ArcMap (desktop version) using either the UI or Python code, to only save one file when saving shapefiles, dBase tables, etc.


texts: different results :Slope analysis in qgis vs ArcMap

texts: I want to make slope analysis on geoTiff. I am using ArcMap and qgis. But i get different results, but both are in degrees. I ArcMap i recive range from 0 to 42 degrees, but i qgis it from 0 to 17 degrees. I ArcMap i use slope function in spatial analysis. Z factor 1.  In qgis i use slope i raster analysis, z factor 1.


texts: Sampling Image collection google earth engine

texts:
Question: How can I get a random sample of images from an image collection?

Background:
I am working with the GRIDMET: University of Idaho Gridded Surface Meteorological Dataset which provides daily weather for the continental US from 1979-2015 as an image collection. I am interested in calculating the 90th percentile value for one of the weather variables (specifically ERC) at a large number (30,000) points in Washington and Oregon. The GRIDMET image Collection is daily from 1979-present about 13,000 images, obviously calculating the 90th percentile from all the images is taking a long time. I have already clipped the image collection to only the points of interest. I am looking to randomly sample the image collection to speed up my code.

Current Code:

// import Fire Features and GRIDMET image collection
var fires = ee.FeatureCollection("users/brycekellogg/Fires_test");
var Idaho = ee.ImageCollection("IDAHO_EPSCOR/GRIDMET");

Map.setCenter(-119.905, 45.836, 6); // Center on Oregon.
Map.addLayer(fires);// add the 5 test fire points to the map.

var firefeature = ee.FeatureCollection(fires); // convert test fires into a
feature collection

var fullERC = Idaho.select('erc')// select only the ERC values

// Create a funtion to clip each image to the feature collection
var clipFeature = function(image) {
  return image.clipToCollection(firefeature);
};

// Map the clip funtion over the Image Collection
var fullERCclip= fullERC.map(clipFeature);

//Find the 90th percentile of the image collection
var ercPrecentile90 = fullERCclip.reduce(ee.Reducer.percentile({
  percentiles:[90]
}));

// define function to extract 90th percentile image at points
var find90piErc = function(feature) {
 var ercpi90 = ercPrecentile90.reduceRegion({
    reducer:ee.Reducer.first(),
    geometry: feature.geometry(),
    scale:400
  });
  return feature.set({erc90pi:ercpi90.get('erc_p90')
  });
};

// map function on to the fire features
var firefeature_pi_added = firefeature.map(find90piErc);

// print first feature
print('Fire Features with 90th Percentile', firefeature_pi_added);

Here is the solution that I came up with to select 1000 random images:

var fullERCselect = ee.ImageCollection(fullERC.randomColumn('random',13).sort('random').limit(1000));


texts:
In my styling file I have a value map with a certain order. But whenever I load it to QGIS the order changes. For example in my styling file I defined this.

<edittype widgetv2type="ValueMap" name="group">
      <widgetv2config fieldEditable="1" constraint="" labelOnTop="0" constraintDescription="" notNull="0">
        <value key="" value=""/>
        <value key="A" value="C"/>
        <value key="B" value="D"/>
        <value key="C" value="A"/>
        <value key="D" value="B"/>
      </widgetv2config>
    </edittype>

But if I load the (*.qml), my value map has another order. If I export it the code changes to this.

<edittype widgetv2type="ValueMap" name="group">
      <widgetv2config fieldEditable="1" constraint="" labelOnTop="0" constraintDescription="" notNull="0">
        <value key="" value=""/>
        <value key="A" value="A"/>
        <value key="B" value="B"/>
        <value key="C" value="C"/>
        <value key="D" value="D"/>
      </widgetv2config>
    </edittype>

Is this a bug or does QGIS automatically change the order to alphabetically?
QGIS version 2.18.14


texts: Getting only time from Python parser of ArcGIS Field Calculator?

texts:
Using ArcGIS 10.4, how I can modify the Date() functions on VB or Python scripts of field Calculator

Now (  )
datetime.datetime.now( )

to get only time like

  12:32


texts: ee.Feature extraction from ee.FeatureCollection by Feature index

texts:
I have a tabular asset in Google Earth Engine that has two features. To extract them one by one in a for-loop I had to introduce site_id property (column) and filter based on that value:
for (var i = 1; i < 2; i++) {
  var feature = myAsset.filterMetadata('site_id', 'equals', i).first();
};

How can one extract one ee.Feature directly using the Feature Index which is present in the asset?
I am looking for an analogy of ee.FeatureCollection.first() but .second(), .third() etc


texts: Which coordinate system is the best for a distances from whole USA?

texts: I have about 100 cities all over the USA and want to calculate some distances. According to that I want use one coordinate projection and need to know which one is the best so I can get the most accurate results. I have compared the results using NAD83 for Texas and NAD 83 for Washington for a same data and the results are not the same.


texts: GeoServer REST API: How to get list of workspaces by username

texts: Is it possible to get list of available workspaces and list of data stores in some workspace for some specific user with GeoServer REST API by sending request with admin user credentials?


texts: Postgis : Geometry vs 3 geom columns in single table

texts:
Requirement:

There a specific requirements where I need a geometry column in table to be flexible ( It can line or point or polygon or all of it )

Justification for requirements :
main concern is to avoid the complexity of the having too many tables for each geometry type as we will be having many categories.

Possible solution I think can be are :

either use geometry column or create 3 different geom columns for point, line, polygon on single table

Has anyone been through this situation before ?

what will be possible  pros and cons of each of above 2 implementation strategy ?


texts: Using Sqlite, Python 2.7 and Spatialite?

texts:
I am trying to spatially enable my sqlite database:

import sqlite3
db = r'c:\temp\mydata.sqlite'
conn = sqlite3.connect(db)
conn.enable_load_extension(True)
cur = conn.cursor()
cur.execute(r"SELECT load_extension("<path>\libspatialite-4.dll");")
conn.commit()
conn.close()
del conn

Python just stops working when I run this code.  It doesn't crash in windows 7 x64, it just closes.  It's very strange.  No error is raised.

When I run the commands directly from python.exe in the command propmpt, I see python crashing:

Am I doing anything wrong?


texts: Image Overlay / Qgis raster projected incorrectly on Leaflet JS

texts:
In leaflet js I added an image to my map (image overlay).
The image is a georeferenced jpeg, and it's a raster generated by Qgis.
In leaflet I setted the image bound points correctly (NE and SW), and I tested this even after the map is rendered.
So, if the extreme points (NE and SW) are corrects, the internal image points do not correspond to the real coordinates, as you can see in the image (see the detail of railway tracks):

Any clue?


texts: Describing input field in ArcMap Field Calculator window?

texts:
A few days ago I was answering Comparing two fields and printing result as 1 (if number matches) or 0 (if not) using Python Parser of ArcGIS Field Calculator?

When describing the process required to get the result using the Codeblock I talked about the "Pre-Logic Script Code", but then didn't know what to call the next box in the window, so I just said

  and in the bottom section enter e

What is that bottom section actually called, so I can use the right term?

[]


texts: How ensure data put into PostGIS gets proper CRS attached?

texts:
I recently had an intractable problem that we finally traced back to the fact that some points in a PostGIS table had a CRS associated with them and others did not.  This resulted in some points being displayed in QGIS while other were not, even when they were they were in different layers and all the coordinates were actually using the same CRS.

I ended up solving this problem by using postgis functions ST_SetSRID and ST_TRANSFORM to set all the point to the same CSR but I am still not sure if these functions correctly set the CRS and whether or not I will not have the same problem after I next add data into the table.

My question now is how do I make sure that all data I put into PostGIS gets the proper CRS attached to the_geom?

Doing this has two advantages as I see it:

you don't get prompted to set the CRS each time you build a layer from the table (I was puzzled why this happened with my 'old' data).
you don't get wierd problems with QGIS when you try and build layers from the data in the tables ;)

The points that had no CRS associated with them were all entered into PostGIS using the SPIT plugin from shape layers that I had saved to disk.  Recently some kind soul on stackGIS pointed out that I could short circuit the process by using the database manager to save a layer to PostGIS.  This appears to save layer CRS with the data into PostGIS.

From now on I intend to use the database manager to save data into PostGIS an the assumption that this will protect me from headaches like this in future.  Is this a reasonable approach?

Apart from getting answers to the question I was also motivated to write this up as an easier to read summary of the question and answers referenced above.  In getting to the answers there we found several red herrings that have confused the basic issue and I have left those out here.


texts: Making a layer read-only in arcgis 9.3.1

texts:
Is there a way to make a layer (a independent shapefile or a layer in a SDE gdb) read-only - so that people will be able to see it but not edit it?

(I'm using 9.3.1)


texts:
I'm pretty new to Earth Engine and I'm trying to export several feature collections.
However for one of those this error keeps occurring whenever I try to run the Export process:
Error: Primary geometry of feature 'null' has 1100253 vertices, above the limit of 1000000 vertices.

Is there any way to reduce the vertices of said feature?


texts:
I am new to QGIS and am trying to use the print composer.

My questions are:

How to change from vertical to horizontal orientation.
How to zoom in on a map that I've added.
How stop a black box from appearing every time I move my map.


texts: Referencing variables inside ST_Distance_Sphere function

texts:
I trying to create an UPDATE function that will ultimately return straight-line distance to the nearest node on an OSM network (which will be later added to pgrouting distances I have already calculated).

This is the basic function (with actual coordinates) that I want to replicate:

SELECT st_distance_sphere/1000
FROM ST_Distance_Sphere((SELECT the_geom
                        FROM ways_vertices_pgr
                        ORDER BY the_geom <-> ST_GeometryFromText('POINT(22.685133 41.554309)')
                        LIMIT 1),
                       (SELECT *
                        FROM ST_GeometryFromText('POINT(22.685133 41.554309 )')));

which returns distance in kilometers perfectly. When I insert variable names into the POINT(X Y) bit, it attempts to parse the variable names. I get the same when trying to calculate the constituent start and end points of the ST_Distance function. Using more or less this code:

UPDATE start_end
SET nearest_dist = (SELECT st_distance_sphere/1000
                    FROM ST_Distance_Sphere(
                        (SELECT the_geom
                         FROM ways_vertices_pgr
                         ORDER BY the_geom <-> ST_GeometryFromText('POINT(start_end.y1 start_end.x1)')
                         LIMIT 1),
                       (SELECT *
                        FROM ST_GeometryFromText('POINT(start_end.y1 start_end.x1)')));

Is there a ready way to adapt the first function using variable names? or does this run deeper than syntax?


texts: Add objects from one layer to another at QGIS

texts:
How can I add some objects from one layer to another one.
I mean, I have one layer with specific data(e.gg polyline)
at the same time I have another layer with different objects .how can I copy/move object between the layers.
Thanks
Hamid


texts: Changing size of window when calculating slope in ArcGIS Pro

texts:
I am running a focal statistic using a 2m resolution raster with a 8x8 window (so 16m x 16m) to calculate the "range" across my tile. I also want to do the same with slope but by default slope uses 3x3 to calculate. The issue is that running slope on my 2m raster is not capturing as much slope as it should be. I could resample the raster to a 5m resolution but then I am losing small details I am trying to capture which is why I am using 2m in the first place.
I want to be able to run slope and focal stat (range) over the same size windows with my 2m data.
Thoughts?


texts: R or Python geometry checker that outputs results as table

texts: Is there a R or Python solution that 1) checks the geometry validity of vectors and 2) outputs the results in form of a table (see table below). QGIS for example has a nice plugin that carries out several geometry checks. But how can I carry out these checks within a R or Python script?


texts: GeometryServer > Project with Datum Transformation

texts:
I am facing the problem that when i use the extent in the url and i go from Amersfoort(WKID 4289) to WGS84(WKID 4326) that without the tranformation the point from longitude and latitude is placed about 100m away from the original point.
In our company this isn't what we want.

In ArcGIS 10.1 it is possible in the rest interface to apply the transformation, but i get a repeated error, when i remove the transformation i get no error.

I use the Transformation: esriSRGeoTransformation_Amersfoort_To_WGS1984

Can you tell me what i am doing wrong, i put a image with this post and a sample url:

https://ourserver.com/arcgis/rest/se...d=false&f=html


texts: Joining attributes from both sides of a line

texts:
I have a layer with street segments that I would like to join with a building + population layer. Desired result: The segment contains population data. Nearest neighbor (NNjoin) finds the closest building and adds its population to the attribute table. Issue: there are buildings on both sides of the segment.
Is there a way I can join attributes from several buildings (within a radius, for example) to the segment, adding their values together?
I use QGIS and GRASS GIS. Any ideas?


texts: Overlaying satellite image tileLayer over OpenStreetMap in LeafletJS?

texts:
I am currently creating a web app for visualizing disaster impact via satellite imagery. I have satellite images (.tif) which I need to overlay over OpenStreetMap. Since .tif file is very large, I used magic slicer to break down image into multiple tiles and used leaflet-deepzoom for loading the tileLayer.

Following is the code for loading OSM layer:

var map = L.map( 'map', {
        center: [28.3949, 84.1240],
        minZoom: 2,
        zoom: 4
    });

 L.tileLayer( 'http://{s}.tile.openstreetmap.org/{z}/{x}/{y}.png', {
        attribution: '&copy; OpenStreetMap',
        subdomains: ['a','b','c']
    }).addTo( map );

This works fine for loading OSM but I need to overlay my satellite image as next tileLayer:

L.tileLayer.deepzoom('tiles/harvey_impact_files/',{
        width: 12800,
        height: 5376,
    }).addTo(map).setOpacity(0.5);

This loads both layers but in a disoriented way. I need to correctly position the satellite image layer. Is there any possible way to provide geocoordinate bounds (topLeft latlng and bottomRight latlng) of this tiles/harvey_impact_files/ image so that this image loads exactly over OSM layer??


texts: Determine the heading or azimuth angle of satellite Meteor M2

texts:
I am working on a project Meteor decoder, where I need to calculate the pixel geolocations of the satellite image to draw map overlays. I have all the angles except the azimuth angle. I tried to calculate it from orbital elements, I took two relatively close coordinates at two different timestamp and calculated the angle between them. But this solution is not perfect.
Sample code for calculating every 10th pixel's geolocation, I hope it is clear:
void PixelGeolocationCalculator::calcPixelCoordinates()
{
    double angle;
    CoordGeodetic satOnGroundPrev;
    CoordGeodetic satOnGround;

    std::vector<Eci> passList;

    DateTime passEnd = DateTime(mPassStart).Add(mPassLength);

    for(DateTime currentTime = mPassStart; currentTime <= passEnd; currentTime = currentTime.AddMicroseconds(PIXELTIME_MS * 10 * 1000))
    {
        passList.push_back(mSgp4.FindPosition(currentTime));
    }

    std::vector<Eci>::const_iterator it = passList.begin();
    std::vector<Eci>::const_iterator prevIt = passList.begin();
    ++it;

    DateTime currentTime = mPassStart;
    for (unsigned int i = 0 ; it != passList.end(); prevIt = it,  ++it, i++)
    {
        satOnGroundPrev = Eci(currentTime, prevIt->Position(), prevIt->Velocity()).ToGeodetic();
        satOnGround = Eci(currentTime, it->Position(), it->Velocity()).ToGeodetic();
        currentTime = currentTime.AddMicroseconds(PIXELTIME_MS * 10 * 1000);

        angle = calculateBearingAngle(satOnGround, satOnGroundPrev);
        angle = degreeToRadian(90) - angle;

        for (int n = 79; n > -79; n--)
        {
            CoordGeodetic coordinate (los_to_earth(satOnGround, degreeToRadian((((mAlfa) / 79.0)) * n), 0, angle));
            mCoordinates.push_back(coordinate);
        }
    }
}

Any better way to calculate this?


texts: Position of points in a QGIS vector layer

texts:
I have created a new layer overlaying a map, then chosen the "add feature" tool to add points, which I have named with numbers, representing the location of customers.
I need to have the coordinate positions of my points. Where do I find them? I have tried Layer-Open attribute table but this only contains my number-names for my points.


texts: Importing TIFF file for use as basemap in ArcGIS Pro

texts:
I have purchased an Ordnance Survey map for use in ArcGIS Pro as a basemap. The file that I received came in a zip file and contained the a TFW File, TIFF File, OVR File and an XML Document. When I open the TIFF file normally in Windows, it is a high resolution file (73 mb). However, when I "Add data" in ArcGIS Pro and import the TIFF file, it is very blurry as a base map and is not high resolution.
Are there any special tips to follow when importing a TIFF file for use as a basemap in ArcGIS Pro?


texts:
When I try to merge two lines selecting them and using the Merge Tool in the Advanced Editing Menu I get an error which says that merging these lines would result in an incompatible geometry for this layer so that it is cancelled.

How can I merge lines?


texts: Catch Execute Error and Continue (arcpy)

texts:
I've got a script that does a number of Date/Time conversions using arcpy.ConvertTimeField_Management. A common issue that I need to handle is when the field is empty and an execute error is raised:

arcgisscripting.ExecuteError: ERROR 000192: Invalid value for Time Format = 'yyyyMMddHHmmss' - was not able to successfully convert any row's time field value

This error occurs when I use the try/except below; I'd like the script to ignore this error and continue. Any ideas?

try:
    results = arcpy.ConvertTimeField_management(SA_TEN_GOV_els, "START_DATE", "yyyyMMddHHmmss", "START_DAT", "DATE", "'Not Used'")
except:
    print results.getMessages()
    pass


texts:
As part of an internal "learn at work month" program, myself and a colleague have volunteered to do a presentation: "The Wonderful World of Web Maps".

To this end, we're now looking for examples of wonderful, weird, functional, beautiful, and just plain awesome web maps that are out there. Any online map that is either interactive or dynamically updates itself would qualify.

Part of this question is already answered with the What Makes a Map Beautiful? question which has some-web maps as does What Makes a Map be classed as Badly Designed? but we're more interested in a diverse range of maps rather than just the beautiful and bad.

I've done a little searching and come up with a couple (that just happen to have a transport theme) like:

The Live Ship Map

Flight Tracker (Heathrow Airport)

What other web-maps are out there that might warrant demonstrating or could use wider exposure?


texts: Using ArcPy to get all database names within server?

texts:
I am trying to populate a folder with ArcSDE database connection files for all the databases within a server instance.  I need to be able to do this for several of our servers.

Although I know within SQL Server how to get a list of DB names from the sys.databases table, using the arcpy function CreateDatabaseConnection I cannot specify "sys" or "master" to make the connection.  Therefore, I need to know the name of a database currently in the server.  I am not sure how to make the check without hard-coding a name.

import arcpy
arcpy.env.overwriteOutput = True

#Input the location where the database connection files will be stored
connectionLocation = "F:/FolderName"

#Get the list of servers
arcpy.CreateDatabaseConnection_management(out_folder_path= arcpy.env.scratchFolder,out_name= "GIS_USERCONFIG" ,database_platform="SQL_SERVER",\
                                          instance = "XXXXXXX" ,account_authentication="DATABASE_AUTH",username="user",password="password",\
                                          save_user_pass="SAVE_USERNAME",database= "GIS_USERCONFIG",schema="#",version_type="TRANSACTIONAL",version="dbo.DEFAULT",date="#")
sdeConn = arcpy.ArcSDESQLExecute(arcpy.env.scratchFolder + "/GIS_USERCONFIG.sde")
serverListSQL = sdeConn.execute("Select name from userservers")
serverList = []
for i in range(0,len(serverListSQL)):
    serverList.append(serverListSQL[i][0])

#Make a folder for each server instance and build the connection files for all DBs in the server instance
for server in serverList:
    folder = connectionLocation + "/" + server
    if arcpy.Exists(folder):
        print server + " folder already created"
    else:
        try:
            arcpy.CreateFolder_management(connectionLocation, server)
            print "Created Folder " + server
        except:
            print "failed to create folder " + server

    #
    arcpy.CreateDatabaseConnection_management(out_folder_path= arcpy.env.scratchFolder,out_name= server ,database_platform="SQL_SERVER",\
                                          instance = server ,account_authentication="OPERATING_SYSTEM_AUTH",\
                                          database= "master",schema="#",version_type="TRANSACTIONAL",version="dbo.DEFAULT",date="#")
    sdeConn = arcpy.ArcSDESQLExecute(arcpy.env.scratchFolder + "/" + server + ".sde")
    dbListSQL = sdeConn.execute("Select name from sys.databases")
    dbList = []
    for i in range(0,len(serverListSQL)):
        dbList.append(dbListSQL[i][0])
    for db in dbList:
        arcpy.CreateDatabaseConnection_management(out_folder_path= folder ,out_name= db,database_platform="SQL_SERVER",\
                                          instance = server ,account_authentication="OPERATING_SYSTEM_AUTH",\
                                          database= db,schema="#",version_type="TRANSACTIONAL",version="dbo.DEFAULT",date="#")

Basically my issue is with the 2nd call for the function "CreateDatabaseConnection_management" since it needs a database that exists in the database (and I could not use "master or sys").  If you open this as a tool in arcmap, the input is a drop-down box that populates choices based on the databases in the chosen server instance.  Is there a way to do this in the code?


texts: Missing field in address locator when geocoding

texts:
I'm a newbie to creating address locators, and I'm hoping someone here can help with a problem I'm having. I'm using ArcMap 10.1.

I created a US - single house address locator using a shapefile of 911-based addresses. I created the address by right-clicking on my geodatabase in the catalog window ---> New---> Address Locator. When creating the locator, I included information for the House Number, Street Name, and Suffix Type fields. The locator appears to be created correctly, as it is showing in the geodatabase in the catalog window, and I can select the address locator through File-->Add Data-->Geocoding-->Geocode Addresses.

After I click on Geocode Addresses and select my address locator, the Geocode Addresses window opens correctly, and I'm able to select my excel table as the Address Table. At that point, though, instead of showing fields for House Number, Street name, and Suffix Type, I'm given only a field/drop-down box for Street. I'd never seen this before, so I tried to geocode my addresses using the concatenated address that contained the Street Number, Street Name, and Suffix Type, which didn't work. I also tried to geocode the addresses using just the Street Name field, which again didn't work (which I expected...). Any ideas as to what might be causing this, or possible fixes?


texts: Can't iterate feature in while nested iteration

texts:
I need to do some operation with every combinations of features in two FeatureClasses. The "for" iterating operator is very very slow-working, so I decided to try "while" operator and using FeatureCursor. So nested iteration I have like below:

    IFeatureCursor fcur1 = fc1.Search(null, false);
    IFeatureCursor fcur2 = fc2.Search(null, false);

    IFeature f1 = fcur1.NextFeature();
    IFeature f2 = fcur2.NextFeature();

        while (f1 != null)
        {
            //external cycle code
            MessageBox.Show("external cycle FID " + f1.Value[0].ToString());

            while (f2 != null)
            {
                //internal cycle code
                MessageBox.Show("internal cycle FID " + f1.Value[0].ToString());
                f2 = fcur2.NextFeature();
            }
            f1 = fcur1.NextFeature();
        }

For example, there are 3 features in each FeatureClass fc1 and fc2.

I use a MessageBoxes to show FID of feature f1 to understand what happens in external and internal iterations.

The expected order of MessageBoxes is:

external cycle FID 0
 internal cycle FID 0
 internal cycle FID 0
 internal cycle FID 0
external cycle FID 1
 internal cycle FID 1
 internal cycle FID 1
 internal cycle FID 1
external cycle FID 2
 internal cycle FID 2
 internal cycle FID 2
 internal cycle FID 2

But result of my code is:

external cycle FID 0
 internal cycle FID 0
 internal cycle FID 0
 internal cycle FID 0
external cycle FID 1
external cycle FID 2

So, external cycle operates with every feature f1, but internal cycle gets the first feature f1 only. Why another features f1 doesn't go to internal cycle?


texts: I am unable to combine Numpy arrays in ArcMap when using FeatureClassToNumPyArray

texts:
I am using FeatureClasstoNumPyArray to get information from two features classes

import arcpy
import numpy

sspipe_fcs = 'ssGravityMain'
sspipe_fl = ('ORIGINAL_COST', 'FEATURE', 'DEVELOPMENT','LENGTH_M4')
ssarr = arcpy.da.FeatureClassToNumPyArray(sspipe_fcs, sspipe_fl)
print ssarr
print ssarr.shape
print ssarr.dtype

wlpipe_fcs = 'Water Mains'
wlpipe_fl = ('ORIGINAL_COST', 'FEATURE', 'DEVELOPMENT','Shape_Length')
wlarr = arcpy.da.FeatureClassToNumPyArray(wlpipe_fcs, wlpipe_fl)
print wlarr
print wlarr.shape
print wlarr.dtype

numpy.concatenate ((ssarr, wlarr))

I would like to combine both arrays but I get an error

TypeError: invalid type promotion

My overall goal is to use the array to eventually prepare some type of pivot table using pandas.

Here are the results of those print statements

[(3087.622879285057, u'Pipe', u'shopping center', 90.81243762603108)
 (4308.758177145577, u'Pipe', u'shopping center', 126.72818168075227)
 (2079.15277816955, u'Pipe', u'shopping center', 61.15155229910441)
 (2055.8024028091863, u'Pipe', u'shopping center', 60.464776553211365)
 (3825.8261621564598, u'Pipe', u'shopping center', 112.5242988869547)
 (5056.951172625341, u'Pipe', u'shopping center', 148.7338580183924)
 (4143.872612491386, u'Pipe', u'shopping center', 121.87860624974664)
 (4545.80563630292, u'Pipe', u'shopping center', 133.70016577361528)
 (3010.085501040413, u'Pipe', u'shopping center', 88.53192650118862)
 (3139.2163487103285, u'Pipe', u'shopping center', 92.32989260912731)
 (3251.2079792976815, u'Pipe', u'shopping center', 95.62376409699063)
 (2104.7324983994613, u'Pipe', u'shopping center', 61.90389701174887)
 (2894.457469764144, u'Pipe', u'shopping center', 85.1311020518866)
 (1922.1176673586654, u'Pipe', u'shopping center', 56.53287256937251)
 (3746.678222477405, u'Pipe', u'shopping center', 110.19641830815897)
 (2629.9600545841795, u'Pipe', u'shopping center', 77.3517663112994)
 (7509.473519849657, u'Pipe', u'shopping center', 220.86686823087229)
 (4927.328877910251, u'Pipe', u'shopping center', 144.92143758559564)
 (5740.15444723423, u'Pipe', u'shopping center', 168.82807197747735)]
(19,)
[('ORIGINAL_COST', '<f8'), ('FEATURE', '<U50'), ('DEVELOPMENT', '<U50'), ('LENGTH_M4', '<f8')]
[(28937.625507077002, u'Pipe', u'shopping center', 826.7893002022)
 (11969.306667602146, u'Pipe', u'shopping center', 341.98019050291845)
 (21051.785700918437, u'Pipe', u'shopping center', 601.4795914548125)]
(3,)
[('ORIGINAL_COST', '<f8'), ('FEATURE', '<U10'), ('DEVELOPMENT', '<U50'), ('Shape_Length', '<f8')]


texts:
I have two LiDAR datasets; one of 1m resolution and another of 2m resolution. I want to stamp 1m LiDAR data on 2m LiDAR.

Similar tool is available in Vertical Mapper in MapInfo called as Splicer where they have an option of stamping two rasters.

Is there any way of stamping two rasters in QGIS?


texts: Save Popup Configurations in ArcGIS Online

texts: Is there a way to save configured pop-ups within ArcGIS Online and then use the saved configuration in other ArcGIS Online web maps?


texts:
I have installed the semi-automatic classification plugin (SCP v. 7.7.1) on QGIS 3.4.11 (Madeira) using Kubuntu 18.04.
I downloaded Sentinel 2 (Level 2A) images and wanted to test one of the semi-automatic classification procedures. However, I failed to create a training input file (*.scp) getting the following information on my QGIS main window:
Information [4]: Select a raster; raster is not loaded

Furthermore, the following information is registered in the log file:

I have the suspicion that there is something not working between QGIS and other libraries such as GDAL. So, how can I inspect the issues more thoroughly and how can I resolve it?


texts: How to cut parallelogram with intersecting line?

texts:
I have a very specific set of circumstances under which I need to break a parallelogram into 2 pieces. The parallelogram obviously always have 4 vertices and the intersecting line intersects the polygon exactly 2 times. In the order in which the intersecting line points are given, the polygon area to the right of the line is the area of interest. I feel this could be solved using a linked list implementation in code but am unsure where to start. Additionally, if there is a simply geometry function to get the job done, that is acceptable too (I've tried intersection and difference but neither applies to the situation). An example is attached below.

Notes:

Polygon vertices are given in clockwise order starting from the upper-right-most point.
The intersecting line can have 2 or more points. (In the case of 2, that means that the same intersecting ling segments intersects 2 of the parallelogram's sides.
It's possible the intersecting line only intersects one side of the polygon.

The area of interest is the upper-right area of the polygon that is closed off by the line. This implies that in this case, the intersecting lines' vertex order is 1) the point to the right, 2) the point in the center, 3) the point at the top.

These are a few more possible cases where the area of interest is shaded (please excuse my poor drawings!).


texts: Edges of Surface Difference tool output are invalid

texts:
See this image, an example output of the ArcGIS Surface Difference tool (gives cut/fill volumes):

I have occasionally come across the illustrated problem where the output does not cover the full extent of where the surfaces overlap. The black line is the edge of the two TINs submitted to the Surface Difference tool. One TIN is the topographic surface, the other is the bottom of oil surface from dip data. In the red circled areas in the image you can see where the tool output polygon does not extend to this edge. Both TINs had been soft clipped with this polygon before being used in the tool. I have tried hard clipping them, same problem. One of them is a topographic surface and so is valid outside of the clipping polygon. I therefore tried it without first clipping this TIN. Same issue. Another example:

When I clipped both surfaces to the same extent the problem was actually reduced but did not disappear, and the volume reported increased, to me indicating that there was indeed a problem introduced by the inconsistencies at the output edge. I ended up being able to fix this by buffering the wet area polygon by 0.01m. Then I softclipped the bottom of oil surface with this. Problem solved:

What is the reason for this solving the problem? I'm guessing it is because the two surfaces intersect, rather than one meeting the other. Through intersecting it is clear what lies between the two surfaces. If they just meet, and the way that this is expressed is influenced by rounding or something else to do with number precision, perhaps the tool will find a smaller area within which the two surfaces overlap?

Arcgis 10.2.0.3348


texts: Shapefile layer cannot be generated - blank warning

texts:
Installed QGIS 3.0 (w/Python 3.6.5) to replace v2.X.
Started a new project - imported shapefiles and everything looks good.

Problem: Click "Create New Shapefile Layer" and I get the pulldown menu (file name, geomeotry type, etc.) Fill everything out and click "ok"
I get this:

  Layer creation failed. Please check message log for further information.

Message log says this:

  2018-04-13T09:04:34 Warning:

That's it. Warning is blank.
Under the OGR tab I've got:

  2018-04-13T09:04:34    WARNING Creation of OGR data source SG sites.shp failed: Failed to create file SG sites.shp: Permission denied.

Why is my permission being denied.  Tried a full reinstall but results are the same.


texts: QGIS Composer: Force loading of all tiles on export

texts:
I face the problem that not all tiles are loaded on a map export via the QGIS composer. This happens randomly with all kinds of tile servers provided by the openlayers-plugin.

Any hint how I can force loading the tiles before the image is saved?

I hope this is not a duplicate - I've been looking through existing threads thoroughly but haven't found anything.

Edit:
Specs: well equipped ThinkPad x230 with Ubuntu 12.04 LTS and QGIS 2.8. Internet-connection not very good at the moment.


texts: Change display name for feature layer when using LayerList widget

texts:
I'm building an app with the ArcGIS JavaScript API, version 3.16. I have a map with two feature layers and I'd like to use the LayerList widget so the user can toggle them on and off. It's functional as I have it written, but not ideal in a couple of ways. The layers appear in the layer list as "graphicsLayer1" and "graphicsLayer2", and I'd like to find a way to give them more descriptive names. Also, the layers are visible when the box is unchecked, and become invisible when I check the box, and I'd like to flip that around.

The relevant section of code is below.

    var mapMain = new Map("mapDiv", {
        basemap: "topo",
        center: [-120, 38],
        zoom: 6
    });

    var resLyrLn = new FeatureLayer("https://data.farwestern.com/arcgis/rest/services/scratch/test_SD_RelatedResources/MapServer/0", {
        mode: FeatureLayer.MODE_ONDEMAND,
        visible: true,
    });
    var resLyrArea = new FeatureLayer("https://data.farwestern.com/arcgis/rest/services/scratch/test_SD_RelatedResources/MapServer/1", {
        mode: FeatureLayer.MODE_ONDEMAND,
        visible: true,
    });

    var lyrs = [resLyrArea, resLyrLn]

    mapMain.addLayers(lyrs);

    var layerListDijit = new LayerList({
        map: mapMain,
        layers: lyrs
    }, "layersDiv");
    layerListDijit.startup();


texts: Appending OS MasterMap feature class to the feature class in PostgreSQL database

texts:
I have my Enterprise geodatabase as PostgreSQL. I want to keep OS MasterMap data on the geodatabase and would like to update it whenever a new tile is downloaded. However, I am struggling to find a way to to append the existing OS MasterMap feature class with new data. I tried OS Translator II in QGIS as the OS Master Map files come with .gz extension but it overwrites the existing feature class in the geodatabase.

Could you give me pointers as to how this can be done? Would I need a Python script or shall I first convert OS MasterMap data into shapefiles every time and then think about adding it to the existing OS MasterMap through ArcGIS.

I would preferably like the solution in QGIS as most of my team-mates have ArcGIS Basic license except me who has ArcGIS Standard license.


texts: Shapefile name specific failure with Append in ArcGIS 10.2

texts:
Can anyone tell me why the append (Desktop Management) function is not working on my 64 bit Windows 7 machine with ArcGIS 10.2.2 under the specific conditions described below?  All tests on my machine were run with Windows 7 10.2.2, except for one test specified below.  Other machines used Windows 7 and 10.1 or 10.2.2  I am unsure if the other machines were 64 or 32 bit.

The specific conditions are:

When the append function is run with two shapefiles, one named sites_p, the append tool runs successfully, but no records are appended to sites_p.  If I rename it to sites_u, sites_t, or anything else, records are appended. This error only happens when sites_p is selected in ArcMap using the dropdown selector in the Append tool.  If I add sites_p as the target dataset by selecting the filepath in ArcMap, the append function works.  It also works correctly when run from ArcCatalog

The append tool was re-run after each of the steps described below.  None of the steps resolved the issue:

The shapefiles were moved to a directory 1 level below my C drive with no success.  All the field types and widths were checked.  A blank .mxd file with only the two shapefiles was used. A new sites_p shapefile was created and tested. .net 3.5 was re-installed. I performed a clean uninstall/re-install as directed by ESRI tech support. As part of the uninstall, I also removed 64 bit background processing.  All tests were repeated after the re-install

The append tool works correctly when:

It is run on other machines with the sites_p shapefile.  The sites_p shapefile is converted to a file geodatabase named sites_p.  It is run with ArcMap 10.2.2 and the sites_p shapefile from Microsoft Windows Technical Preview on a separate partition of my hard drive.


texts:
I want to sum the values of points within each polygon in R, but not exactly sure how to do it. I have computed counts of points, but not a summarization of values of the points within each polygon.
I am also not sure how to generate toy code; sorry for the lack of an example.


texts: Is there any way to fast read pixel value using gdal in C#

texts:
I create with C# a series of lat/lon value pair in the code(using i, j loop) and I want to extract the pixel value based on either lat/lon pair or pixel location (i, j).

I know there is one gdal utility called gdallocationinfo which can read pixel value very fast. But it is written in python.

I am wondering whether there is a way to use this gdal utility in my C# code, or other methods that can extract pixel value in a faster manner.

Final solution:

static void Main(string[] args)
    {
        Console.WriteLine(DateTime.Now.ToString());
        Gdal.AllRegister(); //Register all gdal data drives
        List<int> resolutionId = new List<int>() { 9 };

        Dataset ds = Gdal.Open("D:/SKnoZero.tif", Access.GA_ReadOnly);   //Read raster
        double[] gt = new double[6];
        ds.GetGeoTransform(gt); //Read geo transform info into array
        int Rows = ds.RasterYSize;  //Get the number of rows
        int Cols = ds.RasterXSize;  //Get the number of columns
        Band band = ds.GetRasterBand(1);    //Read band 1
        double startX = gt[0];  //Upper left lon
        double startY = gt[3];  //Upper left lat
        double interval = gt[1];    //Cell size
        double x, y;    //Current lon and lat

        var filePath = "D:/SKcsv.csv";
        using (StreamWriter csv = new StreamWriter(filePath))
        {
            for (int k = 0; k < Rows; k++)  //read one line at a time
            {
                y = startY - k * interval;  //current lat
                int[] buf = new int[Cols];
                //ReadRaster parameters are StartCol, StartRow, ColumnsToRead, RowsToRead, BufferToStoreInto, BufferColumns, BufferRows, 0, 0
                band.ReadRaster(0, k, Cols, 1, buf,Cols,1,0,0);
                //iterate each item in one line
                for (int r = 0; r < Cols; r++)
                {
                    if (buf[r] != -32768)   //if pixel value is not NoData value
                    {
                        x = startX + r * interval;  //current lon
                        var geoidpair = GetVrgId(x, y, resolutionId);   //get geoid
                        var temp = geoidpair;
                        var newLine = string.Format("{0},{1},{2},{3}{4}", x, y, buf[r], temp[0], Environment.NewLine);
                        csv.Write(newLine);
                     }
                 }
         }
         Console.WriteLine(DateTime.Now.ToString());
         Console.WriteLine("Press any key");
         Console.ReadLine();
    }


texts: Changing default or bulk to show Thousands Separators in ArcGIS Pro attribute table?

texts:
I am creating a Map Series in ArcGIS Pro, and have a number of dynamic text elements populated by numeric fields in my map series layer.  The default for numeric values appears to be to not use thousands separators, however I would like to see them, at least in my dynamic text elements.

Now I know I can change the Number Format from the attribute table, by right-clicking a field and then clicking the ellipsis next to "Numeric", and checking on Show thousands separators

But this only changes it for the one field.  I would like this to be the default setting in my project, or at least be able to change it in bulk for the multiple fields I need to display.

How can I change all my (selected) fields to show thousands separators?

I have looked through the text element formatting options and can't see anything in there either.


texts:
I have a fishnet polygon feature class and a land cover (40 categories) raster. Both cover a large area.

I seek to:

a) calculate the area of each land cover type category (40) inside of each Grid from the fishnet,

b) Overlay those two datasets in order to have just 1 dataset. What´s the best option to have just a layer showing the Land cover type inside of each Grid?

c) Clip the Grid for just work with the areas of the grid that match the land cover map. Yet I do not want an exact clipped grid but instead that the fishnet grid from the raster boundaries would not get clipped, thus, maintaining the full grid configuration in overlaid areas.

What is the best approach to reach the goals for a), b) and c)?
Additionally, is it possible to do it without converting the raster to polygon or polygon to raster?


texts: How to create 3D parallel lines with Z coordinates as same as Original 3D polyline using Arcpy?

texts:
I want to calculate cut/fill cost for road alignments. For this reason i need to create a buffer-like polygon from given road centre-line but this polygon has to have Z-values that are taken from corresponding Z coordinates of polyline (check my recent question for more details).
Also my project is to optimize road coordinates to decrease cut/fill cost with genetic algorithm, so manual buffer creation will not be feasible.
After some search for finding solution to the problem i have found this question. The answers from @FelixIP and @jbalk to the question, inspired me the following methodology:

1.Create two 3D parallel polylines with same Z-values

2.Convert the parallel polylines to polygon

The mentioned script is:

import arcpy, math
infc=r'..\SCRARCH\clone.shp'

def CopyParallel(plyP,sLength):
    part=plyP.getPart(0)
    lArray=arcpy.Array();rArray=arcpy.Array()
    for ptX in part:
        dL=plyP.measureOnLine(ptX)
        ptX0=plyP.positionAlongLine (dL-0.01).firstPoint
        ptX1=plyP.positionAlongLine (dL+0.01).firstPoint
        dX=float(ptX1.X)-float(ptX0.X)
        dY=float(ptX1.Y)-float(ptX0.Y)
        lenV=math.hypot(dX,dY)
        sX=-dY*sLength/lenV;sY=dX*sLength/lenV
        leftP=arcpy.Point(ptX.X+sX,ptX.Y+sY)
        lArray.add(leftP)
        rightP=arcpy.Point(ptX.X-sX, ptX.Y-sY)
        rArray.add(rightP)
    array = arcpy.Array([lArray, rArray])
    section=arcpy.Polyline(array)
    return section

with arcpy.da.UpdateCursor(infc,("Shape@","Width")) as cursor:
    for shp,w in cursor:
        twoLines=CopyParallel(shp,w)
        cursor.updateRow((twoLines,w))

However the script in the answer doesn't yield in 3D polyline and i couldn't figure out how to manipulate the script to get my goal.
What changes must be made in the mentioned script to create two 3D parallel polylines?
I have a 3D Analyst license and am using ArcMap/ArcScene/etc with advanced license, however any answer in any platform other than Arcpy, like QGIS could be helpful too.

UPDATE:

Just to clarify my question, i made a 3D model of a simple example in Autodesk AutoCAD and took a screenshot.
In the screenshot image the blue line is road centreline. The white surface is the output of buffer tool from blue line which doesn't help me, also its borders could be generated from mentioned code above. The red polygon, is my desired polygon which i need to calculate cut/fill and yellow lines are subject of this question that i need to create the red polygon.


texts: FME Writing Geo-database relationship classes

texts: I am reading geodatabase using ESRI File reader and this geodatabase contains feature classes with relationship classes. I am trying to write them as it is to another geodatabase. The problem is relationship classes are not created in output geodatabase.


texts: Using FME, Is there a way to transfer Metadata from an SDE FeatureClass to a Shapefile (shp.xml)

texts: Is there a way to transfer, using FME, the metadata from an SDE featureclass to a shapefile?


texts:
I am trying to use difference tool in QGIS to get a mask out the ice coverage from a base vector map. The number of features in ice coverage is around 2000 while my base map has 237 features. I am using

I tried using difference tool to do it but it got unresponsive and stuck at 79%.

I also tried gdifference function to do it in R. But, it doesn't create spatialpolygondataframe but a spatialpolygon when I tried a difference of subset of my vector layer.

I am now trying to do it using a batch process in QGIS differencing 200 features at once. However, it seems to be taking a very long time. The time taken increases with every consequent batch process. Considering that I have to further mask out other layers also, it will probably take a very long time like this.
Is there any other (faster) way to do this?


texts: Why does Print composer draw wrong coordinates?

texts:
QGIS Version: 2.14.3
EPSG: 32630 (WGS 84 UTM Zone 30N)
I´m trying to draw a grid with EPSG: 32630. Layers are in the same CRS. As you see, coordinates are completely off. It does not matter if I use decimals or degree. Doing the same in ArcGIS 10.4 works without problems, so I think the data is fine.
The real coordinates are around 43°N and 3°W. Is this just a bug in the current QGIS version? I did not use QGIS recently, so I do not know if this worked in previous versions like the current LTR (2.8).


texts: Model Builder Error: Python cannot be found

texts:
Building some models in Model Builder for a training session, and it's not working out so well.  Model Builder can't seem to find Python.

Model Builder Error 000574: Python 2.7 is not installed

Here's the kicker...Not only am I looking at Python 2.7 as I type this, but I also ran a more complex model this morning.

Since the first issue this morning I have run a repair on ArcGIS 10.2 and installed Python again, since I was still receiving the same error after the repair.

Windows 7
Data source is read from Database Connection over the network and written locally.
ArcGIS 10.2
Python 2.7.9
I've used several tools and datasets in the model. Both local and netowork sources.


texts: How to record a base layer change event in OpenLayers?

texts: The OpenLayers.Control class possesses an eventListener property which could be used in the layer switcher control. I am not sure how to make use of it to detect a base layer change event so that I can store the variable name of the chosen layer to be used subsequently when the user navigate back to the map page. Can anyone give me some hints?


texts: Converting Degrees Minutes Seconds to Decimal Degrees for coordinates?


texts: What Makes a Map be classed as Badly Designed?

texts:
A poorly designed map can not only look visually unappealing, but can convey the wrong message, which could lead to bad decisions being made.

I would like to ask people to post examples (that are in the public realm) of poorly designed maps, WITH justification on why it is bad design.

Although this 'question' does not have a clear answer, I think it will prove a useful resource to see what merits bad design, so others can learn what NOT to do.  I will let the votes choose the 'right' answer.

I would also like to see examples of bad design around web-mapping.

I would argue that although GIS Professionals generally know how to create a good looking map, I would say that they also have a tendency to over-complicate web-maps, by trying to re-create a GIS desktop application on the web:

Also slow to use, and hard to even get to the map.

I think now that normal people are more used to Google Maps style/simplicity, web-maps should follow a similar approach.

With the explosion of 'NeoGeography' particularly in the web realm, we now have a lot of non-GIS professionals creating maps for the web.  A lot of these developers are often very good at user-interface design, but not trained in cartographic principles.

I think, with web maps, its all about combining the skills of both cartography and user interface design.


texts:
Edit: tl;dr
QGis apparently casts my primary key on postgresql to text when using the identify tool. This causes the form to be slow in startup because it disregards the index on it. A (tested) quick fix is to create an extra index like:
create index on belmap.address USING btree (cast(stableid as text));

I still wonder if there is a proper solution for this.

Original post

I am using a QGIS 3.10.3 (but also tested with QGIS 3.4) directly connected to a PostgreSQL 11 with PostGIS.

I am building a project and data model on which our data operators can do asset management. The project uses several large tables (5.000.000+ rows per table, 4-7 tables). I limit the amount of features fetched via a maximum zoom setting.

I am extremely satisfied with the performance in general, browsing around the map is great. But I am surprised by the performance of the form widget.

If I use the identify tool and drag over a layer to select multiple features, I can see the attributes of the features in the identification result tab in miliseconds. But opening the form, to do actual editing, of a single feature takes +- 5 seconds. Which is a lot, because our data operators will need to open a lot of features.

There is an index on the primary key, there is a spatial index, all foreign key's are indexed, the database is vacuumed....
The only thing that I can think of of doing now is to create an action in QGIS that opens up a custom form.  Oddly enough, I've done a small POC using that workaround, and that seems to be fast. But it seems a lot of work to get the form up just right, and I'd like to limit the amount of code for ease of maintenance purposes. It is slow for all form types (auto generated, drag and drop design or ui file).

Another solution that seemed to work at first, was using a GeoServer + WFS that only fetched the assets in the current screen, but that doesn't work when I try to connect relations, since it does not fetch childeren through the WFS.
(note that with a direct DB connection, it goes slow, even without relations).

Does anyone have any advice?

This is an example ddl of one of the slower tables.

CREATE TABLE belmap.cadastralparcel (
    stableid int8 NOT NULL,
    id int8 NULL,
    versionid int8 NULL,
    validfrom timestamp NULL,
    validto timestamp NULL,
    inuse bool NULL,
    beginlifespan timestamp NULL,
    endlifespan timestamp NULL,
    status varchar NULL,
    nationalreference varchar NULL,
    municipality varchar NULL,
    municipality_stableid int8 NULL,
    geometry geometry(MULTIPOLYGON, 4326) NULL,
    CONSTRAINT cadastralparcel_pkey PRIMARY KEY (stableid)
);

CREATE INDEX idx_cadpar_geom ON belmap.cadastralparcel USING gist (geometry);

This table contains +- 5 million records.

Edit

Screenshots of forms
Basic one (from the table above)

and one of the advanced ones:

The other panels of the advanced form are more or less the name. But contain confidential info.

EDIT 2
I've added the query that is slow here:

SELECT
  st_asbinary("original_geometry",'NDR'),
  "stableid"::text,
  "id"::text,
  "versionid"::text,
  "validfrom"::text,
  "validto"::text,
  boolout("inuse")::text,
  "beginlifespan"::text,
  "endlifespan"::text,
  "type"::text,
  "municipality"::text,
  "municipality_stableid"::text,
  "submunicipality_stableid"::text,
  "road"::text,"subaddressof"::text,
  "status"::text,
  "original_designator"::text,
  "designators"::text,
  "datasources"::text,
  array_out("replace")::text,
  st_asewkt("streetaxis_geometry"),
  st_asewkt("building_geometry"),
  st_asewkt("mailbox_geometry")
FROM
  "belmap"."address"
WHERE (
  (
    ("subaddressof" is null
     and status in ('In use','Not official','Reserved','Suggested')
    )
  )
)
AND ("stableid"::text='136981247')

It is casting the primary key to a text field! Why would it do that? it is clearly disregarding the index!


texts:
As part of a possible solution to  GeoServer WFS Row Level Security? I want create a layer in GeoServer via REST that supplies a SQL statement along with a userid.
I think i would create a feature type and publish it. I need to do this all in REST.
Does anyone have any code examples?


texts: Polymaps + D3+ Polyline with different width and color

texts:
I am actually working on a project which involves polymaps and D3. I was wondering if you could help me on this.

I have 4 locations and I would like to draw a polyline from one location to other 3 location with different widht and color of the polyine. I was wondering if you guys could give me a starting point or maybe a sample code. I think this will be possible by D3.

I also do not understand the locationPoint method, passing in {lat: 0, lon: 0}. or coordinatePoint, passing in {row: 0.5, column: 0.5, zoom: 0},

Also I did try this link http://webcache.googleusercontent.com/search?q=cache:ZJcQP6qOvBwJ:https://github.com/simplegeo/polymaps/issues/91+&cd=1&hl=en&ct=clnk&gl=in
and
http://webcache.googleusercontent.com/search?q=cache:ZJcQP6qOvBwJ:https://github.com/simplegeo/polymaps/issues/91+&cd=1&hl=en&ct=clnk&gl=in but could not understand.

Thanks for your time.


texts: Does study area complexity (number of vertices) affect imagery downloading time in Google Earth Engine?

texts:
Will imagery mapped(clipped) to a polygon ROI consisting of thousands or tens of thousands of vertices be processed and downloaded at the similar rate as a rectangular study area of similar area?
If my study area is a very complex shape, will it be more time-efficient to download a simpler processed image collection and then batch-clip the imagery in GIS software.


texts:
In FUSION you can reproject the final output, for example in case of ASCII to DTM, one can use:
ASCII2DTM canopy.dtm m m 1 10 2 2 canopy_southern.asc

So my question is how can I reproject a LAZ file from NAD83 to WGS_1984_UTM_Zone_11N in FUSION?


texts: Expression returning wrong value in QGIS Field Calculator?

texts:
My QGIS Field Calculator expression:

case when "field_value" = '' or  "field_value" <> NULL then (if( "field_value" = '', 'N', 'Y')) end

returns null value even if "field_value" has something in it.

Note: "field_value" is in string.


texts:
I want to convert some Mapinfo styles to Mapserver.
Here,  I found some examples on how to make symbols. But I don't understand how it works.
For example a dashed line:
SYMBOL
 NAME "mapinfo-pen-73"
 TYPE vector
 FILLED true
 POINTS
  5 25
  25 30
  45 30
  35 25
END # Points
ANCHORPOINT 0.5 0
END # SYMBOL

And i get my dashed line.
What does the POINTS block mean? I can't find any explanation, how to write POINTS by self. Where i can find documentation?
UPDATE
After i'd readed @user30184 and @Andreas answers i'm steal dont understand POINTS block.
I'm looking on my example in map:

i see dashed line.
If POINTS it is a polygons shape it's look like:

This moment dont clear for me. Why shape form different with dashed like form?
Im missed with Y axe. But anyway shape form very strange.


texts: Writing into XMP and EXIF Metadata from attribute table values in PyQGIS

texts:
I've already learnt how to read EXIF and XMP data's from JPG. Now I would like to write from value in attribute table into their XMP and EXIF tags.

This code Import EXIF / XMP , and I would like to do the opposite.

from PIL import Image
from PIL.ExifTags import TAGS, GPSTAGS
from xml.etree import ElementTree as ET

layer = iface.activeLayer()
dpr = layer.dataProvider()

layer.startEditing()
apr_ind = layer.fields().indexFromName("apr")
exp_ind  = layer.fields().indexFromName("exp")
fyd_ind = layer.fields().indexFromName("fyd")

for feature in layer.getFeatures():
    file_path = feature["photo"]
    print(file_path)
    img = Image.open(file_path)
    exif = { TAGS[k]: v for k, v in img._getexif().items() if k in TAGS }
    with open(file_path, "rb") as f:
        s = str(f.read())
    start = s.find('<x:xmpmeta')
    end   = s.find('</x:xmpmeta')
    xmp   = s[start:end+12].replace("\\n","\n")
    tree  = ET.XML(xmp)
    apr = exif["FNumber"][0]/exif["FNumber"][1]
    exp = exif["ExposureBiasValue"][0]
    fyd = tree[0][0].attrib['{http://www.dji.com/drone-dji/1.0/}GimbalPitchDegree']

    dpr.changeAttributeValues({ feature.id():{apr_ind: apr,
                                            exp_ind: exp }})
layer.commitChanges()
print("COMPLETED")


texts: How can I create a radially-symmetric surface from a profile in ArcGIS 10?

texts:
I have created an expected profile for a radially-symmetric feature (an impact crater) using numerical modelling and would like to compare this with the actual topography (a DEM).

My plan is to do this by making a surface/3d polygon based on the profile and then subtracting one from the other. However, I can't figure out a way to create a surface based on the profile. Is there a way to do this in ArcGIS 10?


texts: QGIS 3.10.1: add photos to attribute table using only relative path

texts:
I want to create a shapefile from geotagged photos using only relative path (to display photos from any other GIS software alongside other informations in the shp).
I imported photos using Import Photo plugin and exported them to shp, but when i open the shp from another machine, photos don't show up and an error message appears "No image Path found"
Is there a solution to make it work?


texts: In FME, is there a way to connect multiple input feature types from reader and route to single transformer?

texts:
In FME, is there a way to connect multiple readers to a single transformer? I'm using FME 2015 and want to clip a bunch of feature classes to a single feature, but can't seem to find an easy way to do this.

EDIT: Trying to select multiple features types from input dataset and link them all to a single transformer without having to click individual output arrows on the feature types and link to the input arrows on the transformer.

EDIT 2: Here's a screenshot, hopefully this will help: