texts: British National Grid Sheet Reference Letters in Print Layout

texts:
I have managed to get a custom grid on the print layout to convert northing and eastings (XY) to ESPG:27700 grid references i.e NS00 through to NS99.

The problem I am coming across is when when NS99 crosses into NT00 (see image)

As you can see the x axis will change, but the Y axis on the right hand side should change to NT rather than NS, I know why it isn't - the expression based annotation is reading the NS query first and thus overring the NT query. The same happens on the X axis when you move down from NS to NX (it will read the NS query first).

Is there any way of getting the grid sidings to read independently from each other? Because this has been a real head scratcher for me...

Right now I have 55 different print layouts for each 100km square across the UK, but would love to just have one, and not have to fudge in labels when I have to map across a boundary!

Edit: expressions

case
when @grid_axis = 'x' and
(@grid_number >=200000) and (@grid_number <300000)
then'NS' + left(right(@grid_number,5),2)
when @grid_axis = 'y' and (@grid_number >=600000) and (@grid_number <700000)
then 'NS' + left(right(@grid_number,5),2)
when @grid_axis = 'x' and (@grid_number >=300000) and (@grid_number <400000)   then 'NT' + left(right(@grid_number,5),2)
when @grid_axis = 'y' and (@grid_number >=600000) and (@grid_number <700000)
then'NT' + left(right(@grid_number,5),2)
end


texts: QGIS: modifying each of the attributes in a selection of objects by selecting the attribute value of one of them

texts: In the "simultaneously modify attributes of all selected objects" window, the attributes which have different values for the objects are marked in yellow, but it does not show the list of the different values and it is not possible to pick a value from them to apply it to the whole selection. Is there a way to do that?


texts: Cell size (10m) in Raster properties differs from measured cell size (20m) in multiband (Sentinel) raster?

texts:
I've always thought that if raster resolution is 10m, the cell size in the raster is 10m, and therefore the length of pixel is 10m.

Now I have started to work with Sentinel 2 optical imagery. Sentinel 2 data contains multiple resolutions within a multiband raster: 10m, 20m, 30m and 60m. I have created a raster composite from all 13 bands. When I load back individual single band rasters from raster composite, here I see that all single band rasters have cell size 10 m, even when I can see differences in cell size by naked eye; and while measuring the pixel (cell) length, I see that the cell size is actually 20m.

I don't understand why in Raster properties of one band the cell size is 10 m, but true pixel length is 20 m? If the data have been re-sampled to 10 m (which does not make sense for pixel information, but is computationally possible), the pixel size should be 10 m as well? Or am I misunderstand the concept about what the cell size represent in ArcGIS?

I am working in ArcGIS 10.4, GCS: WGS84


texts: Install Geospatial Modelling Environment

texts: I install Geospatial Modelling (GME) Environment in ArcGIS 10.2. But I can't find where is. I search to ArcGIS and ArcToolbox and as a standalone App but I can't find it.


texts: How to stop openlayers from loading on every zoom?

texts:
I was thinking this sort of stuff was happening in OL by default, but apparently not?

I'm working on an application with a LOT of layers and every time we zoom or pan it starts reloading, which is very annoying as it just adds more and more time to load after the next zoom/pan.

Is there like a way to set some sort of timeout on the zoom/pan events for reloading or am I stuck with a slow application now forever and ever?


texts: Converting an arc-based resolution grid dataset to a fixed aspect-ratio "pixel like" resolution raster image dataset

texts:
I'm looking into the EMAG2 dataset (http://www.geomag.org/models/emag2.html) which claims to be a 2 arc minute resolution. I want to create square pixel tiles where pixels will be of an exact distance like 100m or 100Km. But how can this be done if the grid cells may be of different areas on difference places in the world?

To reiterate, if a dataset is 2 arc minute resolution, does that mean the grid cells change aspect ratio and area size at different coordinates in the world?

I'd prefer to use the GDAL/OGR toolset.

Edited according to clarification: How do you convert an arc-based (e.g. 2 arc-minute) resolution grid dataset to a square aspect-ratio "pixel like" (e.g. 100Km by 100Km pixels) resolution raster image dataset assuming that arc-based resolution grids have grid cells that have different area sizes at different coordinates? (this assumption may not be true, please tell me if so.)


texts: How to divide LINESTRING spatial network into equal length segments using R/sf?

texts:
I'm working on a project to model car crashes within the municipality of Milan using different statistical methods. I have a big dataset of all geocoded car accidents that happenend within the municipality in 2015 and I downloaded data of the street network from OpenStreetMap using the R package osmdata and the following query.

library(sf)
library(osmdata)
# Warning: It downloads a lot of data
geo_milan_highways <- opq("Milan, Italy", timeout = 250) %>%
  add_osm_feature(key = "highway") %>%
  osmdata_sf() %>%
  # Extract LINESTRING geometries
  magrittr::use_series("osm_lines")

The result is a sf object and every street in the municipality is mapped as one or more street segments using a LINESTRING geometry. I excluded the street segments that does not lie within the administrative boundary of the municipality and I filtered several highways since they are not relevant for this analysis (like pedestrian highway or footway).

geo_milan <- getbb("Milan, Italy", format_out = "sf_polygon") %>%
  st_transform(crs = 32632) #proj: UTM WGS84 32N

library(dplyr)
geo_milan_highways <- geo_milan_highways %>%
  st_transform(crs = 32632) %>%
  # filter highways inside the municipality of Milan
  slice(st_contains(geo_milan, .) %>% unlist()) %>%
  # Exclude some highways
  filter(
    !is.na(highway),
    !highway %in% c("motorway", "motorway_link", "trunk", "trunk_link", "pedestrian", "footway", "living_street", "service", "cycleway", "track")
  )

I clustered the geometries that were touching each other using the same approach described here to create a (sort of) highway network.

The result is a sf object where each LINESTRING geometry touches at least another LINESTRING geometry. The problems is that the street segments have very different lengths, from less than 1 meter to more than 3 kilometers. Is it possible to "rebuild" the network in such a way that every segment has (approximately) the same length?

I know that there are some functions to create a spatial network from a sf (like SpatialLinesNetwork of dodgr package) but I don't know how to use it to "refactor" the network with segments of the same length. There is also a function named st_linesubstring of the lwgeom package with an example on how to divide a LINESTRING into LINESTRINGS but it does not consider the neighbouring segments of the spatial network.


texts: "Poor" quality when zooming out a WMS-Layer based on GeoTiff data

texts:
i am having some issues when using the Geoserver Web Mapping Service to display GeoTIFF-Files. When zooming out (in an OpenLayers-Client or in the Geoserver Layer Preview for instance) the quality seems to be very poor. When opening the GeoTIFF in a picture viewer the quality seems to be ok

I am working on a project where i have to display very old maps in an webbased application, so the quality is quite important.
I have alredy tried to alter the compression levels of the wms in Geoserver, in vain.

I'd be grateful if somebody could come up with a hint or an advise.
Thank you.!


texts: Change detection between binary rasters in QGis

texts:
I have two binary rasters, raster1 and raster2, which represent the distribution of grassland habitats inside an area in two different moments in time (raster1 = older, raster2 = recent). In both rasters, grassland areas have the value 255, while everything else has the No data value of 0. The two rasters have the same extent and the cell centers line up.

I would like to do a simple change detection operation and create a new raster whose values are something like:

