texts: How to pass and replace a parameter to a TileServer GL styling JSON-file

texts:
I know it is a very specific question, but I am wondering if it's possible to pass a (key-)parameter like ...?key=1234 to a styling file which will be consumed by the TileServer GL? My styling file looks similar to the following snippet:

{
  "version": 8,
  "name": "basic",
  "metadata": {
    "mapbox:type": "template"
  },
  "sources": {
    "openmaptiles": {
      "type": "vector",
      "url": "http://localhost:8100/data/v3.json?key={key}"
    }
  },
"glyphs": "http://localhost:8100/fonts/{fontstack}/{range}.pbf?key={key}",
"layers": [

...

I want Tileserver GL to actually replace the {key} parameter with 1234 as mentioned above and use the new URL for further requests (further requests shall lead to: http://localhost:8100/data/v3.json?key=1234).

Currently, when I just pass a parameter in the original request (e.g.: http://localhost:8100/styles/basic/style.json?key=1234), Tileserver GL just ignores the parameter even though the server needs the JSON file and its included URL's for the upcoming requests.

Code
Doc

Thanks


texts: What is the most efficient way to create millions of lines using Python from coordinates?

texts:
I am trying to create a flow map and I have 150 million rows of data. I have tried to use QGIS's Virtual Layer (described in Creating virtual layers using Python and export them without using QGIS UI) to create the lines, but it is taking forever (I have an 8GB RAM). How can I create a GeoJSON of the lines as efficiently as possible?
My data looks like this:

TripID
x1
y1
x2
y2

1
103.62183
1.27624
103.62181
1.29713

2
103.62542
1.39132
103.99923
1.37252

3
103.99156
1.27613
104.00284
1.37564

I need to create a line for every "tripID" using the two coordinates, then export as a GeoJSON.


texts: Using wild card selection for value while running feature class iterate in ModelBuilder

texts: I am running iterate in ModelBuilder with an input feature class that has a similar field in another feature class. The basic principal of what I am trying to do seems to work but I cannot figure out how to implement the wild card functionality in the select by attribute while using the Value produced by the iterate function. Below is where I am currently. It's two different processes, the first iterates through a field in one feature class and the other selects by the Value of that iterate but the fields in the second feature class are more complex than the first so I am trying to use FieldName Like '%Value%' to pull all of the features in the second feature class that have at least one instance of the Value in their data. The select function however doesn't seem to recognize the wild card and only pulls the fields that have an exact match for the Value.


texts: is there any relationship between drawing graphics and IIS overloading

texts:
I'm running gis web application built on silverlight but when many users connect to the system the IIS reaches the full memory so could the reason be from many users draw graphics at the same time ,if not what could makes the IIS overloaded and not responding.

regrads.


texts: Move polygon label on geoserver layer by css style

texts:
I wish to change label position of polygon in Geoserver layer (custom local osm map), like that:

but I don't understand how.  label-anchor: was ignored.
my style:
    * {
  font-family: "DejaVu Sans";
  font-size: 8px;
  halo-color: #ffffff;
  halo-radius: 0.4;
  font-weight: normal;
  font-fill: #414142;
  font-style: normal;
}

[@scale < 35000] [@scale > 5000] {
  fill: #ded3be;
  stroke:#C4BBA9;
  stroke-width: 0.5;
  fill-opacity: 0.4, 1;
  geometry: [offset(way, 0.5, -0.5)], [way];
}

[@scale < 5001] [@scale > 2500] {
  fill: #ded3be;
  fill-opacity: 0.4, 1;
  geometry: [offset(way, 0.4, -0.4)], [way];
  stroke: #C4BBA9;
  stroke-width: 0.5;
  font-size: 12;
  font-fill: #414142;
  font-weight: 600;
  label: ["housenumber"];
  label-anchor:  1.0 0.5;
}

[@scale < 2501] {
  fill: #ded3be;
  fill-opacity: 0.4, 1;
  stroke: #C4BBA9;
  stroke-width: 0.5;
  font-size: 12;
  font-fill:#414142;
  label: ["housenumber"];
  label-anchor: 0.5 0.0;
  font-weight: 600;
 }

how to move label to the new label position?


texts: How to permanently rearrange the order of features in a Shapefile with QGIS?

texts:
I have created a trimmed/clipped map of Europe, the Middle East and North Africa in QGIS 1.7.3 by deleting all the other countries from a world map shapefile. I would like to permanently rearrange the countries that are left in the shapefile in a particular sequence in the attribute table. I searched for other similar questions on this website and one answer suggested that the .dbf file be modified using a spreadsheet program, e.g., MS Excel 2003. I tried doing that and it seemed to work; but then I realized that when I pointed at Austria in the shapefile, the name that came up was Egypt. Hence, on the map itself, the names of the countries are all mixed up now.

Can you please suggest a way to arrange the countries in the attribute table according to a specific sequence without its affecting the countries on the map itself?


texts: Can I download Google Map Orthophotos using QGIS?

texts:
Does QGIS have functionality to download Google Map Orthophotos? Or any other open source GIS?

I specifically want to be able to not only view Google Maps on my GIS, but also be able to download the imagery so that I can use it to orthorectify (means: provide georeferences to) other historic digital aerial photos, which is done in a different program.

(I'm inexperienced using QGIS)


texts: Select polyline features and see total length in ArcMap

texts:
I have a layer of polylines representing sewer mains. I want to break the layer down into groups of approximately the same length and area.
Is there a way I can select multiple features in ArcMap and see the total length.  That way I can add or subtract until I get the selection approximately correct.  There are over 3000 lines covering 127km, and I want to break it down into say 10 packages of approximate same total length and in the same general area.  I have a range of similar assets that I want to repeat the process on.
I have seen the summary statistics tool that can do this for the full layer.  However there doesn't seem to be an option to do it for the selection only like QGIS often has.


texts: Combine files in TopoJSON when one is simplified differently to the other

texts:
I have two .shp files, I am separately  converting these to GeoJSON with ogr2ogr.exe and then combining the files with TopoJSON.
Each file independently becomes a ~200kb topojson file, however the combined file is over 2mb.

It appears that one of the shapefiles has already been simplified slightly and so topojson is not finding a good correlation between the coastline arcs.

Other than trying to get a unsimplified source file, is there anything I can do to help topojson combine these files better?


texts: How can I count points from several buffer distances in PostGIS updating the table?

texts:
I have a problem regarding a calulation on PostGIS. I have to calculate how many facilities fall within a given distance. This is not really a problem though, but updating my table is.

I have this code:

CREATE TABLE my_table AS
SELECT buffer1.id, COUNT(facility1) as bf1
FROM buffer1 LEFT JOIN facility1
ON ST_CONTAINS(buffer1.st_buffer,facility.geom)
GROUP BY buffer1.id

Now I have to do this again (35 times), so I thought to update the table by entering:

UPDATE my_table SET buffer2 FROM(
    SELECT /*the count statement*/) as buffer2

But this won't work.

I managed to get a extra column with the Alter Table statement, but it was empty.

I'm able to this manually, but it is easier to this in 1 query, instead of rewriting 36 queries.


texts: Adding separator to ArcGIS python add-in menu?

texts:
In several places (e.g. Making separator and button menu in customized Toolbar using ArcObjects?) I have seen it mentioned that the way to add a separator for add-in menus and toolbars is adding the attribute separator="true" to the button in the config.xml file like: <Button refID="test105_addin.button_1" separator="true"/>.

However, I can only get it to work with ArcGIS 10.3 and not with 10.1/10.2.2.

Does anybody know how a separator is added using ArcGIS 10.1 (or 10.2.2)?


texts: Why is tiff different in epsg 3857 and epsg 900913 in gdal

texts:
When I transform tiff from 4326 to 3857 and from 4326 to 900913 in gdal, I get different result(when display in Qgis) eventhough these two projections(3857 and 900913) should be the same. What is problem for that?

Image shows differance. Brighter(white) layer is from epsg 4326 to 3857 and darker(black) layer is from epsg 4326 to 900913. I know it's not a real standard, but in gdal gdalsrsinfo -e epsg:900913 and gdalsrsinfo -e epsg:3857 are same,I mean both have same proj.4. So I don't understand, why are layers differently displayed in Qgis.


texts: Which is the best way to store roads for the shortest path finding application?

texts:
I would like my application which is simply a shortest path finder in a city to meet the following requirements:

the city roads network should be stored persistently on disk and loaded into memory on demand (to execute the shortest path finding algorithm);
I would like to use something like OpenLayers to view the results that have been yielded by the shortest path algorithm;
the 'OpenLayers view' should contain not only the results, but the map (uninteractive) of the city itself.

As for now I am thinking of the following architecture:

PostGIS enabled database to store roads ('roads' table);
PostGIS enabled database to store the actual map data to render it in OpenLayers;
a table for the results of shortest path algorithm in the PostGIS enabled database (I would like to view a couple of shortest paths simultaneously);
a web app with OpenLayers which shows the map from the 'map' table, shows the roads layer from the 'road' table and the shortest path from the 'shortest path' table.

Does this sound correct or should I change anything?


texts: How to remove Commas from field values in multiple files?

texts:
I'm working in QGIS, presenting data in ArcMap, and Arc is having trouble with the attribute tables.

I have a "species" column which lists the species found in that area and have added commas "," in this field to separate species. Arc can't handle the commas so I would like to replace these in QGIS simply with a space.

Is there an easy way to do this?

I have tried replace and regexp_replace but I'm either doing this wrong or it is just not the way to do it...

I have a large amount of files to work with so an easy fix to this would be welcomed!


texts: How to get Lat/Long between start and end coordinates and spaced at equal distance?

texts:
Suppose I have two GPS points (start and end points) each with lat/lon. I want to get equally spaced GPS points between start and end nodes. That is, I want to find GPS points  for example every 10 meter from start node to end node.

Similar question is asked here.
From the answers given I have referred the one by Vikram Gupta. This solution works fine but if you replace

int dist = (int) d / interval

by

int dist = (int) d

in getLocations(...) method some calculated coordinates lay beyond the end coordinate. How can the code be modified to prevent this?


texts:
I am using jsf2leaf library to add a map to a web application using primefaces. The map and the markers(popup message when clicked) are displayed correctly, however I want to display a primefaces dialog when a marker is clicked instead of the popup message.

any idea of how can I add a click listener to the markers!?


texts: QGIS export circle to KML then import to Google Map

texts:
I am  using QGIS version 3.14.0-Pi.
In QGIS, I display Google Map, with the XYZ tiles. In QGIS, on a Project using EPSG:4326, when I draw a circle, this circle is a real circle (not an ellipse). Then, I export the layer containing this circle into KML. In Google Map ("MyMaps"), I import the KML. Google Map displays this circle as an ellipse; not a circle.
I tried to export the QGIS layer to KML with different projection systems:
EPSG:3857   -- WGS84 Pseudo-Mercator
EPSG:4326   -- WGS84
EPSG:900913 -- Google Maps Global Mercator

But Google Maps is always displaying an ellipse, not a perfect circle.
Does anybody know what is the way to export from QGIS a perfect circle and import it and display it into Google Map as a perfect circle?


texts: Save processing modeler output layers to project folder/relative path?

texts:
I am building an algorithm in the processing modeler which takes a csv file as input and creates a number of vector layers as output. I would like to automatically save these layers in the same folder as the project file. As I understand, it should be possible by entering the path "./filename.shp".

I can set the default save path for each layer in the final "Object Sink Parameter-definition" box in the modeler, but can't get this to work with a relative path.

When setting absolute path, everything works:
C:/Users/Rasmus/Desktop/testfilQgis/målepunkter.shp

When setting only filename (felt.shp), the filename is added to the modelers default output directory path and saved to this path. When running the model on other computers, it will look for the local directory on my machine.

If I use a relative path (./koter.shp) that whole path will be added to the default output directory path (and break since the . folder doesn't exist).

Only the one folder up path seems to be accepted, but returns a "Failed to create: Permission denied".

See below image for the suggested save paths based on these inputs:

I've looked into using the project_folder as suggested here, but the "object sink parameter-definition" box only takes a path, so I'm not sure how to do that.

And setting the save path manually for each layer every time seems like defying the purpose of building a model. Surely, it must be possible to suggest a relative save path?

project > properties > save paths is set to relative (default).

Thanks in advance!


texts: Use Esri PersonalGeodatabase layers in OpenLayers

texts: I would like to know if it is possible to open a Personal Geodatabase of Esri with OpenLayers with FileGDB_API_VS2010_1_2


texts: gdalinfo image.tif>stats.txt - returns empty text file (Windows 7 64 bit)

texts:
I am trying to output the gdalinfo image statistics results into a text file with the following command:

gdalinfo image.tif>stats.txt

The stats.txt file gets created, however it is empty.
It does not contain any statistical information that is printed out in the command window, with the gdalinfo image.tif.


texts:
I am in the process of building a script for my degree capstone project with the use of arcpy and the Planet API module.

Currently it takes a GeoJSON created from a simple spatial join and in turn downloads the required area before cropping and then segmented and classified through the image analysis software eCognition (this last part is not scripted and just uses the Cognitive Neural Network we have developed.

All my parameters are working as expected but I am having trouble with the coverage of the area. I cannot seem to find any information on how to search for 100% of a required area. Some results come back unaffected but others either have partial or totally missing scenes. Below is the script I am using.

#set GEOJSON location
GEOJ = "D:\\Planet_Download\\GeoJSON\\TWP0630804.json"

#open GEOJSON file
with open(GEOJ) as f:
    data = json.loads(f.read())
for feature in data['features']:
    poly_type = feature['geometry']['type']
    coords = feature['geometry']['coordinates']
    geojson_geometry=feature['geometry']
geometry_filter = {
    "type": "GeometryFilter",
    "field_name": "geometry",
    "config": geojson_geometry
}

from planet.api import filters
from datetime import datetime

#set date variables. June to July = most green months
start_date = datetime(year=2018, month=6, day=1)
end_date = datetime(year=2018, month=7, day=31)

date_filter = filters.date_range("acquired", gte=start_date, lte=end_date)
#cloud cover less than 1%
cloud_filter = filters.range_filter("cloud_cover", lte=0.01)
#set requirments togther
and_filter = filters.and_filter(date_filter, cloud_filter, geometry_filter)

#add required platform
item_types = ["REOrthoTile"]
req = filters.build_search_request(and_filter, item_types)
print req
print ""

#open quick search
res = client.quick_search(req)

I am relatively new to Python and scripting so I am not sure if there is an alternative method I could use to receive the total required area or if there is a known process to achieve 100% coverage through the Planet API.


texts: I have been developing a python plugin that uses Processing algorithm for raster calculations. I want to support the gui with progress bar. How can i connect Processing job to progress bar? I found little description in here but no success.


texts: Getting attribute data from one dataset to another

texts:
In the image below you can see two shapefiles: NYC building footprints (outlined), and underneath a PLUTO shapefile rendered with a categorized blue gradient using the "Floor Number" column. The Building Footprint dataset does not have floor number info and the PLUTO dataset doesn't have buildings, just lots.

How can I use QGIS to get the number of floors data from the PLUTO shapefile into the Building Footprints table without having to do it manually, building by building?


texts: Illustrator to QGIS workflow for vectors

texts:
I vectorized a map in Illustrator and would like to import it to QGIS now. I know I must use DXF for that purpose but when I convert the DXF to SHP with the QGIS plugin it is not georeferenced which is actually not surprising. But the QGIS georeferencing tool works for raster images only. So what to do?

The base TIFF image I used as the bottom layer in Illu is georeferenced, so maybe I can use the settings from there but I didn't figure out how. I heard you can use the .prj file from another layer but my TIFF doesn't have this, it has a .points only.


texts: Having spaces in QGIS leaflet web map layer names?

texts: When I export a QGIS project to a leaflet webmap using QGIS2Leaf, is there a way I can have my layer names display with spaces. So for example in the image I would like OffshoreSACS to display as Offshore SACs and LicensedblocksMay2015 as Licensed blocks May 2015 etc.


texts: How to query nearest date in ArcGIS Server REST?

texts:
How to query nearest date (today) in ArcGIS Server REST?

Table:

ID    NAME            DATA
1     test           2013-06-01 13:00:00
2     test1          2013-06-01 16:00:00
3     test2          2013-07-01 18:00:00

Expected result : query  ID = 3


texts: Using arcpy to reclassify "NoData" to a value

texts:
There is a part in my Python script where I convert a Feature to a Raster (arcpy.FeatureToRaster_conversion). The output is a raster with a certain value (based on a field of the Feature). The rest of the extent of the resulting raster has "NoData" values.

However, I want that "NoData" to be 0. Here is the code.

arcpy.env.overwriteOutput = True

raster = arcpy.Raster(export_path)
arcpy.sa.Reclassify(raster, "Value", arcpy.sa.RemapValue([["NoData", 0]]), "").save(export_path)

But this has no effect. I still have "NoData".

How do you reclassify "NoData" values to a numerical value using arcpy?


texts: Determining direction of polygons using QGIS?

texts:
Is it possible to get the direction of the longest side of a rectangle in QGIS?

I know how to get the direction of a line layer. But I need to know the direction of a polygon. One could start at the center of a polygon and find the shortest way to the border by using a line and take the direction of this line, but I don't know how to draw this line automatically for 30000 polygons.


texts: Comparison between two fields of two tables using ArcPy cursors?

texts:
I would like to make a comparison between two fields of two different tables with a Python script.
Fields for comparison are ROAD_NAME and NAME.ROAD.
The idea is to check whether each element of NAME.ROAD field exists in the field ROAD_NAME.
Whether it is necessary to check if it has the right DISTRICT, CITY and good in case there is a mistake we must rectify at the tab1
DATA SOURCE : TAB1 --> FGDB TAB2--> Excel

I'm using ArcGIS 10.2 Advanced.

the code that I can write up to now is the following:

cursor1=arcpy.da.SearchCursor("TAB1",("ROAD_NAME","DISTRICT","CITY"))
cursor2=arcpy.da.SearchCursor("TAB2",("NAME.ROAD","DISTRICT","CITY"))
stockage={}
for row in cursor1:
    stockage[row[0]]=row[0]
for row in cursor2:
    if row[0]==stockage[row[0]]:
        pass
    else:
         ....


texts: Transform Gauss Krüger Zone 3 data into WGS 84 / UTM 33N in QGIS

texts:
I would like to transform this data:
https://opendata.dwd.de/climate_environment/CDC/grids_germany/halfyear/precipitation/ (it is Gauss Krüger Zone 3) into WGS 84 / UTM 33N in QGIS.

How do I do it? Just setting to UTM didn't seem to work.


texts: Sample JSON code to get latitude and longitude values from Google API geolocation

texts:
I'm a fresher in JSON coding.
I got a sample Google API geolocation console from https://developers.google.com/maps/documentation/business/geolocation/ site.

{'cellTowers': [
  {
    'cellId': 21532831,
    'locationAreaCode': 2862,
    'mobileCountryCode': 214,
    'mobileNetworkCode': 7
  }
]}

But, i'm getting below errors.

{
 "error": {
  "errors": [
   {
    "domain": "global",
    "reason": "parseError",
    "message": "Parse Error",
   }
  ],
  "code": 400,
  "message": "Parse Error"
 }
}

Please provide me sample JSON code to get latitude and longitude values from Google API geolocation.


texts: Exporting 2 shp layers to geojson using QGIS?

texts:
I will explain with an example:

I have 2 shp files for Tanzania.

Tanzania inland waters http://www.diva-gis.org/data/wat/TZA_wat.zip

Tanzania administrative boundaries   http://gadm.org/data/shp/TZA_adm.zip

I open both shp files (TZA_wat.shp and TZA_adm0.shp) in Quantum GIS Desktop.
I want to merge the layers : TZA_wat and TZA_adm0 and export/save both layers to one geojson file.

While I can do this for individual layers but not able to do this for both at once. I have written code to process geojson files but I get stuck at this step.

Is there a way to first merge both these .shp files to a single .shp file and then export/save as geojson?

I used the Union option in QGIS and merged the two shp files (here) into one. I uploaded the file here in a zip if anyone wants to take a look. The zip has shp files and geojson file in it.

input shp files : link
output shp file: link


texts: Openlayers 3 ol.extent to ol.geom.Polygon conversion

texts: In openlayers 2 there was a method "toGeometry()" on Openlayers.Bounds() that could convert that bounds into a Geometry. I migrated from ol2 to ol3 and I try to find a solution to easily do the same thing without having to define a ring with all the coordinates and then using the constructor of ol.geom.Polygon(ring). The API is not helping.. Any ideas?


texts: Does anyone know how to make OSM national_park style using GeoServer GeoCSS?

texts:
Does anyone know how to make OSM national_park style using GeoServer's GeoCSS?
At the moment I've succeeded only using stroke-geometry: [buffer(geom, -10)], [geom];. But the result is not satisfying.
This is OSM National Park boundary style (below).

This is my style (below).

As you can see, when I zoom in, distance between original and buffered geometry increases.
Is there any better way? How to make that OSM style in proper way?


texts: Problems installing QGIS 2.14.14 and QGIS 2.18.6 (32-bit) on Windows Vista 32-bit AND QGIS 2.14 and QGIS 2.18.14 64-bit Windows 7 64-bit computer

texts:
I used to have QGIS installed on my computer (functioning) but I decided recently to upgrade.  However, attempts to install QGIS have failed.  My computer is running Windows Vista (upgrading my operating system is not an option in my case) with 32-bit operating system.

During installation (which I have tried with both QGIS 2.18 and QGIS 2.14 several times now), the following error shows up:
"Could not find C:/Program Files/QGIS 2.18/postinstall.log" but I can find a postinstall.txt file when I look manually

When I try to open QGIS, I get the following error message:
"This application failed to start because api.ms.wind.crt.runtime-l1-1-0.dll was not found. Re-installing the application may fix this problem."

If there is a known problem with the latest QGIS and older operating systems, what is the last version of QGIS that is known to work correctly with Windows Vista? [edit: I have successfully installed QGIS 2.10.1 on Windows Vista]

Additional issue: QGIS 2.18.14 (latest release in December 2017) nor QGIS 2.14 (long-term repository release) nor OSGeo4W does not install on my Windows 7 64-bit computer. No error messages were given during the Stand-alone installation but the same error message when I attempt to open QGIS pops up:

"This application failed to start because api.ms.wind.crt.runtime-l1-1-0.dll was not found. Re-installing the application may fix this problem."


texts: Is there a way to write a GDAL raster to a HDF4 or HDF5 dataset using GDAL? I have found many examples that show how to read HDF datasets, but none on how to write them.


texts: Count the number of corners of a polygon

texts: How can I automatically count the number of corners that has a polygon using VBA and ArcObjects?


texts: Projection of Shapefiles for Spain (INE)

texts:
I downloaded the official shapefiles for the municipalities of Spain from INE here.

It has a lot of shapefiles, but I found out that the right one for what I need is esp_muni_0109.shp.  I would like to convert these shapefiles to WGS84 with ogr2ogr.  The problem is that these shapefiles have no .PRJ files, so I don't know from what projection system to convert them from.

Does anyone know the right projection system for these shapefiles?
If so, can anyone tell me what parameters to use with ogr2ogr?  (I'm pretty new in all of this).

Thank you.


texts: How to add GDAL to existing Python installation?

texts:
I would like to download and install GDAL for Python version 2.6 [MCS v.1500 32 bit(Intel)] on win32 that supports reading FileGDB's with Python.  But I am not sure how to go about this.

I found one exe at this link but it was not built with the FileGDB included.  I also tried the OSGeo4w Installer but it currently installs its own 2.7 version of Python rather than allowing the GDAL to be installed with the existing version of Python.

I must work with the 2.6 version installed on our systems.  Can anyone be of assistance in locating a GDAL build that supports FileGDB's, can be installed with my existing Python version 2.6 and provide detailed instructions?


texts: arcpy.sa.SegmentMeanShift() Creates segments with values outside of the range of original raster (stretched?)

texts:
I'm using arcpy.sa.SegmentMeanShift() to produce a segmented raster in ArcGIS 10.7.1.

The segments in the output raster have values outside the range of values (1 to 248) in the input raster (-22 to 92)!

The new range of values in looks like a fit for stretched values (0-255).  But I need the original values, not the stretched values.

How can I get the (mean) original values in my segmented raster?

I have tried arcpy.SetRasterProperties_management(inputRaster, data_type="PROCESSED"), but this has not helped.

UPDATE:  ESRI Rubbish...?

Turns out the problem doesn't always happen.  I've changed my output to be in a geodatabase, and the problem went away!  I haven't tested conclusively if this is why the problem went away... to busy moving on to other things.


texts: Where can I find decent resolution (5m ish) night time orthophotography/satellite imagery data for Ireland?

texts: Anything I've found so far has been incredibly granular. I need a resolution which is fine enough to see general lighting footprints for motorway junctions.


texts:
QGIS 3.2.3: I have a shapefile and an Excel table which is joined to the shapefile based on matching IDs.

I would like to filter the shapefile while using a joined field from the Excel table.

When I select "Filter" at the shapefile, I get the message that filter on joined tables is only possible when creating a "Virtual Layer". I click "Yes" when asked if I want to create a virtual layer.

Nevertheless, when the query builder appears, the joined fields do not appear. Did the creation of the "virtual layer" not work?

Is this not possible or am I doing something wrong?

I have found this older article: is this still valid that it is not possible in QGIS? How do you filter a layer on a joined table in QGIS?

This works in ArcGIS using a "definition query" on a joined layer.


texts: Leaflet - Draw vector layer behind tile layer

texts:
I'm using Leaflet to create a map in which I need to render a tile layer on top of a GeoJSON layer:
L.geoJson(...).addTo(map);
L.tileLayer(...).addTo(map);

The Tile Layer should be render on top.

I have use setZIndex and brinToBack/bringToFront methods without success. I have read that GeoJSON ordering in Leaflet relies on rendering SVG and maybe it's not possible to do that.
So the question is: It's possible to render a GeoJSON layer behind a TileLayer?


texts: How to get the name of location in Google Maps?

texts: I am new in google maps api. I want to get location name after user click to that location. I try to find answers in docs, but it so huge. Maybe someone can help to me or just say where in docs I can read about this?


texts: Printing out list of layers and associated information

texts:
Under Layer/Layer Properties/Information in the menu, there is a table named "Information from Provider" which gives information about the layer such as Name, Path, Geometry, CRS, etc. Is there a menu item somewhere or a plug-in that will print out a list of this information for all layers in a project?

I am using QGIS version LTR 3.10.4.


texts:
I'm performing a simple select by location. I want to select polygons from TQ_Building that fall within rapidair_AOI, a larger polygon. My workspace looks like this:

I'm using QGIS 2.18.13. I've performed this through the Vector > Research Tools > Select by location menu, and through the python console (see code below) but no selection is made and a dialog box with 'unknown See log for more details' appears.

import processing
mapcanvas = iface.mapCanvas()
layers = mapcanvas.layers()
processing.runalg('qgis:selectbylocation',layers[1],layers[0],u'within',0,0)

The layers share the same CRS, they are enetered correctly as input and intersect layers in the code. Any help / suggestions?

UPDATE: following @GeoDatup pointing me to log, here's the log message:

  Uncaught error while executing algorithm          Traceback (most recent
  call last):
              File  "C:/PROGRA~1/QGIS2~1.18/apps/qgis/./python/plugins\processing\core\GeoAlgorithm.py",
  line 203, in execute
                self.processAlgorithm(progress)
              File "C:/PROGRA~1/QGIS2~1.18/apps/qgis/./python/plugins\processing\algs\qgis\SelectByLocation.py",
  line 93, in processAlgorithm
                index = vector.spatialindex(inputLayer)
              File "C:/PROGRA~1/QGIS2~1.18/apps/qgis/./python/plugins\processing\tools\vector.py",
  line 219, in spatialindex
                idx = QgsSpatialIndex(layer.getFeatures(request))           Exception: unknown

UPDATE 2017-11-03

The above was running QGIS 2.18.13 on Windows 2012 r2 server. I have also tried to replicate on my local Dell machine with the following results:

But there's no selection made in the mapviewer (which is what we'd expect given the second 0 parameter in the selectbylocation algorithm) and no extra polygons in the rapidair_AOI layer (used for the intersection enyway).


texts: How to set tolerance value of ST_Intersects query in TileMill?

texts:
I want to place river names in Tilemill. I use "text-dy" to displace label from the center line of a river. But for labeling polygonal rivers, names extracted from lines are used too. So at a small zoom, labels of polygon rivers crawl over river's banks (see attachment_1).

I used ST_intersects to exclude from the selection linear rivers being intersected by polygonal ones, and not to displace labels for them:

( SELECT way,name,
  FROM planet_osm_line
  WHERE waterway IN ('river','canal') AND name IS NOT NULL
  AND ST_Intersects(way,(
    SELECT ST_Collect(way)
    FROM planet_osm_polygon
    WHERE waterway = 'riverbank' AND way && !bbox!))
) AS riv_center

It runs sucessfully, but there are some other problem: line rivers, falling into polygon ones, have short intersection with them and are labelled too (see attachment_2). Maybe one can specify some tolerance for intersection to exclude all line rivers from selection?

Thanks in advance.


texts: Is it possible to create Domains for an Oracle (non SDE) database?

texts:
I have a simple Oracle database (non SDE/Geodatabase) but I want to constrain the values that can be entered into a particular field (ex: the state names that can be stored/entered in a "State" field)

This is exactly what Domains are, in an ESRI geodatabase. Does something like this exist outside of Geodatabases? If so, how can I do it?


texts: Batch open GeoJSON LineString from files with multiple features

texts: I want to open multiple GeoJSON files, each have Points and LineStrings - I only want the LineStrings.  I want to open potentially 1000s of files so don't want the popup asking if I want Points or LineStrings each time.  Is there a tool I can use?


texts: Table join or virtual layer to match multiple points to geometry

texts:
I have:

a non-spatial table of multiple test results with reference to non-unique locations, and
a spatial layer that has all the locations (multiple results per location)

I can do table joins between the locations layer to the results, but I get stuck because each location has multiple results...

I suspect I need to make a virtual layer assigning the geometry to each result, but I get lost at the "query" stage. The QGIS manual goes deep into database SQL queries but I am hoping there is a simpler solution for the non-initiated?

Note: I can create virtual "x" and "y" attribute columns and join them to my table, but still missing a step to use the coordinates in creating the virtual layer


texts: Frontal Detection of sst & Chlorophyll image

texts:
I am working on sea surface temperature satellite images and looking for delineating gradients in open source on the basis of temperature. Is there any tool like Cayula and Cornillon front detection tool of ArcGIS [1] that may help in this.

The one method that is suggested to me is to use high pass filter and contouring. But in both the cases I am not able to draw a line on a temperature specific front. In high pass filter it is detecting all the edge which I am not interested in and in contours it is interval based.

I'm searching for a  tool in QGIS (or any open source) that automatically:

identifies gradients(fronts) based on temperature,
draw the line and
generates shapefile.

Per a commentor's request: the Cayula Cornillon front detection tool automatically detects SST fronts from satellite SST data. The original paper is published here:

  Cayula, J. F., & Cornillon, P. (1992). Edge detection algorithm for
   SST images. Journal of Atmospheric and Oceanic Technology, 9(1), 67-80.

Here's an example dataset: https://www.dropbox.com/s/rtxhln299szonb1/may25.csv?dl=0


texts: Get raster value at specific point (arcpy)

texts:
I have an arcpy script creating thousands of points in a feature layer, based on a river network. At creation time, I'd like to use the value from a (large) elevation raster corresponding to the new point's location.

It needs to be at creation time, because I want to subtract the value from the previous value, and that information will be difficult to get afterwards.

I've tried using sa.ExtractValueToPoints, but it's very slow and it's a lot of extra overhead to create a point feature with my x,y data, run the function, and then extract the value from the newly created feature. There are a few tens of thousands of points...

Is there some way of getting directly at the raster data (while retaining the spatial information) that I don't know about?

Edit: I managed to work around it as explained below in my comments. I'll explore the NumPy option, as I think this is something I might run into again.


texts: How convert .tif to .mbtiles keeping projection using gdal

texts:
I trying convert .tif using:

gdal_translate 4244.tif 4244.mbtiles -of MBTILES
gdaladdo -r average 4244.mbtiles 2 4 8 16

But, o file .mbtiles appear with diferent projection from origin file. Always convert to WG84?

Help me please!!!


texts: geometries with different SRIDs errors (even after setting SRID)

texts:
I've looked at existing questions about this (including this one) and followed the advice that was shared there. I'm still seeing this when I run some queries:

RuntimeError: Postgis Plugin: PSQL error:
ERROR:  Operation on two geometries with different SRIDs
Full sql was: 'SELECT ST_AsBinary("wkb_geometry") AS geom,"ogc_fid" FROM (SELECT ogc_fid,wkb_geometry,ele::integer FROM srtm3_contours_50m) AS "srtm-contour-50m" WHERE "wkb_geometry" && ST_SetSRID('BOX3D(-549661.6930924922 7480280.288652307,-338719.9409595227 7801966.460655085)'::box3d, 900916)'

Here's what I did to try to set the SRID:

UPDATE srtm3_contours_50m SET wkb_geometry = ST_SetSRID(wkb_geometry, 3785);
UPDATE srtm3_contours_50m SET wkb_geometry = ST_Transform(wkb_geometry, 3785);
VACUUM ANALYZE srtm3_contours_50m;

How can this still be happening?! Is my query to blame?


texts: Adding up polygon areas involved in a single area in QGIS

texts:
I have a polygon shapefile comprised of 4 polygons representing landslides (see image). None of them are overlapped.
Number 1 comprises everyone (2, 3 and 4). Number 4 comprises number 3. Finally number 2 and 3 are indivual.

I've calculated the area of every individual polygon but I need to calculate the involved areas as mentioned above.
Is there any possibility of adding up the areas creating a new field using the field calculator?


texts: Divide a polygon in to smaller polygons using JavaScript


texts: Fastest way to transfer data from SQL to ArcMap?

texts: I frequently have to use large tables and results from MS SQL in ArcMap and for that saving the data as .xls is one of the options. But when I add the data to arcmap, processes like joining and plotting the data take a lot of time. My question is that what would be the fastest method of importing SQL data into Arcmap and then working on it?


texts:
I am working with a DEM downloaded from the national park service. The peaks appear like completely white areas, I assumed it was just a subtle change in elevation.

I exported a high resolution (4096 x 4096) region to .png format and brought it into world machine. The resulting model has completely flat peaks.

Could the DEM be bad? Do I need to export to a difference format? Is there something I've missed?


texts:
I work for a government utility company. Our maps are broken down into 500x500 square meter grids. We have scanned images (TIFF) of these maps that we need to overlay and digize. The maps and the scans are saved in different folders with corresponding map and scan being named for example 5049.dwg and 5049.Tiff respectively. In the past we have been Georeferencing in AutoCAD but we migrated to GIS for better accuracy. Each of these scans are the same size but have different coordinates as they are building blocks of a bigger map,  like a jigsaw of sorts. Is there a way to automate the batch Georeferencing of these scans using either python or model builder?.
The question was marked a duplicate, in this case the coordinates of each image is different.


texts: ArcMap - automatically select record based on data driven pages value

texts:
I'm designing a workflow involving 1km grid cells which are referenced by data driven pages. In it the user will use the data driven pages toolbar to skip to the next cell in the index layer and perform a manual interpretation task for that record. This will simply involve clicking on buttons set up in a custom toolbar. So far so good.

What I'd like to do is for ArcMap to automatically select the record in the index layer to which the data driven pages value relates. I have a hunch that this is possible but I can't find a way to make it happen.


texts: Is it possible to reduce the amount of points of polygons from overpass-turbo export?

texts:
I have this query:

[out:json][timeout:25];
{{geocodeArea:Latvia}}->.searchArea;
// gather results
(
  // query part for: “admin_level=6”
  relation["admin_level"="6"](area.searchArea);
);
// print results
out meta;
way(r);
out body;
>;
out skel qt;

http://overpass-turbo.eu/s/cds

Size of output file is near 20mb. I want to reduce size by reducing the amount of points. Is it possible to achieve it?


texts: Keep column titles when exporting SQL query in MapInfo

texts:
Is there a way to keep column titles when exporting a query in MapInfo as a csv or a new tab. This is my result:

And my export looks like this:

It's not a major issue, I jus like saving myself time wherever I can.


texts:
I am using QGIS. I have three different study sites on an island, and have defined these as three polygons using 'clip'. I want to export the geometry data so I could import it into 'Distance' for study design - however, the polygons have an attribute table that is basically empty. Therefore when I try and export, I get .csv file which only contains the column headings "shapeid" "x" and   "y" but no data.

I'm not sure what I have done wrong. The original shapefile that I used as base for clipping has a full attribute table with many points.


texts: Turning one column to multiple columns in Attribute Table?

texts:
This is kind of a two-part question to see which route would be more efficient.

I'm trying to display a label in a shortened version, but the info is all in one column:

Township 1 South Range 20 West

to

T01S-R20W

I tried substr() expression, but since the numbering (for both numbers in the string) vary from 1-29 I cannot use a long substr() expression that joins them together. It moves the start/stop position.

So, I was thinking to somehow do what MS Excel does with "text to columns", but in QGIS. Is there any way to do that? Or, is there an easier way to do what I'm looking for?


texts: Coloring polygons according to total sum of attributes within polygon?

texts:
I have point data about wind turbines across a country with information about the power output. Second shapefile I have is the country state borders (polygon). I now want to color the state polygons according to the total power of the turbines in the polygon.

How can I do that with QGIS?


texts:
I have a simple app that takes a coordinate and determines what GeoJSON polygons it falls inside, then I take the polygon data and put it into a table. I have working code that uses turfjs and Leafletjs. however my task is to do this without a map.

I'm tripping over the part of the code that reads the GeoJSOn and "foreach" poly

var cluster = [];
var point1 = turf.point([-73.79896352764388,42.81585623837577], { });//x,y

$.getJSON(url, function(data){
     ////  replace from here
    geojsonLayer = new L.GeoJSON(data,{style: style, onEachFeature: onEachFeature});
    map.addLayer(geojsonLayer);
    map.fitBounds(geojsonLayer.getBounds(), {padding: [50, 50]});

    geojsonLayer.eachLayer(function (layer) {
    ///// replace to to here....

    //$.each(data, function(layer){  // thought I was on the right track but...

    var isInside = turf.inside(point1,layer.feature);

        if(isInside){
            cluster.push(layer.feature);
        }

    });

    createTable(cluster);

});

Any ideas?


texts: Automatically refreshing ArcGIS93Rest layer in OpenLayers?

texts:
I finally succeeded to import a gtfs-realtime feed into ArcGIS Server (used gtfsrdb.py and the created a script to automatically, every 30 seconds, to update the geometry of a table with the points). It's working ok.

I imported my map service into an OpenLayers map. It's working too, except that I'd like my ArcGIS93Rest layer to redraw every 30 seconds as well.

For now, I have to change zoom level to see the new tile with the new position.

I tried OpenLayers.Strategy.Refresh, but either it's not working with Rest layers or I got it wrong...

var realtime = new OpenLayers.Layer.ArcGIS93Rest(
  "Temps réel",
  "http://xx.xx.xx.x:6080/arcgis/rest/services/trains_rt/MapServer/export",
     {transparent: "true", format: "image/png"},
     {isBaseLayer: false, singleTile: true, visibility: true, format: "image/png"},
     {strategies : [new OpenLayers.Strategy.Refresh({interval: 10000, force: true})]}
  );


texts: Leaflet: How to create toggleable overlays from a single GeoJSON FeatureCollection?

texts:
I am a beginner to Leaflet/Mapbox and I can create a simple map and plot all the points from a GeoJSON FeatureCollection, but I am unsure how to add the ability to toggle (via checkboxes) specific points based on a property named category.

How would I list all the categories layers in my control and give the user the ability to toggle on/off the various layers, as well as a main master layer that shows ALL of the points? Demo FIDDLE here

Here's my example map with 3 example categories in a GeoJSON FeatureCollection:

var collection = [{
"type": "FeatureCollection",
    "features": [
    {"type":"Feature","properties":{"category":"Aviation","Name":"Plant No 1"},"geometry":{"type":"Point","coordinates":[81.73828125,62.59334083012024]}},
    {"type":"Feature","properties":{"category":"Electrical","Name":"Plant No 2"},"geometry":{"type":"Point","coordinates":[94.5703125,58.722598828043374]}},
    {"type":"Feature","properties":{"category":"Military","Name":"Base 1"},"geometry":{"type":"Point","coordinates": [104.4140625,62.91523303947614]}}
  ]}];

The code:

var myStyle = { radius: 10, fillOpacity: 1, stroke: false, weight: 1, opacity: 1, fill: true, clickable: true };

var allPoints = L.geoJson(collection, {
    pointToLayer: function(feature, latlng){
        return L.circleMarker(latlng, myStyle);
    },
    style: function(feature){
        switch(feature.properties.category){
            case 'Aviation' : return { color: "black" };
            case 'Elecrtical' : return { color: "blue" };
            case 'Military' : return { color: "red" };
        }
    },
    onEachFeature: function(feature, layer){
        layer.bindPopup(feature.properties.Name);
    }
});

var basemapsObj = { 'basic': L.mapbox.tileLayer('mapbox.light'), 'satellite': L.mapbox.tileLayer('mapbox.satellite-afternoon') };

var overlaysObj = { 'All points': allPoints.addTo(map) }

L.control.layers(basemapsObj, overlaysObj, { collapsed: false }).addTo(map);


texts:
I am using the ArcGIS Javascript3.22 API routeTask to determine the route around polygon barriers.

How do I properly define the polygon barriers?

What I have tried is:

  routeTask = new RouteTask("https://route.arcgis.com/arcgis/rest/services/World/Route/NAServer/Route_World");
  routeParams = new RouteParameters();
  routeParams.polygonBarriers =  new FeatureSet("http://arcgis-geog865-942335747.us-east-1.elb.amazonaws.com/arcgis/rest/services/RoadLimit_feature_server/FeatureServer/0");

But, when I include the last line above (routeParams.polygonBarriers...) in my code the page stop's working.


texts: Loading temporary layer in processing PyQgis 3.4

texts:
I have been working on processing script in PyQgis 3.4 and i have serious issues with loading temporary layers to project after processing. Since following code works perfectly fine:

entimport processing
layer = iface.activeLayer()
features = layer.getFeatures()

for feature in features:
    print(feature['FIELD_LEN'])
    field_len=feature['FIELD_LEN']
    spacing=feature['SPACING']
    spacing_corr = (field_len/spacing) + 1
    processed_lyr = processing.run('qgis:pointsalonglines', {"INPUT": layer, "DISTANCE": spacing_corr, "START_OFFSET":0, "END_OFFSET":0, "OUTPUT":'d:/test.shp'})

I wanted for output to be temporary layer, so I switched "OUTPUT:"

processed_lyr = processing.run('qgis:pointsalonglines', {"INPUT": layer, "DISTANCE": spacing_corr, "START_OFFSET":0, "END_OFFSET":0, "OUTPUT":'memory:'})
QgsProject.instance().addMapLayers(processed_lyr["OUTPUT"])

but  I have no idea how to acess it since it produces error

TypeError: QgsProject.addMapLayers(): argument 1 has unexpected type 'QgsVectorLayer'

What should be used instead of .addMapLayers?


texts: ArcSDE: Specifying fields to not allow null values

texts:
I have a feature class in SDE that has various fields where the setting AllowNulls is set to FALSE.

When creating a feature in ArcMap in this feature class it still allows nulls to be specified in the attributes, albeit an empty text string.  I can then save my edits without any error messages.

What's going on?

Thanks
Rob


texts: QGIS Densify Geometries based on expression

texts:
I tried using the Expression Based Input in the Densify Geometries tool with the following expression:

attribute(  $currentfeature ,'VertNum')

where VertNum is a field in the attribute table that contains the number of vertices I want to add to each feature. But I just get a "NULL" response in the Output Preview, and when I click "OK" the number of vertices is still set to 1.

I've tried with different variable types (Int, Double) without success. Am I using the Expression Based Input tool for its intended purpose here?


texts: Tiled map image only shows small section of map but has huge file size

texts:
My situation is this: I've purchased some PDF maps that I converted some time ago into some software (MudMap) and I now no longer have the PDFs and I want to convert them back.

My problem is this: From the software, I've got the data file which is about 23MB and the xml file which describes tile information etc. I don't care about the metadata as I just want the map image so my question is with the data file only. As I said, it's 23MB. If I rename it to image extensions, it opens with image viewers but it is only 256x256px and only shows a section of the map - a section that looks like it would only 256px of the whole image. Opening the image in a text editor shows some repeating text elements with "IHDR" so I guess the file is related to a PNG format...??

So my question is: Simply, how can I view the full image. Does the PNG format (or others) have a capability to store additional image information hidden out of view? Or does the file itself set the boundaries of it's viewable area. Finally, is there software out there that reads these files?


texts: Creating Accuracy Assessment Points gives Error 999999 in ArcGIS Desktop?

texts:
I am using ArcGIS Desktop 10.6.1.

I am trying to create 500 stratified random accuracy assessment points for my raster generated using the supervised classification tool (followed by the reclassify tool to combine multiple agricultural classes into a single value).  I open the "Create Accuracy Assessment Points" tool, select my reclassified raster, specify a name and file geodatabase for the output, and leave the default settings for the other parameters (Classified target field, 500 points, Stratified_Random Sampling strategy).

I then get a failed analysis, and under Results: Messages, it states, "Error 999999: Error executing function." and "Failed to execute (CreateAccuracyAssessmentPoints).  The ESRI help page offers two possible reasons for this: the input file does not contain an integer field with the name "Value," or you try to use the Stratified Random sampling strategy with a class that has a count of less than 10 pixels.  I verified that I do have a field in my attribute table named "Value" with the "Long integer" class.  I also checked the counts- none of my five values have pixel counts less than ~55,000.

To try to identify what the problem might be, I tried running the analysis with a Random sampling strategy instead, but got the same error.  I also tried creating a new map document, resetting the only environment variables that I had adjusted (processing extent, snap raster, and raster mask), and processing a different classified raster.  I get the same error each time.

Does anyone have ideas on what might be causing the problem, or even other ways that I might troubleshoot this particular analysis?


texts: I have a flex viewer 3.0 map in wgs 84 (4326) (ESRI basemap) with a wms which is 4326 as well. for some reason the wms doesn't overlay correctly. it is way off. appears on the maps but it is flattened quite a bit.


texts: OpenLayers 3 - Modify Feature - Delete vertex on keypress

texts:
I'm trying to create functionality in OpenLayers 3 that allows me to delete a vertex from a feature that is currently selected using the modify interaction.

I've managed to see that there is a 'deleteCondition' parameter set on the modify interaction when it is created, and that this needs to return an ol.events.condition object. At present I can see from this example (http://openlayers.org/en/master/examples/draw-and-modify-features.html) that is it possible to delete a vertex by pressing SHIFT and mouse click at the same time, however this then creates a new point which seems counter-intuitive to me. If you are deleting a vertex why would you want to create a new vertex in exactly the same place? I may have missed something here though.

So the question is how can I use these ol.events.condition objects to make it so that if the user clicks a specific key on the keyboard the highlighted vertex is deleted?

My initial thoughts were to have something like this:

var modifyInteraction = new ol.interaction.Modify({
    features: selectInteraction.getFeatures(),
    deleteCondition: function(event) {

        var key = event.originalEvent.keycode || event.originalEvent.charCode;
        if(key == 46) {
            return  ol.events.condition.always(event);
        } else {
            return ol.events.condition.never(event);
        }
    }
});

But this produces some very strange results...


texts: I don't see the GRASS toolbar in QGIS 2.10 (Pisa). There are still GRASS commands in the Processing toolbox. I know there was some upgrade of GRASS for 2.10. Is the toolbar not available anymore? If so, how would you make a mapset, etc?


texts: Anyone have recommendations on data collection system for bird point counts?

texts:
We will be conducting point counts this spring and summer, and would like to move away from paper data entry forms (yeah, I know).  I have spent a little time exploring ArcPad (10), there are significant bugs (using "Find" to get to a point doesn't work, data transfer is unreliable, etc).

At each point, we will need to capture a GPS location (one record in a table of points), then will record observations for birds at the associated point (one record per bird in a table of observations).  The once nice thing about ArcPad is that it's easy to set up the one-to-many relationship between points and observations, and enter as many observations as needed.  But I'm just too nervous sending crews out to collect a few weeks' worth of data without feeling 100% confident we'll get the data back.  Also, the Trimble units are too expensive, and the screens are just too small for easy data entry.

I was wondering whether anyone had explored options using Android tablets, or something similar, for data collection in the field...?  In the past, I used a Trimble data collection unit with Pocket Excel.  This kind of system would work, but the one-to-many link depends upon the user entering a point ID correctly in both the point table and the observation table.  I like the idea of a cheap Android tablet with 3 or 4g connectivity running an app that allows data collection in this manner, and would be fine sending crews out with a separate handheld GPS, if accuracy is too poor with a tablet.

Would really appreciate any ideas or suggestions for data collection.


texts: Adjusting resolution of spatstat::pixellate.psp() output

texts:
I am trying to accomplish exactly what is being done in this question: Convert line shapefile to raster, value=total length of lines within cell  using spatstat::pixellate.psp() technique (see third answer) but I cannot figure out how to control and adjust the resolution of the output raster...

Spatstat documentation says to use the "Window" parameter to adjust resolution but I'm getting the same output resolution whether I use it or not.  Also, the third answer in the question cited above uses a dimxy= parameter that is not producing any change in my output resolution no matter what number I input.


texts:
I found some conversation about listing the source paths for multiple mxd files and used the script provided at Crawling directory and listing all datasources using ArcPy?.

It worked great on one folder but on the other two folders I get this error:

  Traceback (most recent call last):   File
  "\de1-s07\Q-Drive\Jobs\Missouri\I_70_PEL_StLouis\AECOM\Map_Docs_MxdCrawler\mxdcrawler.py",
  line 33, in
      main(folderPath, output)   File "\de1-s07\Q-Drive\Jobs\Missouri\I_70_PEL_StLouis\AECOM\Map_Docs_MxdCrawler\mxdcrawler.py",
  line 9, in main
      w.writerows(rows)   File "\de1-s07\Q-Drive\Jobs\Missouri\I_70_PEL_StLouis\AECOM\Map_Docs_MxdCrawler\mxdcrawler.py",
  line 18, in crawlmxds
      dataframes = arcpy.mapping.ListDataFrames(mxd)   File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\utils.py", line 182, in fn_
      return fn(*args, **kw)   File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\mapping.py", line 1496, in
  ListDataFrames
      result = mixins.MapDocumentMixin(map_document).listDataFrames(wildcard)   File
  "C:\Program Files
  (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\arcobjects\mixins.py", line 730,
  in listDataFrames

Each time I run the script it fails at a different point. One time it only read 2 files and this time it read 9 files. The python script is:

import arcpy, os, csv

def main(folder, outputfile):
    with open(outputfile, "wb") as f:
        w = csv.writer(f)
        header = ("Map Document", "MXD Path", "DataFrame Name", "DataFrame Description", "Layer name", "Layer Datasource")
        w.writerow(header)
        rows = crawlmxds(folder)
        w.writerows(rows)

def crawlmxds(folder):
    for root, dirs, files in os.walk(folder):
        for f in files:
            if f.lower().endswith(".mxd"):
                mxdName = os.path.splitext(f)[0]
                mxdPath = os.path.join(root, f)
                mxd = arcpy.mapping.MapDocument(mxdPath)
                dataframes = arcpy.mapping.ListDataFrames(mxd)
                for df in dataframes:
                    dfDesc = df.description if df.description != "" else "None"
                    layers = arcpy.mapping.ListLayers(mxd, "", df)
                    for lyr in layers:
                        lyrName = lyr.name
                        lyrDatasource = lyr.dataSource if lyr.supports("dataSource") else "N/A"
                        seq = (mxdName, mxdPath, df.name, dfDesc, lyrName, lyrDatasource);
                        yield seq
                print "running"
                del mxd

if __name__ == "__main__":
    folderPath = r"\\de1-s07\Q-Drive\Jobs\Missouri\I_70_PEL_StLouis\AECOM\Map_Docs\Restart" # or arcpy.GetParameterAsText(0)
    output = r"\\de1-s07\Q-Drive\Jobs\Missouri\I_70_PEL_StLouis\AECOM\Map_Docs\_mxdcrawler\mxdcrawler_Restart.csv" # or arcpy.GetParameterAsText(1)
    main(folderPath, output)

print "done"

What do I need to fix to get this to work?

I am using python 2.7.8 through ArcGIS 10.3 GUI.


texts: Integrate GRASS tool in custom python script, QGIS 2.8

texts:
I'm trying to integrate "r.null" tool in QGIS 2.8 script editor.
When I try to execute the code I get error message that 'outputs_GRASS7R' is not recognized.

This is the code I'm using:

outputs_GRASS7R.NULL_1=processing.runalg('grass7:r.null', outputs_GDALOGRRASTERCALCULATOR_1['OUTPUT'],None,SUN_ELEVATION,False,False,False,False,False,None,0.0,sun_elevation_raster)

How can I fix that? do I need to import some specific module? (I tried 'import grass, 'import grass.scripts'. I have the GRASS plug-in installed in QGIS.

And I need to be able to send this code to other people that do not understand python.. so hopefully there is a simple solution.


texts:
I'm looking for a tool or a plug-in for QGIS that would allow me to shift rasters to new coordinates. An equivalent of ArcGIS' Shift tool would be the best option.

I know I can edit the world file that comes along with the raster to do this, but some raster formats (like geotiff) contain their metadata inside their binary format.


texts:
I'm observing the following annoying issue after styling layers by categories: All the categories show up on the Layers panel with a very large spacing between categories.
I'm not finding any settings to change this. Anyone seen and fixed this before ?

Please see attached image, for sure will help to understand the issue.

Using
QGIS version 3.10.1-A Coruña


texts: "Sequence size" doesn't match row in Python?

texts:
Essentially I am trying to calculate a value for a couple of fields after adding them. Below you can see how I am adding the fields and then an example on how my "for" statement is looping through the cursors to add values.

I keep getting an error: "TypeError: sequence size must match size of the row"

I'm not sure what to do about this. Is there a better method than to use the UpdateCursor updateRow statements?

fieldName = ["azVEqSet", "azVERis"]
try:
    arcpy.AddField_management(points, fieldName[0] ,"FLOAT")
    arcpy.AddField_management(points, fieldName[1] ,"FLOAT")
except:
    arcpy.DeleteField_management(points, [fieldName[0], fieldName[1]])
    arcpy.AddField_management(points, fieldName[0] ,"FLOAT")
    arcpy.AddField_management(points, fieldName[1] ,"FLOAT")

cursor = arcpy.da.UpdateCursor(points, ["azVEqSet", "azVERis", "POINT_X", "RASTERVALU", "POINT_Y"])
for row in cursor:
...
     numberInsert = ReturnedVariable
     row[0] = numberInsert
     cursor.updateRow([row])


texts:
I'm running code that works with MODIS HDF4 datasets on a Debian Linux distribution.  I am also running this code on a Windows machine.  On both the Windows and the Linux machines, I have the most recent GDAL version (gdal 1.11.1, np19py27_1), however on the Linux machine, GDAL doesn't seem to recognize HDF4 Drivers... it won't load or process any of the HDF4 datasets that are easily processed on the Windows machine (unfortunately this machine doesn't have enough space to run the full code).

Has anyone else encountered this problem?  Do I need to install an additional module on the Linux system?


texts: Certain polygons not obeying JavaScript style rules?

texts:
Styling misbehavior http://catalinagis.org/files/geojson_opacity.jpg

http://catalinagis.org/ownership

You can see that some of the parcels don't follow the fill style rules (opacity: 0.5) even though the same style rule applies to all polygons in the data.  This happens when the data is loaded into http://geojson.io and http://mapbox.github.io/geojson-vt/debug/

Any ideas as to what is happening?


texts: Is there any labeling (by scale) rule of thumb in Cartography?

texts:
I did some research about Map labeling and most of them are about font, spacing, etc.

I am just wondering is there any rule of thumb of labeling in terms of map scale.

For instance, in Google Maps
map zoom level 3, Country Names and Ocean Names
map zoom level 4, Big Cities and Seas
... and so forth.

I tried to do some search and somehow I couldn't find any.

Can you give me some advice or places where I can start to do research?


texts:
I would like to know if in PostGIS is possible to define spatial constraints for avoid topology errors?
I have a polygon table and I need to avoid self-intersection and gaps.
Does anyone know more about that?


texts: Counting how many times a point is inside a set of intersecting polygons in R

texts:
I have a set of points (latitude and longitude) and a set of polygons with the same CRS. Since some of the polygons are interescting, it may be that some of the points are in more than one polygon.

How can I count how many times a point lies in all the set of polygons?

inside <- data.frame(as.numeric(!is.na(over(points, as(get(polygons), "SpatialPolygons")))))

This line of code works only to generate a dummy variable equal to 1 if the point is contained in at least a polygon and 0 otherwise.


texts: Open ArcGIS Business Analyst Vector Tile Service in QGIS 3.14

texts:
Don't even know if this is possible, but I am trying to open up a Vector Tile Service that is subscription content I have through Esri Business Analyst Online. I am able to open in ArcGIS Pro (of course), but I would like to be able to use in QGIS.

I attempted to add as:

ArcGIS Map Service (apeared as an empty raster layer)
ArcGIS Feature service

Adding as vector tiles does not give you the option to enter username and password as the other 2 options do.
Some more info. When I save in ArcGIS as a layer file and open it in a text editor, this is what appears:
{
  "type" : "CIMLayerDocument",
  "version" : "2.5.0",
  "build" : 22081,
  "layers" : [
    "CIMPATH=9062411281594e2a83fcd1acb8bc5817.xml"
  ],
  "layerDefinitions" : [
    {
      "type" : "CIMVectorTileLayer",
      "name" : "CoreLogic Parcel Point Boundaries",
      "uRI" : "CIMPATH=9062411281594e2a83fcd1acb8bc5817.xml",
      "sourceModifiedTime" : {
        "type" : "TimeInstant"
      },
      "useSourceMetadata" : true,
      "layerType" : "Operational",
      "maxScale" : 564.24858800000004,
      "minScale" : 295828763.7957775,
      "showLegends" : true,
      "visibility" : true,
      "displayCacheType" : "Permanent",
      "maxDisplayCacheAge" : 5,
      "showPopups" : true,
      "serviceLayerID" : -1,
      "refreshRate" : -1,
      "refreshRateUnit" : "esriTimeUnitsSeconds",
      "dataConnection" : {
        "type" : "CIMVectorTileDataConnection",
        "uRI" : "https://www.arcgis.com/sharing/rest/content/items/7f3041a9a89f432fbc9d536919dab1ac/resources/styles/root.json"
      }
    }
  ]
}


texts: Comparing two shapefiles based on ID and spatial

texts:
I have two shapefiles that contain position information for individuals, one is their reported location (line segments) and the other is their high-frequency recorded location (point data); they have an ID field that is common between them. What I am looking to do is identify the reported locations in the larger recorded data set with a spatial selection. The problem is my data is so condensed that I cannot do a simple select by location. There is too much overlap of individual locations. I also have over 1,000 individual IDs so I don't want to individually query each ID.

What I am looking for is a way to use possibly a cursor to go through each unique ID, query only the records in each table that share that ID, then do the spatial selection. I am very new to cursors and scripting in general so I don't know if a cursor can be used for two tables. (This may be irrelevant but for this process I can convert my line data to point data if having the same type of shapefiles is important.)

(Removed previous code that I am no longer using)

I am also unaware if maybe Make Query Table will work for this (I have never used it).
 Thanks!

UPDATE: Here is a new code that I worked on with a fellow GIS person. GOOD NEWS! After leaving ArcGIS alone for a while, my code ran! It went as expected, but only did one iteration on one ID. The iteration should be implied in the Cursor as far as I know, or is there a line of code I'm missing that tells the Cursor to repeat for all IDs?

rows = arcpy.SearchCursor("hmb_logbook_2009")
for row in rows:
    whereClause = '"trip_id" = ' + str(row.trip_id)
    arcpy.MakeFeatureLayer_management("hmb_logbook_2009", "currentLines", whereClause)
    arcpy.SelectLayerByLocation_management("hmb2009", "WITHIN_A_DISTANCE", "currentLines", "1 Miles", "NEW_SELECTION")
    arcpy.SelectLayerByAttribute_management("hmb2009", "SUBSET_SELECTION", ' "trip_id" = ' + str(row.trip_id))
    arcpy.CalculateField_management("hmb2009", "match_activity", 1)
    arcpy.SelectLayerByAttribute_management("hmb2009", "CLEAR_SELECTION")
    arcpy.Delete_management("currentLines")
del row, rows

UPDATE 2: So I have troubleshooted, and the reason the code only ran once is this error associated with my Make Feature Layer tool. Here is the error message: ERROR 000622: Failed to execute (Make Feature Layer). Parameters are not valid. ERROR 000628: Cannot set input into parameter in_features.
I tried putting the full address into that tool (only that one) instead of referencing the layer in the map and it works!


texts:
I want to increase icons and texts

icons and texts of Toolbar, Table of Contents and Search window is very small!!!

But, Texts in menu-bar is proper size. I don't understand this situation

Now, I use window 10 and arcmap 10.1


texts: Connecting Access Table to ArcGIS (Replace Function!)

texts: I have an Access Query where the Field names contain symbols like "-", ":" and spaces. When I was connecting my database to ArcCatalogue that specific Query-table was not visible. Then I used Replace function (in Access) to replace those symbols and spaces to "_". Then the table was visible in ArcMap. I did lots of mapping with that and suddenly after few weeks that Table is no longer connecting. My all ArcMap data links are broken. I read somewhere that the Replace function does not work out of Access environment. Then how was it working 2 months ago and now its not connecting? Is anyone experiencing similar problem? any suggestion?


texts: QGIS crashes (dumped crash) while importing SHP layers or SpatialLite

texts:
I meet problems with the last version of QGIS (Lyon 2.12.3) 64bits on Windows 7.

The software crashes as soon as I load a SHP layer or use SpatiaLite layers (created from OSM data).

Besides, each time I open QGIS desktop, I get the following error message:

I tried to uninstall/install again the software, to empty all directories with QGIS stuffs, to empty temp directories, nothing changes.

Moreover, even after unstalling QGIS, when I try to install it again, I get the message that QGIS is still installed on my computer.

I cannot use the app anymore...

Does anyone could help me ?


texts: Get all tileLayers from L.map object in Leaflet

texts:
What is the best way to access all displayed tileLayers from a L.map object?

Traditionally, I have stashed tileLayers in an object that I can access later, like so:

var mapLayers = {}

if (! mapLayers[mapId]){
  mapLayers[mapId] = L.tileLayer(url);
}
map.addLayer(mapLayers[mapId]);

This works well enough, but I keep running into issues with it, and it seems like bad practice to store this data in an object I've created when it's already stored in the map object.

Internally, the map object stores tileLayers in map._layers, along with other layers (grids, feature layers...).  They're indexed according to some internal leaflet scheme, so not sure if this is the best point of entry.

Using leaflet 0.7.3


texts: How can I send a parameter with a RouteTask and get it back in RouteResult

texts:
I'd like to send a parameter with the Route Task and be able to get that same parameter in the RouteResult object.
The goal is to be able to recognize the routes I receive from the server when I perform many solve operations on the same RouteTask.


texts:
I've managed to create a Point shapefile creating a SimpleFeature and outputting it as Test_Point.shp,

public static final SimpleFeatureType POINT = createPointType();

private static SimpleFeature createPoint(){
    Coordinate start = CoordConverterUtils.getStart;
    Coordinate end = CoordConverterUtils.getEnd;

    Coordinate mid = (new LineSegment(start,end)).midpoint());

    GeometryFactory gf = JTSFactoryFinder.getGeometryFactory();
    SimpleFeatureBuilder fb = new SimpleFeatureBuilder(POINT);
    Point p = gf.createPoint(mid);
    fb.add(p);
    return fb.buildFeature(null);
}

private static SimpleFeatureType createPointType(){
    try{
        SimpleFeatureTypeBuilder builder = new SimpleFeatureTypeBuilder();
        builder.setName("Point");
        builder.setCRS(DefaultGeographicCRS.WGS84);
        builder.add("the_geom", Point.class);
        return builder.buildFeatureType();
    catch(Exception e){
        //warning
        return null;
    }
}

In the same directory I have a sld file Test_Point.sld that is from https://docs.geoserver.org/stable/en/user/styling/sld/cookbook/points.html

<?xml version="1.0" encoding="ISO-8859-1"?>
<StyledLayerDescriptor version="1.0.0"
    xsi:schemaLocation="http://www.opengis.net/sld StyledLayerDescriptor.xsd"
    xmlns="http://www.opengis.net/sld"
    xmlns:ogc="http://www.opengis.net/ogc"
    xmlns:xlink="http://www.w3.org/1999/xlink"
    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">
  <NamedLayer>
    <Name>Point</Name>
    <UserStyle>
      <Title>SLD Cook Book: Simple Point</Title>
      <FeatureTypeStyle>
        <Rule>
          <PointSymbolizer>
            <Graphic>
              <Mark>
                <WellKnownName>circle</WellKnownName>
                <Fill>
                  <CssParameter name="fill">#FF0000</CssParameter>
                </Fill>
              </Mark>
              <Size>6</Size>
            </Graphic>
          </PointSymbolizer>
        </Rule>
      </FeatureTypeStyle>
    </UserStyle>
  </NamedLayer>
</StyledLayerDescriptor>

Using Geotools - JMapFrame I can open/ view my Point shapefile, but can't seem to get the styling to be applied. Not sure if this is because

The .sld is not done properly
Geotools JMapFrame and filechooser can only open/ display .shp files


texts: Extract terrain coordinates as seen from a viewpoint + FoV in a digital globe

texts:
Context
Every one working with photogrammetry or 3D knows what a depth image is. For example OpenCV defines a "depth map" from a stereo pair, e.g.:

where the darker a pixel is, the greater the distance to the camera center.
Needs
I need to go a step further.
Let's say I have a triangulated DEM in a digital globe (in my particular case it's a terrain mesh rendered in CesiumJS but I also have the corresponding DEM as a GeoTIFF).
I need to figure out a mean to extract the terrain WGS84 coordinates as it is seen from a simulated camera at a given point of view. This point if view is defined as follows with known values for:

the position,
the 3D orientation,
the focal length,
the image size (width, height and corresponding numbers of pixels, (x,y))

All these values define the 3D cone of visibility, a.k.a. field of view.
As an example, the frame of such an image is depicted here in green (this rectangular shape can be given as a GeoJSON, KML or whatever format in terrain space):

This Cesium screenshot example simply corresponds to a Viewer with a base imageryProvider and a terrainProvider:
var viewer = new Cesium.Viewer('cesiumContainer', {
    imageryProvider : new Cesium.BingMapsImageryProvider({
        url : 'https://dev.virtualearth.net',
        key : `${BING_MAP_TOKEN}`,
        mapStyle : Cesium.BingMapsStyle.AERIAL
    }),

    terrainProvider: new Cesium.CesiumTerrainProvider({
        url : `https://api.maptiler.com/tiles/terrain-quantized-mesh/?key=${MAPTILER_APP_TOKEN}`
    }),
    baseLayerPicker : false,
    infoBox : false,
});

RGB values from terrain should be mapped as R=longitude, G=latitude, B=ellipsoidal height.
All these 3 values corresponds to the terrain coordinates for each defined pixel of the image.
Any 2D tool with some kind of viewshed analysis where one can define a picture frame and its internal parameters are also welcome, but here, I preferred to represent the story using real 3D.
Any hints are welcome.
I'm working on a full open source stack (Ubuntu + GIS tools).


texts: Blank Output From Load raster in ArcCatalog

texts:
I did Load raster data in ArcCatalog, I wait until it finished, but when I opened the Output raster, I just got a White-Blank output, there isn't anything in there.
What happend? And how I fix it?
NB: I did it (with same parameter) to another image and it succeeded, just this one that failure.


texts: Are there open access terrain tiles with a permissive license like OSM?

texts:
Are there open access terrain pyramidal map tiles with a permissive usage license like OpenStreetMap's? I'm building aviation-related maps with Polymaps and need a good raster basemap. The Google terrain tiles are just perfect for my purposes: restrained colours, not too many features, and a clear view of terrain. Unfortunately Google does not allow use of their map tiles in any library but their own Javascript.

Alternatives I've tried.. OpenStreetMap's maps all seem to be road maps. The OpenCycleMap does include terrain, but then it has a lot of other high contrast features on top that make it unsuitable for my needs. USGS has some terrain tiles at the National Map, but they're unstyled and lack any descriptive features. I know it's possible to generate and serve my own tileset using a free datasource like SRTM, but that's an awful lot of work to design and operate.

I'm able to pay a modest cost for bandwidth; mostly I want to avoid spending the time generating my own tileset. Terrain maps seem like a common need, I'm surprised I haven't found one!

Edit since posting this question I've found a few more options. Microsoft's Bing Maps include terrain tiles; with somewhat flexible, but limited terms of service. GeoiQ/GeoCommons has the Acetate tiles which are specifically designed for my purpose.

Edit 2 this post had a small part in inspiring Mike Migurski / Stamen's terrain map based on open data. More at maps.stamen.com.

Edit 3 MapBox also now has a terrain layer. I'm unclear about the licensing, but the company is generally open source friendly.

Edit 4 MapZen now has DEM and raster terrain tiles. Again unclear about the licensing, but the company is generally open source friendly.


texts:
I'm trying to create a buffer

I'm using QGIS 2.4.0
layers and have street departments and rivers with coordinate system EPSG: 4326 - WGS 84

I need to create a buffer layer 500 meter rios

This is my rivers layer

but my result look like this


texts: ArcGIS server 10.1 REST API queries fail inconsistently

texts:
I have a set of map and ESRI feature services running on a 10.1 ArcGIS server with SQL Server back end. When I query the layers for either the geometry or oids using the REST interface, the query only works in certain instances.

If I query a point layer the result works as expected and I get the data
If I query a polygon layer like this, it fails: http://gistmapsdev.itos.uga.edu/arcgis/rest/services/COD_External/SYR_AR_DynamicLabeling/FeatureServer/5/query?returnGeometry=true&f=json
I get the error message:

{"error":{"code":400,"message":"Unable to complete operation.","details":["Unable to perform query operation."]}}

The server log files show:

Database error: Operation without query criteria is not allowed.

However this query (same service, different polygon layer) works as expected and returns the feature geometry:
http://gistmapsdev.itos.uga.edu/arcgis/rest/services/COD_External/SYR_AR_DynamicLabeling/FeatureServer/6/query?returnGeometry=true&f=json

There are 2 polygon layers that I can query and get the expected result but all other polygon layers fail and for the life of me I cannot see any difference in the properties of the feature classes that fail vs those that do not.

If I dirty the query (by adding a where clause 1=1) I get the expected result; so for example, the query in point 2 above works when I add where=1%3D1

What I’ve looked at:

a.  I have not hit the max feature limit set on the server or the service

b.  the behavior is consistent across browsers

c.  the behavior is consistent using either get or post

d.  I have rebuilt the spatial indices and cleared the caches

e.  The behavior is consistent if I ask for json or html response

This is different than the NIM086349 bug described on this site at 52612/arcgis-server-10-1-inconsistent-querying-errors - which fails on alternate queries. In my case the failure is consistent and appears to be tied to particular data types and data layers.

I’d be grateful for any pointers or areas that I should look at.


texts: Is there a way to query Esri's map service rest endpoint but only get the bounding box of each feature (or a single feature)?

texts:
I want to be able to zoom a user to a state of their choosing. But i don't want to load all of the vertices of each state. instead I want to query the rest endpoint and get the 4 bounding box coordinates for each state. Ideally, i would store those client side. Failing that, i would like to be able to query for the bb of a single state so i can zoom a user to that.

I am hitting against an ESRI ags 10.3.1 server


texts: Masking layers by a vector polygon layer in QGIS

texts:
I have a vector layer of buildings:

and a polygon layer coloured by some quantity:

and I want to mask one by the other to produce this:

This can be done by setting the blend mode to "Multiply" with the building polygons on top as long as the polygons are white and the background of that layer is black.

One way of doing that is to use an "Inverted Polygons" fill but that is very slow at low zoom levels since all the buildings are in view and that's a lot of inverted polygons.

The other way I thought I'd cracked this was to create a "Group" layer with the buildings filled in white and a virtual layer consisting of the bounding box of the buldings filled in black. On its own it looks right:

but QGIS (3.6) doesn't let you blend grouped layers (Gimp and Photoshop probably do) so it doesn't work.

I really want to do it in the canvas so solutions involving processing intersections and creating new layers and so on aren't going to work either.

I'm not sure if I'm missing a trick to create a style for a layer to have a black background, or if there's something else in the blending modes that will do this...


texts: Why does mod_tile have a limit of 10 XML files?

texts:
I recently discovered that mod_tile does not, by default, allow more than 10 XML configs to be specified in renderd.conf. Upon closer inspection, I found that I could change that limit in here by setting XMLCONFIGS_MAX to a different value, and then recompiling. I also found that a value larger than ~300 will segfault renderd when started.

Is there a reason this limit is so low? Is it related to renderd crashing with a high XML config count?


texts: How to extract/crop an OSM file to include entries that intersect with a GTiff?

texts:
I am trying to get just the datapoints that intersect with the boundaries of a GeoTiff image for some internal testing on an algorithm.

Are there any tools that can perform this operation?

If not, how would I script this behavior?

I would also like to preserve as much as data as possible from the original OSM file such as tags and other attributes.

I seek a solution that uses GDAL.

Here are the steps I need to perform:

Get the bounding box of the GeoTiff
Find all entries that intersect with those boundaries
Save those entries to a new OSM file.


texts: Hide layers in a workspace from getCapabilities

texts:
I have a Geoserver (2.10.1) with one workspace, some layers and a layer group.

Most of these layers may not make sense alone so I want to allow the users to see only the group layer and a few layers I choose (WMS/WFS services - OpenLayers, getCapabilities, getMap, ...).

I try to read this and try to manage the Users, Groups, and Roles in Geoserver's admin page but I'm very confusing.

According to the user manual (link above):

  Providing access to layers is linked to roles. Layers and roles are
  linked in a file called layers.properties, which is located in the
  security directory in your GeoServer data directory

But the manual says nothing about how to manage the Users, Groups, and Roles in geoserver itself. I don't like to configure anything by editing files when I already have an interface to do it.


texts: Aerial imagery showing in different colour while pan/zoom GeoServer

texts:
I have image/tiff files of a small aerial imagery. I published it under ImageMosaic plugin (GDAL plugin installed), but when I pan or zoom the map in OpenLayers preview I am getting different colour for same position. Screenshots are attached.

What might be the reason for this? Did I need to tune anything?


texts: Oracle MapViewer / Map Builder - Tile algorithm

texts:
Oracle MapBuilder generates three types of layers - theme, base and tiles.
I try to create algorithm for tiled layers. Have you any idea how to decode this?

I mean algorithm for width and height of tiles.

For a tiles is generated XML like this:

<map_tile_layer name="TL_TEST" image_format="PNG" http_header_expires="168.0" concurrent_fetching_threads="3">
   <internal_map_source data_source="10.0.0.57_lrs" base_map="BM_AGC" bgcolor="none"/>
   <tile_storage root_path="/temp"/>
   <coordinate_system srid="2180" minX="-413374.575" minY="193748.248" maxX="1418971.699" maxY="766162.025"/>
   <tile_image width="256" height="256"/>
   <zoom_levels levels="11" min_scale="1.0" max_scale="1.0E7" min_tile_width="0.06773333184979038" min_tile_height="610782.0913333333">
     <zoom_level level="0" name="level0" description="" scale="1.0E7" tile_width="610782.0913333333" tile_height="610782.0913333333">
     </zoom_level>
     <zoom_level level="1" name="level1" description="" scale="1000000.0" tile_width="65440.938357142855" tile_height="65440.938357142855">
     </zoom_level>
     <zoom_level level="2" name="level2" description="" scale="100000.0" tile_width="6761.425365313653" tile_height="6761.425365313653">
     </zoom_level>
     <zoom_level level="3" name="level3" description="" scale="10000.0" tile_width="677.1420081300813" tile_height="677.1420081300813">
     </zoom_level>
     <zoom_level level="4" name="level4" description="" scale="1000.0" tile_width="67.73172195320296" tile_height="67.73172195320296">
     </zoom_level>
     <zoom_level level="5" name="level5" description="" scale="100.0" tile_width="6.773322418713311" tile_height="6.773322418713311">
     </zoom_level>
     <zoom_level level="6" name="level6" description="" scale="10.0" tile_width="0.677333243384311" tile_height="0.677333243384311">
     </zoom_level>
     <zoom_level level="7" name="level7" description="" scale="5.0" tile_width="0.3386666216921555" tile_height="0.3386666216921555">
     </zoom_level>
     <zoom_level level="8" name="level8" description="" scale="4.0" tile_width="0.2709332973537244" tile_height="0.2709332973537244">
     </zoom_level>
     <zoom_level level="9" name="level9" description="" scale="2.0" tile_width="0.13546665869200755" tile_height="0.13546665869200755">
     </zoom_level>
     <zoom_level level="10" name="level10" description="" scale="1.0" tile_width="0.06773333184979038" tile_height="0.06773333184979038">
     </zoom_level>
   </zoom_levels>
</map_tile_layer>


texts: Subset multiband tiff file using Python in ArcGIS for Desktop?

texts:
I have a multiband (1140) raster file in tiff format that I want to subset. Each band represents a prediction of monthly rainfall (12 bands for each year) over a 95 year period (12*95=1140). I need to subset the raster dataset into 10 year increments (120 bands), but cannot figure out how to do this with a multiband file. The output file would simply be a subset of the original raster dataset.

Do I first need to separate into individual bands? I have found a way to create a composite multiband raster file from individual bands (http://resources.arcgis.com/en/help/main/10.1/index.html#//00170000009p000000), but have been unable to find a method for subsetting a multiband raster file in ArcGIS 10.1 for Desktop.


texts:
I have an NDVI raster that I would like to reclassify into some arbitrary bins i.e. all values between 0 - 0.2, 0.2 - 0.4 etc. for spatial analysis. Could anyone advise how to:

1) reclassify a raster in QGIS?

2) then calculate spatial statistics e.g. in the context of NDVI i would like to be able to calculate ground cover by "live" vegetation and then ground cover relating to the "bins" described above and produce some maps and histograms?


texts: Where to get bounding boxes of places / POIs?

texts:
How and where to get bounding boxes of places / POIs? Preferably, for free.
For example, on OSM I can see a place drawing displayed on a may and with zooming-out one can see corners as latitude and longitude.

They are needed to test a coordinate against the bounding box in a query.


texts: Splitting GeoJSON column, what is in one dictionary?

texts:
I'm new to GIS. I can't find anything to solve my problem.
What instrument can I use to do this? (pandas, geopandas)


texts: Finding closest street name from point using MMQGIS?

texts:
I got a list of points spread over a city and their coordinates. I'd like to find out the name of the closest address (or at least street) from each point. I downloaded the plugin MMQGIS which allow to geocode address from coordinate but not the inverse.

Does anyone know a solution to do that ?


texts: How to use, at one-time, an expression to replace encoding errors in QGIS 2.18.14?

texts:
I need to replace several encoding errors in just one task.

I have this:

If I change the encoding to UTF8, the "Label" field is correct, but the error passes to the other field, "name".

Instead of using this expressão several times, how to assemble the expression to do it all at once?

replace( "name" , 'Ã³','ó')
replace( "name" , 'Ã©', 'é')
replace( "name" , 'Ã¢', 'â')
replace( "name" , 'Ãº', 'ú')
replace( "name" , 'Ã¡', 'á')
replace( "name" , 'Ã£', 'ã')
replace( "name" , 'Ã§', 'ç')
replace( "name" , 'Ãµ', 'õ')
replace( "name" , 'Ã?', 'Á')
replace( "name" , 'ÃÂ', 'í')
replace( "name" , 'Ãª', 'ê')
replace( "name" , 'Âº', 'º')


texts: OpenLayers.Format.GML.parseFeature() with GML v3 data

texts:
I have a partial GML object with the geometry type polygon. The GML version must be 3+ since the coordinates are wrapped in the gml:posList tag:

<gml:Polygon srsName="EPSG:25832" xmlns:gml="http://www.opengis.net/gml">
     <gml:exterior>
         <gml:LinearRing srsName="EPSG:25832">
             <gml:posList>
                 N points where N is an even integer.
             </gml:posList>
         </gml:LinearRing>
     </gml:exterior>
</gml:Polygon>

I am using OpenLayers 2.10 and I am trying to convert the GML object into an OpenLayers feature using the OpenLayers.Format.GML class.

When executing the parseFeature method, which later calls the parseGeometry.linestring I encounter this comment:

/**
* Two coordinate variations to consider:
* 1) <gml:posList dimension="d">x0 y0 z0 x1 y1 z1</gml:posList>
* 2) <gml:coordinates>x0, y0, z0 x1, y1, z1</gml:coordinates>
*/

And this line of code which assumes the dimension attribute to be present:

var dim = parseInt(nodeList[0].getAttribute("dimension"));

In OpenLayers.Format.GML.v3 the same line of code looks like this:

var dim = parseInt(node.getAttribute("dimension")) || 2;

The problematic difference for me is the || 2 part which assumes that a missing dimension means 2 dimensional - which is what I want.
Then why not use the v3 edition, but the OpenLayers.Format.GML.v3 does not implement the parseFeature method.

I am recieving GML features from various services and I do not control them.

What is the best approach to handle this problem?


texts: Is there a plugin with CAD tools for editing vector entities such as lines and polygons in QGIS 3.0?


texts: How to export WDPA dataset to text and shapefiles?

texts:
Please have some patience with me as I am not a GIS professional at all, I am merely trying to work with some of the file formats. Let me first explain what I am trying to accomplish:

The WDPA (World Database of Protected Areas) is, as the name says, an attempt to map all protected areas in the world, with an emphasis on biodiversity. Their work can be seen here:

http://www.protectedplanet.net/

You can browse all the areas right there on the site, search for them, etc. In addition, they have the following download options of the entire set (right below the search box):

CSV
KML
SHP
ESRI Web Service

This interests me greatly, since I want to integrate these protected areas into my own website, in the following ways:

Using Google maps, per country (so not all at once) show the polygons as overlay on the current Google Maps canvas (I have to stick with Google Maps here, due to other functionality)
Clicking on an area should show its metadata
Stretch goal: As my website is about photography, it would be incredible if a geotagged photo could automatically be associated with the protected area in which it was taken

To accomplish the above, I figure I need the following:

I need all textual metadata of all areas on the map in a way so that I can convert it into a format that is suitable for MySQL, the internal database I use. For example, each protected area has a unique ID and a country code, I need that data.
I need the shape data per protected area (or per country?). This is where I lack knowledge. Either way it needs to be in such a format so that it can be shown on Google Maps, and I need to know to which protected area (or country) it belongs

Some things I tried:

The first 3 download options provided by WDPA all seem the same, not sure if that is a bug on the site or not. I had no idea how to open them but eventually figured out I need ArcGIS and installed a trial version. I still don't know how to export anything in the format I desire though.

The 4th option seems most promising:
http://ec2-54-204-216-109.compute-1.amazonaws.com:6080/arcgis/rest/services/wdpa/wdpa/MapServer

There the entire set is offered via different web services. I tried opening the layer in Google Earth, but it keeps crashing on me, and I read you can't export layers from Google Earth? The other options are not familiar to me but after trying them I see no easy way to get the data I want in the way I want it.

I'm looking for your guidance in this question. This data set is a treasure and I am very excited to try to bring it into my own site, but I'm failing miserably at the first step.


texts:
I am attempting to rasterize a point layer using the Feature to Raster conversion tool. Each attempt, I receive the following error codes:

ERROR 010429: Error in GRID IO: CellLyrCreateInternal: Grid: n:... alrady exists
ERROR 010067: Error in executing grid expression

This layer does not already exist. Can attribute field type (double vs float vs integer) produce this problem? Each time I try to change field types, or create a new attribute column of a different field type in either ArcMap or ArcCatalog, the column will revert to double, no matter which type I set it as. Any thoughts?


texts:
I have to compute the relative pseudo-range for my student job, unfortunately, nobody can help me because they are not read into the topic enough. So I'm asking here:
My receiver has to have a 1500 bit buffer, so all 5 subframes are found and decoded. This works just fine. I have a buffer counter (num_bits_buffered), which counts after all subframes are decoded until 1500 new bits are found in the buffer.
Another module provides the code phase of the signal.
The Receiver has a function get_TOT(decoder, code_phase) providing the corrected time of transmission.
And here, my first problem occurred: the toc and toe of one satellite were one bit lower - causing the satellite to have a completely different time (16 sec difference). Therefore, it was useless for pseudo-range computing. Do I have to correct the clock time for the pseudo-range computing too?
Now to the main part:
I tried the following formula for pseudo-range computing (All variables unitless):
T_s = TOW * 6 + num_bits_buffered_s / DataFrequency(GPSL1) + code_phase_s / CodeFrequency(GPSL1)
P_s = (T_s - T_1) * c

I get values close to the real ones(+- 2 km), compared to the measure of the distance between my provided position and the computed satellite positions (the one SV with other toc excluded). If I use my PVT algorithm on the pseudo-range data, I'm 1.5km off in a different directions.
Is anybody here seeing my mistake?


texts: Cutting up a shapefile into smaller parts

texts:
I've downloaded the 450MB coastline shapefile from OpenStreetMap, I want to split it up in smaller parts. Is there a way to generate a shapefile with only the polygons contained in a smaller area?

The file comes with the water areas already split up in rectangular shapes.


texts: How to draw lines, dashed lines, routes and add images to the QGIS mapcanvas?

texts:
I need to draw a rectangle in a map canvas and inside it draw routed dashed lines and add small images to it. These actions are to be performed using the QGIS API.

Is there any place with examples that I can start with? All that I need is a direction of implementing this.


texts: How to change the data storage place in Geoserver?

texts: When add a new data storage in geoserver,the url will be a local address.What if I want to direct this url to a remote server address?The data are storing at the amazon server.


texts: How to create pop-ups with graphs on point layers?

texts: I have a PostgreSQL database which is refreshed every 10 minutes and I need to create graphs in a pop-up with the data. When I click, a pop-up (or annotation) should appear. How can I make it possible?


texts:
I have created a geopackage with a spatial view. The first column is a unique rowid that I have aliased as OGC_FID and the second column is a geometry column. When I load this layer into ArcMap 10.4 the query layer user interface appears but no features display after I click OK. I have registered the view in the gpkg_contents table as per the instructions. The same view works for a colleague in ArcGIS 10.6.1.
I even tried simplifying my view by removing any joins. It's a straight select from a single table with geometry. I can add the same table as a layer to ArcMap no problem.
Are spatial views broken or not supported in ArcGIS 10.4?
If not, what was the first version of ArcGIS that supported spatial views?