texts: GeoServer can't determine native CRS from PostGIS db with wkb_geometry column

texts:
I am trying to extract a banks layer from daily OSM extract pbf files for Nepal. After using osmosis and ogr2ogr, I was able to get it into a PostGIS Database.

Here are the commands I used:

osmosis --read-pbf-fast file="latest.pbf" --tf accept-nodes amenity=bank,atm,bureau_de_change --tf reject-ways --tf reject-relations --write-pbf file="banks.pbf"

ogr2ogr -f PGDump banks.sql banks.pbf -lco COLUMN_TYPES=other_tags=hstore --config OSM_CONFIG_FILE conf/banks.ini

Relevant bits of the ini file:

[points]
# common attributes
osm_id=yes
osm_version=no
osm_timestamp=no
osm_uid=no
osm_user=no
osm_changeset=no

# keys to report as OGR fields
attributes=name,barrier,highway,ref,address,is_in,place,man_made
# keys that, alone, are not significant enough to report a node as a OGR point
unsignificant=created_by,converted_by,source,time,ele
# keys that should NOT be reported in the "other_tags" field
ignore=created_by,converted_by,source,time,ele,note,openGeoDB:,fixme,FIXME
# uncomment to avoid creation of "other_tags" field
#other_tags=no

Here's the first column from the PostGIS database:

ogc_fid |                    wkb_geometry                    |  osm_id   |        name         | barrier | highway | ref | address | is_in | place | man_made |                 other_tags
---------+----------------------------------------------------+-----------+---------------------+---------+---------+-----+---------+-------+-------+----------+---------------------------------------------
       1 | 0101000020E61000007D29972169545540B24F00C5C8AE3B40 | 425306730 | Nabil Bank Pulchowk |         |         |     |         |       |       |          | "amenity"=>"bank", "operator"=>"Nabil Bank"

I then performed these steps:

Created a new store in my GeoServer and linked it to this PostGIS db.
Clicked on Layers > Add Layer and chose the new store.
Found the banks layer on it and clicked "Publish".

On the publishing page, the native SRS was not set and the declared SRS was EPSG 404000. The bounding box was also empty.

On clicking Compute from data, I got an error. Compute from native SRS didn't work at all.

Therefore, I set the declared SRS to EPSG 4326 and tried to compute from data. It returns -180 -90 180 90

With this, I published the layer, but I can't preview it on the GeoServer. I don't think the CRS has been applied correctly.

How do I set the Native CRS correctly once the data is in the PostGIS db?


texts: How to project the geographic grid in a shapefile in QGIS

texts: Is it possible to automatically project the geographic grid (View -> Decorations -> Grid) into a shapefile layer in QGIS?


texts: How to calculate the start and end lat/lon of a linestring which is not a part of a ST_buffer?

texts:
I'm trying to find out the lat/lon of points on a linestring that are not contained inside a geometry formed by a ST_buffer. I can find if the linestring intersects the geometry but I'm not able to find how can I get the points.

For example. Consider the following image. The substrings of the linestring that are not covered by the polygon are BC and FG. So, I'm trying to get the lat/lon of the B,C,F,G in order. Is it possible to do this in PostGIS ?

UPDATE:

I've written this query so far but it returns more than expected values. Some of them are inside the ST_buffer.

SELECT ST_AsText(ST_Difference(road::geometry,circle.the_geog::geometry)) FROM (SELECT ST_GeographyFromText('SRID=4326;LINESTRING(85.269150 23.157234,85.278541 23.101098,85.286367 23.047280,85.287737 23.006227)') as road) AS road,towers as circle WHERE ST_Intersects( circle.the_geog,ST_GeographyFromText('SRID=4326;LINESTRING(85.269150 23.157234,85.278541 23.101098,85.286367 23.047280,85.287737 23.006227)')
    );

UPDATE 2:

I also tried ST_SymDifference. Though the results are correct but I'm also getting POLYGON in the GEOMETRYCOLLECTION. Also it returns many rows instead of just one.

SELECT ST_AsText(ST_SymDifference(road::geometry,circle.the_geog::geometry))
FROM (SELECT ST_GeographyFromText('SRID=4326;LINESTRING(85.269150 23.157234,85.278541 23.101098,85.286367 23.047280,85.287737 23.006227)') as road)
AS road,poly as circle WHERE ST_Intersects(circle.the_geog,road);


texts: Geotiffs not displaying colour correctly in MapInfo

texts:
I'm trying to open geotiffs I made in ArcGIS 10.2 in MapInfo 12.5. They open in MapInfo in the right location but displayed on a blue colour scale. These should be colourful RGB images.

Trying to use the Discover enhance image tool just crashes MapInfo.

Any solutions?


texts: Is there a "Disperse Markers" option in Mapbox GL-JS?

texts:
I have used the "Disperse Markers" tool in ArcGIS before, but I am currently working on a project where I am using Mapbox GL-JS to display point data. Many of my points overlap, and I would like them to all be completely visible. This would be possible with "Disperse Markers" in ArcGIS or "Point Displacement" in QGIS.

Here's a link to the ArcGIS tool explanation: https://blogs.esri.com/esri/arcgis/2008/08/01/how-to-scatter-marker-symbols-for-stacked-or-clustered-point-features/


texts:
Still learning my way around QGIS and ran into an odd issue this evening. I'm using a program called Radio Mobile Deluxe to generate RF models.  It exports images as KML ground overlays.  I then bring the KML into QGIS, then georeference the raster.

However, this evening, I noticed that things didn't look quite right.  I should be looking at perfectly square objects (8000x8000 meters and a 800x800 pixel raster).

The overlay imports perfectly into Earth, but when I bring it into QGIS, it's measuring out 10km wide.  However, when I compare visual features in the imagery basemap and roads, everything looks dead on.  I'm really not sure what I'm looking at here, if it's a problem, or how to fix it.

I'm using QGIS 2.18.4, the project CRS is 4326.  Not quite sure what other information I should provide.

Here's a screenshot, showing a very square object in Earth and a not-so-square object in QGIS, along with the N/S/E/W coordinates generated by Radio Mobile.

Here's a sample KML file I created in Earth and loaded into QGIS:

<?xml version="1.0" encoding="UTF-8"?>
<kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2" xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">
<Document>
    <name>foo.kml</name>
    <Style id="s_ylw-pushpin_hl">
        <IconStyle>
            <scale>1.3</scale>
            <Icon>
                <href>http://maps.google.com/mapfiles/kml/pushpin/ylw-pushpin.png</href>
            </Icon>
            <hotSpot x="20" y="2" xunits="pixels" yunits="pixels"/>
        </IconStyle>
        <LineStyle>
            <color>ff0000ff</color>
            <width>3.5</width>
        </LineStyle>
        <PolyStyle>
            <fill>0</fill>
        </PolyStyle>
    </Style>
    <StyleMap id="m_ylw-pushpin">
        <Pair>
            <key>normal</key>
            <styleUrl>#s_ylw-pushpin</styleUrl>
        </Pair>
        <Pair>
            <key>highlight</key>
            <styleUrl>#s_ylw-pushpin_hl</styleUrl>
        </Pair>
    </StyleMap>
    <Style id="s_ylw-pushpin">
        <IconStyle>
            <scale>1.1</scale>
            <Icon>
                <href>http://maps.google.com/mapfiles/kml/pushpin/ylw-pushpin.png</href>
            </Icon>
            <hotSpot x="20" y="2" xunits="pixels" yunits="pixels"/>
        </IconStyle>
        <LineStyle>
            <color>ff0000ff</color>
            <width>3.5</width>
        </LineStyle>
        <PolyStyle>
            <fill>0</fill>
        </PolyStyle>
    </Style>
    <Placemark>
        <name>Untitled Polygon</name>
        <styleUrl>#m_ylw-pushpin</styleUrl>
        <Polygon>
            <tessellate>1</tessellate>
            <outerBoundaryIs>
                <LinearRing>
                    <coordinates>
                        -79.22887,38.08849,0 -79.13737,38.08849,0 -79.13737,38.08849,0 -79.13737,38.16048,0 -79.22887,38.16048,0 -79.22887,38.16048,0 -79.22887,38.16048,0 -79.22887,38.08849,0
                    </coordinates>
                </LinearRing>
            </outerBoundaryIs>
        </Polygon>
    </Placemark>
</Document>
</kml>

EDIT: Here's what the raster looks like (after converting to polygon and removing the background).

When I spot-check the overlay down to the pixel, it's identical between QGIS and Earth.  Accurate enough for modeling RF signals, which tend to get pretty fuzzy in the real world.


texts: Is there a shortcut to create only right angles while editing?

texts:
While digitizing in ArcGIS Pro, I am wondering if there is a way while creating segments to allow only right angles. I understand that this can be done manually by entering direction, deflection, or perpendicular (as well as keyboard shortcuts for those). I want to know if there is a quicker way to achieve this.

Ideally, I would like to have a key-binding that I could hold down while drawing line or polygon segments; it would constrain the next segment to be placed at a perpendicular offset (90 or -90).

I am digitizing a large number of buildings in a city (automatic extraction options not available). It would be too time-consuming to manually enter deflection angles throughout the whole process.


texts: Remove routing-container of leaflet-routing

texts:
How to remove the whole container (see picture below) of leaflet-routing that includes distance and travel time via a button in my index.html?

If i remove line 3508 of leaflet-routing-machine.js: return this._container; the container is not visible. But this way is not an option.


texts: Integer field name changes automatically

texts:
Last week I had imported some shapefiles into file geodatabase. It makes sense but when I checked table of the outputs I had found that some integer field names had been changed to other names. For example original name is 2010 changed to F1.

Please help.


texts:
I'm looking for a C# library that is able to extract from line shapefile all points of the selected line in a sequential manner from the start to the end of the line.

I know it's doable with Python and QGIS, but I'm looking for C# lib since it shall be connected to other modules of the system.


texts: Automatic insert of project name into database

texts:
For every layer I've a field in my postGIS database that should be filled with the name of the current project and today I do this manually.
Is there a way to do this with the Commands that can be defined in the style?


texts: How can I get a window by clicking on a mapsubject in the QGIS web client?

texts:
Is it possible in Qgis-web-client to show a new little window (with custom informations and maybe a little picture) by clicking on a mapsubject? How to modify Qgis-web-client if possible?
Maybe its possible to modify the attribute-window, that i get on the right side by using the identification cursor-mode?

Thanks


texts:
How can i send mail with MapBasic??

I have already tried this code, but it's not working
 i used SenMail Application

sub send
dim cmdstr, subj, msg as string

subj = "Hello from MapInfo-L" '  is the email subject
msg = "This is how you use sendemail to send an email" ' is the message text

cmdstr = "sendemail -f  example@gmail.com -t example@gmail.com -u " + subj + " -m " + msg
run program cmdstr
End Sub


texts:
I've been asked to do a project that maps the bathymetry of a quarry lake. It seems like I've seen a feature in QGIS that would help with this if I took a bunch of depth points in the lake. Any thoughts on where to find this? And how to do it?

And there's not a chance this data could already exist somewhere from remote sensing, soil surveys, etc.?

Thanks


texts: How to generate Report with Unicode Character for Gujarati text?

