texts: Create Space Time Cube From Multidimensional Raster Layer tool missing from ArcGIS Pro

texts: The tool is listed as part of the 'Space Time Pattern Mapping' toolbox, and I can see the other tools within the toolbox, but 'Create Space Time Cube From Multidimensional Raster Layer' does not appear and I cannot search for it.


texts: Display a wms layer automatically along with base layer?

texts:
I have a OpenLayers map that displays a few WMS layers on top of an Open Street Map base layer. I'm would like the county boundary WMS layer to display automatically along with the base layer when the map loads.

var map;
function initMap() {
    var osmlayer;
    var point;
    var dlr_county_boundaries_wms;
    var dlr_libraries_wms;
    var info;

    //create new map object
    map = new OpenLayers.Map("mymap");

    //create new OSM base layer
    osmlayer = new OpenLayers.Layer.OSM("Open Street Map");

    //create WMS layers
    dlr_county_boundaries_wms = new OpenLayers.Layer.WMS("DLR County Boundary (WMS)", "http://localhost:8080/geoserver/egm715/wms",
        {
        layers: "egm715:DLR_County_Boundary",
        transparent: true
        }, {
        projection: "EPSG:3857",
        opacity: 0.5,
        visibility: false
        }
    );

    // DLR Libraries Layer
    dlr_libraries_wms = new OpenLayers.Layer.WMS("DLR Libraries (WMS)", "http://localhost:8080/geoserver/egm715/wms",
        {
        layers: "egm715:DLR_Libraries",
        transparent: true
        }, {
        projection: "EPSG:3857",
        opacity: 0.5,
        visibility: false
        }
    );

//add layers to the map
map.addLayers([osmlayer, dlr_county_boundaries_wms, dlr_libraries_wms]);

//add LayerSwitcher control object
var custom_layerswitcher= new OpenLayers.Control.LayerSwitcher();
map.addControl(custom_layerswitcher);

// change Labels of Base- and Overlay-Section
custom_layerswitcher.dataLbl.innerHTML="DLR Layers";
custom_layerswitcher.baseLbl.innerHTML="Basiskarten";
    // define point as a new LonLat object and transfom
    point = new OpenLayers.LonLat(-6.1832, 53.2560);
    point.transform("EPSG:4326", "EPSG:900913");
    //centre the map and set zoom level
    map.setCenter(point, 11.7);

}


texts:
I posted a similar question before, but I am trying to list all of the element types in the geodatabase, such as feature datasets, networks, rasters, domains, non-spatial tables, etc.

For rasters, I got the following code working.

#Get Raster datasets
rasters = arcpy.ListRasters()

for rs in rasters:
    arcpy.AddMessage(rs)

However, I was told that I should probably use Walk function instead of ListFeature function. The following is the the sample code I found, but was unable to find a way to print the result.

What change do I need to make to the following code to make it print results?

workspace = "c:/data"
rasters = []

walk = arcpy.da.Walk(workspace, topdown=True, datatype="RasterDataset")

for dirpath, dirnames, filenames in walk:
    # Disregard any folder named 'back_up' in creating list of rasters
    if "back_up" in dirnames:
        dirnames.remove('back_up')
    for filename in filenames:
        rasters.append(os.path.join(dirpath, filename))


texts: Retrieving vertices of raster features using python and GDAL

texts:
I have a raster map like the one shown in the picture.
I would like to get the coordinates of the four vertices (A, B, C, D) of the white polygon using python.

The value of the pixels of the white rectangle is 1, black pixels are 0.


texts: Handling user password input in ArcPy scripts?

texts:
I'm currently trying to write a script that has a username & password input which then returns some data if the user has appropriate access.  The script works just fine when I hardcode my own credentials in.

The problem I have is that I don't know how to mask password input using the "Add Script" wizard.  It appears there is an option to accept an Encrypted String as a parameter, which appears to be what I want, but I can't find any documentation on how to actually utilize the GPEncryptedString class within Python.  When I try to extract the parameter with arcpy.getParameterAsText(), it just returns a bunch of asterisks as a string value which obviously does me no good.

I even tried just popping up my own username/password dialog, but scripts have to be ran with the "Always in Foreground" option for the arcpy.mapping.MapDocument("CURRENT")  command to work properly, which constantly steals focus from my dialog.

This seems like a really simple problem, but I'm just not finding any documentation anywhere on this issue and I've now spent more time trying to figure out a way to mask password input than I spent on the actual script itself.


texts: Removing first and third zeros of a fieldname in label expression

texts:
The fieldname only has three characters. I dont want to display zeros in the first and third values. example below.

010 = 1

100 = 10

01Y = 1Y

10Y = 10Y

I have figured out how to strip out zeros except for the second value.


texts:
Where can I download OpenLayers mobile API?
Can I do everything with OpenLayers mobile like OpenLayers API?


texts: Finding azimuths/bearings between many consecutive points in ArcGIS

texts:
I have a ton of consecutive points and I need to know the bearing between each of them (a to b, b to c, c to d,etc).  Ultimately, I will convert the numbers (bearings, azimuths?) to directions: N, NE, NNE, NNW,etc.

I have found a python script for use in the field calculator that comes close to doing this but it returns angles:
https://geonet.esri.com/blogs/dan_patterson/2016/09/01/distance-calculations-using-the-field-calculator

Another option would be a script that can convert the angles returned by the above script to bearings/azimuths. For instance, the above script returns -90 rather than 270 for due west.


texts: Can I use the Search Widget on fields related to a featurelayer?

texts:
Is it possible to point the search widget to related records?  For example, I have a featurelayer of spaces, which is related to a table of names based on the spaceID.  I would like to be able to search for spaces based on names of people related to those spaces.  Something like below:

I'm using ArcGIS Javascript API 4.9.

var searchWidget = new Search({
    view: view,
    includeDefaultSources: false,
    allPlaceholder: "First or Last Name",
    sources:[{
        featureLayer: myFeatureLayer,
        searchFields: ["relationships/0/LASTNAME", "relationships/0/FIRSTNAME"],
        exactMatch: false,
        placeholder: "example: John Doe"
    }]
});

This doesn't work, but I was wondering if there is a similar approach to accomplish this task.  The documentation does not specify.


texts:
So, I have reason to believe that this information is not available. Nonetheless, it involves using Android, and preferably the Google Maps API (though other maps may have APIs which supply this data, so I'm open-minded).

Here is my goal:

Generate a random geospatial coordinate confined to a certain radius of lat/lon (should be easy).
Validate that location, using layering of some API, to determine that it does not require going in the middle of a road, in a body of water, or on private (or otherwise restricted) areas in order to be adjacent to it (within the margin of error of GPS).
If 2 returns invalid, repeat until a coordinate is set.

Now, clearly some companies have access to detailed information about locations (such as Google), but it doesn't appear that they make this information available to developers, having it restricted as a company secret. Does anyone have a workaround?


texts:
Here there are two layers with the representation of the streets and the idea is: We want to update the most recent layer (with the red line) with some of the data storaged in the other one (the light blue), but these layers are similar in the Geographic position that can be very close and a little more distant depending on the place in the map and only one attribute matches in cases that a line is the same in the other one(codlog x codlog5).

For instance: In the middle of the image, there is a square in both layers, but how do I do in the map to recognize the blue one is the same as the green one and update the recent layer?

How can I do it? Is there any plugin or tool thats helps with it?


texts: Simple way to collect the Sentinel data of a given ROI?

texts: If you wanted to collect all Sentinel satellite data for a given region of interest (ROI), say, for a given day or time frame - is there any simple way to do it? That means: Without having to download all the full images manually and cropping the ROI subset manually as well afterwards?


texts: Difference between OpenLayers 2 and OpenLayers 3?

texts:
Can anyone tell me the main difference between OpenLayers 2 and OpenLayers 3?

In other words, what are the main features in OpenLayers 3?


texts: Refreshing open attribute table in ArcMap using ArcPy rather than ArcObjects?

texts:
I am running the Calculate Field GP tool from the arcpy in ArcMap. When a row value is updated, it does not show up untill I use the Reload Cache tool (Table Options > Reload Cache). How do I refresh the attribute table automatically without using ArcObjects?

I have tried updating the definition query as @blah238 suggested Using python, how do you "reload cache" an opened Attribute Table, but this doesn't help. Resetting the data source is not an option because the feature class might be in the editing session and I don't want to lose the connection to the source data.


texts: How to remove some features of a feature layer in the mobile cache

texts:
I have a Hand-Held (HH) device in which I have downloaded to the MobileCache a limited amount of Features from a FeatureLayer. Now I need to remove from the cache all features of a certain type. We were able to delete them from the cache but this places them into a table called fs_table0_deleted_features in the MobileCache.db file. We are only doing a one-way synchronize (download), but if we ever did a two-way synchronize then these features would be deleted from the server too, which is not what we want. Also, if we try to re-synchronize the now deleted features from the server they are not removed from the deleted features table and the deleted features are not placed back into the fs_table0_current_features table; this behavior is causing issues with our application.

We also tried to use the MobileCache.RemoveData method, but this only accepts an Envelope (extent) and a list of FeatureSource. What I want to do is only remove a subset of Features that have been cached leaving the rest cached.  Other solutions we have considered have been:

Copy in memory from the MobileCache the features we want to save, Call RemoveData, then insert into the MobileCache the features that we saved.
RemoveData from the MobileCache and then synchronize with the server for all the features we want
Delete the MobileCache and then copy a backup that was previously made to replace the now deleted MobileCache.
Find the Extent of all the features we want to remove and then call RemoveData with those extents.

Of these, I think #1 or #3 may the better solutions, but they seem overly complex and it could be resource intensive(#1), or may have problems with stale data(#3).

We don't like #2 because the load for synchronizing our devices may be heavy (thousands of devices all must synchronize in a very limited time frame)

The last one, #4 could be a potential problem if there are overlapping features in which we only want to remove one and not all.

Here is how we tried the delete, which is not working and isn't a valid solution.

public void PurgeCacheData(string whereClause)
{
    if (!Cache.IsOpen) return;

    var qFilter = new QueryFilter { WhereClause = whereClause };
    var featureSources = Cache.FeatureSources;

    foreach (var fSource in featureSources
             .Where(fSource => fSource.GetFeatureCount(qFilter) > 0))
    {

        using (var featReader = fSource.GetDataReader(qFilter))
        {
            while (featReader.Read())
            {
                featReader.Delete();
                featReader.Update();
            }
        }

        // Here we tried to delete the entries in the Delete table
        // but this doesn't work. Nothing happens.
        using (FeatureDataTable featureDataTable =
               fSource.GetDataTable(null, EditState.Deleted, null))
        {
            foreach (FeatureDataRow row in featureDataTable.Rows)
            {
                row.Delete();
            }
        }

    }
 }


texts: Coding issue - trying to update single FC using SearchCursor & UpdateCursor

texts:
Fairly new to python and I have a issue with writing new values to an existing text field in a feature class/shapefile. I want to update the existing text values by taking the Shape@ geometry of the polygon and making a series of text values seperated by a comma.
I got it to work to create the text line and I can write the lines to a file - but now I want to take that text string and replace the values in a field that are already there using UpdateCursor. I want to make sure that all the polygon vertices are written to one filed, then move onto the next polygon as well - here I am trying to compare the OID's - if they change then move on.
I am having issues accomplishing that feat - I can get the values to display in the python window - but then I get a
"Runtime error
Traceback (most recent call last):
  File "", line 21, in
TypeError: value #1 - unsupported type: list"

Here is my code so far:

import arcpy

infc = "H:\GIST_8125\Oil_Gas\Gail_test\clip_project.shp"

for row in arcpy.da.SearchCursor(infc, ["OID@", "SHAPE@", "Coord_stri"]):
  # iterate through polygon parts
  partnum = 0

  for part in row [1]:
  # list polygon part points
    for pnt in part:
        print (str(pnt.X) + ',' + str(pnt.Y))
  partnum += 1

# update coord_string

with arcpy.da.UpdateCursor(infc, ["OID@", "Coord_stri"]) as uCursor:
#iterate through rows and update values
  for row in uCursor:
    if row [0] == row [0]:
      row [1] = ([pnt.X, pnt.Y])
      uCursor.updateRow(row)

del uCursor, row

The first part of the polygon is correctly being written to the field - now just need to find a way to write the 4 other points to the same line:

import arcpy

infc = "H:\GIST_8125\Oil_Gas\Gail_test\clip_project.shp"

uCursor = arcpy.da.UpdateCursor(infc, ["OID@", "SHAPE@", "Coord_stri"])

for row in uCursor:
  # iterate through polygon parts
  partnum = 0

  for part in row [1]:
  # list polygon part points
    for pnt in part:
        print (str(pnt.X) + ',' + str(pnt.Y))
        row[2] = (str(pnt.X) + ',' + str(pnt.Y))
        uCursor.updateRow(row)

    partnum += 1

del uCursor, row


texts: Different lines that run on top of each other show as different lines

texts:
Imagine I have some bus lines each with a different colour, and sometimes the different bus lines run on the same road, but I don't want one of them to disappear.
Is it possible to make them show as parallel like in the picture?

If I just simply draw the lines in parallel to each other, they will change appearance a lot as a zoom in or out and I would not have a consistent width between them.


texts:
I recently completed a long bike trip and used a GPS to gather data during the ride, including average speed and elevation loss/gain. I imported the data into QGIS. These data were gathered per mile, and I have them in a point shapefile with the points being the mile marker where the data is stored.

I want to convert these points into a line file that displays different values for each of the mile markers. I am trying to do this using the Points To One plugin but am having trouble, it does not recognize anything from the attribute table.

Is this type of conversion possible, and is it even possible to display different values in a single line shapefile?


texts: How to speed up QGIS browser opening directory?

texts:
If the directory containing many raster files, QGIS browser is very slow?
I think the program is scaning all files which take much times. Are there any setting accelerate this process?


texts: Changing layer color with a button click

texts:
I am trying to change the color of my map layer when filter buttons are selected. So, the goal is: the user sees all of the data in one color. Then, when they click on a filter button, only a subset of the data is shown and in a different color.

I am looking for a command to call the previously defined styles when buttons are selected.

I have several map styles set up like this, each with a different polygon-fill (bluestyle, redstyle, etc):

<style type="cartocss/html" id="bluestyle">
        #blue{
          polygon-fill: #11799F;
          polygon-opacity: 0.7;
          line-color: #FFF;
          line-width: 0.5;
          line-opacity: 1;
        }

    </style>