if raster1 and raster2 = 255 --> assign 255 to new raster (or any other value, i simply want it to represent "no change and remained grassland"

if raster1 and raster2 = 0 --> 0 (no change and remained not grassland)

if raster1 = 0 and raster2 = 255 --> 100 (changed and became grassland from no grassland)

if raster1 = 255 and raster2 = 0 --> 200 (changed and became not grassland from grassland)

Anyone can suggest a way to do this in QGis? (that doesn't require much knowledge in writing codes, if possible)


texts: Smoothing all the intersection for GPX track points?

texts:
I have a.GPX Track file that has too many intersections. I want to smooth out due to all the intersection
Currently, I'm trying to figure out a way to smooth out this track file. The intersections are not smooth and interfering with my perimeter area of these specific sections. Is there any way to smooth this out or will I have to create a new shapefile from scratch?


texts: No qt5widgets.dll installed is reason GIS files can not be loaded

texts:
Original Problem; some GIS files could not be accessed, because no qt5widgets.dll installed.  I downloaded qt5widgets.dll from DLL-Files.com Client. This crashed QGIS program and I get:

The procedure entry point: ??0QAbstractTextDocumentLayoutPrivate@@QEAA@XZ could not be located in the dynamic link library Qt5Gui.dll.
And
Oops, looks like error loading QGIS. Details: Could not load qgis_app.dll.
Windows Error: the specified procedure could not be found.
Help: check C:\Program Files\QGIS 3.4\bin\qgis-g7.env for correct environmental paths.

I downloaded new qt5Gui.dll, but still same error, even with specific directions to download: Qt5Gui.dll Version 5.11.2.0. (Windows 7 64-bit) to \QGIS 3.4\apps\qt5\bin. Also uninstalled and reinstalled QGIS 3.4.3. No luck.

Is there anything I can do?


texts: polygon vertex to vertex distance calculation in arcobjects

texts:
I have polygon with multipart (donuts) I want to calculate vertex to vertex distance. I tried to get outer ring vertex to vertex distance and it is working. Could you please tell me how to get inner ring vertex to outer ring vertex distance?

IPolygon polygon = buld_F.Shape as IPolygon;
IPointCollection ppcol = polygon as IPointCollection;
                    for (int j = 0; j < ppcol.PointCount -1; j++)
                    {
                        IPoint pp = ppcol.get_Point(j);
                        IPoint pp2 = ppcol.get_Point(j + 1);
                        double xvalue2 = pp2.X;
                        double yvalue2 = pp2.Y;
                        double xvalue1 = pp.X;
                        double yvalue1 = pp.Y;
                        double dist = Math.Sqrt(Math.Pow(xvalue2 - xvalue1, 2) + Math.Pow(yvalue2 - yvalue1, 2));
                        double ptop_dist = dist * 100;
                        if (ptop_dist <= 20)
                        {
                            dt.Rows.Add(buld_F.OID, buld_F.Class.AliasName, tgid_val, xvalue2, yvalue2, "Vertex to Vertex distance < 20CM");
                        }
                    }


texts: Clipping tiff raster image using bounding box with GDAL in python

texts:
I would like to know how can I clip a raster image using a bounding box in python.

So far I can do it with gdal_translate:

gdal_translate -projwin -75.3 5.5 -73.5 3.7 -of GTiff original.tif new.tif

I tried:

import gdal
from gdalconst import *
src_filename = 'ori.tif'
dst_filename = 'new.tif'
src_ds = gdal.Open(src_filename, GA_ReadOnly)
format = "GTiff"
driver = gdal.GetDriverByName( format )
dst_ds = driver.CreateCopy( dst_filename, src_ds, 0, [ '-projwin=-75.3 5.5 -73.5 3.7'] )
dst_ds = None
src_ds = None

I know that the -projwin option should not be defined this way. I get the message:

  Driver GTiff does not support -projwin creation option


texts: How do I change a text box to a text area? (ArcGIS API for JavaScript)

texts:
I recently posted this question on the Esri forums and thought I'd give stack exchange a try too.

I'd like to take this sample:
http://help.arcgis.com/en/webapi/javascript/arcgis/demos/ed/ed_multipleAttrInspector.html

and change the text boxes to text areas similar to the "Editor with Toolbar (Union)" sample.

You can find these samples by going to:
http://help.arcgis.com/en/webapi/javascript/arcgis/help/jssamples_start.htm

The one I am trying to customize is listed under Editing - Edit Modifiable Fields
One of the samples that uses text areas instead of text boxes is under Editing - Editor with Toolbar (Union)

Any help would be greatly appreciated.  If you find a line of code that you think I need to add, please tell me exactly where it needs to be added in the "Edit Modifiable Fields" sample.


texts:
My question is about the best way to convert from NAD 83 to WGS 84 in PostGIS.

I have a PostGIS DB containing National Elevation Data rasters. The rasters are in NAD 83 and I would like to convert them to WGS 84 but PostGIS treats NAD 83 and WGS 84 as the same (see here and here) even though they may differ by over a meter.

I figure that I can just add a new spatial reference identifier with non-zero values for the transformation to WGS 84. So basically this below except with TOWGS84[0,0,0,0,0,0,0] replaced with the correct values.

INSERT INTO "spatial_ref_sys" ("srid","auth_name","auth_srid","srtext","proj4text") VALUES (9999,'EPSG',4269,'GEOGCS["NAD83",DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6269"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4269"]]','+proj=longlat +datum=NAD83 +no_defs ');

The problem is that I can't find what those values should be, as in both what the values should refer to and what the value of the values are.


texts: Conversion of Sirgas2000 projection system to wgs84

texts: I am having a bunch of .csv files containing "x" and "y" of some locations (to be specific, Brazil and the coordinate system is SIRGAS 2000). I want to convert them into wgs84 i.e, lat/long using arcpy and rewrite the .csv files having lat/longs to corresponding "x" and"y" using python. Can anyone please help me doing that?


texts: Read a .csv file from a URL with an API call in QGIS 3?

texts:
How do I read a .csv file from a URL with an API call, including key, in QGIS 3?

I've tried creating a .vrt file as outlined here

I've tried the solutions presented here

If I'm doing everything correctly, then the differences between these examples and my case are that my URL is https and that my URL has an API call following the .csv extension.

eg:
https://developer.nrel.gov/api/alt-fuel-stations/v1.csv?fuel_type=HY&state=CA&api_key=G3tUr0wNk33&format=csv

The .vrt file I assembled:

<OGRVRTDataSource>
<OGRVRTLayer name="test">
    <SrcDataSource relativeToVRT="0">/vsicurl/https://developer.nrel.gov/api/alt-fuel-stations/v1.csv?fuel_type=HY&state=CA&api_key=G3tUr0wNk33&format=csv</SrcDataSource>
    <GeometryType>wkbPoint</GeometryType>
    <LayerSRS>WGS84</LayerSRS>
    <GeometryField encoding="PointFromColumns" x="Longitude" y="Latitude"/>
</OGRVRTLayer></OGRVRTDataSource>


texts: Getting (0,0) instead of the intersecting co-ordinates value

texts:
This code works below fine when line(The Circle I mean) intersects more than two new lines but fails when both the intersection falls in the middle of the same line

polygon_path = iface.activeLayer()
        line_path = buffers_dir+"/bufferline"+ str(i.id()) +".shp"
        #Converting the buffer features to lines
        processing.runandload("qgis:polygonstolines",polygon_path,line_path)
        linebuffer=iface.activeLayer()
        selections=[]  # The id of he intersected path of the network layer
        for f in layer.getFeatures():
            for a in linebuffer.getFeatures():
                print(a,'This is A')
                if a.geometry().intersects(f.geometry()):
                    #Finds the intersection of line buffers with network layer
                    intersection = a.geometry().intersection(f.geometry())
                    #same result with .within() and even with a and f switched (the docs aren't that clear on that)
                    print(intersection.asPoint())
                    instersection_point = intersection.asPoint()
                    connecting_points.append(instersection_point) #Contains the (X,Y) co-ordinates of intersection points
                    selections.append( f['WVK_ID'] )

Expected result is something like this:

((114041,400590), ' This is intersection')
([228201037.0, 227201083.0, 228201038.0, 228200076.0, 228201026.0, 227200047.0, 227201082.0], ' This is Selections')

The result i got is:

((0,0), ' This is intersection')
([227200062.0], ' This is Selections')

Error Occured:

  File "<string>", line 186, in <module>
  File "C:\PROGRA~1\QGISWI~1\apps\Python27\lib\random.py", line 274, in choice
    return seq[int(self.random() * len(seq))]  # raises IndexError if seq is empty
IndexError: list index out of range


texts: How to georeference HDF5 file with missing or wrong metadata

texts:
I'm trying to access data in a HDF5 file. I can access the data itself but the metadata, in particular the information to georeference the information, is still a problem. I had partial success in this regard. But the projection claimed in the metadata has to be wrong, because this is data for the Netherlands, not the North pole. I have contacted the provider pointing this out, but they have ignored my point, just briefly answering me with a suggestion to check the product documentation (which I translate with Google Translator).
In the documentation, they just mention two pieces of information regarding georeference:

the size of the image in pixels (which is trivial)
each pixel is 1 km x 1 km
they give the coordinates of the four corners of the image

They do not mention the projection nor any affine transform. Am I missing something? Is this information sufficient to georeference this data? I think not, but I might be misunderstanding this.
Using QGIS I have put these four points in a map and they do not line-up:

Well, the problem is that I do not know the projection. I have tried Mercator (the one used to take this snapshot), some conic and an stereographic (common for radar information). But they still do not line up. Something I find very odd that the points have so different latitudes. I would expect these points to form something very similar to a square, but it seems to be rotated.
Is it possible to georeference this data with the sparse information given in the documentation? It seems currently impossible, and I do not think the providers will help me any further.


texts: Create polygons with blurred edges (SLD (GeoServer))

texts:
I would like to create areas from polygons that will visualize an "uncertain" area. So I would like the polygon to appear a little blurred at the edges. Preferably with SLD in GeoServer so that I could get the layer with WMS. I have OpenLayers in frontend so that could be an option as well if its possible that way..
Does anyone know if this can be done?


texts: I was following a tutorial on time series which showed a screen grab of the properties dialogue with "Temporal" as a tab.  When I folowed the tutorial using the same file Temporal is missing.  Do I need to do anyhing to activate Temporal?  I have reloaded Qgis but still no Temporal.


texts: What column should I use as primary key after importing S-57 into PostgreSQL?

texts:
I've imported data from S-57 format into a PostgreSQL/PostGIS database using ogr2ogr (GDAL 2.1.0) and it resulted in several tables; one for each feature type. I intend to merge these tables into a single table.

However, the current primary keys for the tables are ones which I think was introduced by the importer, ogc_fid, as the values count up from 1 with increments of 1 for each table. That means that the values would not be unique after a merge and can thus not be used as primary key. I intend to drop that column after or before the merge as it doesn't fit my needs.

That leaves three columns which are candidates for being used as primary key, and I'm here to ask which one to pick as I'm not very familiar with the S-57 format.

RCID: This column is not listed here. I've not found information on this column in many places (Googling didn't return many relevant results), though one place referred to it as something like "Record Identification". From what I can see, it contains the same value as FIDN, though I can't say for sure that so is the case for all features.

FIDN: Feature Identification Number. It stands to reason that this one would be unique for each feature, but I'm not so sure. Is it possible that a feature could fall under multiple feature types and thus the same FIDN show up in multiple tables? I'm not familiar enough with S-57 to know.

LNAM: Long Name. It's a combination of AGEN, FIDN and FIDS. So far all features I've seen have had the same values in AGEN and FIDS, so I don't see how this one would be more unique than FIDN if the same FIDN value could appear in multiple tables.

If I stay in doubt I'm probably going to make some combination of FIDN and OBJL into the primary key, because every table have different OBJL values (because it identifies the feature type) and the FIDN should at least be unique within each table.

EDIT: Upon closer inspection I see that the OBJL is already included in the FIDN. So they are already unique between tables.


texts:
I have to tables inside my database (postgis). Now I will do based on QGIS execute a sql-query based on these two tables. How can I do this?
There exists a tool "Execute SQL" but there I can only choose one layer. I will realize the query inside my database in postgis.

I have a Processing model. The first step was to load the two tables inside my database. The next step is to execute the sql query. How can I do this?


texts: Geonames querying for parks in Puerto Rico or US Virgin Islands

texts:
I thought these were US-related entities.

When I make this kind of a query to geoNames, I get nothing back.

http://api.geonames.org/search?featureCode=PRK&username=genadinik&country=US&style=full&adminCode1=CR&maxRows=10

Does anyone know how to query for parks in these territories?


texts: Tile server - problem with tiles loading (pink tiles)

texts:
I have built a tile server by following this instruction: https://switch2osm.org/serving-tiles/building-a-tile-server-from-packages/ and i have a problem.
When i test tile server using http://localhost/osm/slippymap.html, some tiles are pink. After zoom out and zoom in tiles are usually properly displayed.
I tried to run renderd in debugging mode. This is the sample output:

renderd[19572]: DEBUG: Got incoming request with protocol version 2
renderd[19572]: DEBUG: Got command RenderPrio fd(12) xml(default), z(11), x(1156), y(691), mime(image/png), options()
renderd[19572]: DEBUG: Got incoming connection, fd 13, number 6
renderd[19572]: DEBUG: Got incoming request with protocol version 2
renderd[19572]: DEBUG: Got command RenderPrio fd(13) xml(default), z(11), x(1158), y(691), mime(image/png), options()
renderd[19572]: DEBUG: DONE TILE default 11 1160-1167 688-695 in 6.073 seconds
debug: Creating and writing a metatile to /var/lib/mod_tile/default/11/0/0/66/139/128.meta

renderd[19572]: DEBUG: Sending render cmd(3 default 11/1160/689) with protocol version 2 to fd 8
renderd[19572]: DEBUG: Sending render cmd(3 default 11/1160/690) with protocol version 2 to fd 10
renderd[19572]: DEBUG: Sending render cmd(3 default 11/1160/688) with protocol version 2 to fd 9
renderd[19572]: DEBUG: Failed to read cmd on fd 8
renderd[19572]: DEBUG: Connection 0, fd 8 closed, now 5 left
renderd[19572]: DEBUG: Failed to read cmd on fd 10
renderd[19572]: DEBUG: Connection 1, fd 10 closed, now 4 left
renderd[19572]: DEBUG: Failed to read cmd on fd 9
renderd[19572]: DEBUG: Connection 0, fd 9 closed, now 3 left
renderd[19572]: DEBUG: Got incoming connection, fd 8, number 4
renderd[19572]: DEBUG: Got incoming request with protocol version 2
renderd[19572]: DEBUG: Got command RenderPrio fd(8) xml(default), z(11), x(1155), y(691), mime(image/png), options()
renderd[19572]: DEBUG: Got incoming connection, fd 9, number 5
renderd[19572]: DEBUG: Got incoming request with protocol version 2
renderd[19572]: DEBUG: Got command RenderPrio fd(9) xml(default), z(11), x(1159), y(691), mime(image/png), options()
renderd[19572]: DEBUG: Got incoming connection, fd 10, number 6
renderd[19572]: DEBUG: Got incoming request with protocol version 2
renderd[19572]: DEBUG: Got command RenderPrio fd(10) xml(default), z(11), x(1154), y(691), mime(image/png), options()
renderd[19572]: DEBUG: Failed to read cmd on fd 11
renderd[19572]: DEBUG: Connection 0, fd 11 closed, now 5 left
renderd[19572]: DEBUG: Failed to read cmd on fd 12
renderd[19572]: DEBUG: Connection 0, fd 12 closed, now 4 left
renderd[19572]: DEBUG: Failed to read cmd on fd 13
renderd[19572]: DEBUG: Connection 0, fd 13 closed, now 3 left
renderd[19572]: DEBUG: Failed to read cmd on fd 8
renderd[19572]: DEBUG: Connection 0, fd 8 closed, now 2 left
renderd[19572]: DEBUG: Failed to read cmd on fd 9
renderd[19572]: DEBUG: Connection 0, fd 9 closed, now 1 left

I've also tried to pre-render tiles for region. Command for pre-rendering is for example:

render_list -all -s /var/run/renderd/renderd.sock -x 34 -X 36 -y 20 -Y 21 -z 6 -Z 6

Pre-rendering process displays messages pasted below and i don't know if it is working still or not:

debug: init_storage_backend: initialising file storage backend at: /var/lib/mod_tile
Rendering client
Starting 1 rendering threads
Rendering all tiles from zoom 6 to zoom 6
Rendering all tiles for zoom 6 from (34, 20) to (36, 21)
Waiting for rendering threads to finish

I can't see new tiles in /var/lib/mod_tile after that process.

Why does it happen (pink tiles)? How can I prevent this?


texts: Can single user have multiple roles in ArcGIS?

texts:
Is it possible to assign multiple roles to a single user?
I would like to have the following setup:

Role A gives edit/view access to layer A
Role B gives edit/view access to layer B

and

User X has Role A and Role B
User Z has Role A
User Q has Role B


texts:
I need to combine several rasters. One of them has an alphanumeric data in the attribute I need to combine. However, when combined the output raster has a 'LONG' attribute type where should be 'STRING' (I believe it's due the first row in that attribute is a number, not an alphanumeric character, so the combine tool assign a 'LONG' type instead of a 'STRING').

How can I create a sorted new raster, where characters go first and then numbers in the mentioned attribute?

I've tried:

sort the raster (Descending), then export data to a new raster. Didn't work.
sort the raster (D), then lookup tool. Didn't work.

How can I sort a raster attribute, and save the sorted changes?


texts: Is there a way to create an urban extent/population distribution map which only includes areas with a population above 100,000

texts:
I'm looking to do analysis on cities throughout the world. I'm using WorldPop, LandScan and MODIS data to gain urban extent and population density maps. I would like to extract areas which have a population size above 100,000. At the moment I think I may have to manually go through each city which has a population of above 100,000 and clip around this to create a new layer, is there another way around this?
I'm using ArcGIS 10.2.2.


texts: Making darkness/intensity of overlapping buffers additive?

texts:
I don't know much about GIS, but I want to create some maps that show a particular thing.

Specifically, I have a shapefile with several hundred points and I want to create a standard sized buffer around each point (for example, 5km in radius). Secondly it is very important that the buffers for all points are the same color, but vary in intensity, with some being very light and some being very dark. The intensity will depend on a certain attribute of the points.

So my first question is, is it possible to shade the buffers on a continuous gradient, or do I have to use classes?

For example, if the attribute for the points is on a scale from 0 to 1,000, do I have to divide that into classes?

Or can the shading just key directly on the value of the attribute?

Here's my second question: many of the buffers will be overlapping and it is very important that if they overlap that they appear even darker in the overlapping area.

Is that possible?

To be clear, I'm not trying to do any analysis with this job. I just want visualizations.

Alternatively, if that is too difficult, I was thinking of creating a heat map (instead of a map with fixed ring buffers of varying intensity). But in order for the heat map to represent that particular attribute of each point, I would have to multiply the number of points, based on the value of the attribute. For example, if the attribute value for a point was 100, I would want to make 99 more points in that same location. In that way, the number of points at each location would represent the attribute value of each point in the original set of points. Then a heat map would show the density of the points.

But is it difficult to create a script that multiplies points based on an attribute value?


texts:
I'm trying to find a .NET GeoJSON library which I can use to export my SqlServer2008 spatial boundaries (in binary format -> eg. STAsBinary()) from my .NET code to my HTML frontend to display on a google map.

At first I felt that it might be a StackOverflow question, but this site I feel is more specialised to helping me out.

Does anyone know of any?

NOTE: This dll does NOT have to be imported into SqlServer2008. In fact, it will not. I'll use it in my C# MVC WebSite.


texts: What to do with completely bad geometry?

texts:
I have downloaded a shapefile that I need from a website. However, the geographic coordinate system is in a projected coordinate system WGS84 Major Auxiliary Sphere and it is projected in WGS84 Web Mercator (picture attached).
Does anyone know of a way to preserve the polygons within the file but convert this to WGS84 or a geographic coordinate system with only one projection so the polygons are useful?
I have also tried reprojecting to WGS84 (transformation WGS_1984_Major_Auxilary_Sphere_To_WGS_1984) but I get an error saying

ERROR 000365: Invalid geographic transformation.

I understand that I should not define the projection since it is not what the file was originally created in, but did try it just in case and the polygons within the file disappeared. Does anyone know how to get this polygon file into the WGS84 geographic coordinate system? I would prefer it to not have a projection.
Is this problem solvable or do I need a new version from the creator?


texts: How to export tracks of a specific polygon in overpass-turbo?

texts:
I would like to export some features (in this case route=hiking) that are inside a specific boundary. So not the boundary box that I can choose through the pan tool, but all the routes that for example belong to a national park.

Can I do this with the wizard?

I've founded this example in the OSM wiki. I'm trying to perform something like that.


texts: Extract vectors from ArcGIS MapServer into virtual layer with QGIS?

texts:
This is the first time I've used QGIS, so bear with me. I have an ArcGISMapServer data source and added one of its layers to the canvas. It doesn't let me customize any of the styles, because I assume QGIS is just pulling down tile images. All of the features have geometry data, so I'd like to extract the polygons into a virtual layer so that I can style them myself. But I'm not sure what query to use.

Is this even possible and, if so, how would I do it?


texts:
I have a problem loading my shapefile named watershed and link it to my table in postgresql named watershed as mentioned here:

shp2pgsql -I -s 32632 C:\Users\ASUS\Desktop\master_PEER\DATABASE\Shapefile\watershed.shp watershed | psql -U postgres -d Chiba_watershed.

  shp2pgsql : Le terme «shp2pgsql» n'est pas reconnu comme nom d'applet de commande, fonction, fichier de script ou
  programme exécutable. Vérifiez l'orthographe du nom, ou si un chemin d'accès existe, vérifiez que le chemin d'accès
  est correct et réessayez.
  Au caractère Ligne:1 : 1
  + shp2pgsql -I -s 32632 C:\Users\ASUS\Desktop\master_PEER\DATABASE\Shap ...
  + ~~~~~~~~~
      + CategoryInfo          : ObjectNotFound: (shp2pgsql:String) [], CommandNotFoundException
      + FullyQualifiedErrorId : CommandNotFoundException

][2]

[2]: https://i.stack.imgur.com/jETfb.jpgemphasized text

How to resolve this problem?


texts: Changing feature from one feature class to another using ArcGIS Pro?

texts:
I am trying to change a feature from one feature class to another (not copy)

For example, I am working with a town's water system.  We have mainline valves and hydrant valves.  I have come across a number of mainline valves that have accidentally been put in as hydrant valves.  How would I take these valves and change them from hydrant valves to mainline valves.  Another example, we have a feature class for water mains and a feature class for abandoned water mains.  If we recently abandoned a main, how would I move that main from the water main feature class into the abandoned water main feature class? Once last example, we are working with sewer mains and water mains side by side.  A main was attributed to the sewer main feature class, but it is in fact a water main.

How would I take that line and just convert it from a sewer main to a water main?

I am using ArcGIS Pro.


texts: Different output of gdalinfo --version and gdal-config in Xenial

texts:
I have Ubuntu 16.04 with postgresql 9.6 postgis 2.3 QGIS 2.18 installed.

I am trying to install mapview (R library) with no success because of sf defencency.

I had this output:

    configure: GDAL: 1.11.3
checking GDAL version >= 2.0.0... no
configure: error: sf is not compatible with GDAL versions below 2.0.0
ERROR: configuration failed for package ‘sf’
* removing ‘/usr/local/lib/R/site-library/sf’

and this is what happened on the terminal:

$gdalinfo --version
GDAL 2.1.0, released 2016/04/25
$ gdal-config --version
1.11.3

I am wondering the following questions:

Which is the difference between gdalinfo and gdal-config ?

Why does R say me that I have a GDAL version less than 2.0.0 even if output of gdalinfo --version is GDAL 2.1.0?

How can I upgrade Gdal-config to one >= 2.0.0?

If I upgrade, will I have problem with postgresql, postgis Qgis or some other python module that use GDAL?


texts:
I have a layer from one source that was originally in one coordinate system and am attempting to use it with map layers from Google and Open Maps.  I have set the project properties, the map layer and the shapefile all to EPSG 42105 - WGS84.

Is there some other setting or something?

The shapefile is in North America and when I zoom to the layer I see it in an ocean somewhere.    It is still blue water as I zoom to the where I can no longer see the shapefile.


texts:
I have a mbx that creates a layout. The layout consist of a number of .wmf pictures and ofcource a map. During the execution I create .Tab files that reference to the pictures looking like:

Definition Table
  File "picture.wmf"
  Type "RASTER"
  (0,1) (0,0) Label "Pkt 1",
  (1,1) (1,0) Label "Pkt 2",
  (0,0) (0,1) Label "Pkt 3"
  CoordSys NonEarth Units "ft"
  Units "ft"

When something is changed and I make a new layout some pictures are updated and overwriten with new ones. The .Tab files is still the same. The problem is that it's the old pictures that shows up in the layout. When I open the pictures direkt in my filesystem it shows the new one. I think MapInfo caches the old picture in some way and wonder if there is a way to reload the .Tab file or in some way to make MapInfo understand that the information that the .tab file refers to is updated?


texts: Adjust level of detail of Openstreetmap on specific zoom

texts:
I need an overview of a specific large area on QGIS Desktop 2.16.0.
The details disappears with this level of zoom.

Is it possible to control the display level ?

If I zoom too much, I can't see all the concerned area.

If I unzoom, I lose details.

From the answer I was able to show details:

However is it possible to get this bigger and clearer like this ?


texts: Creating weighted list using Select tool in ModelBuilder?

texts:
I've done a raster to point conversion and have a column of the grid codes.

Let's say it is composed of 10 numbers [9,5,7,1,7,3,7,2,8,1].

Right now I am doing a selection on those numbers based on the top 50% of the Max value(Select values greater than 4.5, which gives {9,8,7,7,7,5})

What I would like to do is take the original array and select all of the values that make up greater than 50% of the total(Total is 68/2 = 34, so SUM(9,8,7,7,7) > 34) and have those as my selection points.

I think I can sort the data and use another 1 or 2 calculate fields to get the sum of the column and the 50% value, but how would I go about adding the values to get above 50% and save them as the selection?


texts: How to access both IIS and Geoserver at port 80 from outside the production environment?

texts:
At production server, IIS and GeoServer are installed. I can access IIS at port 80 From out side the production environment. i also want to access GeoServer at the same port 80 from out side of production environment, Like http://92.108.64.13:80/geoserver/web/.
But, i could not access GeoServer at port 80.

Is it possible to access both IIS and Geoserver at port 80 ? Actually I do not want to open any other port for accessing GeoServer. Is there any solution available for accessing both Geoserver and IIS at port 80 ?

need for useful suggestion for solving this problem


texts:
I am using QGIS 2.0.1-Dufouron on Windows 7. When I try to load nviz via GRASS, I get the message "Module nviz not found". Weirdly enough, I can load it via the toolbox and be able to 'run' it (it reaches 100%). But afterwards, nothing happens.

Problem: I merely want to rotate and play around with the view points of layers. This is all for aesthetic purposes, don't want to do anything technical (not yet anyway).

I have read posts on people suggesting to install the plugin SEXTANTE which may solve the issue, however it doesn't show up on my list of plugins (maybe I don't have the correct URL?). I am considering trying ArcMap if I have to but would rather do it in QGIS.