texts:
I am using ArcMap 10.1 and want to generate report using Report Designer.

In Report header i use Label Control for display unicode character (Gujarati) but when i copied gujarati text from word doc and assign to label text property, it display as question marks, 'garbled' or 'nonsense'.

Is there a simple way to display Gujarati characters in an ArcMap label?


texts:
I am new to QGIS.
I want to subtract a sequence of rasters (daily data): one day from the previous.
day100-day99; dayy99-day98; day98-day97 and so on, with file names that appropriately label what each output is.


texts: How to enlarge polygon or a group of polygons by specific amount to each polygon?

texts:
I am working on generalization of Geological maps. As an initial step of it I need to enlarge some of the polygons, but different polygon should have different enlargement amount. It is because minimum size constraint, which dictates to remove, enlarge or aggregate all the polygons which have smaller than the fixed amount. For my case, the amount equals to 2500 m2. I decided to remove some of the polygons and enlarge some of them.

I need to enlarge polygons just above the 2500 m2 so that they meet the minimum size constraint, but as different polygon has different area I need to add  different amount to each of them.

I am working with ArcGIS 10.4, also would like to implement python with arcpy library in order to automate it.

Here are some examples of my polygons and part of the table which shows the area of the polygons and how much should be add to each of them.

Polygons


texts: Tell Dojo to look for minified classes with ArcGIS JavaScript API

texts:
I've got an application build using the ArcGIS JavaScript API. I have several classes built with Dojo define. Each class is in a separate file, and everything is working fine.

However, I've minified my classes (using Visual Studio Web Essentials), and I'd like to tell Dojo to use the minified files. I'm using a dojoConfig like the following:

var dojoConfig = {
  isDebug: false,
  packages: [{
    name: 'js',
    location: '/js'
  }]
};

And dojo loads my files like '/js/MyClass.js'. What do I have to do the get Dojo to load MyClass.min.js ? This would be a successful resolution to my problem.

But, even more, I have gzipped versions of those files, so how can I tell Dojo to download MyClass.min.js.gzip?

I'd even like to take it a step further and bundle all of my classes into one download, say MyBundle.min.js.gzip. That not a problem. But, again, how to explain to the Dojo AMD loader that /js/MyClass1 and /js/MyClass2 are both in MyBundle?

All this, given that I'm not using the Dojo build system.


texts: Census Block GEOid2 conversion to Latitude Longitude

texts:
My question is similar to this question. I have two different data files.

Data1.csv (4500 rows)
GEOid                       GEOid2              Total Population    Urban   Rural
1000000US220019603001020    220019603001020     72                  0       72

Data2.csv (7000 rows)
Latitude    Longitude   ADT
30.48103    -92.61743   234

I want to get GEOid2 for each point (latitude, longitude) which I can do manually by using this API. This task will be then tiresome.

Is there any particular method by which I can use the list of latitude and longitude to get exact GEOid2 for each location point by using that API?

Finally I want my table like this.

Data.csv
Latitude    Longitude   ADT GEOid2            Total Population   Urban  Rural
30.48103    -92.61743   234 220019603001020   72                 0      72


texts: How to mosaic multiple Sentinel-1 images in SNAP software?

texts: I am trying to mosaic multiple sentinel sar images in esa snap software. The scenes are acquired from two adjacent tracks. I used the tool S1 slice wizard in snap under radar tool menu. But my images didn't get mosaiced. I am attaching the screenshot of the output. Kindly advice.


texts: GeoJSON - how to style features of GeometryCollection

texts:
I have a GeoJSON file with a feature as follows:

{
  "type": "Feature",
  "properties": {
    "Start Date": "01-09-2015",
    "End Date": "01-09-2015"
  },
  "geometry": {
    "type": "GeometryCollection",
    "geometries": [
      {
        "type": "Point",
        "name": "test",
        "coordinates": [74.86083984375, -52.94201777829492]
      },
      {
        "type": "LineString",
        "coordinates": [[74.86083984375, -52.94201777829492], [75.5419921875, -53.10721669189339]]
      },
      {
        "type": "Point",
        "coordinates": [75.5419921875, -53.10721669189339]
      }
    ]
  }
},

As you can see it consists of 2 points with a line joining the two points.

I want to style the 2 points with different colors but am stumped as to how to do that. At the moment I have the following ....

        var styleEnd = [
            new ol.style.Style({
                image: new ol.style.Circle({
                    radius: 7,
                    fill: new ol.style.Fill({
                        color: '#FF4500',
                        opacity: 0.8
                    })
                })
            })
        ];

        var styles = {
            'Point': [new ol.style.Style({
                    image: new ol.style.Circle({
                        radius: 7,
                        fill: new ol.style.Fill({
                            color: 'black'
                        })
                    })
                })],
            'LineString': [new ol.style.Style({
                    stroke: new ol.style.Stroke({
                        color: 'green',
                        width: 1
                    })
                })],
            'MultiLineString': [new ol.style.Style({
                    stroke: new ol.style.Stroke({
                        color: 'green',
                        width: 1
                    })
                })],
            'MultiPolygon': [new ol.style.Style({
                    stroke: new ol.style.Stroke({
                        color: 'yellow',
                        width: 1
                    }),
                    fill: new ol.style.Fill({
                        color: 'rgba(255, 255, 0, 0.1)'
                    })
                })],
            'Polygon': [new ol.style.Style({
                    stroke: new ol.style.Stroke({
                        color: 'blue',
                        lineDash: [4],
                        width: 3
                    }),
                    fill: new ol.style.Fill({
                        color: 'rgba(0, 0, 255, 0.1)'
                    })
                })],
            'GeometryCollection': [new ol.style.Style({
                    stroke: new ol.style.Stroke({
                        color: 'magenta',
                        width: 2
                    }),
                    fill: new ol.style.Fill({
                        color: 'magenta'
                    }),
                    image: new ol.style.Circle({
                        radius: 10,
                        fill: null,
                        stroke: new ol.style.Stroke({
                            color: 'magenta'
                        })
                    })
                })],
            'Circle': [new ol.style.Style({
                    stroke: new ol.style.Stroke({
                        color: 'red',
                        width: 2
                    }),
                    fill: new ol.style.Fill({
                        color: 'rgba(255,0,0,0.2)'
                    })
                })]
        };

        var styleFunction = function (feature, resolution) {
            if (feature.getGeometryName() === 'test') {
                return styleEnd;
            } else {
                return styles[feature.getGeometry().getType()];
            }
        };

        var vectorLayer = new ol.layer.Vector({
            source: new ol.source.Vector({
                url: 'map.geojson',
                format: new ol.format.GeoJSON()
            }),
            style: styleFunction
        });


texts: Size of cost raster in Least Cost Path plugin

texts:
When I run the Least Cost Path plugin in QGIS, its success is dependant upon the dimensions of the 'cost raster'. If the original DEM or subsequent 'cost raster' it too large then the operation seems to fail.
With the successful example below, can anyone tell me what the "8765 * 7727" actually is/means?? It is not metres or kilometres. Is it pixels???

QGIS version: 3.8.3-Zanzibar
QGIS code revision: 685d8b15d2
Qt version: 5.11.2
GDAL version: 2.4.1
GEOS version: 3.8.0-CAPI-1.13.1
PROJ version: Rel. 5.2.0, September 15th, 2018
Processing algorithm…
Algorithm 'Least Cost Path' starting…
Input parameters:
{ 'BOOLEAN_OUTPUT_LINEAR_REFERENCE' : False, 'INPUT_COST_RASTER' : 'E:/BU Volunteer/QGIS/Ians-Rift-Valley/10-overlayed-maps/least-cost-stuff/adama.tif', 'INPUT_END_LAYER' : 'E:\\BU Volunteer\\QGIS\\Ians-Rift-Valley\\10-overlayed-maps\\least-cost-stuff\\adama.gpkg|layername=adama', 'INPUT_RASTER_BAND' : 1, 'INPUT_START_LAYER' : 'E:/BU Volunteer/QGIS/Ians-Rift-Valley/10-overlayed-maps/least-cost-stuff/awash-start-4326.gpkg|layername=awash-start-4326', 'OUTPUT' : 'TEMPORARY_OUTPUT' }

The size of cost raster is: 8765 * 7727
Searching least cost path...
Search completed! Saving path...
Execution completed in 4058.56 seconds
Results:
{'OUTPUT': 'Output_least_cost_path_a74cf0f6_b4aa_4ee6_a876_ef6c26f4e164'}

Loading resulting layers
Algorithm 'Least Cost Path' finished


texts: Get Area size using Dragbox in Openlayers 3

texts:
I need to be able to show the square meters of an area using dragbox.  There's a great example using draw polygon, but I need to limit users to a regular rectangle.

The formulas and SHOULD be the same or simular. I've replaced "Draw" with "DragBox" and "drawstart" and "drawend" with "boxstart" and "boxend".  I also commented out the "formatLength" function, as I don't need it.

I'd use this script if I knew how to limit draw polygon to a regular rectangle.


texts: QGIS Server, GetCapabilities not working on Ubuntu 14.04

texts:
I've spent a couple of days now troubleshooting a QGIS Server setup in Ubuntu 14.04.

I am following this guide to install QGIS server:
http://anitagraser.com/2012/03/30/qgis-server-on-ubuntu-step-by-step/

According to the guide above, the installation should respond to the request below.
wget -p -O - "http://localhost/cgi-bin/qgis_mapserv.fcgi?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetCapabilities"

All I am getting is a 404. Apache 2 and FCGID are installed and enabled. I seem to be missing something.

I've tried the setup with 2.8LTR and 2.14LTR both unsuccessfully, but unable to solve this issue.


texts: Setting Python script tool parameter type so that parameter accepts feature class or layer as input?

texts: I would like a few of my parameters in my script tool to be able to accept either a feature class or feature layer as input.  However, I have not been able to figure out a proper parameter type to allow this.  Can anyone help me out on this seemingly simple question?


texts: Convert cluster points into a single point

texts: I need to try and group cluster points by attribute into one point.  So all points with the same attribute need to be grouped onto 1 location.


texts: Can ArcGIS Desktop 9.3.1 open Geodatabase created using ArcGIS for Server 10.1?

texts:
can it possible to use old version of arcgIS desktop, say ver 9.3.1, if we develop a geodatabase using arcGIS 10.1. The above answers were very useful will make me to understand more on arcgis server, since we are planning Thanks

Suplimentary Question: may be question is not understood well..we have old version of arcgis desktop 9.3.1, now we proposing and wanted to develop geodatabase system using new version of arcgis server i.e 10.1, will it compatible.


texts:
I have been using osrm package for a while and didn't ever come across this kind of issue before. It looks so weird and I have been trying to update package and reinstall R as well. But nothing is working. Can someone, guide where am I doing wrong, please.
Code:
library(sf)
data("berlin")
# Get isochones with lon/lat coordinates
iso <- osrmIsochrone(loc = c(13.43,52.47), breaks = seq(0,14,2),
                     returnclass="sf")

Error:
Error: lexical error: invalid char in json text.
                                       <!DOCTYPE HTML PUBLIC "-//IETF/
                     (right here) ------^