Then, the code below (much of which is based on CartoDB's interactivity tutorial). Note I added "**" around the spot where I need to plug in this style change.

  function createSelector(layer) {
            var condition = "";
            var $options = $(".layer_selector").find("li");
            $options.click(function(e) {
                var $li = $(e.target);
                var selected = $li.attr('data');
                var type = $li.data('type');
                if (type === "cartocss") {
                    $options.removeClass('cartocss_selected');
                    if (selected !== "bluestyle") {
                        $li.addClass('cartocss_selected');
                    }
                    condition = $('#'+selected).text();
                    layer.setCartoCSS(condition);
                }  else {
                    $options.removeClass('sql_selected');
                    if (selected !== "") {
                        $li.addClass('sql_selected');
                    }
                    if (selected.indexOf('Button1') !== -1) {
                        **(some code that will change the style to redstyle)** ;
                    } else {

                    }
                    layer.setSQL("SELECT * FROM " + tableName + selected);
                }
            });
        }

I am new to Javascript.


texts: Using a projection other than Mercator with proj4js 2.5 and OpenLayers 5.3

texts:
I am following the OpenLayers FAQ regarding the use of coordinate systems that do not rely on the infamous Mercator projection. I thus produced a very small HTML with a simple example that could be later fed to an OpenLayers View object:

<!DOCTYPE html>

    <title>CRS example</title>
    <link rel="stylesheet" href="https://openlayers.org/en/v5.3.0/css/ol.css" type="text/css">

When this file is loaded in the browser the console shows an error message: "Uncaught EPSG:152160". A stack trace is shown:

This error message is shown both in Chromium and Firefox.

I tried changing the identifier string passed to the defs method, but the result is the same. Without the register method being applied the OpenLayers View object ignores the CRS passed in the projection parameter and uses Mercator instead.

How can I use the desired CRS with this software setup?


texts: Calculating delaunay triangulation with two diagonal lines using ArcGIS Desktop?

texts:
How do I obtain a line output similar to a delaunay triangulation but with  two diagonal lines, from a point layer containing a regular-spaced point grid using ArcGIS Desktop?

The current result is this:

As you can see, this solution only creates a single diagonal line. I've also tried Thiessen polygons, but it merely creates a square grid encompassing each point. I also tried 'Tesselate' after creating a square layers from my points, but it provides an identical solution to the delaunay polygons (a single diagonal line). The aim is to later use this output as a dummy road network (the single diagonal line produces sub-optimal routes currently).


texts: I have a table where points are stored with recorded time. I am storing points with 1s interval. Now if the user stayed approximately at the same spot for 2 hours, a number of useless points are dumped in db. How to avoid this or how to convert delete all points except 2 points. 1, when user entered that location and 2, when he left that location ?


texts: Mapping variation within multiple categories using QGIS?

texts:
I am learning qgis and running into a problem.  Here is what I want to do.

I have county data.  It falls into 4 types.  Each type includes a discrete set of counties, that does not overlap with the other types.  Within each set, counties are assigned continuous values.

I have each type as a separate field in an attribute file.  That is, values for type 1 are in a variable called "type 1" and includes values only for counties within that group.  Remaining counties are missing values.  Then type 2 counties have are in a separate field, with values for those counties, and everything remaining is missing.  ditto for type 3 and 4.

How can I include all four attributes in one map?

Maybe a different way to say it is:  I have four categories, and I want to show variability within each category, but maintain their categorical distinction.  So counties in group 1 could be red, with lower to higher values going from lighter to darker red.  counties in group 2 could be blue, with lower to higher values going from lighter to darker blue.  And so on.

Does that make sense?  if so, how do I do it in QGIS?


texts: Does the TIGER geocoder intersect the polygons to match an address?

texts:
This post has a very good overwiew of how open source geocoders works.

1) Address normalization
2) String matching.
3) Interpolation.

It shows a simple case of how interpolation in made, taking as example a street segment that is supposed to contain houses between 300-400, just
interpolating the position over the line.

Now my question here if the zip code is specified in the address and the street segment covers more than one zip code? is this consider during the interpolation? and if is, How it's done?


texts: QtDesigner/QGIS 3: Which widget to use for multiple selections in ValueRelation?

texts:
I created a field of type ValueRelation allowing multiple selections in PyQGIS. Everything works as it should in the attribute table (see here for example). Now I'm trying to create a custom .ui in QtDesigner. And I can't figure out which widget to use for my field to work properly with the .ui. If I use a listView/listWidget, nothing shows in the list. If I use a comboBox, my list of values shows in the comboBox, but I can't select more than one value (and if I select multiple values directly in the attribute table, then none of them is selected in the comboBox). What is the appropriate QtDesigner Widget to use for multiple selections ValueRelations ?

In case this helps, here's the code I use to populate the list (and that works and poulates the list with the comboBox, but not with the listView/listWidget) :

    layer_ref_milieux=QgsProject.instance().mapLayersByName('REF_grds_types_milieux')[0].id()
    type_grds_types = 'ValueRelation'
    config_grds_types  = {'AllowMulti': True,
    'AllowNull': True,
    'FilterExpression': '',
    'Key': 'grd_type_short',
    'Layer': layer_ref_milieux,
    'LayerName': 'REF_grds_types_milieux',
    'LayerSource': 'dbname=' + bdd + ' host=' + host + ' port=' + port + ' user=' + user + ' key=\'id\' checkPrimaryKeyUnicity=\'1\' table="ref"."grds_types_milieux" (geom) sql=',
    'NofColumns': 1,
    'OrderByValue': True,
    'UseCompleter': False,
    'Value': 'grd_type'}
    grds_types_widget_setup = QgsEditorWidgetSetup(type_grds_types,config_grds_types)
    layer.setEditorWidgetSetup(grds_types_index, grds_types_widget_setup)


texts:
I need to correctly project a TIFF file with specific X/Y coordinates. I am using gdal_translate.

GCP File

-of GTiff -gcp 0 0 -124.5806892 49.9259877 -gcp 2362 0 -65.7337720 49.9259877 -gcp 0 1032 -65.7337720 24.5371839 -gcp 2362 1379 -124.5806892 24.5371839 /home/weather/tmax_plot.tif /home/weather/tmax_plot_proj.tif

Python

gdal_translate --optfile /home/weather/gcp

gdalwarp -r bilinear -t_srs EPSG:4326 /home/weather/tmax_proj.tif /home/weather/tmax_proj_crs.tif

gdalwarp -t_srs EPSG:4326 /home/weather/tmax_proj_crs.tif /home/weather/tmax_proj_final.tif

Expected Result

Actual Result

I think that I have my Uppers/Lower GCP options mixed up but, can't tell for sure.


texts: Exporting layer from Google My Maps as KML has no geometries

texts:
I have a CSV with an address field (no lat/long), so I went to Google My Maps, imported the csv and geocoded of the records.

Now I export the layer as KML and added the KML to QGIS and Google Earth Pro, but in non of them I can see the markers. In fact in QGIS, the KML layer is added as datatable (no geometry).

What's happening? How can I export a geometry layer from google My Maps?