Any suggestions please?


texts: How to search for layers in the layers panel?

texts: I have QGIS projects with a lot of layers. Sometime it's difficult to find a layer into Layers Panel. How to search it?


texts:
I'd like to update records of a GeoPackage containing vector information in R. I've tried this:

library(sf)
library(RSQLite)

nc <- st_read(system.file("shape/nc.shp", package="sf"))

# Write Geo Package
gpkg <- "nc.gpkg"
st_write(nc, gpkg)

# Create connection
con <- dbConnect(SQLite(), dbname = gpkg)

# Update record
res <- dbSendQuery(con, "UPDATE nc SET NAME = 'My new name' WHERE CNTY_ID == 1828")

Which gives me the following error:

Error: no such function: ST_IsEmpty

A bit of digging reveals that ST_IsEmpty is one of several SQL functions that are required for updating a GeoPackage. How would I go about including these? Or are there any other methods of performing this operation?


texts: I'm trying to apply WPS for local shapefile.  Is it possible to process local shapefile directly as a input data?


texts:
Using ArcGIS Desktop 10.7.
I would like to customize the World Terrain Basemap in order to select which roads and topography labels appear.
Is there a way to do this, and if not, then how would I create my own hillshade terrain basemap with a selection of topography features?


texts: Creating voronoi polygons from polygons?