The OSRM server returned an error:
Error: lexical error: invalid char in json text.
                                       <!DOCTYPE HTML PUBLIC "-//IETF/
                     (right here) ------^

The OSRM server returned an error:
Error: lexical error: invalid char in json text.
                                       <!DOCTYPE HTML PUBLIC "-//IETF/
                     (right here) ------^

Error in UseMethod("st_as_sf") :
  no applicable method for 'st_as_sf' applied to an object of class "NULL" ```


texts: Create Boundary Layer on FeatureLayer Geotools

texts:
I want to create a boundary layer on my FeatureLayer in Geotools is there any direct option to create this.
Following Image shows what i want to achieve this.

Please suggest .
we have done this in DotNET by using ConvexHull algorithm .


texts: field calculator python script with coded values

texts:
I am trying to run a field calculator python script that will take a value from one field in the attributes table and create an abbreviation in another field that was designated for it.

Here is a example of the code that was written for this function:

def checkValue(c,d):
    if c == 'ARTERIAL_CRACKSEAL':
        return 'ACS'

The problem is that once I created coded value domains for these attributes the code stopped working (i.e. code:1 description: ARTERIAL_CRACKSEAL ; code:2 description ARTERIAL_microseal ; etc.).  Would anyone know how to adjust the code so the script runs successfully? I need to keep coded value domains for this attribute field.


texts:
I was wondering if anyone would help me on this: I have mapinfo pro and a raster .grd file that contains signals strength levels for different locations. All i want to do is to simply calculate the total area of the .grid file in question. Would there be an easy and quick way to obtain the result?
Thanks in advance,


texts: What are some self-hosted alternatives to Graphhopper's Geocoding API?

texts:
I have a self-hosted instance of Graphhopper. However, it doesn't come with autocomplete like the web demo; from doing some searching, it looks like the autocomplete comes from the Geocoding API, which is a part of GH's Directions API. GH offers the Directions API for business use only, with no self-hosting option available. Are there any geocoding APIs that can:

Integrate with OSM data?
Have a self-hosted option? and
Integrate with Graphhopper, without too much hair-pulling?

If there are self-hosted alternatives for the other components of the overarching Directions API (Routing API, Map Match API, etc.), so much the better.


texts: online map of 242,888 locations: which service to use? (what if temporal heat map?)

texts:
I have a 161.8 MB large CSV file with geocoded locations (with latitude and longitude) for 242,888 lines. I want to have an embeddable online map of all this. I wonder which service could deal with this amount of data, and what you'd recommend.

Relevant details:

— if the data looks too big (162 MB), I can cut down on fields, some are extraneous. But others would be nice to link to the points on the map, even if they are not strictly there for the mapping. (So people could navigate back to the original data source, e.g.)

— I have date field for a temporal map. It would be great to use this to let viewers filter which points to show, or even better, show animated versions if possible.

— one field would be good to use for showing intensity, hot spots on the map (viz. concentration of procurement contracts in space and time)

— the data is the universe of procurement contracts in Hungary since 1998. So good Unicode support is a must, as is global coverage (though I think both are pretty standard these days).


texts: Does r.in.gdal -o reproject the imported raster data to the location projection?

texts: Does r.in.gdal -o actually reproject the imported raster data to the location projection, or do I need to call gdalwarp to reproject prior to importing?


texts: Plotting layer given with Longitudes [0 360] to [-180 180]?

texts: I have a world map in correct Longitudes from -180° to 180°. Additionally, I have a shapefile of data from South America, which shows at 290° instead of -70°. How can I reproject it using QGIS?


texts: QGIS Raster Visualization Problems from a self-produced .asc file

texts:
I have been creating .asc raster files by saving a numpy array with numpy.savetxt and the following header:
header = "ncols        %s\n" % ncols
header += "nrows        %s\n" % nrows
header += "xllcorner    %s\n" % xll
header += "yllcorner    %s\n" % yll
header += "cellsize     %s\n" % cellsize
header += "NODATA_value  -99999"

After producing the .asc files I load them into QGIS and they are by default given a style with min: nan, and max: nan. If I go into the Style dialog, set to Stretch to min/max, and change the min and max to match the appropriate range for the source file (eg 0 to 40), the preview thumbnail shows a pixelated mess and the raster shows as completely black.

Sometimes the raster shows up fine if I just hit apply, but then reverts to all black once I hit OK. If I am zoomed very far in (well within the limits of the raster), then sometimes the raster shows up fine and sometimes it doesn't.

Source and generated .asc files here:
https://github.com/mrlmautner/PrecipRasters


texts: Can I set the max or min zoom on a CartoDB embed?

texts:
I'm using the default iframe embed to display two CartoDB maps, but I'd like to prevent people from zooming out past a certain point, or if that isn't possible, to prevent zooming at all.

Can I do that? I sniffed around but I don't see anything obvious. I tried adding  &maxzoom=5 to the URL, but that doesn't seem to do anything.


texts: Combine columns passed as parameters in QGIS graphical modeller

texts:
I would like to combine two fields in QGIS using graphical modeller. If I were to do it using field calculator, I would do it as follows:

"field1" + ", " + "field2"

Unfortunately, I cannot do it, as I want those fields to be passed as parameters in the graphical modeller.

I have the following model:

As you can see, I have two parameters: "Calle" and "Portal", which are fields of the input vector layer. I would like to combine both fields using the field calculator algorithm. I need to pass the fields' names as parameters as I am looking for a generic solution that may work for different field names.

I have tried with the following formula in the field calculator algorigthm:

Unfortunately, instead of combining the parameters' values, it combines the parameters' names. Therefore, I tried with the parameter() function, which I assume should pass not the parameter itself but its value, as follows:

Unfortunately, when doing so, I get an empty result.

Is there any way to achieve what I am looking for?


texts:
I have standalone application (ArcObjects C# 10.3) uses axSceneControl, I'm trying to get a mouse click point:

private void M_axSceneControl_MouseDown(object sender, ISceneControlEvents_OnMouseDownEvent e)
{
    var m_pSG = m_axSceneControl.SceneGraph;
    var m_pSW = m_pSG.ActiveViewer;
    object a1;
    object b;
    IPoint pt;
    m_pSG.Locate(m_pSW, e.x, e.y, esriScenePickMode.esriScenePickGraphics, true, out pt, out a1, out b);
}

The point is under the mouse pointer with the filled Z coordinate.
But if you get the height at the specified point pt through the
function ISurface.GetElevation(pt.x, pt.y), (IRasterSurface, ZFactor 0.001). I get a value close to Z but not exactly accurate.

X=78.2967, Y=31.8368, Z=2.2807; Elevation=5159.5076
X=77.3918, Y=32.6779, Z=2.2474; Elevation=5732.2082
X=73.3516, Y=30.4610, Z=-0.0559; Elevation=151.3352
X=71.7634, Y=29.0155, Z=-0.0909; Elevation=114.8338

why is this happening ?


texts:
I am trying to split an attribute into its component pieces with FME.

Ex: 1707 WILLINGTON ST, PHILADELPHIA PA 19121

into: 1707 WILLINGTON ST, PHILADELPHIA, PA, and 19121

I am using the AttributeSplitter using a comma Delimiter and I get this:
1707 WILLINGTON ST and PHILADELPHIA PA 19121

How do I get the state - city - ZIP in three separate fields?


texts: Limit the number of lines in a stacked label using python parser in ArcMap labels

texts:
I am labeling names of places by stacking the labels at each space using the python code below. This makes some incoherent labels for text with many spaces and short or abbreviated words.

def FindLabel ( [Label] ):
 S = [Label]
 S = S.replace(' ', '\n')
 return S

I would like to limit the number of stacked lines that are allowed in each label (probably only 2 or 3).

I got it to work fine using the vb code below but due to a linux server it seems that only python will work for us in this case.

left( [Label]   ,instr((len( [Label] ))/2, [Label]  ," ")) & vbnewline &
right( [Label]  ,len( [Label]  )-(instr((len( [Label] ))/2, [Label]  ," ")))

Is there a way to directly define the number of lines allowed using the python parser, or is it best to count the characters like the vb expression did?


texts:
How do I know what projection to use for north-up or 0 rotation? I have a gtiff (ND83-nf) that gives me the error when importing to Grass. Should I re-project or recreate the image with the rotation burned into the image? If I re-project what projection should I use? when  I re-project in fwtools the default comes up as the same as the input. not sure what to choose.
Thanks


texts: Change in ArcGIS locator service causing error

texts:
I have ran into an odd problem the last few days.

I am using ArcGIS for Runtime and ArcGIS for Server (which is located on a remote machine), as well are ArcMap to publish my services. My services are publishing fine and look OK in my server manager, but when I use my C# application containing the URL for the feature layer I am receiving the message:

Locator service failed: System.Net.WebExceptions: The remote server returned an error: (403) Forbidden.
 at System.Net.HttpWebRequest.EndGetResponse(IAsyncResult asyncResult)
 at System.Net.WebClient.GetWebResponse(WebRequest request, IAsyncResult result)
 at System.Net.WebClient.DownloadBitsResponseCallback(IAsyncResult result)

Here is an image of the error:

The odd thing about this problem is that the data is still loaded from the ArcGIS server and I still have all Feature Access capabilities.

My URL also looks like this:

Url="http://**.*.*.***:6080/arcgis/rest/services/DicksonCounty/FeatureServer/0"

So I figured out that the issue was with ArcGIS changing their locator service. I have a ArcGIS online subscription and I can still not figure out what site I am supposed to be using now instead of the one that is listed in all of their documentation.

http://tasks.arcgisonline.com/ArcGIS/rest/services/Locators/TA_Streets_US/GeocodeServer

Has anyone else experienced this?

How should I fix the issue.

Code for locating:

Locator FindlocatorTask = new Locator("http://geocode.arcgis.com/arcgis/rest/services/World/GeocodeServer");
                FindlocatorTask.AddressToLocationsCompleted += FindLocatorTask_AddressToLocationsCompleted;
                FindlocatorTask.Failed += FindLocatorTask_Failed;
                AddressToLocationsParameters addressParams = new AddressToLocationsParameters();
                Dictionary<string, string> address = addressParams.Address;

                if (!string.IsNullOrEmpty(Address.Text))
                    address.Add("Street", Address.Text);
                if (!string.IsNullOrEmpty(City.Text))
                    address.Add("City", City.Text);
                if (!string.IsNullOrEmpty(State.Text))
                    address.Add("State", State.Text);
                if (!string.IsNullOrEmpty(Zip.Text))
                    address.Add("ZIP", Zip.Text);

                FindlocatorTask.AddressToLocationsAsync(addressParams);
                PortalGrid.Visibility = Visibility.Visible;
                AddressGrid.Visibility = Visibility.Collapsed;


texts: GDAL framework does not exist after KyngChaos installation of QGIS

texts:
I'm on a Macbook running Sierra 10.12, so I used the KyngChaos installation of QGIS 3.8. I followed the instructions in the ReadMe, and installed python (though I already had it on my machine), GDAL and then QGIS. I took notice of the 'known issue' as well and created a prepend PATH variable so QGIS could load the GDAL framework, but it did not work. After searching the app's package contents, I realized that there is no "GDAL.framework" folder on my machine even though I downloaded it from the instructions.

I used EasyFind and noticed that I have gdal installed in anaconda, can I use this path for QGIS to be able to load the module? If not, how can I download GDAL so that QGIS can recognize it and it shows up on my machine? I'm not sure why the KyngChaos installation is broken (it's worth noting that the path in the ReadMe file is incorrect as well).

EDIT: Here's a screenshot of the frameworks folder within the app's package contents, alphabetically sorted. As you can see, there is no GDAL.framework.


texts: Anyone who know if there is a method of creating rivers network from digital elevation model in qgis as it does in contour.


texts: How to calculate LS Factor of Watershed for RUSLE Calculation?

texts:
I managed to get subbasins for my catchment (7 km^2) with their areas but I haven't figured out how to get the LS factor (see Revised Universal Soil Loss Equation) from it.

I have the vector file for my subbasins and I have the DEM (digital elevation model) file. I extracted the LS raster file from the DEM file (via r.watershed) but I dont know what to do afterwards. I figured I can isolate the subbasins within the LS raster file and calculate the average LS factor of those areas.
Do you know how to do this?

On a side note, is there a way to split up the subbasins since there are slopes that are steep with slopes that are flat? Thanks in advance.


texts: Issue while configuring DNS with Portal for ArcGIS

texts:
Steps followed:-

Portal for ArcGIS 10.4.1 (Linux SUSE) Installation
ArcGIS Web Adaptor (Linux SUSE) Installation
Configuration of Portal with Web Adaptor (Web adaptor is configured on Apache Tomcat 7.065):

Portal is accessible using fully qualified machine name.
Portal is further configured with SAML: Portal is accessible through SAML successfully

DNS Configuration in Web Adaptor to access portal using DNS name: Portal is accessible using DNS successfully.

Problem:-

One of our web application (build using web app builder developer edition 2.2) is using spatial analysis widget (Buffer) and requires authentication for execution.  Without DNS i.e., using fully qualified machine name (machinename.domain.com) it is giving login option for both SAML and Portal identify store but with DNS (appname.domain.com) the authentication option is only portal identity store and no SAML?

As we have configured DNS after portal installation and ESRI documentation says for assigning friendly portal name, assign the name to the portal machine before you install the software. Then use the friendly name when setting up the portal. Do we have to change portal machine name before Portal for ArcGIS installation?


texts: How is GIS used in landscape architecture?

texts:
I was wandering, if out there somewhere are people that use GIS extensively in the field of landscape architecture (or closely related) and how do they use it?

May be some examples of work, maps, visualisations?

I am landscape architect, but at the time I am using GIS mainly in the field of hydrology.

In the past I did some visibility analyses with use of MapInfo and Vertical Mapper extension and also some landscape character analyses/mapping.

Today I am working mainly with ESRI and OpenSource and I do know that capabilities of GIS are much more than visibility analyses. In example in landscape ecology (which is not exactly landscape architecture but there is some remote relation) there are Fragstat and other software that help computing landscape indices also we can build 3D models expressing actual situation more accurate then ever before. There are some examples of use GIS in landscape architecture that I know, but how in general others landscape professionals do use it? What software do they use?

Landscape architecture is a broad and interdisciplinary field. In different countries understood in different way. I do have this idea, that in order to get better understanding of my own profession it would be great to hear about use GIS in landscape research/analysis/design/planning etc.


texts: Triggering mouse click event within search widget of ArcGIS Online

texts:
I copied this from GeoNet community.
I wrote a function before which is to check whether a point where a user clicks intersects with a layer, like:
view.on("click", function(event)
{
... //my handling of the event inside
queryEarthquakeLayer(event).then(function(y)
...
}
I recently added a search widget so that a user can type in the coordinates/address. How to mimic a mouse
click event inside the function handling the returned results from search:
search.on('search-complete', function(result){
...
//how to trigger a click event here so that I can use my previous function directly by just
//feeding the click event?
}
I found this online but it doesn't work:
map.emit('click', {mapPoint: new esri.geometry.Point(myLong,myLat, esri.SpatialReference({wkid:102100}))
Any idea?


texts: Maxent Fatal error: java.lang.ArrayIndexOutOfBoundExceptions: 44

texts:
I have come across a fatal error for some of the species I am modeling:
Fatal error: java.lang.ArrayIndexOutOfBoundExceptions: 44 (or other numbers).

The error only occurs when I am selecting "Create response curves" and it does not happen for all of
the species I am modelling.

I am using the same climatic and land cover layers for all species.

I have checked the data multiple times in ArcMap and R and all the ocurrence points are within the
extents of the raster data.

The error happens at the very end when Maxent is writing the average files. The single models can be
viewed in the output after the fatal errors shows, but not the averaged data.

If I do not select "create response curves" I get the usual output, minus the response curves
obviously.

Maxent settings are default but with 15 replicates, 10000 max. iterations and ticking "Create response
curves" and "Perform Jackknife tests".

Has anyone come across the same error and what could be the reason for it?


texts: Accessing cell value of a WMS raster served as ol.layer.Image from OpenLayers

texts:
I would like to access the value of the raster cell on clicked location (sort of an identify tool) of a raster displayed in a webapp as ol.layer.Image via MapServer. I have set the .map file of the MapServer to include TEMPLATE link and template.html.

Can anyone help me with the OpenLayers code to access and display the value of the cell?


texts: Concerning "Landsat 5's cloud mask and pixel quality assessment" in Google Earth Engine

texts:
I need assistance for solving a problem with Landsat 5 surface reflectance image's cloud mask.

I am new to this tool/programming and want to try this tool to process cloud mask. Through the provided open source code in the web page of Goole Earth Engine, the processing code is composed as below.

However, I would like to request help to debug this error message "Error: Invalid argument: 'image' must be of type Image."

var cloudy_scene = ee.Image('LANDSAT/LT5_SR/LT5*************');

Map.centerObject(cloudy_scene);

Map.addLayer(cloudy_scene,{},'all_bands_cloud');

// Load a Landsat 5 image collection.
var collection = ee.ImageCollection("LANDSAT/LT5_SR")

// Filter to get a geometry representing an export region.
.filterBounds(ee.Geometry.Rectangle([***.***,**.***,***.***,**.***]))

// Filter to get the years/months/days of data.
.filterDate('1993-01-01', '1993-03-01')  // filter to your favorite period

// Sort the collection in chronological order.
.sort('system:time_start', true);

//Details for cfmask
//0=clear
//1=water
//2=shadow
//3=snow
//4=cloud

var masked = collection.map(function(image){
    return image.updateMask(cloudy_scene.select(['cfmask']).neq(4))
});

// Apply the mask to the image and display the result.
Map.addLayer(masked, {}, 'masked');

var geometry = ee.Geometry.Rectangle([***.***,**.***,***.***,**.***]);

//export the image, specifying scale and region.
Export.image.toDrive({
    image: masked,
    description: 'ImageToDrive',
    scale: 30,
    region: geometry
});


texts: I am trying to load a coverage map onto Arcgis for desktop using .adf files that I downloaded from a county supported filezilla server. I am not experienced with coverage maps, however I have reason to believe that I do not have all the files needed to form a coverage map. When I try to search for the map using catalog nothing shows up and when I tried to import the files using data interoperability and it said that there was no geometry (however I could be using this tool wrong). I have uploaded a picture of the files that I downloaded, am I missing some files, or am I just misunderstanding the use of coverage map format?


texts: ArcMap - Export raster symbology to a raster function template programmatically

texts:
I would like to programmatically export symbology for a raster dataset as a raster function template xml file, as can be done by clicking the export button in the symbology tab of the layer properties dialog in ArcMap 10.2.2.

I wish to do this in a batch environment (outside of ArcMap) using ArcObjects.

Subsequently ESRI Canada support advised that the export required above is not exposed in the ArcObjects API. However I have found that xml serializing the RasterRenderer for the layer, snippet here, does in fact produce the same content as seen as the value for the Renderer variable in the manually / UI exported raster function template (!).

I'm trying to bring the symbology across for individual rasters after they are added to a mosaic dataset. Color maps seem to come across, but stretches don't. So after realizing the above, I'm working on wrapping the RasterRender serialized XML to be a raster function template and then applying it.


texts:
When I try and open psql from the postgis menu in opengeosuite 4 on the Mac I get the following error message.

“psql” can’t be opened because PostGIS is not allowed to open documents in Terminal.

Anyone know how to cure this error.


texts: What Projected or Geographic Coordinate System should I use to calculate km distances between various points on a large scale map?

texts:
I want to measure distances using the distance matrix tool on QGIS between the yellow dots and their closest blue dot. I want my output to be in KM and to be as accurate as possible. I don't want to use multiple different zones like WGS84 /UTM zone 9N etc. Rather just a single system which will give accurate / relativly accurate results.

This page has helped me to some degree (https://support.esri.com/en/technical-article/000006113) but I cant find the PDF to help decide which system to use.

I have tried using the distance matrix tool using WGS84 Psuedo Mercator but get my distances in Decimal Degrees which I've been told is nonsensical.


texts: Value changes while converting Geostatistical layer to grid using ArcGIS Geostatistical Analyst?

texts:
I'm trying to export the GA layer (krigging interpolated temperature data) to raster. While converting, the data range is getting changed in the resultant raster (max and min) as indicated in the image.

Can someone help me out?


texts: Merge function doesn't have output for all FCs, It has just first FC

texts:
The attached file is my script for merging different FCs as they get dissolve and join by the input FC. I don’t know why the output for merge does not include all FCs.

import arcpy, os
from arcpy import env
arcpy.env.overwriteOutput = True
env.workspace = arcpy.GetParameterAsText(0)
output = arcpy.GetParameterAsText(1)
gdbname=arcpy.ListWorkspaces()
for gdb_Name in gdbname:
  env.workspace = gdb_Name
  for fc in arcpy.ListFeatureClasses():
  hh=fc
  gg=arcpy.FeatureClassToFeatureClass_conversion(fc, 'in_memory', '{}'.format(fc) )
  kk=arcpy.Dissolve_management(fc,gg,"Route_Type","","","")
  arcpy.JoinField_management (kk,"Route_Type",hh,"Route_Type",'#')

fClasses = []   #initialize array to hold feature classes
fClasses.append(gg)
arcpy.Merge_management(fClasses, "Merged")
arcpy.FeatureClassToFeatureClass_conversion(fClasses, output, '#' )


texts: Restrict when point labels are displayed using PyQGIS

texts:
I have a number of points that I'm visualising on a points layer, and labelling using PyQGIS.  Unfortunately when zoomed out too much, the resultant map appears a mess of labels.
Can I use PyQGIS to specify that labels should only appear in certain circumstances? I was hoping to test with displaying labels:

Only if there are less than a certain number of points currently being rendered/visible on the screen  (e.g. 50) - I suspect this is the best solution
Failing that, set to only display labels if the scale is less than 1:5000*

Is it possible to do either of the above using PyQGIS?
Edit: * I can do this manually in QGIS by enabling: Rendering > Label options > Scale dependent visibility and setting min/max scale values
Current code is along the lines of:
layer_settings = QgsPalLayerSettings()
layer_settings.drawLabels = True
layer_settings.fieldName = "Vehicle"
layer_settings.isExpression = True

label_settings = QgsVectorLayerSimpleLabeling(layer_settings)

points_layer.setLabeling(label_settings)
points_layer.setLabelsEnabled(True)


texts: Find leaflet id of the latest layer created in Leaflet

texts:
I am adding layers to a layergroup "drawnitems" in leaflet

drawnItems.addLayer(L.polygon(polygon));

I need to find the leaflet id of the layer created using the above code immediately upon its creation. This happens in a loop. I am adding multiple layers.

My goal is to save shapes that have been edited by the user into the DB. These shapes are previously created and then stored in the DB and then rendered on the map. For this I am planning to using leaflet id assigned to each shape and then find the corresponding db entry using some logic in the "draw:edited" event.


texts: Unable to download osm2po from the original site

texts:
I've been trying to download the osm2po tool, in order to import the openstreetmap data into postgis database, from the original site and is impossible for me. All links I find in the page are cyclic, return to the main page and I can't find the way to access the tool. Any body can help me?

Thanks in advance.


texts: Extracting xyz coordinates from DGM file using gdal2xyz (I want the x and y coordinates in degrees not meters)

texts:
I am new to using QGIS and I am trying to compare elevation information between measured points using a total station and elevation given by the DGM files. The elevation comparison will occur in matlab. My problem is that the total station information gives the x and y coordinates in degrees/minutes and using the gdal2xyz in QGIS is only allowing me to extract the x and y coordinates in meters. This prevents me from directly comparing the points taken by the Total Station with the points in the DGM since the x and y coordinate are in degree/minutes and meters respectively.

My question is:

Is there a way to extract the xyz coordinates from the DGM file so that the x and y coordinates are shown in degrees/minutes and not in meters? if not, does anyone have a way around this?


texts: Is it possible to display an OpenLayers map upside-down?

texts:
This is more fun than a serious question, but I'd like to know if it's possible to display an OpenLayers map upside-down as per the uncommon, but interesting upside-down maps.

I guess it's mostly a matter of finding a source of the tiles and turning them 180 degrees, but I figure it's something someone has already thought through.  I'm not the only one who's interested in the stories that are (or are not) being told through maps.


texts: Text in layer toggler in Leaflet is really faint - why?

texts:
I have added a layer toggler using Leaflet and geojson but the text of layer names is too light to see in the toggler - how do I style the text in a darker colour? It's on a wordpress site, if that makes a difference?

I used the usual
    L.control.layers(baseMaps, overlays).addTo(map);


texts: Calculating nearest neighbor within each polygon using QGIS?

texts:
I have 3 layers:

world map divided into countries;
layer of points representing cities;
layer of points representing conflict events (e.g. firefight).

My goal: For each city, calculate the nearest conflict event, conditional on the city and the conflict evenT being in the same country. Assume that all countries have at least one conflict event and one city.

I know this is a "nearest neighbor" calculation, but I could not find an example where this calculation is done within polygons. When I try to do this without taking the world map (country polygons) into account, sometimes it matches cities and conflicts in different countries, which is not what I want.

I use QGIS 3.8.


texts: Can TileStache use external tilestache.cfg file?

texts:
I have a need for an ever evolving tilestache.cfg file.  I would like to be able to load this from an external source.  Is it possible to run a Tilestache server with

tilestache-server.py -c http://s3.amazonaws.com/*myBucket*/tilestache.cfg

or even with

`gunicorn --workers 4 "TileStache:WSGITileServer('http://s3.amazonaws.com/myBucket/tilestache.cfg')"