texts: I have Himawari files in NetCDF format (.nc). But they do not appear to have the conversion factor for transforming albedo to radiance. In MATLAB, I tried ncdisp, but it didn't have the factor. Any advice?


texts: How to import Ordnance Survery raster layer to ArcGIS Desktop?

texts:
I have downloaded raster map data from ONS, the data has come in the form of two files... one is full of  .tlw files and the other is full of .tab files...

How do I import it into a raster layer to ArcCatalog?


texts: Deleting values from raster layer using R

texts:
I am looking for a way to delete values in my raster layer in R programming. Assume I have a raster layer with the following unique values:

0,1,2,3,4. minimum is 0 and maximum is 4.

How do I delete all values of 0 from my raster layer to remain with the following unique values?

1,2,3,4


texts: How to convert a column in an inventory table into a map with graduated colors for polygons

texts:
I am new to GIS, so I am sure I will use some of the wrong terms. I have joined a CSV file to the attribute table of a vector shapefile with polygons (basically all watersheds on earth), using an ID common to both watershed polygons (original vector shapefile) and to the new data linked to watersheds (reliance on blue water, for 92 crops).

So I have 92 columns in the attribute table (plus column for the ID), one for each crop, and in each column I have a single unique data value for each of the watersheds. I would then like to represent the reliance on blue water, by each crop and watershed, with a graduated color ramp. In the end, I would like to have 92 separate shapefile maps, in my ArcMap document.

Somebody gave the advice to click on properties of the shapefile, then symbology, then quantities, then in the value field drop-down list choose the crop - but when I do that, the drop-down list only offers the watershed ID number (despite the data having been joined in the attribute table). I cannot choose categories: unique values, even though the crops show in the drop-down list, because there are hundreds of unique values, one for each watershed.

Any input?


texts: Joining fields with long field attribute names

texts:
I am working with the "Joins" option in QGIS 3.10, which joins the fields from two attribute tables based on a common shared attribute column.

However, this does not work for attribute values with a large length.

E.g. the following field value can be joined:
2404252385

While this one cannot:
551482643342511892418838370

I checked carefully if it may be due to other problems, but it clearly seems to be the problem of the number of characters. For example, I can filter both features by the attribute fields, which returns the respective row. However, the joined fields nevertheless show 'NULL', i.e. the common field name could not be found.

Can somebody confirm this problem and is there maybe a workaround/solution to fix this?

Edit: I am referring to field values, not field names. I tested around and possible issues that might be the problem are:

The columns that I want to join are both called "fid"
The columns that I want to join are of type string, not integer


texts:
I use a PostGIS database with temperature data. I want to use a spatial interpolation method but not the supported IDW method. I need to use kriging, but it is not part of the supported interpolation methods of PostGIS.

I found kriging interpolation method written in R. Now I want to know how can I combine the R code with the data? The problem is that I've never worked with R before.

Edit:

To clarify: I build 4 climate station with Arduino Megas. They send their climate data to PostGIS over wifi/gsm. The received data should now be spatially interpolated and saved to a new table, so I can use in GeoServer for visualisation. Only the interpolated data is shown on a map. So this all should be done on server side.


texts: Create continuous raster from point values in 10.1?

texts:
I have a point dataset of nearly 400 points containing percentage values the point is representing

Is there a way to create a smooth continuous raster surface across the whole of the UK of the point values?
I have tried using Feature to raster and point to raster and I get similar results of the points being surrounded by the raster square cell in a certain colour depending upon the value.
What I am looking for is raster covering the whole of the uk so most of it will be zero where there is no data or no points and the same colour and then then where the point values are they go up in colour showing the difference in values.

I have used kernel density before but I don't particularly want the raster create using magnitude of points near each other.


texts:
German original:

ich beschäfte mich gerade mit QGIS.
Ich habe über die Funktion mein Shapefile hinein geladen. Jetzt ist mein Problem:
In dem Shapefile sind unterschiedliche Punkte, wie Punkte für einen Pfeil oder nur für eine Linie,... und diese sind alle in dem einen Shape.
Ist es irgendwie möglich, dass man den Shape irgendwie auftrennt. Also das es statt einem dann angenommen drei gibt oder so?
Ich hoffe, ihr wisst was ich meine. Könnt ihr mir weiter helfen?

Translation:

I am working currently with QGIS. I have loaded the shapefile into my function. Now my problem is: in the shapefile are different with different kinds of points (one for arrow, one for ...) and these are all in one shapefile.

I would like to know how can I create a separate shape file for each kind of point. So there are three instead of one then accepted?

Can you help me please?


texts: QGIS Treat null as zero in raster calculator

texts:
In QGIS - I have an irregularly shaped raster - I wish to clip an irregular area out of it (call it raster A) with that area set to null.

I masked the area (vector polygon) I want to remove and used warp > reproject to create the second raster (lets call it B) that effectively contains the portion of raster A I want to set to null.

I intended to use raster calculator to subtract raster B from raster A which would leave the areas to be removed with value 0, and I could then reproject this setting 0 to null -  however this doesnt work as outside the masked area 'A - null = null' and I end up with a calculated raster of 0 in only the area I want to remove, and null everywhere else.

How can I treat the areas on raster B that are null as zero in the calculation?  Alternatively, is there any less convoluted way I can complete this process?


texts: Is WGS84 a Coordinate system or projection system?

texts:
WGS84 is often referred to as a "projection system" and sometimes it is confusing to me about what a "projection system" is and where WGS84 falls in one.

Could someone please clear that up for me?


texts: In PostGIS, how do I merge a multilinestring into a single linestring, keeping the pieces in proper order?