texts:
I have a database with buildings and for some further analysis I need to:

aggregate adjacent buildings (preferably without losing the height attribute, but I can work my way around that)
once I have them aggregated, I need to create 'allotment' polygons so basically the Thiessen Polygons based on them - so that they do not intersect an any point aggregated buildings.

I tried to create centroid out of those aggregated buildings and then creating Thiessen Polygons but I failed - they intersect the buildings.
Any idea how to handle it?

Update after fisrt comment


texts: Method for identifying upstream segments

texts:
I need to add to the attribute table of a stream network layer a column that can identify all upstream segments that contribute to the segment (for each segment).

My best idea at the moment is to assign all stream segments with a Strahler order of 1 a prime number.  Then going downstream, I multiply the contributing upstream segments (see graphic below).  From this, I think I can factorise the number into it's primes and then work out the IDs of all upstream segments.

I already have a script that is capable of passing information downstream based on the Strahler stream order so this is not a problem.  I just want to know what a better method would be, my method is going to generate horribly large numbers once I start processing large stream networks!

As an example, in the graphic below, if i was to look at the segment with the 210 attribute, I would factorise this into primes of 2,3,5 and 7.  From these primes I could work out that 6 and 35 were also upstream segements.

I have access to ArcGIS basic 10.2.2 and QGiS.


texts: When importing shapefile gives an error "No dataset found at path"

texts:
I'm trying to follow this tutorial to GeoPandas and am trying to import the "oceans.shp" file, which I downloaded from here (together with the .dbf, .prj, qpj, and .shx files),
When running the code
oceans = gpd.read_file('[...]\oceans.shp')

I get the error
FionaValueError: No dataset found at path '[...]\oceans.shp' using drivers: *

What am I doing wrong - can someone point me into the right direction?