Or am I completely off base here?


texts:
According to the CartoDB tutorial on how to add and edit geometries, using the CartoDB editor, when we click Add Feature, we should get a popup box with three options:

Add point
Add line
Add polygon

But when I click that, I only get the Add point option.
Have the other two options moved somewhere else?
Or is there another way to add a line?


texts: Setting up pre-requisites for creating viewer using suite-sdk

texts:
I am creating a Viewer using the suite-sdk from the following tutorial, which says

The Client SDK should be installed onto a machine where the development will occur. This is typically bundled with the OpenGeo Suite itself (under the sdk subdirectory)

I have installed the OpenGeo and so I wanted to know,

where is the sdk sub directory in the Opengeo suite that i have installed located in Ubuntu?
I have JDK version,

java version "1.6.0_24"
OpenJDK Runtime Environment (IcedTea6 1.11.1) (6b24-1.11.1-4ubuntu3)
OpenJDK 64-Bit Server VM (build 20.0-b12, mixed mode)

ant -version -> 1.8.2

Edit:

I am not sure what the doc meant when it said extract the sdk dir to /opt/opengeo/suite/sdk, but I copied the /usr/share/opengeo-suite/sdk to the /opt/ dir, by creating opengeo, suite dirs. Is this the right method?


texts: Convert resolution of geojson file from 1:10 to 1:110