texts:
I have a multilinestring made up of several pieces, some of which are facing the wrong direction.
(Accessible here: http://pastebin.com/89Zcixv5 )

I think what I'm looking for is some utilization of ST_MakeLine, but when I try to use it, the pieces of the trail get put together wrong. (See: http://i5.minus.com/iJ5QVbaA9rm2g.png )

I have a ton of trails in this format that I need to convert from MultiLineStrings to LineStrings containing the whole trail.

If anybody has any ideas, I'd greatly appreciate all the help I can get!

Thanks!


texts: transformations from a sphere using R for raster and polygon data

texts:
I am planning on using the MODIS burned area products which were produced by Maryland University. These come as either polygons or rasters. In the user guide it says

A shapefile version of the MCD45 product is derived from the monthly Geotiff version by the University of Maryland. The shapefiles are available with the same projection (Plate-Carrée) and geographic extent as the Geotiff sub-continental windows

I used the following code in R to convert the projection from sphere-based datum to WGS 1984 UTM ZONE 38N (into the same projection as the rest of my data, including the dem data).

Is the code below correct for converting these projections?

library(raster)
projection(fire) <- projection(dem)

How would I do this for polygon data?


texts: Creating a new shapefile from points which are overlapped by polygons PyQGIS3

texts:
I want to extract points which are overlapped by polygons to a shapefile with PyQGIS. I typed the code below, I don't get an error, but I get points that I have on the beginning. So my intersection doesn't work, and I can't see why.

points = QgsVectorLayer(points_path, "layer", "ogr")
polygons = QgsVectorLayer(polygon_path, "layerpolygon", "ogr")
int_points = QgsVectorLayer("Point", "inttersected_points", "memory")
pr = int_points.dataProvider()

for polygon in polygons.getFeatures():
    geom_polygon = polygon.geometry()
    features_point = points.getFeatures(QgsFeatureRequest().setFilterRect(geom_polygon.boundingBox()))
    for one_feature_point in features_point:

        if one_feature_point.geometry().intersected(geom_polygon):
             int_points = one_feature_point.getFeatures()

QgsVectorFileWriter.writeAsVectorFormat(int_points, save_int_as_shp, "CP120", points.crs(), "ESRI Shapefile")


texts: Splitting large number of polygons at same time into smaller polygons?

texts:
I would like to split the large polygons in the picture below into smaller polygons using the border. I have about 4000 lots that I would like to divide and find the area of each of the smaller sections of the polygons but every process I have used so far will not allow me to work on more then one polygon at a time.

I am a student that does not have much experience with Python but can try it out.


texts: unique value symbolize layers from a lookup table

texts:
I need to create a large number of Arcmap layer files (v9.3 or 10), perhaps a hundred, from lookup tables. I don't have to set the symbol patch itself, that can be random, but I do need the value and label attributes filled out. The lookup tables look like:

"FC","VALUE_FIELD","VALUE","LABEL"
"Wooded_area","CODE","1240012","Extracted"
"Wooded_area","CODE","1240022","Interpreted"
"Wooded_area","CODE","1240032","FCS-EOSD"
"Wooded_area","CODE","1240042","Land Cover, Circa 2000 - Vector"

And the result would resemble (except the symbols themselves):

I'd prefer a python based solution, via arcobjects+python if required, C# next, but anything that works, even partially or only conceptually, is welcome.


texts: QGIS plugin: set the position of the featureForm

texts:
Im developing QGIS plugin (Python3)

to open the featureForm (depending on some of my conditions (matching geometries,...))I use this code:

            layer = self.iface.activeLayer()
            features = layer.getFeatures()

            for feat in features:
                geom_poly = feat.geometry()
                if self.feature.geometry().intersects(geom_poly):
                    layer.selectByIds([feat.id()])
                    f = feat
            self.iface.openFeatureForm(layer, f, True)

the code works correctly, what I want is to set the position of the featureForm

how can I "catch" featureForm and use something like QDialog.move(x, y)


texts:
I've created a query that creates a line that is perpendicular to my original line using ST_OffsetCurve to draw a line left and right of my original geometry, which I then select the line center of using ST_Interpolate and connect via ST_MakeLine. If I use a small distance (100m offset on a 1km line) it connects as expected and no topology errors are returned. However, if I increase the distance it becomes very prone to topology errors, which are thrown by GEOS. They are always "1st arg isnt a line" or "depth mismatch" errors. I should note that the lines are very 'blocky', as in they have many 90 degree turns, since it's satellite derived data.

I'm not entirely sure what causes it, but ST_OffsetCurve doesn't duplicate my line perfectly, there are always some artefacts to the left and the right, especially near the end of the line, which I believe to be the source of my errors.

I'd explain the code but I'm afraid it's not really easy to explain what's going on. There are five spatial functions jammed into one to form the line (transform, linemerge, line interpolate point, offset curve, remove repeated points). It's below if you care to take a look at it.

SELECT
  ST_MakeLine(
    ST_Line_interpolate_point(
      ST_Transform(
        ST_OffsetCurve(
          ST_RemoveRepeatedPoints(
            ST_LineMerge(
              ST_Transform( cl.geom,3857 )
            )
          ),
          {length},
          'quad_segs=0 join=bevel'
        ), 4326
      ), 0.5
    ),
    (ST_Line_interpolate_point(
      ST_Transform(
        ST_OffsetCurve(
          ST_RemoveRepeatedPoints(
            ST_LineMerge(
              ST_Transform( cl.geom,3857 )
            )
          ),
          -{length},
          'quad_segs=0 join=bevel'
        ), 4326
      ), 0.5
    )
  )
)

Which alternatives to ST_OffsetCurve do I have here, or is there perhaps something I'm doing wrong?

Alternatively, is there any way to perfectly offset the line without any deviation from the input line?


texts: Geodetic Coordinates And Latitude, Longitude

texts: Would it right to say that geodetic coordinates (phi & lambda) are the same as latitude & longitude?


texts: Join CSV to shapefile failing when in saved map document in QGIS

texts:
I am trying to join a CSV file to a shapefile by matching on attributes.

Here is the minimal workflow:

Import shapefile in blank document.
Save attributes of shapefile to CSV.

2.5 eventually, I want to do some calculations in R here and add some new columns to this table, but at this point, I don't even open the CSV outside of QGIS.
import CSV file to map.
join, using add vector join

With this procedure, the join works perfectly if I open a blank map document, and DO NOT save it.  For reference, the files I'm working with only have a few hundred features, they are not large.  When the join works, it's nearly instant on my macbook pro.

However, if I save the map document as something, any time before my step 4 above, then the join simply hangs.  I don't get an error, but QGIS simply stops responding, and I need to force quit.

I'm using QGIS version 2.18.7, on a MacBook Pro, running OS 10.12.5.  I have been trying to solve this problem for a while, so I've updated everything.  For a while, I thought the problem might be related to incompatibilities in the carriage return format, either in my import/export/reimport procedure, or in how QGIS re-imports the file itself, based on some of the comments here:
QGIS 2.0 text delimited layer importing data as single row

My open-a-new-blank-document workaround is liveable for now, but why is this happening?


texts: Cannot read very large data from GeoServer

texts:
I've an application which using Leaflet reads data from GeoServer.

That GeoServer reads data from one underlying PostgreSQL database, that database has a field in one of its tables where I'm storing images via its base64 representation as the PostgreSQL datatype text.

Due to the representation the string to read can be very big, for example, right now I'm attempting to read one which is more than 100000 characters.

By being that big it looks like Leaflet cannot read it, and instead always uses one previosuly cached value. I've checked that if the string is shorter then it has no problem displaying the new image once it reads the base64 data from the database.

Maybe I could use some type of compression for the string and later decompress it, but I'm not sure this is always going to work, so I'd prefer to leave it as it is.

What could I do so it's possible to read such huge strings?

Due to request I'm posting code:

This is the code I use to read:

camping_vectorgrid.on('click', function(e) {

//addInteraction()
displayBase64Image(e.layer.properties.naamnl,e,map);
   })
  .addTo(map);

As long as the string in the naamnl of the database is not too large, e.layer.properties.naamnl gives the proper value (at least in zooms that are near the polygon). But if the string is too long it always uses the cached value.

Images get shown in this popup:

function showPopup(imagePlaceholder,e)
{
  L.popup()
        .bindPopup({maxWidth : 600})
       .setContent(imagePlaceholder)
       .setLatLng(e.latlng)
       .openOn(map);
}

imagePlaceholder contains the image that is generated by setting as its source its base64 representation.


texts: How do I keep the coordinates while calculating tiles from a raster with Rasterio?

texts:
I'm trying to split a georeferenced TIFF image into several tiles for further processing.

I have managed to calculate a single tile and give it the right CRS, evething fine so far.

Once I load the tile up in QGIS, it shows me that the top left corner of the tile is at 0/0. It should be arround 483544,5579377.

I suppose that somehow I am missing to give the right transformation. The code looks like this:

with rasterio.open(imgPath) as src:
    tile = src.read(1, window=Window(1000, 1000, 500, 300))

test = rasterio.open(imgPath)
print(test.crs)# 25832

with rasterio.open(resPath, 'w',driver='GTiff', width=500, height=300, count=1,
        dtype=tile.dtype, crs=src.crs) as dst:
    dst.write(tile, window=Window(0, 0, 500, 300), indexes=1)
    print(dst.crs)#25832

Before anyone asks, I have build this very same tool using GDAL and calling it via subprocess, but it is a nightmare to bundle up with pyinstaller, so I want to try it this way and hope for the best.


texts: Leaflet popup design not working properly on IE11

texts:
I am using Leaflet on my map and I styled popup using Bootstrap 4. Design works good for all browsers expect IE11. Can someone help me to find the problem and fix it?

How popup looks in Firefox:

How popup looks in IE11

Link on site - http://www.rpbv.lv/uploads/Karte_v2/


texts: Inconsistencies in OSM data in AWS Athena

texts:
I have a process that exports OSM data from AWS Athena (https://aws.amazon.com/blogs/big-data/querying-openstreetmap-with-amazon-athena/) to Parquet files on s3 for processing with Spark. While debugging the data set created on 05/05/2020, I found several problems and inconsistencies.

A turn restriction (type='relation', id=10902143) that I know doesn't exist on the ground. It exists in the data I copied from Athena to Parquet, but does not exist in either the planet or planet_history tables.
SELECT * FROM planet
WHERE
type='relation'
AND id=10902143

and
 SELECT * FROM planet_history
 WHERE
 type='relation'
 AND id=10902143

yield zero results. However going to https://www.openstreetmap.org/relation/10902143 shows that it
existed and was cancelled. Why does it not appear in the history table?

Way 91275335 appears in the planet and planet_history tables but not in the Parquet files we created, so it appears that at the moment we ran the query it did not exist. I also noticed that the latest version showing up in the planet_history table for this way id is 17 while the version appearing in planet is 20. To reproduce, run:
 SELECT * FROM planet
 WHERE type='way'
 AND id=91275335

and
 SELECT * FROM planet_history
 WHERE type='way' AND id=91275335
 ORDER BY version DESC

The queries I use to create the parquet are quite straightforward so I doubt they're responsible, but I'll provide them just in case:
CREATE TABLE "temp_db"."osm_ways_temp_to_delete_999" WITH
                (external_location = 's3://my-bucket/osm_ways',
                 format = 'Parquet',
                 parquet_compression = 'SNAPPY')
                 AS
                SELECT id, type, tags, nds, timestamp, version FROM "temp_db"."planet"
                WHERE type='way'
                AND tags['highway'] IN
                ('motorway', 'trunk', 'primary', 'secondary',
                'tertiary', 'unclassified', 'residential',
                'service', 'motorway_link', 'trunk_link',
                'primary_link', 'secondary_link', 'tertiary_link',
                'living_street', 'road', 'construction', 'raceway')

And
CREATE TABLE "temp_db"."osm_relations_temp_to_delete_999" WITH
                    (external_location = 's3://my-bucket/osm_relations',
                     format = 'Parquet',
                     parquet_compression = 'SNAPPY')
                     AS
                    SELECT planet.id, planet.type, planet.tags, planet.members, planet.timestamp AS time_stamp, planet.version, blah.ref AS restriction_ref_id, blah.type AS restriction_ref_type, blah.role AS restriction_ref_role,
                    cast(blah.ref as varchar) || '-' || blah.type || '-' || blah.role AS restriction_string
                    FROM "temp_db"."planet"
                    CROSS JOIN UNNEST(members) as t(blah)
                    WHERE type='relation'
                    AND
                    tags['type']='restriction'

(Please note that anyone with AWS can create the planet and planet_history tables in Athena using the very simple instructions in the link above, but keep in mind that Athena queries are billed based on the amount of data scanned.)
My working theory is that these situations are caused by accidental edits that are eventually cancelled but are somehow released into the datasets provided to Amazon AWS.
If anyone could please shed some light on how to understand, detect and/or avoid such situations, I will be very grateful.


texts:
So, I have this little script (which I would like to build more upon).

I was able to make it work in QGIS3.2 (just change all paths and environments to 3.2) but now, after deleting QGIS 3.2 and installing 3.4 or 3.6 it does not work anymore.

I am able to start QgsApplication. But after testing, isValidLayer(), the result is false. So any processing creates empty output.

Environments and python is set up only for QGIS3.6 (Python37). If I install QGIS 3.2 and set it up back, all works just fine. But i am not able to change it to version 3.6. It is driving me crazy!

Any ideas?

import os, sys, datetime, glob, re, pprint
from qgis.core import *

from qgis.analysis import QgsNativeAlgorithms

sys.path.append(r'c:/SW/QGIS3.6/apps/qgis/python/plugins')

#-------------------------------------------------------------
#  start QGIS processing application

QgsApplication.setPrefixPath(r'C:/SW/QGIS3.6/apps/qgis', True)
print(QgsApplication.prefixPath())
qgs = QgsApplication([], False)
qgs.initQgis()

from processing.core.Processing import Processing
Processing.initialize()
QgsApplication.processingRegistry().addProvider(QgsNativeAlgorithms())
#-------------------------------------------------------------

# folder = r'd:\01_Projects'
# files = glob.glob(folder + '/**/BUDOVY_P.shp', recursive=True)

out = r''
file = r'd:\01_Projects\BUDOVY_P.shp'

layer = QgsVectorLayer(file, 'tmp', 'ogr')
print(layer.isValid())
qgs.exitQgis()

Results is:

C:/SW/QGIS3.6/apps/qgis
False

pth file in site-packages for setting up PyCharm:

C:\SW\QGIS3.6\apps\qgis\python
C:\SW\QGIS3.6\apps\Python37\
C:\SW\QGIS3.6\apps\Python37\Lib\
C:\SW\QGIS3.6\apps\Python37\Lib\site-packages
C:\SW\QGIS3.6\bin
C:\SW\QGIS3.6\include
C:\SW\QGIS3.6\apps\qgis\bin
C:\SW\QGIS3.6\apps\qgis\python\plugins
C:\SW\QGIS3.6\apps\grass\grass-7.6.0
C:\SW\QGIS3.6\apps\Qt5\bin

env variables:

C:\SW\QGIS3.6\apps\qgis\bin
C:\SW\QGIS3.6\bin
C:\SW\QGIS3.6\apps\Python37\Scripts

The same part run in QGIS returns:

Python Console
Use iface to access QGIS API interface or Type help(iface) for more info
Security warning: typing commands from an untrusted source can lead to data loss and/or leak
file = r'd:/01_Projects/BUDOVY_P.shp'
layer = QgsVectorLayer(file, 'tmp_lyr', 'ogr')
print(layer.isValid())
True


texts: Storing brief cartographical tips, information, possible future designs?

texts:
I find inspiring maps and/or brief, 1-line design tips (i.e. "when satellite image is used as a basemap, I like to desaturate the tarnation out of it" that I may want to consider later using in a project or as inspiration as a design for a future map, from a variety of sources (blogs, tweets, forums, here) during my daily internet browsing and online interaction.

I'm curious to hear how each of you, especially cartographers and others specializing in map design store all of this information ?

Right now, I store these bits of information across the sources that I find it on: favorited tweets, bookmarks in my browser (do you use tags ?), pinboard, favorited gis.se questions, and a text file (with the tidbit of information and the URL of where I found it).

I occasionally notice that I ended up bookmarking it in a different place, or not knowing that I had that information already saved somewhere until after I needed it.

Where do you find it efficient to store these tips ? Do you only look at these saved tips when you're looking for something specific or is it beneficial to occasionally review them without a specific goal in mind ?


texts: Cannot enable "Track Edits" on ArcGIS Online Feature Layer

texts:
I have a feature layer that is published to ArcGIS Online and being edited by users there.  I am not able to enable "Track Edits" for this feature layer.  Under the feature layer properties, when I click on "edit", check the "Track Edits" option and click save, the feature layer still shows Track Edits as being disabled.

I have previously been able to enable track edits for feature layers successfully but now it just won't enable.  I've tried this on a couple feature layers and none of them allow me to enable track edits.  Anyone out there experience this issue/know of a solution?

This screenshot shows where I enable "Track Edits" with a check mark.

This screenshot is what I see after clicking save on the previous image.  Track Edits is disabled.


texts: Draw ellipse based on values within an attribute table

texts:
I would like to be able to draw and ellipse on a map based on values within an attributes table. At the moment I have 6 transects at a site with values based on IQI (Infaunal Quality Index) values; these values are used to determine a footprint on the seabed.

I have an attribute table with 6 transects with points along those transects with a know IQI value based invertebrates identified from a sample.
From these values we have generated an IQI value on an S curve (EXCEL) that is demened an exceptable area of effect. From this S curve I have added a line within the attributes table for this value and its distance, where I would like the ellipse boarder to fall.

I was hoping to do it based on an interpolated raster layer, then generate a contour layer and use those contours to generate a polygon based on that contour value, but this just gives individual contours around each transect, rather than what I would like, which is a contour around the whole site.

My other solution was create another layer with manually drawn line features along the transect based on the S curve value then physically draw an ellipse based on these lines. This doesn't seem like the most accurate solution.

Does anyone out there know of a tool that will do this?


texts: Setting up definition query to run through ModelBuilder on ArcMap?

texts:
I am trying to make a model that will symbolize data quickly, but wanted to see if a model could perform a definition query or multiple definition queries on a large data set.

If this is possible, how could one set that up?


texts: Install Orfeo Plugin in QGIS - Plugin is not already in QGIS core

texts:
I am brand new to QGIS (coming from ArcGIS), and am needing to install the Orfeo toolbox to do some image classification. I found a nice tutorial here with an overview on how to install it: https://www.orfeo-toolbox.org/CookBook/QGISInterface.html#open-processing-settings . However, I'm still having trouble with it.
The tutorial specifically states that if you are running a version of QGIS greater than 3.8 (I am running 3.14), the Orfeo Toolbox (OTB) core should already be under the "installed" tab on the plugins manager. However, it is not, and I do not know where to go from here!
I downloaded the OTB from https://www.orfeo-toolbox.org/download/ , but cannot seem to find the OTB core. What can I do to install this toolbox? Like I said, I'm brand new to this program, so I'm having some trouble accomplishing this.


texts: Seeking Mobile GIS applications for Android Tablets?

texts: I know that ArcGIS is available for Android, but does anyone know of any other GIS apps that are available for Android tablets?


texts: Assign more than one button to bindpopup

texts:
I have defined two buttons and want them both to be appeared in my bindpopup. Each of them seperately works ok. But when I try to add both of the buttons in the bindpopup, it shows only one of them on the map and inside popup. I do not know how can I solve this problem.

    let btnDelete = document.createElement('button');
    btnDelete.className = 'ti-trash';
    btnDelete.onclick = async () => {
      await api.deleteZone(zone)
      layer.remove()
      }

    let btnEdit = document.createElement('button');
    btnEdit.className = 'ti-pencil';
    btnEdit.onclick = function () {
      layer.editing.enable();
    }

    layer.bindPopup( btnDelete, btnEdit ).openPopup();


texts: error while trying to compile quantum gis: unresolved external symbol class QgsSearchTreeNode

texts:
i'm trying to compile the quantum gis sources code following the step by step how to compile guide, but i'm facing an error when code is linked.

I'm compiling in Windows 7 x64 in Microsoft Visual Studio 2010 Professional and i'm using the environment variables created by OSGeo4W.bat file.

The exact message error is this one:

1>     Creating library C:/qgis-1.7.4/build/src/core/RelWithDebInfo/qgis_core.lib and object C:/qgis-1.7.4/build/src/core/RelWithDebInfo/qgis_core.exp
1>qgssearchstring.obj : error LNK2019: unresolved external symbol "class QgsSearchTreeNode * __cdecl parseSearchString(class QString const &,class QString &)" (?parseSearchString@@YAPAVQgsSearchTreeNode@@ABVQString@@AAV2@@Z) referenced in function "public: bool __thiscall QgsSearchString::setString(class QString)" (?setString@QgsSearchString@@QAE_NVQString@@@Z)
1>flex_qgssearchstringlexer.obj : error LNK2001: unresolved external symbol "union YYSTYPE yylval" (?yylval@@3TYYSTYPE@@A)
1>C:\qgis-1.7.4\build\src\core\RelWithDebInfo\qgis_core.dll : fatal error LNK1120: 2 unresolved externals

When i compile the qgssearchstringparser.cpp.rule with diagnostic i have the following messages! The problem is from bison compilation? It as bison conflicts shift reduce. How can i resolve them?

1>  C:\Progra~2\GnuWin32\bin\bison.exe -oC:/qgis-1.7.4/build/src/core/qgssearchstringparser.cpp -d -v -t C:/qgis-1.7.4/src/core/qgssearchstringparser.yy
1>  if errorlevel 1 goto :VCEnd (TaskId:11)
1>  Generating qgssearchstringparser.cpp (TaskId:11)
1>  C:\Progra~2\GnuWin32/share/bison/yacc.c:158: m4: Warning: Excess arguments to built-in _m4_popdef' ignored (TaskId:11)
1>  C:\Progra~2\GnuWin32/share/bison/yacc.c:1071: m4: Warning: Excess arguments to built-in_m4_popdef' ignored (TaskId:11)
1>  C:\Progra~2\GnuWin32/share/bison/yacc.c:1740: m4: Warning: Excess arguments to built-in `_m4_popdef' ignored (TaskId:11)
1>  C:\Progra~2\GnuWin32/share/bison/yacc.c:1741: error: m4_divert_pop(0): diversion mismatch:  (TaskId:11)
1>  C:\Progra~2\GnuWin32/share/bison/yacc.c:155: m4_divert_push: 0 (TaskId:11)
1>  stdin:1: m4_divert_push: KILL (TaskId:11)
1>  C:\Progra~2\GnuWin32/share/bison/yacc.c:1741: the top level (TaskId:11)
1>  C:/qgis-1.7.4/src/core/qgssearchstringparser.yy: conflitos: 15 de deslocamento/reduþÒo (TaskId:11)


texts: Dissolve features according to an order set by a field

texts:
I want to dissolve overlapping feature according to a specified order using an attribute field. For instance if I have three overlapping features, like the image below, I'd want to keep the purple feature, then the red one and then what's left the remaining feature

Dissolving would look like this.

I would need to do this with thousands of features so editing manually is not an option.


texts:
I am looking for a simple way to concatenate one geospatial pdf (with multiple layers in a single image) and regular pdf. I have tried tools such as pdfunite, pdftk and pdfconcat, without success. The result is always a regular pdf, losing the layers and geospatial references.

Does it exist a simple tool for merging/concatenate geospatial pdfs? Have anyone else solved this problem using another method?


texts: What does xmin, xmax, ymin, ymax mean on Grid Layer?

texts: I made a grid layer over a country and the data display on it is ID, xmax, xmin, ymax and ymin. It's clear that the last 4 parameters are the dimensions of the square of each ID, but what is the reference? I was expecting a latitude/longitude coord and I dont know how to work with these ones.


texts:
I am building an app that has some GIS requirements that go a little beyond (I think) what I can retrieve from existing online web services resources so I am attempting to build a server of my own to handle these requests.

The app is going to make requests to the GIS server specifying lat, lon and I want the server to return a JSON (or XML) document back containing information such as:
Habitat : Urban|SubUrban|Rural|Forest|Old Growth Forest|Desert|Alpine

Desert|Lake|River|Sea|...
City/Town:
State/Province : Initially just for the US.
Country:
Alias/Also Known As: i.e. "The Big Apple" : "New York City"

I also want the server to respond to such queries as:

"Can the Atlantic Ocean be seen from present location?"
"Can Mount Adams, Washington, USA be seen from present location?"

these will be formatted and POSTed in a well defined JSON document

I would also like the solution's stack to be cloud friendly, so I can utilize Heroku or Joyent or Amazon's cloud services and get any mapping datasets from publicly available sources. I do have the need to add specific location data points to these maps as well and have been able to prototype this out using a mongo db running a simple node.js server script (more out of curiosity than utility, python may be my server code in the end) this serving location information to an iOS client.

So far, I've just been able to insert & index documents with longitude & latitude and query with a given rectangle. I'm interested in performing the more complex tasks listed above. I am not beholden to mongodb or any other technology is something else makes more sense or is more compatible.


texts: QGIS Clip to Hemisphere plugin giving 'NoneType' error

texts:
I'd like to use the Clip to Hemisphere plugin in QGIS to reduce artifacts generated by my Orthographic Projection.

I'm using a custom CRS:

+proj=ortho +lat_0=70.6632 +lon_0=23.6817 +x_0=0 +y_0=0 +a=6371000 +b=6371000 +units=m +no_defs

with a Natural Earth admin shapefile

When I load the plugin, I enter lat and lon, click run, select my custom CRS and get the following error:

'NoneType' object has no attribute 'difference' See log for more detail.

The layer generated clips to random countries instead of clipping the hemisphere like I would like


texts: Changing symbol scale in ArcMap legends?

texts:
I'm using ArcGIS 10.3 and the "Legend Wizard" does not offer me a method to include symbol size within my legend. I have my symbols by color in the legend but I would like them by symbol as well, as seen in the example jpg here.

I saw a similar question from a few years ago, but the solution will not work with 10.3.


texts: how to convert lat and longitude into distance near poles?

texts:
I am trying to convert latitude and longitude into meters in polar regions. I am looking for higher accuracy upto 0.00 meters preferably. I tried the haversine formula but it doesnt seem to be accurate enough for small increments in latitude.
does anyone know a better formula that can be used near poles?


texts:
I have a feature layer found in this map: http://arcg.is/0L8SS5  which is to be used with Collector as a part of a classroom activity.  The feature layer is producing the error message "The shape must have a Z value".

Is the only solution to this error to recreate the feature layer without a z value and republish or can settings be adjusted in AGOL to prevent this error?


texts: After joining a CSV to a shapefile in QGIS, some place names are left with null values. How do I get the original values into the shapefile?

texts:
I'm editing a shapefile for an election map, so I downloaded some county level vote totals and am adding them to a preexisting shapefile. This works for the most part (I can produce most of a new map with the added data) and I have been able to solve most problems I get, but I am stuck on something.
After joining the edited CSV file, most of the new attribute table fills out correctly and can be mapped, but some precinct vote counts are left null. After investigating, I discovered these null precincts are all precinct names that are reused in other counties (counties which have not yet been completed). I have only completed 5 counties so far. In the attribute table of the CSV file, the same precinct fields are not null, and correctly filled out.
The data appears to be organized using two fields, with one field labeled "COUNTYFP", indicating the county using an FIPS code, and the other labeled "NAME", indicating the name of the precinct within the county. Some precinct names are used in multiple counties, and I'm guessing that is the problem, though the original shapefile I started with didn't have any problems with that until I joined more data.

Here is a sample of the attribute table of the CSV file. This is just one small county (Alexander county, IL), but it illustrates the point. Notice how the precinct named "CACHE" has data across all fields.

After joining, this is the resultant portion of attribute table for the same county. The same precinct ("CACHE") now has null values across all fields.

This is the map I get, notice that there is a missing portion. my code is definitely working, bit I think the problem is with the join I am performing. There is another "CACHE" precinct in another county. It has exactly the same name, but a different "COUNTYFPS" number.

Here is snapshot of the the two precincts in an attribute table for the shapefile of a different election (a shapefile which I downloaded). It differentiates the two, despite them having the same name. The resulting map colors both precincts differently in accordance with my color coding.


texts:
My map legend looks like this:

I'd like to indent the circled items farther to the right but can't find any settings to adjust this. Using the item properties I know how to offset the label portion, but not the patch.

This is pretty much the same question I asked as Indenting subitems in QGIS Composer legend?, but now trying to do this with ArcMap instead of QGIS.


texts: How to simplify polygon shapefile and maintain topology?

texts:
I'm looking for a straightforward answer that actually works...not a 'duplicate question' flag. None of the other solutions posted in similar threads actually work...they all have dead-end answers and/or links. I've searched for hours and spent many more actually trying the methods I've found and none seem to work.

As for my question; I'm looking to generalize the polygons of a shapefile to reduce their resolution. Topology needs to be maintained for future geometry functions such as difference and union. There doesn't seem to be any way to preserve topology in QGIS...is there any other way? Any suggestions will help (as long as it doesn't involve ArcGIS).

Previous questions which don't have useful answers:

Simplifying Shapefile using ArcGIS Desktop or alternative? - topology is not considered

Generalizing polygon file while maintaining topology in QGIS? - no accepted answer

How to check Topology in QGIS?

Also another questions of mine:

Simplify by number of points in QGIS


texts: Creating new print composer programmatically?

texts:
I would like to create a print composer (e.g. A3 landscape) containing a single map with a few lines of code. I do not need legends, labels, scale bars and so on.

Does anyone know how to do this?


texts: --flat-nodes osm2pgsql data import and TileMill

texts:
This is more so a question regarding the import process with osm2pgsql and the use of the --flat-nodes option as  have seen that this is very useful especially when it comes to importing the entire planet into a PostgreSQL/PostGIS database. But what exactly is this file used for.

All I really understand about this is that it stores Nodes data in a separate file which is much smaller than the data that would imported into the database, but how is it used and are there any effects to how OSM data is handled in the database if this file is deleted?


texts: Efficiently convert from Shapefile to GeoJSON without GDAL

texts:
Ogr2ogr is the standard way to convert files from Shapefile to GeoJSON on the command line. But it has a few limitations for me:

The output format is not condensed (ie, it contains a lot of needless whitespace, which adds up for huge files)
There doesn't seem to be an option to suppress null values (again, very inefficient on space)
GDAL is an enormous package that does a lot of stuff I don't need (so it's a very heavy dependency for this one tiny job).

Is there any better tool for the job?


texts: Adding points from JSON to ArcGIS for Android?

texts:
I am new to Android, the ArcGIS Platform and ArcGIS for Android.

I need to load points (offline) from a JSON file (GeoJSON) onto the map which displays  a local Tile Package (*.tpk) and shapefile.


texts: QGIS: Derived coordinates not matching delimited text input coordinates

texts:
Have a geo-referenced map and a delimited text file with coordinates of points to plot over top of the map.

Issue is that the points are displaying far down and left of the map. Clicking on them shows the lat and long are correct, but the derived coordinates are not.

Attempted to change the CRS for both the map and points, though did not find a combination that worked. Currently leave both on the same CRS of EPSG:3400 which seems to be reasonable.

Also approached the issue by changing settings during importing the delimited text layer, did not find anything wrong with it.

Adding more points with different coordinates resulted in the same issue, the coordinates being in close vicinity, far down and left of the map.

The map is accurate, and while hovering over the city of Calgary the coordinates displayed on the screen are reasonable and where the points should be.

Hoping to hear what could be causing this issue of the points mapping to the wrong locations, and a solution that they appear where expected.

Using QGIS 3.4


texts: Deleting predefined high values from multi band raster in QGIS?

texts:
I'm working with MOD-16-evapotranspiration and I have a problem. I made a composition 180 bands, but in this images I have  predefined high values(such as 32767, 32765..). So, I would like to hide this values for all 180 bands.
I'm working with QGIS software.


texts: Examples or documentation of the atlasvariables

texts: Are there any examples or a documentation of the atlasexpressions??


texts: Is it a way to use a postGIS/pgRouting enabled database with a PaaS like Windows Azure? If not, do you know alternatives?


texts: How to select lines which are 60% contained in polygon?

texts:
I have two layers : a line layer (red) and a polygon layer (green).

I would like to select the lines which 60% are contained in the polygon.

How can I do that?

Until now, I used "select by location" with the "have their centroid in" method, but in some cases, this method is no longer appropriated.


texts: Convert GPS coordinates to Web Mercator EPSG:3857 using python/pyproj

texts:
I've created a couple of functions using the math package from python 3.5 to transform GPS coordinates to x/y web mercator (pseudo mercator) as explained here: http://wiki.openstreetmap.org/wiki/Mercator#Java

I took the Java version and ported it to Python:

import math

# derived from the Java version explained here: http://wiki.openstreetmap.org/wiki/Mercator
RADIUS = 6378137.0 # in meters on the equator

def lat2y(a):
  return math.log(math.tan(math.pi / 4 + math.radians(a) / 2)) * RADIUS

def lon2x(a):
  return math.radians(a) * RADIUS

Then I can use them to get for example the projection of Trafalgar Square (taking the GPS coordinates from the Google Maps URL):

ts_gm = [51.50809,-0.1285907]
print('latitude web mercator y: {} longitude web mercator x: {}'.format(lat2y(ts_gm[0]), lon2x(ts_gm[1])))

Which gives back:

latitude web mercator y: 6711665.883938471 longitude web mercator x: -14314.651244750603

However I would like to achieve this using pyproj syntax, but I am lost in the details of the string I should pass to build a projection using Web Mercator EPSG:3857.
I suspect the formulas above have some sort of precision issue and I would like to rely on a library like pyproj hoping to minimise projection errors.

Is there a minimum working example I could see to achieve the mapping in python using pyproj?


texts: Removing/hiding items in TOC panel (pie chart symbology)

texts:
I work in ArcMap 10.4. I need to find a way to remove/hide some of the items displayed in the TOC after applying Pie chart symbology to features.
The main idea is the items are somewhat duplicated but nonetheless needed for the correct display of data.
Later the layers will be published for browsing in a map app on Portal for ArcGIS environment. There I need only one set of items to be displayed (as seen in the picture).
Just to clarify - this has nothing to do with legends.

Is there a way to do it?

So basically I need to hide or remove items labelled as '_2'


texts: Dimensional analysis/conversion problem: Map scale to image resolution (dpi)

texts: Given a map scale of 1:17,200 and a raster whose cell size is 6 inches, what resolution (dots per inch) should I choose when exporting the map to an image in order to maintain the raster's native resolution (1 cell = 1 dot). Need a quick refresher on the conversion formulas.


texts: Classifying counties to color bottom / top 25% in area using ArcMap?

texts:
I have a map of my state with county boundaries. What settings should I use to color the counties based on their area but only if its area falls in the bottom or top 25% of county area?

For example, if I have 100 counties, sort them by area then color the top 25 (25%) and bottom 25 (25%).

I found the Symbology tab (Layer properties) and tried various combinations under Quantities, all using the SHAPE.AREA value. I tried changing the Classification. None of  the combinations seem correct or require much manual manipulation. For example, I created 10 classes but set the fill to None for all classes except the bottom and top two classes.

Which classification method should I use and what is the recommended approach to performing this task?


texts: How do I solve Gdalwarp error 'too many points failed to transform' for remapping Geostationary to Lambert conformal?

texts:
I am trying to remap from Geostationary to Lambert conformal using gdalwarp.  My input data is in netcdf, and is in geographic coordinates (degrees) and I would like to output the remapped data to netcdf.  I have created a corresponding vrt file for the input netcdf data.  Gdalwarp will output the netcdf file, but the output data is all zeros and I receive the following error:

Creating output file that is 5120P x 5120L.
Processing input file netcdf.vrt.
ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.
Warning 1: Unable to compute source region for output window 0,0,5120,5120, skipping.
0...10...20...30...40...50...60...70...80...90...100 - done.

I attempted the following command:

/usr/bin/gdalwarp -s_srs "+proj=geos +h=35785831 +lon_0=-75 +x_0=-0.151844 +y_0=0.151844 +a=6378140 +b=6356754.99999591 +units=degrees +no_defs" -t_srs "+proj=lcc +ellps=clrk66 +a=6378137 +b=6378137 +e=0.0818191910435 +lat_0=24.9999 +lon_0=-95 +lat_1=24.9999 +lat_ts=25.0001 +units=meters +no_defs" -te -1952976.3246 -828316.5944 3248431.6754 4373091.4056 -of netCDF -geoloc -overwrite -r bilinear -ts 5120 5120 netcdf.vrt out.nc

Can gdalwarp remap from geographic coordinates to projected?  Or do I need to translate geographic to projected first?  Also, can gdalwarp read projection information straight from netcdf or do you NEED to write to .vrt first?

Here is what gdalinfo outputs from the input file:  (it is a GOES 13 file from CLASS)

Size is 512, 512
Coordinate System is `'
Metadata:
  NC_GLOBAL#Conventions=CF-1.4
  NC_GLOBAL#Satellite Sensor=G-13 IMG
  NC_GLOBAL#Source=McIDAS Area File
Subdatasets:
  SUBDATASET_1_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":auditTrail
  SUBDATASET_1_DESC=[3x80] auditTrail (8-bit character)
  SUBDATASET_2_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":data
  SUBDATASET_2_DESC=[1x665x2036] data (32-bit floating-point)
  SUBDATASET_3_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":lat
  SUBDATASET_3_DESC=[665x2036] lat (32-bit floating-point)
  SUBDATASET_4_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":lon
  SUBDATASET_4_DESC=[665x2036] lon (32-bit floating-point)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)

And additional gdal info on the data variable:

Driver: netCDF/Network Common Data Format
Files: goes13.2013.100.174518.BAND_04.nc
Size is 2036, 665
Coordinate System is `'
Metadata:
  data#coordinates=lon lat
  data#long_name=0-255 Brightness Temperature
  data#type=VISR
  NC_GLOBAL#Conventions=CF-1.4
  NC_GLOBAL#Satellite Sensor=G-13 IMG
  NC_GLOBAL#Source=McIDAS Area File
  NETCDF_DIM_EXTRA={time}
  NETCDF_DIM_time_DEF={1,4}
  NETCDF_DIM_time_VALUES=1365615900
  time#long_name=seconds since 1970-1-1 0:0:0
  time#units=seconds since 1970-1-1 0:0:0
Geolocation:
  LINE_OFFSET=0
  LINE_STEP=1
  PIXEL_OFFSET=0
  PIXEL_STEP=1
  SRS=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"]]
  X_BAND=1
  X_DATASET=NETCDF:"goes13.2013.100.174518.BAND_04.nc":lon
  Y_BAND=1
  Y_DATASET=NETCDF:"goes13.2013.100.174518.BAND_04.nc":lat
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  665.0)
Upper Right ( 2036.0,    0.0)
Lower Right ( 2036.0,  665.0)
Center      ( 1018.0,  332.5)
Band 1 Block=2036x1 Type=Float32, ColorInterp=Undefined
  NoData Value=9.96920996838686905e+36
  Metadata:
    coordinates=lon lat
    long_name=0-255 Brightness Temperature
    NETCDF_DIM_time=1365615900
    NETCDF_VARNAME=data
    type=VISR

Any help will be greatly appreciated!


texts: Route between two point layers following streets

texts:
I have a problem to be solved with shortest path route.

What I have:

layer 1 with startpoints
layer 2 with endpoints
layer 3 with streets (lines)

These two point layers are originally the start and endpoints of a layer of lines, but now i'm looking for an alternate route for these lines, using the shortest path algorithm.

I want to calculate the shortest path between these two point layers. The two pointlayers have both an attribute where the startpoint and endpoint have the same unique value, so they can be connected on that attribute value I suppose.

I already did a try with pgRouting: i followed a tutorial and already did pgr_createTopology on the streetlayer, so the vertices of these streets are already made.

But now i'm stuck, can for example pgr_dijkstra be automised for a whole layer? (because each pointlayer contains +10.000 points)

My idea was to use nearest neighbour to link the start and endpoints to a vertice of the streets. And then use pgr_dijkstra on that streetlayer so that the shortest path is calculated when the values of the pointlayer are the same, so startpoint(UNIQUE_VALUE) = endpoint(UNIQUE_VALUE).

How is this possible, or does someone have another solution?


texts: Finding straight streets going from south to north using QGIS

texts:
How do I find straight streets going from south to north with an inclination of 5° maximum using QGIS.
Here is an example

The red line indicates the street I am looking for.
These streets are very common in the United States, but rare in Europe.


texts:
Is ArcGIS designed in a way that it is possible to open and work in multiple ArcMap instances/documents?

I am wondering, when I have some time-consuming operations running in the ArcMap document 1, but want to work in the meantime on another project where I need to test my arcpy code, which I would like to do in ArcMap document 2, if the operations in ArcMap document 2 influence the performance/duration of the operations in ArcMap document 1?


texts:
I am aware that problems like mine exist (example) but I think, this is slightly different.

I want to use MultiValue to let the user choose various inputs at the same time. The problem is, that it's only possible to use MultiValue in combination with String and not, how I would need it, as SQL.

Because I'am trying to create a toolbox, these are my properties (and on the right side my MultiValue-Field in the Toolbox:

I did it like this:

multival = arcpy.GetParameterAsText(1)
vals = string.split(multival, ";")
where_clause = string.join(vals, "+")

My problem is, that my output now has quotation marks and is not combined to one big SQL-expression.

The expression I am aiming for would look like this:

"Name = 'Store_1' OR Name = 'Store_2' OR Name = 'Store_5' OR Name = 'Store_8'"

But right now it looks like this:

'Name = 'Store_2''+'OR Name = 'Store_3''+'OR Name = 'Store_4''+'OR Name = 'Store_5''

Is there a way to have the outputs without quotations marks and added to one expression?


texts: Creating an online form with intergrated google maps

texts: I am looking at creating an online form that includes a map (google or something similar) where the user can draw a poly and when the user hits submit, i get all the info plus a shapefile or kmz  of the area drawn, have been looking for various solutions to this and not getting very far.


texts: QGis2.8 ValueMap: Why is the first row of my ComboBox empty?

texts:
I define a field of a vector layer as ValueMap. I define the values manually, then the feature attribute dialog appears with the combobox. But the first row in my combobox is empty.
I use QGis2.8. In QGis2.0 the first row was not empty.

Why is the first row of my ComboBox empty?
If I click "ok" immediately without using the combobox I get a strange character () in the attribut tablel.

Regards


texts: OpenLayers custom request for WMS GetFeatureInfo

texts:
In OpenLayers, there is a WMSGetFeatureInfo control that can be attached to a layer and listen to click or hover events on the map, requesting info on the feature being currently clicked (or pointed) at.

However, how to make a custom GetFeatureInfo request, with the following input?

layer
[x,y] mouse position OR [lon,lat] position in the map

Is it possible to do this elegantly with the WMSGetFeatureInfo control, or do I have to ajax the request myself?


texts: Remove gaps between polygons of a shapefile

texts: I have extracted a polygon shapefile using raster to vector tool from TIFF image, and it has produced polygons from the Raster, but it has also generated the gaps between the polygons as attached in below image. What is best way to remove those gaps between polygons?


texts: Iteratively join linestrings that intersect?

texts:
I'm struggling to iteratively join intersecting linestrings of a geometry collection using st_join, meaning I want to join intersecting linestrings, but want to build on that join (i.e. joining all the linestrings that intersect in total).

Doing a reprex would be difficult for this, but curious if anyone has any suggestions.

The net results would be that one linestring joins to all its intersecting linestrings, and remains a linestring, not multiline, then that new line string joins to all intersecting line strings, and so forth until all line strings are consolidated if they intersect.

I've tried for loops and am just spinning wheels at this point. I realize this is not optimal for reproducibility, but figured I'd take a shot in the dark asking here.


texts: Checking for proper administrative location?

texts:
I have two vector layers. The first layer contains points. This layer has a lot of fields with data like name of district, voivodeship etc.

The second layer (polygon) is showing administrative division (commune, district, voivodeship).

Based on the second layer (polygon) I have to check that names in table (first layer - points) like "commune, district , voivodeship" are correct.

Is there any solution?

I am using ArcGIS 10.2 for Desktop.


texts: How to antialiase tiles when seeding a layer from an GeoTiff in GeoServer

texts:
I have an GeoServer running that uses GeoWebCache.
As a data source I use a 1.5 GB big map for a layer.
The outcome of the tiling and scaling are very rough looking tiles.

Is there a way to antialiase the tiles, that they are more pleasant to look at?


texts: Creating layer representing non-intersecting area of two shapefiles in QGIS?

texts: I have a shapefile for a country and a shapefile representing the mobile network coverage of an operator in that country.  I would like to create a shapefile that represents the parts of the country not covered by the mobile operator.  I tried Symmetrical Difference which worked around the edges but that did not get at the internal gaps in the network coverage map.


texts: How to use WMTS url service in my android map application as basemap?

texts:
I am building a android map application and i want to use as basemap some wmts url services of our country.

How can i implement this and get layers from WMTS url services as TileOverlay to google map android api?
Example:

Wmts source:
http://geoportal.asig.gov.al:85/utm34n/wmts/1.0.0/WMTSCapabilities.xml

the tileoverlay url inside xml configuration:

http://geoportal.asig.gov.al/Services/host/host.ashx?url=http://localhost:85/utm34n/ortofoto/{TileMatrixSet}/{TileMatrix}/{TileCol}/{TileRow}.jpg

but Google Map Android Api needed this configuration of url:

http://my.image.server/images/x/y/z.png

x- x coordinate ( start from 0 to maxwidth)
y- y coordinate ( start from 0 to maxheight)
z- Zoom

refer to documentation: https://developers.google.com/maps/documentation/android-api/tileoverlay

How can I format the url from the wmts service ResourceURL( from WMTSCapabilities.xml), to integrate inside google map android api?


texts:
I want to display terrain in OpenLayers.

So, I tried to use dem(.img) and convert it to dem(.jpg) on QGIS.

But, the effect is very bad.

My expected effect is (This picture is captured from google earth)

Which method is better for display terrain in OpenLayers?

This is my website for test:http://www.cv.nctu.edu.tw/shapefile-js-openlayers/ol_simple.html


texts: I'd like to see where commercial districts are, using color-coding. Is there any free map website, other than Google Maps, which can help me?

texts:
Background

Google Maps nowadays shows "areas of interest" in pale orange. (Source.) As far as I can tell, "areas of interest" are basically just high-density commercial districts.

I think it makes sense to show commercial zones in their own color. It's definitely more interesting to walk through a busy area such as Times Square than through a quiet residential street on Staten Island. In fact, even low-density commercial areas are probably more interesting to walk through than residential areas.

You might say, "well, just walk along arterial roads". But this doesn't always help. In my city, large parts of many arterial roads are zoned for residential use (e.g. mid-rise apartments and condos) rather than for commercial or mixed use.

Google Maps does color-coding of "areas of interest" reasonably well. It's true that pale orange is a little bit hard to see at low zoom levels on a low-resolution cellphone; but this is just a minor quibble. Still, I wonder who else shows these areas; maybe one of Google's competitors does a bit better at showing where these districts can be found.

My question

Is there any map website, other than Google Maps, which can use color-coding to show me where commercial areas are?

The website must be free for personal use, since I don't want to pay money to use it.

Edit

On second thought, the official Toronto zoning map is definitely a fine option. (The middle of another post explains how to read the map.) I think I like the fact that it shows all commercial zones: even low-density ones like the ones on Marlee Avenue. Of course, it only covers Toronto, Canada; but it still meets my needs fine.


texts: How to save symbology for raster and also its associated 'value' (using unique values)?

texts: I have a raster in ArcMap (built using sum function in cell statistics), and only want to show the water parts of it (band 128). In doing so I only need to show multiples of 128, however there are over 40 bands which are multiples of 128. I have done 1 of these rasters (but I have 15 more I need to do), however I am wondering if there is a way to save the band value and also the symbology so that I do not need to spend so much time going through each one and deleting the unnecessary bands (of which there are thousands of).


texts: Spatial join of feature class on layer class with PyQGIS

texts:
I have three layers:

One point layer consisting of residents
One polygon layer with buffers of the first layer (one polygon for every point)
Another point layer consisting of containers

I use QGIS 3.4.1

It is possible that multiple container points are within a buffer polygon. My goal is to select all container points within a certain buffer with PyQGIS. Hence, I want to know which containers are near a certain resident. I use the following code:

from osgeo import ogr

#Open datalayers
inw = r"SomePath"
afc = r"SomePath"
inw_buff = r"SomePath"

#set-up for residents
driver_inw = ogr.GetDriverByName("ESRI Shapefile")
dataSource_inw = driver_inw.Open(inw,0)
layer_inw = dataSource_inw.GetLayer()
layerDef_inw = layer_inw.GetLayerDefn()

#set-up for containers
driver_afc = ogr.GetDriverByName("ESRI Shapefile")
dataSource_afc = driver_afc.Open(afc,0)
layer_afc = dataSource_afc.GetLayer()
layerDef_afc = layer_afc.GetLayerDefn()

#set-up for buffers
driver_buff = ogr.GetDriverByName("ESRI Shapefile")
dataSource_buff = driver_buff.Open(inw_buff,0)
layer_buff = dataSource_buff.GetLayer()
layerDef_buff = layer_buff.GetLayerDefn()

#some coordinates from the residents point layer
for f in layer_inw:
    coordx_inw = f[7]
    coordy_inw = f[8]
    ID_inw = f[9]
    break

#find the corresponding buffer polygon
for y in layer_buff:
    if y[9] == ID_inw:
        temp_buff = y

processing.run("native:selectbylocation", {'INPUT':afc,'PREDICATE':[6],'INTERSECT':temp_buff,'METHOD':0})

This gives me the following error:

Traceback (most recent call last):
  File "C:\PROGRA~1\QGIS3~1.4\apps\Python37\lib\code.py", line 90, in runcode
exec(code, self.locals)
File "<input>", line 1, in <module>
File "<string>", line 77, in <module>
File "C:/PROGRA~1/QGIS3~1.4/apps/qgis/./python/plugins\processing\tools\general.py", line 96, in run
return Processing.runAlgorithm(algOrName, parameters, onFinish, feedback, context)
File "C:/PROGRA~1/QGIS3~1.4/apps/qgis/./python/plugins\processing\core\Processing.py", line 139, in runAlgorithm
raise QgsProcessingException(msg)
_core.QgsProcessingException: Niet mogelijk uit te voeren algoritme
Kon bronlaag voor INTERSECT niet laden: ongeldige waarde

Translated from Dutch to English it says that the value for INTERSECT is invalid. I looked it up and it does indeed not except a feature class. I tried to place the attributes of the temp_buff into a freshly made QgsVectorLayer, but so far no success.


texts:
It seems there used to be a GeoServer extension to support "Versioned Web Feature Services" (https://github.com/boundlessgeo/geoserver-2.1.x/tree/master/extension/wfsv).

Interestingly I can find no information about a versioning extension to the WFS standard, or about this GeoServer extension online.

So using GeoServer 2.13+, how do folks implement versioning or at least an audit log for WFS-T transactions? For example how do you track changes to a dataset maintained by multiple members of a team? This could then lead into publishing workflows etc, but I'm trying to wrap my head around how to implement a versioning workflow first ..


texts: Symbolizing a file geodatabase layer using every attribute field?

texts:
I'm trying to export an animation using ArcGIS 10.3.1. I have monthly unemployment rate data dating back to 1990 for every county in the continental US. Each month's unemployment rates are stored in a separate field.

In trying to create an animation, I have only been able to visualize this data in a way that assumes the timestamp is in the attribute, rather than in the field heading.

How can I visualize these data using each attribute field without creating a separate map by hand for every individual month?


texts: After upgrading ArcGIS Pro getting "Failed to Initialize Python Interpreter" error

texts:
I just upgraded ArcGIS Pro to version 2.2 from within the application. After the upgrade, when I open a project, this dialog comes up:

Also, the Python window shows this:

Finally, I cannot seem to see any toolboxes in the geoprocessing window.

Am I going to have to uninstall and reinstall to fix this? Or is there some other way?


texts: Determining if Polygon Contains Negative Space using ArcGIS SOAP Web Service?

texts:
I am running queries against a polygon layer using the ArcGIS SOAP Web Service.  A user selects a point within that layer and the resulting feature set contains PolygonN objects.

Occasionally, the polygons that come back contain negative space.  These are also known as donut or island polygons.  Here is an example:

In this case, I would only want the geometry contained in the purple area, not the green area.

When the user selects a point within the purple area, the resulting PolygonN object contains 2 rings, representing the outer and inner areas, as it should.

The inner ring, as far as I can tell, has no way to indicate that it is actually negative space.  So the resulting polygon looks like it has 2 solid shapes, when in fact it is only one shape with a hole in the center.

Does anyone know how I can determine that a polygon's ring is actually negative space?


texts: Adding Color Map to MrSID Imagery using QGIS

texts:
I am trying to convert a single band (grey) MrSid image to a 3 band (RGB) GeoTiff image. Using QGIS (3.12.3-București) on Windows 7, I saved it as a GeoTiff, which now has the 3 standard bands. Each is 1 to 255 min-max. However, the new GeoTiff image is still grey in color. The source from which I obtained the original image has 2 jpg files (1 is a low-res, grey overview), the other says it's a "layout" file and when viewed, just shows what appears to be the image of a form, though it does have the 3 RGB bands (I'm pretty sure it's not a color map). There are several xml and an shp file as well.
Is there a way to colorize the GeoTiff so that it appears more realistic?


texts: What do pixeltypes in PostGIS raster mean?

texts: I'm trying to build a raster in PostGIS from a vector layer and I was curious as to the meaning of pixeltypes in the raster geometry? How does it effect my raster?


texts:
I am new with QGIS. I know how to generate the x and y coordinates by using the field calculator ($x, $y). I know these are just static values. What I want to do is to change the coordinates by hand so that the values are actually changed.

Is there an easy way to do this?


texts: Style lines differently when intersecting a polygon

texts:
I want to give my contour lines a blue colour on glaciers (inside the polygon that shows the outlines of the glaciers) - and a brown colour everywhere else. What is the easiest way to achieve this? I am using ArcMap 10.5.1.

Here is my result using arctoolbox - analysis tools - extract - clip.


texts:
How can I modify the coordinates of a raster in the QGIS?

i.e. longitude and latitude (geographical positions)

I need to modify

X Min = 144964.482217
X Max = 150226.504610
Y Min = 4398796.699505
Y Max = 4406379.991327

to

min. X : 655295.0176745
max. X : 670395.0176745
min. Y : 4389577.5709399
max. Y : 4401677.5709399


texts: Lat-Lon-Tools -> no external map choosable

texts:
Since today (dont know why!) I can't select an external map in the Lat-Lon-Tools-Settings. I already uninstalled the AddIn, deleted the User-Folder and reinstalled it. Same: there is no option in the drop-down-field... ;(

From the log...
2020-08-15T20:57:28     WARNING    Traceback (most recent call last):
              File "C:/Users/xxxxx/AppData/Roaming/QGIS/QGIS3\profiles\default/python/plugins\latlontools\settings.py", line 552, in showEvent
              for item in settings.userMapProviders:
             TypeError: 'NoneType' object is not iterable

I think it is a bug (because there was a update). I installed an old version (latlontools-3.3.13) and it works)


texts:
First post, so be gentle.

Context - I'm doing some research into healthcare provision and I'd like to visualise the data. I'd like to be able to put these PCT boundaries

as an overlay to google maps and change their colours and so on to show data...  I'm running osx and I'm fairly IT friendly, but I've got no idea where to start on this and google is sending me in all manner of different diagrams with language I'm struggling to get my head around.  So I'm looking for any information people can give me on:

Recommendations of tools that run on OSX that might help
A general outline of the workflow I'm going to be going though
Places where I could put the finished files in an open-sourcy way so that other people might be able to make use of them?

Thanking you all kindly in advance...

EDIT - should point out that I'm happy to trace the lines with a tablet if that's what it takes - doesn't have to be postcode perfect just as long as it's recognizable...


texts: Esri ArcGIS Portal selector window not filtering data view

texts:
I have created a selector box in the header of a web app I am creating in ArcGIS Portal. It is an Ops Dashboard-type web app. I have all of the pie and serial charts responding to the map extent and the category selection (from the selector). But the data view is static. It will not filter for the categories I have selected. I know I should go to the Selector/Activity/select data layer as target and tell it to filter, but (here is the source of the problem) the data layer will not show up in this list unless the Selector is configured for Features (see image)

But when I go to Actions, there is no capability to add a filter:

But even if there was, I need grouped values in the configuration. Why are the right options not showing up and how do I fix this?


texts: Editing two layers topologically in QGIS?

texts:
I have different polygons in layer a_input. The Polygons are adjacent and connected inside the layers and it should be. The Polygons in layer a_input have to fit adjacent to a polygon of layer 2 b_input. I tried v.edit and v.overlay but didn't find the right options to get the result I wanted.

In the Picture, the green Polygon is the b_input and the black lines are the boundaries of the polygons in a_input. They have to fit the boundaries of the polygon in b_input. The red line is the result of v.edit.

How to solve such question in QGIS?


texts: Errors relating to merging point shapefiles

texts:
I have created two different point shapefiles which represent towns within a state. However I want to know the distance between both towns, but first I want to draw a line between them. Is this possible or should I merge both shapefiles together first?

I was trying to merge the point shapefiles together as one but I experienced an error.

  TypeError: list indices must be integers or slices, not str

vicksent = ".vick.shp"
towncent= ".town.shp"

from shapely.geometry import MultiPoint, point, shape, mapping

schema = {'geometry': 'Point',
         'properties': {}, }
with fiona.open('mew.shp', 'w', 'ESRI Shapefile', schema,     crs=from_epsg(29900)) as output:
   prop = {}
   centriod1= vicksent
   centriod2 = towncent
   output.write([centriod1,{'geometry': point, 'properties': prop}])
   output.write([centriod2, {'geometry': point, 'properties': prop}])


texts: Need to convert XY values to features in feature service using ArcGIS Server JavaScript API

texts:
I need to convert XY values in a text file which the user uploads into point features (not graphics) in a feature service using ArcGIS Server JavaScript API.
The ArcGIS server version is 10.1 and API version is 3.5