texts: How to link Number value with Field calculator in QGIS Modeler ?

texts: I want to link 2 number values (RMSE & dosemax) in a field calculator formula. When I run the modeler, I add these values but it doesn't recognize them. I dont want to create a new table. Just use these number in the formula.


texts: How to export Extent as shapefile from multiple shapefiles using OGR?

texts: I want to export EXTENT as Rectangular shp file from multiple shape using OGR2OGR?


texts: Converting trajectory to outer edges polygon

texts:
Using PostGIS, I'm trying to solve the following problem:

I have trajectory data (coordinates from API query), like on this image:

and want to convert it into polygon, having only outer edges:


texts: Transforming coordinate that lies outside a grid transform (e.g. NADCON)?

texts:
As I understand it, a datum transform (e.g. urn:ogc:def:coordinateOperation:EPSG::1241) based on a grid interpolation method, e.g. NADCON, is undefined outside of the grid of points. Therefore the lat/long of Mexico city as expressed in NAD 1927 cannot be converted to WGS 1984 by way of this NADCON conversion -- it's not in the United States.

If you call the ESRI projection engine directly (pe.dll#pe_geog1_to_geog2), it simply doesn't transform the coordinates outside of the NADCON grid, though does return a code reminding you that this has happened. (The return code is the # of points transformed, which if less than the # of points you passed in is a sign of trouble.)

My guess is that (most?) everywhere else in the ESRI stack, this return code is ignored, and the transform is simply left as effectively zero outside of the grid area.

This might be as good an answer as any. What should we be doing instead? In our code at the moment we throw an exception, but this is a contagious and awful thing as it means people looking at studies using this transform can't zoom out past a certain extent.


texts: How to define season by month and average data

texts:
I have 10 year NDVI data i.e. two values per month. I would like to average values within a season e.g Season A (January, February and March) per year and plot a graph. y axis=average ndvi value, x axis=year

I have already given a time vector as an attribute to the stack using teh code below

r<-setZ(r,timeNDVI)

I've also calculated mean NDVI for each raster and converted the array to a data frame using the code below

r1 <- as.data.frame(cellStats(r,mean, na.rm=TRUE))


texts: Splitting lines with polygons that do not intersect using ArcGIS Desktop?

texts:
Is it possible to split polylines based on the boundary of the nearest polygon using ArcGIS Desktop?

I am trying to split a multi-part sidewalk layer based on adjacent parcel boundaries and assign each resulting split sidewalk the address of that parcel.  The problem is that most of the sidewalk does not fall within a parcel boundary.


texts: upload data to geoserver with gxp upload panel

texts:
I installed geoserver in my computer. I use wampserver (php, apache), and when i use the   examble "layer-upload.html" then i have the erro 405:

Failed to load resource: the server responded with a status of 405 (Method Not Allowed)


texts:
So I have a set of 5 Fuzzy Memberships and a polygon feature class, to which I already added a yes/no field. For the final map, I need to combine/overlay all Fuzzy Memberships and the rasterized polygon feature class in order to map suitability, and don't know how to combine it.
Any advice?


texts: Making mouse cursor/pointer of ArcMap HTML Popup tool larger?

texts:
Next week I'm giving a presentation, using the ArcMap application of ArcGIS 10.3.1 for Desktop , where to make things easier for the audience to see what I am doing, I would like to increase the mouse pointer size, in this case for the HTML Popup tool.

I have already gone to ArcMap's Customize | Customize Mode and its Options tab to make sure that Large Icons is checked on, but this seems to only affect icon sizes on toolbars.

I've also gone to Windows 7's Control Panel\All Control Panel Items\Ease of Access Center\Make the mouse easier to use to make sure that I am using Extra Large Inverting.

Is there somewhere else that I can go to make the HTML Popup tool (or all tools) use a mouse pointer icon that is about twice as high and wide as it is now?


texts: Calculate coordinate ( lat / long ) with another coordinate and distance

texts:
I already have a formula that can calculate the distance between 2 geographic coordinates:

public function getGeoDistance(Coordinate $geoCoord1,Coordinate $geoCoord2){
    if(is_object($geoCoord1) && is_object($geoCoord2)){
        $aR = $geoCoord1->getLatitude()*(M_PI/180);
        $bR = $geoCoord2->getLatitude()*(M_PI/180);

        $latD = ($geoCoord2->getLatitude()-$geoCoord1->getLatitude())*(M_PI/180);
        $longD = ($geoCoord2->getLongitude()-$geoCoord1->getLongitude())*(M_PI/180);

        $a = pow(sin($latD/2), 2)+cos($aR)*cos($bR)*pow(sin($longD/2), 2);
        $c = 2*atan2(sqrt($a), sqrt(1-$a));

        return self::getWorldRadius()*$c;
    }

    return -1;
}

Now, i want to create another function that can calculate coordinates from another coordinate pair and distance value. The prototype is:

public function getGeoByGeoAndDistance(Coordinate $geoCoord1,int $distance):Coordinate

Does someone know the math beyond it?

P.S. the code is in PHP.

thanks.

Shlomi.


texts: Finding polygons crossed by line using OGR?

texts:
I'm trying to find all the polygons crossed by a single line (a GPS track).
I'm using the OGR library (from python) for computing this, but it's currently a bit 'brute-force' (and slow). For every point of my track, I call the intersect method with all the polygons. The obvious optimization is to check only with adjacent polygons. But I guess this is a classical problem, with an already known solution (which I can't find...).

I would like to avoid using a dedicated database as I'm trying to write a standalone software (spatialite is an option if the DB is the way to go).

(FYI, the current source code is available here: https://github.com/dkm/airspace-checker )


texts: Directed Shortest Path calculation using Dijkstra ignores one-way streets in pgrouting

texts:
I think I found an interesting bug in pgrouting. When I use the Dijkstra algorithm (no matter if I use the shortest_path() or dijkstra_sp_delta_directed() method) to calculate the directed shortest paths between two particular vertexes (id 72846 and id 85556), it only delivers the correct result for one direction.

I experimented a little bit and then I found out that there are two arcs in the roads_po table with opposite direction which connect two vertexes that are part of both calculated shortest paths (gid 700 and gid 169244).

When I calculate the shortest path from vertex 72846 to vertex 85556, everything works fine as can be seen here:

correct_result:

But when I calculate the shortest path from vertex 85556 to vertex 72846 (the opposite direction), the wrong arc with gid 169244 instead of the correct arc with gid 700 is chosen resulting in a violation of the one-way direction.

wrong_result:

It seems that pgrouting calculates both paths correct, but it has a problem that two vertexes along the paths are connected by two arcs, and chooses the wrong arc for output for one of the paths.

Did someone else had such problems with pgrouting?


texts: Different typesets for labeling within one project

texts:
I'm searching a way to show labels with a Western typeset, as well as a different one.

In my case that's Western (English/German/...) PLUS Russian, which I don't speak at all.  But I'd like to be able to display the native labels - for example names of cities.

My idea was to use a second layer and join both, but I can't find a way :-(

Not that I can't find a way to join, but have Cyrillic displayed :-/

QGIS doesn't seam to save/recall if the typeset is changed.
Layer->Settings->Labelling->Font->Typeset

Please "suggest" a menu path - I'm not a native speaker of either, nor do I use the English version of QGIS.


texts: Configuring geonode proxied by nginx on webfaction

texts:
I'm trying to install geonode on webfaction.
Webfaction doesn't provide root access, so it needs installation from source.
I'm struggling to get authentication working.

Webfaction works like this (simplified a bit):

using their control panel, create an "application". This reserves a port.
also via control panel, create a "website". This maps a domain to this port. This happens via a system nginx instance which proxies all requests to the assigned port.

I created an nginx instance listening on this port, proxying to geonode (i.e. django, served using uwsgi) and geoserver (i.e. a tomcat webapp).

So the total setup is:

webfaction nginx -> my nginx -> / (django)
                             -> /geoserver/ (tomcat)

When loading the geonode homepage, tomcat logs as follows:

2014-05-08 05:49:06,865 WARN [geonode.security] - Error connecting to the GeoNode server for authentication purposes
org.springframework.security.authentication.AuthenticationServiceException: Communication with GeoNode failed
        at org.geonode.security.GeoNodeAuthenticationProvider.authenticate(GeoNodeAuthenticationProvider.java:61)
[...]
Caused by: java.io.IOException: GeoNode communication failed, status report is: 404, Not Found
        at org.geonode.security.HTTPClient.sendGET(HTTPClient.java:73)

Here is my nginx config:

pid /home/.../geonode2_nginx/parts/nginx-geonode/nginx-geonode.pid;
lock_file /home/.../geonode2_nginx/parts/nginx-geonode/nginx-geonode.lock;
daemon off;
events {
    worker_connections 1024;
}
http {
    server {
        listen 9999;
        access_log /home/.../logs/user/access_geonode2_nginx.log;
        error_log /home/.../logs/user/error_geonode2_nginx.log;
        # serve staticfiles
        location /static/ {
            include       /home/.../geonode2_nginx/parts/nginx/conf/mime.types;
            alias /home/.../geonode2_static/;
            expires 7d;
        }
        # serve mediafiles, default 'uploaded' in GeoNode
        location /uploaded/ {
            alias /home/.../geonode2_upload/;
            expires 7d;
        }
        # geoserver proxy
        location /geoserver/ {
            proxy_pass http://127.0.0.1:8888/geoserver/;
            proxy_set_header X-Real-IP $remote_addr;
            proxy_set_header X-Forwarded-For $proxy_add_x_forwarded_for;
            proxy_set_header Host $http_host;
        }
        # gunicorn/uwsgi wsgi proxy
        location / {
            proxy_pass_header Server;
            proxy_set_header Host $http_host;
            proxy_redirect off;
            proxy_set_header X-Real-IP $remote_addr;
            proxy_set_header X-Scheme $scheme;
            proxy_connect_timeout 20;
            proxy_read_timeout 20;
            proxy_pass http://127.0.0.1:7777/;
        }
    }
}

Can anyone see anything wrong?

Any tips on debugging this? What HTTP headers should I look for?

How can I find out what the 404 in the tomcat message above refers to?

I asked this question on the webfaction forums also at
https://community.webfaction.com/questions/16732/geonode-configuration


texts: Converting ESRI SHP to MapInfo TAB via Python QGIS 2.18.1

texts:
I am attempting to convert a few files from ESRI Shapefiles to MapInfo Tab via python. As I understand, the 'gdalogr:convertformat' processing algorithm can be used for this. My code looks something like this at the moment...

def convert_layer(layer,filetype, save_location):
    from processing.core.Processing import Processing
    Processing.initialize()
    print "Converting", layer, "to Mapinfo TAB..."
    processing.tools.general.runalg('gdalogr:convertformat',layer, filetype, "", save_location)

convert_layer("path/to/layer.shp", 5, "path/to/tab")

When I run this, QGIS reports that the conversion was successful, however no files are placed in the output folder.

To make matters more confusing, this same process seems to work fine in QGIS 2.14.

Has something changed between 2.14 and 2.18?


texts:
Where is the Angle Measure Tool in ArcMap 10?

I can't seem to find it.


texts: How to Install QGIS Server on Debian Stretch

texts: I am using QGIS version 3.4. I have downloaded training manual of same version. right now I am try to complete 12 chapter i.e., install QGIS server. Here I dnt understand where to give command or expression.


texts:
I want to export my SHPFile into microstation dgn using qgis. But when, I do click right on layer and save as... This message error shows "L'export du fichier vectoriel a échoué.
Erreur : La création du champ Type a échoué (Erreur OGR : CreateField() not supported by this layer"
That tells me that the export has failed !! :/


texts: What are the pros and cons of PostGIS geography and geometry types?

texts:
My company uses geometry (the_geom) data type for storing geospatial data.

I've recently been acquainted to the concept of geography (the_geog) data type which, as I understand it, stores the SRID along with the geometry.

What are the differences between geography and geometry, and is there any advantage of using one of them in large databases?


texts: Creating arcs between two points in ArcPy?

texts:
I've been working on an python script automating the visualization of live point data. I have a table with two coordinate pairs. I'd like to connect these two point with an arc of a circle or similar (parabola??).

I've been able to connect the two points with a straight line, but it gives me a boring visual.

One of the major hurdles is my license level: ArcView.

Anyone have any idea how I can generate a set of points that represent the path of a line between these two points?

I can then use the points to polyline command.


texts: Priority of street labels over amenity icons

texts:
OK, so I have a CartoCSS project.

These are all of the stylesheets used (as described in my project.mml):

  "Stylesheet": [
    "palette.mss",
    "base.mss",
    "roads.mss",
    "labels.mss",
    "extras.mss"
  ],

In the labels.mss file, I'm styling the labels. Here's the part I'm interested in:

#minorroad_label[zoom>15] {
  text-name:'[name]';
  text-face-name:@sans;
  text-placement:line;
  text-size:9;
  text-fill:@street_text;
  text-halo-fill:@road_halo;
  text-halo-radius:1;
  text-min-distance:60;
  text-size:11;
  text-allow-overlap: false;
}

#minorroad_label[oneway!=0][zoom>=16]{
  marker-placement:line;
  marker-height:1;
  marker-fill: #bbb;
  marker-line-opacity:0;
  marker-spacing: 250;
}

/* ================================================================== */
/* ONE-WAY ARROWS
/* ================================================================== */
#motorway_label[zoom>=16],
#mainroad_label[zoom>=16],
#minorroad_label[zoom>=16] {
  [oneway = 'yes'],
  [oneway='-1'] {
     marker-placement:line;
     marker-max-error: 0.5;
     marker-spacing: 200;
     marker-file: url(img/oneway.svg);
     marker-fill: white;
     [oneway='-1'] { marker-file: url(img/oneway-reverse.svg); }
     [zoom=16] { marker-transform: "scale(0.5)"; }
     [zoom=17] { marker-transform: "scale(0.75)"; }
  }
}

In the extras.mss file, I'm setting up the amenity icons. Here's an example:

.points {

    ::amenity {

        [amenity = 'arts_centre'][zoom >=17] {
            marker-transform: scale(@icon_scale);
            marker-placement: interior;
            marker-fill: @icon_fill;
            marker-file: url('icons/shop/art.svg');
        }

        [amenity = 'atm'][zoom >=17] {
            marker-transform: scale(@icon_scale); marker-placement: interior; marker-fill: @icon_fill;
            marker-file: url('icons/atm.svg');
        }

        //[amenity = 'bank'][zoom >=17] {
        //  marker-placement: interior; marker-fill: @icon_fill;
        //  marker-file: url('icons/bank.svg');
        //}

        [amenity = 'bar'][zoom >= 17] {
            marker-transform: scale(@icon_scale); marker-placement: interior; marker-fill: @icon_fill;
            marker-file: url('icons/bar.svg');
        }

        [amenity = 'biergarten'][zoom >= 17] {
            marker-transform: scale(@icon_scale); marker-placement: interior; marker-fill: @icon_fill;
            marker-file: url('icons/biergarten.svg');
        }

        [amenity = 'cafe'][zoom >= 17] {
            marker-transform: scale(@icon_scale); marker-placement: interior; marker-fill: @icon_fill;
            marker-file: url('icons/cafe.svg');
        }

        [amenity = 'pub'][zoom >= 17] {
            marker-transform: scale(@icon_scale); marker-placement: interior; marker-fill: @icon_fill;
            marker-file: url('icons/pub.svg');
        }

        [amenity = 'restaurant'][zoom >= 17],
        [amenity = 'food_court'][zoom >= 17] {
            marker-transform: scale(@icon_scale); marker-placement: interior; marker-fill: @icon_fill;
            marker-file: url('icons/restaurant.svg');
        }
    }
}

Now the problem is, in places where there are a lot of cafes or restaurants for example, all I can see the related icons -- with the street names, nowhere to be found.

What it looks like with the above code

What it looks like if I remove the icons altogether

What I want to do is: make the street names appear no matter what, and IF there is enough space for the icons, then that's fine. But the street names (#minorroad_label namely) should have priority over the icons.

I've tried moving things around, changing the order of the stylesheets, putting the icons' part before the labels, or after the labels, as if it'd affect the priority, but to no avail. I've also played with text-allow-overlap and marker-allow-overlap. Well, setting them to true does make them appear, but then it's a mess; definitely not what I want.

Any ideas/suggestions?


texts: Can I convert ESRI Grid (DEM) raster to USGS DEM?

texts:
There is a tool for DEM to Raster in ArcGIS but I need to convert my ArcGIS DEM rasters (mostly LiDAR) into USGS DEM.  The idea is to import this DEM into Sketchup and create a "sandbox" TIN surface and for this SketchUp needs the data to be in this particular format.

In the past I mostly used the Sketchup Plugin to export TIN from ArcGIS to Sketchup which worked great but since the plugin no longer works in version 10.  I am trying to create a new workflow but I am not having much success.  The "Collada method" does not work well for complex surfaces. It requires to first create a multipatch and there is the first issue as described here. Besides, it is painfully slow if the conversion even succeeds at all.

There is another option described here but I have not tried it. Can anyone comment on how well the Cloud v8 plugin mentioned there works?


texts: I need to create one mile increment dividers on a river system to separate out data. Is this possible and if so how?


texts: How can I split multiple polygons using a polyline or any other feature in arcgis?

texts: I am trying to split this polygons of the districts by using the polyline that passes through the district into multiple polygons and get the area of each of these polygons. I have tried the analysis toolbar split option but, it did'nt work plus there is only option for selection of one polygon only


texts: Using leaflet with large data - ACS data for all counties in U.S

texts:
I'm plotting ACS data (counts of low-income children) for each county across the entire U.S.  I don't have any markers, just a county shapefile and a value for each county in a popup.

Unsurprisingly this is too large for browsers to handle.  I've seen some similar questions here that deal with clustering, but since I have nothing to cluster I'm starting a new question.

I've tried breaking it up into 4 regional maps, but even those are too much.

Is there some way to lower the footprint, but still have all the county specific data?

I'm using R to generate the map, so R-specific responses are helpful but not necessary.


texts: Map is blurry when mobile viewport is set

texts:
we are creating a map using openlayers 2. It is working fine on PC, but on mobile devices map looks blurry. I also tried using incognito mode but I get same result, so there is no problem with cached tiles.
Viewport is set as:

<meta name="viewport" content="width=device-width, height=device-height, initial-scale=1, minimum-scale=1, maximum-scale=1">

Screenshots from iphone,
this is without viewport set:

and this is using viewport from above:

Do you have any solution?


texts: Free software to publish tiles on Amazon S3?

texts: Does anyone know of free software that would allow me to publish tiles on Amazon S3? I know of Arc2Earth, but it is not free.


texts: What is wrong with this Proj4 transformation?

texts:
I'm implementing support for WMS maps in my company's software.  Servers return a list of projections their maps are available in, but our software only displays data in WGS84 (lat/lon or UTM), and not all servers return the "CRS:84" or "EPSG:4326" CRSes.  As a result I need to transform / reproject the data.  Some servers return quite unusual ones, and I'm having trouble with one of these: EPSG:3112 (GDA94 / Geoscience Australia Lambert) used by Geoscience Australia's WMS server.

In order to display a map onscreen, I need to convert the WGS84 lat/lon coordinates of the required screen or tile bounding box to the server's CRS, request the map image, and then to draw convert the map's bounding box (in the server's CRS) back to WGS-84 lat/lon.  This does not seem to be working quite right:

This image is a portion of the Australian island of Tasmania.  The grey coastline is from Demis, in WGS:84 projection.  The blue outline and grey graticule is from Geoscience Australia, in EPSG:3112, but (theoretically) reprojected through Proj4 to match the grey coastline.  As you can see it's quite inaccurate, to the tune of 50-100 kilometres.  If this display were drawing WGS-84 graticules, they would be perfectly vertical and horizontal.  Note also that the island is stretched, i.e. it's more than a rotation.  This looks to me as though the reprojection is not working at all.

The two Proj4 initialisation strings I'm using are:

WGS-84 (WMS's CRS:84): +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs
EPSG:3112: +proj=lcc +lat_1=-18 +lat_2=-36 +lat_0=0 +lon_0=134 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs

Round-tripping an individual point from WGS-84 through EPSG:3112 and back yields the original point:

WGS:84: (145.0, -42.0)
to EPSG:3112: (931684.430679232, -4764432.89032608)
and back to WGS-84: (145, -41.9999999999999)

All conversions are done using pj_transform.

I don't know enough about map projections or Proj4 to know where to go next, although I have been learning a lot over the past few days, and I'd really appreciate advice from someone who knows a lot more about this library.

Some ideas I've had are:

The blue island appears to be still in Lambert projection.  Is this likely?  It would explain the visual distortion.
We use a version of the Proj4 DLL that was built years ago, and while it appears you can build various data about datums etc into the DLL, I don't know whether this was done.  Proj4 does not flag an error doing the coordinate conversion, though, so at the moment I'm working on the assumption that if the library is happy, it's an external error (i.e. on my part.)  I also don't know the reasoning for having our own, years-old build of the DLL anyway, instead of using a prebuilt binary from the Proj4 website.  However, even the Proj4 binaries are one version old and don't seem to be built to include everything as binary data.
Some versions of Proj4 seem to load datum and other information externally.  I have tried replacing our DLL and configuring it to use some of these files using the PROJ_LIB environment variable, but with no apparent change.
The +towgs84=0,0,0,0,0,0,0 part of the EPSG:3112 string seems odd.  This is from the 'epsg' file in the latest Proj4 download, but could it be wrong?

Perhaps a better formulation of this question is:

What do I need to have or do to convert between these two CRSes?

Thanks for any help!

  I originally posted this question on Stack Overflow, but I've
  re-posted it here.  I think it's more a GIS question than a
  programming question.


texts: Given a set of coordinates, how can I calculate the minimum bounds?

texts:
I have a set of coordinates in lon/lat. There are always at least 3 ordered coordinates that form a poly (if they were projected flat). How can I calculate the minimum bounds for these coordinates as a set of valid longitude and latitude ranges? (by 'valid' I mean ranges that explicitly account for the antimeridian). Its kind of hard to explain what I'm looking for so I have a picture.

The case in the picture is trivial; you just find the absolute minimum and maximum for all the points. This doesn't work for all cases though. Is there a generic solution?

Edit:
To clarify what I mean by 'valid', say I had three longitude values in my data set: -76, -135 and 164. The values cross the antimeridian and I'd like the resulting ranges to be split: -76 to -180 AND 164 to 180.

Some more clarification. The points form a polygon, so in certain cases, the required range could be from -180 to +180 (ie, the full 360 degrees):

The image on the left shows the longitude of four coordinates that occur on one 'half' of the Earth. Imagine it as if you were looking down onto the north pole (black dot). The pink shows the minimum longitudinal range that encompasses the polygon (the polygon is shown between the four points in purple). The case on the left would have two longitudinal ranges: [-180 to -120] and [135 to 180] (just visually estimating it)

The image on the right shows another case where the points go all the way around the Earth. This range would be [-180 to 180].


texts: Converting Sentinel-3 data (netcdf) to GeoTIFF?

texts:
I'm trying to convert Sentinel-3 data to GeoTiff file. The Sentinel-3 data comes as separate nc files (netcdf): a file for each spectral band and a nc file for qualityFlags, geo_coordinates etc' and an xml file. The official why to read Sentinel-3 data is to use the SNAP software and loading to SNAP the xml file that directs SNAP to the different files. I would like to do it using a python script, but I can't figure how. I tried GDAL translate but it does not make the connection between the spectral band file and the coordinate file.  I am able to read each nc file separately and to read the coordinates nc file as text/array.

How can I combine the spectral band data with the coordinate data to create a Geotiff?

My guess is that it is possible doing that by using GDAL but I can't find the way. Is anyone familiar with Sentinel-3 data?
I'm using python 2.7 on win10, I have GDAL 2.1


texts: ST_DWithin uses Geometry or Geography type

texts:
assuming I use the following query to match the locations to a poi within 10km distance:

SELECT
poi.poi_id,
h.loc_id,
(ST_Distance(poi.geom,h.geom)*111.111) as distanceinkm
FROM iut_poi_transl AS poi, hrs_loc as h
WHERE poi.poigr_id != 30 AND ST_DWithin(poi.geom,h.geom,0.09)

Data is SRID 4326.

Does it use the geometry or the geography type?

If uses geometry how do I have to change the query to get use geography type?


texts: How can I change an attribute table field type from string to double in QGIS?

texts:
I work on QGIS 2.2, and I am unable to edit my attribute table to change the field type from string to double.

I have gone through several tutorials.
- I can not run the table manager plugin, probably because its latest version is for qgis 2.0, (well it appears broken on my plugin manager UI).

I am unable to edit my layer so that I can use the field calculator (as several tutorials stipulated).
I am not versed with python to use the python console


texts: Layer files not drawing properly in ArcMap?

texts:
I have a co-worker who had his windows user credentials changed after coming on full time. After this change we noticed an issue with a series of layerfiles not drawing properly when brought into ArcMap.

This only seems to be an issue for the layer files that are referencing ArcGIS Server connections. For example, we have a layerfile that references imagery from the Virginia (VGIN) server and when he brings this in the layer file has no error, it is just completely empty.

No one else in my organization has this issue, only my one co-worker. This makes me believe that it is related to his Windows user credentials being changed. His "C:\Users\UserName\AppData\Roaming\ESRI" folder has been updated with his new user name and I am not able to find any trace of the old user name being referenced anywhere.

Is anyone able to point me in the direction to resolve this issue?


texts: How do I create smaller maps for publishing from multiple tiles representing a large map using QGIS?

texts:
I am completely new to cartography and QGIS - I'm a writer, editor and acceptable designer (Adobe Creative Suite).

I wish to produce PDF tour guides featuring maps that display a route. One version of the PDF will be for print (so high resolution images - 300 dpi) and one version for digital.

I have downloaded Ordnance Survey District Raster data. The data folder I am working with is 75.1MB and contains 68 GeoTIFF Full Colour tiles. The area covered is north west Wales (10,000 square kilometers) but the route I wish to represent passes through a much smaller area and I wish to break the route down into a series of about 10 maps, each covering an area not much greater than 200 square kilometers.

I want to be able to find the areas to be represented, save them as files that I can, if necessary, work on in photoshop or illustrator to show the route i.e. a coloured line overlaid - I have gpx points and am aware they could be loaded via QGIS but perhaps that's a task for later.

To achieve this, my assumption, based on what I can understand from reading QGIS tutorials, is to load all the tiles and then merge them. Once that has been achieved I then clip the area I wish to represent, and save it as a picture file. If that is not correct please tell me what i should be doing.

If it is the correct approach then please can you tell me what the error message "gdal_merge.bat was unexpected at this time" means and what I can do about it? It appears when I try to merge all of the tiles and the merge does not take place.

If I merge half of the tiles then that works so I assume it has something to do with file size.

Also the merged image is a black background (land and sea) and white features (roads, towns etc).  As it only merges some of the tiles then some of the resultant map in the QGIS interface is black and white (like a photographic negative) and the rest (non merged tiles I assume) is colour. In the style section of properties I noted that the unmerged tiles are designated 'palette' and the merged part is designated Singleband gray.

Any assistance to successfully get me out of this quagmire and back onto solid ground would be very welcome.


texts:
Is it possible to make utfgrids from a raster source (e.g. GeoTIFF)?

I am trying to get utfgrids from a .tif with mapnik.xml and TileStache but it seems not to be working. As a workaround I am polygonizing the tif to a shape which works well but is not very efficient.

Has anyone experince with this question? Maybe with other tileservers?

This would be an example for the TileMill generated mapnik.xml I am working with:

<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE Map[]>
<Map srs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over" maximum-extent="-20037508.34,-20037508.34,20037508.34,20037508.34">

<Parameters>
  <Parameter name="bounds">5,43,18,50</Parameter>
  <Parameter name="center">16.7693,48.4365,9</Parameter>
  <Parameter name="format">png8</Parameter>
  <Parameter name="minzoom">0</Parameter>
  <Parameter name="maxzoom">15</Parameter>
  <Parameter name="scale">1</Parameter>
  <Parameter name="metatile">2</Parameter>
  <Parameter name="id"><![CDATA[test]]></Parameter>
  <Parameter name="_updated">1409125380000</Parameter>
  <Parameter name="name"><![CDATA[tmp2m_2_grid]]></Parameter>
  <Parameter name="tilejson"><![CDATA[2.0.0]]></Parameter>
  <Parameter name="scheme"><![CDATA[xyz]]></Parameter>
</Parameters>

<Style name="utf_grid" filter-mode="first">
  <Rule>
    <RasterSymbolizer opacity="1" scaling="near"/>
  </Rule>
</Style>
<Layer name="utf_grid"
  srs="+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over">
    <StyleName>utf_grid</StyleName>
    <Datasource>
       <Parameter name="file"><![CDATA[/path/to/test.tif]]></Parameter>
       <Parameter name="band"><![CDATA[1]]></Parameter>
       <Parameter name="type"><![CDATA[gdal]]></Parameter>
    </Datasource>
  </Layer>

</Map>


texts: QGIS 3.8 could not start ''Qgis_app.dll is missing'' check runtime environment check C:\OSGEO4W64\bin\qgis-dev-bin-g.exe

texts:
I had Qgis 3.4 running well on my PC Windows 10, 64 bit. I tried to upgrade to Qgis 3.8 Zanzibar and I had the error message:
''Qgis_app.dll is missing check runtime environment check C:\OSGEO4W64\bin\qgis-dev-bin-g.exe''.

So far, I tried the folowwing:

install countless times with OSGEO
install from OSGeo and direct download, the 3.8 don't work
go back to the older 3.4 version (from OSGEO or direct donwload), it didn't work
run the from the qgis.bat file
installed qgis 3.4 as administrator
copy and paste the qgis.bat.temp and change the name to qgis.bat
Tried to change the path of OSGeo in the .bat file, but I was a little bit lost...?

OSgeo4W shell does not launch from my start menu either, instead I have a text file opening with a text editor software.
Grass Gis 7.2 doesn't work also, I also have a text file opening with a text editor software.

Funnily, I have two folders created for QGIS (or I believe so, but I'm not sure):

C:\OSGeo4W64\bin

and:

C:\Program Files\QGIS 3.8\bin

Could this be one of the reason? Should I append one of the path file, if yes, how should I do it?


texts:
I have a table with clients containing a field with the name of the districts in which they live.
The client The data I have imported from an Access database and is now as a table in PosgreSQL / PostGIS but no geometry exists.
I have also one shapefile with all districts that is already imported as a Table in POstgreSQL / PostGIS.

Whats is the Sql code so as to have each clients represented as a point based on district names?

Att: There is a field on each table with the names of distrits.


texts: GeoServer not rendering GeoTIFFs after JPEG compression

texts:
I previously asked about some problems I was having converting JP2/MrSid data to GeoTIFF for use in GeoServer. It was suggested that I needed to add overview data to the 4 GB-sized GeoTIFF -- and that helped. GeoServer then rendered and served the raster and was quite fast about it.

I've since been reading about cloud-optimized GeoTIFFs, and tried to follow the three-step process recommended here:

gdal_translate sample.jp2 sample_deflate.tif -co tiled=yes -co BLOCKXSIZE=512 -co BLOCKYSIZE=512 -co COMPRESS=DEFLATE -co PREDICTOR=2

cp sample_deflate.tif sample_addo.tif

gdaladdo -r average sample_addo.tif 2 4 8 16 32 64 128 --config PHOTOMETRIC_OVERVIEW YCBCR --config COMPRESS_OVERVIEW JPEG --config JPEG_QUALITY_OVERVIEW 85 --config INTERLEAVE_OVERVIEW PIXEL

gdal_translate sample_addo.tif sample_compressed.tif -co TILED=YES -co BLOCKXSIZE=512 -co BLOCKYSIZE=512 -co COMPRESS=JPEG -co JPEG_QUALITY=85 -co PHOTOMETRIC=YCBCR -co COPY_SRC_OVERVIEWS=YES --config GDAL_TIFF_OVR_BLOCKSIZE 512

But the resulting file, after JPEG compression with photometric, would not serve at all in GeoServer. The intermediate file -- with just the added overviews -- serves fine in GeoServer, but is 6 GB in size.

Thinking the problem might be the photometric, I backed up and redid the compression with no photometric:

gdal_translate sample_addo.tif sample_compressed.tif -co TILED=YES -co BLOCKXSIZE=512 -co BLOCKYSIZE=512 -co COMPRESS=JPEG -co JPEG_QUALITY=85 -co COPY_SRC_OVERVIEWS=YES --config GDAL_TIFF_OVR_BLOCKSIZE 512

It still wouldn't work. So something to do with my third step, compressing the image, is rendering it unusable by GeoServer. (It could also, I guess, be a problem with my GeoServer settings.)

Here is the gdalinfo for the third (not working) file (with photometric compression):

Driver: GTiff/GeoTIFF
Files: sample_compressed.tif
Size is 43221, 36021
Coordinate System is:
GEOGCS["WGS 84",
    DATUM["WGS_1984",
        SPHEROID["WGS 84",6378137,298.257223563,
            AUTHORITY["EPSG","7030"]],
        AUTHORITY["EPSG","6326"]],
    PRIMEM["Greenwich",0],
    UNIT["degree",0.0174532925199433],
    AUTHORITY["EPSG","4326"]]
Origin = (95.998541666740138,45.001424215510333)
Pixel Size = (0.000138888892252,-0.000138888892252)
Metadata:
  AREA_OR_POINT=Point
Image Structure Metadata:
  COMPRESSION=YCbCr JPEG
  INTERLEAVE=PIXEL
  SOURCE_COLOR_SPACE=YCbCr
Corner Coordinates:
Upper Left  (  95.9985417,  45.0014242) ( 95d59'54.75"E, 45d 0' 5.13"N)
Lower Left  (  95.9985417,  39.9985074) ( 95d59'54.75"E, 39d59'54.63"N)
Upper Right ( 102.0014585,  45.0014242) (102d 0' 5.25"E, 45d 0' 5.13"N)
Lower Right ( 102.0014585,  39.9985074) (102d 0' 5.25"E, 39d59'54.63"N)
Center      (  99.0000001,  42.4999658) ( 99d 0' 0.00"E, 42d29'59.88"N)
Band 1 Block=512x512 Type=Byte, ColorInterp=Red
  Overviews: 21611x18011, 10806x9006, 5403x4503, 2702x2252, 1351x1126, 676x563, 338x282
Band 2 Block=512x512 Type=Byte, ColorInterp=Green
  Overviews: 21611x18011, 10806x9006, 5403x4503, 2702x2252, 1351x1126, 676x563, 338x282
Band 3 Block=512x512 Type=Byte, ColorInterp=Blue
  Overviews: 21611x18011, 10806x9006, 5403x4503, 2702x2252, 1351x1126, 676x563, 338x282

Any ideas what I'm doing wrong or how I can get GeoServer to serve this?


texts:
I'm not very skilled on using postgis but I face a problem that I must solve, i.e. create tables to load geographical data stored as dxf.

I have a thousand dxf files (3d) with polylines, points and polygons and I want to convert them to a postgis/postgres db.

I'm in windows xp using postgres 9 postgis 2.0.

Gdal/OGR are installed using Osgeo4Win installer.
I can write some basic python script to handle data in batch.

Can you please direct me to some useful tutorial/hint on how to face this problem.
Thanks in advance


texts: Histogram from multiple raster in one graph in QGIS?

texts: I have more than 10 rasters which are actually MNDWI index images created from Landsat image. I want to visualize all rasters (tiff) files in one graph to understand the spectral distribution. Is there any plugin in QGIS to do that?


texts: Calculating the area covered by isochrones in QGIS

texts: I created isochrones using Iso-Area algorithm as polygons in QGIS. Is there any tool in QGIS that helps with calculating areas covered by created polygons?


texts: QGIS - Transfer feature selection with new attributes to PostGIS

texts:
I have three shapefiles which represent a conduit network with valves and hydrants. The pipe grid consists of line features in a shapefile and the valves and hydrant each are a point layer in a shapefile.

For the planning of flushing proceedings I need to select some of the pipe elements (the line features) in a certain order which are to be flushed. The underlying attribute table of the pipe elements is also important for the planning process.

Ideally there is a way to create many different flushing selections (e.g. pipe no. 1-5, pipe no. 18-45, etc.) and save them in another layer, let's call it flushings. It is essential that the selection of pipes is assigned a new identifier the flushing ID.

So the flushing layer would contain flushing no. 1,2,3... with each flushing itself containing or refering to a previos made selection of pipe elements.
I then would like to export the generated information of a pipe selection and an assigned flushing ID to either a shapefile or any other way which I can use in my PostGIS database.

Example:

pipe elements attribute table:

id  material  diameter
1   PE        80
2   PE        80
3   PE        100
4   PE        100

I need something like this in my database:

flushings

flushing_id   pipe_id  order
f1            2        1
f1            3        2
f1            4        3
f2            1        1
f2            2        2
f2            3        3

Is there an easy workflow for pipe selection and assigning a flushing id to the selection to import it in a database later? It would also be helpful if each selection could be colored differently to make it easier to work with.


texts: Failed to publish globe service on ArcGIS Server. Error 001369

texts:
I'm trying to publish a globe service on our ArcGIS Server and in ArcGlobe I get this error message

  "Failed to publish service."

---> Went to Results:
[h=1]001369 :

  Failed to create the service.[/h] [h=2]Description[/h]The creation of
  the service failed. This is a generic error. Many different conditions
  could cause this failure. Check the server logs for the specific
  issue.

[h=2]Solution[/h]Check the server logs for the specific issue.

--> Went to ArcGIS Server Logs:

  SEVERE    22.08.2013, 10:15:20    Failed to construct instance of service
  'glob_4.GlobeServer'. Server SEVERE   22.08.2013, 10:15:20    Failed to
  initialize server object 'glob_4': 0x80030003: null   Server
  SEVERE    22.08.2013, 10:15:20    (null)  glob_4.GlobeServer

Unfortunately I really don't know how to solve that. The service is published (I can see it on the server) but apparently it cannot start. I also tried saving a service definition file and then publishing it on the server-> same problem.

So thanks for any advice.

PS: We use ArcGIS Server 10.1 on Linux


texts: Creating universal ArcGIS 'Current tool' property?

texts:
I have a base tool class which handles 90% of my custom logic, but some tools need to set themselves as the current tool. The operations between ArcGIS Engine's IToolbarControl and ArcGIS Desktop's IApplication interface are different.

With ArcGIS Desktop, I need to have a reference to ESRI.ArcGIS.Framework, which is unavailable in a pure ArcEngine deployment.

Is there any way to set the Current Tool property without needing a reference to ESRI.ArcGIS.Framework?

As a side note: Any way of getting the ArcGIS Desktop's window handle without resorting to the Win32 API?


texts: How do I use HTML5 Storage to cache mapservice data?

texts:
I would like to develop a web app that runs on any device (including iPad and iPhone) in a disconnected mode.  It is OK if the user is required to browse around and download mapservice data to a local cache before disconnecting and heading to the field.

I would like to do this with Esri's javascript SDK.  (There is an existing app for which the mapservices have been set up, however it is written in Silverlight.  I've been asked to port it to work with the iPad in a disconnected mode).

Caching ArcGIS Online tiles is a violation of the terms of use.  Are there any restrictions on caching OpenStreetMap?  If not, then how could I cache OSM locally and use that instead of the web service while disconnected.

For the vector data, I would like to provide a button that runs a query and downloads the json locally.  Ideally I'd like to transparently redirect subsequent calls to fetch json from local storage, but I realize this is a bit naive.  Does anyone know any good strategies for doing that?


texts: Tracing polyline between two points in ArcGIS Desktop?

texts:
How i can draw a line between two points...e.g.: I have initial point and endpoint
and i need trace a line between this points...how i can do that?


texts: PostGIS 2.5 error while compiling from source on Ubuntu Server 18.04

texts:
I'm compiling PostGIS 2.5 from source on Ubuntu Server 18.04

The server had no internet connection due to client requirements, so I've been following PostGIS 2.5.1dev Manual, 2.4. Compiling and Install from Source: Detailed Steps ( https://postgis.net/docs/manual-2.5/postgis_installation.html#PGInstall ).

It's required to install (compile from source) the next packages:

Proj4 reprojection library, version 4.9.0 or greater.
GEOS geometry library, version 3.5 or greater.
LibXML2, version 2.5.x or higher.
JSON-C, version 0.9 or higher.
GDAL, version 1.8 or higher.

Everything works fine after typical compiling instructions:

./configure
make
make install

Then I compile PostGIS again with typical instructions:

./configure
make
make install

In order to check if everything is working, I connect to my DB:

sudo -u postgres psql postgres
postgres=# \c testdb
testdb=# \CREATE EXTENSION postgis;

And then...
PostGIS ERROR: could not load library "/usr/lib/postgresql/10/lib/postgis-2.5.so": libproj.so.13: cannot open shared object file: No such file or directory.

As far as I understand, the problem is with libproj so I've compiled again, but had no success.

Did I miss something?


texts: Get and compare area of arcpy geometry from extent object

texts:
I'm trying to get a numerical value for the area of a polygon that I am generating from an extent object.  In this case, I am using the "onRectangle" function in the python addins for arcpy.  The problem starts when using the getArea method on the polygon geometry.  I have no intention of writing this polygon to disk.  Ultimately, I want to be able to compare the extent of a user drawn rectangle to the dataframe extent. Here's what I've got so far.

def onRectangle(self, rectangle_geometry):
    frameExtent = rectangle_geometry
    XMAX = frameExtent.XMax
    XMIN = frameExtent.XMin
    YMAX = frameExtent.YMax
    YMIN = frameExtent.YMin
    pnt1 = arcpy.Point(XMIN, YMIN)
    pnt2 = arcpy.Point(XMIN, YMAX)
    pnt3 = arcpy.Point(XMAX, YMAX)
    pnt4 = arcpy.Point(XMAX, YMIN)
    print XMAX, XMIN, YMAX, YMIN
    array = arcpy.Array()
    array.add(pnt1)
    array.add(pnt2)
    array.add(pnt3)
    array.add(pnt4)
    array.add(pnt1)
    rectPoly = arcpy.Polygon(array)
    global rectArea
    rectArea = rectPoly.getArea()
    print str(rectArea)
    return rectArea

The python window in ArcMap returns this error: TypeError: Invalid geometry type for method.  Any pointers would be appreciated for both getting the area of the polygon and comparing it (percentage wise) to the full data extent would be appreciated.

I'm using ArcGIS 10.2 advanced license.


texts:
I'm trying to create a custom symbol using inkscape and editing the svg file to be able to manage it from the style window in QGIS.

Edited the .svg symbol file deleting the "style" lines and inserting the line (after "path" and "circle"):

fill="param(fill) #FFF" stroke="param(outline) #000" stroke-width="param(outline-width) 1"

I can edit lines colors/thicknesses and the fill color but I do not know how to edit file to get a transparent fill of the symbol. How can i fix it?


texts: Image Collection isn't able to be manipulated after mapping function

texts:
This code works well until the last two blocks. The second to last block multiplies all raster pixel values in an image collection by 10. The last block is supposed to parse the entire collection into separate collections by year. What currently happens is that the image collection seems to work for everything but the last block. When I try to put SAR100 into the parsing block, it returns zero rasters band for each year when there should definitely be more. I receive no error messages when running the code. I will include a GEE link for this code. The study area, named "geometry", is an arbitrary small polygon.
//Changeable Parameters
var startYear = '2016';
var startDate = '-04-01';
var endYear = '2018';
var endDate = '-09-30';
var polarization = 'VH';
var studyRegion = geometry;

//--------------------------------------------------------------------
//          Modify Sentinel-1 collection to desired specifications
//--------------------------------------------------------------------

//Select image collection, filter it to desired polarization,
  //instrument mode, overall start/end dates, and orbital direction
var imgVV = ee.ImageCollection('COPERNICUS/S1_GRD')
        .filter(ee.Filter.listContains('transmitterReceiverPolarisation', polarization))
        .filter(ee.Filter.eq('instrumentMode', 'IW'))
        .select(polarization)
        .filterBounds(studyRegion)
        .filterDate( '2017-01-01', '2020-12-31')
        .filter(ee.Filter.eq('orbitProperties_pass', 'ASCENDING'));

//Clip image collection to study area polygon
var clippedSAR = imgVV.map(function(image) { return image.clip(studyRegion); });
print(clippedSAR, 'clippedSAR');

//--------------------------------------------------------------------
//          Multiply collection by 10
//--------------------------------------------------------------------

var SAR100 = clippedSAR.map(function(image) {
  return image.multiply(10);
});
print(SAR100, 'SAR100');

//--------------------------------------------------------------------
//          Parse complete collection into separate years
//--------------------------------------------------------------------

var collections = {};
for (var year = startYear; year <= endYear; year++)
  collections[year] = SAR100.filterDate(year.toString()
  + startDate, year.toString() + endDate).filterBounds(studyRegion);
print ('collections', collections);

https://code.earthengine.google.com/7f906f531ca5ab8aad20ff2dc3cf41c8


texts: Geocoding South African Addresses in ArcMap?

texts:
I have about 8,000 addresses to geocode in South Africa. I haven't been able to find an address locator in ArcMap that works (trying to use the 'world places' address locator in ArcMap 10 gives me a "resource is not available, regardless of authorization" message.

In the past, I was using an online batch geocode tool from mapdevelopers, however this is limited to 100 records at a time.

Does anyone have info on an address locator for South Africa or another way of doing this?


texts: Can I automate a process in QGIS 3.0 that involves a plug in and grass?

texts:
I would like to automate a process in QGIS 3.0.

The process is:

Open raster calculator and calculate a bands operation (Sentinel 2 bands).
Create dots (centroid of pixel) into several pre-existing polygons.
Open the plug in "point sampling tool" and give to the dots, created in the previous step, the value of the corresponding pixel.
Open GRASS tool v.surf.idw and process the algorithm for every group of dots created in step 3.
Cut the raster created in step 4 to the corresponding polygon.
Select some colours for the cut raster choosing the option singleband pseudocolours.

As I have to do this process for several polygons, it's taking a long time so I would like to automate it; is it possible?


texts: What is the complexity (big / little o) of the ArcGIS watershed delineation algorithm?

texts: I'm curious about the complexity, in terms of Big or little O notation, of watershed delineation algorithms that are implemented in ArcGIS.  I imagine that these will be dependent on either the size of the input raster or on the size of the basin being delineated, but I am not sure.


texts: What value to place for data which has no lat/long value?

texts: i'm trying to create a shapefile from an excel spreadsheet. the spreadsheet contains some data such as "unknown record" for which no lat/long can be obtained while the rest of the data is related to states of a country.when i add spreadsheet to arcgis it displays "null" value under lat long field in front of unknown record and when i create the shapefile the null values are replaced by "0" which causes an inconsistent map extent. i don't want to create a csv as it distorts my rest of the data. i tried placing "0" for the unknown record in lat long fields but in that case lat long wont show in add xy dialog box.please suggest a workaround so that my problem of inconsistent map extent is resolved.


texts: Unable to execute QGIS algorithm while taking Input raster directly from AWS S3 account

texts:
I have a raster file (.TIF) in my AWS S3 account. I am using GDAL to open that .TIF file and to save the reprojected .TIF file to my local computer using QGIS GDAL 'Warp (reproject)' algorithm as given below.

ds = gdal.Open('/vsis3/s3bucket/input.tif')
processing.run("gdal:warpreproject", {'INPUT':ds,'SOURCE_CRS':QgsCoordinateReferenceSystem('EPSG:32643'),'TARGET_CRS':QgsCoordinateReferenceSystem('EPSG:4326'),'RESAMPLING':0,'NODATA':None,'TARGET_RESOLUTION':None,'OPTIONS':'','DATA_TYPE':0,'TARGET_EXTENT':None,'TARGET_EXTENT_CRS':None,'MULTITHREADING':False,'EXTRA':'','OUTPUT':'/home/user/Desktop/reprojected.tif'})

I am able to open the input.tif in my AWS S3 account using gdal.Open.

ds
<osgeo.gdal.Dataset; proxy of <Swig Object of type 'GDALDatasetShadow *' at 0x7f4ea96f1750> >

But, while executing the reprojecting algorithm, I am getting error.

_core.QgsProcessingException: Unable to execute algorithm
Could not load source layer for INPUT: invalid value

Meanwhile, when I give a raster file which is in my local computer as input.tif (instead of directly giving the raster file from my AWS S3 account), this same QGIS processing.run algorithm is getting executed.

I wonder, whether the QGIS processing.run algorithms (including 'Warp (reproject)' can only take an input file only from local computer or whether it can take input file from AWS S3 directly as arguments.

I am able to get all raster statistics and other QGIS functionalities by giving the AWS S3 path of the raster as argument. It is only during the execution of processing.run algorithms, I am getting the error.

I understand that QGIS processing.run algorithms will take only INPUT as the path of raster in terms of 'strings', but ds = ('/vsis3/s3bucket/input.tif') =  >. It is not a string, it is a gdal object.

How can I solve this issue ?


texts: Restricting label/line drawing to 45 deg. lines?

texts:
I've followed this tutorial to be able to move label positions and draw lines to connect them to the label dots, but I'd like to restrict the line drawing to 45˚ angles.

Is there a way to do that? I see this in the changelog for QGIS 3.8 but can't figure out how to implement.


texts:
(base) C:\WINDOWS\system32>conda install -c conda-forge geopandas
Collecting package metadata (current_repodata.json): done
Solving environment: failed with initial frozen solve. Retrying with flexible solve.
Solving environment: |

How can I solve this problem?


texts:
I have defined the following SLD rule for my layer. It says that if the property Weekday is greater than or equal to 0 and less than or equal to 3, then the color of a polygon must be #680000.

<?xml version="1.0" encoding="UTF-8"?>
<StyledLayerDescriptor version="1.0.0" xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc"
  xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"
  xsi:schemaLocation="http://www.opengis.net/sld http://schemas.opengis.net/sld/1.0.0/StyledLayerDescriptor.xsd">
  <NamedLayer>
    <Name>green</Name>
    <UserStyle>
      <Name>Adjustable</Name>
      <Title>Triangle</Title>
      <Abstract>Direction and volume</Abstract>
      <FeatureTypeStyle>
        <Rule>
          <Title>Volume control</Title>
          <Abstract>Volume control</Abstract>
          <Filter>
            <And>
              <PropertyIsGreaterThanOrEqualTo>
                <PropertyName>Weekday</PropertyName>
                <Literal>1</Literal>
              </PropertyIsGreaterThanOrEqualTo>
              <PropertyIsLessThanOrEqualTo>
                <PropertyName>Weekday</PropertyName>
                <Literal>3</Literal>
              </PropertyIsLessThanOrEqualTo>
            </And>
          </Filter>
          <PolygonSymbolizer>
            <Fill>
              <CssParameter name="fill">#680000</CssParameter>
            </Fill>
            <Stroke/>
          </PolygonSymbolizer>
        </Rule>
      </FeatureTypeStyle>
    </UserStyle>
  </NamedLayer>
</StyledLayerDescriptor>

When I validate this rule in Geoserver 2.11, I get the following error:

  line 15: cvc-complex-type.2.4.a: Invalid content was found starting
  with element 'Filter'. One of
  '{"http://www.opengis.net/sld":LegendGraphic,
  "http://www.opengis.net/ogc":Filter,
  "http://www.opengis.net/sld":ElseFilter,
  "http://www.opengis.net/sld":MinScaleDenominator,
  "http://www.opengis.net/sld":MaxScaleDenominator,
  "http://www.opengis.net/sld":Symbolizer}' is expected.

Line 15 is <Filter>. I checked the documentation of SLD and it looks like I defined all mandatory fields. Why do I get this error message and how can I solve this error?


texts: How to reference all the shapefiles within a Windows folder and merge them in ModelBuilder?

texts:
I have a group of shapefiles that are stored in a Windows folder. There are multiple users that will be adding shapefiles to this folder or editing existing shapefiles in this folder on a continual basis. Currently, the shapefiles are manually merged together. I would like to build a very simple model that references all of these shapefiles and merges them on a continuous basis- perhaps 1 time per start of business day. The goal is to have a single layer with all of the shapefiles that is updated automatically as new shapefiles are created. All of the shapefiles created will have the same number of fields, the same field names, and same shape (polygon).

I could use Modelbuilder or possibly create a Python script that runs every morning. I am a beginner with Python.

I am using ArcGIS 10.1.

Update #2- My Python code is as follows:

import arcpy, os, time

arcpy.env.workspace = r'C:\Temp\Planning\Planning_In'
fcs = arcpy.ListFeatureClasses(feature_type='Polygon')

fcOut = r'C:\Temp\Planning\Planning_Out\ProjectMerge.shp'
if fcOut in fcs:fcs.remove(fcOut)

if len(fcs) > 0:
    if arcpy.Exists(fcOut):
        fcArchived = os.path.join(r'C:\Temp\Planning\Planning_Out', 'Archived_{}.shp'.format(time.strftime('%Y%m%d')))
        os.makedirs(os.path.dirname(fcArchived))
        arcpy.Copy_management(fcOut, fcArchived)
        arcpy.Delete_management(fcOut)
    arcpy.Merge_management(fcs, fcOut)

print fcs

My shapefiles are named: 13_012, 13_013, etc.

*** My error message is:

Traceback (most recent call last):
  File "C:\Temp\Planning\import arcpy.py", line 16, in <module>
    arcpy.Merge_management(fcs, fcOut)
  File "C:\Program Files (x86)\ArcGIS\Desktop10.3\ArcPy\arcpy\management.py", line 4221, in Merge
    raise e
ExecuteError: Failed to execute. Parameters are not valid.
ERROR 000732: Input Datasets: Dataset C:\Temp\Planning\Planning_In\13-004.shp;C:\Temp\Planning\Planning_In\1314-004.shp does not exist or is not supported
Failed to execute (Merge).

Additionally, when I remove the merge and just apply the following code, I receive an error message that 'fcArchived' is not defined:

    import arcpy, os, time

    arcpy.env.workspace = r'C:\Temp\Planning\Planning_In'
    fcs = arcpy.ListFeatureClasses(feature_type='Polygon')

    fcOut = r'C:\Temp\Planning\Planning_Out\ProjectMerge.shp'
    if fcOut in fcs:fcs.remove(fcOut)

    if len(fcs) > 0:
        if arcpy.Exists(fcOut):
            fcArchived = os.path.join(r'C:\Temp\Planning\Planning_Out', 'Archived_{}.shp'.format(time.strftime('%Y%m%d')))
            os.makedirs(os.path.dirname(fcArchived))
            arcpy.Copy_management(fcOut, fcArchived)

print fcArchived


texts: Conditional styling and expressions using modulus

texts:
I'm trying to style a contour shapefile in TileMill with all contour lines divisible by 100 bolder than others. Can I use expressions in defining a condition? The following example of course does not work, but it gives an idea of what I am trying to get.

#contours {
  line-width:.4;
  line-color:#222222;

  [(elev % 100 = 0)]{
    line-width: 1;
    line-color: #000;
  }

}


texts:
Is it possible to set style layers in the database? (QGIS) [MySQL database]

I would like to set and share my style with other users, so when they open the layer it's still as my style.


texts: What is the correct color map format (.txt) for a color map I want to load in the Raster Layer Properties/Symbology menu of QGIS (v. 3.6.0)


texts: I am wondering is it possible to have scale bar and north arrow in 3D view when viewing in browser and also when saving image.


texts: define startMonth for jqtimeline plugin in a webmap

texts: I created a webmap with timeline function by using leaflet.js and jqtimeline plugin. Do you know a possibility to not just define the startYear but a startMonth with jqtimeline plugin? The reason is that my datasets start in November and the timeline should begin with Nov and not Jan.


texts: Subtract a reference image from a photographed map

texts:
For an analysis of the use of hand-drawn glyphs (shading, circling, diagrams) in a GIS, I have a set of printed maps of the same area, embellished with participants' pencil and pen marks.

I would like to scan or photograph these maps and then extract the glyphs as shapefiles. A diagram of this schema/workflow is shown above. My printed maps have georeferencing features (gridlines).

What is the best workflow to both georeference the glyphs, and subtract the reference map leaving the glyphs themselves as shapefiles? (The shapefiles would then be imported into ArcGIS Pro).


texts: Can I change how Maplex decides which duplicate labels to remove?

texts:
I have a lot of recreation location points, and most of them are labeled. Sometimes there are many points close together which really only need one label. For example, the recreation location below has a campground, picnic area, and boat ramp.

I don't need/want each of these points to say "Woods Ferry" by them, so using the Maplex labeling, I can turn on "Remove duplicates" and get only one label for the set of points. (Yay.)

But, I would prefer it if the label were put on the boat ramp or picnic table instead, since then it would be covering up a river instead of some roads. In this case, the campground point has the lowest OBJECTID, and apparently that is how priority is determined.

Is there any way I can change how the label priority is chosen?


texts: Where to download high resolution maps or layers for Greece?

texts:
I am trying to create management zones on selected fields in Greece, so I am looking for high-resolution layers (maps) to use as a base (preferably satellite imagery).

Could you suggest me the best sites to download high-resolution maps or layers of Greece?

PS. Check attached picture of already known sites.


texts:
I was using PostgreSQL with PostGIS and Pgrouting. Recently I needed to install SFCGAL for some extra functions for my PostGIS db. After installing SFCGAL, now I don't know how to update my PostGIS db to bring those SFCGAL functions into my db?

Paths:

libSFCGAL.so: /usr/local/lib

sfcgal-config: /usr/local/bin

SFCGAL: /usr/local/include


texts: How to create random points outside polygons?

texts: The tool Create Random Point is able to generate a certain number of points within polygons. I am wondering, given a bounding box, is there any way that I can generate random points outside those polygon?


texts: Add filter on attributes having some special Characters in their attribute name

texts:
We have a Feature layer having multiple attributes. One attribute having attribute name is BSC/RNC . we create some filter on that attribute :-

String strQuery = "BSC/RNC >= 10";
org.opengis.filter.Filter flt = CQL.toFilter(strQuery);

but it gives the result with all records having attribute BSC/RNC as NAN values.


texts: How to calculate distances from points to multiple features at the same time?

texts: I have a point layer and I want to calculate the distance from each point to a number of line features (rivers, roads etc.). Using the "Near" tool I can do this for one feature at a time (eg rivers) and I get 2 new fields in the attribute table indicating the id of the nearest feature and the distance. Fine. However, trying to find distance to the next features (eg roads) simply overwrites the previous ones (for rivers). The problem is I want to get distances to multiple features and get them all added separately to the attribute table. I tried to change the field name to see if this helped but I can only change the "Alias" so it makes no difference. I am working in Arcmap 10.0. Thanks.


texts: Creating a sections polygon layer from tic points

texts:
I am working in Arcmap v10.3.  I'm trying to create a "sections" (PLSS) layer from the section corners (tics).  I know I could convert the points to X,Y coordinates and use the coordinates to create lines.  But I can't figure out what order they will need to be arranged in the table to create squares from the lines.
Another route would be to take the state provided PLSS layer (which in inaccurate) and snap those polygon corners to my tics (the tics have been placed manually using imagery).


texts: How to use perpendicularoffset for PolygonSymbolizer in GeoServer?

texts:
I want to fill the polygon in GeoServer with 2 parallel lines pattern. Example: 1 continuous and 1 dash line. For example like this image:

I can fill the polygon in GeoServer with shape://slash. And combine 2 block , I can create a pattern with 2 lines. But I couldn't offset the second line to make them parallel (tried with perpendicularoffset but always got errors).
<PolygonSymbolizer>
  <Fill>
   <GraphicFill>
     <Graphic>
       <Mark>
         <WellKnownName>shape://slash</WellKnownName>
         <Stroke>
           <CssParameter name="stroke">#e31a1c</CssParameter>
           <CssParameter name="stroke-width">2</CssParameter>
         </Stroke>
       </Mark>
       <Size>20</Size>
     </Graphic>
   </GraphicFill>
 </Fill>
</PolygonSymbolizer>

<PolygonSymbolizer>
  <Fill>
   <GraphicFill>
     <Graphic>
       <Mark>
         <WellKnownName>shape://slash</WellKnownName>
         <Stroke>
           <CssParameter name="stroke">#001cac</CssParameter>
           <CssParameter name="stroke-width">2</CssParameter>
           <CssParameter name="stroke-dasharray">5 2</CssParameter>
         </Stroke>

        <!-- Not sure how to insert Offset correctly -->
         <PerpendicularOffset>5</PerpendicularOffset>

       </Mark>
       <Size>20</Size>
     </Graphic>
   </GraphicFill>
 </Fill>
</PolygonSymbolizer>

The sample use perpendicularoffset seem only work with stroke (boundary of polygons).
So is there anyway I can define a sld file for that purpose?
I can fill with .svg file to meet the requirement but the disadvantage is the difficult to change color of line.
In both ArcGIS & QGIS, I can easily define the pattern like that, so I don't understand why it seems not possible to do in GeoServer.


texts: Automating QGIS project composer printing?

texts:
I have a large project which consists of around 10x QGIS project files.  Each project contains a print composer with an atlas set up to generate up to 10 maps.

What is the easiest way to automate the opening of each project and creating a PDF using the composer atlas?

(Something like this would be ideal - not sure if it's possible)

C:\OSGeo4W\bin\qgis.bat --project MyProject1.qgs --code SaveAtlasAsPdf.py
C:\OSGeo4W\bin\qgis.bat --project MyProject2.qgs --code SaveAtlasAsPdf.py


texts:
I have a polyline shapefile with streams. I would like to turn this to a raster where there would be a cumulative length sum value in each pixel moving downstream. In other words, this would be an downstream increasing cumulative sum based on polyline length.

Unfortunately though, I'm unable to use DEM for calculating flow direction and do this through normal flow accumulation as the stream shapefile and DEM are slightly mismatching. Another problem is that the flow directions are set to both directions of the lines. Thus, I would like to first reset the flow direction of the polylines and then calculate the cumulative sum straight from the polyline file.

I would need 2 tools for this>

Tool that would redirect the line direction based on the DEM value of the two end points
Tool that would calculate cumulative length from the redirected polyline features.

Is anything like these existing?


texts: QGIS 3.4.2 SQL data not displaying on map

texts:
When I connect to SQL Server and tick 'Only look in the geometry_columns metadata table' it speeds up the scan but then when I add the feature class to the map it only appears in the layers list and not spatially on the map.

If I don't have this option ticked these features do display on the map, but the scan is much slower. It is not a problem with our cadastral data which we don't create, so I am wondering if the geometry isn't translating properly. We use Geomedia Professional as our GIS which handles geometries differently.

Josh


texts: OpenLayers seems slow when zooming / panning on vector (geometry line) from GeoJSON. Why?

texts:
I have a layer with 7 lines; when I zoom on some of the lines, the maps seems to be very slow.

I thought it might be because I bring a lot of points to the client (that represents the line).

aa
Line Id     Num of Points As Text
----------- ------------- ---------------------------
1           9971          MULTILINESTRING ((34.797771
2           13180         MULTILINESTRING ((34.808654
3           5465          LINESTRING (34.94746717 31.
4           5397          LINESTRING (34.94746717 31.
5           9550          MULTILINESTRING ((35.104847
6           3770          LINESTRING (34.94746717 31.
7           3770          LINESTRING (34.94746717 31.

Line 1 and 5 are very slow on the map... They are MultiLineString because I use MakeValid on the geometry from MSSQL.

Any Idea why this happens ?


texts: Buffer by field value in the processing modeler in QGIS

texts: I am using the processing modeler in QGIS 3.10. I want to buffer polygons by individual field attributes. This is not a problem when I am not using the processing modeller since one can simply choose the respective attribute field from the loaded layers. However, when I want to automize the workflow in the processing modller, I don't know how to tell QGIS from where it should take the attributes since there are nor preloaded layers in the processing modeller (see picture attached). I guess, I have to use the "Expression string builder", but I don't know which command I should use. Let's assume layer x contains the field y, which has the respective attributes values I want to use for buffering (which differ by polygon, so it's not a constant value). Any Ideas how I express that in the "Expression string builder", when using the processing modeler?


texts: Collector for ArcGIS PDOP and GPS collection configuration settings

texts:
My company just acquired a Geneq SXBlue II GPS receiver to pair with an iPad and the Collector for ArcGIS app.  The functionality is awesome and loved by everyone, but since we have  biologists on our team, they need to know the PDOP per point collected to ensure that it is below 6 to meet minimum standards.

I am basically wondering if/where there is the ability to control GPS collection configuration settings in Collector for ArcGIS, as described in the following support article about ArcGIS Mobile:

http://help.arcgis.com/en/arcgismobile/10.0/help/index.html#//007v0000000v000000.htm

I first started checking with the GPS receiver unit configuration settings to see if I could at least set a threshold to only collect points if the PDOP is under 6; however, I have been unsuccessful.  I next began to ponder that since I want these values per point collected, that maybe the setting was at the shapefile-level.


texts: I can't find where QGIS put my downloaded plugin - Crayfish

texts: I just downloaded the Crayfish plugin. It downloaded successfully but then when I search for it in QGIS, it is nowhere to be found. Can someone suggest where it might be located?


texts: Have a common filter via GET variable for multiple GeoJSON layers in QGIS

texts:
I have a Leaflet page, that fetches various GPS tracks and data from a website and filters the resulting GeoJSON files via HTTP GET variables.
If you want to get the data just for one day, I can use this URL:
https://my.service/geojson/tracker1?day=2020-12-18
I need to view some data now as well in QGIS. I can add the layers fine, but when I want to change the date, I have to manually edit the date in every layer.
Is there any possibility to have a textbox, where I can add the date and it will be applied to all GeoJSON layers.
I have looked at the filtering, but if I export all tracker points, it takes ages to load the data.


texts: Changing table name of command using raster2pgsql?

texts:
I am completely newbie with GIS and PostGIS as well

I'm trying to load raster into postgreSQL with raster2pgsql. It runs well.

However, the table name of the raster is generated automatically. For example, the raster named: landsat.tiff, then the table name is landsat.

Can we change that name of table? and how we can do it?

Here's what I typed in Windows Command line:

\raster2pgsql -s 4326 *.TIFF -t 100x100 > demo.sql

and insert into database

\psql -U postgres -d demo -f demo.sql


texts: Making Default attributes warning popup go away?

texts:
I've got multiple county feature classes within a feature dataset. Within one of the county features I created a new field with coded subtypes. Whenever I'm editing this county in Arc and change this new fields attributes to one of these subtypes I get this pop-up:

  You are changing the subtype for this feature. Do you want to change
  the default attributes for the feature to those defined by the new
  subtype?

Any ideas why this happens and how to make it stop?


texts: Building larger polygons from smaller ones *without* common ID in PostGIS?

texts:
I am looking for a method, using PostGIS, of reconstituting a set of polygons from tens of thousands of smaller polygon 'pieces'.

This is a simplified example, with the pieces shaded:

Unfortunately, the pieces do not have a common ID or other attribute of the larger polygons that would allow me to group them. In some cases, I have the grid that was used to cut the pieces, but in other cases I do not.

The only approach I've been able to devise -- and I haven't been able to make it work yet -- is to use ST_BuildArea on pairs of intersecting pieces in an iterative process as follows (renaming 'composite' to 'pieces' between iterations):

SELECT p1.id,ST_BuildArea(ST_Collect(p1.the_geom,p2.the_geom)) AS the_geom
INTO composite
FROM pieces p1 LEFT OUTER JOIN pieces p2
    ON (ST_DWithin(p1.the_geom,p2.the_geom,1))
WHERE ST_Intersects(p1.the_geom,p2.the_geom) AND
p1.id < p2.id;

This generates larger and larger composite pieces, but many of them (nearly 1,300 after three iterations on the example), and they come increasingly to overlap. It seems to me that to go much further I'd need a means of paring down the results each iteration. (I'm also getting side location conflicts on the fourth iteration, but that may be a separate problem.)


texts: Changing layer's ArcSDE data source user name?

texts:
I am a new hire and I'm in the process up setting up my database connections and default instance. I set up my personal database connection to the server under my name.

I received a MXD from a coworker and saved it on my computer. In the layer properties source tab, the data source shows that the user name is still under my coworkers name.

I would like to change the user name to my name in the source data.

Is there anyway to change the user name for this, and perhaps all of the layers at once?


texts:
Here is my code
How do I add the soil moisture data I am trying to get onto my Google drive?

var dataset = ee.ImageCollection('NASA_USDA/HSL/SMAP_soil_moisture')
                  .filter(ee.Filter.date('2017-04-01', '2017-04-30'));
var soilMoisture = dataset.select('ssm');
var soilMoistureVis = {
  min: 0.0,
  max: 28.0,
  palette: ['0300ff', '418504', 'efff07', 'efff07', 'ff0303'],
};
Map.setCenter(-6.746, 46.529, 2);
Map.addLayer(soilMoisture, soilMoistureVis, 'Soil Moisture');


texts: Rotating vector grid in QGIS 3.10 vs 3.8

texts:
Just wondering if anyone else has had an issue with vector rotation via the advanced editing tools in 3.10?

I have a vertical grid overlaying Vancouver Island that I had tried to rotate using the advanced edit tools in QGIS 3.10. Rotating the grid sends QGIS into spinning wheel of infinite loops.

Doing the same task in 3.8 is not a problem.

Just thought I would check here before submitting a bug.

Here are the GeoPackage files in case anyone wants to test.

https://drive.google.com/file/d/1H8mld8fLWkqnqUqpVa1dlkhnuw1H1kAe/view?usp=sharing


texts: r.viewshed - no output file saved or loaded

texts: I can successfully run r.viewshed (accordingly to the log) but there is an error when loading the output file. Running without the autoload of the output removes the error but doesn't appear to save an output file anywhere.


texts: ArcGIS 10 Locators: TA_Address_NA_10 vs TA_Streets_US_10 -- which to use?

texts:
I see there are two locators on tasks.arcgisonline.com, TA_Address_NA_10 vs TA_Streets_US_10. I'm having difficulty understanding the difference between the two.

From browsing the properties it seems that the "Address" locator is a composite which contains a whole bunch of locators, using Rooftop as first and "Streets" second. I'm guessing that Rooftop usually provides the most accurate location for an address. Is that true?

Does it follow that you should always use "Address"? Why would you ever use the "Streets" locator?

thanks


texts: what kind of accuracy (rmse) can be expected by georeferencing a scanned image?

texts:
I am trying to understand what is the effective accuracy of a georeferenced image and how to attach a standard to the georeferencing process which can be used for future reference. To explain further, I am trying to assess what kind of RMSE would be acceptable when for example georeferencing an accurate architects drawing with an aerial image, or with a small scale landuse plan drawn with thick outlines using parcel data with an accuracy of plus/minus 1cm.

If I am using an inaccurate base map in which to base my passpoints, does the rms error have any meaning whatsoever? Logically, is it not so, that one can only use the rmse as a standard when the accuracy of the basemap can be certified? If the base map is not accurate, any passpoints referenced to it would be even less accurate and an RMS error meaningless.

Am I on the right track?


texts: Zonal Statistics sum of everything touched by polygon

texts:
I'm trying to find the sum of all raster values unweighted, that touch my polygon layer.
In the picture I've below that would give me a sum of 7.4, instead it calculates it as 6.55 because the top left raser value is not added to the selected polygons sum. I've read somewhere that Zonal Statistics uses the center of either raster or vector layer to determine if it is added. I was wondering if there is a way around this.
This is a very large dataset 1.8 gigabyte in a 1.6 x 1.6 m grid covering a whole municipality. Converting everything to vectors would i believe is nearly impossible and using a buffer would probably get false results. The polygon layer represents bluespots with water, and the purpose is to get the sum of the rasters onto the polygon, without weighing the sum. The polygons differ a lot in size, from few cells to hundreds of cells.

I have access to Spatial Analyst, 3D Analyst, and ArcGIS 10.7 Advanced Desktop license.
Any suggestions for an alternativ approach, or way around how Zonal Statistics work?


texts: Getting multipolygon vertexes to points in PostGIS

texts:
I want to extract the coordinates of the vertices of a polygon, and then create a point on each of the vertices.
Sample of Polygon geometry table (multipolygon):

Wished output:


texts: Last moment a layer in Geoserver was viewed

texts:
My company publishes >1,000 layers using several geoservers. We are publishing them on several different websites.
Because there are so many layers, it is impossible for us to know which layer is being used and which one is not.
Therefore, I was wondering whether it is possible to view the last time a layer in Geoserver was requested. Is there any option to know this?


texts: Snap lines to polygon verticies at interscetion

texts:
I'm trying to figure out a way to snap lines to a polygon layer, but only at the intersecting point.  Here is a sample picture of what I'm working with.  I have a layer of polygons and a line string layer.

My goal is to create a topology from the polygons, thin the verticies, and then snap the line layer to polygon verticies within some threshold (i.e. 20 meters).  However, I haven't figured out a way to do this using any of my available software, which is QGIS, Grass 7, and PostGIS.  Ultimately I would like to be able to make changes to say the polygon layer, and have those changes affect the line string layer, but still be able to differentiate the layers and all their attributes.  I have tried various tools in GRASS to no avail, and have only just begun to learn and investigate PostGIS.

Does anyone have suggestions on ways to accomplish this task using these open source tools?  Even multi-step, more complex solutions would be welcome, as I hope to eventually script this work flow as much as possible.


texts: QField: Points aren't straight, they are dirty and blurred

texts:
I've made a project in QGIS, imported a new layer with coordinates (points) from a csv-file. Finally I want to find this marks in the area within 10 cm.
In QGIS everything worked fine.
When I made the QField-Package, it works, I can open it in QField.
But the points are displayed not clear. It looks like, they've got a new dirty, blurred background. I already tried everything, to lose it. But I don't find a possibility.

Does anybody have got an idea?


texts: How should I normalize my data - based on population of an area /neighbourhood or population of the entire city?

texts:
I am studying the driving's forces of crime in a particular city.

I have some independent values, e.g. percentage of the head of the household with low education per zone, the percentage of unemployed people per zone [...]

The percentage is related to the total population in each zone, and not from the total population of the whole city. I think if I use the absolute values it will give me some correlation (multicollinearity) among independent variable since that zone with more population tends to have also a higher number of people with low education, unemployed people and so on.

For example, instead of using the absolute value of the number of people with low education in a particular zone, maybe I should divide that value to the total population of that specific zone. This might help to avoid that such zones with higher population will not necessarily have a higher number of people with low level of education. Because a have difference in the size of zones and consequently the difference in their population.


texts:
Newbie in ArcObjects.
Would like to expedite the learning process by looking at existing code samples.

Please paste links which would provide code samples for common tasks.


texts: Javascript Convert WMS url to coordinate?

texts:
I'm looking for some hints (Javascript custom function or openlayers 4 function) ...) that can help me to convert a wms response URL parameters to get my click coordinates.

Here an example of url to convert to a click coordinate.

https://ws.mapserver.transports.gouv.qc.ca/swtq?SERVICE=WMS
&VERSION=1.3.0
&REQUEST=GetFeatureInfo
&BBOX=-7859605.72457123640924692,5954835.61684618983417749,-7802662.21180228982120752,5988064.99139707256108522
&CRS=EPSG:3857
&WIDTH=1155
&HEIGHT=674
&LAYERS=histo_stations_meteoroutieres
&STYLES=
&FORMAT=image/png
&QUERY_LAYERS=histo_stations_meteoroutieres
&INFO_FORMAT=text/html
&I=722
&J=403
&FEATURE_COUNT=10

My goal is to recreate a gml from this query based on derived wkt Point.


texts: Do you know that is possible to buy commercial license on Geoserver? I know that is under GPL 2.0 license and any changes have to be distributed also under this license. I would like to use Geoserver for commercial project.


texts: Understanding Cost Surface algorithm in ArcGIS Desktop?

texts:
Based on this article "PDF: Extensions to least-cost path algorithms for roadway planning (Chaoqing, Jay, dan Mandy 2002)" (Click on link for PDF).

Whether least cost path tool used in ArcGIS software is already using anisotropic accumulated surface?

If so, why in determining the cost surface, the slope of the data used is the result of the  3D slope analysis?

Whereas based "ArcGIS Help sources" equation used is as shown below:

Whereas the anisotropic equation-accumulated surface, the formula used is the Euclidean Distance formula.

If anyone can explain what alogarithm used for the analysis of least cost path in ArcGIS softwere?

Extension:

In this case, I want to find the shortest path between two points through the terrain surface. Technical criteria that must be fulfilled is, "the slope should be less than 5%" when I compare with other software by hand through contour lines such as AutoCAD, Surpac. the results are very different. see the picture below


texts: Joining or relating multiple data in one single data using QGIS

texts:
I do not know how to solve this problem using QGIS. I have tried several ways but it did not work out.
I have a lot of data, I have one coordinate system for each location, for example

Alor setar with coordinate (6.200,100.400),
batu embun (3.9666,102.3499).

Then for Alor Star there are consist a lot of data based on different year which is for 2012 2011 2010 and other. Then in 2012 2011 and 2010 there also separate data based on month which January February, March, June and other.
How can I solve this data and combine it?
Either use one to many or many to one and how I can relate all the data.


texts: How can I query Postgres tables directly from OpenLayers?

texts:
I have a map that is being hosted with the components of OpenGeo. The map contains WMS layers of millions of parcels and some rail stops. All of our spatial data is on PostgresSQL.

Right now the only interaction I have with Postgres is via OpenLayers/OpenGeo. I'd like to be able to send queries to the database via code (I guess JavaScript or PHP) so I can better interact with the data. For example, I'd like to be able to send a user entered parcel ID number to the database so I can get the latitude and longitude of the parcel, and then I then use OpenLayers to zoom to it. I'd also like to be able to get a list of the rail stop names and populate a drop down menu with them, then zoom to the rail stop when the user selects it in the drop down. I can write the queries I'd need in SQL and execute it in a console, but I am lost as to how to use PHP or JavaScript, etc to actually execute the query and return results. Does anyone have any suggestions or pointers as to where I can start to learn what I need to do.

Thank you very much for your help.


texts: Is there Data Loader for QGIS or Spatialite?

texts:
Is there a tool that works like the data loader of ArcGIS for QGIS (or via Spatialite )?
To integrate data in a database or a shapefile without having the same structure (and the tool allows to propose which column goes with every other column)

For example, we have a layer of trees; two people go on the field and inform each new tree layer (without necessarily having the same structure / field names).
They come back ; I want to integrate their new trees quickly and easily without having to explicitly rename their fields: ArcGIS allows in a geodatabase with the tool "load data" that makes a direct link between the fields of the two layers (linking with a gui).


texts: Plotting lat/long with Tableau averages it out

texts:
I'm trying to plot crimes by lattitude and longitude on Tableau Public. Yet the software averages it out the values:

I've watch this tutorial but, as you can see on the screenshot, I don't have the "Dimensions" feature.


texts: AttributeTable sending an alert that data was not found when clicking on table row

texts:
I have an AttributeTable that I have opening when I start a widget.  When I click on a table row in this instance, I am able to get the information back correctly and not see this message.  I did this by changing the config.json file to have the primary field in the database to not show.

"name": "Current Parties",
"id": "Layer Id",
"layer":{
    "url": "url to Feature Server",
    "fields":[
        {
            "show": false,
            "name": "OBJECTID",
            "alias": "OBJECTID"
        }
    ]
}

With the rest of the fields having the property set to true.  However, when I change the table and have it pulling back from a different database, I see the correct information, but when I click on a record in the table, I get this alert

Which should not be appearing since the new information is being displayed in exactly the same way.  Can anyone explain why this alert is still showing up?  If anything else needs to be added to better help diagnose the problem, let me know and I'll add it.


texts: Calculate neighborhood statistics for polygons?

texts:
I've been searching around for the answer to this (including here) and have yet to find it.  I'm interested in finding out how to calculate a 'neighborhood sum' for vector, polygon data.  I am aware of how to do this with rasters- the equivalent, then, would be to compute a focal sum for each polygon.

I use ArcView 10.

Thanks for your assistance!


texts: Exporting KML/ SHP/ EXCEL Files from Google Maps [Android version]

texts: I want to export my pin locations using Google Maps [ Android version] and save the layer as KML /SHP/EXCEL in my phone. Is this feasible using Google Maps or are there any other applications which can support this functionality?


texts: I am using i.atcorr for atmospheric correction in Grass6.4.4. I would like to add Landsat8 sensor to my project(build). If I copy the iwave.cpp file from grass 7, I get error messages. Do I need to compile the entire Grass7 i.atcorr module again for this or is there any alternate?


texts: Multiple logical filters in a WFS request with OpenLayers 2.13.1

texts:
Is it possible, and if so how, to use multiple logical filters in a wfs request?

For example I have a attribute, type_id that should be 1, 2 or 3.

var wfsTypeFilters = [
    new OpenLayers.Filter.Comparison({
        type: OpenLayers.Filter.Comparison.EQUAL_TO,
        property: "type_id",
        value: 1
    }),
    new OpenLayers.Filter.Comparison({
        type: OpenLayers.Filter.Comparison.EQUAL_TO,
        property: "type_id",
        value: 2
    }),
    new OpenLayers.Filter.Comparison({
        type: OpenLayers.Filter.Comparison.EQUAL_TO,
        property: "type_id",
        value: 3
    })
];

And then I have an attribute called name that should not be empty

var wfsNameIsNotNullFilter = [
    new OpenLayers.Filter.Comparison({
        type: OpenLayers.Filter.Comparison.NOT_EQUAL_TO,
        property: "type_id",
        value :""
    }
)];

After thtis I tried to create two different logical filters. And when combining them I get error:
Uncaught Filter writing not supported for rule type: undefined

var parent_type_filter = new OpenLayers.Filter.Logical({
    type: OpenLayers.Filter.Logical.OR,
    filters: wfsTypeFilters
});

var parent_filter = new OpenLayers.Filter.Logical({
    type: OpenLayers.Filter.Logical.AND,
    filters: [parent_type_filter, wfsNameIsNotNullFilter]
});

If I would have written a sql query this would have been easy by using parentheses

WHERE (type_id = 1 OR type_id = 2 OR type_id = 3 ) AND (name != "")

Br
Mike


texts: Manipulating Transformation parameters with PyQGIS

texts:
Is it possible to add a layer with a local coordinate system, transformed on the fly by a custom transformation entry from this local CRS -towgs84 and furthermore to manipulate the -towgs84 Parameters?

I would like to Helmert-Transform a Polyline from a geodetic survey by parameters I calculate with python based on ground control points.

I am aware of GDALwarp (only works with rasters) and the Georeferencer Plugin, which unfortunately does not support 7-Parameter Helmert Transformation.


texts: How can I define a threshold or outlier for the elevation difference between two different time period taken in rock, that should be zero?

texts:
I am doing the geodetic mass balance of a glacier and I have the elevation difference that was produced in the glacier, but now I have to correct this data but making the mean of the elevation difference in the rock outside the glacier, zero. The difficulty is to know what is rock, because it can get confused with a rocky glacier, glacier, and where there was, for example, a collapse that greatly affects the data. I need only the rock, that is my valid cell to take the mean and the substract this mean to all the dem of glacier difference.
So I extracted everything that was outside the glacier in the DEM that is the elevation difference of
the intersection of two time periods, and I supose everything here is rock, so in this case for example, rock measured in the year 2009 and the rock measured in year 2015 should be the same and the difference should be zero. But in the data there is error between the two measurements or the area I am evaluating, and the values range between -28 m and 25 m. So I need to delete this points that are not valid cells to get the mean that represent the error in valid cells of rock, and that has a low standard deviation, so I need to define a threshold to do so.
I want to know a statistical method to define this threshold or outliers that give me certainty that what I am deleting is not statistically significant. I remember that I did this in university but for the waves in the ocean and we did it with tests: gumbel, rayleigh and frechet. But I wish to do it on Qgis before looking a way that make me export more than 4 million points to R.
I leave an image of the raster that covers the rock area of the elevation difference (black and white), and also the statistic of this raster.


texts: Creating new district layer from parcels layer

texts:
The idea of my script is to create a new layer which will be a layer of districts in the city based on a layer of parcels. The attribute table of parcels has one field that should be used to determine which parcels are from which district.

I managed to do it for one district, but I have no idea how to make it work for all of them, no matter how many parcels or districts are in layer.

import arcpy

arcpy.env.workspace = r"D:\Studia\python'\roboczapy.gdb"

##City shapefile
obreb = r"D:\Studia\python'\P3_Dane\Gdansk_22.shp"

cursor = arcpy.SearchCursor (obreb, fields="IDENTYFIKA")

##Creating a new layer
nowawarstwa = arcpy.MakeFeatureLayer_management(obreb, out_layer="Gdansk_22_Lyr", field_info="FID FID VISIBLE NONE;Shape Shape VISIBLE NONE;OBJECTID OBJECTID VISIBLE NONE;IDENTYFIKA IDENTYFIKA VISIBLE NONE;POWIERZCHN POWIERZCHN VISIBLE NONE;TERYT TERYT VISIBLE NONE;NUMER NUMER VISIBLE NONE;WOJEWODZTW WOJEWODZTW VISIBLE NONE;POWIAT POWIAT VISIBLE NONE;GMINA GMINA VISIBLE NONE;DATA_OD DATA_OD VISIBLE NONE;Shape_Leng Shape_Leng VISIBLE NONE;Shape_Area Shape_Area VISIBLE NONE;obreb obreb VISIBLE NONE;obreby obreby VISIBLE NONE")
arcpy.SaveToLayerFile_management(nowawarstwa, "Gdansk22Lyr", "ABSOLUTE")

## Selecting and exporting to another layer only those parcels from first district using clause [where_clause=""""IDENTYFIKA" LIKE '226101_1.0001.%'"""]
obreb1 = arcpy.FeatureClassToFeatureClass_conversion(obreb, out_path="D:/Studia/python'/roboczapy.gdb", out_name="obreb1", where_clause=""""IDENTYFIKA" LIKE '226101_1.0001.%'""", field_mapping="""OBJECTID "OBJECTID" true true false 11 Double 0 11 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,OBJECTID,-1,-1;IDENTYFIKA "IDENTYFIKA" true true false 100 Text 0 0 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,IDENTYFIKA,-1,-1;POWIERZCHN "POWIERZCHN" true true false 19 Double 15 18 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,POWIERZCHN,-1,-1;TERYT "TERYT" true true false 7 Text 0 0 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,TERYT,-1,-1;NUMER "NUMER" true true false 50 Text 0 0 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,NUMER,-1,-1;WOJEWODZTW "WOJEWODZTW" true true false 254 Text 0 0 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,WOJEWODZTW,-1,-1;POWIAT "POWIAT" true true false 254 Text 0 0 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,POWIAT,-1,-1;GMINA "GMINA" true true false 254 Text 0 0 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,GMINA,-1,-1;DATA_OD "DATA_OD" true true false 8 Date 0 0 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,DATA_OD,-1,-1;Shape_Leng "Shape_Leng" true true false 19 Double 15 18 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,Shape_Leng,-1,-1;Shape_Area "Shape_Area" true true false 19 Double 15 18 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,Shape_Area,-1,-1;obreb "obreb" true true false 5 Short 0 5 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,obreb,-1,-1;obreby "obreby" true true false 5 Short 0 5 ,First,#,D:\Studia\python'\P3_Dane\Gdansk_22.shp,obreby,-1,-1""", config_keyword="")

## Dissolving this one selected district with additional information about average area of parcells, STD of parcels are - but there should be median and highest parcel number in district
arcpy.Dissolve_management(in_features="obreb1", out_feature_class="D:/Studia/python'/roboczapy.gdb/dissolveobreb1", dissolve_field="", statistics_fields="POWIERZCHN MEAN;POWIERZCHN STD;IDENTYFIKA MAX", multi_part="MULTI_PART", unsplit_lines="DISSOLVE_LINES")

print "end"

Some things here look terrible because Copyied Snippets from ArcMap to Python.

Despite dissolving, the attribute table in new layer should contain fields like district number, average parcel area, median of parcel area and highest parcel number, some of these I managed to do in dissolve tool, but that's not enough.

I know I can do it in ModelBuilder or step by step, and it may be easier for me. But unfortunately I HAVE TO do it as a Arcpy/Python code.