texts: After searching a while I finally found a geojson file of the continents. But it is large file (even gzipped) and I do not need the high resolution (1:10). Can I somehow convert it smaller resolution (1:110)?


texts: How to break bond of select and modified feature OpenLayers

texts:
I have select and modify interaction instantiated. Also I have event onChange binded to selected feature.

select.getFeatures().getArray()[0].on('change',function (){
    sources.verteksi.clear()
    koordinate=select.getFeatures().getArray()[0].getGeometry().getCoordinates()
    koordinate.forEach(function(koord){
        koord.forEach(function(ko){
            sources.verteksi.addFeature(new ol.Feature({geometry:new ol.geom.Point(ko,'XY')}))
        })
    })
})

On the other hand I have layer which I modify and select :

arkod_editing=new ol.layer.Vector({
    source: sources.arkod_editing,
    style: new ol.style.Style({
        stroke: new ol.style.Stroke({
            color: 'red',
            width: 2
        }),
        fill: new ol.style.Fill({
            color:'rgb(190,270,80,0.5)'
        }),

        })
  })

The problem is when I modify feature from arkod_editing layer I send it to transactional WFS but before I am setting id of feature which triggers onChange event like feature is coming from select layer. Is it possible to outflank the bond beetween select and modified feature in layer arkod_editing?


texts: Error in function when a class constructor called twice in ArcGIS API for JavaScript?

texts:
I want to create multiple search widgets using the "esri/search module". The second search constructor (for var geosearch) is showing an error that it is not part of function. Is it possible to use a class constructor again in JS?.  Please check the below script:

var map;
require([
"esri/map",
"esri/layers/FeatureLayer",
"esri/dijit/Search",
"dojo/on",
"dojo/dom",
"esri/dijit/editing/Editor",
"esri/toolbars/draw","esri/symbols/SimpleMarkerSymbol",
"esri/symbols/SimpleLineSymbol",
"esri/symbols/SimpleFillSymbol", "dojo/domReady!",
],

function(Map, FeatureLayer, Search,dom, Editor, Draw,SimpleMarkerSymbol,SimpleLineSymbol, SimpleFillSymbol)
 {
 map = new Map("mapDiv", {
 basemap: "streets",
 center: [-56.049, 38.485],
 zoom: 3
 });

 var Search = new Search({
 allPlaceHolder: "Find by address",
 autoNavigate: false,
 autoSelect: false,
 enableHighlight: true,
 enableInfoWindow: true,
 enableSearchingAll: true,
 map:mapDiv,
 showInfoWindowOnSelect: true,
 zoomScale:1000
 },"Address"); });

 var geoSearch = new Search({
 allPlaceHolder: "Find by co-ordinates",
 autoNavigate: false,
 autoSelect: false,
 enableHighlight: true,
 enableInfoWindow: true,

 map:mapDiv,
 showInfoWindowOnSelect: true,
 zoomScale:1000
 },"Coordinates");});


texts: Zoom to a vector feature in a vector layer

texts:
I am using OpenLayers+GeoServer and I have a vector layer in a map with some features
I would like to make a zoom and show the map centered in one of these features that I can get by ID.
How can I do this operation?


texts: GeoServer automate "reload feature type"

texts:
Using GeoServer, if you add a feature type bound to a SQL table that is empty, when you then populate that table no data shows up in layer preview. Clicking the "reload feature type" button corrects the problem. Is there any way to do this through the REST API, or some other way to automate it?

It looks like the resource reset command could achieve this, but it seems to do so for all stores. I would rather only do it for the store in question.


texts: The labels I made on my Contours through the label tab in the properties panel disappear after I save my QGIS project as a shapefile.


texts: Mapbox WMTS layer in QGIS poor print quality

texts:
I've added a Mapbox WMTS layer to my QGIS project.

However when I use the print composer and export my map the result is blurry / not sharp. The image in the print composer itself is sharp, it is only after export that the problems occur. Any idea how to solve this? I've also tried exporting as PDF but the result is the same low quality.

I also tried increasing DPI and total number of pixels but the results are disappointing.


texts: Off panel widget does not appear in web appbuilder panel

texts:
I am using web appbuilder developer edition. And have created an off-panel widget that has a manifest.json like following:
{
  "name": "SomeSearch",
  "platform": "HTML",
  "version": "2.17",
  "wabVersion": "2.17",
  "author": "",
  "description": "This is the widget used in some search.",
  "copyright": "",
  "license": "http://www.apache.org/licenses/LICENSE-2.0",
  "properties": {
    "hasVersionManager": true,
    "inPanel": false,
    "hasConfig": true
  }
}

But it does not appear in side bar red area. Why?


texts: Maintaining number of row and columns in raster after clipping?

texts:
I have a DEM raster (WGS-1984) and MODIS scenes (1200*1200 pixels). I want to clip the DEM raster based on these scenes and after that create a new raster based on MODIS scenes bounds which has value at clipped part and no value at other areas which the DEM can not cover. At the final stage I want to have a raster with 1200*1200 rows and columns. I used clip and extract by mask from arc map (I checked maintain clip extension) but my final output has different number of rows and columns in contrast with MODIS. It is to be noted that MODIS scenes has the following coordinate information:

"Unknown_datum_based_upon_the_custom_spheroid_Sinusoidal
Authority: Custom

Projection: Sinusoidal
false_easting: 0.0
false_northing: 0.0
central_meridian: 0.0
Linear Unit: Meter (1.0)

Geographic Coordinate System: GCS_Unknown_datum_based_upon_the_custom_spheroid
Angular Unit: Degree (0.0174532925199433)
Prime Meridian: Greenwich (0.0)
Datum: D_Not_specified_based_on_custom_spheroid
  Spheroid: Custom_spheroid
    Semimajor Axis: 6371007.181
    Semiminor Axis: 6371007.181
    Inverse Flattening: 0.0"
and DEM has the following coordinate information:"GCS_WGS_1984
WKID: 4326 Authority: EPSG

Angular Unit: Degree (0.0174532925199433)
Prime Meridian: Greenwich (0.0)
Datum: D_WGS_1984
  Spheroid: WGS_1984
    Semimajor Axis: 6378137.0
    Semiminor Axis: 6356752.314245179
    Inverse Flattening: 298.257223563"


texts:
I am quite new to GIS area.
I have a polygon and use concave hull function I could get a boundary like this