texts: Setting up IGeoProcessor2 object in .NET (C#)?

texts:
I'm wondering if I'm instantiating a new class with the IGeoProcessor2 object correctly. It's working the way I have it now, but I feel like there's something else I should be doing so I don't have to turn "Embed Interop Types" in the Properties window to False. I'm working with ArcObjects 10.5, ArcGIS 10.5, and Visual Studio 2015.

This is the code I have to instantiate a new geoprocessing object:

IGeoProcessor2 gp = new GeoProcessorClass();

If I leave the Esri.ArcGIS.Geoprocessing assembly's "Embed Interop Types" property on "True", I get a squiggly red line under GeoProcessorClass. If I set it to "False", the squiggly goes away and everything works as expected. The same thing also happens with the ESRI.ArcGIS.System assembly when I use this code:

IVariantArray parameters = new VarArrayClass();

If I were to set the "Embed Interop Types" to false for the System assembly then the squiggly under VarArrayClass goes away and everything works great.

Is this something everyone has to do to get their code to work, or am I missing something that allows me to keep "Embed Interop Types" as "True"?

I am using IGeoProcessor2 instead of GeoProcessor because I am running a custom tool out of a custom toolbox, so I cannot switch to GeoProcessor. For whatever reason, I cannot get it to work with my code.


texts: Problem running geotools quickstart application?

texts:
package org.geotools.tutorial.quickstart;

import java.io.File;

import org.geotools.data.FileDataStore;
import org.geotools.data.FileDataStoreFinder;
import org.geotools.data.simple.SimpleFeatureSource;
import org.geotools.map.FeatureLayer;
import org.geotools.map.Layer;
import org.geotools.map.MapContent;
import org.geotools.styling.SLD;
import org.geotools.styling.Style;
import org.geotools.swing.JMapFrame;
import org.geotools.swing.data.JFileDataStoreChooser;

/**
 * Prompts the user for a shapefile and displays the contents on the screen in a map frame.
 * <p>
 * This is the GeoTools Quickstart application used in documentationa and tutorials. *
 */
public class Quickstart {

    /**
     * GeoTools Quickstart demo application. Prompts the user for a shapefile and displays its
     * contents on the screen in a map frame
     */
    public static void main(String[] args) throws Exception {
        // display a data store file chooser dialog for shapefiles
        File file = JFileDataStoreChooser.showOpenFile("shp", null);
        if (file == null) {
            return;
        }

        FileDataStore store = FileDataStoreFinder.getDataStore(file);
        SimpleFeatureSource featureSource = store.getFeatureSource();

        // Create a map content and add our shapefile to it
        MapContent map = new MapContent();
        map.setTitle("Quickstart");

        Style style = SLD.createSimpleStyle(featureSource.getSchema());
        Layer layer = new FeatureLayer(featureSource, style);
        map.addLayer(layer);

        // Now display the map
        JMapFrame.showMap(map);
    }
}

  The last line "JMapFrame.showMap(map)";  required:  MapContext   found:
  MapContent

Then I changed

MapContent map = new MapContent();

To

 MapContext map = new MapContext();

Then the application starts but when I select shape file

  Exception in thread "main" java.lang.NoClassDefFoundError:
  javax/media/jai/PropertySource

Am I missing any library?
I have installed both jai-1_1_3_lib-windows and jai_imageio

Exception in thread "main" java.lang.NoClassDefFoundError: javax/media/jai/PropertySource
    at java.lang.ClassLoader.defineClass1(Native Method)
    at java.lang.ClassLoader.defineClass(ClassLoader.java:791)
    at java.security.SecureClassLoader.defineClass(SecureClassLoader.java:142)
    at java.net.URLClassLoader.defineClass(URLClassLoader.java:449)
    at java.net.URLClassLoader.access$100(URLClassLoader.java:71)
    at java.net.URLClassLoader$1.run(URLClassLoader.java:361)
    at java.net.URLClassLoader$1.run(URLClassLoader.java:355)
    at java.security.AccessController.doPrivileged(Native Method)
    at java.net.URLClassLoader.findClass(URLClassLoader.java:354)
    at java.lang.ClassLoader.loadClass(ClassLoader.java:423)


texts:
in QGIS I have to draw a scheme ( kind of CAD) that represents electrical circuits.
The scheme consists of points and lines elements, with alphanumeric attributes associated.

The starting point and the ending point are connected by 1 or more lines, drawn with snap (there are intermediate points).
For example:

CIRCUIT A:
 X--------------------------X-------------X---------X
P1          L1              P2    L2      P3   L3   P4

CIRCUIT B:
 X-------------------X-------------X
P5          L4       P6    L5      P7

This is an examples of my net:

Having selected the starting point only, I would like to determine what is the set of lines that make the connection and eventually what is the ending point.

In QGIS are there "native" network analysis / routing functions or plug-ins for network analysis?

I currently use QGIS 3.0 but I could eventually install a previous version if needed.


texts: Does QGIS Server support QGIS project with Oracle layers?

texts: QGIS recently got native Oracle support. I was wondering if QGIS Server supports generating WMS and WFS from a QGIS project with Oracle Spatial layers. I just tested this but I did not succeed...


texts: When importing a TIFF image into Google Earth Engine, I can't view it as it's completely black

texts:
I'm trying to import a TIFF (Landcover) file into my GEE working folder.
The image I get is completely black and without information about the floor coverings.


texts: Understand gdalwarp and gdalinfo on raster files


texts: Shape comparison language (SCL) - get shapes which intersect extent border

texts:
I am trying to write a SCL string for my ArcGIS query, which should return only shapes which are intersected by a border of a rectangular shape (extent). So far the closest I've got is the following:

"RELATE(G2,G1,'G1 INTERSECT G2.BOUNDARY')"

However this also returns shapes which are fully contained within the rectangular shape and don't intersect the border. Maybe someone with SCL experience could provide a few tips?


texts: Halting QGIS on startup until user is authenticated?

texts:
I have a custom made login form in Python, which launches on QGIS startup. But before the user has a chance to enter the credentials, the program starts normally, bypassing the form. From then on the user can interact with QGIS like the form isn't even there (which defeats the purpose of the form!! ).

Is there a way to halt QGIS startup until the user is authenticated?


texts: Unsplitting Zpolylines without flattening vertical line segments?

texts:
Trying to figure out how to fully keep vertical line segments when using Dissolve/unspit lines...

The first image shows the original Zlines in Arcscene

The second image is the result of a Dissolve /Unsplit lines

Vertical Segment (identical x, y but different z) are missing, the vertex are deleted and not simply stacked.

This is a known Problem (already identified with ArcGis10.0 but apparently unsolved) : http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//002t00000010000000

Looking for some workaround, I found an earlier Q&A titled Vertical lines in a polygon shapefile. It could be one of the step of a larger process.

Has anyone had to deal with dissolving/unsplitting vertical Segment?

Using ArcGis Advanced 10.2 + full Version of Edit Tools, FME, QGIS (AutoCAD is available, though no expert user on site).


texts: QGIS 2.2: Rotated labels not possible any more?

texts:
Using Qgis Valmiera I would like to create feature labels on a layer that are rotated with a rectancular whitish background.

Using "Labels (deprecated)" and entering an angle does the first part of it but I can't make a rectangular backgound. The "buffer labels" option does something similar but it's too fancy and creates giant PDFs.

Using "Labels" I can easily draw a rectangular background but it seems the option to rotate text has been deprecated. Interestingly I can rotate the backgound frame ...

Is there a way to combine both?


texts: Not able to create PostgreSQL file using Osm2po?

texts:
I am new to the GIS community.

After downloading the osm2po jar file, I have run following the instructions from underdark tutorial. Also, from Carsten's answer in gis stack exchange, I have enabled instructions in osm2po.config file but the sql file is not generated in the prefix folder.

I did not get any error while starting the web server.


texts:
I have two point clouds (PCs): (i) PC-A (i.e. Reference PC or ground truth) and
(ii) PC-B (i.e. Bare-Earth PC obtained through filtering)

I want to subtract PC-B from PC-A to calculate commission and omission errors. So, far I tried computing cloud to cloud distances in cloudcompare, and then filtering points that are not within 0.05m threshold distance. However, since the topography is high variable, the aforementioned method is giving a huge spatial variation. I suspect that cloud-to-cloud distance isn't very accurate.

Is there any tool (preferably opensource and not too complex) that I could use to achieve PC-B (-) PC-A?


texts: Stack of time-spatial maps is not read by readmap method

texts:
I am new in PCRaster, so I made a primitive model to test the functionality. I need to read input from precipitation data maps for each day. My code is:

class Model(DynamicModel):
    def __init__(self):
        DynamicModel.__init__(self)
        setclone(src.RasterYSize,src.RasterXSize,src.GetGeoTransform()[1],src.GetGeoTransform()[0],src.GetGeoTransform()[3])
    def initial(self):
        self.dem=self.readmap('dem')
        self.snowmaxTss=TimeoutputTimeseries("snowmax", self, None, noHeader=False)
    def dynamic(self):
        self.snow=self.readmap('/snowdata/snow')
        maxsnow=max(self.snow)
        self.snowmaxTss.sample(maxsnow)
myModel=Model()
dynModel = DynamicFramework(myModel, lastTimeStep=16, firstTimestep=1)
dynModel.run()

My snow maps stack is correct and corresponds to setclone. But I get an error:

RuntimeError                              Traceback (most recent call last)
<ipython-input-188-c7810b190e43> in <module>
     14 myModel=Model()
     15 dynModel = DynamicFramework(myModel, lastTimeStep=16, firstTimestep=1)
---> 16 dynModel.run()

/opt/conda/lib/python3.7/site-packages/pcraster/framework/dynamicFramework.py in run(self)
     63       self._runInitial()
     64
---> 65     self._runDynamic()
     66
     67     # Only execute this section while running filter frameworks.

/opt/conda/lib/python3.7/site-packages/pcraster/framework/frameworkBase.py in _runDynamic(self)
    369         self._incrementIndentLevel()
    370         self._traceIn("dynamic")
--> 371         self._userModel().dynamic()
    372         self._traceOut("dynamic")
    373         self._decrementIndentLevel()

<ipython-input-188-c7810b190e43> in dynamic(self)
      9         self.snowmaxTss=TimeoutputTimeseries("snowmax", self, None, noHeader=False)
     10     def dynamic(self):
---> 11         self.snow=self.readmap('/snowdata/snow')
     12         maxsnow=max(self.snow)
     13         self.snowmaxTss.sample(maxsnow)

/opt/conda/lib/python3.7/site-packages/pcraster/framework/dynamicPCRasterBase.py in readmap(self, name, style)
     78       `style` argument is not used.
     79     """
---> 80     return self._readmapNew(name)
     81
     82   def _setNrTimeSteps(self,

/opt/conda/lib/python3.7/site-packages/pcraster/framework/frameworkBase.py in __call__(self, *args, **kwargs)
     35         if obj:
     36             attr = getattr(obj, self.callback_attr)
---> 37             result = attr(*args, **kwargs)
     38         else:
     39             result = self.default_callback(*args, **kwargs)

/opt/conda/lib/python3.7/site-packages/pcraster/framework/frameworkBase.py in _readmapNew(self, name, style)
    564     assert path is not ""
    565     import pcraster
--> 566     return pcraster.readmap(path)
    567
    568   def _assertAndThrow(self,

RuntimeError: Raster /snowdata/snow0000.001: can not be opened

I tried different formats of file names: 0000.001, 0000.001.map etc - but result is the same. What's the problem with my code? I know, that readmap() method must read stacks by time.
Is an xml file read when reading a map file? I convert from nc layers to map by GDAL and my xml files may contain odd information from the whole nc file.


texts: PDF export from QGIS print layout: line fill pattern vectors shown with pixelation

texts:
I have created grayscale piecharts in QGIS 3.14 Pi on Win10 with geometry generator. So the piecharts are vectors, created with the expression editor, using the wedge_buffer expression. I created a print layout and exported it to a pdf. I opened this pdf in a vector graphic software (Affinity Designer). All elements are separately saved, can be selected and changed accordingly. When I zoom in, most elements behave as expected: as vectors, they don't show a pixel structure, even small dots are shown as perfect circles. A few lines, however, are shown with a pixel structure, see on the left of this screenshot:

I have also tried to export the layout as SVG, but it seems the structure is too complex and the exported image is full of errors. So I have to export as pdf. The problem is not related to Affinity Designer as I have the same effect when opening the exported pdf with InkScape or with a pdf-Viewer. This is strange as I created this pattern with the geometry generator and a line fill pattern, see the next screenshot:

It seems that only lines created with the line fill pattern behave in this way, regardless of the angle. Other straight lines (e.g. borders of the piechart's slices) are rendered correctly and can be zoomed in without showing a pixel pattern.
So the question is: is there a way to avoid this pixelation of vector layers in an exported pdf?


texts: Alignment of the contents of a legend in MapInfo

texts:
Is it possible to change the alignment of the numbers of the ranges in a legend in MapInfo? (without having to "tease" the layout with a text box on top as in "MapInfo legend editing for right to left text" )

I get this:

... where I would like to get this:


texts: Changing Raster Layer Resampling with PyQGIS?

texts: How can I set the zoomedInResampler of a Raster Layer to Cubic and the zoomedOutResampler to Average using PyQGIS?


texts: Is there a way to transform a node OSM point into the reference GPS system of a car?

texts:
I have a car with a GPS (WGS84) and a camera (who shot the road) fixed on it. I would like to know is there is a way to get the transformation from a OSM node to the GPS Cartesian coordinate fixed on the car ? So that I can get easily project the coordinate of the GPS into the picture (from the camera) coordinate.

I don't know if this is a problem of map matching because I don't want to match a trajectory (like for example https://mapmatching.3scale.net/) but I just want to get the node next to the car and then have the transformation between the two.

Do you have an idea ?


texts: Removing slivers after splitting a polygon with a line

texts:
Image 1 shows before editing.
Image 2 shows 'edited line which i used too split with and with labels
Close up in 3rd image.

I often split polygons with a line (ArcGIS 10.1 Split Polygon tool). In this case, and as an example:  The polygons' boundary is coincident with a line (river); and opposite sides of the line/river. I re - digitise (to improve accuracy) the river and smooth the line.  I then use this line to split the adjacent polygons. Often MANY small polygons (100+) are created and these are extremely cumbersome to manually integrate with the correct 'big polygon' they should be merged with. They are so small that it is by labeling them, they can be seen easily.

I usually seect them mannually and merge them with the correct polygon.


texts: Streaming Raster layer does not export with map at large paper size

texts:
I have a layer in an MXD which is a streaming raster layer from a web tile service, it is a mosaic.  If I output the map as a PDF at Arch C size or smaller the raster shows up in the PDF fine, if I choose Arch D or larger the raster does not show up in the exported PDF.  I have set my visibility scales on the raster layer just in case and it is within the range of visibility.  The layer shows fine on my screen in layout and data view of ArcMap.  It is a streaming layer of orthophotos from Planet.  RUnning 10.3.1 but also reproduced in 10.6

Anyone else have this behaviour and know a fix?


texts: How to calculate mean aspect in circular statistics?

texts:
Measurements of angles (e.g. compass bearings, aspect) are circular in nature (e.g. 0 and 360 degrees coincide) and therefore cannot be meaningfully described by ordinary arithmetic means. I have a shapefile of 22 locations containing aspect value for each point. I want to calculate mean aspect of the locations using circular statistics in R. "circstats" is one of the packages that are suggested to me.
Question: Do I need to transform aspect data before feeding into R (aspect data is calculated in GRASS)? How I can calculate the mean direction of aspect in R?


texts: Find polygons that contain a specific number of points

texts:
I want to make a polygon selection based on a specific number of points.
For example, I want to select those polygons that contain more than 30 points. I don't want to select every polygon that contains a point.

I tried to do something with ArcGis's data reviewer and I checked spatial join but I didn't find anything.


texts: OpenLayers measure tool - instant measurement

texts: I need to use measure tool to measure distances (instant measurement!) in my application. Instant measurement example can be found here. But the main problem is that when I use instant measurement tool and try to measure distance, it sums all the pieces and shows sum. So the question: how can I instantly show only the current piece distance (I mean not to sum up all pieces distances)?


texts: Why do hillshades from Image analysis look much better than the hilshade tool in ArcGIS?

texts:
I have been hoping to make my workflow a little more automated, but the hillshades that result from the hillshade tool just don't do it since they have a static look to them. I've tried to play around with the settings but nothing seems to work. I have added pictures from both methods, the only setting changed from default is setting a z-factor of 4.
This is the result from the image analysis method

These are the parameters from the image analysis

This is the result from the hillshade tool

These are the parameters from the hillshade tool. I did not change any of the environment settings.


texts: Drawing markers with custom SQL in javascript with carto

texts:
I am trying to draw markers with custom psql using the js cartobd.createLayer method. The below snippet doesn’t throw any errors, but doesn’t draw any markers on the map either.

window.onload = function() {
  var map = new L.Map('map', {
    center: [40.7128,-74.0059],
    zoom: 14
  });
  var layer = L.mapbox.styleLayer('mapbox://styles/fluffywolf24/cimhxr71f0021qxno7nts1d4w', { detectRetina: true });
  layer.addTo(map);
  console.log('map:', map);
  cartodb.createLayer(map, {
    user_name: 'whayler1',
    type: 'cartodb',
    sublayers: [{
      sql: "SELECT * FROM ny_listings",
      cartocss: '#map {marker-fill: #F0F0F0;}'
    }]
  })
  .addTo(map)
  .done(function(layer) {
    console.log('done', layer);
    // create and add a new sublayer
    layer.createSubLayer({
      sql: "SELECT * FROM ny_listings",
      cartocss: '#map {marker-fill: #000000;}'
    });
  });
}


texts: When copying information into another layers attribute table objects turn invisible

texts:
When I try to copy the information of one layer into another via attribute tables the data is shown in the table but it's not visible (when I push the button in the attribute table to show only visible objects the new copied ones are not shown anymore). What is my mistake? How can I make the "new" objects visible?

First I built a .shp layer and into this one I want to load all my .gpx data. I I already did this (2 years ago), but now it won't work anymore and I don't know where my mistake is. The only difference I could find is that now by saving the first layer as ESRI shapefile is that I now get 6 instead of 5 files earlier on. Now there is also a .cpg file

I'm not able to explain anymore detailed in ways of informatics.


texts: Custom Script tool (ArcPy) Error: Input Features does not exist or is not supported

texts:
I am trying to create a custom tool that take x,y coordinates, buffers them and then puts them into a new feature class. Whenever I try to run the tool I get:

  ERROR000732 Input Features: Dataset (x, y coordinates) does not exist
   or is not supported

Here is my script:

import arcpy
import os

point = arcpy.GetParameterAsText(0)

outputFC = arcpy.GetParameterAsText(1)
buffDist = arcpy.GetParameterAsText(2)

arcpy.Buffer_analysis(point, outputFC, buffDist)

Here is the error from the tool after the script is run:


texts:
I have an issue with my QGIS 3 installation. For some reason, I can make the app be active over other opened windows (other apps or explorers). So, each time I need open QGIS, I need to hide/close all other windows/apps on top of QGIS. Any idea what causes that and how to fix that?
Here is a video that explains the issue with my QGIS installation - https://drive.google.com/file/d/1CDkO1SIt09v0WWnASVZdGxvNylUec0KK/view (Russian voice)


texts:
How can I correctly concatenate the path of a project with a specific shapefile layer and to be able to load it; this is:

QgsProject.instance().readPath("./")   <---path

layer = iface.addVectorLayer(   path   +  capa x "ogr" ) <--example visual

layer = iface.addVectorLayer(QgsProject.instance().readPath("./") "COB.shp", "COB", "ogr")


texts: Mapping land plots which contain two buildings

texts:
I need to show land plots in QGIS that have two buildings.

My query only shows buildings which are placed on land plots.

Create or Replace VIEW dvi_budivli_view AS
Select
"budivly".*,
"dilyanki"
From
public."budivly",
public."dilyanki"
Where
ST_contains ( public."dilyanki".geom, public."budivly".geom);


texts: arcpy Multiple ring buffer: not working for OUTSIDE_ONLY?

texts:
I am trying to create multiple buffers surrounding input polygon fc. My code works well if I run default variables. However, if I would like to keep only outside buffers, I got back and error:

ERROR 000800: The value is not a member of ALL | NONE.
Failed to execute (MultipleRingBuffer).

Do you have any idea, why my script runs in an error for "outside polygons"? It runs well when run as a tool from ArcToolbox. I am using ArcGIS 10.4, python 2.7.

Here is my code:

import arcpy

arcpy.env.overwriteOutput = True

# Define variables
inFC = "C:/Users/inFC"
outLoc = "C:/Users/unique.gdb"
distances = [100, 500, 1000]
bufferUnit = "meters"

outBuff = outLoc + "/" + "b_" + "Koprova"

# Create buffers
arcpy.MultipleRingBuffer_analysis(inFC, outBuff,
                                  distances, bufferUnit,
                                  "", "OUTSIDE_ONLY")  # this works only with "ALL" instead of "OUTSIDE_ONLY"


texts:
I have some polygons with defined colors according to a column in attribute table. I need an example about how to export these polygons into a DGN file with the color but I can't find any after Google Search. Is this possible in QGIS?

I have python background. So if your answer requires some programming, I can take it.


texts: I am trying to delete points from a shapefile with vector information that intersect a buffer that was created. I have tried using the intersect tool, but that only highlights the information that intersects the buffer, it doesn't highlight them in the attribute table of shapefile that I need to delete them from.


texts:
I'm trying to understand the closest_point function in qgis. For that, I created a polygon layer with two features, and through the Geometry Generator I intend to find the closest point between those two.
This is the expression I'm using, but it is not working:
closest_point(geometry (get_feature_by_id('closest', 1)), geometry(get_feature_by_id ('closest', 2)))

Can someone help me to find out why?


texts: Performing spatio-temporal join/filter between two polygon layers using QGIS

texts:
I have been looking for a long time for a good explanation of a spatio-temporal filter or join in QGIS. In my specific case I have two shapefiles (polygons) - one describes the ignitions and the other represents fire perimeter. Both layers also include an attribute called year in which the two  occurred.
My goal is to filter out the fire polygons that do not match the ignition both spatially but also temporally.

I have tried this by converting my ignitions to points, but that did not solve the issue.
How can I easily perform a spatial and temporal filter/join of the two polygon layers in QGIS?
My plan is to look at how the timing of ignitions may change over a time period, therefore, I want to look at how the timing of ignitions is for each year.
I can also look at the average size of the fires that the ignitions have started for each year.


texts: I need to convert GPS coordinates to EPSG:54004 on Android. Are there any libraries or formulas to accomplish the task?


texts:
I'm using GeoXml3 to parse XML files in Google Map and would like to add OSM as custom maptype. I fiddled with this example but it hides completely the original Google layers.

I've found a couple of other examples but unfortunately all of them are using the V2 API.


texts: How to call a package in QGIS python console?

texts:
I'm learning python and I'm also learning GIS (QGIS to be exact).   What I'm trying to do is use the python console in QGIS to practice and get comfortable with python.

Tutorials that I have found are pretty straightforward.  Moving forward, I'm looking at creating a package and calling that package in the QGIS python console.  Specifically, I have created the following on:
E:\OneDrive\Documents\PythonQGis\PythonScripts

File: crsLayers.py

# just a simple function to test accessibility.
def getDefaultProjectCrs():
    return iface.mapCanvas().mapRenderer().destinationCrs().authid()

I have also incuded __init__.py file in the PythonScripts folder.
it has a single line:

__all__ = ["crsLayers"]

In my QGIS python console I issue the following commands:

import sys
sys.path.append("E:\OneDrive\Documents\PythonQGis\PythonScripts")

I have tried the following:

import PythonScripts
import PythonScripts.crsLayers
crsLayers.getDefaultProjectCrs()

They all return error messages:

  Traceback (most recent call last):   File "", line 1, in
     File
  "C:/PROGRA~1/QGISWI~1/apps/qgis-ltr/./python\qgis\utils.py", line 478,
  in _import
      mod = _builtin_import(name, globals, locals, fromlist, level) ImportError: No module named PythonScripts.crsLayers

or

  Traceback (most recent call last):   File "", line 1, in
   NameError: name 'crsLayers' is not defined

What am I missing?


texts: Is there a Batch Add Description option or script?

texts:
I like to make Layer Packages for our group as this keeps things consistent.  One holdup in making my packages is the requirement for all layers to have a layer description in the General Layer Properties section.
Is there a batch process where I could put a blanket "description" in just to speed things up?
Better yet, is there a script I could use to make "Description" the same as "Layer Name" as a batch process.
I have 297 .asc raster that I would ideally like to have the descriptions match the layer name.  However even just adding a generic description like "asc raster" would at least let me make my layer package without having to click on each individual layer and copy and pasting the layer name into description.


texts:
I'm developing a map for a project that has multiple location but they're all in the same county/ city.

How do I limit the user from panning outside the county/ city in Bing Maps?
I know that I can use disablePanning in map options but how do I set it the City boundaries only?


texts: working with hydrography raster the cost array ignores the end of water

texts:
We have our maps laid out using hydrography maps from NOAA.  We had to run the graphics through several steps to create the water feature raster.  I have a process that pulls in that raster, sets a start and end point, and computes a Polyline which provides my distance between points.  But we were getting strange results that had to be figured out.  I added the cost array to the map and it is a regular, circular growth from the target point out.  In a few places it appears to follow rivers but mostly it ignores the existance of land.

It reports values as being between 0 and 17xx.  I'm concerned about the cost value of 0 as I though no raster can have a value of 0.  So is this because I need to convert 0 to NoData?  And if that is the case how do I do the conversion?


texts: Extracting numbers from string field but omit some unnecessary cells using Python Parser of ArcMap Field Calculator?

texts:
I refer to Extracting numbers from string field using Python Parser of ArcMap Field Calculator?

The script above works fine for me, but I have to develop it a little.
What I want is just to omit values without numbers. As you see in the table TEXTSTRING field with ID = 10 got only letters, it's Ok if Python leaves it empty, I just want make it go through all values in table.

target column is test_srednica2

Below you see my attempt  :/

def makestr(TEXTSTRING):

   lis = []
   for i in TEXTSTRING:

      if i.isdigit():
           lis.append(i)
      elif i.isalpha():
           continue

   return ''.join(lis)

Below i put code which has filled field "test_srednica2"
As you see in table it is working until row "13" there is value "woc". There is no number so i think that`s why this function stops.

I want my function leave cells with letters only, and go throught all the others.

If you look on row "11" you see its going throught even if value is "<Null>" its OK for me.

def makestr(TEXTSTRING):

    lis = []
    for i in TEXTSTRING:

    if i.isdigit():
         lis.append(i)

    return ''.join(lis)


texts: Loading OpenStreetMap data to a custom PostgreSQL schema via osm2pgsql

texts:
What is the procedure to load OSM to a different schema than public?

The osm2pgsql documentation mention the parameter --prefix for specifying the tables prefix (planet_osm by default) but does not mention the schema name.
Running osm2pgsql -h -v to list all available commands, it seems there is no parameter for specifying the schema.

How can I load OSM data to a custom schema?


texts:
I'd like to export 127 raster files from grass to .tif in one operation. I already stumbled over this post How to batch export using r.out.gdal command in GRASS but it doesn't work (I don't even know exactly where to enter this -python console or command shell?).

I'm not too familiar with programming, so I hope there's a simple solution.


texts: GeoJSON: Truncate geometry based on circle

texts:
I have a database with a lot of GeoJSON, mostly representing roads and buildings. This GeoJSON is displayed on a map based on a chosen location, i.e. you choose a point on the map and the roads and buildings near to that location are displayed.

At the moment anything which intersects with a circle around the selected point is displayed. However, I'm at the mercy of the underlying dataset in terms of what is actually displayed. For example, if I select a point next to a long road, then that entire road is displayed.

What I would like to do is truncate the geometry displayed based on the circle, so that if there's a long road I only display the portion of the geometry that is contained within the intersection circle.

How can I do this? I'm not even sure where to start.

Ideally I'd like to do this in .NET as my backend is using .NET core and if I can do it server side I won't have to implement it for each client individually.


texts: Converting x, y, zoom for Calculating tile bounds

texts:
i am working on calculating tile bounds lat lon from x,y,z tile coordinate. i have looked gdal2tiles code structure and i have found this good function:

def TileLatLonBounds( tx, ty, zoom ):
"Returns bounds of the given tile in latutude/longitude using WGS84 datum"

bounds = TileBounds(tx, ty, zoom)
minLat, minLon = MetersToLatLon(bounds[0], bounds[1])
maxLat, maxLon = MetersToLatLon(bounds[2], bounds[3])

with this code and its dependies i am getting result of this for tms scheme:

1.TileLatLonBounds(0, 3, 2) # 66.51326044311185, -180.0, 85.0511287798066, -90.0
2.TileLatLonBounds(3, 0, 2) # -85.05112877980659, 89.99999999999999, -66.51326044311185, 180.0

but i know world extent is for upper left -180 90, for lower right is 180 -90. so for first function (minLat, minLon, maxLat, maxLon) maxLat must be 90  instead of 85.0511287798066 and for second function minlat must be -90 instead of -85.05112877980659.

if i want use gdal translate function for rectifying for an image which contains world extent, then web tiling it, and using with openlayers or anything like this i think there will be problem with click into place of its....

and when opening -180 90 180 -90 tif and -180 85.05112877980659 180 -85.05112877980659 tif with global mapper you can see the diffrence of coordinate...

ok. what is your opinion for fixing this problem? or is it not a problem for tiling with gdal2tiles


texts: QGIS Python Discrete ColorBand Classification

texts:
I'm trying to automate a process where we bring in a DEM, extract the contours as a new vector layer, then set the dem to render the band between each contour line a separate colour.

I want to do is the equivalent of setting up the symbology property box like the screenshot below (with differing max, min and classes values calculated from the dem)

Can anyone suggest a good way to do that?


texts: Transforming geometry points coordinates to WGS84 in geodataframe

texts:
I have shapefile that has been brought into my script as a geodataframe. Initially, it shows no crs although the source originally had a datum = NAD83.

I want the coordinates for the geometry (the POINT; not the fields LATITUDE and LONGITUDE) to be in WGS84 and I assume I need to do more than simply set the geodataframe to WGS84 like this : bldg_feat_cent_df.crs = {'init' :'4326'}. It seems the point coordinates need to be transformed to WGS84. I have tried a number of different things using pyproj, shapely and functools but I am not sure how to correctly apply them:

import pyproj
from pyproj import Proj, transform
from functools import partial
from shapely.ops import transform
....
state_feat = state_bldg_centroid
bldg_cents_cols = ['uniqueid', 'LATITUDE', 'LONGITUDE', 'GEOID', 'CensusPop', 'HU_Pop']
bldg_feat_cent_df =
gpd.GeoDataFrame.from_features(records(bldg_cent_cb_intersects_gdb, bldg_cents_cols))
...
//Get a Geoseries of the geometry
geometry=bldg_feat_cent_df.geometry
project = partial(
  pyproj.transform,
  pyproj.Proj(bldg_feat_cent_df.crs), # source coordinate system
  pyproj.Proj(init='epsg:4326'))

This throws an error RuntimeError: b'no arguments in initialization list' and I am not sure what to do from there. Any suggestions on how to accomplish my goal of transforming my points to WGS84?


texts: Creating an annulus sector with mapbasic

texts: Knowing the center of rotation, start and end angles, inner an outer radius, how can this be achieved with MapBasic


texts: Calculating $area in project with multiple layers with different CRS in QGIS

texts:
I have a project with the CRS set to EPSG:31370. In this project there are a couple of layers, with a different CRS than the project CRS namely EPSG:4326 and EPSG:31370. Now I need to calculate the area of the polygons in those layers in square meters.
However, I can't seem to find how I can calculate the area in square meters when using a shapefile with the EPSG:4326.
I've already googled, read the corresponding data in the manual, youtubed and did a search on this StackExchange, and while I'm convinced my problem is an easy one to solve, I can't seem to find the solution.
What steps should I take to be able to calculate the area of shapefiles with CRS's that have degrees as their units?
I am quite new to GIS.


texts: How to covert a GeoTIFF DEM to a standard raster format without losing quality or resolution

texts:
I'm working on a project, a spin-off of my graduate dissertation, to digitize LiDAR data of a portion of Tampa Bay (20x20 km^2 at 1:1 scale) into a Minecraft map. Ideally I want not just accurate heights but land cover, or else the map would be boring and hard to orient in. At first I tried using geo-mapgen, a program for a similar (some may say too similar) game, but it always gave Lua-related errors (which I am completely ignorant of) when I tried to load a map and I've given up unless the programmer can respond to my posts on its message board. Now I've been trying to use WorldPainter for Minecraft and I've had more success: I uploaded in the program a PNG height map converted from a DEM of the land I want to use and got it running in Minecraft but without land cover, so it's just 400 square kilometers of grass, some beaches, and the ocean. In WorldPainter, you can use a transparent layer of a land cover image, and while it can't automatically convert Minecraft land into biomes assigned to numerical values (which geo-mapgen can do), it can assist greatly in helping you manually sculpt biomes which is the program's specialty.

Here's why I'm asking here: I've used QGIS and data from the NOAA Data Access Viewer to convert the Lidar-derived DEM from GeoTIFF to PNG that was able to be read by WorldPainter, but I'm having more trouble converting a land cover image to a PNG or other format (JPEG, GIF, BMP) while preserving the sharp, cell-distinct quality of the land cover GeoTIFF. I can clip it to the extent of my Lidar DEM and lay it over spatially perfectly as a PNG in WorldPainter but it's too blurry after the conversion to help with making biomes. Can you convert TIFFs to a common format without compromising the quality? The DEM wasn't that big so I'm not concerned about file size. I actually used a free online converter for the lander cover TIFF to PNG as I'm not sure which tool to use in QGIS: neither Save as File nor gdal_translate seem to support the file types used by WorldPainter.


texts:
Is there any better way how to find vector nearest to point then this:

select gid,ST_Distance(ST_GeomFromText('POINT($x $y)',4326),cr_2po_4pgr.the_geom)
as d from cr_2po_4pgr order by d limit 1

Sometime i have feeling that this is looking for nearest start/end point and not vector itself.


texts: ARCGIS API GraphicsLayer is not visible on map

texts:
I'm creating an onClick eventhandler to create a point on every click on the map, but the Graphics layer is not visible;
Has anyone else encountered this? Is there something special I'm supposed to do to show it?
Here is some of the code:

        var point = {
            type: "point",
            longitude: event.mapPoint.longitude,
            latitude: event.mapPoint.latitude,
            spatialReference: event.spatialReference,
        };

       var simpleMarkerSymbol = {
         type: "simple-marker",
         size: "50px",
         color: [226, 119, 40], // orange
         outline: {
           color: [255, 255, 255], // white
           width: 1
         }
       };

       var pointGraphic = new Graphic({
         geometry: point,
         symbol: simpleMarkerSymbol
       });

     let pointLayer = this.map.findLayerById("pointLayer") as GraphicsLayer;
     pointLayer.removeAll();
     pointLayer.visible=true;
     pointLayer.add(pointGraphic);
     this.logger.info("nb of graphics : ", pointLayer.graphics.length)

     });
 };


texts: Getting longest distance between multiple points and creating multiple polygons from lines using QGIS?

texts:
I am still a beginner with GIS.

I can't find an easy way to get a maximum distance from multiple points.

I have a layer of up to 100 points (without a unique ID - maybe that is part of the problem? but I could easily add that) and I want to get the length of maximum distance between any two of them.

I tried plenty of line-points plugins, but nothing seems to work, I also tried to make a minimum convex polygon (basically convex hull of all the points) and get a diagonal of that but failed in all tries again.

How do I do this without using a programming language?

I am using QGIS 2.18.5


texts: geoPHP pointOnSurface always returns null

texts:
Using geoPHP with geos installed. pointOnSurface() for a valid polygon always returns null. While other functions such as centroid() work fine.

Example:

geoPHP::geosInstalled();          //true

$geo->geometryType();             //Polygon

$point = $geo->centroid();        //Point {coords: [-73.931519444925, 40.85512387767]}

$geo->covers($point);             //true

$geo->pointOnSurface();           //null

I'm wondering if anyone else has this problem (or I'm missing something)?

Edit: This workaround seems to do the trick:

$geos_geometry = $geo->geos();                    //GEOSGeometry

$geos_POS = $geos_geometry->pointOnSurface();     //GEOSGeometry

$geos_POS->isEmpty();                             //false

$geos_POS->getX();                                //-73.932666589981

$geos_POS->getY();                                //40.85405746088


texts: How to get layer from modify event of insertVertexCondition method?

texts:
For my dev, I want to get the text of the current selected layer from event of modify interaction.
Below is the code snippet I used, but I can't get text of layer.
var modify = new ol.interaction.Modify({
    features,
    pixelTolerance: 1,
    deleteCondition: event => {
      return false;
    },
    insertVertexCondition: event => {
        const coordinates = event.coordinate;

        const layers = event.target.getLayers().getArray();
        layers.forEach(layer => {
          if (
            layer.getProperties().name === 'section-layer' &&
            layer.getSource().getFeatures()[0].getGeometry().getType() === 'MultiPoint'
          ) {
            const layerSource = layer.getSource();

            const layerText = layer.pleaseHelp(); // Please help with this method.

            const calibrationPoint = new ol.Feature({
              geometry: new ol.geom.Point(coordinates),
              text: 'calibration',
              layer: layer.getSource().getFeatures()[0].getProperties().text
            })
            layerSource.addFeatures([calibrationPoint])
          }
        })

        return false;
      }
  });

  map.addInteraction(modify);


texts: GreenGeeks versus WebFaction for Web mapping

texts: Am looking into choosing a hosting service to test out and ultimately deploy a web mapping service. I want to test out Mapserver, Geoserver, Geonode and Ushahidi towards making a decision so looking into a hosting service. Does anyone have any advice/experience with GreenGeeks or WebFaction ?


texts: Zoom to Feature hotkey in ArcMap Table View

texts:
In ArcMap 10 table view, is there a way to set up a hotkey that when a feature is selected, pressing say "z" is equivalent to right-click-->Zoom To?

I would like to be able to highlight a record and just press "z" to zoom to the feature, "p" to pan to the feature and keep the same zoom, etc.  It seems simple enough, but I cannot seem to find a place to set up this function.

Am I missing something simple? I just want to replicate the function shown below without the 2 additional mouse clicks.


texts: how to extract a variable at specific time and depth from netcdf in r

texts:
I have the following netcdf file containing 7 variables and 4 dimensions.

[1] "file data.nc has 4 dimensions:"
[1] "depth   Size: 12"
[1] "time   Size: 37"
[1] "x   Size: 305"
[1] "y   Size: 432"
[1] "------------------------"
[1] "file data.nc has 7 variables:"
[1] "int stereographic[]  Longname:stereographic Missval:NA"
[1] "float longitude[x,y]  Longname:longitude Missval:1e+30"
[1] "float latitude[x,y]  Longname:latitude Missval:1e+30"
[1] "short chla[x,y,depth,time]  Longname:chla Missval:-32767"
[1] "short nitrat[x,y,depth,time]  Longname:nitrat Missval:-32767"
[1] "short phosphat[x,y,depth,time]  Longname:phosphat Missval:-32767"
[1] "short pbiomass[x,y,depth,time]  Longname:pbiomass Missval:-32767"

I would like to extract several months at a specific depth and convert them to a ASCII raster. I am able to use this command myraster<-raster("data.nc", varname="nitrat", lvar=1, band=2, level=20) but I am confused about which depth and time I am actually extracting, it is unclear to me how to use the parameters of the raster function to point to the specific time and depth.
Also I am not able to read a variable using the function get_att and I get this message

test<-open.ncdf("data.nc")
Error in ncatt_get(test, "time") :
  Error, first passed argument must be an object of class ncdf4

I would like to be able to get a list of the different time and depth so I can use the correct parameters in the raster function.I have both "ncdf","ncdf4" packages installed.I am working in Windows with R 3.0


texts: Explanation of reduce.count() function in Google Earth Engine

texts:
Following the indication on Calculating number of pixels in polygon using Google Earthy Engine? I wrote this script:

var gfc =ee.Image('UMD/hansen/global_forest_change_2018_v1_6');
var lossyear = gfc.select(['lossyear']);
//For loop for each years
for (var py=0; py <19; py++){
 var year = py+2000;
 var y = py;
 var table_comp = path for feature collection;
// current and next year
var lossYear= lossyear.eq(y);
var lossYear = lossyear.updateMask(lossYear);
var yp1 =ee.Number(y).add(1);
var lossYear_pl =lossyear.eq(yp1)
var lossYear_pl =lossyear.updateMask(lossYear_pl)
var gfc_count= lossYear.reduceRegions({
   collection: table_comp,
   reducer: ee.Reducer.count(),
   scale: 30});
}

For each geometry of feature collection I would like to have the number of Hansen loss for the reference year and the next.
I'm writing this script to ask information about function of reducer.count() function. I have noticed that the surface that I calculate starting from the pixels extracted from Hansen is often larger than the area of the polygons I previously calculated. Besides the inevitable mistake,I would ask major information about the function reduce.count().

In particular pixels that are on the border of the polygon how are they counted?


texts: Copy feature output is not accessable by next tool in modelbuilder under ArcGis Pro 2.2.2

texts:
Is there a bug in ArcGis Pro 2.2.2 because every tool after the copy feature tool is throwing an error 000229. But when I navigate within the next tool to the output and set it directly as input it is working.

Has anyone experienced the same behaviour or got a solution to this problem?


texts:
How to publish PostGIS Rasters in Geoserver?

  gdalinfo landsat.tif

Driver: GTiff/GeoTIFF
Files: landsat.tif
Size is 1773, 1428
Coordinate System is `'
Origin = (177595.474744186009048,471385.819080783985555)
Pixel Size = (30.000000000000000,-30.000000000000000)
Image Structure Metadata:
 INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (  177595.475,  471385.819)
Lower Left  (  177595.475,  428545.819)
Upper Right (  230785.475,  471385.819)
Lower Right (  230785.475,  428545.819)
Center      (  204190.475,  449965.819)
Band 1 Block=1773x1 Type=Byte, ColorInterp=Red
Band 2 Block=1773x1 Type=Byte, ColorInterp=Green
Band 3 Block=1773x1 Type=Byte, ColorInterp=Blue

  step 1) raster2pgsql

  raster2pgsql -s 4737 -I -M -C landsat.tif -t 100x100 -l 2,4,6 public.land mosaic | psql -U dragonman -d nyc

step 2) MosicPlugIn install

  step 3) creat MOSAIC table

create table {schema}.MOSAIC (NAME varchar(254) not null,
     TileTable varchar(254)not null,
     minX FLOAT8,minY FLOAT8, maxX FLOAT8, maxY FLOAT8,resX FLOAT8, resY FLOAT8, primary key (NAME,TileTable))
   insert into MOSAIC (name, tiletable) values ('oek', 'landmosaic');
   insert into MOSAIC (name, tiletable) values ('oek', 'o_2_landmosaic');
   insert into MOSAIC (name, tiletable) values ('oek', 'o_4_landmosaic');
   insert into MOSAIC (name, tiletable) values ('oek', 'o_6_landmosaic');

  step 4) geoserver create config file (mapping, connenct, config xml file)

  create a directory /data/coverages this needs to contain 3 files: connect.pgraster.xml.inc, mapping.pgraster.xml.inc, oek.pgraster.xml

  step 5) geoserver store create publish

< >Geoserver Error Log

Caused by: java.lang.Exception: Unable to acquire test coverage for format:ImageMosaicJDBC
    at org.geoserver.catalog.CatalogBuilder.buildCoverage(CatalogBuilder.java:970)
    at org.geoserver.catalog.CatalogBuilder.buildCoverage(CatalogBuilder.java:856)
    at org.geoserver.web.data.layer.NewLayerPage.buildLayerInfo(NewLayerPage.java:346)


texts: Looking for help identifying projections and reprojecting

texts:
I'm still struggling with projections in QGIS. The "Vogelschutzgebiete" reprojection worked. But these ones seem to be tougher:

www.projektidee.org/images/extern/Gemeinden.zip

as well as the same URL with Aemter.zip at the end.

Are these the same projections? Does anyone here by any chance know how to reproject these to EPSG:32633 (WGS84 UTM33N)? (shapes without projection information are a pain in the a...)

Referring also to this thread: http://forum.qgis.org/viewtopic.php?f=2&t=10244&sid=52667e4eb2cfdad92b4f88e160c1b577


texts: Is esri-leaflet officially an Esri project on github?

texts:
Is the information found on github for Esri-leaflet, officially part of the ESRI development team for Leaflet?

What I mean is https://github.com/Esri/esri-leaflet officially being developed by Esri developers.

I know this may sound like a pedestrian type of questions since Github is a social network for developers but I am a bit confused, and trying to figure out if Esri-Leaflet has been blessed by Esri, or it is  independent project where someone was able to developed an extention for Esri, but is not really part of the official Esri developer community.


texts: Adding the selected feature on the map in Google Earth Engine

texts:
How can I add the selected county on the map?

var states = ee.FeatureCollection("TIGER/2016/States")
var counties = ee.FeatureCollection("TIGER/2016/Counties")

var statesNames = states.aggregate_array('NAME')

var getCounties = function(state) {
  // Given a state get all counties
  var feat = ee.Feature(states.filterMetadata('NAME', 'equals', state).first())
  var statefp = ee.String(feat.get('STATEFP'))
  var filteredCounties = counties.filterMetadata('STATEFP', 'equals', statefp)
  var filteredCountiesNames = filteredCounties.aggregate_array('NAME')
  return ee.List(filteredCountiesNames)
}

// Empty Dropdowns
var statesDD = ui.Select([], 'Loading..')
var countiesDD = ui.Select([], 'Waiting for a state..')

// Load states
statesNames.evaluate(function(states){
  statesDD.items().reset(states)
  statesDD.setPlaceholder('Select a state')
  statesDD.onChange(function(state){
    // once you select a state (onChange) get all counties and fill the dropdown
    countiesDD.setPlaceholder('Loading...')
    var counties = getCounties(state)
    counties.evaluate(function(countiesNames){
      countiesDD.items().reset(countiesNames)
      countiesDD.setPlaceholder('Select a county')
   Map.addLayer(ee.FeatureCollection(countiesDD))
    })
  })
})

ui.root.insert(0, ui.Panel([statesDD, countiesDD]));

link:
https://code.earthengine.google.com/d4c8625f6723ce1c59098a0e3bf8d3df


texts: Geonode (Geoexplorer) Query Menu Customize

texts:
I am trying to customize query result menu of Geonode. (Layer->Edit Layer->Query)

For example, I want to export query results in a specific format or create a new layer from selected features from right click menu.
I searched some but cannot find required files and it seems time consuming.

How can I do that?


texts: Does anyone know how to solve GDAL problem with batch translating rasters in QGIS?

texts:
I am trying to use the GDAL module in QGIS to convert some NetCDF files into GeoTiff through batch processing. Unfortunately, I keep getting the same error even after several attempts. This is the error's text body:

    2017-01-06T23:57:04 1   Traceback (most recent call last):
          File "C:/PROGRA~1/QGISLA~1/apps/qgis/./python/plugins\GdalTools\tools\widgetBatchBase.py", line 73, in onRun
            self.batchRun()
          File "C:/PROGRA~1/QGISLA~1/apps/qgis/./python/plugins\GdalTools\tools\doTranslate.py", line 306, in batchRun
            self.runItem(self.batchIndex, self.batchTotal)
          File "C:/PROGRA~1/QGISLA~1/apps/qgis/./python/plugins\GdalTools\tools\widgetBatchBase.py", line 107, in runItem
            self.base.refreshArgs(args)
        TypeError: native Qt signal is not callable

Anyone has any idea how to solve this?

N.B: I've tried this earlier in the day on a different laptop and it worked, therfore the files are not a problem.


texts: Polygon of Voronoi cell using scipy.spatial?

texts:
I have points (e.g., lat, lon pairs of cell tower locations) and I need to get the polygon of the Voronoi cells they form.

from scipy.spatial import Voronoi

tower = [[ 24.686 ,  46.7081],
       [ 24.686 ,  46.7081],
       [ 24.686 ,  46.7081]]

c = Voronoi(towers)

Now, I need to get the polygon boundaries in lat,lon coordinates for each cell (and what was the centroid this polygon surrounds). I need this Voronoi to be bounded as well. Meaning that the boundaries don't go to infinity, but rather within a bounding box.


texts: Fix poor quality when using WMTS to import Mapbox to QGIS

texts:
I have a basemap created using Mapbox and have imported it in to QGIS using WMTS. The quality is terribly pixelated on the map and is just as bad when I try export as an image.
Screenshot:

How can I solve this issue?


texts:
I've got an application running www.carstenhogertz.de/c_section.

It has a popup function showing the underlying point attribute. This function shows a strange behaviour on lower zoom levels.

For instance in starting zoom level clicking on the point for Germany the data for Italy is shown -

When I zoom into the map and click on the German circle it correctly show the data for Germany -

What might be the reason?

I use L.TileLayer.BetterWMS.js to do the GetFeatureInfo request.