What I what to do is add a certain-width band along the boundary, so that it could look like this

Is there any way to do this in R?


texts: How can I add properties in WFS layer using OL

texts:
Is there any way I can write in shapefile. As per my knowledge if I load the WFS layer in QGIS and if I have write permission in that layer then I can write in shapefile. Can I add properties in the shapefile using WFS layer using OpenLayers?
My current code is :
fetch(url, {
      method: 'POST',
      body: new XMLSerializer().serializeToString(featureRequest('AB')),
    })
      .then(function (response) {
        return response.json();
      })
      .then(function (json) {
        console.log(json);
        var features = new ol.format.GeoJSON().readFeatures(json);
        //features.set("property", 123);  #Trying to add like this
        vectorSource.addFeatures(features);
        map.addLayer(vector);
        map.getView().fit(vectorSource.getExtent());
      });

I have done modification in code but its still not working:
var pointGML = new ol.format.GML({
        featureNS: 'http://www.openplans.org/uk_portal',
        featureType: ['AB'],
    });

    var pointWFS = function(p,f) {
        var formatWFS2 = new ol.format.WFS();
        switch(p) {
            case 'insert':
                node = formatWFS2.writeTransaction([f],null,null,pointGML);
                break;
            case 'update':
                node = formatWFS2.writeTransaction(null,[f],null,pointGML);
                break;
            case 'delete':
                node = formatWFS2.writeTransaction(null,null,[f],pointGML);
                break;
        }

        var s = new XMLSerializer();
        var str = s.serializeToString(node);
        $.ajax('http://localhost:8080/geoserver/wfs',{
                type: 'POST',
                dataType: 'xml',
                processData: false,
                contentType: 'text/xml',
                data: str,
                success: function(data) {
                    console.log(data);
                }, error: function(data) {
                    console.error(data);
                }
            }).done();

        console.log(" point features were posted to server");
    }

    var interaction = new ol.interaction.Draw({
        type: 'Point',
        source: vector.getSource()
    });
    map.addInteraction(interaction);
    interaction.on('drawend', function(e) {
        var fea = e.feature;
        pointWFS('insert',e.feature.set('ferrets', 'bar'));
    });

I am getting this response:


texts: Drawing a buffer from a POINT and a value for Acreage

texts:
We have a situation where we get a single lat/long point and a value for "acreage", and we need to draw a polygon circle that is roughly the size that the acreage value implies. What we have is a query that comes extremely close, but the total area is still off, for some reason.
create TEMP TABLE g (geom geometry, acres float);
 insert into g VALUES (ST_SETSRID(ST_POINT(-74,34.00), 4326), 100);
 insert into g VALUES (ST_SETSRID(ST_POINT(-74,34.00), 4326), 10000);
 insert into g VALUES (ST_SETSRID(ST_POINT(-74,34.00), 4326), 100000);
create temp table gg  AS select
    ST_TRANSFORM(
        ST_BUFFER(
        -- Build a buffer, calculating radius as sqrt(area/pi)
            geom::geography,
            SQRT(
                 (acres/.00024711)/3.1415
             )
        )::geometry
    , 4326)::geography as b,
    acres,
    acres/.00024711 as area_in_meters
FROM g;

SELECT
    b,
    st_area(b, TRUE), -- This and area_in_meters are not _quite_ equal
    acres,
    area_in_meters
FROM gg;

I believe its one of a few factors, but I am having a hard time deriving the right test.
acres/.00024711 as area_in_meters - is this an accurate way of converting acres to sq meters?
ST_TRANSFORM(buffer_geom, 4326) Will this always yield an accurate geometry for the buffer we want? are we loosing data jumping around between geography and geometry?
st_area(b, TRUE) b is a geography, so this is in meters right? its what the docs say, but perhaps ::geography as b isn't fully performing the cast I think it is?


texts: Unknown coordinates in Shapefile, projecting doesn't give expected results

texts:
I'm reading a Shapefile containing polygons representing certain areas in my neighborhood. For this I use C# and a library called DotSpatial.

I can iterate the so called Features in the Shapefile and each feature has Coordinates with a X and Y property value. They look like this:

[124060,375000004, 487043,270000004]
[124068,295000004, 486998,390000004]
[124071,430000004, 486891,470000004]
[124010,435000004, 486694,295000004]

I need latitude/longitude values so after some reading I found that I needed to reproject the coordinates, either with/without projecting file (see answer: here).

My goal is to get values in the area of Amsterdam (which is what the Shapefile contains): [52.4,4.9]

This is possible by calling (if projection file is included):

  file.Reproject(DotSpatial.Projections.KnownCoordinateSystems.Geographic.World.WGS1984);

Resulting in:

[-37,9115551003779, 23,9336542782764]
[-37,9117191600094, 23,9331853015717]
[-37,9123377534534, 23,9321317059377]
[-37,9144451002488, 23,9304589279205]
[-37,9153029784154, 23,9298198785957]

This looks promising, but nowhere near [52.4,4.9]

Or (without projection file):

  file.Projection = DotSpatial.Projections.KnownCoordinateSystems.Projected.UtmWgs1984.WGS1984UTMZone32N;            file.Reproject(DotSpatial.Projections.KnownCoordinateSystems.Geographic.World.WGS1984);

Resulting in:

[5,61357508612624, 4,39863081202557]
[5,61364817701951, 4,39822583083791]
[5,61368074834023, 4,39726037781928]
[5,61314000884886, 4,39547721608089]
[5,61290560849057, 4,39478423443173]

But both are not the results what I was hoping for. I am in no way any geometry expert and I understand about spheres and cylinders etc. but I'm really unable to figure out how to get my coordinates in lat/long format. Is this possible?

The .prj file looks like this:

PROJCS["RD_New",GEOGCS["GCS_Amersfoort",DATUM["D_Amersfoort",SPHEROID["Bessel_1841",6377397.155,299.1528128]],PRIMEM["Greenwich",0],UNIT["Degree",0.0174532925199432955]],PROJECTION["Double_Stereographic"],PARAMETER["False_Easting",155000],PARAMETER["False_Northing",463000],PARAMETER["Central_Meridian",5.38763888888889],PARAMETER["Scale_Factor",0.9999079],PARAMETER["Latitude_Of_Origin",52.15616055555555],UNIT["Meter",1]]


texts:
I am looking for an alternative tool (proprietary or FOSS) that matches the  functionality of Temporal Analyst for ArcGIS by  DHI Software. The best I could get was the QGIS plugin  Time Manager  whose functionality I perceive is more for visualisation than analytical (similar to esri's Time Slider.)

[Data Analysis - the Temporal Analyst for ArcGIS offers a suite of tools for data analysis and processing, incl. basic statistics and arithmetics, time step resampling, gapfiling, distribution and CDF plots, duration curves, double mass and double scatter.]

I am not a statistician nor hardcore developer so the more the "point and click" I can get the better.


texts: Converting uint16 RGB&NIR geotiff image to a uint8 RGB geotiff?

texts:
I am trying to convert a uint16 geotiff image with 4 bands (RGB&NIR) with LZW compression to a uint8 geotiff image with 3 bands (RGB) without compression. Here is what I have tried:

gdal_translate -ot Byte -b 1 -b 2 -b 3 -r cubic "20170729_182758_0f28_3B_AnalyticMS_SR.tif" "output.tif"

When I drag and drop the image in QGIS, I get meaningful colours. When I generate a tif file using the above command, the colours are off (looks washed out). Any ideas what might be the problem?


texts:
I am trying to install the python library numba with pip from the OSGeo4W shell, but I haven't been able to do it. The only way I could do this was through the anaconda installer. but since Anaconda and OSGeo4W have their own versions of python, this solution is not viable to me.

I need to use this library in a qgis script, Is there a way to do this ?


texts: QGIS: alternative to add/edit virtual layer. With python?

texts:
Following a previous question (Select polygons with other polygons layer?), I select some polygons using the query in the answer to create a virtual layer:

SELECT greens.*  FROM greens, purples
WHERE st_area(st_intersection(greens.geometry, purples.geometry))
      > 0.5 * st_area(greens.geometry)

It works really well, but my layers contain 1.800.000 (and sometimes more) polygons. So that QGIS freezes and virtual layers only appears few days later. When I am trying to record it, it takes again several days to do that. Since I have to do the operation on several layers, it takes an incredible amount of time (that I do not have). Does someone know an alternative to use that kind of query in an other way to select the right polygons?


texts: How to batch merge point, polygon and line feature classes into a new feature class?

texts:
Using ArcGIS 10.3 for desktop.

For our noxious weed surveys, we have some features recorded as polygons, some as lines, and some as points.  I would like to convert everything to points, and then merge it all into one feature class for each species.  I have minimal python experience, but am willing to try.

It is currently set up so that there are three file geodatabases.  NoxiousWeeds_points.gdb, NoxiousWeeds_poly.gdb, and NoxiousWeeds_lines.gdb.  Within each geodatabase is a separate feature class for each species.  E.g., "HimalayanBlackberry_points" or "ReedCanarygrass_poly," etc.  Names are consistent between the three geodatabases, the only thing that varies is the suffix.

What I'm Trying to Do:

I would like to be able to simply point arc to the three geodatabases and have it automatically convert the line and polygon feature classes to points, then merge with the matching point feature class based on the name, and then saved to a fourth, separate geodatabase with a slightly modified name.  For example, "Knotweed_poly" and "knotweed_line" would both be converted to points, merged with the existing "knotweed_points," and then saved as "Knotweed_merged" in a fourth, separate geodatabase named "NoxiousWeeds_merged."  Ideally I could simply point arc to the geodatabases that contain all these files and it would go through this whole process for each species by itself without me having to manually select each file.  Note that I would also like the original source files to remain unaltered.

I suspect this is possible through some combination of python and ModelBuilder.  If so, how would you recommend I go about it?


texts: Using real time data (traffic) for network analyst in ArcGIS for Server?

texts:
I have a problem with using traffic data for routing. I'll describe the scenario:

There is a local traffic service (real time) that has the speed of each segment in the network. This information is converted to a graphic and displayed on a map. I want to use the information in network analyst. For example the user wants to route between point x and point y using one-way, time and traffic restrictions. I want to use these graphics in network analyst such that the results of routing not include the segments that have high traffic.

How can I combine these real time data in network analyst services?


texts: What QGIS line styles work best with Garmin maps?

texts:
I am using QGIS to create some maps to load on to a GPS using the plugin.  When I display the maps on the GPS the lines are very faint and I can not distinguish colours and some colours vanish altogether. But even the black lines don't stand out.

I have tried increasing the thickness of the lines in QGIS but this does not seem to make any difference.

I am displaying the maps on an ETrex 30.


texts: Capitalizing and changing font/color in Label Expression using Python Parser of ArcGIS Pro?

texts:
Currently Using ArcGIS Pro 2.3.2 and Python.

I am currently trying to label Road names, which are all uppercase in the database. I'm unable to edit the database but would like to Capitalize the road names and pick out one street to change the font and color for the road name.

I'm new to python and am having trouble working out how to combine these two expressions.

def FindLabel ( [road_name1] ):
  S = [road_name1]
  S = S.title()
  return S

def FindLabel ( [road_name1] ):
  if str([road_name1]) == 'LOOKOUT ROAD':
    return "<CLR red='0' blue='0' green='0'><FNT size = '22'>" + [road_name1] + "</FNT></CLR>"
  else:
    return "<CLR red='128' blue='128' green='128'><FNT size = '14'>" + [road_name1] + "</FNT></CLR>"


texts: r - Setting custom projection in PROJ >=5 paradigm

texts:
Before PROJ >=5.0, one could set a custom projection centered anywhere in the world by adjusting the +lon_0=and +lat_0= arguments in proj4string specification.

For example, a projection centered on Europe.

proj <- CRS("+proj=laea +lat_0=53 +lon_0=9 +x_0=4321000 +y_0=3210000 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs")

As I understand it, the use of proj4string is being phased out, and this sort of use is discouraged.

So, with the new projection specification methods in R (e.g. WKT2) is there way to set a custom-centered projection? And if not, is anyone aware of a different automated way of selecting a roughly appropriate projection based on a set of X,Y coordinates (e.g. the centroid)?


texts: Adding layers in QGIS with predefined order in TOC and also make some of them to have visibility turned off when added

texts:
I am working on a project and I have a lot of shapefiles with the same name but with different data and when I add them to QGIS they are sorted randomly in TOC (Layers Panel).

Can I configure their order? So when I add them, they are sorted based on some rule or based on type of the layer?

And can I configure that some of them are hidden(visibility turned off) when added?


texts: Why are there multiple symbol icons for the same feature type for raster field in ArcGIS?

texts: Using ArcGIS 10.7 desktop. I manually entered general landuse types for a field I created. For some reason only some of the feature types are consolidating in the symbology while others are creating duplicates. The field names are exact same for each type, copy/pasted, so its not due to additional spaces, etc.


texts: Using CSW for GeoServer PostGIS database

texts:
I am quite new to GeoServer and general OGC fundamentals. I have GeoServer set up and I set up a PostGIS data store which links to a database

I have already been able to make successive calls using WFS to extract data from here. I am looking at other service capabilities offered, like CSW.

I cannot really grasp how it works though. Can someone tell me how I can get records from my database? (using HTTP GET requests)

I can do a GET that returns all the layers I have set up:

http://{host}:{port}/geoserver/csw?service=csw&version=2.0.2&request=GetDomain&propertyName=Title

But using the GetRecords or GetRecordById returns nothing to do with my data in the database though..

Can someone demonstrate a simple REST call to retrieve a record from the PostGIS database?


texts: How to calibrate/put points on a network

texts:
I would like to staled points on a network. I have shape-point and I would like to take them and put/position  them to position correctly on the network.

I think it's easy but I'm lost


texts: 1,000 Feature Limit in AGOL - Find Nearest Tool

texts:
I am trying to run the Find Nearest tool in ArcGIS Online. The tool's documentation does note that it will not run on layers with more than 1,000 features, but there are several Esri blog/Q&A's out there that have addressed this issue to offer work arounds. The solution was to publish the layer(s) as a hosted feature service either through ArcMap or from within Portal. I have done this and I still get the "Cannot process more than 1,000 records" error.

Does anyone know a proper workaround for this? I understand why they limit to 1,000 (because of drawing/rendering), but for a lot of cases, this would be exceeded by many users.

I do not have the Network Analyst extension for my desktop GIS to run this tool there, so I'm looking for alternative ways to compute a closest facility analysis.

My origin layer has 900 records, and my incidents layer has 11,000.


texts: Multiple reduceRegions for time series in GEE js

texts:
I want to gather time series data (mean band values) from MODIS for multiple regions, but I want to filter the data so that I exclude images with low numbers of pixels in my regions.
Below is the code to use reduceRegions to get the pixel count and mean values, but how do I get exclude the images with low pixel counts?
var modis = ee.ImageCollection("MODIS/006/MOD09GA")
  .randomColumn()
  .filter(ee.Filter.lt('random', 0.1))
  .limit(10)

var fc = ee.FeatureCollection([
    ee.Feature(
        ee.Geometry.Polygon(
            [[[-142.36414308128087, 60.44668098352718],
              [-142.36414308128087, 60.171875701856855],
              [-141.62256593284337, 60.171875701856855],
              [-141.62256593284337, 60.44668098352718]]], null, false),
        {"site": "a"}),
    ee.Feature(
        ee.Geometry.Polygon(
            [[[-150.47861099878585, 60.14760911377782],
              [-150.47861099878585, 59.90886167691271],
              [-149.84140396753585, 59.90886167691271],
              [-149.84140396753585, 60.14760911377782]]], null, false),
        {"site": "b"})
]);

// Mask some pixels
var maskModis = function(image){
  image = ee.Image(image)
  var ndsi = image
    .normalizedDifference(['sur_refl_b04', 'sur_refl_b06'])
  var ndsiMask = ndsi.gt(0.8)

  return image.updateMask(ndsiMask)
}

// Apply the mask
var modis = modis.map(maskModis)
// var testImg = ee.Image(modis.first())
// Map.addLayer(testImg)

// 1. Count number of pixels in each feature.
// Return FeatureCollection with date, site, and pixelCount.

var getPixelCount = function(image){
  image = ee.Image(image)
  var count = image
    .select(1)
    .reduceRegions({
      collection: fc,
      reducer: ee.Reducer.count(),
      scale: 500
    })
  return count
}

var counts = modis
  .map(getPixelCount)
  .flatten()
print('counts', counts)

// 2. Get mean band values for each feature, each date
// Same appriach as (1), map reduceRegions

var getMeanVals = function(image){
  image = ee.Image(image)
  var date = ee.Date(image.get('system:time_start')).format('YYYY-MM-dd')
  var mean = image
    .select('sur_refl_b01', 'sur_refl_b04')
    .reduceRegions({
      collection: fc,
      reducer: ee.Reducer.mean(),
      scale: 500
    })
  return mean //returns featureCollection
    // add date to each elt in collection
    .map(function(f){
      return f.set({'date': date})
    })
}

var means = modis
  .map(getMeanVals)
  .flatten()
print('means', means)

This gives me two featureCollections, but I only care about the 'means' collection. How do I filter the 'means' collection by the 'counts' collection?
https://code.earthengine.google.com/?scriptPath=users%2Fcaseyengstrom%2Freprex%3AmultipleReduceRegions


texts: Line Fill Symbol is not shown on legend in ArcGIS Server print service

texts:
I choose Line fill symbol for my polygon layer. When I export map in ArcMap it is ok. But printing with print service this symbol is not shown in the legend.  I add a screenshot.

How can I solve this problem?

I use:

ARCGIS SERVER 10.5.1 (Linux)

Template is desinged in ArcGIS Desktop 10.5.1


texts:
I am still bothered by this and Google has not released another version yet.  Any ideas, please?

Today I upgraded from Google Earth Pro 7.1.7.2602 to the recently released 7.1.8.3036 and I immediately noticed a problem.  The font used by default in all placemark balloons changed from sans-serif to a small serif one.  This affects not only all of my own placemarks but those from the standard Layers.  For my eyes the old font was already marginal and the new one is simply unworkable.  I know I can change the <BalloonStyle> for my own placemarks but for several reasons I don't want to do that, and it still would not fix the font in the standard Layers.

How can I control this font?  It does not originate in myplaces.kml as that file is unchanged.

7.1.7.2602 balloon font:

7.1.8.3036 balloon font:


texts: I have installed QGIS 2.14.1,  I am looking for the OpenLayers Plugin which I could not find under Install Plugins, as well as Web menu.


texts: Union Tool creating slivers in ArcGIS Desktop?

texts:
I am using ArcGIS Desktop 10.4.

I have 2 feature classes in the same GDB. One feature class is a subset of the other one (i.e. it's a direct copy, with a subset of features retained, the others deleted). All 190 features in both FCs are identical.

My assumption then is that a union between these two FCs should yield the same total number of features as the source FCs (i.e. 190) because their geometries are directly coincident. However, it doesn't, and a number of slivers are created in the union shapefile.

Can someone please explain why this happens?


texts: What's the difference between gdal_cachemax and setting the cache size with -wm in gdalwarp?

texts: I'm trying to get better performance out of gdalwarp, and am wondering if there is any difference between setting gdal_cachemax and using the -wm option in gdalwarp?


texts: Python Add In Change tab in Table of Contents

texts:
I'm working on an add in that generates new structures for an addressing system. The tool is working well, but as it runs, at some point it changes tabs in the Table of Contents to the "List By Source" tab. I'm guessing it does this either when I add temporary tables or shapefiles.

The fact that it changes tabs is not a problem, but I would like for the tool to leave the user in the TOC tab "List By Drawing Order" when it is finished.

Is there any way for a python code to change tabs in the TOC ?


texts: Reconditioning and Drainage Line ArcHydro

texts: I cannot get ArcHydro to run Reconditioning - says object reference not set to an instance of an object.  So, not finding what that refers to, I skipped that and used the sequence of Terrain Preprocessing functions up to Drainage Line Delineation, but that won't run either, giving a message that a grid code is dropped. Does someone know how to get past these roadblocks?


texts: How can I get a map image from an ArcGIS mapserver centered on a single feature?

texts:
I have a Python program that is not using arcpy. I need to get an image from a mapservice centered on a single feature using the REST api.

I can call the export method to get a map image:

http://<server>/arcgis/rest/services/<name>/MapServer/export?bbox=<bounds>&f=image

I can also query the mapservice layer using a where clause and get the feature geometry (in my case, it's going to be a polygon feature).

Without using arcpy, is there an easy way to take the polygon's geometry, get a rectangular extent to generate the bounding box necessary to call the export?

What would be ideal is if I could just say "Give me a map image centered on a feature from a layer defined by this where clause" but, unless I'm missing something, you can't do this in one step. Neither can I figure out how to do a Query on a layer and get the extent back instead of just the feature geometry.


texts: Trace upstream/downstream for multiple pairs of points in ArcMap

texts:
I have two point layers, one is set of origin points and one is a set of destination points, each with 608 points. Object IDs match the origins to the destinations. The points are all snapped to the National Hydrography Dataset (NHD) streams. I am working in ArcMap 10.5.
I need to find the distance between the paired points along the NHD stream network, AND state whether the route is upstream or downstream.
Edit 1: Image added by request. In this example, the destination point is upstream of the origin point along the route I highlighted. There are several other destination points which are closer. I need the distance and direction (upstream/downstream) to the matching destination point only. There are 608 of these pairs across the study area.

I have successfully found the distance between origins and destinations along the NHD network, as shown in the image, by generating an OD Cost Matrix in Network Analyst, adding origins and destinations, and selecting the lines where origin ID = destination ID. However, this does not include whether the destination is upstream or downstream.
One potential option would be to find a way to trace upstream of each origin point until reaching the correct destination point. If it does not find the correct destination point upstream, then it could be assumed that the destination point is downstream (or vice-versa). Tracing a geometric network have a way to do this, but I cannot find a way to implement this to find the paired destination point, rather than just finding the closest destination point.
Any ideas for finding the distance and direction along a stream network between pairs of points?


texts: Coordinate systems between QGIS and ArcMAp

texts:
My country (Slovenia) is switching to a newer coordinate system and I've received some sample data in this new coordinate system. I went ahead and searched the definition of the CRS on this website. The CRS is EPSG:3794.

SAMPLE DATA TO REPLICATE THE PROBLEM:

You can access the sample data here
The url for the WMTS is: http://prostor4.gov.si/ows2-gwc-pub/service/wmts?

ARCGIS PROCESS:

Downloaded the option for ArcGIS -> ESRI WKT and imported the downloaded .prj CRS to ArcMap (version 10.2.1 or 10.3).
Restarted the application and reopened a fresh new project.
Firstly I set the newly imported CRS for the project by right clicking the "Add data".
After that I added the provided sample data (you can find it in the link)
Added a new WMTS connection with the service provided by my country in the newest coordinate system, the layer I chose is named "SI.GURS.ZPDZ:DOF050".

I then noticed that the 2 layers do not align.

Here is the picture:

Now this got me abit confused, because that would mean my coordinate systems aren't matching. So to confirm my "theory" I repeated the process in QGIS.

QGIS 2.18 LTS PROCESS:

Copied the .proj4 definition and created my new CRS.
Again added my sample data.
Added a new WMTS connection and inserted the layer: "SI.GURS.ZPDZ:DOF050"

I zoomed in on the data layer  and to add even more confusion the data aligned "perfectly":

Now this made me suspect that ArcMap is behaving incorrectly, but I can't find out why or how.


texts:
I'm trying to import downloaded SVG-Files into QGIS to use them as layer symbols, but nothing works. This question was asked before, but unfortunately, the suggested solutions don't work.

When I do

  "Settings->Options->System->SVG Path"

And open the folder, there are no files in it. The only way I can access those is by opening the folder without any QGIS involved, click on the files and open them in a browser, then it shows me the symbol, but after hours of trying I didn't find any possibility to import them into QGIS.

I only can open them in the browser-window, but not add them as a layer. When right-clicking on properties it says under the headline "Meta Data": Error.

It would be very nice if someone could help me. I may try to figure out the problem more extensive.


texts:
Whenever I run a geoprocess overnight, it seems my computer reboots. Does ArcGIS save some log file information? Where does it save it?

For example, when I run a large dataset or image process I need to perform overnight, does it save the logfile that tells you what it performs from start to finish?


texts: Displaying UTFGrids with Mapnik and TileStache

texts:
I'm new to using Mapnik and TileStache. The reason i wanted to use these two pieces of software was to understand the UTFGrid specification i've seen used all over the place ( examples here, here, here )

I copied the OpenLayers examples linked to above. And i've created my own here.

I have a couple questions about Mapnik and TileStache that maybe someone here can help me resolve:

1) Does anyone know how i can display the UTF Character overlays by themselves? The endgame here is twofold -- a) i want to make a cool ascii-art-like map and b) i want to understand more about how this technology works. I assume the utfgrid overlays are drawing as vector tiles beneath the image tiles but firebug isn't showing me much.

2) It seems that the utfgrid will exactly match the the geometry it is created from in the mapnik xml file -- meaning a direct overlay. If i render custom marker .png(s) for points is there anyway to make the utfgrids overlay for the whole marker instead of just the point? Is there a way to use utfgrids with points or is it meaningless ( pointless even ;-) )?

Thanks
Mazery


texts:
I have a PostGIS 1.5 database and want to update to 2.0. I read that because this is major release I have to follow the 'Hard update' procedure. However the spatial data I have is only test data so I can drop it.

If this is the case, is there easier procedure to update to 2.0 if I don't care about spatial data (I do care about non-spatial tables)?


texts: Using Area Solar Radiation tool for many rasters with different latitudes and Z values?

texts:
I want to count direct radiation for all Europe. As the Area Solar Radiation (ASR) tool allows to work only on small landscapes, I cut my raster-Europe into fragments with respect to latitude. So each piece of raster will be analyzed with a different values of Latitude and Z factor. Everything is working fine, but...

I have over 4000 these fragments of raster and doing it manually does not make sense...

So I used both batch mode and ModelBuilder, but in both cases, both Latitude and Z factor are set to default values, ie, for Lat=45 and for z = 1. If I use the ASR individually, these values are automatically calculated.

What should I do to batch mode or ModelBuilder calculated the value of Latitude and z-factor for each piece raster automatically?


texts: Options for PostGIS and Mapserver outputing SVG, and client with OpenLayers

texts:
This question is related to "Options for displaying PostGIS vectors in OpenLayers?", but I need SVG output on server. OpenLayers must render SVG directly.

We are using Mapserver6, and thinking on WFS server mode (can be other?)... Need TinyOWS to faster output or easier configuration?

There are tutorials (?) for the problem of multiples minscale/maxscale on Mapfile layers -- for use something like st_union and st_simplify PostGIS functions to reduce trafic and complexity on big vectors.

For display we think to use some like,

map.addLayer( new OpenLayers.Layer.Vector(
    "Setores",
    {   strategies: [new OpenLayers.Strategy.BBOX()],
        protocol: new OpenLayers.Protocol.WFS({
            url: conf_url_wfs,
            featureType: "setores_vector",
            featureNS: "http://www.openplans.org/topp"
            }),
        renderers: ["SVG"]
},
    {isBaseLayer:false,visibility: true}
));

There are a Mapserver example like this simple GeoServer example?
Examples with Mapfile also (and connectiontype=postgis).


texts: Change Axis Order when transferring data from WFS to PostGIS DB

texts:
I am using GeoKettle to transfer data from a WFS service to a PostGIS database.

The projection of the geometries is EPSG:4326 :

GEOGCS["WGS 84",
  DATUM["World Geodetic System 1984",
    SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],
    AUTHORITY["EPSG","6326"]],
  PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
  UNIT["degree", 0.017453292519943295],
  AXIS["Longitude", EAST],
  AXIS["Latitude", NORTH],
  AUTHORITY["EPSG","4326"]]

In the data I have latitude/longitude tuples. When it's imported into PostGIS, it appears as longitude/latitude.

I want to specify a CRS WKT parameter so that when data is imported into PostGIS, the geometries will appear where they should.

Is this even possible with an CRS conversion?


texts: Create attribute field on a geopackage with precision using QGIS 3.4

texts:
On a QGIS geopckage when I create a new field for decimal numbers, I can only create it as a decimal (real). So it populates the field with many decimals as shown. I don't mind having more decimal places for the fields. But when I make a map with attributes on print composer, I would like to show only three decimal values.
I don't want to use refactor fields or other processing tools at later stage because I import a lot files so that I can avoid extra steps to link it to the map.


texts: QGIS 'Join attributes by location' algorithm (sometimes?) not working properly in graphical modeler

texts:
In QGIS 2.14.10 I've set up a simple model to assign min and max of attributes of a point layer (adresses, about 2500 features) to a polygon layer (parcels, about 14000 features) to polygons containing the points (hint: polygons were created with polygonize algorithm, so these should have valid topology):

The 'Join attributes by location' algorithm is set up like this:

The resulting bla polygon layer is created, but has zero features:

After some playing around, I found that turning of on-the-fly-projection solved my issue ...

... but unfortunatly not for long. So toggeling on-the-fly-projection sometimes helps, but the behaviour of 'Join attributes by location' algorithm does not seem to depend on whether OTF projection is ON or OFF and is, furthermore, unpredictable when toggling OTF projection.

Am I missing something in modeler, or is there any advice of using this algorithm in the context of OTF projection?


texts: Selecting all objects with specific value from GeoJSON in new variable

texts:
So I have a GeoJSON that I want to pull/filer/select all kitespots with a wind direction that includes "S" in Leaflet. (..and yes, am a beginner at coding..)
var kitespots = {
  "type": "FeatureCollection",
  "features":[
    {
     "type": "Feature",
     "properties": {"name": "STUBBSAND", "windDirection": "S/SW"},
     "geometry": {
       "type": "Point",
       "coordinates": [18.959655761718746, 63.220636696106915]}
    },
    {
     "type": "Feature",
     "properties": {"name": "GULLVIK", "windDirection": "SW/S/SE"},
     "geometry": {
       "type": "Point",
       "coordinates": [ 18.82786273956299, 63.214932082292805]}
    },
    {
     "type": "Feature",
     "properties": {"name": "VANNVIKEN", "windDirection": "SE/S"},
     "geometry": {
       "type": "Point",
       "coordinates": [18.815631866455078, 63.20794965080512]}
    },
    {
     "type": "Feature",
     "properties": {"name": "STORSANDEN", "windDirection": "N/NW"},
     "geometry": {
       "type": "Point",
       "coordinates": [18.528785705566406, 62.99013053076868]}
    }
  ]
}

So far I have tried;
// - Select all kitespots facing south - //
var southSpotCluster = new L.markerClusterGroup();
var windDirectionS = L.geoJson(kitespots,{windDirection: "S"});
southSpotCluster.addLayer(windDirectionS);

// - ...and.. - //

var southSpotCluster = L.geoJson(kitespots)
var windDirection = ["S"]
for(var windDirection in southSpotCluster){
    windDirection.push(southSpotCluster["S"]);
}
southSpotCluster.addLayer(windDirection);

var overLays = {
    'All Kitespots': allSpotsCluster,
    'Kitespots S': southSpotCluster
};

L.control.layers(baseLayers, overLays).addTo(map);


texts:
I'm a new QGIS user and am having trouble adding tif files from OS VectorMap
District.
I get them on to the map but they always appear in the bottom left hand corner, which I would guess indicates a CRS discrepancy. However, I have adjusted this so all layers are using OSGB 1936/BNG EPSG:27700, but nothing happens. Where am I going wrong?


texts: Enabling (python) plugins on QGIS3 (Linux)?

texts:
I'm running qgis 3.0.3 in Fedora 28 and I seem to have several, related problems with plugins.

Plugins --> Manage and Install Plugins: the window that opens shows 7 pre-installed plugins (I think there used to be 8 before, how come?).

Under the Settings tab, it says that "no python support detected, thus no settings available"

What dependency am I missing on my system?

Also, the search from the search field doesn't seem to work, it searches only through the 7 available plugins.

I downloaded a plugin from qgis and extracted it to

/home/myyser/.local/share/QGIS/QGIS3/profiles/default/python/plugins but it's not being loaded.

Any hints?


texts: Multiple Conditionals String with ArcGIS Field Calculator?

texts:
I am trying to use Field Calculator to use a multiple conditional based on a string in another field of the same attribute table.

I have written the code below, but without success.

Pre-Logic Script Code:

def hazardValue( !Group! )
    if !Group! = '1.1':
        return 1
    if !Group! = '1.2':
        return 2
    if !Group! = '1.5':
        return 4
    if !Group! = '2.1':
        return 3
    else
        return 5

Hazard =

hazardValue( !Group! )

It is always returning an error:

  There was a failure during processing, check the Geoprocessing Results window for details


texts:
I am working on a mapping project which includes BCE dates, but I am having trouble getting CartoDB to understand them as dates. So far, the methods I have tried either give me a date in '69 or no date at all.

Any suggestions?